DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

Semiempirical (SE) methods can be derived from either Hartree–Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.

We present a new computational strategy to evaluate the charge-transfer (CT) parameters for hole transfer in DNA. On the basis of a fragment-orbital approach, site energies and coupling integrals for a coarse-grained tight-binding description of the electronic structure of DNA are rapidly calculated using the approximative density functional method SCC-DFTB. The methodology is validated by extensive test calculations in comparison with DFT and ab initio reference data, demonstrating its high accuracy at low computational cost. Environmental effects are captured using a quantum mechanics-molecular mechanics (QM/MM) coupling scheme, and dynamical effects are included by evaluating the CT parameters along classical molecular dynamics simulations. This combined methodology allows for a realistic treatment of CT processes in DNA.

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.

Charge transfer in DNA has received much attention in the last few years due to its role in oxidative damage and repair in DNA and also due to possible applications of DNA in nanoelectronics. Despite intense experimental and theoretical efforts, the mechanism underlying long-range hole transport is still unresolved. This is in particular due to the sensitive dependence of charge transfer on the complex structure and dynamics of DNA and the interaction with the solvent, which could not be addressed adequately in the modeling approaches up to now. In this work, we study the factors governing hole transfer in detail, using a DFT-based fragment-orbital method, which allows to compute the charge transfer parameters along multinanosecond molecular dynamics simulations. Environmental effects are captured using a hybrid quantum mechanics-molecular mechanics (QM/MM) coupling scheme. This methodology allows to analyze several factors responsible for charge transfer in DNA in detail. The fluctuation of counterions, strongly counterbalanced by the surrounding water, leads to large oscillations of onsite energies, which govern the energetics of hole propagation along the DNA strand. In contrast, the electronic couplings depend only on DNA conformation and are not affected by the solvent. In particular, the onsite energies are strongly correlated between neighboring nucleobases, indicating that a conformational-gating type of mechanism may be induced by the collective environmental degrees of freedom.

Altogether eight keto and enol tautomers of guanine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in bulk water. The structures of isolated, as well as mono- and dihydrated tautomers were determined by means of the RI-MP2 method using the extended TZVPP (5s3p2d1f/3s2p1d) basis set. The relative energies of isolated tautomers included the correction to higher correlation energy terms evaluated at the CCSD(T)/aug-cc-pVDZ level. The relative enthalpies at 0 K and relative free energies at 298 K were based on the above-mentioned relative energies and zero-point vibration energies, temperature-dependent enthalpy terms and entropies evaluated at the MP2/6-31G level. The keto form having hydrogen atom at N7 is the global minimum while the canonical form having hydrogen atom at N9 represents the first local minimum at all theoretical levels in vacuo and in the presence of 1 and 2 water molecules. All three unusual rare tautomers having hydrogens at N3 and N7, at N3 and N9, and also at N9 and N7 are systematically considerably less stable and can be hardly detected in the gas phase. The theoretical predictions fully agree with existing theoretical as well as experimental results. The effect of bulk solvent on the relative stability of guanine tautomers was studied by self-consistent reaction field and molecular dynamics free energy calculations using the thermodynamic integration method. Bulk solvent, surprisingly, strongly favored these three rare tautomers over all remaining low-energy tautomers and probably only these forms can exist in water phase. The global minimum (tautomer with hydrogens at N3 and N7) is by 13 kcal/mol more stable than the canonical form (3rd local minimum). Addition of one or two water molecules does not change the relative stability order of isolated guanine tautomers but the respective trend clearly supports the surprising stabilization of three rare forms.

We present a coarse-grained tight-binding method based on density functional theory (DFT) for the simulation of charge transfer in complex materials. The charge-transfer parameters are computed using a fragment-orbital approach combined with the approximative DFT method self-consistent charge density functional tight binding (SCC-DFTB), which allows to follow the dynamics of excess charge along nanosecond MD trajectories, still accounting for the important impact of structural fluctuations and solvent effects. Since DFT suffers from the self-interaction error, which would lead to a delocalization of the hole charge over the entire system, we study the effect of an empirical self-interaction correction in detail. The wave function of the excess charge is propagated within the framework of time-dependent DFT, where the electron (hole) and the atomic system are propagated simultaneously according to the derived coupled equations of motion. In the case of DNA, the solvent polarization effects are a dominant factor affecting the hole transport. The hole charge polarizes the surrounding water, which in turn supports a localization of the hole charge--a water polaron is formed, extended dynamically over several nucleobases. As this polarization of water accompanies the migrating hole, the motion of hole is significantly slowed down due to the solvent reorganization energy involved. The estimated hopping rate between neighboring adenines in poly(A)-DNA is in the order of 100 ns(-1), and our simulations clearly show that the charge transfer occurs in a nonadiabatic fashion, due to the small average electronic coupling of around 0.06 eV.

The photoreceptor protein cryptochrome is thought to host, upon light absorption, a radical pair that is sensitive to very weak magnetic fields, endowing migratory birds with a magnetic compass sense. The molecular mechanism that leads to formation of a stabilized, magnetic field sensitive radical pair has despite various theoretical and experimental efforts not been unambiguously identified yet. We challenge this unambiguity through a unique quantum mechanical molecular dynamics approach where we perform electron transfer dynamics simulations taking into account the motion of the protein upon the electron transfer. This approach allows us to follow the time evolution of the electron transfer in an unbiased fashion and to reveal the molecular driving force that ensures fast electron transfer in cryptochrome guaranteeing formation of a persistent radical pair suitable for magnetoreception. We argue that this unraveled molecular mechanism is a general principle inherent to all proteins of the cryptochrome/photolyase family and that cryptochromes are, therefore, tailored to potentially function as efficient chemical magnetoreceptors.