BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

R. Pakmor

R. Pakmor contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.SRastro-ph.HEastro-ph.COastro-ph.IMgr-qc

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2025arXiv

NLTE spectral modelling for a carbon-oxygen and helium white-dwarf merger as a Ca-rich transient candidate

We carry out NLTE (non local thermodynamic equilibrium) radiative transfer simulations to determine whether explosion during the merger of a carbon-oxygen (CO) white dwarf (WD) with a helium (He) WD can reproduce the characteristic Ca II/[Ca II] and He I lines observed in Ca-rich transients. Our study is based on a 1D representation of a hydrodynamic simulation of a 0.6 $M_{\odot}$ CO + 0.4 $M_{\odot}$ He WD merger. We calculate both photospheric and nebular-phase spectra including treatment for non-thermal electrons, as is required for accurate modelling of He I and [Ca II]. Consistent with Ca-rich transients, our simulation predicts a nebular spectrum dominated by emission from [Ca II] 7291, 7324 angstrom and the Ca II near-infrared (NIR) triplet. The photospheric-phase synthetic spectrum also exhibits a strong Ca II NIR triplet, prominent optical absorption due to He I 5876 angstrom and He I 10830 angstrom in the NIR, as is commonly observed for Ca-rich transients. Overall, our results therefore suggest that CO+He WD mergers are a promising channel for Ca-rich transients. However, the current simulation overpredicts some He I features, in particular both He I 6678 and 7065 angstrom and shows a significant contribution from Ti II, which results in a spectral energy distribution that is substantially redder than most Ca-rich transients at peak. Additionally the Ca II nebular emission features are too broad. Future work should investigate if these discrepancies can be resolved by considering full 3D models and exploring a range of CO+He WD binary configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bipolar planetary nebulae from common envelope evolution of binary stars

Asymmetric shapes and evidence for binary central stars suggest a common-envelope origin for many bipolar planetary nebulae. The bipolar components of the nebulae are observed to expand faster than the rest and the more slowly expanding material has been associated with the bulk of the envelope ejected during the common-envelope phase of a stellar binary system. Common-envelope evolution in general remains one of the biggest uncertainties in binary star evolution and the origin of the fast outflow has not been explained satisfactorily. We perform three-dimensional magnetohydrodynamic simulations of common-envelope interaction with the moving-mesh code AREPO. Starting from the plunge-in of the companion into the envelope of an asymptotic giant branch star and covering hundreds of orbits of the binary star system, we are able to follow the evolution to complete envelope ejection. We find that magnetic fields are strongly amplified in two consecutive episodes. First, when the companion spirals in the envelope and, second, when it forms a contact binary with the core of the former giant star. In the second episode, a magnetically-driven, high-velocity outflow of gas is launched self-consistently in our simulations. The outflow is bipolar and the gas is additionally collimated by the ejected common envelope. The resulting structure reproduces typical morphologies and velocities observed in young planetary nebulae. We propose that the magnetic driving mechanism is a universal consequence of common envelope interaction responsible for a substantial fraction of observed planetary nebulae. Such a mechanism likely also exists in the common-envelope phase of other binary stars that lead to the formation of Type Ia supernovae, X-ray binaries and gravitational-wave merger events.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Introducing the THESAN project: radiation-magnetohydrodynamic simulations of the epoch of reionization

We introduce the THESAN project, a suite of large volume (L = 95.5 cMpc) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the IGM during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter (DM) and baryonic mass resolutions of $3.1 \times 10^6$ $M_\odot$ and $5.8 \times 10^5 $ $M_\odot$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z=5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (AREPO-RT) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying DM models, and numerical convergence. The simulations produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, stellar mass function, star formation rate density, and the mass-metallicity relation, despite the galaxy formation model being mainly calibrated at z=0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21cm experiments to accurately characterise the sources that dominate the ionizing photon budget. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The redshift evolution of the binary black hole merger rate: a weighty matter

Gravitational wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate, $R_{\mathrm{BBH}}(z)$. We make predictions for $R_{\mathrm{BBH}}(z)$ as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations, COMPAS. We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterised by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about $30\rm{M}_{\odot}$ and short delay times ($t_{\rm delay} \lesssim 1$Gyr), while the stable RLOF channel primarily forms systems with BH masses above $30\rm{M}_{\odot}$ and long delay times ($t_{\rm delay} \gtrsim 1$Gyr). We provide a new fit for the metallicity specific star-formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction of $R_{\mathrm{BBH}}(z)$. This leads to a distinct redshift evolution of $R_{\mathrm{BBH}}(z)$ for high and low primary BH masses. We furthermore find that, at high redshift, $R_{\mathrm{BBH}}(z)$ is dominated by the CE channel, while at low redshift it contains a large contribution ($\sim 40\%$) from the stable RLOF channel. Our results predict that, for increasing redshifts, BBHs with component masses above $30\rm{M}_{\odot}$ will become increasingly scarce relative to less massive BBH systems. Evidence of this distinct evolution of $R_{\mathrm{BBH}}(z)$ for different BH masses can be tested with future detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The THESAN project: Lyman-alpha emission and transmission during the Epoch of Reionization

The visibility of high-redshift Lyman-alpha emitting galaxies (LAEs) provides important constraints on galaxy formation processes and the Epoch of Reionization (EoR). However, predicting realistic and representative statistics for comparison with observations represents a significant challenge in the context of large-volume cosmological simulations. The THESAN project offers a unique framework for addressing such limitations by combining state-of-the-art galaxy formation (IllustrisTNG) and dust models with the Arepo-RT radiation-magneto-hydrodynamics solver. In this initial study we present Lyman-alpha centric analysis for the flagship simulation that resolves atomic cooling haloes throughout a (95.5 cMpc)^3 region of the Universe. To avoid numerical artifacts we devise a novel method for accurate frequency-dependent line radiative transfer in the presence of continuous Hubble flow, transferable to broader astrophysical applications as well. Our scalable approach highlights the utility of LAEs and red damping-wing transmission as probes of reionization, which reveal nontrivial trends across different galaxies, sightlines, and frequency bands that can be modelled in the framework of covering fractions. In fact, after accounting for environmental factors influencing large-scale ionized bubble formation such as redshift and UV magnitude, the variation across galaxies and sightlines mainly depends on random processes including peculiar velocities and self-shielded systems that strongly impact unfortunate rays more than others. Throughout the EoR local and cosmological optical depths are often greater than or less than unity such that the exp(-tau) behavior leads to anisotropic and bimodal transmissivity. Future surveys will benefit by targeting both rare bright objects and Goldilocks zone LAEs to infer the presence of these (un)predictable (dis)advantages.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chemodynamics of barred galaxies in cosmological simulations: On the Milky Way's quiescent merger history and in-situ bulge

We explore the chemodynamical properties of a sample of barred galaxies in the Auriga magneto-hydrodynamical cosmological zoom-in simulations, which form boxy/peanut (b/p) bulges, and compare these to the Milky Way (MW). We show that the Auriga galaxies which best reproduce the chemodynamical properties of stellar populations in the MW bulge have quiescent merger histories since redshift $z\sim3.5$: their last major merger occurs at $t_{\rm lookback}>12\,\rm Gyrs$, while subsequent mergers have a stellar mass ratio of $\leq$1:20, suggesting an upper limit of a few percent for the mass ratio of the recently proposed Gaia Sausage/Enceladus merger. These Auriga MW-analogues have a negligible fraction of ex-situ stars in the b/p region ($<1\%$), with flattened, thick disc-like metal-poor stellar populations. The average fraction of ex-situ stars in the central regions of all Auriga galaxies with b/p&#39;s is 3% -- significantly lower than in those which do not host a b/p or a bar. While the central regions of these barred galaxies contain the oldest populations, they also have stars younger than 5Gyrs (>30%) and exhibit X-shaped age and abundance distributions. Examining the discs in our sample, we find that in some cases a star-forming ring forms around the bar, which alters the metallicity of the inner regions of the galaxy. Further out in the disc, bar-induced resonances lead to metal-rich ridges in the $V_ϕ-r$ plane -- the longest of which is due to the Outer Lindblad Resonance. Our results suggest the Milky Way has an uncommonly quiet merger history, which leads to an essentially in-situ bulge, and highlight the significant effects the bar can have on the surrounding disc.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Common-envelope evolution with an asymptotic giant branch star

Common-envelope phases are decisive for the evolution of many binary systems. Of particular interest are cases with asymptotic giant branch (AGB) primary stars, because they are thought to be progenitors of various astrophysical transients. In three-dimensional hydrodynamic simulations with the moving-mesh code AREPO, we study the common-envelope evolution of a $1.0\,M_{\odot}$ early-AGB star with companions of different masses. Although the stellar envelope of the AGB star is less tightly bound than that of a red giant, we find that the release of orbital energy of the core binary is insufficient to eject more than about twenty percent of the envelope mass. Ionization energy released in the expanding envelope, however, can lead to complete envelope ejection. Because recombination proceeds largely at high optical depths in our simulations, it is likely that this effect indeed plays a significant role in the considered systems. The efficiency of mass loss and the final orbital separation of the core binary system depend on the mass ratio between the companion and the primary star. Our results suggest a linear relation between the ratio of final to initial orbital separation and this parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Long-term evolution of a magnetic massive merger product

About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $τ$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

White dwarf deflagrations for Type Iax supernovae: Polarisation signatures from the explosion and companion interaction

Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

WARPFIELD Population Synthesis: The physics of (extra-)Galactic star formation and feedback driven cloud structure and emission from sub-to-kpc scales

We present a novel method to model galactic scale star formation and the resulting emission from star clusters and the multi-phase interstellar medium. We combine global parameters, such as SFR and CMF, with {\sc warpfield} which provides a description of the feedback-driven evolution of individual star-forming regions. Our approach includes stellar evolution, stellar winds, radiation pressure, supernovae, all of which couple to the dynamical evolution of the parental cloud in a highly non-linear fashion. The heating of diffuse gas and dust is calculated self-consistently with the age, mass and density dependent escape fractions of the local star-forming regions. From this we construct the interstellar radiation field at any point in the galaxy, and we employ the multi-frequency Monte Carlo radiative transfer code {\sc polaris} to produce synthetic emission maps for the one-to-one comparison with observational data. We demonstrate the capabilities of our approach by applying the method to a Milky Way like galaxy built-up in a high-resolution cosmological MHD simulation. We give three examples. First, we compute the multi-scale distribution of electron $n_{e^-}$ and $T_{e^{-}}$ and synthesize the MW all-sky H$α$ emission. We use a multipole expansion method to show that the resulting maps are consistent with observations. Second, we predict the expected \SIII 9530~Å emission. This line is a key target of several planned large survey programs. It suffers less extinction than other diagnostic lines and provides information about star formation in very dense environments that are otherwise observationally not readily accessible. Third, we explore the effects of differential extinction as seen by an extra-galactic observer, and discuss the consequences for the correct interpretation of \Ha emission as a star-formation rate tracer at different viewing angles.(abridged)

0 citations0 reviewsNo code linkedOpen access