High-level CCSD Research Articles

Atmospheric Chemistry of Chloroprene Initiated by OH Radicals: Combined Ab Initio/DFT Calculations and Kinetics Analysis.

Chloroprene (CP; CH2═C(Cl)-CH═CH2) is a significant toxic airborne pollutant, often originating from anthropogenic activities. However, the environmental fate of CP is incompletely understood. High level CCSD(T)/aug-cc-pVTZ//M06-2X/aug-cc-pVTZ calculations combined with kinetic modeling were employed here to glean new insight into the reaction mechanism, energies, and kinetics of the reaction of CP with OH radical (•OH). We report the energies of four different addition pathways and six different abstraction pathways. The •OH attack on the terminal C1 atom of the =CH2 group (which is directly attached to the =CCl moiety), leading to the formation of HOCH2-•C(Cl)-CH═CH2, was found to be a major path. The barrier height for the formation of the corresponding transition state was found to be -1.9 kcal mol-1 below that of the starting CP + •OH reactants. Rate coefficients were calculated for addition and abstraction pathways involving the CP + •OH reaction under pre-equilibrium approximation conditions, employing a combination of canonical variational transition state theory and small curvature tunneling. The overall rate coefficient for the reaction of CP + •OH at 298 K was found to be 1.4 × 10-10 cm3 molecule-1 s-1. The thermochemistry of the possible channels and atmospheric implications are provided. In addition, the fate of HOCH2-•C(Cl)-CH═CH2 in the presence of 3O2 was investigated. We found the reaction of the CP-derived peroxy radical adduct with HO2 and NO to make contributions to the formation of products such as formaldehyde, HO2 radical, Cl atom, HOCH2C(OOH)(Cl)CH═CH2, HOCH2C(O)Cl, ClC(O)CH═CH2, HOCH2C(O)CH═CH2, and HC(O) radical.

Theoretical reaction kinetics predictions for acetone and ketene with NH2 radicals: Implications on acetone/ammonia kinetic modeling

The cross-reaction kinetics of acetone/ketene (CH3COCH3/CH2CO) + amino (NH2) radicals are first theoretically reported for a wide range of conditions (T = 300–2500 K and P = 0.1–100 atm) in this work. The high-level electronic structure method CCSD(T)/cc-pVnZ(n = T, Q)//B3LYP-D3BJ/6–311++G(d,p) is used to explore the potential energy profiles on which the temperature- and pressure-dependence kinetic behaviors of the title reactions are characterized using the Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) theory. Corrections of the one-dimensional hindered rotor approximation and asymmetric Eckart tunneling effect are also included in the rate constant calculations. Furthermore, this work delves into the competitive relationships among the H-abstraction, NH2 addition and addition-dissociation reaction pathways. For CH3COCH3 + NH2, the H-abstraction reaction of CH3COCH3 is highly favored over the addition-dissociation reaction and other isomerization channels. For CH2CO + NH2, the reaction channel of addition-dissociation to form CH2NH2 + CO under low temperature conditions is more advantageous, while the H-abstraction reaction plays a dominant role under combustion conditions (T ≥ 1000 K). To reveal the impact of the studied reaction kinetics on model predictions, the rate constants of dominant reaction channels calculated in this work are incorporated into a kinetic model for the auto-ignition, oxidation and pyrolysis of CH3COCH3/NH3 mixtures. The simulated results show that the effect of the updated rate constants on ignition delay time is mainly at low temperatures and the promotion effect of CH3COCH3 addition on the ignition of NH3 presents a nonlinear enhancement. The updated rate constants can also accelerate the consumption of CH3COCH3 and CH2CO formation, while also influence the concentration distributions of other significant species (e.g., NH3, CO). Therefore, the cross-reaction kinetics of CH3COCH3/CH2CO + NH2 are critical in controlling the fuel consumption, important intermediate formation and ignition delay time of CH3COCH3/NH3 mixtures.

CB11S3+: A Triangular Boron-Based Cluster with One Planar Tetracoordinate Carbon at Its Edge.

Boron-based clusters containing planar tetracoordinate carbon (ptC) are unique and scarce. Isoelectronic-replacing after proper vulcanization is an effective strategy to obtain the ptC based on the B12 cluster. Herein we report computational evidence for a ternary CB11S3+ (C2v, 1A1) cluster, which possesses a concentric double-triangle structure containing one ptC atom at the peripheral edge. The unbiased structural explorations of potential energy surfaces and high-level CCSD(T) calculations indicate that the ptC CB11S3+ cluster is a true global minimum. Born-Oppenheimer molecular dynamics (BOMD) simulations reveal that it is dynamically stable against isomerization and decomposition. Chemical bonding analysis reveals that three delocalized π bonds endow the π aromaticity to the CB11 unit of CB11S3+. In addition, the strong S → B π back-bonding is also conducive to the stability of CB11S3+. The current findings offer opportunities for further boron-based ptC clusters.

InnTl4-nH+ (n = 0∼4): Tetracoordinate Hydrogen in a Planar Fashion?

The recent report of planar tetracoordinate hydrogen (ptH) in In4H+ is very intriguing in planar hypercoordinate chemistry. Our high-level CCSD(T) calculations revealed that the proposed D4h-symmetric ptH In4H+ is a first-order saddle point with an imaginary frequency in the out-of-plane mode of the hydrogen atom. In fact, at the CCSD(T)/aug-cc-pV5Z/aug-cc-pV5Z-PP level, the C4v isomer, with the H atom located 0.70 Å above the In4 plane, is 0.5 kcal/mol more stable than the D4h isomer. However, given the small perturbation from planarity and essentially barrierless C4v ↔ D4h ↔ C4v transition, the vibrationally averaged structure can still be considered as a planar. Extending our exploration to the InnTl4-nH+ (n = 0-3) systems, we found all these ptH structures, except for In2Tl2H+, to be the putative global minimum. The single σ-delocalized interaction between the central hydrogen atom and InnTl4-n ligand rings proves pivotal in establishing planarity and aromaticity and conferring substantial stability upon these rule-breaking ptH species.

XB2Bi2 (X = Si, Ge, Sn, Pb): Penta-Atomic Planar Tetracoordinate Si/Ge/Sn/Pb Clusters with 20 Valence Electrons.

Planar tetracoordinate silicon, germanium, tin, and lead (ptSi/Ge/Sn/Pb) species are scarce and exotic. Here, we report a series of penta-atomic ptSi/Ge/Sn/Pb XB2Bi2 (X = Si, Ge, Sn, Pb) clusters with 20 valence electrons (VEs). Ternary XB2Bi2 (X = Si, Ge, Sn, Pb) clusters possess beautiful fan-shaped structures, with a Bi-B-B-Bi chain surrounding the central X core. The unbiased density functional theory (DFT) searches and high-level CCSD(T) calculations reveal that these ptSi/Ge/Sn/Pb species are the global minima on their potential energy surfaces. Born-Oppenheimer molecular dynamics (BOMD) simulations indicate that XB2Bi2 (X = Si, Ge, Sn, Pb) clusters are robust. Bonding analyses indicate that 20 VEs are perfect for the ptX XB2Bi2 (X = Si, Ge, Sn, Pb): two lone pairs of Bi atoms; one 5c-2e π, and three σ bonds (two Bi-X 2c-2e and one B-X-B 3c-2e bonds) between the ligands and X atom; three 2c-2e σ bonds and one delocalized 4c-2e π bond between the ligands. The ptSi/Ge/Sn/Pb XB2Bi2 (X = Si, Ge, Sn, Pb) clusters possess 2π/2σ double aromaticity, according to the (4n + 2) Hückel rule.

Structures and stability of I2 and ICl complexes with pyridine: Ab initio and DFT study.

Theoretical investigation of thermodynamic stability and bonding features of possible isomers of the molecular and ionic complexes of pyridine with molecular iodine and iodine monochloride IX (X = I,Cl) is presented. M06-2X DFT functional is found to provide bond distances and dissociation energies which are close to those obtained at high-level ab initio CCSD(T)/aug-cc-pvtz//CCSD/aug-cc-pvtz benchmark computations for the most stable isomers, formed via donation of a lone pair of nitrogen atom of pyridine to the iodine atom. These isomers are by 23-33 kJ mol-1 (in case of I2) and by 39-56 kJ mol-1 (in case of ICl) more stable than other molecular complexes. T-shaped π-σ* bonded isomers turn out to be energetically comparable with van der Waals bound compounds. Among the ionic isomers, structures featuring [IPy2]+ cation with I3 - or ICl2 - counterions are more stable. Oligomerization favors ionic isomers starting from the tetrameric clusters of the composition (IX)4Py4.

Using single and double laser pulses on the molecular Ni4@C48H36 system to design integrated nanospintronic units.

The accomplishment of long-distance spin transfer scenarios between several magnetic centers is a big challenge for building and supporting spin-logic units for developing future all-optical magnetic unit operations. Using high-level quantum chemistry theory CCSD and EOM-CCSD, we systematically study the ultrafast laser-induced spin-dynamics process on a carbon-based material, to which four magnetic centers are attached. We show that the CCSD method with the 6-31G basis set calculation is sensitive to the C-Ni bond length. The spin density distribution, which is computed using EOM-CCSD with LanL2DZ+ECP calculations, Mulliken population analysis, including spin-orbit-coupling (SOC) and a magnetic field, fulfills the requirements for achieving spin dynamics processes. Different local spin-flip and spin-transfer processes are accomplished within the subpicosecond regime. The impact of the propagation direction of the laser pulse by switching their polar and the azimuthal angles in spherical coordinates on the spin dynamics processes is analyzed. Double laser pulses with time delay δt ≥ 200 × FWHM yield in a realistic magnetic field gradient selectively a lateral resolution, which corresponds to distances smaller than the CMOS scale (2 nm in 2024) while our system size is comparable to the CMOS scale. Here Λ and V processes with two quasi-degenerate intermediate levels are used. We propose a model of an integrated spin-logic processor created from an array of individual spin-logic blocks, which are realized by four magnetic centers Ni. The findings of this study demonstrate the enormous potential of using laser-induced spin dynamics as the fundamental mechanism for future molecular magnetic technology.

Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory.

This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three-five water molecules (Ar(H2O)n=3-5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body-many-body (2b:Mb) and three-body-many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol-1 and by nearly 6 kJ mol-1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm-1 for some hydrogen-bonded OH stretching frequencies.

Ternary XBe4H5- (X = Si, Ge, Sn, Pb) Clusters: Planar Tetracoordinate Si/Ge/Sn/Pb Species with 18 Valence Electrons.

As one of the important probes of chemical bonding, planar tetracoordinate carbon (ptC) compounds have been receiving much attention. Compared with ptC clusters, the heavier planar tetracoordinate silicon, germanium, tin, lead (ptSi/Ge/Sn/Pb) systems are scarcer and more exotic. The 18-valence-electron (ve)-counting is one important guide, though not the only rule, for the design of planar tetra-, penta-coordinate carbon and silicon clusters. The 18ve ptSi/Ge system is very scarce and needs to be expanded. Based on the isoelectronic principle and bonding similarity between the Al atom and the BeH unit, inspired by the previously reported ptSi global minimum (GM) SiAl42-, a series of ternary 18 ve XBe4H5- (X = Si, Ge, Sn, Pb) clusters were predicted with the ptSi/Ge/Sn/Pb centers. Extensive density functional theory (DFT) global minimum searches and high-level CCSD(T) calculations performed herein indicated that these ptSi/Ge/Sn/Pb XBe4H5- (X = Si, Ge, Sn, Pb) clusters were all true GMs on their potential energy surfaces. These GMs of XBe4H5- (X = Si, Ge, Sn, Pb) species possessed the beautiful fan-shaped structures: XBe4 unit can be stabilized by three peripheries bridging H and two terminal H atoms. It should be noted that XBe4H5- (X = Si, Ge, Sn, Pb) were the first ternary 18 ve ptSi/Ge/Sn/Pb species. The natural bond orbital (NBO), canonical molecular orbitals (CMOs) and adaptive natural densitpartitioning (AdNDP) analyses indicated that 18ve are ideal for these ptX clusters: delocalized one π and three σ bonds for the XBe4 core, three Be-H-Be 3c-2e and two Be-H σ bonds for the periphery. Additionally, 2π plus 6σ double aromaticity was found to be crucial for the stability of the ptX XBe4H5- (X = Si, Ge, Sn, Pb) clusters. The simulated photoelectron spectra of XBe4H5- (X = Si, Ge, Sn, Pb) clusters will provide theoretical basis for further experimental characterization.

Computational studies of HCCCCS isomers

The recent identification of interstellar molecules HCCCHCCC and HCCCCS in TMC-1 has prompted an ab initio investigation of these isomers. The isomers of HCCCHCCC have been studied previously but not so much HCCCCS. In this paper we carry out high-level CCSD(T) calculations to characterise the lowest eleven isomers of HCCCCS. We also apply complete basis set extrapolation including core correlation and relativistic effects to further investigate the key isomers and provide high-precision rotational constants which will hopefully aid astronomers in their assignments. Source: This image has been adapted from http://www.eso.org/public/images/eso1209d/ licensed under CC BY-SA 3.0.

CAl4X4 (X = Te, Po): Double Aromatic Molecular Stars Containing Planar Tetracoordinate Carbon Atoms

Planar tetracoordinate carbon (ptC) species are scarce and exotic. Introducing four peripheral Te/Po auxiliary atoms is an effective strategy to flatten the tetrahedral structure of CAl4 (Td, 1A1). Neutral CAl4X4 (X = Te, Po) clusters possess quadrangular star structures containing perfect ptC centers. Unbiased density functional theory (DFT) searches and high-level CCSD(T) calculations suggest that these ptC species are the global minima on the potential energy surfaces. Bonding analyses indicate that 40 valence-electron (VE) is ideal for the ptC CAl4X4 (X = Te, Po): one delocalized π and three σ bonds for the CAl4 core; four lone pairs (LPs) of four X atoms, eight localized Al-X σ bonds, and four delocalized Al-X-Al π bonds for the periphery. Thus, the ptC CAl4X4 (X = Te, Po) clusters possess the stable eight electron structures and 2π + 6σ double aromaticity. Born-Oppenheimer molecular dynamics (BOMD) simulations indicate that neutral ptC CAl4X4 (X = Te, Po) clusters are robust.

Probing the Potential Energy Profile of the I + (H2O)3 → HI + (H2O)2OH Forward and Reverse Reactions: High Level CCSD(T) Studies with Spin-Orbit Coupling Included.

Three different pathways for the atomic iodine plus water trimer reaction I + (H2O)3 → HI + (H2O)2OH were preliminarily examined by the DFT-MPW1K method. Related to previous predictions for the F/Cl/Br + (H2O)3 reactions, three pathways for the I + (H2O)3 reaction are linked in terms of geometry and energetics. To legitimize the results, the "gold standard" CCSD(T) method was employed to investigate the lowest-lying pathway with the correlation-consistent polarized valence basis set up to cc-pVQZ(-PP). According to the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the I + (H2O)3 → HI + (H2O)2OH reaction is predicted to be endothermic by 47.0 kcal mol-1. The submerged transition state is predicted to lie 43.7 kcal mol-1 above the separated reactants. The I···(H2O)3 entrance complex lies below the separated reactants by 4.1 kcal mol-1, and spin-orbit coupling has a significant impact on this dissociation energy. The HI···(H2O)2OH exit complex is bound by 4.3 kcal mol-1 in relation to the separated products. Compared with simpler I + (H2O)2 and I + H2O reactions, the I + (H2O)3 reaction is energetically between them in general. It is speculated that the reaction between the iodine atom and the larger water clusters may be energetically analogous to the I + (H2O)3 reaction. The iodine reaction I + (H2O)3 is connected with the analogous valence isoelectronic bromine/chlorine reactions Br/Cl + (H2O)3 but much different from the F + (H2O)3 reaction. Significant difference with other halogen systems, especially for barrier heights, are seen for the iodine systems.

Hydrogen Bonding versus Halogen Bonding: Spectroscopic Investigation of Gas-Phase Complexes Involving Bromide and Chloromethanes.

Hydrogen bonding and halogen bonding are important non-covalent interactions that are known to occur in large molecular systems, such as in proteins and crystal structures. Although these interactions are important on a large scale, studying hydrogen and halogen bonding in small, gas-phase chemical species allows for the binding strengths to be determined and compared at a fundamental level. In this study, anion photoelectron spectra are presented for the gas-phase complexes involving bromide and the four chloromethanes, CH3 Cl, CH2 Cl2 , CHCl3 , and CCl4 . The stabilisation energy and electron binding energy associated with each complex are determined experimentally, and the spectra are rationalised by high-level CCSD(T) calculations to determine the non-covalent interactions binding the complexes. These calculations involve nucleophilic bromide and electrophilic bromine interactions with chloromethanes, where the binding motifs, dissociation energies and vertical detachment energies are compared in terms of hydrogen bonding and halogen bonding.

Laser-induced ultrafast spin-transfer processes in non-linear zigzag carbon chain systems.

We combine the high-level quantum chemistry theory CCSD and EOM-CCSD together with local and global Λ processes to investigate the details of the laser-induced ultrafast spin manipulation scenarios in non-linear zigzag carbon chain systems Ni2@C32H32 and Ni2@C36H36. The spin density distribution, which is calculated on each many-body state using a Mulliken population analysis, fulfills the requirements to accomplish the spin dynamics processes. Various spin-flip and spin-transfer scenarios are accomplished. All the spin-dynamics processes can be achieved within subpicosecond times. Under the influence of a magnetic field, we find that the spin-transfer scenarios are preserved, while the local spin-flip scenario on a Ni atom can be significantly inhibited depending on the strength of the magnetic field. The impact of the propagation direction of the laser pulse on the spin dynamics processes by varying their polar and azimuthal angles in spherical coordinates is investigated. Additionally, we find that double laser pulses successfully induce the spin-transfer processes. Our outcomes underline the significant potential of carbon chain systems as building blocks for developing future all-optical integrated logic processing units.

Predicting Solvent Effects on SN2 Reaction Rates: Comparison of QM/MM, Implicit, and MM Explicit Solvent Models.

Solvents are one of the key variables in the optimization of a synthesis yield or properties of a synthesis product. In this paper, contemporary solvent models are applied to predict the rates of SN2 reactions in a range of aqueous and non-aqueous solvents. High-level CCSD(T)/CBS//M06-2X/6-31+G(d) gas phase energies were combined with solvation free energies from SMD, SM12, and ADF-COSMO-RS continuum solvent models, as well as molecular mechanics (MM) explicit solvent models with different atomic charge schemes to predict the rate constants of three SN2 reactions in eight protic and aprotic solvents. It is revealed that the prediction of rate constants in organic solvents is not necessarily less challenging than in water and popular solvent models struggle to predict their rate constants to within 3 log units of experimental values. Among the continuum solvent models, the ADF-COSMO-RS model performed the best in predicting absolute rate contants while the SM12 model was best at predicting relative rate constants with an average accuracy of about 1.5 and 0.8 log units, respectively. The use of computationally more demanding MM explicit solvent models did not translate to improvements in absolute rate constants but was quite effective at predicting relative rate constants due to systematic error cancellation. Free energy barriers obtained from umbrella sampling with explicit solvent QM/MM simulations led to excellent agreement with experimental values, provided that a validated level of theory is used to treat the QM region.

Spectroscopic Study of the Br- +CH3 I→I- +CH3 Br SN 2 Reaction.

Mass spectrometry and anion photoelectron spectroscopy have been used to study the gas-phase reaction involving and . The anion photoelectron spectra associated with the reaction intermediates of this reaction are presented. High-level CCSD(T) calculations have been utilised to investigate the reaction intermediates that may form as a result of the reaction along various different reaction pathways, including back-side attack and front-side attack. In addition, simulated vertical detachment energies of each reaction intermediate have been calculated to rationalise the photoelectron spectra.

Halides and the carbon-carbon double bond: Interactions of ethylene with bromide and iodide

A combined experimental and theoretical approach has been used to study the bromide-ethylene and iodide-ethylene gas-phase complexes. Experimental anion photoelectron spectra show the electron stabilisation of the complexes to be 0.19 eV and 0.14 eV associated with the bromide and iodide complexes respectively. High-level CCSD(T) calculations predict two minima with respect to each halide; the halide either appends linearly to a single hydrogen of the ethylene molecule or appends orthogonally to the carbon–carbon double bond. The two structural motifs are found to be close in energy, with the difference in dissociation energy associated with both halides determined to be approximately 1.3 kJ mol−1.

M(BO)k+1- Anions: Novel Superhalogens Based on Boronyl Ligands.

Superhalogens have been a subject of continuous attraction in the last four decades due to their unusual structures and interesting applications. In the quest for new superhalogens based on boronyl (BO) ligands, we investigate M(BO)k+1- anions using a high-level ab initio CCSD(T) and B3LYP density functional method for M = Li, Na, K, Be, Mg, Ca, B, and Al. Their structures are linear for M = Li, Na, K; trigonal planar for M = Be, Mg, Ca; and tetrahedral for M = B, Al. These anions are energetically and thermodynamically stable against various fragmentations. Their superhalogen property has been established due to their higher vertical detachment energy (VDE) than halogen, being in the range 4.44-7.81 eV. For a given k, the stability and VDE of M(BO)k+1- anions decrease with an increase in the size of M, which is due to the decrease in charge delocalization on BO moieties. There exists a linear correlation between frontier orbitals energy gap and VDE with coefficient R2 = 0.93888. It was also noticed that the BO ligand, due to the lower dimerization energy, can be preferably used as an inorganic analogue of CN. However, the structural relaxation in BO-based neutral species significantly affect the VDE and, hence, their superhalogen nature. We believe that these findings should be interesting for researchers working in superatomic chemistry as well as in boron chemistry.

Gas Phase Reaction of Isocyanic Acid: Kinetics, Mechanisms, and Formation of Isopropyl Aminocarbonyl.

Isocyanic acid, HNCO, mainly emitted by combustion processes, is doubted to be detrimental to human health if its concentration surpasses ∼1 ppbv. Very little information has been found regarding the HNCO loss in the gas phase. This study aims to close this knowledge gap by performing a theoretical kinetic study on the reaction of HNCO with the propargyl radical. The potential energy surface of the HNCO + C3H3 reaction was characterized utilizing high-level CCSD(T)/CBS(TQ5)//B3LYP/6-311++G(3df,2p) quantum-chemical approaches, followed by TST and RRKM/ME kinetic computations. The obtained results reveal that the reaction can proceed via H-abstraction, leading to the C3H4 + NCO bimolecular products with energy barriers of 23–25 kcal/mol, and/or addition, resulting in C4H4NO intermediates with 23–26 kcal/mol barrier heights. The C4H4NO adducts when formed can decompose to products and/or return to HNCO + C3H3 in which the reverse decompositions are found to be dominant with a branching ratio that accounts for nearly 100% at 300 K and 760 Torr. The calculated P-independent rate coefficients indicate that at low temperatures, the H-abstraction channels are insignificant. However, at high temperatures (T > 1500 K), the H-abstraction path leading to H3CCCH + NCO prevails with a branching ratio of ∼50–53% in the descending 1800–1500 K temperature range at 760 Torr, while the H-abstraction leading to H2CCCH2 + NCO is favorable at T > 1800 K, with the yield reaching above 50% at 760 Torr. In contrast to the H-abstraction rate constants, the calculated values for the additions and the C4H4NO decompositions show a positive pressure dependence. Both the total rate constants for the reactions HNCO + C3H3 → products and C4H4NO → products, which are, respectively, k_total_bimo(T) = 3.53 × 10–23T3.27 exp[(−21.35 ± 0.06 kcal/mol)/RT] cm3 molecule–1 s–1 and k_total_uni(T) = 1.13 × 1025T–4.02 exp[(−11.77 ± 0.16 kcal/mol)/RT] s–1, increase with the increasing temperature in the 300–2000 K range at 760 Torr. The rate constant of HNCO + C3H3 → products is about 8 orders of magnitude smaller than the value of HCHO + C3H3 → products, showing that HCHO is more reactive toward the C3H3 free radicals than HNCO. The computed heats of formation for several species agree well with the available literature data with the deviation less than 1.0 kcal/mol, indicating that the methods used in this study are extremely reliable. With the given results, it is vigorously suggested that the predicted rate constants, together with the thermodynamic data of the species involved, can be confidently used for modeling HNCO-related systems under atmospheric and combustion conditions.

σ0π2 Singlet Ground State Carbenes Undergo Least-Motion Reactions with H2 and Alkenes.

Ground state singlet carbenes commonly feature σ2π0 orbital occupations and are known for their concerted σ-bond insertion and cycloaddition reactions. Despite the facility of these transformations, orbital symmetry conservation forces them into non-least-motion π-approach reaction pathways. This situation completely changes when the singlet σ0π2 electron configuration becomes the ground state, which we show here by means of high-level CCSD(T) geometry optimizations. Carbenes like the experimentally known 2H-imidazol-2-ylidene react with H2 and ethylene with negligible or no barrier in a σ-fashion, which effectively corresponds to a least-motion reaction trajectory.