We investigate the capability of the J-PAS survey to constrain the primordial power spectrum using a non-parametric Bayesian method. Specifically, we analyze simulated power spectra generated by a local oscillatory primordial feature template motivated by non-standard inflation. The feature is placed within the range of scales where the signal-to-noise ratio is maximized, and we restrict the analysis to $k \in [0.02,0.2] \text{ h} \text{ Mpc}^{-1}$, set by the expected J-PAS coverage and the onset of non-linear effects. Each primordial power spectrum is reconstructed by linearly interpolating $N$ knots in the $\{\log k, \log P_{\mathcal{R}}(k)\}$ plane, which are sampled jointly with the cosmological parameters $\{H_0,\Omega_b h^2, \Omega_c h^2\}$ using PolyChord. To test the primordial features, we apply two statistical tools: the Bayes factor and a hypothesis test that localizes the scales where features are detected. We assess the recovery under different J-PAS specifications, including redshift binning, tracer type, survey area, and filter strategy. Our results show that combining redshift bins and tracers allows the detection of oscillatory features as small as 2\%.

The Universe consists of a variety of objects that formed at different epochs, leading to variations in the formation time which represents the time elapsed from the onset of structure formation until the formation time of a particular object. In this work, we present two approaches to reconstruct and constrain the galaxy formation time $t_f(z)$ using non-parametric reconstruction methods, such as Gaussian Processes (GP) and High-performance Symbolic Regression (SR). Our analysis uses age estimates of 32 old passive galaxies and the Pantheon+ type Ia supernova sample, and considers two different values of the Hubble constant $H_0$ from the SH0ES and Planck Collaborations. When adopting the $\Lambda$CDM model and the GP reconstructions, we find \left=0.72_{-0.16}^{+0.14} Gyr (SH0ES) and \left=1.26_{-0.11}^{+0.10} Gyr (Planck). Without considering a specific cosmological model, we obtain \left=0.71 \pm {0.19} Gyr (SH0ES) and \left = 1.35_{-0.23}^{+0.21} Gyr (Planck). Similar values are obtained from the SR reconstructions, with both methods (GP and SR) indicating the same behavior regarding the time evolution of $t_f(z)$. The results also show significant differences in the formation time from SH0ES and Planck values, highlighting the impact of the $H_0$ tension on the cosmological estimates of $t_f(z)$. In particular, the different approaches used in the analysis agree with each other, demonstrating the robustness and consistency of our results. Overall, this study suggests that galaxies have different evolutionary timescales and that $t_f$ is not constant, with noticeable variations at lower redshifts ($z \lesssim 0.5$).

Eddington ratio is a paramount parameter governing the accretion history and life cycles of Active Galactic Nuclei (AGNs). This short review presents a multi-faceted view of the importance of the Eddington ratio spanning varied AGN studies. We find that the Eddington ratio is crucial for standardizing the Radius-Luminosity (R-L) relation - a necessary step for employing quasars (QSOs) as standardizable cosmological probes to help clarify the standing of the Hubble tension. In this data-driven era, we consolidated disparate aspects by developing novel relations borne out of large datasets, such as the robust, nearly universal anti-correlation between fractional variability and Eddington ratio derived from Zwicky Transient Facility (ZTF) data, which is vital for interpreting forthcoming high-cadence surveys like Rubin Observatory's LSST. Addressing the conundrum where JWST results suggest an overabundance of massive high-redshift black holes, we demonstrate that local AGNs offer clarification: Changing-Look AGNs (CLAGNs), driven by rapid Eddington ratio shifts, cluster in the low-accretion regime, a rate independently confirmed by our integral field spectroscopy and photoionization modeling of a well-known Seyfert 2 galaxy, rich in high-ionization, forbidden, coronal lines. Conversely, for the high-redshift, high-luminosity population where traditional reverberation mapping (RM) is highly impractical, photometric reverberation mapping (PRM) offers a rapid alternative to constrain accretion disk sizes, enabling efficient estimates of black hole masses and Eddington ratios. Finally, we developed tailored semi-empirical spectral energy distributions (SEDs) for extremely high-accretion quasars, successfully validating their characteristic extreme physical conditions.

We investigate a relativistic cosmological model with background rotation, sourced by a non-perfect fluid with anisotropic stress. A modified version of the CLASS Boltzmann code is employed to perform MCMC analyses against Cosmic Microwave Background (CMB) and late-time datasets. The results show that current CMB data constrain the present-day rotation parameter to be negligible. As a consequence, the derived cosmological parameters remain consistent with the standard $\Lambda$CDM values. In contrast, late-time probes such as Type Ia supernovae (SNe) and Baryonic Acoustic Oscillations (BAO) allow for a higher level of rotation and yield an increased Hubble constant. However, this comes at the cost of a higher $\sigma_8$, which remains in tension with DES-Y3 measurement. Combining CMB, SNe and BAO data confirms the preference for non-rotation.

An observational tension on estimates of the Hubble parameter, $H_0$, using early and late Universe information, is being of intense discussion in the literature. Additionally, it is of great importance to measure $H_0$ independently of CMB data and local distance ladder method. In this sense, we analyze 15 measurements of the transversal BAO scale, $\theta_{\rm BAO}$, obtained in a weakly model-dependent approach, in combination with other data sets obtained in a model-independent way, namely, Big Bang Nucleosynthesis (BBN) information, 6 gravitationally lensed quasars with measured time delays by the H0LiCOW team, and measures of cosmic chronometers (CC). We find $H_0 = 74.88_{-2.1}^{+1.9}$ km s${}^{-1}$ Mpc${}^{-1}$ and $H_0 = 72.06_{-1.3}^{+1.2}$ km s${}^{-1}$ Mpc${}^{-1}$ from $\theta_{BAO}$+BBN+H0LiCOW and $\theta_{BAO}$+BBN+CC, respectively, in fully accordance with local measurements. Moreover, we estimate the sound horizon at drag epoch, $r_{\rm d}$, independent of CMB data, and find $r_{\rm d}=144.1_{-5.5}^{+5.3}$ Mpc (from $\theta_{BAO}$+BBN+H0LiCOW) and $r_{\rm d} =150.4_{-3.3}^{+2.7}$ Mpc (from $\theta_{BAO}$+BBN+CC). In a second round of analysis, we test how the presence of a possible spatial curvature, $\Omega_k$, can influence the main results. We compare our constraints on $H_0$ and $r_{\rm d}$ with other reported values. Our results show that it is possible to use a robust compilation of transversal BAO data, $\theta_{BAO}$, jointly with other model-independent measurements, in such a way that the tension on the Hubble parameter can be alleviated.

The Open Cluster Chemical Abundances and Mapping (OCCAM) survey seeks to curate a large, comprehensive, uniform dataset of open clusters and member stars to constrain key Galactic parameters. This eighth entry from the OCCAM survey, based on the newly released SDSS-V/MWM Data Release 19 (DR19), has established a sample of 164 high quality open clusters that are used to constrain the radial and azimuthal gradients of the Milky Way. The DR19 cluster sample [Fe/H] abundances are roughly consistent with measurements from other large-scale spectroscopic surveys. However, the gradients we calculate deviate considerably for some elements. We find an overall linear Galactic radial [Fe/H] gradient of $-0.075 \pm 0.006$ dex kpc$^{-1}$ using the cluster's current Galactocentric Radius ($R_{GC}$) and a gradient of $-0.068 \pm 0.005$ dex kpc$^-1$ with respect to the cluster's guiding center radius. We do not find strong evidence for significant evolution of the differential element gradients ([X/Fe]) investigated here (O, Mg, Si, S, Ca, Ti, Cr, Mn, Fe, Co, Ni, Na, Al, K, Ce, Nd). For the first time using the OCCAM sample we have sufficient numbers of clusters to investigate Galactic azimuthal variations. In this work, we do find evidence of azimuthal variations in the measured radial abundance gradient in the Galactic disk using our open cluster sample.

We present initial results from our JWST NIRSpec program to study the $\alpha$-abundances in the M31 disk. The Milky Way has two chemically-defined disks, the low-$\alpha$ and high-$\alpha$ disks, which are closely related to the thin and thick disks, respectively. The origin of the two populations and the $\alpha$-bimodality between them is not entirely clear, although there are now several models that can reproduce the observed features. To help constrain the models and discern the origin, we have undertaken a study of the chemical abundances of the M31 disk using JWST NIRSpec, in order to determine whether stars in M31's disk also show an $\alpha$-abundance bimodality. Approximately 100 stars were observed in our single NIRSpec field at a projected distance of 18 kpc from the M31 center. The 1-D extracted spectra have an average signal-to-noise ratio of 85 leading to statistical metallicity precision of 0.016 dex, $\alpha$-abundance precision of 0.012 dex, and a radial velocity precision 8 km/s. The initial results indicate that, in contrast to the Milky Way, there is no $\alpha$-bimodality in the M31 disk, and no low-$\alpha$ sequence. The entire stellar population falls along a single chemical sequence very similar to the MW's high-alpha component which had a high star formation rate. While this is somewhat unexpected, the result is not that surprising based on other studies that found the M31 disk has a larger velocity dispersion than the MW and is dominated by a thick component. M31 has had a more active accretion and merger history than the MW which might explain the chemical differences.

Aims. The analysis of variability of astronomical sources is of extraordinary interest, as it allows the study of astrophysical phenomena in real time. This paper presents the Javalambre Variability Survey (J-VAR) which leverages the narrow band filters available at the Javalambre Auxiliary Survey Telescope (JAST80) at the Observatorio Astrofísico de Javalambre (OAJ). Methods. The JAST80 equipped with T80Cam, providing a field of view of 2\,square degrees and a pixel scale of 0.55\,arcsec/pixel has been designed for wide-field studies. The main characteristic is the availability of a variety of narrow band filters strategically located on stellar spectral features (the $J0395$ in correspondence of the Ca H+K doublet, the $J0515$ of the Mg $b$ triplet, the $J0660$ of the H$\alpha$ line, and the $J0861$ of the Ca~triplet). This project combines, for the first time, the wide-field with a variety of narrow band filters for a unique variability survey, observing each field 11~times with a standardised observing sequence. The median limiting magnitude for individual exposures are 19.1 mag in $J0395$ and $J0515$, $19.6$ mag in $J0660$ and $J0861$, $19.8$ mag in $i$, and $20.2$ in $g$ and $r$. The typical FWHM of the $r$-band images is $1.5$ arcsec. Results. This article introduces the first data release of J-VAR including more than 6000 individual asteroids, 10\,detected optical transients (4\,discovered supernovae), and 1.3 million light curves of point-sources. On average, J-VAR delivers an unprecedented $\sim$5000 light curves per square degree of 11\,epochs in 7\,bands, opening research opportunities for theoretical studies and new discoveries alike.

We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analysed in the context of the $\Lambda$CDM model, varying 5 cosmological parameters and 19 nuisance parameters modelling astrophysical and measurement systematics. Our modelling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7\% measurement of the amplitude of fluctuations parameter $S_8\equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.784\pm 0.013$. The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the $\Lambda$CDM model. The addition of the third moment improves the constraints on $S_8$ and $\Omega_{\rm m}$ by $\sim$15\% and $\sim$25\% compared to an analysis that only uses second moments. We compare our results with {\it Planck} constraints from the Cosmic Microwave Background (CMB), finding a $2.2$ \textendash $2.8\sigma$ tension in the full parameter space, depending on the combination of moments considered. The third moment independently is in $2.8\sigma$ tension with {\it Planck}, and thus provides a cross-check on analyses of 2-point correlations.

We use Dark Energy Survey Year 3 (DES Y3) clusters with archival X-ray data from XMM-Newton and Chandra to assess the centering performance of the redMaPPer cluster finder and to measure key richness observable scaling relations. In terms of centering, we find that 10-20% of redMaPPer clusters are miscentered with no significant difference in bins of low versus high richness (20<\lambda<40 and \lambda>40) or redshift ($0.2

Following a quasi model-independent approach we measure the transversal BAO mode at high redshift using the two-point angular correlation function (2PACF). The analyses done here are only possible now with the quasar catalogue from the twelfth data release (DR12Q) from the Sloan Digital Sky Survey, because it is spatially dense enough to allow the measurement of the angular BAO signature with moderate statistical significance and acceptable precision. Our analyses with quasars in the redshift interval z = [2.20,2.25] produce the angular BAO scale theta_BAO = 1.77 +- 0.31 deg with a statistical significance of 2.12 sigma (i.e., 97% confidence level), calculated through a likelihood analysis performed using the theoretical covariance matrix sourced by the analytical power spectra expected in the LCDM concordance model. Additionally, we show that the BAO signal is robust -although with less statistical significance- under diverse bin-size choices and under small displacements of the quasars' angular coordinates. Finally, we also performed cosmological parameter analyses comparing the theta_BAO predictions for wCDM and w(a)CDM models with angular BAO data available in the literature, including the measurement obtained here, jointly with CMB data. The constraints on the parameters Omega_M, w_0 and w_a are in excellent agreement with the LCDM concordance model.

Eddington ratio is a paramount parameter governing the accretion history and life cycles of Active Galactic Nuclei (AGNs). This short review presents a multi-faceted view of the importance of the Eddington ratio spanning varied AGN studies. We find that the Eddington ratio is crucial for standardizing the Radius-Luminosity (R-L) relation - a necessary step for employing quasars (QSOs) as standardizable cosmological probes to help clarify the standing of the Hubble tension. In this data-driven era, we consolidated disparate aspects by developing novel relations borne out of large datasets, such as the robust, nearly universal anti-correlation between fractional variability and Eddington ratio derived from Zwicky Transient Facility (ZTF) data, which is vital for interpreting forthcoming high-cadence surveys like Rubin Observatory's LSST. Addressing the conundrum where JWST results suggest an overabundance of massive high-redshift black holes, we demonstrate that local AGNs offer clarification: Changing-Look AGNs (CLAGNs), driven by rapid Eddington ratio shifts, cluster in the low-accretion regime, a rate independently confirmed by our integral field spectroscopy and photoionization modeling of a well-known Seyfert 2 galaxy, rich in high-ionization, forbidden, coronal lines. Conversely, for the high-redshift, high-luminosity population where traditional reverberation mapping (RM) is highly impractical, photometric reverberation mapping (PRM) offers a rapid alternative to constrain accretion disk sizes, enabling efficient estimates of black hole masses and Eddington ratios. Finally, we developed tailored semi-empirical spectral energy distributions (SEDs) for extremely high-accretion quasars, successfully validating their characteristic extreme physical conditions.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The large-scale structure survey J-PAS is taking data since October 2023. In this work, we present a forecast based on the Fisher matrix method to establish its sensitivity to the sum of the neutrino masses. We adapt the Fisher Galaxy Survey Code (FARO) to account for the neutrino mass under various configurations applied to galaxy clustering measurements. This approach allows us to test the sensitivity of J-PAS to the neutrino mass across different tracers, with and without non-linear corrections, and under varying sky coverage. We perform our forecast for two cosmological models: $\Lambda CDM + \sum m_\nu$ and $w_0w_a CDM + \sum m_\nu$. We combine our J-PAS forecast with Cosmic Microwave Background (CMB) data from the Planck Collaboration and Type Ia supernova (SN) data from Pantheon Plus. Our analysis shows that, for a sky coverage of 8,500 square degrees, J-PAS galaxy clustering data alone will constrain the sum of the neutrino masses to an upper limit at 95\% C.L of \sum m_\nu < 0.32 eV for the $\Lambda CDM + \sum m_\nu$ model, and \sum m_\nu < 0.36 eV for the $w_0w_a CDM + \sum m_\nu$ model. When combined with Planck data, the upper limit improves significantly. For J-PAS+Planck at 95\% C.L, we find \sum m_\nu < 0.061 eV for the $\Lambda CDM + \sum m_\nu$ model, and for J-PAS+Planck+Pantheon Plus, we obtain \sum m_\nu < 0.12 eV for the $w_0w_a CDM + \sum m_\nu$ model. These results demonstrate that J-PAS clustering measurements can play a crucial role in addressing challenges in the neutrino sector, including potential tensions between cosmological and terrestrial measurements of the neutrino mass, as well as in determining the mass ordering.

We report the discovery of eight new ultra-faint dwarf galaxy candidates in the second year of optical imaging data from the Dark Energy Survey (DES). Six of these candidates are detected at high confidence, while two lower-confidence candidates are identified in regions of non-uniform survey coverage. The new stellar systems are found by three independent automated search techniques and are identified as overdensities of stars, consistent with the isochrone and luminosity function of an old and metal-poor simple stellar population. The new systems are faint (Mv > -4.7 mag) and span a range of physical sizes (17 pc < $r_{1/2}$ < 181 pc) and heliocentric distances (25 kpc < D < 214 kpc). All of the new systems have central surface brightnesses consistent with known ultra-faint dwarf galaxies (\mu < 27.5 mag arcsec$^{-2}$). Roughly half of the DES candidates are more distant, less luminous, and/or have lower surface brightnesses than previously known Milky Way satellite galaxies. Most of the candidates are found in the southern part of the DES footprint close to the Magellanic Clouds. We find that the DES data alone exclude (p < 0.001) a spatially isotropic distribution of Milky Way satellites and that the observed distribution can be well, though not uniquely, described by an association between several of the DES satellites and the Magellanic system. Our model predicts that the full sky may hold ~100 ultra-faint galaxies with physical properties comparable to the DES satellites and that 20-30% of these would be spatially associated with the Magellanic Clouds.

Accurate statistical measurement with large imaging surveys has traditionally required throwing away a sizable fraction of the data. This is because most measurements have have relied on selecting nearly complete samples, where variations in the composition of the galaxy population with seeing, depth, or other survey characteristics are small. We introduce a new measurement method that aims to minimize this wastage, allowing precision measurement for any class of stars or galaxies detectable in an imaging survey. We have implemented our proposal in Balrog, a software package which embeds fake objects in real imaging in order to accurately characterize measurement biases. We demonstrate this technique with an angular clustering measurement using Dark Energy Survey (DES) data. We first show that recovery of our injected galaxies depends on a wide variety of survey characteristics in the same way as the real data. We then construct a flux-limited sample of the faintest galaxies in DES, chosen specifically for their sensitivity to depth and seeing variations. Using the synthetic galaxies as randoms in the standard Landy-Szalay correlation function estimator suppresses the effects of variable survey selection by at least two orders of magnitude. With this correction, our measured angular clustering is found to be in excellent agreement with that of a matched sample drawn from much deeper, higher-resolution space-based Cosmological Evolution Survey (COSMOS) imaging; over angular scales of 0.004^{\circ} &lt; \theta &lt; 0.2^{\circ}, we find a best-fit scaling amplitude between the DES and COSMOS measurements of $1.00 \pm 0.09$. We expect this methodology to be broadly useful for extending the statistical reach of measurements in a wide variety of coming imaging surveys.

The Sloan Digital Sky Survey-V (SDSS-V) is pioneering panoptic spectroscopy: it is the first all-sky, multi-epoch, optical-to-infrared spectroscopic survey. SDSS-V is mapping the sky with multi-object spectroscopy (MOS) at telescopes in both hemispheres (the 2.5-m Sloan Foundation Telescope at Apache Point Observatory and the 100-inch du Pont Telescope at Las Campanas Observatory), where 500 zonal robotic fiber positioners feed light from a wide-field focal plane to an optical (R$\sim 2000$, 500 fibers) and a near-infrared (R$\sim 22,000$, 300 fibers) spectrograph. In addition to these MOS capabilities, the survey is pioneering ultra wide-field ($\sim$ 4000~deg$^2$) integral field spectroscopy enabled by a new dedicated facility (LVM-I) at Las Campanas Observatory, where an integral field spectrograph (IFS) with 1801 lenslet-coupled fibers arranged in a 0.5 degree diameter hexagon feeds multiple R$\sim$4000 optical spectrographs that cover 3600-9800 angstroms. SDSS-V's hardware and multi-year survey strategy are designed to decode the chemo-dynamical history of the Milky Way Galaxy and tackle fundamental open issues in stellar physics in its Milky Way Mapper program, trace the growth physics of supermassive black holes in its Black Hole Mapper program, and understand the self-regulation mechanisms and the chemical enrichment of galactic ecosystems at the energy-injection scale in its Local Volume Mapper program. The survey is well-timed to multiply the scientific output from major all-sky space missions. The SDSS-V MOS programs began robotic operations in 2021; IFS observations began in 2023 with the completion of the LVM-I facility. SDSS-V builds upon decades of heritage of SDSS's pioneering advances in data analysis, collaboration spirit, infrastructure, and product deliverables in astronomy.

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large scale structure. This analysis will be carried out on data from the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the tangential shears of background galaxies induced by these foreground lenses. We develop a practical approach to modeling the assumptions and systematic effects affecting small scale lensing, which provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement errors, we study how external priors on different subsets of these parameters affect our growth constraints. Degeneracies within the HOD model, as well as between the HOD and the growth function, are identified as the dominant source of complication, with other systematic effects sub-dominant. The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of the galaxy sample that we employ. We conclude that DES data will provide powerful constraints on the evolution of structure growth in the universe, conservatively/optimistically constraining the growth function to 7.9\%/4.8\% with its first-year data that covered over 1000 square degrees, and to 3.9\%/2.3\% with its full five-year data that will survey 5000 square degrees, including both statistical and systematic uncertainties.

We employ the first two years of data from the near-infrared, high-resolution SDSS-III/APOGEE spectroscopic survey to investigate the distribution of metallicity and alpha-element abundances of stars over a large part of the Milky Way disk. Using a sample of ~10,000 kinematically-unbiased red-clump stars with ~5% distance accuracy as tracers, the [alpha/Fe] vs. [Fe/H] distribution of this sample exhibits a bimodality in [alpha/Fe] at intermediate metallicities, -0.9<[Fe/H]<-0.2, but at higher metallicities ([Fe/H]=+0.2) the two sequences smoothly merge. We investigate the effects of the APOGEE selection function and volume filling fraction and find that these have little qualitative impact on the alpha-element abundance patterns. The described abundance pattern is found throughout the range 5