-
Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics
Authors:
Simon Portegies Zwart,
Inti Pelupessy,
Carmen Martinez-Barbosa,
Arjen van Elteren,
Steve McMillan
Abstract:
Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method…
▽ More
Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method is efficient for solving the equations of motion for a self-gravitating system over a wide range of scales. We adopt a dedicated integrator for solving each particular part of the problem and combine the results to a self-consistent solution. In particular, we explore the possibilities of combining the evolution of one or more microscopic systems that are embedded in a macroscopic system. The here presented generalizations of Bridge include higher-order coupling strategies (from the classic 2nd order up to 10th-order), but we also demonstrate how multiple bridges can be nested and how additional processes can be introduced at the bridge time-step to enrich the physics, for example by incorporating dissipative processes. Such augmentation allows for including additional processes in a classic Newtonian N-body integrator without alterations to the underlying code. These additional processes include for example the Yarkovsky effect, dynamical friction or relativistic dynamics. Some of these processes operate on all particles whereas others apply only to a subset.
The presented method is non-intrusive in the sense that the underlying methods remain operational without changes to the code (apart from adding the get- and set-functions to enable the bridge operator). As a result, the fundamental integrators continue to operate with their internal time step and preserve their local optimizations and parallelism.
... abridged ...
△ Less
Submitted 25 February, 2020;
originally announced February 2020.
-
Survivability of planetary systems in young and dense star clusters
Authors:
Arjen van Elteren,
Simon Portegies Zwart,
Inti Pelupessy,
Maxwell Cai,
Steve McMillan
Abstract:
We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced…
▽ More
We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced. Four, five, or size planets per star were introduced in orbit around the 500 solar-like stars with a maximum orbital separation of 400au. Our study focuses on the production of free-floating planets. A total of 357 become unbound from a total of 2522 planets in the initial conditions of the simulation. Of these, 281 leave the cluster within the crossing timescale of the star cluster; the others remain bound to the cluster as free-floating intra-cluster planets. Five of these free-floating intra-cluster planets are captured at a later time by another star. The two main mechanisms by which planets are lost from their host star, ejection upon a strong encounter with another star or internal planetary scattering, drive the evaporation independent of planet mass of orbital separation at birth. The effect of small perturbations due to slow changes in the cluster potential are important for the evolution of planetary systems. In addition, the probability of a star to lose a planet is independent of the planet mass and independent of its initial orbital separation. As a consequence, the mass distribution of free-floating planets indistinguishable from the mass distribution of planets bound to their host star.
△ Less
Submitted 11 March, 2019; v1 submitted 12 February, 2019;
originally announced February 2019.
-
Gone with the wind: the impact of wind mass transfer on the orbital evolution of AGB binary systems
Authors:
M. I. Saladino,
O. R. Pols,
E. van der Helm,
I. Pelupessy,
S. Portegies Zwart
Abstract:
In low-mass binary systems, mass transfer is likely to occur via a slow and dense stellar wind when one of the stars is in the AGB phase. Observations show that many binaries that have undergone AGB mass transfer have orbital periods of 1-10 yr, at odds with the predictions of binary population synthesis models. We investigate the mass-accretion efficiency and angular-momentum loss via wind mass t…
▽ More
In low-mass binary systems, mass transfer is likely to occur via a slow and dense stellar wind when one of the stars is in the AGB phase. Observations show that many binaries that have undergone AGB mass transfer have orbital periods of 1-10 yr, at odds with the predictions of binary population synthesis models. We investigate the mass-accretion efficiency and angular-momentum loss via wind mass transfer in AGB binary systems. We use these quantities to predict the evolution of the orbit. We perform 3D hydrodynamical simulations of the stellar wind lost by an AGB star using the AMUSE framework. We approximate the thermal evolution of the gas by imposing a simple effective cooling balance and we vary the orbital separation and the velocity of the stellar wind. We find that for wind velocities $v_{\infty}$ larger than the relative orbital velocity of the system $v_\mathrm{orb}$ the flow is described by the Bondi-Hoyle-Lyttleton approximation and the angular-momentum loss is modest, leading to an expansion of the orbit. For low wind velocities an accretion disk is formed around the companion and the accretion efficiency as well as the angular-momentum loss are enhanced, implying that the orbit will shrink. We find that the transfer of angular momentum from the orbit to the outflowing gas occurs within a few orbital separations from the center of mass of the binary. Our results suggest that the orbital evolution of AGB binaries can be predicted as a function of the ratio $v_{\infty}/v_\mathrm{orb}$. Our results can provide insight into the puzzling orbital periods of post-AGB binaries and suggest that the number of stars entering into the common-envelope phase will increase. The latter can have significant implications for the expected formation rates of the end products of low-mass binary evolution, such as cataclysmic binaries, type Ia supernova and double white-dwarf mergers. [ABRIDGED]
△ Less
Submitted 8 May, 2018;
originally announced May 2018.
-
The consequences of a nearby supernova on the early Solar System
Authors:
Simon Portegies Zwart,
Inti Pelupessy,
Arjen van Elteren,
Thomas Wijnen,
Maria Lugaro
Abstract:
If the Sun was born in a relatively compact open cluster, it is quite likely that a massive (10MSun) star was nearby when it exploded in a supernova. The repercussions of a supernova can be rather profound, and the current Solar System may still bear the memory of this traumatic event. The truncation of the Kuiper belt and the tilt of the ecliptic plane with respect to the Sun's rotation axis coul…
▽ More
If the Sun was born in a relatively compact open cluster, it is quite likely that a massive (10MSun) star was nearby when it exploded in a supernova. The repercussions of a supernova can be rather profound, and the current Solar System may still bear the memory of this traumatic event. The truncation of the Kuiper belt and the tilt of the ecliptic plane with respect to the Sun's rotation axis could be such signatures. We simulated the effect of a nearby supernova on the young Solar System using the Astronomical Multipurpose Software Environment. Our calculations are realized in two subsequent steps in which we study the effect of the supernova irradiation on the circumstellar disk and the effect of the impact of the nuclear blast-wave which arrives a few decades later. We find that the blastwave of our adopted supernova exploding at a distance of $0.15$--$0.40$\,pc and at an angle of $35^\circ$--$65^\circ$ with respect to the angular-momentum axis of the circumsolar disk would induce a misalignment between the Sun's equator and its disk to $5^\circ.6\pm1^\circ.2$, consistent with the current value. The blast of a supernova truncates the disk at a radius between $42$ and $55$\,au, which is consistent with the current edge of the Kuiper belt. For the most favored parameters, the irradiation by the supernova as well as the blast wave heat the majority of the disk to $\sim 1200$\,K, which is sufficiently hot to melt chondrules in the circumstellar disk. The majority of planetary system may have been affected by a nearby supernova, some of its repercussions, such as truncation and tilting of the disk, may still be visible in their current planetary system's topology. The amount of material from the supernova blast wave that is accreted by the circumstellar disk is too small by several orders of magnitude to explain the current abundance of the short live radionuclide $^{26}$Al.
△ Less
Submitted 11 April, 2018; v1 submitted 12 February, 2018;
originally announced February 2018.
-
The origin of interstellar asteroidal objects like 1I/2017 U1 'Oumuamua
Authors:
Simon Portegies Zwart,
Santiago Torres,
Inti Pelupessy,
Jeroen Bedorf,
Maxwell Cai
Abstract:
We study the origin of the interstellar object 1I/2017 U1 'Oumuamua by juxtaposing estimates based on the observations with simulations. We speculate that objects like 'Oumuamua are formed in the debris disc as left over from the star and planet formation process, and subsequently liberated. The liberation process is mediated either by interaction with other stars in the parental star-cluster, by…
▽ More
We study the origin of the interstellar object 1I/2017 U1 'Oumuamua by juxtaposing estimates based on the observations with simulations. We speculate that objects like 'Oumuamua are formed in the debris disc as left over from the star and planet formation process, and subsequently liberated. The liberation process is mediated either by interaction with other stars in the parental star-cluster, by resonant interactions within the planetesimal disc or by the relatively sudden mass loss when the host star becomes a compact object. Integrating backward in time in the Galactic potential together with stars from the Gaia-TGAS catalogue we find that about 1.3Myr ago 'Oumuamua passed the nearby star HIP 17288 within a mean distance of $1.3$pc. By comparing nearby observed L-dwarfs with simulations of the Galaxy we conclude that the kinematics of 'Oumuamua is consistent with relatively young objects of $1.1$--$1.7$Gyr. We just met 'Oumuamua by chance, and with a derived mean Galactic density of $\sim 3\times 10^{5}$ similarly sized objects within 100\,au from the Sun or $\sim 10^{14}$ per cubic parsec we expect about 2 to 12 such visitors per year within 1au from the Sun.
△ Less
Submitted 11 May, 2018; v1 submitted 9 November, 2017;
originally announced November 2017.
-
Capture of exocomets and the erosion of the Oort cloud due to stellar encounters in the Galaxy
Authors:
Joris Hanse,
Lucie Jilkova,
Simon F. Portegies Zwart,
F. Inti Pelupessy
Abstract:
The Oort cloud (OC) probably formed more than 4$\,$Gyr ago and has been moving with the Sun in the Galaxy since, exposed to external influences, most prominently to the Galactic tide and passing field stars. Theories suggest that other stars might posses exocomets distributed similarly to our OC. We study the erosion of the OC and the possibility for capturing exocomets during the encounters with…
▽ More
The Oort cloud (OC) probably formed more than 4$\,$Gyr ago and has been moving with the Sun in the Galaxy since, exposed to external influences, most prominently to the Galactic tide and passing field stars. Theories suggest that other stars might posses exocomets distributed similarly to our OC. We study the erosion of the OC and the possibility for capturing exocomets during the encounters with such field stars. We carry out simulations of flybys, where both stars are surrounded by a cloud of comets. We measure how many exocomets are transferred to the OC, how many OC's comets are lost, and how this depends on the other star's mass, velocity and impact parameter. Exocomets are transferred to the OC only during relatively slow ($\lesssim0.5\,$km$\,$s$^{-1}$) and close ($\lesssim10^5\,$AU) flybys and these are expected to be extremely rare. Assuming that all passing stars are surrounded by a cloud of exocomets, we derive that the fraction of exocomets in the OC has been about $10^{-5}$--$10^{-4}$. Finally we simulate the OC for the whole lifetime of the Sun, taking into account the encounters and the tidal effects. The OC has lost 25--65% of its mass, mainly due to stellar encounters, and at most 10% (and usually much less) of its mass can be captured. However, exocomets are often lost shortly after the encounter that delivers them, due to the Galactic tide and consecutive encounters.
△ Less
Submitted 16 October, 2017;
originally announced October 2017.
-
Changes in orientation and shape of protoplanetary discs moving through an ambient medium
Authors:
T. P. G. Wijnen,
F. I. Pelupessy,
O. R. Pols,
S. Portegies Zwart
Abstract:
Misalignments between the orbital planes of planets and the equatorial planes of their host stars have been observed in our solar system, in transiting exoplanets, and in the orbital planes of debris discs. We present a mechanism that causes such a spin-orbit misalignment for a protoplanetary disc due to its movement through an ambient medium. Our physical explanation of the mechanism is based on…
▽ More
Misalignments between the orbital planes of planets and the equatorial planes of their host stars have been observed in our solar system, in transiting exoplanets, and in the orbital planes of debris discs. We present a mechanism that causes such a spin-orbit misalignment for a protoplanetary disc due to its movement through an ambient medium. Our physical explanation of the mechanism is based on the theoretical solutions to the Stark problem. We test this idea by performing self-consistent hydrodynamical simulations and simplified gravitational $N$-body simulations. The $N$-body model reduces the mechanism to the relevant physical processes. The hydrodynamical simulations show the mechanism in its full extent, including gas-dynamical and viscous processes in the disc which are not included in the theoretical framework. We find that a protoplanetary disc embedded in a flow changes its orientation as its angular momentum vector tends to align parallel to the relative velocity vector. Due to the force exerted by the flow, orbits in the disc become eccentric, which produces a net torque and consequentially changes the orbital inclination. The tilting of the disc causes it to contract. Apart from becoming lopsided, the gaseous disc also forms a spiral arm even if the inclination does not change substantially. The process is most effective at high velocities and observational signatures are therefore mostly expected in massive star-forming regions and around winds or supernova ejecta. Our $N$-body model indicates that the interaction with supernova ejecta is a viable explanation for the observed spin-orbit misalignment in our solar system.
△ Less
Submitted 4 July, 2017;
originally announced July 2017.
-
Disc truncation in embedded star clusters: Dynamical encounters versus face-on accretion
Authors:
T. P. G. Wijnen,
O. R. Pols,
F. I. Pelupessy,
S. Portegies Zwart
Abstract:
Observations indicate that the dispersal of protoplanetary discs in star clusters occurs on time scales of about 5 Myr. Several processes are thought to be responsible for this disc dispersal. Here we compare two of these processes: dynamical encounters and interaction with the interstellar medium, which includes face-on accretion and ram pressure stripping. We perform simulations of embedded star…
▽ More
Observations indicate that the dispersal of protoplanetary discs in star clusters occurs on time scales of about 5 Myr. Several processes are thought to be responsible for this disc dispersal. Here we compare two of these processes: dynamical encounters and interaction with the interstellar medium, which includes face-on accretion and ram pressure stripping. We perform simulations of embedded star clusters with parameterisations for both processes to determine the environment in which either of these processes is dominant. We find that face-on accretion, including ram pressure stripping, is the dominant disc truncation process if the fraction of the total cluster mass in stars is $\lesssim 30\,\%$ regardless of the cluster mass and radius. Dynamical encounters require stellar densities $\gtrsim 10^4$ pc$^{-3}$ combined with a mass fraction in stars of $\approx 90\,\%$ to become the dominant process. Our results show that during the embedded phase of the cluster, the truncation of the discs is dominated by face-on accretion and dynamical encounters become dominant when the intra-cluster gas has been expelled. As a result of face-on accretion the protoplanetary discs become compact and their surface density increases. In contrast, dynamical encounters lead to discs that are less massive and remain larger.
△ Less
Submitted 21 June, 2017;
originally announced June 2017.
-
Characterising face-on accretion onto and the subsequent contraction of protoplanetary discs
Authors:
T. P. G. Wijnen,
O. R. Pols,
F. I. Pelupessy,
S. Portegies Zwart
Abstract:
Observations indicate that stars generally lose their protoplanetary discs on a timescale of about 5 Myr. Which mechanisms are responsible for the disc dissipation is still debated. Here we investigate the movement through an ambient medium as a possible cause of disc dispersal. The ram pressure exerted by the flow can truncate the disc and the accretion of material with no azimuthal angular momen…
▽ More
Observations indicate that stars generally lose their protoplanetary discs on a timescale of about 5 Myr. Which mechanisms are responsible for the disc dissipation is still debated. Here we investigate the movement through an ambient medium as a possible cause of disc dispersal. The ram pressure exerted by the flow can truncate the disc and the accretion of material with no azimuthal angular momentum leads to further disc contraction. We derive a theoretical model from accretion disc theory that describes the evolution of the disc radius, mass, and surface density profile as a function of the density and velocity of the ambient medium. We test our model by performing hydrodynamical simulations of a protoplanetary disc embedded in a flow with different velocities and densities. We find that our model gives an adequate description of the evolution of the disc radius and accretion rate onto the disc. The total disc mass in the simulations follows the theoretically expected trend, except at the lowest density where our simulated discs lose mass owing to continuous stripping. This stripping may be a numerical rather than a physical effect. Some quantitative differences exist between the model predictions and the simulations. These are at least partly caused by numerical viscous effects in the disc and depend on the resolution of the simulation. Our model can be used as a conservative estimate for the process of face-on accretion onto protoplanetary discs, as long as viscous processes in the disc can be neglected. The model predicts that in dense gaseous environments, discs can shrink substantially in size and can, in theory, sweep up an amount of gas of the order of their initial mass. This process could be relevant for planet formation in dense environments.
△ Less
Submitted 14 February, 2017;
originally announced February 2017.
-
Constraining cloud parameters using high density gas tracers in galaxies
Authors:
M. V. Kazandjian,
I. Pelupessy,
R. Meijerink,
F. P. Israel,
C. M. Coppola,
M. J. F. Rosenberg,
M. Spaans
Abstract:
Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic…
▽ More
Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic scales. We do this by re-sampling high density gas in a hydrodynamic model of a gas-rich disk galaxy, assuming that the density field of the interstellar medium of the model galaxy follows the probability density function (PDF) inferred from the resolved low density scales. We find that in a narrow gas density PDF, with a mean density of ~10 cm-3 and a dispersion σ= 2.1 in the log of the density, most of the emission of molecular lines, emanates from the 10-1000 cm-3 part of the PDF. We construct synthetic emission maps for the central 2 kpc of the galaxy and fit the line ratios of CO and 13CO up to J = 15-14, as well as HCN, HNC, and HCO+ up to J = 7-6, using one photo-dissociation region (PDR) model. We attribute the goodness of the one component fits for our model galaxy to the fact that the distribution of the luminosity, as a function of density, is peaked at gas densities between 10 and 1000 cm-3.
We explore the impact of different log-normal density PDFs on the distribution of the line-luminosity as a function of density, and we show that it is necessary to have a broad dispersion, corresponding to Mach numbers >~ 30 in order to obtain significant emission from n > 10^4 cm-3 gas. Such Mach numbers are expected in star-forming galaxies, LIRGS, and ULIRGS. By fitting line ratios of HCN(1-0), HNC(1-0), and HCO+(1-0) for a sample of LIRGS and ULIRGS using mechanically heated PDRs, we constrain the Mach number of these galaxies to 29 < M < 77.
△ Less
Submitted 8 September, 2016;
originally announced September 2016.
-
Modelling Mechanical Heating in Star-Forming Galaxies: CO and 13CO Line Ratios as Sensitive Probes
Authors:
M. V. Kazandjian,
I. Pelupessy,
R. Meijerink,
F. P. Israel,
M. Spaans
Abstract:
We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3.
We consider t…
▽ More
We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3.
We consider two galaxy models: a small disk galaxy of solar metallicity and a lighter dwarf galaxy with 0.2 \zsun metallicity. Elevated excitation temperatures for CO(1 - 0) correlate positively with mechanical feedback, that is enhanced towards the central region of both model galaxies. The emission maps of these model galaxies are used to compute line ratios of CO and 13CO transitions. These line ratios are used as diagnostics where we attempt to match them These line ratios are used as diagnostics where we attempt to match them to mechanically heated single component (i.e. uniform density, Far-UV flux, visual extinction and velocity gradient) equilibrium PDR models.
We find that PDRs ignoring mechanical feedback in the heating budget over-estimate the gas density by a factor of 100 and the far-UV flux by factors of ~10 - 1000. In contrast, PDRs that take mechanical feedback into account are able to fit all the line ratios for the central < 2 kpc of the fiducial disk galaxy quite well. The mean mechanical heating rate per H atom that we recover from the line ratio fits of this region varies between $10^{-27}$ -- $10^{-26}$~erg s$^{-1}$.
△ Less
Submitted 8 September, 2016;
originally announced September 2016.
-
Face-on accretion onto a protoplanetary disc
Authors:
T. P. G. Wijnen,
O. R. Pols,
F. I. Pelupessy,
S. Portegies Zwart
Abstract:
Globular clusters (GCs) are known to harbor multiple stellar populations. To explain these observations Bastian et al. suggested a scenario in which a second population is formed by the accretion of enriched material onto the low-mass stars in the initial GC population. The idea is that the low-mass, pre-main sequence stars sweep up gas expelled by the massive stars of the same generation into the…
▽ More
Globular clusters (GCs) are known to harbor multiple stellar populations. To explain these observations Bastian et al. suggested a scenario in which a second population is formed by the accretion of enriched material onto the low-mass stars in the initial GC population. The idea is that the low-mass, pre-main sequence stars sweep up gas expelled by the massive stars of the same generation into their protoplanetary disc as they move through the GC core. We perform simulations with 2 different smoothed particle hydrodynamics codes to investigate if a low-mass star surrounded by a protoplanetary disc can accrete the amount of enriched material required in this scenario. We focus on the gas loading rate onto the disc and star as well as on the lifetime of the disc. We find that the gas loading rate is a factor of 2 smaller than the geometric rate, because the effective cross section of the disc is smaller than its surface area. The loading rate is consistent for both codes, irrespective of resolution. The disc gains mass in the high resolution runs, but loses angular momentum on a time scale of 10^4 yrs. Two effects determine the loss of (specific) angular momentum in our simulations: 1) continuous ram pressure stripping and 2) accretion of material with no azimuthal angular momentum. Our study and previous work suggest that the former, dominant process is mainly caused by numerical rather than physical effects, while the latter is not. The latter process causes the disc to become more compact, increasing the surface density profile at smaller radii. The disc size is determined in the first place by the ram pressure when the flow first hits the disc. Further evolution is governed by the decrease in the specific angular momentum of the disc. We conclude that the size and lifetime of the disc are probably not sufficient to accrete the amount of mass required in Bastian et al.'s scenario.
△ Less
Submitted 4 July, 2016;
originally announced July 2016.
-
The Origin of Chaos in the Orbit of Comet 1P/Halley
Authors:
Tjarda C. N. Boekholt,
F. Inti Pelupessy,
Douglas C. Heggie,
Simon F. Portegies Zwart
Abstract:
According to Munoz-Gutierrez et al. (2015) the orbit of comet 1P/Halley is chaotic with a surprisingly small Lyapunov time scale of order its orbital period. In this work we analyse the origin of chaos in Halley's orbit and the growth of perturbations, in order to get a better understanding of this unusually short time scale. We perform N-body simulations to model Halley's orbit in the Solar Syste…
▽ More
According to Munoz-Gutierrez et al. (2015) the orbit of comet 1P/Halley is chaotic with a surprisingly small Lyapunov time scale of order its orbital period. In this work we analyse the origin of chaos in Halley's orbit and the growth of perturbations, in order to get a better understanding of this unusually short time scale. We perform N-body simulations to model Halley's orbit in the Solar System and measure the separation between neighbouring trajectories. To be able to interpret the numerical results, we use a semi-analytical map to demonstrate different growth modes, i.e. linear, oscillatory or exponential, and transitions between these modes. We find the Lyapunov time scale of Halley's orbit to be of order 300 years, which is significantly longer than previous estimates in the literature. This discrepancy could be due to the different methods used to measure the Lyapunov time scale. A surprising result is that next to Jupiter, also encounters with Venus contribute to the exponential growth in the next 3000 years. Finally, we note an interesting application of the sub-linear, oscillatory growth mode to an ensemble of bodies moving through the Solar System. Whereas in the absence of encounters with a third body the ensemble spreads out linearly in time, the accumulation of weak encounters can increase the lifetime of such systems due to the oscillatory behaviour.
△ Less
Submitted 22 June, 2016;
originally announced June 2016.
-
Stellar winds near massive black holes: The case of the S-stars
Authors:
Nora Lützgendorf,
Edwin van der Helm,
Inti Pelupessy,
Simon Portegies Zwart
Abstract:
The Galactic center provides a unique laboratory to study the interaction of a supermassive black hole (SMBH) with its gaseous and stellar environment. Simulations to determine the accretion of stellar winds from the surrounding O-stars onto the black hole have been performed earlier, but in those the presence of the S-star system was ignored. The S-stars are a group of young massive B-stars in re…
▽ More
The Galactic center provides a unique laboratory to study the interaction of a supermassive black hole (SMBH) with its gaseous and stellar environment. Simulations to determine the accretion of stellar winds from the surrounding O-stars onto the black hole have been performed earlier, but in those the presence of the S-star system was ignored. The S-stars are a group of young massive B-stars in relatively close orbits around the black hole. Here we simulate those stars in order to study their contribution to the accretion rate, without taking the more distant and massive O-stars into account. We use the Astrophysical Multi-purpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics in a single simulation of the S-stars orbiting the supermassive black hole, and use this framework to determine the amount of gas that is accreted onto the black hole. We find that the accretion rate is sensitive to the wind properties of the S-stars (rate of mass-loss and terminal velocity). Our simulations are consistent with the observed accretion rate of the black hole only if the stars exhibit high wind massloss rates that are comparable with those of evolved 7-10 Myr old stars with masses of M=19-25 M_SUN. This is in contrast with observations that have shown that these stars are rather young, main-sequence B-stars. We therefore conclude that the S-stars cannot account for the accretion rate alone.
△ Less
Submitted 10 December, 2015;
originally announced December 2015.
-
The influence of dense gas rings on the dynamics of a stellar disk in the Galactic center
Authors:
Alessandro Alberto Trani,
Michela Mapelli,
Alessandro Bressan,
Federico Inti Pelupessy,
Arjen van Elteren,
Simon Portegies Zwart
Abstract:
The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermo…
▽ More
The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermost parsec. We investigate the perturbations exerted by a gas ring on a nearly-Keplerian stellar disk orbiting a SMBH by means of combined direct N-body and smoothed particle hydrodynamics simulations. We simulate the formation of gas rings through the infall and disruption of a molecular gas cloud, adopting different inclinations between the infalling gas cloud and the stellar disk. We find that a CNR-like ring is not efficient in affecting the stellar disk on a timescale of 3 Myr. In contrast, a gas ring in the innermost 0.5 pc induces precession of the longitude of the ascending node Omega, significantly affecting the stellar disk inclination. Furthermore, the combined effect of two-body relaxation and Omega-precession drives the stellar disk dismembering, displacing the stars from the disk. The impact of precession on the star orbits is stronger when the stellar disk and the inner gas ring are nearly coplanar. We speculate that the warm gas in the inner cavity might have played a major role in the evolution of the clockwise disk.
△ Less
Submitted 8 December, 2015;
originally announced December 2015.
-
Intermediate-mass black holes in globular clusters: observations and simulations
Authors:
Nora Lützgendorf,
Markus Kissler-Patig,
Karl Gebhardt,
Holger Baumgardt,
Diederik Kruijssen,
Eva Noyola,
Nadine Neumayer,
Tim de Zeeuw,
Anja Feldmeier,
Edwin van der Helm,
Inti Pelupessy,
Simon Portegies Zwart
Abstract:
The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. Our group investigated the presence of intermediate-mass black holes for a sample of 10 Galactic globular clusters. We measur…
▽ More
The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. Our group investigated the presence of intermediate-mass black holes for a sample of 10 Galactic globular clusters. We measured the inner kinematic profiles with integral-field spectroscopy and determined masses or upper limits of central black holes in each cluster. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M_bh - sigma) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). We ran N-body simulations of globular clusters containing IMBHs in a tidal field and studied their effects on mass-loss rates and remnant fractions and showed that an IMBH in the center prevents core collapse and ejects massive objects more rapidly. These simulations were further used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighbouring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole.
△ Less
Submitted 29 January, 2015;
originally announced January 2015.
-
Star Formation in the vicinity of Nuclear Black Holes: Young Stellar Objects close to Sgr A*
Authors:
B. Jalali,
F. I. Pelupessy,
A. Eckart,
S. Portegies Zwart,
N. Sabha,
A. Borkar,
J. Moultaka,
K. Mužić,
L. Moser
Abstract:
It is often assumed that the strong gravitational field of a super-massive black hole disrupts an adjacent molecular cloud preventing classical star formation in the deep potential well of the black hole. Yet, young stars have been observed across the entire nuclear star cluster of the Milky Way including the region close ($<$0.5~pc) to the central black hole, Sgr A*. Here, we focus particularly o…
▽ More
It is often assumed that the strong gravitational field of a super-massive black hole disrupts an adjacent molecular cloud preventing classical star formation in the deep potential well of the black hole. Yet, young stars have been observed across the entire nuclear star cluster of the Milky Way including the region close ($<$0.5~pc) to the central black hole, Sgr A*. Here, we focus particularly on small groups of young stars, such as IRS 13N located 0.1 pc away from Sgr A*, which is suggested to contain about five embedded massive young stellar objects ($<$1 Myr). We perform three dimensional hydrodynamical simulations to follow the evolution of molecular clumps orbiting about a $4\times10^6~M_{\odot}$ black hole, to constrain the formation and the physical conditions of such groups. The molecular clumps in our models assumed to be isothermal containing 100 $M_{\odot}$ in $<$0.2 pc radius. Such molecular clumps exist in the circumnuclear disk of the Galaxy. In our highly eccentrically orbiting clump, the strong orbital compression of the clump along the orbital radius vector and perpendicular to the orbital plane causes the gas densities to increase to values higher than the tidal density of Sgr A*, which are required for star formation. Additionally, we speculate that the infrared excess source G2/DSO approaching Sgr A* on a highly eccentric orbit could be associated with a dust enshrouded star that may have been formed recently through the mechanism supported by our models.
△ Less
Submitted 31 July, 2014;
originally announced August 2014.
-
A connected component-based method for efficiently integrating multiscale $N$-body systems
Authors:
Jürgen Jänes,
Federico I. Pelupessy,
Simon F. Portegies Zwart
Abstract:
We present a novel method for efficient direct integration of gravitational N-body systems with a large variation in characteristic time scales. The method is based on a recursive and adaptive partitioning of the system based on the connected components of the graph generated by the particle distribution combined with an interaction-specific time step criterion. It uses an explicit and approximate…
▽ More
We present a novel method for efficient direct integration of gravitational N-body systems with a large variation in characteristic time scales. The method is based on a recursive and adaptive partitioning of the system based on the connected components of the graph generated by the particle distribution combined with an interaction-specific time step criterion. It uses an explicit and approximately time-symmetric time step criterion, and conserves linear and angular momentum to machine precision. In numerical tests on astrophysically relevant setups, the method compares favourably to both alternative Hamiltonian-splitting integrators as well as recently developed block time step-based GPU-accelerated Hermite codes. Our reference implementation is incorporated in the HUAYNO code, which is freely available as a part of the AMUSE framework.
△ Less
Submitted 26 July, 2014;
originally announced July 2014.
-
Diagnostics of the molecular component of PDRs with mechanical heating. II: line intensities and ratios
Authors:
M. V. Kazandjian,
R. Meijerink,
I. Pelupessy,
F. P. Israel,
M. Spaans
Abstract:
CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an…
▽ More
CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates.
△ Less
Submitted 27 March, 2014;
originally announced March 2014.
-
Multi-scale and multi-domain computational astrophysics
Authors:
A. van Elteren,
F. I. Pelupessy,
S. Portegies Zwart
Abstract:
Astronomical phenomena are governed by processes on all spatial and temporal scales, ranging from days to the age of the Universe (13.8,Gyr) as well as from km size up to the size of the Universe. This enormous range in scales is contrived, but as long as there is a physical connection between the smallest and largest scales it is important to be able to resolve them all, and for the study of many…
▽ More
Astronomical phenomena are governed by processes on all spatial and temporal scales, ranging from days to the age of the Universe (13.8,Gyr) as well as from km size up to the size of the Universe. This enormous range in scales is contrived, but as long as there is a physical connection between the smallest and largest scales it is important to be able to resolve them all, and for the study of many astronomical phenomena this governance is present. Although covering all these scales is a challenge for numerical modelers, the most challenging aspect is the equally broad and complex range in physics, and the way in which these processes propagate through all scales. In our recent effort to cover all scales and all relevant physical processes on these scales we have designed the Astrophysics Multipurpose Software Environment (AMUSE). AMUSE is a Python-based framework with production quality community codes and provides a specialized environment to connect this plethora of solvers to a homogeneous problem solving environment.
△ Less
Submitted 19 February, 2014;
originally announced February 2014.
-
Young Stellar Objects close to Sgr A*
Authors:
B. Jalali,
I. Pelupessy,
A. Eckart,
S. Portegies Zwart,
N. Sabha,
A. Borkar,
J. Moultaka,
K. Muzic,
L. Moser
Abstract:
We aim at modelling small groups of young stars such as IRS 13N, 0.1 pc away from Sgr A*, which is suggested to contain a few embedded massive young stellar objects. We perform hydrodynamical simulations to follow the evolution of molecular clumps orbiting about a $4\times10^6 ~ M_{\odot}$ black hole, to constrain the formation and the physical conditions of such groups. We find that, the strong c…
▽ More
We aim at modelling small groups of young stars such as IRS 13N, 0.1 pc away from Sgr A*, which is suggested to contain a few embedded massive young stellar objects. We perform hydrodynamical simulations to follow the evolution of molecular clumps orbiting about a $4\times10^6 ~ M_{\odot}$ black hole, to constrain the formation and the physical conditions of such groups. We find that, the strong compression due to the black hole along the orbital radius vector of clumps evolving on highly eccentric orbits causes the clumps densities to increase to higher than the tidal density of Sgr A*, and required for star formation. This suggests that the tidal compression from the black hole could support star formation. Additionally, we speculate that the infrared excess source G2/DSO approaching Sgr A* on a highly eccentric orbit could be associated with a dust enshrouded star that may have been formed recently through the mechanism supported by our models.
△ Less
Submitted 19 November, 2013;
originally announced November 2013.
-
The Astrophysical Multipurpose Software Environment
Authors:
F. I. Pelupessy,
A. van Elteren,
N. de Vries,
S. L. W. McMillan,
N. Drost,
S. F. Portegies Zwart
Abstract:
We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guaran…
▽ More
We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guarantees that the whole framework is well-suited for distributed computation. It includes facilities for unit handling and data storage. Currently it includes codes for gravitational dynamics, stellar evolution, hydrodynamics and radiative transfer. Within each domain the interfaces to the codes are as similar as possible. We describe the design and implementation of AMUSE, as well as the main components and community codes currently supported and we discuss the code interactions facilitated by the framework. Additionally, we demonstrate how AMUSE can be used to resolve complex astrophysical problems by presenting example applications.
△ Less
Submitted 11 July, 2013;
originally announced July 2013.
-
The formation of planets in circumbinary disks
Authors:
F. I. Pelupessy,
S. Portegies Zwart
Abstract:
We examine the formation of planets around binary stars in light of the recently discovered systems Kepler 16, 34 and 35. We conduct hydrodynamical simulations of self gravitating disks around binary systems. The selected binary and disk parameters are chosen consistent with observed systems. The disks are evolved until they settle in a quasi-equilibrium and the resulting systems are compared with…
▽ More
We examine the formation of planets around binary stars in light of the recently discovered systems Kepler 16, 34 and 35. We conduct hydrodynamical simulations of self gravitating disks around binary systems. The selected binary and disk parameters are chosen consistent with observed systems. The disks are evolved until they settle in a quasi-equilibrium and the resulting systems are compared with the parameters of Kepler 16, 34 and 35. We find a close correspondence of the peak density at the inner disk gap and the orbit of the observed planets. We conclude, based on our simulations, that the orbits of the observed Kepler planets are determined by the size of the inner disk gap which for these systems results from the binary driving. This mediates planet formation either through the density enhancement or through planetary trapping at the density gradient inversion in the inner disk. For all three systems the current eccentricity of the planetary orbit is less than the disk eccentricity in the simulations. This, together with the long term stability of the orbits argues against in situ formation (e.g. a direct collapse scenario of the material in the ring). Conducting additional simulations of systems with a wider range of parameters (taken from a survey of eclipsing binaries), we find that the planet semi-major axis and binary eccentricity in such a scenario should be tightly correlated providing an observational test of this formation mechanism.
△ Less
Submitted 17 October, 2012;
originally announced October 2012.
-
N-body Integrators with Individual Time Steps from Hierarchical Splitting
Authors:
Federico I. Pelupessy,
Jürgen Jänes,
Simon F. Portegies Zwart
Abstract:
We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping, these integrators are momentum conserving and show excellent energy conservation in conjunction with a symmetrized time step criterion. We use an explicit but ap…
▽ More
We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping, these integrators are momentum conserving and show excellent energy conservation in conjunction with a symmetrized time step criterion. We use an explicit but approximate formula for the time symmetrization that is compatible with the use of individual time steps. No iterative scheme is necessary. We implement these ideas in the HUAYNO (available online at www.amusecode.org) code and present tests of the integrators and show that the presented integration schemes shows good energy conservation, with little or no systematic drift, while conserving momentum and angular momentum to machine precision for long term integrations.
△ Less
Submitted 25 May, 2012;
originally announced May 2012.
-
Multi-physics simulations using a hierarchical interchangeable software interface
Authors:
Simon Portegies Zwart,
Steve McMillan,
Arjen van Elteren,
Inti Pelupessy,
Nathan de Vries
Abstract:
We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery f…
▽ More
We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery from errors. The framework can be efficiently implemented and incurs an acceptable overhead. In practice, we measure the time spent in the framework to be less than 1% of the wall-clock time. Due to the unified structure of the interface, incorporating multiple modules addressing the same physics in different ways is relatively straightforward. Different modules may be advanced serially or in parallel. Despite initial concerns, we have encountered relatively few problems with this strict separation between modules, and the results of our simulations are consistent with earlier results using more traditional monolithic approaches. This framework provides a platform to combine existing simulation codes or develop new physical solver codes within a rich "ecosystem" of interchangeable modules.
△ Less
Submitted 16 September, 2012; v1 submitted 24 April, 2012;
originally announced April 2012.
-
Diagnostics of the Molecular Component of PDRs with Mechanical Heating
Authors:
M. V. Kazandjian,
R. Meijerink,
I. Pelupessy,
F. P. Israel,
M. Spaans
Abstract:
Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy.
Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the…
▽ More
Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy.
Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code.
Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude.
Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR.
△ Less
Submitted 22 March, 2012;
originally announced March 2012.
-
High-Performance Distributed Multi-Model / Multi-Kernel Simulations: A Case-Study in Jungle Computing
Authors:
Niels Drost,
Jason Maassen,
Maarten A. J. van Meersbergen,
Henri E. Bal,
F. Inti Pelupessy,
Simon Portegies Zwart,
Michael Kliphuis,
Henk A. Dijkstra,
Frank J. Seinstra
Abstract:
High-performance scientific applications require more and more compute power. The concurrent use of multiple distributed compute resources is vital for making scientific progress. The resulting distributed system, a so-called Jungle Computing System, is both highly heterogeneous and hierarchical, potentially consisting of grids, clouds, stand-alone machines, clusters, desktop grids, mobile devices…
▽ More
High-performance scientific applications require more and more compute power. The concurrent use of multiple distributed compute resources is vital for making scientific progress. The resulting distributed system, a so-called Jungle Computing System, is both highly heterogeneous and hierarchical, potentially consisting of grids, clouds, stand-alone machines, clusters, desktop grids, mobile devices, and supercomputers, possibly with accelerators such as GPUs.
One striking example of applications that can benefit greatly of Jungle Computing Systems are Multi-Model / Multi-Kernel simulations. In these simulations, multiple models, possibly implemented using different techniques and programming models, are coupled into a single simulation of a physical system. Examples include the domain of computational astrophysics and climate modeling.
In this paper we investigate the use of Jungle Computing Systems for such Multi-Model / Multi-Kernel simulations. We make use of the software developed in the Ibis project, which addresses many of the problems faced when running applications on Jungle Computing Systems. We create a prototype Jungle-aware version of AMUSE, an astrophysical simulation framework. We show preliminary experiments with the resulting system, using clusters, grids, stand-alone machines, and GPUs.
△ Less
Submitted 1 March, 2012;
originally announced March 2012.
-
Formation versus destruction: the evolution of the star cluster population in galaxy mergers
Authors:
J. M. Diederik Kruijssen,
F. Inti Pelupessy,
Henny J. G. L. M. Lamers,
Simon F. Portegies Zwart,
Nate Bastian,
Vincent Icke
Abstract:
(Abridged) Interacting galaxies are well-known for their high star formation rates and rich star cluster populations, but the rapidly changing tidal field can also efficiently destroy clusters. We use numerical simulations of merging disc galaxies to investigate which mechanism dominates. The simulations include a model for the formation and dynamical disruption of the entire star cluster populati…
▽ More
(Abridged) Interacting galaxies are well-known for their high star formation rates and rich star cluster populations, but the rapidly changing tidal field can also efficiently destroy clusters. We use numerical simulations of merging disc galaxies to investigate which mechanism dominates. The simulations include a model for the formation and dynamical disruption of the entire star cluster population. We find that the dynamical heating of clusters by tidal shocks is about an order of magnitude higher in interacting galaxies than in isolated galaxies. This is driven by the increased gas density, and is sufficient to destroy star clusters at a higher rate than new clusters are formed: the total number of clusters in the merger remnant is 2-50% of the amount in the progenitor discs, with low-mass clusters being disrupted preferentially. By adopting observationally motivated selection criteria, we find that the observed surplus of star clusters in nearby merging galaxies is caused by the bias to detect young, massive clusters. We provide a general expression for the survival fraction of clusters, which increases with the gas depletion time-scale. Due to the preferential disruption of low-mass clusters, the mass distribution of the surviving star clusters in a merger remnant develops a peak at a mass of about 10^3 Msun, which evolves to higher masses at a rate of 0.3-0.4 dex per Gyr. The peak mass initially depends weakly on the galactocentric radius, but this correlation disappears as the system ages. We discuss the similarities between the cluster populations of the simulated merger remnants and (young) globular cluster systems. Our results suggest that the combination of cluster formation and destruction should be widespread in the dense star-forming environments at high redshifts, which could provide a natural origin to present-day globular cluster systems.
△ Less
Submitted 5 December, 2011;
originally announced December 2011.
-
The evolution of embedded star clusters
Authors:
F. I. Pelupessy,
S. Portegies Zwart
Abstract:
We study the evolution of embedded clusters. The equations of motion of the stars in the cluster are solved by direct N-body integration while taking the effects of stellar evolution and the hydrodynamics of the natal gas content into account. The gravity of the stars and the surrounding gas are coupled self consistently to allow the realistic dynamical evolution of the cluster. While the equation…
▽ More
We study the evolution of embedded clusters. The equations of motion of the stars in the cluster are solved by direct N-body integration while taking the effects of stellar evolution and the hydrodynamics of the natal gas content into account. The gravity of the stars and the surrounding gas are coupled self consistently to allow the realistic dynamical evolution of the cluster. While the equations of motion are solved, a stellar evolution code keeps track of the changes in stellar mass, luminosity and radius. The gas liberated by the stellar winds and supernovae deposits mass and energy into the gas reservoir in which the cluster is embedded. We examine cluster models with 1000 stars, but we varied the star formation efficiency (between 0.05-0.5), cluster radius (0.1-1.0 parsec), the degree of virial support of the initial population of stars (0-100%) and the strength of the feedback. We find that an initial star fraction $M_\star/M_{\rm tot} > 0.05$ is necessary for cluster survival. Survival is more likely if gas is not blown out violently by a supernova and if the cluster has time to approach virial equilibrium during out-gassing. While the cluster is embedded, dynamical friction drives early and efficient mass segregation in the cluster. Stars of $m \gtrsim 2\,M_\odot$ are preferentially retained, at the cost of the loss of less massive stars. We conclude that the degree of mass segregation in open clusters such as the Pleiades is not the result of secular evolution but a remnant of its embedded stage.
△ Less
Submitted 3 November, 2011;
originally announced November 2011.
-
Multi-physics simulations using a hierarchical interchangeable software interface
Authors:
Simon Portegies Zwart,
Stephen McMillan,
Inti Pelupessy,
Arjen van Elteren
Abstract:
We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified software interface. Our framework is intrinsically parallel, and conveniently separates all components in memory. It performs unit conversion between different modules automatically and defines common data structures to communicate across different codes. We use the framework to…
▽ More
We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified software interface. Our framework is intrinsically parallel, and conveniently separates all components in memory. It performs unit conversion between different modules automatically and defines common data structures to communicate across different codes. We use the framework to simulate embedded star clusters. For this purpose we couple solvers for gravitational dynamics, stellar evolution and hydrodynamics to self consistently resolve the dynamical evolution simultaneousy with the internal nuclear evolution of the stars and the hydrodynamic response of the gas. We find, in contrast to earlier studies, that the survival of a young star cluster depends only weakly on the efficiency of star formation. The main reason for this weak dependency is the asymmetric expulsion of the embedding gas from the cluster.
△ Less
Submitted 12 October, 2011;
originally announced October 2011.
-
The escape of ionising radiation from high-redshift dwarf galaxies
Authors:
Jan-Pieter Paardekooper,
Federico I. Pelupessy,
Gabriel Altay,
Chael Kruip
Abstract:
The UV escape fraction from high-redshift galaxies plays a key role in models of cosmic reionisation. Because it is currently not possible to deduce the escape fractions during the epoch of reionisation from observations, we have to rely on numerical simulations. Our aim is to better constrain the escape fraction from high-redshift dwarf galaxies, as these are the most likely sources responsible f…
▽ More
The UV escape fraction from high-redshift galaxies plays a key role in models of cosmic reionisation. Because it is currently not possible to deduce the escape fractions during the epoch of reionisation from observations, we have to rely on numerical simulations. Our aim is to better constrain the escape fraction from high-redshift dwarf galaxies, as these are the most likely sources responsible for reionising the Universe. We employ a N-body/SPH method that includes realistic prescriptions for the physical processes that are important for the evolution of dwarf galaxies. These models are post-processed with radiative transfer to determine the escape fraction of ionising radiation. We perform a parameter study to assess the influence of the spin parameter, gas fraction and formation redshift of the galaxy and study the importance of numerical parameters as resolution, source distribution and local gas clearing. We find that the UV escape fraction from high-redshift dwarf galaxies that have formed a rotationally supported disc lie between 1e-5 and 0.1. The mass and angular momentum of the galaxy are the most important parameters that determine the escape fraction. We compare our results to previous work and discuss the uncertainties of our models. The low escape fraction we find for high-redshift dwarf galaxies is balanced by their high stellar content, resulting in an efficiency parameter for stars that is only marginally lower than the values found by semi-analytic models of reionisation. We therefore conclude that dwarf galaxies play an important role in cosmic reionisation also after the initial starburst phase, when the gas has settled into a disc.
△ Less
Submitted 18 April, 2011;
originally announced April 2011.
-
Modeling the formation and evolution of star cluster populations in galaxy simulations
Authors:
J. M. Diederik Kruijssen,
F. Inti Pelupessy,
Henny J. G. L. M. Lamers,
Simon F. Portegies Zwart,
Vincent Icke
Abstract:
(Abridged) The formation and evolution of star cluster populations are related to the galactic environment. Cluster formation is governed by processes acting on galactic scales, and star cluster disruption is driven by the tidal field. In this paper, we present a self-consistent model for the formation and evolution of star cluster populations, for which we combine an N-body/SPH galaxy evolution c…
▽ More
(Abridged) The formation and evolution of star cluster populations are related to the galactic environment. Cluster formation is governed by processes acting on galactic scales, and star cluster disruption is driven by the tidal field. In this paper, we present a self-consistent model for the formation and evolution of star cluster populations, for which we combine an N-body/SPH galaxy evolution code with semi-analytic models for star cluster evolution. The model includes star formation, feedback, stellar evolution, and star cluster disruption by two-body relaxation and tidal shocks. We apply the model by simulating a suite of 9 isolated disc galaxies and 24 galaxy mergers. The evolutionary histories of individual clusters in these simulations are discussed to illustrate how the environment of clusters changes in time and space. The resulting variability of the disruption rate with time and space affects the properties of star cluster populations. The combined effect of clusters escaping their dense formation sites (`cluster migration') and the preferential disruption of clusters residing in dense environments (`natural selection') implies that the mean disruption rate of the population decreases with cluster age. This affects the slope of the cluster age distribution, which becomes a function of the star formation rate density. The evolutionary histories of clusters in a galaxy merger vary widely and determine which clusters survive the merger. This impacts the age distributions and the locations of the surviving clusters at all times during a merger. We conclude that accounting for the interplay between the formation, disruption, and orbital histories of clusters enables a more sophisticated interpretation of observed properties of cluster populations, thereby extending the role of cluster populations as tracers of galaxy evolution.
△ Less
Submitted 4 February, 2011;
originally announced February 2011.
-
Deviations from the Schmidt-Kennicutt relations during early galaxy evolution
Authors:
Padelis P. Papadopoulos,
Federico I. Pelupessy
Abstract:
We utilize detailed time-varying models of the coupled evolution of stars and the HI, H_2, and CO-bright H_2 gas phases in galaxy-sized numerical simulations to explore the evolution of gas-rich and/or metal-poor systems, expected to be numerous in the Early Universe. The inclusion of the CO-bright H_2 gas phase, and the realistic rendering of star formation as an H_2-regulated process (and the n…
▽ More
We utilize detailed time-varying models of the coupled evolution of stars and the HI, H_2, and CO-bright H_2 gas phases in galaxy-sized numerical simulations to explore the evolution of gas-rich and/or metal-poor systems, expected to be numerous in the Early Universe. The inclusion of the CO-bright H_2 gas phase, and the realistic rendering of star formation as an H_2-regulated process (and the new feedback processes that this entails) allows the most realistic tracking of strongly evolving galaxies, and much better comparison with observations. We find that while galaxies eventually settle into states conforming to Schmidt-Kennicutt (S-K) relations, significant and systematic deviations of their star formation rates (SFRs) from the latter occur, especially pronounced and prolonged for ...
...This indicates potentially serious limitations of (S-K)-type relations as reliable sub-grid elements of star formation physics in simulations of structure formation in the Early Universe. We anticipate that galaxies with marked deviations from the S-K relations will be found at high redshifts as unbiased inventories of total gas mass become possible with ALMA and the EVLA.
△ Less
Submitted 19 May, 2010;
originally announced May 2010.
-
Molecular gas, CO, and star formation in galaxies: emergent empirical relations, feedback, and the evolution of very gas-rich systems
Authors:
Federico I. Pelupessy,
Padelis P. Papadopoulos
Abstract:
We use time-varying models of the coupled evolution of the HI, H_2 gas phases and stars in galaxy-sized numerical simulations to: a) test for the emergence of the Kennicutt-Schmidt (K-S) and the H_2-pressure relation, b) explore a realistic H_2-regulated star formation recipe which brings forth a neglected and potentially significant SF-regulating factor, and c) go beyond typical galactic enviro…
▽ More
We use time-varying models of the coupled evolution of the HI, H_2 gas phases and stars in galaxy-sized numerical simulations to: a) test for the emergence of the Kennicutt-Schmidt (K-S) and the H_2-pressure relation, b) explore a realistic H_2-regulated star formation recipe which brings forth a neglected and potentially significant SF-regulating factor, and c) go beyond typical galactic environments (for which these galactic empirical relations are deduced) to explore the early evolution of very gas-rich galaxies. In this work we model low mass galaxies ($M_{\rm baryon} \le 10^9 \msun$), while incorporating an independent treatment of CO formation and destruction, the most important tracer molecule of H2 in galaxies, along with that for the H2 gas itself. We find that both the K-S and the H_2-pressure empirical relations can robustly emerge in galaxies after a dynamic equilibrium sets in between the various ISM states, the stellar component and its feedback. (abridged)
△ Less
Submitted 3 November, 2009; v1 submitted 26 October, 2009;
originally announced October 2009.
-
Simulations of the Sunyaev-Zeldovich Effect from Quasars
Authors:
Suchetana Chatterjee,
Tiziana Di Matteo,
Arthur Kosowsky,
Inti Pelupessy
Abstract:
Quasar feedback has most likely a substantial but only partially understood impact on the formation of structure in the universe. A potential direct probe of this feedback mechanism is the Sunyaev-Zeldovich effect: energy emitted from quasar heats the surrounding intergalactic medium and induce a distortion in the microwave background radiation passing through the region. Here we examine the for…
▽ More
Quasar feedback has most likely a substantial but only partially understood impact on the formation of structure in the universe. A potential direct probe of this feedback mechanism is the Sunyaev-Zeldovich effect: energy emitted from quasar heats the surrounding intergalactic medium and induce a distortion in the microwave background radiation passing through the region. Here we examine the formation of such hot quasar bubbles using a cosmological hydrodynamic simulation which includes a self-consistent treatment of black hole growth and associated feedback, along with radiative gas cooling and star formation. From this simulation, we construct microwave maps of the resulting Sunyaev-Zeldovich effect around black holes with a range of masses and redshifts. The size of the temperature distortion scales approximately with black hole mass and accretion rate, with a typical amplitude up to a few micro-Kelvin on angular scales around 10 arcseconds. We discuss prospects for the direct detection of this signal with current and future single-dish and interferometric observations, including ALMA and CCAT. These measurements will be challenging, but will allow us to characterize the evolution and growth of supermassive black holes and the role of their energy feedback on galaxy formation.
△ Less
Submitted 2 August, 2008;
originally announced August 2008.
-
SPHRAY: A Smoothed Particle Hydrodynamics Ray Tracer for Radiative Transfer
Authors:
Gabriel Altay,
Rupert A. C. Croft,
Inti Pelupessy
Abstract:
We introduce SPHRAY, a Smoothed Particle Hydrodynamics (SPH) ray tracer designed to solve the 3D, time dependent, radiative transfer (RT) equations for arbitrary density fields. The SPH nature of SPHRAY makes the incorporation of separate hydrodynamics and gravity solvers very natural. SPHRAY relies on a Monte Carlo (MC) ray tracing scheme that does not interpolate the SPH particles onto a grid…
▽ More
We introduce SPHRAY, a Smoothed Particle Hydrodynamics (SPH) ray tracer designed to solve the 3D, time dependent, radiative transfer (RT) equations for arbitrary density fields. The SPH nature of SPHRAY makes the incorporation of separate hydrodynamics and gravity solvers very natural. SPHRAY relies on a Monte Carlo (MC) ray tracing scheme that does not interpolate the SPH particles onto a grid but instead integrates directly through the SPH kernels. Given initial conditions and a description of the sources of ionizing radiation, the code will calculate the non-equilibrium ionization state (HI, HII, HeI, HeII, HeIII, e) and temperature (internal energy/entropy) of each SPH particle. The sources of radiation can include point like objects, diffuse recombination radiation, and a background field from outside the computational volume. The MC ray tracing implementation allows for the quick introduction of new physics and is parallelization friendly. A quick Axis Aligned Bounding Box (AABB) test taken from computer graphics applications allows for the acceleration of the raytracing component. We present the algorithms used in SPHRAY and verify the code by performing all the test problems detailed in the recent Radiative Transfer Comparison Project of Iliev et. al. The Fortran 90 source code for SPHRAY and example SPH density fields are made available on a companion website (www.sphray.org).
△ Less
Submitted 25 February, 2008;
originally announced February 2008.
-
How rapidly do supermassive black hole "seeds" grow at early times?
Authors:
Federico I. Pelupessy,
Tiziana Di Matteo,
Benedetta Ciardi
Abstract:
We investigate the physical conditions for the growth of intermediate mass seed black holes assumed to have formed from remnants of the first generation of massive stars. We follow the collapse of high-sigma halos with Tvir > 1e4 K using cosmological, smooth-particle hydrodynamic (SPH) simulations in the standard LCDM model. During collapse of the parent halo the seed holes are incorporated thro…
▽ More
We investigate the physical conditions for the growth of intermediate mass seed black holes assumed to have formed from remnants of the first generation of massive stars. We follow the collapse of high-sigma halos with Tvir > 1e4 K using cosmological, smooth-particle hydrodynamic (SPH) simulations in the standard LCDM model. During collapse of the parent halo the seed holes are incorporated through mergers into larger systems and accrete mass from the surrounding gas. We include a self-consistent treatment of star formation, black hole accretion and associated feedback processes. Even under optimistic assumptions for the seed black hole mass and for efficient merger rates, we find that seed holes in halos M<1e10 Msun never reach the conditions for critical Eddington growth. Most of the black hole growth in this regime is determined by the initial mass and the merger rates. Critical accretion rates are reached, albeit only after a significant delay, at the time of collapse z~7) for 3-4 sigma halos of M~1e11 Msun. Our results imply M_BH = 5e6 Msun (M_halo/1e11 Msun)^0.78 at the time of collapse. The required conditions of Eddington growth to explain the build-up of supermassive black holes (~1e9 Msun), as implied by Sloan quasars at z>6, are therefore hard to meet in such a scenario. Without a 'jump-start' these conditions may be only achieved in extremely rare halos with M_halo > 1e13 Msun that collapsed before z~6. The sub-Eddington regime in which black holes holes accrete at early time implies a small contribution to the reionization by miniquasar but still sufficient to cause appreciable heating of the IGM at z<15-18.
△ Less
Submitted 29 March, 2007;
originally announced March 2007.
-
Incorporating the molecular gas phase in galaxy-size numerical simulations: first applications in dwarf galaxies
Authors:
F. I. Pelupessy,
P. P. Papadopoulos,
P. P. van der Werf
Abstract:
We present models of the evolution of the gaseous and stellar content of galaxies incorporating the formation of H_2 out of HI gas as part of such a model. We do so by formulating a subgrid model for gas clouds that uses well-known cloud scaling relations and solves for the HI-H_2 balance set by the H_2 formation on dust grains and its FUV-induced photodissociation by the temporally and spatiall…
▽ More
We present models of the evolution of the gaseous and stellar content of galaxies incorporating the formation of H_2 out of HI gas as part of such a model. We do so by formulating a subgrid model for gas clouds that uses well-known cloud scaling relations and solves for the HI-H_2 balance set by the H_2 formation on dust grains and its FUV-induced photodissociation by the temporally and spatially varying interstellar radiation field. This allows the seamless tracking of the evolution of the H_2 gas phase, its precursor Cold Neutral Medium (CNM) HI gas, simultaneously with the star formation. Our most important findings are: a) a significant dependence of the HI-H_2 transition and the resultant H_2 gas mass on the ambient metallicity and the H_2 formation rate, b) the important influence of the characteristic star formation timescale (regulating the ambient FUV radiation field) on the equilibrium H_2 gas mass and c) the possibility of a diffuse H_2 gas phase. Finally, we implement and briefly explore a novel approach of using the ambient H_2 gas mass fraction as a criterion for the onset of star formation. (abridged)
△ Less
Submitted 21 March, 2006;
originally announced March 2006.
-
Periodic bursts of Star Formation in Irregular Galaxies
Authors:
F. I. Pelupessy,
P. P. van der Werf,
V. Icke
Abstract:
We present N-body/SPH simulations of the evolution of an isolated dwarf galaxy including a detailed model for the ISM, star formation and stellar feedback. Depending on the strength of the feedback, the modelled dwarf galaxy shows periodic or quasi-periodic bursts of star formation of moderate strength. The period of the variations is related to the dynamical timescale, of the order of…
▽ More
We present N-body/SPH simulations of the evolution of an isolated dwarf galaxy including a detailed model for the ISM, star formation and stellar feedback. Depending on the strength of the feedback, the modelled dwarf galaxy shows periodic or quasi-periodic bursts of star formation of moderate strength. The period of the variations is related to the dynamical timescale, of the order of $1.5~10^8$ yr. We show that the results of these simulations are in good agreement with recent detailed observations of dwarf irregulars (dIrr) and that the peculiar kinematic and morphological properties of these objects,as revealed by high resolution HI studies, are fully reproduced. We discuss these results in the context of recent surveys of dwarf galaxies and point out that if the star formation pattern of our model galaxy is typical for dwarf irregulars this could explain the scatter of observed properties of dwarf galaxies. Specifically, we show that the time sampled distribution of the ratio between the instanteneous star formation rate (SFR) and the mean SFR is similar to that distribution in observed sample of dwarf galaxies.
△ Less
Submitted 7 April, 2004;
originally announced April 2004.
-
Density Estimators in Particle Hydrodynamics: DTFE versus regular SPH
Authors:
F. I. Pelupessy,
W. E. Schaap,
R. van de Weygaert
Abstract:
We present the results of a study confronting density maps reconstructed by the Delaunay Tessellation Field Estimator (DTFE) and by regular SPH kernel-based techniques. The comparison between the two methods clearly demonstrates the superior performance of the DTFE with respect to conventional SPH methods, in particular at locations where SPH appears to fail. The DTFE is a fully self-adaptive te…
▽ More
We present the results of a study confronting density maps reconstructed by the Delaunay Tessellation Field Estimator (DTFE) and by regular SPH kernel-based techniques. The comparison between the two methods clearly demonstrates the superior performance of the DTFE with respect to conventional SPH methods, in particular at locations where SPH appears to fail. The DTFE is a fully self-adaptive technique for reconstructing continuous density fields from discrete particle distributions, and is based upon the corresponding Delaunay tessellation. Its principal asset is its complete independence of arbitrary smoothing functions and parameters specifying the properties of these. As a result it manages to faithfully reproduce the anisotropies of the local particle distribution and through its adaptive and local nature proves to be optimally suited for uncovering the full structural richness in the density distribution. Through the improvement in local density estimates, calculations invoking the DTFE will yield a much better representation of physical processes which depend on density. The presented results form an encouraging step towards the application and insertion of the DTFE in astrophysical hydrocodes. We describe an outline for the construction of a particle hydrodynamics code in which the DTFE replaces kernel-based methods. Further discussion addresses the issue and possibilities for a moving grid based hydrocode invoking the DTFE, and Delaunay tessellations, in an attempt to combine the virtues of the Eulerian and Lagrangian approaches.
△ Less
Submitted 6 March, 2003; v1 submitted 4 March, 2003;
originally announced March 2003.
-
The radiation driven winds of rotating B[e] supergiants
Authors:
Inti Pelupessy,
Henny J. G. L. M. Lamers,
Jorick S. Vink
Abstract:
We have formulated the momentum equation for sectorial line driven winds from rotating stars including: (a) the oblateness of the star, (b) gravity darkening (von Zeipel effect), (c) conservation of angular momentum, (d) line driving specified by the force multiplier parameters (k, alpha, delta), (e) finite disk correction factors for an oblate star with gravity darkening for both the continuum…
▽ More
We have formulated the momentum equation for sectorial line driven winds from rotating stars including: (a) the oblateness of the star, (b) gravity darkening (von Zeipel effect), (c) conservation of angular momentum, (d) line driving specified by the force multiplier parameters (k, alpha, delta), (e) finite disk correction factors for an oblate star with gravity darkening for both the continuum and the line driving. The equations are solved numerically. We calculated the distribution of the mass flux and the wind velocity from the pole to the equator for the winds of B[e]-supergiants. Rotation decreases the terminal velocity in the equatorial region but hardly affects the wind velocity from the poles; it enhances the mass flux from the poles while the mass flux from the equator remains nearly the same. These effects increase with increasing rotation rates.
We also calculated models with a bi-stability jump around 25 000 K, using force multipliers recently calculated with a Monte Carlo technique. In this case the mass flux increases drastically from the pole to the equator and the terminal velocity decreases drastically from pole to equator. This produces a density contrast in the wind rho(equator/pole) of about a factor 10 independent of the rotation rate of the star. We suggest that the observed density contrast of a factor 100 of the disks of B[e] stars may be reached by taking into account the wind compression due to the trajectories of the gas above the critical point towards the equatorial plane.
△ Less
Submitted 15 May, 2000;
originally announced May 2000.