Procedure - IReNA-NucAstro-QuestionsTools/Supernova-Sunburn GitHub Wiki

Introduction

We took an exsisting simulation in WinNet based on r-process in CCSN winds[1]. We then modified the parameter file and trajectory files to vary the Ye to investigate how sensitive the abundance patternes and pathways with to the Ye value. Below is a step by step guide to how we made these calculations.

Move Working_CCSN_wind_bliss.par file to /WinNet/par/ and all Traj_CPR2_.txt files to /WinNet/data/Example_data/Working_CCSN_wind_bliss/ Each Traj_CPR2_.txt was altered to change the last column to the Ye desired at the start of the calculation.

How to run

  1. Edit lines 87-97 in the Working_CCSN_wind_bliss.par files to match your directory for the reaction data files. Start with Traj_CPR2_Ye_0.20.txt. The first goal is to go to the extremes in order to maximize the impact of the neutrino interactions. You want to see high A value nuclei being produced, to ensure that the neutrinos are in fact being incorporated.

  2. Change line 114 Enuebar to match the value corresponding to the Ye value in the trajectory file. The correct Enuebar value can be found by running Ecalc.py. How this value was calculated can be found in the conclusion page of this wiki.

  3. In main WinNet directory run

    python3 makerun.py -p Working_CCSN_wind_bliss.par -r 'name for output dir'
    
  4. In /WinNet/bin/movie_script run

    python3 winnet_movie.py -i ../../'name for output dir'/ --save --output='name for movie.mp4' --interval=60
    
  5. Repeat steps 1-4 for each Traj_CPR2_Ye_*.txt and Enubar

  6. The Plot_me_Ye_range.py can be used to see the output from these, there is some documentation in that code for how to use the plotting code.

Lessons Learned

The original plan was to use:

   python3 makerun.py -p Working_CCSN_explosive_burning_nup.par -r 'name for output dir'

This file took us ~45 minutes to run per setting. This file was used to see the $\bar{\nu}$ p process, given there was no neutrino interaction incoperated. One day was spent altering the parameter file and exploring the nuflag settings. This taught us that the nuflag 3 or 4 utalize similar reaction rate information. After altering the parameter file to extremes , we did not see the breakthrough to higher A values that is a major attribute of the nutrino interaction in extreme enviroments.

NuFlag 2

NuFlag 3

NuFlag 4

NuFlag 4 had the most neutrino rates utilalized but focuses on neutral charge reaction rather than charged particle reactions. We utalized NuFlag 3 for all the simulations after seeing the very small difference between NuFlag 3 and 4. More information can be found in the WinNet manual. --more time would include risidul plot --

https://github.com/user-attachments/assets/11167400-8144-40ec-b83a-dcb5e4ccbab9

Above is the CCSN_explosive file that was changed to match MRSN parameters to look for the $\bar{\nu}$ p process. Below that will be the MRSN to directly compare the $\bar{\nu}$ p process to what we saw. After some time we discovered that turning the neutrinos on properly for this file would take more time that alloted here. Also the neutrino interactions where not in the scope or design of this file. For this reason we moved to the CCSN wind file.

https://github.com/user-attachments/assets/b7ecd9b0-fc7b-4705-9e6e-bc81009f82ff