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

In the current panorama of large surveys, the vast amount of data obtained with different methods, data types, formats, and stellar samples, is making an efficient use of the available information difficult. The Survey of Surveys is a project to critically compile survey results in a single catalogue, facilitating the scientific use of the available information. In this second release, we present two new catalogs of stellar parameters (Teff, logg, and [Fe/H]). To build the first catalog, SoS-Spectro, we calibrated internally and externally stellar parameters from five spectroscopic surveys (APOGEE, GALAH, Gaia-ESO, RAVE, and LAMOST) and externally on the PASTEL database. The second catalog, SoS-ML catalog, is obtained by using SoS-Spectro as a reference to train a multi-layer perceptron, which predicts stellar parameters based on two photometric surveys, SDSS and SkyMapper. As a novel approach, we build on previous parameters sets, from Gaia DR3 and Andrae et al. (2023), aiming to improve their precision and accuracy. We obtain a catalog of stellar parameters for around 23 millions of stars, which we make publicly available. We validate our results with several comparisons with other machine learning catalogs, stellar clusters, and astroseismic samples. We find substantial improvements in the parameters estimates compared to other Machine Learning methods in terms of precision and accuracy, especially in the metal-poor range, as shown in particular when validating our results with globular clusters. We believe that there are two reasons behind our improved results at the low-metallicity end: first, our use of a reference catalog, the SoS-Spectro, which is calibrated using high-resolution spectroscopic data; and second, our choice to build on pre-existing parameter estimates from em Gaia and Andrae et al., rather than attempting to obtain our predictions from survey data alone.

Euclid is expected to establish new state-of-the-art constraints on extensions beyond the standard LCDM cosmological model by measuring the positions and shapes of billions of galaxies. Specifically, its goal is to shed light on the nature of dark matter and dark energy. Achieving this requires developing and validating advanced statistical tools and theoretical prediction software capable of testing extensions of the LCDM model. In this work, we describe how the Euclid likelihood pipeline, Cosmology Likelihood for Observables in Euclid (CLOE), has been extended to accommodate alternative cosmological models and to refine the theoretical modelling of Euclid primary probes. In particular, we detail modifications made to CLOE to incorporate the magnification bias term into the spectroscopic two-point correlation function of galaxy clustering. Additionally, we explain the adaptations made to CLOE's implementation of Euclid primary photometric probes to account for massive neutrinos and modified gravity extensions. Finally, we present the validation of these CLOE modifications through dedicated forecasts on synthetic Euclid-like data by sampling the full posterior distribution and comparing with the results of previous literature. In conclusion, we have identified in this work several functionalities with regards to beyond-LCDM modelling that could be further improved within CLOE, and outline potential research directions to enhance pipeline efficiency and flexibility through novel inference and machine learning techniques.

We present a high-resolution X-ray spectroscopic study of the Narrow-Line Seyfert 1 galaxy NGC 4051 using two XMM-Newton high-resolution Reflection Grating Spectrometer (RGS) observations. The spectra reveal three distinct layers of photoionized gas flowing outward from the central black hole: a low-ionization phase (LIP), a higher-ionization phase (HIP), and a high-velocity and high ionization phase (HVIP). Each absorber leaves characteristic imprints on the soft X-ray spectrum. While the LIP and HVIP are fully consistent with being in ionization equilibrium with the central radiation field over the course of the $\sim$250 ks spanned by the two observations, the HIP shows a significant change in ionization ($3.8\sigma$), suggesting non-equilibrium. By modeling the two spectra with our time-dependent photoionization code (TEPID), we constrain the density of the HIP gas to $\log n_{\rm H}=7.7^{+0.2}_{-0.9}$ and estimate its distance to be about $R=0.45^{+0.80}_{-0.09}$ light-days from the black hole, corresponding to $R=4000^{+7000}_{-800}$ gravitational radii. In contrast, the narrow soft X-ray emission lines remain constant, consistent with an origin in the more extended narrow-line region. Our results show the value of combining high-resolution and time-resolved spectroscopy to probe the structure, physical conditions, and variability of AGN outflows.

We present deep Chandra X-ray observations of NGC 5005, a LINER-dominated galaxy previously reported to host a broad H$\alpha$ emission line. The diffuse soft X-ray emission (<3 keV) extends out to $\sim$800 pc, while harder emission (>3 keV) is confined to the central $\sim$400 pc. Spatially resolved spectroscopy of the nuclear (r<150 pc) and extended ($150

We report the discovery of X-ray polarization from the X-ray-bright filament. G0.13-0.11 in the Galactic center (GC) region. This filament features a bright, hard X-ray source that is most plausibly a pulsar wind nebula (PWN) and an extended and structured diffuse component. Combining the polarization signal from IXPE with the imaging/spectroscopic data from Chandra, we find that X-ray emission of G0.13-0.11 is highly polarized PD=$57(\pm18)$% in the 3-6 keV band, while the polarization angle is PA=$21^\circ(\pm9^\circ)$. This high degree of polarization proves the synchrotron origin of the X-ray emission from G0.13-0.11. In turn, the measured polarization angle implies that the X-ray emission is polarized approximately perpendicular to a sequence of nonthermal radio filaments that may be part of the GC Radio Arc. The magnetic field on the order of $100\,{\rm\mu G}$ appears to be preferentially ordered along the filaments. The above field strength is the fiducial value that makes our model self-consistent, while the other conclusions are largely model independent.

Recent James Webb Space Telescope (JWST) observations have revealed a population of sources with a compact morphology and a `v-shaped' continuum, namely blue at rest-frame \lambda&lt;4000A and red at longer wavelengths. The nature of these sources, called `little red dots' (LRDs), is still debated, since it is unclear if they host active galactic nuclei (AGN) and their number seems to drastically drop at z<4. We utilise the 63 $deg^2$ covered by the quick Euclid Quick Data Release (Q1) to extend the search for LRDs to brighter magnitudes and to lower z than what has been possible with JWST to have a broader view of the evolution of this peculiar galaxy population. The selection is done by fitting the available photometric data (Euclid, Spitzer/IRAC, and ground-based griz data) with two power laws, to retrieve the rest-frame optical and UV slopes consistently over a large redshift range (i.e, z<7.6). We exclude extended objects and possible line emitters, and perform a visual inspection to remove imaging artefacts. The final selection includes 3341 LRD candidates from z=0.33 to z=3.6, with 29 detected in IRAC. Their rest-frame UV luminosity function, in contrast with previous JWST studies, shows that the number density of LRD candidates increases from high-z down to z=1.5-2.5 and decreases at even lower z. Less evolution is apparent focusing on the subsample of more robust LRD candidates having IRAC detections, which is affected by low statistics and limited by the IRAC resolution. The comparison with previous quasar UV luminosity functions shows that LRDs are not the dominant AGN population at z<4. Follow-up studies of these LRD candidates are key to confirm their nature, probe their physical properties and check for their compatibility with JWST sources, since the different spatial resolution and wavelength coverage of Euclid and JWST could select different samples of compact sources.

The origin of the high-energy emission in astrophysical jets from black holes is a highly debated issue. This is particularly true for jets from supermassive black holes that are among the most powerful particle accelerators in the Universe. So far, the addition of new observations and new messengers have only managed to create more questions than answers. However, the newly available X-ray polarization observations promise to finally distinguish between emission models. We use extensive multiwavelength and polarization campaigns as well as state-of-the-art polarized spectral energy distribution models to attack this problem by focusing on two X-ray polarization observations of blazar BL Lacertae in flaring and quiescent $\gamma$-ray states. We find that regardless of the jet composition and underlying emission model, inverse-Compton scattering from relativistic electrons dominates at X-ray energies.

Long-period variable stars (LPVs) are pulsating red giants, primarily in the asymptotic giant branch phase, and they include both Miras and semi-regular variables (SRVs). Their period-age and period-luminosity relations enable us to trace different stellar populations, as they are intrinsically very bright and cover a wide range in distances and ages. The purpose of this study is to establish a census of LPV stars in a region close to the Galactic center, using the six-year database of the Vista Variables in the V\'ia L\'actea (VVV) ESO Public Survey, as well as to describe the methodology that was employed to search for and characterize LPVs using VVV data. Near-IR surveys such as VVV provide a unique opportunity to probe the high-extinction innermost regions of the Milky Way. The detection and analysis of the intrinsically bright Miras in this region could provide us with an excellent probe of the properties of the Milky Way far behind its bulge. We used point-spread function photometry for all available $K_{s}$-band images in ten VVV tiles, covering $16.4~\deg^2$ in total, overlapping fields observed in the course of the Optical Gravitational Lensing Experiment (OGLE)-III survey. We designed a method to select LPV candidates, and we used the known variables from OGLE-III and other known variables from the literature to test our approach. The reduced $\chi^2$ statistic, along with the flux-independent index $K_{(fi)}$, were used in our analysis. The Lomb-Scargle period search method, Fourier analysis, template fitting, and visual inspection were then performed to refine our sample and characterize the properties of the stars included in our catalog. A final sample of 130 Mira candidates, of which 129 are new discoveries, was thus obtained, with periods in the range between about 80 and 1400~days.

XMM-Newton and NuSTAR multiple exposures allow us to disentangle the different emission components of active galactic nuclei (AGNs) and to study the evolution of their different spectral features. In this work, we present the timing and spectral properties of five simultaneous XMM-Newton and NuSTAR observations of the Narrow Line Seyfert 1 galaxy Mrk 359. We aim to provide the first broadband spectral modeling of Mrk 359 describing its emission spectrum from the UV up to the hard X-rays. To do this, we performed temporal and spectral data analysis, characterising the amplitude and spectral changes of the Mrk 359 time series and computing the 2-10 keV normalised excess variance. The spectral broadband modelling assumes the standard hot Comptonising corona and reflection component, while for the soft excess we tested two different models: a warm, optically thick Comptonising corona (the two-corona model) and a reflection model in which the soft-excess is the result of a blurred reflected continuum and line emission (the reflection model). High and low flux states were observed during the campaign. The former state has a softer spectral shape, while the latter shows a harder one. The photon index is in the 1.75-1.89 range, and only a lower limit to the hot-corona electron temperature can be found. A constant reflection component, likely associated with distant matter, is observed. Regarding the soft excess, we found that among the reflection models we tested, the one providing the better fit (reduced $\chi^2$=1.14) is the high-density one. However, a significantly better fit (reduced $\chi^2$=1.08) is found by modelling the soft excess with a warm Comptonisation model. The present analysis suggests the two-corona model as the best scenario for the optical-UV to X-ray emission spectrum of Mrk 359.

The composition of the solar corona differs from that of the photosphere, with the plasma thought to fractionate in the solar chromosphere according to the First Ionisation Potential (FIP) of the different elements. This produces a FIP bias, wherein elements with a low FIP are preferentially enhanced in the corona compared to their photospheric abundance, but direct observations of this process remain elusive. Here we use a series of spectroscopic observations of Active Region AR 12759 as it transited the solar disc over a period of 6 days from 2-7 April 2020 taken using the Hinode Extreme ultraviolet Imaging Spectrometer (EIS) and Interface Region Imaging Spectrograph (IRIS) instruments to look for signatures of plasma fractionation in the solar chromosphere. Using the Si X/S X and Ca XIV/Ar XIV diagnostics, we find distinct differences between the FIP bias of the leading and following polarities of the active region. The widths of the IRIS Si IV lines exhibited clear differences between the leading and following polarity regions, indicating increased unresolved wave activity in the following polarity region compared to the leading polarity region, with the chromospheric velocities derived using the Mg II lines exhibiting comparable, albeit much weaker, behaviour. These results are consistent with plasma fractionation via resonant/non-resonant waves at different locations in the solar chromosphere following the ponderomotive force model, and indicate that IRIS could be used to further study this fundamental physical process.

We present SABAT ("Sapienza" Balloon Trajectory Simulation Code), a Fortran software for the simulation of zero-pressure scientific balloon trajectories at stratospheric altitudes. Particular attention is made to the correct reproduction of the balloon ascending phase to the floating altitude. In order to obtain trajectories in agreement with those monitored during actual balloon launches and flights, SABAT features a balloon dynamical and geometric model, and a thermal model for the heat exchanges between the vehicle and the surrounding environment. In addition, both a synthetic description of the Earth's atmosphere and real data from radiosondes can be used. The validation tests give results in agreement with the characteristics of real balloon motion. Future developments of the code include optimization for balloons flying at polar latitudes, an improved treatment of the drag force acting on the balloon structure and the development of a user-friendly graphical interface.

An increasing number of pulsar wind nebulae (PWNe) are being identified in the TeV band by ground-based Imaging Air Cherenkov Telescopes such that they constitute the dominant source class of Galactic TeV emitters. However, MeV-GeV PWN counterparts are still largely lacking. To date, only a dozen PWNe are identified by the Fermi-Large Area Telescope (LAT) in the MeV-GeV band. Most PWNe are located along the Galactic plane embedded within the prominent, diffuse Galactic gamma-ray emission, which makes these sources difficult to disentangle from the bright diffuse background. We present a systematic search for gamma-ray counterparts to known PWNe in the 300MeV-2TeV energy band using the Fermi-LAT. We target locations of previously identified PWNe that lack detected Fermi-LAT pulsars to minimize associated pulsar contamination. The sample includes 6 previously identified Fermi-LAT PWNe and 8 Fermi-LAT sources associated with PWNe. We report the analysis of 58 regions of interest and classify detected sources as either a likely PWN or a candidate PWN counterpart based on their morphological and spectral characteristics across the broadband spectrum. There are 9 unidentified Fermi-LAT sources that we consider as likely PWN counterparts, which, if confirmed to be PWNe, would greatly increase the PWN population detected by the Fermi-LAT from 12 to 21. The remaining Fermi-LAT detected sources are considered weaker PWN candidates. A second approach in the systematic search for gamma-ray emitting PWNe will involve studying the off-pulse phases of Fermi-LAT pulsars for the presence of an obscured PWN and will be reported in a subsequent paper.

We present a machine learning method to assign stellar parameters (temperature, surface gravity, metallicity) to the photometric data of large photometric surveys such as SDSS and SKYMAPPER. The method makes use of our previous effort in homogenizing and recalibrating spectroscopic data from surveys like APOGEE, GALAH, or LAMOST into a single catalog, which is used to inform a neural network. We obtain spectroscopic-quality parameters for millions of stars that have only been observed photometrically. The typical uncertainties are of the order of 100K in temperature, 0.1 dex in surface gravity, and 0.1 dex in metallicity and the method performs well down to low metallicity, were obtaining reliable results is known to be difficult.

Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid telescope, using both simulated pre-launch images and the first early-release observations of the Fornax galaxy cluster. The Euclid Wide Survey will provide high-spatial resolution VIS imaging in the broad IE band as well as near-infrared photometry (YE, JE, and HE). We estimate that the galaxies within 100 Mpc in the footprint of the Euclid survey host around 830 000 EGCs of which about 350 000 are within the survey's detection limits. For about half of these EGCs, three infrared colours will be available as well. For any galaxy within 50Mpc the brighter half of its GC luminosity function will be detectable by the Euclid Wide Survey. The detectability of EGCs is mainly driven by the residual surface brightness of their host galaxy. We find that an automated machine-learning EGC-classification method based on real Euclid data of the Fornax galaxy cluster provides an efficient method to generate high purity and high completeness GC candidate catalogues. We confirm that EGCs are spatially resolved compared to pure point sources in VIS images of Fornax. Our analysis of both simulated and first on-sky data show that Euclid will increase the number of GCs accessible with high-resolution imaging substantially compared to previous surveys, and will permit the study of GCs in the outskirts of their hosts. Euclid is unique in enabling systematic studies of EGCs in a spatially unbiased and homogeneous manner and is primed to improve our understanding of many understudied aspects of GC astrophysics.

Galaxy morphologies and shape orientations are expected to correlate with their large-scale environment, since they grow by accreting matter from the cosmic web and are subject to interactions with other galaxies. Cosmic filaments are extracted in projection from the Euclid Quick Data Release 1 (covering 63.1 $\mathrm{deg}^2$) at $0.5 10^{10} M_\odot$) in the projected cosmic web is analysed as a function of morphology measured from VIS data. Specifically, the 2D alignment of galaxy shapes with large-scale filaments is quantified as a function of S\'ersic indices and masses. We find the known trend that more massive galaxies are closer to filament spines. At fixed stellar masses, morphologies correlate both with densities and distances to large-scale filaments. In addition, the large volume of this data set allows us to detect a signal indicating that there is a preferential alignment of the major axis of massive early-type galaxies along projected cosmic filaments. Overall, these results demonstrate our capabilities to carry out detailed studies of galaxy environments with Euclid, which will be extended to higher redshift and lower stellar masses with the future Euclid Deep Survey.

RR Lyrae (RRL) are old, low-mass radially pulsating variable stars in their core helium burning phase. They are popular stellar tracers and primary distance indicators, since they obey to well defined period-luminosity relations in the near-infrared regime. Their photometric identification is not trivial, indeed, RRL samples can be contaminated by eclipsing binaries, especially in large datasets produced by fully automatic pipelines. Interpretable machine-learning approaches for separating eclipsing binaries from RRL are thus needed. Ideally, they should be able to achieve high precision in identifying RRL while generalizing to new data from different instruments. In this paper, we train a simple logistic regression classifier on Catalina Sky Survey (CSS) light curves. It achieves a precision of 87% at 78% recall for the RRL class on unseen CSS light curves. It generalizes on out-of-sample data (ASAS/ASAS-SN light curves) with a precision of 85% at 96% recall. We also considered a L1-regularized version of our classifier, which reaches 90% sparsity in the light-curve features with a limited trade-off in accuracy on our CSS validation set and -- remarkably -- also on the ASAS/ASAS-SN light curve test set. Logistic regression is natively interpretable, and regularization allows us to point out the parts of the light curves that matter the most in classification. We thus achieved both good generalization and full interpretability.

We report on the coordinated observations of the neutron star low-mass X-ray binary (NS-LMXB) \gx in X-rays (IXPE, NICER, Nustar and INTEGRAL), optical (REM and LCO), near-infrared (REM), mid-infrared (VLT VISIR), and radio (ATCA). This Z-source was observed by \IXPE twice in March-April 2023 (Obs. 1 and 2). In the radio band, the source was detected, but only upper-limits to the linear polarization were obtained at a $3\sigma$ level of $6.1\%$ at 5.5 GHz and $5.9\%$ at 9 GHz in Obs.~1 and $12.5\%$ at 5.5~GHz and $20\%$ at 9~GHz in Obs.~2. The mid-IR, near-IR and optical observations suggest the presence of a compact jet which peaks in the mid- or far-IR. The X-ray polarization degree was found to be $3.7\% \pm 0.4 \%$ (at $90\%$ confidence level) during Obs.~1 when the source was in the horizontal branch of the Z-track and $1.8\% \pm 0.4 \%$ during Obs.~2 when the source was in the normal-flaring branch. These results confirm the variation of polarization degree as a function of the position of the source in the color-color diagram as for previously observed Z-track sources (Cyg~X-2 and XTE~1701$-$462). Evidence for a variation of the polarization angle $\sim 20^\circ$ with energy is found in both observations, likely related to the different, non-orthogonal polarization angles of the disk and Comptonization components which peak at different energies.

We present X-ray polarimetry observations from the Imaging X-ray Polarimetry Explorer (IXPE) of three low spectral peak and one intermediate spectral peak blazars, namely 3C 273, 3C 279, 3C 454.3, and S5 0716+714. For none of these objects was IXPE able to detect X-ray polarization at the 3$\sigma$ level. However, we placed upper limits on the polarization degree at $\sim$10-30\%. The undetected polarizations favor models where the X-ray band is dominated by unpolarized photons upscattered by relativistic electrons in the jets of blazars, although hadronic models are not completely eliminated. We discuss the X-ray polarization upper limits in the context of our contemporaneous multiwavelength polarization campaigns.

The Imaging X-ray Polarimetry Explorer (IXPE) measured with high significance the X-ray polarization of the brightest Z-source Scorpius X-1, resulting in the nominal 2-8 keV energy band in a polarization degree of 1.0(0.2)% and a polarization angle of 8(6)° at 90% of confidence level. This observation was strictly simultaneous with observations performed by NICER, NuSTAR, and Insight-HXMT, which allowed for a precise characterization of its broad-band spectrum from soft to hard X-rays. The source has been observed mainly in its soft state, with short periods of flaring. We also observed low-frequency quasi-periodic oscillations. From a spectro-polarimetric analysis, we associate a polarization to the accretion disk at <3.2% at 90% of confidence level, compatible with expectations for an electron-scattering dominated optically thick atmosphere at the Sco X-1 inclination of 44°; for the higher-energy Comptonized component, we obtain a polarization of 1.3(0.4)%, in agreement with expectations for a slab of Thomson optical depth of ~7 and an electron temperature of ~3 keV. A polarization rotation with respect to previous observations by OSO-8 and PolarLight, and also with respect to the radio-jet position angle, is observed. This result may indicate a variation of the polarization with the source state that can be related to relativistic precession or to a change in the corona geometry with the accretion flow.