Many-particle Systems Research Articles

Diffusion enhancement and autoparametric resonance.

The possibility of an autoparametric resonance in an isolated many-particle system induces a specific behavior of the particles in the presence of thermal noise. In particular, the variance associated with a resonant mode, and consequently that of the associated particles, is strongly increased compared to what it would have in the absence of parametric resonance. In this paper we consider a dimer submitted to a periodic potential for which there are only two modes, the center of mass motion and the internal vibration mode. This is the simplest system which is dynamically rich enough to exhibit an autoparametric excitation of the internal vibrations by the center of mass motion. The consequences of this autoparametric excitation on the particles diffusion will be discussed according to the stiffness of the interaction and to the initial energy of the dimer, the relevant parameters which characterize this dynamics.

A variational approach to deriving Hamiltonian functions under non-holonomic constraints — the case of the Gaussian thermostats

The Hamiltonian function (H), if it exists, provides a powerful platform for studying the dynamics of constrained many-particle systems. For systems under holonomic constraints, H can be derived from the principle of least action (PLA) formulated in terms of a suitable set of generalized coordinates. In molecular dynamics (MD) simulations, where the system temperature often prescribes a non-holonomic constraint, the challenge is to find the system’s Hamiltonian. The Gaussian Isokinetic (GIK) thermostat, the first deterministic physics-preserving thermostat in MD, was derived from the principle of least constraint (PLC) in the early 1980s: its H, interestingly, was discovered almost 15 years later. A mathematically rigorous connection of the GIK Hamiltonian with PLA, however, has been missing until now. In this article, we adopt a variational approach to establish the equivalence of PLC with PLA for a class of non-holonomically constrained many-particle systems, and present the general formulation of the Euler–Lagrange equations of motion through time rescaling. Making use of the particular forms of the Lagrange multiplier obtained from the principle of least constraint, we derive H for the GIK and the Gaussian Isoenergetic thermostats.

Controlling inter-particle distances in crowds of motile, cognitive, active particles

Distance control in many-particle systems is a fundamental problem in nature. This becomes particularly relevant in systems of active agents, which can sense their environment and react by adjusting their direction of motion. We employ agent-based simulations to investigate the complex interplay between agent activity, characterized by Péclet number Pe\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox{Pe}$$\\end{document}, reorientation maneuverability Ω\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega$$\\end{document}, vision angle θ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta$$\\end{document} and vision range R0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$R_0$$\\end{document}, and agent density, which determines agent distancing and dynamics. We focus on semi-dense crowds, where the vision range is much larger than the particle size. The minimal distance to the nearest neighbors, exposure time, and persistence of orientation direction are analyzed to characterize the behavior. With increasing particle speed at fixed maneuverability, particles approach each other more closely, and exhibit shorter exposure times. The temporal persistence of motion decreases with increasing Pe\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox{Pe}$$\\end{document}, reflecting the impact of activity and maneuverability on direction changes. For a vision angle θ=π/4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta =\\pi /4$$\\end{document}, we observe the emergence of flocking aggregates with a band-like structure, somewhat reminiscent of the bands in the Vicsek model. Additionally, for vision angles θ≥π/2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta \\ge \\pi /2$$\\end{document}, several quantities are found to display a universal scaling behavior with scaling variable Pe3/2/Ω\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox{Pe}^{3/2}/\\Omega$$\\end{document}. Our results are in good agreement with recent experiments of pedestrians in confined spaces.

Classical and Quantum Mechanical Models of Many-Particle Systems

The workshop focused on the collective behavior of many-particle systems in various application fields: physics (gas dynamics, plasmas, quantum mechanics), mathematical biology (cell mobility, evolution of trait-structured species), and social sciences (wealth distribution). This includes famous models such as the Boltzmann equation of gas dynamics, Vlasov equation for plasmas, Fokker–Planck equations, Smoluchowski and related equations, Keller–Segel system of chemotaxis.

Optimized finite-time work protocols for the Higgs RNA model with external force.

The Higgs RNA model with an added term for a coupling to an external force is studied in regard to finite-time force-driving protocols with a minimal-work requirement. In this paper, RNA sequences which at low temperature exhibit hairpins are considered, which are often cited as typical template systems in stochastic thermodynamics. The optimized work protocols for this glassy many-particle system are determined numerically using the parallel tempering method. The protocols show distinct jumps at the beginning and end, which have been observed for single-particle systems and are proven to be optimal in the fast protocol limit generally. Optimality seems to be achieved by staying close to the equilibrium unfolding transition point, in agreement with experimental and theoretical observations. The change of work distributions, compared to those resulting from a naive linear driving protocol, are discussed generally and in terms of free energy estimation as well as the effect of optimized protocols on rare work process starting conditions.

Optimal evaluation of symmetry-adapted n-correlations via recursive contraction of sparse symmetric tensors

Abstract We present a comprehensive analysis of an algorithm for evaluating high-dimensional polynomials that are invariant (or equi-variant) under permutations and rotations. This task arises in the evaluation linear models as well as equivariant neural network models of many-particle systems. The theoretical bottleneck is the contraction of a high-dimensional symmetric and sparse tensor with a specific sparsity pattern that is directly related to the symmetries imposed on the polynomial. The sparsity of this tensor makes it challenging to construct a highly efficient evaluation scheme. The references [10, 11] introduced a recursive evaluation strategy that relied on a number of heuristics, but performed well in tests. In the present work, we propose an explicit construction of such a recursive evaluation strategy and show that it is in fact optimal in the limit of infinite polynomial degree.

Asymptotic Antipodal Solutions as the Limit of Elliptic Relative Equilibria for the Two- and n-Body Problems in the Two-Dimensional Conformal Sphere

We consider the two- and n-body problems on the two-dimensional conformal sphere MR2, with a radius R>0. We employ an alternative potential free of singularities at antipodal points. We study the limit of relative equilibria under the SO(2) symmetry; we examine the specific conditions under which a pair of positive-mass particles, situated at antipodal points, can maintain a state of relative equilibrium as they traverse along a geodesic. It is identified that, under an appropriate radius–mass relationship, these particles experience an unrestricted and free movement in alignment with the geodesic of the canonical Killing vector field in MR2. An even number of bodies with pairwise conjugated positions, arranged in a regular n-gon, all with the same mass m, move freely on a geodesic with suitable velocities, where this geodesic motion behaves like a relative equilibrium. Also, a center of mass formula is included. A relation is found for the relative equilibrium in the two-body problem in the sphere similar to the Snell law.

Entanglement Phase Transitions in Non-Hermitian Kitaev Chains.

The intricate interplay between unitary evolution and projective measurements could induce entanglement phase transitions in the nonequilibrium dynamics of quantum many-particle systems. In this work, we uncover loss-induced entanglement transitions in non-Hermitian topological superconductors. In prototypical Kitaev chains with onsite particle losses and varying hopping and pairing ranges, the bipartite entanglement entropy of steady states is found to scale logarithmically versus the system size in topologically nontrivial phases and become independent of the system size in the trivial phase. Notably, the scaling coefficients of log-law entangled phases are distinguishable when the underlying system resides in different topological phases. Log-law to log-law and log-law to area-law entanglement phase transitions are further identified when the system switches between different topological phases and goes from a topologically nontrivial to a trivial phase, respectively. These findings not only establish the relationships among spectral, topological and entanglement properties in a class of non-Hermitian topological superconductors but also provide an efficient means to dynamically reveal their distinctive topological features.

HODLR3D: hierarchical matrices for N-body problems in three dimensions

HODLR3D: hierarchical matrices for N-body problems in three dimensions

Perturbing masses: A study of centered co-circular central configurations in power-law n-body problems

This research investigates centered co-circular central configurations in the general power-law potential n-body problem. Firstly, there are no centered co-circular central configurations exist when precisely n−2 of the masses are equal; secondly, unless all masses are equal, no such configurations exist when masses can be divided into two sets of equal masses; thirdly, no such configurations exist when masses is a mix of equal masses and some heavier (lighter) ones, arranged such that the heavier (lighter) ones are non-adjacent.

Virial equation of state for a granular system.

The equation of state for an ideal gas is simple, which is . In the case of imperfect gases where mutual interactions among the constituents are important, pressure P can be expressed as the series expansion of density n with appropriate coefficients, known as virial coefficients . In this paper, we have obtained the first four virial coefficients for a model interaction potential using multidimensional Monte-Carlo integration and importance sampling methods. Next, we perform molecular dynamics simulations with the same for a many-particle system to obtain P as a function of T and n. We compare our numerical data with the virial equation of state.

Effectiveness of educational interventions in minimizing reactance in traffic accident prevention

IntroductionClassical mechanics concepts are vital for understanding traffic safety, encompassing speed, acceleration, stopping distance, impact speed, deformation, kinetic energy, and their effects on vehicle occupants. Statistical physics principles are also employed to study multiple vehicles as a many-particle system. This study examines the effectiveness of using a damaged car and structured instruction in road safety initiatives.MethodsParticipants engaged with a damaged car and discussed deformation caused by excessive speed. They participated in role simulations, envisioning and solving risky traffic scenarios, while considering the potential boomerang effect resulting from excessive fear appeals within the campaign.ResultsCombining emotional appeals with solution-focused, action-oriented, and self-esteem-enhancing interventions is essential to mitigate the boomerang effect and encourage safety-compliant behavior. Educational institutions must implement tailored follow-ups to emotionally charged prevention campaigns. Using Design-Based Research, we developed effective methods for knowledge acquisition and transfer, drawing from cognitive behavioral therapy and established traffic safety programs.DiscussionThis publication explores non-standardized social training within focus groups, focusing on negative emotional states following an emotionalizing campaign similar to “Crash Kurs NRW,” targeting upper-level courses in North Rhine-Westphalian schools. Methods from cognitive behavioral therapy and social learning initiated, observed, and documented changes in participants’ norms, values, and attitudes, discussed herein. Additionally, we investigated whether the presence of police officers in prevention campaigns influences reactive behavior. Comparing pre- and post-surveys confirmed reactance after emotionalizing campaigns, and our structured post-survey effectively influenced stage-event-induced behavior. Surprisingly, police officers’ presence did not significantly impact reactive behavior. Furthermore, our follow-up tool extends its applicability to similar emotionalizing prevention campaigns, broadening its potential reach. Our future research will delve into students’ physics literacy during traffic space analysis of accident hotspots, advancing traffic safety initiatives.

Designer pair statistics of disordered many-particle systems with novel properties.

The knowledge of exact analytical functional forms for the pair correlation function g2(r) and its corresponding structure factor S(k) of disordered many-particle systems is limited. For fundamental and practical reasons, it is highly desirable to add to the existing database of analytical functional forms for such pair statistics. Here, we design a plethora of such pair functions in direct and Fourier spaces across the first three Euclidean space dimensions that are realizable by diverse many-particle systems with varying degrees of correlated disorder across length scales, spanning a wide spectrum of hyperuniform, typical nonhyperuniform, and antihyperuniform ones. This is accomplished by utilizing an efficient inverse algorithm that determines equilibrium states with up to pair interactions at positive temperatures that precisely match targeted forms for both g2(r) and S(k). Among other results, we realize an example with the strongest hyperuniform property among known positive-temperature equilibrium states, critical-point systems (implying unusual 1D systems with phase transitions) that are not in the Ising universality class, systems that attain self-similar pair statistics under Fourier transformation, and an experimentally feasible polymer model. We show that our pair functions enable one to achieve many-particle systems with a wide range of translational order and self-diffusion coefficients D, which are inversely related to one another. One can design other realizable pair statistics via linear combinations of our functions or by applying our inverse procedure to other desirable functional forms. Our approach facilitates the inverse design of materials with desirable physical and chemical properties by tuning their pair statistics.

Energy is not a convex function of particle number for r-k interparticle potentials with k > log34.

The energy of a many-particle system is not convex with respect to particle number for r-k interparticle repulsion potentials if k > log34 ≈ 1.262. With such potentials, some finite electronic systems have ionization potentials that are less than the electron affinity: they have negative band gap (chemical hardness). Although the energy may be a convex function of the number of electrons (for which k = 1), it suggests that finding an analytic proof of convexity will be very difficult. The bound on k is postulated to be tight. An apparent signature of non-convex behavior is that the Dyson orbital corresponding to the lowest-energy mode of electron attachment has a vanishingly small amplitude.

Explicit near-symplectic integrators for post-Newtonian Hamiltonian systems

Explicit symplectic integrators are powerful and widely used for Hamiltonian systems. However, once the post-Newtonian (PN) effect is considered to provide more precise modeling for the N-body problem, explicit symplectic methods cannot be constructed due to the nonseparability of the Hamiltonian. Thus, the available symplectic method is either fully implicit or semi-implicit, which decreases the efficiency because of the implicit iteration used during the evolution. In this paper, we aim to explore efficient explicit methods whose performance is mostly like symplectic methods for PN Hamiltonian systems. Taking the small parameter ε\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varepsilon $$\\end{document} appearing in PN terms into consideration, we replace the implicit symplectic solver with explicit solvers in the mixed symplectic method to solve the PN term and then derive three explicit methods. It is theoretically shown that the proposed methods are respectively second-order, fourth-order, and pseudo-fourth-order, and that their closeness to the corresponding symplectic methods are O(ε3h3),\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathscr {O}}(\\varepsilon ^{3}h^{3}),$$\\end{document}O(ε5h5),\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathscr {O}}(\\varepsilon ^{5}h^{5}),$$\\end{document} and O(ε3h3).\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathscr {O}}(\\varepsilon ^{3}h^{3}).$$\\end{document} That is, they are explicit near-symplectic methods with the presence of the small parameter ε.\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varepsilon .$$\\end{document} Numerical experiments with the Hamiltonian problem of spinning compact binaries show that the energy errors and orbital errors of the proposed explicit near-symplectic methods are indistinguishable from the corresponding mixed semi-implicit symplectic methods. The very small magnitude of the difference between the proposed explicit near-symplectic methods and the mixed symplectic methods confirms our theoretical analysis of their closeness to symplecticity. Finally, the much less CPU time consumed by the proposed methods highlights their most important advantage of high efficiency over the mixed symplectic methods.

Analysis of the geometry of the zero-velocity curves in the N-body ring problem depending on the mass ratio parameter

The purpose of the present study is the investigation of the effect of the mass ratio parameter, β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta$$\\end{document}, on the geometry of the zero-velocity curves, when there are N=3,4,…,100\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=3,4,\\ldots ,100$$\\end{document} peripheral bodies. It is well known that there is a bifurcation value of the β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta$$\\end{document} parameter in the N-body ring problem that produces a change in the number of stationary solutions in the system from 5N to 3N. By examining the behavior of the critical values of the Jacobi constant that define each of the zones of stationary solutions, we have unveiled the existence of other bifurcations or critical values of β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta$$\\end{document} in the scenario with 5N stationary solutions, which cause different changes in the geometry of the zero-velocity curves, which in turn affect the threshold for the total opening of the curves of zero velocity.

No infinite spin for planar total collision

The infinite spin problem is an old problem concerning the rotational behavior of total collision orbits in the n-body problem. It has long been known that when a solution tends to total collision then its normalized configuration curve must converge to the set of normalized central configurations. In the planar n-body problem every normalized configuration determines a circle of rotationally equivalent normalized configurations and, in particular, there are circles of normalized central configurations. It’s conceivable that by means of an infinite spin, a total collision solution could converge to such a circle instead of to a particular point on it. Here we prove that this is not possible, at least if the limiting circle of central configurations is isolated from other circles of central configurations. (It is believed that all central configurations are isolated, but this is not known in general.) Our proof relies on combining the center manifold theorem with the Łojasiewicz gradient inequality.

Evidence of second-order transition and critical scaling for the dynamical ordering transition in current-driven vortices

Dynamical ordering from a disordered plastic flow to an anisotropically ordered smectic flow induced by a dc force has been studied in various many-particle systems, including vortices in type-II superconductors. However, it remains unclear whether the dynamical ordering is a true phase transition because of lack of suitable experimental methods. Here, we study the response of vortex flow to the transverse force using a cross-shaped amorphous Mox\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{x}$$\\end{document}Ge1-x\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{1-x}$$\\end{document} film. From transverse current-voltage (force-velocity) characteristics under various longitudinal currents, we find a change of the transverse response in low voltage (velocity) regions from a nonlinear to linear behavior at a well-defined longitudinal current that marks the dynamical ordering transition. We also find the scaling collapse of the transverse current-voltage curves to a universal scaling function, providing evidence of the second-order transition for the dynamical ordering transition.

MatsubaraFunctions.jl: An equilibrium Green's function library in the Julia programming language

The Matsubara Green’s function formalism stands as a powerful technique for computing the thermodynamic characteristics of interacting quantum many-particle systems at finite temperatures. In this manuscript, our focus centers on introducing MatsubaraFunctions.jl, a Julia library that implements data structures for generalized n-point Green’s functions on Matsubara frequency grids. The package’s architecture prioritizes user-friendliness without compromising the development of efficient solvers for quantum field theories in equilibrium. Following a comprehensive introduction of the fundamental types, we delve into a thorough examination of key facets of the interface. This encompasses avenues for accessing Green’s functions, techniques for extrapolation and interpolation, as well as the incorporation of symmetries and a variety of parallelization strategies. Examples of increasing complexity serve to demonstrate the practical utility of the library, supplemented by discussions on strategies for sidestepping impediments to optimal performance.

Codebase release 0.1 for MatsubaraFunctions.jl

The Matsubara Green’s function formalism stands as a powerful technique for computing the thermodynamic characteristics of interacting quantum many-particle systems at finite temperatures. In this manuscript, our focus centers on introducing MatsubaraFunctions.jl, a Julia library that implements data structures for generalized n-point Green’s functions on Matsubara frequency grids. The package’s architecture prioritizes user-friendliness without compromising the development of efficient solvers for quantum field theories in equilibrium. Following a comprehensive introduction of the fundamental types, we delve into a thorough examination of key facets of the interface. This encompasses avenues for accessing Green’s functions, techniques for extrapolation and interpolation, as well as the incorporation of symmetries and a variety of parallelization strategies. Examples of increasing complexity serve to demonstrate the practical utility of the library, supplemented by discussions on strategies for sidestepping impediments to optimal performance.