We study the effects of Kelvin Helmholtz Instability (KHI) on the cold streams that feed massive halos at high redshift, generalizing our earlier results to include the effects of radiative cooling and heating from a UV background, using analytic models and high resolution idealized simulations. We currently do not consider self-shielding, thermal conduction or gravity. A key parameter in determining the fate of the streams is the ratio of the cooling time in the turbulent mixing layer which forms between the stream and the background following the onset of the instability, $t_{\rm cool,mix}$, to the time in which the mixing layer expands to the width of the stream in the non-radiative case, $t_{\rm shear}$. This can be converted into a critical stream radius, $R_{\rm s,crit}$, such that $R_{\rm s}/R_{\rm s,crit}=t_{\rm shear}/t_{\rm cool,mix}$. If $R_{\rm s}<R_{\rm s,crit}$, the non-linear evolution proceeds similarly to the non-radiative case studied by Mandelker et al. 2019a. If $R_{\rm s}>R_{\rm s,crit}$, which we find to almost always be the case for astrophysical cold streams, the stream is not disrupted by KHI. Rather, background mass cools and condenses onto the stream, and can increase the mass of cold gas by a factor of $\sim 3$ within 10 stream sound crossing times. The mass entrainment induces thermal energy losses from the background and kinetic energy losses from the stream, which we model analytically. Roughly half of the dissipated energy is radiated away from gas with $T<5\times 10^4{\rm K}$, suggesting much of it will be emitted in Ly$\alpha$.

We present an analytic toy model for the radiation produced by the interaction between the cold streams thought to feed massive halos at high redshift and their hot CGM. We begin by deriving cosmologically motivated parameters for the streams as they enter the halo virial radius, $R_{\rm v}$, as a function of halo mass and redshift. For $10^{12}M_{\odot}$ halos at $z=2$, we find the Hydrogen number density in streams to be $n_{\rm H,s}\sim (0.1-5)\times 10^{-2}{\rm cm}^{-3}$, a factor of $\delta \sim (30-300)$ times denser than the hot CGM density, while the stream radii are in the range $R_{\rm s}\sim (0.03-0.50)R_{\rm v}$. As the streams accelerate towards the halo centre, they become denser and narrower. The stream-CGM interaction induces Kelvin-Helmholtz Instability (KHI), which leads to entrainment of CGM mass by the stream and therefore to stream deceleration by momentum conservation. Assuming that the entrainment rates derived by Mandelker et al. 2019 in the absence of gravity can be applied locally at each halocentric radius, we derive equations of motion for the stream in the halo. Using these, we derive the net acceleration, mass growth, and energy dissipation induced by the stream-CGM interaction, as a function of halo mass and redshift, for different CGM density profiles. For the range of model parameters considered, we find that the interaction can induce dissipation luminosities $L_{\rm diss}>10^{42}~{\rm erg~s^{-1}}$ within $\le 0.6 R_{\rm v}$ of halos with $M_{\rm v}>10^{12}M_{\odot}$ at $z=2$, with the emission scaling with halo mass and redshift approximately as $\propto M_{\rm v}\,(1+z)^2$. The magnitude and spatial extent of the emission produced in massive halos at high redshift is consistent with observed Ly$\alpha$ blobs, though better treatment of the UV background and self-shielding is needed to solidify this conclusion.

Self-gravitating gaseous filaments exist on many astrophysical scales, from sub-pc filaments in the interstellar medium to Mpc scale streams feeding galaxies from the cosmic web. These filaments are often subject to Kelvin-Helmotz Instability (KHI) due to shearing against a confining background medium. We study the nonlinear evolution of KHI in pressure-confined self-gravitating gas streams initially in hydrostatic equilibrium, using analytic models and hydrodynamic simulations, not including radiative cooling. We derive a critical line-mass as a function of the stream Mach number and density contrast with respect to the background, $\mu_{cr}(M_b,\delta_c)\le 1$, where $\mu=1$ is normalized to the maximal line mass for which initial hydrostatic equilibrium is possible. For $\mu<\mu_{cr}$, KHI dominates the stream evolution. A turbulent shear layer expands into the background and leads to stream deceleration at a similar rate to the non-gravitating case. However, with gravity, penetration of the shear layer into the stream is halted at roughly half the initial stream radius by stabilizing buoyancy forces, significantly delaying total stream disruption. Streams with $\mu_{cr}<\mu\le 1$ fragment and form round, long-lived clumps by gravitational instability (GI), with typical separations roughly 8 times the stream radius, similar to the case without KHI. When KHI is still somewhat effective, these clumps are below the spherical Jeans mass and are partially confined by external pressure, but they approach the Jeans mass as $\mu\rightarrow 1$ and GI dominates. We discuss potential applications of our results to filaments in the ISM and dense streams feeding galaxies at high redshift.

We study the effects of Kelvin-Helmholtz instability (KHI) on the cold streams that feed high-redshift galaxies through their hot haloes, generalizing our earlier analyses of a 2D slab to a 3D cylinder, but still limiting our analysis to the adiabatic case with no gravity. We combine analytic modeling and numerical simulations in the linear and non-linear regimes. For subsonic or transonic streams with respect to the halo sound speed, the instability in 3D is qualitatively similar to 2D, but progresses at a faster pace. For supersonic streams, the instability grows much faster in 3D and can be qualitatively different due to azimuthal modes, which introduce a strong dependence on the initial width of the stream-background interface. Using analytic toy models and approximations supported by high-resolution simulations, we apply our idealized hydrodynamical analysis to the astrophysical scenario. The upper limit for the radius of a stream that disintegrates prior to reaching the central galaxy is ~70% larger than the 2D estimate; it is in the range (0.5-5)% of the halo virial radius, decreasing with increasing stream density and velocity. Stream disruption generates a turbulent mixing zone around the stream with velocities at the level of ~20% of the initial stream velocity. KHI can cause significant stream deceleration and energy dissipation in 3d, contrary to 2D estimates. For typical streams, up to (10-50)% of the gravitational energy gained by inflow down the dark-matter halo potential can be dissipated, capable of powering Lyman-alpha blobs if most of it is dissipated into radiation.

As part of our long-term campaign to understand how cold streams feed massive galaxies at high redshift, we study the Kelvin-Helmholtz instability (KHI) of a supersonic, cold, dense gas stream as it penetrates through a hot, dilute circumgalactic medium (CGM). A linear analysis (Paper I) showed that, for realistic conditions, KHI may produce nonlinear perturbations to the stream during infall. Therefore, we proceed here to study the nonlinear stage of KHI, still limited to a two-dimensional slab with no radiative cooling or gravity. Using analytic models and numerical simulations, we examine stream breakup, deceleration and heating via surface modes and body modes. The relevant parameters are the density contrast between stream and CGM ($\delta$), the Mach number of the stream velocity with respect to the CGM ($M_{\rm b}$) and the stream radius relative to the halo virial radius ($R_{\rm s}/R_{\rm v}$). We find that sufficiently thin streams disintegrate prior to reaching the central galaxy. The condition for breakup ranges from $R_{\rm s} < 0.03 R_{\rm v}$ for $(M_{\rm b} \sim 0.75, \delta \sim 10)$ to $R_{\rm s} < 0.003 R_{\rm v}$ for $(M_{\rm b} \sim 2.25, \delta \sim 100)$. However, due to the large stream inertia, KHI has only a small effect on the stream inflow rate and a small contribution to heating and subsequent Lyman-$\alpha$ cooling emission.

Cold Fronts and shocks are hallmarks of the complex intra-cluster medium (ICM) in galaxy clusters. They are thought to occur due to gas motions within the ICM and are often attributed to galaxy mergers within the cluster. Using hydro-cosmological simulations of clusters of galaxies, we show that collisions of inflowing gas streams, seen to penetrate to the very centre of about half the clusters, offer an additional mechanism for the formation of shocks and cold fronts in cluster cores. Unlike episodic merger events, a gas stream inflow persists over a period of several Gyrs and it could generate a particular pattern of multiple cold fronts and shocks.

Interstellar gas clouds are often both highly magnetized and supersonically turbulent, with velocity dispersions set by a competition between driving and dissipation. This balance has been studied extensively in the context of gases with constant mean density. However, many astrophysical systems are contracting under the influence of external pressure or gravity, and the balance between driving and dissipation in a contracting, magnetized medium has yet to be studied. In this paper we present three-dimensional (3D) magnetohydrodynamic (MHD) simulations of compression in a turbulent, magnetized medium that resembles the physical conditions inside molecular clouds. We find that in some circumstances the combination of compression and magnetic fields leads to a rate of turbulent dissipation far less than that observed in non-magnetized gas, or in non-compressing magnetized gas. As a result, a compressing, magnetized gas reaches an equilibrium velocity dispersion much greater than would be expected for either the hydrodynamic or the non-compressing case. We use the simulation results to construct an analytic model that gives an effective equation of state for a coarse-grained parcel of the gas, in the form of an ideal equation of state with a polytropic index that depends on the dissipation and energy transfer rates between the magnetic and turbulent components. We argue that the reduced dissipation rate and larger equilibrium velocity dispersion produced by compressing, magnetized turbulence has important implications for the driving and maintenance of turbulence in molecular clouds, and for the rates of chemical and radiative processes that are sensitive to shocks and dissipation.

Observations of soft gamma rays emanating from the Milky Way from SPI/\textit{INTEGRAL} reveal the annihilation of $\sim2\times10^{43}$ positrons every second in the Galactic bulge. The origin of these positrons, which annihilate to produce a prominent emission line centered at 511 keV, has remained mysterious since their discovery almost 50 years ago. A plausible origin for the positrons is in association with the intense star formation ongoing in the Galactic center. Moreover, there is strong evidence for a nuclear outflow in the Milky Way. We find that advective transport and subsequent annihilation of positrons in such an outflow cannot simultaneously replicate the observed morphology of positron annihilation in the Galactic bulge and satisfy the requirement that $90$ per cent of positrons annihilate once the outflow has cooled to $10^4\,\mathrm{K}$.

We report an integral-field spectroscopic (IFS) observation of a gravitationally lensed spiral galaxy A1689B11 at redshift $z=2.54$. It is the most ancient spiral galaxy discovered to date and the second kinematically confirmed spiral at $z\gtrsim2$. Thanks to gravitational lensing, this is also by far the deepest IFS observation with the highest spatial resolution ($\sim$ 400 pc) on a spiral galaxy at a cosmic time when the Hubble sequence is about to emerge. After correcting for a lensing magnification of 7.2 $\pm$ 0.8, this primitive spiral disk has an intrinsic star formation rate of 22 $\pm$ 2 $M_{\odot}$ yr$^{-1}$, a stellar mass of 10$^{9.8 \pm 0.3}$$M_{\odot}$ and a half-light radius of $r_{1/2}=2.6 \pm 0.7$ kpc, typical of a main-sequence star-forming (SF) galaxy at $z\sim2$. However, the H\alpha\ kinematics show a surprisingly tranquil velocity field with an ordered rotation ($V_{\rm c}$ = 200 $\pm$ 12 km/s) and uniformly small velocity dispersions ($V_{\rm \sigma, mean}$ = 23 $\pm$ 4 km/s and $V_{\rm \sigma, outer-disk}$ = 15 $\pm$ 2 km/s). The low gas velocity dispersion is similar to local spiral galaxies and is consistent with the classic density wave theory where spiral arms form in dynamically cold and thin disks. We speculate that A1689B11 belongs to a population of rare spiral galaxies at $z\gtrsim2$ that mark the formation epoch of thin disks. Future observations with JWST will greatly increase the sample of these rare galaxies and unveil the earliest onset of spiral arms.

Virial shocks at edges of cosmic-web structures are a clear prediction of standard structure formation theories. We derive a criterion for the stability of the post-shock gas and of the virial shock itself in spherical, filamentary and planar infall geometries. When gas cooling is important, we find that shocks become unstable, and gas flows uninterrupted towards the center of the respective halo, filament or sheet. For filaments, we impose this criterion on self-similar infall solutions. We find that instability is expected for filament masses between $10^{11}-10^{13}M_\odot Mpc^{-1}.$ Using a simplified toy model, we then show that these filaments will likely feed halos with $10^{10}M_{\odot}\lesssim M_{halo}\lesssim 10^{13}M_{\odot}$ at redshift $z=3$, as well as $10^{12}M_{\odot}\lesssim M_{halo}\lesssim 10^{15}M_{\odot}$ at $z=0$. The instability will affect the survivability of the filaments as they penetrate gaseous halos in a non-trivial way. Additionally, smaller halos accreting onto non-stable filaments will not be subject to ram-pressure inside the filaments. The instreaming gas will continue towards the center, and stop either once its angular momentum balances the gravitational attraction, or when its density becomes so high that it becomes self-shielded to radiation.

Massive galaxies at high redshift are predicted to be fed from the cosmic web by narrow, dense, cold streams. These streams penetrate supersonically through the hot medium encompassed by a stable shock near the virial radius of the dark-matter halo. Our long-term goal is to explore the heating and dissipation rate of the streams and their fragmentation and possible breakup, in order to understand how galaxies are fed, and how this affects their star-formation rate and morphology. We present here the first step, where we analyze the linear Kelvin-Helmholtz instability (KHI) of a cold, dense slab or cylinder flowing through a hot, dilute medium in the transonic regime. The current analysis is limited to the adiabatic case with no gravity and assuming equal pressure in the stream and the medium. By analytically solving the linear dispersion relation, we find a transition from a dominance of the familiar rapidly growing surface modes in the subsonic regime to more slowly growing body modes in the supersonic regime. The system is parameterized by three parameters: the density contrast between the stream and the medium, the Mach number of stream velocity with respect to the medium, and the stream width with respect to the halo virial radius. We find that a realistic choice for these parameters places the streams near the mode transition, with the KHI exponential-growth time in the range 0.01-10 virial crossing times for a perturbation wavelength comparable to the stream width. We confirm our analytic predictions with idealized hydrodynamical simulations. Our linear-KHI estimates thus indicate that KHI may in principle be effective in the evolution of streams by the time they reach the galaxy. More definite conclusions await the extension of the analysis to the nonlinear regime and the inclusion of cooling, thermal conduction, the halo potential well, self-gravity and magnetic fields.

We utilize cosmological simulations of 16 galaxy clusters at redshifts \zeq{0} and \zeq{0.6} to study the effect of inflowing streams on the properties of the X-ray emitting intracluster medium. We find that the mass accretion occurs predominantly along streams that originate from the cosmic web and consist of heated gas. Clusters that are unrelaxed in terms of their X-ray morphology are characterized by higher mass inflow rates and deeper penetration of the streams, typically into the inner third of the virial radius. The penetrating streams generate elevated random motions, bulk flows and cold fronts. The degree of penetration of the streams may change over time such that clusters can switch from being unrelaxed to relaxed over a time-scale of several giga years.

Galaxy evolution depends strongly on the physics of the interstellar medium (ISM). Motivated by the need to incorporate the properties of the ISM in cosmological simulations we construct a simple method to include the contribution of non-thermal components in the calculation of pressure of interstellar gas. In our method we treat three non-thermal components - turbulence, magnetic fields and cosmic rays - and effectively parametrize their amplitude. We assume that the three components settle into a quasi-steady-state that is governed by the star formation rate, and calibrate their magnitude and density dependence by the observed Radio-FIR correlation, relating synchrotron radiation to star formation rates of galaxies. We implement our model in single cell numerical simulation of a parcel of gas with constant pressure boundary conditions and demonstrate its effect and potential. Then, the non-thermal pressure model is incorporated into RAMSES and hydrodynamic simulations of isolated galaxies with and without the non-thermal pressure model are presented and studied. Specifically, we demonstrate that the inclusion of realistic non-thermal pressure reduces the star formation rate by an order of magnitude and increases the gas depletion time by as much. We conclude that the non-thermal pressure can prolong the star formation epoch and achieve consistency with observations without invoking artificially strong stellar feedback.

Many galaxy clusters pose a "cooling-flow problem", where the observed X-ray emission from their cores is not accompanied by enough cold gas or star formation. A continuous energy source is required to balance the cooling rate over the whole core volume. We address the feasibility of a gravitational heating mechanism, utilizing the gravitational energy released by the gas that streams into the potential well of the cluster dark-matter halo. We focus here on a specific form of gravitational heating in which the energy is transferred to the medium thorough the drag exerted on inflowing gas clumps. Using spheri-symmetric hydro simulations with a subgrid representation of these clumps, we confirm our earlier estimates that in haloes >=10^13 solar masses the gravitational heating is more efficient than the cooling everywhere. The worry was that this could overheat the core and generate an instability that might push it away from equilibrium. However, we find that the overheating does not change the global halo properties, and that convection can stabilize the cluster by carrying energy away from the overheated core. In a typical rich cluster of 10^{14-15}solar masses, with ~5% of the accreted baryons in gas clumps of ~10^8 solar masses, we derive upper and lower limits for the temperature and entropy profiles and show that they are consistent with those observed in cool-core clusters. We predict the density and mass of cold gas and the level of turbulence driven by the clump accretion. We conclude that gravitational heating is a feasible mechanism for preventing cooling flows in clusters.

Cold fronts (CFs) - density and temperature plasma discontinuities - are ubiquitous in cool cores of galaxy clusters, where they appear as X-ray brightness edges in the intracluster medium, nearly concentric with the cluster center. We analyze the thermodynamic profiles deprojected across core CFs found in the literature. While the pressure appears continuous across these CFs, we find that all of them require significant centripetal acceleration beneath the front. This is naturally explained by a tangential, nearly sonic bulk flow just below the CF, and a tangential shear flow involving a fair fraction of the plasma beneath the front. Such shear should generate near-equipartition magnetic fields on scales ~<50 pc from the front, and could magnetize the entire core. Such fields would explain the apparent stability of cool-core CFs and the recently reported CF-radio minihalo association.

Cold fronts (CFs) are found in most galaxy clusters, as well as in some galaxies and groups of galaxies. We propose that some CFs are relics of collisions between trailing shocks. Such a collision typically results in a spherical, factor ~1.4-2.7 density/temperature discontinuity. These CFs may be found as far as the virial shock, unlike in other CF formation models. As a demonstration of this effect, we use one dimensional simulations where halo reverberations involving periodic collisions between the virial shock and outgoing secondary shocks exist. These collisions yield a distinctive, concentric geometric sequence of CFs which trace the expansion of the virial shock.

Cold fronts (CFs) are found in most galaxy clusters, as well as in some galaxies and groups of galaxies. We propose that some CFs are relics of merging between two shocks propagating in the same direction. Such shock mergers typically result in a quasi-spherical, factor ~1.4-2.7 discontinuity in density and in temperature. These CFs may be found as far out as the virial shock, unlike what is expected in other CF formation models. As a demonstration of this effect, we use one dimensional simulations of clusters and show that shock induced cold fronts form when perturbations such as explosions or mergers occur near the cluster's centre. Perturbations at a cluster's core induce periodic merging between the virial shock and outgoing secondary shocks. These collisions yield a distinctive, concentric, geometric sequence of CFs which trace the expansion of the virial shock.

There is no shortage of energy around to solve the overcooling problem of cooling flow clusters. AGNs, as well as gravitational energy are both energetic enough to balance the cooling of cores of clusters. The challenge is to couple this energy to the baryons efficiently enough, and to distribute the energy in a manner that will not contradict observational constraints of metalicity and entropy profiles. Here we propose that if a small fraction of the baryons that are accreted to the cluster halo are in the form of cold clumps, they would interact with the hot gas component via hydrodynamic drag. We show that such clumps carry enough energy, penetrate to the center, and heat the core significantly. We then study the dynamic response of the cluster to this kind of heating using a 1D hydrodynamic simulation with sub-grid clump heating, and produce reasonable entropy profile in a dynamic self-consistent way.

Disk galaxies are in hydrostatic equilibrium along their vertical axis. The pressure allowing for this configuration consists of thermal, turbulent, magnetic and cosmic ray components. For the Milky Way(MW) the thermal pressure contributes ~10% of the total pressure near the plane, with this fraction dropping towards higher altitudes. Out of the rest, magnetic fields contribute ~1/3 of the pressure to distances of ~3kpc above the disk plane. In this letter we attempt to extrapolate these local values to high redshift, rapidly accreting, rapidly star forming disk galaxies and study the effect of the extra pressure sources on the accretion of gas onto the galaxies. In particular, magnetic field tension may convert a smooth cold-flow accretion to clumpy, irregular star formation regions and rates. The infalling gas accumulates on the edge of the magnetic fields, supported by magnetic tension. When the mass of the infalling gas exceeds some threshold mass, its gravitational force cannot be balanced by magnetic tension anymore, and it falls toward the disk's plane, rapidly making stars. Simplified estimations of this threshold mass are consistent with clumpy star formation observed in SINS, UDF, GOODS and GEMS surveys. We discuss the shortcomings of pure hydrodynamic codes in simulating the accretion of cold flows into galaxies, and emphasize the need for magneto-hydrodynamic simulations

We calculate the expected imprint of the ionized gas in the Milky-Way halo on the Cosmic Microwave Background (CMB) through the kinetic Sunyaev-Zel'dovich (kSZ) effect. Unlike other Galactic foregrounds, the halo kSZ signature covers the full sky, generates anisotropies on large angular scales, is not accompanied by spectral distortions, and could therefore be confused with primordial CMB anisotropies. We construct theoretical models for various halo components, including smooth diffuse gas, filaments of cold inflowing gas and high velocity clouds. We find that the kSZ effect for all components is above the sensitivity of the Planck satellite, over a range of angular scales. However, the typical halo contribution is well below the cosmic variance noise in the primordial CMB power spectrum. High velocity clouds could dominate the halo contribution and better observational data is required to mask them out. We derive expected kSZ maps based on existing data from tracers of the halo gas distribution, such as 21cm maps of neutral hydrogen and H_alpha maps of recombining gas. The cross-correlation of these maps with the WMAP5 data does not yield any statistically significant signal.

The massive galaxies in the young universe, ten billion years ago, formed stars at surprising intensities. Although this is commonly attributed to violent mergers, the properties of many of these galaxies are incompatible with such events, showing gas-rich, clumpy, extended rotating disks not dominated by spheroids (Genzel et al. 2006, 2008). Cosmological simulations and clustering theory are used to explore how these galaxies acquired their gas. Here we report that they are stream-fed galaxies, formed from steady, narrow, cold gas streams that penetrate the shock-heated media of massive dark matter haloes (Dekel & Birnboim 2006; Keres et al. 2005). A comparison with the observed abundance of star-forming galaxies implies that most of the input gas must rapidly convert to stars. One-third of the stream mass is in gas clumps leading to mergers of mass ratio greater than 1:10, and the rest is in smoother flows. With a merger duy cycle of 0.1, three-quarters of the galaxies forming stars at a given rate are fed by smooth streams. The rarer, submillimetre galaxies that form stars even more intensely are largely merger-induced starbursts. Unlike destructive mergers, the streams are likely to keep the rotating disk configuration intact, although turbulent and broken into giant star-forming clumps that merge into a central spheroid (Noguchi 1999; Genzel et al. 2008, Elmegreen, Bournaud & Elmegreen 2008, Dekel, Sari & Ceverino 2009). This stream-driven scenario for the formation of disks and spheroids is an alternative to the merger picture.

We consider a simple gravitational-heating mechanism for the long-term quenching of cooling flows and star formation in massive dark-matter haloes hosting ellipticals and clusters. The virial shock heating in haloes >10^12 Mo triggers quenching in 10^12-13 Mo haloes (Birnboim, Dekel & Neistein 2007). We show that the long-term quenching in haloes >Mmin~7x10^12 Mo could be due to the gravitational energy of cosmological accretion delivered to the inner-halo hot gas by cold flows via ram-pressure drag and local shocks. Mmin is obtained by comparing the gravitational power of infall into the potential well with the overall radiative cooling rate. The heating wins if the gas inner density cusp is not steeper than r^-0.5 and if the masses in the cold and hot phases are comparable. The effect is stronger at higher redshifts, making the maintenance easier also at later times. Clumps >10^5 Mo penetrate to the inner halo with sufficient kinetic energy before they disintegrate, but they have to be <10^8 Mo for the drag to do enough work in a Hubble time. Pressure confined ~10^4K clumps are stable against their own gravity and remain gaseous once below the Bonnor-Ebert mass ~10^8 Mo. They are also immune to tidal disruption. Clumps in the desired mass range could emerge by thermal instability in the outer halo if the conductivity is not too high. Alternatively, such clumps may be embedded in dark-matter subhaloes if the ionizing flux is ineffective, but they separate from their subhaloes by ram pressure before entering the inner halo. Heating by dynamical friction becomes dominant for massive satellites, which can contribute up to one third of the total gravitational heating. We conclude that gravitational heating by cosmological accretion is a viable alternative to AGN feedback as a long-term quenching mechanism.

We simulate the buildup of galaxies by spherical gas accretion through dark matter haloes, subject to the development of virial shocks. We find that a uniform cosmological accretion turns into a rapidly varying disc buildup rate. The generic sequence of events (Shocked-Accretion Massive Burst & Shutdown: SAMBA) consists of four distinct phases: a) continuous cold accretion while the halo is below a threshold mass Msh~10^12Mo, b) tentative quenching of gas supply for ~2Gyr, starting once the halo is ~Msh and growing a rapidly expanding shock, c) a massive burst due to the big crunch of ~10^11Mo gas in \~0.5Gyr, when the heated gas cools and joins new infalling gas, and d) a long-term shutdown, enhanced by a temporary shock instability in late SAMBAs, those that quench at z~2, burst at z~1 and end up quenched in 10^12-13Mo haloes today. These events occur at all redshifts in galaxies of baryonic mass \~10^11Mo and involve a substantial fraction of this mass. They arise from rather smooth accretion, or minor mergers, which, unlike major mergers, may leave the disc intact while being built in a rapid pace. The early bursts match observed maximum starbursting discs at z>~2, predicted to reside in <~10^13Mo haloes. The late bursts resemble discy LIRGs at z<~1. The tentative quenching gives rise to a substantial population of ~10^11Mo galaxies with a strongly suppressed star-formation rate at z~1-3. The predicted long-term shutdown leads to red & dead galaxies in groups. A complete shutdown in more massive clusters requires an additional quenching mechanism, as may be provided by clumpy accretion. Alternatively, the SAMBA bursts may trigger the AGN activity that couples to the hot gas above Msh and helps the required quenching. The SAMBA phenomenon is yet to be investigated using cosmological simulations.

We address the origin of the robust bi-modality observed in galaxy properties about a characteristic stellar mass ~3x10^{10}Msun. Less massive galaxies tend to be ungrouped blue star-forming discs, while more massive galaxies are typically grouped red old-star spheroids. Colour-magnitude data show a gap between the red and blue sequences, extremely red luminous galaxies already at z~1, a truncation of today's blue sequence above L_*, and massive starbursts at z~2-4. We propose that these features are driven by the thermal properties of the inflowing gas and their interplay with the clustering and feedback processes, all functions of the dark-matter halo mass and associated with a similar characteristic scale. In haloes below a critical shock-heating mass M_shock~10^{12}Msun, discs are built by cold streams, not heated by a virial shock, yielding efficient early star formation. It is regulated by supernova feedback into a long sequence of bursts in blue galaxies constrained to a "fundamental line". Cold streams penetrating through hot media in M>M_shock haloes preferentially at z>2 lead to massive starbursts in L>L_* galaxies. At z<2, in M>M_shock haloes hosting groups, the gas is heated by a virial shock, and being dilute it becomes vulnerable to feedback from energetic sources such as AGNs. This shuts off gas supply and prevents further star formation, leading by passive evolution to "red-and-dead" massive spheroids starting at z~1. A minimum in feedback efficiency near M_shock explains the observed minimum in M/L and the qualitative features of the star-formation history. The cold flows provide a hint for solving the angular-momentum problem. When these processes are incorporated in simulations they recover the main bi-modality features and solve other open puzzles.

We investigate the conditions for the existence of an expanding virial shock in the gas falling within a spherical dark-matter halo. The shock relies on pressure support by the shock-heated gas behind it. When the radiative cooling is efficient compared to the infall rate the post-shock gas becomes unstable; it collapses inwards and cannot support the shock. We find for a monoatomic gas that the shock is stable when the post-shock pressure and density obey gamma effective>10/7, with gamma effective begin the time depended equivalent to the adiabatic index. We express the effective gamma in terms of r, u and rho at the shock to obtain a simple condition for shock stability. This result is confirmed by hydrodynamical simulations, using an accurate spheri-symmetric Lagrangian code. When the stability analysis is applied in cosmology, we find that a virial shock does not develop in most haloes that form before z ~ 2, and it never forms in haloes less massive than a few 10^11 solar masses. In such haloes the infalling gas is never heated to the virial temperature, and it does not need to cool radiatively before it drops into a disc. Instead, the gas collapses at T ~ 10^4K directly into the disc. This should have nontrivial effects on the star-formation rate and on the gas removal by supernova-driven winds. Instead of radiating soft x rays, this gas would emit lyman alpha thus helping explain the low flux of soft x-ray background and the lyman alpha emitters observed at high redshift.

We propose a model for how the buildup of dark halos by merging satellites produces a characteristic inner cusp, of a density profile \rho \prop r^-a with a -> a_as > 1, as seen in cosmological N-body simulations of hierarchical clustering scenarios. Dekel, Devor & Hetzroni (2003) argue that a flat core of a<1 exerts tidal compression which prevents local deposit of satellite material; the satellite sinks intact into the halo center thus causing a rapid steepening to a>1. Using merger N-body simulations, we learn that this cusp is stable under a sequence of mergers, and derive a practical tidal mass-transfer recipe in regions where the local slope of the halo profile is a>1. According to this recipe, the ratio of mean densities of halo and initial satellite within the tidal radius equals a given function psi(a), which is significantly smaller than unity (compared to being 1 according to crude resonance criteria) and is a decreasing function of a. This decrease makes the tidal mass transfer relatively more efficient at larger a, which means steepening when a is small and flattening when a is large, thus causing converges to a stable solution. Given this mass-transfer recipe, linear perturbation analysis, supported by toy simulations, shows that a sequence of cosmological mergers with homologous satellites slowly leads to a fixed-point cusp with an asymptotic slope a_as>1. The slope depends only weakly on the fluctuation power spectrum, in agreement with cosmological simulations. During a long interim period the profile has an NFW-like shape, with a cusp of 1<a<a_as. Thus, a cusp is enforced if enough compact satellite remnants make it intact into the inner halo. In order to maintain a flat core, satellites must be disrupted outside the core, possibly as a result of a modest puffing up due to baryonic feedback.