Search

A wide variety of two-dimensional layered materials are synthesized by liquid-phase exfoliation. Here we examine exfoliation of MoS2 into nanosheets in a mixture of water and isopropanol (IPA) containing cavitation bubbles. Using force fields optimized with experimental data on interfacial energies between MoS2 and the solvent, multimillion-atom molecular dynamics simulations are performed in conjunction with experiments to examine shock-induced collapse of cavitation bubbles and the resulting exfoliation of MoS2. The collapse of cavitation bubbles generates high-speed nanojets and shock waves in the solvent. Large shear stresses due to the nanojet impact on MoS2 surfaces initiate exfoliation, and shock waves reflected from MoS2 surfaces enhance exfoliation. Structural correlations in the solvent indicate that shock induces an ice VII like motif in the first solvation shell of water.

The mechanism of how ice crystal form has been extensively studied by many researchers but remains an open question. Molecular dynamics (MD) simulations are a useful tool for investigating the molecular-scale mechanism of crystal formation. However, the timescale of phenomena that can be analyzed by MD simulations is typically restricted to microseconds or less, which is far too short to explore ice crystal formation that occurs in real systems. In this study, a metadynamics (MTD) method was adopted to overcome this timescale limitation of MD simulations. An MD simulation combined with the MTD method, in which two discrete oxygen–oxygen radial distribution functions represented by Gaussian window functions were used as collective variables, successfully reproduced the formation of several different ice crystals when the Gaussian window functions were set at appropriate oxygen–oxygen distances: cubic ice, stacking disordered ice consisting of cubic ice and hexagonal ice, high-pressure ice VII, layered ice with an ice VII structure, and layered ice with an unknown structure. The free-energy landscape generated by the MTD method suggests that the formation of each ice crystal occurred via high-density water with a similar structure to the formed ice crystal. The present method can be used not only to study the mechanism of crystal formation but also to search for new crystals in real systems.

It has been known for decades that certain aqueous salt solutions of LiCl and LiBr readily form glasses when cooled to below ≈160 K. This fact has recently been exploited to produce a « salty » high-pressure ice form: When the glass is compressed at low temperatures to pressures higher than 4 GPa and subsequently warmed, it crystallizes into ice VII with the ionic species trapped inside the ice lattice. Here we report the extreme limit of salt incorporation into ice VII, using high pressure neutron diffraction and molecular dynamics simulations. We show that high-pressure crystallisation of aqueous solutions of LiCl∙RH2O and LiBr∙RH2O with R = 5.6 leads to solids with strongly expanded volume, a destruction of the hydrogen-bond network with an isotropic distribution of water-dipole moments, as well as a crystal-to-amorphous transition on decompression. This highly unusual behaviour constitutes an interesting pathway from a glass to a crystal where translational periodicity is restored but the rotational degrees of freedom remaining completely random.

We investigated molecular dissociation induced by 10-keV X-ray irradiation in dense ice at pressures up to 40 GPa at 300 K. The dissociation yield estimated from the oxygen K-edge X-ray Raman spectra, showed that the molecular dissociation was enhanced up to 14 GPa and gradually suppressed on further compression to 40 GPa. The molecular dissociation was detected for a rather narrow pressure span of 2–40 GPa by the X-ray spectroscopy. The pressure variation of the dissociation yield was similar to that observed in the electric conductivity of ice VII and likely interpreted in terms of proton mobility.

The neutron diffraction pattern of D2O ice was recently measured at pressures up to 52 GPa by Guthrie et al., who proposed an octahedral interstitial model for ice at pressures above 13 GPa to account for the deviation of the observed crystal structure from that of ice VII. In this article, the octahedral interstitial model was re-examined in terms of the interstitial occupancy and X-ray Raman spectroscopy (XRS) spectra. The interstitial occupancy calculated using first-principles molecular dynamics simulations was negligibly small compared to that of the interstitial model. The oxygen K-edge spectra calculated for the interstitial model exhibited two additional low-energy peaks originating from water molecules and hydroxides that are interacting with interstitial protons, respectively, whereas these low-energy peaks were not observed in the experimentally measured spectra. These results suggest that the interstitial model cannot explain the XRS spectra of ice VII at pressures above 13 GPa and that more precise structure measurements and analyses are necessary to reveal the nature of the pressure-induced transition.

Complex intra-molecular interactions and the hydrogen-bonding network in H2O-volatile mixtures play critical roles in many dynamics processes in physical chemistry, biology and Earth and planetary sciences. We used high pressure Raman spectroscopy to study the pressure evolution of vibrational frequencies and bonding behavior in H2O-CH3OH mixtures. We found that the presence of low CH3OH content in H2O increases the transition pressure where water crystallizes to ice VI, but does not significantly change the pressure where ice VI transforms to ice VII. Furthermore, the stiffening rates of C-H stretching frequencies dω/dP in CH3OH significantly decrease upon the crystallization of water and the softening rates of the O-H stretching frequencies of ice VII are suppressed over a narrow pressure range, after which the frequencies of these modes shift with pressure in ways similar to pure CH3OH and ice VII, respectively. Such complex pressure evolution of Raman frequencies along with pronounced variations in Raman intensities of CH3OH within the sample and the hysteresis of the water-ice VI phase transition suggest pressure-induced segregation of low content CH3OH from ice VII. These findings indicate the significant influence of volatiles on the crystallization of sub-surface ocean and thermal evolution within large icy planets and satellites.

It was discovered that a peak appears near a pressure of Pc = 10 GPa in the electrical conductivity of ice VII as measured through impedance spectroscopy in a diamond anvil cell (DAC) during the process of compression from 2 GPa to 40 GPa at room temperature. The activation energy for the conductivity measured in the cooling/heating process between 278 K and 303 K reached a minimum near Pc. Theoretical modelling and molecular dynamics simulations suggest that the origin of this unique peak is the transition of the major charge carriers from the rotational defects to the ionic defects.

Raman spectra of ice VII and X were investigated using strong plasma shockwave generated by laser-induced breakdown (LIB) in liquid water. Simultaneously, the occurrence of the hydrogen emission lines of 656 nm (Hα), 486 nm (Hβ), 434 nm (Hγ) and 410 nm (Hδ) was observed. At 5 × 10(12) W/cm(2) optical power density, the O-H symmetric stretching, translational and librational modes of ice VII and a single peak at 785 cm(-1) appeared in the spectra. The band was assigned to the Raman-active O-O mode of the monomolecular phase, which was the symmetric hydrogen bond of cuprite ice X. The spectra indicated that ice VII and X structure were formed, as the trajectory of the strong plasma shockwave passes through the stable Pressure-Temperature range of ice VII and X. The shockwave temperature and pressure were calculated by the Grüneisen model.

The structure and phase transition of high-pressure ice are of long-standing interest and challenge, and there is still a huge gap between theoretical and experimental understanding. The quantum nature of protons such as delocalization, quantum tunneling and zero-point motion is crucial to the comprehension of the properties of high-pressure ice. Here we investigated the temperature-induced phase transition and oxygen K-edge x-ray absorption spectra of ice VII, VIII and X using ab initio path-integral molecular dynamics simulations. The tremendous difference between experiments and the previous theoretical predictions is closed for the phase diagram of ice below 300 K at pressures up to 110 GPa. Proton tunneling assists the proton-ordered ice VIII to transform into proton-disordered ice VII where only thermal activated proton-transfer cannot occur. The oxygen K edge with its shift is sensitive to the order-disorder transition, and therefore can be applied to diagnose the dynamics of ice structures.