The chemical composition of diffuse interstellar clouds is not fully established. They host an active chemistry despite their relatively low density and the ubiquitous presence of far-UV radiation. To further explore the chemical composition of diffuse clouds, we performed a spectral scan toward the bright radio source BL Lac in the Q band (from 32 to 50 GHz) using the Yebes 40m telescope. Yebes observations were performed interleaving Frequency Switching and Position Switching integrations toward BL Lac, using a spectral resolution of 38 kHz. The data have been reduced with the CLASS software. We achieved an unprecedented sensitivity on the continuum of 0.02 - 0.07 %, allowing the detection of very faint absorption features. We confirm previous detections of HCS+, C3H, C3H+, CH3CN and HC3N in diffuse clouds and report new detections of CCS, C4H, CH3CHO, H2CCO , HNCO and H2CS along the line of sight to BL Lac, with abundances relative to H2 from a few 10{-11} to a few 10{-10}. We compiled molecular detections toward diffuse clouds to obtain the chemical inventory of a typical diffuse interstellar cloud. The chemical inventory of diffuse interstellar clouds includes complex organic species with up to four heavy atoms. These species are efficiently formed in the diffuse interstellar gas and reach abundances similar to those measured in dense photodissociation regions, pointing to similar gas phase chemical processes.

We present a study of the cold molecular gas kinematics in the inner ~ 4-7 kpc (projected sizes) of three nearby Seyfert galaxies, with AGN luminosities of ~ 10$^{44}$ erg/s, using observations of the CO(2-1) emission line, obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at ~ 0.5-0.8$^{\prime\prime}$ (~ 150-400 pc) spatial resolutions. After modeling the CO profiles with multiple Gaussian components, we detected regions with double-peak profiles that exhibit kinematics distinct from the dominant rotational motion. In NGC 6860, a molecular outflow surrounding the bipolar emission of the [O III] ionized gas is observed extending up to $R_{out}$ ~ 560 pc from the nucleus. There is evidence of molecular inflows along the stellar bar, although an alternative scenario, involving a decoupled rotation in a circumnuclear disk (CND) can also explain the observed kinematics. Mrk 915 shows double-peak CO profiles along one of its spiral arms. Due to its ambiguous disk orientation, part of the CO emission can be interpreted as a molecular gas inflow or an outflow reaching $R_{out}$ ~ 2.8 kpc. MCG -01-24-012 has double-peak profiles associated with a CND, perpendicular to the [O III] bipolar emission. The CO in the CND is rotating while outflowing within $R_{out}$ ~ 3 kpc, with the disturbances possibly being caused by the passage of the ionized gas outflow. Overall, the mass inflow rates are larger than the accretion rate needed to produce the observed luminosities, suggesting that only a fraction of the inflowing gas ends up feeding the central black holes. Although we found signatures of AGN feedback on the cold molecular phase, the mass outflow rates of ~ 0.09-3 M$_\odot$/yr indicate an overall weak impact at these AGN luminosities. Nonetheless, we may be witnessing the start of the depletion and ejection of the molecular gas reservoir that has accumulated over time.

We report here the first detection in the interstellar medium of the cyanomidyl radical (HNCN). Using the Yebes 40m and the IRAM 30m telescopes, we have targeted the doublets of the $N$=2$-$1, 4$-$3, 5$-$4, 6$-$5, and 7$-$6 transitions of HNCN toward the molecular cloud G+0.693-0.027. We have detected three unblended lines of HNCN, these are the $N$=6$-$5 doublet and one line of the $N$=4$-$3 transition. Additionally we present one line of the $N$=5$-$4 transition partially blended with emission from other species. The Local Thermodynamic Equilibrium best fit to the data gives a molecular abundance of (0.91$\pm$0.05)$\times$10$^{-10}$ with respect to H$_2$. The relatively low abundance of this species in G+0.693-0.027, and its high reactivity, suggest that HNCN is possibly produced by gas-phase chemistry. Our work shows that this highly reactive molecule is present in interstellar space, and thus it represents a plausible precursor of larger prebiotic molecules with the NCN backbone such as cyanamide (NH$_2$CN), carbodiimide (HNCNH) and formamidine (NH$_2$CHNH).

We report the first detection in interstellar medium of the 1-cyano propargyl radical, HC$_3$HCN. This species is an isomer of the 3-cyano propargyl radical (CH$_2$C$_3$N), which was recently discovered in TMC-1. The 1-cyano propargyl radical was observed in the cold dark cloud TMC-1 using data from the ongoing QUIJOTE line survey, which is being carried out with the Yebes 40m telescope. A total of seven rotational transitions with multiple hyperfine components were detected in the 31.0-50.4 GHz range. We derived a column density of (2.2$\pm$0.2)$\times$10$^{11}$ cm$^{-2}$ and a rotational temperature of 7$\pm$1\,K. The abundance ratio between HC$_3$HCN and CH$_2$C$_3$N is 1.4. The almost equal abundance of these isomers indicates that the two species may be produced in the same reaction with a similar efficiency, probably in the reaction C + CH$_2$CHCN and perhaps also in the reaction C$_2$ + CH$_3$CN and the dissociative recombination with electrons of CH$_2$C$_3$NH$^+$

(Abridged) We explore the chemistry of the most abundant C, O, S, and N bearing species in molecular clouds, in the context of the IRAM 30 m Large Programme Gas phase Elemental abundances in Molecular Clouds (GEMS). In this work, we aim to assess the limitations introduced in the observational works when a uniform density is assumed along the line of sight for fitting the observations, developing a very simple numerical model of a turbulent box. We perform a MHD simulation in order to reproduce the turbulent steady-state of a turbulent box with properties typical of a molecular filament before collapse. We post-process the results of the MHD simulation with a chemical code to predict molecular abundances, and then post-process this cube with a radiative transfer code to create synthetic emission maps for a series of rotational transitions observed during the GEMS project. From the chemical point of view, we find that turbulence produces variations on the predicted abundances, but they are more or less critical depending on the chosen transition and the chemical age. When compared to real observations, the results from the turbulent simulation provides a better fit than when assuming a uniform gas distribution along the line of sight. In the view of our results, we conclude that taking into account turbulence when fitting observations might significantly improve the agreement with model predictions. This is especially important for sulfur bearing species that are very sensitive to the variations of density produced by turbulence at early times (0.1 Myr). The abundance of CO is also quite sensitive to turbulence when considering the evolution beyond a few 0.1 Myr.

We present a kinematical study of the ionized and molecular gas in the central region (~7-14"~100-200pc) of the nearby radio galaxy Cen A. We used JWST/MIRI MRS 5-28$\mu$m observations taken as part of the MIR Characterization of Nearby Iconic galaxy Centers (MICONIC) of the MIRI EC. The two gas phases present contrasting morphologies and kinematics. The brightest emission from the ionized gas, traced with a range of IP lines ([Fe II] to [Ne VI]), is extended along the direction of the radio jet. We also detected emission from low IP emission lines and H$_2$ transitions in the galaxy disk. Both gas phases present rotational motions but also complex kinematics. The observations reveal several ionized gas kinematical features that are consistent with simulation predictions of a jet-driven bubble and outflow interacting with the galaxy ISM. These include broad components in the nuclear line profiles ($\sigma$~600km/s), high velocities (~ +1000, -1400km/s) confined within the nuclear region, velocities of hundreds of km/s in several directions in the central 2", and enhanced velocity dispersions perpendicular to the radio jet. Moreover, we find evidence of shock excitation in the nuclear region based on MIR line ratios. We compared the ionized gas mass outflow rate with Cen A's AGN luminosity and radio jet power and demonstrate that both mechanisms provide sufficient energy to launch the outflow. The noncircular motions observed in the H$_2$ lines can be reproduced with either a warped rotating disk model or a radial component. The latter might be to related to gas streamers detected in cold molecular gas. There is no clear indication of a fast nuclear H$_2$ outflow, only a weak blueshifted component. This could be due to a relatively low nuclear warm H$_2$ column density and/or the limited geometrical coupling of Cen A's inner radio jet with the circumnuclear disk of the galaxy. (Abridged)

The v=1 and v=2 J=1-0 (43 GHz), and v=1 J=2-1 (86 GHz) SiO masers are intense in AGB stars and have been mapped using VLBI showing ring-like distributions. Those of the v=1, v=2 J=1-0 masers are similar, but the spots are rarely coincident, while the v=1 J=2-1 maser arises from a well separated region farther out. These relative locations can be explained by models tools that include the overlap of two IR lines of SiO and H2O. The v=3 J=1-0 line is not directly affected by any line overlap and its spot structure and position, relative to the other lines, is a good test to the standard pumping models. We present single-dish and simultaneous VLBI observations of the v=1, v=2, and v=3 J=1-0 maser transitions of 28SiO in several AGB stars. The spatial distribution of the SiO maser emission in the v=3 J=1-0 transition from AGB stars is systematically composed of a series of spots that occupy a ring-like structure. The overall ring structure is extremely similar to that found in the other 43 GHz transitions and is very different from the structure of the v=1 J=2-1 maser. The positions of the individual spots of the different 43 GHz lines are, however, very rarely coincident, which in general is separated by about 0.3 AU (between 1 and 5 mas). These results are very difficult to reconcile with standard pumping models, which predict that the masers of rotational transitions within a given vibrational state require very similar excitation conditions, while the transitions of different vibrational states should appear in different positions. However, models including line overlap tend to predict v=1, v=2, v=3 J=1-0 population inversion to occur under very similar conditions, while the requirements for v=1 J=2-1 appear clearly different, and are compatible with the observational results.

Geodetic and astrometric VLBI has entered a new era with the implementation of the VLBI Global Observing System (VGOS). These broadband and dual linear polarization observations aim at an accuracy of station coordinates of 1 mm and a reference frame stability of 0.1 mm/year. Although the extended brightness distribution of many of the radio-loud active galactic nuclei observed during geodetic VLBI sessions is resolved by the interferometer, the established processing chain still treats these objects as point sources. We investigate the impact of source structure on the visibility level and develop tools to remove the structure from the visibility data, right after correlation. Here we present our approach and show results obtained from observational VGOS data.

Post-asymptotic giant branch (post-AGB) binaries are surrounded by massive disks of gas and dust that are similar to protoplanetary disks around young stars. We assembled a catalog of all known Galactic post-AGB binaries with disks. We explore the correlations between the different observables with the aim of learning more about potential disk-binary interactions. We compiled spectral energy distributions of 85 Galactic post-AGB binary systems. We built a color-color diagram to differentiate between the different disk morphologies traced by the infrared excess. We categorized the different disk types and searched for correlations with other observational characteristics of these systems. Between 8 and 12% of our targets are surrounded by transition disks, that is, disks having no or low near-infrared excess. We find a strong link between these transition disks and the depletion of refractory elements on the surface of the post-AGB star. We interpret this link as evidence of a mechanism that stimulates the dust and gas separation within the disk and that also produces the transition disk structure. We propose that such a mechanism is likely to be due to a giant planet carving a hole in the disk, effectively trapping the dust in the outer disk parts. We propose disk evolutionary scenarios, depending on the presence of such a giant planet in the disk. We advocate that giant planets can successfully explain the link between the transition disks and the observed depletion of refractory materials. If the planetary scenario is confirmed, disks around post-AGB binaries could be a unique laboratory for testing planet-disk interactions and their influence on the late evolution of binary stars. The question of whether such planets are first- or second-generation bodies remains to be considered. We argue that these disks are ideal for studying planet formation in unprecedented parameter space.

We present the discovery with the QUIJOTE line survey of the cations HC5N+ and HC7N+ in the direction of TMC-1. Seven lines with half-integer quantum numbers from J=25/2-23/2 to 37/2-35/2 have been assigned to HC5N+ and eight lines from J=55/2-53/2 to 71/2-69/2 to HC7N+. Both species have inverted 2Pi ground electronic states with very good estimates for their B0 and AS0 constants based on optical observations. The lines with the lowest J of HC5N+ exhibit multiple components due to the hyperfine structure introduced by the H and N nuclei. However, these different components collapse for the higher J. No hyperfine structure is found for any of the lines of HC7N+. The derived effective rotational and distortion constants for HC5N+ are Beff = 1336.662+/- 0.001 MHz and Deff = 27.4+/-2.6 Hz, while for HC7N+ they are Beff = 567.85036+/-0.00037 MHz and Deff = 4.01+/-0.19 Hz. From the observed intensities, we derived Trot = 5.5+/-0.5K and N = (9.9+/-1.0)x 1010 cm-2 for HC5N+, while we obtained Trot = 8.5+/-0.5K and N = (2.3+/-0.2)x1010 cm-2 for HC7N+. The HC5N/HC5N+, C5N/HC5N+, C5N-/HC5N+, HC7N/HC7N+, HC5N+/HC7N+, and C7N-/HC7N+ abundance ratios are 670+/-80, 4.8+/-0.8, 1.2+/-0.2, 1000+/-150, 4.2+/-0.5, and 2.2+/-0.2, respectively. We have run chemical modelling calculations to investigate the formation and destruction of these new cations. We find that these species are mainly formed through the reactions of H2 and the cations C5N+ and C7N+, and by the reactions of H+ with HC5N and HC7N, while they are mostly destroyed through a reaction with H2 and a dissociative recombination with electrons. Based on the underestimation of the abundances of HC5N+ and HC7N+ by the chemical model by a factor around 20, we suggest that the rate coefficients currently assumed for the reactions of these cations with H2 could be too high by the same factor, something that will be worth investigating.

Context. The current algorithms used for the calibration and analysis of very long baseline interferometry (VLBI) networks that only use linear polarizers (as is the case of the VLBI Global Observing System, VGOS) do not properly account for instrumental and source-intrinsic polarimetry, which can cause errors in geodetic and astronomical products. Aims. We aim to develop a calibration pipeline for VLBI interferometers that observe in a basis of linear polarization, as is the case of VGOS. The products from this pipeline can be used to obtain valuable full-polarization astronomical information from the observed sources, and they can be used to potentially improve the geodetic results. Methods. We used the algorithm PolConvert to write the correlation products in a basis of circular polarization that is compatible with the standard VLBI calibration procedures. In addition to this, we implemented a wide-band global fringe-fitting algorithm that accounts for dispersive effects (ionospheric delay) and allows us to perform full-polarization imaging of all the observed sources, covering the whole frequency band of VGOS. Results. We present the outcome of our pipeline applied to a global IVS VGOS epoch of observations and show example imaging results in total intensity and polarization. We also discuss issues encountered during the analysis and suggest points of improvement in the VGOS system for an optimum geodetic and astronomical exploitation of this interferometer.

Cosmic rays (CRs) have strong influences on the chemistry of dense molecular clouds (MCs). To study the detailed chemistry induced by CRs, we conducted a Yebes W band line survey towards an unshocked MC (which we named as 3C391:NML) associated with supernova remnant (SNR) 3C391. We detected emission lines of 18 molecular species in total and estimated their column densities with local thermodynamic equilibrium (LTE) and non-LTE analysis. Using the abundance ratio N(HCO+)/N(CO) and an upper limit of N(DCO+)/N(HCO+), we estimated the CR ionization rate of 3C391:NML is $\zeta\gtrsim 2.7\times 10^{-14}\rm \ s^{-1}$ with an analytic method. However, we caution on adopting this value because chemical equilibrium, which is a prerequisite of using the equations, is not necessarily reached in 3C391:NML. We observed lower N(HCO+)/N(HOC+), higher N(HCS+)/N(CS), and higher X($l$-C3H+) by an order of magnitude in 3C391:NML than the typical values in quiescent dense MCs. We found that an enhanced CR ionization rate (of order $\sim 10^{-15}$ or $\sim 10^{-14}\rm \ s^{-1}$) is needed to reproduce the observation with chemical model. This is higher than the values found in typical MCs by 2--3 orders of magnitude.

Utilizing JWST MIRI/MRS IFU observations of the kiloparsec scale central regions, we showcase the diversity of ionized gas distributions and kinematics in six nearby Seyfert galaxies included in the GATOS survey. Specifically, we present spatially resolved flux distribution and velocity field maps of six ionized emission lines covering a large range of ionization potentials ($15.8-97.1$ eV). Based on these maps, we showcase the evidence of ionized gas outflows in the six targets, and find some highly disturbed regions in NGC\,5728, NGC\,5506, and ESO137-G034. We propose AGN-driven radio jets plausibly play an important role in triggering these highly disturbed regions. With the outflow rates estimated based on [Ne~{\footnotesize V}] emission, we find the six targets tend to have ionized outflow rates converged to a narrower range than previous finding. These results have important implication for the outflow properties in AGN of comparable luminosity.

We investigate the capabilities of the mid-infrared instrument (MIRI) of James Webb Space Telescope (JWST) to advance our knowledge of AGN dust using the spectral fitting technique on an AGN collection of 21 nearby (z<0.05) AGN (7 type-1 and 14 type-2) observations obtained with the medium resolution spectroscopy (MRS) mode. This collection includes publicly available AGN and data from the collaboration of Galactic Activity, Torus, and Outflow Survey (GATOS). We developed a tool named MRSPSFisol that decomposes MRS cubes into point-like and extended contributions. We found statistically good fits for 12 targets with current AGN dust models. The model that provides good fits (chi2/dof<2) for {these 12 targets} assumes a combination of clumpy and smooth distribution of dust in a flare-disk geometry where the dust grain size is a free parameter. Still, two and one AGN statistically prefer the disk+wind and the classical clumpy torus model, respectively. However, the currently available models fail to reproduce 40% of the targets, likely due to the extreme silicate features not well reproduced by the models and signatures of water-ice and aliphatic hydrocarbon absorption features in most targets. New models exploring, for instance, new chemistry, are needed to explain the complexity of AGN dust continuum emission observed by JWST.

Molecular gas plays a critical role in regulating star formation and nuclear activity in galaxies. Sulphur bearing molecules, such as H2S, are sensitive to the physical and chemical environments in which they reside and are potential tracers of shocked, dense gas in galactic outflows and active galactic nuclei (AGN). We aim to investigate the origin of H2S emission and its relation to dense gas and outflow activity in the central regions of nearby infrared luminous galaxies. We present ALMA Band 5 observations of the ortho H2S 1(1,0) 1(0,1) transition in three nearby galaxies: NGC 1377, NGC 4418, and NGC 1266. We perform radiative transfer modelling using RADEX to constrain the physical conditions of the H2S emitting gas and compare the results to ancillary CO and continuum data. We detect compact H2S emission in all three galaxies, arising from regions smaller than approximately 150 parsecs. The H2S spectral profiles exhibit broad line wings, suggesting an association with outflowing or shocked gas. In NGC 4418, H2S also appears to be tracing gas that is counterrotating. A peculiar red shifted emission feature may correspond to inflowing gas, or possibly a slanted outflow. RADEX modelling indicates that the H2S emitting gas has high densities (molecular hydrogen density greater than 10^7 cm^-3) and moderately warm temperatures (between 40 and 200 Kelvin). The derived densities exceed those inferred from CO observations, implying that H2S traces denser regions of the interstellar medium.

We aim to study the excitation conditions of the molecular gas in the rotating disk of the Red Rectangle, the only post-Asymptotic-Giant-Branch object in which the existence of an equatorial rotating disk has been demonstrated. For this purpose, we developed a complex numerical code that accurately treats radiative transfer in 2-D, adapted to the study of molecular lines from rotating disks. We present far-infrared Herschel/HIFI observations of the 12CO and 13CO J=6-5, J=10-9, and J=16-15 transitions in the Red Rectangle. We also present our code in detail and discuss the accuracy of its predictions, from comparison with well-tested codes. Theoretical line profiles are compared with the empirical data to deduce the physical conditions in the disk by means of model fitting. We conclude that our code is very efficient and produces reliable results. The comparison of the theoretical predictions with our observations reveals that the temperature of the Red Rectangle disk is typically ~ 100-150 K, about twice as high as previously deduced from mm-wave observations of lower-J lines. We discuss the relevance of these new temperature estimates for understanding the thermodynamics and dynamics of this prototype object, as well as for interpreting observations of other rarely studied post-AGB disks. Despite our sophisticated treatment of the line formation, our model cannot explain the relatively strong line-wing emission for intermediate-J transitions. We argue that a model including a rotating disk only cannot reproduce these data and suggest that there is an additional extended (probably bipolar) structure expanding at about 7--15 km/s.

Type-2 quasars (QSO2s) are active galactic nuclei (AGN) seen through a significant amount of dust and gas that obscures the central supermassive black hole and the broad line region. Here we present new mid-infrared spectra of the central kiloparsec of five optically-selected QSO2s at redshift z~0.1 obtained with JWST/MIRI/MRS. These QSO2s belong to the QSOFEED sample and they have log Lbol=45.5-46.0 erg/s, global SFRs that place them above the main sequence, and practically identical optical spectral shape and [OIII] luminosity, but their nuclear mid-infrared spectra exhibit an unexpected diversity of both continua and features. They show: 1) 9.7 micron silicate features going from emission (strength of S9.7=0.5) to relatively strong absorption (S9.7=-1.0) and 18 and 23 micron silicates either in emission or flat. In addition, two of the QSO2s show absorption bands of CO, H2O, and aliphatic grains, indicating different levels of nuclear obscuration across the sample. 2) [NeV]/[NeII] ratios ranging from 0.1 to 2.1 and [NeIII]/[NeII] from 1.0 to 3.5, indicating different coronal line and ionizing continuum strengths. 3) Warm molecular gas masses of 1-4x10^7 Msun and warm-to-cold gas mass ratios of 1-2%, with molecular gas excitation likely due to jet-induced shocks in J1430+1339, and to UV heating and/or turbulence in J1509+0434. 4) PAH emission features with equivalent widths ranging from <0.002 to 0.075 micron, from which we measure a larger contribution from neutral molecules (PAH 11.3/6.2=1.3-3.4) and SFRs<3-7 Msun/yr. This unprecedented dataset allowed us to start exploring the role of various AGN and galaxy properties including ionizing continuum, obscuration, electron density, and jet-ISM interactions on some of the spectral differences listed above, but larger samples are now required to fully understand the diversity of QSO2s' nuclear mid-infrared spectra.

We present a robust method to estimate the redshift of galaxies using Pan-STARRS1 photometric data. Our method is an adaptation of the one proposed by Beck et al. (2016) for the SDSS Data Release 12. It uses a training set of 2313724 galaxies for which the spectroscopic redshift is obtained from SDSS, and magnitudes and colours are obtained from the Pan-STARRS1 Data Release 2 survey. The photometric redshift of a galaxy is then estimated by means of a local linear regression in a 5-dimensional magnitude and colour space. Our method achieves an average bias of $\overline{\Delta z_{\rm norm}}=-2.01 \times 10^{-4}$, a standard deviation of $\sigma(\Delta z_{\rm norm})=0.0298$, and an outlier rate of $P_o=4.32\%$ when cross-validating on the training set. Even though the relation between each of the Pan-STARRS1 colours and the spectroscopic redshifts is noisier than for SDSS colours, the results obtained by our method are very close to those yielded by SDSS data. The proposed method has the additional advantage of allowing the estimation of photometric redshifts on a larger portion of the sky ($\sim 3/4$ vs $\sim 1/3$). The training set and the code implementing this method are publicly available at this http URL.

Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic environments. As part of the ALCHEMI ALMA large program, the emission of C$_2$H in the central molecular zone of the nearby starburst galaxy NGC 253 was mapped at 1.6 " (28 pc) resolution and characterized to understand its chemical origins. Spectral modelling of the N=1-0 through N=4-3 rotational transitions of C$_2$H was used to derive the C$_2$H column densities towards the dense clouds in NGC 253. Chemical modelling, including PDR, dense cloud, and shock models were then used to investigate the chemical processes and physical conditions that are producing the molecular emission. We find high C$_2$H column densities of $\sim 10^{15} cm^{-3}$ detected towards the dense regions of NGC 253. We further find that these column densities cannot be reproduced by assuming that the emission arises from the PDR regions at the edge of the clouds. Instead, we find that the C$_2$H abundance remains high even in the high visual extinction interior of these clouds and that this is most likely caused by a high cosmic-ray ionization rate.

We report the first detection of the phosphorus monoxide ion (PO$^{+}$) in the interstellar medium. Our unbiased and very sensitive spectral survey towards the G+0.693$-$0.027 molecular cloud covers four different rotational transitions of this molecule, two of which ($J$=1$-$0 and $J$=2$-$1) appear free of contamination from other species. The fit performed, assuming Local Thermodynamic Equilibrium conditions, yields a column density of $N$=(6.0$\pm$0.7)$\times$10$^{11}$ cm$^{-2}$. The resulting molecular abundance with respect to molecular hydrogen is 4.5$\times$10$^{-12}$. The column density of PO$^{+}$ normalised by the cosmic abundance of P is larger than those of NO$^{+}$ and SO$^{+}$, normalised by N and S, by factors of 3.6 and 2.3, respectively. The $N$(PO$^{+}$)/$N$(PO) ratio is 0.12$\pm$0.03, more than one order of magnitude higher than those of $N$(SO$^{+}$)/$N$(SO) and $N$(NO$^{+}$)/$N$(NO). These results indicate that P is more efficiently ionised in the ISM than N and S. We have performed new chemical models that confirm that the PO$^+$ abundance is strongly enhanced in shocked regions with high values of cosmic-ray ionisation rates (10$^{-15}-$10$^{-14}$ s$^{-1}$), as occurs in the G+0.693$-$0.027 molecular cloud. The shocks sputter the interstellar icy grain mantles, releasing into the gas phase most of their P content, mainly in the form of PH$_3$, which is converted into atomic P, and then ionised efficiently by cosmic rays, forming P$^+$. Further reactions with O$_2$ and OH produce PO$^{+}$. The cosmic-ray ionisation of PO might also contribute significantly, which would explain the high $N$(PO$^{+}$)/$N$(PO) observed. The relatively high gas-phase abundance of PO$^{+}$ with respect to other P-bearing species stresses the relevance of this species in the interstellar chemistry of P.