Search NASA⌕ Search

SEARCH · Search NASA

Results for “Real Analytic Functions”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Real-Time Atomic-Scale Structural Analysis Resolves the Amorphous to Crystalline CaCO 3 Mechanism Controversy

Amorphous calcium carbonate (ACC) occurs as a precursor to geological and biogenic calcium carbonate (CaCO 3 ), yet its transformation pathways and reaction mechanisms remain inconsistent and controversial. In this study, we investigated the transformation of ACC to calcite under both solution and dry conditions, in the presence and absence of impurity ions, utilizing operando time-resolved synchrotron X-ray diffraction (TRXRD) and reactive transport modeling. Results demonstrate that TRXRD techniques allow us to differentiate dissolution-reprecipitation versus solid-state transformation mechanisms for amorphous to crystalline phase transitions. Specifically, we observe that in environments with abundant water, ACC transforms to calcite through a dissolution-reprecipitation mechanism. This features an activation energy of 63 ± 2 kJ/mol and unit cell volume contraction during calcite crystal growth. Conversely, under water-limited conditions, ACC to calcite transformation proceeds via a solid-state transformation mechanism, with an activation energy of 210 ± 2 kJ/mol, three times greater than the dissolution-reprecipitation route, and a unit cell expansion during crystalline calcite growth. Further, to illustrate the magnitude of these effects, the rates of calcite growth were similar during dissolution-reprecipitation at 3 °C [0.00207(35) s –1 ] and solid-state transformation at 280 °C [0.00134(11) s –1 ]. Moreover, the incorporation of an impurity, strontium, significantly retards the rate of calcite growth while expanding its unit cell but whose incorporation is history dependent. Reactive transport modeling of the dissolution–precipitation kinetics suggests that ACC must be dissolving as compact aggregates. These various transformation mechanisms drive diverse geological and biological carbonate formations, impacting their use as paleoenvironmental markers and functional materials synthesis.

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

In-Situ Measurement of Hall Thruster Erosion Using a Fiber Optic Regression Probe

One potential life-limiting mechanism in a Hall thruster is the erosion of the ceramic material comprising the discharge channel. This is especially true for missions that require long thrusting periods and can be problematic for lifetime qualification, especially when attempting to qualify a thruster by analysis rather than a test lasting the full duration of the mission. In addition to lifetime, several analytical and numerical models include electrode erosion as a mechanism contributing to enhanced transport properties. However, there is still a great deal of dispute over the importance of erosion to transport in Hall thrusters. The capability to perform an in-situ measurement of discharge channel erosion is useful in addressing both the lifetime and transport concerns. An in-situ measurement would allow for real-time data regarding the erosion rates at different operating points, providing a quick method for empirically anchoring any analysis geared towards lifetime qualification. Erosion rate data over a thruster s operating envelope would also be useful in the modeling of the detailed physics inside the discharge chamber. There are many different sensors and techniques that have been employed to quantify discharge channel erosion in Hall thrusters. Snapshots of the wear pattern can be obtained at regular shutdown intervals using laser profilometry. Many non-intrusive techniques of varying complexity and sensitivity have been employed to detect the time-varying presence of erosion products in the thruster plume. These include the use quartz crystal microbalances, emission spectroscopy, laser induced flourescence, and cavity ring-down spectroscopy. While these techniques can provide a very accurate picture of the level of eroded material in the thruster plume, it is more difficult to use them to determine the location from which the material was eroded. Furthermore, none of the methods cited provide a true in-situ measure of erosion at the channel surface while the thruster is in operation (i.e. none yield a continuous channel erosion measurement). A recent fundamental sensor development effort has led to a novel regression, erosion, and ablation sensor technology (REAST). The REAST sensor allows for measurement of real-time surface erosion rates at a discrete surface location. The sensor was tested using a linear Hall thruster geometry (see Fig. 1), which served as a means of producing plasma erosion of a ceramic discharge chamber. The mass flow rate, discharge voltage, and applied magnetic field strength could be varied, allowing for erosion measurements over a broad thruster operating envelope. Results are presented demonstrating the ability of the REAST sensor to capture not only the insulator erosion rates but also changes in these rates as a function of the discharge parameters.

Polzin, Kurt↗
DOE OSTI2025

Translational diffusion in supercooled water at and near the glass transition temperature—136 K

The properties of amorphous solid water at and near the calorimetric glass transition temperature, T g , of 136 K have been debated for years. One hypothesis is that water turns into a “true” liquid at T g (i.e., it becomes ergodic) and exhibits all the characteristics of an ergodic liquid, including translational diffusion. A competing hypothesis is that only rotational motion becomes active at T g , while the “real” glass transition in water is at a considerably higher temperature. To address this dispute, we have investigated the diffusive mixing in nanoscale water films, with thicknesses up to ∼100 nm, using infrared (IR) spectroscopy. The experiments used films that were composed of at least 90% H 2 O with D 2 O making up the balance and were conducted under conditions where H/D exchange was essentially eliminated. Because the IR spectra of multilayer D 2 O films (e.g., thicknesses of ∼3–6 nm) embedded within thick H 2 O films are distinct from the spectrum of isolated D 2 O molecules within H 2 O, the diffusive mixing of (initially) isotopically layered water films could be followed as a function of annealing time and temperature. The results show that water films with total thicknesses ranging from ∼20 to 100 nm diffusively mixed prior to crystallization for temperatures between 120 and 144 K. The translational diffusion had an Arrhenius temperature dependence with an activation energy of 40.8 ± 3.5 kJ/mol, which indicates that water at and near T g is a strong liquid. The measured diffusion coefficient at 136 K is 6.25 ± 1.4 × 10 −21 m 2 /s.

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

Dry Printing Pure Copper with High Conductivity and Adhesion for Flexible Electronics

Additive manufacturing of functional devices on various rigid and flexible substrates is rising rapidly due to their design flexibility, rapid manufacturing, and lower cost. Current printing technologies are ink-based and focused on printing silver (Ag) as conductive lines due to its matured ink formulation process, low sintering temperature, ease of printing, and low oxidation rate. However, Ag is the 68th most abundant element on Earth, while copper (Cu) is the 25th, making it much cheaper (>100×) while having a comparable conductivity to Ag. Therefore, printing Cu has become technologically and economically more attractive than Ag. Nevertheless, Cu printing is still a significant challenge in ink-based printing methods due to the higher sintering temperature relative to the glass-transition temperature of most flexible substrates, the higher oxidation rate, the challenging ink formulation process, and ink stability concerns. Here, we demonstrate printing highly conductive Cu on flexible polyimide substrates using a dry printing technique. Cu nanoparticles (~3–30 nm) are generated by on-demand laser ablation of a solid Cu target inside the printer head and under argon background gas. These Cu nanoparticles are then transported through a nozzle and onto the substrate, where they are laser-sintered in real time. The argon gas plays three critical roles in laser plume condensation for nanoparticle generation, transport, and sheath gas to avoid oxidation during sintering. The sintered nanoparticles thus show high electrical conductivity and mechanical stability under static and cyclic tests. Our dry printing technique can potentially revolutionize how electronic devices and sensors are additively manufactured for earth and space applications.

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

Combining artificial intelligence and physics-based modeling to directly assess atomic site stabilities: from sub-nanometer clusters to extended surfaces

The performance of functional materials is dictated by chemical and structural properties of individual atomic sites. In catalysts, for instance, the thermodynamic stability of constituting atomic sites is a key descriptor from which more complex properties, such as molecular adsorption energies and reaction rates, can be derived. In this study, we present a widely applicable machine learning (ML) approach to instantaneously compute the stability of individual atomic sites in structurally and electronically complex nano-materials. Conventionally, we determine such site stabilities using computationally intensive first-principles calculations. With our approach, we predict the stability of atomic sites in sub-nanometer metal clusters of 3–55 atoms with mean absolute errors in the range of 0.11–0.14 eV. To extract physical insights from the ML model, we introduce a genetic algorithm (GA) for feature selection. This algorithm distills the key structural and chemical properties governing the stability of atomic sites in size-selected nanoparticles, allowing for physical interpretability of the models and revealing structure–property relationships. The results of the GA are generally model and materials specific. In the limit of large nanoparticles, the GA identifies features consistent with physics-based models for metal–metal interactions. By combining the ML model with the physics-based model, we predict atomic site stabilities in real time for structures ranging from sub-nanometer metal clusters (3–55 atom) to larger nanoparticles (147 to 309 atoms) to extended surfaces using a physically interpretable framework. Finally, we present a proof of principle showcasing how our approach can determine stable and active nanocatalysts across a generic materials space of structure and composition.

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

pyRMG: A framework for high-throughput, large-cell DFT calculations on supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

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

Electron dynamics in complex environments with real-time time dependent density functional theory in a QM-MM framework

This article presents a time dependent density functional theory (TDDFT) implementation to propagate the Kohn-Sham equations in real time, including the effects of a molecular environment through a Quantum-Mechanics Molecular-Mechanics (QM-MM) hamiltonian. The code delivers an all-electron description employing Gaussian basis functions, and incorporates the Amber force-field in the QM-MM treatment. The most expensive parts of the computation, comprising the commutators between the hamiltonian and the density matrix—required to propagate the electron dynamics—, and the evaluation of the exchange-correlation energy, were migrated to the CUDA platform to run on graphics processing units, which remarkably accelerates the performance of the code. The method was validated by reproducing linear-response TDDFT results for the absorption spectra of several molecular species. Two different schemes were tested to propagate the quantum dynamics: (i) a leap-frog Verlet algorithm, and (ii) the Magnus expansion to first-order. These two approaches were confronted, to find that the Magnus scheme is more efficient by a factor of six in small molecules. Interestingly, the presence of iron was found to seriously limitate the length of the integration time step, due to the high frequencies associated with the core-electrons. This highlights the importance of pseudopotentials to alleviate the cost of the propagation of the inner states when heavy nuclei are present. Finally, the methodology was applied to investigate the shifts induced by the chemical environment on the most intense UV absorption bands of two model systems of general relevance: the formamide molecule in water solution, and the carboxy-heme group in Flavohemoglobin. In both cases, shifts of several nanometers are observed, consistently with the available experimental data.

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

Fundamental Impact of Humidity on SOFC Cathode ORR

Although solid oxide fuel cells (SOFC) have demonstrated excellent performance, the durability of SOFCs under real working conditions is still an issue for commercial deployment. In particular cathode exposure to atmospheric air contaminants, such as humidity, can result in long-term performance degradation issues. Therefore, a fundamental understanding of the interaction between water molecules and cathodes is essential to resolve this issue and further enhance cathode durability. In order to study the effects of humidity on the oxygen reduction reaction (ORR), we used in-situ 18 O isotope exchange techniques to probe the exchange of water with two of themost common SOFC cathode materials, (La 0.8 Sr 0.2 ) 0.95 MnO 3±δ (LSM) and La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF). In this experiment, heavy water, D 2 O (with a mass/charge ratio of m/z = 20), is used to avoid the overlapping of H 2 O and the 18 O 2 cracking fraction, which both provide a peak at m/z = 18. A series of temperature programmed isotope exchange measurements were performed to comprehensively study the interaction of water with the cathode surface as a function of temperature, oxygen partial pressure, and water vapor concentration. The results suggest that water and O 2 share the same surface exchange sites, leading to competitive adsorption. Our findings show that water prefers to exchange with LSCF at lower temperatures, around 300–450°C. For LSM, O 2 is more favorable than water to be adsorbed on the surface and the presence of O 2 limits water exchange. The experimental data are summarized in a Temperature-PO 2 diagram to help visualize how the exchange of water on each material depends on the operating conditions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Surface-Functionalized, Two-Dimensional Polymer Electrochromic Layers as Ultrafast, Multi-State Infrared Optical Gates

Electrochromic devices have widespread application potential, but the currently available switching speeds limit broad real-world implementation of this technology. Here, we report surface-engineered two-dimensional polymers with ionophilic pores that offer unprecedented switching speeds in solid-state, two-terminal, electrochromic devices. In particular, we demonstrate that a crystalline donor–acceptor 2DP functionalized with ethylene glycol oligomers exhibits multistate infrared absorption that is 4× faster (t c = 320 ms) with 3× coloration efficiency (491 cm 2 C –1 ) compared to an alkyl functionalized 2DP constructed from the same chromophores. The functionalized nanoporous surfaces enable rapid switching in these materials under either oxidative or reductive conditions, allowing us to access a range of robust, stable optical responses in a single electrochromic layer. These attributes led us to leverage surface-functionalized 2DPs as multistate infrared logic gates. Collectively, this work demonstrates that surface engineering of nanoporous crystalline lattices is a promising approach to co-optimize the electronic and ionic conductivities required to achieve rapidly switchable electrochromic layers. Beyond speed and efficiency, the demonstration of multistate infrared characteristics shows that electrochromic frameworks are useful in integrated optoelectronic circuits. This positions surface-engineered 2DPs as improved electrochromic coatings and a new material platform for photonic information processing and adaptive devices.

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

Adaptive Augmenting Control Flight Characterization Experiment on an F/A-18

This paper summarizes the Adaptive Augmenting Control (AAC) flight characterization experiments performed using an F/A-18 (TN 853). AAC was designed and developed specifically for launch vehicles, and is currently part of the baseline autopilot design for NASA's Space Launch System (SLS). The scope covered here includes a brief overview of the algorithm (covered in more detail elsewhere), motivation and benefits of flight testing, top-level SLS flight test objectives, applicability of the F/A-18 as a platform for testing a launch vehicle control design, test cases designed to fully vet the AAC algorithm, flight test results, and conclusions regarding the functionality of AAC. The AAC algorithm developed at Marshall Space Flight Center is a forward loop gain multiplicative adaptive algorithm that modifies the total attitude control system gain in response to sensed model errors or undesirable parasitic mode resonances. The AAC algorithm provides the capability to improve or decrease performance by balancing attitude tracking with the mitigation of parasitic dynamics, such as control-structure interaction or servo-actuator limit cycles. In the case of the latter, if unmodeled or mismodeled parasitic dynamics are present that would otherwise result in a closed-loop instability or near instability, the adaptive controller decreases the total loop gain to reduce the interaction between these dynamics and the controller. This is in contrast to traditional adaptive control logic, which focuses on improving performance by increasing gain. The computationally simple AAC attitude control algorithm has stability properties that are reconcilable in the context of classical frequency-domain criteria (i.e., gain and phase margin). The algorithm assumes that the baseline attitude control design is well-tuned for a nominal trajectory and is designed to adapt only when necessary. Furthermore, the adaptation is attracted to the nominal design and adapts only on an as-needed basis (see Figure 1). The MSFC algorithm design was formulated during the Constellation Program and reached a high maturity level during SLS through simulation-based development and internal and external analytical review. The AAC algorithm design has three summary-level objectives: (1) "Do no harm;" return to baseline control design when not needed, (2) Increase performance; respond to error in ability of vehicle to track command, and (3) Regain stability; respond to undesirable control-structure interaction or other parasitic dynamics. AAC has been successfully implemented as part of the Space Launch System baseline design, including extensive testing in high-fidelity 6-DOF simulations the details of which are described in [1]. The Dryden Flight Research Center's F/A-18 Full-Scale Advanced Systems Testbed (FAST) platform is used to conduct an algorithm flight characterization experiment intended to fully vet the aforementioned design objectives. FAST was specifically designed with this type of test program in mind. The onboard flight control system has full-authority experiment control of ten aerodynamic effectors and two throttles. It has production and research sensor inputs and pilot engage/disengage and real-time configuration of up to eight different experiments on a single flight. It has failure detection and automatic reversion to fail-safe mode. The F/A-18 aircraft has an experiment envelope cleared for full-authority control and maneuvering and exhibits characteristics for robust recovery from unusual attitudes and configurations aided by the presence of a qualified test pilot. The F/A-18 aircraft has relatively high mass and inertia with exceptional performance; the F/A-18 also has a large thrust-to-weight ratio, owing to its military heritage. This enables the simulation of a portion of the ascent trajectory with a high degree of dynamic similarity to a launch vehicle, and the research flight control system can simulate unstable longitudinal dynamics. Parasitic dynamics such as slosh and bending modes, as well as atmospheric disturbances, are being produced by the airframe via modification of bending filters and the use of secondary control surfaces, including leading and trailing edge flaps, symmetric ailerons, and symmetric rudders. The platform also has the ability to inject signals in flight to simulate structural mode resonances or other challenging dynamics. This platform also offers more test maneuvers and longer maneuver times than a single rocket or missile test, which provides ample opportunity to fully and repeatedly exercise all aspects of the algorithm. Prior to testing on an F/A-18, AAC was the only component of the SLS autopilot design that had not been flight tested. The testing described in this paper raises the Technology Readiness Level (TRL) early in the SLS Program and is able to demonstrate its capabilities and robustness in a flight environment.

VanZwieten, Tannen S.↗
DOE OSTI2021

Connection between the physicochemical characteristics of amorphous carbon thin films and their electrochemical properties

Connecting a material's surface chemistry with its electrocatalytic performance is one of the major questions in analytical electrochemistry. This is especially important in many sensor applications where analytes from complex media need to be measured. Unfortunately, today this connection is still largely missing except perhaps for the most simple ideal model systems. Here we present an approach that can be used to obtain insights about this missing connection and apply it to the case of carbon nanomaterials. In this paper we show that by combining advanced computational techniques augmented by machine learning methods with x-ray absorption spectroscopy (XAS) and electrochemical measurements, it is possible to obtain a deeper understanding of the correlation between local surface chemistry and electrochemical performance. As a test case we show how by computationally assessing the growth of amorphous carbon (a-C) thin films at the atomic level, we can create computational structural motifs that may in turn be used to deconvolute the XAS data from the real samples resulting in local chemical information. Then, by carrying out electrochemical measurements on the same samples from which x-ray spectra were measured and that were further characterized computationally, it is possible to gain insight into the interplay between the local surface chemistry and electrochemical performance. To demonstrate this methodology, we proceed as follows: after assessing the basic electrochemical properties of a-C films, we investigate the effect of short HNO 3 treatment on the sensitivity of these electrodes towards an inner sphere redox probe dopamine to gain knowledge about the influence of altered surface chemistry to observed electrochemical performance. Overall, these results pave the way towards a more general assessment of electrocatalysis in different systems and provide the first steps towards data driven tailoring of electrode surfaces to gain optimal performance in a given application.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

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

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

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

Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT

We present a special symmetric Lanczos algorithm and a kernel polynomial method (KPM) for approximating the absorption spectrum of molecules within the linear response time-dependent density functional theory (TDDFT) framework in the product form. In contrast to existing algorithms, the new algorithms are based on reformulating the original non-Hermitian eigenvalue problem as a product eigenvalue problem and the observation that the product eigenvalue problem is self-adjoint with respect to an appropriately chosen inner product. This allows a simple symmetric Lanczos algorithm to be used to compute the desired absorption spectrum. The use of a symmetric Lanczos algorithm only requires half of the memory compared with the nonsymmetric variant of the Lanczos algorithm. The symmetric Lanczos algorithm is also numerically more stable than the nonsymmetric version. The KPM algorithm is also presented as a low-memory alternative to the Lanczos approach, but the algorithm may require more matrix-vector multiplications in practice. We discuss the pros and cons of these methods in terms of their accuracy as well as their computational and storage cost. Applications to a set of small and medium-sized molecules are also presented.

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

Size-segregated emissions and metal content of vehicle-emitted particles as a function of mileage: Implications to population exposure

The study aims at investigating the characteristics (size distribution, active surface and metal content) of particles emitted by cars as a function of mileage using a novel methodology for characterizing particulate emissions captured by Exhaust Gas Suspension (EGS). EGS was obtained by passing the exhaust gases through a container of deionized water. EGS analysis was performed using laser granulometry, electron scanning microscopy, and high resolution mass spectrometry. Implications of the differences in key features of the emitted particles on population exposure were investigated using numerical simulation for estimating size-segregated PM deposition across human respiratory tract (HRT). It was found that vehicle mileage, age and the respective emissions class have almost no effect on the size distribution of the exhaust gas particulate released into the environment; about half of the examined vehicles with low mileage were found to release particles of aerodynamic diameter above 10 μm. The exhaust gas particulate detected in the EGS of all cars can be classified into three major size classes: (1) 0.1–5 µm – soot and ash particles, metals (Au, Pt, Pd, Ir); (2) 10–30 µm – metal (Cr, Fe, Cu, Zr, Ni) and ash particles; (3) 400–1,000 µm – metal (Fe, Cr, Pb) and ash particles. Newer vehicles with low mileage are substantial sources of soot and metal particles with median diameter of 200 nm with a higher surface area (up to 89,871.16 cm{sup 2}/cm{sup 3}). These tend to deposit in the lower part of the human respiratory tract. - Highlights: • Car mileage has virtually no effect on the size of the solid particles released. • Newer diesel vehicles emit particles of lower aerodynamic diameter. • Particle active surface emitted by newer vehicles is on average 3 times higher. • Real-life emissions were translated into actual internal PM exposure.

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

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

Stretched or noded orbital densities and self-interaction correction in density functional theory

Semilocal approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semilocal approximation makes that approximation exact for all one-electron ground- or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semilocal approximations satisfy, suggesting the need for a generalized SIC. These conclusions are supported by calculations for one-electron densities and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGAs) and meta-GGAs, it reduces and often worsens the atomization energies of molecules. Thus, PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here, in particular, to the strongly constrained and appropriately normed (SCAN) meta-GGA, for which the correlation part is already self-interaction-free. This property makes SCAN a natural first candidate for a generalized SIC.

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

Quantifying Temperature Dependence of Pu(IV) Absorbance Spectra for Advanced Online Monitoring of Nuclear Processes

This article presents a systematic study of Pu(IV) absorbance spectral features as a function of temperature to develop an understanding of this parameter’s effect on chemometric models that can be used as online monitoring tools to support nuclear processing. The descriptive and predictive models that provide real-time feedback of these processes are usually constructed with data collected in conditions typical of a laboratory environment, which can differ drastically from a processing environment. To assess the impact of temperature on Pu(IV) absorbance spectra, 11 samples of Pu(IV) were synthesized with varying HNO 3 concentrations ranging from 0.6 to 9.5 M and heated between 15 and 45 °C. Ultraviolet (UV)–visible (vis)–near-infrared (NIR) absorption spectra collected at different HNO 3 concentrations and temperatures revealed that features associated with Pu(IV) are sensitive to temperature at all HNO 3 concentrations and that changes in features depend on HNO 3 concentration. The contributions of temperature and HNO 3 concentration to variation in Pu(IV) spectral features were evaluated using the principal component analysis of spectra that were baseline-corrected with an asymmetric least-squares method. Furthermore, predictive modeling for HNO 3 concentration with partial least-squares regression of UV–vis–NIR spectra highlighted the importance of accounting for temperature in the calibration set to optimize model performance. This methodology constitutes a new, systematic approach to account for the effect of temperature on the absorption spectra of metal ions and is useful for process monitoring applications in many industries.

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

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗