Amorphous Radial Frustration and Water-Like Anomalies in a Ramp-Shoulder Fluid

Technical Summary: Amorphous Radial Frustration and Water-Like Anomalies in a Ramp-Shoulder Fluid

Problem Statement
The study addresses the microscopic origins of water-like anomalies (non-monotonic variations in density, diffusivity, and structural order under compression or cooling) in soft matter systems. While core-softened (CS) potentials have long been used to model these phenomena, the prevailing view attributes anomalies to the competition between two characteristic interaction distances (shell competition). However, recent evidence suggests this "two-scale" explanation is incomplete, as anomalies can arise from cooperative structural reorganization even without distinct competing scales. Specifically, the role of the detailed shape of the softened interaction region—particularly the presence of shallow attractive shoulders in polymer-grafted nanoparticle systems—and its influence on amorphization and frustration in three dimensions remains an open question.

Methodology
The authors investigate a three-dimensional coarse-grained ramp-shoulder fluid derived from effective interactions between polymer-grafted nanoparticles. The effective potential ($U_{CS}$) combines a Lennard-Jones repulsive core with a series of Gaussian terms representing polymer-mediated effects, resulting in a softened repulsive ramp with a shallow attractive shoulder.

• Simulation: Molecular dynamics simulations were performed using the LAMMPS package with 1024 particles in a cubic box with periodic boundary conditions. Equilibration was conducted in the NVE ensemble, followed by production runs in the NPT ensemble using Nosé–Hoover thermostats and barostats.
• Analysis: The study characterizes thermodynamic, structural, and dynamic properties through:
  • Thermodynamic response functions (isothermal compressibility $\kappa_T$, thermal expansion coefficient $\alpha_P$, specific heat $c_P$).
  • Structural metrics: Radial distribution function $g(r)$, pair excess entropy ($S_2$), cumulative two-body entropy ($C_{s2}$), translational order parameter ($\tau$), and orientational order parameters ($q_6$, $Q_6$).
  • Dynamic metrics: Diffusion coefficient ($D$) calculated via the Einstein relation.
  • Phase identification using PTM analysis to distinguish BCC crystalline, amorphous, and fluid regions.

Key Contributions and Results
The simulations reveal a phase diagram containing crystalline (BCC), amorphous, and fluid regions, exhibiting density, diffusion, and structural anomalies. The primary findings challenge the conventional hierarchy of anomalies in isotropic CS fluids:

1. Partial Decoupling of Anomalies: Unlike standard CS fluids where the structural anomaly typically precedes the diffusion anomaly, which precedes the temperature of maximum density (TMD), this system shows a partial decoupling. The density anomaly (TMD line) extends to significantly higher pressures than the structural anomaly. The diffusion anomaly is closely linked to amorphization and shell migration rather than simple local shell competition.
2. Mechanism of Cooperative Radial Restructuring: The anomalies are not driven solely by the competition between two fixed interaction scales. Instead, they emerge from a regime of cooperative radial restructuring. In this regime, radial correlations increase (particles redistribute between inner and outer shells) without the development of long-range orientational crystalline order.
3. Role of the Attractive Shoulder: The shallow attractive shoulder in the interaction potential is critical. It stabilizes open-shell local arrangements over a broad pressure interval. Consequently, compression does not simply eliminate open shells to form a compact crystal; instead, it forces particles to continuously redistribute between competing local organizations. This leads to amorphous radial frustration.
4. Structural Pathways:
   • Low Pressure: Shell competition is thermally driven; increasing temperature shifts population to shorter scales.
   • Intermediate/High Pressure: Shell competition becomes compression-driven. The system exhibits a crossover where the outer shell population decreases while the inner shell increases.
   • Frustration: Despite increasing radial order (evidenced by rising $g(r)$ peaks and coordination numbers), the orientational order parameter ($Q_6$) remains low throughout the anomalous region. The system remains trapped in a frustrated amorphous state rather than crystallizing.
5. Diffusion Anomaly: The diffusion coefficient exhibits a minimum at intermediate pressures. Further compression destabilizes open-shell arrangements, promoting migration to shorter scales and partially restoring mobility ($\partial D / \partial P > 0$). This dynamic anomaly coincides with the pressure interval where the BCC phase loses stability and amorphous configurations dominate.

Significance and Claims
The paper claims that the detailed shape of the softened interaction region, specifically the inclusion of a shallow attractive shoulder, strongly influences the structural pathways explored under compression. The authors posit that the anomalous behavior in this system cannot be understood merely by the existence of two interaction scales. Instead, the anomalies arise from a specific structural pathway: a regime of amorphous radial frustration where radial organization grows while orientational crystalline order is suppressed.

This work suggests that amorphous radial frustration is a distinct mechanism for generating water-like anomalies in soft-core systems, offering a framework that connects anomalies, amorphization, and frustrated ordering. The results indicate that the competition between scales is a necessary but insufficient condition for anomalies; the system's ability to reorganize cooperatively without crystallizing is the determining factor.