Calculation Logs

= Calculation Log =

NOTE: Figures and data are basically stored in Dropbox.

PPN
- 10/9-10/2019: Calculated overproduction factor (abundance ratio to the case with the Kappeler 1994 rate) as per request of Greg and Grisha. 5 solar mass with solar system metallicity (asplund). Trajectory.inpu was obtained from star.log. See 191010. The results do not agree much with Heshani's NEWTON calculations. But I checked trajectories, etc, carefully, so my calculations should be ok.

- 10/12: Since I found the calculations above are not right (used Z=0.01 Asplund 2005 abundance), I recalculated with Z=0.02 abundance files (GN93). The results still don't match with Heshani's calculations, but they changed quite a bit from Z=0.01 ones.

- 10/13: Run XRB examples.Learned many heavy p-rich nuclei are actually produced by successive (p,g) on stable heavy nuclei (not simply started from light stable nuclei).

- 10/27: Run XRB examples. Tested changing a reaction (59Cu(p,g)) by factor of 100 (in ppn_physics.input) and see if the abundances are changed. I confirmed it does. However, only neighborning nuclei are changed and no effect on lighter or much heavier nuclei. This is because I'm running in single zone and energy generation is not affected by this reaction change (which should be a lot; see Cybert APJ 2016). Also I wrote some scripts in analysis directory.

MPPNP
- 10/10-14/2019 (GCA): Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are (probably) set to Kappeler1994??. Because PPN results didn't agree with Talwar or Heshani either, so I just hoped MPPNP might improve the disagreement.

- 10/17-22/2019 - (WENDI) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Longland2012.

- 10/17-21/2019 - (GCA) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Talwar2016.

- 10/18-19/2019 - (GC-Master) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Massimi2015.

- 10/20/2019 - (GC-Master) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Kappeler1994 (to check 10/10-14 run is truly Kappeler1994).

- 10/24-25/2019 (GCA): Recalculated 5 M Z=0.02 (initially 10/17-21) because some outputs (hdf5) have bad data (missing some cycles; 385001, 435001, ect). So tried by just running these cycles.

- 10/30-31/2019 (GC-MASTER): 25M Z=0.02 from the Canfar hdf5 data. To see s-process impact on massive stars. 22Ne(a,n/g) rates set to Talwar2016.

- 11/1/2019 (GC-MASTER): 25M Z=0.02 from the Canfar hdf5 data. To see s-process impact on massive stars. 22Ne(a,n/g) rates set to Longland2012.

- 11/1/2019 (GC-MASTER): 25M Z=0.02 from the Canfar hdf5 data. To see s-process impact on massive stars. 22Ne(a,n/g) rates set to Kappeler1994.

- 11/2/2019 (GC-MASTER): 25M Z=0.02 from the Canfar hdf5 data. To see s-process impact on massive stars. 22Ne(a,n/g) rates set to Massimi2015.

- 11/3/2019 (GC-MASTER): 25M Z=0.02 from the Canfar hdf5 data. To see s-process impact on massive stars. 22Ne(a,n/g) rates set to TAMU2019 (paper ver. 2019/11/3).

- 11/3/2019 (GC-MASTER): 25M Z=0.02 from the Canfar hdf5 data. To see s-process impact on massive stars. 22Ne(a,n/g) rates set to TAMU2019 (a,n) Upperlimit, (a,g) lowerlimit (paper ver. 2019/11/3).

>>>>>>>>>>>>>>>>>>> 11/3/2019, 25M Z=0.02 completed!

- 11/4/2019 - (GC-Master) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019 (paper ver. 2019/11/3).

- 11/5/2019 - (GC-Master) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019 UL ((a,n) Upperlimit, (a,g) lowerlimit) (paper ver. 2019/11/3).

>>>>>>>>>>>>>>>>>>> 11/5/2019, 5M Z=0.02 completed!

- 11/6-7/2019 - (GC-Master) : Calculated 3M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Talwar2016. '''NOTE Canfar data do not match the .idx file. For example, 70001.hd5 starts with the cycle 70000. So this causes errors (I needed to change the file name to 70000.hd5 and change in .idx file too. Also NOTE this canfar has a mistake in .idx file. It doesn't include until last cycle, so you need to add the whole cycles in .idx file otherwise calculations are interrupted.'''

- 11/6-9/2019 - (GCA) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Talwar2016. '''NOTE Canfar data has a problem in cycle 60000. Both 590001.h5 an 60001.h5 files have cycle 60000, but 59001.h5's one is broken. So I deleted it and saved in hdf5view and also changed the name of 60001.h5 to 60000.h5 and as well as in .idx file. Important things to know are Canfar data are not simply read-only (because cycle 60000 is added if the calculations failed) and name in .idx is important. I.e., if I keep 60001.h5 name then, the error occurs because 60000 is empty and the file says the first cycle is 60001 even 60000 exists in the file. So be careful there's no conflict in file names (file name is read by the mppnp to know about the cycle). '''

- 11/6-8/2019 - (WENDI) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Longland2012.

- 11/8-9/2019 - (GC-Master) : Calculated 3M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Longland2012.

- 11/8-10/2019 - (WENDI) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Massimi2015.'''NOTE: Had the same problem above as GCA. Fixed cycle 60000 and run well. (don't know why longland2012 run with Wendi was ok above...) IMPORTANT!! run stopped before cycle 90000. Almost same reason above (h5 file problem). I tried to fix it, but wget or google drive download didn't work (internet connection??), so I couldn't. I'll fix it later.'''

- 11/9-10/2019 - (GC-Master) : Calculated 3M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to NACRE1999.

- 11/9-12/2019 - (GCA) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to NACRE1999.

- 11/10-12/2019 - (WENDI) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Kp1994.

>>>>>>>>>>>>>>>>>>> '11/11/2019, WENDI data all deleted by Fork mistakenly.... Need to rerun for all WENDI runs up to this day '

- 11/10-11/2019 - (GC-Master) : Calculated 3M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Massimi2015.

- 11/11-12/2019 - (GC-Master) : Calculated 3M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Kp1994.

- 11/12-13/2019 - (GC-Master) : Calculated 3M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019.''This might be wrong (kp1994) because seemingly I didn't change vital.F90... !!Check by analyzing data 112119 and seems ok!''

- 11/12-15/2019 - (GCA) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019.

>>>>>>>>>>>>>>>>>>> '11/13/2019, restart using WENDI after chat with Fork. Now data are saved in /home/user/user/scratch14_wendi3/sota/; do the same calculations for lost data first '

- 11/13-15?/2019 - (WENDI) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Kp1994. ''Cancelled! Too slow (maybe someone's using it).''

- 11/14-15/2019 - (GC-Master) : Calculated 3M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019 UL ((a,n) Upperlimit, (a,g) lowerlimit) (paper ver. 2019/11/3) ''This might be wrong (kp1994) because seemingly I didn't change vital.F90... !!Check by analyzing data 112119 and seems ok!''

>>>>>>>>>>>>>>>>>>> '11/14/2019, start using WENDI Backup server (astrohub.phys.uvic.ca) after chat with Falk (I call this Wendi-b). Now data are saved in /home/user/user/scratch14_wendi3/sota/; do the same calculations for lost data first '

- 11/14-16?/2019 - (WENDI-b) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Massimi2016. ''Cancelled! Too slow (maybe someone's using it).''

>>>>>>>>>>>>>>>>>>> 11/15/2019, 3M Z=0.02 completed!

- 11/15-16/2019 - (GC-Master) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to NACRE1999. This might be wrong (kp1994) because seemingly I didn't change vital.F90... Run again 11/16-17 (surf 000501.surf.h5 has problem. Replaced with Talwar2016 one for now)

- 11/15-18/2019 - (GCA) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019 UL ((a,n) Upperlimit, (a,g) lowerlimit) (paper ver. 2019/11/3).

- 11/16/2019 - (GC-Master) : Calculated 25M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to NACRE1999.

- 11/18/2019 - (GC-Master): Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Longland2012.

- 11/19/2019 - (GC-Master): Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Kp1994.

- 11/19-21/2019 - (GCA) : Calculated 25M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Massimi2016.

- 11/20-21/2019 - (GC-Master): Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Kp1994.

- 11/20-23/2019 - (cycdhcp39+SonyLinux) : Calculated 5M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Massimi2016.

- 11/20-??/2019 - (Wendi-b) : Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to ???.

- 11/21-22/2019 - (GC-Master): Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Massimi2016.

- 11/21-23?/2019 - (Wendi): Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to NACRE1999.

- 11/22-24/2019 - (GCA) : Calculated 25M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to NACRE1999.

- 11/22-23/2019 - (GC-Master) : Calculated 5M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Longland2012.

- 11/23-26/2019 - (cycdhcp39) : Calculated 25M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Kp1994. '''SonyLinux didn't work due to network issue. Stand-alone now. '''

- 11/23-24/2019 - (GC-Master) : Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Longland2012.

- 11/24-26/2019 - (GCA) : Calculated 25M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Talwar2016.

- 11/24-25/2019 - (GC-Master) : Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Talwar2016.

- 11/25-26/2019 - (GC-Master) : Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to NACRE1999.

- 11/26-28/2019 - (GCA) : Calculated 25M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to Longland2012.

- 11/26-29/2019 - (cycdhcp39) : Calculated 25M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019. '''SonyLinux didn't work due to network issue. Stand-alone now. '''

- 11/27-28/2019 - (GC-Master) : Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019.

- 11/28-29/2019 - (GC-Master) : Calculated 3M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019UL.

>>>>>>>>>>>>>>>>>>> 11/29/2019, 3M Z=0.01 completed!

- 11/28-30/2019 - (GCA) : Calculated 25M Z=0.01 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019UL.

>>>>>>>>>>>>>>>>>>> 11/30/2019, 25M Z=0.01 completed!

- 11/29/2019 - (GC-Master) : Calculated 1.65M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019.

- 11/29-12/1/2019 - (cycdhcp39) : Calculated 12M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019.

- 11/30/2019 - (GC-Master) : Calculated 7M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019.

- 11/30-12/1/2019 - (GC-Master) : Calculated 8M Z=0.02 from the Canfar hdf5 data. 22Ne(a,n/g) rates are set to TAMU2019.

MESA
- 6/23/2019: Calculated neutron abundance in 25 solar mass stars. See 190623

- 10/9/: Calculated 5 solar mass AGB (based on 7M_prems_to_AGB). Checked HR diagram and T-rho curves, etc. Confirmed the T-rho trajectory agrees very well with Canfar data's trajectory.

- 10/11: Calculated nova and 25M_pre_ms_to_core_collapse from test_suite on WENDI (just for curiosity). Nova trajectory seems very different from XRB trajectory on PPN examples. Also, 25M trajectory is not similar to 25M_ccsn on PPN examples (but this is probably PPN shows C or He burning's trajectory).

- 10/12: Calculated C13_pocket and He_core_flash from test_suite on WENDI (for curiosity).

- 11/30: Calculation 8M Z=0.02 for creating .se.h5 input files for MPPNP.