XRB generator for helium burst setup

Historic generator file for APPROX-only XRB simulation of $\mathrm{He}$-accretor.

Generator

  1c This generator only runs APPROX
  2c
  3net 1 h1 he3 he4 n14 c12 o16 ne20 mg24
  4net 1 si28 s32 ar36 ca40 ti44 cr48 fe52
  5net 1 ni56 fe54 pn1 nt1
  6m nstar 1.00 fe54
  7c He star abundances
  8m acret 1.0 he4 0.0 c12 0.0 o16
  9c m acret 0.3 he4 0.2 c12 0.5 o16
 10c m acret 0. h1 0.980 he4 0. c12 0.02 n14 0. o16
 11c THIS GRID FOR He ACCRETION
 12g 0   2.0000e25  1 nstar  4.0e+8  1.0e+9
 13g 1   1.9000e25  1 nstar  4.0e+8  1.0e+9
 14g 40  1.0000e22  1 nstar  4.0e+8  1.0e+8
 15g 50  1.0000e21  1 nstar  4.0e+8  1.0e+8
 16g 51  8.0000e20  1 acret  2.0e+8  1.0e+8
 17g 54  2.0000e20  1 acret  1.0e+8  1.0e+6
 18g 55  0.         1 acret  5.0e+7  1.0e+4
 19dstat
 20c genburn  rphebg
 21p dtnew 1.e-4
 22p maxit 40
 23p dtcr .05
 24p dtct .05
 25p dtcd .10
 26p dtcq .1
 27p dtcdt .99
 28p nstop 1000000
 29p nedit 100000
 30p ndump 10
 31p ipup 2
 32p dyemult 50.
 33p dyqmult 2.
 34p dtcp .15
 35p yfloorx 3.e-3
 36p cenu 1.
 37p t7peek 1.e+50
 38p dtcut .1
 39p tfcrbu 2.
 40p tfctbu 10.
 41p fracneut .05
 42p radius0 1.0e+06
 43p summ0 2.8e+33
 44p xlum0 1.6e+34
 45p tnucmin 1.0e+99
 46p fclmax 1.e+99
 47p fclbu 1.e+99
 48p dtsmult 1.e+99
 49p dnratmax .25
 50p rnmin 1.e+6
 51p tnmin 1.e+4
 52p dnmin 2.e-5
 53p izonef 0
 54p idzonef 0
 55p jshell0 51
 56p rnmax 1.e+14
 57p tqsemin 3.e+9
 58p irtype 6
 59p fracrz1 .33
 60p fracrz2 .5
 61p abarratm 1.3
 62p frcsound .0
 63p convlim .0
 64p woversht 0.
 65p fmax0 .01
 66p fmax1 .015
 67p fmax2 .03
 68p nsdump 10
 69p iautoout 5
 70p iflgabar 0
 71p siqselim .02
 72p yfloorbx .003
 73p iwinsize 14001000
 74p wilsonmt -1.
 75p iold 1
 76p kaptab 0
 77p btempmin 1.1e8
 78p ncycqq 100000
 79p mazful -1
 80c p alpth 0.
 81p npixedit 1
 82
 83p h1hdep -1.d0
 84p he4hedep -1.d0
 85p si28dep -1.d0
 86p npixedit 1
 87p lburn 1
 88p iplotb 2
 89p bmasslow 2.80000019999895D33
 90p convlim 1.
 91p frcsound 1.
 92p btempmin 1.d7
 93p tnucmin 1.d7
 94p accrate 1.75d-9
 95p minzone 51
 96p jp0 40
 97p irtype 4
 98p zonemmin 1.5d19
 99p zonemmax 1.d20
100p lcout 10
101p ncnvout 1
102p fcrmax 1.d-8
103p fctmax 1.d-8
104p abunlim 1.d-4
105
106p nadapb 1
107c=======================================================================
108c Now follows the command file
109c=======================================================================
110//*
111c ------- 4U 0614+091 - --------
112c base heating (MeV/nuclen)
113o xqheat 0.45 def
114c accretion rate relative to Eddington
115o xledd 0.1 def
116c hydrogen mass fraction
117o xhfrac 0. def
118c .... accretion rate 1.75D-8 (L/Ledd) * 1.7/(X + 1)
119p accrate 1.75D-8
120c .... substrate luminosity - accrate * 6.0816737e+43 * (Q/MeV)
121c .... 1.0642929e+36 (L/Ledd) * (Q/MeV)
122p xlum0 1.0642929e+36
123c -------------------------
124c ..... SCALE to He/C/O L_Edd accretion
125o xeddf {1.7/(1.+xhfrac)} def
126p accrate {xeddf} *
127p xlum0 {xeddf} *
128c .......
129c substrate L, Q/MeV
130p xlum0 {xqheat} *
131c .......
132c set fraction of Eddington accretion rate
133p accrate {xledd} *
134p xlum0 {xledd} *
135
136c -------------------------
137c get model in equilibrium
138p ncnvout 0
139p nstop 1000000000
140p tnucmin 1.d10
141p tnumin 1.d7
142p accmass 1.d13
143p optconv 0.67
144p iaccadv 0
145p iacceadv 1
146c for APPROX ONLY
147p jp0 0
148p irtype 4
149plot
150c .........................
151c MATCH TO accmass
152@xm(jm)<1.01d18
153p maxbak 20
154p accdepth 1.d99
155p iterbarm 999999
156c .........................
157@time>1.d17
158zerotime
159p toffset 0.
160setcycle 0
161cutbin
162resetacc
163p lburn 0
164p dtnew 1.
165p iaccadv 1
166c use accdepth 5.d20 for He
167c use accdepth 1.d20 for H
168p accdepth 0.5d21
169p tnucmin 1.d7
170p izonef 1
171p idzonef 1
172p iazonef 0
173p zonermax 10.
174p zonemmax 1.d99
175p ddmin 1.d4
176c decretion
177p decrate -1.D0
178p idecmode 1
179p jshell0 0
180c use y coordinate for zoning
181c p fmax0 1.
182c p fmax1 1.
183c p fmax2 1.
184c p zonemmax 1.d99
185c p zonermax 1.d99
186c p zoneymax 0.02
187c
188c some other stuff
189p pbound 5.d18
190c
191c rotation
192c p centmult .667
193c 600 Hz rotation?
194c p angjacc 2.51327d15
195c p nangdis 1
196c p magnet 8
197c p ipixtype 31700
198p irtype 11
199c p xmagfbr 1.d-10
200c p xmagfbt 1.d-5