Polarization interactions and boroxol ring formation in boron oxide: A molecular dynamics study

Many-body effects at the origin of structural transitions in B2O3

The structural properties of glassy diboron trioxide, g-B₂O₃, are investigated from ambient to high pressure conditions using two types of atomic force-field models that account for many-body effects. These models are parameterized by a dipole- and force-fitting procedure of reference datasets created via first-principles calculations on a series of configurations. The predictions of the models are tested against experimental data, where particular attention is paid to the structural transitions in g-B₂O₃ that involve changes to both the short- and medium-range order. The models outperform those previously devised, where improvement originates from the incorporation of two key physical ingredients, namely, (i) the polarizability of the oxide ion and (ii) the ability of an oxide ion to change both size and shape in response to its coordination environment. The results highlight the importance of many-body effects for accurately modeling this challenging system.

Direct Experimental Evidence for Abundant BO4-BO4 Motifs in Borosilicate Glasses From Double-Quantum 11B NMR Spectroscopy

By using double-quantum-single-quantum correlation ¹¹B nuclear magnetic resonance (NMR) experiments and atomistic molecular dynamics (MD) simulations, we resolve the long-standing controversy of whether directly interlinked BO₄-BO₄ groups exist in the technologically ubiquitous class of alkali/alkaline-earth based borosilicate (BS) glasses. Most structural models of Na₂O-B₂O₃-SiO₂ glasses assume the absence of B^[4]-O-B^[4] linkages, whereas they have been suggested to exist in Ca-bearing BS analogs. Our results demonstrate that while B^[4]-O-B^[4] linkages are disfavored relative to their B^[3]-O-B^[3]/B^[4] counterparts, they are nevertheless abundant motifs in Na₂O-B₂O₃-SiO₂ glasses over a large composition space, while the B^[4]-O-B^[4] contents are indeed elevated in Na₂O-CaO-B₂O₃-SiO₂ glasses. We discuss the compositional and structural parameters that control the degree of B^[4]-O-B^[4] bonding.

Structure–composition trends in multicomponent borosilicate-based glasses deduced from molecular dynamics simulations with improved B–O and P–O force fields

Improved B–O and P–O force fields provide accurate molecular dynamics simulations of multicomponent glasses.

Medium-Range Structural Organization of Phosphorus-Bearing Borosilicate Glasses Revealed by Advanced Solid-State NMR Experiments and MD Simulations: Consequences of B/Si Substitutions

The short and intermediate range structures of a large series of bioactive borophosphosilicate (BPS) glasses were probed by solid-state nuclear magnetic resonance (NMR) spectroscopy and atomistic molecular dynamics (MD) simulations. Two BPS glass series were designed by gradually substituting SiO₂ by B₂O₃ in the respective phosphosilicate base compositions 24.1Na₂O-23.3CaO-48.6SiO₂-4.0P₂O₅ ("S49") and 24.6Na₂O-26.7CaO-46.1SiO₂-2.6P₂O₅ ("S46"), the latter constituting the "45S5 Bioglass" utilized for bone grafting applications. The BPS glass networks are built by interconnected SiO₄, BO₄, and BO₃ moieties, whereas P exists mainly as orthophosphate anions, except for a minor network-associated portion involving P-O-Si and P-O-B^[4] motifs, whose populations were estimated by heteronuclear ³¹P{¹¹B} NMR experimentation. The high Na⁺/Ca²⁺ contents give fragmented glass networks with large amounts of nonbridging oxygen (NBO) anions. The MD-generated glass models reveal an increasing propensity for NBO accommodation among the network units according to BO₄ < SiO₄ < BO₃ ≪ PO₄. The BO₄/BO₃ intermixing was examined by double-quantum-single-quantum correlation ¹¹B NMR experiments, which evidenced the presence of all three BO₃-BO₃, BO₃-BO₄, and BO₄-BO₄ connectivities, with B^[3]-O-B^[4] bridges dominating. Notwithstanding that B^[4]-O-B^[4] linkages are disfavored, both NMR spectroscopy and MD simulations established their presence in these modifier-rich BPS glasses, along with non-negligible B^[4]-NBO contacts, at odds with the conventional structural view of borosilicate glasses. We discuss the relative propensities for intermixing of the Si/B/P network formers. Despite the absence of pronounced preferences for Si-O-Si bond formation, the glass models manifest subtle subnanometer-sized structural inhomogeneities, where SiO₄ tetrahedra tend to self-associate into small chain/ring motifs embedded in BO₃/BO₄-dominated domains.

Liquid B2O3 up to 1700 K: x-ray diffraction and boroxol ring dissolution

Using high energy x-ray diffraction, the structure factors of glassy and molten B2O3 were measured with high signal-to-noise, up to a temperature of T  =  1710(20) K. The observed systematic changes with T are shown to be consistent with the dissolution of hexagonal [B3O6] boroxol rings, which are abundant in the glass, whilst the high-T (>~1500 K) liquid can be more closely described as a random network structure based on [BO3] triangular building blocks. We therefore argue that diffraction data are in fact qualitatively sensitive to the presence of small rings, and support the existence of a continuous structural transition in molten B2O3, for which the temperature evolution of the 808 cm−1 Raman scattering band (boroxol breathing mode) has long stood as the most emphatic evidence. Our conclusions are supported by both first-principles and polarizable ion model molecular dynamics simulations which are capable of giving good account of the experimental data, so long as steps are taken to ensure a ring fraction similar to that expected from Raman spectroscopy. The mean thermal expansion of the B-O bond has been measured directly to be αBO  =  3.7(2)  ×  10−6 K−1, which accounts for a few percent of the bulk expansion just above the glass transition temperature, but accounts for greater than one third of the bulk expansion at temperatures in excess of 1673 K.

Study on the polymer translocation induced blockade ionic current inside a nanopore by Langevin dynamics simulation

The blockade ionic current inside a nanopore due to polymer translocation is studied using a three-dimensional Langevin dynamics method. The blockade current IB is dependent on the polymer length N, polymer configuration, polymer-pore interaction, and charge of the polymer. The behavior of IB can be explained using four factors: (1) the volume vacancy fraction fV inside the pore; (2) the conformation of the polymer; (3) the location of the polymer inside the pore; and (4) the total charge Ztot inside the pore. We find that IB increases with fV but decreases with increasing |Ztot|. The influence of the polymer's conformation is complex, dependent on the size of polymer RG and the cross-sectional size of the pore s. A compact conformation can decrease IB when RG > s but increase IB when RG < s. For the latter case, the conformation of the polymer is too small to block the pore, thus providing a broad passage for the ions. At the same fV, monomers will locate close to the surface with a large polymer-pore attraction, which also provides a large IB.

Ionic Current Rectification Through Silica Nanopores

Nanopores immersed in electrolytic solution and under the influence of an electric field can produce ionic current rectification, where ionic currents are higher for one voltage polarity than for the opposite polarity, resulting in an asymmetric current-voltage (I-V) curve. This behavior has been observed in polymer and silicon-based nanopores as well as in theoretically studied continuum models. By means of atomic level molecular dynamics (MD) simulations, we have performed a systematic investigation of KCl conductance in silica nanopores with a total simulation time of 680 ns. We found that ion-binding spots at the silica surfaces, such as dangling atoms, have effects on the ion concentration and electrostatic potential inside the nanopore, producing asymmetric I-V curves. Conversely, silica surfaces without ion-binding spots produce symmetric I-V curves.

Boroxol Rings in Liquid and Vitreous B2O3 from First Principles

We investigate the structural and vibrational properties of glassy B2O3 using first-principles molecular dynamics simulations. In particular, we determine the boroxol rings fraction f for which there is still no consensus in the literature. Two numerical models containing either a low or a high level of boroxol rings are tested against a gamut of experimental probes (static structure factor, Raman, 11B and 17O NMR data). We show that only the boroxol-rich model (f=75%) can reproduce the full set of observables. Total-energy calculations show that at the glass density, boroxol-rich structures are favored by about 6 kcal/(mol boroxol). Finally, the liquid state is explored in the 2,000-4,000 K range and a reduction of f to 10%-20% is obtained.

Structure and stability of small boron and boron oxide clusters

To rationally design and explore a potential energy source based on the highly exothermic oxidation of boron, density functional theory (DFT) was used to characterize small boron clusters with 0-3 oxygen atoms and a total of up to ten atoms. The structures, vibrational frequencies, and stabilities were calculated for each of these clusters. A quantum molecular dynamics procedure was used to locate the global minimum for each species, which proved to be crucial given the unintuitive structure of many of the most stable isomers. Additionally, due to the plane-wave, periodic DFT code used in this study a straightforward comparison of these clusters to the bulk boron and B2O3 structures was possible despite the great structural and energetic differences between the two forms. Through evaluation of previous computational and experimental work, the relevant low-energy structures of all but one of the pure boron clusters can be assigned with great certainty. Nearly all of the boron oxide clusters are described here for the first time, but there are strong indications that the DFT procedure chosen is particularly well suited for the task. Insight into the trends in boron and boron oxide cluster stabilities, as well as the ultimate limits of growth for each, are also provided. The work reported herein provides crucial information towards understanding the oxidation of boron at a molecular level.

Water-silica force field for simulating nanodevices

Amorphous silica is an inorganic material that is central for many nanotechnology applications, such as nanoelectronics, microfluidics, and nanopore sensors. To use molecular dynamics (MD) simulations to study the behavior of biomolecules interacting with silica, we developed a force field for amorphous silica surfaces based on their macroscopic wetting properties that is compatible with the CHARMM force field and TIP3P water model. The contact angle of a water droplet on a silica surface served as a criterion to tune the intermolecular interactions. The resulting force field was used to study the permeation of water through silica nanopores, illustrating the influence of the surface topography and the intermolecular parameters on permeation kinetics. We find that minute modeling of the amorphous surface is critical for MD studies, since the particular arrangement of surface atoms controls sensitively electrostatic interactions between silica and water.

Composition and temperature dependence of cesium-borate glasses by molecular dynamics

The structural aspects of xCs2O-(1-x)B2O3 glasses have been investigated by molecular dynamics as functions of Cs2O content (x=0.2, 0.3, and 0.4) and temperature (T=300 and 1250 K). The tetrahedral (BØ4-) and triangular (BØ3,BØ2O-, and BØO2 (2-)) short-range order borate units were found to be the structure-building entities of the simulated glasses [Ø=bridging oxygen (BO) and O-=nonbridging oxygen (NBO) atom]. The increase of Cs2O content results in the progressive increase of the NBO-containing triangle population at the expense of the BO4- tetrahedral units. The same effect is caused by temperature increase at a fixed Cs2O content, and this was associated with the "fragile" characteristics of alkali borate glasses. A comparison of simulated Cs and Li borates showed very similar structures at x=0.2, but dissimilar ones when the alkali content exceeds this composition. In particular, for x>0.2 Cs borates exhibit a preference for NBO formation relative to Li borates. Differences in the microstructure of sites hosting Cs ions were found, and this permits their classification into bridging (b type) and nonbridging type (nb type) of sites. b-type sites consist exclusively of BO atoms, while both BO and NBO atoms participate in nb-type sites. These differences in Cs-site local bonding characteristics were found to be reflected on the Cs-O(site) vibration frequencies. Also, the computed Cs-O vibrational responses for simulated Cs borates were found to compare well with experimental far-infrared spectra.

A molecular dynamics study of the structural dependence of boron oxide nanoparticles on shape

Molecular dynamics simulation is employed to study the effect of varying nanoparticle shape on the structure of boron oxide nanoparticles. Two nanoshapes are investigated and compared: a sphere of diameter 16 A and a cube of dimension 16 x 16 x 16 A. A many-body polarization model is employed within the simulation, accounting for dipole moments induced by local electric fields. The resulting network is described by a short-range structure consisting of planar BO(3) units, while the intermediate-range structure is described by six-membered planar boroxol rings. Both the fraction of boroxol rings and their locations differ between the two nanoshapes. All planar boroxol rings within the spherical simulation are located on the interior, while planar rings within the cubic simulation aggregate to the cube walls. In addition, structural differences appear between the two shapes at longer ranges, including the formation of "layers" aligned parallel to the walls of the cube, reminiscent of both the low-density crystalline phase and the high-density amorphous form of boron oxide.

A molecular interpretation of vitreous boron oxide dynamics

The mobility of vitreous boron oxide is studied by molecular dynamics simulation. A polarization model that incorporates induced dipoles arising both from charges and from other induced dipoles on atoms with nonzero polarizability is used to simulate boron oxide glass at various temperatures above the glass transition temperature. Particle mobility is investigated through the calculation of the self-intermediate scattering function and the mean-squared displacement. The calculations clearly reveal a two-step relaxation with a plateau at intermediate times for all investigated temperatures. With respect to atomic species, boron atoms are less mobile than oxygen atoms at all temperatures within the plateau region. Through analyzing particle trajectories, it is revealed that BO(3) groups move as one unit and follow each other in a stringlike manner. Three connected BO(3) groups comprise a six-membered boroxol ring, which is shown to move in a collective manner, requiring the simultaneous movement of all ring atoms. The boroxol ring is observed to be confined, or caged, during the plateau region, and jumps to a new location at longer times. This observation is linked to the concept of strong versus fragile glass formers and the potential energy landscape. In addition to the caging feature, an overshoot or dip occurs in the plateau regions of the mean-squared displacement and self-intermediate scattering functions respectively. These features are followed by a ringing pattern, previously associated with finite size effects in other strong glass formers, which persist for the duration of the plateau region. Both features are shown to be consistent with the bending of atomic "cages" from the plane of the boroxol ring, and arise due to the displacement of atoms from local minimum energy configurations.