Search

We present an efficient, open-source formulation for coupled-cluster theory through perturbative triples with domain-based local pair natural orbitals [DLPNO-CCSD(T)]. Similar to the implementation of the DLPNO-CCSD(T) method found in the ORCA package, the most expensive integral generation and contraction steps associated with the CCSD(T) method are linear-scaling. In this work, we show that the t1-transformed Hamiltonian allows for a less complex algorithm when evaluating the local CCSD(T) energy without compromising efficiency or accuracy. Our algorithm yields sub-kJ mol−1 deviations for relative energies when compared with canonical CCSD(T), with typical errors being on the order of 0.1 kcal mol−1, using our TightPNO parameters. We extensively tested and optimized our algorithm and parameters for non-covalent interactions, which have been the most difficult interaction to model for orbital (PNO)-based methods historically. To highlight the capabilities of our code, we tested it on large water clusters, as well as insulin (787 atoms). [ABSTRACT FROM AUTHOR]

A consideration of the point group symmetry of molecules is often advantageous from a computational efficiency standpoint and sometimes necessary for the correct treatment of chemical physics problems. Many modern electronic structure software packages include a treatment of symmetry, but these are sometimes incomplete or unusable outside of that program's environment. Therefore, we have developed the MolSym package for handling molecular symmetry and its associated functionalities to provide a platform for including symmetry in the implementation and development of other methods. Features include point group detection, molecule symmetrization, arbitrary generation of symmetry element sets and character tables, and symmetry adapted linear combinations of real spherical harmonic basis functions, Cartesian displacement coordinates, and internal coordinates. We present some of the advantages of using molecular symmetry as achieved by MolSym, particularly with respect to Hartree–Fock theory, and the reduction of finite difference displacements in gradient/Hessian computations. This package is designed to be easily integrated into other software development efforts and may be extended to further symmetry applications. [ABSTRACT FROM AUTHOR]

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br–O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies −0.2 kcal mol−1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol−1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol−1 barrier. The final product energies lie −37.2, −67.9, and −7.3 kcal mol−1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin–orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol−1. [ABSTRACT FROM AUTHOR]

Multifidelity modeling is a technique for fusing the information from two or more datasets into one model. It is particularly advantageous when one dataset contains few accurate results and the other contains many less accurate results. Within the context of modeling potential energy surfaces, the low-fidelity dataset can be made up of a large number of inexpensive energy computations that provide adequate coverage of the N-dimensional space spanned by the molecular internal coordinates. The high-fidelity dataset can provide fewer but more accurate electronic energies for the molecule in question. Here, we compare the performance of several neural network-based approaches to multifidelity modeling. We show that the four methods (dual, Δ-learning, weight transfer, and Meng–Karniadakis neural networks) outperform a traditional implementation of a neural network, given the same amount of training data. We also show that the Δ-learning approach is the most practical and tends to provide the most accurate model. [ABSTRACT FROM AUTHOR]

Lead is heavily utilised in industry, but it bears potential harm due to its toxic nature. Identifying lead is important in diagnosing and preventing lead poisoning. Hydride generation via inductively-coupled plasma mass spectrometer (HG-ICP-MS) can detect small amounts of lead with high accuracy. Here we investigate seven structures of Pb $ _2 $ 2 H $ _2^+ $ 2 + . The geometries of each structure were optimised using coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] method with the cc-pwCVQZ-PP basis set for lead and the cc-pVQZ basis set for hydrogen. Likewise, harmonic vibrational frequencies were obtained using the same method. Anharmonic frequencies were determined using CCSD(T) with cc-pwCVTZ-PP and cc-pVTZ basis sets for lead and hydrogen, respectively. The focal point approach (FPA), including CCSDT and CCSDT(Q) computations, was applied to the seven structures to find the relative energies with respect to the global minimum butterfly structure. Relative to the global minimum butterfly structure, the energies for the planar transition state, monobridged, vinylidene-like, trans, cis, and linear structures are: 1.1, 18.3, 26.5, 27.8, 28.8, and 93.7 kcal mol $ ^{-1} $ − 1 , respectively. Natural bond orders and dipole moments are reported as well. The predicted vibrational frequencies should be helpful in guiding the laboratory observation of Pb $ _2 $ 2 H $ _2^+ $ 2 + . [ABSTRACT FROM AUTHOR]

Criegee intermediates, formed from the ozonolysis of alkenes, are known to have a role in atmospheric chemistry, including the modulation of the oxidizing capacity of the troposphere. Although studies have been conducted since their discovery, the synthesis of these species in the laboratory has ushered in a new wave of investigations of these structures, both theoretically and experimentally. In some of these theoretical studies, high-order corrections for correlation energy are included to account for the mid multi-reference character found in these systems. Many of these studies include a focus on kinetics; therefore, the calculated energies should be accurate (<1 kcal/mol in error). In this research, we compute the enthalpies of formation for a small set of Criegee intermediates, including higher-order coupled cluster corrections for correlation energy up to coupled cluster with perturbative quintuple excitations. The enthalpies of formation for formaldehyde oxide, anti-acetaldehyde oxide, syn-acetaldehyde oxide, and acetone oxide are presented at 0 K as 26.5, 15.6, 12.2, and 0.1 kcal mol−1, respectively. Additionally, we do not recommend the coupled cluster with perturbative quadruple excitations [CCSDT(Q)] energy correction, as it is approximately twice as large as that of the coupled cluster with full quadruple excitations (CCSDTQ). Half of the CCSDT(Q) energy correction may be included as a reliable, cost-effective estimation of CCSDTQ energies for Criegee intermediates. [ABSTRACT FROM AUTHOR]

While the dithiolene-based N-heterocyclic silane (4) reacts with two equivalents of BX3 (X = Br, I) to give zwitterionic Lewis adducts 5 and 8, respectively, the parallel reaction of 4 with BCl3 results in 10, a dithiolene-substituted N-heterocyclic silane, via the Si-S bond cleavage. Unlike 5, the labile 8 may be readily converted to 9 via BI3-mediated cleavage of the Si-N bond. The formation of 5 and 8 confirms that 4 uniquely possesses dual nucleophilic sites: (a) the terminal sulphur atom of the dithiolene moiety; and (b) the backbone carbon of the N-heterocyclic silane unit. [ABSTRACT FROM AUTHOR]

In 2021 Olaru, Mebs, and Beckmann reported the synthesis of remarkable cationic carbene analogues $ \mathrm{PR_2^+} $ P R 2 + and $ \mathrm{AsR_2^+} $ As R 2 + . This work followed the same group's synthesis of $ \mathrm{SbR_2^+} $ Sb R 2 + and $ \mathrm{BiR_2^+} $ Bi R 2 + . To better understand these important systems, $ \mathrm{SbH_2^+} $ Sb H 2 + and $ \mathrm{BiH_2^+} $ Bi H 2 + have been studied with high level ab initio quantum mechanical methods. Geometries were optimised with the CCSDT(Q) method with the cc-pwCVTZ-PP basis set using small core pseudopotentials. Fundamental vibrational frequencies were computed to provide theoretical predictions for future synthetic studies. Relative energies with respect to $ \mathrm{Pn^+} $ P n + + $ \mathrm{H_2} $ H 2 ( $ \mathrm{Pn} $ Pn = $ \mathrm{Sb} $ Sb , $ \mathrm{Bi} $ Bi ) were determined at the CCSDT/CBS level of theory via the Focal Point Analysis method, and anharmonic zero-point vibrational energy and higher order contributions were also computed. For $ \mathrm{SbH_2^+} $ Sb H 2 + , we obtained R = 1.697 Å and $ \theta =90.8^{\circ } $ θ = 90.8 ∘ with $ \mathrm{Sb^+ + H_2 \rightarrow SbH_2^+} $ S b + + H 2 → Sb H 2 + reaction enthalpy of $ \Delta H = -18.71\,{\rm kcal}\ {\rm mol}^{-1} $ ΔH = − 18.71 kcal mol − 1 . For $ \mathrm{BiH_2^+} $ Bi H 2 + , the analogous results were R = 1.774 Å, $ \theta =89.7^{\circ } $ θ = 89.7 ∘ , and $ \Delta H = -7.64\,{\rm kcal}\,{\rm mol}^{-1} $ ΔH = − 7.64 kcal mol − 1 for $ \mathrm{Bi^+ + H_2 \rightarrow BiH_2^+} $ B i + + H 2 → Bi H 2 + . [ABSTRACT FROM AUTHOR]

Butyl radicals (n-, s-, i-, and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo-tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700 to 3125 cm−1, allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800–3000 cm−1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra. [ABSTRACT FROM AUTHOR]

Small carbon clusters (Cn, n = 2-15) are produced in a molecular beam by pulsed laser vaporization and studied with vacuum ultraviolet (VUV) photoionization mass spectrometry. The required VUV radiation in the 8-12 eV range is provided by the Advanced Light Source (ALS) at the Lawrence Berkeley National Laboratory. Mass spectra at various ionization energies reveal the qualitative relative abundances of the neutral carbon clusters produced. By far the most abundant species is C3. Using the tunability of the ALS, ionization threshold spectra are recorded for the clusters up to 15 atoms in size. The ionization thresholds are compared to those measured previously with charge-transfer bracketing methods. To interpret the ionization thresholds for different cluster sizes, new ab initio calculations are carried out on the clusters for n = 4-10. Geometric structures are optimized at the CCSD(T) level with cc-pVTZ (or cc-pVDZ) basis sets, and focal point extrapolations are applied to both neutral and cation species to determine adiabatic and vertical ionization potentials. The comparison of computed and measured ionization potentials makes it possible to investigate the isomeric structures of the neutral clusters produced in this experiment. The measurements are inconclusive for the n = 4-6 species because of unquenched excited electronic states. However, the data provide evidence for the prominence of linear structures for the n = 7, 9, 11, 13 species and the presence of cyclic C10.

The acetonyl radical (•CH2COCH3) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. No submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol−1 higher in energy than acetonyl radical + O2. The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O2. Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures. [ABSTRACT FROM AUTHOR]

Sulfurous acid (H2SO3) is an infamously elusive molecule. Although some theoretical papers have supposed possible roles for it in more complicated systems, it has yet to be experimentally observed. To aid experiment in detecting this molecule, we have examined the H2O + SO2 potential energy surface at the CCSDT(Q)/CBS//CCSD(T)-F12b/cc-pVTZ-F12b level of theory to resolve standing discrepancies in previous reports and predict the gas-phase vibrational spectrum for H2SO3. We find that sulfurous acid has two potentially detectable rotamers, separated by 1.1 kcal mol−1 ΔH0K with a torsional barrier of 1.6 kcal mol−1. The sulfonic acid isomer is only 6.9 kcal mol−1 above the lowest enthalpy sulfurous acid rotamer, but the barrier to form it is 57.2 kcal mol−1. Error in previous reports can be attributed to misidentified stationary points, the use of density functionals that perform poorly for this system, and, most importantly, the basis set sensitivity of sulfur. Using VPT2+K, we determine that the intense S=O stretch fundamental of each species is separated from other intense peaks by at least 25 cm−1, providing a target for identification by infrared spectroscopy. [ABSTRACT FROM AUTHOR]

The reaction between two OH radicals has been studied in the context of combustion chemistry, atmospheric chemistry, and astrochemistry. The global minimum on the triplet reaction surface is the van der Waals complex between water and O( 3 P) atom. We optimised the geometry of the system using ab initio methods up to CCSDT/aug-cc-pVQZ. We use the focal point approach to determine an accurate dissociation energy of 1.35 kcal mol − 1 for the complex dissociating into water and O( 3 P). We also predict the anharmonic vibrational frequencies of the complex in anticipation of future experiments aimed at detecting this system. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Vinoxy radicals are involved in numerous atmospheric and combustion mechanisms. High-level theoretical methods have recently shed new light on the reaction of the unsubstituted vinoxy radical with O2. The reactions of 1-methylvinoxy radical and 2-methylvinoxy radical with molecular oxygen have experimental high pressure limiting rate constants, k∞, 5–7 times higher than that of the vinoxy plus O2 reaction. In this work, high-level ab initio quantum chemical computations are applied to the 2-methylvinoxy radical plus O2 system, namely, the formation and isomerization of the 1-oxo-2-propylperoxy radical, the immediate product of O2 addition to the 2-methylvinoxy radical. Multireference methods were applied to the entrance channel. No barrier to O2 addition could be located, and more sophisticated treatment of dynamic electron correlation shows that the principal difference between O2 addition to the vinoxy and 2-methylvinoxy radicals is a larger steric factor for 2-methylvinoxy + O2. This is attributed to the favorable interaction between the incoming O2 molecule and the methyl group of the 2-methylvinoxy radical. Via the focal point approach, energetics for this reaction were determined, in most cases, to chemical accuracy. The coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] correlation energy and Hartree–Fock energies were independently extrapolated to the complete basis set limit. A correction for the effect of higher excitations was computed at the CCSDT(Q)/6-31G level. Corrections for the frozen-core approximation, the Born–Oppenheimer approximation, the nonrelativistic approximation, and the zero-point vibrational energy were included. From the 1-oxo-2-propylperoxy radical, dissociation to reactants is competitive with the lowest energy isomerization pathway. The lowest energy isomerization pathway ultimately forms acetaldehyde, CO, and ·OH as the final products. [ABSTRACT FROM AUTHOR]

Benzene and 23 monosubstituted and 32 disubstituted derivatives of benzene were optimized for minimum energy structures using the B3LYP/cc-pVTZ method. The force fields of all the compounds were evaluated at their optimized geometries using the same method and basis set. In order to understand the effect of substitution(s) on the aromaticity of benzene, the aromaticity index based on interaction coordinates (AIBIC) values were computed for each and the change from the benzene value was obtained. This difference, the substituent effect based on interaction coordinates (SEBIC), quantifies the effect of the substituent on the aromaticity of benzene ring satisfactorily. It is found that the AIBIC of disubstituted benzenes (XC6H4Y) could be predicted well by adding the respective SEBIC(C6H5X) and SEBIC(C6H5Y) values to the AIBIC of benzene. The projected force fields of the meta and para fragments of the monosubstituted benzenes when chosen properly contain the information about the directing influence of the substituent in terms of the electron density based on interaction coordinates (EDBIC). When the EDBIC(para) > EDBIC(meta) relative to benzene, the substituent is ortho-para directing, while when the reverse is true, it is meta directing. The effect of conformational changes on aromaticity has been studied using aminophenols and dihydroxybenzenes. The additivity rule and the EDBIC concept work adequately well in that the methods can have several useful practical applications that will benefit various areas of science. A good understanding of the substituent effects and the ability to predict them should add a new dimension to the applications of AIBIC. [ABSTRACT FROM AUTHOR]

The peculiar electronic absorption spectrum of H2CN has been of great interest to experiment. Herein, this system is studied extensively by applying theoretical methods to the ground and low-lying excited electronic states. Employing a large breadth of high-level ab initio computations, including coupled cluster [CCSD(T) and CCSDT(Q)] and multireference configuration interaction [MRCISD+Q] methods, we comprehensively demonstrate that the most recent experimental and theoretical interpretations of the electronic spectrum of H2CN are in error. The previous assignments of the two broad features in the spectrum as the origin 000 (~35 050 cm-1) and 4²0 (~35 600 cm-1) * 4²0 * ²B2 transitions are both found to be incorrect. The presently reported transition energies suggest that the higher energy band near 35 600 cm-1 is the true origin band. Additionally, from the computed anharmonic vibrational frequencies of the * ²B2 and * ²A1 states, we show that this ~550 cm-1 band spacing cannot be attributed to a simple vibronic transition, as claimed by the 4²0 assignment. Possible alternative explanations for the appearance of the lower intensity band near 35 050 cm-1 are discussed. [ABSTRACT FROM AUTHOR]

We use a direct product basis, basis vectors computed by evaluating matrix-vector products, and rank reduction to calculate vibrational energy levels of uracil and naphthalene, with 12 and 18 atoms, respectively. A matrix representing the Hamiltonian in the direct product basis and vectors with as many components as there are direct product basis functions are neither calculated nor stored.We also introduce an improvement of the Hierarchical Intertwined Reduced-Rank Block Power Method (HIRRBPM), proposed previously in Thomas and Carrington, Jr. [J. Chem. Phys. 146, 204110 (2017)]. It decreases the memory cost of the HI-RRBPM and enables one to compute vibrational spectra of molecules with over a dozen atoms with a typical desktop computer. [ABSTRACT FROM AUTHOR]