Filtern
Erscheinungsjahr
Dokumenttyp
- Wissenschaftlicher Artikel (30)
- Dissertation (1)
- Postprint (1)
Sprache
- Englisch (32)
Gehört zur Bibliographie
- ja (32)
Schlagworte
- methods: numerical (32) (entfernen)
Context. The High Energy Stereoscopic System Galactic plane survey (HGPS) is to date the most comprehensive census of Galactic gamma -ray sources at very high energies (VHE; 100 GeV <= E <= 100 TeV). As a consequence of the limited sensitivity of this survey, the 78 detected gamma -ray sources comprise only a small and biased subsample of the overall population. The larger part consists of currently unresolved sources, which contribute to large-scale diffuse emission to a still uncertain amount.Aims. We study the VHE gamma -ray source population in the Milky Way. For this purpose population-synthesis models are derived based on the distributions of source positions, extents, and luminosities.Methods. Several azimuth-symmetric and spiral-arm models are compared for spatial source distribution. The luminosity and radius function of the population are derived from the source properties of the HGPS data set and are corrected for the sensitivity bias of the HGPS. Based on these models, VHE source populations are simulated and the subsets of sources detectable according to the HGPS are compared with HGPS sources.Results. The power-law indices of luminosity and radius functions are determined to range between -1.6 and -1.9 for luminosity and -1.1 and -1.6 for radius. A two-arm spiral structure with central bar is discarded as spatial distribution of VHE sources, while azimuth-symmetric distributions and a distribution following a four-arm spiral structure without bar describe the HGPS data reasonably well. The total number of Galactic VHE sources is predicted to be in the range from 800 to 7000 with a total luminosity and flux of (1.6-6.3) x 10(36) ph s(-1) and (3-15) x 10(-10) ph cm(-2) s(-1), respectively.Conclusions. Depending on the model, the HGPS sample accounts for (68-87)% of the emission of the population in the scanned region. This suggests that unresolved sources represent a critical component of the diffuse emission measurable in the HGPS. With the foreseen jump in sensitivity of the Cherenkov Telescope Array, the number of detectable sources is predicted to increase by a factor between 5 and 9.
Gas flows in galaxy mergers
(2022)
In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magnetohydrodynamical zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated by turbulent pressure, with turbulent Mach numbers peaking at values of 1.6-3.3. This implies that bridges are dominated by supersonic turbulence, and are thus ideal candidates for studying the impact of extreme environments on star formation. We also find that gas accreted from the circumgalactic medium (CGM) during the merger significantly contributes (27-51 percent) to the star formation rate (SFR) at the time of coalescence and drives the subsequent reignition of star formation in the merger remnant. Indeed, 19-53 percent of the SFR at z = 0 originates from gas belonging to the CGM prior the merger. Finally, we investigate the origin of the metallicity-diluted gas at the centre of merging galaxies. We show that this gas is rapidly accreted on to the Galactic Centre with a time-scale much shorter than that of normal star-forming galaxies. This explains why coalescing galaxies are not well-captured by the fundamental metallicity relation.
Galactic winds exhibit a multiphase structure that consists of hot-diffuse and cold-dense phases. Here we present high-resolution idealized simulations of the interaction of a hot supersonic wind with a cold cloud with the moving-mesh code AREPO in setups with and without radiative cooling. We demonstrate that cooling causes clouds with sizes larger than the cooling length to fragment in 2D and 3D simulations. We confirm earlier 2D simulations by McCourt et al. (2018) and highlight differences of the shattering processes of 3D clouds that are exposed to a hot wind. The fragmentation process is quantified with a friends-of-friends analysis of shattered cloudlets and density power spectra. Those show that radiative cooling causes the power spectral index to gradually increase when the initial cloud radius is larger than the cooling length and with increasing time until the cloud is fully dissolved in the hot wind. A resolution of around 1 pc is required to reveal the effect of cooling-induced fragmentation of a 100 pc outflowing cloud. Thus, state-of-the-art cosmological zoom simulations of the circumgalactic medium fall short by orders of magnitudes from resolving this fragmentation process. This physics is, however, necessary to reliably model observed column densities and covering fractions of Lyman alpha haloes, high-velocity clouds, and broad-line regions of active galactic nuclei.
Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix and dissolve in the hot wind. In contrast, stripped cold gas from a large cloud mixes with the hot wind to intermediate temperatures, and then becomes thermally unstable and cools, causing a net accretion of hot gas to the cold tail. Using the magneto-hydrodynamical code AREPO, we perform cloud crushing simulations and test analytical criteria for the transition between the growth and destruction regimes to clarify a current debate in the literature. We find that the hot-wind cooling time sets the transition radius and not the cooling time of the mixed phase. Magnetic fields modify the wind-cloud interaction. Draping of wind magnetic field enhances the field upstream of the cloud, and fluid instabilities are suppressed by a turbulently magnetized wind beyond what is seen for a wind with a uniform magnetic field. We furthermore predict jellyfish galaxies to have ordered magnetic fields aligned with their tails. We finally discuss how the results of idealized simulations can be used to provide input to subgrid models in cosmological (magneto-)hydrodynamical simulations, which cannot resolve the detailed small-scale structure of cold gas clouds in the circumgalactic medium.
The observation of the non-Keplerian behavior of propeller structures in Saturn's outer A ring raises the question: how does the propeller respond to the wandering of the central embedded moonlet? Here, we study numerically how the structural imprint of the propeller changes for a libration of the moonlet. It turns out that the libration induces an asymmetry in the propeller, which depends on the libration period and amplitude of the moonlet. Further, we study the dependence of the asymmetry on the libration period and amplitude for a moonlet with a 400 m Hill radius, which is located in the outer A ring. This allows us to apply our findings to the largest known propeller Blériot, which is expected to be of a similar size. For Blériot, we can conclude that, supposing the moonlet is librating with the largest observed period of 11.1 yr and an azimuthal amplitude of about 1845 km, a small asymmetry should be measurable but depends on the moonlet's libration phase at the observation time. The longitude residuals of other trans-Encke propellers (e.g., Earhart) show amplitudes similar to Blériot, which might allow us to observe larger asymmetries due to their smaller azimuthal extent, allowing us to scan the whole gap structure for asymmetries in one observation. Although the librational model of the moonlet is a simplification, our results are a first step toward the development of a consistent model for the description of the formation of asymmetric propellers caused by a freely moving moonlet.
Context.
Even after the Rosetta mission, some of the mechanical parameters of comet 67P/Churyumov-Gerasimenko's surface material are still not well constrained. They are needed to improve our understanding of cometary activity or for planning sample return procedures.
Aims.
We discuss the physical process dominating the formation of aeolian-like surface features in the form of moats and wind taillike bedforms around obstacles and investigate the mechanical and geometrical parameters involved.
Methods.
By applying the discrete element method (DEM) in a low-gravity environment, we numerically simulated the dynamics of the surface layer particles and the particle stream involved in the formation of aeolian-like morphological features. The material is composed of polydisperse spherical particles that consist of a mixture of dust and water ice, with interparticle forces given by the Hertz contact model, cohesion, friction, and rolling friction. We determined a working set of parameters that enables simulations to be reasonably realistic and investigated morphological changes when modifying these parameters.
Results.
The aeolian-like surface features are reasonably well reproduced using model materials with a tensile strength on the order of 0.1-1 Pa. Stronger materials and obstacles with round shapes impede the formation of a moat and a wind tail. The integrated dust flux required for the formation of moats and wind tails is on the order of 100 kg m(-2), which, based on the timescale of morphological changes inferred from Rosetta images, translates to a near-surface particle density on the order of 10(-6)-10(-4) kg m(-3).
Conclusions.
DEM modeling of the aeolian-like surface features reveals complex formation mechanisms that involve both deposition of ejected material and surface erosion. More numerical work and additional in situ measurements or sample return missions are needed to better investigate mechanical parameters of cometary surface material and to understand the mechanics of cometary activity.
We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at z = 0 for three different phases: 'cold', 'hot' and H i, and compare our results with observations. The amount of material in the hot halo (M-hot a parts per thousand 4-5 x 10(10) M-aS (TM)), and the cold (M-cold(r a parts per thousand(2) 10 kpc) a parts per thousand 10(8) M-aS (TM)) and H i components displays reasonable agreement with observations. We also compute the accretion/ejection rates together with the H i (radial and all-sky) covering fractions. The integrated H i accretion rate within r = 50 kpc gives similar to 0.2-0.3 M-aS (TM) yr(-1), i.e. close to that obtained from high-velocity clouds in the Milky Way. We find that the global accretion rate is dominated by hot material, although ionized gas with T a parts per thousand(2) 10(5) K can contribute significantly too. The net accretion rates of all material at the virial radii are 6-8 M-aS (TM) yr(-1). At z = 0, we find a significant gas excess between the two galaxies, as compared to any other direction, resulting from the overlap of their gaseous haloes. In our simulation, the gas excess first occurs at z similar to 1, as a result of the kinematical evolution of the LG.
We present results of 2D3V particle-in-cell simulations of nonrelativistic plasma collisions with absent or parallel large-scale magnetic field for parameters applicable to the conditions at young supernova remnants. We study the collision of plasma slabs of different density, leading to two different shocks and a contact discontinuity. Electron dynamics play an important role in the development of the system. While nonrelativistic shocks in both unmagnetized and magnetized plasmas can be mediated by Weibel-type instabilities, the efficiency of shock-formation processes is higher when a large-scale magnetic field is present. The electron distributions downstream of the forward and reverse shocks are generally isotropic, whereas that is not always the case for the ions. We do not see any significant evidence of pre-acceleration, neither in the electron population nor in the ion distribution.
We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the timescale of magnetic-field growth depends on the shock strength.
We have investigated via 2D relativistic magnetohydrodynamic simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the post-shock region, magnetic field is strongly amplified by turbulent motions triggered by pre-shock density inhomogeneities. Using a long-simulation box we have followed the magnetic field amplification until it is fully developed and saturated. The turbulent velocity is subrelativistic even for a strong shock. Magnetic field amplification is controlled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic field amplification and saturation depend on the initial strength and direction of the magnetic field in the pre-shock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the post-shock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the post-shock region. If the initial magnetic field in the pre-shock medium is strong, the post-shock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the post-shock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.