36 research outputs found
A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry
We present here a review of the fundamental topics of Hartree-Fock theory in
Quantum Chemistry. From the molecular Hamiltonian, using and discussing the
Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock
equations for the electronic problem. Special emphasis is placed in the most
relevant mathematical aspects of the theoretical derivation of the final
equations, as well as in the results regarding the existence and uniqueness of
their solutions. All Hartree-Fock versions with different spin restrictions are
systematically extracted from the general case, thus providing a unifying
framework. Then, the discretization of the one-electron orbitals space is
reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition
of the basic underlying concepts related to the construction and selection of
Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we
close the review with a section in which the most relevant modern developments
(specially those related to the design of linear-scaling methods) are commented
and linked to the issues discussed. The whole work is intentionally
introductory and rather self-contained, so that it may be useful for non
experts that aim to use quantum chemical methods in interdisciplinary
applications. Moreover, much material that is found scattered in the literature
has been put together here to facilitate comprehension and to serve as a handy
reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and
subeqn package
Calculation of the Free Energy and Cooperativity of Protein Folding
Calculation of the free energy of protein folding and delineation of its pre-organization are of foremost importance for understanding, predicting and designing biological macromolecules. Here, we introduce an energy smoothing variant of parallel tempering replica exchange Monte Carlo (REMS) that allows for efficient configurational sampling of flexible solutes under the conditions of molecular hydration. Its usage to calculate the thermal stability of a model globular protein, Trp cage TC5b, achieves excellent agreement with experimental measurements. We find that the stability of TC5b is attained through the coupled formation of local and non-local interactions. Remarkably, many of these structures persist at high temperature, concomitant with the origin of native-like configurations and mesostates in an otherwise macroscopically disordered unfolded state. Graph manifold learning reveals that the conversion of these mesostates to the native state is structurally heterogeneous, and that the cooperativity of their formation is encoded largely by the unfolded state ensemble. In all, these studies establish the extent of thermodynamic and structural pre-organization of folding of this model globular protein, and achieve the calculation of macromolecular stability ab initio, as required for ab initio structure prediction, genome annotation, and drug design