Comparison of model potentials for molecular-dynamics simulations of silica

Structure, Properties, and Applications of Silica Nanoparticles: Recent Theoretical Modeling Advances, Challenges, and Future Directions

Silica nanoparticles (SNPs), one of the most widely researched materials in modern science, are now commonly exploited in surface coatings, biomedicine, catalysis, and engineering of novel self-assembling materials. Theoretical approaches are invaluable to enhancing fundamental understanding of SNP properties and behavior. Tremendous research attention is dedicated to modeling silica structure, the silica-water interface, and functionalization of silica surfaces for tailored applications. In this review, the range of theoretical methodologies are discussed that have been employed to model bare silica and functionalized silica. The evolution of silica modeling approaches is detailed, including classical, quantum mechanical, and hybrid methods and highlight in particular the last decade of theoretical simulation advances. It is started with discussing investigations of bare silica systems, focusing on the fundamental interactions at the silica-water interface, following with a comprehensively review of the modeling studies that examine the interaction of silica with functional ligands, peptides, ions, surfactants, polymers, and carbonaceous species. The review is concluded with the perspective on existing challenges in the field and promising future directions that will further enhance the utility and importance of the theoretical approaches in guiding the rational design of SNPs for applications in engineering and biomedicine.

Gaussian Process Regression for Materials and Molecules

We provide an introduction to Gaussian process regression (GPR) machine-learning methods in computational materials science and chemistry. The focus of the present review is on the regression of atomistic properties: in particular, on the construction of interatomic potentials, or force fields, in the Gaussian Approximation Potential (GAP) framework; beyond this, we also discuss the fitting of arbitrary scalar, vectorial, and tensorial quantities. Methodological aspects of reference data generation, representation, and regression, as well as the question of how a data-driven model may be validated, are reviewed and critically discussed. A survey of applications to a variety of research questions in chemistry and materials science illustrates the rapid growth in the field. A vision is outlined for the development of the methodology in the years to come.

Modelling the atomic arrangement of amorphous 2D silica: a network analysis

The recent experimental discovery of a semi two-dimensional silica glass has offered a realistic description of the random network theory of a silica glass structure, initially discussed by Zachariasen.

Multiple pathways in pressure-induced phase transition of coesite

High-pressure single-crystal X-ray diffraction method with precise control of hydrostatic conditions, typically with helium or neon as the pressure-transmitting medium, has significantly changed our view on what happens with low-density silica phases under pressure. Coesite is a prototype material for pressure-induced amorphization. However, it was found to transform into a high-pressure octahedral (HPO) phase, or coesite-II and coesite-III. Given that the pressure is believed to be hydrostatic in two recent experiments, the different transformation pathways are striking. Based on molecular dynamic simulations with an ab initio parameterized potential, we reproduced all of the above experiments in three transformation pathways, including the one leading to an HPO phase. This octahedral phase has an oxygen hcp sublattice featuring 2 × 2 zigzag octahedral edge-sharing chains, however with some broken points (i.e., point defects). It transforms into α-PbO₂ phase when it is relaxed under further compression. We show that the HPO phase forms through a continuous rearrangement of the oxygen sublattice toward hcp arrangement. The high-pressure amorphous phases can be described by an fcc and hcp sublattice mixture.

Structure and dynamics in network-forming materials

The study of the structure and dynamics of network-forming materials is reviewed. Experimental techniques used to extract key structural information are briefly considered. Strategies for building simulation models, based on both targeting key (experimentally-accessible) materials and on systematically controlling key model parameters, are discussed. As an example of the first class of materials, a key target system, SiO₂, is used to highlight how the changing structure with applied pressure can be effectively modelled (in three dimensions) and used to link to both experimental results and simple structural models. As an example of the second class the topology of networks of tetrahedra in the MX₂ stoichiometry are controlled using a single model parameter linked to the M-X-M bond angles. The evolution of ordering on multiple length-scales is observed as are the links between the static structure and key dynamical properties. The isomorphous relationship between the structures of amorphous Si and SiO₂ is discussed as are the similarities and differences in the phase diagrams, the latter linked to potential polyamorphic and 'anomalous' (e.g. density maxima) behaviour. Links to both two-dimensional structures for C, Si and Ge and near-two-dimensional bilayers of SiO₂ are discussed. Emerging low-dimensional structures in low temperature molten carbonates are also uncovered.

A molecular dynamics study of ambient and high pressure phases of silica: structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, α-PbO2-type, and pyrite-type - for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase.

Influence of polarizability on metal oxide properties studied by molecular dynamics simulations

We have studied the dependence of metal oxide properties in molecular dynamics (MD) simulations on the polarizability of oxygen ions. We present studies of both liquid and crystalline structures of silica (SiO(2)), magnesia (MgO) and alumina (Al(2)O(3)). For each of the three oxides, two separately optimized sets of force fields were used: (i) long-range Coulomb interactions between oxide and metal ions combined with a short-range pair potential; (ii) extension of force field (i) by adding polarizability to the oxygen ions. We show that while an effective potential of type (i) without polarizable oxygen ions can describe radial distributions and lattice constants reasonably well, potentials of type (ii) are required to obtain correct values for bond angles and the equation of state. The importance of polarizability for metal oxide properties decreases with increasing temperature.

Tuning oxygen packing in silica by nonhydrostatic pressure

The transformation of SiO2 from low pressure tetrahedral phases into denser octahedral phases takes place via the collapse of the oxygen sublattice into a close-packed arrangement. The transition paths and the resulting products are known to be affected by the presence of anisotropic stresses, which are difficult to control, so interpretation of the experimental results is problematic. Based on nonhydrostatic molecular dynamics simulations, we show that the collapse of the oxygen sublattice in the specific case of cristobalite is concomitant with the disappearance of tetrahedral units and that non hydrostatic stresses can be tuned to yield phases with different oxygen close-packed sublattices, including the alpha-PbO2-like phase, for which we provide a microscopic formation path, and phases with a cubic close packing, like anatase, not seen in experiments yet.

Amorphous silica modeled with truncated and screened Coulomb interactions: A molecular dynamics simulation study

We show that finite-range alternatives to the standard long-range pair potential for silica by van Beest et al. [Phys. Rev. Lett. 64, 1955 (1990)] might be used in molecular dynamics simulations. We study two such models that can be efficiently simulated since no Ewald summation is required. We first consider the Wolf method, where the Coulomb interactions are truncated at a cutoff distance rc such that the requirement of charge neutrality holds. Various static and dynamic quantities are computed and compared to results from simulations using Ewald summations. We find very good agreement for rc approximately 10 A. For lower values of rc, the long-range structure is affected which is accompanied by a slight acceleration of dynamic properties. In a second approach, the Coulomb interaction is replaced by an effective Yukawa interaction with two new parameters determined by a force fitting procedure. The same trend as for the Wolf method is seen. However, slightly larger cutoffs have to be used in order to obtain the same accuracy with respect to static and dynamic quantities as for the Wolf method.

Development of an empirical force field for silica. Application to the quartz-water interface

Interactions of pulverized crystalline silica with biological systems, including the lungs, cause cell damage, inflammation, and apoptosis. To allow computational atomistic modeling of these pathogenic processes, including interactions between silica surfaces and biological molecules, new parameters for quartz, compatible with the CHARMM empirical force field were developed. Parameters were optimized to reproduce the experimental geometry of alpha-quartz, ab initio vibrational spectra, and interactions between model compounds and water. The newly developed force field was used to study interactions of water with two singular surfaces of alpha-quartz, (011) and (100). Properties monitored and analyzed include the variation of the density of water molecules in the plane perpendicular to the surface, disruption of the water H-bond network upon adsorption, and space-time correlations of water oxygen atoms in terms of Van Hove self-correlation functions. The vibrational density of states spectra of water in confined compartments were also computed and compared with experimental neutron-scattering results. Both the attenuation and shifting to higher frequencies of the hindered translational peaks upon confinement are clearly reproduced by the model. However, an upshift of librational peaks under the conditions of model confinement still remains underrepresented at the current empirical level.

Infrared and Raman spectra of silica polymorphs from anab initioparametrized polarizable force field

The general aim of this study is to test the reliability of polarizable model potentials for the prediction of vibrational (infrared and Raman) spectra in highly anharmonic systems such as high temperature crystalline phases. By using an ab initio parametrized interatomic potential for SiO2 and molecular dynamics simulations, we calculate the infrared and Raman spectra for quartz, cristobalite, and stishovite at various thermodynamic conditions. The model is found to perform very well in the prediction of infrared spectra. Raman peak positions are also reproduced very well by the model; however, Raman intensities calculated by explicitly taking the derivative of the polarizability with respect to the atomic displacements are found to be in poorer agreement than intensities calculated using a parametrized "bond polarizability" model. Calculated spectra for the high temperature beta phases, where the role of dynamical disorder and anharmonicities is predominant, are found to be in excellent agreement with experiments. For the octahedral phases, our simulations are able to reproduce changes in the Raman spectra across the rutile-to-CaCl2 transition around 50 GPa, including the observed phonon softening.