Photoinduced electron-electron pairing in the extended Falicov-Kimball model
Pith reviewed 2026-05-25 10:18 UTC · model grok-4.3
The pith
Pulse irradiation induces interband electron-electron pair correlations in the extended Falicov-Kimball model while suppressing initial excitonic correlations via an internal SU(2) structure.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Pulse irradiation induces interband electron-electron pair correlation in the photoexcited states of the EFKM while strongly suppressing the excitonic electron-hole pair correlation of the initial ground state; the photoexcited states contain eigenstates with a finite number of interband electron-electron pairs, and this occurs because of the internal SU(2) pairing structure that is preserved even for nonbipartite lattices with direct-type band structure.
What carries the argument
The internal SU(2) pairing structure of the extended Falicov-Kimball model, which maps photoexcited states onto eigenstates carrying interband electron-electron pairs.
Load-bearing premise
The internal SU(2) pairing structure remains intact and governs the dynamics for the pulses and lattices considered.
What would settle it
Measure the time-dependent interband electron-electron pair correlation function after the pulse; if it fails to rise while the electron-hole correlation falls, or if the post-pulse state does not overlap with eigenstates containing a finite number of electron-electron pairs, the claimed induction mechanism does not hold.
Figures
read the original abstract
By employing the time-dependent exact diagonalization method, we investigate the photoexcited states of the excitonic insulator in the extended Falicov-Kimball model (EFKM). We here show that the pulse irradiation can induce the interband electron-electron pair correlation in the photoexcited states, while the excitonic electron-hole pair correlation in the initial ground state is strongly suppressed. We also show that the photoexcited states contains the eigenstates of the EFKM with a finite number of interband electron-electron pairs, which are responsible for the enhancement of the electron-electron pair correlation. The mechanism found here is due to the presence of the internal SU(2) pairing structure in the EFKM and thus it is essentially the same as that for the photoinduced $\eta$-pairing in the repulsive Hubbard model reported recently [T. Kaneko et al., Phys. Rev. Lett. ${\bf 122}$, 077002 (2019)]. This also explains why the nonlinear optical response is effective to induce the electron-electron pairs in the photoexcited states of the EFKM. Furthermore, we show that, unlike the $\eta$-pairing in the Hubbard model, the internal SU(2) structure is preserved even for a nonbipartite lattice when the EFKM has the direct-type band structure, in which the pulse irradiation can induce the electron-electron pair correlation with momentum ${\it {\bf q}}$ = ${\textbf 0}$ in the photoexcited states. We also discuss briefly the effect of a perturbation that breaks the internal SU(2) structure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript employs time-dependent exact diagonalization to study the extended Falicov-Kimball model in its excitonic insulator regime. It claims that optical pulse irradiation induces interband electron-electron pair correlations in the photoexcited states while strongly suppressing the excitonic electron-hole correlations present in the ground state. The photoexcited states are shown to contain eigenstates with a finite number of interband electron-electron pairs. The mechanism is attributed to the model's internal SU(2) pairing symmetry, which remains intact for direct-type bands even on nonbipartite lattices (unlike the Hubbard model), enabling q=0 pairing; a brief qualitative discussion of SU(2)-breaking perturbations is included.
Significance. If the numerical results hold, the work establishes a symmetry-protected mechanism for photoinduced electron-electron pairing in an excitonic insulator, extending the η-pairing analogy from the repulsive Hubbard model to a broader class of lattices and band structures. The exact-diagonalization approach supplies direct, non-perturbative evidence for the role of the internal SU(2) structure in nonlinear optical responses, offering a concrete, falsifiable prediction for light-induced pairing that could be tested in related models or materials.
major comments (1)
- [Discussion of perturbations] Final discussion section (on perturbations breaking the internal SU(2) structure): the robustness of the reported photoinduced electron-electron pair correlations is asserted to fail under SU(2)-breaking terms, yet the manuscript supplies only a qualitative statement without time-dependent ED data, threshold values for perturbation strength (e.g., next-nearest-neighbor hopping or interaction asymmetry), or any quantitative estimate of how small such terms must remain for the q=0 enhancement to survive. Because all quantitative results are obtained strictly inside the symmetric model, this omission leaves the central claim's applicability to realistic extensions unquantified and load-bearing for the mechanism's generality.
minor comments (2)
- [Abstract] The abstract and introduction refer to 'the pulse irradiation' without summarizing the specific pulse parameters (amplitude, frequency, duration) used in the time-dependent simulations; these should be stated explicitly for reproducibility.
- [Results section] Notation for the interband electron-electron pair correlation function is introduced without an explicit equation reference in the main text; adding a numbered equation would improve clarity when comparing ground-state and photoexcited values.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our work and the constructive comment. We address the major point below.
read point-by-point responses
-
Referee: Final discussion section (on perturbations breaking the internal SU(2) structure): the robustness of the reported photoinduced electron-electron pair correlations is asserted to fail under SU(2)-breaking terms, yet the manuscript supplies only a qualitative statement without time-dependent ED data, threshold values for perturbation strength (e.g., next-nearest-neighbor hopping or interaction asymmetry), or any quantitative estimate of how small such terms must remain for the q=0 enhancement to survive. Because all quantitative results are obtained strictly inside the symmetric model, this omission leaves the central claim's applicability to realistic extensions unquantified and load-bearing for the mechanism's generality.
Authors: We agree that the discussion of SU(2)-breaking perturbations is only qualitative, as explicitly noted in the manuscript. The central mechanism relies on the exact internal SU(2) symmetry, which protects the q=0 pairing channel; any breaking term (e.g., next-nearest-neighbor hopping or interaction asymmetry) will lift the degeneracy and suppress the enhancement proportionally to its strength relative to the excitonic gap. While additional time-dependent ED data on perturbed models would be desirable, such calculations are computationally intensive because symmetry breaking enlarges the effective Hilbert space and removes conserved quantities used in the current exact-diagonalization implementation. We will revise the final section to include order-of-magnitude estimates based on the existing ED energy scales, indicating that the perturbation must remain smaller than the photoinduced pairing correlation energy for the effect to survive. revision: partial
Circularity Check
Minor self-citation to Hubbard-model mechanism; central EFKM results from independent TD-ED
full rationale
The paper's core results derive from direct time-dependent exact diagonalization of the EFKM Hamiltonian under pulse driving, showing suppression of excitonic correlations and enhancement of interband electron-electron pair correlations. The cited Kaneko et al. PRL (overlapping authors) is invoked only to identify the shared SU(2) origin of the mechanism, without substituting for or reducing the present numerical evidence. The SU(2) structure is stated as a property of the EFKM itself, preserved under the studied conditions, and the perturbation discussion is brief and qualitative. No self-definitional equations, fitted inputs renamed as predictions, or load-bearing uniqueness theorems appear. This is a standard low-level self-citation that does not force the reported findings.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The extended Falicov-Kimball model possesses an internal SU(2) pairing symmetry that is preserved under the chosen pulse protocol and lattice geometry.
- domain assumption Time-dependent exact diagonalization on finite clusters accurately represents the photoexcited dynamics of the infinite system.
Reference graph
Works this paper leans on
-
[1]
may be written in momentum ( k) space as ˆH0 = ∑ k,α ǫα (k)ˆc† k,α ˆck,α (2) with ǫ1(k) = − 2t(1) h ∑ τ coskτ − D 2 (3) and ǫ2(k) = − 2t(2) h ∑ τ coskτ + D 2, (4) where kτ = k · aτ and aτ is the vector between the nearest-neighbor sites i and j. Here, we implicitly as- sume that the hoppings are finite between sites connected through the primitive translat...
-
[2]
essentially equivalent to the Hubbard model
and (b) an excitonic insulator in the EFKM with t(1) h ·t(2) h < 0. essentially equivalent to the Hubbard model. Therefore, as in the case of the Hubbard model, the EFKM with t(1) h =t(2) h has the internal SU(2) structure defined by the η-pairing operators [ 43–45]. Below, we will show that the EFKM with t(1) h = −t(2) h displays the different internal SU(...
-
[3]
with the ∆-pairing operators that themselves satisfy the SU(2) commutation relations in Eq. ( 12) [46]. Eqs. ( 14) and ( 17) imply that any eigenstate of ˆH is also the eigenstate of ˆ∆2 and ˆ∆z with eigenvalues ∆(∆ +
-
[4]
We denote this eigenstate as |∆, ∆z⟩
and ∆ z, respectively [ 47]. We denote this eigenstate as |∆, ∆z⟩. Assuming that N1 ≥ N2 and L − N1 +N2 is even, |∆, ∆z⟩ can take ∆ = 0 , 1, 2, · · ·, (L − N1 +N2)/ 2 and ∆z = − ∆, − ∆+1, · · ·, ∆− 1, ∆. Note that ∆ z = 0 at half filling with N1 +N2 =L. The state |∆, ∆z = − ∆⟩ is the lowest weight state (L WS) that satisfies ˆ∆− |∆, ∆z = − ∆⟩ = 0 [ 44, 45]....
-
[5]
Repulsive Hubbard model The EFKM with t(2) h = −t(1) h =th can be transformed into the repulsive Hubbard model by the following gauge transformation: ˆcj, 1 → (− 1)j ˆdj, ↑ ˆcj, 2 → ˆdj, ↓ (20) Indeed, the EFKM ˆH is transformed as ˆH → ˆHR = − th ∑ ⟨i,j ⟩,σ ( ˆd† i,σ ˆdj,σ + H. c. ) − D 2 ∑ j (ˆnj, ↑ − ˆnj, ↓ ) +U ∑ j ˆnj, ↑ ˆnj, ↓ (21) 4 provided that t...
-
[6]
can transform the hopping term in the presence of the Peierls phase in Eq. (
-
[7]
Note that here the hoppings are assumed to be finite only between sites on different sub- lattices
as ∑ α t(α ) h e− iA(t)·(Ri− Rj )ˆc† i,α ˆcj,α → th ∑ σ e− iA(t)·(Ri− Rj ) ˆd† i,σ ˆdj,σ , (24) which is exactly the hopping term with the Peierls phase in the Hubbard model. Note that here the hoppings are assumed to be finite only between sites on different sub- lattices. Therefore, even the photoinduced dynamics of the EFKM with t(1) h = −t(2) h is equiv...
-
[8]
Attractive Hubbard model It is also instructive to consider the correspondence between the EFKM and the attractive Hubbard model. Since the repulsive Hubbard model and the attractive Hubbard model are mutually transformed via the so- called Shiba transformation [ 45, 53], it is obvious that the EFKM with t(2) h = −t(1) h = th can also be trans- formed int...
-
[9]
transforms the hopping term with the Peierls phase in Eq. ( 18) as ∑ α t(α ) h e− iA(t)·(Ri− Rj )ˆc† i,α ˆcj,α → the+iA(t)·(Ri− Rj ) ˆd† i, ↑ ˆdj, ↑ +the− iA(t)·(Ri− Rj ) ˆd† i, ↓ ˆdj, ↓, (29) 5 which is different from the hopping term with the Peierls phase in the attractive Hubbard model ˆHA. The differ- ence of the photoexcited dynamics has been discusse...
-
[10]
75th in L = 16, for which the ground state of the 1D EFKM, i.e., the initial state before the pulse irradiation, is the excitonic insulator with N1 = 12 and N2 = 4 (see Fig. 2). Figure 3(a) shows the time evolution of the real-space electron-electron pair correlation function P (j,t ). We confirm the enhancement of P (j,t ) at j = 0, correspond- ing to nd(...
-
[11]
Therefore, the structure factor FIG
and ( 17)]. Therefore, the structure factor FIG. 5. (a) All the eigenenergies εm and P (q = 0) for the eigenstates |ψm⟩ of the half-filled 1D EFKM ˆH with t(2) h = − t(1) h = th for L = 10 under PBC at U = 8 th and D = 0.4th, where N1 = 6 and N2 = 4. The color of each point (diamond) indicates the weight | ⟨ψm|Ψ(t)⟩ |2 of the eigenstate |ψm⟩ in the photoin...
-
[12]
This is ex- pected because, as described in Sec
Square lattice Similarly, ∆ pairs can be photoinduced in the two di- mensional (2D) EFKM in the square lattice. This is ex- pected because, as described in Sec. II D 1, when the sys- tem is bipartite, the 2D EFKM with t(1) h = −t(2) h can be mapped onto the repulsive Hubbard model where η pairs can be induced by the pulse irradiation [ 41]. Since the η pa...
-
[13]
Triangular lattice A nontrivial system is the 2D EFKM in the triangu- lar lattice, for which there is no correspondence to the repulsive Hubbard model, as discussed in Sec. II D 1. In contrast to the case of the η-pairing operators in the Hub- bard model, the ∆-pairing operators in the EFKM sat- isfy [ ˆH, ˆ∆± ] = ±U ˆ∆± , regardless of whether the lattic...
-
[14]
A similar relation (a) (b) FIG
is satisfied even for the EFKM in nonbipartite lattices such as the triangular lattice. A similar relation (a) (b) FIG. 13. (a) Time evolution of the on-site electron-electro n pair correlation function P (j, t) and (b) P (j, t) at t = 0 (blue circles) and t = 30 /th (orange squares). The results are for the 1D EFKM under PBC with t(1) h = t(2) h = th, U =...
- [15]
- [16]
-
[17]
B. I. Halperin and T. M. Rice, Rev. Mod. Phys. 40, 755 (1968)
work page 1968
- [18]
-
[19]
J. Nasu, T. Watanabe, M. Naka, and S. Ishihara, Phys. Rev. B 93, 205136 (2016)
work page 2016
- [20]
-
[21]
H. Cercellier, C. Monney, F. Clerc, C. Battaglia, L. Despont, M. G. Garnier, H. Beck, P. Aebi, L. Patthey, H. Berger, and L. Forr´ o, Phys. Rev. Lett. 99, 146403 (2007)
work page 2007
- [22]
- [23]
- [24]
- [25]
- [26]
-
[27]
Y. Wakisaka, T. Sudayama, K. Takubo, T. Mi- zokawa, M. Arita, H. Namatame, M. Taniguchi, N. Katayama, M. Nohara, and H. Takagi, Phys. Rev. Lett. 103, 026402 (2009)
work page 2009
- [28]
-
[29]
K. Seki, Y. Wakisaka, T. Kaneko, T. Toriyama, T. Konishi, T. Sudayama, N. L. Saini, M. Arita, H. Namatame, M. Taniguchi, N. Katayama, M. No- hara, H. Takagi, T. Mizokawa, and Y. Ohta, Phys. Rev. B 90, 155116 (2014)
work page 2014
- [30]
-
[31]
Y. F. Lu, H. Kono, T. I. Larkin, A. W. Rost, T. Takayama, A. V. Boris, B. Keimer, and H. Takagi, Nat. Commun. 8, 14408 (2017)
work page 2017
-
[32]
K. Sugimoto, S. Nishimoto, T. Kaneko, and Y. Ohta, Phys. Rev. Lett. 120, 247602 (2018)
work page 2018
- [33]
-
[34]
E. M¨ ohr-Vorobeva, S. L. Johnson, P. Beaud, U. Staub, R. De Souza, C. Milne, G. Ingold, J. Demsar, H. Schaefer, and A. Titov, Phys. Rev. Lett. 107, 036403 (2011)
work page 2011
-
[35]
S. Hellmann, T. Rohwer, M. Kall¨ ane, K. Hanff, C. Sohrt, A. Stange, A. Carr, M. M. Murnane, H. C. Kapteyn, L. Kipp, M. Bauer, and K. Rossnagel, Nat. Commun. 3, 1069 (2012)
work page 2012
- [36]
-
[37]
S. Mathias, S. Eich, J. Urbancic, S. Michael, A. V. Carr, S. Emmerich, A. Stange, T. Popmintchev, T. Rohwer, M. Wiesenmayer, A. Ruffing, S. Jakobs, S. Hellmann, P. Matyba, C. Chen, L. Kipp, M. Bauer, H. C. Kapteyn, H. C. Schneider, K. Rossnagel, M. M. Murnane, and M. Aeschlimann, Nat. Commun. 7, 12902 (2016)
work page 2016
- [38]
-
[39]
S. Mor, M. Herzog, D. Goleˇ z, P. Werner, M. Eckstein, N. Katayama, M. Nohara, H. Tak- agi, T. Mizokawa, C. Monney, and J. St¨ ahler, Phys. Rev. Lett. 119, 086401 (2017)
work page 2017
-
[40]
S. Mor, M. Herzog, J. Noack, N. Katayama, M. Nohara, H. Takagi, A. Trunschke, T. Mizokawa, C. Monney, and J. St¨ ahler, Phys. Rev. B 97, 115154 (2018)
work page 2018
-
[41]
D. Werdehausen, T. Takayama, M. H¨ oppner, G. Al- brecht, A. W. Rost, Y. Lu, D. Manske, H. Takagi, and S. Kaiser, Sci. Adv. 4, eaap8652 (2018)
work page 2018
-
[42]
D. Werdehausen, T. Takayama, G. Al- brecht, Y. Lu, H. Takagi, and S. Kaiser, J. Phys.: Condens. Matter 30, 305602 (2018)
work page 2018
-
[43]
K. Okazaki, Y. Ogawa, T. Suzuki, T. Yamamoto, T. Someya, S. Michimae, M. Watanabe, Y. Lu, M. No- hara, H. Takagi, N. Katayama, H. Sawa, M. Fujisawa, T. Kanai, N. Ishii, J. Itatani, T. Mizokawa, and S. Shin, Nat. Commun. 9, 4322 (2018)
work page 2018
- [44]
-
[45]
Y. Murakami, D. Goleˇ z, M. Eckstein, and P. Werner, Phys. Rev. Lett. 119, 247601 (2017)
work page 2017
- [46]
- [47]
-
[48]
C. D. Batista, Phys. Rev. Lett. 89, 166403 (2002) ; Phys. Rev. Lett. 90, 199901 (2003)
work page 2002
-
[49]
D. Ihle, M. Pfafferott, E. Burovski, F. X. Bronold, and H. Fehske, Phys. Rev. B 78, 193103 (2008)
work page 2008
-
[50]
K. Seki, R. Eder, and Y. Ohta, Phys. Rev. B 84, 245106 (2011)
work page 2011
- [51]
- [52]
- [53]
- [54]
- [55]
-
[56]
Direct (indirect) implies that the valence band maximu m and the conduction band minimum are located at the same (different) position(s) in the Brillouin zone
-
[57]
C. N. Yang, Phys. Rev. Lett. 63, 2144 (1989)
work page 1989
-
[58]
F. H. L. Essler, V. E. Korepin, and K. Schoutens, Phys. Rev. Lett. 67, 3848 (1991)
work page 1991
-
[59]
F. H. Essler, H. Frahm, F. G¨ ohmann, A. Kl¨ umper, and V. E. Korepin, The One-Dimensional Hubbard Model (Cambridge University Press, Cambridge, 2005)
work page 2005
-
[60]
( 16), in addition to [ ˆH, ˆ∆z] = 0
In order for ˆH to be fully rotation symmetric with respect to the ∆-pairing operators, ˆH has to satisfy [ ˆH, ˆ∆±] = 0, instead of Eq. ( 16), in addition to [ ˆH, ˆ∆z] = 0. If we add the on-site energy term − U 2 ∑ j(ˆnj,1 + ˆnj,2) to the EFKM in Eq. ( 1), the model satisfies these commutation relations and thus preserves the internal SU(2) symme- try wi...
-
[61]
ˆ∆ and ˆ∆z with the hat are operators and ∆ and ∆ z without the hat are c-numbers
- [62]
-
[63]
G. De Filippis, V. Cataudella, E. A. Nowadnick, T. P. Devereaux, A. S. Mishchenko, and N. Nagaosa, Phys. Rev. Lett. 109, 176402 (2012)
work page 2012
-
[64]
H. Lu, S. Sota, H. Matsueda, J. Bonˇ ca, and T. Tohyama, Phys. Rev. Lett. 109, 197401 (2012)
work page 2012
- [65]
-
[66]
Y. Wang, M. Claassen, B. Moritz, and T. P. Devereaux, Phys. Rev. B 96, 235142 (2017)
work page 2017
- [67]
- [68]
-
[69]
T. J. Park and J. Light, J. Chem. Phys. 85, 5870 (1986)
work page 1986
-
[70]
N. Mohankumar and S. M. Auerbach, Comput. Phys. Commun. 175, 473 (2006)
work page 2006
-
[71]
Z. Lenarˇ ciˇ c, D. Goleˇ z, J. Bonˇ ca, and P. Prelovˇ sek, Phys. Rev. B 89, 125123 (2014)
work page 2014
- [72]
-
[73]
R. M. Fye, M. J. Martins, D. J. Scalapino, J. Wagner, and W. Hanke, Phys. Rev. B 44, 6909 (1991)
work page 1991
-
[74]
E. Jeckelmann, F. Gebhard, and F. H. L. Essler, Phys. Rev. Lett. 85, 3910 (2000)
work page 2000
-
[75]
J. J. Sakurai, Modern Quantum Mechanics (Addison- Wesley, Reading, MA, 1994); M. E. Rose, Elementary Theory of Angular Momentum (Wiley, New York, 1967)
work page 1994
- [76]
-
[77]
W. Hu, S. Kaiser, D. Nicoletti, C. R. Hunt, I. Gierz, M. C. Hoffmann, M. Le Tacon, T. Loew, B. Keimer, and A. Cavalleri, Nat. Mater. 13, 705 (2014)
work page 2014
- [78]
-
[79]
M. Mitrano, A. Cantaluppi, D. Nicoletti, S. Kaiser, A. Perucchi, S. Lupi, P. Di Pietro, D. Pontiroli, M. Ricc` o, S. R. Clark, D. Jaksch, and A. Cavalleri, Nature (London) 530, 461 (2016)
work page 2016
-
[80]
A. Cantaluppi, M. Buzzi, G. Jotzu, D. Nicoletti, M. Mi- trano, D. Pontiroli, M. Ricc` o, A. Perucchi, P. Di Pietro, and A. Cavalleri, Nat. Phys. 14, 837 (2018)
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.