Michael Gaus

Parametrization and Benchmark of DFTB3 for Organic Molecules.

DFTB3 is a recent extension of the self-consistent-charge density-functional tight-binding method (SCC-DFTB) and derived from a third order expansion of the density functional theory (DFT) total energy around a given reference density. Being applied in combination with the parametrization of its predecessor (MIO), DFTB3 improves for hydrogen binding energies, proton affinities, and hydrogen transfer barriers. In the present study, parameters especially designed for DFTB3 are presented, and its performance is evaluated for small organic molecules focusing on thermochemistry, geometries, and vibrational frequencies from our own and several databases from literature. The new parameters remove significant overbinding errors, reduce errors for geometries of noncovalent interactions, and improve the overall performance.

Parameterization of DFTB3/3OB for Sulfur and Phosphorus for Chemical and Biological Applications

We report the parametrization of the approximate density functional tight binding method, DFTB3, for sulfur and phosphorus. The parametrization is done in a framework consistent with our previous 3OB set established for O, N, C, and H, thus the resulting parameters can be used to describe a broad set of organic and biologically relevant molecules. The 3d orbitals are included in the parametrization, and the electronic parameters are chosen to minimize errors in the atomization energies. The parameters are tested using a fairly diverse set of molecules of biological relevance, focusing on the geometries, reaction energies, proton affinities, and hydrogen bonding interactions of these molecules; vibrational frequencies are also examined, although less systematically. The results of DFTB3/3OB are compared to those from DFT (B3LYP and PBE), ab initio (MP2, G3B3), and several popular semiempirical methods (PM6 and PDDG), as well as predictions of DFTB3 with the older parametrization (the MIO set). In general, DFTB3/3OB is a major improvement over the previous parametrization (DFTB3/MIO), and for the majority cases tested here, it also outperforms PM6 and PDDG, especially for structural properties, vibrational frequencies, hydrogen bonding interactions, and proton affinities. For reaction energies, DFTB3/3OB exhibits major improvement over DFTB3/MIO, due mainly to significant reduction of errors in atomization energies; compared to PM6 and PDDG, DFTB3/3OB also generally performs better, although the magnitude of improvement is more modest. Compared to high-level calculations, DFTB3/3OB is most successful at predicting geometries; larger errors are found in the energies, although the results can be greatly improved by computing single point energies at a high level with DFTB3 geometries. There are several remaining issues with the DFTB3/3OB approach, most notably its difficulty in describing phosphate hydrolysis reactions involving a change in the coordination number of the phosphorus, for which a specific parametrization (3OB/OPhyd) is developed as a temporary solution; this suggests that the current DFTB3 methodology has limited transferability for complex phosphorus chemistry at the level of accuracy required for detailed mechanistic investigations. Therefore, fundamental improvements in the DFTB3 methodology are needed for a reliable method that describes phosphorus chemistry without ad hoc parameters. Nevertheless, DFTB3/3OB is expected to be a competitive QM method in QM/MM calculations for studying phosphorus/sulfur chemistry in condensed phase systems, especially as a low-level method that drives the sampling in a dual-level QM/MM framework.

Density functional tight binding: application to organic and biological molecules

In this work, we review recent extensions of the density functional tight binding (DFTB) methodology and its application to organic and biological molecules. DFTB denotes a class of computational models derived from density functional theory (DFT) using a Taylor expansion around a reference density. The first‐ and second‐order models, DFTB1 and DFTB2, have been reviewed recently (WIREs Comput Mol Sci 2012, 2:456–465). Here, we discuss the extension to third order, DFTB3, which in combination with a modification of the Coulomb interactions in the second‐order formalism and a new parametrization scheme leads to a significant improvement of the overall performance. The performance of DFTB2 and DFTB3 for organic and biological molecules are discussed in detail, as well as problems and limitations of the underlying approximations. WIREs Comput Mol Sci 2014, 4:49–61. doi: 10.1002/wcms.1156

Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems

We present an extension to the recent 3OB parametrization of the Density Functional Tight Binding Model DFTB31,2 for biological and organic systems. Parameters for the halogens F, Cl, Br, and I have been developed for use in covalently bound systems and benchmarked on a test set of 106 molecules (the ‘OrgX’ set), using bonding distances, bonding angles, atomization energies, and vibrational frequencies to assess the performance of the parameters. Additional testing has been done with the X40 set of 40 supramolecular systems containing halogens,3 adding a simple correction for the halogen bonds that are strongly overbound in DFTB3. Furthermore, parameters for Ca, K, and Na as counterions in biological systems have been created. To benchmark geometries as well as ligand binding energies a test set ‘BioMe’ of 210 molecules has been created that cover coordination to various functional groups frequently occurring in biological systems. The new DFTB3/3OB parameter set outperforms DFT calculations with a double-ζ basis set in terms of energies and can reproduce DFT geometries, with some minor deviations in bond distances and angles due to the use of a minimal basis set.

Parametrization of DFTB3/3OB for Magnesium and Zinc for Chemical and Biological Applications

We report the parametrization of the approximate density functional theory, DFTB3, for magnesium and zinc for chemical and biological applications. The parametrization strategy follows that established in previous work that parametrized several key main group elements (O, N, C, H, P, and S). This 3OB set of parameters can thus be used to study many chemical and biochemical systems. The parameters are benchmarked using both gas-phase and condensed-phase systems. The gas-phase results are compared to DFT (mostly B3LYP), ab initio (MP2 and G3B3), and PM6, as well as to a previous DFTB parametrization (MIO). The results indicate that DFTB3/3OB is particularly successful at predicting structures, including rather complex dinuclear metalloenzyme active sites, while being semiquantitative (with a typical mean absolute deviation (MAD) of ∼3–5 kcal/mol) for energetics. Single-point calculations with high-level quantum mechanics (QM) methods generally lead to very satisfying (a typical MAD of ∼1 kcal/mol) energetic properties. DFTB3/MM simulations for solution and two enzyme systems also lead to encouraging structural and energetic properties in comparison to available experimental data. The remaining limitations of DFTB3, such as the treatment of interaction between metal ions and highly charged/polarizable ligands, are also discussed.

Parametrization of the SCC-DFTB Method for Halogens

Parametrization of the approximative DFT method SCC-DFTB for halogen elements is presented. The new parameter set is intended to describe halogenated organic as well as inorganic molecules, and it is compatible with the established parametrization of SCC-DFTB for carbon, hydrogen, oxygen, and nitrogen. The performance of the parameter set is tested on a representative set of molecules and discussed.

Extended Polarization in Third-Order SCC-DFTB from Chemical-Potential Equalization

In this work, we augment the approximate density functional method SCC-DFTB (DFTB3) with the chemical-potential equalization (CPE) approach in order to improve the performance for molecular electronic polarizabilities. The CPE method, originally implemented for the NDDO type of methods by Giese and York, has been shown to significantly emend minimal basis methods with respect to the response properties and has been applied to SCC-DFTB recently. CPE allows this inherent limitation of minimal basis methods to be overcome by supplying an additional response density. The systematic underestimation is thereby corrected quantitatively without the need to extend the atomic orbital basis (i.e., without increasing the overall computational cost significantly). The dependency of polarizability as a function of the molecular charge state, especially, was significantly improved from the CPE extension of DFTB3. The empirical parameters introduced by the CPE approach were optimized for 172 organic molecules in order to match the results from density functional theory methods using large basis sets. However, the first-order derivatives of molecular polarizabilities (e.g., required to compute Raman activities) are not improved by the current CPE implementation (i.e., Raman spectra are not improved).

DFTB3 Parametrization for Copper: The Importance of Orbital Angular Momentum Dependence of Hubbard Parameters

We report the parametrization of a density functional tight binding method (DFTB3) for copper in a spin-polarized formulation. The parametrization is consistent with the framework of 3OB for main group elements (ONCHPS) and can be readily used for biological applications that involve copper proteins/peptides. The key to our parametrization is to introduce orbital angular momentum dependence of the Hubbard parameter and its charge derivative, thus allowing the 3d and 4s orbitals to adopt different sizes and responses to the change of charge state. The parametrization has been tested by applying to a fairly broad set of molecules of biological relevance, and the properties of interest include optimized geometries, ligand binding energies, and ligand proton affinities. Compared to the reference QM level (B3LYP/aug-cc-pVTZ, which is shown here to be similar to the B97-1 and CCSD(T) results, in terms of many properties of interest for a set of small copper containing molecules), our parametrization generally gives reliable structural properties for both Cu(I) and Cu(II) compounds, although several exceptions are also noted. For energetics, the results are more accurate for neutral ligands than for charged ligands, likely reflecting the minimal basis limitation of DFTB3; the results generally outperform NDDO based methods such as PM6 and even PBE with the 6-31+G(d,p) basis. For all ligand types, single-point B3LYP calculations at DFTB3 geometries give results very close (∼1–2 kcal/mol) to the reference B3LYP values, highlighting the consistency between DFTB3 and B3LYP structures. Possible further developments of the DFTB3 model for a better treatment of transition-metal ions are also discussed. In the current form, our first generation of DFTB3 copper model is expected to be particularly valuable as a method that drives sampling in systems that feature a dynamical copper binding site.

QM/MM simulations of vibrational spectra of bacteriorhodopsin and channelrhodopsin-2

Channelrhodopsin-2 is a light-gated ion channel, which has been studied intensively over the last decade. Vibrational spectroscopic experiments started to shed light on the structural changes, that occur during the photocycle, especially in the hydrogen-bonded network surrounding the protonated D156 and C128 - the DC gate. However, the interpretation of these experiments was only based on homology models. Since then, an X-ray structure and better computational models became available. In this article, we show that in combination with a recent reparametrization, the approximate DFT method, DFTB, is able to describe the effects of hydrogen bonding on the C=O stretch vibration in carboxylic acids reliably and agrees well with full DFT results. We apply DFTB in a QM/MM framework to perform vibrational analysis of buried aspartic acids in bacteriorhodopsin and channelrhodopsin-2. Using this approach, we can simulate the FTIR spectral difference between D115 in the dark-adapted and K states of bacteriorhodopsin. The FTIR experiments on the DC gate in channelrhodopsin-2 are well described using an indirect model, where D156 and C128 are bridged via a water molecule.

Improved Electronic Properties from Third-Order SCC-DFTB with Cost Efficient Post-SCF Extensions

The present work outlines the implementation and performance of two cost efficient post-SCF extensions into the third-order SCC-DFTB code. The first one, the charge model 3 (CM3), corrects for errors in bond dipoles for an improved description of molecular charge distribution as compared to the standard Mulliken partitioning scheme. The second one focuses on the response of the charge density, that is, the electronic molecular polarizability, described inaccurately from SCC-DFTB due to the usage of a minimal atomic orbital basis. Here, a variational approach, based on scaled dipole integrals, was implemented, which clearly outperforms standard finite electric field approaches for polarizability calculations by approximately 1 order of magnitude. Both extensions in the present work rely on a set of empirical parameters, which were fitted against 112 organic molecules to match a reference data set from full density functional calculations with a large basis. As an achievement, notably improved electronic properties, that is, molecular dipole moments and polarizabilities, result from SCC-DFTB calculations at negligible additional computational cost. Furthermore, the accuracy of infrared and Raman intensities was tested as first-order derivatives of the new dipoles and polarizabilities as a function of normal mode vibrations. As a result, the current implementations cannot contribute to an improved prediction of relative intensity pattern from SCC-DFTB as compared to ab initio reference data.