We report the results of the transit timing variation (TTV) analysis of the extra-solar planet Qatar-1b using thirty eight light curves. Our analysis combines thirty five previously available transit light curves with three new transits observed by us between June 2016 and September 2016 using the 2-m Himalayan Chandra Telescope (HCT) at the Indian Astronomical Observatory (Hanle, India). From these transit data, the physical and orbital parameters of the Qatar-1 system are determined. In addition to this, the ephemeris for the orbital period and mid-transit time are refined to investigate the possible TTV. We find that the null-TTV model provides the better fit to the (O-C) data. This indicates that there is no evidence for TTVs to confirm the presence of additional planets in the Qatar-1 system. The use of the 3.6-m Devasthal Optical Telescope (DOT) operated by the Aryabhatta Research Institute of Observational Sciences (ARIES, Nainital, India) could improve the photometric precision to examine the signature of TTVs in this system with a greater accuracy than in the present work.

We present twelve new transit light curves of the hot-Jupiter TrES-3b observed during $2012-2018$ to probe the transit timing variation (TTV). By combining the mid-transit times determined from these twelve transit data with those re-estimated through uniform procedure from seventy one transit data available in the literature, we derive new linear ephemeris and obtain the timing residuals that suggest the possibility of TTV in TrES-3 system. However, the frequency analysis shows that the possible TTV is unlikely to be periodic, indicating the absence of an additional body in this system. To explore the other possible origins of TTV, the orbital decay and apsidal precession ephemeris models are fitted to the transit time data. We find decay rate of TrES-3b to be $\bf \dot{P_q}= -4.1 \pm 3.1$ $ms$ ${yr}^{-1}$ and the corresponding estimated modified tidal quality factor of ${Q}^{'}_{\ast}$ $\sim 1.11 \times {10}^{5}$ is consistent with the theoretically predicted values for the stars hosting the hot-Jupiters. The shift in the transit arrival time of TrES-3b after eleven years is expected to be ${T}_{shift}\sim 69.55 \ s$, which is consistent with the $RMS$ of the timing residuals. Besides, we find that the apsidal precession ephemeris model is statistically less probable than the other considered ephemeris models. It is also discussed that despite the linear ephemeris model appears to be the most plausible model to represent the transit time data, the possibility of the orbital decay cannot be completely ruled out in TrES-3 system. In order to confirm this, further high-precision and high-cadence follow-up observation of transits of TrES-3b would be important.

We report the detection of type C QPO along with the upper harmonic in the commensurate ratio of 1:2 in the two observations of the low-mass black hole transient H~1743--322 jointly observed by \textit{XMM-Newton} and \textit{NuSTAR} during the 2016 outburst. We find that the QPO and the upper harmonic exhibit shifts in their centroid frequencies in the second observation with respect to the first one. The hardness intensity diagram implies that in contrast to 2008 and 2014 failed outbursts, 2016 outburst was a successful one. We also detect the presence of a broad iron K$\alpha$ line at $\sim$6.5 keV and reflection hump in the energy range of 15--30 keV in both the observations. Along with the shape of the power density spectra, the nature of the characteristic frequencies and the fractional rms amplitude of the timing features imply that the source stayed in the low/hard state during these observations. Moreover, photon index and other spectral parameters also indicate the low/hard state behavior of the source. Unlike the soft lag detected in this source during the 2008 and 2014 failed outbursts, we observe hard time-lag of $0.40\pm0.15$ and $0.32\pm0.07$ s in the 0.07--0.4 Hz frequency range in the two observations during the 2016 outburst. The correlation between the photon index and the centroid frequency of the QPO is consistent with the previous results. Furthermore, the high value of the Comptonized fraction and the weak thermal component indicate that the QPO is being modulated by the Comptonization process.

We study cosmological models with modified Chaplygin gas (in short, MCG) to determine observational constraints on its EoS parameters. The observational data of the background and the growth tests are employed. The background test data namely, H(z)-z data, CMB shift parameter, Baryonic acoustic oscillations (BAO) peak parameter, SN Ia data are considered to study the dynamical aspects of the universe. The growth test data we employ here consists of the linear growth function for the large scale structures of the universe, models are explored assuming MCG as a candidate for dark energy. Considering the observational growth data for a given range of redshift from the Wiggle-Z measurements and rms mass fluctuations from Ly-$\alpha$ measurements, cosmological models are analyzed numerically to determine constraints on the MCG parameters. In this case, the Wang-Steinhardt ansatz for the growth index $\gamma$ and growth function $f$ (defined as $f=\Omega_{m}^{\gamma} (a)$) are also taken into account for the numerical analysis. The best-fit values of the equation of state parameters obtained here are used to study the variation of the growth function ($f$), growth index ($\gamma$), equation of state parameter ($\omega$), squared sound speed $c^2_{s}$ and deceleration parameter with redshift. The observational constraints on the MCG parameters obtained here are then compared with those of the GCG model for viable cosmology. It is noted that MCG models satisfactorily accommodate an accelerating phase followed by a matter dominated phase of the universe. The permitted range of values of the EoS parameters and the associated parameters ($f$, $\gamma$, $\omega$, $\Omega$, $c^2_{s}$, $q$) are compared with those obtained earlier using other observations.

We investigate emergent universe model using recent observational data for the background tests as well as for the growth tests. The background test data comprises of Hubble data, Baryon Acoustic Oscillation (BAO) data, cosmic microwave background (CMB) shift data, Union compilation data. The observational growth data are obtained from Wiggle-Z measurements and rms mass fluctuations data from Lyman-$\alpha$ measurements at different red shifts. The flat emergent universe model obtained by Mukherjee {\it et. al.} is permitted with a non-linear equation of state (in short, EoS) ($p=A\rho-B\rho^{\frac{1}{2}}$), where $A$ and $B$ are constants. The observed cosmological data are used here to estimate the range of allowed values of EoS parameters numerically. The best-fit values of the EoS and growth parameters are determined making use of chi-square minimization technique. The analysis is carried out here considering the Wang-Steinhardt ansatz for growth index ($\gamma$) and growth function ($f$ defined as $f=\Omega_{m}^{\gamma} (a)$). Subsequently, the best-fit values of the EoS parameters are used to study the evolution of the growth function $f$, growth index $\gamma$, state parameter $\omega$ with red shift parameter $z$. The present value of the parameters, namely, $f$, $\gamma$, $\omega$, $\Omega_{m}$ are also estimated. The late accelerating phase of the universe in the model is accommodated satisfactorily.

Five newly observed transit light curves of the TrES-3 planetary system are presented. Together with other light curve data from literature, 23 transit light curves in total, which cover an overall timescale of 911 epochs, have been analyzed through a standard procedure. From these observational data, the system's orbital parameters are determined and possible transit timing variations are investigated. Given that a null transit-timing-variation produces a fit with reduced chi^2=1.52, our results agree with previous work, that transit timing variations might not exist in these data. However, a 1-frequency oscillating transit-timing-variation model, giving a fit with a reduced chi^2=0.93, does possess a statistically higher probability. It is, thus, concluded that future observations and dynamical simulations for this planetary system will be very important.

Through the international collaborators, we recently established a network of existing and working meter-class telescopes to look for planetary transit events. As a first step, we focus on the TrES3 system, and conclude that there could be some level of transit timing variations.

We investigate the linear growth function for the large scale structures of the universe considering modified Chaplygin gas as dark energy. Taking into account observational growth data for a given range of redshift from the Wiggle-Z measurements and rms mass fluctuations from Ly-$\alpha$ measurements we numerically analyze cosmological models to constrain the parameters of the MCG. The observational data of Hubble parameter with redshift ($Z$) is also considered. The Wang-Steinhardt ansatz for growth index $\gamma$ and growth function $f$ (defined as $f=\Omega_{m}^{\gamma} (a)$) are considered for the numerical analysis. The best-fit values of the equation of state parameters obtained here is employed to study the growth function $(f)$, growth index ($\gamma$) and state parameter ($\omega$) with redshift $z$. We note that MCG satisfactorily accommodates an accelerating phase followed by a matter dominated universe.

We present Cosmological models with modified Chaplygin gas (MCG) in the framework of Horava-Lifshitz (HL) theory of gravity both with and without detailed balance. The equation of state (EOS) for a MCG contains three unknown parameters namely, $A$, $\alpha$, $B$. The allowed values of some of these parameters of the EOS are determined using the recent astrophysical and cosmological observational data. Using observational data from $H(z)-z$, BAO peak parameter, CMB shift parameter we study cosmologies in detailed-balance and beyond detailed-balance scenario. In this paper we take up the beyond detailed-balance scenario in totality and contribution of dark radiation in the case of detailed-balance scenario on the parameters of the EOS. We explore the effect of dark radiation on the whole range the of effective neutrino parameter to constrain matter contributing parameter $B$ in both the detailed-balance and the beyond-detailed balance scenario. It has been observed that greater the dark radiation less the matter contribution in the MCG in both the scenario considered here. In order to check the validity of beyond detailed balance scenario we plot supernovae magnitudes ($\mu$) with redshift of Union2 data and then the variation of state parameter with redshift is studied. It has been observed that beyond detailed balance scenario is equally suitable in HL gravity with MCG.

A class of Emergent Universe (EU) model is studied in the light of recent observational data. Significant constraints on model parameters are obtained from the observational data. Density parameter for a class of model is evaluated. Some of the models are in favour of the recent observations. Some models have been found which are not interesting yielding unrealistic present day value of the density parameter.

A specific class of flat Emergent Universe (EU) is considered and its viability is tested in view of the recent observations. Model parameters are constrained from Stern data for Hubble Parameter and Redshift ($H(z)$ vs. $z$) and from a model independent measurement of BAO peak parameter. It is noted that a composition of Exotic matter, dust and dark energy, capable of producing an EU, can not be ruled out with present data. Evolution of other relevant cosmological parameters, viz. density parameter ($\Omega$), effective equation of state (EOS) parameter ($\omega_{eff}$) are also shown.

We study a composition of normal and exotic matter which is required for a flat Emergent Universe scenario permitted by the equation of state (EOS)($p=A\rho-B\rho^{1/2}$) and predict the range of the permissible values for the parameters $A$ and $B$ to explore a physically viable cosmological model. The permitted values of the parameters are determined taking into account the $H(z)-z$ data obtained from observations, a model independent BAO peak parameter and CMB shift parameter (WMAP7 data). It is found that although $A$ can be very close to zero, most of the observations favours a small and negative $A$. As a consequence, the effective Equation of State parameter for this class of Emergent Universe solutions remains negative always. We also compared the magnitude ($\mu (z)$) vs. redshift($z$) curve obtained in the model with that obtained from the union compilation data. According to our analysis the class of Emergent Universe solutions considered here is not ruled out by the observations.

C. Ragoonatha Chary, the first assistant at Madras Observatory during 1864 to 1880 was not only a celebrated observational astronomer but also a person who emphasized the need for incorporating modern observations based improvements into the traditional methods of astronomical calculations. He was one of the first few people who argued for establishment of independent modern Indian observatory for education and training. He was credited with the discovery of two variable stars R Reticuli and another star whose identity is uncertain. The person and his variable star discoveries are discussed.

We study Modified Chaplygin Gas (MCG) as a candidate for dark energy and predict the values of parameters of the gas for a physically viable cosmological model. The equation of state of MCG ($p=B \rho - \frac {A}{\rho^\alpha} $) involves three parameters: $B$, $A$ and $\alpha$. The permitted values of these parameters are determined with the help of dimensionless age parameter ($H_{o}t_{o}$) and $H(z)-z$ Data. Specifically we study the allowed ranges of values of B parameter in terms of $\alpha$ and $A_{s}$ ($A_{s}$ is defined in terms of the constants in the theory). We explore the constraints of the parameters in Cold Dark Matter(CDM) model and UDME(Unified Dark Matter Energy) model respectively.

We have performed smoothed particle hydrodynamics (SPH) simulations to study the response of the central kiloparsec region of a gaseous disk to the imposition of nonaxisymmetric bar potentials. The model galaxies are composed of the three axisymmetric components (halo, disk, and bulge) and a non-axisymmetric bar. These components are assumed to be invariant in time in the frame corotating with the bar. The potential of spherical $\gamma$-models of Dehnen is adopted for the bulge component whose density varies as $r^{-\gamma}$ near the center and $r^{-4}$ at larger radiiand hence, possesses a central density core for $\gamma = 0$ and cusps for $\gamma > 0$. Since the central mass concentration of the model galaxies increases with the cusp parameter $\gamma$, we have examined here the effect of the central mass concentration by varying the cusp parameter $\gamma$ on the mechanism responsible for the formation of the symmetric two-armed nuclear spirals in barred galaxies. Our simulations show that the symmetric two-armed nuclear spirals are formed by hydrodynamic spiral shocks driven by the gravitational torque of the bar for the models with $\gamma = 0$ and 0.5. On the other hand, the symmetric two-armed nuclear spirals in the models with $\gamma=1$ and 1.5 are explained by gas density waves. Thus, we conclude that the mechanism responsible for the formation of the symmetric two-armed nuclear spirals in barred galaxies changes from the hydrodynamic shocks to the gas density waves when the central mass concentration increases from $\gamma = 0$ to 1.5.

We use the usual method of Schwarzschild to construct self-consistent solutions for the triaxial de Zeeuw & Carollo (1996) models with central density cusps. ZC96 models are triaxial generalisations of spherical $\gamma$-models of Dehnen whose densities vary as $r^{-\gamma}$ near the center and $r^{-4}$ at large radii and hence, possess a central density core for $\gamma=0$ and cusps for $\gamma > 0$. We consider four triaxial models from ZC96, two prolate triaxials: $(p, q) = (0.65, 0.60)$ with $\gamma = 1.0$ and 1.5, and two oblate triaxials: $(p, q) = (0.95, 0.60)$ with $\gamma = 1.0$ and 1.5. We compute 4500 orbits in each model for time periods of $10^{5} T_{D}$. We find that a large fraction of the orbits in each model are stochastic by means of their nonzero Liapunov exponents. The stochastic orbits in each model can sustain regular shapes for $\sim 10^{3} T_{D}$ or longer, which suggests that they diffuse slowly through their allowed phase-space. Except for the oblate triaxial models with $\gamma =1.0$, our attempts to construct self-consistent solutions employing only the regular orbits fail for the remaining three models. However, the self-consistent solutions are found to exist for all models when the stochastic and regular orbits are treated in the same way because the mixing-time, $\sim10^{4} T_{D}$, is shorter than the integration time, $10^{5} T_{D}$. Moreover, the ``fully-mixed'' solutions can also be constructed for all models when the stochastic orbits are fully mixed at 15 lowest energy shells. Thus, we conclude that the self-consistent solutions exist for our selected prolate and oblate triaxial models with $\gamma = 1.0$ and 1.5.

The projected properties of triaxial generalization of the modified Hubble mass models are studied. These models are constructed by adding the additional radial functions, each multiplied by a low-order spherical harmonic, to the models of \citet{ct00}. The projected surface density of mass models can be calculated analytically which allows us to derive the analytic expressions of axial ratio and position angles of major axis of constant density elliptical contours at asymptotic radii. The models are more general than those studied earlier in the sense that the inclusions of additional terms in density distribution, allows one to produce varieties of the radial profile of axial ratio and position angle, in particular, their small scale variations at inner radii. Strong correlations are found to exist between the observed axial ratio evaluated at $0.25 R_{e}$ and at $4 R_{e}$ which occupy well-separated regions in the parameter space for different choices of the intrinsic axial ratios. These correlations can be exploited to predict the intrinsic shape of the mass model, independent of the viewing angles. Using Bayesian statistics, the result of a test case launched for an estimation of the shape of a model galaxy is found to be satisfactory.

We have performed smoothed particle hydrodynamics (SPH) simulations for the response of the gaseous disk to the imposed moderately strong non-axisymmetric potentials. The model galaxies are composed of the three stellar components (disk, bulge and bar) and two dark ones (supermassive black hole and halo) whose gravitational potentials are assumed to be invariant in time in the frame corotating with the bar. We found that the torques alone generated by the moderately strong bar that gives the maximum of tangential-to-radial force ratio as $(F_{Tan}/F_{Rad})_{max}= 0.3$ are not sufficient to drive the gas particles close to the center due to the barrier imposed by the inner Lindblad resonances (ILRs). In order to transport the gas particles towards the nucleus ($r<100$ pc), a central supermassive black hole (SMBH) and high sound speed of the gas are required to be present. The former is required to remove the inner inner Lindblad resonance (IILR) that prevents gas inflow close to the nucleus, while the latter provides favourable conditions for the gas particles to lose their angular momentum and to spiral in. Our models that have no IILR show the trailing nuclear spirals whose innermost parts reach close to the center in a curling way when the gas sound speed is $ c_{s} \gtrsim 15$ km s$^{-1}$. They resemble the symmetric two-armed nuclear spirals observed in the central kiloparsec of spiral galaxies. We found that the symmetric two-armed nuclear spirals are formed by the hydrodynamic spiral shocks caused by the gravitational torque of the bar in the presence of a central SMBH that can remove IILR when the sound speed of gas is high enough to drive a large amount of gas inflow deep inside the ILR. However, the detailed morphology of nuclear spirals depends on the sound speed of gas.

We have performed SPH simulations for the response of the gaseous disks to the imposed potentials including those from bars and SMBHs. Evolution of the nuclear regions of gaseous disks depends critically on the masses of SMBHs as well as the sound speeds in the gas.