Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $\alpha_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $\alpha_{\rm CRp}$ down to about $\Delta\alpha_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with \tau_{\chi}&gt;10^{27}s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario.

The Euclid mission of the European Space Agency will deliver weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and extensions thereof. We present forecasts from the combination of these surveys on the sensitivity to cosmological parameters including the summed neutrino mass $M_\nu$ and the effective number of relativistic species $N_{\rm eff}$ in the standard $\Lambda$CDM scenario and in a scenario with dynamical dark energy ($w_0 w_a$CDM). We compare the accuracy of different algorithms predicting the nonlinear matter power spectrum for such models. We then validate several pipelines for Fisher matrix and MCMC forecasts, using different theory codes, algorithms for numerical derivatives, and assumptions concerning the non-linear cut-off scale. The Euclid primary probes alone will reach a sensitivity of $\sigma(M_\nu)=$56meV in the $\Lambda$CDM+$M_\nu$ model, whereas the combination with CMB data from Planck is expected to achieve $\sigma(M_\nu)=$23meV and raise the evidence for a non-zero neutrino mass to at least the $2.6\sigma$ level. This can be pushed to a $4\sigma$ detection if future CMB data from LiteBIRD and CMB Stage-IV are included. In combination with Planck, Euclid will also deliver tight constraints on $\Delta N_{\rm eff}< 0.144$(95$\Lambda$CDM+$M_\nu$+$N_{\rm eff}$model, or$\Delta N_{\rm eff}< 0.063$when future CMB data are included. When floating$(w_0, w_a)$, we find that the sensitivity to$N_{\rm eff}$ remains stable, while that to $M_\nu$ degrades at most by a factor 2. This work illustrates the complementarity between the Euclid spectroscopic and imaging/photometric surveys and between Euclid and CMB constraints. Euclid will have a great potential for measuring the neutrino mass and excluding well-motivated scenarios with additional relativistic particles.

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

The steady-state gamma-ray emission from the Sun is thought to consist of two emission components due to interactions with Galactic cosmic rays: (1) a hadronic disk component, and (2) a leptonic extended component peaking at the solar edge and extending into the heliosphere. The flux of these components is expected to vary with the 11-year solar cycle, being highest during solar minimum and lowest during solar maximum, as it varies with the cosmic-ray flux. No study has yet analyzed the flux variation of each component over solar cycles. In this work, we measure the temporal variations of the flux of each component over 15 years of Fermi Large Area Telescope observations and compare them with the sunspot number and Galactic cosmic-ray flux from AMS-02 near Earth. We find that the flux variation of the disk anticorrelates with the sunspot number and correlates with cosmic-ray protons, as expected, confirming its emission mechanism. In contrast, the extended component exhibits a more complex variation: despite an initial anticorrelation with the sunspot number, we find neither anticorrelation with the sunspot number nor correlation with cosmic-ray electrons over the full 15-year period. This most likely suggests that cosmic-ray transport and modulation in the inner heliosphere are unexpectedly complex and may differ for electrons and protons or, alternatively, that there is an additional, unknown component of gamma rays or cosmic rays. These findings impact space weather research and emphasize the need for close monitoring of Cycle 25 and the ongoing polarity reversal.

El Gordo (ACT-CL J0102$-$4915) is a massive galaxy cluster with two major mass components at redshift $z=0.87$. Using SED fitting results from JWST/NIRCam photometry, the fraction of quenched galaxies in this cluster was measured in two bins of stellar mass: 9&lt;\log{({M_*}/\mathrm{M}_{\odot})}&lt;10 and 10\leq\log{({M_*}/\mathrm{M}_{\odot})}&lt;12. While there is no correlation between the quenched fraction and angular separation from the cluster's overall center of mass, there is a correlation between the quenched fraction and angular separation from the center of the nearest of the two mass components for the less-massive galaxies. This suggests that environmental quenching processes are in place at $z\sim1$, and that dwarf galaxies are more affected by those processes than massive galaxies.

We present new K-band spectroscopy for the giant elliptical galaxy M87 in the Virgo cluster, taken with the LUCI spectrograph at the Large Binocular Telescope (LBT). The new data are used to study line-strengths of K-band absorption features from different chemical species, namely Fe, Mg, Ca, Na, and CO, as a function of galactocentric distance, out to 40arcsec from the center (about half of the galaxy effective radius). The radial trends of spectral indices are compared to those for the bulge of M31, observed with the same instrument. For M87, most K-band indices exhibit flat radial profiles, with the exception of NaI2.21, that decreases outwards, with a negative radial gradient. Significant offsets are found between indices for M87 and those for the bulge of M31, the latter having weaker line-strengths for almost all features, but Fe and Ca, for which we find similar trends in both systems. We find that the behavior of CO features - most prominent in giant stars - is difficult to explain, consistent with previous results for the central regions of massive galaxies. In particular, the CO indices are stronger in M87 than M31, and do not exhibit significant radial gradients in M87, despite its IMF being bottom heavier than M31 especially in its central region. Predictions of state-of-the-art stellar population models, based on results from the optical spectral range, are able to match only the Na and Ca indices of M87, while a significant mismatch is found for all other indices. This shows that state-of-the-art stellar population models should be improved significantly in order to provide reliable constraints on the stellar population content of galaxies in the NIR spectral range.

The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10-year milestone. LAT performance remains well within the specifications defined during the planning phase, validating the design choices and supporting the compelling case to extend the duration of the Fermi mission. The details provided here will be useful when designing the next generation of high-energy gamma-ray observatories.

We present a new catalogue of ~2,400 optically selected quasars with spectroscopic redshifts and X-ray observations from either Chandra or XMM-Newton. The sample can be used to investigate the non-linear relation between the UV and X-ray luminosity of quasars, and to build a Hubble diagram up to redshift z~7.5. We selected sources that are neither reddened by dust in the optical/UV nor obscured by gas in the X-rays, and whose X-ray fluxes are free from flux-limit related biases. After checking for any possible systematics, we confirm, in agreement with our previous works, that (i) the X-ray to UV relation provides distance estimates matching those from supernovae up to z~1.5, and (ii) its slope shows no redshift evolution up to z~5. We provide a full description of the methodology for testing cosmological models, further supporting a trend whereby the Hubble diagram of quasars is well reproduced by the standard flat $\Lambda$CDM model up to z~1.5-2, but strong deviations emerge at higher redshifts. Since we have minimized all non-negligible systematic effects, and proven the stability of the $L_{\rm X}-L_{\rm UV}$ relation at high redshifts, we conclude that an evolution of the expansion rate of the Universe should be considered as a possible explanation for the observed deviation, rather than some systematic (redshift-dependent) effect associated with high-redshift quasars.

ESA's Euclid cosmology mission relies on the very sensitive and accurately calibrated spectroscopy channel of the Near-Infrared Spectrometer and Photometer (NISP). With three operational grisms in two wavelength intervals, NISP provides diffraction-limited slitless spectroscopy over a field of $0.57$ deg$^2$. A blue grism $\text{BG}_\text{E}$ covers the wavelength range $926$--$1366$\,nm at a spectral resolution $R=440$--$900$ for a $0.5''$ diameter source with a dispersion of $1.24$ nm px$^{-1}$. Two red grisms $\text{RG}_\text{E}$ span $1206$ to $1892$\,nm at $R=550$--$740$ and a dispersion of $1.37$ nm px$^{-1}$. We describe the construction of the grisms as well as the ground testing of the flight model of the NISP instrument where these properties were established.

Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Supernova Remnants (SNRs), IC 443, W44, W49B and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 years of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS~I +61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites.

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

In October 2007 a hard X-ray burst was detected by the INTEGRAL satellite from a direction consistent with the position of AX J1818.8-1559, an X-ray source at low Galactic latitude discovered with the ASCA satellite in 1996-1999. The short duration (0.8 s) and soft spectrum (power law photon index of 3.0+/-0.2) of the burst in the 20-100 keV range are typical of Soft Gamma-ray Repeaters and Anomalous X-ray Pulsars. We report on the results of an observation of AX J1818.8-1559 obtained with the Suzaku satellite in October 2011. The source spectrum, a power law with photon index 1.5, and flux 2x10^{-12} erg cm^-2 s^-1 (2-10 keV), do not show significant variations with respect to the values derived from archival data of various satellites (ROSAT, XMM-Newton, Chandra, Swift) obtained from 1993 to 2011. We discuss possible interpretations for AX J1818.8-1559 and, based on its association with the INTEGRAL burst, we propose it as a new member of the small class of magnetar candidates.

We present SIPGI, a spectroscopic pipeline to reduce optical/near-infrared data from slit-based spectrographs. SIPGI is a complete spectroscopic data reduction environment which retains the high level of flexibility and accuracy typical of the standard "by-hand" reduction methods but is characterized by a significantly higher level of efficiency. This is obtained by exploiting three main concepts: $i)$ the instrument model: at the core of the data reduction is an analytic description of the main calibration relations (e.g. spectra location and wavelength calibration) that can be easily checked and adjusted on data using a graphical tool; $ii)$ a built-in data organizer that classifies the data, together with a graphical interface that helps in providing the recipes with the correct input; $iii)$ the design and flexibility of the reduction recipes: the number of tasks required to perform a complete reduction is minimized, while preserving the possibility of verifying the accuracy of the main stages of data-reduction process with provided tools. The current version of SIPGI manages data from the MODS and LUCI spectrographs mounted at the Large Binocular Telescope, and it is our plan to extend SIPGI to support other through-slit spectrographs. Meanwhile, to allow using the same approach based on the instrument model with other instruments, we have developed SpectraPy, a spectrograph independent Python library working on through-slit spectra. In its current version, SpectraPy produces two-dimensional wavelength calibrated spectra corrected by instrument distortions. The current release of SIPGI and its documentation can by downloaded from this http URL, while SpectraPy can be found at this http URL.

GRB 221009A is the brightest gamma-ray burst (GRB) observed to date. Extensive observations of its afterglow emission across the electromagnetic spectrum were performed, providing the first strong evidence of a jet with a nontrivial angular structure in a long GRB. We carried out an extensive observation campaign in very-high-energy (VHE) gamma rays with the first Large-Sized Telescope (LST-1) of the future Cherenkov Telescope Array Observatory (CTAO), starting on 2022 October 10, about one day after the burst. A dedicated analysis of the GRB 221009A data is performed to account for the different moonlight conditions under which data were recorded. We find an excess of gamma-like events with a statistical significance of 4.1$\sigma$ during the observations taken 1.33 days after the burst, followed by background-compatible results for the later days. The results are compared with various models of afterglows from structured jets that are consistent with the published multiwavelength data, but entail significant quantitative and qualitative differences in the VHE emission after one day. We disfavor models that imply VHE flux at one day considerably above $10^{-11}$ erg cm$^{-2}$ s$^{-1}$. Our late-time VHE observations can help disentangle the degeneracy among the models and provide valuable new insight into the structure of GRB jets.

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged]

LensMC is a weak lensing shear measurement method developed for Euclid and Stage-IV surveys. It is based on forward modelling in order to deal with convolution by a point spread function (PSF) with comparable size to many galaxies; sampling the posterior distribution of galaxy parameters via Markov Chain Monte Carlo; and marginalisation over nuisance parameters for each of the 1.5 billion galaxies observed by Euclid. We quantified the scientific performance through high-fidelity images based on the Euclid Flagship simulations and emulation of the Euclid VIS images; realistic clustering with a mean surface number density of 250 arcmin$^{-2}$ (I_{\rm E}&lt;29.5) for galaxies, and 6 arcmin$^{-2}$ (I_{\rm E}&lt;26) for stars; and a diffraction-limited chromatic PSF with a full width at half maximum of $0.^{\!\prime\prime}2$ and spatial variation across the field of view. LensMC measured objects with a density of 90 arcmin$^{-2}$ (I_{\rm E}&lt;26.5) in 4500 deg$^2$. The total shear bias was broken down into measurement (our main focus here) and selection effects (which will be addressed elsewhere). We found measurement multiplicative and additive biases of $m_1=(-3.6\pm0.2)\times10^{-3}$, $m_2=(-4.3\pm0.2)\times10^{-3}$, $c_1=(-1.78\pm0.03)\times10^{-4}$, $c_2=(0.09\pm0.03)\times10^{-4}$; a large detection bias with a multiplicative component of $1.2\times10^{-2}$ and an additive component of $-3\times10^{-4}$; and a measurement PSF leakage of $\alpha_1=(-9\pm3)\times10^{-4}$ and $\alpha_2=(2\pm3)\times10^{-4}$. When model bias is suppressed, the obtained measurement biases are close to Euclid requirement and largely dominated by undetected faint galaxies ($-5\times10^{-3}$). Although significant, model bias will be straightforward to calibrate given the weak sensitivity. LensMC is publicly available at this https URL

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.

Accretion of matter onto a magnetic, rotating object can be strongly affected by the interaction with its magnetic field. This occurs in a variety of astrophysical settings involving young stellar objects, white dwarfs, and neutron stars. As matter is endowed with angular momentum, its inflow toward the star is often mediated by an accretion disc. The pressure of matter and that originating from the stellar magnetic field balance at the magnetospheric radius: at smaller distances the motion of matter is dominated by the magnetic field, and funnelling towards the magnetic poles ensues. However, if the star, and thus its magnetosphere, is fast spinning, most of the inflowing matter will be halted at the magnetospheric radius by centrifugal forces, resulting in a characteristic reduction of the accretion luminosity. The onset of this mechanism, called the propeller, has been widely adopted to interpret a distinctive knee in the decaying phase of the light curve of several transiently accreting X-ray pulsar systems. By comparing the observed luminosity at the knee for different classes of objects with the value predicted by accretion theory on the basis of the independently measured magnetic field, spin-period, mass, and radius of the star, we disclose here a general relation for the onset of the propeller which spans about eight orders of magnitude in spin period and ten in magnetic moment. The parameter-dependence and normalisation constant that we determine are in agreement with basic accretion theory.

The "Magnificent Seven" (M7) are a group of radio-quiet Isolated Neutron Stars (INSs) discovered in the soft X-rays through their purely thermal surface emission. Owing to the large inferred magnetic fields ($B\approx 10^{13}$ G), radiation from these sources is expected to be substantially polarised, independently on the mechanism actually responsible for the thermal emission. A large observed polarisation degree is, however, expected only if quantum-electrodynamics (QED) polarisation effects are present in the magnetised vacuum around the star. The detection of a strongly linearly polarised signal would therefore provide the first observational evidence of QED effects in the strong-field regime. While polarisation measurements in the soft X-rays are not feasible yet, optical polarisation measurements are within reach also for quite faint targets, like the M7 which have optical counterparts with magnitudes $\approx 26$--$28$. Here, we report on the measurement of optical linear polarisation for the prototype, and brightest member, of the class, RX\, J1856.5$-$3754 ($V\sim 25.5$), the first ever for one of the M7, obtained with the Very Large Telescope. We measured a polarisation degree $\mathrm{P.D.} =16.43\% \pm5.26\%$ and a polarisation position angle $\mathrm{P.A.}=145\fdg39\pm9\fdg44$, computed east of the North Celestial Meridian. The $\mathrm{P.D.}$ that we derive is large enough to support the presence of vacuum birefringence, as predicted by QED.