Coordination Number Distribution Research Articles

Comparison of intermediate-range order in GeO2 glass: Molecular dynamics using machine-learning interatomic potential vs reverse Monte Carlo fitting to experimental data.

The short-range order and intermediate-range order in GeO2 glass are investigated by molecular dynamics using machine-learning interatomic potential trained on abinitio calculation data and compared with the reverse Monte Carlo fitting of neutron diffraction data. To characterize the structural differences in each model, the total/partial structure factors, coordination number, ring size and shape distributions, and persistent homology analysis were performed. These results show that although the two approaches yield similar two-body correlations, they can lead to three-dimensional models with different short- and intermediate-range ordering. A clear difference was observed especially in the ring distributions; RMC models exhibit a broad distribution in the ring size distribution, while neural network potential molecular dynamics yield much narrower ring distributions. This confirms that the density functional approximation in the abinitio calculations determines the preferred network assembly more strictly than RMC with simple coordination constraints even when using multiple diffraction data.

Exploring the role of sample size on fracture growth mechanisms in intact rock: Insights from 3D DEM-DFN analysis

The effect of sample size on the deformation characteristic and fracture growth mechanism of intact rock subjected to unconfined compressive stress state is examined in this study via discrete element method coupled with discrete fracture network. Three-dimensional numerical models are first developed based on reported laboratory experiments and verified to realistically replicate the effect of sample size on the macro-mechanical response of the intact rock. This includes strength behavior, fracturing activities and energy budgets. Micromechanical analyses are then performed to understand the role of sample size on the deformational behavior and damage progression in intact rock. Emphasis is placed on the distributions of coordination number, evolutions of crack density and the degree of crack anisotropy with regard to invariants of crack tensors. Results reveal that deformability and strength of intact rock are vastly reliant on the sample size. Increasing sample size facilitate the accumulation of induced microcracks due to the increased probability of interparticle bond failure. It is found that the increase in sample size can result in more tensile openings with strong dilatancy. Finally, the monotonical increase in seismic b-values with the increase in sample size suggests an acceleration in the number of small-magnitude AE events.

An AIMD investigation on the structural properties of Fe-SiO2 melt for converting SiO2 inclusions into the reinforcing particles in Fe matrix

It is essential to understand the structural properties of Fe-SiO2 melt for converting SiO2 inclusions into the reinforcing particles in Fe matrix by liquid metallurgical synthesis. In this work, the structural behaviors of Si and O in Fe27SiO2, Fe24(SiO2)2, Fe21(SiO2)3, and Fe18(SiO2)4 melts at 2800 K and the formation of SiOn clusters in Fe24(SiO2)2 melt (20 at.% SiO2) from 2800 K to 2000 K are investigated by ab-initio molecular dynamics simulations. According to the pseudo-binary mixing enthalpies, partial pair correlation functions, and chemical short-range order, the critical composition range for the formation of SiOn clusters in Fe-SiO2 melt at 2800 K is from 20 at.% SiO2 to 30 at.% SiO2, and the critical temperature range in the quenching process of Fe24(SiO2)2 melt is from 2600 K to 2400 K. By further analyzing the distribution of partial coordination numbers and bond angles, the quenching process of Fe24(SiO2)2 melt from 2800 K to 2000 K can be divided into three stages: the homogenous Fe-SiO2 melt from 2800 K to 2600 K, the dispersive formation of SiOn clusters containing a considerable number of SiO4 tetrahedrons from 2600 K to 2200 K, the gathering of SiOn clusters contributed by the connection of SiO4 tetrahedrons in the pattern of OSi2 from 2200 K to 2000 K.

Structural and Dynamical Effects of the CaO/SrO Substitution in Bioactive Glasses.

Silicate glasses containing silicon, sodium, phosphorous, and calcium have the ability to promote bone regeneration and biodegrade as new tissue is generated. Recently, it has been suggested that adding SrO can benefit tissue growth and silicate glass dissolution. Motivated by these recent developments, the effect of SrO/CaO-CaO/SrO substitution on the local structure and dynamics of Si-Na-P-Ca-O oxide glasses has been studied in this work. Differential thermal analysis has been performed to determine the thermal stability of the glasses after the addition of strontium. The local structure has been studied by neutron diffraction augmented by Reverse Monte Carlo simulation, and the local dynamics by neutron Compton scattering and Raman spectroscopy. Differential thermal analysis has shown that SrO-containing glasses have lower glass transition, melting, and crystallisation temperatures. Moreover, the addition of the Sr2+ ions decreased the thermal stability of the glass structure. The total neutron diffraction augmented by the RMC simulation revealed that Sr played a similar role as Ca in the glass structure when substituted on a molar basis. The bond length and the coordination number distributions of the network modifiers and network formers did not change when SrO (x = 0.125, 0.25) was substituted for CaO (25-x). However, the network connectivity increased in glass with 12.5 mol% CaO due to the increased length of the Si-O-Si interconnected chain. The analysis of Raman spectra revealed that substituting CaO with SrO in the glass structure dramatically enhances the intensity of the high-frequency band of 1110-2000 cm-1. For all glasses under investigation, the changes in the relative intensities of Raman bands and the distributions of the bond lengths and coordination numbers upon the SrO substitution were correlated with the values of the widths of nuclear momentum distributions of Si, Na, P, Ca, O, and Sr. The widths of nuclear momentum distributions were observed to soften compared to the values observed and simulated in their parent metal-oxide crystals. The widths of nuclear momentum distributions, obtained from fitting the experimental data to neutron Compton spectra, were related to the amount of disorder of effective force constants acting on individual atomic species in the glasses.

MOLECULAR DYNAMICS SIMULATION OF STRUCTURAL PROPERTIES OF Al<sub>2</sub>O<sub>3</sub>-SiO<sub>2</sub>-CaO OXIDE DURING THE COOLING PROCESS

The change of structural properties of Al2O3-SiO2-CaO oxide during the cooling process was studied using the molecular dynamic simulation. The microstructures of the sample were investigated through the radial distribution function (RDF), coordination number (CN), and bond angle (BA) distribution. The results show a structural change from the liquid state to the amorphous state when the sample is cooled from 4000 K to 300 K. The glass transition occurs at Tg = 1585 K. The fraction of AlO4, CaO6, and SiO4 units increases with decreasing temperature and dominates at 300 K.

Comparing the MP2, DKH2‐MP2 and SCAN Hybrid Forces Molecular Dynamics Simulation on the Structure and Dynamics of Hydrated Ion Ba2+ in Aqueous Solution

AbstractThe study investigated the structural and dynamical properties of hydrated ion Barium (Ba2+) in aqueous solution using quantum mechanical/molecular mechanical molecular dynamics simulations (QM/MM MD simulation). The primary objective was to assess the performance of the SCAN functional in describing the structure and dynamics of formed Ba2+ hydration shells in aqueous solution. The SCAN/MM MD simulation demonstrated the closest agreement with previously reported experimental data regarding the structural and dynamical properties of hydrated ion Ba2+. Specifically, the Ba2+ and surrounding water molecules distance was determined to be 2.79 Å, with equal distribution of coordination numbers of 8 and 9 in the first hydration shell. Moreover, the stretching frequency of 232.21 cm−1 closely matched the CCSD(T) results. The obtained results from SCAN/MM MD has a better agreement than results from MP2/MM and DKH2 MP2/MM molecular dynamics simulation. These findings underscore the potential of SCAN/MM MD simulation as a promising method for investigating the structural and dynamical properties of hydrated ions.

Effect of Cu addition on the properties of Pd-based metallic glasses by experimentation and molecular dynamics simulations

The effect of trace Cu element additions on the properties of Pd-based amorphous alloys was investigated by a combination of experimental and molecular dynamics (MD) simulations. Thin ribbons of Pd40Ni40P20 and containing different proportions of Cu (1 %, 3 %, and 5 %) were prepared, analyzed by X-ray diffraction (XRD) to determine the phase structure, tested for their thermal properties by differential scanning calorimetry (DSC), and nanoindentation to test their mechanical properties. Meanwhile, the amorphization process and compression behavior of these alloys were simulated by MD simulation, and the microstructures of the alloys were resolved by using the radial distribution function (RDF) and the coordination number distribution function (PDF), and the glass transition temperatures (Tg) were deduced from the energy variations. The Voronoi polyhedral analyses showed that the addition of trace amounts of Cu altered the content of the specific clusters, which affected the macroscopic mechanical properties. Both nanoindentation experimental and simulation results confirmed that the addition of trace Cu reduces the Tg and improved the mechanical properties of the Pd40Ni40P20 alloy while maintaining the amorphous structure.

Effect of basicity on the structure and viscosity properties of HIsmelt slag: A molecular dynamics simulation

ABSTRACT HIsmelt is an emerging and potentially promising non-blast furnace ironmaking process, and the characteristics of its slag are crucial to the smelting process. Currently, there is no atomic-scale research on the slag of the HIsmelt process. In this work, molecular dynamics simulations were employed to investigate the effect of basicity on the structure and viscosity properties of the HIsmelt CaO-SiO2-Al2O3-FeO slag system at 1773 K. The slag structure was characterised by calculating structural parameters, including radial distribution function, coordination number, oxygen type, and bond angle distribution. Furthermore, the viscosity of the slag system was estimated by using the self-diffusion coefficient and compared with mathematical models. The results indicate that basicity has a negligible impact on the short-range ordering of aluminosilicates in the HIsmelt slag system. As slag basicity increases from 0.4 to 1.8, the concentration of bridging oxygen decreases, non-bridging oxygen increases, and the [SiO4]4+ and [AlO4]5+ network structures undergo depolymerisation in the system. Additionally, the self-diffusion coefficient of atoms increases, and the slag viscosity exhibits a declining trend with increasing basicity. In the basicity range of 0.8 −1.8, the current MD simulations align with the results of the previous models.

Numerical Analysis of Non-Sphericity of Particles of Powder Material and their Effect on Packing Structure for Concentrated Solar Power Applications

Concentrated solar power technology represents a novel approach to generating solar power, characterized by high solar radiation density and uninterrupted operation. To store thermal energy and extend system operating hours, concentrated solar power systems rely on thermal energy storage. Selecting the appropriate heat storage media is crucial for designing a cost-efficient and straightforward system capable of withstanding high processing conditions. Among many powder materials, bauxite powder stands out as alternatives to conventional molten salt for thermal storage media due to its widespread availability and suitability for operating conditions. Particle shape, particularly its sphericity, plays a vital role in efficiently packing particles which directly affects to heat transfer properties. To investigate this effect, three samples of particles with different shape were formed: spheres, ellipses, and cylinders using Discrete Element Method with equal volume and material properties, and only variation of their sphericity factors. The study focused on the analysis of the porosity of the packed bed of particles with different shapes and its influence on the packing structure properties. The analysis of local packing fraction revealed that elliptical and cylindrical particles samples exhibited reduced heap formation and a more uniform distribution along z direction during vertical packing. The coordination number and radial distribution function analyses for these non-spherical particles showed a greater number of contacts between particles and disordered distribution of particles due to mechanical interlocking of non-spherical particles. The results indicated that the deviation of sphericity has a positive impact on the denser packing of particles generated more contacts and higher local packing fraction.

Gray correlation analysis between mechanical performance and pore characteristics of permeable concrete

Pervious concrete is a type of sustainable construction material used in pavement engineering due to its advantages in water permeability, skid resistance, and sound absorption performance. This study investigated the pore characteristics of five types of pervious concrete, along with the uniaxial compression behaviors and semi-circular bending (SCB) fracture performance. The critical stress intensity factor was obtained using the equivalent energy approach. Gray relational analysis was conducted between the mechanical performance and the pore characteristics. Pore indicators include the effective porosity, 2D porosity along the specimen height, 2D coordination number (CN) along the specimen height, 2D pore radius, 2D specific surface area, 3D CN, 3D pore radius, and 3D virtual porosity. Results indicated that from M1 to M3, when cement mortar content decreased from 520 kg/m3 to 408 kg/m3, the effective porosity increased from 18.95% to 27.43%, and the compressive strength and fracture toughness decreased obviously. For the same type of pervious concrete mixture, with the increase of notch length, the critical stress intensity factor decreased accordingly. The 2D porosity distribution, 2D pore size distribution, and 2D CN distribution presented a symmetrical trend of decreasing, then increasing, and then decreasing with the increase of specimen height. Aggregate with the gradation within the range of 9.5 mm to 19.0 mm generated a higher compressive strength, larger effective porosity and larger pore size compared to the gradation of 4.75 mm∼9.5 mm. The results of CN and pore radius generated from the 3D pore network model were larger than the corresponding 2D values. The gray correlation analysis revealed that 2D pore radius and 3D pore radius played critical effects on the uniaxial compressive strength and fracture resistance, while 2D CN and 3D CN were the least significant factors affecting the mechanical performance.

Selective laser sintering of amorphous nanoparticles: Molecular dynamics simulations

The paper investigates the process of liquid-phase sintering of amorphous iron-based nanoparticles by the method of molecular dynamics simulations. The classical molecular dynamics package LAMMPS was used for modeling. Visual analysis of the atomic configurations of nanoparticles during their rapid cooling revealed the self-purification effect of the particles. Partial pair correlation functions and coordination number distribution functions were used to analyze the atomic structure of nanoparticles after sintering. As a result of the analysis of the main structural parameters, which were obtained using the specified functions, differences in the atomic composition and structure of the volume and surface of nanoparticles were established.

Molecular dynamics simulation of microstructure and thermophysical properties of LiCl–CaCl2eutectic molten salt

AbstractChloride molten salt systems are widely used as electrolytes for molten salt electrolysis because of their relatively low eutectic temperatures and good thermal stability, but there is a serious lack of data on the thermophysical properties of chloride molten salts at elevated temperatures, whereas the nature of the electrolyte is very important for the electrolysis process. In this paper, the variation of the microstructure and thermophysical properties of the binary mixed LiCl–CaCl2molten salt system with temperature and composition is calculated using molecular dynamics (MD) simulations based on the BMH potential. The microscopic conformations observed in LiCl–CaCl2molten salts are mainly irregular, distorted tetrahedra and octahedra, which dynamically coexist, as analyzed by the radial distribution function, coordination number and angular distribution function. In addition, the effects of temperature and composition on the density, ionic self-diffusion coefficient, shear viscosity, and ionic conductivity of the molten salts were investigated, and the relationships of the thermophysical properties of LiCl–CaCl2molten salts with temperature and composition were obtained, which provide fundamental thermophysical data for the molten salt electrolytes.

Chemdex: quantification and distributions of valence numbers, oxidation numbers, coordination numbers, electron numbers, and covalent bond classes for the elements.

The paper introduces Chemdex, a freely accessible web-based database of over 70 000 compounds characterised by crystallography from across the periodic table. Its software calculates for an atom of interest within each compound classifications including valence number, oxidation number, coordination number, electron number, several covalent bond classifications, and the attached atom set. Users may explore distributions of these classifications by percentages and heat map displays for individual elements or sets of elements, or in several cases for one classification plotted against a second. These properties often display clear periodicity. Based upon distributions across the periodic table of valence numbers, electron numbers, coordination numbers, and attached atom data suggestions are made regarding the placement of hydrogen in the periodic table, membership of group 3 in the periodic table, locations of the early actinoids in the periodic table, and assignments of certain elements as metalloids.

Unveiling the structure of calcium borate glass by neutron diffraction, MAS-NMR, Raman and first-principles calculations

Calcium borate glass is a challenging system, due to the rich variety of borate units and alterable boron coordination number, and hence they are also of particular interest for academic study. In this study, neutron diffraction, solid-state magic angle spinning nuclear magnetic resonance (MAS-NMR) and Raman spectroscopy were combined with first-principles calculations, experience potential structure refinement (EPSR) simulations to investigate the structure of xCaO·(1-x)11B2O3 glass. Property data such as density and glass transition temperature were also measured to study the impact of structural changes on the properties of glasses. The structure of calcium borate glass was analyzed using various methods including structural unit, atomic pair bond length, coordination number, ring structure, chain structure and cavity distribution. Most notably, the atomic labeling method plays a key role in the structural analysis process. Through this method some new conclusions were obtained. For example, the increase of CaO content mainly leads to the conversion of BO to NBO in the glass structure, and has little effect on the content of BO3 and BØ4-. This result further leads to an increase in Ca-NBO interactions. Interestingly, the bond length distribution of Ca-NBO and B-NBO is shorter than that of Ca-BO and B-BO, so the glass structure becomes denser. The structural reason why the density of calcium borate glass is less than the crystal density was also revealed. This study contributes to unveiling the structure of calcium borate glass and exploring the correlation between the structure and properties of calcium borate glass, which provides a valuable reference for developing more valuable amorphous materials.

Insight of temperature and density-driven transition of sawtooth penta-silicene nanoribbons via molecular dynamics study

The density-driven transitions of sawtooth penta-silicene nanoribbon (PSiNR) at different temperature points are investigated using classical Molecular Dynamics (MD). Under the confined z-boundary, we found that the penta-to-tetra and the liquid-to-tetra silicene transitions present a clear dependence on the temperature. Insights of the thermodynamics and structural properties such as ring size (RS) and coordination number (CN) distributions, radial distribution function, buckling, and interatomic distance have been carefully analyzed. Above the melting temperature, the liquid 2D silicene tends to direct aggregate from scattered low-membered ring domains into the tetra-silicene nanoribbon (TSiNR). Meanwhile, the transition in sawtooth penta-SiNR takes place from the boundary into the inner region. A phase diagram indicating a dependence on density and temperature is presented. On the other hand, without the confinement of the z-boundary, the sawtooth penta-SiNRs express the periodic rippling. High temperature and axial compressive strain also could induce the PSiNR roll into the nanotube. Our simulations provided insights into the structural variations involving both the density and temperature of sawtooth PSiNRs.

Multi-scale pore network modelling to evaluate connectivity in ceramic composites

Complex morphologies, such as open or connected feature networks, are present in a wide variety of materials. Characteristics of these networks can impact key performance attributes of the materials themselves, affecting transport properties such as thermal conductivity. Therefore, it is critical to analyze the microstructure of these materials to gain a better understanding of the fundamental characteristics of the morphology. This study utilized pore network modeling as a method to extract morphological information on the solid network formed by boron nitride ceramic flakes in a polymeric resin matrix and uses the characteristics of the model to analyze the connectivity of the flakes. In this work, Micro-CT and FIB/SEM tomography were used in tandem to provide complimentary analyses of the microstructure and nanostructure, respectively, of the flake network to understand how this may contribute to transport properties of the material. Rather than a pore network model (PNM), the flake network model (FNM) was extracted from the tomographic datasets and the coordination number distribution was determined for the flakes detected in each. Micro-CT analysis showed that the flakes had formed a cage-like network around the exterior of the sample with limited connectivity in the interior, likely due to flake agglomeration at the outer surface of the material. A comparison of the full and interior-only Micro-CT FNMs indicated lower connectivity in the interior. This was confirmed by flow rate models generated from the network analysis for the flake contact points. The FNM extracted from the FIB/SEM tomography dataset exhibited similar connectivity compared to the interior-only FNM, indicating that the connectivity of the material was consistent when measured at the micron scale and at the nanometer scale.

Melting Process of the Two-Dimensional Material BN: Insights from Molecular Dynamics Simulations

The structure of the two-dimensional BN containing 9941 atoms has been studied by classical molecular dynamics simulation with Tersoff potential. The periodic boundary condition is applied to the two x and y directions, while the z direction is free. The analysis results of the function of total energy per atom and heat capacity, mean squared displacement, diffusion coefficient, radial distribution function, distribution of coordination number, angle distribution, and ring statistics show that the melting point of the material is about 4600 K. This value is higher than the experimental value as well as the previous simulation results. The observations also show that the melting process begins at the corners and edges and then spreads across the face of the model. The breakage of the B-N bond leads to the formation of clusters of N2 molecules and B with different sizes.

Evolution of Structure and Dynamic of Simulated CdTe Model According to the Melting Process

CdTe model containing 9955 atoms was built and simulated by the method of molecular dynamics (MD) with the Stillinger–Weber (SW) potential. The model was obtained by heating from 300K to 3000K at a rate of 1013 K/s. The structural properties of the model are investigated through the total energy per atom, the partial radial distribution functions (PRDFs), the coordination number distributions and the bond-angle distributions. Calculation results have shown that the model has a melting point of about 1370K and changes from tetrahedral to octahedral structure when the temperature of CdTe model is above the melting point. These results are consistent with the reported experimental and simulation results. In addition, the mobility of Cd and Te atoms is also researched through the mean squared displacement (MSD) and the diffusion coefficients of Cd and Te atoms.

Development of a 3D Computer Simulation Model Using C++ Methods

The article presents modelling of spherical elements based on the developed computer model. We recorded the main combinations of spherical particles during filling, which are formed in the hopper. It was found that the most likely combination that occurs when modelling spherical elements consists of three balls. It should be noted that in the cross-section of such a combination passing through the center of the balls, an equilateral triangle is formed. And in the cross-section of the structure, which consists of four spherical balls, a rhombus is formed, if you connect the centers of these spherical elements. It is worth noting that from this formed combination of spherical elements, it can be seen that the rhombus forms two smaller equilateral triangles that fix the process of pushing the spherical balls apart. In turn, the process of pushing spherical elements apart made it possible to fix the contact between spherical elements, as well as to state the stable position of each (individual) particle. This paper also presents the main fragments of encoding the source text of a 3D computer model for modelling spherical elements, which made it possible to optimize the model parameters. It was found that from the obtained data on the distribution of coordination numbers for different volume fillings of spherical elements, it follows that the largest filling was 72 %, which corresponds to the state when 112 lobules have an average coordination number of 3,92.

Gamma‐ray irradiation effect on microstructure and physical performances of porous silica

AbstractIn this paper, the gamma irradiation effect on the microstructure and physical performances of porous silica, including mechanical, thermal, and optical performances, are systematically investigated by using molecular dynamics and density‐functional theory‐based methods. The study of bond angle distribution, pair distribution function, coordination number distribution, and average ring size distribution show that, after gamma‐ray irradiation, the microstructure of porous silica is obviously modified. The tight packing of SiO2 tetrahedrons in the porous silica network is broken by gamma‐ray irradiation. Defects of three‐coordinated Si and non‐bridging oxygen are induced by gamma‐ray irradiation. Moreover, we find that the defects concentrations rapidly grow as gamma‐ray dose increases. The mechanical, thermal, and optical performances of porous silica are all seriously degenerated by gamma‐ray irradiation. Our results show that, for mechanical performance, Young's modulus, Bulk modulus, and Shear modulus first decrease and then keep stable as gamma‐ray dose increases, but the change of Poisson's ratio is slight. For thermal performance, the thermal conductivity decreases exponentially as gamma‐ray dose increases. For optical performance, light absorption coefficients increase exponentially and light transmittance drops as gamma‐ray dose increases in the working range (photon energy range around 3.5 eV) of inertial confinement fusion. Present work is expected to be valuable for studying the degradation mechanism of silicate materials under gamma radiation and developing gamma‐ray irradiation protection technology.