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.

Application of static charge transfer within an ionic-liquid force field and its effect on structure and dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3-dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are +/- 0.5 e. While the three-dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly-with the positive result that for ionic charges of +/- 0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial-charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic-liquid models in which the ionic charges are assumed to be +/- 1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled-charge sets may be of benefit in the simulation of these systems-especially when looking at properties beyond liquid structure-thus providing on alternative to computationally expensive polarisable force fields.

Thermodynamic properties and atomistic structure of the dry amorphous silica surface from a reactive force field model

A force field model of the Keating type supplemented by rules to break, form, and interchange bonds is applied to investigate thermodynamic and structural properties of the amorphous SiO2 surface. A simulated quench from the liquid phase has been carried out for a silica sample made of 3888 silicon and 7776 oxygen atoms arranged on a slab similar to 40 angstrom thick, periodically repeated along two directions. The quench results into an amorphous sample, exposing two parallel square surfaces of similar to 42 nm(2) area each. Thermal averages computed during the quench allow us to determine the surface thermodynamic properties as a function of temperature. The surface tension turns out to be gamma=310 +/- 20 erg/cm(2) at room temperature and gamma=270 +/- 30 at T=2000 K, in fair agreement with available experimental estimates. The entropy contribution Ts-s to the surface tension is relatively low at all temperatures, representing at most similar to 20% of the surface energy. Almost without exceptions, Si atoms are fourfold coordinated and oxygen atoms are twofold coordinated. Twofold and threefold rings appear only at low concentration and are preferentially found in proximity of the surface. Above the glass temperature T-g=1660 +/- 50 K, the mobility of surface atoms is, as expected, slightly higher than that of bulk atoms. The computation of the height-height correlation function shows that the silica surface is rough in the equilibrium and undercooled liquid phase, becoming smooth below the glass temperature T-g.

Statistical-mechanical description of classical test-particle dynamics in the presence of an external force field: modelling noise and damping from first principles

Aiming to establish a rigorous link between macroscopic random motion (described e.g. by Langevin-type theories) and microscopic dynamics, we have undertaken a kinetic-theoretical study of the dynamics of a classical test-particle weakly coupled to a large heat-bath in thermal equilibrium. Both subsystems are subject to an external force field. From the (time-non-local) generalized master equation a Fokker-Planck-type equation follows as a "quasi-Markovian" approximation. The kinetic operator thus defined is shown to be ill-defined; in specific, it does not preserve the positivity of the test-particle distribution function f(x, v; t). Adopting an alternative approach, previously introduced for quantum open systems, is proposed to lead to a correct kinetic operator, which yields all the expected properties. A set of explicit expressions for the diffusion and drift coefficients are obtained, allowing for modelling macroscopic diffusion and dynamical friction phenomena, in terms of an external field and intrinsic physical parameters.

Field configured assembly: Programmed manipulation and self-assembly at the mesoscale

Field configured assembly is a programmable force field method that permits rapid, "hands-free" manipulation, assembly, and integration of mesoscale objects and devices. In this method, electric fields, configured by specific addressing of receptor and counter electrode sites pre-patterned at a silicon chip substrate, drive the field assisted transport, positioning, and localization of mesoscale devices at selected receptor locations. Using this approach, we demonstrate field configured deterministic and stochastic self-assembly of model mesoscale devices, i.e., 50 mum diameter, 670 nm emitting GaAs-based light emitting diodes, at targeted receptor sites on a silicon chip. The versatility of the field configured assembly method suggests that it is applicable to self-assembly of a wide variety of functionally integrated nanoscale and mesoscale systems.

Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.

The interference effects of non-rotated versus counter-rotated trials in visuomotor adaptation

An isometric torque-production task was used to investigate interference and retention in adaptation to multiple visuomotor environments. Subjects produced isometric flexion-extension and pronation-supination elbow torques to move a cursor to acquire targets as quickly as possible. Adaptation to a 30 degrees counter-clockwise (CCW) rotation (task A), was followed by a period of rest (control), trials with no rotation (task B0), or trials with a 60 degrees clockwise (CW) rotation (task B60). For all groups, retention of task A was assessed 5 h later. With initial training, all groups reduced the angular deviation of cursor paths early in the movements, indicating feedforward adaptation. For the control group, performance at commencement of the retest was significantly better than that at the beginning of the initial learning. For the B0 group, performance in the retest of task A was not dissimilar to that at the start of the initial learning, while for the B60 group retest performance in task A was markedly worse than initially observed. Our results indicate that close juxtaposition of two visuomotor environments precludes improved retest performance in the initial environment. Data for the B60 group, specifically larger angular errors upon retest compared with initial exposures, are consistent with the presence of anterograde interference. Furthermore, full interference occurred even when the visuomotor environment encountered in the second task was not rotated (B0). This latter novel result differs from those obtained for force field learning, where interference does not occur when task B does not impose perturbing forces, i.e., when B consists of a null field (Brashers-Krug et al., Nature 382:252-255, 1996). The results are consistent with recent proposals suggesting different interference mechanisms for visuomotor (kinematic) compared to force field (dynamic) adaptations, and have implications for the use of washout trials when studying interference between multiple visuomotor environments.

Interaction of Room Temperature Ionic Liquid Solutions with a Cholesterol Bilayer

Water solutions of representative (IC(4)mim][Cl] and [C(4)mim][Tf2N] room temperature ionic liquids (ILs) in contact with a neutral lipid bilayer made of cholesterol molecules has been investigated by molecular dynamics simulations based on an empirical force field model. The results show that both ILs display selective adsorption at the water-cholesterol interface, with partial inclusion of ions into the bilayer. In the case Of [C(4)mim][Cl], the adsorption of ions at the water-cholesterol interface is limited by a sizable bulk solubility of the IL, driven by the high water affinity of [Cl](-). The relatively low Solubility Of [C(4)mim][Tf2N], instead, gives rise to a nearly complete segregation of the IL component on the bilayer, altering its volume, compressibility, and electrostatic environment. The computational results display important similarities to the results of recent experimental measurements for ILs in contact with phospholipid model membranes (see Evans, K. O. Int. J. Mol. Sci. 2008, 9, 498-511 and references therein).

Structural motifs of cholesterol nanoparticles

The growth sequence of gas-phase cholesterol clusters (Ch(N)) with up to N=36 molecules has been investigated by atomistic simulation based on an empirical force field model. The results of long annealings from high temperature show that the geometric motifs characterizing the structure of pure cholesterol crystals already appear in nanometric aggregates. In all clusters molecules tend to align along a common direction. For cluster sizes above the smallest ones, dispersion interactions among the hydrocarbon body and tails of cholesterol cooperate with hydrogen bonding to give rise to a bilayer structure. Analysis of snapshots from the annealing shows that the condensation of hydrogen bonds into a connected network of rings and chains is an important step in the self-organization of cholesterol clusters. The effect of solvation on the equilibrium properties of medium-size aggregates is investigated by short molecular dynamics simulations for the N=30 and N=40 clusters in water at near ambient conditions and in supercritical carbon dioxide at T=400 K.

Suggested Binding Mechanism of the HIV-gp120 to its CD4 Receptor

The molecular recognition and attachment of the CD4 molecule and the HIV envelope glycoprotein (gp120) might be described as a consecutive three-step molecular recognition process. 1. (a) Long range interaction: electrostatic pre-orientation, 2. (b) short range interaction: electronic attachment followed by a ‘Locking-in’ (via aromatic ring orientation) and 3. (c) internal interaction (induced fit): conformational readjustment of the protein molecules. On the basis of the preliminary investigations (X-ray structures of CD4 and biological studies of CD4 and gp120 point mutants) we described a computational model. This approach consists of empirical calculations as well as ab initio level of quantum chemistry. The conformational analysis of the wild type and mutant CD4 molecules was supported by molecular mechanics and dynamics (Amber force field). The latter analysis involves the application of a novel method, the Amino Acid Conformation Assignment of Proteins (ACAP) software, developed for the notation of secondary protein structures. According to the cardinal role of the electrostatic factors during this interaction, several ab initio investigations were performed for better understanding of the recognition process on submolecular level. Using the above mentioned computational model, we could interpret the basic behaviours and predict some additional features of CD4-gp120 interaction, in spite of the missing gp120 X-ray structure.

Microscopic theory for random processes in weakly-coupled open systems

In order to relate macroscopic random motion (described e.g. by Langevin-type theories) to microscopic dynamics, we have undertaken the derivation of a Fokker-Planck-type equation from first microscopic principles. Both subsystems are subject to an external force field. Explicit expressions for the diffusion and drift coefficients are obtained, in terms of the field.

A Neutron Diffraction and Molecular Dynamics Investigation of Acetate-based Ionic Liquids as Solvents for Glucose

The liquid state structure of the ionic liquid, 1-ethyl-3-methylimidazolium acetate, and the solute/solvent structure of glucose dissolved in the ionic liquid at a 1: 6 molar ratio have been investigated at 323 K by molecular dynamics simulations and neutron diffraction experiments using H/D isotopically substituted materials. Interactions between hydrogen-bond donating cation sites and polar, directional hydrogen-bond accepting acetate anions are examined. Ion-ion radial distribution functions for the neat ionic liquid, calculated from both MD and derived from the empirical potential structure refinement model to the experimental data, show the alternating shell-structure of anions around the cation, as anticipated. Spatial probability distributions reveal the main anion-to-cation features as in-plane interactions of anions with imidazolium ring hydrogens and cation-cation planar stacking. Interestingly, the presence of the polarised hydrogen-bond acceptor anion leads to increased anion-anion tail-tail structuring within each anion shell, indicating the onset of hydrophobic regions within the anion regions of the liquid.

One-dimensional particle simulation of the filamentation instability: Electrostatic field driven by the magnetic pressure gradient force

Two counterpropagating cool and equally dense electron beams are modeled with particle-in-cell simulations. The electron beam filamentation instability is examined in one spatial dimension, which is an approximation for a quasiplanar filament boundary. It is confirmed that the force on the electrons imposed by the electrostatic field, which develops during the nonlinear stage of the instability, oscillates around a mean value that equals the magnetic pressure gradient force. The forces acting on the electrons due to the electrostatic and the magnetic field have a similar strength. The electrostatic field reduces the confining force close to the stable equilibrium of each filament and increases it farther away, limiting the peak density. The confining time-averaged total potential permits an overlap of current filaments with an opposite flow direction.