diff --git a/config/config_mstar01_C60.yaml b/config/config_mstar01_C60.yaml new file mode 100644 index 0000000..eedc016 --- /dev/null +++ b/config/config_mstar01_C60.yaml @@ -0,0 +1,67 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# used together with jobs/mstar01_C60.yaml +# to generate part of Fig. 1 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.1 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.01 # needed for DiskLabDisk + R0: 50 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.067 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar01_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar01_C60_001.yaml b/config/config_mstar01_C60_001.yaml new file mode 100644 index 0000000..8ea6fba --- /dev/null +++ b/config/config_mstar01_C60_001.yaml @@ -0,0 +1,68 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the disc mass reduced to 1% of the stellar mass +# used together with jobs/mstar01_C60_001.yaml +# to generate part of Fig. A2 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.1 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.001 # needed for DiskLabDisk + R0: 50 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.067 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar01_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar01_C60_moreCO.yaml b/config/config_mstar01_C60_moreCO.yaml new file mode 100644 index 0000000..76a8d1c --- /dev/null +++ b/config/config_mstar01_C60_moreCO.yaml @@ -0,0 +1,69 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the abundance of CH4 reduced in the disc and the excess C/H is partitioned to CO +# naming of the file is weird, i know... +# used together with jobs/mstar01_C60_moreCO.yaml +# to generate part of Fig. A4 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.1 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.01 # needed for DiskLabDisk + R0: 50 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.067 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.29 + CH4_frac: 0.01 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar01_C60_moreCO + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar03_C60.yaml b/config/config_mstar03_C60.yaml new file mode 100644 index 0000000..b62263b --- /dev/null +++ b/config/config_mstar03_C60.yaml @@ -0,0 +1,67 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# used together with jobs/mstar03_C60.yaml +# to generate part of Fig. 1 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.3 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.03 # needed for DiskLabDisk + R0: 110 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.332 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar03_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar03_C60_001.yaml b/config/config_mstar03_C60_001.yaml new file mode 100644 index 0000000..9fb7977 --- /dev/null +++ b/config/config_mstar03_C60_001.yaml @@ -0,0 +1,68 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the disc mass reduced to 1% of the stellar mass +# used together with jobs/mstar03_C60_001.yaml +# to generate part of Fig. A2 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.3 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.003 # needed for DiskLabDisk + R0: 110 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.332 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar03_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar03_C60_moreCO.yaml b/config/config_mstar03_C60_moreCO.yaml new file mode 100644 index 0000000..6f29782 --- /dev/null +++ b/config/config_mstar03_C60_moreCO.yaml @@ -0,0 +1,69 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the abundance of CH4 reduced in the disc and the excess C/H is partitioned to CO +# naming of the file is weird, i know... +# used together with jobs/mstar03_C60_moreCO.yaml +# to generate part of Fig. A4 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.3 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.03 # needed for DiskLabDisk + R0: 110 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.332 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.29 + CH4_frac: 0.01 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar03_C60_moreCO + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar05_C60.yaml b/config/config_mstar05_C60.yaml new file mode 100644 index 0000000..01d7529 --- /dev/null +++ b/config/config_mstar05_C60.yaml @@ -0,0 +1,67 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# used together with jobs/mstar05_C60.yaml +# to generate part of Fig. 1 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.5 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.05 # needed for DiskLabDisk + R0: 155 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.663 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar05_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar05_C60_001.yaml b/config/config_mstar05_C60_001.yaml new file mode 100644 index 0000000..c049e1b --- /dev/null +++ b/config/config_mstar05_C60_001.yaml @@ -0,0 +1,68 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the disc mass reduced to 1% of the stellar mass +# used together with jobs/mstar05_C60_001.yaml +# to generate part of Fig. A2 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.5 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.005 # needed for DiskLabDisk + R0: 155 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.663 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar05_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar05_C60_moreCO.yaml b/config/config_mstar05_C60_moreCO.yaml new file mode 100644 index 0000000..59f90ff --- /dev/null +++ b/config/config_mstar05_C60_moreCO.yaml @@ -0,0 +1,69 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the abundance of CH4 reduced in the disc and the excess C/H is partitioned to CO +# naming of the file is weird, i know... +# used together with jobs/mstar05_C60_moreCO.yaml +# to generate part of Fig. A4 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.5 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.05 # needed for DiskLabDisk + R0: 155 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 0.663 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.29 + CH4_frac: 0.01 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar05_C60_moreCO + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar07_C60.yaml b/config/config_mstar07_C60.yaml new file mode 100644 index 0000000..b44d039 --- /dev/null +++ b/config/config_mstar07_C60.yaml @@ -0,0 +1,67 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# used together with jobs/mstar07_C60.yaml +# to generate part of Fig. 1 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.7 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.07 # needed for DiskLabDisk + R0: 195 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 1.093 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar07_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar07_C60_001.yaml b/config/config_mstar07_C60_001.yaml new file mode 100644 index 0000000..dd26c0d --- /dev/null +++ b/config/config_mstar07_C60_001.yaml @@ -0,0 +1,68 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the disc mass reduced to 1% of the stellar mass +# used together with jobs/mstar07_C60_001.yaml +# to generate part of Fig. A2 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.7 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.007 # needed for DiskLabDisk + R0: 195 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 1.093 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar07_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar07_C60_moreCO.yaml b/config/config_mstar07_C60_moreCO.yaml new file mode 100644 index 0000000..1eb87b2 --- /dev/null +++ b/config/config_mstar07_C60_moreCO.yaml @@ -0,0 +1,69 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the abundance of CH4 reduced in the disc and the excess C/H is partitioned to CO +# naming of the file is weird, i know... +# used together with jobs/mstar07_C60_moreCO.yaml +# to generate part of Fig. A4 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 0.7 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.07 # needed for DiskLabDisk + R0: 195 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 1.093 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.29 + CH4_frac: 0.01 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar07_C60_moreCO + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar10_C60.yaml b/config/config_mstar10_C60.yaml new file mode 100644 index 0000000..693b280 --- /dev/null +++ b/config/config_mstar10_C60.yaml @@ -0,0 +1,67 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# used together with jobs/mstar10_C60.yaml +# to generate part of Fig. 1 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 1.0 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.10 # needed for DiskLabDisk + R0: 250 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 1.933 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar10_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar10_C60_001.yaml b/config/config_mstar10_C60_001.yaml new file mode 100644 index 0000000..0fa2e3b --- /dev/null +++ b/config/config_mstar10_C60_001.yaml @@ -0,0 +1,68 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the disc mass reduced to 1% of the stellar mass +# used together with jobs/mstar10_C60_001.yaml +# to generate part of Fig. A2 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 1.0 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.010 # needed for DiskLabDisk + R0: 250 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 1.933 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.2 + CH4_frac: 0.1 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar10_C60 + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/config/config_mstar10_C60_moreCO.yaml b/config/config_mstar10_C60_moreCO.yaml new file mode 100644 index 0000000..0b31f2c --- /dev/null +++ b/config/config_mstar10_C60_moreCO.yaml @@ -0,0 +1,69 @@ +# Config file of disc C/O vs M_star (Mah et al. 2023) +# with the abundance of CH4 reduced in the disc and the excess C/H is partitioned to CO +# naming of the file is weird, i know... +# used together with jobs/mstar10_C60_moreCO.yaml +# to generate part of Fig. A4 + +# Take care of units here! Otherwise adapt in import_config.scale_units! +# Physical parameters +config_disk: + model: DiskLabDisk + M_STAR: 1.0 # mass of central star in sun massses + ALPHA: 5.0e-4 # Alpha viscosity + ALPHAHEIGHT: 1.0e-4 + M0: 0.10 # needed for DiskLabDisk + R0: 250 # needed for DiskLabDisk + DTG_total: 0.0142 # Dust to gas ratios + static: False + evaporation: True + static_stokes: False + tau_disk: 1.0e4 + begin_photoevap: 10.0 # Myr + temp_evol: False + evap_width: 1.0e-3 + lstar: 1.933 + +chemistry: + FeH: 0.0 + use_FeH: False # use FeH proxy fit according to Bitsch, Battestini 2020 + C_frac: 0.6 # Fraction of C/H in C (will be reduced from CH4) + CO_frac: 0.29 + CH4_frac: 0.01 + CO2_frac: 0.1 + +config_planet: + model: NoAccretion + matter_removal: True + use_heat_torque: True #True + use_dynamical_torque: True + migration: False + M0_fact: 1 + a_p: 5.0 + t_0: 10.5 # time before planets start growing, total time: time_disk_0+t_0 + rho_c: 5.5 # density of the core (g/cm**3) + r_in: 0.2 + keep_peb_iso: True + use_pebiso_diffusion: False + +config_pebble_accretion: + u_frag: 5.0 # m/s + epsilon_p: 0.5 + +config_gas_accretion: + kappa_env: 0.05 + f_machida: 1 + f_disk_max: 1.0 + +# modelling parameters +defaults: + DEF_R_IN: 0.1 # inner r boundary (in AU) + DEF_R_OUT: 1000 # outer r boundary (in AU) + DEF_GRIDSIZE: 500 # radial gridsize + DEF_LIN_SPACING: False # Spacing of radial grid + DEF_TIMESTEP: 10 # Timestep in years + +output: + name: solar_mstar10_C60_moreCO + save_disk: True + save_interval: 5000 # time in years + save_disk_interval: 10 diff --git a/doc/source/Quick-Start.rst b/doc/source/Quick-Start.rst index 3f20deb..c799a5a 100644 --- a/doc/source/Quick-Start.rst +++ b/doc/source/Quick-Start.rst @@ -56,4 +56,10 @@ These examples include: * ``jobs/paper_1_grid_thorngren.yaml`` Job file, used to generate the growthtracks of planets in the first chemcomp paper (Schneider & Bitsch 2021a) * ``jobs/paper_2_Jupiter.yaml`` - Job file, used to generate Jupiter like planets for the second chemcomp paper (Schneider & Bitsch 2021b) \ No newline at end of file + Job file, used to generate Jupiter like planets for the second chemcomp paper (Schneider & Bitsch 2021b) +* ``jobs/mstar01_C60.yaml`` + Job file, used to generate a disc around a 0.1 M_Sun star (Mah, Bitsch et al. 2023). Job files for other stellar masses are also included, check out the ``jobs`` folder! These job files are to be used with the corresponding config files because all the disk parameters are defined in the config files. +* ``jobs/mstar01_C60_moreCO.yaml`` + Job file, used to generate a disc around a 0.1 M_Sun star but with more CO and less CH4 in the disc (Mah, Bitsch et al. 2023). Use this with the corresponding config file. +* ``jobs/mstar01_C60_001.yaml`` + Job file, used to generate a less massive disc around a 0.1 M_Sun star (Mah, Bitsch et al. 2023). Use this with the corresponding config file. diff --git a/jobs/mstar01_C60.yaml b/jobs/mstar01_C60.yaml new file mode 100644 index 0000000..473a840 --- /dev/null +++ b/jobs/mstar01_C60.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar01_C60 +save_disk: True diff --git a/jobs/mstar01_C60_001.yaml b/jobs/mstar01_C60_001.yaml new file mode 100644 index 0000000..4500375 --- /dev/null +++ b/jobs/mstar01_C60_001.yaml @@ -0,0 +1,10 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +output_name: solar_mstar01_C60_001 +save_disk: True diff --git a/jobs/mstar01_C60_moreCO.yaml b/jobs/mstar01_C60_moreCO.yaml new file mode 100644 index 0000000..21fbef7 --- /dev/null +++ b/jobs/mstar01_C60_moreCO.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar01_C60_moreCO +save_disk: True diff --git a/jobs/mstar03_C60.yaml b/jobs/mstar03_C60.yaml new file mode 100644 index 0000000..c771b7f --- /dev/null +++ b/jobs/mstar03_C60.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar03_C60 +save_disk: True diff --git a/jobs/mstar03_C60_001.yaml b/jobs/mstar03_C60_001.yaml new file mode 100644 index 0000000..f23e3a9 --- /dev/null +++ b/jobs/mstar03_C60_001.yaml @@ -0,0 +1,10 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +output_name: solar_mstar03_C60_001 +save_disk: True diff --git a/jobs/mstar03_C60_moreCO.yaml b/jobs/mstar03_C60_moreCO.yaml new file mode 100644 index 0000000..873aa24 --- /dev/null +++ b/jobs/mstar03_C60_moreCO.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar03_C60_moreCO +save_disk: True diff --git a/jobs/mstar05_C60.yaml b/jobs/mstar05_C60.yaml new file mode 100644 index 0000000..15b6365 --- /dev/null +++ b/jobs/mstar05_C60.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar05_C60 +save_disk: True diff --git a/jobs/mstar05_C60_001.yaml b/jobs/mstar05_C60_001.yaml new file mode 100644 index 0000000..fccec43 --- /dev/null +++ b/jobs/mstar05_C60_001.yaml @@ -0,0 +1,10 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +output_name: solar_mstar05_C60_001 +save_disk: True diff --git a/jobs/mstar05_C60_moreCO.yaml b/jobs/mstar05_C60_moreCO.yaml new file mode 100644 index 0000000..600154b --- /dev/null +++ b/jobs/mstar05_C60_moreCO.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar05_C60_moreCO +save_disk: True diff --git a/jobs/mstar07_C60.yaml b/jobs/mstar07_C60.yaml new file mode 100644 index 0000000..7a2f79b --- /dev/null +++ b/jobs/mstar07_C60.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar07_C60 +save_disk: True diff --git a/jobs/mstar07_C60_001.yaml b/jobs/mstar07_C60_001.yaml new file mode 100644 index 0000000..660ed37 --- /dev/null +++ b/jobs/mstar07_C60_001.yaml @@ -0,0 +1,10 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +output_name: solar_mstar07_C60_001 +save_disk: True diff --git a/jobs/mstar07_C60_moreCO.yaml b/jobs/mstar07_C60_moreCO.yaml new file mode 100644 index 0000000..be6d020 --- /dev/null +++ b/jobs/mstar07_C60_moreCO.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar07_C60_moreCO +save_disk: True diff --git a/jobs/mstar10_C60.yaml b/jobs/mstar10_C60.yaml new file mode 100644 index 0000000..a0c8902 --- /dev/null +++ b/jobs/mstar10_C60.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar10_C60 +save_disk: True diff --git a/jobs/mstar10_C60_001.yaml b/jobs/mstar10_C60_001.yaml new file mode 100644 index 0000000..b3b3fe9 --- /dev/null +++ b/jobs/mstar10_C60_001.yaml @@ -0,0 +1,10 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +output_name: solar_mstar10_C60_001 +save_disk: True diff --git a/jobs/mstar10_C60_moreCO.yaml b/jobs/mstar10_C60_moreCO.yaml new file mode 100644 index 0000000..306fef8 --- /dev/null +++ b/jobs/mstar10_C60_moreCO.yaml @@ -0,0 +1,15 @@ +# job file for simulations done in Mah et al. (2023) +# paired with corresponding config file + +par_1: + vals: [ 1e-4,1e-3 ] + name: ALPHA + section: config_disk + +par_2: + vals: [ 1,5,10 ] + name: u_frag + section: config_pebble_accretion + +output_name: solar_mstar10_C60_moreCO +save_disk: True