![]() | This user is a student editor in Purdue_University/Astroparticle_Physics_(Spring_2023). |
Cosmological Simulation
editGeneral Introduction
editIncluding history, methods, and theory === will emphasize on hydrodynamics simulation
Hydrodynamics Simulation
editDark energy and dark matter account for most of the Universe's energy, so it is valid to ignore baryons(N-body simulation) when simulating large-scale structure formation. However, since the visible components of galaxies consist of baryons, it is crucial to include baryons in the simulation to study the detailed structures of galaxies. At first, the baryon component consists of mostly hydrogen and helium gas, which later transforms into stars during the formation of structures.
Euler equations
editIn cosmological simulations, astrophysical gases are typically modeled as inviscid ideal gases that follow the Euler equations, which can be expressed mainly in three different ways: Lagrangian, Eulerian, or arbitrary Lagrange-Eulerian methods. Different methods give specific forms of hydrodynamical equations [1]. When using the Lagrangian approach to specify the field, it is assumed that the observer tracks a specific fluid parcel with its unique characteristics during its movement through space and time. In contrast, the Eulerian approach emphasizes particular locations in space that the fluid passes through as time progresses.
Baryonic Physics
editTo shape the population of galaxies, the hydrodynamical equations must be supplemented by a variety of astrophysical processes mainly governed by baryonic physics.
Gas cooling
editProcesses, such as collisional excitation, ionization, and inverse Compton scattering, can cause the internal energy of the gas to be dissipated. In the simulation, cooling processes are realized by coupling cooling functions to energy equations. Besides the primordial cooling, at high temperature, , heavy elements (metals) cooling dominates [2]. When , the fine structure and molecular cooling also need to be considered to simulate the cold phase of the interstellar medium.
Interstellar medium
editComplex multi-phase structure, including relativistic particles and magnetic field, makes simulation of interstellar medium difficult. In particular, modeling the cold phase of the interstellar medium poses technical difficulties due to the short timescales associated with the dense gas. In the early simulations, the dense gas phase is frequently not modeled directly but rather characterized by an effective polytropic equation of state[3]. More recent simulations use a multimodal distribution[4][5] to describe the gas density and temperature distributions, which directly model the multi-phase structure. However, more detailed physics processes needed to be considered in future simulations, since the structure of the interstellar medium directly affects star formation.
Star formation
editAs cold and dense gas accumulates, it undergoes gravitational collapse and eventually forms stars. To simulate this process, a portion of the gas is transformed into collisionless star particles, which represent coeval, single-metallicity stellar populations and are described by an initial underlying mass function. Observations suggest that star formation efficiency in molecular gas is almost universal, with around 1% of the gas being converted into stars per free fall time[6]. In simulations, the gas is typically converted into star particles using a probabilistic sampling scheme based on the calculated star formation rate. Some simulations seek an alternative to the probabilistic sampling scheme and aim to better capture the clustered nature of star formation by treating star clusters as the fundamental unit of star formation. This approach permits the growth of star particles by accreting material from the surrounding medium[7]. In addition to this, modern models of galaxy formation track the evolution of these stars and the mass they return to the gas component, leading to an enrichment of the gas with metals[8].
Stellar feedback
editStars have an influence on their surrounding gas by injecting energy and momentum. This creates a feedback loop that regulates the process of star formation. To effectively control star formation, stellar feedback must generate galactic-scale outflows that expel gas from galaxies. Various methods are utilized to couple energy and momentum, particularly through supernova explosions, to the surrounding gas. These methods differ in how the energy is deposited, either thermally or kinetically. However, excessive radiative gas cooling must be avoided in the former case. Cooling is expected in dense and cold gas, but it cannot be reliably modeled in cosmological simulations due to low resolution. This leads to artificial and excessive cooling of the gas, causing the supernova feedback energy to be lost via radiation and significantly reducing its effectiveness. In the latter case, kinetic energy cannot be radiated away until it thermalizes. However, using hydrodynamically-decoupled wind particles to inject momentum non-locally into the gas surrounding active star-forming regions may still be necessary to achieve large-scale galactic outflows[9]. Recent models [10]explicitly model stellar feedback. These models not only incorporate supernova feedback but also consider other feedback channels such as energy and momentum injection from stellar winds, photoionization, and radiation pressure resulting from radiation emitted by young, massive stars[11].
Supermassive black holes
editSimulation of supermassive black holes is also considered, numerically seeding them in dark matter haloes, due to their observation in many galaxies[12] and the impact of their mass on the mass density distribution. Their mass accretion rate is frequently modeled by the Bondi-Hoyle model.
Active galactic nuclei
editActive galactic nuclei (AGN) have an impact on the observational phenomena of supermassive black holes, and further have a regulation of black hole growth and star formation. In simulations, AGN feedback is usually classified into two modes, namely quasar and radio mode. Quasar mode feedback is linked to the radiatively efficient mode of black hole growth and is frequently incorporated through energy or momentum injection[13]. The regulation of star formation in massive galaxies is believed to be significantly influenced by radio mode feedback, which occurs due to the presence of highly-collimated jets of relativistic particles. These jets are typically linked to X-ray bubbles that possess enough energy to counterbalance cooling losses[14].
Magnetic fields
editThe ideal magnetohydrodynamics approach is commonly utilized in cosmological simulations since it provides a good approximation for cosmological magnetic fields. The effect of magnetic fields on the dynamics of gas is generally negligible on large cosmological scales. Nevertheless, magnetic fields are a critical component of the interstellar medium since they provide pressure support against gravity[15] and affect the propagation of cosmic rays [16].
Cosmic rays
editCosmic rays play a significant role in the interstellar medium by contributing to its pressure [17], serving as a crucial heating channel [18], and potentially driving galactic gas outflows [19]. The propagation of cosmic rays is highly affected by magnetic fields. So in the simulation, equations describing the cosmic ray energy and flux are coupled to magnetohydrodynamics equations [20].
Radiation Hydrodynamics
editRadiation hydrodynamics simulations are computational methods used to study the interaction of radiation with matter. In astrophysical contexts, radiation hydrodynamics is used to study the epoch of reionization when the Universe had high redshift. There are several numerical methods used for radiation hydrodynamics simulations, including ray-tracing, Monte Carlo, and moment-based methods. Ray-tracing involves tracing the paths of individual photons through the simulation and computing their interactions with matter at each step. This method is computationally expensive but can produce very accurate results.
Results of Hydrodynamical Simulations
editThe simulation results can be compared with observations to test models used in simulations.
References
edit- ^ Gingold, R. A.; Monaghan, J. J. (1 December 1977). "Smoothed particle hydrodynamics: theory and application to non-spherical stars". Monthly Notices of the Royal Astronomical Society. 181 (3): 375–389. doi:10.1093/mnras/181.3.375. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Wiersma, Robert P. C.; Schaye, Joop; Smith, Britton D. (11 February 2009). "The effect of photoionization on the cooling rates of enriched, astrophysical plasmas". Monthly Notices of the Royal Astronomical Society. 393 (1): 99–107. doi:10.1111/j.1365-2966.2008.14191.x. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Springel, V.; Hernquist, L. (21 February 2003). "Cosmological smoothed particle hydrodynamics simulations: a hybrid multiphase model for star formation". Monthly Notices of the Royal Astronomical Society. 339 (2): 289–311. doi:10.1046/j.1365-8711.2003.06206.x. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Hopkins, Philip F.; Quataert, Eliot; Murray, Norman (16 March 2012). "The structure of the interstellar medium of star-forming galaxies". Monthly Notices of the Royal Astronomical Society. 421 (4): 3488–3521. doi:10.1111/j.1365-2966.2012.20578.x. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Agertz, Oscar; Kravtsov, Andrey V.; Leitner, Samuel N.; Gnedin, Nickolay Y. (21 May 2013). "TOWARD A COMPLETE ACCOUNTING OF ENERGY AND MOMENTUM FROM STELLAR FEEDBACK IN GALAXY FORMATION SIMULATIONS". The Astrophysical Journal. 770 (1): 25. doi:10.1088/0004-637X/770/1/25. eISSN 1538-4357. ISSN 0004-637X.
- ^ Bigiel, F.; Leroy, A. K.; Walter, F.; Brinks, E.; de Blok, W. J. G.; Kramer, C.; Rix, H. W.; Schruba, A.; Schuster, K.-F.; Usero, A.; Wiesemeyer, H. W. (3 March 2011). "A CONSTANT MOLECULAR GAS DEPLETION TIME IN NEARBY DISK GALAXIES". The Astrophysical Journal. 730 (2): L13. doi:10.1088/2041-8205/730/2/L13. eISSN 2041-8213. ISSN 2041-8205.
- ^ Li, Hui; Gnedin, Oleg Y.; Gnedin, Nickolay Y.; Meng, Xi; Semenov, Vadim A.; Kravtsov, Andrey V. (3 January 2017). "STAR CLUSTER FORMATION IN COSMOLOGICAL SIMULATIONS. I. PROPERTIES OF YOUNG CLUSTERS". The Astrophysical Journal. 834 (1): 69. doi:10.3847/1538-4357/834/1/69. eISSN 1538-4357.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Vogelsberger, Mark; Genel, Shy; Sijacki, Debora; Torrey, Paul; Springel, Volker; Hernquist, Lars (23 October 2013). "A model for cosmological simulations of galaxy formation physics". Monthly Notices of the Royal Astronomical Society. 436 (4): 3031–3067. doi:10.1093/mnras/stt1789. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Pillepich, Annalisa; Springel, Volker; Nelson, Dylan; Genel, Shy; Naiman, Jill; Pakmor, Rüdiger; Hernquist, Lars; Torrey, Paul; Vogelsberger, Mark; Weinberger, Rainer; Marinacci, Federico (12 October 2017). "Simulating galaxy formation with the IllustrisTNG model". Monthly Notices of the Royal Astronomical Society. 473 (3): 4077–4106. doi:10.1093/mnras/stx2656. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Muratov, Alexander L.; Kereš, Dušan; Faucher-Giguère, Claude-André; Hopkins, Philip F.; Quataert, Eliot; Murray, Norman (13 October 2015). "Gusty, gaseous flows of FIRE: galactic winds in cosmological simulations with explicit stellar feedback". Monthly Notices of the Royal Astronomical Society. 454 (3): 2691–2713. doi:10.1093/mnras/stv2126. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Stinson, G. S.; Brook, C.; Macciò, A. V.; Wadsley, J.; Quinn, T. R.; Couchman, H. M. P. (23 October 2012). "Making Galaxies In a Cosmological Context: the need for early stellar feedback". Monthly Notices of the Royal Astronomical Society. 428 (1): 129–140. doi:10.1093/mnras/sts028. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Moran, Edward C.; Shahinyan, Karlen; Sugarman, Hannah R.; Vélez, Darik O.; Eracleous, Michael (13 November 2014). "BLACK HOLES AT THE CENTERS OF NEARBY DWARF GALAXIES". The Astronomical Journal. 148 (6): 136. doi:10.1088/0004-6256/148/6/136. eISSN 1538-3881.
- ^ Springel, Volker; Di Matteo, Tiziana; Hernquist, Lars (August 2005). "Modelling feedback from stars and black holes in galaxy mergers". Monthly Notices of the Royal Astronomical Society. 361 (3): 776–794. doi:10.1111/j.1365-2966.2005.09238.x. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Sijacki, Debora; Springel, Volker; Di Matteo, Tiziana; Hernquist, Lars (20 August 2007). "A unified model for AGN feedback in cosmological simulations of structure formation". Monthly Notices of the Royal Astronomical Society. 380 (3): 877–900. doi:10.1111/j.1365-2966.2007.12153.x. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link) - ^ Ferrière, Katia M. (5 December 2001). "The interstellar environment of our galaxy". Reviews of Modern Physics. 73 (4): 1031–1066. doi:10.1103/RevModPhys.73.1031. eISSN 1539-0756. ISSN 0034-6861.
- ^ Kotera, Kumiko; Olinto, Angela V. (22 September 2011). "The Astrophysics of Ultrahigh-Energy Cosmic Rays". Annual Review of Astronomy and Astrophysics. 49 (1): 119–153. doi:10.1146/annurev-astro-081710-102620. eISSN 1545-4282. ISSN 0066-4146.
- ^ Cox, Donald P. (1 September 2005). "The Three-Phase Interstellar Medium Revisited". Annual Review of Astronomy and Astrophysics. 43 (1): 337–385. doi:10.1146/annurev.astro.43.072103.150615. eISSN 1545-4282. ISSN 0066-4146.
- ^ Wolfire, M. G.; Hollenbach, D.; McKee, C. F.; Tielens, A. G. G. M.; Bakes, E. L. O. (April 1995). "The neutral atomic phases of the interstellar medium". The Astrophysical Journal. 443: 152. doi:10.1086/175510. eISSN 1538-4357. ISSN 0004-637X.
- ^ Booth, C. M.; Agertz, Oscar; Kravtsov, Andrey V.; Gnedin, Nickolay Y. (18 October 2013). "SIMULATIONS OF DISK GALAXIES WITH COSMIC RAY DRIVEN GALACTIC WINDS". The Astrophysical Journal. 777 (1): L16. doi:10.1088/2041-8205/777/1/L16. eISSN 2041-8213. ISSN 2041-8205.
- ^ Thomas, T; Pfrommer, C (25 January 2019). "Cosmic-ray hydrodynamics: Alfvén-wave regulated transport of cosmic rays". Monthly Notices of the Royal Astronomical Society. 485 (3): 2977–3008. doi:10.1093/mnras/stz263. eISSN 1365-2966. ISSN 0035-8711.
{{cite journal}}
: CS1 maint: unflagged free DOI (link)