In order to study gas evolution in the central region of a barred galaxy, we have performed numerical simulations of gas in the potential of the barred galaxy. We have found that the bar potential produces a gas ring within the central 1 kpc region. In the gas ring, active star and star cluster formations take place. Since the gas ring is dense enough to become self-gravitationally unstable, gas clouds form in the ring. These gas clouds interact gravitationally and collide with the other clouds. Such interaction and collision reduces their angular momentums effectively, and finally gas clouds fall into the galactic center. These processes triggers episodic gas fueling to the galactic center.

Radio-mode feedback from relativistic jets is one of the prominent heating mechanisms in clusters of galaxies. We present a long-term evolution of high-resolution MHD simulation of jets interacting with an environment modeled to represent the Perseus cluster. We investigate the thermodynamics of the ICM due to the gas motion triggered by the action of the jets and show that low-entropy gas is lifted efficiently in the wake of the inflating radio lobe. We look into the uplift mechanism and estimate the energy budget and the rate of thermal conduction. The redistribution of entropy suggests that heat conduction can play a more significant role in the thermal evolution of the cluster core in the presence of jets, which act effectively as a heat pump, thus heating the ICM more efficiently than jets would by themselves in an isentropic cluster.

Binary stars can interact via mass transfer when one member (the primary) ascends onto a giant branch. The amount of gas ejected by the binary and the amount of gas accreted by the secondary over the lifetime of the primary influence the subsequent binary phenomenology. Some of the gas ejected by the binary will remain gravitationally bound and its distribution will be closely related to the formation of planetary nebulae. We investigate the nature of mass transfer in binary systems containing an AGB star by adding radiative transfer to the AstroBEAR AMR Hydro/MHD code.

Intrinsic alignments are believed to be a major source of systematics for future generation of weak gravitational lensing surveys like Euclid or LSST. Direct measurements of the alignment of the projected light distribution of galaxies in wide field imaging data seem to agree on a contamination at a level of a few per cent of the shear correlation functions, although the amplitude of the effect depends on the population of galaxies considered. Given this dependency, it is difficult to use dark matter-only simulations as the sole resource to predict and control intrinsic alignments. We report here estimates on the level of intrinsic alignment in the cosmological hydrodynamical simulation Horizon-AGN that could be a major source of systematic errors in weak gravitational lensing measurements. In particular, assuming that the spin of galaxies is a good proxy for their ellipticity, we show how those spins are spatially correlated and how they couple to the tidal field in which they are embedded. We will also present theoretical calculations that illustrate and qualitatively explain the observed signals.

Galactic disks consist of both stars and gas. The stars gravitationally influence the gas either in disks at large or within spiral arms, leading to the formation of giant clouds and turbulence driving in the gas. In featureless disks as in flocculent galaxies, swing amplification operating in a combined star-gas disk is efficient to form bound condensations and feed a significant level of random gas motions. This occurs when the gaseous Toomre parameter is less than 1.4 for the stellar parameters similar to the solar neighbourhood conditions. In disks with spiral features, on the other hand, spiral-arm spurs and associated giant clouds develop as a consequence of magneto-Jeans instability in which magnetic tension counterbalances the stabilizing Coriolis force. Spiral shocks are inherently unstable when the vertical dimension is taken into account, exhibiting flapping motions of the shock front. This naturally converts the kinetic energy in galaxy rotation into random kinetic energy of the gas. The resulting turbulent motions are supersonic and persist despite strong shock dissipation. Thermal instability occurring in gas flows across spiral arms prompts phases transitions that produce a significant fraction of thermally-unstable, intermediate-temperature gas in the postshock expansion zones.

We discuss recent advances in cloud formation via gravitational instability under the action of self-gravity, magnetic fields, rotational shear, active stars, and/or stellar spiral arms. When shear is strong and the spiral arms are weak, applicable to flocculent galaxies at large, swing amplification exhibits nonlinear threshold behavior such that disks with a Toomre parameter Q < Qc experience gravitational runaway. For most realistic conditions, local models yield Qc ~ 1.4, similar to the observed star formation thresholds. When shear is weak, on the other hand, as in galactic central parts or inside spiral arms, magneto-Jeans instability is very powerful to form spiral-arm substructures including gaseous spurs and giant clouds. The wiggle and Parker instabilities proposed for cloud formation appear to be suppressed by strong non-steady motions inherent in vertically-extended spiral shocks, suggesting that gravitational instability is a primary candidate for cloud formation.

It has been recently shown by several authors that fragmentation of pre-stellar gas (i.e. at densities from 104 to 1010 particles cm−3 and temperatures of order 10-30 K) depends on the gas thermodynamics much more than it was anticipated in earlier studies, in which only an isothermal behaviour has been assumed for the gas. Here we review the results of a number of numerical hydrodynamic simulations (e.g. Li et al. 2003, Jappsen et al. 2005, Bonnell et al. 2006) in which departure from isothermality has been attempted by employing a polytropic equation of state (eos) with exponent different from unity. In particular, in these studies it has been shown that the dominant fragmentation scale of pre-stellar gas, and hence the peak of the initial mass function (IMF), depends on a polytropic exponent that changes value, from below to above unity, at a critical density (Larson 2005). Furthermore, this piecewise polytropic eos depends on the gas metallicity and fundamental constants. Therefore, the peak of the IMF depends, in turn, also on the gas metallicity and fundamental constants rather than on initial conditions, as it has been previously suggested (e.g. Larson 1995). Hence, we are for the first time in a position to infer theoretically the notion of a universal IMF (at least for its low-mass end).

We also present two test cases in which a non-isothermal eos has been used in the context of smoothed particle hydrodynamic (SPH) numerical simulations. In the first case star formation is triggered by means of low-mass clump collisions. These calculations have shown that clump collisions can be a relatively efficient mechanism for the formation of solar-mass protostars and their lower-mass companions (efficiency greater or of order 20-25%; Kitsionas & Whitworth 2006). We have also found that in such collisions protostars form mainly by fragmentation of dense filaments along which it is likely that pairs of protostars capture each other in close binaries surrounded by circumbinary discs. In the second case, the use of a polytropic eos with a varying exponent appropriate for the metallicity of starburst regions (Spaans & Silk 2000, 2005) is shown to be sufficient to obtain a top heavy IMF similar to that observed e.g. in the Galactic centre (Klessen, Spaans & Jappsen 2006). These are preliminary results in the direction of revisiting earlier isothermal calculations that were resolving all densities up to the opacity limit for fragmentation (e.g. Bate et al. 2002ab, 2003), this time also taking into account the thermal properties of the gas in the density range between 104 and 1010 particles cm−3. The next step would be to include self-consistent radiation transport in the calculations, the first attempts for which are already in the making (e.g. Whitehouse & Bate 2004).