We develop a machine learning (ML) framework to populate large dark matter-only simulations with baryonic galaxies. Our ML framework takes input halo properties including halo mass, environment, spin, and recent growth history, and outputs central galaxy and halo baryonic properties including stellar mass ($M_*$), star formation rate (SFR), metallicity ($Z$), neutral ($\rm HI$) and molecular ($\rm H_2$) hydrogen mass. We apply this to the MUFASA cosmological hydrodynamic simulation, and show that it recovers the mean trends of output quantities with halo mass highly accurately, including following the sharp drop in SFR and gas in quenched massive galaxies. However, the scatter around the mean relations is under-predicted. Examining galaxies individually, at $z=0$ the stellar mass and metallicity are accurately recovered ($\sigma\lesssim 0.2$~dex), but SFR and $\rm HI$ show larger scatter ($\sigma\gtrsim 0.3$~dex); these values improve somewhat at $z=1,2$. Remarkably, ML quantitatively recovers second parameter trends in galaxy properties, e.g. that galaxies with higher gas content and lower metallicity have higher SFR at a given $M_*$. Testing various ML algorithms, we find that none perform significantly better than the others, nor does ensembling improve performance, likely because none of the algorithms reproduce the large observed scatter around the mean properties. For the random forest algorithm, we find that halo mass and nearby ($\sim 200$~kpc) environment are the most important predictive variables followed by growth history, while halo spin and $\sim$Mpc scale environment are not important. Finally we study the impact of additionally inputting key baryonic properties $M_*$, SFR and $Z$, as would be available e.g. from an equilibrium model, and show that particularly providing the SFR enables $\rm HI$ to be recovered substantially more accurately.

We extend a data-based model-free multifractal method of exoplanet detection to probe exoplanetary atmospheres. Whereas the transmission spectrum is studied during the primary eclipse, we analyze the emission spectrum during the secondary eclipse, thereby probing the atmospheric limb. In addition to the spectral structure of exoplanet atmospheres, the approach provides information to study phenomena such as atmospheric flows, tidal-locking behavior, and the dayside-nightside redistribution of energy. The approach is demonstrated using Spitzer data for exoplanet HD189733b. The central advantage of the method is the lack of model assumptions in the detection and observational schemes.

Supernova cosmology without spectra will be an important component of future surveys such as LSST. This lack of supernova spectra results in uncertainty in the redshifts which, if ignored, leads to significantly biased estimates of cosmological parameters. Here we present a hierarchical Bayesian formalism -- zBEAMS -- that addresses this problem by marginalising over the unknown or uncertain supernova redshifts to produce unbiased cosmological estimates that are competitive with supernova data with spectroscopically confirmed redshifts. zBEAMS provides a unified treatment of both photometric redshifts and host galaxy misidentification (occurring due to chance galaxy alignments or faint hosts), effectively correcting the inevitable contamination in the Hubble diagram. Like its predecessor BEAMS, our formalism also takes care of non-Ia supernova contamination by marginalising over the unknown supernova type. We illustrate this technique with simulations of supernovae with photometric redshifts and host galaxy misidentification. A novel feature of the photometric redshift case is the important role played by the redshift distribution of the supernovae.

We use state-of-art measurements of the galaxy luminosity function (LF) at z=6, 7 and 8 to derive constraints on warm dark matter (WDM), late-forming dark matter (LFDM) and ultra-light axion dark matter (ULADM) models alternative to the cold dark matter (CDM) paradigm. To this purpose we have run a suite of high-resolution N-body simulations to accurately characterise the low mass-end of the halo mass function and derive DM model predictions of the high-z luminosity function. In order to convert halo masses into UV-magnitudes we introduce an empirical approach based on halo abundance matching which allows us to model the LF in terms of the amplitude and scatter of the ensemble average star formation rate halo mass relation of each DM model, $\langle {\rm SFR}({\rm M_{ h}},z)\rangle$. We find that independent of the DM scenario the average SFR at fixed halo mass increases from z=6 to 8, while the scatter remains constant. At halo mass ${\rm M_{h}}\gtrsim 10^{12}\,{\rm M}_\odot$ h$^{-1}$ the average SFR as function of halo mass follows a double power law trend that is common to all models, while differences occur at smaller masses. In particular, we find that models with a suppressed low-mass halo abundance exhibit higher SFR compared to the CDM results. Using deviance statistics we obtain a lower limit on the WDM thermal relic particle mass, $m_{\rm WDM}\gtrsim 1.5$ keV at $2\sigma$. In the case of LFDM models, the phase transition redshift parameter is bounded to $z_t\gtrsim 8\cdot 10^5$ at $2\sigma$. We find ULADM best-fit models with axion mass $m_a\gtrsim 1.6\cdot 10^{-22}$ eV to be well within $2\sigma$ of the deviance statistics. We remark that measurements at $z=6$ slightly favour a flattening of the LF at faint UV-magnitudes. This tends to prefer some of the non-CDM models in our simulation suite, although not at a statistically significant level to distinguish them from CDM.

Owing to technological advances, the number of exoplanets discovered has risen dramatically in the last few years. However, when trying to observe Earth analogs, it is often difficult to test the veracity of detection. We have developed a new approach to the analysis of exoplanetary spectral observations based on temporal multifractality, which identifies time scales that characterize planetary orbital motion around the host star, and those that arise from stellar features such as spots. Without fitting stellar models to spectral data, we show how the planetary signal can be robustly detected from noisy data using noise amplitude as a source of information. For observation of transiting planets, combining this method with simple geometry allows us to relate the time scales obtained to primary and secondary eclipse of the exoplanets. Making use of data obtained with ground-based and space-based observations we have tested our approach on HD 189733b. Moreover, we have investigated the use of this technique in measuring planetary orbital motion via Doppler shift detection. Finally, we have analyzed synthetic spectra obtained using the SOAP 2.0 tool, which simulates a stellar spectrum and the influence of the presence of a planet or a spot on that spectrum over one orbital period. We have demonstrated that, so long as the signal-to-noise-ratio $\ge$ 75, our approach reconstructs the planetary orbital period, as well as the rotation period of a spot on the stellar surface.

We perform a study of the nonlinear clustering of matter in the late-forming dark matter (LFDM) scenario in which dark matter results from the transition of a nonminimally coupled scalar field from radiation to collisionless matter. A distinct feature of this model is the presence of a damped oscillatory cutoff in the linear matter power spectrum at small scales. We use a suite of high-resolution N-body simulations to study the imprints of LFDM on the nonlinear matter power spectrum, the halo mass and velocity functions and the halo density profiles. The model largely satisfies high-redshift matter power spectrum constraints from Lyman-$\alpha$ forest measurements, while it predicts suppressed abundance of low-mass halos ($\sim 10^{9}-10^{10}$ h$^{-1}$ M$_\odot$) at all redshifts compared to a vanilla $\Lambda$CDM model. The analysis of the LFDM halo velocity function shows a better agreement than the $\Lambda$CDM prediction with the observed abundance of low-velocity galaxies in the local volume. Halos with mass $M\gtrsim 10^{11}$ h$^{-1}$ M$_\odot$ show minor departures of the density profiles from $\Lambda$CDM expectations, while smaller-mass halos are less dense, consistent with the fact that they form later than their $\Lambda$CDM counterparts.

Artificial fragmentation of the matter density field causes the formation of spurious groups of particles in N-body simulations of non-standard Dark Matter (DM) models which are characterized by a small scale cut-off in the linear matter power spectrum. These spurious halos alter the prediction of the mass function in a range of masses where differences among DM models are most relevant to observational tests. Using a suite of high resolution simulations we show that the contamination of artificial groups of particles significantly affect the statistics of halo spin, shape and virial state parameters. We find that spurious halos have systematically larger spin values, are highly elliptical or prolate and significantly deviate from virial equilibrium. These characteristics allow us to detect the presence of spurious halos even in non-standard DM models for which the low-mass end of the mass function remains well behaved. We show that selecting halos near the virial equilibrium provides a simple and effective method to remove the bulk of spurious halos from numerical halo catalogs and consistently recover the halo mass function at low masses.

We present the first report on an innovative new project named "RAD@home", a citizen-science research collaboratory built on free web-services like Facebook, Google, Skype, NASA Skyview, NED, TGSS etc.. This is the first of its kind in India, a zero-funded, zero-infrastructure, human-resource network to educate and directly involve in research, hundreds of science-educated under-graduate population of India, irrespective of their official employment and home-location with in the country. Professional international collaborators are involved in follow up observation and publication of the objects discovered by the collaboratory. We present here ten newly found candidate episodic radio galaxies, already proposed to GMRT, and ten more interesting cases which includes, bent-lobe radio galaxies located in new Mpc-scale filaments, likely tracing cosmological cluster accretion from the cosmic web. Two new Speca-like rare spiral-host large radio galaxies have also been been reported here. Early analyses from our follow up observations with the Subaru and XMM-Newton telescopes have revealed that Speca is likely a new entry to the cluster and is a fast rotating, extremely massive, star forming disk galaxy. Speca-like massive galaxies with giant radio lobes, are possibly remnants of luminous quasars in the early Universe or of first supermassive black holes with in first masssve galaxies. As discoveries of Speca-like galaxies did not require new data from big telescopes, but free archival radio-optical data, these early results demonstrate the discovery potential of RAD@home and how it can help resource-rich professionals, as well as demonstrate a model of academic-growth for resource-poor people in the underdeveloped regions via Internet.

In this paper we introduce PkANN, a freely available software package for interpolating the non-linear matter power spectrum, constructed using Artificial Neural Networks (ANNs). Previously, using Halofit to calculate matter power spectrum, we demonstrated that ANNs can make extremely quick and accurate predictions of the power spectrum. Now, using a suite of 6380 N-body simulations spanning 580 cosmologies, we train ANNs to predict the power spectrum over the cosmological parameter space spanning $3\sigma$ confidence level (CL) around the concordance cosmology. When presented with a set of cosmological parameters ($\Omega_{\rm m} h^2, \Omega_{\rm b} h^2, n_s, w, \sigma_8, \sum m_\nu$ and redshift $z$), the trained ANN interpolates the power spectrum for $z\leq2$ at sub-per cent accuracy for modes up to $k\leq0.9\,h\textrm{Mpc}^{-1}$. PkANN is faster than computationally expensive N-body simulations, yet provides a worst-case error $<1$ per cent fit to the non-linear matter power spectrum deduced through N-body simulations. The overall precision of PkANN is set by the accuracy of our N-body simulations, at 5 per cent level for cosmological models with $\sum m_\nu<0.5$ eV for all redshifts $z\leq2$. For models with $\sum m_\nu>0.5$ eV, predictions are expected to be at 5 (10) per cent level for redshifts $z>1$ ($z\leq1$). The PkANN interpolator may be freely downloaded from http://zuserver2.star.ucl.ac.uk/~fba/PkANN

We calculate the cosmic Mach number M - the ratio of the bulk flow of the velocity field on scale R to the velocity dispersion within regions of scale R. M is effectively a measure of the ratio of large-scale to small-scale power and can be a useful tool to constrain the cosmological parameter space. Using a compilation of existing peculiar velocity surveys, we calculate M and compare it to that estimated from mock catalogues extracted from the LasDamas (a LCDM cosmology) numerical simulations. We find agreement with expectations for the LasDamas cosmology at ~ 1.5 sigma CL. We also show that our Mach estimates for the mocks are not biased by selection function effects. To achieve this, we extract dense and nearly-isotropic distributions using Gaussian selection functions with the same width as the characteristic depth of the real surveys, and show that the Mach numbers estimated from the mocks are very similar to the values based on Gaussian profiles of the corresponding widths. We discuss the importance of the survey window functions in estimating their effective depths. We investigate the nonlinear matter power spectrum interpolator PkANN as an alternative to numerical simulations, in the study of Mach number.

The estimation and analysis of large-scale bulk flow moments of peculiar velocity surveys is complicated by non-spherical survey geometry, the non-uniform sampling of the matter velocity field by the survey objects and the typically large measurement errors of the measured line-of-sight velocities. Previously, we have developed an optimal `minimum variance' (MV) weighting scheme for using peculiar velocity data to estimate bulk flow moments for idealized, dense and isotropic surveys with Gaussian radial distributions, that avoids many of these complications. These moments are designed to be easy to interpret and are comparable between surveys. In this paper, we test the robustness of our MV estimators using numerical simulations. Using MV weights, we estimate the bulk flow moments for various mock catalogues extracted from the LasDamas and the Horizon Run numerical simulations and compare these estimates to the moments calculated directly from the simulation boxes. We show that the MV estimators are unbiased and negligibly affected by non-linear flows.

We investigate the interpolation of power spectra of matter fluctuations using Artificial Neural Network (PkANN). We present a new approach to confront small-scale non-linearities in the power spectrum of matter fluctuations. This ever-present and pernicious uncertainty is often the Achilles' heel in cosmological studies and must be reduced if we are to see the advent of precision cosmology in the late-time Universe. We show that an optimally trained artificial neural network (ANN), when presented with a set of cosmological parameters (Omega_m h^2, Omega_b h^2, n_s, w_0, sigma_8, m_nu and redshift z), can provide a worst-case error <=1 per cent (for z<=2) fit to the non-linear matter power spectrum deduced through N-body simulations, for modes up to k<=0.7 h/Mpc. Our power spectrum interpolator is accurate over the entire parameter space. This is a significant improvement over some of the current matter power spectrum calculators. In this paper, we detail how an accurate interpolation of the matter power spectrum is achievable with only a sparsely sampled grid of cosmological parameters. Unlike large-scale N-body simulations which are computationally expensive and/or infeasible, a well-trained ANN can be an extremely quick and reliable tool in interpreting cosmological observations and parameter estimation. This paper is the first in a series. In this method paper, we generate the non-linear matter power spectra using HaloFit and use them as mock observations to train the ANN. This work sets the foundation for Paper II, where a suite of N-body simulations will be used to compute the non-linear matter power spectra at sub-per cent accuracy, in the quasi-non-linear regime 0.1 h/Mpc <= k <= 0.9 h/Mpc. A trained ANN based on this N-body suite will be released for the scientific community.

The peculiar velocities of galaxies are an inherently valuable cosmological probe, providing an unbiased estimate of the distribution of matter on scales much larger than the depth of the survey. Much research interest has been motivated by the high dipole moment of our local peculiar velocity field, which suggests a large scale excess in the matter power spectrum, and can appear to be in some tension with the LCDM model. We use a composite catalogue of 4,537 peculiar velocity measurements with a characteristic depth of 33 h-1 Mpc to estimate the matter power spectrum. We compare the constraints with this method, directly studying the full peculiar velocity catalogue, to results from Macaulay et al. (2011), studying minimum variance moments of the velocity field, as calculated by Watkins, Feldman & Hudson (2009) and Feldman, Watkins & Hudson (2010). We find good agreement with the LCDM model on scales of k > 0.01 h Mpc-1. We find an excess of power on scales of k < 0.01 h Mpc-1, although with a 1 sigma uncertainty which includes the LCDM model. We find that the uncertainty in the excess at these scales is larger than an alternative result studying only moments of the velocity field, which is due to the minimum variance weights used to calculate the moments. At small scales, we are able to clearly discriminate between linear and nonlinear clustering in simulated peculiar velocity catalogues, and find some evidence (although less clear) for linear clustering in the real peculiar velocity data.

We investigate the impact of massive neutrinos on the distribution of matter in the semi-non-linear regime (0.1<k<0.6 h/Mpc). We present a suite of large-scale N-body simulations quantifying the scale dependent suppression of the total matter power spectrum, resulting from the free-streaming of massive neutrinos out of high-density regions. Our simulations show a power suppression of 3.5-90 per cent at k~0.6 h/Mpc for total neutrino mass, m_nu=0.05-1.9 eV respectively. We also discuss the precision levels that future cosmological datasets would have to achieve in order to distinguish the normal and inverted neutrino mass hierarchies.

We apply the theory of algebraic polynomials to analytically study the transonic properties of general relativistic hydrodynamic axisymmetric accretion onto non-rotating astrophysical black holes. For such accretion phenomena, the conserved specific energy of the flow, which turns out to be one of the two first integrals of motion in the system studied, can be expressed as a 8$^{th}$ degree polynomial of the critical point of the flow configuration. We then construct the corresponding Sturm's chain algorithm to calculate the number of real roots lying within the astrophysically relevant domain of $\mathbb{R}$. This allows, for the first time in literature, to {\it analytically} find out the maximum number of physically acceptable solution an accretion flow with certain geometric configuration, space-time metric, and equation of state can have, and thus to investigate its multi-critical properties {\it completely analytically}, for accretion flow in which the location of the critical points can not be computed without taking recourse to the numerical scheme. This work can further be generalized to analytically calculate the maximal number of equilibrium points certain autonomous dynamical system can have in general. We also demonstrate how the transition from a mono-critical to multi-critical (or vice versa) flow configuration can be realized through the saddle-centre bifurcation phenomena using certain techniques of the catastrophe theory.

By applying the theory of algebraic polynomials and the theory of dynamical systems, we construct the generalized Sturm sequences/chains to investigate the transonic properties of hydrodynamic accretion onto non-rotating astrophysical black holes, to demonstrate, completely analytically, how many critical point such an accretion flow can have. Our work is significantly important, because for the first time in the literature, we provide a purely analytical method, by applying certain powerful theorem of algebraic polynomials in pure mathematics, to check whether certain astrophysical hydrodynamic accretion may undergo more than one sonic transitions. Our work can be generalized to analytically calculate the maximal number of equilibrium points certain autonomous dynamical systems can have in general (Abridged).