Skip to content

Implement rapid quenching into merging satellites#425

Open
aphearin wants to merge 3 commits into
gd_bugfixfrom
gd_bugfix_rapid_q
Open

Implement rapid quenching into merging satellites#425
aphearin wants to merge 3 commits into
gd_bugfixfrom
gd_bugfix_rapid_q

Conversation

@aphearin
Copy link
Copy Markdown
Collaborator

This PR implements a new rapid quenching feature into the gd_phot_kernels.py and gd_phot_kernels_merging.py modules.

aphearin added 3 commits May 18, 2026 15:21
… gd_phot_kernels.py and gd_phot_kernels_merging.py. This changes the API of gd_phot_kernels.py because now we need to know p_merge to compute no-merge photometry
@aphearin aphearin requested a review from zaidikumail May 18, 2026 21:08
@aphearin
Copy link
Copy Markdown
Collaborator Author

The rapid quenching model brought in with this PR implements satellite quenching as a function of p_merge, with the quenching happening directly on the distribution of stellar ages, analogous to the way burstiness is implemented with diffburst in dsps. The figure shows the effect, which is implemented just by dividing the young-star weights by 10,000 and renormalizing.

rapid_quenching_demo_with_burst_0

@aphearin aphearin changed the title Gd bugfix rapid q Implement rapid quenching into merging satellites May 18, 2026
)
t_infall = t_obs - gyr_since_infall
logmp_infall = lgmu_infall + logmhost_infall
p_merge_smooth = merging_model.get_p_merge_from_merging_params(
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yep, I also faced this problem: needing p_merge before running _phot_kern_merging for the model change, in which fq changed as a function of p_merge. I went for a hacky workaround instead.

Copy link
Copy Markdown
Collaborator

@zaidikumail zaidikumail left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This looks good to me @aphearin ! Thanks for rapidly implementing this. I will now run gradient descent with this branch.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants