A dielectric continuum molecular dynamics method

We introduce a novel method to simulate hydrated macromolecules with a dielectric continuum representation of the surrounding solvent. In our approach, the interaction between the solvent and the molecular degrees of freedom is described by means of a polarization density free energy functional which is minimum at electrostatic equilibrium. After a pseudospectral expansion of the polarization and a discretization of the functional, we construct the equations of motion for the system based on a Car-Parrinello technique. In the limit of the adiabatic evolution of the polarization field variables, our method provides the solution of the dielectric continuum problem "on the fly," while the molecular coordinates are propagated. In this first study, we show how our dielectric continuum molecular dynamics method can be successfully applied to hydrated biomolecules, with low cost compared to free energy simulations with explicit solvent. To our knowledge, this is the first time that stable and conservative molecular dynamic simulations of solutes can be performed for a dielectric continuum model of the solvent. (C) 2001 American Institute of Physics.

Charge exchange and dissociative processes in collisions of slow He2+ ions with H2O molecules

Experimental and theoretical studies of one-electron capture in collisions of He2+ ions with H2O molecules have been carried out in the range 0.025-12 keV amu(-1) corresponding to typical solar wind velocities of 70-1523 km s(-1). Translational energy spectroscopy (TES), photon emission spectroscopy (PES), and fragment ion spectroscopy were employed to identify and quantify the collision mechanisms involved. Cross sections for selective single electron capture into n=1, 2, and 3 states of the He+ ion were obtained using TES while PES provided cross sections for capture into the He+(2p) and He+(3p) states. Our model calculations show that He+(n=2) and He+(n=3) formation proceeds via a single-electron process governed by the nucleus-electron interaction. In contrast, the He+(1s) formation mechanism involves an exothermic two-electron process driven by the electron-electron interaction, where the potential energy released by the electron capture is used to remove a second electron thereby resulting in fragmentation of the H2O molecule. This process is found to become increasingly important as the collision energy decreases. The experimental cross sections are found to be in reasonable agreement with cross sections calculated using the Demkov and Landau-Zener models.

Collision mechanisms in one-electron capture by slow H and He- like ions of C, N and O in H and H-2

Translational energy spectroscopy (TES) has been used to study state-selective one-electron capture by H and He-like ions of C, N and O in both H and H-2 within the range 250-900 eV amu(- 1). The main collision mechanisms leading to state-selective electron capture have been identified, their relative importance assessed and compared, where possible, with theoretical predictions and with any previous measurements based on photon emission spectroscopy. For one-electron capture in H-2, the relative importance of contributions from non- dissociative and dissociative capture as well as from two- electron capture into autoionizing states is found to be strikingly different for the cases considered. Our TES measurements in atomic hydrogen provide an important extension of previous measurements to energies below 1000 eV amu(-1) and show that, as the impact energy decreases, electron capture becomes more selective until only a single n product channel is significant. These product main channels are well described by reaction windows calculated using a Landau-Zener approach. However, the same approach applied to the more complex energy- change spectra observed in H-2 is found to be less successful.

Ultrafast ionization study of N-2 in intense linearly and circularly polarized laser fields

A detailed investigation has been carried out of N-2 molecules in intense 55 and 220 fs, linear and circular polarized, 790 nm laser pulses. Using an intensity selective scanning technique, ionization, dissociation, and dissociative ionization channels have been studied. Remarkably similar enhancements of signal with linear polarization observed for double ionization and dissociation channels demonstrate the dominance of dynamic alignment over rescattering effects. Fragmentation energies from dissociative ionization are reasonably well reproduced by classical trajectory calculations, the higher charged fragments displaying evidence of post dissociative ionization.

High harmonic generation by halogen anions and noble gas atoms in a laser field

A comparative study of high harmonic generation (HHG) by atoms and ions with active p-electrons is carried out in the theoretical framework of the rescattering mechanism. The substate with m(l) = 0, i.e. zero orbital momentum projection along the electric vector of a linearly polarized laser wave, is found to give the major contribution to the HHG rate. Our calculations for HHG by an H atom in an excited 2p-state demonstrate that the rate for recombination into a final state with a different value of m(l) (= +/- 1), is higher for lower harmonic orders N, while for higher N (beyond the plateau domain) the difference vanishes. For species with closed electron shells, the m(l)-changing transitions are forbidden by the Pauli exclusion principle. We report absolute HHG rates for halogen ions and noble gas atoms at various intensities. These results demonstrate that the Coulomb binding potential of the atoms considerably enhances both the ionization and recombination steps in the rescattering process. However, the weak binding energy of the anions allows lower orders of HHG to be efficiently produced at relatively low intensities, from which we conclude that observation of HHG by an anion is experimentally feasible.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Solvation effects on equilibria: Triazoles and N-methyl piperidinol

We have performed calculations of the solvation effects on a number of equilibrium constants in water using a recently proposed hybrid quantum classical scheme in which the liquid environment is modelled using classical solvent molecules and the solute electronic structure is computed using modern quantum chemical methods. The liquid phase space is sampled from a fully classical simulation. We find that solvation effects on both triazole tautomeric equilibrium constants and piperidinol conformational equilibrium constants can be interpreted in terms of subtle differences in the local environment which can be seen in probability densities and radial distribution functions. Lower level calculations were performed for comparison and we conclude that the solvation thermodynamics can be predicted from a good classical model of solvent and solute molecules, but the implicit models that we tried are less successful.

Concertedness and solvent effects in multiple proton transfer reactions: The formic acid dimer in solution

The issue of multiple proton transfer (PT) reactions in solution is addressed by performing molecular dynamics simulations for a formic acid dimer embedded in a water cluster. The reactant species is treated quantum mechanically, within a density functional approach, while the solvent is represented by a classical model. By constraining different distances within the dimer we analyze the PT process in a variety of situations representative of more complex environments. Free energy profiles are presented, and analyzed in terms of typical solvated configurations extracted from the simulations. A decrease in the PT barrier height upon solvation is rationalized in terms of a transition state which is more polarized than the stable states. The dynamics of the double PT process is studied in a low-barrier case and correlated with solvent polarization fluctuations. Cooperative effects in the motion of the two protons are observed in two different situations: when the solvent polarization does not favor the transfer of one of the two protons and when the motion of the two protons is not synchronized. This body of observations is correlated with local structural and dynamical properties of the solvent in the vicinity of the reactant. (C) 2000 American Institute of Physics. [S0021-9606(00)51121-0].

Molecular electrostatic properties of ions in an ionic liquid

We have analysed the electronic wave functions from an ab initio simulation of the ionic liquid (room temperature molten salt) dimethyl imidazolium chloride ([dmim][Cl] or [C1mim][Cl]) using localized Wannier orbitals. This allows us to assign electron density to individual ions. The probability distributions of the ionic dipole moments for an isolated ion and for ions in solution are compared. The liquid environment is found to polarize the cation by about 0.7 D and to increase the amplitude of the fluctuations in the dipole moments of both cation and anion. The relative changes in nuclear and electronic contributions are shown. The implications for classical force fields are discussed.

Molecular conduction: Do time-dependent simulations tell you more than the Landauer approach?

A dynamical method for simulating steady-state conduction in atomic and molecular wires is presented which is both computationally and conceptually simple. The method is tested by calculating the current-voltage spectrum of a simple diatomic molecular junction, for which the static Landauer approach produces multiple steady-state solutions. The dynamical method quantitatively reproduces the static results and provides information on the stability of the different solutions. (c) 2006 American Institute of Physics.

Glory and thresholds effects in H+D-2 reactive angular scattering

We analyse H + D-2 reactive angular scattering using the S- matrix elements obtained by Aoiz et al. and Althorpe et al. Enhancement of small angle scattering in the v' = 3

Semiclassical complex angular momentum theory and Pade reconstruction for resonances, rainbows, and reaction thresholds

A semiclassical complex angular momentum theory, used to analyze atom-diatom reactive angular distributions, is applied to several well-known potential (one-particle) problems. Examples include resonance scattering, rainbow scattering, and the Eckart threshold model. Pade reconstruction of the corresponding matrix elements from the values at physical (integral) angular momenta and properties of the Pade approximants are discussed in detail.

Simultaneous histories, path sums, and measurements for noncommuting variables

We analyze von Neumann-like quantum measurements in terms of simultaneous virtual paths constructed for two noncommuting variables. The approach is applied to measurements of operator functions of conjugate variables and to the joint measurements of such variables. The limits of applicability of the restricted phase space path integral are studied. We demonstrate that, for a simple joint measurement, using entangled meter states allows one to manipulate the order in which the measurements are conducted. The effects of '' weakening '' a measurement by choosing unsharp meter states are also discussed.

Superfluid toroidal currents in atomic condensates

The spectrum of collective excitations of oblate toroidal condensates within the Bogoliubov approximation was studied, and the dynamical stability of ring currents around the torus explored. The transition from spheroidal to toroidal geometry of the trap displaced the energy levels into narrow bands. A simple, but accurate, formula was detailed for the lowest angular acoustic modes of excitation, and the splitting energy when a background current is present.