Next Article in Journal
Study of Multi-Pixel Scintillator Detector Configurations for Measuring Polarized Gamma Radiation
Next Article in Special Issue
Reliability of the Ginzburg–Landau Theory in the BCS-BEC Crossover by Including Gaussian Fluctuations for 3D Attractive Fermions
Previous Article in Journal
Probing Electron Properties in ECR Plasmas Using X-ray Bremsstrahlung and Fluorescence Emission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Uniformly Frustrated XY Model: Strengthening of the Vortex Lattice by Intrinsic Disorder

1
Department of Physics, Stockholm University, 106 91 Stockholm, Sweden
2
ISC-CNR, Department of Physics, Sapienza University of Rome, P. le A. Moro 5, 00185 Rome, Italy
*
Author to whom correspondence should be addressed.
Condens. Matter 2021, 6(4), 42; https://doi.org/10.3390/condmat6040042
Submission received: 4 October 2021 / Revised: 3 November 2021 / Accepted: 3 November 2021 / Published: 6 November 2021

Abstract

:
In superconducting films, the role of intrinsic disorder is typically to compete with superconductivity by fragmenting the global phase coherence and lowering the superfluid density. Nonetheless, when a transverse magnetic field is applied to the system and an Abrikosov vortex lattice form, the presence of disorder can actually strengthen the superconducting state against thermal fluctuations. By means of Monte Carlo simulations on the uniformly frustrated XY model in two dimensions, we show that while for weak pinning the superconducting critical temperature T c increases with the applied field H, for strong enough pinning, the experimental decreasing dependence between T c and H is recovered with a resulting more robust vortex lattice.

1. Introduction

The melting transition of two-dimensional (2D) Abrikosov vortex lattice (VL) has always attracted significant experimental and theoretical interest. The interplay between magnetic field, random pinning, and phase fluctuations makes the phase diagram of the system rich of different phases of matter ranging from the 2D Bragg glass to the 2D vortex glass [1,2,3], from the isotropic vortex liquid to the more recently observed hexatic liquid phase [4,5]. In the absence of disorder, at low temperatures and low vortex densities, the vortex lattice exhibits a quasi-long range order. As the temperature or the vortex density increases, topological defects such as dislocations and disclinations form, eventually leading to the melting of the vortex lattice according to the Berezinskii–Kosterlitz–Thouless theory [6,7,8] afterwards refined by Nelson, Halperin and Young [9,10,11]. In this frame, the role of disorder is on the one hand to favour the formation of such lattice defects turning the ground state from a quasi ordered Bragg glass [12,13] to a disordered vortex glass. On the other hand, by acting as a pinning potential for vortices, it also prevents the VL from sliding throughout the system and destroying the global superconducting (SC) phase coherence. Thus, differently from the Meissner state, where disorder competes with superconductivity, in the mixed state, its role is less straightforward. Here, indeed, a true superconducting state exhibiting zero resistivity in the limit of zero current can only settle if a certain degree of disorder is present within the system. Moreover, as recently shown [14], vortex thermal fluctuations responsible for the fast depletion of the superfluid stiffness can be even reduced by the presence of a strong vortex pinning.
From a theoretical perspective, one possible way to investigate such interplay is by studying the classical XY model in the presence of both disorder and transverse magnetic field. As an effective model for superconducting phase fluctuations, the 2D XY model has been applied to the study of disordered and inhomogeneous SC films over the years [15,16,17,18,19,20,21,22,23,24] as well as to the 2D vortex lattice melting for different values of the applied magnetic field [25,26,27,28,29,30,31,32]. The combined effect of quenched disorder and transverse magnetic field has been so far addressed by including in the XY model spatially random quenched magnetic fluxes, but with a resulting zero net magnetic field [2,15,33,34].
In this manuscript, we will instead consider the effect of quenched disorder on the melting of the 2D vortex lattice which forms when a nonzero net transverse magnetic field is included within the XY model. Starting from the consideration that the presence of a lattice introduces a pinning effect at low temperature even in the absence of disorder, our Monte Carlo (MC) simulations show that, apart from disordering the VL ground state, the presence of inhomogeneities makes the vortex lattice stronger against thermal fluctuations with respect the homogenous case according to the experimental observations [14]. The presence of a very weak pinning can indeed lead to a very counter-intuitive observation, namely the increase of the SC critical temperature by the increase of the applied magnetic field. However, when the disorder is sufficiently strong, the experimentally observed dependence between T c and H is recovered.

2. Materials and Methods

We consider a two-dimensional XY model in the presence of random pinning and transverse magnetic field. Its Hamiltonian reads:
H X Y = μ = x ^ , y ^ , i J i μ cos ( θ i θ i + μ + F i μ ) ,
where θ i is the SC phase at site i of a L x × L y lattice, J i μ are the random couplings between the neighbouring sites i and i + μ and the phase shift F i μ accounts via the minimal substitution for the presence of a finite transverse magnetic field H z ^ = × A , being:
F i μ = 2 π Φ 0 r i r i + μ A i μ · d μ ,
with Φ 0 = h c / 2 e the quantum flux, and A i μ is the vector potential along the bond connecting two neighbour spins i and i + μ . The sum of F i μ going counterclockwise around any closed path C of such bonds is 2 π times the number of magnetic flux quanta f C passing trough the relative enclosed area defined by C:
C F i μ = 2 π Φ 0 C A · d l = 2 π Φ C Φ 0 = 2 π f C
Since in the following we will always consider H to be uniform in space, we will refer to the intensity of the applied field in terms of the flux quanta f penetrating through a unitary plaquette P:
2 π f = 2 π Φ 0 P A · d l = 2 π Φ 0 P × H · d l = 2 π Φ 0 H a 2
with a = 1 , so that f = H Φ 0 . In literature, one usually refers to this case as the uniformly frustrated XY model with frustration f. Indeed, the phase shift F i μ in the cosine argument of (1) adds frustration to the system by rendering the ground state no longer ferromagnetic: at T = 0 , the phases θ i ’s instead being all equal will vary from site to site trying to minimise the new gauge-invariant phase ( θ i θ i + μ + F i μ ) . Consequently, the value of f will correspond to the level of such frustration, determining the inhomogeneous space structure of the ground state itself. Specifically, the ground state of the uniformly frustrated XY model will consist of a periodic configuration of vortices in the phase angle θ i , whose number is directly proportional to f. The number of vortices N v of the ground state for a given value of f can be easily derived by rewriting the charge neutrality condition in terms of the new gauge invariant phase:
i P i ( θ A ) d l = 2 π i ( n i f ) = 2 π i q i = 0
where the sum runs over all the L x L y unitary plaquttes P i , so that:
N v = f L x L y .
Since in the present study we will consider periodic boundary conditions, a good gauge choice for the vector potential A is the Coulomb gauge: · A = 0 . For simplicity, we will consider:
A = B ( 0 , x ) ,
so that:
F i μ = 0 if μ = x 2 π f x i if μ = y
Finally, it is important to highlight that, with this choice, not all the values of f will be allowed. Indeed, periodic boundary conditions, together with the Coulomb gauge A y = 2 π f x , give rise to the constraint:
L x · f = 1 , 2 , 3 ,
Hence, for a given value of L x , the smallest frustration we can introduce within the system is: f = 1 / L x .
In the present work, we have studied the model Equation (1) on a square lattice with periodic boundary conditions for a linear size of L x = L y = L = 64 and different values of f. In our MC simulations, we have used a local Metropolis algorithm, needed to probe the correct canonical distribution of the system, combined with a micro-canonical over-relaxation algorithm. Specifically, each MC step consists of five Metropolis spin flips of the whole lattice, followed by 10 over-relaxation sweeps of all the spins. To help the correct thermalization of the system at lower temperatures, we have used a Simulated-Annealing procedure. For each run, we have made 50,000–75,000 MC steps, measuring the main observables with a frequency of five steps, after having discarded the first 25,000. Finally, the averages have been computed over five independent runs for the clean case and over 10 different realizations of quenched disorder for the disordered case.
In the present work, together with the ground state of the vortex lattice, we have studied the SC transition as a function of the frustration f and the level of disorder. To address this issue and measure the SC phase coherence of the system, we have computed the superfluid stiffness J s μ , which accounts for the linear response of the system to an infinitesimal twist Δ μ of the gauge-invariant phase along a given direction μ . As such, it is defined as the second derivative of the free energy with respect to Δ μ at Δ μ = 0 :
J s μ 1 L 2 2 ln Z ( Δ μ ) Δ μ 2 | Δ μ = 0
being finite in the SC phase and zero in the normal phase. Its expression for the model Equation (1) reads:
J s μ = 1 L 2 i J i μ cos ( θ i θ i + μ + F i μ ) + β L 2 i J i μ sin ( θ i θ i + μ + F i μ ) 2 i J i μ sin ( θ i θ i + μ + F i μ ) 2 ,
where β is the inverse temperature and stays both for the MC thermal average and for the average over the independent runs.

3. Results

In what follows, we will present the numerical results obtained via MC simulations both for the clean and for the disordered case, considering different values of f [ 0 , 1 2 ] .

3.1. Clean Case

Let us start with the clean case where J i μ = 1 i , μ so that:
H X Y = J μ = x ^ , y ^ , i cos ( θ i θ i + μ + F i μ ) .
For f = 0 , the model (11) is the classical XY model which undergoes to the well known Berezinskii–Kosterlitz–Thouless [6,7,8] (BKT) transition at T = T B K T . As already mentioned, for finite values of f, the phase transition is no longer driven by the unbinding of vortex–antivortex pairs, but rather by the melting of the vortex lattice which naturally forms when a transverse magnetic field is applied to the system. In the clean case, despite the absence of disorder, the vortex lattice is pinned at low temperatures by the presence of the underlying square grid that defines the array of Josephson junctions. By acting as a periodic pinning potential for vortices, such a square grid can become particularly relevant for large values of f eventually determining the symmetry of the vortex lattice itself. The most emblematic example is the limiting case f = 1 / 2 , wherein the ground state of the vortex lattice assumes a checkerboard ordered pattern. The model Equation (11) with f = 1 / 2 is also known as the fully frustrated XY model (FFXY) and as such it has been extensively studied (see [35] and reference therein) with a particular focus on its critical behavior and the nature of the phase transitions it undergoes.
In the present work, we will instead focus on a smaller value of f, where the periodic pinning potential does not induce such peculiar vortex lattice configuration. The resulting ground state vortex lattices for f = 1 / 64 ; 1 / 32 ; 1 / 8 are shown in Figure 1. For clarity, only a portion of the entire lattice is shown. For all the values of f considered, at low temperatures, the vortices form an ordered and pinned lattice whose symmetry partially depends also on the commensurability with the underlying numerical grid.
At finite but not maximal f, different kinds of phase transitions can occur [27,31], and their nature is still unclear in most of the cases. Without pretending to address this issue, in the present work, we focus on the dependence of the critical temperature T c , at which the superfluid stiffness vanishes, on the applied magnetic field f. The numerical results of J s ( T ) for different values of f are reported in Figure 2.
We can immediately notice that the critical temperature T c at which the system loses its phase coherence, resulting in the vanishing of J s , is strongly suppressed for smaller values of f, completely at odds with the usual experimental observations [2,4,5,36,37,38,39]. The general trend seems, indeed, to be a proportionality between T c and f: lower critical temperatures for smaller frustration.
The observed trend has been already discussed in [31], where Alba et al. have shown that, in the limit of small frustration f = 1 / n (and n 1 ), the critical temperature T c decreases with the increase of n: T c 1 / n 0 as n .
Apart from specific cases where the commensurability with the underlying square grid particularly strengthens or weakens the vortex-lattice structure, by lowering the applied magnetic field, the pinning of the vortex lattice becomes weaker than in the case where the vortex density is higher, with a resulting decrease of the critical temperature with f. By increasing the pinning potential via the introduction of disorder, however, the scenario is reversed and, in agreement with experimental observations, with increasing applied magnetic field, the critical temperature decreases. We will discuss the numerical results of the disordered case in the following section.

3.2. Disordered Case

We now consider the case of a transverse magnetic field applied to a disordered SC film. In the present manuscript, as in some previous works [21,22], we will use as disordered coupling constants in Equation (1) the inhomogeneous local stiffness J i μ derived from the (quantum) XY pseudo-spin 1 / 2 model in random transverse field (RTF) [40,41]:
H P S = i ξ i σ i z + J i , j ( σ i + σ j + h . c ) ,
where the random transverse fields ξ i , extracted from a uniform distribution between W and W, compete with superconductivity by mimicking the Cooper pairs localization due to disorder. The resulting local stiffness that controls the SC phase fluctuations on top of the SC ground state, read [22,41]: J i μ = J σ i x σ i + μ x , with σ i x being the mean field SC order parameter. This disorder has been shown to be appropriate to model disordered superconductors with a non-trivial spatial structure [41,42,43], as well as to account for the experimental observation of a rather broad BKT transition around the critical temperature T B K T [21,22], at which a sharp jump of the superfluid-density would be expected for zero disorder [44]. Previous studies at zero magnetic field [21,22] have shown that such a spatially-correlated disorder with large enough low-stiffness puddles is crucial to induce anomalous vortex nucleation and, consequently, a smearing of the BKT transition. On the other hand, in the present case where vortices are induced by a finite transverse magnetic field, we do not expect that the results will depend crucially on the choice of disorder, so our finding should be general also for different disorder realizations. The level of disorder is here labeled by W / J (see [21] for more details), and it is taken to be quenched in temperature.
Let us start by considering a relatively weak disorder level W / J = 4 . Consistently with the experimental observations of [14], our MC numerical results reveal that the presence of disorder leads to a modification of the ground-state vortex lattice, enlightening further its underlying mechanism. As shown in Figure 3, the core of the vortices is indeed pinned by the inhomogeneity of the local stiffness, which makes them move towards those regions with lower J i μ in order to gain in energy by minimising the Hamiltonian (1). To highlight the correlation in space between low-couplings regions and the vortex lattice deformation, in Figure 3, we have superimposed the vortex lattice to the couplings map, obtained by computing over each plaquette the average value of the local stiffness J i μ . The vortex cores, represented by red dots, are found to be in the darker areas, corresponding to the regions with the lowest local stiffness.
The presence of disorder not only modifies the vortex lattice ground state, but also impacts the superfluid response of the system as a function of the applied magnetic field. In Figure 4a, we report the temperature dependence of the superfluid stiffness J s ( T ) obtained for different values of f. As one can see, the presence of disorder restores the experimentally measured dependence between T c and f by rendering more robust the vortex lattice against thermal fluctuations. In this regard, it is quite impressive to notice that, compared to the homogeneous case (see Figure 3), the critical temperature T c corresponding to the lowest value of f = 1 / 64 in Figure 4a has increased by a factor of ten by the effects of the inhomogeneity.
The strengthening of the vortex lattice due to disorder is even more pronounced when looking at stronger disorder regimes. In Figure 4b, we report the superfluid stiffness trend in temperature for the same values of f but with a disorder level of W / J = 10 . Looking for instance at the lowest value of the field ( f = 1 / 64 ) , the critical temperature is reduced because of the field only by half with respect to the zero-field value, while, at weak disorder ( W / J = 4 ), it was five times smaller.
In order to highlight such increase of robustness, as an effect of the increase of the intrinsic disorder, we have reported in Figure 5a the extrapolated values of the critical temperature as a function of the applied field, for the two disorder regimes considered. In Figure 5b, this effect is made even more evident by rescaling the curves of T c by their value in absence of the magnetic field. The ratio T c ( f ) / T c ( 0 ) decreases with f almost twice as fast at weak disorder W / J = 4 than in the strong disorder case W / J = 10 .

4. Discussion

By means of Monte Carlo simulations on the uniformly frustrated X Y model, we have shown that the presence of disorder, mimicked via random couplings, modifies the ordered vortex lattice of the ground state. For strong enough disorder, the energy of the system is no longer minimized by an ordered vortex-lattice structure, but rather by a disordered structure in which the vortex cores are located where the local stiffness of the superfluid is lower so as to reduce the energy cost associated with a phase twist.
At the same time, we have shown that, in the presence of a very weak pinning as in the case with no disorder, where the only pinning potential is due to the presence of the underlying square grid, the increase of the applied magnetic field does not reduce the superconducting critical temperature, but, in most cases, it contributes to its increase, in stark contrast with the experimental observations [2,4,5,36,37,38,39]. The presence of disorder not only restores the experimentally observed dependence between T c and f, but also acts by making the ground-state vortex lattice more robust against thermal fluctuations. Comparing two different levels of disorder, we have indeed shown that the superconducting critical temperature is suppressed much less with respect to the zero-field case in the strong disorder regime than in the weak disorder regime.
In conclusion, despite separately both the disorder and magnetic fields acting on the SC thin film by suppressing the superfluid stiffness, when a transverse magnetic field is applied to the system, the presence of disorder helps to prevent the destruction of the SC state by thermal fluctuations.

Author Contributions

Conceptualization, I.M., L.B. and C.C.; Data curation, I.M.; Formal analysis, I.M.; Methodology, I.M.; Supervision, L.B. and C.C.; Writing—original draft, I.M.; Writing—review & editing, I.M., L.B. and C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Carl Trygger foundation through grant number CTS 20:75. We also acknowldge financial support under Ateneo 2019 (Grant No. RM11916B56802AFE), Ateneo 2020 (Grant No. RM120172A8CC7CC7) and PRIN 2017 No. 2017Z8TS5B.

Data Availability Statement

The codes needed to assess the presented results are available from the corresponding authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fisher, D.S.; Fisher, M.P.A.; Huse, D.A. Thermal fluctuations, quenched disorder, phase transitions, and transport in type-II superconductors. Phys. Rev. B 1991, 43, 130–159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Blatter, G.; Feigel’man, M.V.; Geshkenbein, V.B.; Larkin, A.I.; Vinokur, V.M. Vortices in high-temperature superconductors. Rev. Mod. Phys. 1994, 66, 1125–1388. [Google Scholar] [CrossRef]
  3. Klein, T.; Joumard, I.; Blanchard, S.; Marcus, J.; Cubitt, R.; Giamarchi, T.; Le Doussal, P. A Bragg glass phase in the vortex lattice of a type II superconductor. Nature 2001, 413, 404–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Guillamón, I.; Suderow, H.; Fernández-Pacheco, A.; Sesé, J.; Córdoba, R.; De Teresa, J.M.; Ibarra, M.R.; Vieira, S. Direct observation of melting in a two-dimensional superconducting vortex lattice. Nat. Phys. 2009, 5, 651–655. [Google Scholar] [CrossRef] [Green Version]
  5. Roy, I.; Dutta, S.; Roy Choudhury, A.N.; Basistha, S.; Maccari, I.; Mandal, S.; Jesudasan, J.; Bagwe, V.; Castellani, C.; Benfatto, L.; et al. Melting of the Vortex Lattice through Intermediate Hexatic Fluid in an a-MoGe Thin Film. Phys. Rev. Lett. 2019, 122, 047001. [Google Scholar] [CrossRef] [Green Version]
  6. Berezinsky, V.L. Destruction of Long-range Order in One-dimensional and Two-dimensional Systems Possessing a Continuous Symmetry Group. II. Quantum Systems. Sov. Phys. JETP 1972, 34, 610. [Google Scholar]
  7. Kosterlitz, J.M.; Thouless, D.J. Ordering, metastability and phase transitions in two-dimensional systems. J. Phys. C Solid State Phys. 1973, 6, 1181–1203. [Google Scholar] [CrossRef]
  8. Kosterlitz, J.M. The critical properties of the two-dimensional xy model. J. Phys. C Solid State Phys. 1974, 7, 1046. [Google Scholar] [CrossRef]
  9. Halperin, B.I.; Nelson, D.R. Theory of Two-Dimensional Melting. Phys. Rev. Lett. 1978, 41, 121–124. [Google Scholar] [CrossRef] [Green Version]
  10. Nelson, D.R.; Halperin, B.I. Dislocation-mediated melting in two dimensions. Phys. Rev. B 1979, 19, 2457–2484. [Google Scholar] [CrossRef]
  11. Young, A.P. Melting and the vector Coulomb gas in two dimensions. Phys. Rev. B 1979, 19, 1855–1866. [Google Scholar] [CrossRef]
  12. Giamarchi, T.; Le Doussal, P. Elastic theory of flux lattices in the presence of weak disorder. Phys. Rev. B 1995, 52, 1242–1270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Le Doussal, P.; Giamarchi, T. Dislocations and Bragg glasses in two dimensions. Phys. C Supercond. 2000, 331, 233–240. [Google Scholar] [CrossRef] [Green Version]
  14. Ganguly, R.; Roy, I.; Banerjee, A.; Singh, H.; Ghosal, A.; Raychaudhuri, P. Magnetic field induced emergent inhomogeneity in a superconducting film with weak and homogeneous disorder. Phys. Rev. B 2017, 96, 054509. [Google Scholar] [CrossRef] [Green Version]
  15. Barabash, S.; Stroud, D.; Hwang, I.J. Conductivity due to classical phase fluctuations in a model for high-Tc superconductors. Phys. Rev. B 2000, 61, R14924–R14927. [Google Scholar] [CrossRef] [Green Version]
  16. Wysin, G.M.; Pereira, A.R.; Marques, I.A.; Leonel, S.A.; Coura, P.Z. Extinction of the Berezinskii-Kosterlitz-Thouless phase transition by nonmagnetic disorder in planar symmetry spin models. Phys. Rev. B 2005, 72, 094418. [Google Scholar] [CrossRef]
  17. Erez, A.; Meir, Y. Effect of amplitude fluctuations on the Berezinskii-Kosterlitz-Thouless transition. Phys. Rev. B 2013, 88, 184510. [Google Scholar] [CrossRef] [Green Version]
  18. Costa, B.V.; Lima, L.S.; Coura, P.Z.; Leonel, S.A.; Lima, A.B. Kosterlitz-Thouless Transition: The Diluted XY model. J. Phys. Conf. Ser. 2014, 487, 012008. [Google Scholar] [CrossRef]
  19. Maccari, I.; Maiorano, A.; Marinari, E.; Ruiz-Lorenzo, J.J. A numerical study of planar arrays of correlated spin islands. Eur. Phys. J. B 2016, 89, 127. [Google Scholar] [CrossRef] [Green Version]
  20. Kumar, M.; Chatterjee, S.; Paul, R.; Puri, S. Ordering kinetics in the random-bond $XY$ model. Phys. Rev. E 2017, 96, 042127. [Google Scholar] [CrossRef] [Green Version]
  21. Maccari, I.; Benfatto, L.; Castellani, C. Broadening of the Berezinskii-Kosterlitz-Thouless transition by correlated disorder. Phys. Rev. B 2017, 96, 060508. [Google Scholar] [CrossRef] [Green Version]
  22. Maccari, I.; Benfatto, L.; Castellani, C. The BKT Universality Class in the Presence of Correlated Disorder. Condens. Matter 2018, 3, 8. [Google Scholar] [CrossRef] [Green Version]
  23. Maccari, I.; Benfatto, L.; Castellani, C. Disordered XY model: Effective medium theory and beyond. Phys. Rev. B 2019, 99, 104509. [Google Scholar] [CrossRef] [Green Version]
  24. Maccari, I.; Defenu, N.; Benfatto, L.; Castellani, C.; Enss, T. Interplay of spin waves and vortices in the two-dimensional XY model at small vortex-core energy. Phys. Rev. B 2020, 102, 104505. [Google Scholar] [CrossRef]
  25. Teitel, S.; Jayaprakash, C. Phase transtions in frustrated two-dimensional $\mathrm{XY}$ models. Phys. Rev. B 1983, 27, 598–601. [Google Scholar] [CrossRef]
  26. Franz, M.; Teitel, S. Vortex lattice melting in 2D superconductors and Josephson arrays. Phys. Rev. Lett. 1994, 73, 480–483. [Google Scholar] [CrossRef] [Green Version]
  27. Franz, M.; Teitel, S. Vortex-lattice melting in two-dimensional superconducting networks and films. Phys. Rev. B 1995, 51, 6551–6574. [Google Scholar] [CrossRef] [Green Version]
  28. Hattel, S.A.; Wheatley, J.M. Flux-lattice melting and depinning in the weakly frustrated two-dimensional XY model. Phys. Rev. B 1995, 51, 11951–11954. [Google Scholar] [CrossRef] [Green Version]
  29. Tanaka, A.; Hu, X. Numerical study of the flux lattice melting transition in 2D superconductors. Phys. C Supercond. 2001, 357-360, 438–441. [Google Scholar] [CrossRef]
  30. Hasenbusch, M.; Pelissetto, A.; Vicari, E. Multicritical behaviour in the fully frustrated XY Model Relat. Syst. J. Stat. Mech. Theory Exp. 2005, 2005, P12002. [Google Scholar] [CrossRef]
  31. Alba, V.; Pelissetto, A.; Vicari, E. The uniformly frustrated two-dimensionalXYmodel in the limit of weak frustration. J. Phys. A Math. Theor. 2008, 41, 175001. [Google Scholar] [CrossRef]
  32. Teitel, S. The Two-Dimensional Fully Frustrated XY Model. In 40 Years of Berezinskii–Kosterlitz–Thouless Theory; World Scientific: Singapore, 2013; pp. 201–235. [Google Scholar] [CrossRef] [Green Version]
  33. Alba, V.; Pelissetto, A.; Vicari, E. Quasi-long-range order in the 2D XY model with random phase shifts. J. Phys. A Math. Theor. 2009, 42, 295001. [Google Scholar] [CrossRef]
  34. Alba, V.; Pelissetto, A.; Vicari, E. Magnetic and glassy transitions in the square-latticeXYmodel with random phase shifts. J. Stat. Mech. Theory Exp. 2010, 2010, P03006. [Google Scholar] [CrossRef] [Green Version]
  35. Hadzibabic, Z.; Dalibard, J. BKT Physics with Two-Dimensional Atomic Gases. In 40 Years of Berezinskii–Kosterlitz–Thouless Theory; World Scientific: Singapore, 2013; pp. 297–323. [Google Scholar] [CrossRef] [Green Version]
  36. Yazdani, A.; White, W.R.; Hahn, M.R.; Gabay, M.; Beasley, M.R.; Kapitulnik, A. Observation of Kosterlitz-Thouless-type melting of the disordered vortex lattice in thin films of a-MoGe. Phys. Rev. Lett. 1993, 70, 505–508. [Google Scholar] [CrossRef]
  37. Chen, B.; Halperin, W.P.; Guptasarma, P.; Hinks, D.G.; Mitrović, V.F.; Reyes, A.P.; Kuhns, P.L. Two-dimensional vortices in superconductors. Nat. Phys. 2007, 3, 239–242. [Google Scholar] [CrossRef] [Green Version]
  38. Misra, S.; Urban, L.; Kim, M.; Sambandamurthy, G.; Yazdani, A. Measurements of the Magnetic-Field-Tuned Conductivity of Disordered Two-Dimensional Mo43Ge57 and InOx Superconducting Films: Evidence for a Universal Minimum Superfluid Response. Phys. Rev. Lett. 2013, 110, 037002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Benyamini, A.; Telford, E.J.; Kennes, D.M.; Wang, D.; Williams, A.; Watanabe, K.; Taniguchi, T.; Shahar, D.; Hone, J.; Dean, C.R.; et al. Fragility of the dissipationless state in clean two-dimensional superconductors. Nat. Phys. 2019, 15, 947–953. [Google Scholar] [CrossRef]
  40. Ma, M.; Lee, P.A. Localized superconductors. Phys. Rev. B 1985, 32, 5658–5667. [Google Scholar] [CrossRef] [PubMed]
  41. Cea, T.; Bucheli, D.; Seibold, G.; Benfatto, L.; Lorenzana, J.; Castellani, C. Optical excitation of phase modes in strongly disordered superconductors. Phys. Rev. B 2014, 89, 174506. [Google Scholar] [CrossRef] [Green Version]
  42. Ioffe, L.B.; Mézard, M. Disorder-Driven Quantum Phase Transitions in Superconductors and Magnets. Phys. Rev. Lett. 2010, 105, 037001. [Google Scholar] [CrossRef] [Green Version]
  43. Lemarié, G.; Kamlapure, A.; Bucheli, D.; Benfatto, L.; Lorenzana, J.; Seibold, G.; Ganguli, S.C.; Raychaudhuri, P.; Castellani, C. Universal scaling of the order-parameter distribution in strongly disordered superconductors. Phys. Rev. B 2013, 87, 184509. [Google Scholar] [CrossRef] [Green Version]
  44. Nelson, D.R.; Kosterlitz, J.M. Universal Jump in the Superfluid Density of Two-Dimensional Superfluids. Phys. Rev. Lett. 1977, 39, 1201–1205. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Vortex lattice ground state for different values of the frustration f = (a) 1 / 64 ; (b) 1 / 32 ; (c) 1 / 8 . Each vortex core is plotted as a blue point. The linear size of the simulated system is L = 64 , while for the sake of clarity in the figure only a portion of the whole VL is shown.
Figure 1. Vortex lattice ground state for different values of the frustration f = (a) 1 / 64 ; (b) 1 / 32 ; (c) 1 / 8 . Each vortex core is plotted as a blue point. The linear size of the simulated system is L = 64 , while for the sake of clarity in the figure only a portion of the whole VL is shown.
Condensedmatter 06 00042 g001
Figure 2. Superfluid stiffness J s as a function of the temperature T for different values of frustration f in the absence of disorder. The linear size of the simulated system is L = 64 .
Figure 2. Superfluid stiffness J s as a function of the temperature T for different values of frustration f in the absence of disorder. The linear size of the simulated system is L = 64 .
Condensedmatter 06 00042 g002
Figure 3. Vortex lattice ground state for different values of the frustration f = (a) 1 / 8 ; (b) 1 / 16 ; (c) 1 / 32 ; (d) 1 / 64 superimposed to a map of a given disorder realization with W / J = 4 . The gray scale refers to the the average value of the local couplings constant J i μ around each single plaquette. Each vortex core is plotted as a red point. The linear size of the simulated system is L = 64 , while, for the sake of clarity in the figure, only a portion of the whole VL is shown.
Figure 3. Vortex lattice ground state for different values of the frustration f = (a) 1 / 8 ; (b) 1 / 16 ; (c) 1 / 32 ; (d) 1 / 64 superimposed to a map of a given disorder realization with W / J = 4 . The gray scale refers to the the average value of the local couplings constant J i μ around each single plaquette. Each vortex core is plotted as a red point. The linear size of the simulated system is L = 64 , while, for the sake of clarity in the figure, only a portion of the whole VL is shown.
Condensedmatter 06 00042 g003
Figure 4. Temperature dependence of the superfluid stiffness in the presence of (a) weak W / J = 4 and (b) strong disorder W / J = 10 , for different values of the frustration f. The continuous gray line is the BKT critical line 2 T / π relative to the case f = 0 .
Figure 4. Temperature dependence of the superfluid stiffness in the presence of (a) weak W / J = 4 and (b) strong disorder W / J = 10 , for different values of the frustration f. The continuous gray line is the BKT critical line 2 T / π relative to the case f = 0 .
Condensedmatter 06 00042 g004
Figure 5. Magnetic field dependence of (a) the critical temperature T c ( f ) and (b) of the rescaled critical temperature T c ( f ) / T c ( f = 0 ) for the two regimes of disorder considered: W / J = 4 and W / J = 10 . For each values of f 0 and disorder level, the critical temperature T c has been estimated from the superfluid stiffness trend in temperature J s ( T ) as the temperature at which J s vanishes starting from low temperatures. For the case f = 0 , we used instead the Nelson–Kosterlitz universal relation [44].
Figure 5. Magnetic field dependence of (a) the critical temperature T c ( f ) and (b) of the rescaled critical temperature T c ( f ) / T c ( f = 0 ) for the two regimes of disorder considered: W / J = 4 and W / J = 10 . For each values of f 0 and disorder level, the critical temperature T c has been estimated from the superfluid stiffness trend in temperature J s ( T ) as the temperature at which J s vanishes starting from low temperatures. For the case f = 0 , we used instead the Nelson–Kosterlitz universal relation [44].
Condensedmatter 06 00042 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maccari, I.; Benfatto, L.; Castellani, C. Uniformly Frustrated XY Model: Strengthening of the Vortex Lattice by Intrinsic Disorder. Condens. Matter 2021, 6, 42. https://doi.org/10.3390/condmat6040042

AMA Style

Maccari I, Benfatto L, Castellani C. Uniformly Frustrated XY Model: Strengthening of the Vortex Lattice by Intrinsic Disorder. Condensed Matter. 2021; 6(4):42. https://doi.org/10.3390/condmat6040042

Chicago/Turabian Style

Maccari, Ilaria, Lara Benfatto, and Claudio Castellani. 2021. "Uniformly Frustrated XY Model: Strengthening of the Vortex Lattice by Intrinsic Disorder" Condensed Matter 6, no. 4: 42. https://doi.org/10.3390/condmat6040042

APA Style

Maccari, I., Benfatto, L., & Castellani, C. (2021). Uniformly Frustrated XY Model: Strengthening of the Vortex Lattice by Intrinsic Disorder. Condensed Matter, 6(4), 42. https://doi.org/10.3390/condmat6040042

Article Metrics

Back to TopTop