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.

Recent extensions of the DMol3 local orbital density functional method for band structure calculations of insulating and metallic solids are described. Furthermore the method for calculating semilocal pseudopotential matrix elements and basis functions are detailed together with other unpublished parts of the methodology pertaining to gradient functionals and local orbital basis sets. The method is applied to calculations of the enthalpy of formation of a set of molecules and solids. We find that the present numerical localized basis sets yield improved results as compared to previous results for the same functionals. Enthalpies for the formation of H, N, O, F, Cl, and C, Si, S atoms from the thermodynamic reference states are calculated at the same level of theory. It is found that the performance in predicting molecular enthalpies of formation is markedly improved for the Perdew–Burke–Ernzerhof [Phys. Rev. Lett. 77, 3865 (1996)] functional.

From molecules to solids with the DMol3 approach

We present the theoretical and technical foundations of the Amsterdam Density Functional (ADF) program with a survey of the characteristics of the code (numerical integration, density fitting for the Coulomb potential, and STO basis functions). Recent developments enhance the efficiency of ADF (e.g., parallelization, near order-N scaling, QM/MM) and its functionality (e.g., NMR chemical shifts, COSMO solvent effects, ZORA relativistic method, excitation energies, frequency-dependent (hyper)polarizabilities, atomic VDD charges). In the Applications section we discuss the physical model of the electronic structure and the chemical bond, i.e., the Kohn–Sham molecular orbital (MO) theory, and illustrate the power of the Kohn–Sham MO model in conjunction with the ADF-typical fragment approach to quantitatively understand and predict chemical phenomena. We review the “Activation-strain TS interaction” (ATS) model of chemical reactivity as a conceptual framework for understanding how activation barriers of various types of (competing) reaction mechanisms arise and how they may be controlled, for example, in organic chemistry or homogeneous catalysis. Finally, we include a brief discussion of exemplary applications in the field of biochemistry (structure and bonding of DNA) and of time-dependent density functional theory (TDDFT) to indicate how this development further reinforces the ADF tools for the analysis of chemical phenomena. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 931–967, 2001

Chemistry with ADF

https://www.zora.uzh.ch/id/eprint/3175/3/PP_VandeVondele_CPCom_2005QuickManVO.pdf

QUICKSTEP: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach

Two kinds of conjugated C3-symmetric perylene dyes, namely, triperylene hexaimides (TPH) and selenium-annulated triperylene hexaimides (TPH-Se), are efficiently synthesized. Both TPH and TPH-Se have broad and strong absorption in the region 300–600 nm together with suitable LUMO levels of about −3.8 eV. Single-crystal X-ray diffraction studies show that TPH displays an extremely twisted three-bladed propeller configuration and a unique 3D network assembly in which three PBI subunits in one TPH molecule have strong π–π intermolecular interactions with PBI subunits in neighboring molecules. The integration of selenophene to TPH endows TPH-Se with a more distorted propeller configuration and a more compact 3D network assembly due to the Se···O interactions. A single-crystal transistor confirms that both TPH and TPH-Se possess good electron-transport ability. TPH and TPH-Se acceptor-based solar cells show high power conversion efficiency of 8.28% and 9.28%, respectively, which mainly results from the combined...

Three-Bladed Rylene Propellers with Three-Dimensional Network Assembly for Organic Electronics.

New generalized gradient approximation (GGA) functionals are reported, using the expansion form of A. D. Becke, J. Chem. Phys. 107, 8554 (1997), with 15 linear parameters. Our original such GGA functional, called HCTH, was determined through a least squares refinement to data of 93 systems. Here, the data are extended to 120 systems and 147 systems, introducing electron and proton affinities, and weakly bound dimers to give the new functionals HCTH/120 and HCTH/147. HCTH/120 has already been shown to give high quality predictions for weakly bound systems. The functionals are applied in a comparative study of the addition reaction of water to formaldehyde and sulfur trioxide, respectively. Furthermore, the performance of the HCTH/120 functional in Car–Parrinello molecular dynamics simulations of liquid water is encouraging.

New generalized gradient approximation functionals

Quantum Simulations of complex many-body systems:from theory to algorithms

An extension of Car-Parrinello (CP) molecular dynamics for efficient treatment of electronically nonadiabatic processes is presented. The current approach couples the ${S}_{1}$ restricted open-shell Kohn-Sham excited state to the ${S}_{0}$ ground state using a surface hopping scheme. Efficient evaluation of the nonadiabatic couplings is achieved by exploiting the available wave function time derivatives. Since the computational cost scales linearly with the number of excited states, the technique makes possible nonadiabatic ab initio simulations of systems of similar complexity to those typically studied by standard CP methods. It is thus ideally suited to study the photochemistry of large molecules, particularly in condensed phases.

Nonadiabatic Car-Parrinello molecular dynamics.

A simple and high yield synthesis of water-soluble arylazopyrazoles (AAPs) featuring superior photophysical properties is reported. The introduction of a carboxylic acid allows the diverse functionalization of AAPs. Based on structural modifications of the switching unit the photophysical properties of the AAPs could be tuned to obtain molecular switches with favorable photostationary states. Furthermore, AAPs form stable and light-responsive host-guest complexes with β-cyclodextrin (β-CD). Our most efficient AAP shows binding affinities comparable to azobenzenes, but more effective switching and higher thermal stability of the Z-isomer. As a proof-of-principle, we investigated two CD-based supramolecular systems, containing either cyclodextrin vesicles (CDVs) or cyclodextrin-functionalized gold nanoparticles (CDAuNPs), which revealed excellent reversible, light-responsive aggregation and dispersion behavior. To conclude, AAPs have great potential to be incorporated as molecular switches in highly demanding and multivalent photoresponsive systems.

Nikos L. Doltsinis

Papers

Three-Bladed Rylene Propellers with Three-Dimensional Network Assembly for Organic Electronics.

New generalized gradient approximation functionals

Quantum Simulations of complex many-body systems:from theory to algorithms

Nonadiabatic Car-Parrinello molecular dynamics.

Arylazopyrazoles as Light-Responsive Molecular Switches in Cyclodextrin-Based Supramolecular Systems