e-eight / chime Goto Github PK
View Code? Open in Web Editor NEWTwo body matrix elements of chiral EFT operators in 3DHO basis.
License: MIT License
Two body matrix elements of chiral EFT operators in 3DHO basis.
License: MIT License
When generating charge radius files, the order defaults to:
lo, n2lo, n3lo, n4lo, nlo
This means that the only way to get the nlo file is to set the target as nlo explicitly; then you get a file for lo, n2lo, n3lo, n4lo, nlo, and an nlo_cumulative file which actually is the sum of all of them (so its really a 'n4lo_cumulative' file). Note also the n2lo and n4lo files have all 0 elements (n2lo because it has no contribution, n4lo since it isn't done yet for 2b).
Note that this means the cumulative files for >= n2lo miss the nlo correction (finite radius).
ChiME's LO elements seem to scale with Nmax, which my best guess is from a norm problem. For example, if I divide ChiME's LO r^2 deuteron channel elements by 4 and subtract ones giving the 'correct' deuteron r^2, I get
Note that the first ME is basically right, but not the others (which makes me think the RME treatment is fine, since the j-hat factors would be a constant scaling factor on all elements otherwise).
Regarding the factor 4, at LO the r operator in relative coordinates is taken to be the distance between nucleons, so half of this is the actual nuclear radius. We should make sure the r operator is consistent everywhere.
It appears a phase error has been introduced in 2641065.
relativecm-gen
input relcm.in
:
1 0 1 1
4 15.0
1.0
mm nlo 2
LENPICchi2bM11C_Nmax04_hw15.0_relcm.dat
Before:
# RELATIVE-CM LSJT
# version
# J0 g0 T0_min T0_max symmetry_phase_mode [P0=(-)^g0]
# Nmax
# T0 Nr' lr' Nc' lc' L' S' J' T' g' Nr lr Nc lc L S J T g JT-RME
1
1 0 1 1 0
4
1 1 1 1 1 0 0 0 0 0 1 1 1 1 0 1 1 1 0 -1.98073763e-01
1 1 1 1 1 0 0 0 0 0 1 1 3 1 0 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 3 1 1 1 0 1 1 1 0 -4.16307108e-02
1 1 1 3 1 0 0 0 0 0 1 1 1 1 0 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 3 1 0 1 1 1 0 -1.98073763e-01
1 1 1 3 1 0 0 0 0 0 3 1 1 1 0 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 1 1 0 1 1 1 0 -4.16307108e-02
1 3 1 1 1 0 0 0 0 0 1 1 3 1 0 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 3 1 1 1 0 1 1 1 0 -1.14243860e-01
1 1 1 1 1 0 0 0 0 0 2 2 2 2 1 0 1 1 0 -0.00000000e+00 <<<
1 1 1 3 1 0 0 0 0 0 2 2 2 2 1 0 1 1 0 -8.73909393e-02 <<<
1 3 1 1 1 0 0 0 0 0 2 2 2 2 1 0 1 1 0 -0.00000000e+00 <<<
1 1 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 1 1 3 1 1 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 3 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 3 1 1 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 3 1 1 1 1 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 3 1 1 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 3 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 1 1 1 1 2 1 1 1 0 2.07579903e-01
1 1 1 1 1 0 0 0 0 0 1 1 3 1 2 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 1 1 3 3 2 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 3 1 1 1 2 1 1 1 0 7.68070656e-02
1 1 1 1 1 0 0 0 0 0 3 3 1 1 2 1 1 1 0 -1.64704540e-01
1 1 1 3 1 0 0 0 0 0 1 1 1 1 2 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 3 1 2 1 1 1 0 2.07579903e-01
1 1 1 3 1 0 0 0 0 0 1 1 3 3 2 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 3 1 1 1 2 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 3 3 1 1 2 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 1 1 2 1 1 1 0 7.68070656e-02
1 3 1 1 1 0 0 0 0 0 1 1 3 1 2 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 3 3 2 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 3 1 1 1 2 1 1 1 0 1.27686495e-01
1 3 1 1 1 0 0 0 0 0 3 3 1 1 2 1 1 1 0 -5.08513652e-02
1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 1 0 0 -2.01831100e-01
1 0 0 0 0 0 0 0 1 0 0 0 2 0 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 2 0 0 0 0 1 1 0 0 2.05509302e-02
1 0 0 0 0 0 0 0 1 0 0 0 4 0 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 2 0 2 0 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 2 2 2 2 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 4 0 0 0 0 1 1 0 0 4.19919025e-02
After:
# RELATIVE-CM LSJT
# version
# J0 g0 T0_min T0_max symmetry_phase_mode [P0=(-)^g0]
# Nmax
# T0 Nr' lr' Nc' lc' L' S' J' T' g' Nr lr Nc lc L S J T g JT-RME
1
1 0 1 1 0
4
1 1 1 1 1 0 0 0 0 0 1 1 1 1 0 1 1 1 0 -1.98073763e-01
1 1 1 1 1 0 0 0 0 0 1 1 3 1 0 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 3 1 1 1 0 1 1 1 0 -4.16307108e-02
1 1 1 3 1 0 0 0 0 0 1 1 1 1 0 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 3 1 0 1 1 1 0 -1.98073763e-01
1 1 1 3 1 0 0 0 0 0 3 1 1 1 0 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 1 1 0 1 1 1 0 -4.16307108e-02
1 3 1 1 1 0 0 0 0 0 1 1 3 1 0 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 3 1 1 1 0 1 1 1 0 -1.14243860e-01
1 1 1 1 1 0 0 0 0 0 2 2 2 2 1 0 1 1 0 0.00000000e+00 <<<
1 1 1 3 1 0 0 0 0 0 2 2 2 2 1 0 1 1 0 8.73909393e-02 <<<
1 3 1 1 1 0 0 0 0 0 2 2 2 2 1 0 1 1 0 0.00000000e+00 <<<
1 1 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 1 1 3 1 1 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 3 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 3 1 1 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 3 1 1 1 1 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 3 1 1 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 3 1 1 1 1 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 1 1 1 1 2 1 1 1 0 2.07579903e-01
1 1 1 1 1 0 0 0 0 0 1 1 3 1 2 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 1 1 3 3 2 1 1 1 0 -0.00000000e+00
1 1 1 1 1 0 0 0 0 0 3 1 1 1 2 1 1 1 0 7.68070656e-02
1 1 1 1 1 0 0 0 0 0 3 3 1 1 2 1 1 1 0 -1.64704540e-01
1 1 1 3 1 0 0 0 0 0 1 1 1 1 2 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 1 1 3 1 2 1 1 1 0 2.07579903e-01
1 1 1 3 1 0 0 0 0 0 1 1 3 3 2 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 3 1 1 1 2 1 1 1 0 -0.00000000e+00
1 1 1 3 1 0 0 0 0 0 3 3 1 1 2 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 1 1 2 1 1 1 0 7.68070656e-02
1 3 1 1 1 0 0 0 0 0 1 1 3 1 2 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 1 1 3 3 2 1 1 1 0 -0.00000000e+00
1 3 1 1 1 0 0 0 0 0 3 1 1 1 2 1 1 1 0 1.27686495e-01
1 3 1 1 1 0 0 0 0 0 3 3 1 1 2 1 1 1 0 -5.08513652e-02
1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 1 0 0 -2.01831100e-01
1 0 0 0 0 0 0 0 1 0 0 0 2 0 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 2 0 0 0 0 1 1 0 0 2.05509302e-02
1 0 0 0 0 0 0 0 1 0 0 0 4 0 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 2 0 2 0 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 2 2 2 2 0 1 1 0 0 -0.00000000e+00
1 0 0 0 0 0 0 0 1 0 4 0 0 0 0 1 1 0 0 4.19919025e-02
Many of the other matrix elements differ in phase, but this is just the first few.
Compiling ChiME failed in:
We should consider generalizing the makefile and use config files for particular systems.
A declarative, efficient, and flexible JavaScript library for building user interfaces.
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. ๐๐๐
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google โค๏ธ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.