## Manifestation of two-channel nonlocal spin transport in the shapes of the Hanle curves

R. C. Roundy , M. C. Prestgard , A. Tiwari , E. G. Mishchenko , and M. E. Raikh 1 2 2 1 1 1 Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112, USA 2 Department of Materials Science and Engineering, University of Utah, Salt Lake City, Utah 84112, USA

Dynamics of charge-density fluctuations in a system of two tunnel-coupled wires contains two diffusion modes with dispersion iω = Dq 2 and iω = Dq 2 + 2 τ t , where D is the diffusion coefficient and τ t is the tunneling time between the wires. The dispersion of corresponding spin-density modes depends on magnetic field as a result of spin precession with Larmour frequency, ω L . The presence of two modes affects the shape of the Hanle curve describing the spin-dependent resistance, R , between ferromagnetic strips covering the non-magnetic wires. We demonstrate that the relative shapes of the R ω ( L )-curves, one measured within the same wire and the other measured between the wires, depends on the ratio τ /τ t s , where τ s is the spin-diffusion time. If the coupling between the wires is local, i.e. only at the point x = 0, then the difference of the shapes of intra-wire and inter-wire Hanle curves reflects the difference in statistics of diffusive trajectories which 'switch' or do not switch near x = 0. When one of the coupled wires is bent into a loop with a radius, a , the shape of the Hanle curve reflects the statistics of random walks on the loop. This statistics is governed by the dimensionless parameter, a √ Dτ s .

PACS numbers: 72.15.Rn, 72.25.Dc, 75.40.Gb, 73.50.-h, 85.75.-d

## I. INTRODUCTION

Spin-orbit interaction is the origin of spin dephasing in semiconductors and metals. On a microscopic level, a finite spin-relaxation time, τ s , results from the momentum-dependent spin-orbit term in the Hamiltonian of a free electron in combination with scatteringinduced momentum relaxation. 1 In metals, the origin of spin dephasing is spin-dependent impurity scattering.

Correspondingly, the second strip, the detector, imitated the analyzer of the emitted light. The characteristic measured in Ref. 6 was the nonlocal resistance, R , which is the ratio of the voltage, generated between the channel and detector, to the current passed through the injector into the channel.

First experimental studies 2-4 of spin relaxation in semiconductors were carried out more than four decades ago. They were based on the notion that if the photoexcited electron has its spin pointing along the x -axis, then in magnetic field, ω L , directed along the z -axis, the projection, S x ( ), t evolves as S x ( ) = cos t ω t L exp( -t/τ s ) which is the result of the Larmour precession. Since the time between generation and recombination is much longer than τ s , so that the spin evolution is completed by the moment of recombination, then the polarization of the luminescence is proportional to ∫ ∞ 0 dtS x ( ), i.e. t

$$\mathcal { P } ( \omega _ { L } ) = \frac { \mathcal { P } ( 0 ) } { 1 + \omega _ { L } ^ { 2 } \tau$$

Numerous experimental measurements reported to date can be fit very accurately with Lorentzian Hanle profile Eq. (1), and when they do not, see e.g. Ref. 5, the deviations reflect the peculiarity of the recombination process.

In the pioneering papers Refs. 6, 7 it was demonstrated that, aside from optics, the underlying physics of spin relaxation manifests itself in transport experiments. The structure fabricated and measured in Ref. 6 represented an aluminum wire with two cobalt ferromagnetic strips on the top. The first strip, injector, played the role of circular-polarized excitation light in optics, in the sense, that it supplied spin-polarized electrons into the wire.

Similar to polarization of luminescence, P ( ω L ), the nonlocal resistance is suppressed with external field, ω L . There is, however, a fundamental difference between the dependencies P ( ω L ) and R ω ( L ). This difference stems from the fact that, in addition to the Larmour precession, formation of nonlocal resistance involves diffusion of carriers over the distance, L , between injector and detector. This diffusion is routinely incorporated into the theory by multiplying cos ω t L exp( -t/τ s ) by a diffusion propagator, P L ( ), t and only subsequently integrating over time. In one dimension, P L ( ) has the form t

$$P _ { L } ( t ) = \frac { 1 } { \left ( 4 \pi D t \right ) ^ { 1 / 2 } } \exp \left ( - \frac { L ^ { 2$$

where D is the diffusion coefficient.

The nonlocal resistance, R ω ( L ), calculated with the help of propagator Eq. (2), is also called the Hanle curve in the literature. The expression for R ω ( L ) contains two unknown parameters, τ s and D . Still it appears that the scores of experimental data accumulated to date can be fitted very accurately with this expression. It is apparent that the shapes of R ω ( L ) is different for 'short' L glyph[lessmuch] ( Dτ s ) 1 / 2 and long L glyph[greatermuch] ( Dτ s ) 1 / 2 samples. This difference in shapes was pointed out already in the seminal paper Ref. 6, where the two samples measured had the lengths L = 50 µ m and L = 300 µ m.

Experimental studies of nonlocal spin transport became a hot topic in 2001 when the measurements of R ω ( L ) were reported 8 for small samples with L ∼ 0 5 . µ m.

FIG. 1: (Color online) (a) Two-channel spin-transport device. The injector is located in the left channel. Two detectors in the left and right channels are located at the same distance, L , from the injector. An electron reaches the first detector by diffusion and the second detector by a combined diffusiontunneling process; (b) In contrast to (a) two channels are coupled locally at the point x = 0; (c) the second wire is bent into a loop. Electron diffusion trajectories encircle the loop several times before the spin polarization is 'forgotten'; (d) two branches, E k ( x ), of the energy spectrum of the tunnel-coupled wires. At small momenta, the splitting, ∆ , is determined by t tunneling, while at large k x the splitting ∆ s is dominated by the spin-orbit coupling in the wires.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0505v1-figure-3b6a31c5-463c-49a6-815b-79cc705cd40e.jpg)

Small structures are appealing for information-storage applications. Indeed, the effect of sign reversal of nonlocal resistance upon reversal of the magnetization of the detector allows one to view the detector as an element of information storage. For this reason, the R ω ( L ) -measurements in spin transport were carried out since 2001 on various structures with L in the micrometer range and with materials of non-magnetic channels ranging from Si and GaAs, see e.g. Refs. 9 and 10, to graphene 11 and organic materials. 12

It turns out that the scope of experimental results on nonlocal resistance are described by the drift-diffusion theory with remarkable precision and including finest details, see e.g. Refs. 13-16. For example, in long samples L glyph[greatermuch] ( Dτ s ) 1 / 2 , the theory predicts several zeros in R ω ( L ) dependence. The origin of these zeros is that during the time L /D 2 of travel between the injector and detector the spin can precess by 2 π , 4 π , and so on. Clearly, the values of R ω ( L ) between these zeros fall off dramatically. Then the number of the zeros observed in experimental R ω ( L ) attests to the accuracy with which the theory captures the process of spin transport. Usually, only the first zero is resolved in experiment. However, the very recent data in Ref. 16 exhibits the second zero as well.

To illustrate the accuracy with which the drift- diffusion theory works for spin-transport devices with variable channel length, in Fig. 2(b), (c) we plot the Hanle shapes measured for two devices fabricated from epitaxial ZnO films. Both devices were fabricated under the same conditions 17 , which included pulsed-laser deposition of ZnO onto a sapphire substrate, deposition of a thin barrier layer of MgO on top, and, finally, the deposition of NiFe electrodes using photolithography and e-beam evaporation. The only difference between the two devices was the distance between the contacts ( L = 90 nm and L = 650 nm) in the four-probe structure, see Fig. 2(a). We see that the seven-fold increase in L changes the shape dramatically, in quantitative agreement with predictions of the drift-diffusion theory.

The fact that the measured R ω ( L ) is so accurately described by the drift-diffusion theory suggests that the shape of the Hanle curve is the characteristics of the spin transport in non-magnetic channel only, and is not affected by the details of injection and detection. It also suggests that description of electron diffusion paths as purely one-dimensional is surprisingly adequate. In fact, in experimental geometries, the length, L , of the channel does not exceed significantly the channel width.

This motivated us to study theoretically the shapes of the Hanle curves in geometries when the transport between the injector and detector does not reduce to a 1D random walk. The results of this study are reported in Sections II and III, where we consider the spin transport along two parallel tunnel-coupled wires and the transport in the special case when one of the wires is bent into a loop. Our objective was to find out whether the statistics of diffusion paths, specific for quasi 1D geometries, can be distinguished in nonlocal spin-transport measurements. Our main message is that, the difference of the Hanle curves measured with two detectors, one located in the same wire as injector and the other located in the neighboring wire, reflects the peculiar statistics of the diffusion paths in a coupled system.

## II. SPIN-DENSITY FLUCTUATIONS IN TUNNEL-COUPLED WIRES

Fluctuations of electron densities, n 1 ( x, t ) and n 2 ( x, t ), in a system of two coupled wires satisfy the system of equations

$$\frac { \partial n _ { 1 } } { \partial t } & = D \frac { \partial ^ { 2 } n _ { 1 } } { \partial x ^ { 2 } } - \frac$$

where τ t is the inter-wire tunneling time. We assume that the wires are disordered so that τ t glyph[greatermuch] τ , where τ is the disorder-scattering time. The latter condition implies that the tunnel splitting, ∆ , of the spectra of the wires t in the absence of disorder is much smaller than τ -1 t . In

FIG. 2: [Color Online] (a) Schematic view of a four-terminal device used for nonlocal spin-transport measurements 17 . Epitaxial ZnO film of a thickness 200 nm, deposited on a sapphire substrate, is spaced from a NiFe layer by a thin MgO barrier; (b) and (c) are the Hanle curves measured for the channel length L = 90 nm and L = 650 nm, respectively. The insets show the theoretical fits plotted from Eqs. (16), (17), (23), and (25). The values τ s used for both fits are the same, so that the dimensionless lengths L = L/ √ 4 Dτ s differ 7 times.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0505v1-figure-3e63f366-9b16-4caf-be50-b52308740982.jpg)

this limit, the expression for τ t reads 18

$$\tau _ { t } = \frac { 1 } { \Delta _ { t } ^ { 2 } \tau }. \quad \quad \quad ( 4 ) \quad \quad \quad$$

Equations Eq. (3) give rise to two diffusion modes with dispersions

$$i \omega = D q ^ { 2 }, \quad i \omega = D q ^ { 2 } + \frac { 2 } { \tau _ { t } }$$

corresponding to symmetric and antisymmetric distributions of densities, respectively. So that the actual distributions n 1 ( x, t ) and n 2 ( x, t ) are linear combinations of the two modes. If electrons are injected into the first wire at x = 0, the combinations satisfying the initial conditions

$$n _ { 1 } ( x, 0 ) = \delta ( x ), \ \ n _ { 2 } ( x, 0 ) = 0. \quad \ \ ( 6 ) \quad \det \epsilon$$

are the sum and the difference of two diffusion modes

$$n _ { 1 } ( x, t ) & = \frac { 1 } { 2 } \ P _ { x } ( t ) \left [ 1 + e ^ { - 2 t / \tau _ { t } } \right ], \quad \text{they} \\ n _ { 2 } ( x, t ) & = \frac { 1 } { 2 } \ P _ { x } ( t ) \left [ 1 - e ^ { - 2 t / \tau _ { t } } \right ], \quad \text{they} \\ \quad \text{$\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\ \quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad$\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $ \quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\ \end{cases}$$

where the diffusion propagator P ( ) is defined by Eq. (2). t x To describe the nonlocal resistance we need the expressions for the spin densities, S 1 ( x, t ) and S 2 ( x, t ), similar to Eq. (7). The system of coupled equations for S 1 ( x, t ), S 2 ( x, t ) has the form

$$\frac { \partial S _ { 1 } } { \partial t } & = \omega _ { L } \times S _ { 1 } - \frac { S _ { 1 } } { \tau _ { s } } + D \frac { \partial ^ { 2 } S _ { 1 } } { \partial x ^ { 2 } } - \frac { 1 } { \tau _ { t } } \left ( S _ { 1 } - S _ { 2 } \right ), \\ \frac { \partial S _ { 2 } } { \partial t } & = \omega _ { L } \times S _ { 2 } - \frac { S _ { 2 } } { \tau _ { s } } + D \frac { \partial ^ { 2 } S _ { 2 } } { \partial x ^ { 2 } } - \frac { 1 } { \tau _ { t } } \left ( S _ { 2 } - S _ { 1 } \right ), \quad ( 8 )$$

and differs from the corresponding equations Eq. (3) describing the charge-density fluctuations in two aspects: both S 1 and S 2 , precess in magnetic field, ω L , and both decay during the spin relaxation time τ s

$$\tau _ { s } = \frac { 1 } { \Delta _ { s } ^ { 2 } \tau },$$

where ∆ s is the spin-orbit splitting of the spectrum in each wire in the absence of disorder. Note that, while both τ t and τ s contain scattering time, the ratio τ /τ s t does not contain disorder, i.e. it is a characteristics of clean wires. The term coupling the wires in the system Eq. (8) has the same form as in the system Eq. (3) since tunneling conserves the spin.

Without boundary conditions, the system Eq. (8) defines four modes

$$\underset { \text{nsel} } { r } \underset { \text{nsel} } { \text{r} } ; \quad i \omega = - D q ^ { 2 } - \frac { 1 } { \tau _ { s } } \pm i \omega _ { L }, \ i \omega = - D q ^ { 2 } - \frac { 1 } { \tau _ { s } } - \frac { 2 } { \tau _ { t } } \pm i \omega _ { L }, \ ( 1 0 )$$

of which the first two correspond to symmetric and the second two to the antisymmetric spin-density fluctuations. With boundary conditions, the solution of the system Eq. (8) can be expressed in terms of the solution Eq. (7) of the system Eq. (3) as follows

$$S _ { 1 } ( x, t ) = s ( t ) n _ { 1 } ( x, t ), \ \ S _ { 2 } ( x, t ) = s ( t ) n _ { 2 } ( x, t ), \ \ ( 1 1 )$$

where the function s ( ) t satisfies the conventional equation of spin dynamics

$$\frac { d s } { d t } = \omega _ { L } \times s - \frac { s } { \tau _ { s } }. \quad \quad \quad ( 1 2 )$$

## III. NONLOCAL RESISTANCES

The initial condition to Eq. (12) is set by the direction of polarization of the injector. We assume that s (0) is directed along the x -axis.

As it is illustrated in Fig. 1, there are two nonlocal resistances: R 11 ( ω L ) is the resistance measured by the detector within the same wire, 1, where polarized electrons are injected, and R 12 is the resistance measured by the detector that covers the wire 2. Within a prefactor they are given by

$$( 7 ) \quad \mathcal { R } _ { 1 1 } = R _ { 0 } L \underset { 0 } { \overset { \cdot } { \lambda } } \overset { \cdot } { \overset { \cdot } { \lambda } } \cos \cdot \underset { 0 } { \lambda } \\ \overset { \cdot } { \mathcal { R } } _ { 1 1 } = R _ { 0 } L \underset { 0 } { \overset { \infty } { \sum } } \overset { d t } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \infty } { \sum } } \overset { d t } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \infty } { \sum } } \overset { d t } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \infty } { \sum } } \overset { d t } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \infty } { \sum } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cdot } { \mathcal { R } } } \overset { \cdot } { \overset { \cd. } { \mathcal { R } } } \overset { \cdot } { \overset { \cd. } { \mathcal { R } } } \overset { \cdot } { \overset { \cd. } { \mathcal { R } } } \overset { \cdot } { \overset { \cd. } { \mathcal { R } } } \overset { \cdot } { \overset { \cd. } { \mathcal { R } } } \overset { \cdot } { \overset { \cd. } { \mathcal { R } } } \overset { \cdot } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } } } \overset { \cd. } { \overset { \cd. } { \mathcal { R } \ } \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } \ } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } \\ \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \ \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \strut } } \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ }{ \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \underset { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } {  \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \overset { \cd. } { \mathcal { R } \ } { \overset { \cd. } { \mathcal { R } \ } { { \mathcal { R } \ } { \overset { \cd. } { \$$

In Eq. (13) it is implicit that the magnetization of the detector is also along the x - axis. In some experiments, say Ref. 19 , the spin transport was studied for polarization of the detector along the y -axis. The corresponding expression for nonlocal resistance reads

$$\xi _ { 1 }, \quad ( 8 ) \quad \widetilde { \mathcal { R } } _ { 1 1 } = R _ { 0 } L \underset { 0 } { \overset { \times } { \int } } \frac { d t } { \tau _ { s } } \ S _ { 1 y } ( L, t ), \ \widetilde { \mathcal { R } } _ { 1 2 } = R _ { 0 } L \underset { 0 } { \overset { \times } { \int } } \frac { d t } { \tau _ { s } } \ S _ { 2 y } ( L, t ).$$

Our goal is to find the expressions for R 11 ( ω L ) and R 12 ( ω L ) for two tunnel-coupled wires. One can see that the coupling strength, τ -1 t , enters into the formulas for R 11 ( ω L ), R 12 ( ω L ) through the last terms in Eq. (3). These terms decay exponentially with time. We also notice that s ( ) t which satisfies Eq. (12) is also an exponential function of time. This observation allows one to express R 11 ( ω L ), R 12 ( ω L ) with tunneling through nonlocal resistance, R ω ( L ), in the absence of tunneling .

Setting τ t = ∞ and substituting

$$s ( t ) = e ^ { - t / \tau _ { s } } \left ( i \, \cos \omega _ { L } t + j \, \sin \omega _ { L } t \right ) \quad ( 1 5 )$$

into Eq. (11) and subsequently into Eq. (13) we restore a standard expression for the Hanle profile of a single channel

$$R ( \omega _ { L }, \tau _ { s } ) = R _ { 0 } F ( \omega _ { L }, \tau _ { s } ) \text{ \quad \ \ } ( 1 6 )$$

where the dimensionless function F ω ( L , τ s ) is defined as

$$F ( \omega _ { L }, \tau _ { s } ) = \frac { L } { \tau _ { s } } \underset { 0 } { \overset { \infty } { \int } } d t \cos \omega _ { L } t \ e ^ { - t / \tau _ { s } } P _ { L } ( t ), \quad ( 1 7 )$$

so that R 0 has the dimensionality of the resistance.

Then, in terms of the function R ω ( L , τ s ), the final result for nonlocal resistances can be presented as

$$\mathcal { R } _ { 1 1 } & = \frac { 1 } { 2 } \left [ R ( \omega _ { L }, \tau _ { s } ) + \frac { \tilde { \tau } _ { s } } { \tau _ { s } } R ( \omega _ { L }, \tilde { \tau } _ { s } ) \right ], \\ \mathcal { R } _ { 1 2 } & = \frac { 1 } { 2 } \left [ R ( \omega _ { L }, \tau _ { s } ) - \frac { \tilde { \tau } _ { s } } { \tau _ { s } } R ( \omega _ { L }, \tilde { \tau } _ { s } ) \right ], \quad \ \ ( 1 8 )$$

where ˜ τ s is an effective spin relaxation time

$$\tilde { \tau } _ { s } = \frac { \tau _ { s } \tau _ { t } } { 2 \tau _ { s } + \tau _ { t } } \quad \quad \quad ( 1 9 ) \quad \text{FI}$$

which includes tunneling an is shorter than τ s . Modifications of the Hanle curves due to tunneling are analyzed below.

## A. Limiting cases

As it was mentioned in the Introduction, the shape of the Hanle curve for a single wire is governed by the dimensionless length

$$\mathcal { L } = \frac { L } { \sqrt { 4 D \tau _ { s } } }. \quad \quad \ \ ( 2 0 ) \quad \ \ t a \, \\ \quad \ \ ( 1 )$$

i . It is apparent that when both wires are long, L glyph[greatermuch] 1, the shapes of the curves R 11 and R 12 do not differ significantly, since a typical electron will have enough time to tunnel before it reaches one of two detectors.

ii . It is also obvious on general grounds that when, the tunneling time is much shorter than the spin-relaxation time, τ t glyph[lessmuch] τ s , the nonlocal resistance R 11 exceeds R 12

a.

FIG. 3: [Color online] Hanle curves R 11 ( ω L ) (blue) and R 12 ( ω L ) (green) measured by two detectors, one located in the same wire as injector and the located in the neighboring wire. Dashed curves are plotted directly from Eq. (18) using the definition Eq. (16); the solid curves are plotted from the assymptotic expansions, R 11 is given by Eq. (17) and R 12 is given by Eq. (21). In all graphs the tunneling time, τ t , is 10 τ s . The three panels correspond to different dimensionless lengths, L = L/ √ 4 Dτ s : L = 3 (a), L = 0 45 (b), and . L = 0 1 . (c). All curves for R 12 are multiplied by τ t τ s = 10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0505v1-figure-5adb0a1c-7875-40a1-aaaf-92255580b31c.jpg)

only slightly. This is because the electron gets equally distributed between the wires before the spin precession takes place. Formally, this follows from Eqs. (18) and (19). In the limit τ t glyph[lessmuch] τ s one has ˜ τ s ≈ τ / t 2. The relative difference, ( R 11 -R 12 ) / R 11 , is of the order of ( τ /τ t s ) 1 / 2 .

iii . The opposite limit of weak tunneling between the wires is most insightful. In this limit, we have τ t glyph[greatermuch] τ s , so that only a small portion of electrons injected in the first wire reach the detector in the second wire. This means that R 12 , is much smaller than R 11 . Formally, two terms in Eq. (18) for R 12 nearly cancel each other. However,

R 12 ( ω L ) possesses a distinctive shape. To find this shape we expand Eq. (18) with respect to τ /τ s t and get

$$\mathcal { R } _ { 1 2 } \approx R _ { 0 } \frac { \tau _ { s } } { \tau _ { t } } G ( \omega _ { L }, \tau _ { s } ),$$

where the function G ω ,τ ( L s ) is defined as

$$G ( \omega _ { L }, \tau _ { s } ) = \frac { L } { \tau _ { s } ^ { 2 } } \int _ { 0 } ^ { \infty } d t \, t \cos \omega _ { L } t \, e ^ { - t / \tau _ { s } } P _ { L } ( t ). \quad ( 2 2 ) \quad \begin{array} { c } \wedge \\ \end{array} \begin{array} { c } \wedge \\ \end{array} \end{array}$$

Analytical expressions for F ω ( L , τ s ) and G ω ,τ ( L s ) for arbitrary length can be found using the identities

$$\int _ { \substack { 0 \\ \infty } } ^ { \infty } \frac { d s } { s ^ { 1 / 2 } } \exp \left [ - \frac { 1 } { s } - y s \right ] = \left ( \frac { \pi } { y } \right ) ^ { 1 / 2 } \exp \left [ - 2 y ^ { 1 / 2 } \right ], \quad ( 2 3$$

$$\underset { 0 } { \overset { \infty } { \int } } d s \ s ^ { 1 / 2 } \exp \left [ - \frac { 1 } { s } - y s \right ] = \frac { \pi ^ { 1 / 2 } } { 2 y ^ { 3 / 2 } } \left ( 1 + 2 y ^ { 1 / 2 } \right ) \exp \left [ - 2 y ^ { 1 / 2 } \right ],$$

$$\int \frac { d s } { s ^ { 1 / 2 } } \exp \left [ - \frac { 1 } { s } - y s \right ] & = \left ( \frac { \pi } { y } \right ) ^ { 1 / 2 } \exp \left [ - 2 y ^ { 1 / 2 } \right ], \quad ( 2 3 ) ^ { \quad ( 2 8 ) \ n } \\ \underset { \overset { \infty } { r } } { \int } \dots \quad \Gamma \ 1 \ 1 \ \mu ^ { 1 / 2 } \ 1 \ \dots \ \ d r \quad \dots \ \ r \quad \dots \ \ r$$

and taking the absolute value and the phase of the complex argument, y , to be

$$| y | = \mathcal { L } ^ { 2 } \left ( 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } \right ) ^ { 1 / 2 }, \quad \phi = \arctan \left ( \omega _ { L } \tau _ { s } \right ). \quad ( 2 5 ) \quad \text{the d}$$

In Fig. 3 we plot these functions which represent the Hanle curves for diagonal and nondiagonal resistances for three domains of L . It is seen, Fig. 3(a), that for large length the shapes of both curves are identical. The smaller is the length the more pronounced is the difference between R 11 and R 12 behaviors. The R 12 ( ω L ) curve is significantly narrower than R 11 ( ω L ) for small length, as it is seen in Fig. 3(c). This narrowing originates from the extra factor, t , in the integrand of Eq. (22) compared to Eq. (17) and can be qualitatively interpreted as follows. In order to reach the detector in the second wire, injected electron diffuses along the first wire, tunnels into the second wire, and diffuses there. Narrower shape indicates that reaching the detector in the second wire takes more time than reaching the detector in the first wire to which electron simply diffuses.

iv . Note that the limit of short wires allows a comprehensive analytical study to which we now turn.

In the limit of small wire length, L glyph[lessmuch] 1, analytical expressions for R 11 and R 12 can be obtained for arbitrary relation between τ t and τ s . In this limit, corresponding to | y | glyph[lessmuch] 1 in Eq. (23), the expression for nonlocal resistances, R ω ( L , τ s ) and ˜ ( R ω L , τ s ) of an isolated wire simplify to

$$R ( \omega _ { L }, \tau _ { s } ) = \frac { R _ { 0 } \mathcal { L } } { \sqrt { 2 } } \frac { \sqrt { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } + 1 } } { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } }, \quad ( 2 6 )$$

$$\tilde { R } ( \omega _ { L }, \tau _ { s } ) = \frac { R _ { 0 } \mathcal { L } } { \sqrt { 2 } } \frac { \stackrel { \bullet } { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } - 1 } } { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } }. \quad ( 2 7 ) \quad \text{the}.$$

Substituting Eq. (26) into Eq. (18), we get

$$& \begin{pmatrix} \mathcal { R } _ { 1 1 } ( \omega _ { L } ) \\ \mathcal { R } _ { 1 2 } ( \omega _ { L } ) \end{pmatrix} = \frac { R _ { 0 } \mathcal { L } } { 2 \sqrt { 2 } } \\ & \quad \times \left [ \frac { \sqrt { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } } + 1 } { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } } \ \pm \sqrt { \frac { \tilde { \tau } _ { s } } { \tau _ { s } } } \ \frac { \sqrt { \sqrt { 1 + \omega _ { L } ^ { 2 } \tilde { \tau } _ { s } ^ { 2 } } } + 1 } { \sqrt { 1 + \omega _ { L } ^ { 2 } \tilde { \tau } _ { s } ^ { 2 } } } \right ]. \\ & \quad \end{pmatrix}.$$

The second term in Eq. (28) is responsible for the difference between R 11 and R 12 . It is apparent that this difference is maximal when ˜ τ s ≈ τ s , i.e. when the tunneling time is long. In the latter case we can simplify Eq. (28) further by expanding with respect to τ /τ s t

$$\mathcal { R } _ { 1 2 } & = \frac { R _ { 0 } \mathcal { L } } { 2 \sqrt { 2 } } \frac { \tau _ { s } } { \tau _ { t } } \frac { 1 } { ( 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } ) ^ { 3 / 2 } } \\ & \times \frac { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } + 1 - \omega _ { L } ^ { 2 } \tau ^ { 2 } } { \sqrt { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } + 1 } }. \quad \quad ( 2 9 )$$

Now R 12 is a function of a single argument, ω τ L s . We see that, while R 11 falls off at large ω τ L s as ( ω τ L s ) -1 / 2 , the decay of R 12 is much faster, as ( ω τ L s ) -3 / 2 . Besides, the distinctive feature of R 12 is that it passes through zero at ω τ L s = √ 3. Overall, see Fig. 3, despite the length is small, the behavior of non-diagonal resistance R 12 resembles the shape of the Hanle curve for a long wire, L glyph[greatermuch] 1.

## IV. LOCAL COUPLING

A different arrangement of two coupled wires is shown in Fig. 1(b). Electron injected into the first wire can cross into the second wire only through a narrow bridge at x = 0. This means that, while R 11 is constituted by all diffusive trajectories in the first wire, the contribution to R 12 comes from a subset of diffusive trajectories which visit the point of contact. More precisely, the bridge serves as 'weak' boundary condition for the diffusion equation. We are going to study how this modification of the diffusion due to crossing into neighboring wire affects the shape of the Hanle curve, R 12 , and compare the result with R 12 ( ω L ) calculated for homogeneous tunneling in the previous Section.

We assume that the coupling via the bridge is weak, so that the concentration n 1 ( x, t ) is given by P x + L 1 ( ). t Presence of the bridge in the diffusion equation Eq. (3) for n 2 ( x, t ) is reflected as a source

$$\frac { \partial n _ { 2 } } { \partial t } - D \frac { \partial ^ { 2 } n _ { 2 } } { \partial x ^ { 2 } } = \frac { l } { \tau _ { t } } \delta ( x ) n _ { 1 } ( 0, t ) \text{ \quad \ \ } ( 3 0 )$$

where l glyph[lessmuch] L represents the width of the bridge. The solution of Eq. (30) can be obtained in a standard way, e.g., by the Fourier expansion of both sides. The expression

a.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0505v1-figure-8f7c78aa-b33f-4874-9fd1-10a637c551c8.jpg)

$$n _ { 2 } ( x, t ) = \frac { l } { \tau _ { t } } \int _ { 0 } ^ { t } d t _ { 1 } P _ { x } ( t - t _ { 1 } ) n _ { 1 } ( 0, t _ { 1 } ). \quad ( 3 1 )$$

Substituting the expression for n 1 (0 , t ) into the integrand, we cast the final result in the form

$$n _ { 2 } ( x, t ) = \frac { l } { \tau _ { t } } \int _ { 0 } ^ { t } d t _ { 1 } \int _ { 0 } ^ { t } d t _ { 2 } P _ { x } ( t _ { 1 } ) P _ { L _ { 1 } } ( t _ { 2 } ) \delta ( t _ { 1 } + t _ { 2 } - t ), \\ \text{$1$-1$} \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{array}$$

which is simply the convolution of two diffusion propagators. Physically, the result Eq. (32) is transparent. It expresses the fact that, to get to the point x in the second wire, electron first diffuses from injector to the bridge and then from the bridge to the point x .

The form Eq. (32) is convenient for the calculation of the nonlocal resistance R 12 ( ω L ). Indeed, for this calculation one has to multiply n 2 ( x, t ) by s ( ), t given by Eq. (15), and integrate over t , which leads to the expression

FIG. 4: [Color online] The difference in the diffusion trajectories in a single wire and in two wires coupled via a bridge manifests itself in the shapes of the Hanle curves. (a) The Hanle curves R 11 ( ω L ) (blue) and R 12 ( ω L ) (green, dashed) are plotted from Eqs. (37) and (38) for a small dimensionless length L = L/ (4 Dτ s ) 1 / 2 = 0 1. . The solid green curve is the Lorentzian asymptote, Eq. (39). The inset shows a comparison of R 12 ( ω L ) for the case of local tunneling (yellow), Eq.(39), and for the case of homogeneous tunneling (brown), Eq. (29); (b) same geometry as in (a). Non-diagonal components of nonlocal resistance corresponding to the perpendicular magnetizations of injector and detector are plotted; (c) the Hanle curve R 11 ( ω L ) (blue) and R 12 ( ω L ) (green) are plotted for the dimensionless length L = 0 45. . The shapes are much closer than in (a).

$$\mathbb { m } _ { 1 } ^ { \cdot } \mathbb { m } _ { 2 } ^ { \cdot } \mathbb { m } _ { 3 } ^ { \cdot } \mathbb { m } _ { 4 } ^ { \cdot } \mathbb { m } _ { 5 } ^ { \cdot } \mathbb { m } _ { 6 } ^ { \cdot } \mathbb { m } _ { 7 } ^ { \cdot } \mathbb { m } _ { 8 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot } \mathbb { m } _ { 9 } ^ { \cdot }. \quad ( 3 3 )$$

Substituting Eq. (32) into Eq. (33) and performing integration over time with the help of the δ -function, we get

$$\mathcal { R } _ { 1 2 } ( \omega _ { L } ) = \frac { R _ { 0 } ( L _ { 1 } + L _ { 2 } ) l } { \tau _ { t } \tau _ { s } } \int _ { 0 } ^ { \infty } d t _ { 1 } \int _ { 0 } ^ { \infty } d t _ { 2 } \ P _ { L _ { 1 } } ( t _ { 1 } ) P _ { L _ { 2 } } ( t _ { 2 } ) e ^ { - ( t _ { 1 } + t _ { 2 } ) / \tau _ { s } } \, ( \cos \omega _ { L } t _ { 1 } \cos \omega _ { L } t _ { 2 } - \sin \omega _ { L } t _ { 1 } \sin \omega _ { L } t _ { 2 } ) \,, \quad ( 3 4 )$$

We now notice that, for both terms in the brackets, the double integral Eq. (34) factorizes into a product of single integrals, which, in turn, can be expressed through the functions R ω ( L ) and ˜ ( R ω L ) for a single wire. The final

expression for R 12 ( ω L ) reads

$$\mathcal { R } _ { 1 2 } = \frac { \tau _ { s } } { \tau _ { t } } \frac { ( L _ { 1 } + L _ { 2 } ) l } { L _ { 1 } L _ { 2 } } \frac { 1 } { R _ { 0 } } \left ( R ( \omega _ { L }, L _ { 1 } ) R ( \omega _ { L }, L _ { 2 } ) - \tilde { R } ( \omega _ { L }, L _ { 1 } ) \tilde { R } ( \omega _ { L }, L _ { 2 } ) \right ).$$

Similarly, for orthogonal magnetizations of injector and detector we get

$$\widetilde { \mathcal { R } } _ { 1 2 } = \frac { \tau _ { s } } { \tau _ { t } } \frac { ( L _ { 1 } + L _ { 2 } ) l } { L _ { 1 } L _ { 2 } } \frac { 1 } { R _ { 0 } } \left ( R ( \omega _ { L }, L _ { 1 } ) \tilde { R } ( \omega _ { L }, L _ { 2 } ) + \tilde { R } ( \omega _ { L }, L _ { 1 } ) R ( \omega _ { L }, L _ { 2 } ) \right ).$$

According to Eq. (35), R 12 depends on both the position of the bridge and the position of the detector. In Appendix A we demonstrate that the dependence of the concentration, n 2 ( L , t 2 ), on the position of the bridge drops out, so that R 12 only depends only on the distance, L = L 1 + L 2 between the injector and the detector. This observation allows L 2 to be set to zero in Eq. (35); correspondingly L 1 should be set equal to the total length L . Technically, this implies that we can use the short-distance asymptotes, Eq. (26) and Eq. (27), while for R ω ( L , L ) and ˜ ( R ω L , L ) the general expressions should be used. These general expressions 8,20 follow from Eq. (23):

$$R ( \omega _ { L }, L ) = \left ( \frac { \pi } { | y | } \right ) ^ { 1 / 2 } \exp \left [ - 2 | y | ^ { 1 / 2 } \cos \frac { \phi } { 2 } \right ] \cos \left ( \frac { \phi } { 2 } + 2 | y | ^ { 1 / 2 } \sin \frac { \phi } { 2 } \right ),$$

where the magnitude, | y | , and phase, φ , are defined by Eq. (25). The corresponding expression for ˜ ( R ω L , L ) differs from Eq. (37) by the replacement of cos with sin in the second factor.

Summarizing, in the geometry of two wires with a bridge, the Hanle curve measured by the first detector is described by R ω ( L , L ), Eq. (37), while the Hanle curve measured by the second detector has the shape given by

$$\mathcal { R } _ { 1 2 } \infty \, \frac { \tau _ { s } } { \tau _ { t } } \, \frac { l } { L } \, \left ( \frac { \sqrt { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } + 1 } } { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } } R ( \omega _ { L }, L ) - \frac { \sqrt { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } - 1 } } { \sqrt { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } } \tilde { R } ( \omega _ { L }, L ) \right ).$$

The most dramatic difference in the shapes of two Hanle curves emerges in the limit of short wires, L glyph[lessmuch] 1. Substituting Eqs. (26) and (27) into Eq. (38) we get the amusingly simple expressions for R 12 , ˜ R 12

the wire while the detector is located in the loop. The spin-transport equation has a form

$$\mathcal { R } _ { 1 2 } = \frac { 1 } { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } }, \quad \tilde { \mathcal { R } } _ { 1 2 } = \frac { \omega _ { L } \tau _ { s } } { 1 + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } }. \quad ( 3 9 ) \quad \text{ \quad \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$ \quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad\text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\ \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$$$

This means that, while the first detector measures the shape Eq. (26), the second detector measures a simple Lorentzian Eq. (1), as in optical measurements. In Fig. 4(a), this difference in shapes is illustrated graphically. We see that R 12 in the second wire is not only narrower, but also possesses a distinctively different shape. Fig. 4(b) illustrates that the difference in the shapes of the two Hanle curves gradually vanished as the wires get longer. Qualitatively this can be understood from Eq. (38). In a long wire the Hanle curve is narrow. This allows us to set ω τ L s glyph[lessmuch] 1 in the prefactors in the brackets. Then the first prefactor close to 1, while the second prefactor is much smaller. Thus we conclude that the ratio of R 12 to R 11 is approximately constant.

## V. COUPLING OF A WIRE TO THE LOOP

As a last example of the modification of the Hanle profile with restricted geometry consider a loop tunnelcoupled to a wire, see Fig. 1. The injector is located in

$$\tilde { \Phi } _ { \mathfrak { m } } \quad \frac { \partial } { \partial t } S ( \theta, t ) = \omega _ { L } \times S ( \theta, t ) + \frac { D } { a ^ { 2 } } \frac { \partial ^ { 2 } } { \partial \theta ^ { 2 } } S ( \theta, t ) - \frac { S ( \theta, t ) } { \tau _ { s } }, \ ( 4 0 )$$

where θ is the azimuthal coordinate and a is the radius of the loop. This equation also allows a factorization: S ( θ, t ) = s ( ) t n θ, t ( ), where s ( ) satisfies Eq. t (12), while the equation for n θ, t ( ) reads

$$\frac { \partial n } { \partial t } = \frac { D } { a ^ { 2 } } \frac { \partial ^ { 2 } n } { \partial \theta ^ { 2 } }.$$

The solution of this equation satisfying the initial condition, n θ, ( 0) = δ θ ( ), can be presented as a sum of angular harmonics

$$\cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ \cdot \\ } } n ( \theta, t ) = \frac { 1 } { 2 \pi } + \frac { 1 } { \pi } \sum _ { k = 1 } ^ { \infty } \exp \left [ - \frac { D k ^ { 2 } t } { a ^ { 2 } } \right ] \cos ( k \theta ). \quad ( 4 2 )$$

Assuming that the detector is located at θ = π , we find the following expression for the spin density

$$\intertext { n l e - 1 } _ { \substack { \text{$el-} \\ \text{in} } } S _ { x } ( t ) = \frac { 1 } { 2 \pi } \sum _ { k = - \infty } ^ { \infty } ( - 1 ) ^ { k } \exp \left [ - \frac { D k ^ { 2 } t } { a ^ { 2 } } - \frac { t } { \tau _ { s } } \right ] \cos \omega _ { L } t.$$

a.

FIG. 5: [Color online] (a) Evolution of the Hanle curves, R ( ω L ), with radius, a . The curves are plotted from Eq. (45) for three values of the dimensionless loop radius λ -1 = a/ Dτ ( s ) 1 / 2 . Dotted, dashed and solid lines correspond to λ = 1 7, 1, and 0 7, respectively. . . (b) same as (a) for hightemperature domain λ = 0 35 (dotted), 0 32 (dashed), . . and 0 27 (solid). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0505v1-figure-f0201704-011b-4c87-af6c-317e967edf43.jpg)

Integration of S x ( ) t over time yields the nonlocal resistance in the form of the infinite sum

$$R ( \omega _ { L } ) = R _ { 0 } \sum _ { k = - \infty } ^ { \infty } \frac { 1 + \frac { D k ^ { 2 } \tau _ { s } } { a ^ { 2 } } } { \left ( 1 + \frac { D k ^ { 2 } \tau _ { s } } { a ^ { 2 } } \right ) ^ { 2 } + \omega _ { L } ^ { 2 } \tau _ { s } ^ { 2 } } }. \quad ( 4 4 ) \quad \text{is for} \quad \text{for} \quad \text{$\nu_{1}$} \\ \text{has a} \quad \text{indic} \quad \text{$\nu$}$$

Closed expression for R ω ( L ) can be obtained with the use of the identity

$$\sum _ { k = - \infty } ^ { \infty } \frac { ( - 1 ) ^ { k } \, ( 1 + k ^ { 2 } \lambda ^ { 2 } ) } { ( k ^ { 2 } \lambda ^ { 2 } + 1 ) ^ { 2 } + z ^ { 2 } } = \frac { 2 } { ( 1 + z ^ { 2 } ) ^ { 1 / 2 } } & & \text{are n} \\ & \times \left [ \frac { x \sinh ( x ) \cos ( y ) - y \sin ( y ) \cosh ( x ) } { \cosh 2 x - \cos 2 y } \right ], \ \ ( 4 5 ) & & \text{the m}$$

where

$$x = \frac { \pi } { \lambda } \sqrt { \frac { \sqrt { 1 + z ^ { 2 } } + 1 } { 2 } }, \ y = \frac { \pi } { \lambda } \sqrt { \frac { \sqrt { 1 + z ^ { 2 } } - 1 } { 2 } }. \quad ( 4 6 ) \quad \bullet$$

Expressing of nonlocal resistance with the help of Eq. (45) requires the following identifications

$$z = \omega _ { L } \tau _ { s }, \quad \lambda = \frac { ( D \tau _ { s } ) ^ { 1 / 2 } } { a }. \quad \quad ( 4 7 )$$

The prefactor π/λ in Eq. (46) is equal to 2 πa/ (4 Dτ s ) 1 / 2 , which is the circumference of the loop in the units of spin-diffusion length. The identity Eq. (45) suggests that nonlocal resistance depends on parameter y in an oscillatory fashion. To clarify the physical meaning of these oscillations consider the limit of weak fields, z glyph[lessmuch] 1, so that the precession angle of spin, δϕ , during the time τ s is small. In this limit the parameter y can be cast in the form

$$y \left | _ { z \ll 1 } \approx \frac { 2 \pi a } { \left ( D \tau _ { s } \right ) ^ { 1 / 2 } } \delta \varphi.$$

The first factor in Eq. (48) can be interpreted as a number of intervals, each having the length equal to the spindiffusion length, covered by electron before it makes a full loop. Since the spin is rotated by δϕ over each interval, the parameter y can be interpreted as a full rotation angle for the whole loop. Then the periodicity of nonlocal resistance corresponds to this full angle being π , 2 π , and so on.

To interpret the oscillations in strong fields, z glyph[greatermuch] 1, we rewrite the parameter y as

$$y \left | _ { z \gg 1 } \approx \frac { 2 \pi a } { \left ( D / \omega _ { L } \right ) ^ { 1 / 2 } }.$$

The denominator in Eq. (49) has the meaning of the length traveled during one Larmour precession period. The fact that R is sensitive to whether the circumference contains integer or half-integer number of these lengths can be interpreted as an effect of finite step-size in the random walk. Naturally, these oscillations are suppressed exponentially, since at strong fields we have x ≈ y .

Suppose now that the temperature is low, so that τ s is long. This means that, before the spin orientation is forgotten, the particle performs many loops, so that density n θ ( ) is nearly homogeneous. This, in turn, suggests that the Hanle profile is unaffected by the diffusion, and has a Lorentzian shape. Plotting R ( ω L ) from Eq. (45) indicates that Lorentzian shape is achieved only for very small loops, such that λ glyph[greaterorsimilar] 40. For moderate values of λ ∼ 1, i.e. for higher temperatures, the Hanle profiles are non-Lorentzian, but rather resemble R ( ω L ) for a long wire, as illustrated in Fig. 5(a). Finally for 'high' temperatures corresponding to λ ∼ 0 3, . see Fig. 5(b), the Hanle curves develop oscillations discussed above, while the magnitude of R ( ω L ) drops rapidly with λ .

## VI. DISCUSSION

- · The fact that experimental Hanle curves are amazingly robust motivated us to investigate whether the charge-transport characteristics could be inferred from their shapes. Namely, in the system of two coupled wires, the tunneling time, τ t , between the wires is a parameter which does not depend on spin. Still, as it is seen in Fig. 3, the Hanle

curves calculated for a given wire, R 11 , and between the wires, R 12 , have visibly different widths. The difference in widths is governed by the ratio. ∆ t / ∆ , s of fundamental 'band-structure' parameters. Thus, this ratio can be inferred from the comparison of these widths. Also, ∆ s can, in principle, be inferred independently from the shape of R 11 . Possibility of such an extraction of tunnel splitting is facilitated by the fact that τ s falls off with increasing temperature dramatically (as T -3 for the Dyakonov-Perel mechanism ), 1 whereas τ t varies slowly. This rapid change of τ s with temperature allows for a 'dimensional crossover' between 0D and 1D statistics of diffusion paths within the same ring-shaped sample, see Fig. 1(c). This crossover manifests itself not only in the width but also in the shape of Hanle curve, which becomes a Lorentzian at low temperatures.

- · In a sense, our quest to reveal the statistics of diffusion paths through the Hanle curves is in line with attempts taken to unravel this statistics from the weak-localization correction to the conductivity, ∆ σ of a 2D sample 21,22 .
- The Hanle profile comes from multiplying the diffusive propagator by cos ω t L and integrating over time. Similarly, the expression for ∆ σ comes from multiplying the diffusive return probability by cos 2 π Φ( ) where Φ( ) is the flux (in the unit of the t t flux quantum) through the area covered by the diffusing particle after time t and integrating over t . Thus, both R and ∆ σ are essentially the Fourier transforms of the diffusion propagator. The role of the spin-flip time τ s in spin transport is played by the phase-breaking time in magnetoresistance.

In fact, a rapid decay of the phase-breaking time with temperature was also exploited previously in the transport studies 23,24 to demonstrate the dimensional crossover from quasi-2D to purely 3D diffusion. It should be noted, however, that while weak-localization relying on the spatial coherence of electron shows up only at low temperatures, the

- 1 M. I. Dyakonov and V. I. Perel, Sov. Phys. Solid State 13, 3023 (1971).
- 2 R. R. Parsons, Phys. Rev. Lett. 23 , 1152 (1969).
- 3 A. I. Ekimov and V. I. Safarov, Pis'ma Zh. Eksp. Teor. Fiz. 12 , 293 (1970)[JETP Lett. 12 , 198 (1970)].
- 4 A. I. Ekimov, D. Z. Garbuzov, and V. I. Safarov, Pis'ma Zh. Eksp. Teor. Fiz. 13 , 36 (1971) [JETP Lett. 13 , 24 (1971)].
- 5 R. I. Dzhioev, V. L. Korenev, B. P. Zakharchenya, D. Gammon, A. S. Bracker, J. G. Tischler, and D. S. Katzer, Phys. Rev. B 66 , 153409 (2002).
- 6 M. Johnson and R. H. Silsbee, Phys. Rev. Lett. 55 , 1790 (1985).

Larmour precession survives at high temperatures and gives rise to the Hanle curve.

- · Throughout the paper we considered a twowire geometry. Another class of structures to which our results might be applicable is tunnelcoupled graphene layers. It was previously demonstrated 11,15,16,25,26 that a single layer of graphene can be used as a channel for nonlocal spintransport measurements. A possibility to fabricate two tunnel-coupled layers was also demonstrated very recently. 27-30 The structures 27-30 were fabricated in order to realize the vertical gate-controlled graphene heterostructures.
- · We considered the spin-current distribution in a loop geometry. Very recently 31 a measurement of nonlocal spin transport in a loop geometry has been reported. The importance of findings of Ref. 31. is that the result of conversion of a spin current into a charge current was revealed not through the voltage buildup in an open-circuit geometry but rather by directly measuring the circulation of current in the loop.
- · The bridge between two channels shown in Fig. 1(b) can be viewed as a boundary condition for the diffusion equation that changes the randomwalk trajectories leading to a modified shape of the Hanle curves. The origin of such a boundary condition can be simply a finite length of the channel. This situation was recently considered theoretically 32 . In accord with our findings, the result of decreasing the channel length is the crossover of the Hanle shape to a Lorentzian.

## VII. ACKNOWLEDGEMENTS

This work was supported by NSF through MRSEC DMR-1121252. E.M. acknowledges support from the Department of Energy, Office of Basic Energy Sciences, Grant No. DE-FG02-06ER46313.

- 7 M. Johnson and R. H. Silsbee, Phys. Rev. B 37 , 5312 (1988).
- 8 F. J. Jedema, A. T. Filip, and B. J. van Wees, Nature (London) 410 , 345 (2001).
- 9 I. Appelbaum, B. Huang, and D. Monsma, Nature (London) 447 , 295 (2007).
- 10 X. Lou, C. Adelmann, S. A. Crooker, E. S. Garlid, J. Zhang, S. M. Reddy, S. D. Flexner, C. J. Palmstrøm, and P. A. Crowell, Nat. Phys. 3 , 197 (2007).
- 11 N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and B. J. van Wees, Nature (London) 448 , 571 (2007).
- 12 S. Watanabe, K. Ando, K. Kang, S. Mooser, Y. Vaynzof, H. Kurebayashi, E. Saitoh, and H. Sirringhaus, Nat. Phys.

10 , 308 (2014).

- 13 K. Kasahara, Y. Fujita, S. Yamada, K. Sawano, M. Miyao, and K. Hamaya, Appl. Phys. Express 7 , 033002 (2014).
- 14 K. Olejn´ ık, J. Wunderlich, A. C. Irvine, R. P. Campion, V. P. Amin, J. Sinova, and T. Jungwirth, Phys. Rev. Lett. 109 , 076601 (2012).
- 15 T. Yamaguchi, Y. Inoue, S. Masubuchi, S. Morikawa, M. Onuki, K. Watanabe, T. Taniguchi, R. Moriya, and T. Machida, Appl. Phys. Express 6 , 073001 (2013).
- 16 A. Dankert, M. V. Kamalakar, J. Bergsten, and S. P. Dash, Appl. Phys. Lett. 104 , 192403 (2014).
- 17 For details of the fabrication, see M. C. Prestgard and A. Tiwari, Appl. Phys. Lett. 104 , 122402 (2014).
- 18 In a superlattice with a period, a , and miniband width, ∆, the diffusion coefficient is equal to a 2 τ t = a 2 ∆ 2 τ , see C. Mauz, A. Rosch, and P. W¨ olfle, Phys. Rev. B 56 , 10953 (1997) the diffusion coefficient in a multilayer system
- 19 S. A. Crooker, M. Furis, X. Lou, C. Adelmann, D. L. Smith, C. J. Palmstrøm, and P. A. Crowell, Science 309 , 2191 (2005).
- 20 T. Sasaki, T. Oikawa, T. Suzuki, M. Shiraishi, Y. Suzuki, and K. Noguchi, IEEE Trans. Magn., 46 , 1436 (2010).
- 21 G. M. Minkov, A. V. Germanenko, V. A. Larionova, S. A. Negashev, and I. V. Gornyi, Phys. Rev. B 61 , 13164 (2000).
- 22 G. M. Minkov, S. A. Negashev, O. E. Rut, A. V. Germanenko, O. I. Khrykin, V. I. Shashkin, and V. M. Daniltsev, Phys. Rev. B 61 , 13172 (2000).

Teor. Fiz. 85 , 2210 (1983).

- 24 A. M. Gilbertson, A. K. M. Newaz, W.-J. Chang, R. Bashir, S. A. Solin1, and L. F. Cohen, Appl. Phys. Lett. 95 , 012113 (2009).
- 25 C. J´zsa, M. Popinciuc, N. Tombros, H. T. Jonkman, and o B. J. van Wees, Phys. Rev. Lett. 100 , 236603 (2008); Phys. Rev. B 79 , 081402 (2009).
- 26 M. H. D. Guimar˜es, A. Veligura, P. J. Zomer, T. Maassen, a I. J. Vera-Marun, N. Tombros, and B. J. van Wees, Nano Lett. 12 , 3512 (2012).
- 27 L. Britnell, R. V. Gorbachev, R. Jalil, B. D. Belle, F. Schedin, A. Mishchenko, T. Georgiou, M. I. Katsnelson, L. Eaves, S. V. Morozov, N. M. R. Peres, J. Leist, A. K. Geim, K. S. Novoselov, and L. A. Ponomarenko, Science 335 , 947 (2012).
- 28 T. Georgiou, R. Jalil, B. D. Belle, L. Britnell, R. V. Gorbachev, S. V. Morozov, Y.-J. Kim, A. Gholinia, S. J. Haigh, O. Makarovsky, L. Eaves, L. A. Ponomarenko, A. K. Geim, K. S. Novoselov, and A. Mishchenko, Nature Nanotechnology 8 , 100 (2013).
- 29 S. Bertolazzi, D. Krasnozhon, and A. Kis, ACS Nano 7 , 3246 (2012).
- 30 S. Kim, D. H. Shin, C. O. Kim, S. S. Kang, J. M. Kim, C. W. Jang, S. S. Joo, J. S. Lee, J. H. Kim, S.-H. Choi, and E. Hwang, ACS Nano, 7 , 5168 (2013).
- 31 Y. Omori, F. Auvray, T. Wakamura, Y. Niimi, A. Fert, and Y. Otani, arXiv:1405.2100.
- 23 A. K. Savchenko, A. S. Rylik, and V. N. Lutskii, Zh. Eksp.
- 32 M. Wojtaszek, I. J. Vera-Marun, and B. J. van Wees, arXiv:1404.6276.

## Appendix A:

In order to substantiate the statement made in Sect. III that the concentration profile in the second wire,

$$n _ { 2 } ( L _ { 2 }, t ) \in \int _ { 0 } ^ { t } d t _ { 1 } \ P _ { L _ { 2 } } ( t - t _ { 1 } ) P _ { L _ { 1 } } ( t _ { 1 } ),$$

does not depend on the position, L 1 , of the bridge, it is convenient to use the Fourier representation of the diffusive propagators P L 1 and P L 2 . In this representation Eq. (A1) acquires the form

$$n _ { 2 } \in \int _ { 0 } ^ { t } d t _ { 1 } \, \left ( \int \frac { d q _ { 2 } } { 2 \pi } \exp \left [ - D q _ { 2 } ^ { 2 } ( t - t _ { 1 } ) + i q _ { 2 } L _ { 2 } \right ] \int \frac { d q _ { 1 } } { 2 \pi } \exp \left [ - D q _ { 1 } ^ { 2 } t _ { 1 } + i q _ { 1 } L _ { 1 } \right ] \right ).$$

Performing the time integration, we get

$$n _ { 2 } \in \int \frac { d q _ { 1 } } { 2 \pi } \int \frac { d q _ { 2 } } { 2 \pi } \exp \left [ i q _ { 1 } L _ { 1 } + i q _ { 2 } L _ { 2 } \right ] \left ( \frac { \exp [ - D q _ { 1 } ^ { 2 } t ] - \exp [ - D q _ { 2 } ^ { 2 } t ] } { q _ { 1 } ^ { 2 } - q _ { 2 } ^ { 2 } } \right ).$$

The sum, L = L 1 + L 2 , which is the total length, does not depend on the position of the bridge, while the difference l = L 1 -L 2 is fully determined by the position of the bridge. In order to decouple L and l , we introduce new variables

$$u = q _ { 1 } + q _ { 2 }, \ \ v = q _ { 1 } - q _ { 2 },$$

so that the integral Eq. (A3) assumes the form

$$n _ { 2 } \in \int \frac { d u } { u } \exp \left [ - \frac { D t u ^ { 2 } } { 4 } + i \frac { L } { 2 } u \right ] \int \frac { d v } { v } \exp \left [ - \frac { D t v ^ { 2 } } { 4 } + i \frac { l } { 2 } v \right ] \sinh \left ( \frac { D t u v } { 2 } \right ).$$

Independence of n 2 on the position of the bridge implies that ∂n /∂l 2 = 0. Differentiating Eq. (A5) with respect to l , we get

$$\frac { \partial n _ { 2 } } { \partial l } \, \infty \, \int \frac { d u } { u } \exp \left [ - \frac { D t u ^ { 2 } } { 4 } + i \frac { L } { 2 } u \right ] \int d v \, \exp \left [ - \frac { D t v ^ { 2 } } { 4 } + i \frac { l } { 2 } v \right ] \sinh \left ( \frac { D t u v } { 2 } \right ).$$

Note that the internal integral in Eq. (A6) can be readily evaluated

$$\int d v \, \exp \left [ - \frac { D t v ^ { 2 } } { 4 } + i \frac { l } { 2 } v \right ] \sinh \left ( \frac { D t u v } { 2 } \right ) \in \exp \left [ \frac { D t u ^ { 2 } } { 4 } \right ] \sin \frac { l u } { 2 }.$$

This allows us to express ∂n /∂l 2 as a single integral

$$\frac { \partial n _ { 2 } } { \partial l } \in \int \frac { d u } { u } \exp \left [ \frac { i L u } { 2 } \right ] \sin \frac { l u } { 2 } = \frac { 1 } { 2 } \int d u \frac { \sin \frac { ( L + l ) u } { 2 } - \sin \frac { ( L - l ) u } { 2 } } { u }.$$

From the identity

$$\int _ { - \infty } ^ { \infty } \frac { d s } { s } \sin \alpha s = \pi \sin ( \alpha ),$$

we conclude that indeed ∂n /∂l 2 is zero as long as l &lt; L . Thus the concentration, n 2 ( L , t 2 ), does not depend on the position of the bridge only when the bridge is located between the injector and the detector.
R. C. Roundy M. C. Prestgard A. Tiwari E. G. Mishchenko M. E. Raikh
Manifestation of two-channel nonlocal spin transport in the shapes of the Hanle curves
Dynamics of charge-density fluctuations in a system of two tunnel-coupled
wires contains two diffusion modes with dispersion iw=Dq^2 and iw
=Dq^2+2/tau_t, where D is the diffusion coefficient and tau_t is the tunneling
time between the wires. The dispersion of corresponding spin-density modes
depends on magnetic field as a result of spin precession with Larmour
frequency, w_L. The presence of two modes affects the shape of the Hanle curve
describing the spin-dependent resistance, R, between ferromagnetic strips
covering the non-magnetic wires. We demonstrate that the relative shapes of the
R(w_L)-curves, one measured within the same wire and the other measured between
the wires, depends on the ratio tau_t/tau_s, where tau_s is the spin-diffusion
time. If the coupling between the wires is local, i.e. only at the point x=0,
then the difference of the shapes of intra-wire and inter-wire Hanle curves
reflects the difference in statistics of diffusive trajectories which "switch"
or do not switch near x=0. When one of the coupled wires is bent into a loop
with a radius, a, the shape of the Hanle curve reflects the statistics of
random walks on the loop. This statistics is governed by the dimensionless
parameter, a/(D tau_s)^(1/2).