Rubens E. G. Machado

Bar formation and evolution in disc galaxies with gas and a triaxial halo: morphology, bar strength and halo properties

We follow the formation and evolution of bars in N-body simulations of disc galaxies with gas and/or a triaxial halo. We find that both the relative gas fraction and the halo shape play a major role in the formation and evolution of the bar. In gas-rich simulations, the disc stays near-axisymmetric much longer than in gas-poor ones, and, when the bar starts growing, it does so at a much slower rate. Due to these two effects combined, large-scale bars form much later in gas-rich than in gas-poor discs. This can explain the observation that bars are in place earlier in massive red disc galaxies than in blue spirals. We also find that the morphological characteristics in the bar region are strongly influenced by the gas fraction. In particular, the bar at the end of the simulation is much weaker in gas-rich cases. In no case did we witness bar destruction. Halo triaxiality has a dual influence on bar strength. In the very early stages of the simulation it induces bar formation to start earlier. On the other hand, during the later, secular evolution phase, triaxial haloes lead to considerably less increase of the bar strength than spherical ones. The shape of the halo evolves considerably with time. The inner halo parts may become more elongated, or more spherical, depending on the bar strength. The main body of initially triaxial haloes evolves towards sphericity, but in initially strongly triaxial cases it stops well short of becoming spherical. Part of the angular momentum absorbed by the halo generates considerable rotation of the halo particles that stay located relatively near the disc for long periods of time. Another part generates halo bulk rotation, which, contrary to that of the bar, increases with time but stays small.

Loss of halo triaxiality due to bar formation

Cosmological N-body simulations indicate that the dark matter haloes of galaxies should be generally triaxial. Yet, the presence of a baryonic disc is believed to alter the shape of the haloes. Here we aim to study how bar formation is affected by halo triaxiality and how, in turn, the presence of the bar influences the shape of the halo. We perform a set of collisionless N-body simulations of disc galaxies with triaxial dark matter haloes, using elliptical discs as initial conditions. Such discs are much closer to equilibrium with their haloes than circular ones, and the ellipticity of the initial disc depends on the ellipticity of the halo gravitational potential. For comparison, we also consider models with initially circular discs, and find that the differences are very important. In all cases, the mass of the disc is grown quasi-adiabatically within the haloes, but the time-scale of growth is not very important. We study models of different halo triaxialities and, to investigate the behaviour of the halo shape in the absence of bar formation, we run models with different disc masses, halo concentrations, disc velocity dispersions and also models where the disc shape is kept artificially axisymmetric. We find that the introduction of a massive disc, even if this is not circular, causes the halo triaxiality to be partially diluted. Once the disc is fully grown, a strong stellar bar develops within the halo that is still non-axisymmetric, causing it to lose its remaining non-axisymmetry. In triaxial haloes in which the parameters of the initial conditions are such that a bar does not form, the halo is able to remain triaxial and the circularization of its shape on the plane of the disc is limited to the period of disc growth. We conclude that part of the circularization of the halo is due to disc growth, but part must be attributed to the formation of a bar. Bars in the halo component, which have already been found in axisymmetric haloes, are also found in triaxial ones. We find that initially circular discs respond excessively to the triaxial potential and become highly elongated. They also lose more angular momentum than the initially elliptical discs and thus form stronger bars. Because of that, the circularization that their bars induce on their haloes is also more rapid. We also analyse halo vertical shapes and observe that their vertical flattenings remain considerable, meaning that the haloes become approximately oblate by the end of the simulations. Finally, we also analyse the kinematics of a subset of halo particles that rotate in disc-like manner. These particles occupy a layer around the plane of the disc and their rotation is more important in the spherical halo than in triaxial ones. We also find that, even though the final shape of the halo is roughly independent of the initial shape, the initially triaxial ones are able to retain the anisotropy of their velocity dispersions.

Chaos and dynamical trends in barred galaxies: bridging the gap between N-body simulations and time-dependent analytical models

Self-consistent N-body simulations are efficient tools to study galactic dynamics. However, using them to study individual trajectories (or ensembles) in detail can be challenging. Such orbital studies are important to shed light on global phase space properties, which are the underlying cause of observed structures. The potentials needed to describe self-consistent models are time-dependent. Here, we aim to investigate dynamical properties (regular/chaotic motion) of a non-autonomous galactic system, whose time-dependent potential adequately mimics certain realistic trends arising from N-body barred galaxy simulations. We construct a fully time-dependent analytical potential, modeling the gravitational potentials of disc, bar and dark matter halo, whose time-dependent parameters are derived from a simulation. We study the dynamical stability of its reduced time-independent 2-degrees of freedom model, charting the different islands of stability associated with certain orbital morphologies and detecting the chaotic and regular regions. In the full 3-degrees of freedom time-dependent case, we show representative trajectories experiencing typical dynamical behaviours, i.e., interplay between regular and chaotic motion for different epochs. Finally, we study its underlying global dynamical transitions, estimating fractions of (un)stable motion of an ensemble of initial conditions taken from the simulation. For such an ensemble, the fraction of regular motion increases with time.

Simulations of the merging galaxy cluster Abell 3376

Observed galaxy clusters often exhibit X-ray morphologies suggestive of recent interaction with an infalling subcluster. Abell 3376 is a nearby (z=0.046) massive galaxy cluster whose bullet-shaped X-ray emission indicates that it may have undergone a recent collision. It displays a pair of Mpc-scale radio relics and its brightest cluster galaxy is located 970 h_70^-1 kpc away from the peak of X-ray emission, where the second brightest galaxy lies. We attempt to recover the dynamical history of Abell 3376. We perform a set of N-body adiabatic hydrodynamical simulations using the SPH code Gadget-2. These simulations of binary cluster collisions are aimed at exploring the parameter space of possible initial configurations. By attempting to match X-ray morphology, temperature, virial mass and X-ray luminosity, we set approximate constraints on some merger parameters. Our best models suggest a collision of clusters with mass ratio in the range 1/6-1/8, and having a subcluster with central gas density four times higher than that of the major cluster. Models with small impact parameter (b<150 kpc), if any, are preferred. We estimate that Abell 3376 is observed approximately 0.5 Gyr after core passage, and that the collision axis is inclined by i~40 degrees with respect to the plane of the sky. The infalling subcluster drives a supersonic shock wave that propagates at almost 2600 km/s, implying a Mach number as high as M~4; but we show how it would have been underestimated as M~3 due to projection effects.

The merger history of the complex cluster Abell 1758: a combined weak lensing and spectroscopic view

We present a weak-lensing and dynamical study of the complex cluster Abell 1758 (A1758, z = 0.278) supported by hydrodynamical simulations. This cluster is composed of two main structures, called A1758N and A1758S. The Northern structure is composed of A1758NW & A1758NE, with lensing determined masses of 7.90_{-1.55}^{+1.89} X 10^{14} M_\odot and 5.49_{-1.33}^{+1.67} X 10^{14} M_\odot, respectively. They show a remarkable feature: while in A1758NW there is a spatial agreement among weak lensing mass distribution, intracluster medium and its brightest cluster galaxy (BCG) in A1758NE the X-ray peak is located 96_{-15}^{+14} arcsec away from the mass peak and BCG positions. Given the detachment between gas and mass we could use the local surface mass density to estimate an upper limit for the dark matter self-interaction cross section: \sigma/m<5.83 cm^2 g^{-1}. Combining our velocity data with hydrodynamical simulations we have shown that A1758 NW \& NE had their closest approach 0.27 Gyr ago and their merger axis is 21+-12 degrees from the plane of the sky. In the A1758S system we have measured a total mass of 4.96_{-1.19}^{+1.08} X 10^{14} M_\odot and, using radial velocity data, we found that the main merger axis is located at 70+-4 degrees from the plane of the sky, therefore closest to the line-of-sight.

Weak lensing and spectroscopic analysis of the nearby dissociative merging galaxy cluster Abell 3376

The galaxy cluster Abell~3376 is a nearby (z=0.046) dissociative merging cluster surrounded by two prominent radio relics and showing an X-ray comet-like morphology. The merger system is comprised of the subclusters A3376W & A3376E. Based on new deep multi-wavelength large-field images and published redshifts, we bring new insights about the history of this merger. Despite the difficulty of applying the weak lensing technique at such low redshift, we successfully recovered the mass distribution in the cluster field. Moreover, with the application of a two-body model, we have addressed the dynamics of these merging system. We have found the individual masses of M_{200}^{W}=3.0_{-1.7}^{+1.3}x10^{14} M_{\odot} and M_{200}^{E}=0.9_{-0.8}^{+0.5}x10^{14} M_{\odot}. The cometary shaped X-ray distribution shows only one peak spatially coincident with both Eastern BCG and the A3376E mass peak whereas the gas content of A3376W seems to be stripped out. Our data allowed us to confirm the existence of a third subcluster located at the North, 1147+-62 kpc apart from the neighbour subcluster A3376E and having a mass M_{200}^{N}=1.4_{-1.0}^{+0.7}x10^{14} M_{\odot}. From our dynamical analysis, we found the merging is taking place very close to the plane of the sky, with the merger axis just 10 deg +-11 deg from it. The application of a two-body analysis code showed that the merging cluster is seen 0.9_{-0.3}^{+0.2} Gyr after the pericentric passage and it is currently going to the point of maximum separation between the subclusters.

Simulations of gas sloshing in galaxy cluster Abell 2052

The intracluster plasma of Abell 2052 exhibits in X-rays a spiral structure extending more than 250 kpc, which is comprised of cool gas. This feature is understood to be the result of gas sloshing caused by the off-axis collision with a smaller subcluster. We aim to recover the dynamical history of Abell 2052 and to reproduce the broad morphology of the spiral feature. To this end, we perform hydrodynamical

Chaotic motion and the evolution of morphological components in a time-dependent model of a barred galaxy within a dark matter halo

N

Uma introdução ao barroquismo de Glauber Rocha: o espaço ambíguo de Terra em transe

-body simulations of cluster collisions. We obtain two regimes that adequately reproduce the desired features. The first scenario is a close encounter and a recent event (0.8 Gyr since pericentric passage), while the second scenario has a larger impact parameter and is older (almost 2.6 Gyr since pericentric passage). In the second case, the simulation predicts that the perturbing subcluster should be located approximately 2 Mpc from the centre of the major cluster. At that position, we are able to identify an observed optical counterpart at the same redshift: a galaxy group with

Simulating the shocks in the dissociative galaxy cluster Abell 1758N

M_{500} = (1.16 \pm 0.43) \times 10^{13} M_{\odot}