|
12 | 12 | "metadata": {}, |
13 | 13 | "source": [ |
14 | 14 | "First, we calculate the overlap `<ψ_i(0)|ψ_f(Q)>` and its derivative. \n", |
| 15 | + "`ψ_i(0)` is a perturbed band-edge state, and `ψ_f(Q)` is a localized defect state at the configuration coordinate `Q` which can be calculated using `get_del_Q.py`.\n", |
15 | 16 | "The electron-phonon coupling matrix element is given by \n", |
16 | 17 | "`W_if = (ε_f - ε_i) d <ψ_i(0)|ψ_f(Q)> / dQ`.\n", |
17 | 18 | "\n", |
|
36 | 37 | "metadata": {}, |
37 | 38 | "source": [ |
38 | 39 | "```bash\n", |
39 | | - "$ python get_wf.py -d 288 -b 291 -D DISP_000 -i DISP_-02 DISP_-01 DISP_001 DISP_002\n", |
| 40 | + "$ python get_wf.py -d <band_index of ψ_i> -b <band_index of ψ_f> -D <path to the initial geometry> -i <path to the deformed geometry 1> [<path to the deformed geometry 2> ...]\n", |
40 | 41 | "\n", |
41 | 42 | "GRID ENCUT 918.2911873931497\n", |
42 | 43 | "finished making projector list\n", |
|
68 | 69 | "source": [ |
69 | 70 | "## W_if calculation\n", |
70 | 71 | "\n", |
71 | | - "Use `get_del_Q.py` to get a distance between structures (`ΔQ`).\n", |
| 72 | + "Use `get_del_Q.py` to get a distance between structures (`ΔQ`) in `amu^0.5Å`. `ΔQ` is a collective variable of the mass-weighted deformation `ΔQ = ΣmΔ𝐑`.\n", |
72 | 73 | "\n", |
73 | | - "`ϵ_i` and `ϵ_f` are KS eigenvalues of the initla and final states.\n" |
| 74 | + "`ϵ_i` and `ϵ_f` are KS eigenvalues of the initial (`ψ_i(0)`) and final states (`ψ_f(Q_f)`).\n" |
74 | 75 | ] |
75 | 76 | }, |
76 | 77 | { |
77 | 78 | "cell_type": "code", |
78 | | - "execution_count": 15, |
| 79 | + "execution_count": 1, |
79 | 80 | "metadata": {}, |
80 | 81 | "outputs": [ |
81 | 82 | { |
|
89 | 90 | "source": [ |
90 | 91 | "# using Plots\n", |
91 | 92 | "# \n", |
92 | | - "ΔQ = 0.18059526482789678 # ΔQ between DISP_000 and DISP_001\n", |
| 93 | + "ΔQ = 0.18059526482789678 # ΔQ between DISP_000 and DISP_001 in [amu^0.5Å]\n", |
93 | 94 | "ϵ_f = 5.405963 # unoccupied trap\n", |
94 | 95 | "ϵ_i = 6.227692 # CBM\n", |
95 | 96 | "\n", |
|
104 | 105 | "println(\"W_if_n = \", W_if_n)\n", |
105 | 106 | "# plot(overlap)" |
106 | 107 | ] |
107 | | - }, |
108 | | - { |
109 | | - "cell_type": "code", |
110 | | - "execution_count": null, |
111 | | - "metadata": {}, |
112 | | - "outputs": [], |
113 | | - "source": [] |
114 | 108 | } |
115 | 109 | ], |
116 | 110 | "metadata": { |
|
0 commit comments