Abstract:Molecular dynamics simulations for liquid CaAl2Si2O8 have been carried out at 72 state points spanning ranges in density (ρ: 2398–4327kg/m3), temperature (T: 3490–6100K) and pressure (P: 0.84–120GPa) relevant to geosystems. The atomic scale structure of the melt is determined by analysis of nearest neighbor coordination statistics as a function of T and P. Dramatic structural change occurs as pressure increases especially for 0< P <20GPa at all temperatures. Changes in structure are encapsulated by examining the coordination of Si, Al, Ca and O around oxygen and vice versa. Si and Al change from predominantly fourfold at low-P to dominantly sixfold for P >∼20GPa. Pentahedrally coordinated Si and Al in distorted trigonal bipyramids attain abundance maxima corresponding to ∼60% of total (Si, Al)O n at 3–5GPa and weakly depend on T. The coordination of Ca by oxygen increases from 7 to 10 for 0< P <20GPa and changes slowly for P >20GPa at 3500K. Similar behavior is seen at 6000K except that the interval of rapid changes occurs at higher pressure. Oxygen with only one nearest Si or Al neighbor (i.e., non-bridging oxygen, NBO) decreases whereas oxygen with two or three nearest neighbors of Si, Al or Ca increases as pressure increases. Changes in melt structure are reflected in the variation of thermodynamic and transport properties of the liquid. Values of the self-diffusivities of Ca, Al, Si and O are fit to a modified Arrhenian expression and compare well to limited laboratory data. Self-diffusivities are best fit using ‘low P’ and ‘high-P’ expressions, identical in form but with different parameters, with activation energies and activation volumes in the range 150–200kJ/mol and +5 to −1cm3/mol, respectively. Green–Kubo calculations for liquid shear viscosity are presented and compare well with limited laboratory results. Application of the Eyring model to determine the characteristic size and number of atoms in the activated cluster based on independently computed D and η suggests that the activated cluster decreases from ∼8 to ∼3 atoms from low to high pressure while its characteristic size shrinks from ∼14Å to ∼3Å providing insight into dynamics of atom mobility and possible cooperative behavior. The equation of state and variation of internal energy with T and V are used in Part II (Ghiorso et al., 2009) to derive a comprehensive thermodynamic description of liquid CaAl2Si2O8. This is best accomplished by allowing for EOS expressions broken into high and low pressure intervals consistent with coordination statistics and MD-derived transport properties.

The Structure and Transport Property of Liquid Al with Different EAM Model

Analysis of semi-empirical interatomic potentials appropriate for simulation of crystalline and liquid Al and Cu

Structural and Thermodynamic Properties of Liquid Transition Metals with Different Embedded-Atom Method Models

The Solid–liquid Interface Free Energy of Al: A Comparison Between Molecular Dynamics Calculations and Experimental Measurements

Molecular Dynamics Computation of the Liquid Structure of Fe50Al50 Alloy

Structural and mechanical properties of Fe–Al compounds: An atomistic study by EAM simulation

Structural and Mechanical Properties of Homogeneous Solid-Liquid Interface of Al Modelled with COMB3 Potential

Ab Initiomolecular Dynamics Simulation of Liquidalxge1−xalloys

A Theoretical and Computational Study of Structure‐Property Relationship in Five Binary Liquid Alloys

Comparative Studies Of Role Of Different Structure Factors In Various Physical Properties Of Some Simple Liquid Metals

Composition- and Temperature-Dependent Liquid Structures in Al–Cu Alloys: Anab Initiomolecular Dynamics and X-Ray Diffraction Study

Correlation Between Viscosity and Local Atomic Structure in Liquid Zr56Co28Al16 Alloy

Structure, thermodynamic and transport properties of CaAl2Si2O8 liquid. Part I: Molecular dynamics simulations

Experimental Andab Initiomolecular Dynamics Simulation Studies of Liquidal60cu40alloy

Liquid Li Structure and Dynamics: A Comparison Between OFDFT and Second Nearest‐neighbor Embedded‐atom Method

Evolution of Atomic Structure in Al75cu25 Liquid from Experimental and Ab Initio Molecular Dynamics Simulation Studies

Self-diffusion in the Intermetallic Compounds NiAl and Ni3Al: an Embedded Atom Method Study

The viscosity re-assessment of the Cu–X (X = Ag, Al, Sn, Mg) and Ag–X (X = Sn, Sb, In, Au) liquid alloys

First-principles study of crystal structure prediction, electronic, thermodynamic and mechanical properties of Al-Li binary system

MEAM interatomic potential for thermodynamic and mechanical properties of lithium allotropes

Transport Properties of Undercooled Liquid Copper: A Molecular Dynamics Study