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.

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.

Beyond Ehrenfest: correlated non-adiabatic molecular dynamics

A method for introducing correlations between electrons and ions that is computationally affordable is described. The central assumption is that the ionic wavefunctions are narrow, which makes possible a moment expansion for the full density matrix. To make the problem tractable we reduce the remaining many-electron problem to a single-electron problem by performing a trace over all electronic degrees of freedom except one. This introduces both one- and two-electron quantities into the equations of motion. Quantities depending on more than one electron are removed by making a Hartree-Fock approximation. Using the first-moment approximation, we perform a number of tight binding simulations of the effect of an electric current on a mobile atom. The classical contribution to the ionic kinetic energy exhibits cooling and is independent of the bias. The quantum contribution exhibits strong heating, with the heating rate proportional to the bias. However, increased scattering of electrons with increasing ionic kinetic energy is not observed. This effect requires the introduction of the second moment.

Conventional free radical polymerization in room temperature ionic liquids: a green approach to commodity polymers with practical advantages

Free-radical polymerization of methyl methacrylate and styrene using conventional organic initiators in the room temperature ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate ([ C(4)mim][PF6]) is rapid and produces polymers with molecular weights up to 10x higher than from benzene; both polymerization and isolation of products were achieved without using VOCs, offering economic as well as environmental advantages.

Operation of a free-electron laser from the extreme ultraviolet to the water window

We report results on the performance of a free-electron laser operating at a wavelength of 13.7 nm where unprecedented peak and average powers for a coherent extreme-ultraviolet radiation source have been measured. In the saturation regime, the peak energy approached 170 J for individual pulses, and the average energy per pulse reached 70 J. The pulse duration was in the region of 10 fs, and peak powers of 10 GW were achieved. At a pulse repetition frequency of 700 pulses per second, the average extreme-ultraviolet power reached 20 mW. The output beam also contained a significant contribution from odd harmonics of approximately 0.6% and 0.03% for the 3rd (4.6 nm) and the 5th (2.75 nm) harmonics, respectively. At 2.75 nm the 5th harmonic of the radiation reaches deep into the water window, a wavelength range that is crucially important for the investigation of biological samples.

On the recollision-free excitation of krypton during ultrafast multi-electron tunnel ionization

The probability of multiple ionization of krypton by 50 femtosecond circularly polarized laser pulses, independent of the optical focal geometry, has been obtained for the first time. The excellent agreement over the intensity range 100 TW cm-2 to 100 PW cm-2 with the recent predictions of Kornev et al (2003 Phys. Rev. A 68 043414) provides the first experimental confirmation that non-recollisional electronic excitation can occur in strong-field ionization. This is particularly true for higher stages of ionization, when the laser intensity exceeds 10 PW cm-2 as the energetic departure of the ionized electron(s) diabatically distorts the wavefunctions of the bound electrons. By scaling the probability of ionization by the focal volume, we discuss why this mechanism was not apparent in previous studies.

Soft x-ray free electron laser microfocus for exploring matter under extreme conditions

We have focused a beam (BL3) of FLASH (Free-electron LASer in Hamburg: lambda = 13.5 nm, pulse length 15 fs, pulse energy 10-40 mu J, 5Hz) using a fine polished off-axis parabola having a focal length of 270 mm and coated with a Mo/Si multilayer with an initial reflectivity of 67% at 13.5 nm. The OAP was mounted and aligned with a picomotor controlled six-axis gimbal. Beam imprints on poly(methyl methacrylate) -PMMA were used to measure focus and the focused beam was used to create isochoric heating of various slab targets. Results show the focal spot has a diameter of

The Y-procedure: How to extract the chemical transformation rate from reaction-diffusion data with no assumption on the kinetic model

The paper presents a new method to extract the chemical transformation rate from reaction–diffusion data with no assumption on the kinetic model (“kinetic model-free procedure”). It is a new non-steady-state kinetic characterization procedure for heterogeneous catalysts. The mathematical foundation of the Y-procedure is a Laplace-domain analysis of the two inert zones in a TZTR followed by transposition to the Fourier domain. When combined with time discretization and filtering the Y-procedure leads to an efficient practical method for reconstructing the concentration and reaction rate in the active zone. Using the Y-procedure the concentration and reaction rate of a non-steady state catalytic process can be determined without any pre-assumption regarding the type of kinetic dependence. The Y-procedure is the basis for advanced software for non-steady state kinetic data interpretation. The Y-procedure can be used to relate changes in the catalytic reaction rate and kinetic parameters to changes in the surface composition (storage) of a catalyst.

Dissociative Recombination of Protonated Formic Acid: Implications for Molecular Cloud and Cometary Chemistry

At the heavy ion storage ring CRYRING in Stockholm, Sweden, we have investigated the dissociative recombination of DCOOD2+ at low relative kinetic energies, from ~1 meV to 1 eV. The thermal rate coefficient has been found to follow the expression k(T) = 8.43 × 10-7 (T/300)^-0.78 cm3 s-1 for electron temperatures, T, ranging from ~10 to ~1000 K. The branching fractions of the reaction have been studied at ~2 meV relative kinetic energy. It has been found that ~87% of the reactions involve breaking a bond between heavy atoms. In only 13% of the reactions do the heavy atoms remain in the same product fragment. This puts limits on the gas-phase production of formic acid, observed in both molecular clouds and cometary comae. Using the experimental results in chemical models of the dark cloud, TMC-1, and using the latest release of the UMIST Database for Astrochemistry improves the agreement with observations for the abundance of formic acid. Our results also strengthen the assumption that formic acid is a component of cometary ices.