Luminescent properties (1D-CCM)

This example shows how to post-process the results of a LumiWork following a one-dimensional configuration coordinate model (1D-CCM). Based on NV- center in diamond (64 atoms supercell).

from abipy.lumi.deltaSCF import DeltaSCF
import abipy.data as abidata
import abipy.abilab as abilab

Read the 4 points netcdf file produced by a LumiWork

NV_center = DeltaSCF.from_four_points_file([
    abidata.ref_file("relaxed_gs_out_GSR.nc"),
    abidata.ref_file("unrelaxed_ex_out_GSR.nc"),
    abidata.ref_file("relaxed_ex_out_GSR.nc"),
    abidata.ref_file("unrelaxed_gs_out_GSR.nc")
])

To draw the configuration coordinates diagram.

plot lumi 1D nv center

To plot the luminescence lineshape following the one effective phonon mode model at 0K

NV_center.plot_lineshape_1D_zero_temp(
    energy_range=[0.5,3],
    max_m=20,
    phonon_width=0.02,
    with_omega_cube=True,
    normalized='Sum');
plot lumi 1D nv center

To get a panda dataframe with the main results:

To plot the magnitude of the displacements/forces, from the N atom.

NV_center.plot_delta_R_distance(defect_symbol="N");
NV_center.plot_delta_F_distance(defect_symbol="N");
  • plot lumi 1D nv center
  • plot lumi 1D nv center

To visualise displacements in 3D

plot lumi 1D nv center

To plot the four band structures associated to each point nscf_files are obtained in a LumiWork by activating the BS computations

nscf_files =[
    abidata.ref_file("relaxed_gs_nscf_GSR.nc"),
    abidata.ref_file("unrelaxed_ex_nscf_GSR.nc"),
    abidata.ref_file("relaxed_ex_nscf_GSR.nc"),
    abidata.ref_file("unrelaxed_gs_nscf_GSR.nc")
]

NV_center.plot_four_BandStructures(nscf_files, ylims=[-4,4]);
$A_g$, $A_g^*$, $A_e^*$, $A_e$

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

Gallery generated by Sphinx-Gallery