Experimental test of universal conductance fluctuations by means of wave-chaotic microwave cavities

Please download to get full document.

View again

of 4
All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.
Document Description
Experimental test of universal conductance fluctuations by means of wave-chaotic microwave cavities
Document Share
Document Tags
Document Transcript
   1 Experimental Test of Universal Conductance Fluctuations by means of Wave-Chaotic Microwave Cavities. Sameer Hemmady 1,2,3,4 ,James Hart 1 , Xing Zheng 5 , Thomas M. Antonsen Jr. 1,2,3 , Edward Ott 1,2,3  and Steven M. Anlage 1,2,4 . 1 Department of Physics, University of Maryland, College Park, MD 20742-4111, U.S.A. [Dated : June 26 th  , 2006] . Abstract: The mathematical equivalence of the time-independent Schrödinger equation and the Helmholtz equation is exploited to provide a novel means of studying universal conductance fluctuations in ballistic chaotic mesoscopic systems using a two-dimensional microwave-cavity. The classically chaotic ray trajectories within a suitably-shaped microwave cavity play a role analogous to that of the chaotic dynamics of non-interacting electron transport through a ballistic quantum dot in the absence of thermal fluctuations. The microwave cavity is coupled through two single-mode ports and the effect of non-ideal coupling between the ports and cavity is removed by a  previously developed method based on the measured radiation impedance matrix. The Landauer-Büttiker formalism is applied to obtain the conductance of a corresponding mesoscopic quantum dot device. We find good agreement for the probability density functions of the experimentally derived surrogate conductance, as well as its mean and variance, with the theoretical predictions of Brouwer and Beenakker. We also observe a linear relation between the quantum dephasing parameter and the cavity ohmic loss parameter. Much attention has been focused on the problem of mesoscopic transport through a quantum dot in which a two-dimensional electron gas system contained within an arbitrarily-shaped potential-well boundary is connected to two electron reservoirs through leads– the source ( s ) and drain ( d  ). Recently it has been possible to fabricate quantum dots with low impurity content where the elastic mean free paths of the enclosed electrons are typically much larger than the physical size of the dot [1] . Electron transport through such “ballistic dots” is governed by elastic collisions off the enclosing potential-well boundaries. It has  been observed that the terminal conductance of such dots, defined as )/(ˆ d ss  V V  I G  −=  where s  I   is the source current flowing into the dot and )(  d s  V V   −  is the potential difference between these two leads, exhibits strong, reproducible fluctuations on the order of the quantum of conductance (  heG / 20  = ) [2, 3, 4] . These fluctuations arise from quantum-interference effects due to the phase-coherent electron transport within such dots and have been explained using the hypothesis that the fluctuations are governed by random matrix theory [5]. Similar universal conductance fluctuations (UCF) have also been observed in other systems such as quasi-one-dimensional metal wires [6, 7] . In a quantum dot, this phase coherence is partly lost by opening the system to the outside world during the process of measurement of the conductance. Quantum phase decoherence (dephasing) can also be induced due to the presence of impurities within the dot, thermal fluctuations, or electron-electron interactions, all of which lead to more classical properties for electron transport [8] . Significant theoretical and experimental effort has been devoted to studying the dephasing of the transport electrons in quantum dots [9, 10, 11] . One class of theoretical dephasing models utilizes a fictitious “voltage probe” attached to the dot that has a number of modes φ   N   each with transparency φ  Γ . Electrons that enter one of the modes of this probe are re-injected into the dot with a phase that is uncorrelated with their initial phase, and there is no net current through the fictitious  probe. An alternative model of electron transport employs a uniform imaginary term in the electron potential [12, 13] , leading to loss of probability density with time. It was shown that [14] , as far as the conductance is concerned, these two models yield equivalent predictions in the limit when the number of channels in the dephasing lead ∞→ φ   N   and 0 →Γ φ  , with the product φ φ  γ   Γ=  N   remaining finite (“the locally weak absorbing limit” )   [15] . A similar idea exists for describing ohmic losses in the microwave cavity in terms of non-ideally coupled “parasitic channels” [16] . Since the ohmic losses in the microwave cavity are to good approximation uniformly distributed, we can make use of the equivalence of the imaginary potential and voltage leads models mentioned above to relate the de-phasing parameter employed by electron-transport theory [14]  to the loss parameter of our microwave cavity ( )/( 22 Qk k  n Δ ) [17, 18, 19, 20] . Here, c f k  /2 π  =  is the wavenumber for the incoming frequency  f  and 2 n k  Δ  is the mean-spacing of the adjacent eigenvalues of the Helmholtz operator, 22 k  +∇ , as predicted by the Weyl Formula [21]  for the closed system. The quantity Q represents the loaded quality factor of the cavity. Using the prescription outlined by ref. [14]  we can directly determine the analog of conductance for the microwave cavity and make detailed comparisons of data to theory. We use an electromagnetic analog of a quantum dot in the form of a two-dimensional chaotic microwave resonator. In the case of a cavity thin in one dimension, Maxwell’s equations reduce to a two-dimensional scalar Helmholtz equation. Owing to the analogy between the scalar Helmholtz equation and the Schrödinger equation [22] , the chaotic microwave cavity is an ideal surrogate for a ballistic quantum dot without the complicating effects of thermal fluctuations [4] , Coulomb interactions, or impurities. The microwave analog also permits detailed measurements of the eigenvalues [23, 24, 25] , eigenfunctions [26, 27] , scattering and reaction matrices [18, 19, 20, 28, 29, 30] , in a system where every detail of the potential and the coupling channels can be controlled. Adopting a variant of the Landaeur-Buttiker formalism, the normalized conductance ( o GGG 2ˆ = ) can be expressed in terms of the scattering matrix ⎟⎟ ⎠ ⎞⎜⎜⎝ ⎛ = 22211211 sssss t  of a chaotic cavity when the leads (ports) are perfectly coupled to the cavity [31] . i.e., 222221212211221222212211212 ||||||||2 )||||1)(||||1( || ssss ssss sG −−−− −−−− += . (1) In the time-reversal symmetric case with single-mode leads, Ref. [14]  has shown that as γ   increases the probability 2 Dept. of Electrical and Computer Engineering. 3 Institute for Research in Electronics and Applied Physics. 4 Center for Superconductivity Research. 5 Now at George Washington University.     2 density function of G  (i.e. );(  γ  GP ) becomes more and more sharply peaked around the classical value of 2/1 = G . In the limit of large γ  , an asymptotic analytic expression for );(  γ  GP is  predicted to be [14],   || )||1( 21);(  x e x xGP  − −+=  γ γ   with )2/1(2  −=  G x  γ  . (2) This yields a large- γ   asymptotic expression for the mean and variance of G  which are predicted to be [14],   )(2121 2 − +−>=<  γ γ  OG , (3), )(43)var( 32 − +=  γ γ  OG . (4). Fig.1: (Color Online) (a) PDFs for the normalized conductance );(  γ  GP  obtained from a chaotic cavity for dry-ice case : 4.1-4.7 GHz (hexagons) ( 1.09.0)/( 22 ±=Δ  Qk k  n ; 1.02.11  ±= 〉〈 T  γ  ) and in (b) Loss Case 0 : 16.8-17.6 GHz (stars) ( 1.08.2)/( 22 ±=Δ  Qk k  n ; 1.01.35  ±= 〉〈 T  γ  ); Loss Case 1 : 8.3-9.5 GHz (circles) ( 1.06.6)/( 22 ±=Δ  Qk k  n ; 1.02.82  ±= 〉〈 T  γ  ) and Loss Case 2 : 16.8-17.6 GHz (squares) ( 1.07.21)/( 22 ±=Δ  Qk k  n ; 1.01.272  ±= 〉〈 T  γ  ). Inset : (i) Schematic (not to scale) diagram for the quarter-bow-tie chaotic cavity. The gray jagged shapes are the two metallic perturbations used to perform ensemble averaging. (ii) Schematic of the coupling geometry of each port shown in profile. (iii) Implementation of the outward radiation boundary condition used for the normalization of the measured cavity impedance. The microwave cavity under study is a metallic, air-filled, quarter bow-tie shaped chaotic resonator (inset (i) of Fig. 1(a)) which is quasi-two-dimensional for frequencies below 19.05 GHz. The cavity is driven by two single-mode, coaxial transmission lines whose inner conductor (diameter 2a=1.27mm for Port 1, 2a=0.635mm for Port 2) extends from the top plate of the cavity and makes contact with the bottom plate (inset (ii) of Fig.1(b)). An ensemble data set of one-hundred similar cavities with different internal field configurations is generated by rotating and translating two metallic perturbations, each of which are roughly the size of a wavelength at 5 GHz (inset (i) of Fig. 1(a)), within the cavity volume. This approach of configuration averaging to approximate a pure ensemble average is similar in principle to deforming the shape of the potential-well boundary of a quantum dot as  performed by [4] although in our case the volume of the system is fixed. In addition to the intrinsic ohmic loss in the cavity, the degree of loss can be further increased in a controlled manner by  partially lining the inner side-walls of the cavity with 2 cm-long strips of microwave absorber having uniform spacing. This results in three experimental Loss Cases - Loss Case 0: no absorbing strips, Loss Case 1: 16 absorbing strips, Loss Case 2: 32 absorbing strips. A fourth experimental Loss Case is created by placing the Loss Case 0 cavity in a bath of dry-ice (solid CO 2  at -78.5 o C). This has the effect of slightly increasing the Loss Case 0 cavity Q value (by ~10%). We refer to this case as the “dry-ice case”. A more detailed explanation of our experimental setup and data analysis can be found in  [18, 20] . For our investigation, we generate a large ensemble of 2x2 cavity scattering matrices ( S  t ) for the four loss cases through measurements for different configurations of the perturbers at many frequencies in a range from 3 to 18 GHz (covering ~ 800 modes of the cavity). Using the “radiation impedance” approach [17, 18, 20]  (inset (iii) of Fig.1(b)), the non-ideal coupling details of the two ports are removed to yield an ensemble of normalized 2x2 scattering matrices ( s t ) from which the conductance statistics are derived using Eq. (1). Prior to reporting results for different “data sets” where each data set corresponds to one of our four loss cases and a frequency range typically spanning about 1 GHz, we estimate the value of γ   for each data set. We derive an analytic expression for 〉〈 T   in terms of γ   from Eq. (17a) of [14] , where T   is the absorption probability and is related to the eigenvalues of ss tt † [14, 20] . )))2/()2))2(2((2 )())2(22(4)1(4(( 41 2/ γ ξ γ γ  γ ξ γ γ γ  γ  γ γ γ γ γ γ  −−−+− −+−−+−−=〉〈  − eeeeeeT   (5) where ∫ ∞−− −=  zt  dt t e z )( ξ   is the exponential integral function. By determining the value of 〉〈 T   from the measured data set, Eq.(5) then uniquely determines the corresponding value of γ  ( 〉〈 ≡  T  γ  ). To determine the cavity loss parameter )/( 22 Qk k  n Δ  for our data sets, we employ two procedures. For data sets with 5)/( 22 ≤Δ  Qk k  n , we numerically generate marginal PDFs of the real and imaginary parts of the normalized impedance ( 1 )1)(1(  − −+=  ttttt ss z ) eigenvalues using random-matrix Monte-Carlo simulations with square matrices of size 1000 =  N  , and the value of )/( 22 Qk k  n Δ  in the simulations ranges from 0.1 to 5 in steps of 0.1. We determine the variance ( 2 σ  ) of these numerically generated PDFs and fit it to a polynomial function ))/(( 222 Qk k  n ΔΘ= σ   of order 20. We then determine the variances of the real and imaginary parts of each experimental data set and solve the inverse polynomial function )()/( 2122 σ  − Θ=Δ  Qk k  n  to obtain a unique estimate of )/( 22 Qk k  n Δ  corresponding to each data set. For data sets with 5)/( 22 >Δ  Qk k  n , we use the relation ))/((1 222 Qk k  n Δ= π σ  [17, 20, 32] which has been validated experimentally in Ref. [18] . To report our results, we begin by examining the relationship between the estimated dephasing parameter 〉〈 T  γ   and the estimated cavity loss parameter )/( 22 Qk k  n Δ . By employing a sliding frequency window 1 GHz wide that runs over each of the   3 three Loss Cases - 0, 1, 2 from 3 to 18 GHz, we estimate the value of 〉〈 T  γ   and the corresponding value of )/( 22 Qk k  n Δ  for each window. The comparison is shown as the black circles in the inset of Fig. 2. A linear fit yields the empirical expression )/()1.05.12( 22 Qk k  nT   Δ±= 〉〈 γ   for ~70 points with values for 〉〈 T  γ   ranging from ~11 to ~300. By comparing the Poynting theorem for the electromagnetic cavity with the continuity equation for the probability density in the quantum system, we find )/(4 22 Qk k  n Δ=  π γ  , with ...56.124  = π  . This result can be considered an empirical confirmation of the proposed equivalence of the imaginary potential (uniform volume losses) and de-phasing lead models in the limit considered in [14] . The 1 GHz width of our sliding window was chosen to be large enough to overcome the effects of short-ray paths (which are not removed by only configuration averaging [17, 20] ), but at the same time small enough that the cavity losses can be assumed to be approximately constant over this frequency range. Fig.2: (Color Online) The universal scaling behavior of the conductance distributions is shown. The vertical-axis represents ]),([ 10 γ γ  GP Log  with the corresponding )2/1(2  −=  G x  γ   along the horizontal-axis for three representative data sets consisting of Loss Case 1 : 5.01-6.08 GHz (stars) ( 1.05.4)/( 22 ±=Δ  Qk k  n ; 1.06.56  ±= 〉〈 T  γ  ); Loss Case 1 : 13.6-14.6 GHz (circles) ( 1.03.7)/( 22 ±=Δ  Qk k  n ; 1.06.91  ±= 〉〈 T  γ  ) and Loss Case 2 : 13.6-14.6 GHz (squares) ( 1.07.17)/( 22 ±=Δ  Qk k  n ; 1.05.220  ±= 〉〈 T  γ  ). Inset: The relation between the experimentally determined 〉〈 T  γ   and )/( 22 Qk k  n Δ . In Fig. 1, the experimentally obtained histogram approximation (symbols) to the PDF of the normalized conductance ( );(  γ  GP ) derived from the normalized scattering matrix s t  and Eq. (1) is shown for four cavity data sets. The colored solid lines (magenta, black, green, red) are the asymptotic analytic expression for ),(  γ  GP  (Eq. (2)) with values of γ   that correspond to the estimated 〉〈 T  γ   values obtained from the four cavity data sets. The purple-colored solid line in Fig.1(a) is a random matrix Monte-Carlo simulation for values of 〉〈 T  γ   corresponding to the data set in Fig.1(a). The red error bars (roughly the size of the symbols) in Fig.1 which are representative of the typical statistical binning error of the experimental histograms show that the agreement between the data (shown by the symbols) and the theoretical predictions (shown by the solid curves) improves as the value of 〉〈 T  γ   increases. This is to be expected as Eq. (2) is valid only in the high dephasing limit ( 1 >> γ  ). Similar good agreement between the data and Eq. (2) is obtained for all of the ~40 data sets that we examined in which the frequency ranges and cavity loss cases resulted in an estimate of the 〉〈 T  γ   parameter to be greater than ~18. Fig.3: Evolution of the variance of the experimentally determined );(  γ  GP  distributions (i.e., )var( G ) for increasing values of 〉〈 T  γ   plotted on a logarithmic scale. Inset: Evolution of the mean of the experimentally determined );(  γ  GP  distributions (i.e.,  〉〈 G ) for increasing values of 〉〈 T  γ  . In order to bring out the universal scaling behavior of the );(  γ  GP  distributions (Eqs. (2)) and also to test that these distributions remain strictly non-Gaussian for increasing values of γ   (as predicted by [14] ), we rescale the );(  γ  GP distributions by  plotting ]);([ 10 γ γ  GP Log  versus )2/1(2  −=  G x  γ   in Fig. 2 for three representative data sets with 〉〈 T  γ   ranging from ~56 to ~220. We observe that the three data sets roughly fall on top of each other. The solid blue curve is the theoretical curve, Eq.(2) which is in good agreement with the data. We observe some deviation of the symbols from the theoretical curve near 2 +≅  x . This is attributed to the lack of adequate statistics in the tails of the experimentally determined histogram approximations to the  probability density functions of the conductance. Overall, for values of  x  ranging from 4 −  to 2 + , the agreement is qualitatively good and applies over other data sets where 〉〈 T  γ   ranges from ~18 to 330. The asymmetric (non-parabolic) nature of the experimental data (symbols) confirms that the experimentally obtained );(  γ  GP  remains strictly non-Gaussian and negatively skewed even for large values of γ  , as predicted by [14] .   4 In Fig. 3, we again employ the sliding frequency window of width 1 GHz to test the asymptotic ( >> 1) relations for the mean ><  G  (Eq.3) and variance )var( G  (Eq.4) of );(  γ  GP  as a function of dephasing (loss) parameter γ  . As before, we determine the value of 〉〈 T  γ   for each frequency window data set that runs from 3 to 18 GHz for the three Loss Cases. We then determine the corresponding values of the mean and variance of the corresponding conductance distributions );(  γ  GP  of each frequency window. In the inset of Fig. 3, each star indicates the experimentally estimated mean value of G (i.e,  ><  G ) for the corresponding value of 〉〈 T  γ  . The standard deviation about the experimentally determined mean is of order 5 10 − . We observe that as 〉〈 T  γ   increases, the stars asymptotically approach the classical value of 2/1 >=<  G . The solid black curve represents the leading terms in Eq. (3) .  The circles in Fig. 3 show a similar analysis for the variance ( )var( G ) of the normalized conductance distributions );(  γ  GP  as a function of γ  . The solid black curve represents the leading term in Eq. (4) .  The mean value of the absolute deviation of the circles about the function )(  〉〈 T   f   γ   is about 0.06. We observe that the circles closely follow the functional approximation for the theoretical curve (Eq. (4)) for the range of 〉〈 T  γ   values from ~ 18 to ~330, with no adjustable parameters. The results discussed in this paper provide experimental evidence in support of the theoretical arguments proposed by [14]  and the hypothesis that random matrix theory provides a good description of the conductance fluctuation statistics in a ballistic chaotic quantum-dot in the presence of dephasing. We have shown that in the “locally weak absorbing limit” as discussed in [14] , the dephasing parameter can be related to the cavity loss parameter. We have derived an empirical linear relation between γ   and the cavity loss-parameter )/( 22 Qk k  n Δ  based on our experimental data. The finite conductivity of the metallic walls of the cavity translates to a minimum-possible experimentally accessible value of γ   ~11 for our experiments (at least for the present cavity geometry and temperatures of -78.5 o C and above). We have shown that our experimentally determined conductance distributions and the asymptotic analytic functional forms for the PDF of G ( )( GP ), its mean value (  ><  G ) and variance ( )var( G ) are in good agreement over a broad range of large γ   values. This also establishes the microwave analog as a method to study detailed theories of non-interacting quantum transport and de-coherence in quantum coherent systems. We acknowledge useful discussions with R. Prange and S. Fishman, as well as comments from P. Brouwer, P. Pereyra and T. Seligman. This work is supported by the DoD MURI for the study of microwave effects under AFOSR Grant F496200110374, AFOSR DURIP Grants FA95500410295 and FA95500510240, and by the Israel/U.S.A. Binational Science Foundation.   References: [1] Y. Alhassid, Rev. Mod. Phys. 72 , 895 (2000). [2] B.L. Altshuler and B. D. Simons, “Mesoscopic Quantum Physics”, edited by E. Akkermans, G. Montambaux, J.-L. Pichard and J. Zinn-Justin. Elsevier, Amsterdam (1995). [3] P.A. Lee and A. D. Stone, Phys. Rev. Lett. 55 , 1622 (1985). [4] A. G. Huibers, et al. , Phys. Rev. Lett. 81 , 1917 (1998). [5] M. L. Mehta, “Random Matrices”, 3 rd   Ed. Academic Press. (2004). [6] P. Mohanty and R. A. Webb, Phys. Rev. Lett. 88 , 146601 (2002). [7] N. Agraït, A. L. Yeyati and J. M. Ruitenbeek, Phys. Rep. 377 , 81 (2003); and references therein. [8] E. Joos, et al. , “Decoherence and the Appearance of a Classical World in Quantum Theory”, Springer (December 1, 1996). [9] M. Büttiker, Phys. Rev. B. 33 , 3020 (1986). [10] H. U. Baranger, Phys. Rev. B, 51 , 4703. [11] A. G. Huibers, et al. , Phys. Rev. Lett. 81 , 200 (1998). [12] K.B. Efetov, Phys. Rev. Lett. 74 , 2299 (1995). [13] E. McCann, I. V. Lerner, J.Phys. Cond. Matt. 8 , 6719 (1996). [14] P. W. Brouwer and C. W. J. Beenakker, Phys. Rev. B, 55 , 4695 (1997). [15] M. R. Zirnbauer, Nucl. Phys. A 560 , 95 (1993). [16] C. H. Lewenkopf, A. Müller and E. Doron, Phys. Rev. A, 45 , 2635 (1992). [17] X. Zheng, T. M. Antonsen and Ott, Electromagnetics, 26 ,3 (2006) ; Electromagnetics, 26 ,37 (2006). [18] S. Hemmady, et al. , Phys. Rev. Lett. 94 , 014102 (2005). [19] S. Hemmady, et al. , Phys. Rev. E 71 , 056215 (2005). [20] S. Hemmady, et al. , cond-mat/ 0512131. [21] E. Ott, “Chaos in Dynamical Systems”, Cambridge Univ. Press (1993). [22] H.-J. Stöckmann, “Quantum Chaos- An Introduction”, Cambridge Univ. Press (1999). [23] H.-J. Stöckmann and J. Stein, Phys. Rev. Lett. 64 , 2215 (1990). [24] H.-D. Gräf, et al. , Phys. Rev. Lett. 69 , 1296 (1992). [25] A. Kudrolli, et al. , Phys. Rev. E 49 , R11 (1994). [26] D. H. Wu, et al.,  Phys. Rev. Lett. 81 , 2890 (1998). [27] S.H. Chung, et al. , Phys. Rev. Lett. 85 , 2482 (2000). [28] U. Kuhl, et al. , Phys. Rev. Lett. 94 , 144101 (2005). [29] R.A. Méndez-Sánchez, et al. , Phys. Rev. Lett. 91 , 174102 (2003). [30] H. Schanze, et al. , Phys. Rev. E 71 , 016223 (2005). [31] C.W.J. Beenakker, Rev. Mod. Phys., 69 , 731 (1997). [32] Since the variance is sensitive to the presence of outliers in our data sets, which arise due to imperfect normalization or calibration errors, we exclude the outlying 1% of the data points while making our estimates of )/( 22 Qk k  n Δ .  
Similar documents
View more...
Search Related
We Need Your Support
Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks