Nuclear repulsion energy is tough to find!
See original GitHub issueWhat is the expected enhancement?
It’s pretty difficult to find the nuclear repulsion energy (and probably some other interesting properties) if one doesn’t already know where to look:
# es_problem is an `ElectronicStructureProblem`
es_problem.grouped_property_transformed.get_property("ElectronicEnergy").nuclear_repulsion_energy
It would be great it the tutorials would mention where they are, e.g. the ground state solver tutorial.
(Originated from this Slack convo)
Issue Analytics
- State:
- Created a year ago
- Comments:5 (4 by maintainers)
Top Results From Across the Web
Hartree Fock Program in Python Part 5 - Nuclear ... - YouTube
In this video, we add the nuclear nuclear repulsion energy term into our HF program. This code is not intended to be an...
Read more >What is Fusion, and Why Is It So Difficult to Achieve? | IAEA
Once the nuclei overcome this repulsion and come within a very close range of each other, the attractive nuclear force between them will ......
Read more >Basic ab initio - Outline of a calculation
Total energy = nuclear repulsion + electronic. The total energy is just evaluated by adding the nuclear repulsion energy to the electronic energy....
Read more >Coulomb Repulsion - an overview ... - ScienceDirect.com
Because of the Coulomb repulsion, too many protons relative to neutrons is bad for nuclear stability. That is why there is no such...
Read more >7.2: Shielding and Effective Nuclear Charge
For an atom or an ion with only a single electron, we can calculate the potential energy by considering only the electrostatic attraction...
Read more >Top Related Medium Post
No results found
Top Related StackOverflow Question
No results found
Troubleshoot Live Code
Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start FreeTop Related Reddit Thread
No results found
Top Related Hackernoon Post
No results found
Top Related Tweet
No results found
Top Related Dev.to Post
No results found
Top Related Hashnode Post
No results found
Top GitHub Comments
You are right! The properties tutorial should cover this but it does not specifically talk about the nuclear repulsion energy. In fact, the string-representation of the
ElectronicEnergy
class should be improved to also include this attribute. This requires implementing theElectronicEnergy.__str__
method.I do not think that this information should go into the GS tutorial since what you are suggesting is specifically designed for extracting the information before running a solver. Instead, a logical link should be made between the electronic structure tutorial and the aforementioned properties tutorial.
To summarize, fixing this issue involves the following tasks:
ElectronicEnergy.__str__
which should:super().__str__()
self.nuclear_repulsion_energy
andself.reference_energy
self.orbital_energies
self.kinetic
andself.overlap
ElectronicStructureProblem
:print(es_problem)
es_problem.get_property("some property name")
A short update on this: the
ElectronicStructureProblem
andElectronicEnergy
are undergoing a major refactoring, of which #796 does a lot of the intended changes. This means, that my outline above is no longer up-to-date. I will make sure, that as part of this refactoring, we will significantly improve the visibility ofnuclear_repulsion_energy
, fixing this issue.But I am unassigning and removing the
good first issue
label from here for the time being,