Angular momentum transport and evolution of lopsided galaxies

The surface brightness distribution in the majority of stellar galactic discs falls off exponentially. Often what lies beyond such a stellar disc is the neutral hydrogen gas whose distribution also follows a nearly exponential profile at least for a number of nearby disc galaxies. Both the stars and gas are commonly known to host lopsided asymmetry especially in the outer parts of a galaxy. The role of such asymmetry in the dynamical evolution of a galaxy has not been explored so far. Following Lindblad's original idea of kinematic density waves, we show that the outer part of an exponential disc is ideally suitable for hosting lopsided asymmetry. Further, we compute the transport of angular momentum in the combined stars and gas disc embedded in a dark matter halo. We show that in a pure star and gas disc, there is a transition point where the free precession frequency of a lopsided mode, Omega - kappa, changes from retrograde to prograde and this in turn reverses the direction of angular momentum flow in the disc leading to an unphysical behaviour. We show that this problem is overcome in the presence of a dark matter halo, which sets the angular momentum flow outwards as required for disc evolution, provided the lopsidedness is leading in nature. This, plus the well-known angular momentum transport in the inner parts due to spiral arms, can facilitate an inflow of gas from outside perhaps through the cosmic filaments.

Lopsided spiral galaxies

The light distribution in the disks of many galaxies is ‘lopsided’ with a spatial extent much larger along one half of a galaxy than the other, as seen in M101. Recent observations show that the stellar disk in a typical spiral galaxy is significantly lopsided, indicating asymmetry in the disk mass distribution. The mean amplitude of lopsidedness is 0.1, measured as the Fourier amplitude of the m=1 component normalized to the average value. Thus, lopsidedness is common, and hence it is important to understand its origin and dynamics. This is a new and exciting area in galactic structure and dynamics, in contrast to the topic of bars and two-armed spirals (m=2) which has been extensively studied in the literature. Lopsidedness is ubiquitous and occurs in a variety of settings and tracers. It is seen in both stars and gas, in the outer disk and the central region, in the field and the group galaxies. The lopsided amplitude is higher by a factor of two for galaxies in a group. The lopsidedness has a strong impact on the dynamics of the galaxy, its evolution, the star formation in it, and on the growth of the central black hole and on the nuclear fuelling. We present here an overview of the observations that measure the lopsided distribution, as well as the theoretical progress made so far to understand its origin and properties. The physical mechanisms studied for its origin include tidal encounters, gas accretion and a global gravitational instability. The related open, challenging problems in this emerging area are discussed.

GMRT observations of the group Holmberg 124: Evolution by tidal forces and ram pressure?

We report new radio continuum and 21 cm HI observations using the Giant Metrewave Radio Telescope (GMRT) of the group Holmberg 124 ( Ho 124) comprising four late-type galaxies, namely NGC 2820, Mrk 108, NGC 2814 and NGC 2805. The three galaxies, NGC 2820, Mrk 108 and NGC 2814 which are closely located in the sky plane have clearly undergone tidal interactions as seen from the various morphological tidal signatures and debris. Moreover we note various features in the group members which we believe might be due to ram pressure. In this paper, we describe four interesting results emerging from our observations: a) detection of the tidal radio continuum bridge at 330 MHz connecting the galaxies NGC 2820+ Mrk 108 with NGC 2814. The radio bridge was discovered at 1465 MHz by van der Hulst & Hummel ( 1985, A& A, 150, 17). We find that the bridge has a fairly steep spectrum with a spectral index alpha(S proportional to nu(alpha)) of - 1.8(-0.2)(+0.3) which is much steeper than the - 0.8 quoted by van der Hulst & Hummel ( 1985); b) detection of other tidal features like the tilted HI and radio continuum disk of NGC 2814, a HI streamer and a radio continuum tail arising from the south of NGC 2814. We also report the detection of a possible tidal dwarf galaxy in HI; c) sharp truncation in the HI distribution in the south of NGC 2820 and in the HI and radio continuum distribution in the north of NGC 2814. The optical disks in both the cases look undisturbed. As pointed out by Davis et al. ( 1997, AJ, 114, 613), ram pressure affects different components of the interstellar medium to varying degrees. Simple estimates of pressure in different components of the interstellar medium ( radio continuum, Ha and HI) in NGC 2820 indicate that ram pressure will significantly influence HI; d) detection of a large one-sided HI loop to the north of NGC 2820. No radio continuum emission or Ha emission is associated with the HI loop. We discuss various scenarios for the origin of this loop including a central starburst, ram pressure stripping and tidal interaction. We do not support the central starburst scenario since the loop is not detected in ionized gas. Using the upper limit on X-ray luminosity of Ho 124 (Mulchaey et al. 2003, ApJS, 145, 39), we estimate an upper limit on the intragroup medium (IGrM) density of 8.8 x 10(-4) cm(-3). For half this electron density, we estimate the ram pressure force of the IGrM to be comparable to the gravitational pull of the disk of NGC 2820. Since tidal interaction has obviously influenced the group, we suggest that the loop could have formed by ram pressure stripping if tidal effects had reduced the surface density of HI in NGC 2820. From the complex observational picture of Ho 124 and the numerical estimates, we suggest that the evolution of the Ho 124 group may be governed by both tidal forces due to the interaction and the ram pressure due to motion of the member galaxies in the IGrM and that the IGrM densities should not be too low (i.e. >= 4 x 10(-4)). However this needs to be verified by further observations.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

Lopsidedness in WHISP galaxies I. Rotation curves and kinematic lopsidedness

The frequently observed lopsidedness of the distribution of stars and gas in disc galaxies is still considered as a major problem in galaxy dynamics. It is even discussed as an imprint of the formation history of discs and the evolution of baryons in dark matter haloes. Here, we analyse a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. The Hi data allow us to follow the morphology and the kinematics out to very large radii. In the present paper, we present the rotation curves and study the kinematic asymmetry. We extract the rotation curves of the receding and approaching sides separately and show that the kinematic behaviour of disc galaxies can be classified into five different types: symmetric velocity fields where the rotation curves of the receding and approaching sides are almost identical; global distortions where the rotation velocities of the receding and approaching sides have an offset that is constant with radius; local distortions leading to large deviations in the inner and negligible deviations in the outer parts (and vice versa); and distortions that divide the galaxies into two kinematic systems that are visible in terms of the different behaviour of the rotation curves of the receding and approaching sides, which leads to a crossing and a change in side. The kinematic lopsidedness is measured from the maximum rotation velocities, averaged over the plateau of the rotation curves. This gives a good estimate of the global lopsidedness in the outer parts of the sample galaxies. We find that the mean value of the perturbation parameter denoting the lopsided potential as obtained from the kinematic data is 0.056. Altogether, 36% of the sample galaxies are globally lopsided, which can be interpreted as the disc responding to a halo that was distorted by a tidal encounter. In Paper II, we study the morphological lopsidedness of the same sample of galaxies.

Lopsidedness in WHISP galaxies II. Morphological lopsidedness

The distribution of stars and gas in many galaxies is asymmetric. This so-called lopsidedness is expected to significantly affect the dynamics and evolution of the disc, including the star formation activity. Here, we measure the degree of lopsidedness for the gas distribution in a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. This complements our earlier work (Paper I) where the kinematic lopsidedness was derived for the same galaxies. The morphological lopsidedness is measured by performing a harmonic decomposition of the surface density maps. The amplitude of lopsidedness A(1), the fractional value of the first Fourier component, is typically quite high (about 0.1) within the optical disc and has a constant phase. Thus, lopsidedness is a common feature in galaxies and indicates a global mode. We measure A(1) out to typically one to four optical radii, sometimes even further. This is, on average, four times larger than the distance to which lopsidedness was measured in the past using near-IR as a tracer of the old stellar component, and therefore provides a new, more stringent constraint on the mechanism for the origin of lopsidedness. Interestingly, the value of A(1) saturates beyond the optical radius. Furthermore, the plot of A(1) versus radius shows fluctuations that we argue are due to local spiral features. We also try to explain the physical origin of this observed disc lopsidedness. No clear trend is found when the degree of lopsidedness is compared to a measure of the isolation or interaction probability of the sample galaxies. However, this does not rule out a tidal origin if the lopsidedness is long-lived. In addition, we find that the early-type galaxies tend to be more morphologically lopsided than the late-type galaxies. Both results together indicate that lopsidedness has a tidal origin.

Lopsidedness and sloshing in centres of advanced mergers of galaxies

We measure the non-axisymmetry in the luminosity distribution in the central few kpc of a sample of advanced mergers of galaxies, by analyzing their 2MASS images. All mergers show a high central asymmetry: the centres of isophotes show a striking sloshing pattern with a spatial variation of upto 30% within the central 1 kpc; and the Fourier amplitude for lopsidedness (m = 1) shows high values upto 0.2 within the central 5 kpc. The central asymmetry is estimated to be long-lived, lasting for ~ a few Gyr or ~ 100 local dynamical timescales. This will significantly affect the dynamical evolution of this region, by helping fuel the central active galactic nucleus, and also by causing the secular growth of the bulge driven by lopsidedness.

Superbubble breakout and galactic winds from disc galaxies

We study the conditions for disc galaxies to produce superbubbles that can break out of the disc and produce a galactic wind. We argue that the threshold surface density of supernovae rate for seeding a wind depends on the ability of superbubble energetics to compensate for radiative cooling. We first adapt Kompaneets formalism for expanding bubbles in a stratified medium to the case of continuous energy injection and include the effects of radiative cooling in the shell. With the help of hydrodynamic simulations, we then study the evolution of superbubbles evolving in stratified discs with typical disc parameters. We identify two crucial energy injection rates that differ in their effects, the corresponding breakout ranging from being gentle to a vigorous one. (a) Superbubbles that break out of the disc with a Mach number of the order of 2-3 correspond to an energy injection rate of the order of 10(-4) erg cm(-2) s(-1), which is relevant for disc galaxies with synchrotron emitting gas in the extra-planar regions. (b) A larger energy injection threshold, of the order of 10(-3) erg cm(-2) s(-1), or equivalently, a star formation surface density of similar to 0.1 M-circle dot yr(-1) kpc(-2), corresponds to superbubbles with a Mach number similar to 5-10. While the milder superbubbles can be produced by large OB associations, the latter kind requires super-starclusters. These derived conditions compare well with observations of disc galaxies with winds and the existence of multiphase halo gas. Furthermore, we find that contrary to the general belief that superbubbles fragment through Rayleigh-Taylor (RT) instability when they reach a vertical height of the order of the scaleheight, the superbubbles are first affected by thermal instability for typical disc parameters and that RT instability takes over when the shells reach a distance of approximately twice the scaleheight.

Narrow escape: how ionizing photons escape from disc galaxies

In this paper we calculate the escape fraction (f(esc)) of ionizing photons from starburst galaxies. Using 2D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the escape fraction of ionizing photons from the centre of the disc along different angles through the superbubble and the gas disc. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of similar to 40 degrees, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scaleheights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed similar to1 - cos (1 rad)] = 0.5 from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time-and angle-averaged escape fraction on the mid-plane disc gas density (in the range n(0) = 0.15-50 cm(-3)) and the disc scaleheight (between z(0) = 10 and 600 pc). We find that the escape fraction is related to the disc parameters (the mid-plane disc density and scaleheight) roughly so that f(esc)(alpha)n(0)(2)z(0)(3) (with alpha approximate to 2.2) is a constant. For discs with a given warm neutral medium temperature, massive discs have lower escape fraction than low-mass galaxies. For Milky Way ISM parameters, we find f(esc) similar to 5 per cent, and it increases to approximate to 10 per cent for a galaxy 10 times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the Universe.

Evolution of multivariant microstructures with anisotropic misfit: A phase field study

We study, in two dimensions, the effect of misfit anisotropy on microstructural evolution during precipitation of an ordered beta phase from a disordered alpha matrix; these phases have, respectively, 2- and 6-fold rotation symmetries. Thus, precipitation produces three orientational variants of beta phase particles, and they have an anisotropic (and crystallographically equivalent) misfit strain with the matrix. The anisotropy in misfit is characterized using a parameter t = epsilon(yy)/epsilon(xx), where epsilon(xx) and epsilon(yy) are the principal components of the misfit strain tensor. Our phase field, simulations show that the morphology of beta phase particles is significantly influenced by 1, the level of misfit anisotropy. Particles are circular in systems with dilatational misfit (t = 1), elongated along the direction of lower principal misfit when 0 < t < 1 and elongated along the invariant direction when - 1 <= t <= 0. In the special case of a pure shear misfit strain (t = - 1), the microstructure exhibits star, wedge and checkerboard patterns; these microstructural features are in agreement with those in Ti-Al-Nb alloys.

Nonlinear Enzyme Kinetics Can Lead to High Metabolic Flux Control Coefficients: Implications for the Evolution of Dominance

In a classic study, Kacser & Burns (1981, Genetics 97, 639-666) demonstrated that given certain plausible assumptions, the flux in a metabolic pathway was more or less indifferent to the activity of any of the enzymes in the pathway taken singly. It was inferred from this that the observed dominance of most wild-type alleles with respect to loss-of-function mutations did not require an adaptive, meaning selectionist, explanation. Cornish-Bowden (1987, J. theor. Biol. 125, 333-338) showed that the Kacser-Burns inference was not valid when substrate concentrations were large relative to the relevant Michaelis constants. We find that in a randomly constructed functional pathway, even when substrate levels are small, one can expect high values of control coefficients for metabolic flux in the presence of significant nonlinearities as exemplified by enzymes with Hill coefficients ranging from two to six, or by the existence of oscillatory loops. Under these conditions the flux can be quite sensitive to changes in enzyme activity as might be caused by inactivating one of the two alleles in a diploid. Therefore, the phenomenon of dominance cannot be a trivial ''default'' consequence of physiology but must be intimately linked to the manner in which metabolic networks have been moulded by natural selection.

Evolution of electronic structure with dimensionality in divalent nickelates

We investigate the evolution of electronic structure with dimensionality (d) of Ni-O-Ni connectivity in divalent nickelates, NiO (3-d), La2NiO4, Pr2NiO4 (2-d), Y2BaNiO5 (1-d) and Lu2BaNi5 (0-d), by analyzing the valence band and the Ni 2p core-level photoemission spectra in conjunction with detailed many-body calculations including full multiplet interactions. Experimental results exhibit a reduction in the intensity of correlation-induced satellite features with decreasing dimensionality. The calculations based on the cluster model, but evaluating both Ni 3d and O 2p related photoemission processes on the same footing, provide a consistent description of both valence-band and core-level spectra in terms of various interaction strengths. While the correlation-induced satellite features in NiO is dominated by poorly screened d(8) states as described in the existing literature, we find that the satellite features in the nickelates with lower dimensional Ni-O-Ni connectivity are in fact dominated by the over-screened d(10)L(2) states. It is found that the changing electronic structure with the dimensionality is primarily driven by two factors: (i) a suppression of the nonlocal contribution to screening; and (ii) a systematic decrease of the charge-transfer energy Delta driven by changes in the Madelung potential. [S0163-1829(99)09619-8].

Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

Microstructural evolution and mechanical characteristics in nanocrystalline nickel with a bimodal grain-size distribution

An attractive microstructural possibility for enhancing the ductility of high-strength nanocrystals is to develop a bimodal grain-size distribution, in which the fine grains provide strength, and the coarser grains enable strain hardening. Annealing of nanocrystalline Ni over a range of temperatures and times led to microstructures with varying volume fractions of coarse grains and a change in texture. Tensile tests revealed a drastic reduction in ductility with increasing volume fraction of coarse grains. The reduction in ductility may be related to the segregation of sulphur to grain boundaries.

Study of texture evolution in metastable beta-Ti alloy as a function of strain path and its effect on alpha transformation texture

Texture evolution in a low cost beta titanium alloy was studied for different modes of rolling and heat treatments. The alloy was cold rolled by unidirectional and multi-step cross rolling. The cold rolled material was either aged directly or recrystallized and then aged. The evolution of texture in alpha and beta phases were studied. The rolling texture of beta phase that is characterized by the gamma fiber is stronger for MSCR than UDR; while the trend is reversed on recrystallization. The mode of rolling affects alpha transformation texture on aging with smaller alpha lath size and stronger alpha texture in UDR than in MSCR. The defect structure in beta phase influences the evolution of a texture on aging. A stronger defect structure in beta phase leads to variant selection with the rolled samples showing fewer variants than the recrystallized samples.