The 136 year long optical light curve of OJ~287 is explained by a binary black hole model where the secondary is in a 12 year orbit around the primary. Impacts of the secondary on the accretion disk of the primary generate a series of optical flares which follow a quasi-Keplerian relativistic mathematical model. The orientation of the binary in space is determined from the behavior of the primary jet. Here we ask how the jet of the secondary black hole projects onto the sky plane. Assuming that the jet is initially perpendicular to the disk, and that it is ballistic, we follow its evolution after the Lorentz transformation to the observer's frame. Since the orbital speed of the secondary is of the order of one-tenth of the speed of light, the result is a change in the jet direction by more than a radian during an orbital cycle. We match the theoretical jet line with the recent 12 $\mu$as-resolution RadioAstron map of OJ~287, and determine the only free parameter of the problem, the apparent speed of the jet relative to speed of light. It turns out that the Doppler factor of the jet, $\delta\sim5$, is much lower than in the primary jet. Besides following a unique shape of the jet path, the secondary jet is also distinguished by a different spectral shape than in the primary jet. The present result on the spectral shape agrees with the huge optical flare of 2021 November 12, also arising from the secondary jet.

We investigate the time evolution of sub-Keplerian transonic accretion flow onto a non-rotating black hole using axisymmetric viscous hydrodynamic simulations. We simulate the accretion flow using boundary values from semi-analytical analysis and set up three different models. Two of the models do not predict accretion shocks from the semi-analytic analysis, while one of them does. We also consider radiative cooling along with viscosity in the simulation. Our two-dimensional simulation deviated from the one-dimensional semi-analytical solution and admitted shocks in all three models. Viscous dissipation tends to push the shock front outward, and radiative cooling will push it in. Additionally, gravity is attractive. Depending on the competing strengths of all three processes, it may trigger shock oscillation. Different rates of angular-momentum transport in various layers may trigger eddies, which will enhance the shock oscillation. We show that any simple power law cannot approximate these solutions. We find that hot and higher angular-momentum flow requires higher viscosity to produce oscillatory shocks. From the temporal variation of the luminosity, shock oscillations generate QPOs in the range of sub-Hertz to a few Hertz frequencies if a ten solar mass black hole is assumed.

The Seyfert 1 galaxy J1626+5120 is estimated to host a $10^8 M_{\odot}$ black hole (BH) accreting at Eddington ratio $\dot{m}_{\text{Edd}} \approx 0.043$. Its long-term multi-band light curve data show flicker-like variations, but in a well-sampled $g$-band light curve, we are able to determine a $\simeq 329$\,d quasi-periodic oscillation (QPO) at a $\sim$4.53$\sigma$ significance. Six optical spectra were obtained for the source, three of which were taken by us. The spectra show that the variations were mainly because of flux changes blueward of 4000\,Å. We also analyze X-ray and ultraviolet (UV) data obtained with {\it the Neil Gehrels Swift Observatory (Swift)}, which targeted the source in the past two years. X-ray and UV emissions of the source show variations correlated with optical. Time lags of four UV bands and four optical bands are determined with respect to the X-ray emission, which are consistent with a continuum reprocessing disk model. These properties point out a disk origin for the QPO, likely due to Lense-Thirring (LT) precession of the accretion flow at $\sim$20 gravitational radii of the BH. This QPO could be a key case linking sub-year long QPOs in jets, which have more cases reported, to LT precession.

Radiative feedback from massive stars plays a central role in the evolution of molecular clouds and the interstellar medium. This paper presents a multi-wavelength analysis of the bright-rimmed cloud, BRC 44, which is located at the periphery of the Hii region Sh2-145 and is excited by the massive stars in the region. We use a combination of archival and newly obtained infrared data, along with new optical observations, to provide a census of young stellar objects (YSOs) in the region and to estimate stellar parameters such as age, mass etc. The spatial distribution of YSOs visible in the optical wavelength suggests that they are distributed in separate clumps compared to the embedded YSOs and are relatively older. Near-Infrared (NIR) spectroscopy of four YSOs in this region using the TANSPEC mounted on the 3.6m Devasthal Optical Telescope (DOT) confirms their youth. From Spectral Energy Distribution (SED) fitting, most of the embedded YSO candidates are in their early stage of evolution, with the majority of them in their Class II and some in Class I stage. The relative proper motions of the YSOs with respect to the ionizing source are indicative of the rocket effect in the BRC. The 12CO, 13CO, and C18O observations with the Purple Mountain Observatory are used to trace the distribution of molecular gas in the region. A comparison of the cold molecular gas distribution with simple analytical model calculations shows that the cloud is in the compression stage, and massive stars may be influencing the formation of young embedded stars in the BRC region due to radiative feedback.

At 66 Mpc, AT2019qiz is the closest optical tidal disruption event (TDE) to date, with a luminosity intermediate between the bulk of the population and iPTF16fnl. Its proximity allowed a very early detection and triggering of multiwavelength and spectroscopic follow-up well before maximum light. The velocity dispersion of the host galaxy and fits to the TDE light curve indicate a black hole mass $\approx 10^6$ M$_\odot$, disrupting a star of $\approx 1$ M$_\odot$. Comprehensive UV, optical and X-ray data shows that the early optical emission is dominated by an outflow, with a luminosity evolution $L \propto t^2$, consistent with a photosphere expanding at constant velocity ($\gtrsim 2000$ km s$^{-1}$), and a line-forming region producing initially blueshifted H and He II profiles with $v=3000-10000$ km s$^{-1}$. The fastest optical ejecta approach the velocity inferred from radio detections (modelled in a forthcoming companion paper from K.~D.~Alexander et al.), thus the same outflow may be responsible for both the fast optical rise and the radio emission -- the first time this connection has been observed in a TDE. The light curve rise begins $29 \pm 2$ days before maximum light, peaking when the photosphere reaches the radius where optical photons can escape. The photosphere then undergoes a sudden transition, first cooling at constant radius then contracting at constant temperature. At the same time, the blueshifts disappear from the spectrum and Bowen fluorescence lines (N III) become prominent, implying a source of far-UV photons, while the X-ray light curve peaks at $\approx 10^{41}$ erg s$^{-1}$. Assuming that these X-rays are from prompt accretion, the size and mass of the outflow are consistent with the reprocessing layer needed to explain the large optical to X-ray ratio in this and other optical TDEs, possibly favouring accretion-powered over collision-powered outflow models.

Star clusters, including high-mass stars, form within hub-filament systems (HFSs). Observations of HFSs that remain unaffected by feedback from embedded stars are rare yet crucial for understanding the mass inflow process in high-mass star formation. Using the JWST NIRCAM images, Dewangan et al. 2024, reported that the high-mass protostar G11P1 is embedded in a candidate HFS (G11P1-HFS; <0.6 pc). Utilizing ALMA N$_{2}$H$^{+}$(1-0) data, we confirm the presence of G11P1-HFS and study the dense gas kinematics. We analyzed the position-position-velocity (PPV) map and estimated on-sky velocity gradient ($V_g$) and gravity ($\mathcal{F}_{g}$) vectors. The spatial distribution of gas velocity and H$_2$ column density was examined. The steep $V_g$ of 5 km s$^{-1}$ pc$^{-1}$ and $-$7 km s$^{-1}$ pc$^{-1}$ toward either side of G11P1-hub, and the decreasing $V_g$ toward the hub, identify G11P1-HFS as a small-scale HFS in its nascent phase. $V_g$ and $\mathcal{F}_{g}$ align along the filaments, indicating gravity-driven flows. This work highlights the wiggled, funnel-shaped morphology of a HFS in PPV space, suggesting the importance of subfilaments or transverse gas flows in mass transportation to the hub.

PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required.

In this work, we employ two publicly available analysis tools to study four hydrogen(H)--stripped core--collapse supernovae (CCSNe) namely, SN 2009jf, iPTF13bvn, SN 2015ap, and SN 2016bau. We use the Modular Open-Source Fitter for Transients ({\tt MOSFiT}) to model the multi band light curves. {\tt MOSFiT} analyses show ejecta masses (log M$_{ej}$) of $0.80_{-0.13}^{+0.18}$ M$_{\odot}$, $0.15_{-0.09}^{+0.13}$ M$_{\odot}$, $0.19_{-0.03}^{+0.03}$ M$_{\odot}$, and $0.19_{+0.02}^{-0.01}$ M$_{\odot}$ for SN 2009jf, iPTF13vn, SN 2015ap, and SN 2016au, respectively. Later, Modules for Experiments in Stellar Astrophysics ({\tt MESA}), is used to construct models of stars from pre-main sequence upto core collapse which serve as the possible progenitors of these H-stripped CCSNe. Based on literature, we model a 12 M$_{\odot}$ ZAMS star as the possible progenitor for iPTF13vn, SN 2015ap, and SN 2016bau while a 20 M$_{\odot}$ ZAMS star is modeled as the possible progenitor for SN 2009jf. Glimpses of stellar engineering and the physical properties of models at various stages of their lifetime have been presented to demonstrate the usefulness of these analysis threads to understand the observed properties of several classes of transients in detail.

The highly variable blazar OJ~287 is commonly discussed as an example of a binary black hole system. The 130 year long optical light curve is well explained by a model where the central body is a massive black hole of 18.35$\times$10$^9$ solar mass that supports a thin accretion disc. The secondary black hole of 0.15$\times$10$^9$ solar mass impacts the disc twice during its 12 year orbit, and causes observable flares. Recently, it has been argued that an accretion disc with a typical AGN accretion rate and above mentioned central body mass should be at least six magnitudes brighter than OJ~287's host galaxy and would therefore be observationally excluded. Based on the observations of OJ~287's radio jet, detailed in Marscher and Jorstad (2011), and up-to-date accretion disc models of Azadi et al. (2022), we show that the V-band magnitude of the accretion disc is unlikely to exceed the host galaxy brightness by more than one magnitude, and could well be fainter than the host. This is because accretion power is necessary to launch the jet as well as to create electromagnetic radiation, distributed across many wavelengths, and not concentrated especially on the optical V-band. Further, we note that the claimed V-band concentration of accretion power leads to serious problems while interpreting observations of other Active Galactic Nuclei. Therefore, we infer that the mass of the primary black hole and its accretion rate do not need to be smaller than what is determined in the standard model for OJ~287.

Low Mass X-ray binaries (LMXBs) are binary systems where one of the components is either a black hole or a neutron star and the other is a less massive star. It is challenging to unambiguously determine whether a LMXB hosts a black hole or a neutron star. In the last few decades, multiple observational works have tried, with different levels of success, to address this problem. In this paper, we explore the use of machine learning to tackle this observational challenge. We train a random forest classifier to identify the type of compact object using the energy spectrum in the energy range 5-25 keV obtained from the Rossi X-ray Timing Explorer archive. We report an average accuracy of 87+/-13 in classifying the spectra of LMXB sources. We further use the trained model for predicting the classes for LMXB systems with unknown or ambiguous classification. With the ever-increasing volume of astronomical data in the X-ray domain from present and upcoming missions (e.g., SWIFT, XMM-Newton, XARM, ATHENA, NICER), such methods can be extremely useful for faster and robust classification of X-ray sources and can also be deployed as part of the data reduction pipeline.

Hot Jupiters should initially form at considerable distances from host stars and subsequently migrate towards inner regions, supported directly by transit timing variation (TTV). We report the TTV of K2-237b, using reproduced timings fitted from \textit{Kepler} K2 and \textit{TESS} data. The timings span from 2016 to 2021, leading to an observational baseline of 5 years. The timing evolution presents a significant bias to a constant period scenario. The model evidence is evaluated utilizing the Bayesian Information Criterion (BIC), which favours the scenario of period decay with a $\Delta$BIC of 14.1. The detected TTV induces a period decay rate ($\dot{P}$) of -1.14$\pm$0.28$\times$10$^{-8}$ days per day ($-$0.36 s/year). Fitting the spectral energy distribution, we find infrared excess at the significance level of 1.5 $\sigma$ for WISE W1 and W2 bands, and 2 $\sigma$ level for W3 and W4 bands. This potentially reveals the existence of a stellar disk, consisting of hot dust at 800$\pm$300 K, showing a $L_{dust}/L_{\ast}$ of 5$\pm$3$\times$10$^{-3}$. We obtain a stellar age of 1.0$^{+1.4}_{-0.7}$$\times$10$^{9}$ yr from isochrone fitting. The properties of K2-237b potentially serve as a direct observational support to the planet disk migration though more observation are needed.

The Galactic magnetar SGR J1935+2154 has undergone another outburst since 2022 October 10. We present the results of searching for an optical/NIR counterpart of SGR J1935+2154 before and during this outburst. No counterpart was detected at the magnetar's position in ${r'}$ and ${z'}$ bands, providing stringent upper limits of $r'\gtrsim 28.65$ and $z'\gtrsim 26.27$. Using archival X-ray data from NICER, we investigated the properties of the bursts and the spectral evolution of persistent emission. The burst flux $F$ showed a power-law distribution of $N\propto F^{-0.76\pm0.10}$ for flux $\gtrsim 2.6\times 10^{-9}\rm{\ erg\ cm^{-2}\ s^{-1}}$, while the temperature and radius followed a lognormal distribution with $kT=1.63^{+0.73}_{-0.50}\ \rm{keV}$ and $R_{\rm bb}=4.35_{-1.35}^{+1.95}\ \rm{km}$, respectively. The persistent flux evolution experienced a quick decay and an enhancement $\sim 27$ days after the BAT trigger. Using the near-infrared (NIR) and X-ray emission, together with the GTC optical/NIR upper limits, we discussed the origin of the NIR emission from the magnetar based on the fallback disk model and magnetosphere model. We found that either model cannot be ruled out with currently available data. Further mid-infrared observations are needed to find out the mechanism for producing the NIR emission from SGR J1935+2154.

The Vanishing & Appearing Sources during a Century of Observations (VASCO) project investigates astronomical surveys spanning a time interval of 70 years, searching for unusual and exotic transients. We present herein the VASCO Citizen Science Project, which can identify unusual candidates driven by three different approaches: hypothesis, exploratory, and machine learning, which is particularly useful for SETI searches. To address the big data challenge, VASCO combines three methods: the Virtual Observatory, user-aided machine learning, and visual inspection through citizen science. Here we demonstrate the citizen science project and its improved candidate selection process, and we give a progress report. We also present the VASCO citizen science network led by amateur astronomy associations mainly located in Algeria, Cameroon, and Nigeria. At the moment of writing, the citizen science project has carefully examined 15,593 candidate image pairs in the data (ca. 10% of the candidates), and has so far identified 798 objects classified as "vanished". The most interesting candidates will be followed up with optical and infrared imaging, together with the observations by the most potent radio telescopes.

In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ( |A| < 2\%) and precision (< 6\% per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

We analyze a prominence-like cool plasma structure as observed by Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). We perform the Differential Emission Measure (DEM) analysis using various filters of AIA, and also deduce the temperature and density structure in and around the observed flux-tube. In addition to deducing plasma parameters, we also find an evidence of multiple harmonics of fast magnetoacoustic kink waves in the observed prominence-like magnetic structure. Making use of estimated plasma parameters and observed wave parameters, under the baseline of MHD seismology, we deduce magnetic field in the flux-tube. The wave period ratio P1/P2 = 2.18 is also observed in the flux-tube, which carries the signature of magnetic field divergence where we estimate the tube expansion factor as 1.27. We discuss constraints in the estimation of plasma and magnetic field properties in such a structure in the current observational perspective, which may shed new light on the localized plasma dynamics and heating scenario in the solar atmosphere.

We present a detailed comparative systematic study using a sample of 221 Narrow-line Seyfert 1 (NLSy1) galaxies in comparison to a redshift matched sample of 154 Broad-line Seyfert 1 (BLSy1) galaxies based on their observations using ROSAT and/or XMM-Newton telescopes in soft X-ray band (0.1-2.0 keV). A homogeneous analysis is carried out to estimate their soft X-ray photon indices ($\Gamma^{s}_{X}$) and its correlations with other parameters of nuclear activities such as Eddington ratios (R$_\mathrm{Edd}$), bolometric luminosities (L$_\mathrm{bol}$), black hole masses (M$_\mathrm{BH}$) and the widths of the broad component of H$\beta$ lines (FWHM(H$\beta$)). In our analysis, we found clear evidence of the difference in the $\Gamma^{s}_{X}$ and R$_\mathrm{Edd}$ distributions among NLSy1 and BLSy1 galaxies, with steeper $\Gamma^{s}_{X}$ and higher R$_\mathrm{Edd}$ for the former. Such a difference also exists in the spectral indices distribution in hard X-ray ($\Gamma^{h}_{X}$), based on the analysis of 53 NLSy1 and 46 BLSy1 galaxies in the 2-10 keV energy band. The difference in R$_\mathrm{Edd}$ distributions does exist even after applying the average correction for the difference in the inclination angle of NLSy1 and BLSy1 galaxies. We also estimated R$_\mathrm{Edd}$, based on SED fitting of 34 NLSy1 and 30 BLSy1 galaxies over the 0.3-10 keV energy band and found that results are still consistent with R$_\mathrm{Edd}$ estimates based on the optical bolometric luminosity. Our analysis suggests that the higher R$_\mathrm{Edd}$ in NLSy1 is responsible for its steeper X-ray spectral slope compared to the BLSy1, consistent with the disc-corona model as proposed for the luminous AGNs.

The Devasthal Fast Optical Telescope (DFOT) is a 1.3 meter aperture optical telescope, recently installed at Devasthal, Nainital. We present here the first results using an \Ha filter with this telescope on a Wolf-Rayet dwarf galaxy Mrk 996. The instrumental response and the \Ha sensitivity obtained with the telescope are $(3.3 \pm 0.3)\times 10^{-15} (erg s^{-1} cm^{-2})/(counts s^{-1})$and$7.5\times10^{-17} erg s^{-1} cm^{-2} arcsec^{-2}$ respectively. The \Ha flux and the equivalent width for Mrk~996 are estimated as $(132 \pm 37)\times 10^{-14} erg s^{-1} cm^{-2}$and$\sim$96$\mathring{A}$ respectively. The star formation rate is estimated as $0.4\pm0.1 M_\odot$yr$^{-1}$. Mrk~996 deviates from the radio-FIR correlation known for normal star forming galaxies with a deficiency in its radio continuum. The ionized gas as traced by \Ha emission is found in a disk shape which is misaligned with respect to the old stellar disk. This misalignment is indicative of a recent tidal interaction in the galaxy. We believe that galaxy-galaxy tidal interaction is the main cause of the WR phase in Mrk~996.

A very unique strength of the Devasthal Observatory is its capability of detecting optical transients with the 4-m International Liquid Mirror Telescope (ILMT) and to rapidly follow them up using the 1.3-m Devasthal Fast Optical Telescope (DFOT) and/or the 3.6-m Devasthal Optical Telescope (DOT), installed right next to it. In this context, we have inspected 20 fields observed during 9 consecutive nights in October-November 2022 during the first commissioning phase of the ILMT. Each of these fields has an angular extent of $22^\prime$ in declination by $9 \times 22^\prime$ in right ascension. Combining both a visual search for optical transients and an automatic search for these using an image subtraction technique (see the ILMT poster paper by Pranshu et al.), we report a total of 232 significant transient candidates. After consulting the Minor Planet Center database of asteroids, we could identify among these 219 positions of known asteroids brighter than $V=22$. These correspond to the confirmed positions of 78 distinct known asteroids. Analysis of the remaining CCD frames covering 19 more fields (out of 20) should lead to an impressive number of asteroids observed in only 9 nights. The conclusion is that in order to detect and characterize new supernovae, micro-lensing events, highly variable stars, multiply imaged quasars, etc. among the ILMT optical transients, we shall first have to identify all known and new asteroids. Thanks to its large diameter and short focal length (f/D $\sim$ 2.4), the ILMT turns out to be an excellent asteroid hunter.

We present B-, V-, R-, and I-bands linear polarimetric observations of 73 stars in the direction of open star cluster Casado Alessi 1 (hereafter Alessi 1). We aim to use polarimetry as a tool to investigate the properties and distribution of dust grains toward the direction of the cluster. The polarimetric observations were carried out using the ARIES IMaging POLarimeter mounted at the 104 cm telescope of ARIES, Nainital (India). Using the Gaia photometric data the age and distance of the cluster are estimated to be $0.8\pm0.1$ Gyr and $673\pm98$ pc, respectively. A total of 66 stars with a 26 arcmin radius from the cluster are identified as members of the cluster using the astrometric approach. Out of these 66 members, 15 stars were observed polarimetrically and found to have the same value of polarization. The majority of the stars in the region follow the general law of the polarization for the interstellar medium, indicating that polarization toward the cluster Alessi 1 is dominated by foreground dust grains. The average values of the maximum polarization ($P_{max}$) and the wavelength corresponding to the maximum polarization ($\lambda_{max}$) toward the cluster are found to be $0.83\pm0.03$% and $0.59\pm0.04$ $\mu$m, respectively. Also, dust grains toward the cluster appear to be aligned, possibly due to the galactic magnetic field.

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.