CCSD Method Research Articles

A kinetic and mechanistic study of the self-reaction between two propargyl radicals.

The propargyl radical plays a critical role in various chemical processes, including hydrocarbon combustion, flame synthesis, and interstellar chemistry. Its unique stability arises from the delocalization of π-electrons, allowing it to participate in a wide range of reactions despite being a radical. The self-reaction of propargyl radicals is a fundamental step in synthesizing polycyclic aromatic hydrocarbons. In this work, therefore, a computational study into the C3H3 + C3H3 potential energy surface has been carefully characterized. The calculated results indicate that the reaction can occur by H-abstraction or addition of two propargyl radicals together. The H-abstraction mechanism can create the products P3 (H2CCC + H3CCCH) and P4 (H2CCCH2 + HCCCH) but the energy barriers of the two H-abstraction channels are very high (from 12 to 22kcal/mol). In contrast, the addition mechanism of two propargyl radicals forming the intermediates (I1, I5, I12) and the bimolecular products (P1, P2, P7, P11, P12) are dominant. Among the bimolecular products, the P11 (C6H4 + H2) product is the most energetically favorable, and the channel leading to this product is also the most preferred path compared to all other paths throughout the PES. The calculated enthalpy changes of various reaction paths in this study are in good agreement with the available literature data, indicating that the CCSD(T) method is suitable for the title reaction. The overall rate constant of the reaction depends on both temperature and pressure, reducing with temperature but rising with pressure. The calculated results agree closely with the available experimental values ​​and previous calculated data. Thus, it can be affirmed that in addition to the CASPT2 method as applied in the study of Georgievskii et al. (Phys. Chem. Chem. Phys., 2007, 9, 4259-4268), the CCSD(T) method is also very good for the self-reaction of two propargyl radicals. The M06-2X and CCSD(T) methods with the aug-cc-pVTZ basis set were used to optimize and calculate single-point energies for all species of the reaction. The bimolecular rate constants of the dominant reaction paths were predicted in the temperature and pressure ranges of 300-1800K and 0 - 76,000Torr, respectively, using the VTST and RRKM models with Eckart tunneling correction for the H-shift steps.

Extending the Perturbative Triples Correction in Coupled-Cluster Theory to Account for Inactive Orbitals.

Coupled-cluster singles and doubles calculations with perturbative triples [CCSD(T)] offer high accuracy with significant cost. The standard way to reduce the cost of a CCSD(T) calculation is by excluding molecular orbitals from the correlated calculation. While this speeds up the calculation, it also throws away the contribution from the inactive orbitals to the correlation energy. Here, we propose extending the standard CCSD(T) method to account for the effects of these inactive orbitals. This approach is based on a perturbation expansion of the similarity-transformed Hamiltonian, and the final method includes external singles and doubles corrections along with a semi-internal triples term. Compared to all-electron CCSD(T) calculations with the cc-pCVTZ basis set for a set of small molecules, we recovered, on average, 98% of the total correlation energy while using only 30% of the molecular orbitals. Using 72% of the molecular orbitals gave 99.5% of the correlation energy.

Cooperativity and halonium transfer in the ternary NCI···CH3I···-CN halogen-bonded complex: An ab initio gas phase study.

The strength and nature of the two halogen bonds in the NCI···CH3I···-CN halogen-bonded ternary complex are studied in the gas phase via ab initio calculations. Different indicators of halogen bond strength were employed to examine the interactions including geometries, complexation energies, Natural Bond Order (NBO) Wiberg bond indices, and Atoms in Molecules (AIM)-based charge density topological properties. The results show that the halogen bond is strong and partly covalent in nature when CH3I donates the halogen bond, but weak and noncovalent in nature when CH3I accepts the halogen bond. Significant halogen bond cooperativity emerges in the ternary complex relative to the corresponding heterodimer complexes, NCI···CH3I and CH3I···-CN, respectively. For example, the CCSD(T) complexation energy of the ternary complex (-18.27kcal/mol) is about twice the sum of the complexation energies of the component dimers (-9.54kcal/mol). The halonium transfer reaction that converts the ternary complex into an equivalent one was also investigated. The electronic barrier for the halonium transfer was calculated to be 6.70kcal/mol at the CCSD(T) level. Although the MP2 level underestimates and the MP3 overestimates the barrier, their calculated MP2.5 average barrier (6.44kcal/mol) is close to that of the more robust CCSD(T) level. Insights on the halonium ion transfer reaction was obtained by examining the reaction energy and force profiles along the intrinsic reaction coordinate, IRC. The corresponding evolution of other properties such as bond lengths, Wiberg bond indices, and Mulliken charges provides specific insight on the extent of structural rearrangements and electronic redistribution throughout the entire IRC space. The MP2 method was used for geometry optimizations. Energy calculations were performed using the CCSD(T) method. The aug-cc-pVTZ basis set was employed for all atoms other than iodine for which the aug-cc-pVTZ-PP basis set was used instead.

A Critical Evaluation of the Hybrid KS DFT Functionals Based on the KS Exchange-Correlation Potentials.

We have developed a critical methodology for the evaluation of the quality of hybrid exchange-correlation (XC) density functional approximations (DFAs) based on very fundamental quantities, i.e., Kohn-Sham (KS) XC potentials, self-consistent electron densities, first ionization potentials (IPs), and total energies. Since the XC potentials, the primary objects in the current study, are not directly accessible for the hybrids, we calculate them by inverting the KS electron densities. Utilizing this methodology, we tested 155 hybrid DFAs available in the LIBXC library using FCI and CCSD(T) methods as a reference. We have found that a group of functionals produces very decent XC potentials, mainly those with a large mixture of Hartree-Fock exchange. Moreover, the value of IP strongly depends on the XC potential quality. On the other hand, we show that the XC energy is dominated by functional-driven error, which in some cases leads to substantial errors in electronic densities. The study shows new directions for constructing more accurate XC functionals within the KS-DFT framework.

Complete Active Space Iterative Coupled Cluster Theory.

In this work, we investigate the possibility of improving multireference-driven coupled cluster (CC) approaches with an algorithm that iteratively combines complete active space (CAS) calculations with tailored CC and externally corrected CC. This is accomplished by establishing a feedback loop between the CC and CAS parts of a calculation through a similarity transformation of the Hamiltonian with those CC amplitudes that are not encompassed by the active space. We denote this approach as the complete active space iterative coupled cluster (CASiCC) ansatz. We investigate its efficiency and accuracy in the singles and doubles approximation by studying the prototypical molecules H4, H8, H2O, and N2. Our results demonstrate that CASiCC systematically improves on the single-reference CCSD and the externally corrected CCSD methods across entire potential energy curves while retaining modest computational costs. However, the tailored coupled cluster method shows superior performance in the strong correlation regime, suggesting that its accuracy is based on error compensation. We find that the iterative versions of externally corrected and tailored coupled cluster methods converge to the same results.

Electric Field Dependence of EPR Hyperfine Coupling Constants.

In this article, we introduce a simple but reliable method to calculate the electric field dependence of the isotropic hyperfine coupling tensor Aiso for free radicals. This dependence, also referred to as the Bloembergen effect, can be of interest in analyzing EPR experiments for solid-state materials but is rarely studied for isolated radicals in the gaseous phase. The proposed method uses the numerical differentiation of the field-perturbed A tensor and, consequently, as a purely numerical method, does not depend on a quantum chemical method used to determine the hyperfine tensor A. To test the performance of the proposed method, we used a set of 28 systems, including seven organic radicals, the majority of which were taken from a recent benchmark study by Bartlett's group. We employed the well-tested and robust density functional theory (DFT) functional, namely CAM-B3LYP, and the variant of the CCSD method based on local pair natural orbitals, namely DLPNO-CCSD to calculate the Aiso tensor itself, and all first derivatives or Bloembergen effect constants ai(1).

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.

Accurate and efficient open-source implementation of domain-based local pair natural orbital (DLPNO) coupled-cluster theory using a t1-transformed Hamiltonian.

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.1kcal 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).

Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry.

Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)Λ contributions to reaction energies as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Λ/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol-1 and an L2D of 0.52 kcal mol-1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably "ballpark" how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.

Energetics of the OH radical H-abstraction reactions from simple aldehydes and their geminal diol forms.

Carbonyl compounds, especially aldehydes, emitted to the atmosphere, may suffer hydration in aerosols or water droplets in clouds. At the same time, they can react with hydroxyl radicals which may add or abstract hydrogen atoms from these species. The interplay between hydration and hydrogen abstraction is studied using density functional and quantum composite theoretical methods, both in the gas phase and in simulated bulk water. The H-abstraction from the aldehydic and geminal diol forms of formaldehyde, acetaldehyde, glycolaldehyde, glyoxal, methylglyoxal, and acrolein is studied to determine whether the substituent has any noticeable effect in the preference for the abstraction of one form or another. It is found that abstraction of the H-atom adjacent to the carbonyl group gives a more stable radical than same abstraction from the geminal diol in the case of formaldehyde, acetaldehyde, and glycolaldehyde. The presence of a delocalizing group in the Cα (a carbonyl group in glyoxal and methylglyoxal, and a vinyl group in acrolein), reverts this trend, and now the abstraction of the H-atom from the geminal diol gives more stable radicals. A further study was conducted abstracting hydrogen atoms from the other different positions in the species considered, both in the aldehydic and geminal diol forms. Only in the case of glycolaldehyde, the radical formed by H-abstraction from the -CH2OH group is more stable than any of the other radical species. Abstraction of the hydrogen atom in one of the hydroxyl groups in the geminal diol is equivalent to the addition of the •OH radical to the aldehyde. It leads, in some cases, to decomposition into a smaller radical and a neutral molecule. In these cases, some interesting theoretical differences are observed between the results in gas phase and (simulated) bulk solvent, as well as with respect to the method of calculation chosen. DFT (M06-2X, B2PLYP, PW6B95), CCSD(T), and composite (CBS-QB3, jun-ChS, SCVECV-f12) methods using Dunning basis sets and extrapolation to the CBS limit were used to study the energetics of closed shell aldehydes in their keto and geminal-diol forms, as well as the radical derived from them by hydrogen abstraction. Both gas phase and simulated bulk solvent calculations were performed, in the last case using the Polarizable Continuum Model.

Quantum chemical calculations of electron affinities of alkaline earth metal atoms (Ca, Sr, Ba, and Ra)

We performed high-level ab initio quantum chemical calculations, incorporating higher-order excitations, spin–orbit coupling (SOC), and the Gaunt interaction, to calculate the electron affinities (EAs) of alkaline earth (AE) metal atoms (Ca, Sr, Ba, and Ra), which are notably small. The coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method is insufficient to accurately calculate the EAs of AE metal atoms. Higher-order excitations proved crucial, with the coupled-cluster singles, doubles, and triples with perturbative quadruples [CCSDT(2)Q] method effectively capturing dynamic electron correlation effects. The contributions of SOC (ΔESOs) to the EAs calculated using the multireference configuration interaction method with the Davidson correction, including SOC, positively enhance the EAs; however, these contributions are overestimated. The Dirac–Hartree–Fock (DHF)-CCSD(T) method addresses this overestimation and provides reasonable values for ΔESO (ΔESO−D). Employing additional sets of diffuse and core–valence correlation basis sets is critical for accurately calculating the EAs of AE metal atoms. The contributions of the Gaunt interaction (ΔEGaunt) to the EAs of AE metal atoms are negligible. Notably, the CCSDT(2)Q with the complete basis set limit + ΔESO−D + ΔEGaunt produced EA values for Ca, Sr, and Ba that closely aligned with experimental data and achieved accuracy exceeding the chemical accuracy. Based on our findings, the accurately proposed EA for Ra is 9.88 kJ/mol.

A new "gold standard": Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing.

A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the "gold standard" CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.

State-Specific Coupled-Cluster Methods for Excited States.

We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.

Kinetic Computation of Cyclization Reactions of Large Keto-Hydroperoxide Radicals in Low Temperature Combustion.

The cyclization reactions of keto-hydroperoxide (KHP) radicals leading to the formation of keto cyclic ethers and OH radicals play an important role in low temperature combustion for hydrocarbon fuels or oxygenated hydrocarbon fuels. However, due to the lack of kinetic data of cyclization reactions of KHP radicals, researchers often derive high-pressure-limit rate constants of cyclization reactions of KHP radicals from analogous cyclization reactions of hydroperoxyl alkyl radicals during construction of the combustion mechanism. This study aims to systematically investigate the kinetics of cyclization reactions of KHP radicals involving short-to-large-sized radicals. The studied reactions are divided into 7 reaction classes, according to the size of the cyclic transition state, the conjugative effect (whether KHP radicals are resonance-stabilized or not), and the position of the carbonyl group (whether the carbonyl group is inside or outside of the reaction center). The isodesmic reaction method, in conjunction with transition state theory, is utilized for each reaction class to compute the energy barriers and high-pressure-limit rate constants at the DFT level. The study revealed that energy barriers calculated at the DFT level with correction by the isodesmic reaction method are close to the results from the benchmark CCSD(T) method. To develop more accurate rate rules, these reaction classes are further divided into subclasses based on the relative site of the OOH group with the carbonyl group, the type of carbon atoms where the OOH group is located, and the type of carbon atoms where the radical site is located. For each subclass, high-pressure-limit rate rules are derived by averaging the rate constants of reactions in the subclass, and it is found that the maximum absolute deviation of the energy barrier and the ratio of the largest rate constant to the smallest rate constant among reactions in each subclass are within chemical accuracy limits, indicating acceptable use of the developed rate rules. A comparison of the rate constants for cyclization reactions of KHP radicals with the values of analogous cyclization reactions of hydroperoxyalkyl radicals as provided in reported mechanisms is made. Additionally, a comparison is drawn between our developed rate rules for subclasses of the cyclization reactions of KHP radicals and the rate rules for analogous subclasses of cyclization reactions of hydroperoxyl alkyl radicals. These comparisons demonstrate significant differences and highlight the necessity for improved rate rules for cyclization reactions of KHP radicals to enhance the automatically generated combustion mechanisms for hydrocarbon and oxygenated hydrocarbon fuels.

Thermal rate constants and kinetic isotope effects of the H + H2O2 reactions: barrier height and reaction energy from single- and multireference methods.

One of the more significant sub-mechanisms of H2/O2 combustion involves the reaction of hydrogen peroxide with hydrogen atoms (H + H2O2), resulting in the production of OH + H2O (R1) and H2 + HO2 (R2) paths. Previous experimental and ab initio calculations reveal some variations in the barrier height for (R1). To improve the energetics of both (R1) and (R2), single reference and multireference ab initio methods are employed, and the rate constants and H/D kinetic isotope effects (KIEs) are calculated as a function of temperature. For (R1), the best results for the barrier height and reaction energies computed with the CASPT2(15,11)/aug-cc-pV6Z are 5.2 and - 70.3kcal.mol-1, respectively. CCSD(T)/aug-cc-pV5Z + CV (core-valence) calculations for (R2) give 9.7 and - 15.6kcal.mol-1 to those parameters. The CVT/SCT rate constants of both paths agree well with the fitted rate constants from uncertainty-weighted statistical analysis of the 14-mechanism of H2/O2. The kinetic isotopic effect (kH/kD) for the reaction D + H2O2 → DH + HO2 was found to be 0.47, which is in excellent agreement with the experimental value of 0.43. The structures of reactants, transition state, and products of (R1) and (R2) are calculated with the aug-cc-pVTZ basis set and M062X DFT, CCSD(T), and CASSCF methods. The barrier heights and reaction energies of (R1) and (R2) are computed using the M06-2X, CCSD(T), MRCI, and CASPT2 methods and various basis sets. The rate constants are calculated with the variational transition state theory including multidimensional tunneling corrections (VTST-MT), with potential energy surfaces built by the M06-2X/aug-cc-pVTZ approach.

Structure, dissociation pathways and thermochemistry of some nitrogen-containing astrophysical molecules: a theoretical study

The structure, relative energies, dissociation pathways and thermochemistry of acetonitrile (CH3CN or methyl cyanide), acrylonitrile (CH2CHCN or vinyl cyanide), propionitrile (CH3CH2CN or ethyl cyanide) and their stable isomers have been studied in detail using density functional B3LYP, B3PW91 and ab initio MP2, QCISD and CCSD(T) methods. Several space missions have already confirmed the abundance of these molecules in different interstellar mediums (ISMs). The existence of related radicals, atoms and molecules are also confirmed using various IR-spectrometric techniques and that may be evolved in solar irradiation in the cold dark core TMC-1 or massive star-forming regions such as Sagittarius B2 or Oion-KL. Possible ways of fragmentation of ground state CH3NC, CH2CHNC and CH3CH2NC through different dissociation channels and their link to other isomers have been considered to justify the presence of such observed radicals in the interstellar cloud. Cyanide(-C≡N), iso-cyanide (-N = C) and imine (-C = N-H) isomeric forms have been considered here and out of these, cyanide isomers are more stable species compared to its iso-nitrile and imine form. Favourable product molecules, such as CH2 radical, C2H2, C2H4 and HCN/HNC, are considered in the dissociation channels.

Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes.

Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.

A Simulation of Hydrazine Molecule’s Potential Energy Surface using Variational Quantum Eigensolver Algorithm

Quantum computing is a technology that utilizes the principles of quantum mechanics to perform complex computational processes. In this work, we use Qiskit Module from IBM to do quantum computational calculation using Variational Quantum Eigensolver (VQE) algorithm. VQE is a hybrid quantum-classical method that combines a quantum computer to measure energies and a classical computer to process the measurement results and update the parameters of the quantum computer. The purpose of VQE is to find the ground state energy of a chemical system. In the previous study, many of the VQE calculations have been done on simple molecules. So, in this study, we would like to use Hydrazine molecule as our object of VQE calculation. Furthermore, these results will be compared with the results from the classical calculation (MP2, CCSD(T), QCISD(T), and CASSCF) methods for testing the effectiveness of VQE using Unitary Coupled-Cluster Single and Double excitations (UCCSD) Ansatz. The quantum algorithm based on the UCCSD Ansatz led to a simplification of the algorithm by reducing the circuit depth. Then, the possibility to use active space approximation, can be used to reduce the quantum gates while trying to keep a good level of accuracy. In this study, we chose (2,2) and (4,4) active spaces. Based on the results, as we increase the size of the active space during the evaluation of the single-point energy, the estimated ground states obtained from the VQE algorithm yield nearly identical values. Conversely, in CASSCF calculations, expanding the active space introduces more energy corrections, thus making it more sensitive. Additionally, when examining potential energy surfaces, VQE demonstrates results that gradually align with CCSD(T) and QCISD(T) methods.

Benchmarking model chemistry composite calculations for vertical ionization potential of molecular systems

In chemical science, the vertical ionization potential (VIP) is a crucial metric for understanding the electronegativity, hardness and softness of chemical material systems as well as the electronic structure and stability of molecules. Ever since the last century, the model chemistry composite methods have witnessed tremendous developments in computing the thermodynamic properties as well as the barrier heights. However, their performance in realm of the vertical electron processes of molecular systems has been rarely explored. In this study, we for the first time benchmarked the model chemistry composite methods (e.g., CBS-QB3, G4 and W1BD) in comparison with the commonly used Koopmans's theorem (KT), electron propagator theory (e.g., OVGF, D2, P3 and P3+) and CCSD(T) methods in calculating the VIP for up to 613 molecular systems with available experimental measurements. The large-scale test calculations strongly showed that the CBS-QB3 model chemistry composite technique can be well recommended to calculate VIP from the perspectives of accuracy, economy and applicability. Notably, the VIP values of up to 7 molecules were identified to have the absolute errors of larger than 0.3 eV at all calculation levels, which have strong hints that their VIP experimental values should be re-investigated.

Review and Theoretical Analysis of Fluorinated Radicals in Direct CAr−H Functionalization of (Hetero)arenes

AbstractWe highlight key contributions in the field of direct radical CAr− H (hetero)aromatic functionalization involving fluorinated radicals. A compilation of Functional Group Transfer Reagents and their diverse activation mechanisms leading to the release of radicals are discussed. The substrate scope for each radical is analyzed and classified into three categories according to the electronic properties of the substrates. Density functional theory computational analysis provides insights into the chemical reactivity of several fluorinated radicals through their electrophilicity and nucleophilicity parameters. Theoretical analysis of their reduction potentials also highlights the remarkable correlation between electrophilicity and oxidizing ability. It is also established that highly fluorinated radicals (e.g. ⋅OCF3) are capable of engaging in single‐electron transfer (SET) processes rather than radical addition, which is in good agreement with experimental literature data. A reactivity scale, based on activation barrier of addition of these radicals to benzene is also elaborated using the high accuracy DLPNO‐(U)CCSD(T) method.