A key uncertainty in galaxy evolution is the physics regulating star formation, ranging from small-scale processes related to the life-cycle of molecular clouds within galaxies to large-scale processes such as gas accretion onto galaxies. We study the imprint of such processes on the time-variability of star formation with an analytical approach tracking the gas mass of galaxies ("regulator model"). Specifically, we quantify the strength of the fluctuation in the star-formation rate (SFR) on different timescales, i.e. the power spectral density (PSD) of the star-formation history, and connect it to gas inflow and the life-cycle of molecular clouds. We show that in the general case the PSD of the SFR has three breaks, corresponding to the correlation time of the inflow rate, the equilibrium timescale of the gas reservoir of the galaxy, and the average lifetime of individual molecular clouds. On long and intermediate timescales (relative to the dynamical timescale of the galaxy), the PSD is typically set by the variability of the inflow rate and the interplay between outflows and gas depletion. On short timescales, the PSD shows an additional component related to the life-cycle of molecular clouds, which can be described by a damped random walk with a power-law slope of $\beta\approx2$ at high frequencies with a break near the average cloud lifetime. We discuss star-formation "burstiness" in a wide range of galaxy regimes, study the evolution of galaxies about the main sequence ridgeline, and explore the applicability of our method for understanding the star-formation process on cloud-scale from galaxy-integrated measurements.

We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($\lambda_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $\lambda_{\rm EDD}$ distribution, and with a positive shift of the break $\lambda^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($\lambda_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $\lambda_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.

To gain insights into long-term Active Galactic Nuclei (AGN) variability, we analyze an AGN sample from the Sloan Digital Sky Survey (SDSS) and compare their photometry with observations from the Hyper Suprime-Cam survey (HSC) observed $\langle 14.85 \rangle$ years after SDSS. On average, the AGN are fainter in HSC than SDSS. We demonstrate that the difference is not due to subtle differences in the SDSS versus HSC filters or photometry. The decrease in mean brightness is redshift dependent, consistent with expectations for a change that is a function of the rest-frame time separation between observations. At a given redshift, the mean decrease in brightness is stronger for more luminous AGN and for objects with longer time separation between measurements. We demonstrate that the dependence on redshift and luminosity of measured mean brightness decrease is consistent with simple models of Eddington ratio variability in AGN on long (Myr, Gyr) timescales. We show how our results can be used to constrain the variability and demographic properties of AGN populations.

We present a numerical framework for the variability of active galactic nuclei (AGN), which links the variability of AGN over a broad range of timescales and luminosities to the observed properties of the AGN population as a whole, and particularly the Eddington ratio distribution function (ERDF). We have implemented our framework on GPU architecture, relying on previously published time series generating algorithms. After extensive tests that characterise several intrinsic and numerical aspects of the simulations, we describe some applications used for current and future time domain surveys and for the study of extremely variable sources (e.g., "changing look" or flaring AGN). Specifically, we define a simulation setup which reproduces the AGN variability observed in the PTF/iPTF survey, and use it to forward model longer light curves of the kind that may be observed within the LSST main survey. Thanks to our effcient implementations, these simulations are able to cover for example over 1 Myr with a roughly weekly cadence. We envision that this framework will become highly valuable to prepare for, and best exploit, data from upcoming time domain surveys, such as for example LSST.

We present a framework for modelling the star-formation histories of galaxies as a stochastic process. We define this stochastic process through a power spectrum density with a functional form of a broken power-law. Star-formation histories are correlated on short timescales, the strength of this correlation described by a power-law slope, $\alpha$, and they decorrelate to resemble white noise over a timescale that is proportional to the timescale of the break in the power spectrum density, $\tau_{\rm break}$. We use this framework to explore the properties of the stochastic process that, we assume, gives rise to the log-normal scatter about the relationship between star-formation rate and stellar mass, the so-called galaxy star-forming main sequence. Specifically, we show how the measurements of the normalisation and width ($\sigma_{\rm MS}$) of the main sequence, measured in several passbands that probe different timescales, give a constraint on the parameters of the underlying power spectrum density. We first derive these results analytically for a simplified case where we model observations by averaging over the recent star-formation history. We then run numerical simulations to find results for more realistic observational cases. As a proof of concept, we use observational estimates of the main sequence scatter at $z\sim0$ and $M_{\star}\approx10^{10}~M_{\odot}$ measured in H$\alpha$, UV+IR and the u-band, and show that combination of these point to $\tau_{\rm break}=178^{+104}_{-66}$ Myr, when assuming $\alpha=2$. This implies that star-formation histories of galaxies lose "memory" of their previous activity on a timescale of $\sim200$ Myr, highlighting the importance of baryonic effects that act over the dynamical timescales of galaxies.

In order to relate the observed evolution of the galaxy stellar mass function and the luminosity function of active galactic nuclei (AGN), we explore a co-evolution scenario in which AGN are associated only with the very last phases of the star-forming life of a galaxy. We derive analytically the connections between the parameters of the observed quasar luminosity functions and galaxy mass functions. The $(m_{\rm bh}/m_{*})_{Qing}$ associated with quenching is given by the ratio of the global black hole accretion rate density (BHARD) and star-formation rate density (SFRD) at the epoch in question. Observational data on the SFRD and BHARD suggests $(m_{\rm bh}/m_{*})_{Qing} \propto (1+z)^{1.5}$ below redshift 2. This evolution reproduces the observed mass-luminosity plane of SDSS quasars, and also reproduces the local $m_{\rm bh}/m_{*}$ relation in passive galaxies. The characteristic Eddington ratio, $\lambda^*$, is derived from both the BHARD/SFRD ratio and the evolving $L^*$ of the AGN population. This increases up to $z \sim 2$ as $\lambda^* \propto (1+z)^{2.5}$ but at higher redshifts, $\lambda^*$ stabilizes at the physically interesting Eddington limit, $\lambda^* \sim 1$. The new model may be thought of as an opposite extreme to our earlier co-evolution scenario in Caplar et al. 2015. The main observable difference between the two co-evolution scenarios, presented here and in Caplar et al. 2015, is in the active fraction of low mass star-forming galaxies. We compare the predictions with the data from deep multi-wavelength surveys and find that the "quenching" scenario developed in the current paper is much to be preferred.

We present a framework to link and describe AGN variability on a wide range of timescales, from days to billions of years. In particular, we concentrate on the AGN variability features related to changes in black hole fuelling and accretion rate. In our framework, the variability features observed in different AGN at different timescales may be explained as realisations of the same underlying statistical properties. In this context, we propose a model to simulate the evolution of AGN light curves with time based on the probability density function (PDF) and power spectral density (PSD) of the Eddington ratio ($L/L_{\rm Edd}$) distribution. Motivated by general galaxy population properties, we propose that the PDF may be inspired by the $L/L_{\rm Edd}$ distribution function (ERDF), and that a single (or limited number of) ERDF+PSD set may explain all observed variability features. After outlining the framework and the model, we compile a set of variability measurements in terms of structure function (SF) and magnitude difference. We then combine the variability measurements on a SF plot ranging from days to Gyr. The proposed framework enables constraints on the underlying PSD and the ability to link AGN variability on different timescales, therefore providing new insights into AGN variability and black hole growth phenomena.

We use stellar mass functions to study the properties and the significance of quenching through major galaxy mergers. In addition to SDSS DR7 and Galaxy Zoo 1 data, we use samples of visually selected major galaxy mergers and post merger galaxies. We determine the stellar mass functions of the stages that we would expect major merger quenched galaxies to pass through on their way from the blue cloud to the red sequence: 1: major merger, 2: post merger, 3: blue early type, 4: green early type and 5: red early type. Based on the similar mass function shapes we conclude that major mergers are likely to form an evolutionary sequence from star formation to quiescence via quenching. Relative to all blue galaxies, the major merger fraction increases as a function of stellar mass. Major merger quenching is inconsistent with the mass and environment quenching model. At z~0 major merger quenched galaxies are unlikely to constitute the majority of galaxies that transition the green valley. Furthermore, between z~0-0.5 major merger quenched galaxies account for 1-5% of all quenched galaxies at a given stellar mass. Major galaxy mergers are therefore not a significant quenching pathway, neither at z~0 nor within the last 5 Gyr. The majority of red galaxies must have been quenched through an alternative quenching mechanism which causes a slow blue to red evolution.

We use a phenomenological model to show that black hole growth in the local Universe (z < 0.1) can be described by two separate, mass independent Eddington ratio distribution functions (ERDFs). We assume that black holes can be divided into two independent groups: those with radiatively efficient accretion, primarily hosted by optically blue and green galaxies, and those with radiatively inefficient accretion, which are mainly found in red galaxies. With observed galaxy stellar mass functions as input, we show that the observed AGN luminosity functions can be reproduced by using mass independent, broken power law shaped ERDFs. We use the observed hard X-ray and 1.4 GHz radio luminosity functions to constrain the ERDF for radiatively efficient and inefficient AGN, respectively. We also test alternative ERDF shapes and mass dependent models. Our results are consistent with a mass independent AGN fraction and AGN hosts being randomly drawn from the galaxy population. We argue that the ERDF is not shaped by galaxy-scale effects, but by how efficiently material can be transported from the inner few parsecs to the accretion disc. Our results are incompatible with the simplest form of mass quenching where massive galaxies host higher accretion rate AGN. Furthermore, if reaching a certain Eddington ratio is a sufficient condition for maintenance mode, it can occur in all red galaxies, not just the most massive ones.

We characterize the optical variability of quasars in the intermediate Palomar Transient Factory (iPTF) and Palomar Transient Factory (PTF) surveys. We re-calibrate the $r$-band light curves for $\sim$28,000 luminous, broad-line AGNs from the SDSS, producing a total of $\sim$2.4 million photometric data points. We utilize both the structure function (SF) and power spectrum density (PSD) formalisms to search for links between the optical variability and the physical parameters of the accreting supermassive black holes that power the quasars. The excess variance (SF$^{2}$) of the quasar sample tends to zero at very short time separations, validating our re-calibration of the time-series data. We find that the the amplitude of variability at a given time-interval, or equivalently the time-scale of variability to reach a certain amplitude, is most strongly correlated with luminosity with weak or no dependence on black hole mass and redshift. For a variability level of SF($\tau$)=0.07 mag, the time-scale has a dependency of $\tau \propto L^{0.4}$. This is broadly consistent with the expectation from a simple Keplerian accretion disk model, which provides $\tau \propto L^{0.5}$. The PSD analysis also reveals that many quasar light curves are steeper than a damped random walk. We find a correlation between the steepness of the PSD slopes, specifically the fraction of slopes steeper than 2.5, and black hole mass, although we cannot exclude the possibility that luminosity or Eddington ratio are the drivers of this effect. This effect is also seen in the SF analysis of the (i)PTF data, and in a PSD analysis of quasars in the SDSS Stripe 82.

We analyze the role of first (leading) author gender on the number of citations that a paper receives, on the publishing frequency and on the self-citing tendency. We consider a complete sample of over 200,000 publications from 1950 to 2015 from five major astronomy journals. We determine the gender of the first author for over 70% of all publications. The fraction of papers which have a female first author has increased from less than 5% in the 1960s to about 25% today. We find that the increase of the fraction of papers authored by females is slowest in the most prestigious journals such as Science and Nature. Furthermore, female authors write 19$\pm$7% fewer papers in seven years following their first paper than their male colleagues. At all times papers with male first authors receive more citations than papers with female first authors. This difference has been decreasing with time and amounts to $\sim$6% measured over the last 30 years. To account for the fact that the properties of female and male first author papers differ intrinsically, we use a random forest algorithm to control for the non-gender specific properties of these papers which include seniority of the first author, number of references, total number of authors, year of publication, publication journal, field of study and region of the first author's institution. We show that papers authored by females receive 10.4$\pm$0.9% fewer citations than what would be expected if the papers with the same non-gender specific properties were written by the male authors. Finally, we also find that female authors in our sample tend to self-cite more, but that this effect disappears when controlled for non-gender specific variables.

We explore the connections between the evolving galaxy and AGN populations. We present a simple phenomenological model that links the evolving galaxy mass function and the evolving quasar luminosity function, which makes specific and testable predictions for the distribution of host galaxy masses for AGN of different luminosities. We show that the $\phi^{*}$ normalisations of the galaxy mass function and of the AGN luminosity function closely track each other over a wide range of redshifts, implying a constant "duty cycle" of AGN activity. The strong redshift evolution in the AGN $L^*$ can be produced by either an evolution in the distribution of Eddington ratios, or in the $m_{bh}/m_{*}$ mass ratio, or both. To try to break this degeneracy we look at the distribution of AGN in the SDSS ($m_{bh},L$) plane, showing that an evolving ratio $m_{bh}/m_{*} \propto (1+z)^2$ reproduces the observed data and also reproduces the local relations which connect the black hole population with the host galaxies for both quenched and star-forming populations. We stress that observational studies that compare the masses of black holes in active galaxies at high redshift with those in quiescent galaxies locally will always see much weaker evolution. Evolution of this form would produce, or could be produced by, a redshift-independent $m_{bh} - \sigma$ relation and could explain why the local $m_{bh} - \sigma$ relation is tighter than $m_{bh} - m_{*}$ even if $\sigma$ is not directly linked to black hole growth. Irrespective of the evolution of $m_{bh}/m_{*}$, the model reproduces both the appearance of "downsizing" and the so-called "sub-Eddington boundary" without any mass-dependence in the evolution of black hole growth rates.

A model of unification of dark matter and dark energy based on the modeling of the speed of sound as a function of the parameter of the equation of state is introduced. It is found that the model in which the speed of sound depends on the power of the parameter of the equation of state, $c_s^2=\alpha (-w)^{\gamma}$, contains the generalized Chaplygin gas models as its subclass. An effective scalar field description of the model is obtained in a parametric form which in some cases can be translated into a closed form solution for the scalar field potential. A constraint on model parameters is obtained using the observational data on the Hubble parameter at different redshifts.