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:1911.02817 · Physical review. B./Physical review. B

Superconductor-insulator transition in Josephson junction chains by quantum Monte Carlo calculations

D. M. Basko, Frank Pfeiffer, P. Adamus, Markus Holzmann, F. W. J. Hekking

Published Jan 28, 2020• 18 citations
arXivHTMLPDF

Abstract

We study the zero-temperature phase diagram of a dissipationless and disorder-free Josephson junction chain. Namely, we determine the critical Josephson energy below which the chain becomes insulating as a function of the ratio of two capacitances: the capacitance of each Josephson junction and the capacitance between each superconducting island and the ground. We develop an imaginary-time path integral quantum Monte Carlo algorithm in the charge representation, which enables us to efficiently handle the electrostatic part of the chain Hamiltonian. We find that a large part of the phase diagram is determined by anharmonic corrections which are not captured by the standard Kosterlitz-Thouless renormalization group description of the transition.

DOI
10.1103/physrevb.101.024518

Related papers

Field-induced superconductor-to-insulator transition in Josephson-junction ladders

Enzo Granato

The superconductor-to-insulator transition is studied in a self-charging model for a ladder of Josephson junctions in the presence of an external magnetic field. Path integral Monte Carlo simulations of the equivalent $(1+1)$-dimensional classical model are used to study the phase diagram and critical behavior. In addition to a superconducting (vortex-free) phase, a vortex phase can also occur for an increasing magnetic field and small charging energy. It is found that an intervening insulating phase separates the superconducting from the vortex phases. Surprisingly, a finite-size scaling analysis shows that the field-induced superconducting-to-insulator transition is in the KT universality class, even though the external field breaks time-reversal symmetry.

arXiv:cond-mat/0509800•Physical Review B•Sep 30, 2005•14 citations

Phase transitions in the one-dimensional frustrated quantum XY model and Josephson-junction ladders

Enzo Granato

A one-dimensional quantum version of the frustrated XY (planar rotor) model is considered which can be physically realized as a ladder of Josephson junctions at half a flux quantum per plaquette. This system undergoes a superconductor to insulator transition at zero temperature as a function of charging energy. The critical behavior is studied using a Monte Carlo transfer matrix applied to the path-integral representation of the model and a finite-size-scaling analysis of data on small system sizes. Depending on the ratio between the interchain and intrachain couplings the system can have single or double transitions which is consistent with the prediction that its critical behavior should be described by the two-dimensional classical XY-Ising model.

arXiv:cond-mat/9403035•Journal of Applied Physics•May 15, 1994•7 citations

Phase diagram of Josephson junction arrays with capacitive disorder

Francesco Paolo Mancini, Pasquale Sodano, Andrea Trombettoni

We study the phase diagram at finite temperature of Josephson junction arrays with capacitive disorder (i.e., random offset charges and/or random charging energies): in the limit of large particle numbers per junction, this is a remarkable physical realization of the disordered boson Hubbard model. By using a mean-field approximation, we compute the average free energy and the equation for the phase boundary line between the insulating and superconducting phases. We find that the Mott-insulating lobe structure disappears for large variance $(\ensuremath{\sigma}\ensuremath{\gtrsim}e)$ of the offset charges probability distribution. Further, with nearest-neighbor interactions, the insulating lobe around $q=e$ is destroyed even for small values of $\ensuremath{\sigma}.$ In the case of random charging energies, until the variance of the distribution reaches some critical value the superconducting phase increases in comparison to the situation in which all self-capacitances are equal.

arXiv:cond-mat/0205568•Physical review. B, Condensed matter•Jan 21, 2003•11 citations

Response of Josephson-junction arrays near the quantum phase transition

Anne van Otterlo, Karl-Heinz Wagenblast, Rosario Fazio +1 more

We study the superconductor-insulator transition of a two-dimensional Bose-Hubbard model, considering as a specific example, an array of Josephson junctions. Within a coarse-graining approach we derive an effective free-energy functional from which we determine the phase diagram. At zero temperature it consists of a superconducting phase and Mott-insulating lobes. The phase boundaries of some of these lobes display reentrant behavior as a function of temperature. Next, we evaluate the electromagnetic-response functions of the system. The real part of the longitudinal conductivity is characterized by an excitation gap, whereas the imaginary part describes a capacitor. In an ideal system, under certain conditions a universal conductance is found at the transition. If we add low-frequency dissipation to the model a different value of the universal conductance is found, but still it is independent of the strength of the dissipation. Qualitatively differing results are obtained for frustrated and unfrustrated systems. We also discuss the Hall conductance of the system.

Physical review. B, Condensed matter•Aug 1, 1993•80 citations

Quantum Vortex Melting and Superconductor Insulator Transition in a 2D Josephson Junction Array in a Perpendicular Magnetic Field via Diffusion Monte Carlo

Pragalv Karki, Yen Lee Loh

In this study, we simulated a quantum rotor model describing a Josephson junction array (JJA) in the phase representation at zero temperature in a perpendicular magnetic field [Formula: see text] (in units of [Formula: see text]) on a [Formula: see text] square lattice with spacing a for [Formula: see text]. The superconductor-insulator transition (SIT) is tuned by the ratio of charging energy to Josephson coupling, U/J. Abrupt drops in the magnetization values were observed in the bigger lattices at certain values of B and U/J caused by the formation of vortices. Increasing U/J at a fixed B field causes quantum vortex melting. The magnetization drops to zero around [Formula: see text] indicating SIT. For B  =  0.1 the SIT occurs without an intermediate vortex state and the magnetization scales as [Formula: see text], whereas for B  =  0.4 the scaling is [Formula: see text] during the vortex melting. For B between 0.1 and 0.4 the scaling is not clear. We used the diffusion Monte Carlo (DMC) method with a guiding wavefunction optimized using the variational Monte Carlo (VMC) method. The ground state energy is calculated easily in DMC and its error estimates were generally smaller than [Formula: see text], both with and without the guiding wavefunction. Quantities like magnetization and vorticity that do not commute with the Hamiltonian were calculated using an efficient forward walking algorithm. Their estimates are affected severely in absence of the guiding wavefunction. With the guiding wavefunction, errors for the magnetization were generally less than [Formula: see text] and going up to [Formula: see text] percent around the phase transition from the Meissner to the vortex state, and without the guiding wavefunction errors were generally higher than [Formula: see text] and going up to [Formula: see text] around the critical point.

APS•Jan 1, 2019•0 citations

Quantum vortex melting and superconductor insulator transition in a 2D Josephson junction array in a perpendicular magnetic field via diffusion Monte Carlo

Pragalv Karki, Yen Lee Loh

In this study, we simulated a quantum rotor model describing a Josephson junction array (JJA) in the phase representation at zero temperature in a perpendicular magnetic field [Formula: see text] (in units of [Formula: see text]) on a [Formula: see text] square lattice with spacing a for [Formula: see text]. The superconductor-insulator transition (SIT) is tuned by the ratio of charging energy to Josephson coupling, U/J. Abrupt drops in the magnetization values were observed in the bigger lattices at certain values of B and U/J caused by the formation of vortices. Increasing U/J at a fixed B field causes quantum vortex melting. The magnetization drops to zero around [Formula: see text] indicating SIT. For B = 0.1 the SIT occurs without an intermediate vortex state and the magnetization scales as [Formula: see text], whereas for B = 0.4 the scaling is [Formula: see text] during the vortex melting. For B between 0.1 and 0.4 the scaling is not clear. We used the diffusion Monte Carlo (DMC) method with a guiding wavefunction optimized using the variational Monte Carlo (VMC) method. The ground state energy is calculated easily in DMC and its error estimates were generally smaller than [Formula: see text], both with and without the guiding wavefunction. Quantities like magnetization and vorticity that do not commute with the Hamiltonian were calculated using an efficient forward walking algorithm. Their estimates are affected severely in absence of the guiding wavefunction. With the guiding wavefunction, errors for the magnetization were generally less than [Formula: see text] and going up to [Formula: see text] percent around the phase transition from the Meissner to the vortex state, and without the guiding wavefunction errors were generally higher than [Formula: see text] and going up to [Formula: see text] around the critical point.

Journal of Physics Condensed Matter•Aug 17, 2018•1 citation
We study the zero-temperature phase diagram of a dissipationless and disorder-free Josephson junction chain. Namely, we determine the critical Josephson energy below which the chain becomes insulating as a function of the ratio of two capacitances: the capacitance of each Josephson junction and the capacitance between each superconducting island and the ground. We develop an imaginary-time path integral quantum Monte Carlo algorithm in the charge representation, which enables us to efficiently handle the electrostatic part of the chain Hamiltonian. We find that a large part of the phase diagram is determined by anharmonic corrections which are not captured by the standard Kosterlitz-Thouless renormalization group description of the transition.
The superconductor-to-insulator transition is studied in a self-charging model for a ladder of Josephson junctions in the presence of an external magnetic field. Path integral Monte Carlo simulations of the equivalent (1+1)(1+1)(1+1)-dimensional classical model are used to study the phase diagram and critical behavior. In addition to a superconducting (vortex-free) phase, a vortex phase can also occur for an increasing magnetic field and small charging energy. It is found that an intervening insulating phase separates the superconducting from the vortex phases. Surprisingly, a finite-size scaling analysis shows that the field-induced superconducting-to-insulator transition is in the KT universality class, even though the external field breaks time-reversal symmetry.
A one-dimensional quantum version of the frustrated XY (planar rotor) model is considered which can be physically realized as a ladder of Josephson junctions at half a flux quantum per plaquette. This system undergoes a superconductor to insulator transition at zero temperature as a function of charging energy. The critical behavior is studied using a Monte Carlo transfer matrix applied to the path-integral representation of the model and a finite-size-scaling analysis of data on small system sizes. Depending on the ratio between the interchain and intrachain couplings the system can have single or double transitions which is consistent with the prediction that its critical behavior should be described by the two-dimensional classical XY-Ising model.
We study the phase diagram at finite temperature of Josephson junction arrays with capacitive disorder (i.e., random offset charges and/or random charging energies): in the limit of large particle numbers per junction, this is a remarkable physical realization of the disordered boson Hubbard model. By using a mean-field approximation, we compute the average free energy and the equation for the phase boundary line between the insulating and superconducting phases. We find that the Mott-insulating lobe structure disappears for large variance $(\ensuremath{\sigma}\ensuremath{\gtrsim}e)$ of the offset charges probability distribution. Further, with nearest-neighbor interactions, the insulating lobe around q=eq=eq=e is destroyed even for small values of $\ensuremath{\sigma}.$ In the case of random charging energies, until the variance of the distribution reaches some critical value the superconducting phase increases in comparison to the situation in which all self-capacitances are equal.
We study the superconductor-insulator transition of a two-dimensional Bose-Hubbard model, considering as a specific example, an array of Josephson junctions. Within a coarse-graining approach we derive an effective free-energy functional from which we determine the phase diagram. At zero temperature it consists of a superconducting phase and Mott-insulating lobes. The phase boundaries of some of these lobes display reentrant behavior as a function of temperature. Next, we evaluate the electromagnetic-response functions of the system. The real part of the longitudinal conductivity is characterized by an excitation gap, whereas the imaginary part describes a capacitor. In an ideal system, under certain conditions a universal conductance is found at the transition. If we add low-frequency dissipation to the model a different value of the universal conductance is found, but still it is independent of the strength of the dissipation. Qualitatively differing results are obtained for frustrated and unfrustrated systems. We also discuss the Hall conductance of the system.
In this study, we simulated a quantum rotor model describing a Josephson junction array (JJA) in the phase representation at zero temperature in a perpendicular magnetic field [Formula: see text] (in units of [Formula: see text]) on a [Formula: see text] square lattice with spacing a for [Formula: see text]. The superconductor-insulator transition (SIT) is tuned by the ratio of charging energy to Josephson coupling, U/J. Abrupt drops in the magnetization values were observed in the bigger lattices at certain values of B and U/J caused by the formation of vortices. Increasing U/J at a fixed B field causes quantum vortex melting. The magnetization drops to zero around [Formula: see text] indicating SIT. For B  =  0.1 the SIT occurs without an intermediate vortex state and the magnetization scales as [Formula: see text], whereas for B  =  0.4 the scaling is [Formula: see text] during the vortex melting. For B between 0.1 and 0.4 the scaling is not clear. We used the diffusion Monte Carlo (DMC) method with a guiding wavefunction optimized using the variational Monte Carlo (VMC) method. The ground state energy is calculated easily in DMC and its error estimates were generally smaller than [Formula: see text], both with and without the guiding wavefunction. Quantities like magnetization and vorticity that do not commute with the Hamiltonian were calculated using an efficient forward walking algorithm. Their estimates are affected severely in absence of the guiding wavefunction. With the guiding wavefunction, errors for the magnetization were generally less than [Formula: see text] and going up to [Formula: see text] percent around the phase transition from the Meissner to the vortex state, and without the guiding wavefunction errors were generally higher than [Formula: see text] and going up to [Formula: see text] around the critical point.
In this study, we simulated a quantum rotor model describing a Josephson junction array (JJA) in the phase representation at zero temperature in a perpendicular magnetic field [Formula: see text] (in units of [Formula: see text]) on a [Formula: see text] square lattice with spacing a for [Formula: see text]. The superconductor-insulator transition (SIT) is tuned by the ratio of charging energy to Josephson coupling, U/J. Abrupt drops in the magnetization values were observed in the bigger lattices at certain values of B and U/J caused by the formation of vortices. Increasing U/J at a fixed B field causes quantum vortex melting. The magnetization drops to zero around [Formula: see text] indicating SIT. For B = 0.1 the SIT occurs without an intermediate vortex state and the magnetization scales as [Formula: see text], whereas for B = 0.4 the scaling is [Formula: see text] during the vortex melting. For B between 0.1 and 0.4 the scaling is not clear. We used the diffusion Monte Carlo (DMC) method with a guiding wavefunction optimized using the variational Monte Carlo (VMC) method. The ground state energy is calculated easily in DMC and its error estimates were generally smaller than [Formula: see text], both with and without the guiding wavefunction. Quantities like magnetization and vorticity that do not commute with the Hamiltonian were calculated using an efficient forward walking algorithm. Their estimates are affected severely in absence of the guiding wavefunction. With the guiding wavefunction, errors for the magnetization were generally less than [Formula: see text] and going up to [Formula: see text] percent around the phase transition from the Meissner to the vortex state, and without the guiding wavefunction errors were generally higher than [Formula: see text] and going up to [Formula: see text] around the critical point.