Influence of Device Geometry and Imperfections on the Interpretation of Transverse Magnetic Focusing Experiments

Spatially separating electrons of different spins and efficiently generating spin currents are crucial steps towards building practical spintronics devices. Transverse magnetic focusing is a potential technique to accomplish both those tasks. In a material where there is significant Rashba spin–orbit interaction, electrons of different spins will traverse different paths in the presence of an external magnetic field. Experiments have demonstrated the viability of this technique by measuring conductance spectra that indicate the separation of spin-up and spin-down electrons. However the effect that the geometry of the leads has on these measurements is not well understood. By simulating an InGaAs-based transverse magnetic focusing device, we show that the resolution of features in the conductance spectra is affected by the shape, separation and width of the leads. Furthermore, the number of subbands occupied by the electrons in the leads affects the ratio between the amplitudes of the spin-split peaks in the spectra. We simulated devices with random onsite potentials and observed that transverse magnetic focusing devices are sensitive to disorder. Ultimately we show that careful choice and characterisation of device geometry are crucial for correctly interpreting the results of transverse magnetic focusing experiments. Supplementary Information The online version contains supplementary material available at 10.1186/s11671-022-03671-x.


Introduction
Unlike conventional electronics, spintronics manipulates the spin degree of freedom of an electron to transfer and store information [1]. One of the challenges in spintronics is to generate a spin current efficiently [2], and another is to spatially separate spin currents [3]. Transverse magnetic focusing (TMF) offers a solution to both these challenges by exploiting the spin-momentum coupling induced by the Rashba spin-orbit interaction [4]. Electrons injected into a two-dimensional electron gas (2DEG) with a significant Rashba spin-orbit coupling are separated by their spins and travel in different cyclotron paths as shown in Fig. 1a when a perpendicular magnetic field is applied to the device. Measuring the conductance between the points of electron injection and collection, while varying the strength of this magnetic field, gives a spectrum like the one shown in Fig. 1b, with peaks occurring at magnetic field strengths that focus an electron at the point of collection.
TMF in a metal was first demonstrated by Tsoi in 1974 using a Bi single crystal [5]. TMF in a 2DEG was shown by van Houten et al. in 1988 using a GaAs/AlGaAs heterostructure at temperatures ranging from 0.05 K to 7 K [6]. In 1992, Heremans et al. observed ballistic hole transport by performing a TMF experiment in a twodimensional hole gas (2DHG) using a GaAs/AlGaAs heterostructure [7,8]. More recently, TMF has been shown to be a viable spin filter in a wide range of experiments [3,[9][10][11][12][13][14][15][16][17][18]. There are also indications that TMF need not Open Access *Correspondence: yik.kheng.lee@student.rmit.edu.au 2 ARC Centre of Excellence in Future Low-Energy Electronics Technologies, RMIT University, Melbourne, Australia Full list of author information is available at the end of the article be limited to semiconductor heterostructures, as Taychatanapat et al. [19] demonstrated it can be performed in graphene devices up to a temperature of 300 K.
In a TMF experiment, spatial separation of the spin currents is indicated by the presence of a pair of spin-split peaks in the conductance spectrum that each correspond to a different spin, such as the peaks at B z = 0.061 T and 0.089 T in Fig. 1b. In materials without Rashba spin-orbit coupling, the spin currents would not spatially separate, resulting in a single higher peak, similar to the peak at B z = 0.151 T in Fig. 1b. The spin-split peaks do not have equal amplitudes, which has previously been attributed to spin polarisation in the injected electrons [11,20]. The wave-like nature of electrons suggests that there should be wave interference, which can be observed as interference fringes in the spectrum [21,22]. Although transport in these devices is considered ballistic, there are many factors that could affect measurements in experiments, such as temperature [9,13,15] and disorder. The effects of physical device parameters such as the shape and width of the leads used remain unclear. Square-shaped [3], trumpet-shaped [12] and funnel-shaped [9] leads have all been used in experiments.
In this paper, we focus on computing TMF measurements in an InGaAs-based 2DEG with Rashba spin-orbit coupling with realistic device geometries. We simulate TMF using wave function matching as implemented through the program KWANT [23] using a finite difference method with a square grid discretisation scheme. A brief introduction of TMF is given in the "Semi-classical Model" section. We then present our numerical approach in the "Numerical Model" section. The effects that the shape, separation distance and width of the leads have on the TMF response are presented in the "Influence of device geometry" section. The effects of varying the number of subbands occupied by the injected electrons are presented in the "Subbands" section. We simulate the leads as quantum point contacts (QPCs) in the "Simulation of realistic QPCs" section and show that effects seen in previous sections are linked in this more realistic model. We investigate the influence that disorder in the surface potential of a 2DEG has on the TMF results in the "Disorder" section and conclude in the "Conclusions" section.

Semi-classical Model
A typical TMF experiment consists of a two-dimensional focusing area connected to two leads which serve as injector and collector for electrons (see Fig. 1a). An outof-plane magnetic field B z is applied such that an electron injected into the system experiences a Lorentz force which causes it to move in a cyclotron orbit with a Larmor radius R L given by [6]: where is the reduced Planck constant, k is the magnitude of the wavevector of the electron, and e is the charge of an electron.
At particular values of B z , where the distance between the leads coincides with integer multiples of 2R L , the probability of transmission for an electron between these two leads is nonzero as long as the reflections from the wall are specular. This is observed as sharp peaks in the conductance (as shown in Fig. 1b). Experiments usually measure current and resistance, which are proportional to the conductance spectra we calculate in this paper.
(1) R L = k |e|B z Fig. 1 a Schematic of a typical transverse magnetic focusing device. Electrons injected at lead L 1 are focused and collected at lead L 2 by applying an out-of-plane magnetic field B z . If the Rashba effect is present, spin-up (red) and spin-down (blue) electrons have different trajectories due to spin-orbit interaction. Distance between the leads (D) is defined from the middle of each lead. b Example conductance spectrum for a transverse magnetic focusing device showing conductance as a function of the out-of-plane magnetic field calculated using the method described in Sect. Numerical Model. The solid line (purple) shows the total conductance, the dashed (red) and dotted (blue) lines show the spin projected conductance for spin-up and spin-down in x, respectively When a strong Rashba spin-orbit interaction is present in the system, spin-up and spin-down electrons entering the 2DEG with the same energy have different crystal momenta. Consequently, spin-up and spin-down electrons have orbits with different Larmor radii, which results in the spatial separation of spin currents as indicated in Fig. 1a and the presence of spin-split peaks in the conductance spectrum as seen in Fig. 1b. The Rashba interaction is described by the Hamiltonian [4]: where α is the Rashba coupling parameter, σ is the Pauli vector, k is the wavevector of the electron, and ẑ is the unit vector normal to the plane of the 2DEG.
The dispersion relation of a 2DEG with Rashba interaction is then given by [4]: where E ± are the energies of the spin-up(+)/down(−) electrons, k = |k| , and m * is the effective mass of the electron.
The magnetic field strengths, B + z (B − z ), required to focus spin-up(spin-down) electrons of the same Fermi energy, E F , into orbits of the same Larmor radius will be [11]: The difference between these two magnetic field strengths is the separation between the spin-up and spin-down peak pairs in the conductance spectrum. This semi-classical approximation ignores quantum effects such as the Berry curvature correction to the velocity of the electron [24]. However, for the parameter regime we consider, such corrections are several orders of magnitude smaller than the velocities. Thus we do not consider such corrections in this paper.
It is also worth noting that classically the Larmor radius (Eq. 1) is written in terms of the electron's tangential velocity. In a semi-classical treatment, one might therefore be tempted to substitute the group velocity ( v g = 1 ∂E ∂k ) for the tangential velocity, instead of substituting the crystal momentum ( k ) for its counterpart in the classical theory. There is a subtlety in the semi-classical treatment, however, where both the Larmor radius and the cyclotron frequency depend on the spin state of the electron via the Rashba coupling, and these effects cancel such that the electron's velocity is independent of spin [25]. In other words, the fact that the velocities are the same for the two spin states does not mean their Larmor radii are equal, as their cyclotron frequencies also differ. Our numerical approach, which is fully quantum mechanical, explicitly solves for the wave propagation of the electron via the system's Hamiltonian and therefore does not depend on the definition of the Larmor radius. Nevertheless we find our results correspond well to the radii predicted by equations 1 and 4.
We can now discuss the interpretation of a typical TMF response using this semi-classical model. The conductance between the leads varies with magnetic field strength, as illustrated in Fig. 1b. The first pair of peaks located at approximately 0.061 T and 0.089 T correspond to when the Larmor radii of the spin-down and spin-up electrons coincide with the distance between the leads, respectively, resulting in almost pure spin currents. The effect can be seen more clearly by plotting the probability density of conduction electrons injected into a TMF device. In Fig. 2a, b we see unfiltered spin-down and spin-up electrons spatially separate into distinct trajectories in the focusing area.
The peak at B z ≈ 0.151 T in Fig. 1b shows where the spin currents combine into a single current, giving a much higher total conductance. If we continue increasing the magnetic field strength, we eventually get another pair of peaks around 0.225 T that are less distinct compared to the first pair. Our numerical results show good agreement with the semi-classical approximation of the peak positions calculated using Eq. 4. More sophisticated semi-classical approximations can also allow the finer features in the conductance spectrum to be analytically calculated [22].
The recombination of the spatially separated spin currents upon reflection from a wall of the device is a known phenomenon [3,21], but deserves more discussion as the reason it happens may not immediately be clear. The probability density plot in Fig. 2c corresponds to the highest peak at approximately 0.151 T in Fig. 1b. We see that the spatially separated trajectories of the spin-up and spin-down electrons overlap at the collector after reflecting off the bottom edge of the device. This suggests the crystal momenta of the electrons change after reflection from a hard wall boundary as the Larmor radius of each trajectory has changed.
The picture of the spin-up and spin-down electrons moving in their own trajectories as shown in Fig. 1a is greatly simplified and commonly used in the TMF literature. In this picture, spin-up and down are defined to be parallel or anti-parallel to the Rashba effective magnetic field, which is always perpendicular to the direction of the instantaneous velocity of the electron. In contrast, if a stationary frame is chosen for the definition of spinup and spin-down (e.g. Cartesian coordinates in the lab frame), it becomes more obvious that the spin of the electron precesses as it travels through the focusing area [25]. Figure 3a illustrates this; the spin of the electron is always aligned tangentially to the Fermi surface in k-space, and at points where the electron is moving only in the x or y direction, its spin is aligned along the y or x direction, respectively.
If we follow the path of an electron that has spin-down in x when injected into the focusing area as shown by the solid line in Fig. 3b, we see how the spin precession leads to the spin currents recombining after reflecting off a hard wall boundary in the device. At the injection point, labelled 1 , the electron is moving in the positive y direction, which corresponds to point 1 on the outer Fermi surface in k-space. It has a larger crystal momentum than an electron with spin-up in x (which starts on the inner Fermi surface at 4 ) and therefore a larger Larmor radius (Eq. 1). Here we see that the group velocity is proportional to ∂E ∂k , and the gradients of the inner and outer Fermi surfaces are equal in magnitude, so the spinup and spin-down electrons have equal group velocities. Although the the group velocities have the same signs as the corresponding crystal momenta at this particular energy, this is not necessarily always the case, e.g. when the Fermi surfaces are no longer concentric near the bottom of the band.
As the electron moves through the focusing area, it flips to spin-up in x after point 2 . When it reaches the boundary of the device at point 3 , the electron is moving in the negative y direction before it is reflected. After reflection, spin is conserved but velocity changes sign, i.e. it switches from point 3 on the outer Fermi surface to point 4 on the inner Fermi surface, similar to the effect observed by Chen et al. [26]. The electron now has a lower crystal momentum and a shorter Larmor radius as it continues moving through point 5 , where it flips back to spin-down in S x , and exits the device at point 6 . An electron which is spin-up in S x at injection would go through a similar process, but starts with a lower k at injection and changes to a higher k after reflection. Thus the separate spin currents recombine after reflecting off the boundary of the system once (as seen in Fig. 2c), resulting in a sharp peak in the conductance spectrum.

Numerical Model
For all calculations in this paper, we simulated a TMF device with a geometry as shown in Fig. 1a. The injection and collection leads of the device were separated by 1500 nm as measured from the centre of each lead. The left, right and top edges of the device are simulated as open boundaries by attaching semi-infinite leads to them. The bottom edge is a hard wall boundary except for the injection and collection leads. An additional 312.5 nm The colour in the plots indicate the direction of spin as measured along the x axis. The ratio of red to blue in the RGB colour values at each point is set to be proportional to the expected spin states at that point. Pure red indicates a pure spin-up state, and pure blue indicates a pure spin-down state, while purple indicates a mixed state. The transparency values are on a normalised tanh scale, given by tanh(2d i /d max )/ tanh(2) , where d i is the density at each point and d max is the maximum density in each plot. Points with the maximum density in each plot will be completely opaque, while points with zero density will be completely transparent was added between the side edges and the leads to serve as a buffer. The size of this buffer was chosen to be large enough that it did not significantly affect the results. The height of the device (length in y direction) was chosen to be the equal to the length of the separation between the contacts, 1500 nm, so that the focusing area was large enough to accommodate the cyclotron trajectory of an electron focused from the injection lead to the collection lead.
The device was described as a 2DEG confined to the xy plane, with an external magnetic field B z applied along the z axis and Rashba spin-orbit coupling with the following Hamiltonian [21]: where g is the effective g-factor, µ B is the Bohr magneton, and σ x,y,z are the Pauli spin matrices. To approximate the effect of the magnetic field on the electron we used the Peierls substitution [27]: where A x is the x-component of the magnetic vector potential. The y-components were similarly substituted. The Landau gauge in ŷ , A = −xB zŷ , was used everywhere except in the semi-infinite leads on the left and right edges and a thin strip of the focusing area next to those leads, where the Landau gauge in x , A = −yB zx , was used instead. We used the effective mass m * = 0.04m e where m e is the free electron mass, and an effective g-factor g = 9 for InGaAs [28]. For the Rashba coupling parameter, we used the value for an InGaAs-based 2DEG, [3]. We note that other spin-orbit interactions such as the Dresselhaus effect are present in InGaAs [29], however they are not as strong as the Rashba effect, therefore we only included the Rashba effect in this paper for simplicity. The external magnetic field and Rashba spin-orbit coupling were applied everywhere in the system, including within the leads. Note that while we refer to the injection and collection points as leads, they are still part of the 2DEG, and typically correspond to QPCs in experiments. Our calculations used a finite difference method to discretise the Hamiltonian given in Eq. 5, resulting in a real space grid as shown in Fig. 4. Each vertex in the grid represents a discretisation site that has an associated on-site potential, and the lines represent hoppings that couple neighbouring sites. Neighbouring sites in the grid are separated by a length of 3.125 nm. This grid spacing was chosen by balancing computational efficiency and numerical accuracy, while ensuring that results converged. Energies in our calculations are given in terms of a "hopping energy" t, defined as: where a is the grid spacing. The Peierls substitution modifies the hopping terms by a Peierls phase [27], e.g. if the hopping term between sites r 1 and r 2 is t, then: We calculated our results using the software package KWANT, which uses a wave function matching method that is mathematically equivalent to the non-equilibrium Green's function method [23,30]. This model assumes ballistic transport (the electrons have mean free paths longer than the device length) and the zero temperature limit (no thermal broadening). It is a fully quantum description that captures the wave-like nature of electrons as can be seen in the density plots in Fig. 2 where the trajectories are clearly not classical cyclotron trajectories. We did not apply any bias between the leads in all of our calculations, so conductance was calculated using the Landauer-Büttiker formalism with zero-bias voltage, giving [31]: where h is Planck's constant, T(E) is the transmission probability between the injector and collector, and f 0 (E) is the Fermi-Dirac distribution function. In the zero temperature limit, the derivative of the Fermi-Dirac distribution function is a Dirac delta function, thus the integral in Eq. 9 is only evaluated at the exact energy chosen. The conductance is then directly proportional to the transmission probability, with a proportionality constant given by the conductance quantum, 2e 2 /h. The spin-projected probability densities in Fig. 2 are calculated by: where for each site i, ρ i,± is the spatially-resolved probability density for spin-up/down in x, ψ i is the wavefunction element corresponding to that site, and χ ± are the eigenvectors for the x Pauli matrix.

Influence of Device Geometry
An important consideration in TMF experiments is the influence of device geometry on the resulting signal. Whilst the leads in actual devices are QPCs formed by smoothly varying potentials, here we model them as being confined by simple hard walls such that there are well-defined widths. We present a more realistic model of the leads in the "Simulation of realistic QPCs" section. We define the width of a lead as the width of its narrowest section, and its shape by the position of that section relative to the contact point between the lead and the focusing area. We find three different shapes used in TMF experiments in the literature. First, there is the funnel-shaped lead (Fig. 4a) which has its narrowest section at the contact point and becomes wider towards the electron source [6,9,16]. Second, the square-shaped lead (Fig. 4b) maintains the same width throughout its entire length [3]. Third, the trumpet-shaped lead (Fig. 4c) is at its widest at the contact point and becomes narrower towards the electron source [12]. Figure 4a-c shows the grids we used to simulate each of these shapes. Finally, we also show a grid in Fig. 4d that represents an idealised point (one-dimensional) lead, the end of which is connected to a number of sites in the focusing area, giving the lead an effective width. For example, a 50 nm wide point lead refers to a 1D lead that is connected to 17 sites in the channel, which for a grid spacing of 3.125 nm is equivalent to a width of 50 nm (see Fig. 4d).
This idealised point lead will only have one transmission mode [21]. It is a simple theoretical model that ignores the effects of wave interference in the lead itself. Similar grids for the point lead [21] and square-shaped lead [32] have been used in the literature for numerical simulations of TMF experiments.
The shape of the lead influences the impedance matching between the lead and the focusing area, and thus we expect that there is an effect on the measured conductance spectra. The separation distance between the leads has also been shown to have an effect on the conductance spectrum [14]. We calculated conductance as a function of the out-of-plane magnetic field for TMF devices using injector and collector leads of each of the aforementioned shapes with widths of 50 nm. The separation distance (10)  between the leads was set to 500 nm, 700 nm, 1500 nm and 2500 nm. Here a Fermi energy of E F = 7 meV was used. Figure 5 shows the resulting conductance spectra with the conductance normalised to the maximum value in each data series. Note that each column has a different scale for the x-axis since the peak positions shift towards lower magnetic field strengths as the lead separation is increased, which agrees with Eq. 1. We see that for each separation distance, the spectra for all lead shapes have features similar to the spectrum for the idealised lead (point).
The main peaks are easily resolvable in all the spectra, but finer structures such as the interference fringes are only clearly observed at the larger separations of 1500 nm and 2500 nm, when the peaks become sharper. The shapes of the spectra change as the separation distance is changed, but stay relatively similar when the lead shape is changed for the same separation. Interestingly, at a separation distance of 750 nm, a splitting appears in the second focusing peak, which should not have spin-splitting as discussed in Sec. 2.1. This splitting appears to be dependent on the separation distance as it is not noticeable at other separation distances.
The spectra for the trumpet shaped lead appear to have less fluctuations compared to the square and funnel shaped leads for all separation distances. While the interference fringes before the first and second peaks are better resolved in the funnel and square shaped leads, we note that the magnitude of the fringes are comparable to the fluctuations seen in the peaks at higher magnetic field strengths. Thus it would not be easy to differentiate between path interference fringes [22] and fluctuations from coincidental interference effects. We chose to use the trumpet-shaped lead for subsequent calculations in this paper as it appeared to consistently give the best combination of low noise and clear resolution of features in the conductance spectrum, while still being an experimentally realistic geometry.
Another geometric aspect that affects the measured conductance spectra is lead width. Intuitively we would expect that the wider the leads, the wider the peaks would be in the spectra, since a wider injector would create a larger wavefront while a wider collector would collect more of that wavefront at each magnetic field strength. Figure 6 shows the transmission spectra around the first peak pair in these spectra at B z ≈ 0.061 T and B z ≈ 0.089 T for TMF devices with the injector and collector widths set to combinations of 50 nm, 100 nm, 150 nm, 200 nm, and 250 nm. Here a Fermi energy of 7 meV was used. We see that resolving the spin-split peaks depends on both the injector and collector widths. The resolvability of the peaks is reduced as we move across each row and down each column of Fig. 6, with the clearest peaks observed when both leads are at their minimum widths (here 50 nm).
If we look at the conductance spectra of a TMF device using idealised leads of the same widths, as shown in Fig. 7, we can see a similar degradation in the resolvability of the peaks as the widths of the leads increase. The individual peaks broaden as the lead widens, and we observe that the interference fringes before the spin-split peaks are not observable at widths above 50 nm. In general, the spin-split peaks broaden as the leads widen and the interference fringes become less observable.

Subbands
Thus far we have not considered the number of subbands occupied by the electrons when we choose the Fermi energy in our calculations. In experiments this can be controlled by changing the gate voltage at the leads to have conductance values that correspond to the number of subbands occupied [28]. Adjusting the gate voltage in experiments effectively shifts the band structure within the leads, whilst the Fermi energy stays constant.
In our calculations, we chose to shift the Fermi energy instead for numerical simplicity. While this will shift the positions of the peaks in the conductance spectra for different Fermi energies, the effect this has on the overall features of the spectra (e.g. peak shapes) is equivalent to shifting the band structure. We will present a more realistic model of the leads as QPCs controlled by gate voltages in Sec. 3.3. We examine in Fig. 8 how the features of the peaks at the first focusing magnetic field strength change as the the Fermi energy shifts through the lowest four odd numbers of subbands. The TMF device simulated has trumpet-shaped leads with widths of 100 nm separated by 1500 nm. To interpret these results, Fig. 9 shows the 1D band structure calculated for a nanowire with a width of 100 nm and the dotted lines indicate the corresponding Fermi energies shown in Fig. 8.
In Fig. 8, we see that the amplitude of the peaks increases as the Fermi energy increases, as more conducting modes are available when more subbands are occupied. The ratio of the first peak amplitude to the second peak amplitude changes from low:high to high:low as more subbands are occupied, but the peak shapes and amplitudes do not change significantly when the Fermi energy changes within the same subband. In the spectra that correspond to one occupied subband (blue), the first spin-split peak is greatly suppressed. This strong spin polarisation is explained by adiabatic transitions Fig. 7 Conductance spectra calculated with the same parameters as in Fig. 6 using idealised point leads (solid lines) and trumpet leads (dashed lines) with both injector and collector widths as labelled and set 1500 nm apart. Conductance values are normalised to the maximum value in each spectrum and offset for clarity Fig. 8 Conductance as a function of the out-of-plane magnetic field for a transverse magnetic focusing device with 100 nm wide trumpet shaped leads set 1500 nm apart, calculated with Fermi energies set between 8 and 75 meV. The Fermi energy crosses one (1 meV and 2 meV), three (8 meV and 11 meV), five (21 meV and 26 meV), and seven (40 meV and 48 meV) subbands as indicated, corresponding to the band structure shown in Fig. 9. The shape of the peaks changes as more subbands are occupied, but stays the same for energies within the same subband  Fig. 8 between anti-crossed subbands as the lead widens [33]. However, note that we do not observe the suppression of the first spin-split peak in the spectra along the first row and column of Fig. 6, even though there is only one occupied subband at 7 meV for a 50 nm wide nanowire. For a 50 nm wide nanowire, the anti-crossing between the first and second subbands occurs at an energy close to the minimum of the third subband, a much higher energy compared to the case for a 100 nm wide nanowire (see Fig. 9). This suggests the 50 nm wide lead does not widen enough for the electron to transition through the anti-crossing point, thus the first spin-split peak is not suppressed.
The band structure significantly affects the features of the TMF spectra, therefore it is vitally important to ensure that the number of subbands occupied and the spin compositions of the subbands are known when interpreting TMF spectra.

Simulation of Realistic QPCs
In experiments, the injector and collector leads in Fig. 1a are QPCs formed by applying voltage to surface gates, depleting electrons within the 2DEG under the gates [6]. Changing the gate voltages effectively changes the width of the QPC [6] and shifts the band structure within the QPC. We can approximate the potential field within the 2DEG applied by the surface gates using the method derived by Davies et al. [34]. Adding this potential field as a position dependent potential term to the Hamiltonian in Eq. 5 allows us to model the TMF device more realistically.
For our calculations we used the potential generated from surface gates with shapes illustrated by the orange lines in Fig. 10.
The injector QPC is formed by gates 3 and 4, while the collector QPC is formed by gates 1 and 2. The gates forming each QPC are separated by 300 nm. To avoid asymmetric biasing in the QPCs, we apply equal voltages on gates 1 and 2, as well as gates 3 and 4. We define the injector and collector gate voltages ( V g ) as the voltage applied to the gates forming the respective QPCs. The semi-infinite leads attached below each QPC have the same potential profile as the sites adjacent to the semiinfinite leads.
We calculated the transmission spectra for a TMF device with a lead separation of 1500 nm and a Fermi energy of 7 meV. The injector and collector gate voltages were set to combinations of 10 mV, 20 mV, 30 mV, 40 mV, and 50 mV (Fig. 11). The conductance values in Fig. 11 are normalised for better comparison of the peak features. However, we note that the conductance values for gate voltages of 50 mV (bottom row and right-most column in Fig. 11) are actually several orders of magnitude smaller than the conductance values in the other spectra. Also, the conductance values for gate voltages of 40 mV are close to one conductance quanta. This suggests that the QPC formed by gate voltages of 40 mV is only allowing transmission in the lowest subband, and at 50 mV, the only transmission occurring is from tunneling electrons. We see that similar to Fig. 6, the peaks broaden as the effective width of the QPCs  Conductance as a function of the out-of-plane magnetic field for transverse magnetic focusing devices with various combinations of gate voltages at the injector and collector. All leads are set 1500 nm apart, with the Fermi energy set to 7 meV. Conductance values are normalised to the maximum value in each spectrum increase (i.e. the voltage decreases). However, unlike Fig. 6, reducing the effective width of the QPC does not allow the interference fringes to be resolved even with the narrowest QPCs.
This suggests that to observe the interference fringes experimentally, the QPCs will need to have sharp hard wall-like potential profiles.
We also see that the first peak is greatly suppressed at gate voltages of 40 mV, similar to the single subband spectra in Fig. 8. Again, this spin polarisation can be explained by adiabatic transitions between anti-crossed subbands as the QPC narrows and widens [33]. Thus, we can see that the effect of changing the gate voltage is a combination of changing the width and Fermi energy in our previous simpler model. While the hard wall model of leads is less realistic, it allows us to explore the effects of width and Fermi energy independently.

Disorder
Devices used in experiments are usually not perfect, they exhibit defects, impurities, surface roughness and other imperfections. This causes electrons to scatter as they travel through the device. Here we simulate disorder using an Anderson localisation model [35], where we add a random onsite potential with values chosen from a uniform distribution within the range [−U /2, U /2] to each site within the focusing area, but not within the semiinfinite leads (the red sites in Fig. 4). U is defined as the disorder strength and can be related to the mean free path ( l mfp ) of the electrons using Fermi's golden rule [36] (Derivation shown in Additional file 1): where a is the length of the grid spacing, E F is the Fermi energy, and t is the hopping energy as defined in Eq. 7. For our calculations, we set a to be 3.125 nm and E F to be 7 meV. The mean free path lengths at varying disorder strengths calculated from Eq. 11 is shown in Fig. 12. The two inset probability density plots show that the spin-split trajectories become less well-defined as the mean free path length becomes shorter than the distance between the injector and collector leads. Figure 13 shows the conductance spectra calculated for a TMF device for various disorder strengths U. We performed conductance calculations for 201 different disorder realisations at each disorder strength. The mean of the conductance values at each magnetic field strength and disorder strength are shown as solid lines in Fig. 13, and the shaded areas show the variance of the conductance values for each disorder strength. We observe that (11) as disorder increases, the amplitude of the conductance peaks decrease as expected since scattering reduces the probability that the electron reaches the collector lead. The spin-split peaks at B z ≈ 0.061 T and B z ≈ 0.089 T are greatly suppressed and indistinguishable from signal noise at disorder strengths of 0.25t and above.
This agrees with our calculation that the mean free path of the electron at U = 0.25t is less than half the distance between the injector and collector leads (Fig. 12), thus the electron is likely to be scattered without reaching the collector.
We only simulated disorder within the weak localisation regime where the mean free path length is much greater than the Fermi wavelength. However, even within this regime, the characteristic features in the conductance spectra rapidly became indistinguishable as the devices were increasingly disordered. Therefore, this implies that to observe the finer structure in the spectra such as the interference fringes, the device used has to be relatively free from disorder.

Conclusions
We have shown the geometry of the device used in a TMF experiment plays a significant role in the measured characteristics of the conductance spectra. The shape of the leads used does not significantly affect the separation of the spin-split peaks, but does affect the level of noise or fluctuations in the spectra. The separation distance between the injector and collector leads can significantly Fig. 12 Mean free path length ( l mfp ) of an electron in a device with disorder strength U calculated using Eq. 11. The lattice constant (a) and Fermi energy ( E F ) were set to be 3.125 nm and 7 meV, respectively. The inset plots show the probability density for devices at a magnetic field strength of B z = 0.089 T, with particular realisations of disorder strengths 0.10t and 0.25t. The horizontal black dashed line at l mfp = 1500 nm indicates the distance D between the injector and collector leads in the device used in Fig. 13. The vertical dashed lines correspond to the disorder strengths used in Fig. 13 alter the features of the spectra. The width of the leads used significantly affects the resolution of the finer features in the TMF spectra, e.g. the interference fringes. Furthermore, the spin composition of the subbands and the number of subbands occupied both affect the ratio of the spin-split peak amplitudes. Thus, accurate knowledge of the device geometry and band structure is essential when interpreting TMF results.
As a result of these effects, the potentials forming QPCs in a TMF device need to be very well-defined in order to observe finer structures in the TMF spectra, e.g. the interference fringes. In addition, we show that the TMF device is sensitive to disorder and features in the measured conductance spectrum quickly become obfuscated as the device becomes disordered. The use of real-space numerical simulation of TMF experiments allows us to disentangle geometric effects and device imperfections from the spin-split structure and the interference fringes. This in turn allows for unambiguous interpretation of TMF experiments in the future.

Fig. 13
Conductance as a function of the out-of-plane magnetic field for a transverse magnetic focusing device with injector and collector QPCs set 1500 nm apart, and both gate voltages set to a 30 mV and b 40 mV. Fermi energies were set to 7 meV. Disorder was implemented in the systems by applying onsite potentials chosen at random from a uniform distribution within the range [−U/2, U/2] . Solid lines show the mean values for 201 realisations of that particular disorder strength, and the corresponding shaded areas show the peak-to-peak variance