Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Accurate computational methods for two-electron atom-laser interactions

We discuss the application of quantitatively accurate computational methods to the study of laser-driven two-electron atoms in short intense laser pulses. The fundamental importance of such calculations to the subject area is emphasized. Calculations of single- and double-electron ionization rates at 390 nm are presented. (C) 2001 Optical Society of America.

Molecular dynamics simulation of ionic liquids: The effect of electronic polarizability

An electronically polarizable model has been developed for the ionic liquid 1-ethyl-3-methylimidazolium nitrate (EMIM+/NO3-), Molecular dynamics simulation studies were then performed on both the polarizable and nonpolarizable versions of the model. Comparisons of shear viscosity and diffusion constants at 400 K show that the effects of polarizability are quite substantial and the polarizable model results are in better agreement with the experimental values.

A reactive force field simulation of liquid-liquid phase transitions in phosphorus

A force field model of phosphorus has been developed based on density functional (DF) computations and experimental results, covering low energy forms of local tetrahedral symmetry and more compact (simple cubic) structures that arise with increasing pressure. Rules tailored to DF data for the addition, deletion, and exchange of covalent bonds allow the system to adapt the bonding configuration to the thermodynamic state. Monte Carlo simulations in the N-P-T ensemble show that the molecular (P-4) liquid phase, stable at low pressure P and relatively low temperature T, transforms to a polymeric (gel) state on increasing either P or T. These phase changes are observed in recent experiments at similar thermodynamic conditions, as shown by the close agreement of computed and measured structure factors in the molecular and polymer phases. The polymeric phase obtained by increasing pressure has a dominant simple cubic character, while the polymer obtained by raising T at moderate pressure is tetrahedral. Comparison with DF results suggests that the latter is a semiconductor, while the cubic form is metallic. The simulations show that the T-induced polymerization is due to the entropy of the configuration of covalent bonds, as in the polymerization transition in sulfur. The transition observed with increasing P is the continuation at high T of the black P to arsenic (A17) structure observed in the solid state, and also corresponds to a semiconductor to metal transition. (C) 2004 American Institute of Physics.

High-pressure deformation mechanism in scolecite: A combined computational-experimental study

Pressure-induced structural modifications in scolecite were studied by means of in situ synchrotron X-ray powder diffraction and density functional computations. The experimental cell parameters were refined up to 8.5 GPa. Discontinuities in the slope of the unit-cell parameters vs. pressure dependence were observed; as a consequence, an increase in the slope of the linear pressure-volume dependence is observed at about 6 GPa, suggesting an enhanced compressibility at higher pressures. Weakening and broadening of the diffraction peaks reveals increasing structural disorder with pressure, preventing refinement of the lattice parameters above 8.5 GPa. Diffraction patterns collected during decompression show that the disorder is irreversible. Atomic coordinates within unit cells of different dimensions were determined by means of Car-Parrinello simulations. The discontinuous rise in compressibility at about 6 GPa is reproduced by the computation, allowing us to attribute it to re-organization of the hydrogen bonding network, with the formation of water dimers. Moreover we found that, with increasing pressure, the tetrahedral chains parallel to c rotate along their elongation axis and display an increasing twisting along a direction perpendicular to c. At the same time, we observed the compression of the channels. We discuss the modification of the Ca polyhedra under pressure, and the increase in coordination number (from 4 to 5) of one of the two Al atoms, resulting from the approach of a water molecule. We speculate that this last transformation triggers the irreversible disordering of the system.

Numerical simulation of the charge balance in an EBIT

The electron beam ions traps (EBITs) are widely used to study highly charged ions (HCIs). In an EBIT, a high energy electron beam collides with atoms and ions to generate HCIs in the trap region. It is important to study the physics in the trap. The atomic processes, such as electron impact ionisation (EI), radiative recombination (RR), dielectronic recombination (DR) and charge exchange (CX), occur in the trap and numerical simulation can give some parameters for design, predict the composition and describe charge state evolution in an EBIT [Phys. Rev. A 43 (199 1) 4861]. We are presently developing a new code, which additionally includes a description of the overlaps between the ion clouds of the various charge-states. It has been written so that it can simulate experiments where various machine parameters (e.g. beam energy and current) can vary throughout the simulation and will be able to use cross- sections either based on scaling laws or derived from atomic structure calculations. An object-oriented method is used in developing the new software, which is an efficient way to organize and write code. (C) 2003 Elsevier Science B.V. All rights reserved.