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

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.

Energy-effective Grinding of Inorganic Solids Using Organic Additives.

We present our research findings related to new formulations of the organic additives (grinding aids) needed for the efficient grinding of inorganic solids. Even though the size reduction phenomena of the inorganic solid particles in a ball mill is purely a physical process, the addition of grinding aids in milling media introduces a complex physicochemical process. In addition to further gain in productivity, the organic additive helps to reduce the energy needed for grinding, which in the case of cement clinker has major environmental implications worldwide. This is primarily due to the tremendous amounts of cement produced and almost 30% of the associated electrical energy is consumed for grinding. In this paper, we examine the question of how to optimize these grinding aids linking molecular insight into their working mechanisms, and also how to design chemical additives of improved performance for industrial comminution.

A force field for tricalcium aluminate to characterize surface properties, initial hydration, and organically modified interfaces in atomic resolution.

Tricalcium aluminate (C3A) is a major phase of Portland cement clinker and some dental root filling cements. An accurate all-atom force field is introduced to examine structural, surface, and hydration properties as well as organic interfaces to overcome challenges using current laboratory instrumentation. Molecular dynamics simulation demonstrates excellent agreement of computed structural, thermal, mechanical, and surface properties with available experimental data. The parameters are integrated into multiple potential energy expressions, including the PCFF, CVFF, CHARMM, AMBER, OPLS, and INTERFACE force fields. This choice enables the simulation of a wide range of inorganic-organic interfaces at the 1 to 100 nm scale at a million times lower computational cost than DFT methods. Molecular models of dry and partially hydrated surfaces are introduced to examine cleavage, agglomeration, and the role of adsorbed organic molecules. Cleavage of crystalline tricalcium aluminate requires approximately 1300 mJ m(-2) and superficial hydration introduces an amorphous calcium hydroxide surface layer that reduces the agglomeration energy from approximately 850 mJ m(-2) to 500 mJ m(-2), as well as to lower values upon surface displacement. The adsorption of several alcohols and amines was examined to understand their role as grinding aids and as hydration modifiers in cement. The molecules mitigate local electric fields through complexation of calcium ions, hydrogen bonds, and introduction of hydrophobicity upon binding. Molecularly thin layers of about 0.5 nm thickness reduce agglomeration energies to between 100 and 30 mJ m(-2). Molecule-specific trends were found to be similar for tricalcium aluminate and tricalcium silicate. The models allow quantitative predictions and are a starting point to provide fundamental understanding of the role of C3A and organic additives in cement. Extensions to impure phases and advanced hydration stages are feasible.

Freezing of fluids confined between mica surfaces.

Using grand ensemble simulations, we show that octamethyl-cyclo-tetra-siloxane (OMCTS) confined between two mica surfaces can form a variety of frozen phases which undergo solid-solid transitions as a function of the separation between the surfaces. For atomically smooth mica surfaces, the following sequence of transitions 1[triangle up] --> 1[triangle up]b --> 2B --> 2 square --> 2[triangle up] are observed in the one- and two-layered regimes, where n[triangle up], n[square], and nB denote triangular, square, and buckled phases, respectively, with the prefix n denoting the number of confined layers. The presence of potassium on mica is seen to have a strong influence on the degree of order induced in the fluid. The sequence of solid-solid transitions that occurs with the smooth mica surface is no longer observed. When equilibrated with a state point near the liquid-solid transition, a counterintuitive freezing scenario is observed in the presence of potassium. Potassium disrupts in-plane ordering in the fluid in contact with the mica surface, and freezing is observed only in the inner confined layers. The largest mica separations at which frozen phases were observed ranged from separations that could accommodate six to seven fluid layers. The extent of freezing and the square-to-triangular lattice transition was found to be sensitive to the presence of potassium as well as the thermodynamic conditions of the bulk fluid. The implications of our results on interpretation of surface force experiments as well as the generic phase behavior of confined soft spheres is discussed.