Jul 25, 2014 - For electronic circuits, the standard Johnson-. Nyquist noise is ..... by âim, where m is the fourth Majorana fermion (needed anyway to construct the ... can be obtained as a saddle-point solution in the path- ..... physics, repr. ed
May 20, 2005 - 840. 845. 850. 0.1. 1. 10. 100. 1.510 1.512 1.514 1.516 1.518 1.520. -5. 0. 5. 10. 15 ....  B. Dong, H. L. Cui, and X. L. Lei, Phys. Rev. Lett. 94,.
Mar 25, 2011 - In this review, we first introduce the basic principles of SNS and the ... 17. âCorresponding author. Email address: [email protected] (Jens HÃ¼bner). 5 Applications. 18 ..... function of the probe laser spot size A for a
Jul 3, 2013 - The major advantage of spin noise spectroscopy4â8 to all other probes of semiconductor spin dynamics (e.g., op- tical pump-probe technique9 ...
May 10, 2016 - (Dated: May 1, 2016). We report quantum enhancement of Faraday rotation spin noise spectroscopy by polarization squeezing of the probe beam. Using natural abundance Rb in 100 Torr of N2 buffer gas, and squeezed light from a sub-thresho
Jan 25, 2015 - 1Department of Physics, Graduate School of Science, Osaka University, .... FIG. 1. (color online). (a) Current noise at T = 0 for âÂµ = 400 ÂµeV, when a non-magnetic detector electrode is used . The inset shows the schematic ener
Dec 22, 2015 - International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China and. Collaborative Innovation Center of ... CRANN and School of Physics, Trinity College, Dublin 2, Ireland. X.-G. Zhangâ¡ .... (c
Feb 3, 2016 - to the 4th order spin correlators is contained in a simple function H ... correlated events that have a Poisson distribution, which contributes to the ...
Oct 6, 2015 - We demonstrate a method to extend the range of pulsed laser spin noise measurements to long spin lifetimes. We use an ...... trep âfÎ²,Î±. (B2) where we make use of the fact that Î± + Î² = trep. The weighting factor for each series is
Jul 26, 2013 - ing quantum dots (QDs) represent a building block for realization of a ......  H. J. Krenner, M. Sabathil, E. C. Clark, A. Kress, D. Schuh, M.
of the Angstrom scale. The bulk of the experiments per- formed to date have been on magnetically ordered states, such as antiferromagnets. Here we propose to ...
Jul 10, 2014 - spin fluctuations from those singly-charged QDs hav- ing charged exciton (trion) optical transition energies close to the laser energy. Consider ...
Mar 28, 2017 - The theoretical description takes into account generation and relaxation of charge carriers in the quantum dot system. The spin noise is ... system by non-resonant excitation affects the noise signal monitoring the ground state ... As
coupling (DD) [7, 12]. It is based on the application of a sequence of control pulses to the system to effectively isolate it from the environment. Different DD sequences were developed [7, 12â14] and tested experimentally [15â. 20]. Different so
Jul 18, 2017 - They can be thought of as a layer-0 level of protection designed ... Very re- cently, for example, a 10-minute record-breaking coher- ... eating your free lunch. Linearization. Sec. 5. Sec. 4. Figure 1: Estimation demands data. At some
Jan 23, 2017 - wide frequency range, which gives significant advantage in performance and time resolution with respect to ... form (FFT/DFT) [23, 24] of EV and current response is used[25â27]. After which, the investigated .... Eq. (3) with integra
Feb 4, 2015 - 2Spin Optics Laboratory, St. Petersburg State University,. 1 Ul'anovskaya, Peterhof, St. ... ijkl are the second and third order nonlinear sus- ceptibilities, and dots ..... Optics: Principles and Applications (Taylor & Francis,. 2014).
Sep 20, 2011 - taneously as both, probing and pumping beam, respec- tively. The experiments ... and calculations reveal that the observed spin noise un- der resonant probing .... (Color online) Spin noise power density in color coding as a function o
Feb 4, 2010 - ultrafast and depth resolved SNS measurements  to ex- plore the physical origin of the inhomogeneous spin de- phasing in n-doped GaAs ...
Nov 8, 2016 - Optical spin noise spectroscopy: application for study of gyrotropy spatial correlations ... what we call Gaussian beam and introduce the model of polarimetric detector used in further analysis. ... of independent paramagnetic centers a
Jan 1, 2006 - Ultracold alkali atoms [1, 2, 3, 4] provide experimentally accessible model systems for probing quantum states that ... to estimate the expected signal magnitudes for spin noise measurements in ultracold atom systems and to show that ..
Sep 19, 2008 - structures are extremely long even at room temperature. .... which is based on the change of the refractive index due to the fluctuating imbalance of electron ... experimental data. Figure 2 shows a more detailed com- parison of the me
Jun 20, 2013 - We attract attention to the fact that the ultimate (shot-noise-limited) polarimetric sensitivity can be enhanced by orders of magnitude leaving the photon flux incident onto the photodetector on the same low level. This opportunity is
Jul 28, 2017 - the photon-polariton interactions, in the strong light-matter coupling regime, and of the quantum ... short, intermediate and long timescales of spin and pho- ...... tum phase switch between a single solid-state spin and a photon ...
Nonequilibrium spin noise spectroscopy Fuxiang Li Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, NM 87545 USA and Department of Physics, Texas A&M University, College Station, TX 77845 USA
Yuriy V. Pershin Department of Physics and Astronomy, University of South Carolina, Columbia, SC, 29208 USA
arXiv:1307.4426v1 [cond-mat.mes-hall] 16 Jul 2013
Valeriy A. Slipko Department of Physics and Astronomy, University of South Carolina, Columbia, SC, 29208 and Department of Physics and Technology, V. N. Karazin Kharkov National University, Kharkov 61077, Ukraine
N. A. Sinitsyn∗ Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 USA (Dated: November 9, 2018) Spin Noise Spectroscopy (SNS) is an experimental approach to obtain correlators of mesoscopic spin fluctuations in time by purely optical means. We explore the information that this technique can provide when it is applied to a weakly non-equilibrium regime when an electric current is driven through a sample by an electric field. We find that the noise power spectrum of conducting electrons experiences a shift, which is proportional to the strength of the spin-orbit coupling for electrons moving along the electric field direction. We propose applications of this effect to measurements of spin orbit coupling anisotropy and separation of spin noise of conducting and localized electrons.
Introduction. Optical spin noise spectroscopy (SNS)  has been recently introduced as a promising approach for probing local fluctuations of spins in semiconducting materials [2–8], atomic gases [9–14] and quantum dots [15–17]. In particular, if Sz (t) is the time-dependent polarization of spins in a mesoscopic region of a semiconductor, SNS can be used to obtain the spin noise power spectrum [18–23]: Z ∞ P (ω) = 2 dt cos (ωt)hSz (t)Sz (0)i. (1) 0
The advantage of the optical SNS over the other measurement approaches (e.g., optical pump-probe [24, 25] or STM measurements of a single spin ) is usually associated with the minimal energy dissipated, i.e., it can probe spin dynamics at thermodynamic equilibrium. In addition, the SNS allows accumulation of a large statistics that smoothes out the statistical noise in the data, so one can study subtle details at the tails of spin-spin correlators . Sensitivity of this approach is continuously improving, e.g. the recently introduced Ultrahigh Bandwidth SNS can resolve spin correlations with picosecond resolution . At the thermodynamic equilibrium, the fluctuationdissipation theorem predicts that the knowledge of the spin correlator hSz (t)Sz (0)i, which is obtained by SNS, is formally equivalent to the information that can be obtained from a linear response pump-probe measurements. So, for systems at the thermodynamic equilibrium, SNS probes characteristics that, at least in prin-
ciple, can be obtained from more traditional approaches based on characterizing system’s linear response. In this letter we explore the possibility to apply SNS to semiconductors in a non-equilibrium steady state. The behavior of the spin-spin correlator in a non-equilibrium regime may no longer be the subject of the fluctuationdissipation theorem. Hence, even if a perturbation from the equilibrium is weak (in our case it will be a weak electric field that induces an electric current), identifying non-equilibrium contribution to the spin correlator may provide the information about the system that cannot be obtained from the linear response characteristics. For a demonstration, we consider the effect of an electric field on spin fluctuations of conducting elections in 2D electron gas with Rashba and Dresselhaus spin orbit couplings and spin-independent scatterings. We develop an approach that allows us not only to derive the equations for mean spin polarization  but also to relate parameters of spin fluctuations to the shot noise at microscopic scattering events. Stochastic dynamics of spin fluctuations. Consider spin fluctuations from the mean steady state of an electron system. Let ρˆk be the spin density matrix in the momentum space, which is a 2×2 matrix in spin indexes. We assume that the observation region is much larger than the spin diffusion length, so that we can consider dynamics only in the momentum space . The evolution of the spin density matrix in a momentum space volume k is described by the quantum Boltzmann equation: ˆ 0 , ρˆk ] = Iˆcoll , ρˆ˙ k − eE · ∇k ρˆk + i[H
where σ is the vector of Pauli matrices, H is the inplane magnetic field (with absorbed Bohr magneton and g-factor). The spin orbit coupling field can be written as a sum of two parts Ωk = kx Ω1 + ky Ω2 , with:
ωkk0 fk (1 − fk0 )δ(t − t0 )δµν (δkk1 δk0 k01 − δkk01 δk0 k1 ).(8) Let’s introduce variables describing coarse-grained spin characteristics: X X S0 = sk , S1,2 = sk kx,y , (9) k
Ω1 = 2(β x ˆ − αˆ y ), Ω2 = 2(αˆ x − β yˆ),
where α and β are strengths of, respectively, Rashba and Dresselhaus couplings. It is convenient to introduce the spin density sµk = Tr[σµ ρˆk ]/2, where µ = x, y, z, so that eq.(2) can be rewritten as s˙ k − eE · ∇k sk = sk × H + sk × Ωk +
Jk0 →k , (5)
where Jk0 →k is the stochastic spin current, in the momentum space, due to scattering between states k0 and k. For spin conserving scatterings, Jk0 →k = −Jk→k0 . Consider a scattering channel that connects states at k0 and k. Let ak,↑ , ak,↓ be the annihilation operators of an electron with momentum k for spin up and down, respectively. For such scatterings, without spin-flipping, the evolution of spin-up and spin-down electrons during a small time interval t is described by a scattering matrix : ak,↑/↓ (t) ak,↑/↓ (0) Rkk Tkk0 = , (6) Tk0 k Rk0 k0 ak0 ,↑/↓ (t) ak0 ,↑/↓ (0) where Tkk0 and Rkk are, respectively, time-dependent transmission and reflection amplitudes. The spin operator is defined as ˆsk = 12 Ψ†k σΨk , with Ψk = (ak,↑ , ak,↓ )T . If we assume that scattering is weak, then |Tkk0 |2 = |Tk0 k |2 1, |Rkk |2 = 1 − |Tkk0 |2 . The spin current R t 0 operator0 due to such scatterings is defined by: dt Jˆk0 →k (t ) ≡ ˆsk (t) − ˆsk (0). 0 To determine first two cumulants of the spin current , one can take the trace of the first and second powers of ˆsk (t) − ˆsk (0) with the density matrix ρˆst ˆk , k + ρ where ρˆst is the density matrix at the steady state and k ρˆk is due to the currently present spin fluctuation . The steady state density matrix is approximated by a ˆ spin-diagonal matrix ρˆst k ≈ fk 1k , where fk = f (k ) is the Dirac-Fermi distribution over energy k and ˆ1k is a unit matrix in the spin space of fermions with momentum k. This is equivalent to disregarding terms O(E 2 ) and higher order corrections in Ω/F , where F is the Fermi energy, in the final expression for the spin correlator. We should also assume that there are no initial correlations between different phase space volumes and between spin currents at different time moments. Scattering probability, |Tkk0 |2 , is linearly growing with time when energies of k0 and k are the same. Introducing the scattering rate
and the transport life time τtr : X 1 = ωk,k0 [1 − cos(ϕ − ϕ0 )], τtr 0
with ϕ and ϕ0 are angles of, respectively, k and k0 taken at the Fermi surface. By summing over k in Eq. (5) we then obtain S˙ 0 = S0 × H + S1 × Ω1 + S2 × Ω2 . (11) Due to kB TP F , where T is temperature, we can ap2 proximate k kx,y sk ≈ (kF2 /2)S0 , where kF is the Fermi momentum. Multiplying Eq. (5) by kx and ky and summing over k we then find k2 S˙ 1 = −eEx S0 + F S0 × Ω1 − 2 k2 S˙ 2 = −eEy S0 + F S0 × Ω2 − 2
S1 + η1 , τtr S2 + η2 , τtr
where we neglect the terms proportional to H because H 1/τtr . The relaxation terms in (12) and (13) originate from the mean value of the stochastic spin current, and the spin noise terms are defined as: X η1,2 = kx,y (Jk0 →k − hJk0 →k i), (14) k,k0
with their averages zero and correlations: hηiµ ηjν i = kF2 δµν δij DkB T δ(t − t0 )/(2τtr ), i, j = 1, 2, (15) where D is the density of states per spin in the full observation region. Linear dependence on temperature T appears in (15) after integration of fk (1 − fk ) over energy, i.e. it can be traced to the Dirac-Fermi statistics of electrons. Due to short correlation time of spin currents and due to fast relaxation, first harmonics S1µ and S2µ will change at fast time-scales, at which variables S0µ can be considered constant. This allows us to express first harmonics, e.g. S1 , as functions of S0 : S1 = −eEx τtr S0 +
kF2 τtr S0 × Ω1 + κ1 (t), 2
where κ1 is the solution of the equation, κ˙ 1 = − τκtr1 + η1 (t), which describes a noise with correlators hκµ1 (t)i = 0,
hκµ1 (t)κν1 (t0 )i =
kF2 DkB T −|t−t0 |/τtr e δµ,ν . 4 (17)
3 We also note that the correlator (17) has very short decay time τtr (less than 10ps in Rashba 2D electron gas), so at time scales of spin relaxation it can be safely approximated by a white noise hκµ1 (t)κν1 (t0 )i ≈
kF2 DkB T τtr δ(t − t0 )δµ,ν . 2
By similarly working out equations for S2 and substituting results into (11) we obtain equations for slowly changing spin density. In order to simplify notation, we introduce new parameters HSO = 2λeEτtr with λ = q p 2 2 2 2 α + β and E = Ex + Ey , and τ1s = (2λkF )2 τtr as well as the noise variables ξ: ξ = κ1 × Ω1 + κ2 × Ω2 .
Note that hξ x ξ y i = 6 0. Equations for dynamics of S0µ are: Sz S˙ 0z = (S0 × Heff )z − 0 + ξz , τs S y sin(2φ) Sx + ξx , S˙ 0x = (S0 × Heff )x − 0 − 0 2τs 2τs
S x sin(2φ) Sy S˙ 0y = (S0 × Heff )y − 0 − 0 + ξy , 2τs 2τs with Heff = H + HSO , in which HSO = HSO (− sin(θ + φ), cos(θ − φ), 0),
where θ denotes the angle that the in-plane electric field makes with x-axis, and tan φ = β/α. Eqs. (20) and (21) show that τs is a characteristic relaxation time for the out-of-plane spin component of the fluctuation, and HSO is the effective magnetic field which is induced by the electric current. By taking Fourier transform of (20), e.g. S0z (ω) = R Tm /2 dteiωt S0z (t) with Tm being the measurement time, −Tm /2 which is much larger than τs , we obtain the correlator in the frequency domain. When |Heff | 1/τs and a magnetic field is in-plane of the sample, we find that the noise power is given by P (ω) ≡
where ωL = |Heff | and the effective spin relaxation time is τe−1 = τs−1 (3 − sin 2ϕ sin 2φ)/4 in which ϕ is the angle between the direction of Heff and the x-axis. Eq. (22) shows that the shape of the noise power spectrum is Lorentzian but the Larmor frequency is influenced by the electric field and the peak width is renormalized by a factor, which depends on the direction of the magnetic field and the spin orbit coupling anisotropy. Applications. A straightforward experimentally testable prediction of Eqs. (21) and (22) is that the electric field shifts the position of the maximum of the Lorentzian peak by the amount δωL ∼ 2eEλτtr , as we
FIG. 1. (a) The spin noise power spectrum without and with an in-plane electric field (black and red curves, respectively) in strained bulk GaAs. Maximum of the peak is normalized to 1. The peak is shifted by an amount of δωL ≡ ωL − H. Here, H = 100MHz, HSO = 15MHz, τs−1 = 10MHz. (b) The polar plot of the peak shift |δωL | in a 2D electron gas as a function of θ, i.e. the angle of an in-plane electric field with x-axis. Here E = 12V/cm, and the magnetic field is always perpendicular to the electric field. Red, black and green curves correspond to α = λ (β = 0), α = β and α = 0 (β = λ), respectively.
show in Fig. 1(a). The sign of the shift is proportional to e, i.e. depends on the sign of the carriers. Taking the values from  for Rashba coupling in 2D electron gas at GaAs/AlGaAs interface, λR = 1.5 · 10−13 eV · m, relaxation time 1/τtr = 10−3 eV, and assuming the electric field E = 12V/cm, we find HSO ∼ 800MHz. which is comparable to the spin relaxation rate in such systems . Alternatively, a linear Rashba-type spin orbit coupling is induced in bulk 3D GaAs samples at imposed strains, i.e. α ∼ ε ≡ (εxx − εyy ), where εαβ are components of the strain tensor. The effective spin orbit field, HSO , induced by an electric field E = 9V/cm in a 3D GaAs sample with a strain ε = 0.015% was previously determined to be about ∼ 1Gauss by a local Hanle measurement approach . Considering that strains can be increased by an order of magnitude, the field E ∼ 25V/cm should produce the shift of the conducting peak by 15MHz, which would be larger than its width (∼ 10MHz) and hence clearly observable (Fig. 1(a)). The magnitude of the peak-shift effect is sensitive to the anisotropy of the spin orbit coupling, and hence, in a 2D electron gas, depends on the direction of the electric field, as illustrated in Fig. 1(b). In fact, measuring the shift of the Larmor frequency at two transverse directions of the external electric field, one can determine strengths of the Rashba and Dresselhaus couplings separately. Another application of the peak-shift effect can be in studies of localized states at the presence of conducting electrons. At low doping, below the conductinginsulating phase transition [33, 34], there can be donor impurities that are well separated from each other. If the distance between impurities exceeds some critical value R ∼ 200nm, electron hopping between them will be strongly suppressed and localized electron states near such impurities become akin to localized states in quantum dots . Here we predict that spin noise from con-
FIG. 2. Estimate of the noise power spectrum (a) at zero electric field and (b) at an in-plane electric field E = 30V/cm for a GaAs sample with comparable numbers of conducting and strongly localized electrons (n ∼ 3 · 1014 cm−3 ). Spin relaxation times for conducting electrons is τs = 100ns and for localized electrons τloc = 500ns. Electric field E = 30V/cm at strain ε = 0.15% splits the peak at zero frequency into two peaks at ±HSO ∼ 15MHz due to conducting electrons and a peak that remains at zero frequency, which is produced by localized electrons.
ducting and localized electrons can be distinguished due to their different behavior in an applied electric field. Namely, the noise power spectrum for mobile electrons will experience a displacement in the electric field, as explained above, but the power spectrum of localized electron spins will stay intact, as is shown in Fig. 2. In order to estimate this effect, consider a 3D strained GaAs sample with an arbitrary donor concentration n. Let V = 4πR3 /3 be the volume that is needed for a localized electron of a donor impurity to be well separated from other impurities. The probability that a localized state is not overlapped with any other electron state is given by e−nV , so the concentration of impurities whose electron wave functions remain well separated from other donors is given by nloc (n) = ne−nV .
localized and thermally activated conducting electrons, and hence their contribution to the spin noise power, can be made comparable. The physics of spin relaxation of isolated localized electrons is expected to be dominated by the hyperfine coupling in essentially the same way as in the spin of electron-doped InGaAs quantum dots, which was discussed in . The theory  predicts that if the magnetic field is set to zero the localized states of a single donor impurities produce a sharp noise power peak. At low temperatures (below 7-10K), this peak has a nonLorentzian power-law shape at frequencies below 1MHz with a broader shoulder, whose width is determined by the typical strength of the quadrupolar coupling of nuclear spins (Fig. 5b in ). For GaAs, the latter is in the order of several megahertz. At moderately large temperatures (7-30 K), phonon mediated mechanisms of localized spin relaxation make this peak shape Lorentzian . Fig. 2 shows that a reasonably strong electric field is sufficient to shift the peak of conducting electrons and distinguish it from the peak of localized states. Conclusion. We predict that measuring the spin noise power spectrum at steady non-equilibrium conditions is a promising research direction with applications to parameter estimation and uncovering new phenomena. We showed that an electric field leads to a measurable shift of the noise power peak of conducting electrons, which can be used for characterizing the anisotropy of the spin orbit coupling and separating the spin noise of localized states from the spin noise of conducting electrons. Future research directions on the non-equilibrium SNS may include effects of an AC electric field, spin noise measured from optically polarized electrons, studies of high order fluctuation-dissipation relations , and spin noise in the non-Ohmic regime at strong electric fields .
One can find that for R = 200nm and a doping above the metal-insulator transition (n ∼ 1016 cm−3 ), the number of such well localized states is negligibly small. However, at lower doping (n ∼ 1014 cm−3 ), the number of well
Authors thank S. Crooker, D. Smith, A. Saxena, and Yan Li for useful discussions. This work was funded by DOE under Contract No. DE-AC52-06NA25396.
 E. B. Aleksandrov and V. S. Zapasskii, Zh. Eksp. Teor. Fiz. 81, 132 (1981) [JETP 54, 64 (1981)].  M. Oestreich, M. R¨ omer, R. J. Haug, and D. H¨ agele, Phys. Rev. Lett. 95, 216603(2005).  G. M. M¨ uller et al., Phys. Rev. Lett. 101, 206601 (2008).  S. A. Crooker, L. Cheng, and D. L. Smith, Phys.Rev. B 79, 035208 (2009).  Georg M. M¨ uller, Michael Oestreich, Michael R¨ omer, and Jens H¨ ubner, Physica E 43, 569 (2010).  M. R¨ omer, et al, Phys. Rev. B 81 075216 (2010).  Q. Huang and D. S. Steel, Phys. Rev. B 83, 155204 (2011).
 F. Berski, H. Kuhn, J. G. Lonnemann, J. H¨ ubner, and M. Oestreich, Preprint arXiv:1207.0081 (2012); S. Starosielec, and D. H¨ agele, Appl. Phys. Lett. 93, 051116 (2008).  S. A. Crooker, D. G. Rickel, A. V. Balatsky, and D. L. Smith, Nature (London) 431, 49 (2004).  B. Mihaila et al., Phys. Rev. A 74, 043819 (2006).  B. Mihaila et al., Phys. Rev. A 74, 063608 (2006).  G. E. Katsoprinakis, A. T. Dellis, and I. K. Kominis, Phys. Rev. A 75, 042502 (2007).  V. Shah, G. Vasilakis, and M.V. Romalis, Phys. Rev. Lett. 104, 013601 (2010).  V. S. Zapasskii, A. Greilich, S. A. Crooker, Yan Li, G. G.
        
Kozlov, D. R. Yakovlev, D. Reuter, A. D. Wieck, and M. Bayer, Phys. Rev. Lett. 110, 176601 (2013). S. A. Crooker, J. Brandt, C. Sandfort, A. Greilich, D. R. Yakovlev, D. Reuter, A. D. Wieck, and M. Bayer, Phys. Rev. Lett.104, 036601 (2010). Y. Li et al., Phys. Rev. Lett. 108, 186603 (2012). Andreas V. Kuhlmann, et al, Preprint arXiv:1301.6381v1 (2013). M. Braun and J. K¨ onig, Phys. Rev. B 75, 085310 (2007). S. Kos, A. V. Balatsky, P. B. Littlewood, and D. L. Smith, Phys. Rev. B 81, 064407(2010). M. M. Glazov and E. Y. Sherman, Phys. Rev. Lett. 107, 156602(2011). Y. V. Pershin, V. A. Slipko, D. Roy, N. Sinitsyn, Appl. Phys. Lett. 102, 202405 (2013). N. A. Sinitsyn, Yan Li, S. A. Crooker, A. Saxena, D. L. Smith, Phys. Rev. Lett. 109, 166605 (2012). D.Roy, et al, Phys. Rev. B (in press). R.-B. Liu, W. Yao, and L. J. Sham, Adv. Phys. 59, 703 (2010); I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 (2004); R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2007); D. Loss, and D.P. DiVincenzo, Phys. Rev. A 57, 120 (1998).
 M. Furis, D. L. Smith, S. Kos, E. S. Garlid, K. S. M. Reddy, C. J. Palmstrom, P. A. Crowell, and S. A. Crooker, New J. Phys. 9, 347 (2007).  A. V. Balatsky, J. Fransson, D. Mozyrsky, and Y. Manassen, Phys Rev B. 73,184429 (2006); A. V. Balatsky, Y. Manassen, and R. Salem, Phys. Rev. B 66, 195416 (2002).  O. Bleibaum, Phys. Rev. B 72, 075366 (2005).  W. Kohn and J. M. Luttinger, Phys. Rev. 108, 590 (1957).  Ya. M. Blanter, M. Buttiker, Phys. Rep. 336, 1 (2000).  The procedure to obtain higher order correlators is more complex, see e.g. I. Klich, “Quantum Noise in Mesoscopic Systems,” ed. Yu V Nazarov (Kluwer, 2003).  L. Meier, et al, Nature Physics, 3, 650 (2007).  A. Balocchi, et al, Phys. Rev. Lett. 107, 136604 (2011).  B. I. Shklovskii, A. L. Efros, “Electronic properties of doped semiconductors.” Springer, Heidelberg (1984).  F. Gebhardt, “The Mott metal-insulator transition: Models and methods.” Springer Tracts in Modern Physics (1997).  H. Forster and M. Buttiker, Phys. Rev. Lett. 101, 136805 (2008).  M. Furis, D. L. Smith, S. A. Crooker, and J. L. Reno, Appl. Phys. Lett. 89, 102102 (2006).