cosmic-popsynth 3.6.2__cp313-cp313-macosx_14_0_arm64.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- cosmic/.dylibs/libgcc_s.1.1.dylib +0 -0
- cosmic/.dylibs/libgfortran.5.dylib +0 -0
- cosmic/.dylibs/libquadmath.0.dylib +0 -0
- cosmic/Match.py +191 -0
- cosmic/__init__.py +32 -0
- cosmic/_commit_hash.py +1 -0
- cosmic/_evolvebin.cpython-313-darwin.so +0 -0
- cosmic/_version.py +1 -0
- cosmic/bse_utils/__init__.py +18 -0
- cosmic/bse_utils/zcnsts.py +570 -0
- cosmic/bse_utils/zdata.py +596 -0
- cosmic/checkstate.py +128 -0
- cosmic/data/cosmic-settings.json +1635 -0
- cosmic/evolve.py +607 -0
- cosmic/filter.py +214 -0
- cosmic/get_commit_hash.py +15 -0
- cosmic/output.py +466 -0
- cosmic/plotting.py +680 -0
- cosmic/sample/__init__.py +26 -0
- cosmic/sample/cmc/__init__.py +18 -0
- cosmic/sample/cmc/elson.py +411 -0
- cosmic/sample/cmc/king.py +260 -0
- cosmic/sample/initialbinarytable.py +251 -0
- cosmic/sample/initialcmctable.py +449 -0
- cosmic/sample/sampler/__init__.py +25 -0
- cosmic/sample/sampler/cmc.py +418 -0
- cosmic/sample/sampler/independent.py +1252 -0
- cosmic/sample/sampler/multidim.py +882 -0
- cosmic/sample/sampler/sampler.py +130 -0
- cosmic/test_evolve.py +108 -0
- cosmic/test_match.py +30 -0
- cosmic/test_sample.py +580 -0
- cosmic/test_utils.py +198 -0
- cosmic/utils.py +1574 -0
- cosmic_popsynth-3.6.2.data/scripts/cosmic-pop +544 -0
- cosmic_popsynth-3.6.2.dist-info/METADATA +55 -0
- cosmic_popsynth-3.6.2.dist-info/RECORD +38 -0
- cosmic_popsynth-3.6.2.dist-info/WHEEL +6 -0
|
@@ -0,0 +1,1635 @@
|
|
|
1
|
+
[
|
|
2
|
+
{
|
|
3
|
+
"category": "filters",
|
|
4
|
+
"category_label": "Filters",
|
|
5
|
+
"category_description": "Settings that filter the data returned by COSMIC simulations",
|
|
6
|
+
"docs-colour": "cadetblue",
|
|
7
|
+
"settings": [
|
|
8
|
+
{
|
|
9
|
+
"name": "binary_state",
|
|
10
|
+
"description": "Filter for the final state of the binaries you wish to retain",
|
|
11
|
+
"type": "checkbox",
|
|
12
|
+
"options-preface": "Each binary system will end its evolution in one of three states. Use these options to retain only certain endstates.",
|
|
13
|
+
"options": [
|
|
14
|
+
{
|
|
15
|
+
"name": 0,
|
|
16
|
+
"description": "Retain only binaries still alive today",
|
|
17
|
+
"default": true
|
|
18
|
+
},
|
|
19
|
+
{
|
|
20
|
+
"name": 1,
|
|
21
|
+
"description": "Retains binaries that merged",
|
|
22
|
+
"default": true
|
|
23
|
+
},
|
|
24
|
+
{
|
|
25
|
+
"name": 2,
|
|
26
|
+
"description": "Retains binaries that were disrupted"
|
|
27
|
+
}
|
|
28
|
+
]
|
|
29
|
+
},
|
|
30
|
+
{
|
|
31
|
+
"name": "timestep_conditions",
|
|
32
|
+
"description": "Pick specific time resolutions to print at targeted stages of the binary evolution.",
|
|
33
|
+
"type": "string",
|
|
34
|
+
"options-preface": "This is used in conjunction with the [bse] section value dtp to determine the timestep resolution for printing to the bcm array. See the <a href='evolve/resolution.html'>related guide</a> for more information.",
|
|
35
|
+
"options": [
|
|
36
|
+
{
|
|
37
|
+
"name": "'dtp=None'",
|
|
38
|
+
"description": "Only the final time step is printed to the bcm array",
|
|
39
|
+
"default": true
|
|
40
|
+
},
|
|
41
|
+
{
|
|
42
|
+
"name": "'dtp=1.0'",
|
|
43
|
+
"description": "Use 1 Myr timesteps for all evolutionary stages"
|
|
44
|
+
},
|
|
45
|
+
{
|
|
46
|
+
"name": "[['binstate==0', 'dtp=1.0']]",
|
|
47
|
+
"description": "Use 1 Myr timesteps for binaries until they merge or disrupt"
|
|
48
|
+
}
|
|
49
|
+
]
|
|
50
|
+
}
|
|
51
|
+
]
|
|
52
|
+
},
|
|
53
|
+
{
|
|
54
|
+
"category": "sampling",
|
|
55
|
+
"category_label": "Sampling",
|
|
56
|
+
"category_description": "Settings that change how the initial binary population is sampled",
|
|
57
|
+
"docs-colour": "darkseagreen",
|
|
58
|
+
"settings": [
|
|
59
|
+
{
|
|
60
|
+
"name": "sampling_method",
|
|
61
|
+
"description": "Select which models to use to generate an initial sample of binary parameters at Zero Age Main Sequence",
|
|
62
|
+
"type": "dropdown",
|
|
63
|
+
"options-preface": "",
|
|
64
|
+
"options": [
|
|
65
|
+
{
|
|
66
|
+
"name": "independent",
|
|
67
|
+
"description": "Initialize binaries with independent parameter distributions for the primary mass, mass ratio, eccentricity, separation, and binary fraction",
|
|
68
|
+
"default": true
|
|
69
|
+
},
|
|
70
|
+
{
|
|
71
|
+
"name": "multidim",
|
|
72
|
+
"description": "Initialize binaries with multidimensional parameter distributions according to <a href='http://adsabs.harvard.edu/abs/2017ApJS..230...15M'>Moe & Di Stefano 2017</a>"
|
|
73
|
+
}
|
|
74
|
+
]
|
|
75
|
+
},
|
|
76
|
+
{
|
|
77
|
+
"name": "primary_model",
|
|
78
|
+
"description": "Model for sampling primary masses <br><i>[Only used when <code>sampling_method = independent</code>]</i>",
|
|
79
|
+
"type": "dropdown",
|
|
80
|
+
"options-preface": "",
|
|
81
|
+
"options": [
|
|
82
|
+
{
|
|
83
|
+
"name": "salpeter55",
|
|
84
|
+
"description": "Use the Salpeter 1955 IMF"
|
|
85
|
+
},
|
|
86
|
+
{
|
|
87
|
+
"name": "kroupa93",
|
|
88
|
+
"description": "Use the Kroupa 1993 IMF"
|
|
89
|
+
},
|
|
90
|
+
{
|
|
91
|
+
"name": "kroupa01",
|
|
92
|
+
"description": "Use the Kroupa 2001 IMF",
|
|
93
|
+
"default": true
|
|
94
|
+
}
|
|
95
|
+
]
|
|
96
|
+
},
|
|
97
|
+
{
|
|
98
|
+
"name": "porb_model",
|
|
99
|
+
"description": "Model for sampling orbital periods <br><i>[Only used when <code>sampling_method = independent</code>]</i>",
|
|
100
|
+
"type": "dropdown",
|
|
101
|
+
"options-preface": "",
|
|
102
|
+
"options": [
|
|
103
|
+
{
|
|
104
|
+
"name": "sana12",
|
|
105
|
+
"description": "Sample from power law orbital period between 0.15 < log(P/day) < 5.5 following <a href='https://ui.adsabs.harvard.edu/abs/2012Sci...337..444S/abstract'>Sana+2012</a>",
|
|
106
|
+
"default": true
|
|
107
|
+
},
|
|
108
|
+
{
|
|
109
|
+
"name": "log_uniform",
|
|
110
|
+
"description": "Sample semi-major axis flat in log space from RRLO < 0.5 up to \\(10^{5} {\\rm R_{\\odot}}\\) according to <a href='http://adsabs.harvard.edu/abs/1983ARA%26A..21..343A'>Abt (1983)</a> and consistent with Dominik+2012,2013 - then convert to orbital period in days using Kepler's third law."
|
|
111
|
+
},
|
|
112
|
+
{
|
|
113
|
+
"name": "renzo19",
|
|
114
|
+
"description": "Uses <code>sana12</code> for massive binaries (\\(m_1 > 15 {\\rm M_{\\odot}}\\)) and flat in log otherwise (following <a href='https://ui.adsabs.harvard.edu/abs/2019A%26A...624A..66R/abstract'>Renzo+19</a>)."
|
|
115
|
+
},
|
|
116
|
+
{
|
|
117
|
+
"name": "raghavan10",
|
|
118
|
+
"description": "Sample log normal orbital periods in days with mean_logP = 4.9 and sigma_logP = 2.3 between \\(0 < \\log_{10}(P / {\\rm day}) < 9\\) following <a href='https://ui.adsabs.harvard.edu/abs/2010ApJS..190....1R/abstract'>Raghavan+2010 </a>"
|
|
119
|
+
},
|
|
120
|
+
{
|
|
121
|
+
"name": "moe19",
|
|
122
|
+
"description": "As <code>raghavan10</code> but with different close binary fractions following <a href='https://ui.adsabs.harvard.edu/abs/2019ApJ...875...61M/abstract'>Moe+2019</a>"
|
|
123
|
+
},
|
|
124
|
+
{
|
|
125
|
+
"name": "martinez26",
|
|
126
|
+
"description": "Uses sana12 for massive binaries (\\(m_1 ≥ 8 {\\rm M_{\\odot}}\\)) with upper limit at 3000 days and raghavan10 for low-mass binaries. Used in <a href='https://ui.adsabs.harvard.edu/abs/2025arXiv251123285M/abstract'>Martinez+2026</a>."
|
|
127
|
+
},
|
|
128
|
+
{
|
|
129
|
+
"name": "custom",
|
|
130
|
+
"description": "Sample from a custom power law. The user provides a dictionary of min, max and slope values for the power law."
|
|
131
|
+
}
|
|
132
|
+
]
|
|
133
|
+
},
|
|
134
|
+
{
|
|
135
|
+
"name": "ecc_model",
|
|
136
|
+
"description": "Model for sampling eccentricity <br><i>[Only used when <code>sampling_method = independent</code>]</i>",
|
|
137
|
+
"type": "dropdown",
|
|
138
|
+
"options-preface": "",
|
|
139
|
+
"options": [
|
|
140
|
+
{
|
|
141
|
+
"name": "thermal",
|
|
142
|
+
"description": "Samples from a thermal eccentricity distribution following <a href='http://adsabs.harvard.edu/abs/1975MNRAS.173..729H'>Heggie (1975)</a>"
|
|
143
|
+
},
|
|
144
|
+
{
|
|
145
|
+
"name": "uniform",
|
|
146
|
+
"description": "Samples from a uniform eccentricity distribution"
|
|
147
|
+
},
|
|
148
|
+
{
|
|
149
|
+
"name": "sana12",
|
|
150
|
+
"description": "Samples from the eccentricity distribution from <a href='https://ui.adsabs.harvard.edu/abs/2012Sci...337..444S/abstract'>Sana+2012</a>",
|
|
151
|
+
"default": true
|
|
152
|
+
},
|
|
153
|
+
{
|
|
154
|
+
"name": "circular",
|
|
155
|
+
"description": "Assumes zero eccentricity for all systems"
|
|
156
|
+
}
|
|
157
|
+
]
|
|
158
|
+
},
|
|
159
|
+
{
|
|
160
|
+
"name": "qmin",
|
|
161
|
+
"description": "Minimum mass ratio for sampling the secondary mass <br><i>[Only used when <code>sampling_method = independent</code>]</i>",
|
|
162
|
+
"type": "number",
|
|
163
|
+
"options-preface": "The assumed mass ratio distribution is flat in \\(q \\equiv m_2 / m_1\\). NOTE: only one of <code>qmin</code> and <code>m2_min</code> should be specified.",
|
|
164
|
+
"options": [
|
|
165
|
+
{
|
|
166
|
+
"name": "values in [0, 1]",
|
|
167
|
+
"description": "Sets the minimum mass ratio"
|
|
168
|
+
},
|
|
169
|
+
{
|
|
170
|
+
"name": -1,
|
|
171
|
+
"description": "Set the minimum mass ratio such that the pre-MS lifetime of the secondary is not longer than the full lifetime of the primary if it were to evolve as a single star",
|
|
172
|
+
"default": true
|
|
173
|
+
}
|
|
174
|
+
]
|
|
175
|
+
},
|
|
176
|
+
{
|
|
177
|
+
"name": "m2_min",
|
|
178
|
+
"description": "Minimum secondary mass for sampling <br><i>[Only used when <code>sampling_method = independent</code>]</i>",
|
|
179
|
+
"type": "string",
|
|
180
|
+
"options-preface": "NOTE: only one of <code>qmin</code> and <code>m2_min</code> should be specified.",
|
|
181
|
+
"options": [
|
|
182
|
+
{
|
|
183
|
+
"name": "positive values",
|
|
184
|
+
"description": "Sample the secondary mass uniformly between <code>m2_min</code> and <code>mass_1</code>"
|
|
185
|
+
},
|
|
186
|
+
{
|
|
187
|
+
"name": "None",
|
|
188
|
+
"description": "Default value (ignore, use <code>qmin</code>)",
|
|
189
|
+
"default": true
|
|
190
|
+
}
|
|
191
|
+
]
|
|
192
|
+
},
|
|
193
|
+
{
|
|
194
|
+
"name": "binfrac_model",
|
|
195
|
+
"description": "Model for sampling binary fraction <br><i>[Only used when <code>sampling_method = independent</code>]</i>",
|
|
196
|
+
"type": "string",
|
|
197
|
+
"options-preface": "",
|
|
198
|
+
"options": [
|
|
199
|
+
{
|
|
200
|
+
"name": "values between [0, 1]",
|
|
201
|
+
"description": "Fixed binary fraction"
|
|
202
|
+
},
|
|
203
|
+
{
|
|
204
|
+
"name": "vanHaaften",
|
|
205
|
+
"description": "Primary mass dependent binary fraction following van Haaften+05"
|
|
206
|
+
},
|
|
207
|
+
{
|
|
208
|
+
"name": "offner23",
|
|
209
|
+
"description": "Primary mass dependent binary fraction following Offner+23"
|
|
210
|
+
},
|
|
211
|
+
{
|
|
212
|
+
"name": 0.5,
|
|
213
|
+
"description": "Default value",
|
|
214
|
+
"default": true
|
|
215
|
+
}
|
|
216
|
+
]
|
|
217
|
+
},
|
|
218
|
+
{
|
|
219
|
+
"name": "SF_start",
|
|
220
|
+
"description": "Sets the time in the past when star formation initiates in Myr.",
|
|
221
|
+
"type": "number",
|
|
222
|
+
"options-preface": "",
|
|
223
|
+
"options": [
|
|
224
|
+
{
|
|
225
|
+
"name": "positive values",
|
|
226
|
+
"description": "Start time of star formation in Myr"
|
|
227
|
+
},
|
|
228
|
+
{
|
|
229
|
+
"name": 13700.0,
|
|
230
|
+
"description": "For example, this specifies a start time at the beginning of a Hubble time",
|
|
231
|
+
"default": true
|
|
232
|
+
}
|
|
233
|
+
]
|
|
234
|
+
},
|
|
235
|
+
{
|
|
236
|
+
"name": "SF_duration",
|
|
237
|
+
"description": "Sets the duration of constant star formation from ``SF_start`` in Myr.",
|
|
238
|
+
"type": "number",
|
|
239
|
+
"options-preface": "",
|
|
240
|
+
"options": [
|
|
241
|
+
{
|
|
242
|
+
"name": "positive values",
|
|
243
|
+
"description": "Duration of star formation in Myr"
|
|
244
|
+
},
|
|
245
|
+
{
|
|
246
|
+
"name": 0.0,
|
|
247
|
+
"description": "A single burst of star formation",
|
|
248
|
+
"default": true
|
|
249
|
+
},
|
|
250
|
+
{
|
|
251
|
+
"name": 13700.0,
|
|
252
|
+
"description": "For example, this specifies a constant star formation rate over a Hubble time"
|
|
253
|
+
}
|
|
254
|
+
]
|
|
255
|
+
|
|
256
|
+
},
|
|
257
|
+
{
|
|
258
|
+
"name": "metallicity",
|
|
259
|
+
"description": "Sets the metallicity of the stellar population.",
|
|
260
|
+
"type": "number",
|
|
261
|
+
"options-preface": "COSMIC expects an <b>absolute</b> metallicity (i.e., <b>NOT</b> units of <i>zsun</i>)",
|
|
262
|
+
"options": [
|
|
263
|
+
{
|
|
264
|
+
"name": "positive values",
|
|
265
|
+
"description": "Absolute metallicity"
|
|
266
|
+
},
|
|
267
|
+
{
|
|
268
|
+
"name": 0.02,
|
|
269
|
+
"description": "For example, this sets the metallicity to approximately solar metallicity",
|
|
270
|
+
"default": true
|
|
271
|
+
}
|
|
272
|
+
]
|
|
273
|
+
}
|
|
274
|
+
]
|
|
275
|
+
},
|
|
276
|
+
{
|
|
277
|
+
"category": "convergence",
|
|
278
|
+
"category_label": "Convergence",
|
|
279
|
+
"category_description": "Settings that control the convergence of the simulation run with cosmic-pop",
|
|
280
|
+
"docs-colour": "#a873dd",
|
|
281
|
+
"settings": [
|
|
282
|
+
{
|
|
283
|
+
"name": "convergence_params",
|
|
284
|
+
"description": "A list of parameters you would like to verify have converged to a single distribution shape when running cosmic-pop from the command line.",
|
|
285
|
+
"type": "checkbox",
|
|
286
|
+
"options-preface": "",
|
|
287
|
+
"options": [
|
|
288
|
+
{
|
|
289
|
+
"name": "mass_1",
|
|
290
|
+
"description": "Primary mass",
|
|
291
|
+
"default": true
|
|
292
|
+
},
|
|
293
|
+
{
|
|
294
|
+
"name": "mass_2",
|
|
295
|
+
"description": "Secondary mass",
|
|
296
|
+
"default": true
|
|
297
|
+
},
|
|
298
|
+
{
|
|
299
|
+
"name": "sep",
|
|
300
|
+
"description": "Separation"
|
|
301
|
+
},
|
|
302
|
+
{
|
|
303
|
+
"name": "porb",
|
|
304
|
+
"description": "Orbital period",
|
|
305
|
+
"default": true
|
|
306
|
+
},
|
|
307
|
+
{
|
|
308
|
+
"name": "ecc",
|
|
309
|
+
"description": "Eccentricity",
|
|
310
|
+
"default": true
|
|
311
|
+
},
|
|
312
|
+
{
|
|
313
|
+
"name": "massc_1",
|
|
314
|
+
"description": "Primary core mass"
|
|
315
|
+
},
|
|
316
|
+
{
|
|
317
|
+
"name": "massc_2",
|
|
318
|
+
"description": "Secondary core mass"
|
|
319
|
+
},
|
|
320
|
+
{
|
|
321
|
+
"name": "rad_1",
|
|
322
|
+
"description": "Primary radius"
|
|
323
|
+
},
|
|
324
|
+
{
|
|
325
|
+
"name": "rad_2",
|
|
326
|
+
"description": "Secondary radius"
|
|
327
|
+
}
|
|
328
|
+
]
|
|
329
|
+
},
|
|
330
|
+
{
|
|
331
|
+
"name": "convergence_limits",
|
|
332
|
+
"description": "Specifies limits for parameters included in the <code>convergence_params</code> list.",
|
|
333
|
+
"type": "string",
|
|
334
|
+
"options-preface": "For each parameter specified <code>convergence_limits</code>, the lower and upper limit must be included.",
|
|
335
|
+
"options": [
|
|
336
|
+
{
|
|
337
|
+
"name": "{}",
|
|
338
|
+
"description": "No limits specified",
|
|
339
|
+
"default": true
|
|
340
|
+
},
|
|
341
|
+
{
|
|
342
|
+
"name": "{'mass_1' : [5, 10], 'sep' : [0, 10]}",
|
|
343
|
+
"description": "For example, this specifies that the primary mass must be between 5 and 10 solar masses, and the separation must be between 0 and 10 Rsun."
|
|
344
|
+
}
|
|
345
|
+
]
|
|
346
|
+
},
|
|
347
|
+
{
|
|
348
|
+
"name": "pop_select",
|
|
349
|
+
"description": "Selects the stage of the evolution at which you would like to check for convergence.",
|
|
350
|
+
"type": "dropdown",
|
|
351
|
+
"options-preface": "This will filter for systems that satisfy the <code>final_kstar1</code> and <code>final_kstar2</code> selections from the command line call of cosmic-pop at the following states:",
|
|
352
|
+
"options": [
|
|
353
|
+
{
|
|
354
|
+
"name": "formation",
|
|
355
|
+
"description": "At binary formation",
|
|
356
|
+
"default": true
|
|
357
|
+
},
|
|
358
|
+
{
|
|
359
|
+
"name": "1_SN",
|
|
360
|
+
"description": "Just before the first supernova"
|
|
361
|
+
},
|
|
362
|
+
{
|
|
363
|
+
"name": "2_SN",
|
|
364
|
+
"description": "Just before the second supernova"
|
|
365
|
+
},
|
|
366
|
+
{
|
|
367
|
+
"name": "disruption",
|
|
368
|
+
"description": "Just before binary disruption"
|
|
369
|
+
},
|
|
370
|
+
{
|
|
371
|
+
"name": "final_state",
|
|
372
|
+
"description": "After the full evolution specified by the user-supplied evolution time"
|
|
373
|
+
},
|
|
374
|
+
{
|
|
375
|
+
"name": "XRB_form",
|
|
376
|
+
"description": "At the start of RLOF following the first supernova"
|
|
377
|
+
}
|
|
378
|
+
]
|
|
379
|
+
},
|
|
380
|
+
{
|
|
381
|
+
"name": "apply_convergence_limits",
|
|
382
|
+
"description": "Whether to filter the binary population (including the bcm, bpp, initC, and kick_info DataFrames) to only contain the binaries that satisfy the constraints from convergence_limits.",
|
|
383
|
+
"type": "dropdown",
|
|
384
|
+
"options-preface": "",
|
|
385
|
+
"options": [
|
|
386
|
+
{
|
|
387
|
+
"name": "True",
|
|
388
|
+
"description": "Filter the binary population to only contain the binaries that satisfy the constraints from convergence_limits"
|
|
389
|
+
},
|
|
390
|
+
{
|
|
391
|
+
"name": "False",
|
|
392
|
+
"description": "Do not filter the binary population",
|
|
393
|
+
"default": true
|
|
394
|
+
}
|
|
395
|
+
]
|
|
396
|
+
},
|
|
397
|
+
{
|
|
398
|
+
"name": "match",
|
|
399
|
+
"description": "Provides the tolerance for the convergence calculation and is calculated as \\({\\rm match} = \\log_{\\rm 10} (1 - {\\rm convergence})\\)",
|
|
400
|
+
"type": "number",
|
|
401
|
+
"options-preface": "",
|
|
402
|
+
"options": [
|
|
403
|
+
{
|
|
404
|
+
"name": -5.0,
|
|
405
|
+
"description": "For example, this specifies a tolerance of \\(10^{-5}\\) for convergence",
|
|
406
|
+
"default": true
|
|
407
|
+
}
|
|
408
|
+
]
|
|
409
|
+
}
|
|
410
|
+
]
|
|
411
|
+
},
|
|
412
|
+
{
|
|
413
|
+
"category": "rand_seed",
|
|
414
|
+
"category_label": "Random Seed",
|
|
415
|
+
"category_description": "Settings that control the random number generation used",
|
|
416
|
+
"docs-colour": "#f0ba80",
|
|
417
|
+
"settings": [
|
|
418
|
+
{
|
|
419
|
+
"name": "seed",
|
|
420
|
+
"description": "Sets the seed for the random number generator (for <code>numpy.random.seed()</code>)",
|
|
421
|
+
"type": "number",
|
|
422
|
+
"options-preface": "",
|
|
423
|
+
"options": [
|
|
424
|
+
{
|
|
425
|
+
"name": "integer values",
|
|
426
|
+
"description": "Random seed"
|
|
427
|
+
},
|
|
428
|
+
{
|
|
429
|
+
"name": 42,
|
|
430
|
+
"description": "For example, this sets the seed to 42",
|
|
431
|
+
"default": true
|
|
432
|
+
}
|
|
433
|
+
]
|
|
434
|
+
}
|
|
435
|
+
]
|
|
436
|
+
},
|
|
437
|
+
{
|
|
438
|
+
"category": "bse",
|
|
439
|
+
"category_label": "Binary physics",
|
|
440
|
+
"category_description": "Settings that control the binary physics in the simulation",
|
|
441
|
+
"docs-colour": "lightcoral",
|
|
442
|
+
"settings": [
|
|
443
|
+
{
|
|
444
|
+
"name": "pts1",
|
|
445
|
+
"description": "Sets the timestep modifier for main sequence stars (<code>dtp *= pts1</code>)",
|
|
446
|
+
"type": "number",
|
|
447
|
+
"settings-section": "Timesteps",
|
|
448
|
+
"options-preface": "",
|
|
449
|
+
"options": [
|
|
450
|
+
{
|
|
451
|
+
"name": "positive values",
|
|
452
|
+
"description": "Timestep modifier"
|
|
453
|
+
},
|
|
454
|
+
{
|
|
455
|
+
"name": 0.001,
|
|
456
|
+
"description": "Recommended value from <a href='https://ui.adsabs.harvard.edu/abs/2019arXiv190207718B/abstract'>Bannerjee+2019</a> for NS/BH progenitors",
|
|
457
|
+
"default": true
|
|
458
|
+
},
|
|
459
|
+
{
|
|
460
|
+
"name": 0.05,
|
|
461
|
+
"description": "Recommended value from <a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Hurley+2000</a>"
|
|
462
|
+
}
|
|
463
|
+
]
|
|
464
|
+
},
|
|
465
|
+
{
|
|
466
|
+
"name": "pts2",
|
|
467
|
+
"description": "Sets the timestep modifier for Giant Branch (GB, CHeB, AGB, HeGB) stars (<code>dtp *= pts2</code>)",
|
|
468
|
+
"type": "number",
|
|
469
|
+
"options-preface": "",
|
|
470
|
+
"options": [
|
|
471
|
+
{
|
|
472
|
+
"name": "positive values",
|
|
473
|
+
"description": "Timestep modifier"
|
|
474
|
+
},
|
|
475
|
+
{
|
|
476
|
+
"name": 0.01,
|
|
477
|
+
"description": "Recommended value from <a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Hurley+2000</a>",
|
|
478
|
+
"default": true
|
|
479
|
+
}
|
|
480
|
+
]
|
|
481
|
+
},
|
|
482
|
+
{
|
|
483
|
+
"name": "pts3",
|
|
484
|
+
"description": "Sets the timestep modifier for HG, HeMS stars (<code>dtp *= pts3</code>)",
|
|
485
|
+
"type": "number",
|
|
486
|
+
"options-preface": "",
|
|
487
|
+
"options": [
|
|
488
|
+
{
|
|
489
|
+
"name": "positive values",
|
|
490
|
+
"description": "Timestep modifier"
|
|
491
|
+
},
|
|
492
|
+
{
|
|
493
|
+
"name": 0.02,
|
|
494
|
+
"description": "Recommended value from <a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Hurley+2000</a>",
|
|
495
|
+
"default": true
|
|
496
|
+
}
|
|
497
|
+
]
|
|
498
|
+
},
|
|
499
|
+
{
|
|
500
|
+
"name": "zsun",
|
|
501
|
+
"description": "Sets the metallicity of the Sun which primarily affects stellar winds.",
|
|
502
|
+
"type": "number",
|
|
503
|
+
"settings-section": "Metallicity",
|
|
504
|
+
"options-preface": "Note that the wind prescriptions for OB stars are calibrated to zsun = 0.019 as described in <a href='https://ui.adsabs.harvard.edu/abs/2001A%26A...369..574V/abstract'>Vink+2001.</a>",
|
|
505
|
+
"options": [
|
|
506
|
+
{
|
|
507
|
+
"name": "positive values",
|
|
508
|
+
"description": "Set the solar metallicity"
|
|
509
|
+
},
|
|
510
|
+
{
|
|
511
|
+
"name": 0.014,
|
|
512
|
+
"description": "Following <a href='https://ui.adsabs.harvard.edu/abs/2009ARA%26A..47..481A/abstract'>Asplund+2009</a>",
|
|
513
|
+
"default": true
|
|
514
|
+
}
|
|
515
|
+
]
|
|
516
|
+
},
|
|
517
|
+
{
|
|
518
|
+
"name": "windflag",
|
|
519
|
+
"description": "Selects the model for wind mass loss for each star",
|
|
520
|
+
"type": "dropdown",
|
|
521
|
+
"settings-section": "Stellar Winds",
|
|
522
|
+
"options-preface": "",
|
|
523
|
+
"options": [
|
|
524
|
+
{
|
|
525
|
+
"name": 0,
|
|
526
|
+
"description": "Standard SSE/BSE (<a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Hurley+2000</a>)"
|
|
527
|
+
},
|
|
528
|
+
{
|
|
529
|
+
"name": 1,
|
|
530
|
+
"description": "StarTrack (<a href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>)"
|
|
531
|
+
},
|
|
532
|
+
{
|
|
533
|
+
"name": 2,
|
|
534
|
+
"description": "Metallicity dependence for O/B stars and Wolf Rayet stars (<a href='https://ui.adsabs.harvard.edu/abs/2001A%26A...369..574V/abstract'>Vink+2001</a>, <a href='https://ui.adsabs.harvard.edu/abs/2005A%26A...442..587V/abstract'>Vink+2005</a>)"
|
|
535
|
+
},
|
|
536
|
+
{
|
|
537
|
+
"name": 3,
|
|
538
|
+
"description": "Same as 2, but LBV-like mass loss for giants and non-degenerate stars beyond the Humphreys-Davidson limit",
|
|
539
|
+
"default": true
|
|
540
|
+
}
|
|
541
|
+
]
|
|
542
|
+
},
|
|
543
|
+
{
|
|
544
|
+
"name": "eddlimflag",
|
|
545
|
+
"description": "Adjusts the dependence of mass loss on metallicity for stars near the Eddington limit (see <a href='https://ui.adsabs.harvard.edu/abs/2011A%26A...535A..56G/abstract'>Grafener+2011</a>, <a href='https://ui.adsabs.harvard.edu/abs/2018MNRAS.474.2959G/abstract'>Giacobbo+2018</a>).",
|
|
546
|
+
"type": "dropdown",
|
|
547
|
+
"options-preface": "",
|
|
548
|
+
"options": [
|
|
549
|
+
{
|
|
550
|
+
"name": 0,
|
|
551
|
+
"description": "does not adjust metallicity dependence for stars near the Eddington limit",
|
|
552
|
+
"default": true
|
|
553
|
+
},
|
|
554
|
+
{
|
|
555
|
+
"name": 1,
|
|
556
|
+
"description": "adjusts metallicity dependence for stars near the Eddington limit as in <a href='https://ui.adsabs.harvard.edu/abs/2018MNRAS.474.2959G/abstract'>Giacobbo+2018</a>."
|
|
557
|
+
}
|
|
558
|
+
]
|
|
559
|
+
},
|
|
560
|
+
{
|
|
561
|
+
"name": "neta",
|
|
562
|
+
"description": "Reimers mass-loss coefficient (<a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Equation 106 of SSE</a>).",
|
|
563
|
+
"type": "number",
|
|
564
|
+
"options-preface": "Note: this equation has a typo. There is an extra \\(\\eta\\) out front; the correct rate is directly proportional to \\(\\eta\\). See also <a href='https://ui.adsabs.harvard.edu/abs/1978A%26A....70..227K/abstract'>Kurdritzki+1978, Section Vb</a> for discussion.",
|
|
565
|
+
"options": [
|
|
566
|
+
{
|
|
567
|
+
"name": "positive values",
|
|
568
|
+
"description": "Set \\(\\eta\\) value"
|
|
569
|
+
},
|
|
570
|
+
{
|
|
571
|
+
"name": 0.5,
|
|
572
|
+
"description": "Sets \\(\\eta\\) for Reimers mass loss to 0.5",
|
|
573
|
+
"default": true
|
|
574
|
+
}
|
|
575
|
+
]
|
|
576
|
+
},
|
|
577
|
+
{
|
|
578
|
+
"name": "bwind",
|
|
579
|
+
"description": "Binary enhanced mass loss parameter",
|
|
580
|
+
"type": "number",
|
|
581
|
+
"options-preface": "This parameter is used to enhance mass loss in binaries. See <a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Hurley+2000, Eq. 12</a> for more information.",
|
|
582
|
+
"options": [
|
|
583
|
+
{
|
|
584
|
+
"name": "positive values",
|
|
585
|
+
"description": "Sets the binary enhanced mass loss parameter, \\(B_w\\), from <a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Hurley+2000, Eq. 12</a>"
|
|
586
|
+
},
|
|
587
|
+
{
|
|
588
|
+
"name": 0.0,
|
|
589
|
+
"description": "Default value (no effect on single stars)",
|
|
590
|
+
"default": true
|
|
591
|
+
}
|
|
592
|
+
]
|
|
593
|
+
},
|
|
594
|
+
{
|
|
595
|
+
"name": "hewind",
|
|
596
|
+
"description": "Helium star mass loss parameter",
|
|
597
|
+
"type": "number",
|
|
598
|
+
"options-preface": "\\( 10^{-13} {\\rm \\ \\texttt{hewind} \\ } L^{2/3}\\) gives He star mass-loss. Equivalent to \\(1 - \\mu\\) in the last equation on <a href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Hurley+2000, page 19</a>.",
|
|
599
|
+
"options": [
|
|
600
|
+
{
|
|
601
|
+
"name": "range [0, 1]",
|
|
602
|
+
"description": "Sets the helium star mass loss parameter"
|
|
603
|
+
},
|
|
604
|
+
{
|
|
605
|
+
"name": 0.5,
|
|
606
|
+
"description": "Default value",
|
|
607
|
+
"default": true
|
|
608
|
+
}
|
|
609
|
+
]
|
|
610
|
+
},
|
|
611
|
+
{
|
|
612
|
+
"name": "beta",
|
|
613
|
+
"description": "Wind velocity factor. \\( v_{\\rm wind}^2 \\propto \\beta\\), see <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002, Eq. 9</a>.",
|
|
614
|
+
"type": "number",
|
|
615
|
+
"options-preface": "",
|
|
616
|
+
"options": [
|
|
617
|
+
{
|
|
618
|
+
"name": -1,
|
|
619
|
+
"description": "Follow StarTrack prescription for wind velocity factor, \\(\\beta_w\\), from <a href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>"
|
|
620
|
+
},
|
|
621
|
+
{
|
|
622
|
+
"name": "positive values",
|
|
623
|
+
"description": "Sets the wind velocity factor, \\(\\beta_w\\)"
|
|
624
|
+
},
|
|
625
|
+
{
|
|
626
|
+
"name": 0.125,
|
|
627
|
+
"description": "Default value",
|
|
628
|
+
"default": true
|
|
629
|
+
}
|
|
630
|
+
]
|
|
631
|
+
},
|
|
632
|
+
{
|
|
633
|
+
"name": "xi",
|
|
634
|
+
"description": "Wind accretion efficiency factor, which gives the fraction of angular momentum lost via winds from the primary that transfers to the spin angular momentum of the companion.",
|
|
635
|
+
"type": "number",
|
|
636
|
+
"options-preface": "Corresponds to \\(\\mu_w\\) in <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002, Eq. 11</a>",
|
|
637
|
+
"options": [
|
|
638
|
+
{
|
|
639
|
+
"name": "range [0, 1]",
|
|
640
|
+
"description": "Sets the wind accretion efficiency factor"
|
|
641
|
+
},
|
|
642
|
+
{
|
|
643
|
+
"name": 0.5,
|
|
644
|
+
"description": "Default value",
|
|
645
|
+
"default": true
|
|
646
|
+
}
|
|
647
|
+
]
|
|
648
|
+
},
|
|
649
|
+
{
|
|
650
|
+
"name": "acc2",
|
|
651
|
+
"description": "Bondi-Hoyle wind accretion factor where the mean wind accretion rate onto the secondary is proportional to <code>acc2</code>. See <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002, Eq. 6</a>.",
|
|
652
|
+
"type": "number",
|
|
653
|
+
"options-preface": "This value is equivalent to \\(\\alpha_w\\) in <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002, Eq. 6</a>.",
|
|
654
|
+
"options": [
|
|
655
|
+
{
|
|
656
|
+
"name": "positive values",
|
|
657
|
+
"description": "Sets the Bondi-Hoyle wind accretion factor, \\(\\alpha_w\\)"
|
|
658
|
+
},
|
|
659
|
+
{
|
|
660
|
+
"name": 1.5,
|
|
661
|
+
"description": "Default value",
|
|
662
|
+
"default": true
|
|
663
|
+
}
|
|
664
|
+
]
|
|
665
|
+
},
|
|
666
|
+
{
|
|
667
|
+
"settings-section": "Common-envelope",
|
|
668
|
+
"settings-section-description": "Note: there are cases where a common envelope is forced regardless of the critical mass ratio for unstable mass transfer. In the following cases, a common envelope occurs regardless of the choices below:<ul><li><b>contact</b>: the stellar radii go into contact (common for similar ZAMS systems)</li><li><b>periapse contact</b>: the periapse distance is smaller than either of the stellar radii (common for highly eccentric systems)</li><li><b>core Roche overflow</b>: either of the stellar radii overflow their component's Roche radius (in this case, mass transfer from the convective core is always dynamically unstable)</li></ul>",
|
|
669
|
+
"name": "alpha1",
|
|
670
|
+
"description": "Common-envelope efficiency parameter which scales the efficiency of transferring orbital energy to the envelope. See <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002, Eq. 71</a>.",
|
|
671
|
+
"type": "number",
|
|
672
|
+
"options-preface": "This value is equivalent to \\(\\alpha\\) in <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002, Eq. 71</a>.",
|
|
673
|
+
"options": [
|
|
674
|
+
{
|
|
675
|
+
"name": "positive values",
|
|
676
|
+
"description": "Sets the common-envelope efficiency parameter, \\(\\alpha\\)"
|
|
677
|
+
},
|
|
678
|
+
{
|
|
679
|
+
"name": 1.0,
|
|
680
|
+
"description": "Default value",
|
|
681
|
+
"default": true
|
|
682
|
+
}
|
|
683
|
+
]
|
|
684
|
+
},
|
|
685
|
+
{
|
|
686
|
+
"name": "lambdaf",
|
|
687
|
+
"description": "Binding energy factor for common envelope evolution. The initial binding energy of the stellar envelope is proportional to \\(1 / \\lambda\\). See <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002, Eq. 69</a>.",
|
|
688
|
+
"type": "number",
|
|
689
|
+
"options-preface": "",
|
|
690
|
+
"options": [
|
|
691
|
+
{
|
|
692
|
+
"name": "positive values",
|
|
693
|
+
"description": "uses variable lambda prescription detailed in appendix of <a href='https://ui.adsabs.harvard.edu/abs/2014A%26A...563A..83C/abstract'>Claeys+2014</a> where lambdaf is the fraction of the ionization energy that can go into ejecting the envelope; to use this prescription without extra ionization energy, set lambdaf = 0"
|
|
694
|
+
},
|
|
695
|
+
{
|
|
696
|
+
"name": 0.0,
|
|
697
|
+
"description": "As above, this is the default choice",
|
|
698
|
+
"default": true
|
|
699
|
+
},
|
|
700
|
+
{
|
|
701
|
+
"name": "negative values",
|
|
702
|
+
"description": "Uses a fixed value (i.e. fixes \\( \\lambda \\) to a value of -lambdaf)"
|
|
703
|
+
}
|
|
704
|
+
]
|
|
705
|
+
},
|
|
706
|
+
{
|
|
707
|
+
"name": "ceflag",
|
|
708
|
+
"description": "Selects the model to set the initial orbital energy (choose whether to use total mass of the stars instead of the core masses)",
|
|
709
|
+
"type": "dropdown",
|
|
710
|
+
"options-preface": "",
|
|
711
|
+
"options": [
|
|
712
|
+
{
|
|
713
|
+
"name": 0,
|
|
714
|
+
"description": "Use the core masses of the stars to set the initial orbital energy (as in <a href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002 Eq. 70</a>)"
|
|
715
|
+
},
|
|
716
|
+
{
|
|
717
|
+
"name": 1,
|
|
718
|
+
"description": "Use the total mass of the stars to set the initial orbital energy (as in <a href='https://ui.adsabs.harvard.edu/abs/1990ApJ...358..189D/abstract'>de Kool 1990</a>)",
|
|
719
|
+
"default": true
|
|
720
|
+
}
|
|
721
|
+
]
|
|
722
|
+
},
|
|
723
|
+
{
|
|
724
|
+
"name": "cekickflag",
|
|
725
|
+
"description": "Selects which mass and separation values to use when a supernova occurs during the CE and a kick needs to be applied.",
|
|
726
|
+
"type": "dropdown",
|
|
727
|
+
"options-preface": "",
|
|
728
|
+
"options": [
|
|
729
|
+
{
|
|
730
|
+
"name": 0,
|
|
731
|
+
"description": "uses pre-CE mass and post-CE sep (BSE default)"
|
|
732
|
+
},
|
|
733
|
+
{
|
|
734
|
+
"name": 1,
|
|
735
|
+
"description": "uses pre-CE mass and sep values"
|
|
736
|
+
},
|
|
737
|
+
{
|
|
738
|
+
"name": 2,
|
|
739
|
+
"description": "uses post-CE mass and sep",
|
|
740
|
+
"default": true
|
|
741
|
+
}
|
|
742
|
+
]
|
|
743
|
+
},
|
|
744
|
+
{
|
|
745
|
+
"name": "cemergeflag",
|
|
746
|
+
"description": "Determines whether stars that begin a CE without a distinct core-envelope boundary automatically lead to merger in a CE. These systems include: kstars = [0,1,2,7,8,10,11,12].",
|
|
747
|
+
"type": "dropdown",
|
|
748
|
+
"options-preface": "Note that while the optimal choice is <code>cemergeflag=1</code> according to <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>, <code>cemergeflag=0</code> allows for both options to be explored, since it is trivial to remove these systems from a population in post processing.",
|
|
749
|
+
"options": [
|
|
750
|
+
{
|
|
751
|
+
"name": 0,
|
|
752
|
+
"description": "allows the CE to proceed (optimistic CE)"
|
|
753
|
+
},
|
|
754
|
+
{
|
|
755
|
+
"name": 1,
|
|
756
|
+
"description": "forces the stars to merge in the CE (pessimistic CE)",
|
|
757
|
+
"default": true
|
|
758
|
+
}
|
|
759
|
+
]
|
|
760
|
+
},
|
|
761
|
+
{
|
|
762
|
+
"name": "cehestarflag",
|
|
763
|
+
"description": "Uses fitting formulae from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2015MNRAS.451.2123T/abstract'>Tauris+2015</a> for evolving RLO systems with a helium star donor and compact object accretor. NOTE: this flag will override <code>cekickflag</code> if set",
|
|
764
|
+
"type": "dropdown",
|
|
765
|
+
"options-preface": "",
|
|
766
|
+
"options": [
|
|
767
|
+
{
|
|
768
|
+
"name": 0,
|
|
769
|
+
"description": "Do not use the fitting formulae from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2015MNRAS.451.2123T/abstract'>Tauris+2015</a>",
|
|
770
|
+
"default": true
|
|
771
|
+
},
|
|
772
|
+
{
|
|
773
|
+
"name": 1,
|
|
774
|
+
"description": "Use the fitting formulae from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2015MNRAS.451.2123T/abstract'>Tauris+2015</a> for final period only"
|
|
775
|
+
},
|
|
776
|
+
{
|
|
777
|
+
"name": 2,
|
|
778
|
+
"description": "Use the fitting formulae from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2015MNRAS.451.2123T/abstract'>Tauris+2015</a> for both final mass and final period"
|
|
779
|
+
}
|
|
780
|
+
]
|
|
781
|
+
},
|
|
782
|
+
{
|
|
783
|
+
"name": "qcflag",
|
|
784
|
+
"description": "Selects model to determine critical mass ratios for the onset of unstable mass transfer and/or a common envelope during RLO. NOTE: this is overridden by <code>qcrit_array</code> if any of its values are non-zero.",
|
|
785
|
+
"type": "dropdown",
|
|
786
|
+
"options-preface": "The table below shows the values for <code>qcflag</code> across the choices of this flag and the kstar type. The equations in the table correspond to:<br> - Eq.1: <code class='docutils literal notranslate'><span class='pre'>qc = 0.362 + 1.0/(3.0*(1.0 - massc(j1)/mass(j1)))</span></code> (from Hjellming & Webbink 1983)<br> - Eq 2: <code class='docutils literal notranslate'><span class='pre'>qc = (1.67d0-zpars(7)+2.d0*(massc(j1)/mass(j1))**5)/2.13d0</span></code> (from Claeys+ 2014)",
|
|
787
|
+
"options": [
|
|
788
|
+
{
|
|
789
|
+
"name": 0,
|
|
790
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Section 2.6 of Hurley+2002</a> (Default BSE)"
|
|
791
|
+
},
|
|
792
|
+
{
|
|
793
|
+
"name": 1,
|
|
794
|
+
"description": "same as 0 but with <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/1987ApJ...318..794H/abstract'>Hjellming & Webbink 1987</a> for GB/AGB stars"
|
|
795
|
+
},
|
|
796
|
+
{
|
|
797
|
+
"name": 2,
|
|
798
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2014A%26A...563A..83C/abstract'>Table 2 of Claeys+2014</a>"
|
|
799
|
+
},
|
|
800
|
+
{
|
|
801
|
+
"name": 3,
|
|
802
|
+
"description": "same as 2 but with <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/1987ApJ...318..794H/abstract'>Hjellming & Webbink 1987</a> for GB/AGB stars"
|
|
803
|
+
},
|
|
804
|
+
{
|
|
805
|
+
"name": 4,
|
|
806
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Section 5.1 of Belcyznski+2008</a> except for WD donors which follow BSE"
|
|
807
|
+
},
|
|
808
|
+
{
|
|
809
|
+
"name": 5,
|
|
810
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2019MNRAS.490.3740N/abstract'>Section 2.3 of Neijssel+2020</a>; mass transfer from stripped stars is always assumed to be dynamically stable",
|
|
811
|
+
"default": true
|
|
812
|
+
}
|
|
813
|
+
]
|
|
814
|
+
},
|
|
815
|
+
{
|
|
816
|
+
"name": "qcrit_array",
|
|
817
|
+
"description": "User-defined values for critical mass ratios for the onset of unstable mass transfer and a common envelope during RLOF.",
|
|
818
|
+
"type": "string",
|
|
819
|
+
"options-preface": "Array of dimensions (1,16) specifying user-input values for the critical mass ratios that govern the onset of unstable mass transfer and a common envelope. Each item is set individually for its associated kstar, and a value of 0.0 will apply the prescription specified qcflag for that kstar.<br>NOTE: Recall from the start of this section that there are cases where a common envelope is forced regardless of the critical mass ratio for unstable mass transfer.",
|
|
820
|
+
"options": [
|
|
821
|
+
{
|
|
822
|
+
"name": "[0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0]",
|
|
823
|
+
"description": "Default value, which applies the prescription specified by <code>qcrit</code> for each kstar type",
|
|
824
|
+
"default": true
|
|
825
|
+
}
|
|
826
|
+
]
|
|
827
|
+
},
|
|
828
|
+
{
|
|
829
|
+
"settings-section": "Natal kicks",
|
|
830
|
+
"name": "kickflag",
|
|
831
|
+
"description": "Sets the particular natal kick prescription to use. Note that <code class='docutils literal notranslate'><span class='pre'>sigmadiv</span></code>, <code class='docutils literal notranslate'><span class='pre'>bhflag</span></code>, <code class='docutils literal notranslate'><span class='pre'>bhsigmafrac</span></code>, <code class='docutils literal notranslate'><span class='pre'>aic</span></code>, and <code class='docutils literal notranslate'><span class='pre'>ussn</span></code>, which are described below, are only used when <code class='docutils literal notranslate'><span class='pre'>abs(kickflag)=1</span></code>. Positive values use the Pfahl+2002 prescription for handling natal kicks.",
|
|
832
|
+
"type": "dropdown",
|
|
833
|
+
"options-preface": "",
|
|
834
|
+
"options": [
|
|
835
|
+
{
|
|
836
|
+
"name": 1,
|
|
837
|
+
"description": "The standard COSMIC kick prescription, where kicks are drawn from a bimodal distribution with standard FeCCSN getting a kick drawn from a Maxwellian distribution with dispersion parameter <code class='docutils literal notranslate'><span class='pre'>sigma</span></code> and ECSN/USSN are drawn according to <code class='docutils literal notranslate'><span class='pre'>sigmadiv</span></code>. This setting has additional possible options for <code class='docutils literal notranslate'><span class='pre'>bhflag</span></code>, <code class='docutils literal notranslate'><span class='pre'>bhsigmafrac</span></code>, <code class='docutils literal notranslate'><span class='pre'>aic</span></code> and <code class='docutils literal notranslate'><span class='pre'>ussn</span></code>.",
|
|
838
|
+
"default": true
|
|
839
|
+
},
|
|
840
|
+
{
|
|
841
|
+
"name": 2,
|
|
842
|
+
"description": "Natal kicks are drawn according to <code class='docutils literal notranslate'><span class='pre'>sigma</span></code> and scaled by the ejecta mass and remnant mass following Eq. 1 of <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2020ApJ...891..141G/abstract'>Giacobbo & Mapelli 2020</a></p> with their default parameters (\\(m_{\\rm NS = 1.2 {\\rm M_\\odot}\\), \\(m_{\\rm ej = 9 {\\rm M_\\odot}\\))"
|
|
843
|
+
},
|
|
844
|
+
{
|
|
845
|
+
"name": 3,
|
|
846
|
+
"description": "Natal kicks are drawn according to <code class='docutils literal notranslate'><span class='pre'>sigma</span></code> and scaled by just the ejecta mass following Eq. 2 of <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2020ApJ...891..141G/abstract'>Giacobbo & Mapelli 2020</a>, which does not scale the kick by (\\(m_{\\rm NS\\)"
|
|
847
|
+
},
|
|
848
|
+
{
|
|
849
|
+
"name": 4,
|
|
850
|
+
"description": "Natal kicks are drawn according to Eq. 1 of <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2016MNRAS.461.3747B/abstract'>Bray & Eldridge 2016</a>, with their default parameters (\\(\\alpha=70 \\, {\\rm km/s}, \\beta = 120 \\, {\\rm km/s)}"
|
|
851
|
+
},
|
|
852
|
+
{
|
|
853
|
+
"name": 5,
|
|
854
|
+
"description": "Follows the same prescription as 1, but uses the kick prescription described in <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2025arXiv250522102D/abstract'>Disberg & Mandel 2025</a> for CCSN."
|
|
855
|
+
},
|
|
856
|
+
{
|
|
857
|
+
"name": -1,
|
|
858
|
+
"description": "As 1, but using the old Kiel & Hurley 2009 prescription for changing the orbital configuration of the binary, available for reproducibility purposes but not recommended for new work"
|
|
859
|
+
},
|
|
860
|
+
{
|
|
861
|
+
"name": -2,
|
|
862
|
+
"description": "As 2, but using the old Kiel & Hurley 2009 prescription for changing the orbital configuration of the binary, available for reproducibility purposes but not recommended for new work"
|
|
863
|
+
},
|
|
864
|
+
{
|
|
865
|
+
"name": -3,
|
|
866
|
+
"description": "As 3, but using the old Kiel & Hurley 2009 prescription for changing the orbital configuration of the binary, available for reproducibility purposes but not recommended for new work"
|
|
867
|
+
},
|
|
868
|
+
{
|
|
869
|
+
"name": -4,
|
|
870
|
+
"description": "As 4, but using the old Kiel & Hurley 2009 prescription for changing the orbital configuration of the binary, available for reproducibility purposes but not recommended for new work"
|
|
871
|
+
}
|
|
872
|
+
]
|
|
873
|
+
},
|
|
874
|
+
{
|
|
875
|
+
"name": "sigma",
|
|
876
|
+
"description": "Sets the dispersion in the Maxwellian for the SN kick velocity in km/s",
|
|
877
|
+
"type": "number",
|
|
878
|
+
"options-preface": "",
|
|
879
|
+
"options": [
|
|
880
|
+
{
|
|
881
|
+
"name": "positive values",
|
|
882
|
+
"description": "Sets the dispersion in the Maxwellian for the SN kick velocity in km/s"
|
|
883
|
+
},
|
|
884
|
+
{
|
|
885
|
+
"name": 0.0,
|
|
886
|
+
"description": "No natal kicks for core-collapse when using Hobbs kicks"
|
|
887
|
+
},
|
|
888
|
+
{
|
|
889
|
+
"name": 265.0,
|
|
890
|
+
"description": "Default choice",
|
|
891
|
+
"default": true
|
|
892
|
+
}
|
|
893
|
+
]
|
|
894
|
+
},
|
|
895
|
+
{
|
|
896
|
+
"name": "bhflag",
|
|
897
|
+
"description": "Sets the model for how SN kicks are applied to BHs, where bhflag != 0 allows for velocity kick at BH formation",
|
|
898
|
+
"type": "dropdown",
|
|
899
|
+
"options-preface": "",
|
|
900
|
+
"options": [
|
|
901
|
+
{
|
|
902
|
+
"name": 0,
|
|
903
|
+
"description": "No BH kick"
|
|
904
|
+
},
|
|
905
|
+
{
|
|
906
|
+
"name": 1,
|
|
907
|
+
"description": "fallback-modulated kicks following <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2012ApJ...749...91F/abstract'>Fryer+2012</a>",
|
|
908
|
+
"default": true
|
|
909
|
+
},
|
|
910
|
+
{
|
|
911
|
+
"name": 2,
|
|
912
|
+
"description": "kicks decreased by ratio of BH mass to NS mass (1.44 Msun); conserves linear momentum"
|
|
913
|
+
},
|
|
914
|
+
{
|
|
915
|
+
"name": 3,
|
|
916
|
+
"description": "BH natal kicks are not decreased compared to NS kicks and are drawn from the same Maxwellian distribution with dispersion = <code>sigma</code> set above"
|
|
917
|
+
}
|
|
918
|
+
]
|
|
919
|
+
},
|
|
920
|
+
{
|
|
921
|
+
"name": "bhsigmafrac",
|
|
922
|
+
"description": "Sets a fractional modification which scales down sigma for BHs. This works in addition to whatever is chosen for bhflag, and is applied to sigma before the bhflag prescriptions are applied",
|
|
923
|
+
"type": "number",
|
|
924
|
+
"options-preface": "",
|
|
925
|
+
"options": [
|
|
926
|
+
{
|
|
927
|
+
"name": "range (0, 1]",
|
|
928
|
+
"description": "reduces sigma by bhsigmafrac for BHs"
|
|
929
|
+
},
|
|
930
|
+
{
|
|
931
|
+
"name": 1.0,
|
|
932
|
+
"description": "Default choice",
|
|
933
|
+
"default": true
|
|
934
|
+
}
|
|
935
|
+
]
|
|
936
|
+
},
|
|
937
|
+
{
|
|
938
|
+
"name": "sigmadiv",
|
|
939
|
+
"description": "Sets the modified ECSN kick strength",
|
|
940
|
+
"type": "number",
|
|
941
|
+
"options-preface": "",
|
|
942
|
+
"options": [
|
|
943
|
+
{
|
|
944
|
+
"name": "positive values",
|
|
945
|
+
"description": "divide sigma (defined above) by sigmadiv"
|
|
946
|
+
},
|
|
947
|
+
{
|
|
948
|
+
"name": "negative values",
|
|
949
|
+
"description": " sets ECSN kicks to be drawn from a Maxwellian distribution with dispersion given by sigmadiv"
|
|
950
|
+
},
|
|
951
|
+
{
|
|
952
|
+
"name": -20.0,
|
|
953
|
+
"description": "Default choice",
|
|
954
|
+
"default": true
|
|
955
|
+
}
|
|
956
|
+
]
|
|
957
|
+
},
|
|
958
|
+
{
|
|
959
|
+
"name": "ecsn",
|
|
960
|
+
"description": "Allows for electron capture SNe and sets the maximum He-star mass (at core helium depletion) that will result in an ECSN",
|
|
961
|
+
"type": "number",
|
|
962
|
+
"options-preface": "",
|
|
963
|
+
"options": [
|
|
964
|
+
{
|
|
965
|
+
"name": 0,
|
|
966
|
+
"description": "turns off ECSN"
|
|
967
|
+
},
|
|
968
|
+
{
|
|
969
|
+
"name": "positive values",
|
|
970
|
+
"description": "sets maximum He-star mass for ECSN; <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>BSE (Hurley+2002)</a> and <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>StarTrack (Belczynski+2008)</a> use ecsn = 2.25, while <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2004ApJ...612.1044P/abstract'>Podsiadlowksi+2004</a> argues that binarity can increase this to ecsn = 2.5"
|
|
971
|
+
},
|
|
972
|
+
{
|
|
973
|
+
"name": 2.25,
|
|
974
|
+
"description": "Default choice",
|
|
975
|
+
"default": true
|
|
976
|
+
}
|
|
977
|
+
]
|
|
978
|
+
},
|
|
979
|
+
{
|
|
980
|
+
"name": "ecsn_mlow",
|
|
981
|
+
"description": "Sets the low end of the ECSN mass range",
|
|
982
|
+
"type": "number",
|
|
983
|
+
"options-preface": "",
|
|
984
|
+
"options": [
|
|
985
|
+
{
|
|
986
|
+
"name": "positive values",
|
|
987
|
+
"description": "sets maximum He-star mass for ECSN; <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>BSE (Hurley+2002)</a> uses ecsn_mlow = 1.6, <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>StarTrack (Belczynski+2008)</a> uses ecsn_mlow = 1.85, <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2004ApJ...612.1044P/abstract'>Podsiadlowksi+2004</a> argues that binarity can decrease this to ecsn_mlow = 1.4"
|
|
988
|
+
},
|
|
989
|
+
{
|
|
990
|
+
"name": 0.0,
|
|
991
|
+
"description": "no lower limit on He-star mass for ECSN"
|
|
992
|
+
},
|
|
993
|
+
{
|
|
994
|
+
"name": 1.6,
|
|
995
|
+
"description": "Default choice",
|
|
996
|
+
"default": true
|
|
997
|
+
}
|
|
998
|
+
]
|
|
999
|
+
},
|
|
1000
|
+
{
|
|
1001
|
+
"name": "aic",
|
|
1002
|
+
"description": "Sets the model for accretion induced collapse SN natal kicks",
|
|
1003
|
+
"type": "dropdown",
|
|
1004
|
+
"options-preface": "Applies even if <code>ecsn = 0<.code>",
|
|
1005
|
+
"options": [
|
|
1006
|
+
{
|
|
1007
|
+
"name": 0,
|
|
1008
|
+
"description": "AIC SN receive kicks drawn from Maxwellian with dispersion = sigma defined above"
|
|
1009
|
+
},
|
|
1010
|
+
{
|
|
1011
|
+
"name": 1,
|
|
1012
|
+
"description": "sets AIC SN kick strength according to sigmadiv; NOTE that this will apply even if ecsn = 0.0",
|
|
1013
|
+
"default": true
|
|
1014
|
+
}
|
|
1015
|
+
]
|
|
1016
|
+
},
|
|
1017
|
+
{
|
|
1018
|
+
"name": "ussn",
|
|
1019
|
+
"description": "Reduces kicks according to the sigmadiv selection for ultra-stripped supernovae, assumed to happen if a He-star undergoes a CE with a compact companion",
|
|
1020
|
+
"type": "dropdown",
|
|
1021
|
+
"options-preface": "",
|
|
1022
|
+
"options": [
|
|
1023
|
+
{
|
|
1024
|
+
"name": 0,
|
|
1025
|
+
"description": "USSN receive kicks drawn from Maxwellian with dispersion = sigma defined above"
|
|
1026
|
+
},
|
|
1027
|
+
{
|
|
1028
|
+
"name": 1,
|
|
1029
|
+
"description": "sets USSN kick strength according to sigmadiv",
|
|
1030
|
+
"default": true
|
|
1031
|
+
}
|
|
1032
|
+
]
|
|
1033
|
+
},
|
|
1034
|
+
{
|
|
1035
|
+
"name": "pisn",
|
|
1036
|
+
"description": "Allows for (pulsational) pair instability supernovae and sets either the model to use or the maximum mass of the remnant.",
|
|
1037
|
+
"type": "number",
|
|
1038
|
+
"options-preface": "",
|
|
1039
|
+
"options": [
|
|
1040
|
+
{
|
|
1041
|
+
"name": 0,
|
|
1042
|
+
"description": "no pulsational pair instability SN"
|
|
1043
|
+
},
|
|
1044
|
+
{
|
|
1045
|
+
"name": -1,
|
|
1046
|
+
"description": " uses the formulae from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2017MNRAS.470.4739S/abstract'>Spera & Mapelli 2017</a>"
|
|
1047
|
+
},
|
|
1048
|
+
{
|
|
1049
|
+
"name": -2,
|
|
1050
|
+
"description": "uses a polynomial fit to <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2018arXiv181013412M/abstract'>Table 1 in Marchant+2018</a>",
|
|
1051
|
+
"default": true
|
|
1052
|
+
},
|
|
1053
|
+
{
|
|
1054
|
+
"name": -3,
|
|
1055
|
+
"description": "uses a polynomial fit to <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2019ApJ...878...49W/abstract'>Table 5 in Woosley 2019</a>"
|
|
1056
|
+
},
|
|
1057
|
+
{
|
|
1058
|
+
"name": "positive values",
|
|
1059
|
+
"description": "turns on pulsational pair instability and pair instability SNe, and sets the maximum mass of the allowed remnant (i.e., the bottom of the pair instability mass gap). He core masses between pisn and 65 Msun are assumed to go through pulsational pair instability and limit the He core mass to pisn, while He core masses from 65-135 Msun are assumed have a pair instability SN and leave no remnant."
|
|
1060
|
+
}
|
|
1061
|
+
]
|
|
1062
|
+
},
|
|
1063
|
+
{
|
|
1064
|
+
"name": "polar_kick_angle",
|
|
1065
|
+
"description": "Sets the opening angle of the SN kick relative to the pole of the exploding star",
|
|
1066
|
+
"type": "number",
|
|
1067
|
+
"options-preface": "",
|
|
1068
|
+
"options": [
|
|
1069
|
+
{
|
|
1070
|
+
"name": "range [0, 90]",
|
|
1071
|
+
"description": "Sets the opening angle of the SN kick relative to the pole of the exploding star"
|
|
1072
|
+
},
|
|
1073
|
+
{
|
|
1074
|
+
"name": 0.0,
|
|
1075
|
+
"description": "Strictly polar kicks"
|
|
1076
|
+
},
|
|
1077
|
+
{
|
|
1078
|
+
"name": 90.0,
|
|
1079
|
+
"description": "Fully isotropic kicks (default choice)",
|
|
1080
|
+
"default": true
|
|
1081
|
+
}
|
|
1082
|
+
]
|
|
1083
|
+
},
|
|
1084
|
+
{
|
|
1085
|
+
"name": "natal_kick_array",
|
|
1086
|
+
"description": "Array of dimensions (2,5) which takes user input values for the SN natal kick, where the first row corresponds to the first star and the second row corresponds to the second star and columns are: [vk, phi, theta, mean_anomaly, rand_seed].",
|
|
1087
|
+
"type": "string",
|
|
1088
|
+
"options-preface": "NOTE: any numbers outside the ranges below will be sampled in the standard ways detailed above.",
|
|
1089
|
+
"options": [
|
|
1090
|
+
{
|
|
1091
|
+
"name": "vk",
|
|
1092
|
+
"description": "SN kick velocity in km/s, valid on the range [0, inf]"
|
|
1093
|
+
},
|
|
1094
|
+
{
|
|
1095
|
+
"name": "phi",
|
|
1096
|
+
"description": "SN kick co-lateral polar angle in degrees, valid on the range [-90, 90]"
|
|
1097
|
+
},
|
|
1098
|
+
{
|
|
1099
|
+
"name": "theta",
|
|
1100
|
+
"description": "SN kick azimuthal angle in degrees, valid on the range [0, 360]"
|
|
1101
|
+
},
|
|
1102
|
+
{
|
|
1103
|
+
"name": "mean_anomaly",
|
|
1104
|
+
"description": "SN kick mean anomaly in degrees, valid on the range [0, 360]"
|
|
1105
|
+
},
|
|
1106
|
+
{
|
|
1107
|
+
"name": "rand_seed",
|
|
1108
|
+
"description": "supplied if restarting evolution after a supernova has already occurred"
|
|
1109
|
+
},
|
|
1110
|
+
{
|
|
1111
|
+
"name": "[[-100.0, -100.0, -100.0, -100.0, 0.0], [-100.0, -100.0, -100.0, -100.0, 0.0]]",
|
|
1112
|
+
"description": "Default choice, which applies the standard natal kick prescription",
|
|
1113
|
+
"default": true
|
|
1114
|
+
}
|
|
1115
|
+
]
|
|
1116
|
+
},
|
|
1117
|
+
{
|
|
1118
|
+
"settings-section": "Remnant mass",
|
|
1119
|
+
"name": "remnantflag",
|
|
1120
|
+
"description": "Determines the remnant mass prescription used for NSs and BHs.",
|
|
1121
|
+
"type": "dropdown",
|
|
1122
|
+
"options-preface": "",
|
|
1123
|
+
"options": [
|
|
1124
|
+
{
|
|
1125
|
+
"name": 0,
|
|
1126
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2000MNRAS.315..543H/abstract'>Section 6 of Hurley+2000</a> (default BSE)"
|
|
1127
|
+
},
|
|
1128
|
+
{
|
|
1129
|
+
"name": 1,
|
|
1130
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002ApJ...572..407B/abstract'>Belczynski+2002</a>"
|
|
1131
|
+
},
|
|
1132
|
+
{
|
|
1133
|
+
"name": 2,
|
|
1134
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>"
|
|
1135
|
+
},
|
|
1136
|
+
{
|
|
1137
|
+
"name": 3,
|
|
1138
|
+
"description": "follows the rapid prescription from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2012ApJ...749...91F/abstract'>Fryer+2012</a>, with updated proto-core mass from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2020ApJ...891..141G/abstract'>Giacobbo & Mapelli 2020</a>. This leads to a mass gap between neutron stars and black holes."
|
|
1139
|
+
},
|
|
1140
|
+
{
|
|
1141
|
+
"name": 4,
|
|
1142
|
+
"description": "follows the delayed prescription from <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2012ApJ...749...91F/abstract'>Fryer+2012</a>. This fills the mass gap between neutron stars and black holes.",
|
|
1143
|
+
"default": true
|
|
1144
|
+
}
|
|
1145
|
+
]
|
|
1146
|
+
},
|
|
1147
|
+
{
|
|
1148
|
+
"name": "mxns",
|
|
1149
|
+
"description": "Sets the boundary between the maximum NS mass and the minimum BH mass",
|
|
1150
|
+
"type": "number",
|
|
1151
|
+
"options-preface": "",
|
|
1152
|
+
"options": [
|
|
1153
|
+
{
|
|
1154
|
+
"name": "positive values",
|
|
1155
|
+
"description": "Sets the boundary between the maximum NS mass and the minimum BH mass"
|
|
1156
|
+
},
|
|
1157
|
+
{
|
|
1158
|
+
"name": 3.0,
|
|
1159
|
+
"description": "Default choice",
|
|
1160
|
+
"default": true
|
|
1161
|
+
}
|
|
1162
|
+
]
|
|
1163
|
+
},
|
|
1164
|
+
{
|
|
1165
|
+
"name": "rembar_massloss",
|
|
1166
|
+
"description": "Determines the prescriptions for mass conversion due to neutrino emission during the collapse of the proto-compact object",
|
|
1167
|
+
"type": "number",
|
|
1168
|
+
"options-preface": "",
|
|
1169
|
+
"options": [
|
|
1170
|
+
{
|
|
1171
|
+
"name": "positive values",
|
|
1172
|
+
"description": "sets the maximum amount of mass loss, which should be about 10% of the maximum mass of an iron core (\\({\\sim 5 \\mathrm{M}_\\odot}\\) Fryer, private communication)"
|
|
1173
|
+
},
|
|
1174
|
+
{
|
|
1175
|
+
"name": "range [-1, 0]",
|
|
1176
|
+
"description": "assumes that proto-compact objects lose a constant fraction of their baryonic mass when collapsing to a black hole, such that \\(M_{\\rm rem} = (1 + \\texttt{rembar\\_massloss}) M_{\\rm rem}\\) (e.g., rembar_massloss = -0.1 gives the black hole a gravitational mass that is 90% of the proto-compact object's baryonic mass)"
|
|
1177
|
+
},
|
|
1178
|
+
{
|
|
1179
|
+
"name": 0.5,
|
|
1180
|
+
"description": "Default choice",
|
|
1181
|
+
"default": true
|
|
1182
|
+
}
|
|
1183
|
+
]
|
|
1184
|
+
},
|
|
1185
|
+
{
|
|
1186
|
+
"name": "wd_mass_lim",
|
|
1187
|
+
"description": "Determines if the maximum white dwarf mass is limited to the chandraekhar mass during merger induced collapse",
|
|
1188
|
+
"type": "dropdown",
|
|
1189
|
+
"options-preface": "",
|
|
1190
|
+
"options": [
|
|
1191
|
+
{
|
|
1192
|
+
"name": 0,
|
|
1193
|
+
"description": "Do not apply the limit"
|
|
1194
|
+
},
|
|
1195
|
+
{
|
|
1196
|
+
"name": 1,
|
|
1197
|
+
"description": "Apply the limit",
|
|
1198
|
+
"default": true
|
|
1199
|
+
}
|
|
1200
|
+
]
|
|
1201
|
+
},
|
|
1202
|
+
{
|
|
1203
|
+
"settings-section": "Remnant spin",
|
|
1204
|
+
"name": "bhspinflag",
|
|
1205
|
+
"description": "Uses different prescriptions for BH spin after formation",
|
|
1206
|
+
"type": "dropdown",
|
|
1207
|
+
"options-preface": "",
|
|
1208
|
+
"options": [
|
|
1209
|
+
{
|
|
1210
|
+
"name": 0,
|
|
1211
|
+
"description": "sets all BH spins to bhspinmag",
|
|
1212
|
+
"default": true
|
|
1213
|
+
},
|
|
1214
|
+
{
|
|
1215
|
+
"name": 1,
|
|
1216
|
+
"description": "draws a random BH spin between 0 and bhspinmag for every BH"
|
|
1217
|
+
},
|
|
1218
|
+
{
|
|
1219
|
+
"name": 2,
|
|
1220
|
+
"description": "core-mass dependent BH spin (based on <a class='reference external' href='https://arxiv.org/abs/1706.07053v1'>Belczynski+2017 v1</a>)"
|
|
1221
|
+
}
|
|
1222
|
+
]
|
|
1223
|
+
},
|
|
1224
|
+
{
|
|
1225
|
+
"name": "bhspinmag",
|
|
1226
|
+
"description": "Sets either the spin of all BHs or the upper limit of the uniform distribution for BH spins (see bhspinflag)",
|
|
1227
|
+
"type": "number",
|
|
1228
|
+
"options-preface": "",
|
|
1229
|
+
"options": [
|
|
1230
|
+
{
|
|
1231
|
+
"name": "range [0, 1]",
|
|
1232
|
+
"description": "Sets either the spin of all BHs or the upper limit of the uniform distribution for BH spins (see bhspinflag)"
|
|
1233
|
+
},
|
|
1234
|
+
{
|
|
1235
|
+
"name": 0.0,
|
|
1236
|
+
"description": "Default choice",
|
|
1237
|
+
"default": true
|
|
1238
|
+
}
|
|
1239
|
+
]
|
|
1240
|
+
},
|
|
1241
|
+
{
|
|
1242
|
+
"settings-section": "GR Orbital Decay",
|
|
1243
|
+
"name": "grflag",
|
|
1244
|
+
"description": "Turns on or off orbital decay due to gravitational wave emission",
|
|
1245
|
+
"type": "dropdown",
|
|
1246
|
+
"options-preface": "",
|
|
1247
|
+
"options": [
|
|
1248
|
+
{
|
|
1249
|
+
"name": 0,
|
|
1250
|
+
"description": "No orbital decay due to gravitational wave emission"
|
|
1251
|
+
},
|
|
1252
|
+
{
|
|
1253
|
+
"name": 1,
|
|
1254
|
+
"description": "Orbital decay due to gravitational wave emission is turned on",
|
|
1255
|
+
"default": true
|
|
1256
|
+
}
|
|
1257
|
+
]
|
|
1258
|
+
},
|
|
1259
|
+
{
|
|
1260
|
+
"settings-section": "Mass transfer",
|
|
1261
|
+
"name": "eddfac",
|
|
1262
|
+
"description": "Eddington limit factor for mass transfer.",
|
|
1263
|
+
"options-preface": "",
|
|
1264
|
+
"type": "number",
|
|
1265
|
+
"options": [
|
|
1266
|
+
{
|
|
1267
|
+
"name": "positive values",
|
|
1268
|
+
"description": "restrict accretion limit to fraction of Eddington (sub- or super-Eddington accretion)"
|
|
1269
|
+
},
|
|
1270
|
+
{
|
|
1271
|
+
"name": 0,
|
|
1272
|
+
"description": "no mass transfer onto compact objects"
|
|
1273
|
+
},
|
|
1274
|
+
{
|
|
1275
|
+
"name": 1,
|
|
1276
|
+
"description": "mass transfer rate is limited by the Eddington rate following Equation 67 in <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a>"
|
|
1277
|
+
},
|
|
1278
|
+
{
|
|
1279
|
+
"name": 10,
|
|
1280
|
+
"description": "mass transfer rate is limited to 10x the Eddington rate",
|
|
1281
|
+
"default": true
|
|
1282
|
+
}
|
|
1283
|
+
]
|
|
1284
|
+
},
|
|
1285
|
+
{
|
|
1286
|
+
"name": "gamma",
|
|
1287
|
+
"description": "Angular momentum prescriptions for mass lost during Roche-lobe overflow at super-Eddington mass transfer rates",
|
|
1288
|
+
"type": "number",
|
|
1289
|
+
"options-preface": "",
|
|
1290
|
+
"options": [
|
|
1291
|
+
{
|
|
1292
|
+
"name": -1,
|
|
1293
|
+
"description": "assumes the lost material carries away the specific angular momentum of the primary"
|
|
1294
|
+
},
|
|
1295
|
+
{
|
|
1296
|
+
"name": -2,
|
|
1297
|
+
"description": "assumes material is lost from the system as if it is a wind from the secondary",
|
|
1298
|
+
"default": true
|
|
1299
|
+
},
|
|
1300
|
+
{
|
|
1301
|
+
"name": -3,
|
|
1302
|
+
"description": "assumes mass is lost through the outer Lagrangian point, forming a circumbinary disk. See Zapartas+17 Eq. 9 and Artymowicz & Lubow (1994)."
|
|
1303
|
+
},
|
|
1304
|
+
{
|
|
1305
|
+
"name": "positive values",
|
|
1306
|
+
"description": "assumes that the lost material takes away a fraction <code>gamma</code> of the orbital angular momentum"
|
|
1307
|
+
},
|
|
1308
|
+
{
|
|
1309
|
+
"name": 0.0,
|
|
1310
|
+
"description": "no angular momentum loss from the system due to mass loss"
|
|
1311
|
+
}
|
|
1312
|
+
]
|
|
1313
|
+
},
|
|
1314
|
+
{
|
|
1315
|
+
"name": "don_lim",
|
|
1316
|
+
"description": "Determines the rate of mass loss through Roche-lobe overflow mass transfer from the donor star",
|
|
1317
|
+
"type": "dropdown",
|
|
1318
|
+
"options-preface": "",
|
|
1319
|
+
"options": [
|
|
1320
|
+
{
|
|
1321
|
+
"name": -1,
|
|
1322
|
+
"description": "donor mass loss rate is calculated following <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a>",
|
|
1323
|
+
"default": true
|
|
1324
|
+
},
|
|
1325
|
+
{
|
|
1326
|
+
"name": -2,
|
|
1327
|
+
"description": "donor mass loss rate is calculated following <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2014A%26A...563A..83C/abstract'>Claeys+2014</a>"
|
|
1328
|
+
}
|
|
1329
|
+
]
|
|
1330
|
+
},
|
|
1331
|
+
{
|
|
1332
|
+
"name": "acc_lim",
|
|
1333
|
+
"description": "Limits the amount of mass accreted during Roche-lobe overflow",
|
|
1334
|
+
"type": "number",
|
|
1335
|
+
"options-preface": "",
|
|
1336
|
+
"options": [
|
|
1337
|
+
{
|
|
1338
|
+
"name": -1,
|
|
1339
|
+
"description": "limited to 10x the thermal rate of the accretor for MS/HG/CHeB and unlimited for GB/EAGB/AGB stars",
|
|
1340
|
+
"default": true
|
|
1341
|
+
},
|
|
1342
|
+
{
|
|
1343
|
+
"name": -2,
|
|
1344
|
+
"description": "limited to 1x the thermal rate of the accretor for MS/HG/CHeB and unlimited for GB/EAGB/AGB stars"
|
|
1345
|
+
},
|
|
1346
|
+
{
|
|
1347
|
+
"name": -3,
|
|
1348
|
+
"description": "limited to 10x the thermal rate of the accretor for all stars"
|
|
1349
|
+
},
|
|
1350
|
+
{
|
|
1351
|
+
"name": -4,
|
|
1352
|
+
"description": "limited to 1x the thermal rate of the accretor for all stars"
|
|
1353
|
+
},
|
|
1354
|
+
{
|
|
1355
|
+
"name": "positive values",
|
|
1356
|
+
"description": "sets overall fraction of donor material that is accreted, with the rest being lost from the system (<code>acc_lim = 0.5</code> assumes 50% accretion efficiency as in <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>)"
|
|
1357
|
+
},
|
|
1358
|
+
{
|
|
1359
|
+
"name": 0.0,
|
|
1360
|
+
"description": "assume all mass is lost from the system (fully nonconservative)"
|
|
1361
|
+
}
|
|
1362
|
+
]
|
|
1363
|
+
},
|
|
1364
|
+
{
|
|
1365
|
+
"settings-section": "Tides",
|
|
1366
|
+
"name": "tflag",
|
|
1367
|
+
"description": "Activates tidal circularization following <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a>",
|
|
1368
|
+
"type": "dropdown",
|
|
1369
|
+
"options-preface": "",
|
|
1370
|
+
"options": [
|
|
1371
|
+
{
|
|
1372
|
+
"name": 0,
|
|
1373
|
+
"description": "No tidal circularization"
|
|
1374
|
+
},
|
|
1375
|
+
{
|
|
1376
|
+
"name": 1,
|
|
1377
|
+
"description": "Tidal circularization is turned on",
|
|
1378
|
+
"default": true
|
|
1379
|
+
}
|
|
1380
|
+
]
|
|
1381
|
+
},
|
|
1382
|
+
{
|
|
1383
|
+
"name": "ST_tide",
|
|
1384
|
+
"description": "Activates StarTrack setup for tides following <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>",
|
|
1385
|
+
"type": "dropdown",
|
|
1386
|
+
"options-preface": "",
|
|
1387
|
+
"options": [
|
|
1388
|
+
{
|
|
1389
|
+
"name": 0,
|
|
1390
|
+
"description": "Follows BSE (<a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a>)"
|
|
1391
|
+
},
|
|
1392
|
+
{
|
|
1393
|
+
"name": 1,
|
|
1394
|
+
"description": "Follow StarTrack (<a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>). Note StarTrack method does not use a better integration scheme (yet) but simply ; follows similar set up to startrack (including initial vrot, using roche-lobe check ; at periastron, and circularisation and synchronisation at start of MT).",
|
|
1395
|
+
"default": true
|
|
1396
|
+
}
|
|
1397
|
+
]
|
|
1398
|
+
},
|
|
1399
|
+
{
|
|
1400
|
+
"name": "fprimc_array",
|
|
1401
|
+
"description": "Controls the scaling factor for convective tides. Each value in the array is set individually for its associated kstar. The relevant equation is <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002 Eq. 21</a>.",
|
|
1402
|
+
"type": "string",
|
|
1403
|
+
"options-preface": "",
|
|
1404
|
+
"options": [
|
|
1405
|
+
{
|
|
1406
|
+
"name": "positive values",
|
|
1407
|
+
"description": "sets scaling factor of Equation 21 referenced above"
|
|
1408
|
+
},
|
|
1409
|
+
{
|
|
1410
|
+
"name": "[2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0,2.0/21.0]",
|
|
1411
|
+
"description": "Default choice, which applies the standard scaling factor",
|
|
1412
|
+
"default": true
|
|
1413
|
+
}
|
|
1414
|
+
]
|
|
1415
|
+
},
|
|
1416
|
+
{
|
|
1417
|
+
"settings-section": "White dwarfs",
|
|
1418
|
+
"name": "ifflag",
|
|
1419
|
+
"description": "Activates the initial-final white dwarf mass relation from Han+1995 <a class='reference external' href='http://adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1995MNRAS.272..800H&link_type=ARTICLE&db_key=AST&high=#page=4'>Equations 3, 4, and 5</a>.",
|
|
1420
|
+
"type": "dropdown",
|
|
1421
|
+
"options-preface": "",
|
|
1422
|
+
"options": [
|
|
1423
|
+
{
|
|
1424
|
+
"name": 0,
|
|
1425
|
+
"description": "No modifications to BSE"
|
|
1426
|
+
},
|
|
1427
|
+
{
|
|
1428
|
+
"name": 1,
|
|
1429
|
+
"description": "Initial-final white dwarf mass relation is turned on",
|
|
1430
|
+
"default": true
|
|
1431
|
+
}
|
|
1432
|
+
]
|
|
1433
|
+
},
|
|
1434
|
+
{
|
|
1435
|
+
"name": "wdflag",
|
|
1436
|
+
"description": "Activates an alternate cooling law found in the description immediately following <a class='reference external' href='http://iopscience.iop.org/article/10.1086/374637/pdf#page=3'>Equation 1</a> in Hurley & Shara 2003. Equation 1 gives the BSE default Mestel cooling law.",
|
|
1437
|
+
"type": "dropdown",
|
|
1438
|
+
"options-preface": "",
|
|
1439
|
+
"options": [
|
|
1440
|
+
{
|
|
1441
|
+
"name": 0,
|
|
1442
|
+
"description": "No modifications to BSE"
|
|
1443
|
+
},
|
|
1444
|
+
{
|
|
1445
|
+
"name": 1,
|
|
1446
|
+
"description": "Alternate cooling law is turned on",
|
|
1447
|
+
"default": true
|
|
1448
|
+
}
|
|
1449
|
+
]
|
|
1450
|
+
},
|
|
1451
|
+
{
|
|
1452
|
+
"name": "epsnov",
|
|
1453
|
+
"description": "Fraction of accreted matter retained in a nova eruption. This is relevant for accretion onto degenerate objects; see Section 2.6.6.2 in <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a>.",
|
|
1454
|
+
"type": "number",
|
|
1455
|
+
"options-preface": "",
|
|
1456
|
+
"options": [
|
|
1457
|
+
{
|
|
1458
|
+
"name": "range [0, 1]",
|
|
1459
|
+
"description": "Retains epsnov fraction of accreted matter"
|
|
1460
|
+
},
|
|
1461
|
+
{
|
|
1462
|
+
"name": 0.001,
|
|
1463
|
+
"description": "Default choice",
|
|
1464
|
+
"default": true
|
|
1465
|
+
}
|
|
1466
|
+
]
|
|
1467
|
+
},
|
|
1468
|
+
{
|
|
1469
|
+
"settings-section": "Pulsars",
|
|
1470
|
+
"name": "bdecayfac",
|
|
1471
|
+
"description": "Activates different models for accretion induced field decay; see <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008MNRAS.388..393K/abstract'>Kiel+2008</a>.",
|
|
1472
|
+
"type": "dropdown",
|
|
1473
|
+
"options-preface": "",
|
|
1474
|
+
"options": [
|
|
1475
|
+
{
|
|
1476
|
+
"name": 0,
|
|
1477
|
+
"description": "uses an exponential decay"
|
|
1478
|
+
},
|
|
1479
|
+
{
|
|
1480
|
+
"name": 1,
|
|
1481
|
+
"description": "uses an inverse decay",
|
|
1482
|
+
"default": true
|
|
1483
|
+
}
|
|
1484
|
+
]
|
|
1485
|
+
},
|
|
1486
|
+
{
|
|
1487
|
+
"name": "bconst",
|
|
1488
|
+
"description": "Sets the magnetic field decay timescale for pulsars following Section 3 of <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008MNRAS.388..393K/abstract'>Kiel+2008</a>.",
|
|
1489
|
+
"type": "number",
|
|
1490
|
+
"options-preface": "",
|
|
1491
|
+
"options": [
|
|
1492
|
+
{
|
|
1493
|
+
"name": "positive values",
|
|
1494
|
+
"description": "sets \\(k\\) in Myr from Equation 8 to <code>bconst</code>"
|
|
1495
|
+
},
|
|
1496
|
+
{
|
|
1497
|
+
"name": 3000,
|
|
1498
|
+
"description": "Default choice",
|
|
1499
|
+
"default": true
|
|
1500
|
+
}
|
|
1501
|
+
]
|
|
1502
|
+
},
|
|
1503
|
+
{
|
|
1504
|
+
"name": "ck",
|
|
1505
|
+
"description": "Sets the magnetic field decay timescale for pulsars following Section 3 of <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008MNRAS.388..393K/abstract'>Kiel+2008</a>.",
|
|
1506
|
+
"type": "number",
|
|
1507
|
+
"options-preface": "",
|
|
1508
|
+
"options": [
|
|
1509
|
+
{
|
|
1510
|
+
"name": "positive values",
|
|
1511
|
+
"description": "sets \\(\\tau_b\\) in Myr from Equation 2 to <code>ck</code>"
|
|
1512
|
+
},
|
|
1513
|
+
{
|
|
1514
|
+
"name": 1000,
|
|
1515
|
+
"description": "Default choice",
|
|
1516
|
+
"default": true
|
|
1517
|
+
}
|
|
1518
|
+
]
|
|
1519
|
+
},
|
|
1520
|
+
{
|
|
1521
|
+
"settings-section": "Mixing variables",
|
|
1522
|
+
"name": "rejuv_fac",
|
|
1523
|
+
"description": "Sets the mixing factor in main sequence star collisions. This is hard coded to 0.1 in the original BSE release and in Equation 80 of <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a> but can lead to extended main sequence lifetimes in some cases.",
|
|
1524
|
+
"type": "number",
|
|
1525
|
+
"options-preface": "",
|
|
1526
|
+
"options": [
|
|
1527
|
+
{
|
|
1528
|
+
"name": "range [0, 1]",
|
|
1529
|
+
"description": "sets the mixing factor in main sequence star collisions"
|
|
1530
|
+
},
|
|
1531
|
+
{
|
|
1532
|
+
"name": 1.0,
|
|
1533
|
+
"description": "Default choice",
|
|
1534
|
+
"default": true
|
|
1535
|
+
}
|
|
1536
|
+
]
|
|
1537
|
+
},
|
|
1538
|
+
{
|
|
1539
|
+
"name": "rejuvflag",
|
|
1540
|
+
"description": "Sets whether to use the orginal prescription for mixing of main-sequence stars (based on equation 80 of <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a>) or whether to use the ratio of the pre-merger He core mass at the base of the giant branch to the merger product's He core mass at the base of the giant branch",
|
|
1541
|
+
"type": "dropdown",
|
|
1542
|
+
"options-preface": "",
|
|
1543
|
+
"options": [
|
|
1544
|
+
{
|
|
1545
|
+
"name": 0,
|
|
1546
|
+
"description": "No modifications to BSE",
|
|
1547
|
+
"default": true
|
|
1548
|
+
},
|
|
1549
|
+
{
|
|
1550
|
+
"name": 1,
|
|
1551
|
+
"description": "modified mixing times"
|
|
1552
|
+
}
|
|
1553
|
+
]
|
|
1554
|
+
},
|
|
1555
|
+
{
|
|
1556
|
+
"name": "bhms_coll_flag",
|
|
1557
|
+
"description": "If set to 1, then the star is not destroyed in a BH+star collision if \\(M_{\\rm star} > M_{\\rm BH}\\)",
|
|
1558
|
+
"type": "dropdown",
|
|
1559
|
+
"options-preface": "",
|
|
1560
|
+
"options": [
|
|
1561
|
+
{
|
|
1562
|
+
"name": 0,
|
|
1563
|
+
"description": "Star is destroyed in a BH+star collision even if \\(M_{\\rm star} > M_{\\rm BH}\\)",
|
|
1564
|
+
"default": true
|
|
1565
|
+
},
|
|
1566
|
+
{
|
|
1567
|
+
"name": 1,
|
|
1568
|
+
"description": "Star is not destroyed in a BH+star collision if \\(M_{\\rm star} > M_{\\rm BH}\\)"
|
|
1569
|
+
}
|
|
1570
|
+
]
|
|
1571
|
+
},
|
|
1572
|
+
{
|
|
1573
|
+
"settings-section": "Magnetic Braking",
|
|
1574
|
+
"name": "htpmb",
|
|
1575
|
+
"description": "Activates different models for magnetic braking",
|
|
1576
|
+
"type": "dropdown",
|
|
1577
|
+
"options-preface": "",
|
|
1578
|
+
"options": [
|
|
1579
|
+
{
|
|
1580
|
+
"name": -1,
|
|
1581
|
+
"description": "No magnetic braking"
|
|
1582
|
+
},
|
|
1583
|
+
{
|
|
1584
|
+
"name": 0,
|
|
1585
|
+
"description": "No modifications to BSE"
|
|
1586
|
+
},
|
|
1587
|
+
{
|
|
1588
|
+
"name": 1,
|
|
1589
|
+
"description": "follows <a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2003ApJ...599..516I/abstract'>Ivanona and Taam 2003</a>",
|
|
1590
|
+
"default": true
|
|
1591
|
+
}
|
|
1592
|
+
]
|
|
1593
|
+
},
|
|
1594
|
+
{
|
|
1595
|
+
"settings-section": "Miscellaneous",
|
|
1596
|
+
"name": "ST_cr",
|
|
1597
|
+
"description": "Activates different convective vs radiative boundaries",
|
|
1598
|
+
"type": "dropdown",
|
|
1599
|
+
"options-preface": "",
|
|
1600
|
+
"options": [
|
|
1601
|
+
{
|
|
1602
|
+
"name": 0,
|
|
1603
|
+
"description": "Follows BSE (<a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2002MNRAS.329..897H/abstract'>Hurley+2002</a>)"
|
|
1604
|
+
},
|
|
1605
|
+
{
|
|
1606
|
+
"name": 1,
|
|
1607
|
+
"description": "Follow StarTrack (<a class='reference external' href='https://ui.adsabs.harvard.edu/abs/2008ApJS..174..223B/abstract'>Belczynski+2008</a>)",
|
|
1608
|
+
"default": true
|
|
1609
|
+
}
|
|
1610
|
+
]
|
|
1611
|
+
},
|
|
1612
|
+
{
|
|
1613
|
+
"name": "rtmsflag",
|
|
1614
|
+
"description": "Flag for calculating the radius at the end of the main sequence.",
|
|
1615
|
+
"type": "number",
|
|
1616
|
+
"options-preface": "",
|
|
1617
|
+
"options": [
|
|
1618
|
+
{
|
|
1619
|
+
"name": 0,
|
|
1620
|
+
"description": "Use the original SSE calculation <b>except</b> if M>200Msun and Z<0.0008 (~0.04 Zsun). For the exception, an extrapolation is used. This ad hoc extrapolation works well for stars with M < 4000 Msun for Z >= 0.01 Zsun. For lower metallicities, you may still run into issues of negative radii for very massive stars - use at your own risk.",
|
|
1621
|
+
"default": true
|
|
1622
|
+
},
|
|
1623
|
+
{
|
|
1624
|
+
"name": 1,
|
|
1625
|
+
"description": "Calculate using BoOST simulation data. BoOST metallicities = [1.1e-4, 2.1e-4, 1e-3, 2e-3] corresponding to [dwarfD, IZw18, dwarfA, SMC] models (Szecsi et al. (2022))"
|
|
1626
|
+
},
|
|
1627
|
+
{
|
|
1628
|
+
"name": 2,
|
|
1629
|
+
"description": "Calculating using BPASS simulation data. BPASS metallicities = [1e-4, *2e-4*, 1e-3, 2e-3]. NOTE : For BPASS, we used a power law to fit the rtms v/s mzams values for each metallicity. We have coded the best fit power-laws for the above BPASS metallicities here. Since, Z = 2e-4 model is not available in the BPASS tracks, we assume the same rtms v/s mzams power law as Z=1e-4 for Z=2e-4."
|
|
1630
|
+
}
|
|
1631
|
+
]
|
|
1632
|
+
}
|
|
1633
|
+
]
|
|
1634
|
+
}
|
|
1635
|
+
]
|