Time-dependent density functional theory with the orthogonal projector augmented wave method
Explore the source record for details and available documents.
Nguyen, Minh [Sandia National Laboratories (SNL-NM↗
SEARCH · Search NASA
Results for “TIME FUNCTION”
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.
Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion
Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.
36 MATERIALS SCIENCE↗
Codebase release r1.1 for SmoQyDEAC.jl
We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl package. It also provides detailed benchmarks of the package’s output against the popular maximum entropy and stochastic analytic continuation methods. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDEAC.jl or installed using the Julia package manager. The online documentation, including examples, can be accessed at https://smoqysuite.github.io/SmoQyDEAC.jl/stable/.
Neuhaus, James (ORCID:0000000169048510)↗
Nanoscale Editing of Multi and Single Layer Tungsten Disulfide via Gas‐Assisted Focused Electron Beam Induced Etching for Device Prototyping
Focused electron beam induced etching (FEBIE) with XeF 2 (xenon difluoride) precursor is conducted on multi-layer exfoliated WS 2 (tungsten disulfide) and monolayer WS 2 grown by chemical vapor deposition (CVD). The films are characterized by atomic force microscopy (AFM) and Raman and photoluminescence (PL) spectroscopy post-etching. The etch rates/efficiencies are reported as a function of electron beam energy, current, dwell time, and XeF 2 pressure. Bulk film Raman spectra are unchanged post-FEBIE, indicating minimal subsurface damage. Monolayer WS 2 shows a decrease in Raman and PL intensity post-FEBIE, with a dose-to-clear of ≈2 nC µm −2 . The study reveals regimes affected by the various mass transport contributions such as refresh time and the ratio of electrons/XeF 2 . Spontaneous etching was discovered during FEBIE of large patterned areas due to the long frame/refresh times. Density functional theory and ab initio molecular dynamics simulations compares desorption of SF x and WF x molecules from pristine WS 2 basal planes and pore edges, revealing the spontaneous etching is consistent with etching of partially etched monolayers during each frame. Single-line etching width of 21 nm, and patterning flakes into 100 nm wide channels are demonstrated. In conclusion, this work demonstrates the possibility of editing WS 2 flakes into electronic devices of arbitrary dimensions for semiconductor applications.
2D materials↗
Amorphization and siliconization of silicon carbide as a first wall material
The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Post-irradiation Heating Tests of As-Irradiated AGR-3/4 TRISO Fuel Compacts
Four post-irradiation heating tests of fuel compacts from the U.S. Advanced Gas Reactor (AGR)-3/4 irradiation experiment were completed. In addition to tristructural isotropic (TRISO)-coated driver fuel, each compact contained designed-to-fail (DTF) particles with fuel kernels coated only in pyrocarbon so as to simulate exposed kernels. Tests at 1600/1700°C, 1400°C, and 1200°C were performed to measure fission product release as a function of time and temperature. Silver release was highest in the 1200°C test, supporting the observation that silver release rates are highest in the 1100–1300°C range. Compared to tests of AGR-1 compacts with no exposed kernels, the Cs-134 and Kr-85 releases were noticeably higher in AGR-3/4. The exposed kernels’ contributions to Eu and Sr release are inconclusive, due to the difficulty in distinguishing among the combined effects of higher irradiation temperatures in these particular AGR-3/4 compacts, the presence of the DTF particles, and the Fuel Accident Condition Simulator (FACS) test temperatures. These data can be used to make inferences about fission product retention in exposed kernels as a function of time and temperature.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Novel application of neutrinos to evaluate U.S. nuclear weapons performance
There is a growing realization that neutrinos can be used as a diagnostic tool to better understand the inner workings of a nuclear weapon. Robust estimates demonstrate that an Inverse Beta Decay (IBD) neutrino scintillation detector built at the Nevada Test Site with a 1000-ton active target mass at a standoff distance of 500 m would detect thousands of antineutrino events per nuclear test. This would provide less than 4% statistical error on the measured antineutrino rate and 5% error on antineutrino energy. Extrapolating this to an error on the test device explosive yield requires knowledge from evaluated nuclear databases, non-equilibrium fission rates, and assumptions on internal neutron fluxes. Initial calculations demonstrate that the total number of neutrinos emitted per fission in the first 10 3 s after a short pulse of 239 Pu fission is about a factor of two less than that from Pu fissioning under steady state conditions. Furthermore, there are significant energy spectral differences as a function of time after the pulse that must be considered. These and other model dependencies will be discussed in the paper. In the absence of nuclear weapons testing, many of the technical and theoretical challenges of a full nuclear test could be mitigated with a low cost smaller scale 20 ton fiducial mass IBD demonstration detector placed near a pulsed reactor. Potential reactors include the Texas A&M University TRIGA 1 GW–10 ms pulsed facility or the Sandia Annular Core Research Reactor. The short duty cycle and repeatability of pulses would provide critical real environment testing and measurements, which would be valuable for planning a possible real test shot in the future. Furthermore, the antineutrino rate as a function of time data would provide unique constraints on fission databases and model assumptions. Finally, there are impactful science drivers such as sensitive searches for ∼1 eV 2 sterile neutrinos and ∼MeV scale axions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Lignin Extraction and Condensation as a Function of Temperature, Residence Time, and Solvent System in Flow-through Reactors
Solvolytic extraction of lignin from biomass is a critical step in lignin-first biorefining, including the reductive catalytic fractionation (RCF) process. Key to optimal RCF processing is the ability to rapidly extract lignin from biomass at high delignification extents and transfer the lignin molecules to a catalyst surface in a time frame that minimizes lignin condensation reactions. Here, we use a flow-through reactor to study the effects of temperature (175-250 °C), residence time (9 to 36 min), and solvent composition (methanol and methanol-water) on lignin extraction and condensation. We evaluated three metrics at each condition: total delignification, delignification rate, and extent of condensation, the latter measured by a decrease in monomer yield for batch hydrogenolysis reactions of solvolysis liquor compared to batch RCF reactions. We observe that delignification is predominantly determined by temperature, while residence time dictates the lignin condensation extent. Moreover, the extent of both extraction and condensation increased in the methanol-water solvent system compared to that in the methanol system. Lignin extracted in methanol is stable up to 18-min residence times at or below 225 °C, while a majority of the lignin extracted in methanol-water is condensed with a 9-min residence time at 200 °C. These results can inform reactor designs and solvent selection for lignin-first biorefining processes that aim to physically separate the biomass and catalyst.
09 BIOMASS FUELS↗
Lignin Extraction and Condensation as a Function of Temperature, Residence Time, and Solvent System in Flow-through Reactors
Solvolytic extraction of lignin from biomass is a critical step in lignin-first biorefining, including the reductive catalytic fractionation (RCF) process. Key to optimal RCF processing is the ability to rapidly extract lignin from biomass at high delignification extents and transfer the lignin molecules to a catalyst surface in a time frame that minimizes lignin condensation reactions. Here, we use a flow-through reactor to study the effects of temperature (175–250 °C), residence time (9 to 36 min), and solvent composition (methanol and methanol–water) on lignin extraction and condensation. We evaluated three metrics at each condition: total delignification, delignification rate, and extent of condensation, the latter measured by a decrease in monomer yield for batch hydrogenolysis reactions of solvolysis liquor compared to batch RCF reactions. We observe that delignification is predominantly determined by temperature, while residence time dictates the lignin condensation extent. Moreover, the extent of both extraction and condensation increased in the methanol–water solvent system compared to that in the methanol system. Lignin extracted in methanol is stable up to 18-min residence times at or below 225 °C, while a majority of the lignin extracted in methanol–water is condensed with a 9-min residence time at 200 °C. These results can inform reactor designs and solvent selection for lignin-first biorefining processes that aim to physically separate the biomass and catalyst.
biorefining↗
Lignin Extraction and Condensation as a Function of Temperature, Residence Time, and Solvent System in Flow-through Reactors -Figure data
Data for Figures from Lignin extraction paper, https://pubs.acs.org/doi/10.1021/acssuschemeng.5c04198
biorefining↗
Lignin Extraction and Condensation as a Function of Temperature, Residence Time, and Solvent System in Flow-through Reactors -Tables
tabulated data for Lignin extraction paper, https://pubs.acs.org/doi/10.1021/acssuschemeng.5c04198
biorefining↗
Rotational memory function of SPC/E water
Memory effects are essential for the dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for water have never been calculated directly from molecular dynamics simulations. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The normalized memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, ≲1 fs, justifies the use of the mathematics of rotational diffusion to describe the dynamics of a single molecular dipole moment in bulk water. Here, an analytical equation for the rotational memory time is derived.
Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
Fouling behavior of zwitterionic membranes compared to polyamide membranes
Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.
Cross-linked polyamide↗
Time-domain theory of transient heat conduction in the local limit
Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Elastic Bayesian Model Calibration
Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.
97 MATHEMATICS AND COMPUTING↗
Correlation functions from tensor network influence functionals: The case of the spin-boson model
We investigate the application of matrix product state (MPS) representations of the influence functionals (IFs) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We examine three different IF approaches: one based on the Kadanoff–Baym contour targeting correlation functions at all times, one based on a complex contour targeting the correlation function at a single time, and a steady state formulation, which avoids imaginary or complex times, while providing access to correlation functions at all times. We show that within the IF language, the steady state formulation provides a powerful approach to evaluate equilibrium correlation functions.
Chemistry↗
Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes
We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Utah FORGE CoreFlooding Experimental Results
This submission contains associated data from 100C and 200C core-flooding experiments conducted by Lawrence Livermore National Laboratory. The samples used were sourced from 16A(78)-32 well core. The primary objectives of these tests were to determine the change in calculated hydraulic fracture and permeability over time, under constant confining pressure, in predominantly constant flowrate conditions, while monitoring effluent chemistry as a function of time.
15 GEOTHERMAL ENERGY↗