Hexatic phase and water-like anomalies in a two-dimensional fluid of particles with a weakly softened core

Core-softened colloid under extreme geometrical confinement

Geometrical constraints offer a promising strategy for assembling colloidal crystal structures that are not typically observed in bulk or under 2D conditions. Core-softened colloids, in particular, have emerged as versatile chemical building blocks with applications across various scientific and technological areas. In this study, we investigate the behavior of a core-softened model confined between two parallel walls. Employing molecular dynamics simulations, we analyze the system's response under extreme confinement, where only one or two layers of colloids are permitted. The system comprises particles modeled by a ramp-like potential confined within slit nanoslits created by two flat, purely repulsive walls with a lateral side L separated by a distance Lz. Through a systematic analysis of the phase behavior as Lz increases, or as the system undergoes decompression, for different values of L, we identified a mono-to-bilayer transition associated with changes in the colloidal structure. In the monolayer regime, we observed solid phases at lower densities than those observed in the 2D case. Importantly, we demonstrated that confinement at specific Lz values, allowing particle arrangement into two layers, can lead to the emergence of the square phase, which was not observed under monolayer or 2D conditions. By correlating thermodynamic, translational, and orientational ordering, as well as the dynamics of this confined colloidal system, our findings offer valuable insights into the utilization of geometrical constraints to induce and manipulate structural changes.

Signatures of continuous hexatic-liquid transition in two-dimensional melting

Recent studies have shown that the melting of two-dimensional crystals can be either continuous or discontinuous, relying on multiple parameters such as particle stiffness, density, and particle size dispersity. However, what determines the continuity or discontinuity of the two-dimensional melting remains elusive. Here we study the two-dimensional melting of binary mixtures of soft-core particles. The two particle species are different in either particle size or particle stiffness. Starting with the mono-component systems which exhibit discontinuous hexatic-liquid transition, we gradually increase the particle size or stiffness dispersity and find that the hexatic-liquid coexistent region shrinks and eventually vanishes above a critical dispersity. Therefore, the growth of disorder caused by the particle size or stiffness dispersity leads to the discontinuous-continuous transition of the two-dimensional melting. We further find that as long as the melting is continuous the defect concentrations on the boundary between hexatic and liquid phases remain almost constant, accompanied by an almost constant correlation length. These characteristic defect concentrations and correlation length are universal and independent of particle interactions, temperature, and type of particle dispersity, which act as signatures of the continuous two-dimensional melting.

Phase behavior of a two-dimensional core-softened system: new physical insights

In this work, we report results of extensive computer simulations regarding the phase behavior of a core-softened system. By using structural and thermodynamic descriptors, as well as self-diffusion coefficients, we provide a comprehensive view of the rich phase behavior displayed by the particular instance of the model studied in here. Our calculations agree with previously published results focused on a smaller region in the temperature-density parameter space (Dudalov et al 2014 Soft Matter 10 4966). In this work, we explore a broader region in this parameter space, and uncover interesting fluid phases with low-symmetry local order, that were not reported by previous works. Solid phases were also found, and have been previously characterized in detail by (Kryuchkov et al 2018 Soft Matter 14 2152). Our results support previously reported findings, and provide new physical insights regarding the emergence of order as disordered phases transform into solids by providing radial distribution function maps and specific heat data. Our results are summarized in terms of a phase diagram.

Event-chain Monte Carlo simulations of the liquid to solid transition of two-dimensional decagonal colloidal quasicrystals

In event-chain Monte Carlo simulations, we model colloidal particles in two dimensions that interact according to an isotropic short-ranged pair potential which supports the two typical length scales present in decagonal quasicrystals. We investigate the assembled structures as we vary the density and temperature. Our special interest is related to the transition from quasicrystal to liquid. In contrast to the KTHNY melting theory for quasicrystals which predicts an intermediate pentahedratic phase, we find a one-step first-order melting transition. However, we discover that the slow relaxation of phasonic flips, i.e. rearrangements of the particles due to additional degrees of freedom in quasicrystals, changes the positional correlation functions, to the extent that structures with long-range orientational correlations, but exponentially decaying positional correlations, are observed.

Hysteresis in the temperature dependence of the IR bending vibration of deeply cooled confined water

Using Fourier Transform Infrared (FTIR) spectroscopy, we investigate the temperature dependence of the bending vibrations of water confined in the pores of a silica hydrogel in the temperature interval of 270-180 K. We also investigate the presence of thermal hysteresis by cooling and reheating temperature scans. The results clearly show the presence, at about 230 K, of a crossover in the temperature dependence of the IR spectra; moreover, the presence of hysteresis is clearly demonstrated. By comparing FTIR data with neutron diffraction data and previous calorimetric data on the same samples, we conclude that the crossover and the hysteretical behavior do not involve a water glass transition or crystallization but are related to a first-order-like liquid-liquid transition.

Hierarchy of anomalies in the two-dimensional Mercedes-Benz model of water

We investigate by Monte Carlo simulations density, diffusion, and structural anomalies of the simple two-dimensional Mercedes-Benz (MB) model of water, which is a very simple toy model for explaining the origin of water properties. MB water molecules are modeled as two-dimensional Lennard-Jones disks, with three orientation-dependent hydrogen-bonding arms, arranged as in the MB logo. The model is in a way also a variance of silica-like models. Beside the known thermodynamic anomaly for the model we also found diffusion and structural anomalies and map out the cascade of density, structural, pair entropy, and diffusivity anomalies for MB model. The orientational order parameters with three and six-fold symmetry were determined and maximum for each one observed. The anomalies occur in hierarchy order, which is a slight variation of the hierarchy order in real water. The diffusion anomaly region is the innermost in the hierarchy while for water it is the density anomaly region.

Waterlike anomalies in a two-dimensional core-softened potential

We investigate the structural, thermodynamic, and dynamic behavior of a two-dimensional (2D) core-corona system using Langevin dynamics simulations. The particles are modeled by employing a core-softened potential which exhibits waterlike anomalies in three dimensions. In previous studies in a quasi-2D system a new region in the pressure versus temperature phase diagram of structural anomalies was observed. Here we show that for the two-dimensional case two regions in the pressure versus temperature phase diagram with structural, density, and diffusion anomalies are observed. Our findings indicate that, while the anomalous region at lower densities is due the competition between the two length scales in the potential at higher densities, the anomalous region is related to the reentrance of the melting line.

Density Affects the Nature of the Hexatic-Liquid Transition in Two-Dimensional Melting of Soft-Core Systems

We find that both continuous and discontinuous hexatic-liquid transitions can happen in the melting of two-dimensional solids of soft-core disks. For three typical model systems, Hertzian, harmonic, and Gaussian-core models, we observe the same scenarios. These systems exhibit reentrant crystallization (melting) with a maximum melting temperature T_{m} happening at a crossover density ρ_{m}. The hexatic-liquid transition at a density smaller than ρ_{m} is discontinuous. Liquid and hexatic phases coexist in a density interval, which becomes narrower with increasing temperature and tends to vanish approximately at T_{m}. Above ρ_{m}, the transition is continuous, in agreement with the Kosterlitz-Thouless-Halperin-Nelson-Young theory. For these soft-core systems, the nature of the hexatic-liquid transition depends on density (pressure), with the melting at ρ_{m} being a plausible transition point from discontinuous to continuous hexatic-liquid transition.

Grain-boundary-induced melting in quenched polycrystalline monolayers

Melting in two dimensions can successfully be explained with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario which describes the formation of the high-symmetry phase with the thermal activation of topological defects within an (ideally) infinite monodomain. With all state variables being well defined, it should hold also as freezing scenario where oppositely charged topological defects annihilate. The Kibble-Zurek mechanism, on the other hand, shows that spontaneous symmetry breaking alongside a continuous phase transition cannot support an infinite monodomain but leads to polycrystallinity. For any nonzero cooling rate, critical fluctuations will be frozen out in the vicinity of the transition temperature. This leads to domains with different director of the broken symmetry, separated by a defect structure, e.g., grain boundaries in crystalline systems. After instantaneously quenching a colloidal monolayer from a polycrystalline to the isotropic fluid state, we show that such grain boundaries increase the probability for the formation of dislocations. In addition, we determine the temporal decay of defect core energies during the first few Brownian times after the quench. Despite the fact that the KTHNY scenario describes a continuous phase transition and phase equilibrium does not exist, melting in polycrystalline samples starts at grain boundaries similar to first-order phase transitions.

Random pinning changes the melting scenario of a two-dimensional core-softened potential system

In experiments two-dimensional systems are realized mainly on solid substrates, which introduce quenched disorder due to some inherent defects. The defects of substrates influence the melting scenario of the systems and have to be taken into account in the interpretation of experimental results. We present the results of molecular dynamics simulations of a two-dimensional system with a core-softened potential in which a small fraction of the particles is pinned, inducing quenched disorder. Ppotentials of this type are widely used for the qualitative description of systems with waterlike anomalies. In our previous publications it was shown that the system demonstrates an anomalous melting scenario: at low densities the system melts through two continuous transitions in accordance with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory with an intermediate hexatic phase, while at high densities the conventional first-order melting transition takes place. We find that the well-known disorder-induced widening of the hexatic phase occurs at low densities, while in the high-density part of the phase diagram random pinning transforms the first-order melting into two transitions: a continuous KTHNY-like solid-hexatic transition and a first-order hexatic-isotropic liquid transition.

Specific heat in two-dimensional melting

We report the specific heat cN around the melting transition(s) of micrometer-sized superparamagnetic particles confined in two dimensions, calculated from fluctuations of positions and internal energy, and corresponding Monte Carlo simulations. Since colloidal systems provide single particle resolution, they offer the unique possibility to compare the experimental temperatures of the peak position of cN(T) and symmetry breaking, respectively. While order parameter correlation functions confirm the Kosterlitz-Thouless-Halperin-Nelson-Young melting scenario where translational and orientational order symmetries are broken at different temperatures with an intermediate so called hexatic phase, we observe a single peak of the specific heat within the hexatic phase, with excellent agreement between experiment and simulation. Thus, the peak is not associated with broken symmetries but can be explained with the total defect density, which correlates with the maximum increase of isolated dislocations. The absence of a latent heat strongly supports the continuous character of both transitions.

How dimensionality changes the anomalous behavior and melting scenario of a core-softened potential system?

We present a computer simulation study of the phase diagram and anomalous behavior of two-dimensional (2D) and three-dimensional (3D) classical particles repelling each other through an isotropic core-softened potential. As in the analogous three-dimensional case, in 2D a reentrant-melting transition occurs upon compression under not too high pressure, along with a spectrum of thermodynamic and dynamic anomalies in the fluid phase. However, in two dimensions the order of the region of anomalous diffusion and the region of structural anomaly is inverted in comparison with the 3D case, where there exists a water-like sequence of anomalies, and has a silica-like sequence. In the low density part of the 2D phase diagram, melting is a continuous two-stage transition, with an intermediate hexatic phase. All available evidence supports the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario for this melting transition. On the other hand, at high density part of the phase diagram one first-order transition takes place.

Anomalous melting scenario of the two-dimensional core-softened system

We present a computer simulation study of the phase behavior of two-dimensional (2D) classical particles repelling each other through an isotropic core-softened potential. As in the analogous three-dimensional (3D) case, a reentrant-melting transition occurs upon compression for not too high pressures, along with a spectrum of waterlike anomalies in the fluid phase. However, in two dimensions in the low density part of the phase diagram melting is a continuous two-stage transition, with an intermediate hexatic phase. All available evidence supports the Kosterlitz-Thouless-Halperin-Nelson-Young scenario for this melting transition. On the other hand, at the high density part of the phase diagram one first-order transition takes place.

Volume crossover in deeply supercooled water adiabatically freezing under isobaric conditions

The irreversible return of a supercooled liquid to stable thermodynamic equilibrium often begins as a fast process which adiabatically drives the system to solid-liquid coexistence. Only at a later stage will solidification proceed with the expected exchange of thermal energy with the external bath. In this paper we discuss some aspects of the adiabatic freezing of metastable water at constant pressure. In particular, we investigated the thermal behavior of the isobaric gap between the molar volume of supercooled water and that of the warmer ice-water mixture which eventually forms at equilibrium. The available experimental data at ambient pressure, extrapolated into the metastable region within the scheme provided by the reference IAPWS-95 formulation, show that water ordinarily expands upon (partially) freezing under isenthalpic conditions. However, the same scheme also suggests that, for increasing undercoolings, the volume gap is gradually reduced and eventually vanishes at a temperature close to the currently estimated homogeneous ice nucleation temperature. This behavior is contrasted with that of substances which do not display a volumetric anomaly. The effect of increasing pressures on the alleged volume crossover from an expanded to a contracted ice-water mixture is also discussed.