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

다중혈관 관상동맥 환자에서 y-문합을 이용하여 양쪽 내흉동맥만을 사용한 우회술의 조기 성적

It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a standard "zero-damping" formula and rational damping to finite values for small interatomic distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coefficients is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interatomic forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramolecular dispersion in four representative molecular structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermolecular distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of corrected GGAs for non-covalent interactions. According to the thermodynamic benchmarks BJ-damping is more accurate especially for medium-range electron correlation problems and only small and practically insignificant double-counting effects are observed. It seems to provide a physically correct short-range behavior of correlation/dispersion even with unmodified standard functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying density functional.

Effect of the damping function in dispersion corrected density functional theory

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.

TDDFT and RI-CC2 calculations have been performed on the excited-state intramolecular proton transfer in malonaldehyde, o-hydroxybenzaldehyde, salicylic acid, 7-hydroxy-1-indanone, and 2-(2'-hydroxyphenyl)benzothiazole. Vertical and adiabatic excitation energies have been computed for the npi and pipi states. Overall, we have found that both RI-CC2 and TDDFT methods are good candidates for the description of ESIPT potential energy surfaces. Proton transfer (PT) curves have been computed for both excited states. An essentially barrierless and very shallow energy profile has been found for the pi pi* state. For the n pi* state the keto minimum is more pronounced than for the pi pi* state and, depending on the case, energy barriers ranging from values <0.1 eV up to 0.5 eV were found. From the computed PT curves we conclude that extended crossing regions between the two excited states will occur.

Excited-State Intramolecular Proton Transfer: A Survey of TDDFT and RI-CC2 Excited-State Potential Energy Surfaces

In the present study a benchmark set of medium-sized and large aromatic organic molecules with 10–78 atoms is presented. For this test set 0–0 transition energies measured in supersonic jets are compared to those calculated with DFT and the B3LYP functional, ADC(2), CC2 and the spin-scaled CC2 variants SOS-CC2 and SCS-CC2. Geometries of the ground and excited states have been optimized with these methods in polarized triple zeta basis sets. Zero-point vibrational corrections have been calculated with the same methods and basis sets. In addition the energies have been corrected by single point calculations with a triple zeta basis augmented with diffuse functions, aug-cc-pVTZ. The deviations of the theoretical results from experimental electronic origins, which have all been measured in the gas phase with high-resolution techniques, were evaluated. The accuracy of SOS-CC2 is comparable to that of unscaled CC2, whereas ADC(2) has slightly larger errors. The lowest errors were found for SCS-CC2. All correlated wave function methods provide significantly better results than DFT with the B3LYP functional. The effects of the energy corrections from the augmented basis set and the method-consistent calculation of the zero-point vibrational corrections are small. With this benchmark set reliable reference data for 0–0 transition energies for larger organic chromophores are available that can be used to benchmark the accuracy of other quantum chemical methods such as new DFT functionals or semi-empirical methods for excitation energies and structures and thereby augments available benchmark sets augments present benchmark sets which include mainly smaller molecules.

Benchmarks for 0–0 transitions of aromatic organic molecules: DFT/B3LYP, ADC(2), CC2, SOS-CC2 and SCS-CC2 compared to high-resolution gas-phase data

Distributed polarizabilities using the topological theory of atoms in molecules

We present a detailed study of the effects of electron correlation on two-photon absorption calculated by coupled cluster quadratic response theory. The hierarchy of coupled cluster models CCS, CC2, CCSD, and CC3 has been used to investigate the effects of electron correlation on the two-photon absorption cross sections of formaldehyde (CH2O), diacetylene (C4H2), and water (H2O). In particular, the effects of triple excitations on two-photon transition cross sections are determined for the first time. In addition, we present a detailed comparison of the coupled cluster results with those obtained from Hartree-Fock and density-functional response theories. We have investigated the local-density approximation, the pure Becke-Lee-Yang-Parr (BLYP) functional, the hybrid Becke-3-parameter-Lee-Yang-Parr (B3LYP), and the Coulomb-attenuated B3LYP (CAM-B3LYP) functionals. Our results show that the CAM-B3LYP functional, when used in conjuction with a one-particle basis-set containing diffuse functions, has much promise; however, care must still be exercised for diffuse Rydberg-type states.

/pdf/benchmarking-two-photon-absorption-with-cc3-quadratic-jwk70g600n.pdf

Christof Hättig

Papers

Excited-State Intramolecular Proton Transfer: A Survey of TDDFT and RI-CC2 Excited-State Potential Energy Surfaces

Benchmarks for 0–0 transitions of aromatic organic molecules: DFT/B3LYP, ADC(2), CC2, SOS-CC2 and SCS-CC2 compared to high-resolution gas-phase data

Distributed polarizabilities using the topological theory of atoms in molecules

Benchmarking two-photon absorption with CC3 quadratic response theory, and comparison with density-functional response theory.

Frequency-dependent first hyperpolarizabilities using coupled cluster quadratic response theory