Nonperturbative regime of low-order harmonic generation in intense low-frequency laser field

This paper demonstrates that while standard perturbation theory fails for atomic responses in intense low-frequency laser fields above $0.6 \times 10^{14}$ W/cm$^2$, a Padé expansion fitted to numerical TDSE solutions effectively describes the nonperturbative regime up to $1.4 \times 10^{14}$ W/cm$^2$ and successfully models the intensity-dependent growth of third and fifth harmonic generation and optical rectification, despite limitations in predicting the nonlinear refractive index.

The Gist

Imagine an atom as a tiny, spring-loaded trampoline. When you shine a light (a laser) on it, the light's electric field pushes and pulls on the trampoline, making it bounce. This bouncing creates a new kind of light called a "harmonic," which is like a higher-pitched echo of the original laser.

For a long time, scientists thought they could predict exactly how this trampoline would bounce using a simple rule: the harder you push, the more it bounces, in a perfectly straight line. This is called "perturbation theory." It works great for gentle pushes (weak lasers).

However, this paper investigates what happens when you push that trampoline really hard with an intense laser. The authors, S. A. Bondarenko and V. V. Strelkov, found that the simple straight-line rule breaks down completely.

Here is a breakdown of their findings using everyday analogies:

1. The "Straight Line" Breaks (The Problem)

When the laser gets too strong (specifically, above a certain intensity), the trampoline stops behaving like a simple spring.

• The Old Way: Scientists tried to fix the broken rule by just adding more terms to their math, like saying, "Okay, maybe the bounce isn't just straight; maybe it curves a little bit, then a lot, then a tiny bit more." They kept adding these "curves" (higher-order non-linearities) to their equations.
• The Reality: No matter how many extra curves they added, the math still didn't match what was actually happening in the computer simulation. The trampoline was doing something the straight-line logic simply couldn't predict. It was entering a "non-perturbative" regime—a fancy way of saying the rules of the game had changed, and the old playbook was useless.

2. The New Map (The Padé Solution)

Instead of trying to force the trampoline into a straight line or a series of curves, the authors tried a different approach. They looked at the actual data from their super-computer simulations (solving the Schrödinger equation, which is the master rulebook for how quantum particles move).

They found that the trampoline's behavior looked like it was heading toward a "cliff" or a singularity at a specific strength of push. To describe this, they used a Padé approximation.

• The Analogy: Imagine trying to draw a map of a winding mountain road. A polynomial series (the old way) tries to draw it using only straight lines and gentle curves, which eventually fails to capture the sharp turns. The Padé approximation is like using a flexible, stretchy rubber band that can snap into the exact shape of the road, even if it has a sharp cliff or a loop.
• The Result: This new "rubber band" map fit the computer data perfectly, even when the laser was very strong (up to about $1.4 \times 10^{14}$ W/cm²). It worked for both weak pushes and strong ones.

3. The "Nonlinear Oscillator" Model

Once they had this perfect map of how the trampoline behaves in a static (non-moving) field, they wanted to see if they could use it to predict what happens when the laser is actually oscillating (wiggling back and forth).

They built a Nonlinear Oscillator Model.

• The Analogy: Think of a child on a swing. If you push the swing gently, it moves back and forth predictably. If you push it hard, the swing's chains might stretch, or the seat might tilt, changing how it moves. The authors created a mathematical "swing" where the restoring force (the pull back to the center) was defined by their new "rubber band" map.
• What it Got Right: This model successfully predicted how the efficiency of creating new light (harmonics) grows as the laser gets stronger. It worked well for:
  • Creating the 3rd and 5th "echoes" (harmonics) of the light in infrared fields.
  • Creating a steady "rectified" field (like turning AC power into DC) using two different colored lasers.
• What it Got Wrong: The model failed to predict the behavior of the refractive index (how much the light bends or slows down) in the non-perturbative zone.
  • Why? The model treats the atom as a perfect, closed system. In reality, when the laser is that strong, it starts ripping electrons off the atom (photoionization). These free electrons act like a crowd of people running around the trampoline, messing up the bounce. The model didn't account for these "runaway" electrons, nor did it account for specific resonances (when the laser frequency accidentally matches the atom's natural vibration).

Summary

The paper is essentially a story about knowing when to stop using old maps.

1. Old Map (Perturbation Theory): Works for weak lasers, fails for strong ones. Adding more details to the map didn't help.
2. New Map (Padé Approximation): A flexible mathematical tool that fits the actual data perfectly for strong lasers, up to the point where the atom starts breaking apart (ionizing).
3. The Simulation (The Oscillator Model): Using this new map, they built a model that correctly predicts how efficiently new light is created in strong fields. However, it cannot predict how the light bends (refractive index) because it ignores the messy reality of electrons being ripped out of the atom.

In short: They found a better way to describe how atoms react to intense light, but only up to the point where the light becomes so strong it starts destroying the atom.

Technical Summary

Technical Summary: Nonperturbative Regime of Low-Order Harmonic Generation in Intense Low-Frequency Laser Fields

Problem Statement
The paper addresses the limitations of perturbation theory in describing the atomic response to intense, low-frequency laser fields. While perturbation theory successfully models harmonic generation at moderate intensities, it fails to accurately describe the system when laser intensities exceed approximately $0.6 \times 10^{14}$ W/cm$^2$. In this regime, the response deviates significantly from the linear dependence on intensity predicted by low-order nonlinearities, and standard polynomial expansions (Taylor series) fail to converge or provide satisfactory accuracy, even when higher-order terms are included. The authors aim to characterize this nonperturbative behavior in the quasi-static limit (before photoionization becomes the dominant mechanism) and develop a model capable of describing the atomic response beyond the perturbative domain.

Methodology
The study employs a multi-faceted approach combining numerical simulation, analytical approximation, and physical modeling:

1. Numerical Solution of TDSE: The authors numerically solve the three-dimensional time-dependent Schrödinger equation (TDSE) for a model argon atom in a linearly polarized external field. The simulation uses a single-active electron approximation with an effective potential. To isolate the bound-state response, the numerical box is sized to absorb detached wave packets, and the dipole moment is corrected for wave-function depletion due to photoionization.
2. Quasi-Static Analysis: By analyzing the TDSE results for very low-frequency pulses, the authors reconstruct the atomic dipole moment $d_{qs}(E)$ as a function of instantaneous field strength. They compare this numerical data against:
   • Polynomial series expansions (using least squares and Chebyshev approximations).
   • A Padé approximation, which fits the numerical data with a rational function designed to capture the singularity behavior near the barrier-suppression ionization limit.
3. Nonlinear Oscillator Model: To extend the quasi-static findings to oscillating fields, the authors propose a nonlinear oscillator model. The restoring force $F(d)$ in the oscillator equation is defined implicitly by the Padé-fitted quasi-static response. This model is solved numerically to predict the dispersion of nonlinear polarizabilities for finite-frequency fields.

Key Contributions and Results

• Breakdown of Perturbation Theory: The study confirms that for intensities above $\sim 0.6 \times 10^{14}$ W/cm$^2$, polynomial series approximations of the dipole moment lose precision. Adding higher-order terms does not improve the fit, indicating a fundamental breakdown of the perturbative approach in this regime.
• Padé Approximation: The authors successfully fit the numerical TDSE results for the quasi-static dipole moment using a Padé expansion. This approximation accurately describes the response from weak fields up to intensities of approximately $1.4 \times 10^{14}$ W/cm$^2$, effectively bridging the gap between the perturbative and nonperturbative regimes before photoionization dominates.
• Nonlinear Oscillator Model Performance:
  • Harmonic Generation: The model successfully predicts the nonperturbative growth of efficiency for third- and fifth-harmonic generation in infrared (IR) fields. The results align well with TDSE calculations for IR frequencies.
  • Optical Rectification: The model accurately describes the generation of quasi-static fields (optical rectification) in two-color fields, correctly reproducing the rapid growth of polarizability with intensity and its dependence on the phase difference between the two colors.
  • Limitations (Refractive Index): The model fails to predict the behavior of the nonlinear refractive index (first-order polarizability) in the nonperturbative domain. Specifically, it predicts a nonperturbative increase in the refractive index, whereas TDSE results show a decrease or saturation. The authors attribute this discrepancy to the model's exclusion of photoelectron contributions (dominant in IR) and multiphoton resonances (dominant in UV), which are critical for the refractive index but not for the higher-order harmonic generation efficiencies in the studied regimes.
  • Phase Behavior: The model predicts zero harmonic phases under steady-state conditions, whereas TDSE results show a non-zero phase delay at high intensities. The authors note that the oscillator model only captures non-zero phases in a narrow range near the ionization threshold.

Significance and Claims
The paper claims to provide a robust alternative to perturbation theory for describing low-order harmonic generation in intense laser fields. By utilizing a Padé expansion derived from numerical TDSE solutions, the authors establish a method to describe the atomic response in the nonperturbative regime up to the onset of significant photoionization.

The primary significance lies in the demonstration that while the nonlinear oscillator model fails to capture the complex physics of the nonlinear refractive index (due to missing photoelectron and resonance effects), it remains a valid and effective tool for describing the efficiency growth of third- and fifth-harmonic generation and optical rectification in IR and two-color fields. The work delineates the specific boundaries where perturbative descriptions fail and where simplified physical models can still offer accurate predictions for specific nonlinear optical processes, provided the dominant mechanisms (photoionization and resonances) are accounted for or are negligible.