Galactic Winds Research Articles

Warped diffusive radio halo around the quiescent spiral edge-on galaxy NGC 4565

Context. Cosmic rays play a pivotal role in launching galactic winds, particularly in quiescently star-forming galaxies where the hot gas alone is not sufficient to drive a wind. Except for the Milky Way, not much is known about the transport of cosmic rays in galaxies. Aims. In this Letter, we present low-frequency observations of the nearby edge-on spiral galaxy NGC 4565 using the LOw-Frequency ARray (LOFAR). With our deep 144 MHz observations, we obtain a clean estimate of the emission originating from old cosmic-ray electrons (CRe), which is almost free from contamination by thermal emission. Methods. We measured vertical profiles of the non-thermal radio continuum emission that we fitted with Gaussian and exponential functions. The different profile shapes correspond to 1D cosmic-ray transport models of pure diffusion and advection, respectively. Results. We detect a warp in the radio continuum that is reminiscent of the previously known H I warp. Because the warp is not seen at GHz-frequencies in the radio continuum, its minimum age must be about 100 Myr. The warp also explains the slight flaring of the thick radio disc that can otherwise be well described by a Gaussian profile with an FWHM of 65 arcsec (3.7 kpc). Conclusions. The diffusive radio halo together with the extra-planar X-ray emission may be remnants of enhanced star-forming activity in the past where the galaxy had a galactic wind, as GHz-observations indicate only a weak outflow in the last 40 Myr. NGC 4565 could be in transition from an outflow- to an inflow-dominated phase.

Probing gaseous halos of galaxies with radio jets

Context.Gaseous halos play a key role in understanding inflow, feedback, and the overall baryon budget in galaxies. Literature models predict transitions of the state of the gaseous halo between cold and hot accretion, winds, fountains, and hydrostatic halos at certain galaxy masses. Since luminosities of radio AGN are sensitive to halo densities, any significant transition would be expected to show up in the radio luminosities of large samples of galaxies. The LOw Frequency ARray (LOFAR) Two-Metre Sky Survey (LoTSS) has identified a galaxy stellar mass scale, 1011 M⊙, above which the radio luminosities increase disproportionately.Aims.We investigate if radio luminosities of galaxies, especially the marked rise at galaxy masses around 1011 M⊙, can be explained with standard assumptions regarding jet powers, scaling between black hole mass and galaxy mass, and gaseous halos.Methods.Based on observational data and theoretical constraints, we developed models for the radio luminosity of radio AGN in halos under infall, galactic wind, and hydrostatic conditions. We compared these models to LoTSS data for a large sample of galaxies in the mass range between 108.5 M⊙and 1012 M⊙.Results.Under the assumption that the same characteristic upper limit to jet powers known from high galaxy masses holds at all masses, we find the maximum radio luminosities for the hydrostatic gas halos to lie close to the upper envelope of the distribution of the LOFAR data. The marked rise in radio luminosity at 1011 M⊙is matched in our model and is related to a significant change in halo gas density around this galaxy mass, which is a consequence of lower cooling rates at a higher virial temperature. Wind and infall models overpredict the radio luminosities for small galaxy masses and have no particular steepening of the run of the radio luminosities predicted at any galaxy mass.Conclusions.Radio AGN could have the same characteristic Eddington-scaled upper limit to jet powers in galaxies of all masses in the sample if the galaxies have hydrostatic gas halos in phases when radio AGN are active. We find no evidence of a change of the type of galaxy halo with the galaxy mass. Galactic winds and quasi-spherical cosmological inflow phases cannot frequently occur at the same time as powerful jet episodes unless the jet properties in these phases are significantly different from what we assumed in our model.

Simulating the effect of photoheating feedback during reionization

Abstract We present self-consistent radiation hydrodynamic simulations of hydrogen reionization performed with arepo-rt complemented by a state-of-the-art galaxy formation model. We examine how photoheating feedback, due to reionization, shapes the galaxies properties. Our fiducial model completes reionization by z ≈ 6 and matches observations of the Ly α forest, the cosmic microwave background electron scattering optical depth, the high-redshift ultraviolet (UV) luminosity function, and stellar mass function. Contrary to previous works, photoheating suppresses star formation rates by more than $50{{\ \rm per\ cent}}$ only in haloes less massive than ∼108.4 M⊙ (∼108.8 M⊙) at z = 6 (z = 5), suggesting inefficient photoheating feedback from photons within galaxies. The use of a uniform UV background that heats up the gas at z ≈ 10.7 generates an earlier onset of suppression of star formation compared to our fiducial model. This discrepancy can be mitigated by adopting a UV background model with a more realistic reionization history. In the absence of stellar feedback, photoheating alone is only able to quench haloes less massive than ∼109 M⊙ at z ≳ 5, implying that photoheating feedback is sub-dominant in regulating star formation. In addition, stellar feedback, implemented as a non-local galactic wind scheme in the simulations, weakens the strength of photoheating feedback by reducing the amount of stellar sources. Most importantly, photoheating does not leave observable imprints in the UV luminosity function, stellar mass function, or the cosmic star formation rate density. The feasibility of using these observables to detect imprints of reionization therefore requires further investigation.

Kinematics of Circumgalactic Gas: Feeding Galaxies and Feedback

Abstract We present observations of 50 pairs of redshift z ≈ 0.2 star-forming galaxies and background quasars. These sightlines probe the circumgalactic medium (CGM) out to half the virial radius, and we describe the circumgalactic gas kinematics relative to the reference frame defined by the galactic disks. We detect halo gas in Mg ii absorption, measure the equivalent-width-weighted Doppler shifts relative to each galaxy, and find that the CGM has a component of angular momentum that is aligned with the galactic disk. No net counter-rotation of the CGM is detected within 45° of the major axis at any impact parameter. The velocity offset of the circumgalactic gas correlates with the projected rotation speed in the disk plane out to disk radii of roughly 70 kpc. We confirm previous claims that the Mg ii absorption becomes stronger near the galactic minor axis, and we show that the equivalent width correlates with the velocity range of the absorption. We cannot directly measure the location of any absorber along the sightline, but we explore the hypothesis that individual velocity components can be associated with gas orbiting in the disk plane or flowing radially outward in a conical outflow. We conclude that centrifugal forces partially support the low-ionization gas and galactic outflows kinematically disturb the CGM producing excess absorption. Our results firmly rule out schema for the inner CGM that lack rotation and suggest that angular momentum as well as galactic winds should be included in any viable model for the low-redshift CGM.

IGM‐Vis: Analyzing Intergalactic and Circumgalactic Medium Absorption Using Quasar Sightlines in a Cosmic Web Context

AbstractWe introduce IGM‐Vis, a novel astrophysics visualization and data analysis application for investigating galaxies and the gas that surrounds them in context with their larger scale environment, the Cosmic Web. Environment is an important factor in the evolution of galaxies from actively forming stars to quiescent states with little, if any, discernible star formation activity. The gaseous halos of galaxies (the circumgalactic medium, or CGM) play a critical role in their evolution, because the gas necessary to fuel star formation and any gas expelled from widely observed galactic winds must encounter this interface region between galaxies and the intergalactic medium (IGM). We present a taxonomy of tasks typically employed in IGM/CGM studies informed by a survey of astrophysicists at various career levels, and demonstrate how these tasks are facilitated via the use of our visualization software. Finally, we evaluate the effectiveness of IGM‐Vis through two in‐depth use cases that depict real‐world analysis sessions that use IGM/CGM data.

Large reservoirs of turbulent diffuse gas around high-z starburst galaxies

AbstractStarburst galaxies at z ∼ 2 – 4 are among the most intensely star-forming galaxies in the universe. The way they accrete their gas to form stars at such high rates is still a controversial issue. ALMA has detected the CH+ (J = 1-0) line in emission and/or absorption in all the gravitationally lensed starburst galaxies targeted so far at z ∼ 3. Its unique spectroscopic and chemical properties enable CH+ to highlight the sites of most intense dissipation of mechanical energy. The absorption lines reveal highly turbulent, massive reservoirs of low-density molecular gas. The broad emission lines, arising in myriad UV-irradiated molecular shocks, reveal powerful galactic winds. The CH+ lines therefore probe the fate of prodigious energy releases, due to infall and/or outflows, and primarily stored in turbulence before being radiated by cool molecular gas. The turbulent reservoirs act as mass and energy buffers over the duration of the starburst phase.

Magnetic Fields and Halos in Spiral Galaxies

Radio continuum and polarization observations reveal best the magnetic field structure and strength in nearby spiral galaxies. They show a similar magnetic field pattern, which is of spiral shape along the disk plane and X-shaped in the halo, sometimes accompanied by strong vertical fields above and below the central region of the disk. The strength of the total halo field is comparable to that of the disk. The small- and large-scale dynamo action is discussed to explain the observations with special emphasis on the rôle of star formation on the α − Ω dynamo and the magnetic field strength and structure in the disk and halo. Recently, with RM-synthesis of the CHANG-ES observations, we obtained the first observational evidence for the existence of regular magnetic fields in the halo. The analysis of the radio scale heights indicate escape-dominated radio halos with convective cosmic ray propagation for many galaxies. These galactic winds may be essential for an effective dynamo action and may transport large-scale magnetic field from the disk into the halo.

Mufasa: Time-scales for H i consumption and SFR depletion of satellite galaxies in groups

We investigate the connection between the HI~content, SFR and environment of galaxies using a hydrodynamic simulation that incorporates scaling relations for galactic wind and a heuristic halo mass-based quenching prescription. We run two zoom-in simulations of galaxy groups with $M_{halo}>10^{13}M_\odot$ at z=0, selected to have quiet merger histories. We track galaxies as they become satellites, and compute the delay time $\tau_{d}$ during which the satellites are similar to central galaxies at a given stellar mass, and a fading time $\tau_{f}$ during which satellites go from gas-rich and star-forming to gas-poor and quiescent. We find $0.7< \tau_{d}< 3$ Gyr at $z=0$, and depends inversely on the satellite halo mass at infall. At z ~ 1 we find ~$0.3< \tau_{d} < 2$ Gyr, broadly consistent with a positive correlation with the Hubble time. For a given halo mass, lower stellar mass galaxies at infall time have higher $\tau_{d}$. We generally find $\tau_{f}\ll \tau_{d}$, ranging between ~ 150 Myr at z~0 and ~ 80 Myr at z ~ 1 based on linear interpolation, with some uncertainty because they are smaller than our simulation output frequency ($200-300$ Myr). $\tau_{f}$ has no obvious dependency on infall halo mass. Both timescales show little difference between HI depletion and SF quenching, indicating that using up the gas reservoir by star formation without refilling is the main mechanism to transform satellite galaxies at these halo masses. At a given physical distance from the center of the main halo of interest, higher redshift galaxies have on average higher cold gas content, but the ratio of gas (HI or H$_2$) to star formation rate is similar, indicating that star formation is consistently fed through reservoirs of HI then H$_2$. For a given amount of HI, galaxies have shorter consumption times in more massive halo structures.

On the dynamics and survival of fractal clouds in galactic winds

Abstract Recent observations suggest that dense gas clouds can survive even in hot galactic winds. Here we show that the inclusion of turbulent densities with different statistical properties has significant effects on the evolution of wind-swept clouds. We investigate how the initial standard deviation of the lognormal density field influences the dynamics of quasi-isothermal clouds embedded in supersonic winds. We compare uniform, fractal solenoidal, and fractal compressive cloud models in both 3D and 2D hydrodynamical simulations. We find that the processes of cloud disruption and dense gas entrainment are functions of the initial density distribution in the cloud. Fractal clouds accelerate, mix, and are disrupted earlier than uniform clouds. Within the fractal cloud sample, compressive clouds retain high-density nuclei, so they are more confined, less accelerated, and have lower velocity dispersions than their solenoidal counterparts. Compressive clouds are also less prone to Kelvin–Helmholtz and Rayleigh–Taylor instabilities, so they survive longer than solenoidal clouds. By comparing the cloud properties at the destruction time, we find that dense gas entrainment is more effective in uniform clouds than in either of the fractal clouds, and it is more effective in solenoidal than in compressive models. In contrast, mass loading into the wind is more efficient in compressive cloud models than in uniform or solenoidal models. Overall, wide density distributions lead to inefficient entrainment, but they facilitate mass loading and favour the survival of very dense gas in hot galactic winds.

The insignificance of Seyfert 2 activity in driving cold-gas galactic winds

We study the interstellar Na I $\lambda \lambda 5890, 5895$ (Na D) absorption-line doublet in a nearly-complete sample of $\sim$9900 nearby Seyfert 2 galaxies, in order to quantify the significance of optical AGN activity in driving kpc-scale outflows that can quench star formation. Comparison to a carefully matched sample of $\sim$44,000 control objects indicates that the Seyfert and control population have similar Na D detection rates ($\sim 5-6%$). Only 53 Seyferts (or 0.5% of the population) are found to potentially display galactic-scale winds, compared to 0.8% of the control galaxies. While nearly a third of the Na D outflows observed in our Seyfert 2 galaxies occur around the brightest AGN, both radio and infrared data indicate that star formation could play the dominant role in driving cold-gas outflows in an even higher fraction of the Na D-outflowing Seyfert 2s. Our results indicate that galactic-scale outflows at low redshift are no more frequent in Seyferts than they are in their non-active counterparts, that optical AGN are not significant contributors to the quenching of star formation in the nearby Universe, and that star-formation may actually be the principal driver of outflows even in systems that do host an AGN.

Interpreting the Relation between the Gamma-Ray and Infrared Luminosities of Star-forming Galaxies

Abstract It has been found that there is a quasi-linear scaling relationship between the gamma-ray luminosity in GeV energies and the total infrared luminosity of star-forming galaxies, i.e., with α ≃ 1. However, the origin of this linear slope is not well understood. Although extreme starburst galaxies can be regarded as calorimeters for hadronic cosmic-ray interaction and thus a quasi-linear scaling may hold, it may not be the case for low star formation rate (SFR) galaxies, as the majority of cosmic rays in these galaxies are expected to escape. We calculate the gamma-ray production efficiency in star-forming galaxies by considering realistic galaxy properties, such as the gas density and galactic wind velocity in star-forming galaxies. We find that the slope for the relation between gamma-ray luminosity and the infrared luminosity gets steeper for low infrared luminosity galaxies, i.e., , due to increasingly lower efficiency for the production of gamma-ray emission. We further find that the measured data of the gamma-ray luminosity is compatible with such a steepening. The steepening in the slope suggests that cosmic-ray escape is very important in low-SFR galaxies.

Cosmic magnetism evolution using cosmological simulations

AbstractThe Intergalactic Medium (IGM) is the region comprising the environment between the galaxies. Gamma-ray observations have provided lower limits to IGM magnetic fields of the order of ≳10–16 G. Magnetic fields are continuously ejected from galaxies by jets and galactic winds. However, the origin and evolution of cosmic magnetic fields in the more diffuse regions, like voids, is still debated. The difficulties in directly measuring magnetic fields and their coherent scales, make hydrodynamic and magnetohydrodynamic (MHD) cosmological simulations useful tools to shed light on this debate. As a first approach, we have performed hydrodynamic cosmological simulations assuming energy equipartition as an initial condition between the baryonic gas and the magnetic field, starting at z = 8, to track the evolution of magnetic fields, and compare with results of MHD simulations. We have found that for halos and cores, our results are comparable to the MHD description. For the less dense regions, the equipartition condition clearly overestimates the observed limits. In forthcoming work, we will investigate MHD simulations of cosmological evolution and amplification of seed magnetic fields, considering all relevant feedback processes and exploring turbulent dynamo amplification versus primordial mechanisms across cosmological timescales.

CHANG-ES

Context. The observation of total and linearly polarized synchrotron radiation of spiral galaxies in the radio continuum reveals the distribution and structure of their magnetic fields. By observing these, information about the proposed dynamo processes that preserve the large-scale magnetic fields in spiral galaxies can be gained. Additionally, by analyzing the synchrotron intensity, the transport processes of cosmic rays into the halo of edge-on spiral galaxies can be investigated. Aims. We analyze the magnetic field geometry and the transport processes of the cosmic rays of the edge-on spiral starburst galaxy NGC 4666 from CHANG-ES radio data in two frequencies; 6 GHz (C-band) and 1.5 GHz (L-band). Supplementary X-ray data are used to investigate the hot gas in NGC 4666. Methods. We determine the radio scale heights of total power emission at both frequencies for this galaxy. We show the magnetic field orientations derived from the polarization data. Using rotation measure (RM) synthesis we further study the behavior of the RM values along the disk in C-band to investigate the large-scale magnetic-field pattern. We use the revised equipartition formula to calculate a map of the magnetic field strength. Furthermore, we model the processes of cosmic-ray transport into the halo with the 1D SPINNAKER model. Results. The extended radio halo of NGC 4666 is box-shaped and is probably produced by the previously observed supernova-driven superwind. This is supported by our finding of an advective cosmic-ray transport such as that expected for a galactic wind. The scale-height analysis revealed an asymmetric halo above and below the disk as well as between the two sides of the major axis. A central point source as well as a bubble structure is seen in the radio data for the first time. Our X-ray data show a box-shaped hot halo around NGC 4666 and furthermore confirm the AGN nature of the central source. NGC 4666 has a large-scale X-shaped magnetic field in the halo, as has been observed in other edge-on galaxies. The analysis furthermore revealed that the disk of NGC 4666 shows hints of field reversals along its radius, which is the first detection of this phenomenon in an external galaxy.

NIHAO XIX: how supernova feedback shapes the galaxy baryon cycle

We have used the NIHAO simulations to explore how supernovae (SNe) affect star formation in galaxies. We find that SN feedback operates on all scales from the interstellar medium (ISM) to several virial radii. SNe regulate star formation by preventing condensation of HI into H$_2$ and by moving cold neutral gas to the hot HII phase. The first effect explains why the cold neutral gas in dwarf galaxies forms stars inefficiently. The second maintains the hot ISM of massive galaxies (HII vents out at lower masses). At $v_{\rm vir}\simeq 67{\rm\,km\,s}^{-1}$, the outflow rate follows the relation: $\dot{M}_{\rm out}=23\,(v_{\rm vir}/67{\rm\,km\,s}^{-1})^{-4.6}\,{\rm SFR}$. $20\%$ to $70\%$ of the gas expelled from galaxies escapes from the halo (ejective feedback) but outflows are dominated by cold swept-up gas, most of which falls back onto the galaxy on a $\sim 1\,$Gyr timescale. This `fountain feedback' reduces the masses of galaxies by a factor of two to four, since gas spends half to three quarter of its time in the fountain. Less than $10\%$ of the ejected gas mixes with the hot circumgalactic medium and this gas is usually not reaccreted. On scales as large as $6r_{\rm vir}$, galactic winds divert the incoming gas from cosmic filaments and prevent if from accreting onto galaxies (pre-emptive feedback). This process is the main reason for the low baryon content of ultradwarves.

Correlation between SFR Surface Density and Thermal Pressure of Ionized Gas in Local Analogs of High-redshift Galaxies

Abstract We explore the relation between the star formation rate surface density (ΣSFR) and the interstellar gas pressure for nearby compact starburst galaxies. The sample consists of 17 green peas and 19 Lyman break analogs (LBAs). Green peas are nearby analogs of Lyα emitters at high redshift and LBAs are nearby analogs of Lyman break galaxies at high redshift. We measure the sizes of green peas using Hubble Space Telescope Cosmic Origins Spectrograph near-UV images with a spatial resolution of ∼0.″05. We estimate the gas thermal pressure in H ii regions by P = N total Tk B ≃ 2n e Tk B . The electron density is derived using the [S ii] doublet at 6716,6731 Å and the temperature is calculated from the [O iii] lines. The correlation is characterized by ΣSFR = . Green peas and LBAs have high ΣSFR up to 1.2 M ⊙ yr−1 kpc−2 and high thermal pressure in the H ii region up to P/k B ∼ 107.2 K cm−3. These values are at the highest end of the range seen in nearby starburst galaxies. The high gas pressure and the correlation are in agreement with those found in star-forming galaxies at z ∼ 2.5. These extreme pressures are shown to be responsible for driving galactic winds in nearby starbursts. These outflows may be crucial in enabling Lyα and Lyman-continuum to escape.

CHANG-ES XII

Context. Low-frequency radio continuum studies of star-forming edge-on galaxies can help to further understand how cosmic-ray electrons (CRe) propagate through the interstellar medium into the halo and how this is affected by energy losses and magnetic fields. Aims. Observations with the Very Large Array (VLA) from Continuum Haloes in Nearby Galaxies – an EVLA Survey (CHANG-ES) are combined with those with the LOw Frequency ARray (LOFAR) from the LOFAR Two-metre Sky Survey (LoTSS ) to identify the prevailing mode of cosmic-ray transport in the edge-on spiral galaxy NGC 3556. Methods. We mapped the radio spectral index, magnetic field strength, and orientation using VLA 1.5 and 6 GHz and LOFAR 144 MHz data, and we fit 1D cosmic-ray propagation models to these maps using SPINNAKER (Spectral Index Numerical Analysis of K(c)osmic-ray electron radio emission) and its interactive wrapper SPINTERACTIVE. Results. We find that the spectral index in the galactic midplane is, as expected for young CRe, α ≈ −0.7 and steepens towards the halo of the galaxy as a consequence of spectral ageing. The intensity scale heights are about 1.4 and 1.9 kpc for the thin disc, and 3.3 and 5.9 kpc for the thick disc at 1.5 GHz and 144 MHz, respectively. While pure diffusion cannot explain our data, advection can, particularly if we assume a linearly accelerating wind. Our best-fitting model has an initial speed of 123 km s−1 in the galactic midplane and reaches the escape velocity at heights between 5 kpc and 15 kpc above the disc, depending on the assumed dark matter halo of the galaxy. This galactic wind scenario is corroborated by the existence of vertical filaments seen both in the radio continuum and in H α in the disc-halo interface and of a large-scale reservoir of hot, X-ray emitting gas in the halo. Conclusions. Radio haloes show the existence of galactic winds, possibly driven by cosmic rays, in typical star-forming spiral galaxies.

Two-face(s): ionized and neutral gas winds in the local Universe

We present a comprehensive study of the Na Iλ5890, 5895 (Na I D) resonant lines in the Sloan Digital Sky Survey (SDSS, DR7) spectroscopic sample to look for neutral gas outflows in the local galaxies. Individual galaxy spectra are stacked in bins of stellar mass (M⋆) and star formation rate (SFR) to investigate the dependence of galactic wind occurrence and velocity as a function of the galaxy position in the SFR-M⋆plane. While in most cases the interstellar medium (ISM) absorption and emission lines are at the galaxy systemic velocity, in massive galaxies (M⋆ &gt; 5 × M⊙) at the high SFR tail (SFR &gt; 10−12M⊙yr−1) we find evidence of a significant blue-shifted Na I D absorption, which we interpret as evidence of neutral outflowing gas. The occurrence of the blue-shifted absorption in the stacked spectra does not depend on the contribution of the nuclear activity as it is observed at the same significance for purely star-forming (SF) galaxies, active galactic nuclei (AGN), and composite systems at fixed SFR. We confirm, instead, for all classes of objects a clear dependence on the galaxy disc inclination: the blueshift is the largest and the Na I D equivalent width the smallest for face-on galaxies, while the absorption feature is at the systemic velocity for edge-on systems. This indicates that the neutral outflow is mostly perpendicular or biconical with respect to the galactic disc. We also compare the kinematics of the neutral gas with the ionized gas phase as traced by the [OIII]λ5007, Hα, [NII]λ6548, and [NII]λ6584 emission lines in the same galaxy spectra. Differently for the neutral gas phase, all the emission lines show evidence of perturbed kinematics only in galaxies with a significant level of nuclear activity, and they are independent of the galactic disc inclination. This would suggest that, while neutral winds originate from the galactic disc and are powered by SF feedback, ionized outflows are instead due to AGN feedback originating from the black hole accretion disc. In both the neutral and ionized gas phases, the observed wind velocities (of the order of 100−200 kms−1) suggest that the outflowing gas remains bound to the galaxy with no definitive effect on the gas reservoir.

Ly α properties of simulated galaxies in overdense regions: effects of galactic winds at $z \gtrsim 6$

We perform Monte-Carlo radiative transfer calculations to model the Lyman alpha properties of galaxies in high-resolution, zoom-in cosmological simulations at z ~ 6.6. The simulations include both constrained and unconstrained runs, representing respectively a highly overdense region and an average field. Different galactic wind models are used in the simulations in order to investigate the effects of these winds on the apparent Ly alpha properties of galaxies. We find that, for models including galactic winds, the Ly alpha properties of massive galaxies residing in the overdense region match well recent observations of luminous Ly alpha emitters (LAEs) at z ~ 6-7, in terms of apparent Ly alpha luminosity, Ly alpha line width and Ly alpha equivalent width distributions. Without winds, the same galaxies appear less Ly alpha bright as a result of both differences in the line profile emerging from galaxies themselves, and, in the distributions of neutral gas in the circumgalactic (CGM) and intergalactic medium (IGM). We also study the relations between apparent Ly alpha luminosity and various galaxy properties: stellar mass, star formation rate (SFR) and host halo mass. At fixed halo mass, the apparent Ly alpha luminosity of galaxies appears to depend on the large-scale environment while this is no longer true for galaxies at a given stellar mass or SFR. We provide simple linear fits to these relations that can be used for quickly constructing mock LAE samples from N-body simulations. Our results suggest that the observed luminous LAEs at z ~ 6.6 are hosted by ~10^{12} h^{-1} Mo, dark matter haloes, residing in large, overdense ionized regions.

Gamma-Ray Production in the Extended Halo of the Galaxy and Possible Implications for the Origin of Galactic Cosmic Rays

Various studies have implied the existence of a gaseous halo around the Galaxy extending out to 100 kpc. Galactic cosmic rays (CRs) that propagate to the halo, either by diffusion or by convection with the possibly existing large-scale Galactic wind, can interact with the gas therein and produce gamma-rays via proton-proton collision. We calculate the cosmic ray distribution in the halo and the gamma-ray flux, and explore the dependence of the result on model parameters such as diffusion coefficient, CR luminosity, CR spectral index. We find that the current measurement of isotropic gamma-ray background at $\lesssim$TeV with Fermi Large Area Telescope already approaches a level that can provide interesting constraints on the properties of Galactic cosmic ray (e.g., with CR luminosity $L_{CR}\leq 10^{41}$erg/s). We also discuss the possibilities of the Fermi bubble and IceCube neutrinos originating from the proton-proton collision between cosmic rays and gas in the halo, as well as the implication of our results for the baryon budget of the hot circumgalactic medium of our Galaxy. Given that the isotropic gamma-ray background is likely to be dominated by unresolved extragalactic sources, future telescopes may extract more individual sources from the IGRB, and hence put even more stringent restriction on the relevant quantities (such as Galactic cosmic ray luminosity and baryon budget in the halo) in the presence of a turbulent halo that we consider.

Osaka feedback model: isolated disc galaxy simulations

We study various implementations of supernova feedback model and present the results of our `Osaka feedback model' using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {\small GADGET-3}. Our model is a modified version of Stinson et al.'s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical properties of SN remnants from the results of Chevalier and McKee \& Ostriker, such as the effective radius of SN bubble and the remnant life-time, in the form of Sedov-Taylor (ST)-like solutions with the effect of radiative cooling. Our model utilizes the local, physical parameters such as density and temperature of the ISM rather than galactic or halo properties to determine the galactic wind velocity or mass-loading factor. The Osaka model succeeds in self-regulating star formation, and naturally produces galactic outflow with variable velocities depending on the local environment and available SN energy as a function of time.An important addition to our previous work by Aoyama et al. is the implementation of the {\small CELib} chemistry library which allows us to deal with the time-dependent input of energy and metal yields for type Ia \& II supernovae (SNe) and asymptotic giant branch (AGB) stars. As initial tests of our model, we apply it to isolated galaxy simulations, and examine various galactic properties and compare with observational data including metal abundances.