Researcher profile

Debora Sijacki

Debora Sijacki contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
7topics
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

9 published item(s)

preprint2022arXiv

A disturbing FABLE of mergers, feedback, turbulence, and mass biases in simulated galaxy clusters

The use of galaxy clusters as cosmological probes often relies on understanding the properties and evolution of the intracluster medium (ICM). However, the ICM is a complex plasma, regularly stirred by mergers and feedback, with non-negligible bulk and turbulent motions and a non-thermal pressure component, making it difficult to construct a coherent and comprehensive picture. To this end, we use the FABLE simulations to investigate how the hydrostatic mass bias is affected by mergers, turbulence, and feedback. Following in detail a single, massive cluster we find the bias varies significantly over cosmic time, rarely staying at the average value found at a particular epoch. Variations of the bias at a given radius are contemporaneous with periods where outflows dominate the mass flux, either due to mergers or interestingly, at high redshift, AGN feedback. The $z=0$ ensemble median mass bias in FABLE is $\sim\!13$ per cent at $R_\mathrm{500}$ and $\sim\!15$ per cent at $R_\mathrm{200}$, but with a large scatter in individual values. In halo central regions, we see an increase in temperature and a decrease in non-thermal pressure support with cosmic time as turbulence thermalises, leading to a reduction in the mass bias within $\sim\!0.2 \, R_\mathrm{200}$. When using a fitted pressure profile, instead of the simulation data, to estimate the bias, we find there can be significant differences, particularly at larger radii and higher redshift. We therefore caution over the use of such fits in future work when comparing with the next generation of X-ray and SZ observations.

preprint2022arXiv

Blandford-Znajek jets in galaxy formation simulations: exploring the diversity of outflows produced by spin-driven AGN jets in Seyfert galaxies

Recent observations of Seyfert galaxies indicate that low power, misaligned jets can undergo significant interaction with the gas in the galactic disc and may be able to drive large-scale, multiphase outflows. We apply our novel sub-grid model for Blandford-Znajek jets to simulations of the central regions of Seyferts, in which a black hole is embedded in a dense, sub-kpc circumnuclear disc (CND) and surrounded by a dilute circumgalactic medium (CGM). We find that the variability of the accretion flow is highly sensitive both to the jet power and to the CND thermodynamics and, ultimately, is determined by the complex interplay between jet-driven outflows and backflows. Even at moderate Eddington ratios, AGN jets are able to significantly alter the thermodynamics and kinematics of CNDs and entrain up to 10% of their mass in the outflow. Mass outflow rates and kinetic powers of the warm outflowing component are in agreement with recent observations for black holes with similar bolometric luminosities, with outflow velocities that are able to reach 500 km/s. Depending on their power and direction, jets are able to drive a wide variety of large-scale outflows, ranging from light, hot and collimated structures to highly mass-loaded, multiphase, bipolar winds. This diversity of jet-driven outflows highlights the importance of applying physically motivated models of AGN feedback to realistic galaxy formation contexts. Such simulations will play a crucial role in accurately interpreting the wealth of data that next generation facilities such as JWST, SKA and Athena will provide.

preprint2022arXiv

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy.

preprint2022arXiv

Two can play at that game: constraining the role of supernova and AGN feedback in dwarf galaxies with cosmological zoom-in simulations

There is growing observational evidence for dwarf galaxies hosting active galactic nuclei (AGN), including hints of AGN-driven outflows in dwarfs. However, in the common theoretical model of galaxy formation, efficient supernova (SN) feedback is the tool of choice for regulating star formation in the low-mass regime. In this paper, we present a suite of high-resolution cosmological dwarf zoom-in simulations relaxing the assumption of strong SN feedback, with the goal to determine whether more moderate SN feedback in combination with an efficient AGN could be a suitable alternative. Importantly, we find that there are sufficient amounts of gas to power brief Eddington-limited accretion episodes in dwarfs. This leads to a variety of outcomes depending on the AGN accretion model: from no additional suppression to moderate regulation of star formation to catastrophic quenching. Efficient AGN can drive powerful outflows, depleting the gas reservoir of their hosts via ejective feedback and then maintaining a quiescent state through heating the circumgalactic medium. Moderate AGN outflows can be as efficient as the strong SN feedback commonly employed, leading to star formation regulation and HI gas masses in agreement with observations of field dwarfs. All efficient AGN set-ups are associated with overmassive black holes (BHs) compared to the (heavily extrapolated) observed BH mass - stellar mass scaling relations, with future direct observational constraints in this mass regime being crucially needed. Efficient AGN activity is mostly restricted to high redshifts, with hot, accelerated outflows and high X-ray luminosities being the clearest tell-tale signs for future observational campaigns.

preprint2021arXiv

Momentum deposition of supernovae with cosmic rays

The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We explore multiple ambient densities, magnetic fields and fractions of supernova energy deposited as CRs ($χ_{\rm CR}$), accounting for cosmic ray anisotropic diffusion and streaming. All our runs have comparable evolutions until the end of the Sedov-Taylor phase. However, our CRMHD simulations experience an additional CR pressure-driven snowplough phase once the CR energy dominates inside the SNR. We present a model for the final momentum deposited by supernovae that captures this new phase: $p_{\rm SNR} = 2.87\times 10^{5} (χ_{\text{CR}} + 1)^{4.82}\left(\frac{n}{\text{cm}^{-3}}\right)^{-0.196} M_{\odot}$ km s$^{-1}$. Assuming a 10% fraction of SN energy in CRs leads to a 50% boost of the final momentum, with our model predicting even higher impacts at lower ambient densities. The anisotropic diffusion of CRs assuming an initially uniform magnetic field leads to extended gas and cosmic ray outflows escaping from the supernova poles. We also study a tangled initial configuration of the magnetic field, resulting instead in a quasi-isotropic diffusion of CRs and earlier momentum deposition. Finally, synthetic synchrotron observations of our simulations using the POLARIS code show that the local magnetic field configuration in the interstellar medium modifies the overall radio emission morphology and polarisation.

preprint2020arXiv

Metal enrichment in the circumgalactic medium and Lyα haloes around quasars at z $\sim$3

Deep observations have detected extended Ly$α$ emission nebulae surrounding tens of quasars at redshift 2 to 6. However, the metallicity of such extended haloes is still poorly understood. We perform a detailed analysis on a large sample of 80 quasars at $z\sim3$ based on MUSE-VLT data. We find clear evidence of extended emission of the UV nebular lines such as CIV or HeII for about 20$\%$ of the sample, while CIII] is only marginally detected in a few objects. By stacking the cubes we detect emission of CIV, HeII and CIII] out to a radius of about 45 kpc. CIV and HeII show a radial decline much steeper than Ly$α$, while CIII] shows a shallower profile similar to Ly$α$ in the inner 45 kpc. We infer that the average metallicity of the circumgalactic gas within the central 30-50~kpc is $\sim$0.5 solar, or even higher. However, we also find evidence of a component of the Ly$α$ haloes, which has much weaker metal emission lines relative to Ly$α$. We suggest that the high metallicity of the circumgalactic medium within the central 30-50 kpc is associated with chemical pre-enrichment by past quasar-driven outflows and that there is a more extended component of the CGM that has much lower metallicity and likely associated with near-pristine gas accreted from the intergalactic medium. We show that our observational results are in good agreement with the expectations of the FABLE zoom-in cosmological simulations.

preprint2013arXiv

Following the flow: tracer particles in astrophysical fluid simulations

We present two numerical schemes for passive tracer particles in the hydrodynamical moving-mesh code AREPO, and compare their performance for various problems, from simple setups to cosmological simulations. The purpose of tracer particles is to allow the flow to be followed in a Lagrangian way, tracing the evolution of the fluid with time, and allowing the thermodynamical history of individual fluid parcels to be recorded. We find that the commonly-used `velocity field tracers', which are advected using the fluid velocity field, do not in general follow the mass flow correctly, and explain why this is the case. This method can result in orders-of-magnitude biases in simulations of driven turbulence and in cosmological simulations, rendering the velocity field tracers inappropriate for following these flows. We then discuss a novel implementation of `Monte Carlo tracers', which are moved along with fluid cells, and are exchanged probabilistically between them following the mass flux. This method reproduces the mass distribution of the fluid correctly. The main limitation of this approach is that it is more diffusive than the fluid itself. Nonetheless, we show that this novel approach is more reliable than what has been employed previously and demonstrate that it is appropriate for following hydrodynamical flows in mesh-based codes. The Monte Carlo tracers can also naturally be transferred between fluid cells and other types of particles, such as stellar particles, so that the mass flow in cosmological simulations can be followed in its entirety.

preprint2013arXiv

Moving mesh cosmology: tracing cosmological gas accretion

We investigate the nature of gas accretion onto haloes and galaxies at z=2 using cosmological hydrodynamic simulations run with the moving mesh code AREPO. Implementing a Monte Carlo tracer particle scheme to determine the origin and thermodynamic history of accreting gas, we make quantitative comparisons to an otherwise identical simulation run with the smoothed particle hydrodynamics (SPH) code GADGET-3. Contrasting these two numerical approaches, we find significant physical differences in the thermodynamic history of accreted gas in haloes above 10^10.5 solar masses. In agreement with previous work, GADGET simulations show a cold fraction near unity for galaxies forming in massive haloes, implying that only a small percentage of accreted gas heats to an appreciable fraction of the virial temperature during accretion. The same galaxies in AREPO show a much lower cold fraction, <20% in haloes above 10^11 solar masses. This results from a hot gas accretion rate which, at this same halo mass, is an order of magnitude larger than with GADGET, while the cold accretion rate is also lower. These discrepancies increase for more massive systems, and we explain both as due to numerical inaccuracies in the standard formulation of SPH. We also observe that the relatively sharp transition from cold to hot mode dominated accretion, at a halo mass of ~10^11, is a consequence of comparing past gas temperatures to a constant threshold value independent of virial temperature. Examining the spatial distribution of accreting gas, we find that gas filaments in GADGET tend to remain collimated and flow coherently to small radii, or artificially fragment and form a large number of purely numerical &#34;blobs&#34;. Similar gas streams in AREPO show increased heating and disruption at 0.25-0.5 virial radii and contribute to the hot gas accretion rate in a manner distinct from classical cooling flows.

preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.