Bar dissolution in non-spherical halos

Dynamical evolution of N-body bars embedded in spherical and prolate dark matter halos is investigated. In particular, the configuration such that galactic disks are placed in the plane perpendicular to the equatorial plane of the prolate halos is considered. Such a configuration is frequently found in cosmological simulations. N-body disks embedded in a fixed external halo potential were simulated, so that the barred structure was formed via dynamical instability in initially cool disks. In the subsequent evolution, bars in prolate halos dissolved gradually with time, while the bar pattern in spherical halos remained almost unchanged until the end of simulations. The e-folding time of bars suggest that they could be destroyed in a time smaller than a Hubble time. This revised version was published online in August 2006 with corrections to the Cover Date.     

Escape from planetary neighbourhoods

In this paper we use recently developed phase-space transport theory coupled with a so-called classical spectral theorem to develop a dynamically exact and computationally efficient procedure for studying escape from a planetary neighbourhood. The 'planetary neighbourhood' is a bounded region of phase space where entrance and escape are only possible by entering or exiting narrow 'bottlenecks' created by the influence of a saddle point. The method therefore immediately applies to, for example, the circular restricted three-body problem and Hill's lunar problem (which we use to illustrate the results), but it also applies to more complex, and higher-dimensional, systems possessing the relevant phase-space structure. It is shown how one can efficiently compute the mean passage time through the planetary neighbourhood, the phase-space flux in, and out, of the planetary neighbourhood, the phase-space volume of initial conditions corresponding to trajectories that escape from the planetary neighbourhood, and the fraction of initial conditions in the planetary neighbourhood corresponding to bound trajectories. These quantities are computed for Hill's problem. We study the dependence of the proportions of these quantities on energy and dimensionality (two-dimensional planar and three-dimensional spatial Hill's problem). The methods and quantities presented are of central interest for many celestial and stellar dynamical applications such as, for example, the capture and escape of moons near giant planets, the formation of binaries in the Kuiper belt and the escape of stars from star clusters orbiting about a galaxy.     

Anisotropic Inhomogeneous Riemann S-Ellipsoids in Spheroidal Halos. I.

The classical Riemann S-ellipsoids are generalized to account for a Ferrers inhomogeneity in the mass distribution, a three-dimensional anisotropy of the velocity dispersions, and the gravitation of a spheroidal halo. The Ferrers inhomogeneity does not affect the equilibrium and stability conditions for the ellipsoids, but only changes the numerical coefficients in the equilibrium parameters. The anisotropy of the velocity dispersions changes both the equilibrium conditions and the criterion for stability of the ellipsoids, either extending or limiting the region in which the ellipsoids exist. A halo extends this region in all cases and stabilizes the ellipsoids. The observed of the velocity dispersions of stars in the neighborhood of the sun are obtained using the inhomgeneous model. __________ Translated from Astrofizika, Vol. 48, No. 4, pp. 613–632 (November 2005).     

Correlations, spectra, and instability of phase-space density fluctuations in open-cluster models

The dynamical evolution of six open star cluster models is analyzed using the correlation and spectral analysis of phase-space density fluctuations. The two-time and mutual correlation functions are computed for the fluctuations of the phase-space density of cluster models. The data for two-time and two-particle correlations are used to determine the correlation time for phase-space density fluctuations ((0.1–1) τ _v.r., where τ _v.r. is the violent relaxation time of the model) and the average phase velocities of the propagation of such fluctuations in cluster models. These velocities are 2–20 times smaller than the root mean square velocities of the stars in the cluster core. The power spectra and dispersion curves of phase-space density fluctuations are computed using the Fourier transform of mutual correlation functions. The results confirm the presence of known unstable phase-space density fluctuations due to homologous fluctuations of the cluster cores. The models are found to exhibit a number of new unstable phase-space density fluctuations (up to 32–41 pairs of fluctuations with different complex conjugate frequencies in each model; the e-folding time of the amplitude growth of such fluctuations is (0.4–10) τ _v.r. and their phases are distributed rather uniformly). Astrophysical applications of the obtained results (irregular structure of open star clusters, formation and decay of quasi-stationary states in such clusters) are discussed.     

Anisotropic and inhomogeneous S-type Riemann ellipsoids inside spheroidal halos. II

Classical S-type Riemann ellipsoids are generalized taking a Ferrers-type inhomogeneity in the mass distribution, a three-dimensional anisotropy in the dispersion in the velocities, and the gravitation of a spherical halo into account. A Ferrers inhomogeneity in the mass distribution does not affect the equilibrium conditions for the ellipsoids, but only changes the numerical coefficients in the equilibrium parameters. An anisotropy in the dispersion in the velocities does change the equilibrium conditions, expanding or contracting the region in which embedded ellipsoids exist. A halo extends this region in all cases. __________ Translated from Astrofizika, Vol. 49, No. 3, pp.359–373 (August 2006).     

Light subsystem with a linear velocity field inside a gravitating sphere

We have obtained possible figures of equilibrium for a light subsystem with internal material flows of constant vorticity inside a gravitating homogeneous sphere. Along with ellipsoidal figures of equilibrium (light analogs of Riemann ellipsoids), we have obtained figures of equilibrium in the form of one-sheet and two-sheeted hyperboloids undergoing tilted rotation.Translated from Astrofizika, Vol. 38, No. 1, pp. 55–73, January–March, 1995.     

Signatures of galaxy mergers in the Milky Way: Here, there and everywhere

The hierarchical paradigm predicts that large galaxies like the Milky Way formed through mergers of smaller systems, which are expected to leave behind substructures in the halo of the final product. Recently the first tests of this prediction on galaxies other than the Milky Way have been made (eg NGC 5907 by Sackett et al.; M31 by Ibata et al.), but one should bare in mind that it is extremely difficult to detect halos in external galaxies let alone substructures in those halos. On the other hand, the multi-dimensional phase-space information available for our Galaxy (6d for stars in the vicinity of the Sun, and 4d for more distant ones) enables us to directly search for merger signatures. This revised version was published online in September 2006 with corrections to the Cover Date.     

Evolution of Polytropic Reimann Disks in Halos

The equilibrium of elliptical Riemann disks with a polytropic equation of state and their evolution under the influence of viscosity and gravitational radiation inside spheroidal halos with a relative surface mass density k. The evolutionary trajectory of a disk inside a halo with k<0.5, which is analogous to the evolution of an isolated disk, differs from that of a disk inside a denser halo.     

Non-zero rest mass neutrinos and the cosmological constant

Recently it was pointed out that a non-zero cosmological constant can play a role in the formation of neutrino halos only in the case of neutrinos of very low rest mass (m _v <-0.1eV). However, phase-space considerations would requirem _v >50 eV if neutrinos dominate the missing mass in halos of large spiral galaxies and moreoverm _v >200 eV is implied in the case of dwarf spheroidals. These larger neutrino masses would be in conflict with observed constraints on the age of the Universe unless a cosmological constant is invoked.     

Galactic models in the light of collision dynamics

A simple and flexible model of a spiral galaxy is developed in the light of the study of the dynamics of colliding galaxies. As such interactions are strongly characterised by the binding energy distributions of the systems involved; hence, the main criterion used in constructing the model is its compatibility with the observed shapes and intensity distribution of light of spirals. The model basically consists of thick exponential disk component and a spherical polytropic bulge component. The relevant quantities pertaining to the model, significant for dynamical studies, are computed and analysed.The study of the dynamics of colliding galaxies require galactic models with and without halos, depending upon the problem under study and the effects to be accessed. Consequently, we have preferred to tackle the problem of galactic models with massive halos separately in the next paper.     

Stability of a Riemann S Ellipsoid with a Spheroidal Halo

The stability of Riemann S ellipsoids inside an oblate halo with respect to the second form of oscillations is investigated. It is shown that some ellipsoids with reverse internal circulation of matter, which are stable inside a spherical halo or in its absence, become unstable with respect to second odd forms of oscillation inside an oblate halo. Here there is asymmetry between conjugate ellipsoids from the standpoint of their stability. Only those conjugate ellipsoids that correspond to higher frequencies of reverse circulation of matter than their corresponding angular rotation rates are unstable. The domains of instability of light and heavy conjugate embedded ellipsoids are obtained as a function of the oblateness measure and relative density of the halo.     

MOND plus classical neutrinos are not enough for cluster lensing

Clusters of galaxies offer a robust test bed for probing the nature of dark matter that is insensitive to the assumption of the gravity theories. Both Modified Newtonian Dynamics (MOND) and General Relativity (GR) would require similar amounts of non-baryonic matter in clusters as MOND boosts the gravity only mildly on cluster scales. Gravitational lensing allows us to estimate the enclosed mass in clusters on small (∼20–50 kpc) and large (∼several 100 kpc) scales independent of the assumptions of equilibrium. Here, we show for the first time that a combination of strong and weak gravitational lensing effects can set interesting limits on the phase-space density of dark matter in the centres of clusters. The phase-space densities derived from lensing observations are inconsistent with neutrino masses ranging from 2–7 eV, and hence do not support the 2 eV-range particles required by MOND. To survive, the most plausible modification for MOND may be an additional degree of dynamical freedom in a covariant incarnation.     

The Connection between Radio Halos and Cluster Mergers and the Statistical Properties of the Radio Halo Population

We discuss the statistical properties of the radio halo population in galaxy clusters. Radio bi-modality is observed in galaxy clusters: a fraction of clusters host giant radio halos while a majority of clusters do not show evidence of diffuse cluster-scale radio emission. The radio bi-modality has a correspondence in terms of dynamical state of the hosting clusters showing that merging clusters host radio halos and follow the well-known radio-X-ray correlation, while more relaxed clusters do not host radio halos and populate a region well separated from that correlation. These evidences can be understood in the framework of a scenario where merger-driven turbulence re-accelerate the radio emitting electrons. We discuss the main statistical expectations of this scenario underlining the important role of the upcoming LOFAR surveys to test present models.     

Varying Fundamental Constants: A Dynamical Systems Approach

Theory of dynamical systems offers a possibility of investigating the space of all possible solutions. In the context of simple cosmological models such like Varying Speed of Light Friedman-Robertson-Walker (VSL FRW) models there exists a systematic method of reducing field equations to certain two-dimensional dynamical system. One of the features of this reduction is the possibility of representing the model as a Hamiltonian system in which the properties of the potential function V(X) can serve as a tool for qualitative classification of possible evolutions of a(t). Some important features like resolution of the flatness problem, existence of event horizons near the singularity can be visualized as domains on the phase-space. Then one is able to see how large is the class of solutions (labelled by the initial conditions) leading to the desired property. This revised version was published online in July 2006 with corrections to the Cover Date.     

Interacting galaxies and mergers

A study of the merger time-scales of various types of interacting galaxies is conducted on the basis of the collisional theory. The results indicate that in the absence of halos, violently interacting galaxies merge in a time-scale of ~ 10⁸ years; but the mildly interacting ones have merger time-scales from ~ 10⁹ to 10¹⁰ years. However, in the presence of halos, all types of interacting galaxies are likely to merge in a time-scale of 10⁸ years (as indicated by preliminary calculations). Galaxy evolution by mutual interactions is likely to have its reflection on the fundamental plane, as during the process the dynamical structures of the progenitors change and dissipation occurs.     

Self-consistent Modelling of the interstellar medium

The dynamical evolution of hot optically thin plasmas in the ISMcrucially depends on the heating and cooling processes. It isessential to realize that all physical processes that contributeoperate on different time scales. In particular detailedbalancing is often violated since the statistically inverseprocess of e.g. collisional ionization is recombination of an ionwith two electrons, which as a three-body collision is usuallydominated by radiative recombination, causing a departure fromcollisional ionization equilibrium. On top of these differences inatomic time scales, hot plasmas are often in a dynamical state,thereby introducing another time scale, which canbe the shortest one.The non-equilibrium effects will be illustrated and discussed inthe case of galactic outflows. It will be shown, that spectralanalyses of X-ray data of edge-on galaxies show a clear signaturein the form of ‘multi-temperature’ halos, which can mostnaturally be explained by the ‘freezing-in’ of highly ionizedspecies in the outflow, which contribute to the overall spectrumby delayed recombination. This naturally leads to anon-equilibrium cooling function, which modifies the dynamics,which in turn changes the plasma densities and thermal energybudget, thus feeding back on the ionization structure. Thereforeself-consistent modelling is needed.     

Gravitational Mesoscopic Constraints in Cosmological Dark Matter Halos

We present an analysis of the behaviour of the ‘coarse-grained’ (‘mesoscopic’) rank partitioning of the mean energy of collections of particles composing virialized dark matter halos in a Λ-CDM cosmological simulation. We find evidence that rank preservation depends on halo mass, in the sense that more massive halos show more rank preservation than less massive ones. We find that the most massive halos obey Arnold’s theorem (on the ordering of the characteristic frequencies of the system) more frequently than less massive halos. This method may be useful to evaluate the coarse-graining level (minimum number of particles per energy cell) necessary to reasonably measure signatures of ‘mesoscopic’ rank orderings in a gravitational system.     

Does Deuterium Enable the Formation of Primordial Brown Dwarfs?

The time evolution of barred structures is examined under the influence of the external forces exerted by a spherical halo and by prolate halos. In particular, galaxy disks are placed in the plane including the major axis of prolate halos, whose configuration is often found in cosmological simulations. N-body disks in fixed external halo fields are simulated, so that bars are formed via dynamical instability. In the subsequent evolution, the bars in prolate halos dissolve gradually with time, while the bar pattern in a spherical halo remains almost unchanged to the end of the simulation. The decay times of the bars suggest that they can be destroyed in a time smaller than a Hubble time. Our results indicate that this dissolution process could occur in real barred galaxies, if they are surrounded by massive dark prolate halos, and the configuration adopted here is not unusual from the viewpoint of galaxy formation. For a prolate halo model, an additional simulation that is restricted to two-dimensional in-plane motions has also ended up with similar bar dissolution. This means that the vertical motions of disk stars do not play an essential role in the bar dissolution demonstrated here.     

A simple approach to the high-redshift sub-millimeter galaxies

We aim at understanding the statistical properties of luminous sub-millimeter (submm) galaxies (SMGs) in the context of cosmological structure formation. By utilizing a cosmological N-body simulation to calculate the distribution of dark halos in the Universe, we consider the dust enrichment in individual halos by Type II supernovae (SNe II). The SN II rate is estimated under a star formation activity which is assumed to occur on a dynamical timescale in the dark matter potential. Our simple framework successfully explains the luminosity function, the typical star formation rate, and the typical dust mass of an observational SMG sample at z～3. We also examine the clustering properties of SMGs, since a positive cross correlation between SMGs and Lyα emitters (LAEs) is indeed observed by a recent observation. In the simulation, we select SMGs by FIR dust luminosity >10¹² L _⊙, while LAEs are chosen such that the age and the virial mass are consistent with the observed LAE properties. The SMGs and LAEs selected in this way show a spatial cross correlation whose strength is consistent with the observation. This confirms that the SMGs really trace the most clustered regions at z～3 and that their luminosities can be explained by the dust accumulation as a result of their star formation activities. We extend our prediction to higher redshifts, finding that a statistical sample of submm galaxies at z≥6 can be obtained by ALMA with a 100 arcmin² survey. With the same survey, a few submm galaxies at z～10 may be detected.     

Prompt Iron Enrichment, Two r-Process Components, and Abundances in Very Metal-Poor Stars

We compare the structure and substructure of dark matter halos in model universes dominated by collisional, strongly self-interacting dark matter (SIDM) and collisionless, weakly interacting dark matter (CDM). While SIDM virialized halos are more nearly spherical than CDM halos, they can be rotationally flattened by as much as 20% in their inner regions. Substructure halos suffer ram-pressure truncation and drag, which are more rapid and severe than their gravitational counterparts tidal stripping and dynamical friction. Lensing constraints on the size of galactic halos in clusters are a factor of 2 smaller than predicted by gravitational stripping, and the recent detection of tidal streams of stars escaping from the satellite galaxy Carina suggests that its tidal radius is close to its optical radius of a few hundred parsecs-an order of magnitude smaller than predicted by CDM models but consistent with SIDM models. The orbits of SIDM satellites suffer significant velocity bias, sigmaSIDM&solm0;sigmaCDM=0.85, and are more circular than CDM satellites, betaSIDM approximately 0.5, in agreement with the inferred orbits of the Galaxy's satellites. In the limit of a short mean free path, SIDM halos have singular isothermal density profiles; thus, in its simplest incarnation SIDM, is inconsistent with galactic rotation curves.