Accurate CCSD Research Articles

Accurate Prediction of Nuclear Magnetic Resonance Parameters via the XYG3 Type of Doubly Hybrid Density Functionals.

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful and versatile tools in elucidating molecular structures. To eliminate ambiguities of experimental assignments, accurate calculations of NMR spectra are of great importance. Here, a method for theoretical evaluation of the NMR shielding constants by analytic derivatives using gauge including atomic orbitals (GIAO) has been implemented for the XYG3 type of doubly hybrid density functionals (xDH), namely, the GIAO-xDH method. Benchmark calculations on shielding constants and chemical shifts demonstrate the remarkable accuracy of the GIAO-xDH method, compared to the accurate CCSD(T) references. It is shown here that the XYGJ-OS functional is able to give a mean absolute deviation (MAD) of ∼3.0 ppm in the calculated shielding constants for 13C, 15N, 17O, 19F, while both XYGJ-OS and xDH-PBE0 functionals are able to provide a satisfactory estimation of chemical shifts with MADs of ∼0.03 and 1.0 ppm for 1H and 13C, respectively. The basis set influence upon the method has been examined and a computational scheme considering both accuracy and efficiency has been proposed and tested to predict the experimental 13C chemical shifts of five medium-sized natural product molecules, yielding a MAD of ∼1.0 ppm, which demonstrates the practical feasibility of the GIAO-xDH method.

Non-covalent interactions atlas benchmark data sets 4: σ-hole interactions.

The SH250×10 dataset presented here extends the Non-Covalent Interactions Atlas database (https://www.nciatlas.org) to complexes bound by σ-hole interactions - halogen, chalcogen and pnictogen bonds. It comprises 250 complexes where Cl, Br, I, S, Se, P and As interact with diverse electron donors. An accurate CCSD(T)/CBS benchmark is provided for ten points along a dissociation curve of each complex. The SH250×10 set is used in testing a wide variety of DFT functionals and semiempirical quantum-mechanical methods. In DFT calculations, the new data set exposes large errors of some functionals related to exaggerated charge transfer. The size and diversity of the data set have also been exploited in the reparametrization of a halogen-bond correction for the PM6 semiempirical method.

Vibrational Properties of CO Adsorbed on Au Single Atom Catalysts on TiO2(101), ZrO2(101), CeO2(111), and LaFeO3(001) Surfaces: A DFT Study

The nature and local environment of Au single atoms supported and stabilized on four different oxides is studied by means of DFT + U calculations using CO as probe molecule and its stretching frequency, ωe, as a fingerprint of the site where the Au atom is bound. Four oxides are considered, anatase TiO2, tetragonal ZrO2, cubic CeO2, and a perovskite LaFeO3. In this latter case a recently reported experimental study has detected a stretching mode for CO adsorbed on Au1/LaFeO3 of 2215 cm−1, with a large blue shift, ∆ω(CO) = 72 cm−1 with respect to free CO. In order to identify the Au adsorption site that can give rise to this large blue-shift we have considered five cases: (a) Au replacing a lattice cation, (Au)subM; (b) Au replacing a lattice O anion, (Au)subO; (c) Au adsorbed on the surface, (Au)ads; (d) Au bound to an extra O atom on the surface, (AuO)ads, or (e) Au bound to two extra O atoms on the surface, (AuO2)ads. It turns out that the correct reproduction of ∆ω for CO adsorbed on positively charged gold, Auδ+, is challenging for DFT. Therefore, we have performed a comparative study of Auδ+-CO molecular compounds for which ωe(CO) is known experimentally using various kinds of DFT functionals and accurate CCSD and CCSD(T) quantum chemistry methods. Also based on this comparison we propose a tentative assignment for the observed frequency of CO adsorbed on Au1/LaFeO3 single atom catalyst.Graphic

Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data.

Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH2, Ac-Ala-NH2, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.

High-Level Systematic Ab Initio Comparison of Carbon- and Silicon-Centered SN2 Reactions.

We characterize the stationary points along the Walden inversion, front-side attack, and double-inversion pathways of the X– + CH3Y and X– + SiH3Y [X, Y = F, Cl, Br, I] SN2 reactions using chemically accurate CCSD(T)-F12b/aug-cc-pVnZ [n = D, T, Q] levels of theory. At the carbon center, Walden inversion dominates and proceeds via prereaction (X–···H3CY) and postreaction (XCH3···Y–) ion-dipole wells separated by a usually submerged transition state (X–H3C–Y)−, front-side attack occurs over high barriers, double inversion is the lowest-energy retention pathway for X = F, and hydrogen- (F–···HCH2Y) and halogen-bonded (X–···YCH3) complexes exist in the entrance channel. At the silicon center, Walden inversion proceeds through a single minimum (X–SiH3–Y)−, the front-side attack is competitive via a usually submerged transition state separating pre- and postreaction minima having X–Si–Y angles close to 90°, double inversion occurs over positive, often high barriers, and hydrogen- and halogen-bonded complexes are not found. In addition to the SN2 channels (Y– + CH3X/SiH3X), we report reaction enthalpies for proton abstraction (HX + CH2Y–/SiH2Y–), hydride substitution (H– + CH2XY/SiH2XY), XH···Y– complex formation (XH···Y– + 1CH2/1SiH2), and halogen abstraction (XY + CH3–/SiH3– and XY– + CH3/SiH3).

Coupled Cluster Benchmark of New DFT and Local Correlation Methods: Mechanisms of Hydroarylation and Oxidative Coupling Catalyzed by Ru(II, III) Chloride Carbonyls.

We have evaluated a set of accurate canonical CCSD(T) energies for stationary points on the potential energy surface for Ru(II, III) chloride carbonyl catalysis of two competing reactions between benzene and methyl acrylate (MA), namely, hydroarylation and oxidative coupling. We have then applied this set to evaluate the performance of localized orbital coupled-cluster methods and several new and common density functionals. We find that (a) DLPNO-CCSD(T) with TightPNO cutoffs is an acceptable substitute for full canonical CCSD(T) calculations on this system; (b) for the closed-shell systems where it could be applied, LNO-CCSD(T) with tight convergence criteria is very close to the canonical results; (c) the recent ωB97X-V and ωB97M-V functionals exhibit superior performance to commonly used DFT functionals in both closed- and open-shell calculations; (d) the revDSD-PBEP86 revision of the DSD-PBEP86 double hybrid represents an improvement over the original, even though transition metals were not involved in its parametrization; and (e) DSD-SCAN and DOD-SCAN show comparable efficiency. Most tested (meta)-GGA and hybrid density functionals perform better for open-shell than for closed-shell complexes; this is not the case for the double hybrids considered.

A Spectroscopic Validation of the Improved Lennard-Jones Model.

The Lennard–Jones (LJ) and Improved Lennard–Jones (ILJ) potential models have been deeply tested on the most accurate CCSD(T)/CBS electronic energies calculated for some weakly bound prototype systems. These results are important to plan the correct application of such models to systems at increasing complexity. CCSD(T)/CBS ground state electronic energies were determined for 21 diatomic systems composed by the combination of the noble gas atoms. These potentials were employed to calculate the rovibrational spectroscopic constants, and the results show that for 20 of the 21 pairs the ILJ predictions agree more effectively with the experimental data than those of the LJ model. The CCSD(T)/CBS energies were also used to determine the parameter of the ILJ form, related to the softness/hardness of the interacting partners and controlling the shape of the potential well. This information supports the experimental finding that suggests the adoption of for most of the systems involving noble gas atoms. The He-Ne and He-Ar molecules have a lifetime of less than 1ps in the 200–500 K temperature range, indicating that they are not considered stable under thermal conditions of gaseous bulks. Furthermore, the controversy concerning the presence of a “virtual” or a “real” vibrational state in the He molecule is discussed.

Non-Covalent Interactions Atlas Benchmark Data Sets 3: Repulsive Contacts.

The new R739×5 data set from the Non-Covalent Interactions Atlas series (www.nciatlas.org) focuses on repulsive contacts in molecular complexes, covering organic molecules, sulfur, phosphorus, halogens, and noble gases. Information on the repulsive parts of the potential energy surface is crucial for the development of robust empirically parametrized computational methods. We use the new data set of highly accurate CCSD(T)/CBS interaction energies to test selected density functional theory (DFT) and semiempirical quantum-mechanical methods. The double-hybrid functionals were the best performing, with the revDSD-PBEP86-D3 being the most accurate DFT method, followed by the range-separated ωB97X functionals. Out of semiempirical methods, GFN2-xTB yielded the best results. On the example of the PM6 method, we analyze the source of error and its relation to the difficulties in the description of conformational energies, and we also devise an immediately applicable correction that fixes the most serious uncorrected issues previously encountered in practical calculations.

Accurate Reduced-Cost CCSD(T) Energies: Parallel Implementation, Benchmarks, and Large-Scale Applications.

The accurate and systematically improvable frozen natural orbital (FNO) and natural auxiliary function (NAF) cost-reducing approaches are combined with our recent coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] implementations. Both of the closed- and open-shell FNO-CCSD(T) codes benefit from OpenMP parallelism, completely or partially integral-direct density-fitting algorithms, checkpointing, and hand-optimized, memory- and operation count effective implementations exploiting all permutational symmetries. The closed-shell CCSD(T) code requires negligible disk I/O and network bandwidth, is MPI/OpenMP parallel, and exhibits outstanding peak performance utilization of 50–70% up to hundreds of cores. Conservative FNO and NAF truncation thresholds benchmarked for challenging reaction, atomization, and ionization energies of both closed- and open-shell species are shown to maintain 1 kJ/mol accuracy against canonical CCSD(T) for systems of 31–43 atoms even with large basis sets. The cost reduction of up to an order of magnitude achieved extends the reach of FNO-CCSD(T) to systems of 50–75 atoms (up to 2124 atomic orbitals) with triple- and quadruple-ζ basis sets, which is unprecedented without local approximations. Consequently, a considerably larger portion of the chemical compound space can now be covered by the practically “gold standard” quality FNO-CCSD(T) method using affordable resources and about a week of wall time. Large-scale applications are presented for organocatalytic and transition-metal reactions as well as noncovalent interactions. Possible applications for benchmarking local CCSD(T) methods, as well as for the accuracy assessment or parametrization of less complete models, for example, density functional approximations or machine learning potentials, are also outlined.

Can density functional theory ‘Cope’ with highly fluxional shapeshifting molecules?

Shape-shifting molecules such as bullvalene undergo rapid structural reorganizations via degenerate Cope rearrangements. Here, we obtain accurate CCSD(T)/CBS barrier heights and reaction energies for a wide range of Cope rearrangements in substituted bullvalenes (C10H9R, R = NH3, OH, CH3, H, F, Cl, SH, and CN). We use this benchmark dataset to evaluate the performance of DFT and ab initio methods for the kinetics and thermodynamics of these reactions. The reaction barrier heights pose a significant challenge for DFT methods – the best methods attain root-mean-square deviations of 4.9 (BMK), 4.5 (PBE0), 4.2 (PW6B95), and 3.8 (B1B95) kJ mol−1. Overall, only three DFT functionals (BMK, PW6B95, and MN12-SX) are able to surpass (or attain near) chemical accuracy for both barrier heights and reaction energies. In contrast, the double-hybrid DFT procedures ωB97X-2(LP), ωB97X-2(TQZ), PWPB95-D3, PBEQI-DH, and DSD-PBEB95-D3 give good-to-excellent performance.

Thermal Stability of Bis-Tetrazole and Bis-Triazole Derivatives with Long Catenated Nitrogen Chains: Quantitative Insights from High-Level Quantum Chemical Calculations.

Azobis tetrazole and triazole derivatives containing long catenated nitrogen atom chains are of great interest as promising green energetic materials. However, these compounds often exhibit poor thermal stability and high impact sensitivity. Kinetics and mechanism of the primary decomposition reactions are directly related to these issues. In the present work, with the aid of highly accurate CCSD(T)-F12 quantum chemical calculations, we obtained reliable bond dissociation energies and activation barriers of thermolysis reactions for a number of N-rich heterocycles. We studied all existing 1,1'-azobistetrazoles containing an N10 chain, their counterparts with the 5,5'-bridging pattern, and the species with hydrazo- and azoxy-bridges, which are often present energetic moieties. The N8-containing azobistriazole was considered as well. For all compounds studied, the radical decomposition channel was found to be kinetically unfavorable. All species decompose via the ring-opening reaction yielding a transient azide (or diazo) intermediate followed by the N2 elimination. In the case of azobistetrazole derivatives, the calculated effective activation barriers of decomposition are ∼26-33 kcal mol-1, which is notably lower than that of tetrazole (∼40 kcal mol-1). This fact agrees well with the low thermal stability and high impact sensitivities of the former species. The activation barriers of the N2 elimination were found to be almost the same for the azobis compounds and the parent tetrazole, and the effective decomposition barrier is determined by the thermodynamics of the tetrazole-azide rearrangement. In comparison with 1,1'-azobistetrazole, the hydrazo-bridged compound is more stable kinetically due to the lack of pi-conjugation in the azide intermediate. In turn, the azoxy-bridged compounds are entirely unstable due to tremendous azide stabilization by the hydrogen bond formation. In general, the 5,5'-bridged species are more thermally stable than their 1,1'-counterparts due to a much higher barrier of the N2 elimination. Apart from this, the highly accurate gas-phase formation enthalpies were calculated at the W1-F12 level of theory for all species studied.

Pressure DSC for energetic materials. Part 2. Switching between evaporation and thermal decomposition of 3,5-dinitropyrazole

Differential scanning calorimetry at an elevated external pressure (Pressure DSC) allows for shifting the vaporization of a sample to higher temperatures, thus often facilitating the direct observation of its thermal decomposition. In the present work, the thermolysis of a promising energetic material 3,5-dinitropyrazole was studied under pressures up to 10 MPa. The thermokinetic analysis of the datasets acquired at 2 and 5 MPa exhibited a pronounced kinetic compensation, thus allowing to build the joint formal kinetic model. The final kinetic scheme is comprised of the two parallel reactions, viz., the first-order process (Ea = 198 kJ mol−1) along with a first-order autocatalysis (Ea = 127 kJ mol−1). The experiment was complemented with the highly accurate CCSD(T)-F12 quantum chemical calculations. Theory revealed an unusual primary decomposition channel, viz., a sigmatropic [1,5]H-shift followed by the pyrazole ring opening yielding a molecular nitrogen and a nitro radical as simple primary products. Apart from this, the comparative thermogravimetry at a normal pressure yielded the vapor pressure of 3,5-dinitropyrazole along with the internally consistent set of phase change enthalpies. In general, the pressure DSC is a facile technique to study the true decomposition kinetics of the compounds that vaporize/sublime in conventional DSC experiments.

The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules

Maximum diversification of data is a central theme in building generalized and accurate machine learning (ML) models. In chemistry, ML has been used to develop models for predicting molecular properties, for example quantum mechanics (QM) calculated potential energy surfaces and atomic charge models. The ANI-1x and ANI-1ccx ML-based general-purpose potentials for organic molecules were developed through active learning; an automated data diversification process. Here, we describe the ANI-1x and ANI-1ccx data sets. To demonstrate data diversity, we visualize it with a dimensionality reduction scheme, and contrast against existing data sets. The ANI-1x data set contains multiple QM properties from 5 M density functional theory calculations, while the ANI-1ccx data set contains 500 k data points obtained with an accurate CCSD(T)/CBS extrapolation. Approximately 14 million CPU core-hours were expended to generate this data. Multiple QM calculated properties for the chemical elements C, H, N, and O are provided: energies, atomic forces, multipole moments, atomic charges, etc. We provide this data to the community to aid research and development of ML models for chemistry.

Computational investigation on the reaction of dimethyl ether with nitric dioxide. II. Detailed chemical kinetic modeling

Hydrogen abstraction from dimethyl ether (DME) by nitric dioxide (NO2) has been considered to be very sensitive to the DME ignition behavior. Thermal rate coefficients for the reaction of DME with NO2 are computed by employing the multi-structural canonical variational transition-state theory with multidimensional tunneling (MS-CVT/MT) at temperatures from 200 to 3000 K. The B2PLYP/TZVP theoretical level has been verified in the first part of this study against experiment for geometries of stationary points and against accurate CCSD(T)/CBS and CCSD(T)/6-311+G(2df,2p) computations for the energetics related to the title reaction, and then the compound CCSD(T)/6-311+G(2df,2p)//B2PLYP/TZVP computations are performed here to obtain the properties of all stationary points and the geometries, energies, gradients, and Hessians of non-stationary points along each reaction minimum-energy pathway. It has been proposed that the internal rotations in the transition structures are anharmonic and strongly coupled to each other to produce multiple conformations whose contributions are included in their partition function approximations. The previous evaluations for these rate coefficients based on transition-state theory with one-dimensional tunneling corrections and one-dimensional hindered rotor treatments of torsions are essentially distinct from our MS-CVT/SCT results where the critical effects of corner cutting due to small reaction-path curvature are seriously treated via multidimensional small-curvature tunneling (SCT) treatment and torsions are explicitly considered by multistructural treatments, and thus the MS-CVT/SCT rate coefficients are proposed particularly to refine the current kinetic model for DME combustion.

Impact of Approximate DFT Density Delocalization Error on Potential Energy Surfaces in Transition Metal Chemistry.

For approximate density functional theory (DFT) to be useful in catalytic applications of transition metal complexes, modeling strategies must simultaneously address electronic, geometric, and energetic properties of the relevant species. We show that for representative transition metal triatomics (MO2, where M = Cr, Mn, Fe, Co, or Ni) and related diatomics the incorporation of Hartree-Fock (HF) exchange in most cases improves the properties of the Born-Oppenheimer potential energy surface (PES) with respect to accurate experimental or CCSD(T) references. We rationalize this observation by noting reduced delocalization obtained with hybrid functionals (20-40% HF exchange), as evidenced by reduced hybridization of non-bonding orbitals and increases in metal partial charges. Although we show that the optimal exchange fraction is both property and system specific, incorporating HF exchange synergistically improves properties of density, structure, and energetics within a single PES characterized by moderately covalent bonding. The same improvement is not observed in the ordering of MO2 spin states, as good agreement of semi-local DFT spin state ordering is worsened by over-stabilization of higher spin states when HF exchange is added. More work is needed to understand minimal functional forms capable of improving multiple properties with respect to semi-local DFT descriptions of transition metal chemistry.

Theoretical study of H-atom abstraction reactions from CH3CH2OCH2CH3, CHF2CF2OCH2CF3 and CF3CH2OCH3 by NO3 radical & subsequent degradation

The nocturnal reactions of CH3CH2OCH2CH3, CHF2CF2OCH2CF3 and CF3CH2OCH3 initiated by NO3 radicals are important sources of alkyl radicals and nitric acids. In this paper, the thermodynamics and kinetics of CH3CH2OCH2CH3, CHF2CF2OCH2CF3 and CF3CH2OCH3 induced by NO3 radical in gas phase are studied in detail by BHandHLYP method combined with 6-311G(d,p) basis set, and the single point correction is calculated by relatively accurate CCSD(T) method. In the temperature range of 200–400 K, the rate constants of title reactions are fitted to the three-parameter Arrhenius formula:k1 = 1.13 × 10−40T9.24exp(1675.99/T)k2 = 2.23 × 10−23T2.81exp(-4476.24/T)k3 = 5.63 × 1043T−19.20exp(-9344.12/T)All the rate constants calculated by the canonical variational transition state theory and the small curvature tunneling are basically consistent with the limited experimental data. By comparing the reaction rate constants of ethyl ether and its isomer methyl propyl ether with NO3 radical at 293 ± 2 K, the higher the symmetry is, the faster the reaction rate of ether is.Thermodynamic calculations and kinetic data of the title reactions indicate that the H-abstraction reactions at the –OCH2- sites are the main reaction pathways. The thermodynamic and kinetic data of the reaction CH3CH2OCH2CH3, CHF2CF2OCH2CF3 with NO3 radical, showing that the reaction activity could be reduced due to the addition of fluorine atoms, which is further verified by the enthalpies, Gibbs free energies of the title reactions and C–H bond dissociation energies of the CH3CH2OCH2CH3, CHF2CF2OCH2CF3 and CF3CH2OCH3 molecules. The reaction thermodynamics and kinetics are determined, and the formation mechanisms of the products are proposed, which are crucial to determine the influence of CH3CH2OCH2CH3, CHF2CF2OCH2CF3 and CF3CH2OCH3 on air quality, as well as its atmospheric lifetime and durability.The atmospheric lifetimes of CH3CH2OCH2CH3, CHF2CF2OCH2CF3 and CF3CH2OCH3 are evaluated in the NO3-concentration range of 5 × 108-2 × 109 molecule cm−3 to fully consider the effects of different regions on their nocturnal migration. The radiation efficiency and global warming potentials (GWPs) have been reported. The products of title reaction CH3CH2OCHCH3, CF3CHOCF2CHF2 and CF3CHOCH3 are further oxidized into organic nitrates in the presence of O2 and NO. Organic nitrites can be isomerized into organic nitrates or degraded to form CH3C(O•)HOCH2CH3, CF3C(O•)HOCF2CHF2 and CH3C(O•)HOCF3 alkoxy radicals and NO2. This work provides deep insight into the night migration and transformation mechanism of the three ethers.

Description of noncovalent interactions involving π-system with high precision: An assessment of RPA, MP2, and DFT-D methods.

Efficient approaches with high precision are essential for understanding the formation and stability of noncovalent interaction complexes. Here, 21 noncovalent interaction complexes involving π-system are selected and grouped in three subsets according to ETS-NOCV method: dispersion-dominated, electrostatic-dominated, and mixed. We mainly focus on examining the performance of random-phase approximation (RPA) on these π systems. The tested RPA-based method includes standard RPA and its variants including the related single excitations (SEs), renormalized single excitations (rSEs), second-order screened exchange (SOSEX), and the renormalized second-order perturbation theory (rPT2). The routine second-order Møller-Plesset perturbation theory (MP2) and three popular DFT-D functionals (M06-2X-D3, ωB97XD, and PBE-D3(BJ)) are also assessed for comparison. In this work, besides the calculation of interaction energies at Dunning-type aug-cc-pVDZ and aug-cc-pVTZ basis set, we also present a larger database of interaction energies calculated using MP2 and RPA methods with Dunning-type aug-cc-pVQZ basis set. An accurate CCSD(T)/CBS scheme is used to provide benchmark database. In addition to the high-level results, we also provide potential energy surfaces (PES) of different interaction type. Among all the tested methods, MP2 has a satisfactory performance on electrostatic-dominated and mixed-type systems, except for dispersion-dominated systems. DFT-D functionals, especially ωB97XD functional, has a balanced performance across all the tested systems. Importantly, for RPA-based methods, the calculation accuracy can be dramatically improved by taking into account SE or exchange effects, especially in the mixed complexes. We conclude that rPT2 among all the test RPA-based methods gives an overall satisfactory performance across different interaction types. © 2019 Wiley Periodicals, Inc.

Understanding the Nature of Weak Interactions between Functionalized Boranes and N2/O2, Promising Functional Groups for Gas Separations.

The separation of nitrogen and oxygen gases is considered as a very challenging process, since both O2 and N2 are nonpolar molecules with similar kinetic diameters. Electronic structure theory can provide a fundamental understanding of effects that can lead to selective binding of nitrogen or oxygen gas for the development of novel separation processes. Boranes can bind dinitrogen through a dative bond, where the boron acts as a σ acceptor and back-donates through π orbitals. To better understand these interactions, we have performed highly accurate CCSD(F12)(T) and CCSDT(Q) computations for the BH3-N2 and BH3-O2 complexes. The coupled-cluster binding energies were used as reference for benchmarking different density functionals, and larger functionalized boranes were examined at the M05/def2-TZVPPD level. Symmetry adapted perturbation theory (SAPT) calculations were performed for the elucidation of the nature of the interaction between nitrogen and substituted boranes and how direct or distal functionalizations affect the strength of the weak dative bonds. By use of these methods, several boranes were found to bind N2 over O2. These molecular species are promising functional groups for incorporation into the next generation of advanced materials for efficient N2/O2 separations.

Basis Set Effects in the Description of the Cl-O Bond in ClO and XClO/ClOX Isomers (X = H, O, and Cl) Using DFT and CCSD(T) Methods

The performance of a group of density functional methods of progressive complexity for the description of the ClO bond in a series of chlorine oxides was investigated. The simplest ClO radical species and the two isomeric structures XClO/ClOX for each X = H, Cl, and O were studied using the PW91, TPSS, B3LYP, PBE0, M06, M06-2X, BMK, and B2PLYP functionals. Geometry optimizations and reaction enthalpies and enthalpies of formation for each species were calculated using Pople basis sets and the (aug)-cc-pVnZ Dunning sets, with n = D, T, Q, 5, and 6. For the calculation of enthalpies of formation, atomization and isodesmic reactions were employed. Both the precision of the methods with respect to the increase of the basis sets, as well as their accuracy, were gauged by comparing the results with the more accurate CCSD(T) calculations, performed using the same basis sets as for the DFT methods. The results obtained employing composite chemical methods (G4, CBS-QB3, and W1BD) were also used for the comparisons, as well as the experimental results when they are available. The results obtained show that error compensation is the key for successful description of molecular properties (geometries and energies) by carefully selecting the method and basis sets. In general, expansion of the one-electron basis set to the limit of completeness does not improve results at the DFT level, but just the opposite. The enthalpies of formation calculated at the CCSD(T)/aug-cc-pV6Z for the species considered are generally in agreement with experimental determinations and the most accurate theoretical values. Different sources of error in the calculations are discussed in detail.

Benchmarking DFT methods on linear and nonlinear electric properties of spatially confined molecules

AbstractThe accurate evaluation of molecular electric properties still constitutes a challenge for density functional theory (DFT) methods. Moreover, some new issues appear, within this area of research, that require analysis. One of them concerns the electric properties of confined matter. This study aims at a quantitative evaluation of the performance of DFT framework in the description of dipole moment, polarizability, and first hyperpolarizability of the isolated and spatially confined molecules. The calculations were performed using a set of 35 exchange‐correlation functionals. The reference data were determined employing the highly accurate CCSD(T) method. To mimic the effect of spatial confinement the cylindrically symmetric harmonic oscillator potential was applied. Based on the conducted research it was found that the behavior of the considered functionals substantially differs depending on the studied electric property and molecule. Moreover, it was observed that the strength of spatial confinement strongly influences the performance of DFT methods.