Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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 OSTI2020

Modelling of the laser amplification process with allowance for the effect of the temperature distribution in an Yb : YAG gain element on the thermophysical and lasing characteristics of the medium

A time-dependent three-dimensional model for the laser amplification process has been constructed with allowance for the effect of the temperature distribution on the thermophysical and lasing characteristics of gain media. We have performed numerical modelling of the laser amplification process in the gain elements of a two-stage subjoule-level cryogenic laser amplifier operating at a pulse repetition rate of up to 1 kHz. It has been shown that taking into account the temperature distribution is of critical importance in calculation of cryogenically cooled laser amplifiers pumped with high-power diodes. We have found optimal diode pump parameters at which the maximum achievable pulse energy at the amplifier output can reach 300 and 570 mJ at pulse repetition rates of 1000 and 500 Hz, respectively. (extreme light fields and their interaction with matter)

42 ENGINEERING↗
DOE OSTI2025

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploring the exact limits of the real-time equation-of-motion coupled cluster cumulant Green’s functions

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green’s function approach [Rehr et al., J. Chem. Phys. 152, 174113 (2020)]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz, ensuring the exactness in the expansion limit. In addition, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green’s function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green’s function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green’s function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green’s function approach and other Green’s function approaches in the numerical exact limit.

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

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M –1 cm –1 ) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Benchmarking the multiconfigurational Hartree method by the exact wavefunction of two harmonically trapped bosons with contact interaction

Highlights: • Comparison of exact and numerical single-particle density matrices. • Benchmark for multiconfigurational Hartree. • Convergence for large number of orbitals. We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, as well as overestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Understanding the chemical contribution to the enhancement mechanism in SERS: Connection with Hammett parameters

The enhancement mechanism due to the molecule-surface chemical interaction in surface-enhanced Raman scattering (SERS) has been characterized using a theoretical approach based on time dependent density functional theory. This includes a systematic study of the chemical mechanism (CM) to the SERS enhancement for halogen substituted benzenethiols interacting with a silver cluster. Changing the halogen on benzenethiol enables us to systematically modulate interactions between the benzenethiol ring and the metal cluster. We observe a decrease in the CM enhancement factor with an increase in atomic number of the halogen for para- substitutions. For meta- substitutions, there is no such trend. However, the results scale linearly with the Hammett parameters for both meta and para halogens, which provides an important predictive tool for interpreting chemical enhancements. We also study the effect of solvation on the CM, showing that there is a systematic increase in enhancement with increasing solvent dielectric constant. As a result, the correlation of CM with other properties, such as the amount of charge transfer between adsorbate and metal, or the excitation energies of charge transfer states, is much less predictive than the Hammett parameter correlation.

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

Effect of Substituent Location on the Relationship between the Transition Dipole Moments, Difference Static Dipole, and Hydrophobicity in Squaraine Dyes for Quantum Information Devices

Aggregates of organic dyes that exhibit excitonic coupling have a wide array of applications, including medical imaging, organic photovoltaics, and quantum information devices. The optical properties of a dye monomer, as a basis of dye aggregate, can be modified to strengthen excitonic coupling. Squaraine (SQ) dyes are attractive for those applications due to their strong absorbance peak in the visible range. While the effects of substituent types on the optical properties of SQ dyes have been previously examined, the effects of various substituent locations have not yet been investigated. In this study, density functional theory (DFT) and time-dependent density functional theory (TD-DFT) were used to investigate the relationships between SQ substituent location and several key properties of the performance of dye aggregate systems, namely, difference static dipole (Δd), transition dipole moment (μ), hydrophobicity, and the angle (θ) between Δd and μ. We found that attaching substituents along the long axis of the dye could increase μ while placement off the long axis was shown to increase Δd and reduce θ. The reduction in θ is largely due to a change in the direction of Δd as the direction of μ is not significantly affected by substituent position. Hydrophobicity decreases when electron-donating substituents are located close to the nitrogen of the indolenine ring. These results provide insight into the structure–property relationships of SQ dyes and guide the design of dye monomers for aggregate systems with desired properties and performance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Plasmonic Hot-Carrier Generation and Catalysis in Ti 3 C 2 O 2 from Real-Time TDDFT Simulations

Photoinduced hot electrons are central to plasmon-driven catalysis. Atomically thin Ti 3 C 2 O 2 , with high carrier density and broad optical absorption, offers a promising platform for plasmon-driven reactions. However, comprehensive investigations of its plasmon resonance, hot-carrier generation, and plasmonic catalytic performance remain limited. In this work, real-time time-dependent density functional theory (rt-TDDFT) was employed to study Ti 3 C 2 O 2 ’s plasmon excitation and hot-carrier generation from nonradiative plasmon damping. The temporal evolution of the dipole moment reveals plasmon resonance in Ti 3 C 2 O 2 , followed by strong plasmon damping that redistributes the stored energy to generate hot carriers. Ti 3 C 2 O 2 with low oxygen vacancy concentration (O v -Ti 3 C 2 O 2 ) exhibits plasmonic behavior resembling the pristine surface, and the plasmon-generated hot electrons can markedly reduce the dissociation barrier of CO 2 at the oxygen vacancy. These findings provide fundamental insights into the plasmonic properties of Ti 3 C 2 O 2 and how they drive its catalytic performance in surface reactions, which is valuable for advancing plasmon-driven catalysis.

CO2 reduction↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Tuned and screened range-separated hybrid density functional theory for describing electronic and optical properties of defective gallium nitride

We apply a hybrid density functional theory approach, based on a tuned and screened range-separated hybrid (SRSH) exchange-correlation functional, to describe the optoelectronic properties of defective gallium nitride (GaN). SRSH and time-dependent SRSH (TDSRSH) are tuned to produce accurate energetics for the pristine material and applied to the study of a series of point defects in bulk GaN, a blue light-emitting material that degrades in the presence of defects. We first establish the accuracy of the method by comparing the predicted quasiparticle gap and low-energy excitation spectra of (TD)SRSH and many-body perturbation theory for both pristine GaN and GaN containing a single nitrogen vacancy. Aided by the reduced computational cost of (TD)SRSH, we then report on three additional technologically relevant point defects and defect complexes in GaN: the gallium vacancy, the carbon interstitial, and the carbon-silicon complex. We compute the low-energy optical absorption spectra for these defects and show the presence of defect-centered transitions. Furthermore, by estimating the Stokes shift, we predict, in agreement with previous studies, that the carbon substitutional defect is a candidate for the detrimental yellow luminescence in GaN. Furthermore, this study indicates that TDSRSH is a promising and computationally feasible approach for quantitatively accurate, first-principles modeling of defective semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

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

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

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

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗
DOE OSTI2024

Site-specific electronic structure of covalently linked bimetallic dyads from nitrogen K-edge x-ray absorption spectroscopy

A nitrogen K-edge x-ray absorption near-edge structure (XANES) survey is presented for tetrapyrido[3,2-a:',3'-c:3'',2''-h:2''',3'''-j]phenazine (tpphz)-bridged bimetallic assemblies that couple chromophore and catalyst transition metal complexes for light driven catalysis, as well as their individual molecular constituents. We demonstrate the high N site sensitivity of the N pre-edge XANES features, which are energetically well-separated for the phenazine bridge N atoms and for the individual metal-bound N atoms of the inner coordination sphere ligands. By comparison with the time-dependent density functional theory calculated spectra, we determine the origins of these distinguishable spectral features. We find that metal coordination generates large shifts toward higher energy for the metal-bound N atoms, with increasing shift for 3d < 4d < 5d metal bonding. This is attributed to increasing ligand-to-metal σ donation that increases the effective charge of the bound N atoms and stabilizes the N 1s core electrons. In contrast, the phenazine bridge N pre-edge peak is found at a lower energy due to stabilization of the low energy electron accepting orbital localized on the phenazine motif. While no sensitivity to ground state electronic coupling between the individual molecular subunits was observed, the spectra are sensitive to structural distortions of the tpphz bridge. Here these results demonstrate N K-edge XANES as a local probe of electronic structure in large bridging ligand motifs, able to distinctly investigate the ligand-centered orbitals involved in metal-to-ligand and ligand-to-ligand electron transfer following light absorption.

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

Isotopic discrepancy in microscopic fusion of 16 O + 40,48 Ca

We present a detailed comparison of fusion in 16 O + 40,48 Ca within a fully microscopic time-dependent Hartree-Fock framework over sub-, near-, and above-barrier energies. While 16 O + 40 Ca fusion cross sections are reproduced within the available experimental energy range, calculations for the neutron-rich 16 O + 48 Ca system systematically overestimate the measured near-barrier cross sections. The discrepancy persists across the tested energy density functionals, indicating that it is not tied to a particular parametrization. Analysis of the dynamical interaction potentials reveals an outward shift and softening of the barrier region at the lowest collision energies. These results are then compared to a frozen density description where the barrier modifications are absent and the predicted cross sections are brought closer to experiment. Here, this demonstrates that the observed discrepancy is closely linked to dynamical surface rearrangement and isovector effects, and highlights the sensitivity of low-energy fusion to both the structure and dynamics of neutron-rich nuclei. This suggests that limitations in the isovector sector of conventional Skyrme EDFs in describing the structure of 48 Ca, as highlighted by CREX, are also manifested in the fusion dynamics of the 16 O + 48 Ca system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Spectroscopic and Photometric Monitoring of a Poorly Known Highly Luminous OH/IR Star: IRAS 18278+0931

We present the time-dependent properties of a poorly known OH/IR star, IRAS 18278+0931 (hereafter IRAS 18+09), toward the Ophiuchus constellation. We have carried out long-term optical/near-infrared photometric and spectroscopic observations to study the object. From optical R- and I-band light curves, the period of IRAS 18+09 is estimated to be 575 ± 30 days and the variability amplitudes range from ΔR ∼ 4.0 mag to ΔI ∼ 3.5 mag. From the standard period–luminosity relations, the distance (D) to the object, 4.0 ± 1.3 kpc, is estimated. Applying this distance in the radiative transfer model, the spectral energy distribution is constructed from multiwavelength photometric and IRAS-LRS spectral data, which provide the luminosity, optical depth, and gas mass-loss rate of the object to be 9600 ± 500 L {sub ⊙}, 9.1 ± 0.6 at 0.55 μm, and 1.0 × 10{sup −6} M {sub ⊙} yr{sup −1}, respectively. The current mass of the object is inferred to be in the range 1.0−1.5 M {sub ⊙} assuming solar metallicity. Notably, the temporal variation of atomic and molecular features (e.g., TiO, Na i, Ca i, CO, H{sub 2}O) over the pulsation cycle of the OH/IR star illustrates the sensitivity of the spectral features to the dynamical atmosphere as observed in pulsating AGB stars.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator

The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3–4 orders of magnitude fewer Hamiltonian evaluations and 2–3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2–3/1–2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ~1.3×10 9 parameters to obtain accurate dynamics.

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

X-ray absorption spectroscopy and theoretical investigations of the effect of extended ligands in potassium organic matter interaction

Potassium (K) is an essential nutrient for plant growth, and despite its abundance in soil, most of the K is structurally bound in minerals, limiting its bioavailability and making this soil K reservoir largely inaccessible to plants. Microbial biochemical weathering has been shown to be a promising pathway to sustainably increase plant available K. However, the mechanisms underpinning microbial K uptake, transformation, storage, and sharing are poorly resolved. Here, to better understand the controls on microbial K transformations, we performed K K-edge x-ray absorption near-edge structure (XANES) spectroscopy on K-organic salts, including acetate, citrate, nitrate, oxalate, and tartrate, which are frequently observed as low molecular weight organic acids secreted by soil microbes, as well as humic acid, which acts as a proxy for higher molecular weight organic acids. The organic salts display feature-rich K XANES spectra, each demonstrating numerous unique features spanning ~13 eV range across the absorption edge. In contrast, the spectra for humic acid have one broad, wide feature across the same energy range. We used a combination of time-dependent density functional theory and the Bethe–Salpeter equation based approach within the OCEAN code to simulate the experimental spectra for K-nitrate (KNO 3 ) and K-citrate [K 3 (C 6 H 5 O 7 )·H 2 O] to identify the electronic transitions that give rise to some of the outlying and unique spectral features in the organic salts. KNO 3 has both the lowest and highest lying energy features, and K 3 (C 6 H 5 O 7 )·H 2 O is produced by several soil microbes and is effective at mineral weathering. Our results analyze the K-organic salt bonding in detail to elucidate why the spectral shapes differ and indicate that the K K-edge XANES spectra are associated with the entire ligand despite similar first-shell bonding environments around the K center. The improved understanding of K bonding environments with organic ligands and their use for interpretation of the K-XANES spectra provides an important toolkit to understand how K is transformed by microbial processes and made bioavailable for plant uptake.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗