Molecular dynamics simulation of amorphous SiO2 nanoparticles

Tuning the Properties of Halide Nanocomposite Solid Electrolytes with Diverse Oxides for All-Solid-State Batteries

Herein, we report halide nanocomposite solid electrolytes (HNSEs) that integrate diverse oxides with alterations that allow tuning of their ionic conductivity, (electro)chemical stability, and specific density. A two-step mechanochemical process enabled the synthesis of multimetal (or nonmetal) HNSEs, MO2-2Li2ZrCl6, as verified by pair distribution function and X-ray diffraction analyses. The multimetal (or nonmetal) HNSE strategy increases the ionic conductivity of Li2ZrCl6 from 0.40 to 0.82 mS cm-1. Additionally, cyclic voltammetry test findings corroborated the enhanced passivating properties of the HNSEs. Notably, incorporating SiO2 into HNSEs leads to a substantial reduction in the specific density of HNSEs, demonstrating their strong potential for achieving a high energy density and lowering costs. Fluorinated SiO2-2Li2ZrCl5F HNSEs exhibited enhanced interfacial compatibility with Li6PS5Cl and LiCoO2 electrodes. Cells employing SiO2-2Li2ZrCl5F with LiCoO2 exhibit superior electrochemical performance delivering the initial discharge capacity of 162 mA h g-1 with 93.7% capacity retention at the 100th cycle at 60 °C.

Modeling Colloidal Particle Aggregation Using Cluster Aggregation with Multiple Particle Interactions

In this study, we investigate the aggregation dynamics of colloidal silica by generating simulated structures and comparing them to experimental data gathered through scanning transmission electron microscopy (STEM). More specifically, diffusion-limited cluster aggregation and reaction-limited cluster aggregation models with different functions for the probability of particles sticking upon contact were used. Aside from using a constant sticking probability, the sticking probability was allowed to depend on the masses of the colliding clusters and on the number of particles close to the collision between clusters. The different models of the sticking probability were evaluated based on the goodness-of-fit of spatial summary statistics. Furthermore, the models were compared to the experimental data by calculating the structures' fractal dimension and mass transport properties from simulations of flow and diffusion. The sticking probability, depending on the interaction with multiple particles close to the collision site, led to structures most similar to the STEM data.

Developments in Atomistic and Nano Structure Evolution Mechanisms of Molten Slag Using Atomistic Simulation Methods

Molten slag has different properties depending on its composition. The relationship between its composition, structure, and properties has been the focus of attention in industrial manufacturing processes. This review describes the atomistic scale mechanisms by which oxides of different compositions affect the properties and structure of slag, and depicts the current state of research in the atomic simulation of molten slag. At present, the research on the macroscopic properties of molten slag mainly focuses on viscosity, free-running temperature, melting point, and desulphurization capacity. Regulating the composition has become the most direct and effective way to control slag properties. Analysis of the microevolution mechanism is the fundamental way to grasp the macroscopic properties. The microstructural evolution mechanism, especially at the atomic and nanoscale of molten slag, is reviewed from three aspects: basic oxides, acidic oxides, and amphoteric oxides. The evolution of macroscopic properties is analyzed in depth through the evolution of the atomic structure. Resolution of the macroscopic properties of molten slag by the atomic structure plays a crucial role in the development of fundamental theories of physicochemistry.

First-Principles Prediction of Amorphous Silica Nanoparticle Surface Charge: Effect of Size, pH, and Ionic Strength

The quantification of surface charge properties of silica nanoparticles is essential for several applications. To determine these properties, many experimental and theoretical methods have been introduced, which are time-consuming and/or challenging to use. In this study, a first-principles approach is developed to determine the surface charge properties of amorphous silica nanoparticles against the nanoparticle size, pH, and ionic strength without relying on experimental data. An amorphous silica nanoparticle of 1.34 nm diameter is simulated by using integrated molecular dynamics and Monte Carlo methods. A detailed analysis of the nanoparticle structure is provided by analyzing the types of silanol groups on the surface. Moreover, a model is developed to estimate the probability distribution of the surface silanol groups based on the nearest neighbor distances and the diameter of the nanoparticle to determine the number of surface silanols on larger nanoparticles. Thereafter, a computational chemistry approach is used to calculate the acid dissociation constants of the corresponding deprotonation reactions. The calculated constants and the point of zero charge value are in excellent agreement with experiments. The surface charge properties of the nanoparticle with various diameters are then estimated by using a mean-field model at different pH and ionic strength values. The results of the developed model are compared to the Poisson-Boltzmann equation as a reference model. The developed model predictions agree well with the reference model for low and mid-electrolyte concentrations (1 and 10 mM) and small nanoparticles (smaller than 100 nm). However, the developed model seems to qualitatively predict the surface charge properties more accurately than the Poisson-Boltzmann model for high electrolyte concentrations.

Atomistic Simulation of Physical Vapor Deposition of Optical Thin Films

A review of the methods and results of atomistic modeling of the deposition of thin optical films and a calculation of their characteristics is presented. The simulation of various processes in a vacuum chamber, including target sputtering and the formation of film layers, is considered. Methods for calculating the structural, mechanical, optical, and electronic properties of thin optical films and film-forming materials are discussed. The application of these methods to studying the dependences of the characteristics of thin optical films on the main deposition parameters is considered. The simulation results are compared with experimental data.

Antioxidant Effects of Bioactive Glasses (BGs) and Their Significance in Tissue Engineering Strategies

Elevated levels of oxidative stress are usually observed following injuries, leading to impaired tissue repair due to oxidation-related chronic inflammation. Several attempts have been made to manage this unfavorable situation, and the use of biomaterials with antioxidant activity is showing great promise in tissue engineering and regenerative medicine approaches. Bioactive glasses (BGs) are a versatile group of inorganic substances that exhibit an outstanding regenerative capacity for both hard and soft damaged tissues. The chemical composition of BGs provides a great opportunity for imparting specific biological activities to them. On this point, BGs may easily become antioxidant substances through simple physicochemical modifications. For example, particular antioxidant elements (mostly cerium (Ce)) can be added to the basic composition of the glasses. On the other hand, grafting natural antioxidant substances (e.g., polyphenols) on the BG surface is feasible for making antioxidant substitutes with promising results in vitro. Mesoporous BGs (MBGs) were demonstrated to have unique merits compared with melt-derived BGs since they make it possible to load antioxidants and deliver them to the desired locations. However, there are actually limited in vivo experimental studies on the capability of modified BGs for scavenging free radicals (e.g., reactive oxygen species (ROS)). Therefore, more research is required to determine the actual potential of BGs in decreasing oxidative stress and subsequently improving tissue repair and regeneration. The present work aims to highlight the potential of different types of BGs in modulating oxidative stress and subsequently improving tissue healing.

Gas Adsorption and Diffusion Behaviors in Interfacial Systems Composed of a Polymer of Intrinsic Microporosity and Amorphous Silica: A Molecular Simulation Study

We investigate the adsorption and diffusion behaviors of CO₂, CH₄, and N₂ in interfacial systems composed of a polymer of intrinsic microporosity (PIM-1) and amorphous silica using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations. We build model systems of mixed matrix membranes (MMMs) with PIM-1 chains sandwiched between silica surfaces. Gas adsorption analysis using GCMC simulations shows that gas molecules are preferentially adsorbed in microcavities distributed near silica surfaces, resulting in an increase in the solubility coefficients of CO₂, CH₄, and N₂ compared to bulk PIM-1. In contrast, diffusion coefficients obtained from MD simulations and then calibrated using the dual-mode sorption model show different tendencies depending on gas species: CO₂ diffusivity decreases in MMMs compared to PIM-1, whereas CH₄ and N₂ diffusivities increase. These differences are attributed to competing effects of silica surfaces: the emergence of larger pores as a result of chain packing disruption, which enhances gas diffusion, and a quadrupole-dipole interaction between gas molecules and silica surface hydroxyl groups, which retards gas diffusion. The former has a greater impact on CH₄ and N₂ diffusivities, whereas the latter has a greater impact on CO₂ diffusivity due to the strong quadrupole-dipole interaction between CO₂ and surface hydroxyls. These findings add to our understanding of gas adsorption and diffusion behaviors in the vicinity of PIM-1/silica interfaces, which are unobtainable in experimental studies.

SiO2-Ag Composite as a Highly Virucidal Material: A Roadmap that Rapidly Eliminates SARS-CoV-2

COVID-19, as the cause of a global pandemic, has resulted in lockdowns all over the world since early 2020. Both theoretical and experimental efforts are being made to find an effective treatment to suppress the virus, constituting the forefront of current global safety concerns and a significant burden on global economies. The development of innovative materials able to prevent the transmission, spread, and entry of COVID-19 pathogens into the human body is currently in the spotlight. The synthesis of these materials is, therefore, gaining momentum, as methods providing nontoxic and environmentally friendly procedures are in high demand. Here, a highly virucidal material constructed from SiO2-Ag composite immobilized in a polymeric matrix (ethyl vinyl acetate) is presented. The experimental results indicated that the as-fabricated samples exhibited high antibacterial activity towards Escherichia coli (E. coli) and Staphylococcus aureus (S. aureus) as well as towards SARS-CoV-2. Based on the present results and radical scavenger experiments, we propose a possible mechanism to explain the enhancement of the biocidal activity. In the presence of O2 and H2O, the plasmon-assisted surface mechanism is the major reaction channel generating reactive oxygen species (ROS). We believe that the present strategy based on the plasmonic effect would be a significant contribution to the design and preparation of efficient biocidal materials. This fundamental research is a precedent for the design and application of adequate technology to the next-generation of antiviral surfaces to combat SARS-CoV-2.

Computational Study of APTES Surface Functionalization of Diatom-like Amorphous SiO2 Surfaces for Heavy Metal Adsorption

The amorphous silica (SiO₂) shell on diatom frustules is a highly attractive biomaterial for removing pollutants from aquatic ecosystems. The surface activity of silica can be enhanced by modification with organosilanes. In this work, we present an atomic-level theoretical study based on molecular dynamics and dispersion-corrected density functional theory calculations on the surface stability and adsorption of heavy metal (HM) compounds on silane- and 3-aminopropyltriethoxysilane (APTES)-covered SiO₂ surfaces. Our simulations show that at low APTES coverage, the molecular adsorption of Cd(OH)₂ and HgCl₂ is more favorable near the modifier, compared to As(OH)₃ that binds at the hydroxylated region on silica. At higher coverages, the metallic compounds are preferentially adsorbed by the terminating amino group on the surface, whereas the adsorption in the region between APTES and the oxide surface is also spontaneous. The adsorption is strongly driven by van der Waals interactions at the highly covered surface, where the consideration of dispersion corrections reduces the modifier-adsorbate interatomic distances and increases the adsorption energy by ca. 0.4-0.7 eV. The adsorption of water is favorable, although it is generally weaker than for the HM compounds. Based on our results, we conclude that the addition of APTES modifiers on silica increases the adsorption strength and provides extra binding sites for the adsorption of HM pollutants. These outcomes can be used for the design of more efficient structures of biomaterials for depollution of HMs.

Earthquake lubrication and healing explained by amorphous nanosilica

During earthquake propagation, geologic faults lose their strength, then strengthen as slip slows and stops. Many slip-weakening mechanisms are active in the upper-mid crust, but healing is not always well-explained. Here we show that the distinct structure and rate-dependent properties of amorphous nanopowder (not silica gel) formed by grinding of quartz can cause extreme strength loss at high slip rates. We propose a weakening and related strengthening mechanism that may act throughout the quartz-bearing continental crust. The action of two slip rate-dependent mechanisms offers a plausible explanation for the observed weakening: thermally-enhanced plasticity, and particulate flow aided by hydrodynamic lubrication. Rapid cooling of the particles causes rapid strengthening, and inter-particle bonds form at longer timescales. The timescales of these two processes correspond to the timescales of post-seismic healing observed in earthquakes. In natural faults, this nanopowder crystallizes to quartz over 10s–100s years, leaving veins which may be indistinguishable from common quartz veins.

Tectonic faults weaken during slip in order to accelerate and produce earthquakes. Here the authors show a mechanism for weakening faults through the transformation of quartz to amorphous nanoparticulate wear powders that lubricate friction experiments, and transform back to quartz under geologic conditions.

On the Structure of Amorphous Mesoporous Silica Nanoparticles by Aberration-Corrected STEM

Mesoporous silica materials have demonstrated a vast spectrum of applications, stimulating an intensive field of study due to their potential use as nanocarriers. Nonetheless, when produced at the nanoscale, their structural characterization is hindered due to the re-arrangement of the pores. To address this issue, this work combines molecular dynamics simulations with electron microscopy computer simulations and experimental results to provide an insight into the structure of amorphous mesoporous silica nanoparticles. The amorphous silica model is prepared using a simple melt-quench molecular dynamics method, while the reconstruction of the mesoporous nanoparticles is carried out using a methodology to avoid false symmetry in the final model. Simulated scanning transmission electron microscopy images are compared with experimental images, revealing the existence of structural domains, created by the misalignment of the pores to compensate the surface tension of these spherical nanoparticles.

Propane-Water Mixtures Confined within Cylindrical Silica Nanopores: Structural and Dynamical Properties Probed by Molecular Dynamics

Despite the multiple length and time scales over which fluid-mineral interactions occur, interfacial phenomena control the exchange of matter and impact the nature of multiphase flow, as well as the reactivity of C-O-H fluids in geologic systems. In general, the properties of confined fluids, and their influence on porous geologic phenomena are much less well understood compared to those of bulk fluids. We used equilibrium molecular dynamics simulations to study fluid systems composed of propane and water, at different compositions, confined within cylindrical pores of diameter ∼16 Å carved out of amorphous silica. The simulations are conducted within a single cylindrical pore. In the simulated system all the dangling silicon and oxygen atoms were saturated with hydroxyl groups and hydrogen atoms, respectively, yielding a total surface density of 3.8 -OH/nm2. Simulations were performed at 300 K, at different bulk propane pressures, and varying the composition of the system. The structure of the confined fluids was quantified in terms of the molecular distribution of the various molecules within the pore as well as their orientation. This allowed us to quantify the hydrogen bond network and to observe the segregation of propane near the pore center. Transport properties were quantified in terms of the mean square displacement in the direction parallel to the pore axis, which allows us to extract self-diffusion coefficients. The diffusivity of propane in the cylindrical pore was found to depend on pressure, as well as on the amount of water present. It was found that the propane self-diffusion coefficient decreases with increasing water loading because of the formation of water bridges across the silica pores, at sufficiently high water content, which hinder propane transport. The rotational diffusion, the lifespan of hydrogen bonds, and the residence time of water molecules at contact with the silica substrate were quantified from the simulated trajectories using the appropriate autocorrelation functions. The simulations contribute to a better understanding of the molecular phenomena relevant to the behavior of fluids in the subsurface.

SUePDF: a program to obtain quantitative pair distribution functions from electron diffraction data

SUePDF is a graphical user interface program written in MATLAB to achieve quantitative pair distribution functions (PDFs) from electron diffraction data. The program facilitates structural studies of amorphous materials and small nanoparticles using electron diffraction data from transmission electron microscopes. It is based on the physics of electron scattering as well as the total scattering methodology. A method of background modeling is introduced to treat the intensity tail of the direct beam, inelastic scattering and incoherent multiple scattering. Kinematical electron scattering intensity is scaled using the electron scattering factors. The PDFs obtained after Fourier transforms are normalized with respect to number density, nanoparticle form factor and the non-negativity of probability density. SUePDF is distributed as free software for academic users.

Study of thermal conductivity of ice clusters after impact deposition on the silica surfaces using the ReaxFF reactive force field

During aircraft or spacecraft missions, ice accumulates on different parts of their surface elements. An important parameter affecting the ability to remove this ice from the surface is the heat transfer characteristics of the accumulated ice. The ice heat transfer is related to the process of ice formation and its density and internal structure. In this study we investigate the effects of the ice and silica structure and the ice cluster attachment mechanism to the silica surface on the thermal conductivity (TC) of the attached ice cluster using the ReaxFF reactive force field. The purpose of this study is to investigate the thermal transport in amorphous and crystalline ice after high-velocity deposition on the silica surfaces. A dual thermostat method has been applied for the calculation of TC values. The validity of this method has been verified by comparing the calculated values of TC for crystal and amorphous ice with available experimental values. Our calculations show that the TC values of both crystal and amorphous ice drop after deposition on the silica surfaces. This decrease in the TC is more significant for the ice deposition on suboxide silica surfaces. Furthermore, crystal ice shows higher TC values than amorphous ice after accumulation. However, when crystal ice impacts on the silica surface at 1 km s(-1) impact speed, the crystalline shape of the ice cluster is lost to a considerable level and the TC values obtained for the ice clusters in such cases are closer to amorphous ice TC values. We observed a decrease in the TC values when ionic species are added inside the ice clusters.

The antioxidant properties of Ce-containing bioactive glass nanoparticles explained by Molecular Dynamics simulations

Molecular dynamics simulations of two glass nanoparticles with composition 25Na2O·25CaO 50SiO2 mol% (Ce-K NP) and 46.1SiO2·24.4Na2O·26.9CaO· 2.6P2O5 mol.% (Ce-BG NP) doped with 3.6 mol% of CeO2 have been carried out in order to explain the enhanced antioxidant properties of the former glass with respect to the latter.

The present models show that the different catalase mimetic activity of the two NPs is related to the Ce3+/Ce4+ ratio exposed at their surface. In fact, this ratio is about 3.5 and 13 in the bulk and at the surface of the Ce-BG NP, and 1.0 and 2.1 in the bulk and at the surface of the Ce-K NPs, respectively. Since both oxidation states are necessary for the catalysis of the dismutation reaction of hydrogen peroxides, NPs with a very high Ce3+/Ce4+ ratio possess poorer antioxidant properties.

Moreover, our simulations reveal that the already low silicate connectivity found in the bulk glasses examined here is further reduced on the nanoparticle surface, whereas the Na+/Ca2+ ratio rapidly increases. Sodium, calcium and cerium sites in proximity of the surface are found to be under-coordinated, prone to quickly react with water present in physiological environments, thus accelerating the glass biodegradation

Stress in titania nanoparticles: an atomistic study

The size-dependent surface and bulk stresses intrinsic to titania nanoparticles are investigated using atomistic simulation. Surface charge is also shown to induce a significant tensile stress.

Influence of capillary bridge formation onto the silica nanoparticle interaction studied by grand canonical Monte Carlo simulations

Adhesion forces between nanoparticles strongly depend on the amount of adsorbed condensed water from ambient atmosphere. Liquid water forms bridges in the cavities separating the particles, giving rise to the so-called capillary forces which in most cases dominate the van der Waals and long-range electrostatic interactions. Capillary forces promote the undesirable agglomeration of particles to large clusters, thereby hindering the flowability of dry powders in process containers. In process engineering macroscopic theories based on the Laplace pressures are used to estimate the strength of the capillary forces. However, especially for low relative humidity and when the wetting of rough or small nanoparticles is studied, those theories can fail. Molecular dynamic simulations can help to give better insight into the water-particle interface. The simulated force versus distance curve as well as adhesion forces and the adsorption isotherm for silica nanoparticles at varying relative humidity will be discussed in comparison to experiments, theories, and simulations.

Calculation of noncontact forces between silica nanospheres

Quantification of the interactions between nanoparticles is important in understanding their dynamic behaviors and many related phenomena. In this study, molecular dynamics simulation is used to calculate the interaction potentials (i.e., van der Waals attraction, Born repulsion, and electrostatic interaction) between two silica nanospheres of equal radius in the range of 0.975 to 5.137 nm. The results are compared with those obtained from the conventional Hamaker approach, leading to the development of modified formulas to calculate the van der Waals attraction and Born repulsion between nanospheres, respectively. Moreover, Coulomb's law is found to be valid for calculating the electrostatic potential between nanospheres. The developed formulas should be useful in the study of the dynamic behaviors of nanoparticle systems under different conditions.

Simulation of Forces between Humid Amorphous Silica Surfaces: A Comparison of Empirical Atomistic Force Fields

Atmospheric humidity strongly influences the interactions between dry granular particles in process containers. To reduce the energy loss in industrial production processes caused by particle agglomeration, a basic understanding of the dependence of particle interactions on humidity is necessary. Hence, in this study, molecular dynamic simulations were carried out to calculate the adhesion between silica surfaces in the presence of adsorbed water. For a realistic description, the choice of force field is crucial. Because of their frequent use and transferability to biochemical systems, the Clay and CWCA force fields were investigated with respect to their ability to describe the water-silica interface in comparison to the more advanced Reax force field, ab initio calculations, and experiments.

Molecular dynamics simulation of liquid and amorphous Fe nanoparticles

Liquid and amorphous Fe nanoparticles have been studied in spherical models with sizes of 3, 4 and 5 nm. We use the molecular dynamics (MD) method with the Pak-Doyam interatomic pair potential. Structural properties of Fe nanoparticles have been studied via radial distribution function (RDF), interatomic distances, coordination number distribution and radial density profile. Local icosahedral order was detected via using Honeycutt-Andersen (HA) analysis compared with that observed in the bulk counterparts. Surface and core structures of Fe nanoparticles have been studied in detail. Potential and surface energies of liquid and amorphous Fe nanoparticles have been calculated and discussed in addition to the size dependence of a glass transition temperature.

Static and thermodynamic properties of liquid and amorphous Fe(2)O(3) nanoparticles

Molecular dynamics (MD) simulation of the structure and thermodynamic properties of liquid and amorphous Fe(2)O(3) nanoparticles, with different sizes ranging from 2 to 5 nm, has been carried out based on the Born-Mayer pair potentials under non-periodic boundary conditions. The partial radial distribution functions (PRDFs), coordination number distributions, bond-angle distributions and interatomic distances of liquid nanoparticles at 3500 K were calculated in detail. In addition, we show the radial density profile and stoichiometry in nanoparticles. The temperature dependence of the surface structure and surface energy of nanoparticles has been obtained and presented. Moreover, the size dependence of the glass transition temperature has been found and discussed.

Surface structure and structural point defects of liquid and amorphous aluminosilicate nanoparticles

The surface structure of liquid and amorphous aluminosilicate nanoparticles of composition Al(2)O(3)·2SiO(2) has been investigated in a model of different sizes ranging from 2.0 to 5.0 nm with the Born-Mayer type pair potential under non-periodic boundary conditions. Models have been obtained by cooling from the melts at a constant density of 2.6 g cm(-3) via molecular dynamics (MD) simulation. The surface structure has been investigated via the coordination number, bond-angle distributions and structural point defects. Calculations show that surface effects on surface static and thermodynamic properties of models are significant according to the change in the number of Al atoms in the surface layers. Evolution of the local environment of oxygen in the surface shell of nanoparticles upon cooling from the melt toward the glassy state was also found and discussed. In addition, the nanosize dependence of the glass transition temperature was presented.