Dataset Information

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets.

ABSTRACT: We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (∼1.2 × 10^-3E_h/a₀) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

SUBMITTER: Shpiro B

PROVIDER: S-EPMC8908760 | biostudies-literature | 2022 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets.

Shpiro Ben B Fabian Marcel David MD Rabani Eran E Baer Roi R

Journal of chemical theory and computation 20220131 3

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. <i>Wiley Interdiscip. Rev.: Comput. Mol. Sci.</i> 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is t ...[more]

PMID: 35099187

Similar Datasets

Project description:Electron and charge transfers are part of many vital processes in nature and technology. Ab initio descriptions of these processes provide useful insights that can be utilized for applications. A combination of the embedded cluster material model and nonorthogonal configuration interaction (NOCI), in which the cluster wave functions are expanded in many-electron basis functions (MEBFs) consisting of spin-adapted, antisymmetrized products of multiconfigurational wave functions of fragments (which are usually molecules) in the cluster, appears to provide a compromise between accuracy and calculation time. Additional advantages of this NOCI-Fragments approach are the chemically convenient interpretation of the wave function in terms of molecular states, and the direct accessibility of electronic coupling between diabatic states to describe energy and electron transfer processes. Bottlenecks in this method are the large number of two-electron integrals that have to be handled for the calculation of an electronic coupling matrix element and the enormous number of matrix elements over determinant pairs that have to be evaluated for the calculation of one matrix element between the MEBFs. We show here how we created a reduced common molecular orbital basis that is utilized to significantly reduce the number of two-electron integrals that need to be handled. The results obtained with this basis do not show any loss of accuracy in relevant quantities like electronic couplings and vertical excitation energies. We also show a significant reduction in computation time without loss in accuracy when matrix elements over determinant pairs with small weights are neglected in the NOCI. These improvements in the methodology render NOCI-Fragments to be also applicable to treat clusters of larger molecular systems with larger atomic basis sets and larger active spaces, as the computation time becomes dependent on the number of occupied orbitals and less dependent on the size of the active space.

Project description:Recent advances in local electron correlation approaches have enabled the relatively routine access to CCSD(T) [that is, coupled cluster (CC) with single, double, and perturbative triple excitations] computations for molecules of a hundred or more atoms. Here, approaching their complete basis set (CBS) limit becomes more challenging due to extensive basis set superposition errors, often necessitating the use of large atomic orbital (AO) basis sets with diffuse functions. Here, we study a potential remedy in the form of non-atom-centered or floating orbitals (FOs). FOs are still rarely employed even for small molecules due to the practical complication of defining their position, number, exponents, etc. The most frequently used FO method thus simply places a single FO center with a large number of FOs toward the middle of noncovalent dimers; however, a single FO center for larger complexes can soon become insufficient. A recent alternative uses a grid of FO centers around the monomers with a single s function per center, which is currently applicable only for H, C, N, and O atoms. Here, we build on the above advantages and mitigate some drawbacks of previous FO approaches by using a layer of FO centers and 4-9 FOs/center for each monomer. Thus, a double layer of FOs is placed between the interacting subsystems. When extending the double-ζ AO basis with this double layer of FOs, the quality of conventional augmented double-ζ or conventional triple-ζ AO bases can be reached or surpassed with less orbitals, leading to few tenths of a kcal/mol basis set errors for medium-sized dimers. This good performance extends to larger molecules (shown here up to 72 atoms), as efficient local natural orbital (LNO) CCSD(T) computations with only double-ζ AO and 4 FOs/center FO bases match our LNO-CCSD(T)/CBS reference within ca. 0.1 kcal/mol. These developments introduce FO methods to the accurate modeling of large molecular complexes without limitations to atom types by further accelerating efficient correlation calculations, like LNO-CCSD(T).

Dataset Information

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets.

Publications

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets