Symmetric Hydrogen Bonds Research Articles

Effect of temperature on the symmetrization of AlOOH hydrogen bonds

Hydrogen bonding symmetrization of water-bearing minerals plays a key role in the water cycle of the deep Earth. The lower mantle is a high-temperature and high-pressure environment, and in this study a classical molecular dynamics simulation of the δ-AlOOH hydrogen bond symmetrization process is carried out using a deep learning potential model. Calculation the length of the HO bond reveals that AlOOH undergoes hydrogen bond symmetrization when the pressure reaches 9 GPa. The volume thermal expansion profile of AlOOH shows an anomalous phenomenon of negative correlation between phase transition pressure and temperature. The radial distribution functions gOH(r) and gHH(r) show the asymmetry of hydrogen bonding at 8 GPa and indicate the critical conditions for hydrogen bonding symmetry at 10 GPa. Calculations of the O1O2 bond lengths show that the critical pressure required for hydrogen bond symmetry increases with temperature. Our simulation results provide some insight into the hydrogen bonding symmetry of water-bearing minerals in the deep Earth’s.

Probing the Effects of Size and Charge on the Monohydration and Dihydration of SiF5- and SiF62- via Comparisons with BF4- and PF6.

This study systematically examines the interactions of the trigonal bipyramidal silicon pentafluoride and octahedral silicon hexafluoride anions with either one or two water molecules, (SiF5-(H2O)n and SiF62-(H2O)n, respectively, where n = 1, 2). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using the CCSD(T) ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for Si, O, and F; denoted as haTZ). Two monohydrate and six dihydrate minima have been identified for the SiF5-(H2O)n systems, whereas one monohydrate and five dihydrate minima have been identified for the SiF62-(H2O)n systems. Both monohydrated anions have a minimum in which the water molecule adopts a symmetric double ionic hydrogen bond (DIHB) motif with C2v symmetry. However, a second unique monohydrate minimum has been identified for SiF5- in which the water molecule adopts an asymmetric DIHB motif along the edge of the trigonal bipyramidal anion between one axial and one equatorial F atom. This Cs structure is more than 2 kcal mol-1 lower in energy than the C2v local minimum at the CCSD(T)/haTZ level of theory. While the interactions between the solvent and ionic solute are quite strong for the monohydrated anions (electronic dissociation energies of ≈12 and ≈24 kcal mol-1 for the SiF5-(H2O)1 and SiF62-(H2O)1 global minima, respectively), these values are nearly perfectly doubled for the dihydrates, with the lowest-energy SiF5-(H2O)2 and SiF62-(H2O)2 minima exhibiting dissociation energies of ≈24 and ≈47 kcal mol-1, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -100 cm-1 for the SiF5-(H2O)2 minimum and -300 cm-1 for the SiF62-(H2O)2 minimum relative to an isolated H2O molecule at the CCSD(T)/haTZ level of theory. This work also corrects the OH stretching frequency shifts for two dihydrate minima of PF6- that were previously erroneously reported ( J. Phys. Chem. A 2020, 124, 8744-8752, DOI: 10.1021/acs.jpca.0c06466).

Quantum symmetrization transition in superconducting sulfur hydride from quantum Monte Carlo and path integral molecular dynamics

AbstractWe study the structural phase transition, originally associated with the highest superconducting critical temperature Tc measured in high-pressure sulfur hydride. A quantitative description of its pressure dependence has been elusive for any ab initio theory attempted so far, raising questions on the actual mechanism leading to the maximum of Tc. Here, we estimate the critical pressure of the hydrogen bond symmetrization in the Im$$\bar{3}$$ 3 ¯ m structure, by combining density functional theory and quantum Monte Carlo simulations for electrons with path integral molecular dynamics for quantum nuclei. We find that the Tc maximum corresponds to pressures where local dipole moments dynamically form on the hydrogen sites, as precursors of the ferroelectric Im$$\bar{3}$$ 3 ¯ m-R3m transition, happening at lower pressures. For comparison, we also apply the self-consistent harmonic approximation, whose ferroelectric critical pressure lies in between the ferroelectric transition estimated by path integral molecular dynamics and the local dipole formation. Nuclear quantum effects play a major role in a significant reduction (≈50 GPa) of the classical ferroelectric transition pressure at 200 K and in a large isotope shift (≈25 GPa) upon hydrogen-to-deuterium substitution of the local dipole formation pressure, in agreement with the corresponding change in the Tc maximum location.

Breaking hydrogen-bonds in aqueous electrolytes towards highly reversible zinc-ion batteries

An electrolyte additive with uniquely symmetric structure breaks water's hydrogen bonds and modifies the primary solvation structure of Zn2+ in aqueous zinc batteries, enhancing battery performance and stability.

Competition between Solvent···Solvent and Solvent···Solute Interactions in the Microhydration of the Tetrafluoroborate Anion, BF4-(H2O)n=1,2,3,4.

This study systematically examines the interactions of the tetrafluoroborate anion (BF4-) with up to four water molecules (BF4-(H2O)n=1,2,3,4). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using a variety of density functional theory (DFT) methods (B3LYP, B3LYP-D3BJ, and M06-2X) and the MP2 ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for B, O, and F; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the CCSD(T) ab initio method and the haTZ basis set for the mono- and dihydrate (n = 1, 2) structures. The 2-body:Many-body (2b:Mb) technique, in which CCSD(T) computations capture the 1- and 2-body contributions to the interactions and MP2 computations recover all higher-order contributions, was used to extend these demanding computations to the tri- and tetrahydrate (n = 3, 4) systems. Four, five, and eight new stationary points have been identified for the di-, tri-, and tetrahydrate systems, respectively. The global minimum of the monohydrate adopts a symmetric double ionic hydrogen bond motif with C2v symmetry and an electronic dissociation energy of 13.17 kcal mol-1 at the CCSD(T)/haTZ level of theory. This strong solvent···solute interaction, however, competes with solute···solute interactions in the lowest-energy BF4-(H2O)n=2,3,4 minima that are not seen in the other di-, tri-, or tetrahydrate minima. The latter interactions help increase the 2b:Mb dissociation energies to more than 26, 41, and 51 kcal mol-1 for n = 2, 3, and 4, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -280 cm-1 relative to an isolated H2O molecule at the 2b:Mb/haTZ level of theory.

Magnetic Spin-Flop Transition of ε-FeOOH at 8 GPa

Mössbauer spectra of polycrystalline ε-57FeOOH were measured up to 14.6 GPa at room temperature using the energy-domain synchrotron 57Fe Mössbauer source at BL11XU in SPring-8. Based on measurements with/without an external magnetic field, ε-FeOOH appears to be antiferromagnetic at ambient pressure and temperature. On compression, the quadrupole shift 2ε shifted from negative to positive at ∼8 GPa, which is consistent with a magnetic spin rotation. Previous XRD studies did not detect any transitions at this pressure; thus we infer the change in 2ε to be the third pressure-induced transition without major structural change in addition to the hydrogen-bond symmetrization and electronic spin transition.

Equation of State and Spin Crossover of (Al, Fe)‐Phase H

AbstractThe transport of hydrogen into Earth's deep interior may have an impact on lower mantle dynamics as well as on the seismic signature of subducted material. Due to the stability of the hydrous phases δ‐AlOOH (delta phase), MgSiO2(OH)2 (phase H), and ε‐FeOOH at high temperatures and pressures, their solid solutions may transport significant amounts of hydrogen as deep as the core‐mantle boundary. We have constrained the equation of state, including the effects of a spin crossover in the Fe3+ atoms, of (Al, Fe)‐phase H: Al0.84Fe3+ 0.07Mg0.02Si0.06OOH, using powder X‐ray diffraction measurements to 125 GPa, supported by synchrotron Mössbauer spectroscopy measurements on (Al, Fe)‐phase H and δ‐(Al, Fe)OOH. The changes in spin state of Fe3+ in (Al, Fe)‐phase H results in a significant decrease in bulk sound velocity and occurs over a different pressure range (48–62 GPa) compared with δ‐(Al, Fe)OOH (32–40 GPa). Changes in axial compressibilities indicate a decrease in the compressibility of hydrogen bonds in (Al, Fe)‐phase H near 30 GPa, which may be associated with hydrogen bond symmetrization. The formation of (Al, Fe)‐phase H in subducted oceanic crust may contribute to scattering of seismic waves in the mid‐lower mantle (∼1,100–1,550 km). Accumulation of 1–4 wt.% (Al, Fe)‐phase H could reproduce some of the seismic signatures of large, low seismic‐velocity provinces. Our results suggest that changes in the electronic structure of phases in the (δ‐AlOOH)‐(MgSiO2(OH)2)‐(ε‐FeOOH) solid solution are sensitive to composition and that the presence of these phases in subducted oceanic crust could be seismically detectable throughout the lower mantle.

Hydrogen bond symmetrization and high-spin to low-spin transition of ε-FeOOH at the pressure of Earth’s lower mantle

Abstract We focus on the ferric end-member of phase H: ε-FeOOH using density functional theory at the PBEsol+U level. At 300 K, we find that ε-FeOOH undergoes a hydrogen bond symmetrization at 37 GPa and a sharp high-spin to low-spin transition at 45 GPa. We find excellent agreement with experimental measurements of the equation of state, lattice parameters, atomic positions, vibrational frequencies, and optical properties as related to the band gap, which we find to be finite and small, decreasing with pressure. The hydrogen bond symmetrization transition is neither first-nor second-order, with no discontinuity in volume or any of the elastic moduli. Computed IR and Raman frequencies and intensities show that vibrational spectroscopy may provide the best opportunity for locating the hydrogen bond symmetrization transition experimentally. We find that ε-FeOOH is highly anisotropic in both longitudinal- and shear-wave velocities at all pressures, with the shear wave velocity varying with propagation and polarization direction by as much as 24% at zero pressure and 43% at 46 GPa. The shear and bulk elastic moduli increase by 18% across the high-spin to low-spin transition.

Polyether phases of formic acid revealed under high pressure.

Evolutionary structure searches reveal that at 20 GPa, formic acid (FA) transforms to an orthorhombic ether-chain polymer phase, ruling out hydrogen bond symmetrization, followed by a novel crowned cyclic-ether tetragonal phase above 60 GPa. Emergence of characteristic polyether and new OH stretching modes in infrared experiments validate the findings. Resemblance of polymer chain with a cosmopolymer polyoxymethylene, shows an engrossing multifaceted evolution of FA.

Calculated Elasticity of Al-Bearing Phase D

Using first-principles calculations, this study evaluates the structure, equation of state, and elasticity of three compositions of phase D up to 75 GPa: (1) the magnesium endmember [MgSi2O4(OH)2], (2) the aluminum endmember [Al2SiO4(OH)2], and (3) phase D with 50% Al-substitution [AlMg0.5Si1.5O4(OH)2]. We find that the Mg-endmember undergoes hydrogen-bond symmetrization and that this symmetrization is linked to a 22% increase in the bulk modulus of phase D, in agreement with previous studies. Al2SiO4(OH)2 also undergoes hydrogen-bond symmetrization, but the concomitant increase in bulk modulus is only 13%—a significant departure from the 22% increase of the Mg-endmember. Additionally, Al-endmember phase D is denser (2%–6%), less compressible (6%–25%), and has faster compressional (6%–12%) and shear velocities (12%–15%) relative to its Mg-endmember counterpart. Finally, we investigated the properties of phase D with 50% Al-substitution [AlMg0.5Si1.5O4(OH)2], and found that the hydrogen-bond symmetrization, equation of state parameters, and elastic constants of this tie-line composition cannot be accurately modeled by interpolating the properties of the Mg- and Al-endmembers.

Superconducting LaP2H2 with graphenelike phosphorus layers

Novel structural building blocks in compounds could induce interesting physical and chemical properties. Although phosphorus tends to form very different motifs, the existence of lone pair electrons has always prevented the formation of graphenelike structures. Here, the application of first-principles swarm structural calculations has allowed us to predict the stability of pressure-induced hexagonal LaP2H2 containing graphenelike phosphorus, which derives from the trigonal bipyramid configuration of P atoms regulated by symmetric hydrogen bonds. LaP2 in LaP2H2 has the same configuration as MgB2, and P and H atoms form a three-dimensional framework as H3S. Interestingly, LaP2H2 shows a superconductivity dominated by the graphenelike phosphorus layer and its coupling with La atoms. On the other hand, LaP2H2 is not only superconducting at a lower pressure than the H-rich LaPH6, but it also shows a superconducting transition temperature three times higher. Our work provides an example which extends the landscape of conventional superconductors at lower pressures.

Potential energy barrier for proton transfer in compressed benzoic acid.

Benzoic acid (BA) is a model system for studying proton transfer (PT) reactions. The properties of solid BA subject to high pressure (exceeding 1 kbar = 0.1 GPa) are of particular interest due to the possibility of compression-tuning of the PT barrier. Here we present simulations aimed at evaluating the value of this barrier in solid BA in the 1 atm – 15 GPa pressure range. We find that pressure-induced shortening of O⋯O contacts within the BA dimers leads to a decrease in the PT barrier, and subsequent symmetrization of the hydrogen bond. However, this effect is obtained only after taking into account zero-point energy (ZPE) differences between BA tautomers and the transition state. The obtained results shed light on previous experiments on compressed benzoic acid, and indicate that a common scaling behavior with respect to the O⋯O distance might be applicable for hydrogen-bond symmetrization in both organic and inorganic systems.

Absence of proton tunneling during the hydrogen-bond symmetrization in δ−AlOOH

$\delta$-AlOOH is of significant crystallochemical interest due to a subtle structural transition near 10 GPa from a $P2_1nm$ to a $Pnnm$ structure, the nature and origin of hydrogen disorder, the symmetrization of the O-H$\cdots$O hydrogen bond and their interplay. We perform a series of density functional theory based simulations in combination with high-pressure nuclear magnetic resonance experiments on $\delta$-AlOOH up to 40 GPa with the goal to better characterize the hydrogen potential and therefore the nature of hydrogen disorder. Simulations predict a phase transition in agreement with our nuclear magnetic resonance experiments at $10-11$ GPa and hydrogen bond symmetrization at $14.7$ GPa. Calculated hydrogen potentials do not show any double-well character and there is no evidence for proton tunneling in our nuclear magnetic resonance data.

Pressure-Induced Synthesis and Properties of an H2S-H2Se-H2 Molecular Alloy.

The chalcogens are known to react with one another to form interchalcogens, which exhibit a diverse range of bonding and conductive behavior due to the difference in electronegativity between the group members. Through a series of high-pressure diamond anvil experiments combined with density functional theory calculations, we report the synthesis of an S-Se hydride. At pressures above 4 GPa we observe the formation of a single solid composed of both H2Se and H2S molecular units. Further compression in a hydrogen medium leads to the formation of an alloyed compound (H2SxSe1-x)2H2, after which there is a sequence of pressure-induced phase transitions associated with the arrested rotation of molecules. At pressures above 50 GPa, there is a symmetrization of hydrogen bonds concomitantly with a closing band gap and increased reflectivity of the compound, indicative of a transition to a metallic state.

Nb6Cl12(CH3OH)4(OCH3)2] ⋅ DABCO ⋅ 1.66 CH2Cl2: Cluster Units Wrapped in a Network of Hydrogen Bonds

AbstractAn easy‐to‐do synthesis for the hexanuclear niobium cluster compound [Nb6Cl12(CH3OH)4(OCH3)2] ⋅ DABCO ⋅ 1.66 CH2Cl2 has been developed. An one‐pot reaction between the cluster precursor [Nb6Cl14(H2O)4] ⋅ 4H2O and methanol with the addition of DABCO leads to the crystallization of the title compound in high yield within a few minutes. The single‐crystal X‐ray structure of this cluster compound has been determined. Very strong, nearly symmetric intercluster hydrogen bonds Nb6−MeO⋅⋅⋅H⋅⋅⋅OMe−Nb6 are present between the cluster units. A bridging co‐crystalline DABCO molecule is also involved in a three‐dimensional hydrogen‐bonding network.

Quantum Symmetrization of Hydrogen Bonds in Ice

Quantum Symmetrization of Hydrogen Bonds in Ice

High-pressure vibrational spectra of humite-group minerals: Fluorine effect on thermodynamic properties and hydrogen bonds

The humite-group minerals could be potential water carriers as the dehydration products of serpentine in subduction slabs. Fluorine is an important element incorporated into the crystal structures through the substitution mechanism of OH− = F−, and has significant impacts on the thermodynamic properties. We conducted in situ high-pressure Raman and Fourier transform infrared (FTIR) measurements on natural F-bearing clinohumite, humite and norbergite samples, up to approximately 19GPa. The isothermal mode Grüneisen and intrinsic anharmonic parameters were calculated for both the lattice and OH-stretching modes, and the anharmonic contribution to the heat capacities was further evaluated for these hydrous phases. As compared with the OH-clinohumite sample, the anharmonic contribution is larger in magnitude for the F-bearing clinohumite, while such anharmonic contributions are even more significant in the humite and norbergite samples with higher fluorine concentrations. Both the Raman and FTIR spectra indicate that the OH-stretching modes above 3450 cm−1, which are controlled by the neighboring HH, systematically show positive pressure dependences. On the other hand, the OH bands below 3450 cm−1 shift to lower frequencies at elevated pressure, due to F− or Ti4+ substitutions. The F− substitution of OH− group could alleviate local HH repulsion effect, and create new H sites with shorter hydrogen bond lengths. These hydrogen bond would get shorter (negative pressure dependences), and are likely to participate in the hydrogen bond symmetrization at high pressure.

Phase transitions in ε-FeOOH at high pressure and ambient temperature

AbstractConstraining the accommodation, distribution, and circulation of hydrogen in the Earth's interior is vital to our broader understanding of the deep Earth due to the significant influence of hydrogen on the material and rheological properties of minerals. Recently, a great deal of attention has been paid to the high-pressure polymorphs of FeOOH (space groups P21nm and Pnnm). These structures potentially form a hydrogen-bearing solid solution with AlOOH and phase H (MgSiO4H2) that may transport water (OH–) deep into the Earth's lower mantle. Additionally, the pyrite-type polymorph (space group Pa3 of FeOOH), and its potential dehydration have been linked to phenomena as diverse as the introduction of hydrogen into the outer core (Nishi et al. 2017), the formation of ultralow-velocity zones (ULVZs) (Liu et al. 2017), and the Great Oxidation Event (Hu et al. 2016). In this study, the high-pressure evolution of FeOOH was re-evaluated up to ~75 GPa using a combination of synchrotron-based X-ray diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), and optical absorption spectroscopy. Based on these measurements, we report three principal findings: (1) pressure-induced changes in hydrogen bonding (proton disordering or hydrogen bond symmetrization) occur at substantially lower pressures in ε-FeOOH than previously reported and are unlikely to be linked to the high-spin to low-spin transition; (2) ε-FeOOH undergoes a 10% volume collapse coincident with an isostructural Pnnm → Pnnm transition at approximately 45 GPa; and (3) a pressure-induced band gap reduction is observed in FeOOH at pressures consistent with the previously reported spin transition (40 to 50 GPa).

Electrical conductivity of diaspore, δ-AlOOH and ε-FeOOH

AbstractElectrical conductivities of diaspore (α-AlOOH), δ-AlOOH, and ε-FeOOH were measured by impedance spectroscopy with a frequency range from 10–1 to 106 Hz at pressures from 8 to 20 GPa and temperatures from 500 to 1200 K, well below the dehydration temperatures of these phases at the relevant pressures. For diaspore, the relationship between electrical conductivity and reciprocal temperature can be well fitted by the Arrhenius formula:σ = σ 0 exp ⁡ [ − ( Δ E + P Δ V ) k T ] ,where σ0 is the pre-exponential factor, ∆E is the activation energy, and ∆V is activation volume of 56.0 ± 1.2 S/m, 0.55 ± 0.02 eV, and 1.68 ± 0.12 cm3/mol, respectively. The electrical conductivity of diaspore decreases with increasing pressure ranging from 8 to 12 GPa by a half order of magnitude, whereas the conductivity becomes almost constant in a pressure range above 12 GPa. δ-AlOOH and ε-FeOOH show one and two orders of magnitude higher electrical conductivity than diaspore. Electrical conductivities of δ-AlOOH and ε-FeOOH, which have isostructural CaCl2-type hydroxide structure, show the nearly identical activation enthalpies (0.38 ± 0.01, 0.33 ± 0.05 eV), which are relatively lower than that of diaspore. The dominant conduction mechanism of AlOOH phases can be regarded as proton conduction. The conductivity difference between diaspore and δ-AlOOH attributes to result in the different O1H bond lengths of each phase. The reduction of O1H bond length with increasing pressure could enhance the proton migration by reducing the potential barrier, thereby raising the electrical conductivity. Small polaron conduction may contribute to the conductivity of ε-FeOOH to generate higher conductivity than δ-AlOOH. Furthermore, hydrogen bond symmetrization will also play an important role in the conductivity discrepancy of these hydrous minerals with CaCl2-type hydroxide structure. For subducted sedimentary rocks, polymorphs of AlOOH and FeOOH are representative hydrous phases. Al-rich sediments show conductivity reduction with increasing depth until phase transformation occurs because diaspore represents negative pressure dependence of conductivity. After transformation to δ-AlOOH, the conductivity will jump up around 18 GPa. If ε-FeOOH is stable above 5 GPa in an iron-rich lithology, such as banded iron formation (BIF), a high conductivity zone with positive pressure dependence could be observed to the deep lower mantle.

The caesium phosphates Cs3(H1.5PO4)2(H2O)2, Cs3(H1.5PO4)2, Cs4P2O7(H2O)4, and CsPO3

The caesium phosphates Cs3(H1.5PO4)2(H2O)2 and Cs3(H1.5PO4)2 were obtained from aqueous solutions, and Cs4P2O7(H2O)4 and CsPO3 from solid state reactions, respectively. Cs3(H1.5PO4)2, Cs4P2O7(H2O)4, and CsPO3 were fully structurally characterized for the first time on basis of single-crystal X-ray diffraction data recorded at − 173 °C. Monoclinic Cs3(H1.5PO4)2 (Z = 2, C2/m) represents a new structure type and comprises hydrogen phosphate groups involved in the formation of a strong non-symmetrical hydrogen bond (accompanied by a disordered H atom over a twofold rotation axis) and a very strong symmetric hydrogen bond (with the H atom situated on an inversion centre) with symmetry-related neighbouring anions. Triclinic Cs4P2O7(H2O)4 (Z = 2, Pbar{1}) crystallizes also in a new structure type and is represented by a diphosphate group with a P–O–P bridging angle of 128.5°. Although H atoms of the water molecules were not modelled, O···O distances point to hydrogen bonds of medium strengths in the crystal structure. CsPO3 is monoclinic (Z = 4, P21/n) and belongs to the family of catena-polyphosphates (MPO3)n with a repetition period of 2. It is isotypic with the room-temperature modification of RbPO3. The crystal structure of Cs3(H1.5PO4)2(H2O)2 was re-evaluated on the basis of single-crystal X-ray diffraction data at − 173 °C, revealing that two adjacent hydrogen phosphate anions are connected by a very strong and non-symmetrical hydrogen bond, in contrast to the previously described symmetrical bonding situation derived from room temperature X-ray diffraction data. In the four title crystal structures, coordination numbers of the caesium cations range from 7 to 12.Graphic abstract