The Dalton quantum chemistry program system.

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

Insights into the dynamics of evaporation and proton migration in protonated water clusters from large-scale Born-Oppenheimer direct dynamics.

Large-scale on-the-fly Born-Oppenheimer molecular dynamics simulations using recent advances in linear scaling electronic structure theory and trajectory integration techniques have been performed for protonated water clusters around the magic number (H(2)O)(n)H(+) , for n = 20 and 21. Besides demonstrating the feasibility and efficiency of the computational approach, the calculations reveal interesting dynamical details. Elimination of water molecules is found to be fast for both cluster sizes but rather insensitive to the initial geometry. The water molecules released acquire velocities compatible with thermal energies. The proton solvation shell changes between the well-known Eigen and Zundel motifs and is characterized by specific low-frequency vibrational modes, which have been quantified. The proton transfer mechanism largely resembles that of bulk water but one interesting variation was observed.

An efficient density-functional-theory force evaluation for large molecular systems.

An efficient, linear-scaling implementation of Kohn-Sham density-functional theory for the calculation of molecular forces for systems containing hundreds of atoms is presented. The density-fitted Coulomb force contribution is calculated in linear time by combining atomic integral screening with the continuous fast multipole method. For higher efficiency and greater simplicity, the near-field Coulomb force contribution is calculated by expanding the solid-harmonic Gaussian basis functions in Hermite rather than Cartesian Gaussians. The efficiency and linear complexity of the molecular-force evaluation is demonstrated by sample calculations and applied to the geometry optimization of a few selected large systems.

Synthesis, structure, and antimycobacterial activity of 6-[1(3H)-isobenzofuranylidenemethyl]purines and analogs.

6-Benzofuryl-, styryl, benzyl, and furfurylpurines as well as 6-[1(3H)-isobenzofuranylidenemethyl]purines have been synthesized and their activities against Mycobacterium tuberculosis (Mtb) determined. Several compounds displayed profound antimycobacterial activity in combination with low toxicity towards mammalian cells. NMR and X-ray crystallography were employed to determine the detailed structures and the results were supported by quantum chemical calculations.

Atomistic and electronic structure of bimetallic cobalt/rhenium clusters from density functional theory calculations.

We have carried out computational density functional investigations of Co I Re J (J=0,1,2; I+J=14) metal atom clusters. Through thorough optimization of geometry, spin polarization, and electronic configuration, the most stable structures for each cluster have been identified. While the global minima are found to be well defined and energetically well separated from other local minima, the study reveals a plethora of different structures and symmetries only moderately higher in energy. A key point of interest is the effect of doping the cobalt clusters with rhenium. Aside from significant structural reorganizations, rhenium is found to stabilize the clusters and couple down the spin. Furthermore, the most stable clusters comprise highly coordinated rhenium and, in the case of Co 12 Re 2, Re-Re bonding. Our results are compared to earlier experimental and computational data.

Models of fragmentations induced by electron attachment to protonated peptides.

Invoking a number of theoretical levels ranging from HF/STO-3G to CCSD(T)/aug-cc-pVQZ, we have made a detailed survey of six potential energy surfaces (NH4+, NH4*, [CH3CONHCH3]H+, [CH3CONHCH3]H*, [HCONHCH2CONH2]H+, [HCONHCH2CONH2]H*). In conjunction with this, ab inito direct dynamics calculations have been conducted, tracing out several hundred reaction trajectories to reveal details of the electron-capture dissociation mechanism. The model calculations suggest the possibility of a bimodal pattern where some of the radicals, formed upon recombination, dissociate almost directly within one picosecond, the remaining radicals being subject to complete energy redistribution with a subsequent dissociation occurring at the microsecond timescale. Both processes give rise to c and z backbone fragments, resulting from cleavage of N-Calpha bonds of the peptide chain.

Decomposition of protonated formic acid: one transition state--two product channels.

The unimolecular chemistry of protonated formic acid, [HCOOH]H(+), has been investigated by analyzing the fragmentation of metastable ions (MI) during flight in a sector mass spectrometer, and by proton transfer to formic acid in a Fourier-transform ion cyclotron resonance (FT-ICR) mass spectrometer. High level ab initio calculations have been used to model the relevant parts of the potential energy surface (PES). In addition, ab initio direct dynamics calculations have been conducted, tracing out 60 different reaction trajectories. The only stable isomer in the mass spectrometric experiments is HC(OH)(2)(+), which is the precursor to both observed ionic products, HCO(+) and H(3)O(+), via the same saddle point of the potential energy surface. The detailed motion of the dissociating molecule during passage of the post-transition state region of the PES therefore determines which product ion is formed. After passing the TS a transient HC(O)OH(2)(+) molecule is first formed. High total energy increases the probability that the nascent water molecule will have sufficient speed to escape the HCO(+) moiety. Otherwise, typically at low energies, the two units recombine, upon which intra-complex proton transfer is very likely. Eventually, this will give the more stable H(3)O(+).

Calculation of electric dipole hypershieldings at the nuclei in the Hellmann-Feynman approximation.

The third-rank electric hypershieldings at the nuclei of four small molecules have been evaluated at the Hartree-Fock level of theory in the Hellmann-Feynman approximation. The nuclear electric hypershieldings are closely related to molecular vibrational absorption intensities and a generalization of the atomic polar tensors (expanded in powers of the electric field strength) is proposed to rationalize these intensities. It is shown that the sum rules for rototranslational invariance and the constraints imposed by the virial theorem provide useful criteria for basis-set completeness and for near Hartree-Fock quality of nuclear shieldings and hypershieldings evaluated in the Hellmann-Feynman approximation. Twelve basis sets of different size and quality have been employed for the water molecule in an extended numerical test on the practicality of the proposed scheme. The best results are obtained with the R12 and R12+ basis sets, designed for the calculation of electronic energies by the explicitly correlated R12 method. The R12 basis set is subsequently used to investigate three other molecules, CO, N2, and NH3, verifying that the R12 basis consistently performs very well.