This paper has been prepared by the Symphony collaboration (University of Warsaw, Uniwersytet Jagielloński, DESY/CNR and ICFO) on the occasion of the 25th anniversary of the ‘simple man’s models’ which underlie most of the phenomena that occur when intense ultrashort laser pulses interact with matter. The phenomena in question include high-harmonic generation (HHG), above-threshold ionization (ATI), and non-sequential multielectron ionization (NSMI). ‘Simple man’s models’ provide both an intuitive basis for understanding the numerical solutions of the time-dependent Schrödinger equation and the motivation for the powerful analytic approximations generally known as the strong field approximation (SFA). In this paper we first review the SFA in the form developed by us in the last 25 years. In this approach the SFA is a method to solve the TDSE, in which the non-perturbative interactions are described by including continuum–continuum interactions in a systematic perturbation-like theory. In this review we focus on recent applications of the SFA to HHG, ATI and NSMI from multi-electron atoms and from multi-atom molecules. The main novel part of the presented theory concerns generalizations of the SFA to: (i) time-dependent treatment of two-electron atoms, allowing for studies of an interplay between electron impact ionization and resonant excitation with subsequent ionization; (ii) time-dependent treatment in the single active electron approximation of ‘large’ molecules and targets which are themselves undergoing dynamics during the HHG or ATI processes. In particular, we formulate the general expressions for the case of arbitrary molecules, combining input from quantum chemistry and quantum dynamics. We formulate also theory of time-dependent separable molecular potentials to model analytically the dynamics of realistic electronic wave packets for molecules in strong laser fields. We dedicate this work to the memory of Bertrand Carré, who passed away in March 2018 at the age of 60.

Topological materials are of interest to both fundamental science and advanced technologies because topological states are robust with respect to perturbations and dissipation. Experimental detection of topological invariants is thus in great demand but remains extremely challenging. Ultrafast laser-matter interactions, and in particular high-harmonic generation, meanwhile, were proposed several years ago as tools to explore the structural and dynamical properties of various matter targets. Here we show that the high-harmonic emission signal produced by a circularly polarized laser contains signatures of topological phases and transitions in the paradigmatic Haldane model. In addition to clear shifts of the overall emissivity and harmonic cutoff, the high-harmonic emission shows a unique circular dichroism, which exhibits clear changes in behavior at the topological phase boundary. Our findings pave the way to understand fundamental questions about the ultrafast electron-hole pair dynamics in topological materials via nonlinear high-harmonic generation spectroscopy.

We investigate theoretically the strong-field regime of light-matter interactions in the topological-insulator class of quantum materials. In particular, we focus on the process of nonperturbative high-order harmonic generation from the paradigmatic three-dimensional topological insulator bismuth selenide (${\mathrm{Bi}}_2{\mathrm{Se}}_3$) subjected to intense midinfrared laser fields. We analyze the contributions from the spin-orbit-coupled bulk states and the topological surface bands separately and reveal a major difference in how their harmonic yields depend on the ellipticity of the laser field. Bulk harmonics show a monotonic decrease in their yield as the ellipticity increases, in a manner reminiscent of high harmonic generation in gaseous media. However, the surface contribution exhibits a highly nontrivial dependence, culminating with a maximum for circularly polarized fields. We attribute the observed anomalous behavior to (i) the enhanced amplitude and the circular pattern of the interband dipole and the Berry connections in the vicinity of the Dirac point and (ii) the influence of the higher-order, hexagonal warping terms in the Hamiltonian, which are responsible for the hexagonal deformation of the energy surface at higher momenta. The latter are associated directly with spin-orbit-coupling parameters. Our results thus establish the sensitivity of strong-field-driven high harmonic emission to the topology of the band structure as well as to the manifestations of spin-orbit interaction.

High-order harmonic generation (HHG) is a fundamental process which can be simplified as the production of high energetic photons from a material subjected to a strong driving laser field. This highly nonlinear optical process contains rich information concerning the electron structure and dynamics of matter, for instance, gases, solids and liquids. Moreover, the HHG from solids has recently attracted the attention of both attosecond science and condensed matter physicists, since the HHG spectra can carry information of electron-hole dynamics in bands and inter- and intra-band current dynamics. In this paper, we study the effect of interlayer coupling and symmetry in two-dimensional (2D) material by analyzing high-order harmonic generation from monolayer and two differently stacked bilayer hexagonal boron nitrides (hBNs). These simulations reveal that high-order harmonic emission patterns strongly depend on crystal inversion symmetry (IS), rotation symmetry and interlayer coupling.

Progress in attosecond science has relied on advancements in few-cycle pulse generation technology and its application to high-order harmonic generation. Traditionally, self-phase modulation in bulk solids has been used for the compression of moderate-energy pulses, additionally exhibiting favorable dispersion properties for mid-infrared (mid-IR) pulses. Here, we use the anomalous dispersion of ${\mathrm{Y}}_3{\mathrm{Al}}_5{\mathrm{O}}_{12}$ (YAG) to self-compress many-cycle pulses from a 50 kHz mid-IR OPA down to produce sub-three-cycle 10 μJ pulses and further use them to generate high-order harmonics in a ZnO crystal. In agreement with theoretical predictions, we observe a boost in the harmonic yield by a factor of two, and spectral broadening of above-gap harmonics, compared to longer driving pulses. The enhanced yield results from an increase in the intensity for the self-compressed pulses.

Two-dimensional hexagonal materials such as transition metal dichalcogenides exhibit valley degrees of freedom, offering fascinating potential for valley-based quantum computing and optoelectronics. In nonlinear optics, the K and K’ valleys provide excitation resonances that can be used for ultrafast control of excitons, Bloch oscillations, and Floquet physics. Under intense laser fields, however, the role of coherent carrier dynamics away from the K/K’ valleys is largely unexplored. In this study, we observe quantum interferences in high harmonic generation from monolayer WS₂ as laser fields drive electrons from the valleys across the full Brillouin zone. In the perturbative regime, interband resonances at the valleys enhance high harmonic generation through multi-photon excitations. In the strong-field regime, the high harmonic generation is sensitively controlled by quantum interferences of laser-field-driven electrons occupying various points in the Brillouin zone, including regions far from the K/K’ valleys. Our experimental observations are in strong agreement with quantum simulations, validating their interpretation. This work proposes new routes for harnessing laser-driven quantum interference in two-dimensional hexagonal systems and all-optical techniques to occupy and read-out electronic structures in the full Brillouin zone via strong-field nonlinear optics, advancing quantum technologies.