We present the novel wide & deep neural network GalaxyNet, which connects the properties of galaxies and dark matter haloes, and is directly trained on observed galaxy statistics using reinforcement learning. The most important halo properties to predict stellar mass and star formation rate (SFR) are halo mass, growth rate, and scale factor at the time the mass peaks, which results from a feature importance analysis with random forests. We train different models with supervised learning to find the optimal network architecture. GalaxyNet is then trained with a reinforcement learning approach: for a fixed set of weights and biases, we compute the galaxy properties for all haloes and then derive mock statistics (stellar mass functions, cosmic and specific SFRs, quenched fractions, and clustering). Comparing these statistics to observations we get the model loss, which is minimised with particle swarm optimisation. GalaxyNet reproduces the observed data very accurately ($\chi_\mathrm{red}=1.05$), and predicts a stellar-to-halo mass relation with a lower normalisation and shallower low-mass slope at high redshift than empirical models. We find that at low mass, the galaxies with the highest SFRs are satellites, although most satellites are quenched. The normalisation of the instantaneous conversion efficiency increases with redshift, but stays constant above $z\gtrsim0.7$. Finally, we use GalaxyNet to populate a cosmic volume of $(5.9~\mathrm{Gpc})^3$ with galaxies and predict the BAO signal, the bias, and the clustering of active and passive galaxies up to $z=4$, which can be tested with next-generation surveys, such as LSST and Euclid.

We explore the galaxy-galaxy merger rate with the empirical model for galaxy formation, EMERGE. On average, we find that between $2$ per cent and $6$ per cent of massive galaxies ($\log_{10}(m_{*}/M_{\odot}) \geq 10.3$) will experience a major merger per Gyr. Our model predicts galaxy merger rates that do not scale as a power-law with redshift, and exhibit a clear stellar mass and mass-ratio dependence, unlike the underlying halo-halo merger rate. Specifically, major mergers are more frequent at high masses and at low redshift. We show mergers are significant for the stellar mass growth of galaxies $\log_{10}(m_{*}/M_{\odot}) \gtrsim 11.0$. For the most massive galaxies major mergers dominate the accreted mass fraction, contributing as much as $90$ per cent of the total accreted stellar mass. We argue that these phenomena are a direct result of the stellar-to-halo mass relation, which results in massive galaxies having a higher likelihood of experiencing major mergers than low mass galaxies. Furthermore, we find no statistical evidence for mergers as a driving mechanism for quenching, with more than $58$ per cent of the most massive galaxies experiencing no major merger within $2$ Gyr of quenching. Our model produces a galaxy pair fraction consistent with resent observations, exhibiting a form best described by a power-law exponential function. Translating these pair fractions into merger rates results in an over prediction compared to the model intrinsic values. We find the pair fraction can be best mapped to the intrinsic merger rate by adopting a constant observation timescale $T_{obs}=2.56$ Gyr, assuming all observed pairs merge by $z=0$.