pith. sign in

arxiv: 2604.25093 · v1 · submitted 2026-04-28 · ❄️ cond-mat.mtrl-sci

Ultrafast Energy Absorption in Silicon Controlled by Two-Color Double Pulses

Eiyu S. Gushiken , Mizuki Tani , Hiroki Katow , Kenichi L. Ishikawa This is my paper

Pith reviewed 2026-05-07 16:14 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords siliconenergy absorptiontwo-color pulsesfemtosecond lasersultrafast dynamicssemiconductorsnonequilibrium stateslaser control
0
0 comments X

The pith

Two-color double pulses control how much energy silicon absorbs depending on intensity and pulse order.

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

The paper shows that energy absorption in crystalline silicon can be tuned by sending two femtosecond laser pulses of different wavelengths one after the other. A sympathetic reader would care because this offers a way to manage ultrafast energy transfer in a common semiconductor simply by choosing the wavelengths, intensities, and sequence of the pulses. Simulations across low to high intensities reveal three regimes: shorter wavelengths absorb more at low power through multiphoton processes, longer wavelengths work better at high power where tunneling dominates, and in the middle range a short pulse first followed by a long one boosts total absorption by raising the energy each electron gains from the second pulse.

Core claim

Energy absorption in crystalline silicon can be controlled by two-color femtosecond double-pulse irradiation, in which two temporally separated pulses with different wavelengths interact sequentially with the system. The governing mechanism and optimal wavelength combination depend on intensity: multiphoton interband absorption dominates at low intensity with shorter wavelengths preferred, while tunneling ionization and intraband acceleration favor longer wavelengths at high intensity. In the intermediate regime, a short-wavelength pulse preceding a long one enhances absorption because the nonequilibrium state left by the first pulse increases the energy gain per excited electron from the se

What carries the argument

The nonequilibrium electronic state prepared by the first pulse, which alters the excitation process and energy gain per carrier induced by the second pulse of different wavelength.

If this is right

  • At low intensities, pairs of shorter wavelengths enhance absorption via multiphoton interband transitions.
  • At high intensities, pairs of longer wavelengths enhance absorption through increased tunneling ionization and intraband acceleration.
  • At intermediate intensities, a short-wavelength pulse followed by a long-wavelength pulse maximizes absorption by raising energy gain per excited electron.
  • The total absorbed energy is governed by both the number of excited carriers and the energy each gains, both of which can be adjusted by pulse wavelength and order.
  • This tunability holds over the examined range of peak intensities from 2 times 10 to the 11 to 10 to the 13 watts per square centimeter and wavelengths of 515, 1030, and 2060 nanometers.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The pulse-sequence control could be applied to design more precise laser protocols for modifying silicon surfaces or creating defects at specific depths.
  • Similar intensity-dependent regime shifts might appear in other group-IV semiconductors, suggesting the effect is not unique to silicon.
  • Extending the approach to shaped pulses or more than two colors could further refine the energy deposition for applications like optical data storage or ultrafast switching.

Load-bearing premise

The simulations correctly capture the electron behavior and energy absorption in silicon for these laser conditions without major errors from the model approximations.

What would settle it

An experiment measuring absorbed energy for short-first versus long-first pulse pairs at intermediate intensities and finding no extra absorption for the short-first order would show the claimed enhancement is absent.

Figures

Figures reproduced from arXiv: 2604.25093 by Eiyu S. Gushiken, Hiroki Katow, Kenichi L. Ishikawa, Mizuki Tani.

Figure 1
Figure 1. Figure 1: FIG. 1. Temporal profile of the electric field E( view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Temporal profile of the calculated number of view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Mean absorbed energy per excited electron for each view at source ↗
Figure 5
Figure 5. Figure 5: shows snapshots of the DDOS at t = 35 and 65 fs, corresponding to the times immediately af￾ter the first and second laser pulses, respectively, for I1 = I2 = 3.5 × 1012 W/cm2 , Tdelay = 35 fs, for the 515+2060 nm [panel (a)] and 2060+515 nm [panel (b)]. We first discuss the 515+2060 nm case shown in view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Number of electrons excited to states within and be view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Relative phase dependence of the absorbed energy view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. (a) Temporal profile of the calculated number of ex view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Mean absorbed energy per excited electron for each view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. (b) Total number of electrons excited to energy view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. (a) Temporal profile and (b) final number of excited view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. Mean absorbed energy per excited electron for each view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. DDOS at view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17. (a) DDOS at the end of the simulation and (b) view at source ↗
Figure 19
Figure 19. Figure 19: FIG. 19. The calculated density of states of Si. 0 eV is set to view at source ↗
Figure 20
Figure 20. Figure 20: FIG. 20. The calculated dielectric function of crystalline Si. view at source ↗
read the original abstract

We theoretically show that energy absorption in crystalline silicon can be controlled by two-color femtosecond double-pulse irradiation, in which two temporally separated pulses with different wavelengths interact sequentially with the system. Using time-dependent density functional theory, we systematically examine the wavelength and intensity dependence of the absorbed energy over peak intensities of $2\times10^{11}$-$10^{13}$ W/cm$^2$ and wavelengths of 515, 1030, and 2060 nm. We find that the mechanism governing energy absorption and the optimal wavelength combination strongly depend on the intensity regime. In the low-intensity regime, multiphoton interband absorption is dominant, and energy absorption is enhanced for pulse pairs composed of shorter wavelengths. In contrast, in the high-intensity regime, the contributions of tunneling ionization and intraband acceleration become significant, leading to enhanced absorption for longer-wavelength combinations. In the intermediate-intensity regime, a pronounced enhancement is observed when a short-wavelength pulse precedes a long-wavelength pulse. Our analysis reveals that the nonequilibrium electronic state prepared by the first pulse modifies the excitation process induced by the second pulse, thereby enhancing the absorbed energy through an increased energy gain per excited electron. In this regime, the energy absorption is governed not only by the number of excited carriers but also by the energy gain per excited electron, which can be strongly modified by the pulse sequence. These results indicate that ultrafast energy transfer in semiconductors is tunable by appropriately designing the wavelength and intensity combination of the two pulses, and provide microscopic insight into two-color strong-field excitation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript uses time-dependent density functional theory (TDDFT) simulations to show that energy absorption in crystalline silicon can be controlled via two-color femtosecond double-pulse irradiation. It reports a systematic scan over wavelengths (515, 1030, 2060 nm) and peak intensities (2×10^11 to 10^13 W/cm²), finding intensity-dependent optimal pulse sequences: shorter-wavelength pairs enhance absorption at low intensity via multiphoton interband processes, longer-wavelength pairs at high intensity via tunneling and intraband acceleration, and short-then-long sequences at intermediate intensity due to nonequilibrium states prepared by the first pulse that increase energy gain per excited electron.

Significance. If the simulation results are robust, the work delivers microscopic insight into tunable nonequilibrium excitation in semiconductors under strong fields, particularly the role of pulse sequencing in modifying carrier energies beyond simple density changes. The forward, parameter-free nature of the TDDFT scan (no fitting to prior outputs) and the clear regime partitioning provide a useful theoretical framework for designing ultrafast laser interactions, though quantitative reliability depends on addressing standard TDDFT limitations.

major comments (3)
  1. [Methods section] Methods section: The calculations rely on standard adiabatic TDDFT (likely LDA/GGA) without any reported scissor correction, hybrid-functional benchmark, or explicit test of bandgap sensitivity. This directly affects the low- and intermediate-intensity regime boundaries, as the ~0.5 eV Si gap underestimate alters multiphoton absorption thresholds and the intensity at which tunneling/intraband contributions dominate, undermining the claimed wavelength optima and the 'energy gain per excited electron' mechanism.
  2. [Results section (intensity-regime analysis)] Results section (intensity-regime analysis): No convergence tests are provided for critical numerical parameters (k-point sampling, real-space grid, time step, or supercell size) despite quantitative claims on absorbed-energy ratios and enhancements across three orders of magnitude in intensity. Without these, it is unclear whether the reported differences between pulse sequences are numerically converged or sensitive to discretization choices.
  3. [Abstract and regime-partitioning discussion] Abstract and regime-partitioning discussion: The central claim partitions behavior into low-, intermediate-, and high-intensity regimes with distinct optimal combinations, but the manuscript lacks any direct comparison to experiment or higher-level methods (e.g., TDDFT with meta-GGA or real-time GW). This leaves the attribution of mechanisms and the intermediate-regime enhancement vulnerable to functional errors.
minor comments (2)
  1. [Figures] Figure captions and legends should explicitly state the pulse ordering (e.g., '515 nm first, 2060 nm second') and the definition of absorbed energy per unit volume for clarity.
  2. [Discussion] The phrase 'energy gain per excited electron' is used repeatedly but never given an explicit formula; adding an equation (e.g., E_abs / N_excited) would improve reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments, which have helped us identify areas for improvement. We address each major comment point by point below, providing our responses and indicating planned revisions to the manuscript.

read point-by-point responses

  1. Referee: [Methods section] The calculations rely on standard adiabatic TDDFT (likely LDA/GGA) without any reported scissor correction, hybrid-functional benchmark, or explicit test of bandgap sensitivity. This directly affects the low- and intermediate-intensity regime boundaries, as the ~0.5 eV Si gap underestimate alters multiphoton absorption thresholds and the intensity at which tunneling/intraband contributions dominate, undermining the claimed wavelength optima and the 'energy gain per excited electron' mechanism.

    Authors: We acknowledge that the LDA bandgap underestimate (~0.5 eV) is a standard limitation and can shift absolute multiphoton thresholds. Our analysis, however, emphasizes relative trends and mechanisms across pulse sequences rather than precise boundary values. We will revise the Methods section to include a brief discussion of this effect and add results from scissor-corrected calculations for the low-intensity regime, confirming that the preference for shorter-wavelength pairs persists. In the intermediate regime, the per-electron energy gain arises primarily from intraband dynamics in the nonequilibrium state, which remains robust to moderate gap shifts. We will also reference prior TDDFT studies on silicon to contextualize the functional choice. revision: partial

  2. Referee: [Results section (intensity-regime analysis)] No convergence tests are provided for critical numerical parameters (k-point sampling, real-space grid, time step, or supercell size) despite quantitative claims on absorbed-energy ratios and enhancements across three orders of magnitude in intensity. Without these, it is unclear whether the reported differences between pulse sequences are numerically converged or sensitive to discretization choices.

    Authors: We thank the referee for highlighting this omission. Our simulations employed a 4×4×4 k-grid, 0.3 a.u. real-space spacing, 0.02 a.u. time step, and 8-atom supercell, with internal checks showing absorbed-energy variations below 5% under moderate parameter changes. To address the concern transparently, we will add a dedicated paragraph in the Methods or Results section (with supporting data in the SI) demonstrating convergence of the key ratios and sequence-dependent enhancements with respect to these parameters. revision: yes

  3. Referee: [Abstract and regime-partitioning discussion] The central claim partitions behavior into low-, intermediate-, and high-intensity regimes with distinct optimal combinations, but the manuscript lacks any direct comparison to experiment or higher-level methods (e.g., TDDFT with meta-GGA or real-time GW). This leaves the attribution of mechanisms and the intermediate-regime enhancement vulnerable to functional errors.

    Authors: We agree that higher-level benchmarks would be desirable, but real-time GW or meta-GGA scans over the full intensity-wavelength grid remain computationally prohibitive. Our TDDFT approach follows established practice for strong-field semiconductor dynamics, and the identified mechanisms (multiphoton, tunneling, intraband) align with prior literature. The intermediate-regime enhancement is substantiated by explicit time-dependent occupation and kinetic-energy analyses showing pulse-order effects on per-carrier energy gain. We will expand the discussion section to cite relevant experimental two-color studies, explicitly note functional limitations, and qualify the regime boundaries as TDDFT-specific while preserving the qualitative insights. revision: partial

Circularity Check

0 steps flagged

No circularity: results from direct TDDFT forward simulations

full rationale

The paper's claims derive from solving the time-dependent Kohn-Sham equations under two-color laser fields for crystalline silicon, computing absorbed energy as a function of intensity and wavelength combinations. All reported regime boundaries, optimal sequences, and mechanisms (multiphoton dominance at low I, tunneling/intraband at high I, nonequilibrium modification at intermediate I) are outputs of these explicit propagations rather than inputs. No parameters are fitted to data subsets and renamed as predictions, no self-definitional loops exist in the equations, and no load-bearing self-citations reduce the central result to prior author outputs. The derivation chain remains independent of the paper's own results.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The study rests on the domain assumption that TDDFT faithfully reproduces strong-field electron dynamics in silicon; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Time-dependent density functional theory provides a sufficiently accurate description of electron dynamics in solids under intense laser fields for the purposes of this study.
    All reported results are obtained from TDDFT simulations as stated in the abstract.

pith-pipeline@v0.9.0 · 5593 in / 1276 out tokens · 66948 ms · 2026-05-07T16:14:40.914389+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

72 extracted references · 72 canonical work pages

  1. [1]

    Fig- ure 2(a) shows the temporal evolution of the absorbed energy

    Energy transfer to electrons We now examine how the energy transfer depends on the combination of wavelengths under the fixed condition ofI 1 =I 2 = 3.5×10 12 W/cm2 andT delay = 35 fs. Fig- ure 2(a) shows the temporal evolution of the absorbed energy. The energy deposited by the first pulse (around 30 fs) is found to be nearly independent of its wave- len...

    work page 2060

  2. [2]

    Number of excited electrons and mean absorbed energy The absorbed energy can be decomposed into two con- tributions: the number of generated carriers (electrons promoted to the conduction band) and their mean energy absorption. Then, is the difference in energy deposition observed in the previous subsection primarily due to that in the former or the latte...

    work page 2060

  3. [3]

    To this end, we compute the differ- ence density of states (DDOS), which reflects changes in electron occupancy relative to the ground state

    Excitation dynamics analysis The findings in the previous subsection motivate a more detailed analysis of the time- and energy-resolved excitation behavior. To this end, we compute the differ- ence density of states (DDOS), which reflects changes in electron occupancy relative to the ground state. DDOS is calculated by projecting the time-dependent Kohn-S...

    work page 2060

  4. [4]

    A positive de- lay means that the 515 nm pulse precedes the 2060 nm pulse

    Delay-time and relative-phase dependence Figure 7 shows the total absorbed energy and the num- ber of excited electrons as functions of the delay time Tdelay atI 1 =I 2 = 3.5×10 12 W/cm2. A positive de- lay means that the 515 nm pulse precedes the 2060 nm pulse. When the pulses overlap (−20≲T delay ≲20), the energy absorption is increased, as we have prev...

    work page 2060

  5. [5]

    78 ∝ I 2.02 ∝ I 1.90 ∝ I 3.47 FIG. 18. Intensity dependence of total absorbed energy for se- lected wavelength combinations: 515+515 nm, 515+2060 nm, 2060+515 nm, and 2060+2060 nm. The time delay is fixed atT delay = 35 fs. The dashed lines represent power-law fits W∝I s obtained from the two lowest-intensity points. estimated using the impact ionization ...

    work page 2060

  6. [6]

    Ghimire, A

    S. Ghimire, A. D. DiChiara, E. Sistrunk, P. Agostini, L. F. DiMauro, and D. A. Reis, Observation of high-order harmonic generation in a bulk crystal, Nature Physics7, 138 (2011)

    work page 2011

  7. [7]

    Schubert, M

    O. Schubert, M. Hohenleutner, F. Langer, B. Urbanek, C. Lange, U. Huttner, D. Golde, T. Meier, M. Kira, S. W. Koch, and R. Huber, Sub-cycle control of terahertz high- harmonic generation by dynamical bloch oscillations, Na- ture Photonics8, 119 (2014)

    work page 2014

  8. [8]

    Hohenleutner, F

    M. Hohenleutner, F. Langer, O. Schubert, M. Knorr, U. Huttner, S. W. Koch, M. Kira, and R. Huber, Real- time observation of interfering crystal electrons in high- harmonic generation, Nature523, 572 (2015)

    work page 2015

  9. [9]

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

    work page 2015

  10. [10]

    Sugioka and Y

    K. Sugioka and Y. Cheng, Ultrafast lasers—reliable tools for advanced materials processing, Light: Science & Ap- plications3, e149 (2014)

    work page 2014

  11. [11]

    R. R. Gattass and E. Mazur, Femtosecond laser micro- machining in transparent materials, Nature photonics2, 219 (2008)

    work page 2008

  12. [12]

    Malinauskas, A

    M. Malinauskas, A. ˇZukauskas, S. Hasegawa, Y. Hayasaki, V. Mizeikis, R. Buividas, and S. Juodkazis, Ultrafast laser processing of materials: from science to industry, Light: Science & Applications5, e16133 (2016)

    work page 2016

  13. [13]

    C. Yang, C. Ji, S. Feng, Y. Liu, W. Wei, and Y. Long, Ultrafast laser-matter interaction mechanisms and appli- cations in functional device fabrication: Recent advances and perspectives, Applied Physics Reviews12, 031325 (2025)

    work page 2025

  14. [14]

    J. Ma, J. Wu, Z. Lin, J. Wang, W. Yao, Y. Zhang, X. Zhang, L. Zhu, Y. Hayasaki, and H. Zhang, Femtosecond-laser preparation of hydrogel with micro/nano-structures and their biomedical appli- cations, Small Science5, 2400400 (2025)

    work page 2025

  15. [15]

    B. N. Chichkov, C. Momma, S. Nolte, F. Von Al- vensleben, and A. T¨ unnermann, Femtosecond, picosec- ond and nanosecond laser ablation of solids, Applied physics A63, 109 (1996)

    work page 1996

  16. [16]

    Sokolowski-Tinten, J

    K. Sokolowski-Tinten, J. Bialkowski, A. Cavalleri, D. von der Linde, A. Oparin, J. Meyer-ter Vehn, and S. I. Anisimov, Transient states of matter during short pulse laser ablation, Phys. Rev. Lett.81, 224 (1998)

    work page 1998

  17. [17]

    B. C. Stuart, M. D. Feit, S. Herman, A. M. 13 Rubenchik, B. W. Shore, and M. D. Perry, Nanosecond- to-femtosecond laser-induced breakdown in dielectrics, Phys. Rev. B53, 1749 (1996)

    work page 1996

  18. [18]

    Kerse, H

    C. Kerse, H. Kalaycıo˘ glu, P. Elahi, B. C ¸ etin, D. K. Kesim, ¨O. Ak¸ caalan, S. Yava¸ s, M. D. A¸ sık, B.¨Oktem, H. Hoogland, R. Holzwarth, and F. ¨O. Ilday, Ablation- cooled material removal with ultrafast bursts of pulses, Nature537, 84 (2016)

    work page 2016

  19. [19]

    X. Liu, D. Du, and G. Mourou, Laser ablation and mi- cromachining with ultrashort laser pulses, IEEE Journal of Quantum Electronics33, 1706 (1997)

    work page 1997

  20. [20]

    Mishchik, G

    K. Mishchik, G. Bonamis, J. Qiao, J. Lopez, E. Au- douard, E. Mottay, C. H¨ onninger, and I. Manek- H¨ onninger, High-efficiency femtosecond ablation of sili- con with ghz repetition rate laser source, Opt. Lett.44, 2193 (2019)

    work page 2019

  21. [21]

    Bonamis, K

    G. Bonamis, K. Mishchick, E. Audouard, C. H¨ onninger, E. Mottay, J. Lopez, and I. Manek-H¨ onninger, High ef- ficiency femtosecond laser ablation with gigahertz level bursts, Journal of Laser Applications31, 022205 (2019)

    work page 2019

  22. [22]

    Caballero-Lucas, K

    F. Caballero-Lucas, K. Obata, and K. Sugioka, Enhanced ablation efficiency for silicon by femtosecond laser micro- processing with ghz bursts in mhz bursts(biburst), In- ternational Journal of Extreme Manufacturing4, 015103 (2022)

    work page 2022

  23. [23]

    Zubauskas, E

    L. Zubauskas, E. Markauskas, P. ˇSleini¯ ut˙ e, and P. Geˇ cys, Classical top-down fused silica milling with a femtosec- ond laser using different laser pulse burst regimes, Jour- nal of Manufacturing Processes155, 820 (2025)

    work page 2025

  24. [24]

    Gaudfrin, J

    K. Gaudfrin, J. Lopez, K. Mishchik, L. Gemini, R. Kling, and G. Duchateau, Fused silica ablation by double fem- tosecond laser pulses: influence of polarization state, Opt. Express28, 15189 (2020)

    work page 2020

  25. [25]

    Zukerstein, V

    M. Zukerstein, V. P. Zhukov, T. J.-Y. Derrien, O. Fedo- tova, and N. M. Bulgakova, Double-pulse-laser volumet- ric modification of fused silica: the effect of pulse delay on light propagation and energy deposition, Opt. Express 32, 12882 (2024)

    work page 2024

  26. [26]

    Xiaoming, C

    Y. Xiaoming, C. Zenghu, C. P. B., and L. Shuting, Fab- ricating nanostructures on fused silica using femtosecond infrared pulses combined with sub-nanojoule ultraviolet pulses, Opt. Lett.39, 5638 (2014)

    work page 2014

  27. [27]

    Gedvilas, J

    M. Gedvilas, J. Mikˇ sys, J. Berzinˇ s, V. Stankeviˇ c, and G. Raˇ ciukaitis, Multi-photon absorption enhancement by dual-wavelength double-pulse laser irradiation for ef- ficient dicing of sapphire wafers, Scientific Reports7, 10.1038/s41598-017-05548-x (2017)

    work page doi:10.1038/s41598-017-05548-x 2017

  28. [28]

    D. Chu, P. Yao, X. Sun, K. Yin, and C. Huang, Ablation enhancement of fused silica glass by femtosecond laser double-pulse bessel beam, J. Opt. Soc. Am. B37, 3535 (2020)

    work page 2020

  29. [29]

    S. I. Kudryashov, A. A. Samokhvalov, Y. D. Golubev, D. S. Ivanov, M. E. Garcia, V. P. Veiko, B. Rethfeld, and V. Y. Mikhailovskii, Dynamic all-optical control in ultrashort double-pulse laser ablation, Applied Surface Science537, 147940 (2021)

    work page 2021

  30. [30]

    Zhang, S

    J. Zhang, S. Wang, M. Wang, and Z. Chu, Femtosecond laser double pulses nanofabrication on silicon, IOP Con- ference Series: Materials Science and Engineering565, 012018 (2019)

    work page 2019

  31. [31]

    Zoppel, J

    S. Zoppel, J. Zehetner, and G. A. Reider, Two color laser ablation: Enhanced yield, improved machining, Applied Surface Science253, 7692 (2007), photon-Assisted Syn- thesis and Processing of Functional Materials

    work page 2007

  32. [32]

    Duchateau, A

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

    work page 2022

  33. [33]

    Y. Ito, R. Yoshizaki, N. Miyamoto, and N. Sugita, Ultra- fast and precision drilling of glass by selective absorption of fiber-laser pulse into femtosecond-laser-induced fila- ment, Applied Physics Letters113, 061101 (2018)

    work page 2018

  34. [34]

    Yoshizaki, Y

    R. Yoshizaki, Y. Ito, S. Yoshitake, C. Wei, A. Shibata, I. Nagasawa, K. Nagato, and N. Sugita, Mechanism of material removal through transient and selective laser ab- sorption into excited electrons in fused silica, Journal of Applied Physics130, 053102 (2021)

    work page 2021

  35. [35]

    Zhang, T

    Y. Zhang, T. Koike, R. Yoshizaki, G. Ren, A. Shi- bata, S. Kiriake, R. Hasegawa, I. N. nd Keisuke Nagato, N. Sugita, and Y. Ito, Ultrahigh-speed laser drilling of transparent materials via transient electronic excitation, Science Advances11, eadv4436 (2025)

    work page 2025

  36. [36]

    Schultze, K

    M. Schultze, K. Ramasesha, C. Pemmaraju, S. Sato, D. Whitmore, A. Gandman, J. S. Prell, L. J. Borja, D. Prendergast, K. Yabana, D. M. Neumark, and S. R. Leone, Attosecond band-gap dynamics in silicon, Science 346, 1348 (2014)

    work page 2014

  37. [37]

    Wachter, C

    G. Wachter, C. Lemell, J. Burgd¨ orfer, S. A. Sato, X.-M. Tong, and K. Yabana, Ab initio simulation of electrical currents induced by ultrafast laser excitation of dielectric materials, Phys. Rev. Lett.113, 087401 (2014)

    work page 2014

  38. [38]

    Sugioka, S

    K. Sugioka, S. Wada, A. Tsunemi, T. Sakai, H. Takai, H. Moriwaki, A. Nakamura, H. T. Hideo Tashiro, and K. T. Koichi Toyoda, Micropatterning of quartz sub- strates by multi-wavelengthvacuum-ultraviolet laser ab- lation, Japanese Journal of Applied Physics32, 6185 (1993)

    work page 1993

  39. [39]

    Zhang, K

    J. Zhang, K. Sugioka, S. Wada, H. Tashiro, and K. Toy- oda, Dual-beam ablation of fused quartz using 266 nm and vuv lasers with different delay-times, Applied Physics A64, 477 (1997)

    work page 1997

  40. [40]

    Zhang, K

    J. Zhang, K. Sugioka, T. Takahashi, K. Toyoda, and K. Midorikawa, Dual-beam ablation of fused silica by multiwavelength excitation process using krf excimer and f2 lasers, Applied Physics A71, 23 (2000)

    work page 2000

  41. [41]

    Obata, K

    K. Obata, K. Sugioka, T. Akane, N. Aoki, K. Toyoda, and K. Midorikawa, Influence of laser fluence and irradiation timing of f2 laser on ablation properties of fused silica in f2-krf excimer laser multi-wavelength excitation process, Applied Physics A73, 755 (2001)

    work page 2001

  42. [42]

    Zoppel, R

    S. Zoppel, R. Merz, J. Zehetner, and G. Reider, Enhance- ment of laser ablation yield by two color excitation, Ap- plied Physics A81, 847 (2005)

    work page 2005

  43. [43]

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

    work page 2013

  44. [44]

    Yang, C.-H

    C.-S. Yang, C.-H. Lin, A. Zaytsev, K.-C. Teng, T.-H. Her, and C.-L. Pan, Femtosecond laser ablation of poly- methylmethacrylate via dual-color synthesized waveform, Applied Physics Letters106, 051902 (2015)

    work page 2015

  45. [45]

    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. B106, 195141 (2022)

    work page 2022

  46. [46]

    Runge and E

    E. Runge and E. K. Gross, Density-functional theory for time-dependent systems, Physical Review Letters52, 997 (1984)

    work page 1984

  47. [47]

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

    work page 2021

  48. [48]

    Koz´ ak, M

    M. Koz´ ak, M. Mart´ ınek, T. Otobe, F. Troj´ anek, and P. Mal´ y, Observation of ultrafast impact ionization in di- amond driven by mid-infrared femtosecond pulses, Jour- nal of Applied Physics128, 015701 (2020)

    work page 2020

  49. [49]

    Miyamoto, Direct treatment of interaction between laser-field and electrons for simulating laser processing of metals, Scientific Reports11, 14626 (2021)

    Y. Miyamoto, Direct treatment of interaction between laser-field and electrons for simulating laser processing of metals, Scientific Reports11, 14626 (2021)

    work page 2021

  50. [50]

    Otobe, High-harmonic generation inα-quartz by electron-hole recombination, Phys

    T. Otobe, High-harmonic generation inα-quartz by electron-hole recombination, Phys. Rev. B94, 235152 (2016)

    work page 2016

  51. [51]

    Tancogne-Dejean, O

    N. Tancogne-Dejean, O. D. M¨ ucke, F. X. K¨ artner, and A. Rubio, Impact of the electronic band structure in high- harmonic generation spectra of solids, Phys. Rev. Lett. 118, 087403 (2017)

    work page 2017

  52. [52]

    Tancogne-Dejean, M

    N. Tancogne-Dejean, M. A. Sentef, and A. Rubio, Ul- trafast modification of hubbarduin a strongly corre- lated material: Ab initio high-harmonic generation in nio, Phys. Rev. Lett.121, 097402 (2018)

    work page 2018

  53. [53]

    Atomic-like high-harmonic generation from two-dimensional materials

    N. Tancogne-Dejean and A. Rubio, Atomic-like high-harmonic generation from two-dimensional materials, Science Advances4, eaao5207 (2018), https://www.science.org/doi/pdf/10.1126/sciadv.aao5207

    work page doi:10.1126/sciadv.aao5207 2018

  54. [54]

    Floss, C

    I. Floss, C. Lemell, K. Yabana, and J. Burgd¨ orfer, In- corporating decoherence into solid-state time-dependent density functional theory, Phys. Rev. B99, 224301 (2019)

    work page 2019

  55. [55]

    Yamada and K

    S. Yamada and K. Yabana, Determining the optimum thickness for high harmonic generation from nanoscale thin films: An ab initio computational study, Phys. Rev. B103, 155426 (2021)

    work page 2021

  56. [56]

    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 field, Phys. Rev. B77, 165104 (2008)

    work page 2008

  57. [57]

    Chu and G

    X. Chu and G. C. Groenenboom, Time-dependent density-functional-theory calculation of high-order- harmonic generation of h 2, Phys. Rev. A85, 053402 (2012)

    work page 2012

  58. [58]

    S. A. Sato, Y. Shinohara, T. Otobe, and K. Yabana, Di- electric response of laser-excited silicon at finite electron temperature, Phys. Rev. B90, 174303 (2014)

    work page 2014

  59. [59]

    S. A. Sato, K. Yabana, Y. Shinohara, T. Otobe, and G. F. Bertsch, Numerical pump-probe experiments of laser-excited silicon in nonequilibrium phase, Phys. Rev. B89, 064304 (2014)

    work page 2014

  60. [60]

    Tancogne-Dejean, O

    N. Tancogne-Dejean, O. D. M¨ ucke, F. X. K¨ artner, and A. Rubio, Ellipticity dependence of high-harmonic gen- eration in solids originating from coupled intraband and interband dynamics, Nature Communications8, 745 (2017)

    work page 2017

  61. [61]

    C. Lian, M. Guan, S. Hu, J. Zhang, and S. Meng, Pho- toexcitation in solids: First-principles quantum simula- tions by real-time tddft, Advanced Theory and Simula- tions1, 1800055 (2018)

    work page 2018

  62. [62]

    C. Yu, S. Jiang, and R. Lu, High order harmonic generation in solids: a review on recent numerical methods, Advances in Physics: X4, 1562982 (2019), https://doi.org/10.1080/23746149.2018.1562982

    work page doi:10.1080/23746149.2018.1562982 2019

  63. [63]

    Yamada and K

    A. Yamada and K. Yabana, Energy transfer from intense laser pulse to dielectrics in time-dependent density func- tional theory, The European Physical Journal D73, 87 (2019)

    work page 2019

  64. [64]

    S. A. Sato, H. H¨ ubener, U. De Giovannini, and A. Rubio, Technical review: Time-dependent density functional theory for attosecond physics ranging from gas-phase to solids, npj Computational Materials11, 233 (2025)

    work page 2025

  65. [65]

    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 nanoscience, Computer Physics Communications235, 356 (2019)

    work page 2019

  66. [66]

    Fuchs and M

    M. Fuchs and M. Scheffler, Ab initio pseudopotentials for electronic structure calculations of poly-atomic sys- tems using density-functional theory, Computer Physics Communications119, 67 (1999)

    work page 1999

  67. [67]

    Tran and P

    F. Tran and P. Blaha, Accurate band gaps of semiconduc- tors and insulators with a semilocal exchange-correlation potential, Physical Review Letters102, 226401 (2009)

    work page 2009

  68. [68]

    S. A. Sato, Y. Taniguchi, Y. Shinohara, and K. Yabana, Nonlinear electronic excitations in crystalline solids us- ing meta-generalized gradient approximation and hybrid functional in time-dependent density functional theory, Journal of Chemical Physics143, 224116 (2015)

    work page 2015

  69. [69]

    Massa, G

    E. Massa, G. Mana, U. Kuetgens, and L. Ferroglio, Mea- surement of the lattice parameter of a silicon crystal, New Journal of Physics11, 053013 (2009)

    work page 2009

  70. [70]

    Kotani and M

    T. Kotani and M. van Schilfgaarde, Impact ionization rates for si, gaas, inas, zns, and gan in thegwapproxi- mation, Phys. Rev. B81, 125201 (2010)

    work page 2010

  71. [71]

    Schinke, P

    C. Schinke, P. C. Peest, J. Schmidt, R. Brendel, K. Bothe, M. R. Vogt, I. Kr¨ oger, S. Winter, A. Schir- macher, S. Lim, H. T. Nguyen, and D. MacDonald, Un- certainty analysis for the coefficient of band-to-band ab- sorption of crystalline silicon, AIP Advances5, 067168 (2015)

    work page 2015

  72. [72]

    D. E. Aspnes and A. A. Studna, Dielectric functions and optical parameters of si, ge, GaP, GaAs, GaSb, InP, InAs, and InSb from 1.5 to 6.0 ev, Physical Review B 27, 985 (1983)

    work page 1983