Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time-dependent Schrödinger equation”

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.

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. Here, we use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear Optimal Control of Electron Dynamics Within Hartree-Fock Theory

Consider the problem of determining the optimal applied electric field to drive a molecule from an initial state to a desired target state. For even moderately sized molecules, solving this problem directly using the exact equations of motion—the time-dependent Schrödinger equation (TDSE)—is numerically intractable. Here, we present a solution of this problem within time-dependent Hartree-Fock (TDHF) theory, a mean field approximation of the TDSE. Optimality is defined in terms of minimizing the total control effort while maximizing the overlap between desired and achieved target states. We frame this problem as an optimization problem constrained by the nonlinear TDHF equations; we solve it using trust region optimization with gradients computed via a custom-built adjoint state method. For three molecular systems, we show that with very small neural network parametrizations of the control, our method yields solutions that achieve desired targets within acceptable constraints and tolerances.

97 MATHEMATICS AND COMPUTING↗

Control simulations of many-body quantum systems by a synergism of discrete real-time learning and optimal control theory

We present a self-consistent algorithm for optimal control simulations of many-body quantum systems. The algorithm features a two-step synergism that combines discrete real-time machine learning (DRTL) with Quantum Optimal Control Theory (QOCT) using the time-dependent Schrödinger equation. Specifically, in step (1), DRTL is employed to identify a compact working space (i.e., the important portion of the Hilbert space) for the time evolution of the many-body quantum system in the presence of a control field (i.e., the initial or previously updated field), and in step (2), QOCT utilizes the DRTL-determined working space to find a newly updated control field for a chosen objective. Steps 1 and 2 are iterated until a self-consistent control objective value is reached such that the resulting optimal control field yields the same targeted objective value when the corresponding working space is systematically enlarged. Furthermore, to demonstrate this two-step self-consistent DRTL-QOCT synergistic algorithm, we perform optimal control simulations of strongly interacting 1D as well as 2D Heisenberg spin systems. In both scenarios, only a single spin (at the left end site for 1D and the upper left corner site for 2D) is driven by the time-dependent control fields to create an excitation at the opposite site as the target. It is found that, starting from all spin-down zero excitation states, the synergistic method is able to identify working spaces and convergence of the desired controlled dynamics with just a few iterations of the overall algorithm. In the cases studied, the dimensionality of the working space scales only quasi-linearly with the number of spins.

Artificial neural networks↗

Ultrafast x-ray imaging of coherently controlled molecular dynamics in real space and time

Coherent control aims to manipulate chemical processes on the latent length and timescales of atoms and bonds, scales that are intrinsic to molecular dynamics but not directly resolved by most experimental probes. As a result, existing coherent-control experiments, which overwhelmingly rely on spectroscopic observables, leave a crucial blind spot: the direct, real-space recovery of all atomic and molecular rearrangements in response to coherently controlled excitations. Here, we overcome this limitation by integrating a Tannor-Kosloff-Rice pump-control-probe scheme with ultrafast X-ray scattering to capture snapshots of the wavepacket motion in a benchmark molecular system. We demonstrate this by photoexciting diatomic iodine vapor with a visible pump pulse, selectively steering the wavepacket toward ground-state recombination or dissociative pathways with a time-delayed near-IR control pulse, and recording the dynamics with angstrom and femtosecond precision with an ultrashort hard X-ray probe pulse. By comparing these structural observations with numerical solutions of the time-dependent Schrödinger equation, we reveal how coherent control actively reshapes the molecular charge density distribution. Our results pave the way for leveraging structural feedback as a control handle and provide a fundamental microscopic visualization of quantum decoherence and energy redistribution at the atomic level.

Hopper, Thomas R. [SLAC National Accelerator Labor↗

Attosecond Optical Orientation

Circularly polarized light offers opportunities to probe symmetry-dependent properties of matter such as chirality and spin. Circular-dichroic measurements typically require further intrinsic or extrinsic breaking of symmetry by, e.g., enantiomeric excess, orientation, magnetic fields, or direction-sensitive detectors. Here we introduce circular-dichroic attosecond transient absorption spectroscopy by leveraging the angular momentum of two circularly polarized pulses, both pump and probe, in an isotropic medium, optically orienting the angular momentum of excited states on an attosecond timescale. We investigate a circular-dichroic measurement of the attosecond transient absorption of He Rydberg states. By limiting the allowed pathways via dipole selection rules for co- and counterrotating circularly polarized near-infrared and extreme ultraviolet (XUV) pulses, different spectral reshapings of the XUV transient absorption due to the AC Stark effect are observed. Furthermore, paired with time-dependent Schrödinger equation calculations, the results show the role of selection and propensity rules and open up new opportunities to study coupling pathways of excited states as well as spin-dependent dynamics in atoms and beyond via attosecond optical orientation.

Angular momentum of light↗

Photoelectron spectra in circularly and elliptically polarized laser pulses

We present results of numerical simulations of the time-dependent Schrödinger equation and theoretical analysis concerning the interaction of a rare gas atom with circularly and elliptically polarized laser pulses. In agreement with recent observation in circularly polarized fields, the photoelectron energy spectra for counterrotating electrons are peaked at lower kinetic energy than those for corotating electrons. We show that this difference can be interpreted as being due to the additional pathways to ionization that are available for the counterrotating electrons only. Furthermore, our results show, in agreement with earlier work, that the offset angle by which the emission of electrons is rotated in an elliptically polarized field increases with each successive above-threshold ionization (ATI) order and is larger for the emission from counterrotating states as compared to that from corotating states. As a result, a simple model based on the interference and relative phase difference between just three continuum states provides remarkable agreement and once again emphasizes the importance of the additional ionization pathways available for the counterrotating electrons.

74 ATOMIC AND MOLECULAR PHYSICS↗

Theory of resonant x-ray scattering with ultrafast intense pulses

Here, we present a time-dependent Schrödinger equation approach within a nonrelativistic quantum electrodynamics framework to investigate resonant x-ray scattering driven by intense x-ray pulses. This method enables us to explore how coherent x-ray electron dynamics influence scattering signals from Ne + . We account for both resonance fluorescence and elastic scattering channels, while also considering competing photoionization and inner-shell decay processes. By computing the angular distribution and energy spectrum of scattered photons, we uncover interference effects between elastic scattering and resonance fluorescence pathways. Notably, this interference results in a small asymmetry in the energy spectrum. We discuss the experimental potential for detecting signatures of interference. Our findings demonstrate that the x-ray Rabi dynamics can be used to control scattering responses and provide insights into interference mechanisms and scattering efficiency in high-intensity x-ray regimes.

Venkatesh, Akilesh [Argonne National Laboratory (A↗

Local enhancement in transient absorption spectroscopy by gating the resonance in the time domain

The concept of resonant perfect absorption, enabled by the combined action of pulse propagation and an auxiliary gate pulse, was recently proposed and demonstrated in a group of two-level systems [Y. He , ]. Here we exploit this method in a more realistic scenario by solving the coupled time-dependent Schrödinger equation and the Maxwell wave equation in helium. Through emptying the population of the 1 s 2 p excited state after its excitation, we explore the evolution of the spectral profile with time delay and propagation distance and link the observations to the controlled interference between the original field and the gated new field. We find that resonant absorption for higher-lying states can also be strongly enhanced, in spite of the congestion of multiple resonances and the presence of complex laser-induced couplings. Our results show that interferometric control of absorption using intense laser fields can be applied selectively in both the temporal and spatial domains. Published by the American Physical Society 2024

He, Yu (ORCID:0000000152120176)↗

Wave-packet manipulation of He Rydberg states by a seeded free-electron laser

We report a two-dimensional pump-control-probe spectroscopic study of the dynamics of singly excited He Rydberg-state wave packets with a seeded extreme ultraviolet (XUV) free-electron laser (FEL) source. A pair of coherent XUV pulses, defined by their coarse time separation and relative phase, created and manipulated the wave packets. The He atoms were postionized by infrared (IR) pulses, and the ion yield was measured as a function of XUV phase and IR arrival time. We tagged and sorted the relative phase of the XUV pulse pair on a single-shot basis by fitting each FEL spectrum with a suitable function that accounts for nonidealities of the XUV pulse pairs, associated with the seeding process; more generally, the fit returns the time-dependent electric field of the FEL spectra. The experimental two-dimensional maps of ion yields, measured as a function of IR (probe) delay and of XUV (pump-control) phase, were compared with the solution of the first-order time-dependent Schrödinger equation for this field. Despite the fact that the experimental conditions imply strong excitation, beyond the approximations of first-order perturbation theory, the simulated map satisfactorily reproduces the experimental one for temporally well-separated pulses. We show that by selecting data at appropriate values of pump-control phase, we enhance or suppress the amplitude of chosen wave-packet components consisting of two or more Rydberg states. When the temporal overlap of the pulse pair cannot be neglected, the phase reconstruction is underdetermined, and we provide a simplified comparison between data and simulations. Published by the American Physical Society 2024

Dumergue, Mathieu (ORCID:0000000236984406)↗

Resonant propagation of extreme-ultraviolet pulses through strongly driven high-density media

We show that by combining strong-field dressing and resonant propagation of XUV pulses, the transition of absorption lines from their natural Lorentzian profiles through Fano and complex multipeak shapes all the way back to broadened near-Lorentzian profiles can be achieved, in the limit of optically thick samples. The final stage of this spectral modification can be understood in terms of a significant temporal stretching and delay of the resonant XUV pulse as it propagates through the dense gas, which alters the ultrafast absorption that is modified by the time-synchronized few-femtosecond laser pulse. We first demonstrate this concept in numerical calculations, both using a model system and through a fully coupled solution of the time-dependent Schrödinger equation and the Maxwell wave equation. The applicability and generality of the underlying mechanism is then illustrated in proof-of-principle attosecond transient absorption measurements in a helium gas with a widely varying atomic density. These results provide insights into the interaction of ultrashort laser pulses with dense media and its coherent control. Published by the American Physical Society 2025

He, Yu (ORCID:0000000152120176)↗

New correlated numerical methods for attosecond molecular single and double ionization

This project resulted in the development of ASTRA, a new molecular ionization code capable of describing single and double ionization processes in polyatomic molecules with correlated electrons and time-dependent fields. Key achievements include: (i) the development of the PRISM hybrid-integral library enabling core-ionization and high-angular-momentum convergence; (ii) implementation of time-dependent Schrödinger equation solvers for pump–probe spectroscopies; (iii) validation against experimental and theoretical benchmarks for atoms and molecules; and (iv) extension to two-photon double ionization using the finite-pulse virtual-sequential model. These developments establish ASTRA as a versatile platform for attosecond molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Comments on the paper “Eliminating beam-induced depolarizing effects in the hydrogen jet target for high-precision proton beam polarimetry at the Electron-Ion Collider” (Part I)

A critical review of the methodology used in F. Rathmann et al., Phys. Rev. Accel. Beams 29, 021001 (2026), to evaluate beam-induced depolarization of the Atomic Polarized Hydrogen Gas Jet (HJET) target at the Electron–Ion Collider (EIC) is presented. It is shown that several key assumptions underlying that analysis—including the introduction of a photon emission threshold, the application of Fermi’s Golden Rule to coherent hyperfine transitions, the interpretation of power broadening as a physical linewidth increase, and the treatment of spatial magnetic fields—are either incorrect or internally inconsistent. As a consequence, the predicted large depolarization effects are demonstrated to be artifacts of the adopted methodology rather than genuine physical phenomena. A consistent quantum-mechanical treatment based on the time-dependent Schrödinger equation shows that beam-induced depolarization probabilities at the EIC are negligibly small.

43 PARTICLE ACCELERATORS↗

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗

Conformal duality of the nonlinear Schrödinger equation: Theory and applications to parameter estimation

The nonlinear Schrödinger equation (NLSE) in one spatial dimension has stationary solutions similar to those of the linear Schrödinger equation (LSE) as well as more exotic solutions such as solitary waves and quantum droplets. Here, we present a newly discovered conformal duality which unifies the stationary and time-dependent traveling-wave solutions of the one-dimensional cubic-quintic NLSE, the cubic NLSE and LSE. Any two systems that are classified by the same single number called the cross ratio are related by this symmetry. Notably, the conformal duality can also be adapted in Newtonian mechanics and serves as a powerful tool for investigating physical systems that otherwise cannot be directly accessed in experiments. Further, we show that the conformal symmetry is a valuable resource to substantially improve NLSE parameter estimation from noisy empirical data by introducing an optimization afterburner. The new method therefore has far reaching practical applications for nonlinear physical systems. Published by the American Physical Society 2025

Reinhardt, David B. (ORCID:0009000409812838)↗