The study of the differential rotation in the chromosphere of the Sun is of significant importance as it provides valuable insights into the rotational behaviour of the solar atmosphere at higher altitudes and the coupling mechanism between the various layers of the solar atmosphere. In this work, we employed the image correlation technique, explicitly focusing on plages, intending to estimate the chromospheric differential rotation. For this purpose, we have utilized Ca ii K spectroheliograms (1907 – 2007) from the Kodaikanal Solar Observatory (KoSO), recently calibrated with a better technique to ensure accuracy. Our analysis indicates that plages in the chromosphere exhibit faster rotation and a smaller latitudinal gradient when compared to the rotation rate obtained through sunspot tracking. Furthermore, we investigate the temporal analysis of the chromospheric differential rotation parameters across various solar cycles.

Cosmological simulations fail to reproduce realistic galaxy populations without energy injection from active galactic nuclei (AGN) into the interstellar medium (ISM) and circumgalactic medium (CGM); a process called ‘AGN feedback’. Consequently, observational work searches for evidence that luminous AGN impact their host galaxies. Here, we review some of this work. Multi-phase AGN outflows are common, some with potential for significant impact. Additionally, multiple feedback channels can be observed simultaneously; e.g., radio jets from ‘radio quiet’ quasars can inject turbulence on ISM scales, and displace CGM-scale molecular gas. However, caution must be taken comparing outflows to simulations (e.g., kinetic coupling efficiencies) to infer feedback potential, due to a lack of comparable predictions. Furthermore, some work claims limited evidence for feedback because AGN live in gas-rich, star-forming galaxies. However, simulations do not predict instantaneous, global impact on molecular gas or star formation. The impact is expected to be cumulative, over multiple episodes.

Helioseismology has discovered a thin layer beneath the solar surface where the rotation rate increases rapidly with depth. The normalized rotational shear in the upper 10 Mm of the layer is constant with latitude. Differential rotation theory explains such a rotational state by a radial-type anisotropy of the near-surface convection and a short correlation time of convective turbulence compared to the rotation period. The shear layer is the main driver of the global meridional circulation.

The origin of the sub-terahertz (sub-THz) component of radio emission from solar flares, which is characterized by the increase flux with frequency in the 100-400 GHz range, is considered. On the basis of equations of 1D non-LTE radiation hydrodynamics we simulated the altitude distribution of the plasma density and temperature inside the flare loop caused by the interaction of non-stationary beam of accelerated electrons in the form of a triangular pulse with the chromospheric plasma. The FLARIX numerical code was used to calculate the dynamics of the flare plasma parameters at different heights which are compared with the RADYN numerical code. We found that the characteristic heights of the formation of sub-THz emission vary over a wide range with time for both codes. The main contribution to the sub-THz emission comes from the chromospheric and transition region plasma with temperatures of 104–105K.

The first quasars at the Reionisation Epoch, z ∼ 6-7.5, probe the early stages of supermassive black holes and host galaxy assembly. In this paper, we present recent results, exploiting VLT, ALMA and NOEMA observations, that allow us to constrain the onset of strong black hole feedback, the dust properties and star formation rates in high redshift quasars with unprecedented accuracy. These results highlight the strategic importance of ALMA high frequency (i.e. Band 9 and 8) observations to obtain a reliable overview of the host galaxy and supermassive black hole growth out to the highest redshifts.

The Sun and solar-type stars exhibit irregular cyclic variations in their magnetic activity over long time scales. To understand this irregularity, we employed the flux transport dynamo models to investigate the behavior of one solar mass star at various rotation rates. To achieve this, we have utilized a mean-field hydrodynamic model to specify differential rotation and meridional circulation, and we have incorporated stochastic fluctuations in the Babcock–Leighton source of the poloidal field to capture inherent fluctuations in the stellar convection. Our simulations successfully demonstrated consistency with the observational data, revealing that rapidly rotating stars exhibit highly irregular cycles with strong magnetic fields and no Maunder-like grand minima. On the other hand, slow rotators produce smoother cycles with weaker magnetic fields, long-term amplitude modulation, and occasional extended grand minima. We observed that the frequency and duration of grand minima increase with the decreasing rotation rate. These results can be understood as the tendency of a less supercritical dynamo in slow rotators to be more prone to produce extended grand minima. We further explore the possible existence of the dynamo in the subcritical regime in a Babcock–Leighton-type framework and in the presence of a small-scale dynamo.

We here present 0.02–0.04″ resolution ALMA observation of the compact obscured nucleus (CON) of IRAS 17578-0400. A dusty torus within the nucleus, approximately 4 pc in radius, has been uncovered, exhibiting a usually flat spectral index at ALMA band 3, likely due to the millimeter corona emission from the central supermassive black hole (SMBH). The dense gas disk, traced by 13CO(1-0), spans 7 pc in radius and suggests an outflow driven by a disk wind due to its asymmetrical structure along the minor axis. Collimated molecular outflows (CMO), traced by the low-velocity components of the HCN(3-2) and HCO+(3-2) lines, align with the minor axis gas disk. Examination of position-velocity plots of HCN(3-2) and HCO+(3-2) reveals a flared dense gas disk extended a radius of ∼60 pc, infalling and rotating at speeds of about 200 km s−1 and 300 km s−1 respectively. A centrifugal barrier, located around 4 pc from the dynamical center, implies an SMBH mass of approximately 108 Mȯ, consistent with millimeter corona emission estimates. The CMO maintains a steady rotation speed of 200 km s−1 over the 100 pc scale along the minor axis. The projected speed of the CMO is about 80 km s−1, corresponding to around ∼500 km s−1, assuming an inclination angle of 80°. Such a kinematics structure of disk-driven collimated rotating molecular outflow with gas supplies from a falling rotating disk indicates that the feedback of the compact obscured nucleus is likely regulated by the momentum transfer of the molecular gas that connects to both the feeding of the nuclear starburst and supermassive black hole.

The central parsec of AGN is a key region for the launching of winds, and near-infrared interferometry is a unique tool for its study. With GRAVITY at the VLT interferometer, we can now spatially resolve not just the hot dust continuum on milliarcsecond ‘torus’ scales through imaging but also the broad-line region (BLR) on microarcsecond scales through spectro-astrometry. We have mapped the kinematics of the BLR in seven nearby AGN, measured sizes of the hot dust for seventeen AGN, and reconstructed dust images for two AGN. BLR kinematics has allowed us to measure the BLR size and supermassive black hole mass independent of reverberation mapping. The ongoing GRAVITY+ upgrade will greatly enhance the sensitivity and sky coverage of GRAVITY, and first results demonstrate its power for AGN science at z∼2 and beyond.

We present analysis of the evolution of subsurface flows in and around active regions with peculiar magnetic configurations and compare their characteristics with the normal active regions. We also study the zonal and meridional components of subsurface flows separately in different polarity regions separately to better understand their role in flux migration. We use the techniques of local correlation tracking and ring diagrams for computing surface and subsurface flows, respectively. Our study manifests an evidence that the meridional component of the flows near anti-Hale active regions is predominantly equatorward which disagrees with the poleward flow pattern seen in pro-Hale active regions. We also find clockwise or anti-clockwise flows surrounding the anti-Joy active regions depending on their locations in the Southern or Northern hemispheres, respectively.

The radiative mode of AGN feedback, operated through outflows, plays an essential role in the evolution of galaxies. Quasar outflows are detected as blue-shifted broad absorption lines in the UV/optical spectra of quasars. Thanks to the Sloan digital sky survey, 100,000 broad absorption line quasars are available now for ensemble statistical studies. This rich dataset has also enabled us to identify some peculiar cases of these sources. By quantifying the BAL fraction in radio-loud BAL quasars, our studies demonstrate a clear trend of increasing BAL fraction as the viewing angle approaches an edge-on orientation, favoring the orientation model of BAL quasars. Also, by contrasting the properties of BAL quasars with appearing and disappearing BAL troughs, our analysis suggests that the extreme variations in BAL troughs are driven by ionization changes.

We show that, contrary to simple predictions, most AGNs show at best only a small increase of lags with increasing wavelength in the J, H, K, and L bands. We suggest that a possible cause of this near simultaneity from the near-IR to the mid-IR is that the hot dust is in a hollow bi-conical outflow of which we preferentially see the near side. In the proposed model sublimation or re-creation of dust (with some delay relative luminosity variations) along our line of sight in the hollow cone as the flux varies could be a factor in explaining the AGN changing-look phenomenon (CL). Variations in the dust obscuration can help explain changes in relationship of Hβ time delay on Luv variability. The relative wavelength independence of IR lags simplifies the use of IR lags for estimating cosmological parameters.

The canonical undestanding of stellar convection has recently been put under doubt due to helioseismic results and global 3D convection simulations. This “convective conundrum” is manifested by much higher velocity amplitudes in simulations at large scales in comparison to helioseismic results, and the difficulty in reproducing the solar differential rotation and dynamo with global 3D simulations. Here some aspects of this conundrum are discussed from the viewpoint of hydrodynamic Cartesian 3D simulations targeted at testing the rotational influence and surface forcing on deep convection. More specifically, the dominant scale of convection and the depths of the convection zone and the weakly subadiabatic – yet convecting – Deardorff zone are discussed in detail.

Tidal forces in close binaries and multiple systems that contain magnetically active component are supposed to influence the operation of magnetic dynamo. Through synchronization the tidal effect of a close companion helps maintain fast rotation, thus supporting an efficient dynamo. At the same time, it can also suppress the differential rotation of the convection zone, or even force the formation of active longitudes at certain phases fixed to the orbit. V815 Her is a four-star system consisting of two close binaries orbiting each other, one of which contains an active G-type main-sequence star. Therefore, the system offers an excellent opportunity to investigate the influence of gravitational effects on solar-type magnetic activity using different methods.

In the late 80s of the 20th century, Crimean astronomers, studying the structure of transverse magnetic fields in active regions (ARs), discovered signs of the presence of large-scale vertical electric currents – global electric currents (Abramenko, Gopasyuk 1987). In 2018–2020, we finalized and adapted the method for detecting large-scale electric currents to the data of modern instruments for studying the Sun, and began studying their dynamics on time scales of 3–5 days (Fursyak et. al 2020). Our researches carried out during 2020–2023 showed that: 1) Large-scale electric currents with values of the order of ~ 1013 A exist in ARs with nonzero flare activity. 2) Large-scale electric currents extend to the upper layers of the solar atmosphere in one part of the AR, and close through the chromosphere and corona in the remaining part of the AR. This assumption for the AR NOAA 12192 is confirmed by the results of numerical simulations performed in 2016 (Jiang et al. 2016). 3) The greater the magnitude of the large-scale electric current, the higher the probability of occurrence of M- and X- class solar flares in the AR. 4) At the final stages of AR evolution, a nonzero large-scale electric current can have a stabilizing effect on the sunspot, preventing its decay by its own magnetic field. 5) Large-scale electric currents are involved in coronal heating processes. Ohmic dissipation of a large-scale electric current is one of the mechanisms of quasi-stationary heating of coronal plasma above the AR. Our research on large-scale electric currents and the processes in which they take part continues.

We present a theoretical model of the near-surface shear layer (NSSL) of the Sun. Convection cells deeper down are affected by the Sun’s rotation, but this is not the case in a layer just below the solar surface due to the smallness of the convection cells there. Based on this idea, we show that the thermal wind balance equation (the basic equation in the theory of the meridional circulation which holds inside the convection zone) can be solved to obtain the structure of the NSSL, matching observational data remarkably well.

Planetary influence on a stellar convective shell can result in a periodic modulation of stellar dynamo drivers. Similar modulation can arise in stellar binary systems. Using the Parker low-mode dynamo model we investigate the properties of nonlinear parametric resonance. This model is a system of four ordinary differential equations and, in the first approximation, describes the processes of generation and oscillation of large-scale magnetic fields in stellar systems. In the absence of nonlinear suppression effects, the problem, by analogy with a system of harmonic oscillations, allows an asymptotic selection of multiple resonant frequencies. Despite the fact that at first glance at these frequencies it is reasonable to expect an increase in the amplitude, the behavior of the system can be just the opposite. All this stuff deserves a systematic analysis of swing excitation in the dynamo sistems in comparison with classical swing excitation in the framework of the Mathieu equation.

An accurate description of the center-to-limb variation (CLV) of stellar spectra is becoming an increasingly critical factor in both stellar and exoplanet characterization. In particular, the CLV of spectral lines is extremely challenging as its characterization requires highly detailed knowledge of the stellar physical conditions. To this end, we present the Numerical Empirical Sun-as-a-Star Integrator (NESSI) as a tool for translating high-resolution solar observations of a partial field of view into disk-integrated spectra that can be used to test common assumptions in stellar physics.

We have measured zonal and meridional components of subsurface flows up to a depth of 30 Mm below the solar surface by applying the technique of ring diagram on Dopplergrams which are constructed from the spherical harmonic (SH) coefficients. The SH coefficients are obtained from the Helioseismic and Magnetic Imager (HMI) full-disk Dopplergrams. We find a good agreement and some differences between the flows obtained in this study with those from the traditional methods using direct Dopplergrams.

The Sun’s meridional circulation is a crucial component for understanding the Sun’s dynamo and its interior dynamics. However, the determination of meridional circulation is affected by a systematic center-to-limb (CtoL) effect, which introduces systematic errors 5–10 times stronger than the meridional-flow-induced travel-time shifts in deep-flow measurements. Recently, it was found that the CtoL effect has a significant acoustic-frequency dependence, while flow-induced travel-time shifts show little frequency dependence (Chen & Zhao 2018). This discovery forms the basis for designing a new method to remove the CtoL effect. We therefore propose a frequency-dependent approach to measure the CtoL effect and the flow-induced signals in the Fourier domain. In this work, we present this new method and compare time–distance measurements in different frequency bands with those obtained by previous time-domain methods. The results demonstrate consistency with conventional time-domain fitting methods in the dominant frequency range, promising the potential for conducting meridional flow inversion across a broader frequency spectrum.