A chemical dynamics study on the gas-phase formation of triplet and singlet C5H2 carbenes

Atomistic dynamics of elimination and substitution driven by entrance channel

E2 elimination and SN2 substitution reactions are of central importance in preparative organic synthesis due to their stereospecificity. Herein, atomistic dynamics of a prototype reaction of ethyl chloride with hydroxide ion are uncovered that show strikingly distinct features from the case with fluoride anion. Chemical dynamics simulations reproduce the experimental reaction rate and reveal that the E2 proceeding through a direct elimination mechanism dominates over SN2 for the hydroxide ion reaction. This unexpected finding of a pronounced contribution of direct reaction dynamics, even at a near-thermal energy, is in strong contrast to the complex-mediated indirect mechanism for the fluoride case that characterizes the low-energy ion-molecule reactions. The entrance channel structures are found to be crucial and the differences are attributed to subtle changes in the hydrogen-bonding interaction of the approaching reactants. This effect presents in E2/SN2 reactions of different bases and alkyl halides and might play a role in complex chemical networks and environments.

Theoretical Study of the Subsequent Decomposition Mechanisms of 1,1-Diamino-2,2-dinitroethene (FOX-7)

This computational study focuses on the mechanism of the consecutive decomposition of FOX-7 and compares the results with recent experimental study [J. Phys. Chem. A 2023, 127, 7707] under 202 nm photolysis (592 kJ/mol). The mechanisms of forming these compounds, including cyanamide variants (HNCNH and NH2CN), hydroxylamine (NH2OH), nitrosamine (NH2NO), diaminoacetylene (H2NCCNH2), cyanogen (NCCN), water (H2O), ammonia (NH3), urea ((NH2)2CO), hydroxyurea (NH2C(O)NHOH), and formamide (NH2CHO), have only been speculated on without any energetic information previously. This study employed an unsupervised potential energy profile search protocol and ab initio molecular dynamics (AIMD) simulations to identify reaction pathways leading to these compounds. The calculations reveal that although some products (e.g., HNCNH, NH2CN, H2NCCNH2, and NCCN) can be formed via unimolecular decomposition, other products (e.g., NH2OH, NH2NO, H2O, NH3, (NH2)2CO, NH2C(O)NHOH, and NH2CHO) are energetically favored if they are formed via bimolecular recombination between unimolecular decomposition products or a product and a FOX-7 molecule.

Energy and spectroscopic parameters of neutral and cations isomers of the CnH2 (n = 2-6) families using high-level ab-initio approaches

Cationic species, previously detected from ion-induced desorption of solid methane by plasma desorption mass spectrometry (PDMS), and neutral species, are investigated using high-level ab-initio approaches. From a set of 25 cationic and 26 neutral structures belonging to CnH2 (n = 2-6) families, it was obtained the energy, rotational constants, harmonic vibrational frequency, charge distribution and excitation energies. The ZPVE-corrected energies, at CCSD(T)-F12; CCSD(T)-F12/RI/(cc-pVTZ-F12, cc-pVTZ-F12-CABS, cc-pVQZ/C) (n = 2-5) and CCSD(T)/cc-pVTZ (n = 6) levels, reveal that the topology of the most stable isomer vary with n and the charge. Out of 674 harmonic frequencies, those with maximum intensity are generally in the 3000-3500 cm-1 range. Analysis of 169 vertical transition energies calculated with the EOM-CCSD approach, suggest three C6H2 species as potential carriers of the diffuse interstellar bands (DIB). Systematic comparison of properties between neutral and cationic species can assist in the structural description of complex matrices.

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

Ion molecule reactions in the HBr+ + CH4 system: a combined experimental and theoretical study

Reactions in the system HBr+ + CH4 have been investigated inside a guided ion-beam apparatus under single-collision conditions. The HBr+ is vibrational and rotational state selected in the electronic X2Π1/2 state created by (2+1)-REMPI. Due to the exitation scheme employed different rotational states of the HBr+ are accessible. Four reaction channels have been observed. The cross section, σ, for the exothermic proton transfer channel (PT) decreases with increasing collision energy, steeper than predicted by the Langevin model. The cross section also decreases with increasing rotational energy in the HBr+, with the effect of the rotational energy being stronger than that of translational energy. The cross section for the endothermic charge transfer (CT) increased with increasing collision energy. The energy dependence is well reproduced by a simple line of center (loc) model. Although the bromine transfer (BT) is exothermic the observed cross section increased with increasing collision energy due to an activation barrier on the potential energy surface (PES). Analysis by a modified loc model suggest the relevance of an angle dependence of σ. The cross section for the endothermic hydrogen atom abstraction (HA) exhibits a maximum at 2 eV Ecm. The measured cross sections are rationalized by means of reaction dynamics simulations which show good agreement with the experimental cross sections. The dynamics simulations are carried out with a machine learning potential that is developed and benchmarked with ab initio molecular dynamics simulation. The absolute cross sections predicted by reaction dynamics simulations are well within the same order of magnitude while reproducing the trends over three different collision energies for all four reaction channels. Furthermore, the simulations demonstrate various reaction mechanisms for these reaction channels, including a very interesting HBr+ orientation selectivity for the BT reaction channel.

CCSD(T) Rotational Constants for Highly Challenging C5H2 Isomers-A Comparison between Theory and Experiment

We evaluate the accuracy of CCSD(T) and density functional theory (DFT) methods for the calculation of equilibrium rotational constants (Ae, Be, and Ce) for four experimentally detected low-lying C5H2 isomers (ethynylcyclopropenylidene (2), pentatetraenylidene (3), ethynylpropadienylidene (5), and 2-cyclopropen-1-ylidenethenylidene (8)). The calculated rotational constants are compared to semi-experimental rotational constants obtained by converting the vibrationally averaged experimental rotational constants (A0, B0, and C0) to equilibrium values by subtracting the vibrational contributions (calculated at the B3LYP/jun-cc-pVTZ level of the theory). The considered isomers are closed-shell carbenes, with cumulene, acetylene, or strained cyclopropene moieties, and are therefore highly challenging from an electronic structure point of view. We consider both frozen-core and all-electron CCSD(T) calculations, as well as a range of DFT methods. We find that calculating the equilibrium rotational constants of these C5H2 isomers is a difficult task, even at the CCSD(T) level. For example, at the all-electron CCSD(T)/cc-pwCVTZ level of the theory, we obtain percentage errors ≤0.4% (Ce of isomer 3, Be and Ce of isomer 5, and Be of isomer 8) and 0.9-1.5% (Be and Ce of isomer 2, Ae of isomer 5, and Ce of isomer 8), whereas for the Ae rotational constant of isomers 2 and 8 and Be rotational constant of isomer 3, high percentage errors above 3% are obtained. These results highlight the challenges associated with calculating accurate rotational constants for isomers with highly challenging electronic structures, which is further complicated by the need to convert vibrationally averaged experimental rotational constants to equilibrium values. We use our best CCSD(T) rotational constants (namely, ae-CCSD(T)/cc-pwCVTZ for isomers 2 and 5, and ae-CCSD(T)/cc-pCVQZ for isomers 3 and 8) to evaluate the performance of DFT methods across the rungs of Jacob's Ladder. We find that the considered pure functionals (BLYP-D3BJ, PBE-D3BJ, and TPSS-D3BJ) perform significantly better than the global and range-separated hybrid functionals. The double-hybrid DSD-PBEP86-D3BJ method shows the best overall performance, with percentage errors below 0.5% in nearly all cases.

Ion-molecule reactions in the HBr+ + HCl (DCl) system: a combined experimental and theoretical study

Reactions in the system HBr⁺ + HCl (DCl) were investigated inside a guided ion-beam apparatus under single-collision conditions. In the HBr⁺ + HCl system, the proton transfer (PT_HCl) and charge transfer (CT) are observable. In the HBr⁺ + DCl system, proton transfer (PT_DCl) and deuterium abstraction (DA) are accessible. The cross sections for all reaction channels were measured as a function of the collision energy E_cm and of the rotational energy E_rot of the ion. The rotationally state-selective formation of the ionic species was realized by resonance-enhanced multiphoton ionization (REMPI). As expected, the PT-channels are exothermic, and the cross section decreases with increasing collision energy for both PT_HCl and PT_DCl. The cross section for DA also decreases with an increasing E_c.m.. In the case of a considerably endothermic CT-channel, the reaction efficiency increases with increasing collision energy but has an overall much smaller cross sections compared to PT and DA reactions. Both PT-reactions are hindered by ion rotation, whereas DA is independent of E_rot. The CT-channel shows a rotational enhancement near the thermochemical threshold. The experiment is complemented by theory, using ab initio molecular dynamics (AIMD, also known as direct dynamics) simulations and taking the rotational enhancement of HBr⁺ into account. The simulations show good agreement with the experimental results. The cross section of PT_HCl decreases with an increase of the rotational energy. Furthermore, the absolute cross sections are in the same order of magnitude. The CT channel shows no reactions in the simulation.

Why Are MgC3H Isomers Missing in the Interstellar Medium?

Considering the recent findings of linear doublet (2Σ+) MgCnH isomers (n = 2, 4, and 6) in the evolved carbon star IRC+10216, various structural isomers of MgC3H and MgC3H+ are theoretically investigated here. For MgC3H, 11 doublet and 8 quartet stationary points ranging from 0.0 to 71.8 and 0.0 to 110.1 kcal mol-1, respectively, have been identified initially at the UωB97XD/6-311++G(2d,2p) level. To get accurate relative energies, further energy evaluations are carried out for all isomers with coupled cluster methods and thermochemical modules such as G3//B3LYP, G4MP2, and CBS-QB3 methods. Unlike the even series, where the global minima are linear molecules with a Mg atom at one end, in the case of MgC3H, the global minimum geometry turns out to be a cyclic isomer, 2-magnesabicyclo[1.1.0]but-1,3,4-triyl (1, C2v, 2A1). In addition, five low-lying isomers, magnesium-substituted cyclopropenylidene (2, Cs, 2A'), 1-magnesabut-2,3-dien-1-yl-4-ylidene (3, Cs, 2A″), 1-magnesabut-2-yn-1-yl-4-ylidene (4, Cs, 2A″), 2λ3-magnesabicyclo[1.1.0]but-1,3-diyl-4-ylidene (5, C2v;, 2A1), and 1-magnesabut-2,3-dien-2-yl-4-ylidene (6, C∞v, 2Σ+), were also identified. The doublet linear isomer of MgC3H, 1-magnesabutatrienyl (10, C∞v, 2Σ+) turns out to be a minimum but lies 54.1 kcal mol-1 above 1 at the ROCCSD(T)/cc-pVTZ level. The quartet (4Σ+) electronic state of 10 was also found to be a minimum, but it lies 8.0 kcal mol-1 above 1 at the same level. Among quartets, isomer 10 is the most stable molecule. The next quartet electronic state (of isomer 11) is 34.4 kcal mol-1 above 10, and all other quartet electronic states of other isomers are not energetically close to low-lying doublet isomers 2 to 6. Overall, the chemical space of MgC3H contains more cyclic isomers (1, 2, and 3) on the low-energy side unlike their even-numbered MgCnH counterparts (n = 2, 4, and 6). Though the quartet electronic state of 10 is linear, it is not the global minimum geometry on the MgC3H potential energy surface. Isomerization pathways among the low-lying isomers (doublets of 1-4 and a quartet of 10) reveal that these molecules are kinetically stable. For the cation, MgC3H+, the cyclic isomers (1+, 2+, and 3+) are on the low-energy side. The singlet linear isomer, 10+, is a fourth-order saddle point. The low-lying cations are quite polar, with dipole moment values of >7.00 D. The current theoretical data would be helpful to both laboratory astrophysicists and radioastronomers for further studies on the MgC3H0/+ isomers.

From Molecules with a Planar Tetracoordinate Carbon to an Astronomically Known C5H2 Carbene

Ethynylcyclopropenylidene (2), an isomer of C₅H₂, is a known molecule in the laboratory and has recently been identified in Taurus Molecular Cloud-1 (TMC-1). Using high-level coupled-cluster methods up to the CCSDT(Q)/CBS level of theory, it is shown that two isomers of C₅H₂ with a planar tetracoordinate carbon (ptC) atom, (SP-4)-spiro[2.2]pent-1,4-dien-1,4-diyl (11) and (SP-4)-spiro[2.2]pent-1,4-dien-1,5-diyl (13), serve as the reactive intermediates for the formation of 2. Here, a theoretical connection has been established between molecules containing ptC atoms (11 and 13) and a molecule (2) that is present nearly 430 light years away, thus providing evidence for the existence of ptC species in the interstellar medium. The reaction pathways connecting the transition states and the reactants and products have been confirmed by intrinsic reaction coordinate calculations at the CCSDT(Q)/CBS//B3LYP-D3BJ/cc-pVTZ level. While isomer 11 is non-polar (μ = 0), isomers 2 and 13 are polar, with dipole moment values of 3.52 and 5.17 Debye at the CCSD(T)/cc-pVTZ level. Therefore, 13 is also a suitable candidate for both laboratory and radioastronomical studies.

Directed gas phase preparation of ethynylallene (H2CCCHCCH; X1A′) via the crossed molecular beam reaction of the methylidyne radical (CH; X2Π) with vinylacetylene (H2CCHCCH; X1A′)

The elementary reaction of the methylidyne radical with vinylacetylene leading to the predominant formation of ethynylallene and atomic hydrogen via indirect scattering dynamics.

A chemical dynamics study of the reaction of the methylidyne radical (CH, X2Π) with dimethylacetylene (CH3CCCH3, X1A1g)

The gas-phase reaction of the methylidyne (CH; X²Π) radical with dimethylacetylene (CH₃CCCH₃; X¹A_1g) was studied at a collision energy of 20.6 kJ mol^-1 under single collision conditions with experimental results merged with ab initio calculations of the potential energy surface (PES) and ab initio molecule dynamics (AIMD) simulations. The crossed molecular beam experiment reveals that the reaction proceeds barrierless via indirect scattering dynamics through long-lived C₅H₇ reaction intermediate(s) ultimately dissociating to C₅H₆ isomers along with atomic hydrogen with atomic hydrogen predominantly released from the methyl groups as verified by replacing the methylidyne with the D1-methylidyne reactant. AIMD simulations reveal that the reaction dynamics are statistical leading predominantly to p28 (1-methyl-3-methylenecyclopropene, 13%) and p8 (1-penten-3-yne, 81%) plus atomic hydrogen with a significant amount of available energy being channeled into the internal excitation of the polyatomic reaction products. The dynamics are controlled by addition to the carbon-carbon triple bond with the reaction intermediates eventually eliminating a hydrogen atom from the methyl groups of the dimethylacetylene reactant forming 1-methyl-3-methylenecyclopropene (p28). The dominating pathways reveal an unexpected insertion of methylidyne into one of the six carbon-hydrogen single bonds of the methyl groups of dimethylacetylene leading to the acyclic intermediate, which then decomposes to 1-penten-3-yne (p8). Therefore, the methyl groups of dimethylacetylene effectively 'screen' the carbon-carbon triple bond from being attacked by addition thus directing the dynamics to an insertion process as seen exclusively in the reaction of methylidyne with ethane (C₂H₆) forming propylene (CH₃C₂H₃). Therefore, driven by the screening of the triple bond, one propynyl moiety (CH₃CC) acts in four out of five trajectories as a spectator thus driving an unexpected, but dominating chemistry in analogy to the methylidyne - ethane system.

Directed Gas-Phase Formation of Aminosilylene (HSiNH2; 1A'): The Simplest Silicon Analogue of an Aminocarbene, under Single-Collision Conditions

The aminosilylene molecule (HSiNH₂, X¹A')-the simplest representative of an unsaturated nitrogen-silylene-has been formed under single collision conditions via the gas phase elementary reaction involving the silylidyne radical (SiH) and ammonia (NH₃). The reaction is initiated by the barrierless addition of the silylidyne radical to the nonbonding electron pair of nitrogen forming an HSiNH₃ collision complex, which then undergoes unimolecular decomposition to aminosilylene (HSiNH₂) via atomic hydrogen loss from the nitrogen atom. Compared to the isovalent aminomethylene carbene (HCNH₂, X¹A'), by replacing a single carbon atom with silicon, a profound effect on the stability and chemical bonding of the isovalent methanimine (H₂CNH)-aminomethylene (HNCH₂) and aminosilylene (HSiNH₂)-silanimine (H₂SiNH) isomer pairs is shown; i.e., thermodynamical stabilities of the carbene versus silylene are reversed by 220 kJ mol^-1. Hence, the isovalency of the main group XIV element silicon was found to exhibit little similarities with the atomic carbon revealing a remarkable effect not only on the reactivity but also on the thermochemistry and chemical bonding.

Six Low-Lying Isomers of C11H8 Are Unidentified in the Laboratory-A Theoretical Study

Isomers of C₁₁H₈ have been theoretically examined using density functional theory and coupled-cluster methods. The current investigation reveals that 2aH-cyclopenta[cd]indene (2), 7-ethynyl-1H-indene (6), 4-ethynyl-1H-indene (7), 6-ethynyl-1H-indene (8), 5-ethynyl-1H-indene (9), and 7bH-cyclopenta[cd]indene (10) remain elusive till date in the laboratory. The puckered low-lying isomer 2 lies at 9 kJ mol^-1 below the experimentally known molecule, cyclobuta[de]naphthalene (3), at the fc-CCSD(T)/cc-pVTZ//fc-CCSD(T)/cc-pVDZ level of theory. 2 lies at 36 kJ mol^-1 above the thermodynamically most stable and experimentally known isomer, 1H-cyclopenta[cd]indene (1), at the same level. It is identified that 1,2-H transfer from 1 yields 2H-cyclopenta[cd]indene (14) and subsequent 1,2-H shift from 14 yields 2. Appropriate transition states have been identified, and intrinsic reaction coordinate calculations have been carried out at the B3LYP/6-311+G(d,p) level of theory. Recently, 1-ethynyl-1H-indene (11) has been detected using synchrotron-based vacuum ultraviolet ionization mass spectrometry. 2-Ethynyl-1H-indene (4) and 3-ethynyl-1H-indene (5) have been synthetically characterized in the past. While the derivatives of 7bH-cyclopenta[cd]indene (10) have been isolated elsewhere, the parent compound remains unidentified till date in the laboratory. Although C₁₁H₈ is a key elemental composition of astronomical interest for the formation of polycyclic aromatic hydrocarbons in the interstellar medium, none of its low-lying isomers have been characterized by rotational spectroscopy though they are having a permanent dipole moment (μ ≠ 0). Therefore, energetic and spectroscopic properties have been computed, and the present investigation necessitates new synthetic studies on C₁₁H₈, in particular 2, 6-10, and also rotational spectroscopic studies on all low-lying isomers.