We extend the basic formalism of mimetic-metric-torsion gravity theory, in a way that the mimetic scalar field can manifest itself geometrically as the source of not only the trace mode of torsion, but also its axial (or, pseudo-trace) mode. Specifically, we consider the mimetic field to be (i) coupled explicitly to the well-known Holst extension of the Riemann-Cartan action, and (ii) identified with the square of the associated Barbero-Immirzi field, which is presumed to be a pseudo-scalar. The conformal symmetry originally prevalent in the theory would still hold, as the associated Cartan transformations do not affect the torsion pseudo-trace, and hence the Holst term. Demanding the theory to preserve the spatial parity symmetry as well, we focus on a geometric unification of the cosmological dark sector, and show that a super-accelerating regime in the course of evolution of the universe is always feasible. From the observational perspective, assuming the cosmological evolution profile to be very close to that for $\L$CDM, we further show that there could be a smooth crossing of the so-called phantom barrier at a low red-shift, however for a very restricted parametric domain. The extent of the super-acceleration have subsequently been ascertained by examining the evolution of the relevant torsion parameters.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

This work presents the broad-band time-averaged spectral analysis of neutron star low-mass X-ray binary, XTE J1710-281 by using the Suzaku archival data. The source was in a hard or an intermediate spectral state during this observation. This is the first time that a detailed spectral analysis of the persistent emission spectra of XTE J1710-281 has been done up to 30 keV with improved constraints on its spectral parameters. By simultaneously fitting the XIS (0.6-9.0 keV) and the HXD-PIN (15.0-30.0 keV) data, we have modelled the persistent spectrum of the source with models comprising a soft component from accretion disc and/or neutron star surface/boundary layer and a hard Comptonizing component. The 0.6-30 keV continuum with neutral absorber can be described by a multi-colour disc blackbody with an inner disc temperature of $kT_{\rm disc} = 0.28$ keV, which is significantly Comptonized by the hot electron cloud with electron temperature of $kT_{\rm e} \approx 5$ keV and described by photon index $\Gamma = 1.86$. A more complex three-component model comprising a multi-colour disc blackbody $\approx 0.30$ keV, single temperature blackbody $\approx 0.65$ keV and Comptonization from the disc, partially absorbed (about 38 per cent) by an ionized absorber (log($\xi$) $\approx$ 4) describes the broad-band spectrum equally well.

SAX J1748.9-2021 is a transient accretion powered millisecond X-ray pulsar located in the Globular cluster NGC 6440. We report on the spectral and timing analysis of SAX J1748.9-2021 performed on AstroSat data taken during its faint and short outburst of 2017. We derived the best-fitting orbital solution for the 2017 outburst and obtained an average local spin frequency of 442.361098(3) Hz. The pulse profile obtained from 3-7 keV and 7-20 keV energy bands suggest constant fractional amplitude ~0.5% for fundamental component, contrary to previously observed energy pulse profile dependence. Our AstroSat observations revealed the source to be in a hard spectral state. The 1-50 keV spectrum from SXT and LAXPC on-board AstroSat can be well described with a single temperature blackbody and thermal Comptonization. Moreover, we found that the combined spectra from XMM-Newton (EPIC-PN) and AstroSat (SXT+LAXPC) indicated the presence of reflection features in the form of iron (Fe K${\alpha}$) line that we modeled with the reflection model xillvercp. One of the two X-ray burst observed during the AstroSat/LAXPC observation showed hard X-ray emission (>30 keV) due to Compton up-scattering of thermal photons by the hot corona. Time resolved analysis performed on the bursts revealed complex evolution in emission radius of blackbody for second burst suggestive of mild photospheric radius expansion.

The Chevalier & Clegg (1985) spherical wind is widely used to model galactic outflows. Efforts are underway to understand the production of multiphase gas in such outflows. Since the important scales for the multiphase gas are very small compared to the size of the wind, it is very useful to have a local model of the cold gas interacting with a volume-filling expanding wind. In this note we explicitly derive the governing hydrodynamic equations in such a local box for the first time. We highlight that the expansion in a wind occurs only in the transverse direction, and this results in anisotropic source terms in the momentum equation. Our equations will be very useful in realistic modelling of the small scale structure of multiphase galactic outflows. [original note does not include an abstract]

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events.

We generalize the basic theory of mimetic gravity by extending its purview to the general metric-compatible geometries that admit torsion, in addition to curvature. This essentially implies reinstating the mimetic principle of isolating the conformal degree of freedom of gravity in presence of torsion, by parametrizing both the physical metric and torsion in terms of the scalar `mimetic' field and the metric and torsion of a fiducial space. We assert the requisite torsion parametrization from an inspection of the fiducial space Cartan transformation which, together with the conformal transformation of the fiducial metric, preserve the physical metric and torsion. In formulating the scalar-tensor equivalent Lagrangian, we consider an explicit contact coupling of the mimetic field with torsion, so that the former can manifest itself geometrically as the source of a torsion mode, and most importantly, give rise to a viable `dark universe' picture from a mimicry of an evolving dust-like cosmological fluid with a non-zero pressure. A further consideration of higher derivatives of the mimetic field in the Lagrangian leads to physical bounds on the mimetic-torsion coupling strength, which we determine explicitly.

The role of the symmetry energy on the internal constitution and the global structure of a cold nonaccreted neutron star is studied using a set of unified equations of state. Based on the nuclear energy-density functional theory, these equations of state provide a thermodynamically consistent treatment of all regions of the star and were calculated using the four different Brussels-Montreal functionals BSk22, BSk24, BSk25 and BSk26. Our predictions are compared to various constraints inferred from astrophysical observations including the recent detection of the gravitational wave signal GW170817 from a binary neutron-star merger.

The theory of the nuclear energy-density functional is used to provide a unified and thermodynamically consistent treatment of all regions of cold non-accreting neutron stars. In order to assess the impact of our lack of complete knowledge of the density dependence of the symmetry energy on the constitution and the global structure of neutron stars, we employ four different functionals. All of them were precision fitted to essentially all the nuclear-mass data with the Hartree-Fock-Bogoliubov method and two different neutron-matter equations of state based on realistic nuclear forces. For each functional, we calculate the composition, the pressure-density relation, and the chemical potentials throughout the star. We show that uncertainties in the symmetry energy can significantly affect the theoretical results for the composition and global structure of neutron stars. To facilitate astrophysical applications, we construct analytic fits to our numerical results.

We report the spectral analysis of accretion powered millisecond X-ray pulsar SAX J1748.9-2021 from the NuSTAR observation made during its 2015 outburst. The spectra displayed a broad emission line at $\sim 6.5$ keV with line width of $\sim 0.5$ keV and excess above $\sim 20$ keV due to Compton hump. The persistent emission of SAX J1748.9-2021 is described by a combination of soft thermal component with $kT =0.62^{+0.03}_{-0.05}$ keV and thermally Comptonized component with $kT_e=2.50^{+0.06}_{-0.03}$ keV reflected from the disc with reflection fraction of $0.30^{+0.08}_{-0.04}$. We have modeled the reflection spectrum with self-consistent model relxillCP and have found the inclination of the accretion disc to be $32.3^{^\circ +4.8}_{-4.7}$ and log $\xi$ is equal to $3.05^{+0.40}_{-0.34}$. We have obtained an upper limit of 57 km for inner disc radius; and the estimated magnetic field strength at the poles is less than $3.4 \times 10^9$ G. This kind of detailed investigation of reflection spectrum of SAX J1748.9-2021, especially above 10 keV, has been achieved only because of high statistics NuSTAR data of the source.

We report results from a broadband spectral analysis of the low mass X-ray binary MXB 1658-298 with the Swift/XRT and NuSTAR observations made during its 2015-2017 outburst. The source showed different spectral states and accretion rates during this outburst. The source was in low/hard state during 2015; and it was in high/soft state during the 2016 NuSTAR observations. This is the first time that a comparison of the soft and hard spectral states during an outburst is being reported in MXB 1658-298. Compared with the observation of 2015, the X-ray luminosity was about four times higher during 2016. The hard state spectrum can be well described with models consisting of a single-temperature blackbody component along with Comptonized disc emission, or three component model comprising of multi-colour disc, single-temperature blackbody and thermal Comptonization components. The soft state spectrum can be described with blackbody or disc blackbody and Comptonization component, where the neutron star surface (or boundary layer) is the dominant source for the Comptonization seed photons. We have also found a link between the spectral state of the source and the Fe K absorber. The absorption features due to highly ionized Fe were observed only in the soft state. This suggests a probable connection between accretion disc wind (or atmosphere) and spectral state (or accretion state) of MXB 1658-298.

In this paper, we will deform the usual energy-momentum dispersion relation of a photon gas at an intermediate scale between the Planck and electroweak scales. We will demonstrate that such a deformation can have non-trivial effects on the physics of dumb holes. So, motivated by the physics of dumb holes, we will first analyse the effect of such a deformation on thermodynamics. Then we observe that the velocity of sound also gets modified due to such a modification of the thermodynamics. This changes the position of the horizon of a dumb hole, and the analogous Hawking radiation from a dumb hole. Therefore, dumb holes can be used to measure the deformation of the usual energy-momentum dispersion relation.

We characterise the distribution of quasars within dark matter halos using a direct measurement technique for the first time at redshifts as high as $z \sim 1$. Using the Planck Sunyaev-Zeldovich (SZ) catalogue for galaxy groups and the Sloan Digital Sky Survey (SDSS) DR12 quasar dataset, we assign host clusters/groups to the quasars and make a measurement of the mean number of quasars within dark matter halos as a function of halo mass. We find that a simple power-law fit of $\log\left <N\right> = (2.11 \pm 0.01) \log (M) -(32.77 \pm 0.11)$ can be used to model the quasar fraction in dark matter halos. This suggests that the quasar fraction increases monotonically as a function of halo mass even to redshifts as high as $z\sim 1$.

We present an X-ray timing analysis of the transient X-ray binary MXB 1658-298, using data obtained from the RXTE and XMM-Newton observatories. We have made 27 new mid eclipse time measurements from observations made during the two outbursts of the source. These new measurements have been combined with the previously known values to study long term changes in orbital period of the binary system. We have found that the mid-eclipse timing record of MXB 1658-298 is quite unusual. The long term evolution of mid-eclipse times indicates an overall orbital period decay with a time scale of -- 6.5(7) x 10^7 year. Over and above this orbital period decay, the O-C residual curve also shows a periodic residual on shorter timescales. This sinusoidal variation has an amplitude of ~9 lt-sec and a period of ~760 d. This is indicative of presence of a third body around the compact X-ray binary. The mass and orbital radius of the third body are estimated to lie in the range, 20.5-26.9 Jupiter mass and 750-860 lt-sec, respectively. If true, then it will be the most massive circumbinary planet and also the smallest period binary known to host a planet.

4U 1822-37 is a Low Mass X-ray Binary (LMXB) system with an Accretion Disk Corona. We have obtained 16 new mid-eclipse time measurements of this source during the last 13 years using X- ray observations made with the RXTE-PCA, RXTE-ASM, Swift-XRT, XMM-Newton and Chandra observatories. These, along with the earlier known mid-eclipse times have been used to accurately determine the timescale for a change in the orbital period of 4U 1822-37. We have derived an orbital period Porb = 0.23210887(15) d, which is changing at the rate of \cdot Porb = 1.3(3) x 10-10 d d-1 (at T0 = MJD 45614). The timescale for a change in the orbital period is Porb/ \cdot Porb of 4.9(1.1) x 106 yr. We also report the detection of 0.59290132(11) s (at T0 = MJD 51975) X-ray pulsations from the source with a long term average \cdot Pspin of -2.481(4) x 10-12 s s-1, i.e., a spin-up time scale (Pspin/ \cdot Pspin) of 7578(13) yr. In view of these results, we have discussed various mechanisms that could be responsible for the orbital evolution in this system. Assuming the extreme case of conservative mass transfer, we have found that the measured \cdot Porb requires a large mass transfer rate of (4.2 - 5.2) x 10-8 M\odot yr-1 which together with the spin up rate implies a magnetic field strength in the range of (1-3) x 108 G. Using the long term RXTE-ASM light curve, we have found that the X-ray intensity of the source has decreased over the last 13 years by ? 40% and there are long term fluctuations at time scales of about a year. In addition to the long term intensity variation, we have also observed significant variation in the intensity during the eclipse.

4U 1626-67 is an accretion powered X-ray pulsar that shows remarkably stable X-ray luminosity above hours timescale and gradual intensity variation on a few years timescale. Here we report a significant increase in the X-ray intensity in the long term RXTE-ASM light curve of 4U 1626-67. Similar enhancement in the X-ray flux has also been detected in the Swift-BAT light curve. The increase in the X-ray flux took place over a long period of about 100 days and there appears to be two episodes of flux enhancement. We have investigated the spectral and timing features of 4U 1626-67 during its current state of enhanced flux emission with data obtained from the Proportional Counter Array and the High-Energy X-ray Timing Explorer on board the Rossi X-ray Timing Explorer. The source has entered a new spin-up phase with a spin-up rate of 4.02(5) $\times$ 10$^{-13}$ Hz s$^{-1}$. The present spin-up rate is almost half of the earlier spin-up and spin-down trends. A significant excess in soft X-ray photon emission is observed during the enhanced flux state, which is similar to the energy spectrum obtained during the spin-up era of the pulsar before 1990. We also report detection of a significant broadening in the wings of the 130 mHz peak and a change in the shape of the continuum of the power spectrum.

We present the orbital X-ray modulation study of three high mass X- ray binary systems, IGR J18027-2016, IGR J18483-0311 and IGR J16318-4848 using data obtained with RXTE-ASM, Swift-BAT and INTEGRAL-ISGRI. Using the long term light curves of the eclipsing HMXB IGR J18027-2016, obtained with Swift-BAT in the energy range 15-50 keV and INTEGRAL-ISGRI in the energy range 22-40 keV, we have determined three new mid eclipse times. The newly determined mid eclipse times together with the known values were used to derive an accurate value of the orbital period of 4.5693(4) d at MJD 52168 and an upper limit of 3.9(1.2)x10^-7 d d^-1 on the period derivative. We have also accurately determined an orbital period of 18.5482(88) d for the intermediate system IGR J18483-0311, which displays an unusual behaviour and shares many properties with the known SFXTs and persistent supergiant systems. This is a transient source and the outbursts occur intermittently at intervals of 18.55 d. Similarly, in the third supergiant system, IGR J16318-4848, we have found that the outbursts are separated by intervals of 80 d or its multiples, suggesting a possible orbital period.

We report here discovery of a 3.32 day orbital period in the Supergiant Fast X-ray Transient (SFXT) source IGR J16479-4514. Using the long term light curve of this source obtained with Swift-BAT in the energy range of 15-50 keV, we have clearly detected an orbital modulation including a full eclipse of duration ~0.6 day. In the hard X-ray band of the BAT instrument, the eclipse ingress and egress are rapid. We have also used the long term light curve obtained with the RXTE -ASM in the energy range of 1.5-12 keV. Taken independently, the detection of orbital modulation in the RXTE -ASM light curve is not significant. However, considering a clear detection of orbital modulation in the BAT light curve, we have used the ASM light curve for a more precise determination of the orbital period. IGR J16479-4514 has the shortest orbital period among the three SFXTs with measured/known orbital period. We discuss the implication of a short orbital period with the various mechanisms proposed to explain the transient nature of this class of sources.

An approach to the equation of state for the inner crust of neutron stars based on Skyrme-type forces is presented. Working within the Wigner-Seitz picture, the energy is calculated by the TETF (temperature-dependent extended Thomas-Fermi) method, with proton shell corrections added self-consistently by the Strutinsky-integral method. Using a Skyrme force that has been fitted to both neutron matter and to essentially all the nuclear mass data, we find strong proton shell effects: proton numbers $Z$ = 50, 40 and 20 are the only values possible in the inner crust, assuming that nuclear equilibrium is maintained in the cooling neutron star right down to the ambient temperature. Convergence problems with the TETF expansion for the entropy, and our way of handling them, are discussed. Full TETF expressions for the specific heat of inhomogeneous nuclear matter are presented. Our treatment of the electron gas, including its specific heat, is essentially exact, and is described in detail.