Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “species dynamics”

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 199 records · Page 11

Operando characterization and regulation of metal dissolution and redeposition dynamics near battery electrode surface

Mn dissolution has been a longstanding, ubiquitous issue that negatively impacts the performance of Mn-based battery materials. Mn dissolution involves complex chemical and structural transformations at the electrode-electrolyte interface. The continuously evolving electrodeelectrolyte interface has posed great challenges for characterizing the dynamic interfacial process and quantitatively establishing the correlation with battery performance. In this study, we visualize and quantify the temporally and spatially resolved Mn dissolution/redeposition (D/R) dynamics of electrochemically operating Mn-containing cathodes. The particle-level and electrode-level analyses reveal that the D/R dynamics is associated with distinct interfacial degradation mechanisms at different states of charge. Our results statistically differentiate the contributions of surface reconstruction and Jahn-Teller distortion to the Mn dissolution at different operating voltages. Introducing sulfonated polymers (Nafion) into composite electrodes can modulate the D/R dynamics through trapping the dissolved Mn species and rapidly establishing the local Mn D/R equilibrium. This work represents an inaugural effort to pinpoint the chemical and structural transformations responsible for Mn dissolution via an operando synchrotron study and develops an effective method to regulate Mn interfacial dynamics for improving battery performance.

25 ENERGY STORAGE↗

Dynamic behavior of molecular Pd-acetate trimers and dimers in heterogeneous vinyl acetate synthesis

Vinyl acetate monomer (VAM) is a crucial intermediate in the production of various polymers. While molecular Pd-acetate trimers and dimers, such as Pd 3 (OAc) 6 and K 2 Pd 2 (OAc) 6 , are known to form on potassium acetate (KOAc)-promoted PdAu catalysts during heterogeneous VAM synthesis, their mechanistic role remains unclear. Here, we study the dynamics of different Pd-acetate species by utilizing in situ and operando crystallographic and spectroscopic characterizations combined with computational modeling on monometallic Pd model catalysts. The promoter-free catalyst expectedly shows low catalytic activity and VAM selectivity, corresponding to the complete reduction of Pdn(OAc) 2n species to form Pd 0 and PdC x nanoparticles. Conversely, noticeable quantities of K n Pd 2 (OAc) n+4 species remain on the KOAc-promoted catalyst, leading to smaller nanoparticle formation with 10 times the activity and double the selectivity for VAM. This study reveals that molecular Pd-acetate trimers and dimers are significant indicators of catalytic performance and highlights their structurally dynamic nature in heterogeneous vinyl acetate chemistry.

Jacobs, Hunter P. [Rice Univ., Houston, TX (United↗

The transition from resistance to acceptance: Managing a marine invasive species in a changing world

Abstract Marine invasive species can transform coastal ecosystems, yet mitigating their effects can be difficult, and even impractical. Often, marine invasive species are managed at poorly matched spatial scales, and at the same time, rates of spread and establishment are increasing under climate change and can outpace resources available for population suppression. These circumstances challenge traditional conservation goals of maintaining a historic environmental state, especially for a species like the European green crab ( Carcinus maenas ), a formidable invader with few examples of successful long‐term removal programs. A management paradigm where decision alternatives include resisting or accepting a new ecological trajectory may be needed. We apply mathematical concepts from decision theory to develop a quantitative framework for navigating management decisions in this new resist‐accept paradigm. We develop a model of European green crab growth, removal and colonization, and we find optimal levels of removal effort that minimize both ecological change and removal cost. We establish a benchmark of colonization pressure at which green crab density becomes decoupled from a decision maker's actions, such that population control can no longer shape the invasion trajectory. For informing the decision boundary between resistance and acceptance, our results highlight that a decision maker's understanding of how removal cost scales with removal effort is more important than understanding the density‐impact relationship. We show that assuming stationary system dynamics can result in sub‐optimal levels of species removal effort, highlighting the importance of developing anticipatory management strategies by accounting for non‐stationary dynamics. Policy implications . For marine invasive species that can disperse across long distances and recolonize rapidly after removal, the focus of conservation policy should shift away from understanding how to resist change to understanding when to stop resisting change. Navigating this decision problem involves trade‐offs among competing objectives, highlighting the need for structured approaches to elicit objective weights that reflect the values of the decision maker. For natural resource managers facing possible ecosystem transformation, this decision framework can enable proactive and strategic decisions made under uncertainty in a changing world.

Keller, Abigail G. [Department of Environment Scie↗

On Kinetics Modeling of Vibrational Energy Transfer

Two models of vibrational energy exchange are compared at equilibrium to the elementary vibrational exchange reaction for a binary mixture. The first model, non-linear in the species vibrational energies, was derived by Schwartz, Slawsky, and Herzfeld (SSH) by considering the detailed kinetics of vibrational energy levels. This model recovers the result demanded at equilibrium by the elementary reaction. The second model is more recent, and is gaining use in certain areas of computational fluid dynamics. This model, linear in the species vibrational energies, is shown not to recover the required equilibrium result. Further, this more recent model is inconsistent with its suggested rate constants in that those rate constants were inferred from measurements by using the SSH model to reduce the data. The non-linear versus linear nature of these two models can lead to significant differences in vibrational energy coupling. Use of the contemporary model may lead to significant misconceptions, especially when integrated in computer codes considering multiple energy coupling mechanisms.

Gilmore, John O.↗

Electron–Ion Covariance Reveals Geometry-Resolved Inner-Shell Spectra in CS 2 Photodissociation

The chemical shifts of inner-shell atomic orbitals are highly sensitive to chemical bonding and molecular structure. In evolving systems, however, the spectra of distinct molecular species and geometries overlap, obscuring the underlying chemical dynamics. Here we demonstrate the use of electron–ion covariance analysis to combine the structural sensitivity of Coulomb explosion imaging with inner-shell spectroscopy, yielding geometry-specific spectra of transient and product species. We apply this approach to the excited state dynamics of CS 2 probed by ionization above the S 2p edge. Electron–ion covariance with time- and momentum-selected S + and S 2+ ions isolates distinct S 2p photoelectron spectra for ground-state CS 2 , bent photoexcited CS 2 , the CS photoproduct, and bare atomic sulfur ─ species whose spectra overlap strongly in the channel-averaged measurement. Clear chemical shifts are observed in the covariance photoelectron spectrum for each of these species, all of which are consistent with high-level calculations. Here, by extracting the atomic S contribution to the photoelectron spectrum in a finely time-resolved manner we can disentangle this contribution to the overall time-resolved photoelectron spectrum as the photodissociation proceeds. These results demonstrate the promise of electron–ion covariance as a general approach to geometry-resolved inner-shell spectroscopy, opening a route to tracking structural evolution through chemical shifts in complex photoexcited molecules.

Ionization↗

Dynamically Stable Active Sites from Surface Evolution of Perovskite Materials during the Oxygen Evolution Reaction

Perovskite oxides are an important class of oxygen evolution reaction (OER) catalysts in alkaline media, despite the elusive nature of their active sites. In this work, we demonstrate that the origin of the OER activity in a La 1-x SrxCoO 3 model perovskite arises from a thin surface layer of Co hydr(oxy)oxide (CoO x H y ) that interacts with trace-level Fe species present in the electrolyte, creating dynamically stable active sites. Generation of the hydr(oxy)oxide layer is a consequence of a surface evolution process driven by the A-site dissolution and O-vacancy creation. In turn, this imparts a 10-fold improvement in stability against Co dissolution and a 3-fold increase in the activity-stability factor for CoO x H y / LSCO when compared to nanoscale Co-hydr(oxy)oxides clusters. Our results suggest new design rules for active and stable perovskite oxide-based OER materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-resolution infrared spectroscopy of jet cooled cyclobutyl in the α-CH stretch region: large-amplitude puckering dynamics in a 4-membered ring radical

We report gas-phase cyclobutyl radical (c-C 4 H 7 ) is generated at a rotational temperature of T rot = 26(1) K in a slit-jet discharge mixture of 70% Ne/30% He and 0.5–0.6% cyclobromobutane (c-C 4 H 7 Br). A fully rovibrationally resolved absorption spectrum of the α-CH stretch fundamental band between 3062.9 cm −1 to 3075.7 cm −1 is obtained and analyzed, yielding first precision structural and dynamical information for this novel radical species. The α-CH stretch band origin is determined to be 3068.7887(4) cm −1 , which implies only a modest (≈0.8 cm −1 ) blue shift from rotationally unresolved infrared spectroscopic studies of cyclobutyl radicals in liquid He droplets [ A. R. Brown, P. R. Franke and G. E. Douberly, J. Phys. Chem. A , 2017, 121 , 7576–7587]. Of particular dynamical interest, a one-dimensional potential energy surface with respect to the ring puckering coordinate is computed at CCSD(T)/ANO2 level of theory and reveals a double minimum C s puckered geometry, separated by an exceedingly shallow planar C 2v transition state barrier (E barr ≈ 1 cm −1 ). Numerical solutions on this double minimum potential yield a zero-point energy for the ground state (E zero-point ≈ 27 cm −1 ) greatly in excess of the interconversion barrier. This is indicative of highly delocalized, large amplitude motion of the four-membered ring structure, for which proper vibrationally averaging of the moment of inertia tensor reproduces the experimentally determined inertial defect remarkably well. Finally, intensity alternation in the experimental spectrum due to nuclear spin statistics upon exchange of three indistinguishable H atom pairs (IH = ½) matches K a + K c = even : odd = 36 : 28 predictions, implying that the unpaired electron in the radical center lies in an out-of-plane p π orbital. Thus, high-resolution infrared spectroscopy provides first experimental confirmation of a shallow double minimum ring puckering potential with a highly delocalized ground state wave function peaked at a planar C 2v transition state geometry consistent with a cyclobutyl π radical.

Chemistry↗

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗

Ultrafast Dynamics and Rearrangement of the EUV Photoacid Generator Phenyl Triflate

Extreme UV photolithography is a high-energy process crucial for modern computer chip fabrication, where photoacid generators (PAGs) enhance radiation exposure by producing reactive chemical species. This study investigates the ultrafast dissociative dynamics of phenyl triflate, a widely used PAG, under high-energy ionizing conditions. The dissociative ionization of the phenyl triflate cation, which includes a molecular rearrangement that releases SO 2 and forms the phenyl trifluoromethyl ether cation, is measured over the first 5 ps. These dynamics reveal vibrational coherence corresponding to a torsional mode of the cation, involving the twisting of the phenyl group into the O–S–C plane. Electronic structure calculations for the radical cation are in good agreement with the experimentally observed vibrational coherence. These findings provide valuable insights into the behavior of phenyl triflate under ionizing conditions similar to its industrial usage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potential vorticity structure of Titan's polar vortices from Cassini CIRS observations

The Cassini mission has provided the best opportunity to date to extensively study the seasonal variation in Titan’s atmosphere, with observations spanning almost half a Titan year (𝐿𝑠 = 293 − 93°). An important feature in the Titan middle-atmosphere is the formation of a polar vortex. Observations have shown that an initially well-developed northern vortex enriched with trace gas species gradually breaks down after spring equinox as a new vortex emerges in southern winter. Here we use Cassini CIRS observations to derive the temperature and composition of the middle-atmosphere. We use the gradient wind equation to first estimate the mean zonal winds, and then the Potential Vorticity (PV) throughout Titan’s atmosphere over the timespan of the Cassini mission. PV is a useful diagnostic quantity for studying the dynamics of polar vortices because it is materially conserved for adiabatic and frictionless flows, and can be inverted to find all other dynamical fields. Our results show the formation of a strong zonal jet in the winter hemisphere, with wind velocities reaching 220 m/s, which is consistent with previous studies. An annular PV structure is also observed over the winter poles, whereby a ring of PV encircles a local minima over the pole. Such distributions are often found to be unstable without a restoring force, yet they are seen here in numerous observations in both the northern and southern hemispheres. A comparison with the annular Martian vortices shows that latent heat release from condensation or subsidence-induced adiabatic heating may explain the origin and stability of the annulus. Finally, we investigate the evolution of the size of the vortices and the role of strong PV gradients as a dynamical mixing barrier for trace gas species across the vortex edge. We find that longer lived gases are less confined to the vortex than those with shorter photochemical lifetimes.

Conor Andrew Nixon↗

Extreme elevational migration spurred cryptic speciation in giant hummingbirds

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how “outlier” species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world’s largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

Science & Technology - Other Topics↗

DEVELOP Poster: Medicine Bow Disasters

The Medicine Bow National Forest (MBNF) consists of approximately 1,383,790 acres of forested land, grassland, and sagebrush steppe in southeastern Wyoming. Cheatgrass (Bromus tectorum), an invasive plant species in the Western US, occurs in the grasslands throughout MBNF. Cheatgrass is known to rapidly colonize disturbed sites and dramatically alter historic fire regimes and nutrient/water dynamics as well as outcompete native plant species that are important forage for mule deer (Odocoileus hemionus) and elk (Cervus canadensis). In 2012, the Squirrel Creek Fire burned approximately 10,587 acres of land within MBNF, exacerbating the spread of cheatgrass. In 2015, the Wyoming Ecological Forecasting DEVELOP team identified areas of high cheatgrass abundance within the fire boundary in order to guide US Forest Service (USFS) herbicide spraying efforts to reduce cheatgrass in 2016. This research used Landsat 8 Operational Land Imager (OLI) and Sentinel-2 MultiSpectral Instrument (MSI) data to create a 2019 probabilistic cheatgrass occurrence map. This map allowed an analysis of the effectiveness of aerial spraying to inform future land management techniques for the USFS. Based on the results of the Generalized Linear Model, we found that treated areas decreased in cheatgrass cover by 36% while untreated areas increased in cheatgrass cover by 6%, suggesting that herbicide treatment has been effective.

Disasters↗

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗

Landscape fragmentation overturns classical metapopulation thinking

Habitat loss and isolation caused by landscape fragmentation represent a growing threat to global biodiversity. Existing theory suggests that the process will lead to a decline in metapopulation viability. However, since most metapopulation models are restricted to simple networks of discrete habitat patches, the effects of real landscape fragmentation, particularly in stochastic environments, are not well understood. To close this major gap in ecological theory, we developed a spatially explicit, individual-based model applicable to realistic landscape structures, bridging metapopulation ecology and landscape ecology. This model reproduced classical metapopulation dynamics under conventional model assumptions, but on fragmented landscapes, it uncovered general dynamics that are in stark contradiction to the prevailing views in the ecological and conservation literature. Notably, fragmentation can give rise to a series of dualities: a) positive and negative responses to environmental noise, b) relative slowdown and acceleration in density decline, and c) synchronization and desynchronization of local population dynamics. Furthermore, counter to common intuition, species that interact locally (“residents”) were often more resilient to fragmentation than long-ranging “migrants.” This set of findings signals a need to fundamentally reconsider our approach to ecosystem management in a noisy and fragmented world.

54 ENVIRONMENTAL SCIENCES↗

Migrating Populus with climate change: Phenology, coppice management, cold spell susceptibility, leaf dynamics, and biomass production

Understanding the phenology and productivity of Populus species is crucial for effective management and conservation strategies amid climate change. We investigated leaf budbreak timing, susceptibility to cold damage, leaf dynamics, and biomass production of 168 Populus genotypes with diverse provenances in the southeastern United States. Our study revealed significant variation in budbreak timing across different taxa and years, with genotypes inheriting traits adapted to their parents’ local climates. Temperature emerged as a key factor triggering budbreak, while leaf development depended on other environmental cues such as photoperiod. Notably, budbreak occurred approximately 20 days earlier in 2023 compared to 2022 due to higher accumulated degree days (ADDs). Short-rotation-coppice (SRC) management delayed budbreak by five to ten days. Cold damage was significant in 2023, particularly for genotypes from northern provenances and those with P. maximowiczii parentage. Severe damage was also observed in eastern cottonwood ( Populus deltoides ​× ​ Populus deltoides (D ​× ​D)) genotypes, despite most having southeastern US parentages. Leaf dynamics, including leaf duration and leaf area index (LAI), varied across taxa and sites, with earlier budbreak correlating with extended growing seasons and increased LAI. Biomass production was intricately linked to phenological events, with earlier budbreak leading to increased biomass production and greater susceptibility to cold damage. Our findings highlight the importance of genetics, environment, and coppicing management in understanding and managing Populus phenology and biomass production. These insights provide valuable guidance for developing effective breeding, conservation, and management strategies for Populus species in the context of climate change.

Accumulated degree days (ADDs)↗

Assessment of the impact of reactor residence time distribution on non-equilibrium product selectivity of polypropylene pyrolysis using reactive molecular dynamics simulations

Mass residence time distribution (RTD) is considered to be an important factor controlling the product selectivity in the pyrolysis of biomass and plastic wastes along with the pyrolysis chemistry. However, due to the complex pyrolysis chemistry of biomass and plastic waste, the coupling between the reaction chemistry, RTD, and product selectivity is challenging to understand. Here, we introduce a reaction molecular dynamics-based method to examine pyrolysis chemistry and species timescales to assess the impact of RTD on product selectivity and yield. To validate this method, reactive molecular dynamics simulations were conducted for polypropylene pyrolysis and its non-equilibrium product selectivity using 6 different RTDs. We find that the RTD and the reaction chemistry control the peak non-equilibrium product concentrations. The peak monomer (C 3 H 6 ) concentration during pyrolysis can be increased by up to 25 % by using a narrow RTD in the case of polypropylene pyrolysis. We also find that product selectivity is strongly affected by the average residence time and RTD. This coupling between the reaction chemistry, RTD, and product selectivity highlights the need to understand detailed reaction chemistry to control RTD and optimize non-equilibrium product selectivity during polymer and biomass pyrolysis. The present method provides a new way to design RTD for reactors to reach maximized product selectivity of plastic waste and biomass.

33 ADVANCED PROPULSION SYSTEMS↗

First Principles Based Reactive Atomistic Simulations to Understand the Effects of Molecular Hypervelocity Impact on Cassini's Ion and Neutral Mass Spectrometer

We report here on the predicted impact of species such as ice-water, CO2, CH4, and NH3, on oxidized titanium, as well as HC species on diamond surfaces. These simulations provide the dynamics of product distributions during and after a hypervelocity impact event, ionization fractions, and dissociation probabilities for the various species of interest as a function of impact velocity (energy). We are using these results to determine the relevance of the fragmentation process to Cassini INMS results, and to quantify its effects on the observed spectra.

hypervelocity collisions↗

Kinetic isotope effects for dissociative recombination of tritiated ketenyl ion ( 3 HCCO + ): A surface-hopping ab initio molecular dynamics study

Dissociative recombination (DR) reactions are important when modeling charged species in the presence of free electrons. While experimental measurements of DR reaction rates are challenging, surface hopping ab initio molecular dynamics (SH-AIMD) simulations provide an attractive alternative. SH-AIMD is especially well-suited for estimating branching ratios, i.e., the relative rates of competing production channels, for DR reactions. Although the radiolysis of diatomic tritium has been studied experimentally, previous attempts to model these systems have failed to account for isotope effects in DR reactions. Previous SH-AIMD studies have also not investigated tritium isotope effects for the branching ratios of DR reactions. In this study, we compute the DR branching ratios of the protiated and tritiated ketenyl ion. Comparison with literature values for the protiated branching ratios provides confidence in the reliability of our SH-AIMD results. Our simulations predict a significant increase of the HC + CO branching ratio for the tritiated system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗