Papers for Wednesday, Jan 15 2025

Observations suggest an SMBH-host coevolution. We consider the mass and energy flow in a bulge suffused by gases of varying temperatures. By assuming the rate of energy flow independent of the distance from the bulge center and the local virial equilibrium for permeated gases, a key parameter $\varepsilon_b$ was identified that quantifies the rate of mass and energy flow in gases and the efficiency of gas cooling and thus regulates the coevolution of SMBHs and hosts. Using Illustris simulations, we found $\varepsilon_b\propto (1+z)^{5/2}$. A higher $\varepsilon_b$ in the early Universe means a more efficient gas cooling that allows initial rapid growth of SMBHs and hosts. This simple theory, characterized by $\varepsilon_b$, provides the dominant mean cosmic evolution of SMBHs and hosts. All other transient phenomena may only contribute to the dispersion around mean evolution. Relevant scaling laws involving $\varepsilon_b$ were identified. For host galaxies, the mass-size relation $M_b\propto \varepsilon_b^{2/3}r_b^{5/3}G^{-1}$, dispersion-size relation $\sigma_b^2\propto(\varepsilon_b r_b)^{2/3}\propto (1+z)$, or the mass-dispersion relation $M_b\propto \varepsilon_b^{-1}G^{-1}\sigma_b^5$ were identified, where size $r_b\propto (1+z)^{-1}$. For SMBHs, three evolution phases were found involving an initial rapid growth stage with a rising luminosity $L_B\propto (\varepsilon_b M_{BH})^{4/5}$, a transition stage with a declining $L_B\propto \varepsilon_b^2 M_{BH} \propto (1+z)^5$, and a dormant stage with $L_B\propto (\varepsilon_b M_{BH})^{4/3}$. Results suggest a rapid initial super-Eddington growth with a new redshift-dependent luminosity limit $L_X\propto\varepsilon_b^{4/5}M_{BH}^{4/5}G^{-1/5}c$, in contrast to the Eddington limit. Analytical solutions are formulated for the BH and AGN mass functions and AGN duty cycle and predict a slope of -1/5 for the faint-end luminosity function.

Highly collimated outflows (jets) are observed across a wide range of astrophysical systems involving the accretion of material onto central objects, from supermassive black holes in active galaxies to proto-brown dwarfs and stellar-mass black holes. Despite the diversity of their driving sources, it is believed that all jets are different manifestations of a single universal phenomenon. However, a unified explanation for their ejection and collimation remains elusive. In this study we present the first rotation measure analysis of the polarized synchrotron emission ever performed in a protostellar radio jet, which allows us to reveal its true 3D magnetic structure. Unlike extragalactic radio jets, which often exhibit faint counterjets, protostellar radio jets allow both the jet and the counterjet to be analyzed. This exceptional circumstance allows us to unveil the magnetic field structure of both components. Our findings provide the first solid evidence for a helical magnetic field within a protostellar jet, supporting the universality of the jet collimation mechanism.

The kinetic Sunyaev-Zel'dovich effect (kSZ) and patchy screening effect are two complementary cosmic microwave background (CMB) probes of the reionization era. The kSZ effect is a relatively strong signal, but is difficult to disentangle from other sources of temperature anisotropy, whereas patchy screening is weaker but can be reconstructed using the cleaner polarization channel. Here, we explore the potential of using upcoming CMB surveys to correlate a reconstructed map of patchy screening with (the square of) the kSZ map, and what a detection of this cross-correlation would mean for reionization science. To do this, we use simulations and theory to quantify the contributions to this signal from different redshifts. We then use the expected survey properties for CMB-S4 and CMB-HD to make detection forecasts before exploring how the signal depends on reionization parameters. We forecast that CMB-S4 will obtain a hint of this signal at up to 1.8$\sigma$ and CMB-HD will detect it at up to 14$\sigma$. We find that this signal contains unique information about the timing and morphology of reionization, particularly sensitivity to the first half of reionization and to the bispectrum of the ionized gas distribution.

The early evolution of protostellar, star-forming discs, including their density structure, turbulence, magnetic dynamics, and accretion variability, remains poorly understood. We present high-resolution magnetohydrodynamic simulations, using adaptive mesh refinement to capture detailed disc dynamics down to sub-AU scales. Starting from initial conditions derived from a molecular cloud simulation, we model the collapse of a dense core into a protostellar disc over 10,000 yr following sink particle (star) formation, achieving a maximum effective resolution of 0.63 AU. This simulation traces the evolution of the disc density, accretion rates, turbulence, and magnetic field structures. We find that the protostellar disc grows to a diameter of approximately 100 AU, with mass accretion occurring in episodic bursts influenced by the turbulence of the core from which the disc builds up. The disc is highly turbulent with a sonic Mach number of $\sim2$. The magnetic field is wound up by rotational motions, generating a magneto-centrifugal outflow, while episodic accretion events within the disc cause intermittent increases in mass and magnetic energy density, resulting in an equipartition of the thermal and magnetic pressure and an Alfvén Mach number of ~ 2. The disc density profiles steepen over time, following a power law consistent with observed young stellar discs and the minimum mass solar nebula. These results underscore the role of turbulence in early accretion variability and offer new insights into the physical and magnetic structure of young protostellar discs.

Whether the mass of supermassive black hole ($M_\mathrm{BH}$) is directly linked to the quasar radio luminosity remains a long-debated issue, and understanding the role of $M_\mathrm{BH}$ in the evolution of quasars is pivotal to unveiling the mechanism of AGN feedback. In this work, based on a two-component Bayesian model, we examine how $M_\mathrm{BH}$ affects the radio emission from quasars, separating the contributions from host galaxy star formation (SF) and AGN activity. By modelling the radio flux density distribution of Sloan Digital Sky Survey (SDSS) quasars from the LOFAR Two-metre Sky Survey Data Release 2, we find no correlation between $M_\mathrm{BH}$ and SF rate (SFR) at any mass for quasars at a given redshift and bolometric luminosity. The same holds for AGN activity across most $M_\mathrm{BH}$ values; however, quasars with the top 20\% most massive SMBHs are 2 to 3 times more likely to host strong radio jets than those with lower-mass SMBHs at similar redshift and luminosity. We suggest defining radio quasar populations by their AGN and SF contributions instead of radio loudness; our new definition unifies previously divergent observational results on the role of $M_\mathrm{BH}$ in quasar radio emissions. We further demonstrate that this radio enhancement in quasars with the 20\% most massive SMBHs affects only the $\sim5\%$ most radio bright quasars at a given redshift and bolometric luminosity. We discuss possible physical origins of this radio excess in the most massive and radio-bright quasar population, which remains an interest for future study.

We present numerical simulations of dark matter and stellar dynamics in Ultra Light Dark Matter halos tailored to mimic dwarf galaxies. For dark matter particle mass $m\approx 1\times 10^{-22}$ eV, dynamical heating causes the half-light radius to over-shoot surface brightness data of the Fornax galaxy. For $m\approx 1\times 10^{-21}$ eV, soliton core formation leads to a velocity dispersion peak incompatible with kinematics data. Extending the analysis to the Carina and Leo II galaxies, the tension persists up to $m\approx 5\times 10^{-21}$ eV. A caveat in our analysis is the omission of stellar self-gravity. This would not change dynamics today, but could affect extrapolation back in time if the stellar body was more compact in the past.

The vertical shear instability and the streaming instability are two robust sources of turbulence in protoplanetary disks. The former has been found to induce anisotropic turbulence that is stronger in the vertical than in the radial dimension and to be overall stronger compared to the largely isotropic turbulence caused by the streaming instability. In this study, we shed light on the dust diffusion by the vertical shear instability and the streaming instability separately and together, and in particular on the direction- and scale-dependence of the diffusion. To this end, we employ two-dimensional global models of the two instabilities either in isolation or in combination. The vertical shear instability in isolation diffuses dust more strongly in the vertical direction than the streaming instability in isolation, resulting in a wave-shaped dust layer in our two-dimensional simulations. Compared with this large-scale diffusion, though, our study highlights that the vertical shear instability causes substantially weaker or even negligible small-scale diffusion. We validate this result using previously published three-dimensional simulations. In particular when simulating centimetre-sized dust, the undulating dust layer becomes internally razor-thin. In contrast, the diffusion owing to the streaming instability exhibits only a marginal scale-dependence, with the dust layer possessing a Gaussian shape. In models including both instabilities, the undulating mid-plane layer is broadened to a width set by the intrinsic diffusion level caused by the streaming instability.

Our current knowledge of the thermodynamic properties of galaxy clusters comes primarily from detailed studies of clusters selected by their minority components: hot baryons. Most of these studies select the clusters using the component that is being investigated, the intracluster medium (ICM), making the sample choice prone to selection effects. Weak-gravitational lensing allows us to select clusters by the total mass component and, being independent of the type of matter, makes the sample choice unbiased with respect to the baryon content. In this paper, we study four galaxy clusters at intermediate redshift ($0.25<z<0.61$), selected from the weak-lensing survey of Miyazaki et al. (2018). We derive core-excised X-ray luminosities, richness-based masses, Compton parameters, and profiles of mass, pressure and electron densities. These quantities are derived from shear data, Compton maps, and our own X-ray and SZ follow-up. When compared to ICM-selected clusters of the same mass, in the range $2$ to $5 \ 10^{14}$ M$_\odot$, our small sample of four clusters is expected to have on average 0.2 rare ($>2\sigma$) features, while we observed on average two rare features in each one of the seven explored properties: richness, core-excised luminosity, Compton parameter, pressure and electron pressure profiles, and central values of them. The abundance of rare and unique features in such a small sample indicates a fundamental bias in our knowledge of the thermodynamic properties of clusters when derived from ICM-selected samples.

The dynamic properties of molecular clouds are set by the interplay of their self-gravity, turbulence, external pressure and magnetic fields. Extended surveys of Galactic molecular clouds typically find that their kinetic energy ($E_{\rm k}$) counterbalances their self-gravitational energy ($E_{\rm g}$), setting their virial parameter $\alpha_{\rm vir}=2E_{\rm k}/|E_{\rm g}|\approx1$. However, past studies either have been biased by the use of optically-thick lines or have been limited within the solar neighborhood and the inner Galaxy (Galactocentric radius $R_{\rm gc}<R_{\rm gc,\odot} \approx 8$ kpc). Here we present sensitive mapping observations of optically thin $^{13}$CO lines towards molecular clouds in the low-metallicity Galactic outer disk ($R_{\rm gc}\sim9-24$ kpc). By combining archival data from the inner Galaxy and four nearby metal-poor dwarf galaxies, we reveal a systematic trend of $\alpha_{\rm vir}$, which declines from supervirial dynamic states in metal-rich clouds to extremely subvirial dynamic states in metal-poor clouds. In these metal-poor environments, turbulence alone is insufficient to counterbalance the self-gravity of a cloud. A cloud-volumetric magnetic field may replace turbulence as the dominant cloud-supporting mechanism in low-metallicity conditions, for example, the outermost galactic disks, dwarf galaxies and galaxies in the early Universe, which would then inevitably impact the initial conditions for star formation in such environments.

With enough X-ray flux to be detected in a 160s scan by SRG/eROSITA, the $z = 6.19$ quasar CFHQS J142952+544717 is, by far, the most luminous X-ray source known at $z > 6$. We present deep (245 ks) NuSTAR observations of this source; with $\sim180$ net counts in the combined observations, CFHQS J142952+544717 is the most distant object ever observed by the observatory. Fortuitously, this source was independently observed by Chandra $\sim110$ days earlier, enabling the identification of two nearby (30'' and 45'' away), fainter X-ray sources. We jointly fit both Chandra and NuSTAR observations--self-consistently including interloper sources--and find that, to greater than 90% confidence, the observed 3-7 keV flux varied by a factor of $\sim2.6$ during that period, corresponding to approximately two weeks in the quasar rest-frame. This brightening is one the most extreme instances of statistically significant X-ray variability seen in the Epoch of Reionization. We discuss possible scenarios that could produce such rapid change, including X-ray emission from jets too faint at radio frequencies to be observed.

Relics are diffuse, highly-polarized radio sources that trace merger-driven shocks at the periphery of merging galaxy clusters. The LOFAR survey recently discovered a rare example of double relics in the low-mass cluster PSZ2 G181.06+48.47. Through a detailed exploration of new Chandra and XMM-Newton observations, we reveal that PSZ2 G181.06+48.47 has a lower mass ($M_{500,X}=2.57^{+0.37}_{-0.38}\times10^{14}$ M$_{\odot}$) than previously thought. Despite its cool global temperature of $kT_{500}=3.86^{+0.27}_{-0.29}$ keV, PSZ2 G181.06+48.47 is one of the most disturbed clusters in the Planck sample, with a complex morphological and thermodynamic structure. We discover a set of three discontinuities within <500 kpc of the cluster center, and, from a surface brightness analysis, place $5\sigma$ upper limits of $M_{NE}<1.43$ and $M_{SW}<1.57$ for any shock associated with the relic locations. We also revise established scaling relations for double radio-relics by adding 12 new systems not included in previous work. The PSZ2 G181.06+48.47 relics have the widest separation (scaled for $r_{500}$) of all known double-relic systems. The exceptional distance from the cluster center ($>r_{200}$), indicates the relics may be associated with shocks in the ``run-away" phase. We propose that this late-stage, post-apocenter merger is captured as the two subclusters with a mass ratio of 1.2-1.4 fall back into each other. The outer relic shocks were likely produced at the first core passage, while the inner discontinuities are associated with the second infall.

Thermodynamical explosions of White Dwarfs (WD)are one of the keys to high precision cosmology. Nebular spectra, namely mid-infrared (MIR) with JWST are an effective tool to probe for the multi-dimensional imprints of the explosion physics of WDs and their progenitor systems but also pose a challenge for simulations. What we observe as SNe Ia are low-energy photons, namely light curves, and spectra detected some days to years after the explosion. The light is emitted from a rapidly expanding envelope consisting of a low-density and low-temperature plasma with atomic population numbers far from thermodynamical equilibrium. SNe Ia are powered radioactive decays which produce hard X- and gamma-rays and MeV leptons which are converted within the ejecta to low-energy photons. We find that the optical and IR nebular spectra depend sensitively on the proper treatment of the physical conversion of high to low energies. The low-energy photons produced by forbidden line transitions originate from a mostly optically thin envelope. However, the UV is optically thick because of a quasi-continuum formed by allowed lines and bound-free transitions even several years after the explosion. We find that stimulated recombination limits the over-ionization of high ions with populations governed by the far UV. The requirements to simulate nebular spectra are well beyond both classical stellar atmospheres and nebulae. Using our full non-LTE HYDrodynamical RAdiation code (HYDRA) as a test-bed, the sensitivity on the physics on synthetic spectra are demonstrated using observations as a benchmark. At some examples, we establish the power of high-precision nebular spectroscopy as quantitative tool. Centrally ignited, off-center delayed-detonation near Chandrasekhar-mass models can reproduce line-ratios and line profiles of Branch-normal and underluminous SNe Ia observed with JWST.

X-ray observations collected over the last decades have revealed a strongly variable X-ray signal within the Milky Way's Galactic center, interpreted as X-ray echoes from its supermassive black hole, Sgr A*. These echoes are traced by the strong Fe K$\alpha$ fluorescent line at 6.4 keV of which its intensity is proportional to the density of the illuminated molecular gas. Over time, the echo scans through molecular clouds (MCs) in our Galactic center, revealing their 3D structure and highlighting their densest parts. While previous studies have utilized spectral line doppler shifts along with kinematic models to constrain the geometry of the CMZ or to study the structure of individual clouds, these methods have limitations, particularly in the turbulent region of the CMZ. We use archival Chandra X-ray data to construct one of the first 3D representations of one prominent MC, the Stone Cloud, located at (l = 0.068°, b = -0.076°) at a distance of ~20pc from Sgr A* in projection. Using the Chandra X-ray Observatory, we followed the X-ray echo in this cloud from 2008 to 2017. We combine this data with 1.3 mm dust continuum emission observed with the Submillimeter Array (SMA) and the Herschel Space Observatory to re-construct the 3D structure of the cloud and estimate the column densities for each year's observed slice. The analysis of the X-ray echoes along with velocities from SMA molecular line data indicate that the structure of the Stone cloud can be described as a very diffuse background with multiple dense clumps throughout.

Galaxy formation and evolution models, such as semi-analytic models, are powerful theoretical tools for predicting how galaxies evolve across cosmic time. These models follow the evolution of galaxies based on the halo assembly histories inferred from large $N$-body cosmological simulations. This process requires codes to identify halos ("halo finder") and to track their time evolution ("tree builder"). While these codes generally perform well, they encounter numerical issues when handling dense environments. In this paper, we present how relevant these issues are in state-of-the-art cosmological simulations. We characterize two major numerical artefacts in halo assembly histories: (i) the non-physical swapping of large amounts of mass between subhalos, and (ii) the sudden formation of already massive subhalos at late cosmic times. We quantify these artefacts for different combinations of halo finder (SUBFIND, VELOCIRAPTOR, HBT+HERONS) and tree builder codes (D-TRESS+DHALO, TREEFROG, HBT+HERONS), finding that in general more than $50\%$ ($80\%$) of the more massive subhalos with $>10^{3}$ ($>10^{4}$) particles at $z=0$ inherit them in most cases. However, HBT+HERONS, which explicitly incorporates temporal information, effectively reduces the occurrence of these artefacts to $5\%$ ($10\%$). We then use the semi-analytic models SHARK and GALFORM to explore how these artefacts impact galaxy formation predictions. We demonstrate that the issues above lead to non-physical predictions in galaxies hosted by affected halos, particularly in SHARK where the modelling of baryons relies on subhalo information. Finally, we propose and implement fixes for the numerical artefacts at the semi-analytic model level, and use SHARK to show the improvements, especially at the high-mass end, after applying them.

PLAnetary Transits and Oscillations of stars (PLATO) is an ESA M-class mission to be launched by the end of 2026 to discover and characterize transiting planets around bright and nearby stars, and in particular habitable rocky planets hosted by solar-like stars. Over the mission lifetime, an average of 8% of the science data rate will be allocated to Guest Observer programs (GOs) selected by ESA through public calls, hence it is essential for the community to know in advance where the observing fields will be located. In a previous paper, we identified two preliminary long-pointing fields (LOPN1 and LOPS1) for PLATO, respectively in the northern and southern hemisphere. Here we present LOPS2, a slightly adjusted version of the southern field that has recently been selected by the PLATO Science Working Team as the first field to be observed by PLATO for at least two continuous years, following the scientific requirements. In this paper, we describe the astrophysical content of LOPS2 in detail, including known planetary systems, bright/variable/binary stars, clusters and synergies with other current and future facilities.

Using optical and near-infrared images of the Cassiopeia A (Cas A) supernova remnant covering the time period 1951 to 2022, together with optical spectra of selected filaments, we present an investigation of Cas A's reverse shock velocity and the effects it has on the remnant's metal-rich ejecta. We find the sequence of optical ejecta brightening and the appearance of new optical ejecta indicating the advancement of the remnant's reverse shock in the remnant's main shell has velocities typically between 1000 and 2000 km/s, which is ~1000 km/s less than recent measurements made in X-rays. We further find the reverse shock appears to move much more slowly and is nearly even stationary in the sky frame along the remnant's western limb. However, we do not find the reverse shock to move inward at velocities as large as ~2000 km/s as has been reported. Optical ejecta in Cas A's main emission shell have proper motions indicating outward tangential motions ~3500 - 6000 km/s, with the smaller values preferentially along the remnant's southern regions which we speculate may be partially the cause of the remnant's faint and more slowly evolving southern sections. Following interaction with the reverse shock, ejecta knots exhibit extended mass ablated trails 0.2" - 0.5" in length leading to extended emission indicating reverse shock induced decelerated velocities as large as 1000 km/s. Such ablated material is most prominently seen in higher ionization line emissions, whereas denser parts of ejecta knots show surprisingly little deceleration.

The solar tachocline is an internal boundary layer in the Sun located between the differentially-rotating convection zone and the uniformly-rotating radiative interior beneath. Spiegel and Zahn (1992) proposed the first hydrodynamical model, which here we call SZ92, arguing that the tachocline is essentially in a steady state of thermal-wind balance, angular-momentum balance, and thermal equilibrium. Angular momentum transport in their model is assumed to be dominated by strongly anisotropic turbulence, primarily horizontal owing to the strong stable stratification of the radiative interior. By contrast, the heat transport is assumed to be dominated by a predominantly vertical diffusive heat flux owing to the thinness of the tachocline. In this paper, we demonstrate that these assumptions are not consistent with the new model of stratified turbulence recently proposed by Chini et al. (2022) and Shah et al. (2024), which has been numerically validated by Garaud et al. (2024). We then propose a simple self-consistent alternative to the SZ92 model, namely, a scenario wherein angular momentum and heat transport are both dominated by horizontal turbulent diffusion. The thickness of the tachocline in the new model scales as $\Omega_\odot / N_m$, where $\Omega_\odot$ is the mean angular velocity of the Sun, and $N_m$ the buoyancy frequency in the tachocline region. We discuss other properties of the model, and show that it has several desirable features, but does not resolve some of the other well-known problems of the SZ92 model.

Astronomers have used observations of the Galactic gas and dust via infrared, microwave, and radio to study molecular clouds in extreme environments such as the Galactic center. More recently, X- ray telescopes have opened up a new wavelength range in which to study these molecular clouds. Previous flaring events from SgrA* propagate X-rays outwards in all directions, and these X-rays interact with the surrounding molecular gas, illuminating different parts of the clouds over time. We use a combination of X-ray observations from Chandra and molecular gas tracers (line data from Herschel and the Submillimeter Array) to analyze specific features in the Sticks cloud, one of three clouds in the Three Little Pigs system in the Central Molecular Zone (Galactic longitude and latitude of 0.106 and -0.082 degrees respectively). We also present a novel X-ray tomography method we used to create 3-D map of the Sticks cloud. By combining X-ray and molecular tracer observations, we are able to learn more about the environment inside the Sticks cloud.

Millisecond pulsars (MSPs) in binary systems are precise laboratories for tests of gravity and the physics of dense matter. Their orbits can show relativistic effects that provide a measurement of the neutron star mass and the pulsars are included in timing array experiments that search for gravitational waves. Neutron star mass measurements are key to eventually solving the neutron star equation of state and these can be obtained by a measure of the Shapiro delay if the orbit is viewed near edge-on. Here we report on the timing and noise analysis of five MSPs observed with the MeerKAT radio telescope: PSRs J0900$-$3144, J0921$-$5202, J1216$-$6410, J1327$-$0755 and J1543$-$5149. We searched for the Shapiro delay in all of the pulsars and obtain weak detections for PSRs J0900$-$3144, J1216$-$6410, and J1327$-$0755. We report a higher significance detection of the Shapiro delay for PSR J1543$-$5149, giving a precise pulsar mass of $M_{\rm p} = 1.349^{+0.043}_{-0.061}\,$M$_\odot$ and companion white-dwarf mass $M_{\rm c} = 0.223^{+0.005}_{-0.007}\,$M$_\odot$. This is an atypically low mass measurement for a recycled MSP. In addition to these Shapiro delays, we also obtain timing model parameters including proper motions and parallax constraints for most of the pulsars.

The advent of both multi-messenger and time-domain astrophysics over the last decade has seen a large interest in the development of small-scale, cheap, and robust gamma-ray detectors. This has been further encouraged by the availability of CubeSat platforms. Of particular interest are detectors capable of producing spectral and localization measurements of X and gamma-ray transients to allow for accurate follow-up measurements at different wavelengths. A vast number of the instruments developed for such purposes in the last years use a combination of scintillators and Silicon Photomultipliers (SiPMs) for photon detection. Here, we present the design, performance, and space qualification of a readout system capable of reading out 64 SiPM channels. This low-power and low-cost system was originally designed for the POLAR-2 mission, a large scale gamma-ray polarimeter. However, its flexible design makes it equally suitable for use on various CubeSat missions. The system was found to perform well when reading out both plastic and high Z scintillators using a total of 1.8~W. The space qualified design furthermore relies on commercial off-the-shelf components, thereby also removing most international export issues. In this paper, we will present the overall design, the performance of the electronics, its performance when reading out various scintillators and the successful space-qualification of this design.

We measure the void size function (VSF) from the Baryon Oscillation Spectroscopic Survey (BOSS DR16) and perform the cosmological constraints. The BOSS DR16 galaxy sample is selected in the redshift range from $z = 0.2$ to 0.8, considering the selection criteria based on galaxy number density. We identify non-spherical voids from this galaxy catalog using the Voronoi tessellation and watershed algorithm without assuming any void shape. We select the void samples based on the void ellipticity, and derive the VSFs in two redshift bins, i.e. $z=0.2-0.5$ and $0.5-0.8$. The VSF model we use is based on the excursion-set theory, including the void linear underdensity threshold $\delta_{\rm v}$ and the redshift space distortion (RSD) parameter $\beta$. The Markov Chain Monte Carlo (MCMC) method is applied to perform the joint constraints on the cosmological and void parameters. We find that the VSF measurement from BOSS DR16 gives $w = -1.214_{-0.375}^{+0.293}$, $\Omega_{\rm m} = 0.280_{-0.047}^{+0.056}$, and $\sigma_8 = 0.874_{-0.210}^{+0.203}$, which can be a good complementary probe to galaxy clustering measurements. Our method demonstrates the potential of using the VSF to study cosmological models, and it can provide a reference for future VSF analysis in the upcoming galaxy spectroscopic surveys.

Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we have performed an Ursa Major cluster HI Survey (FUMaS) covering the entire UMa region centered at RA=11$^h$59$^m$28$^s$.3, DEC=49\degr05\arcmin18\arcsec with a radius of 7.5\degr. We have obtained the most complete catalog of HIsources in the UMa cluster, containing 179 HI sources with velocities in the range 625-1213.4 km~s$^{-1}$ and masses in the range 10$^{6.0}$-10$^{10.1}$ M$_{\odot}$ assuming a distance of 17.4 Mpc. Among them, 55 HI sources were detected for the first time. 32 HI sources do not have optical counterparts with known redshifts, and we found 25 possible counterparts in the multicolor optical image, and another 7 may be HI clouds. We also detected HI distributions in some interacting systems, e.g. the overlapping gas disks between NGC 3992 and its three companion galaxies, filaments around NGC 4026 and NGC 4111, and debris-like gas around NGC 3998. We computed the HIMF of the UMa cluster using the 1/V$_\mathrm{max}$ method and fitting it with the non-linear least squares (NLLS) and modified maximum likelihood (MML) method for the parameters: log$_{10}$($\phi_*$/Mpc$^{-3}$) = -0.86 $\pm$ 0.18, $\alpha$ = -1.10 $\pm$ 0.08 and log$_{10}$($M_*$/$M_{\odot}$) = 9.92 $\pm$ 0.23 for the NLLS method, and log$_{10}$($\phi_*$/Mpc$^{-3}$) = -0.78 $\pm$ 0.11, $\alpha$ = -1.10 $\pm$ 0.05 and log$_{10}$($M_*$/$M_{\odot}$) = 9.88 $\pm$ 0.14 for the MML method. This result is similar to that derived with VLA data, but the slope is steeper at the low-mass end because we detected more low-mass galaxies. The slope is flatter than that of the global HIMF, which is in agreement with the theoretical prediction that galaxies in high-density regions are stripped of gas due to interactions.

The $H_{0}$ tension stands as a prominent challenge in cosmology, serving as a primary driver for exploring alternative models of dark energy. Another tension arises from measurements of the $ S_{8} $ parameter, which is characterize the amplitude of matter fluctuations in the universe. In this study, we address the alleviation of both the Hubble tension and $ S_{8} $ tension by incorporating Kaniadakis horizon entropy. We investigate two scenarios to explore the impact of this entropy on cosmological parameters. In the first scenario, utilizing modified Friedmann equations through Kaniadakis entropy, we estimate the values of $H_{0}$ and $ S_{8} $. In the subsequent scenario, we introduce the neutrino term and assess its effect on mitigating the Hubble and $ S_{8} $ tensions. Our findings reveal that when considering the first scenario, the results closely align with Planck's 2018 outcomes for Hubble and $ S_{8} $ tensions. Moreover, with the inclusion of neutrinos, these tensions are alleviated to approximately 2$\sigma$, and the $ S_{8} $ value is in full agreement with the results from the KiDS and DES survey. Furthermore, we impose a constraint on the parameter $K$ in each scenario. Our analysis yields $K = 0.12\pm 0.41$ for Kaniadakis entropy without neutrinos and $K = 0.39\pm 0.4$ for the combined dataset considering Kaniadakis entropy in the presence of neutrinos. We demonstrate that the value of K may be affected by neutrino mass, which can cause energy transfer between different parts of the universe and alter the Hubble parameter value.

Planetary nebulae are formed by the matter ejected by low-to-intermediate mass stars (~0.8-8 times the mass of the Sun) towards the end of their lives. As hydrogen and then helium fuel sources run out, stars expand. During these giant phases stars also lose sizable amounts of mass. During the second giant phase, after the exhaustion of core helium, the mass loss is so great that stars lose a large fraction of their mass (50 - 90%), leaving behind a small, hot core, known as a white dwarf, surrounded by a nebula. Planetary nebulae are the result of many processes that shape and alter their ionization structure and chemical composition. The resulting nebula, illuminated by the ultraviolet-rich spectrum of the remnant very hot stellar core, is a spectacle of beauty and science. In this chapter, we show that these objects are invaluable laboratories for astrophysics, astrochemistry, and astromineralogy studies, with impact in many areas of Astronomy.

We present the results of searching for new dwarf galaxies in the Local Volume. We found 40 satellite candidates in the double-virial-radius regions of 20 Milky Way-like and LMC-like galaxies in the southern sky using DESI Legacy Imaging Surveys, 10 of which were known but not clearly associated with the Local Volume previously. Among the 40 satellite candidates, 8 are supposed members of the NGC6744 group and 13 are located in the vicinity of the Sombrero galaxy. Based on seven companions to the giant spiral galaxy NGC\,6744 with measured radial velocities, we estimate that the total mass of the group is $M_T = (1.88\pm0.71)\times 10^{12} M_{\odot}$ and the total mass-to-K-luminosity ratio $M_T/L_K = (16.1\pm6.0) M_{\odot}/L_{\odot}$. We reproduce a distribution of 68 early- and late-type galaxies in the Local Volume situated around the Sombrero, noting their strong morphological segregation and also the presence of a foreground diffuse association of dwarf galaxies at 8 degrees to SE from the Sombrero.

This chapter provides a comprehensive overview of weak gravitational lensing and its current applications in cosmology. We begin by introducing the fundamental concepts of gravitational lensing and derive the key equations for the deflection angle, lensing potential, convergence, and shear. We explore how weak lensing can be used as a cosmological probe, discussing cosmic shear, galaxy-galaxy lensing, and their combination with galaxy clustering in the 3$\times$2pt analysis. The chapter covers the theoretical framework for modeling lensing observables, shear estimation techniques, and major systematic effects such as intrinsic alignments and baryonic feedback. We review the current results of weak lensing cosmology from major surveys and outline prospects for future advancements in the field.

So far, most of the about 5700 exoplanets have been discovered mainly with radial velocity and transit methods. These techniques are sensitive to planets in close orbits, not being able to probe large star--planet separations. $\mu$-lensing is the indirect method that allows us to probe the planetary systems at the snow-line and beyond, but it is not a repeatable observation. On the contrary, direct imaging (DI) allows for the detection and characterization of low mass companions at wide separation (\mbox{$\leq$ 5--6 au}). The main challenge of DI is that a typical planet--star contrast ranges from $10^{-6}$, for a young Jupiter in emitted light, to $10^{-9}$ for Earth in reflected light. In the last two decades, a lot of efforts have been dedicated to combining large (D $\geq$ 5 m) telescopes (to reduce the impact of diffraction) with coronagraphs and high-order adaptive optics (to correct phase errors induced by atmospheric turbulence), with sophisticated image post-processing, to reach such a contrast between the star and the planet in order to detect and characterize cooler and closer companions to nearby stars. Building on the first pioneering instrumentation, the second generation of high-contrast imagers, SPHERE, GPI, and SCExAO, allowed us to probe hundreds of stars (e.g., 500--600 stars using SHINE and GPIES), contributing to a better understanding of the demography and the occurrence of planetary systems. The DI offers a possible clear vision for studying the formation and physical properties of gas giant planets and brown dwarfs, and the future DI (space and ground-based) instruments with deeper detection limits will enhance this vision. In this paper, we briefly review the methods, the instruments, the main sample of targeted stars, the remarkable results, and the perspective of this rising technique.

The James Webb Space Telescope (JWST) is discovering an increasing number of galaxies well into the early stages of the Epoch of Reionization (EoR). Many of these galaxies are clustered with strong Lyman alpha (Ly$\alpha$) emission, motivating the presence of surrounding cosmic HII regions that would facilitate Ly$\alpha$ transmission through the intergalactic medium (IGM). Detecting these HII "bubbles" would allow us to connect their growth to the properties of the galaxies inside them. Here we develop a new forward-modeling framework to estimate the local HII region size and location from Ly$\alpha$ spectra of galaxy groups in the early stages of the EoR. Our model takes advantage of the complementary information provided by neighboring sightlines through the IGM. Our forward models sample the main sources of uncertainty, including: (i) the global neutral fraction; (ii) EoR morphology; (iii) emergent Ly$\alpha$ emission; and (iv) NIRSpec instrument noise. Depending on the availability of complementary nebular lines, $\sim$ 0.006 $\unicode{x2013}$ 0.01 galaxies per cMpc$^3$, are required to be $\gtrsim$95\% confident that the HII bubble location and size recovered by our method is accurate to within $\sim$ 1 comoving Mpc. This corresponds roughly to tens of galaxies at $z\sim7\unicode{x2013}8$ in $\sim$2x2 tiled pointing with JWST NIRSpec. Such a sample is achievable with a targeted survey with completeness down to $M_{\rm UV}^{\rm min}\lesssim$ -19 $\unicode{x2013}$ -17, depending on the over-density of the field. We test our method on 3D EoR simulations as well as misspecified equivalent width distributions, in both cases accurately recovering the HII region surrounding targeted galaxy groups.

Remote observations are often limited by user interfaces, which seem frozen to another computer era: low performances, outdated programming languages, command-line scripting, high version-dependent software. Instead, web-based tools are standard: using nothing more than a browser, astronomers can interact with a generic observatory in a native cross-platform, remote-born way. We used this approach while advancing in the remotization of the 1m-class OARPAF telescope, located in Northern Italy. The web-based control software provides easy and integrated management of its components. This solution can be exported not only to similar hardware/software facilities, but also to large instruments such as SHARK-NIR at LBT; not only for single operations, but also for procedure scripting. In this contribution we describe our best practices and present two recent, orthogonal use cases: an in-place professional use for exoplanetary transit follow-ups, and the first remote control of the telescope from a Japanese high school, allowing students to independently observe, in their daytime, globular clusters

In the context of SHARK-NIR (System for coronagraphy with High Order adaptive optics in Z and H band), we present the development of SHINS, the SHARK-NIR INstrument control Software, in particular focusing on the changes introduced during the Assembly, Integration, and Test (AIT) phase. SHARK-NIR observing sessions will be carried out with "ESO-style" Observation Blocks (OBs) based on so-called Templates scripts that will be prepared by observers. We decided to develop Templates also for the large number of AIT tests (flexures, coronagraphic mask alignment, scientific camera performances...). Here we present the adopted HTTP API for the OBs generation and a web-based frontend that implements it. Taking advantage of this approach, we decided to expose APIs also for individual device movement and monitoring, as well as for general status. These APIs are then used in the web-based instrument control and synoptic panels. During the recent AIT phase, a potential collision issue between two motorized components emerged. While we are exploring the possibility of a hardware interlock, we present a software solution developed at the Observation Software level, that is also available while using other software such as engineering panels. The system is based on three protection layers and it has been successfully tested.

A 2.1-year periodic oscillation of the gamma-ray flux from the blazar PG 1553+113 has previously been tentatively identified in almost 7 year of data from the Fermi Large Area Telescope. After 15 years of Fermi sky-survey observations, doubling the total time range, we report >7 cycle gamma-ray modulation with an estimated significance of 4 sigma against stochastic red noise. Independent determinations of oscillation period and phase in the earlier and the new data are in close agreement (chance probability <0.01). Pulse timing over the full light curve is also consistent with a coherent periodicity. Multiwavelength new data from Swift X-Ray Telescope, Burst Alert Telescope, and UVOT, and from KAIT, Catalina Sky Survey, All-Sky Automated Survey for Supernovae, and Owens Valley Radio Observatory ground-based observatories as well as archival Rossi X-Ray Timing Explorer satellite-All Sky Monitor data, published optical data of Tuorla, and optical historical Harvard plates data are included in our work. Optical and radio light curves show clear correlations with the gamma-ray modulation, possibly with a nonconstant time lag for the radio flux. We interpret the gamma-ray periodicity as possibly arising from a pulsational accretion flow in a sub-parsec binary supermassive black hole system of elevated mass ratio, with orbital modulation of the supplied material and energy in the jet. Other astrophysical scenarios introduced include instabilities, disk and jet precession, rotation or nutation, and perturbations by massive stars or intermediate-mass black holes in polar orbit.

OARPAF (Osservatorio Astronomico Regionale Parco Antola Comune di Fascia, Italy) hosts an 80cm Astelco telescope with a Gambato Dome, SBIG-STX camera, Davis weather station, and SBIG AllSky camera. We present a layer-structured python3 framework to control these devices. Layer 1 provides straightforward getter/setter interface for "atomic" operations on devices. Layer 2 wraps the above mentioned atomic operations into "ESO-style" Templates, to perform sequences of common pointing, observation, and calibration operations called "Observation Blocks" (OBs) that are run by a sequencer. Layer 3 is a REST API based on HTTP verbs to expose methods that control Layer 1 devices and Layer 2. We also present a web interface built on top of this layer. The work is part of the frame for remoting and robotizing the observatory.

The discovery of the solar activity cycle was linked from the outset to the observation of the temporal variability of sunspots, which we know to be the result of complex processes associated with the dynamics of inner layers. Numerous recent studies have highlighted changes in the Sun's Near-Surface Shear Layer (NSSL), pointing to the role of the leptocline, a shallow and sharp rotational shear layer in the top around 8 Mm. The leptocline, mainly characterized by a strong radial rotational gradient at middle latitudes and self-organized meridional flows, is the cradle of numerous phenomena: opacity, superadiabaticity, and turbulent pressure changes; the hydrogen and helium ionization processes; a sharp decrease in the sound speed; and, probably, variations of the seismic radius associated with a nonmonotonic expansion of subsurface layers with depth. In addition, the leptocline may play a key role in forming the magnetic butterfly diagram. Such results are a starting point for further systematic investigations of the structure and dynamics of this layer, which will lead to a better understanding of solar activity.

The work investigates the features of galactic cosmic ray density and anisotropy behavior and their relation to solar sources, interplanetary and geomagnetic disturbances from May 8 to May 13, 2024. During this time, powerful solar flares and fast CMEs were recorded, leading to registration of an extreme geomagnetic storm along with one of the most significant Forbush effects for the entire observation period. All the calculations of cosmic ray characteristics are made using the data of global neutron monitor network and unique methods maintained at IZMIRAN: the Global Survey Method and the Ring of Stations Method. It is determined that the magnitude of Forbush effect under study was 15.7% (for particles with 10 GV rigidity) and as an extreme geomagnetic storm was recorded there was a significant magnetospheric effect observed in the data of neutron monitors (~4%).

After the first multi-messenger observation of a binary neutron star (BNS) merger powering a short-duration gamma-ray burst (GRB), GW170817-GRB 170817A, remarkable effort is ongoing to unravel the evolution of the collimated, relativistic outflow (or jet) that was launched during the merger and fed the GRB event, imprinting its angular structure onto the follow-up afterglow signal. Current theoretical models, based on relativistic magneto-hydrodynamic (RMHD) simulations, offer detailed insights into the launch and propagation processes that govern jet evolution. Notably, these simulations point out that jet injection parameters, such as luminosity, magnetization, power decay time scale, and launch time relative to merger, play a crucial role. However, the impact of these parameters is typically investigated within simplified jet propagation environments, lacking a direct connection with a realistic BNS merger aftermath. In this work, we present the first suite of 3D RMHD simulations exploring the influence of such parameters on the propagation of magnetized incipient GRB jets injected into magnetized environments directly imported from the outcome of a general-relativistic MHD BNS merger simulation. Our results demonstrate that, alongside the injection parameters, the BNS merger environment has a central role in shaping the overall jet evolution. Specifically, under identical jet parameters, the fate of an incipient jet (whether it successfully breaks out or becomes choked) depends strongly on the properties of such environment. Further quantitative comparison between realistic and simplified environments reveals major differences, emphasizing the importance of incorporating the former for accurate modeling.

AR Aur A+B is a close binary of astrophysical interest, because dissimilar surface compositions are reported between similar late B-type dwarfs. A new spectroscopic study on this system was carried out based on the disentangled spectra, in order to determine their atmospheric parameters and elemental abundances, The effective temperature and microturbulence (determined from the equivalent widths of Fe II lines) turned out (11150K, 0.9km/s) and (10650K, 0.1km/s) for A and B. The chemical abundances of 28 elements were derived while taking into account the non-LTE effect for Z<=15 elements (Z: atomic number). The following trends were elucidated for [X/H] (abundance of X relative to the Sun): (1) Qualitatively, [X/H] shows a rough global tendency of increasing with Z, with the gradient steeper for A than for B. (2) However, considerable dispersion is involved for A, since prominently large peculiarities are seen in specific elements reflecting the characteristics of HgMn stars (e.g., very deficient N, Al, Sc, Ni; markedly overabundant P, Mn). (3) In contrast, the Z-dependence of [X/H] for B tends to be nearly linear with only a small dispersion. These observational facts may serve as a key to understanding the critical condition for the emergence of chemical anomaly.

Our study presents a time-resolved X-ray spectral analysis of Mkn 421, using AstroSat observations taken during different epochs between 2016 and 2019. The variability of the source in X-rays is utilized to investigate the evolution of its spectral properties. Each observation period was divided into segments of about 10 ks, and we employed three forms of particle distributions: broken-power law (BPL), log-parabola (LP), and power-law with maximum electron energy (xi-max model) undergoing synchrotron losses to fit the broad X-ray spectrum in each segment. We observed that all of these models provided good fits to the spectra. In the case of the broken-power law model, we investigated the relationship between normalized particle density at an energy less than the break energy and the index before the break. The results revealed an inverse correlation between the index and particle density with no time delay. Additionally, correlations between spectral parameters were used to determine the pivot energy. We observed that the pivot energy remained the same across the observations. For xi-max and LP models, we define analogous pivot energies and show that they also do not vary, indicating the model-independent nature of the result. The constant pivot energy suggests that the source's variability arises from index variations and not due to changes in the normalization. Consequently, parameters such as magnetic field strength, Doppler factor, etc., do not contribute to the source's variability. Instead, variations are primarily associated with the acceleration or escape timescales of emitted particles within the source.

We search for reflection-dominated Compton-thick active galactic nuclei (CT AGN) candidates in the Lockman Hole region using the data of SRG/eROSITA Lockman Hole survey. We selected sources with anomalously hard photon indices in the $0.3 - 8.0$ keV band, untypical for type I AGN. In particular, we required that the upper end of the $90\%$ error interval did not exceed a fiducial boundary of $\Gamma=1.3$. We found 291 sources which constitute a rare subpopulation among extragalactic X-ray sources detected by eROSITA in the Lockman Hole field, $\approx 5\%$. These sources constitute the eROSITA sample of CT AGN candidates in the Lockman Hole field. We further divide the sources into three categories depending on the availability of reliable redshift and statistically significant detection of intrinsic absorption. We present two catalogues: the bright sample (37 sources) and the faint one (254). We estimate the fraction and sky density of reflection-dominated CT AGN candidates. We show examples of individual spectra and use stacking analysis to search for possible redshift evolution of their properties with redshift. We analyse combined eROSITA spectra of bright sources of different categories with a physically motivated spectral model UXCLUMPY and find them fully consistent with the fits to the about $\sim 1$ Msec XMM-Newton data for one of our reflection-dominated CT candidates, Type 2 galaxy $\text{SRGe J105348.6+573032}$. The catalogues of CT AGN candidates could be a good starting point for planning future studies and follow-ups at all wavelengths.

Open star clusters are dynamic systems whose evolution is critically influenced by initial conditions such as star formation efficiency and orbital parameters. Understanding their dissolution mechanisms provides insight into stellar population dynamical mixing in the Milky Way. We aim to investigate the dynamical evolution and dissolution of initially non-virialised open clusters by examining how different global star formation efficiencies and orbital characteristics impact the cluster longevity and structural changes. We followed the evolution of the clusters up to their dissolution time on the basis of our calculations. We compare our open cluster dynamical evolutionary models with the observed open clusters in our Galaxy's solar vicinity. Using high-order direct N-body simulations, we modelled cluster evolution across different Galactic orbits, systematically varying initial star formation efficiencies to comprehensively explore dissolution mechanisms. Our simulations reveal that open clusters typically survive approximately ten orbital periods, with cluster lifetime being strongly dependent on global star formation efficiency and only marginally influenced by orbital eccentricity. We estimate gas expulsion timescales of 0.9 Myr, with initial supernova explosions efficiently removing gaseous components from the cluster. The expected lifetime of the cluster (in units of orbital periods) strongly depends on the cluster global star-formation efficiency and only slightly on the orbital eccentricities of the cluster. The theoretical models demonstrate a remarkable agreement of the Roche-volobe filling parameter with the recent observed Gaia DR3 cluster catalogues in the solar vicinity. By incorporating a mixed sample of clusters with varying star formation efficiencies, we provide a more nuanced understanding of open cluster evolution in the Galactic disc.

We study the mass growth histories of the halos of Milky Way and M31 analogues formed in constrained cosmological simulations of the Local Group. These simulations constitute a fair and representative set of $\Lambda$CDM realisations conditioned on properties of the main Local Group galaxies, such as their masses, relative separation, dynamics and environment. Comparing with isolated analogues extracted from the TNG dark-matter-only simulations, we find that while our M31 halos have a comparable mass growth history to their isolated counterparts, our Milky Ways typically form earlier and their growth is suppressed at late times. Mass growth associated to major and minor mergers is also biased early for the Milky Way in comparison to M31, with most accretion occurring 1 - 4 Gyr after the Big Bang, and a relatively quiescent history at later times. 32% of our Milky Ways experienced a Gaia-Enceladus/Sausage (GES)-like merger, while 13% host an LMC-like object at the present day, with 5% having both. In one case, an SMC- and a Sagittarius-analogue are also present, showing that the most important mergers of the Milky Way in its Local Group environment can be reproduced in $\Lambda$CDM. We find that the material that makes up the Milky Way and M31 halos at the present day first collapsed onto a plane roughly aligned with the Local Sheet and Supergalactic plane; after $z \sim 2$, accretion occurred mostly within this plane, with the tidal effects of the heavier companion, M31, significantly impacting the late growth history of the Milky Way.

Global 3-D magnetohydrodynamical simulations have been conducted to study magnetospheric accretion around stars with various spin rates. For slow rotators, characterized by a fastness parameter $\omega_s\lesssim 0.78$, the disk's inner edge at the magnetospheric truncation radius becomes unstable to the interchange instability, leading to intruding filaments which produce hot spots closer to the stellar equator. Depending on spin rate, slow rotators can be in ``chaotic'' or ``ordered'' unstable regimes. For fast rotators, the interchange instability is suppressed by the super-Keplerian rotation beyond the corotation radius, and hot spots are generated only through polar accretion. Low- and mid-energy flux hot spots cover $\lesssim20\%$ and $\lesssim3\%$ of the surface, with faster rotators tending to produce hotter spots. Beyond the truncation radius, angular momentum transfers from the disk surface to the midplane, resulting in surface accretion and midplane outflow. The midplane outflow may transport thermally processed materials (e.g. those in chondrites) to the outer disk. Field inflation generates episodic winds with mass-loss rates $\sim 1\%-40\%$ of the accretion rate, depending on stellar spin. Frequent magnetic reconnections lead to efficient star-disk coupling. We derive the torque exerted by the disk on the star as a function of stellar spin. For fast rotators/propellers, both spin-down torque and disk wind rate increase dramatically with stellar spin. The equilibrium spin state occurs at $\omega_s\sim0.7$, with wind/jet speeds ($\sim$500 km/s) and mass loss rates ($\sim10\%$ accretion rate) aligning with observations. Most results are insensitive to disk thickness. Finally, we present testable predictions for how observables vary with stellar spin.

We investigate the behavior of multiple stellar populations in globular clusters (GCs) using photometric diagrams constructed with Euclid photometry. By employing synthetic spectra and isochrones that incorporate the chemical differences between first-population (1P) stars, resembling field stars, and second-population (2P) stars, enriched in helium and nitrogen but depleted in carbon and oxygen, we identify, from a theoretical perspective, the color-magnitude diagrams and the chromosome maps most effective at distinguishing these populations within GCs. Euclid photometry proves to be a powerful tool for identifying multiple populations among M-dwarfs, as 1P and 2P stars form distinct sequences in well-chosen photometric diagrams, driven by differences in the strength of oxygen-based molecular features, such as water vapor. To validate our theoretical findings, we analyzed Euclid photometry and astrometry of the GC NGC 6397, complemented by photometric and astrometric data from the Hubble Space Telescope and James Webb Space Telescope, enabling a comprehensive study of its stellar populations across a wide field of view. We find that the 1P constitutes about 30% of the M-dwarfs in NGC 6397, with the fraction of 1P stars remaining consistent across different stellar masses and throughout the entire field of view. 2P stars exhibit an [O/Fe] depletion of about 0.3 dex relative to 1P stars, and both populations display isotropic proper motions. This study represents the first comprehensive analysis of multiple populations among M-dwarfs across a wide field of view, demonstrating that Euclid photometry is a powerful instrument for investigating multiple populations in GCs.

We present radio polarimetric images of 12 Seyfert and Low-Ionization Nuclear Emission-line Region (LINER) galaxies belonging to the Centre for Astrophysics (CfA)+12 micron sample exhibiting kiloparsec-scale radio outflows (KSRs). These observations have been carried out at 10 GHz with Karl G. Jansky Very Large Array (VLA) in D-array and at 1.4 GHz with the BnA$\rightarrow$A array configurations. We find signatures of organized magnetic (B-) field structures in the cores, jets and lobes of these galaxies. The linear polarization fraction varies from a few per cent in the cores to $47\pm18$ per cent in the lobes. The inferred B-fields are toroidal in the cores of several sources making them consistent with the presence of either a sheath-like or a wind-like component surrounding the jet. The in-band spectral index images typically show the presence of flat/inverted spectrum cores and steep spectrum lobes. Radio cores with flatter spectra are found to have lower Eddington ratios while the steeper ones have higher. A strong correlation is observed between the Seyfert/LINER radio outflow properties and the mass of the supermassive black holes (SMBHs); correlations with Eddington ratios are weaker. We find signatures of jet-medium interaction and both positive and negative AGN feedback in these sources. Overall, our study indicates that radio-quiet (RQ) AGN with KSRs possess radio outflows driven by magnetic fields anchored to their black holes - accretion disks, which significantly impact their environments.

Planck data to study the Virgo Cluster's galaxy M 90 and its surroundings. We find, as in the case of certain galaxies of the Local Group and its vicinity, a substantial temperature asymmetry that probably arises from the rotation of the M 90 halo and extends up to about one degree from its centre. This temperature asymmetry is particularly remarkable as M 90 is a rare blueshifted galaxy of the Virgo Cluster, and it thus has possible implications for the cluster internal dynamics versus the galactic halo's formation and structure. Possible explanations for the observed effect are discussed.

We present results from fitting $p$-mode spectra derived from 7-d segments of Sun-as-a-star helioseismic observations from the Birmingham Solar Oscillations Network covering 32 yr. The results show a clear dependence of the mode frequencies on solar activity, and the frequency dependence of the sensitivity to activity can also be seen. Because we use data segments that cover less than half of a solar rotation, we are able to test for the effect of activity on the solar far side. By fitting with a model that takes into account activity on the far side of the Sun, we show that the frequency shifts are sensitive to activity from the whole Sun, not just the side facing the observer. Our results suggest that there is potential to investigate activity-related asteroseismic frequency shifts in solar-like oscillators using short time series of observations.

This paper presents a classification framework for the architectures of planetary systems based on a complete survey of the confirmed exoplanet population. With nearly 6000 confirmed exoplanets discovered, including more than 300 multiplanet systems with three or more planets, the current observational sample has reached the point where it is both feasible and useful to build a classification system that divides the observed population into meaningful categories. This framework provides a criterion to split planetary systems into inner and outer regimes, and then further divides inner systems into dynamical classes. The resulting categories include "peas-in-a-pod systems" with uniformly small planets and "warm Jupiter systems" with a mix of large and small planets, as well as "closely-spaced systems" and "gapped systems," with further subdivisions based on the locations of gaps and other features. These categories can classify nearly all of the confirmed systems with three or more planets with minimal ambiguity. We qualitatively examine the relative prevalence of each type of system, subject to observational selection effects, as well as other notable features such as the presence of hot Jupiters. A small number of outlier systems are also discussed. Potential additional classes of systems yet to be discovered are proposed.

We used a supervised machine learning algorithm (probabilistic random forest) to classify ~130 million sources in the VISTA Survey of the Magellanic Clouds (VMC). We used multi-wavelength photometry from optical to far-infrared as features to be trained on, and spectra of Active Galactic Nuclei (AGN), galaxies and a range of stellar classes including from new observations with the Southern African Large Telescope (SALT) and SAAO 1.9m telescope. We also retain a label for sources that remain unknown. This yielded average classifier accuracies of ~79% (SMC) and ~87% (LMC). Restricting to the 56,696,719 sources with class probabilities (P$_{class}$) > 80% yields accuracies of ~90% (SMC) and ~98% (LMC). After removing sources classed as 'Unknown', we classify a total of 707,939 (SMC) and 397,899 (LMC) sources, including > 77,600 extragalactic sources behind the Magellanic Clouds. The extragalactic sources are distributed evenly across the field, whereas the Magellanic sources concentrate at the centres of the Clouds, and both concentrate in optical/IR colour-colour/magnitude diagrams as expected. We also test these classifications using independent datasets, finding that, as expected, the majority of X-ray sources are classified as AGN (554/883) and the majority of radio sources are classed as AGN (1756/2694) or galaxies (659/2694), where the relative AGN-galaxy proportions vary substantially with radio flux density. We have found: > 49,500 hitherto unknown AGN candidates, likely including more AGN dust dominated sources which are in a critical phase of their evolution; > 26,500 new galaxy candidates and > 2800 new Young Stellar Object (YSO) candidates.

Radio searches for single pulses provide the opportunity to discover one-off events, fast transients and some pulsars that might otherwise be missed by conventional periodicity searches. The MeerTRAP real-time search pipeline operates commensally to observations with the MeerKAT telescope. Here, we report on 26 new Galactic radio transients, mostly rotating radio transients (RRATs) and also the detection of one RRAT and two pulsars that were independently discovered by other surveys. The dispersion measures of two of the new sources marginally exceed the Galactic contribution depending on the electron density model used. Using a simple method of fitting a Gaussian function to individual pulses, and obtaining positions of arcsecond accuracy from image-based localisations using channelised voltage data from our transient buffer, we have derived timing solutions spanning multiple years for five sources. The timing parameters imply ages of several Myr and low surface magnetic field strengths which is characteristic of RRATs. We were able to measure spin periods for eight more transients, including one source which appears to rotate every 17.5 seconds. A majority of the sources have only been seen in one observation, sometimes despite multiple return visits to the field. Some sources exhibit complex emission features like component switching and periodic microstructure.

Compact, star-forming galaxies with high star formation rate surface densities ($\Sigma_{\text{SFR}}$) are often efficient Lyman continuum (LyC) emitters at $z\leq 4.5$, likely as intense stellar feedback creates low-density channels that allow photons to escape. Irregular or disturbed morphologies, such as those resulting from mergers, can also facilitate LyC escape by creating anisotropic gas distributions. We investigate the influence of galaxy morphology on LyC production and escape at redshifts $5 \leq z \leq 7$ using observations from various \textit{James Webb Space Telescope} (JWST) surveys. Our sample consists of 436 sources, which are predominantly low-mass ($\sim 10^{8.15} M_\odot$), star-forming galaxies with ionizing photon efficiency ($\xi_{\rm ion}$) values consistent with canonical expectations. Since direct measurements of $f_{\rm esc}$ are not possible during the Epoch of Reionization (EoR), we predict $f_{\rm esc}$ for high-redshift galaxies by applying survival analysis to a subsample of LyC emitters from the Low-Redshift Lyman Continuum Survey (LzLCS), selected to be direct analogs of reionization-era galaxies. We find that these galaxies exhibit on average modest predicted escape fractions ($\sim 0.04$). Additionally, we assess the correlation between morphological features and LyC emission. Our findings indicate that neither $\xi_{\rm ion}$ nor the predicted $f_{\rm esc}$ values show a significant correlation with the presence of merger signatures. This suggests that in low-mass galaxies at $z \geq 5$, strong morphological disturbances are not the primary mechanism driving LyC emission and leakage. Instead, compactness and star formation activity likely play a more pivotal role in regulating LyC escape.

As the sample of mid-inclination disks with measured CO emission surfaces grows, a fundamental unanswered question is how these vertical profiles connect to their host properties. This project aims to relate the vertical extent of protoplanetary disks as traced by $^{12}$CO $2-1$ to key stellar and physical parameters. In order to produce a result that is applicable towards an observational analysis, we benchmark our results with ALMA observations of CO emission from nineteen disks. We produce a grid of disk models using the physical-chemical code DALI, for a template T Tauri and Herbig star. Our models use an iterative solver to calculate the hydrostatic equilibrium equations and determine a physically-motivated density structure. Key stellar and disk parameters such as stellar luminosity and temperature, total disk mass, carbon abundance and critical radius are varied to determine their effect on the CO emitting surface. Each vertical profile is fitted by an exponentially tapered power-law and characterized by the $z/r$ value that represents the structure inwards of 80% of the tapering radius. The CO emission surface location is primarily determined by the disk mass ($M_d$) and the level of volatile carbon depletion. T Tauri and Herbig systems show different vertical profiles, with disks around T Tauri stars being more vertically extended. We derive a $z/r$-$M_d$ relationship, which for each stellar type has a degeneracy with the volatile carbon abundance. In order to reconcile total disk mass estimates from the characteristic $z/r$ and the values obtained based on dust continuum analysis, a volatile carbon depletion of 10-100 (with respect to the ISM) is needed for the majority of our sources. Our carbon depletion values are in agreement with previous literature estimates, highlighting the potential of this method to rapidly calculate key disk parameters.

The disintegrating ultra-short period rocky exoplanet K2-22b periodically emits dusty clouds in a dynamically chaotic process resulting in a variable transit depth from 0-1.3%. The effluents that sublimate off the surface and condense out in space are probably representative of the formerly interior layers convectively transported to the molten surface. Transmission spectroscopy of these transiting clouds reveal spectral fingerprints of the interior composition of this rocky world. We used JWST's Mid-Infrared Instrument (MIRI) as a low-resolution slitless spectrograph to observe four predicted transit windows for K2-22b. For each observation, we extracted a transmission spectrum over the spectral range of 4.3-11.8 $\mu$m. We detect one transit at high significance and two at low significance. We find that the data 1) disfavor featureless, iron-dominated core material, 2) are consistent with some form of magnesium silicate minerals, likely from mantle material, and 3) show a distinct and unexpected feature at $\sim$5 $\mu$m. The unexpected feature, also seen weakly in the low-significance transits, is consistent with some gas features, possibly NO and/or CO$_2$. These findings warrant further study to improve the constraints on the composition of this disintegrating rocky world.

The largest electric fields between 18 and 30 solar radii are in narrow band waves simultaneously observed at a few Hz (somewhat above the local proton gyrofrequency) and a few hundred Hz (far below the lower hybrid frequency), with the higher frequency wave triggered at specific phases of the lower frequency wave. This wave pair, called triggered ion acoustic waves (TIAW), has been shown to both be physical and to occur at times of electron heating. A theory of electron heating and acceleration by the low frequency wave has been presented. While this theory and the TIAW results strongly suggest the presence of low frequency electric fields that are parallel to the local magnetic field, such fields have not been directly observed. In this paper, such parallel electric field observations are reported and TIAW are further described to conclude that they occur during about 75 percent of the Parker Solar Probe passes through 18 to 30 solar radii and, when present, they are the dominant wave signal, lasting for hours. In the presence of these parallel electric fields, electrons are heated while, in their absence, there is no electron heating. That there is no heating between 18 and 30 solar radii in the absence of TIAW is a most significant result because it invalidates other proposed mechanisms that predict heating in this radial range all of the time.

Transition-edge sensor (TES) is a highly sensitive device that is capable of detecting extremely low levels of energy. It is characterised by low noise performance and high energy resolution. A mature method for reading out TES signal is through time-division multiplexing (TDM) direct current superconducting quantum interference device (SQUID). In a TDM system, the signal readout chain represents a significant source of noise other than the TES intrinsic noise. The noise generated by TES is in the range of several tens to several hundreds of $pA/\sqrt{Hz}$. In order to ensure the high energy resolution of TES, it is necessary to ensure that the noise contribution from the room temperature readout electronics is less than $10$ $pA/\sqrt{Hz}$ above 100 $Hz$. In this work, we have designed a low-noise, high-resolution room temperature readout circuit for TDM. The equivalent current noise contribution of ADC is about $0.05$ $pA/\sqrt{Hz}$ above 100 $Hz$ and $0.46$ $pA\sqrt{Hz}$ under 30:1 multiplexing. The resolution of the analog to digital converter (ADC) is larger than 11.5 bits, which can reconstruct the TES signal without distortion. The readout board, which has eight channels, has JESD204B serial ports, which has greatly simplified the space of room temperature electronics. The readout chain is based on multi-threaded CPU processing and can transfer data at 2 $Gbps$ for each channel in real time. This readout board can be used in a TDM system with smaller size for large TES arrays.

Neutrino self-interactions beyond the standard model have profound implications in astrophysics and cosmology. In this work, we study an uncharted scenario in which one of the three neutrino species has a mass much smaller than the temperature of the cosmic neutrino background. This results in a relativistic component that significantly broadens the absorption feature on the astrophysical neutrino spectra, in contrast to the sharply peaked absorption expected in the extensively studied scenarios assuming a fully nonrelativistic cosmic neutrino background. By solving the Boltzmann equations for neutrino absorption and regeneration, we demonstrate that this mechanism provides novel sensitivity to sub-keV mediator masses, well below the traditional $\sim 1$--100 MeV range. Future observations of the diffuse supernova neutrino background with Hyper-Kamiokande could probe coupling strengths down to $g \sim 10^{-8}$, surpassing existing constraints by orders of magnitude. These findings open new directions for discoveries and offer crucial insights into the interplay between neutrinos and the dark sector.

We study the non-Gaussianities generated by light axions, or compact scalar fields, during inflation. To correctly calculate their impact on primordial statistics, we will argue that it is necessary to account for the periodicity, or gauge symmetry, of the compact scalars. We illustrate this point by comparing the predictions for the squeezed kinematic limit of the primordial bispectrum generated by two cases: a non-compact scalar $\sigma$ and a compact scalar $\varphi$. We demonstrate that while a light non-compact scalar predicts a bispectrum of the so-called local shape, the light compact scalar predicts a qualitatively different shape characterised by the ratio of the Hubble scale to its field-space circumference. In doing so, we show that ignoring the gauge symmetry of the compact scalar during inflation leads to spurious infrared enhancements which are softened by working with appropriate gauge-invariant operators. In addition, we connect our results for the primordial bispectrum with late-time cosmological observables and show that it is possible to measure the decay constant of the compact scalar using galaxy clustering measurements.

Core-collapse supernovae (SNe) provide a unique environment to study Feebly Interacting Particles (FIPs) such as Axion-Like Particles (ALPs), sterile neutrinos, and Dark Photons (DPs). This paper focuses on heavy FIPs produced in SNe, whose decay produces electrons and positrons, generating observable secondary signals during their propagation and annihilation. We focus on the In-flight Annihilation (IA) of positrons, which emerge as the most significant contribution to the resulting gamma-ray spectrum. Using data from COMPTEL and EGRET we derive the most stringent bounds on the FIP-electron couplings for heavy ALPs, sterile neutrinos, and DPs. These results strenghten existing bounds of one or two orders of magnitude, depending on the FIP model.

The search for axions sits at the intersection of solving critical problems in fundamental physics, including the strong CP problem in QCD, uncovering the nature of dark matter, and understanding the origin of the universe's matter-antimatter asymmetry. The measurement of axion-mediated spin-dependent interactions offers a powerful approach for axion detection. However, it has long been restricted to regions outside the 'axion window' due to a significant trade-off: the need to effectively suppress the magnetic leakage from highly polarized spin sources while simultaneously detecting sub-femtotesla level exotic physics signals at sub-decimeter-scale distances. In this work, we report new experimental results on axion-mediated exotic spin-spin interactions using an iron-shielded SmCo$_5$ spin source in combination with a specially designed self-compensation comagnetometer. Employing a composite shielding structure, we achieved a suppression of the magnetic field by up to $10^{11}$. This enabled us to establish new constraints on the coupling between electrons and neutrons, surpassing previous experimental limits by more than 10000 times within the axion window. Furthermore, we also set strongest constraints on the coupling between electrons and protons. The proposed method holds substantial potential not only for advancing the search for new physics beyond the Standard Model but also for enabling transformative applications in biological and chemical research.

Self-consistent mean-field methods with Skyrme-type effective interactions and semiclassical approximations, such as the Thomas-Fermi approach and its extensions are particularly well-suited for describing in a thermodynamically consistent way the various phases of the dense matter present in the interior of neutron stars. These methods have been applied to predict the composition of the different regions, including the inner crust constituted by nuclear clusters coexisting with free neutrons and electrons. Because of the computational cost, the energy is typically calculated for a few selected average baryon number densities, and the results are interpolated to obtain the pressure numerically. However, this may introduce systematic errors in the calculations of the global structure of a neutron star and its dynamical evolution. In this paper, we show how the full equation of state can be consistently calculated within the same framework by deriving exact formulas for the chemical potentials and for the pressure that can be easily implemented in existing computer codes. These formulas are applicable to both catalyzed and accreted crusts. We discuss in each case the suitable conditions to impose to determine the composition. Numerical examples are also presented and discussed. Results from refined calculations of the BSk24 equation of state for the inner crust of nonaccreted neutron stars and the corresponding adiabatic index are provided.

New limits on the axion-nucleon coupling over the axion mass region $10^{-24} \leq m_a \leq 5 \times 10^{-21}$ eV are derived by reanalyzing data from laboratory measurements on Lorentz and $CPT$ violation. These results establish the first laboratory constraints on the axion-nucleon coupling for axion masses below $10^{-22}$ eV. For $10^{-22} \leq m_a \leq 5 \times 10^{-21}$ eV, the results improve upon previous laboratory limits by more than 3 orders of magnitude, exceeding for the first time the astrophysical limits from supernova SN1987A cooling. For the axion mass range of interest corresponding to ultralow frequencies, the crucial initial phase of the axion field is considered. Furthermore, the obtained limits are nearly equivalent to those projected for a recently proposed experiment employing high-intensity neutron beams at the European Spallation Source [Phys. Rev. Lett. 133, 181001 (2024)]. For an alternative type of axion-nucleon interaction, the quadratic wind coupling, the constraints exceed the current best results by approximately 2 orders of magnitude.

Within the framework of Chern-Simons gravity, a theory that dynamically violates parity, we analyze the power spectrum of gravitational waves in light of the damping effect due to the free streaming relativistic neutrinos and dark fermions. The power spectrum is expressed terms of right- and left-handed polarizations, and the evolution of the gravitational waves is studied numerically. Birefringence is explicitly shown in the power spectrum, though the difference in the amplitudes is small. Specific features of peaks and dips appear gravitational wave power spectrum mirroring chiral gravitational wave mediated parametric resonance during reheating. Our result represents a useful tool to test Chern-Simons gravity and enables to constrain mechanisms of inflation and reheating related to this theory. We predict a falsifiable pattern of observable peaks and dips in the chiral independent gravitational power spectrum, eventually observable in next space-borne gravitational interferometers, including LISA, Taiji and Tianqin.

We conduct the first full Bayesian inference analysis for LISA parameter estimation incorporating the effects of subdominant harmonics and spin-precession through a full time domain response. The substantial computational demands of using time domain waveforms for LISA are significantly mitigated by implementing a novel Python version of the IMRPhenomT family of waveform models and the LISA response with GPU acceleration. This time domain response alleviates the theoretical necessity of developing specific transfer functions to approximate the LISA response in the Fourier domain for each specific type of system and allows for the use of unequal arms configurations and realistic LISA orbits. Our analysis includes a series of zero-noise injections for a Massive Black Hole Binary with aligned and precessing spins. We investigate the impact of including subdominant harmonics, compare equal and unequal arm configurations, and analyze different Time-Delay-Interferometry (TDI) configurations. We utilize full and uniform priors, with a lower frequency cutoff of 0.1mHz, and a signal duration of approximately two months, sampled every 5 seconds. The sampler is initialized based on Fisher estimates. Our results demonstrate LISA capability to measure the two spin magnitudes and the primary spin tilt angle, alongside sky localization, with percent-level precision, while component masses are determined with sub-percent accuracy.