Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1007.1609 · Physical Review Letters

Efficient Numerical Approach to Inhomogeneous Superconductivity: The Chebyshev-Bogoliubov–de Gennes Method

Lucian Covaci, François Maria Peeters, Mona Berciu

Published Oct 12, 2010• 112 citations
arXivHTMLPDF

Abstract

We propose a highly efficient numerical method to describe inhomogeneous superconductivity by using the kernel polynomial method in order to calculate the Green's functions of a superconductor. Broken translational invariance of any type (impurities, surfaces, or magnetic fields) can be easily incorporated. We show that limitations due to system size can be easily circumvented and therefore this method opens the way for the study of scenarios and/or geometries that were unaccessible before. The proposed method is highly efficient and amenable to large scale parallel computation. Although we only use it in the context of superconductivity, it is applicable to other inhomogeneous mean-field theories.

DOI
10.1103/physrevlett.105.167006

Related papers

Efficient Numerical Self-Consistent Mean-Field Approach for Fermionic Many-Body Systems by Polynomial Expansion on Spectral Density

Yuki Nagai, Yukihiro Ota, Masahiko Machida

We propose an efficient numerical algorithm to solve Bogoliubov de Gennes equations self-consistently for inhomogeneous superconducting systems with a reformulated polynomial expansion scheme. This proposed method is applied to typical issues such as a vortex under randomly distributed impurities and a normal conducting junction sandwiched between superconductors. With various technical remarks, we show that its efficiency becomes remarkable in large-scale parallel performance.

arXiv:1105.4939•Journal of the Physical Society of Japan•Jan 20, 2012•50 citations

Spectral-density polynomial expansion for two-particle Green's function in BdG framework

Yuki Nagai, Y. Ota, Masahiko Machida

We develop an effcient numerical method on the basis of the spectral-density polynomial expansion to calculate two particle Green's functions in Bogoliubov-de Gennes framework for inhomogeneous superconducting systems. Our method allows systematic evaluation such as spin or charge susceptibility, nuclear magnetic relaxation rate, optical or thermal conductivity even in magnetic fields. In addition, its effciencybecomes remarkable in large-scale parallel performance.

Physics Procedia•Jan 1, 2012•4 citations

N-independent Localized Krylov–Bogoliubov-de Gennes Method: Ultra-fast Numerical Approach to Large-scale Inhomogeneous Superconductors

Yuki Nagai

We propose the ultra-fast numerical approach to large-scale inhomogeneous superconductors, which we call the Localized Krylov Bogoliubov-de Gennes method (LK-BdG). In the LK-BdG method, the computational complexity of the local Green's function, which is used to calculate the local density of states and the mean-fields, does not depend on the system size $N$. The calculation cost of self-consistent calculations is ${\cal O}(N)$, which enables us to open a new avenue for treating extremely large systems with millions of lattice sites. To show the power of the LK-BdG method, we demonstrate a self-consistent calculation on the 143806-site Penrose quasicrystal lattice with a vortex and a calculation on 1016064-site two-dimensional nearest-neighbor square-lattice tight-binding model with many vortices. We also demonstrate that it takes less than 30 seconds with one CPU core to calculate the local density of states with whole energy range in 100-millions-site tight-binding model.

arXiv:2001.02362•Journal of the Physical Society of Japan•Jun 8, 2020•26 citations

Numerical Construction of a Low-Energy Effective Hamiltonian in a Self-Consistent Bogoliubov–de Gennes Approach of Superconductivity

Yuki Nagai, Yasushi Shinohara, Yasunori Futamura +2 more

We propose a fast and efficient approach for solving the Bogoliubov-de Gennes (BdG) equations in superconductivity, with a numerical matrix-size reduction procedure proposed by Sakurai and Sugiura [J. Comput. Appl. Math. 159, 119 (2003)]. The resultant small-size Hamiltonian contains the information of the original BdG Hamiltonian in a given energy domain. In other words, the present approach leads to a numerical construction of a low-energy effective theory in superconductivity. The combination with the polynomial expansion method allows a self-consistent calculation of the BdG equations. Through numerical calculations of quasi-particle excitations in a vortex lattice, thermal conductivity, and nuclear magnetic relaxation rate, we show that our approach is suitable for evaluating physical quantities in a large-size superconductor and a nano-scale superconducting device, with the mean-field superconducting theory.

arXiv:1303.3683•Journal of the Physical Society of Japan•Aug 2, 2013•21 citations

Efficient numerical method for evaluating normal and anomalous time-domain equilibrium Green's functions in inhomogeneous systems

Tomas Löthman, Christopher Triola, Jorge Cayao +1 more

In this work we develop EPOCH (equilibrium propagator by orthogonal polynomial chain), a computationally efficient method to calculate the time-dependent equilibrium Green's functions, including the anomalous Green's functions of superconductors, to capture the time evolution in large inhomogeneous systems. The EPOCH method generalizes the Chebyshev wave-packet propagation method from quantum chemistry and efficiently incorporates the Fermi-Dirac statistics that is needed for equilibrium quantum condensed matter systems. The computational cost of EPOCH scales only linearly in the system degrees of freedom, generating an extremely efficient algorithm also for very large systems. We demonstrate the power of the EPOCH method by calculating the time evolution of an excitation near a superconductor--normal metal interface in two and three dimensions, capturing transmission as well as normal and Andreev reflections.

arXiv:2004.10783•Physical review. B./Physical review. B•Sep 2, 2021•13 citations

Self-consistent calculations of the effects of disorder in d-wave and s-wave superconductors

Long He, Yun Song

The kernel polynomial method, as an approach to deriving highly-accurate spectral properties of large sparse matrices, is extended to investigate the effects of disorder in two-dimensional inhomogeneous superconductors. Being able to solve the Bogoliubov-de Gennes equations for quite large square lattices self-consistently, we gain new microscopic insights and make intuitive observations of the localization of low-energy quasiparticles on the nanometer scale. We also compare the behaviors of the optical conductivities of superconductors with different pairing symmetries. We find that the Drude weight turns out to be finite only for the weakly-disordered d -wave superconductor while obvious suppression of the optical conductivity in the infrared range can be observed in the disordered s -wave superconductor.

Journal of the Korean Physical Society•Jun 1, 2013•4 citations
We propose a highly efficient numerical method to describe inhomogeneous superconductivity by using the kernel polynomial method in order to calculate the Green's functions of a superconductor. Broken translational invariance of any type (impurities, surfaces, or magnetic fields) can be easily incorporated. We show that limitations due to system size can be easily circumvented and therefore this method opens the way for the study of scenarios and/or geometries that were unaccessible before. The proposed method is highly efficient and amenable to large scale parallel computation. Although we only use it in the context of superconductivity, it is applicable to other inhomogeneous mean-field theories.
We propose the ultra-fast numerical approach to large-scale inhomogeneous superconductors, which we call the Localized Krylov Bogoliubov-de Gennes method (LK-BdG). In the LK-BdG method, the computational complexity of the local Green's function, which is used to calculate the local density of states and the mean-fields, does not depend on the system size NNN. The calculation cost of self-consistent calculations is O(N){\cal O}(N)O(N), which enables us to open a new avenue for treating extremely large systems with millions of lattice sites. To show the power of the LK-BdG method, we demonstrate a self-consistent calculation on the 143806-site Penrose quasicrystal lattice with a vortex and a calculation on 1016064-site two-dimensional nearest-neighbor square-lattice tight-binding model with many vortices. We also demonstrate that it takes less than 30 seconds with one CPU core to calculate the local density of states with whole energy range in 100-millions-site tight-binding model.
We propose a fast and efficient approach for solving the Bogoliubov-de Gennes (BdG) equations in superconductivity, with a numerical matrix-size reduction procedure proposed by Sakurai and Sugiura [J. Comput. Appl. Math. 159, 119 (2003)]. The resultant small-size Hamiltonian contains the information of the original BdG Hamiltonian in a given energy domain. In other words, the present approach leads to a numerical construction of a low-energy effective theory in superconductivity. The combination with the polynomial expansion method allows a self-consistent calculation of the BdG equations. Through numerical calculations of quasi-particle excitations in a vortex lattice, thermal conductivity, and nuclear magnetic relaxation rate, we show that our approach is suitable for evaluating physical quantities in a large-size superconductor and a nano-scale superconducting device, with the mean-field superconducting theory.
In this work we develop EPOCH (equilibrium propagator by orthogonal polynomial chain), a computationally efficient method to calculate the time-dependent equilibrium Green's functions, including the anomalous Green's functions of superconductors, to capture the time evolution in large inhomogeneous systems. The EPOCH method generalizes the Chebyshev wave-packet propagation method from quantum chemistry and efficiently incorporates the Fermi-Dirac statistics that is needed for equilibrium quantum condensed matter systems. The computational cost of EPOCH scales only linearly in the system degrees of freedom, generating an extremely efficient algorithm also for very large systems. We demonstrate the power of the EPOCH method by calculating the time evolution of an excitation near a superconductor--normal metal interface in two and three dimensions, capturing transmission as well as normal and Andreev reflections.
The kernel polynomial method, as an approach to deriving highly-accurate spectral properties of large sparse matrices, is extended to investigate the effects of disorder in two-dimensional inhomogeneous superconductors. Being able to solve the Bogoliubov-de Gennes equations for quite large square lattices self-consistently, we gain new microscopic insights and make intuitive observations of the localization of low-energy quasiparticles on the nanometer scale. We also compare the behaviors of the optical conductivities of superconductors with different pairing symmetries. We find that the Drude weight turns out to be finite only for the weakly-disordered d -wave superconductor while obvious suppression of the optical conductivity in the infrared range can be observed in the disordered s -wave superconductor.