On the microcanonical solution of a system of fully coupled particles

We study the Hamiltonian mean field (HMF) model, a system of N fully coupled particles, in the microcanonical ensemble. We use the previously obtained free energy in the canonical ensemble to derive entropy as a function of energy, using Legendre transform techniques. The temperature–energy relation is found to coincide with the one obtained in the canonical ensemble and includes a metastable branch which represents spatially homogeneous states below the critical energy. “Water bag” states, with removed tails momentum distribution, lying on this branch, are shown to relax to equilibrium on a time which diverges linearly with N in an energy region just below the phase transition.     

Minimum Energy Problem on the Hypersphere

We consider the minimum energy problem on the unit sphere \(\mathbb {S}^{d-1}\) in the Euclidean space \(\mathbb {R}^{d}\), d = 3, in the presence of an external field Q, where the charges are assumed to interact according to Newtonian potential 1/r ^d-2, with r denoting the Euclidean distance. We solve the problem by finding the support of the extremal measure, and obtaining an explicit expression for the density of the extremal measure. We then apply our results to an external field generated by a point charge of positive magnitude, placed at the North Pole of the sphere, and to a quadratic external field.     

Localization of the essential spectrum of the energy operators ofn-particle quantum systems in a magnetic field

Theorems on the localization of the essential spectrum of givenSO(2) symmetry forn-particle Hamiltonians in an external magnetic field are proved. The energy operators are studied for systems of arbitrary particles in the absence of an external potential field and for systems of identical particles in an external potential field. Some of the results were announced in earlier papers [1,2].Radiophysics Research Institute. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 97, No.1, pp. 94–112, October, 1993.     

Equivalence of the Duffin-Kemmer-Petiau and Klein-Gordon-Fock equations

A strict proof of the equivalence of the Duffin-Kemmer-Petiau and Klein-Gordon-Fock theories is presented for physical S-matrix elements in the case of charged scalar particles minimally interacting with an external or quantized electromagnetic field. The Hamiltonian canonical approach to the Duffin-Kemmer-Petiau theory is first developed in both the component and the matrix form. The theory is then quantized through the construction of the generating functional for the Green's functions, and the physical matrix elements of the S-matrix are proved to be relativistic invariants. The equivalence of the two theories is then proved for the matrix elements of the scattered scalar particles using the reduction formulas of Lehmann, Symanzik, and Zimmermann and for the many-photon Green's functions. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 124, No. 3, pp. 445–462, September, 2000.     

Can Soliton Attractors Exist in Realistic 3+1-D Conservative Systems?

High-energy physicists already know that stable attractors (solitons) can exist in 3+1-dimensional conservative Lagrangian systems, so long as the definition of an attractor is based on weak notions of stability and the fields admit topological charge. This paper explores the possibility of attractors in Lagrangian field theories without topological charge, using a new, stronger concept of stability—Convective quantized Asymptotic Orbital Stability (ChAOS) . Under certain conditions, ChAOS is related to additive Liapunov stability or energetic stability. Russian physicists have argued that such stability tends to require topological charge; however, this paper describes systems which avoid those arguments, and suggests how numerical examples might be constructed. Solitons have been proposed to explain the existence and nature of elementary particles within the Feynman version of quantum theory; Section 6cites this literature, as well as new possibilities for alternative versions with testable nuclear implications.     

Universality properties of Gelfand–Tsetlin patterns

A standard Gelfand–Tsetlin pattern of depth n is a configuration of particles in ${\{1,\ldots,n\} \times \mathbb{R}}$ . For each ${r \in \{1, \ldots, n\}, \{r\} \times \mathbb{R}}$ is referred to as the rth level of the pattern. A standard Gelfand–Tsetlin pattern has exactly r particles on each level r, and particles on adjacent levels satisfy an interlacing constraint. Probability distributions on the set of Gelfand–Tsetlin patterns of depth n arise naturally as distributions of eigenvalue minor processes of random Hermitian matrices of size n. We consider such probability spaces when the distribution of the matrix is unitarily invariant, prove a determinantal structure for a broad subclass, and calculate the correlation kernel. In particular we consider the case where the eigenvalues of the random matrix are fixed. This corresponds to choosing uniformly from the set of Gelfand–Tsetlin patterns whose nth level is fixed at the eigenvalues of the matrix. Fixing ${q_n \in \{1,\ldots,n\}}$ , and letting n → ∞ under the assumption that ${\frac{q_n}n \to \alpha \in (0, 1)}$ and the empirical distribution of the particles on the nth level converges weakly, the asymptotic behaviour of particles on level q _n is relevant to free probability theory. Saddle point analysis is used to identify the set in which these particles behave asymptotically like a determinantal random point field with the Sine kernel.     

A Large-Deformation Phase-Field Approach for the Modeling of Magneto-Sensitive Elastomers

Magneto-sensitive materials show magneto-mechanical coupled response and are thus of increasing interest in the recent age of smart functional materials. Ferromagnetic particles suspended in an elastomeric matrix show realignment under the influence of an external applied field, in turn causing large deformations of the substrate material. The magneto-mechanical coupling in this case is governed by the magnetic properties of the inclusion and the mechancial properties of the matrix. The magnetic phenomenon in ferromagnetic materials is governed by the formation and evolution of domains on the micro scale. A better understanding of the behavior of these particles under the influence of an external applied field is required to accurately predict the behavior of such materials. In this context it is of particular importance to model the macro scopic magneto-mechanically coupled behavior based on the micro-magnetic domain evolution. The key aspect of this work is to develop a large-deformation micro-magnetic model that can accurately capture the microscopic response of such materials. Rigorous exploitation of appropriate rate-type variational principles and consequent incremental variational principles directly give us field equations including the time evolution equation of the magnetization, which acts as the order parameter in our formulation. The theory presented here is the continuation of the work presented in [1, 7] for small deformations. A summary of magneto-mechanical theories spanning over multiple scales has been presented in [4]. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Riesz External Field Problems on the Hypersphere and Optimal Point Separation

We consider the minimal energy problem on the unit sphere ?? ^d in the Euclidean space ? ^d+1 in the presence of an external field Q, where the energy arises from the Riesz potential 1/r ^s (where r is the Euclidean distance and s is the Riesz parameter) or the logarithmic potential log(1/r). Characterization theorems of Frostman-type for the associated extremal measure, previously obtained by the last two authors, are extended to the range d ? 2 ≤ s < d ? 1. The proof uses a maximum principle for measures supported on ?? ^d . When Q is the Riesz s-potential of a signed measure and d ? 2 ≤ s < d, our results lead to explicit point-separation estimates for (Q,s)-Fekete points, which are n-point configurations minimizing the Riesz s-energy on ?? ^d with external field Q. In the hyper-singular case s > d, the short-range pair-interaction enforces well-separation even in the presence of more general external fields. As a further application, we determine the extremal and signed equilibria when the external field is due to a negative point charge outside a positively charged isolated sphere. Moreover, we provide a rigorous analysis of the three point external field problem and numerical results for the four point problem.     

PDEs for the Gaussian ensemble with external source and the Pearcey distribution

The present paper studies a Gaussian Hermitian random matrix ensemble with external source, given by a fixed diagonal matrix with two eigenvalues ±a. As a first result, the probability that the eigenvalues of the ensemble belong to an interval E satisfies a fourth‐order PDE with quartic nonlinearity; the variables are the eigenvalue a and the boundary of E. This equation enables one to find a PDE for the Pearcey distribution. The latter describes the statistics of the eigenvalues near the closure of a gap, i.e., when the support of the equilibrium measure for large‐size random matrices has a gap that can be made to close. The Gaussian Hermitian random matrix ensemble with external source, described above, has this feature. The Pearcey distribution is shown to satisfy a fourth‐order PDE with cubic nonlinearity. This also gives the PDE for the transition probability of the Pearcey process, a limiting process associated with nonintersecting Brownian motions on ℝ. © 2006 Wiley Periodicals, Inc.     

Edge scaling limits for a family of non-Hermitian random matrix ensembles

A family of random matrix ensembles interpolating between the Ginibre ensemble of n × n matrices with iid centered complex Gaussian entries and the Gaussian unitary ensemble (GUE) is considered. The asymptotic spectral distribution in these models is uniform in an ellipse in the complex plane, which collapses to an interval of the real line as the degree of non-Hermiticity diminishes. Scaling limit theorems are proven for the eigenvalue point process at the rightmost edge of the spectrum, and it is shown that a non-trivial transition occurs between Poisson and Airy point process statistics when the ratio of the axes of the supporting ellipse is of order n ^?1/3. In this regime, the family of limiting probability distributions of the maximum of the real parts of the eigenvalues interpolates between the Gumbel and Tracy–Widom distributions.     

Riesz extremal measures on the sphere for axis-supported external fields

We investigate the minimal Riesz s-energy problem for positive measures on the d-dimensional unit sphere Sd in the presence of an external field induced by a point charge, and more generally by a line charge. The model interaction is that of Riesz potentials |x−y|^−s with d−2?s<d. For a given axis-supported external field, the support and the density of the corresponding extremal measure on Sd is determined. The special case s=d−2 yields interesting phenomena, which we investigate in detail. A weak^∗ asymptotic analysis is provided as s→⁺(d−2).     

Gauge invariance and implementability of the S-operator for spin-0 and spin-12 particles in time-dependent external fields

It is proved that the classical S-operator for relativistic spin-0 and spin-$\text{1}\text{2}$ particles in time-dependent external fields is gauge invariant, and that S_+- and S_?+ are entire functions of the coupling constant in the Hilbert-Schmidt norm. As a result the Fock space S-operator exists for any real value of the coupling constant, and is gauge invariant. The external fields and the gauge function are assumed to be real-valued resp. complex-valued functions in S(R⁴).     

Higher‐Dimensional Coulomb Gases and Renormalized Energy Functionals

We consider a classical system of n charged particles in an external confining potential in any dimension d ≥ 2. The particles interact via pairwise repulsive Coulomb forces and the coupling parameter is of order n^?1 (mean‐field scaling). By a suitable splitting of the Hamiltonian, we extract the next‐to‐leading‐order term in the ground state energy beyond the mean‐field limit. We show that this next order term, which characterizes the fluctuations of the system, is governed by a new “renormalized energy” functional providing a way to compute the total Coulomb energy of a jellium (i.e., an infinite set of point charges screened by a uniform neutralizing background) in any dimension. The renormalization that cuts out the infinite part of the energy is achieved by smearing out the point charges at a small scale, as in Onsager's lemma. We obtain consequences for the statistical mechanics of the Coulomb gas: next‐to‐leading‐order asymptotic expansion of the free energy or partition function, characterizations of the Gibbs measures, estimates on the local charge fluctuations, and factorization estimates for reduced densities. This extends results of Sandier and Serfaty to dimension higher than 2 by an alternative approach. © 2016 Wiley Periodicals, Inc.     

The Limiting Spectral Measure for Ensembles of Symmetric Block Circulant Matrices

Given an ensemble of N×N random matrices, a natural question to ask is whether or not the empirical spectral measures of typical matrices converge to a limiting spectral measure as N→∞. While this has been proved for many thin patterned ensembles sitting inside all real symmetric matrices, frequently there is no nice closed form expression for the limiting measure. Further, current theorems provide few pictures of transitions between ensembles. We consider the ensemble of symmetric m-block circulant matrices with entries i.i.d.r.v. These matrices have toroidal diagonals periodic of period m. We view m as a “dial” we can “turn” from the thin ensemble of symmetric circulant matrices, whose limiting eigenvalue density is a Gaussian, to all real symmetric matrices, whose limiting eigenvalue density is a semi-circle. The limiting eigenvalue densities f _m show a visually stunning convergence to the semi-circle as m→∞, which we prove. In contrast to most studies of patterned matrix ensembles, our paper gives explicit closed form expressions for the densities. We prove that f _m is the product of a Gaussian and a certain even polynomial of degree 2m?2; the formula is the same as that for the m×m Gaussian Unitary Ensemble (GUE). The proof is by derivation of the moments from the eigenvalue trace formula. The new feature, which allows us to obtain closed form expressions, is converting the central combinatorial problem in the moment calculation into an equivalent counting problem in algebraic topology. We end with a generalization of the m-block circulant pattern, dropping the assumption that the m random variables be distinct. We prove that the limiting spectral distribution exists and is determined by the pattern of the independent elements within an m-period, depending not only on the frequency at which each element appears, but also on the way the elements are arranged.     

Convergence of Metropolis-type algorithms for a large canonical ensemble

In this paper, we study the convergence of Metropolis-type algorithms used in modeling statistical systems with a fluctuating number of particles located in a finite volume. We justify the use of Metropolis algorithms for a particular class of such statistical systems. We prove a theorem on the geometric ergodicity of the Markov process modeling the behavior of an ensemble with a fluctuating number of particles in a finite volume whose interaction is described by a potential bounded below and decreasing according to the law r ^?3?α, α ≥ 0, as r → 0.     

Spectra of Large Random Trees

We analyze the eigenvalues of the adjacency matrices of a wide variety of random trees. Using general, broadly applicable arguments based on the interlacing inequalities for the eigenvalues of a principal submatrix of a Hermitian matrix and a suitable notion of local weak convergence for an ensemble of random trees that we call probability fringe convergence, we show that the empirical spectral distributions for many random tree models converge to a deterministic (model-dependent) limit as the number of vertices goes to infinity. Moreover, the masses assigned by the empirical spectral distributions to individual points also converge in distribution to constants. We conclude for ensembles such as the linear preferential attachment models, random recursive trees, and the uniform random trees that the limiting spectral distribution has a set of atoms that is dense in the real line. We obtain lower bounds on the mass assigned to zero by the empirical spectral measures via the connection between the number of zero eigenvalues of the adjacency matrix of a tree and the cardinality of a maximal matching on the tree. In particular, we employ a simplified version of an algorithm due to Karp and Sipser to construct maximal matchings and understand their properties. Moreover, we show that the total weight of a weighted matching is asymptotically equivalent to a constant multiple of the number of vertices when the edge weights are independent, identically distributed, nonnegative random variables with finite expected value, thereby significantly extending a result obtained by Aldous and Steele in the special case of uniform random trees. We greatly generalize a celebrated result obtained by Schwenk for the uniform random trees by showing that if any ensemble converges in the probability fringe sense and a very mild further condition holds, then, with probability converging to one, the spectrum of a realization is shared by at least one other (nonisomorphic) tree. For the linear preferential attachment model with parameter a>?1, we show that for any fixed k, the k largest eigenvalues jointly converge in distribution to a nontrivial limit when rescaled by $n^{1/2\gamma_{a}}$ , where ?? _a =a+2 is the Malthusian rate of growth parameter for an associated continuous-time branching process.     

Classification of (v, 3)-configurations

A (v, 3)-configuration is a nondegenerate matrix of dimension v over the field GF(2) considered up to permutation of rows and columns and containing exactly three 1’s in the rows and columns, while the inverse matrix has also exactly three 1’s in the rows and columns. It is proved that, for each even v ≥ 4, there is only one indecomposable (v, 3)-configuration, while, for odd v, there are no such configurations, the only exception being the unique (5, 3)-configuration.     

Existence of an A-optimal model for a regression experiment

The existence of certain A-optimal models for a regression experiment is equivalent to the solvability of certain matrix equations. It is proved that for any m × m matrix V (over the real field), there exists an orthogonal matrix Q such that QVQ′ has equal diagonal elements. This result is used to solve the equations mentioned above. The connection of this result to Hadamard matrices is discussed.     

Regularity of the solution and well-posedness of a mixed problem for an elliptic system with quadratic nonlinearity in gradients

A boundary value problem for an elliptic system of equations is studied that arises in the analysis of a new hydrodynamic model describing charge transport in a planar semiconductor MESFET (metal semiconductor field effect transistor). The problem has a number of features, specifically, the equations of the system involve squared components of the gradients of the unknown functions; the boundary conditions are of a mixed character, i.e., Dirichlet and Neumann conditions are set on different portions of the boundary; and the boundary of the domain is a nonsmooth curve, namely, a rectangle. Under a certain optimal condition, the C ^1,α-regularity of a weakened solution of the problem is justified and its existence is proved, while its uniqueness is shown under additional constraints. The results are used to justify the stabilization method as applied to finding approximate stationary solutions of the hydrodynamic model.     

Simultaneous local exact controllability of 1D bilinear Schrödinger equations

We consider N independent quantum particles, in an infinite square potential well coupled to an external laser field. These particles are modelled by a system of linear Schrödinger equations on a bounded interval. This is a bilinear control system in which the state is the N  -tuple of wave functions. The control is the real amplitude of the laser field. For N=1$N=1$, Beauchard and Laurent proved local exact controllability around the ground state in arbitrary time. We prove, under an extra generic assumption, that their result does not hold in small time if N?2$N?2$. Still, for N=2$N=2$, we prove that local controllability holds either in arbitrary time up to a global phase or exactly up to a global delay. This is proved using Coron's return method. We also prove that for N?3$N?3$, local controllability does not hold in small time even up to a global phase. Finally, for N=3$N=3$, we prove that local controllability holds up to a global phase and a global delay.