This paper proposes k-medoids clustering method to reveal the distinct groups of 1,318 variable stars in the Galaxy based on their light curves, where each light curve represents the graph of brightness of the star against time. To overcome the deficiencies of subjective traditional classification, we separate the stars more scientifically according to their geometrical configuration and show that our approach outperforms the existing classification schemes in astronomy. It results in two optimum groups of eclipsing binaries corresponding to bright, massive systems and fainter, less massive systems.

Formation mechanism of present day population of elliptical galaxies have been revisited in the context of hierarchical cosmological models accompanied by accretion and minor mergers through cross correlation function including spatial effect. The present work investigates the formation and evolution of several components of nearby massive early type galaxies (ETGs) through cross-correlation in the spatial coordinates, right ascension and declination (RA, DEC) and mass-size parameter space with high redshift $(0.5\leq z\leq2.7)$ ETGs. It is found that innermost components of nearby ETGs are highly correlated with ETGs in the redshift range $(2\leq z\leq2.7)$ known as 'red nuggets'. The intermediate and outermost parts have moderate correlations with ETGs in the redshift range $(0.5\leq z\leq0.75)$. The quantitative measures are highly consistent with the two phase formation scenario of massive nearby early type galaxies as suggested by various authors and resolves the conflict raised in a previous work suggesting other possibilities for the formation of outermost part of nearby massive ETGs. The improvement is expected to be due to inclusion of spatial effects in addition to other linear parameters.

We consider the problem related to clustering of gamma-ray bursts (from "BATSE" catalogue) through kernel principal component analysis in which our proposed kernel outperforms results of other competent kernels in terms of clustering accuracy and we obtain three physically interpretable groups of gamma-ray bursts. The effectivity of the suggested kernel in combination with kernel principal component analysis in revealing natural clusters in noisy and nonlinear data while reducing the dimension of the data is also explored in two simulated data sets.

Clustering objects into synthetic groups is a natural activity of any science. Astrophysics is not an exception and is now facing a deluge of data. For galaxies, the one-century old Hubble classification and the Hubble tuning fork are still largely in use, together with numerous mono-or bivariate classifications most often made by eye. However, a classification must be driven by the data, and sophisticated multivariate statistical tools are used more and more often. In this paper we review these different approaches in order to situate them in the general context of unsupervised and supervised learning. We insist on the astrophysical outcomes of these studies to show that multivariate analyses provide an obvious path toward a renewal of our classification of galaxies and are invaluable tools to investigate the physics and evolution of galaxies.

Theoretical and indirect observational evidences suggest that stellar initial mass function (IMF) increases with redshift. On the other hand star formation rates (SFR) may be as high as 100 $M_{\odot}$ yr$^{-1}$ in star burst galaxies. These may lead to formation of massive clusters hence massive stars to make the integrated galactic stellar initial mass function (IGIMF) top heavy (i.e. proportion of massive stars is higher than less massive stars). We investigate the joint effect of evolving IMF and several measures of SFR in dependence of galaxy wide IMF. The resulting IGIMF have slopes $\alpha_{2,IGIMF}$ in the high mass regime, which is highly dependent on the minimum mass of the embedded cluster ($M_{ecl,min}$), star formation rates and mass spectrum indices of embedded clusters (viz. $\beta$). It is found that for z $\sim$ 0 - 2, $\alpha_{2,IGIMF}$ becomes steeper (i.e. bottom heavy), for z $\sim$ 2 - 4, $\alpha_{2,IGIMF}$ becomes flatter (i.e. top heavy ) and from z $\sim$ 4 onwards $\alpha_{2,IGIMF}$ becomes again steeper. The effects are faster for higher values of $\beta$. $\alpha_{2,IGIMF}$ is flatter also for higher values of $M_{ecl,min}$. All these effects might be counted for the joint effect of increasing temperature of the ambient medium as well as varying SFR with increasing redshift.

Origin of linear instability resulting in rotating sheared accretion flows has remained a controversial subject for long. While some explanations of such non-normal transient growth of disturbances in the Rayleigh stable limit were available for magnetized accretion flows, similar instabilities in absence of magnetic perturbations remained unexplained. This dichotomy was resolved in two recent publications by Chattopadhyay, {\it et al} where it was shown that such instabilities, especially for non-magnetized accretion flows, were introduced through interaction of the inherent stochastic noise in the system (even a \enquote{cold} accretion flow at 3000K is too \enquote{hot} in the statistical parlance and is capable of inducing strong thermal modes) with the underlying Taylor-Couette flow profiles. Both studies, however, excluded the additional energy influx (or efflux) that could result from nonzero cross-correlation of a noise perturbing the velocity flow, say, with the noise that is driving the vorticity flow (or equivalently the magnetic field and magnetic vorticity flow dynamics). In this article we show that nonzero noise cross-correlations essentially renormalize the strength of temporal correlations. Apart from an overall boost in the energy rate (both for spatial and temporal correlations, and hence in the ensemble averaged energy spectra), this results in mutual competition in growth rates of affected variables often resulting in suppression of oscillating Alfven waves at small times while leading to faster saturations at relatively longer time scales. The effects are seen to be more pronounced with magnetic field fluxes where the noise cross-correlation magnifies the strength of the field concerned.

Spatial clustering nature of galaxies have been studied previously through auto correlation function. The same type of cross correlation function has been used to investigate parametric clustering nature of galaxies e.g. with respect to masses and sizes of galaxies. Here formation and evolution of several components of nearby massive early type galaxies have been envisaged through cross correlations, in the mass-size parametric plane, with high redshift early type galaxies (hereafter ETG).It is found that the inner most components of nearby ETGs have significant correlation with ETGs in the highest redshift range called red nuggets whereas intermediate components are highly correlated with ETGs in the redshift range with z value greater than 0.5 and less than 0.75. The outer most part has no correlation in any range, suggesting a scenario through in situ accretion. The above formation scenario is consistent with the previous results obtained for NGC5128 (Chattopadhyay et al. (2009), Chattopadhyay et al. (2013)) and to some extent for nearby elliptical galaxies (Huang et al. (2013)) after considering a sample of ETGs at high redshift with stellar masses greater than or equal to 108.73 M-Sun. So the present work indicates a three phase formation of massive nearby elliptical galaxies instead of two as discussed in previous works.

Cluster analysis is the distribution of objects into different groups or more precisely the partitioning of a data set into subsets (clusters) so that the data in subsets share some common trait according to some distance measure. Unlike classi cation, in clustering one has to rst decide the optimum number of clusters and then assign the objects into different clusters. Solution of such problems for a large number of high dimensional data points is quite complicated and most of the existing algorithms will not perform properly. In the present work a new clustering technique applicable to large data set has been used to cluster the spectra of 702248 galaxies and quasars having 1540 points in wavelength range imposed by the instrument. The proposed technique has successfully discovered ve clusters from this 702248X1540 data matrix.

We investigate the evolution of magnetohydrodynamic/hydromagnetic perturbations in the presence of stochastic noise in rotating shear flows. The particular emphasis is the flows whose angular velocity decreases but specific angular momentum increases with increasing radial coordinate. Such flows, however, are Rayleigh stable, but must be turbulent in order to explain astrophysical observed data and, hence, reveal a mismatch between the linear theory and observations/experiments. The mismatch seems to have been resolved, atleast in certain regimes, in the presence of weak magnetic field revealing magnetorotational instability. The present work explores the effects of stochastic noise on such magnetohydrodynamic flows, in order to resolve the above mismatch generically for the hot flows. We essentially concentrate on a small section of such a flow which is nothing but a plane shear flow supplemented by the Coriolis effect, mimicking a small section of an astrophysical accretion disk around a compact object. It is found that such stochastically driven flows exhibit large temporal and spatial auto-correlations and cross-correlations of perturbation and hence large energy dissipations of perturbation, which generate instability. Interestingly, auto-correlations and cross-correlations appear independent of background angular velocity profiles, which are Rayleigh stable, indicating their universality. This work, to the best of our knowledge, is the first attempt to understand the evolution of three-dimensional hydromagnetic perturbations in rotating shear flows in the presence of stochastic noise.

Origin of hydrodynamic turbulence in rotating shear flows is investigated. The particular emphasis is the flows whose angular velocity decreases but specific angular momentum increases with increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to explain observed data. Such a mismatch between the linear theory and observations/experiments is more severe when any hydromagnetic/magnetohydrodynamic instability and then the corresponding turbulence therein is ruled out. The present work explores the effect of stochastic noise on such hydrodynamic flows. We essentially concentrate on a small section of such a flow which is nothing but a plane shear flow supplemented by the Coriolis effect. This also mimics a small section of an astrophysical accretion disk. It is found that such stochastically driven flows exhibit large temporal and spatial correlations of perturbation velocities, and hence large energy dissipations of perturbation, which presumably generate instability. A range of angular velocity (\Omega) profiles of background flow, starting from that of constant specific angular momentum (\lambda = \Omega r^2 ; r being the radial coordinate) to that of constant circular velocity (v_\phi = \Omega r), is explored. However, all the background angular velocities exhibit identical growth and roughness exponents of perturbations, revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear flows. This work, to the best of our knowledge, is the first attempt to understand origin of instability and turbulence in the three-dimensional Rayleigh stable rotating shear flows by introducing additive noise to the underlying linearized governing equations. This has important implications to resolve the turbulence problem in astrophysical hydrodynamic flows such as accretion disks.

Galaxy diversification proceeds by transforming events like accretion, interaction or mergers. These explain the formation and evolution of galaxies that can now be described with many observables. Multivariate analyses are the obvious tools to tackle the datasets and understand the differences between different kinds of objects. However, depending on the method used, redundancies, incompatibilities or subjective choices of the parameters can void the usefulness of such analyses. The behaviour of the available parameters should be analysed before an objective reduction of dimensionality and subsequent clustering analyses can be undertaken, especially in an evolutionary context. We study a sample of 424 early-type galaxies described by 25 parameters, ten of which are Lick indices, to identify the most structuring parameters and determine an evolutionary classification of these objects. Four independent statistical methods are used to investigate the discriminant properties of the observables and the partitioning of the 424 galaxies: Principal Component Analysis, K-means cluster analysis, Minimum Contradiction Analysis and Cladistics. (abridged)

We present a statistical analysis of the properties of a large sample of dynamically hot old stellar systems, from globular clusters to giant ellipticals, which was performed in order to investigate the origin of ultra-compact dwarf galaxies. The data were mostly drawn from Forbes et al. (2008). We recalculated some of the effective radii, computed mean surface brightnesses and mass-to-light-ratios, estimated ages and metallicities. We completed the sample with globular clusters of M31. We used a multivariate statistical technique (K-Means clustering), together with a new algorithm (Gap Statistics) for finding the optimum number of homogeneous sub-groups in the sample, using a total of six parameters (absolute magnitude, effective radius, virial mass-to-light ratio, stellar mass-to-light ratio and metallicity). We found six groups. FK1 and FK5 are composed of high- and low-mass elliptical galaxies respectively. FK3 and FK6 are composed of high-metallicity and low-metallicity objects, respectively, and both include globular clusters and ultra-compact dwarf galaxies. Two very small groups, FK2 and FK4, are composed of Local Group dwarf spheroidals. Our groups differ in their mean masses and virial mass-to-light ratios. The relations between these two parameters are also different for the various groups. The probability density distributions of metallicity for the four groups of galaxies is similar to that of the globular clusters and UCDs. The brightest low-metallicity globular clusters and ultra-compact dwarf galaxies tend to follow the mass-metallicity relation like elliptical galaxies. The objects of FK3 are more metal-rich per unit effective luminosity density than high-mass ellipticals.

The fundamental plane of early-type galaxies is a rather tight three-parameter correlation discovered more than twenty years ago. It has resisted a both global and precise physical interpretation despite a consequent number of works, observational, theoretical or using numerical simulations. It appears that its precise properties depend on the population of galaxies in study. Instead of selecting a priori these populations, we propose to objectively construct homologous populations from multivariate analyses. We have undertaken multivariate cluster and cladistic analyses of a sample of 56 low-redshift galaxy clusters containing 699 early-type galaxies, using four parameters: effective radius, velocity dispersion, surface brightness averaged over effective radius, and Mg2 index. All our analyses are consistent with seven groups that define separate regions on the global fundamental plane, not across its thickness. In fact, each group shows its own fundamental plane, which is more loosely defined for less diversified groups. We conclude that the global fundamental plane is not a bent surface, but made of a collection of several groups characterizing several fundamental planes with different thicknesses and orientations in the parameter space. Our diversification scenario probably indicates that the level of diversity is linked to the number and the nature of transforming events and that the fundamental plane is the result of several transforming events. We also show that our classification, not the fundamental planes, is universal within our redshift range (0.007 - 0.053). We find that the three groups with the thinnest fundamental planes presumably formed through dissipative (wet) mergers. In one of them, this(ese) merger(s) must have been quite ancient because of the relatively low metallicity of its galaxies, Two of these groups have subsequently undergone dry mergers to increase their masses. In the k-space, the third one clearly occupies the region where bulges (of lenticular or spiral galaxies) lie and might also have formed through minor mergers and accretions. The two least diversified groups probably did not form by major mergers and must have been strongly affected by interactions, some of the gas in the objects of one of these groups having possibly been swept out. The interpretation, based on specific assembly histories of galaxies of our seven groups, shows that they are truly homologous. They were obtained directly from several observables, thus independently of any a priori classification. The diversification scenario relating these groups does not depend on models or numerical simulations, but is objectively provided by the cladistic analysis. Consequently, our classification is more easily compared to models and numerical simulations, and our work can be readily repeated with additional observables.

An objective classification of the globular clusters of NGC 5128 has been carried out by using a model-based approach of cluster analysis. The set of observable parameters includes structural parameters, spectroscopically determined Lick indices and radial velocities from the literature. The optimum set of parameters for this type of analysis is selected through a modified technique of Principal Component Analysis, which differs from the classical one in the sense that it takes into consideration the effects of outliers present in the data. Then a mixture model based approach has been used to classify the globular clusters into groups. The efficiency of the techniques used is tested through the comparison of the misclassification probabilities with those obtained using the K-means clustering technique. On the basis of the above classification scheme three coherent groups of globular clusters have been found. We propose that the clusters of one group originated in the original cluster formation event that coincided with the formation of the elliptical galaxy, and that the clusters of the two other groups are of external origin, from tidally stripped dwarf galaxies on random orbits around NGC 5128 for one group, and from an accreted spiral galaxy for the other.

Two different multivariate clustering techniques, the K-means partitioning method and the Dirichlet process of mixture modeling, have been applied to the BATSE Gamma-ray burst (GRB) catalog, to obtain the optimum number of coherent groups. In the standard paradigm, GRB are classified in only two groups, the long and short bursts. However, for both the clustering techniques, the optimal number of classes was found to be three, a result which is consistent with previous statistical analysis. In this classification, the long bursts are further divided into two groups which are primarily differentiated by their total fluence and duration and hence are named low and high fluence GRB. Analysis of GRB with known red-shifts and spectral parameters suggests that low fluence GRB have nearly constant isotropic energy output of 10^{52} ergs while for the high fluence ones, the energy output ranges from 10^{52} to 10^{54} ergs. It is speculated that the three kinds of GRBs reflect three different origins: mergers of neutron star systems, mergers between white dwarfs and neutron stars, and collapse of massive stars.