Developing efficient path integral (PI) methods for atomistic simulations of vibrational spectra in heterogeneous condensed phases and interfaces has long been a challenging task. Here, we present the h-CMD method, short for hybrid centroid molecular dynamics, which combines the recently introduced fast quasi-CMD (f-QCMD) method with fast CMD (f-CMD). In this scheme, molecules that are believed to suffer more seriously from the curvature problem of CMD, e.g., water, are treated with f-QCMD, while the rest, e.g., solid surfaces, are treated with f-CMD. To test the accuracy of the newly introduced scheme, the infrared spectra of the interfacial D2O confined in the archetypal ZIF-90 framework are simulated using h-CMD compared to a variety of other PI methods, including thermostatted ring-polymer molecular dynamics (T-RPMD) and partially adiabatic CMD as well as f-CMD and experiment as reference. Comparisons are also made with classical MD, where nuclear quantum effects are neglected entirely. Our detailed comparisons at different temperatures of 250–600 K show that h-CMD produces O–D stretches that are in close agreement with the experiment, correcting the known curvature problem and redshifting of the stretch peaks of CMD. h-CMD also corrects the known issues associated with too artificially dampened and broadened spectra of T-RPMD, which leads to missing the characteristic doublet feature of the interfacial confined water, rendering it unsuitable for these systems. The new h-CMD method broadens the applicability of f-QCMD to heterogeneous condensed phases and interfaces, where defining curvilinear coordinates for the entire system is not feasible.
more »
« less
A local mode study of ring puckering effects in the infrared spectra of cyclopentane
We report and interpret recently recorded high-resolution infrared spectra for the fundamentals of the CH 2 scissors and CH stretches of gas phase cyclopentane at −26.1 and −50 °C, respectively. We extend previous theoretical studies of this molecule, which is known to undergo barrierless pseudorotation due to ring puckering, by constructing local mode Hamiltonians of the stretching and scissor vibrations for which the frequencies, couplings, and linear dipoles are calculated as functions of the pseudorotation angle using B3LYP/6-311++(d,p) and MP2/cc-pVTZ levels of theory. Symmetrization ( D 5 h ) of the vibrational basis sets leads to simple vibration/pseudorotation Hamiltonians whose solutions lead to good agreement with the experiment at medium resolution, but which miss interesting line fractionation when compared to the high-resolution spectra. In contrast to the scissor motion, pseudorotation leads to significant state mixing of the CH stretches, which themselves are Fermi coupled to the scissor overtones.
more »
« less
- Award ID(s):
- 1900095
- PAR ID:
- 10350692
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 21
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 214305
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The photodissociation dynamics of the dimethyl-substituted acetone oxide Criegee intermediate [(CH 3 ) 2 COO] is characterized following electronic excitation on the π*←π transition, which leads to O ( 1 D) + acetone [(CH 3 ) 2 CO, S0] products. The UV action spectrum of (CH 3 ) 2 COO recorded with O ( 1 D) detection under jet-cooled conditions is broad, unstructured, and essentially unchanged from the corresponding electronic absorption spectrum obtained using a UV-induced depletion method. This indicates that UV excitation of (CH 3 ) 2 COO leads predominantly to the O ( 1 D) product channel. A higher energy O ( 3 P) + (CH3)2CO (T1) product channel is not observed, although it is energetically accessible. This is attributed to the relatively weak absorption cross section at UV excitation energies above the threshold. In addition, complementary MS-CASPT2 trajectory surface-hopping (TSH) simulations indicate minimal population leading to the O ( 3 P) channel and non-unity overall probability for dissociation (within 100 fs). Velocity map imaging of the O ( 1 D) products is utilized to reveal the total kinetic energy release (TKER) distribution upon photodissociation of (CH 3 ) 2 COO at various UV excitation energies. Simulation of the TKER distributions is performed using a hybrid model that combines an impulsive model with a statistical component, the latter reflecting the longer-lived (> 100 fs) trajectories identified in the TSH calculations. The impulsive model accounts for vibrational activation of (CH 3 ) 2 CO arising from geometrical changes between the Criegee intermediate and the carbonyl product, indicating the importance of CO stretch, CCO bend, and CC stretch along with activation of hindered rotation and rock of the methyl groups in the (CH 3 ) 2 CO product. Detailed comparison is also made with the TKER distribution arising from photodissociation dynamics of CH 2 OO upon UV excitation.more » « less
-
Sum frequency generation (SFG) * Equal contributors. spectroscopy was used to deduce the orientation of the terminal methyl (CH3) group of self-assembled monolayers (SAMs) at the air–solid and air–liquid interfaces at surface concentrations as low as 1% protonated molecules in the presence of 99% deuterated molecules. The SFG spectra of octadecanethiol (ODT) and deuterated octadecanethiol (d37ODT) SAMs on gold were used for analysis at the air–solid interface. However, the eicosanoic acid (EA) and deuterated EA (d39EA) SAMs on the water were analyzed at the air–liquid interface. The tilt angle of the terminal CH3group was estimated to be ∼39 ° for a SAM of 1% ODT : 99% d37ODT, whereas the tilt angle of the terminal CH3group of the 1% EA : 99% d39EA monolayer was estimated to be ∼32 °. The reliability of the orientational analysis at low concentrations was validated by testing the sensitivity of the SFG spectroscopy. A signal-to-noise (S/N) ratio of ∼60 and ∼45 was obtained for the CH3symmetric stretch (SS) of 1% ODT : 99% d37ODT and 1% EA : 99% d39EA, respectively. The estimated increase in S/N ratio values, as a measure of the sensitivity of the SFG spectroscopy, verified the capacity to acquire the SFG spectra at low concentrations of interfacial molecules under ambient conditions. Overall, the orientational analysis of CH3SS vibrational mode was feasible at low concentrations of protonated molecules due to increased S/N ratio. In support, the improved S/N ratio on varying incident power density of the visible beam was also experimentally demonstrated.more » « less
-
Context. At present, there are strong indications that white dwarf (WD) stars with masses well below the Chandrasekhar limit ( M Ch ≈ 1.4 M ⊙ ) contribute a significant fraction of SN Ia progenitors. The relative fraction of stable iron-group elements synthesized in the explosion has been suggested as a possible discriminant between M Ch and sub- M Ch events. In particular, it is thought that the higher-density ejecta of M Ch WDs, which favours the synthesis of stable isotopes of nickel, results in prominent [Ni II ] lines in late-time spectra (≳150 d past explosion). Aims. We study the explosive nucleosynthesis of stable nickel in SNe Ia resulting from M Ch and sub- M Ch progenitors. We explore the potential for lines of [Ni II ] in the optical an near-infrared (at 7378 Å and 1.94 μm) in late-time spectra to serve as a diagnostic of the exploding WD mass. Methods. We reviewed stable Ni yields across a large variety of published SN Ia models. Using 1D M Ch delayed-detonation and sub- M Ch detonation models, we studied the synthesis of stable Ni isotopes (in particular, 58 Ni) and investigated the formation of [Ni II ] lines using non-local thermodynamic equilibrium radiative-transfer simulations with the CMFGEN code. Results. We confirm that stable Ni production is generally more efficient in M Ch explosions at solar metallicity (typically 0.02–0.08 M ⊙ for the 58 Ni isotope), but we note that the 58 Ni yield in sub- M Ch events systematically exceeds 0.01 M ⊙ for WDs that are more massive than one solar mass. We find that the radiative proton-capture reaction 57 Co( p , γ ) 58 Ni is the dominant production mode for 58 Ni in both M Ch and sub- M Ch models, while the α -capture reaction on 54 Fe has a negligible impact on the final 58 Ni yield. More importantly, we demonstrate that the lack of [Ni II ] lines in late-time spectra of sub- M Ch events is not always due to an under-abundance of stable Ni; rather, it results from the higher ionization of Ni in the inner ejecta. Conversely, the strong [Ni II ] lines predicted in our 1D M Ch models are completely suppressed when 56 Ni is sufficiently mixed with the innermost layers, which are rich in stable iron-group elements. Conclusions. [Ni II ] lines in late-time SN Ia spectra have a complex dependency on the abundance of stable Ni, which limits their use in distinguishing among M Ch and sub- M Ch progenitors. However, we argue that a low-luminosity SN Ia displaying strong [Ni II ] lines would most likely result from a Chandrasekhar-mass progenitor.more » « less
-
Electrocatalytic Methane Functionalization with d 0 Early Transition Metals Under Ambient ConditionsAbstract The undesirable loss of methane (CH4) at remote locations welcomes approaches that ambiently functionalize CH4on‐site without intense infrastructure investment. Recently, we found that electrochemical oxidation of vanadium(V)‐oxo with bisulfate ligand leads to CH4activation at ambient conditions. The key question is whether such an observation is a one‐off coincidence or a general strategy for electrocatalyst design. Here, a general scheme of electrocatalytic CH4activation with d0early transition metals is established. The pre‐catalysts’ molecular structure, electrocatalytic kinetics, and mechanism were detailed for titanium (IV), vanadium (V), and chromium (VI) species as model systems. After a turnover‐limiting one‐electron electrochemical oxidation, the yielded ligand‐centered cation radicals activate CH4with low activation energy and high selectivity. The reactivities are universal among early transition metals from Period 4 to 6, and the reactivities trend for different early transition metals correlate with their d orbital energies across periodic table. Our results offer new chemical insights towards developing advanced ambient electrocatalysts of natural gas.more » « less
An official website of the United States government

