We show how wet major mergers can create disc galaxies in a cosmological context, using the Illustris simulation. We select a sample of 38 disc galaxies having experienced a major merger in their history with no subsequent significant minor merger, and study how they transform into discs after the merger. In agreement with what was previously found in controlled simulations of such mergers, we find that their disc is built gradually from young stars formed after the merger in the disc region, while the old stars born before the merger form an ellipsoidal component. Focusing on one fiducial case from our sample, we show how the gas was initially dispersed in the halo region right after the merger, but is then accreted onto a disc to form stars, and builds the disc component. We then select a sample of major mergers creating elliptical galaxies, to show that those cases correspond mainly to dry mergers, where the lack of star formation prevents the formation of a disc in the remnant galaxy. The amount of gas in the remnant galaxy after the merger is therefore essential to determine the final outcome of a major merger.

We study barred galaxies selected from the Illustris cosmological simulation, focusing on tidally induced bars formed from flyby interactions. To guarantee high enough resolution we focus on high mass disc galaxies ($M_*>8.3 \times 10^{10} M_{\odot}$). We find that the fraction of barred galaxies among those (21% at redshift $z=0$) is lower in Illustris than observed in the local Universe, and the fraction grows slightly with redshift. The bar fraction also increases with the stellar mass and decreases with the amount of gas in the disc. Only very few bars at redshift $z=0$ are formed in secular evolution ($\sim$ 7%) and most of them are triggered by external perturbers in mergers or flybys. Many of these bars disappear over time, mostly during secular evolution, which leads to a lower fraction of bars at redshift $z=0$. We then focus on the effect of flyby interactions on the disc and look at tidally induced bars created by a flyby, or pre-existing bars influenced by the passage of a perturber. In the latter case, the interaction can enhance or weaken the bar. During the interaction, the change in the bar strength occurs right after the pericentre passage. The resulting tidally induced bars tend to be stronger than the overall bar sample in Illustris. The preferred scenario to create or enhance a bar seems to be with a strong interaction involving a perturber on a prograde orbit. Furthermore, the strength of the created bar grows with the strength of the interaction.

We study the effect of angular momentum on the surface density profiles of disc galaxies, using high resolution simulations of major mergers whose remnants have downbending radial density profiles (type II). As described in the previous papers of this series, in this scenario, most of the disc mass is acquired after the collision via accretion from a hot gaseous halo. We find that the inner and outer disc scalelengths, as well as the break radius, correlate with the total angular momentum of the initial merging system, and are larger for high angular momentum systems. We follow the angular momentum redistribution in our simulated galaxies, and find that, like the mass, the disc angular momentum is acquired via accretion, i.e. to the detriment of the gaseous halo. Furthermore, high angular momentum systems give more angular momentum to their discs, which affects directly their radial density profile. Adding simulations of isolated galaxies to our sample, we find that the correlations are valid also for disc galaxies evolved in isolation. We show that the outer part of the disc at the end of the simulation is populated mainly by inside-out stellar migration, and that in galaxies with higher angular momentum, stars travel radially further out. This, however, does not mean that outer disc stars (in type II discs) were mostly born in the inner disc. Indeed, generally the break radius increases over time, and not taking this into account leads to overestimating the number of stars born in the inner disc.

Context: In a series of papers, we study the major merger of two disk galaxies in order to establish whether or not such a merger can produce a disc galaxy. Aims: Our aim here is to describe in detail the technical aspects of our numerical experiments. Methods: We discuss the initial conditions of our major merger, which consist of two protogalaxies on a collision orbit. We show that such merger simulations can produce a non-realistic central mass concentration, and we propose simple, parametric, AGN-like feedback as a solution to this problem. Our AGN-like feedback algorithm is very simple: at each time-step we take all particles whose local volume density is above a given threshold value and increase their temperature to a preset value. We also compare the GADGET3 and GIZMO codes, by applying both of them to the same initial conditions. Results: We show that the evolution of isolated protogalaxies resembles the evolution of disk galaxies, thus arguing that our protogalaxies are well suited for our merger simulations. We demonstrate that the problem with the unphysical central mass concentration in our merger simulations is further aggravated when we increase the resolution. We show that our AGN-like feedback removes this non-physical central mass concentration, and thus allows the formation of realistic bars. Note that our AGN-like feedback mainly affects the central region of a model, without significantly modifying the rest of the galaxy. We demonstrate that, in the context of our kind of simulation, GADGET3 gives results which are very similar to those obtained with the PSPH (density independent SPH) flavor of GIZMO. Moreover, in the examples we tried, the differences between the results of the two flavors of GIZMO, namely PSPH, and MFM (mesh-less algorithm) are similar to and, in some comparisons, larger than the differences between the results of GADGET3 and PSPH.

Using three fiducial Nbody+SPH simulations, we follow the merging of two disk galaxies with a hot gaseous halo component each, and examine whether the merger remnant can be a spiral galaxy. The stellar progenitor disks are destroyed by violent relaxation during the merging and most of their stars form a classical bulge, while the remaining form a thick disk and its bar. A new stellar disk forms subsequently and gradually in the remnant from the gas accreted mainly from the halo. It is vertically thin and well extended in its equatorial plane. A bar starts forming before the disk is fully in place, contrary to what is assumed in idealised simulations of isolated bar-forming galaxies. It has morphological features such as ansae and boxy/peanut bulges. Stars of different ages populate different parts of the box/peanut. A disky pseudobulge forms also, so that by the end of the simulation, all three types of bulges coexist. The oldest stars are found in the classical bulge, followed by those of the thick disk, then by those in the thin disk. The youngest stars are in the spiral arms and the disky pseudobulge. The disk surface density profiles are of type II (exponential with downbending), and the circular velocity curves are flat and show that the disks are submaximum in these examples: two clearly so and one near-borderline between maximum and submaximum. On average, only roughly between 10 and 20% of the stellar mass is in the classical bulge of the final models, i.e. much less than in previous simulations.