Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sensitivity coefficients”

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 361 records · Page 20

Slow-Light Mid-IR Silicon Photonic Chips for NO 2 and CH 4 Gas Detection

A compact, chip-scale mid-infrared gas sensor is demonstrated, leveraging a two-dimensional photonic crystal waveguide (PCW) fabricated on a silicon-on-insulator (SOI) platform. The PCW comprises a hexagonal lattice with lattice constant a = 860 nm and hole radius r = 0.22a, incorporating a central line defect of reduced-radius holes (r s = 0.7r) to induce slow-light propagation near the photonic band edge with a group index of approximately 73, thereby enhancing light-matter interaction. The sensor operates at fundamental absorption wavelengths of 3.42 μm for nitrogen dioxide (NO 2 ) and 3.40 μm for methane (CH 4 ), utilizing the strongest molecular vibrational transitions for maximum sensitivity. Experimental validation was conducted using dynamically diluted gas mixtures generated by mass flow controllers, with signal acquisition performed by a liquid nitrogen-cooled InSb detector. For NO 2 , the sensor exhibited excellent linear response over 5–25 ppm (part per million) with coefficient of determination R 2 = 0.9934, achieving a detection limit of 210 ppb (part per billion)─representing the first reported silicon photonic-based NO 2 detection. For CH 4 , exposure to 25 ppm resulted in a 6.4% decrease in transmitted intensity, demonstrating multigas sensing capability. The CMOS-compatible fabrication process and compact 3 mm device footprint establish this SOI-PCW platform as a scalable, low-power solution for integrated mid-infrared gas sensing, with significant potential for environmental monitoring and industrial safety applications.

Crystals↗

Sensitivity study on the fin effect of thermocouple mounted on the heated surface under film boiling condition

In the reactivity-initiated accident (RIA) scenario, right after the power burst with the occurrence of departure from nucleate boiling (DNB), the cladding temperature quickly reaches its maximum and is maintained in a steady film boiling regime. Then, under limited heat transfer performance of film boiling, the local cooling effect by the thermocouple wire may not be neglected. A computational model is developed to estimate the local cooling effect of the thermocouple wire mounted on the cladding surface under film boiling condition. The model predicts the temperature profile along the cladding surface as physically expected and is verified against the benchmark problem. Sensitivity studies reveal the significant effects of the thickness of vapor film layer and the diameter of thermocouple on the temperature drop near the thermocouple, for they increase the heat transfer area to the liquid water region. The impact of film boiling heat transfer coefficient, however, is not significant, as the primary heat transfer mechanism near the thermocouple is the heat conduction to the wire. Future study will include the contribution of radiative heat transfer between the cladding surface and the vapor film layer, as well as a transient power history which represents an actual RIA scenario.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessing free tropospheric quasi-equilibrium for different GCM resolutions using a cloud-resolving model simulation of tropical convection

Abstract This study examines the free-tropospheric quasi-equilibrium at different global climate model (GCM) resolutions using the simulation of tropical convection by a cloud-resolving model during the Tropical Western Pacific International Cloud Experiment. The simulated dynamic and thermodynamic fields within the model domain are averaged over subdomains of different sizes equivalent to different GCM resolutions. These coarse-grained fields are then used to compute CAPE and its change with time, and their relationships with simulated convection. Results show that CAPE change with time is controlled predominantly by variations of thermodynamic properties in the planetary boundary layer for all subdomain sizes ranging from 64 to 4 km. Lag correlation analysis shows that CAPE generation by the free-tropospheric dynamical advection (dCAPE ls ) leads convective precipitation but is in phase with convective mass flux at 600 mb and 500 mb vertical velocity for all subdomain sizes. However, the correlation coefficients and regression slopes decrease as the subdomain size decreases for subdomain sizes smaller than 16 km. This is probably due to increased randomness of convection and more scale-dependence of the relationships when the subdomain size reaches the grey zone. By examining the sensitivity of the relationships of convection with dCAPE ls to temporal scales in different subdomain size, it shows that the quasi-equilibrium between dCAPE ls and convection holds well for timescales of 30 min or longer at all subdomain sizes. These results suggest that the free tropospheric quasi-equilibrium assumption may still be useable even for GCM resolutions in the grey zone.

54 ENVIRONMENTAL SCIENCES↗

Characterizing how size distribution and concentration affect echogenicity of ultrasound contrast agents

Here, we investigated the effects of UCA gas bubble size distribution and concentration on the generated ultrasound echogenicity signal. Gas bubble size characterization using Coulter Counter and cryogenic-SEM revealed the hollow structure and rare presence of microbubbles >10 µm in a commercial UCA product, Lumason™. Volume-weighed size and concentration were observed to be more sensitive to changes in UCA bubble stability than number-weighted size and concentration. Size distribution measurements showed that the force (e.g., shaking/agitation energy) used to redisperse the sample did not affect the size distribution, concentration, or echogenicity of the UCA sample. The ultrasound backscattering coefficient (BSC) of size fractionated and serial diluted microbubbles showed that the echogenicity signal correlates most with UCA bubble concentration, especially volume-weighted concentration. Findings from this study may be used to support demonstrating the equivalence of a generic UCA product to the reference listed drug.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Anomalous thermal effect in Zr Te 5 observed via photothermal measurements

In this study, we explore the magneto-thermoelectric power (MTP) of $ZrTe$ 5 , a canonical Dirac semimetal, through a novel photothermal technique. Unlike conventional thermoelectric studies that rely on on-chip heaters and are limited by fabrication processes, especially for stress-sensitive materials, our approach utilizes photothermal effects to induce temperature gradients. Our experiments, applying a magnetic field approximately parallel and transverse to the photocurrent detection direction, reveal that the photothermal method efficiently and reliably extracts both diagonal and off-diagonal components of the thermoelectric coefficient of $ZrTe$ 5 . Here, we observe that the longitudinal MTP reproduces features previously reported in thermal transport studies, while the photoinduced transverse MTP confirms the anomalous Nernst effect. This photothermal measurement technique opens new avenues for investigating transport properties in a wide range of quantum materials, both in 3D and 2D systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reanalysis of the 𝛽− $\overline{v}$ 𝑒 Angular Correlation Measurement from the aSPECT Experiment with New Constraints on Fierz Interference

On the basis of revisions of some of the systematic errors, we reanalyzed the electron-antineutrino angular correlation (𝑎 coefficient) in free neutron decay inferred from the recoil energy spectrum of the protons which are detected in 4⁢𝜋 by the aSPECT spectrometer. With 𝑎=−0.104 02⁢(82), the new value differs only marginally from the one published in 2020. Here, the experiment also has sensitivity to 𝑏, the Fierz interference term. From a correlated (𝑏,𝑎) fit to the proton recoil spectrum, we derive a limit of 𝑏 =−0.0098⁢(193) which translates into a somewhat improved 90% confidence interval region of −0.041 ≤ 𝑏 ≤ 0.022 on this hypothetical term. Tighter constraints on 𝑏 can be set from a combined analysis of the PERKEO III (𝛽 asymmetry) and aSPECT measurement which suggests a finite value of 𝑏 with 𝑏 (𝑐) =−0.0181 ± 0.0065 deviating by 2.82⁢𝜎 from the standard model.

Beta decay↗

Electronic 1/f noise as a probe of dimensional effects on spin-glass dynamics

Over the past decade, spin-glass simulations have improved to the point that they now access time- and length-scales comparable to experiments at the mesoscale. A recent series of thin-film field-cooled/zero-field-cooled magnetization (FC/ZFC) experiments demonstrated activated spin dynamics, with a temperature-independent activation energy proportional to the logarithm of the film thickness and with coefficients in remarkable agreement with the simulation. These measurements require the application of small magnetic fields, which has been shown to affect the spin-glass energy landscape. Measurements of the 1/f noise in metallic spin-glasses have been previously shown to be a sensitive probe of the spin dynamics, and the measurements can be made without applying a magnetic field. In this mini-review, we review these techniques and discuss how transport measurements can fit into the current landscape of spin-glass measurements. We compare previous measurements to more recent measurements on similar films, made with ostensibly different cooling protocols, and compare both the previous and recent measurements to the magnetometry. The transport measurements—taken over a wider range of temperature than magnetometry—suggest that the maximum spin-glass energy barrier height is temperature-dependent, not fixed, possibly due to two-dimensional dynamics. We discuss this possibility, along with future measurements, which may be able to resolve this mystery.

1/f noise↗

Four Isotope-Labeled Recombination Pathways of Ozone Formation

A theoretical approach is developed for the description of all possible recombination pathways in the ozone forming reaction, without neglecting any process a priori, and without decoupling the individual pathways one from another. These pathways become physically distinct when a rare isotope of oxygen is introduced, such as 18 O, which represents a sensitive probe of the ozone forming reaction. Each isotopologue of O 3 contains two types of physically distinct entrance channels and two types of physically distinct product wells, creating four recombination pathways. Calculations are done for singly and doubly substituted isotopologues of ozone, eight rate coefficients total. Two pathways for the formation of asymmetric ozone isotopomer exhibit rather different rate coefficients, indicating large isotope effect driven by ΔZPE-difference. Rate coefficient for the formation of symmetric isotopomer of ozone (third pathway) is found to be in between of those two, while the rate of insertion pathway is smaller by two orders of magnitude. These trends are in good agreement with experiments, for both singly and doubly substituted ozone. The total formation rates for asymmetric isotopomers are found to be somewhat larger than those for symmetric isotopomers, but not as much as in the experiment. Altogether, the distribution of lifetimes is found to be very similar for the metastable states in symmetric and asymmetric ozone isotopomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

Lorentz invariance violation (LIV) is often described by dispersion relations of the form E i 2 = m i 2 +p i 2 +δ i,n E 2+n with delta different based on particle type i , with energy E, momentum p and rest mass m. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients δ i,n tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond 10 19 eV, we obtain δ γ,0 > -10 -21 , δ γ,1 > -10 -40 eV -1 and δ γ,2 > -10 -58 eV -2 in the case of a subdominant proton component up to 10 20 eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as δ had,0 < 10 -19 , δ had,1 < 10 -38 eV -1 and δ had,2 < 10 -57 eV -2 at 5σ CL.

79 ASTRONOMY AND ASTROPHYSICS↗

Iterative Bragg peak removal on X-ray absorption spectra with automatic intensity correction

This study introduces a novel iterative Bragg peak removal with automatic intensity correction (IBR-AIC) methodology for X-ray absorption spectroscopy (XAS), specifically addressing the challenge of Bragg peak interference in the analysis of crystalline materials. The approach integrates experimental adjustments and sophisticated post-processing, including an iterative algorithm for robust calculation of the scaling factor of the absorption coefficients and efficient elimination of the Bragg peaks, a common obstacle in accurately interpreting XAS data, particularly in crystalline samples. The method was thoroughly evaluated on dilute catalysts and thin films, with fluorescence mode and large-angle rotation. The results underscore the technique's effectiveness, adaptability and substantial potential in improving the precision of XAS data analysis. While demonstrating significant promise, the method does have limitations related to signal-to-noise ratio sensitivity and the necessity for meticulous angle selection during experimentation. Overall, IBR-AIC represents a significant advancement in XAS, offering a pragmatic solution to Bragg peak contamination challenges, thereby expanding the applications of XAS in understanding complex materials under diverse experimental conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generation of the Thermal Scattering Law of Uranium Dioxide with Ab Initio Lattice Dynamics to Capture Crystal Binding Effects on Neutron Interactions

Scattering of thermal neutrons and Doppler broadening of epithermal neutron resonances in uranium and its compounds may be sensitive to crystal binding. The thermal scattering law (TSL) for uranium dioxide, which captures crystal binding effects, has been reevaluated for ENDF/B-VIII.0. Phonon spectra were generated using ab initio lattice dynamics for the paramagnetic phase and validated against experiment. Improved agreement with the Debye-Waller coefficient as a function of temperature is found relative to the spectrum used for the ENDF/B-VII.1 evaluation. The TSL was generated using the phonon expansion method within the NJOY nuclear data processing package and was found to be in reasonable agreement with inelastic neutron scattering measurements. The present evaluation predicts a reduction in the inelastic scattering cross section relative to ENDF/B-VII.1 and a total scattering cross section consistent with neutron transmission experiments

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploring the influence of bulk viscosity of QCD on dilepton tomography

The collective behavior of hadrons and of electromagnetic radiation in heavy-ion collisions has been widely used to study the properties of the quark-gluon plasma (QGP). Indeed this collectivity, as measured by anisotropic flow coefficients, can be used to constrain the transport properties of QGP. The goal of this contribution is to investigate the influence of the specific bulk viscosity (ζ/s) on dilepton production, both at Relativistic Heavy-Ion Collider and Large Hadron Collider energies. We explore the sensitivity of dileptons to dynamical features that bulk viscosity induces on the evolution of a strongly interacting medium, and highlight what makes them a valuable probe in the pursuit to also constrain ζ/s.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Landau-Devonshire Approach of Elastic Strain Effects on Ferroelectric Al1-xScxN

Landau-Devonshire thermodynamic modeling enables the phenomenological description of ferroelectrics. 1 The approach has been validated across a wide range of ferroelectric materials from the first ferroelectric discovered, Rochelle Salt, to recent HfO2-based materials. 2-4 Extension to what is now referred to as the Landau-Ginsburg-Devonshire model enables rigorous investigation of polarization boundaries and reorientation, phase transformations, and domain formation/evolution. We present a thermodynamic analysis of the recently discovered nitride ferroelectric materials 5-7 using the classic Landau-Devonshire approach. The electrostrictive and dielectric stiffness coefficients of Al1-xScxN with wurtzite structure (6mm) are determined using a free energy density function assuming a hexagonal parent phase (6/mmm), with the first order phase transition based on the dielectric stiffness relationships (Fig. 1). The results of this analysis show that the strain sensitivity of the energy barrier is one order of magnitude larger than that of the spontaneous polarization in these wurtzite ferroelectrics, yet both are less sensitive to strain compared to classic perovskite ferroelectrics (Fig. 2). These analysis results reported here explain experimentally reported sensitivity of coercive field to elastic strain/stress in Al1-xScxN films, and would enable further thermodynamic analysis via phase field simulation and related methods.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗

Estimating the color lifetime of energetic quarks

Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.

Electron-ion collider↗

Plasma decay in hydrocarbons and hydrocarbon- and H 2 O-containing mixtures excited by high-voltage nanosecond discharge at elevated gas temperatures

Plasma decay after a high-voltage nanosecond discharge was experimentally and numerically studied in pure hydrocarbons (C 2 H 6 and C 3 H 8 ), and H 2 O:N 2 and C 3 H 8 :O 2 mixtures for pressures in the range 2–4 Torr and gas temperatures from 300 to 600 K. In a stoichiometric C 3 H 8 :O 2 mixture, plasma decay was also studied in a repetitively pulsed discharge for varing numbers of discharge pulses (varying degrees of fuel oxidation). The rate of plasma decay was determined from the temporal evolution of electron density measured using the microwave interferometer. It was observed that gas heating to 600 K leads to a decrease in the rate of plasma decay in all cases. The effect of heating on plasma decay was most profound in the H 2 O:N 2 mixture (after a single discharge pulse) and in the C 3 H 8 :O 2 mixture for high degrees of fuel oxidation. A kinetic scheme was developed to numerically simulate the plasma decay in hydrocarbons and combustible mixtures. Numerical analysis showed that, under the conditions studied, plasma decay was controlled by dissociative electron recombination with simple molecular and cluster ions. Gas heating led to a decrease in the rate of the electron-ion recombination and the rate of conversion of molecular ions to cluster ions. As a result, the gas temperature increase caused a decrease in the fraction of cluster ions for which the recombination coefficients are an order of magnitude higher than the recombination coefficients for molecular ions. The influence of gas heating on the decrease of the amount of cluster ions was more important when the ion composition was dominated by hydrated H 3 O + (H 2 O) k ions. The rates of the formation of these ions are extremely sensitive to any variations in gas temperature. A result, in agreement with our observations, gas heating led to an anomalous decrease in the rates of plasma decay in the H 2 O:N 2 mixture, as well as in the C 3 H 8 :O 2 mixture when H 2 O molecules were produced due to fuel oxidation.

20 FOSSIL-FUELED POWER PLANTS↗

New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles

Abstract. The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.

Harp, Dylan R. (ORCID:0000000197778000)↗

Predicting Device Parameters for Dye-Sensitized Solar Cells from Electronic Structure Calculations to Reproduce Experiment.

Given that improvements to the power-conversion efficiency (PCE) of dye-sensitized solar cells (DSSCs) have slowed in recent years, a means to accurately predict the device parameters yielded by trial dyes in silico, without having to synthesize them, would be extremely valuable to speed up the design process. Currently, the best-performing methods of calculating device parameters rely on a set of experimentally determined kinetic coefficients. In practice, it is very difficult to measure these kinetic parameters accurately, limiting the overall accuracy of such predictive methods. This work proposes a model to obtain key parameters such as J(SC), V-OC, and PCE using only the results from density functional theory (DFT) and time-dependent DFT calculations, noting that rates of electron-transfer steps are ultimately linked to the electronic structure of the dye center dot center dot center dot TiO2 working electrode. Six organic DSSC dyes from dissimilar chemical classes (L0, L1, L2, WS-2, WS-92, and C281) were chosen to demonstrate the power of this approach. Their a priori known experimentally determined device performance metrics served to validate our predictions. The greatest absolute error in our predicted PCE values was 0.36% relative to the experiment, while the greatest fractional error was 0.042. This indicates that the proposed model offers a dramatic improvement on previous predictive methods for DSSC device parameters, both in accuracy and in consistency. Moreover, such a predictive model has great potential to be applied to other photovoltaic applications, further enabling the design of novel, highly efficient photoactive materials.

density functional theory↗