Density functional theories with simulation, combination

The main purpose of this chapter is to present the basics of ab initio molecular dynamics, focusing on the practical aspects of the simulations, and in particular, on modeling chemical reactions. Although CP-MD is a general molecular dynamics scheme which potentially can be applied in combination with any electronic structure method, the Car-Parinello MD is usually implemented within the framework of density functional theory with plane-waves as the basis set. Such an approach is conceptually quite distant from the commonly applied static approaches of quantum-chemistry with atom-centered basis sets. Therefore, a main... 

For calculation of the equilibrium compositions of the liquid phase either the equilibrium constants of the dissociation and polycondensation reactions have to be known or they can be computed by methods which use the approach of minimizing Gibbs free energy [200-202]. In addition, ab initio modeling techniques such as density functional theory (DFT) in combination with reactive molecular dynamic (MD) simulations could be used. Once the liquid phase system is modeled, there are in principle two options to describe the vapor-liquid equilibrium. Either equations of state (EOS) or excess Gibbs free energy models (g -models) may be used to describe the thermodynamics of the liquid... 

In 1985 Car and Parrinello invented a method [111-113] in which molecular dynamics (MD) methods are combined with first-principles computations such that the interatomic forces due to the electronic degrees of freedom are computed by density functional theory [114-116] and the statistical properties by the MD method. This method and related ab initio simulations have been successfully applied to carbon [117], silicon [118-120], copper [121], surface reconstruction [122-128], atomic clusters [129-133], molecular crystals [134], the epitaxial growth of metals [135-140], and many other systems for a review see Ref. 113. 

Further studies were carried out with halocarbene amides 34 and 357 Although again no direct spectroscopic signatures for specifically solvated carbenes were found, compelling evidence for such solvation was obtained with a combination of laser flash photolysis (LFP) with UV-VIS detection via pyridine ylides, TRIR spectroscopy, density functional theory (DFT) calculations, and kinetic simulations. Carbenes 34 and 35 were generated by photolysis of indan-based precursors (Scheme 4.7) and were directly observed by TRIR spectroscopy in Freon-113 at 1635 and 1650 cm , respectively. The addition of small amounts of dioxane or THF significantly retarded the rate of biomolecular reaction with both pyridine and TME in Freon-113. Also, the addition of dioxane increased the observed lifetime of carbene 34 in Freon-113. These are both unprecedented observations. 

Ab initio molecular orbital methodology or density functional theory [158-160] would be suited for this combined QM/MM approach. However, in order to be able to compute the QM energies along the Monte Carlo simulation, nowadays a semiempirical Hamiltonian, like AMI [161], is a much more computationally efficient method. Before using AMI, the goodness of the semiempirical results in gas phase in comparison with the ab initio ones has to be tested. For systems in which the semiempirical results are poor, the relation... 

Theoretical calculations and simulations using ah initio and density function theory (DFT) methodologies are also seeing increasing use. Combining these theoretical calculations with spectroscopic data can assist in the interpretation of the observed spectral features and an improved understanding of how a probe molecule interacts with the various types of sites in zeolitic systems. 

The electronic coupling of donor and acceptor sites, connected via a t-stack, can either be treated by carrying out a calculation on the complete system or by employing a divide-and-conquer (DC) strategy. With the Hartree-Fock (HF) method or a method based on density functional theory (DFT), full treatment of a d-a system is feasible for relatively small systems. Whereas such calculations can be performed for models consisting of up to about ten WCPs, they are essentially inaccessible even for dimers when one attempts to combine them with MD simulations. Semiempirical quantum chemical methods require considerably less effort than HF or DFT methods also, one can afford application to larger models. However, standard semiempirical methods, e.g., AMI or PM3, considerably underestimate the electronic couplings between r-stacked donor and acceptor sites and, therefore, a special parameterization has to be invoked (see below). 

In the present chapter, we will focus on the simulation of the dynamics of photoexcited nucleobases, in particular on the investigation of radiationless decay dynamics and the determination of associated characteristic time constants. We use a nonadiabatic extension of ab initio molecular dynamics (AIMD) [15, 18, 21, 22] which is formulated entirely within the framework of density functional theory. This approach couples the restricted open-shell Kohn-Sham (ROKS) [26-28] first singlet excited state, Su to the Kohn-Sham ground state, S0, by means of the surface hopping method [15, 18, 94-97], The current implementation employs a plane-wave basis set in combination with periodic boundary conditions and is therefore ideally suited to condensed phase applications. Hence, in addition to gas phase reference simulations, we will also present nonadiabatic AIMD (na-AIMD) simulations of nucleobases and base pairs in aqueous solution. 

First principles approaches are important as they avoid many of the pitfalls associated with using parameterized descriptions of the interatomic interactions. Additionally, simulation of chemical reactivity, reactions and reaction kinetics really requires electronic structure calculations [108]. However, such calculations were traditionally limited in applicability to rather simplistic models. Developments in density functional theory are now broadening the scope of what is viable. Car-Parrinello first principles molecular dynamics are now being applied to real zeolite models [109,110], and the combined use of classical and quantum mechanical methods allows quantum chemical methods to be applied to cluster models embedded in a simpler description of the zeoUte cluster environment [105,111]. 

In the first part of this work, a brief overview over several strategies to combine such time domain transport simulations with first principles electronic structure theory is given. For the latter, we restrict ourselves to a discussion of time dependent density functional theory (TDDFT) only. This method is by far the most employed many body approach in this field and provides an excellent ratio of accuracy over computational cost, allowing for the treatment of realistic molecular devices. This digest builds on the earlier excellent survey by Koentopp and co-workers on a similar topic [13]. Admittedly and inevitably, the choice of the covered material is biased by the authors interests and background. 

The geometric and electronic structures of the mononuclear [Cu(02)(TpAd,pr)] complex have been evaluated using Cu K- and L-edge X-ray absorption spectroscopy (XAS) studies in combination with valence bond configuration interaction (VBCI) simulations and spin-unrestricted broken symmetry density functional theory (DFT) calculations.205... 

A promising method, developed in recent years, is the use of first principles molecular dynamics as exemplified by the Car-Parrinello technique (8]. In these calculations the interatomic potentials are explicitly derived from the electronic ground-state within the density functional theory in local or non-local approximation. It combines quantum mechanical calculations with molecular dynamics simulations and, therefore, overcomes the limitations of both methods. Actual computers allow only simulations of aqueous solutions of about 60 water molecules for several ps (10 s). This limit is still at least one order of magnitude shorter than the fastest directly measured water exchange rate, k = 3.5 x 10 s for [Eu(H20)8], i.e. one exchange event every (8 x 3.5 x lO s ) = 36 ps [9]. Nevertheless, several publications appeared in the late 1990s on solvated Be [10], K+ [11] and Cu + [12] presenting mainly structural results. 

---