We generalize the analytic solutions presented in Pantoni et al. (2019) by including a simple yet effective description of wind recycling and galactic fountains, with the aim of self-consistently investigating the spatially-averaged time evolution of the gas, stellar, metal, and dust content in disc-dominated late-type galaxies (LTGs). Our analytic solutions, when supplemented with specific prescriptions for parameter setting and with halo accretion rates from $N-$body simulations, can be exploited to reproduce the main statistical relationships followed by local LTGs; these involve, as a function of the stellar mass, the star formation efficiency, the gas mass fraction, the gas/stellar metallicity, the dust mass, the star formation rate, the specific angular momentum, and the overall mass/metal budget. Our analytic solutions allow to easily disentangle the diverse role of the main physical processes ruling galaxy formation in LTGs; in particular, we highlight the crucial relevance of wind recycling and galactic fountains in efficiently refurnishing the gas mass, extending the star-formation timescale, and boosting the metal enrichment in gas and stars. All in all, our analytic solutions constitute a transparent, handy, and fast tool that can provide a basis for improving the (subgrid) physical recipes presently implemented in more sophisticated semi-analytic models and numerical simulations, and can offer a benchmark for interpreting and forecasting current and future spatially-averaged observations of local and higher redshift LTGs.

We present a set of new analytic solutions aimed at self-consistently describing the spatially-averaged time evolution of the gas, stellar, metal, and dust content in an individual starforming galaxy hosted within a dark halo of given mass and formation redshift. Then, as an application, we show that our solutions, when coupled to specific prescriptions for parameter setting (inspired by in-situ galaxy-black hole coevolution scenarios) and merger rates (based on numerical simulations), can be exploited to reproduce the main statistical relationships followed by early-type galaxies and by their high-redshift starforming progenitors. Our analytic solutions allow to easily disentangle the diverse role of the main physical processes regulating galaxy formation, to quickly explore the related parameter space, and to make transparent predictions on spatially-averaged quantities. As such, our analytic solutions may provide a basis for improving the (subgrid) physical recipes presently implemented in theoretical approaches and numerical simulations, and can offer a benchmark for interpreting and forecasting current and future broadband observations of high-redshift starforming galaxies.

[ABRIDGED] We aim to provide a holistic view on the typical size and kinematic evolution of massive early-type galaxies (ETGs), that encompasses their high-$z$ star-forming progenitors, their high-$z$ quiescent counterparts, and their configurations in the local Universe. Our investigation covers the main processes playing a relevant role in the cosmic evolution of ETGs. Specifically, their early fast evolution comprises: biased collapse of the low angular momentum gaseous baryons located in the inner regions of the host dark matter halo; cooling, fragmentation, and infall of the gas down to the radius set by the centrifugal barrier; further rapid compaction via clump/gas migration toward the galaxy center, where strong heavily dust-enshrouded star-formation takes place and most of the stellar mass is accumulated; ejection of substantial gas amount from the inner regions by feedback processes, which causes a dramatic puffing up of the stellar component. In the late slow evolution, passive aging of stellar populations and mass additions by dry merger events occur. We describe these processes relying on prescriptions inspired by basic physical arguments and by numerical simulations, to derive new analytical estimates of the relevant sizes, timescales, and kinematic properties for individual galaxies along their evolution. Then we obtain quantitative results as a function of galaxy mass and redshift, and compare them to recent observational constraints on half-light size $R_e$, on the ratio $v/\sigma$ between rotation velocity and velocity dispersion (for gas and stars) and on the specific angular momentum $j_\star$ of the stellar component; we find good consistency with the available multi-band data in average values and dispersion, both for local ETGs and for their $z\sim 1-2$ star-forming and quiescent progenitors.