Dos

class py4vasp.data.Dos(raw_data)

Bases: py4vasp.data._base.DataBase, py4vasp.data._export.Image

The electronic density of states (DOS).

You can use this class to extract the DOS data of a Vasp calculation. Typically you want to run a non self consistent calculation with a denser mesh for a smoother DOS, but the class will work independent of it. If you generated orbital decomposed DOS, you can use this class to select which subset of these orbitals to read or plot.

Parameters

raw_dos (RawDos) – Dataclass containing the raw data necessary to produce a DOS.

Methods Summary

from_dict(dict_)

Initialize refinement class from data dictionary

from_file([file])

Read the data dictionary from the given file.

plot([selection])

Read the data and generate a plotly figure.

print()

Print a string representation of this class to standard output.

read([selection])

Read the data into a dictionary.

to_dict([selection])

Read the data into a dictionary.

to_frame([selection])

Read the data into a pandas DataFrame.

to_image(*args[, filename])

Read the data and generate an image writing to the given filename.

to_plotly([selection])

Read the data and generate a plotly figure.

Methods Documentation

classmethod from_dict(dict_)

Initialize refinement class from data dictionary

Parameters

data_dict (dict) – Data dictionary that contains one or more different raw data sources.

classmethod from_file(file=None)

Read the data dictionary from the given file.

You want to use this method if you want to avoid using the Calculation wrapper, for example because you renamed the output of the VASP calculation.

Parameters

file (str or Path or raw.File) – Filename from which the data is extracted. If not present the default filename is used. Alternatively, you can open the file yourself and pass the File object. In that case, you need to take care the file is properly closed again and be aware the generated instance of this class becomes unusable after the file is closed.

Returns

The returned instance handles opening and closing the file for every function called on it, unless a File object in which case this is left to the user.

Return type

DataBase

plot(selection=None)

Read the data and generate a plotly figure.

Parameters
  • selection (str) –

    A string specifying the projection of the orbitals. There are three distinct possibilities:

    • To specify the atom, you can either use its element name (Si, Al, …) or its index as given in the input file (1, 2, …). For the latter option it is also possible to specify ranges (e.g. 1:4).

    • To select a particular orbital you can give a string (s, px, dxz, …) or select multiple orbitals by their angular momentum (s, p, d, f).

    • For the spin, you have the options up, down, or total.

    You separate multiple selections by commas or whitespace and can nest them using parenthesis, e.g. Sr(s, p) or s(up), p(down). The order of the selections does not matter, but it is case sensitive to distinguish p (angular momentum l = 1) from P (phosphorus).

  • source (str, optional) – If you used a KPOINTS_OPT file to use a second k-point mesh, you can provide a keyword argument source=”kpoints_opt” to use the k-points defined in that file instead of the one specified in KPOINTS.

Returns

plotly figure containing the total DOS. If the calculation was spin polarized, the resulting DOS is spin resolved and the spin-down DOS is plotted towards negative values. If a selection the orbital resolved DOS is given for the specified projectors.

Return type

plotly.graph_objects.Figure

Examples

Select the p orbitals of the first atom in the POSCAR file:

>>> calc.dos.plot(selection="1(p)")

Select the d orbitals of Mn, Co, and Fe:

>>> calc.dos.plot("d(Mn, Co, Fe)")

Select the spin-up contribution of the first three atoms combined

>>> calc.dos.plot("up(1:3)")
print()

Print a string representation of this class to standard output.

read(selection=None)

Read the data into a dictionary.

Parameters
  • selection (str) –

    A string specifying the projection of the orbitals. There are three distinct possibilities:

    • To specify the atom, you can either use its element name (Si, Al, …) or its index as given in the input file (1, 2, …). For the latter option it is also possible to specify ranges (e.g. 1:4).

    • To select a particular orbital you can give a string (s, px, dxz, …) or select multiple orbitals by their angular momentum (s, p, d, f).

    • For the spin, you have the options up, down, or total.

    You separate multiple selections by commas or whitespace and can nest them using parenthesis, e.g. Sr(s, p) or s(up), p(down). The order of the selections does not matter, but it is case sensitive to distinguish p (angular momentum l = 1) from P (phosphorus).

  • source (str, optional) – If you used a KPOINTS_OPT file to use a second k-point mesh, you can provide a keyword argument source=”kpoints_opt” to use the k-points defined in that file instead of the one specified in KPOINTS.

Returns

Contains the energies at which the DOS was evaluated aligned to the Fermi energy and the total DOS or the spin-resolved DOS for spin-polarized calculations. If available and a selection is passed, the orbital resolved DOS for the selected orbitals is included.

Return type

dict

Examples

Select the p orbitals of the first atom in the POSCAR file:

>>> calc.dos.read(selection="1(p)")

Select the d orbitals of Mn, Co, and Fe:

>>> calc.dos.read("d(Mn, Co, Fe)")

Select the spin-up contribution of the first three atoms combined

>>> calc.dos.read("up(1:3)")
to_dict(selection=None)

Read the data into a dictionary.

Parameters
  • selection (str) –

    A string specifying the projection of the orbitals. There are three distinct possibilities:

    • To specify the atom, you can either use its element name (Si, Al, …) or its index as given in the input file (1, 2, …). For the latter option it is also possible to specify ranges (e.g. 1:4).

    • To select a particular orbital you can give a string (s, px, dxz, …) or select multiple orbitals by their angular momentum (s, p, d, f).

    • For the spin, you have the options up, down, or total.

    You separate multiple selections by commas or whitespace and can nest them using parenthesis, e.g. Sr(s, p) or s(up), p(down). The order of the selections does not matter, but it is case sensitive to distinguish p (angular momentum l = 1) from P (phosphorus).

  • source (str, optional) – If you used a KPOINTS_OPT file to use a second k-point mesh, you can provide a keyword argument source=”kpoints_opt” to use the k-points defined in that file instead of the one specified in KPOINTS.

Returns

Contains the energies at which the DOS was evaluated aligned to the Fermi energy and the total DOS or the spin-resolved DOS for spin-polarized calculations. If available and a selection is passed, the orbital resolved DOS for the selected orbitals is included.

Return type

dict

Examples

Select the p orbitals of the first atom in the POSCAR file:

>>> calc.dos.read(selection="1(p)")

Select the d orbitals of Mn, Co, and Fe:

>>> calc.dos.read("d(Mn, Co, Fe)")

Select the spin-up contribution of the first three atoms combined

>>> calc.dos.read("up(1:3)")
to_frame(selection=None)

Read the data into a pandas DataFrame.

Parameters
  • selection (str) –

    A string specifying the projection of the orbitals. There are three distinct possibilities:

    • To specify the atom, you can either use its element name (Si, Al, …) or its index as given in the input file (1, 2, …). For the latter option it is also possible to specify ranges (e.g. 1:4).

    • To select a particular orbital you can give a string (s, px, dxz, …) or select multiple orbitals by their angular momentum (s, p, d, f).

    • For the spin, you have the options up, down, or total.

    You separate multiple selections by commas or whitespace and can nest them using parenthesis, e.g. Sr(s, p) or s(up), p(down). The order of the selections does not matter, but it is case sensitive to distinguish p (angular momentum l = 1) from P (phosphorus).

  • source (str, optional) – If you used a KPOINTS_OPT file to use a second k-point mesh, you can provide a keyword argument source=”kpoints_opt” to use the k-points defined in that file instead of the one specified in KPOINTS.

Returns

Contains the energies at which the DOS was evaluated aligned to the Fermi energy and the total DOS or the spin-resolved DOS for spin-polarized calculations. If available and a selection is passed, the orbital resolved DOS for the selected orbitals is included.

Return type

pd.DataFrame

Examples

Select the p orbitals of the first atom in the POSCAR file:

>>> calc.dos.to_frame(selection="1(p)")

Select the d orbitals of Mn, Co, and Fe:

>>> calc.dos.to_frame("d(Mn, Co, Fe)")

Select the spin-up contribution of the first three atoms combined

>>> calc.dos.to_frame("up(1:3)")
to_image(*args, filename=None, **kwargs)

Read the data and generate an image writing to the given filename.

The filetype is automatically deduced from the filename; possible are common raster (png, jpg) and vector (svg, pdf) formats. If no filename is provided a default filename is deduced from the name of the class and the picture has png format.

Note that the filename must be a keyword argument, i.e., you explicitly need to write filename=”name_of_file” because the arguments are passed on to the plot() function. Please check the documentation of that function to learn which arguments are allowed.

to_plotly(selection=None)

Read the data and generate a plotly figure.

Parameters
  • selection (str) –

    A string specifying the projection of the orbitals. There are three distinct possibilities:

    • To specify the atom, you can either use its element name (Si, Al, …) or its index as given in the input file (1, 2, …). For the latter option it is also possible to specify ranges (e.g. 1:4).

    • To select a particular orbital you can give a string (s, px, dxz, …) or select multiple orbitals by their angular momentum (s, p, d, f).

    • For the spin, you have the options up, down, or total.

    You separate multiple selections by commas or whitespace and can nest them using parenthesis, e.g. Sr(s, p) or s(up), p(down). The order of the selections does not matter, but it is case sensitive to distinguish p (angular momentum l = 1) from P (phosphorus).

  • source (str, optional) – If you used a KPOINTS_OPT file to use a second k-point mesh, you can provide a keyword argument source=”kpoints_opt” to use the k-points defined in that file instead of the one specified in KPOINTS.

Returns

plotly figure containing the total DOS. If the calculation was spin polarized, the resulting DOS is spin resolved and the spin-down DOS is plotted towards negative values. If a selection the orbital resolved DOS is given for the specified projectors.

Return type

plotly.graph_objects.Figure

Examples

Select the p orbitals of the first atom in the POSCAR file:

>>> calc.dos.plot(selection="1(p)")

Select the d orbitals of Mn, Co, and Fe:

>>> calc.dos.plot("d(Mn, Co, Fe)")

Select the spin-up contribution of the first three atoms combined

>>> calc.dos.plot("up(1:3)")