Vibrational properties of nanometric AB(2) ionic clusters

A broad survey of harmonic dynamics in AB(2) clusters with up to N = 3000 atoms is performed using a simple rigid ion model, with ionic radii selected to give rutile as the ground state structure for the corresponding extended crystal. The vibrational density of states is already close to its bulk counterpart for N similar to 500, with characteristic differences due to surfaces, edges and vertices. Two methods are proposed and tested to map the cluster vibrational states onto the rutile crystal phonons. The net distinction between infrared (IR) active and Raman active modes that exists for bulk rutile becomes more and more blurred as the cluster size is reduced. It is found that, in general, the higher the IR activity of the mode, the more this is affected by the system size. IR active modes are found to spread over a wide frequency range for the finite clusters. Simple models based on either a crude confinement constraint or surface pressure arguments fail to reproduce the results of the calculations. The effects of the stoichiometry and dielectric properties of the surrounding medium on the vibrational properties of the clusters are also investigated.

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.

Density functional and Monte Carlo studies of sulfur. I. Structure and bonding in S-n rings and chains (n=2-18) (vol 118, pg 9257, 2003)

Density functional calculations have been performed for ring isomers of sulfur with up to 18 atoms, and for chains with up to ten atoms. There are many isomers of both types, and the calculations predict the existence of new forms. Larger rings and chains are very flexible, with numerous local energy minima. Apart from a small, but consistent overestimate in the bond lengths, the results reproduce experimental structures where known. Calculations are also performed on the energy surfaces of S8 rings, on the interaction between a pair of such rings, and the reaction between one S8 ring and the triplet diradical S8 chain. The results for potential energies, vibrational frequencies, and reaction mechanisms in sulfur rings and chains provide essential ingredients for Monte Carlo simulations of the liquid–liquid phase transition. The results of these simulations will be presented in Part II.