Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Al(HO)3”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Thermodynamic Understanding of Impurity Phase Segregation in a PdCrO 2 /CuCrO 2 Heterostructure

Abstract PdCrO 2 films are synthesized on CuCrO 2 buffer layers on Al 2 O 3 substrates. This synthesis is accompanied by impurity phase segregation, which hampers the synthesis of high quality PdCrO 2 films. Potential causes ofth impurity phase segregation are studied by using a combination of experiments and ab initio calculations. X‐ray diffraction and scanning transmission electron microscopy experiments reveal impurity phases of Cu x Pd 1 − x alloy and chromium oxides, Cr 2 O 3 and Cr 3 O 4 , in PdCrO 2 . Calculations determine that oxygen deficiency can cause the impurity phase segregation. Therefore, preventing oxygen release from delafossites can suppress the impurity phase segregation. The amounts of Cr 2 O 3 and Cr 3 O 4 depend on temperature and oxygen partial pressure. A reasonable theory‐based explanation for this experimental observation is provided.

2D metal↗
DOE OSTI2023

Effect of nitrogen doping and pressure on the stability of $\mathrm{LuH_3}$

A report on the near-ambient superconductivity in a nitrogen-doped lutetium hydride has stimulated great interest in this material [Dasenbrock-Gammon et al. , Nature (London) 615, 244 (2023)]. While its superconductivity is still a subject of debate, the structure of the claimed cubic phase remains uncertain. In this work, we study the effect of nitrogen doping and pressure on the energetic and dynamic stability of cubic $\mathrm{LuH_3}$. Our findings indicate that both pressure and nitrogen doping can enhance the stability of the cubic $\mathrm{LuH_3}$ phase. Here, we propose a $\mathrm{Lu_8}$$\mathrm{H_{21}}$$\mathrm{Lu}$ structure that exhibits a stable phonon, reasonable thermodynamic stability at 1 GPa, and an x-ray diffraction pattern similar to the experimental data. However, we do not observe electron-phonon coupling in the zone-center phonon modes of these crystal structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Reversing the Irreversible: Thermodynamic Stabilization of LiAlH 4 Nanoconfined Within a Nitrogen-Doped Carbon Host

A general problem when designing functional nanomaterials for energy storage is the lack of control over the stability and reactivity of metastable phases. Using the high-capacity hydrogen storage candidate LiAlH4 as an exemplar, we demonstrate an alternative approach to the thermodynamic stabilization of metastable metal hydrides by coordination to nitrogen binding sites within the nanopores of N-doped CMK-3 carbon (NCMK-3). The resulting LiAlH 4 @NCMK-3 material releases H 2 at temperatures as low as 126 °C with full decomposition below 240 °C, bypassing the usual Li 3 AlH 6 intermediate observed in bulk. Moreover, >80% of LiAlH 4 can be regenerated under 100 MPa H 2 , a feat previously thought to be impossible. Nitrogen sites are critical to these improvements, as no reversibility is observed with undoped CMK-3. Density functional theory predicts a drastically reduced Al–H bond dissociation energy and supports the observed change in the reaction pathway. Finally, the calculations also provide a rationale for the solid-state reversibility, which derives from the combined effects of nanoconfinement, Li adatom formation, and charge redistribution between the metal hydride and the host.

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

A Sobering Examination of the Feasibility of Aqueous Aluminum Batteries

Recent reports on aqueous aluminum (Al) batteries boast encouraging metrics including high energy density, nonflammability, and low raw material cost. However, the poorly understood solvation environment around the trivalent aluminum ions (Al3+) raises uncertainty regarding their feasibility. Furthermore, the corrosiveness of the acidic solutions and degree of hydration of aluminum trifluoromethanesulfonate (Al(OTf)3) salt is routinely mischaracterized. This study provides a rigorous examination of aqueous Al electrolytes, with the first compelling evidence for a dynamic octahedral solvation structure around Al3+ dominated by labile water and OH-, without Al-OTf contact ion pairs, at high concentrations. This solvation behavior and the concomitant electrostatic hydrolysis differs greatly from previously reported water-in-salt electrolytes, and occurs due to the high charge density of the Lewis acidic Al3+. Nuclear magnetic resonance spectroscopy and other physicochemical measurements quantitatively reveal how species activity evolves with concentration. The high proton activity observed in transport and electrochemical measurements constitutes a prohibitive barrier to the feasibility of rechargeable aluminum batteries with aqueous electrolyte.

Pastel, Glenn R.↗
DOE OSTI2020

Remote heteroepitaxy of GaN microrod heterostructures for deformable light-emitting diodes and wafer recycle

There have been rapidly increasing demands for flexible lighting apparatus, and micrometer-scale light-emitting diodes (LEDs) are regarded as one of the promising lighting sources for deformable device applications. Herein, we demonstrate a method of creating a deformable LED, based on remote heteroepitaxy of GaN microrod (MR) p - n junction arrays on c -Al 2 O 3 wafer across graphene. The use of graphene allows the transfer of MR LED arrays onto a copper plate, and spatially separate MR arrays offer ideal device geometry suitable for deformable LED in various shapes without serious device performance degradation. Moreover, remote heteroepitaxy also allows the wafer to be reused, allowing reproducible production of MR LEDs using a single substrate without noticeable device degradation. The remote heteroepitaxial relation is determined by high-resolution scanning transmission electron microscopy, and the density functional theory simulations clarify how the remote heteroepitaxy is made possible through graphene.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction of superconductivity in metallic boron–carbon compounds from 0 to 100 GPa by high-throughput screening

Boron–carbon compounds have been shown to have feasible superconductivity. In our earlier paper [Zheng et al., Phys. Rev. B, 2023, 107, 014508], we identified a new conventional superconductor of LiB 3 C at 100 GPa. Here, we aim to extend the investigation of possible superconductivity in this structural framework by replacing Li atoms with 27 different cations from periods 3, 4, and 5 under pressures ranging from 0 to 100 GPa. Using the high-throughput screening method of zone-center electron–phonon interaction, we found that ternary compounds like CaB 3 C, SrB 3 C, TiB 3 C, and VB 3 C are promising candidates for superconductivity. The consecutive calculations using the full Brillouin zone confirm that they have a T c of <31 K at moderate pressures. In conclusion, our study demonstrates that fast screening of superconductivity by calculating zone-center electron–phonon coupling strength is an effective strategy for high-throughput identification of new superconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Spatial variation of the physical conditions of molecular gas in galaxies

Multi-line studies of CO-12, CO-13, C-18O, HCN, and HCO(+) at 3 mm, 1.3 mm, and 0.8 mm using the Institute for Radio Astronomy in the Millimeter range (IRAM) 30 m telescope, with the IRAM superconductor insulator superconductor (SIS) receivers and the Max Planck Institute for External Physics (MPE) 350 GHz SIS receiver, show that the densities and temperatures of molecular gas in external galaxies change significantly with position. CO-12 measures the densities and temperature of diffuse interclump molecular gas, but not the bulk of the molecular gas. Simple one-component models, with or without external heating, cannot account for the weakness of the CO-12 J = 3 to 2 line relative to J = 2 to 1 and J = 1 to 0. CO-12 does not trace the bulk of the molecular gas, and optical depth effects obviate a straightforward interpretation of CO-12 data. Instead, researchers turned to the optically thin CO isotopes and other molecular species. Isotopic CO lines measure the bulk of the molecular gas, and HCN and HCO(+) pick out denser regions. Researchers find a warm ridge of gas in IC 342 (Eckart et al. 1989), denser gas in the starburst nucleus of IC 342, and a possible hot-spot in NGC 2903. In IC 342, NGC 2146, and NGC 6764, the CO-13 J = 2 to 1 line is subthermally populated, implying gas densities less than or equal to 10(exp 4) cm(-3).

Jackson, James M.↗
NASA NTRS2021

Constraining the Neutron Star Mass–Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes

We describe the X-ray pulse profile models we use and how we use them to analyze Neutron Star Interior Composition Explorer(NICER)observations of rotation-powered millisecond pulsars to obtain information about the mass–radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we analyzed in Miller et al. and Riley et al. and describe a cross-verification of computations of the pulse profiles of a star with R/M 3, in case stars this compact need to be considered in future analyses. We also present our early cross-verification efforts of the parameter estimation procedures used by Miller et al. and Riley et al. by analyzing two distinct synthetic data sets. Both codes yielded credible regions in the mass–radius plane that are statistically consistent with one another, and both gave posterior distributions for model parameter values consistent with the values that were used to generate the data. We also summarize the additional tests of the parameter estimation procedure of Miller et al. that used synthetic pulse profiles and the NICER pulse profile of PSR J0030+0451. We then illustrate how the precision of mass and radius estimates depends on the pulsar’s spin rate and the size of its hot spot by analyzing four different synthetic pulse profiles. Finally, we assess possible sources of systematic error in the estimates made using this technique, some of which may warrant further investigation.

Slavko Bogdanov↗
DOE OSTI2021

Strong Fermi-level pinning at metal contacts to halide perovskites

The performance of halide perovskite-based electronic and optoelectronic devices is often related to interfacial charge transport. To shed light on the underlying physical and chemical properties of CH 3 NH 3 PbI 3 (MAPbI 3 ) in direct contact with common electrodes Al, Ti, Cr, Ag, and Au, the evolution of interfacial properties and Fermi level pinning is systematically studied. Given a unique experimental facility, pristine interfaces without any exposure to ambient air were prepared. We observe aggregation of substantial amounts of metallic lead (Pb 0 ) at the metal/MAPbI 3 interface, resulting from the interfacial reaction between the deposited metal and iodine ions from MAPbI 3 . It is found that the Schottky barrier height at the metal/MAPbI 3 interface is independent of the metal work function due to strong Fermi level pinning, possibly due to the metallic Pb 0 aggregates, which act as interfacial trap sites. The charge neutrality level of MAPbI 3 is consistent with the energy level of Pb 0 -related defects, indicating that Pb 0 interfacial trap states can be nonradiative recombination sites. Here this work underlines that control of chemical bonding at interfaces is a key factor for designing future halide perovskite-based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Yb Substitution and Ultralow Thermal Conductivity of the Ca 3– x Yb x AlSb 3 (0 ≤ x ≤ 0.81(1)) System

Here, a series of Yb-substituted Zintl phases in the Ca 3–x Yb x AlSb 3 (0 ≤ x ≤ 0.81(1)) system has been synthesized by initial arc melting and post-heat treatment, and their isotypic crystal structures were characterized by both powder and single crystal X-ray diffraction analysis. All four title compounds adopted the Ca 3 AlAs 3 -type structure (space group Pnma, Pearson code oP28, Z = 4). The overall structure can be described as a combination of the 1-dimensional (1D) infinite chain of ∞ 1 [Al(Sb 2 Sb 2/2 )] formed by two vertices sharing [AlSb 4 ] tetrahedral moieties and three Ca 2+ /Yb 2+ mixed sites located in between these 1D chains. The charge balance and the resultant independency of the 1D chains in the title system were explained by the Zintl-Klemm formalism [Ca 2+ /Yb 2+ ] 3 [(4b-Al 1– )(1b-Sb 2– ) 2 (2b-Sb 1– ) 2/2 ]. A series of DFT calculations proved that (1) the band overlap between the d-orbital states from two types of cations and the p-orbital states from Sb at the high symmetry Γ point implied a heavily doped degenerate semiconducting behavior of the quaternary Ca 2 YbAlSb 3 model and (2) the site preference of Yb for the M1 site was due to the electronic-factor criterion based on the Q values of each atomic site. The electron localization function calculations also proved that the two different shapes of lone pairs of the Sb atoms—the “umbrella-shape” and the “C-shape”—are determined by local geometry and the coordination environment on the anionic frameworks. Thermoelectric measurements of the quaternary title compound Ca 2.19(1) Yb 0.81 AlSb 3 showed an approximately two times larger ZT value than that of ternary Ca 3 AlSb 3 at 623 K due to increased electrical conductivity and ultralow thermal conductivity originated from Yb substitution for Ca.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

DE-1 observations of polar O(+) stream bulk parameters and comparison with a model of the centrifugally-accelerated polar wind

A survey of bulk parameters of analyzable O(+) outward streams in the mid-altitude (3-4.7 R(sub E) geocentric distance) polar cap magnetosphere is obtained from measurements by the Retarding Ion Mass Spectrometer (RIMS) aboard the Dynamics Explorer-1 (DE-1) spacecraft. There is wide scatter in the obtained densities, but they do display discernible trends: the average O(+) density in these streams decreases from about 60 ions/cc at 3.5 R(sub E) to about 1 ion/cc at 4.6 R(sub E). The streaming velocities are somewhat more defined, and their average increases from about 8 km/s at 3.5 R(sub E) to about 12 km/s at 4.6 R(sub E). The densities and bulk velocities are inversely correlated. We have further compared these observational trends with model profiles for the centrifugally-accelerated polar wind as recently described by Horwitz et al. (1994). The large outflow velocities observed can be understood in part as centrifugally-driven by convection with ionospheric electric field magnitudes of the order 50-70 mV/m, perhaps including plasma expansion effects.

Ho, C. Wing↗
NASA NTRS1994

DE1 observations of polar O(+) stream bulk parameters and comparison with a model of the centrifugally-accelerated polar wind

A survey of bulk parameters of analyzable O(+) outward streams in the midaltitude (3-4.7 R(sub E) geocentric distance) polar cap magnetosphere is obtained from measurments by the Retarding Ion Mass Spectrometer (RIMS) aboard the Dynamics Explorer 1 (DE1) spacecraft. There is wide scatter in the obtained densities, but they do display discernible trends: the average O(+) density in these decreases from about 30 ions/cu cm at 3.5 R(sub E) to about 1 ion/cu cm at 4.6 R(sub E). The streaming velocities are somewhat more defined, and their average increases from about 8 km/s at 3.5 R(sub E) to about 12 km/s at 4 R(sub E). The outward ion flux increases by a factor of 4 with magnetic activity from K(sub p) = 2 to 5. We have further compared these observational trends with model profiles for the centrifugally-accelerated polar wind as recently described by Horwitz et al. (1994). The large outflow velocities observed can be understood in part as centrifugally-driven by convection with ionospheric electric field magnitudes of the order 50-70 mV/m, purhaps including plasma expansion effects.

Wing Ho, C.↗
DOE OSTI2026

Mechanism of ion pairing–induced co-ion penetration into the stern layer

In classical electric double layer (EDL) theory, co-ions—ions carrying the same charge as the electrode—are assumed to be excluded from the Stern layer because of electrostatic repulsion. However, ion pairing with adsorbed counterions may enable co-ions to access the EDL, a phenomenon that remains underexplored. Key questions therefore remain unresolved, including how deeply co-ions penetrate, how ion-specific properties govern them, and how ion pairing within the EDL differs from bulk solution. Here, in this study, we employ molecular dynamics simulations combined with well-tempered metadynamics to quantify co-ion penetration and ion-pairing thermodynamics at positively charged Mg-Al layered double hydroxide interfaces (σ = +20.6 μC/cm 2 ) in aqueous chloride electrolytes across 0.005–3 M range. We reveal a concentration-driven crossover: at dilute conditions, counterion–co-ion contact ion pairing (CIP) is thermodynamically unfavorable or only weakly favorable (ΔG CIP = +0.12 kcal/mol for Li + to −0.42 kcal/mol for K + ), leading to co-ion exclusion. As concentration increases (≥0.5 M), enhanced screening reduces surface repulsion and stabilizes CIP thermodynamics (ΔG CIP = −0.89 kcal/mol for Li + to −1.05 kcal/mol for K + at 3 M), enabling monovalent co-ions to penetrate the Stern layer following the hydration hierarchy K + > Na + > Li + , while divalent co-ions (Mg 2+ , Ca 2+ ) remain in the diffuse layer, forming only solvent-separated pairs. Comparison with bulk solution shows interfacial ion pairing is suppressed at low concentrations but converges to bulk-like thermodynamics at high concentrations, indicating surface effects vanish under strong screening. These results establish a quantitative framework linking ion-pairing thermodynamics to co-ion penetration and providing strategies to design electrochemical interfaces.

Co-ion penetration↗
DOE OSTI2026

Crystal Structures, Optical Behavior, and Magnetic Properties in Hydrated Lanthanide Iron Sulfates

Single crystals of LnFe(SO 4 ) 3 (H 2 O) 2 (Ln = La, Ce, Pr, Nd, Sm, Eu, Gd, Dy, Ho, Er, Tm; compounds 1–11) and LnFe(SO 4 ) 3 (H 2 O) (Ln = Tm, Yb, Lu; compounds 12–14) were synthesized under hydrothermal conditions. Single-crystal X-ray diffraction (SCXRD) analysis revealed that the dihydrated compounds (1–11) crystallize in centrosymmetric (CS) structures, with the lanthanide ions adopting eight-coordinate geometries. In contrast, the monohydrated compounds (12–14) exhibit noncentrosymmetric (NCS) structures, where the lanthanide ions are seven-coordinated. Vibrating sample magnetometry (VSM) confirmed that compounds 2, 4, and 7–11 are paramagnetic below 400 K, with compound 8 displaying the highest magnetic susceptibility. Compounds 1, 5, 6, 13, and 14 show a sharp increase in magnetic susceptibility at Néel temperatures ( T N ) of approximately 72, 76, 70, 58, and 56 K, respectively, indicating antiferromagnetic ordering. High-temperature magnetic susceptibility measurements further support the presence of antiferromagnetic transitions in these compounds. Second harmonic generation (SHG) measurements showed that the noncentrosymmetric Yb (compound 13) and Lu (compound 14) compounds exhibit SHG intensities of 0.3× and 1.6× that of potassium dihydrogen phosphate (KDP), respectively.

anions↗
NASA NTRS1994

Three-dimensional structure of Schistosoma japonicum glutathione S-transferase fused with a six-amino acid conserved neutralizing epitope of gp41 from HIV

The 3-dimensional crystal structure of glutathione S-transferase (GST) of Schistosoma japonicum (Sj) fused with a conserved neutralizing epitope on gp41 (glycoprotein, 41 kDa) of human immunodeficiency virus type 1 (HIV-1) (Muster T et al., 1993, J Virol 67:6642-6647) was determined at 2.5 A resolution. The structure of the 3-3 isozyme rat GST of the mu gene class (Ji X, Zhang P, Armstrong RN, Gilliland GL, 1992, Biochemistry 31:10169-10184) was used as a molecular replacement model. The structure consists of a 4-stranded beta-sheet and 3 alpha-helices in domain 1 and 5 alpha-helices in domain 2. The space group of the Sj GST crystal is P4(3)2(1)2, with unit cell dimensions of a = b = 94.7 A, and c = 58.1 A. The crystal has 1 GST monomer per asymmetric unit, and 2 monomers that form an active dimer are related by crystallographic 2-fold symmetry. In the binding site, the ordered structure of reduced glutathione is observed. The gp41 peptide (Glu-Leu-Asp-Lys-Trp-Ala) fused to the C-terminus of Sj GST forms a loop stabilized by symmetry-related GSTs. The Sj GST structure is compared with previously determined GST structures of mammalian gene classes mu, alpha, and pi. Conserved amino acid residues among the 4 GSTs that are important for hydrophobic and hydrophilic interactions for dimer association and glutathione binding are discussed.

Schistosoma japonicum/enzymology/genetics/immunolo↗
DOE OSTI2021

Controlling quantum many-body dynamics in driven Rydberg atom arrays

Dynamic stabilization of an array Large-scale systems comprising one-dimensional chains and two-dimensional arrays of excited atoms held in a programmable optical lattice are a powerful platform with which to simulate emergent phenomena. Bluvstein et al. built an array of up to 200 Rydberg atoms and subjected the system to periodic excitation. Under such driven excitation, they found that the array of atoms stabilized, freezing periodically into what looked like time crystals. Understanding and controlling the dynamic interactions in quantum many-body systems lies at the heart of contemporary condensed matter physics and the exotic phenomena that can occur. Science , this issue p. 1355

Science & Technology - Other Topics↗
NASA NTRS1999

Cepheid Calibration of the Peak Brightness of SNe Ia.: SN 1989B in NGC 3627 - 9

Repeated imaging observations have been made of NGC 3627 with the Hubble Space Telescope in 1997/98, over an interval of 58 days. Images were obtained on 12 epochs in the F555W band and on five epochs in the F8141,V band. The galaxy hosted the prototypical, "Branch normal", type la supernova SN 1989B. A total of 83 variables have been found, of which 68 are definite Cepheid variables with periods ranging from 75 days to 3.85 days. The de-reddened distance modulus is determined to be (m - M)(sub 0) = 30.22 +/- 0.12 (internal uncertainty) using a subset of the Cepheid data whose reddening and error parameters are secure. The photometric data of Wells et al. (1994), combined with the Cepheid data for NGC 3627 give MB(max) = -19.36 +/- 0.18 and M(sub V)(max) = -19.34 +/- 0.16 for SN 1989B. Combined with the previous six calibrations in this program, plus two additional calibrations determined by others gives the mean absolute magnitudes at maximum of (M(sub B)) = -19.48 +/- 0.07 for "Brunch normal" SNe Ia at this interim stage in the calibration program. Using the argument by Wells et al. (1994) that SN 1989B here is virtually identical in decay rate and colors at maximum with SN 198ON in NGC 1316 in the Fornax cluster, and that such identity means nearly identical absolute magnitude, it follows that the difference in the distance modulus of NGC 3627 and NGC 1316 is 1.62 +/- 0.03 mag. Thus the NGC 3627 modulus implies that (m - M)(sub 0) = 31.84 for NGC 1316. The second parameter correlations of M(max) of blue SNe la with decay rate, color at maximum, and Hubble type are re-investigated. The dependence of (M(max)) on decay rate is non-linear, showing a minimum for decay rates between 1.0 less than ADelta(sub m)15 less than 1.6. Magnitudes corrected for decay rate show no dependence on Hubble type, but a dependence on color remains. Correcting both the fiducial sample of 34 SNe la with decay-rate data and the current eight calibrating SNe la for the correlation with decay rate as well as color gives H(sub 0) = 60 +/- 2 (internal) km/s Mpc, in both B and V. The same value to within 4% is obtained if only the SNe la in spirals (without second parameter corrections) are considered. The correlation of SNe la color at maximum with M(max) cannot be due to internal absorption because the slope coefficients in B, V, and I with the change in magnitude are far from or even opposite to the canonical reddening values. The color effect must be intrinsic to the supernova physics. "Absorption" corrections of distant blue SNe la will lead to incorrect values of H(sub 0). The Cepheid distances used in this series are insensitive to metallicity differences. The zeropoint of the P-L relation is based on an assumed LMC modulus of (m - M)(sub 0) = 18.50. As this may have to be increased by 0(sup m).06 to 0(sup m).08, all distances in this paper will follow and Ho will decrease by 3 - 4%.

Saha, A.↗
DOE OSTI2020

Searching for correlations between vibrational spectral features and structural parameters of silicate glass network

Infrared (IR) and Raman spectroscopic features of silicate glasses are usually interpreted based on the analogy with those of smaller molecules, molecular clusters, or crystalline counterparts; this study tests the accuracy and validity of these widely cited peak assignment schemes by comparing vibrational spectral features with bond parameters of the glass network created by molecular dynamics (MD) simulations. A series of sodium silicate glasses with compositions of [Na 2 O] x [Al 2 O 3 ] 2 [SiO 2 ] 98- x with x = 7, 12, 17, and 22 were synthesized and analyzed with IR and Raman. A silica glass substrate and a crystalline quartz were also analyzed for comparison. Glass structures with the same compositions were generated with MD simulations using three types of potentials: fixed partial charge pairwise (Teter), partial diffuse charge potential (MGFF), and bond order-based charge transfer potential (ReaxFF). The comparison of simulated and experimental IR spectra showed that, among these three potentials tested, ReaxFF reproduces the concentration dependence of spectral features closest to the experimentally observed trend. Hence, the bond length and angle distributions as well as Si– Q n species and ring size distributions of silica and sodium silicate glasses were obtained from ReaxFF-MD simulations and further compared with the peak assignment or deconvolution schemes—which have been widely used since 1970s and 1980s—(a) correlation between the IR peak position in the Si–O stretch region (1050-1120 cm -1 ) and the Si–O–Si bond angle; (b) deconvolution of the Raman bands in the Si–O stretch region with the Q n speciation; and (c) assignment of the Raman bands in the 420-600 cm -1 region to the bending modes of (SiO) n rings with different sizes (typically, n = 3-6). The comparisons showed that none of these widely used methods is congruent with the bond parameters or structures of silicate glass networks produced via ReaxFF-MD simulations. This result invokes that the adequacy of these spectral interpretation methods must be questioned. Alternative interpretations are proposed, which are to be tested independently in future studies.

36 MATERIALS SCIENCE↗