XRB generator for helium burst setup¶
Todo
add details on simulation
c This generator only runs APPROX
c
net 1 h1 he3 he4 n14 c12 o16 ne20 mg24
net 1 si28 s32 ar36 ca40 ti44 cr48 fe52
net 1 ni56 fe54 pn1 nt1
m nstar 1.00 fe54
c He star abundances
m acret 1.0 he4 0.0 c12 0.0 o16
c m acret 0.3 he4 0.2 c12 0.5 o16
c m acret 0. h1 0.980 he4 0. c12 0.02 n14 0. o16
c THIS GRID FOR He ACCRETION
g 0 2.0000e25 1 nstar 4.0e+8 1.0e+9
g 1 1.9000e25 1 nstar 4.0e+8 1.0e+9
g 40 1.0000e22 1 nstar 4.0e+8 1.0e+8
g 50 1.0000e21 1 nstar 4.0e+8 1.0e+8
g 51 8.0000e20 1 acret 2.0e+8 1.0e+8
g 54 2.0000e20 1 acret 1.0e+8 1.0e+6
g 55 0. 1 acret 5.0e+7 1.0e+4
dstat
c genburn rphebg
p 1 1.e-4
p 5 40
p 6 .05
p 7 .05
p 8 .10
p 9 .1
p 10 .99
p 14 1000000
p 16 100000
p 18 10
p 28 2
p 39 50.
p 40 2.
p 46 .15
p 47 3.e-3
p 48 1.
p 49 1.e+50
p 53 .1
p 54 2.
p 55 10.
p 59 .05
p 60 1.0e+06
p 61 2.8e+33
p 62 1.6e+34
p 65 1.0e+99
p 70 1.e+99
p 73 1.e+99
p 75 1.e+99
p 80 .25
p 82 1.e+6
p 83 1.e+4
p 84 2.e-5
p 86 0
p 87 0
p 93 51
p 88 1.e+14
p 105 3.e+9
p 132 6
p 138 .33
p 139 .5
p 144 1.3
p 146 .0
p 147 .0
p 148 0.
p 150 .01
p 151 .015
p 152 .03
p 156 10
p 159 5
p 160 0
p 189 .02
p 206 .003
p 42 14001000
p 199 -1.
p 388 1
p 377 0
p 233 1.1e8
p 299 100000
p 265 -1
c p 425 0.
p 64 1
p 405 -1.d0
p 406 -1.d0
p 420 -1.d0
p 64 1
p 434 1
p 443 2
p 419 2.80000019999895D33
p 147 1.
p 146 1.
p 233 1.d7
p 65 1.d7
p 211 1.75d-9
p 444 51
p 119 40
p 132 4
p 336 1.5d19
p 445 1.d20
p 437 10
p 376 1
p 11 1.d-8
p 12 1.d-8
p 128 1.d-4
p 137 1
c=======================================================================
c Now follows the command file
c=======================================================================
//*
c ------- 4U 0614+091 - --------
c base heating (MeV/nuclen)
o xqheat 0.45 def
c accretion rate relative to Eddington
o xledd 0.1 def
c hydrogen mass fraction
o xhfrac 0. def
c .... accretion rate 1.75D-8 (L/Ledd) * 1.7/(X + 1)
p accrate 1.75D-8
c .... substrate luminosity - accrate * 6.0816737e+43 * (Q/MeV)
c .... 1.0642929e+36 (L/Ledd) * (Q/MeV)
p xlum0 1.0642929e+36
c -------------------------
c ..... SCALE to He/C/O L_Edd accretion
o xeddf {1.7/(1.+xhfrac)} def
p accrate {xeddf} *
p xlum0 {xeddf} *
c .......
c substrate L, Q/MeV
p xlum0 {xqheat} *
c .......
c set fraction of Eddington accretion rate
p accrate {xledd} *
p xlum0 {xledd} *
c -------------------------
c get model in equilibrium
p ncnvout 0
p nstop 1000000000
p tnucmin 1.d10
p tnumin 1.d7
p accmass 1.d13
p optconv 0.67
p 521 0
p 520 1
c for APPROX ONLY
p jp0 0
p 132 4
plot
c .........................
c MATCH TO accmass
@xm(jm)<1.01d18
p 52 20
p accdepth 1.d99
p iterbarm 999999
c .........................
@time>1.d17
zerotime
p toffset 0.
setcycle 0
cutbin
resetacc
p lburn 0
p 1 1.
p 521 1
c use accdepth 5.d20 for He
c use accdepth 1.d20 for H
p accdepth 0.5d21
p tnucmin 1.d7
p 86 1
p 87 1
p 452 0
p 458 10.
p zonemmax 1.d99
p ddmin 1.d4
c decretion
p decrate -1.D0
p idecmode 1
p jshell0 0
c use y coordinate for zoning
c p 150 1.
c p 151 1.
c p 152 1.
c p 445 1.d99
c p 458 1.d99
c p 513 0.02
c
c some other stuff
p 69 5.d18
c
c rotation
c p 480 .667
c 600 Hz rotation?
c p 512 2.51327d15
c p 551 1
c p magnet 8
c p 113 31700
p 132 11
c p 552 1.d-10
c p 553 1.d-5