Papers for Thursday, Mar 28 2024

Simulations of structure formation in the standard cold dark matter cosmological model quantify the dark matter halos of galaxies. Taking into account dynamical friction between the dark matter halos, we investigate the past orbital dynamical evolution of the Magellanic Clouds in the presence of the Galaxy. Our calculations are based on a three-body model of rigid Navarro-Frenk-White profiles for the dark matter halos, but were verified in a previous publication by comparison to high-resolution N-body simulations of live self-consistent systems. Under the requirement that the LMC and SMC had an encounter within 20 kpc between 1 and 4 Gyr ago, in order to allow the development of the Magellanic Stream, and using the latest astrometric data, the dynamical evolution of the MW/LMC/SMC system is calculated backwards in time. With the employment of the genetic algorithm and a Markov-Chain Monte-Carlo method, the present state of this system is unlikely with a probability of <10^{-9} (6 sigma complement), because solutions found do not fit into the error bars for the observed plane-of-sky velocity components of the Magellanic Clouds. This implies that orbital solutions that assume dark matter halos according to cosmological structure formation theory to exist around the Magellanic Clouds and the Milky Way are not possible with a confidence of more than 6 sigma

Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify the systematic errors that arise from the kinematic distance method in the presence of a Galactic potential that is non-axisymmetric. We investigate how these errors differ in certain regions of the Galaxy and how they relate to the underlying dynamics. We perform 2D isothermal hydrodynamical simulation of the gas disk with the moving-mesh code Arepo, adding the capability of using an external potential provided by the Agama library for galactic dynamics. We introduce a new analytic potential of the Milky Way, taking elements from existing models and adjusting parameters to match recent observational constraints. We find significant errors in the kinematic distance estimate for gas close to the Sun, along sight lines towards the Galactic centre and anti-centre, and significant deviations associated with the Galactic bar. Kinematic distance errors are low within the spiral arms as gas resides close to local potential minima and the resulting line-of-sight velocity is close to what is expected for an axisymmetric potential. Interarm regions exhibit large deviations at any given Galactic radius. This is caused by the gas being sped up or slowed down as it travels into or out of the spiral arm. We are able to define 'zones of avoidance' in the lv-diagram, where the kinematic distance method is particularly unreliable and should only be used with caution. We report a power law relation between the kinematic distance error and the deviation of the project line-of-sight velocity from circular motion.

We present the results of a search for ammonia maser emission in 119 Galactic high-mass star-forming regions (HMSFRs) known to host 22 GHz H$_2$O maser emission. Our survey has led to the discovery of non-metastable NH$_3$ inversion line masers toward 14 of these sources. This doubles the number of known non-metastable ammonia masers in our Galaxy, including nine new very high excitation ($J,K$)~=~(9,6) maser sources. These maser lines, including NH$_3$ (5,4), (6,4), (6,5), (7,6), (8,6), (9,6), (9,8), (10,8), and (11,9), arise from energy levels of 342 K, 513 K, 465 K, 606 K, 834 K, 1090 K, 942 K, 1226 K, and 1449 K above the ground state. Additionally, we tentatively report a new metastable NH$_3$ (3,3) maser in G048.49 and an NH$_3$ (7,7) maser in G029.95. Our observations reveal that all of the newly detected NH$_3$ maser lines exhibit either blueshifted or redshifted velocities with respect to the source systemic velocities. Among the non-metastable ammonia maser lines, larger velocity distributions, offset from the source systemic velocities, are found in the ortho-NH$_3$ ($K=3n$) than in the para-NH$_3$ ($K\neq3n$) transitions.

Spins play a crucial role in the appearance, evolution and occupation fraction of massive black holes. To date observational estimates of massive black hole spins are scarce, and the assumptions commonly made in such estimates have been recently questioned. Similarly, theoretical models for massive black hole spin evolution, while reproducing the few observational constraints, are based on possibly oversimplified assumptions. New independent constraints on massive black hole spins are therefore of primary importance. We present a rigorous statistical analysis of the relative orientation of radio jets and mega-maser discs in ten low-redshift galaxies. We find a strong preference for (partial) alignment between jets and mega-masers, that can be attributed to two different causes: "coherent" and "selective" accretion. In the first case the partial alignment is due to an anisotropy in the gas reservoir fueling the growth of MBHs. In the second case the spin-dependent anisotropic feedback allows for long-lived accretion only if the orbits of the gas inflows are almost aligned to the MBH equatorial plane. A discussion of the implications of the two accretion scenarios on the evolution of MBHs is presented, together with an outlook on future observational tests aiming at discriminating between the two scenarios and at checking if any of the two apply to different redshifts and black hole mass regimes.

Only indirect evidence of the role of magnetic braking in regulating gravitational collapse and the formation of circumstellar disks was found from observational work, such as compact disk sizes and the launching of high-velocity collimated protostellar jets. More direct tests of the magnetic braking shaping the angular momentum (AM) of the gas in Class 0 protostars are crucially needed. In the present work we have used non-ideal MHD models of protostellar collapse and synthetic observations of molecular gas spectral emission that we analyze to test whether possible kinematic signatures of the magnetic braking in the gas velocity field can be captured from maps of the molecular gas emission in protostellar envelopes. By comparing the 3D Specific AM of models with varying turbulent energy and magnetization, we show that, in the numerical models of protostellar evolution explored, the increase in magnetization and its consequences on the spatial redistribution of SAM modifies the shapes of the radial profiles of SAM. We show that widely used observational methods fail to quantitatively capture the magnitude of SAM of the gas in protostellar envelopes, and that no method allows to measure the differences in radial evolution of SAM due to different magnetization at all envelope radii. This is especially true in the more magnetized cases. However, our analysis suggests that the detection of symmetric patterns and organized velocity fields, in the moment-1 maps of the molecular line emission as well as monotonous radial profiles of the SAM showing a power-law decline, should be suggestive of a less magnetized scenario. Protostellar cores where efficient magnetic braking is at work are more likely to present a highly asymmetric velocity field, and more prone to show complex radial profiles of their specific angular momentum measured in the equatorial plane.

PDS 456 is the most luminous RQQ at z<0.3 and can be regarded as a local counterpart of the powerful QSOs shining at Cosmic Noon. It hosts a strong nuclear X-ray ultra-fast outflow, and a massive and clumpy CO(3-2) molecular outflow extending up to 5 kpc from the nucleus. We analyzed the first MUSE WFM and AO-NFM optical integral field spectroscopic observations of PDS456. The AO-NFM observations provide an unprecedented spatial resolution, reaching up to 280 pc. Our findings reveal a complex circumgalactic medium around PDS 456, extending up to a maximum projected size of ~46 kpc. This includes a reservoir of gas with a mass of ~1e7-1e8 Modot, along with eight companion galaxies, and a multi-phase outflow. WFM and NFM MUSE data reveal an outflow on a large scale (~12 kpc from the quasar) in [OIII], and on smaller scales (within 3 kpc) with higher resolution (about 280 pc) in Halpha, respectively. The [OIII] outflow mass rate is 2.3 +/- 0.2 Modot/yr which is significantly lower than those typically found in other luminous quasars. Remarkably, the Ha outflow shows a similar scale, morphology, and kinematics to the CO(3-2) molecular outflow, with the latter dominating in terms of kinetic energy and mass outflow rate by two and one orders of magnitude, respectively. Our results therefore indicate that mergers, powerful AGN activity, and feedback through AGN-driven winds will collectively contribute to shaping the host galaxy evolution of PDS 456, and likely, that of similar objects at the brightest end of the AGN luminosity function across all redshifts. Moreover, the finding that the momentum boost of the total outflow deviates from the expected energy-conserving expansion for large-scale outflows highlights the need of novel AGN-driven outflow models to comprehensively interpret these phenomena.

The BL Lac object 3C 371 has been observed by the Transiting Exoplanet Survey Satellite (\textit{TESS}) for approximately a year, between July 2019 and July 2020, with an unmatched 2-minute observing cadence. In parallel, the Whole Earth Blazar Telescope (WEBT) Collaboration organized an extensive observing campaign, providing three years of continuous optical monitoring between 2018 and 2020. These datasets allow for a thorough investigation of the variability of the source. The goal of this study is evaluating the optical variability of 3C 371. Taking advantage of the remarkable cadence of \textit{TESS} data, we aim to characterize the intra-day variability (IDV) displayed by the source and identify its shortest variability timescale. With this estimate, constraints on the size of the emitting region and black hole mass can be calculated. Moreover, WEBT data is used to investigate long-term variability (LTV), including understanding spectral behaviour of the source and the polarization variability. Based on the derived characteristics, information on the origin of the variability on different timescales is extracted. We evaluated the variability applying the variability amplitude tool that quantifies how variable the emission is. Moreover, we employed common tools like ANOVA (ANalysis Of VAariance) tests, wavelet and power spectral density (PSD) analyses to characterize the shortest variability timescales present in the emission and the underlying noise affecting the data. Short- and long-term colour behaviours have been evaluated to understand the spectral behaviour of the source. The polarized emission was analyzed, studying its variability and possible rotation patterns of the electric vector position angle (EVPA). Flux distributions of IDV and LTV were also studied with the aim of linking the flux variations to turbulent and/or accretion disc related processes.

The ``Residual-to-Residual DNN series for high-Dynamic range imaging'' (R2D2) approach was recently introduced for Radio-Interferometric (RI) imaging in astronomy. R2D2's reconstruction is formed as a series of residual images, iteratively estimated as outputs of Deep Neural Networks (DNNs) taking the previous iteration's image estimate and associated data residual as inputs. In this work, we investigate the robustness of the R2D2 image estimation process, by studying the uncertainty associated with its series of learned models. Adopting an ensemble averaging approach, multiple series can be trained, arising from different random DNN initializations of the training process at each iteration. The resulting multiple R2D2 instances can also be leveraged to generate ``R2D2 samples'', from which empirical mean and standard deviation endow the algorithm with a joint estimation and uncertainty quantification functionality. Focusing on RI imaging, and adopting a telescope-specific approach, multiple R2D2 instances were trained to encompass the most general observation setting of the Very Large Array (VLA). Simulations and real-data experiments confirm that: (i) R2D2's image estimation capability is superior to that of the state-of-the-art algorithms; (ii) its ultra-fast reconstruction capability (arising from series with only few DNNs) makes the computation of multiple reconstruction samples and of uncertainty maps practical even at large image dimension; (iii) it is characterized by a very low model uncertainty.

We present the semi-analytical light curve modelling of 13 supernovae associated with gamma-ray bursts (GRB-SNe) along with two relativistic broad-lined (Ic-BL) SNe without GRBs association (SNe 2009bb and 2012ap), considering millisecond magnetars as central-engine-based power sources for these events. The bolometric light curves of all 15 SNe in our sample are well-regenerated utilising a $\chi^2-$minimisation code, $\texttt{MINIM}$, and numerous parameters are constrained. The median values of ejecta mass ($M_{\textrm{ej}}$), magnetar's initial spin period ($P_\textrm{i}$) and magnetic field ($B$) for GRB-SNe are determined to be $\approx$ 5.2 M$_\odot$, 20.5 ms and 20.1 $\times$ 10$^{14}$ G, respectively. We leverage machine learning (ML) algorithms to comprehensively compare the 3-dimensional parameter space encompassing $M_{\textrm{ej}}$, $P_\textrm{i}$, and $B$ for GRB-SNe determined herein to those of H-deficient superluminous SNe (SLSNe-I), fast blue optical transients (FBOTs), long GRBs (LGRBs), and short GRBs (SGRBs) obtained from the literature. The application of unsupervised ML clustering algorithms on the parameters $M_{\textrm{ej}}$, $P_\textrm{i}$, and $B$ for GRB-SNe, SLSNe-I, and FBOTs yields a classification accuracy of $\sim$95%. Extending these methods to classify GRB-SNe, SLSNe-I, LGRBs, and SGRBs based on $P_\textrm{i}$ and $B$ values results in an accuracy of $\sim$84%. Our investigations show that GRB-SNe and relativistic Ic-BL SNe presented in this study occupy different parameter spaces for $M_{\textrm{ej}}$, $P_\textrm{i}$, and $B$ than those of SLSNe-I, FBOTs, LGRBs and SGRBs. This indicates that magnetars with different $P_\textrm{i}$ and $B$ can give birth to distinct types of transients.

Galaxy clusters are the largest gravitationally-bound structures in the Universe. As such, during merger events with similar systems, they release an enormous amount of energy that is dissipated through the formation of shock waves and turbulence in the intracluster medium (ICM), the hot ionised plasma permeating the cluster volume. These shock waves are believed to be ideal sites for electron acceleration that, in the presence of ubiquitous magnetic fields in the ICM, are capable of producing elongated non-thermal radio features typically observed in the outskirts of dynamically disturbed clusters, also known as radio relics. In this work, we analyse a hydrodynamical re-simulation of merging galaxy clusters, extracted from a large set of zoom-in cosmological simulations of The Three Hundred Project, to study the evolution and diversity of merger shocks and their associated diffuse radio emission within the framework of the diffusive shock acceleration scenario.

The origin of the chemical diversity observed around low-mass protostars probably resides in the earliest history of these systems. We aim to investigate the impact of protostellar feedback on the chemistry and grain growth in the circumstellar medium of multiple stellar systems. In the context of the ALMA Large Program FAUST, we present high-resolution (50 au) observations of CH$_3$OH, H$_2$CO, and SiO and continuum emission at 1.3 mm and 3 mm towards the Corona Australis star cluster. Methanol emission reveals an arc-like structure at $\sim$1800 au from the protostellar system IRS7B along the direction perpendicular to the major axis of the disc. The arc is located at the edge of two elongated continuum structures that define a cone emerging from IRS7B. The region inside the cone is probed by H$_2$CO, while the eastern wall of the arc shows bright emission in SiO, a typical shock tracer. Taking into account the association with a previously detected radio jet imaged with JVLA at 6 cm, the molecular arc reveals for the first time a bow shock driven by IRS7B and a two-sided dust cavity opened by the mass-loss process. For each cavity wall, we derive an average H$_2$ column density of $\sim$7$\times$10$^{21}$ cm$^{-2}$, a mass of $\sim$9$\times$10$^{-3}$ M$_\odot$, and a lower limit on the dust spectral index of $1.4$. These observations provide the first evidence of a shock and a conical dust cavity opened by the jet driven by IRS7B, with important implications for the chemical enrichment and grain growth in the envelope of Solar System analogues.

Observations of cosmic microwave background polarisation, essential for probing a potential phase of inflation in the early universe, suffer from contamination by polarised emission from the Galactic interstellar medium. This work combines existing observations from the WMAP and Planck space missions to make a low-noise map of polarised synchrotron emission that can be used to clean forthcoming CMB observations. We combine WMAP K, Ka and Q maps with Planck LFI 30~GHz and 44~GHz maps using weights that near-optimally combine the observations as a function of sky direction, angular scale, and polarisation orientation. We publish well-characterised maps of synchrotron Q and U Stokes parameters at nu = 30GHz and 1 degree angular resolution. A statistical description of uncertainties is provided with Monte-Carlo simulations of additive and multiplicative errors. Our maps are the most sensitive full-sky maps of synchrotron polarisation to date, and are made available to the scientific community on a dedicated web site.

Relativistic jets are observed from accreting and cataclysmic transients throughout the Universe, and have a profound affect on their surroundings. Despite their importance, their launch mechanism is not known. For accreting neutron stars, the speed of their compact jets can reveal whether the jets are powered by magnetic fields anchored in the accretion flow or in the star itself, but to-date no such measurements exist. These objects can display bright explosions on their surface due to unstable thermonuclear burning of recently accreted material, called type-I X-ray bursts, during which the mass accretion rate increases. Here, we report on bright flares in the jet emission for a few minutes after each X-ray burst, attributed to the increased accretion rate. With these flares, we measure the speed of a neutron star compact jet to be $v=0.38^{+0.11}_{-0.08}$c, much slower than those from black holes at similar luminosities. This discovery provides a powerful new tool in which we can determine the role that individual system properties have on the jet speed, revealing the dominant jet launching mechanism.

The main sequence (MS) of star-forming galaxies (SFGs) is the tight relation between the galaxy stellar mass and its star formation rate (SFR) and was observed up to z ~ 6. The MS relation can be used as a reference for understanding the differences among galaxies, characterised by different rates of stellar production (starbursts, SFGs, and passive galaxies), and those inside a galaxy made up of different components (bulge, disk, and halo). To investigate peculiar features found in our sample galaxies, we focus here on their star formation history (SFH). We performed a spectral energy distribution fitting procedure that accounted for the energetic balance between UV and far-IR radiation on a sample of eight nearby face-on spiral galaxies from the DustPedia sample. This approach allowed us to study the spatially resolved MS of the sample and to recover the past SFH. By exploiting the BAGPIPES code, we constrained the SFHs for each galaxy with a delayed exponentially declining model to derive their mass-weighted age (tMW). A central old region (tMW up to~7Gyr, consistent with the presence of a bulge for various systems) is followed by younger regions in which the disks are still forming stars (tMW~4Gyr). At larger distances, tMW increases mildly in general. Strikingly, in two galaxies (NGC4321 and NGC5194), we found a steep increase in tMW that reached levels similar to those of the bulge. These old stellar populations in the very galaxy outskirts are unexpected. We discuss their potential origin by considering the different gas phases of the source with the most prominent quenched ring, NGC4321, and argue for two main possibilities: 1) some environmental effect (e.g. starvation) or 2) the circumgalactic medium of sources outside of high-density clusters might have stopped to supply pristine gas to the galaxy (e.g. if its specific angular moment is too high for being accreted).

Recent studies indicate that some Galactic open clusters (OCs) exhibit extended main-sequence turnoff (eMSTO) in their colour-magnitude diagrams (CMDs). However, the number of Galactic OCs with eMSTO structures detected so far is limited, and the reasons for their formation are still unclear. This work identifies 26 Galactic OCs with undiscovered eMSTOs and investigates the causes of these features. Stellar population types and fundamental parameters of cluster samples are acquired using CMD fitting methods. Among them, the results of 11 OCs are reliable as the observed CMDs are well-reproduced. We propose the crucial role of stellar binarity and confirm the importance of stellar rotation in reproducing eMSTO morphologies. The results also show that the impact of age spread is important, as it can adequately explain the structure of young OCs and fit the observed CMDs of intermediate-age OCs better.

While substantial progress has been made in studying the pre-reconstruction galaxy bispectrum, investigations of the post-reconstruction bispectrum are still in nascent stages. In this paper, we present a bispectrum model that incorporates one-loop corrections in the Standard Perturbation Theory (SPT), while simultaneously addressing infrared (IR) effects in the post-reconstruction density fluctuations in a non-perturbative approach. This model accurately captures the nonlinear behavior of the Baryon Acoustic Oscillation (BAO) signal within the post-reconstruction bispectrum framework. Furthermore, the incorporation of the one-loop correction extends its applicability to smaller scales. Given that the approach to addressing IR effects remains invariant before and after the reconstruction process, the resulting IR-resummed bispectrum model exhibits uniformity across both scenarios. Throughout this analysis, we achieve a comprehensive and unified description of the bispectrum that seamlessly spans the pre- and post-reconstruction phases. An explanatory video is available at https://youtu.be/bHdlIi3jQ4A.

In this study, we explore the characteristics of carbon stars within our Galaxy through a comprehensive analysis of observational data spanning visual and infrared (IR) bands. Leveraging datasets from IRAS, ISO, Akari, MSX, 2MASS, WISE, Gaia DR3, AAVSO, and the SIMBAD object database, we conduct a detailed comparison between the observational data and theoretical models. To facilitate this comparison, we introduce various IR two-color diagrams (2CDs), IR color-magnitude diagrams (CMDs), and spectral energy distributions (SEDs). We find that the CMDs, which utilize the latest distance and extinction data from Gaia DR3 for a substantial number of carbon stars, are very useful to distinguish carbon-rich asymptotic giant branch (CAGB) stars from extrinsic carbon stars that are not in the AGB phase. To enhance the accuracy of our analysis, we employ theoretical radiative transfer models for dust shells around CAGB stars. These theoretical dust shell models demonstrate a commendable ability to approximate the observations of CAGB stars across various SEDs, 2CDs, and CMDs. We present the infrared properties of known pulsating variables and explore the infrared variability of the sample stars by analyzing WISE photometric data spanning the last 14 yr. Additionally, we present a novel catalog of CAGB stars, offering enhanced reliability and a wealth of additional information.

Recent high-angular resolution ALMA observations have revealed rich information about protoplanetary disks, including ubiquitous substructures and three-dimensional gas kinematics at different emission layers. One interpretation of these observations is embedded planets. Previous 3-D planet-disk interaction studies are either based on viscous simulations, or non-ideal magnetohydrodynamics (MHD) simulations with simple prescribed magnetic diffusivities. This study investigates the dynamics of gap formation in 3-D non-ideal MHD disks using non-ideal MHD coefficients from the look-up table that is self-consistently calculated based on the thermo-chemical code. We find a concentration of the poloidal magnetic flux in the planet-opened gap (in agreement with previous work) and enhanced field-matter coupling due to gas depletion, which together enable efficient magnetic braking of the gap material, driving a fast accretion layer significantly displaced from the disk midplane. The fast accretion helps deplete the gap further and is expected to negatively impact the growth of planetary embryos. It also affects the corotation torque by shrinking the region of horseshoe orbits on the trailing side of the planet. Together with the magnetically driven disk wind, the fast accretion layer generates a large, persistent meridional vortex in the gap, which breaks the mirror symmetry of gas kinematics between the top and bottom disk surfaces. Finally, by studying the kinematics at the emission surfaces, we discuss the implications of planets in realistic non-ideal MHD disks on kinematics observations.

Image super-resolution has been an important subject in image processing and recognition. Here, we present an attention-aided convolutional neural network (CNN) for solar image super-resolution. Our method, named SolarCNN, aims to enhance the quality of line-of-sight (LOS) magnetograms of solar active regions (ARs) collected by the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO). The ground-truth labels used for training SolarCNN are the LOS magnetograms collected by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). Solar ARs consist of strong magnetic fields in which magnetic energy can suddenly be released to produce extreme space weather events, such as solar flares, coronal mass ejections, and solar energetic particles. SOHO/MDI covers Solar Cycle 23, which is stronger with more eruptive events than Cycle 24. Enhanced SOHO/MDI magnetograms allow for better understanding and forecasting of violent events of space weather. Experimental results show that SolarCNN improves the quality of SOHO/MDI magnetograms in terms of the structural similarity index measure (SSIM), Pearson's correlation coefficient (PCC), and the peak signal-to-noise ratio (PSNR).

We present a progress report of our project aiming to increase the number of known Cepheids in double-lined binary (SB2) systems from six to 100 or more. This will allow us, among other goals, to accurately measure masses for a large sample of Cepheids. Currently, only six accurate Cepheid masses are available, which hinders our understanding of their physical properties and renders the Cepheid mass--luminosity relation poorly constrained. At the same time, Cepheids are widely used for essential measurements (e.g., extragalactic distances, the Hubble constant). To examine Cepheid period--luminosity relations, we selected as binary candidates Cepheids that are too bright for their periods. To date, we have confirmed 56 SB2 systems, including the detection of significant orbital motions of the components for 32. We identified systems with orbital periods up to five times shorter than the shortest reported period to date, as well as systems with mass ratios significantly different from unity (suggesting past merger events). Both features are essential to understand how multiplicity affects the formation and destruction of Cepheid progenitors and what effect this has on global Cepheid properties. We also present eight new systems composed of two Cepheids (only one such system was known before). Among confirmed SB2 Cepheids, there are also several wide-orbit systems. In the future, these may facilitate independent accurate geometric distance measurements to the Large and Small Magellanic Clouds.

The detection and analysis of the solar coronal holes (CHs) is an important field of study in the domain of solar physics. Mainly, it is required for the proper prediction of the geomagnetic storms which directly or indirectly affect various space and ground-based systems. For the detection of CHs till date, the solar scientist depends on manual hand-drawn approaches. However, with the advancement of image processing technologies, some automated image segmentation methods have been used for the detection of CHs. In-spite of this, fast and accurate detection of CHs are till a major issues. Here in this work, a novel quantum computing-based fast fuzzy c-mean technique has been developed for fast detection of the CHs region. The task has been carried out in two stages, in first stage the solar image has been segmented using a quantum computing based fast fuzzy c-mean (QCFFCM) and in the later stage the CHs has been extracted out from the segmented image based on image morphological operation. In the work, quantum computing has been used to optimize the cost function of the fast fuzzy c-mean (FFCM) algorithm, where quantum approximate optimization algorithm (QAOA) has been used to optimize the quadratic part of the cost function. The proposed method has been tested for 193 \AA{} SDO/AIA full-disk solar image datasets and has been compared with the existing techniques. The outcome shows the comparable performance of the proposed method with the existing one within a very lesser time.

Aims. Solar active regions (ARs), which are formed by flux emergence, serve as the primary sources of solar eruptions. However, the specific physical mechanism that governs the emergence process and its relationship with flare productivity remains to be thoroughly understood. Methods. We examined 136 emerging ARs, focusing on the evolution of their magnetic helicity and magnetic energy during the emergence phase. Based on the relation between helicity accumulation and magnetic flux evolution, we categorized the samples and investigated their flare productivity. Results. The emerging ARs we studied can be categorized into three types, Type-I, Type-II, and Type-III, and they account for 52.2%, 25%, and 22.8% of the total number in our sample, respectively. Type-I ARs exhibit a synchronous increase in both the magnetic flux and magnetic helicity, while the magnetic helicity in Type-II ARs displays a lag in increasing behind the magnetic flux. Type-III ARs show obvious helicity injections of opposite signs. Significantly, 90% of the flare-productive ARs (flare index > 6) were identified as Type-I ARs, suggesting that this type of AR has a higher potential to become flare productive. In contrast, Type-II and Type-III ARs exhibited a low and moderate likelihood of becoming active, respectively. Our statistical analysis also revealed that Type-I ARs accumulate more magnetic helicity and energy, far beyond what is found in Type-II and Type-III ARs. Moreover, we observed that flare-productive ARs consistently accumulate a significant amount of helicity and energy during their emergence phase. Conclusions. These findings provide valuable insight into the flux emergence phenomena, offering promising possibilities for early-stage predictions of solar eruptions.

An essential factor in determining the flow characteristics of an accretion flow is its angular momentum. According to the angular momentum of the flow, semi-analytical analysis suggests various types of accretion solutions. It is critical to test it with numerical simulations using the most advanced framework available (general relativistic magnetohydrodynamics) to understand how flow changes with different angular momentum. By changing the initial condition of the accretion torus minimally, we can simulate steady, low angular momentum accretion flow around a Kerr black hole. We focus primarily on the lower limits of angular momentum and come upon that accretion flow with an intermediate range of angular momentum differs significantly from high or very low angular momentum flow. The intermediate angular momentum accretion flow has the highest density, pressure, and temperature near the black hole, making it easier to observe. We find that the density and pressure have power-law scalings $\rho\propto r^{n-3/2}$ and $p_g\propto r^{n-5/2}$ which only hold for very low angular momentum cases. With the increase in flow angular momentum, it develops a non-axisymmetric nature. In this case, simple self-similarity does not hold. We also find that the sonic surface moves away from the innermost stable circular orbit as its angular momentum decreases. Finally, we emphasize that intermediate angular momentum flow could provide a possible solution to explain the complex observation features of the supermassive black hole Sgr~A$^*$ at our galactic center.

At present, there is no doubt that relativistic jets observed in active galactic nuclei pass from highly magnetized to weakly magnetized stage, which is observed as a break in the dependence on their width $d_{\rm jet}(z)$ on the distance $z$ to the central engine. In this paper, we discuss the possibility of observing another break, which should be located at shorter distances. The position of this break can be associated with the region of formation of the dense central core near the jet axis which was predicted both analytically and numerically more than a decade ago, but has not yet received sufficient attention. In this case, the observed width should be determined by the dense core, and not by the total transverse size of the jet. The calculations carried out in this paper, which took into account both the transverse electromagnetic structure of the jet and the change in the spectrum of emitting particles along its axis, indeed showed such behaviour. We also found the evidence of the predicted break in the jet expansion profile using stacked 15 GHz VLBA image of M87 radio jet and constrain the light cylinder radius.

We investigate non-ideal magnetohydrodynamical (MHD) effects in the chromosphere on the solar wind by performing MHD simulations for Alfv\'en-wave driven winds with explicitly including Ohmic and ambipolar diffusion. We find that MHD waves are significantly damped in the chromosphere by ambipolar diffusion so that the Alfv\'enic Poynting flux that reaches the corona is substantially reduced. As a result, the coronal temperature and the mass loss rate of the solar wind are considerably reduced, compared with those obtained from an ideal MHD case, which is indicative of a great importance of the non-ideal MHD effects in the solar atmosphere. However, the temperature and the mass loss rate are recovered by a small increase in the convection-originated velocity perturbation at the photosphere because of the sensitive dependence of the ambipolar diffusion and reflection of Alfv\'en waves on the physical properties of the chromosphere. We also find that density perturbations in the corona are reduced by the ambipolar diffusion of Alfv\'en waves in the chromosphere because the nonlinear generation of compressible perturbations is suppressed.

Context. Fully fledged penumbrae have been widely studied both observationally and theoretically. Yet the relatively fast process of penumbra formation has not been studied closely with high spatial resolution. Aims. We investigate the stages previous to and during the formation of penumbral filaments in a developing sunspot. Methods. We analysed Milne-Eddington inversions from spectro-polarimetric data of the leading sunspot of NOAA 11024 during the development of its penumbra. We focused on selected areas of this protospot in which segments of penumbra develop. Results. We find that few types of distinctive flow patterns develop at the protospot limb and centre sides previous to penumbra formation. The flow in the centre side is often characterised by a persistent (>20 min) inflow-outflow pattern extending radially over 4 arcsec at the direct periphery of the protospot umbra. This inflow-outflow system often correlates with elongated granules, as seen in continuum intensity maps, and is also coupled with magnetic bipolar patches at its edges, as seen in magnetograms. The field is close to horizontal between the bipolar patches, which is indicative of its possible loop configuration. All of these aspects are analogous to observations of magnetic flux emergence. In the protospot limb side, however, we observed a mostly regular pattern associated with small granules located near the protospot intensity boundary. Locally, an inflow develops adjacent to an existing penumbral segment, and this inflow is correlated with a single bright penumbral filament that is brighter than filaments containing the Evershed flow. All investigated areas at the centre and limb side eventually develop penumbral filaments with an actual Evershed flow that starts at the umbral boundary and grows outwards radially as the penumbral filaments become longer in time

We present the results of a machine learning study to measure the dust content of galaxies observed with JWST at z > 6 through the use of trained neural networks based on high-resolution IllustrisTNG simulations. Dust is an important unknown in the evolution and observability of distant galaxies and is degenerate with other stellar population features through spectral energy fitting. As such, we develop and test a new SED-independent machine learning method to predict dust attenuation and sSFR of high redshift (z > 6) galaxies. Simulated galaxies were constructed using the IllustrisTNG model, with a variety of dust contents parameterized by E(B-V) and A(V) values, then used to train Convolutional Neural Network (CNN) models using supervised learning through a regression model. We demonstrate that within the context of these simulations, our single and multi-band models are able to predict dust content of distant galaxies to within a 1$\sigma$ dispersion of A(V) $\sim 0.1$. Applied to spectroscopically confirmed z > 6 galaxies from the JADES and CEERS programs, our models predicted attenuation values of A(V) < 0.7 for all systems, with a low average (A(V) = 0.28). Our CNN predictions show larger dust attenuation but lower amounts of star formation compared to SED fitted values. Both results show that distant galaxies with confirmed spectroscopy are not extremely dusty, although this sample is potentially significantly biased. We discuss these issues and present ideas on how to accurately measure dust features at the highest redshifts using a combination of machine learning and SED fitting.

The Fermi and eROSITA bubbles, large diffuse structures in our Galaxy, can be the by-products of the steady star formation activity. To simultaneously explain the star formation history of the Milky Way and the metallicity of $\sim$ Z$_\odot$ at the Galactic disk, a steady Galactic wind driven by cosmic-rays is required. For tenuous gases with a density of $\lesssim$10$^{-3}$ cm$^{-3}$, the cosmic-ray heating dominates over radiative cooling, and the gas can maintain the virial temperature of $\sim$0.3 keV ideal for escape from the Galactic system as the wind. A part of the wind falls back onto the disk like a galactic fountain flow. We model the wind dynamics according to the Galactic evolution scenario and find that the scale height and surface brightness of the X-ray and the hadronic gamma-ray emissions from such fountain flow region can be consistent with the observed properties of the Fermi and eROSITA bubbles. This implies that the bubbles are persistent structures of the Milky Way existing over (at least) the last $\sim$1 Gyr, rather than evanescent structures formed by non-trivial, $\sim$10 Myr past Galactic Center transient activities.

Free-floating planets are a new class of planets recently discovered. These planets don't orbit within stellar systems, instead living a nomadic life within the galaxy. How such objects formed remains elusive. Numerous works have explored mechanisms to form such objects, but have not yet provided predictions on their distributions that could differentiate between formation mechanisms. In this work we form these objects within circumbinary systems, where these planets are readily formed and ejected through interactions with the central binary stars. We find significant differences between planets ejected through planet-planet interactions and those by the binary stars. The main differences that arise are in the distributions of excess velocity, where binary stars eject planets with faster velocities. These differences should be observable amongst known free-floating planets in nearby star-forming regions. We predict that targeted observations of directly imaged free-floating planets in these regions should be able to determine their preferred formation pathway, either by planet formation in single or multiple stellar systems, or through processes akin to star formation. Additionally the mass distributions of free-floating planets can yield important insights into the underlying planet populations. We find that for planets more massive than 20 $\rm M_{\oplus}$, their frequencies are similar to those planets remaining bound and orbiting near the central binaries. This similarity allows for effective and informative comparisons between mass distributions from microlensing surveys, to those of transit and radial velocities. Ultimately, by observing the velocity dispersion and mass distribution of free-floating planets, it will be possible to effectively compare with predictions from planet formation models, and to further understand the formation and evolution of these exotic worlds.

Cosmic rays are highly energetic messengers propagating in magnetized plasma, which are, possibly but not exclusively, accelerated at astrophysical shocks. Amongst the variety of astrophysical objects presenting shocks, the huge circumstellar stellar wind bubbles forming around very massive stars, are potential non thermal emitters. We present the 1D magnetohydrodynamical simulation of the evolving magnetized surroundings of a single, OB type main sequence 60 Mo star, which is post processed to calculate the re-acceleration of preexisting non-thermal particles of the Galactic cosmic ray background. It is found that the forward shock of such circumstellar bubble can, during the early phase (1 Myr) of its expansion, act as a substantial reaccelerator of pre existing interstellar cosmic rays. This results in an increasing excess emission flux by a factor of 5, the hadronic component producing gamma-rays by pion0 decay being more important than those by synchrotron and inverse Compton radiation mechanisms. We propose that this effect is at work in the circumstellar environments of massive stars in general and we conjecture that other nebulae such as the stellar wind bow shocks of runaway massive stars also act as Galactic cosmic-ray re-accelerators. Particularly, this study supports the interpretation of the enhanced hadronic emission flux measured from the surroundings of kappa Ori as originating from the acceleration of pre-existing particles at the forward shock of its wind bubble.

This research note concerns the application of deep-learning-based multi-view-imaging techniques to data from the H.E.S.S. Imaging Atmospheric Cherenkov Telescope array. We find that the earlier the fusion of layer information from different views takes place in the neural network, the better our model performs with this data. Our analysis shows that the point in the network where the information from the different views is combined is far more important for the model performance than the method used to combine the information.

Dark matter (DM) constitutes the predominant portion of matter in our universe. Despite compelling evidence, the precise characteristics of DM remain elusive. Among the leading DM candidates are weakly interacting massive particles, which may aggregate into steep concentrations around the central black holes of galaxies, forming dense spikes. Employing Schwarzschild geometry, we assess the density and velocity distribution of DM within such spikes. Through variations in black hole masses and dark halo parameters, we identify universal features in the density profile of DM and fit them with Gaussian distributions. Additionally, we investigate the impact of dynamical friction on gravitational waves (GWs) generated by extreme-mass-ratio inspirals (EMRIs) within DM spikes. Our findings uncover phase shifts in the time-domain waveform, potentially providing significant insights for the GW-based detection of DM in galactic centers.

Future observations of the Sunyaev-Zeldovich (SZ) effect promise ever improving measurements in terms of both sensitivity and angular resolution. As such, it is increasingly relevant to model `higher-order' contributions to the SZ effect. This work examines the effects of high-energy non-thermal electron distributions and those of anisotropic electron and photon distributions on the SZ signals. Analytic forms of the anisotropic scattering kernels for photons and electrons have been derived and investigated. We present a method for determining the anisotropic contributions through a spherical harmonic decomposition to arbitrary angular multipoles, and discuss the behaviour of these scattering kernels. We then carry out an exploration of various simplistic models of high energy non-thermal electron distributions, and examine their anisotropic behaviour. The kinematic SZ in the relativistic regime is studied using the kernel formulation allowing us to clarifying the role of kinematic corrections to the scattering optical depth. We finally present a release of an updated and refined version of SZpack including a new integrated Python interface and new modules for the calculation of various SZ signals, including those described in this paper.

Red (S > 10%/0.1 micron) spectral slopes are common among Centaurs and trans-Neptunian objects (TNOs) in the outer solar system. Interior to and co-orbital with Jupiter, the red (S approx. = 10%/0.1 micron) slopes of D-type main-belt and Jupiter Trojan asteroids are thought to reflect their hypothesized shared origin with TNOs beyond the orbit of Jupiter. In order to quantify the abundance of red-sloped asteroids within the main belt, we conducted a survey using the NASA Infrared Telescope Facility and the Lowell Discovery Telescope. We followed up on 32 candidate red objects identified via spectrophotometry from the Sloan Digital Sky Survey's Moving Object Catalog to confirm their steep spectral slopes and determine their taxonomic classifications. We find that our criteria for identifying candidate red objects from the Moving Object Catalog result in an approx. 50% confirmation rate for steeply red-sloped asteroids. We also compare our observations of main-belt asteroids to existing literature spectra of the Jupiter Trojans and steeply red-sloped main-belt asteroids. We show that some red-sloped asteroids have linearly increasing reflectance with increasing wavelength, while other red-sloped asteroids show a flattening in slope at longer near-infrared wavelengths, indicating a diversity among the population of spectrally red main-belt asteroids suggestive of a variety of origins among the population of steep-sloped asteroids.

Aims: we propose that the condition of relative motion between us and the objects that we observe in the Universe should generate relativistic aberration on the photons that such objects emit, varying the observed flux similarly to the cases of blazars and neutron stars, with instead a decrease of this radiative flux. Methods: we follow some important papers and textbook used in modern cosmology, adding the effects of relativistic aberration in the theoretical description of the Cosmos. Results: New definitions for the luminosity distance and the angular distance arise, with the consequence of changing the cosmological parameters measured in the last years. A qualitative description of all of this is also performed on the multipole analysis of the CMB.

The phenomenon of the Gnevyshev gap was first identified in the solar-corona irradiance data (green line). Later, it was studied in the sunspot, coronal, and heliospheric data. We have investigated the Gnevyshev gap in the magnetic field data and have arrived at the conclusion that it reflects the behavior of the large-scale magnetic field. The Gnevyshev gap occurs at the polarity reversal of the solar magnetic field at the photosphere level. The presence of the Gnevyshev gap in sunspot data at the photosphere level is disguised by non-global structures that retain dependence on both latitude and longitude (the accepted mathematical term is tessaral, see below for more details). However, it is clearly visible in the magnetic field data at the photosphere level and is even more pronounced at the source surface (i.e., in the corona).

Multiple nations are planning activity on the Moon's surface, and to deconflict lunar operations we must understand the sandblasting damage from rocket exhaust blowing soil. Prior research disagreed over the scaling of the erosion rate, which determines the magnitude of the damage. Reduced gravity experiments and two other lines of evidence now indicate that the erosion rate scales with the kinetic energy flux at the bottom of the laminar sublayer of the gas. Because the rocket exhaust is so fast, eroded particles lifted higher in the boundary layer do not impact the surface for kilometers (if at all; some leave the Moon entirely), so there is no saltation in the vicinity of the gas. As a result, there is little transport of gas kinetic energy from higher in the boundary layer down to the surface, so the emission of soil into the gas is a surprisingly low energy process. In low lunar gravity, a dominant source of resistance to this small energy flux turns out to be the cohesive energy density of the lunar soil, which arises primarily from particles in the 0.3 to 3 micron size range. These particles constitute only a tiny fraction of the mass of lunar soil and have been largely ignored in most studies, so they are poorly characterized.

In the companion paper ("Erosion rate of lunar soil under a landing rocket, part 1: identifying the rate-limiting physics", this issue) an equation was developed for the rate that lunar soil erodes under the exhaust of a landing rocket. That equation has only one parameter that is not calibrated from first principles, so here it is calibrated by the blowing soil's optical density curve during an Apollo landing. An excellent fit is obtained, helping validate the equation. However, when extrapolating the erosion rate all the way to touchdown on the lunar surface, a soil model is needed to handle the increased resistance to erosion as the deeper, more compacted soil is exposed. Relying on models derived from Apollo measurements and from Lunar Reconnaissance Orbiter (LRO) Diviner thermal inertia measurements, only one additional soil parameter is unknown: the scale of increasing cohesive energy with soil compaction. Treating this as an additional fitting parameter results in some degeneracy in the solutions, but the depth of erosion scour in the post-landing imagery provides an additional constraint on the solution. The results show that about 4 to 10 times more soil was blown in each Apollo landing than previously believed, so the potential for sandblasting damage is worse than prior estimates. This also shows that, with further development, instruments to measure the soil erosion during lunar landings can constrain the soil column's density profile complementary to the thermal inertia measurements, providing insight into the landing site's geology.

We discuss an algorithm whereby the massive galaxy clusters detected in the SRG/eROSITA all-sky survey are identified and their photometric redshifts are estimated. For this purpose, we use photometric redshift estimates for galaxies and WISE forced photometry. To estimate the algorithm operation quality, we used a sample of 634 massive galaxy clusters from the Planck survey with known spectroscopic redshifts in the range $0.1 < z_{spec} < 0.6$. The accuracy of the photometric redshift estimates for this sample is $\delta z_{phot}/(1+z_{phot}) \approx 0.5$%, the fraction of large deviations is 1.3%. We show that these large deviations arise mainly from the projections of galaxy clusters or other large-scale structures at different redshifts in the X-ray source field. Measuring the infrared (IR) luminosities of galaxy clusters allows one to estimate the reliability of the optical identification of the clusters detected in the SRG/eROSITA survey and to obtain an additional independent measurement of their total gravitational masses, $M_{500}$. We show that the masses $M_{500}$ of the galaxy clusters estimated from their IR luminosity measurements have an accuracy $\sigma_{\lg\,M_{500}} = 0.124$, comparable to the accuracy of the mass estimation for the galaxy clusters from their X-ray luminosities.

We present the results of the optical identification and spectroscopic redshift measurements of 216 galaxy clusters detected in the SRG/eROSITA all-sky X-ray survey. The spectroscopic observations were performed in 2020-2023 with the 6-m BTA telescope at the Special Astrophysical Observatory of the Russian Academy of Sciences, the 2.5-m telescope at the Caucasus Mountain Observatory of the Sternberg Astronomical Institute of the Moscow State University, the 1.6-m AZT-33IK telescope at the Sayan Solar Observatory of the Institute of Solar-Terrestrial Physics of the Siberian Branch of the Russian Academy of Sciences, and the 1.5-m Russian-Turkish telescope (RTT-150) at the T\"{U}B\.{I}TAK Observatory. For all of the galaxy clusters presented here the spectroscopic redshift measurements have been obtained for the first time. Of these, 139 galaxy clusters have been detected for the first time in the SRG/eROSITA survey and 22 galaxy clusters are at redshifts $z_{spec} \gtrsim 0.7$, including three at $z_{spec} \gtrsim 1$. Deep direct images with the rizJK filters have also been obtained for four distant galaxy clusters at $z_{spec} > 0.7$. For these observations the most massive clusters are selected. Therefore, most of the galaxy clusters presented here most likely will be included in the cosmological samples of galaxy clusters from the SRG/eROSITA survey.

The splashback radius was proposed as a physically motivated boundary of clusters as it sets the limit between the infalling and the orbitally dominated regions. However, galaxy clusters are complex objects connected to filaments of the cosmic web from which they accrete matter that disturbs them and modifies their morphology. In this context, estimating the splashback radius and the cluster boundary becomes challenging. In this work, we use a constrained hydrodynamical simulation of the Virgo cluster's replica embedded in its large-scale structure to investigate the impact of its local environment on the splashback radius estimate. We identify the splashback radius from 3D radial profiles of dark matter density, baryons density, and pressure in three regions representative of different dynamical states: accretion from spherical collapse, filaments, and matter outflow. We also identify the splashback radius from 2D-projected radial profiles of observation-like quantities: mass surface density, emission measure, and Compton-y. We show that the splashback radius mainly depends on the dynamics in each region and the physical processes traced by the different probes. We find multiple values for the splashback radius ranging from 3.3$\pm$0.2 to 5.5$\pm$0.3 Mpc. Particularly, in the regions of collapsing and outflowing material, the splashback radii estimated from baryon density and pressure radial profiles overestimate that of the dark matter density profiles, which is considered the reference value originally defined from dark matter simulations. Consequently, caution is required when using the splashback radius as a boundary of clusters, particularly in the case of highly disturbed clusters like Virgo. We also discuss the detection of the splashback radius from pressure radial profiles, which could be more related to an accretion shock, and its detection from stacked radial profiles.

Black holes long-lived enough to be the dark matter have temperatures below the MeV. Since Hawking evaporation is a quasi-thermal process, no GeV emission is predicted to be produced by black holes if they are part, or all, of the cosmological dark matter. However, black holes could be ``spawned'' at late times with masses that correspond to short lifetimes, and as such be significantly hotter and produce particles well in excess of the GeV. Here, we show that such late-forming black holes could, at once, explain the tantalizing excesses found in the gamma radiation from the Galactic center, in the flux of cosmic-ray antiproton, and in the few tentative antihelium events reported by the anti-matter spectrometer AMS-02. We compute accurate predictions for the anti-deuteron, high-energy neutrino, and positron fluxes if this scenario is realized in nature. We find that while the neutrino and positron fluxes are too small compared to the expected background, a significant number of anti-deuteron events is expected both at AMS-02 and at the future General AntiParticle Spectrometer (GAPS).

Supernovae are responsible for the elemental enrichment of the galaxy and some are postulated to leave behind a black hole. In a stellar binary system the supernova pollutes its companion, and the black hole can accrete back its own debris and emit X-rays. In this sequence of events, which is only poorly understood, winds are ejected, and observed through X-ray absorption lines. Measuring abundances of elements in the wind can lead to inferences about the historical explosion and possibly identify the long-gone progenitor of the compact object. Here, we re-analyze the uniquely rich X-ray spectrum of the 2005 outburst of GRO J1655-40. We reconstruct the absorption measure distribution (AMD) of the wind, and find that it increases sharply with ionization from H-like O up to H-like Ca, and then flattens out. The AMD is then used to measure relative abundances of 18 different elements. The present abundances are in partial agreement with a previous work with discrepancies mostly for low-Z elements. The overabundance of odd-Z elements hints at a high-metallicity, high-mass ($\simeq25\,M_\odot$) progenitor. Interestingly, the abundances are different from those measured in the companion atmosphere, indicating that the wind entrains lingering ambient supernova debris. This can be expected since the current total stellar mass of the binary ($<10\,M_\odot$) is much less than the progenitor mass.

Capturing high resolution imagery of the Earth's surface often calls for a telescope of considerable size, even from Low Earth Orbits (LEO). A large aperture often requires large and expensive platforms. For instance, achieving a resolution of 1m at visible wavelengths from LEO typically requires an aperture diameter of at least 30cm. Additionally, ensuring high revisit times often prompts the use of multiple satellites. In light of these challenges, a small, segmented, deployable CubeSat telescope was recently proposed creating the additional need of phasing the telescope's mirrors. Phasing methods on compact platforms are constrained by the limited volume and power available, excluding solutions that rely on dedicated hardware or demand substantial computational resources. Neural Network (NN) are known for their computationally efficient inference and reduced on board requirements. Therefore we developed a NN based method to measure co phasing errors inherent to a deployable telescope. The proposed technique demonstrates its ability to detect phasing error at the targeted performance level (typically a wavefront error (WFE) below 15 nm RMS for a visible imager operating at the diffraction limit) using a point source. The robustness of the NN method is verified in presence of high order aberrations or noise and the results are compared against existing state of the art techniques. The developed NN model ensures its feasibility and provides a realistic pathway towards achieving diffraction limited images.

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will catalogue the light-curves of up to 100 million quasars. Among these there can be up to approximately 100 ultra-compact massive black hole (MBH) binaries, which 5-15 years later can be detected in gravitational waves (GWs) by the Laser Interferometer Space Antenna (LISA). Here we assume that GWs from a MBH binary have been detected by LISA, and we assess whether or not its electromagnetic (EM) counterpart can be uniquely identified in archival LSST data as a periodic quasar. We use the binary's properties derived from the LISA waveform, such as the past evolution of its orbital frequency, its total mass, distance and sky localization, to predict the redshift, magnitude and historical periodicity of the quasar expected in the archival LSST data. We then use Monte Carlo simulations to compute the false alarm probability, i.e. the number of quasars in the LSST catalogue matching these properties by chance, based on the (extrapolated) quasar luminosity function, the sampling cadence of LSST, and intrinsic ``damped random walk (DRW)" quasar variability. We perform our analysis on four fiducial LISA binaries, with total masses and redshifts of $(M_{\rm bin}/{\rm M_{\odot}},z) = (3\times10^5,0.3)$, $(3\times10^6,0.3)$, $(10^7,0.3)$ and $(10^7,1)$. While DRW noise and aliasing due to LSST's cadence can produce false periodicities by chance, we find that the frequency chirp of the LISA source during the LSST observations washes out these noise peaks and allows the genuine source to stand out in Lomb-Scargle periodograms. We find that all four fiducial binaries yield excellent chances to be uniquely identified, with false alarm probabilities below $10^{-5}$, a week or more before their merger. This then enables deep follow-up EM observations targeting the individual EM counterparts during their inspiral stage.

Context. The complexity and variety exhibited by the light curves of long gamma-ray bursts (GRBs) enclose a wealth of information that still awaits being fully deciphered. Despite the tremendous advance in the knowledge of the energetics, structure, and composition of the relativistic jet that results from the core collapse of the progenitor star, the nature of the inner engine, how it powers the relativistic outflow, and the dissipation mechanisms remain open issues. Aims. A promising way to gain insights is describing GRB light curves as the result of a common stochastic process. In the Burst And Transient Source Experiment (BATSE) era, a stochastic pulse avalanche model was proposed and tested through the comparison of ensemble-average properties of simulated and real light curves. Here we aim to revive and further test this model. Methods. We apply it to two independent data sets, BATSE and Swift/BAT, through a machine learning approach: the model parameters are optimised using a genetic algorithm. Results. The average properties are successfully reproduced. Notwithstanding the different populations and passbands of both data sets, the corresponding optimal parameters are interestingly similar. In particular, for both sets the dynamics appears to be close to a critical state, which is key to reproduce the observed variety of time profiles. Conclusions. Our results propel the avalanche character in a critical regime as a key trait of the energy release in GRB engines, which underpins some kind of instability.

We present a new perspective on the symmetries that govern the formation of large-scale structures across the Universe, particularly focusing on the transition from the seeds of galaxy clusters to the seeds of galaxies themselves. We address two main features of cosmological fluid dynamics pertaining to both the linear and non-linear regimes. The linear dynamics of cosmological perturbations within the Hubble horizon is characterized by the Jeans length, which separates stable configurations from unstable fluctuations due to the gravitational effect on sufficiently large (and therefore, massive enough) overdensities. On the other hand, the non-linear dynamics of the cosmological fluid is associated with a turbulent behavior once the Reynolds numbers reach a sufficiently high level. This turbulent regime leads to energy dissipation across smaller and smaller scales, resulting in a fractal distribution of eddies throughout physical space. The proposed scenario suggests that the spatial scale of eddy formation is associated with the Jeans length of various levels of fragmentation from an original large-scale structure. By focusing on the fragmentation of galaxy cluster seeds versus galaxy seeds, we arrived at a phenomenological law that links the ratio of the two structure densities to the number of galaxies in each cluster and to the Hausdorff number of the Universe matter distribution. Finally, we introduced a primordial magnetic field and studied its influence on the Jeans length dynamics. The resulting anisotropic behavior of the density contrast led us to infer that the main features of the turbulence could be reduced to a 2D Euler equation. Numerical simulations showed that the two lowest wavenumbers contained the major energy contribution of the spectrum.

Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by fNL, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spectrum. We propose a framework to measure fNL using the configuration space two-point correlation function (2pcf) monopole and three-point correlation function (3pcf) monopole of survey tracers. Our model estimates the effect of the scale-dependent bias induced by the presence of PNG on the 2pcf and 3pcf from the clustering of simulated dark matter halos. We describe how this effect may be scaled to an arbitrary tracer of the cosmological matter density. The 2pcf and 3pcf of this tracer are measured to constrain the value of fNL. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our model, and forecast the ability to constrainfNL to a precision of sigma(fNL) = 22 with one year of DESI survey data.

The Solar Orbiter (SO) mission provides the opportunity to study the evolution of solar wind turbulence. We use SO observations of nine extended intervals of homogeneous turbulence to determine when turbulent magnetic field fluctuations may be characterized as: (i) wave-packets and (ii) coherent structures (CS). We perform the first systematic scale-by-scale decomposition of the magnetic field using two wavelets known to resolve wave-packets and discontinuities, the Daubechies 10 (Db10) and Haar respectively. The probability distributions (pdfs) of turbulent fluctuations on small scales exhibit stretched tails, becoming Gaussian at the outer scale of the cascade. Using quantile-quantile plots, we directly compare the wavelet fluctuations pdfs, revealing three distinct regimes of behaviour. Deep within the inertial range (IR) both decompositions give essentially the same fluctuation pdfs. Deep within the kinetic range (KR) the pdfs are distinct as the Haar wavelet fluctuations have larger variance and more extended tails. On intermediate scales, spanning the IR-KR break, the pdf is composed of two populations: a core of common functional form containing $\sim97\%$ of fluctuations, and tails which are more extended for Haar fluctuations than Db10 fluctuations. This establishes a crossover between wave-packet (core) and CS (tail) phenomenology in the IR and KR respectively. The range of scales where the pdfs are $2$-component is narrow at $0.9$ au ($4-16$ s) and broader ($0.5-8$ s) at $0.4$ au. As CS and wave-wave interactions are both candidates to mediate the turbulent cascade, these results offer new insights into the distinct physics of the IR and KR.

Discovering exoplanets in orbit around distant stars via direct imaging is fundamentally impeded by the high dynamic range between the star and the planet. Coronagraphs strive to increase the signal-to-noise ratio of exoplanet signatures by optically rejecting light from the host star while leaving light from the exoplanet mostly unaltered. However it is unclear whether coronagraphs constitute an optimal strategy for attaining fundamental limits relevant exoplanet discovery. In this work, we report the quantum information limits of exoplanet detection and localization specified by the Quantum Chernoff Exponent (QCE) and the Quantum Fisher Information Matrix (QFIM) respectively. In view of these quantum limits, we assess and compare several high-performance coronagraph designs that theoretically achieve total rejection of an on-axis point source. We find that systems which exclusively eliminate the fundamental mode of the telescope without attenuating higher-order orthogonal modes are quantum-optimal in the regime of high star-planet contrasts. Importantly, the QFIM is shown to persist well below the diffraction-limit of the telescope, suggesting that quantum-optimal coronagraphs may further expand the domain of accessible exoplanets.

We consider here natural inflation in the low energy (two-derivative) metric-affine theory containing only the minimal degrees of freedom in the inflationary sector, i.e. the massless graviton and the pseudo-Nambu-Goldstone boson (PNGB). This theory contains the Ricci-like and parity-odd Holst invariants together with non-minimal couplings between the PNGB and the above-mentioned invariants. The Palatini and Einstein-Cartan realizations of natural inflation are particular cases of our construction. Explicit models of this type are shown to admit an explicit UV completion in a QCD-like theory with a Planckian confining scale. Moreover, for these models, we find regions of the parameter space where the inflationary predictions agree with the most recent observations at the $2\sigma$ level. We find that in order to enter the $1\sigma$ region it is necessary (and sufficient) to have a finite value of the Barbero-Immirzi parameter and a sizable non-minimal coupling between the inflaton and the Holst invariant (with sign opposite to the Barbero-Immirzi parameter). Indeed, in this case the potential of the canonically normalized inflaton develops a plateau as shown analytically.

The percent-level precision attained by modern cosmic ray (CR) observations motivates reaching a comparable or better control of theoretical uncertainties. Here we focus on energy-loss processes affecting low-energy CR protons ($\sim 0.1-5$ GeV), where the experimental errors are small and collisional effects play a comparatively larger role with respect to collisionless transport ones. We study three aspects of the problem: i) We quantitatively assess the role of the nuclear elastic cross-section, for the first time, providing analytical formulae for the stopping power and inelasticity. ii) We discuss the error arising from treating both elastic and pion production inelastic interactions as continuous energy loss processes, as opposed to catastrophic ones. The former is the approximation used in virtually all modern numerical calculations. iii) We consider sub-leading effects such as relativistic corrections, radiative and medium processes in ionization energy-losses. Our analysis reveals that neglecting i) leads to errors close to 1%, notably around and below 1 GeV; neglecting ii) leads to errors reaching about 3% within the considered energy range; iii) contributes to a minor effect, gauged at the level of 0.1%. Consequently, while iii) can currently be neglected, ii) warrants consideration, and we also recommend incorporating i) into computations. We conclude with some perspectives on further steps to be taken towards a high-precision goal of theoretical CR predictions regarding the treatment of energy-losses.

Historically, astronomy has prioritized visuals to present information, with scientists and communicators overlooking the critical need to communicate astrophysics with blind or low-vision audiences and provide novel channels for sighted audiences to process scientific information. This study sonified NASA data of three astronomical objects presented as aural visualizations, then surveyed blind or low-vision and sighted individuals to elicit feedback on the experience of these pieces as it relates to enjoyment, education, and trust of the scientific data. Data analyses from 3,184 sighted or blind or low-vision survey participants yielded significant self-reported learning gains and positive experiential responses. Results showed that astrophysical data engaging multiple senses could establish additional avenues of trust, increase access, and promote awareness of accessibility in sighted and blind or low-vision communities.

A unified analytical solution is presented for constructing the phase space near collinear libration points in the Circular Restricted Three-body Problem (CRTBP), encompassing Lissajous orbits and quasihalo orbits, their invariant manifolds, as well as transit and non-transit orbits. Traditional methods could only derive separate analytical solutions for the invariant manifolds of Lissajous orbits and halo orbits, falling short for the invariant manifolds of quasihalo orbits. By introducing a coupling coefficient {\eta} and a bifurcation equation, a unified series solution for these orbits is systematically developed using a coupling-induced bifurcation mechanism and Lindstedt-Poincar\'e method. Analyzing the third-order bifurcation equation reveals bifurcation conditions for halo orbits, quasihalo orbits, and their invariant manifolds. Furthermore, new families of periodic orbits similar to halo orbits are discovered, and non-periodic/quasi-periodic orbits, such as transit orbits and non-transit orbits, are found to undergo bifurcations. When {\eta} = 0, the series solution describes Lissajous orbits and their invariant manifolds, transit, and non-transit orbits. As {\eta} varies from zero to non-zero values, the solution seamlessly transitions to describe quasihalo orbits and their invariant manifolds, as well as newly bifurcated transit and non-transit orbits. This unified analytical framework provides a more comprehensive understanding of the complex phase space structures near collinear libration points in the CRTBP.

We present new, simple analytical formulas to accurately describe light rays in spherically symmetric static spacetimes. These formulas extend those introduced by Beloborodov and refined by Poutanen for the Schwarzschild metric. Our enhanced formulas are designed to be applicable to a broader range of spacetimes, making them particularly valuable for describing phenomena around compact objects like neutron stars and black holes. As an illustration of their application, we present analytical studies of images of thin accretion disks surrounding black holes and explore their associated polarimetry.

The scaling of the relativistic reconnection outflow speed is studied in the presence of both shear flows parallel to the reconnecting magnetic fields and guide fields pointing out of the reconnection plane. In nonrelativistic reconnection, super-Alfv\'enic shear flows have been found to suppress reconnection. We extend the analytical model of this phenomenon to the relativistic regime and find similar behavior, which is confirmed by particle-in-cell simulations. Unlike the nonrelativistic limit, the addition of a guide field lowers the in-plane Alfv\'en velocity, contributing to slower outflow jets and the more efficient suppression of reconnection in strongly magnetized plasmas.

High frequency gravitational waves (HFGWs) are predicted in various exotic scenarios involving both cosmological and astrophysical sources. These elusive signals have recently sparked the interest of a diverse community of researchers, due to the possibility of HFGW detection in the laboratory through graviton-photon conversion in strong magnetic fields. Notable examples include the redesign of the resonant cavities currently under development to detect the cosmic axion. In this work, we derive the sensitivities of some existing and planned resonant cavities to detect a HFGW background. As a concrete scenario, we consider the collective signals that originate from the merging of compact objects, such as two primordial black holes (PBHs) in the asteroid mass window. Our findings improve over existing work by explicitly discussing and quantifying the loss in the experimental reach due to the actual coherence of the source. We elucidate on the approach we adopt in relation with recent literature on the topic. Most notably, we give a recipe for the estimate of the stochastic background that focuses on the presence of the signal in the cavity at all times and showing that, in the relevant PBH mass region, the signal is dominated by coherent binary mergers.

The promise of multi-messenger astronomy relies on the rapid detection of gravitational waves at very low latencies ($\mathcal{O}$(1\,s)) in order to maximize the amount of time available for follow-up observations. In recent years, neural-networks have demonstrated robust non-linear modeling capabilities and millisecond-scale inference at a comparatively small computational footprint, making them an attractive family of algorithms in this context. However, integration of these algorithms into the gravitational-wave astrophysics research ecosystem has proven non-trivial. Here, we present the first fully machine learning-based pipeline for the detection of gravitational waves from compact binary coalescences (CBCs) running in low-latency. We demonstrate this pipeline to have a fraction of the latency of traditional matched filtering search pipelines while achieving state-of-the-art sensitivity to higher-mass stellar binary black holes.

The aim of this study is to investigate the effect of dark matter (DM) on $f$-mode oscillations in DM admixed neutron stars (NSs). We consider hadronic matter modeled by the relativistic mean field model and the DM model based on the neutron decay anomaly. We study the non-radial $f$-mode oscillations for such DM admixed NS in a full general relativistic framework. We investigate the impact of DM, DM self-interaction, and DM fraction on the $f$-mode characteristics. We derive relations encoding the effect of DM on $f$-mode parameters. We then perform a systematic study by varying all the model parameters within their known uncertainty range and obtain a universal relation for the DM fraction based on the total mass of the star and DM self-interaction strength. We also perform a correlation study among model parameters and NS observables, in particular, $f$-mode parameters. Finally, we check the $f$-mode universal relations (URs) for the case of DM admixed NSs and demonstrate the existence of a degeneracy between purely hadronic NSs and DM admixed NSs.

In this paper, we explore the effects of General Relativity modification on the Mass-Radius relations of Neutron Stars induced by the presence of the Quintessence field. We consider, in particular, the Kiselev model, according to which the Quintessence field, being present in the entire Universe, might also be present around massive objects. Considering the Equation of State (EoS) for Baryonic matter BSk22 derived by A. Y. Potekhin et al., we infer the upper limit for NS masses in the presence of Quintessence. The presence of Quintessence generates a peculiar effect for which the Mass-Radius relation is unvaried and therefore the presence of Quintessence is indistinguishable from ordinary matter, at least for the Kiselev model studied in this paper.

Neutron stars (NSs) are compact objects that host the densest forms of matter in the observable universe, providing unique opportunities to study the behaviour of matter at extreme densities. While precision measurements of NS masses through pulsar timing have imposed effective constraints on the equation of state (EoS) of dense matter, accurately determining the radius or moment of inertia (MoI) of a NS remains a major challenge. This article presents a detailed review on measuring the Lense-Thirring (LT) precession effect in the orbit of binary pulsars, which would give access to the MoI of NSs and offer further constraints on the EoS. We discuss the suitability of certain classes of binary pulsars for measuring the LT precession from the perspective of binary star evolution, and highlight five pulsars that exhibit properties promising to realise these goals in the near future. Finally, discoveries of compact binaries with shorter orbital periods hold the potential to greatly enhance measurements of the MoI of NSs. The MoI measurements of binary pulsars are pivotal to advancing our understanding of matter at supranuclear densities as well as improving the precision of gravity tests, such as the orbital decay due to gravitational wave emission and of tests of alternative gravity theories.