Heat and charge transport in H2O at ice-giant conditions from ab initio molecular dynamics simulations

Combined experiments and molecular simulations for understanding the thermo-responsive behavior and gelation of methylated glucans with different glycosidic linkages

HYPOTHESIS

Elucidation of the micro-mechanisms of sol-gel transition of gelling glucans with different glycosidic linkages is crucial for understanding their structure-property relationship and for various applications. Glucans with distinct molecular chain structures exhibit unique gelation behaviors. The disparate gelation phenomena observed in two methylated glucans, methylated (1,3)-β-d-glucan of curdlan (MECD) and methylated (1,4)-β-d-glucan of cellulose (MC), notwithstanding their equivalent degrees of substitution, are intricately linked to their unique molecular architectures and interactions between glucan and water.

EXPERIMENTS

Density functional theory and molecular dynamics simulations focused on the electronic property distinctions between MECD and MC, alongside conformational variations during thermal gelation. Inline attenuated total reflection Fourier transform infrared spectroscopy tracked secondary structure alterations in MECD and MC. To corroborate the simulation results, additional analyses including circular dichroism, rheology, and micro-differential scanning calorimetry were performed.

FINDINGS

Despite having similar thermally induced gel networks, MECD and MC display distinct physical gelation patterns and molecular-level conformational changes during gelation. The network of MC gel was formed via a "coil-to-ring" transition, followed by ring stacking. In contrast, the MECD gel comprised compact irregular helices accompanied by notable volume shrinkage. These variations in gelation behavior are ascribed to heightened hydrophobic interactions and diminished hydrogen bonding in both systems upon heating, resulting in gelation. These findings provide valuable insights into the microstructural changes during gelation and the thermo-gelation mechanisms of structurally similar polysaccharides.

Seebeck Coefficient of Ionic Conductors from Bayesian Regression Analysis

We propose a novel approach to evaluating the ionic Seebeck coefficient in electrolytes from relatively short equilibrium molecular dynamics simulations, based on the Green-Kubo theory of linear response and Bayesian regression analysis. By exploiting the probability distribution of the off-diagonal elements of a Wishart matrix, we develop a consistent and unbiased estimator for the Seebeck coefficient, whose statistical uncertainty can be arbitrarily reduced in the long-time limit. We assess the efficacy of our method by benchmarking it against extensive equilibrium molecular dynamics simulations conducted on molten CsF using empirical force fields. We then employ this procedure to calculate the Seebeck coefficient of molten NaCl, KCl, and LiCl using neural network force fields trained on ab initio data over a range of pressure-temperature conditions.

Thermal Conductivity of Water at Extreme Conditions

Measuring the thermal conductivity (κ) of water at extreme conditions is a challenging task, and few experimental data are available. We predict κ for temperatures and pressures relevant to the conditions of the Earth mantle, between 1,000 and 2,000 K and up to 22 GPa. We employ close to equilibrium molecular dynamics simulations and a deep neural network potential fitted to density functional theory data. We then interpret our results by computing the equation of state of water on a fine grid of points and using a simple model for κ. We find that the thermal conductivity is weakly dependent on temperature and monotonically increases with pressure with an approximate square-root behavior. In addition, we show how the increase of κ at high pressure, relative to ambient conditions, is related to the corresponding increase in the sound velocity. Although the relationships between the thermal conductivity, pressure and sound velocity established here are not rigorous, they are sufficiently accurate to allow for a robust estimate of the thermal conductivity of water in a broad range of temperatures and pressures, where experiments are still difficult to perform.

Magnesium oxide-water compounds at megabar pressure and implications on planetary interiors

Magnesium Oxide (MgO) and water (H₂O) are abundant in the interior of planets. Their properties, and in particular their interaction, significantly affect the planet interior structure and thermal evolution. Here, using crystal structure predictions and ab initio molecular dynamics simulations, we find that MgO and H₂O can react again at ultrahigh pressure, although Mg(OH)₂ decomposes at low pressure. The reemergent MgO-H₂O compounds are: Mg₂O₃H₂ above 400 GPa, MgO₃H₄ above 600 GPa, and MgO₄H₆ in the pressure range of 270-600 GPa. Importantly, MgO₄H₆ contains 57.3 wt % of water, which is a much higher water content than any reported hydrous mineral. Our results suggest that a substantial amount of water can be stored in MgO rock in the deep interiors of Earth to Neptune mass planets. Based on molecular dynamics simulations we show that these three compounds exhibit superionic behavior at the pressure-temperature conditions as in the interiors of Uranus and Neptune. Moreover, the water-rich compound MgO₄H₆ could be stable inside the early Earth and therefore may serve as a possible early Earth water reservoir. Our findings, in the poorly explored megabar pressure regime, provide constraints for interior and evolution models of wet planets in our solar system and beyond.

Ultralow Melting Temperature of High-Pressure Face-Centered Cubic Superionic Ice

Superionic ice with oxygen in a face-centered cubic (fcc) sublattice is ascribed to the origin of magnetic fields of Uranus and Neptune, since the melting temperature (T_m) of fcc-superionic ice is believed to be higher than the isentropes of ice giants. However, precisely measuring the fcc-superionic phase experimentally remains a difficult task. The majority of the systematic investigations of its T_m were performed using perfect oxygen fcc-sublattice computations, which could result in superheating and overestimation of T_m. On the basis of the ab initio molecular dynamics method and the model with H₂O vacancy, we avoid superheating and obtain a much lower T_m than previous reports, indicating that fcc-superionic ice cannot exist in the interiors of Uranus and Neptune. Further simulations with the two-phase method justify the conclusion. The results suggest that superheating should be seriously treated when simulating the phase diagram of other hydrogen-related superionic states, which are widely used to understand the properties of ice giants, Earth, and Venus.

Investigating finite-size effects in molecular dynamics simulations of ion diffusion, heat transport, and thermal motion in superionic materials

The effects of the finite size of the simulation box in equilibrium molecular dynamics simulations are investigated for prototypical superionic conductors of different types, namely, the fluorite-structure materials PbF₂, CaF₂, and UO₂ (type II), and the α phase of AgI (type I). Largely validated empirical force-fields are employed to run ns-long simulations and extract general trends for several properties, at increasing size and in a wide temperature range. This work shows that, for the considered type-II superionic conductors, the diffusivity dramatically depends on the system size and that the superionic regime is shifted to larger temperatures in smaller cells. Furthermore, only simulations of several hundred atoms are able to capture the experimentally observed, characteristic change in the activation energy of the diffusion process, occurring at the order-disorder transition to the superionic regime. Finite-size effects on ion diffusion are instead much weaker in α-AgI. The thermal conductivity is found generally smaller for smaller cells, where the temperature-independent (Allen-Feldman) regime is also reached at significantly lower temperatures. The finite-size effects on the thermal motion of the non-mobile ions composing the solid matrix follow the simple law that holds for solids.

Metallic Aluminum Suboxides with Ultrahigh Electrical Conductivity at High Pressure

Aluminum, as the most abundant metallic elemental content in the Earth's crust, usually exists in the form of alumina (Al₂O₃). However, the oxidation state of aluminum and the crystal structures of aluminum oxides in the pressure range of planetary interiors are not well established. Here, we predicted two aluminum suboxides (Al₂O, AlO) and two superoxides (Al₄O₇, AlO₃) with uncommon stoichiometries at high pressures using first-principle calculations and crystal structure prediction methods. We find that the P4/nmm Al₂O becomes stable above ~765 GPa and may survive in the deep mantles or cores of giant planets such as Neptune. Interestingly, the Al₂O and AlO are metallic and have electride features, in which some electrons are localized in the interstitials between atoms. We find that Al₂O has an electrical conductivity one order of magnitude higher than that of iron under the same pressure-temperature conditions, which may influence the total conductivity of giant planets. Our findings enrich the high-pressure phase diagram of aluminum oxides and improve our understanding of the interior structure of giant planets.

Invariance principles in the theory and computation of transport coefficients

ABSTRACT

In this work, we elaborate on two recently discovered invariance principles, according to which transport coefficients are, to a large extent, independent of the microscopic definition of the densities and currents of the conserved quantities being transported (energy, momentum, mass, charge). The first such principle, gauge invariance, allows one to define a quantum adiabatic energy current from density-functional theory, from which the heat conductivity can be uniquely defined and computed using equilibrium ab initio molecular dynamics. When combined with a novel topological definition of atomic oxidation states, gauge invariance also sheds new light onto the mechanisms of charge transport in ionic conductors. The second principle, convective invariance, allows one to extend the analysis to multi-component systems. These invariance principles can be combined with new spectral analysis methods for the current time series to be fed into the Green-Kubo formula to obtain accurate estimates of transport coefficients from relatively short molecular dynamics simulations.

First principles study of dense and metallic nitric sulfur hydrides

Studies of molecular mixtures containing hydrogen sulfide (H₂S) could open up new routes towards hydrogen-rich high-temperature superconductors under pressure. H₂S and ammonia (NH₃) form hydrogen-bonded molecular mixtures at ambient conditions, but their phase behavior and propensity towards mixing under pressure is not well understood. Here, we show stable phases in the H₂S-NH₃ system under extreme pressure conditions to 4 Mbar from first-principles crystal structure prediction methods. We identify four stable compositions, two of which, (H₂S) (NH₃) and (H₂S) (NH₃)₄, are stable in a sequence of structures to the Mbar regime. A re-entrant stabilization of (H₂S) (NH₃)₄ above 300 GPa is driven by a marked reversal of sulfur-hydrogen chemistry. Several stable phases exhibit metallic character. Electron-phonon coupling calculations predict superconducting temperatures up to 50 K, in the Cmma phase of (H₂S) (NH₃) at 150 GPa. The present findings shed light on how sulfur hydride bonding and superconductivity are affected in molecular mixtures. They also suggest a reservoir for hydrogen sulfide in the upper mantle regions of icy planets in a potentially metallic mixture, which could have implications for their magnetic field formation.