Time evolution as an optimization problem: The hydrogen atom in strong laser fields in a basis of time-dependent Gaussian wave packets

Simon Elias Schrader,Håkon Emil Kristiansen,Thomas Bondo Pedersen,Simen Kvaal

2024-04-11

Abstract:Recent advances in attosecond science have made it increasingly important to develop stable, reliable and accurate algorithms and methods to model the time evolution of atoms and molecules in intense laser fields. A key process in attosecond science is high-harmonic generation, which is challenging to model with fixed Gaussian basis sets, as it produces high-energy electrons, with a resulting rapidly varying and highly oscillatory wave function that extends over dozens of ångström. Recently, Rothe's method, where time evolution is rephrased as an optimization problem, has been applied to the one-dimensional Schrödinger equation. Here, we apply Rothe's method to the hydrogen wave function and demonstrate that complex-valued Gaussian wave packets with time-dependent width, center, and momentum parameters are able to reproduce spectra obtained from essentially exact grid calculations for high-harmonic generation with only 50-181 Gaussians for field strengths up to $5\times 10^{14}$W/cm$^2$. This paves the way for the inclusion of continuum contributions into real-time, time-dependent electronic-structure theory with Gaussian basis sets for strong fields, and eventually accurate simulations of the time evolution of molecules without the Born-Oppenheimer approximation.

Chemical Physics,Atomic Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the time - evolution problem of hydrogen atoms in a strong laser field, especially the challenges in the high - harmonic - generation (HHG) process. Specifically, the paper focuses on how to accurately simulate the electron - dynamics process of hydrogen atoms using a time - dependent Gaussian wave - packet basis set under the action of a strong laser field. Traditional methods such as the fixed Gaussian basis set face difficulties when dealing with high - energy electrons because these methods cannot effectively describe rapidly changing and highly oscillating wave functions. To solve these problems, the paper adopts Rothe's method, reformulating the time - evolution problem as an optimization problem, and achieving an accurate description of the wave function by adjusting the width, central position, and momentum parameters of the Gaussian wave packets. This method can reproduce the HHG spectra that match the exact grid - calculation results with a relatively small number of Gaussian wave packets (50 - 181), and is applicable to cases where the field strength is as high as \(5\times10^{14}\,\text{W/cm}^2\). This paves the way for incorporating continuous - state contributions into real - time, time - dependent electron - structure theory under strong - field conditions, and ultimately enables the accurate simulation of molecular time - evolution without using the Born - Oppenheimer approximation.