Single-simulation determination of phase boundaries: A dynamic Clausius–Clapeyron integration method

Atomistic study of the solid state inside graphene nanobubbles

A two-dimensional (2D) material placed on an atomically flat substrate can lead to the formation of surface nanobubbles trapping different types of substances. In this paper graphene nanobubbles of the radius of 7-34 nm with argon atoms inside are studied using molecular dynamics (MD). All modeled graphene nanobubbles except for the smallest ones exhibit an universal shape, i.e., a constant ratio of a bubble height to its footprint radius, which is in an agreement with experimental studies and their interpretation using the elastic theory of membranes. MD simulations reveal that argon does exist in a solid close-packed phase, although the internal pressure in the nanobubble is not sufficiently high for the ordinary crystallization that would occur in a bulk system. The smallest graphene bubbles with a radius of 7 nm exhibit an unusual "pancake" shape. Previously, nanobubbles with a similar pancake shape were experimentally observed in completely different systems at the interface between water and a hydrophobic surface.

Addressing the discrepancy of finding the equilibrium melting point of silicon using molecular dynamics simulations

We performed molecular dynamics simulations to study the equilibrium melting point of silicon using (i) the solid-liquid coexistence method and (ii) the Gibbs free energy technique, and compared our novel results with the previously published results obtained from the Monte Carlo (MC) void-nucleated melting method based on the Tersoff-ARK interatomic potential (Agrawal et al. Phys. Rev. B72, 125206. (doi:10.1103/PhysRevB.72.125206)). Considerable discrepancy was observed (approx. 20%) between the former two methods and the MC void-nucleated melting result, leading us to question the applicability of the empirical MC void-nucleated melting method to study a wide range of atomic and molecular systems. A wider impact of the study is that it highlights the bottleneck of the Tersoff-ARK potential in correctly estimating the melting point of silicon.

The Uhlenbeck-Ford model: Exact virial coefficients and application as a reference system in fluid-phase free-energy calculations

The Uhlenbeck-Ford (UF) model was originally proposed for the theoretical study of imperfect gases, given that all its virial coefficients can be evaluated exactly, in principle. Here, in addition to computing the previously unknown coefficients B₁₁ through B₁₃, we assess its applicability as a reference system in fluid-phase free-energy calculations using molecular simulation techniques. Our results demonstrate that, although the UF model itself is too soft, appropriately scaled Uhlenbeck-Ford (sUF) models provide robust reference systems that allow accurate fluid-phase free-energy calculations without the need for an intermediate reference model. Indeed, in addition to the accuracy with which their free energies are known and their convenient scaling properties, the fluid is the only thermodynamically stable phase for a wide range of sUF models. This set of favorable properties may potentially put the sUF fluid-phase reference systems on par with the standard role that harmonic and Einstein solids play as reference systems for solid-phase free-energy calculations.

Vibrational contribution to the thermodynamics of nanosized precipitates: vacancy-copper clusters in bcc-Fe

The effects of lattice vibration on the thermodynamics of nanosized coherent clusters in bcc-Fe consisting of vacancies and/or copper are investigated within the harmonic approximation. A combination of on-lattice simulated annealing based on Metropolis Monte Carlo simulations and off-lattice relaxation by molecular dynamics is applied to obtain the most stable cluster configurations at T = 0 K. The most recent interatomic potential built within the framework of the embedded-atom method for the Fe-Cu system is used. The total free energy of pure bcc-Fe and fcc-Cu as well as the total formation free energy and the total binding free energy of the vacancy-copper clusters are determined for finite temperatures. Our results are compared with the available data from previous investigations performed using many-body interatomic potentials and first-principles methods. For further applications in rate theory and object kinetic Monte Carlo simulations, the vibrational effects evaluated in the present study are included in the previously developed analytical fitting formulae.

Theoretical investigation of methane under pressure

We present computer simulations of liquid and solid phases of condensed methane at pressures below 25 GPa, between 150 and 300 K, where no appreciable molecular dissociation occurs. We used molecular dynamics (MD) and metadynamics techniques and empirical potentials in the rigid molecule approximation, whose validity was confirmed a posteriori by carrying out ab initio MD simulations for selected pressure and temperature conditions. Our results for the melting line are in satisfactory agreement with existing measurements. We find that the fcc crystal transforms into a hcp structure with four molecules per unit cell (B phase) at about 10 GPa and 150 K, and that the B phase transforms into a monoclinic high pressure phase above 20 GPa. Our results for solid/solid phase transitions are consistent with those of Raman studies but the phase boundaries estimated in our calculations are at higher pressure than those inferred from spectroscopic data.

Molecular dynamics simulations of surface-initiated melting of nitromethane

The melting of nitromethane initiated at solid-vacuum interfaces has been investigated using molecular dynamics nvt simulations with a realistic force field [D. C. Sorescu et al., J. Phys. Chem. B 104, 8406 (2000)]. The calculated melting point (251+/-5 K) is in good agreement with experiment (244.73 K) and values obtained previously (approximately 255.5 and 266.5+/-8 K) using other simulation methods [P. M. Agrawal et al., J. Chem. Phys. 119, 9617 (2003)]. Analyses of the molecular orientations and diffusion during the simulations as functions of the distance from the exposed surfaces show that the melting is a direct crystal-to-liquid transition, in which the molecules first gain rotational freedom, then mobility. There is a slight dependence of the melting temperature on the exposed crystallographic face.

Molecular dynamics studies of melting and solid-state transitions of ammonium nitrate

Molecular dynamics simulations are used to calculate the melting point and some aspects of high-temperature solid-state phase transitions of ammonium nitrate (AN). The force field used in the simulations is that developed by Sorescu and Thompson [J. Phys. Chem. A 105, 720 (2001)] to describe the solid-state properties of the low-temperature phase-V AN. Simulations at various temperatures were performed with this force field for a 4 x 4 x 5 supercell of phase-II AN. The melting point of AN was determined from calculations on this supercell with voids introduced in the solid structure to eliminate superheating effects. The melting temperature was determined by calculating the density and the nitrogen-nitrogen radial distribution functions as functions of temperature. The melting point was predicted to be in the range 445 +/- 10 K, in excellent agreement with the experimental value of 442 K. The computed temperature dependences of the density, diffusion, and viscosity coefficient for the liquid are in good agreement with experiment. Structural changes in the perfect crystal at various temperatures were also investigated. The ammonium ions in the phase-II structure are rotationally disordered at 400 K. At higher temperatures, beginning at 530 K, the nitrate ions are essentially rotationally unhindered. The density and radial distribution functions in this temperature range show that the AN solid is superheated. The rotational disorder is qualitatively similar to that observed in the experimental phase-II to phase-I solid-state transition.

Molecular dynamics studies of melting and some liquid-state properties of 1-ethyl-3-methylimidazolium hexafluorophosphate [emim][PF6]

Molecular dynamics simulations are used to study the liquid-state properties and melting of 1-ethyl-3-methylimidazolium hexafluorosphosphate [emim][PF6] using the force field of Canongia Lopes et al. [J. Phys. Chem. B 108, 2038 (2004)] and geometric constants from crystallographic data. The structures of the solid and liquid states are characterized by carbon-carbon, carbon-phosphorous, and phosphorous-phosphorous radial distribution functions. Spatial correlations among the ions are strong in the liquid state. The cohesive energy density and the temperature dependences of the molar volume and density of the liquid have been computed. The melting point is determined by equilibrating the solid-state supercells in which void defects have been introduced to eliminate the free-energy barrier for the formation of a solid-liquid interface. The computed melting point is 375+/-10 K, which is approximately 10% higher than the experimental value of 333 K.

Phase diagram of silicon from atomistic simulations

In this Letter we present a calculation of the temperature-pressure phase diagram of Si in a range of pressures covering from -5 to 20 GPa and temperatures up to the melting point. The phase boundaries and triple points between the diamond, liquid, beta-Sn, and Si34 clathrate phases are reported. We have employed efficient simulation techniques to calculate free energies and to numerically integrate the Clausius-Clapeyron equation, combined with a tight-binding model capable of an accuracy comparable to that of first-principles methods. The resulting phase diagram agrees well with the available experimental data.

Toward a robust and general molecular simulation method for computing solid-liquid coexistence

A rigorous and generally applicable method for computing solid-liquid coexistence is presented. The method overcomes some of the technical difficulties associated with other solid-liquid simulation procedures and can be implemented within either a molecular dynamics or Monte Carlo framework. The method consists of three steps: First, relative Gibbs free energy curves are created for the solid and liquid phases using histogram reweighting. Next, the free energy difference between the solid and liquid phases is evaluated at a single state point by integrating along a pseudosupercritical transformation path that connects the two phases. Using this result, the solid and liquid free energy curves are referenced to a common point, allowing a single coexistence point to be determined. Finally, Gibbs-Duhem integration is used to determine the full coexistence curve. To evaluate its utility, this method is applied to the Lennard-Jones and NaCl systems. Results for solid-liquid coexistence agree with previous calculations for these systems. In addition, it is shown that the NaCl model does not correctly describe solid-liquid coexistence at high pressures. An analysis of the accuracy of the method indicates that the results are most sensitive to the transformation free energy calculation.

Transitions between disordered phases in supercooled liquid silicon

We have investigated the transitions between disordered phases in supercooled liquid silicon using computer simulations. The thermodynamic properties were directly obtained from the free energy, which was computed using the recently proposed reversible scaling method. The calculated free energies of the crystalline and liquid phases of silicon at zero pressure, obtained using the environment dependent interatomic potential, are in excellent agreement with the available experimental data. The results show that, at zero pressure, a weak first-order liquid-liquid transition occurs at 1135 K and a continuous liquid-amorphous transition takes place at 843 K. These results are consistent with the existence of a second critical point for the liquid-liquid transition at a negative pressure.

Synergistic science

Should computational materials science be recognized as a field with a role in the community comparable to computational physics or chemistry? With the emergence of multiscale modelling, the answer is a resounding 'yes'.