Florent Renaud

A sub-parsec resolution simulation of the Milky Way: global structure of the interstellar medium and properties of molecular clouds

We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy at a resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through photoionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at sub-parsec resolution for a few cloud lifetimes and at 0.05 pc for about a cloud-crossing time. The turbulence cascade and gravitation from kpc scales are de facto included in smaller structures like molecular clouds. We show that the formation of a bar influences the dynamics of the central ˜100 pc by creating resonances. At larger radii, the spiral arms host the formation of regularly spaced clouds: beads on a string and spurs. These instabilities pump turbulent energy into the gas, generally in the supersonic regime. Because of asymmetric drift, the supernovae explode outside their gaseous nursery, which diminishes the effect of feedback on the structure of clouds. The evolution of clouds is thus mostly due to fragmentation and gas consumption, regulated mainly by supersonic turbulence. The transition from turbulence-supported to self-gravitating gas is detected in the gas density probability distribution function at ˜2000 cm-3. The power-spectrum density suggests that gravitation governs the hierarchical organization of structures from the galactic scale down to a few pc.

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 (\textless= 50 Myr), like the molecular clouds in todays galaxies, or whether they survive stellar feedback for more than a disk orbital time (approximate to 300 Myr) in which case they can migrate inward and help building the central bulge. We present 3.5-7 pc resolution adaptive mesh refinement simulations of high-redshift disks including photoionization, radiation pressure, and supernovae feedback. 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 the galaxy. The clumps also lose mass, especially old stars, by tidal stripping, and the stellar populations contained in the clumps hence remain relatively young (\textless= 200 Myr), as observed. The clumps survive gaseous outflows and stellar loss, because they are wandering in gas-rich turbulent disks from which they can reaccrete gas at high rates compensating for outflows and tidal stripping, overall keeping realistic and self-regulated gaseous and stellar masses. The outflow and accretion rates have specific timescales of a few 10(8) yr, as opposed to rapid and repeated dispersion and reformation of clumps. 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. These long-lived clumps are gas-dominated and contain a moderate mass fraction of stars; they drive inside-out disk evolution, thickening, spheroid growth, and fueling of the central black hole.

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 starburst 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.

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.

One moment in time - modeling star formation in the antennae

We present a new high-resolution N-body/smoothed particle hydrodynamics simulation of an encounter of two gas-rich disk galaxies that closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation, and feedback from Type II supernovae. 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 (Δt ≈ 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. (Less)

The interplay between a galactic bar and a supermassive black hole: nuclear fuelling in a subparsec 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.

STAR FORMATION LAWS AND THRESHOLDS FROM INTERSTELLAR MEDIUM 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.

Evolution of the mass, size, and star formation rate in high redshift merging galaxies - MIRAGE – A new sample of simulations with detailed stellar feedback

We aim at addressing the questions related to galaxy mass assembly through major and minor wet merging processes in the redshift range 1<z<2. A consequent fraction of Milky Way like galaxies are thought to have undergone an unstable clumpy phase at this early stage. Using the adaptive mesh refinement code RAMSES, with a recent physically-motivated implementation of stellar feedback, we build the Merging and Isolated high-Redshift Adaptive mesh refinement Galaxies (MIRAGE) sample. It is composed of 20 mergers and 3 isolated idealized disks simulations with global physical properties in accordance with the 1<z<2 mass complete sample MASSIV. The numerical hydrodynamical resolution reaches 7 parsecs in the smallest Eulerian cells. Our simulations include: star formation, metal line cooling, metallicity advection, and a recent implementation of stellar feedback which encompasses OB-type stars radiative pressure, photo-ionization heating, and supernovae. The initial conditions are set to match the z~2 observations, thanks to a new public code DICE. The numerical resolution allows us to follow the formation and evolution of giant clumps formed in-situ from Jeans instabilities triggered by high initial gas fraction. The star formation history of isolated disks shows stochastic star formation rate, which proceeds from the complex behavior of the giant clumps. Our minor and major gas-rich merger simulations do not trigger starbursts, suggesting a saturation of the star formation in a turbulent and clumpy interstellar medium fed by substantial accretion from the circum-galactic medium. Our simulations are close to the normal regime of the disk-like star formation on a Schmidt-Kennicutt diagram. The mass-size relation and its rate of evolution matches observations, suggesting that the inside-out growth mechanisms of the stellar disk do not necessarily require to be achieved through a cold accretion.

The role of galaxy mergers on the evolution of star clusters

Interacting galaxies favour 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 time-scale (∼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.

Modeling CO emission from hydrodynamic simulations of nearby spirals, starbursting mergers, and high-redshift galaxies

We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and highredshift galaxies with very high resolutions (3 pc and 10 3 M� ) and detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes, and galactic winds. The simulations are analyzed with a large velocity gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer, opacity effect, and the intensity emerging from galaxies to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially concentrated stellar feedback sources, although CO(1−0) to CO(3−2) emission is dominated overall by low-excitation gas around the densest clumps. These results generally highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn affect the CO emission properties.