Cavity-Modified Unimolecular Dissociation Reactions via Intramolecular Vibrational Energy Redistribution
Not Available
SEARCH · Search NASA
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.
Not Available
Not Available
Results of a study showing that a metastable argon-carbon dioxide reaction results in dissociation of carbon dioxide and electronic excitation of one of the products, carbon monoxide or oxygen. A flow system using a 2450-MHz discharge was used to produce metastable argon atoms. Metastable argon in the afterglow was confirmed by adding nitrogen to the afterglow. Without addition of carbon dioxide no argon line emission, or any other emission, is observed from the reaction zone. Absence of argon line emission produced by recombination indicates the absence of charged species.
Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.
Abstract Field evaporation from ionic or covalently bonded materials often leads to the emission of molecular ions. The metastability of these molecular ions, particularly under the influence of the intense electrostatic field (10 10 Vm −1 ), makes them prone to dissociation with or without an exchange of energy amongst them. These processes can affect the analytical performance of atom probe tomography (APT). For instance, neutral molecules formed through dissociation may not be detected at all or with a time of flight no longer related to their mass, causing their loss from the analysis. Here, we evaluated the changes in the measured composition of FeO, Fe 2 O 3 and Fe 3 O 4 across a wide range of analysis conditions. Possible dissociation reactions are predicted by density-functional theory calculations considering the spin states of the molecules. The energetically favoured reactions are traced on to the multi-hit ion correlation histograms, to confirm their existence within experiments, using an automated Python-based routine. The detected reactions are carefully analyzed to reflect upon the influence of these neutrals from dissociation reactions on the performance of APT for analysing iron oxides.
Spontaneous chemical reactivity at multivalent (Mg, Ca, Zn, Al) electrode surfaces is critical to solid electrolyte interphase (SEI) formation, and hence, directly affects the longevity of batteries. Here, we report an investigation of the reactivity of 0.5 M Mg(TFSI)2 in 1,2-dimethoxyethane (DME) solvent at a Mg(0001) surface using ab initio molecular dynamics (AIMD) simulations and detailed Bader charge analysis. Based on the simulations, the initial degradation reactions of the electrolyte strongly depend on the structure of the Mg(TFSI)2 species near the anode surface. At the surface, the dissociation of Mg(TFSI)2 species occurs via cleavage of the N-S bond for the solvent separated ion pair (SSIP) and via cleavage of the C-S bond for the contact ion pair (CIP) configuration. In the case of the CIP, both TFSI anions undergo spontaneous bond dissociation reactions to form atomic O, C, S, F, and N species adsorbed on the surface of the Mg anode. These products indicate that the initial SEI layer formed on the surface of the pristine Mg anode consists of a complex mixture of multiple components such as oxides, carbides, sulfides, fluorides, and nitrides. We believe that the atomic level insights gained from these simulations will lay the groundwork for the rational design of tailored and functional interphases that are critical for the success of multivalent battery technology.
Laser flash photolysis of Cl2/CO/M mixtures (M = N2, CO, Ar, CO2) has been employed in conjunction with Cl((sup 2)P(sub J)) detection by time-resolved resonance fluorescence spectroscopy to investigate equilibration kinetics in the reactions Cl((sup 2)P(sub J)) + CO ClCO as a function of temperature (185-260 K) and pressure (14-200 Torr). The association and dissociation reactions are found to be in the low-pressure limit over the range of experimental conditions investigated. In N2 and/or CO buffer gases, the temperature dependences of the ClCO formation and dissociation reaction rate constants are described by the Arrhenius expressions k(sub 1) = (1.05 +/- 0.36) x 10(exp -34) exp[(810 +/- 70)/T] cm(exp 6)/molecules(exp 2).s and k(sub -1) = (4.1 +/- 3.1) x 10(exp -10) exp[(-2960 +/- 60)/T]cu cm/(molecule.s) (errors are 2 sigma). Second- and third-law analyses of the temperature dependence of the equilbrium constant (k/k-1) lead to the following thermodynamic parameters for the association reaction: Delta-H(sub 298) = -7.7 +/- 0.6 kcal/mol, Delta-H(sub 0) = -6.9 +/- 0.7 kcal/mol, Delta-S(sub 298) = -23.8 +/- 2.0 cal/mole.K, Delta-H(sub f,298)(ClCO) = 5.2 +/- 0.6 kcal/mol (errors are 2 sigma). The results repported in this study significantly reduce the uncertainties in all reported kinetic and thermodynamic parameters.
Conspectus (Photo)electrochemical energy conversion is important in the development of a carbon-neutral energy economy because it can provide a pathway for mitigating the intermittency of renewable energy sources such as wind and solar. In order to operate efficiently, these technologies, which include photoelectrochemical cells, water and CO 2 electrolyzers, fuel cells, and redox flow batteries, require fast charge transfer kinetics at the electrode/electrolyte interface as well as robust ion management in the electrolyte. In conventional electrolyzers and fuel cells, the electrolyte is strongly acidic or basic and ionic current is carried by H + or OH – ions. In contrast, photoelectrodes and electrocatalysts for water splitting are often studied in buffered solutions. The question of ion balance in these systems led us to analyze the polarization losses due to ion concentration gradients in cells that employed various buffer–membrane combinations. Continuously driving the buffer ions across an ionomer membrane not only lowers the buffer capacity of an aqueous electrolyte but also introduces pH gradients that result in significant energy losses. To address the problem, we and other groups have studied the use of reverse-biased bipolar membranes (BPMs) in (photo)electrolytic cells. BPMs consist of an anion exchange layer (AEL) laminated with a cation exchange layer (CEL) and are usually equipped with a catalytic layer in between to accelerate the water dissociation reaction. At the AEL/CEL interface, water dissociates into protons and hydroxide ions, which replenish those consumed at the cathode and anode. Compared to conventional water electrolyzers with proton/anion exchange membranes (PEM/AEM), BPM electrolyzers provide the unique advantage of continuously operating the cathode and anode under different pH conditions, which is desirable when the two electrode reactions have different pH requirements. BPMs also enable the use of buffered electrolytes at pH values that are optimized for electrode stability and product selectivity in applications such as CO 2 electrolysis. Product crossover losses and CO 2 pumping can be dramatically reduced in BPM-based CO 2 electrolyzers, relative to conventional alkaline membranes, by electrostatic repulsion (of anionic products) and electroosmotic drag (for neutral products). BPM-based gas fed CO 2 electrolyzers can achieve high current density, but they suffer from low Faradaic efficiency (FE) due to the acidic local environment of the CEL. This problem can be mitigated by adding an aqueous buffering layer or by creating a weak acid cation exchange film on the CEL face of the membrane. The use of BPMs in fuel cells and redox flow batteries offers some interesting advantages. Configurations with both reverse and forward bias have been studied, but forward bias has been favored due to material compatibility, reaction kinetics, and thermodynamic considerations. The net reaction at the AEL/CEL interface is the acid–base neutralization reaction, which has a high inherent reaction rate constant, but in the BPM is limited to a nanometric space-charge layer and requires efficient catalysis to achieve high current density. Understanding the mechanism of the acid–base neutralization and the opposite process, the water dissociation reaction, will be essential for improving the performance of forward-biased BPMs. In conclusion, this Account reviews our current understanding of the working mechanisms of BPMs and discusses how we can use them to effectively manage ions for various (photo)electrochemical applications.
In this report, we focus on the calculations of the energetics and chemical kinetics of heterogeneous reactions for Organometallic vapor phase epitaxy (OMVPE). The work described in this report builds upon our own previous thermochemical and chemical kinetics studies. The first of these articles refers to the prediction of thermochemical properties, and the latter one deals with the prediction of rate constants for gaseous homolytic dissociation reactions. The calculations of this investigation are at the microscopic level. The systems chosen consisted of a gallium nitride (GaN) substrate, and molecular nitrogen (N2) and ammonia (NH3) as adsorbants. The energetics for the adsorption and the adsorbant dissociation processes were estimated, and reaction rate constants for the dissociation reactions of free and adsorbed molecules were predicted. The energetics for substrate decomposition was also computed. The ONIOM method, implemented in the Gaussian98 program, was used to perform the calculations. This approach has been selected since it allows dividing the system into two layers that can be treated at different levels of accuracy. The atoms of the substrate were modeled using molecular mechanics6 with universal force fields, whereas the adsorbed molecules were approximated using quantum mechanics, based on density functional theory methods with B3LYP functionals and 6-311G(d,p) basis sets. Calculations for the substrate were performed in slabs of several unit cells in each direction. The N2 and NH3 adsorbates were attached to a central location at the Ga-lined surface.
The equilibrium energy, enthalpy, entropy, specific heat at constant volume and constant pressure, and the equation of state of the gas are all derived for chemically reacting gas mixtures in terms of the compressibility, the mol fractions, the thermodynamic properties of the pure gas components, and the change in zero point energy due to reaction. Results are illustrated for a simple diatomic dissociation reaction and nitrogen is used as an example. Next, a gas mixture resulting from combined diatomic dissociation and atomic ionization reactions is treated and, again, nitrogen is used as an example. A short discussion is given of the additional complexities involved when precise solutions for high-temperature air are desired, including effects caused by NO produced in shuffle reactions and by other trace species formed from CO2, H2O and Ar found in normal air.
Chemical reactions at high temperatures have been considered extensively because of their importance to the heating effects on re-entry of space vehicles. Data on these reactions however, are not abundant and even when found there are discrepancies in data collected by various investigators. In particular, data for recombination reactions are calculated from the dissociation reactions or vice versa through the equilibrium constant. This involves the use of the principle of detailed balancing. This principle is discussed in reference to conditions where it is valid as well as to those where it is not valid. Related topics that merit further study or for which applicable information was available are briefly mentioned in an appendix to this report.
In the nonequilibrium atmospheric chemistry regime that will be encountered by the proposed Aeroassisted Orbital Transfer Vehicle in the upper atmosphere, where air density is too low for thermal and chemical equilibrium to be maintained, the detailed high temperature air chemistry plays a critical role in defining radiative and convective heating loads. Although vibrational and electronic temperatures remain low (less than 15,000 K), rotational and translational temperatures may reach 50,000 K. Attention is presently given to the effects of multiple temperatures on the magnitudes of various chemical reaction rate constants, for the cases of both bimolecular exchange reactions and collisional excitation and dissociation reactions.
Organic compounds found in meteorites often show isotopic signatures of their interstellar/protosolar heritage as enrichments in D and 15N. Meteoritic organics found to be enriched in D include amino acids, hydroxy and dicarboxylic acids, as well as polycyclic aromatic hydrocarbons (PAHs). Processes that can produce isotopic enrichments in presolar/protosolar materials include gas-phase ion-molecule reactions, gas-grain surface reactions, and unimolecular photo-dissociation reactions involving PAHs. Because many molecules in interstellar clouds are enriched in D, the presence of D anomalies in meteorites is thought to originate from preserved or slightly altered interstellar/protostellar materials. However, the link between isotopic enrichments seen in space and those in meteoritic compounds and their relationship remain unclear. In this work, we present results of hydrogen isotopic fractionation for compounds in organic residues produced from the UV irradiation using an H2- discharge UV lamp of H2O:pyrimidine = 20:1 and H2O:NH3:pyrimidine = 20:2:1 ice mixtures at low temperature (is less than 20 K). After irradiation, the resulting residues are dissolved in H2O and analyzed with gas chromatography-mass spectrometry coupled with isotope ratio mass spectrometry (GC-MS/IRMS) [1], following a protocol similar to that used for previous analyses of comparable samples [2,3]. We used this technique to measure compound-specific D/H isotopic ratios for the initial pyrimidine and for two photo-products present in the residues, namely, 2,2'-bipyrimidine and an unidentified bipyrimidine isomer [2-4]. Measuring D enrichments in bipyrimidines has the advantage that the H atoms on these molecules are not easily exchangeable with other compounds, in particular the H2O and NH3 present in the ices or the solvents used to extract the samples for GC-MS/IRMS measurements. The δD value for the initial pyrimidine, measured with a high-temperature conversion elemental analyzer connected to the IRMS, was found to be -30% per mille. Preliminary measurements made on a residue produced from the UV irradiation of an H2O:NH3:pyrimidine = 20:2:1 ice mixture indicate δD values of +118% per mille for 2,2'- bipyrimidine and +92% per mille for the other bipyrimidine isomer, and therefore show a significant D enrichment during the photo-processing and warm-up that lead to their formation [5]. New measurements are currently being performed on a number of residues produced from simpler H2O:pyrimidine = 20:1 ice mixtures under different experimental conditions and will be presented here.
The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.
Here, photodissociation of pyruvic acid (PA) was studied in the gas-phase at 193 nm using two complementary techniques. The time-sliced velocity map imaging arrangement was used to determine kinetic energy release distributions of fragments and estimate dissociation timescales. The multiplexed photoionization mass spectrometer setup was used to identify and quantify photoproducts, including isomers and free radicals, by their mass-to-charge ratios, photoionization spectra, and kinetic time profiles. Using these two techniques, it is possible to observe the major dissociation products of PA photodissociation: CO 2 , CO, H, OH, HCO, CH 2 CO, CH 3 CO, and CH 3 . Acetaldehyde and vinyl alcohol are minor primary photoproducts at 193 nm, but products that are known to arise from their unimolecular dissociation, such as HCO, H 2 CO, and CH 4 , are identified and quantified. A multivariate analysis that takes into account the yields of the observed products and assumes a set of feasible primary dissociation reactions provides a reasonable description of the photoinitiated chemistry of PA despite the necessary simplifications caused by the complexity of the dissociation. These experiments offer the first comprehensive description of the dissociation pathways of PA initiated on the S 3 excited state. Most of the observed products and yields are rationalized on the basis of three reaction mechanisms: (i) decarboxylation terminating in CO 2 + other primary products (~50%); (ii) Norrish type I dissociation typical of carbonyls (~30%); and (iii) O—H and C—H bond fission reactions generating the H atom (~10%). The analysis shows that most of the dissociation reactions create more than two products. This observation is not surprising considering the high excitation energy (~51 800 cm –1 ) and fairly low energy required for dissociation of PA. We find that two-body fragmentation processes yielding CO 2 are minor, and the expected, unstable primary co-fragment, methylhydroxycarbene, is not observed because it probably undergoes fast secondary dissociation and/or isomerization. Norrish type I dissociation pathways generate OH and only small yields of CH 3 CO and HOCO, which have low dissociation energies and further decompose via three-body fragmentation processes. Experiments with d 1 -PA (CH 3 COCOOD) support the interpretations. The dissociation on S 3 is fast, as indicated by the products’ recoil angular anisotropy, but the roles of internal conversion and intersystem crossing to lower states are yet to be determined.
Understanding the time-dependent effects of charged vacancies on the electromechanical responses of materials is at the forefront of research for designing materials exhibiting metal-insulator transitions and memristive behavior. A Rayleighian approach is used to develop a model for studying the nonlinear kinetics of the reaction leading to generation of vacancies and electrons via the dissociation of vacancy-electron pairs. Also, diffusion and elastic effects of charged vacancies are considered to model polarization-electric potential and strain-electric potential hysteresis loops. The model captures multiphysics phenomena by introducing couplings among polarization, the electric potential, stress, strain, and concentrations of charged (multivalent) vacancies and electrons (treated as classical negatively charged particles), where the concentrations can vary due to association-dissociation reactions. A derivation of coupled time-dependent equations based on the Rayleighian approach is presented. Three limiting cases of the governing equations are considered, highlighting the effects of (1) nonlinear reaction kinetics on the generation of charged vacancies and electrons, (2) Vegard's law (i.e., the concentration-dependent local strain) on asymmetric strain-electric potential relations, and (3) coupling between a fast component and the slow component of the net polarization on the polarization-electric-field relations. The Rayleighian approach discussed in this work should pave the way for developing a multiscale modeling framework in a thermodynamically consistent manner while capturing multiphysics phenomena in ferroelectric materials. Published by the American Physical Society 2025
Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.
BCl 3 is a promising acceptor precursor for atomic-precision δ-doping of silicon, as it has been observed to rapidly dissociate into boron doped into the silicon surface and surface chlorine, which can be removed upon annealing. The chemical pathway and the resulting kinetics, through which BCl 3 adsorbs and dissociates on silicon, however, have only been partially explained. Here, in this work, we use density functional theory to expand the dissociation reactions of BCl 3 to include reactions that take place across multiple silicon dimer rows and reactions which end in a bare B atom either at the surface, substituted for a surface silicon, or in a subsurface position. We further simulate the resulting scanning tunneling microscopy images for each of these BCl x dissociation fragments, demonstrating that they often display distinct features that may allow for relatively confident experimental identification. Finally, we input the full dissociation pathway for BCl 3 into a kinetic Monte Carlo model, which simulates realistic reaction pathways as a function of environmental conditions, such as the pressure and temperature of dosing. We find that BCl 2 is broadly dominant at low temperatures, while high temperatures and ample space on the silicon surface for dissociation encourage the formation of bridging BCl fragments and B substitutions on the surface. This work provides the chemical mechanisms for understanding atomic-precision doping of Si with B, enabling a number of relevant quantum applications, such as bipolar nanoelectronics, acceptor-based qubits, and superconducting Si.