The discovery of HTS [1], excitingly opened up a new field of research in solid state physics experimentalists as well as theorists. The effect of anharmonicity, impurities as electron-phonon problem in HTS is least studied [2-5] due to its very complicated nature, because the unit cell of these compounds contain, beside two superconducting planes, the chain elements connected by bridge with plain through the apical oxygen ions . The interaction of electrons with anharmonic lattice vibrations is a longstanding problem that is not yet fully understood. The effects of anharmonicity on the electron-phonon problem and more specifically on superconductivity are relatively unknown. The present investigation deals with the impurity induced anharmonic phonon-electron problem, in which the contribution due to anharmonicities, isotropic impurities and interference has been dealt via an almost complete Hamiltonian. Having developed the electron Green’s function the expressions for energy spectrum and electron density of states have been obtained with few new features for the High temperature superconductors.

Let us consider the double-time thermodynamic electron retarded Green’s function [6,7]

via an almost complete Hamiltonian [5-7]

where, , , and, respectively are unperturbed electron-, unperturbed phonon-, electronphonon-, anharmonic (upto quartic terms)-, and defect contributions to the Hamiltonian and are expressible in the form [8]

In the above equations and are the electron creation (annihilation) and phonon field and momentum operators, respectively. (and are phonon and electron wave vectors) and stands for electron-phonon coupling coefficient., and are anharmonic coupling coefficients mass and force constant change parameters, respectively. Following the equation of motion technique [8-10] of quantum dynamics via Hamiltonian (2) and using Dyson’s equation approach we can obtain the Green’s function as

where and are electron and pairon frequencies.

The excitation spectrum i.e. a response function can be expressed as

Here and are the electron energy shift and electron line width respectively. Higher order Green’s function appearing in response function is decoupled using an appropriate decoupling scheme and the remaining Green’s function is evaluated via a renormalized electron and phonon Hamiltonian

(5)

Equation (3) can be written in simplified form as

with electron perturbed modeand renormalized electron mode frequencies

Now electron energy shift and line widths are obtainable as

where the superscript “D”, “3A”, “4A” and “ep” stand for the contributions due to defects, anharmonicities (cubic and quartic) and electron-phonon interactions, respectively.

The electron density of states (EDOS) in Lehman representation can be expressed as

Imaginary part of is given by

Using the imaginary part of Green’s function from Equation (23) in Equation (22) in Lehman representation we can write result 

Equation (24) can be reasonably approximated for small values of line width in the form

The various contributions to EDOS appears in Equation (26) can be summarized as

All the above expressions of density of states depend on the temperature. Let us examine EDOS in the following two regions.

Case-1 At the temperature very close to critical temperature, most of the electrons in superconducting state are paired. The collection of pairons (cooper pair, bipolarons) constitutes the condense or super fluid i.e. the pairons dominate over the normal electron

Case-2 In this case , the contribution of normal energy becomes least than the pairon energy, so that the pairon energy dominates over the unpaired energy of the system and transits to the superconducting state. The EDOS in this situation becomes as

The various symbols in above equations can be expressed as

is pairon distribution function.

Rest of the symbols are defined in references elsewhere [8].

The variation of EDOS for and is depicted in the Figures 1 and 2. These graphics are similar to the experimental imaging of quasi particle density of states for impurity induced EDOS and anharmonic EDOS. The central peak in Figure 1 exhibits the enhanced peak in the close vicinity of impurity site. This also reveals the evidence of four fold symmetric quasi particle cloud intensity peaks aligned with the nodes of the d-wave superconducting gap which is believed to characterize the HTSC and well supports the experimental observations of Pan et al. [11-12]. Figure 2 describes the effects of cubic anharmonicities on electron density of states. Four almost sharp peaks are found symmetrically distributed from both and axes. It also notable that the pairon distribution function, energy and temperature functions heavily influence the anharmonic contribution to EDOS as well as the other terms appearing in the EDOS. It emerges from the present investigations that EDOS not only depends on energy but also on renormalized frequency, anharmonicities, impurity concentration and

temperature and well support to d-wave superconductivity.

The authors (HS) and (AS) are thankful to Ministry of Human Resource development (MHRD), New Delhi and Council of Scientific and Industrial Research (CSIR), New Delhi for the financial support to carry out this research work.