In this paper, we use photometric data from the S-PLUS DR4 survey to identify isolated galaxy pairs and analyse their characteristics and properties. Our results align with previous spectroscopic studies, particularly in luminosity function parameters, suggesting a consistent trait among galaxy systems. Our findings reveal a high fraction of red galaxies across all samples, irrespective of projected distance, velocity difference, or luminosity ratio. We found that the proximity of a neighbour to its central galaxy influences its colour due to environmental effects. We also found that central and neighbour have different behaviours: central galaxies maintain a stable red colour regardless of luminosity, while neighbour colours vary based on luminosity ratios. When the central is significantly brighter, the neighbour tends to be less red. According to our division in red, blue and mixed pairs, we found evidence of galactic conformity. Red pair fractions increase in closer pairs and in pairs of similar luminosity, indicating shared environments promoting red galaxy formation. Analysing local density, the expected colour-density relation is of course recovered, but it is strongly determined by the stellar mass of the pair. In denser environments, the red pair fractions increase, blue pairs decrease and for mixed pairs it depends on their stellar mass: more massive mixed pairs decrease their fraction whereas the lower massive ones increase it. These results shed light on the intricate relationship between galaxy pairs, their characteristics, and environmental influences on colour, providing insights into their evolutionary histories.

The measurement of the structure of stellar populations in the Milky Way disk places fundamental constraints on models of galaxy formation and evolution. Previously, the disk's structure has been studied in terms of populations defined geometrically and/or chemically, but a decomposition based on stellar ages provides a more direct connection to the history of the disk, and stronger constraint on theory. Here, we use positions, abundances and ages for 31,244 red giant branch stars from the SDSS-APOGEE survey, spanning 3 < R_{\mathrm{gc}} < 15 kpc, to dissect the disk into mono-age and mono-[Fe/H] populations at low and high [$\alpha$/Fe]. For each population, with \Delta \mathrm{age} < 2 Gyr and \Delta \mathrm{[Fe/H]} < 0.1 dex, we measure the structure and surface-mass density contribution. We find that low [$\alpha$/Fe] mono-age populations are fit well by a broken exponential, which increases to a peak radius and decreases thereafter. We show that this profile becomes broader with age, interpreted here as a new signal of disk heating and radial migration. High [$\alpha$/Fe] populations are well fit as single exponentials within the radial range considered, with an average scale length of $1.9\pm 0.1$ kpc. We find that the relative contribution of high to low [$\alpha$/Fe] populations at $R_0$ is $f_\Sigma = 18\% \pm 5\%$; high [$\alpha$/Fe] contributes most of the mass at old ages, and low [$\alpha$/Fe] at young ages. The low and high [$\alpha$/Fe] populations overlap in age at intermediate [Fe/H], although both contribute mass at $R_{0}$ across the full range of [Fe/H]. The mass weighted scale height $h_Z$ distribution is a smoothly declining exponential function. High [$\alpha$/Fe] populations are thicker than low [$\alpha$/Fe], and the average $h_Z$ increases steadily with age, between 200 and 600 pc.

We investigate and update observational constraints on cosmological parameters within the $\Lambda$CDM and dynamical dark energy frameworks, using a new compilation of the transverse (or 2D) BAO data, measurements that provide a relatively model-independent estimate of the BAO angular scale at a given redshift. Firstly, we assess the consistency of this compilation with CMB-Planck data and recent BAO results from the DESI collaboration. After confirming minimal tension with CMB data, we perform a series of joint analyses combining CMB data with the 2D~BAO compilation, as well as with several recent Type Ia supernova (SNIa) samples. In all cases, we compare the constraining power of the 2D~BAO data with that of DESI~DR2 samples. Our results indicate that combining 2D~BAO with CMB and SNIa data provides observational constraints that are competitive with those obtained using DESI~DR2. Although the precision of DESI~DR2 results remains higher, as expected due to the more accurate 3D measurements, the 2D~BAO compilation yields strong constraints. For example, in the $\Lambda$CDM context, we find $H_0 = 68.16^{+0.41}_{-0.37} \,\,\, \text{km s}^{-1}\,\text{Mpc}^{-1}$ (CMB + 2D~BAO) and \Sigma m_{\nu} < 0.081~\mathrm{eV} (95\%~CL). These results are comparable to analogous analyses using DESI~DR2. Several other cases are analyzed and presented in the main text. Due to these results, we conclude that this new 2D~BAO compilation is both robust and competitive in constraining cosmological parameters, and, importantly, it does not exhibit significant tension with CMB measurements.

We perform a search for stellar streams around the Milky Way using the first three years of multi-band optical imaging data from the Dark Energy Survey (DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of g > 23.5 with a relative photometric calibration uncertainty of < 1 \%. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a matched-filter in color-magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of eleven new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extra-tidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, as well as the large- and small-scale distribution of dark matter around the Milky Way.

This is the second paper in the HOWLS (higher-order weak lensing statistics) series exploring the usage of non-Gaussian statistics for cosmology inference within \textit{Euclid}. With respect to our first paper, we develop a full tomographic analysis based on realistic photometric redshifts which allows us to derive Fisher forecasts in the ($\sigma_8$, $w_0$) plane for a \textit{Euclid}-like data release 1 (DR1) setup. We find that the 5 higher-order statistics (HOSs) that satisfy the Gaussian likelihood assumption of the Fisher formalism (1-point probability distribution function, $\ell$1-norm, peak counts, Minkowski functionals, and Betti numbers) each outperform the shear 2-point correlation functions by a factor $2.5$ on the $w_0$ forecasts, with only marginal improvement when used in combination with 2-point estimators, suggesting that every HOS is able to retrieve both the non-Gaussian and Gaussian information of the matter density field. The similar performance of the different estimators\inlinecomment{, with a slight preference for Minkowski functionals and 1-point probability distribution function,} is explained by a homogeneous use of multi-scale and tomographic information, optimized to lower computational costs. These results hold for the $3$ mass mapping techniques of the \textit{Euclid} pipeline: aperture mass, Kaiser--Squires, and Kaiser--Squires plus, and are unaffected by the application of realistic star masks. Finally, we explore the use of HOSs with the Bernardeau--Nishimichi--Taruya (BNT) nulling scheme approach, finding promising results towards applying physical scale cuts to HOSs.

We employ a model-independent approach in both the correlation function estimation and the angular BAO feature estimation by computing the angular two-point correlation function. First, we conducted a series of tests to the available DESI tracers to check their representativeness to angular clustering; the result was that, considering the completeness of the first data release across the footprint, we could only make use of the BGS sample for the effective redshifts 0.21 (BGS1) and 0.25 (BGS2). For a reliable analysis in such low redshift, we consider Lagrangian Perturbation Theory at first order on our mocks, which approximately reproduces the expected non-linearities, and generate the corresponding random catalogues. We use a purely statistical method to correct the projection effects and find that our results show reasonable agreement with the $\theta_{\rm BAO}$ expected by the CPL parameters obtained by DESI DR1, being BGS1 $11.78 \pm 1.12$ degrees and BGS2 $11.81 \pm 1.20$ degrees. This means a tension at the $1.5\sigma$ ($2.6\sigma$) level for BGS1 (BGS2) CPL parametrization, while a $2\sigma$ ($3.3\sigma$) discrepancy within the predicted by $\Lambda$CDM. We conclude that, with the current sample available, the use of an angular correlation function serving as the BAO probe, although prefers the CPL parametrization, does not provide conclusive results regarding the best cosmological model.

Fast Radio Bursts (FRBs) are millisecond-duration radio transients with an observed dispersion measure ($DM$) greater than the expected Milky Way contribution, which suggests that such events are of extragalactic origin. Although some models have been proposed to explain the physics of the pulse, the mechanism behind the FRBs emission is still unknown. From FRBs data with known host galaxies, the redshift is directly measured and can be combined with estimates of the $DM$ to constrain the cosmological parameters, such as the baryon number density and the Hubble constant. However, the poor knowledge of the fraction of baryonic mass in the intergalactic medium ($f_{IGM}$) and its degeneracy with the cosmological parameters impose limits on the cosmological application of FRBs. In this work we present a cosmological model-independent method to determine the evolution of $f_{IGM}$ combining the latest FRBs observations with localized host galaxy and current supernovae data. We consider constant and time-dependent $f_{IGM}$ parameterizations and show, through a Bayesian model selection analysis, that a conclusive answer about the time-evolution of $f_{IGM}$ depend strongly on the $DM$ fluctuations due to the spatial variation in cosmic electron density ($\delta$). In particular, our analysis show that the evidence varies from strong (in favor of a growing evolution of $f_{IGM}$ with redshift) to inconclusive, as larger values of $\delta$ are considered.

Precise measurements of black hole masses are essential to understanding the coevolution of these sources and their host galaxies. We develop a novel approach for computing black hole virial masses using measurements of continuum luminosities and emission line widths from partially overlapping, narrow-band observations of quasars; we refer to this technique as single-epoch photometry. This novel method relies on forward-modelling quasar observations for estimating emission line widths, which enables unbiased measurements even for lines coarsely resolved by narrow-band data. We assess the performance of this technique using quasars from the Sloan Digital Sky Survey (SDSS) observed by the miniJPAS survey, a proof-of-concept project of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) collaboration covering $\simeq1\,\mathrm{deg}^2$ of the northern sky using the 56 J-PAS narrow-band filters. We find remarkable agreement between black hole masses from single-epoch SDSS spectra and single-epoch miniJPAS photometry, with no systematic difference between these and a scatter ranging from 0.4 to 0.07 dex for masses from $\log(M_\mathrm{BH})\simeq8$ to 9.75, respectively. Reverberation mapping studies show that single-epoch masses present approximately 0.4 dex precision, letting us conclude that our novel technique delivers black hole masses with only mildly lower precision than single-epoch spectroscopy. The J-PAS survey will soon start observing thousands of square degrees without any source preselection other than the photometric depth in the detection band, and thus single-epoch photometry has the potential to provide details on the physical properties of quasar populations that do not satisfy the preselection criteria of previous spectroscopic surveys.

The Southern Photometric Local Universe Survey (S-PLUS) is a project to map $\sim9300$ sq deg of the sky using twelve bands (seven narrow and five broadbands). Observations are performed with the T80-South telescope, a robotic telescope located at the Cerro Tololo Observatory in Chile. The survey footprint consists of several large contiguous areas, including fields at high and low galactic latitudes, and towards the Magellanic Clouds. S-PLUS uses fixed exposure times to reach point source depths of about $21$ mag in the $griz$ and $20$ mag in the $u$ and the narrow filters. This paper describes the S-PLUS Data Release 4 (DR4), which includes calibrated images and derived catalogues for over 3000 sq deg, covering the aforementioned area. The catalogues provide multi-band photometry performed with the tools \texttt{DoPHOT} and \texttt{SExtractor} -- point spread function (\PSF) and aperture photometry, respectively. In addition to the characterization, we also present the scientific potential of the data. We use statistical tools to present and compare the photometry obtained through different methods. Overall we find good agreement between the different methods, with a slight systematic offset of 0.05\,mag between our \PSF and aperture photometry. We show that the astrometry accuracy is equivalent to that obtained in previous S-PLUS data releases, even in very crowded fields where photometric extraction is challenging. The depths of main survey (MS) photometry for a minimum signal-to-noise ratio $S/N = 3$ reach from $\sim19.5$ for the bluer bands to $\sim21.5$ mag on the red. The range of magnitudes over which accurate \PSF photometry is obtained is shallower, reaching $\sim19$ to $\sim20.5$ mag depending on the filter. Based on these photometric data, we provide star-galaxy-quasar classification and photometric redshift for millions of objects.

We describe the photometric data set assembled from the full six years of observations by the Dark Energy Survey (DES) in support of static-sky cosmology analyses. DES Y6 Gold is a curated data set derived from DES Data Release 2 (DR2) that incorporates improved measurement, photometric calibration, object classification and value added information. Y6 Gold comprises nearly $5000~{\rm deg}^2$ of $grizY$ imaging in the south Galactic cap and includes 669 million objects with a depth of $i_{AB} \sim 23.4$ mag at S/N $\sim 10$ for extended objects and a top-of-the-atmosphere photometric uniformity < 2~{\rm mmag}. Y6 Gold augments DES DR2 with simultaneous fits to multi-epoch photometry for more robust galaxy shapes, colors, and photometric redshift estimates. Y6 Gold features improved morphological star-galaxy classification with efficiency $98.6\%$ and contamination $0.8\%$ for galaxies with 17.5 < i_{AB} < 22.5. Additionally, it includes per-object quality information, and accompanying maps of the footprint coverage, masked regions, imaging depth, survey conditions, and astrophysical foregrounds that are used for cosmology analyses. After quality selections, benchmark samples contain 448 million galaxies and 120 million stars. This paper will be complemented by online data access and documentation.

(119951) 2002 KX14 is a large classical TNO with limited previous observations and unresolved questions regarding its physical properties. Five stellar occultations by 2002 KX14 were observed from 2020 to 2023, involving multiple telescopes across different locations in Europe and the Americas. The five occultations resulted in 15 positive chords, accurately measuring the 2002 KX14's shape and size. The projected ellipse has semi-major and semi-minor axes of $241.0 \pm 7.2$ km and $157.1 \pm 5.2$ km, respectively, corresponding to an average area-equivalent diameter of $389.2 \pm 8.7$ km. The geometric albedo was estimated at $11.9 \pm 0.7\%$.

Studying the SDSS-DR16 quasar catalog, we detect a baryon acoustic oscillation (BAO) signal in the two-point angular correlation function with a statistical significance of $3\sigma$, at an effective redshift of $z_{\rm eff}=1.725$. Using a simple parameterization-comprising a polynomial plus a Gaussian function-we measure the transverse BAO scale as $\theta_{\rm BAO}=1.928^{\circ}\pm0.094^{\circ}$. This measurement is obtained from a narrow redshift shell, $z \in [1.72, 1.73]$ (i.e., $\Delta z=0.01$), thin enough that projection-effect corrections are negligible, making it only weakly dependent on the assumed fiducial cosmology. The only assumption adopted is isotropy in the computation of the correlation function, further ensuring that the result depends only weakly on specific cosmological-model hypotheses. We also investigate possible systematics that could affect the detection or significance of the BAO signal and find them to be subdominant or implausible. When combined with other transverse BAO measurements from the literature, our result shows good concordance-within the $1\sigma$ confidence level-with the cosmological parameter values reported by the Planck and DESI collaborations. This new measurement of the transverse BAO scale, obtained from the SDSS quasar sample with minimal cosmological-model assumptions, provides an additional independent constraint for updated statistical studies aimed at probing the nature of dark energy.

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.

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the Universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (3$\times$2pt) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining $\Lambda$ cold dark matter ($\Lambda$CDM) parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure $\Omega_\mathrm{m}=0.300\pm0.017$ and $\sigma_8=0.797\pm0.026$. Compared to constraints from Planck primary cosmic microwave background (CMB) anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 ($1.2\sigma$) for the two-parameter difference. We further obtain $S_8\equiv\sigma_8(\Omega_\mathrm{m}/0.3)^{0.5}=0.796\pm0.013$ which is lower than the Planck measurement at the $1.6\sigma$ level. The combined SPT cluster, DES 3$\times$2pt, and Planck datasets mildly prefer a nonzero positive neutrino mass, with a 95% upper limit \sum m_\nu<0.25~\mathrm{eV} on the sum of neutrino masses. Assuming a $w$CDM model, we constrain the dark energy equation of state parameter $w=-1.15^{+0.23}_{-0.17}$ and when combining with Planck primary CMB anisotropies, we recover $w=-1.20^{+0.15}_{-0.09}$, a $1.7\sigma$ difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology.

J-PAS (Javalambre Physics of the Accelerating Universe Astrophysical Survey) will present a groundbreaking photometric survey covering 8500 deg$^2$ of the visible sky from Javalambre, capturing data in 56 narrow band filters. This survey promises to revolutionize galaxy evolution studies by observing $\sim$10$^8$ galaxies with low spectral resolution. A crucial aspect of this analysis involves predicting stellar population parameters from the observed galaxy photometry. In this study, we combine the exquisite J-PAS photometry with state-of-the-art single stellar population (SSP) libraries to accurately predict stellar age, metallicity, and dust attenuation with a neural network (NN) model. The NN is trained on synthetic J-PAS photometry from different SSP librares (E-MILES, Charlot & Bruzual, XSL), to enhance the robustness of our predictions against individual SSP model variations and limitations. To create mock samples with varying observed magnitudes we add artificial noise in the form of random Gaussian variations within typical observational uncertainties in each band. Our results indicate that the NN can accurately estimate stellar parameters for SSP models without evident degeneracies, surpassing a bayesian SED-fitting method on the same test set. We obtain median bias, scatter and percentage of outliers $\mu$ = (0.01 dex, 0.00 dex, 0.00 mag), $\sigma_{NMAD}$ = (0.23 dex, 0.29 dex, 0.04 mag), f$_{o}$ = (17 %, 24 %, 1 %) at $i \sim$17 mag for age, metallicity and dust attenuation, respectively. The accuracy of the predictions is highly dependent on the signal-to-noise (S/N) ratio of the photometry, achieving robust predictions up to $i$ $\sim$ 20 mag.

We present the point-spread function (PSF) modeling for weak lensing shear measurement using the full six years of the Dark Energy Survey (DES Y6) data. We review the PSF estimation procedure using the PIFF (PSFs In the Full FOV) software package and describe the key improvements made to PIFF and modeling diagnostics since the DES year three (Y3) analysis: (i) use of external Gaia and infrared photometry catalogs to ensure higher purity of the stellar sample used for model fitting, (ii) addition of color-dependent PSF modeling, the first for any weak lensing analysis, and (iii) inclusion of model diagnostics inspecting fourth-order moments, which can bias weak lensing measurements to a similar degree as second-order modeling errors. Through a comprehensive set of diagnostic tests, we demonstrate the improved accuracy of the Y6 models evident in significantly smaller systematic errors than those of the Y3 analysis, in which all $g$ band data were excluded due to insufficiently accurate PSF models. For the Y6 weak lensing analysis, we include $g$ band photometry data in addition to the $riz$ bands, providing a fourth band for photometric redshift estimation. Looking forward to the next generation of wide-field surveys, we describe several ongoing improvements to PIFF, which will be the default PSF modeling software for weak lensing analyses for the Vera C. Rubin Observatory's Legacy Survey of Space and Time.

According to the perturbed Friedmann model, the difference between Hubble constant measurements in two rest frames, at leading order in velocity, is determined solely by the relative motion of the observers and remains unaffected by the peculiar velocities of the sources. This implies that, when averaging over a sufficiently large and distant set of sources where local nonlinear inhomogeneities are diminished, such a difference should vanish, so that the Hubble flow is statistically uniform, as predicted by the Cosmological Principle -- a core assumption of the standard cosmological paradigm. In previous works, distance measurement compilations, e.g. the CosmicFlows-3 catalogue, were used for this purpose, as it comprises a large number ($\sim 10^4$) of sources of different types. Due to the increasing amount of precise luminosity distance measurements of Type Ia Supenovae (SNe) in the last few years, in this work we investigate whether we can confirm the uniformity of the Hubble flow with low-$z$ SN distances only. By means of the Pantheon+ and SH0ES compilation, we find that the results align well with previous works based on the CF3 catalogue, and are in good agreement with the expected Hubble variance in the standard model across cosmic scales of $20-150$ Mpc. Notably, the Hubble constant difference $\Delta H_0 \approx 0$ is observed at around $85$ Mpc. Despite the smaller sample size ($\sim 10^2$ versus $\sim 10^4$) relative to CF3 at those scales, our analysis show that the Pantheon+ and SH0ES dataset supports the standard model paradigm, which indicates that the Hubble flow becomes statistically uniform at around $70-100$ Mpc, which is compatible with independent determinations of the homogeneity scale based on galaxy number counts.

We introduce redMaGiC, an automated algorithm for selecting Luminous Red Galaxies (LRGs). The algorithm was specifically developed to minimize photometric redshift uncertainties in photometric large-scale structure studies. redMaGiC achieves this by self-training the color-cuts necessary to produce a luminosity-thresholded LRG sample of constant comoving density. We demonstrate that redMaGiC photozs are very nearly as accurate as the best machine-learning based methods, yet they require minimal spectroscopic training, do not suffer from extrapolation biases, and are very nearly Gaussian. We apply our algorithm to Dark Energy Survey (DES) Science Verification (SV) data to produce a redMaGiC catalog sampling the redshift range $z\in[0.2,0.8]$. Our fiducial sample has a comoving space density of $10^{-3}\ (h^{-1} Mpc)^{-3}$, and a median photoz bias ($z_{spec}-z_{photo}$) and scatter $(\sigma_z/(1+z))$ of 0.005 and 0.017 respectively. The corresponding $5\sigma$ outlier fraction is 1.4%. We also test our algorithm with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) and Stripe 82 data, and discuss how spectroscopic training can be used to control photoz biases at the 0.1% level.

We measure the transverse baryon acoustic oscillations (BAO) signal in the local Universe using a sample of blue galaxies from the Sloan Digital Sky Survey (SDSS) survey as a cosmological tracer. The method is weakly dependent on a cosmological model and is suitable for 2D analyses in thin redshift bins to investigate the SDSS data in the interval $z {\in} [0.105, 0.115]$. We detect the transverse BAO signal ${\theta}_{BAO} = 19.8^° {\pm} 1.05^°$ at $z_{eff} = 0.11$, with a statistical significance of $2.2 {\sigma}$. Additionally, we perform tests that confirm the robustness of this angular BAO signature. Supported by a large set of log-normal simulations, our error analyses include statistical and systematic contributions. In addition, considering the sound horizon scale calculated by the Planck Collaboration, $r_{s}^{Planck}$, and the ${\theta}_{BAO}$ value obtained here, we obtain a measurement of the angular diameter distance $D_{A}(0.11) = 258.31 {\pm} 13.71 \,Mpc/h$. Moreover, combining this ${\theta}_{BAO}$ measurement at low redshift with other BAO angular scale data reported in the literature, we perform statistical analyses for the cosmological parameters of some Lambda cold dark matter (${\Lambda}$CDM) type models.

We present a suite of 18 synthetic sky catalogs designed to support science analysis of galaxies in the Dark Energy Survey Year 1 (DES Y1) data. For each catalog, we use a computationally efficient empirical approach, ADDGALS, to embed galaxies within light-cone outputs of three dark matter simulations that resolve halos with masses above ~5x10^12 h^-1 m_sun at z <= 0.32 and 10^13 h^-1 m_sun at z~2. The embedding method is tuned to match the observed evolution of galaxy counts at different luminosities as well as the spatial clustering of the galaxy population. Galaxies are lensed by matter along the line of sight --- including magnification, shear, and multiple images --- using CALCLENS, an algorithm that calculates shear with 0.42 arcmin resolution at galaxy positions in the full catalog. The catalogs presented here, each with the same LCDM cosmology (denoted Buzzard), contain on average 820 million galaxies over an area of 1120 square degrees with positions, magnitudes, shapes, photometric errors, and photometric redshift estimates. We show that the weak-lensing shear catalog, redMaGiC galaxy catalogs and redMaPPer cluster catalogs provide plausible realizations of the same catalogs in the DES Y1 data by comparing their magnitude, color and redshift distributions, angular clustering, and mass-observable relations, making them useful for testing analyses that use these samples. We make public the galaxy samples appropriate for the DES Y1 data, as well as the data vectors used for cosmology analyses on these simulations.