The class of generating functions for completely monotone sequences (moments of finite positive measures on $[0,1]$) has an elegant characterization as the class of Pick functions analytic and positive on $(-\infty,1)$. We establish this and another such characterization and develop a variety of consequences. In particular, we characterize generating functions for moments of convex and concave probability distribution functions on $[0,1]$, and we provide a simple analytic proof that for any real $p,r>0$, the Fuss-Catalan or Raney numbers $\frac{r}{pn+r}\binom{pn+r}{n}$, $n=0,1,\ldots$ are the moments of a probability distribution on some interval $[0,\tau]$ if and only if $p\ge1$ and $p\ge r$.

From a variational perspective, we derive a series of magnetization and quantum spin current systems coupled via an “s-d” potential term, including the Schrödinger--Landau--Lifshitz--Maxwell system, the Pauli--Landau--Lifshitz system, and the Schrödinger--Landau--Lifshitz system with successive simplifications. For the latter two systems, we propose using the time splitting spectral method for the quantum spin current and the Gauss--Seidel projection method for the magnetization. Accuracy of the time splitting spectral method applied to the Pauli equation is analyzed and verified by numerous examples. Moreover, behaviors of the Schrödinger--Landau--Lifshitz system in different “s-d” coupling regimes are explored numerically.

Abstract We investigate a system of N Brownian particles with the Coulomb interaction in any dimension \(d\ge 2\), and we assume that the initial data are independent and identically distributed with a common density \(\rho _0\) satisfying \(\int _{\mathbb {R}^{d}}\rho _0\ln \rho _0\,\hbox {d}x<\infty \) and \(\rho _0\in L^{\frac{2d}{d+2}} (\mathbb {R}^{d}) \cap L^1(\mathbb {R}^{d}, (1+|x|^2)\,\hbox {d}x)\). We prove that there exists a unique global strong solution for this interacting partsicle system and there is no collision among particles almost surely. For \(d=2\), we rigorously prove the propagation of chaos for this particle system globally in time without any cutoff in the following sense. When \(N\rightarrow \infty \), the empirical measure of the particle system converges in law to a probability measure and this measure possesses a density which is the unique weak solution to the mean-field Poisson–Nernst–Planck equation of single component.

In this paper, we discuss error estimates associated with three different aggregation-diffusion splitting schemes for the Keller-Segel equations. We start with one algorithm based on the Trotter product formula, and we show that the convergence rate is CΔt, where Δt is the time-step size. Secondly, we prove the convergence rate CΔt2 for the Strang's splitting. Lastly, we study a splitting scheme with the linear transport approximation, and prove the convergence rate CΔt.

In this paper we study the existence assertion of the initial boundary value problem for the equation $\frac{\partial u}{\partial t} = \Delta e^{-\Delta u}$. This problem arises in the mathematical description of the evolution of crystal surfaces. Our investigations reveal that the exponent in the equation can have a singular part in the sense of the Lebesgue decomposition theorem, and the exponential nonlinearity somehow “cancels” it out. The net result is that we obtain a solution $u$ that satisfies the equation and the initial boundary conditions in the almost everywhere (a.e.) sense.