Automatic Clustering of Excited-State Trajectories: Application to Photoexcited Dynamics
Not Available
SEARCH · Search NASA
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.
Not Available
Multiconfiguration pair-density functional theory (MC-PDFT) is a computationally efficient method that computes the energies of electronic states in a state specific or state average framework via an on-top functional. However, MC-PDFT does not include state interaction among these states since the final energies do not come from the diagonalization of an effective model-space Hamiltonian. Recently, multistate extensions such as linearized PDFT (L-PDFT) have been developed to accurately model the potentials near conical intersections and avoided crossings. However, there has not been any systematic study evaluating their performance for predicting vertical excitations at the equilibrium geometry of a molecule, when the excited states are generally well separated. In this paper, we report the performance of L-PDFT on the extensive QUESTDB data set of vertical excitations using a database of automatically selected active spaces. We show that L-PDFT performs well on all these excitations and successfully reproduces the performance of MC-PDFT. These results further demonstrate the potential of L-PDFT, as its scaling is constant with the number of states included in the state-average manifold, whereas MC-PDFT scales linearly in this regard.
Hydropower units often rely on legacy excitation systems with fixed automatic voltage regulator (AVR) settings, but parameter drift, changing operating conditions, and evolving grid-service demands can degrade terminal-voltage regulation over time. This problem should be addressed because poor voltage regulation affects reactive-power support, system stability, and reliable plant operation, and replacing installed AVR and limiter logic is costly and disruptive in practice. To address this issue, this paper proposes a retrofit-friendly enhancement for hydropower generator excitation systems that combines adaptive tuning of the AVR middle lead–lag branch with a small proportional-integral correction for removing steady bias and slow drift. The adaptive laws are normalized, projection-bounded, and equipped with leakage, while adaptation is enabled only in the middle branch. Thus, the intent of the underexcitation and overexcitation limiters, as well as legacy current and rate limits, remains unchanged. Stability is analyzed under incremental performance assumptions, yielding boundedness and asymptotic tracking in the absence of saturation. The method is evaluated using real operational data from a utility-scale hydropower unit at Rocky Reach Dam, Unit C-8. Compared with the fixed AVR baseline, the proposed method reduces terminal-voltage mean squared error by about 70%, removes long periods of steady-state bias, and improves transient tracking while coexisting with limiter actions. These results show that substantial voltage-control improvement can be achieved on existing hydropower units with minimal integration effort and without replacing the installed AVR structure.
Aims. We present an improved database of temperature-dependent rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions. The database includes 231 transitions between the lower para-states of H 2 O and 210 transitions between its lower ortho-states (up toj= 7) and can be employed in cometary and planetary applications up to the temperature of 1000 K. Methods. We developed and applied a new general method that allows the generation of rate coefficients for excitation and quenching processes that automatically satisfy the principle of microscopic reversibility and also helps to cover the range of low collision energies by interpolation of cross sections between the process threshold and the computed data points. Results. We find that in the range of intermediate temperatures, 150
Abstract Background Estimation of the accuracy (quality) of protein structural models is important for both prediction and use of protein structural models. Deep learning methods have been used to integrate protein structure features to predict the quality of protein models. Inter-residue distances are key information for predicting protein’s tertiary structures and therefore have good potentials to predict the quality of protein structural models. However, few methods have been developed to fully take advantage of predicted inter-residue distance maps to estimate the accuracy of a single protein structural model. Result We developed an attentive 2D convolutional neural network (CNN) with channel-wise attention to take only a raw difference map between the inter-residue distance map calculated from a single protein model and the distance map predicted from the protein sequence as input to predict the quality of the model. The network comprises multiple convolutional layers, batch normalization layers, dense layers, and Squeeze-and-Excitation blocks with attention to automatically extract features relevant to protein model quality from the raw input without using any expert-curated features. We evaluated DISTEMA’s capability of selecting the best models for CASP13 targets in terms of ranking loss of GDT-TS score. The ranking loss of DISTEMA is 0.079, lower than several state-of-the-art single-model quality assessment methods. Conclusion This work demonstrates that using raw inter-residue distance information with deep learning can predict the quality of protein structural models reasonably well. DISTEMA is freely at https://github.com/jianlin-cheng/DISTEMA
The emergence of variational quantum applications has led to the development of automatic differentiation techniques in quantum computing. Existing work has formulated differentiable quantum programming with bounded loops, providing a framework for scalable gradient calculation by quantum means for training quantum variational applications. However, promising parameterized quantum applications, e.g., quantum walk and unitary implementation, cannot be trained in the existing framework due to the natural involvement of unbounded loops. To fill in the gap, we provide the first differentiable quantum programming framework with unbounded loops, including a newly designed differentiation rule, code transformation, and their correctness proof. Technically, we introduce a randomized estimator for derivatives to deal with the infinite sum in the differentiation of unbounded loops, whose applicability in classical and probabilistic programming is also discussed. We implement our framework with Python and Q# and demonstrate a reasonable sample efficiency. Through extensive case studies, we showcase an exciting application of our framework in automatically identifying close-to-optimal parameters for several parameterized quantum applications.
Apparatus and methods are disclosed for control systems for improving stability of electrical grids by temporarily reducing voltage output of electrical generators responsive to transient events on an electrical grid. In one example of the disclosed technology, a controller is coupled is to an automatic voltage regulator, which in turn adjusts excitation current of an electrical generator responsive to changes in frequency detected for the electrical grid. Reducing the output voltage temporarily allows for smaller generators to provide power to the microgrid. Methods for selecting parameters determining how the controller generates a regulation signal used to adjust the excitation current are further disclosed.
We show that linear superpositions of plane waves involving a single-valued, covariantly stable dispersion relation ω ( k ) always propagate outside the light cone unless ω ( k ) = a + b k . This implies that there is no notion of causality for individual dispersion relations since no mathematical condition on the function ω ( k ) (such as the front velocity or the asymptotic group velocity conditions) can serve as a sufficient condition for subluminal propagation in dispersive media. Instead, causality can only emerge from a careful cancellation that occurs when one superimposes all the excitation branches of a physical model. This happens automatically in local theories of matter that are covariantly stable. Hence, we find that the need for nonhydrodynamic modes in relativistic fluid mechanics is analogous to the need for antiparticles in relativistic quantum mechanics. Published by the American Physical Society 2024
Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.
Understanding how neurons interact across the brain to control animal behaviors is one of the central goals in neuroscience. Recent developments in fluorescent microscopy and genetically-encoded calcium indicators led to the establishment of whole-brain imaging methods in zebrafish, which record neural activity across a brain-wide volume with single-cell resolution. Pioneering studies of whole-brain imaging used custom light-sheet microscopes, and their operation relied on commercially developed and maintained software not available globally. Hence it has been challenging to disseminate and develop the technology in the research community. Here, we present PyZebrascope, an open-source Python platform designed for neural activity imaging in zebrafish using light-sheet microscopy. PyZebrascope has intuitive user interfaces and supports essential features for whole-brain imaging, such as two orthogonal excitation beams and eye damage prevention. Its camera module can handle image data throughput of up to 800 MB/s from camera acquisition to file writing while maintaining stable CPU and memory usage. Its modular architecture allows the inclusion of advanced algorithms for microscope control and image processing. As a proof of concept, we implemented a novel automatic algorithm for maximizing the image resolution in the brain by precisely aligning the excitation beams to the image focal plane. PyZebrascope enables whole-brain neural activity imaging in fish behaving in a virtual reality environment. Thus, PyZebrascope will help disseminate and develop light-sheet microscopy techniques in the neuroscience community and advance our understanding of whole-brain neural dynamics during animal behaviors.
We introduce an extension to the PySCF package, which makes it automatically differentiable. The implementation strategy is discussed, and example applications are presented to demonstrate the automatic differentiation framework for quantum chemistry methodology development. These include orbital optimization, properties, excited-state energies, and derivative couplings, at the mean-field level and beyond, in both molecules and solids. We also discuss some current limitations and directions for future work.
Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.
Abstract Ferroelectrics are being increasingly called upon for electronic devices in extreme environments. Device performance and energy efficiency is highly correlated to clock frequency, operational voltage, and resistive loss. To increase performance it is common to engineer ferroelectric domain structure with highly‐correlated electrical and elastic coupling that elicit fast and efficient collective switching. Designing domain structures with advantageous properties is difficult because the mechanisms involved in collective switching are poorly understood and difficult to investigate. Collective switching is a hierarchical process where the nano‐ and mesoscale responses control the macroscopic properties. Using chemical solution synthesis, epitaxially nearly‐relaxed (100) BaTiO 3 films are synthesized. Thermal strain induces a strongly‐correlated domain structure with alternating domains of polarization along the [010] and [001] in‐plane axes and 90° domain walls along the [011] or [01] directions. Simultaneous capacitance–voltage measurements and band‐excitation piezoresponse force microscopy revealed strong collective switching behavior. Using a deep convolutional autoencoder, hierarchical switching is automatically tracked and the switching pathway is identified. The collective switching velocities are calculated to be ≈500 cm s −1 at 5 V (7 kV cm −1 ), orders‐of‐magnitude faster than expected. These combinations of properties are promising for high‐speed tunable dielectrics and low‐voltage ferroelectric memories and logic.
Apparatus (10) and methods for making simultaneous measurements of composition (water-cut), fluid flow, and sound attenuation in a multiphase fluid flowing (15) through a pipe (12) in real-time, using the same apparatus (10) are described. Additionally, the apparatus (10) provides real-time pipe wall thickness monitoring for observing pipe corrosion or internal deposition. Knowledge of wall thickness is necessary to correct for water-cut (oil-water composition) automatically by adjusting the liquid path length internal to the pipe (spool). The use of short duration frequency chirp excitation signals (24) enables the apparatus to provide information that can be used to extract multiple levels of information from the same measurement in multiphase fluids including the presence of a significant quantity of gas (˜60% gas volume fraction) in different flow regimes. Besides measuring steady flow, this device is useful for measurements during fast changing flows, such as for a rod-pumped well. Measurements up to about 1000 times a second can reliably be made.
This paper investigates the potential of using brushless excitation (BLE) for not only riding through the fault but also to damp inter-area power oscillation with the help of a synchronous generator (SG) used in Type 5 Wind Power Plant (WPP). In BLE, an auxiliary synchronous generator (ASG) behaving like an exciter is coupled and driven by the rotor of the main SG. The BLE's field current's ASG is fed by two separate loops of the automatic voltage regulator (AVR) and power system stabilizer (PSS). The AVR system implements a control loop to regulate the generator terminal voltage, VT. For the PSS, the kinetic energy of the wind turbine is utilized according to the estimated rotational speed of the synchronous generator shaft of the SG. It may mitigate the necessity of any curtailment of active power for damping. The effectiveness of the proposed control scheme is verified with a three-phase short circuit fault in a two-area power system.
This paper investigates the potential of using brushless excitation (BLE) for not only riding through the fault but also to damp inter-area power oscillation with the help of a synchronous generator (SG) used in Type 5 Wind Power Plant (WPP). In BLE, an auxiliary synchronous generator (ASG) behaving like an exciter is coupled and driven by the rotor of the main SG. The BLE's field current's ASG is fed by two separate loops of the automatic voltage regulator (AVR) and Power System Stabilizer (PSS). The AVR system implements a control loop to regulate the generator terminal voltage, VT. For the PSS, the kinetic energy of the wind turbine is utilized according to the estimated rotational speed of the synchronous generator shaft of the SG. It mitigates the necessity of any curtailment of active power for damping. The effectiveness of the proposed control scheme is verified with a three-phase short circuit fault in a two-area power system.
Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.
The isotope group at Jefferson Lab is carrying out R&D for producing medically interesting radioisotopes, especially those with theranostic (therapeutic and diagnostic) attributes. Here the search for viable production mechanisms has been expanded to multi-step reactions where a daughter is produced from the target and decays into a medically interesting granddaughter radioisotope. It is difficult to find efficient production routes when investigating both the initial excitation reaction as well as the decay routes leading to medically interesting isotopes. The overall goal of this project is to create a structured code in Python to find these decay routes by automatically exploring the large number of isotopes and their possible decay modes. The program structure is described, and preliminary results are presented.