Understanding how galaxies form, interact, and evolve comes largely from comparing theory predictions with observational data. Numerical simulations of galaxies provide the most accurate approach... Show moreUnderstanding how galaxies form, interact, and evolve comes largely from comparing theory predictions with observational data. Numerical simulations of galaxies provide the most accurate approach to testing the theory, as they follow the non-linear evolution of gas and dark matter in great detail and incorporate numerous baryonic processes, among which are energy feedback from supernovae (SNe) and Active Galactic Nuclei (AGN). In this thesis, we show the results of the development of the new model COLIBRE for cosmological simulations of galaxy formation that include a cold interstellar medium. First, we present a new SN feedback recipe developed for COLIBRE, whereby SN energy is injected into the gas in thermal and kinetic forms, and the total energy and momentum of the system of gas and stars are exactly conserved. Second, we conduct a detailed comparison of different ways in which SN energy is distributed in the gas environment around young stellar populations. Third, by using our simulation setup originally developed to test COLIBRE’s SN feedback, we show that the radioactive isotope Fe60 that has been detected on Earth is likely of SN origin. Finally, we present the calibration of the SN and AGN feedback of the COLIBRE model using machine learning. Show less
In this thesis, we examine various systems through the lens of several numerical methods. We delve into questions concerning thermalization in closed unitary systems, lattice gauge theories, and... Show moreIn this thesis, we examine various systems through the lens of several numerical methods. We delve into questions concerning thermalization in closed unitary systems, lattice gauge theories, and the intriguing properties of deep neural network phase spaces. Leveraging modern advancements in both software and hardware, we scrutinize these systems in greater detail, accessing previously unreachable regimes. Show less
Circumstellar discs are the reservoirs of gas and dust that surround young stars and have the potential to become planetary systems. Their evolution will determine the time and material available... Show moreCircumstellar discs are the reservoirs of gas and dust that surround young stars and have the potential to become planetary systems. Their evolution will determine the time and material available to form planets. Studying the evolution of circumstellar discs can then help us understand planet formation and the diversity of observed planetary systems. These discs develop almost immediately after star formation, as a direct consequence of the collapse of a molecular cloud and angular momentum conservation. Their surroundings are rich in gas and neighbouring stars, which can be hostile to the discs and affect their evolution in different ways: dynamical encounters with nearby stars can truncate the discs; stellar winds and supernovae explosions can truncate, tilt, or completely destroy the discs; and the presence of bright, massive stars in the vicinity of circumstellar discs can heat their surface enough to evaporate mass from them. This process, known as external photoevaporation, is arguably one of the most important environmental mechanisms in depleting mass from young circumstellar discs. The work performed for this thesis consisted of simulating the early evolution of circumstellar discs in star clusters and the effects of the environment, in particular, truncations due to close encounters and photoevaporation. The results show that photoevaporation is extremely efficient in removing mass from the discs, greatly limiting the amount of material and time available to form planets. Show less
Ever since Isaac Newton in 1687 posed the N-body problem, astronomers have been looking for its solutions in order to understand the evolution of dynamical systems, such as our own solar system,... Show moreEver since Isaac Newton in 1687 posed the N-body problem, astronomers have been looking for its solutions in order to understand the evolution of dynamical systems, such as our own solar system, star clusters and galaxies. The main difficulty is that small errors grow exponentially, so that numerical solutions diverge easily from the mathematical solution. This thesis presents two new state of the art N-body algorithms, one of which is designed for high precision (Brutus) and the other for speed (Sakura). The assumption that N-body results are accurate in a statistical sense, is put to the test for three-body configurations. Finally, a new mathematical model is constructed that describes the origin of chaos in a dynamical systems, and explains the short Liapounov time of Comet Halley's orbit. Show less
The appearance of the first luminous sources in the Universe marks the beginning of the epoch of reionisation. The cosmic reionisation of hydrogen has been the last major phase transition in the... Show moreThe appearance of the first luminous sources in the Universe marks the beginning of the epoch of reionisation. The cosmic reionisation of hydrogen has been the last major phase transition in the evolution of the Universe, turning cool neutral gas into the hot, ionised plasma we observe today. Up to now, direct observations of reionisation have been impossible, most of our current knowledge therefore stems from numerical studies. However, even with current computing facilities the modelling of the reionisation epoch presents some serious challenges. Radiative transfer is perhaps the biggest challenge, due to the high dimensionality and non-locality of the problem. In this thesis the latest improvements of the SimpleX method, that was designed to overcome some of the problems related to radiative transfer, are presented. We have used this method to study the escape fraction of ionising photons from galaxies. Our results show that high-redshift dwarf galaxies can be important contributors to cosmic reionisation. In addition, we have studied the effect of the inclusion of diffuse recombination radiation in reionisation simulations on large scales, which has been neglected by recent studies. We show that diffuse recombination radiation does not have a significant influence on the reionisation process. Show less
Planets form in disks that are commonly found around young stars. The intimate relationship that exists between planet and disk can account for a lot of the exotic extrasolar planetary systems... Show morePlanets form in disks that are commonly found around young stars. The intimate relationship that exists between planet and disk can account for a lot of the exotic extrasolar planetary systems known today. In this thesis we explore disk-planet interaction using numerical hydrodynamical simulations. We study the growth and migration of embedded planets, as well as the condition for gap formation in the disk. These planetary gaps provide an important link to future observations of circumstellar disks. Show less
