Ultrashort laser pulses on the attosecond timescale are typically achieved via high-order harmonic generation (HHG), a nonlinear process in which atoms interact with intense light fields to emit a broad spectrum of harmonics. HHG is commonly described in terms of a `quantum orbits' model based on several interfering electron trajectories, thereby incorporating both quantum-mechanical effects and an intuitive picture of classical dynamics. By tuning the parameters of the driving laser field, the interplay between these trajectories can be controlled, shaping the emitted light. Mathematically, this model expresses the harmonic response as a highly oscillatory integral. Applying saddle-point methods to this integral allows it to be decomposed into contributions from distinct saddle points of the semi-classical action, thereby linking quantum dynamics to classical trajectories. However, a general framework for applying these methods across arbitrary parameters and laser configurations has been missing.
In this thesis, we introduce Picard--Lefschetz theory and develop practical numerical methods for its application. These enable the evaluation of oscillatory integrals and identification of contributions from individual critical points.
We apply these techniques to strong-field ionisation and HHG, focusing on caustics -- enhancement features where trajectories coalesce and standard approximations fail. Our methods remain valid in these regions, allowing systematic analysis of parameter regimes and revealing previously inaccessible features. This work improves the understanding and control of ultrafast light--matter interactions.
Comments
PhD thesis, first submitted version (31 January 2026)
Ultrashort laser pulses on the attosecond timescale are typically achieved via high-order harmonic generation (HHG), a nonlinear process in which atoms interact with intense light fields to emit a broad spectrum of harmonics. HHG is commonly described in terms of a `quantum orbits' model based on several interfering electron trajectories, thereby incorporating both quantum-mechanical effects and an intuitive picture of classical dynamics. By tuning the parameters of the driving laser field, the interplay between these trajectories can be controlled, shaping the emitted light. Mathematically, this model expresses the harmonic response as a highly oscillatory integral. Applying saddle-point methods to this integral allows it to be decomposed into contributions from distinct saddle points of the semi-classical action, thereby linking quantum dynamics to classical trajectories. However, a general framework for applying these methods across arbitrary parameters and laser configurations has been missing.
In this thesis, we introduce Picard--Lefschetz theory and develop practical numerical methods for its application. These enable the evaluation of oscillatory integrals and identification of contributions from individual critical points.
We apply these techniques to strong-field ionisation and HHG, focusing on caustics -- enhancement features where trajectories coalesce and standard approximations fail. Our methods remain valid in these regions, allowing systematic analysis of parameter regimes and revealing previously inaccessible features. This work improves the understanding and control of ultrafast light--matter interactions.
When matter is exposed to strong fields, electrons are ionized and a number of processes such as high harmonic generation (HHG), above threshold ionization (ATI), non-sequential ionization happen. These processes are usually investigated semiclassically, i.e. matter is treated by quantum mechanics, and radiation is treated classically. In particular HHG which is the generation of ultrashort coherent extreme ultraviolet (XUV) pulses, has enabled a wealth of novel ultrashort spectroscopic methods in atomic, molecular, and condensed matter physics. Further, ionization has applications in micro-machining. Strong field physics continues to have numerous open questions. As an illustration, the mechanisms underlying dominant HHG pathways in solids remain insufficiently understood. It is not clear which mechanism dominates HHG in various materials. Possible pathways are the interband and intraband contributions. Furthermore, many of these processes in strong field physics have been described by single active electron (SAE) approximation. This approach relies on the assumption that the interaction between a single electron and the intense laser field dominates all other interactions. However, this is clearly not the case especially in solids where interactions with the lattice and other electrons can be important. In addition, quantum optical aspects of strong field physics have generated increasing interest, yet they remain partially explored. This chapter involves treating electromagnetic fields using the formalism of second quantization. Quantum optical properties, such as squeezing, entanglement and the negativity of the Wigner function, are of fundamental importance for the field of quantum information and quantum computation. The central theme of my thesis is the development of quantum optical and statistical theoretical frameworks for describing intense laser field processes such as ionization and HHG in atomic and molecular gases and in solids. The second chapter focuses on gaining a detailed understanding of the mathematical steps involved in the Keldysh theory of ionization to establish a solid theoretical foundation. Through this analysis we identified an additional factor of two compared to Keldysh’s original derivation of atomic ionization rates. The third chapter addresses the inadequately understood mechanisms that govern dominant high harmonic generation (HHG) pathways in solids. One approach to clarify these mechanisms involves introducing real (resonant) and virtual processes. We developed the strong field adiabatic following (SFAF) formalism which is based on Dyson expansion using the von-Neumann equation of density matrix. Using SFAF a diagnostic method is obtained to separate virtual and resonant channels. Through this separation, and by comparing with experimental results, we identified the need to incorporate many-body effects. The fourth chapter explores the fact that a solid is a complex many body system in which the SAE approximation is very crude. The electron interacts with holes, other electrons, collective excitations such as plasmons, and phonons. Our work is based on the idea that all theses effects can be treated as a quantum statistical heat bath of bosonic harmonic oscillators. We apply this to our SFAF formalism. In the ionization case, this work has settled a long standing issue which is the fact that simple phenomenological approaches such as the relaxation time approximation 𝑇2 result in nonphysical enhancement of ionization. In chapter five, we have investigated ways to transfer quantum optical properties on the otherwise classical high harmonic radiation. The goal is to use HHG to scale quantum sources to smaller wavelengths. Quantum properties can be imprinted by perturbing HHG with a quantum field, such as bright squeezed vacuum (BSV). In this chapter, the theory of quantum sideband high harmonic generation (QSHHG) in atoms and solids is derived to find ways by which to transfer quantum properties from the perturbative BSV to the harmonic sideband. The theoretical framework is a quantum generalization of the semi-classical Lewenstein model of HHG. It gives closed-form solutions for the HHG and QSHHG wavefunctions. Knowing the wavefunction, we can identify the quantum properties of QSHHG. The additional photons absorbed and emitted from the quantum perturbation, here BSV, create entanglement between individual harmonic sidebands and between the harmonic sidebands and the BSV. We show how this entanglement can be used to create a variety of non-classical states commonly used in quantum information science, such as high-purity single-photon states, Schrödinger cat states, and photon-added squeezed vacuum states. Some of these non-Gaussian states with negative Wigner functions have been shown to provide a quantum computational advantage over their Gaussian counterparts. Additionally, they play a significant role in quantum metrology, enhancing precision measurements beyond classical limits. This chapter has opened a path to quantum engineering of HHG.
The generation and control of extreme ultraviolet (XUV) radiation by high harmonic generation (HHG) have advanced ultrafast science, providing direct insights into electron dynamics on their natural time scale. Attosecond science has established the capability to resolve ultrafast quantum phenomena in matter by characterizing and controlling the classical properties of the high harmonics. Recent theoretical proposals have introduced novel schemes for generating and manipulating XUV HHG with distinct quantum features, paving the way to attosecond quantum optics. In this work, we transfer fundamental concepts in quantum optics into attosecond science. By driving the HHG process with a combination of an infrared bright squeezed vacuum (BSV, a non-classical state of light), and a strong coherent field, we imprint the quantum correlations of the input BSV onto both the ultrafast electron wavefunction and the harmonics' field. Performing in-situ HHG interferometry provides an insight into the underlying sub-cycle dynamics, revealing squeezing in the statistical properties of one of the most fundamental strong-field phenomena -- field induced tunneling. Our measurement allows the reconstruction of the quantum state of the harmonics through homodyne-like tomography, resolving correlated fluctuations in the harmonic field that mirror those of the input BSV. By controlling the delay between the two driving fields, we manipulate the photon statistics of the emitted attosecond pulses with sub-cycle accuracy. The ability to measure and control quantum correlations in both electrons and XUV attosecond pulses establishes a foundation for attosecond electrodynamics, manipulating the quantum state of electrons and photons with sub-cycle precision.
High harmonic generation (HHG) is an extremely nonlinear effect generating coherent broadband radiation and pulse durations reaching attosecond timescales. Conventional models of HHG that treat the driving and emitted fields classically are usually very successful but inherently cannot capture the quantum-optical nature of the process. Although prior work considered quantum HHG, it remains unknown in what conditions the spectral and statistical properties of the radiation depart considerably from the known phenomenology of HHG. The discovery of such conditions could lead to novel sources of attosecond light having squeezing and entanglement. Here, we present a fully-quantum theory of extreme nonlinear optics, predicting quantum effects that alter both the spectrum and photon statistics of HHG, thus departing from all previous approaches. We predict the emission of shifted frequency combs and identify spectral features arising from the breakdown of the dipole approximation for the emission. Our results show that each frequency component of HHG can be bunched and squeezed and that each emitted photon is a superposition of all frequencies in the spectrum, i.e., each photon is a comb. Our general approach is applicable to a wide range of nonlinear optical processes, paving the way towards novel quantum phenomena in extreme nonlinear optics.
The generation of higher-order harmonic radiation originating from the interaction of intense laser pulses with matter is typically described semiclassically: While the electronic structure and dynamics of matter is described quantum mechanically, the intense light field is described classically---and accordingly the generated harmonic radiation. However, pioneering experiments on quantum optical properties of high-order harmonic generation (HHG) in atomic gases from the group of P. Tzallas [I. Gonoskov et al., Sci. Rep. 6, 32821 (2016) and N. Tsatrafyllis et al., Nat. Commun. 8, 15170 (2017)], and theoretical investigations from the group of M. Lewenstein [M. Lewenstein et al., Nat. Phys. 17, 1104 (2021) and P. Stammer et al., Phys. Rev. Lett. 128, 123603 (2022)] have impressively demonstrated that light's quantum properties are observable in the strong-field realm. In this paper, we develop a quantum optical description of HHG in a bulk semiconductor originating from the nonlinear Bloch current. This mechanism of HHG, known as the intraband current, constitutes the major contribution to the emission spectrum for harmonic with energies of quanta roughly below the bandgap. Under certain approximations, employing a quantum description of both light and matter, we obtain analytical solutions, which allow us to analyze the classical and quantum optical characteristics of the fundamental mode of light and the harmonic modes. We find intricate but sufficiently mild modifications of the photon statistics of the fundamental mode and coherent displacements depending on the parameters of the driving laser field. Similar to high-harmonic generation in atoms, the fundamental and emitted harmonic field modes are entangled. Moreover, our analytical model predicts parameter ranges where these quantum optical properties will be most pronounced, allowing protocols for quantum information processing with high photon numbers over a large range of frequencies.
High-order harmonic generation (HHG) provides a broad spectral bandwidth for synthesizing attosecond pulses. However, in the current HHG schemes, only part of the harmonics can be phase-locked, which limits the ability to achieve shorter attosecond pulses. Here, we study attosecond pulse synthesis from HHG of an atom driven by an intense quantum light, i.e., squeezed light. It is interestingly found that the harmonics in the whole spectrum can be phase-locked and, by using these harmonics, the width of the synthesized attosecond pulse is greatly reduced. By developing strong-field approximation theory in squeezed light, the physics of the phase-locked harmonic generation throughout the HHG spectrum is revealed and is found to be independent of the target system. Furthermore, we uncover the dependence of the synthesized attosecond pulse width on the squeezing parameter of the squeezed light. Our findings provide a robust tool for obtaining phase-locked harmonics throughout the HHG spectrum for synthesizing short attosecond pulses.
Intense light-matter interaction largely relies on the use of high-power light sources, creating fields comparable to, or even stronger than, the field keeping the electrons bound in atoms. Under such conditions, the interaction induces highly nonlinear processes such as high harmonic generation, in which the low-frequency photons of a driving laser field are upconverted into higher-frequency photons. These processes have enabled numerous groundbreaking advances in atomic, molecular, and optical physics, and they form the foundation of attosecond science. Until recently, however, such processes were typically described using semi-classical approximations, since the quantum properties of the light field were not required to explain the observables. This has changed in the recent past. Ongoing theoretical and experimental advances show that fully quantized descriptions of intense light-matter interactions, which explicitly incorporate the quantum nature of the light field, open new avenues for both fundamental research and technological applications at the fully quantized level. These advances emerge from the convergence of quantum optics with strong-field physics and ultrafast science. Together, they have given rise to the field of quantum optics and quantum electrodynamics of strong-field processes.