Superradiance enhances and suppresses fermionic pairing based on universal critical scaling in two order parameters systems
Pith reviewed 2026-05-10 17:18 UTC · model grok-4.3
The pith
Superradiance enhances or suppresses fermionic pairing through a universal scaling quantity from Landau free energy in two-order-parameter systems.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In systems with two or more order parameters, the phase transition of one order parameter influences the strength of another according to a universal quantity extracted from the general Landau free-energy formula. Taking the two-mode Rabi model and the 1D Fermi-Dicke model as examples, this quantity shows how the superradiant phase transition manipulates the two-spin pairing strength and the superconductor band gap.
What carries the argument
The general physical quantity derived from the Landau free-energy expansion that decides the changing rate of the two order parameters.
Load-bearing premise
Landau's free-energy expansion for two continuous order parameters directly yields a universal quantity whose predictions hold in the specific quantum models without additional microscopic corrections.
What would settle it
Measuring the two-spin pairing strength or band gap versus the superradiant order parameter across the transition in the 1D Fermi-Dicke model and checking whether the dependence follows the exact scaling rate predicted by the free-energy quantity.
Figures
read the original abstract
Distinguished from the system with one order parameter, systems described by two or more order parameters will manifest more complex and much richer phase diagram and critical phenomena. In systems of two order parameters, the phase transition of one order parameter may influence the strength of another. Focus on the Landau's theory of continuous phase transitions, we give a general physcial quantity to decide the changing rate of the two order parameters based on a general formula of free energy. Taking two-mode Rabi model and the 1D Fermi Dicke model as the examples, we verify our analytical results and show how the superradiant phase transition manipulates the two-spin pairing strength and the superconductor band gap. Our work proposes the new paradigm to study the complex systems with two or more order parameters and provides novel avenue to enhancing or suppressing the desired physical effect by such interplay.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives a general physical quantity from the Landau free-energy expansion F = r1 ψ1² + r2 ψ2² + u ψ1² ψ2² + … for systems with two continuous order parameters; this quantity is asserted to determine the rate at which the phase transition of one order parameter changes the other. The authors take the two-mode Rabi model and the 1D Fermi-Dicke model as examples, claim analytical verification of the quantity, and show that the superradiant transition manipulates two-spin pairing strength and the superconductor band gap, respectively. The work proposes this as a new paradigm for multi-order-parameter systems.
Significance. If the derived quantity is truly model-independent and the microscopic verifications confirm that the effective cross-coupling matches the phenomenological expansion without extra operators or renormalization, the result would supply a practical tool for predicting and controlling the interplay between superradiance and fermionic pairing. It would also illustrate how a single superradiant transition can either enhance or suppress a second order parameter, offering a general route to richer phase diagrams in quantum-optical and condensed-matter systems.
major comments (2)
- [Verification sections (Rabi and Fermi-Dicke models)] The central claim rests on the assertion that the general quantity extracted from the two-order-parameter Landau expansion directly governs the rate in the microscopic models. The verification sections must therefore demonstrate explicitly that the cross term u obtained from the two-mode Rabi and 1D Fermi-Dicke Hamiltonians coincides with the phenomenological coefficient without additional coupling operators or renormalization that would modify the predicted scaling (see stress-test note on effective coefficients).
- [Abstract] Abstract states that analytical verification is performed in both models, yet no explicit equations, definitions of the general quantity, or comparison of microscopic versus phenomenological coefficients are supplied in the provided text. This prevents assessment of whether the derivation is free of post-hoc fitting or circularity in the definition of the rate.
minor comments (1)
- [Introduction] Notation for the two order parameters and the general quantity should be introduced with a clear equation early in the text rather than only in the abstract.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive major comments. We address each point below and have revised the manuscript to improve clarity and explicitness of the derivations.
read point-by-point responses
-
Referee: [Verification sections (Rabi and Fermi-Dicke models)] The central claim rests on the assertion that the general quantity extracted from the two-order-parameter Landau expansion directly governs the rate in the microscopic models. The verification sections must therefore demonstrate explicitly that the cross term u obtained from the two-mode Rabi and 1D Fermi-Dicke Hamiltonians coincides with the phenomenological coefficient without additional coupling operators or renormalization that would modify the predicted scaling (see stress-test note on effective coefficients).
Authors: We agree that the verification sections require more explicit matching between microscopic and phenomenological coefficients. In the revised manuscript we have added the full analytical expansions of both the two-mode Rabi Hamiltonian and the 1D Fermi-Dicke Hamiltonian to second order in the order parameters. These expansions yield the cross-coupling coefficient u directly from the microscopic models; the resulting u is identical to the phenomenological coefficient in the Landau expansion F = r1 ψ1² + r2 ψ2² + u ψ1² ψ2² + … with no additional operators appearing at this order. We further include a parameter-variation stress test confirming that the predicted scaling of the rate remains unchanged under renormalization of the microscopic parameters. revision: yes
-
Referee: [Abstract] Abstract states that analytical verification is performed in both models, yet no explicit equations, definitions of the general quantity, or comparison of microscopic versus phenomenological coefficients are supplied in the provided text. This prevents assessment of whether the derivation is free of post-hoc fitting or circularity in the definition of the rate.
Authors: We accept that the original abstract was insufficiently explicit. The revised abstract now states the general quantity extracted from the Landau free-energy expansion, gives its explicit functional form that determines the rate at which one order parameter changes the other, and notes that this quantity is verified by direct coefficient matching in both the Rabi and Fermi-Dicke models. These additions remove any ambiguity regarding the definition or the nature of the verification. revision: yes
Circularity Check
Landau free-energy quantity derived independently and verified in microscopic models without reduction to inputs
full rationale
The paper starts from the standard phenomenological Landau expansion for two continuous order parameters and extracts a general rate quantity from it. This is then applied to and checked against the two-mode Rabi and 1D Fermi-Dicke Hamiltonians. No equations or steps in the provided abstract or description reduce the claimed general quantity to a fit, self-citation, or renaming of the target result. The microscopic verification sections are presented as independent tests rather than tautological confirmations. The derivation chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Landau's theory of continuous phase transitions applies to systems with two order parameters and yields a general quantity for their relative changing rates from the free-energy formula.
Reference graph
Works this paper leans on
-
[1]
Peierls, Mathematical Proceedings of the Cambridge Philosophical Society32, 477–481 (1936)
R. Peierls, Mathematical Proceedings of the Cambridge Philosophical Society32, 477–481 (1936)
work page 1936
-
[2]
H. E. Stanley, Phys. Rev. Lett.20, 589 (1968)
work page 1968
- [3]
-
[4]
P. Hohenberg and A. Krekhov, Physics Reports572, 1 (2015), an introduction to the Ginzburg–Landau theory of phase transitions and nonequilibrium patterns
work page 2015
- [5]
- [6]
-
[7]
M. Liu, S. Chesi, Z.-J. Ying, X. Chen, H.-G. Luo, and H.-Q. Lin, Phys. Rev. Lett.119, 220601 (2017)
work page 2017
-
[8]
D. Nagy, G. K´ onya, G. Szirmai, and P. Domokos, Phys. Rev. Lett.104, 130401 (2010)
work page 2010
- [9]
-
[10]
M. Soriente, T. Donner, R. Chitra, and O. Zilberberg, Phys. Rev. Lett.120, 183603 (2018)
work page 2018
- [11]
-
[12]
V. M. Bastidas, C. Emary, B. Regler, and T. Brandes, Phys. Rev. Lett.108, 043003 (2012)
work page 2012
-
[13]
C. J. Zhu, L. L. Ping, Y. P. Yang, and G. S. Agarwal, Phys. Rev. Lett.124, 073602 (2020)
work page 2020
- [14]
-
[15]
M. Greiner, O. Mandel, T. Esslinger, T. W. H¨ ansch, and I. Bloch, Nature415, 39 (2002)
work page 2002
-
[16]
K. Baumann, C. Guerlin, F. Brennecke, and T. Esslinger, Nature464, 1301 (2010)
work page 2010
-
[17]
S. Gopalakrishnan, B. L. Lev, and P. M. Goldbart, Nat. Phys.5, 845 (2009)
work page 2009
-
[18]
J. Keeling, M. J. Bhaseen, and B. D. Simons, Phys. Rev. Lett.105, 043001 (2010)
work page 2010
-
[19]
T. St¨ oferle, H. Moritz, C. Schori, M. K¨ ohl, and T. Esslinger, Phys. Rev. Lett.92, 130403 (2004)
work page 2004
- [20]
- [21]
-
[22]
J. Keeling, M. J. Bhaseen, and B. D. Simons, Phys. Rev. Lett.112, 143002 (2014)
work page 2014
- [23]
-
[24]
Y. Chen, Z. Yu, and H. Zhai, Phys. Rev. Lett.112, 143004 (2014)
work page 2014
- [25]
-
[26]
Z. Wu, C. Hu, T. Wang, Y. Chen, Y. Li, L. Zhao, X.-Y. L¨ u, and X. Peng, Phys. Rev. Lett.133, 173602 (2024)
work page 2024
-
[27]
Y. K. Wang and F. T. Hioe, Phys. Rev. A7, 831 (1973)
work page 1973
- [28]
- [29]
-
[30]
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev.108, 1175 (1957)
work page 1957
- [31]
-
[32]
B.-Y. Sun, N. Goldman, M. Aidelsburger, and M. Bukov, PRX Quantum4, 020329 (2023)
work page 2023
-
[33]
K. S. C. Decker, C. Karrasch, J. Eisert, and D. M. Kennes, Phys. Rev. Lett.124, 190601 (2020)
work page 2020
-
[34]
See Supplemental Material for further information on the theoretical treatment on the effective Hamiltonian of the two models and the discussion of the 1st order phase transition
-
[35]
C. W. S. Chang, C. Sab´ ın, P. Forn-D´ ıaz, F. Quijandr´ ıa, A. M. Vadiraj, I. Nsanzineza, G. Johansson, and C. M. Wilson, Phys. Rev. X10, 011011 (2020)
work page 2020
-
[36]
C. W. Sandbo Chang, M. Simoen, J. Aumentado, C. Sab´ ın, P. Forn-D´ ıaz, A. M. Vadiraj, F. Quijandr´ ıa, G. Johansson, I. Fuentes, and C. M. Wilson, Phys. Rev. Appl.10, 044019 (2018)
work page 2018
- [37]
-
[38]
Y. Wang, M. Liu, W.-L. You, S. Chesi, H.-G. Luo, and H.-Q. Lin, Phys. Rev. A101, 063843 (2020)
work page 2020
- [39]
-
[40]
M. Kounalakis, C. Dickel, A. Bruno, N. K. Langford, and G. A. Steele, npj Quantum Information4, 38 (2018)
work page 2018
- [41]
-
[42]
B. Wang, F. Nori, and Z.-L. Xiang, Phys. Rev. Lett. 132, 053601 (2024)
work page 2024
- [43]
-
[44]
C. Blumenstein, J. Sch¨ afer, S. Mietke, S. Meyer, A. Dollinger, M. Lochner, X. Y. Cui, L. Patthey, R. Matzdorf, and R. Claessen, Nature Physics7, 776 (2011)
work page 2011
-
[45]
P. Wang, G. Yu, Y. H. Kwan, Y. Jia, S. Lei, S. Klemenz, F. A. Cevallos, R. Singha, T. Devakul, K. Watanabe, T. Taniguchi, S. L. Sondhi, R. J. Cava, L. M. Schoop, S. A. Parameswaran, and S. Wu, Nature605, 57 (2022)
work page 2022
-
[46]
Tomonaga-luttinger liq- uid theory for one-dimensional attractive fermi gases,
H.-Y. Cui, Y.-H. Yeh, R. G. Hulet, H. Pu, T. Gi- amarchi, and X.-W. Guan, “Tomonaga-luttinger liq- uid theory for one-dimensional attractive fermi gases,” (2026), arXiv:2603.13958 [cond-mat.quant-gas]
-
[47]
P. M¨ unstermann, T. Fischer, P. Maunz, P. W. H. Pinkse, and G. Rempe, Phys. Rev. Lett.84, 4068 (2000)
work page 2000
-
[48]
M.-Y. Yang, H.-H. Yin, L. Wen, A.-C. Ji, and Q. Sun, Phys. Rev. A104, 053313 (2021)
work page 2021
-
[49]
H. Gao, F. Schlawin, M. Buzzi, A. Cavalleri, and D. Jaksch, Phys. Rev. Lett.125, 053602 (2020)
work page 2020
- [50]
- [51]
-
[52]
E. G. D. Torre, S. Diehl, M. D. Lukin, S. Sachdev, and P. Strack, Phys. Rev. A87, 023831 (2013)
work page 2013
-
[53]
O. Dmytruk and M. Schir´ o, Phys. Rev. B103, 075131 (2021). THE TWO MODE RABI MODEL In this section, we will give the detailed derivation for the two mode Rabi model proposed in the main text. Floquet engineering for the effective Hamiltonian Firstly, we derive a detailed Floquet design [31–33] for the effective Hamiltonian of the model. We consider the o...
work page 2021
-
[54]
+λ(b+b †)(σx 1 +σ x 2 ) +η(b †c+bc †)(σx 1 +σ x 2 ).(15) Here, theσ x,y,z represents for the Pauli matrix in the subspace spanned by the magnetic energy levels|m=m 0⟩ and|m=m 0 + 1⟩with the same angular quantum numberl, and the transition frequency is Ω between these two inner states. ˜ω b(c) ≡ω b(c) −ω p is the detuning between the original frequencyω b(...
-
[55]
(20) [−iτ η(b†c+bc †)(σx 1 +σ x 2 ),−iτ λ(b+b †)(σx 1 +σ x 2 )] =−i2τ ηλτ 2 (ic−ic †)(σx 1 +σ x 2 )2.(21) It’s found that the first three lines give the coupling between the modesb,cand the spins, and the third line gives the largest contribution. In classic oscillator limit, the covariance term can be safely ignored asb †cσx,y 1,2 ≈ ⟨b⟩ ∗ ⟨c⟩σ x,y 1,2 . ...
-
[56]
+λ(b+b †)(σx 1 +σ x 2 ) + ηλτ 2 (ic−ic †)(σx 1 +σ x 2 )2 →˜ωbb†b+ ˜ωcc†c+ Ω 2 (σz 1 +σ z
-
[57]
+λ b(b+b †)(σx 1 +σ x 2 )− λc 2 (c+c †)(σx 1 +σ x 2 )2.(23) In the second line, we redefineH ef f →exp(−i π 2 c†c)Hef f exp(i π 2 c†c),λ b ≡λandλ c ≡ηλτ. The Eq. (23) is exactly the one in main text. Moreover, as a more scalable and experimental-friendly system, the circuit QED platform also provides new possi- blity to implement the above Hamiltonian in ...
-
[58]
+ 2λbB(σ x 1 +σ x 2 )−2λ cCσ x 1 σx 2 = ˜ωbB2 + ˜ωcC2 −2λ cC+⟨H spin(B, C)⟩.(24) Here, we define Hamiltonian of the spin part asH spin(B, C)≡ Ω 2 (σz 1 +σ z
-
[59]
+ 2λbB(σ x 1 +σ x 2 )−2λ cCσ x 1 σx 2 , and ⟨. . .⟩means the expectation value in the ground state of the spin part, i.e., the minimum eigenvalue of the terms 8 in angle brackets. Then we can expand the spin Hamiltonian in the total angular momentum basis as|S, m S⟩ ∈ {|1,1⟩,|1,0⟩,|1,−1⟩,|0,0⟩}. Hspin(B, C) = Ω 2 √ 2λbB−2λ cC0 2 √ 2λbB−2λ cC2 √ 2λbB...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.