Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nonadiabatic couplings from a variational excited state method based on constrained DFT

Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.

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

Interaction models and configurational entropies of binary MoTa and the MoNbTaW high entropy alloy

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. Furthermore, these interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's cluster variation method entropy formalism is applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included. A high temperature series expansion aids in our understanding of the convergence. Computer codes implementing these methods, and supporting data, are freely available on the internet.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Path integral Monte Carlo calculation of the momentum distribution of the homogeneous electron gas at finite temperature

Path integral Monte Carlo (PIMC) simulations are employed to calculate the momentum distribution of the homogeneous electron gas at finite temperature. Open paths are introduced to sample off-diagonal elements of the real-space density matrix. Here it is demonstrated how the restricted PIMC method can be extended to incorporate open paths in order to allow for the simulations in fermionic systems where a sign problem is present. The computed momentum distribution shows significant deviations from behavior of free fermions when strong correlations are present but agrees with predictions from variational methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Size Effect of Local Current-Voltage Characteristics of MX 2 Nanoflakes: Local Density of States Reconstruction from Scanning Tunneling Microscopy Experiments

Local current-voltage characteristics for low-dimensional transition-metal dichalcogenides (LDTMD), as well as the reconstruction of their local density of states (LDOS) from scanning tunneling microscopy (STM) experiments, are of fundamental interest and can be useful for advanced applications. Most of the existing models either have limited applicability for complex-shaped LDTMDs (e.g., those based on the Simmons approach) or require solving of an ill-defined integral equation to deconvolute the unknown LDOS (e.g., those based on the Tersoff approach). Using a serial expansion of the Tersoff formulas, we propose a flexible method to reconstruct the LDOS from local current-voltage characteristics measured in STM experiments. We establish a set of key physical parameters, which characterize the tunneling current of a STM-probe–sample contact and the sample LDOS expanded in Gaussian functions. Using a direct variational method coupled with probabilistic analysis, we determine these parameters from the STM experiments for MoS2 nanoflakes with different numbers of layers. The main result is the reconstruction of the LDOS in a relatively wide energy range around the Fermi level, which allows us to gain insight into the local band structure of LDTMDs. The reconstructed LDOS reveal pronounced size effects for the single-layer, two-layer, and three-layer MoS 2 nanoflakes, which we relate to the low dimensionality and strong bending or corrugation of the nanoflakes. We hope that the proposed elaboration of the Tersoff approach, allowing LDOS reconstruction, will be of critical interest for the quantitative description of STM experiments and also of use to better understand the microscopic physical aspects of the surface, strain, and bending contributions to the LDTMDs’ electronic properties.

42 ENGINEERING↗
DOE OSTI2015

Resonances in Coupled-Channel Scattering of Light Mesons from Lattice QCD

Recent calculations of light meson scattering will be presented. Several advances have made it possible to extract coupled-channel hadron-hadron scattering amplitudes in a systematically improvable way using lattice QCD. The methods used in some of these first calculations will be explained in detail, and specific applications will be described. I will focus on channels involving pairs of ?, K, ?,?? mesons. Finite volume energy levels describing these systems are obtained from large correlation matrices that are then diagonalized using the variational method to obtain a spectrum. These finite volume spectra are then used to constrain infinite volume scattering amplitudes that may in turn be analyzed for their pole content to extract resonance properties.

Wilson, David↗
DOE OSTI2017

Variational principles for dissipative (sub)systems, with applications to the theory of linear dispersion and geometrical optics

Applications of variational methods are typically restricted to conservative systems. Some extensions to dissipative systems have been reported too but require ad hoc techniques such as the artificial doubling of the dynamical variables. We propose a different approach. Here, we show that for a broad class of dissipative systems of practical interest, variational principles can be formulated using constant Lagrange multipliers and Lagrangians nonlocal in time, which allow treating reversible and irreversible dynamics on the same footing. A general variational theory of linear dispersion is formulated as an example. Particularly, we present a variational formulation for linear geometrical optics in a general dissipative medium, which is allowed to be nonstationary, inhomogeneous, anisotropic, and exhibit both temporal and spatial dispersion simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Few-nucleon matrix elements in pionless effective field theory in a finite volume

Here, pionless effective field theory in a finite volume (FVEFT $\notπ$ ) is investigated as a framework for the analysis of multinucleon spectra and matrix elements calculated in lattice QCD (LQCD). By combining FVEFT π with the stochastic variational method, the spectra of nuclei with atomic number A ∈ {2, 3} are matched to existing finite-volume LQCD calculations at heavier-than-physical quark masses corresponding to a pion mass m $\notπ$ = 806 MeV , thereby enabling infinite-volume binding energies to be determined using infinite-volume variational calculations. Based on the variational wave functions that are constructed in this approach, the finite-volume matrix elements of various local operators are computed in FVEFT $\notπ$ and matched to LQCD calculations of the corresponding QCD operators in the same volume, thereby determining the relevant one- and two-body effective field theory counterterms and enabling an extrapolation of the LQCD matrix elements to infinite volume. As examples, the scalar, tensor, and axial matrix elements are considered, in addition to the magnetic moments and the isovector longitudinal momentum fraction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Measurement of neutron and gamma-ray noise in large subcritical condition for criticality safety - 14485

Various experimental methods have been proposed for subcriticality measurement and neutron noise analysis method is one of important experimental techniques for this purpose. The variance to mean value method, called Feynman-alpha method, is one of well-known noise analysis methods. In this method, variation of neutron count rate distribution from the Poisson distribution caused by fission chain reactions, which is called Y-value, is measured, however, it is difficult to measure the Y-value in large subcritical state such as nuclear fuel facility, because number of emitted neutrons by fission chain reaction is small, which causes large statistical error. We developed new experimental technique based on Feynman-alpha method combined with pulsed neutron source, called Delayed Neutron Noise Analysis (DNNA) method, to monitor the change of subcriticality. In this new method, neutron and gamma-ray counts measured only appeared in the delayed neutron region after injection of pulsed neutron are analyzed, and the prompt neutron region is completely ignored. The advantage of DNNA compared with conventional method is that spatial dependency of detector positions on measured results can be greatly decreased, which result in increase of accuracy of measured subcriticality. The validity of the DNNA method was experimentally investigated in the Kyoto University Critical Assembly (KUCA) at Kyoto University Research Reactor Institute. In the experiments, pulsed neutrons by D-T reactions are injected into the core and we used {sup 3}He detectors to measure neutron noise signals and BGO scintillation detectors to measure gamma-ray noise signals, and all data in the delayed neutron region were analyzed based on Feynman-alpha method. From the experimental results, it was found that the Y-values are clearly larger than zero in the case of existence of fuel assemblies and this method can be applied to detection of hidden special nuclear materials for example in a cargo container. (authors)

42 ENGINEERING↗
DOE OSTI2016

Detecting sulphate aerosol geoengineering with different methods

Sulphate aerosol injection has been widely discussed as a possible way to engineer future climate. Monitoring it would require detecting its effects amidst internal variability and in the presence of other external forcings. Here, we investigate how the use of different detection methods and filtering techniques affects the detectability of sulphate aerosol geoengineering in annual-mean global-mean near-surface air temperature. This is done by assuming a future scenario that injects 5 Tg yr -1 of sulphur dioxide into the stratosphere and cross-comparing simulations from 5 climate models. 64% of the studied comparisons would require 25 years or more for detection when no filter and the multi-variate method that has been extensively used for attributing climate change are used, while 66% of the same comparisons would require fewer than 10 years for detection using a trend-based filter. This then highlights the high sensitivity of sulphate aerosol geoengineering detectability to the choice of filter. With the same trend-based filter but a non-stationary method, 80% of the comparisons would require fewer than 10 years for detection. This does not imply sulphate aerosol geoengineering should be deployed, but suggests that both detection methods could be used for monitoring geoengineering in global, annual mean temperature should it be needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes

Small-angle neutron scattering (SANS) is a powerful tool for studying biological membranes and model lipid bilayer membranes. The length scales probed by SANS, being from 1 nm to over 100 nm, are well-matched to the relevant length scales of the bilayer, particularly when it is in the form of a vesicle. However, it is the ability of SANS to differentiate between isotopes of hydrogen as well as the availability of deuterium labeled lipids that truly enable SANS to reveal details of membranes that are not accessible with the use of other techniques, such as small-angle X-ray scattering. In this work, an overview of the use of SANS for studying unilamellar lipid bilayer vesicles is presented. The technique is briefly presented, and the power of selective deuteration and contrast variation methods is discussed. Approaches to modeling SANS data from unilamellar lipid bilayer vesicles are presented. Finally, recent examples are discussed. While the emphasis is on studies of unilamellar vesicles, examples of the use of SANS to study intact cells are also presented.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Establishing reference ranges for circulating biomarkers of drug‐induced liver injury in healthy human volunteers 1

Aims The potential of mechanistic biomarkers to improve prediction of drug‐induced liver injury (DILI) and hepatic regeneration is widely acknowledged. We sought to determine reference intervals for new biomarkers of DILI and regeneration, as well as to characterize their natural variability and impact of diurnal variation. Methods Serum samples from 227 healthy volunteers were recruited as part of a cross‐sectional study; of these, 25 subjects had weekly serial sampling over 3 weeks, while 23 had intensive blood sampling over a 24h period. Alanine aminotransferase (ALT), MicroRNA‐122 (miR‐122), High Mobility Group Box‐1 (HMGB1), total Keratin‐18 (K18), caspase‐cleaved Keratin‐18 (ccK18), Glutamate Dehydrogenase (GLDH) and Macrophage Colony‐Stimulating Factor‐1 (CSF‐1) were assayed. Results Reference intervals were established for each biomarker based on the 97.5% quantile (90% CI) following the assessment of fixed effects in univariate and multivariable models. Intra‐individual variability was found to be non‐significant, and there was no significant impact of diurnal variation. Conclusion Reference intervals for novel DILI biomarkers have been described. An upper limit of a reference range might represent the most appropriate mechanism to utilize these data. These data can now be used to interpret data from exploratory clinical DILI studies and to assist their further qualification as required by regulatory authorities.

Jorgensen, Andrea L. [Department of Health Data Sc↗
DOE OSTI2014

Adaptive Replanning to Account for Lumpectomy Cavity Change in Sequential Boost After Whole-Breast Irradiation

Purpose: To evaluate the efficiency of standard image-guided radiation therapy (IGRT) to account for lumpectomy cavity (LC) variation during whole-breast irradiation (WBI) and propose an adaptive strategy to improve dosimetry if IGRT fails to address the interfraction LC variations. Methods and Materials: Daily diagnostic-quality CT data acquired during IGRT in the boost stage using an in-room CT for 19 breast cancer patients treated with sequential boost after WBI in the prone position were retrospectively analyzed. Contours of the LC, treated breast, ipsilateral lung, and heart were generated by populating contours from planning CTs to boost fraction CTs using an auto-segmentation tool with manual editing. Three plans were generated on each fraction CT: (1) a repositioning plan by applying the original boost plan with the shift determined by IGRT; (2) an adaptive plan by modifying the original plan according to a fraction CT; and (3) a reoptimization plan by a full-scale optimization. Results: Significant variations were observed in LC. The change in LC volume at the first boost fraction ranged from a 70% decrease to a 50% increase of that on the planning CT. The adaptive and reoptimization plans were comparable. Compared with the repositioning plans, the adaptive plans led to an improvement in target coverage for an increased LC case (1 of 19, 7.5% increase in planning target volume evaluation volume V{sub 95%}), and breast tissue sparing for an LC decrease larger than 35% (3 of 19, 7.5% decrease in breast evaluation volume V{sub 50%}; P=.008). Conclusion: Significant changes in LC shape and volume at the time of boost that deviate from the original plan for WBI with sequential boost can be addressed by adaptive replanning at the first boost fraction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Latent Twins

Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.

Latent Twins↗
DOE OSTI2020

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Multileaf collimator performance monitoring and improvement using semiautomated quality control testing and statistical process control

Purpose: High-quality radiation therapy using highly conformal dose distributions and image-guided techniques requires optimum machine delivery performance. In this work, a monitoring system for multileaf collimator (MLC) performance, integrating semiautomated MLC quality control (QC) tests and statistical process control tools, was developed. The MLC performance monitoring system was used for almost a year on two commercially available MLC models. Control charts were used to establish MLC performance and assess test frequency required to achieve a given level of performance. MLC-related interlocks and servicing events were recorded during the monitoring period and were investigated as indicators of MLC performance variations. Methods: The QC test developed as part of the MLC performance monitoring system uses 2D megavoltage images (acquired using an electronic portal imaging device) of 23 fields to determine the location of the leaves with respect to the radiation isocenter. The precision of the MLC performance monitoring QC test and the MLC itself was assessed by detecting the MLC leaf positions on 127 megavoltage images of a static field. After initial calibration, the MLC performance monitoring QC test was performed 3–4 times/week over a period of 10–11 months to monitor positional accuracy of individual leaves for two different MLC models. Analysis of test results was performed using individuals control charts per leaf with control limits computed based on the measurements as well as two sets of specifications of ±0.5 and ±1 mm. Out-of-specification and out-of-control leaves were automatically flagged by the monitoring system and reviewed monthly by physicists. MLC-related interlocks reported by the linear accelerator and servicing events were recorded to help identify potential causes of nonrandom MLC leaf positioning variations. Results: The precision of the MLC performance monitoring QC test and the MLC itself was within ±0.22 mm for most MLC leaves and the majority of the apparent leaf motion was attributed to beam spot displacements between irradiations. The MLC QC test was performed 193 and 162 times over the monitoring period for the studied units and recalibration had to be repeated up to three times on one of these units. For both units, rate of MLC interlocks was moderately associated with MLC servicing events. The strongest association with the MLC performance was observed between the MLC servicing events and the total number of out-of-control leaves. The average elapsed time for which the number of out-of-specification or out-of-control leaves was within a given performance threshold was computed and used to assess adequacy of MLC test frequency. Conclusions: A MLC performance monitoring system has been developed and implemented to acquire high-quality QC data at high frequency. This is enabled by the relatively short acquisition time for the images and automatic image analysis. The monitoring system was also used to record and track the rate of MLC-related interlocks and servicing events. MLC performances for two commercially available MLC models have been assessed and the results support monthly test frequency for widely accepted ±1 mm specifications. Higher QC test frequency is however required to maintain tighter specification and in-control behavior.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-E-J-262: Variability in Texture Analysis of Gynecological Tumors in the Context of An 18F-FDG PET Adaptive Protocol

Purpose: This study examines the effect on texture analysis due to variable reconstruction of PET images in the context of an adaptive FDG PET protocol for node positive gynecologic cancer patients. By measuring variability in texture features from baseline and intra-treatment PET-CT, we can isolate unreliable texture features due to large variation. Methods: A subset of seven patients with node positive gynecological cancers visible on PET was selected for this study. Prescribed dose varied between 45–50.4Gy, with a 55–70Gy boost to the PET positive nodes. A baseline and intratreatment (between 30–36Gy) PET-CT were obtained on a Siemens Biograph mCT. Each clinical PET image set was reconstructed 6 times using a TrueX+TOF algorithm with varying iterations and Gaussian filter. Baseline and intra-treatment primary GTVs were segmented using PET Edge (MIM Software Inc., Cleveland, OH), a semi-automatic gradient-based algorithm, on the clinical PET and transferred to the other reconstructed sets. Using an in-house MATLAB program, four 3D texture matrices describing relationships between voxel intensities in the GTV were generated: co-occurrence, run length, size zone, and neighborhood difference. From these, 39 textural features characterizing texture were calculated in addition to SUV histogram features. The percent variability among parameters was first calculated. Each reconstructed texture feature from baseline and intra-treatment per patient was normalized to the clinical baseline scan and compared using the Wilcoxon signed-rank test in order to isolate variations due to reconstruction parameters. Results: For the baseline scans, 13 texture features showed a mean range greater than 10%. For the intra scans, 28 texture features showed a mean range greater than 10%. Comparing baseline to intra scans, 25 texture features showed p <0.05. Conclusion: Variability due to different reconstruction parameters increased with treatment, however, the majority of texture features showed significant changes during treatment independent of reconstruction effects.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Solid-State Solvation and Enhanced Exciton Diffusion in Doped Organic Thin Films under Mechanical Pressure

Direct modification of exciton energy has been previously used to optimize the operation of organic optoelectronic devices. One demonstrated method for exciton energy modification is through the use of the solvent dielectric effects in doped molecular films. To gain a deeper appreciation of the underlying physical mechanisms, we test the solid-state solvation effect in molecular thin films under applied external pressure. We observe that external mechanical pressure increases dipole–dipole interactions, leading to shifts in the Frenkel exciton energy and enhancement of the time-resolved spectral red shift associated with the energy-transfer-mediated exciton diffusion. Measurements are performed on host:dopant molecular thin films, which show bathochromic shifts in photoluminescence (PL) under increasing pressure. This is in agreement with a simple solvation theory model of exciton energetics with a fitting parameter based on the mechanical properties of the host matrix material. We measure no significant change in exciton lifetime with increasing pressure, consistent with unchanged aggregation in molecular films under compression. However, we do observe an increase in exciton spectral thermalization rate for compressed molecular films, indicating enhanced exciton diffusion for increased dipole–dipole interactions under pressure. The results highlight the contrast between molecular energy landscapes obtained when dipole–dipole interactions are increased by the pressure technique versus the conventional dopant concentration variation methods, which can lead to extraneous effects such as aggregation at higher doping concentrations. The present work demonstrates the use of pressure-probing techniques in studying energy disorder and exciton dynamics in amorphous molecular thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Evaluation of molecular volume change of block copolymer depending on temperature: A SANS study

Amphiphilic Pluronic triblock copolymers form various self-assembled structures such as sphere, cylinder, lamellae and so on, depending on temperature, leading to the increase of hydrophobicity of block copolymers. However, the effective molecular volume change of the block copolymer has not been fully exploited yet, when temperature increases. Here in this paper, we have investigated the effective molecular volume change of the block copolymer upon heating by using the contrast variation small angle neutron scattering. The scattering length densities (SLDs) of the block copolymer were experimentally obtained from the neutron scattering contrast variation method between the solvent and the block copolymer at varying temperature. Even though the SLD, which is the intrinsic property of the material, should not be changed by temperature elevation, it was dependent on temperature, indicating that the molecular volume is changed. Therefore, we obtained the increase rate of the molecular volume change of the block copolymer (the effective molecular volume change) from the comparison of the calculated SLD and the standard SLD, which is evaluated by plotting the SANS intensity at the first order Bragg peak as the function of temperature at each volume fraction of D 2 O and H 2 O that is about 25.5%–51.3% depending on temperature.

36 MATERIALS SCIENCE↗