Beyond Eliashberg superconductivity in MgB2: anharmonicity, two-phonon scattering, and multiple gaps

Density-functional calculations of the phonon spectrum and electron-phonon coupling in MgB (2) are presented. The E(2g) phonons, which involve in-plane B displacements, couple strongly to the p(x,y) electronic bands. The isotropic electron-phonon coupling constant is calculated to be about 0.8. Allowing for different order parameters in different bands, the superconducting lambda in the clean limit is calculated to be significantly larger. The E(2g) phonons are strongly anharmonic, and the nonlinear contribution to the coupling between the E(2g) modes and the p(x,y) bands is significant.     

Phonon dispersion and lifetimes in MgB2

We measure phonon dispersion and linewidth in a single crystal of MgB2 along the Gamma-A, Gamma-M, and A-L directions using inelastic x-ray scattering. We use density functional theory to compute the effect of both electron-phonon coupling and anharmonicity on the linewidth, obtaining excellent agreement with experiment. Anomalous broadening of the E(2g) phonon mode is found all along Gamma-A. The dominant contribution to the linewidth is always the electron-phonon coupling.     

Observation of interband pairing interaction in a two-band superconductor: MgB2

The recently discovered anisotropic superconductor MgB2 is the first of its kind showing the intriguing properties of two-band superconductivity. By tunneling experiments using thin film tunnel junctions, electron-coupled phonon spectra were determined showing that superconductivity in MgB2 is phonon mediated. In a further analysis, which involves first principles calculations, the strongest feature in these spectra could be traced back to the key quantity of two-band superconductivity, the interband pairing interaction. For the phonons, this interaction turns out quite selective. It involves mainly low-energy optical phonon modes, where the boron atoms move perpendicular to the boron planes.     

Nonlocal screening, electron-phonon coupling, and phonon renormalization in metals

A method for calculating the phonon self-energy in metals arising from the coupling between phonons and electrons near the Fermi surface is developed. The essence of this scheme is the separation of the inter- and intraband parts of the electron polarizability. The intraband contribution provides extra screenings and is closely related to the electron-phonon coupling and phonon softening in metals. Applications of this scheme to phonons in MgB2 give excellent results when compared with experiments and previous theoretical work. In addition, both electron and hole dopings are found to reduce the renormalization effect of the E(2g) phonon mode, which indicates a weakened electron-phonon coupling in the doped systems.     

Raman scattering from a superconductivity-induced bound state in MgB2

It is shown that the sharp peak in the E(2g) Raman spectrum of superconducting MgB2 is due to a bound state caused by the electron-phonon coupling. Our theory explains why this peak appears only in the spectra with E(2g) symmetry and only in the sigma but not the pi bands. The properties of the bound state and the Raman spectrum are investigated, also in the presence of impurity scattering.     

Enhancement of the superconducting transition temperature of MgB2 by a strain-induced bond-stretching mode softening

We report a systematic increase of the superconducting transition temperature T(c) with a biaxial tensile strain in MgB2 films to well beyond the bulk value. The tensile strain increases with the MgB2 film thickness, caused primarily by the coalescence of initially nucleated discrete islands (the Volmer-Weber growth mode.) The T(c) increase was observed in epitaxial films on SiC and sapphire substrates, although the T(c) values were different for the two substrates due to different lattice parameters and thermal expansion coefficients. We identified, by first-principles calculations, the underlying mechanism for the T(c) increase to be the softening of the bond-stretching E(2g) phonon mode, and we confirmed this conclusion by Raman scattering measurements. The result suggests that the E(2g) phonon softening is a possible avenue to achieve even higher T(c) in MgB2-related material systems.     

Superconductivity of metallic boron in MgB2

Boron in MgB2 forms stacks of honeycomb layers with magnesium as a space filler. Band structure calculations indicate that Mg is substantially ionized, and the bands at the Fermi level derive mainly from B orbitals. Strong bonding with an ionic component and considerable metallic density of states yield a sizable electron-phonon coupling. Together with high phonon frequencies, which we estimate via zone-center frozen phonon calculations to be between 300 and 700 cm(-1), this produces a high critical temperature, consistent with recent experiments. Thus MgB2 can be viewed as an analog of the long sought, but still hypothetical, superconducting metallic hydrogen.     

Lattice dynamics of FeSb2

The lattice dynamics of FeSb(2) is investigated by first-principles density functional theory calculations and Raman spectroscopy. All Raman- and infrared-active phonon modes are properly assigned. The calculated and measured phonon energies are in good agreement. We have observed strong mixing of the A(g) symmetry modes, with the intensity exchange in the temperature range 210 and 260 K. The A(g) mode repulsion increases by doping FeSb(2) with Co, with no signatures of the electron-phonon interaction for these modes.     

Phonon anharmonicities in graphite and graphene

We determine from first principles the finite-temperature properties-linewidths, line shifts, and lifetimes-of the key vibrational modes that dominate inelastic losses in graphitic materials. In graphite, the phonon linewidth of the Raman-active E(2g) mode is found to decrease with temperature; such anomalous behavior is driven entirely by electron-phonon interactions, and does not appear in the nearly degenerate infrared-active E(1u) mode. In graphene, the phonon anharmonic lifetimes and decay channels of the A(1)' mode at K dominate over E(2g) at Gamma and couple strongly with acoustic phonons, highlighting how ballistic transport in carbon-based interconnects requires careful engineering of phonon decays and thermalization.     

Phonon dispersion and electron-phonon coupling in MgB2 and AlB2.

We present a first principles investigation of the lattice dynamics and electron-phonon coupling of the superconductor MgB2 and the isostructural AlB2 within the framework of density functional perturbation theory using a mixed-basis pseudopotential method. Complete phonon dispersion curves and Eliashberg functions alpha2F are calculated for both systems. The main differences are related to high frequency in-plane boron vibrations, which are strongly softened in MgB2 and exhibit an exceptionally strong electron-phonon coupling. We also report on Raman measurements, which support the theoretical findings. Implications for the superconducting transition temperature are briefly discussed.     

Three-dimensional MgB2-type superconductivity in hole-doped diamond

We substantiate by numerical and analytical calculations that the recently discovered superconductivity below 4 K in 3% boron-doped diamond is caused by electron-phonon coupling of the same type as in MgB2, albeit in three dimensions. Holes at the top of the zone-centered, degenerate sigma-bonding valence-band couple strongly to the optical bond-stretching modes. The increase from two to three dimensions reduces the mode softening crucial for T(c) reaching 40 K in MgB2. Even if diamond had the same bare coupling constant as MgB2, which could be achieved with 10% doping, T(c) would be only 25 K. Superconductivity above 1 K in Si (Ge) requires hole doping beyond 5% (10%).     

Superconductivity in boron-doped diamond

Superconductivity of boron-doped diamond, reported recently at T(c)=4 K, is investigated exploiting its electronic and vibrational analogies to MgB2. The deformation potential of the hole states arising from the C-C bond-stretch mode is 60% larger than the corresponding quantity in MgB2 that drives its high T(c), leading to very large electron-phonon matrix elements. The calculated coupling strength lambda approximately 0.5 leads to T(c) in the 5-10 K range and makes phonon coupling the likely mechanism. Higher doping should increase T(c) somewhat, but the effects of three dimensionality primarily on the density of states keep doped diamond from having a T(c) closer to that of MgB2.     

Temperature and layer number dependence of the G and 2D phonon energy and damping in graphene

We have studied the temperature and size dependence of the G and 2D phonon modes in graphene. It is shown that in a graphene monolayer the phonon energy decreases whereas the phonon damping increases with increasing temperature. The electron-phonon interaction leads to hardening whereas the fourth-order anharmonic phonon-phonon processes lead to softening of the phonon energy with increasing temperature. We have shown that the electron-phonon interaction plays an important role also by the dispersion dependence of the phonon G mode, by the observation of the Kohn anomaly. The G mode frequency decreases and damping increases, whereas the 2D phonon frequency and damping increase with increasing layer number. The temperature and size effects of the 2D mode are much stronger than those of the G mode.     

Specific heat and thermal conductivity in the mixed state of MgB2

The specific heat C and the electronic and phononic thermal conductivities kappa(e) and kappa(ph) are calculated in the mixed state for magnetic fields H near H(c2), including the effects of supercurrent flow and Andreev scattering. The resulting function C(H) is nearly linear while kappa(e)(H) exhibits an upward curvature near H(c2). The slopes decrease with impurity scattering which improves the agreement with the data on MgB2. The ratio of phonon relaxation times tau(n)/tau(s)=g(omega(0),H) for phonon energy omega(0) is smeared out around omega(0)=2Delta and tends to one for increasing H. This leads to a rapid reduction of kappa(ph)(H) in MgB2 for relatively small fields due to the rapid suppression of the smaller energy gap.     

Phonon States and Dispersive Spectra of Polar Optical Phonons in Quasi-One-Dimensional Nanowires of Wurtzite ZnO and Zinc-Blend MgO Semiconductors

Within the framework of the macroscopic dielectric continuum model and Loudon's uniaxial crystal model, the phonon modes of a wurtzite/zinc-blende one-dimensional (1D) cylindrical nanowire (NW) are derived and studied. The analytical phonon states of phonon modes are given. It is found that there exist two types of polar phonon modes, i.e. interface optical (IO) phonon modes and the quasi-confined (QC) phonon modes existing in 1D wurtzite/zinc-blende NWs. Via the standard procedure of field quantization, the Fröhlich electron-phonon interaction Hamiltonians are obtained. Numerical calculations of dispersive behavior of these phonon modes on a wurtzite/zinc-blende ZnO/MgO NW are performed. The frequency ranges of the IO and QC phonon modes of the ZnO/MgO NWs are analyzed and discussed. It is found that the IO modes only exist in one frequency range, while QC modes may appear in three frequency ranges. The dispersive properties of the IO and QC modes on the free wave-number k_z and the azimuthal quantum number m arediscussed. The analytical Hamiltonians of electron-phonon interaction obtained here are quite useful for further investigating phonon influence on optoelectronics properties of wurtzite/zinc-blende 1D NW structures.     

Phonon density of states and phonon dispersion of superconducting MgB2

For the superconductor MgB2, we have calculated the phonon density of states (DOS), phonon dispersion and Eliashberg function throughout the Brillouin zone (BZ), using an empirical potential model. The calculated values are consistent with the theoretical and experimental values. The calculated results show our empirical potential model is available for MgB2.     

Structure and superconductivity in the ternary silicide CaAlSi

Using the linear response-linearized Muffin-tin orbital (LR-LMTO) method, we study the electronic band structure, phonon spectra, electron-phonon coupling and superconductivity for c-axis ferromagnetic-like (F-like) and antiferromagnetic-like (AF-like) structures in ternary silicide CaAlSi. The following conclusions are drawn from our calculations. If Al and Si atoms are assumed to arrange along the c axis in an F-like long-range ordering (-Al-Al-Al-and-Si-Si-Si-), one could obtain the ultrasoft B_1g phonon mode and thus very strong electron-phonon coupling in CaAlSi. However, the appearance of imaginary frequency phonon modes indicates the instability of such a structure. For Al and Si atoms arranging along the c axis in an AF-like long-range ordering (-Al-Si-Al-), the calculated electron-phonon coupling constant is equal to 0.8 and the logarithmically averaged frequency is 146.8 K. This calculated result can correctly yield the superconducting transition temperature of CaAlSi by the standard BCS theory in the moderate electron-phonon coupling strength. We propose that an AF-like superlattice model for Al (or Si) atoms along the c direction may mediate the inconsistency estimated from theory and experiment, and explain the anomalous superconductivity in CaAlSi.      

Temperature dependence of the band gap of semiconducting carbon nanotubes

The temperature dependence of the band gap of semiconducting single-wall carbon nanotubes (SWNTs) is calculated by direct evaluation of electron-phonon couplings within a "frozen-phonon" scheme. An interesting diameter and chirality dependence of E(g)(T) is obtained, including nonmonotonic behavior for certain tubes and distinct "family" behavior. These results are traced to a strong and complex coupling between band-edge states and the lowest-energy optical phonon modes in SWNTs. The E(g)(T) curves are modeled by an analytic function with diameter- and chirality-dependent parameters; these provide a valuable guide for systematic estimates of E(g)(T) for any given SWNT. The magnitudes of the temperature shifts at 300 K are smaller than 12 meV and should not affect (n,m) assignments based on optical measurements.     

铁基氟化物超导体SrFe_(1-x)Co_xAsF(x=0,0.125)声子特性的第一性原理计算研究

采用基于第一性原理的平面波赝势方法,计算了铁基氟化物及其钴掺杂超导体SrFe1-xCoxAsF(x=0,0.125)在四方非磁态与正交条纹反铁磁态下的声子谱(声子色散曲线、声子态密度)及电-声子耦合常数.计算发现:条纹反铁磁相互作用下的自旋-声子耦合效应强于电-声子耦合作用使声子谱的宽度减小;自旋效应使声子的有效质量增加导致条纹反铁磁态下Fe原子与As原子的耦合振动频率减小.另外,掺杂和自旋效应是提高电-声子耦合常数的两个有效方法,但计算所得超导转变温度远小于实验测量值,表明铁基超导电性非简单的电-声子耦合配对机理.     

Kohn anomalies and electron-phonon interactions in graphite

We demonstrate that graphite phonon dispersions have two Kohn anomalies at the Gamma-E(2g) and K-A'1 modes. The anomalies are revealed by two sharp kinks. By an exact analytic derivation, we show that the slope of these kinks is proportional to the square of the electron-phonon coupling (EPC). Thus, we can directly measure the EPC from the experimental dispersions. The Gamma-E(2g) and K-A'1 EPCs are particularly large, while they are negligible for all the other modes at Gamma and K.