No Arabic abstract
In this study we present three-dimensional radiative cooling hydrodynamical simulations of galactic winds generated particularly in M82-like starburst galaxies. We have considered intermittent winds induced by SNe explosions within super star clusters randomly distributed in the central region of the galaxy and were able to reproduce the observed M82 wind conditions with its complex morphological outflow structure. We have found that the environmental conditions in the disk in nearly recent past are crucial to determine whether the wind will develop a large scale rich filamentary structure, as in M82 wind, or not. Also, the numerical evolution of the SN ejecta have allowed us to obtain the abundance distribution over the first 3 kpc extension of the wind and we have found that the SNe explosions change significantly the metallicity only of the hot, low-density wind component. Moreover, we have found that the SN-driven wind transports to outside the disk large amounts of energy, momentum and gas, but the more massive high-density component reaches only intermediate altitudes smaller than 1.5 kpc. Therefore, no significant amounts of gas mass are lost to the IGM and the mass evolution of the galaxy is not much affected by the starburst events occurring in the nuclear region.
We investigate the properties of the galaxies that reionized the Universe and the history of cosmic reionization using the Evolution and Assembly of GaLaxies and their environments (EAGLE) cosmological hydrodynamical simulations. We obtain the evolution of the escape fraction of ionizing photons in galaxies assuming that galactic winds create channels through which 20~percent of photons escape when the local surface density of star formation is greater than $0.1$ M$_odot$ yr$^{-1}$ kpc$^{-2}$. Such threshold behaviour for the generation of winds is observed, and the rare local objects which have such high star formation surface densities exhibit high escape fractions of $sim 10$ percent. In our model the luminosity-weighted mean escape fraction increases with redshift as $bar f_{rm esc}=0.045~((1+z)/4)^{1.1}$ at $z>3$, and the galaxy number weighted mean as $langle f_{rm esc} rangle=2.2times10^{-3}~((1+z)/4)^4$, and becomes constant $approx0.2$ at redshift $z>10$. The escape fraction evolves as an increasingly large fraction of stars forms above the critical surface density of star formation at earlier times. This evolution of the escape fraction, combined with that of the star formation rate density from EAGLE, reproduces the inferred evolution of the filling factor of ionized regions during the reionization epoch ($6<z<8$), the evolution of the post-reionization ($0leq z<6$) hydrogen photoionisation rate, and the optical depth due to Thomson scattering of the cosmic microwave background photons measured by the Planck satellite.
In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae.
We study galaxy super-winds driven in major mergers, using pc-resolution simulations with detailed models for stellar feedback that can self-consistently follow the formation/destruction of GMCs and generation of winds. The models include molecular cooling, star formation at high densities in GMCs, and gas recycling and feedback from SNe (I&II), stellar winds, and radiation pressure. We study mergers of systems from SMC-like dwarfs and Milky Way analogues to z~2 starburst disks. Multi-phase super-winds are generated in all passages, with outflow rates up to ~1000 M_sun/yr. However, the wind mass-loading efficiency (outflow rate divided by SFR) is similar to that in isolated galaxy counterparts of each merger: it depends more on global galaxy properties (mass, size, escape velocity) than on the dynamical state of the merger. Winds tend to be bi- or uni-polar, but multiple events build up complex morphologies with overlapping, differently-oriented bubbles/shells at a range of radii. The winds have complex velocity and phase structure, with material at a range of speeds up to ~1000 km/s, and a mix of molecular, ionized, and hot gas that depends on galaxy properties and different feedback mechanisms. These simulations resolve a problem in some sub-grid models, where simple wind prescriptions can dramatically suppress merger-induced starbursts. But despite large mass-loading factors (>~10) in the winds, the peak SFRs are comparable to those in no wind simulations. Wind acceleration does not act equally, so cold dense gas can still lose angular momentum and form stars, while blowing out gas that would not have participated in the starburst in the first place. Considerable wind material is not unbound, and falls back on the disk at later times post-merger, leading to higher post-starburst SFRs in the presence of stellar feedback. This may require AGN feedback to explain galaxy quenching.
The space-borne missions have provided us with a wealth of high-quality observational data that allows for seismic inferences of stellar interiors. This requires the computation of precise and accurate theoretical frequencies, but imperfect modeling of the uppermost stellar layers introduces systematic errors. To overcome this problem, an empirical correction has been introduced by Kjeldsen et al. (2008, ApJ, 683, L175) and is now commonly used for seismic inferences. Nevertheless, we still lack a physical justification allowing for the quantification of the surface-effect corrections. We used a grid of these simulations computed with the CO$^5$BOLD code to model the outer layers of solar-like stars. Upper layers of the corresponding 1D standard models were then replaced by the layers obtained from the horizontally averaged 3D models. The frequency differences between these patched models and the 1D standard models were then calculated using the adiabatic approximation and allowed us to constrain the Kjeldsen et al. power law, as well as a Lorentzian formulation. We find that the surface effects on modal frequencies depend significantly on both the effective temperature and the surface gravity. We further provide the variation in the parameters related to the surface-effect corrections using their power law as well as a Lorentzian formulation. Scaling relations between these parameters and the elevation (related to the Mach number) is also provided. The Lorentzian formulation is shown to be more robust for the whole frequency spectrum, while the power law is not suitable for the frequency shifts in the frequency range above $ u_{rm max}$.
NGC 1097 is a nearby barred spiral galaxy believed to be interacting with the elliptical galaxy NGC 1097A located to its northwest. It hosts a Seyfert 1 nucleus surrounded by a circumnuclear starburst ring. Two straight dust lanes connected to the ring extend almost continuously out to the bar. The other ends of the dust lanes attach to two main spiral arms. To provide a physical understanding of its structural and kinematical properties, two-dimensional hydrodynamical simulations have been carried out. Numerical calculations reveal that many features of the gas morphology and kinematics can be reproduced provided that the gas flow is governed by a gravitational potential associated with a slowly rotating strong bar. By including the self-gravity of the gas disk in our calculation, we have found the starburst ring to be gravitationally unstable which is consistent with the observation in citet{hsieh11}. Our simulations show that the gas inflow rate is 0.17 M$_sun$ yr$^{-1}$ into the region within the starburst ring even after its formation, leading to the coexistence of both a nuclear ring and a circumnuclear disk.