Electron Excitation Probability in Dielectrics under Two-color Intense Laser Fields

Kenichi L. Ishikawa; Mizuki Tani; Tomohito Otobe

REVIEW 4 major objections 5 minor 56 references

Electron Excitation Probability in Dielectrics under Two-color Intense Laser Fields

T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper derives an analytical formula for the excitation rate in dielectrics under collinear two-color intense laser fields and shows that it reproduces the intensity and relative-phase trends of TDDFT simulations for α-quartz.

desk verdict New two-color Bessel formalism with a solid φ=0 derivation, but the φ≠0 result rests on an invalid modified-Bessel identity and the phase-dependence claim is unsupported. read the letter →

arxiv 2507.01023 v1 pith:WHK54JW2 submitted 2025-06-13 physics.optics physics.atom-phphysics.plasm-ph

Mizuki Tani , Kenichi L. Ishikawa , Tomohito Otobe This is my paper · ORCID

classification physics.opticsphysics.atom-phphysics.plasm-ph

keywords two-colorlaserfieldsstrong-fieldionizationdielectricsparabolictwo-bandmodelgeneralizedBesselfunctionsrelativephasecontrolTDDFTalpha-quartz

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper derives a closed analytical formula for the rate at which an intense, collinear two-color laser field ($\omega$ plus its second harmonic $2\omega$) excites electrons across the band gap of a dielectric. The formula comes from a parabolic two-band model in the Houston basis, with the transition amplitude expressed through generalized and 'twisted' multivariable Bessel functions; the relative phase $\varphi$ enters analytically through those functions. The authors benchmark Eq. (30) against real-time time-dependent density functional theory (TDDFT) simulations for $\alpha$-quartz and report that the intensity dependence and the $\varphi$-dependence of the excitation probability reproduce the TDDFT trends, including enhanced two-color excitation at high intensity and a minimum near $\varphi \simeq 2\pi$ with a maximum near $\varphi \simeq 1.4\pi$. The purpose is to give a fast, physically transparent substitute for expensive first-principles calculations in strong-field dielectrics.

What carries the argument

The workhorse is Eq. (30), the total interband transition probability per unit time and volume. It is built by writing the time-dependent orbital in the Houston basis, assuming a parabolic two-band gap $\Delta^{\mathrm{vc}}_{\mathbf{k}} = E_g + k^2/(2\mu)$, and approximating the momentum matrix element as $|P_{\mathrm{vc}}|^2 \simeq E_g/(4\mu)$. The time integral of the accumulated phase is expanded into a sum over multiphoton channels using 4-variable generalized Bessel functions $J_l(\eta_1,\eta_2,\eta_3,\eta_4)$; for nonzero relative phase, the non-vanishing oscillatory terms are collected into a 'twisted' multivariable Bessel function $J^\varphi_l$ whose extra arguments involve $\sin\varphi$ and $\cos 2\varphi$. Each channel is weighted by Bessel-index shifts $\pm 1$ (for the $\omega$ component), $\pm 2$ (for the $2\omega$ component), and their cross terms, while the Dirac delta $\delta(\xi_l)$ with $\xi_l = k^2/\mu + U_p + E_g + l\omega$ enforces energy conservation including the ponderomotive shift $U_p = (A_1^2 + A_2^2)/(4\mu)$. This structure is what yields both the channel-closing dips and the phase-dependent interference.

What would settle it

A decisive test is a relative-phase scan of the excitation probability in $\alpha$-quartz at equal 800 nm and 400 nm intensities near $I_{\mathrm{tot}}=10^{13}$ W/cm$^2$: the paper predicts a minimum near $\varphi \simeq 2\pi$ and a maximum near $\varphi \simeq 1.4\pi$. A measurement or a fully ab initio TDDFT calculation that places the extrema elsewhere, or reverses their ordering, would falsify the central claim. Because the predicted channel-closing dips are washed out by pulse bandwidth, the intensity-dependence part of the claim should instead be tested with narrow-band or long-pulse fields.

Watch

Extended reading notes

Core claim

The central claim is that the analytical transition probability, Eq. (30), captures the essential physics of two-color strong-field ionization in a real dielectric. When applied to $\alpha$-quartz with $\hbar\omega=1.55$ eV, equal $\omega$ and $2\omega$ intensities, $E_g=9$ eV and $\mu=0.28$ a.u., the formula reproduces the TDDFT results within roughly a factor of three and, more importantly, tracks the same qualitative trends: at total intensities above about $10^{13}$ W/cm$^2$ the two-color field excites more electrons than either single-color field, while at lower intensities the $2\omega$ field dominates; the relative-phase scan shows a minimum at $\varphi\simeq 2\pi$ and a maximum at $\varphi\simeq 1.4\pi$. The authors interpret the agreement as evidence that band dispersion near the gap, rather than detailed band-structure features, determines the dominant ionization dynamics in strong fields, while attributing residual discrepancies to the parabolic, isotropic band assumption and the lack of Brillouin-zone periodicity.

Load-bearing premise

The entire formula rests on treating the material as a parabolic, isotropic two-band system with one reduced mass and a momentum matrix element fixed by the band gap; if a target's real bands are strongly nonparabolic or anisotropic, the predicted intensity and phase dependence can shift or fail.

Editorial extensions

If this is right

The formula predicts that above roughly $10^{13}$ W/cm$^2$ an equally mixed $\omega+2\omega$ field excites more electrons than either frequency alone, matching the TDDFT enhancement and offering a target for two-color machining experiments.
It gives an explicit relative-phase dependence with a minimum near $\varphi \simeq 2\pi$ and a maximum near $\varphi \simeq 1.4\pi$ at $I_{\mathrm{tot}}=10^{13}$ W/cm$^2$, so phase-locked two-color pulses can be chosen to maximize or minimize excitation.
Channel-closing dips appear at specific intensities in the analytical scan (for example near $1.5\times10^{13}$ W/cm$^2$ for the $\omega+2\omega$ case in $\alpha$-quartz), marking where the ponderomotive shift pushes a multiphoton channel out of resonance.
Because the formula is analytic, parameter scans over $E_g$, $\mu$, intensity ratio, and phase can be evaluated almost instantly, making systematic material surveys feasible without TDDFT.
The cross terms between the $\omega$ and $2\omega$ amplitudes are isolated as the source of two-color mixing, explaining mechanistically why two-color excitation exceeds the sum of the single-color rates.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The channel-closing dips rely on the monochromatic, continuous-wave assumption, so finite-bandwidth experimental pulses will smear them; the qualitative phase dependence is the more robust prediction to test.
The same Houston-basis machinery could be pushed to non-collinear polarization geometries or to three-color fields, though the order of the multivariable Bessel functions would grow; the paper does not carry out that extension.
If Eq. (30) remains accurate in other wide-gap materials, two-color phase control could be used to steer ionization thresholds, damage morphology, or LIPSS formation without expensive first-principles runs; this practical consequence is left implicit.
A sharper stress test would apply the formula to a narrow-gap, strongly nonparabolic material such as silicon, where the authors' own error attribution predicts the quantitative agreement should degrade.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Two things to know. The paper extends the single-color Houston/Bessel formalism to a collinear ω+2ω field, and the φ=0 derivation in Sec. II.C.1 is explicit and internally consistent. The benchmark against TDDFT for α-quartz is honest: it shows a factor-of-three offset, and the authors flag the pulse-bandwidth and band-structure assumptions. That part is a useful, if approximate, analytical tool.

The problem is the φ≠0 result. Equations (22)-(24) define a 'twisted' Bessel function using real-argument modified Bessel functions I_n. But Eq. (22) expands a pure phase exp(i∫Δ dt): the generating function Σ_l J^φ_l(z) z^l must have modulus |z|=1. Modified Bessel functions with real arguments generate exp(η cos nωt), which is a real exponential, not an oscillatory phase. The correct Jacobi-Anger coefficient for exp(i u cos mωt) is i^n J_n(u), not I_n(u). So Eq. (30) does not follow from Eq. (22) as printed, and the phase-dependence comparison in Fig. 2 is built on an invalid identity. The φ=0 limit is safe because η'_5–η'_8 vanish, but the central qualitative claim about relative-phase control is currently unsupported.

I checked the stress-test note against the manuscript and it holds up. This is not a matter of taste; the equations are internally inconsistent. The factor-of-three intensity offset and the continuous-wave versus pulsed-field issue are secondary and well acknowledged. The k-resolved analysis in Figs. 3-4 is a good idea.

The paper is for people who want a quick analytical estimate of two-color excitation trends. It won't replace TDDFT for quantitative work, but if the φ≠0 identity is corrected, it could be a serviceable design tool. As it stands, only the φ=0 formula is usable.

Recommendation: send it to peer review, but demand a derivation or verification of Eq. (24) before acceptance. This is a serious within-subfield contribution that needs a substantive revision, not a desk reject.

Referee Report

4 major / 5 minor

Summary. The manuscript derives an analytical expression for the electron excitation probability in a parabolic two-band dielectric driven by a collinear two-color (ω+2ω) laser field, starting from the Houston-basis form of the time-dependent Schrödinger equation. The zero-relative-phase (φ=0) case is derived explicitly in terms of multivariable Bessel functions, while the nonzero-phase case is handled by introducing a 'twisted multivariable Bessel function.' The resulting formula is benchmarked against real-time TDDFT calculations for α-quartz for both the intensity dependence and the relative-phase dependence. The authors report qualitative agreement, including enhanced two-color excitation at high intensity and a phase dependence with extrema near φ≈1.4π and φ≈2π.

Significance. If the formulas were correct, this paper would provide a compact analytical tool for two-color strong-field excitation in dielectrics, complementing costly first-principles simulations and extending the single-color theory of JPSJ 88, 024706 (2019). The φ=0 derivation is explicit and internally coherent, and the TDDFT benchmarking is a useful attempt to validate the model against an independent method. The paper is also honest about the factor-of-three quantitative discrepancy and about the limitations of the parabolic two-band approximation. However, the φ≠0 part, which is the main novelty of the paper, contains mathematical errors that currently invalidate the reported relative-phase dependence. The paper has the right scope and a plausible framework, but the load-bearing equations need to be re-derived and the numerical results re-examined before the central claim can be accepted.

major comments (4)

[Sec. II.C.2, Eqs. (22)–(24)] The 'twisted' Bessel expansion is not a valid Fourier expansion of the pure phase factor exp(i∫Δ dt). A real-argument modified Bessel function I_n has generating function exp(z cos x), not exp(i z cos x), so inserting I_n(η'_6), I_n(η'_7), I_n(η'_8) into Eq. (24) cannot reproduce the oscillatory phase. Concretely, consider A1=0, kcosθ=0, and φ=π/4; then the phase contains the term i η'_8 cos4ωt with η'_8 = A2²/(16μω). The exact coefficient of e^{i4ωt} at first order in η'_8 is +i η'_8/2, whereas Eq. (24) gives −i I_1(η'_8) ≈ −i η'_8/2 after using J_0(0,η'_2,0,0)=1. The printed expansion therefore violates unitarity (Σ_l |c_l|²=1 for a pure phase), and Eq. (22) is not the phase factor (13) for φ≠0. Since Eq. (30) and the phase-dependence results in Fig. 2 are generated from this expansion, the central phase-dependence claim is currently unsupported.
[Sec. II.C.2, Eqs. (25)–(29)] The coefficient definitions for the φ≠0 phase factor are inconsistent with the explicit φ=0 result. Generalizing the Appendix A derivation to φ≠0 gives the coefficient of sin2ωt in the exponent as kA2 cosθ/(2μω) cosφ + A1²/(8μω). Equation (26), however, reads η'_2 = A1²/(8μω) + A2²/(2μω) cosφ, which omits the kA2 cosφ term and introduces an A2² term that belongs to the 4ω component (via sin(4ωt+2φ)). Consequently, setting φ=0 in Eq. (30) does not reduce to Eq. (21), even after correcting the Bessel-function issue. This is a separate, concrete algebraic error that must be fixed.
[Sec. II.C.1 and Sec. II.C.2, Eqs. (21) and (30)] The prefactors of the φ=0 total rate in Eq. (21) and in the φ=0 limit of Eq. (30) differ by a factor of 4π³: Eq. (21) has π²|Pvc|²μ^{3/2}/(4√2), while Eq. (30) has |Pvc|²μ^{3/2}/(16√2π). Since Eq. (30) should reduce to Eq. (21) when φ=0, at least one of these prefactors is wrong. The absolute scale of the analytical curves in Figs. 1 and 2 depends on which expression is used, so this inconsistency must be resolved before the quantitative benchmark can be assessed.
[Sec. II.C.2, Eqs. (22)–(30)] The nonzero-phase result is asserted without derivation: unlike the φ=0 case, there is no appendix or intermediate step showing how Eq. (11) becomes Eqs. (22)–(24). Given the mathematical errors identified above, this omission is not merely a presentation issue. The authors should provide a complete derivation of the φ≠0 phase-factor expansion, including the generating-function identity for the 'twisted' Bessel function, or clearly state any additional approximation that could justify the printed form.

minor comments (5)

[Abstract and Sec. III.B] The claim of 'remarkable qualitative accuracy' overstates the agreement shown in Fig. 1, where the analytical rate differs from TDDFT by a factor of about three and the channel-closing structure is largely absent in the TDDFT data. Suggest tempering the wording to 'qualitative agreement' and quantifying the discrepancy in the abstract.
[Sec. II.B, Eq. (18)] The factor |Pvc| appears inside the second line after A2(...), which is dimensionally inconsistent; it should be factored out of the entire bracket as in the first line.
[Fig. 3 caption] The word 'pannel' appears three times and should be 'panel'.
[Sec. II.C.2, Eq. (23)] The 'twisted multivariable Bessel function' is introduced with no reference or generating-function identity; if it is a new construct, its convergence properties and reduction to known cases should be stated explicitly.
[Sec. III.B] The analytical model uses Eg=9 eV while the TDDFT optical gap is reported as 8.75 eV; the 0.25 eV difference can shift channel-closing intensities and should be listed as a source of quantitative discrepancy.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: Eq. (30) is derived from the TDSE with explicit approximations and benchmarked against independent TDDFT, with material parameters taken from prior literature rather than fitted to the two-color data.

full rationale

The derivation chain is self-contained: the paper starts from the time-dependent Schrodinger equation in the Houston basis (Eqs. 1-6), makes the two-band parabolic approximation (Eq. 10), and obtains transition probabilities via Fourier/Bessel expansions of the phase factor (Eqs. 13-30). The material parameters Eg = 9 eV and mu = 0.28 a.u. are stated as assumptions taken from the earlier JPSJ work [25], not optimized against the TDDFT results for two-color pulses. The TDDFT benchmark (Eqs. 31-36) is an independent ab initio calculation, and the paper reports qualitative agreement with known discrepancies due to band-structure and periodicity assumptions. Self-citation of [25] supplies the single-color formalism and parameter values, but it does not by itself force the new nonzero-phase result; no fitted parameter is relabeled as a prediction, and no uniqueness theorem is imported. The skeptical concern about Eq. (24) (modified Bessel functions appearing in a pure-phase expansion) is a potential mathematical error or missing verification in the derivation, which is a correctness issue, not a reduction of the prediction to its inputs by construction. Therefore no circular step is exhibited, and the circularity score is 0.

Assumptions & free parameters 2 free parameters · 6 assumptions · 0 invented entities

The derivation is a first-principles perturbative calculation with material parameters inherited from the literature; no new physical entity is introduced. The main unproved mathematical input is the twisted Bessel identity used for the relative-phase case.

free parameters (2)

Bandgap E_g = 9 eV
Chosen from the measured range 8.6-9.5 eV for α-SiO2; the TDDFT-computed gap is 8.75 eV. Not fitted to the benchmark data, but directly affects the resonance condition in Eq. (21).
Reduced mass μ = 0.28 a.u.
Taken from prior work [25]; sets the parabolic curvature and the phase-space factor μ^{3/2} in Eq. (21). Not fitted to the two-color TDDFT data.

assumptions (6)

domain assumption Parabolic two-band dispersion Δ_vc^k = E_g + k^2/(2μ) (Eq. 10).
Replaces the real α-quartz band structure; the authors argue it works best in the strong-field regime.
domain assumption Independent-electron approximation and dipole approximation (Eq. 1).
Standard for strong-field solid theories; many-body and non-dipole effects are not included.
domain assumption Population remains in valence band: |C_vv| ≈ 1 >> |C_vc| (Sec. II.B).
Weak-excitation assumption used to truncate the Houston-basis equations to a first-order transition amplitude (Eq. 6).
domain assumption Momentum matrix element is uniform in k-space: |P_vc|^2 ≈ E_g/(4μ) (Eq. 19).
Kane-type approximation; ignores k-dependence and interband coupling variation.
domain assumption Long-time limit with continuous-wave field; oscillatory terms and terms containing Δ_vc^{k+A(t)} are neglected (Sec. II.B, Eq. 6).
Converts finite-pulse TDDFT dynamics to a CW rate with δ-function energy conservation; the authors note this causes discrepancies at low intensity.
ad hoc to paper The phase factor for φ≠0 equals the asserted twisted multivariable Bessel expansion (Eqs. 22-24).
Stated without derivation in Sec. II.C.2, this identity underlies the phase-dependence result; it is not verified numerically within the paper.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Electron Excitation Probability in Dielectrics under Two-color Intense Laser Fields." pith.science (2026). https://pith.science/paper/WHK54JW2

@misc{pith2026250701023,
  author       = {Pith},
  title        = {Pith review of: Electron Excitation Probability in Dielectrics under Two-color Intense Laser Fields},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WHK54JW2}},
  note         = {Machine review of arXiv:2507.01023}
}

abstract

Two-color laser fields offer significantly enhanced control over electron excitation dynamics under ultrashort intense laser pulses compared to monochromatic fields. However, their strong nonlinearity necessitates computationally expensive first-principles calculations to accurately predict ionization dynamics. To overcome this challenge, we derive an analytical expression for the ionization rate in dielectrics subjected to intense two-color laser fields, refining the theoretical framework introduced in JPSJ {\bf 88}, 024706 (2019). By benchmarking our formula against first-principles calculations based on time-dependent density functional theory (TDDFT) for $\alpha$-quartz, we demonstrate that our model captures the essential physics of ionization dynamics with remarkable qualitative accuracy, despite employing certain approximations. This analytical approach not only provides deeper physical insight but also offers a computationally efficient alternative for predicting strong-field interactions in dielectrics.

Figures

Figures reproduced from arXiv: 2507.01023 by the authors.

**Figure 2.** FIG. 2. Relative-phase dependence of excitation probabili [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 3.** (b) and the dotted line in [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

Figures from the paper (1 more)

**Figure 4.** Figure 4: FIG. 4. Angle-resolved transition probability [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

56 extracted references · 54 canonical work pages

[1]

Zero relative phase φ = 0 First, let us examine the φ = 0 case. Equation (11) is transformed into, exp ( i ∫ t ∆k+A(t′) vc dt′ ) = ∑ l exp ( i [k2 2µ +Up +Eg +lω ] t ) Jl(η1,η 2,η 3,η 4), (13) where Up = (A2 1 +A2 2)/4µ is the ponderomotive energy, Jl(η1,η 2,η 3,η 4) stands for the l-th order 4-variable Bessel function [48], deﬁned by, Jl(η1,η 2,η 3,η 4) ...

work page
[2]

Non-zero relative phase φ ⁄= 0 Next, turn to the φ ⁄= 0 case in Eq. (12). The phase factor [Eq. (13)] is modiﬁed as, exp ( i ∫ t ∆k+A(t′) vc dt′ ) = ∑ l exp ( i [k2 2µ +Up +Eg +lω ] t ) ×J φ l (η′ 1,η ′ 2,η ′ 3,η ′ 4,η ′ 5,η ′ 6,η ′ 7,η ′ 8), (22) whereJ φ l is the l-th order “twisted” multivariable Bessel function deﬁned as, J φ l (η′ 1,η ′ 2,η ′ 3,η ′ 4...

work page
[3]

Temple and M

P. Temple and M. Soileau, Polarization charge model for laser-induced ripple patterns in dielectric materials , IEEE J. Quantum Electron. 17, 2067 (1981)

work page 1981
[4]

(30) Again, the ﬁrst two terms J φ l±1 and the others J φ l±2 cor- respond to excitation by ω and 2ω, respectively, whereas their cross-terms of them represent the two-color eﬀect

˜Jl′ (η′ 5,η ′ 6,η ′ 7,η ′ 8), (23) where ˜Jl′ is deﬁned as, ˜Jl′ (η′ 5,η ′ 6,η ′ 7,η ′ 8) = ∑ l6,l7,l8 i−l6−2l7−3l8Jl′−2l6−3l7−4l8 (η′ 5)Il6 (η′ 6)Il7 (η′ 7)Il8 (η′ 8), (24) and, η′ 1 = A1k cosθ µω + A1A2 2µω cosφ (25) η′ 2 = A2 1 8µω + A2 2µω cosφ (26) η′ 3 = A1A2 6µω cosφ,η ′ 4 = A2 2 16µω cos 2φ (27) η′ 5 = A1A2 2µω sinφ,η ′ 6 = A2k cosθ 2µω sinφ (28)...

work page
[5]

The results obtained by the analytical formula derived in the previous section [Fig

Laser intensity dependence Figure 1 illustrates the intensity dependence of the ex- citation probability for the φ = 0 case for monochro- matic and two-color laser ﬁelds. The results obtained by the analytical formula derived in the previous section [Fig. 1 (a)] reproduce the overall trend of the TDDFT results [Fig. 1 (b)] surprisingly well, albeit diﬀere...

work page
[6]

Figure 2 compares our analytical expression with TDDFT calcu- lation results for an equally mixed two-color laser ﬁeld

Relative phase dependence Let us next discuss relative phase φ dependence of ex- citation probability for Itot = 10 13W/cm2 case. Figure 2 compares our analytical expression with TDDFT calcu- lation results for an equally mixed two-color laser ﬁeld. Here, the domain of φ is set to [ π, 2π]. The vector po- tential waveform is invariant under the transforma...

work page
[7]

B. C. Stuart, M. D. Feit, A. M. Rubenchik, B. W. Shore, and M. D. Perry, Laser-induced damage in dielectrics with nanosecond to subpicosecond pulses, Phys. Rev. Lett. 74, 2248 (1995)

work page 1995
[8]

T. T. Luu, M. Garg, S. Y. Kruchinin, A. Moulet, M. T. Hassan, and E. Goulielmakis, Extreme ultraviolet high- harmonic spectroscopy of solids, Nature 521, 498 (2015)

work page 2015

Show all 56 references

[9]

R. R. Gattass and E. Mazur, Femtosecond laser micro- machining in transparent materials, Nat. Photonics 2, 219 (2008)

2008
[10]

Sugioka and Y

K. Sugioka and Y. Cheng, Femtosecond laser three- dimensional micro- and nanofabrication, Appl. Phys. 9 Rev. 1, 041303 (2014)

2014
[11]

Haug and S

H. Haug and S. Schmitt-Rink, Electron theory of the optical properties of laser-excited semiconductors, Prog . Quantum Electron. 9, 3 (1984)

1984
[12]

Golde, T

D. Golde, T. Meier, and S. W. Koch, High harmonics generated in semiconductor nanostructures by the cou- pled dynamics of optical inter- and intraband excitations, Phys. Rev. B Condens. Matter 77, 075330 (2008)

2008
[13]

Kaneshima, Y

K. Kaneshima, Y. Shinohara, K. Takeuchi, N. Ishii, K. Imasaka, T. Kaji, S. Ashihara, K. L. Ishikawa, and J. Itatani, Polarization-resolved study of high harmonics from bulk semiconductors, Phys. Rev. Lett. 120, 243903 (2018)

2018
[14]

Yue and M

L. Yue and M. B. Gaarde, Introduction to theory of high- harmonic generation in solids: tutorial, J. Opt. Soc. Am. B 39, 535 (2022)

2022
[15]

Iwata, K

J.-I. Iwata, K. Yabana, and G. F. Bertsch, Real-space computation of dynamic hyperpolarizabilities, J. Chem. Phys. 115, 8773 (2001)

2001
[16]

Otobe, M

T. Otobe, M. Yamagiwa, J.-I. Iwata, K. Yabana, T. Nakatsukasa, and G. F. Bertsch, First-principles elec- tron dynamics simulation for optical breakdown of di- electrics under an intense laser ﬁeld, Phys. Rev. B Con- dens. Matter 77, 165104 (2008)

2008
[17]

M. Noda, S. A. Sato, Y. Hirokawa, M. Uemoto, T. Takeuchi, S. Yamada, A. Yamada, Y. Shinohara, M. Yamaguchi, K. Iida, I. Floss, T. Otobe, K.-M. Lee, K. Ishimura, T. Boku, G. F. Bertsch, K. Nobusada, and K. Yabana, SALMON: Scalable ab-initio Light–Matter simulator for optics and...

2019
[18]

Tancogne-Dejean, M

N. Tancogne-Dejean, M. J. T. Oliveira, X. Andrade, H. Appel, C. H. Borca, G. Le Breton, F. Buchholz, A. Castro, S. Corni, A. A. Correa, U. De Giovannini, A. Delgado, F. G. Eich, J. Flick, G. Gil, A. Gomez, N. Helbig, H. H¨ ubener, R. Jest¨ adt, J. Jornet-Somoza, A. H. Larsen, ...

2020
[19]

Fennel, G

T. Fennel, G. F. Bertsch, and K.-H. Meiwes-Broer, Ion- ization dynamics of simple metal clusters in intense ﬁelds by the thomas-fermi-vlasov method, The European Phys- ical Journal D - Atomic, Molecular, Optical and Plasma Physics 29, 367 (2004)

2004
[20]

K¨ ohn, R

J. K¨ ohn, R. Redmer, K.-H. Meiwes-Broer, and T. Fennel, Non-resonant absorption enhancement in laser-excited simple metal clusters through electron-electron collisio ns, Phys. Rev. A 77, 033202 (2008)

2008
[21]

Fennel, K.-H

T. Fennel, K.-H. Meiwes-Broer, J. Tiggesb¨ aumker, P.- G. Reinhard, P. M. Dinh, and E. Suraud, Laser-driven nonlinear cluster dynamics, Rev. Mod. Phys. 82, 1793 (2010)

2010
[22]

M. Tani, T. Otobe, Y. Shinohara, and K. L. Ishikawa, Semiclassical description of electron dynamics in ex- tended systems under intense laser ﬁelds, Phys. Rev. B Condens. Matter 104, 075157 (2021)

2021
[23]

L. V. Keldysh and Others, Diagram technique for nonequilibrium processes, Sov. Phys. JETP 20, 1018 (1965)

1965
[24]

Lewenstein, P

M. Lewenstein, P. Balcou, M. Y. Ivanov, A. L’Huillier, and P. B. Corkum, Theory of high-harmonic generation by low-frequency laser ﬁelds, Phys. Rev. A 49, 2117 (1994)

1994
[25]

S. M. Golin, S. E. Kirkwood, D. D. Klug, D. M. Vil- leneuve, D. M. Rayner, C. A. T. Herrero, and P. B. Corkum, Strong ﬁeld processes inside gallium arsenide, J. Phys. B At. Mol. Opt. Phys. 47, 204025 (2014)

2014
[26]

Koz´ ak, T

M. Koz´ ak, T. Otobe, M. Zukerstein, F. Troj´ anek, and P. Mal´ y, Anisotropy and polarization dependence of mul- tiphoton charge carrier generation rate in diamond, Phys. Rev. B. 99, 104305 (2019)

2019
[27]

F. H. M. Faisal, Multiple absorption of laser photons by atoms, J. Phys., B Proc. Phys. Soc., At. Mol. Phys. 6, L89 (1973)

1973
[28]

H. D. Jones and H. R. Reiss, Intense-ﬁeld eﬀects in solid s, Phys. Rev. 16, 2466 (1977)

1977
[29]

H. R. Reiss, Eﬀect of an intense electromagnetic ﬁeld on a weakly bound system, Phys. Rev. A Gen. Phys. 22, 1786 (1980)

1980
[30]

Otobe, Y

T. Otobe, Y. Shinohara, S. A. Sato, and K. Yabana, The- ory for electron excitation in dielectrics under an intense linear and circularly polarized laser ﬁelds, J. Phys. Soc. Jpn. 88, 024706 (2019)

2019
[31]

Venkat and T

P. Venkat and T. Otobe, Electron excitation rate in di- electrics under an intense elliptically polarized laser ﬁe ld, J. Phys. Soc. Jpn. 91 (2022)

2022
[32]

V. E. Gruzdev, Photoionization rate in wide band-gap crystals, Phys. Rev. B Condens. Matter Mater. Phys. 75 (2007)

2007
[33]

H. Taya, T. Fujimori, T. Misumi, M. Nitta, and N. Sakai, Exact wkb analysis of the vacuum pair production by time-dependent electric ﬁelds, Phys. Rev. A 2021, 82 (2021)

2021
[34]

Sugioka, S

K. Sugioka, S. Wada, A. Tsunemi, T. Sakai, H. Takai, H. Moriwaki, A. Nakamura, H. T. H. Tashiro, and K. T. K. Toyoda, Micropatterning of quartz substrates by multi-WavelengthVacuum-ultraviolet laser ablation, Jpn. J. Appl. Phys. 32, 6185 (1993)

1993
[35]

X. Yu, Q. Bian, Z. Chang, P. B. Corkum, and S. Lei, Femtosecond laser nanomachining initiated by ultraviolet multiphoton ionization, Opt. Express 21, 24185 (2013)

2013
[36]

Gedvilas, J

M. Gedvilas, J. Mikˇ sys, J. Berzinˇ s, V. Stankeviˇ c, an d G. Raˇ ciukaitis, Multi-photon absorption enhancement by dual-wavelength double-pulse laser irradiation for eﬃ- cient dicing of sapphire wafers, Sci. Rep. 7, 5218 (2017)

2017
[37]

V. P. Zhukov and N. M. Bulgakova, Volumetric mod- iﬁcation of transparent materials with Two-Color laser irradiation: Insight from numerical modeling, Materials 17 (2024)

2024
[38]

H. Yao, M. Chen, L. Jiang, B. Guo, Z. Yu, G. Wang, and S. Song, Selective processing of CFRP by two-color double-pulse femtosecond laser, Opt. Express 33, 3101 (2025)

2025
[39]

Ishikawa, Photoemission and ionization of he+ un- der simultaneous irradiation of fundamental laser and high-order harmonic pulses, Phys

K. Ishikawa, Photoemission and ionization of he+ un- der simultaneous irradiation of fundamental laser and high-order harmonic pulses, Phys. Rev. Lett. 91, 043002 (2003)

2003
[40]

E. J. Takahashi, T. Kanai, K. L. Ishikawa, Y. Nabekawa, and K. Midorikawa, Dramatic enhancement of high-order harmonic generation, Phys. Rev. Lett. 99, 053904 (2007)

2007
[41]

X. Song, S. Yang, R. Zuo, T. Meier, and W. Yang, En- hanced high-order harmonic generation in semiconduc- tors by excitation with multicolor pulses, Phys. Rev. A 10 101, 033410 (2020)

2020
[42]

Duchateau, A

G. Duchateau, A. Yamada, and K. Yabana, Electron dy- namics in α-quartz induced by two-color 10-femtosecond laser pulses, Phys. Rev. B Condens. Matter 105, 165128 (2022)

2022
[43]

M. Tani, K. Sasaki, Y. Shinohara, and K. L. Ishikawa, Enhanced energy deposition and carrier generation in silicon induced by two-color intense femtosecond laser pulses, Phys. Rev. B Condens. Matter 106, 195141 (2022)

2022
[44]

Sheehy, B

B. Sheehy, B. Walker, and L. F. DiMauro, Phase control in the two-color photodissociation of HD+, Phys. Rev. Lett. 74, 4799 (1995)

1995
[45]

T. Endo, H. Fujise, Y. Kawachi, A. Ishihara, A. Mat- suda, M. Fushitani, H. Kono, and A. Hishikawa, Selective bond breaking of CO2 in phase-locked two-color intense laser ﬁelds: laser ﬁeld intensity dependence, Phys. Chem. Chem. Phys. 19, 3550 (2017)

2017
[46]

Di Fraia, O

M. Di Fraia, O. Plekan, C. Callegari, K. C. Prince, L. Gi- annessi, E. Allaria, L. Badano, G. De Ninno, M. Trov` o, B. Diviacco, D. Gauthier, N. Mirian, G. Penco, P. R. Ribiˇ c, S. Spampinati, C. Spezzani, G. Gaio, Y. Orimo, O. Tugs, T. Sato, K. L. Ishikawa, P. A. Carpeggiani,...

2019
[47]

D. You, K. Ueda, E. V. Gryzlova, A. N. Grum- Grzhimailo, M. M. Popova, E. I. Staroselskaya, O. Tugs, Y. Orimo, T. Sato, K. L. Ishikawa, P. A. Carpeg- giani, T. Csizmadia, M. F¨ ule, G. Sansone, P. K. Maroju, A. D’Elia, T. Mazza, M. Meyer, C. Callegari, M. Di Fraia, O. Plekan, ...

2020
[48]

Orimo, O

Y. Orimo, O. Tugs, T. Sato, D. You, K. Ueda, and K. L. Ishikawa, Interferometric extraction of photoionization - path amplitudes and phases from time-dependent multi- conﬁguration self-consistent-ﬁeld simulations, J. Phys. B At. Mol. Opt. Phys. 54, 074001 (2021)

2021
[49]

E. V. Gryzlova, P. Carpeggiani, M. M. Popova, M. D. Kiselev, N. Douguet, M. Reduzzi, M. Negro, A. Comby, H. Ahmadi, V. Wanie, M. C. Castrovilli, A. Fischer, P. Eng-Johnsson, M. Meyer, K. Bartschat, S. M. Burkov, T. Csizmadia, M. Dumergue, S. K¨ uhn, N. G. Harshitha, M. Fule, F...

2022
[50]

K. L. Ishikawa, K. C. Prince, and K. Ueda, Control of ion-photoelectron entanglement and coherence via rabi oscillations, J. Phys. Chem. A 127, 10638 (2023)

2023
[51]

Runge and E

E. Runge and E. K. U. Gross, Density-functional the- ory for time-dependent systems, Phys. Rev. Lett. 52, 997 (1984)

1984
[52]

Yabana and G

K. Yabana and G. F. Bertsch, Time-dependent local- density approximation in real time, Phys. Rev. B Con- dens. Matter 54, 4484 (1996)

1996
[53]

Dattoli, C

G. Dattoli, C. Chiccoli, S. Lorenzutta, G. Maino, M. Richetta, and A. Torre, Generating functions of mul- tivariable generalized bessel functions and jacobi-ellip tic functions, J. Math. Phys. 33, 25 (1992)

1992
[54]

Tran and P

F. Tran and P. Blaha, Accurate band gaps of semiconduc- tors and insulators with a semilocal exchange-correlation potential, Phys. Rev. Lett. 102, 226401 (2009)

2009
[55]

F. Bart, M. Gautier, J. P. Duraud, and M. Henriot, (011
[56]

α-quartz surface: a LEED, XANES and ELS study, Surf. Sci. 274, 317 (1992)

1992

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [1]

Zero relative phase φ = 0 First, let us examine the φ = 0 case. Equation (11) is transformed into, exp ( i ∫ t ∆k+A(t′) vc dt′ ) = ∑ l exp ( i [k2 2µ +Up +Eg +lω ] t ) Jl(η1,η 2,η 3,η 4), (13) where Up = (A2 1 +A2 2)/4µ is the ponderomotive energy, Jl(η1,η 2,η 3,η 4) stands for the l-th order 4-variable Bessel function [48], deﬁned by, Jl(η1,η 2,η 3,η 4) ...

work page

[2] [2]

Non-zero relative phase φ ⁄= 0 Next, turn to the φ ⁄= 0 case in Eq. (12). The phase factor [Eq. (13)] is modiﬁed as, exp ( i ∫ t ∆k+A(t′) vc dt′ ) = ∑ l exp ( i [k2 2µ +Up +Eg +lω ] t ) ×J φ l (η′ 1,η ′ 2,η ′ 3,η ′ 4,η ′ 5,η ′ 6,η ′ 7,η ′ 8), (22) whereJ φ l is the l-th order “twisted” multivariable Bessel function deﬁned as, J φ l (η′ 1,η ′ 2,η ′ 3,η ′ 4...

work page

[3] [3]

Temple and M

P. Temple and M. Soileau, Polarization charge model for laser-induced ripple patterns in dielectric materials , IEEE J. Quantum Electron. 17, 2067 (1981)

work page 1981

[4] [4]

(30) Again, the ﬁrst two terms J φ l±1 and the others J φ l±2 cor- respond to excitation by ω and 2ω, respectively, whereas their cross-terms of them represent the two-color eﬀect

˜Jl′ (η′ 5,η ′ 6,η ′ 7,η ′ 8), (23) where ˜Jl′ is deﬁned as, ˜Jl′ (η′ 5,η ′ 6,η ′ 7,η ′ 8) = ∑ l6,l7,l8 i−l6−2l7−3l8Jl′−2l6−3l7−4l8 (η′ 5)Il6 (η′ 6)Il7 (η′ 7)Il8 (η′ 8), (24) and, η′ 1 = A1k cosθ µω + A1A2 2µω cosφ (25) η′ 2 = A2 1 8µω + A2 2µω cosφ (26) η′ 3 = A1A2 6µω cosφ,η ′ 4 = A2 2 16µω cos 2φ (27) η′ 5 = A1A2 2µω sinφ,η ′ 6 = A2k cosθ 2µω sinφ (28)...

work page

[5] [5]

The results obtained by the analytical formula derived in the previous section [Fig

Laser intensity dependence Figure 1 illustrates the intensity dependence of the ex- citation probability for the φ = 0 case for monochro- matic and two-color laser ﬁelds. The results obtained by the analytical formula derived in the previous section [Fig. 1 (a)] reproduce the overall trend of the TDDFT results [Fig. 1 (b)] surprisingly well, albeit diﬀere...

work page

[6] [6]

Figure 2 compares our analytical expression with TDDFT calcu- lation results for an equally mixed two-color laser ﬁeld

Relative phase dependence Let us next discuss relative phase φ dependence of ex- citation probability for Itot = 10 13W/cm2 case. Figure 2 compares our analytical expression with TDDFT calcu- lation results for an equally mixed two-color laser ﬁeld. Here, the domain of φ is set to [ π, 2π]. The vector po- tential waveform is invariant under the transforma...

work page

[7] [7]

B. C. Stuart, M. D. Feit, A. M. Rubenchik, B. W. Shore, and M. D. Perry, Laser-induced damage in dielectrics with nanosecond to subpicosecond pulses, Phys. Rev. Lett. 74, 2248 (1995)

work page 1995

[8] [8]

T. T. Luu, M. Garg, S. Y. Kruchinin, A. Moulet, M. T. Hassan, and E. Goulielmakis, Extreme ultraviolet high- harmonic spectroscopy of solids, Nature 521, 498 (2015)

work page 2015