Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Optical Performance Modeling of FUSE Telescope Mirror

We describe the Metrology Data Processor (METDAT), the Optical Surface Analysis Code (OSAC), and their application to the image evaluation of the Far Ultraviolet Spectroscopic Explorer (FUSE) mirrors. The FUSE instrument - designed and developed by the Johns Hopkins University and launched in June 1999 is an astrophysics satellite which provides high resolution spectra (lambda/Delta(lambda) = 20,000 - 25,000) in the wavelength region from 90.5 to 118.7 nm The FUSE instrument is comprised of four co-aligned, normal incidence, off-axis parabolic mirrors, four Rowland circle spectrograph channels with holographic gratings, and delay line microchannel plate detectors. The OSAC code provides a comprehensive analysis of optical system performance, including the effects of optical surface misalignments, low spatial frequency deformations described by discrete polynomial terms, mid- and high-spatial frequency deformations (surface roughness), and diffraction due to the finite size of the aperture. Both normal incidence (traditionally infrared, visible, and near ultraviolet mirror systems) and grazing incidence (x-ray mirror systems) systems can be analyzed. The code also properly accounts for reflectance losses on the mirror surfaces. Low frequency surface errors are described in OSAC by using Zernike polynomials for normal incidence mirrors and Legendre-Fourier polynomials for grazing incidence mirrors. The scatter analysis of the mirror is based on scalar scatter theory. The program accepts simple autocovariance (ACV) function models or power spectral density (PSD) models derived from mirror surface metrology data as input to the scatter calculation. The end product of the program is a user-defined pixel array containing the system Point Spread Function (PSF). The METDAT routine is used in conjunction with the OSAC program. This code reads in laboratory metrology data in a normalized format. The code then fits the data using Zernike polynomials for normal incidence systems or Legendre-Fourier polynomials for grazing incidence systems. It removes low order terms from the metrology data, calculates statistical ACV or PSD functions, and fits these data to OSAC models for the scatter analysis. In this paper we briefly describe the laboratory image testing of FUSE spare mirror performed in the near and vacuum ultraviolet at John Hopkins University and OSAC modeling of the test setup performed at NASA/GSFC. The test setup is a double-pass configuration consisting of a Hg discharge source, the FUSE off-axis parabolic mirror under test, an autocollimating flat mirror, and a tomographic imaging detector. Two additional, small fold flats are used in the optical train to accommodate the light source and the detector. The modeling is based on Zernike fitting and PSD analysis of surface metrology data measured by both the mirror vendor (Tinsley) and JHU. The results of our models agree well with the laboratory imaging data, thus validating our theoretical model. Finally, we predict the imaging performance of FUSE mirrors in their flight configuration at far-ultraviolet wavelengths.

Saha, Timo T.↗
DOE OSTI2025

Interplay between electron localization, magnetic order, and Jahn-Teller distortion dictates LiMnO2 phase stability

The development of manganese (Mn)-rich cathodes for Li-ion batteries promises to alleviate potential supply chain bottlenecks in battery manufacturing. Fundamental challenges in Mn-rich cathodes arise from phenomena such as structural changes due to cooperative Jahn-Teller (JT) distortions of in octahedral environments, Mn migration, and phase transformations to spinel-like order, all of which affect the electrochemical performance. These physically complex phenomena motivate an re-examination of the Li-Mn-O rock-salt space, with a focus on the thermodynamics of the prototypical, polymorphs. It is found that the generalized gradient approximation (GGA-PBEsol) and meta-GGA ( ) density functionals with empirically fitted on-site Hubbard corrections yield spurious stable phases for , such as predicting a phase with -like order ( ) to be the ground state instead of the orthorhombic (Pmmn) phase, which is the experimentally known ground state. Accounting for antiferromagnetic order in each structure is shown to have a substantial effect on the total energies and resulting phase stability. By using hybrid-GGA (HSE06) and GGA with self-consistent Hubbard parameters (on-site and inter-site ) calculated from linear response theory, the experimentally observed phase stability trends are recovered. The calculated on-site between Mn- states in the experimentally observed orthorhombic, layered, and spinel phases are significantly smaller than in and disordered layered structures, by within GGA. The smaller values of are shown to be correlated with a collinear ordering of JT distortions, in which all orbitals are oriented in the same direction. This cooperative JT effect can lead to greater electron delocalization from Mn along the states due to increased Mn-O covalency, which contributes to the greater electronic stability compared to the phases with noncollinear JT arrangements. The structures with collinear ordering of JT distortions also generate greater vibrational entropy, which helps stabilize these phases at high temperature. These phases are shown to be strongly insulating with large calculated band gaps , which are computed using HSE06 and .

Kam, Ronald L↗
DOE OSTI2022

Machine learning for design principles for single atom catalysts towards electrochemical reactions

Machine learning (ML) integrated density functional theory (DFT) calculations have recently been used to accelerate the design and discovery of heterogeneous catalysts such as single atom catalysts (SACs) through the establishment of deep structure–activity relationships. Here, this review provides recent progress in the ML-aided rational design of heterogeneous catalysts with the focus on SACs in terms of structure–activity relationships, feature importance analysis, high-throughput screening, stability, and metal–support interactions for electrochemistry. Support vector machine (SVM), random forest regression (RFR), and deep neural networks (DNN) along with atomic properties are mainly used for the design of SACs. The ML results have shown that the number of electrons in the d orbital, oxide formation enthalpy, ionization energy, Bader charge, d-band center, and enthalpy of vaporization are mainly the most important parameters for the defining of the structure–activity relationships for electrochemistry. However, the black-box nature of ML techniques occasionally makes a physical interpretation of descriptors, such as the Bader charge, d-band center, and enthalpy of vaporization, non-trivial. At the current stage, ML application is limited by the lack of a large and high-quality database. Future prospects for the development of a large database and a generalized ML algorithm for SAC design are discussed to give insights for further studies in this field.

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

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2023

Probing lithium mobility at a solid electrolyte surface

Abstract Solid-state electrolytes overcome many challenges of present-day lithium ion batteries, such as safety hazards and dendrite formation 1,2 . However, detailed understanding of the involved lithium dynamics is missing due to a lack of in operando measurements with chemical and interfacial specificity. Here we investigate a prototypical solid-state electrolyte using linear and nonlinear extreme-ultraviolet spectroscopies. Leveraging the surface sensitivity of extreme-ultraviolet-second-harmonic-generation spectroscopy, we obtained a direct spectral signature of surface lithium ions, showing a distinct blueshift relative to bulk absorption spectra. First-principles simulations attributed the shift to transitions from the lithium 1 s state to hybridized Li- s /Ti- d orbitals at the surface. Our calculations further suggest a reduction in lithium interfacial mobility due to suppressed low-frequency rattling modes, which is the fundamental origin of the large interfacial resistance in this material. Our findings pave the way for new optimization strategies to develop these electrochemical devices via interfacial engineering of lithium ions.

25 ENERGY STORAGE↗
DOE OSTI2021

COMPASO : A new halo finder for competitive assignment to spherical overdensities

We describe a new method (compaso) for identifying groups of particles in cosmological N-body simulations. compaso builds upon existing spherical overdensity (SO) algorithms by taking into consideration the tidal radius around a smaller halo before competitively assigning halo membership to the particles. In this way, the compaso finder allows for more effective deblending of haloes in close proximity as well as the formation of new haloes on the outskirts of larger ones. This halo-finding algorithm is used in the abacussummit suite of N-body simulations, designed to meet the cosmological simulation requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. compaso is developed as a highly efficient on-the-fly group finder, which is crucial for enabling good load-balancing between the GPU and CPU and the creation of high-resolution merger trees. In this paper, we describe the halo-finding procedure and its particular implementation in abacus, accompanying it with a qualitative analysis of the finder. We test the robustness of the compaso catalogues before and after applying the cleaning method described in an accompanying paper and demonstrate its effectiveness by comparing it with other validation techniques. We then visualize the haloes and their density profiles, finding that they are well fit by the NFW formalism. Finally, we compare other properties such as radius–mass relationships and two-point correlation functions with that of another widely used halo finder, rockstar.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

O VI Emission from the Interstellar Medium

This project's primary goal was to examine the Local Bubble, a large hot bubble surrounding the solar neighborhood. In order to do this, we observed the 1032 and 1038 A resonance line emission from O VI in the bubble and used the,results to comment on models for the Local Bubble and its embedded clouds. In order to maximize the signal to noise of our spectrum, we combined the awarded guest investigator observation with unpublished FUSE In Orbit Checkout observations. The resulting spectrum was sufficiently good as to enable us to place tight 2 sigma upper limits on the intensities of the 1032 and the 1038 A resonance lines. We also measured or placed upper limits on the other cosmic lines in the bandpass, including C III and C II. These are the first known ultraviolet emission line measurements and/or upper limits for the gas in the Local Bubble (as opposed to gas anywhere along long lines of sight). With the O VI upper limits, we were able to quantitatively evaluate competing theories for the origins of the Local Bubble. The upper limits are well below those expected in the Breitschwerdt model (which proposes that during its its early development, the Local Bubble rapidly expanded beyond its nascent cloud and, as a result, is now vastly underionized). The upper limits on the O VI resonance line doublet intensity and the measurement of the C III intensity, garnered from this project, combined with measurements of the O VI column density, garnered from another project, are so far below the predictions, that they make a good case for eliminating the Breitschwerdt model from the field of possibilities. Thus, instead of being vastly underionized, the Local Bubble is near ionizational equilibrium. In addition, the upper limits challenge the other well-known model for the Local Bubble. In that model, the Local Bubble was blown by a series of supernova explosions and winds and contains a myriad of evaporating clouds. The intensity of the O VI resonance line doublet predicted from models exceeds the observationally determined 2 sigma upper limits. The discrepancy prompts the exciting suggestion that there may be some astrophysics yet to be considered in standard hydrodynamic models. These measurements performed a second useful function. By subtract- ing them from the O VI intensity observed on long sight-lines through the Galactic halo and Local Bubble, we were able to estimate the intensity originating the Galactic halo.

Shelton, Robin L.↗
DOE OSTI2025

Understanding the Active Site Structures and Achieving Catalytic Activity Tuning of Atomically Dispersed FeN 4 Sites for Oxygen Reduction Reaction

Atomically dispersed Fe–—N—C catalysts with high oxygen reduction reaction (ORR) activity have attracted great attention since the last decade. Due to comparable ORR activity and low material cost, they are promising platinum group metal (PGM)-free catalysts that can replace the commercialized Pt/C materials; furthermore, it can facilitate the efficiency of the fuel cell technologies and mitigate dependence on fossil fuels. Great advancements have been made to experimentally optimize the synthesis approach of the Fe–—N—C catalysts, enhance the ORR activity, and improve the catalyst stability. Similarly, recent theoretical studies also provide enriched understanding of the active site structures, properties, and reaction mechanisms. In this review, discussions are made upon utilizing combined experimental and computational spectroscopy to reveal the active site structures, employing mechanistic studies to investigate reaction thermodynamics and kinetics, as well as developing scaling relationships to assist the design and development of future PGM-free catalyst materials. Furthermore, recent advances in studying Fe–—N—C catalysts utilizing electrified surface models and explicit solvation models are also discussed. Not only can these aspects improve the accuracy of theoretical simulation and predictions but also deepen the understanding of the catalyst properties and reaction mechanisms under the effect of surface charges and solvent molecules.

Fe single-atom catalysts↗
DOE OSTI2026

Stoichiometry dependent properties of cerium hydride: An active learning developed interatomic potential study

Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

0D and 2D: The Cases of Phenylethylammonium Tin Bromide Hybrids

Tin halide perovskites and perovskite-related materials have emerged as promising lead-free hybrid materials for various optoelectronic applications. While remarkable progress has been achieved in the development of organic tin halide hybrids with diverse structures and controlled dimensionalities at the molecular level, some controversial results that have been reported recently need to be addressed. For instance, different photophysical properties have been reported for two-dimensional (2D) (PEA) 2 SnBr 4 (PEA = phenylethylammonium) by several groups with distinct emission peaks at around 468 and 550 nm. Here we report our efforts in the synthesis of phenylethylammonium tin bromide hybrids with zero-dimensional (0D) and 2D structures, and characterizations of their structural and photophysical properties. 0D [(PEA) 4 SnBr 6 ][(PEA)Br] 2 [CCl 2 H 2 ] 2 was found to exhibit strong yellow emission peak at 566 nm with a photoluminescence quantum efficiency (PLQE) of ~90%, while 2D (PEA) 2 SnBr 4 had weak emission peak at 470 nm with a PLQE of <0.1%. Interestingly, 0D [(PEA) 4 SnBr 6 ][(PEA)Br] 2 [CCl 2 H 2 ] 2 can be converted into 2D (PEA) 2 SnBr 4 upon drying, which would return to 0D [(PEA) 4 SnBr 6 ][(PEA)Br] 2 [CCl 2 H 2 ] 2 upon addition of dichloromethane. Powder X-ray diffraction results confirmed the reversible transformation between 0D and 2D structures. Density functional theory calculations showed that excitons in 0D [(PEA) 4 SnBr 6 ][(PEA)Br] 2 [CCl 2 H 2 ] 2 are highly localized, resulting in a strongly Stokes shifted broadband emission, while delocalized electronic states in 2D (PEA) 2 SnBr 4 result in weaker exciton binding, a higher exciton mobility, and a higher nonradiative decay.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pt Nanoparticles on Atomic-Metal-Rich Carbon for Heavy-Duty Fuel Cell Catalysts: Durability Enhancement and Degradation Behavior in Membrane Electrode Assemblies

Proton exchange membrane fuel cells (PEMFCs) are a promising zero-emission power source for heavy-duty vehicles (HDVs). However, long-term durability of up to 25,000 h is challenging because current carbon support, catalyst, membrane, and ionomer developed for traditional light-duty vehicles cannot meet the stringent requirement. Therefore, understanding catalyst degradation mechanisms under the HDV condition is crucial for rationally designing highly active and durable platinum group metal (PGM) catalysts for high-performance membrane electrode assemblies (MEAs). Herein, we report a PGM catalyst consisting of platinum nanoparticles with a high content (40 wt %) on atomic-metal-site (e.g., MnN 4 )-rich carbon support. MEAs with the Pt (40 wt %)/Mn–N–C cathode catalyst achieved significantly enhanced performance and durability, generating 1.41 A cm –2 at 0.7 V under HDV conditions (0.25 mgPt cm –2 and 250 kPa abs pressure) and retaining 1.20 A cm –2 after an extended and accelerated stress test up to 150,000 voltage cycles. Electron microscopy studies indicate that most fine Pt nanoparticles are retained on or/and in the carbon support covered with the ionomer throughout the catalyst layer at the end of life. During the long-term stability test, the observed electrochemical active surface area reduction and performance loss primarily result from Pt depletion in the catalyst layer due to Pt dissolution and redeposition at the interface of the cathode and membrane. Importantly, the first-principle density functional theory calculations further reveal a support entrapment effect of the Mn–N–C, in which the MnN 4 site can specifically adsorb the Pt atom and further retard the Pt dissolution and migration, therefore enhancing long-term MEA durability.

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

Bridging the Gap between the X-ray Absorption Spectroscopy and the Computational Catalysis Communities in Heterogeneous Catalysis: A Perspective on the Current and Future Research Directions

X-ray absorption spectroscopy (XAS) [extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES)] is a key technique within the heterogeneous catalysis community to probe the structure and properties of the active site(s) for a diverse range of catalytic materials. Furthermore, the interpretation of the raw experimental data to derive an atomistic picture of the catalyst requires modeling and analysis; the EXAFS data are compared to a model, and a goodness of fit parameter is used to judge the best fit. This EXAFS modeling can often be nontrivial and time-consuming; overcoming or improving these limitations remains a central challenge for the community. Considering these limitations, this Perspective highlights how recent developments in analysis software, increased availability of reliable computational models, and application of data science tools can be used to improve the speed, accuracy, and reliability of EXAFS interpretation. In particular, we emphasize the advantages of combining theory and EXAFS as a unified technique that should be treated as a standard (when applicable) to identify catalytic sites and not two separate complementary methods. Building on the recent trends in the computational catalysis community, we also present a community-driven approach to adopt FAIR Guiding Principles for the collection, analysis, dissemination, and storage of XAS data. Written with both the experimental and theory audience in mind, we provide a unified roadmap to foster collaborations between the two communities.

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

Crystalline solutions of the Kohn-Sham equations in the fractional quantum Hall regime

A Kohn-Sham density functional approach has recently been developed for the fractional quantum Hall effect, which maps the strongly interacting electrons into a system of weakly interacting composite fermions subject to an exchange correlation potential as well as a density dependent gauge field that mimics the "flux quanta" bound to composite fermions. To get a feel for the role of various terms, we study the behavior of the self-consistent solution as a function of the strength of the exchange correlation potential, which is varied through an ad hoc multiplicative factor. Here, we find that a crystal phase is stabilized when the exchange correlation interaction is sufficiently strong relative to the composite-fermion cyclotron energy. Various properties of this crystal are examined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Computational and Experimental Mechanistic Insights into the Ethanol-to-Butanol Upgrading Reaction over MgO

The mechanism of ethanol upgrading to higher products is still under debate, especially regarding intermediate species and hydrogenation and dehydrogenation steps. In this work, we conducted a combined theoretical and experimental approach to contribute to this discussion. For such, detailed electronic structure density functional theory calculations (aiming at probing density of states, infrared spectra, geometric parameters, charge densities, and reaction energetics) and diffuse reflectance infrared Fourier transform spectroscopy experiments were carried out revealing the relevance of an appropriate combination of reactive surface sites to support the formation of several intermediates that are formed in the C-C coupling over MgO. The roles of Mg and O sites were also studied under an electronic perspective and different geometrical arrangements. We found that a kink configuration was the most adequate for ethanol to 1-butanol upgrading. Our calculations also gave us arguments to propose distinct reaction routes, whose mutual predominance would depend upon reaction temperature. At temperatures up to 573 K, the so-called β-route, which goes through scission of a Cβ-H bond and formation of an oxametallacycle-like intermediate, would dominate the coupling, whereas at higher temperatures, up to 673 K, a more usual Guerbet mechanism, via an aldol coupling step and then consecutive hydrogenations, would be expected. Here, the theoretical conclusions were followed by a careful experimental strategy using sequential experimental planning techniques in order to estimate accurate parameters with the lowest possible experimental load. Information from these different sources were coupled to develop a mathematical model for the rate of the ethanol upgrading reaction, using a Langmuir-Hinshelwood-Hougen-Watson approach. The developed and statistically validated model adequately described the experimental data at 673 K and 1.1 bar total pressure for ethanol partial pressures in the range from 0 to 20 kPa.

09 BIOMASS FUELS↗
DOE OSTI2021

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ion Mobility Spectrometry Characterization of the Intermediate Hydrogen-Containing Gold Cluster Au7(PPh3)7H52+

We employ ion mobility spectrometry and density functional theory to determine the structure of Au7(PPh3)7H52+ (PPh3 = triphenylphosphine) which was recently identified by high mass resolution mass spectrometry. Experimental ion-neutral collision cross sections represent the momentum transfer between the ionic clusters and gas molecules averaged over the relative thermal velocities of the colliding pair, thereby providing structural insights. Theoretical calculations indicate the geometry of Au7(PPh3)7H52+ is similar to Au7(PPh3)7+, with three hydrogen atoms bridging two gold atoms and two hydrogen atoms forming single Au–H bonds. Collision induced dissociation products observed during IMS experiments reveal that smaller hydrogen containing clusters may be produced through fragmentation of Au7(PPh3)7H52+. Our findings indicate that hydrogen containing species like Au7(PPh3)7H52+ act as intermediates in the formation of larger phosphine ligated gold clusters. These results advance the understanding and ability to control the mechanisms of size selective cluster formation which is necessary for scalable synthesis of clusters with tailored properties. This work was supported by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. M.A.H. acknowledges support from the DOE Science Undergraduate Laboratory Internship (SULI) program. H.H. acknowledges support from the DOE Office of Workforce Development for Teachers and Scientist (WDTS) under the Visiting Faculty Program (VFP). The efforts of K.A.M. and B.H.C. were supported by the Defense Threat Reduction Agency (Grant No. HDTRA1-14-1-0023). The research was performed using EMSL, a national scientific user facility sponsored by the DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is operated by Battelle for DOE under Contract DE-AC05-76RL01830.

Ligare, Marshall R.↗
NASA NTRS2001

Retrieval of Polar Stratospheric Cloud Microphysical Properties From Lidar Measurements: Dependence on Particle Shape Assumptions

Knowledge of particle sizes and number densities of polar stratospheric clouds (PSCs) is highly important, because they are critical parameters for the modeling of the ozone chemistry of the stratosphere. In situ measurements of PSC particles are rare. the main instrument for the accumulation of PSC data are lidar systems. Therefore the derivation of some microphysical properties of PSCS from the optical parameters measured by lidars would be highly beneficial for ozone research. Inversion of lidar data obtained in the presence of PSCs formed from crystalline particles type 11 and the various nitric acid tri Ydrrate (NAT) types cannot be easily accomplished, because a suitable scattering theory for small faceted crystals has not been readily available tip to now. As a consequence, the T-matrix method is commonly used for the interpretation of these PSC lidar data. Here the assumption is made that the optical properties of an ensemble of spheroids resemble those of crystalline PSCs, and microphysical properties of the PSC are inferred from the optical signatures of the PSC at two or more wavelengths. The problem with the T-matrix approach is that the assumption of spheroidal instead of faceted particles can lead to dramatically wrong results: Usually cloud particle properties are deduced from analysis of lidar profiles of backscatter ratio and depolarization ratio. The particle contribution to the backscatter ratio is given by the product of the particle number density and the backscattering cross section. The latter is proportional to the value of the particle's scattering phase function at 180 degrees scattering angle. At 180 degrees however, the phase functions of rough, faceted crystals and of spheroids with same maximum dimension differ by a factor of 6. From this it follows that for a PSC consisting of faceted crystals, the particle number density is underestimated by roughly the same factor if spheroidal particles are unrealistically assumed. We are currently developing a retrieval technique for determining the microphysical parameters of crystalline PSCs that takes into account the faceted shape of the PSC particles. This approach utilizes finite-difference time-domain (FDTD) calculations of particle optical properties. The accuracy and the free choice of the shape of the scattering particle make the FDTD technique a promising tool for the inversion of PSC lidar data. A first comparison of FDTD and T-matrix calculations will be presented.

Reichardt, Susanne↗
DOE OSTI2024

QuTree: A tree tensor network package

Here we present QuTree, a C++ library for tree tensor network approaches. QuTree provides class structures for tensors, tensor trees, and related linear algebra functions that facilitate the fast development of tree tensor network approaches such as the multilayer multiconfigurational time-dependent Hartree approach or the density matrix renormalization group approach and its various extensions. We investigate the efficiency of relevant tensor and tensor network operations and show that the overhead for managing the network structure is negligible, even in cases with a million leaves and small tensors. QuTree focuses on providing simple, high-level routines while retaining easy access to the backend to facilitate novel developments. We demonstrate the capabilities of the package by computing the eigenstates of coupled harmonic oscillator Hamiltonians and performing random circuit simulations on a virtual quantum computer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗