Master’s Thesis
Galaxy clusters assemble hierarchically through mergers whose shocks, turbulence, and core disruption leave rich thermal (X-ray) and non-thermal (radio) imprints. Inferring the underlying merger parameters—time since/to pericenter, collision velocity, mass ratio, pericenter distance, and component masses—is scientifically valuable but much of the relevant information is lost or distorted in observations. Moreover, observations provide only single-epoch snapshots rather than time-resolved tracks over Gyr, so reliable ground-truth labels are scarce. Hence we use TNG–Cluster simulations and its 352 zoom-in halos in three orthogonal projections across 0 ≤ z ≤ 1 to create intrinsic maps and their merger parameters. Building on these data, this thesis develops an end-to-end, simulation-based route from images to calibrated posteriors over merger physics. First, SimCLR encoder, a contrastive learning framework, are trained separately for X-ray, radio, and for paired X-ray+radio maps to learn morphology-aware representation that distill high-resolution maps into feature vectors that capture the relevant structure. A physics-aware augmentation suite (rotations, flips, Gaussian blur, additive noise, and affine zoom/shift; applied channel-consistently in the joint case) promotes invariance to observing nuisances without erasing merger signatures. Second, a conditional invertible neural network (cINN) with rational–quadratic spline couplings and a mixture-of-experts partition in representation space enabling inference of the p(x|c), where x is the merger parameter and c is the conditioner. Two conditioning modes are evaluated: (i) representation-conditioned, where c is the learned embedding (X-ray, radio, or joint); and (ii) scalar-conditioned, where c is a vector of scalar observable parameters that are fed directly to the cINN, without contrastive learning.
Across four conditioners schemes; scalars, X-ray embedding, radio embedding, and joint X-ray+radio, the method is evaluated on both last- and next-merger targets. Radio conditioning delivers the sharpest and best-calibrated posteriors, with typical MAP ranges of collision time ∼[−5, 10]%, velocity ≤ ±1%, main cluster's M500c, subcluster mass ∼ ±0.5%, and pericenter distance ∼ ±3%, reflecting the close coupling of radio emission and its geometry to recent cluster dynamics. X-ray conditioning remains informative, but produces broader posteriors and larger scatter in MAP estimates. Joint conditioning is consistently intermediate: it narrows and stabilizes inferences relative to X-ray alone, but it does not surpass radio conditioning. Scalar conditioning on the other hand, performs competitively only for collision time and the main cluster mass, while elsewhere it is not calibrated and suffers from heteroscedastic dispersions, pointing to the information bottle neck that could be introduced by only using scalar parameters. In all cases, the cINN not only recovers ΛCDM-consistent cross-parameter correlations, but also uses the learned correlations encoded in the joint density to propagate information from well constrained coordinated to weaker ones. Forecasts for the next merger mirror last-merger performance with modest, physically expected broadening (most visible for timing). Methodologically, this thesis demonstrates that label-free contrastive representation used directly as cINN condition enable survey-scale, uncertainty-aware inference of cluster assembly from imaging alone. Practically, the results argue for prioritizing high-fidelity radio imaging. In future research, this framework can be extended toward instrument-aware mock observables for sim-to-real transfer, multi-resolution fusion including SZ and weak-lensing, and scalable pipelines suitable for eROSITA/Chandra/XMM paired with LOFAR/MeerKAT/VLA.
Master's Projects
Potential-Density Pairs for Galaxy Discs with Exponential or Vertical Profile
The main results of this study are the analytic mass models for thin discs with exponential or sech2 vertical profile. These are obtained by modifying the razor-thin Kuzmin (1956) disc, very similar to how the Miyamoto & Nagai (1975; MN) model is constructed. In fact, many properties of the new models closely follow those of the MN model (with appropriately chosen parameters), so that the main difference is the vertical structure. Our approach can be used to construct modified Kuzmin models with arbitrary vertical profiles. The central density of an exponential vertical profile is significantly higher than for the MN or the sech2 models (at the same surface density and exponential scale height). This higher central density results in a stronger vertical force at small |z|, and therefore also in higher vertical orbital frequencies Ωz for stars with small maximal orbital height zmax. This in turn affects various dynamical effects, such as the presence or absence of orbital resonances (between vertical and radial motion) and the degree of phase-mixing in the vertical phase space. Thus, the modelling of phenomena related to the vertical structure of galactic discs, such as z–vz phase spirals and breathing or bending waves, are likely to be affected by the assumed vertical profile. Another result of this study are the cored-exponential profiles, for which we also give analytic mass models constructed by modifying the Kuzmin (1956) disc. These models share with the exponential and sech2 profiles the exponential decay with scale height h at large |z|, but differ at small |z|, where they possess a core of near-constant density with adjustable width w ≤ h (the sech2 profile has core width w = 2h). As the precise vertical profiles of galactic discs near z = 0 are difficult to assess observationally and no theoretical foundations for either the exponential or sech2 (or any other profile) exist yet, these cored-exponential profiles are a useful addition to the dynamicist’s toolbox and allow study of the effect of a density core of any width w ≤ h. Like the MN models, our new models are suitable bases for the construction of non-axisymmetric bar-shaped models by convolving them with a function f(x) (Long & Murali 1992) and in fact, we have added them to the suite of mass models provided by the discBar code (Dehnen & Aly 2023).
Inferring Halo Gas Fractions with Mock X-ray Observations of Cosmological Hydrodynamical Simulations
This paper investigates the baryonic gas content in galaxy clusters by using mock X-ray analysis based on the IllustrisTNG simulation and comparing the results to those from an observational study by Chiu et al. The primary focus is on estimating gas fractions within simulated halos and comparing them to observational trends. Following the methodology used by Chiu et al., we applied a modified beta-model to calculate gas fractions and compared the mock results with both simulation data and observational estimates. Our findings reveal that the mock data consistently overestimate the gas content compared to direct simulation data. Additionally, while Chiu et al.'s observational data show a trend of increasing gas fraction with halo mass, our mock data show no such trend, with gas fractions remaining relatively flat across different masses. This suggests that the feedback mechanisms within the IllustrisTNG simulation, particularly AGN feedback, are not strong enough to expel gas beyond the halo’s virial radius, leading to inflated gas fractions. These results indicate that current estimates of halo gas content, both in simulations and potentially in observational studies, may be overestimated. The research recommends future improvements, including stronger feedback models and more accurate X-ray luminosity equations, to better match observed gas distributions. Extending the analysis to the TNG-Cluster Simulation is also suggested to explore gas fractions in more massive halos.
Dynamics of Spiral Structure
This study investigates the mechanisms responsible for the formation and evolution of spiral structures in galaxies. Swing amplification is highlighted as a key process that converts leading density waves into trailing ones, with amplification occurring at the corotation radius, as supported by Toomre’s (1981) simulations. Tidal interactions are analyzed using N-body simulations from Semczuk et al. (2017), demonstrating how gravitational interactions with nearby satellite galaxies can generate grand-design spirals through repeated encounters. Bar-driven spirals are examined through hydrodynamical simulations from Huntley et al. (1978), showing the alignment of spiral structures with the rotation of bars, leading to the formation of stable, bisymmetric spiral arms. The results indicate that swing amplification, facilitated by differential rotation and self-gravity, effectively produces transient, multi-armed, and grand-design spiral patterns. Meanwhile, tidal interactions contribute to the persistence of grand-design spirals, and bar-driven mechanisms lead to well-defined, long-lived spiral structures. Observational evidence, such as radial mixing and gas distribution in spirals, supports the transient nature of these patterns, while simulations confirm the emergence of multi-armed spirals under certain conditions. The study emphasizes that swing amplification is the primary mechanism for forming transient spiral patterns, with tidal and bar-driven processes playing key roles in the structure of grand-design galaxies.
Bachelor’s Thesis
Black Hole Information Paradox
In this project, we studied the paradox of black hole information and a few possible solutions that have been suggested. We began with a brief look at the history of black holes and the introduction to some of the useful mathematical tools for studying them and then moved on to the similarities between the laws of classical and black holes’ thermodynamics. In the next chapter, we obtained Hawking radiation using the Unruh effect, and the reason for using this method was to demonstrate the dependence of the concepts of particles and vacuum on the observer. Next, information and its relationship with entropy were introduced. In the last section, we discussed the leading cause of the paradox, which is the vacuum without information in the event horizon predicted by general relativity. As the black holes evaporate, this vacuum stretches to produce a Hawking pair.
The first solution was the idea of remnants; a relatively old idea that the evaporation of a black hole will stop when it reaches Planck size. Following that, we analyzed Page’s average information and proved it mathematically. Moreover, we looked at Page’s theory which suggests that the information starts to come out only after the black hole has reached its Page time. Finally, we studied the complementarity theory and the interesting article on entanglement and wormholes. The proposed solutions each had contradictions: infinite malice in the remnant theory, loss of information in Page theory, and the existence of a firewall at the black hole’s event horizon in the complementarity approach, which conflicts with general relativity’s equivalence principle.