Search NASA⌕ Search

SEARCH · Search NASA

Results for “CASE METHOD”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Continuity of the maximum-entropy inference: Convex geometry and numerical ranges approach

We study the continuity of an abstract generalization of the maximum-entropy inference—a maximizer. It is defined as a right-inverse of a linear map restricted to a convex body which uniquely maximizes on each fiber of the linear map a continuous function on the convex body. Using convex geometry we prove, amongst others, the existence of discontinuities of the maximizer at limits of extremal points not being extremal points themselves and apply the result to quantum correlations. Further, we use numerical range methods in the case of quantum inference which refers to two observables. One result is a complete characterization of points of discontinuity for 3 × 3 matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Solution of D dimensional Dirac equation for hyperbolic tangent potential using NU method and its application in material properties

The analytical solution of D-dimensional Dirac equation for hyperbolic tangent potential is investigated using Nikiforov-Uvarov method. In the case of spin symmetry the D dimensional Dirac equation reduces to the D dimensional Schrodinger equation. The D dimensional relativistic energy spectra are obtained from D dimensional relativistic energy eigen value equation by using Mat Lab software. The corresponding D dimensional radial wave functions are formulated in the form of generalized Jacobi polynomials. The thermodynamically properties of materials are generated from the non-relativistic energy eigen-values in the classical limit. In the non-relativistic limit, the relativistic energy equation reduces to the non-relativistic energy. The thermal quantities of the system, partition function and specific heat, are expressed in terms of error function and imaginary error function which are numerically calculated using Mat Lab software.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Locally coupled open subsystems: A formalism for affordable electronic structure calculations featuring fractional charges and size consistency

This manuscript introduces a methodology (within the Born-Oppenheimer picture) to compute electronic ground-state properties of molecules and solids/surfaces with fractionally occupied components. Given a user-defined division of the molecule into subsystems, our theory uses an auxiliary global Hamiltonian that is defined as the sum of subsystem Hamiltonians, plus the spatial integral of a second-quantized local operator that allows the electrons to be transferred between subsystems. This electron transfer operator depends on a local potential that can be determined using density functional approximations and/or other techniques such as machine learning. The present framework employs superpositions of tensor-product wave functions, which can satisfy size consistency and avoid spurious fractional charges at large bond distances. The electronic population of each subsystem is in general a positive real number and is obtained from wave-function amplitudes, which are calculated by means of ground-state matrix diagonalization (or matrix propagation in the time-dependent case). Furthermore, our method can provide pathways to explore charge-transfer effects in environments where dividing the molecule into subsystems is convenient and to develop computationally affordable electronic structure algorithms.

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

Methods to accelerate high-throughput screening of atomic qubit candidates in van der Waals materials

The discovery of atom-like spin emitters associated with defects in two-dimensional (2D) wide-bandgap (WBG) semiconductors presents new opportunities for highly tunable and versatile qubits. So far, the study of such spin emitters has focused on defects in hexagonal boron nitride (hBN). However, hBN necessarily contains a high density of nuclear spins, which are expected to create a strong incoherent spin-bath that leads to poor coherence properties of spins hosted in the material. Therefore, identification of new qubit candidates in other 2DWBG materials is necessary. Given the time demands of ab initio methods, new approaches for rapid screening and calculations of identifying properties of suitable atom-like qubits are required. In this work, we present two new methods for rapid estimation of the zero-phonon line (ZPL), a key property of atomic qubits in WBG materials. First, the ZPL is calculated by exploiting Janak’s theorem. For finite changes in occupation, we provide the leading-order estimate of the correction to the ZPL obtained using Janak’s theorem, which is more rapid than the standard method (ΔSCF). Next, we demonstrate an approach to converging excited states that is faster for systems with small strain than the standard approach used in the ΔSCF method. We illustrate these methods using the case of the singly negatively charged calcium vacancy in SiS2, which we are the first to propose as a qubit candidate. This work has the potential to assist in accelerating the high-throughput search for quantum defects in materials, with applications in quantum sensing and quantum computing.

Kuate Defo, Rodrick (ORCID:0000000242788196)↗
DOE OSTI2021

Integrated atomic force microscopy and x-ray irradiation for in situ characterization of radiation-induced processes

Understanding radiation-induced chemical and physical transformations at material interfaces is important across diverse fields, but experimental approaches are often limited to either ex situ observations or in situ electron microscopy or synchrotron-based methods in which cases the radiation type and dose are inextricably tied to the imaging basis itself. In this work we overcome this limitation by demonstrating integration of an X-ray source with an atomic force microscope to directly monitor radiolytically-driven interfacial chemistry at the nanoscale. We illustrate the value of in situ observations by examining effects of radiolysis on material adhesion forces in aqueous solution, as well as examining the production of alkali nitrates at the interface between an alkali halide crystal surface and air. For the examined salt-air interface, direct visualization under flexible experimental conditions greatly extends prior observations by enabling the transformation process to be followed comprehensively from source-to-sink with mass balance quantitation. Overall, our novel rad-AFM opens doors into understanding the dynamics of radiolytically-driven mass transfer and surface alteration at the nanoscale in real-time.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

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

In situ X-ray studies of adlayer-induced crystal nucleation at the liquid-liquid interface

Crystal nucleation and growth at a liquid–liquid interface is studied on the atomic scale by in situ Å-resolution X-ray scattering methods for the case of liquid Hg and an electrochemical dilute electrolyte containing Pb 2+ , F - , and Br - ions. In the regime negative of the Pb amalgamation potential Φ rp = -0.70 V, no change is observed from the surface-layered structure of pure Hg. Upon potential-induced release of Pb 2+ from the Hg bulk at Graphic, the formation of an intriguing interface structure is observed, comprising a well-defined 7.6-Å–thick adlayer, decorated with structurally related 3D crystallites. Both are identified by their diffraction peaks as PbFBr, preferentially aligned with their Graphic axis along the interface normal. X-ray reflectivity shows the adlayer to consist of a stack of five ionic layers, forming a single-unit-cell–thick crystalline PbFBr precursor film, which acts as a template for the subsequent quasiepitaxial 3D crystal growth. This growth behavior is assigned to the combined action of electrostatic and short-range chemical interactions.

Elsen, Annika↗
DOE OSTI2020

Radiation Source Localization Using Surrogate Models Constructed from 3-D Monte Carlo Transport Physics Simulations

Recent research has focused on the development of surrogate models for radiation source localization in a simulated urban domain. We employ the Monte Carlo N-Particle (MCNP) code to provide high- delity simulations of radiation transport within an urban domain. The model is constructed to employ a source location (x, y, z) as input and return the estimated count rate for a set of speci ed detector locations. Because MCNP simulations are computationally expensive, we develop e cient and accurate surrogate models of the detector responses. We construct surrogate models using Gaussian processes (GP) and neural networks (NN) that we train and verify using the MCNP simulations. The trained surrogate models provide an e cient framework for Bayesian inference and experimental design. We employ Delayed Rejection Adaptive Metropolis (DRAM), a Markov Chain Monte Carlo (MCMC) algorithm, to infer the location and intensity of an unknown source. The DRAM results yield a posterior probability distribution for the source's location conditioned on the observed detector count rates. The posterior distribution exhibits regions of high and low probability within the simulated environment identifying potential source locations. In this manner, we can quantify the source location to within at least one of these regions of high probability in the considered cases. Employing these methods, we are able to reduce the space of potential source locations by at least 60%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Prolate spheroidal harmonic expansion of gravitational field

As a modification of the oblate spheroidal case, a recursive method is developed to compute the point value and a few low-order derivatives of the prolate spheroidal harmonics of the second kind, Q{sub nm} (y), namely the unnormalized associated Legendre function (ALF) of the second kind with its argument in the domain, 1 < y < ∞. They are required in evaluating the prolate spheroidal harmonic expansion of the gravitational field in addition to the point value and the low-order derivatives of P-bar {sub nm}(t), the 4π fully normalized ALF of the first kind with its argument in the domain, |t| ≤ 1. The new method will be useful in the gravitational field computation of elongated celestial objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

COLLISIONAL EXCITATION OF FAR-INFRARED LINE EMISSIONS FROM WARM INTERSTELLAR CARBON MONOXIDE (CO)

Motivated by recent observations with Herschel/PACS, and the availability of new rate coefficients for the collisional excitation of CO, the excitation of warm astrophysical CO is revisited with the use of numerical and analytic methods. For the case of an isothermal medium, results have been obtained for a wide range of gas temperatures (100-5000 K) and H{sub 2} densities (10{sup 3}-10{sup 9} cm{sup -3}), and presented in the form of rotational diagrams, in which the logarithm of the column density per magnetic substate, log (N{sub J} /g{sub J} ), is plotted for each state, as a function of its energy, E{sub J} . For rotational transitions in the wavelength range accessible to Herschel/PACS, such diagrams are nearly linear when n(H{sub 2}) {>=} 10{sup 8} cm{sup -3}. When n(H{sub 2}) {approx} 10{sup 6.8}-10{sup 8} cm{sup -3}, they exhibit significant negative curvature, whereas when n(H{sub 2}) {<=} 10{sup 4.8} cm{sup -3}, the curvature is uniformly positive throughout the PACS-accessible range. Thus, the observation of a positively curved CO rotational diagram does not necessarily require the presence of multiple temperature components. Indeed, for some sources observed with Herschel/PACS, the CO rotational diagrams show a modest positive curvature that can be explained by a single isothermal component. Typically, the required physical parameters are densities in the 10{sup 4}-10{sup 5} cm{sup -3} range and temperatures close to the maximum at which CO can survive. Other sources exhibit rotational diagrams with more curvature than can be accounted for by a single temperature component. For the case of a medium with a power-law distribution of gas temperatures, dN/dT{proportional_to}T{sup -b}, results have been obtained for H{sub 2} densities 10{sup 3}-10{sup 9} cm{sup -3} and power-law indices, b, in the range 1-5; such a medium can account for a CO rotational diagram that is more positively curved than any resulting from an isothermal medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE NEGATIVE EFFECTIVE MAGNETIC PRESSURE IN STRATIFIED FORCED TURBULENCE

To understand the basic mechanism of the formation of magnetic flux concentrations, we determine by direct numerical simulations the turbulence contributions to the mean magnetic pressure in a strongly stratified isothermal layer with large plasma beta, where a weak uniform horizontal mean magnetic field is applied. The negative contribution of turbulence to the effective mean magnetic pressure is determined for strongly stratified forced turbulence over a range of values of magnetic Reynolds and Prandtl numbers. Small-scale dynamo action is shown to reduce the negative effect of turbulence on the effective mean magnetic pressure. However, the turbulence coefficients describing the negative effective magnetic pressure phenomenon are found to converge for magnetic Reynolds numbers between 60 and 600, which is the largest value considered here. In all these models, the turbulent intensity is arranged to be nearly independent of height, so the kinetic energy density decreases with height due to the decrease in density. In a second series of numerical experiments, the turbulent intensity increases with height such that the turbulent kinetic energy density is nearly independent of height. Turbulent magnetic diffusivity and turbulent pumping velocity are determined with the test-field method for both cases. The vertical profile of the turbulent magnetic diffusivity is found to agree with what is expected based on simple mixing length expressions. Turbulent pumping is shown to be down the gradient of turbulent magnetic diffusivity, but it is twice as large as expected. Corresponding numerical mean-field models are used to show that a large-scale instability can occur in both cases, provided the degree of scale separation is large enough and hence the turbulent magnetic diffusivity small enough.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ON DETECTING TRANSIENT PHENOMENA

Transient phenomena are interesting and potentially highly revealing of details about the processes under observation and study that could otherwise go unnoticed. It is therefore important to maximize the sensitivity of the method used to identify such events. In this article, we present a general procedure based on the use of the likelihood function for identifying transients which is particularly suited for real-time applications because it requires no grouping or pre-processing of the data. The method makes use of all the information that is available in the data throughout the statistical decision-making process, and is suitable for a wide range of applications. Here we consider those most common in astrophysics, which involve searching for transient sources, events or features in images, time series, energy spectra, and power spectra, and demonstrate the use of the method in the case of a weak X-ray flare in a time series and a short-lived quasi-periodic oscillation in a power spectrum. We derive a fit statistic that is ideal for fitting arbitrarily shaped models to a power density distribution, which is of general interest in all applications involving periodogram analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Oscillator strength measurement for the A(0–6)–X(0), C(0)–X(0), and E(0)–X(0) transitions of CO by the dipole (γ, γ) method

The oscillator strengths of CO, especially for the excitations to the predissociating levels C{sup 1}Σ{sup +} and E{sup 1}Π, have significant applications in the studies of interstellar gases. A large number of investigations, applying a variety of experimental and calculational methods, have focused on the determination of these oscillator strengths. However, deviations of more than 10% have been observed, even for the same experimental method. In this case, we introduce a new technique, the dipole (γ, γ), to measure the oscillator strengths. The obtained oscillator strengths of the strong bands C{sup 1}Σ{sup +}(ν{sup ′}=0)−X{sup 1}Σ{sup +} and E{sup 1}Π(ν{sup ′}=0)−X{sup 1}Σ{sup +} in this work reach high accuracy and are in excellent agreement with some previous results, although they are different from the recommended values in other works. In addition, we measured the oscillator strengths of the A{sup 1}Π(ν{sup ′}=0−6)−X{sup 1}Σ{sup +} bands, and the electronic transition moment dependence on r-centroid was derived.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Exact restoration of Galilei invariance in density functional calculations with quantum Monte Carlo

We investigate some aspects of the connection between mean-field and ab-initio calculations for any-nucleon systems employing the Gogny effective interaction. In particular, we present variational Monte Carlo calculations of 4 He and 16 O, for which we extended the formalism of the linear method to the case of a density-dependent Hamiltonian. Monte Carlo calculations using an uncorrelated basis reduce to Hartree-Fock ones, but Galilean invariance can be preserved, contrary to the usual treatment. An analysis of the results obtained with the Gogny force when used in conjunction with correlated wave functions has also been performed. Here, the calculations show the collapse of the deuteron, leading to a set of unphysical predictions for the structure of heavier nuclei. This behavior remarks the different role of the Gogny effective interaction as compared to bare interactions typically used in fully microscopic calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Reduced-order atomistic cascade method for simulating radiation damage in metals

Atomistic modeling of radiation damage through displacement cascades is deceptively non-trivial. Because of the high energy and stochastic nature of atomic collisions, individual primary knock-on atom (PKA) cascade simulations are computationally expensive and ill-suited for length and dose upscaling. Here, we propose a reduced-order atomistic cascade model capable of predicting and replicating radiation events in metals across a wide range of recoil energies. Our methodology approximates cascade and displacement damage production by modeling the cascade as a core-shell atomic structure composed of two damage production estimators, namely an athermal recombination corrected displacements per atom (arc-dpa) in the shell and a replacements per atom (rpa) representing atomic mixing in the core. These estimators are calibrated from explicit PKA simulations and a standard displacement damage model that incorporates cascade defect production efficiency and mixing effects. We highlight the predictability and accuracy of our reduced-order atomistic cascade method for the cases of copper and niobium by comparing its results with those from full PKA simulations in terms of defect production as well as the resulting cascade evolution and structure. We offer examples for simulating high energy cascade fragmentation and large dose ion-bombardment to demonstrate its possible applicability. Finally, we discuss the various practical considerations and challenges associated with this methodology especially when simulating subcascade formation and dose effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Application of State Quantization-Based Methods in HEP Particle Transport Simulation

Simulation of particle-matter interactions in complex geometries is one of the main tasks in high energy physics (HEP) research. An essential aspect of it is an accurate and efficient particle transportation in a non-uniform magnetic field, which includes the handling of volume crossings within a predefined 3D geometry. Quantized State Systems (QSS) is a family of numerical methods that provides attractive features for particle transportation processes, such as dense output (sequences of polynomial segments changing only according to accuracy-driven discrete events) and lightweight detection and handling of volume crossings (based on simple root-finding of polynomial functions). In this work we present a proof-of-concept performance comparison between a QSS-based standalone numerical solver and an application based on the Geant4 simulation toolkit, with its default Runge-Kutta based adaptive step method. In a case study with a charged particle circulating in a vacuum (with interactions with matter turned off), in a uniform magnetic field, and crossing up to 200 volume boundaries twice per turn, simulation results showed speedups of up to 6 times in favor of QSS while it being 10 times slower in the case with zero volume boundaries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Linear point and sound horizon as purely geometric standard rulers

The baryon acoustic oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the purely geometric BAO methods. To date two different purely geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable (e.g., the BAO). We define the linear point and sound horizon and we fully characterize and compare the two rulers’ cosmological-parameter dependence. We find that they are both geometrical (i.e., independent of the primordial cosmological parameters) within the required accuracy, and that they have the same parameter dependence for a wide range of parameter values. We estimate the rulers’ best-fit values and errors, given the cosmological constraints obtained by the Planck Satellite team from their measurements of the cosmic microwave background temperature and polarization anisotropies. We do this for three different cosmological models encompassed by the purely geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by 0.5σ when we do not fix the spatial geometry to be flat in Λ CDM. Finally, this points toward a sensitivity of the rulers to different cosmological models when they are estimated from the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

On the effect of Ga and In substitutions in the Ca 11 Bi 10 and Yb 11 Bi 10 bismuthides crystallizing in the tetragonal Ho 11 Ge 10 structure type

The Ga- and In-substituted bismuthides Ca 11 Ga x Bi 10– x , Ca 11 In x Bi 10– x , Yb 11 Ga x Bi 10– x , and Yb 11 In x Bi 10– x ( x < 2) can be readily synthesized employing molten Ga or In metals as fluxes. They crystallize in the tetragonal space group I 4/ mmm and adopt the Ho 11 Ge 10 structure type (Pearson code tI 84; Wyckoff sequence n 2 m j h 2 e 2 d ). The structural response to the substitution of Bi with smaller and electron-poorer In or Ga has been studied by single-crystal X-ray diffraction methods for the case of Ca 11 In x Bi 10– x [ x = 1.73 (2); octabismuth undecacalcium diindium]. The refinements show that the In atoms substitute Bi only at the 8 h site. The refined interatomic distances show an unconventional – for this structure type – bond-length distribution within the anionic sublattice. The latter can be viewed as consisting of isolated Bi 3- anions and [In 4 Bi 8 20- ] clusters for the idealized Ca 11 In 2 Bi 8 model. Formal electron counting and first-principle calculations show that the peculiar bonding in this compound drives the system toward an electron-precise state, thereby stabilizing the observed bond-length pattern.

Chemistry↗