# Effect of tilted magnetic field on magnetosubbands and conductance of bi-layer quantum wire

## Abstract

The single electron magnetotransport in a vertical bi-layer semiconductor nanowire made of InAlAs/InGaAs and AlGaAs/GaAs heterostructure is theoretically studied. The magnetic field is directed perpendicularily to the main (transport) axis of the quantum wire and both non-zero components of magnetic field, that is the transverse and the vertical ones, allow to change the magnitude of intra-layer and inter-layer subbands mixing, respectively. We analyze in detail the changes introduced to energy dispersion relation E(k) by strong titled magnetic field up to several teslas for a symmetric and an asymmetric confining potential in the growth direction. These calculated energy dispersion relations are thereafter used to show that the value of conductance of bi-layer nanowire may jump as well as drop by few conductance quanta when the Fermi energy is changed what in conjunction with spin Zeeman effect may give a moderately spin polarized current.

###### pacs:

72.25.Dc,73.21.Hb## I Introduction

The electron transport properties of a nanosystem consisting of two quantum wires being placed close to each other depends largely on the magnitude of their tunnel coupling.Eugster and del Alamo (1991); Bierwagen et al. (2003); Lyo (1996); Mourokh et al. (2007); Fischer et al. (2006a) For strong and moderate coupling strength the electron’s wave functions originating from separate wires can hybridize and the magnitude of hybridization depends naturally on the energy splitting of the electron’s subbands and the electron’s wave functions overlap. The last two quantities are particularily sensitive to the strength of magnetic field and its direction.Huang et al. (2008); Fischer et al. (2005) If magnetic field is directed along the main axis of a coupled wires system, it squeezes the electron’s wave functions what leads to their stronger localization within particular nanowire and in consequence it weakens the effect of subbands hybridization.Thomas et al. (1999); Fischer et al. (2005) However, if direction of magnetic field is set perpendicularily to the axis of a nanosystem, then the Lorentz force pushes the electrons to the edges of the nanowires. The wave functions being pushed to the central barrier separating the wires can be in such case easily mixed. Such inter-wire subbands mixing by perpendicular magnetic field transforms the crossings of hybridized subbands in the energy dispersion relation E(k) into anticrossings.Shi and Gu (1997); Lyo (1996) Inside an anticrossing a pseudogap is formed what means that the conductance of two coupled nanowires is lowered when the Fermi energy level is raised and enters this region. Qualitatively mechanism of the pseudogaps formation in E(k) spectrum in magnetic field is the same for two laterally or vertically coupled nanowires systems. However in practical applications a nanosystem built of two vertically coupled wires has advantage over a system with two wires aligned laterally since: i) it allows to tune the Fermi energy in particular wire independently of its value in the second one by means of the top, bottom and side gates Bielejec et al. (2005); Kumar et al. (2014); Buchholz et al. (2008); Hew et al. (2009); Fischer et al. (2006a) and, ii) each of the three components of magnetic field modifies in a different manner the conductance of a nanosystemFischer S. F. et al. (2006); Mourokh et al. (2007) giving thus an opportunity to control over a single subband.Fischer S. F. et al. (2006)

In this paper, we theoretically analyze the effect of tilted magnetic field on a single electron transport in a bi-layer quantum wire made of a double inverted heterostructure like InAlAs/InGaAs and AlGaAs/GaAs. We assume the electron current flows in the wire (along x axis) without scattering i.e. the electron transport is ballistic, layers are vertically aligned in z (growth) direction one over another and the surrounding confining potential is formed by the rectangular barriers. Direction of magnetic field is perpendicular to the wire axis, that is only its vertical () and transverse () components have non-zero values and both can be changed independently with precision.Tarasov et al. (2010) Such approach give us an ability to tune an interlayer and intralayer modes coupling by changing and , respectively. First, we discuss how the simultaneous mixing of the vertical and the transverse eigenstates for the assumed double-well confining potential modifies the energy dispersion relation of bi-layer nanowire. Next the oscillating behaviour of magnetosubbands in vicinity of is considered in context of the conductance variations as function of Fermi energy. We show that for a nanosystem working in a ballistic regime these oscillations give contribution to the wire conductance which may jump as well as drop by few conductance quanta when Fermi energy level is successively raised or lowered between two neighbouring pseudogaps. In the last part of this paper we discuss a possibility of application of bilayer quantum wire as a source of partly spin polarized current for moderate Fermi energies.

Paper is organized as follows. In Sec.II we present theoretical model used in calculations, properties of magnetosubbands for tilted magnetic field are discussed in Sec.III while in Sec.IV we present the potential application of pseudogaps appearing in energy spectrum for partial spin polarization of conductance. We end up with conclusions given in Sec.V.

## Ii Theoretical model

The confining potential in a conventional semiconductor quantum wire can be formed electrostatically by gating the 2DEGHew et al. (2009), etching of nano-groovesFischer et al. (2006a) on the layered nanostructure that holds 2DEG few tens of nanometers beneath the surface or by cleaved-edge overgrowth.Huber et al. (2005) In all cases the 2DEG is formed within a square well created by double inverted heterojunction.Fischer et al. (2006a) The first three methods give soft lateral confining potential which was widely used in theoretical works beforeAbdullah et al. (2010); Ihnatsenka and Zozoulenko (2006) while the last one generates rectangular confinement that we have adopted for this work. We consider a quantum wire in which the electrons can move freely along x axis but their motion in y-z plane is quantized due to the rectangle shapes of external barriers. The quantum well is defined for and with high confining potential outside this region. Throughout this paper we use and . For simplicity we assume the barrier surrounding the wire is infinite while the confining potential inside the channel depends only on position in the growth direction (z-axis) i.e. . We model the confining potential by formula:

(1) |

which describes the potential with maximum localized in a central region of original well [non-zero parameter breaks the symmetry in ]. Such potential is formed within a wide quantum well when one delta doping layer is placed below and another one above the well. Then, the positively ionized dopants effectively lower the confining potential near both edges of the well giving thus a bi-layer coupled system within a single nanowire.Fischer et al. (2005); Schnez et al. (2011) Depth of an upper well, or in other words the value of parameter , can be adjusted by changing the voltage applied to the central top gateThomas et al. (1999) which may cover the whole structureHackens et al. (2002) or to the top split gate.Fischer et al. (2005)

An example of a such confining potential is showed in Fig.1(a). A bi-layer system can be also formed by stacking two quantum wires one above the other during the epitaxial growth with very narrow tunnel barrier separating them.Fischer et al. (2006a); Fischer S. F. et al. (2006); Buchholz et al. (2008) Since all effects, we investigate here, depend mainly on the energy difference between two lowest eigenstates of vertical quantization, an actual shape of vertical confinement is of little importance and the results presented below are representative for both types of confinements. Both layers, the upper and the lower one are pierced by magnetic field which has whereas the values of two other components can be freely changed. Since magnetic field lifts spin degeneracy, we limit our considerations to subbands with spins set parallel to magnetic field until otherwise stated. The energies for subbands with antiparallel spins can be simply obtaned by adding spin-Zeeman splitting energy . In calculations we use a non-symmetric vector potential for which the single electron Hamiltonian reads:

(2) |

where: is an imagnary unit, is an electron’s charge and is its effective mass ( for GaAs and for InGaAs). The eigenstates of Hamiltonian given by Eq. 2) can be expressed as the linear combination of products of the plane waves for x direction and the eigenstates for the transverse (y-axis) and vertical (z-axis) directions. We define the electron’s wave function for p-th subband and the wave vector k as follows:

(3) |

In Eq. 3, and are the basis functions for quantization in z and y directions, respectively. Due to infinite barriers surrounding the quantum wire, the normalized basis functions for y direction are simply:

(4) |

The basis functions have been found by solving the eigenproblem for z direction i.e. with Hamiltonian . The coefficients of linear combination appearing in Eq. 3 are real numbers which are strictly 0 or 1 when magnetic mixing is absent () otherwise they fulfill condition .

The maximum in the confining potential divides the original well into two narrower tunnel-coupled wells which for have the same eigenenergies. Figure 1(b) shows the probability densities for the ground [] and the first excited [] eigenstates of Hamiltonian . We notice that for large amplitude of [see Eq.1], the densities are localized mainly in two narrower wells but mutually overlap to some extent. Since energy difference for these eigenstates [] is small, even a small distortion in the confinement [ in Figs.1(b)] can significantly mix them what may result in their spatial separation [red curves in Fig.1(b)]. Therefore, these two thin wells form two transport layers and the nanowire becomes in fact a bi-layer system. To get deeper insight into nature of this process, we calculate as function of . This dependency is showed in Fig.2(a). In spite of value, decreases if grows but the lowest we get always for . For , the ground state is localized in the upper well what minimizes its energy while the first excited one occupies the lower well because of orthogonality constrictions but for the price of its increased energy.Hackens et al. (2002) Due to a large energy separation of higher eigenstates () we have neglected them in calculations and use in Eq.3 while the number of transverse modes was limited to .

Having the wave functions defined by Eq.3, we may used them to transform an original Hamiltonian given by Eq.2 to much simpler algebraic form. For this puropse, we first eliminate x variable from Hamiltonian by averaging it over this variable: . Reduced Hamiltonian then reads:

(5) |

Next, we calculate and thereafter . The effective Hamiltonian has the following form:

(6) |

(7) | |||||

(8) | |||||

(9) | |||||

(10) |

In above array of equations we have assumed: ; , ; ; ; ; ; and is a Kronecker’s delta. Indices stand for transverse modes (defined in Eq.4) while pair of mark the vertical ones. The effective Hamiltonian has matrix form which for only two-element basis becomes block matrix with real elements:

(11) |

Since all terms appearing in Eqs.7-10 are real, we immediately get condition . The only term depending on is and therefore it contributes to diagonal submatrices and (their rank equals N that is the number of basis states ). In consequence does not mix z-eigenstates. If confining potential is symmetric () then the wave function has even parity (ground or bounding state) while has odd parity (first excited or antibonding state). Since the term has itself even parity, the matrix elements give contributions only to diagonal submatrices, otherwise () they contribute also to off-diagonal submatrix . On the other hand, term has odd parity for . It mixes vertical states and and in consequence gives contribution to due to and terms. Notice, that for these Hamiltonians give contributions to diagonal matrices and .

## Iii Results

### iii.1 Two-state model for

For there is no mixing between basis states since the term disappears in while term reduces to due to [see Eqs.7-10]. This gives us possibility to limit our considerations for a moment to the case with keeping in mind that energy dispersion for transverse modes with indcies are strictly replicas of that with . Then, only diagonal elements of are shifted towards higher energy in the same manner. The two-state Hamiltonian now reads:

(12) | |||||

(13) |

Eigenvalues of can be written as:

(14) |

where we have used following abbrievations: ; ; . Energy dispersion relation for these two subbands are displayed in Fig.3(a) (black color). In first subband there are three extremums: two minimums separated by maximum localized at . Localization of minimums can be found by imposing condition on dispersion relation what gives:

(15) |

Value of depends on and and therefore, for a fixed geometry and values of material parameters of nanowire, the non-zero value of depends only on due to requirement of non-negativity of expression in square root in Eq.15. Minimal value of that gives can be estimated from the following formula:

(16) |

For , when energy difference is determinated by value of , the increase of decreases value of and consequently smaller is needed for lateral minima to appear. Such dependence is displayed in Fig.2(b) which shows a difference between maximum and minimum of energy for first subband in function of value. Simply, the lower is, the lower is needed for this difference to have non-zero value. In close vicinity of the squared term in Eq.15 can be expanded in power series of k. Neglecting the terms with exponents greater than 3 which are small, we get the parabolic shape of energy dispersion:

(17) |

where sign corresponds to upper (lower) subband. For , the energy difference between both subbands equals and depends on the sum of two differences: i) and, ii) . Let us notice that, for increasing value of both basis states for vertical direction i.e. and becomes degenerated [Fig.2(a)] with similar densities [Fig.1(b)]. For this reason, the energy difference between two lowest subbands gets smaller when is increased for .

Appearance of two additional energy minimums in magnetosubbands brings severe consequences for electron transport. If is large enough to create energy minima then the kinetic wave vector becomes negative for . Similarly, due to symmetry of a confinig potential, for negative k in the range , the effective wave vector becomes positive. Any oscillations in energy depending on canonical wave vector k produce thus negative energy dispersion relations. To study this problem in detail we have plotted the energies of two lowest subbands in function of in Fig.3 (red color). After transformation , the dependence of electron’s energy on becomes ambiguous for the first subband when is small. It consists of two curves: the closed loop surrounded by a parabolic branch. The horizontal arrows in Fig.3 indicates directions of the wave vector’s transformations. First, let us notice that for the closed loop, transforms to and due to symmetry transforms to . Second remark concerns the scalability of . In Fig.3(a) we see that value of kinetic wave vector is compressed for the lower subband and expanded for the second one with respect to canonical wave vector value. And last, the two lowest subbands become degenerated for much larger Fermi energies what is showed in Fig.3(b). It means that electrons in both subbands move with the same group velocity . This fact can be easily explained if energy dispersion relation given by Eq.17 will be expressed as function of instead of k. For this purpose, we first calculate an expectation value of kinetic wave vector for p-th subband:

(18) |

where coefficents and are components of two-element effective Hamiltonian () given by Eq.13. Using the components of two orthogonal eigenvectors, after some algebra we obtain a formula for the kinetic wave vectors:

(19) |

where stands for first (lower) subband while for the second (upper) one. Assuming very large k value we may expand expression with square root leaving only the first term and then rearrange equation to get k: . Substitution this approximate expression for k into Eq.14 gives:

(20) |

Now, it is easily to notice that if the kinetic wave vectors for the first and second subbands have the same value then these subbands are degenerated [Fig.3(b)].

### iii.2 Inter-layer subbands mixing for a non-symmetric confining potential

In this section we consider an effect of pure inter-layer subbands mixing () on the energy dispersion relation for and on the conductance of a bi-layer nanowire. In Fig.4 we have plotted the low energy spectra calculated within our model for and . If there is no in-plane magnetic field [first row in Fig.4], the vertical eigenmodes are not mixed and the symmetry of energy dispersion i.e. is kept independently of value. For the energy branches of subsequent subbands are only shifted upwards on energy scale in comparison to the case with . When value of is increased to [second row in Fig.4], then for the negative energy dispersion relation appears () in the three lowest energy subbands which correspond to the ground, the first and the second excited states for y direction. However, the difference between the maximum and minimum of energy for each of these subbands is very small since it equals only . If potential becomes slightly nonsymmetric for , it destroys the symmetry of subbands too. In such case, the upper quantum well is wider than the lower one what leads to a larger energy separation between the basis states and which for is . Although, this growth is not large, it is sufficient to suppress the negative energy dispersion relation for the moderate value of an in-plane magnetic field () since it is too small to overcome defined in Eq.16. In this case, the electrons which have and are localized in the lower layer due to an action of the magnetic force, have higher energies than those with [see localization of electron densities in Fig.1(b)] which move in a wider upper well. For this reason, the right parts of the energy spectra in Figs.4(d) and (f) are shifted upwards with respect to their left parts. Stronger magnetic field () enhances the negative energy dispersion in relation to case. In Fig.4(e) we may notice that the difference between maximum and minimum of energy equals 5 meV for three lowest subbands and is much larger when compared to fraction of meV we have got for . Mixing of vertical subbands is now so strong that the negative dispersion relation of energy is reconstructed also for [cf. Figs.4(d) and (f)].

Since mixing of the vertical modes and an appearance of additional lateral energy minima in magnetosubbands depend on , we have repeated calculations for its smaller value. The left column in Fig.5 display the energy spectra for and . Again, in the absence of an in-plane magnetic field component (), subbands have parabolic shape but with lower energy spacings between neighbours. Now however, in contrary to the previous case, a moderate in-plane magnetic field () effectively mixes the vertical eigenmodes what leads to a formation of two equally deep minima for and one deep (the left one) and one shallow (the right one) for . Note also that the subbands with negative energy dispersion relation which are lying higher on energy scale, cross with subbands of parabolic shape. These parabolic branches are formed when an electron occupies first excited state in vertical direction [in Eq.3 value of m is then reduced to ]. These subbands are not mixed by magnetic field and have therefore different parity than those which were already mixed. These energy crossings survive also for [cf. Figs.5(d) and 5(f)]. Let us notice here, that such crossings are absent in the low energy spectra presented in Fig.4(c) due to the larger value of as well as due to a limited range of the energy scale in Figs. 4(e) and 4(f). For stronger magnetic field () the energy minima for become more than five times deeper than for moderate field and almost two times deeper than those for [cf. Figs. 5(e) and 4(e)].

In Fig.6 we have plotted the conductance of a bi-layer wire for spin-up electrons and that of a single-layer wire for comparison. Conductance has been calculated for temperature by counting the crossings between the horizontal line at Fermi energy level with those parts of subbands which have . Results obtained for a single-layer wire [] show standard, well-known step-like raising function (grey colour) with hights of steps equal to the conductance unit [the spin degeneracy is lifted]. On the other hand, the conductance for a bi-layer wire [, red colour stands for while the black for ] still exhibits a step-like character but now two additional features appear: i) the heights of the steps may be equal to or and, ii) the value of conductance drops by when the Fermi energy exceeds the height of the central maximum in a particular subband. Due to a stronger coupling of the vertical modes for , the rising and the falling steps become better separated than for . This is easily noticeable if we compare e.g. the changes in conductance for in Figs.6 (a) and (b). When a confining potential looses spatial symmetry (), the conductance of bi-layer nanosystem is generally lower than for and approaches in some points a lower limit established by the value of conductance of a single-layer wire. In both cases, even though the conductance is significantly larger for bi-layer wire than for a single-layer case, their ratio very rarely reaches its upper limit which equals 2.

### iii.3 Mixing of the vertical and the transverse modes in tilted magnetic field

If we account the vertical component of magnetic field () in our considerations, then the term appearing in the diagonal part of effective Hamiltonian [Eq.6] becomes responsible for mixing of the transverse eigenstates while the transverse and vertical eigenmodes can be simultaneously mixed by matrix elements defined in Eq.10 if . The latter term contribute to the off-diagonal elements in Eq.11 for due to non-zero value of or simultaneously to the diagonal and off-diagonal submatrices in Eq.11 for .

Magnetosubbands for and are plotted in Figs.7(a-c) for and in Fig.7(d) for . For , subbands have parabolic-like shapes and do not cross each other. Since , all the off-diagonal elements [see Eqs.8-10] in effective Hamiltonian [Eq.11] vanish and only the diagonal elements given by Eq.7 survive but they cannot mix the vertical modes. Therefore, the subbands corresponding to vertical excitation of an electron [dashed lines in Fig.7(a)] are simply the replicas of those subbands in which the electron occupies the ground state in vertical direction. They are only shifted upwards on energy scale by . The wave functions of these subbands have different symmetries with respect to reflection and for this reason may mutually intersect in higher part of energy spectrum for [not shown in Fig.7(a)]. If both and have non-zero values, the off-diagonal elements in effective Hamiltonian which are defined by Eq.10 may mix the vertical and transverse modes since the products of terms and have non-zero values. In Fig.7(b) we see that for there are two deep lateral minima in the lowest subband similarily as in the case for and [cf. Figs. 4(c) and 5(c)]. However, for , due to mixing of the transverse modes, the subbands lying higher on energy scale are not the replicas of the lowest one any longer. In vicinty of , crossings are replaced by avoided crossings at positions where they maximally mix with their neighbours what in consequence opens additional pseudogaps in the energy spectrum. The widths of antycrossings are dependent on strength. If value of is increased, then the spatial localizations of subbands in the lower and upper layers are enhanced. It diminishes the overlaps between transverse modes belonging to different layers what in consequence decreases the widths of anticrossings [cf. Figs. 7(b) and 7(c)]. Generally, the pattern of these avoiding crossings does not change much when th upper and lower quantum wells are slightly different [] besides an accuracy of additional positive slope in energy [cf. Figs. 7(c) and 7(d)]. Even though, for the vertical modes shall be effectively mixed, the two deep lateral minima are present only in two lowest subbands. In upper subbands they are replaced by bending points due to a comparable mixing of the vertical and transverse modes. However, when the transverse component of magnetic field becomes two times stronger (), the coupling of vertical modes dominates in the system and these minima are again visible in the lowest subbands [cf. Figs. 7(c) and 7(d)]. In addition, the minima in question are localized at points that is exactly as in the lowest subband [Fig.7(c)].

An appearance of additional pseudogaps in energy spectrum significantly modifies the conductane of a bi-layer wire. An example dependence of wire’s conductance on Fermi energy is shown in Fig.8 for and . For symmetric confinement () the conductance raises by when the energy exceeds the bottoms of the first three subbands [i.e. for , and in Fig.8] which are determined by their lateral energy minima and it drops by for an energy exceeding next the central maximum in particular subband. When the energy is increased and passes through the pseudogap and then through the local energy minima localized in vicinity of , the value of conductance again jumps. Now however, the heights of G steps depend on the number of local minima belonging to a particular subband which equals the subband’s index plus one. For example, in Fig.8 we see that the conductance grows by when energy increases from to and even by when it is changed between and . On the other hand, it may also significantly drop if an electron’s energy is shifted through the set of local maxima localized in proximity of because it enters a pseudogap region then. We notice such a large drop, being equal to , when Fermi energy is increased from to . Generally, the conductance steps for are lower than for mainly due to a considerable shift of energy levels ladder in the narrower layer. This effectively suppresses the mixing of vertical modes because the value of significantly grows then. The conductance of a single layer wire constitutes a lower bound for conductance of bi-layer wire. Both can be equal only if Fermi energy is localized within a pseudogap. Such case is visible in Fig.8 for both values of and energies: and .

## Iv Spin polarization of conductance in bi-layer nanowire

In previous sections we have discussed the process of formation of pseudogaps in electron’s energy spectra and its influence on the wire’s conductance but have neglected the spin Zeeman effect contribution to energy. An interaction of electron’s spin with strong magnetic field, what is the case considered here, splits the spin-down and spin-up subbands by . For this reason, the conductance of bi-layer wire may be partly spin polarized Das et al. (2012) i.e. and this polarization shall be dependent not only on a number of active spin-up and spin-down subbands as it is in the case of a single-layer wire but also on that whether the Fermi energy is pinned within the pseudogap or not. In the later case one may expect a larger value of . In Fig.9 we have plotted the conductance, its derivative with respect to energy () and value of spin polarization of conductance () for bi-layer wire in function of and energy for and . These outcomes were obtained for the wire made of GaAs [first and second columns for and ] and of InGaAsRamos et al. (2007) [third column for and ].

In three figures 9(a,d,g) which show the conductance, we may notice two
characteristic regions lying above and under the anti-diagonal. In first region (above
the anti-diagonal), the changes of conductance values are frequent and can be increased as well as
decreased when electron’s energy grows. The second
characteristic region appears rather for strong magnetc field (under anti-diagonal). It has more
regular pattern resembling very much that of a single layer quantum wire as the value of
conductance increases by when subsequent subband becomes active.
Besides the conductance, also a transconductance is very often measured in experiment as it
directly reveals the dynamical properties of nanosystem being sensitive to the variations of
voltages applied to metallic gatesFischer
et al. (2006b) used e.g. to tune the Fermi energy in
the wire. Figures 9(b,e,h) show dependence of similar quantity i.e.
on and energy values. For GaAs, that has low value of g
factor, subsequent spin-up and spin-down subbands are gathered in pairs even for strong due
to a small energy splitting [Figs.9(b),(e)]. However,
because of small Zeeman energy splitting, conductance becomes partly polarized only for a very
narrow energy stripes what show Figs.9(c) and 9(f). Despite this fact, we
have found that even for GaAs wire, conductance can to a large extent be spin polarized. For
example, for
[Fig.9(c)] polarization may reach and
for pairs of parameters: , and ,
, respectively. Regions with similar values of polarization we have also found
in Fig.9(f) for .

If a bi-layer nanowire is made of InGaAs which has much larger g factor than
GaAs, then an energy splitting due to spin Zeeman effect becomes even comparible with an energy
difference between the bottoms of two neighbouring subbands [see the right part of
Fig.9(h)]. In such case, subbands with negative energy dispersion relation (blue curves)
are shifted significantly on energy scale even for moderate magnetic field [e.g. in Fig.9(h)].
For this reason, the regions of partly spin polarized conductance appearing for a wire
with low g factor in form of narrow stripes, now become much wider [see two distinct reddish stripes
appearing near the
central part of Fig.9(i), which mark and conductance polarization,
respectively].
This example show the advantage of bilayer quantum wire over e.g. Y-shaped
nanostructuresCummings et al. (2008); Wócik et al. (2015) in preparing partially spin polarized current. Bi-layer
nanowire enables one to get partially spin polarized current not only for the
lowest subband but also for those lying higher on energy scale giving thus larger conductances and
currents. Drawback of this solution is however that, it still requires a strong magnetic
field to work.

We have repeated calculations for spin polarization of conductance for InGaAs wire for temperature . Results are displayed in Fig.10. As expected, the temperature smearing of subbands makes the originally moderately spin polarized regions smaller and additionally, it lowers their polarization to about . This unvafourable effect can however be limited if semiconductor materials with much larger g factor like e.g. InSb Nadj-Perge et al. (2012) or InAsNadj-Perge et al. (2010) are to be used for the nanowire fabrication process.

## V Conclusions

We have theoretically investigated an effect of the magnetic field on the inter-layer and intra-layer subbands mixing for two vertically aligned nanowires with a rectangular-like external confining potential. For this purpose, a simple semi-analytical method was developed which has enabled us to calculate the energy subbands in dependence on electron’s wave vector value. It has been showed, that the transverse component of magnetic field, which is perpendicular to the wire’s axis but parallel to the layers, can effectively mix two lowest vertical eigenmodes what transforms the low energy subbands’ parabolas into slowly oscillating curves with two deep lateral energy minima. If besides the transverse component of magnetic field, the vertical one is also taken into account, then both the vertical and transverse modes are mixed simultaneously, crossings between subbands are replaced by avoided crossings what lifts the degeneracy between subbands in vicinity of and additional small pseudogaps appear in the energy spectrum. A qualitatively similar behaviour of magnetosubbands were predicted for two laterally aligned wires by Shi and Gu.Shi and Gu (1997) They showed that only one component of magnetic field, namely the perpendicular one, is needed for an effective mixing of all pairs of magnetosubbands which were originally localized in the same and in different wires. Consequently, such unrestricted hybridization of magnetosubbands gives then simultaneously both types of oscillations in energy spectrum E(k) that is, two deep lateral minima and the small-amplitude oscillations near . In contrary to a nanosystem with two laterally aligned wires, an interlayer subbands hybridization in the vertically aligned bi-layer wire depends on the transverse component of magnetic field only while the perpendicular one is responsible for the intralayer modes mixing.

Irrespective of the coupling direction, if these small-amplitude oscillations appear in energy dispersion relation E(k) then the conductance of a bi-layer wire may jump as well as drop by a few conductance quanta provided that the confining potential has low number of defects. When Fermi energy is shifted through these oscillations, value of conductance first jumps and then falls, within a thin energy region, even by several units of conductance. Results of our simulations show that the conductance of bi-layer quantum wire can be spin-polarized up to at zero temperature. Magnitude of spin polarization can be tuned by changing the strength of magnetic field and the value of Fermi energy.

Although there are no direct experimental results confirming our predictions, a number of experiments were performed for similar bi-layer nanosystems. Thomas et al.Thomas et al. (1999) have measured the conductance of two vertically coupled wires in dependence of the top gate voltages for parallel and perpendicular magnetic fields. Also in work of Fischer et al.Fischer et al. (2005), the conductance of similar bi-layer nanowire with nonsymmetric vertical confining potential was experimentally investigated. In both experiments however, the lateral confinement was smooth whereas our predictions concern the nanowires with rectangular-like lateral confinement and therefore the outcomes of calculations and the experimental data can not be directly compared.

## Acknowledgements

The work was financed by Polish Ministry of Science and Higher Education (MNiSW)

## References

### References

- C. C. Eugster and J. A. del Alamo, Phys. Rev. Lett. 67, 3586 (1991).
- O. Bierwagen, C. Walther, W. T. Masselink, and K.-J. Friedland, Phys. Rev. B 67, 195331 (2003).
- S. K. Lyo, J. Phys.: Condens. Matter 8, L703 (1996).
- L. G. Mourokh, A. Y. Smirnov, and S. F. Fischer, Appl. Phys. Lett. 90, 132108 (2007).
- S. F. Fischer, G. Apetrii, U. Kunze, D. Schuh, and G. Abstreiter, Phys. Rev. B 74, 115324 (2006a).
- D. Huang, S. K. Lyo, K. J. Thomas, and M. Pepper, Phys. Rev. B 77, 085320 (2008).
- S. F. Fischer, G. Apetrii, U. Kunze, D. Schuh, and G. Abstreiter, Phys. Rev. B 71, 195330 (2005).
- K. J. Thomas, J. T. Nicholls, M. Y. Simmons, W. R. Tribe, A. G. Davies, and M. Pepper, Phys. Rev. B 59, 12252 (1999).
- J.-R. Shi and B.-Y. Gu, Phys. Rev. B 55, 9941 (1997).
- E. Bielejec, J. A. Seamons, J. L. Reno, and M. P. Lilly, App. Phys. Lett. 86, 083101 (2005).
- S. Kumar, K. J. Thomas, L. W. Smith, M. Pepper, G. L. Creeth, I. Farrer, D. Ritchie, G. Jones, and J. Griffiths, Phys. Rev. B 90, 201304 (2014).
- S. Buchholz, S. Fischer, U. Kunze, D. Schuh, and G. Abstreiter, Physica E 40, 1448 (2008).
- W. K. Hew, K. J. Thomas, M. Pepper, I. Farrer, D. Anderson, G. A. C. Jones, and D. A. Ritchie, Phys. Rev. Lett. 102, 056804 (2009).
- Fischer S. F., Apetrii G., Kunze U., Schuh D., and Abstreiter G., Nat. Phys. 2, 91 (2006).
- A. Tarasov, S. Hugger, H. Xu, M. Cerchez, T. Heinzel, I. V. Zozoulenko, U. Gasser-Szerer, D. Reuter, and A. D. Wieck, Phys. Rev. Lett. 104, 186801 (2010).
- M. Huber, M. Grayson, M. Rother, W. Biberacher, W. Wegscheider, and G. Abstreiter, Phys. Rev. Lett. 94, 016805 (2005).
- N. R. Abdullah, C.-S. Tang, and V. Gudmundsson, Phys. Rev. B 82, 195325 (2010).
- S. Ihnatsenka and I. V. Zozoulenko, Phys. Rev. B 73, 075331 (2006).
- S. Schnez, C. Rössler, T. Ihn, K. Ensslin, C. Reichl, and W. Wegscheider, Phys. Rev. B 84, 195322 (2011).
- B. Hackens, F. Delfosse, S. Faniel, C. Gustin, H. Boutry, X. Wallart, S. Bollaert, A. Cappy, and V. Bayot, Phys. Rev. B 66, 241305 (2002).
- P. P. Das, N. K. Bhandari, J. Wan, J. Charles, M. Cahay, K. B. Chetry, R. S. Newrock, and S. T. Herbert, Nanotechnology 23, 215201 (2012).
- A. V. Ramos, M.-J. Guittet, J.-B. Moussy, R. Mattana, C. Deranlot, F. Petroff, and C. Gatel, App. Phys. Lett. 91, 122107 (2007).
- S. F. Fischer, G. Apetrii, U. Kunze, D. Schuh, and G. Abstreiter, Phys. Rev. B 74, 115324 (2006b).
- A. W. Cummings, R. Akis, D. K. Ferry, J. Jacob, T. Matsuyama, U. Merkt, and G. Meier, J. Appl. Phys. 104, 066106 (2008).
- P. Wócik, J. Adamowski, M. Wołoszyn, and B. J. Spisak, J. Appl. Phys. 118, 014302 (2015).
- S. Nadj-Perge, V. S. Pribiag, J. W. G. van den Berg, K. Zuo, S. R. Plissard, E. P. A. M. Bakkers, S. M. Frolov, and L. P. Kouwenhoven, Phys. Rev. Lett. 108, 166801 (2012).
- S. Nadj-Perge, S. M. Frolov, J. W. W. van Tilburg, J. Danon, Y. V. Nazarov, R. Algra, E. P. A. M. Bakkers, and L. P. Kouwenhoven, Phys. Rev. B 81, 201305 (2010).