Horizon-AGN Simulation Research Articles

How complex are galaxies? A non-parametric estimation of the intrinsic dimensionality of wide-band photometric data

Abstract Galaxies are complex objects, yet the number of independent parameters to describe them remains unknown. We present here a non-parametric method to estimate the intrinsic dimensionality of large datasets. We apply it to wide-band photometric data drawn from the COSMOS2020 catalogue and a comparable mock catalogue from the Horizon-AGN simulation. Our galaxy catalogues are limited in signal-to-noise ratio in all optical and NIR bands. Our results reveal that most of the variance in the wide-band photometry of this galaxy sample can be described with at most 4.3 ± 0.5 independent parameters for star-forming galaxies and 2.9 ± 0.2 for passive ones, both in the observed and simulated catalogues. We identify one of these parameters to be noise-driven, and recover that stellar mass and redshift are two key independent parameters driving the magnitudes. Our findings support the idea that wide-band photometry does not provide more than one additional independent parameter for star-forming galaxies. Although our sample is not mass-limited and may miss some passive galaxies due to our cut in SNR, our work suggests that dimensionality reduction techniques may be effectively used to explore and analyse wide-band photometric data, provided the used latent space is at least four-dimensional.

Assembly of the intracluster light in the Horizon-AGN simulation

ABSTRACT The diffuse stellar component of galaxy clusters made up of intergalactic stars is termed the intracluster light (ICL). Although there is a developing understanding of the mechanisms by which the ICL is formed, no strong consensus has yet been reached on which objects the stars of the ICL are primarily sourced from. We investigate the assembly of the ICL starting approximately 10 Gyr before $z=0$ in 11 galaxy clusters (halo masses between $\sim 1\times 10^{14}$ and $\sim 7\times 10^{14}$ M$_{\odot }$ at $z\approx 0$) in the Horizon-AGN simulation. By tracking the stars of galaxies that fall into these clusters past cluster infall, we are able to link almost all of the $z\approx 0$ ICL back to progenitor objects. Satellite stripping, mergers, and pre-processing are all found to make significant contributions to the ICL, but any contribution from in situ star formation directly into the ICL appears negligible. Even after compensating for resolution effects, we find that approximately 90 per cent of the stacked ICL of the 11 clusters that is not pre-processed should come from galaxies infalling with stellar masses above $10^{9}$ M$_{\odot }$, with roughly half coming from infalling galaxies with stellar masses within half a dex of $10^{11}$ M$_{\odot }$. The fact that the ICL appears largely sourced from such massive objects suggests that the ICL assembly of any individual cluster may be principally stochastic.

BayeSED-GALAXIES. I. Performance Test for Simultaneous Photometric Redshift and Stellar Population Parameter Estimation of Galaxies in the CSST Wide-field Multiband Imaging Survey

The forthcoming Chinese Space Station Telescope (CSST) wide-field multiband imaging survey will produce seven-band photometric spectral energy distributions (SEDs) for billions of galaxies. The effective extraction of astronomical information from these massive data sets of SEDs relies on the techniques of SED synthesis (or modeling) and SED analysis (or fitting). We evaluate the performance of the latest version of the BayeSED code combined with SED models with increasing complexity for simultaneously determining the photometric redshifts and stellar population parameters of galaxies in this survey. By using an empirical statistics–based mock galaxy sample without SED modeling errors, we show that the random observational errors in photometries are more important sources of errors than the parameter degeneracies and Bayesian analysis method and tool. By using a Horizon-AGN hydrodynamical simulation–based mock galaxy sample with SED modeling errors about the star formation histories (SFHs) and dust attenuation laws (DALs), the simple typical assumptions lead to significantly worse parameter estimation with CSST photometries only. SED models with more flexible (or complicated) forms of SFH/DAL do not necessarily lead to better estimation of redshift and stellar population parameters. We discuss the selection of the best SED model by means of Bayesian model comparison in different surveys. Our results reveal that Bayesian model comparison with Bayesian evidence may favor SED models with different complexities when using photometries from different surveys. Meanwhile, the SED model with the largest Bayesian evidence tends to give the best performance of parameter estimation, which is clearer for photometries with higher discriminative power.

Alignments between Galaxies and the Cosmic Web at z ∼ 1–2 in the IllustrisTNG Simulations

Galaxy formation theories predict that galaxy shapes and angular momenta have nonrandom alignments with the cosmic web. This leads to so-called intrinsic alignment between pairs of galaxies, which is important to quantify as a nuisance parameter for weak lensing. We study galaxy–cosmic web alignment in the IllustrisTNG suite of hydrodynamical simulations at redshifts 1 and 2, finding that alignment trends are consistent with previous studies. However, we find that the magnitude of the spin alignment signal is ∼2.4× weaker than seen in previous studies of the Horizon-AGN simulation, suggesting that this signal may have a significant dependence on subgrid physics. Based on IllustrisTNG, we then construct mock observational spectroscopic surveys that can probe shape–cosmic web alignment at z ∼ 1–2, modeled on the low-z galaxy redshift and IGM tomography surveys on the upcoming Subaru Prime Focus Spectrograph Galaxy Evolution (PFS GE) survey. However, even over box sizes of L = 205 h −1 Mpc, we find that global anisotropies induce a sample variance in the 2D projected alignment signal that depend on the projected direction; this induces significant errors in the observed alignment. We predict a 5.3σ detection of IllustrisTNG’s shape alignment signal at z ∼ 1 from Subaru PFS GE, although a detection would be challenging at z ∼ 2. However, a rough rescaling of the relative alignment signal strengths between the TNG and Horizon-AGN simulations suggests that PFS GE should be able to more easily constrain the latter’s stronger signal.

Evidence for non-merger co-evolution of galaxies and their supermassive black holes

ABSTRACT Recent observational and theoretical studies have suggested that supermassive black holes (SMBHs) grow mostly through non-merger (‘secular’) processes. Since galaxy mergers lead to dynamical bulge growth, the only way to observationally isolate non-merger growth is to study galaxies with low bulge-to-total mass ratio (e.g. $B/T\lt 10~{{\ \rm per\ cent}}$). However, bulge growth can also occur due to secular processes, such as disc instabilities, making disc-dominated selections a somewhat incomplete way to select merger-free systems. Here we use the Horizon-AGN simulation to select simulated galaxies which have not undergone a merger since z = 2, regardless of bulge mass, and investigate their location on typical black hole-galaxy scaling relations in comparison to galaxies with merger dominated histories. While the existence of these correlations has long been interpreted as co-evolution of galaxies and their SMBHs driven by galaxy mergers, we show here that they persist even in the absence of mergers. We find that the correlations between SMBH mass and both total mass and stellar velocity dispersion are independent of B/T ratio for both merger-free and merger-dominated galaxies. In addition, the bulge mass and SMBH mass correlation is still apparent for merger-free galaxies, the intercept for which is dependent on B/T. Galaxy mergers reduce the scatter around the scaling relations, with merger-free systems showing broader scatter. We show that for merger-free galaxies, the co-evolution is dominated by radio-mode feedback, and suggest that the long periods of time between galaxy mergers make an important contribution to the co-evolution between galaxies and SMBHs in all galaxies.

Supermassive black holes in merger-free galaxies have higher spins which are preferentially aligned with their host galaxy

ABSTRACT Here, we use the Horizon–active galactic nucleus (AGN) simulation to test whether the spins of supermassive black hole (SMBH) in merger-free galaxies are higher. We select samples using an observationally motivated bulge-to-total mass ratio of <0.1, along with two simulation-motivated thresholds selecting galaxies which have not undergone a galaxy merger since z = 2, and those SMBHs with $\lt 10~{{\ \rm per\ cent}}$ of their mass due to SMBH mergers. We find higher spins (>5σ) in all three sample compared to the rest of the population. In addition, we find that SMBHs with their growth dominated by BH mergers following galaxy mergers are less likely to be aligned with their galaxy spin than those that have grown through accretion in the absence of galaxy mergers (3.4σ). We discuss the implications this has for the impact of active galactic nucleus (AGN) feedback, finding that merger-free SMBHs spend on average 91 per cent of their lifetimes since z = 2 in a radio mode of feedback (88 per cent for merger-dominated galaxies). Given that previous observational and theoretical works have concluded that merger-free processes dominate SMBH-galaxy co-evolution, our results suggest that this co-evolution could be regulated by radio mode AGN feedback.

Star Formation History and Transition Epoch of Cluster Galaxies Based on the Horizon-AGN Simulation

Cluster galaxies exhibit substantially lower star formation rates than field galaxies today, but it is conceivable that clusters were sites of more active star formation in the early universe. Herein, we present an interpretation of the star formation history (SFH) of group/cluster galaxies based on the large-scale cosmological hydrodynamic simulation, Horizon-AGN. We find that massive galaxies in general have small values of e-folding timescales of star formation decay (i.e., “mass quenching”) regardless of their environment, while low-mass galaxies exhibit prominent environmental dependence. In massive host halos (i.e., clusters), the e-folding timescales of low-mass galaxies are further decreased if they reside in such halos for a longer period of time. This “environmental quenching” trend is consistent with the theoretical expectation from ram pressure stripping. Furthermore, we define a “transition epoch” as where cluster galaxies become less star-forming than field galaxies. The transition epoch of group/cluster galaxies varies according to their stellar and host-cluster halo masses. Low-mass galaxies in massive clusters show the earliest transition epoch of ∼7.6 Gyr ago in lookback time. However, this decreases to ∼5.2 Gyr for massive galaxies in low-mass clusters. Based on our findings, we can describe a cluster galaxy’s SFH with regard to the cluster halo-to-stellar mass ratio.

Dual AGN in the Horizon-AGN simulation and their link to galaxy and massive black hole mergers, with an excursus on multiple AGN

ABSTRACT The occurrence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on massive black holes (MBHs) and to derive clues on MBH mergers, using dual AGN as a parent population of precursors. We investigate the properties of dual AGN in the cosmological simulation Horizon-AGN. We create catalogs of dual AGN selected with distance and luminosity criteria, plus sub-catalogs where further mass cuts are applied. We divide the sample into dual AGN hosted in different galaxies, on the way to a merger, and into those hosted in one galaxy, after the galaxy merger has happened. We find that the relation between MBH and galaxy mass is similar to that of general AGN population and we compare the properties of dual AGN also with a control sample, discussing differences and similarities in masses and Eddington ratios. The typical mass ratio of galaxy mergers associated to dual AGN is 0.2, with mass loss in the smaller galaxy decreasing the mass ratio as the merger progresses. Between 30 and 80 per cent of dual AGN with separations between 4 and 30 kpc can be matched to an ensuing MBH merger. The dual AGN fraction increases with redshift and with separation threshold, although above 50 kpc the increase of multiple AGN limits that of duals. Multiple AGN are generally associated with massive haloes, and mass loss of satellites shapes the galaxy–halo relation.

Optical Spectroscopy of Dual Quasar Candidates from the Subaru HSC-SSP program

Abstract We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars, which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest-redshift (z = 3.1) quasar pair with a separation <10 kpc. Based on the sample acquired to date, the success rate of identifying physically associated dual quasars is 19% when excluding stars based on their HSC colors. Using the full sample of six spectroscopically confirmed dual quasars, including three previously published, we find that the black holes in these systems have black hole masses (M BH ∼ 108−9 M ⊙), bolometric luminosities (log L bol ∼ 44.5–47.5 erg s–1) and Eddington ratios (0.01–0.3) similar to single SDSS quasars. We measure the stellar mass of their host galaxies based on 2D image decomposition of the five-band (grizy) optical emission and assess the mass relation between supermassive black holes (SMBHs) and their hosts. Dual SMBHs appear to have elevated masses relative to their host galaxies. Thus, mergers may not necessarily align such systems onto the local mass relation, as suggested by the Horizon-AGN simulation. This study suggests that dual luminous quasars are triggered by mergers prior to the final coalescence of the two SMBHs, resulting in early mass growth of the black holes relative to their host galaxies.

Calibrating galaxy formation effects in galactic tests of fundamental physics

Galactic scale tests have proven to be powerful tools in constraining fundamental physics in previously under-explored regions of parameter space. The astrophysical regime which they probe is inherently complicated, and the inference methods used to make these constraints should be robust to baryonic effects. Previous analyses have assumed simple empirical models for astrophysical noise without detailed calibration or justification. We outline a framework for assessing the reliability of such methods by constructing and testing more advanced baryonic models using cosmological hydrodynamical simulations. As a case study, we use the Horizon-AGN simulation to investigate warping of stellar disks and offsets between gas and stars within galaxies, which are powerful probes of screened fifth forces. We show that the degree of `U'-shaped warping of galaxies is well modelled by Gaussian random noise, but that the magnitude of the gas-star offset is correlated with the virial radius of the host halo. By incorporating this correlation we confirm recent results ruling out astrophysically relevant Hu-Sawicki $f(R)$ gravity, and identify a $\sim 30\%$ systematic uncertainty due to baryonic physics. Such an analysis must be performed case-by-case for future galactic tests of fundamental physics.

Star–Gas Misalignment in Galaxies. II. Origins Found from the Horizon-AGN Simulation

Abstract There have been many studies aiming to reveal the origins of the star–gas misalignment found in galaxies, but there still is a lack of understanding of the contribution from each formation channel candidate. We aim to answer the question by investigating the misaligned galaxies in the Horizon-AGN simulation. There are 27,903 galaxies of stellar mass M * > 1010 M ⊙ in our sample, of which 5984 are in a group in the halo mass of M 200 > 1012 M ⊙. We have identified four main formation channels of misalignment and quantified their levels of contribution: mergers (35%), interaction with nearby galaxies (23%), interaction with dense environments or their central galaxies (21%), and secular evolution, including smooth accretion from neighboring filaments (21%). We found in the simulation that the gas, rather than stars, is typically more vulnerable to dynamical disturbances; hence, misalignment formation is mainly due to the change in the rotational axis of the gas rather than stars, regardless of the origin. We have also inspected the lifetime (duration) of the misalignment. The decay timescale of the misalignment shows a strong anticorrelation with the kinematic morphology (V/σ) and the cold gas fraction of the galaxy. The misalignment has a longer lifetime in denser regions, which is linked with the environmental impact on the host galaxy. There is a substantial difference in the length of the misalignment lifetime depending on the origin, and it can be explained by the magnitude of the initial position angle offset and the physical properties of the galaxies.

Euclid preparation

The analysis of weak gravitational lensing in wide-field imaging surveys is considered to be a major cosmological probe of dark energy. Our capacity to constrain the dark energy equation of state relies on an accurate knowledge of the galaxy mean redshift ⟨z⟩. We investigate the possibility of measuring ⟨z⟩ with an accuracy better than 0.002 (1 + z) in ten tomographic bins spanning the redshift interval 0.2 < z < 2.2, the requirements for the cosmic shear analysis of Euclid. We implement a sufficiently realistic simulation in order to understand the advantages and complementarity, as well as the shortcomings, of two standard approaches: the direct calibration of ⟨z⟩ with a dedicated spectroscopic sample and the combination of the photometric redshift probability distribution functions (zPDFs) of individual galaxies. We base our study on the Horizon-AGN hydrodynamical simulation, which we analyse with a standard galaxy spectral energy distribution template-fitting code. Such a procedure produces photometric redshifts with realistic biases, precisions, and failure rates. We find that the current Euclid design for direct calibration is sufficiently robust to reach the requirement on the mean redshift, provided that the purity level of the spectroscopic sample is maintained at an extremely high level of > 99.8%. The zPDF approach can also be successful if the zPDF is de-biased using a spectroscopic training sample. This approach requires deep imaging data but is weakly sensitive to spectroscopic redshift failures in the training sample. We improve the de-biasing method and confirm our finding by applying it to real-world weak-lensing datasets (COSMOS and KiDS+VIKING-450).

Beyond halo mass: quenching galaxy mass assembly at the edge of filaments

ABSTRACT We examine how the mass assembly of central galaxies depends on their location in the cosmic web. The Horizon-AGN simulation is analysed at z ∼ 2 using the DisPerSE code to extract multi-scale cosmic filaments. We find that the dependency of galaxy properties on large-scale environment is mostly inherited from the (large-scale) environmental dependency of their host halo mass. When adopting a residual analysis that removes the host halo mass effect, we detect a direct and non-negligible influence of cosmic filaments. Proximity to filaments enhances the build-up of stellar mass, a result in agreement with previous studies. However, our multi-scale analysis also reveals that, at the edge of filaments, star formation is suppressed. In addition, we find clues for compaction of the stellar distribution at close proximity to filaments. We suggest that gas transfer from the outside to the inside of the haloes (where galaxies reside) becomes less efficient closer to filaments, due to high angular momentum supply at the vorticity-rich edge of filaments. This quenching mechanism may partly explain the larger fraction of passive galaxies in filaments, as inferred from observations at lower redshifts.

Spatially offset black holes in the Horizon-AGN simulation and comparison to observations

ABSTRACT We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a subgrid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that (i) the fraction of spatially offset BHs increases with cosmic time, (ii) BHs live on prograde orbits in the plane of the galaxy with an orbital radius that decays with time but stalls near z = 0, and (iii) the magnitudes of offsets from the galaxy centres are substantially larger in the simulation than in observations. We attribute the stalling of the infall and excessive offset magnitudes to the fact that dynamical friction from stars and dark matter is not modelled in the simulation, and hence provide a way to improve the BH dynamics of future simulations.

Dual Supermassive Black Holes at Close Separation Revealed by the Hyper Suprime-Cam Subaru Strategic Program

Abstract The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known Sloan Digital Sky Survey quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.″6–4.″0 correspond to projected separations of 3–30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck I/Low Resolution Imaging Spectrometer and Gemini-N/Near-Infrared Integral Field Spectrometer, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [O iii]λ5007 emission lines have mild velocity offsets, thus the joint [O iii] line profile is not double-peaked. We find a dual quasar fraction of 0.26 ± 0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

The impact of AGN feedback on the 1D power spectra from the Ly α forest using the Horizon-AGN suite of simulations

ABSTRACT The Lyman-α forest is a powerful probe for cosmology, but it is also strongly impacted by galaxy evolution and baryonic processes such as active galactic nucleus (AGN) feedback, which can redistribute mass and energy on large scales. We constrain the signatures of AGN feedback on the 1D power spectrum of the Lyman-α forest using a series of eight hydro-cosmological simulations performed with the adaptive mesh refinement code ramses. This series starts from the Horizon-AGN simulation and varies the subgrid parameters for AGN feeding, feedback, and stochasticity. These simulations cover the whole plausible range of feedback and feeding parameters according to the resulting galaxy properties. AGNs globally suppress the Lyman-α power at all scales. On large scales, the energy injection and ionization dominate over the supply of gas mass from AGN-driven galactic winds, thus suppressing power. On small scales, faster cooling of denser gas mitigates the suppression. This effect increases with decreasing redshift. We provide lower and upper limits of this signature at nine redshifts between z = 4.25 and 2.0, making it possible to account for it at post-processing stage in future work given that running simulations without AGN feedback can save considerable amounts of computing resources. Ignoring AGN feedback in cosmological inference analyses leads to strong biases with 2 per cent shift on σ8 and 1 per cent shift on ns, which represents twice the standards deviation of the current constraints on ns.

Star–Gas Misalignment in Galaxies. I. The Properties of Galaxies from the Horizon-AGN Simulation and Comparisons to SAMI

Abstract Recent integral field spectroscopy observations have found that about 11% of galaxies show star–gas misalignment. The misalignment possibly results from external effects such as gas accretion, interaction with other objects, and other environmental effects, hence providing clues to these effects. We explore the properties of misaligned galaxies using Horizon-AGN, a large-volume cosmological simulation, and compare the results with those of the Sydney-AAO Multi-object integral field spectrograph (SAMI) Galaxy Survey. Horizon-AGN can match the overall misalignment fraction and reproduces the distribution of misalignment angles found by observations surprisingly closely. The misalignment fraction is found to be highly correlated with galaxy morphology both in observations and in the simulation: early-type galaxies are substantially more frequently misaligned than late-type galaxies. The gas fraction is another important factor associated with misalignment in the sense that misalignment increases with decreasing gas fraction. However, there is a significant discrepancy between the SAMI and Horizon-AGN data in the misalignment fraction for the galaxies in dense (cluster) environments. We discuss possible origins of misalignment and disagreement.

Kinematic unrest of low mass galaxy groups

In an effort to better understand the formation of galaxy groups, we examine the kinematics of a large sample of spectroscopically confirmed X-ray galaxy groups in the Cosmic Evolution Survey with a high sampling of galaxy group members up to z = 1. We compare our results with predictions from the cosmological hydrodynamical simulation of HORIZON-AGN. Using a phase-space analysis of dynamics of groups with halo masses of M200c ∼ 1012.6 − 1014.50 M⊙, we show that the brightest group galaxies (BGG) in low mass galaxy groups (M200c < 2 × 1013 M⊙) have larger proper motions relative to the group velocity dispersion than high mass groups. The dispersion in the ratio of the BGG proper velocity to the velocity dispersion of the group, σBGG/σgroup, is on average 1.48 ± 0.13 for low mass groups and 1.01 ± 0.09 for high mass groups. A comparative analysis of the HORIZON-AGN simulation reveals a similar increase in the spread of peculiar velocities of BGGs with decreasing group mass, though consistency in the amplitude, shape, and mode of the BGG peculiar velocity distribution is only achieved for high mass groups. The groups hosting a BGG with a large peculiar velocity are more likely to be offset from the Lx − σv relation; this is probably because the peculiar motion of the BGG is influenced by the accretion of new members.

Real galaxy mergers from galaxy pair catalogues

ABSTRACT Mergers of galaxies are extremely violent events shaping their evolution. Such events are thought to trigger starbursts and, possibly, black hole accretion. Nonetheless, it is still not clear how to know the fate of a galaxy pair from the data available at a given time, limiting our ability to constrain the exact role of mergers. In this paper we use the light-cone of the horizon-agn simulation, for which we know the fate of each pair, to test three selection processes aiming at identifying true merging pairs. We find that the simplest one (selecting objects within two thresholds on projected distance d and redshift difference Δz) gives similar results than the most complex one (based on a neural network analysing d, Δz, redshift of the primary, masses/star formation rates/aspect ratio of both galaxies). Our best thresholds are $d_\mathrm{th}\sim 100\, \mathrm{kpc}$ and Δzth ∼ 10−3, in agreement with recent results.

Cosmological evolution of orientations of cluster-sized dark matter haloes and their central galaxies in the Horizon-AGN simulation

Abstract It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than 5 × 1013 M⊙ and their central galaxies (CGs) at z ≈ 0 in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between 0 < z < 5. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, z > 1. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within ∼30○ in the three dimensional space and ∼20○ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the slowest collapsing (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch.