Transcription

1 Accuracy of the coherent potential approximation for a onedimensional array with a Gaussian distribution of fluctuations in the on-site potential I. Avgin Department of Electrical and Electronics Engineering, Ege University, Bornova 35100, Izmir, Turkey A. Boukahil Physics Department, University of Wisconsin-Whitewater, Whitewater, WI D. L. Huber Physics Department, University of Wisconsin-Madison, Madison, WI We investigate the accuracy of the coherent potential approximation (CPA) for a one-dimensional array with nearest-neighbor interactions and a Gaussian distribution of fluctuations in the on-site potential. The CPA values of the integrated density of states and the inverse localization length are compared with the results of mode-counting studies carried out on arrays of sites. Good agreement is obtained suggesting that the CPA may be exact for this model. We also consider the asymptotic behavior of the inverse localization length and show that it can be approximated by the reciprocal of the decay length of a state localized about a single, strongly perturbed site in an otherwise perfect lattice. PACS: An, Aa, Dj 1

2 I. Introduction Since its introduction more than forty years ago 1,2, the coherent potential approximation, or CPA, as it is commonly known, has proven to be a successful method of determining the distribution and properties of one-electron states in disordered materials. Due to its success, efforts have been made to determine if there are models for which the CPA gives exact results. Up to this point, the CPA has proven to be exact only for the Lloyd model where there is a Cauchy (or Lorentzian) distribution of the fluctuations in the on-site potential. 3,4 Early applications of the CPA to one- dimensional arrays having a Gaussian distribution of fluctuations have shown that the approximation is unusually effective for this model as well. 5,6 Because of this effectiveness, we have undertaken a systematic study of the accuracy of the CPA for the one-dimensional, nearest-neighbor model with Gaussian disorder. As will be discussed below, we found good agreement between the CPA predictions for integrated density of states and the inverse localization length and the corresponding results obtained by mode-counting in arrays of sites. The major source for the discrepancy between CPA and mode counting appears to be errors in the numerical solution of the CPA equation, leading us to conjecture that the CPA may be exact for this model as well. As a by-product of this study, we investigate the asymptotic behavior of the inverse localization length and show that it is the characteristic of the decay rate of a wave function localized at a strongly perturbed site in an otherwise perfect array. 2

3 II. Analysis The Hamiltonian for the one-dimensional array takes the form n (1) n n H = V n >< n [ n >< n+ 1 + n >< n 1 ] where the on-site potential, V n, has a Gaussian distribution with zero mean and variance σ 2. Note that the nearest-neighbor coupling is equal to 1, so that in the absence of disorder, the band of eigenstates, symmetric about 0, ranges from 2 to 2. In the coherent potential approximation, as applied to the system characterized by CPA Eq.(1), the Green s function, G0 ( E ) is expressed as π 1 1 π c 0 CPA G0 ( E) = dk[ E + 2cos( k) V ( E)] (2) where the coherent potential, V c (E), satisfies the equation CPA dvp ( V )[ V V ( E c )]/[1 ( V V c( E ) G 0 ( E )] = 0 (3) with P(V) having the Gaussian form. The main focus of this paper is on assessing the accuracy of the coherent potential approximation for the single-particle Green s function by comparing the integral of CPA G0 ( E) with accurate numerical results obtained from large arrays. A straightforward way to do this would be to diagonalize the Hamiltonian and calculate the Green s function directly from the eigenvalues. This approach is limited, however, by number of sites in the array (typically, ). In one dimensional arrays with nearest-neighbor interactions, one can make use of mode-counting techniques to establish accuracy of the CPA in arrays of 10 8 (or more) sites. 3

4 Information about the real and imaginary parts of the Green s function can be obtained from an analysis of the sequence of amplitude ratios, R n = a n /a n-1, associated with the eigenstates of the Hamiltonian 7 R ( E) = E V 1/ R ( E) (4) n n n 1 with N being the number of sites, the number of sign changes in the sequence R 1 (E) (= E V 1 ), R 2 (E),, R N (E) corresponds to the number of modes with energies less than E. When N >> 1, which is the case here, we identify the number of sign changes, divided by N, with the integrated density of states per site, IDOS(E). The integrated density of states can also be obtained from the imaginary part of the Green s function using the equation E ( ) (0) 1 π 'Im 0( ') 0 IDOS E IDOS = de G E (5) Information about the spatial extent of the wave functions follows from a consideration of the logarithm of R n (E). 8,9 The inverse localization length, ILL(E), which is identified with the reciprocal of the average fall-off distance for eigenstates with energy E, is expressed as 8 ILL E = N < a a > 1 ( ) ln N / 1 E N N 1 1 n n N (6) n= 1 n= 1 = N ln R ( E) = N ln R ( E) The equivalent expression for ILL(E) involving the Green s function utilizes the real part of G 0 (E). It can be derived by integrating the real part of the spectral representation for G 0 (E) and making use of the expression for the inverse localization length derived in Ref. 8 4

5 ILL( E) = dxρ( x)ln E x (7) in which ρ(x) denotes the density of states. The equation takes the form E ILL( E) ILL(0) = de 'Re( G ( E ') (8) 0 0 In assessing the accuracy of the CPA, we focus on IDOS(E) IDOS(0) and ILL(E) ILL(0) whose derivatives with respect to E yield the real and imaginary parts of G 0 (E). In Figs. 1 and 2, we compare the mode-counting and CPA results for the integrals of the imaginary (Fig. 1 and Eq. (5)) and real (Fig. 2 and Eq. (8)) parts of the Green s functions calculated with σ 2 = 0.25, 1.0, and 4.0. The numerical results for the IDOS involved counting the number of sign changes in the sequence, R 1, R 2,, R N, divided by N, whereas the values for the ILL were obtained from Eq. (6). In both cases, N = The values of IDOS(0) and ILL(0) were and for σ 2 = 0.25, and for σ 2 = 1.0, and and for σ 2 = 4.0. The CPA results were CPA obtained from by solving Eqs. (2) and (3) for G0 ( E) which is then used in the evaluation of the right hand side of Eqs. (5) and (8). Standard MatLab programs in double precision were used in the solution of the self-consistent equation for the coherent potential and in the evaluation of the integrals. We note that the IDOS(E) IDOS(0) curves approach 1/2 as E, consistent with the fact that the number of modes is equal to the number of sites. In the limit N, IDOS(0) = ½ since there are equal numbers of positive and negative energy modes. The asymptotic behavior of the ILL is more interesting. In the limit of large E, ILL approaches the reciprocal of the decay length of a localized state with energy E 5

6 associated with a perturbed site in an otherwise unperturbed lattice. The connection between the energy of the localized state and the shift in the on-site potential, V 0, is 10, 11 through the equation E 2 1/2 =± 2(1 + ( V0 / 2) ) (9) where the sign corresponds to the sign of V 0. The reciprocal of the decay length, δ (E), (in units of the inverse lattice constant) associated with the localized state is given by 2 1/2 1/ ( E) ln{( E / 2 )[1 (1 4 / E ) ]} δ = (10) which is appropriate only for E > 2 and has the limiting behavior 1/δ (E) ln E for large E. In Fig. 3 we compare the mode-counting values of the ILL(E) with the results obtained by approximating ILL(E) as 1/δ (E), which we refer to as the single-site approximation. In contrast to IDOS(0), there appears to be no simple analytical expression for ILL(0). We postpone consideration of these results until the following section. III. Discussion It is evident from Figs. 1 and 2 that there is good agreement between the numerical results and the results obtained using the coherent potential approximation. The extent of the agreement for σ 2 = 1.0 is shown in greater detail in Fig. 4 where we plot the relative difference ΔIDOS = [IDOS(E) mode-counting IDOS(E) CPA ]/IDOS(E) mode-counting and in Fig. 5 where we plot ΔILL = [ILL(E) mode-counting ILL(E) CPA ]/ILL(E) mode-counting 6

7 where in both cases N = In Figs. 4 and 5 it is apparent that ΔIDOS(E) < whereas ΔILL(E) < Since we have carried out similar computations with a different sequence of potential fluctuations and obtained results that were indistinguishable from those shown, we suspect that the qualitative difference between the IDOS and the ILL is connected with numerical solution of the self-consistent equation for the coherent potential. The results displayed in Fig. 3, along with Eq. (9), indicate that the single-site approximation works well when the energy of the localized state lies in a region where 2 the probability of the corresponding potential fluctuation V0 =± 2[( E / 2) 1], which is expressed as PV = E E (11) 2 1/ /2 2 2 ( 0) (2 πσ ) exp{ 2[( / 2) 1] / σ } (2 πσ ) exp[ / 2 σ ] is extremely small, i.e. P(V 0 )/P(0) exp[ E 2 /2σ 2 ] << 1. It should be noted that the single-site approximation is not limited to the Gaussian distribution, but is applicable to other distributions in the large-e limit, e.g. the Cauchy distribution, as can be seen from the exact expression for the ILL of the Cauchy distribution given in Ref. 8. Taken together, the results presented in this paper are evidence that the CPA is a good approximate theory for the single-particle Green s function for the one-dimensional Gaussian model and lead us to conjecture that it may be exact. To strengthen the argument, it is necessary to make significant improvements in the accuracy of the numerical solutions of the CPA equations, although only analytical studies can establish that the coherent potential approximation is exact. 7

Chapter 4 Genomic correlations and sequence alignment statistics Long-range correlations in genomic base composition are a ubiquitous statistical feature among many eukaryotic genomes. As we have shown

The Quantum Harmonic Oscillator Stephen Webb The Importance of the Harmonic Oscillator The quantum harmonic oscillator holds a unique importance in quantum mechanics, as it is both one of the few problems

Feshbach molecules in a one-dimensional Fermi gas Abstract We consider the binding energy and the wave function of Feshbach molecules confined in a one-dimensional matter waveguide. We compare the binding

II.D Many Random Variables With more than one random variable, the set of outcomes is an -dimensional space, S x = { < x 1, x,, x < }. For example, describing the location and velocity of a gas particle

Transport properties of a ladder with two random dimer chains Hu Dong-Sheng( ), Zhu Chen-Ping( ), and Zhang Yong-Mei( ) College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing 211100,

Continuous Random Variables The probability that a continuous random variable, X, has a value between a and b is computed by integrating its probability density function (p.d.f.) over the interval [a,b]:

Teoretisk Fysik Hand-in problem B, SI1142, Spring 2010 In 1955 Fermi, Pasta and Ulam 1 numerically studied a simple model for a one dimensional chain of non-linear oscillators to see how the energy distribution

Chapter 1 Thomas-Fermi Theory The Thomas-Fermi theory provides a functional form for the kinetic energy of a non-interacting electron gas in some known external potential V (r) (usually due to impurities)

Eigenvalues, Eigenvectors, and Differential Equations William Cherry April 009 (with a typo correction in November 05) The concepts of eigenvalue and eigenvector occur throughout advanced mathematics They

Landau Theory o Second Order Phase Transitions Order Parameter Second order phase transitions occur when a new state o reduced symmetry develops continuously rom the disordered (high temperature) phase.

A MODEL OF PULSED EDDY CURRENT CRACK DETECTION R.E. Beissner and J.L. Fisher Southwest Research Institute San Antonio, Texas 78284 INTRODUCTION As its name implies, the pulsed eddy current (PEC) method

Chapter 4 The Landau-Ginzburg Free Energy as a Hamiltonian c 2012 by W. Klein, Harvey Gould, and Jan Tobochnik 1 October 2012 4.1 Introduction In this chapter we will follow the paradigm discussed in Chapters

A2.PS.1 Use a variety of problem solving strategies to understand new mathematical content A2.PS.2 Recognize and understand equivalent representations of a problem situation or a mathematical concept A2.PS.3

Review of Fundamental Mathematics As explained in the Preface and in Chapter 1 of your textbook, managerial economics applies microeconomic theory to business decision making. The decision-making tools

PHY464 Introduction to Quantum Mechanics Fall 4 Practice Test 3 November, 4 These problems are similar but not identical to the actual test. One or two parts will actually show up.. Short answer. (a) Recall

for all values of x, greater than 1 and less than 1 Andrew York May 7, 2016 In my previous article "Throwing a Curve at Logarithmic Spirals" I discussed a new way to generate logarithmic or equiangular

RANDOM VIBRATION AN OVERVIEW by Barry Controls, Hopkinton, MA ABSTRACT Random vibration is becoming increasingly recognized as the most realistic method of simulating the dynamic environment of military

Chapter 4 Tunneling and the semiclassical propagator 4.1 A quantum-semiclassical approach 4.1.1 Motivation In the foregoing Chapter 3, tunneling was described in the energy domain by means of classical

Improvement of a System for Prime Factorization Based on Optical Interferometer Kouichi Nitta, Nobuto Katsuta, and Osamu Matoba Department of Computer Science and Systems Engineering, Graduate School of

FYSA21 Mathematical Tools in Science Signal analysis in astronomy Lennart Lindegren Department of Astronomy and Theoretical Physics, Lund University 1 Aim of the project This project will illustrate an

CALCULATIONS & STATISTICS CALCULATION OF SCORES Conversion of 1-5 scale to 0-100 scores When you look at your report, you will notice that the scores are reported on a 0-100 scale, even though respondents

The ensity of States When describing very large systems in which an exact treatment of all particles is impractical or undesirable, we often statistically treat the influence of these many degrees of freedom,

Numerical methods for Ordinary Differential Equations Didier Gonze November 18, 2013 I have come to believe that ones knowledge of any dynamical system is deficient unless one knows a valid way to numerically

Chapter 4 Green s function formalism Having a numerical representation of the system in terms of a tight-binding model, one still needs a mathematical framework within which different physical properties

Time Independent Perturbation Theory Contd. Keep in mind that if you can come up with arguments that make the matrix element is zero, it will make the calculation easier as we had in the HO example. If