Structure and dynamics of amorphous silica surfaces

Generation of Amorphous Silica Surfaces with Controlled Roughness

Amorphous silica (a-SiO2) surfaces, when grafted with select metals on the active sites of the functionalized surfaces, can act as useful heterogeneous catalysts. From a molecular modeling perspective, one challenge has been generating a-SiO2 slab models with controllable surface roughness to facilitate the study of the effect of surface morphology on the material properties. Previous computational methods either generate relatively flat surfaces or periodically corrugated surfaces that do not mimic the full range of potential surface roughness of the amorphous silica material. In this work, we present a new method, inspired by the capillary fluctuation theory of interfaces, in which rough silica slabs are generated by cleaving a bulk amorphous sample using a cleaving plane with Fourier components randomly generated from a Gaussian distribution. The width of this Gaussian distribution (and thus the degree of surface roughness) can be tuned by varying the surface roughness parameter α. Using the van Beest, Kramer, and van Santen (BKS) force field, we create a large number of silica slabs using cleaving surfaces of varying roughness (α) and using two different system sizes. These surfaces are then characterized to determine their roughness (mean-squared displacement), density profile, and ring size distribution. This analysis shows a higher concentration of surface defects (under-/overcoordinated atoms and strained rings) as the surface roughness increases. To examine the effect of the roughness on surface reactivity, we re-equilibriate a subset of these slabs using the reactive force field ReaxFF and then expose the slabs to water and observe the formation of surface silanols. We observe that the rougher surfaces exhibit higher silanol concentrations as well as bimodal acidity.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Surface Diffusion Is Controlled by Bulk Fragility across All Glass Types

Surface diffusion is vastly faster than bulk diffusion in some glasses, but only moderately enhanced in others. We show that this variation is closely linked to bulk fragility, a common measure of how quickly dynamics is excited when a glass is heated to become a liquid. In fragile molecular glasses, surface diffusion can be a factor of 10^{8} faster than bulk diffusion at the glass transition temperature, while in the strong system SiO_{2}, the enhancement is a factor of 10. Between these two extremes lie systems of intermediate fragility, including metallic glasses and amorphous selenium and silicon. This indicates that stronger liquids have greater resistance to dynamic excitation from bulk to surface and enables prediction of surface diffusion, surface crystallization, and formation of stable glasses by vapor deposition.

Engineering the Surface Pattern of Microparticles: From Raspberry-like to Golf Ball-like

Control of the shape and uniformity of colloid particles is essential for realizing their functionality in various applications. Herein, we report a facile approach for the synthesis of narrowly dispersed anisotropic microparticles with well-defined raspberry-like and golf ball-like surface patterns. First, we demonstrate that hybrid raspberry-like particles can be achieved through a one-pot polymerization method using glycidyl polyhedral oligomeric silsesquioxane (GPOSS) and pentaerythritol tetra(3-mercaptopropionate) (PETMP) as monomers. Varying the polymerization parameters such as catalyst loading, monomer concentration, and the molar ratio of monomers, we are able to regulate the sizes and surface protrusion numbers of these raspberry-like microparticles. The formation mechanism is attributed to a competition balance between thiol-epoxy reaction and thiol-thiol coupling reaction. The former promotes rapid formation of large core particles between PETMP and GPOSS droplets (which can serve as core particles), while the latter allows for generation of surface protrusions by PETMP self-polymerization, leading to the formation of raspberry-like surface patterns. Based on the different POSS contents in the surface protrusions and cores of the raspberry-like microparticles, we demonstrate that they can be used as precursors to produce microporous silica (sub)microparticles with golf ball-like morphology via pyrolysis subsequently. Overall, this work provides a facile yet controllable approach to synthesize narrowly dispersed anisotropic microparticles with diverse surface patterns.

Factors correlating to enhanced surface diffusion in metallic glasses

The enhancement of surface diffusion (D_S) over the bulk (D_V) in metallic glasses (MGs) is well documented and likely to strongly influence the properties of glasses grown by vapor deposition. Here, we use classical molecular dynamics (MD) simulations to identify different factors influencing the enhancement of surface diffusion in MGs. MGs have a simple atomic structure and belong to the category of moderately fragile glasses that undergo pronounced slowdown of bulk dynamics with cooling close to the glass transition temperature (T_g). We observe that D_S exhibits a much more moderate slowdown compared to D_V when approaching T_g, and D_S/D_V at T_g varies by two orders of magnitude among the MGs investigated. We demonstrate that both the surface energy and the fraction of missing bonds for surface atoms show good correlation to D_S/D_V, implying that the loss of nearest neighbors at the surface directly translates into higher mobility, unlike the behavior of network-bonded and hydrogen-bonded organic glasses. Fragility, a measure of the slowdown of bulk dynamics close to T_g, also correlates to D_S/D_V, with more fragile systems having larger surface enhancement of mobility. The deviations observed in the fragility-D_S/D_V relationship are shown to be correlated to the extent of segregation or depletion of the mobile element at the surface. Finally, we explore the relationship between the diffusion pre-exponential factor (D₀) and the activation energy (Q) and compare it to a ln(D₀)-Q correlation previously established for bulk glasses, demonstrating similar correlations from MD as in the experiments and that the surface and bulk have very similar ln(D₀)-Q correlations.

Structure and vibrational properties of sodium silicate glass surfaces

Using molecular dynamics simulations, we investigate how the structural and vibrational properties of the surfaces of sodo-silicate glasses depend on the sodium content as well as the nature of the surface. Two types of glass surfaces are considered: A melt-formed surface (MS) in which a liquid with a free surface has been cooled down into the glass phase and a fracture surface (FS) obtained by tensile loading of a glass sample. We find that the MS is more abundant in Na and non-bridging oxygen atoms than the FS and the bulk glass, whereas the FS has higher concentration of structural defects such as two-membered rings and under-coordinated Si than the MS. We associate these structural differences to the production histories of the glasses and the mobility of the Na ions. It is also found that for Na-poor systems, the fluctuations in composition and local atomic charge density decay with a power-law as a function of distance from the surface, while Na-rich systems show an exponential decay with a typical decay length of ≈2.3 Å. The vibrational density of states shows that the presence of the surfaces leads to a decrease in the characteristic frequencies in the system. The two-membered rings give rise to a pronounce band at ≈880 cm^-1, which is in good agreement with experimental observations.

Kinetic and Binding Studies Reveal Cooperativity and Off-Cycle Competition for H-Bonding Catalysis with Silsesquioxane Silanols

The catalytic activity, kinetics, and quantification of H-bonding ability of incompletely condensed polyhedral oligomeric silsesquioxane (POSS) silanols are reported. POSS-triols, a homogeneous model for vicinal silica surface sites, exhibit enhanced H-bonding compared with other silanols and alcohols as quantified using a ³¹ P NMR probe. Evaluation of a Friedel-Crafts addition reaction shows that phenyl-POSS-triol is active as an H-bond donor catalyst whereas other POSS silanols studied are not. An in-depth kinetic study (using RPKA and VTNA) highlights the concentration-dependent H-bonding behavior of POSS-triols, which is attributed to intermolecular association forming an off-cycle dimeric species. Binding constants provide additional support for reduced H-bond ability at higher concentrations, which is attributed to competitive association. POSS-triol self-association disrupts H-bond donor abilities relevant for catalysis by reducing the concentration of active monomeric catalyst.

The Relevance of Point Defects in Studying Silica-Based Materials from Bulk to Nanosystems

The macroscopic properties of silica can be modified by the presence of local microscopic modifications at the scale of the basic molecular units (point defects). Such defects can be generated during the production of glass, devices, or by the environments where the latter have to operate, impacting on the devices’ performance. For these reasons, the identification of defects, their generation processes, and the knowledge of their electrical and optical features are relevant for microelectronics and optoelectronics. The aim of this manuscript is to report some examples of how defects can be generated, how they can impact device performance, and how a defect species or a physical phenomenon that is a disadvantage in some fields can be used as an advantage in others.

Atomic-scale description of interfaces in rutile/sodium silicate glass-crystal composites

In this work interfaces between (Na2O)x(SiO2)1-x glasses (for x = 0.0, 0.1 and 0.2) and TiO2 crystals are simulated using molecular dynamics and empirical potentials. Interfaces are presented for the distinct terminating surfaces of TiO2 with Miller indices ≤2, the properties of which have been investigated using atomistic models. Simulations showed that partially ordered layers had been induced in the glass close to the interfaces, with successive oxygen-rich and cation-rich planes being noted. The first silicate layer in contact with the crystal tended to be highly-structured, with Si ions occupying well-defined positions that depend on the orientation of the crystal at the interface, and showing 2-dimensional ordering depending on glass composition. Finally, interface energies were calculated. These indicated that the interface formation may stabilise a crystal surface in comparison to maintaining a free surface. Results are presented suggesting that the structural flexibility of the glass network allows it to conform to the crystal, thereby providing charge compensation and avoiding large relaxation of the crystal structure close to the interfaces. Such interfacial properties could be crucial to improving phenomenological models of glass-crystal composite properties.

Molecular Dynamics Simulation of Water Confinement in Disordered Aluminosilicate Subnanopores

The porous structure and mass transport characteristics of disordered silicate porous media were investigated via a geometry based analysis of water confined in the pores. Disordered silicate porous media were constructed to mimic the dissolution behavior of an alkali aluminoborosilicate glass, i.e., soluble Na and B were removed from the bulk glass, and then water molecules and Na were introduced into the pores to provide a complex porous structure filled with water. This modelling approach revealed large surface areas of disordered porous media. In addition, a number of isolated water molecules were observed in the pores, despite accessible porous connectivity. As the fraction of mobile water was approximately 1%, the main water dynamics corresponded to vibrational motion in a confined space. This significantly reduced water mobility was due to strong hydrogen-bonding water-surface interactions resulting from the large surface area. This original approach provides a method for predicting the porous structure and water transport characteristics of disordered silicate porous media.

Modelling the surface of amorphous dehydroxylated silica: the influence of the potential on the nature and density of defects

The nature and density of defects on the amorphous dehydroxylated silica surface are studied by molecular dynamics for information on the silanol groups of pretreated silica.

Structural Origin of Enhanced Dynamics at the Surface of a Glassy Alloy

The enhancement of mobility at the surface of an amorphous alloy is studied using a combination of molecular dynamic simulations and normal mode analysis of the nonuniform distribution of Debye-Waller factors. The increased mobility at the surface is found to be associated with the appearance of Arrhenius temperature dependence. We show that the transverse Debye-Waller factor exhibits a peak at the surface. Over the accessible temperature range, we find that the bulk and surface diffusion coefficients obey the same empirical relationship with the respective Debye-Waller factors. Extrapolating this relationship to lower T, we argue that the observed decrease in the constraint at the surface is sufficient to account for the experimentally observed surface enhancement of mobility.

Influence of Hydrogen Bonding on the Surface Diffusion of Molecular Glasses: Comparison of Three Triazines

Surface grating decay measurements have been performed on three closely related molecular glasses to study the effect of intermolecular hydrogen bonds on surface diffusion. The three molecules are derivatives of bis(3,5-dimethyl-phenylamino)-1,3,5-triazine and differ only in the functional group R at the 2-position, with R being C₂H₅, OCH₃, and NHCH₃, and referred to as "Et", "OMe", and "NHMe", respectively. Of the three molecules, NHMe forms more extensive intermolecular hydrogen bonds than Et and OMe and was found to have slower surface diffusion. For Et and OMe, surface diffusion is so fast that it replaces viscous flow as the mechanism of surface grating decay as temperature is lowered. In contrast, no such transition was observed for NHMe under the same conditions, indicating significantly slower surface diffusion. This result is consistent with the previous finding that extensive intermolecular hydrogen bonds slow down surface diffusion in molecular glasses and is attributed to the persistence of hydrogen bonds even in the surface environment. This result is also consistent with the lower stability of the vapor-deposited glass of NHMe relative to those of Et and OMe and supports the view that surface mobility controls the stability of vapor-deposited glasses.

Non-monotonic effect of confinement on the glass transition

The relaxation dynamics of glass forming liquids and their structure are influenced in the vicinity of confining walls. This effect has mostly been observed to be a monotonic function of the slit width. Recently, a qualitatively new behaviour has been uncovered by Mittal and coworkers, who reported that the single particle dynamics in a hard-sphere fluid confined in a planar slit varies in a non-monotonic way as the slit width is decreased from five to roughly two particle diametres (Mittal et al 2008 Phys. Rev. Lett. 100 145901). In view of the great potential of this effect for applications in those fields of science and industry, where liquids occur under strong confinement (e.g. nano-technology), the number of researchers studying various aspects and consequences of this non-monotonic behaviour has been rapidly growing. This review aims at providing an overview of the research activity in this newly emerging field. We first briefly discuss how competing mechanisms such as packing effects and short-range attraction may lead to a non-monotonic glass transition scenario in the bulk. We then analyse confinement effects on the dynamics of fluids using a thermodynamic route which relates the single particle dynamics to the excess entropy. Moreover, relating the diffusive dynamics to the Widom's insertion probability, the oscillations of the local dynamics with density at moderate densities are fairly well described. At high densities belonging to the supercooled regime, however, this approach breaks down signaling the onset of strongly collective effects. Indeed, confinement introduces a new length scale which in the limit of high densities and small pore sizes competes with the short-range local order of the fluid. This gives rise to a non-monotonic dependence of the packing structure on confinement, with a corresponding effect on the dynamics of structural relaxation. This non-monotonic effect occurs also in the case of a cone-plate type channel, where the degree of confinement varies with distance from the apex. This is a very promising issue for future research with the possibility of uncovering the existence of alternating glassy and liquid-like domains.

Molecular dynamics study of coagulation in silica-nanocolloid-water-NaCl systems based on the atomistic model

In the present work, large scale molecular dynamics (MD) simulations of nanocolloidal silica in aqueous NaCl solutions were performed using a fully atomistic model to study the microscopic structures and dynamics of the systems that lead to aggregation or gelation. Our attention is focused on the self-organizations that occur in the structures of the colloidal silica and water for various concentrations of NaCl. As the salt concentration increased, coagulation developed through the direct bonding of SiO4 units. The trend was explained by the systematic changes in the pair correlation functions related to the barrier height in the potential of mean force [J. G. Kirkwood, J. Chem. Phys., 1935, 3, 300]. Network structures of silica were visualised, and their fractal dimensions were examined by computing the running coordination numbers of Si-Si pairs and also by the analysis of two dimensional images. The calculated dimension by the former method was comparable to the experimental observations for the aggregation of silica colloids, and at longer length scales, super-aggregation was evident in the gelation process. The result from the 2D images is found to be insensitive to the differences in the structure. Clear changes in both the structure and mobility of the water were observed as the NaCl concentration increased, suggesting the importance of the solvent structures to these processes, although such a feature is lacking in the conventional models and most simulations of colloids.

Brittle to ductile transition in densified silica glass

Current understanding of the brittleness of glass is limited by our poor understanding and control over the microscopic structure. In this study, we used a pressure quenching route to tune the structure of silica glass in a controllable manner, and observed a systematic increase in ductility in samples quenched under increasingly higher pressure. The brittle to ductile transition in densified silica glass can be attributed to the critical role of 5-fold Si coordination defects (bonded to 5 O neighbors) in facilitating shear deformation and in dissipating energy by converting back to the 4-fold coordination state during deformation. As an archetypal glass former and one of the most abundant minerals in the Earth's crest, a fundamental understanding of the microscopic structure underpinning the ductility of silica glass will not only pave the way toward rational design of strong glasses, but also advance our knowledge of the geological processes in the Earth's interior.

Probing of 2 dimensional confinement-induced structural transitions in amorphous oxide thin film

Whereas the atomic structure of surface of crystals is known to be distinct from that of bulk, experimental evidence for thickness-induced structural transitions in amorphous oxides is lacking. We report the NMR result for amorphous alumina with varying thickness from bulk up to 5 nm, revealing the nature of structural transitions near amorphous oxide surfaces/interfaces. The coordination environments in the confined amorphous alumina thin film are distinct from those of bulk, highlighted by a decrease in the fractions of high-energy clusters (and thus the degree of disorder) with thickness. The result implies that a wide range of variations in amorphous structures may be identified by controlling its dimensionality.

Glassy nature of water in an ultraconfining disordered material: the case of calcium-silicate-hydrate

We present the structural and dynamic nature of water ultraconfined in the quasi-two-dimensional nanopores of the highly disordered calcium-silicate-hydrate (C-S-H), the major binding phase in cement. Our approach is based on classical molecular simulations. We demonstrate that the C-S-H nanopore space is hydrophilic, particularly because of the nonbridging oxygen atoms on the disordered silicate chains which serve as hydrogen-bond acceptor sites, directionally orienting the hydrogen atoms of the interfacial water molecules toward the calcium-silicate layers. The water in this interlayer space adopts a unique multirange structure: a distorted tetrahedral coordination at short range up to 2.7 Å, a disordered structure similar to that of dense fluids and supercooled phases at intermediate range up to 4.2 Å, and persisting spatial correlations through dipole-dipole interactions up to 10 Å. A three-stage dynamics governs the mean square displacement (MSD) of water molecules, with a clear cage stage characteristic of the dynamics in supercooled liquids and glasses, consistent with its intermediate-range structure identified here. At the intermediate time scales corresponding to the β-relaxation of glassy materials, coincident with the cage stage in MSD, the non-Gaussian parameter indicates a significant heterogeneity in the translational dynamics. This dynamic heterogeneity is induced primarily because of the heterogeneity in the distribution of hydrogen bond strengths. The strongly attractive interactions of water molecules with the calcium silicate walls serve to constrain their motion. Our findings have important implications on describing the cohesion and mechanical behavior of cement from its setting to its aging.

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.

Mechanical heterogeneity in nanoscale films of liquid silica

Molecular dynamics simulations are carried out for slabs of silica liquid with thicknesses between 1 and 3 nm . A local analysis of the Born contribution to the elastic modulus, CB, shows that the elasticity is not uniform throughout the slabs--CB is identical to that of bulk silica in the slab interior, but CB is larger at the slab edges. The larger CB at the slab edges is due to a distinct atomic level structure characterized by larger density, larger concentration of more highly coordinated ions, and smaller silica rings.

The glass transition and thermodynamics of liquid and amorphous TiO(2) nanoparticles

The glass transition and thermodynamics of spherical liquid TiO(2) nanoparticles, with different sizes ranging from 2 to 5 nm, have been studied in a model under non-periodic boundary conditions. We use the pairwise interatomic potentials proposed by Matsui and Akaogi. Models have been obtained by cooling from the melt via molecular dynamics (MD) simulation. The structural properties of liquid nanoparticles at 3500 K have been analyzed in detail through the partial radial distribution functions (PRDFs), coordination number distributions, bond-angle distributions and interatomic distances. Moreover, we also show the radial density profile in nanoparticles. Calculations show that size effects on the structure of a model are significant and that liquid TiO(2) nanoparticles have a distorted pentahedral network structure with the mean coordination numbers Z(Ti-O)≈5.0 and Z(O-Ti)≈2.5, while amorphous TiO(2) nanoparticles have an octahedral network structure. The temperature dependence of the surface structure and surface energy of the nanoparticles has been obtained and is presented. In addition, the size dependence of the glass transition temperature and the temperature dependence of the diffusion constant of atomic species have been found and are discussed.

Estimating thermal conductivity of amorphous silica nanoparticles and nanowires using molecular dynamics simulations

In several recent applications, including those aimed at developing thermal interface materials, nanoparticulate systems have been proposed to improve the effective behavior of the system. While nanoparticles by themselves may have low conductivities relative to larger particles owing to interfacial resistance, their use along with larger particles is believed to enhance the percolation threshold leading to better effective behavior overall. One critical challenge in using nanoparticulate systems is the lack of knowledge regarding their thermal conductivity. In this paper, the thermal conductivity of silica clusters (or nanoparticles) as well as nanowires is determined using molecular dynamics (MD) simulations. The equilibrium MD simulations of nanoparticles using Green-Kubo relations are demonstrated to be computationally very expensive and unsuitable for such nanoscaled systems. A nonequilibrium MD method adapted from the study of Müller-Plathe is shown to be faster and more accurate. The method is first demonstrated on bulk amorphous silica (using both cubic and orthorhombic simulation cells) and silica nanowires. The thermal conductivity values are compared to those reported in the literature. The mean thermal conductivity values for bulk silica and silica nanowire were estimated to be 1.2 W/mK and 1.435 W/mK, respectively. To model nanoparticles, the Müller-Plathe technique is adapted by dividing the cluster into concentric shells so as to capture the naturally radial mode of heat transfer. The mean thermal conductivity value of a 600-atom silica nanoparticle obtained using this approach was 0.589 W/mK. This value is approximately 50-60% lower than those of bulk silica or silica nanowire.

Molecular dynamics simulation of amorphous SiO2 nanoparticles

Molecular dynamics simulation of amorphous SiO2 spherical nanoparticles has been carried out in a model with different sizes, 2, 4, and 6 nm, under non-periodic boundary conditions. We use the pair interatomic potentials which have weak Coulomb interaction and Morse type short-range interaction. Models have been obtained by cooling from the melt via molecular dynamics (MD) simulation. Structural properties of amorphous nanoparticles obtained at 350 K have been studied via partial radial distribution functions (PRDFs), mean interatomic distances, coordination numbers, and bond-angle distributions, which are compared with those observed in the bulk. Calculations of the radial density profile in nanoparticles show the tendency of oxygen to concentrate at the surface as observed previously in other amorphous clusters or thin films. Size effects on structure of nanosized models are significant. The calculations show that if the size is larger than 4 nm, amorphous SiO2 nanoparticles have a distorted tetrahedral network structure with the mean coordination number ZSi-O approximately 4.0 and ZO-Si approximately 2.0 like those observed in the bulk. Surface structure, surface energy, and glass transition temperature of SiO2 nanoparticles have been obtained and presented.

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.

Model for the Water−Amorphous Silica Interface: The Undissociated Surface

The physical and chemical properties of the amorphous silica-water interface are of crucial importance for a fundamental understanding of electrochemical and electrokinetic phenomena, and for various applications including chromatography, sensors, metal ion extraction, and the construction of micro- and nanoscale devices. A model for the undissociated amorphous silica-water interface reported here is a step toward a practical microscopic model of this important system. We have extended the popular BKS and SPC/E models for bulk silica and water to describe the hydrated, hydroxylated amorphous silica surface. The parameters of our model were determined using ab initio quantum chemical studies on small fragments. Our model will be useful in empirical potential studies, and as a starting point for ab initio molecular dynamics calculations. At this stage, we present a model for the undissociated surface. Our calculated value for the heat of immersion, 0.3 J x m(-2), falls within the range of reported experimental values of 0.2-0.8 J x m(-2). We also study the perturbation of water properties near the silica-water interface. The disordered surface is characterized by regions that are hydrophilic and hydrophobic, depending on the statistical variations in silanol group density.

Density functional theory model of adsorption on amorphous and microporous silica materials

We present a novel quenched solid density functional theory (QSDFT) model of adsorption on heterogeneous surfaces and porous solids, which accounts for the effects of surface roughness and microporosity. Within QSDFT, solid atoms are considered as quenched component(s) of the solid-fluid system with given density distribution(s). Solid-fluid intermolecular interactions are split into hard-sphere repulsive and mean-field attractive parts. The former are treated with the multicomponent fundamental measure density functional. Capabilities of QSDFT are demonstrated by drawing on the example of adsorption on amorphous silica materials. We show that, using established intermolecular potentials and a realistic model for silica surfaces, QSDFT quantitatively describes adsorption/desorption isotherms of Ar and Kr on reference MCM-41, SBA-15, and LiChrosphere materials in a wide range of relative pressures. QSDFT offers a systematic approach to the practical problems of characterization of microporous, mesoporous, and amorphous silica materials, including an assessment of microporosity, surface roughness, and adsorption deformation. Predictions for the pore diameter and the extent of pore surface roughness in MCM-41 and SBA-15 materials are in very good agreement with recent X-ray diffraction studies.

Large melting-point hysteresis of Ge nanocrystals embedded in SiO2

The melting behavior of Ge nanocrystals embedded within SiO2 is evaluated using in situ transmission electron microscopy. The observed melting-point hysteresis is large (+/-17%) and nearly symmetric about the bulk melting point. This hysteresis is modeled successfully using classical nucleation theory without the need to invoke epitaxy.

Silica Nanoarchitectures with Tailored Pores Based on the Hybrid Three- and Four-Membered Rings

Inspired by the recent developments in the controlled synthesis of porous materials, we present herein the structural prediction of silica nanoarchitectures by using the three- (3MRs) and four-membered rings (4MRs), which are more frequently found in the nanometer-sized particles than in the bulk form, as building blocks. The proposed models include the active molecular rings, thin nanowires, hollow nanotubes, discrete fullerene-like cages, and porous zeolite-like three-dimensional networks. Their geometrical and electronic structures and properties were studied by performing density functional calculations. These silica nanostructures were proved, using molecular dynamics simulations, to possess intrinsic structural stabilities with highly symmetrical geometries and regular nanochannels. These atomically well-defined clusters, (SiO)(n), are chemically more reactive than those proposed earlier and are energetically more favorable for n > 20 in high-level density functional calculations over the corresponding two-membered ring (2MR) chains and rings as well as the pure 3MR networks. The nanoparticles and nanodevices based on them are expected to have potential technological applications that mainly make use of their characteristic geometrical structures (nanosized pores) and novel electronic properties.

Structural Model of Silica Nanowire Assembled from a Highly Stable (SiO2)8 Unit

The ground-state structures of silica clusters (SiO2)n for n = 1-8 were studied by performing calculations at the B3LYP/6-311+G(d) level of density functional theory. The results indicate that the growth mode of a silica nanowire based on small silica clusters may change at different wire lengths. A linear chain might be assembled from the smallest clusters of rhombic two-membered ring (2MR) with n < or = 5, while the growth motif changes at n = 6 into a more compact form composed of three-membered-rings (3MRs). The 3MR-containing structures become energetically favorable configurations for even longer silica clusters. In particular, the closed molecular ring consisting of 3MRs at n = 8 (i.e., (SiO2)8) with a high symmetry shows extreme energetic stability and relatively high chemical reactivity and thus is considered to be an important building block to assemble into silica nanowires. The relative stability of so-assembled silica nanowires were evaluated and compared with the models of silica nanowires in the literature.

Theoretical Study of the Interaction between Selected Adhesives and Oxide Surfaces

We investigate the competition of the various organic components of two representative adhesive systems for reactive defect sites at model surfaces of both SiO2 and Al2O3. The reaction energies of resin monomers, curing agents, and in some cases also of additional adhesion promoters with the defects are calculated. We applied a density-functional based tight-binding method including a self-consistent correction of the Mulliken charges, which has already proven to be a useful tool for computational materials science, delivering reliable structural and energetic information.