Quantum-Matter Lewenstein Model

• Quantum-Matter Lewenstein Model is a quantum-optical generalization of the classical Lewenstein HHG theory that quantizes the emitted harmonic field from an atom-field Hamiltonian.
• It decomposes harmonic emission into two channels—interband-like recombination and intraband-like ionization nonlinearity—linking gas-phase and solid-state phenomena.
• The model introduces a vacuum field scaling factor that alters the spectral roll-off, enabling a unified description of high-harmonic generation across different media.

Searching arXiv for related HHG Lewenstein-model papers to support disambiguation and context. The quantum-matter Lewenstein model is a quantum-optical generalization of the standard Lewenstein strong-field description of high-harmonic generation (HHG) in gases. It starts from the full atom plus quantized-field Hamiltonian and treats the emitted harmonics as quantized-field excitations rather than only as a classical radiated current. In this formulation, the harmonic field amplitude decomposes into two channels, one interband-like and one intraband-like, and both are multiplied by the vacuum electric field amplitude. The resulting framework recovers the familiar gas-phase ionization–continuum propagation–recombination picture while simultaneously establishing a direct correspondence with interband and intraband HHG in solids (Thorpe et al., 19 Aug 2025).

1. Definition and formal scope

In the quantum-matter formulation, the starting Hamiltonian is

$\hat{H}=\hat{H}_m+\hat{H}_f+\hat{H}_i,$

with

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$

The electric field is quantized mode by mode as

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$

with vacuum field amplitude

$$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$

The proportionality $\eta_\kappa\propto \omega_\kappa^{-1/2}$ is identified as the origin of the modified frequency scaling and of the faster roll-off of harmonic power with order (Thorpe et al., 19 Aug 2025).

This construction differs from semiclassical HHG theories in a specific way: the emitted field is quantized from the outset, while the strong driving field is separated out as a classical coherent component. A plausible implication is that the model is not merely a reformulation of the standard gas-phase theory, but a change in what is taken as the primary observable. Instead of inferring radiation solely from a classical current, the theory computes photon-number observables in a mixed matter–field Hilbert space.

2. Matter–field reduction and Lewenstein-type ansatz

The model uses the standard strong-field, dipole-level Lewenstein reduction for a hydrogen-like atom with a ground state $|0\rangle$ and continuum states $|\mathbf p\rangle$. The wavefunction ansatz is

$$|\Psi(t)\rangle = |0\rangle\otimes |\phi_0(t)\rangle + \int d^3p\,|\mathbf p\rangle\otimes |\phi_{\mathbf p}(t)\rangle.$$

This decomposes the state into a bound component and a continuum component labeled by momentum $\mathbf p$. The HHG signal is extracted from the emitted photon number expectation value,

$$\langle \hat n\rangle(t)\approx \langle \phi_0(t)|\hat n|\phi_0(t)\rangle.$$

To separate the strong classical pump from the quantized emitted harmonics, the derivation uses the interaction picture, a velocity-gauge transformation, and a displacement operator for the strong coherent driver. The field thereby splits as

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$0

where $$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$1 is classical and $$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$2 is the quantized harmonic field. The resulting dynamics reduce to a Volterra integro-differential equation for the ground-state photon wavefunction,

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$3

with coupling operator

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$4

and classical action

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$5

Under the same approximation spirit as Lewenstein—slowly varying photon wavefunction, first-order operator expansion, and neglect of strong depletion—the solution becomes

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$6

The emitted harmonic field is therefore a coherent state. This preserves the Lewenstein computational architecture while shifting the interpretation of the emitted radiation into an explicitly quantum-optical form (Thorpe et al., 19 Aug 2025).

3. Two-channel structure of harmonic emission

The central structural result is the decomposition

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$7

The first contribution,

$$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$8

is the interband-like channel. Here $$\hat{H}_m=\frac{\hat{\mathbf p}^2}{2m}+V(\mathbf x), \qquad \hat{H}_f=\sum_\kappa \hbar\omega_\kappa \hat a_\kappa^\dagger \hat a_\kappa, \qquad \hat{H}_i = |e|\,\mathbf r\cdot \hat{\mathbf F}(\mathbf x).$$9 is the Lewenstein time-dependent dipole,

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$0

This is the familiar gas-phase strong-field picture of ionization, propagation in the continuum, and recombination to the ground state.

The second contribution,

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$1

with

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$2

is the intraband-like channel. In the formulation of (Thorpe et al., 19 Aug 2025), this second channel was not emphasized in the original semiclassical Lewenstein discussion. It arises from the ionization nonlinearity and from the dressing of the continuum electron by coherent radiation modes.

A common simplification is to identify the gas-phase Lewenstein model solely with recombination radiation. The quantum-matter construction rejects that reduction. It preserves the recombination channel, but also exposes a second source term whose structure is analogous to intraband HHG in solids. This suggests that the standard gas/solid distinction can be reorganized around a common two-channel architecture rather than around separate phenomenological mechanisms.

4. Current representation and gas–solid correspondence

Both channels can be rewritten in a semiclassical current form: $$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$3

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$4

Here

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$5

is the interband current in the gas-phase mapping, while

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$6

is the intraband current, with continuum velocity

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$7

and continuum amplitude

$$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$8

In both cases the current is multiplied by the vacuum field factor $$\hat{\mathbf F}(\mathbf x) = i\sum_\kappa \eta_\kappa \omega_\kappa \mathbf e_\kappa \left( \hat a_\kappa e^{i\mathbf k\cdot \mathbf x} - \hat a_\kappa^\dagger e^{-i\mathbf k\cdot \mathbf x} \right),$$9, a feature absent from the corresponding semiclassical treatments (Thorpe et al., 19 Aug 2025).

The gas/solid mapping is explicit.

Sector	Gas-phase identification	Solid-state identification
Momentum variable	$$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$0 is the free-electron momentum	$$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$1 is replaced by crystal momentum $$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$2
Interband channel	Lewenstein recombination dipole	Valence–conduction coherence
Intraband channel	Ionization nonlinearity and continuum dressing; related to Brunel HHG	Carrier motion within bands with $$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$3

In gases, the continuum is described as basically linear, so intraband-like HHG is not due to band curvature. In solids, the intraband contribution is associated with carrier acceleration within bands. The unified perspective is therefore not based on identical microscopic carriers, but on a shared decomposition into coherence-generated emission and population-transport emission. A plausible implication is that the distinction between gases and solids shifts from channel taxonomy to the detailed origin of the corresponding currents.

5. Spectral observables, scaling, and crossover behavior

The total photon number is written as

$$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$4

After the continuum limit and angular simplifications, the spectral density becomes

$$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$5

and the yield around harmonic $$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$6 is

$$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$7

The $$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$8 dependence cancels because $$\eta_\kappa=\sqrt{\frac{\hbar}{2\omega_\kappa \varepsilon_0 V}}.$$9, so the yield is volume-independent (Thorpe et al., 19 Aug 2025).

The principal quantitative correction follows from $\eta_\kappa\propto \omega_\kappa^{-1/2}$0. As a result, harmonic amplitudes are smaller than in previous semiclassical quantum-matter theories by a frequency-dependent factor, and harmonic power falls off faster with harmonic order. The model is therefore presented as a quantitative correction, not just a conceptual reinterpretation. In the numerical results for argon, the interband-like channel dominates most of the HHG plateau and determines the cutoff behavior, whereas the intraband-like channel dominates at low harmonic orders.

The crossover or switchover harmonic is defined as the highest odd harmonic where the intraband spectrum is comparable to or exceeds the interband spectrum. When plotted versus driving wavelength, the crossover follows a linear trend in harmonic order corresponding to a constant photon energy,

$\eta_\kappa\propto \omega_\kappa^{-1/2}$1

Since $\eta_\kappa\propto \omega_\kappa^{-1/2}$2, a fixed crossover energy implies

$\eta_\kappa\propto \omega_\kappa^{-1/2}$3

This provides a wavelength-independent criterion in photon-energy space for the transition between the two mechanisms. The same framework is also checked against the standard gas HHG cutoff law

$\eta_\kappa\propto \omega_\kappa^{-1/2}$4

showing that the interband-like channel reproduces the usual plateau and cutoff structure (Thorpe et al., 19 Aug 2025).

6. Conceptual significance and terminological disambiguation

The significance of the quantum-matter Lewenstein model lies in the combination of three claims: it preserves the intuitive Lewenstein three-step picture, it quantizes the emitted field so that vacuum fluctuations enter explicitly, and it reveals a second HHG channel that maps naturally to intraband emission. Because the model is still computationally efficient, especially after saddle-point treatment, it can be incorporated into macroscopic propagation and phase-matching simulations for realistic HHG experiments (Thorpe et al., 19 Aug 2025).

One conceptual misconception addressed by this framework is that quantizing the emitted field would merely reproduce the classical emission formulas with operator notation. The explicit $\eta_\kappa\propto \omega_\kappa^{-1/2}$5 factor prevents such a reduction, because it changes the scaling of harmonic amplitudes and hence the spectral roll-off. A second misconception is that the gas-phase Lewenstein model contains only the recombination channel; the two-channel decomposition shows that, in this quantum-optical treatment, the gas phase also contains an intraband-like source term related to ionization nonlinearity and continuum dressing.

The term “Lewenstein model” also has an unrelated usage in quantum game theory. In that context, it denotes the Eisert–Wilkens–Lewenstein scheme, treated as a quantum extension of a standard $\eta_\kappa\propto \omega_\kappa^{-1/2}$6 strategic-form game and, in one interpretation, as a nontrivial extension of the mixed extension of a classical game (Frąckiewicz, 2020). The quantum-matter Lewenstein model is distinct from that usage: it concerns HHG, atom–field dynamics, quantized emitted radiation, and the unification of gas and solid HHG channels.