The current attempt is aimed to honor the first centennial of Johannes Diderik van der Waals (VDW) awarding Nobel Prize in Physics. The VDW theory of ordinary fluids is reviewed in the first part of the paper, where special effort is devoted to the equation of state and the law of corresponding states. In addition, a few mathematical features involving properties of cubic equations are discussed, for appreciating the intrinsic beauty of the VDW theory. A theory of astrophysical fluids is shortly reviewed in the second part of the paper, grounding on the tensor virial theorem for two-component systems, and an equation of state is formulated with a convenient choice of reduced variables. Additional effort is devoted to particular choices of density profiles, namely a simple guidance case and two cases of astrophysical interest. The related macroisothermal curves are found to be qualitatively similar to VDW isothermal curves below the critical threshold and, for sufficiently steep density profiles, a critical macroisothermal curve exists, with a single horizontal inflexion point. Under the working hypothesis of a phase transition (assumed to be gas-stars) for astrophysical fluids, similar to the vapour-liquid phase transition in ordinary fluids, the location of gas clouds, stellar systems, galaxies, cluster of galaxies, on the plane scanned by reduced variables, is tentatively assigned. A brief discussion shows how van der Waals' two great discoveries, namely a gas equation of state where tidal interactions between molecules are taken into account, and the law of corresponding states, related to microcosmos, find a counterpart with regard to macrocosmos. In conclusion, after more than a century since the awarding of the Nobel Prize in Physics, van der Waals' ideas are still valid and helpful to day for a full understanding of the universe.

The occurrence of a desert between neutron star (NS) and black hole (BH) remnants is reviewed. The dependence of stellar remnants on the zero age main sequence (ZAMS) progenitor mass for solar metallicity is taken from a recent investigation and further effort is devoted to NS and BH remnants. In particular, a density parameter is defined and related to NS mass and radius. Spinning BHs in Kerr metrics are considered as infinitely thin, homogeneous, rigidly rotating disks in Newtonian mechanics. Physical parameters for nonrotating (TOV) and equatorial breakup (EQB) configurations are taken or inferred from a recent investigation with regard to 4 NS and 3 quark star (QS) physically motivated equation of state (EOS) kinds. A comparison is performed with counterparts related to nonrotating and maximally rotating BHs. The results are also considered in the light of empirical relations present in literature. With regard to $J$-$M$ relation, EQB configurations are placed on a sequence of similar slope in comparison to maximally rotating BHs, but shifted downwards due to lower angular momentum by a factor of about 3.6. Under the assumption heavy baryons are NS constituents and instantaneously undergo quark-level reactions, the energy released (or adsorbed) is calculated using results from a recent investigation. Even if NSs exclusively host heavy baryons of the kind considered, the total amount cannot exceed about 10\% of the binding energy, which inhibits supernova explosions as in supramassive white dwarf (WD) remnants or implosions into BH. Alternative channels for submassive ($2<=M/m_\odot<=4$) BH formation are shortly discussed.

A classical model of lunar recession is reviewed and updated, where input parameters are taken or inferred from earlier astronomical computation of insolation quantities on Earth back to 0.25 Gyr. Free parameters are (i) the transition age, where mean lunar recession is suddenly lowered, and (ii) the merging age, where Earth-Moon distance (EMD), drops to zero. Predicted mean EMD and length of day (LOD) slightly overstimate their counterparts from the above mentioned astronomical computation. Predicted LOD, where the effect of atmospheric tides and nontidal processes is considered and supposed to be time independent, slightly understimates a linear interpolation from paleontological data related to Phanerozoic, inferred in earlier investigations. If short-period ($\Delta t\approx10^{-3}$ Gyr) EMD fluctuations ($\Delta a\approx\mp0.5R_\oplus$), resulting from the above mentioned astronomical computation back to 0.25 Gyr, occur along the whole evolution of Earth-Moon system (EMS), then the predicted mean EMD is consistent with values inferred from three well studied paleontological data sets: Elatina-Reynella (ER), Big Cottonwood (BC), Weeli Wolli (WW), for a suitable choice of input parameters. The same holds for LOD where, in addition, computed values are consistent with a linear interpolation from paleontological data related to Proterozoic, inferred in earlier investigations. The effect of atmospheric tides and nontidal processes is also discussed. In conclusion, the current model can be considered as a useful zeroth-order approximation to the evaluation of lunar recession and LOD. More accurate paleontological data, expecially in connection with error evaluation, would be desirable in view of improved results and further constraints involving both Astronomy and Paleontology.

Concerning bivariate least squares linear regression, the classical results obtained for extreme structural models in earlier attempts are reviewed using a new formalism in terms of deviation (matrix) traces which, for homoscedastic data, reduce to usual quantities leaving aside an unessential (but dimensional) multiplicative factor. Within the framework of classical error models, the dependent variable relates to the independent variable according to the usual additive model. The classes of linear models considered are regression lines in the limit of uncorrelated errors in X and in Y. For homoscedastic data, the results are taken from earlier attempts and rewritten using a more compact notation. For heteroscedastic data, the results are inferred from a procedure related to functional models. An example of astronomical application is considered, concerning the [O/H]-[Fe/H] empirical relations deduced from five samples related to different stars and/or different methods of oxygen abundance determination. For low-dispersion samples and assigned methods, different regression models yield results which are in agreement within the errors for both heteroscedastic and homoscedastic data, while the contrary holds for large-dispersion samples. In any case, samples related to different methods produce discrepant results, due to the presence of (still undetected) systematic errors, which implies no definitive statement can be made at present. Asymptotic expressions approximate regression line slope and intercept variance estimators, for normal residuals, to a better extent with respect to earlier attempts. Related fractional discrepancies are not exceeding a few percent for low-dispersion data, which grows up to about 10% for large-dispersion data. An extension of the formalism to generic structural models is left to a forthcoming paper.

With regard to generic two-component systems, the theory of first variations of global quantities is reviewed and explicit expressions are inferred for subsystem potential energies and potential-energy tensors. Performing a conceptual experiment, a physical interpretation of subsystem potential energies and potential-energy tensors is discussed. Subsystem tidal radii are defined by requiring an unbound component in absence of the other one. To this respect, a few guidance examples are presented as: (i) an embedding and an embedded homogeneous sphere; (ii) an embedding and an embedded truncated, singular isothermal sphere where related centres are sufficiently distant; (iii) a homogeneous sphere and a Roche system i.e. a mass point surrounded by a vanishing atmosphere. The results are discussed and compared with the findings of earlier investigations.

The concept of reduced variables is revisited with regard to van der Waals' theory and an application is made to polytropic spheres, where the reduced radial coordinate is ${\rm red}(r)=r/R=\xi/\Xi$, $R$ radius, and the reduced density is ${\rm red}(\rho)=\rho/\lambda=\theta^n$, $\lambda$ central density. Reduced density profiles are plotted for several polytropic indexes within the range, $0\le n\le5$, disclosing two noticeable features. First, any point of coordinates, $({\rm red}(r),{\rm red}(\rho))$, $0\le{\rm red}(r)\le1$, $0\le{\rm red}(\rho)\le1$, belongs to a reduced density profile of the kind considered. Second, sufficiently steep i.e. large $n$ reduced density profiles exhibit an oblique inflection point, where the threshold is found to be located at $n=n_{\rm th}=0.888715$. Reduced pressure profiles, ${\rm red}(P)=P/\varpi=\theta^{n+1}$, $\varpi$ central pressure, Lane-Emden fucntions, $\theta=(\rho/\lambda)^{1/n}$, and polytropic curves, ${\rm red}(P)={\rm red}(P)({\rm red}(\rho))$, are also plotted. The method can be extended to nonspherical polytropes with regard to a selected direction, ${\rm red}(r)(\mu)=r(\mu)/R(\mu)=\xi(\mu)/\Xi(\mu)$. The results can be extended to polytropic spheres made of collisionless particles, for polytropic index within a more restricted range, $1/2\le n\le5$.

With respect to earlier investigations, the theory of multi-component, concentric, copolar, axisymmetric, rigidly rotating polytropes is improved and extended, including subsystems with nonzero density on the boundary and subsystems with intersecting boundaries. The formulation is restricted to two subsystems for simplicity but, in principle, can be extended to $N$ subsystems. Equilibrium configurations are independent of the nature of the fluid i.e. collisional or collisionless, provided the polytropic index lies within the range, $1/2\le n\le5$, as in one-component systems. The solution of the equilibrium equations is expanded in power series, which can be continued up to the boundary and outside via starting points placed at increasingly larger distance from the centre of mass. A detailed analysis is devoted to special cases where the solution of the equilibrium equations can be expressed analytically. Finally a guidance example is shown, involving homogeneous subsystems with intersecting boundaries, where a substantially flattened component extends outside a slightly flattened one.

Within the context of Newton's theory of gravitation, restricted to point-like test particles and central bodies, stable circular orbits in ordinary space are related to stable circular paths on a massless, unmovable, undeformable vortex-like surface, under the action of a tidal gravitational field along the symmetry axis. An interpretation is made in the light of a holographic principle, in the sense that motions in ordinary space are connected with motions on a selected surface and vice versa. Then ordinary space is conceived as a 3-hypersurface bounding a $n$-hypervolume where gravitation takes origin, within a $n$-hyperspace. The extension of the holographic principle to extra dimensions implies the existence of a minimum distance where test particles may still be considered as distinct from the central body. Below that threshold, it is inferred test particles lose theirs individuality and "glue" to the central body via unification of the four known interactions and, in addition, (i) particles can no longer be conceived as point-like but e.g., strings or membranes, and (ii) quantum effects are dominant and matter turns back to a pre-big bang state. A more detailed formulation including noncircular motions within the context of general relativity, together with further knowledge on neutron stars, quark stars and black holes, would provide further insight on the formulation of quantum gravity.

Pure power-law density profiles, $\rho(r)\propto r^{b-3}$, are classified in connection with the following reference cases: (i) isodensity, $b=3$, $\rho=$ const; (ii) isogravity, $b=2$, $g=$ const; (iii) isothermal, $b=1$, $v=[GM(r)/r]^{1/2}=$ const; (iv) isomass, $b=0$, $M=$ const. A restricted number of different families of density profiles including, in addition, cored power-law, generalized power-law, polytropes, are studied in detail with regard to both one-component and two-component systems. Considerable effort is devoted to the existence of an extremum point (maximum absolute value) in the gravitational acceleration within the matter distribution. Predicted velocity curves are compared to the data inferred from observations. Tidal effects on an inner subsystem are investigated and an application is made to globular clusters within the Galaxy. To this aim, the tidal radius is defined by balancing the opposite gravitational forces from the Galaxy and the selected cluster on a special point of the cluster boundary, lying between related centres of mass. The position of 17 globular clusters with respect to the stability region, where the tidal radius exceeds the observed radius, is shown for assigned dark-to-visible mass ratios and density profiles, among those considered, which are currently used for the description of galaxies and/or dark matter haloes.

A linear [Fe/H]-[O/H] relation is found for different stellar populations in the Galaxy (halo, thick disk, thin disk) from a data sample obtained in a recent investigation (Ram{\'\i}rez et al. 2013). These correlations support previous results inferred from poorer samples: stars display a "main sequence" expressed as [Fe/H] = $a$[O/H$]+b\mp\Delta b$ where a unit slope, $a=1$, implies a constant [O/Fe] abundance ratio. Oxygen and iron empirical abundance distributions are then determined for different subsamples, which are well explained by the theoretical predictions of multistage closed-(box+reservoir) (MCBR) chemical evolution models by taking into account the found correlations. The interpretation of these distributions in the framework of MCBR models gives us clues about inflow/outflow rates in these different Galactic regions and their corresponding evolution. Outflow rate for the thick and the thin disks are lower than the halo outflow rate. Moreover if the thin disk built up from the thick disk, both systems result of comparable masses. Besides that, the iron-to-oxygen yield ratio and the primary to not primary contribution ratio for the iron production are obtained from the data, resulting consistent with SNII progenitor nucleosynthesis and with the iron production from SNIa supernova events.

The computed $\log(\tau/{\rm y})$-$\log(m/m_\odot)$ relation for the stellar initial mass range, 0.6-120.0, and the stellar initial metallicity range, 0.0004-0.0500, tabulated in an earlier attempt (Portinari et al. 1998) is fitted to a good extent by a four-parameter curve, expressed by a double power-law, for assigned stellar initial metallicity, which can be reduced to a three-parameter curve, expressed by a single power-law, for the whole set of stellar initial metallicities. The relative errors do not exceed about 2% and 4%, respectively. The extent to which the interpolation curve, expressed by a single power-law, can be extrapolated towards both high-mass and low-mass stars, is also investigated. High-mass star lifetimes are understimated by a factor less than 2 up to $m/m_\odot= 1000$ and by a fiducial factor less than 4 up to infinite. Low-mass star lifetimes are overstimated by a factor of about 3 down to $m/m_\odot=0.25$ and by an unacceptably large factor down to $m/m_\odot=0.08$. The interpolation curve, expressed by a single power-law, is related (in differential form) to a generalization of the equation of the exponential decay, which could be the starting point for a theoretical interpretation. As a simple application, the star mass fraction of a single star generation with stellar initial mass function defined by a power-law, is plotted vs. the logarithmic stellar lifetime. The star mass fraction declines in time at a decreasing rate for mild stellar initial mass function and at an increasing rate for steep stellar initial mass function, where a linear trend is exhibited for a value of the exponent close to the Salpeter's value, equal to -2.35.

Axisymmetric, rigidly rotating polytropes are considered in the framework of both the original Chandrasekhar (C33) approximation and a different version (extended C33 approximation). Special effort is devoted to two specific points, namely (i) a contradiction between the binomial series evaluation and the vanishing density on the boundary, which affects the self-consistency of the above mentioned approximations, and (ii) the continuity of selected parameters as a function of the polytropic index. Concerning (i), it is shown the Emdem-Chandrasekhar (EC) associated functions are defined at any internal point even if related EC associated equations hold only for a particular subvolume, in the framework of both the C33 and the extended C33 approximation. Concerning (ii), the continuity may safely be established, on the boundary of the domain, for part of the parameters, while additional data are needed for the remaining part. Simple fitting curves, valid to a good extent for a wide range of polytropic index, involve exponential functions and, in a single case, two straight lines joined by a parabolic segment. The expression of physical parameters in terms of the polytropic index can be used in building up sequences of configurations with changing density profile for assigned mass and angular momentum.

Linear [Q/H]-[O/H] relations, Q = Na, Mg, Si, Ca, Ti, Cr, Fe, Ni, are inferred from a sample (N=67) of recently studied FGK-type dwarf stars in the solar neighbourhood including different populations. Regression line slope and intercept estimators and related variance estimators are determined. With regard to the straight line, [Q/H] = aQ [O/H] + bQ, sample stars display along a "main sequence", [Q,O] = [aQ,bQ, Delta bQ], leaving aside the two OL stars which, in most cases (e.g., Na), lie outside. A unit slope, aQ = 1, implies Q is a primary element synthesised via SNII progenitors in presence of universal stellar initial mass function (defined as simple primary element). To this respect, Mg, Si, Ti, show aQ = 1 within 2 sigma(aQ); Cr, Fe, Ni, within 3 sigma(aQ); Na, Ca, exceeding 3 sigma(aQ). The empirical, differential element abundance distributions are inferred from different subsamples, where related regression lines represent their theoretical counterparts within the framework of simple MCBR (multistage closed box + reservoir) chemical evolution models. Hence the fractional yields are determined and (as an example) a comparison is shown with their theoretical counterparts inferred from SNII progenitor nucleosynthesis under the assumption of a power-law stellar initial mass function. The generalized fractional yields are determined regardless from the chemical evolution model. The ratio of outflow to star formation rate is compared for different populations, in the framework of simple MCBR models. The opposite situation of element abundance variation entirely due to cosmic scatter is also considered under reasonable assumptions. The related differential element abundance distribution fits to the data as well as its counterpart inferred in the opposite limit of instantaneous mixing in presence of chemical evolution.

Starting from recent determination of Fe, O, Na abundance on a restricted sample (N=67) of halo and thick disk stars, a natural and well motivated selection criterion is defined for the classification of globular cluster stars. An application is performed to M13 using a sample (N=112) for which Fe, O, Na abundance has been recently inferred from observations. A comparison is made between current and earlier M13 star classification. Both O and Na empirical differential abundance distributions are determined for each class and the whole sample (with the addition of Fe in the last case) and compared with their theoretical counterparts due to cosmic scatter obeying a Gaussian distribution whose parameters are inferred from related subsamples. The procedure is repeated with regard to four additional classes according if oxygen and sodium abundance is above (stage CE) or below (stage AF) a selected threshold. Both O and Na empirical differential abundance distributions, related to the whole sample, exhibit a linear fit for AF and CE stage. Within the errors, the oxygen slope for CE stage is equal and opposite in sign with respect of the sodium slope for AF stage, while the contrary holds in dealing with the oxygen slope for AF stage with respect to the sodium slope for CE stage. In the light of simple models of chemical evolution applied to M13, the O-Na anticorrelation and the primordial to current M13 mass ratio are discussed. The proposed classification of globular cluster stars may be extended in a twofold manner, namely to (i) elements other than Na and Fe and (ii) globular clusters other than M13.

Simple MCBR models of chemical evolution are extended to the limit of dominant gas inflow or outflow with respect to gas locked up into long-lived stars and remnants. For an assigned empirical differential oxygen abundance distribution, which can be linearly fitted, a family of theoretical curves is built up with assigned prescriptions. For curves with increasing cut parameter, the gas mass fraction locked up into long-lived stars and remnants is found to attain a maximum and then decrease towards zero as the flow tends to infinity, while the remaining parameters show a monotonic trend. The theoretical integral oxygen abundance distribution is also expressed. An application is performed to the empirical distribution deduced from two different samples of disk stars, for both the thin and the thick disk. The constraints on formation and evolution are discussed in the light of the model. The evolution is tentatively subdivided into four stages, A, F, C, E. The empirical distribution related to any stage is fitted by all curves for a wide range of the cut parameter. The F stage may safely be described by a steady inflow regime, implying a flat theoretical distribution, in agreement with the results of hydrodynamical simulations. Finally, (1) the change of fractional mass due to the extension of the linear fit to the empirical distribution, towards both the (undetected) low-metallicity and high-metallicity tail, is evaluated and (2) the idea of a thick disk-thin disk collapse is discussed, in the light of the model.

The observed relationship between stellar mass and effective radius for early type galaxies, pointed out by many authors, is interpreted in the context of Clausius' virial maximum theory. In this view, it is strongly underlined that the key of the above mentioned correlation is owing to the presence of a deep link between cosmology and the existence of the galaxy Fundamental Plane. Then the ultimate meaning is: understanding visible mass - size correlation and/or Fundamental Plane means understanding how galaxies form. The mass - size relationship involves baryon (mainly stellar) mass and its typical dimension related to the light, but it gets memory of the cosmological mass variance at the equivalence epoch. The reason is that the baryonic component virializes by sharing virial energy in about equal amount between baryons and dark matter, this sharing depending, in turn, on the steepness of the dark matter distribution. The general strategy consists in using the two-component tensor virial theorem for determining the virialized baryonic configurations. A King and a Zhao density profile are assumed for the inner baryonic and the outer dark matter component, respectively, at the end of the relaxation phase. All the considerations are restricted to spherical symmetry for simplicity. The effect of changing the dark-to-baryon mass ratio, m, is investigated inside a LambdaCDM scenario. A theoretical mass - size relation is expressed for the baryonic component, which fits fairly well to the data from a recently studied galaxy sample. Finally, the play of intrinsic dispersion on the mass ratio, m, is discussed in the light of the cusp/core problem and some consequences are speculated about the existence of a limit, m_l, expected by the theory.

Concerning bivariate least squares linear regression, the classical approach pursued for functional models in earlier attempts is reviewed using a new formalism in terms of deviation (matrix) traces. Within the framework of classical error models, the dependent variable relates to the independent variable according to the usual additive model. Linear models of regression lines are considered in the general case of correlated errors in X and in Y for heteroscedastic data. The special case of (C) generalized orthogonal regression is considered in detail together with well known subcases. In the limit of homoscedastic data, the results determined for functional models are compared with their counterparts related to extreme structural models. While regression line slope and intercept estimators for functional and structural models necessarily coincide, the contrary holds for related variance estimators even if the residuals obey a Gaussian distribution, with a single exception. An example of astronomical application is considered, concerning the [O/H]-[Fe/H] empirical relations deduced from five samples related to different stars and/or different methods of oxygen abundance determination. For selected samples and assigned methods, different regression models yield consistent results within the errors for both heteroscedastic and homoscedastic data. Conversely, samples related to different methods produce discrepant results, due to the presence of (still undetected) systematic errors, which implies no definitive statement can be made at present. A comparison is also made between different expressions of regression line slope and intercept variance estimators, where fractional discrepancies are found to be not exceeding a few percent, which grows up to about 20% in presence of large dispersion data.

Simple closed-box (CB) models of chemical evolution are extended on two respects: (i) simple closed-(box+reservoir) (CBR) models allowing gas outflow from the box into the reservoir or gas inflow into the box from the reservoir with same composition as the preexisting gas and rate proportional to the star formation rate, and (ii) simple multistage closed-(box+reservoir) (MCBR) models allowing different stages of evolution characterized by different inflow or outflow rates. The stellar initial mass function is assumed to be universal, and mass conservation holds for the whole system (box+reservoir) while it is violated for each subsystem (box and reservoir). The theoretical differential oxygen abundance distribution (TDOD) predicted by the model, under the assumption of instantaneous recycling, is a continuous broken line, where different slopes are related to different inflow rates. For an application of the model (a) a fictitious sample is built up from two distinct samples and taken as representative of the inner Galactic halo, and (b) different [O/H]-[Fe/H] empirical relations are deduced from five different samples related to different methods, and two of them are selected for determining the empirical differential oxygen abundance distribution (EDOD) with regard to the fictitious sample. In both cases the EDOD is represented, to an acceptable extent, as a continuous broken line. If the inner halo and the metal-poor bulge (after the inflow stage) are represented by the box and the reservoir, respectively, then the inner halo fractional mass (normalized to the halo) is comparable with, or exceeds by a factor up to 4, the metal-poor bulge fractional mass (normalized to the bulge), for current estimates of the halo-to-bulge mass ratio of about 0.05-0.10.

Macrogases are defined as two-component, large-scale celestial objects where the subsystems interact only via gravitation. The macrogas equation of state is formulated and compared to the van der Waals equation of state for ordinary gases. By analogy, it is assumed that real macroisothermal curves in macrogases occur as real isothermal curves in ordinary gases, where a phase transition takes place along a horisontal line in the macrovolume-macropressure (Mv-Mp) plane. A simple guidance case and two density profiles which satisfactorily fit to observations or simulations, are studied in detail. For sufficiently steep density profiles, a critical macroisothermal curve exists as shown by ordinary gases, where the critical point coincides with the horisontal inflexion point. By analogy with ordinary gases, the first quadrant of the (Mv-Mp) plane may be divided into three parts, namely (i) the G region, where only gas exists; (ii) the S region, where only stars exist; (iii) the GS region, where both gas and stars exist. An application is made to a subsample of elliptical galaxies investigated within the SAURON project. Different models characterized by equal subsystem mass ratio and different scaled truncation radii, are considered and the related position of sample objects on the (Mv-Mp) plane is determined. Tipically, fast rotators are found to lie within the S region, while slow rotators are close (from both sides) to the boundary between the S and the GS region. The net effect of the uncertainty affecting observed quantities, on the position of sample objects on the (Mv-Mp) plane, is also investigated. Finally, a principle of corresponding states is formulated for macrogases with assigned density profiles and scaled truncation radii.

Tidal interactions between neighboring objects span across the whole admissible range of lengths in nature: from, say, atoms to clusters of galaxies i.e. from micro to macrocosms. According to current cosmological theories, galaxies are embedded within massive non-baryonic dark matter (DM) halos, which affects their formation and evolution. It is therefore highly rewarding to understand the role of tidal interaction between the dark and luminous matter in galaxies. The current investigation is devoted to Early-Type Galaxies (ETGs), looking in particular at the possibility of establishing whether the tidal interaction of the DM halo with the luminous baryonic component may be at the origin of the so-called "tilt" of the Fundamental Plane (FP). The extension of the tensor virial theorem to two-component matter distributions implies the calculation of the self potential energy due to a selected subsystem, and the tidal potential energy induced by the other one. The additional assumption of homeoidally striated density profiles allows analytical expressions of the results for some cases of astrophysical interest. The current investigation raises from the fact that the profile of the (self + tidal) potential energy of the inner component shows maxima and minima, suggesting the possible existence of preferential scales for the virialized structure, i.e. a viable explanation of the so called "tilt" of the FP. It is found that configurations related to the maxima do not suffice, by themselves, to interpret the FP tilt, and some other relation has to be looked for.

Different sequences of ellipsoids are represented on the ellipticity-rotation plane. The rotation parameter is defined as the ratio of kinetic energy related to the mean tangential equatorial velocity component to kinetic energy related to tangential equatorial component velocity dispersion and residual motions. Systems with isotropic stress tensor are considered as adjoint configurations to their counterparts with anisotropic stress tensor, different angular momentum, and equal remaining parameters. Both nonequilibrium figures and figures elongated by imaginary rotation are represented on the ellipticity-rotation plane. An application is made to a reduced sample of elliptical galaxies. The position on the ellipticity-rotation plane of both sample objects and related adjoint configurations with isotropic stress tensor is inferred from existing observations within the SAURON project. With a single exception, slow rotators are characterized by low ellipticities, low anisotropy parameters, and low rotation parameters, while the contrary holds for fast rotators. A possible interpretation of slow rotators as nonrotating at all and elongated due to negative anisotropy parameters, is exploited. Finally, the elliptical side of the Hubble sequence is interpreted as a sequence of equilibrium (adjoint) configurations where the ellipticity is an increasing function of the rotation parameter, slow rotators correspond to early classes and fast rotators to late classes. In this view, boundaries are rotationally distorted regardless of angular momentum and stress tensor, where rotation has to be intended as due to additional kinetic energy of tangential equatorial velocity components, with respect to spherical configurations with isotropic stress tensor.

The empirical differential oxygen abundance distribution (EDOD) is deduced from subsamples related to two different samples involving solar neighbourhood (SN) thick disk, thin disk, halo, and bulge stars. The EDOD of the SN thick + thin disk is determined by weighting the mass, for assumed SN thick to thin disk mass ratio within the range, 0.1-0.9. Inhomogeneous models of chemical evolution for the SN thick disk, the SN thin disk, the SN thick + thin disk, the SN halo, and the bulge, are computed assuming the instantaneous recycling approximation. The EDOD data are fitted, to an acceptable extent, by their TDOD counterparts provided (i) still undetected, low-oxygen abundance thin disk stars exist, and (ii) a single oxygen overabundant star is removed from a thin disk subsample. In any case, the (assumed power-law) stellar initial mass function (IMF) is universal but gas can be inhibited from, or enhanced in, forming stars at different rates with respect to a selected reference case. Models involving a strictly universal IMF (i.e. gas neither inhibited from, nor enhanced in, forming stars with respect to a selected reference case) can also reproduce the data. The existence of a strictly universal IMF makes similar chemical enrichment within active (i.e. undergoing star formation) regions placed in different environments, but increasing probability of a region being active passing from SN halo to SN thick + thin disk, SN thin disk, SN thick disk, and bulge. On the basis of the results, it is realized that the chemical evolution of the SN thick + thin disk as a whole cannot be excluded.

Two-component, ideal, self-gravitating fluids are conceived as macrogases, and the related equation of state is expressed using the virial theorem for subsystems, under the restriction of homeoidally striated density profiles. Shallower density profiles are found to yield an equation of state, \phi=\phi(y,m), characterized (for assigned values of the fractional mass, m=M_j/ M_i) by the occurrence of two extremum points, a minimum and a maximum. Steeper density profiles produce a similar equation of state, which implies that a special value of m is related to a critical curve where the above mentioned extremum points reduce to a single horizontal inflexion point, and curves below the critical one show no extremum points. The similarity of the isofractional mass curves to van der Waals' isothermal curves, suggests the possibility of a phase transition in a bell-shaped region of the (O y \phi) plane, where the fractional truncation radius along a selected direction is y=R_j/R_i, and the fractional virial potential energy is \phi=(E_{ji})_{vir}/(E_{ij})_{vir}. Further investigation is devoted to mass distributions described by Hernquist (1990) density profiles, for which an additional relation can be used to represent a sample of N=16 elliptical galaxies (EGs) on the (O y \phi) plane, under the assumption that the fractional mass related to EGs and their hosting dark matter (DM) haloes, has a universal value. In the light of the model, the evolution of isolated EGs appears to be other than strictly homologous, due to different values of fractional truncation radii, y, or fractional scaling radii, y^\dagger, deduced from sample objects.

The tensor virial theorem for subsystems is formulated for three-component systems and further effort is devoted to a special case where the inner subsystems and the central region of the outer one are homogeneous, the last surrounded by an isothermal homeoid. The virial equations are explicitly written under additional restrictions. An application is made to hole + vortex + bulge systems, in the limit of flattened inner subsystems. Using the Faber-Jackson relation, the standard $M_{\rm H}$-$\sigma_0$ form is deduced from qualitative considerations. The projected bulge velocity dispersion to projected vortex velocity ratio, $\eta$, as a function of the fractional radius, y_{\rm BV}, and the fractional masses, $m_{\rm BH}$, and $m_{\rm VH}$, is plotted for several cases. It is shown that a fixed value of $\eta$ below the maximum corresponds to two different configurations: a compact bulge on the left and an extended bulge on the right. In addition, for fixed $m_{\rm BH}$ or $m_{\rm BV}$, and $y_{\rm BV}$, more massive bulges are related to larger $\eta$ and vice versa. The model is applied to NGC 4374 and NGC 4486, and the bulge mass is inferred and compared with results from different methods. In presence of a massive vortex $(m_{\rm VH}=5)$, the hole mass has to be reduced by a factor 2-3 with respect to the case of a massless vortex, to get the fit.

The empirical differential metallicity distribution (EDMD) is deduced for (i) local thick disk stars; (ii) likely metal-weak thick disk stars; (iii) chemically selected local G dwarfs, with the corrections performed in order to take into account the stellar scale height; in addition to previous results related to (iv) solar neighbourhood halo subdwarfs; and (v) K-giant bulge stars. The thick disk is conceived as made of two distinct regions: the halo-like and the bulge-like thick disk, and the related EDMD is deduced. Under the assumption that each distribution is typical for the corresponding subsystem, the EDMD of the thick disk, the thick + thin disk, and the Galaxy, is determined by weighting the mass. Models of chemical evolution are computed for each subsystem assuming the instantaneous recycling approximation. The EDMD data are reasonably fitted by simple models implying both homogeneous and inhomogeneous star formation, provided that star formation is inhibited during thick disk evolution. The initial mass function (IMF) is assumed to be a universal power law, which implies an unchanged true yield in different subsystems. The theoretical differential metallicity distribution (TDMD) is first determined for the halo-like thich disk, the bulge-like thick disk, and the thin disk separately, and then for the Galaxy by weighting the mass. An indicative comparison is performed between the EDMD deduced for the disk both in presence and in absence of [O/Fe] plateau, and its counterpart computed for (vi) nearby stars for which the oxygen abundance has been determined both in presence and in absence of the local thermodynamical equilibrium (LTE) approximation. Both distributions are found to exhibit a similar trend, though systematic differences exist.

A theory of collisionless fluids is developed in a unified picture, where nonrotating figures with anisotropic random velocity component distributions and rotating figures with isotropic random velocity component distributions, make adjoints configurations to the same system. R fluids are defined and mean and rms angular velocities and mean and rms tangential velocity components are expressed, by weighting on the moment of inertia and the mass, respectively. The definition of figure rotation is extended to R fluids. The generalized tensor virial equations are formulated for R fluids and further attention is devoted to axisymmetric configurations where, for selected coordinate axes, a variation in figure rotation has to be counterbalanced by a variation in anisotropy excess and vice versa. A microscopical analysis of systematic and random motions is performed under a few general hypotheses, by reversing the sign of tangential or axial velocity components of an assigned fraction of particles, leaving the distribution function and other parameters unchanged (Meza 2002). The application of the reversion process to tangential velocity components, implies the conversion of random motion rotation kinetic energy into systematic motion rotation kinetic energy. The application of the reversion process to axial velocity components, implies the conversion of random motion translation kinetic energy into systematic motion translation kinetic energy, and the loss related to a change of reference frame is expressed in terms of systematic (imaginary) motion rotation kinetic energy. A procedure is sketched for deriving the spin parameter distribution (including imaginary rotation) from a sample of observed or simulated large-scale collisionless fluids i.e. galaxies and galaxy clusters.

Unsustained matter distributions unescapely collapse unless fragmentation and centrifugal or pressure support take place. Starting from the above evidence, supermassive compact objects at the centre of large-mass galaxies are conceived as the end-product of the gravitational collapse of local density maxima around which overdensities are located. At the beginning of evolution, local density maxima are idealized as homogeneous peaks, while the surrounding envelopes are described by a power-law density profile. The dependence of the density profile on a second parameter, chosen to be the ratio between peak and total mass, is analysed. Overdensity evolution is discussed in the context of quintessence cosmological models and further investigation is devoted to a special case with the aim to describe the central collapse. An empirical relation between hole and dark halo mass is translated into a dependence of the fractional hole mass on the overdensity mass. Computations are performed up to the end of central collapse, and density profiles of related configurations are determined together with additional parameters. The central collapse is completed in early times, no longer than a few hundredths of Gyr, which implies hole formation when proto-haloes, proto-bulges, and proto-disks are still expanding or contracting. No appreciable change in evolution is found with regard to different mean peak heights related to equal masses. On the other hand, it is recognized that homogeneous peaks collapse ``faster'' with respect to surroundings envelopes, in low-mass than in large-mass overdensities. In conclusion, it is inferred that gravitational collapse of homogeneous peaks within overdensities may be a viable mechanism for hole generation.

With regard to homeoidally striated Jacobi ellipsoids, a unified theory of systematically rotating and peculiar motions is developed, where both real and imaginary rotation are considered. The effect of positive or negative residual motion excess, is shown to be equivalent to an additional real or imaginary rotation, respectively. Then it is realized that a homeoidally striated Jacobi ellipsoid with assigned velocity distribution, always admits an adjoint configuration i.e. a classical Jacobi ellipsoid of equal mass and axes. In addition, further constraints are established on the amount of residual velocity anisotropy along the principal axes, for triaxial configurations. Special effort is devoted to investigating sequences of virial equilibrium configurations in terms of normalized parameters, including the effects of both density and velocity profile. In particular, it is shown that bifurcation points from axisymmetric to triaxial configurations occur as in classical Jacobi ellipsoids, contrary to earlier results. The reasons of the above mentioned discrepancy are also explained. An interpretation of recent results from numerical simulations on stability (Meza 2002), is provided in the light of the model. A physical interpretation of the early Hubble sequence is shortly reviewed and discussed from the standpoint of the model, according if elliptical galaxies are considered as isolated systems or embedded within dark matter haloes. In any case, a lower limit to the flattening of oblate-like configurations is established. On the other hand, it is found no lower limit to the elongation of prolate-like configurations, and observed lack of elliptical galaxies more elongated than E7 needs a different physical interpretation, such as the fully investigated bending instabilities.

Elliptical galaxies are modelled as homeoidally striated Jacobi ellipsoids where the peculiar velocity distribution is anisotropic, or equivalently as their adjoints configurations i.e. classical Jacobi ellipsoids of equal mass and axes, in real or imaginary rotation. Reasons for the coincidence of bifurcation points from axisymmetric to triaxial configurations in both the sequences, contrary to earlier findings, are presented and discussed. The effect of centrifugal support at the ends of the major equatorial axis, is briefly outlined. The existence of a lower limit to the flattening of elliptical galaxies is investigated in dealing with a number of limiting situations. More specifically, (i) elliptical galaxies are considered as isolated systems, and an allowed region within Ellipsoidland, related to the occurrence of bifurcation points from ellipsoidal to pear-shaped configurations, is shown to be consistent with observations; (ii) elliptical galaxies are considered as embedded within dark matter haloes and, under reasonable assumptions, it is shown that tidal effects from hosting haloes have little influence on the above mentioned results; (iii) dark matter haloes and embedded elliptical galaxies, idealized as a single homeoidally striated Jacobi ellipsoid, are considered in connection with the cosmological transition from expansion to relaxation, by generalizing an earlier model, and the existence of a lower limit to the flattening of relaxed (oblate-like) configurations, is established. On the other hand, no lower limit is found to the elongation of relaxed (prolate-like) configurations, and the observed lack of elliptical galaxies more elongated than E7 needs a different physical interpretation, such as the fully investigated bending instabilities.

With regard to large-scale astrophysical systems, the current paper deals with (i) formulation of tensor virial equations from the standpoint of analytical mechanics; (ii) investigation on the role of systematic and random motions for virial equilibrium configurations; (iii) extent to which systematic and random motions are equivalent in changing a fluid shape. The tensor virial equations are formulated using analytical mechanics, and the self potential-energy tensor is shown to be symmetric. The role of systematic and random motions in collisionless, ideal, self-gravitating fluids, is analysed in detail including radial and tangential velocity dispersion on the equatorial plane. R3 fluids are defined as ideal, self-gravitating fluids in virial equilibrium, with systematic rotation around a principal axis of inertia, and ihe related virial equations are formulated. A unified theory of systematic and random motions is developed for R3 fluids, taking into consideration imaginary rotation. The effect of random motion excess is shown to be equivalent to an additional real or imaginary rotation, respectively, inducing flattening or elongation. R3 fluids are found to admit adjoint configurations with isotropic random velocity distribution. Further constraints are established on the amount of random velocity anisotropy along the principal axes, for triaxial configurations. A necessary condition is formulated for the occurrence of bifurcation points from axisymmetric to triaxial configurations in virial equilibrium, which is independent of the anisotropy parameters. In the special case of homeoidally striated Jacobi ellipsoid, some previously known results are reproduced.

We present a panoramic review of several observational and theoretical aspects of the modern astrophysical research about the origin of the Fundamental Plane (FP) relation for Early-Type Galaxies (ETGs). The discussion is focused on the problem of the tilt and the tightness of the FP, and on the attempts to derive the luminosity evolution of ETGs with redshift. Finally, a number of observed features in the FP are interpreted from the standpoint of a new theoretical approach based on the two-component tensor virial theorem.

Using two alternative [O/H]-[Fe/H] dependences, the empirical oxygen abundance distribution (EGD) is deduced from two different samples of old objects, available in the literature. Under the assumption that each distribution (with the addition of a third one previously known) is typical for the corresponding subsystem, the EGD of the Galactic spheroid is determined weighting by mass. A fit to the data is provided by simple models of chemical evolution implying both homogeneous and inhomogeneous mixing, where star formation is inhibited during halo formation and enhanced during bulge formation, with respect to the disk. The theoretical oxygen abundance distribution (TGD) is first determined for the halo and the bulge separately, and then for the Galactic spheroid weighting by mass. Simple models implying inhomogeneous mixing are able to fit an empirical age-metallicity relation (EAMR) recently determined from a homogeneous sample of globular clusters, which shows a non monotonic trend characterized by large dispersion. Some possible implications related to the formation of Galactic spheroid and disk are also discussed.

The formulation of the tensor virial equations is generalized to unrelaxed configurations, where virial equilibrium does not coincide with dynamical (or hydrostatic) equilibrium. Further investigation is devoted to special classes of homeoidally striated ellipsoids, defined as homeoidally striated, Jacobi ellipsoids. In particular, virial equilibrium configurations with constant anisotropy parameters are studied with more detail, including both flattened and elongated, triaxial configurations, and the determination of the related bifurcation points. The explicit expression of different rotation parameters is also determined. An application is made to dark matter haloes hosting giant, galaxies, with regard to assigned initial and final configuration, following and generalizing to many respects a procedure conceived by Thuan & Gott (1975). The dependence of the limiting axis ratios, below which no configuration is allowed for the sequence under consideration, on the change in mass, total energy, and angular momentum, during the evolution, is illustrated in some representative situations. The dependence of axis ratios and rotation parameters on an additional parameter, related to the initial conditions of the density perturbation, is analysed in connection with a few special cases. Within the range of Peebles (1969) rotation parameter, inferred from high-resolution numerical simulations, the shape of dark matter haloes is mainly decided by the amount of anisotropy in residual velocity distribution. On the other hand, the contribution of rotation has only a minor effect on the meridional plane, and no effect on the equatorial plane, as bifurcation points occur for larger values of Peebles (1969) rotation parameter. To this respect, dark matter haloes are found to resemble giant elliptical galaxies.

Analytical and geometrical properties of generalized power-law (GPL) density profiles are investigated in detail. In particular, a one-to-one correspondence is found between mathematical parameters and geometrical parameters. Then GPL density profiles are compared with simulated dark haloes (SDH) density profiles, and nonlinear least-absolute values and least-squares fits involving the above mentioned five parameters (RFSM5 method) are prescribed. More specifically, the sum of absolute values or squares of absolute logarithmic residuals is evaluated on a large number of points making a 5-dimension hypergrid, through a few iterations. The size is progressively reduced around a fiducial minimum, and superpositions on nodes of earlier hypergrids are avoided. An application is made to a sample of 17 SDHs on the scale of cluster of galaxies, within a flat $\Lambda$CDM cosmological model (Rasia et al. 2004). In dealing with the mean SDH density profile, a virial radius, averaged over the whole sample, is assigned, which allows the calculation of the remaining parameters. Using a RFSM5 method provides a better fit with respect to other methods. No evident correlation is found between SDH dynamical state (relaxed or merging) and asymptotic inner slope of the logarithmic density profile or (for SDH comparable virial masses) scaled radius. Mean values and standard deviations of some parameters are calculated, and a comparison with previous results is made with regard to the scaled radius. A certain degree of degeneracy is found in fitting GPL to SDH density profiles. If it is intrinsic to the RFSM5 method or it could be reduced by the next generation of high-resolution simulations, still remains an open question.

Generalized power-law (GPL) density profiles are compared with simulated dark haloes (SDH) density profiles, and nonlinear least-squares fits are prescribed, involving five parameters which specify the fitting profile (RFSM5 method). More specifically, the validity of a necessary condition for the occurrence of an extremal point, is related to the existence of an intersection between three surfaces in a 3-dimension space. An application is made to a sample of 17 SDHs on the scale of cluster of galaxies, within a flat /\CDM cosmological model. In dealing with the averaged SDH density profile (ADP), a virial radius, equal to the mean over the whole sample, is assigned, which allows the calculation of the remaining parameters. The exponents of both the best fitting GPL density profile to ADP and related averages calculated over the whole halo sample, are far from their NFW counterparts. Values of asymptotic inner slope of fitting logarithmic density profiles, are consistent with results from recent high-resolution simulations.

Simulated dark matter haloes are fitted by self-similar, universal density profiles, where the scaled parameters depend only on a scaled (truncation) radius which, in turn, is supposed to be independent on the mass and the formation redshift. A criterion for the choice of the best fitting density profile is proposed, with regard to a set of high-resolution simulations, where some averaging procedure on scaled density profiles has been performed, in connection with a number of fitting density profiles. An application is made to a pair of sets each made of a dozen of high-resolution simulations, which are available in literature, in connection with two currently used fitting density profiles, where the dependence of the scaled radius on the mass and the formation redshift, may be neglected to a first extent. Some features of the early evolution of dark matter haloes represented by fitting density profiles, are discussed in the limit of the spherical top-hat model.

A general theory of homeoidally striated density profiles where no divergence occurs, is adapted to cuspy density profiles, with a suitable choice of the scaling density and the scaling radius. A general formulation of some physical parameters, such as angular-momentum vector, rotational-energy tensor (both calculated in connection with a special class of rotational velocity fields), inertia tensor, and self potential-energy tensor, is performed. Other potential-energy tensors involving two density profiles where the boundaries are similar and similarly placed, are also expressed. Explicit results are attained for three special cases of physical interest which fit to a good extent the results of high-resolution simulations for dark matter haloes, and a class of density profiles, which closely approximate the de Vaucouleurs $r^{1/4}$ law for elliptical galaxies. The virial theorem in tensor form for two-component systems is written for each subsystem, and applied to giant elliptical galaxies. The predicted velocity dispersion along the line of sight, in the limiting case where a principal axis points towards the observer, is found to be consistent with observations except for (intrinsic) $E7$ configurations where the major axis points towards the observer. If dark matter haloes host an amount of undetected baryons about twice as massive as the stellar subsystem, and undetected baryons trace non baryonic matter therein, two main consequences arise, namely (i) velocity dispersions along the line of sight are lower than in absence of undetected baryons, and (ii) dark matter haloes are dynamically ``hotter'' than stellar ellipsoids, the transition occurring when the amount of undetected baryons is about one and a half times that of the stellar subsystem.