\omega}Centauri ({\omega}Cen) is the most enigmatic Galactic globular cluster (GC), with unmatched chemical complexity. We combine photometric and spectroscopic catalogs to identify its distinct stellar populations and to investigate their spatial distribution and chemical properties, uncovering new insights into the cluster's formation history. We identify the iron-poor stars commonly found in GCs: the first population (1P), with halo-like chemical composition, and the second population (2P), enriched in elements produced by p-capture processes. Similarly, we divided the iron-rich stars (the anomalous stars) into two groups: the AI and the AII, exhibiting light-element abundance distributions similar to 1P and 2P stars, respectively. The wide extension of our dataset (five times the half-light radius) allowed us to directly and unambiguously compare the fraction of these populations at different radii. We find that 2P and AII stars are more centrally concentrated than the 1P and AI. The remarkable similarities between the 1P-2P and AI-AII radial distributions strongly suggest that these two groups of stars originated from similar mechanisms. Our chemical analysis indicates that the 1P and AI stars (the lower stream) developed their inhomogeneities through core-collapse supernova (and possibly other massive stars') self-enrichment, and that these populations contributed p-capture-processed material to the intracluster medium, from which the chemically extreme 2P and AII stars (the upper stream) formed. Additional polluters, such as intermediate-mass asymptotic giant branch stars and Type Ia supernovae, likely played a role in shaping the AII. Finally, we propose that 2P and AII stars with intermediate light-element abundances (the middle stream) formed via dilution between the pure ejecta that created the upper stream and lower-stream material.

We present the first measurement of the cross-correlation between anisotropic birefringence and galaxy number counts, utilizing polarization data from Planck NPIPE and the Quaia quasar catalog. By employing a QML/pseudo-$C_\ell$ combined estimator, we compute the angular power spectrum up to $\ell=191$ from birefringence and clustering maps at $N_{\rm side}=64$. Our analysis indicates that the observed spectrum is well consistent with the null-hypothesis, with a probability to exceed of 37% and an estimated scale-invariant amplitude of $A^{\mathcal{D}_\ell}=(2.22\pm2.09)\times10^{-4}\,\text{deg}$, at the 68% confidence level. Finally, we derive constraints on the axion-parameters within an early dark energy model of birefringence. Our findings reveal an unprecedented upper bound on the axion-photon coupling down to $g_{\phi\gamma}=10^{-15}\,\text{GeV}^{-1}$ for masses around $10^{-32}\,\text{eV}$ and high initial misalignment angles. This result opens a previously unexplored window in parameter space, providing the first constraint in this ultra-light mass regime.

JWST has identified a large population of faint, broad-line active galactic nuclei (AGN) in the early universe that are powered by black holes (BHs) that often appear overmassive relative to their host galaxies. In this study, we examine the relationship between BH mass and galaxy stellar mass at $33\sigma$ above the relationship measured for local broad-line AGN. We derive an intrinsic scatter in this relationship of $0.9$ dex, which does not vary over the redshift range of our sample. We also find that the $M_{\rm BH}/M_{\star}$ ratio increases by $2.3$ dex from $z = 3.5$ and $z = 6.5$ with a confidence level of > 3\sigma. We attribute this trend with the increasing fraction of LRDs in our sample at z>4 as their host masses are $\sim1$ dex lower than the non-LRD AGN in our sample. These results support a picture in which the BHs powering JWST's broad-line AGN are genuinely overmassive and become increasingly so with redshift. We discuss the implications of our findings on early BH growth relative to that of their host galaxies and the constraints it places on BH seeding models.

We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This enables the detection of LRDs over a wider redshift range and with less contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at z&lt;8 and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are $\sim1$ dex more numerous than X-ray and UV selected AGN at z~5-7. Within our sample, we have identified the first two X-ray detected LRDs. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec observations from the RUBIES survey of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is 71\% for sources with F444W<26.5. In addition, we find three LRDs with blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low-ionization gas from near the central engine of these faint, red AGN.

The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The i&lt;21 sources have sub-percent photometric redshift accuracy and even the faintest sources at $25

Context. Low- and intermediate-mass giants undergo a complex chemical evolution that has yet to be observationally probed. The influence of core helium flash on the chemical composition of stellar atmospheres has been an open question since its theoretical prediction 60 years ago. Aims. Based on high-resolution spectral observations of 44 open star clusters in the Gaia-ESO survey, our aim is to perform the first large-scale homogeneous investigation into the carbon and nitrogen photospheric content of low- and intermediate-mass giant stars in different phases of evolution. Methods. We determined carbon and nitrogen abundances using spectral synthesis of the C2 Swan (1,0) band head at 5135 Å and C2 Swan (0,1) band head at 5635.5 Å, 12C14N bands in the interval 6470 - 6490 Å, and the forbidden [O i] line at 6300.31 Å. Results. We revealed differences in C/N abundance ratios between pre- and post-core-He-flash stars. The lower C/N ratios in core He-burning red clump stars are mainly due to the enhancement of nitrogen abundances. We presented calibrations of the relationship between [C/N] and stellar age for solar metallicity low- and intermediate-mass giants taking into account different evolutionary stages. Conclusions. The C/N abundance ratios in the investigated first-ascent giant stars are slightly less affected by the first dredge-up than predicted by the theoretical models. The rotation-induced extra mixing is not as efficient as theoretically predicted. The core helium flash may trigger additional alterations in carbon and nitrogen abundances that are not yet theoretically modelled. We found that the evolutionary stage of stars must be taken into account when using [C/N] as an age indicator.

We present a statistical study of spatially resolved chemical enrichment in 18 main-sequence galaxies at $z=4$--6, observed with \jwst/NIRSpec IFU as part of the ALPINE-CRISTAL-\jwst\ survey. Performing an optimized reduction and calibration procedure, including local background subtraction, light-leakage masking, stripe removal, and astrometry refinement, we achieve robust emission-line mapping on kiloparsec scales. Although line-ratio distributions vary across galaxies in our sample, we generally find mild central enhancements in [O\,\textsc{iii}]/H$\beta$, [O\,\textsc{ii}]/[O\,\textsc{iii}], [S\,\textsc{ii}]$_{6732}$/[S\,\textsc{ii}]$_{6718}$, H$\alpha$/H$\beta$, and $L_{\rm H\alpha}/L_{\rm UV}$, consistent with elevated electron density, dust obscuration, and bursty star formation accompanied by reduced metallicity and ionization parameter. These features point to inside-out growth fueled by recent inflows of pristine gas. Nevertheless, the median metallicity gradient is nearly flat over a few kpc scale, $\Delta \log({\rm O/H}) = 0.02 \pm 0.01$ dex kpc$^{-1}$, implying efficient chemical mixing through inflows, outflows, and mergers. From pixel-by-pixel stellar and emission-line characterizations, we further investigate the resolved Fundamental Metallicity Relation (rFMR). Metallicity is described by a fundamental plane with stellar mass and SFR surface densities, but with a stronger dependence on $\Sigma_{\rm SFR}$ than seen in local galaxies. Our results indicate that the regulatory processes linking star formation, gas flows, and metal enrichment were already vigorous $\sim$1 Gyr after the Big Bang, producing the nearly flat metallicity gradient and a stronger coupling between star formation and metallicity than observed in evolved systems in the local universe.

We present the first measurement of the cross-correlation between anisotropic birefringence and galaxy number counts, utilizing polarization data from Planck NPIPE and the Quaia quasar catalog. By employing a QML/pseudo-$C_\ell$ combined estimator, we compute the angular power spectrum up to $\ell=191$ from birefringence and clustering maps at $N_{\rm side}=64$. Our analysis indicates that the observed spectrum is well consistent with the null-hypothesis, with a probability to exceed of 37% and an estimated scale-invariant amplitude of $A^{\mathcal{D}_\ell}=(2.22\pm2.09)\times10^{-4}\,\text{deg}$, at the 68% confidence level. Finally, we derive constraints on the axion-parameters within an early dark energy model of birefringence. Our findings reveal an unprecedented upper bound on the axion-photon coupling down to $g_{\phi\gamma}=10^{-15}\,\text{GeV}^{-1}$ for masses around $10^{-32}\,\text{eV}$ and high initial misalignment angles. This result opens a previously unexplored window in parameter space, providing the first constraint in this ultra-light mass regime.

We investigate the X-shaped radio galaxies (XRGs) with optical double-peaked narrow emission (DPNEL) as potential hosts of dual or binary supermassive black holes (SMBHs). Using a sample of 187 XRGs selected from SDSS and DESI optical spectroscopic surveys, we check the AGN nature of both emission components using the BPT diagnostics of multiple emission lines, namely {[O III]$\lambda\lambda$4959,5007}, H$\alpha$, {[N II]$\lambda\lambda$6548,6584}, {[S II]$\lambda\lambda$6716,6731}, and H$\beta$, and mid-infrared colors. We find that the detection rate of [O III] DPNEL features in XRGs is 30% compared to just 1% in the general galaxy population (mostly radio quiet). The dual AGN fraction in DPNEL galaxies is found to depend strongly on the radio luminosity, increasing from $\sim$25% for radio-undetected to $\sim$58% in the radio-detected sample of general DPNEL galaxies. In contrast, the DPNEL XRGs and FR-II radio galaxies having higher radio power show a $\sim$95% likelihood of hosting a dual AGN. Secondly, the detection of companion galaxies in more than 30% of DPNEL XRGs suggests a vital role of mergers in the XRG formation. We also investigate the parsec-scale radio structure of the nuclei of several XRGs using Very Long Baseline Array (VLBA) maps at 1.4 GHz, 4.3 GHz or 7.6 GHz and find a resolved core for only one of the XRGs. However, the flat spectral indices of the VLBA cores along with the DPNEL components exhibiting AGN characteristics, together with the detection of radio-optical offsets between the VLBA and Gaia position, are strongly indicative of XRGs being likely candidates for hosting dual/binary AGNs.

Jets from accreting black hole (BH) X-ray binaries (XRBs) are powerful outflows that release a large fraction of the accretion energy to the surrounding environment, providing a feedback mechanism that may alter the interstellar medium (ISM) properties. Studying accretion and feedback together enables estimates of matter and energy input/output around accreting BHs. We focus on the extended jet structures of the BH-XRB GRS1758-258. First seen in VLA data, these parsec-scale jets arise from jet-ISM interaction and show a Z-shaped morphology. Using the MeerKAT telescope we observed GRS1758-258 in L-band for a total exposure of 7 hr. Applying a calorimetry-based method developed for AGN and later used for XRBs, we estimated the properties of the jets and of the surrounding ISM. We detect a jet and counter-jet terminating in bow-shocks. Within the northern jet lobe we identify synchrotron and bremsstrahlung emission, while the southern lobe is dominated by thermal emission. We measure ISM densities between 10-40 cm-3 across both jets, slightly lower in the northern region. The estimated ages of the two lobes range from 6-51 kyr. The time-averaged jet power lies between 4.4x10^33 and 3.3x10^36 erg/s, with differences between north and south likely due to different local ISM conditions. Comparing new MeerKAT with archival VLA data, we measured a proper motion of 130 mas/yr in a portion of the northern jet. Jet-ISM interaction structures on both sides of GRS1758-258 reveal different ISM properties. The comparison between these structures and those from other XRBs suggests that the lobes in GRS1758-258 are younger and may result from different jet activity phases. The time-averaged energy transferred to the environment is slightly lower than in other XRBs, consistent with the younger age of the lobes in GRS1758-258 relative to those of other systems.

Euclid is expected to establish new state-of-the-art constraints on extensions beyond the standard LCDM cosmological model by measuring the positions and shapes of billions of galaxies. Specifically, its goal is to shed light on the nature of dark matter and dark energy. Achieving this requires developing and validating advanced statistical tools and theoretical prediction software capable of testing extensions of the LCDM model. In this work, we describe how the Euclid likelihood pipeline, Cosmology Likelihood for Observables in Euclid (CLOE), has been extended to accommodate alternative cosmological models and to refine the theoretical modelling of Euclid primary probes. In particular, we detail modifications made to CLOE to incorporate the magnification bias term into the spectroscopic two-point correlation function of galaxy clustering. Additionally, we explain the adaptations made to CLOE's implementation of Euclid primary photometric probes to account for massive neutrinos and modified gravity extensions. Finally, we present the validation of these CLOE modifications through dedicated forecasts on synthetic Euclid-like data by sampling the full posterior distribution and comparing with the results of previous literature. In conclusion, we have identified in this work several functionalities with regards to beyond-LCDM modelling that could be further improved within CLOE, and outline potential research directions to enhance pipeline efficiency and flexibility through novel inference and machine learning techniques.

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

We present the Blue Jay survey, a Cycle-1 JWST program aimed at studying the stellar and gas content of galaxies at Cosmic Noon. The survey consists of deep spectroscopy for 153 targets observed over two pointings in the COSMOS field using the NIRSpec micro-shutter assembly (MSA). We employ the three medium-resolution gratings G140M, G235M, and G395M, with exposure times of 13 hours, 3.2 hours, and 1.6 hours, respectively. We thus obtain full coverage of the 1-5 micron range, corresponding to the entire rest-frame optical wavelength range. The sample is carefully selected to provide a census of galaxies over the redshift range 1.7 < z < 3.5 above a redshift-dependent minimum stellar mass that ranges from 10^8.7 Msun to 10^9.3 this http URL Blue Jay sample is representative of the entire galaxy population at these redshifts, without strong biases in color, star formation rate, or other properties. The sizes of massive galaxies at these redshifts are comparable to the NIRSpec shutters, which requires custom strategies for designing and reducing the observations. Since the standard A-B nod subtraction leads to flux self-subtraction, we construct a master background from empty shutters and subtract it from each of the science spectra. This, in turn, allows for the use of shorter slitlets consisting of only two shutters per galaxy instead of the usual three, with a substantial increase in the multiplexing of the NIRSpec MSA. We measure multi-band photometry using archival JWST and HST observations in two different ways: in a large elliptical aperture encompassing the entire source and from the exact area in the sky where the NIRSpec 1D spectrum is extracted. This enables self-consistent fits of spectroscopic and photometric data. The Blue Jay dataset, which we publicly release, represents the ideal sample for studying the stellar populations, neutral gas, and ionized gas in Cosmic Noon galaxies.

We present reconstructed convergence maps, \textit{mass maps}, from the Dark Energy Survey (DES) third year (Y3) weak gravitational lensing data set. The mass maps are weighted projections of the density field (primarily dark matter) in the foreground of the observed galaxies. We use four reconstruction methods, each is a \textit{maximum a posteriori} estimate with a different model for the prior probability of the map: Kaiser-Squires, null B-mode prior, Gaussian prior, and a sparsity prior. All methods are implemented on the celestial sphere to accommodate the large sky coverage of the DES Y3 data. We compare the methods using realistic $\Lambda$CDM simulations with mock data that are closely matched to the DES Y3 data. We quantify the performance of the methods at the map level and then apply the reconstruction methods to the DES Y3 data, performing tests for systematic error effects. The maps are compared with optical foreground cosmic-web structures and are used to evaluate the lensing signal from cosmic-void profiles. The recovered dark matter map covers the largest sky fraction of any galaxy weak lensing map to date.

We compare the performance of the flat-sky approximation and Limber approximation for the clustering analysis of the photometric galaxy catalogue of Euclid. We study a 6 bin configuration representing the first data release (DR1) and a 13 bin configuration representative of the third and final data release (DR3). We find that the Limber approximation is sufficiently accurate for the analysis of the wide bins of DR1. Contrarily, the 13 bins of DR3 cannot be modelled accurately with the Limber approximation. Instead, the flat-sky approximation is accurate to below $5\%$ in recovering the angular power spectra of galaxy number counts in both cases and can be used to simplify the computation of the full power spectrum in harmonic space for the data analysis of DR3.

Astronomy is a discipline primarily reliant on visual data. However, alternative data representation techniques are being explored, in particular ''sonification'', namely, the representation of data into sound. While there is increasing interest in the astronomical community in using sonification in research and educational contexts, its full potential is still to be explored. This study measured the performance of astronomers and non-astronomers to detect a transit-like feature in time series data (i.e., light curves), that were represented visually or auditorily, adopting different data-to-sound mappings. We also assessed the bias that participants exhibited in the different conditions. We simulated the data of 160 light curves with different signal-to-noise ratios (SNR). We represented them as visual plots or auditory streams with different sound parameters to represent brightness: pitch, duration, or the redundant duration & pitch. We asked the participants to identify the presence of transit-like features in these four conditions in a session that included an equal number of stimuli with and without transit-like features. With auditory stimuli, participants detected transits with performances above the chance level. However, visual stimuli led to overall better performances compared to auditory stimuli and astronomers outperformed non-astronomers. Visualisations led to a conservative response bias (reluctance to answer ''yes, there is a transit''), whereas sonifications led to more liberal responses (proneness to respond ''yes, there is a transit''). Overall, this study contributes to understanding how different representations (visual or auditory) and sound mappings (pitch, duration, duration & pitch) of time series data affect detection accuracy and biases.

Characterizing dust attenuation is crucial for revealing the intrinsic physical properties of galaxies. We present an analysis of dust attenuation in 18 spectroscopically confirmed star-forming main-sequence galaxies at $z = 4.4-5.7$ observed with JWST/NIRSpec IFU and NIRCam, selected from the ALPINE and CRISTAL ALMA large programs. We fit the emission line fluxes from NIRSpec and the broad-band photometry from NIRCam with Prospector, using both spatially integrated emission and $\sim0.6$ kpc pixel-by-pixel measurements. We derive the stellar-to-nebular dust attenuation ratio ($f=E(B-V)_{\mathrm{star}}/E(B-V)_{\mathrm{neb}}$) from the SED fits and the Balmer decrement with H$\alpha$ and H$\beta$. Although individual galaxies show large scatter, the best-fit value is $f = 0.51^{+0.04}_{-0.03}$, slightly higher than that measured for local starburst galaxies. We find weak correlations of $f$ with galaxy properties, increasing with higher specific star-formation rates, younger stellar ages, and more recent star-formation. For the range of $E(B-V)_{\mathrm{star}} = 0.009-0.15$ mag for in our sample, assuming $f = 1$ (often adopted in high-redshift studies) instead of $f = 0.51$ underestimate line luminosities and ionizing photon production efficiency $\xi_\text{ion}$ by $\sim3-36\%$ and $\sim4-46\%$, respectively. We also find that the total stellar masses estimated from spatially-integrated SED fits with a delayed-$\tau$ star-formation histories are systematically smaller than the sum of pixel-by-pixel SED fits, with a median offset of $\sim 0.26$ dex, likely because the integrated fits are biased toward luminous young stellar populations.

BL Lac objects are an extreme type of active galactic nuclei (AGNs) that belong to the largest population of $\gamma$-ray sources: blazars. This class of AGNs shows a double-bumped spectral energy distribution that is commonly described in terms of a synchrotron self-Compton (SSC) emission process, whereas the low-energy component that dominates their emission between the infrared and the X-ray band is tightly connected to the high-energy component that peaks in the $\gamma$-rays. Two strong connections that link radio and mid-infrared emission of blazars to the emission in the $\gamma$-ray band are well established. They constitute the basis for associating $\gamma$-ray sources with their low-energy counterparts. We searched for a possible link between X-ray and $\gamma$-ray emissions for the subclass of BL Lacs using all archival Swift/XRT observations combined with Fermi data for a selected sample of 351 sources. Analyzing $\sim$2400 ks of Swift/XRT observations that were carried out until December 2018, we discovered that above the $\gamma$-ray flux threshold $F_{\gamma}\approx3\times10^{-12}\,\rm{erg}\,\rm{cm}^{-2}\,\rm{s}^{-1}$, 96\% of all \emph{Fermi} BL Lacs have an X-ray counterpart that is detected with signal-to-noise ratio higher than 3. We did not find any correlation or clear trend between X-ray and $\gamma$-ray fluxes and/or spectral shapes, but we discovered a correlation between the X-ray flux and the mid-infrared color. Finally, we discuss on a possible interpretation of our results in the SSC framework.

A solar active region can significantly disrupt the Sun Earth space environment, often leading to severe space weather events such as solar flares and coronal mass ejections. As a consequence, the automatic classification of active region groups is the crucial starting point for accurately and promptly predicting solar activity. This study presents our results concerned with the application of deep learning techniques to the classification of active region cutouts based on the Mount Wilson classification scheme. Specifically, we have explored the latest advancements in image classification architectures, from Convolutional Neural Networks to Vision Transformers, and reported on their performances for the active region classification task, showing that the crucial point for their effectiveness consists in a robust training process based on the latest advances in the field.

Using deep JWST/NIRSpec spectra from the Blue Jay survey, we perform the first systematic investigation of neutral gas content in massive galaxies at Cosmic Noon based on the Ca II H, K absorption lines. We analyze a sample of 9 galaxies at 1.8 < z < 2.8 with stellar masses > 10.6, for which we detect neutral gas absorption both in Ca II and in Na I. After removing the stellar continuum using the best-fit model obtained with Prospector, we fit the excess absorption due to neutral gas in the Ca II H, K doublet and in the Na I D doublet, together with nearby emission lines produced by ionized gas. We measure covering fractions between 0.2 and 0.9 from the Ca II H and K lines, which are spectrally well resolved in the NIRSpec R ~ 1000 observations, unlike the absorption lines in the Na I D doublet. We measure the velocity shift, velocity dispersion, and column density separately for Ca II and Na I. About half of the galaxies present blueshifted Ca II, indicative of an outflow of neutral gas, consistent with previous results based on Na I. The velocity shift and the column density measured from Ca II are correlated with those measured from Na I, implying that these absorption lines trace gas in similar physical conditions. However, the column densities are not in a 1:1 relation, meaning that the relative amount of Ca II and Na I atoms along the line of sight varies with the gas column density. After discussing possible reasons for this behavior, we derive an empirical relation between the column density of Ca II and the column density of Na I and, in a more indirect way, of neutral hydrogen H I. This calibration offers a new way to estimate the outflow mass and the mass outflow rate for the neutral phase from current and future JWST observations of massive galaxies at Cosmic Noon and beyond