Phonon Spectral Function Research Articles

Large impact of phonon lineshapes on the superconductivity of solid hydrogen

Phonon anharmonicity plays a crucial role in determining the stability and vibrational properties of high-pressure hydrides. Furthermore, strong anharmonicity can render phonon quasiparticle picture obsolete questioning standard approaches for modeling superconductivity in these material systems. In this work, we show the effects of non-Lorentzian phonon lineshapes on the superconductivity of high-pressure solid hydrogen. We calculate the superconducting critical temperature TC ab initio considering the full phonon spectral function and show that it overall enhances the TC estimate. The anharmonicity-induced phonon softening exhibited in spectral functions increases the estimate of the critical temperature, while the broadening of phonon lines due to phonon-phonon interaction decreases it. Our calculations also reveal that superconductivity emerges in hydrogen in the Cmca − 12 molecular phase VI at pressures between 450 and 500 GPa and explain the disagreement between the previous theoretical results and experiments.

Effects of anisotropy and disorder on the superconducting properties of niobium

We report results for the superconducting transition temperature and anisotropic energy gap for pure niobium based on Eliashberg’s equations and electron and phonon band structures computed from density functional theory. The electronic band structure is used to construct the Fermi surface and calculate the Fermi velocity at each point on the Fermi surface. The phonon bands are in excellent agreement with inelastic neutron scattering data. The corresponding phonon density of states and electron–phonon coupling define the electron–phonon spectral function, α2F(p, p′; ω), and the corresponding electron–phonon pairing interaction, which is the basis for computing the superconducting properties. The electron–phonon spectral function is in good agreement with existing tunneling spectroscopy data except for the spectral weight of the longitudinal phonon peak at ℏωLO = 23 meV. We obtain an electron–phonon coupling constant of λ = 1.057, renormalized Coulomb interaction μ⋆ = 0.218, and transition temperature Tc = 9.33 K. The corresponding strong-coupling gap at T = 0 is modestly enhanced, Δ0 = 1.55 meV, compared to the weak-coupling BCS value Δ0wc=1.78kBTc=1.43meV. The superconducting gap function exhibits substantial anisotropy on the Fermi surface. We analyze the distribution of gap anisotropy and compute the suppression of the superconducting transition temperature using a self-consistent T-matrix theory for quasiparticle-impurity scattering to describe niobium doped with non-magnetic impurities. We compare these results with experimental results on niobium SRF cavities doped with nitrogen impurities.

Local lattice distortions and the structural instabilities in bcc Nb–Ta–Ti–Hf high-entropy alloys: An ab initio computational study

Local lattice distortions (LLD) and structural stability of body-centered cubic (bcc) Nb–Ta–Ti–Hf high-entropy alloys (HEAs) are studied as functions of composition employing ab initio density-functional theory calculations, with specific focus on the role of the relative concentrations of group IV (Ti and Hf) versus group V (Nb and Ta) elements. Calculated results are presented as a function of composition x in NbxTa0.25Ti(0.75−x)/2Hf(0.75−x)/2 alloys, for elastic moduli, phonon spectral functions, LLD and structural energy differences for the bcc and competing hexagonal close-packed (hcp) and ω phases. The results highlight the important role of group V elements and LLD in stabilizing the bcc structure. They further reveal how composition x can be tuned to alter both the magnitude of the LLD and structural energy differences. Specifically, the magnitude of the structural energy differences, and elastic and dynamic stability of the bcc phase, are enhanced with increasing x, while the LLD increase in magnitude as this concentration is decreased. The results also show evidence of correlated LLD at lower values of x, reflecting local structural distortions towards the ω phase, but not hcp. The degree of ω-collapse is nevertheless partial i.e., transformation towards this phase is not observed to be complete due to the presence of Ta and Nb. At lower values of x we further find an energy landscape characterized by multiple, nearly degenerate local energy minima for different values of the LLD.

Dynamical Phonons Following Electron Relaxation Stages in Photoexcited Graphene.

Ultrafast electron-phonon relaxation dynamics in graphene hides many distinct phenomena, such as hot phonon generation, dynamical Kohn anomalies, and phonon decoupling, yet it still remains largely unexplored. Here, we unravel intricate mechanisms governing the vibrational relaxation and phonon dressing in graphene at a highly nonequilibrium state by means of first-principles techniques. We calculate dynamical phonon spectral functions and momentum-resolved line widths for various stages of electron relaxation and find photoinduced phonon hardening, overall increase of relaxation rate and nonadiabaticity, as well as phonon gain. Namely, the initial stage of photoexcitation is found to be governed by strong phonon anomalies of finite-momentum optical modes along with incoherent phonon production. The population inversion state, on the other hand, allows the production of coherent and strongly coupled phonon modes. Our research provides vital insights into the electron-phonon coupling phenomena in graphene and serves as a foundation for exploring nonequilibrium phonon dressing in materials where ordered states and phase transitions can be induced by photoexcitation.

Front Cover: Temperature‐Dependent Anharmonic Phonons in Quantum Paraelectric KTaO3 by First Principles and Machine‐Learned Force Fields (Adv. Quantum Technol. 4/2023)

The stability of the KTaO3 quantum paraelectric perovskite at low T (ordered structure in the centre) is studied beyond the harmonic approximation for the nuclear motion, with the inclusion of the zero-point quantum energy. Combining first-principles methods, quantum stochastic schemes, and machine learning force fields employing local radial and angular descriptors (green and purple arcs) it is possible to asses key boson quantities like the phonon spectral function (dispersion in the background) with great accuracy and reasonable time. More details can be found in article number 2200131 by Cesare Franchini and co-workers.

Dynamics and lifetimes of resonant phonons in the quasiparticle and nonquasiparticle regimes

The dynamics and lifetimes of a phonon with model resonant interactions from the quasiparticle to nonquasiparticle regime were investigated employing Green's function and the Green-Kubo method. In the weak-coupling case, the dynamics of the resonant phonon are analogous to a damped harmonic oscillator and its lifetime $\ensuremath{\tau}$ is in accordance with the standard phonon transport theory $\ensuremath{\tau}=\frac{1}{2\mathrm{\ensuremath{\Gamma}}}$ ($\mathrm{\ensuremath{\Gamma}}$ being the imaginary part of phonon self-energy). In the strong-coupling nonquasiparticle regime, however, the resonant phonon ``propagates'' in a complex form of wave packets and the actual phonon lifetime ${\ensuremath{\tau}}_{\mathrm{GK}}$ as determined from the Green-Kubo formula significantly deviates from the standard $\ensuremath{\tau}=\frac{1}{2\mathrm{\ensuremath{\Gamma}}}$ relation. Taking the four-phonon resonant ${\mathrm{AgCrSe}}_{2}$ and three-phonon resonant PbSe model systems as examples, the phonon nonquasiparticle dynamics and their lifetimes in real materials are further investigated by first-principle calculations. Substantial discrepancies between ${\ensuremath{\tau}}_{\mathrm{GK}}$ and $\ensuremath{\tau}$ are found for the strongly resonant phonons at high-symmetrical points in both materials. Meanwhile, the lifetime ${\ensuremath{\tau}}_{\mathrm{GK}}$ of the phonons that are not subjected to resonant phonon interactions almost recovers to the conventional theoretical result $\ensuremath{\tau}$ despite the non-Lorentzian spectral feature of these phonons. It is suggested that $\ensuremath{\omega}{\ensuremath{\tau}}_{\mathrm{GK}}$, instead of the conventional $\ensuremath{\omega}\ensuremath{\tau}$, can be applied as a criterion to distinguish phonon quasiparticles and nonquasiparticles in addition to the phonon spectral functions.

Lattice dynamics and elastic properties of α-U at high-temperature and high-pressure by machine learning potential simulations

Studying the physical properties of materials under high pressure and temperature through experiments is difficult. Theoretical simulations can compensate for this deficiency. Currently, large-scale simulations using machine learning force fields are gaining popularity. As an important nuclear energy material, the evolution of the physical properties of uranium under extreme conditions is still unclear. Herein, we trained an accurate machine learning force field on α-U and predicted the lattice dynamics and elastic properties at high pressures and temperatures. The force field agrees well with the ab initio molecular dynamics (AIMD) and experimental results and it exhibits higher accuracy than classical potentials. Based on the high-temperature lattice dynamics study, we first present the temperature-pressure range in which the Kohn anomalous behavior of the Σ4 optical mode exists. Phonon spectral function analysis showed that the phonon anharmonicity of α-U is very weak. We predict that the single-crystal elastic constants C44, C55, C66, polycrystalline modulus (E, G), and polycrystalline sound velocity (CL, CS) have strong heating-induced softening. All the elastic moduli exhibited compression-induced hardening behavior. The Poisson's ratio shows that it is difficult to compress α-U at high pressures and temperatures. Moreover, we observed that the material becomes substantially more anisotropic at high pressures and temperatures. The accurate predictions of α-U demonstrate the reliability of the method. This versatile method facilitates the study of other complex metallic materials.

Off-diagonal electron-phonon coupling effects on two-level system dynamics

Electromagnetically coupled two-level systems play a central role in several condensed matter components being considered for quantum information processing applications. If the states couple to phonon excitations, their electromagnetic response is altered via phonon-assisted transitions and lifetime broadening. The former has been treated extensively for a number of specific two-level systems (e.g., excitons in artificial quantum dots, localized states associated with impurities or defects, etc.), but the latter has received less attention. Here we study a microscopic model of the dipole transition of a two-level system under the influence of both diagonal and nondiagonal interactions with a bath of phonons. Our results capture both the influence of the frequency distribution of phonons on the relative spectral weight of the zero-phonon transition and the phonon sidebands, and the broadening of the zero-phonon line due to nondiagonal electron-phonon coupling. We use a formalism that includes non-Markovian effects related to the feedback mechanism between the two-level system and the phonon bath. For simplified forms of the phonon spectral functions we provide analytical expressions up to second order in the coupling strength that demonstrate the importance of including both forms of electron-phonon coupling in studies of these systems. Our formalism can be generalized to higher orders of coupling and for realistic phonon spectral functions.

Detection of squeezed phonons in pump-probe spectroscopy

Robust engineering of phonon squeezed states in optically excited solids has emerged as a promising tool to control and manipulate their properties. However, in contrast to quantum optical systems, detection of phonon squeezing is subtle and elusive, and an important question is what constitutes an unambiguous signature of it. The state of the art involves observing oscillations at twice the phonon frequency in time resolved measurements of the out of equilibrium phonon fluctuation. Using Keldysh formalism we show that such a signal is a necessary but not a sufficient signature of a squeezed phonon, since we identify several mechanisms that do not involve squeezing and yet which produce similar oscillations. We show that a reliable detection requires a time and frequency resolved measurement of the phonon spectral function.

The Pseudopotential Approach within Density-Functional Theory: The Case of Atomic Metallic Hydrogen

Internal energies, enthalpies, phonon dispersion curves, and superconductivity of atomic metallic hydrogen are calculated. The standard use of pseudopotentials in density-functional theory are compared with full Coulomb-potential, all-electron linear muffin-tin orbital calculations. Quantitatively similar results are found as far as internal energies are concerned. Larger differences are found for phase-transition pressures; significant enough to affect the phase diagram. Electron–phonon spectral functions α2F(ω) also show significant differences. Against expectation, the estimated superconducting critical-temperature Tc of the first atomic metallic phase I41/amd (Cs-IV) at 500 GPa is actually higher.

Finite-temperature density-matrix renormalization group method for electron-phonon systems: Thermodynamics and Holstein-polaron spectral functions

We investigate the thermodynamics and finite-temperature spectral functions of the Holstein polaron using a density-matrix renormalization group method. Our method combines purification and local basis optimization (LBO) as an efficient treatment of phonon modes. LBO is a scheme which relies on finding the optimal local basis by diagonalizing the local reduced density matrix. By transforming the state into this basis, one can truncate the local Hilbert space with a negligible loss of accuracy for a wide range of parameters. In this work, we focus on the crossover regime between large and small polarons of the Holstein model. Here, no analytical solution exists and we show that the thermal expectation values at low temperatures are independent of the phonon Hilbert space truncation provided the basis is chosen large enough. We then demonstrate that we can extract the electron spectral function and establish consistency with results from a finite-temperature Lanczos method. We additionally calculate the electron emission spectrum and the phonon spectral function and show that all the computations are significantly simplified by the local basis optimization. We observe that the electron emission spectrum shifts spectral weight to both lower frequencies and larger momenta as the temperature is increased. The phonon spectral function experiences a large broadening and the polaron peak at large momenta gets significantly flattened and merges almost completely into the free-phonon peak.

Anharmonicity and Ultralow Thermal Conductivity in Lead-Free Halide Double Perovskites.

The lead-free halide double perovskite class of materials offers a promising venue for resolving issues related to toxicity of Pb and long-term stability of the lead-containing halide perovskites. We present a first-principles study of the lattice vibrations in Cs_{2}AgBiBr_{6}, the prototypical compound in this class and show that the lattice dynamics of Cs_{2}AgBiBr_{6} is highly anharmonic, largely in regards to tilting of AgBr_{6} and BiBr_{6} octahedra. Using an energy- and temperature-dependent phonon spectral function, we then show how the experimentally observed cubic-to-tetragonal phase transformation is caused by the collapse of a soft phonon branch. We finally reveal that the softness and anharmonicity of Cs_{2}AgBiBr_{6} yield an ultralow thermal conductivity, unexpected of high-symmetry cubic structures.

Polaron transformations in the realistic model of the strongly correlated electron system

Electron-phonon coupling, diagonal in a real space formulation, leads to polaron paradigm of smoothly varying properties. However, fundamental changes, namely the singular behavior of polarons, occur if non-diagonal pairing is involved into consideration. The study of polaron transformations and related properties of matter is of particular interest for realistic models, since competition between diagonal and non-diagonal electron-phonon contributions in the presence of other strong interactions can result in unconventional behavior of the system. Here we consider the multiband pd-model of cuprate superconductors with electron-phonon interaction and analyze the features of the systems that are caused by the competition of diagonal and non-diagonal electron-phonon contributions in the limit of strong electron correlations. Using the polaronic version of the generalized tight-binding method, we describe the evolution of the band structure, Fermi surface, density of states at Fermi level, and phonon spectral function in the space of electron-phonon parameters ranging from weak to strong coupling strength of the adiabatic limit. On the phase diagram of polaron properties we reveal two quantum phase transitions and show how electron-phonon interaction gives rise to Fermi surface transformation (i) from hole pockets to Fermi arcs and (ii) from hole to electron type of conductivity. We also demonstrate the emergence of new states in the phonon spectral function of the polaron and discuss their origin.

Quantifying phonon-induced non-Markovianity in color centers in diamond

The degree of non-Markovianity of a continuous bath can be quantified by means of the coherence. This simple measure is experimentally accessible through Ramsey spectroscopy, but it is limited to incoherent dynamical maps. We propose an extension of this measure and discuss its application to color centers in diamond, where the optical coherence between two orbital states is affected by interactions with a structured phonon bath. By taking realistic phonon spectral density functions into account, we show that this measure is well-behaved at arbitrary temperatures and that it provides additional insights about how non-Markoviantiy is affected by the presence of both bulk and quasi-localized phonon modes. Importantly, with only a little overhead the measure can be adapted to eliminate the false signs of non-Markovianity from coherent dynamical maps and is thus applicable for a large class of systems modeled by the spin-boson Hamiltonian.

Electronic and phononic properties of a two-dimensional electron gas coupled to dipolar phonons via small-momentum-transfer scattering

We study the electron and phonon spectral functions as well as superconductivity for a two dimensional electron gas couped to dipolar phonons via small-momentum transfer scattering. The results reported here are obtained through sign-problem-free quantum Monte-Carlo simulation. Hence aside from modeling there is no further approximations. Our results can be relevant to the high-temperature superconductor at the interface between a unit-cell-thick FeSe and SrTiO3.

High-Tc superconductivity in H3S: pressure effects on the superconducting critical temperature and Cooper pair distribution function

We use first-principles calculations to study the effects of pressure on the vibrational and superconducting properties of H3S in the cubic phase for the pressure range where the superconducting critical temperature (Tc) was measured (155–225 GPa). The pressure effects were incorporated using the functional derivative method (FDM). In this paper we present for the first time the Cooper pair distribution functions Dcp(ω, Tc) for H3S, which will allow us to identify the spectral regions where Cooper pair formation at temperature Tc is more favorable. We analyze in detail the effects of pressure on the electron–phonon spectral density function α2F(ω) and the phonon density of states and their relationship with Tc. The FDM manages to reproduce the trend of the pressure dependence of critical temperature, in good agreement with experimental data in the range of 155 to 190 GPa. Dcp(ω, Tc) suggests that the low-frequency vibration region is where Cooper pairs are possible, which means that S-vibrations have an important role in the superconductivity properties of H3S.

High-pressure phase diagram of hydrogen and deuterium sulfides from first principles: Structural and vibrational properties including quantum and anharmonic effects

We study the structural and vibrational properties of the high-temperature superconducting sulfur trihydride and trideuteride in the high-pressure $Im\bar{3}m$ and $R3m$ phases by first-principles density-functional-theory calculations. On lowering pressure, the rhombohedral transition $Im\bar{3}m \rightarrow R3m$ is expected, with hydrogen bond desymmetrization and occurrence of trigonal lattice distortion. In hydrostatic conditions we find that, contrary to what suggested in some recent experiments, if the rhombohedral distortion exists it affects mainly the hydrogen-bonds, whereas the resulting cell distortion is minimal. We estimate that the occurrence of a stress anisotropy of approximately $10\%$ could explain this discrepancy. Assuming hydrostatic conditions, we calculate the critical pressure at which the rhombohedral transition occurs. Quantum and anharmonic effects, which are relevant in this system, are included at nonperturbative level with the stochastic self-consistent harmonic approximation (SSCHA). Within this approach, we determine the transition pressure by calculating the free energy Hessian. We find that quantum anharmonic effects are responsible for a strong reduction of the critical pressure with respect to the one obtained with the classical harmonic approach. Moreover, we observe a prominent isotope effect, as we estimate higher pressure transition for D${}_3$S than for H${}_3$S. Finally, within SSCHA we calculate the anharmonic phonon spectral functions in the $Im\bar{3}m$ phase. The strong anharmonicity of the system is confirmed by the occurrence of very large anharmonic broadenings leading to complex non-Lorentzian line shapes. However, for the vibrational spectra at zone center, accessible e.g. by infrared spectroscopy, the broadenings are very small (linewidth at most around 2~meV) and anharmonic phonon quasiparticles are well defined.

Influence of phonon-phonon coupling on superconducting state in honeycomb-type crystal lattice

We have taken into account the superconducting state inducing in the crystal lattice of the honeycomb-type. In the framework of the Eliashberg theory, we have determined the thermodynamic properties of the system. The phonon spectral function, which is the input parameter in the Eliashberg equations, has been calculated by using the thermodynamic Green functions. We have considered the model of the coupled Einstein oscillators with frequency ω0=100 meV. We have shown that the increasing inter-phonon coupling constant (f) causes the rapid growth of the critical temperature ([TC]max=36.2 K) just below the maximum value of f equal to 0.25ω0. Simultaneously, the order parameter and the thermodynamic critical field take the values increasingly distant from the predictions of the BCS theory, which results from the strong-coupling and the retardation effects.

Phonon broadening in high entropy alloys

Refractory high entropy alloys feature outstanding properties making them a promising materials class for next-generation high-temperature applications. At high temperatures, materials properties are strongly affected by lattice vibrations (phonons). Phonons critically influence thermal stability, thermodynamic and elastic properties, as well as thermal conductivity. In contrast to perfect crystals and ordered alloys, the inherently present mass and force constant fluctuations in multi-component random alloys (high entropy alloys) can induce significant phonon scattering and broadening. Despite their importance, phonon scattering and broadening have so far only scarcely been investigated for high entropy alloys. We tackle this challenge from a theoretical perspective and employ ab initio calculations to systematically study the impact of force constant and mass fluctuations on the phonon spectral functions of 12 body-centered cubic random alloys, from binaries up to 5-component high entropy alloys, addressing the key question of how chemical complexity impacts phonons. We find that it is crucial to include both mass and force constant fluctuations. If one or the other is neglected, qualitatively wrong results can be obtained such as artificial phonon band gaps. We analyze how the results obtained for the phonons translate into thermodynamically integrated quantities, specifically the vibrational entropy. Changes in the vibrational entropy with increasing the number of elements can be as large as changes in the configurational entropy and are thus important for phase stability considerations. The set of studied alloys includes MoTa, MoTaNb, MoTaNbW, MoTaNbWV, VW, VWNb, VWTa, VWNbTa, VTaNbTi, VWNbTaTi, HfZrNb, HfMoTaTiZr.

Critical temperature of metallic hydrogen sulfide at 225-GPa pressure

The Eliashberg theory generalized for electron—phonon systems with a nonconstant density of electron states and with allowance made for the frequency behavior of the electron mass and chemical potential renormalizations is used to study T c in the SH3 phase of hydrogen sulfide under pressure. The phonon contribution to the anomalous electron Green’s function is considered. The pairing within the total width of the electron band and not only in a narrow layer near the Fermi surface is taken into account. The frequency and temperature dependences of the complex mass renormalization ReZ(ω), the density of states N(e) renormalized by the electron—phonon interactions, and the electron—phonon spectral function obtained computationally are used to calculate the anomalous electron Green’s function. A generalized Eliashberg equation with a variable density of electron states has been solved. The frequency dependence of the real and imaginary parts of the order parameter in the SH3 phase has been obtained. The value of T c ≈ 177 K in the SH3 phase of hydrogen sulfide at pressure P = 225 GPa has been determined by solving the system of Eliashberg equations.