I'm an astrophysicist and research software engineer working on stellar dynamics, population synthesis, fluid dynamics, and efficient simulation frameworks.
I lead the development of the GPU-based meshless computational fluid dynamics codebase at Morpheus Fluid Ltd and maintain several astrophysical research tools including binary_c-python, SSPC, and Astrotalks.
Welcome to my new Org-mode powered website!
A while back I decided I was fed up with my HUGO-academic themed website, and a simpler approach that was fully emacs lisp/org-mode driven seemed like a fun idea.
On this website I showcase my academic research and make write-ups about research, coding, and other topics that interest me.
My posts are generally categorized into the following types:
My most recent posts are listed below
About the research projects I did during my Bachelor and Master degree at the UvA, Amsterdam
About the astrophysics research code `binary_c-python` which I've written to simulate populations of stars
About the astrophysics research code `SSPC` which I've written to combine stellar population-synthesis data with star formation histories.
An overview of all my publications is available on the publications page, but I list the most recent highlighted publications and submitted papers.
Current observations of binary black-hole (BBH) merger events show support for a feature in the primary BH-mass distribution at $\sim \,35 \, solarmass$, previously interpreted as a signature of pulsational pair-instability (PPISN) supernovae. Such supernovae are expected to map a wide range of pre-supernova carbon-oxygen (CO) core masses to a narrow range of BH masses, producing a peak in the BH mass distribution. However, recent numerical simulations place the mass location of this peak above $50\,\solarmass$. Motivated by uncertainties in the progenitor's evolution and explosion mechanism, we explore how modifying the distribution of BH masses resulting from PPISN affects the populations of gravitational-wave (GW) and electromagnetic (EM) transients. To this end, we simulate populations of isolated {BBH} systems and combine them with cosmic star-formation rates. Our results are the first cosmological BBH-merger predictions made using the binary_c rapid population synthesis framework. We find that our fiducial model does not match the observed GW peak. We can only explain the $35\,\solarmass$ peak with PPISNe by shifting the expected CO core-mass range for PPISN downwards by $\sim,15\,\solarmass$. Apart from being in tension with state-of-the art stellar models, we also find that this is likely in tension with the observed rate of hydrogen-less super-luminous supernovae. Conversely, shifting the mass range upward, based on recent stellar models, leads to a predicted third peak in the BH mass function at $\sim,64\,\solarmass$. Thus we conclude that the $\sim,35\,\solarmass$ feature is unlikely to be related to PPISN.
Mass-transfer interactions in binary stars can lead to accretion disk formation, mass loss from the system and spin-up of the accretor. To determine the trajectory of the mass-transfer stream, and whether it directly impacts the accretor, or forms an accretion disk, requires numerical simulations. The mass-transfer stream is approximately ballistic, and analytic approximations based on such trajectories are used in many binary population synthesis codes as well as in detailed stellar evolution codes. We use binary population synthesis to explore the conditions under which mass transfer takes place. We then solve the reduced three-body equations to compute the trajectory of a particle in the stream for systems with varying system mass ratio, donor synchronicity and initial stream velocity. Our results show that on average both more mass and more time is spent during mass transfer from a sub-synchronous donor than from a synchronous donor. Moreover, we find that at low initial stream velocity the asynchronous rotation of the donor leads to self-accretion over a large range of mass ratios, especially for super-synchronous donors. The stream (self-)intersects in a narrow region of parameter space where it transitions between accreting onto the donor or the accretor. Increasing the initial stream velocity leads to larger areas of the parameter space where the stream accretes onto the accretor, but also more (self-)intersection. The radii of closest approach generally increase, but the range of specific angular momenta that these trajectories carry at the radius of closest approach gets broader. Our results are made publicly available.
Binary stars evolve into chemically-peculiar objects and are a major driver of the Galactic enrichment of heavy elements. During their evolution they undergo interactions, including tides, that circularize orbits and synchronize stellar spins, impacting both individual systems and stellar populations. Using Zahn's tidal theory and Mesa main-sequence model grids, we derive the governing parameters $lambda_{lm}$ and $E_2$, and implement them in the new MINT library of the stellar population code binary_c. Our MINT equilibrium tides are 2 to 5 times more efficient than the ubiquitous bse prescriptions while the radiative-tide efficiency drops sharply with increasing age. We also implement precise initial distributions based on bias-corrected observations. We assess the impact of tides and initial orbital-parameter distributions on circularization and synchronization in eight open clusters, comparing synthetic populations and observations through a bootstrapping method. We find that changing the tidal prescription yields no statistically-significant improvement as both calculations typically lie within 0.5 $sigma$. The initial distribution, especially the primordial concentration of systems at $log_{10}(P/{ m d}) approx 0.8, e approx 0.05$ dominates the statistics even when artificially increasing tidal strength. This confirms the inefficiency of tides on the main sequence and shows that constraining tidal-efficiency parameters using the $e-log_{10}(P/{ m d})$ distribution alone is difficult or impossible. Orbital synchronization carries a more striking age-dependent signature of tidal interactions. In M35 we find twice as many synchronized rotators in our MINT calculation as with se. This measure of tidal efficiency is verifiable with combined measurements of orbital parameters and stellar spins.