When immiscible becomes miscible-Methane in water at high pressures

At low pressures, the solubility of gases in liquids is governed by Henry’s law, which states that the saturated solubility of a gas in a liquid is proportional to the partial pressure of the gas. As the pressure increases, most gases depart from this ideal behavior in a sublinear fashion, leveling off at pressures in the 1- to 5-kbar (0.1 to 0.5 GPa) range with solubilities of less than 1 mole percent (mol %). This contrasts strikingly with the well-known marked increase in solubility of simple gases in water at high temperature associated with the critical point (647 K and 212 bar). The solubility of the smallest hydrocarbon, the simple gas methane, in water under a range of pressure and temperature is of widespread importance, because it is a paradigmatic hydrophobe and occurs widely in terrestrial and extraterrestrial geology. We report measurements up to 3.5 GPa of the pressure dependence of the solubility of methane in water at 100°C—well below the latter’s critical temperature. Our results reveal a marked increase in solubility between 1 and 2 GPa, leading to a state above 2 GPa where the maximum solubility of methane in water exceeds 35 mol %

Efficacious calculation of Raman spectra in high pressure hydrogen

We present and evaluate an efficient method for simulating Raman spectra from molecular dynamics (MD) calculations {\it without} defining normal modes. We apply the method to high pressure hydrogen in the high-temperature "Phase IV": a plastic crystal in which the conventional picture of fixed phonon eigenmodes breaks down. Projecting trajectories onto in-phase molecular stretches is shown to be many orders of magnitude faster than polarisability calculations, allowing statistical averaging at high-temperature. The simulations are extended into metastable regimes and identify several regimes associated with symmetry-breaking on different timescales, which are shown to exhibit features in the Raman spectra at the current experimental limit of resolvability. In this paper we have concentrated on the methodology, a fuller description of the structure of Phase IV hydrogen is given in a previous paperComment: EHPRG conference 2013, High Pressure Research: Volume 34, Issue 2, 201

Twinning hierarchy, shape memory, and superelasticity demonstrated by molecular dynamics

A martensitic phase transition exhibiting shape memory, transformation-induced plasticity, or superelasticity typically involves a transformation between a high temperature, high symmetry phase and a low temperature, low symmetry phase. There have been numerous attempts using molecular dynamics to simulate the shape memory behavior, where the memory is stored in a twinned martensite and deformation occurs by motion of twin boundaries. However, the 3D case has always proved elusive, because suitable interatomic potentials to produce a unique low temperature phase are difficult to obtain. Here we present a study in which the binary Morse potential is tuned specifically to maximize the difference between L1_{0} and B19 (Strukturbericht notation, spacegroups P4/mmm and Pmma) structures. The twinned structure of martensite has been induced by gradually cooling the sample below the transition temperature. A bar-shaped sample was plastically deformed in the martensite phase, and on reheating above the transition temperature the initial shape was recovered. The effect of the shear-induced phase transition on the nanostructure of resulting martensite has also been investigated. An unusual discovery is that of a hierarchy of twins: nanotwins accommodate the mismatch between austenite and martensite at the habit plane, while dynamically created macrotwins are responsible for the deformation behavior and shape memory

Krypton and the fundamental flaw of the Lennard-Jones Potential

We have performed a series of neutron scattering experiments on supercritical krypton. Our data and analysis allow us to characterize the Frenkel line crossover in this model monatomic fluid. The data from our measurements was analyzed using Empirical Potential Structure Refinement to determine the short- and medium-range structure of the fluids. We find evidence for several shells of neighbors which form approximately concentric rings of density about each atom. The ratio of second to first shell radius is significantly larger than in any crystal structure. Modeling krypton using a Lennard-Jones potential is shown to give significant errors, notably that the liquid is overstructured. The true potential appears to be longer ranged and with a softer core than the 6–12 powerlaws permit

Crystal Structures of Dense Lithium: A Metal-Semiconductor-Metal Transition

Ab initio random structure searching and single-crystal x-ray diffraction have been used to determine the full structures of three phases of lithium, recently discovered at low temperature above 60 GPa. A structure with C2mb symmetry, calculated to be a poor metal, is proposed for the oC88 phase (60-65 GPa). The oC40 phase (65-95 GPa) is found to have a lowest-enthalpy structure with C2cb symmetry, in excellent agreement with the x-ray data. It is calculated to be a semiconductor with a band gap of similar to 1 eV at 90 GPa. oC24, stable above 95 GPa, has the space group Cmca, and refined atomic coordinates are in excellent agreement with previous calculations

Understanding high pressure hydrogen with a hierarchical machine-learned potential

The hydrogen phase diagram has a number of unusual features which are generally well reproduced by density functional calculations. Unfortunately, these calculations fail to provide good physical insights into why those features occur. In this paper, we parameterize a model potential for molecular hydrogen which permits long and large simulations. The model shows excellent reproduction of the phase diagram, including the broken-symmetry Phase II, an efficiently-packed phase III and the maximum in the melt curve. It also gives an excellent reproduction of the vibrational frequencies, including the maximum in the vibrational frequency $\nu(P)$ and negative thermal expansion. By detailed study of lengthy molecular dynamics, we give intuitive explanations for observed and calculated properties. All solid structures approximate to hexagonal close packed, with symmetry broken by molecular orientation. At high pressure, Phase I shows significant short-ranged correlations between molecular orientations. The turnover in Raman frequency is due to increased coupling between neighboring molecules, rather than weakening of the bond. The liquid is denser than the close-packed solid because, at molecular separations below 2.3\AA, the favoured relative orientation switches from quadrupole-energy-minimising to steric-repulsion-minimising. The latter allows molecules to get closer together, without atoms getting closer but this cannot be achieved within the constraints of a close-packed layer

Inelastic response of silicon to shock compression

The elastic and inelastic response of [001] oriented silicon to laser compression has been a topic of considerable discussion for well over a decade, yet there has been little progress in understanding the basic behaviour of this apparently simple material. We present experimental x-ray diffraction data showing complex elastic strain profiles in laser compressed samples on nanosecond timescales. We also present molecular dynamics and elasticity code modelling which suggests that a pressure induced phase transition is the cause of the previously reported 'anomalous' elastic waves. Moreover, this interpretation allows for measurement of the kinetic timescales for transition. This model is also discussed in the wider context of reported deformation of silicon to rapid compression in the literature