There is, I think, something ethereal about i —the square root of minus one. I remember first hearing about it at school. It seemed an odd beast at that time—an intruder hovering on the edge of reality.

Usually familiarity dulls this sense of the bizarre, but in the case of i it was the reverse: over the years the sense of its surreal nature intensified. It seemed that it was impossible to write mathematics that described the real world in …

I and i

The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.

http://dns2.asia.edu.tw/~ysho/YSHO-English/1000%20WC/PDF/J%20Che%20Phy132,%20154104.pdf

A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu

Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

Semiempirical GGA-type density functional constructed with a long-range dispersion correction.

6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

/pdf/quantum-mechanical-continuum-solvation-models-3edv963g14.pdf

Quantum mechanical continuum solvation models.

Kalte Ionenspektroskopie zur Lösung der Gasphasenstruktur eines Decapeptids

Ab initio mol. dynamics (AIMD) and combined Hybrid/AIMD simulations appear as promising candidates for an in situ modeling of enzymic reactions. In order to probe the capabilities of these methods for the characterization of enzymic processes we have chosen the enzyme Human Carbonic Anhydrase II (HCAII) as a test case system. Several models of the active site have been studied in this work including ab initio cluster models of different sizes, ranging from about 30-90 atoms. We have also extended these quantum mech. (QM) models by taking into account the electrostatic effects of the surrounding protein through a hybrid scheme where the external field is described via charges from a classical MD force field. By comparing the structural and electronic properties of these different models we can probe the importance of size effects of the QM part as well as the influence of the protein environment. For the largest size cluster model we have been able to directly observe part of the enzymic reaction cycle, namely the initial proton transfer steps from the zinc-bound water towards the Nd-atom of the proton acceptor group His 64. [on SciFinder (R)]

Ab initio and hybrid molecular dynamics simulations of the active site of human carbonic anhydrase II: a test case study

The implementation of local control theory using nonadiabatic molecular dynamics within the framework of linear-response time-dependent density functional theory is discussed. The method is applied to study the photoexcitation of lithium fluoride, for which we demonstrate that this approach can efficiently generate a pulse, on-the-fly, able to control the population transfer between two selected electronic states. Analysis of the computed control pulse yields insights into the photophysics of the process identifying the relevant frequencies associated to the curvature of the initial and final state potential energy curves and their energy differences. The limitations inherent to the use of the trajectory surface hopping approach are also discussed.

/pdf/local-control-theory-using-trajectory-surface-hopping-and-xjg9men0m8.pdf

Local Control Theory using Trajectory Surface Hopping and Linear-Response Time-Dependent Density Functional Theory

From a week to less than a day: Speedup and scaling of coordinate-scaled exact exchange calculations in plane waves

Conversion of adenosine triphosphate (ATP) to the second messenger cyclic adenosine monophosphate (cAMP) is an essential reaction mechanism that takes place in eukaryotes, triggering a variety of signal transduction pathways. ATP conversion is catalyzed by the enzyme adenylyl cyclase (AC), which can be regulated by binding inhibitory, Gαi, and stimulatory, Gαs subunits. In the past twenty years, several crystal structures of AC in isolated form and complexed to Gαs subunits have been resolved. Nevertheless, the molecular basis of the inhibition mechanism of AC, induced by Gαi, is still far from being fully understood. Here, classical molecular dynamics simulations of the isolated holo AC protein type 5 and the holo binary complex AC5:Gαi have been analyzed to investigate the conformational impact of Gαi association on ATP-bound AC5. The results show that Gαi appears to inhibit the activity of AC5 by preventing the formation of a reactive ATP conformation.

/pdf/gai1-inhibition-mechanism-of-atp-bound-adenylyl-cyclase-type-3j13ckn7u2.pdf

Ursula Rothlisberger

Papers

Kalte Ionenspektroskopie zur Lösung der Gasphasenstruktur eines Decapeptids

Ab initio and hybrid molecular dynamics simulations of the active site of human carbonic anhydrase II: a test case study

Local Control Theory using Trajectory Surface Hopping and Linear-Response Time-Dependent Density Functional Theory

From a week to less than a day: Speedup and scaling of coordinate-scaled exact exchange calculations in plane waves

Gαi1 inhibition mechanism of ATP-bound adenylyl cyclase type 5.