Nonharmonic transverse vibration of the H-bonded molecules and the THz spectra in ice and water

A nonlinear equation of motion is found for the dimer comprising two charged H2O molecules. The THz dielectric response to nonharmonic vibration of a nonrigid dipole, forming the hydrogen bond (HB), is found in the direction transverse to this bond. An explicit expression is derived for the autocorrelator that governs the spectrum generated by transverse vibration (TV) of such a dipole. This expression is obtained by analytical solution of the truncated set of recurrence equations. The far infrared (FIR) spectra of ice at the temperature - 7 degrees C are calculated. The wideband, in the wavenumber (frequency) v range 0... 100.0 cm(-1), spectra are obtained for liquid water at room temperature and for supercooled water at -5.6 degrees C. All spectra are represented in terms of the complex permittivity epsilon(v) and the absorption coefficient alpha(v). The obtained analytical formula for epsilon comprises the term epsilon(perpendicular to) pertinent to the studied TV mechanism with three additional terms Delta epsilon(q), Delta epsilon(mu), and epsilon(or) arising, respectively, from: elastic harmonic vibration of charged molecules along the H-bond; elastic reorientation of HB permanent dipoles; and rather free libration of permanent dipoles in 'defects' of water/ice structure. The suggested TV-dielectric relaxation mechanism allows us: (a) to remove the THz 'deficit' of loss epsilon" inherent in previous theoretical studies; (b) to explain the THz loss and absorption spectra in supercooled (SC) water; and (c) to describe, in agreement with the experiment, the low- and high-frequency tails of the two bands of ice H2O located in the range 10...300 cm(-1). Specific THz dielectric properties of SC water are ascribed to association of water molecules, revealed in our study by transverse vibration of HB charged molecules. (C) 2006 Published by Elsevier B.V.

Density functional and Monte Carlo studies of sulfur. I. Structure and bonding in S-n rings and chains (n=2-18) (vol 118, pg 9257, 2003)

Density functional calculations have been performed for ring isomers of sulfur with up to 18 atoms, and for chains with up to ten atoms. There are many isomers of both types, and the calculations predict the existence of new forms. Larger rings and chains are very flexible, with numerous local energy minima. Apart from a small, but consistent overestimate in the bond lengths, the results reproduce experimental structures where known. Calculations are also performed on the energy surfaces of S8 rings, on the interaction between a pair of such rings, and the reaction between one S8 ring and the triplet diradical S8 chain. The results for potential energies, vibrational frequencies, and reaction mechanisms in sulfur rings and chains provide essential ingredients for Monte Carlo simulations of the liquid–liquid phase transition. The results of these simulations will be presented in Part II.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Equilibrium polymerization of cyclic carbonate oligomers. III. Chain branching and the gel transition

Ring-opening polymerization of cyclic polycarbonate oligomers, where monofunctional active sites act on difunctional monomers to produce an equilibrium distribution of rings and chains, leads to a "living polymer." Monte Carlo simulations [two-dimensional (2D) and three-dimensional (3D)] of the effects of single [J. Chem. Phys. 115, 3895 (2001)] and multiple active sites [J. Chem. Phys. 116, 7724 (2002)] are extended here to trifunctional active sites that lead to branching. Low concentrations of trifunctional particles c(3) reduce the degree of polymerization significantly in 2D, and higher concentrations (up to 32%) lead to further large changes in the phase diagram. Gel formation is observed at high total density and sizable c(3) as a continuous transition similar to percolation. Polymer and gel are much more stable in 3D than in 2D, and both the total density and the value of c(3) required to produce high molecular weight aggregates are reduced significantly. The degree of polymerization in high-density 3D systems is increased by the addition of trifunctional monomers and reduced slightly at low densities and low c(3). The presence of branching makes equilibrium states more sensitive (in 2D and 3D) to changes in temperature T. The stabilities of polymer and gel are enhanced by increasing T, and-for sufficiently high values of c(3)-there is a reversible polymer-gel transformation at a density-dependent floor temperature. (C) 2002 American Institute of Physics.

Hydrogen bonding and collective proton modes in clusters and periodic layers of squaric acid: A density functional study

Hydrogen bonding in clusters and extended layers of squaric acid molecules has been investigated by density functional computations. Equilibrium geometries, harmonic vibrational frequencies, and energy barriers for proton transfer along hydrogen bonds have been determined using the Car-Parrinello method. The results provide crucial parameters for a first principles modeling of the potential energy surface, and highlight the role of collective modes in the low-energy proton dynamics. The importance of quantum effects in condensed squaric acid systems has been investigated, and shown to be negligible for the lowest-energy collective proton modes. This information provides a quantitative basis for improved atomistic models of the order-disorder and displacive transitions undergone by squaric acid crystals as a function of temperature and pressure. (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.

Symmetry, delocalization, and molecular conductance

Molecules bonded between two metal contacts form the simplest possible molecular devices. Coupled by the molecule, the left and right contact-based states form symmetric and antisymmetric pairs near the Fermi level. We relate the size of the resulting energy splitting DeltaE to the symmetry and degree of delocalization of the coupling molecular orbital. Qualitative trends in molecular conductances are then estimated from the variations in DeltaE. We examine benzenedithiol and other molecules of interest in transport. (C) 2005 American Institute of Physics.

Accurate and efficient non-adiabatic quantum molecular dynamics approach for laser-matter interactions

A non-adiabatic quantum molecular dynamics approach for treating the interaction of matter with intense, short-duration laser pulses is developed. This approach, which is parallelized to run on massively-parallel supercomputers, is shown to be both accurate and efficient. Illustrative results are presented for harmonic generation occurring in diatomic molecules using linearly polarized laser pulses.

Ionization dynamics of laser-driven H-2(+)

We set out aspects of a numerical algorithm used in solving the full-dimensionality time-dependent Schrodinger equation describing the electronic motion of the hydrogen molecular ion driven by an intense, linearly polarized laser pulse aligned along the molecular axis. This algorithm has been implemented within the fixed inter-nuclear separation approximation in a parallel computer code, a brief summary of which is given. Ionization rates are calculated and compared with results from other methods, notably the time-independent Floquet method. Our results compare very favourably with the precise predictions of the Floquet method, although there is some disagreement with other wavepacket calculations. Visualizations of the electron dynamics are also presented in which electron rescattering is observed.