Grassmann time-evolving matrix product operators for fermionic impurities coupled to a superconducting bath
Pith reviewed 2026-05-08 07:23 UTC · model grok-4.3
The pith
The Bogoliubov transformation enables direct adaptation of Grassmann time-evolving matrix product operators to fermionic impurities in superconducting baths via the Nambu formalism.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By employing the Bogoliubov transformation for the superconducting bath, one obtains the analytic expression of the Feynman-Vernon influence functional in a similar form to the case of a normal bath, after which the core algorithms of GTEMPO can be straightforwardly adapted to the Nambu formalism.
What carries the argument
The Bogoliubov transformation on the superconducting bath that yields a Feynman-Vernon influence functional with structure preserved for temporal matrix product operator encoding in GTEMPO.
If this is right
- The method matches exact diagonalization results in several solvable cases for imaginary- and real-time evolution.
- It converges with CT-QMC in DMFT iterations for non-integrable interacting impurity models.
- Both equilibrium and non-equilibrium calculations are supported within the same adapted framework.
- GTEMPO becomes usable as an impurity solver in Nambu-formalism DMFT and its dynamical extensions.
Where Pith is reading between the lines
- The preserved structure suggests GTEMPO can now address pairing correlations at the impurity without redesigning the algorithm.
- Real-time results may be used to study quench dynamics in superconducting environments.
- Connections to other symmetry-broken baths could be explored using analogous transformations.
Load-bearing premise
The influence functional after the Bogoliubov transformation retains the exact structural properties required by the GTEMPO algorithm, without introducing any loss of accuracy or need for further approximations.
What would settle it
A mismatch between the adapted GTEMPO results and exact diagonalization in an exactly solvable model with a superconducting bath would falsify the claim that the adaptation works without accuracy loss.
Figures
read the original abstract
The Grassmann time-evolving matrix product operator (GTEMPO) method, which represents the Feynman-Vernon influence functional as a temporal matrix product state, has been shown to be a flexible and potentially scalable solution for fermionic quantum impurity problems. In this work, we extend GTEMPO to solve fermionic impurity problems in the Nambu formalism, in which the impurity is coupled to a superconducting bath. A key insight is that by employing the Bogoliubov transformation for the superconducting bath, one could obtain the analytic expression of the Feynman-Vernon influence functional in a similar form to the case of a normal bath, after which the core algorithms of GTEMPO can be straightforwardly adapted. We demonstrate the accuracy of our method by benchmarking it against exact diagonalization in several exactly solvable cases, and against the continuous-time quantum Monte Carlo method using converged dynamical mean field theory (DMFT) iterations on the imaginary contour in the non-integrable case. In all cases, we perform both imaginary- and real-time calculations to illustrate the flexibility of our method. These results illustrate that our method could be potentially useful as an impurity solver in DMFT as well as its non-equilibrium extension for fermionic impurity problems in the Nambu formalism.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript extends the Grassmann time-evolving matrix product operator (GTEMPO) method to fermionic impurity problems coupled to superconducting baths within the Nambu formalism. By applying the Bogoliubov transformation to the quadratic bath, an analytic expression for the Feynman-Vernon influence functional is obtained that is structurally similar to the normal-metal case. This allows direct adaptation of the GTEMPO algorithms for representing the influence functional as a temporal matrix product state. The method is benchmarked against exact diagonalization for solvable models and against converged CT-QMC DMFT for interacting cases, with calculations performed on both the imaginary and real-time axes.
Significance. If the central claim holds, the work supplies a flexible and potentially scalable impurity solver for DMFT and its non-equilibrium extensions in the presence of superconductivity. The exact analytic influence functional (arising because the bath remains Gaussian after the unitary transformation) together with explicit benchmarks on both time axes against ED and CT-QMC constitute concrete strengths that could enable studies of unconventional superconductors.
Simulated Author's Rebuttal
We thank the referee for their positive and accurate summary of our manuscript, which correctly identifies the extension of GTEMPO to the Nambu formalism for superconducting baths via the Bogoliubov transformation, the resulting analytic influence functional, and the benchmarks against exact diagonalization and CT-QMC on both imaginary and real-time axes. We appreciate the recognition of the method's potential utility as an impurity solver for DMFT and non-equilibrium extensions in superconducting systems. No specific major comments were provided in the report.
Circularity Check
No significant circularity; derivation uses standard external transformation and prior independent method
full rationale
The paper's central derivation applies the established Bogoliubov transformation (a standard technique from superconductivity literature, independent of this work) to the quadratic bath, yielding an exact analytic Feynman-Vernon influence functional whose 2-point correlator structure remains Gaussian and thus directly compatible with the existing GTEMPO MPO representation. This adaptation is presented as a straightforward substitution of Nambu-space correlators rather than a re-derivation or redefinition of GTEMPO itself. The prior GTEMPO framework is cited as an external algorithmic foundation, not as a self-referential load-bearing premise. Independent validation is supplied via benchmarks against exact diagonalization on solvable models and converged CT-QMC DMFT on interacting cases, confirming that the extension preserves accuracy without hidden reductions to fitted inputs or self-definitional loops. No steps in the chain reduce by construction to the paper's own outputs or self-citations.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Bogoliubov transformation diagonalizes the superconducting bath and produces an influence functional of the same form as a normal bath
Reference graph
Works this paper leans on
-
[1]
P. W. Anderson, Localized magnetic states in metals, Phys. Rev.124, 41 (1961)
work page 1961
-
[2]
A. Georges, The beauty of impurities: Two revivals of friedel’s virtual bound-state concept, Comptes Rendus Physique17, 430 (2016), physique de la mati` ere con- dens´ ee au XXIe si` ecle: l’h´ eritage de Jacques Friedel
work page 2016
-
[3]
G. D. Mahan,Many-Particle Physics(Springer; 3nd edi- tion, 2000)
work page 2000
-
[4]
Weiss,Quantum Dissipative Systems(World Scien- tific, Singapore, 1993)
U. Weiss,Quantum Dissipative Systems(World Scien- tific, Singapore, 1993)
work page 1993
-
[5]
W. Metzner and D. Vollhardt, Correlated lattice fermions ind=∞dimensions, Phys. Rev. Lett.62, 324 (1989)
work page 1989
-
[6]
A. Georges and G. Kotliar, Hubbard model in infinite dimensions, Phys. Rev. B45, 6479 (1992)
work page 1992
-
[7]
A. Georges, G. Kotliar, W. Krauth, and M. J. Rozen- berg, Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions, Rev. Mod. Phys.68, 13 (1996)
work page 1996
-
[8]
G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Tran- sition from metallic to tunneling regimes in supercon- ducting microconstrictions: Excess current, charge im- balance, and supercurrent conversion, Phys. Rev. B25, 4515 (1982)
work page 1982
-
[9]
J. C. Cuevas, A. Mart´ ın-Rodero, and A. L. Yeyati, Hamil- tonian approach to the transport properties of supercon- ducting quantum point contacts, Phys. Rev. B54, 7366 (1996)
work page 1996
-
[10]
B. M. Andersen, K. Flensberg, V. Koerting, and J. Paaske, Nonequilibrium transport through a spinful quantum dot with superconducting leads, Phys. Rev. Lett.107, 256802 (2011)
work page 2011
-
[11]
A. Y. Kitaev, Unpaired majorana fermions in quantum wires, Physics-Uspekhi44, 131 (2001)
work page 2001
-
[12]
J. Alicea, New directions in the pursuit of majorana fermions in solid state systems, Reports on Progress in Physics75, 076501 (2012)
work page 2012
-
[13]
A. Sakurai, Comments on superconductors with mag- netic impurities, Progress of Theoretical Physics 44, 1472 (1970), https://academic.oup.com/ptp/article- pdf/44/6/1472/5378530/44-6-1472.pdf
work page 1970
-
[14]
E. M¨ uller-Hartmann and J. Zittartz, Kondo effect in su- perconductors, Phys. Rev. Lett.26, 428 (1971)
work page 1971
-
[15]
T. Matsuura, The effects of impurities on superconduc- tors with kondo effect, Progress of Theoretical Physics 57, 1823 (1977), https://academic.oup.com/ptp/article- pdf/57/6/1823/5461826/57-6-1823.pdf
work page 1977
-
[16]
A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Impurity- induced states in conventional and unconventional su- perconductors, Rev. Mod. Phys.78, 373 (2006)
work page 2006
-
[17]
A. Yazdani, B. A. Jones, C. P. Lutz, M. F. Crommie, and D. M. Eigler, Probing the local effects of magnetic im- purities on superconductivity, Science275, 1767 (1997), https://www.science.org/doi/pdf/10.1126/science.275.5307.1767
-
[18]
K. J. Franke, G. Schulze, and J. I. Pascual, Com- petition of superconducting phenomena and kondo screening at the nanoscale, Science332, 940 (2011), https://www.science.org/doi/pdf/10.1126/science.1202204
-
[19]
L. P. Gor’kov, On the energy spectrum of superconduc- tors, Sov. Phys. JETP7, 505 (1958)
work page 1958
-
[20]
Nambu, Quasi-particles and gauge invariance in the theory of superconductivity, Phys
Y. Nambu, Quasi-particles and gauge invariance in the theory of superconductivity, Phys. Rev.117, 648 (1960)
work page 1960
-
[21]
A. I. Lichtenstein and M. I. Katsnelson, Antiferromag- netism and d-wave superconductivity in cuprates: A clus- ter dynamical mean-field theory, Phys. Rev. B62, R9283 (2000)
work page 2000
-
[22]
M. Caffarel and W. Krauth, Exact diagonalization ap- proach to correlated fermions in infinite dimensions: Mott transition and superconductivity, Phys. Rev. Lett. 72, 1545 (1994)
work page 1994
-
[23]
S. S. Kancharla, B. Kyung, D. S´ en´ echal, M. Civelli, M. Capone, G. Kotliar, and A.-M. S. Tremblay, Anoma- lous superconductivity and its competition with antifer- romagnetism in doped mott insulators, Phys. Rev. B77, 184516 (2008)
work page 2008
-
[24]
Civelli, Doping-driven evolution of the superconduct- ing state from a doped mott insulator, Phys
M. Civelli, Doping-driven evolution of the superconduct- ing state from a doped mott insulator, Phys. Rev. B79, 195113 (2009)
work page 2009
- [25]
- [26]
- [27]
- [28]
- [29]
-
[30]
E. Gull and A. J. Millis, Superconductivity and the pseu- dogap in the two-dimensional hubbard model, Phys. Rev. Lett.110, 216405 (2013)
work page 2013
-
[31]
X. Dong, E. Gull, and A. J. Millis, Quantifying the role of antiferromagnetic fluctuations in the superconductivity of the doped hubbard model, Nat. Phys.18, 1293 (2022)
work page 2022
-
[32]
C. Yue, H. Aoki, and P. Werner, Superconductivity en- hanced by pair fluctuations between wide and narrow bands, Phys. Rev. B106, L180506 (2022)
work page 2022
- [33]
- [34]
- [35]
- [36]
-
[37]
C. Yue, S. Hoshino, A. Koga, and P. Werner, Unconven- tional pairing from local orbital fluctuations in strongly correlatedA 3c60, Phys. Rev. B104, 075107 (2021)
work page 2021
-
[38]
C. Yue, Y. Nomura, K. Prassides, and P. Werner, In- stability of thepa 3 fulleride cs 3c60 at ambient pressure and superconducting state of the fcc phase, Phys. Rev. B108, L220508 (2023). 14
work page 2023
-
[39]
C. Yue, Y. Nomura, and P. Werner, Doping asymmetry and layer-selective metal-insulator transition in trilayer k3+xc60, Phys. Rev. Lett.129, 066403 (2022)
work page 2022
- [40]
-
[41]
B. Bacq-Labreuil, B. Lacasse, A.-M. S. Tremblay, D. S´ en´ echal, and K. Haule, Toward an ab initio theory of high-temperature superconductors: A study of multi- layer cuprates, Phys. Rev. X15, 021071 (2025)
work page 2025
- [42]
-
[43]
Y.-H. Tian, Y. Chen, J.-M. Wang, R.-Q. He, and Z.-Y. Lu, Correlation effects and concomitant two-orbitals ±- wave superconductivity in la 3ni2o7 under high pressure, Phys. Rev. B109, 165154 (2024)
work page 2024
-
[44]
S.-C. Mo, Y.-Y. Zheng, and W. W´ u, Intertwined elec- tron pairing in the bilayer two-orbital kanamori-hubbard model: A unified picture of two superconductivities, arXiv preprint (2025), 2508.04554
- [45]
- [46]
-
[47]
S. Sakai and N. Takemori, Doped mott insulator on a penrose tiling, Phys. Rev. B105, 205138 (2022)
work page 2022
-
[48]
T. Fukushima, N. Takemori, S. Sakai, M. Ichioka, and A. Jagannathan, Supercurrent distribution in real space and anomalous paramagnetic response in a supercon- ducting quasicrystal, Phys. Rev. Research5, 043164 (2023)
work page 2023
-
[49]
D. J. Luitz and F. F. Assaad, Weak-coupling continuous- time quantum monte carlo study of the single impurity and periodic anderson models withs-wave superconduct- ing baths, Phys. Rev. B81, 024509 (2010)
work page 2010
-
[50]
D. J. Luitz, F. F. Assaad, T. Novotn´ y, C. Karrasch, and V. Meden, Understanding the josephson current through a kondo-correlated quantum dot, Phys. Rev. Lett.108, 227001 (2012)
work page 2012
-
[51]
K. Satori, H. Shiba, O. Sakai, and Y. Shimizu, Numerical renormalization group study of magnetic impurities in su- perconductors, Journal of the Physical Society of Japan 61, 3239 (1992), https://doi.org/10.1143/JPSJ.61.3239
-
[52]
O. Sakai, Y. Shimizu, H. Shiba, and K. Satori, Numer- ical renormalization group study of magnetic impuri- ties in superconductors. ii. dynamical excitation spec- tra and spatial variation of the order parameter, Jour- nal of the Physical Society of Japan62, 3181 (1993), https://doi.org/10.1143/JPSJ.62.3181
- [53]
-
[54]
C. Karrasch, A. Oguri, and V. Meden, Josephson current through a single anderson impurity coupled to bcs leads, Phys. Rev. B77, 024517 (2008)
work page 2008
- [55]
- [56]
-
[57]
M. M. Wauters, C.-M. Chung, L. Maffi, and M. Burrello, Topological kondo model out of equilibrium, Phys. Rev. B108, L220302 (2023)
work page 2023
-
[58]
A. N. Rubtsov and A. I. Lichtenstein, Continuous-time quantum monte carlo method for fermions: Beyond aux- iliary field framework, J. Exp. Theor. Phys. Lett.80, 61 (2004)
work page 2004
-
[59]
A. N. Rubtsov, V. V. Savkin, and A. I. Lichten- stein, Continuous-time quantum monte carlo method for fermions, Phys. Rev. B72, 035122 (2005)
work page 2005
- [60]
-
[61]
E. Gull, P. Werner, O. Parcollet, and M. Troyer, Continuous-time auxiliary-field monte carlo for quantum impurity models, EPL82, 57003 (2008)
work page 2008
-
[62]
C.-K. Chan, P. Werner, and A. J. Millis, Magnetism and orbital ordering in an interacting three-band model: A dynamical mean-field study, Phys. Rev. B80, 235114 (2009)
work page 2009
-
[63]
E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. Troyer, and P. Werner, Continuous-time monte carlo methods for quantum impurity models, Rev. Mod. Phys. 83, 349 (2011)
work page 2011
-
[64]
I. de Vega, U. Schollw¨ ock, and F. A. Wolf, How to dis- cretize a quantum bath for real-time evolution, Phys. Rev. B92, 155126 (2015)
work page 2015
-
[65]
L. Kohn and G. E. Santoro, Efficient mapping for ander- son impurity problems with matrix product states, Phys. Rev. B104, 014303 (2021)
work page 2021
-
[66]
L. Kohn and G. E. Santoro, Quench dynamics of the an- derson impurity model at finite temperature using ma- trix product states: entanglement and bath dynamics, J. Stat. Mech. Theory Exp.2022, 063102 (2022)
work page 2022
-
[67]
D. E. Makarov and N. Makri, Path integrals for dissipa- tive systems by tensor multiplication. condensed phase quantum dynamics for arbitrarily long time, Chem. Phys. Lett.221, 482 (1994)
work page 1994
-
[68]
Makri, Numerical path integral techniques for long time dynamics of quantum dissipative systems, J
N. Makri, Numerical path integral techniques for long time dynamics of quantum dissipative systems, J. Math. Phys.36, 2430 (1995)
work page 1995
-
[69]
A. Strathearn, P. Kirton, D. Kilda, J. Keeling, and B. W. Lovett, Efficient non-markovian quantum dynamics using time-evolving matrix product operators, Nat. Commun. 9, 3322 (2018)
work page 2018
-
[70]
Gu, Efficient simulation of grassmann tensor prod- uct states, Phys
Z.-C. Gu, Efficient simulation of grassmann tensor prod- uct states, Phys. Rev. B88, 115139 (2013)
work page 2013
-
[71]
Y. Yoshimura, Y. Kuramashi, Y. Nakamura, S. Takeda, and R. Sakai, Calculation of fermionic green func- tions with grassmann higher-order tensor renormaliza- tion group, Phys. Rev. D97, 054511 (2018)
work page 2018
-
[72]
S. Akiyama and D. Kadoh, More about the grassmann tensor renormalization group, J. High Energy Phys. 2021, 188
work page 2021
-
[73]
R. Chen, X. Xu, and C. Guo, Grassmann time-evolving matrix product operators for equilibrium quantum impu- rity problems, New J. Phys.26, 013019 (2024)
work page 2024
-
[74]
R. Chen, X. Xu, and C. Guo, Grassmann time-evolving matrix product operators for quantum impurity models, Phys. Rev. B109, 045140 (2024). 15
work page 2024
-
[75]
R. Chen and C. Guo, Solving equilibrium quantum im- purity problems on the l-shaped kadanoff-baym contour, Phys. Rev. B110, 165114 (2024)
work page 2024
-
[76]
R. Chen, L. Gu, and C. Guo, Tensor network algorithm to solve polaron impurity problems, Chin. Phys. Lett.42, 120701 (2025)
work page 2025
-
[77]
Z. Sun, R. Chen, Z. Li, and C. Guo, Scalable tensor net- work algorithm for quantum impurity problems, Phys. Rev. B112, 155115 (2025)
work page 2025
-
[78]
M. Troyer and U.-J. Wiese, Computational complexity and fundamental limitations to fermionic quantum monte carlo simulations, Phys. Rev. Lett.94, 170201 (2005)
work page 2005
-
[79]
M. Eckstein, M. Kollar, and P. Werner, Interaction quench in the hubbard model: Relaxation of the spec- tral function and the optical conductivity, Phys. Rev. B 81, 115131 (2010)
work page 2010
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.