Search NASA⌕ Search

SEARCH · Search NASA

Results for “isotope effect”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electrochemical Kinetics of Stable Isotopes (Final Report)

Stable isotopes are sensitive indicators to the mass transport and chemical reactions that occur during potentiostat-controlled electroplating of metals from aqueous (and some non-aqueous) metal-salt solutions. The goal of this research program is to measure the stable isotope fractionation during electrochemical processes, how the isotope-dependent rates are sensitive to chemistry, temperature, and electrochemical kinetics, and to use these measurements to develop a theoretical framework for a kinetic isotope effect that is able to predict kinetic isotope fractionation factors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Profiles of alpha 13 C and alpha D in methane from the lower stratosphere

Methane is an important greenhouse gas of biogenic and anthropogenic origin for which global budgets are being constructed from a variety of data. One approach to the global methane budget is the use of the stable isotopes C-13 and D, and the radionuclide C-14 as tracers. The authors measured the isotopic composition of methane from various sources and in tropospheric air for a number of locations. Here, the authors report on the isotopic composition of methane from the lower stratosphere. Measurements of this concentration in the stratosphere can yield estimates for the kinetic isotope effects in the methane destruction reactions. These effects have to be known for quantitative isotopic methane budgets.

Wahlen, Martin↗
NASA NTRS1999

Technical Considerations for Advanced Propulsion Systems

This presentation reviews concerns involving advanced propulsion systems. The problems involved with the use of Am-242m, is that it has a high "eta" plus an order of magnitude larger fission cross section than other fissionable materials, and that it is extremely rare. However other americium isotopes are much more common, but extremely effective isotopic separation is required. Deuterium-Tritium fusion is also not attractive for space propulsion applications. Because the pulsed systems cannot breed adequate amounts of tritium and it is difficult and expensive to bring tritium from Earth. The systems that do breed tritium have severely limited performance. However, other fusion processes should still be evaluated. Another problem with advanced propellants is that inefficiencies in converting the total energy generated into propellant energy can lead to tremendous heat rejection requirements. Therefore Many. advanced propulsion concepts benefit greatly from low-mass radiators.

Houts, Michael G.↗
DOE OSTI2022

Aldol Condensation and Esterification over Ti-Substituted *BEA Zeolite: Mechanisms and Effects of Pore Hydrophobicity

Aldol condensation and esterification reactions provide paths to upgrade ethanol and acetaldehyde to higher- value molecules useful as fuels or intermediates for the synthesis of polymers. Transition-metal-substituted BEA zeolites (M-BEA) catalyze these reactions; however, the mechanisms for these processes in M-BEA and the effects of incidental or purposefully included silanol groups are not reported. Here, we combine kinetic and spectroscopic measurements obtained during catalytic reactions of acetaldehyde (CH 3 CHO), ethanol (C 2 H 5 OH), and hydrogen (H 2 ) mixtures over a series of Ti-BEA catalysts that possess a known range of silanol group densities to examine the kinetic relevance of intervening steps and the impact of silanol groups on catalytic rates. Across all Ti-BEA, rates for aldol condensation and esterification increase with the pressure of CH 3 CHO; however, C 2 H 5 OH and H 2 O weakly inhibit the rates of these reactions. The substitution of CD 3 CDO for CH 3 CHO decreases aldol condensation rates slightly (~10%) but leads to greater esterification rates (2- to 5-fold). The kinetic isotope effects together with the measured dependence of rates on reactant pressures suggest that aldol condensation and esterification occur on unoccupied Ti sites and involve multiple kinetically relevant steps. CH 3 CHO deprotonates irreversibly, and the kinetically relevant nucleophilic attack of the enolate to CH 3 CHO* (i.e., adsorbed CH 3 CHO on Ti sites) leads to aldol products, while the nucleophilic attack of the enolate to C 2 H 5 OH* gives esters. Selectivities toward aldol condensation increase with the ratio of CH 3 CHO to C 2 H 5 OH pressure and with increases in the silanol density of the as-synthesized Ti-BEA. During catalysis, in situ infrared spectroscopy demonstrates that these silanol groups react with C 2 H 5 OH to form ethoxysilane groups (i.e., SiOC 2 H 5 ) that modify the polarity of the environment near Ti active sites. As initial silanol densities increase, steady-state turnover rates for aldol condensation and esterification increase by factors of 5 and 2, respectively. The changes in rates and selectivities among Ti-BEA catalysts likely reflect changes in excess free energies of transition states for enolization and nucleophilic attack of the enolate to adsorbed coreactants. The differences in excess stability report on the interactions among reactive intermediates at framework Ti atoms and the ethoxysilane and remaining silanol groups present. Furthermore, the in situ modification of these pore environments confers changes in the stability of reactive species in a manner that contradicts intuition when considering the initial state of the catalyst but can be reconciled after accounting for the formation of persistent alkoxy surface moieties in the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Four Isotope-Labeled Recombination Pathways of Ozone Formation

A theoretical approach is developed for the description of all possible recombination pathways in the ozone forming reaction, without neglecting any process a priori, and without decoupling the individual pathways one from another. These pathways become physically distinct when a rare isotope of oxygen is introduced, such as 18 O, which represents a sensitive probe of the ozone forming reaction. Each isotopologue of O 3 contains two types of physically distinct entrance channels and two types of physically distinct product wells, creating four recombination pathways. Calculations are done for singly and doubly substituted isotopologues of ozone, eight rate coefficients total. Two pathways for the formation of asymmetric ozone isotopomer exhibit rather different rate coefficients, indicating large isotope effect driven by ΔZPE-difference. Rate coefficient for the formation of symmetric isotopomer of ozone (third pathway) is found to be in between of those two, while the rate of insertion pathway is smaller by two orders of magnitude. These trends are in good agreement with experiments, for both singly and doubly substituted ozone. The total formation rates for asymmetric isotopomers are found to be somewhat larger than those for symmetric isotopomers, but not as much as in the experiment. Altogether, the distribution of lifetimes is found to be very similar for the metastable states in symmetric and asymmetric ozone isotopomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Evaluating the Origins of Aerobic Oxidation Catalysis with TAM-3, a MOF with Accessible Co(II) Sites and Large Pores

Metal-organic frameworks (MOFs) are attractive platforms that merge concepts of homogeneous and heterogeneous catalysis. Catalyst design and optimization are enabled by an array of synthetic methods that offer independent control over the local chemical structure of lattice-embedded metal ions (i.e., ligand identity and geometry) and the long-range materials properties (i.e., porosity). Establishing the origin of catalytic activity in MOF-promoted reactions remains a significant challenge: The relative rates of catalyst turnover and substrate diffusion dictate the extent to which interstitial sites are accessible and operational in catalysis. To minimize the contributions of surface sites in catalysis, materials with large pore dimensions are often sought, however, the impact of pore expansion on the origins of catalytic activity is similarly challenging to establish. Here, we describe TAM-3, a Co(II) based MOF with accessible metal sites supported by a facially coordinating tris-tetrazole ligand set. TAM-3 features large channel-like pores (17 × 23 Å) and promotes aerobic C−H oxidation and olefin epoxidation. Using a set of simple kinetics experiments, based on the analysis of kinetic isotope effects and olefin oxidation diastereoselectivities, we demonstrate that despite the large pores, interstitial metal ions do not significantly contribute to the observed substrate oxidation. This study highlights the importance of conducting kinetic experiments to assess the origin of apparent catalytic activity with MOFs and the challenge of harnessing reactive oxidants with microporous catalyst materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1978

Calcium isotopic anomalies in the Allende meteorite

Isotopic anomalies in Ca are reported which were found in two Ca-Al-rich inclusions of the Allende meteorite. The corrected Ca data show nonlinear isotopic effects in Ca-48 of +13.5 per mil and in Ca-42 of +1.7 per mil for one sample. The second sample shows a Ca-48 depletion of -2.9 per mil, but all other isotopes are normal. Samples with large excesses in Mg-26 show no Ca anomalies. The effects demonstrate that isotopic anomalies exist for higher-atomic-number refractory elements in solar-system materials and do not appear to be readily explainable by a simple model. Observed anomalous Mg and Ca compositions for coexisting mineral phases are uniform within each inclusion and require initial isotopic homogeneity within an inclusion but the preservation of wide variations between inclusions. Assuming formation of these inclusions by condensation from a gaseous part of the solar nebula, this implies isotopic heterogeneity on a scale of 10 to 100 km within the nebula.

Lee, T.↗
DOE OSTI2022

Kinetic Probes of the Origin of Activity in MOF-Based C–H Oxidation Catalysis

The development of homogeneous catalysis is enabled by the availability of a rich toolkit of kinetics experiments, such as the Hgdrop test, that differentiate catalytic activity at ligand-supported metal complexes from potential heterogeneous catalysts derived from decomposition of molecular species. Metal-organic frameworks (MOFs) have garnered significant attention as platforms for catalysis at site-isolated, interstitial catalyst sites. Unlike homogeneous catalysis, relatively few strategies have been advanced to evaluate the origin of catalytic activity in MOF-promoted reactions. Many of the MOFs that have been examined as potential catalysts are comprised of molecular constituents that represent viable catalysts in the absence of the extended MOF lattice, and thus interfacial sites and leached homogeneous species represent potential sources of catalyst activity. Here, we demonstrate that analysis of deuterium kinetic isotope effects (KIEs) and olefin epoxidation diastereoselectivity provides probes of the origin of catalytic activity in MOF-promoted oxidation reactions. These analyses support involvement of lattice-based Fe sites in the turnover-limiting step of C–H activation with Fe-MOF-74-based materials (i.e., the MOF functions as a bona fide catalyst) and evidence that Cu 2 -based MOF MIL-125-Cu 2 O 2 functions as a solid-state initiator for solution-phase oxidation chemistry and is not involved in the turnover limiting step (i.e., the MOF does not function as a catalyst for substrate functionalization). furthermore, we anticipate that the simple experiments described here will provide a valuable tool for clarifying the role of MOFs in C–H oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Experimental Investigation of the Thermal Decomposition Pathways and Kinetics of TATB by Isotopic Substitution

Real-time measurements of the product gases arising from the thermal decomposition of triamino-trinitro benzene (TATB), its deuterated analogue, and plastically bonded TATB (LX-17) are presented in this study. Gas-phase decomposition products are identified by IR absorption spectroscopy. The frequency shifts in rovibrational spectra due to isotopic substitution and the change in rate of formation of decomposition products due to the kinetic-isotope-effect (KIE) help elucidate the decomposition pathways. The formation of H 2 O precedes other molecules (e.g., HCN, HNCO) during decomposition. After the concentrations of HCN and HNCO molecules reach a peak, their amounts gradually decrease. The concentrations of the other decomposition products (e.g., NH 3 and CO 2 ) rapidly rise after an induction period, which is attributed to the presence of autocatalytic reactions. The trends of chemical evolution are similar for all the samples, but their kinetic behaviors are different. This indicates the rates of consistent pathways are changed during thermal decomposition. Additionally, the kinetics of deuterated TATB decomposition is slower than that of unsubstituted TATB due to the KIE (kH/kD~1.41). The rate of LX-17 decomposition is slightly lower than unsubstituted TATB (k TATB /k LX-17 ~1.15). The KIE is more pronounced during the early stage of decomposition, which is attributed to the first steps of TATB decomposition involving water formation (i.e., H vs D transfer). The KIE slows down the formation of all gases, including those lacking hydrogen (e.g., CO 2 ). These results suggest the TATB thermal decomposition mechanism might involve a series of pathways rather than a set of independent and parallel reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Mechanistic Studies of an Iron-Catalyzed Intermolecular C–H Amination Reaction under Catalytic Conditions and Having a Large KIE

The conversion of C–H bonds into amines by nitrene insertion is an attractive transformation since it is both atom- and step-economical, and provides a direct route to functionalizing hydrocarbons. Using an iron catalyst [{( tBu pyrr) 2 pyr}Fe(OEt 2 )] (1-OEt 2 ) (( tBu pyrr) 2 pyr 2– = 3,5- t Bu 2 -bis(pyrrolyl)pyridine), we recently demonstrated the catalytic conversion of weak C–H bonds into secondary amines using aryl azides as the nitrene source [Zars, E.; Angew. Chem., Int. Ed. 2023, 62, e202311749]. Here, we describe detailed mechanistic studies of this intermolecular C–H amination reaction under catalytic conditions. We find by Variable Time Normalization Analysis (VTNA) that the conversion of xanthene (2-H 2 ) and 2,4,6-trimethyl-phenyl azide ( Me 3) catalyzed by 1-OEt 2 is an overall 3/2 order process, being 1 st order in 2-H 2 and half order in Me 3. A kinetic isotope effect study (KIE) using 2-d 2 results in a significant decrease in the rate (KIE = 61(15)), which clearly implicates the C–H insertion step as rate-determining. Furthermore, treatment of 1-OEt 2 with one equivalent of N 3 -2,6- i Pr 2 –C 6 H 3 yields the mixed-valence C–N coupled product [( tBu pyrr) 2 pyrFe-N═C(2,6 i Pr 2 –Ph)═N-(2,6 i Pr 2 –Ph))Fe tBu pyrrpyr(2-H-pyrr)] (5 iPr ). Quantum chemical calculations confirm the electronic structure of the mixed-valence dimer in 5 iPr and rationalize the Hammett correlation by a delicate balance in the dinuclearization of the catalytically active monomers. Calculations further indicate significant tunneling for the pivotal H atom abstraction by the iron-imidyl complex. Combining all these results allows us to propose a mechanism consisting of imido formation in equilibrium with a radical-coupled diiron system, followed by stepwise C–H insertion via a linear H atom abstraction transition state and subsequent radical rebound.

azides↗
DOE OSTI2026

Progress in Extending the High-Temperature Plasma Regime in the LHD

In this paper, we show the expansion of the high-temperature regime and related physical phenomena in the Large Helical Device (LHD). At the LHD, we have developed a high-temperature operational regime for both ion (T i ) and electron (T e ) temperatures through understanding the physical characteristics of high-temperature plasma, the enhancement of plasma heating, and the development of technical methods such as wall recycling control and impurity control. In particular, the formation of the ion internal transport barrier (i-ITB) and electron internal transport barrier (e-ITB) in the plasma core region is the key to improving plasma performance. In addition, the deuterium experiment that began in 2017 showed that, due to the isotope effect, the thermal transport of both ions and electrons is suppressed in deuterium plasma compared to light hydrogen plasma, and that a high central temperature can be achieved efficiently. Furthermore, the combination of neutral beam injection (NBI) and electron cyclotron heating (ECH) has extended the operating conditions that simultaneously maintain high T i and T e , and comprehensive research on plasma confinement characteristics with an eye to future fusion reactor conditions has progressed.

Large Helical Device↗
DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 2: Experimental determination of hydrogen isotopic equilibrium for ethane-H2 from 30 to 200 °C and propane-H2 from 75 to 200 °C

The stable isotopic compositions of light n-alkanes, including methane, ethane, and propane, are often used to identify the sources and thermal maturity of natural gas samples. Though stable isotopic compositions of these molecules are commonly assumed to be controlled by kinetic isotope effects, recent studies have proposed both carbon and hydrogen isotopic equilibrium may also occur in some samples. Assessing whether samples are in isotopic equilibrium requires knowledge of light alkane equilibrium fractionation factors over geologically relevant temperatures for formation and storage (up to ∼300 °C). In this study, we report experimental results of hydrogen isotopic equilibrium between ethane and H2 from 30 to 200 °C and propane and H2 from 75 to 200 °C. We compare these results with high-level theoretical calculations and provide a preferred polynomial fit to describe equilibrium fractionation factors. Comparison of these fractionation factors with a compilation of ∼500 compiled environmental gas samples supports the proposal that many (∼50%) of these natural gas samples exhibit hydrogen isotopic compositions consistent with having formed in or attained methane-ethane-propane hydrogen isotopic equilibrium over geologically relevant temperatures for formation and storage (50–300 °C).

Turner, Andrew C↗
DOE OSTI2021

Fossil records of early solar irradiation and cosmolocation of the CAI factory: A reappraisal

Calcium-aluminum–rich inclusions (CAIs) in meteorites carry crucial information about the environmental conditions of the nascent Solar System prior to planet formation. Based on models of 50V– 10 Be co-production by in-situ irradiation, CAIs are considered to have formed within ~0.1 AU from the proto-Sun. Here, we present vanadium (V) and strontium (Sr) isotopic co-variations in fine- and coarse-grained CAIs and demonstrate that kinetic isotope effects during partial condensation and evaporation best explain V isotope anomalies previously attributed to solar particle irradiation. We also report initial excesses of 10 Be and argue that CV CAIs possess essentially a homogeneous level of 10 Be, inherited during their formation. Based on numerical modeling of 50V–10Be co-production by irradiation, we show that CAI formation during protoplanetary disk build-up likely occurred at greater heliocentric distances than previously considered, up to planet-forming regions (~1AU), where solar particle fluxes were sufficiently low to avoid substantial in-situ irradiation of CAIs.

58 GEOSCIENCES↗
NASA NTRS2013

Carbonate Mineral Formation on Mars: Clues from Stable Isotope Variation Seen in Cryogenic Laboratory Studies of Carbonate Salts

The geologic history of water on the planet Mars is intimately connected to the formation of carbonate minerals through atmospheric CO2 and its control of the climate history of Mars. Carbonate mineral formation under modern martian atmospheric conditions could be a critical factor in controlling the martian climate in a means similar to the rock weathering cycle on Earth. The combination of evidence for liquid water on the martian surface and cold surface conditions suggest fluid freezing could be very common on the surface of Mars. Cryogenic calcite forms readily when a rise in pH occurs as a result of carbon dioxide degassing quickly from freezing Ca-bicarbonate-rich water solutions. This is a process that has been observed in some terrestrial settings such as arctic permafrost cave deposits, lakebeds of the Dry Valleys of Antarctica, and in aufeis (river icings) from rivers of N.E. Alaska. We report here the results of a series of laboratory experiments that were conducted to simulate potential cryogenic carbonate formation on the planet Mars. These results indicate that carbonates grown under martian conditions (controlled atmospheric pressure and temperature) show enrichments from starting bicarbonate fluids in both carbon and oxygen isotopes beyond equilibrium values with average delta13C(DIC-CARB) values of ~20.5%0 which exceed the expected equilibrium fractionation factor of [10(sup 3) ln alpha = ~13%0] at 0 degC. Oxygen isotopes showed a smaller enrichment with delta18O(H2O-CARB) values of ~35.5%0, slightly exceeding the equilibrium fractionation factor of [10(sup 3) ln alpha = ~34%0 ] at 0degC. Large kinetic carbon isotope effects during carbonate precipitation could substantially affect the carbon isotope evolution of CO2 on Mars allowing for more efficient removal of 13C from the Noachian atmosphere enriched by atmospheric loss. This mechanism would be consistent with the observations of large carbon isotope variations in martian materials despite the relative paucity of carbonate minerals in the martian crust.

Socki, Richard↗
DOE OSTI2021

Effect of hydrogen isotope on plasma impedance and thermal pinching induced by single aerosol droplets injected into an inductively coupled plasma

The mass of the hydrogen isotope in micrometer-sized aqueous droplets introduced into an ICP has been found to affect the degree of resulting plasma perturbation. Studied plasma perturbations include local plasma cooling, plasma shrinkage due to thermal pinching, and changes in plasma impedance. The local cooling effect caused by a vaporizing aerosol droplet was found to be similar between H 2 O and D 2 O. This finding is likely a result of the similar heat capacities and vaporization rates of H 2 O and D 2 O droplets and of the comparable thermal conductivities and dissociation kinetics of H 2 O and D 2 O vapor. In contrast, a clear isotope effect was observed for plasma shrinkage and impedance change, with H 2 O causing a considerably stronger effect. Plasma shrinkage, gauged by the reduction in Ar emission at the edge of the plasma, was detected at all observation heights within the load-coil region, with D 2 O droplets causing as much as 30% less of an effect than H 2 O. Likewise, D 2 O droplets induce significantly less change in plasma impedance than H 2 O droplets. Depending on experimental conditions and the probing method, changes in plasma impedance caused by monodisperse D 2 O droplets were only 52%–80% that of H 2 O. The lower influence on plasma shrinkage and impedance change for D 2 O droplets is believed to be related to the much lower thermal conductivity of atomic deuterium, around 70% that of atomic hydrogen (protium). The time required for the plasma impedance to return to its original steady-state value was comparatively long (~10 ms) and similar for H 2 O and D 2 O droplets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Isotopic inferences of ancient biochemistries - Carbon, sulfur, hydrogen, and nitrogen

In processes of biological incorporation and subsequent biochemical processing sizable isotope effects occur as a result of both thermodynamic and kinetic fractionations which take place during metabolic and biosynthetic reactions. In this chapter a review is provided of earlier work and recent studies on isotope fractionations in the biogeochemical cycles of carbon, sulfur, hydrogen, and nitrogen. Attention is given to the biochemistry of carbon isotope fractionation, carbon isotope fractionation in extant plants and microorganisms, isotope fractionation in the terrestrial carbon cycle, the effects of diagenesis and metamorphism on the isotopic composition of sedimentary carbon, the isotopic composition of sedimentary carbon through time, implications of the sedimentary carbon isotope record, the biochemistry of sulfur isotope fractionation, pathways of the biogeochemical cycle of nitrogen, and the D/H ratio in naturally occurring materials.

Schidlowski, M.↗
DOE OSTI2021

Unimolecular and bimolecular formic acid decomposition routes on dispersed Cu nanoparticles

The elementary steps and site requirements in formic acid (HCOOH) dehydrogenation on Cu surfaces remain of keen interest because formate species act as intermediates or spectators in methanol synthesis and water–gas shift reactions. Steady-state and transient kinetic data, isotopic effects, infrared spectra during catalytic and stoichiometric reactions, and theoretical treatments based on density functional theory (DFT) provide evidence for bimolecular reactions, in which saturated bidentate formate (*HCOO*) adlayers, present at 0.25 ML (0.25 *HCOO* per surface Cu atom), react with undissociated species (HCOOH$^\square$) bound at interstices within formate adlayers ($^\square$) to form H-bonded bimolecular HCOOH$^\square$-*HCOO* adducts. The co-existence of vicinal HCOOH$^\square$ and *HCOO* moieties is evident from antisymmetric infrared bands for *HCOO* that become stronger as a result of their H-bonding that perturbs the induced dipole moment of *HCOO* upon vibration, consistent with DFT-derived vibrational frequencies and intensities for such perturbed species. The *HCOO* moiety in this complex undergoes C-H activation via a transition state that is preferentially stabilized through H-bonding with the vicinal HCOOH$^\square$ relative to its *HCOO* precursor. DFT-derived HCOOH dehydrogenation activation barriers and those determined from the evolution of CO 2 from pre-adsorbed *HCOO* species are about 10 kJ mol -1 smaller in the presence of gaseous HCOOH reactants (because of HCOOH$^\square$-*HCOO* interactions) than those for the unimolecular decomposition of bound *HCOO* species. Such bimolecular routes are consistent with measured effects of HCOOH, H 2 , and CO pressures and of H/D isotopic substitution on dehydrogenation turnover rates and represent the predominant channel for the formation of CO 2 and H 2 during catalytic HCOOH dehydrogenation on Cu nanoparticles. A saturated *HCOO* adlayer that retains binding interstices and the presence of HCOOH(g) enable a sequence of elementary steps unavailable for *HCOO* species, thus circumventing unassisted unimolecular routes that exhibit higher activation barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Excited State Intramolecular Proton Transfer with Nuclear-Electronic Orbital Ehrenfest Dynamics

The recent development of the Ehrenfest dynamics approach in the nuclear-electronic orbital (NEO) framework provides a promising way to simulate coupled nuclear-electronic dynamics. Our previous study showed that the NEO-Ehrenfest approach with a semi- classical traveling proton basis method yields accurate predictions of molecular vibrational frequencies. In this work, we provide a more thorough analysis of the semiclassical traveling proton basis method to elucidate its validity and convergence behavior. We also conduct NEO-Ehrenfest dynamics simulations to study an excited state intramolecular proton transfer process. These simulations reveal that nuclear quantum effects influence the predictions of proton transfer reaction rates and kinetic isotope effects due to the intrinsic delocalized nature of the quantum nuclear wave function. This work illustrates the importance of nuclear quantum effects in coupled nuclear-electronic dynamical processes and shows that the NEO-Ehrenfest approach can be a powerful tool to provide insights and predictions for these processes.

Zhao, Luning↗