JWST has revealed an abundance of supermassive black holes (BHs) in the early Universe, and yet the lowest mass seed black holes that gave rise to these populations remain elusive. Here we present a systematic search for broad-line Active Galactic Nuclei (AGNs) in some of the faintest high-$z$ galaxies surveyed yet by combining ultra-deep JWST/NIRSpec G395M spectroscopy with the strong lensing aid in Abell S1063. By employing the profile of the [OIII]$\lambda 5007$ emission lines as a template for narrow-line components and carefully cross-validating with mock observations, we identify a sample of ten broad-line AGNs at $4.5

The detection of strong Balmer breaks and absorption features in Little Red Dots (LRDs) suggests they host AGN embedded within dense gas envelopes, potentially powered by super-Eddington accretion. We present GLIMPSE-17775, a luminous ($L_{\rm bol}\sim10^{45}$ erg s$^{-1}$) LRD at $z=3.501$ behind Abell S1063 ($\mu\sim2$), observed with deep JWST/NIRCam and a $\sim$20 hr (80 hr de-lensed) NIRSpec/G395M spectrum. The data reveal 40+ emission and absorption features, including a rich forest of low-ionization FeII lines and numerous broad hydrogen recombination transitions. We use this depth to test the dense-gas interpretation through five independent diagnostics. Nearly all permitted lines show exponential wings with consistent FWHM, the signature of Thomson scattering requiring $n_e\gtrsim10^8$ cm$^{-3}$. Adopting this width yields $M_{\rm BH}\sim10^{6.7}M_\odot$, a factor of ten lower than Gaussian fits, and $\lambda_{\rm Edd}\sim1.8$. Additional diagnostics support the same picture: a pronounced Balmer break ($f_{\nu,4050}/f_{\nu,3670}=2.0\pm0.1$), enhanced HeI $\lambda7065$ and $\lambda10830$ with P-Cygni absorption, Bowen-fluorescent OI $\lambda8446$-$\lambda11290$ emission requiring Ly$\beta$ pumping, and 16 FeII lines matching fluorescence models. These features indicate a dense ($n\sim10^8$ cm$^{-3}$), partially ionized cocoon where scattering and fluorescence dominate line formation, providing strong evidence that at least some LRDs are powered by super-Eddington black-hole growth in the early Universe.

We present an overview of the JWST GLIMPSE program, highlighting its survey design, primary science goals, gravitational lensing models, and first results. GLIMPSE provides ultra-deep JWST/NIRCam imaging across seven broadband filters (F090W, F115W, F200W, F277W, F356W, F444W) and two medium-band filters (F410M, F480M), with exposure times ranging from 20 to 40 hours per filter. This yields a 5$\sigma$ limiting magnitude of 30.9 AB (measured in a 0.2 arcsec diameter aperture). The field is supported by extensive ancillary data, including deep HST imaging from the Hubble Frontier Fields program, VLT/MUSE spectroscopy, and deep JWST/NIRSpec medium-resolution multi-object spectroscopy. Exploiting the strong gravitational lensing of the galaxy cluster Abell S1063, GLIMPSE probes intrinsic depths beyond 33 AB magnitudes and covers an effective source-plane area of approximately 4.4 arcmin$^2$ at $z \sim 6$. The program's central aim is to constrain the abundance of the faintest galaxies from $z \sim 6$ up to the highest redshifts, providing crucial benchmarks for galaxy formation models, which have so far been tested primarily on relatively bright systems. We present an initial sample of $\sim 540$ galaxy candidates identified at 6 < z < 16, with intrinsic UV magnitudes spanning $M_{\mathrm UV}$ = $-$20 to $-$12. This enables unprecedented constraints on the extreme faint end of the UV luminosity function at these epochs. In addition, GLIMPSE opens new windows for spatially resolved studies of star clusters in early galaxies and the detection and characterization of faint high-$z$ active galactic nuclei. This paper accompanies the first public data release, which includes reduced JWST and HST mosaics, photometric catalogs, and gravitational lensing models.

One of the first exoplanet hosts discovered thirty years ago, the star 55 Cnc has been constantly observed ever since. It is now known to host at least five planets with orbital periods ranging from 17 hours to 15 years. It is also one of the most extreme metal rich stars in the neighbourhood and it has a low-mass secondary star. In this article, we present data obtained at the Canada-France-Hawai'i Telescope with the SPIRou spectropolarimeter on both components of the 55 Cnc stellar system. We revisit the long-period radial-velocity signals of 55 Cnc A, with a focus on the role of the magnetic cycle, and propose the existence of a sixth planet candidate, whose period falls close to that of the magnetic cycle, or half of it. The other massive outer planet has a revised period of 13.15 years and a minimum mass of 3.8 MJup. Although some uncertainty remains on these outer planets, the characterization of the four inner planets is very robust through the combination of many different data sets, and all signals are consistent in the nIR and optical domains. In addition, the magnetic topology of the solar-type primary component of the system is observed by SPIRou at the minimum of its activity cycle, characterized by an amplitude ten times smaller than observed during its maximum in 2017. For the low-mass component 55 Cnc B, we report the discovery of two exoplanets in the system, with a period of 6.799+-0.0014 and 33.75+-0.04 days and a minimum mass of 3.5+-0.8 and 5.3+-1.4 MEarth, respectively. The secondary magnetic field is very weak and the current data set does not allow its precise characterization, setting an upper limit of 10 G. The system 55 Cnc stands out as the sixth binary system with planetary systems around both components, and the first one with non equal-mass stellar components.

The formation process of galaxy groups is not yet fully understood. In particular, that of fossil groups (FGs) is still under debate. Due to the relative rarity of FGs, large samples of such objects are still missing. The present paper aims to analyse the properties of groups in various evolutionary stages (FGs, "almost" FGs, and non-FGs), and to increase the sample of FG candidates. We have spectroscopically observed galaxies in four groups and ten candidate FGs detected in the Canada France Hawaii Telescope Legacy Survey. We searched for substructures by applying the Serna-Gerbal dendrogram method to analyse the dynamical structure of each group. By applying the FIREFLY software to the continuum and PIPE_VIS to the emission lines, we derived the stellar population properties in various regions for each group. A roughly continuous variation in properties is found between a group that is still building up (XCLASS 1330), a well-formed massive group (MCG+00-27-023), a dynamically complex non-FG (NGC 4065), and a near-FG (NGC 4104). We also optically confirm two FGs in the Canada France Hawaii Telescope Legacy Survey, but their X-ray luminosity is still unknown. We observe that the lower the mass of the substructure, the more recent the stellar population in the considered groups. We also show an apparent lack of high-mass substructures for low-metallicity systems. These results are consistent with the generally adopted model of energy transfer during interactions of the galaxies with the group and cluster potential wells. Furthermore, the fossil status of a group might be related to the large-scale environment. Therefore, studying the positions of non-FGs, near-FGs, and FGs within the cosmic web can provide insights into the process of how fossil systems come into being in the Universe.

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

Cosmic microwave background (CMB) photons are deflected by large-scale structure through gravitational lensing. This secondary effect introduces higher-order correlations in CMB anisotropies, which are used to reconstruct lensing deflections. This allows mapping of the integrated matter distribution along the line of sight, probing the growth of structure, and recovering an undistorted view of the last-scattering surface. Gravitational lensing has been measured by previous CMB experiments, with $\textit{Planck}$'s $42\,\sigma$ detection being the current best full-sky lensing map. We present an enhanced $\textit{LiteBIRD}$ lensing map by extending the CMB multipole range and including the minimum-variance estimation, leading to a $49$ to $58\,\sigma$ detection over $80\,\%$ of the sky, depending on the final complexity of polarized Galactic emission. The combination of $\textit{Planck}$ and $\textit{LiteBIRD}$ will be the best full-sky lensing map in the 2030s, providing a $72$ to $78\,\sigma$ detection over $80\,\%$ of the sky, almost doubling $\textit{Planck}$'s sensitivity. Finally, we explore different applications of the lensing map, including cosmological parameter estimation using a lensing-only likelihood and internal delensing, showing that the combination of both experiments leads to improved constraints. The combination of $\textit{Planck}$ + $\textit{LiteBIRD}$ will improve the $S_8$ constraint by a factor of 2 compared to $\textit{Planck}$, and $\textit{Planck}$ + $\textit{LiteBIRD}$ internal delensing will improve $\textit{LiteBIRD}$'s tensor-to-scalar ratio constraint by $6\,\%$. We have tested the robustness of our results against foreground models of different complexity, showing that improvements remains even for the most complex foregrounds.

Context: The relation between Galactic cosmic-ray electrons, magnetic fields and synchrotron radiation. Aims: We exploit synchrotron radiation to constrain the low-energy interstellar electron spectrum, using various radio surveys and connecting with electron data from Fermi-LAT and other experiments. Methods: The GALPROP programme for cosmic-ray propagation, gamma-ray and synchrotron radiation is used. Secondary electrons and positrons are included. Propagation models based on cosmic-ray and gamma-ray data are tested against synchrotron data from 22 MHz to 94 GHz. Results: The synchrotron data confirm the need for a low-energy break in the cosmic-ray electron injection spectrum. The interstellar spectrum below a few GeV has to be lower than standard models predict, and this suggests less solar modulation than usually assumed. Reacceleration models are more difficult to reconcile with the synchrotron constraints. We show that secondary leptons are important for the interpretation of synchrotron emission. We also consider a cosmic-ray propagation origin for the low-energy break. Conclusions: Exploiting the complementary information on cosmic rays and synchrotron gives unique and essential constraints on electrons, and has implications for gamma rays. This connection is especially relevant now in view of the ongoing PLANCK and Fermi missions.

We present a joint cosmological analysis combining data from the Planck satellite, the Atacama Cosmology Telescope, and the South Pole Telescope, constructing a unified likelihood that reproduces the measured temperature and polarisation power spectra by jointly modelling the cosmic microwave background (CMB) signal, galactic and extragalactic foregrounds, and instrumental systematics across all datasets. This approach reduces reliance on external priors and improves the robustness of parameter estimation. Within this joint analysis, $\Lambda$CDM parameters exhibit remarkable stability with respect to variations in foreground modelling. Extended cosmological parameters are more sensitive to these assumptions, with uncertainties increasing by up to 35%. Despite this, the combined constraints show no significant deviation from $\Lambda$CDM expectations, and several previously reported tensions -- such as the preference for non-zero curvature or the excess of lensing amplitude A_L -- are significantly reduced or resolved. In contrast, the determination of foreground parameters more severely depends on the assumptions made about the underlying models. Overall, this work demonstrates the feasibility and reliability of a fully joint analysis of current CMB experiments, and emphasizes the importance of consistent and accurate foreground modelling for the scientific goals of next-generation, high-sensitivity CMB surveys.

In this paper, we present extended gas kinematic maps of the Perseus cluster by combining five new XRISM/Resolve pointings observed in 2025 with four Performance Verification datasets from 2024, totaling 745 ks net exposure. To date, Perseus remains the only cluster that has been extensively mapped out to ~0.7$r_{2500}$ by XRISM/Resolve, while simultaneously offering sufficient spatial resolution to resolve gaseous substructures driven by mergers and AGN feedback. Our observations cover multiple radial directions and a broad dynamical range, enabling us to characterize the intracluster medium kinematics up to the scale of ~500 kpc. In the measurements, we detect high velocity dispersions ($\simeq$300 km/s) in the eastern region of the cluster, corresponding to a nonthermal pressure fraction of $\simeq$7-13%. The velocity field outside the AGN-dominant region can be effectively described by a single, large-scale kinematic driver based on the velocity structure function, which statistically favors an energy injection scale of at least a few hundred kpc. The estimated turbulent dissipation energy is comparable to the gravitational potential energy released by a recent merger, implying a significant role of turbulent cascade in the merger energy conversion. In the bulk velocity field, we observe a dipole-like pattern along the east-west direction with an amplitude of $\simeq\pm$200-300 km/s, indicating rotational motions induced by the recent merger event. This feature constrains the viewing direction to ~30$^\circ$-50$^\circ$ relative to the normal of the merger plane. Our hydrodynamic simulations suggest that Perseus has experienced at least two energetic mergers since redshift z~1, the latest associated with the radio galaxy IC310. This study showcases exciting scientific opportunities for future missions with high-resolution spectroscopic capabilities (e.g., HUBS, LEM, and NewAthena).

We investigate the impact of AGN feedback, on the entropy and characteristic temperature measurements of galaxy groups detected in the SRG/eROSITA's first All-Sky Survey (eRASS1) to shed light on the characteristics of the feedback mechanisms. We analyze deeper eROSITA observations of 1178 galaxy groups detected in eRASS1. We divide the sample into 271 subsamples and extract average thermodynamic properties, including electron density, temperature, and entropy at three characteristic radii along with the integrated temperature by jointly analyzing X-ray images and spectra following a Bayesian approach. We present the tightest constraints on the impact of AGN feedback through our average entropy and characteristic temperature measurements of the largest group sample used in X-ray studies, incorporating major systematics in our analysis. We find that entropy shows an increasing trend with temperature in the form of a power-law-like relation at the higher intra-group medium temperatures, while for the low mass groups, a slight flattening is observed on the average entropy. Overall, the observed entropy measurements agree well with the earlier measurements in the literature. The comparisons with the state-of-the-art cosmological hydrodynamic simulations (MillenniumTNG, Magneticum, OWL simulations) after the applications of the selection function calibrated for our galaxy groups reveal that observed entropy profiles in the cores are below the predictions of simulations. At the mid-region, the entropy measurements agree well with the Magneticum simulations, whereas the predictions of MillenniumTNG and OWL simulations fall below observations. At the outskirts, the overall agreement between the observations and simulations improves, with Magneticum simulations reproducing the observations the best. Our measurements will pave the way for more realistic AGN feedback implementations in simulations.

The XMM-Newton X-ray observatory has played a prominent role in astrophysics, conducting precise and thorough observations of the X-ray sky for the past two decades. The most recent iteration of the 4XMM catalogue and one of its latest data releases DR11 mark significant improvements over previous XMM-Newton catalogues, serving as a cornerstone for comprehending the diverse inhabitants of the X-ray sky. We employ detections and spectra extracted from the 4XMM-DR11 catalogue, subjecting them to fitting procedures using simple models. Our study operates within the framework of the XMM2ATHENA project, which focuses on developing state-of-the-art methods that exploit existing XMM-Newton data. We introduce and publicly release four catalogues containing measurements derived from X-ray spectral modelling of sources. The first catalogue encompasses outcomes obtained by fitting an absorbed power law model to all the extracted spectra for individual detections within the 4XMM-DR11 dataset. The second catalogue presents results obtained by fitting both an absorbed power law and an absorbed blackbody model to all unique physical sources listed in the 4XMM-DR11s catalogue, which documents source detection results from overlapping XMM-Newton observations. For the third catalogue we use the five band count rates derived from the pipe line detection of X-ray sources to mimic low resolution spectra to get a rough estimate of the spectral shape (absorbed power-law) of all 4XMM-DR11 detections. In the fourth catalogue, we conduct spectral analyses for the subset of identified sources with extracted spectra, employing various models based on their classification into categories such as AGN, stars, X-ray binaries, and cataclysmic variables. The scientific potential of these catalogues is highlighted by discussing the capabilities of optical and mid-infrared colours for selecting absorbed AGN. (abridged)

JWST observations have revealed rare galaxies with UV spectra exhibiting intense lines of nitrogen, indicative of super-solar N/O abundances at low metallicity. To better understand these enigmatic objects and provide new constraints on proposed scenarios, we have undertaken a systematic search for galaxies with UV emission lines of nitrogen. Using public JWST NIRSpec data, we have identified 45 N-emitters with robust NIII] or NIV] detections, including 4 previously known objects. We find N-emitters from redshift $z\sim 3-11$ among a broad diversity of galaxies, in terms of morphology, UV magnitude, stellar mass, SFR, metallicity, and rest-optical line strengths. The UV nitrogen lines show typical equivalent widths between $\sim 5-50$ Å. Carbon lines are generally fainter than the N lines. Using strong line calibrations established at high-redshift, we find metallicities $12+log(O/H)\sim 7.15-8.5$, including thus also high metallicities. The H$\beta$ equivalent width of N-emitters varies strongly, and sources with low EWs show clear signs of a Balmer break, indicative of composite stellar populations combining both young (< 10 Myr) stars responsible of the UV emission lines and an older population contributing to the rest-optical spectrum. Supersolar N/O ratios are found in all N-emitters. C/O abundances are comparable to those of galaxies at the same metallicity, and all N-emitters show high N/C ratios or lower limits (\log(N/C) &gt; 0.5), independently of metallicity. The observed abundance ratios are compatible with ejecta from H-burning and do not show signs of Carbon enhancements, even at higher metallicities. Finally, we find that the fraction of N-emitters increases with redshift, and we quantify this evolution. Our study increases the sample of known N-emitters by a factor $\sim 3$, reveals a diversity of properties among N-emitters, and provides new constraints on their nature.

The magnetorotational (MRI) dynamo has long been considered one of the possible drivers of turbulent angular momentum transport in astrophysical accretion disks. However, various numerical results suggest that this dynamo may be difficult to excite in the astrophysically relevant regime of magnetic Prandtl number (Pm) significantly smaller than unity, for reasons currently not well understood. The aim of this article is to present the first results of an ongoing numerical investigation of the role of both linear and nonlinear dissipative effects in this problem. Combining a parametric exploration and an energy analysis of incompressible nonlinear MRI dynamo cycles representative of the transitional dynamics in large aspect ratio shearing boxes, we find that turbulent magnetic diffusion makes the excitation and sustainment of this dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for decreasing Pm. This results in an increase in the critical Rm of the dynamo for increasing kinematic Reynolds number (Re), in agreement with earlier numerical results. Given its very generic nature, we argue that turbulent magnetic diffusion could be an important determinant of MRI dynamo excitation in disks, and may also limit the efficiency of angular momentum transport by MRI turbulence in low Pm regimes.

We present a comprehensive analysis of the presence of very massive stars (VMS > $100 M_{\odot}$) in the integrated spectra of 13 UV-bright star-forming galaxies at $2.2 \lesssim z \lesssim 3.6$ taken with the Gran Telescopio Canarias (GTC). These galaxies have very high UV absolute magnitudes ($M_{\rm UV} \simeq -24$), intense star formation (SFR $\simeq 100-1000$ $M_{\odot}$ yr$^{-1}$), and metallicities in the range of 12+log(O/H) $\simeq8.10-8.50$ inferred from strong rest-optical lines. The GTC rest-UV spectra reveal spectral features indicative of very young stellar populations with VMS, such as strong P-Cygni line profiles in the wind lines N~{\sc v} $\lambda 1240$ and C~{\sc iv} $\lambda 1550$ along with intense and broad He~{\sc ii} $\lambda 1640$ emission with $EW_{0}$ $\simeq 1.40-4.60$ Å, and FWHM $\simeq 1150-3170$ $km \ s^{-1}$. A Comparison with known VMS-dominated sources and typical galaxies without VMS reveals that some UV-bright galaxies closely resemble VMS-dominated clusters (e.g., R136 cluster). The presence of VMS is further supported by a quantitative comparison of the observed strength of the He~{\sc ii} emission with population synthesis models with and without VMS, where models with VMS are clearly preferred. Employing an empirical threshold for $EW_{0}$ (\heii) $\geq 3.0$ Å, along with the detection of other VMS-related spectral profiles (N~{\sc iv} $\lambda 1486, 1719$), we classify nine out of 13 UV-bright galaxies as VMS-dominated sources. This high incidence of VMS-dominated sources in the UV-bright galaxy population ($\approx 70\%$) contrasts significantly with the negligible presence of VMS in typical $L_{\rm UV}^{*}$ LBGs at similar redshifts (&lt;1\%). Our results thus indicate that VMS are common in UV-bright galaxies, suggesting a different initial mass function (IMF) with upper mass limits between $175 M_{\odot}$ and $475 M_{\odot}$.

We present ionizing photon production efficiencies (${\xi}_{\rm ion}$) for 63 z=1.5-6.9 star-forming galaxies using precise nebular dust attenuation corrections from the JWST/AURORA survey. A subset of objects within AURORA have individually-determined nebular dust attenuation curves, which vary significantly in shape and normalization, resulting in reduced systematic uncertainty when constraining the total attenuation of H${\alpha}$ luminosity, and thus the intrinsic ionizing output within our sample. We find evidence for positive correlations between ${\xi}_{\rm ion}$ and redshift, equivalent width of [OIII]${\lambda}$5007, and O32=[OIII]${\lambda}$5007/[OII]${\lambda}$3726,3729, and negative correlations between ${\xi}_{\rm ion}$ and stellar attenuation, UV luminosity (L$_{\rm UV}$), stellar mass, and direct-method metallicity. We test alternate dust prescriptions within this sample, and find that the total attenuation is lower when using the commonly-assumed Galactic extinction curve or when assuming that stellar attenuation is equal to nebular attenuation. We also find that assuming either of these alternate dust prescriptions can change the slope of relationships between ${\xi}_{\rm ion}$ and galaxy property, notably inducing a flat trend between ${\xi}_{\rm ion}$ and L$_{\rm UV}$ within AURORA. While the novel nebular dust curves derived from AURORA spectroscopy reveal obscured ionizing photon production within star-forming galaxies at these redshifts, a more complete understanding of stellar attenuation is required to fully reduce dust systematics on ${\xi}_{\rm ion}$ for inclusion in reionization models.

Herschel/HIFI observations toward the compact HII region W51 has revealed the presence of a cold dense core along its line of sight in a high-velocity stream located just in front of W51. This detection has been made possible through absorption measurements of low-energy transitions of HDO, NH3, and C3 against the bright background emitted by the star-forming region. We present a follow-up study of this core using the high sensitivity and high spectral resolution provided by the IRAM 30-meter telescope. We report new detections of this core in absorption for DCO+ (2-1, 3-2), H13CO+ (1-0), DNC (3-2), HN13C (1-0), p-H2CO (202-101, 303-202), and in emission for o-NH2D. We also report interferometric observation of this last species using the IRAM/NOEMA telescope, revealing the fragmented nature of the source through the detection of two cores, separated by 0.19 - 0.24 pc, with average sizes of less than 0.16 - 0.19 pc. From a non-LTE analysis, we are able to estimate the density (~2.5e4 cm-3) and temperature (~10 K) of this component, typical of what is found in dark clouds. This component (called W51-core) has the same DCO+/HCO+ ratio (0.02) as TMC-1 and a high DNC/HNC ratio (0.14). Detection of these deuterated species indicates that W51-core is similar to an early-phase low-mass star-forming region, formed from the interaction between the W51 giant molecular cloud and the high-velocity stream in front of it. The W51 complex being at about 5 kpc, these findings lead to what is the first detection of the earliest phase of low-mass star-forming region at such a large distance.

Cherenkov Telescope Array Observatory (CTAO) is the next-generation ground-based gamma-ray observatory operating in the energy range from 20 GeV up to 300 TeV, with two sites in La Palma (Spain) and Paranal (Chile). It will consist of telescopes of three sizes, covering different parts of the large energy range. We report on the performance of Large-Sized Telescope prototype (LST-1) in the detection and characterization of extragalactic gamma-ray sources, with a focus on the reconstructed gamma-ray spectra and variability of classical bright BL Lacertae objects, which were observed during the early commissioning phase of the instrument. LST-1 data from known bright gamma-ray blazars - Markarian 421, Markarian 501, 1ES 1959+650, 1ES 0647+250, and PG 1553+113 - were collected between July 10, 2020, and May 23, 2022, covering a zenith angle range of 4 deg to 57 deg. The reconstructed light curves were analyzed using a Bayesian block algorithm to distinguish the different activity phases of each blazar. Simultaneous Fermi-LAT data were utilized to reconstruct the broadband $\gamma$-ray spectra for the sources during each activity phase. High-level reconstructed data in a format compatible with gammapy are provided together with measured light curves and spectral energy distributions (SEDs) for several bright blazars and an interpretation of the observed variability in long and short timescales. Simulations of historical flares are generated to evaluate the sensitivity of LST-1. This work represents the first milestone in monitoring bright BL Lacertae objects with a CTAO telescope.