代码拉取完成,页面将自动刷新
nohup: ignoring input
p0n1 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p0n1
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 0 pairs of protons, and 1 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=0 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -4.751733, it took 0 s elapsed time
iter = 2, f = -5.045849, it took 1 s elapsed time
iter = 3, f = -5.606877, it took 0 s elapsed time
iter = 4, f = -6.559887, it took 0 s elapsed time
iter = 5, f = -7.605635, it took 0 s elapsed time
iter = 6, f = -7.762406, it took 1 s elapsed time
iter = 7, f = -8.596781, it took 0 s elapsed time
iter = 8, f = -9.193417, it took 1 s elapsed time
iter = 9, f = -9.204772, it took 0 s elapsed time
iter = 10, f = -9.448427, it took 0 s elapsed time
iter = 11, f = -9.880032, it took 1 s elapsed time
iter = 12, f = -10.496259, it took 0 s elapsed time
iter = 13, f = -11.013029, it took 0 s elapsed time
iter = 14, f = -11.098524, it took 1 s elapsed time
iter = 15, f = -11.105704, it took 0 s elapsed time
iter = 16, f = -11.251052, it took 1 s elapsed time
iter = 17, f = -11.472975, it took 0 s elapsed time
n=132, Minimum=-11.681159 found after 18 iterations
************************************************
Found the Minimum!
************************************************
iter = 18, f = -11.681159, it took 0 s elapsed time
before normalization, pove = 1 nove = 30.019
after normalization, pove = 1 nove = 1
---------- variation ends, it took 6 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -11.6812
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 6 s.
p0n1 variation ends
p1n1 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p1n1
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 1 pairs of protons, and 1 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -18.309542, it took 0 s elapsed time
iter = 2, f = -18.861803, it took 1 s elapsed time
iter = 3, f = -19.929807, it took 0 s elapsed time
iter = 4, f = -21.840212, it took 1 s elapsed time
iter = 5, f = -24.363021, it took 0 s elapsed time
iter = 6, f = -24.698510, it took 1 s elapsed time
iter = 7, f = -25.093332, it took 1 s elapsed time
iter = 8, f = -26.574568, it took 1 s elapsed time
iter = 9, f = -28.288623, it took 0 s elapsed time
iter = 10, f = -28.653653, it took 1 s elapsed time
iter = 11, f = -28.755030, it took 1 s elapsed time
iter = 12, f = -29.386270, it took 0 s elapsed time
iter = 13, f = -30.111195, it took 1 s elapsed time
iter = 14, f = -30.268877, it took 1 s elapsed time
iter = 15, f = -30.965819, it took 0 s elapsed time
iter = 16, f = -31.903622, it took 1 s elapsed time
iter = 17, f = -32.289647, it took 1 s elapsed time
iter = 18, f = -32.313511, it took 1 s elapsed time
iter = 19, f = -32.683732, it took 0 s elapsed time
iter = 20, f = -33.296518, it took 1 s elapsed time
iter = 21, f = -34.090746, it took 0 s elapsed time
iter = 22, f = -34.630807, it took 1 s elapsed time
iter = 23, f = -34.652011, it took 1 s elapsed time
iter = 24, f = -35.352565, it took 0 s elapsed time
iter = 25, f = -36.083737, it took 1 s elapsed time
iter = 26, f = -36.218859, it took 1 s elapsed time
n=132, Minimum=-36.599254 found after 27 iterations
************************************************
Found the Minimum!
************************************************
iter = 27, f = -36.599254, it took 1 s elapsed time
before normalization, pove = 70.4542 nove = 67.8365
after normalization, pove = 1 nove = 1
---------- variation ends, it took 19 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -36.5993
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 19 s.
p1n1 variation ends
p0n2 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p0n2
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 0 pairs of protons, and 2 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=0 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -9.766303, it took 0 s elapsed time
iter = 2, f = -10.485808, it took 0 s elapsed time
iter = 3, f = -11.832961, it took 1 s elapsed time
iter = 4, f = -13.842860, it took 0 s elapsed time
iter = 5, f = -15.028157, it took 0 s elapsed time
iter = 6, f = -15.028676, it took 1 s elapsed time
iter = 7, f = -15.082838, it took 0 s elapsed time
iter = 8, f = -15.189658, it took 1 s elapsed time
iter = 9, f = -15.397181, it took 0 s elapsed time
iter = 10, f = -15.787120, it took 0 s elapsed time
iter = 11, f = -16.463844, it took 1 s elapsed time
iter = 12, f = -17.414397, it took 0 s elapsed time
iter = 13, f = -18.058046, it took 0 s elapsed time
iter = 14, f = -18.059949, it took 1 s elapsed time
iter = 15, f = -18.207637, it took 0 s elapsed time
iter = 16, f = -18.483303, it took 0 s elapsed time
iter = 17, f = -18.953994, it took 1 s elapsed time
iter = 18, f = -19.582909, it took 0 s elapsed time
iter = 19, f = -19.891457, it took 0 s elapsed time
iter = 20, f = -19.901256, it took 1 s elapsed time
iter = 21, f = -20.066388, it took 0 s elapsed time
iter = 22, f = -20.355645, it took 1 s elapsed time
iter = 23, f = -20.773961, it took 0 s elapsed time
iter = 24, f = -21.080943, it took 0 s elapsed time
iter = 25, f = -21.082523, it took 1 s elapsed time
iter = 26, f = -21.173589, it took 0 s elapsed time
iter = 27, f = -21.337893, it took 0 s elapsed time
iter = 28, f = -21.598483, it took 1 s elapsed time
iter = 29, f = -21.883619, it took 0 s elapsed time
iter = 30, f = -21.936358, it took 1 s elapsed time
iter = 31, f = -22.273222, it took 0 s elapsed time
iter = 32, f = -22.531925, it took 0 s elapsed time
iter = 33, f = -22.534542, it took 1 s elapsed time
iter = 34, f = -22.608407, it took 0 s elapsed time
iter = 35, f = -22.723531, it took 1 s elapsed time
iter = 36, f = -22.837996, it took 0 s elapsed time
iter = 37, f = -22.846734, it took 0 s elapsed time
iter = 38, f = -22.923148, it took 1 s elapsed time
n=132, Minimum=-23.004985 found after 39 iterations
************************************************
Found the Minimum!
************************************************
iter = 39, f = -23.004985, it took 0 s elapsed time
before normalization, pove = 1 nove = 409.352
after normalization, pove = 1 nove = 1
---------- variation ends, it took 14 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -23.005
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 14 s.
p0n2 variation ends
p1n2 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p1n2
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 1 pairs of protons, and 2 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -31.339486, it took 0 s elapsed time
iter = 2, f = -32.159113, it took 1 s elapsed time
iter = 3, f = -33.720067, it took 0 s elapsed time
iter = 4, f = -36.242268, it took 1 s elapsed time
iter = 5, f = -38.370091, it took 0 s elapsed time
iter = 6, f = -38.448205, it took 2 s elapsed time
iter = 7, f = -39.693140, it took 0 s elapsed time
iter = 8, f = -41.503639, it took 1 s elapsed time
iter = 9, f = -42.554174, it took 0 s elapsed time
iter = 10, f = -42.570432, it took 1 s elapsed time
iter = 11, f = -42.929890, it took 1 s elapsed time
iter = 12, f = -43.597596, it took 1 s elapsed time
iter = 13, f = -44.711772, it took 0 s elapsed time
iter = 14, f = -46.058718, it took 1 s elapsed time
iter = 15, f = -46.308448, it took 0 s elapsed time
iter = 16, f = -46.453316, it took 1 s elapsed time
iter = 17, f = -47.599158, it took 1 s elapsed time
iter = 18, f = -49.104626, it took 0 s elapsed time
iter = 19, f = -49.773353, it took 1 s elapsed time
iter = 20, f = -49.802230, it took 1 s elapsed time
iter = 21, f = -50.201542, it took 1 s elapsed time
iter = 22, f = -50.833445, it took 0 s elapsed time
iter = 23, f = -51.490903, it took 1 s elapsed time
iter = 24, f = -51.568404, it took 1 s elapsed time
iter = 25, f = -52.133335, it took 0 s elapsed time
iter = 26, f = -52.629071, it took 1 s elapsed time
iter = 27, f = -52.641270, it took 1 s elapsed time
iter = 28, f = -52.839521, it took 1 s elapsed time
iter = 29, f = -53.119998, it took 0 s elapsed time
iter = 30, f = -53.266537, it took 1 s elapsed time
iter = 31, f = -53.276713, it took 1 s elapsed time
iter = 32, f = -53.379609, it took 1 s elapsed time
iter = 33, f = -53.515497, it took 0 s elapsed time
iter = 34, f = -53.527615, it took 1 s elapsed time
iter = 35, f = -53.562659, it took 1 s elapsed time
n=132, Minimum=-53.653775 found after 36 iterations
************************************************
Found the Minimum!
************************************************
iter = 36, f = -53.653775, it took 0 s elapsed time
before normalization, pove = 11.012 nove = 375.873
after normalization, pove = 1 nove = 1
---------- variation ends, it took 25 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -53.6538
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 25 s.
p1n2 variation ends
p2n2 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p2n2
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 2 pairs of protons, and 2 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -48.757851, it took 1 s elapsed time
iter = 2, f = -49.808821, it took 0 s elapsed time
iter = 3, f = -51.855828, it took 1 s elapsed time
iter = 4, f = -55.358214, it took 0 s elapsed time
iter = 5, f = -58.613705, it took 1 s elapsed time
iter = 6, f = -58.712967, it took 1 s elapsed time
iter = 7, f = -59.989070, it took 1 s elapsed time
iter = 8, f = -61.881989, it took 0 s elapsed time
iter = 9, f = -63.178934, it took 1 s elapsed time
iter = 10, f = -63.179039, it took 1 s elapsed time
iter = 11, f = -63.209594, it took 1 s elapsed time
iter = 12, f = -63.270404, it took 0 s elapsed time
iter = 13, f = -63.390802, it took 1 s elapsed time
iter = 14, f = -63.626558, it took 0 s elapsed time
iter = 15, f = -64.076718, it took 1 s elapsed time
iter = 16, f = -64.883683, it took 0 s elapsed time
iter = 17, f = -66.087182, it took 1 s elapsed time
iter = 18, f = -66.930417, it took 1 s elapsed time
iter = 19, f = -66.930528, it took 1 s elapsed time
iter = 20, f = -66.962138, it took 0 s elapsed time
iter = 21, f = -67.024985, it took 1 s elapsed time
iter = 22, f = -67.149164, it took 1 s elapsed time
iter = 23, f = -67.391281, it took 0 s elapsed time
iter = 24, f = -67.849284, it took 1 s elapsed time
iter = 25, f = -68.652870, it took 0 s elapsed time
iter = 26, f = -69.789841, it took 1 s elapsed time
iter = 27, f = -70.452671, it took 0 s elapsed time
iter = 28, f = -70.460176, it took 2 s elapsed time
iter = 29, f = -70.715852, it took 0 s elapsed time
iter = 30, f = -71.181908, it took 1 s elapsed time
iter = 31, f = -71.934126, it took 0 s elapsed time
iter = 32, f = -72.777076, it took 1 s elapsed time
iter = 33, f = -72.935691, it took 1 s elapsed time
iter = 34, f = -73.986910, it took 1 s elapsed time
iter = 35, f = -75.001451, it took 0 s elapsed time
iter = 36, f = -75.076035, it took 1 s elapsed time
iter = 37, f = -75.882879, it took 1 s elapsed time
iter = 38, f = -76.962506, it took 1 s elapsed time
iter = 39, f = -77.503377, it took 0 s elapsed time
iter = 40, f = -77.521375, it took 1 s elapsed time
iter = 41, f = -77.868477, it took 1 s elapsed time
iter = 42, f = -78.434618, it took 1 s elapsed time
iter = 43, f = -79.109776, it took 0 s elapsed time
iter = 44, f = -79.136091, it took 1 s elapsed time
iter = 45, f = -79.290423, it took 1 s elapsed time
iter = 46, f = -79.994440, it took 1 s elapsed time
iter = 47, f = -80.351263, it took 0 s elapsed time
iter = 48, f = -80.373914, it took 1 s elapsed time
iter = 49, f = -80.572034, it took 1 s elapsed time
iter = 50, f = -80.767103, it took 0 s elapsed time
iter = 51, f = -80.776506, it took 2 s elapsed time
iter = 52, f = -80.849755, it took 0 s elapsed time
iter = 53, f = -80.921292, it took 1 s elapsed time
iter = 54, f = -80.923968, it took 1 s elapsed time
n=132, Minimum=-80.945895 found after 55 iterations
************************************************
Found the Minimum!
************************************************
iter = 55, f = -80.945895, it took 1 s elapsed time
before normalization, pove = 1165.87 nove = 673.537
after normalization, pove = 1 nove = 1
---------- variation ends, it took 40 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -80.9459
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 40 s.
p2n2 variation ends
p0n3 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p0n3
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 0 pairs of protons, and 3 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=0 nNN=3 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -15.873766, it took 0 s elapsed time
iter = 2, f = -16.774428, it took 1 s elapsed time
iter = 3, f = -18.322549, it took 0 s elapsed time
iter = 4, f = -19.961001, it took 0 s elapsed time
iter = 5, f = -20.110183, it took 1 s elapsed time
iter = 6, f = -21.405903, it took 0 s elapsed time
iter = 7, f = -22.736544, it took 1 s elapsed time
iter = 8, f = -22.874802, it took 0 s elapsed time
iter = 9, f = -23.931006, it took 1 s elapsed time
iter = 10, f = -25.023011, it took 0 s elapsed time
iter = 11, f = -25.188678, it took 1 s elapsed time
iter = 12, f = -26.757793, it took 0 s elapsed time
iter = 13, f = -28.656497, it took 0 s elapsed time
iter = 14, f = -28.863399, it took 1 s elapsed time
iter = 15, f = -29.068410, it took 0 s elapsed time
iter = 16, f = -30.210214, it took 0 s elapsed time
iter = 17, f = -30.945952, it took 1 s elapsed time
iter = 18, f = -30.951043, it took 0 s elapsed time
iter = 19, f = -31.080232, it took 1 s elapsed time
iter = 20, f = -31.298284, it took 0 s elapsed time
iter = 21, f = -31.573591, it took 0 s elapsed time
iter = 22, f = -31.646914, it took 1 s elapsed time
iter = 23, f = -31.906232, it took 0 s elapsed time
iter = 24, f = -32.120034, it took 1 s elapsed time
iter = 25, f = -32.120352, it took 0 s elapsed time
iter = 26, f = -32.138048, it took 1 s elapsed time
iter = 27, f = -32.170933, it took 0 s elapsed time
iter = 28, f = -32.226744, it took 0 s elapsed time
iter = 29, f = -32.299353, it took 1 s elapsed time
iter = 30, f = -32.322190, it took 0 s elapsed time
iter = 31, f = -32.402822, it took 1 s elapsed time
iter = 32, f = -32.436764, it took 0 s elapsed time
iter = 33, f = -32.443692, it took 1 s elapsed time
iter = 34, f = -32.481549, it took 0 s elapsed time
iter = 35, f = -32.522339, it took 0 s elapsed time
iter = 36, f = -32.525964, it took 1 s elapsed time
iter = 37, f = -32.546715, it took 0 s elapsed time
iter = 38, f = -32.565608, it took 1 s elapsed time
iter = 39, f = -32.565866, it took 0 s elapsed time
iter = 40, f = -32.570863, it took 1 s elapsed time
iter = 41, f = -32.579799, it took 0 s elapsed time
n=132, Minimum=-32.593441 found after 42 iterations
************************************************
Found the Minimum!
************************************************
iter = 42, f = -32.593441, it took 0 s elapsed time
before normalization, pove = 1 nove = 3321.1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 17 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -32.5934
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 17 s.
p0n3 variation ends
p1n3 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p1n3
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 1 pairs of protons, and 3 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=3 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -42.828251, it took 1 s elapsed time
iter = 2, f = -43.876887, it took 0 s elapsed time
iter = 3, f = -45.876831, it took 1 s elapsed time
iter = 4, f = -48.823761, it took 0 s elapsed time
iter = 5, f = -50.196944, it took 1 s elapsed time
iter = 6, f = -50.210354, it took 1 s elapsed time
iter = 7, f = -50.564627, it took 1 s elapsed time
iter = 8, f = -51.235878, it took 0 s elapsed time
iter = 9, f = -52.412597, it took 1 s elapsed time
iter = 10, f = -54.068479, it took 0 s elapsed time
iter = 11, f = -55.199556, it took 1 s elapsed time
iter = 12, f = -55.208832, it took 1 s elapsed time
iter = 13, f = -55.676694, it took 1 s elapsed time
iter = 14, f = -56.496280, it took 0 s elapsed time
iter = 15, f = -57.665643, it took 1 s elapsed time
iter = 16, f = -58.501826, it took 1 s elapsed time
iter = 17, f = -58.510905, it took 1 s elapsed time
iter = 18, f = -58.950688, it took 0 s elapsed time
iter = 19, f = -59.736719, it took 1 s elapsed time
iter = 20, f = -60.905908, it took 0 s elapsed time
iter = 21, f = -61.709184, it took 1 s elapsed time
iter = 22, f = -61.709952, it took 1 s elapsed time
iter = 23, f = -61.784716, it took 1 s elapsed time
iter = 24, f = -61.931231, it took 0 s elapsed time
iter = 25, f = -62.211584, it took 1 s elapsed time
iter = 26, f = -62.717049, it took 1 s elapsed time
iter = 27, f = -63.480319, it took 0 s elapsed time
iter = 28, f = -63.991126, it took 1 s elapsed time
iter = 29, f = -63.994783, it took 1 s elapsed time
iter = 30, f = -64.103992, it took 1 s elapsed time
iter = 31, f = -64.304073, it took 0 s elapsed time
iter = 32, f = -64.631408, it took 1 s elapsed time
iter = 33, f = -65.009299, it took 0 s elapsed time
iter = 34, f = -65.078626, it took 1 s elapsed time
iter = 35, f = -65.524791, it took 1 s elapsed time
iter = 36, f = -65.937179, it took 1 s elapsed time
iter = 37, f = -65.953249, it took 1 s elapsed time
iter = 38, f = -66.148278, it took 0 s elapsed time
iter = 39, f = -66.367190, it took 1 s elapsed time
iter = 40, f = -66.398191, it took 1 s elapsed time
iter = 41, f = -66.559447, it took 1 s elapsed time
iter = 42, f = -66.679584, it took 0 s elapsed time
iter = 43, f = -66.679733, it took 2 s elapsed time
iter = 44, f = -66.689307, it took 0 s elapsed time
iter = 45, f = -66.707494, it took 1 s elapsed time
iter = 46, f = -66.740027, it took 0 s elapsed time
iter = 47, f = -66.789744, it took 1 s elapsed time
iter = 48, f = -66.828285, it took 0 s elapsed time
iter = 49, f = -66.828300, it took 2 s elapsed time
iter = 50, f = -66.830591, it took 0 s elapsed time
iter = 51, f = -66.835099, it took 1 s elapsed time
iter = 52, f = -66.843829, it took 0 s elapsed time
iter = 53, f = -66.860141, it took 1 s elapsed time
iter = 54, f = -66.888177, it took 1 s elapsed time
iter = 55, f = -66.925988, it took 0 s elapsed time
iter = 56, f = -66.929746, it took 1 s elapsed time
iter = 57, f = -66.939771, it took 1 s elapsed time
iter = 58, f = -66.984176, it took 1 s elapsed time
iter = 59, f = -67.012285, it took 0 s elapsed time
iter = 60, f = -67.013187, it took 1 s elapsed time
iter = 61, f = -67.025172, it took 1 s elapsed time
n=132, Minimum=-67.045046 found after 62 iterations
************************************************
Found the Minimum!
************************************************
iter = 62, f = -67.045046, it took 1 s elapsed time
before normalization, pove = 7.20872 nove = 4910.09
after normalization, pove = 1 nove = 1
---------- variation ends, it took 45 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -67.045
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 45 s.
p1n3 variation ends
p2n3 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p2n3
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 2 pairs of protons, and 3 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=3 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -67.604260, it took 1 s elapsed time
iter = 2, f = -68.839692, it took 0 s elapsed time
iter = 3, f = -71.091317, it took 1 s elapsed time
iter = 4, f = -74.112956, it took 1 s elapsed time
iter = 5, f = -75.000240, it took 0 s elapsed time
iter = 6, f = -75.167254, it took 1 s elapsed time
iter = 7, f = -76.471225, it took 1 s elapsed time
iter = 8, f = -78.354464, it took 1 s elapsed time
iter = 9, f = -79.609499, it took 0 s elapsed time
iter = 10, f = -79.610298, it took 2 s elapsed time
iter = 11, f = -79.733060, it took 0 s elapsed time
iter = 12, f = -79.971740, it took 1 s elapsed time
iter = 13, f = -80.421334, it took 0 s elapsed time
iter = 14, f = -81.207891, it took 1 s elapsed time
iter = 15, f = -82.338928, it took 1 s elapsed time
iter = 16, f = -83.120008, it took 0 s elapsed time
iter = 17, f = -83.120099, it took 1 s elapsed time
iter = 18, f = -83.156318, it took 1 s elapsed time
iter = 19, f = -83.228045, it took 1 s elapsed time
iter = 20, f = -83.368641, it took 0 s elapsed time
iter = 21, f = -83.638265, it took 1 s elapsed time
iter = 22, f = -84.130477, it took 0 s elapsed time
iter = 23, f = -84.925287, it took 1 s elapsed time
iter = 24, f = -85.805235, it took 1 s elapsed time
iter = 25, f = -85.956055, it took 1 s elapsed time
iter = 26, f = -87.025236, it took 1 s elapsed time
iter = 27, f = -88.013710, it took 0 s elapsed time
iter = 28, f = -88.049850, it took 1 s elapsed time
iter = 29, f = -88.563409, it took 1 s elapsed time
iter = 30, f = -89.352787, it took 1 s elapsed time
iter = 31, f = -89.965773, it took 0 s elapsed time
iter = 32, f = -89.966239, it took 1 s elapsed time
iter = 33, f = -90.031587, it took 1 s elapsed time
iter = 34, f = -90.160474, it took 1 s elapsed time
iter = 35, f = -90.410328, it took 0 s elapsed time
iter = 36, f = -90.873503, it took 1 s elapsed time
iter = 37, f = -91.627813, it took 0 s elapsed time
iter = 38, f = -92.426370, it took 1 s elapsed time
iter = 39, f = -92.537364, it took 1 s elapsed time
iter = 40, f = -93.580618, it took 1 s elapsed time
iter = 41, f = -94.394846, it took 1 s elapsed time
iter = 42, f = -94.409030, it took 1 s elapsed time
iter = 43, f = -94.773320, it took 0 s elapsed time
iter = 44, f = -95.236440, it took 1 s elapsed time
iter = 45, f = -95.284277, it took 1 s elapsed time
iter = 46, f = -95.385335, it took 1 s elapsed time
iter = 47, f = -95.818091, it took 0 s elapsed time
iter = 48, f = -96.177211, it took 1 s elapsed time
iter = 49, f = -96.180731, it took 1 s elapsed time
iter = 50, f = -96.279317, it took 1 s elapsed time
iter = 51, f = -96.442027, it took 0 s elapsed time
iter = 52, f = -96.635718, it took 1 s elapsed time
iter = 53, f = -96.674326, it took 1 s elapsed time
iter = 54, f = -96.888716, it took 1 s elapsed time
iter = 55, f = -97.150889, it took 1 s elapsed time
iter = 56, f = -97.209346, it took 1 s elapsed time
iter = 57, f = -97.463944, it took 0 s elapsed time
iter = 58, f = -97.606236, it took 1 s elapsed time
iter = 59, f = -97.614843, it took 1 s elapsed time
iter = 60, f = -97.713010, it took 1 s elapsed time
iter = 61, f = -97.873306, it took 0 s elapsed time
iter = 62, f = -98.054891, it took 1 s elapsed time
iter = 63, f = -98.080807, it took 1 s elapsed time
iter = 64, f = -98.242280, it took 1 s elapsed time
iter = 65, f = -98.301319, it took 0 s elapsed time
iter = 66, f = -98.318593, it took 2 s elapsed time
iter = 67, f = -98.409298, it took 0 s elapsed time
iter = 68, f = -98.490489, it took 1 s elapsed time
iter = 69, f = -98.491685, it took 1 s elapsed time
iter = 70, f = -98.515615, it took 1 s elapsed time
iter = 71, f = -98.554757, it took 0 s elapsed time
iter = 72, f = -98.598783, it took 1 s elapsed time
iter = 73, f = -98.604063, it took 1 s elapsed time
iter = 74, f = -98.639864, it took 1 s elapsed time
n=132, Minimum=-98.679397 found after 75 iterations
************************************************
Found the Minimum!
************************************************
iter = 75, f = -98.679397, it took 0 s elapsed time
before normalization, pove = 1243.97 nove = 112641
after normalization, pove = 1 nove = 1
---------- variation ends, it took 57 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -98.6794
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 57 s.
p2n3 variation ends
p3n3 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p3n3
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 3 pairs of protons, and 3 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=3 nNN=3 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -88.734116, it took 1 s elapsed time
iter = 2, f = -90.152596, it took 0 s elapsed time
iter = 3, f = -92.870659, it took 1 s elapsed time
iter = 4, f = -97.175394, it took 1 s elapsed time
iter = 5, f = -99.697672, it took 0 s elapsed time
iter = 6, f = -99.737322, it took 2 s elapsed time
iter = 7, f = -100.509275, it took 0 s elapsed time
iter = 8, f = -101.933554, it took 1 s elapsed time
iter = 9, f = -104.214305, it took 0 s elapsed time
iter = 10, f = -106.488762, it took 1 s elapsed time
iter = 11, f = -106.752781, it took 1 s elapsed time
iter = 12, f = -109.781155, it took 1 s elapsed time
iter = 13, f = -113.541555, it took 1 s elapsed time
iter = 14, f = -114.027854, it took 0 s elapsed time
iter = 15, f = -114.501463, it took 2 s elapsed time
iter = 16, f = -117.170133, it took 0 s elapsed time
iter = 17, f = -119.589794, it took 1 s elapsed time
iter = 18, f = -119.698065, it took 1 s elapsed time
iter = 19, f = -121.187902, it took 1 s elapsed time
iter = 20, f = -122.863692, it took 0 s elapsed time
iter = 21, f = -123.130052, it took 2 s elapsed time
iter = 22, f = -124.569299, it took 0 s elapsed time
iter = 23, f = -125.626343, it took 1 s elapsed time
iter = 24, f = -125.626534, it took 1 s elapsed time
iter = 25, f = -125.668822, it took 1 s elapsed time
iter = 26, f = -125.750937, it took 0 s elapsed time
iter = 27, f = -125.905322, it took 1 s elapsed time
iter = 28, f = -126.174784, it took 1 s elapsed time
iter = 29, f = -126.558995, it took 0 s elapsed time
iter = 30, f = -126.753906, it took 1 s elapsed time
iter = 31, f = -126.771451, it took 1 s elapsed time
iter = 32, f = -126.947725, it took 1 s elapsed time
iter = 33, f = -127.208293, it took 0 s elapsed time
iter = 34, f = -127.373361, it took 1 s elapsed time
iter = 35, f = -127.377661, it took 2 s elapsed time
iter = 36, f = -127.450403, it took 0 s elapsed time
iter = 37, f = -127.577720, it took 1 s elapsed time
iter = 38, f = -127.760038, it took 0 s elapsed time
iter = 39, f = -127.843087, it took 1 s elapsed time
iter = 40, f = -127.856620, it took 1 s elapsed time
iter = 41, f = -127.971786, it took 1 s elapsed time
iter = 42, f = -128.140997, it took 1 s elapsed time
iter = 43, f = -128.232094, it took 0 s elapsed time
iter = 44, f = -128.240333, it took 1 s elapsed time
iter = 45, f = -128.329968, it took 1 s elapsed time
iter = 46, f = -128.479837, it took 1 s elapsed time
iter = 47, f = -128.661294, it took 0 s elapsed time
iter = 48, f = -128.694352, it took 2 s elapsed time
iter = 49, f = -128.882424, it took 0 s elapsed time
iter = 50, f = -129.072097, it took 1 s elapsed time
iter = 51, f = -129.079597, it took 1 s elapsed time
iter = 52, f = -129.167920, it took 1 s elapsed time
iter = 53, f = -129.310870, it took 0 s elapsed time
iter = 54, f = -129.460426, it took 1 s elapsed time
iter = 55, f = -129.469514, it took 1 s elapsed time
iter = 56, f = -129.562449, it took 1 s elapsed time
iter = 57, f = -129.671249, it took 1 s elapsed time
iter = 58, f = -129.688778, it took 1 s elapsed time
iter = 59, f = -129.778757, it took 0 s elapsed time
iter = 60, f = -129.839967, it took 1 s elapsed time
iter = 61, f = -129.840705, it took 1 s elapsed time
iter = 62, f = -129.857362, it took 1 s elapsed time
iter = 63, f = -129.884759, it took 1 s elapsed time
iter = 64, f = -129.916174, it took 0 s elapsed time
iter = 65, f = -129.920379, it took 2 s elapsed time
iter = 66, f = -129.945370, it took 0 s elapsed time
iter = 67, f = -129.969521, it took 1 s elapsed time
iter = 68, f = -129.970354, it took 1 s elapsed time
iter = 69, f = -129.979788, it took 1 s elapsed time
n=132, Minimum=-129.991966 found after 70 iterations
************************************************
Found the Minimum!
************************************************
iter = 70, f = -129.991966, it took 0 s elapsed time
before normalization, pove = 7755.63 nove = 9353.9
after normalization, pove = 1 nove = 1
---------- variation ends, it took 55 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -129.992
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 55 s.
p3n3 variation ends
p0n4 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p0n4
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 0 pairs of protons, and 4 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=0 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = 16.969763, it took 0 s elapsed time
iter = 2, f = 16.183126, it took 0 s elapsed time
iter = 3, f = 14.654661, it took 1 s elapsed time
iter = 4, f = 12.130539, it took 0 s elapsed time
iter = 5, f = 9.796232, it took 0 s elapsed time
iter = 6, f = 9.555080, it took 1 s elapsed time
iter = 7, f = 7.420973, it took 0 s elapsed time
iter = 8, f = 5.779309, it took 1 s elapsed time
iter = 9, f = 5.774869, it took 0 s elapsed time
iter = 10, f = 5.621348, it took 1 s elapsed time
iter = 11, f = 5.338565, it took 0 s elapsed time
iter = 12, f = 4.876649, it took 0 s elapsed time
iter = 13, f = 4.368604, it took 0 s elapsed time
iter = 14, f = 4.293867, it took 1 s elapsed time
iter = 15, f = 3.660630, it took 0 s elapsed time
iter = 16, f = 2.666962, it took 1 s elapsed time
iter = 17, f = 1.730224, it took 0 s elapsed time
iter = 18, f = 1.707148, it took 0 s elapsed time
iter = 19, f = 1.621911, it took 1 s elapsed time
iter = 20, f = 0.887799, it took 0 s elapsed time
iter = 21, f = 0.382089, it took 1 s elapsed time
iter = 22, f = 0.382054, it took 0 s elapsed time
iter = 23, f = 0.371588, it took 1 s elapsed time
iter = 24, f = 0.350943, it took 0 s elapsed time
iter = 25, f = 0.310796, it took 0 s elapsed time
iter = 26, f = 0.235049, it took 1 s elapsed time
iter = 27, f = 0.101483, it took 0 s elapsed time
iter = 28, f = -0.096589, it took 0 s elapsed time
iter = 29, f = -0.245078, it took 0 s elapsed time
iter = 30, f = -0.245088, it took 1 s elapsed time
iter = 31, f = -0.248936, it took 0 s elapsed time
iter = 32, f = -0.256508, it took 1 s elapsed time
iter = 33, f = -0.271155, it took 0 s elapsed time
iter = 34, f = -0.298468, it took 0 s elapsed time
iter = 35, f = -0.345245, it took 1 s elapsed time
n=132, Minimum=-0.408098 found after 36 iterations
************************************************
Found the Minimum!
************************************************
iter = 36, f = -0.408098, it took 0 s elapsed time
before normalization, pove = 1 nove = 704.682
after normalization, pove = 1 nove = 1
---------- variation ends, it took 14 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -0.408098
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 14 s.
p0n4 variation ends
p1n4 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p1n4
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 1 pairs of protons, and 4 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -14.225662, it took 0 s elapsed time
iter = 2, f = -15.053828, it took 1 s elapsed time
iter = 3, f = -16.651497, it took 0 s elapsed time
iter = 4, f = -19.364579, it took 1 s elapsed time
iter = 5, f = -22.449738, it took 0 s elapsed time
iter = 6, f = -23.658808, it took 1 s elapsed time
iter = 7, f = -23.663766, it took 1 s elapsed time
iter = 8, f = -23.948178, it took 1 s elapsed time
iter = 9, f = -24.507655, it took 0 s elapsed time
iter = 10, f = -25.582008, it took 1 s elapsed time
iter = 11, f = -27.499503, it took 0 s elapsed time
iter = 12, f = -30.152579, it took 1 s elapsed time
iter = 13, f = -31.364402, it took 0 s elapsed time
iter = 14, f = -31.407283, it took 2 s elapsed time
iter = 15, f = -32.017530, it took 0 s elapsed time
iter = 16, f = -33.026830, it took 1 s elapsed time
iter = 17, f = -34.221247, it took 0 s elapsed time
iter = 18, f = -34.300576, it took 1 s elapsed time
iter = 19, f = -34.515393, it took 1 s elapsed time
iter = 20, f = -35.368696, it took 0 s elapsed time
iter = 21, f = -35.859834, it took 1 s elapsed time
iter = 22, f = -35.871881, it took 1 s elapsed time
iter = 23, f = -36.013532, it took 1 s elapsed time
iter = 24, f = -36.240936, it took 0 s elapsed time
iter = 25, f = -36.478977, it took 1 s elapsed time
iter = 26, f = -36.498564, it took 1 s elapsed time
iter = 27, f = -36.660087, it took 1 s elapsed time
iter = 28, f = -36.905790, it took 0 s elapsed time
iter = 29, f = -37.098825, it took 1 s elapsed time
iter = 30, f = -37.098954, it took 1 s elapsed time
iter = 31, f = -37.114581, it took 0 s elapsed time
iter = 32, f = -37.144865, it took 1 s elapsed time
iter = 33, f = -37.201540, it took 0 s elapsed time
iter = 34, f = -37.299302, it took 1 s elapsed time
iter = 35, f = -37.433234, it took 1 s elapsed time
iter = 36, f = -37.472543, it took 0 s elapsed time
iter = 37, f = -37.490985, it took 1 s elapsed time
iter = 38, f = -37.616548, it took 1 s elapsed time
iter = 39, f = -37.789562, it took 0 s elapsed time
iter = 40, f = -37.865748, it took 1 s elapsed time
iter = 41, f = -37.875817, it took 1 s elapsed time
iter = 42, f = -37.970072, it took 1 s elapsed time
iter = 43, f = -38.095283, it took 0 s elapsed time
iter = 44, f = -38.123222, it took 1 s elapsed time
iter = 45, f = -38.144225, it took 1 s elapsed time
iter = 46, f = -38.253050, it took 0 s elapsed time
iter = 47, f = -38.344880, it took 1 s elapsed time
iter = 48, f = -38.345724, it took 1 s elapsed time
iter = 49, f = -38.368458, it took 1 s elapsed time
n=132, Minimum=-38.406368 found after 50 iterations
************************************************
Found the Minimum!
************************************************
iter = 50, f = -38.406368, it took 0 s elapsed time
before normalization, pove = 16.4109 nove = 653.177
after normalization, pove = 1 nove = 1
---------- variation ends, it took 34 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -38.4064
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 34 s.
p1n4 variation ends
p2n4 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p2n4
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 2 pairs of protons, and 4 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -43.175452, it took 1 s elapsed time
iter = 2, f = -44.215662, it took 0 s elapsed time
iter = 3, f = -46.243090, it took 1 s elapsed time
iter = 4, f = -49.833224, it took 0 s elapsed time
iter = 5, f = -54.318499, it took 1 s elapsed time
iter = 6, f = -55.761699, it took 1 s elapsed time
iter = 7, f = -55.878357, it took 1 s elapsed time
iter = 8, f = -57.131856, it took 0 s elapsed time
iter = 9, f = -59.135310, it took 1 s elapsed time
iter = 10, f = -61.088529, it took 0 s elapsed time
iter = 11, f = -61.222237, it took 2 s elapsed time
iter = 12, f = -62.519959, it took 0 s elapsed time
iter = 13, f = -64.315498, it took 1 s elapsed time
iter = 14, f = -65.051200, it took 0 s elapsed time
iter = 15, f = -65.134987, it took 2 s elapsed time
iter = 16, f = -65.982222, it took 0 s elapsed time
iter = 17, f = -67.310836, it took 1 s elapsed time
iter = 18, f = -68.558825, it took 0 s elapsed time
iter = 19, f = -68.648688, it took 1 s elapsed time
iter = 20, f = -69.580765, it took 1 s elapsed time
iter = 21, f = -70.241774, it took 1 s elapsed time
iter = 22, f = -70.241782, it took 1 s elapsed time
iter = 23, f = -70.248234, it took 0 s elapsed time
iter = 24, f = -70.261074, it took 1 s elapsed time
iter = 25, f = -70.286498, it took 1 s elapsed time
iter = 26, f = -70.336323, it took 0 s elapsed time
iter = 27, f = -70.431884, it took 1 s elapsed time
iter = 28, f = -70.606744, it took 0 s elapsed time
iter = 29, f = -70.892674, it took 1 s elapsed time
iter = 30, f = -71.226330, it took 1 s elapsed time
iter = 31, f = -71.294035, it took 1 s elapsed time
iter = 32, f = -71.694250, it took 0 s elapsed time
iter = 33, f = -72.208602, it took 1 s elapsed time
iter = 34, f = -72.315804, it took 0 s elapsed time
iter = 35, f = -72.390625, it took 2 s elapsed time
iter = 36, f = -72.795047, it took 0 s elapsed time
iter = 37, f = -73.093605, it took 1 s elapsed time
iter = 38, f = -73.093656, it took 1 s elapsed time
iter = 39, f = -73.105981, it took 1 s elapsed time
iter = 40, f = -73.130281, it took 0 s elapsed time
iter = 41, f = -73.177493, it took 1 s elapsed time
iter = 42, f = -73.266387, it took 0 s elapsed time
iter = 43, f = -73.422298, it took 1 s elapsed time
iter = 44, f = -73.649454, it took 0 s elapsed time
iter = 45, f = -73.798592, it took 1 s elapsed time
iter = 46, f = -73.800042, it took 1 s elapsed time
iter = 47, f = -73.855935, it took 1 s elapsed time
iter = 48, f = -73.961621, it took 0 s elapsed time
iter = 49, f = -74.148680, it took 1 s elapsed time
iter = 50, f = -74.427645, it took 0 s elapsed time
iter = 51, f = -74.640300, it took 1 s elapsed time
iter = 52, f = -74.640429, it took 1 s elapsed time
iter = 53, f = -74.663881, it took 1 s elapsed time
iter = 54, f = -74.709713, it took 0 s elapsed time
iter = 55, f = -74.797102, it took 1 s elapsed time
iter = 56, f = -74.954940, it took 1 s elapsed time
iter = 57, f = -75.204739, it took 0 s elapsed time
iter = 58, f = -75.463681, it took 1 s elapsed time
iter = 59, f = -75.489271, it took 1 s elapsed time
iter = 60, f = -75.767127, it took 0 s elapsed time
iter = 61, f = -76.097426, it took 1 s elapsed time
iter = 62, f = -76.176316, it took 1 s elapsed time
iter = 63, f = -76.526114, it took 1 s elapsed time
iter = 64, f = -76.648969, it took 0 s elapsed time
iter = 65, f = -76.679205, it took 2 s elapsed time
n=132, Minimum=-76.853235 found after 66 iterations
************************************************
Found the Minimum!
************************************************
iter = 66, f = -76.853235, it took 0 s elapsed time
before normalization, pove = 7621.64 nove = 10601.5
after normalization, pove = 1 nove = 1
---------- variation ends, it took 47 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -76.8532
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 47 s.
p2n4 variation ends
p3n4 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p3n4
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 3 pairs of protons, and 4 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=3 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -75.024556, it took 0 s elapsed time
iter = 2, f = -76.199234, it took 1 s elapsed time
iter = 3, f = -78.404882, it took 1 s elapsed time
iter = 4, f = -81.958265, it took 0 s elapsed time
iter = 5, f = -85.498737, it took 1 s elapsed time
iter = 6, f = -85.952676, it took 0 s elapsed time
iter = 7, f = -86.106214, it took 2 s elapsed time
iter = 8, f = -88.086739, it took 0 s elapsed time
iter = 9, f = -91.277084, it took 1 s elapsed time
iter = 10, f = -94.193535, it took 0 s elapsed time
iter = 11, f = -94.389166, it took 2 s elapsed time
iter = 12, f = -97.174539, it took 0 s elapsed time
iter = 13, f = -100.558630, it took 1 s elapsed time
iter = 14, f = -101.201430, it took 1 s elapsed time
iter = 15, f = -103.375125, it took 1 s elapsed time
iter = 16, f = -103.674659, it took 1 s elapsed time
iter = 17, f = -104.784063, it took 1 s elapsed time
iter = 18, f = -104.953468, it took 0 s elapsed time
iter = 19, f = -105.160748, it took 1 s elapsed time
iter = 20, f = -105.933441, it took 1 s elapsed time
iter = 21, f = -106.688183, it took 1 s elapsed time
iter = 22, f = -106.733478, it took 1 s elapsed time
iter = 23, f = -107.196218, it took 0 s elapsed time
iter = 24, f = -107.773647, it took 1 s elapsed time
iter = 25, f = -107.818948, it took 1 s elapsed time
iter = 26, f = -107.945992, it took 1 s elapsed time
iter = 27, f = -108.397890, it took 0 s elapsed time
iter = 28, f = -108.494693, it took 2 s elapsed time
iter = 29, f = -108.700518, it took 0 s elapsed time
iter = 30, f = -108.714108, it took 2 s elapsed time
iter = 31, f = -108.793740, it took 0 s elapsed time
iter = 32, f = -108.903039, it took 1 s elapsed time
iter = 33, f = -108.922937, it took 0 s elapsed time
iter = 34, f = -108.946498, it took 2 s elapsed time
iter = 35, f = -109.046158, it took 0 s elapsed time
iter = 36, f = -109.139654, it took 1 s elapsed time
iter = 37, f = -109.141570, it took 1 s elapsed time
iter = 38, f = -109.173662, it took 1 s elapsed time
iter = 39, f = -109.228339, it took 0 s elapsed time
iter = 40, f = -109.298887, it took 1 s elapsed time
iter = 41, f = -109.317404, it took 1 s elapsed time
iter = 42, f = -109.394951, it took 1 s elapsed time
iter = 43, f = -109.441876, it took 0 s elapsed time
iter = 44, f = -109.444137, it took 2 s elapsed time
iter = 45, f = -109.471359, it took 0 s elapsed time
iter = 46, f = -109.513531, it took 1 s elapsed time
iter = 47, f = -109.548916, it took 0 s elapsed time
iter = 48, f = -109.548981, it took 2 s elapsed time
iter = 49, f = -109.553456, it took 0 s elapsed time
iter = 50, f = -109.561972, it took 1 s elapsed time
iter = 51, f = -109.577282, it took 1 s elapsed time
iter = 52, f = -109.601094, it took 0 s elapsed time
iter = 53, f = -109.622120, it took 1 s elapsed time
iter = 54, f = -109.622316, it took 1 s elapsed time
iter = 55, f = -109.629200, it took 1 s elapsed time
iter = 56, f = -109.642079, it took 0 s elapsed time
iter = 57, f = -109.664298, it took 1 s elapsed time
iter = 58, f = -109.694743, it took 0 s elapsed time
iter = 59, f = -109.700772, it took 1 s elapsed time
iter = 60, f = -109.707046, it took 1 s elapsed time
iter = 61, f = -109.742494, it took 1 s elapsed time
iter = 62, f = -109.778525, it took 0 s elapsed time
iter = 63, f = -109.780226, it took 2 s elapsed time
iter = 64, f = -109.798950, it took 0 s elapsed time
n=132, Minimum=-109.826185 found after 65 iterations
************************************************
Found the Minimum!
************************************************
iter = 65, f = -109.826185, it took 1 s elapsed time
before normalization, pove = 28127.8 nove = 400.542
after normalization, pove = 1 nove = 1
---------- variation ends, it took 51 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -109.826
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 51 s.
p3n4 variation ends
p4n4 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p4n4
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 4 pairs of protons, and 4 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=2 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = 151.288626, it took 1 s elapsed time
iter = 2, f = 150.168286, it took 1 s elapsed time
iter = 3, f = 147.990449, it took 0 s elapsed time
iter = 4, f = 144.227536, it took 1 s elapsed time
iter = 5, f = 140.342335, it took 0 s elapsed time
iter = 6, f = 140.024340, it took 2 s elapsed time
iter = 7, f = 137.304795, it took 0 s elapsed time
iter = 8, f = 133.684258, it took 1 s elapsed time
iter = 9, f = 132.484583, it took 0 s elapsed time
iter = 10, f = 132.309290, it took 1 s elapsed time
iter = 11, f = 130.761391, it took 1 s elapsed time
iter = 12, f = 128.411278, it took 1 s elapsed time
iter = 13, f = 126.632823, it took 0 s elapsed time
iter = 14, f = 126.618582, it took 1 s elapsed time
iter = 15, f = 126.036683, it took 1 s elapsed time
iter = 16, f = 125.023972, it took 1 s elapsed time
iter = 17, f = 123.591001, it took 0 s elapsed time
iter = 18, f = 122.585523, it took 1 s elapsed time
iter = 19, f = 122.578879, it took 1 s elapsed time
iter = 20, f = 122.259435, it took 0 s elapsed time
iter = 21, f = 121.733623, it took 1 s elapsed time
iter = 22, f = 121.111247, it took 1 s elapsed time
iter = 23, f = 120.975760, it took 1 s elapsed time
iter = 24, f = 120.599180, it took 0 s elapsed time
iter = 25, f = 120.542926, it took 1 s elapsed time
iter = 26, f = 120.386573, it took 1 s elapsed time
iter = 27, f = 120.330296, it took 1 s elapsed time
iter = 28, f = 120.312895, it took 1 s elapsed time
iter = 29, f = 120.244397, it took 0 s elapsed time
iter = 30, f = 120.194088, it took 1 s elapsed time
iter = 31, f = 120.193969, it took 1 s elapsed time
iter = 32, f = 120.188755, it took 1 s elapsed time
iter = 33, f = 120.178792, it took 0 s elapsed time
iter = 34, f = 120.160717, it took 1 s elapsed time
n=132, Minimum=120.131920 found after 35 iterations
************************************************
Found the Minimum!
************************************************
iter = 35, f = 120.131920, it took 1 s elapsed time
before normalization, pove = 112.061 nove = 113.132
after normalization, pove = 1 nove = 1
---------- variation ends, it took 26 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 120.132
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 26 s.
p4n4 variation ends
p0n5 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p0n5
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 0 pairs of protons, and 5 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=0 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = 8.331228, it took 0 s elapsed time
iter = 2, f = 7.971652, it took 1 s elapsed time
iter = 3, f = 7.266919, it took 0 s elapsed time
iter = 4, f = 5.997284, it took 0 s elapsed time
iter = 5, f = 4.356524, it took 0 s elapsed time
iter = 6, f = 3.854712, it took 1 s elapsed time
iter = 7, f = 3.798076, it took 0 s elapsed time
iter = 8, f = 3.330537, it took 1 s elapsed time
iter = 9, f = 2.546446, it took 0 s elapsed time
iter = 10, f = 1.569516, it took 0 s elapsed time
iter = 11, f = 1.020201, it took 1 s elapsed time
iter = 12, f = 1.011274, it took 0 s elapsed time
iter = 13, f = 0.775533, it took 0 s elapsed time
n=132, Minimum=0.598326 found after 14 iterations
************************************************
Found the Minimum!
************************************************
iter = 14, f = 0.598326, it took 1 s elapsed time
before normalization, pove = 1 nove = 9.82308
after normalization, pove = 1 nove = 1
---------- variation ends, it took 5 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 0.598326
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 5 s.
p0n5 variation ends
p1n5 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p1n5
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 1 pairs of protons, and 5 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -27.302706, it took 1 s elapsed time
iter = 2, f = -27.837065, it took 0 s elapsed time
iter = 3, f = -28.892633, it took 1 s elapsed time
iter = 4, f = -30.884905, it took 1 s elapsed time
iter = 5, f = -34.014574, it took 0 s elapsed time
iter = 6, f = -36.443784, it took 1 s elapsed time
iter = 7, f = -36.497027, it took 1 s elapsed time
iter = 8, f = -37.459107, it took 0 s elapsed time
iter = 9, f = -38.956444, it took 1 s elapsed time
iter = 10, f = -40.436437, it took 0 s elapsed time
iter = 11, f = -40.446441, it took 1 s elapsed time
iter = 12, f = -40.694431, it took 2 s elapsed time
iter = 13, f = -42.102641, it took 0 s elapsed time
iter = 14, f = -42.306855, it took 1 s elapsed time
iter = 15, f = -42.485520, it took 1 s elapsed time
iter = 16, f = -43.491494, it took 1 s elapsed time
iter = 17, f = -44.181012, it took 0 s elapsed time
iter = 18, f = -44.182896, it took 1 s elapsed time
iter = 19, f = -44.295887, it took 1 s elapsed time
iter = 20, f = -44.501106, it took 1 s elapsed time
iter = 21, f = -44.832890, it took 0 s elapsed time
iter = 22, f = -45.228913, it took 1 s elapsed time
iter = 23, f = -45.341265, it took 0 s elapsed time
iter = 24, f = -45.367509, it took 1 s elapsed time
iter = 25, f = -45.574469, it took 1 s elapsed time
n=132, Minimum=-45.748480 found after 26 iterations
************************************************
Found the Minimum!
************************************************
iter = 26, f = -45.748480, it took 0 s elapsed time
before normalization, pove = 28.41 nove = 30.7967
after normalization, pove = 1 nove = 1
---------- variation ends, it took 18 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -45.7485
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 18 s.
p1n5 variation ends
p2n5 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p2n5
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 2 pairs of protons, and 5 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -64.293587, it took 0 s elapsed time
iter = 2, f = -65.047589, it took 1 s elapsed time
iter = 3, f = -66.480283, it took 1 s elapsed time
iter = 4, f = -68.830234, it took 0 s elapsed time
iter = 5, f = -71.281425, it took 1 s elapsed time
iter = 6, f = -71.960786, it took 0 s elapsed time
iter = 7, f = -71.997234, it took 1 s elapsed time
iter = 8, f = -72.789434, it took 1 s elapsed time
iter = 9, f = -74.215462, it took 1 s elapsed time
iter = 10, f = -76.425483, it took 0 s elapsed time
iter = 11, f = -78.584752, it took 1 s elapsed time
iter = 12, f = -78.854569, it took 1 s elapsed time
iter = 13, f = -81.043020, it took 0 s elapsed time
iter = 14, f = -82.319607, it took 1 s elapsed time
iter = 15, f = -82.332913, it took 1 s elapsed time
iter = 16, f = -82.560104, it took 1 s elapsed time
iter = 17, f = -82.949241, it took 0 s elapsed time
iter = 18, f = -83.460812, it took 1 s elapsed time
iter = 19, f = -83.517629, it took 0 s elapsed time
iter = 20, f = -83.632504, it took 1 s elapsed time
iter = 21, f = -84.057282, it took 1 s elapsed time
iter = 22, f = -84.469701, it took 1 s elapsed time
iter = 23, f = -84.487211, it took 1 s elapsed time
iter = 24, f = -84.695180, it took 0 s elapsed time
iter = 25, f = -85.014532, it took 1 s elapsed time
iter = 26, f = -85.286922, it took 0 s elapsed time
iter = 27, f = -85.290519, it took 2 s elapsed time
iter = 28, f = -85.399686, it took 0 s elapsed time
iter = 29, f = -85.595524, it took 1 s elapsed time
iter = 30, f = -85.896312, it took 0 s elapsed time
iter = 31, f = -86.160892, it took 1 s elapsed time
iter = 32, f = -86.167454, it took 1 s elapsed time
iter = 33, f = -86.322141, it took 0 s elapsed time
iter = 34, f = -86.580630, it took 1 s elapsed time
iter = 35, f = -86.911976, it took 1 s elapsed time
iter = 36, f = -87.014936, it took 0 s elapsed time
iter = 37, f = -87.044145, it took 1 s elapsed time
iter = 38, f = -87.318169, it took 1 s elapsed time
iter = 39, f = -87.654328, it took 0 s elapsed time
iter = 40, f = -87.684956, it took 1 s elapsed time
iter = 41, f = -87.754753, it took 1 s elapsed time
iter = 42, f = -87.997874, it took 1 s elapsed time
iter = 43, f = -88.020212, it took 1 s elapsed time
n=132, Minimum=-88.132405 found after 44 iterations
************************************************
Found the Minimum!
************************************************
iter = 44, f = -88.132405, it took 0 s elapsed time
before normalization, pove = 2113.94 nove = 39.7327
after normalization, pove = 1 nove = 1
---------- variation ends, it took 31 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -88.1324
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 31 s.
p2n5 variation ends
p3n5 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p3n5
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 3 pairs of protons, and 5 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=3 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -99.861031, it took 1 s elapsed time
iter = 2, f = -101.132105, it took 0 s elapsed time
iter = 3, f = -103.402581, it took 1 s elapsed time
iter = 4, f = -106.058260, it took 0 s elapsed time
iter = 5, f = -106.099052, it took 1 s elapsed time
iter = 6, f = -106.349683, it took 1 s elapsed time
iter = 7, f = -108.097929, it took 1 s elapsed time
iter = 8, f = -110.808441, it took 0 s elapsed time
iter = 9, f = -113.369565, it took 1 s elapsed time
iter = 10, f = -113.662911, it took 1 s elapsed time
iter = 11, f = -113.804085, it took 1 s elapsed time
iter = 12, f = -115.388948, it took 0 s elapsed time
iter = 13, f = -117.553864, it took 1 s elapsed time
iter = 14, f = -117.784686, it took 0 s elapsed time
iter = 15, f = -118.200352, it took 2 s elapsed time
iter = 16, f = -119.507322, it took 0 s elapsed time
iter = 17, f = -120.723489, it took 1 s elapsed time
iter = 18, f = -120.761645, it took 1 s elapsed time
iter = 19, f = -121.377169, it took 1 s elapsed time
iter = 20, f = -122.425638, it took 0 s elapsed time
iter = 21, f = -123.646673, it took 1 s elapsed time
iter = 22, f = -123.858364, it took 1 s elapsed time
iter = 23, f = -125.066514, it took 1 s elapsed time
iter = 24, f = -125.679912, it took 0 s elapsed time
iter = 25, f = -125.699721, it took 1 s elapsed time
iter = 26, f = -125.981218, it took 1 s elapsed time
iter = 27, f = -126.410832, it took 1 s elapsed time
iter = 28, f = -126.779209, it took 0 s elapsed time
iter = 29, f = -126.784490, it took 1 s elapsed time
iter = 30, f = -126.920301, it took 1 s elapsed time
iter = 31, f = -127.153843, it took 1 s elapsed time
iter = 32, f = -127.476132, it took 0 s elapsed time
iter = 33, f = -127.619442, it took 1 s elapsed time
iter = 34, f = -127.638463, it took 1 s elapsed time
iter = 35, f = -127.809240, it took 0 s elapsed time
iter = 36, f = -128.005534, it took 1 s elapsed time
iter = 37, f = -128.034922, it took 1 s elapsed time
iter = 38, f = -128.176465, it took 1 s elapsed time
iter = 39, f = -128.301052, it took 0 s elapsed time
iter = 40, f = -128.302358, it took 2 s elapsed time
iter = 41, f = -128.330244, it took 0 s elapsed time
iter = 42, f = -128.375773, it took 1 s elapsed time
iter = 43, f = -128.426123, it took 0 s elapsed time
iter = 44, f = -128.431298, it took 2 s elapsed time
iter = 45, f = -128.466503, it took 0 s elapsed time
iter = 46, f = -128.512552, it took 1 s elapsed time
iter = 47, f = -128.526361, it took 1 s elapsed time
n=132, Minimum=-128.567605 found after 48 iterations
************************************************
Found the Minimum!
************************************************
iter = 48, f = -128.567605, it took 0 s elapsed time
before normalization, pove = 22536.6 nove = 9.82908
after normalization, pove = 1 nove = 1
---------- variation ends, it took 36 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -128.568
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 36 s.
p3n5 variation ends
p4n5 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p4n5
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 4 pairs of protons, and 5 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = 118.767192, it took 1 s elapsed time
iter = 2, f = 117.818844, it took 0 s elapsed time
iter = 3, f = 115.920620, it took 1 s elapsed time
iter = 4, f = 112.490433, it took 0 s elapsed time
iter = 5, f = 108.356993, it took 1 s elapsed time
iter = 6, f = 106.952725, it took 0 s elapsed time
iter = 7, f = 106.907295, it took 2 s elapsed time
iter = 8, f = 106.112748, it took 0 s elapsed time
iter = 9, f = 104.674163, it took 1 s elapsed time
iter = 10, f = 102.461366, it took 0 s elapsed time
iter = 11, f = 100.499035, it took 1 s elapsed time
iter = 12, f = 100.385148, it took 1 s elapsed time
iter = 13, f = 99.206273, it took 1 s elapsed time
iter = 14, f = 98.065888, it took 0 s elapsed time
iter = 15, f = 98.000197, it took 1 s elapsed time
iter = 16, f = 97.385033, it took 1 s elapsed time
iter = 17, f = 96.539226, it took 0 s elapsed time
iter = 18, f = 96.224826, it took 1 s elapsed time
iter = 19, f = 96.167064, it took 1 s elapsed time
iter = 20, f = 95.702443, it took 1 s elapsed time
iter = 21, f = 95.084938, it took 0 s elapsed time
iter = 22, f = 94.829910, it took 1 s elapsed time
iter = 23, f = 94.802578, it took 1 s elapsed time
iter = 24, f = 94.530533, it took 0 s elapsed time
iter = 25, f = 94.125209, it took 1 s elapsed time
iter = 26, f = 93.789111, it took 1 s elapsed time
iter = 27, f = 93.783178, it took 1 s elapsed time
iter = 28, f = 93.637306, it took 0 s elapsed time
iter = 29, f = 93.435793, it took 1 s elapsed time
iter = 30, f = 93.357001, it took 0 s elapsed time
iter = 31, f = 93.339514, it took 2 s elapsed time
iter = 32, f = 93.245783, it took 0 s elapsed time
iter = 33, f = 93.151794, it took 1 s elapsed time
iter = 34, f = 93.146391, it took 1 s elapsed time
iter = 35, f = 93.101012, it took 0 s elapsed time
n=132, Minimum=93.043283 found after 36 iterations
************************************************
Found the Minimum!
************************************************
iter = 36, f = 93.043283, it took 1 s elapsed time
before normalization, pove = 469.831 nove = 10.6003
after normalization, pove = 1 nove = 1
---------- variation ends, it took 26 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 93.0433
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 26 s.
p4n5 variation ends
p5n5 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p5n5
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 5 pairs of protons, and 5 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=1 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = 78.467562, it took 1 s elapsed time
iter = 2, f = 77.852522, it took 1 s elapsed time
iter = 3, f = 76.638248, it took 0 s elapsed time
iter = 4, f = 74.356947, it took 1 s elapsed time
iter = 5, f = 70.862018, it took 0 s elapsed time
iter = 6, f = 68.540372, it took 1 s elapsed time
iter = 7, f = 68.539094, it took 1 s elapsed time
iter = 8, f = 68.401074, it took 0 s elapsed time
iter = 9, f = 68.132131, it took 1 s elapsed time
iter = 10, f = 67.623761, it took 1 s elapsed time
iter = 11, f = 66.730735, it took 0 s elapsed time
iter = 12, f = 65.439153, it took 1 s elapsed time
iter = 13, f = 64.425562, it took 0 s elapsed time
iter = 14, f = 64.389017, it took 1 s elapsed time
iter = 15, f = 63.701320, it took 1 s elapsed time
iter = 16, f = 62.757664, it took 0 s elapsed time
iter = 17, f = 62.176362, it took 1 s elapsed time
iter = 18, f = 62.176340, it took 1 s elapsed time
iter = 19, f = 62.168173, it took 0 s elapsed time
iter = 20, f = 62.151944, it took 1 s elapsed time
iter = 21, f = 62.119905, it took 1 s elapsed time
iter = 22, f = 62.057506, it took 0 s elapsed time
iter = 23, f = 61.939401, it took 1 s elapsed time
iter = 24, f = 61.729674, it took 0 s elapsed time
iter = 25, f = 61.411941, it took 1 s elapsed time
iter = 26, f = 61.129456, it took 0 s elapsed time
iter = 27, f = 61.118942, it took 1 s elapsed time
iter = 28, f = 60.945253, it took 1 s elapsed time
iter = 29, f = 60.695825, it took 0 s elapsed time
iter = 30, f = 60.532738, it took 1 s elapsed time
iter = 31, f = 60.531948, it took 1 s elapsed time
iter = 32, f = 60.507401, it took 1 s elapsed time
iter = 33, f = 60.464700, it took 0 s elapsed time
n=132, Minimum=60.404223 found after 34 iterations
************************************************
Found the Minimum!
************************************************
iter = 34, f = 60.404223, it took 1 s elapsed time
before normalization, pove = 12.159 nove = 11.0875
after normalization, pove = 1 nove = 1
---------- variation ends, it took 22 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 60.4042
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 22 s.
p5n5 variation ends
p0n6 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p0n6
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 0 pairs of protons, and 6 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=0 nNN=0 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
before normalization, pove = 1 nove = 1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 0 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 0
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 0 s.
p0n6 variation ends
p1n6 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p1n6
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 1 pairs of protons, and 6 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=0 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -41.564958, it took 0 s elapsed time
iter = 2, f = -41.943032, it took 1 s elapsed time
iter = 3, f = -42.671068, it took 0 s elapsed time
iter = 4, f = -43.934172, it took 0 s elapsed time
iter = 5, f = -45.410319, it took 1 s elapsed time
iter = 6, f = -45.480844, it took 0 s elapsed time
iter = 7, f = -45.729090, it took 1 s elapsed time
iter = 8, f = -46.860100, it took 0 s elapsed time
iter = 9, f = -47.286964, it took 0 s elapsed time
iter = 10, f = -47.313737, it took 1 s elapsed time
iter = 11, f = -47.489899, it took 0 s elapsed time
iter = 12, f = -47.728294, it took 1 s elapsed time
iter = 13, f = -47.839323, it took 0 s elapsed time
iter = 14, f = -47.847463, it took 1 s elapsed time
iter = 15, f = -47.936451, it took 0 s elapsed time
iter = 16, f = -48.102584, it took 0 s elapsed time
iter = 17, f = -48.381421, it took 0 s elapsed time
iter = 18, f = -48.734714, it took 1 s elapsed time
iter = 19, f = -48.953123, it took 0 s elapsed time
iter = 20, f = -48.961763, it took 1 s elapsed time
iter = 21, f = -49.297883, it took 0 s elapsed time
n=132, Minimum=-49.480323 found after 22 iterations
************************************************
Found the Minimum!
************************************************
iter = 22, f = -49.480323, it took 0 s elapsed time
before normalization, pove = 34.5253 nove = 1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 8 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -49.4803
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 8 s.
p1n6 variation ends
p2n6 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p2n6
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 2 pairs of protons, and 6 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=0 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -83.360996, it took 0 s elapsed time
iter = 2, f = -84.113410, it took 0 s elapsed time
iter = 3, f = -85.537747, it took 1 s elapsed time
iter = 4, f = -87.723161, it took 0 s elapsed time
iter = 5, f = -89.147402, it took 0 s elapsed time
iter = 6, f = -89.147884, it took 1 s elapsed time
iter = 7, f = -89.220622, it took 0 s elapsed time
iter = 8, f = -89.364167, it took 0 s elapsed time
iter = 9, f = -89.643201, it took 1 s elapsed time
iter = 10, f = -90.166691, it took 0 s elapsed time
iter = 11, f = -91.061305, it took 0 s elapsed time
iter = 12, f = -92.206474, it took 1 s elapsed time
iter = 13, f = -92.488538, it took 0 s elapsed time
iter = 14, f = -92.599049, it took 1 s elapsed time
iter = 15, f = -93.395718, it took 0 s elapsed time
iter = 16, f = -94.450758, it took 0 s elapsed time
iter = 17, f = -94.687425, it took 1 s elapsed time
iter = 18, f = -94.814501, it took 0 s elapsed time
iter = 19, f = -95.274381, it took 1 s elapsed time
iter = 20, f = -95.387426, it took 0 s elapsed time
iter = 21, f = -95.576954, it took 1 s elapsed time
iter = 22, f = -95.599488, it took 0 s elapsed time
iter = 23, f = -95.661585, it took 1 s elapsed time
iter = 24, f = -95.677063, it took 0 s elapsed time
iter = 25, f = -95.718866, it took 1 s elapsed time
iter = 26, f = -95.760444, it took 0 s elapsed time
iter = 27, f = -95.762161, it took 0 s elapsed time
iter = 28, f = -95.782198, it took 1 s elapsed time
iter = 29, f = -95.817927, it took 0 s elapsed time
iter = 30, f = -95.871331, it took 0 s elapsed time
iter = 31, f = -95.904984, it took 1 s elapsed time
iter = 32, f = -95.905974, it took 0 s elapsed time
iter = 33, f = -95.926655, it took 1 s elapsed time
iter = 34, f = -95.961525, it took 0 s elapsed time
iter = 35, f = -96.006171, it took 0 s elapsed time
iter = 36, f = -96.007413, it took 1 s elapsed time
iter = 37, f = -96.019957, it took 0 s elapsed time
iter = 38, f = -96.074268, it took 1 s elapsed time
iter = 39, f = -96.105180, it took 0 s elapsed time
iter = 40, f = -96.106366, it took 0 s elapsed time
iter = 41, f = -96.123416, it took 1 s elapsed time
iter = 42, f = -96.149251, it took 0 s elapsed time
iter = 43, f = -96.170256, it took 0 s elapsed time
iter = 44, f = -96.170309, it took 1 s elapsed time
iter = 45, f = -96.173751, it took 0 s elapsed time
iter = 46, f = -96.180304, it took 1 s elapsed time
n=132, Minimum=-96.192100 found after 47 iterations
************************************************
Found the Minimum!
************************************************
iter = 47, f = -96.192100, it took 0 s elapsed time
before normalization, pove = 109.405 nove = 1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 19 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -96.1921
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 19 s.
p2n6 variation ends
p3n6 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p3n6
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 3 pairs of protons, and 6 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=3 nNN=0 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = -123.141246, it took 1 s elapsed time
iter = 2, f = -124.155035, it took 0 s elapsed time
iter = 3, f = -126.037987, it took 0 s elapsed time
iter = 4, f = -128.628477, it took 1 s elapsed time
iter = 5, f = -129.148295, it took 0 s elapsed time
iter = 6, f = -129.390077, it took 1 s elapsed time
iter = 7, f = -130.774132, it took 0 s elapsed time
iter = 8, f = -132.848007, it took 0 s elapsed time
iter = 9, f = -134.410430, it took 1 s elapsed time
iter = 10, f = -134.465066, it took 0 s elapsed time
iter = 11, f = -135.582313, it took 1 s elapsed time
iter = 12, f = -136.347941, it took 0 s elapsed time
iter = 13, f = -136.352818, it took 1 s elapsed time
iter = 14, f = -136.455955, it took 0 s elapsed time
iter = 15, f = -136.646355, it took 0 s elapsed time
iter = 16, f = -136.958086, it took 1 s elapsed time
iter = 17, f = -137.283150, it took 0 s elapsed time
iter = 18, f = -137.304414, it took 1 s elapsed time
iter = 19, f = -137.553195, it took 0 s elapsed time
iter = 20, f = -138.017916, it took 0 s elapsed time
iter = 21, f = -138.753221, it took 1 s elapsed time
iter = 22, f = -139.376289, it took 0 s elapsed time
iter = 23, f = -139.394288, it took 1 s elapsed time
iter = 24, f = -139.845905, it took 0 s elapsed time
iter = 25, f = -140.510846, it took 0 s elapsed time
iter = 26, f = -141.027453, it took 1 s elapsed time
iter = 27, f = -141.034942, it took 0 s elapsed time
iter = 28, f = -141.309477, it took 1 s elapsed time
iter = 29, f = -141.744437, it took 0 s elapsed time
iter = 30, f = -142.207993, it took 0 s elapsed time
iter = 31, f = -142.275901, it took 1 s elapsed time
iter = 32, f = -142.742544, it took 1 s elapsed time
iter = 33, f = -142.952921, it took 0 s elapsed time
iter = 34, f = -142.976021, it took 1 s elapsed time
iter = 35, f = -143.153181, it took 0 s elapsed time
iter = 36, f = -143.364480, it took 0 s elapsed time
iter = 37, f = -143.403187, it took 1 s elapsed time
iter = 38, f = -143.525674, it took 0 s elapsed time
iter = 39, f = -143.558386, it took 1 s elapsed time
n=132, Minimum=-143.635512 found after 40 iterations
************************************************
Found the Minimum!
************************************************
iter = 40, f = -143.635512, it took 0 s elapsed time
before normalization, pove = 13481.3 nove = 1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 17 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = -143.636
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 17 s.
p3n6 variation ends
p4n6 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p4n6
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 4 pairs of protons, and 6 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=2 nNN=0 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = 86.150036, it took 0 s elapsed time
iter = 2, f = 85.275044, it took 0 s elapsed time
iter = 3, f = 83.581779, it took 1 s elapsed time
iter = 4, f = 80.915479, it took 0 s elapsed time
iter = 5, f = 78.916331, it took 0 s elapsed time
iter = 6, f = 78.852715, it took 1 s elapsed time
iter = 7, f = 77.652179, it took 0 s elapsed time
iter = 8, f = 75.836853, it took 0 s elapsed time
iter = 9, f = 74.296819, it took 1 s elapsed time
iter = 10, f = 74.229506, it took 0 s elapsed time
iter = 11, f = 73.157085, it took 1 s elapsed time
iter = 12, f = 71.874295, it took 0 s elapsed time
iter = 13, f = 71.837896, it took 0 s elapsed time
iter = 14, f = 71.580779, it took 1 s elapsed time
iter = 15, f = 70.484935, it took 1 s elapsed time
iter = 16, f = 70.285186, it took 0 s elapsed time
iter = 17, f = 69.719961, it took 1 s elapsed time
iter = 18, f = 69.637247, it took 0 s elapsed time
iter = 19, f = 69.361401, it took 1 s elapsed time
iter = 20, f = 69.289099, it took 0 s elapsed time
iter = 21, f = 69.251596, it took 1 s elapsed time
iter = 22, f = 69.114871, it took 0 s elapsed time
iter = 23, f = 69.053646, it took 0 s elapsed time
iter = 24, f = 69.044633, it took 1 s elapsed time
n=132, Minimum=69.006790 found after 25 iterations
************************************************
Found the Minimum!
************************************************
iter = 25, f = 69.006790, it took 0 s elapsed time
before normalization, pove = 151.492 nove = 1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 11 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 69.0068
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 11 s.
p4n6 variation ends
p5n6 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p5n6
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 5 pairs of protons, and 6 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=1 nNN=0 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
iter = 1, f = 42.850317, it took 1 s elapsed time
iter = 2, f = 42.447670, it took 0 s elapsed time
iter = 3, f = 41.665004, it took 0 s elapsed time
iter = 4, f = 40.279250, it took 1 s elapsed time
iter = 5, f = 38.566597, it took 0 s elapsed time
iter = 6, f = 38.233471, it took 0 s elapsed time
iter = 7, f = 38.110743, it took 1 s elapsed time
iter = 8, f = 37.379989, it took 0 s elapsed time
iter = 9, f = 36.285801, it took 1 s elapsed time
iter = 10, f = 35.286869, it took 0 s elapsed time
iter = 11, f = 35.166490, it took 0 s elapsed time
iter = 12, f = 35.093965, it took 1 s elapsed time
iter = 13, f = 34.564683, it took 0 s elapsed time
iter = 14, f = 34.472054, it took 1 s elapsed time
iter = 15, f = 34.189365, it took 0 s elapsed time
iter = 16, f = 33.888625, it took 1 s elapsed time
iter = 17, f = 33.836348, it took 0 s elapsed time
n=132, Minimum=33.539922 found after 18 iterations
************************************************
Found the Minimum!
************************************************
iter = 18, f = 33.539922, it took 0 s elapsed time
before normalization, pove = 22.5978 nove = 1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 7 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 33.5399
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 7 s.
p5n6 variation ends
p6n6 variation starts
--------------------------------------------------------------
PVPC3.0
--------------------------------------------------------------
Mode Options
(v) variation: optimize a pair condensate |PC> = | (P_pi)^{Np/2} (P_nu)^{Nn/2} >, so that < PC | H | PC > / <PC|PC> gets a minimum.
(CrankLinear) cranked variation with a linear constraint: with given omega, optimize a pair condensate, so that <PC|H|PC> - omega * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz.
(CrankQuadratic) cranked variation with a quadratic constraint: with given ALPHA, cons, optimize a pair condensate, so that <PC|H|PC> + ALPHA * <PC|ConsOP|PC> gets a minimum, ConsOP = Jx or Jz
(projection-quadrature) project states out of condensates by quadrature.
(LAP) Linear Algebraic Projection.
(pv) Projection after Variation of pair condensate(s) .
(DGTSR) Double Gamow-Teller transition sum rule.
(u) unit tests of basic components.
(ramdom-M0-pair) gadget: make a M=0 random pair.
(CG) gadget: calculate a Clebsch-Gorden coefficient.
(M-orbits) gadget: with given j-scheme orbits, show m-scheme orbits.
(Pandya) gadget: with given particle numbers and valence space, and interactions in pn format, do pandya transformation.
(projection-wo-Pandya) project states out of one condensate, without doing implicit Pandya particle-hole transformation.
(pv-wo-Pandya) projection after variation of pair condensates, without doing implicit Pandya particle-hole transformation.
(tensor-decomposition) With a given pair condensate, decompose it into parts with good angular momentum (J,M), and print percentages.
(EMtrans) With given one-body pp/nn transition operator, 1 pair condensate, and projected wavefunctions, calculate the reduced E&M transition probability.
So mode = v
Which nucleus do we study? (this affects nothing, only to mark) nucleus =
nucleus = p6n6
How many pairs are in proton condensate? pN =
How many pairs are in neutron condensate? nN =
So there are 6 pairs of protons, and 6 pairs of neutrons.
Please input the file name for single particle orbit info.
So j-scheme single particle orbits will be read from file: example/sd/sd.jsp
Detected j-scheme space info:
3 proton j-scheme orbits, and 3 neutron j-scheme orbits
Acore = 16 Zcore = 8
j-scheme orbs :
0-th jorb: n=1 l=2 2j=3 tz=-1
1-th jorb: n=1 l=2 2j=5 tz=-1
2-th jorb: n=2 l=0 2j=1 tz=-1
3-th jorb: n=1 l=2 2j=3 tz=1
4-th jorb: n=1 l=2 2j=5 tz=1
5-th jorb: n=2 l=0 2j=1 tz=1
proton : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=-1
1-th morb: n=1 l=2 j=3 m=-1 tz=-1
2-th morb: n=1 l=2 j=3 m=1 tz=-1
3-th morb: n=1 l=2 j=3 m=3 tz=-1
4-th morb: n=1 l=2 j=5 m=-5 tz=-1
5-th morb: n=1 l=2 j=5 m=-3 tz=-1
6-th morb: n=1 l=2 j=5 m=-1 tz=-1
7-th morb: n=1 l=2 j=5 m=1 tz=-1
8-th morb: n=1 l=2 j=5 m=3 tz=-1
9-th morb: n=1 l=2 j=5 m=5 tz=-1
10-th morb: n=2 l=0 j=1 m=-1 tz=-1
11-th morb: n=2 l=0 j=1 m=1 tz=-1
neutron : 12 m-scheme orbits:
0-th morb: n=1 l=2 j=3 m=-3 tz=1
1-th morb: n=1 l=2 j=3 m=-1 tz=1
2-th morb: n=1 l=2 j=3 m=1 tz=1
3-th morb: n=1 l=2 j=3 m=3 tz=1
4-th morb: n=1 l=2 j=5 m=-5 tz=1
5-th morb: n=1 l=2 j=5 m=-3 tz=1
6-th morb: n=1 l=2 j=5 m=-1 tz=1
7-th morb: n=1 l=2 j=5 m=1 tz=1
8-th morb: n=1 l=2 j=5 m=3 tz=1
9-th morb: n=1 l=2 j=5 m=5 tz=1
10-th morb: n=2 l=0 j=1 m=-1 tz=1
11-th morb: n=2 l=0 j=1 m=1 tz=1
Totally there are 12 proton m-scheme orbits and 12 neutron m-scheme orbits.
Interaction format = ? (jPQ/xpn)
Please input file name for interactions.
Scale the 2body interactions, or not? (y/n/h)
(y) automatical scaling : pow( (Acore+2)/A, 0.3 )
(n) no scaling
(h) handly input the scaling factor, which will be asked for later
So, shell model interactions in pn format will be read from example/sd/sd.xpn with scaling
How many pair condensates would you like to optimize?
kk = So totally we'll optimize 1 pair condensates.
Does the variation start from a random start? (y/n)
---------------------------------------------------------------------
--------------- variation2 starts --------------------------------
xpnint::read has read single particle energies : 2.1117 -3.9257 -3.2079 2.1117 -3.9257 -3.2079
Detected 159 lines of 2body interaction.
xpnint::read has read 30 Vpp's, 30 Vnn's, 99 Vpn's
finished reduction of Vpn terms, there are 98 Vpn's left.
Sequences of j-orbits are already sorted, so that a <= b, c <=d, (a,b) <= (c,d)
After scaling and optionally Pandya, new interactions are generated in pn format, written into output/2use.pn
before GnrPQint
after generating GnrPQint
P+Q interactions are generated, written into output/temp.jPQ
Finished pre-processing of the xpn interactions, xpn interactions after Pandya-transformation and scaling are written into output/2use.pn, equivalent J-scheme P+Q interactions are written into output/temp.jPQ
pNN=0 nNN=0 pNorb=12 nNorb=12 kk=1
-----------------------------------------------------------
varying the 0-th pair
-----------------------------------------------------------
The variation will start from a random position, with pseudo random numbers initialized with seconds from 1970/01/01/0:00
finished setting pystart, nystart
passed py[0-k], ny[0-k] to params
passed Ove, S, B, H to params
parameter number = 132
I jump into main loop of GSL_fdf_minimizer
before normalization, pove = 1 nove = 1
after normalization, pove = 1 nove = 1
---------- variation ends, it took 1 s. ------
--------------------------------------------------------------------------
The 0-th optimized orthogonal ket has < H > = 0
The optimized pairs are written into:
output/best-proton.pair0, output/best-neutron.pair0,
time = 1 s.
p6n6 variation ends
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。