-
Notifications
You must be signed in to change notification settings - Fork 436
Bethe Salpeter equation for molecules #3308
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
Conversation
Added Bethe Salpeter on top of GW for the calculation of electronic excitations in molecules.
src/bse_iterative.F
Outdated
| WRITE (unit_nr, *) "z_space_energy_cutoff", z_space_energy_cutoff | ||
| END IF | ||
|
|
||
| !MG to del |
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.
Please remove this debug code and similar code below
src/bse_iterative.F
Outdated
| WRITE (unit_nr, *) "MG: Resetting dimension in i_iter=", i_iter | ||
| END IF | ||
| IF (unit_nr > 0) THEN | ||
| OPEN (10, file=res_filename, status="old", access="sequential", position="append") |
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.
Please ask CP2K for a unit number (close_file/open_file/get_unit_number from the module "cp_files") or use the provide ones (default_input_unit/default_output_unit from the module "machine")
src/bse_util.F
Outdated
| DO i = 1, nrow_local | ||
| DO j = 1, ncol_local | ||
| IF (ABS(fm%local_data(i, j)) > thresh) THEN | ||
| WRITE (*, "(A7,T10,I5,T20,I5,T30,F13.5)") header, row_indices(i), col_indices(j), & |
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.
Add suitable unit numbers
src/input_cp2k_mp2.F
Outdated
| lone_keyword_l_val=.TRUE.) | ||
| CALL section_add_keyword(section, keyword) | ||
| CALL keyword_release(keyword) | ||
| ! CALL keyword_create(keyword, __LOCATION__, name="BSE", & |
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.
Why do you have a keyword and a section? Remove it if it is obsolete and not expected to be reused.
fstein93
left a comment
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.
I have left a few formal comments. The CI will complain about most of them anyways.
|
Thank you for the review! I incorporated all your comments in the new commit. |
|
There is still the |
Added Bethe Salpeter on top of GW for the calculation of electronic excitations in molecules.