We develop a Lagrangian Perturbation Theory (LPT) framework to study the clustering of cold dark matter (CDM) in cosmologies with massive neutrinos. We follow the trajectories of CDM particles with Lagrangian displacements fields up to third order in perturbation theory. Once the neutrinos become non-relativistic, their density fluctuations are modeled as being proportional to the CDM density fluctuations, with a scale-dependent proportionality factor. This yields a gravitational back-reaction that introduces additional scales to the linear growth function, which is accounted for in the higher order LPT kernels. Through non-linear mappings from Eulerian to Lagrangian frames, we ensure that our theory has a well behaved large scale behavior free of unwanted UV divergences, which are common when neutrino and CDM densities are not treated on an equal footing, and in resummation schemes that manifestly break Galilean invariance. We use our theory to construct correlation functions for both the underlying matter field, as well as for biased tracers using Convolution-LPT. Redshift-space distortions effects are modeled using the Gaussian Streaming Model. When comparing our analytical results to simulated data from the Quijote simulation suite, we find good accuracy down to $r=20 \,\text{Mpc} \, h^{-1}$ at redshift $z=0.5$, for the real space and redshift space monopole particle correlation functions with no free parameters. The same accuracy is reached for the redshift space quadrupole if we additionally consider an effective field theory parameter that shifts the pairwise velocity dispersion. For modeling the correlation functions of tracers we adopt a simple Lagrangian biasing scheme with only density and curvature operators, which we find sufficient to reach down to $r=20 \,\text{Mpc} \, h^{-1}$ when comparing to simulated halos.

The Galactic microquasar GRS 1915+105 shows rich variability which is categorized into different classes. In this paper, we report temporal and spectral analysis of GRS 1915+105 to study the properties of the accretion flow when the light curve is showing $\theta$ class variability. For this purpose, we use the LAXPC data from the Target of Opportunity observations of India's first multi-wavelength astronomy satellite AstroSat. The $\theta$ class is marked by the recurrent appearance of U-shaped regions in the light curve, where the photon count rate first decreases rapidly and then increases slowly. For our analysis, we use U-shaped regions of first two orbits (02345 and 02346) on 2016 March 04. In both of the cases, the dynamic Power Density Spectra (PDS) showed significant power at around $4-5$ Hz, suggesting the presence of a low-frequency Quasi-Periodic Oscillation (QPO) around that frequency interval. The QPO frequency is found to increase with time when the enhancement of the energy flux also takes place. From the evolution of the spectra, we determine the evolution of the accretion flow parameters in both of these observations. Fitting the spectra with the transonic flow solution based Two-Component Advective Flow (TCAF) model in $4-25$ keV energy band shows that the Keplerian disk accretion rate rises with the rise in the radiation intensity, while the location of the centrifugal pressure driven shock front decreases. In both these data, a gradual increment of power-law photon index with intensity is observed suggesting the progressive softening of the source.

We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, $w$, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of $w$ that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of $w$. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of $\sim 200\ \rm{km\ s}^{-1}$, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

Galactic transient black hole candidate (BHC) MAXI~J1535-571 was discovered on 2017 September 02 simultaneously by {\it MAXI}/GSC and {\it Swift}/BAT instruments. It has also been observed by India's first multi-wavelength astronomy-mission satellite {\it AstroSat}, during the rising phase of its 2017-18 outburst. We make both the spectral and the temporal analysis of the source during 2017 September 12-17 using data of {\it AstroSat}'s Large Area X-ray Proportional Counter (LAXPC) in the energy range of $3-40$~keV to infer the accretion flow properties of the source. Spectral analysis is done with the physical two-component advective flow (TCAF) solution-based {\it fits} file. From the nature of the variation of the TCAF model fitted physical flow parameters, we conclude and confirm that during our analysis period source was in the soft-intermediate spectral state. We observe sharp type-C quasi-periodic oscillations (QPOs) in the frequency range of $\sim 1.74-2.81$~Hz, which are not common in the soft-intermediate spectral state. For a better understanding of the nature and evolution of these type-C QPOs, a dynamic study of the power density spectra is done. We also investigate the origin of these QPOs from the shock oscillation model. We find that as is expected in the soft-Intermediate spectral states, non-satisfaction of Rankine-Hugoniot conditions for non-dissipative shocks and not their resonance oscillations may be the cause of the origin of type-C QPOs.

The class variable source GRS 1915+105 exhibits a wide range of time variabilities in time scales of a few seconds to a few days. Depending on the count rates in different energy bands and the nature of the conventional color-color diagram, the variabilities were classified into sixteen classes, which were later sequenced in ascending order of Comptonization Efficiency (CE). The CE is the ratio of the phenomenological power-law and disk blackbody model fitted spectral photon fluxes. For a better understanding of the flow dynamics in one of the variability classes, namely $\chi$, we fitted spectra with the physical two-component advective flow (TCAF) model. In $\chi$ class, X-ray flux is steady with no significant variation, however, various $\chi$ subclasses are observed at different X-ray flux and variations of count rates across different $\chi$ subclasses must be linked to the variation of flow parameters such as the accretion rates, be it the Keplerian disc rate and/or the low angular momentum halo rate. We find that in the $\chi_{2,4}$ classes, which are reportedly devoid of significant outflows, the spectra could be fitted well using TCAF solution alone. In the $\chi_{1,3}$ class, which are always linked with outflows, a cutoff power-law model is needed in addition to the TCAF solution. At the same time, the normalization required by this model along with the variation of photon index and exponential roll-off factor provides us with the information on the relative dominance of the outflow in the latter two classes. TCAF fit also provides us with the size and location of the Compton cloud along with its optical depth. Thus by fitting with TCAF, a physical understanding of the flow geometry in different $\chi$ classes of GRS 1915+105 has been obtained.

Local determinations of the Hubble constant $H_0$ favour a higher value than Planck based on CMB and $\Lambda$CDM. Through a model-independent expansion, we show that low redshift ($z \lesssim 0.7$) data comprising baryon acoustic oscillations (BAO), cosmic chronometers and Type Ia supernovae has a preference for Quintessence models that lower $H_0$ relative to $\Lambda$CDM. In addition, we confirm that an exponential coupling to dark matter cannot alter this conclusion in the same redshift range. Our results leave open the possibility that a coupling in the matter-dominated epoch, potentially even in the dark ages, may yet save $H_0$ from sinking in the string theory Swampland.

We present the results of a detailed investigation of the prompt and afterglow emission in the HESS detected GRB 190829A. Swift and Fermi observations of the prompt phase of this GRB reveal two isolated sub-bursts or episodes, separated by a quiescent phase. The energetic and the spectral properties of the first episode are in stark contrast to the second. The first episode, which has a higher spectral peak of $\sim 120\:\text{keV}$ and a low isotropic energy $\sim 10^{50}\:\text{erg}$ is an outlier to the Amati correlation and marginally satisfies the Yonetoku correlation. However, the energetically dominant second episode has lower peak energy and is consistent with the above correlations. We compared this GRB to other low luminosity GRBs (LLGRBs). Prompt emission of LLGRBs also indicates a relativistic shock breakout origin of the radiation. For GRB 190829A, some of the properties of a shock breakout origin are satisfied. However, the absence of an accompanying thermal component and energy above the shock breakout critical limit precludes a shock breakout origin. In the afterglow, an unusual long-lasting late time flare of duration $\sim 10^4\:\text{s}$ is observed. We also analyzed the late-time \fermi-LAT emission that encapsulates the H.E.S.S. detection. Some of the LAT photons are likely to be associated with the source. All the above observational facts suggest GRB 190829A is a peculiar low luminosity GRB that is not powered by a shock breakout, and with an unusual rebrightening due to a patchy emission or a refreshed shock during the afterglow. Furthermore, our results show that TeV energy photons seem common in both high luminosity GRBs and LLGRBs.

Low surface brightness (LSBs) superthins constitute classic examples of very late-type galaxies, with their disc dynamics strongly regulated by their dark matter halos. In this work we consider a gravitational origin of dark matter in the brane world scenario, where the higher dimensional Weyl stress term projected onto the 3-brane acts as the source of dark matter. In the context of the braneworld model, this dark matter is referred to as the \emph{`dark mass'}.This model has been successful in reproducing the rotation curves of several low surface brightness and high surface brightness galaxies. Therefore it is interesting to study the prospect of this model in explaining the vertical structure of galaxies which has not been explored in the literature so far. Using our 2-component model of gravitationally-coupled stars and gas in the external force field of this \emph{dark mass}, we fit the observed scale heights of stellar and atomic hydrogen (HI) gas of superthin galaxy `UGC7321' using the Markov Chain Monte Carlo approach. We find that the observed scaleheights of `UGC7321' can be successfully modelled in the context of the braneworld scenario. In addition, the model predicted rotation curve also matches the observed one. The implications on the model parameters are discussed.

We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the non linear regime. We also investigate extensions of the Schr\"odinger-Poisson system that employ multiple stacked cold streams, and the von Neumann-Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. We present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.

Superthin galaxies are a class of low surface brightness galaxies with strikingly high values of planar-to-vertical axes ratio $\rm(b/a> 10 - 20)$ with little or no bulge component, possibly indicating the presence of an ultra-cold stellar disc, the origin, and evolution of which is not well understood. Using the 2-component model of gravitationally-coupled stars and gas in the external force field of a dark matter halo and assumed to be in vertical hydrostatic equilibrium, we determine the vertical velocity dispersion of stars and gas as a function of galactocentric radius for a sample of five superthin galaxies (UGC 7321, IC 5249, FGC 1540, IC2233 and UGC711) using observed stellar and atomic hydrogen gas (HI) scaleheights as constraints, and employing the Markov Chain Monte Carlo (MCMC) method. The mass models for our sample galaxies, constructed using stellar photometry and HI 21cm radio-synthesis observations were already available in the literature. We find that the central vertical velocity dispersion for the stellar disc in the optical band varies between ${\sigma}_0$ $\sim$ $10.2 - 18.4 km/s$ and falls off with an exponential scale length of $2.6$ to $3.2 R_{D}$ where $R_{D}$ is the exponential stellar disc scale length. Interestingly, in the $3.6 \mu$m, the same, averaged over the two components of the stellar disc, varies between $5.9$ to $11.8 km/s$, mainly representative of the denser, thinner and smaller of the two-disc components. However, the dispersion of the more massive disc component varies between 15.9 - 24.7 with a scalelength of $\sim$ 2.2 $R_{D}$. Our calculated values of the N-component disc stability parameter $Q_{min}$ lies between 1.7 - 5.7, thus confirming the dynamic stability of our model. Further, our model results are consistent with AGAMA(Action-based Galaxy Modelling Architecture).

The origin of the bulk of the astrophysical neutrinos detected by the IceCube Observatory remains a mystery. Previous source-finding analyses compare the directions of IceCube events and individual sources in astrophysical catalogs. The source association method is technically challenging when the number of source candidates is much larger than the number of the observed astrophysical neutrinos. We show that in this large source number regime, a cross-correlation analysis of neutrino data and source catalog can instead be used to constrain potential source populations for the high-energy astrophysical neutrinos, and provide spatial evidence for the existence of astrophysical neutrinos. We present an analysis of the cross-correlation of the IceCube 2010-2012 point-source data and a WISE-2MASS galaxy sample. While we find no significant detection of cross-correlation with the publicly available neutrino dataset, we show that, when applied to the full IceCube data, which has a longer observation time and higher astrophysical neutrino purity, our method has sufficient statistical power to detect a cross-correlation signal if the neutrino sources trace the Large Scale Structure of the Universe.

Constructing dynamical models for interacting pair of galaxies as constrained by their observed structure and kinematics crucially depends on the correct choice of the values of the relative inclination ($i$) between their galactic planes as well as the viewing angle ($\theta$), the angle between the line of sight and the normal to the plane of their orbital motion. We construct Deep Convolutional Neural Network (DCNN) models to determine the relative inclination ($i$) and the viewing angle ($\theta$) of interacting galaxy pairs, using N-body $+$ Smoothed Particle Hydrodynamics (SPH) simulation data from the GALMER database for training the same. In order to classify galaxy pairs based on their $i$ values only, we first construct DCNN models for a (a) 2-class ( $i$ = 0 $^{\circ}$, 45$^{\circ}$ ) and (b) 3-class ($i = 0^{\circ}, 45^{\circ} \text{ and } 90^{\circ}$) classification, obtaining $F_1$ scores of 99% and 98% respectively. Further, for a classification based on both $i$ and $\theta$ values, we develop a DCNN model for a 9-class classification ($(i,\theta) \sim (0^{\circ},15^{\circ}) ,(0^{\circ},45^{\circ}), (0^{\circ},90^{\circ}), (45^{\circ},15^{\circ}), (45^{\circ}, 45^{\circ}), (45^{\circ}, 90^{\circ}), (90^{\circ}, 15^{\circ}), (90^{\circ}, 45^{\circ}), (90^{\circ},90^{\circ})$), and the $F_1$ score was 97$\%$. Finally, we tested our 2-class model on real data of interacting galaxy pairs from the Sloan Digital Sky Survey (SDSS) DR15, and achieve an $F_1$ score of 78%. Our DCNN models could be further extended to determine additional parameters needed to model dynamics of interacting galaxy pairs, which is currently accomplished by trial and error method.

Recently constructed high-redshift Hubble diagrams of supernovae, quasars and gamma-ray bursts show $\sim4 \, \sigma$ deviations from flat $\Lambda$CDM, where the discrepancies are observed in cosmographic expansions. In this note, we show that existing approaches to high-redshift cosmography can lead to "tensions" that are simply an artifact of the expansions, even if the underlying model is fully consistent with flat $\Lambda$CDM. Thus, in this note we revisit the $\sim 4 \, \sigma$ tensions to ascertain if results have been impacted by the cosmography. By directly fitting respective datasets to flat $\Lambda$CDM we find that the tension is real, but the quasar data, which has been calibrated by external Type Ia supernovae has a 7:3 preference for a Universe with no dark energy over one described by $\Lambda$CDM with canonical (Planck) values. Results based on cosmography mask this feature.

We discuss the sensitivity of the present and near-future axion dark matter experiments to a halo of axions or axion-like particles gravitationally bound to the Earth or the Sun. The existence of such halos, assuming they are formed, renders a significant gain in the sensitivity of axion searches while satisfying all the present experimental bounds. The structure and coherence properties of these halos also imply novel signals, which can depend on the latitude or orientation of the detector. We demonstrate this by analysing the sensitivity of several distinct types of axion dark matter experiments.

The cosmological relaxion can address the hierarchy problem, while its coherent oscillations can constitute dark matter in the present universe. We consider the possibility that the relaxion forms gravitationally bound objects that we denote as relaxion stars. The density of these stars would be higher than that of the local dark matter density, resulting in enhanced signals in table-top detectors, among others. Furthermore, we raise the possibility that these objects may be trapped by an external gravitational potential, such as that of the Earth or the Sun. This leads to formation of relaxion halos of even greater density. We discuss several interesting implications of relaxion halos, as well as detection strategies to probe them.

We comprehensively analyse the cosmology dependence of counts-in-cell statistics. We focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly nonlinear scales. Based on large-deviation statistics, we parametrise the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the nonlinear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass $M_\nu$ and can disentangle it from the clustering amplitude $\sigma_8$. Using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the $\nu\Lambda$CDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first $\nu\Lambda$CDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly nonlinear scales. The matter PDF is particularly powerful for constraining the matter density $\Omega_m$, clustering amplitude $\sigma_8$ and the total neutrino mass $M_\nu$. Adding the mildly nonlinear matter PDF to the mildly nonlinear matter power spectrum improves constraints on $\Omega_m$ by a factor of 5 and $\sigma_8$ by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc/h)$^3$. We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

We show for the first time that the loop-driven kinetic mixing between visible and dark Abelian gauge bosons can facilitate dark matter production in the early Universe by creating a 'dynamic' portal, which depends on the energy of the process. The required smallness of the strength of the portal interaction, suited for freeze-in, is justified by a suppression arising from the mass of a heavy vector-like fermion. The strong temperature sensitivity associated with the interaction is responsible for most of the dark matter production during the early stages of reheating.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

Models of the spatial distribution of dark matter halos must achieve new levels of precision and accuracy in order to satisfy the requirements of upcoming experiments. In this work, we present a halo bias emulator for modeling the clustering of halos on large scales. It incorporates the cosmological dependence of the bias beyond the mapping of halo mass to peak height. The emulator makes substantial improvements in accuracy compared to the widely used Tinker et al. (2010) model. Halos in this work are defined using an overdensity criteria of 200 relative to the mean background density. Halo catalogs are produced for 40 N-body simulations as part of the Aemulus project at snapshots from z=3 to z=0. The emulator is trained over the mass range $6\times10^{12}-7\times10^{15}\ h^{-1}M_{\odot}$. Using an additional suite of 35 simulations, we determine that the precision of the emulator is redshift dependent, achieving sub-percent levels for a majority of the redshift range. Two additional simulation suites are used to test the ability of the emulator to extrapolate to higher and lower masses. Our high-resolution simulation suite is used to develop an extrapolation scheme in which the emulator asymptotes to the Tinker et al. (2010) model at low mass, achieving ~3% accuracy down to $10^{11}\ h^{-1}M_{\odot}$. Finally, we present a method to propagate emulator modeling uncertainty into an error budget. Our emulator is made publicly available at \url{https://github.com/AemulusProject/bias_emulator}.

Nonlinear objects like halos and voids exhibit a scale-dependent bias on linear scales in massive neutrino cosmologies. The shape of this scale-dependent bias is a unique signature of the neutrino masses, but the amplitude of the signal is generally small, of the order of $f_\nu$, the contribution of neutrinos to the total matter content ($\lesssim 1\%$). In this paper, we demonstrate for the first time how the strength of this signal can be substantially enhanced by using information about the halo environment at a range of scales. This enhancement is achieved by using certain combinations of the large scale Cold Dark Matter and total matter environments of halos, both of which are measurable from galaxy clustering and weak lensing surveys.

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($\sigma/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.

Thawing and freezing quintessence models are compared thermodynamically. Both of them are found to disobey the Generalized Second Law of Thermodynamics. However, for freezing models, there is still a scope as this breakdown occurs in the past, deep inside the radiation dominated era, when a standard scalar field model with a pressureless matter is not a correct description of the matter content. The thawing model has a pathological breakdown in terms of thermodynamics in a finite future.

Superthin galaxies are low surface brightness (LSB) bulgeless disc galaxies having stellar discs with unusually high planar-to-vertical axes ratio $b/a$ $>$ 10 - 20, the formation and evolution of which is not well-understood. We calculate the specific angular momenta of a sample of six superthins and nine other bulgeless LSBs using stellar photometry, atomic hydrogen (HI) surface density and high resolution HI rotation curves available in the literature. We find that the stellar specific angular mometum $j_{\rm{s}}$, and hence the stellar disc size given by the exponential stellar disc scale length R$_D$, of three superthins and seven LSBs lie above the 95.4 $\%$ confidence band of the $j_{\rm{s}}$ - $V_{\rm{rot}}$ regression line for ordinary bulgeless disc galaxies, $V_{\rm{rot}}$ being the asymptotic rotational velocity. Further, we find that superthins and LSBs have higher $j_{\rm{s}}$ and R$_D$ values for a given value of stellar mass $M_{\rm{s}}$ at high values of statistical significance, compared to ordinary disc galaxies. Therefore, we conclude, a superthin is may be distinguished by a characteristically larger disc size which could possibly explain the origin of its large planar-to-vertical axes ratio. Interestingly, we find that the median spin parameter $\lambda = \frac{ j_{\rm{stars}} } { {\sqrt(2) V_{\rm{vir}} R_{\rm{vir}} } }$, $V_{\rm{vir}}$ and $R_{\rm{vir}}$ being the virial velocity and virial radius of the galaxy respectively, is 0.13 $\pm$ 0.01 for superthin galaxies which is an order of magnitude higher than those of LSBs and ordinary disc galaxies, which may have important implications for the existence of superthin stellar discs in these low surface brightness galaxies.

GRB 190114C is a famous Gamma-Ray Burst (GRB) due to its detection at sub-TeV energies by MAGIC, seen at redshift z = 0.42. This burst is one of the brightest GRB detected by \fermi. We present a detailed analysis of GRB 190114C prompt emission, using the two \fermi detectors: GBM and LAT. The \emph{LAT low energy events} (LLE) data is also considered. A joint GBM-LAT analysis reveals a sub-$GeV$ spectral cutoff. A similar high energy cutoff was reported in GRB 160509A and GRB 100724B earlier, and a handful of other sources. The cutoff can be explained by the intrinsic opacity due to pair production within the emitting region. Such morphology in these GRBs suggests that they belong to one specific class having a similar source of the radiation mechanism. GRB 190114C shows a transition from non-thermal to a quasi-thermal-like spectrum along with radiation due to external shock. From spectrum analysis and $Lorentz$ factor evolution from the trigger time to late emission, considering the fact that sub-TeV photons are detected in MAGIC, we are able to draw an emission mechanism picture, where the prompt emission spectrum is more consistent with spectrum via photospheric dissipation with presence of external shock emission simultaneously.

We present the UNIT $N$-body cosmological simulations project, designed to provide precise predictions for nonlinear statistics of the galaxy distribution. We focus on characterizing statistics relevant to emission line and luminous red galaxies in the current and upcoming generation of galaxy surveys. We use a suite of precise particle mesh simulations (FastPM) as well as with full $N$-body calculations with a mass resolution of $\sim 1.2\times10^9\,h^{-1}$M$_{\odot}$ to investigate the recently suggested technique of Angulo & Pontzen 2016 to suppress the variance of cosmological simulations We study redshift space distortions, cosmic voids, higher order statistics from $z=2$ down to $z=0$. We find that both two- and three-point statistics are unbiased. Over the scales of interest for baryon acoustic oscillations and redshift-space distortions, we find that the variance is greatly reduced in the two-point statistics and in the cross correlation between halos and cosmic voids, but is not reduced significantly for the three-point statistics. We demonstrate that the accuracy of the two-point correlation function for a galaxy survey with effective volume of 20 ($h^{-1}$Gpc)$^3$ is improved by about a factor of 40, indicating that two pairs of simulations with a volume of 1 ($h^{-1}$Gpc)$^3$ lead to the equivalent variance of $\sim$150 such simulations. The $N$-body simulations presented here thus provide an effective survey volume of about seven times the effective survey volume of DESI or Euclid. The data from this project, including dark matter fields, halo catalogues, and their clustering statistics, are publicly available at http://www.unitsims.org.

The galactic black hole candidate Cygnus X-1, one of the brightest sources in the sky, is the first ever black hole candidate to be discovered. Despite being a very well-studied object due to its persistent brightness in X-rays, there has been much difficulty in determining its mass since its discovery. In spite of continuous efforts since the early nineteen seventies, there is yet no concensus on its mass for nearly four decades. The uncertainties in the mass measurements are due to the high degree of error involved in the measurement of its distance. In 2011, Orosz et al. constrained the mass of this object to be M = 14.8 +/- 1.0 M_Sun using dynamical methods. In this paper, we obtained the mass of Cygnus X-1, using a completely independent method, namely, carrying out the spectral analysis using Two Component Advective Flow (TCAF) solution based fits file and the archival data of RXTE PCA instrument. Our result does not require the distance of the source or the information about the companion. Each spectral fit with the TCAF gives one best fitted mass. Averaging fitted masses of Cygnus X-1 over a span of five months of observation during its persistent hard phase, mass of the source comes out to be M_avg = 14.20 +/- 0.36 M_Sun, which is consistent with the dynamically estimated mass.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The black hole transient H1743-322 exhibited several outbursts with temporal and spectral variability since its discovery in 1977. These outbursts occur at a quasi-regular recurrence period of around $0.5-2$ years, since its rediscovery in March 2003. We investigate accretion flow dynamics around the Low Mass X-ray Binary H1743-322 during its 2004 outburst using the RXTE/PCA archival data. We use Two Component Advective Flow (TCAF) solution to analyse the spectral data. From the fits with TCAF solution, we obtain day to day variation of physical accretion rates of Keplerian and sub-Keplerian components, size of the Compton cloud and its other properties. Analysis of the spectral properties of the 2004 outburst by keeping fitted normalization to be in a narrow range and its timing properties in terms of the presence and absence of QPOs, enable us to constrain the mass of the black hole in a range of $10.31 M_{\odot} - 14.07 M_{\odot}$ which is consistent with other estimates reported in the literature.

We present high resolution H{\sc i} 21cm Giant Meterwave Radio Telescope (GMRT) observations of the superthin galaxy FGC1540 with a spatial resolution of 10$''$ $\times$ 8$''$ and a spectral resolution of 1.73 kms$^{-1}$ and an rms noise of 0.9 mJy per beam. We obtain its rotation curve as well as deprojected radial H{\sc i} surface density profile by fitting a 3-dimensional tilted ring model directly to the H{\sc i} data cubes by using the publicly-available software, Fully Automated Tirrific (FAT). We also present the rotation curve of FGC1540 derived from its optical spectroscopy study using the 6-m BTA telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. We use the rotation curve, the H{\sc i} surface density profile together with Spitzer 3.6 $\mu$m and the SDSS $i$--band data to construct the mass models for FGC1540. We find that both the Pseudo-isothermal (PIS), as well as Navarro-Frenk-White (NFW) dark matter (DM) halos, fit the observed rotation curve equally well. The PIS model indicates a compact dark matter halo ($R_{\rm C}/R_{\rm D}$ < 2), with the best-fitting core radius ($R_{\rm C}$) approximately half the exponential stellar disc scale length ($R_{\rm D}$), which is in agreement with the mass models of superthin galaxies studied earlier in the literature. Since the vertical thickness of the galactic stellar disc is determined by a balance between the net gravitational field and the velocity dispersion in the vertical direction, the compact dark matter halo may be primarily responsible in regulating the superthin vertical structure of the stellar disc in FGC1540 as was found in case of the superthin galaxy UGC7321.

Self-Interacting Dark Matter (SIDM) has long been proposed as a solution to small scale problems posed by standard Cold Dark Matter (CDM). We use numerical simulations to study the effect of dark matter interactions on the morphology of disk galaxies falling into galaxy clusters. The effective drag force on dark matter leads to offsets of the stellar disk with respect to the surrounding halo, causing distortions in the disk. For anisotropic scattering cross-sections of 0.5 and 1.0$\,\textrm{cm}^{2}\textrm{g}^{-1}$, we show that potentially observable warps, asymmetries, and thickening of the disk occur in simulations. We discuss observational tests of SIDM with galaxy surveys and more realistic simulations needed to obtain detailed predictions.

In this article we perform a detailed theoretical analysis for a class of new exact solutions with anisotropic fluid distribution of matter for compact objects in hydrostatic equilibrium. To achieve this we call the relation between the metric functions, namely, embedding class one condition. The investigation is carried out by generalising the properties of a spherical star with an emphasis on hydrostatic equilibrium equation, i.e., the generalised Tolman-Oppenheimer-Volkoff equation, in our understanding of these compact objects. We match the interior solution to an exterior Reissner-Nordstrom solution, and study some physical features of this models, such as the energy conditions, speeds of sound, and mass - radius relation of the star. We also show that obtained solution is compatible with observational data for compact object Her X-1.

Late-type galaxies falling into a cluster would evolve being influenced by the interactions with both the cluster and the nearby cluster member galaxies. Most numerical studies, however, tend to focus on the effects of the former with little work done on those of the latter. We thus perform a numerical study on the evolution of a late-type galaxy interacting with neighboring early-type galaxies at high speed, using hydrodynamic simulations. Based on the information obtained from the Coma cluster, we set up the simulations for the case where a Milky Way-like late-type galaxy experiences six consecutive collisions with twice as massive early-type galaxies having hot gas in their halos at the closest approach distances of 15-65 kpc/h at the relative velocities of 1500-1600 km/s. Our simulations show that the evolution of the late-type galaxy can be significantly affected by the accumulated effects of the high-speed multiple collisions with the early-type galaxies, such as on cold gas content and star formation activity of the late-type galaxy, particularly through the hydrodynamic interactions between cold disk and hot gas halos. We find that the late-type galaxy can lose most of its cold gas after the six collisions and have more star formation activity during the collisions. By comparing our simulation results with those of galaxy-cluster interactions, we claim that the role of the galaxy-galaxy interactions on the evolution of late-type galaxies in clusters could be comparable with that of the galaxy-cluster interactions, depending on the dynamical history.

We present a new method for generating initial conditions for numerical cosmological simulations in which massive neutrinos are treated as an extra set of N-body (collisionless) particles. It allows us to accurately follow the density field for both Cold Dark Matter (CDM) and neutrinos at both high and low redshifts. At high redshifts, the new method is able to reduce the shot noise in the neutrino power spectrum by a factor of more than $10^7$ compared to previous methods, where the power spectrum was dominated by shot noise at all scales. We find that our new approach also helps to reduce the noise on the total matter power spectrum on large scales, whereas on small scales the results agree with previous simulations. Our new method also allows for a systematic study of clustering of the low velocity tail of the distribution function of neutrinos. This method also allows for the study of the evolution of the overall velocity distribution as a function of the environment determined by the CDM field.

We investigate for an Effective Field Theory (EFT) framework that can consistently explain inflation to Large Scale Structures (LSS). With the development of the construction algorithm of EFT, we arrive at a properly truncated action for the entire scenario. Using this, we compute the two-point correlation function for quantum fluctuations from Goldstone modes and related inflationary observables in terms of coefficients of relevant EFT operators, which we constrain using Planck 2015 data. We then carry forward this primordial power spectrum with the same set of EFT parameters to explain the linear and non-linear regimes of LSS by loop-calculations of the matter overdensity two-point function. For comparative analysis, we make use of two widely accepted transfer functions, namely, BBKS and Eisenstein-Hu, thereby making the analysis robust. We finally corroborate our results with LSS data from SDSS-DR7 and WiggleZ. The analysis thus results in a consistent, model-independent EFT framework for inflation to structures.

(abridged) We investigate the signatures left by the cosmic neutrino background on the clustering of matter, CDM+baryons and halos in redshift-space using a set of more than 1000 N-body and hydrodynamical simulations with massless and massive neutrinos. We find that the effect neutrinos induce on the clustering of CDM+baryons in redshift-space on small scales is almost entirely due to the change in $\sigma_8$. Neutrinos imprint a characteristic signature in the quadrupole of the matter (CDM+baryons+neutrinos) field on small scales, that can be used to disentangle the effect of $\sigma_8$ and $M_\nu$. We show that the effect of neutrinos on the clustering of halos is very different, on all scales, to the one induced by $\sigma_8$. We find that the effects of neutrinos of the growth rate of CDM+baryons ranges from $\sim0.3\%$ to $2\%$ on scales $k\in[0.01, 0.5]~h{\rm Mpc}^{-1}$ for neutrinos with masses $M_\nu \leqslant 0.15$ eV. We compute the bias between the momentum of halos and the momentum of CDM+baryon and find it to be 1 on large scales for all models with massless and massive neutrinos considered. This point towards a velocity bias between halos and total matter on large scales that it is important to account for in order to extract unbiased neutrino information from velocity/momentum surveys such as kSZ observations. We show that baryonic effects can affect the clustering of matter and CDM+baryons in redshift-space by up to a few percent down to $k=0.5~h{\rm Mpc}^{-1}$. We find that hydrodynamics and astrophysical processes, as implemented in our simulations, only distort the relative effect that neutrinos induce on the anisotropic clustering of matter, CDM+baryons and halos in redshift-space by less than $1\%$. Thus, the effect of neutrinos in the fully non-linear regime can be written as a transfer function with very weak dependence on astrophysics.

Low Surface Brighness Galaxies (LSBs), inspite of being gas rich, have low star formation rates and are, therefore, low surface brightness in nature. We calculate Q$_{\rm{RW}}$, the 2-component disc stability parameter as proposed by Romeo \& Wiegert (2011), as a function of galactocentric radius $R$ for a sample of five LSBs, for which mass models, as obtained from HI 21cm radio-sythesis observations and R-band photometry, were available in the literature. We find that the median value of Q$_{\rm{RW}}^{\rm{min}}$, the minimum of Q$_{\rm{RW}}$ over $R$, lies between 2.6 and 3.1 for our sample LSBs, which is higher than the median value of 1.8 $\pm$ 0.3 for Q$_{\rm{RW}}^{\rm{min}}$ for a sampleof high surface brightness galaxies (HSBs) as obtained in earlier studies. This clearly shows that LSBs have more stable discs than HSBs, which could explain their low star formation rates and, possibly, their low surface brightness nature. Interestingly, the calculated values of Q$_{\rm{RW}}$ decrease only slightly (median Q$_{\rm{RW}}^{\rm{min}}$ $\sim$ 2.3 - 3) if the discs were taken to respond to the gravitational potential of the dark matter halo only, but reduce by $\sim$ a factor of 2-3 (median Q$_{\rm{RW}}^{\rm{min}}$ $\sim$ 0.7 - 1.5) if they respond to their self-gravity alone. This implies that the dark matter halo is crucial in regulating disc stability in LSBs, which may have important implications for models of galaxy formation and evolution.

Superthin galaxies are low surface brightness disc galaxies, characterised by optical discs with strikingly high values of planar-to-vertical axes ratios ($>$ 10), the physical origin and evolution of which continue to be a puzzle. We present mass models for three superthin galaxies: IC5249, UGC7321 and IC2233. We use high resolution rotation curves and gas surface density distributions obtained from HI 21 cm radio-synthesis observations, in combination with their two-dimensional structural surface brightness decompositions at Spitzer 3.6 $\mu$m band, all of which were available in the literature. We find that while models with the pseudo-isothermal (PIS) and the Navarro-Frenk-White (NFW) dark matter density profiles fit the observed rotation curves of IC5249 and UGC7321 equally well, those with the NFW profile does not comply with the slowly-rising rotation curve of IC2233. Interestingly, for all of our sample galaxies, the best-fitting mass models with a PIS dark matter density profile indicate a {\it compact} dark matter halo i.e., $R_c/R_D$ $<$ 2 where $R_c$ is the core radius of the PIS dark matter halo, and $R_D$ is the radial scale-length of the exponential stellar disc. The compact dark matter halo may be fundamentally responsible for the superthin nature of the stellar disc, and therefore our results may have important implications for the formation and evolution models of superthin galaxies in the universe.

We analyze the ability of galaxy and CMB lensing surveys to constrain massive neutrinos and new models of dark radiation. We present a Fisher forecast analysis for neutrino mass constraints with the LSST galaxy survey and the CMB S4 survey. A joint analysis of the three galaxy and shear 2-point functions, along with key systematics parameters and Planck priors, constrains the neutrino masses to $\sum m_\nu = 0.041\,$eV at 1-$\sigma$ level, comparable to constraints expected from Stage 4 CMB lensing. If low redshift information from upcoming spectroscopic surveys like DESI is included, the constraint becomes $\sum m_\nu = 0.032\,$eV. These constraints are derived having marginalized over the number of relativistic species ($N_{\rm eff}$), which is somewhat degenerate with the neutrino mass. We also explore the gain by combining LSST and CMB S4, that is, using the five relevant auto- and cross-correlations of the two datasets. We conclude that advances in modeling the nonlinear regime and the measurements of other parameters are required to ensure a neutrino mass detection. Using the same datasets, we explore the ability of LSST-era surveys to test "nonstandard" models with dark radiation. We find that if evidence for dark radiation is found from $N_{\rm eff}$ measurements, the mass of the dark radiation candidate can be measured at a 1-$\sigma$ level of $0.162\,$eV for fermionic dark radiation, and $0.137\,$eV for bosonic dark radiation, for $\Delta N_{\rm eff} = 0.15$. We also find that the NNaturalness model of Arkani-Hamed et al 2016, with extra light degrees of freedom, has a sub-percent effect on the power spectrum: even more ambitious surveys than the ones considered here will be needed to test such models.

We present a new method for simulating cosmologies that contain massive particles with thermal free streaming motion, such as massive neutrinos or warm/hot dark matter. This method combines particle and fluid descriptions of the thermal species to eliminate the shot noise known to plague conventional N-body simulations. We describe this method in detail, along with results for a number of test cases to validate our method, and check its range of applicability. Using this method, we demonstrate that massive neutrinos can produce a significant scale-dependence in the large-scale biasing of deep voids in the matter field. We show that this scale-dependence may be quantitatively understood using an extremely simple spherical expansion model which reproduces the behavior of the void bias for different neutrino parameters.

The Cosmic Microwave Background can provide information regarding physics of the very early universe, more specifically, of the matter-radiation distribution of the inflationary era. Starting from the effective field theory of inflation, we use the Goldstone action to calculate the three point correlation function for the Goldstone field, whose results can be directly applied to the field describing the curvature perturbations around a de Sitter solution for the inflationary era. We then use the data from the recent Planck mission for the parameters $f_{NL}^{equil}$ and $f_{NL}^{orthog}$ which parametrize the size and shape of non-Gaussianities generated in single field models of inflation. Using these known values, we calculate the parameters relevant to our analysis, $f_{NL}^{\dot{\pi}^3}$, $f_{NL}^{\dot{\pi}(\partial _i \pi)^2}$ and the speed of sound $c_s$ which parametrize the non-Gaussianities arising from two different kinds of generalized interactions of the scalar field in question.

We model the observed vertical distribution of the neutral hydrogen (HI) in the faint (M_B ~ -13.7 mag) edge-on dwarf irregular galaxy KK250. Our model assumes that the galaxy consists of axi-symmetric, co-planar gas and stellar discs in the external force-field of a spherical dark matter halo, and in vertical hydrostatic equilibrium. The velocity dispersion of the gas is left as a free parameter in the model. Our best fit model is able to reproduce the observed vertical distribution of the HI gas, as well as the observed velocity profiles. The best fit model has a large velocity dispersion (~ 22 km/s) at the centre of the galaxy, which falls to a value of ~ 8 km/s by a galacto-centric radius of ~ 1 kpc, which is similar to both the scale-length of the stellar disc, as well as the angular resolution of the data along the radial direction. Similarly we find that the thickness of the HI disc is also minimum at ~ 1 kpc, and increases by about a factor of ~ 2 as one goes to the centre of the galaxy or out to ~ 3 kpc. The minimum intrinsic HWHM of the HI vertical distribution in KK250 is ~ 350 pc. For comparison the HWHM of the vertical distribution of the HI in the solar neighbourhood is ~ 70-140 pc. Our results are hence consistent with other observations which indicate that dwarf galaxies have significantly puffier gas discs than spirals.

We report the timing and spatial resolution from the Muon Telescope Detector (MTD) installed in the STAR experiment at RHIC. Cosmic ray muons traversing the STAR detector have an average transverse momentum of 6 GeV/c. Due to their very small multiple scattering, these cosmic muons provide an ideal tool to calibrate the detectors and measure their timing and spatial resolution. The values obtained were ~100 ps and ~1-2 cm, respectively. These values are comparable to those obtained from cosmic-ray bench tests and test beams.

Using the Tremaine-Weinberg method, we measure the speed of the HI bar seen in the disk of NGC 3741. NGC 3741 is an extremely gas rich galaxy with an {H\,{\sc i}} disk which extends to about 8.3 times its Holmberg radius. It is also highly dark matter-dominated. Our calculated value of the pattern speed $\Omega_p$ is 17.1 $\pm$ 3.4 km $\textrm{s}^{-1}\textrm{kpc}^{-1}$. We also find the ratio of the co-rotation radius to the bar semi-major axis to be (1.6 $\pm$ 0.3), indicating a slow bar. This is consistent with bar models in which dynamical friction results in a slow bar in dark matter dominated galaxies.

Some low surface brightness galaxies are known to have extremely thin stellar disks with the vertical to planar axes ratio 0.1 or less, often referred to as superthin disks. Although their existence is known for over three decades, the physical origin for the thin distribution is not understood. We model the stellar thickness for a two-component (gravitationally coupled stars and gas) disk embedded in a dark matter halo, for a superthin galaxy UGC 7321 which has a dense, compact halo, and compare with a typical dwarf galaxy HoII which has a non-compact halo. We show that while the presence of gas does constrain the disk thickness, it is the compact dark matter halo which plays the decisive role in determining the superthin disk distribution in low-mass disks. Thus the compact dark matter halo significantly affects the disk structure and this could be important for the early evolution of galaxies.

Neutrino-neutrino interactions in dense neutrino streams, like those emitted by a core-collapse supernova, can lead to self-induced neutrino flavor conversions. While this is a nonlinear phenomenon, the onset of these conversions can be examined through a standard stability analysis of the linearized equations of motion. The problem is reduced to a linear eigenvalue equation that involves the neutrino density, energy spectrum, angular distribution, and matter density. In the single-angle case, we reproduce previous results and use them to identify two generic instabilities: The system is stable above a cutoff density ("cutoff mode"), or can approach an asymptotic instability for increasing density ("saturation mode"). We analyze multi-angle effects on these generic types of instabilities and find that even the saturation mode is suppressed at large densities. For both types of modes, a given multi-angle spectrum typically is unstable when the neutrino and electron densities are comparable, but stable when the neutrino density is much smaller or much larger than the electron density. The role of an instability in the SN context depends on the available growth time and on the range of affected modes. At large matter density, most modes are off-resonance even when the system is unstable.

A galactic disk in a spiral galaxy is generally believed to be embedded in an extended dark matter halo, which dominates its dynamics in the outer parts. However, the shape of the halo is not clearly understood. Here we show that the dark matter halo in the Milky Way Galaxy is prolate in shape. Further, it is increasingly more prolate at larger radii, with the vertical-to-planar axis ratio monotonically increasing to 2.0 at 24 kpc. This is obtained by modeling the observed steeply flaring atomic hydrogen gas layer in the outer Galactic disk, where the gas is supported by pressure against the net gravitational field of the disk and the halo. The resulting prolate-shaped halo can explain several long-standing puzzles in galactic dynamics, for example, it permits long-lived warps thus explaining their ubiquitous nature.

In this paper, we model dwarf galaxies as a two-component system of gravitationally coupled stars and atomic hydrogen gas in the external force field of a pseudo-isothermal dark matter halo, and numerically obtain the radial distribution of {H\,{\sc i}} vertical scale heights. This is done for a group of four dwarf galaxies (DDO\,154, Ho\,II, IC\,2574 and NGC\,2366) for which most necessary input parameters are available from observations. The formulation of the equations takes into account the rising rotation curves generally observed in dwarf galaxies. The inclusion of self-gravity of the gas into the model at par with that of the stars results in scale heights that are smaller than what was obtained by previous authors. This is important as the gas scale height is often used for deriving other physical quantities. The inclusion of gas self-gravity is particularly relevant in the case of dwarf galaxies where the gas cannot be considered a minor perturbation to the mass distribution of the stars. We find that three out of four galaxies studied show a flaring of their {H\,{\sc i}} disks with increasing radius, by a factor of a few within several disk scale lengths. The fourth galaxy has a thick {H\,{\sc i}} disk throughout. This arises as a result of the gas velocity dispersion remaining constant or decreasing only slightly while the disk mass distribution declines exponentially as a function of radius.

We model the shape and density profile of the dark matter halo of the low surface brightness, superthin galaxy UGC 7321, using the observed rotation curve and the HI scale height data as simultaneous constraints. We treat the galaxy as a gravitationally coupled system of stars and gas, responding to the gravitational potential of the dark matter halo. An isothermal halo of spherical shape with a core density in the range of 0.039 - 0.057 M_sun/pc^3 and a core radius between 2.5 - 2.9 kpc, gives the best fit to the observations for a range of realistic gas parameters assumed. We find that the best-fit core radius is only slightly higher than the stellar disc scale length (2.1 kpc), unlike the case of the high surface brightness galaxies where the halo core radius is typically 3-4 times the disc scale length of the stars. Thus our model shows that the dark matter halo dominates the dynamics of the low surface brightness, superthin galaxy UGC 7321 at all radii, including the inner parts of the galaxy.

In this paper, we use the outer-galactic HI scale height data as well as the observed rotation curve as constraints to determine the halo density distribution of the Andromeda galaxy (M31). We model the galaxy as a gravitationally-coupled system of stars and gas, responding to the external force-field of a known Hernquist bulge and the dark matter halo, the density profile of the latter being characterized by four free parameters. The parameter space of the halo is optimized so as to match the observed HI thickness distribution as well as the rotation curve on an equal footing, unlike the previous studies of M31 which were based on rotation curves alone. We show that an oblate halo, with an isothermal density profile, provides the best fit to the observed data. This gives a central density of 0.011 M_sun /pc^3, a core radius of 21 kpc, and an axis ratio of 0.4. The main result from this work is the flattened dark matter halo for M31, which is required to match the outer galactic HI scale height data. Interestingly, such flattened halos lie at the most oblate end of the distribution of halo shapes found in recent cosmological simulations.

Over the past two decades observations have revealed that the vertical density distribution of stars in galaxies near the mid-plane is substantially steeper than the sech^2 function that is expected from an isothermal approximation. However, the physical origin for this has not been explained so far. Here we show that such steep profiles result naturally even within the isothermal regime, on taking into account the gravitational force due to the gas. Due to its low velocity dispersion the gas is concentrated closer to the galactic mid-plane than the stars, and hence it strongly affects the vertical stellar distribution even though its contribution to the total surface density is small. We apply a three-component galactic disk model consisting of gravitationally coupled stars and the HI and H_2 gas, embedded in the dark matter halo, and calculate the vertical density distribution of stars for the Galaxy. The resulting vertical density distribution of stars is shown to be steeper than the sech^2 function, and lies between the sech and an exponential function, in good agreement with observations of galaxies. We also show that a multi-component stellar disk consisting of coupled dwarfs and two populations of giants does not explain the observed steep stellar profiles.