A Minimum Quantum Chemistry CCSD(T)/CBS Data Set of Dimeric Interaction Energies for Small Organic Functional Groups: Heterodimers
Ontology highlight
ABSTRACT:
SUBMITTER: Huang H
PROVIDER: S-EPMC9201891 | biostudies-literature | 2022 May
REPOSITORIES: biostudies-literature
Similar Datasets
Project description: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+AT1/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.
| S-EPMC11008106 | biostudies-literature
Project description:A thorough comparison is carried out for explicitly correlated and density-based basis-set correction approaches, which were primarily developed to mitigate the basis-set incompleteness error of wave function methods. An efficient implementation of the density-based scheme is also presented, utilizing the density-fitting approximation. The performance of these approaches is comprehensively tested for the second-order Møller-Plesset (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] methods with respect to the corresponding complete basis set references. It is demonstrated that the density-based correction together with complementary auxiliary basis set (CABS)-corrected Hartree-Fock energies is highly robust and effectively reduces the error of the standard approaches; however, it does not outperform the corresponding explicitly correlated methods. Nevertheless, what still makes the density-corrected CCSD and CCSD(T) methods competitive is that their computational costs are roughly half of those of the corresponding explicitly correlated variants. Additionally, an incremental approach for standard CCSD and CCSD(T) is introduced. In this simple scheme, the total energies are corrected with the CABS correction and explicitly correlated MP2 contributions. As demonstrated, the resulting methods yield surprisingly good results, below 1 kcal/mol for thermochemical properties even with a double-ζ basis, while their computational expenses are practically identical to those of the density-based basis-set correction approaches.
| S-EPMC10688194 | biostudies-literature
Project description:The calculation of density-based basis-set correction (DBBSC), which remedies the basis-set incompleteness (BSI) error of the correlation energy, is combined with local approximations. Aiming at large-scale applications, the procedure is implemented in our efficient local natural orbital-based coupled-cluster singles and doubles with perturbative triples [LNO-CCSD(T)] scheme. To this end, the range-separation function, which characterizes the one-electron BSI in space, is decomposed into the sum of contributions from individual localized molecular orbitals (LMOs). A compact domain is constructed around each LMO, and the corresponding contributions are evaluated only within these restricted domains. Furthermore, for the calculation of the complementary auxiliary basis set (CABS) correction, which significantly improves the Hartree-Fock (HF) energy, the local density fitting approximation is utilized. The errors arising from the local approximations are examined in detail, efficient prescreening techniques are introduced to compress the numerical quadrature used for DBBSC, and conservative default thresholds are selected for the truncation parameters. The efficiency of the DBBSC-LNO-CCSD(T) method is demonstrated through representative examples of up to 1000 atoms. Based on the numerical results, we conclude that the corrections drastically reduce the BSI error using double-ζ basis sets, often to below 1 kcal/mol compared to the reliable LNO-CCSD(T) complete basis set references, while significant improvements are also achieved with triple-ζ basis sets. Considering that the calculation of the DBBSC and CABS corrections only moderately increases the wall-clock time required for the post-HF steps in practical applications, the proposed DBBSC-LNO-CCSD(T) method offers a highly efficient and robust tool for large-scale calculations.
| S-EPMC11391584 | biostudies-literature
Project description:Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).
| S-EPMC9835832 | biostudies-literature
Project description: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.
| S-EPMC7884001 | biostudies-literature
Project description:Membranes fuse by forming highly curved intermediates, culminating in structures described as fusion pores. These hourglass-like figures that join two fusing membranes have high bending energies, which can be estimated using continuum elasticity models. Fusion pore bending energies depend strongly on shape, and the present study developed a method for determining the shape that minimizes bending energy. This was first applied to a fusion pore modeled as a single surface and then extended to a more realistic model treating a bilayer as two monolayers. For the two-monolayer model, fusion pores were found to have metastable states with energy minima at particular values of the pore diameter and bilayer separation. Fusion pore energies were relatively insensitive to membrane thickness but highly sensitive to spontaneous curvature and membrane asymmetry. With symmetrical bilayers and monolayer spontaneous curvatures of -0.1 nm(-1) (a typical value) separated by 6 nm (closest distance determined by repulsive hydration forces), fusion pore formation required 43-65 kT. The pore radius of approximately 2.25 nm fell within the range estimated from conductance measurements. With bilayer separation >6 nm, fusion pore formation required less energy, suggesting that protein scaffolds can promote fusion by bending membranes toward one another. With nonzero spontaneous monolayer curvature, the shape that minimized the energy change during fusion pore formation differed from the shape that minimized its energy after it formed. Thus, a nascent fusion pore will relax spontaneously to a new shape, consistent with the experimentally observed expansion of nascent fusion pores during viral fusion.
| S-EPMC2833281 | biostudies-literature
Project description:A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids "on-the-fly" QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.
| S-EPMC2736613 | biostudies-other
Project description:A vast number of non-covalent interaction energies at the counterpoise corrected CCSD(T) level have been collected from the literature to build a diverse new dataset. The whole dataset, which consists of 2027 CCSD(T) energies, includes most of the published data at this level. A large subset of the data was then used to train a novel, B3LYP specific, empirical correction scheme for non-covalent interactions and basis set superposition error (abbreviated as B3LYP-MM). Results obtained with our new correction scheme were directly compared to benchmark results obtained with B3LYP-D3(1) and M06-2X(2) (two popular density functions designed specifically to accurately model non-covalent interactions). For non-covalent complexes dominated by dispersion or dipole-dipole interactions all three tested methods give accurate results with the medium size aug-cc-pVDZ(3-6) basis set with MUE's of 0.27 (B3LYP-MM), 0.32 (B3LYP-D3) and 0.47 kcal/mol (M06-2X) (with explicit counterpoise corrections). These results validate both B3LYP-D3 and M06-2X for interactions of this type using a much larger data set than was presented in prior work. However, our new dispersion correction scheme shows some clear advantages for dispersion and dipole-dipole dominated complexes with the small LACVP* basis set, which is very popular in use due to its low associated computational cost: The MUE for B3LYP-MM with the LACVP* basis set for this subset of complexes (without explicit counterpoise corrections) is only 0.28 kcal/mol, compared to 0.65 kcal/mol for M06-2X or 1.16 kcal/mol for B3LYP-D3. Additionally, our new correction scheme also shows major improvements in accuracy for hydrogen-bonded systems and for systems involving ionic interactions, for example cation-π interactions. Compared to B3LYP-D3 and M06-2X, we also find that our new B3LYP-MM correction scheme gives results of higher or equal accuracy for a large dataset of conformer energies of di- and tripeptides, sugars, and cysteine.
| S-EPMC3206731 | biostudies-literature
Project description:We have investigated the title question for both a subset of the W4-11 total atomization energies benchmark, and for the A24x8 noncovalent interactions benchmark. Overall, counterpoise corrections to post-CCSD(T) contributions are about two orders of magnitude less important than those to the CCSD(T) interaction energy. Counterpoise corrections for connected quadruple substitutions (Q) are negligible, and QΛ-Q$$ {(Q)}_{\Lambda}-(Q) $$ or T4-Q$$ {T}_4-(Q) $$ especially so. In contrast, for atomization energies, the T3-T$$ {T}_3-(T) $$ counterpoise correction can reach about 0.05 kcal/mol for small basis sets like cc-pVDZ, thought it rapidly tapers off with cc-pVTZ and especially aug-cc-pVTZ basis sets. It is reduced to insignificance by the extrapolation of T3-T$$ {T}_3-(T) $$ applied in both W4 and HEAT thermochemistry protocols. In noncovalent dimers, the differential BSSE on post-CCSD(T) correlation contributions is negligible even in basis sets as small as the unpolarized split-valence cc-pVDZ(no d).
| S-EPMC11668233 | biostudies-literature
Project description:The non-Boltzmann Bennett (NBB) free energy estimator method is applied to 21 molecules from the blind subset of the SAMPL4 challenge. When NBB is applied with the SMD implicit solvent model, and the OLYP/DZP level of quantum chemistry, highly accurate hydration free energy calculations are obtained with respect to experiment (RMSD=0.89kcal·mol-1). Other quantum chemical methods are also tested, and the effects of solvent model, density functional, basis set are explored in this benchmarking study, providing a framework for improvements in calculating hydration free energies. We provide a practical guide for using the best QM-NBB protocols that are consistently more accurate than either pure QM or pure MM alone. In situations where high accuracy hydration free energy predictions are needed, the QM-NBB method with SMD implicit solvent should be the first choice of quantum chemists.
| S-EPMC5068830 | biostudies-literature