Numerical modelling of the vertical structure and dark halo parameters in disc galaxies

The non‐linear dynamics of bending instability and vertical structure of a galactic stellar disc embedded into a spherical halo are studied with N‐body numerical modelling. Development of the bending instability in stellar galactic disc is considered as the main factor that increases the disc thickness. Correlation between the disc vertical scale height and the halo‐to‐disc mass ratio is predicted from the simulations. The method of assessment of the spherical‐to‐disc mass ratio for edge‐on spiral galaxies with a small bulge is considered. Modelling of eight edge‐on galaxies: NGC 891, NGC 4738, NGC 5170, UGC 6080, UGC 7321, UGC 8286, UGC 9422 and UGC 9556 is performed. Parameters of stellar discs, dark haloes and bulges are estimated. The lower limit of the dark‐to‐luminous mass ratio in our galaxies is of the order of one within the limits of their stellar discs. The dark haloes dominate by mass in the galaxies with very thin stellar discs (NGC 5170, UGC 7321 and UGC 8286) (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic fields and the dynamics of spiral galaxies

We investigate the dynamics of magnetic fields in spiral galaxies by performing 3D magnetohydrodynamics simulations of galactic discs subject to a spiral potential using cold gas, warm gas and a two-phase mixture of both. Recent hydrodynamic simulations have demonstrated the formation of interarm spurs as well as spiral arm molecular clouds, provided the interstellar medium model includes a cold H  i phase. We find that the main effect of adding a magnetic field to these calculations is to inhibit the formation of structure in the disc. However, provided a cold phase is included, spurs and spiral arm clumps are still present if β≳ 0.1 in the cold gas. A caveat to the two-phase calculations though is that by assuming a uniform initial distribution, β≳ 10 in the warm gas, emphasizing that models with more consistent initial conditions and thermodynamics are required. Our simulations with only warm gas do not show such structure, irrespective of the magnetic field strength.
Furthermore, we find that the introduction of a cold H  i phase naturally produces the observed degree of disorder in the magnetic field, which is again absent from simulations using only warm gas. Whilst the global magnetic field follows the large-scale gas flow, the magnetic field also contains a substantial random component that is produced by the velocity dispersion induced in the cold gas during the passage through a spiral shock. Without any cold gas, the magnetic field in the warm phase remains relatively well ordered apart from becoming compressed in the spiral shocks. Our results provide a natural explanation for the observed high proportions of disordered magnetic field in spiral galaxies and we thus predict that the relative strengths of the random and ordered components of the magnetic field observed in spiral galaxies will depend on the dynamics of spiral shocks.     

Do LINER 2 galaxies harbour low-luminosity active galactic nuclei?

We use ROSAT HRI spatial data and ASCA spectral measurements for a sample of seven nearby, early-type spiral galaxies, to address the question of whether a low-luminosity active galactic nucleus (LLAGN) is present in galaxies that have a LINER 2 classification. The brightest discrete X-ray source in the ROSAT HRI observations is invariably found to be positionally coincident with the optical galactic nucleus, and in most cases its flux dominates the X-ray emission from the central region of the galaxy. All seven galaxies have X-ray spectra consistent with a two-component, soft thermal plus hard power-law, spectral form. If we exclude the two galaxies with relatively hard X-ray spectra, NGC 3628 and NGC 4594, for which there is supporting evidence for a LLAGN (or alternatively in the case of NGC 3628 a dominant ultraluminous X-ray binary), then the remaining galaxies show surprisingly similar X-ray spectral properties. Specifically the flux ratio F _X(0.5–1)/ F _X(2–5) , which measures the relative strengths of the thermal and non-thermal emission components, shows little scatter about a mean of 0.66, a value very similar to that measured in the classic starburst galaxy NGC 253. As there is no obvious reason why the luminosity of the hard power-law continuum emanating from a putative LLAGN should be very closely correlated with the thermal emission of the surrounding region, this suggests that that the broad-band (0.5–5 keV) X-ray emission from these LINER 2 galaxies may originate in a common set of processes probably associated with the starburst phenomenon. Conversely, it appears that in many LINER 2 galaxies and perhaps the majority, the nuclear X-ray luminosity does not derive directly from the presence of a LLAGN.     

New active galactic nuclei among the INTEGRAL and SWIFT X-ray sources

We present the results of our optical identifications of a set of X-ray sources from the INTEGRAL and SWIFT all-sky surveys. The optical data have been obtained with the 1.5-m Russian-Turkish Telescope (RTT-150). Nine X-ray sources have been identified with active galactic nuclei (AGNs). Two of them are located in the nearby spiral galaxies MCG-01-05-047 and NGC 973 seen almost edge-on. One source, IGR J16562-3301, is probably a BL Lac object (blazar). The remaining AGNs are observed as the starlike nuclei of spiral galaxies whose spectra exhibit broad emission lines. The relation between the hard X-ray (17–60 keV) luminosity and the [O III] 5007 line luminosity, log L _x/L _{[O III]} ≈ 2.1, holds good for most of the AGNs detected in hard X rays. However, the luminosities of some AGNs deviate from this relation. The fraction of such objects can reach ～20%. In particular, the [O III] line flux is lower for two nearby edge-on spiral galaxies. This can be explained by the effect of absorption in the galactic disks.     

Spiral-Vortex Structure in the Gaseous Disks of Galaxies

General ideas, as well as experimental and theoretical efforts concerning the prediction and discovery of new structures in the disks of spiral galaxies – giant anticyclones - are reviewed. A crucial point is the development of a new method to restore the full vector velocity field of the galactic gas from the line-of-sight velocity field. This method can be used to get self-consistent solutions for the following problems: 1) determination of non-circular velocities associated with spiral-vortex structure; 2) determination of fundamental parameters of this structure: pattern speed, corotation radius, location of giant anticyclones; 3) refinement of galactic rotation curves taking into account regular non-circular motion in the spiral density wave, which makes it possible to build more accurate models of the mass distribution in the galaxy; 4) refinement of parameters of the rotating gaseous disk: inclination angle, center of rotation and position angle of the major dynamical axis, systematic velocity. The method is demonstrated using the restoration of the velocity field of the galaxy NGC 157 as an example. Results for this and some other spiral galaxies suggest that giant anticyclones are a universal property of galaxies with grand design structure.     

A ROSAT survey of Wolf–Rayet galaxies — II. The extended sample

We present results from an ongoing X-ray survey of Wolf–Rayet (WR) galaxies, a class of objects believed to be very young starbursts. This paper extends the first X-ray survey of WR galaxies by Stevens &38; Strickland by studying WR galaxies identified subsequent to the original WR galaxy catalogue of Conti.   Out of a sample of 40 new WR galaxies a total of 10 have been observed with the ROSAT PSPC, and of these seven have been detected (NGC 1365, NGC 1569, I Zw 18, NGC 3353, NGC 4449, NGC 5408 and a marginal detection of NGC 2366). Of these, all are dwarf starbursts except for NGC 1365, which is a barred spiral galaxy possibly with an active nucleus. We also report on observations of the related emission-line galaxy IRAS 0833+6517.   The X-ray properties of these galaxies are broadly in line with those found for the original sample; they are X-ray overluminous compared with their blue luminosity and have thermal spectra with typically kT  ∼ 0.4 − 1.0 keV. There are some oddities: NGC 5408 is very overluminous in X-rays, even compared with other WR galaxies; I Zw 18 has a harder X-ray spectrum; NGC 1365, although thought to contain an active nucleus, has X-ray properties that are broadly similar to other WR galaxies, and we suggest that the X-ray emission from NGC 1365 is due to starburst activity.   A good correlation between X-ray and blue luminosity is found for the WR galaxy sample as a whole. However, when just dwarf galaxies are considered there is little evidence of correlation. We discuss the implications of these results on our understanding of the X-ray emission from WR galaxies and suggest that the best explanation for the X-ray activity is starburst activity from a young starburst region.     

Observations and modelling of the interacting galaxies NGC 3395 and 3396

We present both observations and modelling of the atomic hydrogen in the closely interacting galaxies NGC 3395 and 3396. The observations were made with the VLA in both C- and D-arrays. We detect a large 'tail' of H  i extending to a projected distance of 63 kpc (10 arcmin) south-west of the pair, as well as two smaller galaxies, IC 2604, 14 arcmin to the south-west, and IC 2608, 14 arcmin to the south-east. However, these galaxies appear to have had at most a minor influence on the dynamics of NGC 3395/6. By means of N -body simulation we show that the tail is gas that has been stripped from NGC 3395 during a prograde encounter with NGC 3396, and that the pair has had one previous close approach. It is shown that the galaxies are within 5 × 10⁷ yr of their second perigalactic passage. Comparison of the time-scales for starburst activity with those from the simulations shows that the current starbursts are a result of the current close approach and not the first one. The interaction between NGC 3395 and 3396 has flattened the rotation curve of NGC 3396 owing to the parameters of the interaction. This naturally explains the more nucleated radio continuum structure observed in this galaxy, as significant infall and a subsequent central starburst would be expected in this scenario. The velocity structure and line profiles of the H  i are best explained if both cloud–cloud collisions between the two gas discs and tidal forces have been important.     

H i and CO observations of Arp 104: a spiral–elliptical interacting pair

We present data probing the spatial and kinematical distribution of both the atomic (H  i ) and molecular (CO) gas in NGC 5218, the late-type barred spiral galaxy in the spiral–elliptical interacting pair, Arp 104. We consider these data in conjunction with far-infrared and radio-continuum data, and N -body simulations, to study the galaxies interactions, and the star formation properties of NGC 5218. We use these data to assess the importance of the bar and tidal interaction on the evolution of NGC 5218, and the extent to which the tidal interaction may have been important in triggering the bar. The molecular gas distribution of NGC 5218 appears to have been strongly affected by the bar; the distribution is centrally condensed with a very large surface density in the central region. The N -body simulations indicate a time-scale since perigalacticon of  ∼3 × 10⁸ yr  , which is consistent with the interaction having triggered or enhanced the bar potential in NGC 5218, leading to inflow and the large central molecular gas density observed. Whilst NGC 5218 appears to be undergoing active star formation, its star formation efficiency is comparable to a 'normal' SBb galaxy. We propose that this system may be on the brink of a more active phase of star formation.     

Constraining the star formation histories of spiral bulges

Stellar populations in spiral bulges are investigated using the Lick system of spectral indices. Long-slit spectroscopic observations of line strengths and kinematics made along the minor axes of four spiral bulges are reported. Comparisons are made between central line strengths in spiral bulges and those in other morphological types [elliptical, spheroidal (Sph) and S0]. The bulges investigated are found to have central line strengths comparable to those of single stellar populations of approximately solar abundance or above. Negative radial gradients are observed in line strengths, similar to those exhibited by elliptical galaxies. The bulge data are also consistent with correlations between Mg₂, Mg₂ gradient and central velocity dispersion observed in elliptical galaxies. In contrast to elliptical galaxies, central line strengths lie within the loci defining the range of 〈Fe〉 and Mg₂ achieved by Worthey's solar abundance ratio, single stellar populations (SSPs). The implication of solar abundance ratios indicates significant differences in the star formation histories of spiral bulges and elliptical galaxies. A 'single zone with infall' model of galactic chemical evolution, using Worthey's SSPs, is used to constrain the possible star formation histories of our sample. We show that the 〈Fe〉, Mg₂ and H β line strengths observed in these bulges cannot be reproduced using primordial collapse models of formation but can be reproduced by models with extended infall of gas and star formation (2–17 Gyr) in the region modelled. One galaxy (NGC 5689) shows a central population with a luminosity-weighted average age of ∼5 Gyr, supporting the idea of extended star formation. Kinematic substructure, possibly associated with a central spike in metallicity, is observed at the centre of the Sa galaxy NGC 3623.     

The X-ray spectrum of NGC 7213 and the Seyfert–LINER connection

We present an XMM–Newton observation of the Seyfert–LINER (low-ionization nuclear emission-line region) galaxy NGC 7213. The RGS soft X-ray spectrum is well fitted with a power law plus soft X-ray collisionally ionized thermal plasma  ( kT = 0.18^+0.03_−0.01 keV)  . We confirm the presence of Fe  i , Fe  xxv and Fe  xxvi Kα emission in the EPIC spectrum and set tighter constraints on their equivalent widths of  82⁺¹⁰₋₁₃, 24⁺⁹₋₁₁  and 24⁺¹⁰₋₁₃ eV, respectively. We compare the observed properties together with the inferred mass accretion rate of NGC 7213 with those of other Seyfert and LINER galaxies. We find that NGC 7213 has intermediate X-ray spectral properties lying between those of the weak active galactic nucleus found in the LINER M81 and higher-luminosity Seyfert galaxies. There appears to be a continuous sequence of X-ray properties from the Galactic Centre through LINER galaxies to Seyferts, probably determined by the amount of material available for accretion in the central regions.     

Vertical distribution of Galactic disc stars and gas constrained by a molecular cloud complex

We investigate the dynamical effects of a molecular cloud complex with a mass ∼ 10⁷ M_⊙ and a size ∼ a few 100 pc on the vertical distribution of stars and atomic hydrogen gas in a spiral galactic disc. Such massive complexes have now been observed in a number of spiral galaxies. The extended mass distribution in a complex, with an average mass density 6 times higher than the Oort limit, is shown to dominate the local gravitational field. This results in a significant redistribution or clustering of the surrounding disc components towards the mid-plane, with a resulting decrease in their vertical scaleheights.
The modified, self-consistent stellar density distribution is obtained by solving the combined Poisson equation and the force equation along the z -direction for an isothermal stellar disc on which the complex is imposed. The effect of the complex is strongest at its centre, where the stellar mid-plane density increases by a factor of 2.6 and the vertical scaleheight decreases by a factor of 3.4 compared with the undisturbed stellar disc. A surprising result is the large radial distance of ∼ 500 pc from the complex centre over which the complex influences the disc; this is due to the extended mass distribution in a complex. The complex has a comparable effect on the vertical distribution of the atomic hydrogen gas in the galactic disc. This 'pinching' or constraining effect should be detectable in the nearby spiral galaxies, as for example has been done for NGC 2403 by Sicking. Thus the gravitational field of a complex results in local corrugations of the stellar and H  i vertical scaleheights, and the galactic disc potential is highly non-uniform on scales of the intercomplex separation of ∼ 1 kpc.     

Galaxy secular mass flow rate determination using the potential-density phase shift approach: Application to six nearby spiral galaxies

Using the potential-density phase shift approach developed by the present authors in earlier publications, we estimate the magnitude of radial mass accretion/excretion rates across the disks of six nearby spiral galaxies (NGC 628, NGC 3351, NGC 3627, NGC 4321, NGC 4736, and NGC 5194) having a range of Hubble types. Our goal is to examine these rates in the context of bulge building and secular morphological evolution along the Hubble sequence. Stellar surface density maps of the sample galaxies are derived from SINGS 3.6 μm and SDSS i-band images using colors as an indicator of mass-to-light ratios. Corresponding molecular and atomic gas surface densities are derived from published CO (1-0) and HI interferometric observations of the BIMA SONG, THINGS, and VIVA surveys. The mass flow rate calculations utilize a volume-type torque integral to calculate the angular momentum exchange rate between the basic state disk matter and what we assume to be density wave modes in the observed galaxies. This volume-type integral contains the contributions from both the gravitational surface torque couple and the advective surface torque couple at the nonlinear, quasi-steady state of the wave modes, in sharp contrast to its behavior in the linear regime, where it contains only the contribution from the gravitational surface torque couple used by Lynden-Bell & Kalnajs in 1972. The potential-density phase shift approach yields angular momentum transport rates several times higher than those estimated using the Lynden-Bell and Kalnajs approach. And unlike Lynden-Bell and Kalnajs, whose approach predicts zero mass redistribution across the majority of the disk surface (apart from the isolated locations of wave-particle resonances) for quasi-steady waves, the current approach leads to predictions of significant mass redistribution induced by the quasi-steady density wave modes, enough for the morphological types of disks to evolve substantially within its lifetime. This difference with the earlier conclusions of Lynden-Bell and Kalnajs reflects the dominant role played by collisionless shocks in the secular evolution of galaxies containing extremely non-linear, quasi-steady density wave modes, thus enabling significant morphological transformation along the Hubble sequence during a Hubble time. We show for the first time also, using observational data, that stellar mass accretion/excretion is just as important, and oftentimes much more important, than the corresponding accretion/excretion processes in the gaseous component, with the latter being what had been emphasized in most of the previous secular evolution studies.     

Study of the Structure of the Coma Cluster Based on a Hierarchical Force Clustering Method

Six subclusters in the Coma cluster have been selected on the basis of a hierarchical clustering method that takes the gravitational interaction among galaxies into account. Of these, 3 central subclusters around the galaxies NGC 4889, NGC 4874, and NGC 4839 have been singled out. We have used the objective statistical criterion applied by Vennik and Anosova in studies of close groups of galaxies to evaluate each member included in a subcluster with a high probability. Galaxies with a significant deficit of hydrogen HI, including objects from the Bravo-Alfaro list, have been identified with members of the subclusters, with the greatest number of them in the subclusters around NGC 4874 and NGC 4839. A quantitative estimate of the hydrogen deficit using the HI index in the RCG3 catalog reveals a statistically significant excess value for those galaxies that are members of the subclusters compared to galaxies with a hydrogen deficit in the overall Coma cluster field. A substantial number of the spiral galaxies with a hydrogen deficit in the subclusters turned out to be radio galaxies as well. __________ Translated from Astrofizika, Vol. 48, No. 4, pp. 515–528 (November 2005).     

Numerical studies on the structure of the cosmic ray electron halo in starburst galaxies

The structure of the cosmic ray electron halo of a starburst galaxy depends strongly on the nature of galactic wind and the configuration of the magnetic field. We have investigated these dependencies by solving numerically the propagation of electrons originating in starburst galaxies, most likely in supernova remnants. The calculations are made for several models for the galactic winds and for the configuration of the magnetic fields for comparison with observations. Our simulation of a quasi-radio halo reproduces both the extended structure of ∼9 kpc and the subtle hollow structure near the polar region of the radio halo that are observed in the starburst galaxy NGC 253. These findings suggest the existence of strong galactic wind in NGC 253.      

Tidal disruption model for the interacting pair of galaxies VV 117 (NGC 2444/45)

A series of N-body simulations have been performed to interpret the interacting pair of galaxies VV 117 (NGC 2444/45). The galaxies have been modelled assuming a mass ratio two. The simulations use various values for the distance of closest approach and the eccentricity of the relative orbit of the pair. A plausible scenario for the tidal disruption of the less massive galaxy is proposed. NGC2444, having double the mass of NGC 2445, has undergone penetrating collision with the latter in a hyperbolic or a parabolic orbit. After the first collision, the orbit has become bound. Our results show that VV 117 has either just emerged from the first collision or are on the verge of a second collision. NGC 2445suffers considerable disruption and mass loss. NGC 2444 is not affected much. The second collision is expected to culminate in the merger of the two galaxies. This revised version was published online in July 2006 with corrections to the Cover Date.     

Acceleration of cosmic-ray electrons by large-scale galactic spiral shocks — An observational report

The hypothesis that large-scale shocks accompanying spiral density waves are the main sources of cosmic-ray electrons for whole galactic disks has been tested utilizing results of 10.7 GHz observations of spiral galaxies. It is suggested that the 10.7 GHz emission of galactic disks is largely controlled by star-forming processes, and that large-scale shocks apparently play a minor role.     

Dark haloes of spiral galaxies: ISO photometry

We exclude hydrogen-burning stars, of any mass above the hydrogen-burning limit and any metallicity, as significant contributors to the massive haloes deduced from rotation curves to dominate the outer parts of spiral galaxies. We present and analyse images of four nearly edge-on bulgeless spiral galaxies (UGC 711, NGC 2915, UGC 12426, UGC 1459) obtained with ISOCAM (The CAMera instrument on board the Infrared Space Observatory ) at 14.5 and 6.75 μm. Our sensitivity limit for detection of any diffuse infrared emission associated with the dark haloes in these galaxies is a few tens of μJy per 6 × 6 arcsec² pixel, with this limit currently set by remaining difficulties in modelling the non-linear behaviour of the detectors. All four galaxies show zero detected signal from extended non-disc emission, consistent with zero halo-like luminosity density distribution. The 95 per cent upper limit on any emission, for NGC 2915 in particular, allows us to exclude very low mass main-sequence stars ( M  > 0.08 M⊙) and young brown dwarfs (≲1 Gyr) as significant contributors to dark matter in galactic haloes. Combining our results with those of the Galactic microlensing surveys, which exclude objects with M  < 0.01 M⊙, excludes almost the entire possible mass range of compact baryonic objects from contributing to Galactic dark matter.     

Vertical Scale Parameter Estimates for 48 Non-edge-on Spiral Galaxies

In the first paper of this series, we directly studied the mathematical forms, symmetry of spiral structure, and the projection of galactic discs on the images, and measured the pitch angles of the spiral arms and inclination angles of the galactic discs for 60 spiral galaxies. In this second paper, we estimate the vertical scale parameters of 48 non-edge-on spiral galaxies based on the method proposed by Peng et al. and on the results given in Paper Ⅰ. As we know, for edge-on disc galaxies we can obtain the vertical scale parameter from the photometry, once a mathematical form is specified for the vertical light distribution. For non-edgeon galaxies, some other methods have to be used. The statistical result was that the vertical scale parameter is comparable for edge-on and non-edge-on galaxies, although it is obtained from two very different methods.     

Simulations of stellar disks in galaxies with non-Newtonian interactions

We describe gravitationalN-body simulations to investigate whether various non-Newtonian interactions between the stars of a system could explain the flat rotational curves which are characteristic of actual isolated spiral galaxies. It is shown that replacing the standard Newtonian interaction by the models of Sanders (1984), Kuhn and Kruglyak (1987) and Milgrom (1983), no massive halo (or dark matter) is required to produce the flat rotational curves of the systems under consideration. All models also generate the exponential surface mass density distribution which is in agreement with that observed in disk-shaped galaxies. In relation to the spiral structure of galaxies, we present the evidence that the non-Newtonian interactions can reproduce the multiple armed patterns in stellar disks without dark matter.     

Exploring spiral galaxy potentials with hydrodynamical simulations

We study how well the complex gas velocity fields induced by massive spiral arms are modelled by the hydrodynamical simulations that we used recently to constrain the dark matter fraction in nearby spiral galaxies. More specifically, we explore the dependence of the positions and amplitudes of features in the gas flow on the temperature of the interstellar medium (assumed to behave as a one-component isothermal fluid), the non-axisymmetric disc contribution to the galactic potential, the pattern speed  Ω_p  , and finally the numerical resolution of the simulation. We argue that, after constraining the pattern speed reasonably well by matching the simulations to the observed spiral arm morphology, the amplitude of the non-axisymmetric perturbation (the disc fraction) is left as the primary parameter determining the gas dynamics. However, owing to the sensitivity of the positions of the shocks to modelling parameters, one has to be cautious when quantitatively comparing the simulations to observations. In particular, we show that a global least-squares analysis is not the optimal method for distinguishing different models, as it tends to slightly favour low disc fraction models. Nevertheless, we conclude that, given observational data of reasonably high spatial resolution and an accurate shock-resolving hydro-code, this method tightly constrains the dark matter content within spiral galaxies. We further argue that, even if the perturbations induced by spiral arms are weaker than those of strong bars, they are better suited for this kind of analysis because the spiral arms extend to larger radii where effects like inflows due to numerical viscosity and morphological dependence on gas sound speed are less of a concern than they are in the centres of discs.