CP2K: default Hirshfeld charges (Gaussian shape functions) are unphysical for ionic solids; generate SHAPE_FUNCTION DENSITY or label them
Nobody has claimed this yet.
Assessment
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Newbie friendliness
- 42/100
Research direction
Start by tracing how ECCE builds CP2K inputs when charge analysis is requested, then find how the Viewer labels or displays CP2K Hirshfeld charges. The issue identifies the desired density-based SHAPE_FUNCTION option and an approximate label or hiding behavior for default Gaussian-shape results. Done means charge requests use the specified density-based analysis and default-run charges are labeled or hidden; the payload names no files or tests.
Written by the indexing model from the issue text.
Description
Andy's CP2K run 03_SOMS_ord1_final_PBE-D3BJ_700Ry_rel100_forces (sodium niobate, GTH-q9 Na, PBE-D3BJ): CP2K's own output gives Mulliken Na +0.82 but Hirshfeld Na -1.02, Nb -0.41, O +0.24 (Hirshfeld table: Ref Charge 9.000, Population 10.02). The input has no &HIRSHFELD section, so CP2K used its default analysis (SHAPE_FUNCTION GAUSSIAN, radii from tables), which is a crude promolecule and fails for ions with semicore pseudopotentials. ECCE parses the table correctly; the numbers are CP2K's. Wanted: when ECCE's CP2K input asks for charges, write &PRINT &HIRSHFELD SHAPE_FUNCTION DENSITY &END (atomic reference densities; optionally SELF_CONSISTENT T for iterative Hirshfeld), and in the Viewer label CP2K Hirshfeld charges from a default (Gaussian-shape) run as approximate, or hide them.
- Dominant language
- C
- Stars
- 10
- Forks
- 8
- Avg merge
- 2h 39m
- Merged PRs (30d)
- 9
Getting set up
This project ships no dev container, Dockerfile or contributing guide, so setting up is up to you: start from its README, and see our first-contribution guide for the general steps.
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
More from FriendsofECCE/ECCE
-
Difficulty 2/5 1-3 hours Newbie friendliness 70/100
FriendsofECCE/ECCE#336 · 1 comment ·
-
Difficulty 4/5 3-5 days Newbie friendliness 55/100
FriendsofECCE/ECCE#712 ·
-
Launch Info: wall time left for a running job (requested while queued, used once ended)Possibly taken @ohlincha claimed this today. Open
FriendsofECCE/ECCE#711 · 1 assignee ·
-
Recipes: create a workflow's linked boxes from a recipe, with the level of theory chosen in a dialogOpen
Difficulty 5/5 Over a week Newbie friendliness 20/100
FriendsofECCE/ECCE#710 · 2 comments ·
-
Difficulty 5/5 Over a week Newbie friendliness 25/100
FriendsofECCE/ECCE#709 ·
All issues in FriendsofECCE/ECCE
Similar issues
-
vim.filetype.match() prefixes a URL buffer name with the cwd before passing it to filetype functionsOpen
Difficulty 2/5 1-3 hours Newbie friendliness 73/100
Maintainers usually reply within 1 day
-
compilation errorOpen
Difficulty 2/5 1-3 hours Newbie friendliness 68/100
COVESA/dlt-daemon#925 · 1 reaction ·
-
stream_curl: retry_count never resets after a successful reconnect, so long sequential streams die after 5 dropsPossibly taken A pull request linked to this issue is open or already merged. Open
Difficulty 2/5 1-3 hours Newbie friendliness 70/100
mpv-player/mpv#18728 ·
Maintainers usually reply within 1 day
-
followup Prio 📌
Difficulty 2/5 1-3 hours Newbie friendliness 66/100
Maintainers usually reply within 1 day
-
Difficulty 2/5 1-3 hours Newbie friendliness 72/100
cyrusimap/cyrus-imapd#6640 ·
Maintainers usually reply within 2 days