question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

Nuclear repulsion energy is tough to find!

See original GitHub issue

What 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:closed
  • Created a year ago
  • Comments:5 (4 by maintainers)

github_iconTop GitHub Comments

1reaction
mrossinekcommented, May 31, 2022

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 the ElectronicEnergy.__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:

  • implement ElectronicEnergy.__str__ which should:
    • call super().__str__()
    • append information on self.nuclear_repulsion_energy and self.reference_energy
    • append information on self.orbital_energies
    • append information on self.kinetic and self.overlap
  • add a section to the Electronic Structure tutorial highlighting how to extract properties from the ElectronicStructureProblem:
    • showcase print(es_problem)
    • i.e. es_problem.get_property("some property name")
    • link to the properties tutorial for more information
0reactions
mrossinekcommented, Aug 18, 2022

A short update on this: the ElectronicStructureProblem and ElectronicEnergy 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 of nuclear_repulsion_energy, fixing this issue.

But I am unassigning and removing the good first issue label from here for the time being,

Read more comments on GitHub >

github_iconTop 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 >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found