Lennard-Jones Ideal Glass

• Lennard-Jones ideal glass is an amorphous solid with a flat potential energy landscape and local order resembling FCC, HCP, BCC, or icosahedral packings.
• The system exhibits distinctive thermodynamic markers and cooperative atomic rearrangements that pinpoint the glass transition via techniques like Monte Carlo and molecular dynamics.
• Advanced simulation methods and integral-equation theories precisely capture its energy landscape, indicating minimized defect networks and robust quasi-crystalline domains.

A Lennard-Jones ideal glass is a theoretically well-defined amorphous solid phase generated from atoms interacting via the canonical 12–6 Lennard-Jones pair potential. Distinguished by a unique flat potential energy landscape, cooperative atomic rearrangements, and a distinct local order dominated by quasi-crystalline motifs, the Lennard-Jones ideal glass serves as a central reference system for computational and theoretical investigations into glass formation, transition thermodynamics, and microscopic topological constraints. It exhibits neither long-range translational order nor true crystallinity, but is characterized by pervasive short-range order with local domains closely resembling face-centered cubic (FCC), hexagonal close-packed (HCP), body-centered cubic (BCC), or icosahedral atomic packing.

1. Simulation Frameworks and Potential Specifications

The formation of a Lennard-Jones ideal glass employs either single-component or binary mixtures of atoms subject to the pairwise 12–6 Lennard-Jones potential: $$V_{ab}(r) = 4\varepsilon_{ab}\left[ \left(\frac{\sigma_{ab}}{r}\right)^{12} - \left(\frac{\sigma_{ab}}{r}\right)^{6} \right],$$ with $\varepsilon_{ab}$ and $\sigma_{ab}$ specifying the energy and length scales, respectively (Abraham, 2015, Sun et al., 2017, Derlet, 2020, Jean-Marc et al., 27 Dec 2025). Typical simulation protocols include:

• Isothermal–isobaric (NPT) Monte Carlo: Sequential quenching from high temperature ($T^* = 1.0$) to low temperature ($T^* = 0.1$) at fixed $P^* = 1.0$, monitoring acceptance ratio, density, and enthalpy for equilibrium (Abraham, 2015).
• Classical molecular dynamics: Velocity-Verlet integration with Nosé–Hoover thermostat, long time annealing up to $14 \,\mu$s, and steepest-descent or conjugate-gradient minimizations to obtain inherent structures (Sun et al., 2017, Derlet, 2020).
• Integral-equation theory: Replicated hybrid mean spherical approximation (R-HMSA) with the Optimized Division Scheme (ODS) for the potential split, solving Ornstein–Zernike equations for real and reciprocal space correlation functions (Jean-Marc et al., 27 Dec 2025).

Periodic boundary conditions and cutoffs ($r_c \approx 2.5 \sigma$) are employed to accurately represent multiple coordination shells.

2. Identification of the Glass Transition and Thermodynamic Markers

The transition from supercooled liquid to glass is pinpointed by complementary thermodynamic and structural diagnostics:

• Density $\rho^*$, enthalpy $H^*$: Subtle kinks in quasi-linear temperature dependence near $T^* \approx 0.4$ indicate the onset of glassy behavior (Abraham, 2015).
• Pair distribution function $g(r)$: Bimodal splitting and emergence of additional peaks between first and second neighbor shells upon supercooling, with the extra peak at $r_{peak} \approx 1.45\,\sigma$ in the glass, interpreted as fcc-like local order (Jean-Marc et al., 27 Dec 2025).
• Structural metric $R^*=g_{min}/g_{max}$: Develops a distinct kink at the transition temperature.
• Dynamical transitions: In replicated liquid theory, the crossover temperature $T_{cr}$ is precisely identified where glass and liquid branches of excess free energy per particle intersect; values depend on density and typically, $T_{cr} \in [0.13, 0.30]$ for $\rho^* \in [0.90,1.00]$ (Jean-Marc et al., 27 Dec 2025).
• Configurational entropy $S_{conf}$: Emerges above the Kauzmann temperature $T_s \simeq 0.45\,\varepsilon/k_B$; $S_{conf}\to0$ marks the onset of the ideal glass regime (Sun et al., 2017).

3. Microscopic Structural Motifs and Bond-Topology Analysis

Detailed real-space analysis indicates that the Lennard-Jones ideal glass is not randomly disordered, but exhibits pervasive local environments closely matching crystalline reference geometries:

• Ackland–Jones method: Classifies each atom as quasi-FCC, HCP, BCC, icosahedral, or "other" via angular distribution functions and bond-angle cosines within the first neighbor shell. The root-mean-square deviation $\Delta_M(i)$ from ideal structural motifs is used for assignment (Abraham, 2015).
• Temperature-dependent populations: At high $T^*$ ($>0.6$), most atoms are "other" (>50%). Upon cooling below $T^*\approx0.4$, over 70% are classified quasi-crystalline (see Table), with icosahedral motifs remaining minor (1–3%) (Abraham, 2015).

$T^*$	Other	FCC	HCP	BCC	ICO
1.0	86.5	0.3	7.4	5.4	0.4
0.4	30.9	7.1	33.7	25.5	2.7
0.1	21.6	15.6	38.9	21.6	2.3

• SU(2) bond-topology (Nelson): In binary LJ glasses (Wahnström mixture), modified radical Voronoi tessellation coupled with SU(2) connectivity rules shows $>95\%$ of atoms satisfy the SU(2) closure law, indicating minimal disclination network frustration. Defect bonds are nearly eliminated (defect density $\leq0.3\%$) and mean bond order approaches the Frank–Kasper crystal limit (Derlet, 2020).

4. Potential Energy Landscape and Dynamics

A defining characteristic of the Lennard-Jones ideal glass is a flat and homogeneous potential energy surface (PES):

• Inherent structures: Minimization yields a landscape of local minima with a narrow distribution of transition barriers, $P(\Delta E)$ peaking at $\Delta E \approx 0.1\,\varepsilon$ (Sun et al., 2017).
• Cooperative diffusion: Mean-square displacement is almost arrested, but collective jumps between PES minima occur in brief bursts. The glass differs from the liquid, wherein atom motion is decoupled and barriers are high and variable (Sun et al., 2017).
• Thermally activated structural excitations: Localized events ("string" or "loop"-like), involving atoms violating SU(2) topology, tend to restore allowed geometric packings, lowering local frustration (Derlet, 2020).
• Frustration metrics: Excess bond energy due to deviation from equilibrium bond length ($f$) is minimized in the relaxed ideal glass.

5. Statistical and Integral-Equation Descriptions

Integral equation theory provides quantitative predictions for structural and thermodynamic properties in the ideal-glass phase:

• Replicated HMSA with ODS: Solves Ornstein–Zernike equations for both one-component and two-replica mixtures. The overlap order parameter $Q$ is large and persistent in the glass. The ideal-glass solution features a well-resolved additional peak in $g(r)$ between the main and second peaks, signaling local fcc-like structure despite absence of long-range periodicity or Bragg peaks (Jean-Marc et al., 27 Dec 2025).
• Transition quantification: Density-dependent threshold temperatures $T_0(\rho^*)$ mark the development of the structural peak (e.g., $T_0=0.10$ for $\rho^*=0.90$). The glass transition in excess free energy $T_{cr}$ closely tracks these changes.
• Free energy: In the ideal glass, configurational entropy vanishes below $T_s$ and the free-energy curve parallels that of the solid. Above $T_s$, $S_{conf}$ increases, lowering the free energy and leading to crossover with the liquid phase (Sun et al., 2017).

6. Physical Interpretation and Implications

The Lennard-Jones ideal glass exemplifies an amorphous solid phase where thermodynamic stability and mechanical rigidity are attributed to a percolated mosaic of nanoscopic quasi-crystalline domains. Key physical implications are:

• The glass transition correlates with the emergence of connectivity among quasi-crystalline clusters, manifested as sharp increases in viscosity and arrest of diffusive motion (Abraham, 2015, Sun et al., 2017).
• Compared to quasicrystals (nonperiodic long-range order) and single crystals (global translational symmetry), the Lennard-Jones ideal glass features exclusively short-range crystalline-like order, with local motifs nearly indistinguishable from ideal FCC/HCP/BCC geometries within disorder (Abraham, 2015, Jean-Marc et al., 27 Dec 2025).
• Minimization of bond-length frustration and defect network density in binary glasses maximizes the degree of topological order, driving the system toward lower cohesive energy and suppressed rearrangement activity (Derlet, 2020).
• The success of replica-based integral-equation approaches at ultra-low temperatures supports rigorous thermodynamic treatment and facilitates systematic study of ideal glass phases in atomistic systems without extrinsic disorder (Jean-Marc et al., 27 Dec 2025).
• Cooperative atomic diffusion and flat PES are the intrinsic features enabling configurational entropy generation and kinetic arrest, fundamental to glass transition theory (Sun et al., 2017).

The Lennard-Jones ideal glass thus represents a paradigm for understanding the interplay of local order, topology, energetics, and phase transitions in amorphous materials, bridging mean-field statistical theory and real-space structural analysis across simulation and theoretical frameworks.