Nuclear Quantum Effects on the Electronic Structure of Water and Ice.

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

An experimental and computational view of the photoionization of diol-water clusters.

In the interstellar medium, diols and other prebiotic molecules adsorb onto icy mantles surrounding dust grains. Water in the ice may affect the reactivity and photoionization of these diols. Ethylene glycol (EG), 1,2-propylene glycol, and 1,3-propylene glycol clusters with water clusters were used as a proxy to study these interactions. The diol-water clusters were generated in a continuous supersonic molecular beam, photoionized by synchrotron-based vacuum ultraviolet light from the Advanced Light Source, and subsequently detected by reflectron time-of-flight mass spectrometry. The appearance energies for the detected clusters were determined from the mass spectra, collected at increasing photon energy. Clusters of both diol fragments and unfragmented diols with water were detected. The lowest energy geometry optimized conformers for the observed EG-water clusters and EG fragment-water clusters have been visualized using density functional theory (DFT), providing insight into hydrogen bonding networks and how these affect fragmentation and appearance energy. As the number of water molecules clustered around EG fragments (m/z 31 and 32) increased, the appearance energy for the cluster decreased, indicating a stabilization by water. This trend was supported by DFT calculations. Fragment clusters from 1,2-propylene glycol exhibited a similar trend, but with a smaller energy decrease, and no trend was observed from 1,3-propylene glycol. We discuss and suggest that the reactivity and photoionization of diols in the presence of water depend on the size of the diol, the location of the hydroxyl group, and the number of waters clustered around the diol.

New Trick for an Old Dog: From Prediction to Properties of "Hidden Clathrates" Ba2Zn5As6 and Ba2Zn5Sb6.

The zinc-antimony phase space has been heavily investigated due to the structural complexity and abundance of high-performing thermoelectric materials. Consequentially, the desire to use zinc and antimony as framework elements to encage rattling cations and achieve phonon-glass-electron-crystal-type properties has remained an enticing goal with only two alkali metal clathrates to date, Cs8Zn18Sb28 and K58Zn122Sb207. Guided by Zintl electron-counting predictions, we explored the Ba-Zn-Pn (Pn = As, Sb) phase space proximal to the expected composition of the type-I clathrate. In situ powder X-ray diffraction studies revealed two "hidden" compounds which can only be synthesized in a narrow temperature range. The ex situ synthesis and crystal growth unveiled that instead of type-I clathrates, compositionally close but structurally different new clathrate-like compounds formed, Ba2Zn5As6 and Ba2Zn5Sb6. These materials crystallize in a unique structure, in the orthorhombic space group Pmna with the Wyckoff sequence i2h6gfe. Single-phase synthesis enabled the exploration of their transport properties. Rattling of the Ba cations in oversized cages manifested low thermal conductivity, which, coupled with the high Seebeck coefficients observed, are prerequisites for a promising thermoelectric material. Potential for further optimization of the thermoelectric performance by aliovalent doping was computationally analyzed.

UV-Visible Absorption Spectra of Solvated Molecules by Quantum Chemical Machine Learning.

Predicting UV-visible absorption spectra is essential to understand photochemical processes and design energy materials. Quantum chemical methods can deliver accurate calculations of UV-visible absorption spectra, but they are computationally expensive, especially for large systems or when one computes line shapes from thermal averages. Here, we present an approach to predict UV-visible absorption spectra of solvated aromatic molecules by quantum chemistry (QC) and machine learning (ML). We show that a ML model, trained on the high-level QC calculation of the excitation energy of a set of aromatic molecules, can accurately predict the line shape of the lowest-energy UV-visible absorption band of several related molecules with less than 0.1 eV deviation with respect to reference experimental spectra. Applying linear decomposition analysis on the excitation energies, we unveil that our ML models probe vertical excitations of these aromatic molecules primarily by learning the atomic environment of their phenyl rings, which align with the physical origin of the π →π* electronic transition. Our study provides an effective workflow that combines ML with quantum chemical methods to accelerate the calculations of UV-visible absorption spectra for various molecular systems.

Effect of sodium chloride adsorption on the surface premelting of ice.

We characterise the structural properties of the quasi-liquid layer (QLL) at two low-index ice surfaces in the presence of sodium chloride (Na+/Cl-) ions by molecular dynamics simulations. We find that the presence of a high surface density of Na+/Cl- pairs changes the surface melting behaviour from step-wise to gradual melting. The ions lead to an overall increase of the thickness and the disorder of the QLL, and to a low-temperature roughening transition of the air-ice interface. The local molecular structure of the QLL is similar to that of liquid water, and the differences between the basal and primary prismatic surface are attenuated by the presence of Na+/Cl- pairs. These changes modify the crystal growth rates of different facets and the solvation environment at the surface of sea-water ice with a potential impact on light scattering and environmental chemical reactions.

Unprecedented superstructure in the type I family of clathrates.

The first arsenic-based clathrate exhibiting superstructural ordering due to optimization of Au-As, As-As, and Ba-Au bonding is reported. Ba8Au16As30 crystallizes in a unique P21/c monoclinic clathrate structure. The synthesis, crystal and electronic structure, and transport properties are discussed.

Engineering Thermal Transport across Layered Graphene-MoS2 Superlattices.

Layering two-dimensional van der Waals materials provides a high degree of control over atomic placement, which could enable tailoring of vibrational spectra and heat flow at the sub-nanometer scale. Here, using spatially resolved ultrafast thermoreflectance and spectroscopy, we uncover the design rules governing cross-plane heat transport in superlattices assembled from monolayers of graphene (G) and MoS2 (M). Using a combinatorial experimental approach, we probe nine different stacking sequences, G, GG, MG, GGG, GMG, GGMG, GMGG, GMMG, and GMGMG, and identify the effects of vibrational mismatch, interlayer adhesion, and junction asymmetry on thermal transport. Pure G sequences display evidence of quasi-ballistic transport, whereas adding even a single M layer strongly disrupts heat conduction. The experimental data are described well by molecular dynamics simulations, which include thermal expansion, accounting for the effect of finite temperature on the interlayer spacing. The simulations show that an increase of ∼2.4% in the layer separation of GMGMG, relative to its value at 300 K, can lead to a doubling of the thermal resistance. Using these design rules, we experimentally demonstrate a five-layer GMGMG superlattice "thermal metamaterial" with an ultralow effective cross-plane thermal conductivity comparable to that of air.

Solid Solution Yb2-xCaxCdSb2: Structure, Thermoelectric Properties, and Quality Factor.

Solid solutions of Yb2-xAxCdSb2 (A = Ca, Sr, Eu; x ≤ 1) are of interest for their promising thermoelectric (TE) properties. Of these solid solutions, Yb2-xCaxCdSb2 has end members with different crystal structures. Yb2CdSb2 crystallizes in the polar space group Cmc21, whereas Ca2CdSb2 crystallizes in the centrosymmetric space group Pnma. Other solid solutions, Yb2-xAxCdSb2 (A = Sr, Eu), crystallize in the polar space group for x ≤ 1, and compositions with x ≥ 1 have not been reported. Both structure types are composed of corner-sharing CdSb4 tetrahedra condensed into sheets that differ by the stacking of the layers. Single crystals of the solid solution Yb2-xCaxCdSb2 (x = 0-1) were studied to elucidate the structural transition between the Yb2CdSb2 and Ca2CdSb2 structure types. For x ≤ 1, the structures remain in the polar space group Cmc21. As the Ca content is increased, a positional disorder arises in the intralayer cation sites (Yb2/Ca2) and the Cd site, resulting in inversion of the CdSb4 tetrahedral chain. This phenomenon could be indicative of an intergrowth of the opposing space group. The TE properties of polycrystalline samples of Yb2-xCaxCdSb2 (x ≤ 1) were measured from 300 to 525 K. The lattice thermal conductivity is extremely low (0.3-0.4 W/m·K) and the Seebeck coefficients are high (100-180 µV/K) across the temperature range. First-principles calculations show a minimum in the thermal conductivity for the x = 0.3 composition, in good agreement with experimental data. The low thermal conductivity stems from the acoustic branches being confined to low frequencies and a large number of phonon scattering channels provided by the localized optical branches. The TE quality factor of the Yb1.7A0.3CdSb2 (A = Ca, Sr, Eu) series has been calculated and predicts that the A = Ca and Sr solid solutions may not improve with carrier concentration optimization but that the Eu series is worthy of additional modifications. Overall, the x = 0.3 compositions provide the highest zT because they provide the best electronic properties with the lowest thermal conductivity.

Direct observation of ultrafast hydrogen bond strengthening in liquid water.

Water is one of the most important, yet least understood, liquids in nature. Many anomalous properties of liquid water originate from its well-connected hydrogen bond network1, including unusually efficient vibrational energy redistribution and relaxation2. An accurate description of the ultrafast vibrational motion of water molecules is essential for understanding the nature of hydrogen bonds and many solution-phase chemical reactions. Most existing knowledge of vibrational relaxation in water is built upon ultrafast spectroscopy experiments2-7. However, these experiments cannot directly resolve the motion of the atomic positions and require difficult translation of spectral dynamics into hydrogen bond dynamics. Here, we measure the ultrafast structural response to the excitation of the OH stretching vibration in liquid water with femtosecond temporal and atomic spatial resolution using liquid ultrafast electron scattering. We observed a transient hydrogen bond contraction of roughly 0.04 Å on a timescale of 80 femtoseconds, followed by a thermalization on a timescale of approximately 1 picosecond. Molecular dynamics simulations reveal the need to treat the distribution of the shared proton in the hydrogen bond quantum mechanically to capture the structural dynamics on femtosecond timescales. Our experiment and simulations unveil the intermolecular character of the water vibration preceding the relaxation of the OH stretch.

Ultrahigh Convergent Thermal Conductivity of Carbon Nanotubes from Comprehensive Atomistic Modeling.

Anomalous heat transport in one-dimensional nanostructures, such as nanotubes and nanowires, is a widely debated problem in condensed matter and statistical physics, with contradicting pieces of evidence from experiments and simulations. Using a comprehensive modeling approach, comprised of lattice dynamics and molecular dynamics simulations, we proved that the infinite length limit of the thermal conductivity of a (10,0) single-wall carbon nanotube is finite but this limit is reached only for macroscopic lengths due to a thermal phonon mean free path of several millimeters. Our calculations showed that the extremely high thermal conductivity of this system at room temperature is dictated by quantum effects. Modal analysis showed that the divergent nature of thermal conductivity, observed in one-dimensional model systems, is suppressed in carbon nanotubes by anharmonic scattering channels provided by the flexural and optical modes with polarization in the plane orthogonal to the transport direction.

Mn-intercalated MoSe2 under pressure: Electronic structure and vibrational characterization of a dilute magnetic semiconductor.

Intercalation offers a promising way to alter the physical properties of two-dimensional (2D) layered materials. Here, we investigate the electronic and vibrational properties of 2D layered MoSe2 intercalated with atomic manganese at ambient and high pressure up to 7 GPa by Raman scattering and electronic structure calculations. The behavior of optical phonons is studied experimentally with a diamond anvil cell and computationally through density functional theory calculations. Experiment and theory show excellent agreement in optical phonon behavior. The previously Raman inactive A2u mode is activated and enhanced with intercalation and pressure, and a new Raman mode appears upon decompression, indicating a possible onset of a localized structural transition, involving the bonding or trapping of the intercalant in 2D layered materials. Density functional theory calculations reveal a shift of the Fermi level into the conduction band and spin polarization in MnxMoSe2 that increases at low Mn concentrations and low pressure. Our results suggest that intercalation and pressurization of van der Waals materials may allow one to obtain dilute magnetic semiconductors with controllable properties, providing a viable route for the development of new materials for spintronic applications.

Bathochromic Shift in the UV-Visible Absorption Spectra of Phenols at Ice Surfaces: Insights from First-Principles Calculations.

Some organic pollutants in snowpacks undergo faster photodegradation than in solution. One possible explanation for such effect is that their UV-visible absorption spectra are shifted toward lower energy when the molecules are adsorbed at the air-ice interface. However, such bathochromic shift is difficult to measure experimentally. Here, we employ a multiscale/multimodel approach that combines classical and first-principles molecular dynamics, quantum chemical methods, and statistical learning to compute the light absorption spectra of two phenolic molecules in different solvation environments at the relevant thermodynamic conditions. Our calculations provide an accurate estimate of the bathochromic shift of the lowest-energy UV-visible absorption band when these molecules are adsorbed at the air-ice interface, and they shed light into its molecular origin.

Photodecay of guaiacol is faster in ice, and even more rapid on ice, than in aqueous solution.

Snowpacks contain a wide variety of inorganic and organic compounds, including some that absorb sunlight and undergo direct photoreactions. How the rates of these reactions in, and on, ice compare to rates in water is unclear: some studies report similar rates, while others find faster rates in/on ice. Further complicating our understanding, there is conflicting evidence whether chemicals react more quickly at the air-ice interface compared to in liquid-like regions (LLRs) within the ice. To address these questions, we measured the photodegradation rate of guaiacol (2-methoxyphenol) in various sample types, including in solution, in ice, and at the air-ice interface of nature-identical snow. Compared to aqueous solution, we find modest rate constant enhancements (increases of 3- to 6-fold) in ice LLRs, and much larger enhancements (of 17- to 77-fold) at the air-ice interface of nature-identical snow. Our computational modeling suggests the absorption spectrum for guaiacol red-shifts and increases on ice surfaces, leading to more light absorption, but these changes explain only a small portion (roughly 2 to 9%) of the observed rate constant enhancements in/on ice. This indicates that increases in the quantum yield are primarily responsible for the increased photoreactivity of guaiacol on ice; relative to solution, our results suggest that the quantum yield is larger by a factor of roughly 3-6 in liquid-like regions and 12-40 at the air-ice interface.

Carbon dioxide, bicarbonate and carbonate ions in aqueous solutions under deep Earth conditions.

We investigate the effect of pressure, temperature and acidity on the composition of water-rich carbon-bearing fluids under thermodynamic conditions that correspond to the Earth's deep crust and upper mantle. Our first-principles molecular dynamics simulations provide mechanistic insight into the hydration shell of carbon dioxide, bicarbonate and carbonate ions, and into the pathways of the acid/base reactions that convert these carbon species into one another in aqueous solutions. At temperatures of 1000 K and higher, our simulations can sample the chemical equilibrium of these acid/base reactions, thus allowing us to estimate the chemical composition of diluted carbon dioxide and (bi)carbonate ions as a function of acidity and thermodynamic conditions. We find that, especially at the highest temperature, the acidity of the solution is essential to determine the stability domain of CO2vs. HCO3-vs. CO32-.

Thermal transport in amorphous small organic materials: a mechanistic study.

Understanding the thermal transport mechanisms in amorphous organic materials is of great importance to solve hot-spot issues in organic-electronics nanodevices. Here we studied thermal transport in two popular molecular electronic materials, N,N-dicarbazolyl-3,5-benzene (mCP) and N,N'-diphenyl-N,N'-di(3-methylphenyl)-(1,1'-biphenyl)-4,4'diamine (TPD), in the amorphous state by molecular dynamics simulations. We found that due to the softness of organic materials, the low thermal conductivity of both systems can be greatly enhanced under pressure. Notably, in such systems, the convective term of heat flux provides an important contribution to thermal transport as it cross-correlates with the Virial term in the Green-Kubo formula. Mode diffusivity calculations reveal that low-frequency modes can contribute significantly to thermal transport in both mCP and TPD. By increasing the pressure, the sound velocity and relaxation time of such low-frequency modes can be enhanced, and a part of these modes converts from diffusons to propagons. The cooperation of these three effects is responsible for the strong pressure dependence of thermal transport in amorphous organic systems. Molecular pair heat flux calculations demonstrate that heat transfer mainly happens between pairs of molecules with distances below 1.4 nm. This work paves the way for the optimization of thermal transport in amorphous organic materials widely used in opto-electronics, e.g. as OLED and OPV.

III-V Clathrate Semiconductors with Outstanding Hole Mobility: Cs8In27Sb19 and A8Ga27Sb19 (A = Cs, Rb).

Three novel unconventional clathrates with unprecedented III-V semiconducting frameworks have been synthesized: Cs8In27Sb19, Cs8Ga27Sb19, and Rb8Ga27Sb19. These clathrates represent the first examples of tetrel-free clathrates that are completely composed of main group elements. All title compounds crystallize in an ordered superstructure of clathrate-I in the Ia3̅ space group (No. 206; Z = 8). In the clathrate framework, a full ordering of {Ga or In} and Sb is observed by a combination of high-resolution synchrotron single-crystal and powder X-ray diffraction techniques. Density functional theory (DFT) calculations show that all three clathrates are energetically stable with relaxed lattice constants matching the experimental data. Due to the complexity of the crystal structure composed of heavy elements, the reported clathrates exhibit ultralow thermal conductivities of less than 1 W·m-1·K-1 at room temperature. All compounds are predicted and experimentally confirmed to be narrow-bandgap p-type semiconductors with high Seebeck thermopower values, up to 250 µV·K-1 at 300 K for Cs8In27Sb19. The latter compound shows carrier concentrations and mobilities, 1.42 × 1015 cm-3 and 880 cm2 ·V-1·s-1, which are on par with the values for parent binary InSb, one of the best electronic semiconductors. The high hole carrier mobility is uncommon for complex bulk materials and a highly desirable trait, opening ways to design semiconducting materials based on tunable III-V clathrates.

Influence of thermostatting on nonequilibrium molecular dynamics simulations of heat conduction in solids.

Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here, we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green's function method in the ballistic regime and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it into the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification.

Publisher Correction: An electrochemical thermal transistor.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Modeling heat transport in crystals and glasses from a unified lattice-dynamical approach.

We introduce a novel approach to model heat transport in solids, based on the Green-Kubo theory of linear response. It naturally bridges the Boltzmann kinetic approach in crystals and the Allen-Feldman model in glasses, leveraging interatomic force constants and normal-mode linewidths computed at mechanical equilibrium. At variance with molecular dynamics, our approach naturally and easily accounts for quantum mechanical effects in energy transport. Our methodology is carefully validated against results for crystalline and amorphous silicon from equilibrium molecular dynamics and, in the former case, from the Boltzmann transport equation.