Search NASA⌕ Search

SEARCH · Search NASA

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,225 records · Page 68

Sm 2 Ru 3 Sn 5 : A Noncentrosymmetric Cubic Member of the Ln 2 M 3 X 5 Family

An optimized synthetic method is presented for Sm 2 Ru 3 Sn 5 and investigate its physical properties and electronic structure. Sm 2 Ru 3 Sn 5 is prepared by arc-melting stoichiometric ratios of the elements and is confirmed by single crystal and powder X-ray diffraction. An antiferromagnetic transition is observed at T N = 3.8 K. A modified Curie-Weiss fit to the data in the range 50–150 K yields a Curie-Weiss temperature: θ CW = −36.6 K and an effective magnetic moment: μ eff = 0.83 μ B , in agreement with a Sm 3+ oxidation state. Field-dependent magnetization up to H = 7 T at 2 K shows a maximum response of 0.06 μ B , which is significantly lower than the expected Sm 3+ saturation moment (0.71 μ B ). Resistivity measurements indicate metallic behavior, and analysis of the magnetic entropy from the heat capacity reveals a doublet ground state due to crystal electric field splitting. The electronic structure and density of states are calculated with density function theory and further supported by the local density approximation with dynamical mean-field theory. Finally, the experimental and computational results highlight localized Sm 3+ moments and suggest a possible interplay between Ruddelman–Kitel–Kasuya–Yosida and Kondo interactions, positioning Sm 2 Ru 3 Sn 5 as a promising material for studying topology and complex physical phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Orientation-dependent structural properties during growth and growth mechanism of CoO films

The orientation-dependent local structural properties of CoO films on sapphire substrates during growth were investigated through linearly-polarized extended X-ray absorption fine structure (EXAFS) measurements. Specifically, CoO(1 1 1) and (100) crystals with a rock-salt structure (Fm3m) were epitaxially grown on alpha-Al2O3(0001) and (1012) substrates, respectively, at 700 degrees C using a radio-frequency sputtering system. The local structural properties of CoO films in the in-plane and out-of-plane orientations were quantitatively determined using linearly-polarized EXAFS at the Co K-edge during growth. The EXAFS analysis revealed that during the initial stages of growth, the local structural properties exhibit significant differences compared to thick films, with short atomic distances and large (small) Debye-Waller factors observed in the out-of-plane (in-plane) orientations. The local structural strain mostly diminished when approximately 20 CoO layers accumulated on the substrate. Density functional theory (DFT) calculations further supported these findings by confirming that cobalt atoms initially form stable bonds with the sapphire surface, leading to the simultaneous growth of oxygen and cobalt layers in a coordinated manner through layer-by-layer growth.

36 MATERIALS SCIENCE↗

Formulation and solution approach for calibrating activity-based travel demand model-system via microsimulation

This study addresses the problem of calibrating utility-maximizing nested logit activity-based travel demand model-systems. After estimation, it is common practice to use aggregate measurements to calibrate the estimated model-system’s parameters prior to their application in transportation planning, policy making, and operations. However, calibration of activity-based model-systems has received much less attention. Existing calibration approaches are myopic heuristics in the sense that they do not consider the fundamental inter-dependencies among choice-models and do not have a systematic way to adjust model parameters. Also, other purely simulation-based approaches do not perform well in large-scale applications. In this study, we focus on utility-maximizing nested logit activity-based model-systems and calibrating aggregate statistics such as activity shares, mode shares, time-dependent & mode-specific OD flows, and time-dependent & mode-specific sensor counts. We formulate the calibration problem as a simulation-based optimization problem and propose a stochastic gradient-based solution procedure to solve it. The solution procedure relies on microsimulation to calculate expectations of the aggregate statistics of interest to the calibration problem. Additionally, we derive approximate analytical expressions for the gradient of the objective function —that are evaluated through microsimulation on mini-batches of the population. The proposed solution procedure is sensitive to the fundamental structure of the activity-based model-system and is non-myopic in considering the dependencies across its model components. The formulated optimization problem is non-convex, highly nonlinear, and potentially has multiple-minima. Lastly, we show —through a real-world application— that the proposed solution procedure outperforms other state-of-the-art purely simulation-based optimization approaches in terms of computational efficiency, stability, and convergence. We also compare various gradient-based solution algorithms to determine the best algorithm to update the parameters. This work has the potential to facilitate wider and easier application of activity-based model-systems.

97 MATHEMATICS AND COMPUTING↗

Optimization problems governed by systems of PDEs with uncertainties

This paper reviews current theoretical and numerical approaches to optimization problems governed by partial differential equations (PDEs) that depend on random variables or random fields. Such problems arise in many engineering, science, economics and societal decision-making tasks. This paper focuses on problems in which the governing PDEs are parametrized by the random variables/fields, and the decisions are made at the beginning and are not revised once uncertainty is revealed. Examples of such problems are presented to motivate the topic of this paper, and to illustrate the impact of different ways to model uncertainty in the formulations of the optimization problem and their impact on the solution. A linear–quadratic elliptic optimal control problem is used to provide a detailed discussion of the set-up for the risk-neutral optimization problem formulation, study the existence and characterization of its solution, and survey numerical methods for computing it. Different ways to model uncertainty in the PDE-constrained optimization problem are surveyed in an abstract setting, including risk measures, distributionally robust optimization formulations, probabilistic functions and chance constraints, and stochastic orders. Furthermore, approximation-based optimization approaches and stochastic methods for the solution of the large-scale PDE-constrained optimization problems under uncertainty are described. Some possible future research directions are outlined.

Heinkenschloss, Matthias [Rice Univ., Houston, TX ↗

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New Dimension in Ab Initio Electronic Structure Theory: Temperature, Pressure, and Chemical Potential

Ab initio electronic structure theory has transformed gas-phase molecular science with its predictive ability. In the attempt to bring such predictive ability to macroscopic systems and condensed matter, the theory must integrate quantum mechanics with statistical thermodynamics, so that thermodynamic functions such as free energy, internal energy, entropy, and chemical potentials are computed as functions of temperature in a systematically converging series of approximations. Here, a general, versatile strategy of elevating ab initio electronic structure theory to nonzero temperatures is introduced and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Metal Halide Perovskites Demonstrate Radiation Hardness and Defect Healing in Vacuum

Herein, we subject formamidinium lead iodide films to oxygen-containing gases (flowing O 2 or free diffusion of lab atmosphere), inert gases (flowing He, Ar, or N 2 ), and vacuum. Our films are irradiated by Cu Ka X-rays and held at 75 °C while X-ray diffraction is recorded. Under all gas conditions, we observe a reproducible 1.1 ± 0.5 A 3 perovskite lattice contraction from an initial unit cell volume of 256.5 ± 0.8 A 3 concurrent with continuous perovskite loss and lead iodide growth. Oxygen-containing gases increase the reaction rates without materially altering perovskite structural changes. Under the same temperature and irradiation conditions in vacuo, a self-healing reaction is observed, exhibited by a reproducible (0.9 ± 0.3 A 3 ) lattice expansion and stabilization of the perovskite. Furthermore, interactions between the perovskite, defects, and minority phases are simulated by generalized gradient approximation Perdew-Burke-Ernzerhof (GGA-PBE) density functional theory. Lattice contraction indicates an increase in the concentration of Schottky defects-pairs of formamidinium and iodine vacancies. Under irradiation in every atmospheric condition, a solid solution of Schottky defects with a concentration of several percent diffuses and precipitates forming lead iodide and consuming the defects. In the presence of ionized gases, this framework is modified to include the continual loss of formamidinium and iodine ions from the perovskite forming Schottky defects.

14 SOLAR ENERGY↗

Charge Density Wave and Narrow Energy Gap at Room Temperature in 2D Pb 3– x Sb 1+ x S 4 Te 2-δ with Square Te Sheets

We report a new two-dimensional compound Pb 3– x Sb 1+ x S 4 Te 2-δ has a charge density wave (CDW) at room temperature. The CDW is incommensurate with q-vector of 0.248(6)a* + 0.246(8)b* + 0.387(9)c* for x = 0.29(2) and d = 0.37(3) due to positional and occupational long range ordering of Te atoms in the sheets. The modulated structure was refined from the single crystal X-ray diffraction data with a superspace group P$\bar{1}$(αβγ)0 using (3 + 1)-dimensional crystallography. The resistivity increases with decreasing temperature, suggesting semiconducting behavior. The transition temperature (T CDW ) of the CDW is ~ 345 K above which the Te square sheets become disordered with no q-vector. Lastly, first-principles density functional theory calculations on the undistorted structure and an approximate commensurate supercell reveal that the gap is due to the structure modulation.

36 MATERIALS SCIENCE↗

Interpreting Tafel behavior of consecutive electrochemical reactions through combined thermodynamic and steady state microkinetic approaches

Assessing the reaction pathway of multi-electron-transfer reactions is an essential yet difficult task for the rational design of electrocatalysts. In this work, we develop a heuristic approach that combines thermodynamic adsorption energetics calculated through density functional theory with microkinetic modeling using the steady state approximation to interpret the potential-dependent Tafel behavior of consecutive electrochemical reactions. In doing so, we introduce a kinetic framework for ab initio calculations that ensures self-consistent adsorption energetics based on kinetically limited adsorbate coverages. The approach is applied to experimental results on CoO x (OH) 2-x single crystal electrocatalyst particles yielding coverage dependent mechanistic information and identification of the rate-limiting step with standard rate constants for the oxygen evolution reaction on the ($11\bar{2}0$) surfaces of the β-Co(OH) 2 , β-CoOOH, and CoO 2 bulk phases. Finally, this generalizable method enables catalyst benchmarking based on determining the active species involved and associated intrinsic reaction rate constants in consecutive multi-electron-transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An experimental and theoretical study of the C 1s ionization satellites in CH 3 I

We report the C 1s ionization spectrum of CH 3 I has been studied both experimentally and theoretically. Synchrotron radiation has been employed to record polarization dependent photoelectron spectra at a photon energy of 614 eV. These spectra encompass the main-line due to the C 1s single-hole state and the peaks associated with the shake-up satellites. Vertical ionization energies and relative photoelectron intensities have been computed using the fourth-order algebraic-diagrammatic construction approximation scheme for the one-particle Green’s function and the 6-311++G** basis set. The theoretical spectrum derived from these calculations agrees qualitatively with the experimental results, thereby allowing the principal spectral features to be assigned. According to our calculations, two 2 A 1 shake-up states of the C 1s -1 σ CI → σ$*\atop{CI}$ type with singlet and triplet intermediate coupling of the electron spins (S' = 0, 1) play an important role in the spectrum and contribute significantly to the overall intensity. Both of these states are expected to have dissociative diabatic potential energy surfaces with respect to the C–I separation. Whereas the upper of these states perturbs the manifold of Rydberg states, the lower state forms a band which is characterized by a strongly increased width. Our results indicate that the lowest shake-up peak with significant spectral intensity is due to the pair (S' = 0, 1) of 2 E (C 1s -1 I 5p → σ$*\atop{CI}$) states. We predict that these 2 E states acquire photoelectron intensity due to spin-orbit interaction. Lastly, such interactions play an important role here due to the involvement of the I 5p orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recent advances in the quantification of uncertainties in reaction theory

Uncertainty quantification has become increasingly more prominent in nuclear physics over the past several years. In few-body reaction theory, there are four main sources that contribute to the uncertainties in the calculated observables: the effective potentials, approximations made to the few-body problem, structure functions, and degrees of freedom left out of the model space. In this work, we illustrate some of the features that can be obtained when modern statistical tools are applied in the context of nuclear reactions. This work consists of a summary of the progress that has been made in quantifying theoretical uncertainties in this domain, focusing primarily on those uncertainties coming from the effective optical potential as well as their propagation within various reaction theories. We use, as the central example, reactions on the doubly-magic stable nucleus 40 Ca, namely neutron and proton elastic scattering and single-nucleon transfer 40 Ca(d,p) 41 Ca. First, we show different optimization schemes used to constrain the optical potential from differential cross sections and other experimental constraints; we then discuss how these uncertainties propagate to the transfer cross section, comparing two reaction theories. Finally, we finish by laying out our future plans.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MOSAIC: a joint modeling methodology for combined circadian and non-circadian analysis of multi-omics data

Abstract Motivation Circadian rhythms are approximately 24-h endogenous cycles that control many biological functions. To identify these rhythms, biological samples are taken over circadian time and analyzed using a single omics type, such as transcriptomics or proteomics. By comparing data from these single omics approaches, it has been shown that transcriptional rhythms are not necessarily conserved at the protein level, implying extensive circadian post-transcriptional regulation. However, as proteomics methods are known to be noisier than transcriptomic methods, this suggests that previously identified arrhythmic proteins with rhythmic transcripts could have been missed due to noise and may not be due to post-transcriptional regulation. Results To determine if one can use information from less-noisy transcriptomic data to inform rhythms in more-noisy proteomic data, and thus more accurately identify rhythms in the proteome, we have created the Multi-Omics Selection with Amplitude Independent Criteria (MOSAIC) application. MOSAIC combines model selection and joint modeling of multiple omics types to recover significant circadian and non-circadian trends. Using both synthetic data and proteomic data from Neurospora crassa, we showed that MOSAIC accurately recovers circadian rhythms at higher rates in not only the proteome but the transcriptome as well, outperforming existing methods for rhythm identification. In addition, by quantifying non-circadian trends in addition to circadian trends in data, our methodology allowed for the recognition of the diversity of circadian regulation as compared to non-circadian regulation. Availability and implementation MOSAIC’s full interface is available at https://github.com/delosh653/MOSAIC. An R package for this functionality, mosaic.find, can be downloaded at https://CRAN.R-project.org/package=mosaic.find. Supplementary information Supplementary data are available at Bioinformatics online.

De los Santos, Hannah↗

Detecting deviations from Gaussianity in high-redshift CMB lensing maps

ABSTRACT While the probability density function of the cosmic microwave background (CMB) convergence field approximately follows a Gaussian distribution, primordial non-Gaussianities, and small contributions from structures at low redshifts make the overall distribution slightly non-Gaussian. Some of the late-time component can be modelled using the distribution of galaxies and subtracted off from the original CMB lensing map to produce a map of matter distribution at high redshifts. Using this high-redshift mass map, we are able to directly study the early phases of structure formation. In this work, we forecast the detectability of signatures of non-Gaussianity due to non-linear structure formation at z > 1.2. Assuming the optimal case of no systematics, we find that it is challenging to detect such signatures in current surveys, but future experiments such as the deep field of CMB-S4 will be able to make detections of ∼7σ.

Astronomy & Astrophysics↗

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↗

Ground-state properties of the Hubbard model in one and two dimensions from the Gutzwiller conjugate gradient minimization theory

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Influence of hydrostatic pressure on hidden order, the Kondo lattice, and magnetism in URu 2 Si 2-x P x

Within the chemical substitution series URu 2 Si 2-x P x , there is an evolution in the ground-state behavior from hidden ordered (HO) for x ≲ 0.03, to Kondo lattice behavior with no ordering (NO) for 0.03 ≲ x ≲ 0.26, to antiferromagnetism (AFM-2) for 0.26 ≲ x ≲ 0.5 [A. Gallagher et al., Nat. Commun. 7, 10712 (2016); A. Gallagher et al., J. Phys.: Condens. Matter 29, 024004 (2016)]. To better understand what factors control this behavior, temperature-dependent electrical resistivity measurements are performed for this series under applied pressures P up to 20.5 kbar. Specimens in the HO x region show similarities to the parent compound, where HO transforms into antiferromagnetism (AFM-1) at a critical pressure (P c ). P c decreases with increasing x and collapses towards P = 0 near x ≈ 0.03, suggesting that AFM-1 occurs at ambient pressure for this concentration. No pressure-induced phase transitions are observed in the NO x region and the AFM-2 state is only weakly suppressed by P. Measurements further reveal that AFM-1 and AFM-2 are distinct from each other. Calculations of the wave functions using the tight-binding Hartree-Fock approximation are performed and show (i) that the radial probability distributions for the phosphorus ions are more tightly bound than those for the silicon and (ii) that the energy difference between the orbitals decreases with increasing x. The cumulative effect of these two factors is that Si → P substitution decreases the hybridization strength, which correlates with the weakening of HO. At large x, additional effects such as electrical charge tuning also play an important role in determining the ground-state

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Compositional patterning in irradiated alloys: Effective potentials and effective interfacial energy

Compositional patterning (CP) in binary alloys during energetic particle irradiation is studied using a kinetic model that considers two competing kinetic processes, a thermally activated one promoting macroscopic phase separation (MPS) of the concentration field c(r,t) and a forced one resulting in finite-range random atomic mixing. The forced mixing is modeled by a Gaussian relocation distribution with a characteristic distance R. A series of approximate kinetic models are introduced by expanding the mixing function into a series of n terms, thus replacing the non-local evaluations of the concentration field c(r'-r,t) by local derivatives of c(r,t). This approach makes it possible to obtain exact effective potentials and build steady-state diagrams for each order-n model. Phase field (PF) simulations using these order-n models reveal that near the onset of patterning, phase evolution is accurately described using an order-3 model, which changes smoothly from an extended Cahn-Hilliard free energy in the MPS regime to a one-mode Swift-Hohenberg functional in the CP regime. Deeper into the patterning regime, higher-order models are required to achieve convergence, yielding square-like concentration profiles characteristic of a strong segregation regime. These higher-order effective free energies are analogous to multimodal Swift-Hohenberg functionals. Here, a new definition for the effective interfacial energy is proposed in the CP regime, since the interfacial area is no longer an excess quantity in that regime, precluding the use of the standard thermodynamic definition of interfacial energy.

36 MATERIALS SCIENCE↗

Model of chiral spin liquids with Abelian and non-Abelian topological phases

In this article, we present a two-dimensional lattice model for quantum spin-1/2 for which the low-energy limit is governed by four flavors of strongly interacting Majorana fermions. We study this low-energy effective theory using two alternative approaches. The first consists of a mean-field approximation. The second consists of a random phase approximation (RPA) for the single-particle Green's functions of the Majorana fermions built from their exact forms in a certain one-dimensional limit. The resulting phase diagram consists of two competing chiral phases, one with Abelian and the other with non-Abelian topological order, separated by a continuous phase transition. Remarkably, the Majorana fermions propagate in the two-dimensional bulk, as in the Kitaev model for a spin liquid on the honeycomb lattice. We identify the vison fields, which are mobile (they are static in the Kitaev model) domain walls propagating along only one of the two space directions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗