From 125e0cb5b2147cd675eac0efb0d8c9a68791cb1f Mon Sep 17 00:00:00 2001 From: Khanak Bhargava Date: Fri, 1 Dec 2023 15:52:08 -0500 Subject: [PATCH 1/6] Adding a readme for wdmergers --- Exec/science/wdmerger/README | 17 +++++++++++++++++ 1 file changed, 17 insertions(+) create mode 100644 Exec/science/wdmerger/README diff --git a/Exec/science/wdmerger/README b/Exec/science/wdmerger/README new file mode 100644 index 0000000000..d73d40cc0b --- /dev/null +++ b/Exec/science/wdmerger/README @@ -0,0 +1,17 @@ +# This contains information about the input parameters that can be used for running WD mergers + +# Easiest way to start a merger +# Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E +problem.roche_radius_factor = 1.0e0 + +# To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun,use +problem.hybrid_wd_c_frac = 0.00e0 +problem.hybrid_wd_o_frac = 0.00e0 +problem.hybrid_wd_he_shell_mass = 0.00e0 + +# To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use +problem.co_wd_c_frac = 0.00e0 +problem.co_wd_o_frac = 0.00e0 +problem.co_wd_he_shell_mass = 0.00e0 + + From e74f58cec3e3e5e3c78309b09a282eec4f290694 Mon Sep 17 00:00:00 2001 From: Khanak Bhargava Date: Fri, 1 Dec 2023 17:32:57 -0500 Subject: [PATCH 2/6] Changed to markdown --- Exec/science/wdmerger/README.md | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) create mode 100644 Exec/science/wdmerger/README.md diff --git a/Exec/science/wdmerger/README.md b/Exec/science/wdmerger/README.md new file mode 100644 index 0000000000..b8e15ff0c1 --- /dev/null +++ b/Exec/science/wdmerger/README.md @@ -0,0 +1,20 @@ +# Input parameters to simulate WD mergers + +### Easiest way to start a merger +### Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E +``` +problem.roche_radius_factor = 1.0e0 +``` +### To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun,use +``` +problem.hybrid_wd_c_frac = 0.00e0 +problem.hybrid_wd_o_frac = 0.00e0 +problem.hybrid_wd_he_shell_mass = 0.00e0 +``` +### To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use +``` +problem.co_wd_c_frac = 0.00e0 +problem.co_wd_o_frac = 0.00e0 +problem.co_wd_he_shell_mass = 0.00e0 +``` + From 43bfaa89a5666ff45141d2c27b1276767a4b4f3c Mon Sep 17 00:00:00 2001 From: Khanak Bhargava Date: Fri, 1 Dec 2023 19:27:51 -0500 Subject: [PATCH 3/6] Updated WD Readme --- Exec/science/wdmerger/README.md | 12 ++++++++---- 1 file changed, 8 insertions(+), 4 deletions(-) diff --git a/Exec/science/wdmerger/README.md b/Exec/science/wdmerger/README.md index b8e15ff0c1..2ad7c332f5 100644 --- a/Exec/science/wdmerger/README.md +++ b/Exec/science/wdmerger/README.md @@ -1,17 +1,21 @@ # Input parameters to simulate WD mergers -### Easiest way to start a merger -### Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E +### Easiest way to start a merger: + ``` problem.roche_radius_factor = 1.0e0 ``` -### To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun,use +### This is used to calculate the orbital distance and the rotational period. (Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E) + +### Castro either uses the inital WD masses or the initial central densities to generate desired stars. Users can specify one of these and leave the other parameter negative. If both are provided, it will use the central densities to generate the model and assume masses as an estimate for the purpose of determining the envelope mass boundary. This is established [here](https://github.com/AMReX-Astro/Castro/blob/10534968cf09ba73f1044186389740559adc4188/Util/model_parser/model_parser.H#L208). + +### To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun, use: ``` problem.hybrid_wd_c_frac = 0.00e0 problem.hybrid_wd_o_frac = 0.00e0 problem.hybrid_wd_he_shell_mass = 0.00e0 ``` -### To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use +### To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use: ``` problem.co_wd_c_frac = 0.00e0 problem.co_wd_o_frac = 0.00e0 From 9f13cb0a2d9af0cdb9ae5897252249212ba25188 Mon Sep 17 00:00:00 2001 From: Khanak Bhargava Date: Sat, 2 Dec 2023 09:41:27 -0500 Subject: [PATCH 4/6] Fixed format --- Exec/science/wdmerger/README | 17 ----------------- Exec/science/wdmerger/README.md | 10 +++++----- 2 files changed, 5 insertions(+), 22 deletions(-) delete mode 100644 Exec/science/wdmerger/README diff --git a/Exec/science/wdmerger/README b/Exec/science/wdmerger/README deleted file mode 100644 index d73d40cc0b..0000000000 --- a/Exec/science/wdmerger/README +++ /dev/null @@ -1,17 +0,0 @@ -# This contains information about the input parameters that can be used for running WD mergers - -# Easiest way to start a merger -# Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E -problem.roche_radius_factor = 1.0e0 - -# To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun,use -problem.hybrid_wd_c_frac = 0.00e0 -problem.hybrid_wd_o_frac = 0.00e0 -problem.hybrid_wd_he_shell_mass = 0.00e0 - -# To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use -problem.co_wd_c_frac = 0.00e0 -problem.co_wd_o_frac = 0.00e0 -problem.co_wd_he_shell_mass = 0.00e0 - - diff --git a/Exec/science/wdmerger/README.md b/Exec/science/wdmerger/README.md index 2ad7c332f5..fa00c49f5c 100644 --- a/Exec/science/wdmerger/README.md +++ b/Exec/science/wdmerger/README.md @@ -1,21 +1,21 @@ # Input parameters to simulate WD mergers -### Easiest way to start a merger: +Easiest way to start a merger: ``` problem.roche_radius_factor = 1.0e0 ``` -### This is used to calculate the orbital distance and the rotational period. (Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E) +This is used to calculate the orbital distance and the rotational period. (Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E) -### Castro either uses the inital WD masses or the initial central densities to generate desired stars. Users can specify one of these and leave the other parameter negative. If both are provided, it will use the central densities to generate the model and assume masses as an estimate for the purpose of determining the envelope mass boundary. This is established [here](https://github.com/AMReX-Astro/Castro/blob/10534968cf09ba73f1044186389740559adc4188/Util/model_parser/model_parser.H#L208). +Castro either uses the inital WD masses or the initial central densities to generate desired stars. Users can specify one of these and leave the other parameter negative. If both are provided, it will use the central densities to generate the model and assume masses as an estimate for the purpose of determining the envelope mass boundary. This is established [here](https://github.com/AMReX-Astro/Castro/blob/10534968cf09ba73f1044186389740559adc4188/Util/model_parser/model_parser.H#L208). -### To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun, use: +To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun, use: ``` problem.hybrid_wd_c_frac = 0.00e0 problem.hybrid_wd_o_frac = 0.00e0 problem.hybrid_wd_he_shell_mass = 0.00e0 ``` -### To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use: +To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use: ``` problem.co_wd_c_frac = 0.00e0 problem.co_wd_o_frac = 0.00e0 From 61305ad1131a8fad3c852bac8ca1b29e162ff78c Mon Sep 17 00:00:00 2001 From: Khanak Bhargava Date: Wed, 14 Feb 2024 14:01:01 -0500 Subject: [PATCH 5/6] Added info to readme --- Exec/science/wdmerger/README.md | 58 ++++++++++++++++++++++++++++----- 1 file changed, 49 insertions(+), 9 deletions(-) diff --git a/Exec/science/wdmerger/README.md b/Exec/science/wdmerger/README.md index fa00c49f5c..082cf6ee26 100644 --- a/Exec/science/wdmerger/README.md +++ b/Exec/science/wdmerger/README.md @@ -1,24 +1,64 @@ # Input parameters to simulate WD mergers -Easiest way to start a merger: +### Problem set up +``` +problem.problem = 1 +``` +### Easiest way to make the two stars merge ``` problem.roche_radius_factor = 1.0e0 ``` This is used to calculate the orbital distance and the rotational period. (Ref: https://articles.adsabs.harvard.edu/pdf/1983ApJ...268..368E) +### Desired masses or densities + Castro either uses the inital WD masses or the initial central densities to generate desired stars. Users can specify one of these and leave the other parameter negative. If both are provided, it will use the central densities to generate the model and assume masses as an estimate for the purpose of determining the envelope mass boundary. This is established [here](https://github.com/AMReX-Astro/Castro/blob/10534968cf09ba73f1044186389740559adc4188/Util/model_parser/model_parser.H#L208). -To set the CO mass fractions and He shell mass for WDs with masses between 0.45 M_sun and 0.60 M_sun, use: +**Primary** +``` +problem.mass_P = 1.0e0 +problem.central_density_P = -1.0e0 +``` + +**Secondary** +``` +problem.mass_S = 1.0e0 +problem.central_density_S = -1.0e0 +``` + +### Desired compositions + +Inspired by Dan et al. 2012, we classify the WDs using their masses and the maximum values that they can hold, and set their compositions accordingly. + +Stars with masses less than 0.45 $M_{\odot}$ are given a pure He composition. This maximum limit can be changed using: ``` -problem.hybrid_wd_c_frac = 0.00e0 -problem.hybrid_wd_o_frac = 0.00e0 -problem.hybrid_wd_he_shell_mass = 0.00e0 +problem.max_he_wd_mass = 0.45e0 ``` -To set the CO mass fractions and He shell mass for WDs with masses between 0.6 M_sun and 1.05 M_sun, use: + +To simulate stars with pure CO cores and a thin He shell between 0.45 $M_{\odot}$ and 0.6 $M_{\odot}$ we use: + ``` -problem.co_wd_c_frac = 0.00e0 -problem.co_wd_o_frac = 0.00e0 -problem.co_wd_he_shell_mass = 0.00e0 +problem.max_he_wd_mass = 0.45e0 +problem.max_hybrid_wd_mass = 0.6e0 + +problem.hybrid_wd_c_frac = 0.50e0 +problem.hybrid_wd_o_frac = 0.50e0 +problem.hybrid_wd_he_shell_mass = 0.10e0 ``` +On the other hand, to simlute stars with pure CO cores and a He shell between 0.6 $M_{\odot}$ and 1.05 $M_{\odot}$, and with different CO mass fractions, we use: +``` +problem.max_hybrid_wd_mass = 0.6e0 +problem.max_co_wd_mass = 1.05e0 + +problem.co_wd_c_frac = 0.40e0 +problem.co_wd_o_frac = 0.60e0 +problem.co_wd_he_shell_mass = 0.0e0 +``` +**NOTE:** To create stars with masses higher than the default maximum values, we need too increase the limits otherwise the conditions will not be satisfied and hence, we won't get the desired compositions. +### Stellar temperature +We can set the initial stellar temperature to the desired value using +``` +problem.stellar_temp = 1.0e7 +``` \ No newline at end of file From 63281fa6a328ff39c09f1537e1a84aa00dd15157 Mon Sep 17 00:00:00 2001 From: Khanak Bhargava Date: Wed, 14 Feb 2024 16:18:52 -0500 Subject: [PATCH 6/6] Codespell fix --- Exec/science/wdmerger/README.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Exec/science/wdmerger/README.md b/Exec/science/wdmerger/README.md index 082cf6ee26..5c7dce023f 100644 --- a/Exec/science/wdmerger/README.md +++ b/Exec/science/wdmerger/README.md @@ -13,7 +13,7 @@ This is used to calculate the orbital distance and the rotational period. (Ref: ### Desired masses or densities -Castro either uses the inital WD masses or the initial central densities to generate desired stars. Users can specify one of these and leave the other parameter negative. If both are provided, it will use the central densities to generate the model and assume masses as an estimate for the purpose of determining the envelope mass boundary. This is established [here](https://github.com/AMReX-Astro/Castro/blob/10534968cf09ba73f1044186389740559adc4188/Util/model_parser/model_parser.H#L208). +Castro either uses the initial WD masses or the initial central densities to generate desired stars. Users can specify one of these and leave the other parameter negative. If both are provided, it will use the central densities to generate the model and assume masses as an estimate for the purpose of determining the envelope mass boundary. This is established [here](https://github.com/AMReX-Astro/Castro/blob/10534968cf09ba73f1044186389740559adc4188/Util/model_parser/model_parser.H#L208). **Primary** ```