Search NASA⌕ Search

SEARCH · Search NASA

Results for “Real Functions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multiscale simulations of electron and ion dynamics in self-irradiated silicon

The interaction of energetic ions with the electronic and ionic system of target materials is an interesting but challenging multiscale problem, and understanding of the early stages after impact of heavy, initially charged ions is particularly poor. At the same time, energy deposition during these early stages determines later formation of damage cascades. In this work, we address the multiscale character by combining real-time time-dependent density functional theory for electron dynamics with molecular dynamics simulations of damage cascades. Our first-principles simulations prove that core electrons affect electronic stopping and have an unexpected influence on the charge state of the projectile. We show that this effect is absent for light projectiles, but dominates the stopping physics for heavy projectiles. By parametrizing an inelastic energy loss friction term in the molecular dynamics simulations using our first-principles results, we also show a qualitative influence of electronic stopping physics on radiation-damage cascades.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dynamic short-range correlation in photoinduced disorder phase transitions

Ultrafast photoexcitation can induce a nonequilibrium dynamic with electron-lattice interaction, offering an effective way to study photoinduced phase transitions (PIPTs) in solids. The issue that atomic displacements after photoexcitation belong to a coherent change or disordered process has become controversial in the PIPT community. Using real-time, time-dependent density functional theory (rt-TDDFT) simulations, we were able to obtained both the coherent and the disordered PIPTs (dimer dissociation) in IrTe 2 with different electronic occupations. More importantly, we found that in the disordered phase transition there exists a local correlation between different dimers regarding their dissociation status. We define these Ir-Ir dimers directly connected by Te atoms, including intralayer and vertically across the layers, as a group (group I). Other Ir-Ir dimers separated by five Ir atoms from Ir-Ir dimers in group I are divided into another group (group II). The dimers in the same group will dissociate in a correlated fashion; they either all dissociate or all do not dissociate. On the other hand, the dimers in neighboring groups will have an anticorrelation: If the dimers in one group dissociate, the dimers in the neighboring group tend not to be dissociated, and vice versa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Electronic Excitation Response of DNA to High-Energy Proton Radiation in Water

The lack of molecular-level understanding for the electronic excitation response of DNA to charged particle radiation, such as high-energy protons, remains a fundamental scientific bottleneck in advancing proton and other ion beam cancer therapies. In particular, the dependence of different types of DNA damage on high-energy protons represents a significant knowledge void. Here, in this work, we employ first-principles real-time time-dependent density functional theory simulation, using a massively parallel supercomputer, to unravel the quantum-mechanical details of the energy transfer from high-energy protons to DNA in water. The calculations reveal that protons deposit significantly more energy onto the DNA sugar-phosphate side chains than onto the nucleobases, and greater energy transfer is expected onto the DNA side chains than onto water. As a result of this electronic stopping process, highly energetic holes are generated on the DNA side chains as a source of oxidative damage.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Dielectric screening and vacancy formation for large neutral and charged Si n H m ( n > 1500 ) nanocrystals using real-space pseudopotentials

A commonly used procedure for computing the properties of defects in crystalline materials is to consider a large supercell that includes the defect of interest. This is a straightforward technique as standard energy band codes can be used for such computations. For neutral defects, the only impediment of such an approach is to avoid defect-defect interactions between adjoining cells. However, this procedure can be complex if the defect of interest is charged as the system at large contains Coulombic divergences. Moreover, some have recently argued that the conventional definition of formation energies for charged defects cannot be reconciled with statistical mechanics. Here, we focus on an alternative approach. We consider large nanocrystals wherein a charged defect can be placed. Since the system is confined, a charged defect within the nanocrystal does not result in a Coulombic divergence. The chief impediment is computational, i.e., while no defect-defect or Coulombic divergences are present, the nanocrystal must be sufficiently large to allow the system to properly replicate a bulklike configuration. With the development of new algorithms and hardware advances, computations for systems of sufficient size to address this issue are feasible. In particular, we solve the Kohn-Sham equation in real space using pseudopotential-density-functional theory for large silicon nanocrystals, which contain thousands of atoms. Further, we focus on (i) the screening of a point charge and (ii) the formation of a charged vacancy in hydrogen-terminated silicon nanocrystals. This approach allows us to examine the role of quantum confinement in addition to exploring the bulk limit. Comparisons to other methods confirm the viability of this approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Shape-induced superstructure formation in concentrated ferrofluids under applied magnetic fields

The field-induced ordering of concentrated ferrofluids based on spherical and cuboidal maghemite nanoparticles is studied using small-angle neutron scattering, revealing a qualitative effect of the faceted shape on the interparticle interactions as shown in the structure factor and correlation lengths. Whereas a spatially disordered hard-sphere interaction potential with a short correlation length is found for ~9 nm spherical nanoparticles, nanocubes of a comparable particle size exhibit a more pronounced interparticle interaction and the formation of linear arrangements. Analysis of the anisotropic two-dimensional pair distance correlation function gives insight into the real-space arrangement of the nanoparticles. On the basis of the short interparticle distances found here, oriented attachment, i.e. a face-to-face arrangement of the nanocubes, is likely. The unusual field dependence of the interparticle correlations suggests a field-induced structural rearrangement.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The seeds and homogeneous nucleation of photoinduced nonthermal melting in semiconductors due to self-amplified local dynamic instability

Laser-induced nonthermal melting in semiconductors has been studied over the past four decades, but the underlying mechanism is still under debate. Here, by using an advanced real-time time-dependent density functional theory simulation, we reveal that the photoexcitation-induced ultrafast nonthermal melting in silicon occurs via homogeneous nucleation with random seeds originating from a self-amplified local dynamic instability. Because of this local dynamic instability, any initial small random thermal displacements of atoms can be amplified by a charge transfer of photoexcited carriers, which, in turn, creates a local self-trapping center for the excited carriers and yields the random nucleation seeds. Because a sufficient amount of photoexcited hot carriers must be cooled down to band edges before participating in the self-amplification of local lattice distortions, the time needed for hot carrier cooling is the response for the longer melting time scales at shorter laser wavelengths. This finding provides fresh insights into photoinduced ultrafast nonthermal melting.

42 ENGINEERING↗
DOE OSTI2021

Directional dependency of electronic stopping in nickel, projectile’s excited charge state and momentum transfer

Abstract We studied the directional dependency of electronic stopping power of swift light ions in nickel using real-time time-dependent density functional theory. We report a variation of electronic stopping for moving ions as the projectile probes different electronic densities of the host material. These results show that while the predicted magnitude stays in reasonable agreement with experiment, for $$v > 2$$ v > 2 . a.u. simulating only low index crystallographic directions is not enough to sample the experimental average values. The ab initio simulations give us access to microscopic quantities, such as non-adiabatic forces, momentum transfer and transient excited state charges of the projectile and host ions, which are not available through other methods. We report these quantities for the first time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Quantum Non-Equilibrium Dynamics Prediction of Electronic Transport Coefficients in Non-Linear Regimes

The goal of this work was to develop unique LLNL capabilities for the computational prediction of electrical, optical and thermal transport properties from quantum simulations. Using a recently developed theory of microscopic quantum transport, we were able to model compressed metallic hydrogen and it’s non-linear and out-of-equilibrium contributions in electrical and thermal transport. This theory was implemented in the open-source, real-space time-dependent density functional theory code Octopus. This model is general and widely applicable to high energy density systems that are important to LLNL and NNSA mission areas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Tau-Jet Matching and Tag-and-Probe Tau Tagging Efficiencies for HH→bbτ+τ− Searches Using Run 3 CMS Scouting Data

We are studying tau-jet matching efficiency and tau tagging efficiency for the search for Higgs boson pair production in the HH→bbτ+τ− decay channel using Run 3 CMS scouting data. Collision events in the CMS detector produce many particle candidates, so jets are clustered using the anti-kT algorithm to identify possible tau signatures. In Monte Carlo simulation, tau-jet matching efficiency can be measured by comparing generated taus to reconstructed jets. In actual collision data, however, generated taus are not available, so we use a tag-and-probe method with Z→τ+τ− candidate events. In this method, one tau candidate is used as the tag and the other as the probe. The invariant mass of the tau-pair candidates is plotted and fitted with signal and background functions to estimate the number of real tau events and determine the tagging efficiency. By comparing efficiencies in data and Monte Carlo simulation, this work helps evaluate the performance of scouting-based tau reconstruction and tau tagging for future HH→bbτ+τ− searches.

Bellot, Annella [North Central Coll.; Fermilab]↗
DOE OSTI2026

Tau-Jet Matching and Tag-and-Probe Tau Tagging Efficiencies for HH→bbτ+τ− Searches Using Run 3 CMS Scouting Data

We are studying tau-jet matching efficiency and tau tagging efficiency for the search for Higgs boson pair production in the HH→bbτ+τ− decay channel using Run 3 CMS scouting data. Collision events in the CMS detector produce many particle candidates, so jets are clustered using the anti-kT algorithm to identify possible tau signatures. In Monte Carlo simulation, tau-jet matching efficiency can be measured by comparing generated taus to reconstructed jets. In actual collision data, however, generated taus are not available, so we use a tag-and-probe method with Z→τ+τ− candidate events. In this method, one tau candidate is used as the tag and the other as the probe. The invariant mass of the tau-pair candidates is plotted and fitted with signal and background functions to estimate the number of real tau events and determine the tagging efficiency. By comparing efficiencies in data and Monte Carlo simulation, this work helps evaluate the performance of scouting-based tau reconstruction and tau tagging for future HH→bbτ+τ− searches.

Bellot, Annella [North Central Coll.; Fermilab]↗
DOE OSTI2020

Optimal Renewable Resource Allocation and Load Scheduling of Resilient Communities

This paper presents a methodology for enhancing community resilience through optimal renewable resource allocation and load scheduling in order to minimize unserved load and thermal discomfort. The proposed control architecture distributes the computational effort and is easier to be scaled up than traditional centralized control. The decentralized control architecture consists of two layers: The community operator layer (COL) allocates the limited amount of renewable energy resource according to the power flexibility of each building. The building agent layer (BAL) addresses the optimal load scheduling problem for each building with the allowable load determined by the COL. Both layers are formulated as a model predictive control (MPC) based optimization. Simulation scenarios are designed to compare different combinations of building weighting methods and objective functions to provide guidance for real-world deployment by community and microgrid operators. The results indicate that the impact of power flexibility is more prominent than the weighting factor to the resource allocation process. Allocation based purely on occupancy status could lead to an increase of PV curtailment. Further, it is necessary for the building agent to have multi-objective optimization to minimize unserved load ratio and maximize comfort simultaneously.

model predictive control↗
DOE OSTI2022

Reverse-mode differentiation in arbitrary tensor network format: with application to supervised learning.

This paper describes an efficient reverse-mode differentiation algorithm for contraction operations of tensor networks that may have arbitrary and unconventional network topologies. The approach leverages the tensor contraction tree of Evenbly and Pfeifer (2014), which provides an instruction set for the contraction sequence of a network. We show that this tree can be efficiently leveraged for differentiation of a full tensor network contraction using a recursive scheme that exploits (1) the bilinear property of contraction and (2) the property that trees have single path from root to leaves. While differentiation of tensor-tensor contraction is already possible in most automatic differentiation packages, we show that exploiting these two additional properties in the specific context of contraction sequences can improve efficiency. Following a description of the algorithm and computational complexity analysis, we investigate its utility for gradient-based supervised learning for low-rank function recovery and for fitting real-world unstructured datasets. We demonstrate improved performance over alternating least-squares optimization approaches and the capability to handle heterogeneous and arbitrary tensor network formats. When compared to alternating minimization algorithms, we find that the gradient-based approach requires a smaller oversampling ratio (number of samples compared to number model parameters) for recovery. This increased efficiency extends to fitting unstructured data of varying dimensionality and when employing a variety of tensor network formats. Here, we show improved learning using the hierarchical Tucker method over the tensor-train in high-dimensional settings on a number of benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Efficient generalized boundary detection

Fast, efficient, and robust compression-based methods for detecting boundaries in arbitrary datasets, including sequences (1D datasets), are desired. The methods, each employing three simple algorithms, approximate the information distance between two adjacent sliding windows within a dataset. One of the algorithms calculates an initial ordered list of subsequences; while a second algorithm updates the ordered list of subsequences by dropping a first entry and appending a last entry rather than calculating completely new ordered lists with each iteration. Large values in the distance metric are indicative of boundary locations. A smoothed z-score or a wavelet-based algorithm may then be used to locate peaks in the distance metric, thereby identifying boundary locations. An adaptive version of the method employs a collection of window sizes and corresponding weighting functions, making it more amenable to real datasets with unknown, complex, and changing structures.

Ting, Christina↗
DOE OSTI2023

The Association Between Baseline Hepatic or Renal Function and Clinical Outcomes for Patients With Non–Small Cell Lung Cancer Treated With a PD –1/PD– L1 Blocking Antibody Using Real–World and Trial Data

Clinical trials have demonstrated the benefit of PD-1/PD-L1 blocking antibodies for the treatment of patients with advanced non-small cell lung cancer (NSCLC) in defined patient populations that often exclude patients with moderate or severe hepatic or renal impairment. We assessed the association between overall survival (OS) and baseline organ function in patients with advanced NSCLC treated with PD-1/PD-L1 blocking antibodies in real-world data (RWD; patient-level data from electronic health records) and pooled clinical trial data submitted to the US Food and Drug Administration (FDA). The Kaplan–Meier estimator was used to estimate OS in different subgroups based on organ function. Unadjusted and adjusted Cox proportional hazards models were used to estimate the association between OS and organ function. In this hypothesis-generating study, baseline renal impairment did not appear to be associated with OS, while patients with baseline liver impairment had shorter OS. RWD provided information on a broader range of renal and hepatic function than was evaluated in clinical trials and hold promise to complement trial data in better understanding populations not represented in clinical trials.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

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

Are grid cells used for navigation? On local metrics, subjective spaces, and black holes

The symmetric, lattice-like spatial pattern of grid-cell activity is thought to provide a neuronal global metric for space. This view is compatible with grid cells recorded in empty boxes but inconsistent with data from more naturalistic settings. Here, we review evidence arguing against the global-metric notion, including the distortion and disintegration of the grid pattern in complex and three-dimensional environments. We argue that deviations from lattice symmetry are key for understanding grid-cell function. We propose three possible functions for grid cells, which treat real-world grid distortions as a feature rather than a bug. First, grid cells may constitute a local metric for proximal space rather than a global metric for all space. Second, grid cells could form a metric for subjective action-relevant space rather than physical space. Third, distortions may represent salient locations. Finally, we discuss mechanisms that can underlie these functions. These ideas may transform our thinking about grid cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

All-electron, density-functional-based method for angle-resolved tunneling ionization in the adiabatic regime

Here we develop and test a method that integrates many-electron weak-field asymptotic theory (ME-WFAT) [O. I. Tolstikhin, L. B. Madsen, and T. Morishita, Phys. Rev. A 89, 013421 (2014)] in the integral representation (IR) into the density-functional-theory (DFT) framework. In particular, we present modifications of the integral formula in the IR ME-WFAT to incorporate the potential terms unique to DFT. By solving an adiabatic rate equation for the angle-resolved ionization yield in our DFT-based ME-WFAT method, we show that the results are in excellent agreement with those of real-time time-dependent density-functional-theory (TDDFT) simulations for NO, OCS, C⁢H 3 ⁢Br, and C⁢H 3 ⁢Cl interacting with one- and two-color laser fields with a fundamental wavelength of 800 nm. This agreement is significant because the WFAT calculations take only a small fraction of the time of full TDDFT calculations. These results suggest that in the wavelength region commonly used in strong-field experiments (800 nm and longer), our DFT-based WFAT treatment can be used to rapidly screen for the ionization properties of a large number of molecules as a function of alignment or orientation between the molecule and the strong field.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Logarithmic enhancements in conformal perturbation theory and their real time interpretation

We study various corrections of correlation functions to leading order in conformal perturbation theory, both on the cylinder and on the plane. Many problems on the cylinder are mathematically equivalent to those in the plane if we give the perturbations a position dependent scaling profile. The integrals to be done are then similar to those in the study of correlation functions with one additional insertion at the center of the profile. We will be primarily interested in the divergence structure of these corrections when computed in dimensional regularization. In particular, we show that the logarithmic divergences (enhancements) that show up in the plane under these circumstances can be understood in terms of resonant behavior in time dependent perturbation theory, for a transition between states that is induced by an oscillatory perturbation on the cylinder.

Physics↗