DOCSLIB.ORG

T-Matrix Theories of the Attractive Hubbard Model

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
T-Matrix Theories of the Attractive Hubbard Model Diagrammatic Analysis of T-matrix Theories of the Attractive Hubbard Model Iiwin S. D.Beach .A r tiesis submitted to the Department of Physics in conformity wit 11 the rccpirernents for the degree of Master of Science Queen's University Iiingstori. Ontario. Canada September. 1999 @ Kevin S. D.Beach. 1999 National Library Bibliotheque nationale I*I ofCanada du Canada Acquisitions and Acquisitions et Bibliographic Services services bibliographiques 395 Wellmgton Street 395, rue Wellington Ottawa ON K 1A ON4 Ottawa ON K 1A ON4 Canada Canada Your h~eVow reference Our fm Mtre reldwnce The author has granted a non- L'auteur a accorde une licence non exclusive licence allowing the exclusive permettant a la National Library of Canada to Bibliotheque nationale du Canada de reproduce, loan, distribute or sell reproduire, preter, distribuer ou copies of this thesis in microform, vendre des copies de cette these sous paper or electronic formats. la forme de microfiche/filrn. de reproduction sur papier ou sur format electronique. The author retains ownership of the L'auteur conserve la propriete du copyright in this thesis. Neither the droit d'auteur qui protege cette these. thesis nor substantial extracts from it Ni la these ni des extraits substantiels may be printed or otherwise de celle-ci ne doivent ttre imprimes reproduced without the author's ou autrement reproduits sans son permission. autorisation. Abstract 111 r his uork, r~~~iny-bocfre(-hniqu~s are used to st.udy rho effwt of ~lectroripair c-orr~latious on th(1 nor111;d~tate prop~rti~s of th~t~vo-clirn~nsi~naI attractiv~ Hubbarci ntociel npar t tlr c.ritical tt~r~iperstrirc( T.) of the sr~p~rcon~iuctingtransition. The role clt' .;elf-t.onsistrn(*y ir~t hr st ,zndxr[i T-matrix npprosirnation is pxarninrd ivi t h particular at t~ntionpaid to its c.orltriLution to the suppression of the r-ritical temprrature and to the depletion of spertrai 1i.cxiqht. Ind~ed.ivt) show that the ~~r~forc.crnrntof .;elf-consistency ~irarnatic-allyxlttlrs t hcl valu~oS I'. 1 from its r~on-self-(-onsistont value) and that t he choice of self-c-onsisttwc.y~h~tw cmrnpl~telydetprrriines the predicted density of states at T,,. Mmovcr, at low tenlpt.rnturw Iitaitr T-, non-self-consistent cornpu tational results suggest that the depletion of ~pc!c*trid ivvight at the fern~ilevel is prunounced ard leads to the opening of a pseudoqap. a prwursor to t he cl!icr%y gap of the superconducting state. In xdditior~,we critique L Pad6 analytic continucztion met hod whereby a rational poly- nornial function is fit to a set of' input points by means of a single matrix inversion. The r~let,llo(lis applied to the problerri of determining the spectral function of a one-particle propagator known only st a finite number of llstsubara frequencies. This is dono for prop- agators constructed from two example self-energies, drawn from the T-matrix theory of the attractive Hubbard model. which are calculated to extremely high accuracy using a novel symbolic corn putation algorithm. Cl'e present a systematic analysis of the effects of error in the input points and discuss how the convergence of the polynomial coefficients to allowed values can serve as a quantitative measure of the goodness-of-fit. Statement of originality Acknowledgements Ihis tvork ivas carried out 11ncic.rthe supervision of Dr. Robert d. Gooding (Qsrtri'.q I'rlirrr- .i,t// ) to LV horn the a11t hor o~~sa great ~lebt of gratitude for his guidance and cncouragernrrit. I'hx~ksare illdo doe to Dr. Frank llnrsiglio f-nr~.rrsit!jof .-1lb~rtnlfor many helpful tiis- ~-l~ssions.rnd for qrariously pro\*i(iing nnm~ri(.illrwdts ~vhichivr.ere r~quir~dto ti~bug.in11 vrrify the computer code used in this project. Fillally, the author wishes to nckno~vlniqe the unconditional support. of his family and of his partner, Diana Drappel. Notational convent ions .\ v~ctorcluarltity is (lenoted by an arrow arrd its c:ornponents by supt\rsc.ript..;: rb.g. t.he s!*rrrbol i: = ( rl,rl.. uD) represents a vector of dimension D. :In operator in rht: ~wcupationr~unibcfr forr~idism is marked by a, caret. except for the fernlion wrni hilation ;L% possible). P.E. ti = cfv. .-Ijrrmronrr operator consists of an odd number uf c nnql r."i ivhortw .a br).wr~tcoperator consists of an even nurnher. - The cornrnutntor of two operators .iand b is given by [.-i.B] = .4B - H.4 and rhc ~lr~trt.or~lrrl~rtntorby {.-i.b} = ..id + 6.i. \Ye shall use brackets wirh a ii%a iut~scripr. to ~l~nut,cc'ittl~r the commutator or anticornnlut.ator according to [..i.B]* = .ib + fj.-i. -4 useful consequence of this nition is the identity = - .41 def [.-iB,i'] .-i[d. i'] [i'. _ 8. ivhich holds for both kinds of statistics. Generally. whenever we treat fernli and hose c;tses jimdtnneously using + and 7. the fernii case will be represented by the upper sign. The herrnitian conjugate (or adjoint) is denoted by a dagger (t) and the con~plexcon- jugate by an asterisk (*). An operator .iis called hermitinn (or self-adjoint) if .-it = .Lie Every matrix element of s herrnitian operator satisfies Sotice that if we choose lo) = 10') to be an eigenstate of .i.then the equation above tells 11s that the eigenvalues of a hermitian operator are real. rl is inlplicitlv cv;rhlat.ed in the limit rl -. 0'. l\'e occasionally employ the shorthand -r, - r + 11. In t.ern~sof the infinitesimal rl. the Heaviside and Dirac t'unctions can be represented by set union set intersection set subtraction csrtesian product direct sum absolute value norm NSC' non-self-consistent pp 1 partic~e-prticI~ I I particle-hole (JhIC' quantum 1Iorltc C'arlo TCL Thouless criterion line Table of contents Abstract Statement of originality Acknowledgements Notational conventions 1v ... Table of contents Vlll List of figures 1 Attractive Hubbard model 1.1 [ntroduction . I Uoclel harniltonian . L .?.I Real space reprcsentrtt,ion . 1 \C',zvevector representation . 1 Spin representation . , . 1. Exact solutions . 1.-I Thermody narnics of the Hubbard model . l .j S~mmetriesof the haniiltonian . - . 1. Limiting cases . - . L.6.1 Band limit . , . 1.6 Atomic limit . 1.6.3 >on-trivial values of IV/t. , , . 2 Formalism of diagrammatic field theory 2.1 Approaching the strongly correlated electron problem . 2. Functional differentiation . 2. Equation of motion method . 2.4 Fourier representation . 2.5 T-matrix approximation . 2.6 Transition to the superconducting state . 3 rllt~alysisof the T-matrix approximation 67 1 O\.~rviciv ..................................... ci: -I A new approach to analytic continuation 115 1. 1 Ovt1rvil\w ..................................... 11 .: 1.2 Green's t'unctiori forrnalisrn ........................... llti 1.3 Pati6 approxinlmts ................................ 118 !. Sunlerical rcsdts ................................. 123 t.5 Summary ..................................... 131 5 Conclusions and future work 13-1 A Basic formalism 137 .\. 1 Occupation number formalism (second quantization) ............. 1;jT A.2 Real space representation ............................ LXl .4 .3 It'avevector representation ............................ 142 :\..1 Spin operator ................................... 1-15 ..I..? Fluxoperator ................................... 143 .A .rj Local number conservation ........................... 1-52 .A .; Density of levels ................................. 1.34 B Many-body formalism 161 B.1 Green's function formalism ........................... 161 B . 1.1 Time ordering operator ......................... 161 8.1.2 Equations of motion ........................... 162 8.1.3 Green's functions ............................. 163 B . 1.4 Fourier series for temperature functions ................ 166 8.1..5 Fourier respresentations ......................... 167 B.1.6 Analytic continuation .......................... 171 B.1.7 1Iornents of the spectral function .................... 172 B.2 Occupation functions ............................... 171 B.2.1 Basic occupation functions ....................... 174 B.2.2 Extended occupation functions ..................... 179 8.2.3 Partial occupation functions ....................... 132 C Functional differentiation 193 '1 F111lct.ionaltiitfcr~ntiiitiorl . 193 ('.2 5kclcto11 clspanslorl . I07 D Equation of motion method 300 D.1 Lirlt!ar t'unc.t.ionnls . 200 2 Evaluating t tie corllrnutntor ul the ~ranti~~~nonical harniltoniarl ~vit11 all ;in- !lihilation operator . 204 D.3 Equation of motiur~. 205 L). I Terrilinating the hierarch! 01' equations . .)()*- D.5 C;cneraiized intt>ractiori . 2 I I E h1APLE code 215 E.l Pndr. package . 215 I:.? €luhbi1rd2ci package . 2 l(i Vitae 254 List of figures [ortic xr~erlirlgof r.1t.c trorls ........................... .) (I;eneric loiv.tertlperaturt.i plla.~~liagrnrn of t.he h~gh.1.. I-ripr:itcs ........1 Levrl ~rrri'acwof the dispersion relation in tile Brillouin zone ........ '1 ~lomcnturr~transfer in a singlcl dectrorl-electron x;lttcrinq ovetlt....... 10 Li~lerrysp~ctrurn for ;in ~lcctrnn pair ...................... 1ri t3h;tscl ~iiagrarn(A' tllc t'rcle c:lectron ga ..................... -'q- Density of levtals it\ une . t.tvo. and three ltirn~nsior\s .............. :10 Phase tliagranl of the attractive Hubbard nlmlel in the atomic limit .... 3-4 Doubie occupancy in the atomic limit ..................... 3.5 2.1 Sktlleturl vspa~~sionof t.he Hubbard sdfen~rg,y ................ 1-1 2.2 Skt.l~torltlupansion of t.he opposit~-.; pin t~vo-particlerorrelator ....... li 2.3 Skeleton expansion of the equal-spin two-particle corrrlator ......... lfj 2 . I Dccon~positionof
Load more

Recommended publications
  • The Digamma Function and Explicit Permutations of the Alternating Harmonic Series
    The digamma function and explicit permutations of the alternating harmonic series. Maxim Gilula February 20, 2015 Abstract The main goal is to present a countable family of permutations of the natural numbers that provide explicit rearrangements of the alternating harmonic series and that we can easily define by some closed expression. The digamma function presents its ubiquity in mathematics once more by being the key tool in computing explicitly the simple rearrangements presented in this paper. The permutations are simple in the sense that composing one with itself will give the identity. We show that the count- able set of rearrangements presented are dense in the reals. Then, slight generalizations are presented. Finally, we reprove a result given originally by J.H. Smith in 1975 that for any conditionally convergent real series guarantees permutations of infinite cycle type give all rearrangements of the series [4]. This result provides a refinement of the well known theorem by Riemann (see e.g. Rudin [3] Theorem 3.54). 1 Introduction A permutation of order n of a conditionally convergent series is a bijection φ of the positive integers N with the property that φn = φ ◦ · · · ◦ φ is the identity on N and n is the least such. Given a conditionally convergent series, a nat- ural question to ask is whether for any real number L there is a permutation of order 2 (or n > 1) such that the rearrangement induced by the permutation equals L. This turns out to be an easy corollary of [4], and is reproved below with elementary methods. Other \simple"rearrangements have been considered elsewhere, such as in Stout [5] and the comprehensive references therein.
    [Show full text]
  • The Riemann and Hurwitz Zeta Functions, Apery's Constant and New
    The Riemann and Hurwitz zeta functions, Apery’s constant and new rational series representations involving ζ(2k) Cezar Lupu1 1Department of Mathematics University of Pittsburgh Pittsburgh, PA, USA Algebra, Combinatorics and Geometry Graduate Student Research Seminar, February 2, 2017, Pittsburgh, PA A quick overview of the Riemann zeta function. The Riemann zeta function is defined by 1 X 1 ζ(s) = ; Re s > 1: ns n=1 Originally, Riemann zeta function was defined for real arguments. Also, Euler found another formula which relates the Riemann zeta function with prime numbrs, namely Y 1 ζ(s) = ; 1 p 1 − ps where p runs through all primes p = 2; 3; 5;:::. A quick overview of the Riemann zeta function. Moreover, Riemann proved that the following ζ(s) satisfies the following integral representation formula: 1 Z 1 us−1 ζ(s) = u du; Re s > 1; Γ(s) 0 e − 1 Z 1 where Γ(s) = ts−1e−t dt, Re s > 0 is the Euler gamma 0 function. Also, another important fact is that one can extend ζ(s) from Re s > 1 to Re s > 0. By an easy computation one has 1 X 1 (1 − 21−s )ζ(s) = (−1)n−1 ; ns n=1 and therefore we have A quick overview of the Riemann function. 1 1 X 1 ζ(s) = (−1)n−1 ; Re s > 0; s 6= 1: 1 − 21−s ns n=1 It is well-known that ζ is analytic and it has an analytic continuation at s = 1. At s = 1 it has a simple pole with residue 1.
    [Show full text]
  • COMPLETE MONOTONICITY for a NEW RATIO of FINITE MANY GAMMA FUNCTIONS Feng Qi
    COMPLETE MONOTONICITY FOR A NEW RATIO OF FINITE MANY GAMMA FUNCTIONS Feng Qi To cite this version: Feng Qi. COMPLETE MONOTONICITY FOR A NEW RATIO OF FINITE MANY GAMMA FUNCTIONS. 2020. hal-02511909 HAL Id: hal-02511909 https://hal.archives-ouvertes.fr/hal-02511909 Preprint submitted on 19 Mar 2020 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. COMPLETE MONOTONICITY FOR A NEW RATIO OF FINITE MANY GAMMA FUNCTIONS FENG QI Dedicated to people facing and fighting COVID-19 Abstract. In the paper, by deriving an inequality involving the generating function of the Bernoulli numbers, the author introduces a new ratio of finite many gamma functions, finds complete monotonicity of the second logarithmic derivative of the ratio, and simply reviews complete monotonicity of several linear combinations of finite many digamma or trigamma functions. Contents 1. Preliminaries and motivations 1 2. A lemma 3 3. Complete monotonicity 4 4. A simple review 5 References 7 1. Preliminaries and motivations Let f(x) be an infinite differentiable function on (0; 1). If (−1)kf (k)(x) ≥ 0 for all k ≥ 0 and x 2 (0; 1), then we call f(x) a completely monotonic function on (0; 1).
    [Show full text]
  • Euler-Maclaurin and Euler-Boole Formulas
    Appendix A Euler-MacLaurin and Euler-Boole Formulas A.1 A Taylor Formula The classical Taylor formula Z Xm xk x .x t/m f .x/ D @kf .0/ C @mC1f .t/dt kŠ 0 mŠ kD0 xk can be generalized if we replace the polynomial kŠ by other polynomials (Viskov 1988; Bourbaki 1959). Definition If is a linear form on C0.R/ such that .1/ D 1,wedefinethe polynomials .Pn/ by: P0 D 1 @Pn D Pn1 , .Pn/ D 0 for n 1 P . / k The Generating Function k0 Pk x z We have formally X X X k k k @x. Pk.x/z / D Pk1.x/z D z. Pk.x/z / k0 k1 k0 © Springer International Publishing AG 2017 175 B. Candelpergher, Ramanujan Summation of Divergent Series, Lecture Notes in Mathematics 2185, DOI 10.1007/978-3-319-63630-6 176 A Euler-MacLaurin and Euler-Boole Formulas thus X k xz Pk.x/z D C.z/e k0 To evaluate C.z/ we use the notation x for and by definition of .Pn/ we can write X X k k x. Pk.x/z / D x.Pk.x//z D 1 k0 k0 X k xz xz x. Pk.x/z / D x.C.z/e / D C.z/x.e / k0 1 this gives C.z/ D xz . Thus the generating function of the sequence .Pn/ is x.e / X n xz Pn.x/z D e =M.z/ n xz where the function M is defined by M.z/ D x.e /: Examples P xn n xz .
    [Show full text]
  • On Some Series Representations of the Hurwitz Zeta Function Mark W
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector Journal of Computational and Applied Mathematics 216 (2008) 297–305 www.elsevier.com/locate/cam On some series representations of the Hurwitz zeta function Mark W. Coffey Department of Physics, Colorado School of Mines, Golden, CO 80401, USA Received 21 November 2006; received in revised form 3 May 2007 Abstract A variety of infinite series representations for the Hurwitz zeta function are obtained. Particular cases recover known results, while others are new. Specialization of the series representations apply to the Riemann zeta function, leading to additional results. The method is briefly extended to the Lerch zeta function. Most of the series representations exhibit fast convergence, making them attractive for the computation of special functions and fundamental constants. © 2007 Elsevier B.V. All rights reserved. MSC: 11M06; 11M35; 33B15 Keywords: Hurwitz zeta function; Riemann zeta function; Polygamma function; Lerch zeta function; Series representation; Integral representation; Generalized harmonic numbers 1. Introduction (s, a)= ∞ (n+a)−s s> a> The Hurwitz zeta function, defined by n=0 for Re 1 and Re 0, extends to a meromorphic function in the entire complex s-plane. This analytic continuation to C has a simple pole of residue one. This is reflected in the Laurent expansion ∞ n 1 (−1) n (s, a) = + n(a)(s − 1) , (1) s − 1 n! n=0 (a) (a)=−(a) =/ wherein k are designated the Stieltjes constants [3,4,9,13,18,20] and 0 , where is the digamma a a= 1 function.
    [Show full text]
  • A New Entire Factorial Function
    A new entire factorial function Matthew D. Klimek Laboratory for Elementary Particle Physics, Cornell University, Ithaca, NY, 14853, USA Department of Physics, Korea University, Seoul 02841, Republic of Korea Abstract We introduce a new factorial function which agrees with the usual Euler gamma function at both the positive integers and at all half-integers, but which is also entire. We describe the basic features of this function. The canonical extension of the factorial to the complex plane is given by Euler's gamma function Γ(z). It is the only such extension that satisfies the recurrence relation zΓ(z) = Γ(z + 1) and is logarithmically convex for all positive real z [1, 2]. (For an alternative function theoretic characterization of Γ(z), see [3].) As a consequence of the recurrence relation, Γ(z) is a meromorphic function when continued onto − the complex plane. It has poles at all non-positive integer values of z 2 Z0 . However, dropping the two conditions above, there are infinitely many other extensions of the factorial that could be constructed. Perhaps the second best known factorial function after Euler's is that of Hadamard [4] which can be expressed as sin πz z z + 1 H(z) = Γ(z) 1 + − ≡ Γ(z)[1 + Q(z)] (1) 2π 2 2 where (z) is the digamma function d Γ0(z) (z) = log Γ(z) = : (2) dz Γ(z) The second term in brackets is given the name Q(z) for later convenience. We see that the Hadamard gamma is a certain multiplicative modification to the Euler gamma.
    [Show full text]
  • Transcendence of Various Infinite Series Applications of Baker's
    Transcendence of Various Infinite Series and Applications of Baker's Theorem by Chester Jay Weatherby A thesis submitted to the Department of Mathematics and Statistics in conformity with the requirements for the degree of Doctor of Philosophy Queen's University Kingston, Ontario, Canada April 2009 Copyright c Chester Jay Weatherby, 2009 Library and Archives Bibliothèque et Canada Archives Canada Published Heritage Direction du Branch Patrimoine de l’édition 395 Wellington Street 395, rue Wellington Ottawa ON K1A 0N4 Ottawa ON K1A 0N4 Canada Canada Your file Votre référence ISBN: 978-0-494-65340-1 Our file Notre référence ISBN: 978-0-494-65340-1 NOTICE: AVIS: The author has granted a non- L’auteur a accordé une licence non exclusive exclusive license allowing Library and permettant à la Bibliothèque et Archives Archives Canada to reproduce, Canada de reproduire, publier, archiver, publish, archive, preserve, conserve, sauvegarder, conserver, transmettre au public communicate to the public by par télécommunication ou par l’Internet, prêter, telecommunication or on the Internet, distribuer et vendre des thèses partout dans le loan, distribute and sell theses monde, à des fins commerciales ou autres, sur worldwide, for commercial or non- support microforme, papier, électronique et/ou commercial purposes, in microform, autres formats. paper, electronic and/or any other formats. The author retains copyright L’auteur conserve la propriété du droit d’auteur ownership and moral rights in this et des droits moraux qui protège cette thèse. Ni thesis. Neither the thesis nor la thèse ni des extraits substantiels de celle-ci substantial extracts from it may be ne doivent être imprimés ou autrement printed or otherwise reproduced reproduits sans son autorisation.
    [Show full text]
  • POCKETBOOKOF MATHEMATICAL FUNCTIONS Abridged Edition of Handbook of Mathematical Functions Milton Abramowitz and Irene A
    POCKETBOOKOF MATHEMATICAL FUNCTIONS Abridged edition of Handbook of Mathematical Functions Milton Abramowitz and Irene A. Stegun (eds.) Material selected by Michael Danos and Johann Rafelski 1984 Verlag Harri Deutsch - Thun - Frankfurt/Main CONTENTS Forewordtothe Original NBS Handbook 5 Pref ace 6 2. PHYSICAL CONSTANTS AND CONVERSION FACTORS 17 A.G. McNish, revised by the editors Table 2.1. Common Units and Conversion Factors 17 Table 2.2. Names and Conversion Factors for Electric and Magnetic Units 17 Table 2.3. AdjustedValuesof Constants 18 Table 2.4. Miscellaneous Conversion Factors 19 Table 2.5. FactorsforConvertingCustomaryU.S. UnitstoSIUnits 19 Table 2.6. Geodetic Constants 19 Table 2.7. Physical andNumericalConstants 20 Table 2.8. Periodic Table of the Elements 21 Table 2.9. Electromagnetic Relations 22 Table 2.10. Radioactivity and Radiation Protection 22 3. ELEMENTARY ANALYTICAL METHODS 23 Milton Abramowitz 3.1. Binomial Theorem and Binomial Coefficients; Arithmetic and Geometrie Progressions; Arithmetic, Geometrie, Harmonie and Generalized Means 23 3.2. Inequalities 23 3.3. Rules for Differentiation and Integration 24 3.4. Limits, Maxima and Minima 26 3.5. Absolute and Relative Errors 27 3.6. Infinite Series 27 3.7. Complex Numbers and Functions 29 3.8. Algebraic Equations 30 3.9. Successive Approximation Methods 31 3.10. TheoremsonContinuedFractions 32 4. ELEMENTARY TRANSCENDENTAL FUNCTIONS 33 Logarithmic, Exponential, Circular and Hyperbolic Functions Ruth Zucker 4.1. Logarithmic Function 33 4.2. Exponential Function 35 4.3. Circular Functions 37 4.4. Inverse Circular Functions 45 4.5. Hyperbolic Functions 49 4.6. InverseHyberbolicFunctions 52 5. EXPONENTIAL INTEGRAL AND RELATED FUNCTIONS .... 56 Walter Gautschi and William F.
    [Show full text]
  • A Note on the Zeros and Local Extrema of Digamma Related Functions
    A note on the zeros and local extrema of Digamma related functions Istv´an Mez˝o 1 School of Mathematics and Statistics, Nanjing University of Information Science and Technology, No.219 Ningliu Rd, Pukou, Nanjing, Jiangsu, P. R. China Abstract Little is known about the zeros of the Digamma function. Establishing some Weier- strassian infinite product representations for a given regularization of the Digamma function we find interesting sums of its zeros. In addition, we study the same ques- tions for the zeros of the logarithmic derivative of the Barnes G function. At the end of the paper we provide rather accurate approximations of the hyperfactorial where, rather interestingly, the Lambert function appears. Key words: Euler Gamma function; Digamma function; Barnes G function; zeros; Weierstrass Product Theorem 1991 MSC: 33B15, 30C15, 30D99 arXiv:1409.2971v2 [math.CV] 9 Feb 2016 1 Introduction The Euler gamma function is defined by the improper integral ∞ Γ(z)= e−ztz−1dt (ℜ(z) > 0). Z0 The logarithmic derivative of Γ is the Digamma function Γ′(z) ψ(z)= (z ∈ C \{0, −1, −2,... }). Γ(z) Email address: [email protected] (Istv´an Mez˝o). 1 The research of Istv´an Mez˝owas supported by the Scientific Research Foundation of Nanjing University of Information Science & Technology, the Startup Foundation for Introducing Talent of NUIST. Project no.: S8113062001, and the National Natural Science Foundation for China. Grant no. 11501299. Preprint submitted to 10 February 2016 This function is analytic everywhere except the non positive integers, and it has first order poles in these points. It is known that the Digamma function has only real and simple zeros, and all of them except only one are negative.
    [Show full text]
  • Clustering, Features Selection, and Outliers Rejection
    Statistical Modeling for Simultaneous Data Clustering, Features Selection, and Outliers Rejection Khaled Almakadmeh A Thesis in The Concordia Institute for Information Systems Engineering Presented in Partial Fulfillment of the Requirements for the Degree of Master of Applied Science in Information Systems Security at Concordia University Montréal, Québec, Canada May 2010 © Khaled Almakadmeh, 2010 Library and Archives Bibliothèque et ?F? Canada Archives Canada Published Heritage Direction du Branch Patrimoine de l'édition 395 Wellington Street 395, rue Wellington Ottawa ON K1A 0N4 OttawaONK1A0N4 Canada Canada Your We Votre référence ISBN: 978-0-494-71008-1 Our file Notre référence ISBN: 978-0-494-71008-1 NOTICE: AVIS: The author has granted a non- L'auteur a accordé une licence non exclusive exclusive license allowing Library and permettant à la Bibliothèque et Archives Archives Canada to reproduce, Canada de reproduire, publier, archiver, publish, archive, preserve, conserve, sauvegarder, conserver, transmettre au public communicate to the public by par télécommunication ou par l'Internet, prêter, telecommunication or on the Internet, distribuer et vendre des thèses partout dans le loan, distribute and sell theses monde, à des fins commerciales ou autres, sur worldwide, for commercial or non- support microforme, papier, électronique et/ou commercial purposes, in microform, autres formats. paper, electronic and/or any other formats. The author retains copyright L'auteur conserve la propriété du droit d'auteur ownership and moral rights in this et des droits moraux qui protège cette thèse. Ni thesis. Neither the thesis nor la thèse ni des extraits substantiels de celle-ci substantial extracts from it may be ne doivent être imprimés ou autrement printed or otherwise reproduced reproduits sans son autorisation.
    [Show full text]
  • Some Inequalities for the Trigamma Function in Terms of the Digamma
    SOME INEQUALITIES FOR THE TRIGAMMA FUNCTION IN TERMS OF THE DIGAMMA FUNCTION FENG QI AND CRISTINEL MORTICI Abstract. In the paper, the authors establish three kinds of double inequal- ities for the trigamma function in terms of the exponential function to powers of the digamma function. These newly established inequalities extend some known results. The method in the paper utilizes some facts from the asymp- totic theory and is a natural way to solve problems for approximating some quantities for large values of the variable. 1. Motivations and main results The classical Euler gamma function may be defined for (z) > 0 by ∞ ℜ − − Γ(z)= tz 1e t d t. Z0 The logarithmic derivative of the gamma function Γ(z) is denoted by d Γ′(z) ψ(z)= [ln Γ(z)] = d z Γ(z) and called the digamma function. The derivatives ψ′(z) and ψ′′(z) are called the trigamma and tetragamma functions respectively. As a whole, the functions ψ(k)(z) for k 0 N are called the polygamma functions. These functions are widely used in theoretical∈{ }∪ and practical problems in all branches of mathematical science. Con- sequently, many mathematicians were preoccupied to establish new results about the gamma function, polygamma functions, and other related functions. 1.1. The first main result. In 2007, Alzer and Batir [3, Corollary] discovered that the double inequality 1 1 √2π xx exp x ψ(x + α) < Γ(x) < √2π xx exp x ψ(x + β) (1.1) arXiv:1503.03020v1 [math.CA] 7 Mar 2015 − − 2 − − 2 1 holds for x> 0 if and only if α 3 and β 0.
    [Show full text]
  • Arxiv:1705.06547V1 [Math.CA] 18 May 2017 Hysatisfy They Am Ucin Oyam Ucin,Inequalities
    INEQUALITIES FOR THE INVERSES OF THE POLYGAMMA FUNCTIONS NECDET BATIR Abstract. We provide an elementary proof of the left side in- equality and improve the right inequality in 1 n+1 n! − − < ((−1)n 1ψ(n)) 1(x) x − (x−1/n + α)−n 1 n! n+1 < , x − (x−1/n + β)−n where α = [(n − 1)!]−1/n and β = [n!ζ(n + 1)]−1/n, which was proved in [6], and we prove the following inequalities for the inverse of the digamma function ψ. 1 − 1 < ψ 1(x) <ex + , x ∈ R. log(1 + e−x) 2 The proofs are based on nice applications of the mean value theo- rem for differentiation and elementary properties of the polygamma functions. 1. introduction As it is well known for a positive real number x the gamma function Γ is defined to be ∞ Γ(x)= ux−1e−udu. Z0 It is a common knowledge that the gamma function plays a special role arXiv:1705.06547v1 [math.CA] 18 May 2017 in the theory of special functions. The most important function related to the gamma function is the digamma or psi function ψ(x), which is defined by logarithmic derivative of the gamma function Γ(x), that is, Γ′(x) ψ(x)= . Γ(x) The digamma function ψ is closely related with the Euler-Mascheroni 1 ··· 1 constant γ (= 0.57721...) and harmonic numbers Hn =1+ 2 + + n . They satisfy ψ(n +1) = −γ + Hn (n ∈ N). In [14] it is proved that 2000 Mathematics Subject Classification. Primary: 33B15; Secondary: 26D07.
    [Show full text]