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 mechanisms leading to the formation of disks around young stellar objects (YSOs) and to the launching of the associated jets are crucial to the understanding of the earliest stages of star and planet formation. HH 212 is a privileged laboratory to study a pristine jet-disk system. Therefore we investigate the innermost region (<100 AU) around the HH 212-MM1 protostar through ALMA band\,7 observations of methanol. The 8 GHz bandwidth spectrum towards the peak of the continuum emission of the HH 212 system reveals at least 19 transitions of methanol. Several of these lines (among which several vibrationally excited lines in the v$_{\rm t}=1,2$ states) have upper energies above 500 K. They originate from a compact (<135 AU in diameter), hot ($\sim 295$ K) region elongated along the direction of the SiO jet. We performed a fit in the $uv$ plane of various velocity channels of the strongest high-excitation lines. The blue- and red-shifted velocity centroids are shifted roughly symmetrically on either side of the jet axis, indicating that the line-of-sight velocity beyond 0.7 km s$^{-1}$ from systemic is dominated by rotational motions. The velocity increases moving away from the protostar further indicating that the emission of methanol is not associated with a Keplerian disk or rotating-infalling cavity, and it is more likely associated with outflowing gas. We speculate that CH$_3$OH traces a disk wind gas accelerated at the base. The launching region would be at a radius of a few astronomical units from the YSO.

Recent JWST/NIRCam imaging taken for the ultra-deep UNCOVER program reveals a very red dropout object at $z_{\mathrm{phot}}\simeq7.6$, triply imaged by the galaxy cluster Abell 2744 ($z_{\mathrm{d}}=0.308$). All three images are very compact, i.e. unresolved, with a de-lensed size upper-limit of $r_{e}\lesssim35$ pc. The images have apparent magnitudes of $m_{\mathrm{F444W}}\sim25-26$ AB, and the magnification-corrected absolute UV magnitude of the source is $M_{\mathrm{UV},1450}=-16.81\pm0.09$. From the sum of observed fluxes and from a spectral energy distribution (SED) analysis, we obtain estimates of the bolometric luminosities of the source of $L_{\mathrm{bol}}\gtrsim10^{43} \frac{\mathrm{erg}}{\mathrm{s}}$ and $L_{\mathrm{bol}}\sim10^{44}-10^{46} \frac{\mathrm{erg}}{\mathrm{s}}$, respectively. Based on its compact, point-like appearance, its position in color-color space and the SED analysis, we tentatively conclude that this object is a UV-faint dust-obscured quasar-like object, i.e. an active galactic nucleus (AGN) at high redshift. We also discuss other alternative origins for the object's emission features, including a massive star cluster, Population III, supermassive, or dark stars, or a direct-collapse black hole. Although populations of red galaxies at similar photometric redshifts have been detected with JWST, this object is unique in that its high-redshift nature is corroborated geometrically by lensing, that it is unresolved despite being magnified -- and thus intrinsically even more compact -- and that it occupies notably distinct regions in both size-luminosity and color-color space. Planned UNCOVER JWST/NIRSpec observations, scheduled in Cycle 1, will enable a more detailed analysis of this object.

Studying molecular gas in nearby galaxies using hydrogen cyanide (HCN) as a tracer for higher densities than CO emission still poses a significant challenge. Even though several galaxies have HCN maps on a few kpc scales, higher-resolution maps are still required. Our goal is to examine the contrast in intensity between two tracers that probe different density regimes - HCN(1-0)/CO(2-1) ratio - and their kinematics across NGC 253. By utilizing the advanced capabilities of the Atacama Large Millimeter/submillimeter Array (ALMA), we can map these features at high resolution across a large field of view and uncover the nature of such dense gas in extragalactic systems. We present new ALMA Atacama Compact Array and Total Power (ACA+TP) observations of the HCN emission across NGC 253, covering the inner 8.6' of the galaxy disk at 300 pc scales. We analyze the integrated intensity and mean velocity of HCN and CO along each line of sight and use SCOUSE software to perform spectral decomposition, which considers each velocity component separately. Molecular gas traced by HCN piles up in a ring-like structure at a radius of 2 kpc. The HCN emission is enhanced by 2 orders of magnitude in the central 2 kpc regions, beyond which its intensity decreases with increasing galactocentric distance. The number of components in the HCN spectra shows a robust environmental dependence, with multiple velocity features across the center and bar. We have identified an increase in the HCN/CO ratio in these regions, corresponding to a velocity component likely associated with a molecular outflow. We have also discovered that the ratio between the total infrared luminosity and dense gas mass, which indicates the star formation efficiency of dense gas, is anti-correlated with the molecular gas surface density up to approximately 200 Msul/pc^2. In contrast, beyond this point, the ratio starts to increase.

We compare the properties of giant molecular clouds (GMCs) in M51 identified by the Plateau de Bure Interferometer Whirlpool Arcsecond Survey (PAWS) with GMCs identified in wide-field, high resolution surveys of CO emission in M33 and the Large Magellanic Cloud (LMC). We find that GMCs in M51 are larger, brighter and have higher velocity dispersions relative to their size than equivalent structures in M33 and the LMC. These differences imply that there are genuine variations in the average mass surface density of the different GMC populations. To explain this, we propose that the pressure in the interstellar medium surrounding the GMCs plays a role in regulating their density and velocity dispersion. We find no evidence for a correlation between size and linewidth in any of M51, M33 or the LMC when the CO emission is decomposed into GMCs, although moderately robust correlations are apparent when regions of contiguous CO emission (with no size limitation) are used. Our work demonstrates that observational bias remains an important obstacle to the identification and study of extragalactic GMC populations using CO emission, especially in molecule-rich galactic environments.

We present results from a program of Atacama Large Millimeter Array (ALMA) 1.3 mm (Band 6) molecular line mapping of a sample of nearby, bipolar/pinched-waist, molecule-rich PNe (NGC 6302, Hubble 5, NGC 2440, NGC 6445, NGC 2899, and NGC 2818). Maps of $^{12}$CO(2$-$1) and $^{13}$CO(2$-$1) emission as well as emission lines of HCN, HNC, HCO$^+$, CN, and CS $-$ many of these detected in these PNe for the first time $-$ reveal the molecular mass distributions, compositions, and velocity fields of the equatorial and, in some cases, polar regions of the sample PNe. In each case, the bulk of the molecular gas traces an expanding equatorial torus, with torus expansion velocities ranging from $\sim$15 to $\sim$50 km s$^{-1}$ and molecular masses from $\sim$0.002 to $\sim$0.1 $M_\odot$. The inferred molecular torus dynamical ages, which span the range $\sim$500 yr (Hb 5) to $\sim$11000 yr (NGC 2818), provide support for a model wherein molecular torus ejection precedes bipolar lobe formation. Collectively, these ALMA survey results provide insight into the rapid structural evolution as well as the zones of irradiated molecular gas within bipolar PNe that are descended from relatively massive progenitors, likely residing in interacting binary systems, over $\sim$10 kyr of the post-AGB evolution of such systems.

Stellar population synthesis is a crucial methodology in astrophysics, enabling the interpretation of the integrated light of galaxies and stellar clusters. By combining empirical and/or theoretical libraries of the spectral energy distribution emitted by simple stellar populations (SSPs) with models of the star formation history (SFH) and chemical evolution, population synthesis facilitates the estimation of essential galaxy properties, such as total stellar mass, star formation rate, mass-weighted age and metallicity, etc. The Population Synthesis Toolkit (PST) is a Python library that offers a comprehensive and flexible framework for stellar population synthesis. Its main goal is to compute composite spectra using different galaxy evolution models and SSP libraries with ease and efficiency. It incorporates additional effects, such as cosmic redshift and dust extinction, and it computes several observable quantities derived from the spectra, including broadband photometric fluxes and equivalent widths.

We performed a detailed analysis of the isotopologues with 13C, 34S, 33S, and 36S of the sulphur-bearing molecules CS, CCS, CCCS,HCS+, HCCS+, and H2CS towards TMC-1 using the QUIJOTE1. The observations were obtained with the Yebes radio telescope. Observations with the IRAM 30m of the most abundant isotopologues of these species are also presented and used to estimate volume densities and to constrain the excitation conditions. We report the first detection in space of C13C34S, CC33S, CCC33S, HC33S+, and HCC34S+. C36S is also detected for the first time in a cold object. We also complemented with maps that provide the spatial distribution of most of these species. Using the available collisional rate coefficients for each species, we modeled the observed line intensities using the large velocity gradient method for the radiative transfer. We report the most complete analysis of the column densities of the CnS family and to compare the abundance ratios of all detected isotopologues. Adopting a T_k for TMC-1 of 9K, we found that n(H2)=0.9-1.5X10^4cm-3 can explain the observed decline in intensity with increasing J. We derived the rot. constants for the C13C34S, CC33S, CCC33S, HC33S+, and HCC34S+ isotopologues from new laboratory data and complemented them with the frequencies of the observed lines. We find that all S isotopologues are consistent with solar isotopic abundance ratios. Accurate 12C/13C abundances were derived and, as previously suggested, the 13C isotopologues of CCS and CCCS show strong abundance anomalies depending on the position of the substituted carbon. Nevertheless, the 12C/13C abundance ratio is practically identical to the solar value for CS, HCS+, and H2CS. We also searched for the isotopologues of other S-bearing molecules. The expected intensities for their 34S and 13C isotopologues are too low to be detected with the present sensitivity of the QUIJOTE, however.

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

AGN feedback is invoked as one of the most relevant mechanisms that shape the evolution of galaxies. Our goal is to understand the interplay between AGN feedback and the interstellar medium in M51, a nearby spiral galaxy with a modest AGN and a kpc-scale radio jet expanding through the disc of the galaxy. For that purpose, we combine molecular gas observations in the CO(1-0) and HCN(1-0) lines from the Plateau de Bure interferometer with archival radio, X-ray, and optical data. We show that there is a significant scarcity of CO emission in the ionisation cone, while molecular gas emission tends to accumulate towards the edges of the cone. The distribution and kinematics of CO and HCN line emission reveal AGN feedback effects out to r~500pc, covering the whole extent of the radio jet, with complex kinematics in the molecular gas which displays strong local variations. We propose that this is the result of the almost coplanar jet pushing on molecular gas in different directions as it expands; the effects are more pronounced in HCN than in CO emission, probably as the result of radiative shocks. Following previous interpretation of the redshifted molecular line in the central 5" as caused by a molecular outflow, we estimate the outflow rates to be Mdot_H2~0.9Msun/yr and Mdot_dense~0.6Msun/yr, which are comparable to the molecular inflow rates (~1Msun/yr); gas inflow and AGN feedback could be mutually regulated processes. The agreement with findings in other nearby radio galaxies suggests that this is not an isolated case, and probably the paradigm of AGN feedback through radio jets, at least for galaxies hosting low-luminosity active nuclei.

The enormous radiative and mechanical luminosities of massive stars impact a vast range of scales and processes, from the reionization of the universe, to the evolution of galaxies, to the regulation of the interstellar medium, to the formation of star clusters, and even to the formation of planets around stars in such clusters. Two main classes of massive star formation theory are under active study, Core Accretion and Competitive Accretion. In Core Accretion, the initial conditions are self-gravitating, centrally concentrated cores that condense with a range of masses from the surrounding, fragmenting clump environment. They then undergo relatively ordered collapse via a central disk to form a single star or a small-N multiple. In this case, the pre-stellar core mass function has a similar form to the stellar initial mass function. In Competitive Accretion, the material that forms a massive star is drawn more chaotically from a wider region of the clump without passing through a phase of being in a massive, coherent core. In this case, massive star formation must proceed hand in hand with star cluster formation. If stellar densities become very high near the cluster center, then collisions between stars may also help to form the most massive stars. We review recent theoretical and observational progress towards understanding massive star formation, considering physical and chemical processes, comparisons with low and intermediate-mass stars, and connections to star cluster formation.

We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density of Companions method we find a characteristic spatial scale at ~560 AU, and we use this spatial scale to define the boundary of 19 `cores' in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ~2900 AU, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation towards the southern filament. In addition, the cores of the southern filament are also the densest (within 1100 AU) cores in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 AU seems to be governed by thermal Jeans processes in OMC-1S.

We present the analysis of the multi-phase gas properties in the Seyfert II galaxy NGC 424, using spatially resolved spectroscopic data from JWST/MIRI, part of the Mid-InfraRed Activity of Circumnuclear Line Emission (MIRACLE) program, as well as VLT/MUSE and ALMA. We trace the properties of the multi-phase medium, from cold and warm molecular gas to hot ionised gas, using emission lines such as CO(2-1), H2 S(1), [OIII]5007, [NeIII]15, and [NeV]14. These lines reveal the intricate interplay between the different gas phases within the circumnuclear region, spanning approximately 1.4x1.4 kpc^2, with a resolution of 10 pc. Exploiting the multi-wavelength and multi-scale observations of gas emission we model the galaxy disc rotation curve from scales of a few parsec up to 5 kpc from the nucleus and infer a dynamical mass of 1.09\pm0.08x10^10 M_{\odot} with a disc scale radius of 0.48\pm0.02 kpc. We detect a compact ionised outflow with velocities up to 10^3 km/s, traced by the [OIII], [NeIII], and [NeV] transitions, with no evidence of cold or warm molecular outflows. We suggest that the ionised outflow might be able to inject a significant amount of energy into the circumnuclear region, potentially hindering the formation of a molecular wind, as the molecular gas is observed to be denser and less diffuse. The combined multi-band observations also reveal, in all gas phases, a strong enhancement of the gas velocity dispersion directed along the galaxy minor axis, perpendicular to the high-velocity ionised outflow, and extending up to 1 kpc from the nucleus. Our findings suggest that the outflow might play a key role in such enhancement by injecting energy into the host disc and perturbing the ambient material.

The shapes of only 12 trans-Neptunian objects have been directly measured, offering crucial insights into their internal structure. These properties are strongly connected to the processes that shaped the early Solar System, and provide important clues about its evolution. The aim of the present work is to characterise the size, shape, geometric albedo, and beaming parameter of the TNO (470316) 2007 OC10 . We compared these values to the effective diameter and geometric albedo obtained from thermal data by the TNOs are Cool survey. We also combined occultation and thermal data to constrain the size of a putative unresolved satellite. We predicted an occultation of the star Gaia DR3 2727866328215869952 by 2007 OC10 on 2022 August 22. Four stations detected the occultation. We implemented an elliptical shape model for the projection of 2007 OC10. Following a Bayesian approach, we obtained the posterior probability density in the model parameter space using a Markov chain Monte Carlo method. The elliptical limb of 2007 OC10 has semi-axes of $ 215^{+10}_{-7} \times 141 ^{+24}_{-23}$km, and thus the projected axis ratio is$b/a = 0.58^{+0.16}_{-0.16}$. The area-equivalent diameter is$330^{+56}_{-55}$,km. From our own absolute magnitude value of $H_V = 5.40 \pm 0.02$, the geometric albedo is $p_V = 11.2 ^{+2.1}_{-5.0}$ %. Combining the occultation results with thermal data, we constrain the beaming parameter to $\eta = 1.42^{+0.75}_{-0.58}$. Occultation data reveal that the star is double. The secondary star has a position angle with respect to the primary of $56^{+3}_{-17}$ degrees, has an angular separation of $57^{+4}_{-11}$ mas, and is $1.18^{+0.07}_{-0.07}$ magnitudes fainter than the primary.

During the last years and decades several individual studies and large-scale spectroscopic surveys significantly improved our knowledge of the Galactic metallicity distribution based on open clusters. The availability of Gaia data provided a further step forward in our knowledge. However, still some open issues remain, for example the influence of radial migration on the interpretation of the observed gradients. We used spectroscopic metallicities from individual studies and from the APOGEE survey to compile a sample of 136 open clusters, with a membership verification based on Gaia DR2. Additionally, we present photometric metallicity estimates of 14 open clusters in a somewhat outer Galactic region. Eight age groups allow us to study the evolution of the metallicity gradient in detail, showing within the errors an almost constant gradient of about $-$0.06 dex kpc$^{-1}$. Furthermore, using the derived gradients and an analysis of the individual objects, we estimate a mean migration rate of 1 kpc Gyr$^{-1}$ for objects up to about 2 Gyr. Here, the change of the guiding radius is clearly the main contributor. For older and dynamically hotter objects up to 6 Gyr we infer a lower migration rate of up to 0.5 kpc Gyr$^{-1}$. The influence of epicyclic excursions increases with age and contributes already about 1 kpc to the total migration distance after 6 Gyr. A comparison of our results with available models shows good agreement. However, there is still a lack of a suitable coverage of older objects, future studies are still needed to provide a better sampling in this respect.

In this paper I will describe a new software package developed using the Java programming language, aimed to compute the positions of any Solar System body (among asteroids, comets, planets, and satellites) to help to perform cross-matches of them in observations taken from earth- and space-based observatories. The space telescopes supported are Hubble, James Webb, Euclid, XMM-Newton, Spitzer, Herschel, Gaia, Kepler, Chandra, and TESS, although the flexibility of the software allows to support any other mission without the need to change a single line of code. The orbital elements can be selected among the asteroid database from the Lowell observatory (completed with the cometpro database of comets maintained by the LTE), and the JPL database of minor bodies. The software does not depend on external tools, and performs its own numerical integration of minor bodies. The dynamical model implemented for the Solar System includes the gravity effects of all major bodies, including the Earth, Moon, and Pluto as individual bodies, 16 perturbing asteroids as in other tools, the General Relativity effects, the oblateness of the Sun, Earth, and Moon, and the non-gravitational forces for both comets and asteroids. A complete set of web services allow to compute the cross-matches (that are later to be confirmed, for instance by visual inspection of the images) and also ephemerides of specific bodies. The code is highly optimized and follows the highest standards in terms of software quality and documentation.

The circumstellar environments of PNe are wonderful chemically rich astrophysical laboratories in which the ionization of atoms and the formation of simple and complex molecules can be studied. The new high-sensitivity receivers open the possibility to carry out deep observations, essential to unveil weak atomic and molecular spectra in the mm range. The main goal of this work is to study the emission lines detected in the spectra of the bright C-rich PNe IC 418 and NGC 7027 and to identify all those emission features associated with radio recombination lines (RRLs) of light elements. We aim to analyze the RRLs detected on each source, and to model the sources and derive their physical parameters. This work allows us to provide the most complete and updated catalog of RRLs in space, carried out at 2, 3 and 7 mm with the IRAM 30m and the Yebes 40m radio telescopes. We compare these observational data sets with synthetic models produced with the radiation transfer code Co3RaL. Our observations reveal the presence of several H and He I RRLs at mm wavelengths in the spectra of IC 418 and NGC 7027 and also of HeIIi RRLs in the spectrum of NGC 7027. Many of these lines had remained undetected until now due to their weakness and the lack of high-sensitivity observations at these frequencies. The data also confirm the absence of molecular emission towards IC 418, above a detection level of ~3 mK [Tmb]. These mm observations represent the most extended RRL line survey of two C-rich PNe carried out so far, with most of the lines never reported before. These extremely complete catalogs evidence the importance of high-sensitivity observations and are expected to be very helpful in the line identification process in mm observations, in particular for still unknown or poorly characterized molecular species existing in the vicinity of ionized environments.

We present a multi-phase study of the gas in the circumnuclear region (~1.1x1.0 kpc^2) of the nearby Seyfert 1.8 galaxy NGC 1365, observed in the context of the Mid-IR Activity of Circumnuclear Line Emission (MIRACLE) program. We combined spatially resolved spectroscopic observations from JWST/MIRI, VLT/MUSE, and ALMA to investigate the ionized atomic gas and the warm and cold molecular phases. MIRI data revealed over 40 mid-IR emission lines from ionized and warm molecular gas. Moment maps show that both cold and warm molecular gas follow the rotation of the stellar disk along the circumnuclear ring. The ionized gas displays flux and kinematic patterns that depend on ionization potential (IP): low-IP species (<25 eV) trace the disk, while higher-IP lines (up to ~120 eV) trace outflowing material. The [O III]5700 and [Ne V]14 lines both trace the southeast nuclear outflow cone. Additionally, [Ne V]14 detects the northwest counter-cone, obscured in the optical and thus invisible in [O III]5700. Mid-IR diagnostics, unlike optical ones, clearly reveal the AGN as the primary ionization source in the nucleus. Emission from high-IP species is spatially coincident with the ionization cones and not with star-forming regions. Using the [Ne V]24/[Ne V]14 ratio, we derive an electron density of (750+-440) cm^(-3), in agreement with values from the [S II] optical doublet. For the first time, we apply a fully self-consistent approach combining advanced photoionization and kinematic models (HOMERUN+MOKA3D) to constrain intrinsic outflow properties, overcoming the limitations of simplified classical methods. Exploiting the synergy of JWST/MIRI and VLT/MUSE, HOMERUN reproduces fluxes of over 60 emission lines from optical to mid-IR, disentangling AGN and star formation contributions and yielding robust estimates of outflow mass, geometry, and energetics.

We detected the linear $^3\Sigma^-$ radicals SiC$_3$ and SiC$_5$ toward IRC+10216 using an ultrasensitive line survey gathered with the Yebes 40\,m radio telescope. The derived column densities of $l$-SiC$_3$ and $l$-SiC$_5$ are (3.6$\pm$0.4)$\times$10$^{12}$ cm$^{-2}$ and (1.8$\pm$0.2)$\times$10$^{12}$ cm$^{-2}$, respectively. The linear SiC$_3$ radical is $\sim$2 times less abundant that its singlet rhomboidal prolate isomer, for which we provide a new analysis based on recent sensitive observations in the Q band (7\,mm), and at 3 and 2\,mm with the IRAM 30m telescope. The emission detected from these species arises from the cool external layers of the circumstellar envelope. We speculate whether ion-neutral routes involving SiC$_n$H$_m$$^+$ cations or neutral-neutral reactions involving Si and SiC$_2$ could efficiently synthesize these species.

We present the discovery of the unsubstituted polycyclic aromatic hydrocarbon (PAH) phenalene ($c$-C$_{13}$H$_{10}$) in TMC-1 as part of the QUIJOTE line survey. In spite of the low dipole moment of this three-ring PAH we have found a total of 267 rotational transitions with quantum numbers $J$ and $K_a$ up to 34 and 14, respectively, corresponding to 100 independent frequencies. The identification of this new PAH from our survey was based on the agreement between the rotational parameters derived from the analysis of the lines and those obtained by quantum chemical calculations. Subsequent chemical synthesis of this PAH and the investigation of its laboratory microwave spectrum unequivocally support our identification. The column density of phenalene in TMC-1 is (2.8$\pm$1.6)$\times$10$^{13}$ cm$^{-2}$.