-
Notifications
You must be signed in to change notification settings - Fork 22
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
OptaDOS photoemission module #38
base: develop
Are you sure you want to change the base?
Changes from 1 commit
File filter
Filter by extension
Conversations
Jump to
Diff view
Diff view
There are no files selected for viewing
Large diffs are not rendered by default.
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,29 @@ | ||
%BLOCK LATTICE_ABC | ||
2.517885 2.517885 50 | ||
90.000000 90.000000 90.000000 | ||
%ENDBLOCK LATTICE_ABC | ||
|
||
%BLOCK POSITIONS_FRAC | ||
Cu -0.250000 0.250000 0.266523 | ||
Cu 0.250000 -0.250000 0.231863 | ||
Cu -0.250000 0.250000 0.196014 | ||
Cu 0.250000 -0.250000 0.160204 | ||
Cu -0.250000 0.250000 0.124622 | ||
Cu 0.250000 -0.250000 0.089005 | ||
Cu -0.250000 0.250000 0.053413 | ||
Cu 0.250000 -0.250000 0.017804 | ||
Cu -0.250000 0.250000 -0.017804 | ||
Cu 0.250000 -0.250000 -0.053413 | ||
Cu -0.250000 0.250000 -0.089005 | ||
Cu 0.250000 -0.250000 -0.124622 | ||
Cu -0.250000 0.250000 -0.160204 | ||
Cu 0.250000 -0.250000 -0.196014 | ||
Cu -0.250000 0.250000 -0.231863 | ||
Cu 0.250000 -0.250000 -0.266523 | ||
%ENDBLOCK POSITIONS_FRAC | ||
|
||
symmetry_generate | ||
|
||
fix_all_cell : true | ||
|
||
kpoint_mp_grid 32 32 1 |
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,20 @@ | ||
TASK : photoemission | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Might be worth commenting this file more, to explain what the new keywords do here. |
||
work_function : 4.556 | ||
s_area : 6.339744873 | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. What does s_area mean? Slab area? If so, it would be better to call it slab_area to match slab_volume. |
||
slab_volume : 190.942338 | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Why isn't the slab volume calculated from the lattice vectors ? |
||
elec_field : 0.001 | ||
imfp_const : 19.0 | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. I think I can work out all of the other keywords, this one seems very mysterious. Is there a better name for it? There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. I changed s_area to surface_area, that should be more clear. I will checked the name of the keywords to make them more clear. There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. There's an argument that photoemission keywords relevant to that module /only/ might be prefixed by |
||
JDOS_SPACING : 0.1 | ||
JDOS_MAX_ENERGY : 10 | ||
BROADENING : linear | ||
OPTICS_GEOM : unpolar | ||
optics_qdir : 1 1.000 1.0000 | ||
photon_energy : 4.7 | ||
linear_smearing : 0.026 | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Why do we have both smearing set here? |
||
fixed_smearing : 0.026 | ||
optics_intraband : true | ||
photo_model : 3step | ||
momentum : crystal | ||
hybrid_linear : true | ||
temp : 300 | ||
photo_output = be | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. What is be? Can we have a better name for this please? There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. I changed it to binding_energy. It is a keyword to generate a plot of the photoemission intensities as a function of the binding energy. |
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,17 @@ | ||
TASK : SPECTRAL | ||
SPECTRAL_TASK : ALL | ||
PDOS_CALCULATE_WEIGHTS : TRUE | ||
xc_functional : PBEsol | ||
opt_strategy : speed | ||
smearing_width : 300 K | ||
cut_off_energy : 600 eV | ||
elec_energy_tol : 1e-8 eV | ||
PERC_EXTRA_BANDS : 100 | ||
MAX_SCF_CYCLES : 1000 | ||
FIX_OCCUPANCY : FALSE | ||
METALS_METHOD : DM | ||
SPIN_POLARIZED : FALSE | ||
WRITE_FORMATTED_POTENTIAL : TRUE | ||
WRITE_FORMATTED_DENSITY : TRUE | ||
MIXING_SCHEME : Pulay | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Are all of these parameters necessary? It would be helpful not to have red-herrings in here. There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. You are right, they are not a requirement. They has to do more with other results in my work. |
||
|
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,46 @@ | ||
# CELL file written 10:16:14 (GMT+2.0) 18th July 2020 from run Si2 | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. This shouldn't be checked in. |
||
|
||
%BLOCK lattice_cart | ||
ANG | ||
0.00000000000000 2.73000000000000 2.73000000000000 | ||
2.73000000000000 0.00000000000000 2.73000000000000 | ||
2.73000000000000 2.73000000000000 0.00000000000000 | ||
%ENDBLOCK lattice_cart | ||
|
||
%BLOCK positions_frac | ||
Si 0.500000000000000 0.500000000000000 0.500000000000000 | ||
Si 0.750000000000000 0.750000000000000 0.750000000000000 | ||
%ENDBLOCK positions_frac | ||
|
||
FIX_COM : true | ||
|
||
%BLOCK species_pot | ||
Si 3|1.8|5|6|7|30:31:32 | ||
%ENDBLOCK species_pot | ||
|
||
%BLOCK species_lcao_states | ||
Si 2 | ||
%ENDBLOCK species_lcao_states | ||
|
||
SYMMETRY_TOL : 0.001000 | ||
|
||
%BLOCK symmetry_ops | ||
# Symm. op. 1 E | ||
1.000000000000000 0.000000000000000 0.000000000000000 | ||
0.000000000000000 1.000000000000000 0.000000000000000 | ||
0.000000000000000 0.000000000000000 1.000000000000000 | ||
0.000000000000000 0.000000000000000 0.000000000000000 | ||
%ENDBLOCK symmetry_ops | ||
|
||
kpoint_mp_grid : 10 10 10 | ||
|
||
bs_kpoint_mp_grid : 10 10 10 | ||
|
||
optics_kpoint_mp_grid : 10 10 10 | ||
|
||
elnes_kpoint_mp_grid : 10 10 10 | ||
|
||
magres_kpoint_mp_grid : 10 10 10 | ||
|
||
phonon_kpoint_mp_grid : 10 10 10 | ||
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is this important for the calculation? If not, please can we have a minimal working example ?