Abstract:The simulation of metals, oxides, and hydroxides can accelerate the design of therapeutics, alloys, catalysts, cement-based materials, ceramics, bioinspired composites, and glasses. Here we introduce the INTERFACE force field (IFF) and surface models for α-Al2O3, α-Cr2O3, α-Fe2O3, NiO, CaO, MgO, β-Ca(OH)2, β-Mg(OH)2, and β-Ni(OH)2. The force field parameters are nonbonded, including atomic charges for Coulomb interactions, Lennard-Jones (LJ) potentials for van der Waals interactions with 12-6 and 9-6 options, and harmonic bond stretching for hydroxide ions. The models outperform DFT calculations and earlier atomistic models (Pedone, ReaxFF, UFF, CLAYFF) up to 2 orders of magnitude in reliability, compatibility, and interpretability due to a quantitative representation of chemical bonding consistent with other compounds across the periodic table and curated experimental data for validation. The IFF models exhibit average deviations of 0.2% in lattice parameters, <10% in surface energies (to the extent known), and 6% in bulk moduli relative to experiments. The parameters and models can be used with existing parameters for solvents, inorganic compounds, organic compounds, biomolecules, and polymers in IFF, CHARMM, CVFF, AMBER, OPLS-AA, PCFF, and COMPASS, to simulate bulk oxides, hydroxides, electrolyte interfaces, and multiphase, biological, and organic hybrid materials at length scales from atoms to micrometers. The nonbonded character of the models also enables the analysis of mixed oxides, glasses, and certain chemical reactions, and well-performing nonbonded models for silica phases, SiO2, are introduced. Automated model building is available in the CHARMM-GUI Nanomaterial Modeler. We illustrate applications of the models to predict the structure of mixed oxides, and energy barriers of ion migration, as well as binding energies of water and organic molecules in outstanding agreement with experimental data and calculations at the CCSD(T) level. Examples of model building for hydrated, pH-sensitive oxide surfaces to simulate solid-electrolyte interfaces are discussed.

A polarizable interatomic force field for TiO$_2$ parameterized using density functional theory

Theoretical Investigation of the High Pressure Structure, Lattice Dynamics, Phase Transition, and Thermal Equation of State of Titanium Metal

Pressure‐induced phase transition and elastic properties of TiO2 polymorphs

Pressure-induced phase transition and elastic properties of TiO<sub>2</sub>polymorphs

A systematic DFT study of structure and electronic properties of titanium dioxide

Rutile, anatase, brookite and titania thin film from Hubbard corrected and hybrid DFT

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Investigating Polaron Formation in Anatase and Brookite TiO2 by Density Functional Theory with Hybrid-Functional and DFT + U Methods

DFT study with inclusion of the Grimme potential on anatase TiO2: structure, electronic, and vibrational analyses

Using Density Functional Theory to Model Realistic TiO2 Nanoparticles, Their Photoactivation and Interaction with Water

Phase Stability of TiO$_2$ Polymorphs from Diffusion Quantum Monte Carlo

Structural Stability of TiO2 at High Pressure in Density-Functional Theory Based Calculations.

Ab Initio, Physically Motivated Force Fields for Co2 Adsorption in Zeolitic Imidazolate Frameworks

Ab Initio Study of Phase Stability in Doped TiO2

Investigation on the stability, mechanism and electronic properties of rutile TiO2 via first-principle calculations

Ab initio multi-scale modeling of ferroelectrics: The case of PbTiO3

Phase behaviour of empirical potentials of titanium dioxide

Phase transitions in titanium with an analytic bond-order potential

Accurate Force Fields for Atomistic Simulations of Oxides, Hydroxides, and Organic Hybrid Materials up to the Micrometer Scale

Atomic-Scale Modeling of the Dynamics of Titanium Oxidation

Electronic and Optical Properties of TiO2 Thin Films: Combined Experimental and Theoretical Study