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.

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.

Post-depositional modification of atmospheric dust on a granite building in central Rio de Janeiro: implications for surface induration and subsequent stone decay

Extensive contour scaling of a 200 year old granite church is associated with the breaching of an apparently iron-rich crust and the widespread deposition of atmospheric dust within the canyon-like streetscape of Rio de Janeiro. Contemporary dust, accumulated dust from within the a depression on the building surface, the surface crust and the underlying granite are examined by a combination of total element analysis and sequential extraction, X-ray diffraction and energy dispersive X-ray fluorescence. Results indicate an increase in total organic carbon and a marked decrease in pH within the accumulated dust, and a rapid mobilization of anions and cations from the water-soluble and carbonate phases. It is considered that the latter is linked to salt accumulation within and eventual salt weathering of the granite. Post-depositional alteration of the dust is also linked with the de-silicification of clay minerals (Illite to kaolinite) and the loss of silica from the amorphous Fe/Mn phase of the accumulated dust under the initially saline and progressively more acidic conditions experienced at the stone - atmosphere interface. This mobilization of silica is associated with the formation of what is, in effect, a thin silica-rich surface crust or glaze. Within the glaze, assessory amounts of extractable iron are concentrated within the amorphous and crystalline Fe/Mn phases at levels that are significantly elevated with respect to the underlying granite, but much lower than the equivalent phases of the accumulated dust from which it is principally assumed to derive. The protection afforded to the stone work by the crust is not, however, permanent and within the last 15 years it has been possible to observe a rapid increase in the surface delamination of the church close to street level.