We present a sample of 33 Balmer Break Galaxies (BBGs) selected as HST/F160W dropouts in the deepest CANDELS/GOODS fields ($H\gtrsim27.3$~mag) but relatively bright in {\it Spitzer}/IRAC ($[3.6],[4.5]<24.5$~mag), implying red colors (median and quartiles: $\langle H-[3.6]\rangle=3.1^{3.4}_{2.8}$\,mag). Half of these BBGs are newly identified sources. Our BBGs are massive ($\langle \log(\rm{M}/\rm{M}_\odot)\rangle=10.8$) high redshift ($\langle z\rangle=4.8$) dusty ($\langle \rm{A(V)}\rangle=2.0$~mag) galaxies. The SEDs of half of our sample indicate that they are star-forming galaxies with typical specific SFRs 0.5-1.0~Gyr$^{-1}$, qualifying them as main sequence (MS) galaxies at $3<z<6$. One third of those SEDs indicates the presence of prominent emission lines (H$\beta$+$[OIII]$, H$\alpha$$+$[NII]) boosting the IRAC fluxes and red colors. Approximately 20\% of the BBGs are very dusty ($\rm{A(V)}\sim2.5$~mag) starbursts with strong mid-to-far infrared detections and extreme SFRs ($\rm{SFR}>10^{3}\,\rm{M}_\odot/yr$) that place them above the MS. The rest, 30\%, are post-starbursts or quiescent galaxies located $>2\sigma$ below the MS with mass-weighted ages older than 700~Myr. Only 2 of the 33 galaxies are X-ray detected AGN with optical/near-infrared SEDs dominated by stellar emission, but the presence of obscured AGN in the rest of sources cannot be discarded. Our sample accounts for 8\% of the total number density of $\log(\rm{M}/\rm{M}_\odot)>10$ galaxies at $z>3$, but it is a significant contributor (30\%) to the general population of red $\log(\rm{M}/\rm{M}_\odot)>11$ galaxies at $4<z<6$. Finally, our results point out that 1 of every 30 massive $\log(\rm{M}/\rm{M}_\odot)>11$ galaxies in the local Universe was assembled in the first 1.5~Gyr after the Big Bang, a fraction that is not reproduced by state-of-the-art galaxy formation simulations.

The Hubble Ultra Deep field (HUDF) is the deepest region ever observed with the Hubble Space Telescope. With the main objective of unveiling the nature of galaxies up to $z \sim 7-8$, the observing and reduction strategy have focused on the properties of small and unresolved objects, rather than the outskirts of the largest objects, which are usually over-subtracted. We aim to create a new set of WFC3/IR mosaics of the HUDF using novel techniques to preserve the properties of the low surface brightness regions. We created ABYSS: a pipeline that optimises the estimate and modelling of low-level systematic effects to obtain a robust background subtraction. We have improved four key points in the reduction: 1) creation of new absolute sky flat fields, 2) extended persistence models, 3) dedicated sky background subtraction and 4) robust co-adding. The new mosaics successfully recover the low surface brightness structure removed on the previous HUDF published reductions. The amount of light recovered with a mean surface brightness dimmer than $\overline{\mu}=26$ mar arcsec$^{-2}$ is equivalent to a m=19 mag source when compared to the XDF and a m=20 mag compared to the HUDF12. We present a set of techniques to reduce ultra-deep images ($\mu>32.5$ mag arcsec$^{-2}$, $3\sigma$ in $10\times10$ arcsec boxes), that successfully allow to detect the low surface brightness structure of extended sources on ultra deep surveys. The developed procedures are applicable to HST, JWST, EUCLID and many other space and ground-based observatories. We will make the final ABYSS WFC3/IR HUDF mosaics publicly available at http://www.iac.es/proyecto/abyss/.

Major mergers are popularly considered too destructive to produce the relaxed regular structures and the morphological inner components (ICs) usually observed in lenticular (S0) galaxies. We aim to test if major mergers can produce remnants with realistic S0 morphologies. We have selected a sample of relaxed discy remnants resulting from the dissipative merger simulations of the GalMer database and derived their properties mimicking the typical conditions of current observational data. We compare their global morphologies, visual components, and merger relics in mock photometric images with their real counterparts. Only $\sim$1-2 Gyr after the full merger, we find that: 1) many remnants (67 major and 29 minor events) present relaxed structures and typical S0 or E/S0 morphologies, for a wide variety of orbits and even in gas-poor cases. 2) Contrary to popular expectations, most of them do not exhibit any morphological traces of their past merger origin under typical observing conditions and at distances as nearby as 30 Mpc. 3) The merger relics are more persistent in minor mergers than in major ones for similar relaxing time periods. 4) No major-merger S0-like remnant develops a significant bar. 5) Nearly 58% of the major-merger S0 remnants host visually detectable ICs, such as embedded inner discs, rings, pseudo-rings, inner spirals, nuclear bars, and compact sources, very frequent in real S0s too. 6) All remnants contain a lens or oval, identically ubiquitous in local S0s. 7) These lenses and ovals do not come from bar dilution in major merger cases, but are associated with stellar halos or embedded inner discs instead (thick or thin). We conclude that the relaxed morphologies, lenses, ovals, and other ICs of real S0s do not necessarily come from internal secular evolution, gas infall or environmental mechanisms, as traditionally assumed, but they can result from major mergers as well.

We have undertaken a comprehensive search for both Lyman Alpha Emitters (LAEs) and Lyman Break Galaxies (LBGs) in the SHARDS Survey of the GOODS-N field. SHARDS is a deep imaging survey, made with the 10.4 m Gran Telescopio Canarias (GTC), employing 25 medium band filters in the range from 500 to 941 nm. This is the first time that both LAEs and LBGs are surveyed simultaneously in a systematic way in a large field. We draw a sample of 1558 sources; 528 of them are LAEs. Most of the sources (1434) show rest-frame UV continua. A minority of them (124) are pure LAEs with virtually no continuum detected in SHARDS. We study these sources from $z\sim3.35$ up to $z\sim6.8$, well into the epoch of reionization. Note that surveys done with just one or two narrow band filters lack the possibility to spot the rest-frame UV continuum present in most of our LAEs. We derive redshifts, Star Formation Rates (SFRs), Ly$\alpha$ Equivalent Widths (EWs) and Luminosity Functions (LFs). Grouping within our sample is also studied, finding 92 pairs or small groups of galaxies at the same redshift separated by less than 60 comoving kpc. In addition, we relate 87 and 55 UV-selected objects with two known overdensities at $z=4.05$ and $z=5.198$, respectively. Finally, we show that surveys made with broad band filters are prone to introduce many unwanted sources ($\sim20$% interlopers), which means that previous studies may be overestimating the calculated LFs, specially at the faint end.

The physical processes driving the evolution of star formation (SF) in galaxies over cosmic time still present many open questions. Recent galaxy surveys allow now to study these processes in great detail at intermediate redshift. In this work, we build a complete sample of star-forming galaxies and analyze their properties, reaching systems with low stellar masses and low star formation rates (SFRs) at intermediate-to-low redshift. We use data from the SHARDS multiband survey in the GOODS-North field. Its depth (up to magnitude $\langle m_{3\sigma}\rangle\sim26.5$) and its spectro-photometric resolution ($R\sim50$) provides us with an ideal dataset to search for emission line galaxies (ELGs). We develop a new algorithm to identify low-redshift ($z$<0.36) ELGs by detecting the [OIII]5007 and $H\alpha$ emission lines simultaneously. We fit the spectral energy distribution (SED) of the selected sample, using a model with two single stellar populations. We find 160 star-forming galaxies with equivalent widths (EWs) as low as 12 {\AA}, with median values for the sample of $\sim$ 35 {\AA} in [OIII]5007 and $\sim$ 56 {\AA} in $H\alpha$, respectively. Results from the SED fitting show a young stellar population with low median metallicity (36% of the solar value) and extinction ($A_V \sim$ 0.37), with median galaxy stellar mass $\sim$ 10$^{8.5}$ M$_{\odot}$. Gas-phase metallicities measured from available spectra are also low. ELGs in our sample present bluer colors in the UVJ plane than the median color-selected star-forming galaxy in SHARDS. We suggest a new (V-J) color criterion to separate ELGs from non-ELGs in blue galaxy samples. In addition, several galaxies present high densities of O-type stars. Robust fits to the full SEDs can only be obtained including an old stellar population, suggesting the young component is built up by a recent burst of SF in an otherwise old galaxy.

Type-III S0 galaxies present tight scaling relations between their surface brightness photometric and structural parameters. Several evolutionary models have been proposed for the formation of Type-III S0 galaxies but the observations of are usually limited to the local Universe. We study the evolution of the photometric and structural scaling relations found between the parameters of the surface brightness profiles in the rest-frame R-band of Type-III S0 galaxies with z and the possible differences between the rest-frame (B-R) colours of the inner and outer disc profiles. We make use of a sample of 14 Type-III E/S0--S0 galaxies at 0.2<z<0.6 to study if the correlations found in local Type-III S0 galaxies were present ~6 Gyr ago. We analyse the distribution of the surface brightness characteristic parameters as a function of the stellar mass and if there is a significant change with z. We also derive their rest-frame (B-R) colour profiles and we compare these results with the predictions from a grid of SSP models. We find that the inner and outer scale-lengths of Type-III S0 galaxies at 0.4<z<0.6 follow compatible trends and scaling relations with those observed in local S0 galaxies. We do not detect any significant differences between the location of Rbreak between z~0.6 and z=0 for a fixed stellar mass of the object, whereas the surface brightness at the break radius is ~1.5 mag arcsec-2 dimmer in the local Universe than at z~0.6 for a fixed stellar mass. In contrast to Type-II profiles, the Type-III surface brightness profiles of S0 galaxies present compatible Rbreak values and scaling relations during the last 6 Gyr. This result and the similarity of the colours of the inner and outer discs point to a highly scalable and stable formation process, probably more related to gravitational and dynamical processes than to the evolution of stellar populations (abridged).

The intrinsic shape of galactic bulges in nearby galaxies provides crucial information to separate bulge types. We intended to derive accurate constraints to the intrinsic shape of bulges to provide new clues on their formation mechanisms and set new limitations for future simulations. We retrieved the intrinsic shape of a sample of CALIFA bulges using a statistical approach. Taking advantage of GalMer numerical simulations of binary mergers we estimated the reliability of the procedure. Analyzing the $i$-band mock images of resulting lenticular remnants, we studied the intrinsic shape of their bulges at different galaxy inclinations. Finally, we introduced a new ($B/A$, $C/A$) diagram to analyze possible correlations between the intrinsic shape and the properties of bulges. We tested the method on simulated lenticular remnants, finding that for galaxies with inclinations $25^{\circ} < \theta < 65^{\circ}$ we can safely derive the intrinsic shape of their bulges. We found that our CALIFA bulges tend to be nearly oblate systems (66%), with a smaller fraction of prolate spheroids (19%) and triaxial ellipsoids (15%). The majority of triaxial bulges are in barred galaxies (75%). Moreover, we found that bulges with low S\'ersic indices or in galaxies with low bulge-to-total luminosity ratios form a heterogeneous class of objects; additionally, also bulges in late-type galaxies or in less massive galaxies have no preference in being oblate, prolate, or triaxial. On the contrary, bulges with high S\'ersic index, in early-type galaxies, or in more massive galaxies are mostly oblate systems. We concluded that various evolutionary pathways may coexist in galaxies, with merging events and dissipative collapse being the main mechanisms driving the formation of the most massive oblate bulges and bar evolution reshaping the less massive triaxial bulges.

(Abridged version) We explore whether a scenario that combines an origin by mergers at $z\sim$1.8-1.5 with a subsequent passive evolution of the resulting S0 remnants since $z \sim$0.8-1 is compatible with observational data of S0s in the Tully-Fisher relation (TFR). We studied a set of major and minor merger experiments from the GalMer database that generate massive S0 remnants. We analysed the location of these remnants in the photometric and stellar TFRs assuming that they correspond to $z\sim0.8$ galaxies. We then estimated their evolution in these planes over the last 7 Gyr. The results were compared with data of real S0s and spirals at different redshifts. We also tested how the use of Vcirc or Vrot,max affects the results. We found that just after $\sim$1-2 Gyr of coalescence, major mergers generate S0 remnants that are outliers of the local photometric and stellar TFRs at $z\sim0.8$. After $\sim$4-7 Gyr of passive evolution in isolation, the S0 remnants move towards the local TFR, although the initial scatter among them persists. This scatter is sensitive to the indicator used for the rotation velocity: Vcirc values yield a lower scatter than when Vrot,max values are considered instead. In the planes involving Vrot,max, a clear segregation of the S0 remnants in terms of the spin-orbit coupling of the model is observed, in which the remnants of retrograde encounters overlap with local S0s hosting counter-rotating discs. The location of the S0 remnants at $z\sim 0$ agrees well with the observed distribution of local S0 galaxies in the $\sigma_0$-$M_K$, Vcirc-$\sigma_0$ and Vrot,max-$\sigma_0$ planes. Thus, massive S0 galaxies may have been formed through major mergers that occurred at high redshift and have later evolved towards the local TFR through passive evolution in relative isolation, a mechanism that would also contribute to the scatter observed in this relation.

The controversy about the origin of the structure of S0--E/S0 galaxies may be due to the difficulty of comparing surface brightness profiles with different depths, photometric corrections and PSF effects (almost always ignored). We aim to quantify the properties of Type-III (anti-truncated) discs in a sample of S0 galaxies at 0.2<z<0.6. In this paper, we present the sample selection and describe in detail the methods to robustly trace the structure in their outskirts and correct for PSF effects. We have selected and classified a sample of 150 quiescent galaxies at 0.2<z<0.6 in the GOODS-N field. We perform a quantitative structural analysis of 44 S0-E/S0 galaxies. We corrected their surface brightness profiles for PSF distortions and analysed the biases in the structural and photometric parameters when the PSF correction is not applied. Additionally, we have developed Elbow, an automatic statistical method to determine whether a possible break is significant - or not - and its type and made it publicly available. We found 14 anti-truncated S0-E/S0 galaxies in the range 0.2<z<0.6 (~30% of the final sample). This fraction is similar to the those reported in the local Universe. In our sample, ~25% of the Type-III breaks observed in PSF-uncorrected profiles are artifacts, and their profiles turn into a Type I after PSF correction. PSF effects also soften Type-II profiles. We found that the profiles of Type-I S0 and E/S0 galaxies of our sample are compatible with the inner profiles of the Type-III, in contrast with the outer profiles. We have obtained the first robust and reliable sample of 14 anti-truncated S0--E/S0 galaxies beyond the local Universe, in the range 0.2<z<0.6. PSF effects significantly affect the shape of the surface brightness profiles in galaxy discs even in the case of the narrow PSF of HST/ACS images, so future studies on the subject should make an effort to correct them.

Previous numerical studies had apparently ruled out the possibility that flares in galaxy discs could give rise to the apparent breaks in their luminosity profiles when observed edge-on. However the studies have not, until now, analyzed this hypothesis systematically using realistic models for the disc, the flare, and the bulge. We revisit this theme by analyzing a series of models which sample a wide range of observationally based structural parameters for these three components. We have considered realistic distributions of bulge to disc ratios, morphological parameters of bulges and discs, vertical scale heights of discs and their radial gradients defining the flare for different morphological types and stellar mass bins, based on observations. The surface brightness profiles for the face-on and edge-on views of each model were simulated to find out whether the flared disc produces a Type-II break in the disc profile when observed edge-on, and if so under what conditions. Contrary to previous claims, we find that discs with realistic flares can produce significant breaks in discs when observed edge-on. Specifically a flare with the parameters of that of the Milky Way would produce a significant break of the disc at a Rbreak of ~8.6 kpc if observed edge-on. Central bulges have no significant effects on the results. These simulations show that flared discs can explain the existence of many Type-II breaks observed in edge-on galaxies, in a range of galaxies with low-to-intermediate break strength values of -0.25<S<-0.1.

Lenticular galaxies (S0s) represent the majority of early-type galaxies in the local Universe, but their formation channels are still poorly understood. While galaxy mergers are obvious pathways to suppress star formation and increase bulge sizes, the marked parallelism between spiral and lenticular galaxies (e.g. photometric bulge-disc coupling) seemed to rule out a potential merger origin. Here, we summarise our recent work in which we have shown, through N-body numerical simulations, that disc-dominated lenticulars can emerge from major mergers of spiral galaxies, in good agreement with observational photometric scaling relations. Moreover, we show that mergers simultaneously increase the light concentration and reduce the angular momentum relative to their spiral progenitors. This explains the mismatch in angular momentum and concentration between spirals and lenticulars recently revealed by CALIFA observations, which is hard to reconcile with simple fading mechanisms (e.g. ram-pressure stripping).

We present Star Formation Histories (SFHs) for a sample of 104 massive (stellar mass M > 10$^{10}$ M$_{\odot}$) quiescent galaxies (MQGs) at $z$=1.0-1.5 from the analysis of spectro-photometric data from the SHARDS and HST/WFC3 G102 and G141 surveys of the GOODS-N field, jointly with broad-band observations from ultraviolet (UV) to far-infrared (Far-IR). The sample is constructed on the basis of rest-frame UVJ colours and specific star formation rates (sSFR=SFR/Mass). The Spectral Energy Distributions (SEDs) of each galaxy are compared to models assuming a delayed exponentially declining SFH. A Monte Carlo algorithm characterizes the degeneracies, which we are able to break taking advantage of the SHARDS data resolution, by measuring indices such as MgUV and D4000. The population of MQGs shows a duality in their properties. The sample is dominated (85%) by galaxies with young mass-weighted ages, tM $<$ 2 Gyr, short star formation timescales, $\langle \tau \rangle$ $\sim$ 60-200 Myr, and masses log(M/M$_{\odot}$) $\sim$ 10.5. There is an older population (15%) with tM$=$2 - 4 Gyr, longer star formation timescales, $\langle \tau \rangle$$\sim$ 400 Myr, and larger masses, log(M/M$_{\odot}$) $\sim$ 10.7. The SFHs of our MQGs are consistent with the slope and the location of the Main Sequence (MS) of star-forming galaxies at $z$ $>$ 1.0, when our galaxies were 0.5-1.0 Gyr old. According to these SFHs, all the MQGs experienced a Luminous Infrared Galaxy (LIRG) phase that lasts for $\sim$ 500 Myr, and half of them an Ultra Luminous Infrared Galaxy (ULIRG) phase for $\sim$ 100 Myr. We find that the MQG population is almost assembled at $z$ $\sim$ 1, and continues evolving passively with few additions to the population.

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermediate redshifts (z~1). We combine the ultra-deep optical spectro-photometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) with deep UV-to-FIR observations in the GOODS-N field. Exploiting two of the 25 SHARDS medium-band filters, F687W17 and F823W17, we select [OII] emission line galaxies at z~0.84 and z~1.23 and characterize their physical properties. Their rest-frame equivalent widths (EW$_{\mathrm{rf}}$([OII])), line fluxes, luminosities, star formation rates (SFRs) and dust attenuation properties are investigated. The evolution of the EW$_{\mathrm{rf}}$([OII]) closely follows the SFR density evolution of the universe, with a trend of EW$_{\mathrm{rf}}$([OII])$\propto$(1+z)$^3$ up to redshift z~1, followed by a possible flattening. The SF properties of the galaxies selected on the basis of their [OII] emission are compared with complementary samples of SFGs selected by their MIR and FIR emission, and also with a general mass-selected sample of galaxies at the same redshifts. We demonstrate observationally that the UVJ diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent galaxies from the SFGs. The SFR-M$_*$ relation is investigated for the different samples, yelding a logarithmic slope ~1, in good agreement with previous results. The dust attenuations derived from different SFR indicators (UV(1600), UV(2800), [OII], IR) are compared and show clear trends with respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies being affected by stronger dust attenuation.

A number of simulators have argued that major mergers can sometimes preserve discs (e.g. Springel & Hernquist 2005), but the possibility that they could explain the emergence of lenticular galaxies (S0s) has been generally neglected. In fact, observations of S0s reveal a strong structural coupling between their bulges and discs, which seems difficult to reconcile with the idea that they come from major mergers. However, in Querejeta et al. (2015a) we have used N-body simulations of binary mergers to show that, under favourable conditions, discs are first destroyed but soon regrow out of the leftover debris, matching observational photometric scaling relations (e.g. Laurikainen et al. 2010). Additionally, in Querejeta et al. (2015b) we have shown how the merger scenario agrees with the recent discovery that S0s and most spirals are not compatible in an angular momentum--concentration plane. This important result from CALIFA constitutes a serious objection to the idea that spirals transform into S0s mainly by fading (e.g. via ram-pressure stripping, as that would not explain the observed simultaneous change in $\lambda_\mathrm{Re}$ and concentration), but our simulations of major mergers do explain that mismatch. From such a 3D comparison we conclude that mergers must be a relevant process in the build-up of the current population of S0s.

The CALIFA team has recently found that the stellar angular momentum and concentration of late-type spiral galaxies are incompatible with those of lenticular galaxies (S0s), concluding that fading alone cannot satisfactorily explain the evolution from spirals into S0s. Here we explore whether major mergers can provide an alternative way to transform spirals into S0s by analysing the spiral-spiral major mergers from the GalMer database that lead to realistic, relaxed S0-like galaxies. We find that the change in stellar angular momentum and concentration can explain the differences in the $\lambda_\mathrm{Re}$--$R_{90}/R_{50}$ plane found by the CALIFA team. Major mergers thus offer a feasible explanation for the transformation of spirals into S0s.

It has been recently found that the characteristic photometric parameters of antitruncated discs in S0 galaxies follow tight scaling relations. We investigate if similar scaling relations are satisfied by galaxies of other morphological types. We have analysed the trends in several photometric planes relating the characteristic surface brightness and scalelengths of the breaks and the inner and outer discs of local antitruncated S0-Scd galaxies, using published data and fits performed to the surface brightness profiles of two samples of Type-III galaxies in the R and Spitzer 3.6 microns bands. We have performed linear fits to the correlations followed by different galaxy types in each plane, as well as several statistical tests to determine their significance. We have found that: 1) the antitruncated discs of all galaxy types from Sa to Scd obey tight scaling relations both in R and 3.6 microns, as observed in S0s; 2) the majority of these correlations are significant accounting for the numbers of the available data samples; 3) the trends are clearly linear when the characteristic scalelengths are plotted on a logarithmic scale; and 4) the correlations relating the characteristic surface brightnesses of the inner and outer discs and the breaks with the various characteristic scalelengths significantly improve when the latter are normalized to the optical radius of the galaxy. The results suggest that the scaling relations of Type-III discs are independent of the morphological type and the presence (or absence) of bars within the observational uncertainties of the available datasets, although larger and deeper samples are required to confirm this. The tight structural coupling implied by these scaling relations impose strong constraints on the mechanisms proposed for explaining the formation of antitruncated stellar discs in the galaxies across the whole Hubble Sequence (Abridged).

K-band galaxy number counts (GNCs) exhibit a slope change at K~17.5 mag not present in optical bands. To unveil the nature of this feature, we have derived the contribution of different galaxy types to the total K-band GNCs at 0.3<z<1.5 by redshift bins and compared the results with expectations from several galaxy evolutionary models. We show that the slope change is caused by a sudden swap of the galaxy population that numerically dominates the total GNCs (from quiescent E-S0's at K<17.5 mag to blue star-forming discs at fainter magnitudes), and that it is associated with a flattening of the contribution of the E-S0's at 0.6<z<1 to the total GNCs. We confirm previous studies showing that models in which the bulk of massive E-S0's have evolved passively since z>2 cannot predict the slope change, whereas those imposing a relatively late assembly on them (z<1.5) can reproduce it. The K-band GNCs by redshift bins and morphological types point to a progressively definitive build-up of ~50% of this galaxy population at 0.8<z<1.5, which can be explained only through the major mergers reported by observations. We conclude that the slope change in total K-band GNCs is a vestige of the definitive assembly of a substantial fraction of present-day massive E-S0's at 0.8<z<1.5.

Lenticular galaxies (S0s) are more likely to host antitruncated (Type-III) stellar discs than galaxies of later Hubble types. Major mergers are popularly considered too violent mechanisms to form these breaks. We have investigated whether major mergers can result into S0-like remnants with realistic antitruncated stellar discs or not. We have analysed 67 relaxed S0 and E/S0 remnants resulting from dissipative N-body simulations of major mergers from the GalMer database. We have simulated realistic R-band surface brightness profiles of the remnants to identify those with antitruncated stellar discs. Their inner and outer discs and the breaks are quantitatively characterized to compare with real data. Nearly 70% of our S0-like remnants are antitruncated, meaning that major mergers that result in S0s have a high probability of producing Type-III stellar discs. Our remnants lie on top of the extrapolations of the observational trends (towards brighter magnitudes and higher break radii) in several photometric diagrams. In scale-free photometric diagrams, simulations and observations overlap and the remnants reproduce the observational trends, so the physical mechanism after antitruncations is highly scalable. We have found novel photometric scaling relations between the characteristic parameters of the antitruncations in real S0s, which are also reproduced by our simulations. The trends in all the photometric planes can be derived from three basic scaling relations that real and simulated Type-III S0s fulfill: h_i \prop R_brkIII, h_o \prop R_brkIII, and mu_brkIII \prop R_brkIII, where h_i and h_o are the scalelenghts of the inner and outer discs, and mu_brkIII and R_brkIII are the surface brightness and radius of the breaks. Mayor mergers provide a feasible mechanism to form realistic antitruncated S0 galaxies (abridged).

Observations reveal a strong structural coupling between bulge and disc in S0 galaxies, which seems difficult to explain if they have formed from supposedly catastrophic events such as major mergers. We face this question by quantifying the bulge-disc coupling in dissipative simulations of major and minor mergers that result in realistic S0s. We have studied the dissipative N-body binary merger simulations from the GalMer database that give rise to realistic, relaxed E/S0 and S0 remnants (67 major and 29 minor mergers). We simulate surface brightness profiles of these S0-like remnants in the K-band, mimicking typical observational conditions, to perform bulge-disc decompositions analogous to those carried out in real S0s. The global bulge-disc structure of these remnants has been compared with real data, and they distribute in the B/T - r_e - h_d parameter space consistently with real bright S0s, where B/T is the bulge-to-total luminosity ratio, r_e is the bulge effective radius, and h_d is the disc scalelength. Major mergers can rebuild a bulge-disc coupling in the remnants after having destroyed the structures of the progenitors, whereas minor mergers directly preserve them. Remnants exhibit B/T and r_e/h_d spanning a wide range of values, and their distribution is consistent with observations. Many remnants have bulge Sersic indices ranging 1<n<2, flat appearance, and contain residual star formation in embedded discs, a result which agrees with the presence of pseudobulges in real S0s. Contrary to the popular view, mergers (and in particular, major events) can result in S0 remnants with realistically coupled bulge-disc structures in less than ~3 Gyr. In conclusion, the bulge-disc coupling and the presence of pseudobulges in real S0s cannot be used as an argument against the possible major-merger origin of these galaxies.

We analysed collisionless N-body simulations of intermediate and minor dry mergers onto S0s to test whether these mergers can generate S0 galaxies with kinematics intermediate between fast and slow rotators. We find that minor mergers induce a lower decrease of the global rotational support than encounters of lower mass ratios, which results in S0s with properties intermediate between fast and slow rotators. The resulting remnants are intrinsically more triaxial, less flattened, and span the whole range of apparent ellipticities up to $\epsilon_\mathrm{e} \sim 0.8$. They do not show lower apparent ellipticities in random projections than initially; on the contrary, the formation of oval distortions and the disc thickening increase the percentage of projections at $0.4 < \epsilon_\mathrm{e} < 0.7$. In the experiments with S0b progenitor galaxies, minor mergers tend to spin up the bulge and to decrease slightly its intrinsic ellipticity, whereas in the cases of primary S0c galaxies they keep the rotational support of the bulge nearly constant and decrease significantly its intrinsic ellipticity. The remnant bulges remain nearly spherical ($B/A \sim C/A > 0.9$), but exhibit a wide range of triaxialities ($0.20 < T < 1.00$). In the plane of global anisotropy of velocities ($\delta$) vs. intrinsic ellipticity ($\epsilon_\mathrm{e,intr}$), some of our models extend the linear trend found in previous major merger simulations towards higher $\epsilon_\mathrm{e,intr}$ values, while others depart from it. This is consistent with the wide dispersion exhibited by real S0s in this diagram compared with ellipticals, which follow the linear trend drawn by major merger simulations. The different trends exhibited by ellipticals and S0 galaxies in the $\delta$ - $\epsilon_\mathrm{e}$ diagram may be pointing to the different role played by major mergers in the build-up of each morphological type.

Galaxy mergers are considered as questionable mechanisms for the evolution of lenticular galaxies (S0's), on the basis that even minor ones induce structural changes that are difficult to reconcile with the strong bulge-disk coupling observed in the photometric scaling relations of S0's. We check if the evolution induced onto S0's by dry intermediate and minor mergers can reproduce their photometric scaling relations, analysing the bulge-disk decompositions of the merger simulations presented in Eliche-Moral et al. (2012). The mergers induce an evolution in the photometric planes compatible with the data of S0's, even in those ones indicating a strong bulge-disk coupling. The mergers drive the formation of the observed photometric relation in some cases, whereas they induce a slight dispersion compatible with data in others. Therefore, this evolutionary mechanism tends to preserve these scaling relations. In those photometric planes where the morphological types segregate, the mergers always induce evolution towards the region populated by S0's. The structural coupling of the bulge and the disk is preserved or reinforced because the mergers trigger internal secular processes in the primary disk that induce significant bulge growth, even although these models do not induce bars. Intermediate and minor mergers can thus be considered as plausible mechanisms for the evolution of S0's attending to their photometric scaling relations, as they can preserve and even strengthen any pre-existing structural bulge-disk coupling, triggering significant internal secular evolution (even in the absence of bars or dissipational effects). This means that it may be difficult to isolate the effects of pure internal secular evolution from those of the merger-driven one in present-day early-type disks (abridged).

Recent studies have argued that galaxy mergers are not important drivers for the evolution of S0's, on the basis that mergers cannot preserve the coupling between the bulge and disk scale-lengths observed in these galaxies and the lack of correlation of their ratio with the S0 Hubble type. We investigate whether the remnants resulting from collision-less N-body simulations of intermediate and minor mergers onto S0 galaxies evolve fulfilling global structural relations observed between the bulges and disks of these galaxies. Different initial bulge-to-disk ratios of the primary S0 have been considered, as well as different satellite densities, mass ratios, and orbits of the encounter. We have analysed the final morphology of the remnants in images simulating the typical observing conditions of S0 surveys. We derive bulge+disk decompositions of the final remnants to compare their global bulge-to-disk structure with observations. We show that all remnants present undisturbed S0 morphologies according to the prescriptions of specialized surveys. The dry intermediate and minor mergers induce noticeable bulge growth (S0c --> S0b and S0b --> S0a), but affect negligibly to the bulge and disk scale-lengths. Therefore, if a coupling between these two components exists prior to the merger, the encounter does not break this coupling. This fact provides a simple explanation for the lack of correlation between the ratio of bulge and disk scale-lengths and the S0 Hubble type reported by observations. These models prove that dry intermediate and minor mergers can induce global structural evolution within the sequence of S0 Hubble types compatible with observations, meaning that these processes should not be discarded from the evolutionary scenarios of S0's just on the basis of the strong coupling observed between the bulge and disk scale-lengths in these galaxies (abridged).

Some recent observations seem to disagree with hierarchical theories of galaxy formation about the role played by major mergers in the late buildup of massive E-S0's. We re-address this question by analysing the morphology, structural distortion level, and star formation enhancement of a sample of massive galaxies (M_* > 5 * 10^10 Msun) lying on the Red Sequence and its surroundings at 0.3 < z < 1.5. We have used an initial sample of ~1800 sources with K_s < 20.5 mag over an area ~155 arcmin^2 on the Groth Strip, combining data from the Rainbow Extragalactic Database and the GOYA Survey. Red galaxy classes that can be directly associated to intermediate stages of major mergers and to their final products have been defined. We report observational evidence of the existence of a dominant evolutionary path among massive red galaxies at 0.6<z<1.5, consisting in the conversion of irregular disks into irregular spheroids, and of these ones into regular spheroids. This result implies that: 1) the massive red regular galaxies at low redshifts derive from the irregular ones populating the Red Sequence and its neighbourhood at earlier epochs up to z~1.5; 2) the progenitors of the bulk of present-day massive red regular galaxies have been blue disks that have migrated to the Red Sequence mostly through major mergers at 0.6 < z < 1.2 (these mergers thus starting at z~1.5); and 3) the formation of E-S0's that end up with M_*> 10^11 Msun at z=0 through gas-rich major mergers has frozen since z~0.6. All these facts support that major mergers have played the dominant role in the definitive buildup of present-day E-S0's with M_*> 10^11 Msun at 0.6<z<1.2, in good agreement with hierarchical scenarios of galaxy formation.

Recent observations show that inner discs and rings (IDs and IRs, henceforth) are not preferably found in barred galaxies, a fact that points to the relevance of formation mechanisms different to the traditional bar-origin scenario. In contrast, the role of minor mergers in the formation of these inner components (ICs), while often invoked, is still poorly understood. We have investigated the capability of minor mergers to trigger the formation of IDs and IRs in spiral galaxies through collisionless N-body simulations. We have run a battery of minor mergers in which both primary and secondary are modelled as disc-bulge-halo galaxies with realistic density ratios. A detailed analysis of the morphology, structure, and kinematics of the ICs resulting from the minor merger has been carried out. All the simulated minor mergers develop thin ICs out of satellite material, supported by rotation. A wide morphological zoo of ICs has been obtained (including IDs, IRs, pseudo-rings, nested IDs, spiral patterns, and combinations of them), but all with structural and kinematical properties similar to observations. The existence of the resulting ICs can be deduced through the features that they imprint in the isophotal profiles and kinemetric maps of the final remnant, as in many real galaxies. The realistic density ratios used in the present models make the satellites to experience more efficient orbital circularization and disruption than in previous studies. Combined with the disc resonances induced by the encounter, these processes give place to highly aligned co- and counter-rotating ICs in the remnant centre. Therefore, minor mergers are an efficient mechanism to form rotationally-supported stellar ICs in spiral galaxies, neither requiring strong dissipation nor the development of noticeable bars (abridged).

Hierarchical models predict that massive early-type galaxies (mETGs) derive from the most massive and violent merging sequences occurred in the Universe. However, the role of wet, mixed, and dry major mergers in the assembly of mETGs is questioned by some recent observations. We have developed a semi-analytical model to test the feasibility of the major-merger origin hypothesis for mETGs, just accounting for the effects on galaxy evolution of the major mergers strictly reported by observations. The model proves that it is feasible to reproduce the observed number density evolution of mETGs since z~1, just accounting for the coordinated effects of wet/mixed/dry major mergers. It can also reconcile the different assembly redshifts derived by hierarchical models and by mass downsizing data for mETGs, just considering that a mETG observed at a certain redshift is not necessarily in place since then. The model predicts that wet major mergers have controlled the mETGs buildup since z~1, although dry and mixed mergers have also played an essential role in it. The bulk of this assembly took place at 0.7<z<1, being nearly frozen at z<~0.7 due to the negligible number of major mergers occurred per existing mETG since then. The model suggests that major mergers have been the main driver for the observational migration of mass from the massive end of the blue galaxy cloud to that of the red sequence in the last ~8 Gyr.

Hierarchical models predict that massive early-type galaxies (mETGs) are the latest systems to be in place into the cosmic scenario (at z<~0.5), conflicting with the observational phenomenon of galaxy mass downsizing, which poses that the most massive galaxies have been in place earlier that their lower-mass counterparts (since z~0.7). We have developed a semi-analytical model to test the feasibility of the major-merger origin hypothesis for mETGs, just accounting for the effects on galaxy evolution of the major mergers strictly reported by observations. The most striking model prediction is that very few present-day mETGs have been really in place since z~1, because ~90% of the mETGs existing at z~1 are going to be involved in a major merger between z~1 and the present. Accounting for this, the model derives an assembly redshift for mETGs in good agreement with hierarchical expectations, reproducing observational mass downsizing trends at the same time.

Several studies have tried to ascertain whether or not the increase in abundance of the early-type galaxies (E-S0a's) with time is mainly due to major mergers, reaching opposite conclusions. We have tested it directly through semi-analytical modelling, by studying how the massive early-type galaxies with log(M_*/Msun)>11 at z~0 (mETGs) would have evolved backwards-in-time, under the hypothesis that each major merger gives place to an early-type galaxy. The study was carried out just considering the major mergers strictly reported by observations at each redshift, and assuming that gas-rich major mergers experience transitory phases of dust-reddened, star-forming galaxies (DSFs). The model is able to reproduce the observed evolution of the galaxy LFs at z<~1, simultaneously for different rest-frame bands (B, I, and K) and for different selection criteria on color and morphology. It also provides a framework in which apparently-contradictory results on the recent evolution of the luminosity function (LF) of massive, red galaxies can be reconciled, just considering that observational samples of red galaxies can be significantly contaminated by DSFs. The model proves that it is feasible to build up ~50-60% of the present-day mETG population at z<~1 and to reproduce the observational excess by a factor of ~4-5 of late-type galaxies at 0.8<z<1 through the coordinated action of wet, mixed, and dry major mergers, fulfilling global trends that are in general agreement with mass-downsizing. The bulk of this assembly takes place during ~1 Gyr elapsed at 0.8<z<1. The model suggests that major mergers have been the main driver for the observational migration of mass from the massive-end of the blue galaxy cloud to that of the red sequence in the last ~8 Gyr.(Abridged)

Hierarchical models predict that present-day massive early-type galaxies (mETGs) have finished their assembly at a quite late cosmic epoch (z~0.5), conflicting directly with galaxy mass-downsizing. In Eliche-Moral et al. (2010), we presented a semi-analytical model that predicts the increase by a factor of ~2.5 observed in the number density of mETGs since z~1 to the present, just accounting for the effects of the major mergers strictly-reported by observations. Here, we describe the relative, coordinated role of wet, mixed, and dry major mergers in driving this assembly. Accordingly to observations, the model predicts that: 1) wet major mergers have controlled the mETGs buildup since z~1, although dry and mixed mergers have also contributed significantly to it; 2) the bulk of this assembly takes place during the ~1.4 Gyr time-period elapsed at 0.7<z<1, being nearly frozen at z<~0.7; 3) this frostbite can be explained just accounting for the observational decrease of the major merger fraction since z~0.7, implying that major mergers (and, in particular, dry events) have contributed negligibly to the mETGs assembly during the last ~6.3 Gyr; and 4) major mergers are responsible for doubling the stellar mass at the massive-end of the red sequence since z~1. The most striking model prediction is that at least ~87% of the mETGs existing at z~1 are not the passively-evolved, high-z counterparts of present-day mETGs, but their gas-poor progenitors instead. This implies that <~5% of present-day mETGs have been really in place since z~1. The model derives a redshift of final assembly for present-day mETGs in agreement with hierarchical models (z~0.5), reproducing at the same time the observed buildup of mETGs at z<~1.(Abridged)

Nuclear disks and rings are frequent galaxy substructures, for a wide range of morphological types (from S0 to Sc). We have investigated the possible minor-merger origin of inner disks and rings in spiral galaxies through collisionless N-body simulations. The models confirm that minor mergers can drive the formation of thin, kinematically-cold structures in the center of galaxies out of satellite material, without requiring the previous formation of a bar. Satellite core particles tend to be deposited in circular orbits in the central potential, due to the strong circularization experienced by the satellite orbit through dynamical friction. The material of the satellite core reaches the remnant center if satellites are dense or massive, building up a thin inner disk; whereas it is fully disrupted before reaching the center in the case of low-mass satellites, creating an inner ring instead.

FRIDA (inFRared Imager and Dissector for the Adaptative optics system of the GTC) will be a NIR (1-2.5microns) imager and Integral Field Unit spectrograph to operate with the Adaptative Optics system of the 10.4m GTC telescope. FRIDA will offer broad and narrow band diffraction-limited imaging and integral field spectroscopy at low, intermediate and high spectral resolution. The Extragalactic Astrophysics and Astronomical Instrumentation group of the Universidad Complutense de Madrid (GUAIX) is developing the Data Reduction Pipeline for FRIDA. Specific tools for converting output, reduced datacubes to the standard Euro3D FITS format will be developed, in order to allow users to exploit existing VO applications for analysis. FRIDA is to be commissioned on the telescope in 2011.

(Abridged) We measure the fraction of galaxies undergoing disk-disk major mergers (f_m) at intermediate redshifts (0.35 <= z < 0.85) by studying the asymmetry index A of galaxy images. Results are provided for B- and Ks-band absolute magnitude selected samples from the Groth strip in the GOYA photometric survey. Three sources of systematic error are carefully addressed: (i) we avoid morphological K-corrections, (ii) we measure asymmetries in artificially redshifted to z_d = 0.75 galaxies to lead with loss of morphological information with redshift, and (iii) we take into account the observational errors in z and A, that tend to overestimate the merger fraction, by maximum likelihood techniques. We find: (i) our data allow for a robust merger fraction to be provided for a single redshift bin centered at z=0.6. (ii) Merger fractions have low values: f_m = 0.045 for M_B <= -20 galaxies, and f_m = 0.031 for M_Ks <= -23.5 galaxies. And, (iii) failure to address the effects of the observational errors leads to overestimating f_m by factors of 10%-60%. Combining our results with those on literature, and parameterizing the merger fraction evolution as f_m(z) = f_m(0)(1+z)^m, we obtain that m = 2.9 +- 0.8, and f_m(0) = 0.012 +- 0.004$. Assuming a Ks-band mass-to-light ratio not varying with luminosity, we infer that the merger rate of galaxies with stellar mass M >= 3.5x10^10 M_Sun is R_m = 1.6x10^-4 Mpc^-3 Gyr^-1. When we compare with previous studies at similar redshifts, we find that the merger rate decreases when mass increases.

We aim to define a sample of intermediate-z disk galaxies harbouring central bulges, and a complementary sample of disk galaxies without measurable bulges. We intend to provide colour profiles for both samples, as well as measurements of nuclear, disk, and global colours, which may be used to constrain the relative ages of bulges and disks. We select a diameter-limited sample of galaxies in images from the HST/WFPC2 Groth Strip survey, which is divided into two subsamples of higher and lower inclination to assess the role of dust in the measures quantities. Mergers are visually identified and excluded. We take special care to control the pollution by ellipticals. The bulge sample is defined with a criterion based on nuclear surface brightness excess over the inward extrapolation of the exponential law fitted to the outer regions of the galaxies. We extract colour profiles on the semi-minor axis least affected by dust in the disk, and measure nuclear colours at 0.85 kpc from the centre over those profiles. Disk colours are measured on major axis profiles; global colours are obtained from 2.6"-diameter apertures. We obtain a parent sample containing 248 galaxies with known redshifts, spectroscopic or photometric, spanning 0.1 < z < 1.2. The bulge subsample comprises 54 galaxies (21.8% of the total), while the subsample with no measureable bulges is 55.2% of the total (137 galaxies). The remainder (23%) is composed of mergers. We list nuclear, disk, and global colours (observed and restframe) and magnitudes (apparent and absolute), as well as galaxy colour gradients for the samples with and without bulges. We also provide images, colour maps, plots of spectral energy distributions, major-axis surface brightness profiles, and minor-axis colour profiles for both samples.

Satellite accretion events have been invoked for mimicking the internal secular evolutionary processes of bulge growth. However, N-body simulations of satellite accretions have paid little attention to the evolution of bulge photometric parameters, to the processes driving this evolution, and to the consistency of this evolution with observations. We want to investigate whether satellite accretions indeed drive the growth of bulges, and whether they are consistent with global scaling relations of bulges and discs. We perform N-body models of the accretion of satellites onto disc galaxies. A Tully-Fisher (M \propto V_{rot}^ {alpha_TF}) scaling between primary and satellite ensures that density ratios, critical to the outcome of the accretion, are realistic. We carry out a full structural, kinematic and dynamical analysis of the evolution of the bulge mass, bulge central concentration, and bulge-to-disc scaling relations. The remnants of the accretion have bulge-disc structure. Both the bulge-to-disc ratio (B/D) and the Sersic index (n) of the remnant bulge increase as a result of the accretion, with moderate final bulge Sersic indices: n = 1.0 to 1.9. Bulge growth occurs no matter the fate of the secondary, which fully disrupts for alpha_TF=3 and partially survives to the remnant center for alpha_TF = 3.5 or 4. Global structural parameters evolve following trends similar to observations. We show that the dominant mechanism for bulge growth is the inward flow of material from the disc to the bulge region during the satellite decay. The models confirm that the growth of the bulge out of disc material, a central ingredient of secular evolution models, may be triggered externally through satellite accretion.

We present U and B galaxy differential number counts from a field of ~900 arcmin^2, based on GOYA Survey imaging of the HST Groth-Westphal strip. Source detection efficiency corrections as a function of the object size have been applied. A variation of the half-exposure image method has been devised to identify and remove spurious detections. Achieved 50% detection efficiencies are 24.8 mag in U and 25.5 mag in B in the Vega system. Number count slopes are d log (N)/dm = 0.50 +/- 0.02 for B=21.0-24.5, and d log (N)/dm = 0.48 +/- 0.03 for U=21.0-24.0. Simple number count models are presented that simultaneously reproduce the counts over 15 mag in U and B,and over 10 mag in K_s, using a Lambda-dominated cosmology and SDSS local luminosity functions. Only by setting a recent z_f ~ 1.5 formation redshift for early-type, red galaxies do the models reproduce the change of slope observed at Ks = 17.5 in NIR counts. A moderate optical depth (tau_B = 0.6) for all galaxy types ensures that the recent formation for ellipticals does not leave a signature in the U or B number counts, which are featureless at intermediate magnitudes. No ad-hoc disappearing populations are needed to explain the counts if number evolution is introduced using an observationally-based z-evolution of the merger fraction.

For bulges of spiral galaxies, the concentration, or Sersic index, increases with bulge luminosity and bulge-to-disk ratio B/D (Andredakis, Peletier, & Balcells 1995, MNRAS, 275, 874). Does this trend trace the growth of bulges via satellite accretion? And, is satellite infall consistent with this trend? Aguerri, Balcells, & Peletier (2001, A&A, 367, 428) investigated this question with N-body simulations of the accretion of dense, spheroidal satellites. Here, we expand on that work by running N-body simulations of the accretion of satellites that have realistic densities. Satellites are modelled as disk-bulge structures with their own dark-matter halo. A realistic density scaling with the primary galaxy is ensured by using the Tully-Fisher relation. Our merger models show that most satellites disrupt before reaching the center. However, a bulge-disk decomposition of the surface density profile after the accretion shows an increase of both the B/D and the Sersic index n of the bulge. The increase in the mass and concentration of the inner Sersic component is due to inward piling up of disk material due to transient bars during the satellite orbital decay. This research is described in Eliche-Moral, Balcells, Aguerri, & Gonzalez-Garcia, 2005 (in preparation).