Source author record

Florent Renaud

Florent Renaud 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

32works
4topics
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

32 published item(s)

preprint2023arXiv

Pre-supernova feedback sets the star cluster mass function to a power law and reduces the cluster formation efficiency

The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrain star formation and as fundamental units for stellar feedback models. Detailed galaxy models must therefore accurately capture the basic properties of observed clusters to be considered predictive. We use hydrodynamical simulations of a dwarf galaxy as a laboratory to study star cluster formation. We test different combinations of stellar feedback mechanisms, including stellar winds, ionizing radiation, and supernovae. Each feedback mechanism affects the cluster formation efficiency and cluster mass function. Increasing the feedback budget by combining the different types of feedback decreases the cluster formation efficiency by reducing the number of massive clusters. Ionizing radiation is found to be especially influential. This effect depends on the timing of feedback initiation, as shown by comparing early and late feedback. Early feedback occurs from ionizing radiation and stellar winds with onset immediately after a massive star is formed. Late feedback occurs when energy injection only starts after the main-sequence lifetime of the most massive SN progenitor, a timing that is further influenced by the choice of the most massive SN progenitor. Late feedback alone results in a broad, flat mass function, approaching a log-normal shape in the complete absence of feedback. Early feedback, on the other hand, produces a power-law cluster mass function with lower formation efficiency, albeit with a steeper slope than that usually observed.

preprint2022arXiv

EDGE: The sensitivity of ultra-faint dwarfs to a metallicity-dependent initial mass function

Motivated by the observed bottom-light initial mass function (IMF) in faint dwarfs, we study how a metallicity-dependent IMF affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions, which we complement with an empirically motivated subgrid prescription for systematic IMF variations. In this framework, at the low gas metallicities typical of faint dwarfs, the IMF of newborn stellar populations becomes top-heavy, increasing the efficiency of supernova and photoionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is nonetheless met by increased metal production from more numerous massive stars, leading to nearly constant iron content at $z=0$. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star}\sim 10^3 \rm M_{\odot}$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass--metallicity relation.

preprint2022arXiv

From giant clumps to clouds -- III. The connection between star formation and turbulence in the ISM

Supersonic gas turbulence is a ubiquitous property of the interstellar medium. The level of turbulence, quantified by the gas velocity dispersion ($σ_{\rm g}$), is observed to increase with the star formation rate (SFR) of a galaxy, but it is yet not established whether this trend is driven by stellar feedback or gravitational instabilities. In this work we carry out hydrodynamical simulations of entire disc galaxies, with different gas fractions, to understand the origins of the SFR-$σ_{\rm g}$ relation. We show that disc galaxies reach the same levels of turbulence regardless of the presence of stellar feedback processes, and argue that this is an outcome of the way disc galaxies regulate their gravitational stability. The simulations match the SFR-$σ_{\rm g}$ relation up to SFRs of the order of tens of M$_\odot$ yr$^{-1}$ and $σ_{\rm g}\sim 50$ km s$^{-1}$ in neutral hydrogen and molecular gas, but fail to reach the very large values ($>100$ km s$^{-1}$) reported in the literature for rapidly star-forming galaxies. We demonstrate that such high values of $σ_{\rm g}$ can be explained by 1) insufficient beam smearing corrections in observations, and 2) stellar feedback being coupled to the ionised gas phase traced by recombination lines. Given that the observed SFR-$σ_{\rm g}$ relation is composed of highly heterogeneous data, with $σ_{\rm g}$ at high SFRs almost exclusively being derived from H$α$ observations of high redshift galaxies with complex morphologies, we caution against analytical models that attempt to explain the SFR-$σ_{\rm g}$ relation without accounting for these effects.

preprint2022arXiv

Stellar feedback in M83 as observed with MUSE -- I. Overview, an unprecedented view of the stellar and gas kinematics and evidence of outflowing gas

We present a large VLT/MUSE mosaic (3.8 x 3.8 kpc) of the nearby spiral galaxy M83, with a spatial resolution ~20 pc. We obtained the kinematics of the stars and ionised gas, and compared them with molecular gas kinematics from ALMA CO(2-1). We separated the ionised gas into HII regions and diffuse ionised gas (DIG) and determined the fraction of Ha luminosity originating from the DIG (f_DIG). We observe that both stars and gas trace the galactic disk rotation, as well as a fast-rotating nuclear component, likely connected to secular processes driven by the galactic bar. In the gas kinematics, we observe a stream east of the nucleus, redshifted with respect to the disk. The stream is surrounded by an extended ionised gas region with enhanced velocity dispersion and a high ionisation state, which is largely consistent with being ionised by slow shocks. We interpret this feature as either the superposition of the disk and an extraplanar layer of DIG, or as a bar-driven inflow of shocked gas. A double Gaussian component fit to the Ha line also reveals the presence of a nuclear biconic structure whose axis of symmetry is perpendicular to the bar. The two cones appear blue- and redshifted along the line of sight and stand out for having an Ha emission separated by up to 200 km s-1 from that of the disk, and a high velocity dispersion ~80-200 km s-1. At the far end of the cones, we observe that the gas is consistent with being ionised by shocks. These features had never been observed before in M83; we postulate that they are tracing a starburst-driven outflow shocking into the surrounding ISM. Finally, we obtain f_DIG ~ 13% in our field of view. We inspect the emission of the HII regions and DIG in `BPT' diagrams, finding that in HII regions photoionisation accounts for 99.8% of the Ha flux, whereas the DIG has a mixed contribution from photoionisation (94.9%) and shocks (5.1%). [abridged]

preprint2021arXiv

From giant clumps to clouds -- II. The emergence of thick disc kinematics from the conditions of star formation in high redshift gas rich galaxies

High redshift disc galaxies are more gas rich, clumpier, and more turbulent than local Universe galaxies. This early era of galaxy formation imprints the distribution and kinematics of the stars that we observe today, but it is not yet well established how. In this work, we use simulations of isolated Milky Way-mass disc galaxies to study how kinematic properties of stars change when varying the gas fraction. This allows us to quantify the roles played by internal processes, e.g. gas turbulence and gravitational scattering off massive gas clumps, in establishing the observed stellar velocity dispersions and orbital eccentricities. We find that models with gas fractions $>20$ per cent feature a turbulent and clumpy interstellar medium (ISM), leading to zero-age stellar velocity dispersions $\sim 20-30~{\rm km\, s}^{-1}$ and high mean orbital eccentricities. Low eccentricities cannot arise from these physical conditions. For gas fractions below $20$ per cent, the ISM becomes less turbulent, with stellar velocity dispersions $<10~{\rm km\, s}^{-1}$, and nearly circular orbits for young stars. The turbulence present in gas-rich high redshift galaxies hence acts as a `barrier' against the formation of thin discs. We compare our findings to the Milky Way's age-velocity dispersion relation and argue that velocity dispersions imprinted already at star formation by the ISM contribute significantly at all times. Finally, we show that observed orbital eccentricities in the Milky Way's thick and thin discs can be explained entirely as imprints by the star-forming ISM, rather than by mergers or secular processes.

preprint2020arXiv

How runaway stars boost galactic outflows

Roughly ten per cent of OB stars are kicked out of their natal clusters before ending their life as supernovae. These so-called runaway stars can travel hundreds of parsecs into the low-density interstellar medium, where momentum and energy from stellar feedback is efficiently deposited. In this work we explore how this mechanism affects large scale properties of the galaxy, such as outflows. To do so we use a new model which treats OB stars and their associated feedback processes on a star-by-star basis. With this model we compare two hydrodynamical simulations of Milky Way-like galaxies, one where we include runaways, and one where we ignore them. Including runaway stars leads to twice as many supernovae explosions in regions with gas densities ranging from 1e-5 cm^-3 to 1e-3 cm^-3. This results in more efficient heating of the inter-arm regions, and drives strong galactic winds with mass loading factors boosted by up to one order of magnitude. These outflows produce a more massive and extended multi-phase circumgalactic medium, as well as a population of dense clouds in the halo. Conversely, since less energy and momentum is released in the dense star forming regions, the cold phase of the interstellar medium is less disturbed by feedback effects.

preprint2019arXiv

3 things they don't tell you about star clusters

Dense stellar systems in general and star clusters in particular have recently regained the interest of the extragalactic and even cosmology communities, due to the role they could play as actors and probes of re-ionization, galactic archeology and the dark matter content of galaxies, among many others. In the era of the exploitation and the preparation of large stellar surveys (Gaia, APOGEE, 4MOST, WEAVE), of the detection of gravitational waves mostly originating from dense regions like the cores of clusters (Ligo, LISA), and in an always more holistic view of galaxy formation (HARMONI, Euclid, LSST, soon to be known as the Vera Rubin Survey Telescope, VRST), a complete theory on the formation and evolution of clusters is needed to interpret the on-going and forthcoming data avalanche. In this context, the community carries an effort to model the aspects of star cluster formation and evolution in galactic and even cosmological context. However, it is not always easy to understand the caveats and the shortcuts taken in theories and simulations, and their implications on the conclusions drawn. I take the opportunity of this document to highlight three of these topics and discuss why some shortcuts taken by the community are or could be misleading.

preprint2016arXiv

If it does not kill them, it makes them stronger: collisional evolution of star clusters with tidal shocks

The radii of young (<100 Myr) star clusters correlate only weakly with their masses. This shallow relation has been used to argue that impulsive tidal perturbations, or `shocks', by passing giant molecular clouds (GMCs) preferentially disrupt low-mass clusters. We show that this mass-radius relation is in fact the result of the combined effect of two-body relaxation and repeated tidal shocks. Clusters in a broad range of environments including those like the solar neighbourhood evolve towards a typical radius of a few parsecs, as observed, independent of the initial radius. This equilibrium mass-radius relation is the result of a competition between expansion by relaxation and shrinking due to shocks. Interactions with GMCs are more disruptive for low-mass clusters, which helps to evolve the globular cluster mass function (GCMF). However, the properties of the interstellar medium in high-redshift galaxies required to establish a universal GCMF shape are more extreme than previously derived, challenging the idea that all GCs formed with the same power-law mass function.

preprint2016arXiv

New insights on the formation of nuclear star clusters

Nuclear Clusters (NCs) are common stellar systems in the centres of galaxies. Yet, the physical mechanisms involved in their formation are still debated. Using a parsec-resolution hydrodynamical simulation of a dwarf galaxy, we propose an updated formation scenario for NCs. In this 'wet migration scenario', a massive star cluster forms in the gas-rich disc, keeping a gas reservoir, and growing further while it migrates to the centre via a combination of interactions with other substructures and dynamical friction. A wet merger with another dense cluster and its own gas reservoir can occur, although this is not a pre-requisite for the actual formation of the NC. The merging process does significantly alter the properties of the NC (mass, morphology, star formation history), also quenching the on-going local star formation activity, thus leading to interesting observational diagnostics for the physical origin of NCs. A population of lower mass clusters co-exist during the simulation, but these are either destroyed via tidal forces, or have high angular momentum preventing them to interact with the NC and contribute to its growth. The proposed updated scenario emphasises the role of gas reservoirs associated with the densest star clusters formed in a gas-rich low-mass galaxy.

preprint2016arXiv

Star formation triggered by galaxy interactions in modified gravity

Together with interstellar turbulence, gravitation is one key player in star formation. It acts both at galactic scales in the assembly of gas into dense clouds, and inside those structures for their collapse and the formation of pre-stellar cores. To understand to what extent the large scale dynamics govern the star formation activity of galaxies, we present hydrodynamical simulations in which we generalise the behaviour of gravity to make it differ from Newtonian dynamics in the low acceleration regime. We focus on the extreme cases of interacting galaxies, and compare the evolution of galaxy pairs in the dark matter paradigm to that in the Milgromian Dynamics (MOND) framework. Following up on the seminal work by Tiret & Combes, this paper documents the first simulations of galaxy encounters in MOND with a detailed Eulerian hydrodynamical treatment of baryonic physics, including star formation and stellar feedback. We show that similar morphologies of the interacting systems can be produced by both the dark matter and MOND formalisms, but require a much slower orbital velocity in the MOND case. Furthermore, we find that the star formation activity and history are significantly more extended in space and time in MOND interactions, in particular in the tidal debris. Such differences could be used as observational diagnostics and make interacting galaxies prime objects in the study of the nature of gravitation at galactic scales.

preprint2016arXiv

The origin of the Milky Way globular clusters

We present a cosmological zoom-in simulation of a Milky Way-like galaxy used to explore the formation and evolution of star clusters. We investigate in particular the origin of the bimodality observed in the colour and metallicity of globular clusters, and the environmental evolution through cosmic times in the form of tidal tensors. Our results self-consistently confirm previous findings that the blue, metal-poor clusters form in satellite galaxies which are accreted onto the Milky Way, while the red, metal-rich clusters form mostly in situ or, to a lower extent in massive, self-enriched galaxies merging with the Milky Way. By monitoring the tidal fields these populations experience, we find that clusters formed in situ (generally centrally concentrated) feel significantly stronger tides than the accreted ones, both in the present-day, and when averaged over their entire life. Furthermore, we note that the tidal field experienced by Milky Way clusters is significantly weaker in the past than at present-day, confirming that it is unlikely that a power-law cluster initial mass function like that of young massive clusters, is transformed into the observed peaked distribution in the Milky Way with relaxation-driven evaporation in a tidal field.

preprint2015arXiv

A flexible method to evolve collisional systems and their tidal debris in external potentials

We introduce a numerical method to integrate tidal effects on collisional systems, using any definition of the external potential as a function of space and time. Rather than using a linearisation of the tidal field, this new method follows a differential technique to numerically evaluate the tidal acceleration and its time derivative. Theses are then used to integrate the motions of the components of the collisional systems, like stars in star clusters, using a predictor-corrector scheme. The versatility of this approach allows the study of star clusters, including their tidal tails, in complex, multi-components, time-evolving external potentials. The method is implemented in the code nbody6 (Aarseth 2003).

preprint2015arXiv

The effect of secular galactic growth on the evolution of star clusters

The growth of galaxies through adiabatic accretion of dark matter is one of the main drivers of galaxy evolution. By isolating it from other processes like mergers, we analyse how it affects the evolution of star clusters. Our study comprises a fast and approximate exploration of the orbital and intrinsic cluster parameter space, and more detailed monitoring of their evolution, through N-body simulations for a handful of cases. We find that the properties of present-day star clusters and their tidal tails differ very little, whether the clusters are embedded in a growing galactic halo for 12 Gyr, or in a static one.

preprint2015arXiv

The incorrect rotation curve of the Milky Way

In the fundamental quest of the rotation curve of the Milky Way, the tangent-point (TP) method has long been the simplest way to infer velocities for the inner, low latitude regions of the Galactic disk from observations of the gas component. We test the validity of the method on realistic gas distribution and kinematics of the Milky Way, using a numerical simulation of the Galaxy. We show that the resulting velocity profile strongly deviates from the true rotation curve of the simulation, as it overstimates it in the central regions, and underestimates it around the bar corotation. Also, its shape strongly depends on the orientation of the stellar bar. The discrepancies are caused by highly non-uniform azimuthal velocities, and the systematic selection by the TP method of high-velocity gas along the bar and spiral arms, or low-velocity gas in less dense regions. The velocity profile is in good agreement with the rotation curve only beyond corotation, far from massive asymmetric structures. Therefore the observed velocity profile of the Milky Way inferred by the TP method is expected to be very close to the true Galactic rotation curve for 4.5<R<8 kpc. Another consequence is that the Galactic velocity profile for R<4-4.5 kpc is very likely flawed by the non-uniform azimuthal velocities, and does not represent the true Galactic rotation curve, but instead local motions. The real shape of the innermost rotation curve is probably shallower than previously thought. Using a wrong rotation curve has a dramatic impact on the modelling of the mass distribution, in particular for the bulge component of which derived enclosed mass within the central kpc and scale radius are, respectively, twice and half of the actual values. We thus strongly argue against using terminal velocities or the velocity curve from the TP method for modelling the mass distribution of the Milky Way. (abridged)

preprint2014arXiv

A parsec-resolution simulation of the Antennae galaxies: Formation of star clusters during the merger

We present a hydrodynamical simulation of an Antennae-like galaxy merger at parsec resolution, including a multi-component model for stellar feedback and reaching numerical convergence in the global star formation rate for the first time. We analyse the properties of the dense stellar objects formed during the different stages of the interaction. Each galactic encounter triggers a starburst activity, but the varying physical conditions change the triggering mechanism of each starburst. During the first two pericenter passages, the starburst is spatially extended and forms many star clusters. However, the starburst associated to the third, final passage is more centrally concentrated: stars form almost exclusively in the galactic nucleus and no new star cluster is formed. The maximum mass of stars clusters in this merger is more than 30 times higher than those in a simulation of an isolated Milky Way-like galaxy. Antennae-like mergers are therefore a formation channel of young massive clusters possibly leading to globular clusters. Monitoring the evolution of a few clusters reveals the diversity of formation scenarios including the gathering and merger of gas clumps, the monolithic formation and the hierarchical formation in sub-structures inside a single cloud. Two stellar objects formed in the simulation yield the same properties as ultra-compact dwarf galaxies. They share the same formation scenario than the most massive clusters, but have a larger radius either since birth, or get it after a violent interaction with the galactic center. The diversity of environments across space and time in a galaxy merger can account for the diversity of the stellar objects formed, both in terms of mass and size.

preprint2014arXiv

Starbursts triggered by inter-galactic tides and interstellar compressive turbulence

Using parsec-resolution simulations of a typical galaxy merger, we study the triggering of starbursts by connecting the (inter-)galactic dynamics to the structure of the interstellar medium. The gravitational encounter between two galaxies enhances tidal compression over large volumes, which increases and modifies the turbulence, in particular its compressive mode with respect to the solenoidal one. This generates an excess of dense gas leading to intense star formation activity. Along the interaction, the compressive turbulence modifies the efficiency of gas-to-star conversion which, in the Schmidt-Kennicutt diagram, drives the galaxies from the sequence of discs to that of starbursts.

preprint2014arXiv

The inefficiency of satellite accretion in forming extended star clusters

The distinction between globular clusters and dwarf galaxies has been progressively blurred by the recent discoveries of several extended star clusters, with size (20-30 pc) and luminosity (-6 < Mv < -2) comparable to the one of faint dwarf spheroidals. In order to explain their sparse structure, it has been suggested that they formed as star clusters in dwarf galaxy satellites that later accreted onto the Milky Way. If these clusters form in the centre of dwarf galaxies, they evolve in a tidally-compressive environment where the contribution of the tides to the virial balance can become significant, and lead to a super-virial state and subsequent expansion of the cluster, once removed. Using N-body simulations, we show that a cluster formed in such an extreme environment undergoes a sizable expansion, during the drastic variation of the external tidal field due to the accretion process. However, we show that the expansion due to the removal of the compressive tides is not enough to explain the observed extended structure, since the stellar systems resulting from this process are always more compact than the corresponding clusters that expand in isolation due to two-body relaxation. We conclude that an accreted origin of extended globular clusters is unlikely to explain their large spatial extent, and rather favor the hypothesis that such clusters are already extended at the stage of their formation.

preprint2014arXiv

The interplay between a galactic bar and a supermassive black hole: nuclear fueling in a sub-parsec resolution galaxy simulation

We study the connection between the large-scale dynamics and the gas fueling toward a central black hole via the analysis of a Milky Way-like simulation at sub-parsec resolution. This allows us to follow a set of processes at various scales (e.g., the triggering of inward gas motion towards inner resonances via the large-scale bar, the connection to the central black hole via mini spirals) in a self-consistent manner. This simulation provides further insights on the role of shear for the inhibition of star formation within the bar in regions with significant amount of gas. We also witness the decoupling of the central gas and nuclear cluster from the large-scale disc, via interactions with the black hole. This break of symmetry in the mass distribution triggers the formation of gas clumps organised in a time-varying 250 pc ring-like structure, the black hole being offset by about 70 pc from its centre. Some clumps form stars, while most get disrupted or merge. Supernovae feedback further creates bubbles and filaments, some of the gas being expelled to 100 pc or higher above the galaxy plane. This helps remove angular momentum from the gas, which gets closer to the central dark mass. Part of the gas raining down is being accreted, forming a 10~pc polar disc-like structure around the black hole, leading to an episode of star formation. This gives rise to multiple stellar populations with significantly different angular momentum vectors, and may lead to a natural intermittence in the fueling of the black hole.

preprint2014arXiv

The role of turbulence in star formation laws and thresholds

The Schmidt-Kennicutt relation links the surface densities of gas to the star formation rate in galaxies. The physical origin of this relation, and in particular its break, i.e. the transition between an inefficient regime at low gas surface densities and a main regime at higher densities, remains debated. Here, we study the physical origin of the star formation relations and breaks in several low-redshift galaxies, from dwarf irregulars to massive spirals. We use numerical simulations representative of the Milky Way, the Large and the Small Magellanic Clouds with parsec up to subparsec resolution, and which reproduce the observed star formation relations and the relative variations of the star formation thresholds. We analyze the role of interstellar turbulence, gas cooling, and geometry in drawing these relations, at 100 pc scale. We suggest in particular that the existence of a break in the Schmidt- Kennicutt relation could be linked to the transition from subsonic to supersonic turbulence and is independent of self-shielding effects. This transition being connected to the gas thermal properties and thus to the metallicity, the break is shifted toward high surface densities in metal-poor galaxies, as observed in dwarf galaxies. Our results suggest that together with the collapse of clouds under self-gravity, turbulence (injected at galactic scale) can induce the compression of gas and regulate star formation.

preprint2013arXiv

The long lives of giant clumps and the birth of outflows in gas-rich galaxies at high redshift

Star-forming disk galaxies at high redshift are often subject to violent disk instability, characterized by giant clumps whose fate is yet to be understood. The main question is whether the clumps disrupt within their dynamical timescale (<50Myr), like molecular clouds in today's galaxies, or whether they survive stellar feedback for more than a disk orbital time (~300Myr) in which case they can migrate inward and help building the central bulge. We present 3.5-7pc resolution AMR simulations of high-redshift disks including photo-ionization, radiation pressure, and supernovae feedback (Renaud et al. 2013, and Perret et al., this astro-ph issue). Our modeling of radiation pressure determines the mass loading and initial velocity of winds from basic physical principles. We find that the giant clumps produce steady outflow rates comparable to and sometimes somewhat larger than their star formation rate, with velocities largely sufficient to escape galaxy. The clumps also lose mass, especially old stars, by tidal stripping, and the stellar populations contained in the clumps hence remain relatively young (<=200Myr), as observed. The clumps survive gaseous outflows and stellar loss, because they are wandering in gas-rich turbulent disks from which they can re-accrete gas at high rates compensating for outflows and tidal stripping, overall keeping realistic and self-regulated gaseous and stellar masses. Our simulations produce gaseous outflows with velocities, densities and mass loading consistent with observations, and at the same time suggest that the giant clumps survive for hundreds of Myr and complete their migration to the center of high-redshift galaxies, without rapid dispersion and reformation of clumps.

preprint2013arXiv

The role of galaxy mergers on the evolution of star clusters

Interacting galaxies favor the formation of star clusters but are also suspected to affect their evolution through an intense and rapidly varying tidal field. Treating this complex behaviour remains out-of-reach of (semi-)analytical studies. By computing the tidal field from galactic models and including it into star-by-star N-body simulations of star clusters, we monitor the structure and mass evolution of a population of clusters in a galaxy major merger, taking the Antennae galaxies (NGC 4038/39) as a prototype. On the long timescale (~ 10^9 yr), the merger only indirectly affects the evolution of clusters by modifying their orbits in or around the galaxies: the mass-loss of clusters in the merger remnant is faster, while clusters ejected in the tidal debris survive much longer, compared to in an isolated galaxy. The tidal perturbations of the galactic collisions themselves are too short lived and not strong enough to significantly influence the structure and dissolution of realistically dense/massive star clusters.

preprint2013arXiv

Volume Rendering of AMR Simulations

High-resolution simulations often rely on the Adaptive Mesh Resolution (AMR) technique to optimize memory consumption versus attainable precision. While this technique allows for dramatic improvements in terms of computing performance, the analysis and visualization of its data outputs remain challenging. The lack of effective volume renderers for the octree-based AMR used by the RAMSES simulation program has led to the development of the solutions presented in this paper. Two custom algorithms are discussed, based on the splatting and the ray-casting techniques. Their usage is illustrated in the context of the visualization of a high-resolution, 6000-processor simulation of a Milky Way-like galaxy. Performance obtained in terms of memory management and parallelism speedup are presented.

preprint2012arXiv

Star formation laws and thresholds from ISM structure and turbulence

We present an analytical model of the relation between the surface density of gas and star formation rate in galaxies and clouds, as a function of the presence of supersonic turbulence and the associated structure of the interstellar medium. The model predicts a power-law relation of index 3/2, flattened under the effects of stellar feedback at high densities or in very turbulent media, and a break at low surface densities when ISM turbulence becomes too weak to induce strong compression. This model explains the diversity of star formation laws and thresholds observed in nearby spirals and their resolved regions, the Small Magellanic Cloud, high-redshift disks and starbursting mergers, as well as Galactic molecular clouds. While other models have proposed interstellar dust content and molecule formation to be key ingredients to the observed variations of the star formation efficiency, we demonstrate instead that these variations can be explained by interstellar medium turbulence and structure in various types of galaxies.

preprint2011arXiv

Evolution of star clusters in arbitrary tidal fields

We present a novel and flexible tensor approach to computing the effect of a time-dependent tidal field acting on a stellar system. The tidal forces are recovered from the tensor by polynomial interpolation in time. The method has been implemented in a direct-summation stellar dynamics integrator (NBODY6) and test-proved through a set of reference calculations: heating, dissolution time and structural evolution of model star clusters are all recovered accurately. The tensor method is applicable to arbitrary configurations, including the important situation where the background potential is a strong function of time. This opens up new perspectives in stellar population studies reaching to the formation epoch of the host galaxy or galaxy cluster, as well as for star-burst events taking place during the merger of large galaxies. A pilot application to a star cluster in the merging galaxies NGC 4038/39 (the Antennae) is presented.

preprint2011arXiv

Models of Stephan's Quintet: Hydrodynamic Constraints on the Group's Evolution

We present smoothed particle hydrodynamic models of the interactions in the compact galaxy group, Stephan's Quintet. This work is extension of the earlier collisionless N-body simulations of Renaud et al. in which the large-scale stellar morphology of the group was modeled with a series of galaxy-galaxy interactions in the simulations. Including thermohydrodynamic effects in this work, we further investigate the dynamical interaction history and evolution of the intergalactic gas of Stephan's Quintet. The major features of the group, such as the extended tidal features and the group-wide shock, enabled us to constrain the models reasonably well, while trying to reproduce multiple features of the system. We found that reconstructing the two long tails extending from NGC 7319 toward NGC 7320c one after the other in two separate encounters is very difficult and unlikely, because the second encounter usually destroys or distorts the already-generated tidal structure. Our models suggest the two long tails may be formed simultaneously from a single encounter between NGC 7319 and 7320c, resulting in a thinner and denser inner tail than the outer one. The tails then also run parallel to each other as observed. The model results support the ideas that the group-wide shock detected in multi-wavelength observations between NGC 7319 and 7318b and the starburst region north of NGC 7318b are triggered by the high-speed collision between NGC 7318b and the intergalactic gas. Our models show that a gas bridge is formed by the high-speed collision and clouds in the bridge continue to interact for some tens of millions of years after the impact. This produces many small shocks in that region, resulting a much longer cooling time than that of a single impact shock.

preprint2011arXiv

Tides in colliding galaxies

Long tails and streams of stars are the most noticeable upshots of galaxy collisions. Their origin as gravitational, tidal, disturbances has however been recognized only less than fifty years ago and more than ten years after their first observations. This Review describes how the idea of galactic tides emerged, in particular thanks to the advances in numerical simulations, from the first ones that included tens of particles to the most sophisticated ones with tens of millions of them and state-of-the-art hydrodynamical prescriptions. Theoretical aspects pertaining to the formation of tidal tails are then presented. The third part of the review turns to observations and underlines the need for collecting deep multi-wavelength data to tackle the variety of physical processes exhibited by collisional debris. Tidal tails are not just stellar structures, but turn out to contain all the components usually found in galactic disks, in particular atomic / molecular gas and dust. They host star-forming complexes and are able to form star-clusters or even second-generation dwarf galaxies. The final part of the review discusses what tidal tails can tell us (or not) about the structure and content of present-day galaxies, including their dark components, and explains how tidal tails may be used to probe the past evolution of galaxies and their mass assembly history. On-going deep wide-field surveys disclose many new low-surface brightness structures in the nearby Universe, offering great opportunities for attempting galactic archeology with tidal tails.

preprint2010arXiv

Dynamics of the Tidal Fields and Formation of Star Clusters in Galaxy Mergers

In interacting galaxies, strong tidal forces disturb the global morphology of the progenitors and give birth to the long stellar, gaseous and dusty tails often observed. In addition to this destructive effect, tidal forces can morph into a transient, protective setting called compressive mode. Such modes then shelter the matter in their midst by increasing its gravitational binding energy. This thesis focuses on the study of this poorly known regime by quantifying its properties thanks to numerical and analytical tools applied to a spectacular merging system of two galaxies, commonly known as the Antennae galaxies. N-body simulations of this pair yield compressive modes in the regions where observations reveal a burst of star formation. Furthermore, characteristic time- and energy scales of these modes match well those of self-gravitating substructures such as star clusters and tidal dwarf galaxies. These results suggest that the compressive modes of tidal fields plays an important role in the formation and evolution of young clusters, at least in a statistical sense, over a lapse of ~10 million years. Preliminary results from simulations of stellar associations highlight the importance of embedding the clusters in the evolving background galaxies to account precisely for their morphology and internal evolution.

preprint2010arXiv

N-body simulation of the Stephan's Quintet

The evolution of compact groups of galaxies may represent one of the few places in the nearby universe in which massive galaxies are being forged through a complex set of processes involving tidal interaction, ram-pressure stripping, and perhaps finally "dry-mergers" of galaxies stripped of their cool gas. Using collisionless N-body simulations, we propose a possible scenario for the formation of one of the best studied compact groups: Stephan's Quintet. We define a serial approach which allows us to consider the history of the group as sequence of galaxy-galaxy interactions seen as relatively separate events in time, but chained together in such a way as to provide a plausible scenario that ends in the current configuration of the galaxies. By covering a large set of parameters, we claim that it is very unlikely that both major tidal tails of the group have been created by the interaction between the main galaxy and a single intruder. We propose instead a scenario based on two satellites orbiting the main disk, plus the recent involvement of an additional interloper, coming from the background at high speed. This purely N-body study of the quintet will provide a parameter-space exploration of the basic dynamics of the group that can be used as a basis for a more sophisticated N-body/hydrodynamic study of the group that is necessary to explain the giant shock structure and other purely gaseous phenomena observed in both the cold, warm and hot gas in the group.

preprint2010arXiv

One moment in time - modeling star formation in the Antennae

We present a new high-resolution N-body/SPH simulation of an encounter of two gas-rich disk galaxies which closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation and feedback from SNII. The large-scale morphology and kinematics are determined by the internal structure and the orbit of the progenitor disks. The properties of the central region, in particular the starburst in the overlap region, only match the observations for a very short time interval of ~20 Myr after the second encounter. This indicates that the Antennae galaxies are in a special phase only about 40 Myr after the second encounter and 50 Myr before their final collision. This is the only phase in the simulation when a gas-rich overlap region between the nuclei is forming accompanied by enhanced star formation. The star formation rate as well as the recent star formation history in the central region agree well with observational estimates. For the first time this new model explains the distributed extra-nuclear star formation in the Antennae galaxies as a consequence of the recent second encounter. The proposed model predicts that the Antennae are in a later merger stage than the Mice (NGC 4676) and would therefore lose their first place in the classical Toomre sequence.

preprint2009arXiv

Fully compressive tides in galaxy mergers

The disruptive effect of galactic tides is a textbook example of gravitational dynamics. However, depending on the shape of the potential, tides can also become fully compressive. When that is the case, they might trigger or strengthen the formation of galactic substructures (star clusters, tidal dwarf galaxies), instead of destroying them. We perform N-body simulations of interacting galaxies to quantify this effect. We demonstrate that tidal compression occurs repeatedly during a galaxy merger, independently of the specific choice of parameterization. With a model tailored to the Antennae galaxies, we show that the distribution of compressive tides matches the locations and timescales of observed substructures. After extending our study to a broad range of parameters, we conclude that neither the importance of the compressive tides (~15% of the stellar mass) nor their duration (~ 10 Myr) are strongly affected by changes in the progenitors' configurations and orbits. Moreover, we show that individual clumps of matter can enter compressive regions several times in the course of a simulation. We speculate that this may spawn multiple star formation episodes in some star clusters, through e.g., enhanced gas retention.

preprint2009arXiv

The mass-to-light ratio of rich star clusters

We point out a strong time-evolution of the mass-to-light conversion factor eta commonly used to estimate masses of unresolved star clusters from observed cluster spectro-photometric measures. We present a series of gas-dynamical models coupled with the Cambridge stellar evolution tracks to compute line-of-sight velocity dispersions and half-light radii weighted by the luminosity. We explore a range of initial conditions, varying in turn the cluster mass and/or density, and the stellar population's IMF. We find that eta, and hence the estimated cluster mass, may increase by factors as large as 3 over time-scales of 50 million years. We apply these results to an hypothetic cluster mass distribution function (d.f.) and show that the d.f. shape may be strongly affected at the low-mass end by this effect. Fitting truncated isothermal (Michie-King) models to the projected light profile leads to over-estimates of the concentration parameter c of delta c ~ 0.3 compared to the same functional fit applied to the projected mass density.