In recent years, strong-field processes such as high-order harmonic generation (HHG) and above-threshold ionization driven by nonclassical states of light have become an increasingly popular field of study. The theoretical modeling of these processes often applies an approximate phase-space expansion of the nonclassical driving field in terms of coherent states, which has been shown to accurately predict the harmonic spectrum. However, its accuracy for the computation of quantum optical observables like the degree of squeezing and photon statistics has not been thoroughly considered. In this work, we introduce this approximative phase-space description and discuss its accuracy, and we find that it mischaracterizes the quantum optical properties of the driving laser by making it an incoherent mixture of classical states. We further show that this error in the driving field description maps onto the light emitted from HHG, as neither sub-Poissonian photon statistics nor quadrature squeezing below vacuum fluctuations can be captured by the approximative phase-space description. Lastly, to benchmark the approximative phase-space description, we consider the quantum HHG from a one-band model, which yields an exact analytical solution. Using the approximative phase-space representation with this specific model, we find a small quantitative error in the quadrature variance of the emitted field that scales with pulse duration and emitter density. Our results show that using this approximative phase-space description can mischaracterize quantum optical observables. Attributing physical meaning to such results should therefore be accompanied by a quantitative analysis of the error.
In recent years, strong-field processes such as high-order harmonic generation (HHG) and above-threshold ionization driven by nonclassical states of light have become an increasingly popular field of study. The theoretical modeling of these processes often applies an approximate phase-space expansion of the nonclassical driving field in terms of coherent states, which has been shown to accurately predict the harmonic spectrum. However, its accuracy for the computation of quantum optical observables like the degree of squeezing and photon statistics has not been thoroughly considered. In this work, we introduce this approximative phase-space description and discuss its accuracy, and we find that it mischaracterizes the quantum optical properties of the driving laser by making it an incoherent mixture of classical states. We further show that this error in the driving field description maps onto the light emitted from HHG, as neither sub-Poissonian photon statistics nor quadrature squeezing below vacuum fluctuations can be captured by the approximative phase-space description. Lastly, to benchmark the approximative phase-space description, we consider the quantum HHG from a one-band model, which yields an exact analytical solution. Using the approximative phase-space representation with this specific model, we find a small quantitative error in the quadrature variance of the emitted field that scales with pulse duration and emitter density. Our results show that using this approximative phase-space description can mischaracterize quantum optical observables. Attributing physical meaning to such results should therefore be accompanied by a quantitative analysis of the error.
Quantum-optical descriptions of strong-field processes have attracted significant attention in recent years. Typically, the theoretical modeling has been conducted in the Schrödinger picture, where results are only obtainable under certain approximations, while, in contrast, the Heisenberg picture has remained relatively unexplored. In this work, we develop an accurately controlled perturbative expansion of the time-evolution operator in the Heisenberg picture and derive beyond-semiclassical corrections to the emitter dynamics due to the coupling to the quantized electromagnetic field, capturing effects of the quantum fluctuations present in the latter. We focus on high-order harmonic generation (HHG), where the approach is accurate in parameter regimes of current interest and it gives closed-form expressions for key observables. This formulation not only simplifies numerical calculations compared to the Schrödinger-picture approach but also provides a clear correspondence between nonclassical features of the emitted light and the underlying induced dynamics of the generating medium including quantum fluctuations. Moreover, the Heisenberg framework naturally yields scaling relations with the number of independent emitters, enabling us to assess whether nonclassical behavior should persist under typical experimental conditions involving large emitter ensembles. Interestingly, we find that the degree of squeezing increases with the number of emitters, whereas the photon statistics approaches a classical Poissonian distribution in the many-emitter limit. We also find that the beyond-semiclassical emitter dynamics significantly enhances the degree of squeezing of the emitted light. Our work advances the theoretical understanding of quantum-optical HHG and introduces an accessible and well-controlled framework to describe realistic experiments.
Strong-field high-harmonic generation (HHG) is conventionally described semiclassically, with the driving laser treated as a classical field. This approximation becomes insufficient in nanoscale interaction geometries, where extreme spatial confinement raises the vacuum-field amplitude to the ~10^-2 level relative to the driving-field amplitude. When the quantum fluctuations of the driving field are redistributed between conjugate quadratures by squeezing, they can be directly imprinted onto macroscopic HHG observables. To model this interaction, we employ a Wigner phase-space approach that maps the quantum-optical driver onto a stochastic ensemble of time-dependent Schrodinger equation realizations. Although each realization remains classically simulable, the sub-vacuum quadrature covariance structure of squeezed states cannot be reproduced by any field admitting a non-negative Glauber-Sudarshan P-representation. Within this single-mode Gaussian framework, we show that amplitude squeezing suppresses the shot-to-shot variance of the HHG cutoff below the standard quantum limit (SQL). To leading order in the vacuum-to-driver ratio, the normalized cutoff variance decays exponentially with the squeezing parameter, independent of the absolute vacuum-field amplitude and therefore robust against uncertainties in the effective nanoscale mode volume. A subleading phase-noise contribution from the anti-squeezed quadrature produces a variance minimum near r_opt ~ 1.6, providing a concrete experimental target. These results establish the HHG cutoff variance ratio as a nonlinear, self-calibrating Heisenberg witness in which sub-SQL scaling directly reflects the redistribution of quantum uncertainty in the driving field.
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 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.
High-harmonic generation (HHG) driven by intense quantum light has recently emerged as a new frontier in strong-field physics. Yet, the mechanism underlying HHG with nonclassical fields such as bright squeezed vacuum (BSV) light—characterized by a vanishing mean electric field—remains elusive. In this work, we develop a fully quantum theory for HHG under arbitrary quantum driving fields. We show that the fluctuation component of the dipole correlation operator plays a central role in HHG with quantum light and is essential for BSV-driven HHG. Furthermore, our framework reveals the photon statistics distribution within individual harmonics and across different harmonics, demonstrating that the statistics inherit the statistical properties of the driving field. This theory overcomes key limitations of existing approaches and uncovers new quantum features of HHG in the presence of nonclassical light.
High-harmonic generation (HHG) is a prototypical strong-field process in which intense light drives matter to emit radiation at integer multiples of the driving frequency. Extending HHG into the quantum-optical regime offers new opportunities to probe and control strongly nonlinear light-matter interactions using nonclassical states of light. Yet describing this regime requires a fully quantum treatment of the correlated electron-photon dynamics, which becomes computationally challenging for broadband, strongly squeezed fields. Here we solve the quantum-electrodynamical dynamics of a two-level system driven by bright squeezed vacuum in a converged multimode Hilbert space. Both the driving field and emitted harmonics are fully quantized, with the light-matter interaction treated nonperturbatively. This enables direct access to the multimode quantum state and its higher-order correlations beyond semiclassical sampling or perturbative descriptions. We show that squeezed-vacuum driving produces harmonic emission with qualitatively distinct second- and third-order photon correlations compared with coherent excitation. Moreover, back-action from the driven emitter strongly reshapes the incident squeezed field, generating pronounced Wigner-function negativities that evolve on attosecond timescales. Our results establish a fully quantum framework for broadband strong-field dynamics with squeezed light and provide a route to predicting and interpreting quantum-HHG experiments and their extension to more complex emitters.