Markov processes conditioned on their location at large exponential times

Suppose that ${\left(X_t\right)}_{t\ge 0}$ is a one-dimensional Brownian motion with negative drift $?\mu$. It is possible to make sense of conditioning this process to be in the state 0 at an independent exponential random time and if we kill the conditioned process at the exponential time the resulting process is Markov. If we let the rate parameter of the random time go to 0, then the limit of the killed Markov process evolves like $X$ conditioned to hit 0, after which time it behaves as $X$ killed at the last time $X$ visits 0. Equivalently, the limit process has the dynamics of the killed “bang–bang” Brownian motion that evolves like Brownian motion with positive drift $+\mu$ when it is negative, like Brownian motion with negative drift $?\mu$ when it is positive, and is killed according to the local time spent at 0.An extension of this result holds in great generality for a Borel right process conditioned to be in some state $a$ at an exponential random time, at which time it is killed. Our proofs involve understanding the Campbell measures associated with local times, the use of excursion theory, and the development of a suitable analogue of the “bang–bang” construction for a general Markov process.As examples, we consider the special case when the transient Borel right process is a one-dimensional diffusion. Characterizing the limiting conditioned and killed process via its infinitesimal generator leads to an investigation of the $h$-transforms of transient one-dimensional diffusion processes that goes beyond what is known and is of independent interest.     

Existence of strong minimizers for the Griffith static fracture model in dimension two

We consider the Griffith fracture model in two spatial dimensions, and prove existence of strong minimizers, with closed jump set and continuously differentiable deformation fields. One key ingredient, which is the object of the present paper, is a generalization to the vectorial situation of the decay estimate by De Giorgi, Carriero, and Leaci. This is based on replacing the coarea formula by a method to approximate $SB{D}^p$ functions with small jump set by Sobolev functions, and is restricted to two dimensions. The other two ingredients will appear in companion papers and consist respectively in regularity results for vectorial elliptic problems of the elasticity type and in a method to approximate in energy $GSB{D}^p$ functions by $SB{V}^p$ ones.     

Optimizing and improving of the C-to-R method for solving complex symmetric linear systems

For complex symmetric linear systems, Axelsson et al. (2014) proposed the C-to-R method. In this paper, by further studying the C-to-R method with $W$ and $T$ being symmetric positive semidefinite, the optimal iteration parameter for the C-to-R method ${\alpha }_{opt}=\frac{\sqrt{2\sqrt{2}}}{2}$ is obtained and the C-to-R method is optimized. Furthermore, based on the optimized C-to-R method, we further propose an optimized preconditioner. Eigenvalue properties of the optimized preconditioned matrix are analyzed, which show that all the eigenvalues of the preconditioned matrix are located in tighter interval. Numerical results are presented, not only confirm the validity of the theoretical analysis, but also demonstrate the feasibility and effectiveness of the proposed optimized C-to-R method.     

Collision and escape orbits in a generalized Hénon–Heiles model

The motion of a material point of unit mass in a field determined by a generalized Hénon–Heiles potential $U=A{q}_1^2+B{q}_2^2+C{q}_1^2{q}_2+D{q}_2^3$, with $(q_1,q_2)=$ standard Cartesian coordinates and $(A,B,C,D)\in {(0,\infty )}^2\times R^2$, is addressed for two limit situations: collision and escape. Using McGehee-type transformations, the corresponding collision and infinity boundary manifolds pasted on the phase space are determined. These are fictitious manifolds, but, due to the continuity with respect to initial data, their flow determines the near by orbit behaviour.The dynamics on the collision and infinity manifolds is fully described. The topology of the flow on the collision manifold is independent of the parameters. In the full phase space, while spiraling collision orbits are present, most of the orbits avoid collision. The topology of the flow on the infinity manifold changes as the ratio between $C$ and $D$ varies. More precisely, there are two symmetric pitchfork bifurcations along the line $2C?3D=0$, due to the reshaping of the potential along the bifurcation line. Besides rectilinear and spiraling orbits, the near-escape dynamics includes oscillatory orbits, for which angular momentum alternates sign.     

Global existence and time decay of the non-cutoff Boltzmann equation with hard potential

This paper is concerned with the Cauchy problem on the Boltzmann equation without angular cutoff assumption for hard potential in the whole space. When the initial data is a small perturbation of a global Maxwellian, the global existence of solution to this problem is proved in unweighted Sobolev spaces $H^N\left({\mathbb{R}}_{x,v}^6\right)$ with $N\ge 2$. But if we want to obtain the optimal temporal decay estimates, we need to add the velocity weight function, in this case the global existence and the optimal temporal decay estimate of the Boltzmann equation are all established. Meanwhile, we further gain a more accurate energy estimate, which can guarantee the validity of the assumption in Chen et al. (0000).     

Lattice BGK simulations of double diffusive natural convection in a rectangular enclosure in the presences of magnetic field and heat source

In this paper, we develop a temperature-concentration lattice Bhatnagar-Gross-Krook (TCLBGK) model, with a robust boundary scheme for simulating the two-dimensional, hydromagnetic, double-diffusive convective flow of a binary gas mixture in a rectangular enclosure, in which the upper and lower walls are insulated, while the left and right walls are at a constant temperature and concentration and a uniform magnetic field is applied in the $x$-direction. In the model, the velocity, temperature and concentration fields are solved by three independent LBGK equations which are combined into a coupled equation for the whole system. In our simulations, we take the Prandtl number $Pr=1$, the Lewis number $Le=2$, the thermal Rayleigh number $R{a}_T={10}^5,{10}^6$, the Hartmann number $Ha=0,10,25,50$, the dimensionless heat generation or absorption $?=0.0,?1.0$, the buoyancy ratio $N=0.8,1.3$, and the aspect ratio $A=2$ for the enclosure. The numerical results are found to be in good agreement with those of previous studies [A.J. Chamkha, H. Al-Naser, Hydromagnetic double-diffusive convection in a rectangular enclosure with opposing temperature and concentration gradients, Int. J. Heat Mass Transfer 45 (2002) 2465C2483].     

Solving bivariate systems using Rational Univariate Representations

Given two coprime polynomials $P$ and $Q$ in $\mathbb{Z}[x,y]$ of degree bounded by $d$ and bitsize bounded by $\tau$, we address the problem of solving the system $\{P,Q\}$. We are interested in certified numerical approximations or, more precisely, isolating boxes of the solutions. We are also interested in computing, as intermediate symbolic objects, rational parameterizations of the solutions, and in particular Rational Univariate Representations (RURs), which can easily turn many queries on the system into queries on univariate polynomials. Such representations require the computation of a separating form for the system, that is a linear combination of the variables that takes different values when evaluated at the distinct solutions of the system.We present new algorithms for computing linear separating forms, RUR decompositions and isolating boxes of the solutions. We show that these three algorithms have worst-case bit complexity ${\tilde{O}}_B(d^6+d^5\tau )$, where $\tilde{O}$ refers to the complexity where polylogarithmic factors are omitted and $O_B$ refers to the bit complexity. We also present probabilistic Las Vegas variants of our two first algorithms, which have expected bit complexity ${\tilde{O}}_B(d^5+d^4\tau )$. A key ingredient of our proofs of complexity is an amortized analysis of the triangular decomposition algorithm via subresultants, which is of independent interest.     

Optimal lower eigenvalue estimates for Hodge-Laplacian and applications

We consider the eigenvalue problem for Hodge-Laplacian on a Riemannian manifold M isometrically immersed into another Riemannian manifold $\overline{M}$. We first assume the pull back Weitzenböck operator of $\overline{M}$ bounded from below, and obtain an extrinsic lower bound for the first eigenvalue of Hodge-Laplacian. As applications, we obtain some rigidity results. Second, when the pull back Weitzenböck operator of $\overline{M}$ bounded from both sides, we give a lower bound of the first eigenvalue by the Ricci curvature of M and some extrinsic geometry. As a consequence, we prove a weak Ejiri type theorem, that is, if the Ricci curvature bounded from below pointwisely by a function of the norm square of the mean curvature vector, then M is a homology sphere. In the end, we give an example to show that all the eigenvalue estimates are optimal when $\overline{M}$ is the space form.     

Degree-equipartite graphs

A graph $G$ of order $2n$ is called degree-equipartite if for every $n$-element set $A?V\left(G\right)$, the degree sequences of the induced subgraphs $G\left[A\right]$ and $G[V\left(G\right)?A]$ are the same. In this paper, we characterize all degree-equipartite graphs. This answers Problem 1 in the paper by Grünbaum et al. [B. Grünbaum, T. Kaiser, D. Král, and M. Rosenfeld, Equipartite graphs, Israel J. Math. 168 (2008) 431–444].     

The Cauchy problem for a generalized Camassa–Holm equation

In this paper, we are concerned with the Cauchy problem of the generalized Camassa–Holm equation. Using a Galerkin-type approximation scheme, it is shown that this equation is well-posed in Sobolev spaces $H^s$, $s>3/2$ for both the periodic and the nonperiodic case in the sense of Hadamard. That is, the data-to-solution map is continuous. Furthermore, it is proved that this dependence is sharp by showing that the solution map is not uniformly continuous. The nonuniform dependence is proved using the method of approximate solutions and well-posedness estimates. Moreover, it is shown that the solution map for the generalized Camassa–Holm equation is Hölder continuous in $H^r$-topology. Finally, with analytic initial data, we show that its solutions are analytic in both variables, globally in space and locally in time.     

An Edgeworth expansion for functionals of Gaussian fields and its applications

This paper is concerned with the rate of convergence in the normal approximation of the sequence $\left\{F_n\right\}$, where each $F_n$ is a functional of an infinite-dimensional Gaussian field. We develop new and powerful techniques for computing the exact rate of convergence in distribution with respect to the Kolmogorov distance. As a tool for our works, the Edgeworth expansion of general orders, with an explicitly expressed remainder, will be obtained, and this remainder term will be controlled to find upper and lower bounds of the Kolmogorov distance in the case of an arbitrary sequence $\left\{F_n\right\}$. As applications, we provide the optimal fourth moment theorem of the sequence $\left\{F_n\right\}$ in the case when $\left\{F_n\right\}$ is a sequence of random variables living in a fixed Wiener chaos or a finite sum of Wiener chaoses. In the former case, our results show that the conditions given in this paper seem more natural and minimal than ones appeared in the previous works.