Fitting equations of state User deﬁnes a start guess for mineral properties, the parameters they wish to change the pressure, temperature and volume data and optional uncertainties. Function returns ﬁtted parameters and the associated covariance matrix. The authors are thankful for input from (alphabetically by surname): Wolfgang Bangerth, Zack Geballe, Marc Hirschmann, Yu Huang, Jessica Irving, JiaChao Liu, Valentina Magni, Bill McDonough, Motohiko Murakami, Wendy Panero, Barbara Romanowicz and Quentin Williams. We thank Lars Stixrude for providing benchmarking calculations. This project was initiated at the CIDER workshop in 2012 (www.deepearth.org). The project is funded and supported by: BurnMan: Towards a multidisciplinary toolkit for reproducible deep Earth science Sanne Cottaar1, Robert Myhill2, Timo Heister3, Ian Rose4, Cayman Unterborn5, Juliane Dannberg6 1 University of Cambridge, 2 University of Bristol, 3 Clemson University, 4 University of California, Berkeley, 5 Arizona State University, 6 Colorado State University/Texas A&M University We present BurnMan[1], an open-source mineral physics python toolbox to determine thermodynamic and elastic properties for Earth minerals and rocks. BurnMan currently includes: Several equations of state (including Vinet, Birch Murnaghan, Mie-Debye-Grueneisen, Modiﬁed Tait) Popular solid solution models (ideal, symmetric, asymmetric, subregular) Published databases of minerals (Holland and Powell, 2011 [2]; Stixrude & Lithgow-Bertelloni, 2005; 2011[3]) Different averaging schemes for seismic velocities (Reuss, Voigt, VRH, Hashin-Shtrikman) Reference seismic models- 1D and soon 3D Reference geotherms Repository of work published using BurnMan A wide range of examples, benchmarks and unit tests BurnMan is designed to be simple and intuitive to use. For example, calculating and printing the properties of Mg0.9Fe0.1SiO3-bridgmanite at 30 GPa and 2000 K is just four lines of code: BurnMan is available at www.burnman.org and is being openly developed at https://github.com/geodynamics/burnman. Latest mineral physics additions Summary Interface with seismology Interface with planetary science Interface with geodynamics $ burnman.tools.fit_PVT_data(mineral, fit_params, PT, V, V_sigma) Hugoniot calculations (for shock experiments) User deﬁnes the reference conditions at the start of the shock, the set of pressures along which to calculate the Hugoniot, and an optional reference mineral. Function returns the volumes and temperatures along the Hugoniot. $ burnman.tools.hugoniot(mineral, P_ref, T_ref, pressures, reference_mineral) Equilibrium assemblage calculations (work in progress) User deﬁnes a bulk composition, a set of phases (endmembers or solutions), two constraints (compositional or P-T), and optional compositional guesses. Function updates assemblage with equilibrium compositions and conditions. [right] Olivine phase diagram constructed using BurnMan. $ burnman.equilibriumassemblage.equilibrate(composition, assemblage, constraints, guesses) $ bdg = burnman.minerals.SLB_2011.mg_fe_al_bridgmanite() $ bdg.set_composition([0.9, 0.1, 0.0]) $ bdg.set_state(30.e9, 2000.) $ print bdg.gibbs, bdg.S, bdg.K_T, bdg.heat_capacity_p Example applications • High pressure experimental data analysis • Geodynamic model interrogation for seismic properties • Forward-modelling for seismic proﬁles using bulk/mineral composition and temperature constraints • Self-consistent planetary mass/moment of inertia/pressure-depth calculations An example of using BurnMan in an MCMC inversion of outer core elastic properties will be presented at DI43A-2660 3D reference models So far, BurnMan has provided a range of seismically derived 1D radial models of the elastic properties of the Earth. An average 1D model however is not a true representation of the Earth. Constraints from fully 3D seismic models are needed to test, which thermo-chemical models ﬁt the velocity variations observed. Currently, we are working on provi