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:cond-mat/0611297 · Physical Review B

Cooper-pair transport through a Hubbard chain sandwiched between two superconductors: Density matrix renormalization group calculations

Adrian E. Feiguin, Steven R. White, Douglas J. Scalapino

Published Jan 11, 2007• 11 citations
arXivHTMLPDF

Abstract

We present a numerical approach to study the coherent transport of Cooper pairs through a Hubbard chain, and study the role of the contacts in achieving perfect Andreev reflection. We calculate the pair transport using the density matrix renormalization group by measuring the response of the system to quantum pair fields with complex phases on the two ends of an open system. This approach gives an effective superfluid weight which is in close agreement with the Bethe ansatz results for the superfluid weight for closed Hubbard rings.

DOI
10.1103/physrevb.75.024505

Related papers

Andreev scattering and Josephson current in a one-dimensional electron liquid

Ian K. Affleck, Jean-Sébastien Caux, Alexandre M. Zagoskin

Andreev scattering and the Josephson current through a one-dimensional interacting electron liquid sandwiched between two superconductors are reexamined. We first present some apparently new results on the noninteracting case by studying an exactly solvable tight-binding model rather than the usual continuum model. We show that perfect Andreev scattering (i.e., zero normal scattering) at the Fermi energy can only be achieved by fine-tuning junction parameters, a fine-tuning which is possible even with bandwidth mismatch between superconductor and normal metal. We also obtain exact results for the Josephson current, which is generally a smooth function of the superconducting phase difference except when the junction parameters are adjusted to give perfect Andreev scattering, in which case it becomes a sawtooth function. We then observe that, even when interactions are included, all low-energy properties of a junction $(E\ensuremath{\ll}\ensuremath{\Delta},$ the superconducting gap) can be obtained by ``integrating out'' the superconducting electrons to obtain an effective Hamiltonian describing the metallic electrons only with a boundary pairing interaction. This boundary model provides a suitable starting point for bosonization-renormalization group-boundary conformal field theory analysis. We argue that total normal reflection and total Andreev reflection correspond to two fixed points of the boundary renormalization group. For repulsive bulk interactions the Andreev fixed point is unstable and the normal one stable. However, the reverse is true for attractive interactions. This implies that a generic junction Hamiltonian (without fine-tuned junction parameters) will renormalize to the normal fixed point for repulsive interations but to the Andreev one for attractive interations. An exact mapping of our tight-binding model to the Hubbard model with a transverse magnetic field is used to help understand this behavior. We calculate the critical exponents, which are different at these two different fixed points.

arXiv:cond-mat/0001375•Physical review. B, Condensed matter•Jul 1, 2000•116 citations

Analytic approach to transport in Josephson junctions beyond the Andreev approximation: General theory and applications to the BEC-BCS crossover

F. Setiawan, Johannes Hofmann

Transport in Josephson junctions is commonly described using a simplifying assumption called the Andreev approximation, which assumes that excitations are fixed at the Fermi momentum and only Andreev reflections occur at interfaces (with no normal reflections). This approximation is appropriate for BCS-type superconductors, where the chemical potential vastly exceeds the pairing gap, but it breaks down for superconductors with low carrier density, such as topological superconductors, doped semiconductors, or superfluid quantum gases. Here, we present a generic $analytical$ framework for calculating transport in Josephson junctions that lifts up the requirement of the Andreev approximation. Using this general framework, we study in detail transport in Josephson junctions across the BCS-BEC crossover, which describes the evolution from a BCS-type superconductor with loosely-paired Cooper pairs to a BEC of tighly-paired dimers. As the interaction is tuned from the BCS to the BEC regime, we find that the overall subgap current caused by multiple Andreev reflections decreases, but nonlinearities in the current-voltage characteristic called the subharmonic gap structure become more pronounced near the intermediate unitary limit, giving rise to sharp peaks and dips in the differential conductance with even $negative$ conductance at specific voltages.

•Aug 23, 2021•0 citations

Entanglement and quantum phase transition in the asymmetric Hubbard chain: density-matrix renormalization group calculations

Chan Wen-Ling, Shi-Jian Gu

We study the ground state quantum phase transition by means of entanglement in the one-dimensional asymmetric Hubbard model with open boundary conditions. The local entanglement between the middle two sites and the rest of the system, and the block entanglement between the left and right portions of the system, are calculated using the density-matrix renormalization group (DMRG) method. We find that the entanglement shows interesting scaling and singular behavior around the phase transition line.

arXiv:0805.2021•Journal of Physics Condensed Matter•Aug 1, 2008•10 citations

Friedel oscillations in the open Hubbard chain

Gerald Bedürftig, B. Brendel, Holger Frahm +1 more

Using the density-matrix renormalization-group (DMRG) technique, we calculate critical exponents for the one-dimensional Hubbard model with open boundary conditions with and without additional boundary potentials at both ends. A direct comparison with open boundary condition Bethe ansatz calculations provides a good check for the DMRG calculations on large system sizes. On the other hand, the DMRG calculations provide an independent check of the predictions of conformal field theory, which are needed to obtain the critical exponents from the Bethe ansatz. From the Bethe ansatz we predict the behavior of the $1/L$-corrected mean value of the Friedel oscillations (for the density and the magnetization) and the characteristic wave vectors, and show numerically that these conjectures are fulfilled with and without boundary potentials. The quality of the numerical results allows us to determine the behavior of the coefficients of the Friedel oscillations as a function of the Hubbard interaction.

arXiv:cond-mat/9805123•Physical review. B, Condensed matter•Oct 15, 1998•74 citations

Symmetry of the Gap Function in High-$T_c$ Superconductors

Mihir Arjunwadkar, D. G. Kanhere, G. Baskaran +1 more

One of the most debated issues related to high-$T_c$ superconductivity is the symmetry of the Cooper pair or the gap function. In this report, we present numerical results regarding the gap function in strongly correlated electron systems using $t-J$ and Hubbard models in one and two dimensions. To this end, we use exact diagonalization to study the ground states of 8- and 16-site clusters consisting of single or coupled layers. We calculate a reduced two-particle density matrix in momentum space which is a measure of the gap function. We then analyze the eigenvectors of this density matrix, which display the possible Cooper pair symmetries. The eigenvector corresponding to the largest eigenvalue indicates a vanishing gap on the Fermi surface (which is in favour of odd-gap pairing) although $d_{x^2-y^2}$ symmetry is seen to be a very close contestant in many of the cases.

arXiv:cond-mat/9502008•Feb 2, 1995•0 citations

Transport properties of the one-dimensional Hubbard model at finite temperature

Christoph Karrasch, Dante Marvin Kennes, Joel E. Moore

We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is twofold: First, we compute both the charge- and the spin-current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ in the presence of small integrability-breaking next-nearest-neighbor interactions (the extended Hubbard model). The dc conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a test bed for our numerics, we compute ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ for the integrable XXZ spin chain in its gapped phase.

arXiv:1408.1891•Physical Review B•Oct 3, 2014•56 citations
We present a numerical approach to study the coherent transport of Cooper pairs through a Hubbard chain, and study the role of the contacts in achieving perfect Andreev reflection. We calculate the pair transport using the density matrix renormalization group by measuring the response of the system to quantum pair fields with complex phases on the two ends of an open system. This approach gives an effective superfluid weight which is in close agreement with the Bethe ansatz results for the superfluid weight for closed Hubbard rings.
Andreev scattering and the Josephson current through a one-dimensional interacting electron liquid sandwiched between two superconductors are reexamined. We first present some apparently new results on the noninteracting case by studying an exactly solvable tight-binding model rather than the usual continuum model. We show that perfect Andreev scattering (i.e., zero normal scattering) at the Fermi energy can only be achieved by fine-tuning junction parameters, a fine-tuning which is possible even with bandwidth mismatch between superconductor and normal metal. We also obtain exact results for the Josephson current, which is generally a smooth function of the superconducting phase difference except when the junction parameters are adjusted to give perfect Andreev scattering, in which case it becomes a sawtooth function. We then observe that, even when interactions are included, all low-energy properties of a junction $(E\ensuremath{\ll}\ensuremath{\Delta},$ the superconducting gap) can be obtained by ``integrating out'' the superconducting electrons to obtain an effective Hamiltonian describing the metallic electrons only with a boundary pairing interaction. This boundary model provides a suitable starting point for bosonization-renormalization group-boundary conformal field theory analysis. We argue that total normal reflection and total Andreev reflection correspond to two fixed points of the boundary renormalization group. For repulsive bulk interactions the Andreev fixed point is unstable and the normal one stable. However, the reverse is true for attractive interactions. This implies that a generic junction Hamiltonian (without fine-tuned junction parameters) will renormalize to the normal fixed point for repulsive interations but to the Andreev one for attractive interations. An exact mapping of our tight-binding model to the Hubbard model with a transverse magnetic field is used to help understand this behavior. We calculate the critical exponents, which are different at these two different fixed points.
Transport in Josephson junctions is commonly described using a simplifying assumption called the Andreev approximation, which assumes that excitations are fixed at the Fermi momentum and only Andreev reflections occur at interfaces (with no normal reflections). This approximation is appropriate for BCS-type superconductors, where the chemical potential vastly exceeds the pairing gap, but it breaks down for superconductors with low carrier density, such as topological superconductors, doped semiconductors, or superfluid quantum gases. Here, we present a generic analyticalanalyticalanalytical framework for calculating transport in Josephson junctions that lifts up the requirement of the Andreev approximation. Using this general framework, we study in detail transport in Josephson junctions across the BCS-BEC crossover, which describes the evolution from a BCS-type superconductor with loosely-paired Cooper pairs to a BEC of tighly-paired dimers. As the interaction is tuned from the BCS to the BEC regime, we find that the overall subgap current caused by multiple Andreev reflections decreases, but nonlinearities in the current-voltage characteristic called the subharmonic gap structure become more pronounced near the intermediate unitary limit, giving rise to sharp peaks and dips in the differential conductance with even negativenegativenegative conductance at specific voltages.
We study the ground state quantum phase transition by means of entanglement in the one-dimensional asymmetric Hubbard model with open boundary conditions. The local entanglement between the middle two sites and the rest of the system, and the block entanglement between the left and right portions of the system, are calculated using the density-matrix renormalization group (DMRG) method. We find that the entanglement shows interesting scaling and singular behavior around the phase transition line.
Using the density-matrix renormalization-group (DMRG) technique, we calculate critical exponents for the one-dimensional Hubbard model with open boundary conditions with and without additional boundary potentials at both ends. A direct comparison with open boundary condition Bethe ansatz calculations provides a good check for the DMRG calculations on large system sizes. On the other hand, the DMRG calculations provide an independent check of the predictions of conformal field theory, which are needed to obtain the critical exponents from the Bethe ansatz. From the Bethe ansatz we predict the behavior of the 1/L1/L1/L-corrected mean value of the Friedel oscillations (for the density and the magnetization) and the characteristic wave vectors, and show numerically that these conjectures are fulfilled with and without boundary potentials. The quality of the numerical results allows us to determine the behavior of the coefficients of the Friedel oscillations as a function of the Hubbard interaction.
Symmetry of the Gap Function in High-TcT_cTc​ Superconductors
One of the most debated issues related to high-TcT_cTc​ superconductivity is the symmetry of the Cooper pair or the gap function. In this report, we present numerical results regarding the gap function in strongly correlated electron systems using t−Jt-Jt−J and Hubbard models in one and two dimensions. To this end, we use exact diagonalization to study the ground states of 8- and 16-site clusters consisting of single or coupled layers. We calculate a reduced two-particle density matrix in momentum space which is a measure of the gap function. We then analyze the eigenvectors of this density matrix, which display the possible Cooper pair symmetries. The eigenvector corresponding to the largest eigenvalue indicates a vanishing gap on the Fermi surface (which is in favour of odd-gap pairing) although dx2−y2d_{x^2-y^2}dx2−y2​ symmetry is seen to be a very close contestant in many of the cases.
We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is twofold: First, we compute both the charge- and the spin-current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ in the presence of small integrability-breaking next-nearest-neighbor interactions (the extended Hubbard model). The dc conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a test bed for our numerics, we compute ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ for the integrable XXZ spin chain in its gapped phase.