Researcher profile

Thomas Berlok

Thomas Berlok contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
4topics
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

4 published item(s)

preprint2026arXiv

Characterizing turbulence in galaxy clusters: defining turbulent energies and assessing multi-scale versus fixed-scale filters

Disentangling turbulence and bulk motions in the intracluster medium (ICM) of galaxy clusters is inherently ambiguous, as the plasma is continuously stirred by different processes on disparate scales. This poses a serious problem in the interpretation of both observations and numerical simulations. In this paper, we use filtering operators in real space to separate bulk motion from turbulence at different scales. We show how filters can be used to define consistent kinetic and magnetic energies for the bulk and turbulent component. We apply our GPU-accelerated filtering pipeline to a simulation of a major galaxy cluster merger, which is part of the PICO-Clusters suite of zoom-in cosmological simulations of massive clusters using the moving mesh code Arepo and the IllustrisTNG galaxy formation model. We find that during the merger the turbulent pressure fraction on physical scales $\lesssim$160 kpc reaches a maximum of 5%, before decreasing to 2% after $\sim$1.3 Gyr from the core passage. These low values are consistent with recent observations of clusters with XRISM, and suggest that unless a cluster was recently perturbed by a major merger, turbulence levels are low. We then re-examine the popular multiscale iterative filter method. In our tests, we find that its use can introduce artifacts, and that it does not reliably disentangle fluctuations living on widely separated length scales. Rather, we believe it is more fruitful to use fixed-scale filters and turbulent energies to compare between simulations and observations. This work significantly improves our understanding of turbulence generation by major mergers in galaxy clusters, which can be probed by XRISM and next-generation X-ray telescopes, allowing us to connect high-resolution cosmological simulations to observations.

preprint2020arXiv

Local Simulations of Instabilities Driven by Composition Gradients in the ICM

The distribution of Helium in the intracluster medium (ICM) permeating galaxy clusters is not well constrained due to the very high plasma temperature. Therefore, the plasma is often assumed to be homogeneous. A non-uniform Helium distribution can however lead to biases when measuring key cluster parameters. This has motivated one-dimensional models that evolve the ICM composition assuming that the effects of magnetic fields can be parameterized or ignored. Such models for non-isothermal clusters show that Helium can sediment in the cluster core leading to a peak in concentration offset from the cluster center. The resulting profiles have recently been shown to be linearly unstable when the weakly-collisional character of the magnetized plasma is considered. In this paper, we present a modified version of the MHD code Athena, which makes it possible to evolve a weakly-collisional plasma subject to a gravitational field and stratified in both temperature and composition. We thoroughly test our implementation and confirm excellent agreement against several analytical results. In order to isolate the effects of composition, in this initial study we focus our attention on isothermal plasmas. We show that plasma instabilities, feeding off gradients in composition, can induce turbulent mixing and saturate by re-arranging magnetic field lines and alleviating the composition gradient. Composition profiles that increase with radius lead to instabilities that saturate by driving the average magnetic field inclination to roughly $45^{\circ}$. We speculate that this effect may alleviate the core insulation observed in homogeneous settings, with potential consequences for the associated cooling flow problem.

preprint2020arXiv

On Helium Mixing in Quasi-global Simulations of the Intracluster Medium

The assumption of a spatially uniform helium distribution in the intracluster medium can lead to biases in the estimates of key cluster parameters if composition gradients are present. The helium concentration profile in galaxy clusters is unfortunately not directly observable. Current models addressing the putative sedimentation are one-dimensional and parametrize the presence of magnetic fields in a crude way, ignoring the weakly-collisional, magnetized nature of the medium. When these effects are considered, a wide variety of instabilities can play an important role in the plasma dynamics. In a series of recent papers, we have developed the local, linear theory of these instabilities and addressed their non-linear development with a modified version of Athena. Here, we extend our study by developing a quasi-global approach that we use to simulate the mixing of helium as induced by generalizations of the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI), which feed off thermal and composition gradients. In the inner region of the ICM, mixing can occur on few Gyrs, after which the average magnetic field inclination angle is $\sim 45^{\circ}$ resulting in an averaged Spitzer parameter higher by about 20 % than the value obtained in homogeneous simulations. In the cluster outskirts the instabilities are rather inefficient, due to the shallow gradients. This suggests that compositions gradients in cluster cores might be shallower than one-dimensional models predict. More quantitative statements demand more refined models that can incorporate the physics driving the sedimentation process and simultaneously account for the weakly-collisional nature of the plasma.

preprint2020arXiv

The impact of magnetic fields on cold streams feeding galaxies

High redshift, massive halos are observed to have sustained, high star formation rates, which require that the amount of cold gas in the halo is continuously replenished. The cooling time scale for the hot virialized halo gas is too long to provide the source of cold gas. Supersonic, cold streams have been invoked as a mechanism for feeding massive halos at high redshift and deliver the cold gas required for continued star formation at the rates observed. This mechanism for replenishing the cold gas reservoir is motivated by some cosmological simulations. However, the cold streams are likely to be subject to the supersonic version of the Kelvin-Helmholtz instability (KHI), which eventually leads to stream disruption. Cosmological simulations have yet to obtain the spatial resolution required for understanding the detailed stability properties of cold streams. In this paper, we consider instead an idealized model of magnetized cold streams that we spatially resolve. Using linear theory we show how magnetic fields with dynamically important field strengths do not inhibit the KHI but rather enhance its growth rate. We perform nonlinear simulations of magnetized stream disruption and find that magnetic fields can nevertheless increase stream survival times by suppressing the mixing rate of cold gas with the circumgalactic medium. We find that magnetic fields can allow streams to survive $\sim 2-8$ times longer and, consequently, that streams $\sim 2-8$ times thinner can reach the central galaxy if the magnetic field strength is $\sim 0.3-0.8 μ$G.