Spatially-resolved emission line kinematics are invaluable to investigating fundamental galaxy properties and have become increasingly accessible for galaxies at $z\gtrsim0.5$ through sensitive near-infrared imaging spectroscopy and millimeter interferometry. Kinematic modeling is at the core of the analysis and interpretation of such data sets, which at high-z present challenges due to lower signal-to-noise ratio (S/N) and resolution compared to data of local galaxies. We present and test the 3D fitting functionality of DysmalPy, examining how well it recovers intrinsic disk rotation velocity and velocity dispersion, using a large suite of axisymmetric models, covering a range of galaxy properties and observational parameters typical of $z\sim1$-$3$ star-forming galaxies. We also compare DysmalPy's recovery performance to that of two other commonly used codes, GalPak3D and 3DBarolo, which we use in turn to create additional sets of models to benchmark DysmalPy. Over the ranges of S/N, resolution, mass, and velocity dispersion explored, the rotation velocity is accurately recovered by all tools. The velocity dispersion is recovered well at high S/N, but the impact of methodology differences is more apparent. In particular, template differences for parametric tools and S/N sensitivity for the non-parametric tool can lead to differences up to a factor of 2. Our tests highlight and the importance of deep, high-resolution data and the need for careful consideration of: (1) the choice of priors (parametric approaches), (2) the masking (all approaches) and, more generally, evaluating the suitability of each approach to the specific data at hand. This paper accompanies the public release of DysmalPy.

In the standard model of cosmic structure formation, dark matter haloes form by gravitational instability. The process is hierarchical: smaller systems collapse earlier, and later merge to form larger haloes. The galaxy clusters, hosted by the largest dark matter haloes, are at the top of this hierarchy representing the largest as well as the last structures formed in the universe, while the smaller and first haloes are those Earth-sized dark subhaloes which have been both predicted by theoretical considerations and found in numerical simulations, though it does not exist any observational hints of their existence. The probability that a halo of mass $m$ at redshift $z$ will be part of a larger halo of mass $M$ at the present time can be described in the frame of the extended Press & Schecter theory making use of the progenitor (conditional) mass function. Using the progenitor mass function we calculate analytically, at redshift zero, the distribution of subhaloes in mass, formation epoch and rarity of the peak of the density field at the formation epoch. That is done for a Milky Way-size system, assuming both a spherical and an ellipsoidal collapse model. Our calculation assumes that small progenitors do not lose mass due to dynamical processes after entering the parent halo, and that they do not interact with other subhaloes. For a $\mathrm{\Lambda}$CDM power spectrum we obtain a subhalo mass function $\mathrm{d}n/\mathrm{d}m$ proportional to $m^{- \alpha}$ with a model-independent $\alpha \sim 2$. Assuming the dark matter is a weakly interacting massive particle, the inferred distributions is used to test the feasibility of an indirect detection in the $\gamma$-rays energy band of such a population of subhaloes with a GLAST-like satellite.

A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio $r$ parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., $\approx 150$ GHz) and of few tens of arcmin at the lowest (i.e., $\approx 40$ GHz) and highest (i.e., $\approx 400$ GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are $\approx 5$ times less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on $r$ due to systematics below the limit established by the LiteBIRD collaboration.

In the third APOKASC catalog, we present data for the complete sample of 15,808 evolved stars with APOGEE spectroscopic parameters and Kepler asteroseismology. We used ten independent asteroseismic analysis techniques and anchor our system on fundamental radii derived from Gaia $L$ and spectroscopic $T_{\rm eff}$. We provide evolutionary state, asteroseismic surface gravity, mass, radius, age, and the spectroscopic and asteroseismic measurements used to derive them for 12,418 stars. This includes 10,036 exceptionally precise measurements, with median fractional uncertainties in \nmax, \dnu, mass, radius and age of 0.6\%, 0.6\%, 3.8\%, 1.8\%, and 11.1\% respectively. We provide more limited data for 1,624 additional stars which either have lower quality data or are outside of our primary calibration domain. Using lower red giant branch (RGB) stars, we find a median age for the chemical thick disk of $9.14 \pm 0.05 ({\rm ran}) \pm 0.9 ({\rm sys})$ Gyr with an age dispersion of 1.1 Gyr, consistent with our error model. We calibrate our red clump (RC) mass loss to derive an age consistent with the lower RGB and provide asymptotic GB and RGB ages for luminous stars. We also find a sharp upper age boundary in the chemical thin disk. We find that scaling relations are precise and accurate on the lower RGB and RC, but they become more model dependent for more luminous giants and break down at the tip of the RGB. We recommend the usage of multiple methods, calibration to a fundamental scale, and the usage of stellar models to interpret frequency spacings.

We analyze the evolution of massive (log$_{10}$ [$M_\star/M_\odot$] $>10$) galaxies at $z \sim$ 4--8 selected from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam wide-field slitless spectroscopy. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei (AGN). We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for $z \sim$ 1--4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology.

(abridged) We present a new near-infrared survey covering the 2 deg sq COSMOS field. Combining our survey with Subaru B and z images we construct a deep, wide-field optical-infrared catalogue. At Ks<23 (AB magnitudes) our survey completeness is greater than 90% and 70% for stars and galaxies respectively and contains 143,466 galaxies and 13,254 stars. At z~2 our catalogues contain 3931 quiescent and 25,757 star-forming BzK-selected galaxies representing the largest and most secure sample of these objects to date. Our counts of quiescent galaxies turns over at Ks~22 an effect which we demonstrate cannot be due to sample incompleteness. In our survey both the number of faint and bright quiescent objects exceeds the predictions of a semi-analytic galaxy formation model, indicating potentially the need for further refinements in the amount of merging and AGN feedback at z~2 in these models. We measure the angular correlation function for each sample and find that at small scales the correlation function for passive BzK galaxies exceeds the clustering of dark matter. We use 30-band photometric redshifts to derive the spatial correlation length and the redshift distributions for each object class. At Ks<22 we find r_0^{\gamma/1.8}=7.0 +/-0.5h^{-1} Mpc for the passive BzK candidates and 4.7+/-0.8h^{-1} Mpc for the star-forming BzK galaxies. Our pBzK galaxies have an average photometric redshift of z_p~1.4, in approximate agreement with the limited spectroscopic information currently available. The stacked Ks image will be made publicly available from IRSA.

The study of resonant oscillation modes in low-mass red giant branch stars enables their ages to be inferred with exceptional ($\sim$10%) precision, unlocking the possibility to reconstruct the temporal evolution of the Milky Way at early cosmic times. Ensuring the accuracy of such a precise age scale is a fundamental yet difficult challenge. Since the age of red giant branch stars primarily hinges on their mass, an independent mass determination for an oscillating red giant star provides the means for such assessment. We analyze the old eclipsing binary KIC10001167, which hosts an oscillating red giant branch star and is a member of the thick disk of the Milky Way. Of the known red giants in eclipsing binaries, this is the only member of the thick disk that has asteroseismic signal of high enough quality to test the seismic mass inference at the 2% level. We measure the binary orbit and obtain fundamental stellar parameters through combined analysis of light curve eclipses and radial velocities, and perform a detailed asteroseismic, photospheric, and Galactic kinematic characterization of the red giant and binary system. We show that the dynamically determined mass $0.9337\pm0.0077 \rm\ M_{\odot}$ (0.8%) of this 10 Gyr-old star agrees within 1.4% with the mass inferred from detailed modelling of individual pulsation mode frequencies (1.6%). This is now the only thick disk stellar system, hosting a red giant, where the mass has been determined both asteroseismically with better than 2% precision, and through a model-independent method at 1% precision, and we hereby affirm the potential of asteroseismology to define an accurate age scale for ancient stars to trace the Milky Way assembly history.

We analyse deep images from the VISTA survey of the Magellanic Clouds in the YJKs filters, covering 14 sqrdeg (10 tiles), split into 120 subregions, and comprising the main body and Wing of the Small Magellanic Cloud (SMC). We apply a colour--magnitude diagram reconstruction method that returns their best-fitting star formation rate SFR(t), age-metallicity relation (AMR), distance and mean reddening, together with 68% confidence intervals. The distance data can be approximated by a plane tilted in the East-West direction with a mean inclination of 39 deg, although deviations of up to 3 kpc suggest a distorted and warped disk. After assigning to every observed star a probability of belonging to a given age-metallicity interval, we build high-resolution population maps. These dramatically reveal the flocculent nature of the young star-forming regions and the nearly smooth features traced by older stellar generations. They document the formation of the SMC Wing at ages <0.2 Gyr and the peak of star formation in the SMC Bar at 40 Myr. We clearly detect periods of enhanced star formation at 1.5 Gyr and 5 Gyr. The former is possibly related to a new feature found in the AMR, which suggests ingestion of metal-poor gas at ages slightly larger than 1 Gyr. The latter constitutes a major period of stellar mass formation. We confirm that the SFR(t) was moderately low at even older ages.

This paper provides an update of our previous scaling relations (Genzel et al.2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M_solar)=9.0-11.8, and star formation rates relative to that on the MS, delta_MS=SFR/SFR(MS), from 10^{-1.3} to 10^{2.2}. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t_depl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z)^{-0.6}x(delta_MS)^{-0.44}, and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass mu_gas depends on (1+z)^{2.5}x (delta_MS)^{0.52}x(M*)^{-0.36}, which tracks the evolution of the specific star formation rate. The redshift dependence of mu_gas requires a curvature term, as may the mass-dependences of t_depl and mu_gas. We find no or only weak correlations of t_depl and mu_gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.

We explore the multiplicity of exoplanet host stars with high-resolution images obtained with VLT/SPHERE. Two different samples of systems were observed: one containing low-eccentricity outer planets, and the other containing high-eccentricity outer planets. We find that 10 out of 34 stars in the high-eccentricity systems are members of a binary, while the proportion is 3 out of 27 for circular systems. Eccentric-exoplanet hosts are, therefore, significantly more likely to have a stellar companion than circular-exoplanet hosts. The median magnitude contrast over the 68 data sets is 11.26 and 9.25, in H and K, respectively, at 0.30 arcsec. The derived detection limits reveal that binaries with separations of less than 50au are rarer for exoplanet hosts than for field stars. Our results also imply that the majority of high-eccentricity planets are not embedded in multiple stellar systems (24 out of 34), since our detection limits exclude the presence of a stellar companion. We detect the low-mass stellar companions of HD 7449 and HD 211847, both members of our high-eccentricity sample. HD 7449B was already detected by Rodigas et al (2016) and our independent observation is in agreement with this earlier work. HD 211847's substellar companion, previously detected by the radial velocity method, is actually a low-mass star seen face-on. The role of stellar multiplicity in shaping planetary systems is confirmed by this work, although it does not appear as the only source of dynamical excitation.

In this work, we study the impact of an imperfect knowledge of the instrument bandpasses on the estimate of the tensor-to-scalar ratio $r$ in the context of the next-generation LiteBIRD satellite. We develop a pipeline to include bandpass integration in both the time-ordered data (TOD) and the map-making processing steps. We introduce the systematic effect by having a mismatch between the ``real'', high resolution bandpass $\tau$, entering the TOD, and the estimated one $\tau_s$, used in the map-making. We focus on two aspects: the effect of degrading the $\tau_s$ resolution, and the addition of a Gaussian error $\sigma$ to $\tau_s$. To reduce the computational load of the analysis, the two effects are explored separately, for three representative LiteBIRD channels (40 GHz, 140 GHz and 402 GHz) and for three bandpass shapes. Computing the amount of bias on $r$, $\Delta r$, caused by these effects on a single channel, we find that a resolution $\lesssim 1.5$ GHz and $\sigma \lesssim 0.0089$ do not exceed the LiteBIRD budget allocation per systematic effect, \Delta r &lt; 6.5 \times 10^{-6}. We then check that propagating separately the uncertainties due to a resolution of 1 GHz and a measurement error with $\sigma = 0.0089$ in all LiteBIRD frequency channels, for the most pessimistic bandpass shape of the three considered, still produces a \Delta r &lt; 6.5 \times 10^{-6}. This is done both with the simple deprojection approach and with a blind component separation technique, the Needlet Internal Linear Combination (NILC). Due to the effectiveness of NILC in cleaning the systematic residuals, we have tested that the requirement on $\sigma$ can be relaxed to $\sigma \lesssim 0.05$. (Abridged)

A unique galaxy at z = 2.2, zC406690, has a striking clumpy large-scale ring structure that persists from rest UV to near-infrared, yet has an ordered rotation and lies on the star-formation main sequence. We combine new JWST/NIRCam and ALMA band 4 observations, together with previous VLT/SINFONI integral field spectroscopy and HST imaging to re-examine its nature. The high-resolution H$\alpha$ kinematics are best fitted if the mass is distributed within a ring with total mass $M_{\rm{ring}} = 2 \times 10^{10} M_\odot$ and radius $R_{ring}$ = 4.6 kpc, together with a central undetected mass component (e.g., a "bulge") with a dynamical mass of $M_{bulge} = 8 \times 10^{10} M_\odot$. We also consider a purely flux emitting ring superposed over a faint exponential disk, or a highly "cuspy" dark matter halo, both disfavored against a massive ring model. The low-resolution CO(4-3) line and 142GHz continuum emission imply a total molecular and dust gas masses of $M_{mol,gas} = 7.1 \times 10^{10}M_\odot$and$M_{dust} = 3 \times 10^8 M_\odot$ over the entire galaxy, giving a dust-to-mass ratio of 0.7%. We estimate that roughly half the gas and dust mass lie inside the ring, and that $\sim 10\%$ of the total dust is in a foreground screen that attenuates the stellar light of the bulge in the rest-UV to near-infrared. Sensitive high-resolution ALMA observations will be essential to confirm this scenario and study the gas and dust distribution.

We analyze Ha or CO rotation curves (RCs) extending out to several galaxy effective radii for 100 massive, large, star-forming disk galaxies (SFGs) across the peak of cosmic galaxy star formation (z~0.6-2.5), more than doubling the previous sample presented by Genzel et al. (2020) and Price et al. (2021). The observations were taken with SINFONI and KMOS integral-field spectrographs at ESO-VLT, LUCI at LBT, NOEMA at IRAM, and ALMA. We fit the major axis kinematics with beam-convolved, forward models of turbulent rotating disks with bulges embedded in dark matter (DM) halos, including the effects of pressure support. The fraction of dark to total matter within the disk effective radius ($R_e ~ 5 kpc$), $f_DM (R_e)=V_{DM}^2 (R_e)/V_{circ}^2 (R_e)$, decreases with redshift: At z~1 (z~2) the median DM fraction is $0.38\pm 0.23$ ($0.27\pm 0.18$), and a third (half) of all galaxies are "maximal" disks with f_{DM} (R_e)&lt;0.28. Dark matter fractions correlate inversely with the baryonic surface density, and the low DM fractions require a flattened, or cored, inner DM density distribution. At z~2 there is ~40% less dark matter mass on average within $R_e$ compared to expected values based on cosmological stellar-mass halo-mass relations. The DM deficit is more evident at high star formation rate (SFR) surface densities (\Sigma_{SFR}&gt;2.5 M_{\odot} yr^{-1} kpc^{-2}) and galaxies with massive bulges (M_{bulge}&gt;10^{10} M_{\odot}). A combination of stellar or active galactic nucleus (AGN) feedback, and/or heating due to dynamical friction, either from satellite accretion or clump migration, may drive the DM from cuspy into cored mass distributions. The observations plausibly indicate an efficient build-up of massive bulges and central black holes at z~2 SFGs.

(abridged) We present a new near-infrared survey covering the 2 deg sq COSMOS field. Combining our survey with Subaru B and z images we construct a deep, wide-field optical-infrared catalogue. At Ks<23 (AB magnitudes) our survey completeness is greater than 90% and 70% for stars and galaxies respectively and contains 143,466 galaxies and 13,254 stars. At z~2 our catalogues contain 3931 quiescent and 25,757 star-forming BzK-selected galaxies representing the largest and most secure sample of these objects to date. Our counts of quiescent galaxies turns over at Ks~22 an effect which we demonstrate cannot be due to sample incompleteness. In our survey both the number of faint and bright quiescent objects exceeds the predictions of a semi-analytic galaxy formation model, indicating potentially the need for further refinements in the amount of merging and AGN feedback at z~2 in these models. We measure the angular correlation function for each sample and find that at small scales the correlation function for passive BzK galaxies exceeds the clustering of dark matter. We use 30-band photometric redshifts to derive the spatial correlation length and the redshift distributions for each object class. At Ks<22 we find r_0^{\gamma/1.8}=7.0 +/-0.5h^{-1} Mpc for the passive BzK candidates and 4.7+/-0.8h^{-1} Mpc for the star-forming BzK galaxies. Our pBzK galaxies have an average photometric redshift of z_p~1.4, in approximate agreement with the limited spectroscopic information currently available. The stacked Ks image will be made publicly available from IRSA.

We present the updated version of the code used to compute stellar evolutionary tracks in Padova. It is the result of a thorough revision of the major input physics, together with the inclusion of the pre-main sequence phase, not present in our previous releases of stellar models. Another innovative aspect is the possibility of promptly generating accurate opacity tables fully consistent with any selected initial chemical composition, by coupling the OPAL opacity data at high temperatures to the molecular opacities computed with our AESOPUS code (Marigo & Aringer 2009). In this work we present extended sets of stellar evolutionary models for various initial chemical compositions, while other sets with different metallicities and/or different distributions of heavy elements are being computed. For the present release of models we adopt the solar distribution of heavy elements from the recent revision by Caffau et al. (2011), corresponding to a Sun's metallicity Z=0.0152. From all computed sets of stellar tracks, we also derive isochrones in several photometric systems. The aim is to provide the community with the basic tools to model star clusters and galaxies by means of population synthesis techniques.

Accurate physical parameters of exoplanet systems are essential for further exploration of planetary internal structure, atmospheres, and formation history. We aim to use simultaneous multicolour transit photometry to improve the estimation of transit parameters, to search for transit timing variations (TTVs), and to establish which of our targets should be prioritised for follow-up transmission spectroscopy. We performed time series photometric observations of 12 transits for the hot Jupiters HAT-P-19b, HAT-P-51b, HAT-P-55b, and HAT-P-65b using the simultaneous four-colour camera MuSCAT2 on the Telescopio Carlos Sánchez. We collected 56 additional transit light curves from TESS photometry. To derive transit parameters, we modelled the MuSCAT2 light curves with Gaussian processes to account for correlated noise. To derive physical parameters, we performed EXOFASTv2 global fits to the available transit and radial velocity data sets, together with the Gaia DR3 parallax, isochrones, and spectral energy distributions. To assess the potential for atmospheric characterisation, we compared the multicolour transit depths with a flat line and a clear atmosphere model. We consistently refined the transit and physical parameters. We improved the orbital period and ephemeris estimates, and found no evidence for TTVs or orbital decay. The MuSCAT2 broadband transmission spectra of HAT-P-19b and HAT-P-65b are consistent with previously published low-resolution transmission spectra. We also found that, except for HAT-P-65b, the assumption of a planetary atmosphere can improve the fit to the MuSCAT2 data. In particular, we identified HAT-P-55b as a priority target among these four planets for further atmospheric studies using transmission spectroscopy.

The assembly history of the Milky Way (MW) is a rapidly evolving subject, with numerous small accretion events and at least one major merger proposed in the MW's history. Accreted alongside these dwarf galaxies are globular clusters (GCs), which act as spatially coherent remnants of these past events. Using high precision differential abundance measurements from our recently published study, we investigate the likelihood that the MW clusters NGC 362 and NGC 288 are galactic siblings, accreted as part of the Gaia-Sausage-Enceladus (GSE) merger. To do this, we compare the two GCs at the 0.01 dex level for 20+ elements for the first time. Strong similarities are found, with the two showing chemical similarity on the same order as those seen between the three LMC GCs, NGC 1786, NGC 2210 and NGC 2257. However, when comparing GC abundances directly to GSE stars, marked differences are observed. NGC 362 shows good agreement with GSE stars in the ratio of Eu to Mg and Si, as well as a clear dominance in the r- compared to the s-process, while NGC 288 exhibits only a slight r-process dominance. When fitting the two GC abundances with a GSE-like galactic chemical evolution model, NGC 362 shows agreement with both the model predictions and GSE abundance ratios (considering Si, Ni, Ba and Eu) at the same metallicity. This is not the case for NGC 288. We propose that the two are either not galactic siblings, or GSE was chemically inhomogeneous enough to birth two similar, but not identical clusters with distinct chemistry relative to constituent stars.

The Mid-Infrared Instrument (MIRI) aboard the James Webb Space Telescope (JWST) provides the observatory with a huge advance in mid-infrared imaging and spectroscopy covering the wavelength range of 5 to 28 microns. This paper describes the performance and characteristics of the MIRI imager as understood during observatory commissioning activities, and through its first year of science operations. We discuss the measurements and results of the imager's point spread function, flux calibration, background, distortion and flat fields as well as results pertaining to best observing practices for MIRI imaging, and discuss known imaging artefacts that may be seen during or after data processing. Overall, we show that the MIRI imager has met or exceeded all its pre-flight requirements, and we expect it to make a significant contribution to mid-infrared science for the astronomy community for years to come.

We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey (VIKING). This is the first time that a full optical to near-infrared data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning $450~$deg$^2$, allows us to improve significantly the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and - most importantly - solidify our knowledge of the redshift distributions of the sources. Based on a flat $\Lambda$CDM model we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.737_{-0.036}^{+0.040}$ in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with $S_8$ differing by $2.3\sigma$. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.

In this paper we present a new chemical evolution model for the Galaxy which assumes two main infall episodes for the formation of halo-thick disk and thin disk, respectively. We do not try to take into account explicitly the evolution of the halo but we implicitly assume that the timescale for the formation of the halo was of the same order as the timescale for the formation of the thick disk. The formation of the thin-disk is much longer than that of the thick disk, implying that the infalling gas forming the thin-disk comes not only from the thick disk but mainly from the intergalactic medium. The timescale for the formation of the thin-disk is assumed to be a function of the galactocentric distance, leading to an inside-out picture for the Galaxy building. The model takes into account the most up to date nucleosynthesis prescriptions and adopts a threshold in the star formation process which naturally produces a hiatus in the star formation rate at the end of the thick disk phase, as suggested by recent observations. The model results are compared with an extended set of observational constraints. Among these constraints, the tightest one is the metallicity distribution of the G-dwarf stars for which new data are now available. Our model fits very well these new data. We show that in order to reproduce most of these constraints a timescale $\le 1$ Gyr for the (halo)-thick-disk and of 8 Gyr for the thin-disk formation in the solar vicinity are required. We predict that the radial abundance gradients in the inner regions of the disk (R&lt; R_{\odot}) are steeper than in the outer regions, a result confirmed by recent abundance determinations, and that the inner ones steepen in time during the Galactic lifetime.