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.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (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.

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.