Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pair distribution function analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Multitemporal Cross-Calibration of the Terra MODIS and Landsat 7 ETM+ Reflective Solar Bands

In recent years, there has been a significant increase in the use of remotely sensed data to address global issues. With the open data policy, the data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Enhanced Thematic Mapper Plus (ETM+) sensors have become a critical component of numerous applications. These two sensors have been operational for more than a decade, providing a rich archive of multispectral imagery for analysis of mutitemporal remote sensing data. This paper focuses on evaluating the radiometric calibration agreement between MODIS and ETM+ using the near-simultaneous and cloud-free image pairs over an African pseudo-invariant calibration site, Libya 4. To account for the combined uncertainties in the top-of-atmosphere (TOA) reflectance due to surface and atmospheric bidirectional reflectance distribution function (BRDF), a semiempirical BRDF model was adopted to normalize the TOA reflectance to the same illumination and viewing geometry. In addition, the spectra from the Earth Observing-1 (EO-1) Hyperion were used to compute spectral corrections between the corresponding MODIS and ETM+ spectral bands. As EO-1 Hyperion scenes were not available for all MODIS and ETM+ data pairs, MODerate resolution atmospheric TRANsmission (MODTRAN) 5.0 simulations were also used to adjust for differences due to the presence or lack of absorption features in some of the bands. A MODIS split-window algorithm provides the atmospheric water vapor column abundance during the overpasses for the MODTRAN simulations. Additionally, the column atmospheric water vapor content during the overpass was retrieved using the MODIS precipitable water vapor product. After performing these adjustments, the radiometric cross-calibration of the two sensors was consistent to within 7%. Some drifts in the response of the bands are evident, with MODIS band 3 being the largest of about 6% over 10 years, a change that will be corrected in Collection 6 MODIS processing.

MODIS↗
DOE OSTI2024

RMCProfile7 : reverse Monte Carlo for multiphase systems

This work introduces a completely rewritten version of the programRMCProfile(version 7), big-box, reverse Monte Carlo modelling software for analysis of total scattering data. The major new feature ofRMCProfile7is the ability to refine multiple phases simultaneously, which is relevant for many current research areas such as energy materials, catalysis and engineering. Other new features include improved support for molecular potentials and rigid-body refinements, as well as multiple different data sets. An empirical resolution correction and calculation of the pair distribution function as a back-Fourier transform are now also available.RMCProfile7is freely available for download at https://rmcprofile.ornl.gov/.

Chemistry↗
DOE OSTI2021

Powder Diffraction Crystallography : 98 Years as Plan B?

A historical essay written in 2014 to accompany a talk in the Transactions Symposium "100 Years of Crystallography" at the 2014 American Crystallographic Association Annual Meeting (talk at https://history.amercrystalassn.org/h-brian-toby-videos). Abstract: Powder diffraction patterns were used for structure solution almost contemporaneously with the Bragg’s initial single-crystal determinations. Structures have continued to be determined from powders, in the subsequent century. There has been a nearly explosive growth in this technique since the 1970’s. This was initially driven by neutron diffraction, which prompted the invention of Rietveld analysis, but synchrotrons have gone on to have at least as large an impact. There is no apparent upper limit to the scale of systems that can be structurally studied with powder diffraction; novel protein structures have been determined. When short-range ordering or disordered regions of a material are of interest, the Fourier transform of a powder diffraction pattern, the pair distribution function can be modeled.

Toby, Brian H.↗
DOE OSTI2023

Impedance-Based Detection of NO 2 Using Ni-MOF-74: Influence of Competitive Gas Adsorption

Chemically robust, low-power sensors are needed for the direct electrical detection of toxic gases. Metal–organic frameworks (MOFs) offer exceptional chemical and structural tunability to meet this challenge, though further understanding is needed regarding how coadsorbed gases influence or interfere with the electrical response. To probe the influence of competitive gases on trace NO 2 detection in a simulated flue gas stream, a combined structure–property study integrating synchrotron powder diffraction and pair distribution function analyses was undertaken, to elucidate how structural changes associated with gas binding inside Ni-MOF-74 pores correlate with the electrical response from Ni-MOF-74-based sensors. Data were evaluated for 16 gas combinations of N 2 , NO 2 , SO 2 , CO 2 , and H 2 O at 50 °C. Fourier difference maps from a rigid-body Rietveld analysis showed that additional electron density localized around the Ni-MOF-74 lattice correlated with large decreases in Ni-MOF-74 film resistance of up to a factor of 6 × 10 3 , observed only when NO 2 was present. These changes in resistance were significantly amplified by the presence of competing gases, except for CO 2 . Without NO 2 , H 2 O rapidly (<120 s) produced small (1–3×) decreases in resistance, though this effect could be differentiated from the slower adsorption of NO 2 by the evaluation of the MOF’s capacitance. Furthermore, samples exposed to H 2 O displayed a significant shift in lattice parameters toward a larger lattice and more diffuse charge density in the MOF pore. Evaluating the Ni-MOF-74 impedance in real time, NO 2 adsorption was associated with two electrically distinct processes, the faster of which was inhibited by competitive adsorption of CO 2 . Together, this work points to the unique interaction of NO 2 and other specific gases (e.g., H 2 O, SO 2 ) with the MOF’s surface, leading to orders of magnitude decrease in MOF resistance and enhanced NO2 detection. Finally, understanding and leveraging these coadsorbed gases will further improve the gas detection properties of MOF materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Search for the critical point of strongly-interacting matter in $^{40}$Ar + $^{45}$Sc collisions at 150A Ge V /c using scaled factorial moments of protons

The critical point of dense, strongly interacting matter is searched for at the CERN SPS in $^{40}$Ar + $^{45}$Sc collisions at 150A Ge V /c. The dependence of second-order scaled factorial moments of proton multiplicity distribution on the number of subdivisions of transverse momentum space is measured. The intermittency analysis is performed using both transverse momentum and cumulative transverse momentum. For the first time, statistically independent data sets are used for each subdivision number. The obtained results do not indicate any statistically significant intermittency pattern. An upper limit on the fraction of correlated proton pairs and the power of the correlation function is obtained based on a comparison with the Power-law Model developed for this purpose.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Models to Incorporate Reaction Mechanisms into DG-OSPREY: Fixed-Bed Simulations for Organic Iodides Capture Using Ag 0 Z

As one of the most potent radioisotopes released during spent nuclear fuel reprocessing, 129 I is strictly regulated and must be removed before discharge. Organic iodides (primarily alkyl iodides with different chain lengths, i.e., CH 3 I, C 4 H 9 I, and C 12 H 25 I) comprise ~2% of the total iodine in the reprocessing off-gases and are primarily present in vessel off-gas (VOG). Reduced silver mordenite (Ag 0 Z) is predominantly considered for the removal of radioiodine; however, its capture performance and underlying interaction processes with long-chain organic iodides are not fully understood. Two major tasks were accomplished in this study. First, to improve upon the previous experimental studies where Ag 0 Z was used to capture CH 3 I, C 4 H 9 I, and C 12 H 25 I at different concentrations, we comprehensively investigated the corresponding capture mechanisms by characterizing fully loaded Ag 0 Z samples. Second, computational codes were implemented to perform fixed-bed simulations that account for transport and reaction mechanisms. Scanning electron microscopy with energy dispersive X-ray analysis (SEM-EDX), powder X-ray diffraction (PXRD), UV-visible diffuse reflectance spectroscopy (UV-vis DRS), and thermogravimetric analysis (TGA) were conducted on Ag 0 Z samples that are saturated with I 2 , CH 3 I, C 4 H 9 I, and C 12 H 25 I), respectively. Results indicate that AgI is the predominant adsorption product regardless of the adsorbed iodine species, yet alkyl iodides with different carbon chain lengths may have different compositions of α- and γ-AgI. Synchrotron pair distribution function (PDF) measurements and TGA coupled with a Fourier transformed infrared detector (TGA FTIR) have been performed, and experimental data are currently analyzed. Results are expected to provide further insights into the adsorption mechanisms. The fixed-bed capture performance for CH 3 I was successfully simulated by solving mathematical equations that describe the underlying transport processes and adsorption reactions. The computational framework, Catalytic After Treatment System (CATS), that was originated in our research group, was used to solve the governing equations. Kinetic parameters, including the pore diffusivity and reaction rate constant were obtained by optimization techniques using data from thin-bed experiments performed at Oak Ridge National Laboratory. The performance of a deep bed predicted using the optimized parameters showed promising agreement with the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1988

Modeling and measurement of error propagation in a multimodule computing system

An error propagation model has been developed for multimodule computing systems in which the main parameters are the distribution functions of error propagation times. A digraph model is used to represent a multimodule computing system, and error propagation in the system is modeled by general distributions of error propagation times between all pairs of modules. Two algorithms are developed to compute systematically and efficiently the distributions of error propagation times. Experiments are also conducted to measure the distributions of error propagation times with the fault-tolerant microprocessor (FTMP). Statistical analysis of experimental data shows that the error propagation times in FTMP do not follow a well-known distribution, thus justifying the use of general distributions in the present model.

Shin, Kang G.↗
DOE OSTI2021

Measurement of the inclusive and differential Higgs boson production cross sections in the leptonic WW decay mode at $\sqrt{s} =$ 13 TeV

Measurements of the fiducial inclusive and differential production cross sections of the Higgs boson in proton-proton collisions at $ \sqrt{s} $ = 13 TeV are performed using events where the Higgs boson decays into a pair of W bosons that subsequently decay into a final state with an electron, a muon, and a pair of neutrinos. The analysis is based on data collected with the CMS detector at the LHC during 2016–2018, corresponding to an integrated luminosity of 137 fb$^{−1}$. Production cross sections are measured as a function of the transverse momentum of the Higgs boson and the associated jet multiplicity. The Higgs boson signal is extracted and simultaneously unfolded to correct for selection efficiency and resolution effects using maximum-likelihood fits to the observed distributions in data. The integrated fiducial cross section is measured to be 86.5 ± 9.5 fb, consistent with the Standard Model expectation of 82.5 ± 4.2 fb. No significant deviation from the Standard Model expectations is observed in the differential measurements.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quantifying a partial polyamorphic transition in a cerium-based metallic glass during cooling

Cerium-based metallic glasses are prototype polyamorphous systems with pressure-induced polyamorphic transitions extensively reported. Cooling typically has a similar effect on materials as compression with regard to reducing volume. However, previous studies show dramatically different behavior of Ce-based metallic glasses between cooling and compression, whose origin remains unclear. Here, using in situ low-temperature synchrotron high-energy x-ray diffraction, the structural evolution of a Ce68Al10Cu20Co2 metallic glass is accurately determined and analyzed by a structure factor and a reduced pair distribution function (PDF) during cooling from 298 to 83 K. An unusually large linear thermal expansion coefficient is revealed, which is associated with both continuous but inconsistent structural changes between the two subpeaks of the first atomic shell in terms of average bond lengths and coordination numbers. These phenomena are suggested to be attributed to a gradual 4f electron delocalization of only a minimal amount (∼2.6% at 83 K) of Ce atoms by quantitative analysis of the PDF data. However, a previously expected global polymorphic transition from a low-density amorphous state to a high-density amorphous state with an abrupt volume collapse is not observed. Moreover, electrical resistivity also shows a continuous increase during cooling without any sharp change. It is clarified that cryogenic temperatures could facilitate but are not powerful enough alone to trigger a global polymorphic transition in the Ce68Al10Cu20Co2 metallic glass, suggesting a wide distribution of its local atomic environment.

Chen, Zhi↗
DOE OSTI2025

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2025

Conditional Pseudo-Reversible Normalizing Flow for Surrogate Modeling in Quantifying Uncertainty Propagation

We introduce a conditional pseudo-reversible normalizing flow (PR-NF) that directly learns conditional probability distributions from noisy physical models to efficiently quantify both forward and inverse uncertainty propagation. Traditional surrogate modeling approaches approximate only the deterministic component of physical models, requiring separate noise characterization and computationally expensive sampling methods for inverse problems. Here, in this work, we develop the conditional PR-NF model to directly learn and efficiently generate samples from the conditional probability density functions (PDFs). The training process utilizes dataset consisting of input-output pairs without requiring prior knowledge about the noise and the function. Once trained, our model efficiently generates samples from conditional PDFs for any input within the training domain. Moreover, the pseudo-reversibility feature allows for the use of fully connected neural network architectures, which simplifies the implementation and enables theoretical analysis. We provide a rigorous convergence analysis of the conditional PR-NF model, showing its ability to converge to the target conditional PDF using the Kullback−Leibler divergence. To demonstrate the effectiveness of our method, we apply it to several benchmark tests and a real-world geologic carbon storage problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Measurement of the W-boson angular coefficients and transverse momentum in pp collisions at s=13 TeV with the ATLAS detector

The angular distributions of Drell–Yan lepton pairs provide sensitive probes of the underlying dynamics of quantum chromodynamics (QCD) effects in vector-boson production. This paper presents for the first time the measurement of the full set of angular coefficients together with the differential cross-section as a function of the transverse momentum of the W boson, in the full phase space of the decay leptons. The measurements are performed separately for the W-$$W^-$$ and W+$$W^+$$ channels. The analysis uses proton–proton collision data recorded by the ATLAS experiment at the Large Hadron Collider in 2017 and 2018, during special low-luminosity runs with a reduced number of interactions per bunch crossings (pile-up). The data correspond to an integrated luminosity of 338 pb-1$$^{-1}$$ at a centre-of-mass energy of s=13$$\sqrt{s} = 13$$ TeV. The low pile-up environment provides excellent experimental conditions for high-precision measurements of W-boson production. All results agree with theoretical predictions incorporating finite-order QCD corrections up to order αS2$$\alpha _S^2$$.

Aad, G↗
DOE OSTI2021

Characterization of zirconium carbide microspheres synthesized via internal gelation

Microspheres of zirconium carbide with grain size 40-320 nm were synthesized using internal gelation techniques and characterized by scanning electron microscopy, synchrotron X-ray diffraction, and neutron total scattering. Compared with a polycrystalline benchmark ZrC sample prepared by plasma vapor phase deposition, the microspheres display lower variance in grain size and considerably lower microstrain representative of a more homogenous internal microstructure. However, excess carbon was present in the microspheres, evident in both the X-ray and neutron diffraction data as well as a corresponding hypostoichiometric ZrC phase. Furthermore, the excess carbon phase is assumed to be pockets of carbon that remained unreacted through the combined internal gelation and subsequent carbothermic reduction synthesis process. Atomic-scale structural characterization with neutron PDF analysis confirmed the presence of localized nano-sized domains of graphite-like carbon material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Capture and Decomposition of the Nerve Agent Simulant, DMCP, Using the Zeolitic Imidazolate Framework (ZIF-8)

Understanding mechanisms of decontamination of chemical warfare agents (CWA) is an area of intense research aimed at developing new filtration materials to protect soldiers and civilians in case of state-sponsored or terrorist attack. In this study, we employed complementary structural, chemical, and dynamic probes and in situ data collection, to elucidate the complex chemistry, capture, and decomposition of the CWA simulant, dimethyl chlorophosphonate (DMCP). Our work reveals key details of the reactive adsorption of DMCP and demonstrates the versatility of zeolitic imidazolate framework (ZIF-8) as a plausible material for CWA capture and decomposition. Here, the in situ synchrotron-based powder X-ray diffraction (PXRD) and pair distribution function (PDF) studies, combined with Fourier transform infrared (FT-IR) spectroscopy, X-ray photoelectron spectroscopy (XPS), zinc K-edge X-ray absorption near edge structure (XANES), and Raman spectroscopies, showed that the unique structure, chemical state, and topology of ZIF-8 enable accessibility, adsorption, and hydrolysis of DMCP into the pores and revealed the importance of linker chemistry and Zn 2+ sites for nerve agent decomposition. DMCP decontamination and decomposition product(s) formation were observed by thermogravimetric analysis, FT-IR spectroscopy, and phosphorus (P) K-edge XANES studies. Differential PDF analysis indicated that the average structure of ZIF-8 (at the 30 Å scale) remains unchanged after DMCP dosing and provided information on the dynamics of interactions of DMCP with the ZIF-8 framework. Using in situ PXRD and diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), we showed that nearly 90% regeneration of the ZIF-8 structure and complete liberation of DMCP and decomposition products occur upon heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

From phase decomposition to evaporation: A multi-modal evaluation of thermally degraded model lightweight high-entropy alloy

Lightweight high-entropy alloys (LHEAs) have the potential to replace conventional lightweight materials due to their superior mechanical properties and thermal stability. However, the thermal degradation pattern of LHEAs from phase decomposition to evaporation is not clear. We develop a new Al-based dual phase (FCC + HCP) LHEA—AlTi 0.45 CuZn, and further investigate its thermal degradation behavior for potential high-temperature structural applications. Using multimodal advanced characterization techniques such as differential scanning calorimetry/thermogravimetric analysis, scanning/transmission electron microscopy, and synchrotron X-ray diffraction/pair distribution function (XRD/PDF), a sequence of thermal degradation events beyond the thermal phase stability limit—between 250 and 360 °C—is observed. These include phase decomposition at ~360 °C, Zn evaporation at ~750 °C, and LHEA melting at 880 °C which results in ~25% cumulative weight loss. The formation of Al-Ti phase off the AlTi 0.45 CuZn matrix is due to the largest negative mixing enthalpy for Al-Ti than other binary pairs. Similarly, Zn evaporation from AlTi 0.45 CuZn LHEA is due to its faster evaporation rate than other constituent elements. The high-resolution synchrotron XRD and PDF results support the aforementioned observations; in addition, they reveal local atomic arrangements, local strain, and sluggish grain growth in the LHEA. Among other LHEAs of close density range (5.55 ≤ ρ ≤ 5.85 g/cc), the investigated LHEA exhibits outstanding nano-indentation hardness values due to the coupled grain size effect and HCP phase strengthening of the FCC matrix. As the search for LHEAs for lightweight applications grows, this study shows the potential use of AlTi 0.45 CuZn LHEA for structural applications even at elevated temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

A deconvolution technique for Hubble Space Telescope FGS fringe analysis

A technique has been developed for directly transforming interferometer fringe visibility functions ('S curves') from the Hubble Space Telescope (HST) fine guidance sensors (FGSs) into intensity profiles of the program object. In a process analogous to Fourier transform image deconvolution, an S curve from a double star yields a pair of narrow profiles containing information on the separation (in one coordinate), and relative brightness of the two stars. The procedure has yielded high internal precision for the separation and relative brightness in tests with HST data. Simulations indicate that it can also deconvolve S curves from multiple stars or continuous intensity distributions such as resolvable stellar or extragalactic objects. Similar deconvolution analyses may be useful in other types of interferometry.

Hershey, John L.↗
NASA NTRS2003

Examination of the Atomic Pair Distribution Function (PDF) of SiC Nanocrystals by In-situ High Pressure Diffraction

Key properties of nanocrystals are determined by their real atomic structure, therefore a reasonable understanding and meaningful interpretation of their properties requires a realistic model of the structure. In this paper we present an evidence of a complex response of the lattice distances to external pressure indicating a presence of a complex structure of Sic nanopowders. The experiments were performed on nanocrystalline Sic subjected to hydrostatic or isostatic pressure using synchrotron and neutron powder diffraction. Elastic properties of the samples were examined based on X-ray diffraction data using a Diamond Anvil Cell (DAC) in HASYLAB at DESY. The dependence'of the lattice parameters and of the Bragg reflections width with pressure exhibits a ha1 nature of the properties (compressibilities) of the powders and indicates a complex structure of the grains. We interpreted tws behaviour as originating from different elastic properties of the grain interior and surface. Analysis of the dependence of individual interatomic distances on pressure was based on in-situ neutron diffraction measurements done with HbD diffractometer at LANSCE in Los Alamos National Laboratory with the Paris-Edinburgh cell under pressures up to 8 GPa (Qmax = 26/A). Interatomic distances were obtained by PDF analysis using the PDFgetN program. We have found that the interatomic distances undergo a complex, non-monotonic changes. Even under substantial pressures a considerable relaxation of the lattice may take place: some interatomic distances increase with an increase in pressure. We relate this phenomenon to: (1), changes of the microstructure of the densified material, in particular breaking of its fractal chain structure and, (2), its complex structure resembling that of a material composed of two phases, each with its distinct elastic properties.

Grzanka, E.↗
DOE OSTI2024

Insights on the structure and properties of sodium iron phosphate glasses from molecular dynamics simulations

Iron phosphate glasses are promising nuclear waste forms while more detailed understanding of their structures and structure-property relations are still needed to better design waste glass compositions. In this work we report studies of three series of sodium iron phosphate (NFP) glasses: 60P 2 O 5 -(40-x)Fe 2 O 3 -xNa 2 O (x = 0→40), (100–2x)P 2 O 5 -xFe 2 O 3 -xNa 2 O (x = 5→17.5) and one with different iron redox ratio, to understand the composition as well as the iron redox effects on the structure and properties of these glasses using molecular dynamics simulations with effective two-body and three-body potentials. Structural analyses, including pair distribution function, bond angle distribution, Q n distribution, and polyhedral connectivity, were performed to obtain in-depth information on short-range and medium-range structural features. The P-O pair distributions showed a first peak splitting with phosphorus-bridging and non-bridging oxygen contributions. This and the average P-O and other cation-oxygen bond distances are in excellent agreement with experiments. The coordination number of P 5+ remained four while that of Fe 3+ increased from 4.30 to 4.72 with decreasing Fe/Na ratio. Polyhedral linkage analysis showed [PO 4 ] units linked with [PO 4 ] and [FeO x ] through corner-sharing while the [PO 4 ]-[FeO x ] linkages become dominant for compositions with Fe 2 O 3 larger than 15 mol%. The effect of iron redox ratio on the structure of NFP glasses was also studied and it was found that bond lengths and coordination numbers were not strongly affected, while the reduction of iron introduced higher network distortions, as evident by O-P-O bond angle and Q n distribution. The glass transition temperature (T g ) showed a monotonic increase with Fe 2 O 3 in the first series, in good agreement with experiments, while those of the second series showed a maximum at P 2 O 5 = 82 mol%. Here, calculated elastic moduli were found to increase with Fe 2 O 3 in the first glass series, which was be explained by the increase of network connectivity, while those of the second series decrease with Fe 2 O 3 due to decrease of P 2 O 5 .

36 MATERIALS SCIENCE↗