Franz Saija

Phase diagram of the Gaussian-core model.

We trace with high numerical accuracy the phase diagram of the Gaussian-core model, a classical system of point particles interacting via a Gaussian-shaped, purely repulsive potential. This model, which provides a reliable qualitative description of the thermal behavior of interpenetrable globular polymers, is known to exhibit a polymorphic fcc-bcc transition at low densities and reentrant melting at high densities. Extensive Monte Carlo simulations, carried out in conjunction with accurate calculations of the solid free energies, lead to a thermodynamic scenario that is partially modified with respect to previous knowledge. In particular, we find that: (i) the fluid-bcc-fcc triple-point temperature is about one third of the maximum freezing temperature; (ii) upon isothermal compression, the model exhibits a fluid-bcc-fcc-bcc-fluid sequence of phases in a narrow range of temperatures just above the triple point. We discuss these results in relation to the behavior of star-polymer solutions and of other softly repulsive systems.

Ab initio molecular dynamics study of dissociation of water under an electric field.

The behavior of liquid water under an electric field is a crucial phenomenon in science and engineering. However, its detailed description at a microscopic level is difficult to achieve experimentally. Here we report on the first ab initio molecular-dynamics study on water under an electric field. We observe that the hydrogen-bond length and the molecular orientation are significantly modified at low-to-moderate field intensities. Fields beyond a threshold of about 0.35  V/Å are able to dissociate molecules and sustain an ionic current via a series of correlated proton jumps. Upon applying even more intense fields (∼1.0  V/Å), a 15%-20% fraction of molecules are instantaneously dissociated and the resulting ionic flow yields a conductance of about 7.8  Ω-1 cm-1, in good agreement with experimental values. This result paves the way to quantum-accurate microscopic studies of the effect of electric fields on aqueous solutions and, thus, to massive applications of ab initio molecular dynamics in neurobiology, electrochemistry, and hydrogen economy.

Phase diagram of softly repulsive systems: The Gaussian and inverse-power-law potentials

We redraw, using state-of-the-art methods for free-energy calculations, the phase diagrams of two reference models for the liquid state: the Gaussian and inverse-power-law repulsive potentials. Notwithstanding the different behaviors of the two potentials for vanishing interparticle distances, their thermodynamic properties are similar in a range of densities and temperatures, being ruled by the competition between the body-centered-cubic (bcc) and face-centered-cubic (fcc) crystalline structures and the fluid phase. We confirm the existence of a reentrant bcc phase in the phase diagram of the Gaussian-core model, just above the triple point. We also trace the bcc-fcc coexistence line of the inverse-power-law model as a function of the power exponent n and relate the common features in the phase diagrams of such systems to the softness degree of the interaction.

Miller experiments in atomistic computer simulations

Significance In 1953, Stanley Miller reported on the spontaneous formation of glycine when applying an electric discharge on a mixture of simple molecules, giving birth to modern research on the origins of life. The effect of electric fields on mixtures of simple molecules is presently studied in computer simulations at the quantum level, and Miller results are reproduced for the first time, to our knowledge, in atomistic simulations, as glycine forms spontaneously only in the presence of electric fields. However, this occurs through reaction pathways more complex than believed, identifying formamide as a key compound in prebiotic chemistry. Moreover, electric fields are naturally present at mineral surfaces, suggesting a potentially crucial role in the biogeochemistry of both the primordial and the modern Earth. The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.

Hexatic phase in the two-dimensional Gaussian-core model

We present a Monte Carlo simulation study of the phase behavior of two-dimensional classical particles repelling each other through an isotropic Gaussian potential. As in the analogous three-dimensional case, a reentrant-melting transition occurs upon compression for not too high temperatures, along with a spectrum of waterlike anomalies in the fluid phase. However, in two dimensions melting is a continuous two-stage transition, with an intermediate hexatic phase which becomes increasingly more definite as pressure grows. All available evidence supports the Kosterlitz-Thouless-Halperin-Nelson-Young scenario for this melting transition. We expect that such a phenomenology can be checked in confined monolayers of charge-stabilized colloids with a softened core.

Evaluation of phenomenological one-phase criteria for the melting and freezing of softly repulsive particles

We test the validity of some widely used phenomenological criteria for the localization of the fluid-solid transition thresholds against the phase diagrams of particles interacting through the exp-6, inverse-power-law, and Gaussian potentials. We find that one-phase rules give, on the whole, reliable estimates of freezing/melting points. The agreement is ordinarily better for a face-centered-cubic solid than for a body-centered-cubic crystal, even more so in the presence of a pressure-driven reentrant transition of the solid into a denser fluid phase, as found in the Gaussian-core model.

Anomalous phase behavior in a model fluid with only one type of local structure

We present evidence that the concurrent existence of two populations of particles with different effective diameters is not a prerequisite for the occurrence of anomalous phase behaviors in systems of particles interacting through spherically symmetric unbounded potentials. Our results show that an extremely weak softening of the interparticle repulsion, which yields a single nearest-neighbor separation, is able to originate a wide spectrum of unconventional features including reentrant melting, solid polymorphism, as well as thermodynamic, dynamic, and structural anomalies. These findings extend the possibility of anomalous phase behavior to a class of systems much broader than currently assumed.

Anomalous phase behavior of a soft-repulsive potential with a strictly monotonic force

We study the phase behavior of a classical system of particles interacting through a strictly convex soft-repulsive potential which, at variance with the pairwise softened repulsions considered so far in the literature, lacks a region of downward or zero curvature. Nonetheless, such interaction is characterized by two length scales, owing to the presence of a range of interparticle distances where the repulsive force increases, for decreasing distance, much more slowly than in the adjacent regions. We investigate, using extensive Monte Carlo simulations combined with accurate free-energy calculations, the phase diagram of the system under consideration. We find that the model exhibits a fluid-solid coexistence line with multiple re-entrant regions, an extremely rich solid polymorphism with solid-solid transitions, and waterlike anomalies. In spite of the isotropic nature of the interparticle potential, we find that among the crystal structures in which the system can exist, there are also a number of non-Bravais lattices, such as cI16 and diamond.

Anomalous melting behavior under extreme conditions: Hard matter turning “soft”

We show that a system of particles interacting through the exp-6 pair potential, commonly used to describe effective interatomic forces under high compression, exhibits anomalous melting features such as reentrant melting and a rich solid polymorphism, including a stable BC8 crystal. We relate this behavior to the crossover, with increasing pressure, between two different regimes of local order that are associated with the two repulsive length scales of the potential. Our results provide a unifying picture for the high-pressure melting anomalies observed in many elements and point out that, under extreme conditions, atomic systems may reveal surprising similarities with soft matter.

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

We study a two-dimensional fluid of particles interacting through a spherically symmetric and marginally soft two-body repulsion. This model can exist in three different crystal phases, one of them with square symmetry and the other two triangular. We show that, while the triangular solids first melt into a hexatic fluid, the square solid is directly transformed on heating into an isotropic fluid through a first-order transition, with no intermediate tetratic phase. In the low-pressure triangular and square crystals, melting is reentrant provided the temperature is not too low, but without the necessity of two competing nearest-neighbor distances over a range of pressures. A whole spectrum of water-like fluid anomalies completes the picture for this model potential.