Author Topic: phonon spectral function  (Read 2449 times)

0 Members and 1 Guest are viewing this topic.

Offline Roc2019

  • Heavy QuantumATK user
  • ***
  • Posts: 53
  • Country: cn
  • Reputation: 0
    • View Profile
phonon spectral function
« on: May 1, 2020, 01:56 »
Dear Sir,
Is it possible to obtain the phonon spectral function (includinng the weight of each phonon mode) like the attachment?

Thank you so much.

Offline Petr Khomyakov

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 1290
  • Country: dk
  • Reputation: 25
    • View Profile
Re: phonon spectral function
« Reply #1 on: May 5, 2020, 11:36 »
Please clarify how you define the weight of phonon bands and what is your system(s) of interest.

Offline Roc2019

  • Heavy QuantumATK user
  • ***
  • Posts: 53
  • Country: cn
  • Reputation: 0
    • View Profile
Re: phonon spectral function
« Reply #2 on: May 7, 2020, 05:50 »
Dear Dr.  Khomyakov,

I mean the spectral function of phonon like  the Fig. 6 in the following paper

PHYSICAL REVIEW MATERIALS 4, 045403 (2020)

So, how can one obtain that picture by using the newest version 2019.12.sp1. ?  (for example, http://ollehellman.github.io/page/workflows/minimal_example_1.html), is it beyond the newest version?

Another, how can one  obtain the electron and phonon linewidth? (for example, the Fig.3 to Fig. 6 obtained by EPW  http://epw.org.uk/Documentation/B-dopedDiamond)

Thank you so much.
« Last Edit: May 7, 2020, 11:44 by Roc2019 »

Offline Troels-Markussen

  • QuantumATK Staff
  • Heavy QuantumATK user
  • *****
  • Posts: 41
  • Country: dk
  • Reputation: 3
  • QuantumATK staff
    • View Profile
Re: phonon spectral function
« Reply #3 on: May 12, 2020, 09:23 »
Hi, It is currently not possible in QuantumATK to calculate the phonon linewidth due to anharmonic phonon-phonon interactions as  done here http://ollehellman.github.io/page/workflows/minimal_example_1.html . In principle you can calculate the electron linewidth due to electron-phonon coupling.  It will require some scripting by yourself, but we might make an actual implementation next year. The way to do it would be something like:
Code
# Calculate electron-phonon coupling.
# This requires a HamiltonianDerivatives and a DynamicalMatrix object.
# Also, specify the k-points at which to calculate the broadening.
# Manually specify the bands to consider (speed up computation time).

electron_phonon_coupling = ElectronPhononCoupling(
    hamiltonian_derivatives,
    dynamical_matrix,
    kpoints_fractional=kpoint,
    qpoints_fractional=MonkhorstPackGrid(30,30,30),
    electron_bands=[3, 4, 5],
    energy_tolerance=1.0*eV,
    initial_state_energy_range=[-100, 100]*eV,
    )

# Calculate the inverse life-time (broadening)
inverse_life_time = electron_phonon_coupling.inverseLifeTime()
The calculated inverse life time corresponds to the linewidth calculated with the EPW code. Regards, Troels

Offline Roc2019

  • Heavy QuantumATK user
  • ***
  • Posts: 53
  • Country: cn
  • Reputation: 0
    • View Profile
Re: phonon spectral function
« Reply #4 on: May 12, 2020, 14:24 »
Thank you Dr. Markussen .