Limit Of Infinite Dimensions Research Articles

Transverse forces and glassy liquids in infinite dimensions.

We explore the dynamics of a simple liquid whose particles, in addition to standard potential-based interactions, are also subjected to transverse forces preserving the Boltzmann distribution. We derive the effective dynamics of one and two tracer particles in the infinite-dimensional limit. We determine the amount of acceleration of the dynamics caused by the transverse forces, in particular in the vicinity of the glass transition. We analyze the emergence and evolution of odd transport phenomena induced by the transverse forces.

Frustrated magnets in the limit of infinite dimensions: Dynamics and disorder-free glass transition

Frustrated magnets in the limit of infinite dimensions: Dynamics and disorder-free glass transition

Glass transition of quantum hard spheres in high dimensions.

We study the equilibrium thermodynamics of quantum hard spheres in the infinite-dimensional limit, determining the boundary between liquid and glass phases in the temperature-density plane by means of the Franz-Parisi potential. We find that as the temperature decreases from high values, the effective radius of the spheres is enhanced by a multiple of the thermal de Broglie wavelength, thus increasing the effective filling fraction and decreasing the critical density for the glass phase. Numerical calculations show that the critical density continues to decrease monotonically as the temperature decreases further, suggesting that the system will form a glass at sufficiently low temperatures for any density. The methods used in this paper can be extended to more general potentials, and also to other transitions such as the Kauzman/Replica Symmetry Breaking (RSB) transition, the Gardner transition, and potentially even jamming.

Chaos in Stochastic 2d Galerkin-Navier–Stokes

We prove that all Galerkin truncations of the 2d stochastic Navier–Stokes equations in vorticity form on any rectangular torus subjected to hypoelliptic, additive stochastic forcing are chaotic at sufficiently small viscosity, provided the frequency truncation satisfies N≥392\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N\\ge 392$$\\end{document}. By “chaotic” we mean having a strictly positive Lyapunov exponent, i.e. almost-sure asymptotic exponential growth of the derivative with respect to generic initial conditions. A sufficient condition for such results was derived in previous joint work with Alex Blumenthal which reduces the question to the non-degeneracy of a matrix Lie algebra implying Hörmander’s condition for the Markov process lifted to the sphere bundle (projective hypoellipticity). The purpose of this work is to reformulate this condition to be more amenable for Galerkin truncations of PDEs and then to verify this condition using (a) a reduction to genericity properties of a diagonal sub-algebra inspired by the root space decomposition of semi-simple Lie algebras and (b) computational algebraic geometry executed by Maple in exact rational arithmetic. Note that even though we use a computer assisted proof, the result is valid for all aspect ratios and all sufficiently high dimensional truncations; in fact, certain steps simplify in the formal infinite dimensional limit.

Low-frequency spin dynamics in diluted magnetic semiconductors

The low-frequency dynamical magnetic properties in paramagnetic diluted magnetic semiconductors are addressed in the framework of the dynamical mean-field theory applied for the Kondo lattice model. In the infinite-dimensional limit, a set of self-consistent equations is derived so the single-particle Green's function and its self-energy can be evaluated numerically. In terms of the Green's function and self-energies, the local dynamical spin susceptibility function and then the spin-relaxation rate are explicitly expressed based on the Baym-Kadanoff approach. It is found that the spin fluctuations become dominated, indicated by the sharp peak appearing at the low frequency of the spin dynamical susceptibility function in the case of large magnetic coupling and temperature close to the paramagnetic-ferromagnetic transition point. The low-frequency spin dynamic in the systems is also addressed in the signatures of the spin-relaxation process. In the case of large temperature and small magnetic coupling, the spin-relaxation rate releases the scenario of the Korringa process specifying the weak correlation systems likely normal metals. Otherwise, i.e., at small temperature and large magnetic coupling, we find exponential behavior of the spin-relaxation rate versus temperature. Moreover, at a temperature approaching the paramagnetic-ferromagnetic transition point, one finds sharp suppression of the spin-relaxation rate or speeding up of the spin-relaxation time. These scenarios are attributed to the appearance of the magnetic coherence bound state or the spin clusters in diluted magnetic semiconductors due to the strongly magnetic correlations.

Hard-sphere jamming through the lens of linear optimization.

The jamming transition is ubiquitous. It is present in granular matter, foams, colloids, structural glasses, and many other systems. Yet, it defines a critical point whose properties still need to be fully understood. Recently, a major breakthrough came about when the replica formalism was extended to build a mean-field theory that provides an exact description of the jamming transition of spherical particles in the infinite-dimensional limit. While such theory explains the jamming critical behavior of both soft and hard spheres, investigating the transition in finite-dimensional systems poses very difficult and different problems, in particular from the numerical point of view. Soft particles are modeled by continuous potentials; thus, their jamming point can be reached through efficient energy minimization algorithms. In contrast, the latter methods are inapplicable to hard-sphere (HS) systems since the interaction energy among the particles is always zero by construction. To overcome these difficulties, here we recast the jamming of hard spheres as a constrained optimization problem and introduce the CALiPPSO algorithm, capable of readily producing jammed HS packings without including any effective potential. This algorithm brings a HS configuration of arbitrary dimensions to its jamming point by solving a chain of linear optimization problems. We show that there is a strict correspondence between the force balance conditions of jammed packings and the properties of the optimal solutions of CALiPPSO, whence we prove analytically that our packings are always isostatic and in mechanical equilibrium. Furthermore, using extensive numerical simulations, we show that our algorithm is able to probe the complex structure of the free-energy landscape, finding qualitative agreement with mean-field predictions. We also characterize the algorithmic complexity of CALiPPSO and provide an open-source implementation of it.

Contracting differential equations in weighted Banach spaces

Geodesic contraction in vector-valued differential equations is readily verified by linearized operators which are uniformly negative-definite in the Riemannian metric. In the infinite-dimensional setting, however, such analysis is generally restricted to norm-contracting systems. We develop a generalization of geodesic contraction rates to Banach spaces using a smoothly-weighted semi-inner product structure on tangent spaces. We show that negative contraction rates in bijectively weighted spaces imply asymptotic norm-contraction, and apply recent results on asymptotic contractions in Banach spaces to establish the existence of fixed points. We show that contraction in surjectively weighted spaces verifies non-equilibrium asymptotic properties, such as convergence to finite- and infinite-dimensional subspaces, submanifolds, limit cycles, and phase-locking phenomena. We use contraction rates in weighted Sobolev spaces to establish existence and continuous data dependence in nonlinear PDEs, and pose a method for constructing weak solutions using vanishing one-sided Lipschitz approximations. We discuss applications to control and order reduction of PDEs.

Extraction of ergotropy: free energy bound and application to open cycle engines

The second law of thermodynamics uses change in free energy of macroscopic systems to set a bound on performed work. Ergotropy plays a similar role in microscopic scenarios, and is defined as the maximum amount of energy that can be extracted from a system by a unitary operation. In this analysis, we quantify how much ergotropy can be induced on a system as a result of system's interaction with a thermal bath, with a perspective of using it as a source of work performed by microscopic machines. We provide the fundamental bound on the amount of ergotropy which can be extracted from environment in this way. The bound is expressed in terms of the non-equilibrium free energy difference and can be saturated in the limit of infinite dimension of the system's Hamiltonian. The ergotropy extraction process leading to this saturation is numerically analyzed for finite dimensional systems. Furthermore, we apply the idea of extraction of ergotropy from environment in a design of a new class of stroke heat engines, which we label open-cycle engines. Efficiency and work production of these machines can be completely optimized for systems of dimensions 2 and 3, and numerical analysis is provided for higher dimensions.

On the Infinite Dimension Limit of Invariant Measures and Solutions of Zeitlin’s 2D Euler Equations

In this work we consider a finite dimensional approximation for the 2D Euler equations on the sphere {mathbb {S}}^2, proposed by V. Zeitlin, and show their convergence towards a solution to Euler equations with marginals distributed as the enstrophy measure. The method relies on nontrivial computations on the structure constants of {mathbb {S}}^2, that appear to be new. In the last section we discuss the problem of extending our results to Gibbsian measures associated with higher Casimirs.

On the empirical spectral distribution for certain models related to sample covariance matrices with different correlations

Given [Formula: see text], we study two classes of large random matrices of the form [Formula: see text] where for every [Formula: see text], [Formula: see text] are iid copies of a random variable [Formula: see text], [Formula: see text], [Formula: see text] are two (not necessarily independent) sets of independent random vectors having different covariance matrices and generating well concentrated bilinear forms. We consider two main asymptotic regimes as [Formula: see text]: a standard one, where [Formula: see text], and a slightly modified one, where [Formula: see text] and [Formula: see text] while [Formula: see text] for some [Formula: see text]. Assuming that vectors [Formula: see text] and [Formula: see text] are normalized and isotropic “in average”, we prove the convergence in probability of the empirical spectral distributions of [Formula: see text] and [Formula: see text] to a version of the Marchenko–Pastur law and the so-called effective medium spectral distribution, correspondingly. In particular, choosing normalized Rademacher random variables as [Formula: see text], in the modified regime one can get a shifted semicircle and semicircle laws. We also apply our results to the certain classes of matrices having block structures, which were studied in [G. M. Cicuta, J. Krausser, R. Milkus and A. Zaccone, Unifying model for random matrix theory in arbitrary space dimensions, Phys. Rev. E 97(3) (2018) 032113, MR3789138; M. Pernici and G. M. Cicuta, Proof of a conjecture on the infinite dimension limit of a unifying model for random matrix theory, J. Stat. Phys. 175 (2) (2019) 384–401, MR3968860].

Global well-posedness and long-time behavior of the fractional NLS

In this paper, our discussion mainly focuses on equations with energy supercritical nonlinearities. We establish probabilistic global well-posedness (GWP) results for the cubic Schrödinger equation with any fractional power of the Laplacian in all dimensions. We consider both low and high regularities in the radial setting, in dimension ge 2. In the high regularity result, an Inviscid - Infinite dimensional (IID) limit is employed while in the low regularity global well-posedness result, we make use of the Skorokhod representation theorem. The IID limit is presented in details as an independent approach that applies to a wide range of Hamiltonian PDEs. Moreover we discuss the adaptation to the periodic settings, in any dimension, for smooth regularities.

Mean-field dynamics of infinite-dimensional particle systems: global shear versus random local forcing

In infinite dimensions, many-body systems of pairwise interacting particles provide exact analytical benchmarks for the features of amorphous materials, such as the stress–strain curve of glasses under quasistatic shear. Here, instead of global shear, we consider an alternative driving protocol, as recently introduced by Morse et al 2020 (arXiv:2009.07706), which consists of randomly assigning a constant local displacement on each particle, with a finite spatial correlation length. We show that, in the infinite-dimensional limit, the mean-field dynamics under such a random forcing are strictly equivalent to those under global shear, upon a simple rescaling of the accumulated strain. Moreover, the scaling factor is essentially given by the variance of the relative local displacements of interacting pairs of particles, which encodes the presence of a finite spatial correlation. In this framework, global shear is simply a special case of a much broader family of local forcing, which can be explored by tuning its spatial correlations. We discuss the specific implications for the quasistatic driving of glasses—initially prepared at a replica-symmetric equilibrium—and how the corresponding ‘stress–strain’-like curves and elastic moduli can be rescaled onto their quasistatic-shear counterparts. These results hint at a unifying framework for establishing rigorous analogies, at the mean-field level, between different driven disordered systems.

Effect of spatial dimension on a model of fluid turbulence

A numerical study of the $d$ -dimensional eddy damped quasi-normal Markovian equations is performed to investigate the dependence on spatial dimension of homogeneous isotropic fluid turbulence. Relationships between structure functions and energy and transfer spectra are derived for the $d$ -dimensional case. Additionally, an equation for the $d$ -dimensional enstrophy analogue is derived and related to the velocity derivative skewness. Comparisons are made to recent four-dimensional direct numerical simulation results. Measured energy spectra show a magnified bottleneck effect which grows with dimension whilst transfer spectra show a varying peak in the nonlinear energy transfer as the dimension is increased. These results are consistent with an increased forward energy transfer at higher dimensions, further evidenced by measurements of a larger asymptotic dissipation rate with growing dimension. The enstrophy production term, related to the velocity derivative skewness, is seen to reach a maximum at around five dimensions and may reach zero in the limit of infinite dimensions, raising interesting questions about the nature of turbulence in this limit.

Mean-Field Predictions of Scaling Prefactors Match Low-Dimensional Jammed Packings

No known analytic framework precisely explains all the phenomena observed in jamming. The replica theory for glasses and jamming is a mean-field theory which attempts to do so by working in the limit of infinite dimensions, such that correlations between neighbors are negligible. As such, results from this mean-field theory are not guaranteed to be observed in finite dimensions. However, many results in mean field for jamming have been shown to be exact or nearly exact in low dimensions. This suggests that the infinite dimensional limit is not necessary to obtain these results. In this Letter, we perform precision measurements of jamming scaling relationships between pressure, excess packing fraction, and number of excess contacts from dimensions 2-10 in order to extract the prefactors to these scalings. While these prefactors should be highly sensitive to finite dimensional corrections, we find the mean-field predictions for these prefactors to be exact in low dimensions. Thus the mean-field approximation is not necessary for deriving these prefactors. We present an exact, first-principles derivation for one, leaving the other as an open question. Our results suggest that mean-field theories of critical phenomena may compute more for d≥d_{u} than has been previously appreciated.

Surface growth on tree-like lattices and the upper critical dimension of the KPZ class

Aiming to investigate the upper critical dimension, d u , of the KPZ class, in Saberi A. A., EPL, 103 (2013) 10005, some growth models were numerically analyzed using Cayley trees (CTs) as substrates, as a way to access their behavior in the infinite-dimensional limit, and some unexpected results were reported: logarithmic roughness scaling, differing for EW and KPZ models (indicating that even at the KPZ nonlinearity is still relevant); beyond asymptotically rough EW surfaces above the upper critical dimension of the EW class. Motivated by these strange findings, I revisit these growth models here to show that such results are simple consequences of boundary effects, inherent to systems defined on CTs. In fact, I demonstrate that the anomalous boundary of the CT leads the growing surfaces to develop curved shapes, which explains the strange behaviors previously found for these systems, once the global “roughness” was analyzed for non-flat surfaces in the study above. Importantly, by measuring the height fluctuations at the central site of the CT, which can be seen as an approximation for the Bethe lattice, smooth surfaces are found for both EW and KPZ classes, consistently with the behavior expected for growing systems in dimensions . Interesting features of the 1-point height fluctuations, such as the possibility of non-saturation in the steady state regime, are also discussed for substrates in general.

Marginally jammed states of hard disks in a one-dimensional channel.

We have studied a class of marginally jammed states in a system of hard disks confined in a narrow channel-a quasi-one-dimensional system-whose exponents are not those predicted by theories valid in the infinite dimensional limit. The exponent γ which describes the distribution of small gaps takes the value 1 rather than the infinite dimensional value 0.41269⋯. Our work shows that there exist jammed states not found within the tiling approach of Ashwin and Bowles. The most dense of these marginal states is an unusual state of matter that is asymptotically crystalline.

A numerical approach to Kolmogorov equation in high dimension based on Gaussian analysis

For Kolmogorov equations associated to finite dimensional stochastic differential equations (SDEs) in high dimension, a numerical method alternative to Monte Carlo simulations is proposed. The structure of the SDE is inspired by stochastic Partial Differential Equations (SPDE) and thus contains an underlying Gaussian process which is the key of the algorithm. A series development of the solution in terms of iterated integrals of the Gaussian process is given, it is proved to converge - also in the infinite dimensional limit - and it is numerically tested in a number of examples.

Extending the Gutzwiller approximation to intersite interactions

We develop an extension of the Gutzwiller approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping, and Coulomb-assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multiband Hubbard models in the presence of only local interactions. This is accomplished by combining the $1/z$ expansion---where $z$ is the coordination number, and only the leading-order terms contribute in the limit of infinite dimensions---with a ${P}_{R}^{\ifmmode\dagger\else\textdagger\fi{}}{P}_{R}\ensuremath{-}I$ expansion, where ${P}_{R}$ is the Gutzwiller projector on the site $R$. The method is conveniently formulated in terms of a Gutzwiller Lagrange function. We apply our theory to the extended single-band Hubbard model. Similarly to the usual Brinkman-Rice mechanism, we find a Mott transition. A valence skipping transition is observed, where the occupation of the empty and doubly occupied states for the Gutzwiller wave function is enhanced with respect to the uncorrelated Slater determinant wave function.

Random geometric graphs in high dimension.

Many machine learning algorithms used for dimensional reduction and manifold learning leverage on the computation of the nearest neighbors to each point of a data set to perform their tasks. These proximity relations define a so-called geometric graph, where two nodes are linked if they are sufficiently close to each other. Random geometric graphs, where the positions of nodes are randomly generated in a subset of R^{d}, offer a null model to study typical properties of data sets and of machine learning algorithms. Up to now, most of the literature focused on the characterization of low-dimensional random geometric graphs whereas typical data sets of interest in machine learning live in high-dimensional spaces (d≫10^{2}). In this work, we consider the infinite dimensions limit of hard and soft random geometric graphs and we show how to compute the average number of subgraphs of given finite size k, e.g., the average number of k cliques. This analysis highlights that local observables display different behaviors depending on the chosen ensemble: soft random geometric graphs with continuous activation functions converge to the naive infinite-dimensional limit provided by Erdös-Rényi graphs, whereas hard random geometric graphs can show systematic deviations from it. We present numerical evidence that our analytical results, exact in infinite dimensions, provide a good approximation also for dimension d≳10.

Edge of spiked beta ensembles, stochastic Airy semigroups and reflected Brownian motions

Nous accedons a l’extremite du spectre des ensembles beta gaussiens avec perturbation de rang un par l’entremise de grandes puissances des matrices tridiagonales qui y sont associees. Pour les valeurs traditionnelles $\beta =1,2,4$, ceci correspond a l’etude des fluctuations des valeurs propres maximales des matrices aleatoires GOE/GUE/GSE assujetties a une perturbation additive de rang un. En dimensions infinies, nos resultats nous menent vers une famille de semi-groupes de type Feynman–Kac qui, dans un cas extreme, correspond au stochastic Airy semigroup introduit par Gorin et Shkolnikov (Ann. Probab. 46 (2018) 2287–2344). De plus, nos resultats ont pour corollaire des formules de Feynman–Kac pour les spiked stochastic Airy operators introduits par Bloemendal et Virag (Probab. Theory Related Fields 156 (2013) 795–825). Ces formules sont exprimees a l’aide de certaines fonctionnelles du mouvement brownien reflechi et de ses temps locaux. Ce faisant, les operateurs en question peuvent etre etudies a l’aide du calcul stochastique. Nous obtenons un premier resultat dans cette lignee en demontrant une nouvelle identite decrivant la distribution du mouvement brownien reflechi ayant ete conditionne sur son temps local a zero. La principale innovation de notre demonstration consiste en la preuve d’un nouveau resultat sur l’approximation forte du mouvement brownien reflechi et de son temps local par une marche aleatoire non negative en utilisant la methode de reflexion de Skorokhod.