AlAs Phonon bandsΒΆ

This example shows how to plot the phonon band structure of AlAs. See tutorial/lesson_rf2.html



Warning: file /Users/gmatteo/git_repos/abipy/abipy/data/refs/alas_phonons/ does not contain atomic_numbers.
Particular methods need them!
/Users/gmatteo/git_repos/pymatgen/pymatgen/util/ UserWarning: Matplotlib is currently using agg, which is a non-GUI backend, so cannot show the figure.

from abipy.abilab import abiopen
import as abidata

# Open PHBST file produced by anaddb and extract the phonon bands object.
# (alternatively one can use the shell and ` -nb`
# to open the file in a jupyter notebook.
with abiopen(abidata.ref_file("")) as ncfile:
    phbands = ncfile.phbands

# Plot the phonon frequencies. Note that the labels for the q-points
# are found automatically in an internal database.
phbands.plot(units="cm-1", title="AlAs without LO-TO splitting")

# Alternatively you can use the optional argument qlabels
# that defines the mapping reduced_coordinates --> name of the q-point.
#qlabels = {
#    (0,0,0): "$\Gamma$",
#    (0.375, 0.375, 0.7500): "K",
#    (0.5, 0.5, 1.0): "X",
#    (0.5, 0.5, 0.5): "L",
#    (0.5, 0.0, 0.5): "X",
#    (0.5, 0.25, 0.75): "W",

# and pass it to the plot method:
#phbands.plot(title="AlAs without LO-TO splitting", qlabels=qlabels)

Total running time of the script: ( 0 minutes 0.877 seconds)

Gallery generated by Sphinx-Gallery