Add Yang+2024 parametric SIDM gravothermal-sequence profile

autogalaxy/config/priors/mass/dark/yang24.yaml

@@ -0,0 +1,51 @@
+YangSIDMSph:
+  centre_0:
+    type: Gaussian
+    mean: 0.0
+    sigma: 0.1
+    width_modifier:
+      type: Absolute
+      value: 0.05
+    limits:
+      lower: -inf
+      upper: inf
+  centre_1:
+    type: Gaussian
+    mean: 0.0
+    sigma: 0.1
+    width_modifier:
+      type: Absolute
+      value: 0.05
+    limits:
+      lower: -inf
+      upper: inf
+  kappa_s:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 1.0
+    width_modifier:
+      type: Relative
+      value: 0.2
+    limits:
+      lower: 0.0
+      upper: inf
+  scale_radius:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 30.0
+    width_modifier:
+      type: Relative
+      value: 0.2
+    limits:
+      lower: 0.0
+      upper: inf
+  tau:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 1.0
+    width_modifier:
+      type: Relative
+      value: 0.2
+    limits:
+      lower: 0.0
+      upper: 1.0

autogalaxy/config/priors/mass/dark/yang24_mcr.yaml

@@ -0,0 +1,91 @@
+YangSIDMMCRLudlowSph:
+  centre_0:
+    type: Gaussian
+    mean: 0.0
+    sigma: 0.1
+    width_modifier:
+      type: Absolute
+      value: 0.05
+    limits:
+      lower: -inf
+      upper: inf
+  centre_1:
+    type: Gaussian
+    mean: 0.0
+    sigma: 0.1
+    width_modifier:
+      type: Absolute
+      value: 0.05
+    limits:
+      lower: -inf
+      upper: inf
+  mass_at_200:
+    type: LogUniform
+    lower_limit: 100000000.0
+    upper_limit: 1000000000000000.0
+    width_modifier:
+      type: Relative
+      value: 0.5
+    limits:
+      lower: 0.0
+      upper: inf
+  sigma_over_m:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 10.0
+    width_modifier:
+      type: Relative
+      value: 0.5
+    limits:
+      lower: 0.0
+      upper: inf
+  velocity_exponent:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 4.0
+    width_modifier:
+      type: Relative
+      value: 0.5
+    limits:
+      lower: 0.0
+      upper: 8.0
+  velocity_ref:
+    type: Uniform
+    lower_limit: 1.0
+    upper_limit: 100.0
+    width_modifier:
+      type: Relative
+      value: 0.5
+    limits:
+      lower: 0.0
+      upper: inf
+  t_age:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 13.8
+    width_modifier:
+      type: Relative
+      value: 0.2
+    limits:
+      lower: 0.0
+      upper: inf
+  redshift_object:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 1.0
+    width_modifier:
+      type: Relative
+      value: 0.5
+    limits:
+      lower: 0.0
+      upper: inf
+  redshift_source:
+    type: Uniform
+    lower_limit: 0.0
+    upper_limit: 1.0
+    width_modifier:
+      type: Relative
+      value: 0.5
+    limits:
+      lower: 0.0
+      upper: inf

autogalaxy/profiles/mass/__init__.py

@@ -48,6 +48,8 @@
     cNFWMCRScatterLudlowSph,
     KaplinghatCoredNFWSph,
     KaplinghatCoredNFWMCRLudlowSph,
+    YangSIDMSph,
+    YangSIDMMCRLudlowSph,
 )
 from .stellar import (
     Gaussian,

autogalaxy/profiles/mass/dark/__init__.py

@@ -14,3 +14,5 @@
 from .cnfw_mcr_scatter import cNFWMCRScatterLudlow, cNFWMCRScatterLudlowSph
 from .kaplinghat import KaplinghatCoredNFWSph
 from .kaplinghat_mcr import KaplinghatCoredNFWMCRLudlowSph
+from .yang24 import YangSIDMSph
+from .yang24_mcr import YangSIDMMCRLudlowSph