In weakly collisional stellar systems such as some globular clusters, partial energy equipartition and mass segregation are expected to develop as a result of the cumulative effect of stellar encounters even in systems initially characterized by star-mass independent density and energy distributions. In parallel, numerical simulations have demonstrated that radially-biased pressure anisotropy slowly builds up in realistic models of globular clusters from initial isotropic conditions, leading to anisotropy profiles that, to some extent, mimic those resulting from incomplete violent relaxation known to be relevant to elliptical galaxies. In this paper we consider a set of realistic simulations realized by means of Monte Carlo methods and analyze them by means of self-consistent two-component models. For the purpose, we refer to an underlying distribution function, originally conceived to describe elliptical galaxies, that has been recently truncated and adapted to the context of globular clusters. The two components are meant to represent light stars (combining all main sequence stars) and heavy stars (giants, dark remnants, and binaries). We show that this conceptually simple family of two-component truncated models provides a reasonable description of simulated density, velocity dispersion, and anisotropy profiles, especially for the most relaxed systems, with the capability to express quantitatively the attained levels of energy equipartition and mass segregation. In contrast, two-component isotropic models based on the King distribution function do not offer a comparably satisfactory representation of the simulated globular clusters. With this work we provide a new reliable diagnostic tool applicable to nonrotating globular clusters that are characterized by significant gradients in the local value of the mass-to-light ratio, beyond the commonly used one-component dynamical models.

As a result of the slow action of two-body encounters, globular clusters develop mass segregation and attain a condition of only partial energy equipartition even in their central, most relaxed regions. Realistic numerical simulations show that, during the process, a radially-biased anisotropy profile slowly builds up, mimicking that resulting from incomplete violent relaxation. Commonly used dynamical models, such as the one-component King models, cannot describe these properties. Here we show that simple two-component models based on a distribution function originally conceived to describe elliptical galaxies, recently truncated and adapted to the context of globular clusters, can describe in detail what is observed in complex and realistic numerical simulations.

With the help of high-resolution long-slit and integral-field spectroscopy observations, the number of confirmed cases of galaxies with counterrotation is increasing rapidly. The evolution of such counterrotating galaxies remains far from being well understood. In this paper we study the dynamics of counterrotating collisionless stellar disks by means of $N$-body simulations. We show that, in the presence of counterrotation, an otherwise gravitationally stable disk can naturally generate bending waves accompanied by strong disk heating across the disk plane, that is in the vertical direction. Such conclusion is found to hold even for dynamically warm systems with typical values of the initial vertical-to-radial velocity dispersion ratio $\sigma_{\rm z}/\sigma_{\rm R} \approx 0.5$, for which the role of pressure anisotropy should be unimportant. We note that, during evolution, the $\sigma_{\rm z}/\sigma_{\rm R}$ ratio tends to rise up to values close to unity in the case of locally Jeans-stable disks, whereas in disks that are initially Jeans-unstable it may reach even higher values, especially in the innermost regions. This unusual behavior of the $\sigma_{\rm z}/\sigma_{\rm R}$ ratio in galaxies with counterrotation appears not to have been noticed earlier. Our investigations of systems made of two counterrotating components with different mass-ratios suggest that even apparently normal disk galaxies (i.e., with a minor counterrotating component so as to escape detection in current observations) might be subject to significant disk heating especially in the vertical direction.

Recently, a class of non-truncated radially-anisotropic models (the so-called $f^{(\nu)}$-models), originally constructed in the context of violent relaxation and modeling of elliptical galaxies, has been found to possess interesting qualities in relation to observed and simulated globular clusters. In view of new applications to globular clusters, we improve this class of models along two directions. To make them more suitable for the description of small stellar systems hosted by galaxies, we introduce a 'tidal' truncation (by means of a procedure that guarantees full continuity of the distribution function). The new $f_T^{(\nu)}$-models are shown to provide a better fit to the observed photometric and spectroscopic profiles for a sample of 13 globular clusters studied earlier by means of non-truncated models; interestingly, the best-fit models also perform better with respect to the radial-orbit instability. Then we design a flexible but simple two-component family of truncated models, to study the separate issues of mass segregation and of multiple populations. We do not aim at a fully realistic description of globular clusters, to compete with the description currently obtained by means of dedicated simulations. The goal here is to try to identify the simplest models, that is, those with the smallest number of free parameters, but still able to provide a reasonable description for clusters that are evidently beyond the reach of one-component models: with this tool we aim at identifying the key factors that characterize mass segregation or the presence of multiple populations. To reduce the relevant parameter space, we formulate a few physical arguments (based on recent observations and simulations). A first application to two well-studied globular clusters is briefly described and discussed.

By means of 3D hydrodynamical simulations, in a separate paper we have discussed the properties of non-axisymmetric density wave trains in the outermost regions of galaxy disks, based on the picture that self-excited global spiral modes in the bright optical stellar disk are accompanied by low-amplitude short trailing wave signals outside corotation; in the gas, such wave trains can penetrate through the outer Lindblad resonance and propagate outwards, forming prominent spiral patterns. In this paper we present the synthetic 21~cm velocity maps expected from simulated models of the outer gaseous disk, focusing on the case when the disk is dominated by a two-armed spiral pattern, but considering also other more complex situations. We discuss some aspects of the spiral pattern in the gaseous periphery of galaxy disks noted in our simulations that might be interesting to compare with specific observed cases.

Deep HI observations of the outer parts of disc galaxies demonstrate the frequent presence of extended, well-developed spiral arms far beyond the optical radius. To understand the nature and the origin of such outer spiral structure, we investigate the propagation in the outer gaseous disc of large-scale spiral waves excited in the bright optical disc. Using hydrodynamical simulations, we show that non-axisymmetric density waves, penetrating in the gas through the outer Lindblad resonance, can exhibit relatively regular spiral structures outside the bright optical stellar disc. For low-amplitude structures, the results of numerical simulations match the predictions of a simple WKB linear theory. The amplitude of spiral structure increases rapidly with radius. Beyond $\approx 2$ optical radii, spirals become nonlinear (the linear theory becomes quantitatively and qualitatively inadequate) and unstable to Kelvin-Helmholtz instability. In numerical simulations, in models for which gas is available very far out, spiral arms can extend out to 25 disc scale-lengths. A comparison between the properties of the models we have investigated and the observed properties of individual galaxies may shed light into the problem of the amount and distribution of dark matter in the outer halo.

Internal rotation is considered to play a major role in the dynamics of some globular clusters. However, in only few cases it has been studied by quantitative application of realistic and physically justified global models. Here we present a dynamical analysis of the photometry and three-dimensional kinematics of omega Cen, 47 Tuc, and M15, by means of a recently introduced family of self-consistent axisymmetric rotating models. The three clusters, characterized by different relaxation conditions, show evidence of differential rotation and deviations from sphericity. The combination of line-of-sight velocities and proper motions allows us to determine their internal dynamics, predict their morphology, and estimate their dynamical distance. The well-relaxed cluster 47 Tuc is very well interpreted by our model; internal rotation is found to explain the observed morphology. For M15, we provide a global model in good agreement with the data, including the central behavior of the rotation profile and the shape of the ellipticity profile. For the partially relaxed cluster omega Cen, the selected model reproduces the complex three-dimensional kinematics; in particular the observed anisotropy profile, characterized by a transition from isotropy, to weakly-radial anisotropy, and then to tangential anisotropy in the outer parts. The discrepancy found for the steep central gradient in the observed line-of-sight velocity dispersion profile and for the ellipticity profile is ascribed to the condition of only partial relaxation of this cluster and the interplay between rotation and radial anisotropy.

We present new number density and surface brightness profiles for the globular cluster M92 (NGC 6341). These profiles are calculated from optical images collected with the CCD mosaic camera MegaCam at the Canada-France-Hawaii-Telescope and with the Advanced Camera for Surveys on the Hubble Space Telescope. The ground-based data were supplemented with the Sloan Digital Sky Survey photometric catalog. Special care was taken to discriminate candidate cluster stars from field stars and to subtract the background contamination from both profiles. By examining the contour levels of the number density, we found that the stellar distribution becomes clumpy at radial distances larger than about 13 arcminutes, and there is no preferred orientation of contours in space. We performed detailed fits of King and Wilson models to the observed profiles. The best-fit models underestimate the number density inside the core radius. Wilson models better represent the observations, in particular in the outermost cluster regions: the good global agreement of these models with the observations suggests that there is no need to introduce an extra-tidal halo to explain the radial distribution of stars at large radial distances. The best-fit models for the number density and the surface brightness profiles are different, even though they are based on the same observations. Additional tests support the evidence that this fact reflects the difference in the radial distribution of the stellar tracers that determine the observed profiles (main sequence stars for the number density, bright evolved stars for the surface brightness).

Internal rotation is considered to play a major role in determining the structure and dynamics of some globular clusters. We present a dynamical analysis of the photometry and three-dimensional kinematics of 47 Tuc and omega Cen, by means of a new family of self-consistent axisymmetric rotating models. The combined use of line-of-sight velocities and proper motions allows us to obtain a global description of the internal dynamical structure of the objects together with an estimate of their dynamical distances. The well-relaxed cluster 47 Tuc is very well interpreted by our dynamical models; in particular, internal rotation is found to explain the observed morphology. For the partially relaxed cluster omega Cen, the selected model provides a good representation of its complex three-dimensional kinematics, in general qualitative agreement with the observed anisotropy profile, which is characterized by tangential anisotropy in the outer parts; discrepancies are found between the observed and the expected ellipticity profile and are ascribed to the presence of a high degree of radial anisotropy in the intermediate region and to its interplay with rotation.

A re-investigation of the gravothermal catastrophe is presented. By means of a linear perturbation analysis, we study the dynamical stability of a spherical self-gravitating isothermal fluid of finite volume and find that the conditions for the onset of the gravothermal catastrophe, under different external conditions, coincide with those obtained from thermodynamical arguments. This suggests that the gravothermal catastrophe may reduce to Jeans instability, rediscovered in an inhomogeneous framework. We find normal modes and frequencies for the fluid system and show that instability develops on the dynamical time scale. We then discuss several related issues. In particular: (1) For perturbations at constant total energy and constant volume, we introduce a simple heuristic term in the energy budget to mimic the role of binaries. (2) We outline the analysis of the two-component case and show how linear perturbation analysis can be carried out also in this more complex context in a relatively straightforward way. (3) We compare the behavior of the fluid model with that of the collisionless sphere. In the collisionless case the instability seems to disappear, which is at variance with the linear Jeans stability analysis in the homogeneous case; we argue that a key ingredient to understand the difference (a spherical stellar system is expected to undergo the gravothermal catastrophe only in the presence of some collisionality, which suggests that the instability is dissipative and not dynamical) lies in the role of the detailed angular momentum in a collisionless system. Finally, we briefly comment on the meaning of the Boltzmann entropy and its applicability to the study of the dynamics of self- gravitating inhomogeneous gaseous systems.

Two new families of self-consistent axisymmetric truncated equilibrium models for the description of quasi-relaxed rotating stellar systems are presented. The first extends the spherical King models to the case of solid-body rotation. The second is characterized by differential rotation, designed to be rigid in the central regions and to vanish in the outer parts, where the energy truncation becomes effective. The models are constructed by solving the nonlinear Poisson equation for the self-consistent mean-field potential. For rigidly rotating configurations, the solutions are obtained by an asymptotic expansion on the rotation strength parameter. The differentially rotating models are constructed by means of an iterative approach based on a Legendre series expansion of the density and the potential. The two classes of models exhibit complementary properties. The rigidly rotating configurations are flattened toward the equatorial plane, with deviations from spherical symmetry that increase with the distance from the center. For models of the second family, the deviations from spherical symmetry are strongest in the central region, whereas the outer parts tend to be quasi-spherical. The relevant parameter spaces are explored and the intrinsic and projected structural properties are described. Special attention is given to the effect of different options for the truncation of the distribution function in phase space. Models in the moderate rotation regime are best suited to applications to globular clusters. For general interest in stellar dynamics, at high values of the rotation strength the differentially rotating models exhibit a toroidal core embedded in a quasi-spherical configuration. Physically simple analytical models of the kind presented here provide insights into dynamical mechanisms and may be a basis for more realistic investigations with the help of N-body simulations.

We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic f^(\nu) models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the first time to the study of globular clusters. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, f^(\nu) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of accurate kinematic data-points, well distributed over the cluster field.

Weak lensing is one of the best available diagnostic tools to measure the total density profiles of distant clusters of galaxies. Unfortunately, it suffers from the well-known mass-sheet degeneracy, so that weak lensing analyses cannot lead to fully reliable determinations of the total mass of the clusters. One possible way to set the relevant scale of the density profile would be to make a direct measurement of the magnification produced by the clusters as gravitational lenses. In this paper we revisit a suggestion made a few years ago for this general purpose, based on the use of the Fundamental Plane as a standard rod for early-type galaxies. Here we move one step further, beyond the simple outline of the idea given earlier, and quantify some statistical properties of this innovative diagnostic tool, with the final goal of identifying clear guidelines for a future observational test on concrete cases, which turns out to be well within the current instrument capabilities. The study is carried out by discussing the statistical properties of Fundamental Plane measurements for a sample of early-type source galaxies behind a massive cluster, for which a weak lensing analysis is assumed to be available. We proceed to study the best strategy to use Fundamental Plane measurements to determine the mass scale of a given cluster and find that the optimal choice is that of a sample of early-type galaxies behind the cluster distributed approximately uniformly in the sky. We also show that substructures do not contribute much to the magnification signal that we are looking for, but only add a modest amount of scatter. We find that for a massive cluster (M > 10^15 M_Sun) located at redshift 0.3 +/- 0.1, a set of about 20 Fundamental Plane measurements, combined with a good weak lensing analysis, should be able to lead to a mass determination with a precision of 20% or better.

Disk-halo decompositions of galaxy rotation curves are generally performed in a parametric way. We construct self-consistent models of nonspherical isothermal halos embedding a zero-thickness disk, by assuming that the halo distribution function is a Maxwellian. The method developed here can be used to study other physically-based choices for the halo distribution function and the case of a disk accompanied by a bulge. In a preliminary investigation we note the existence of a fine tuning between the scalelengths R_{\Omega} and h, respectively characterizing the rise of the rotation curve and the luminosity profile of the disk, which surprisingly applies to both high surface brightness and low surface brightness galaxies. This empirical correlation identifies a much stronger conspiracy than the one required by the smoothness and flatness of the rotation curve (disk-halo conspiracy). The self-consistent models are characterized by smooth and flat rotation curves for very different disk-to-halo mass ratios, hence suggesting that conspiracy is not as dramatic as often imagined. For a typical rotation curve, with asymptotically flat rotation curve at V_{\infty} (the precise value of which can also be treated as a free parameter), and a typical density profile of the disk, self-consistent models are characterized by two dimensionless parameters, which correspond to the dimensional scales (the disk mass-to-light ratio M/L and the halo central density) of standard disk-halo decompositions. We show that if the rotation curve is decomposed by means of our self-consistent models, the disk-halo degeneracy is removed and typical rotation curves are fitted by models that are below the maximum-disk prescription. Similar results are obtained from a study of NGC 3198. Finally, we quantify the flattening of the spheroidal halo, which is significant, especially on the scale of the visible disk.

In a previous paper we found that the Globular Clusters of our Galaxy lie around a line in the log(Re), SBe, log(sigma) parameter space, with a moderate degree of scatter and remarkable axi-symmetry. This implies the existence of a purely photometric scaling law obtained by projecting such a line onto the log(Re), SBe plane. Such photometric quantities are readily available for large samples of clusters, as opposed to stellar velocity dispersion data. We study a sample of 129 Galactic and extragalactic clusters on such photometric plane in the V-band. We look for a linear relation between SBe and log(Re) and study how the scatter around it is influenced by age and dynamical environment. We interpret our results as a test on the evolutionary versus primordial origin of the Fundamental Line. We perform a detailed analysis of surface brightness profiles, which allows us to present a catalogue of structural properties, without relying on a given dynamical model. We find a linear relation between SBe and log(Re), in the form SBe = (5.25 +- 0.44) log(Re) + (15.58 +- 0.28), where SBe is measured in mag/arcsec^2 and Re in parsec. Both young and old clusters lie on the scaling law, with a scatter of approximately 1 mag in SBe. However, young clusters display more scatter and a clear trend of such scatter with age, which old clusters do not. Such trend becomes tighter if cluster age is measured in units of the cluster half-light relaxation time. Two-body relaxation therefore plays a major role, together with passive stellar population evolution, in shaping the relation between SBe, log(Re), and cluster age. We argue that the log(Re)-SBe relation and hence the Fundamental Line scaling law is not primordially set at cluster formation, but rather is the result of combined stellar evolution and collisional dynamical evolution.

Context: Several spiral galaxies, as beautifully exhibited by the case of NGC 6946, display a prominent large-scale spiral structure in their gaseous outer disk. Such structure is often thought to pose a dynamical puzzle, because grand-design spiral structure is traditionally interpreted as the result of density waves carried mostly in the stellar disk. Aims. Here we argue that the outer spiral arms in the cold gas outside the bright optical disk actually have a natural interpretation as the manifestation of the mechanism that excites grand-design spiral structure in the main, star-dominated body of the disk: the excitation is driven by angular momentum transport to the outer regions, through trailing density waves outside the corotation circle that can penetrate beyond the Outer Lindblad Resonance in the gaseous component of the disk. Methods: Because of conservation of the density wave action, these outgoing waves are likely to become more prominent in the outer disk and eventually reach non-linear amplitudes. To calculate the desired amplitude profiles, we make use of the theory of dispersive waves. Results: If the conditions beyond the optical radius allow for an approximate treatment in terms of a linear theory, we show that fitting the observed amplitude profiles leads to a quantitative test on the density of the disk material and thus on the dark matter distribution in the outer parts of the galaxy. Conclusions: This study is thus of interest to the general problem of the disk-halo decomposition of rotation curves.

In a previous paper, we have constructed a family of self-consistent triaxial models of quasi-relaxed stellar systems, shaped by the tidal field of the hosting galaxy, as an extension of the well-known spherical King models. For a given tidal field, the models are characterized by two physical scales (such as total mass and central velocity dispersion) and two dimensionless parameters (the concentration parameter and the tidal strength). The most significant departure from spherical symmetry occurs when the truncation radius of the corresponding spherical King model is of the order of the tidal radius, which, for a given tidal strength, is set by the maximum concentration value admitted. For such maximally extended (or "critical") models the outer boundary has a generally triaxial shape, given by the zero-velocity surface of the relevant Jacobi integral, which is basically independent of the concentration parameter. In turn, the external tidal field can give rise to significant global departures from spherical symmetry (as measured, for example, by the quadrupole of the mass distribution of the stellar system) only for low-concentration models, for which the allowed maximal value of the tidal strength can be relatively high. In this paper we describe in detail the intrinsic and the projected structure and kinematics of the models, covering the entire parameter space, from the case of sub-critical (characterized by "underfilling" of the relevant Roche volume) to that of critical models. The intrinsic properties can be a useful starting point for numerical simulations and other investigations that require initialization of a stellar system in dynamical equilibrium. The projected properties are a key step in the direction of a comparison with observed globular clusters and other candidate stellar systems.

Spherical models of collisionless but quasi-relaxed stellar systems have long been studied as a natural framework for the description of globular clusters. Here we consider the construction of self-consistent models under the same physical conditions, but including explicitly the ingredients that lead to departures from spherical symmetry. In particular, we focus on the effects of the tidal field associated with the hosting galaxy. We then take a stellar system on a circular orbit inside a galaxy represented as a "frozen" external field. The equilibrium distribution function is obtained from the one describing the spherical case by replacing the energy integral with the relevant Jacobi integral in the presence of the external tidal field. Then the construction of the model requires the investigation of a singular perturbation problem for an elliptic partial differential equation with a free boundary, for which we provide a method of solution to any desired order, with explicit solutions to two orders. We outline the relevant parameter space, thus opening the way to a systematic study of the properties of a two-parameter family of physically justified non-spherical models of quasi-relaxed stellar systems. The general method developed here can also be used to construct models for which the non-spherical shape is due to internal rotation. Eventually, the models will be a useful tool to investigate whether the shapes of globular clusters are primarily determined by internal rotation, by external tides, or by pressure anisotropy.

The globular clusters (GC) of our Galaxy have been found to lie close to a plane in the log(R_e), log(sigma), SB_e space, on the continuation of the Fundamental Plane (FP) known to characterize the properties of early-type galaxies. We reexamine the issue on a sample of 48 GCs selected in terms of homogeneity criteria for the photometric data available from the literature and perform a model-independent analysis of surface brightness profiles (SBP) and distance moduli, estimating error bars and studying selection effects with non-parametric statistical tests. We determine the coefficients of the FP and their error bars. The scatter from the FP relation is likely to be intrinsic, i.e. not due to measurement errors only. We find that in the standard FP coordinates our sample occupies a slim, axisymmetric region of parameter space, suggesting that the scaling relation might be around a Fundamental Line, rather than a plane, confirming a result noted earlier. This is likely to be the origin of the difficulties in the fit by a plane mentioned in previous investigations. Such FL relation would imply a pure photometric scaling law, which might be tested on wider samples and on extra-galactic GC systems. We find a correlation of the residuals from the FP relation with the central slope of the SBP. No other correlations are found. Finally, we reconstruct the distribution of the values of the quantity log(K_V/(M/L)) (virial coefficient divided by the mass-to-light ratio) through kernel density estimation and find evidence for bimodality, which suggests that the galactic GC system may be composed of at least two dynamically different populations. Yet, these populations do not reflect the standard dichotomy between disk and halo clusters. (abridged).

Gravitational lens modeling is presented for the first discovered example of a three-component source for which each component is quadruply imaged. The lens is a massive galaxy member of the cluster Cl J0152.7-1357 at z ~ 0.84. Taking advantage of this exceptional configuration and of the excellent angular resolution of the HST/ACS, we measure the properties of the lens. Several parametric macroscopic models were developed for the lens galaxy, starting from pointlike to extended image models. For a lens model in terms of a singular isothermal sphere with external shear, the Einstein radius is found to be R_{E} = (9.54 +/- 0.15) kpc. The external shear points to the cluster's northern mass peak. The unknown redshift of the source is determined to be higher than 1.9 and lower than 2.9. Our estimate of the lensing projected total mass inside the Einstein radius, M_{len}(R < 9.54 kpc), depends on the source distance and lies between 4.6 and 6.2 x 10^{11} M_{Sun}. This result turns out to be compatible with the dynamical estimate based on an isothermal model. By considering the constraint on the stellar mass-to-light ratio that comes from the evolution of the Fundamental Plane, we can exclude the possibility that at more than 4 sigma level the total mass enclosed inside the Einstein ring is only luminous matter. Moreover, the photometric-stellar mass measurement within the Einstein radius gives a minimum value of 50% (1 sigma) for the dark-to-total matter fraction. The lensing analysis has allowed us to investigate the distribution of mass of the deflector, also providing some interesting indications on scales that are larger and smaller than the Einstein radius of the lens galaxy. The combination of different diagnostics has proved to be essential in determining quantities that otherwise would have not been directly measurable with the current data.

[Abridged] We consider two clusters (A496 and Coma) that are representative of the two classes of cool-core and non-cool-core clusters. We first refer to a two-component dynamical model that ignores the contribution from the galaxy density distribution and study the condition of hydrostatic equilibrium for the hot intracluster medium (ICM) under the assumption of spherical symmetry, in the presence of dark matter. We model the ICM density distribution in terms of a standard $\beta$-model with $\beta=2/3$, i.e. with a distribution similar to that of a regular isothermal sphere (RIS), and fit the observed X-ray brightness profiles. With the explicit purpose of ignoring cosmological arguments, we na\"ively assume that dark matter, if present, has an analogous density distribution, with the freedom of two different density and length scales. The relative distribution of visible and dark matter is then derived by fitting the temperature data for the ICM under conditions of hydrostatic equilibrium. For both clusters, we find that dark matter is more concentrated with respect to visible matter. We then test whether the conclusion changes significantly when dark matter is taken to be distributed according to cosmologically favored density profiles and when the contribution of the mass contained in galaxies is taken into account. Although the qualitative conclusions remain unchanged, we find that the contribution of galaxies to the mass budget is more important than generally assumed. We also show that, without resorting to additional information on the small scale, it is not possible to tell whether a density cusp is present or absent in these systems. [Abridged]

We show how the combination of observations related to strong gravitational lensing and stellar dynamics in ellipticals offers a new way to measure the cosmological matter and dark-energy density parameters. A gravitational lensing estimate of the mass enclosed inside the Einstein circle can be obtained by measuring the Einstein angle, once the critical density of the system is known. A model-dependent dynamical estimate of this mass can also be obtained by measuring the central velocity dispersion of the stellar component. By assuming the well-tested homologous 1/r^{2} profile for the total density distribution in the lens elliptical galaxies, these two mass measurements can be properly compared. Thus, a relation between the Einstein angle and the central stellar velocity dispersion is derived, and the cosmological matter and the dark-energy density parameters can be estimated from this. We determined the accuracy of the cosmological parameter estimates by means of simulations that include realistic measurement uncertainties on the relevant quantities. Interestingly, the expected constraints on the cosmological parameter plane are complementary to those coming from other observational techniques. Then, we applied the method to the data sets of the Sloan Lens ACS and the Lenses Structure and Dynamics Surveys, and showed that the concordance value between 0.7 and 0.8 for the dark-energy density parameter is included in our 99% confidence regions. The small number of lenses available to date prevents us from precisely determining the cosmological parameters, but it still proves the feasibility of the method. When applied to samples made of hundreds of lenses that are expected to become available from forthcoming surveys, this technique will be an important tool for measuring the geometry of the Universe.

In this paper we study the infall, by dynamical friction, of rigid "satellites" taken in a variety of initial configurations and models. We thus measure how dynamical friction depends on the density concentration and on the pressure anisotropy of the host galaxy and, separately, the response of the galaxy during the process. The investigation is carried out by using a numerical laboratory set up in previous papers of this series. We find that the pressure anisotropy present in the host galaxy has little effect on dynamical friction. Instead, the shape of the galaxy density profile is very important. The classical idealized description appears to be applicable to galaxy models characterized by a small density gradient but not to concentrated models. In contrast to the behavior found in small density gradient models, the orbits of satellites captured in concentrated models are not circularized by dynamical friction. The induced evolution in the host galaxy reflects the initial conditions adopted for the captured satellite. Satellites spiraling in on quasi-circular orbits tend to modify the pressure tensor of the host galaxy in the tangential direction and the galaxy shape from spherical to oblate, while satellites captured along quasi-radial orbits tend to induce pressure anisotropy in the radial direction and to generate a prolate shape. This result suggests that the shape of early-type galaxies may result from occasional mergers rather than being directly related to the effectiveness of the radial-orbit instability during the process of formation via collisionless collapse, as often argued in the past.

Context. We consider a simple self-gravitating disk, made of two fluid components characterized by different effective thermal speeds and interacting with one another only through gravity; two-component models of this type have often been considered in order to estimate the impact of the cold interstellar medium on gravitational instabilities in star-dominated galaxy disks. Aims. This simple model allows us to produce a unified description of instabilities in non-viscous self-gravitating disks, some originating from Jeans collapse, and others from the relative motion between the two components. In particular, the model suggests that the small streaming velocity between the two components associated with the so-called asymmetric drift may be the origin of instability for suitable non-axisymmetric perturbations. Methods. The result is obtained by examining the properties of a local, linear dispersion relation for tightly wound density waves in such two-component model. The parameters characterizing the equilibrium model and the related dispersion relation allow us to recover as natural limits the cases, known in the literature, in which the relative drift between the two components is ignored. Results. Dynamically, the instability is similar to (although gentler than) that known to affect counter-rotating disks. However, in contrast to the instability induced by counter-rotation, which is a relatively rare phenomenon, the mechanism discussed in this paper is likely to be rather common in nature. Conclusions. We briefly indicate some consequences of the instability on the evolution of galaxy disks and possible applications to other astrophysical systems, in particular to protostellar disks and accretion disks.

We consider the Fundamental Plane of elliptical galaxies lensed by the gravitational field of a massive deflector (typically, a cluster of galaxies). We show that the Fundamental Plane relation provides a straightforward measurement of the projected mass distribution of the lens with a typical accuracy of ~0.15 in the dimensionless column density kappa. The proposed technique breaks the mass-sheet degeneracy completely and is thus expected to serve as an important complement to other lensing-based analyses. Moreover, its ability to measure directly the mass distribution on the small pencil beams that characterize the size of background galaxies may lead to crucial tests for current scenarios of structure formation.

Many astrophysical problems, ranging from structure formation in cosmology to dynamics of elliptical galaxies, refer to slow processes of evolution of essentially collisionless self-gravitating systems. In order to determine the relevant quasi-equilibrium configuration at time t from given initial conditions, it is often argued that such slow evolution may be approximated in terms of adiabatic evolution, for the calculation of which efficient semi--analytical techniques are available. Here we focus on the slow process of evolution, induced by dynamical friction of a host stellar system on a minority component of "satellites", that we have investigated in a previous paper, to determine to what extent an adiabatic description might be applied. The study is realized by comparing directly N--body simulations of the stellar system evolution (in two significantly different models) from initial to final conditions in a controlled numerical environment. We demonstrate that for the examined process the adiabatic description is going to provide incorrect answers, not only quantitatively, but also qualitatively. The two classes of models considered exhibit generally similar trends in evolution, with one exception noted in relation to the evolution of the total density profile. This simple conclusion should be taken as a warning against the indiscriminate use of adiabatic growth prescriptions in studies of structure of galaxies.

It is well known that the simple criterion proposed originally by Polyachenko and Shukhman (1981) for the onset of the radial orbit instability, although being generally a useful tool, faces significant exceptions both on the side of mildly anisotropic systems (with some that can be proved to be unstable) and on the side of strongly anisotropic models (with some that can be shown to be stable). In this paper we address two issues: Are there processes of collisionless collapse that can lead to equilibria of the exceptional type? What is the intrinsic structural property that is responsible for the sometimes noted exceptional stability behavior? To clarify these issues, we have performed a series of simulations of collisionless collapse that start from homogeneous, highly symmetrized, cold initial conditions and, because of such special conditions, are characterized by very little mixing. For these runs, the end-states can be associated with large values of the global pressure anisotropy parameter up to 2K_r/K_T \approx 2.75. The highly anisotropic equilibrium states thus constructed show no significant traces of radial anisotropy in their central region, with a very sharp transition to a radially anisotropic envelope occurring well inside the half-mass radius (around 0.2 r_M). To check whether the existence of such almost perfectly isotropic "nucleus" might be responsible for the apparent suppression of the radial orbit instability, we could not resort to equilibrium models with the above characteristics and with analytically available distribution function; instead, we studied and confirmed the stability of configurations with those characteristics by initializing N-body approximate equilibria (with given density and pressure anisotropy profiles) with the help of the Jeans equations.

We present 21-cm line observations of the spiral galaxy NGC 4559, made with the Westerbork Synthesis Radio Telescope. We have used them to study the HI distribution and kinematics, the relative amount and distribution of luminous and dark matter in this galaxy and, in particular, the presence of extra-planar gas. Our data do reveal the presence of such a component, in the form of a thick disk, with a mass of 5.9 x 10^8 Mo (one tenth of the total HI mass) and a mean rotation velocity 25-50 km/s lower than that of the thin disk. The extra-planar gas may be the result of galactic fountains but accretion from the IGM cannot be ruled out. With this study we confirm that lagging, thick HI layers are likely to be common in spiral galaxies.

We start from a study of the density-potential relation for classical homeoids in terms of an asymptotic expansion for small deviations from spherical symmetry. We then show that such expansion is a useful device that allows us to construct a variety of exact density-potential pairs with spheroidal, toroidal, or triaxial shapes for which the deviation from spherical symmetry is finite. As concrete analytical applications, we describe: (1) The construction of a family of toroidal axisymmetric density-potential pairs one of which is associated with a perfectly flat rotation curve (for a member of this family, the supporting two-integral phase-space distribution function is obtained in closed form); (2) The determination of the aperture velocity dispersion in a wide class of two-integral axisymmetric models not stratified on homeoids with central black hole, which may be useful for the discussion of the dynamical contributions to the characteristics of the Fundamental Plane of early-type galaxies; and (3) For such class of models, the construction of the v/sigma-ellipticity relation, often considered to assess the role of rotation in the structure of elliptical galaxies.

[Abridged] N-body simulations of collisionless collapse have offered important clues to the construction of realistic stellar dynamical models of elliptical galaxies. Understanding this idealized and relatively simple process, by which stellar systems can reach partially relaxed equilibrium configurations (characterized by isotropic central regions and radially anisotropic envelopes), is a prerequisite to more ambitious attempts at constructing physically justified models of elliptical galaxies in which the problem of galaxy formation is set in the generally accepted cosmological context of hierarchical clustering. In a previous paper, we have discussed the dynamical properties of a family of models of partially relaxed stellar systems (the f^{(\nu)} models), designed to incorporate the qualitative properties of the products of collisionless collapse at small and at large radii. Here we revisit the problem of incomplete violent relaxation, by making a direct comparison between the detailed properties of such family of models and those of the products of collisionless collapse found in N-body simulations that we have run for the purpose. Surprisingly, the models thus identified are able to match the simulated density distributions over nine orders of magnitude and also to provide an excellent fit to the anisotropy profiles and a good representation of the overall structure in phase space.

[Abridged] Recently we have found that a family of models of partially relaxed, anisotropic stellar systems, inspired earlier by studies of incomplete violent relaxation, exhibits some interesting thermodynamic properties. Here we present a systematic investigation of its dynamical characteristics, in order to establish the basis for a detailed comparison with simulations of collisionless collapse, planned for a separate paper. For a full comparison with the observations of elliptical galaxies, the models should be extended to allow for the presence a sizable dark halo and of significant rotation. In the spherical limit, the family is characterized by two dimensionless parameters, i.e. \Psi, measuring the depth of the galaxy potential, and \nu, defining the form of a third global quantity Q, which is argued to be approximately conserved during collisionless collapse (in addition to the total energy and the total number of stars).

N-body simulations of collisionless collapse have offered important clues to the construction of realistic stellar dynamical models of elliptical galaxies. Such simulations confirm and quantify the qualitative expectation that rapid collapse of a self-gravitating collisionless system, initially cool and significantly far from equilibrium, leads to incomplete relaxation. In this paper we revisit the problem, by comparing the detailed properties of a family of distribution functions derived from statistical mechanics arguments to those of the products of collisionless collapse found in N-body simulations.

We present preliminary results on a new family of distribution functions that are able to generate axisymmetric, truncated (i.e., finite size) stellar dynamical models characterized by toroidal shapes. The relevant distribution functions generalize those that are known to describe polytropic spheres, for which all the dynamical and structural properties of the system can be expressed in explicit form as elementary functions of the system gravitational potential. The model construction is then completed by a numerical study of the associated Poisson equation. We note that our axisymmetric models can also include the presence of an external gravitational field, such as that produced by a massive disk or by a central mass concentration (e.g., a supermassive black hole).

N-body simulations of collisionless collapse have offered important clues to the construction of realistic stellar dynamical models of elliptical galaxies. Such simulations confirm and quantify the qualitative expectation that rapid collapse of a self-gravitating collisionless system, initially cool and significantly far from equilibrium, leads to incomplete relaxation, that is to a quasi-equilibrium configuration characterized by isotropic, quasi-Maxwellian distribution of stellar orbits in the inner regions and by radially biased anisotropic pressure in the outer parts. In earlier studies, as illustrated in a number of papers several years ago (see Bertin et al. 1993 and references therein), the attention was largely focused on the successful comparison between the models (constructed under the qualitative clues offered by the N-body simulations mentioned above) and the observations. In this paper we revisit the problem of incomplete violent relaxation, by making a direct comparison between the detailed properties of a family of distribution functions and those of the products of collisionless collapse found in N-body simulations.

Studies of the Spectral Energy Distribution of Young Stellar Objects suggest that the outer disk of FU Orionis objects might be self-gravitating. In this paper we propose a method to test directly whether, in these objects, significant deviations from Keplerian rotation occur. In a first approach, we have used a simplified model of the disk vertical structure that allows us to quickly bring out effects related to the disk self-gravity. We find that the often studied optical and near-infrared line profiles are produced too close to the central object to provide significant evidence for non-Keplerian rotation. Based on parameters relevant for the case of FU Ori, we show that high-resolution long-wavelength spectroscopy, of the far-infrared H$_2$ pure rotational lines (sometimes observed in ``passive'' protostellar disks) and sub-mm CO lines, should be well suited to probe the rotation curve in the outer disk, thus measuring to what extent it is affected by the disk self-gravity. The results of the present exploratory paper should be extended soon to a more realistic treatment of the disk vertical structure.

The main goal of this paper is to set up a numerical laboratory for the study of the slow evolution of the density and of the pressure tensor profiles of an otherwise collisionless stellar system, as a result of the interactions with a minority component of heavier "particles".The effects that we would like to study are those attributed to slow collisional relaxation and generically called "dynamical friction". We start by revisiting the problem of the sinking of a satellite inside an initially isotropic, non-rotating, spherical galaxy, which we follow by means of N-body simulations using about one million particles. We then consider a quasi-spherical problem, in which the satellite is fragmented into a set of many smaller masses with a spherically symmetric initial density distribution. In a wide set of experiments, designed in order to bring out effects genuinely associated with dynamical friction, we are able to demonstrate the slow evolution of the density profile and the development of a tangentially biased pressure in the underlying stellar system, while we briefly address the issue of the circularization of orbits induced by dynamical friction on the population of fragments. The results of the simulations presented here and others planned for future investigations allow us to study the basic mechanisms of slow relaxation in stellar systems and thus may be of general interest for a variety of problems, especially in the cosmological context. Here our experiments are conceived with the specific goal to clarify some mechanisms that may play a role in the evolution of an elliptical galaxy as a result of the interaction between the stars and a significant population of globular clusters or of the merging of a large number of small satellites.

Studies of the Fundamental Plane of early-type galaxies, from small to intermediate redshifts, are often carried out under the guiding principle that the Fundamental Plane reflects the existence of an underlying mass-luminosity relation for such galaxies, in a scenario where elliptical galaxies are homologous systems in dynamical equilibrium. Here I will re-examine the issue of whether empirical evidence supports the view that significant systematic deviations from strict homology occur in the structure and dynamics of bright elliptical galaxies. In addition, I will discuss possible mechanisms of dynamical evolution for these systems, in the light of some classical thermodynamical arguments and of recent N-body simulations for stellar systems under the influence of weak collisionality.

The nucleus of the Seyfert galaxy NGC 1068 is believed to host a supermassive black hole. Evidence for the presence of a massive central object is provided by water maser emission, which displays a linear pattern in the sky, suggestive of a rotating disk. The rotating disk hypothesis is further strengthened by the declining shape of the derived rotation curve. Similar maser emission from NGC 4258 has led to a reliable estimate of the mass of the central black hole, because in this case the rotation curve is Keplerian. In the case of NGC 1068 the rotation curve traced by the water maser is non-Keplerian. In this paper we provide an interpretation of the non-Keplerian rotation in NGC 1068 by means of a self-gravitating accretion disk model. We obtain a good fit to the available data and derive a black hole mass M_{bh}=(8.0\pm 0.3) 10^6M_{sun}. The resulting disk mass is comparable to the black hole mass. As an interesting by-product of our fitting procedure, we are able to estimate the viscosity parameter, which turns out to be \alpha\approx 10^{-2}, in line with some theoretical expectations.

We examine the thermodynamical properties of a family of partially relaxed, anisotropic stellar systems, derived earlier from the Boltzmann entropy under the assumption that a third quantity Q is conserved in addition to the total energy and the total number of stars. We now show that the family of models conforms to the paradigm of the gravothermal catastrophe, which is expected to occur (in the presence of adequate energy transport mechanisms) when the one-parameter equilibrium sequence attains sufficiently high values of the concentration parameter; these are the values for which the models are well fitted by the R^(1/4) law. In the intermediate concentration regime the models belonging to the sequence exhibit significant deviations from the R^(1/4) law. Curiously, in the low-concentration regime, the global thermodynamical temperature associated with the models becomes negative when the models become too anisotropic so that they are unstable against the radial orbit instability; this latter behavior, while offering a new clue to the physical interpretation of the radial orbit instability, is at variance with respect to the low-concentration limit of the classical case of the isotropic, isothermal sphere investigated by Bonnor (1956) and Lynden-Bell & Wood (1968).

We start by studying a small set of objects characterized by photometric profiles that have been pointed out to deviate significantly from the standard R^{1/4} law. For these objects we confirm that a generic R^{1/n} law, with n a free parameter, can provide superior fits (the best-fit value of n can be lower than 2.5 or higher than 10), better than those that can be obtained by a pure R^{1/4} law, by an R^{1/4}+exponential model, and by other dynamically justified self--consistent models. Therefore, strictly speaking, elliptical galaxies should not be considered homologous dynamical systems. Still, a case for "weak homology", useful for the interpretation of the Fundamental Plane of elliptical galaxies, could be made if the best-fit parameter n, as often reported, correlates with galaxy luminosity L, provided the underlying dynamical structure also follows a systematic trend with luminosity. We demonstrate that this statement may be true even in the presence of significant scatter in the correlation n(L). Preliminary indications provided by a set of "data points" associated with a sample of 14 galaxies suggest that neither the strict homology nor the constant stellar mass--to--light solution are a satisfactory explanation of the observed Fundamental Plane (abridged).

We have measured the Fundamental Plane (FP) parameters for a sample of 30 field early-type galaxies (E/S0) in the redshift range 0.1<z<0.66. We find that: i) the FP is defined and tight out to the highest redshift bin; ii) the intercept \gamma evolves as d\gamma/dz=0.58+0.09-0.13 (for \Omega=0.3, \Omega_{\Lambda}=0.7), or, in terms of average effective mass to light ratio, as d\log(M/L_B)/dz=-0.72+0.11-0.16, i.e. faster than is observed for cluster E/S0 -0.49+-0.05. In addition, we detect [OII] emission >5\AA in 22% of an enlarged sample of 42 massive E/S0 in the range 0.1<z<0.73, in contrast with the quiescent population observed in clusters at similar z. We interpret these findings as evidence that a significant fraction of massive field E/S0 experiences secondary episodes of star-formation at z<1.

We have carried out a search for Lyman alpha emission from galaxies at z~=2.4 over a field of 1200 sq. arcmin using the CFH12K camera at the CFHT and a custom medium band filter. The search has uncovered 58 candidates, corresponding to a completeness-corrected source density of ~0.07 sq. arcmin^-2. Our sources have red colors (B-I~=1.8) which imply either that a large fraction of the light is highly reddened and we are detecting Lyman alpha through special lines of sight, or that these objects contain an underlying older stellar population. While for each individual object we cannot discriminate between these alternatives, we conclude that most of the objects actually contain an older component because the star formation rates inferred from the picture based on reddening, applied to all candidates, would imply an exceedingly high star formation rate, i.e. more than two orders of magnitude above the peak cosmic star formation rate (e.g. Lilly et al. 1996).

The long wavelength emission of protostellar objects is commonly attributed to a disk of gas and dust around the central protostar. In the first stages of disk accretion or in the case of high mass protostars, the disk mass is likely to be sufficiently large, so that the disk self-gravity may have an impact on the dynamics and the emission properties of the disk. In this paper we describe the spectral energy distribution (SED) produced by a simple, non-flaring, self-gravitating accretion disk model. Self-gravity is included in the calculation of the rotation curve of the disk and in the energy balance equation, as a term of effective heating related to Jeans instability. In order to quantify in detail the requirements on the mass of the disk and on the accretion rate posed on the models by realistic situations, we compare the SEDs produced by these models with the observed SEDs of a small sample of well-studied protostellar objects. We find that relatively modest disks - even lighter than the central star - can lead to an interesting fit to the infrared SED of the FU Orionis objects considered, while in the case of T Tauri stars the required parameters fall outside the range suggested as acceptable by the general theoretical and observational scenario. On the basis of the present results, we may conclude that the contribution of a self-gravitating disk is plausible in several cases (in particular, for FU Orionis objects) and that, in the standard irradiation dominated disk scenario, it would help softening the requirements encountered by Keplerian accretion models.

We have re-examined the classical problem of the macroscopic equation of state for a hydrostatic isothermal self-gravitating gas cloud bounded by an external medium at constant pressure. We have obtained analytical conditions for its equilibrium and stability without imposing any specific shape and symmetry to the cloud density distribution. The equilibrium condition can be stated in the form of an upper limit to the cloud mass; this is found to be inversely proportional to the power 3/2 of a form factor \mu characterizing the shape of the cloud. In this respect, the spherical solution, associated with the maximum value of the form factor, \mu = 1, turns out to correspond to the shape that is most difficult to realize. Surprisingly, the condition that defines the onset of the Bonnor instability (or gravothermal catastrophe) can be cast in the form of an upper limit to the density contrast within the cloud that is independent of the cloud shape. We have then carried out a similar analysis in the two-dimensional case of infinite cylinders, without assuming axisymmetry. The results obtained in this paper generalize well-known results available for spherical or axisymmetric cylindrical isothermal clouds that have had wide astrophysical applications, especially in the study of the interstellar medium.

We report on the study of a sample of 25 field early-type galaxies, in the redshift range z~0.1-0.5, selected on the basis of colours and morphology from the HST-MDS. Field early-type galaxies define a tight Fundamental Plane (FP) out to z~0.4, with scatter unchanged with respect to local samples, within the observational errors. The intermediate redshift FP is offset with respect to the FP of the Coma Cluster. The offset of the FP is found to increase with redshift. The evolution of the FP is studied quantitatively with a Bayesian-Montecarlo technique. By applying this technique, we find that the offset of the intercept of the FP (\Delta \gamma) with respect to the local FP increases as \Delta \gamma = \tau z with the following 68 per cent limits: 0.33<\tau<0.44 (for \Omega=1, \Omega_{\Lambda}=0) or 0.44<\tau<0.56 (for \Omega=0.3,\Omega_{\Lambda}=0.7). In addition, we interpret the results in terms of the evolution of the stellar populations, under the assumption of passive evolution. In a single-burst scenario, the observed properties are consistent with those of a stellar population formed at z>2 (for \Omega=1, \Omega_{\Lambda}=0, H_0=50 \kms Mpc^{-1}) or 0.8<z<1.6 (for \Omega=0.3, \Omega_{\Lambda}=0.7, H_0=65 \kms Mpc^{-1}). If a small fraction of the stellar mass is formed in a secondary burst, the primary burst may have occurred at higher z. Finally, the intercept and scatter of the FP found for field early-type galaxies and for cluster data at z~0.4 are mutually consistent, within the observational errors. If higher redshift (up to z=0.83) cluster data are considered, the ages of the stellar populations of field early-type galaxies inferred from a single-burst scenario are found to be marginally smaller than the ages derived for the cluster galaxies. [shortened]

A sample of field early-type galaxies (E/S0) at intermediate redshift ($z\sim0.1-0.6$) is selected, based on morphology and colours from HST-WFPC2 parallel images. Photometric structural parameters (effective radius $R_{\tx{e}}$ and effective surface brightness $SB_{\tx{e}}$) are derived through the F606W and F814W filters, using luminosity profile fitting and two-dimensional fitting techniques. The combined parameter that enters the Fundamental Plane ($\log R_{\tx{e}}-\beta SB_{\tx{e}}$, with $\beta\approx0.32$) is shown to suffer from significantly smaller uncertainties (r.m.s. 0.03) than the individual structural parameters (e.g. $\sim 15$ per cent r.m.s. on the effective radius). High signal-to-noise intermediate resolution spectra, taken at the ESO-3.6m, yield redshifts for 35 galaxies and central velocity dispersions for 22 galaxies. Central velocity dispersions are derived using a library of stellar templates covering a wide range of spectral types, in order to study the effects of templates mismatches. The average random error on the central velocity dispersion is found to be 8 per cent and the average systematic error due to template mismatch is found to be 5 per cent. The errors on velocity dispersion measurement and the effects of template mismatches are studied by means of extensive Montecarlo simulations. In addition, we investigate whether the determination of the velocity dispersion is sensitive to the spectral range used, finding that the value of velocity dispersion is unchanged when the spectral regions that include the absorption features Ca HK and NaD are masked out during the fit.

In the dynamics of accretion disks, the presence of collective effects associated with the self-gravity of the disk is expected to affect not only the momentum transport, but also the relevant energy balance equations, which could differ substantially from the non-self-gravitating case. Here we follow the picture that, when the disk is sufficiently cold, the stirring due to Jeans-related instabilities acts as a source of effective heating. The corresponding reformulation of the energy equations allows us to: ({\it i}) demonstrate how self-regulation can be established, so that the stability parameter $Q$ is maintained close to a threshold value, with weak dependence on radius; ({\it ii}) rediscuss the opacity properties in the self-gravitating regime. In particular, we show that, if cooling is dominated by {\it bremsstrahlung}, an optically thin stationary accretion solution is thermally stable, even in the non-advective case, provided the disk is self-gravitating. The details of the cooling function have little effect on the structure of such accretion disk, which is in any case brought by self-gravity to self-regulate. This condition of self-gravitating accretion is expected to be appropriate for the outer regions of many disks of astrophysical interest. With the reformulation of the energy equations described in this paper we have also secured: ({\it iii}) a starting point for the study of the emission properties of self-gravitating accretion disks; ({\it iv}) a tool to analyze the structure of the transition region, where the disk becomes self-gravitating.

We present preliminary results on the evolution of the stellar populations of field early-type galaxies (E/S0) from z=0.4 to z=0. The diagnostic tool used in this study is the Fundamental Plane (FP), a tight empirical correlation between their central velocity dispersion (sigma), effective radius (Re), and effective surface brightness (SBe), which is observed to hold in the local Universe. Using HST-WFPC2 archive images and spectra obtained at the ESO-3.6m telescope we measured the FP parameters for a sample of ~30 field E/S0s at z=0.2-0.4. Remarkably, field E/S0s at intermediate redshift also define a tight FP, with scatter unchanged with respect to that of local samples. The intermediate redshift FP is offset from the local one, in the sense that, for given RE and sigma, galaxies are brighter at z=0.4 than at z=0. The implication of the offset of the FP in terms of passive evolution of the stellar population depends on its star formation history. In a single burst scenario, the stellar populations of field E/S0s were formed at z=0.8-1.6 (Omega=0.3; Omega_{Lambda}=0.7; H_0=50 km/s/Mpc). Alternatively, the bulk of stars (90% in mass) can be formed at high redshift (z~3), and the rest in a secondary burst occurred more recently (z~0.5-0.8).

The analysis of the shear induced by a single cluster on the images of a large number of background galaxies is all centered around the curl-free character of a well-known vector field that can be derived from the data. Such basic property breaks down when the source galaxies happen to be observed through two clusters at different redshifts, partially aligned along the line of sight. In this paper we address the study of double lenses and obtain five main results. (i) First we generalize the procedure to extract the available information, contained in the observed shear field, from the case of a single lens to that of a double lens. (ii) Then we evaluate the possibility of detecting the signature of double lensing given the known properties of the distribution of clusters of galaxies. (iii) As a different astrophysical application, we demonstrate how the method can be used to detect the presence of a dark cluster that might happen to be partially aligned with a bright cluster studied in terms of statistical lensing. (iv) In addition, we show that the redshift distribution of the source galaxies, which in principle might also contribute to break the curl-free character of the shear field, actually produces systematic effects typically two orders of magnitude smaller than the double lensing effects we are focusing on. (v) Remarkably, a discussion of relevant contributions to the noise of the shear measurement has brought up an intrinsic limitation of weak lensing analyses, since one specific contribution, associated with the presence of a non-vanishing two-galaxy correlation function, turns out not to decrease with the density of source galaxies (and thus with the depth of the observations).

In this paper we describe some analytical properties of the R^{1/m} law proposed by Sersic (1968) to categorize the photometric profiles of elliptical galaxies. In particular, we present the full asymptotic expansion for the dimensionless scale factor b(m) that is introduced when referring the profile to the standard effective radius. Surprisingly, our asymptotic analysis turns out to be useful even for values of m as low as unity, thus providing a unified analytical tool for observational and theoretical investigations based on the R^{1/m} law for the entire range of interesting photometric profiles, from spiral to elliptical galaxies.

In this paper we describe some analytical properties of the R^{1/m} law proposed by Sersic to categorize the photometric profiles of elliptical galaxies. In particular, we present the full asymptotic expansion for the dimensionless scale factor b(m) that is introduced when referring the profile to the standard effective radius. Surprisingly, our asymptotic analysis turns out to be useful even for values of m as low as unity, thus providing a unified analytical tool for observational and theoretical investigations based on the R^{1/m} law for the entire range of interesting photometric profiles, from spiral to elliptical galaxies.

We consider a class of steady-state self-gravitating accretion disks for which efficient cooling mechanisms are assumed to operate so that the disk is self-regulated at a condition of approximate marginal Jeans stability. In an earlier paper, this scenario had been shown to lead naturally, in the absence of a central point mass, to a self-similar solution characterized by a flat rotation curve. In this article we investigate the entire parameter space available for such self-regulated accretion disks and provide two non-trivial extensions of the model. The first extension is that of a bimodal disk, obtained by partially relaxing the self-regulation constraint, so that full matching with an inner "standard" Keplerian accretion disk takes place. The second extension is the construction of self-regulated accretion disks embedded in a diffuse spherical "halo". The analysis is further strengthened by a careful discussion of the vertical structure of the disk, in such a way that the transition from self-gravity dominated to non-gravitating disks is covered uniformly.

Statistical analyses of observed galaxy distortions are often used to reconstruct the mass distribution of an intervening cluster responsible for gravitational lensing. In current projects, distortions of thousands of source galaxies have to be handled efficiently; much larger data bases and more massive investigations are envisaged for new major observational initiatives. In this article we present an efficient mass reconstruction procedure, a direct method that solves a variational principle noted in an earlier paper, which, for rectangular fields, turns out to reduce the relevant execution time by a factor from 100 to 1000 with respect to the fastest methods currently used, so that for grid numbers N = 400 the required CPU time on a good workstation can be kept within the order of 1 second. The acquired speed also opens the way to some long-term projects based on simulated observations (addressing statistical or cosmological questions) that would be, at present, practically not viable for intrinsically slow reconstruction methods.

We present measurements of the Fundamental Plane (FP) parameters (the effective radius, the mean effective surface brightness, and the central velocity dispersion) of six field elliptical galaxies at intermediate redshift. The imaging is taken from the Medium Deep Survey of the Hubble Space Telescope, while the kinematical data are obtained from long-slit spectroscopy using the 3.6-m ESO telescope. The Fundamental Plane appears well defined in the field even at redshift $\approx$ 0.3. The data show a shift in the FP zero point with respect to the local relation, possibly indicating modest evolution, consistent with the result found for intermediate redshift cluster samples. The FP slopes derived for our field data, plus other cluster ellipticals at intermediate redshift taken from the literature, differ from the local ones, but are still consistent with the interpretation of the FP as a result of homology, of the virial theorem and of the existence of a relation between luminosity and mass, $L \propto M^{\eta}$. We also derive the surface brightness vs. effective radius relation for nine galaxies with redshift up to $z \approx0.6$, and data from the literature; the evolution that can be inferred is consistent with what is found using the FP.

Recently, it has been shown that it is possible to reconstruct the projected mass distribution of a cluster from weak lensing provided that both the geometry of the universe and the probability distribution of galaxy redshifts are known; actually, when additional photometric data are taken to be available, the galaxy redshift distribution could be determined jointly with the cluster mass from the weak lensing analysis. In this paper we develop, in the spirit of a ``thought experiment,'' a method to constrain the geometry of the universe from weak lensing, provided that the redshifts of the source galaxies are measured. The quantitative limits and merits of the method are discussed analytically and with a set of simulations, in relation to point estimation, interval estimation, and test of hypotheses for homogeneous Friedmann-Lemaitre models. The constraints turn out to be significant when a few thousand source galaxies are used.

In this paper we provide a statistical analysis of the parameter-free method often used in weak lensing mass reconstructions. It is found that a proper assessment of the errors involved in such a non-local analysis requires the study of the relevant two-point correlation functions. After calculating the two-point correlation function for the reduced shear, we determine the expected error on the inferred mass distribution and on other related quantities, such as the total mass, and derive the error power spectrum. This allows us to optimize the reconstruction method, with respect to the kernel used in the inversion procedure. In particular, we find that curl-free kernels are bound to lead to more accurate mass reconstructions. Our analytical results clarify the arguments and the numerical simulations by Seitz & Schneider (1996).

Different options can be used in order to measure the shear from observations in the context of weak lensing. Here we introduce new methods where the isotropy assumption for the distribution of the source galaxies is implemented directly on the observed quadrupole moments. A quantitative analysis of the error associated with the finite number of source galaxies and with their ellipticity distribution is provided, applicable even when the shear is not weak. Monte Carlo simulations based on a realistic sample of source galaxies show that our procedure generally leads to errors ~30% smaller than those associated with the standard method of Kaiser and Squires (1993).

We consider a class of fully self-gravitating accretion disks, for which efficient cooling mechanisms are assumed to maintain the disk close to the margin of Jeans instability. For such self-regulated disks the equations become very simple in the outer regions, where the angular momentum convective transport approximately balances the viscous transport. These latter equations are shown to lead naturally to a self-similar solution with flat rotation curve, with circular velocity proportional to Mdot^(1/3) and essentially fixed opening angle.

A simple spherically-symmetric, steady-state, cooling-flow description with gas loss (following Sarazin \& Ashe 1989), within galaxy models constrained by radially extended stellar dynamical data, is shown to provide generally reasonable fits to the existing data on X-ray emission profiles and temperatures for a set of bright elliptical galaxies in Virgo and Fornax. Three free parameters are needed to specify the model: the external mass flux, the external pressure, and a dimensionless factor, which regulates the mass deposition rate along the flow. Three different assumptions on the supernova rate have been considered. A moderate value for the supernova rate in elliptical galaxies is found to be preferred. Confining pressures of $p_{ext}\sim4\div15\times10^3\kelvin$ cm$^{-3}$ and significant accretion rates of external material, up to $4\msolar/\yr$, are suggested by our models. A possible correlation between $L_X/L_B$ and the iron abundance in the gas inside ellipticals is pointed out.