OH-stretching Region Research Articles

Tracking the Ti4+ substitution in phlogopite by spectroscopic imaging: A tool for unravelling the growth of micas at HP-HT conditions

Phlogopite solid-solutions have a wide pressure–temperature (P-T) stability field and are ubiquitous in a wide variety of geological settings, from deep lithosphere magmatic environments to upper crust metamorphic domains. Phlogopite composition represents therefore a valuable physical–chemical archive and may provide important information regarding its crystallization and the petrogenesis of the host-rock. In this paper we examine the phlogopite phenocrysts from the well-known Fort Regent mica-bearing lamprophyre minette from St. Helier (Island of Jersey, UK). Phlogopite phenocrystals from lamprophyres generally show normal-step and continuous compositional zoning, however those from the Fort Regent minette show a peculiar texture characterized by dark brown high-Ti (average TiO2 ≈ 8.5 wt.%) cores enveloped by euhedral low- to mid-amplitude zonation due to oscillatory contents in Ti, Fe and Mg. Thermo-barometry modelling based on biotite-only composition yields relatively high P-T estimates (T ≈ 970 ± 54 °C at P ≈ 0.73 ± 0.13 GPa) for cores whereas lower values (T ≈ 790 ± 54 °C at P ≈ 0.29 ± 0.13 GPa) are obtained for the outer rims. Comparable temperatures (T ≈ 1075 ± 54 °C) but extremely high and anomalous pressure values (P ≈ 1.82 ± 0.13 GPa) are obtained for the yellowish inner rims. The combination of electron micro probe (EMP) analysis and single-crystal infra-red (FTIR) imaging in the OH-stretching region shows that the exceptional and oscillatory Ti contents are due to the Ti-vacancy substitution, typical of crystallization and growth processes of HP/HT environments. Raman imaging provides additional insight for this process, confirming the dominant dioctahedral nature for the Ti-Fe-rich cores and outer rims. Interpretation of thermobaric estimates obtained from the phlogopite composition-only model, based on the fine-scale compositional evolution, shows that pressure–temperature values from low-Ti high-Mg domains should be carefully evaluated because the substitution mechanisms during the dark mica growth are not univocally related to pressure–temperature variation of the crystallizing environment. Our results demonstrate how a multidisciplinary approach based on the combination of chemical investigations and vibrational spectroscopies could represent a valuable tool to evaluate pressure–temperature estimates from biotite composition-only thermo-barometry models and therefore to correctly unravel HP/HT petrogenetic processes at a very fine scale.

Predicting the Raman Spectra of Liquid Water with a Monomer-Field Model.

The Raman spectrum of liquid water is quite complex, reflecting its strong sensitivity to the local environment of the individual waters. The OH-stretch region of the spectrum, which captures the influence of hydrogen bonding, has only just begun to be unraveled. Here we develop a model for predicting the Raman spectra of the OH-stretch region by considering how local electric fields distort the energy surface of each water monomer. We find that our model is capable of reproducing the bimodal nature of the main peak, with the shoulder at 3250 cm-1 resulting almost entirely from Fermi resonance. Furthermore, we capture the temperature and polarization dependence of the shoulder, which has proven to be difficult to obtain with previous methods, and analyze the origin of this dependence. We expect our model to be generally useful for understanding and predicting how Raman spectra change under different conditions and with different probe reporters beyond water.

Investigating Intramolecular H Atom Transfer Dynamics in β-Diketones with Ultrafast Infrared Spectroscopies and Theoretical Modeling.

The vibrational signatures and ultrafast dynamics of the intramolecular H-bond in a series of β-diketones are investigated with 2D IR spectroscopy and computational modeling. The chosen β-diketones exhibit a range of H atom donor-acceptor distances and asymmetry along the H atom transfer coordinate that tunes the intramolecular H-bond strength. The species with the strongest H-bonds are calculated to have very soft H atom potentials, resulting in highly red-shifted OH stretch fundamental frequencies and dislocation of the H atom upon vibrational excitation. These soft potentials lead to significant coupling to the other normal mode coordinates and give rise to the very broad vibrational signatures observed experimentally. The 2D IR spectra in both the OH and OD stretch regions of the light and deuterated isotopologues reveal broadened and long-lived ground-state bleach signatures of the vibrationally hot molecules. Polarization-sensitive transient absorption measurements in the OH and OD stretch regions reveal notable isotopic differences in orientational dynamics. Orientational relaxation was measured to occur on ∼600 fs and ∼2 ps time scales for the light and deuterated isotopologues, respectively. The orientational dynamics are interpreted in terms of activated H/D atom transfer events driven by collective intramolecular structural rearrangements.

Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region.

We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.

Ion-Specific Water–Macromolecule Interactions at the Air/Aqueous Interface: An Insight into Hofmeister Effect

The specificity of ions in inducing conformational changes in macromolecules is introduced as the Hofmeister series; however, the detailed underlying mechanism is not comprehensible yet. We utilized surface-specific sum frequency generation (SFG) vibrational spectroscopy to explore the Hofmeister effect at the air/polyvinylpyrrolidone (PVP)/water interface. The spectral signature observed from the ssp polarization scheme reveals ion-specific ordering of water molecules following the Hofmeister series attributed to the ion-macromolecule interactions. Along with this, the presence of ions does not reflect any significant influence on the structure of the PVP macromolecule. However, the ppp-SFG spectra in the CH-stretch region reveal the impact of ions on the orientation angle of vinyl chain CH2-groups, which follows the Hofmeister series: SO42- > Cl- > NO3- > Br- > ClO4- > SCN-. The minimal orientation angle of CH2-groups indicates significant reordering in PVP vinyl chains in the presence of chaotropic anions ClO4-, and SCN-. The observation is attributed to the ion-specific water-macromolecule interactions at the air/aqueous interface. It is compelling to observe the signature of spectral blue shifts in the OH-stretch region in the ppp configuration in the presence of chaotropic anions. The origin of spectral blue shifts has been ascribed to the existence of weaker interactions between the interfacial water molecules and the backbone CH- and CH2-moieties of the PVP macromolecules. The ion-specific modulation in water-macromolecule interactions is endorsed by the relative propensity of anion's adsorption toward the air/aqueous interface. The experimental findings highlight the existence and cooperative participation of ion-specific water-macromolecule interactions in the mechanism of the Hofmeister effect, along with the illustrious ion-water and ion-macromolecule interactions.

Elucidation of the pH-Dependent Electric Double Layer Structure at the Silica/Water Interface Using Heterodyne-Detected Vibrational Sum Frequency Generation Spectroscopy

The silica/water interface is one of the most abundant charged interfaces in natural environments, and the elucidation of the water structure at the silica/water interface is essential. In the present study, we measured the interface-selective vibrational (χ(2)) spectra in the OH stretch region of the silica/water interface in a wide pH range of pH 2.0-12.0 while changing the salt concentration by heterodyne-detected vibrational sum-frequency generation spectroscopy. With the help of singular value decomposition analysis, it is shown that the imaginary part of the χ(2) (Imχ(2)) spectra can be decomposed into the spectra of the diffuse Gouy-Chapman layer (DL) and the compact Stern layer (SL), which enables us to quantitatively analyze the spectra of DL and SL separately. The salt-concentration dependence of the DL spectra at different pH values is analyzed using the modified Gouy-Chapman theory, and the pH-dependent surface charge density and the pKa value (4.8 ± 0.2) of the silica/water interface are evaluated. Furthermore, it is found that the pH-dependent change of the SL spectra is quantitatively explained by three spectral components that represent the three characteristic water species appearing in different pH regions in SL. The quantitative understanding obtained from the analysis of each spectral component in the Imχ(2) spectra provides a clear molecular-level picture of the electric double layer at the silica/water interface.

Optimal Raman-scattering signal for estimating the Fe3+ content on the clinozoisite–epidote join

Abstract. To provide a tool for fast estimation of the Fe3+ content in Ca2(Al, Fe3+)3Si3O12(OH) epidote grains, including in thin sections and crude-rock samples, we applied Raman spectroscopy to 33 areas from 15 natural samples with Fe3+ ranging from 0.22 to 1.13 atoms per formula unit (apfu), the chemistry of which was independently determined by wavelength-dispersive electron microprobe analysis (WD-EPMA). The Raman spectra were collected from the very areas subjected to WD-EPMA. We have analysed both the OH-stretching region (3215–3615 cm−1) and the spectral range generated by the framework vibrations (15–1215 cm−1). Similarly to the IR spectra, the Raman peaks in the OH-stretching region shift toward higher wavenumbers with increasing Fe. However, the quantification of Fe3+ based on OH-stretching Raman peaks can be hindered by the multicomponent overlapping and significant intensity variations with the crystal orientation. Among the Raman signals generated by framework vibrations, the position of four peaks (near 250, 570, 600, and 1090 cm−1) exhibit a steady linear regression with the increase in Fe content (in apfu). However, the peak near 250 cm−1 attributed to MO6 vibrations also depends on the crystal orientation and therefore is not always well resolved, which worsens the accuracy in Fe-content determination based on its position. The peaks near 570, 600, and 1090 cm−1 arise from Si2O7 vibrational modes, and although their intensities also vary with the crystal orientation, all three signals are well resolved in a random orientation. However, among the three Si2O7-related signals, the 570 cm−1 peak is the sharpest (peak width <10 cm−1) and is easily recognized as a separate peak. Hence, we propose to use the position of this peak as a highly reliable parameter to estimate the Fe content, via the linear trend given as ω570=577.1(3)-12.7(4)x, where ω is the wavenumber (cm−1) and x is Fe content (apfu), with accuracy ± 0.04 Fe3+ apfu. The peaks near 600 and 1090 cm−1 may be complementarily used for the Fe estimate, based on the following relations: ω600=611.6(2)-13.8(4)x and ω1090=1098.8(3)-13.5(5)x. Analyses of the effect of Sr as a substitution for Ca and Cr at the octahedral sites indicate that contents of Sr <0.12 apfu do not interfere with the quantification of Fe via the ω570 (x) relation, whereas Cr >0.16 apfu leads to overestimation of Fe; Cr presence can be recognized however by the broadening of the peaks near 95 and 250 cm−1.

Transferability of vibrational spectroscopic map from TIP4P to TIP4P-like water models.

We computed the IR, Raman, and sum frequency generation spectra of water in the OH-stretch region by employing the quantum/classical mixed approach that consists of a vibrational spectroscopic map and molecular dynamics (MD) simulation. We carried out the MD simulation with the TIP4P, TIP4P/2005, and TIP4P/Ice models and applied the map designed for TIP4P by Skinner et al. to each MD trajectory. Although the map is not tuned for TIP4P-like models, TIP4P/2005 and TIP4P/Ice provide the best reproduction of the experimental vibrational spectra of liquid water and crystalline ice, respectively. This result demonstrates the transferability of the map from TIP4P to TIP4P/2005 and TIP4P/Ice, meaning that one can choose an appropriate TIP4P-like model to calculate the vibrational spectra of an aqueous system without rebuilding the map.

High-resolution rovibrational spectroscopy of trans-formic acid in the v1 OH stretching fundamental: Dark state coupling and evidence for charge delocalization dynamics

High-resolution infrared (IR) reduced-Doppler absorption spectra of jet-cooled gas phase trans-formic acid in the v1 OH stretching fundamental region are reported for the first time, obtained by supersonically expanding trans-formic acid/Ar mixtures through a slit jet nozzle source and rotationally cooling to Trot ≈ 10.9(5) K, with absorption signals recorded by high-resolution difference-frequency IR absorption spectroscopy. Two a/b-type rovibrational bands of comparable intensity, one ∼10-fold weaker b-type band, and one ∼6-fold weaker a-type band are observed, with vibrational band origins at 3570.493(5), 3566.793(5), 3560.032(9), and 3534.6869(2) cm−1, respectively. Based on previous Raman jet spectroscopic work by Nejad and Sibert [A. Nejad, E.L. Sibert III, The Raman jet spectrum of trans-formic acid and its deuterated isotopologs: Combining theory and experiment to extend the vibrational database, J. Chem. Phys. 154(6) (2021) 064301.], these four rovibrational bands have been assigned to v1, (v2 + v7), (v6 + 2v7 + 2v9), and 2v3, respectively. Specifically, two of the three upper dark states (2171 (a′) and 617292 (a′)) are close enough to the “bright” 11 (a′) state to facilitate strong anharmonic resonance interactions, which results in intensity mixing into the two zero-order bands that would otherwise be “dark”. Furthermore, our high-resolution spectral analysis reveals that there are local rotational crossings between these zero-order 11 and 2171 states resulting in extra lines (i.e., some upper levels in the nominally v1 band have majority zero-order 2171 state character). This motivates development of a 3 coupled state (11, 2171, and 617292) picture to aid in the spectral analysis, which is able to match all 3 observed band origins and relative band intensities, as well as indicate the necessity of multistate (> 2) coupling. Though limited by range of J and Ka levels (J' ≤ 9 and Ka' ≤ 3) populated at supersonic jet temperatures, this work offers first precision spectroscopic analysis of trans-formic acid in the v1 OH stretch region, which should aid in assignment of the more complete yet highly congested room temperature FTIR spectra [D. Hurtmans, F. Herregodts, M. Herman, J. Liévin, A. Campargue, A. Garnache, A. Kachanov, Spectroscopic and ab initio investigation of the νOH overtone excitation in trans-formic acid, J. Chem. Phys. 113(4) (2000) 1535–1545.]. Finally, and in sharp contrast to the spectral complexity in the three predominantly b-type bands, the lone a-type 2v3 rovibrational band at 3534.6869(2) cm−1 is well described by a simple, rigid asymmetric top Hamiltonian.

Vibrational Sum Frequency Generation Spectra of Water-Vapor Interfaces Covered by Alcohols: Effects of Surface Coverage and Coupling between Oscillators.

The present study deals with the effects of varying coverage of water surface by alcohols on the vibrational sum frequency generation (VSFG) spectrum of interfacial water. We have considered two different alcohols: Tertiary butyl alcohol (TBA) whose alkyl part is fully branched and stearyl alcohol (STA) which has a long linear alkyl chain with larger hydrophobic surface area than that of TBA. With increase of the alcohol concentration, the hydrogen bonded OH stretch region of the VSFG spectrum is found to change following a regular trend for the STA-water system, whereas non-monotonic variation of the VSFG spectrum is observed for the TBA-water system which can be correlated with the presence of very different interactions of TBA molecules at different concentrations. On increasing the concentration of TBA, the hydrophobic groups get more tilted towards the water phase and significant hydrophobic interactions are introduced at higher concentrations. Whereas, for STA, there is a gradual increase in the hydrophilic interaction. Because of stacking interactions between the long chain alkyl groups, the hydrophobic parts stay outward from the water phase at higher concentrations and a regular change in the VSFG spectrum is observed. We have also presented a computationally efficient scheme to calculate the VSFG spectrum of interfacial systems for coupled oscillators which is expected to be beneficial for the treatment of coupling where the interfacial system size is inherently large.

Hydrothermal genesis and growth of the banded agates from the Allumiere-Tolfa volcanic district (Latium, Italy)

In this work, we studied the hydrothermal agates from the Neogene–Quaternary volcanic district of Allumiere-Tolfa, north-west of Rome (Latium, Italy) using a combination of micro-textural, spectroscopic, and geochemical data. The examined sample consists of (1) an outer cristobalite layer deposited during the early stages of growth, (2) a sequence of chalcedonic bands (including i.e., length-fast, zebraic, and minor length-slow chalcedony) with variable moganite content (up to ca. 48 wt%), (3) an inner layer of terminated hyaline quartz crystals. The textures of the various SiO2 phases and their trace element content (Al, Li, B, Ti, Ga, Ge, As), as well as the presence of mineral inclusions (i.e., Fe-oxides and sulfates), is the result of physicochemical fluctuations of SiO2-bearing fluids. Positive correlation between Al and Li, low Al/Li ratio, and low Ti in hyaline quartz points to low-temperature hydrothermal environment. Local enrichment of B and As in chalcedony-rich layers are attributed to pH fluctuations. Analysis of the FT-IR spectra in the principal OH-stretching region (2750–3750 cm−1) shows that the silanol and molecular water signals are directly proportional. Strikingly, combined Raman and FT-IR spectroscopy on the chalcedonic bands reveals an anticorrelation between the moganite content and total water (SiOH + molH2O) signal. The moganite content is compatible with magmatic-hydrothermal sulfate/alkaline fluids at a temperature of 100–200 °C, whereas the boron-rich chalcedony can be favored by neutral/acidic conditions. The final Bambauer quartz growth lamellae testifies diluted SiO2-bearing solutions at lower temperature. These findings suggest a genetic scenario dominated by pH fluctuations in the circulating hydrothermal fluid.

Stabilization Effects Induced by Trehalose on Creatine Aqueous Solutions Investigated by Infrared Spectroscopy.

Creatine is a very popular amino acid widely utilized in the sports world due to its functions mainly related to muscle building and increasing performance. The present work investigates the behavior of creatine aqueous solutions and of creatine aqueous in the presence of trehalose as a function of time changes by means of Infrared spectroscopy. Infrared spectra have been gathered and studied over time for both the full spectrum and the intramolecular OH-stretching region for the two mixtures. This latter region was studied more specifically using a cutting-edge technique called Spectral Distance (SD). From this analysis of the spectral features of the investigated samples, it emerges that trehalose has a significant stabilizing effect on creatine aqueous solutions.

Picosecond infrared laser driven sample delivery for simultaneous liquid-phase and gas-phase electron diffraction studies.

Here, we report on a new approach based on laser driven molecular beams that provides simultaneously nanoscale liquid droplets and gas-phase sample delivery for femtosecond electron diffraction studies. The method relies on Picosecond InfraRed Laser (PIRL) excitation of vibrational modes to strongly drive phase transitions under energy confinement by a mechanism referred to as Desorption by Impulsive Vibrational Excitation (DIVE). This approach is demonstrated using glycerol as the medium with selective excitation of the OH stretch region for energy deposition. The resulting plume was imaged with both an ultrafast electron gun and a pulsed bright-field optical microscope to characterize the sample source simultaneously under the same conditions with time synchronization equivalent to sub-micrometer spatial resolution in imaging the plume dynamics. The ablation front gives the expected isolated gas phase, whereas the trailing edge of the plume is found to consist of nanoscale liquid droplets to thin films depending on the excitation conditions. Thus, it is possible by adjusting the timing to go continuously from probing gas phase to solution phase dynamics in a single experiment with 100% hit rates and very low sample consumption (<100 nl per diffraction image). This approach will be particularly interesting for biomolecules that are susceptible to denaturation in turbulent flow, whereas PIRL–DIVE has been shown to inject molecules as large as proteins into the gas phase fully intact. This method opens the door as a general approach to atomically resolving solution phase chemistry as well as conformational dynamics of large molecular systems and allow separation of the solvent coordinate on the dynamics of interest.

Ultrafast vibrational dynamics of the free OD at the air/water interface: Negligible isotopic dilution effect but large isotope substitution effect.

Vibrational relaxation dynamics of the OH stretch of water at the air/water interface has been a subject of intensive research, facilitated by recent developments in ultrafast interface-selective nonlinear spectroscopy. However, a reliable determination of the vibrational relaxation dynamics in the OD stretch region at the air/D2O interface has not been yet achieved. Here, we report a study of the vibrational relaxation of the free OD carried out by time-resolved heterodyne-detected vibrational sum frequency generation spectroscopy. The results obtained with the aid of singular value decomposition analysis indicate that the vibrational relaxation (T1) time of the free OD at the air/D2O interface and air/isotopically diluted water (HOD-H2O) interfaces show no detectable isotopic dilution effect within the experimental error, as in the case of the free OH in the OH stretch region. Thus, it is concluded that the relaxation of the excited free OH/OD predominantly proceeds with their reorientation, negating a major contribution of the intramolecular energy transfer. It is also shown that the T1 time of the free OD is substantially longer than that of the free OH, further supporting the reorientation relaxation mechanism. The large difference in the T1 time between the free OD and the free OH (factor of ∼2) may indicate the nuclear quantum effect on the diffusive reorientation of the free OD/OH because this difference is significantly larger than the value expected for a classical rotational motion.

Crystallographic orientation mapping of lizardite serpentinite by Raman spectroscopy

Abstract. The serpentine mineral lizardite displays strong Raman anisotropy in the OH-stretching region, resulting in significant wavenumber shifts (up to ca. 14.5 cm−1) that depend on the orientation of the impinging excitation laser relative to the crystallographic axes. We quantified the relationship between crystallographic orientation and Raman wavenumber using well-characterised samples of Monte Fico lizardite by applying Raman spectroscopy and electron backscatter diffraction (EBSD) mapping on thin sections of polycrystalline samples and grain mounts of selected single crystals, as well as by a spindle stage Raman study of an oriented cylinder drilled from a single crystal. We demonstrate that the main band in the OH-stretching region undergoes a systematic shift that depends on the inclination of the c-axis of the lizardite crystal. The data are used to derive an empirical relationship between the position of this main band and the c-axis inclination of a measured lizardite crystal: y=14.5cos 4 (0.013x+0.02)+(3670±1), where y is the inclination of the c-axis with respect to the normal vector (in degrees), and x is the main band position (wavenumber in cm −1) in the OH-stretching region. This new method provides a simple and cost-effective technique for measuring and quantifying the crystallographic orientation of lizardite-bearing serpentinite fault rocks, which can be difficult to achieve using EBSD alone. In addition to the samples used to determine the above empirical relationship, we demonstrate the applicability of the technique by mapping the orientations of lizardite in a more complex sample of deformed serpentinite from Elba Island, Italy.

Vibrational Dynamics of the Intramolecular H-Bond in Acetylacetone Investigated with Transient and 2D IR Spectroscopy.

Acetylacetone (AcAc) has proven to be a fruitful but highly challenging model system for the experimental and computational interrogation of strong intramolecular hydrogen bonds. Key questions remain, however, regarding the identity of the minimum-energy structure of AcAc and the dynamics of intramolecular proton transfer. Here, we investigate the OH/OD stretch and bend regions of the enol tautomer of AcAc and its deuterated isotopologue with transient absorption and 2D IR spectroscopy. The OH bend region reveals a single dominant diagonal transition near 1625 cm-1 with intense cross peaks to lower-frequency modes, demonstrating highly mixed fingerprint transitions that contain OH bend character. The anharmonic coupling of the OH bend results in a highly elongated OH bend excited-state absorption transition that indicates a large manifold of OH bend overtone/combination bands in the OH stretch region that leads to strong bend-stretch Fermi resonance interactions. The OH and OD stretch regions consist of broad ground-state bleach signals, but there is no clear evidence of ω21 excited-state absorptions due to rapid population relaxation arising from strong intramolecular coupling to bending, fingerprint, and low-frequency H-bond modes. Orientational relaxation dynamics persist for timescales longer than the vibrational lifetimes, with polarization anisotropy components decaying within approximately 2 and 10 periods of the O-O oscillation for the OH and OD stretch, respectively. The significant isotopic dependence of the orientational dynamics is discussed in the context of intramolecular mode coupling, diffusional processes, and contributions from proton/deuteron transfer dynamics.

Hydrolysis and Condensation of Tetraethyl Orthosilicate at the Air–Aqueous Interface: Implications for Silica Nanoparticle Formation

The silica (SiO2) nanoparticles of a well-known silica precursor tetraethyl orthosilicate (TEOS) are generally synthesized via a promising solution-gelation inorganic polymerization process. The monodisperse silica nanoparticles have potential applications ranging from the formation of dental nanocomposites to antireflective coatings. In the present study, we have systematically investigated the in situ interfacial molecular structure of TEOS and its impact on the interfacial water structure during the processes of hydrolysis and condensation at the air–aqueous interface using sum-frequency generation (SFG) vibrational spectroscopy. With the presence of water, a gradual decrease in the SFG intensity in the CH-stretch region for each concentration of TEOS with time reflects the elimination of ethoxy groups, which is a signature of the hydrolysis process. Further, the impact of the hydrolysis process is revealed from the significant enhancement in the SFG signal in the OH-stretch region. The hydrolysis of TEOS is then followed by condensation in which the ≡Si─O– charged species are replaced by forming the ≡Si─O─Si≡ bridging network. The signature of the condensation process is reflected with the gradual decrease in the observed enhanced SFG signal in the OH-stretch region. The formation of monodispersed silica nanoparticles as an end product of size variation from 1.75 to 4.67 nm with the increase in TEOS concentration is confirmed with the DLS measurements. We have also probed the pH-dependent SFG studies at three different pH values (2.0, 5.8, and 9.0). The dominant pH-dependent hydrolysis process is revealed from the observed molecular structure of TEOS at the air–aqueous interface.

Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD).

We present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-molecular dynamics (MD) classical and nuclear quantum dynamical simulations, in which the energies and forces of the full system are approximated by the two-, three-, and four-body terms of the MBE, and compare the average potential and the vibrational density of states with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the four-body term converges with near-identical behavior to the one from the full simulation. The three-body makes a substantial contribution (∼20%) to the energy, whereas the four-body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (∼2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm-1) at the four-body level and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a four-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N4, N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that the MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.

Raman spectroscopy of high salinity brines and ices

Raman spectroscopy is an ideal tool to analyze the geochemistry and mineralogy of heterogenous mixtures of solids, liquid, and gases in situ, while maintaining planetary protection protocols. Here we characterize saturated CaCl2, MgCl2, MgSO4, Na2SO4, NaCl, and NaClO4 brines, as well as ultrapure water, and mixed MgSO4-NaCl, MgSO4-NaClO4, Na2SO4-NaCl, Na2SO4-NaClO4, and NaCl-NaClO4 brines from 200 K to 295 K to determine how changes in temperature affect spectral signatures of planetary analogue brines. The resulting reference dataset can be used to interpret spectra from future samples analyzed in situ on planetary bodies. Sulfate and perchlorate brines produced clear, distinct peaks associated with each polyatomic anion. While chloride brines did not produce anion peaks, subtle changes were observed in the OH-stretching region, suggesting changes to the molecular water vibration states due to complexation. Solid-liquid phase transitions were clearly observed in each of the solutions using both 785 nm (red) and 532 nm (green) excitation lasers, particularly in the OH-stretching region between 3000 and 4000 cm−1 with the 532 nm laser. Differences observed in the spectra of frozen sulfate brines suggest that cooling rates may influence the hydration state and/or crystallinity of the solid magnesium and sodium- sulfate salts. These experiments and the resulting spectral library will allow future researchers to use Raman spectroscopy to look for in situ melting, freezing, evaporation, and deliquescence as well as identify the composition of high salinity brines and their frozen products in a range of planetary environments, including permafrost and recurring slope lineae on Mars, potential ice and salt-rich regolith on asteroids such as Ceres, and ice shells and possible seeps or geysers on icy moons and other bodies.

Two-Dimensional Adiabatic Model for Calculating Progressions Resulting from Stretch–Rock Coupling in Vibrational Spectra of Anion–Water Complexes

Several anion-water dimers feature a distinct progression in the OH stretch region of their vibrational spectra. This progression arises from strong anharmonic couplings between the OH stretch and low-frequency intermolecular modes. In this work, we introduce a two-dimensional adiabatic model accounting explicitly for the water and anion rock degrees of freedom and use it to calculate the vibrational spectra of HCO2-·(H2O) and NO3-·(H2O). The spectra calculated by using this model are in excellent agreement with experiment, in terms of both peak spacings and lengths of the progressions, and represent a substantial improvement over earlier models.