Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure

Superconductivity induced phase-controlled mesoscopic magnetic effects in a two-dimensional electron gas that bridges two superconducting reservoirs are investigated. Giant paramagnetic response of the junction, occuring at certain values of the phase difference of the order parameter, is predicte...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2008
Автори: Romanovsky, I.A., Bogachek, E.N., Krive, I.V., Landman, U.
Формат: Стаття
Мова:English
Опубліковано: Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України 2008
Назва видання:Физика низких температур
Теми:
Онлайн доступ:https://nasplib.isofts.kiev.ua/handle/123456789/117553
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure / I.A. Romanovsky, E.N. Bogachek, I.V. Krive, U. Landman // Физика низких температур. — 2008. — Т. 34, № 10. — С. 1098–1106. — Бібліогр.: 32 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-117553
record_format dspace
spelling nasplib_isofts_kiev_ua-123456789-1175532025-02-10T01:33:36Z Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure Romanovsky, I.A. Bogachek, E.N. Krive, I.V. Landman, U. Carbon nanotubes, quantum wires and Luttinger liquid Superconductivity induced phase-controlled mesoscopic magnetic effects in a two-dimensional electron gas that bridges two superconducting reservoirs are investigated. Giant paramagnetic response of the junction, occuring at certain values of the phase difference of the order parameter, is predicted. A geometrically similar system, consisting of a graphene ribbon stretched between two superconducting leads, is also considered. The magnetic effects in this system are found to be small and the difference between the magnetic properties of the two systems is discussed. The research of I.A.R., E.N.B., and U.L. was supported by the U.S. Department of Energy, grant No.FG-05-86ER 45234. I.V.K. acknowledges the financial support from the grant «Effects of electronic, magnetic and elastic properties in strongly inhomogeneous nanostructures» by the National Academy of Sciences of Ukraine and the financial support from the Royal Swedish Academy of Sciences 2008 Article Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure / I.A. Romanovsky, E.N. Bogachek, I.V. Krive, U. Landman // Физика низких температур. — 2008. — Т. 34, № 10. — С. 1098–1106. — Бібліогр.: 32 назв. — англ. 0132-6414 PACS: 73.63.Nm;74.78.Na;74.45.+c;74.25.Bt;74.25.Ha https://nasplib.isofts.kiev.ua/handle/123456789/117553 en Физика низких температур application/pdf Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
topic Carbon nanotubes, quantum wires and Luttinger liquid
Carbon nanotubes, quantum wires and Luttinger liquid
spellingShingle Carbon nanotubes, quantum wires and Luttinger liquid
Carbon nanotubes, quantum wires and Luttinger liquid
Romanovsky, I.A.
Bogachek, E.N.
Krive, I.V.
Landman, U.
Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure
Физика низких температур
description Superconductivity induced phase-controlled mesoscopic magnetic effects in a two-dimensional electron gas that bridges two superconducting reservoirs are investigated. Giant paramagnetic response of the junction, occuring at certain values of the phase difference of the order parameter, is predicted. A geometrically similar system, consisting of a graphene ribbon stretched between two superconducting leads, is also considered. The magnetic effects in this system are found to be small and the difference between the magnetic properties of the two systems is discussed.
format Article
author Romanovsky, I.A.
Bogachek, E.N.
Krive, I.V.
Landman, U.
author_facet Romanovsky, I.A.
Bogachek, E.N.
Krive, I.V.
Landman, U.
author_sort Romanovsky, I.A.
title Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure
title_short Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure
title_full Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure
title_fullStr Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure
title_full_unstemmed Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure
title_sort giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure
publisher Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
publishDate 2008
topic_facet Carbon nanotubes, quantum wires and Luttinger liquid
url https://nasplib.isofts.kiev.ua/handle/123456789/117553
citation_txt Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure / I.A. Romanovsky, E.N. Bogachek, I.V. Krive, U. Landman // Физика низких температур. — 2008. — Т. 34, № 10. — С. 1098–1106. — Бібліогр.: 32 назв. — англ.
series Физика низких температур
work_keys_str_mv AT romanovskyia giantmagnetizationofasuperconductortwodimensionalelectrongassuperconductorstructure
AT bogacheken giantmagnetizationofasuperconductortwodimensionalelectrongassuperconductorstructure
AT kriveiv giantmagnetizationofasuperconductortwodimensionalelectrongassuperconductorstructure
AT landmanu giantmagnetizationofasuperconductortwodimensionalelectrongassuperconductorstructure
first_indexed 2025-12-02T12:18:25Z
last_indexed 2025-12-02T12:18:25Z
_version_ 1850398899934593024
fulltext Fizika Nizkikh Temperatur, 2008, v. 34, No. 10, p. 1098–1106 Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure I.A. Romanovsky1, E.N. Bogachek1, I.V. Krive2,3, and Uzi Landman1 1 School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332-0430 E-mail: igor.romanovsky@gmail.com 2 B. Verkin Institute for Low Temperature Physics and Engineering of the National Academy of Sciences of Ukraine, 47 Lenin Ave., Kharkov 61103, Ukraine 3 Department of Physics, Göteborg University, SE 41296 Göteborg , Sweden Received April 23, 2008 Superconductivity induced phase-controlled mesoscopic magnetic effects in a two-dimensional electron gas that bridges two superconducting reservoirs are investigated. Giant paramagnetic response of the junc- tion, occuring at certain values of the phase difference of the order parameter, is predicted. A geometrically similar system, consisting of a graphene ribbon stretched between two superconducting leads, is also consi- dered. The magnetic effects in this system are found to be small and the difference between the magnetic properties of the two systems is discussed. PACS: 73.63.Nm Quantum wires; 74.78.Na Mesoscopic and nanoscale systems; 74.45.+c Proximity effects; Andreev effect; SN and SNS junctions; 74.25.Bt Thermodynamic properties; 74.25.Ha Magnetic properties. Keywords: supercurrent, magnetization, Andreev levels, graphene. 1. Introduction Low-dimensional mesoscopic systems exhibit a num- ber of properties which are of fundamental scientific in- terest, as well as some that may offer opportunities for the use of such systems in future miniaturized electronic de- vices. The new behavior exhibited by systems in this size regime made them the subject of extensive theoretical and experimental research over the past two decades. Quantum wires (QW), are electric conductors with lat- eral (transverse) dimensions reduced such that the motion of the charge carriers in the wire becomes quantized. Quantum wires can be fabricated in many different ways. One of the common methods for the creation of QWs is the gate voltage method applied to semiconducting heterostructures with a two-dimensional electron gas (2DEG) at the interface [1,2]. Another often used method is the controlled break junction technique (CBJ). In this case a contact is pulled apart (for example, separating a contact formed between a tip and a surface) [3], or bent [4], in such a way that at the point where it breaks one gets a very narrow junction with a diameter comparable to the wavelength of the electrons at the Fermi level, resulting in conductance quantization [1,2,5,6] and/or force oscil- lations [3,4] that emerge upon continuous pulling of the contact. Quantum wires can also be made from quasi- two-dimensional conductors, such as graphene, by re- stricting the planar motion of the electrons in one of the two planar directions through the formation of sufficient- ly narrow ribbons [7]. Magnetic field effects in three-dimensional quantum wires connecting two superconductors have been studied previously in [8], where the phenomenon of giant magne- tization oscillations was predicted. Here we focus on me- soscopic magnetic effects of superconducting-normal-su- perconducting (SNS) constrictions, shown schematically in Fig. 1, where the narrow (normal) part that connects the two bulk superconductors is a quasi-two-dimensional conductor. We consider two cases: (i) a superconductor– two-dimensional electron gas–superconductor (S/2DEG/S) junction and (ii) a superconductor–graphene mono- layer–superconductor (S/GM/S) junction. Both types of � I.A. Romanovsky, E.N. Bogachek, I.V. Krive, and Uzi Landman, 2008 systems can be realized experimentally (see, e.g., [9–12]). In the following section we consider the S/2DEG/S system. We investigate the influence of the barriers and the carrier effective mass differences at the interfaces be- tween the semiconductor and superconductor on the Andreev levels and calculate the magnetization. In Sec. 3 we discuss the differences between the S/2DEG/S junc- tion and S/GM/S junction, and explain the smallness of the magnetic field effects in the S/GM/S junction. We summarize our results in Sec. 4. 2. Magnetic effects in S/2DEG/S junctions In this section we consider a S/N/S contact made from a semiconductor heterostructures, that bridges two bulk superconductors. The difference in the electron band structure of the InAs, In Ga Asx x1− and In Al Asx x1− semi- conductors leads to the formation of a narrow (∼ λ F ) po- tential well in the InAs layer [13]. As a result, some of the conduction electrons of the semiconductors are trapped in this well, forming a two-dimensional electron gas. The quantum wires are fabricated by restricting the lateral di- mensions of the two-dimensional electron gas in the semi- conductor heterostructure through the use of electrostatic voltage gates. The potential of the electric field of the gates forms a bottleneck whose size in the narrowest part can be reduced to reach the order of several Fermi wavelengths. In- side this narrow region motion of the electrons in the direc- tions that are transverse to the constriction axis is highly quantized and, consequently, at low temperatures only few transverse modes are populated and can conduct a current. Since the motion of the electrons in the transverse di- rection is quantized, the motion of the electrons in each transverse mode inside the 2DEG is effectively one-di- mensional. Assuming that reflection at the interfaces be- tween the semiconducting bridge and the superconduct- ing reservoirs is specular (i.e., electrons from different modes do not mix and, therefore, can be treated indepen- dently), we can describe each mode with a simplified one-dimensional model [15]. We assume that the left and right superconductors have uniform order parameters of the same magnitude Δ, but with different phases, φ1 and φ2. The order parameter inside the semiconducting bridge vanishes, i.e. Δ Δ Δ ( ) / , , / / , / . x x L L x L x L i i = < − − < < > + ⎧ ⎨ ⎪ ⎪ ⎩ ⎪ ⎪ φ φ e e 1 2 2 0 2 2 2 (1) In our model we neglect interelectron interactions inside the 2DEG. For such a quasi-1D constriction we could have used a Luttinger liquid model to describe the one-di- mensional transverse channels between the superconduc- tors (see, e.g., Refs. 16 and 17). It turns out, however, that for an adiabatic constriction both models produce the same final result. The influence of a magnetic field on the electrons in the 2DEG is twofold. First, the magnetic field acts on the magnetic moments of the electrons and holes, and second it affects the spatial motion of the charge carriers via the Lorentz force. However, if the magnetic field is applied parallel to the plane of the 2D electron gas, orbital effects can be neglected. In this article we consider only the ef- fects of Zeeman interactions. We assume that the properties of the 2DEG/S interface are very similar to the properties of a normal (N) metal/su- perconductor (S), N–S, interface, namely, the electrons and holes from the nonsuperconducting part of the con- striction, whose energies are lower than the energy gap Δ in the spectrum of the superconductors, cannot leave the normal part of our S/2DEG/S junction. Instead they are being Andreev reflected [18] at the 2DEG/S interfaces. The interference of the incident electron and the reflected hole produces a set of discrete Andreev–Kulik (AK) states [19], that are responsible for the occurence of a Josephson current through the system. The spectrum of the AK states can be found by solving the Bogolyubov– De Gennes equations [15] Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1099 Fig. 1. A schematic description of the superconductor–2DEG–su- perconductor junction. Two superconductors with different pha- ses of the order parameter are connected via a two-dimensional normal electron gas, or, alternatively, a narrow ribbon of gra- phene. A magnetic field is applied locally between the super- conductors, parallel to the plane of the 2D conductor and is taken to be negligibly weak near the SN interfaces. A similar system has been considered in Ref. 14, where the effects of Zeeman splitting and spin-orbit interaction on the Josephson current were studied. � ( ) � ( ) � * ( ) � ( ) ( ) ( ) ( H D D H u v u v x x x x x x x − ⎛ ⎝ ⎜ ⎞ ⎠ ⎟ ⎛ ⎝⎜ ⎞ ⎠⎟ = ( ) ε x ) ⎛ ⎝⎜ ⎞ ⎠⎟ , (2) where the two-component spinors u( ) ( ( ), ( ))x u x u x= ↑ ↓ and v( ) ( ( ), ( ))x v x v x= ↑ ↓ are, respectively, electron and hole probability amplitudes, and the ↑ and ↓ subscripts determine the quasiparticle spin projection. The Hamil- tonian �H is a 2 2× matrix � ( )H x he= ⎛ ⎝⎜ ⎞ ⎠⎟ + − ⎛ ⎝⎜ ⎞ ⎠⎟ 1 0 0 1 1 0 0 1 μ . (3) Here h p m x U xe F= + −2 2/ ( ) ( ) ε and μ μ= g B / 2. The first term in Eq. (3) describes the nonmagnetic part of the Hamiltonian with εF being the Fermi energy, and m x( ) the effective mass of the quasiparticles. The effective mass ms in the superconductors can be different from the effective mass mn in the semiconductor. The function U x( ) describes the external potential. We assume that the S/2DEG/S junction is clean and, therefore, U x( ) = 0 everywhere except at the interfaces. At the interfaces ( / )x L= ± 2 the difference in material composition of su- perconductor and semiconductor commonly leads to for- mation of potential barriers. We will model these poten- tial barriers by Dirac delta-functions of equal strength W , i.e. U x W x L x L( ) ( ( / ) ( / ))= + + −δ δ2 2 . (4) The second term in Eq. (3) is the Zeeman coupling of the magnetic momenta of the electrons to the external field B x( ). We assume that the magnetic field is weak ( ( )g B xBμ << Δ, where μB is the Bohr magneton and g is the Lande g-factor), and that it acts only in the nonsuper- conducting part of the junction B x B x L ( ) , | | / , , = <⎧ ⎨ ⎩ if otherwise . 2 0 These requirements are important in order to guarantee that the field does not create currents (or equivalently, phase gradients) along the S/2DEG interfaces. Finally, the off-diagonal potential �D( )x , which enters Eq. (2) is equal to � ( ) ( ) ,D x x= ⎛ ⎝⎜ ⎞ ⎠⎟ Δ 1 0 0 1 (5) where Δ( )x is the order parameter, Eq. (1). The use of the Bogolyubov–de Gennes equations to describe a S/2DEG/S structure is valid only when the nonsuperconducting seg- ment of the system is shorter than the electron–hole pair correlation decay length l v k TT F B= � / ( ) (see, e.g., dis- cussion in Ref. 20). The Andreev–Kulik levels are ob- tained by matching the solutions for the uniform regions at the S/2DEG interfaces. The result of the matching pro- cedure is quite cumbersome. To simplify it we will use the Andreev approximation [18,19]. That is, to the lowest nonvanishing order in max ( , ) /Δ E Fε we (approximate- ly) replace the wave vectors of the electrons and holes in- side the S and 2DEG parts of the junction with the Fermi wave vectors k k ke h F n≈ ≈ ( ) and q q ke h F s≈ ≈ ( ) , where k k m m F s F n s n ( ) ( ) / /= ; and for the differences between the wave vectors we write k k m ke h n F n− ≈ ε / ( ) ( ) � 2 and q q m ke h s F s− ≈ Δ / ( ) ( ) � 2 . In a long junction (L vF>> =ξ 0 � / Δ) the bound state energies close to the Fermi level, ε << Δ, are ε πω π σ μσn L Bn g B, ( ) ( ) ± φ = ± +⎛ ⎝⎜ ⎞ ⎠⎟ +Θ 2 2 , (6) where the new phase Θ is defined as Θ ≈ − φ + +⎛ ⎝ ⎜ ⎞ ⎠ ⎟π arccos R R Lk R Lk D F F Z 1 2 32 2cos ( ) cos ( ) sin ( ) , [ ( ) ] , , [( ) D m m m m Z R m m R m m m Z s n s n n s n s s = + + = = − − 4 42 2 2 1 3 2 2 2 + − + = − + 8 3 16 8 4 2 2 4 4 3 m m m m Z m Z R m m Z m m m m n n s s n n s n s s ( ) )] , [( ) nZ2 2] . (7) Here Z Wm k F= / ( )� is the parameter which describes the barrier strength and ωL Fv L= � / is the energy level spac- ing. If the effective masses of the charge carriers in the superconducting leads and in the 2DEG are equal, the re- sults obtained in Refs. 21–23 are recovered. Equation (6) describes two sets (±) of discrete levels, labeled by inte- ger indices n = ± ± …0 1 2, , , and an additional index σ = ±1 which characterizes the splitting of the energy levels in a magnetic field. Since the scattering properties of the interfaces are assumed to be spin-independent, spin is conserved in both the Andreev and normal reflections. The magnetization of each populated AK level is, accord- ing Eq. (6), g Bμ / 2. The mesoscopic harmonic factors, sin ( )2Lk F and cos ( )2Lk F in Eq. (7), are associated with interference of the incident and normally reflected quasi- particle waves. They appear due to a large change in the quasimomentum of the quasiparticle at the interface in normal reflection. In a transparent SNS junction (Z = 0) with equal effective masses, these oscillations are absent since in pure Andreev reflection processes momentum is 1100 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 I.A. Romanovsky, E.N. Bogachek, I.V. Krive, and Uzi Landman approximately conserved. As a result, Eq. (6) reduces to the well known Kulik spectrum [19]. For a transparent SNS junction in zero magnetic field, B, the AK levels, corresponding to two different sets ( ± φ), intersect at φ = +r rπ π2 (r = ± ± …0 1 2, , , ). At these special points the levels are four-fold degenerate for a single channel junction. The situation is changed if the junction is not transparent. Even if the barriers are very small, levels do not intersect and they oscillate periodi- cally with phase, approaching each other at the points φ = = ± ±πk k( , , )0 1 2 . Every AK level is now two-fold degenerate. An external magnetic field removes the re- maining degeneracy. The magnetization generated by the AK state nσ( )± of the transverse mode l is given by M n l l n , ( ) , , ( ) ( ) exp ( ( )) σ σ μσ βε ± ± = − + − φ1 . (8) The subscript l which enters this formula implies that the transverse modes are different. The total magnetization of a single transverse mode is given by a sum over all the AK states. It is useful to define the magnetizations produced by pairs of levels with opposite directions of the mag- netic moment M n M n M n l l l ( ) , ( ) ,( ) ( ) ( ) ( )± + ± −= + ± 1 1 and then sum up the pair contributions M n l ( ) ( ) ± over the index n using the Poisson summation formula M n M d M l n l l k ( ) ( ) ( ) ( ) ( ) Re ± = −∞ +∞ ± −∞ ∞ ± −∞ ∞ = +∞ ∑ ∫ ∫∑ = = +ν ν 2 1 ( ) .ν νπ νe 2 ik d ⎧ ⎨ ⎪ ⎩⎪ ⎫ ⎬ ⎪ ⎭⎪ (9) Such ordering guarantees convergence of the first and second terms in Eq. (9). The total magnetization of the single transverse mode is [8] M g B Tg k k kT l B L B L Lk = + = 2 2 1 4 2 μ πω μ ω χ π ω cos ( ) sin ( ) ( / ) Θ sinh ∞ ∑ , (10) where χ μ ω= g BB L/ . The first term in Eq. (10) does not depend on the phase difference and the temperature. This term describes magnetization of the junction in the absence of Andreev reflection; that is, for nonsuperconducting reservoirs (Pauli magnetization). It is the second term which is re- sponsible for superconductivity-induced oscillations of the magnetization. At high temperatures (T L>> ω ), when the Fermi distribution is smeared over many AK levels the amplitude of these oscillations is small in comparison with the first term in Eq. (10). The sum in Eq. (10) can be truncated at k = 1 so that the magnetization of each trans- verse mode oscillates as A cos Θ with an amplitude A g BT TL L= −( / ) exp ( / )4 22 2 2μ ω π ω . In the opposite limiting case, at low temperatures ( )T g BB< μ , only a small number of states near the Fermi energy contributes to the magnetization. The AK levels are shifted in energy because of the change of the phase difference φ between the superconductors. At certain val- ues of the phase difference some AK states approach very close to the Fermi energy. At these phase differences the magnetization of the junction will increase. Taking into account only a single pair of states which is the closest (at the given φ = φr) to the Fermi energy and neglecting the contribution of the other states, we can ap- proximate the previous formula by a simple expression M g T g T B g T g l B L B B L (osc) ≈ −⎛ ⎝⎜ ⎞ ⎠⎟ + − − + μ ω μ μ ω μ exp exp 2 2 1 2 Θ Θ B T B 2 1 ⎛ ⎝⎜ ⎞ ⎠⎟ + , (11) where the angle Θ is defined by Eq. (7). The same result can be achieved if we use the Euler–Maclaurin summa- tion formula to approximate Eq. (10) in the vicinity of the points φ = +r rπ( )1 2 . Formally Eq. (10) describes only the part of the mag- netization that is due to the discrete spectrum of AK states. The continuous spectrum (| |ε > Δ) also contributes to the magnetization [24]. At temperatures T << Δ the ef- fect of continuous spectrum is to compensate the nonsmooth contributions of the discrete spectrum arising when summing up to the finite number of bound states in- side the energy interval | |ε < Δ. In our calculation method (see Eqs. (9) and (10)) the superconducting gap was for- mally put equal to infinity, and as a result only the discrete spectrum survives. This method reproduces the correct results for thermodynamic properties of a long SNS junc- tion at temperatures T << Δ (see the corresponding discussion in a recent review [25]). In order to evaluate the total magnetization of the junc- tion we have to sum up the contributions of all the open transverse modes (channels). A channel is open if the en- ergy of its lowest longitudinal mode is smaller than the Fermi energy of the superconducting leads that it is con- nected to. Since the barriers at the S/2DEG interfaces are assumed small the motion in the longitudinal direction is almost unperturbed. This situation corresponds to the case of junctions fabricated in InGaAs heterostructures with Nb electrodes [9]. Confinement in the lateral direc- tion is usually created by an electrostatic potential gene- rated by gate electrodes etched on the surface of the semiconductor. Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1101 We obtained results for two simple models of the con- fining potential: i) a hard wall potential, and ii) a para- bolic potential. For more accurate description one needs to solve the Schr�dinger and Poisson equations self-con- sistently as described in Refs. 26 and 27. For hard wall boundary conditions, we can express the velocity of the electrons and the holes, with energies near the Fermi energy in the lth channel, in the form v m m l dF l n F n ( ) = − ⎛ ⎝⎜ ⎞ ⎠⎟ ⎛ ⎝ ⎜ ⎜ ⎞ ⎠ ⎟ ⎟ 2 1 2 2 ε π� . (12) The velocities of quasiparticles in a multichannel junc- tion with parabolic confinement U y m y /n( ) = ⊥Ω 2 2 2 are v m l F l n F ( ) = − +⎛ ⎝⎜ ⎞ ⎠⎟ ⎛ ⎝⎜ ⎞ ⎠⎟⊥ 2 1 2 ε �Ω , (13) where the lateral frequency Ω⊥ = 8 2εF nm d/ ( ) was chosen such that U y d F( / )= ± =2 ε . From this formula we see that the velocities and, con- sequently, the separation between the AK levels ω ω= =l F l v L� ( ) / are different for different transverse modes. If the 2DEG junction is wide the number of modes is large and we can replace the summation over them by an integration and write the magnetization of the S/2DEG/S junction as M N g D dB≈ ⊥ ∫μ ξ α β λ ξ ξ( , , ) ( ) 0 1 . (14) The function D f f( , , ) ( ) ( )ξ α β αξ β αξ β= − − + is the dif- ference between the occupation of the two states with op- posite momenta, f ( ) / ( )η η= +1 1e is the Fermi distribu- tion, α = �v LTFΘ / ( )2 and β μ= g B TB / ( )2 . The weight function λ ξ( ) is equal to 2ξ for the «soft walls» potential and to ξ ξ/ 1 2− for the «hard walls» potential. The integral on the right hand side of Eq. (14) cannot be calculated analytically, but we can find the asymptotic behavior of the magnetization in the vicinity of several values of the of phase difference φ. In the vicinity of the resonance points φr where Θ << LT vF/ ( )� , the magneti- zation is the highest, and it is approximately equal to M N g g B T C v LT g B B F B ≈ × × − ⎛ ⎝⎜ ⎞ ⎠⎟ ⊥ μ μ μ tanh sech ( / ) ( 4 1 12 1 2 2Θ� B T/ ) .4 ⎡ ⎣ ⎢ ⎢ ⎤ ⎦ ⎥ ⎥ (15) The constant C1 is equal to 1 2/ in the «hard wall» po- tential and to 2 3/ in the «soft wall» potential model. On the other hand, if the phase is far away from resonance ( / ( ))Θ >> LT vF� the magnetization approaches slowly the value M C N g BT L v B F ≈ ⎛ ⎝ ⎜ ⎞ ⎠ ⎟⊥2 2 2 2 2 2μ � Θ ln . (16) The constant C 2 is equal to 1 for the «hard wall» model and it takes the value 2 for the «soft wall» model. While Eq. (15) works well for transparent S/2DEG/2 junction, it may be inaccurate if the barriers at the interfaces are not small enough. In this case the condition Θ << 1 cannot be satisfied at any phase φ. In Fig. 2 we show the result of numerical calculations for the magnetization as a function of the phase differ- ence, and compare this numerical result with the asymp- totic behavior in Eqs. (15) and (16). Figure 3 shows the behavior of the magnetization at several different temper- atures. Higher temperatures lead to smearing of the Fermi distribution and the resonance peaks become smaller and broader. Since for transparent clean junction the resonance con- dition φ = φr is the same for all transverse modes, magne- tization at the resonances is proportional to the number of transverse modes. This effect is analogous to the giant oscillations of the conductance considered in Ref. 28. Therefore, the wider the junction, the higher would be the magnetization at the resonances. This behavior is illus- 1102 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 I.A. Romanovsky, E.N. Bogachek, I.V. Krive, and Uzi Landman 0.9 1.0 1.1 0 4 8 12 φ π/ M /μ B M MA MB Fig. 2. Magnetization of a transparent (Z = 0) S/2DEG/S junc- tion with a harmonic (parabolic) lateral confining potential in a magnetic field B = 10 Oe, plotted as a function of the phase difference. Also shown are its asymptotes MA (Eq. (15)) and MB (Eq. (16)). For this demonstration we choose materials with the same effective masses. We assume that the superconduc- tors are made from niobium, with a Fermi energy εF ( ) .Nb erg= ⋅ −8 52 10 12 , and effective masses m ms n= . The length of the 2DEG part of the junction is L = −10 4 cm, and it's width is d = ⋅ −6 10 5 cm. Results are shown for a temperature T = 0 1. K, which is much lower than the critical temperature of niobium Tc( ) .Nb K= 9 2 . Note that the resonance peaks are ex- tremely narrow, their width δ μ ωφ ∼ g BB L/ is approximately equal to the phase change necessary to shift the spectrum of AK levels by an amount equal to Zeeman splitting. trated in Fig. 4. On the other hand, varying the length of the junction will affect only the width of the peaks, whereas their height will remain unchanged (see Fig. 5). In the case of nontransparent junctions, the electrons and holes incident on the surface may be reflected nor- mally from the S/2DEG interface. This process modifies the spectrum of the AK levels. The AK levels do not inter- sect and do not cross the Fermi energy at any value of the phase φ. As a consequence, the amplitude of the magneti- zation peaks decreases. Differences between the effective masses of the 2DEG and the superconductors have a simi- lar effect on the magnetization as the presence of barriers at the interfaces. Namely, because of the mismatch be- tween the Fermi velocities some of the electrons are being normally reflected at the 2DEG/S surface. For larger dif- ference between the effective masses, larger fraction of all of the incident electrons will be reflected from the interfaces normally. In clean S/2DEG/S or S/N/S junctions every Andreev ref lected electron picks up an addi t ional phase δφ = ± φ +1 2, ( / )arccos E Δ at the interface, whereas nor- mal electrons do not pick up such a phase. The wave func- tion of the reflected electron-hole pair is a mixture of Andreev and normally reflected electrons and holes. The resonance conditions for this mixture are different from the resonance conditions for the wave function in the clean S/2DEG/S junction, where only Andreev reflection can occur. The resonance conditions in the junction with barriers, or with different effective masses, depend not only on the phase difference, but also on the length of the junction. Although it is always possible to achieve reso- nance for one channel, it is difficult to satisfy the reso- nance conditions for many channels simultaneously, since different transverse channels have different longitu- dinal velocities. The resonant peaks for magnetization of each single channel are very narrow, with the majority of the channels being off resonance even for small barriers or slightly different materials. As a result, the total mag- netization of the junction for both barrier models (hard and soft walls) and for different masses, is strongly sup- pressed. Therefore, for observation of strong resonance behavior the use of materials with similar effective masses is recommended. The magnetization, as well as the superconducting current in a S/2DEG/S junction, result from the difference in the population of the different AK levels. Josephson current through the junction is possible only when the en- Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1103 0.95 1.00 1.05 0 5 10 15 20 T = 0.05 K T = 0.10 K T = 0.15 K M /μ B φ π/ Fig. 3. Magnetization of the transparent S/2DEG/S junction in a magnetic field B = 10 Oe at several temperatures plotted versus the phase difference. The width of the junction is d = ⋅ −5 10 5 cm, and all the other parameters of the junction are as in Fig. 2. 14 12 10 8 6 4 2 0 M /μ B φ π/ 0.95 1.00 1.05 d = 2·10 cm –5 d = 6·10 cm –5 d = 4·10 cm –5 Fig. 4. A magnetization as a function of the width of the 2DEG part of a transparent S/2DEG junction. The length of the junc- tion is L = −10 4 cm and the temperature T = 0 1. K. The strength of the magnetic field and all other parameters are the same as in Fig. 2. 0.9 1.0 1.1 0 2 4 6 8 10 12 M /μ B φ π/ L = 2·10 cm 4– L = 3·10 cm 4– L = 1·10 cm –4 Fig. 5. The magnetization of a transparent S/2DEG/S junction versus the phase difference φ, for different lengths of the 2DEG part. The width of the junction is d = ⋅ −5 10 5 cm, temper- ature T = 0 1. K and the magnetic field B = 10 Oe. ergy levels of the two sets (dE dn / φ > 0 and dE dn / φ < 0) are unequally populated. Similarly, magnetization of the S/2DEG/S junction is a result of different population of states with opposite directions of the magnetic moments. It is most interesting to note that the magnetization is more sensitive to the barriers and mass differences than the Josephson current. Indeed, since we are constrained to use weak external magnetic fields in order to avoid de- struction of superconductivity, Zeeman splitting of the AK levels is small (μ ωB Z LB = <<Δ ). In this situation, impurities, barriers, or effective mass difference that modify the spectrum of the AK levels can move them away from the region (∼ T ) where the gradient of the Fermi distribution is high. Even a small shift Δ b of the en- ergy levels from the Fermi level may result in the situa- tion T Z b<< <<Δ Δ when both levels with opposite direc- tion of the magnetic moment are almost equally populated, combining to give a small magnetization M l B<< μ . This situation does not occur for levels which belong to different sets. These levels are well separated ∼ Δ L and small change Δ Δb L<< in their positions cannot affect their population significantly. 3. Differences between properties of S/2DEG/S and S/GM/S junctions We turn now to analysis of a junction with the geome- try similar to that of the S/2DEG/S junction dicussed above, but where the two superconducting leads are con- nected to each other through a monolayer of graphene. Graphene attracted much recent attention. In this material the dispersion relation for the low-energy electrons and holes is similar to the dispersion relation of relativistic massless fermions. This quasi-relativistic behavior of the electrons in graphene has many interesting consequences. One of these is the existence of two types of Andreev re- flections at the graphene/superconductor boundary [29]. Another consequence is the effect of Klein tunneling through the a potential barrier [30]. One might expect that Klein tunneling and specular Andreev reflection would enhance magnetic effects by reducing the strong destruc- tive interference between normally and Andreev reflected electrons and holes in the nonsuperconducting (graphene) part of the junction. In addition to specular Andreev re- flection and Klein tunneling one should also take into ac- count the fact that reflection of the electrons from the edges of the graphene ribbon depends on the orientation of the crystallographic axes with respect to the line of the edge at the sides of the graphene ribbon — this affects the quantization of the transverse motion within the ribbon [31]; in our calculations we assumed infinite-mass boundary conditions at the edges, and as a result our transverse momentum is quantized as q l Wl = +( / ) /1 2 π . Similar to the S/2DEG/S junction, by matching the so- lutions of the Schr�dinger equation for each part of the contact (see Ref. 32) one can find the spectrum of the AK levels in the S/GM/S contact and use this spectrum to cal- culate the total magnetization of the S/GM/S junction un- der the same conditions as for the S/2DEG/S one. It turns out, however, that the magnetization of the S/GM/S junc- tion is very small even in the absence of potential barriers at the interfaces between the superconductor leads and 1104 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 I.A. Romanovsky, E.N. Bogachek, I.V. Krive, and Uzi Landman 10 8 6 4 2 0 10 8 6 4 2 0 ε, ar b . u n it s ε, ar b . u n it s –1 –10 01 12 23 3 φ π/ φ π/ Fig. 6. The spectrum of AK levels in a S/2DEG/S junction (left) and in a S/GM/S junction (right). Both figures are schematic, that is, these figures are plotted for junctions with different geometrical parameters with the main purpose to illustrate the basic charac- teristic features of the energy levels in these systems. Note, that unlike the S/2DEG/S junction where the lowest AK levels of all transverse modes approach the Fermi level and become degenerate at the resonant phases (φ = +r rπ π2 , r = ± ± …0 1 2, , , ), the lowest AK levels of the S/GM/S junction, belonging to different modes, never become degenerate at the resonant phases, and all (except one level of the lowest transverse mode) stay away from the Fermi level. the graphene. In fact the magnetization is much smaller than that found for the S/2DEG/S junction (see above). There are several reasons for such a small effect. First is the fact that the density of states of graphene near the Fermi level (for undoped graphene) is very small and therefore the number of AK levels in the range of energies near the Fermi level for the S/GM/S junction should be much smaller than for the S/2DEG/S junction. The other reason for the smallness of magnetization is that the AK levels in different transverse modes of the S/GM/S junc- tion behave differently (see Fig. 6) and they do not con- verge to a single degenerate state at the Fermi level for the resonant phase differences, as happens for S/2DEG/S junctions (see above). Since the energy levels of the trans- verse modes in the S/GM/S junction never approach the Fermi level and always stay away from the region with maximal gradients of the Fermi distribution, their contri- bution to the total magnetization is very small. We conclude that the different behavior of the AK le- vels in the semiconductor 2DEG and in graphene origi- nate from the different type of dispersion relations of the electrons in the two cases: a quasi-relativistic linear dis- persion relation for the electrons in graphene versus the regular quadratic dispersion relation in the two-dimen- sional electron gas. 4. Summary In summary, we considered magnetic effects in a two-dimensional electron gas bridging two superconduct- ing reservoirs. We demonstrated that this system can ex- hibit interesting superconductivity-induced magnetic res- onance effects. Namely, we predict sharp increases in the magnetic susceptibility of the junction at special values φ = +r rπ π2 of the phase difference φ = φ −φ1 2 of the or- der parameter between the two superconductors. This effect results from a change of the population of the AK levels near the Fermi energy. In general, magneti- zation of a single transverse mode due to the Andreev lev- els is very small but since for transparent junctions the resonance conditions are identical for all transverse modes, the magnetic response at the resonances (at low temperatures) is proportional to the number of transverse modes. Consequently, when the number of transverse mo- des is large, the total magnetization of the junction may become large enough to allow experimental detection. We also considered a junction made from a graphene ribbon bridging the two superconductors. We found that the giant magnetization oscillations that we predict for the S/2DEG/S junction are absent in the Josephson con- tact made with the monolayer graphene ribbon. Reasons for the differences in magnetization response between the S/2DEG/S and S/GM/S systems were discussed. The research of I.A.R., E.N.B., and U.L. was sup- ported by the U.S. Department of Energy, grant No. FG-05-86ER 45234. I.V.K. acknowledges the financial support from the grant «Effects of electronic, magnetic and elastic proper- ties in strongly inhomogeneous nanostructures» by the National Academy of Sciences of Ukraine and the fi- nancial support from the Royal Swedish Academy of Sciences. 1. B.J. van Wees, H. van Houten, C.W.J. Beenakker, J.G. Williamson, L.P. Kouwenhoven, D. van der Marel, and C.T. Foxon, Phys. Rev. Lett. 60, 848 (1988). 2. D.A. Wharam, T.J. Thornton, R. Newbury, M. Pepper, H. Ahmed, J.E.F. Frost, D.G. Hasko, D.C. Peacock, D.A. Rit- chie, and G.A.C. Jones, J. Phys. C21, L209 (1988). 3. U. Landman, W.D. Luedtke, N. Burnham, and R.L. Colton, Science 248, 454 (1990). 4. N. Agrait, A.L. Yeyati, and J.M. van Ruitenbeek, Phys. Rep. 377, 8103 (2003). 5. E.N. Bogachek, A.M. Zagoskin, and I.O. Kulik, Fiz. Nizk. Temp. 16, 1404 (1990) [Sov. J. Low Temp. Phys. 16, 796 (1990)]. 6. A.G. Scherbakov, E.N. Bogachek, and U. Landman, Phys. Rev. B53, 4054 (1996). 7. Xiaolin Li, Xinran Wang, Li Zhang, Sangwon Lee, and Hongjie Dai, Science 29, 1229 (2008). 8. I.V. Krive, I.A. Romanovsky, E.N. Bogachek, and U. Landman, Phys. Rev. Lett. 92, 126802 (2004). 9. H. Takayanagi, T. Akazaki, and J. Nitta, Phys. Rev. Lett. 75, 3533 (1995). 10. S.G. den Hartog, B.J. van Wees, Yu.V. Nazarov, T.M. Klapwijk, and G. Borghs, Phys. Rev. Lett. 79, 3250 (1997). 11. H.B. Heersche, P. Jarillo-Herrero, J.B. Oostinga, L.M.K. Vandersypen, and A.F. Morpurgo, Nature 446, 56 (2007). 12. F. Miao, S. Wijeratne, Y. Zhang, U.C. Coskun, W. Bao, and C.N. Lau, Science 317, 1530 (2007). 13. Supriyo Datta, Electronic Transport in Mesoscopic Phy- sics, Cambridge University Press, New York (1995). 14. E.V. Bezuglyi, A.S. Rozhavsky, I.D. Vagner, and P. Wy- der, Phys. Rev. B66, 052508 (2002). 15. P.G. de Gennes, Superconductivity of Metals and Alloys, Addison-Wesley, New York (1992). 16. D.L. Maslov, M. Stone, P.M. Goldbart, and D. Loss, Phys. Rev. B53, 1548 (1996). 17. J.S. Caux, H. Saleur, and F. Siano, Phys. Rev. Lett. 88, 106402 (2002). 18. A.F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 19. I.O. Kulik, Sov. Phys. JETP 30, 944 (1970). 20. A.M. Zagoskin, Quantum Theory of Many-Body Systems, Springer, New York (1998). 21. G.A. Gogadze and A.M. Kosevich, Fiz. Nizk. Temp. 24, 716 (1998) [Low Temp. Phys. 24, 540 (1998)]. 22. O. Cakir and I.O. Kulik, Phys. Rev. B67, 174514 (2003). 23. P.F. Bagwell, Phys. Rev. B46, 12573 (1992). 24. C.W. Beenakker and H. van Houten, Phys. Rev. Lett. 66, 3056 (1991). Giant magnetization of a superconductor–two-dimensional electron gas–superconductor structure Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1105 25. I.V. Krive, S.I. Kulinich, R.I. Shekhter, and M. Jonson, Fiz. Nizk. Temp. 30, 738 (2004) [Low Temp. Phys. 30, 554 (2004)]. 26. S.E. Laux, D.J. Frank, and F. Stern, Surface Science 196, 101 (1988). 27. J.H. Davies and J.A. Nixon, Phys. Rev. B39, 3423 (1989). 28. A. Kadigrobov, A.M. Zagoskin, R.I. Shekhter, and M. Jon- son, Phys. Rev. B52, R8662 (1995). 29. C.W. Beenakker, Phys. Rev. Lett. 97, 067007 (2006). 30. M.I. Katsnelson, K.S. Novoselov, and A.K. Geim, Nature Physics 2, 620 (2006). 31. J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and C.W.J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 32. M. Titov and C.W. Beenakker, Phys. Rev. B74, 041401(R) (2006). 1106 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 I.A. Romanovsky, E.N. Bogachek, I.V. Krive, and Uzi Landman