Source author record

Hugo Martel

Hugo Martel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2026arXiv

Integral field spectroscopy and numerical simulations of the NGC 2207/IC 2163 system

We present integral field spectroscopy of the interacting galaxy system NGC 2207/IC 2163 obtained with the imaging Fourier Transform Spectrometer SITELLE. Approximately 1000 HII region complexes are detected in both galaxies and analyzed using their strong optical emission lines. Their properties were studied via BPT diagrams and their luminosity function. We conducted a detailed study of the distribution of oxygen abundance across the system using a series of strong-line O/H indicators and calibrations. Both galaxies exhibit negative galactocentric abundance gradients with a slope -0.015 dex kpc$^{-1}$. There are marginal signs of discontinuities in the O/H gradients with some indicators while no significant azimuthal variations are seen. A shallower slope in the HII region luminosity function between the arm and inter-arm regions in IC 2163 is observed, supporting previous conclusion that the star formation process in this galaxy eyelids has been altered during the interaction. The kinematics of the ionised gas reveal disturbed velocity fields, AGN-like features in the nucleus of NGC 2207, and elevated velocity dispersion in turbulent or feedback-driven regions. To interpret these findings, we modeled the collision using the numerical algorithm GCD+. The simulation reproduces key features of the system and demonstrates how close passages drive enhanced star formation and localized chemical enrichment. Finally, two dwarf galaxies in the field are found to have very similar systemic velocities as their larger counterparts, and could well play a minor role in the global interaction based on their morphology and position.

preprint2022arXiv

The interacting pair of galaxies Arp 82: Integral field spectroscopy and numerical simulations

Spectral data cubes of the interacting pair of galaxies NGC 2535 and NGC 2536 (the Arp 82 system) targeting bright emission lines in the visible band, obtained with the imaging Fourier transform spectrometer (iFTS) SITELLE attached to the Canada-France-Hawaii Telescope (CFHT), are presented. Analysis of H$\upalpha$ velocity maps reveals a bar in $\rm NGC\,2536$. In $\rm NGC\,2535$, we find strong non-circular motions outside the ocular ring, in the elliptical arc and tidal tails of $\rm NGC\,2535$ and a misalignment between the kinematic and photometric position angles. We detect 155 HII region complexes in the interacting pair of galaxies and determine oxygen abundances for 66 of them using different calibrators. We find, regardless of the indicator used, that the oxygen abundance distribution in $\rm NGC\,2536$ is shallow whereas, in $\rm NGC\,2535$, it is best fitted by two slopes, the break occurring beyond the ocular ring. The inner slope is comparable to the one observed in isolated normal star-forming galaxies but the outer slope is shallow. We present a numerical simulation of the interaction that reproduces the observed tidal features, kinematics, and metallicity distribution, to investigate the effect of the interaction on the galaxies. The model indicates that the galaxies have undergone a close encounter, strongly prograde for the primary, and are half way in their course to a second close encounter.

preprint2020arXiv

The Evolution of Magellanic-like Galaxy Pairs and The Production of Magellanic Stream Analogues in Simulations with Tides, Ram Pressure, and Stellar Feedback

We present a series of chemodynamical simulations of Magellanic-like systems consisting of two interacting, equal-mass dwarf galaxies orbiting a massive host galaxy, including feedback and star formation, tides, and ram pressure. We study the star formation and chemical enrichment history of the dwarfs, and the production of a Magellanic Stream analogue. The dwarfs interact with each other through tidal forces, distorting their morphologies and triggering star formation. A stream is naturally produced as outflows, induced by feedback and interactions, are stretched by tidal forces. Counter to some recent simulations, we find that the leading arm persists even in the presence of ram pressure from the host galaxy. Interactions between the dwarfs and the host galaxies produce multiple kinematic components in the stream, as observed. A combination of ongoing star-formation and entrained low-metallicity gas causes the stream to have a complex chemical structure, with an average metallicity that is generally lower than that of the dwarfs.

preprint2016arXiv

Chemodynamic evolution of dwarf galaxies in tidal fields

The mass-metallicity relation shows that the galaxies with the lowest mass have the lowest metallicities. As most dwarf galaxies are in group environments, interaction effects such as tides could contribute to this trend. We perform a series of smoothed particle hydrodynamics (SPH) simulations of dwarf galaxies in external tidal fields to examine the effects of tides on their metallicities and metallicity gradients. In our simulated galaxies, gravitational instabilities drive gas inwards and produce centralized star formation and a significant metallicity gradient. Strong tides can contribute to these instabilities, but their primary effect is to strip the outer low-metallicity gas, producing a truncated gas disk with a large metallicity. This suggests that the role of tides on the mass-metallicity relation is to move dwarf galaxies to higher metallicities.

preprint2016arXiv

Metal diffusion in smoothed particle hydrodynamics simulations of dwarf galaxies

We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows, and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths, by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]-[Fe/H] distribution in the gas, and cuts off the gas metallicity distribution function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between "metal mass-loading" (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present.

preprint2016arXiv

The Mass Dependence of Star Formation Histories in Barred Spiral Galaxies

We performed a series of 29 gasdynamical simulations of disc galaxies, barred and unbarred, with various stellar masses, to study the impact of the bar on star formation history. Unbarred galaxies evolve very smoothly, with a star formation rate (SFR) that varies by at most a factor of three over a period of 2 Gyr. The evolution of barred galaxies is much more irregular, especially at high stellar masses. In these galaxies, the bar drives a substantial amount of gas toward the centre, resulting in a high SFR, and producing a starburst in the most massive galaxies. Most of the gas is converted into stars, and gas exhaustion leads to a rapid drop of star formation after the starburst. In massive barred galaxies (stellar mass M* > 2x10^10 Msun) the large amount of gas funnelled toward the centre is completely consumed by the starburst, while in lower-mass barred galaxies it is only partially consumed. Gas concentration is thus higher in lower-mass barred galaxies than it is in higher-mass ones. Even though unbarred galaxies funnelled less gas toward their centre, the lower SFR allows this gas to accumulate. At late times, the star formation efficiency is higher in barred galaxies than unbarred ones, enabling these galaxies to maintain a higher SFR with a smaller gas supply. Several properties, such as the global SFR, central SFR, or central gas concentration, vary monotonically with time for unbarred galaxies, but not for barred galaxies. Therefore one must be careful when comparing barred and unbarred galaxies that share one observational property, since these galaxies might be at very different stages of their respective evolution.

preprint2014arXiv

Major Cluster Mergers and the Location of the Brightest Cluster Galaxy

Using a large N-body cosmological simulation combined with a subgrid treatment of galaxy formation, we study the formation and evolution of the galaxy and cluster population in a comoving volume (100 Mpc)^3 in a LCDM universe. At z = 0, our computational volume contains 1788 clusters with mass M_cl > 1.1x10^12 Msun, including 18 massive clusters with M_cl > 10^14 Msun. It also contains 1 088 797 galaxies with mass M_gal > 2x10^9 Msun and luminosity L > 9.5x10^5 Lsun. For each cluster, we identified the brightest cluster galaxy (BCG). We then computed the fraction f_BNC of clusters in which the BCG is not the closest galaxy to the center of the cluster in projection, and the ratio Dv/s, where Dv is the difference in radial velocity between the BCG and the whole cluster, and s is the radial velocity dispersion of the cluster. f_BNC increases from 0.05 for low-mass clusters (M_cl ~ 10^12 Msun) to 0.5 for high-mass ones (M_cl > 10^14 Msun), with no dependence on cluster redshift. The values of Dv/s vary from 0 to 1.8. These results are consistent with previous observational studies, and indicate that the central galaxy paradigm, which states that the BCG should be at rest at the center of the cluster, is usually valid, but exceptions are too common to be ignored. Analysis of the merger trees for the 18 most massive clusters in the simulation reveals that 16 of these clusters have experienced major mergers in the past. These mergers leave each cluster in a non-equilibrium state, but eventually the cluster settles into an equilibrium configuration, unless it is disturbed by another major merger. We found evidence that these mergers are responsible for the off-center positions and peculiar velocities of some BCGs. Our results thus support the merging-group scenario, in which some clusters form by the merger of smaller groups in which the galaxies have already formed.

preprint2013arXiv

Investigating Bar Structure of Disc Galaxies via PRIMAL: A PaRtIcle-by-particle M2M ALgorithm

We have modified our particle-by-particle adaptation of the made-to-measure (M2M) method, with the aim of modelling the Galactic disc from upcoming Galactic stellar survey data. In our new particle-by-particle M2M algorithm, PRIMAL, the observables of the target system are compared with those of the model galaxy at the position of the target stars, i.e. particles. The mass of the model particles are adjusted to reproduce the observables of the target system, and the gravitational potential is automatically adjusted by the changing mass of the particles. This paper builds upon our previous work, introducing likelihood-based velocity constraints in PRIMAL. In this paper we apply PRIMAL to barred disc galaxies created by a N-body simulation in a known dark matter potential, with no error in the observables. This paper demonstrates that PRIMAL can recover the radial profiles of the surface density, velocity dispersion in the radial and perpendicular directions, and the rotational velocity of the target discs, along with the apparent bar structure and pattern speed of the bar, especially when the reference frame is adjusted so that the bar angle of the target galaxy is aligned to that of the model galaxy at every timestep.

preprint2013arXiv

The Connection between Star Formation and Metallicity Evolution in Barred Spiral Galaxies

We have performed a series of chemodynamical simulations of barred disc galaxies. Our goal is to determine the physical processes responsible for the increase in the central gas-phase metallicity and of the central SFR observed in SDSS. We define a 2kpc diameter central aperture to approximate the integrated spectroscopic fibre measurements from the SDSS. The chemical evolution observed within this central region depends critically upon the relative size of the bar and the aperture, which evolves strongly with time. At t~0.5Gyr, a strong bar forms, whose length is considerably longer than the 2kpc aperture. The stars and gas lose angular momentum and follow elongated orbits that cause an intense mixing of the gas between the central region and its surroundings. During the next 1.5Gyr, the orbits of the gas contract significantly until the entire gas bar is contained in the 2kpc aperture, resulting in a net flux of gas into the central region. During this period, the metallicity in the central region increases steadily, and this enrichment is dominated by metal-rich gas that is flowing into the central region. The main result of this work is therefore that the observed enrichment in the centres of barred galaxies is not dominated by in-situ enrichment by stars formed in the centre. Rather, star formation occurs along the full length of the bar, much of which occurs initially outside the 2kpc aperture. About 50% of the metals that end up in the central region originate from this extended bar-long star formation, but flow into the central region due to loss of angular momentum. We conclude that there is no direct connection between central SFR and central metallicity. The central metallicity does not originate exclusively from central stars. Instead, the global SFR (especially along the bar) and the large-scale flow of enriched gas play a major role.

preprint2012arXiv

Fragmentation and Evolution of Molecular Clouds. III: The Effect of Dust and Gas Energetics

Dust and gas energetics are incorporated into a cluster-scale simulation of star formation in order to study the effect of heating and cooling on the star formation process. We build on our previous work by calculating separately the dust and gas temperatures. The dust temperature is set by radiative equilibrium between heating by embedded stars and radiation from dust. The gas temperature is determined using an energy-rate balance algorithm which includes molecular cooling, dust-gas collisional energy transfer, and cosmic-ray ionization. The fragmentation proceeds roughly similarly to simulations in which the gas temperature is set to the dust temperature, but there are differences. The structure of regions around sink particles have properties similar to those of Class 0 objects, but the infall speeds and mass accretion rates were, on average, higher than those seen for regions forming only low-mass stars. The gas and dust temperature have complex distributions not well modeled by approximations that ignore the detailed thermal physics. There is no simple relationship between density and kinetic temperature. In particular, high density regions have a large range of temperatures, determined by their location relative to heating sources. The total luminosity underestimates the star formation rate at these early stages, before ionizing sources are included, by an order of magnitude. As predicted in our previous work, a larger number of intermediate mass objects form when improved thermal physics is included, but the resulting IMF still has too few low mass stars. However, if we consider recent evidence on core-to-star efficiencies, the match to the IMF is improved.

preprint2012arXiv

The Fate of Dwarf Galaxies in Clusters and the Origin of Intracluster Stars. II. Cosmological Simulations

We combine a N-body simulation algorithm with a subgrid treatment of galaxy formation, mergers, and tidal destruction, and an observed conditional luminosity function Phi(L|M), to study the origin and evolution of galactic and extragalactic light inside a cosmological volume of size (100 Mpc)^3, in a concordance LCDM model. This algorithm simulates the growth of large-scale structures and the formation of clusters, the evolution of the galaxy population in clusters, the destruction of galaxies by mergers and tides, and the evolution of the intracluster light. We find that destruction of galaxies by mergers dominates over destruction by tides by about an order of magnitude at all redshifts. However, tidal destruction is sufficient to produce intracluster light fractions f_ICL that are sufficiently high to match observations. The bulk of the intracluster light (60%) is provided by intermediate galaxies of total masses 10^11 Msun-10^12 Msun and stellar masses 6x10^8 Msun-3x10^10 Msun that were tidally destroyed by even more massive galaxies. The contribution of low-mass galaxies to the intracluster light is small and the contribution of dwarf galaxies is negligible, even though, by numbers, most galaxies that are tidally destroyed are dwarfs. Tracking clusters back in time, we find that their values of f_ICL tend to increase over time, but can experience sudden changes that are sometimes non-monotonic. These changes occur during major mergers involving clusters of comparable masses but very different intracluster luminosities. Most of the tidal destruction events take place in the central regions of clusters. As a result, the intracluster light is more centrally concentrated than the galactic light. Our results support tidal destruction of intermediate-mass galaxies as a plausible scenario for the origin of the intracluster light.

preprint2011arXiv

Chemical signature of gas-rich disc-disc mergers at high redshift

We performed numerical simulations of mergers between gas-rich disc galaxies, which result in the formation of late-type galaxies. Stars formed during the merger end up in a thick disc that is partially supported by velocity dispersion and has high [alpha/Fe] ratios at all metallicities. Stars formed later end up in a thin, rotationally supported disc which has lower [alpha/Fe] ratios. While the structural and kinematical properties of the merger remnants depend strongly upon the orbital parameters of the mergers, we find a clear chemical signature of gas-rich mergers.

preprint2011arXiv

Galactic Outflows and Evolution of the Interstellar Medium

We present a model to self-consistently describe the joint evolution of starburst galaxies and the galactic wind resulting from this evolution. We combine the population synthesis code Starburst99 with a semi-analytical model of galactic outflows and a model for the distribution and abundances of chemical elements inside the outflows. Starting with a galaxy mass, formation redshift, and adopting a particular form for the star formation rate, we describe the evolution of the stellar populations in the galaxy, the evolution of the metallicity and chemical composition of the interstellar medium (ISM), the propagation of the galactic wind, and the metal-enrichment of the intergalactic medium (IGM). In this paper, we study the properties of the model, by varying the mass of the galaxy, the star formation rate, and the efficiency of star formation. Our main results are the following: (1) For a given star formation efficiency f*, a more extended period of active star formation tends to produce a galactic wind that reaches a larger extent. If f* is sufficiently large, the energy deposited by the stars completely expels the ISM. Eventually, the ISM is being replenished by mass loss from supernovae and stellar winds. (2) For galaxies with masses above 10^11 Msun, the material ejected in the IGM always falls back onto the galaxy. Hence lower-mass galaxies are the ones responsible for enriching the IGM. (3) Stellar winds play a minor role in the dynamical evolution of the galactic wind, because their energy input is small compared to supernovae. However, they contribute significantly to the chemical composition of the galactic wind. We conclude that the history of the ISM enrichment plays a determinant role in the chemical composition and extent of the galactic wind, and therefore its ability to enrich the IGM.

preprint2011arXiv

Numerical Simulations of Galactic Outflows and Evolution of the IGM

Galactic outflows play a major role in the evolution of galaxies and the intergalactic medium (IGM). The energy deposited into the interstellar medium by supernovae and active galactic nuclei can accelerate the gas past the escape velocity, and eject it into the IGM. This will affect the subsequent evolution of the galaxy, by reducing or eliminating star formation, and quenching the accretion of matter onto the central AGN. Galactic outflows is the main process by which energy and processed interstellar matter is transported into the IGM. This affects the subsequent formation of other galaxies. The energy carried by outflows can strip protogalactic halos of their gas, preventing galaxies from forming. Conversely, the metals carried by outflows can modify the composition and cooling rates of the gas in protogalactic halos, favoring the formation of galaxies. In this paper, I review the various techniques used to simulate galactic outflows and their impact on galaxy and IGM evolution.

preprint2010arXiv

Anisotropic AGN Outflows and Enrichment of the Intergalactic Medium. II. Metallicity

We investigate the large-scale influence of outflows from AGNs in enriching the IGM with metals in a cosmological context. We combine cosmological simulations of large scale structure formation with a detailed model of metal enrichment, in which outflows expand anisotropically along the direction of least resistance, distributing metals into the IGM. The metals carried by the outflows are generated by two separate stellar populations: stars located near the central AGN, and stars located in the greater galaxy. Using this algorithm, we performed a series of 5 simulations of the propagation of AGN-driven outflows in a cosmological volume of size (128/h Mpc)^3 in a Lambda-CDM universe, and analyze the resulting metal enrichment of the IGM. We found that the metallicity induced in the IGM is greatly dominated by AGNs having bolometric luminosity L > 10^9 L_sun, sources with 10^8 < L / L_sun < 10^9 having a negligible contribution. Our simulations produced an average IGM metallicity of [O/H] = -5 at z = 5.5, which then rises gradually, and remains relatively flat at a value [O/H] = -2.8 between z = 2 and z = 0. The ejection of metals from AGN host galaxies by AGN-driven outflows is found to enrich the IGM to > 10 - 20% of the observed values, the number dependent on redshift. The enriched IGM volume fractions are small at z > 3, then rise rapidly to the following values at z = 0: 6 - 10% of the volume enriched to [O/H] > -2.5, 14 - 24% volume to [O/H] > -3, and 34 - 45% volume to [O/H] > -4. At z > 2, there is a gradient of the induced enrichment, the metallicity decreasing with increasing IGM density, enriching the underdense IGM to higher metallicities, a trend more prominent with increasing anisotropy of the outflows. This can explain observations of metal-enriched low-density IGM at z = 3 - 4.

preprint2010arXiv

Anisotropic Galactic Outflows and Enrichment of the Intergalactic Medium. II. Numerical Simulations

We combine an analytic model for anisotropic outflows and galaxy formation with numerical simulations of large-scale structure and halo formation to study the impact of galactic outflows on the evolution of the IGM. We have simulated the evolution of a comoving volume (15 Mpc)^3 in the LCDM universe. We follow the formation of 20000-60000 galaxies and simulate the galactic outflows produced by these galaxies, for five outflow opening angles, alpha=60, 90, 120, 150, and 180 degrees (isotropic outflows). Anisotropic outflows follow the path of least resistance and thus travel preferentially into low-density regions, away from cosmological structures where galaxies form. These anisotropic outflows are less likely to overlap with one another, or to hit pre-galactic collapsing halos and strip them of their gas, preventing a galaxy from forming. Going from 180 deg to 60 deg, the number of galaxies that actually form doubles, producing twice as many outflows, and these outflows overlap to a lesser extent. As a result, the metal volume filling factor of the IGM goes from 8% for isotropic outflows up to 28% for anisotropic ones. High density regions are more efficiently enriched than low density ones (~80% compared to ~20% by volume), even though most enriched regions are low densities. Increasing the anisotropy of outflows increases the extent of enrichment at all densities, low and high. This is in part because anisotropic outflows are more numerous. When this effect is factored-out, we find that the probability a galaxy will enrich systems at densities up to 10 rho_mean is higher for increasingly anisotropic outflows. This is an effect of the dynamical evolution of the IGM. Anisotropic outflows expand preferentially into underdense gas, but that gas can later accrete onto overdense structures.

preprint2009arXiv

Fragmentation and Evolution of Molecular Clouds. II: The Effect of Dust Heating

We investigate the effect of heating by luminosity sources in a simulation of clustered star formation. Our heating method involves a simplified continuum radiative transfer method that calculates the dust temperature. The gas temperature is set by the dust temperature. We present the results of four simulations, two simulations assume an isothermal equation of state and the two other simulations include dust heating. We investigate two mass regimes, i.e., 84 Msun and 671 Msun, using these two different energetics algorithms. The mass functions for the isothermal simulations and simulations which include dust heating are drastically different. In the isothermal simulation, we do not form any objects with masses above 1 Msun. However, the simulation with dust heating, while missing some of the low-mass objects, forms high-mass objects (~20 Msun) which have a distribution similar to the Salpeter IMF. The envelope density profiles around the stars formed in our simulation match observed values around isolated, low-mass star-forming cores. We find the accretion rates to be highly variable and, on average, increasing with final stellar mass. By including radiative feedback from stars in a cluster-scale simulation, we have determined that it is a very important effect which drastically affects the mass function and yields important insights into the formation of massive stars.

preprint2009arXiv

Structure, Kinematics, and Chemical Enrichment Patterns after Major Gas-Rich Disc-Disc Mergers

We used an N-body smoothed particle hydrodynamics algorithm, with a detailed treatment of star formation, supernovae feedback, and chemical enrichment, to perform eight simulations of mergers between gas-rich disc galaxies. We vary the mass ratio of the progenitors, their rotation axes, and their orbital parameters and analyze the kinematic, structural, and chemical properties of the remnants. Six of these simulations result in the formation of a merger remnant with a disc morphology as a result of the large gas-fraction of the remnants. We show that stars formed during the merger (a sudden starburst occur in our simulation and last for 0.2-0.3 Gyr) and those formed after the merger have different kinematical and chemical properties. The first ones are located in thick disc or the halo. They are partially supported by velocity dispersion and have high [alpha/Fe] ratios even at metallicities as high as [Fe/H]=-0.5. The former ones -- the young component -- are located in a thin disc rotationally supported and have lower [alpha/Fe] ratios. The difference in the rotational support of both components results in the rotation of the thick disc lagging that of the thin disc by as much as a factor of two, as recently observed.We find that, while the kinematic and structural properties of the merger remnant depends strongly upon the orbital parameters of the mergers, there is a remarkable uniformity in the chemical properties of the mergers. This suggests that general conclusions about the chemical signature of gas-rich mergers can be drawn.

preprint2006arXiv

The Spatial Distribution of the Galactic First Stars II: SPH Approach

We use cosmological, chemo-dynamical, smoothed particle hydrodynamical simulations of Milky-Way-analogue galaxies to find the expected present-day distributions of both metal-free stars that formed from primordial gas and the oldest star populations. We find that metal-free stars continue to form until z~4 in halos that are chemically isolated and located far away from the biggest progenitor of the final system. As a result, if the Population III initial mass function allows stars with low enough mass to survive until z=0 (< 0.8 Msol), they would be distributed throughout the Galactic halo. On the other hand, the oldest stars form in halos that collapsed close to the highest density peak of the final system, and at z=0 they are located preferentially in the central region of the Galaxy, i.e., in the bulge. According to our models, these trends are not sensitive to the merger histories of the disk galaxies or the implementation of supernova feedback. Furthermore, these full hydrodynamics results are consistent with our N-body results in Paper I, and lend further weight to the conclusion that surveys of low-metallicity stars in the Galactic halo can be used to directly constrain the properties of primordial stars. In particular, they suggest that the current lack of detections of metal-free stars implies that their lifetimes were shorter than a Hubble time, placing constraints on the metal-free initial mass function.

preprint2005arXiv

Gas Rich Mergers in Disk Formation

In order to explain disk galaxy formation within the hierarchical structure formation, it seems that gas rich mergers must play an important role. We review here our previous studies which have shown the importance of mergers at high redshift being gas rich, in the formation of both the stellar halo and thick disk components of disk galaxies. Regulation of star formation in the building blocks of our galaxy is required to form a low mass low metallicity stellar halo. This regulation results in high redshift, gas rich mergers during which the thick disk forms. In these proceedings, we categorise stars from our simulated disk galaxy into thin and thick disk components by using the Toomre diagram.

preprint2005arXiv

Understanding the Equilibrium Structure of CDM Halos

N-body simulations find a universal structure for the halos which result from the nonlinear growth of Gaussian-random-noise density fluctuations in the CDM universe. This talk summarized our attempts to derive and explain this universal structure by analytical approximation and simplified models. As an example, we show here that a 1D spherical infall model involving a fluid approximation derived from the Boltzmann equation can explain not only the halo density profile but its phase-space density profile, as well.

preprint1998arXiv

Light Propagation in inhomogeneous Universes

Using a multi-plane lensing method that we have developed, we follow the evolution of light beams as they propagate through inhomogeneous universes. We use a P3M code to simulate the formation and evolution of large-scale structure. The resolution of the simulations is increased to sub-Megaparsec scales by using a Monte Carlo method to locate galaxies inside the computational volume according to the underlying particle distribution. The galaxies are approximated by isothermal spheres, with each morphological type having its own distribution of masses and core radii. The morphological types are chosen in order to reproduce the observed morphology-density relation. This algorithm has an effective resolution of 9 orders of magnitudes in length, from the size of superclusters down to the core radii of the smallest galaxies. We consider cold dark matter models normalized to COBE, and perform a large parameter survey by varying the cosmological parameters Omega_0, lambda_0, H_0, and n (the tilt of the primordial power spectrum). The values of n are chosen by imposing particular values or sigma_8, the rms mass fluctuation at a scale of 8/h Mpc. We use the power spectrum given by Bunn & White. This is the largest parameter survey ever done is this field.

preprint1997arXiv

Adaptive Smoothed Particle Hydrodynamics: Methodology II

This paper presents an alternative formulation of the ASPH algorithm for evolving anisotropic smoothing kernels, in which the geometric approach of Shapiro et al. (1996; Paper I) is replaced by an approach involving a local transformation of coordinates to those in which the underlying anisotropic volume changes appear to be isotropic. The ASPH method is presented in 2D and 3D, including a number of details not previously included in Paper I, some of which represent either advances or different choices with respect to Paper I. Among the advances included here are an asynchronous time-integration scheme with different time steps for different particles and the generalization of the ASPH method to 3D. The shock-tracking algorithm described in Paper I for locally adapting the artificial viscosity to restrict viscous heating just to particles encountering shocks, is not included here. Instead, we adopt a different interpolation kernel for use with the artificial viscosity, which has the effect of spatially localizing effects of the artificial viscosity. This version of the ASPH method in 2D and 3D is then applied to a series of 1D, 2D, and 3D test problems, and the results are compared to those of standard SPH applied to the same problems. These include the problem of cosmological pancake collapse, the Riemann shock tube, cylindrical and spherical Sedov blast waves, the collision of two strong shocks, and problems involving shearing disks intended to test the angular momentum conservation properties of the method. These results further support the idea that ASPH has significantly better resolving power than standard SPH for a wide range of problems, including that of cosmological structure formation. (Abridged)