Linear-Scaling Implementation of Multilevel Hartree-Fock Theory.
Ontology highlight
ABSTRACT:
SUBMITTER: Goletto L
PROVIDER: S-EPMC8675138 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The ability to perform ab initio electronic structure calculations that scales linearly with the system size is one of the central aims in theoretical chemistry. In this study, the implementation of the divide-and-conquer (DC) algorithm, an algorithm with the potential to aid the achievement of true linear scaling within Hartree-Fock (HF) theory is revisited. Standard HF calculations solve the Roothaan-Hall equations for the whole system; in the DC-HF approach, the diagonalization of the Fock matrix is carried out on smaller subsystems. The DC algorithm for HF calculations was validated on polyglycines, polyalanines and eleven real three-dimensional proteins of up to 608 atoms in this work. We also found that a fragment-based initial guess using molecular fractionation with conjugated caps (MFCC) method significantly reduces the number of SCF cycles and even is capable of achieving convergence for some globular proteins where the simple superposition of atomic densities (SAD) initial guess fails.
| S-EPMC2853773 | biostudies-literature
Project description:The X-ray constrained wavefunction (XC-WF) method proposed by Jayatilaka [Jayatilaka & Grimwood (2001) ▸, Acta Cryst. A57, 76-86] has attracted much attention because it represents a possible third way of theoretically studying the electronic structure of atoms and molecules, combining features of the more popular wavefunction- and DFT-based approaches. In its original formulation, the XC-WF technique extracts statistically plausible wavefunctions from experimental X-ray diffraction data of molecular crystals. A weight is used to constrain the pure Hartree-Fock solution to the observed X-ray structure factors. Despite the wavefunction being a single Slater determinant, it is generally assumed that its flexibility could guarantee the capture, better than any other experimental model, of electron correlation effects, absent in the Hartree-Fock Hamiltonian but present in the structure factors measured experimentally. However, although the approach has been known for long time, careful testing of this fundamental hypothesis is still missing. Since a formal demonstration is impossible, the validation can only be done heuristically and, to accomplish this task, X-ray constrained Hartree-Fock calculations have been performed using structure factor amplitudes computed at a very high correlation level (coupled cluster) for selected molecules in isolation, in order to avoid the perturbations due to intermolecular interactions. The results show that a single-determinant XC-WF is able to capture the electron correlation effects only partially. The largest amount of electron correlation is extracted when: (i) a large external weight is used (much larger than what has normally been used in XC-WF calculations using experimental data); and (ii) the high-order reflections, which carry less information on the electron correlation, are down-weighted (or even excluded), otherwise they would bias the fitting towards the unconstrained Hartree-Fock wavefunction.
| S-EPMC5330524 | biostudies-literature
Project description:Automatic differentiation (AD) is a powerful tool that allows calculating derivatives of implemented algorithms with respect to all of their parameters up to machine precision, without the need to explicitly add any additional functions. Thus, AD has great potential in quantum chemistry, where gradients are omnipresent but also difficult to obtain, and researchers typically spend a considerable amount of time finding suitable analytical forms when implementing derivatives. Here, we demonstrate that AD can be used to compute gradients with respect to any parameter throughout a complete quantum chemistry method. We present DiffiQult, a Hartree-Fock implementation, entirely differentiated with the use of AD tools. DiffiQult is a software package written in plain Python with minimal deviation from standard code which illustrates the capability of AD to save human effort and time in implementations of exact gradients in quantum chemistry. We leverage the obtained gradients to optimize the parameters of one-particle basis sets in the context of the floating Gaussian framework.
| S-EPMC5968443 | biostudies-literature
Project description:Linear-scaling density functional theory (DFT) is an efficient method to describe the electronic structures of molecules, semiconductors, and insulators to avoid the high cubic-scaling cost in conventional DFT calculations. Here, we present a parallel implementation of linear-scaling density matrix trace correcting (TC) purification algorithm to solve the Kohn-Sham (KS) equations with the numerical atomic orbitals in the HONPAS package. Such a linear-scaling density matrix purification algorithm is based on the Kohn's nearsightedness principle, resulting in a sparse Hamiltonian matrix with localized basis sets in the DFT calculations. Therefore, sparse matrix multiplication is the most time-consuming step in the density matrix purification algorithm for linear-scaling DFT calculations. We propose to use the MPI_Allgather function for parallel programming to deal with the sparse matrix multiplication within the compressed sparse row (CSR) format, which can scale up to hundreds of processing cores on modern heterogeneous supercomputers. We demonstrate the computational accuracy and efficiency of this parallel density matrix purification algorithm by performing large-scale DFT calculations on boron nitrogen nanotubes containing tens of thousands of atoms.
| S-EPMC7726133 | biostudies-literature
Project description:The counterpoise (CP) correction by Boys and Bernardi has been well accepted as a reliable strategy to account for basis set superposition error (BSSE) in intermolecular complexes. The behavior of the CP correction was thoroughly studied in individual molecules of molecular complexes. This work studies the performance of the CP correction in many-body clusters including three-body clusters of organic compounds in the 3B-69 dataset. Additionally, we used crystal structures of polymorphs of benzene, aspirin, and oxalyl dihydrazide (ODH) to construct a many-body cluster dataset, abbreviated as the MBC-36 dataset, consisting of two, four and eight molecules, and 16 molecules in the case of benzene. A series of Dunning's basis sets-cc-pXZ and aug-cc-pXZ (X = D and T)-were used to predict CP-corrected Hartree-Fock (HF) interaction energies of the 3B-69 and MBC-36 datasets. The CP-corrected interaction energies were found to be basis-set independent, whereas the non-CP corrected interaction energies were found not to a follow a smooth exponential fitting as previously found for electronic energies of individual molecules. This observation was attributed to the presence of non-additive induction forces in some clusters. Two 2 × 2 × 2 supercells of benzene polymorphs were constructed to explore the local nature of BSSE effects. A cut-off radius of 10 Å was demonstrated to be sufficient to fully recover these effects. Although the behavior of CP correction was found to be non-conventional in many-body clusters of organic compounds, the use of a small basis set such as cc-pVDZ showed excellent performance in the prediction of HF interaction energies.
| S-EPMC9303541 | biostudies-literature
Project description:ObjectiveThe present study aims to establish the relationship between the reported radical-scavenging activities of flavonoids and some enthalpy changes that may occur during flavonoids' reactions with free radicals.MethodEight flavonoids were chosen for the study on the basis of their structural merits and reported 1,1-diphenyl-2-picryl-hydrazyl scavenging activities. Enthalpy changes accompanying interconversions between selected conformations (including spin multiplicities) and homolytic dissociations were estimated.ResultsA novel relationship exists between the total enthalpy of reaction for the abstraction of two hydrogen atoms from flavonoids, their reported radical-scavenging activities and the enthalpy of the homolytic dissociation of hydrogen molecule (104.206 kcal mol(-1)). Only those flavonoids which could give up two hydrogen atoms with total enthalpy changes well below 104.206 kcal mol(-1) were active radical scavengers.DiscussionBy appealing to equilibrium dynamics, we demonstrated that, for flavonoids to be able to donate hydrogen atoms, the change in enthalpy accompanying the abstraction of two hydrogen atoms needs to be less than 104.206 kcal mol(-1). This condition does not seem to be restricted to flavonoids only but rather generally applicable to chian-breaking antioxidants.ConclusionThermodynamical relationships may be the most important factors governing the radical-scavenging reactions of flavonoids and possibly other compounds as well. Nevertheless, a more complete characterization of antioxidants would necessitate kinetic analysis.
| S-EPMC6837546 | biostudies-literature
Project description:The full configuration interaction (full-CI) method is capable of providing the numerically best wave functions and energies of atoms and molecules within basis sets being used, although it is intractable for classical computers. Quantum computers can perform full-CI calculations in polynomial time against the system size by adopting a quantum phase estimation algorithm (QPEA). In the QPEA, the preparation of initial guess wave functions having sufficiently large overlap with the exact wave function is recommended. The Hartree-Fock (HF) wave function is a good initial guess only for closed shell singlet molecules and high-spin molecules carrying no spin-? unpaired electrons, around their equilibrium geometry, and thus, the construction of multiconfigurational wave functions without performing post-HF calculations on classical computers is highly desired for applying the method to a wide variety of chemistries and physics. In this work, we propose a method to construct multiconfigurational initial guess wave functions suitable for QPEA-based full-CI calculations on quantum computers, by utilizing diradical characters computed from spin-projected UHF wave functions. The proposed approach drastically improves the wave function overlap, particularly in molecules with intermediate diradical characters.
| S-EPMC6346626 | biostudies-literature
Project description:Polyphenolic compounds are now widely studied using computational chemistry approaches, the most popular of which is Density Functional Theory. To ease this process, it is critical to identify the optimal level of theory in terms of both accuracy and resource usage-a challenge we tackle in this study. Eleven DFT functionals with varied Hartree-Fock exchange values, both global and range-separated hybrids, were combined with 14 differently augmented basis sets to calculate the reactivity indices of caffeic acid, a phenolic acid representative, and compare them to experimental data or a high-level of theory outcome. Aside from the main course, a validation of the widely used Janak's theorem in the establishment of vertical ionization potential and vertical electron affinity was evaluated. To investigate what influences the values of the properties under consideration, linear regression models were developed and thoroughly discussed. The results were utilized to compute the scores, which let us determine the best and worst combinations and make broad suggestions on the final option. The study demonstrates that M06-2X/6-311G(d,p) is the best fit for such research, and, curiously, it is not necessarily essential to include a diffuse function to produce satisfactory results.
| S-EPMC8398206 | biostudies-literature
Project description:Following recent developments in multilevel embedding methods, we introduce a novel density matrix-based multilevel approach within the framework of density functional theory (DFT). In this multilevel DFT, the system is partitioned in an active and an inactive fragment, and all interactions are retained between the two parts. The decomposition of the total system is performed upon the density matrix. The orthogonality between the two parts is maintained by solving the Kohn-Sham equations in the MO basis for the active part only, while keeping the inactive density matrix frozen. This results in the reduction of computational cost. We outline the theory and implementation and discuss the differences and similarities with state-of-the-art DFT embedding methods. We present applications to aqueous solutions of methyloxirane and glycidol.
| S-EPMC7880574 | biostudies-literature
Project description:Networks distribute energy, materials and information to the components of a variety of natural and human-engineered systems, including organisms, brains, the Internet and microprocessors. Distribution networks enable the integrated and coordinated functioning of these systems, and they also constrain their design. The similar hierarchical branching networks observed in organisms and microprocessors are striking, given that the structure of organisms has evolved via natural selection, while microprocessors are designed by engineers. Metabolic scaling theory (MST) shows that the rate at which networks deliver energy to an organism is proportional to its mass raised to the 3/4 power. We show that computational systems are also characterized by nonlinear network scaling and use MST principles to characterize how information networks scale, focusing on how MST predicts properties of clock distribution networks in microprocessors. The MST equations are modified to account for variation in the size and density of transistors and terminal wires in microprocessors. Based on the scaling of the clock distribution network, we predict a set of trade-offs and performance properties that scale with chip size and the number of transistors. However, there are systematic deviations between power requirements on microprocessors and predictions derived directly from MST. These deviations are addressed by augmenting the model to account for decentralized flow in some microprocessor networks (e.g. in logic networks). More generally, we hypothesize a set of constraints between the size, power and performance of networked information systems including transistors on chips, hosts on the Internet and neurons in the brain.
| S-EPMC2607348 | biostudies-literature