Accurate pulsar astrometric estimates play an essential role in almost all high-precision pulsar timing experiments. Traditional pulsar timing techniques refine these estimates by including them as free parameters when fitting a model to observed pulse time-of-arrival measurements. However, reliable sub-milliarcsecond astrometric estimations require years of observations and, even then, power from red noise can be inadvertently absorbed into astrometric parameter fits, biasing the resulting estimations and reducing our sensitivity to red noise processes, including gravitational waves (GWs). In this work, we seek to mitigate these shortcomings by using pulsar astrometric estimates derived from Very Long Baseline Interferometry (VLBI) as priors for the timing fit. First, we calibrated a frame tie to account for the offsets between the reference frames used in VLBI and timing. Then, we used the VLBI-informed priors and timing-based likelihoods of several astrometric solutions consistent with both techniques to obtain a maximum-posterior astrometric solution. We found offsets between our results and the timing-based astrometric solutions, which, if real, would lead to absorption of spectral power at frequencies of interest for single-source GW searches. However, we do not find significant power absorption due to astrometric fitting at the low-frequency domain of the GW background.

We present \textit{JWST} NIRSpec and MIRI MRS observations of the central kiloparsec of M58 (NGC 4579), a nearby LINER galaxy hosting a low-luminosity AGN (LLAGN; $L_\mathrm{bol} \sim 10^{42}$ erg s$^{-1}$) with a low-power jet. These data provide an unprecedented view of the warm molecular gas phase and reveal clear signatures of feedback. We detect 44 H$_2$ lines, including bright pure rotational lines (S(1)-S(18)) and rovibrational lines up to $\nu=2$, probing a wide range of excitation conditions. Excitation diagrams show that rotational lines follow a power-law temperature distribution with an exponential cutoff, consistent with heating by low-velocity shocks. H$_2$ rovibrational lines deviate from thermal models primarily because of sub-thermal excitation at low density. Additionally, there may be a 10% contribution by AGN X-ray heating in the nucleus. The dust lanes associated with the spiral inflow appear dynamically undisturbed but show signs of shock heating, while the inner $\sim$200 pc exhibits turbulent kinematics produced by outflowing molecular gas. These results reveal the subtle yet measurable impact of LLAGN feedback on the interstellar medium, demonstrating that even weak, vertically oriented jets and low radiative accretion rates can perturb molecular gas and regulate nuclear reservoirs. This study highlights JWST's transformative ability to uncover hidden modes of AGN feedback.

The intrinsic properties of galaxies are influenced by their environments, underscoring the environment's critical role in galaxy formation and evolution. Traditionally, these environments are categorized into four fixed classifications: knots, filaments, walls, and voids, which collectively describe the complex organization of galaxies within large-scale structures. We propose an alternative description that complements the traditional quadripartite categorization by introducing a continuous framework, allowing for a more nuanced examination of the relationship between the intrinsic properties of galaxies and their environments. This complementary description is applied using one of the most prevalent methodologies: categorization using the eigenvalues of the Hessian matrix extracted from the matter density field. We integrated our findings into a semi-analytical model of galaxy formation, combined with cosmological numerical simulations, to analyze how the intrinsic properties of galaxies are influenced by environmental changes. In our study, we find a continuous distribution of eigenvalue ratios, revealing a clear dependence of galaxy properties on their surrounding environments. This method allowed us to identify critical values at which transitions in the behavior of key astrophysical galaxy properties become evident.

We present a framework for dark matter (DM) halo formation based on a kinetic theory of self-gravitating fermions together with a solid connection to thermodynamics. Based on maximum entropy arguments, this approach predicts a most likely phase-space distribution which takes into account the Pauli exclusion principle, relativistic effects, and particle evaporation. The most general equilibrium configurations depend on the particle mass and develop a degenerate compact core embedded in a diluted halo, both linked by their fermionic nature. By applying such a theory to the Milky Way we analyze the stability of different families of equilibrium solutions with implications on the DM distribution and the mass of the DM candidate. We find that stable core-halo profiles, which explain the DM distribution in the Galaxy, exist only in the range $mc^2 \approx 194 - 387\,\rm{keV}$. The lower bound is a consequence of imposing thermodynamical stability on the core-halo solutions having a $4.2\times 10^6 M_\odot$ quantum core mass alternative to the black hole hypothesis at the Galaxy center. The upper bound is solely an outcome of general relativity when the quantum core reaches the Oppenheimer-Volkoff limit and undergoes gravitational collapse towards a black hole. We demonstrate that there exists a set of stable core-halo profiles which are astrophysically relevant in the sense that their total mass is finite, do not suffer from the gravothermal catastrophe, and agree with observations. The morphology of the outer halo tail is described by a polytrope of index $5/2$, developing a sharp decline of the density beyond $25\,\rm{kpc}$ in excellent agreement with the latest Gaia DR3 rotation curve data. Moreover, we obtain a total mass of about $2\times 10^{11} M_\odot$ including baryons and a local DM density of about $0.4\,\rm{GeV}\,c^{-2}\,\rm{cm}^{-3}$ in line with recent independent estimates.

The radius valley separating super-Earths from mini-Neptunes is a fundamental benchmark for theories of planet formation and evolution. Observations show that the location of the radius valley decreases with decreasing stellar mass and with increasing orbital period. Here, we build from our previous pebble-based formation model, which, combined with photoevaporation after disc dispersal, unveiled the radius valley as a separator between rocky- and water-worlds. We expand our models for a range of stellar masses spanning from 0.1 to 1.5 $M_\odot$. We find that the location of the radius valley is well described by a power-law in stellar mass as $R_{\rm valley} = 1.8197 \, M_{\star}^{\!0.14({+0.02}/{-0.01})}$, which is in excellent agreement with observations. We also find very good agreement with the dependence of the radius valley on orbital period, both for FGK- and M-dwarfs. Additionally, we note that the radius valley gets filled towards low stellar masses, particularly at 0.1-0.4 $M_\odot$, yielding a rather flat slope in $R_{\rm valley} - P_{\rm orb}$. This is the result of orbital migration occurring at lower planet mass for less massive stars, which allows for low-mass water-worlds to reach the inner regions of the system, blurring the separation in mass (and size) between rocky- and water-worlds. Furthermore, we find that for planetary equilibrium temperatures above 400 K, the water in the volatile layer exists fully in the form of steam, puffing the planet radius up compared to condensed-water worlds. This produces an increase in planet radii of $\sim 30\%$at 1$M_\oplus$, and of$\sim 15\%$at 5$M_\oplus$, compared to condensed-water-worlds. As with Sun-like stars, pebble accretion leaves its imprint on the overall exoplanet population as a depletion of planets with intermediate compositions, carving a valley in planet density for all spectral types (abridged).

The recent detection in archival HST images of an object at the the location of supernova (SN) iPTF13bvn may represent the first direct evidence of the progenitor of a Type Ib SN. The object's photometry was found to be compatible with a Wolf-Rayet pre-SN star mass of ~11 Msun. However, based on hydrodynamical models we show that the progenitor had a pre-SN mass of ~3.5 Msun and that it could not be larger than ~8 Msun. We propose an interacting binary system as the SN progenitor and perform evolutionary calculations that are able to self-consistently explain the light-curve shape, the absence of hydrogen, and the pre-SN photometry. We further discuss the range of allowed binary systems and predict that the remaining companion is a luminous O-type star of significantly lower flux in the optical than the pre-SN object. A future detection of such star may be possible and would provide the first robust identification of a progenitor system for a Type Ib SN.

Although dust constitutes only about 1% of the mass of a protoplanetary disk, recent studies demonstrate that it can exert a significant torque on low- and intermediate-mass planetary cores. We compute and quantify for the first time the influence of the dust torque on the evolution of growing planetary embryos as they move in a protoplanetary disk while growing via gas and pebble accretion. Our global model evolves the gaseous disk via viscous accretion and X-ray photoevaporation, while accounting for dust growth and evolution including coagulation, drift, and fragmentation. Our research indicates that dust torque significantly influences planetary migration, particularly driving substantial outward migration for planets forming within the water ice-line. This effect occurs due to an increased dust-to-gas mass ratio in the inner disk, resulting from inward pebble drift from outer regions. In contrast, for planets initially located beyond the water ice-line, the dust torque mitigates inward migration but does not significantly alter their paths, as the dust-to-gas ratio diminishes rapidly due to rapid pebble drift and the brief timescales of planet formation in these areas. These findings underscore the pivotal role of dust torque in shaping the migration patterns of low- and intermediate-mass planets, especially when enhanced dust concentrations in the inner disk amplify its effects

Surrounding Sgr A*, a cluster of young and massive stars coexist with a population of dust-enshrouded objects, whose astrometric data can be used to scrutinize the nature of Sgr A*. An alternative to the black hole (BH) scenario has been recently proposed in terms of a supermassive compact object composed of self-gravitating fermionic dark matter (DM). Such horizon-less configurations can reproduce the relativistic effects measured for S2 orbit, while being part of a single continuous configuration whose extended halo reproduces the latest GAIA-DR3 rotation curve. In this work, we statistically compare different fermionic DM configurations aimed to fit the astrometric data of S2, and five G-sources, and compare with the BH potential when appropriate. We sample the parameter spaces via Markov Chain Monte Carlo statistics and perform a quantitative comparison estimating Bayes factors for models that share the same likelihood function. We extend previous results of the S2 and G2 orbital fits for 56 keV fermions (low core-compactness) and show the results for 300 keV fermions (high core-compactness). For the selected S2 dataset, the former model is slightly favoured over the latter. However, more precise S2 datasets, as obtained by the GRAVITY instrument, remain to be analysed in light of the fermionic models. For the G-objects, no conclusive preference emerges between models. For all stellar objects tested, the BH and fermionic models predict orbital parameters that differ by less than 1%. More accurate data, particularly from stars closer to Sgr A*, is necessary to statistically distinguish between the models considered.

For decades, the theoretical understanding of planetary nebulae (PNe) has remained in tension with the observed universal bright-end cutoff of the PN luminosity function (PNLF). While the brightest younger PN populations are expected to be brighter in their [OIII] emission than observed, recent studies have proposed circumnebular extinction to be a key ingredient for bringing their brightness down to the observed bright end. In this work we use the recently introduced PICS (PNe In Cosmological Simulations) framework to investigate the impact of different circumnebular extinction treatments on the modeled PNe and their PNLF for a large range of stellar ages and metallicities. We test how different slopes in the observed relation of extinction versus central star mass modify the bright-end cutoffs of the PNLF, finding that steeper slopes lead to large changes for young stellar populations. In contrast, the differences for older PNe are much smaller. However, for individual PNe, the extinctions observed in nearby galaxies appear to be much higher than the models predict, showing that improvements on both the modeling and observational sides are needed to gain a better understanding of the brightest and strongly extincted PNe. These findings further advance the theoretical foundation for interpreting observed extragalactic PN populations coming from more complex composite stellar populations in the future.

We investigate the structural properties, distribution and abundance of LCDM dark matter subhaloes using the Phi-4096 and Uchuu suite of N-body cosmological simulations. Thanks to the combination of their large volume, high mass resolution and superb statistics, we are able to quantify -- for the first time consistently over more than seven decades in ratio of subhalo-to-host-halo mass -- dependencies of subhalo properties with mass, maximum circular velocity, Vmax, host halo mass and distance to host halo centre. We also dissect the evolution of these dependencies over cosmic time. We provide accurate fits for the subhalo mass and velocity functions, both exhibiting decreasing power-law slopes in the expected range of values and with no significant dependence on redshift. We also find subhalo abundance to depend weakly on host halo mass. We explore the distribution of subhaloes within their hosts and its evolution over cosmic time for subhaloes located as deep as ~0.1 per cent of the host virial radius. Subhalo structural properties are codified via a concentration parameter, cV, that does not depend on any specific, pre-defined density profile and relies only on Vmax. We derive the cV-Vmax relation in the range 7-1500 km/s and find an important dependence on distance of the subhalo to the host halo centre, as already described in Moliné et al. (2017). Interestingly, we also find subhaloes of the same mass to be significantly more concentrated into more massive hosts. Finally, we investigate the redshift evolution of cV, and provide accurate fits that take into account all mentioned dependencies. Our results offer an unprecedented detailed characterization of the subhalo population, consistent over a wide range of subhalo and host halo masses, as well as cosmic times. Our work enables precision work in any future research involving dark matter halo substructure.

We present time-series photometry of 31 massive DA white dwarfs with $M\gtrsim 0.9~M_\odot$ within the ZZ Ceti instability strip from the Montreal White Dwarf Database 100 pc sample. The majority of the targets had no previous time-series photometry available, though several were classified as non-variable or potential pulsators in the literature. Out of the 31 candidates, we confirm 16 as pulsating. Our observations at three observatories have led us to discover the most massive pulsating white dwarf currently known, J0959$-$1828 ($M=1.32$ or $1.27~M_\odot$ for a CO versus ONe core), which is slightly more massive than the previous record holder J0049$-$2525. We study the sample properties of massive ZZ Ceti white dwarfs, and find several trends with their weighted mean periods. As predicted by theory, we see an increase in the weighted mean periods with decreasing effective temperature, and a decrease in pulsation amplitudes at the red edge of the instability strip. Furthermore, the weighted mean periods decrease with increasing stellar mass. Our observations show that the ZZ Ceti instability strip may not be pure at high masses. This is likely because the non-variable white dwarfs in the middle of the strip may be weakly magnetic, which could escape detection in the available low-resolution spectroscopy data, but may be sufficient to suppress pulsations. Extensive follow-up observations of the most massive white dwarfs in our sample have the potential to probe the interior structures and core-compositions of these white dwarfs with significantly crystallized cores.

We present a comprehensive study of the X-ray binary system XTE~J1550-564, with the primary objective of analyzing the evolution of the black hole's spin parameter. To achieve this objective, we embarked on the necessary step of identifying a plausible progenitor for the system. Using a set of models covering various parameter combinations, we were able to replicate the system's observed characteristics within acceptable error margins, including fundamental parameters such as component masses, orbital period, donor luminosity, and effective temperature. The model results indicate the possibility of diverse evolutionary pathways for the system, highlighting the significant role played by the initial mass of the donor star and the efficiency of mass transfer episodes. While some models are well-aligned with estimates of the mass transfer rate, they all fall short of explaining the black hole's observed moderate spin ($a^* = 0.49$). We also explored alternative magnetic braking prescriptions, finding that only an extreme and fully conservative scenario, based on the convection and rotation boosted prescription, can reproduce the observed spin and only in a marginal way. Our study attempts to shed light on the complex dynamics of black hole X-ray binaries and the challenges of explaining their observed properties with theoretical models.

We characterize the overdensity (spike) of fermionic dark matter (DM) particles around a supermassive black hole (SMBH) within a general relativistic analysis. The initial DM halo distribution is obtained by solving the equilibrium equations of a self-gravitating system of massive fermions at a finite temperature, according to the Ruffini-Argüelles-Rueda (RAR) model. The final fermionic DM spike is calculated around a Schwarzschild SMBH. We explore two possible interpretations for the origin and evolution of the SMBH seed. One corresponds to the traditional scenario proposed by Gondolo & Silk (1999), where a small BH mass of baryonic origin sits at the halo's center and grows adiabatically. The other one, from DM origin, where the dense and degenerate fermion core predicted by the RAR model grows adiabatically by capturing baryons until its gravitational collapse, providing a heavy SMBH seed, whose specific value depends on the fermion mass. We study different initial fermionic DM profiles that the theory allows. We show that overall dilute (i.e., Boltzmannian) fermionic DM develops the well-known spike with density profile $\rho \sim r^{-3/2}$. Instead, for semi-degenerate fermions with a dense and compact core surrounded by a diluted halo, we find novel spike profiles that depend on the particle mass and nature. In the more general case, fermionic spikes do not develop a simple power-law profile. Furthermore, the SMBH does not always imply an enhancement of the surrounding DM density; it might also deplete it. Thus, the self-consistent inclusion of the DM candidate nature and mass in determining the structure and distribution of DM in galaxies, including the DM spikes around SMBHs, is essential for the specification of DM astrophysical probes such as BH mergers, gravitational waves, or stellar orbits.

The metal content of a galaxy's interstellar medium reflects the interplay between different evolutionary processes such as feedback from massive stars and the accretion of gas from the intergalactic medium. Despite the expected abundance of low-luminosity galaxies, the low-mass and low-metallicity regime remains relatively understudied. Since the properties of their interstellar medium resemble those of early galaxies, identifying such objects in the Local Universe is crucial to understand the early stages of galaxy evolution. We used the DR3 catalog of the Southern Photometric Local Universe Survey (S-PLUS) to select low-metallicity dwarf galaxy candidates based on color selection criteria typical of metal-poor, star-forming, low-mass systems. The final sample contains approximately 50 candidates. Spectral energy distribution fitting of the 12 S-PLUS bands reveals that $\sim$ 60% of the candidates are best fit by models with low stellar metallicities. We obtained long-slit observations with the Gemini Multi-Object Spectrograph to follow-up a pilot sample and confirm whether these galaxies have low metallicities. We find oxygen abundances in the range 7.28< 12 + log(O/H) < 7.82 (4% to 13% of the solar value), confirming their metal-poor nature. Most targets are outliers in the mass-metallicity relation, i.e. they display a low metal content relative to their observed stellar masses. In some cases, perturbed optical morphologies might give evidence of dwarf-dwarf interactions or mergers. These results suggest that the low oxygen abundances may be associated with an external event causing the accretion of metal-poor gas, which dilutes the oxygen abundance in these systems.

The Red-Giant Branch Bump (RGBB) is one of the most noteworthy features in the red-giant luminosity function of stellar clusters. It is caused by the passage of the hydrogen-burning shell through the composition discontinuity left at the point of the deepest penetration by the convective envelope. When crossing the discontinuity the usual trend in increasing luminosity reverses for a short time before it increases again, causing a zig-zag in the evolutionary track. In spite of its apparent simplicity the actual physical reason behind the decrease in luminosity is not well understood and several different explanations have been offered. Here we use a recently proposed simple toy-model for the structure of low-mass red giants, together with previous results, to show beyond reasonable doubt that the change in luminosity at the RGBB can be traced to the change in the mean molecular weight of the layers on top of the burning shell. And that these changes happen on a nuclear timescale. The change in the effective mean molecular weight, as the burning shell approaches the discontinuity, causes a drop in the temperature of the burning shell which is attenuated by the consequent feedback contraction of the layers immediately below the burning shell. Our work shows that, when applied correctly, including the feedback on the structure of the core together with of the increase in the mass of the core, shell-source homology relations do a great quantitative job in explaining the properties of full evolutionary models at the RGBB.

It has been known for decades that the observed number of baryons in the local universe falls about 30-40% short of the total number of baryons predicted by Big-Bang Nucleosynthesis, as inferred from density fluctuations of the Cosmic Microwave Background and seen during the first 2-3 billion years of the universe in the so called Lyman-alpha Forest. A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas - hydrogen - is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios. Indeed, despite the large observational efforts, only a few marginal claims of detection have been made so far. Here we report observations of two absorbers of highly ionized oxygen (OVII) in the high signal-to-noise-ratio X-ray spectrum of a quasar at redshift >0.4. These absorbers show no variability over a 2-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The OVII systems lie in regions characterized by large (x4 compared to average) galaxy over-densities and their number (down to the sensitivity threshold of our data), agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium (WHIM). We conclude that the missing baryons have been found.

We present tests of a new method to simultaneously estimate stellar population and emission line (EL) properties of galaxies out of S-PLUS photometry. The technique uses the AlStar code, updated with an empirical prior which greatly improves its ability to estimate ELs using only the survey's 12 bands. The tests compare the output of (noise-perturbed) synthetic photometry of SDSS galaxies to properties derived from previous full spectral fitting and detailed EL analysis. For realistic signal-to-noise ratios, stellar population properties are recovered to better than 0.2 dex in masses, mean ages, metallicities and $\pm 0.2$ mag for the extinction. More importantly, ELs are recovered remarkably well for a photometric survey. We obtain input $-$ output dispersions of 0.05--0.2 dex for the equivalent widths of $[\mathrm{O}\,\rm{II}]$, $[\mathrm{O}\,\rm{III}]$, H$\beta$, H$\alpha$, $[\mathrm{N}\,\rm{II}]$, and $[\mathrm{S}\,\rm{II}]$, and even better for lines stronger than $\sim 5$ $\mathring{A}$. These excellent results are achieved by combining two empirical facts into a prior which restricts the EL space available for the fits: (1) Because, for the redshifts explored here, H$\alpha$ and $[\mathrm{N}\,\rm{II}]$ fall in a single narrow band (J0660), their combined equivalent width is always well recovered, even when $[\mathrm{N}\,\rm{II}]$/H$\alpha$ is not. (2) We know from SDSS that $W_{H\alpha+[\mathrm{N}\,\rm{II}]}$ correlates with $[\mathrm{N}\,\rm{II}]$/H$\alpha$, which can be used to tell if a galaxy belongs to the left or right wings in the classical BPT diagnostic diagram. Example applications to integrated light and spatially resolved data are also presented, including a comparison with independent results obtained with MUSE-based integral field spectroscopy.

We present observations of ASASSN-13dn, one of the first supernovae discovered by ASAS-SN, and a new member of the rare group of Luminous Type II Supernovae (LSNe II). It was discovered near maximum light, reaching an absolute magnitude of M$_{v}$ $\sim$ -19 mag, placing this object between normal luminosity type II SNe and superluminous SNe A detailed analysis of the photometric and spectroscopic data of ASASSN-13dn is performed. The spectra are characterized by broad lines, in particular the H$\alpha$ lines where we measure expansion velocities ranging between 14000 - 6000 km s$^{-1}$ over the first 100 days. H$\alpha$ dominates the nebular spectra, and we detect a narrow P-Cygni absorption within the broader emission line with an expansion velocity of 1100 km s$^{-1}$. Photometrically, its light curve shows a re-brightening of $\sim$ 0.6 mag in the $gri$ bands starting at 25$\pm$2 days after discovery, with a secondary peak at $\sim 73$d, followed by an abrupt and nearly linear decay of 0.09 mag d$^{-1}$ for the next 35 days. At later times, after a drop of 4 magnitudes from the second maximum, the light curves of ASASSN-13dn shows softer undulations from 125 to 175 days. We compare ASASSN-13dn with other LSNe II in the literature, finding no match to both light curve and spectroscopic properties. We discuss the main powering mechanism and suggest that interaction between the ejecta and a dense CSM produced by eruptions from an LBV-like progenitor could potentially explain the observations.

Ultra-massive hydrogen-rich (DA spectral type) white dwarf (WD) stars ($M_{\star} > 1M_{\odot}$) coming from single-star evolution are expected to harbor cores made of $^{16}$O and $^{20}$Ne, resulting from semi-degenerate carbon burning when the progenitor star evolves through the super asymptotic giant branch (S-AGB) phase. These stars are expected to be crystallized by the time they reach the ZZ Ceti instability strip ($T_{\rm eff} \sim 12\,500$ K). Theoretical models predict that crystallization leads to a separation of $^{16}$O and $^{20}$Ne in the core of ultra-massive WDs, which impacts their pulsational properties. This property offers a unique opportunity to study the processes of crystallization. Here, we present the first results of a detailed asteroseismic analysis of the best-studied ultra-massive ZZ Ceti star BPM~37093. As a second step, we plan to repeat this analysis using ultra-massive DA WD models with C/O cores in order to study the possibility of elucidating the core chemical composition of BPM~37093 and shed some light on its possible evolutionary origin. We also plan to extend this kind of analyses to other stars observed from the ground and also from space missions like Kepler and TESS.

Integral Field Spectroscopy (IFS) is well known for providing detailed insight of extended sources thanks to the possibility of handling space resolved spectroscopic information. Simple and straightforward analysis such as single line fitting yield interesting results, although it might miss a more complete picture in many cases. Violent star forming regions, such as starburst galaxies, display very complex emission line profiles due to multiple kinematic components superposed in the line of sight. We perform a spatially resolved kinematical study of a single Green Pea (GP) galaxy, SDSSJ083843.63+385350.5, using a new method for analyzing Integral Field Unit (IFU) observations of emission line spectra. The method considers the presence of multiple components in the emission-line profiles and makes use of a statistical indicator to determine the meaningful number of components to fit the observed profiles. We are able to identify three distinct kinematic features throughout the field and discuss their link with a rotating component, a strong outflow and a turbulent mixing layer. We also derive an updated star formation rate for \ourobj and discuss the link between the observed signatures of a large scale outflow and of the Lyman continuum (LyC) leakage detected in GP galaxies.