G. Franzese,G. Malescio,A. Skibinsky,S. V. Buldyrev,H. E. Stanley
Abstract:We investigate the phase behavior of a single-component system in 3 dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature {\bf 409}, 692 (2001)] that, even with no evidences of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas--low-density liquid (LDL) critical point, and the other in a gas--high-density liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the 3-parameter space of the soft-core potential and we perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to explore the behavior of liquid - liquid phase transitions and its relationship with the crystal structure under specific soft - core potential conditions in a single - component system. Specifically:
1. **Research Background**: The paper focuses on a single - component system with spherically symmetric, pairwise additive, soft - core interactions, which include long - range attractions, medium - range repulsive shoulders, and short - range hard - core repulsions. This type of potential is often used to describe liquid systems such as colloids, protein solutions, or liquid metals.
2. **Known Phenomena**: Previous studies have shown that even without evidence of density anomalies, there are two first - order fluid - fluid phase transitions in the phase diagram of this system, ending at the gas - low - density liquid (LDL) critical point and the gas - high - density liquid (HDL) critical point respectively, and there is also an LDL - HDL phase transition at low temperatures.
3. **Research Objectives**:
- **Detailed Exploration of the Phase Diagram**: Explore the three - parameter space of the soft - core potential through integral equation calculations and molecular dynamics simulations, especially the fluid metastable phases in the supercooled phase.
- **Analysis of the Crystal Structure**: Analyze the structure of the crystal phase and find that it does not depend on density within the considered density range.
- **Study of Density Anomalies**: Prove the non - existence of density anomalies in this system through explicit thermodynamic calculations and propose that this phenomenon may be related to the existence of only one stable crystal structure in the system.
4. **Main Contributions**:
- **Detailed Analysis of the Phase Diagram**: Through integral equations and molecular dynamics simulations, the phase diagram of the system is analyzed in detail, especially the fluid metastable phases in the supercooled phase.
- **Confirmation of the Crystal Structure**: Through molecular dynamics simulations, it is confirmed that there is only one stable crystal structure in the system.
- **Explanation of Density Anomalies**: Through thermodynamic calculations, the existence of density anomalies is excluded, and the relationship between density anomalies and crystal structure diversity is proposed.
In summary, this paper aims to deeply explore the liquid - liquid phase transition behavior of a single - component system under specific soft - core potential conditions through theoretical calculations and simulations, especially the relationship between these phase transitions and the crystal structure of the system.