Radiative transfer (RT) problems in which the source function includes a scattering-like integral are typical two-points boundary problems. Their solution via differential equations implies to make hypotheses on the solution itself, namely the specific intensity I(tau;n) of the radiation field. On the contrary, integral methods require to make hypotheses on the source function S(tau). It looks of course more reasonable to make hypotheses on the latter because one can expect that the run of S(tau) with depth be smoother than that of I(tau;n). In previous works we assumed a piece-wise parabolic approximation for the source function, which warrants the continuity of S(tau) and its first derivative at each depth point. Here we impose the continuity of the second derivative S"(tau). In other words, we adopt a cubic spline representation to the source function, which highly stabilize the numerical processes.

The structural parameters of a magnitude-limited sample of 148 unbarred S0-Sb galaxies were derived using the GASP2D algorithm and then analyzed to derive the intrinsic shape of their bulges. We developed a new method to derive the intrinsic shape of bulges based on photometric data and on the geometrical relationships between the apparent and intrinsic shapes of bulges and disks. The method is conceived as completely independent of the studied class of objects, and it can be applied whenever a triaxial ellipsoid embedded in an axisymmetric component is considered. We found that the intrinsic shape is well constrained for a subsample of 115 bulges with favorable viewing angles. A large fraction of them is characterized by an elliptical section (B/A<0.9). This fraction is 33%, 55%, and 43% if using their maximum, mean, or median equatorial ellipticity, respectively. Most of them are flattened along their polar axis. The distribution of triaxiality is strongly bimodal. This bimodality is driven by bulges with Sersic index n>2, or equivalently, by the bulges of galaxies with a bulge-to-total ratio B/T>0.3. Bulges with n<2 and with B/T<0.3 follow a similar distribution, which is different from that of bulges with n>2 and with B/T>0.3. In particular, bulges with n<2 and with B/T<0.3 show a larger fraction of oblate axisymmetric (or nearly axisymmetric) bulges, a smaller fraction of triaxial bulges, and fewer prolate axisymmetric (or nearly axisymmetric) bulges with respect to bulges with n>2 and with B/T>0.3, respectively. According to the numerical simulations of bulge formation, bulges with n<2, which show a high fraction of oblate axisymmetric (or nearly axisymmetric) shapes and have B/T<0.3, could be the result of dissipational minor mergers. Both major dissipational and dissipationless mergers seem to be required to explain the variety of shapes found for bulges with n>2 and B/T>0.3.

(Abridged) The structural parameters of a magnitude-limited sample of 148 unbarred S0-Sb galaxies were analyzed to derive the intrinsic shape of their bulges. We developed a new method to derive the intrinsic shape of bulges based on the geometrical relationships between the apparent and intrinsic shapes of bulges and disks. The equatorial ellipticity and intrinsic flattening of bulges were obtained from the length of the apparent major and minor semi-axes of the bulge, twist angle between the apparent major axis of the bulge and the galaxy line of nodes, and galaxy inclination. We found that the intrinsic shape is well constrained for a subsample of 115 bulges with favorable viewing angles. A large fraction of them is characterized by an elliptical section (B/A<0.9). This fraction is 33%, 55%, and 43% if using their maximum, mean, or median equatorial ellipticity, respectively. Most are flattened along their polar axis (C<(A+B)/2). The distribution of triaxiality is strongly bimodal. This bimodality is driven by bulges with Sersic index n>2, or equivalently, by the bulges of galaxies with a bulge-to-total ratio B/T>0.3. In particular, bulges with n\leq2 and with B/T\leq0.3 show a larger fraction of oblate axisymmetric (or nearly axisymmetric) bulges, a smaller fraction of triaxial bulges, and fewer prolate axisymmetric (or nearly axisymmetric) bulges with respect to bulges with n>2 and with B/T>0.3, respectively. According to predictions of the numerical simulations of bulge formation, bulges with n\leq2, which show a high fraction of oblate axisymmetric (or nearly axisymmetric) shapes and have B/T\leq0.3, could be the result of dissipational minor mergers. Both major dissipational and dissipationless mergers seem to be required to explain the variety of shapes found for bulges with n>2 and B/T>0.3.

The structural parameters of a magnitude-limited sample of 148 unbarred S0-Sb galaxies were derived to study the correlations between bulge and disk parameters as well as the probability distribution function (PDF) of the intrinsic equatorial ellipticity of bulges. A new algorithm (GASP2D) was used to perform the bidimensional bulge-disk decomposition of the J-band galaxy images extracted from the archive of the 2MASS survey. The PDF of intrinsic ellipticities was derived from the distribution of the observed ellipticities of the bulges and misalignments between the the bulges and disks. About 80% of the observed bulges are not oblate but triaxial ellipsoids. Their mean axial ratio in the equatorial plane is <B/A>=0.85. There is not significant dependence of their PDF on morphology, light concentration or luminosity. This has to be explained by the different scenarios of bulge formation.

(Abridged) A variety of formation scenarios was proposed to explain the diversity of properties observed in bulges. Studying their intrinsic shape can help in constraining the dominant mechanism at the epochs of their assembly. The structural parameters of a magnitude-limited sample of 148 unbarred S0--Sb galaxies were derived in order to study the correlations between bulges and disks as well as the probability distribution function (PDF) of the intrinsic equatorial ellipticity of bulges. It is presented a new fitting algorithm (GASP2D) to perform the two-dimensional photometric decomposition of galaxy surface-brightness distribution. This was assumed to be the sum of the contribution of a bulge and disk component characterized by elliptical and concentric isophotes with constant (but possibly different) ellipticity and position angles. Bulge and disk parameters of the sample galaxies were derived from the J-band images which were available in the Two Micron All Sky Survey. The PDF of the equatorial ellipticity of the bulges was derived from the distribution of the observed ellipticities of bulges and misalignments between bulges and disks. Strong correlations between the bulge and disk parameters were found. About 80% of bulges in unbarred lenticular and early-to-intermediate spiral galaxies are not oblate but triaxial ellipsoids. Their mean axial ratio in the equatorial plane is <B/A> = 0.85. There is not significant dependence of their PDF on morphology, light concentration, and luminosity. The interplay between bulge and disk parameters favors scenarios in which bulges assembled from mergers and/or grew over long times through disk secular evolution. But all these mechanisms have to be tested against the derived distribution of bulge intrinsic ellipticities.

A mayor problem that arises in the computation of stellar atmosphere models is the self consistent determination of the temperature distribution via the constraint of energy conservation. The energy balance includes the gains due to the absorption of radiation and the losses due to emission. It is well known that within each one of the two above integrals the part corresponding to spectral ranges whose opacity X(nu) is huge can overcome by many orders of magnitude the part that corresponds to the remaining frequencies. On the other hand, at those frequencies where X(nu) is very large, the mean intensity J(nu) of the radiation field shall be equal, up to many significant digits, to the source function S(nu) and consequently to the Planck function B(nu,T). Then their net share to the energy balance shall be null, albeit separately their contributions to the gain and loss integrals are the most important numerically. Thus the spectral range whose physical contribution to the overall balance is null will dominate numerically both sides of the energy balance equation, and consequently the errors on the determination of J(nu) and S(nu) at these frequencies will falsify the balance. It is possible to circumvent the numerical problem brought about by the foregoing circumstances by solving the radiative transfer equation for the variable I(n,nu) - S(nu), instead of the customary intensity I(n,nu). We present here a novel iterative algorithm, based on iteration factors already employed by us with success, which makes it possible a fast correction of the temperature by computing directly the difference between the radiative losses and gains at each step of the iterations.

We study the density profiles of collapsed galaxy-size dark matter halos with masses 1e11-5e12 Msun focusing mostly on the halo outer regions from the formal virial radius Rvir up to 5-7Rvir. We find that isolated halos in this mass range extend well beyond Rvir exhibiting all properties of virialized objects up to 2-3Rvir: relatively smooth density profiles and no systematic infall velocities. The dark matter halos in this mass range do not grow as one naively may expect through a steady accretion of satellites, i.e., on average there is no mass infall. This is strikingly different from more massive halos, which have large infall velocities outside of the virial radius. We provide accurate fit for the density profile of these galaxy-size halos. For a wide range (0.01-2)Rvir of radii the halo density profiles are fit with the approximation rho=rho_s exp(-2n[x^{1/n}-1])+rho_m, where x=r/r_s, rho_m is the mean matter density of the Universe, and the index n is in the range n=6-7.5. These profiles do not show a sudden change of behavior beyond the virial radius. For larger radii we combine the statistics of the initial fluctuations with the spherical collapse model to obtain predictions for the mean and most probable density profiles for halos of several masses. The model give excellent results beyond 2-3 formal virial radii.

If dark matter (DM), which is considered to constitute most of the mass of galaxies, is made of supersymmetric (SUSY) particles, the centers of galaxies should emit gamma-rays produced by their self-annihilation. We present accurate estimates of continuum gamma-ray fluxes due to neutralino annihilation in the central regions of the Milky Way. We use detailed models of our Galaxy, which satisfy available observational data, and include some important physical processes, which were previously neglected. Our models predict that spatially extended annihilation signal should be detected at high confidence levels by incoming experiments assuming that neutralinos make up most of the DM in the Universe and that they annihilate according to current SUSY models.

A method based on Lucy (1974, AJ 79, 745) iterative algorithm is developed to invert the equation of stellar statistics for the Galactic bulge and is then applied to the K-band star counts from the Two-Micron Galactic Survey in a number of off-plane regions (10 deg.>|b|>2 deg., |l|<15 deg.). The top end of the K-band luminosity function is derived and the morphology of the stellar density function is fitted to triaxial ellipsoids, assuming a non-variable luminosity function within the bulge. The results, which have already been outlined by Lopez-Corredoira et al.(1997, MNRAS 292, L15), are shown in this paper with a full explanation of the steps of the inversion: the luminosity function shows a sharp decrease brighter than M_K=-8.0 mag when compared with the disc population; the bulge fits triaxial ellipsoids with the major axis in the Galactic plane at an angle with the line of sight to the Galactic centre of 12 deg. in the first quadrant; the axial ratios are 1:0.54:0.33, and the distance of the Sun from the centre of the triaxial ellipsoid is 7860 pc. The major-minor axial ratio of the ellipsoids is found not to be constant. However, the interpretation of this is controversial. An eccentricity of the true density-ellipsoid gradient and a population gradient are two possible explanations. The best fit for the stellar density, for 1300 pc<t<3000 pc, are calculated for both cases, assuming an ellipsoidal distribution with constant axial ratios, and when K_z is allowed to vary. From these, the total number of bulge stars is ~ 3 10^{10} or ~ 4 10^{10}, respectively.

By considering that the distribution of measured light along any galactocentric radius of an elliptical galaxy can be well represented by a R^{1/n} Sersic law, we propose in this paper a ``discrete ordenate'' method, which, for any value of n, allows an explicit expression for the luminous density, rho(r), to be found that can be numerically evaluated to any required precision. Once we have this semi-analytical expression for the spatial density, rho(r), the mass distribution, M(r), the potential, phi(r), and the velocity dispersions, sigma_s^2(r), in the space and in the observational plane, sigma_p^2(R), can be computed in a straightforward manner.