Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “excited state methods”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S 2 /S 3 (D 2d ) → S 1 (D 8h ) → S o (D 4h ) → S o (D 2d ) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. Here, the obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

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

Spectral Signatures of Ultrafast Excited-State Intramolecular Proton Transfer from Computational Multi-edge Transient X-ray Absorption Spectroscopy

Excited state intramolecular proton transfer (ESIPT) is a fundamental chemical process with several applications. Ultrafast ESIPT involves coupled electronic and atomic motions and has been primarily studied using femtosecond optical spectroscopy. X-ray spectroscopy is particularly useful because it is element-specic and enables direct, individual probes of the proton donating and accepting atoms. Here, we report a computational study to resolve the ESIPT in 10-Hydroxybenzo[h]quinoline (HBQ), an intramolecularly hydrogen bonded compound, using time-dependent density functional theory combined with ab initio molecular dynamics (AIMD) and time-resolved X-ray absorption spectroscopy (XAS) computations to track the ultrafast excited-state dynamics. Overall, our results reveal clear X-ray spectral signatures of coupled electronic and atomic motions during and following ESIPT at the oxygen and nitrogen K-edge paving the way for future experiments at X-ray free electron lasers.

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

Excited-state uncertainties in lattice-QCD calculations of multi-hadron systems

Excited-state effects lead to hard-to-quantify systematic uncertainties in lattice quantum chromodynamics (LQCD) spectroscopy calculations when computationally accessible imaginary times are smaller than inverse excitation gaps, as often arises for multi-hadron systems with signal-to-noise problems. Lanczos residual bounds address this by providing two-sided constraints on energies that do not require assumptions beyond Hermiticity, but often give very conservative systematic uncertainty estimates. Here, a more-constraining set of gap bounds is introduced for hadron spectroscopy. These bounds provide tighter constraints whose validity requires an explicit assumption about an energy gap. Exactly solvable lattice field theory correlators are used to test the utility of residual and gap bounds at finite and infinite statistics. Two-sided bounds and other analysis methods are then applied to a high-statistics LQCD calculation of nucleon-nucleon scattering at $m_π\sim 800$ MeV. Generalized eigenvalue problem (GEVP) and Lanczos energy estimators are compatible when applied to the same correlator data, but analyses including different interpolating operators show statistically significant inconsistencies. However, two-sided bounds from all operators are consistent. Under the assumption that the number of energy levels below $NΔ$ and $ΔΔ$ thresholds is the same as for non-interacting nucleons, gap bounds are sufficient to constrain nucleon-nucleon scattering amplitudes at phenomenologically relevant precision. Lanczos methods further reveal that energy-eigenstate estimates from previously studied asymmetric correlators have not converged over accessible imaginary times. Nevertheless, data-driven examples demonstrate why assumptions are required to draw conclusions about the natures of two-nucleon ground states at these masses.

Detmold, William [MIT, Cambridge, CTP]↗
DOE OSTI2022

Colloidal Nanoribbons: From Infrared to Visible

Using the cation-exchange method, colloidal PbS nanoribbons are converted completely into CdS nanoribbons. This process expands the emission spectrum of the nanoribbons from infrared to visible. Here, the morphology of nanoribbons remains the same after cation exchange, but the crystal structure changes from rock salt to zincblende. CdS nanoribbons exhibit blue band-edge photoluminescence under ultraviolet-light excitation. Cathodoluminescence spectroscopy of the CdS nanoribbons shows multicolor (blue, green, and red) emissions. Further time-resolved photoluminescence spectroscopy studies show that the lifetime of the midgap states is more than 2 orders of magnitude longer than that of the band-edge states.

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

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Establishing computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: Here, we propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: 44,46,48 Ti, 48,50 Cr, 52 Fe, 60,62,64 Zn, 66,68 Ge, 68 Se, and 108,110 Xe. We predict low-lying states of 112–114 Ba and 116–120 Ce, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal

Highlights: • We look at the role played by generic scale invariance (GSI) in the phase transition between a U(1) spin liquid and a Fermi liquid. • We identify dangerously irrelevant operators in the renormalization group transformation, involved with the GSI. • We highlight novel mean-field dynamics of this Mott quantum criticality, resulting from the GSI. We investigate the role of generic scale invariance in a Mott transition from a U(1) spin-liquid insulator to a Landau Fermi-liquid metal, where there exist massless degrees of freedom in addition to quantum critical fluctuations. Here, the Mott quantum criticality is described by critical charge fluctuations, and additional gapless excitations are U(1) gauge-field fluctuations coupled to a spinon Fermi surface in the spin-liquid state, which turn out to play a central role in the Mott transition. An interesting feature of this problem is that the scaling dimension of effective leading local interactions between critical charge fluctuations differs from that of the coupling constant between U(1) gauge fields and matter-field fluctuations in the presence of a Fermi surface. As a result, there appear dangerously irrelevant operators, which can cause conceptual difficulty in the implementation of renormalization group (RG) transformations. Indeed, we find that the curvature term along the angular direction of the spinon Fermi surface is dangerously irrelevant at this spin-liquid Mott quantum criticality, responsible for divergence of the self-energy correction term in U(1) gauge-field fluctuations. Performing the RG analysis in the one-loop level based on the dimensional regularization method, we reveal that such extremely overdamped dynamics of U(1) gauge-field fluctuations, which originates from the emergent one-dimensional dynamics of spinons, does not cause any renormalization effects to the effective dynamics of both critical charge fluctuations and spinon excitations. However, it turns out that the coupling between U(1) gauge-field fluctuations and both matter-field excitations still persists at this Mott transition, which results in novel mean-field dynamics to explain the nature of the spin-liquid Mott quantum criticality. We discuss physical implications of effective one-dimensional spin dynamics and extremely overdamped gauge dynamics at the Mott quantum criticality.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigation of the isoscalar monopole response in the proton-rich nucleus 14 O

Deuteron inelastic scattering on 14 O was measured in inverse kinematics using an active-target time projection chamber and a magnetic spectrograph. The experimental technique enabled precise measurements of deuteron recoiling particles in coincidence with beam-like fragments detected in the spectrograph focal plane. The double differential cross section was reconstructed for scattering angles of 3–6 degrees and excitation energies up to 26 MeV. The monopole strength distribution was obtained from the data using a multipole decomposition analysis. The results were compared to quasiparticle random-phase approximation (QRPA) and generator coordinate method (GCM) calculations. The QRPA calculation accurately describes experimental data in the energy range of 13 to 26 MeV. GCM calculations assuming a 12 C (g . s .) + p + p cluster configuration were used to determine the 0 + strength in 14 O below 13 MeV. The monopole transition strength of these cluster states provides a good description of the experimental distribution in the 9–11 MeV region, while the 0$^{+}_{2}$ state accounts for only a small fraction of the experimental strength around 6 MeV.

Active target↗
DOE OSTI2026

Optimisation of the Kaplan hydropower system via PID 2 and digital twin

Here, this paper proposes a proportional–integral-double–derivative (PID 2 ) optimisation method for the Kaplan hydropower system by building a digital twin. The study first uses one multilayer perceptron (MLP) to model the hydroturbine dynamic and then adopts three connected MLPs to model the generator dynamic, both in an open-loop fashion. Inspired by stochastic distribution control (SDC) theory, we regard the training of the turbine's neural network model as a process control problem, and we propose minimising entropy loss to update the network parameters. The next step is to build the digital twin by connecting the neural network models with a PID 2 controller and a lead-lag exciter and run the whole model in a closed-loop fashion. After that, a binary search approach is applied to optimise the PID 2 parameters based on the obtained digital twin model. The simulation results show that the proposed method can reduce the mean square tracking error by more than 90%. Furthermore, the method is extended to jointly optimise the PID 2 controller and excitation system gains through multiobjective optimisation, leveraging Pareto frontier analysis to balance active power and voltage tracking performance. Simulation results confirm the effectiveness of the proposed method, achieving a 83.46% reduction in relative mean square error of active power, a 47.13% reduction in terminal voltage tracking error, and an 82.78% improvement in the overall scalarized objective.

Hydropower system↗
DOE OSTI2024

A molecular ground electronic state with an occupied 5g spinor—The superheavy (E125)F molecule

Fully relativistic calculations, primarily at the 4-component coupled-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory with the Dirac–Coulomb (DC) Hamiltonian, have been carried out for the superheavy (E125)F molecule using large Gaussian basis sets. The electronic ground state is determined to have an [Og]8s 2 5g 1 6f 3 configuration on E125 with an Ω = 6 ground state and an 8p electron largely donated to F. A Mulliken population analysis indicates that the ground state is mainly ionic with a partial charge of +0.79 on E125 and a single sigma bond involving the F 2p and E125 8p spinors. The occupied g spinor is not involved in the bonding. With the largest basis set used in this work, the (0 K) dissociation energy was calculated at the DC-CCSD(T) level of theory to be 7.02 eV. Analogous calculations were also carried out for the E125 atom, both the neutral and its cation. The lowest energy electron configuration of E125 + , [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{3}_{5/2}$ with a J = 6 ground state, was found to be similar to that in (E125)F, while the neutral E125 atom has an [Og]$8s$$^{2}_{1/2}$$5g$$^{1}_{7/2}$$6f$$^{2}_{5/2}$$7d$$^{1}_{3/2}$$8p$$^{1}_{1/2}$ ground state electron configuration with a J = 17/2 ground state. The ionization energy (IE) of E125 is reported for the first time and is calculated to be 4.70 eV at the DC-CCSD(T) level of theory. Non-relativistic calculations were also carried out on the E125 atom and the (E125)F molecule. Here, the non-relativistic ground state of the E125 atom was calculated to have a 5g 5 ground state with an IE of just 3.4 eV. The net effect of relativity on (E125)F is to stabilize its bonding.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Source of Bright Near-Infrared Luminescence in Gold Nanoclusters

Gold nanoclusters with near-infrared (NIR) photoluminescence (PL) have great potential as sensing and imaging materials in biomedical and bioimaging applications. In this work, Au 21 (S-Adm) 15 and Au 38 S 2 (S-Adm) 20 are used to unravel the underlying mechanisms for the improved quantum yields (QY), large Stokes shifts and long PL lifetimes in gold nanoclusters. Both nanoclusters show decent PL QY. In particular, the Au 38 S 2 (S-Adm) 20 nanocluster shows a bright NIR PL at 900 nm with QY up to 15% in normal solvents (such as toluene) at ambient conditions. The relatively lower QY for Au 21 (S-Adm) 15 (4%) compared to Au 38 S 2 (S-Adm) 20 is attributed to the lowest-lying excited state being symmetry-disallowed, as evidenced by the pressure-dependent anti-spectral shift of the absorption spectra compared to PL. Yet, Au 21 (S-Adm) 15 maintains some emissive properties due to a nearby symmetry-allowed excited state. Furthermore, our results show that suppression of non-radiative decay due to the surface “lock rings” which encircle the Au kernel and the surface “lock atoms” which bridge the fundamental Au-kernel units (e.g., tetrahedra, icosahedra, etc.) is the key to obtain high QYs in gold nanoclusters. Here, the complicated excited-state processes and the small absorption coefficient of the band-edge transition lead to the large Stokes shifts and the long PL lifetimes that are widely observed in gold nanoclusters.

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

Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language

The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.

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

Search for an s -wave resonance in Li 7 just above the proton-decay threshold

Near-threshold resonances play an important role in nucleosynthesis and applied nuclear science. The study of nuclei removed from stability has greatly extended the list of resonances close to decay thresholds. The no core shell model with continuum (NCSMC) recently predicted an s-wave resonance just above the proton-decay threshold of 7 Li at an excitation energy of 10.2 MeV. This potential case of a near-threshold resonance is dependent on the quantum mechanics of the p + 6 He fragments extended into the continuum. Further, the 6 He(d,n) 7 Li* reaction was employed in an attempt to populate this resonance and search for its proton decay via the invariant-mass method. No evidence of this resonance was found. However, the data collected in this search led to better constraints on the energy and width of the isobaric analog state (IAS) with (J π ,T) = (3/2 - , 3/2) and revealed a new weak resonance just above the IAS in energy, predicted by the NCSMC as a (J π , T ) = (3/2 – , 1/2) state and potentially a part of the antianalog structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Roper Resonance Structure and Exploration of Emergent Hadron Mass from CLAS Electroproduction Data

The N(1440)1/2 + nucleon resonance, first identified in 1964 by Roper and collaborators in analyses of πN hadroproduction data, has continued to provide pivotal insights that serve to advance our understanding of nucleon excited states. In this contribution, we present results from studies of the structure of the Roper resonance based on exclusive πN and π + π − p electro production data measured with the CLAS detector at Jefferson Lab. These analyses have revealed the Roper resonance as a complex interplay between an inner core of three dressed quarks and an external meson–baryon cloud. Analyses of the CLAS results on the evolution of the Roper resonance electroexcitation amplitudes with photon virtuality Q 2 , within the framework of the Continuum Schwinger Method, have conclusively demonstrated the capability to gain insight into the strong interaction dynamics responsible for generating more than 98% of hadron mass. Further extension of such studies to higher Q 2 , through experiments currently underway with the CLAS12 detector and in the future with a potential CEBAF energy upgrade to 22 GeV, offers the only foreseeable opportunity to explore the full range of distances where the dominant portion of hadron mass and resonance structure emerge.

Mokeev, Victor I. [Thomas Jefferson National Accel↗
DOE OSTI2022

Hadron spectroscopy of pseudoscalar and vector meson photoproduction using linearly polarised photons at CLAS

The field of hadron spectroscopy is informed by experiments which study excited states of the nucleon by detecting the decay products following excitation with a hadronic, leptonic or electromagnetic probe. Experiments using a polarised beam of photons provide an important contribution to the world data. This thesis describes the analysis of data from one such experiment, the CEBAF Large Acceptance Spectrometer at Jefferson Laboratory, Virginia, USA. In the experiment, an incident beam of electrons passes through a diamond radiator producing a linearly polarised photon beam of energy from 1.1 to 2.1 GeV on a liquid hydrogen target. The analysis method and results for two reaction channels are presented. In the first analysis, the polarisation observables {?, P, T, Ox, Oz} for the reaction ??? p ? K0 S ?+ are extracted simultaneously using likelihood sampling with Markov-Chain Monte-Carlo. The results are presented as angular distributions (cos ?K0S) for bins in photon energy E? for a total of 21 bins. The values extracted for T , Ox and Oz are a first measurement of these quantities for this reaction and the beam asymmetry (?) measurement supplements the one previous measurement for this quantity and extends the energy range over which it has been extracted. The measurement of the recoil polarisation, P, is consistent with previous measurements. In a second analysis, the spin density matrix elements (SDMEs) for ? p ? p ? are extracted using similar analysis techniques. The results are presented as angular distributions (cos ?? ) for bins in photon energy E? for a total of 24 bins. This work is the first to extract all nine SDMEs for the full angular range of the ? meson production. A comparison is made to the predictions of the Vector Meson Dominance (VMD) model and the results provide evidence for contributing processes other than VMD, particularly at more backward angles of the ? meson production. The data also indicate non-helicity conserving process in both the s-channel and the t-channel. The results from the first analysis have been passed to the Juelich-Bonn theory group and the effect of including the new data within their dynamical coupled channel model is described. Several nucleon resonances are affected and the impact on the pole positions, widths and photon decay amplitudes of the affected resonances is described.

Clark, Louise↗
DOE OSTI2021

Wavelength conversion through plasmon-coupled surface states

Surface states generally degrade semiconductor device performance by raising the charge injection barrier height, introducing localized trap states, inducing surface leakage current, and altering the electric potential. We show that the giant built-in electric field created by the surface states can be harnessed to enable passive wavelength conversion without utilizing any nonlinear optical phenomena. Photo-excited surface plasmons are coupled to the surface states to generate an electron gas, which is routed to a nanoantenna array through the giant electric field created by the surface states. The induced current on the nanoantennas, which contains mixing product of different optical frequency components, generates radiation at the beat frequencies of the incident photons. We utilize the functionalities of plasmon-coupled surface states to demonstrate passive wavelength conversion of nanojoule optical pulses at a 1550 nm center wavelength to terahertz regime with efficiencies that exceed nonlinear optical methods by 4-orders of magnitude.

42 ENGINEERING↗
DOE OSTI2023

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

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

Machine-Learning X-Ray Absorption Spectra to Quantitative Accuracy

Simulations of excited state properties, such as spectral functions, are often computationally expensive and therefore not suitable for high-throughput modeling. As a proof of principle, here we demonstrate that graph-based neural networks can be used to predict the x-ray absorption near-edge structure spectra of molecules to quantitative accuracy. Specifically, the predicted spectra reproduce nearly all prominent peaks, with 90% of the predicted peak locations within 1 eV of the ground truth. Besides its own utility in spectral analysis and structure inference, our method can be combined with structure search algorithms to enable high-throughput spectrum sampling of the vast material configuration space, which opens up new pathways to material design and discovery.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Locating the Source of Forced Oscillations in Transmission Power Grids

A forced oscillation event in power grids refers to a state where malfunctioning or abnormally operating equipment causes persisting periodic disturbances in the system. While power grids are designed to damp most perturbations during standard operations, some of them can excite normal modes of the system and cause significant energy transfers across the system, creating large oscillations thousands of miles away from the source. Localization of the source of such disturbances remains an outstanding challenge due to a limited knowledge of the system parameters outside of the zone of responsibility of system operators. Here, we propose a new method for locating the source of forced oscillations that addresses this challenge by performing a simultaneous dynamic model identification using a principled maximum likelihood approach. We illustrate the validity of the algorithm on a variety of examples where forcing leads to resonance conditions in the system dynamics. Our results establish that an accurate knowledge of system parameters is not required for a successful inference of the source and frequency of a forced oscillation. We anticipate that our method will find a broader application in general dynamical systems that can be well described by their linearized dynamics over short periods of time.

24 POWER TRANSMISSION AND DISTRIBUTION↗