The gravitational lens SDSS J1004+4112 was the first discovered system where a background quasar is lensed by a galaxy cluster instead of a single galaxy. We use the 14.5-year r-band light curves together with the recently measured time delay of the fourth brightest quasar image (Munõz et al. (2022)) and the mass model from Forés-Toribio et al. (2022) to study the microlensing effect in this system. We constrain the quasar accretion disk size to light-days at 2407Å in the restframe which is compatible with most previous estimates. We also infer the fraction of mass in stars at the positions of the quasar images: $${\alpha _A} = 0.058_{ - 0.032}^{ + 0.024},{\alpha _B} = 0.048_{ - 0.014}^{ + 0.032},{\alpha _C} = 0.018_{ - 0.018}^{ + 0.015}$$ and $${\alpha _D} = 0.008_{ - 0.008}^{ + 0.033}$$. The stellar fraction estimates are reasonable for intracluster medium although the stellar fractions at images A and B are slightly larger, suggesting the presence of a near undetected galaxy.

Recent ALMA observations detected protostellar outflows in 70-μm dark infrared dark clouds (IRDCs). These sources are candidates for the initial stages of high-mass star formation. We launched a new survey for free-free emission from outflow shocks using the Yamaguchi Interferometer (YI) at 8 GHz. We aim to catalog “proto-high-mass protostar” candidates that are still in the low to intermediate-mass phase. We selected starless-like clumps without any 70-μm point source from Traficante et al. (2015). We currently detected 82 sources from 167 clumps. 37 of them are fainter than 20 mJy (down to a few mJy). They tend to associate with colder and denser clumps that are suitable for star formation. This fact suggests that, at least, some of them trace star-formation activities. The highest-density clumps are, in fact, associated with several masers and molecular outflows. Furthermore, some of them have already shown a signature of ongoing cluster formation.

The 6.7 GHz methanol maser emitted by a high-mass star forming region G33.641-0.228 is known to exhibit fast flux density variability on timescales of less than one day. The mechanism of this variability, called burst, has not been known. We observed the circular polarization of this maser. As a result, we found that only the spectral components representing burst exhibit strong circular polarization exceeding 10%. This suggests that the two phenomena of the burst and circular polarization are related to each other.

We report on a direct comparison of VLBI maser data and ALMA thermal-emission data for the high-mass protostar G353.273+0.641. We detected a gravitationally-unstable disk by dust and a high-velocity jet traced by a thermal CO line by ALMA long-baselines (LB). 6.7 GHz CH3OH masers trace infalling streamlines inside the disk. The innermost maser ring indicates another compact accretion disk of 30 au. Such a nested system could be caused by angular momentum transfer by the spiral arms. 22 GHz H2O masers trace the jet-accelerating region, which are directly connecting the CO jet and the protostar. The recurrent maser flares imply episodic jet ejections per 1–2 yr, while typical separation of CO knots indicates a variation of outflow rate per 100 yr. Our study demonstrates that VLBI maser observations are still a powerful tool to explore detailed structures nearby high-mass protostars by combining ALMA LB.

In this contribution, we explore the question on the formation of multiple massive stellar systems via disk fragmentation with the help of the highest-resolution simulations to date of a fragmenting disk in the context of massive star formation. The simulations start from a collapsing cloud of 200 solar masses, followed by the formation of an accretion disk that develops spiral arms and fragments. Due to the high resolution of our grid, we are able to self-consistently form the fragments without the need for a subgrid module such as sink particles. We track the formed fragments into the first stages of companion formation, which allows us to give an estimate of the multiplicity of the final system due to disk fragmentation. We find in total around ∼6 fragments, some at orbits of ∼ 1000 au, and some close (possibly spectroscopic) companions.

We model evolution of exoplanets of S-type in close binary systems at the stage when the companion starts to lose mass via a slow stellar wind. At this stage an accretion disc is formed around the planets’ host. Detailed structure of such discs is calculated in quasi-stationary and non-stationary approaches. We model migration of planets embedded in these discs.

Situation with highly magnetized neutron stars in binary systems is not yet certain. On the one hand, all best studied magnetars seem to be isolated objects. On the other, there are many claims based on model-dependent analysis of spin properties or/and luminosity of neutron stars in X-ray binaries in favour of large fields. In addition, there are a few results suggesting a magnetar-like activity of neutron stars in close binary systems. Most of theoretical considerations do not favour even existence, not speaking about active decay, of magnetar-scale fields in neutron stars older than ∼106 yrs. However, alternative scenarios of the field evolution exist. I provide a brief review of theoretical and observational results related to the presence of neutron stars with large magnetic field in binaries and discuss perspectives of future studies.

Many simulations have been performed to elucidate the formation process of first stars. In first star formation, radiative feedback is a key process in determining stellar masses. However, previous simulations which follow the feedback process don’t resolve the small scale ( 10 AU) to realize long-term calculation, and the structure near massive protostars is still unknown. To clarify how the radiation from the protostar works, we need to resolve small scale and calculate the interaction between the radiation and the dense gas in such a region. As a first step towards understanding the phenomenon in this region, we perform the high-resolution simulation around the massive protostar without radiative transfer. We find that dense gas covers the protostar even in the polar direction and the HII region cannot expand. Solving the radiative transfer for getting accurate results is our future work. We are currently developing the new radiation hydrodynamics code for that.

We introduce the HandOff, a set of computational tools to transfer general relativistic magnetohydrodynamic and spacetime data from the numerical code IllinoisGRMHD to Harm3d. While the former simulates binary neutron star (BNS) mergers in full dynamical General Relativity, the latter specializes in modeling accretion disks around black holes (BH) over long timescales. The HandOff, then, enables us to transfer the outcome of BNS mergers after BH formation from IllinoisGRMHD to Harm3d and to evolve the post-merger system efficiently over long timescales. We show our latest results in this respect, for matter modeled as a magnetized ideal-fluid, and discuss our future plans which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.

We analyzed 39 ks NuSTAR data of Cen X-3 through both orbital- and pulse-phase resolved spectroscopy. Orbital-phase resolved spectra show extrinsic fluctuations due to absorption by surrounding plasma, as the spectral fluctuation mainly emerges below 10 keV. Pulse-phase resolved spectra, on the other hand, show intrinsic fluctuations depending on effectiveness of Comptonization, since the spectrum becomes hard above 10 keV at the pulse peak.

One of the models explaining the high luminosity of pulsing ultra-luminous X-ray sources (pULXs) was suggested by Mushtukov et al. (2015). They showed that the accretion columns on the surfaces of highly magnetized neutron stars can be very luminous due to opacity reduction in the high magnetic field. However, a strong magnetic field leads also to amplification of the electron-positron pairs creation. Therefore, increasing of the electron and positron number densities compensates the cross-section reduction, and the electron scattering opacity does not decrease with the magnetic field magnification. As a result, the maximum possible luminosity of the accretion column does not increase with the magnetic field. It ranges between 1040 − 1041 erg s−1 depending only slightly on the magnetic field strength.

Hα emission is one of the most prominent features of young stellar objects in the optical range, and importantly, the equivalent width (EW) of Hα emission [EW(Hα)] is used to characterise an evolutionary stage of young stars. The aim of this work is to identify and study the stellar objects with variable EW(Hα) in the young stellar cluster IC 348. We performed photometric and slit-less observations at several epochs in order to reveal the variable objects. Significant variability of EW(Hα) was found in 90 out of 127 examined stars. From all epochs of observations, 32 objects were classified as CTT (classical T Tauri) and 69 as WTT (weak-line T Tauri) objects. The fraction of the variables in these samples is ~60%. We also identified 20 stellar objects, which showed not only a significant variability of the EW, but which also change their apparent evolutionary stage (CTT ⇆ WTT). For six stars, Hα line was observed in both emission and absorption.

The analysis of data obtained over a wide wavelength range (from X-ray to mid-infrared) has shown that Hα activity and the measure of its variability are in good agreement with the activity of stellar objects measured with its other parameters, such as X-ray radiation and the mass accretion rate. The EW(Hα) differs not only between objects at different evolutionary stages, but also between variable and non-variable objects. The variables in the CTT and WTT samples are more active than non-variables although they have almost the same evolutionary age. Another distinct difference between these variables and non-variables is their average masses. The variables from both CTT and WTT samples are noticeably more massive than non-variables. Our data confirm the assumption made for other star formation regions that the decay of accretion activity occurs more slowly for more massive CTT objects. Apparently, a similar trend is also present in WTT objects, which are at a later stage of evolution. The variability of the stellar objects, which change their evolutionary classes (CTT ⇆ WTT), at least in a fraction of them, is due to the fact that they are close binaries, which affects and modulates their Hα emission activity.

In Panda et al.2018a, we constructed a refined sample from the original Shen et al.(2011) QSO catalog. Based on our hypothesis — the main driver of the Quasar Main Sequence is the maximum of the accretion disk temperature (TBBB) defined by the Big Blue Bump on the Spectral Energy Distribution (Panda et al.2017; Panda et al.2018b). We select the four extreme sources that have RFeII ⩾ 4.0 and use {CIGALE (Boquien et al.2018) to fit their multi—band photometric data. We also perform detailed spectral fitting including the Fe II pseudo—continuum (based on Śniegowska et al.2018)) to estimate and compare the value of RFEII. We show the dependence of FeII strength on changing metallicity.

We address the effect of orientation of the accretion disk plane and the geometry of the broad line region (BLR) in the context of understanding the distribution of quasars along their Main Sequence. We utilize the photoionization code CLOUDY to model the BLR, incorporating the ‘un-constant’ virial factor. We show the preliminary results of the analysis to highlight the co-dependence of the Eigenvector 1 parameter, RFeII on the broad HβFWHM (i.e. the line dispersion) and the inclination angle (θ), assuming fixed values for the Eddington ratio (Lbol/ LEdd), black hole mass (MBH), spectral energy distribution (SED) shape, cloud density (nH) and composition.†

The Event Horizon Telescope (EHT) provides a unique opportunity to probe the physics of supermassive black holes through Very Large Baseline Interferometry (VLBI), such as the existence of the event horizon, the accretion processes as well as jet formation in Low Luminosity AGN (LLAGN). We build a theoretical model which includes an Advection Dominated Accretion Flow (ADAF) and a simple radio jet outflow. The predicted spectral energy distribution (SED) of this model can be compared to observations to get the best estimates of the model parameters. The model-predicted radial emission profiles at different frequency bands can be used to predict whether the inflow can be resolved by the EHT or other telescopes. We have applied this method to some nearby LLAGN such as M84, NGC 4594, NGC 4278 and NGC 3998. We also estimate the model parameters for each of them using high resolution data from different surveys.

In previous works, the radio-X-ray slope in FRI radio galaxies is found to be steeper compared with that in low-luminosity AGNs, indicating different origin of the X-ray emission. Here we reinvestigate this point by compiling a sample of 13 low-excitation radio galaxies (LERG) from 3CR radio catalog of galaxies, where the central engine in LERG is accepted to be a radiatively inefficient accretion flow (RIAF). The core radio and X-ray emissions in all the objects of our sample are detected by VLA/VLBI/VLBA and Chandra/XMM-Newton, respectively. Surprisingly, a shallower slope of Lr – Lx relation () is given by our sample, which demonstrates that the X-ray emission in LERG may come from accretion disk rather than a jet as suggested by previous works. In addition, the slope in the fundamental plane ((log LR = 0.52 log LX + 0.84 log MBH + 10.84) of LERG is found to be well consistent with that reported by Merloni et al. (2003).

A strong quasi-periodic modulation has recently been revealed in the X-ray flux of the X-ray source XMMUJ134736.6+173403. The two observed twin-peak quasiperiodic oscillations (QPOs) exhibit a 3:1 frequency ratio and strongly support the evidence for the presence of an active galactic nucleus black hole (AGN BH). It has been suggested that detections of twin-peak QPOs with commensurable frequency ratios and scaling of their periods with BH mass could provide the basis for a method intended to determine the mass of BH sources, such as AGNs. Assuming the orbital origin of QPOs, we calculate the upper and lower limit on the AGN BH mass M, reaching M ≍ 107–109M⊙. Compared to mass estimates of other sources, XMMUJ134736.6+173403 appears to be the most massive source with commensurable QPO frequencies, and its mass represents the current observational upper limit on the AGN BH mass obtained from the QPO observations.

The interpretation of the main sequence of quasars has become a frontier subject in the last years. This considers the effect of a highly flattened, axially symmetric geometry for the broad line region (BLR) on the parameters related to the distribution of quasars along their main sequence. We utilize the photoionization code CLOUDY to model the BLR, assuming ‘un-constant’ virial factor with a strong dependence on the viewing angle. We show the preliminary results of the analysis to highlight the co-dependence of the Eigenvector 1 parameter, RFeii on the broad Hβ FWHM (i.e. the line dispersion) and the inclination angle (θ), assuming fixed values for the Eddington ratio (Lbo1/LEdd), black hole mass (MBH) and spectral energy distribution (SED) shape. We consider four cases with changing cloud density (nH) and composition. Understanding the emitting region is crucial as this knowledge can be extended to the use of quasars as distance indicators for Cosmology.‡