Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse modeling”

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

Neutrino and gamma-ray emissions from NGC 1068

IceCube has recently reported the detection of ∼1–10 TeV neutrinos from the nearby active galaxy, NGC 1068. The lack of TeV-scale emission from this source suggests that these neutrinos are generated in the dense corona that surrounds NGC 1068’s supermassive black hole. In this paper, we present a physical model for this source, including the processes of pair production, pion production, synchrotron, and inverse Compton scattering. We have also performed a new analysis of Fermi-LAT data from the direction of NGC 1068, finding that the gamma-ray emission from this source is very soft but bright at energies below ∼1 GeV . Our model can predict a gamma-ray spectrum that is consistent with Fermi-LAT observations when the magnetic field within the corona of this active galactic nucleus (AGN) is quite high, namely 𝐵 ≳6 kG . To explain the observed neutrino emission, this source must accelerate protons with a total power that is comparable to its intrinsic x-ray luminosity. In this context, we consider two additional nearby active galaxies, NGC 4151 and NGC 3079, which have been identified as promising targets for IceCube.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Parameterization of Turbulent Dissipation and Pressure Damping Time Scales in Stably Stratified Inversions, and its Effects on Low Clouds in Global Simulations

Abstract It is difficult for coarse‐resolution global models of the atmosphere to accurately simulate the observed distribution of low clouds. In particular, it is difficult for moist turbulence closure models to simulate sufficiently bright near‐coastal stratocumulus (Sc) without simulating overly bright marine shallow cumuli (Cu). To parameterize bright Sc, a turbulence parameterization must damp the turbulent fluxes of heat and moisture above cloud top in order to prevent excessive entrainment of dry air into cloud top. To parameterize dim shallow Cu, the subgrid variances of temperature and moisture must remain large, in order to permit partial cloudiness. However, damping the fluxes but not the variances just above cloud top is difficult if a parameterization uses a single “master” time scale to damp both. In nature, the above‐cloud fluxes are damped by pressure fluctuations, whereas scalar variances are damped by a different process, namely, turbulent dissipation. In a stably stratified inversion above cloud, pressure damping is large but turbulent dissipation is small. To avoid this problem, a multitime scale parameterization for damping has been developed. The damping parameterization has been implemented in a global model and evaluated. The parameterization is capable of dimming shallow Cu while producing adequately bright Sc.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Deep Learning for Simultaneous Inference of Hydraulic and Transport Properties

Abstract Identification of a heterogeneous conductivity field and reconstruction of a contaminant release history are key aspects of subsurface remediation. These two goals are achieved by combining model predictions with sparse and noisy hydraulic head and concentration measurements. Solution of this inverse problem is notoriously difficult due to, in part, high dimensionality of the parameter space and high computational cost of repeated forward solves. We use a convolutional adversarial autoencoder (CAAE) to parameterize a heterogeneous non‐Gaussian conductivity field via a low‐dimensional latent representation. A three‐dimensional dense convolutional encoder‐decoder (DenseED) network serves as a forward surrogate of the flow and transport model. The CAAE‐DenseED surrogate is fed into the ensemble smoother with multiple data assimilation (ESMDA) algorithm to sample from the Bayesian posterior distribution of the unknown parameters, forming a CAAE‐DenseED‐ESMDA inversion framework. The resulting CAAE‐DenseED‐ESMDA inversion strategy is used to identify a three‐dimensional contaminant source and conductivity field. A comparison of the inversion results from CAAE‐ESMDA with physical flow and transport simulator and from CAAE‐DenseED‐ESMDA shows that the latter yields accurate reconstruction results at the fraction of the computational cost of the former.

Zhou, Zitong↗
DOE OSTI2023

Effects of Size and Shape on the Tolerances for Misalignment and Probabilities for Successful Oriented Attachment of Nanoparticles

Oriented attachment (OA) of nanoparticles is an important pathway of crystal growth, but tools for quantitatively modeling OA are lacking. Here we present several simple models that relate the probability of achieving OA to basic geometric parameters such as particle size, shape, and lattice periodicity. A Moiré-domain model is applied to understand twist-misorientations between parallel surfaces, and it predicts that the range of twist angles yielding perfect OA is inversely related to the width of the contact area. This is confirmed and further developed using a surface functional model, which predicts how crystallographic registration forces drive the emergence of complex orientational energy landscapes. The energy landscapes are predicted to possess local minima that can trap particles in imperfect alignments, and these local minima become deeper and more numerous as the contact area increases, making OA more challenging for large particles. Further, a second set of models is presented to understand the sequence of events by which two crystallographic faces become co-planer after collision. We use a ‘central force approximation’ to quantitatively predict the odds of attaining coalignment between various faces when particles collide with random misalignments, and we show that in the absence of biasing forces, the probability of attaining alignment on a given face is roughly proportional to its solid angle as viewed from the center of the particle. The model predicts that OA is most favorable between well-faceted particles and becomes exceedingly unlikely for large spherical particles that express many microfacets.

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

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

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

On the Neutrino and Gamma-Ray Emission from NGC 1068

IceCube has recently reported the detection of $\sim 1-10 \,{\rm TeV}$ neutrinos from the nearby active galaxy, NGC 1068. The lack of TeV-scale emission from this source suggests that these neutrinos are generated in the dense corona that surrounds NGC 1068's supermassive black hole. In this paper, we present a physical model for this source, including the processes of pair production, pion production, synchrotron, and inverse Compton scattering. We have also performed a new analysis of Fermi-LAT data from the direction of NGC 1068, finding that the gamma-ray emission from this source is very soft but bright at energies below $\sim 1 \, {\rm GeV}$. Our model can predict a gamma-ray spectrum that is consistent with Fermi-LAT observations, but only if the magnetic field within the corona of this active galactic nucleus (AGN) is quite high, namely $B\gtrsim 6 \, {\rm kG}$. To explain the observed neutrino emission, this source must accelerate protons with a total power that is comparable to its intrinsic X-ray luminosity. In this context, we consider two additional nearby active galaxies, NGC 4151 and NGC 3079, which have been identified as promising targets for IceCube.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Utah FORGE: Seismic Velocity Models, February 2021

This dataset contains a map, showing the Utah FORGE seismic stations, and seismic velocity model data. There are 61 1-D velocity models which are in a compressed TAR file. A paper is referenced at the end of this description which discusses the use of these data in 3D modelling. The paper summary follows: We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2-5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ~22 degrees to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity. Reference: Zhang, H. and K. L. Pankow (2021). High-resolution Bayesian spatial auto-correlation (SPAC) pseudo-3D Vs model of Utah FORGE site with a dense geophone array, Geophys. Res. Int, https://doi.org/10.1093/gji/ggab049

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

RSPECT: A PROSPECT-based model incorporating the real structure of rice leaves

Radiative transfer models (RTMs) describe how light is absorbed, scattered, and transmitted within leaves by simulating mechanistic light propagation processes. The PROSPECT model is based on measurable parameters (the leaf biochemical content) and a non-measurable parameter (the leaf anatomical structure represented by the leaf structure parameter (N)). The effect of N on the optical properties of leaves has been investigated through a number of local and global sensitivity analyses. Other studies have directly evaluated the effect of the leaf anatomical structure on spectral reflectance, particularly in the near infrared region. However, the relationship between N and the anatomical structure is unclear. Here, in this study, we leveraged eLeaf, a ray tracing-based 3D rice leaf simulator, to establish relationships between leaf anatomical features and spectral properties, enabling us to replace N in the PROSPECT-4 model with measurable leaf anatomical parameters and develop the RSPECT model. The leaf thickness at minor vein, leaf thickness at bulliform cells, mesophyll thickness at minor vein, and distance between two minor veins could be used to predict N effectively. The RSPECT model achieved spectral simulation accuracy comparable to PROSPECT-4 and was more suitable for parameter inversion of the physical and chemical properties of rice leaves, with relative root mean square errors of 7.4% for chlorophyll content, 5.6% for equivalent water thickness, and 7.5% for dry matter content. In conclusion, RSPECT improves radiative transfer modeling by integrating measurable anatomical features and provides a framework for extending this approach to other vegetation types.

hyperspectral↗
DOE OSTI2022

Theoretical and numerical studies of inverse source problem for the linear parabolic equation with sparse boundary measurements

We consider the inverse source problem in the parabolic equation, where the unknown source possesses the semi-discrete formulation. Theoretically, we prove that the flux data from any nonempty open subset of the boundary can uniquely determine the semi-discrete source. This means the observed area can be extremely small, and that is the reason we call it sparse boundary data. For the numerical reconstruction, we formulate the problem from the Bayesian sequential prediction perspective and conduct the numerical examples which estimate the space-time-dependent source state by state. To better demonstrate the method’s performance, we solve two common multiscale problems from two models with a long source sequence. The numerical results illustrate that the inversion is accurate and efficient.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Investigating effects of soil chemicals on density of small mammal bioindicators using spatial capture-recapture models

Monitoring the ecological impacts of environmental pollution and the effectiveness of remediation efforts requires identifying relationships between contaminants and the disruption of biological processes in populations, communities, or ecosystems. Wildlife are useful bioindicators, but traditional comparative experimental approaches rely on a staunch and typically unverifiable assumption that, in the absence of contaminants, reference and contaminated sites would support the same densities of bioindicators, thereby inferring direct causation from indirect data. We demonstrate the utility of spatial capture-recapture (SCR) models for overcoming these issues, testing if community density of common small mammal bioindicators was directly influenced by soil chemical concentrations. By modeling density as an inhomogeneous Poisson point process, we found evidence for an inverse spatial relationship between Peromyscus density and soil mercury concentrations, but not other chemicals, such as polychlorinated biphenyls, at a site formerly occupied by a nuclear reactor. Although the coefficient point estimate supported Peromyscus density being lower where mercury concentrations were higher (β = –0.44), the 95% confidence interval overlapped zero, suggesting no effect was also compatible with our data. Estimated density from the most parsimonious model (2.88 mice/ha; 95% CI = 1.63–5.08), which did not support a density-chemical relationship, was within the range of reported densities for Peromyscus that did not inhabit contaminated sites elsewhere. Environmental pollution remains a global threat to biodiversity and ecosystem and human health, and our study provides an illustrative example of the utility of SCR models for investigating the effects that chemicals may have on wildlife bioindicator populations and communities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Sensitivity of Joint Inversions of Seismic and Geodynamic Data to Mantle Viscosity

Seismic tomography has revealed the existence of large-scale velocity heterogeneities in the mantle. The interpretation of seismic velocity anomalies in terms of temperature and chemical composition is nonunique. We use geodynamic observations including gravity, plate motions, dynamic topography, and excess ellipticity of the core-mantle boundary combined with seismic observations to investigate the thermo-chemical structure of the mantle through joint inversions. An outstanding issue, however, is the physical connection between mantle density anomalies and the surface geodynamic observations, which requires knowledge of the mantle viscosity structure. Here we perform joint inversions assuming different viscosity profiles and examine the dependence of the results on the viscosity. We first assume that mantle heterogeneity is due to thermal variations, which places a constraint on the relation between seismic velocity and density, and we subsequently relax the constraint to allow for potential nonthermal effects. In all of our joint inversions, a nonthermal origin of density anomalies is required to explain the geodynamic data, though the amount varies with the assumed viscosity structure. A common observation is a high-density chemical signal in the center of the large low-shear-velocity provinces at the base of the mantle resulting in a near neutral or slightly dense overall buoyancy there. Using the derived density models and their corresponding viscosity profiles, we also calculate instantaneous mantle flow fields. The predicted flow fields derived from joint inversions are generally similar but are quite different from flow fields using density models derived from a posteriori scaling of pure seismic tomography models.

58 GEOSCIENCES↗
DOE OSTI2020

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗
DOE OSTI2024

A Mechanism-Based Reaction–Diffusion Model for Accelerated Discovery of Thermoset Resins Frontally Polymerized by Olefin Metathesis

Frontal ring-opening metathesis polymerization (FROMP) involves a self-perpetuating exothermic reaction, which enables the rapid and energy-efficient manufacturing of thermoset polymers and composites. Current state-of-the-art reaction–diffusion FROMP models rely on a phenomenological description of the olefin metathesis kinetics, limiting their ability to model the governing thermo-chemical FROMP processes. Furthermore, the existing models are unable to predict the variations in FROMP kinetics with changes in the resin composition and as a result are of limited utility toward accelerated discovery of new resin formulations. In this work, we formulate a chemically meaningful model grounded in the established mechanism of ring-opening metathesis polymerization (ROMP). Our study aims to validate the hypothesis that the ROMP mechanism, applicable to monomer-initiator solutions below 100 °C, remains valid under the nonideal conditions encountered in FROMP, including ambient to >200 °C temperatures, sharp temperature gradients, and neat monomer environments. Through extensive simulations, we demonstrate that our mechanism-based model accurately predicts the FROMP behavior across various resin compositions, including polymerization front velocities and thermal characteristics (e.g., T max ). Additionally, we introduce a semi-inverse workflow that predicts FROMP behavior from a single experimental data point. Notably, the physiochemical parameters utilized in our model can be obtained through DFT calculations and minimal experiments, highlighting the model’s potential for rapid screening of new FROMP chemistries in pursuit of thermoset polymers with superior thermo-chemo-mechanical properties.

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

Mercury and Radiocesium Accumulation and Associations With Sublethal Endpoints in the Florida Green Watersnake ( Nerodia floridana )

Mercury (Hg) and radiocesium ( 137 Cs) are well-known environmental contaminants with the potential to impact the health of humans and wildlife. Snakes have several characteristics conducive to studying environmental contamination but have rarely been included in the monitoring of polluted sites. We investigated the bioaccumulation of Hg and 137 Cs and associations with sublethal effects (standard metabolic rate [SMR] and hemoparasite infections) in Florida green watersnakes (Nerodia floridana). We captured 78 snakes from three former nuclear cooling reservoirs on the US Department of Energy's Savannah River Site in South Carolina (USA). For captured snakes, we (1) determined whole-body 137 Cs, (2) quantified total Hg (THg) using snake tail clips, (3) conducted hemoparasite counts, and (4) measured the SMR. We used multiple regression models to determine associations among snake body size, capture location, sex, tail THg, whole-body 137 Cs, Hepatozoon spp. prevalence and parasitemia, and SMR. Average whole-body 137 Cs (0.23 ± 0.08 Becquerels [Bq]/g; range: 0.00–1.02 Bq/g) was correlated with snake body size and differed significantly by capture site (Pond B: 0.67 ± 0.05 Bq/g; Par Pond: 0.10 ± 0.02 Bq/g; Pond 2: 0.03 ± 0.02 Bq/g). Tail THg (0.33 ± 0.03 mg/kg dry wt; range: 0.16–2.10 mg/kg) was significantly correlated with snake body size but did not differ by capture site. We found no clear relationship between SMR and contaminant burdens. However, models indicated that the prevalence of Hepatozoon spp. in snakes was inversely related to increasing whole-body 137 Cs burdens. Our results indicate the bioaccumulation of Hg and 137 Cs in N. floridana and further demonstrate the utility of aquatic snakes as bioindicators. Furthermore, our results also suggest a decrease in Hepatozoon spp. prevalence related to increased burdens of 137 Cs. Although the results are intriguing, further research is needed to understand the dynamics between 137 Cs and Hepatozoon spp. infections in semiaquatic snakes.

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

Enhanced neutrino polarizability

We point out that neutrinos can have enhanced couplings to photons, if light (pseudo)scalar mediators are present, resulting in a potentially measurable neutrino polarizability. We show that the expected suppression from small neutrino masses can be compensated by the light mediator mass, generating dimension 7 Rayleigh operators at low scales. We explore the rich phenomenology of such models, computing in detail the constraints on the viable parameter space, spanned by the couplings of the mediator to neutrinos and photons. Finally, we build several explicit models that lead to an enhanced neutrino polarizability by modifying the inverse see-saw majoron, i.e., the pseudo-Nambu-Goldstone boson of the U(1) L global lepton number responsible for generating small neutrino masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Impact of Changing Winds on the Mauna Loa CO 2 Seasonal Cycle in Relation to the Pacific Decadal Oscillation

Abstract Long‐term measurements at the Mauna Loa Observatory (MLO) show that the CO 2 seasonal cycle amplitude (SCA) increased from 1959 to 2019 at an overall rate of 0.22 0.034 ppm decade −1 while also varying on interannual to decadal time scales. These SCA changes are a signature of changes in land ecological CO 2 fluxes as well as shifting winds. Simulations with the TM3 tracer transport model and CO 2 fluxes from the Jena CarboScope CO 2 Inversion suggest that shifting winds alone have contributed to a decrease in SCA of −0.10 0.022 ppm decade −1 from 1959 to 2019, partly offsetting the observed long‐term SCA increase associated with enhanced ecosystem net primary production. According to these simulations and MIROC‐ACTM simulations, the shorter‐term variability of MLO SCA is nearly equally driven by varying ecological CO 2 fluxes (49%) and varying winds (51%). We also show that the MLO SCA is strongly correlated with the Pacific Decadal Oscillation (PDO) due to varying winds, as well as with a closely related wind index (U‐PDO). Since 1980, 44% of the wind‐driven SCA decrease has been tied to a secular trend in the U‐PDO, which is associated with a progressive weakening of westerly winds at 700 mbar over the central Pacific from 20°N to 40°N. Similar impacts of varying winds on the SCA are seen in simulations at other low‐latitude Pacific stations, illustrating the difficulty of constraining trend and variability of land CO 2 fluxes using observations from low latitudes due to the complexity of circulation changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Optically induced umklapp shift currents in striped cuprates

Motivated by recent experiments that observed low-frequency second-order optical responses in doped striped superconductors, here we investigate the nonlinear electrodynamics of systems exhibiting a charge density wave (CDW) order parameter. Due to the Bragg scattering off the CDW order, an incoming spatially homogeneous electric field in addition to zero momentum current generates umklapp currents that are modulated in space at momenta of the reciprocal CDW lattice. In particular, here we predict and microscopically evaluate the umklapp shift current, a finite momentum analog of the regular shift current which represents the second-order optical process that downconverts homogeneous AC electric field into low-frequency, zero momentum current. Specifically, we evaluate real-time response functions within mean-field theory via the Keldysh technique and use the Peierls substitution to compute observables at finite momenta in lattice models. We find that systems with certain lattice symmetries (such as inversion symmetry), where the regular shift current is disallowed, may give rise to the umklapp one. We apply our framework to investigate lattice symmetries in layered materials with helical-like stripes and show that both types of shift currents provide insight into the nature of intertwined phases of matter. Finally, we discuss the relation of our findings to recent experiments in striped superconductors. Published by the American Physical Society 2024

Materials Science↗