Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function approximation”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2015

Relativistic Coulomb excitation within the time dependent superfluid local density approximation

Within the framework of the unrestricted time-dependent density functional theory, we present for the first time an analysis of the relativistic Coulomb excitation of the heavy deformed open shell nucleus 238 U. The approach is based on the superfluid local density approximation formulated on a spatial lattice that can take into account coupling to the continuum, enabling self-consistent studies of superfluid dynamics of any nuclear shape. We compute the energy deposited in the target nucleus as a function of the impact parameter, finding it to be significantly larger than the estimate using the Goldhaber-Teller model. The isovector giant dipole resonance, the dipole pygmy resonance, and giant quadrupole modes are excited during the process. As a result, the one-body dissipation of collective dipole modes is shown to lead a damping width Γ↓≈0.4 MeV and the number of preequilibrium neutrons emitted has been quantified.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Phase‐Dependent Band Gap Engineering in Alloys of Metal‐Semiconductor Transition Metal Dichalcogenides

Abstract Bandgap engineering plays a critical role in optimizing the electrical, optical and (photo)‐electrochemical applications of semiconductors. Alloying has been a historically successful way of tuning bandgaps by making solid solutions of two isovalent semiconductors. In this work, a novel form of bandgap engineering involving alloying non‐isovalent cations in a 2D transition metal dichalcogenide (TMDC) is presented. By alloying semiconducting MoSe 2 with metallic NbSe 2 , two structural phases of Mo 0.5 Nb 0.5 Se 2 , the 1T and 2H phases, are produced each with emergent electronic structure. At room temperature, it is observed that the 1T and 2H phases are semiconducting and metallic, respectively. For the 1T structure, scanning tunneling microscopy/spectroscopy (STM/STS) is used to measure band gaps in the range of 0.42–0.58 at 77 K. Electron diffraction patterns of the 1T structure obtained at room temperature show the presence of a nearly commensurate charge density wave (NCCDW) phase with periodic lattice distortions that result in an uncommon 4 × 4 supercell, rotated approximately 4° from the lattice. Density‐functional‐theory calculations confirm that local distortions, such as those in a NCCDW, can open up a band gap in 1T ‐Mo 0.5 Nb 0.5 Se 2 , but not in the 2H phase. This work expands the boundaries of alloy‐based bandgap engineering by introducing a novel technique that facilitates CDW phases through alloying.

Wang, Shuxi↗
DOE OSTI2015

Comparison of multigrid algorithms for high‐order continuous finite element discretizations

Summary We present a comparison of different multigrid approaches for the solution of systems arising from high‐order continuous finite element discretizations of elliptic partial differential equations on complex geometries. We consider the pointwise Jacobi, the Chebyshev‐accelerated Jacobi, and the symmetric successive over‐relaxation smoothers, as well as elementwise block Jacobi smoothing. Three approaches for the multigrid hierarchy are compared: (1) high‐order h ‐multigrid, which uses high‐order interpolation and restriction between geometrically coarsened meshes; (2) p ‐multigrid, in which the polynomial order is reduced while the mesh remains unchanged, and the interpolation and restriction incorporate the different‐order basis functions; and (3) a first‐order approximation multigrid preconditioner constructed using the nodes of the high‐order discretization. This latter approach is often combined with algebraic multigrid for the low‐order operator and is attractive for high‐order discretizations on unstructured meshes, where geometric coarsening is difficult. Based on a simple performance model, we compare the computational cost of the different approaches. Using scalar test problems in two and three dimensions with constant and varying coefficients, we compare the performance of the different multigrid approaches for polynomial orders up to 16. Overall, both h ‐multigrid and p ‐multigrid work well; the first‐order approximation is less efficient. For constant coefficients, all smoothers work well. For variable coefficients, Chebyshev and symmetric successive over‐relaxation smoothing outperform Jacobi smoothing. While all of the tested methods converge in a mesh‐independent number of iterations, none of them behaves completely independent of the polynomial order. When multigrid is used as a preconditioner in a Krylov method, the iteration number decreases significantly compared with using multigrid as a solver. Copyright © 2015 John Wiley & Sons, Ltd.

Sundar, Hari↗
DOE OSTI2021

Subleading corrections to the S 3 free energy of necklace quiver theories dual to massive IIA

We investigate the S 3 free energy of N = 3 Chern-Simons-matter quiver gauge theories with gauge group U(N) r (r ≥ 2) where the sum of Chern-Simons levels does not vanish, beyond the leading order in the large-N limit. We take two different approaches to explore the sub-leading structures of the free energy. First we evaluate the matrix integral for the partition function in the ’t Hooft limit using a saddle point approximation. Second we use an ideal Fermi-gas model to compute the same partition function, but in the limit of fixed Chern-Simons levels. The resulting expressions for the free energy F = – log Z are then compared in the overlapping parameter regime. The Fermi-gas approach also hints at a universal 1/6 log N correction to the free energy. Since the quiver gauge theories we consider are dual to massive Type IIA theory, we expect the sub-leading correction of the planar free energy in the large ’t Hooft parameter limit to match higher-derivative corrections to the tree-level holographic dual free energy, which have not yet been fully investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Outcomes of Temporary Partially Covered Stent Placement for Benign Tracheobronchial Stenosis

PurposeTo evaluate the intermediate outcomes of temporary partially covered tracheobronchial stenting in patients with benign tracheobronchial stenosis.Materials and MethodsWe conducted a retrospective study of patients with benign tracheobronchial stenosis who underwent stent placement. All stents were removed approximately 3 months after placement. Respiratory function was assessed using the visual analogue scale (VAS) and Karnofsky performance status scale (KPS) before and after stent placement. The lumen diameters of the stenotic lesions were measured using chest computed tomography (CT) and compared between before stent placement and after stent removal.ResultsA total of 51 stents were successfully placed in consecutive 51 patients with benign tracheobronchial stenosis. No serious complications occurred. The mean VAS and KPS scores significantly improved after stent removal (6.291 ± 0.495 and 25.352 ± 10.533, respectively) compared with those before stent placement (1.493 ± 0.504 and 60.140 ± 16.344, respectively; P < 0.05). The mean lumen transverse diameters of the stenotic site in trachea and main bronchus after stent removal (17.235 ± 3.457 and 8.993 ± 0.961 mm; 1 month post-removal; 16.353 ± 4.132 and 8.357 ± 1.082 mm; 6 months post-removal) were significantly larger than those before stent placement (7.876 ± 2.351 and 2.143 ± 0.770 mm, respectively; P < 0.05). However, the mean lumen diameters between 1 and 6 months after stent removal had no significant difference (P > 0.05).ConclusionTemporary partially covered stenting may be a safe and effective treatment for benign tracheobronchial stenosis.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Stochastic derivative-free optimization using a trust region framework

This study presents a trust region algorithm to minimize a function f when one has access only to noise-corrupted function values f¯. The model-based algorithm dynamically adjusts its step length, taking larger steps when the model and function agree and smaller steps when the model is less accurate. The method does not require the user to specify a fixed pattern of points used to build local models and does not repeatedly sample points. If f is sufficiently smooth and the noise is independent and identically distributed with mean zero and finite variance, we prove that our algorithm produces iterates such that the corresponding function gradients converge in probability to zero. As a result, we present a prototype of our algorithm that, while simplistic in its management of previously evaluated points, solves benchmark problems in fewer function evaluations than do existing stochastic approximation methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

GGA+U method from first principles: application to reduction–oxidation properties in ceria-based oxides

We implement DFT calculations by a Hubbard-like correction for localized strongly correlated electrons, coupled with a generalized gradient approximation to the exchange-correlation functional to elucidate the role of the ceria based catalytically active supports for the chemical reactions involving reduction–oxidation processes. These catalytic processes are relevant for many industrial applications, such as catalytic converters in automotive applications, solid oxide fuel cells, and hydrogen production from biomass. The Hubbard-like correction U is computed from first principles as physical property of the system. We find that the high performance of ceria-based oxides as an active support for noble metals in catalysis relies on an efficient supply of lattice oxygen at reaction sites governed by oxygen vacancy formation.

Materials Science↗
DOE OSTI2019

Effect of V{sub Zn}/V{sub O} on Stability, Magnetism, and Electronic Characteristic of Oxygen Ions for Li-Doped ZnO

The effects of Zn/O vacancy (V{sub Zn}/V{sub O}) and different proportions of Li/V{sub Zn} on the magnetism of Li-doped ZnO are analyzed through first-principle calculation by using generalized gradient approximation+ U (GGA + U) under density functional theory. Results reveal that Li-doped ZnO with V{sub Zn} can realize ferromagnetic long-range order and has a Curie temperature above room temperature. Under different proportions of Li/V{sub Zn} (i.e., 1:1, 1:2, and 2:2) in ZnO (2 × 2 × 4), the doping system containing 2Li/2V{sub Zn} (Zn{sub 28}Li{sub 2}O{sub 32}) shows the greatest magnetic moment and the smallest differential charge density. These characteristics pave the way for enhancing the magnetic properties of dilute magnetic semiconductors. The oxygen atoms in Zn{sub 28}Li{sub 2}O{sub 32} show acceptor and donor characteristics and exist in the forms of itinerant electrons (O{sup 1−}) and local electrons (O{sup 2−}), which have different effects on Zn{sub 28}Li{sub 2}O{sub 32} magnetism. The spin-polarization double-exchange effect among the unpaired itinerant electron (O{sup 1−}) orbit, local electron (O{sup 2−}) orbit, and unpaired Zn-3d electron orbit is the origin of magnetism for Li-doped ZnO with V{sub Zn}. By contrast, the doping systems of Li-doped ZnO with V{sub O} are nonmagnetic, rendering such systems inapplicable.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Optimal carbon storage reservoir management through deep reinforcement learning

Model-based optimization plays a central role in energy system design and management. The complexity and high-dimensionality of many process-level models, especially those used for geosystem energy exploration and utilization, often lead to formidable computational costs when the dimension of decision space is also large. This work adopts elements of recently advanced deep learning techniques to solve a sequential decision-making problem in applied geosystem management. Specifically, a deep reinforcement learning framework was formed for optimal multiperiod planning, in which a deep Q-learning network (DQN) agent was trained to maximize rewards by learning from high-dimensional inputs and from exploitation of its past experiences. To expedite computation, deep multitask learning was used to approximate high-dimensional, multistate transition functions. Both DQN and deep multitask learning are pattern based. As a demonstration, the framework was applied to optimal carbon sequestration reservoir planning using two different types of management strategies: monitoring only and brine extraction. Both strategies are designed to mitigate potential risks due to pressure buildup. Results show that the DQN agent can identify the optimal policies to maximize the reward for given risk and cost constraints. Finally, experiments also show that knowledge the agent gained from interacting with one environment is largely preserved when deploying the same agent in other similar environments.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Deep neural operators can predict the real-time response of floating offshore structures under irregular waves

The use of neural operators in a digital twin model of an offshore floating structure holds the potential for a significant shift in the prediction of structural responses and health monitoring, offering valuable real-time control insights. In this work, we investigate the effectiveness of three neural operators, namely the deep operator network (DeepONet), the Fourier neural operator (FNO), and the Wavelet neural operator (WNO), to accurately capture the responses of a floating structure under six different sea state codes (3 − 8) based on the wave characteristics described by the World Meteorological Organization (WMO). To further enhance the accuracy of the vanilla architecture of the neural operators, novel extensions, such as wavelet-DeepONet and self-adaptive WNO, are proposed in this paper. The results demonstrate that these high-precision neural operators can deliver structural responses more efficiently, up to two orders of magnitude faster than a dynamic analysis using conventional numerical solvers. Additionally, compared to gated recurrent units (GRUs), a commonly used recurrent neural network for time-series estimation, neural operators are both more accurate and efficient, especially in situations with limited data availability. Taken together, our study shows that FNO outperforms all other operators for approximating the mapping of one input functional space to the output space as well as for responses that have small bandwidth of the frequency spectrum. Conversely, DeepONet, with historical states, proves most accurate in learning the mapping of multiple input functions to the output space and capturing responses within a broad frequency spectrum.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A systematic approach to resolve high impedance of multifunctional energy storage composites

Previous work has proposed and characterized the structural and electrical performance of Multifunctional Energy Storage Composite (MESC) structures: structural elements with embedded lithium-ion batteries which were developed by the Structures and Composite Laboratory (SACL) at Stanford University. This work conducts a comprehensive parametric study over the design variables in the manufacturing process to determine the root cause of increased DC impedance compared to commercial pouch cell batteries. Several MESC samples were assembled under various controlled environmental temperatures and pressures. The samples underwent Hybrid Peak Pulse Characterization (HPPC) to approximate the DC impedance as a function of state of charge (SOC). Failure limits of the battery with respect to assembly temperature and pressure were tested, and an MESC manufacturing simulation model was constructed in Abaqus to validate the manufacturing limits. Samples assembled without enough temperature or pressure failed due to a lack of structural integrity or hermetic seal. Conversely, samples which underwent elevated temperatures and pressures resulted in high DC impedance due to nonuniform electrolyte wetting. This work resolves the primary concern of MESC by achieving comparable DC impedance to commercial pouch cells of the same chemistry and testing procedure.

25 ENERGY STORAGE↗
DOE OSTI2014

Neurochemical Evidence of Potential Neurotoxicity After Prophylactic Cranial Irradiation

Purpose: To examine whether cerebrospinal fluid biomarkers for neuroaxonal damage, neuroglial activation, and amyloid β–related processes could characterize the neurochemical response to cranial radiation. Methods and Materials: Before prophylactic cranial irradiation (PCI) of patients with small cell lung cancer, each patient underwent magnetic resonance imaging of the brain, lumbar puncture, and Mini-Mental State Examination of cognitive function. These examinations were repeated at approximately 3 and 12 months after radiation. Results: The major findings were as follows. (1) Cerebrospinal fluid markers for neuronal and neuroglial injury were elevated during the subacute phase after PCI. Neurofilament and T-tau increased 120% and 50%, respectively, after PCI (P<.05). The same was seen for the neuroglial markers YKL-40 and glial fibrillary acidic protein, which increased 144% and 106%, respectively, after PCI (P<.05). (2) The levels of secreted amyloid precursor protein-α and -β were reduced 44% and 46%, respectively, 3 months after PCI, and the levels continued to decrease as long as 1 year after treatment (P<.05). (3) Mini-Mental State Examination did not reveal any cognitive decline, indicating that a more sensitive test should be used in future studies. Conclusion: In conclusion, we were able to detect radiation therapy–induced changes in several markers reflecting neuronal injury, inflammatory/astroglial activation, and altered amyloid precursor protein/amyloid β metabolism, despite the low number of patients and quite moderate radiation doses (20-30 Gy). These changes are hypothesis generating and could potentially be used to assess the individual risk of developing long-term symptoms of chronic encephalopathy after PCI. This has to be evaluated in large studies with extended clinical follow-up and more detailed neurocognitive assessments.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Adjoint DSMC for nonlinear Boltzmann equation constrained optimization

Applications for kinetic equations such as optimal design and inverse problems often involve finding unknown parameters through gradient-based optimization algorithms. Based on the adjoint-state method, we derive two different frameworks for approximating the gradient of an objective functional constrained by the nonlinear Boltzmann equation. While the forward problem can be solved by the DSMC method, it is difficult to efficiently solve the high-dimensional continuous adjoint equation obtained by the “optimize-then-discretize” approach. This challenge motivates us to propose an adjoint DSMC method following the “discretize-then-optimize” approach for Boltzmann-constrained optimization. We also analyze the properties of the two frameworks and their connections. Here, several numerical examples are presented to demonstrate their accuracy and efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Adversarial sampling of unknown and high-dimensional conditional distributions

Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximations need to be made when the distribution is unknown. In this paper the sampling method, as well as the inference of the underlying distribution, are both handled with a data-driven method known as generative adversarial networks (GAN), which trains two competing neural networks to produce a network that can effectively generate samples from the training set distribution. In practice, it is often necessary to draw samples from conditional distributions. When the conditional variables are continuous, only one (if any) data point corresponding to a particular value of a conditioning variable may be available, which is not sufficient to estimate the conditional distribution. This work handles this problem using an a priori estimation of the conditional moments of a PDF. Herein, two approaches, stochastic estimation, and an external neural network are compared for computing these moments; however, any preferred method can be used. The algorithm is demonstrated in the case of the deconvolution of a filtered turbulent flow field. It is shown that all the versions of the proposed algorithm effectively sample the target conditional distribution with minimal impact on the quality of the samples compared to state-of-the-art methods. Additionally, the procedure can be used as a metric for the diversity of samples generated by a conditional GAN (cGAN) conditioned with continuous variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

First-Principles Search for Efficient Activators for LaI3

First-principles calculations were performed using density functional theory with Hubbard corrections or hybrid exchange-correlation functionals, as well as the GW approximation, to identify activator dopants that could maximize the light yield of LaI3, a potentially very bright scintillator. The dopants considered in this work included a series of lanthanide ions (Ce, Pr, Nd, Eu, Gd, and Tb) and several ns2 ions (Tl, Pb, Sb, and Bi). Based on both ground-state calculations and constrained DFT calculations to simulate excited states, the trivalent lanthanide dopants were shown not to constitute an improvement over Ce3+, for which experimental data exist, as they showed occupied 4f states below the valence band maximum (VBM) and 5d states above the conduction band maximum (CBM). In contrast, the only divalent lanthanide considered, Eu2+, displayed occupied 4f states within the band gap of the host, but its 5d states were calculated to be above the CBM. Eu2+ could nonetheless be exploited as an activator by increasing the band gap energy slightly through substitution of iodide ions by bromide ions. Similar results were obtained for Tl+. Finally, Bi3+ and Sb3+ were predicted to be efficient activators in LaI3 by virtue of having unoccupied p states below the CBM, which can serve as electron traps and can combine with holes at the VBM to form localized emission centers. Comparison with empirical relationships of the relative positions of the energy levels of rare-earth dopants and implications for efficient doping schemes and scintillation mechanisms in LaI3 are also discussed.

Wu, Dangxin↗