Search

We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.

The oxidative cyclization of 1,5-dienes by metal oxo species is a powerful method for stereocontrolled synthesis of tetrahydrofuran diols (THF-diols), structural motifs present in many bioactive natural products. Oxidative cyclization of (2E,6E)-octa-2,6-diene catalyzed by OsO4/NMO as has been studied using density functional theory (DFT) calculations (M06-2X/aug-cc-pVDZ/Hay-Wadt VDZ (n+1) ECP), highlighting the remarkable effect of acid on the fate of the first intermediate, an Os(VI) dioxoglycolate. Strong acid promotes cyclization of the Os(VI) dioxoglycolate, or its NMO complex, through protonation of an oxo ligand to give more electrophilic species. By contrast, in absence of acid reoxidation may occur to afford the Os(VIII) trioxoglycolate, which is shown to favor conventional “second cycle” dihydroxylation reactivity rather than cyclization. The results of the calculations are consistent with experimental results for reactions of OsO4/NMO with 1,5-dienes with acid (oxidative cyclization) and without acid (second cycle osmylation / dihydroxylation). Detailed evaluation of potential catalytic cycles support oxidation of the cyclized Os(IV) THF diolate intermediate to the corresponding Os(VI) species followed by slow hydrolysis, and finally regeneration of OsO4.

The charge distribution of NO2 groups within the crystalline polymorphs of energetic materials strongly affects their explosive properties. We use the recently introduced basis-space iterated stockholder atom partitioning of high-quality charge distributions to examine the approximations that can be made in modeling polymorphs and their physical properties, using 1,3,5-trinitroperhydro-1,3,5-triazine, trinitrotoluene, 1-3-5-trinitrobenzene, and hexanitrobenzene as exemplars. The NO2 charge distribution is strongly affected by the neighboring atoms, the rest of the molecules, and also significantly by the NO2 torsion angle within the possible variations found in observed crystal structures. Thus, the proposed correlations between the molecular electrostatic properties, such as trigger-bond potential or maxima in the electrostatic potential, and impact sensitivity will be affected by the changes in conformation that occur on crystallization. We establish the relationship between the NO2 torsion angle and the likelihood of occurrence in observed crystal structures, the conformational energy, and the charge and dipole magnitude on each atom, and how this varies with the neighboring groups. We examine the effect of analytically rotating the atomic multipole moments to model changes in torsion angle and establish that this is a viable approach for crystal structures but is not accurate enough to model the relative lattice energies. This establishes the basis of transferability of the NO2 charge distribution for realistic nonempirical model intermolecular potentials for simulating energetic materials.

First-principles quantum mechanical calculations with methods such as density functional theory (DFT) allow the accurate calculation of interaction energies between molecules. These interaction energies can be dissected into chemically relevant components such as electrostatics, polarization, and charge transfer using energy decomposition analysis (EDA) approaches. Typically EDA has been used to study interactions between small molecules; however, it has great potential to be applied to large biomolecular assemblies such as protein–protein and protein–ligand interactions. We present an application of EDA calculations to the study of ligands that bind to the thrombin protein, using the ONETEP program for linear-scaling DFT calculations. Our approach goes beyond simply providing the components of the interaction energy; we are also able to provide visual representations of the changes in density that happen as a result of polarization and charge transfer, thus pinpointing the functional groups between the ligand and protein that participate in each kind of interaction. We also demonstrate with this approach that we can focus on studying parts (fragments) of ligands. The method is relatively insensitive to the protocol that is used to prepare the structures, and the results obtained are therefore robust. This is an application to a real protein drug target of a whole new capability where accurate DFT calculations can produce both energetic and visual descriptors of interactions. These descriptors can be used to provide insights for tailoring interactions, as needed for example in drug design.

We report the development and implementation of an energy decomposition analysis (EDA) scheme in the ONETEP linear-scaling electronic structure package. Our approach is hybrid as it combines the localized molecular orbital EDA (Su, P.; Li, H. J. Chem. Phys., 2009, 131, 014102) and the absolutely localized molecular orbital EDA (Khaliullin, R. Z.; et al. J. Phys. Chem. A, 2007, 111, 8753–8765) to partition the intermolecular interaction energy into chemically distinct components (electrostatic, exchange, correlation, Pauli repulsion, polarization, and charge transfer). Limitations shared in EDA approaches such as the issue of basis set dependence in polarization and charge transfer are discussed, and a remedy to this problem is proposed that exploits the strictly localized property of the ONETEP orbitals. Our method is validated on a range of complexes with interactions relevant to drug design. We demonstrate the capabilities for large-scale calculations with our approach on complexes of thrombin with an inhibitor comprised of up to 4975 atoms. Given the capability of ONETEP for large-scale calculations, such as on entire proteins, we expect that our EDA scheme can be applied in a large range of biomolecular problems, especially in the context of drug design.

The partitioning of the energy in ab initio quantum mechanical calculations into its chemical origins (e.g., electrostatics, exchange-repulsion, polarization, and charge transfer) is a relatively recent development; such concepts of isolating chemically meaningful energy components from the interaction energy have been demonstrated by variational and perturbation based energy decomposition analysis approaches. The variational methods are typically derived from the early energy decomposition analysis of Morokuma [Morokuma, J. Chem. Phys., 1971, 55, 1236], and the perturbation approaches from the popular symmetry-adapted perturbation theory scheme [Jeziorski et al., Methods and Techniques in Computational Chemistry: METECC-94, 1993, ch. 13, p. 79]. Since these early works, many developments have taken place aiming to overcome limitations of the original schemes and provide more chemical significance to the energy components, which are not uniquely defined. In this review, after a brief overview of the origins of these methods we examine the theory behind the currently popular variational and perturbation based methods from the point of view of biochemical applications. We also compare and discuss the chemical relevance of energy components produced by these methods on six test sets that comprise model systems that display interactions typical of biomolecules (such as hydrogen bonding and pi-pi stacking interactions) including various treatments of the dispersion energy.