Search NASA⌕ Search

SEARCH · Search NASA

Results for “Southern California”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Atacama Cosmology Telescope DR6 and DESI: Structure growth measurements from the cross-correlation of DESI legacy imaging galaxies and CMB lensing from ACT DR6 and 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 PR4

We measure the growth of cosmic density fluctuations on large scales and across the redshift range 0.3 < 𝑧 < 0.8 through galaxy clustering and the cross-correlation of the ACT data release 6 cosmic microwave background (CMB) lensing map and galaxies from the Dark Energy Spectroscopic Instrument Legacy Survey, using three galaxy samples spanning the redshifts of 0.3 ≲ 𝑧 ≲ 0.45, 0.45 ≲ 𝑧 ≲ 0.6, 0.6 ≲ 𝑧 ≲ 0.8. We adopt a scale cut where nonlinear effects are negligible, so that the cosmological constraints are derived from the linear regime. We determine the amplitude of matter fluctuations over all three redshift bins using Atacama Cosmology Telescope (ACT) data alone to be 𝑆 8 ≡ 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 =0.772 ± 0.040 in a joint analysis combining the three redshift bins and ACT lensing alone. Using a combination of ACT and Planck data we obtain 𝑆 8 = 0.765 ± 0.032. The lowest redshift bin used is the least constraining and exhibits a ∼2⁢𝜎 tension with the other redshift bins; thus we also report constraints excluding the first redshift bin, giving 𝑆 8 = 0.785 ± 0.033 for the combination of ACT and Planck. This result is in excellent agreement at the 0.3⁢𝜎 level with measurements from galaxy lensing, but is 1.8⁢𝜎 lower than predictions based on Planck primary CMB data. Understanding whether this hint of discrepancy in the growth of structure at low redshifts arises from a fluctuation, from systematics in data, or from new physics is a high priority for forthcoming CMB lensing and galaxy cross-correlation analyses.

cosmic microwave background↗
DOE OSTI2025

Neuromorphic ionic computing in droplet interface synapses

Ionic devices with memory capabilities can emulate neural functionality, enabling neuromorphic computing and biomedical applications. In this study, we report an ionic spiking synapse based on aqueous droplet interface bilayer assembly. Under stepwise triangular voltages, the device displays coupled memcapacitive-memristive behavior, showing noncrossing pinched hysteretic I-V loops. This hysteretic ion dynamics can be regulated by modifying bilayer components, reconstituting protein channels, or adjusting droplet assembly configuration. Droplet interface synapses (DIS) exhibit fundamental neuromorphic behaviors such as paired-pulse facilitation/depression, spike rate–dependent plasticity, Hebbian learning, and short-term associative learning under classical conditioning. We also used reservoir computing with DIS to implement two learning algorithms: a classification algorithm that recognizes handwritten digits and a reinforcement learning algorithm that learns to play a board game of tic-tac-toe.

Li, Zhongwu [Lawrence Livermore National Laborator↗
DOE OSTI2025

A terminology for scientific workflow systems

The term “scientific workflow” has evolved over the last two decades to encompass a broad range of compositions of interdependent compute tasks and data movements. It has also become an umbrella term for processing in modern scientific applications. Today, many scientific applications can be considered as workflows made of multiple dependent steps, and hundreds of workflow systems have been developed to manage and run these scientific workflows. However, no turnkey solution has emerged from the field to address the diversity of scientific processes and the infrastructure on which they are supposed to be implemented. Instead, new research problems requiring the execution of scientific workflows with some novel feature often lead to the development of an entirely new workflow system. A direct consequence of this situation is that many existing workflow management systems (WMSs) share some salient features, offer similar functionalities, and can manage the same categories of workflows but at the same time also have some distinct capabilities that can be important for specific applications. This situation makes researchers who develop workflows face the complex question of selecting a WMS. This selection can be driven by technical considerations, to find the system that is the most appropriate for their application and for the computing and storage resources available to them, or other factors such as reputation, adoption, strong community support, or long-term sustainability. To address this problem, a group of WMS developers and practitioners joined their efforts to produce a community-based terminology of WMSs. This paper summarizes their findings and introduces this new terminology to characterize WMSs. Furthermore, this terminology is composed of fives axes: workflow structure and characteristics, composition, orchestration, data management, and metadata capture. Each axis comprises several concepts that capture the prominent features of WMSs. Based on this terminology, this paper also presents a classification of 23 existing WMSs according to the proposed axes and terms.

Community-based terminology↗
DOE OSTI2022

Lead isotopic compositions of airborne particulates collected from snow and from dry deposition in the East River Watershed, Colorado and at the Southern Sierra Critical Zone Observatory, California

The data set consists of lead (Pb) isotopic ratios (206Pb/207Pb and 208Pb/207Pb, as well as 206Pb/204Pb, 207Pb/204Pb, 208Pb/204Pb) for particulate samples collected from snow and dry deposition. This was part of a study (O’Day et al. 2020) that was looking at the speciation of phosphorous in splits of these samples to evaluate the flux and bio-availability of phosphorous to the mountain sites in the study. The Pb isotopic data allow sourcing of long-range transported Pb, specifically the fraction of Pb sourced in Asia, and so a comparison of the relative distances of PM sources between the high and low elevation sites. Samples were collected at the East River Watershed, CO (dry deposition collected over Sept.-Oct. 2016, show sampling winter 2017) and at the Southern Sierra Critical Zone Observatory (dry deposition collected over Sept-Nov 2016, July-Sept 2017, June-July 2017). The data includes Pb concentration data, as well as calculated fraction of the Pb sourced in Asia. This data set was published in part in O’Day et al. (2020) Phosphorus Speciation in Atmospherically Deposited Particulate Matter and Implications for Terrestrial Ecosystem Productivity. ES&T, https://dx.doi.org/10.1021/acs.est.9b06150

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Programmable simulations of molecules and materials with reconfigurable quantum processors

Simulations of quantum chemistry and quantum materials are believed to be among the most important applications of quantum information processors. However, realizing practical quantum advantage for such problems is challenging because of the prohibitive computational cost of programming typical problems into quantum hardware. Here we introduce a simulation framework for strongly correlated quantum systems represented by model spin Hamiltonians that uses reconfigurable qubit architectures to simulate real-time dynamics in a programmable way. Our approach also introduces an algorithm for extracting chemically relevant spectral properties via classical co-processing of quantum measurement results. We develop a digital–analogue simulation toolbox for efficient Hamiltonian time evolution using digital Floquet engineering and hardware-optimized multi-qubit operations to accurately realize complex spin–spin interactions. As an example, we propose an implementation based on Rydberg atom arrays. In addition, we show how detailed spectral information can be extracted from the dynamics through snapshot measurements and single-ancilla control, enabling the evaluation of excitation energies and finite-temperature susceptibilities from a single dataset. To illustrate the approach, we show how to use the method to compute key properties of a polynuclear transition-metal catalyst and two-dimensional magnetic materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Quantum-State-Controlled Collisions of Ultracold Polyatomic Molecules

Collisions between ultracold calcium monohydroxide (CaOH) molecules are realized and studied. Inelastic collision rate constants are measured for CaOH prepared in ground and excited vibrational states, and the electric field dependence of these rates is measured for molecules in single quantum states of the parity-doubled bending mode. Theoretical calculations of collision rate coefficients are performed and found to agree with measured values. The lowest collisional loss rates are for states with repulsive long-range potentials that shield ultracold molecules from loss channels at short distance. These results unveil the collisional behavior of parity-doublet molecules in the ultracold regime and lay the foundation for future experiments to evaporatively cool polyatomic molecules to quantum degeneracy.

Vilas, Nathaniel B. [Harvard University, Cambridge↗
DOE OSTI2025

Materials laboratories of the future for alloys, amorphous, and composite materials

In alignment with the Materials Genome Initiative and as the product of a workshop sponsored by the US National Science Foundation, we define a vision for materials laboratories of the future in alloys, amorphous materials, and composite materials; chart a roadmap for realizing this vision; identify technical bottlenecks and barriers to access; and propose pathways to equitable and democratic access to integrated toolsets in a manner that addresses urgent societal needs, accelerates technological innovation, and enhances manufacturing competitiveness. Spanning three important materials classes, this article summarizes the areas of alignment and unifying themes, distinctive needs of different materials research communities, key science drivers that cannot be accomplished within the capabilities of current materials laboratories, and open questions that need further community input. Here, we provide a broader context for the workshop, synopsize the salient findings, outline a shared vision for democratizing access and accelerating materials discovery, highlight some case studies across the three different materials classes, and identify significant issues that need further discussion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Test of the Gravitational Force Law on Cosmological Scales Using the Kinematic Sunyaev-Zeldovich Effect

The mean pairwise velocity of massive halos reflects the gravitational force law on cosmic scales. For this work, we combine cosmic microwave background intensity maps from the Atacama Cosmology Telescope and a galaxy catalog from the Sloan Digital Sky Survey to estimate the mean pairwise velocity using the kinematic Sunyaev-Zeldovich (kSZ) effect. On scales from 30 to 230 megaparsecs, we constrain the gravitational acceleration between pairs of halos at separation 𝑟 to be 𝑔 ∝ 1/𝑟 𝑛 with 𝑛 = 2.1 ± 0.3, which is consistent with Newtonian gravity in an expanding spacetime (i.e., the standard Λ⁢ CDM model). This constraint shows agreement with an inverse quadratic radial dependence over the large distances separating galaxy halos, as expected in standard cosmology. Upcoming surveys have the potential to rule out 𝑛 = 1 at 10⁢𝜎 significance. Our results establish the kSZ effect as a powerful tool for testing gravity on cosmological scales.

alternative gravity theories↗
DOE OSTI2026

Developing Scenario‐Based Strategies for Health, Climate, and Environmental Preparedness: The One Health, One Earth Approach

Climate change amplifies many threats to human health. Despite advances in understanding climate change dynamics and impacts, there remains a critical gap in translating scientific knowledge into equitable, and community-driven health interventions. The inaugural One Earth, One Health workshop sought to explore this gap through human-centered design exercises involving interdisciplinary researchers from climate and Earth sciences, engineering, epidemiology, microbiology, and environmental health. Although participants did not co-develop solutions with affected communities, they used stakeholder role-playing to guide ideation and lay groundwork for actionable plans. Through these methods, participants identified community needs and proposed prototype solutions to alleviate health threats exacerbated by global environmental change. Prototypes were organized around infectious diseases, extreme weather, and air quality, as illustrative themes rather than an exhaustive set of risks. Key solutions included strategies for anticipatory systems and early warning (e.g., integrating environmental signals with health data), inclusive communication and infrastructure needs for responding to extreme weather events, and integrated platforms visualizing air quality trends to support tailored, context-aware guidance beyond one-size-fits-all alerts. The workshop highlighted opportunities such as leveraging machine learning, Earth observation, and real-time surveillance to protect communities, but also noted barriers including data quality, technological redundancy, privacy, and governance challenges. Additionally, participants emphasized the need for interdisciplinary teams capable of collaborating across sectors, breaking down silos and addressing gaps in training and education. Overall, the workshop illustrates how process-driven, human-centered approaches can help surface user needs and generate testable prototype concepts, while underscoring the importance of direct community partnership for implementation.

Abadi, Azar M. [University of Alabama, Birmingham,↗
DOE OSTI2025

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2025

A2SD: Accelerating Scientific Innovation Through Autonomous Discovery Systems

The 2025 Advancing Autonomous Scientific Discovery (A2SD) workshop convened researchers from academia, national laboratories, and industry to explore the transformative role of autonomy in scientific discovery. The workshop highlighted a convergence of artificial intelligence, robotics, and computational workflows into autonomous systems capable of accelerating the scientific process. Presentations and discussions spanned autonomous experimentation, intelligent workflow orchestration, digital twins, and agent-based systems for managing complex research ecosystems. Key challenges discussed included interoperability across heterogeneous infrastructures, near real-time data management under FAIR principles, reproducibility, and the integration of human oversight. The workshop also emphasized the need for modular software interfaces, federated learning models, and education initiatives to support a next-generation scientific workforce.

Taufer, Michela [University of Tennessee, Knoxvill↗
DOE OSTI2024

Nanotwinned alloys under high pressure

Nanotwinned alloys are of interest due to their high strength and ductility, but twin boundaries may not be stable under shear. Computational studies indicate that high hydrostatic pressure may suppress detwinning mechanisms. Here, in this study, we investigate the microstructural changes of nanotwinned-nanocrystalline copper-nickel and Inconel 725 alloys under quasi-hydrostatic pressures up to 50 gigapascals (GPa). The alloys are compressed in a diamond anvil cell. In-situ x-ray diffraction (XRD) and ex-situ transmission electron microscopy (TEM) were employed to monitor microstructural changes. Twin boundary deformation and grain growth occur at 11.4 GPa quasi-hydrostatic pressure in the copper-nickel alloy. Molecular dynamics (MD) simulations reveal that hydrostatic pressure causes elevated local shear stress at grain boundaries, which leads to atomic rearrangements. A superposition of hydrostatic and deviatoric pressures lead to partial dislocation mediated twin boundary migration. In contrast, the Inconel 725 alloy showed stable twin and grain boundaries up to a quasi-hydrostatic pressure of 12.7 GPa. Texture, high solid solution strengthening, and low stacking fault energy are hypothesized to the enhanced microstructural stability in Inconel 725.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗
DOE OSTI2025

Environmental photosensitizers can exhibit enhanced actinic absorption in microhydrated clusters compared to solution

Brown carbon chromophores at environmental air-water interfaces often act as photosensitizers that absorb sunlight and subsequently transfer energy to nearby molecules, initiating a wide variety of chemical reactions. Despite their importance to understanding daytime chemistry at these airwater interfaces, little is known about the role of the solvation environment on the photophysical properties of these photosensitizers. In this work, we present a joint experimental-theoretical study of the vibrational and photophysical properties of microhydrated protonated and deprotonated 4-benzoylbenzoic acid (4-BBA), a key model system for environmental photosensitizers. We find that for protonated 4-BBAH + ∙(H 2 O) 0-1 , representing photosensitizers in very acidic conditions, a single bright state dominates the UV-Vis spectrum between 280 and 400 nm. Comparing the experimental UV-Vis spectra and quantum chemistry-predicted spectra of 4-BBA + ∙(H 2 O) 0-2 , we find that the degree of microhydration has little effect on the UV-Vis spectra or the orbitals of the dominant feature. For deprotonated 4-BBA ‒ , representing photosensitizers in basic conditions, quantum chemistry calculations predict that the UV-Vis spectra are ~3x weaker in intensity than the brightest 4-BBAH + ∙(H 2 O) 0-1 features and were not observed experimentally. Quantum chemistry calculations predict one dominant UV-Vis feature is present in the spectra of 4-BBA ‒ ∙(H 2 O) 0-2 , which exhibit minor shifts with degree of microhydration. While 4-BBA in bulk solution over a range of pH values has relatively weak absorption within the solar actinic region, we show that microhydrated 4-BBA has bright transitions within the actinic region. This indicates that the complex structure of environmental air-water interfaces can shift the absorption maximum of photosensitizers into the actinic region for enhanced absorption of sunlight and subsequent enhancement of photosensitizer-driven reactions.

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

Spectroscopy and dynamics of the v = 1 levels derived from the OH(D) stretching modes in the I‾∙HDO complex using CW and time-resolved, resonant two-photon infrared excitation of the cryogenically cooled ions

The vibrational energy levels of the two isotopomers adopted by the I‾∙HDO ion-molecule complex occur such that the OH(D) stretching fundamentals span its dissociation energy, thus enabling a spectroscopic investigation of the dynamics displayed by a system prepared in the vicinity of the dissociation threshold. This regime is explored using infrared photoexcitation of the mass-selected complexes cooled in a cryogenic radiofrequency (Paul) ion trap. Survey spectra are obtained at modest resolution using two-color, IR-IR photodissociation with nanosecond lasers to establish the level structure and unimolecular decay dynamics of the v = 1 and 2 levels of the bound OH(D) oscillator. The v = 1 levels are prepared by fixed frequency excitation in the trap and the absorption spectra arising from this excited state are probed by photofragmentation of the complex with a second pulsed IR laser after a variable delay time (0 to 20 ms). The transitions to levels above the dissociation threshold for I‾ + HDO formation are observed to be sharp (~5 cm -1 FWHM). At low pressure, the bound OD (v = 1) population relaxes very slowly (~3 ms), consistent with resonant fluorescence in the IR. The collisional quenching rate constants of this level by the He buffer gas were estimated to be on the order 3 x 10 -12 cm 3 /s based on a crude Stern-Volmer analysis. Here, the rotational fine structure and linewidths (≲0.01 cm -1 FWHM) arising from transitions of the non-bonded OH stretch fundamental of the OD-bound isotopomer that lies just above the dissociation energy are determined using single photon photodissociation by excitation of the 10 K ions with a single-frequency, CW IR laser in the ion trap.

Cryogenic ion infrared spectroscopy↗
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗