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

: The unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio. Harmonic force fields are obtained using Density Functional Theory (DFT), MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set. DFT calculations use the Local Spin Density Approximation (LSDA), BLYP, and Becke3LYP (B3LYP) density functionals. Mid-IR spectra predicted using LSDA, BLYP, and B3LYP force fields are of significantly different quality, the B3LYP force field yielding spectra in clearly superior, and overall excellent, agreement with experiment. The MP2 force field yields spectra in slightly worse agreement with experiment than the B3LYP force field. The SCF force field yields spectra in poor agreement with experiment.The basis set dependence of B3LYP force fields is also explored: the 6-31G* and TZ2P basis sets give very similar results while the 3-21G basis set yields spectra in substantially worse agreements with experiment. jg

Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields

In undergraduate level organometallic chemistry courses students are usually taught that organometallic compounds are toxic and unstable in air and water. While this is true of many complexes, some are also non-toxic and stable in air and water. Indeed, bioorganometallic chemistry, the study of biomolecules or biologically active molecules that contain at least one carbon directly bound to a metal, is a thriving subject, and air and water stability is a general pre-requisite. This interdisciplinary field is located at the borderline between chemistry, biochemistry, biology and medicine. In this tutorial review, various aspects of bioorganometallic chemistry are introduced, with the main emphasis on medicinal organometallic compounds. Organometallic therapeutics for cancer, HIV and malaria and other medicinal applications are described. It is also shown how rational ligand design has led to new improved therapies much in the same way that an organometallic chemist working in catalysis will design new ligands for improved activities.

Bioorganometallic chemistry—from teaching paradigms to medicinal applications

Application of chitosan and its derivatives as adsorbents for dye removal from water and wastewater: a review.

In the past decade, attosecond technology has opened up the investigation of ultrafast electronic processes in atoms, simple molecules, and solids. Here, we report the application of isolated attosecond pulses to prompt ionization of the amino acid phenylalanine and the subsequent detection of ultrafast dynamics on a sub–4.5-femtosecond temporal scale, which is shorter than the vibrational response of the molecule. The ability to initiate and observe such electronic dynamics in polyatomic molecules represents a crucial step forward in attosecond science, which is progressively moving toward the investigation of more and more complex systems.

/pdf/ultrafast-electron-dynamics-in-phenylalanine-initiated-by-4cvrok3xg0.pdf

Ultrafast electron dynamics in phenylalanine initiated by attosecond pulses

Convergence of the multicenter B-spline DFT approach for the continuum

The accurate crystal structure determinations of MeCbl (1), CNCbl·2LiCl (2), and CNCbl·KCl (3), based on synchrotron diffraction data collected at 100 K and using high-quality single crystals, are reported. Refinements gave R1 indices of 0.0834 (1), 0.0434 (2), and 0.0773 (3). The influence of the water of crystallization and ion content on the crystal packing of these and other cobalamins (XCbl) is discussed, and a relationship between the crystal packing and the corrin side chain conformations is presented. An analysis of the bond lengths within the corrin moiety, based on 13 accurate structures with several X groups, shows that the trend of the C−C and C−N distances can be interpreted in terms of electronic and steric factors. The variation in structural, NMR and IR spectroscopic, and electrochemical properties are compared with those of cobaloximes, the B12 model, when X is varied. This comparison indicates that the π-back-donation from metal to the CN axial ligand and the transmission of the trans in...

Similarities and differences between cobalamins and cobaloximes. Accurate structural determination of methylcobalamin and of LiCl- and KCl-containing cyanocobalamins by synchrotron radiation.

We present an adaptive density-guided approach for the construction of Born–Oppenheimer potential energy surfaces (PES) in rectilinear normal coordinates for use in vibrational structure calculations. The procedure uses one-mode densities from vibrational structure calculations for a dynamic sampling of PESs. The implementation of the procedure is described and the accuracy and versatility of the method is tested for a selection of model potentials, water, difluoromethane and pyrimidine. The test calculations illustrate the advantage of local basis sets over harmonic oscillator basis sets in some important aspects of our procedure.

An adaptive density-guided approach for the generation of potential energy surfaces of polyatomic molecules

A new implementation of the vibrational self-consistent field (VSCF) method is presented on the basis of a second quantization formulation. A so-called active terms algorithm is shown to be a significant improvement over a standard implementation reducing the computational effort by one order in the number of degrees of freedom. Various types of screening provide even further reductions in computational scaling and absolute CPU time. VSCF calculations on large polyaromatic hydrocarbon model systems are presented. Further, it is demonstrated that in cases where distant modes are not directly coupled in the Hamiltonian, down to linear scaling of the required CPU time with respect to the number of vibrational modes can be obtained. This is illustrated with calculations on simple model systems with up to 1 million degrees of freedom.

New Formulation and Implementation of Vibrational Self-Consistent Field Theory.

A procedure for the automatic construction of Born-Oppenheimer (BO) potential energy and molecular property surfaces in rectilinear normal coordinates is presented and its suitability and accuracy when combined with vibrational structure calculations are assessed. The procedure relies on a hierarchical n-mode representation of the BO potential energy or molecular property surface, where the n-mode term of the sequence of potentials/molecular properties includes only the couplings between n or less vibrational degrees of freedom. Each n-mode cut of the energy/molecular property surface is first evaluated in a grid of points with ab initio electronic structure methods. The ab initio data are then spline interpolated and a subsequent polynomial fitting provides an analytical semiglobal representation for use in vibrational structure programs. The implementation of the procedure is outlined and the accuracy of the method is tested on water and difluoromethane. Strategies for improving the proposed algorithm are also discussed.

Daniele Toffoli

Papers

Convergence of the multicenter B-spline DFT approach for the continuum

Similarities and differences between cobalamins and cobaloximes. Accurate structural determination of methylcobalamin and of LiCl- and KCl-containing cyanocobalamins by synchrotron radiation.

An adaptive density-guided approach for the generation of potential energy surfaces of polyatomic molecules

New Formulation and Implementation of Vibrational Self-Consistent Field Theory.

Automatic generation of potential energy and property surfaces of polyatomic molecules in normal coordinates.