Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA

The p-block challenge: assessing quantum chemistry methods for inorganic heterocycle dimerizations

The elements of the p-block of the periodic table are of high interest in various chemical and technical applications like frustrated Lewis-pairs (FLP) or opto-electronics. However, high-quality benchmark data to assess approximate density functional theory (DFT) for their theoretical description are sparse. In this work, we present a benchmark set of 604 dimerization energies of 302 "inorganic benzenes" composed of all non-carbon p-block elements of main groups III to VI up to polonium. This so-called IHD302 test set comprises two classes of structures formed by covalent bonding and by weaker donor-acceptor (WDA) interactions, respectively. Generating reliable reference data with ab initio methods is challenging due to large electron correlation contributions, core-valence correlation effects, and especially the slow basis set convergence. To compute reference values for these dimerization reactions, after thorough testing, we applied a computational protocol using state-of-the-art explicitly correlated local coupled cluster theory termed PNO-LCCSD(T)-F12/cc-VTZ-PP-F12(corr.). It includes a basis set correction at the PNO-LMP2-F12/aug-cc-pwCVTZ level. Based on these reference data, we assess 26 DFT methods in combination with three different dispersion corrections and the def2-QZVPP basis set, five composite DFT approaches, and five semi-empirical quantum mechanical methods. For the covalent dimerizations, the r2SCAN-D4 meta-GGA, the r2SCAN0-D4 and ωB97M-V hybrids, and the revDSD-PBEP86-D4 double-hybrid functional are found to be the best-performing methods among the evaluated functionals of the respective class. However, since def2 basis sets for the 4th period are not associated to relativistic pseudo-potentials, we obtained significant errors in the covalent dimerization energies (up to 6 kcal mol-1) for molecules containing p-block elements of the 4th period. Significant improvements were achieved for systems containing 4th row elements by using ECP10MDF pseudopotentials along with re-contracted aug-cc-pVQZ-PP-KS basis sets introduced in this work with the contraction coefficients taken from atomic DFT (PBE0) calculations. Overall, the IHD302 set represents a challenge to contemporary quantum chemical methods. This is due to a large number of spatially close p-element bonds which are underrepresented in other benchmark sets, and the partial covalent bonding character for the WDA interactions. The IHD302 set may be helpful to develop more robust and transferable approximate quantum chemical methods in the future.

Improved CPS and CBS Extrapolation of PNO-CCSD(T) Energies: The MOBH35 and ISOL24 Data Sets

Computation of heats of reaction of large molecules is now feasible using the domain-based pair natural orbital (PNO)-coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] theory. However, to obtain agreement within 1 kcal/mol of experiment, it is necessary to eliminate basis set incompleteness error, which comprises both the AO basis set error and the PNO truncation error. Our investigation into the convergence to the canonical limit of PNO-CCSD(T) energies with the PNO truncation threshold T shows that errors follow the model E ( T ) = E + A T 1 / 2 . Therefore, PNO truncation errors can be eliminated using a simple two-point CPS extrapolation to the canonical limit so that subsequent CBS extrapolation is not limited by the residual PNO truncation error. Using the ISOL24 and MOBH35 data sets, we find that PNO truncation errors are larger for molecules with significant static correlation and that it is necessary to use very tight thresholds of T = 10 - 8 to ensure that errors do not exceed 1 kcal/mol. We present a lower-cost extrapolation scheme that uses information from small basis sets to estimate the PNO truncation errors for larger basis sets. In this way, the canonical limit of CCSD(T) calculations on sizable molecules with large basis sets can be reliably estimated in a practical way. Using this approach, we report near complete basis set (CBS)-CCSD(T) reaction energies for the full ISOL24 and MOBH35 data sets.

CREST-A program for the exploration of low-energy molecular chemical space

Conformer-rotamer sampling tool (CREST) is an open-source program for the efficient and automated exploration of molecular chemical space. Originally developed in Pracht et al. [Phys. Chem. Chem. Phys. 22, 7169 (2020)] as an automated driver for calculations at the extended tight-binding level (xTB), it offers a variety of molecular- and metadynamics simulations, geometry optimization, and molecular structure analysis capabilities. Implemented algorithms include automated procedures for conformational sampling, explicit solvation studies, the calculation of absolute molecular entropy, and the identification of molecular protonation and deprotonation sites. Calculations are set up to run concurrently, providing efficient single-node parallelization. CREST is designed to require minimal user input and comes with an implementation of the GFNn-xTB Hamiltonians and the GFN-FF force-field. Furthermore, interfaces to any quantum chemistry and force-field software can easily be created. In this article, we present recent developments in the CREST code and show a selection of applications for the most important features of the program. An important novelty is the refactored calculation backend, which provides significant speed-up for sampling of small or medium-sized drug molecules and allows for more sophisticated setups, for example, quantum mechanics/molecular mechanics and minimum energy crossing point calculations.

Conformational energies of biomolecules in solution: Extending the MPCONF196 benchmark with explicit water molecules

A prerequisite for the computational prediction of molecular properties like conformational energies of biomolecules is a reliable, robust, and computationally affordable method usually selected according to its performance for relevant benchmark sets. However, most of these sets comprise molecules in the gas phase and do not cover interactions with a solvent, even though biomolecules typically occur in aqueous solution. To address this issue, we introduce a with explicit water molecules solvated version of a gas-phase benchmark set containing 196 conformers of 13 peptides and other relevant macrocycles, namely MPCONF196 [J. Řezáč et al., JCTC 2018, 14, 1254-1266], and provide very accurate PNO-LCCSD(T)-F12b/AVQZ' reference values. The novel solvMPCONF196 benchmark set features two additional challenges beyond the description of conformers in the gas phase: conformer-water and water-water interactions. The overall best performing method for this set is the double hybrid revDSDPBEP86-D4/def2-QZVPP yielding conformational energies of almost coupled cluster quality. Furthermore, some (meta-)GGAs and hybrid functionals like B97M-V and ω B97M-D with a large basis set reproduce the coupled cluster reference with an MAD below 1 kcal mol- 1  . If more efficient methods are required, the composite DFT-method r2 SCAN-3c (MAD of 1.2 kcal mol- 1 ) is a good alternative, and when conformational energies of polypeptides or macrocycles with more than 500-1000 atoms are in the focus, the semi-empirical GFN2-xTB or the MMFF94 force field (for very large systems) are recommended.

Confined Lewis Pairs: Investigation of the X- →Si20 Interaction in Halogen-Encapsulating Silafulleranes

A joint theoretical and experimental study on 32 endohedral silafullerane derivatives [X@Si20 Y20 ]- (X=F-I; Y=F-I, H, Me, Et) andT h ${T_h }$ -[Cl@Si20 H12 Y8 ]- (Y=F-I) is presented. First, we evaluated the structure-determining template effect of Cl- in a systematic series of concave silapolyquinane model systems. Second, we investigated the X- →Si20 interaction energy (E int ${E_{{\rm{int}}} }$ ) as a function of X- and Y and found the largestE int ${E_{{\rm{int}}} }$ values for electron-withdrawing exohedral substituents Y. Given that X- ions can be considered as Lewis bases and empty Si20 Y20 clusters as Lewis acids, we classify our inseparable host-guest complexes [X@Si20 Y20 ]- as "confined Lewis pairs". Third, 35 Cl NMR spectroscopy proved to be highly diagnostic for an experimental assessment of the Cl- →Si20 interaction as the paramagnetic shielding and, in turn, δ ${\delta }$ (35 Cl) of the endohedral Cl- ion correlate inversely withE int ${E_{{\rm{int}}} }$ . Finally, we disclose the synthesis of [PPN][Cl@Si20 Y20 ] (Y=Me, Et, Br) and provide a thorough characterization of these new silafulleranes.

Influence of Steric and Dispersion Interactions on the Thermochemistry of Crowded (Fluoro)alkyl Compounds

ConspectusAlkanes play a pivotal role in industrial, environmental, and biological processes. They are characterized by their carbon-carbon single-bond structure, remarkable stability, and conformational diversity. Fluorination of such compounds imparts unique physicochemical properties that often enhance pharmacokinetic profiles, metabolic stability, and receptor interactions while keeping beneficial properties. However, such per- and polyfluoroalkyl substances (PFAS) show a persistent presence in the environment and potential adverse health effects, which propelled them to the forefront of global environmental and health discussions. Alkyl compounds are also prototypical for stereoelectronic (SE) effects that are widely applied in chemistry. Substituents are typically described as electron-density-donating/withdrawing and/or responsible for sterically interacting with reagents or strategic groups in the molecule. That alkane branching can result in higher stability compared to less-branched isomers has been investigated in detail also by testing quantum chemical methods, in particular density functional theory (DFT). Alkane branching results in spatially compact structures with close intramolecular contacts so that at a specific size the detailed balance of attractive London dispersion and covalent versus repulsive Pauli exchange interactions shifts to new, chemically unfragile situations. This may lead to dissociation at room temperature and opens the central question: what is the smallest crowed alkane that cannot be made synthetically? In this Account, we try to shed light on the interplay among the various (free) energy components for crowded (fluoro)alkane dissociation. In this context, homolytic cleavage of the central C-C bond in a series of model alkanes of increasing size with tert-butyl (tBu), adamantyl (Ad), and [1.1.1]propellanyl (Prop) substituents is investigated. Reference energies are calculated at the PNO-LCCSD(T)-F12b level and used to benchmark the performance of contemporary DFT functionals. In line with previous conclusions, the application of dispersion corrections to density functionals is mandatory. For crowed structures, the accurate description of the midrange correlation effects, specifically repulsive van der Waals interactions, is crucial, and we observed that the density-dependent VV10 correction is superior to D4 in this context, although the asymptotic region is better described by the latter. The best available dispersion-inclusive functionals show systematic and reasonably small residual errors and can be safely applied to large systems (>100 atoms), for which coupled cluster methods with large basis sets are not computationally feasible anymore. For qualitatively correct predictions of synthetic accessibility under equilibrium conditions (free energy), the inclusion of thermostatistical (entropy) contributions is also essential. According to our results, tetra-tert-butylmethane (C17tBu) is the largest and most crowded system with a positive dissociation free energy and should be synthesizable. The difference between hydrogenated and perfluorinated systems originates from the increase in the steric repulsion of spatially close substituents, which is not compensated to the same extent by attractive orbital and dispersion interactions. A sometimes-assumed similar steric demand for fluorine and hydrogen atoms is not corroborated by our investigations on crowded systems. Perfluorination is found to substantially decrease thermal stability, rendering perfluorinated hexamethylethane (C8tBuF) the last potentially stable representative.

Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules

In this article, we present a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework, including explicitly correlated second-order Møller-Plesset perturbation (MP2), coupled-cluster singles and doubles (CCSD), domain-based local pair natural orbital CCSD (DLPNO-CCSD), and DLPNO-CCSD with perturbative triples (DLPNO-CCSD(T)). In these methods, F12 correction is decomposed into contributions from each occupied local molecular orbital and then evaluated independently in a given cluster, which consists of a subset of localized orbitals. These newly developed methods allow F12 calculations of large molecules (up to 145 atoms for quasi-one-dimensional systems) on a single node. We use these methods to investigate the relative stability between extended and folded alkane C30H62, the relative stability of four secondary structures of a polyglycine Ace(Gly)10NH2, and the binding energies of two host-guest complexes. The results demonstrate that the combination of CIM with F12 methods is a promising way to investigate large molecules with small basis set errors.

Assessment of DLPNO-MP2 Approximations in Double-Hybrid DFT

The unfavorable scaling (N5) of the conventional second-order Møller-Plesset theory (MP2) typically prevents the application of double-hybrid (DH) density functionals to large systems with more than 100 atoms. A prominent approach to reduce the computational demand of electron correlation methods is the domain-based local pair natural orbital (DLPNO) approximation that is successfully used in the framework of DLPNO-CCSD(T). Its extension to MP2 [Pinski P.; Riplinger, C.; Valeev, E. F.; Neese, F. J. Chem. Phys. 2015, 143, 034108.] paved the way for DLPNO-based DH (DLPNO-DH) methods. In this work, we assess the accuracy of the DLPNO-DH approximation compared to conventional DHs on a large number of 7925 data points for thermochemistry and 239 data points for structural features, including main-group and transition-metal systems. It is shown that DLPNO-DH-DFT can be applied successfully to perform energy calculations and geometry optimizations for large molecules at a drastically reduced computational cost. Furthermore, PNO space extrapolation is shown to be applicable, similar to its DLPNO-CCSD(T) counterpart, to reduce the remaining error.

Accurate Thermochemical and Kinetic Parameters at Affordable Cost by Means of the Pisa Composite Scheme (PCS)

A new strategy for the computation at an affordable cost of geometrical structures, thermochemical parameters, and rate constants for medium-sized molecules in the gas phase is proposed. The most distinctive features of the new model are the systematic use of cc-pVnZ-F12 basis sets, the addition of MP2 core-valence correlation in geometry optimizations by a double-hybrid functional, the separate extrapolation of MP2 and post-MP2 contributions, and the inclusion of anharmonic contributions in zero-point energies and thermodynamic functions. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme outperforms the most well-known model chemistries without the need for any empirical parameter. Additional tests show that the computed zero-point energies and thermal contributions can be confidently used for obtaining accurate thermochemical and kinetic parameters. Since the whole computational workflow is translated in a black-box procedure, which can be followed with standard electronic structure codes, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.