Git Product home page Git Product logo

Comments (4)

benroberts999 avatar benroberts999 commented on July 27, 2024

Some of the Breit corrections to Hyperfine constants differ significantly from Derevianko 10.1103/PhysRevA.65.012106.
The Breit corrections at Hartree-Fock level agree perfectly.

This is in contrast to E1, where Breit corrections at HF level match exactly, and agree very well at RPA level. Also, Breit corrections to E1 match exactly with diagram and TDHF method for RPA.

The difference is at most a factor of two for the RPA correction part; the result is only noticable when there is large cancellation between HF and RPA contributions, so in the end, doesn't seem like major issue.

HF ampsci Dzuba Derevianko
6s+ 0.011 0.010 0.011
6p- -0.680 -0.680 -0.680
6p+ -0.062 -0.062 -0.060
5d- 0.099 0.099 0.099
5d+ 0.053 0.053
HF+RPA ampsci Dzuba Derevianko
6s+ 4.383 4.010 4.111
6p- 0.012 -0.270 -0.29
6p+ 0.114 0.072 0.04
5d- 0.266 0.233 0.199
5d+ -0.169 -0.183
RPA only ampsci Dzuba Derevianko
6s+ 4.372 4.000 4.1
6p- 0.692 0.410 0.39
6p+ 0.176 0.135 0.1
5d- 0.168 0.134 0.1
5d+ -0.221 -0.236

Full code outputs:
rpa-breit.txt

from ampsci.

benroberts999 avatar benroberts999 commented on July 27, 2024

Note: get exactly the same Breit corrections with diagram and "TDHF Basis" methods

  • This essentially means issue is not in RPA equations, which are implemented very differently between these two methods
  • Possibly issue in Breit equations; or, possibly no issue at all!
  • "Basis" method is equivalent to TDHF, but instead of solving TDHF equations, expands solutions using a basis
    • nb: "basis" method is extremely slow, and should only be used for testing
HF+RPA RPA Only
Diagram Basis Diagram Basis
6s+ 4.383 4.383 4.372 4.372
6p- 0.012 0.015 0.692 0.694
6p+ 0.114 0.115 0.176 0.177
5d- 0.266 0.269 0.168 0.170
5d+ -0.169 -0.167 -0.221 -0.220

Full code input/output: rpa-breit-basis.txt

I'll leave this issue up for now as a note that this probably requires extra testing, comparison with other calculations etc.
But essentially, it seems there is not an error in ampsci

from ampsci.

benroberts999 avatar benroberts999 commented on July 27, 2024

As an extra note - also less than perfect agreement for two-body Breit correction to energies at second-order compared to Derevianko.
This correction is complicated, and depends strongly on Brueckner convergence, and overall correction still agrees well. Therefore, not really major issue, but something that would be good to check:

(Nb: it's not non-linear in Breit terms which make the difference, which can be easily check by scaling Breit potential by small number)

Breit corrections at Hartree-Fock and Brueckner levels, accounting for one-body (1) and two-body (2) Breit contributions. Derevianko results 10.1103/PhysRevA.65.012106 in italics.

HF BO(1) BO(2) Total
6s_1/2 3.2 3.2 -6.1 -5.0 -0.3 -0.8 -3.2 -2.6
7s_1/2 1.1 1.1 -1.0 -1.1 -0.1 -0.3 0.02 -0.26
6p_1/2 7.5 7.5 -0.2 -0.1 -0.4 -0.3 6.9 7.1
7p_1/2 2.7 2.7 -0.1 -0.1 -0.1 -0.1 2.5 2.5
6p_3/2 2.9 2.9 -2.3 -1.8 0.02 -0.25 0.59 0.84
7p_3/2 1.0 1.0 -0.6 -0.6 0.01 -0.09 0.44 0.38
5d_3/2 -10.2 -10.2 -15.8 -12 -1.2 -0.4 -27.1 -22
5d_5/2 -11.7 -11.8 -19.2 -14 0.4 -0.3 -30.6 -26

Full code output: 240125-Breit-BO2.txt

from ampsci.

benroberts999 avatar benroberts999 commented on July 27, 2024

Note: issue is now resolved (at least, the RPA part) - issue was with diagram method, TDHF works very well.

See commits: fefa5bc
and: 39b721e

Updated output:

HF ampsci Dzuba Derevianko
6s+ 0.011 0.010 0.011
6p- -0.680 -0.680 -0.680
6p+ -0.062 -0.062 -0.060
5d- 0.099 0.099 0.099
5d+ 0.053 0.053
HF+RPA ampsci Dzuba Derevianko
6s+ 4.025 4.010 4.111
6p- -0.237 -0.270 -0.29
6p+ 0.075 0.072 0.04
5d- 0.233 0.233 0.199
5d+ -0.179 -0.183
RPA only ampsci Dzuba Derevianko
6s+ 4.014 4.000 4.1
6p- 0.442 0.410 0.39
6p+ 0.137 0.135 0.1
5d- 0.134 0.134 0.1
5d+ -0.232 -0.236

Full code outputs: rpa-breit-updated.txt

from ampsci.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.