To predict phase transitions in ceramics and minerals from molecular dynamics simulations, we have developed a force field in which the charges are allowed to readjust instantaneously to the atomic configurations. These charges are calculated using the charge equilibration (QEq) method. In addition to electrostatics, a two-body Morse stretch potential is included to account for shortrange nonelectrostatic interactions. This MS-Q potential is applied herein to SiO 2 , where we find that it describes well the fourfold coordinated and sixfold coordinated systems (such as quartz and stishovite), silica glass, and the pressure-induced phase transition from quartz to stishovite.[ S0031-9007(99) In principle, the phase transitions in minerals and ceramics can be predicted from first principles quantum mechanics (QM), and significant progress is being made along this line [1][2][3]. However, most QM studies have considered only static conditions since dynamical QM simulations are usually not practical for long times (at least nanoseconds) and large system sizes of interest for ceramics. Thus, we need to use classical molecular dynamics (MD) for predicting such phenomena. The problem here is that standard approaches to force fields (FF) [4][5][6][7] for ceramics and oxides use simplifications (fixed charges, three-body potentials) or valence terms [8] which may not be appropriate for describing phase transitions, where the ligancy and structure may change dramatically. For this reason, we have developed an alternative procedure more likely to describe phase transitions of ceramics. Herein we outline the methodology and report the first results for SiO 2 systems.Because electrostatics play an essential role in determining the structure and properties of ceramics, we consider that the first priority of the FF is to produce plausible charges. Since the charges may depend on the distances, angles, and ligancy, we consider that the charge must be allowed to readjust to the instantaneous configuration of the atoms. To do this we use the charge equilibration (QEq) procedure developed by Rappé and Goddard [9] and used in many applications on organic and inorganic systems [10]. Rather than keeping the charges fixed during the dynamics, as in previous calculations, we now allow the charges to adjust to the instantaneous geometric configuration of the atoms.In QEq the charges are determined by requiring that the chemical potential x A be equal on all atoms, and x A is a function of the charges on all of the atoms of the system:Here, the atomic parameters x Thus, J AB describes simple Coulomb for large separations, but it is shielded for short distances. The QEq parameters used for Si and O are given in Table Ib [9].The nonelectrostatic interactions (short-range Pauli repulsion, covalency, dispersion, etc.) are included via a simple two-body Morse-Stretch (MS) term,The MS parameters for Si-O, O-O, and Si-Si, were optimized to describe the properties (density, cohesive energy,