Thermodynamics data of Alkali Feldspars from FPMD simulations

<p>We computed the thermodynamic properties (pressure, temperature, internal energy, heat capacity) and thermoelastic coefficients (isobaric expansivity, isothermal compressibility, thermal pressure coefficient) on&nbsp;the two alkali feldspars end-members&nbsp;using&nbsp;<em>ab initio</em> molecular dynamics simulation in the&nbsp;2000-20000&nbsp;K temperature range and in the 0.5-6 g.cm^-3 density range.</p> <p>Simulations are performed using the&nbsp;Vienna Ab Initio Simulation Package (VASP) (Kresse and Furthmuller, 1996) in the canonical (NVT) ensemble with a timestep of 0.5-2 fs for 5-20 ps depending on the temperature and density. We model the feldspar end-members in a cubic cell containing 208 atoms (16 formula units) and 1024 or 1152 electrons for the Na- and K-feldspars respectively. For simulations at low density we used pseudopotentials which require a lower plane wave energy cutoff, set to 370 eV.&nbsp;For Na-end-member, we also used&nbsp;&nbsp;hard pseudopotentials at high density in order to reduce the overlap of electronic spheres, in particular for Na-Na pairs. The energy cutoff for this set of pseudopotentials is&nbsp;950 eV.&nbsp;Additional details can be found in the manuscript.</p> <p>&nbsp;</p> <p>There is one dataset for each different composition and set of pseudopotentials used:</p> <ol> <li>kobsch-ds01.txt&nbsp; &nbsp; --&gt;&nbsp; &nbsp; NaAlSi₃O₈</li> <li>kobsch-ds02.txt&nbsp; &nbsp; --&gt;&nbsp; &nbsp; NaAlSi₃O_8&nbsp;and set of pseudopotentials with a lower plane wave energy cutoff than in 01</li> <li>kobsch-ds03.txt&nbsp; &nbsp; --&gt;&nbsp; &nbsp; NaAlSi₃O₈&nbsp;and harder pseudopotentials than in 01</li> <li>kobsch-ds04.txt&nbsp; &nbsp; --&gt;&nbsp; &nbsp; KAlSi₃O₈</li> <li>kobsch-ds05.txt&nbsp; &nbsp; --&gt;&nbsp; &nbsp; KAlSi₃O₈&nbsp;and set of pseudopotentials with a lower plane wave energy cutoff than in 04</li> </ol> <p>&nbsp;</p> <p>Each file present the arithmetic time averages of the pressure (P), temperature (T) and internal energy (E). The standard deviation of the data to the mean is indicated by stdev_X, where X is P, T or E. The statistical error to the mean (err_X) is computed using the blocking method as described by Flyvbjerg and Petersen (1989). The sign &#39;&gt;&#39; is indicated before the value of the statistical error when no convergence &nbsp;was reached during the estimation of this error. The heat capacity Cv is computed using fluctuations on both potential and kinetic energies (Allen and Tildesley, 1989) and its statistical error stdev_Cv&nbsp;is computed using the bootstrap method.&nbsp;</p> <p>We computed the thermoelastic coefficients only for densities (<span class="math-tex">\(\rho\)</span>) above 1.5 g.cm^-3. The thermal pressure coefficient (TPC = <span class="math-tex">\(\frac{\partial P}{\partial T}\big|_V\)</span>) is &nbsp;the slope of linear fit of P vs. T isochores. The isothermal compressibility (<span class="math-tex">\(\beta = -\frac{1}{\rho} \frac{\partial \rho}{\partial P}\big|_T\)</span>) is computed using central finite differences on our P vs. <span class="math-tex">\(\rho\)</span>&nbsp;isotherms. The isobaric expansivity (<span class="math-tex">\(\alpha = \frac{1}{\rho} \frac{\partial \rho}{\partial T}\big|_P\)</span>) is computed using the previously computed <span class="math-tex">\(\beta\)</span>&nbsp;and TPC.</p>