We present stellar_wind.py, a module that provides multiple methods of simulating stellar winds using smoothed particle hydrodynamics codes (SPH) within the astrophysical multipurpose software environment (AMUSE) framework. With the simple wind mode, we create SPH wind particles in a spherically symmetric shell. We inject the wind particles with a velocity equal to their terminal velocity. The accelerating wind mode is similar, but with this method particles can be injected with a lower initial velocity than the terminal velocity and they are accelerated away from the star according to an acceleration function. With the heating wind mode, SPH particles are created with zero initial velocity with respect to the star, but instead wind particles are given an internal energy based on the integrated mechanical luminosity of the star. This mode is designed to be used on longer timescales and larger spatial scales compared to the other two modes and assumes that the star is embedded in a gas cloud. For fast winds, we find that both the simple and accelerating mode can reproduce the desired velocity, density and temperature profiles. For slow winds, the simple wind mode is insufficient due to dominant hydrodynamical effects that change the wind velocities. The accelerating mode, with additional options to account for these hydrodynamical effects, can still reproduce the desired wind profiles. We test the heating mode by simulating both a normal wind and a supernova explosion of a single star in a uniform density medium. The stellar wind simulation results matches the analytical solution for an expanding wind bubble. The supernova simulation gives qualitatively correct results, but the simulated bubble expands faster than the analytical solution predicts. We conclude with an example of a triple star system which includes the colliding winds of all three stars.

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]

We model the intermediate-mass black hole HLX-1, using the Hubble Space Telescope, XMM-Newton and Swift. We quantify the relative contributions of a bluer component, function of X-ray irradiation, and a redder component, constant and likely coming from an old stellar population. We estimate a black hole mass of about (2^{+2}_{-1}) x 10^4 M_{sun}, a spin parameter a/M ~ 0.9 for moderately face-on view, and a peak outburst luminosity of about 0.3 times the Eddington luminosity. We discuss the discrepancy between the characteristic sizes inferred from the short X-ray timescale (R ~ a few 10^{11} cm) and from the optical emitter (R sqrt[cos theta] ~ 2.2 x 10^{13} cm). One possibility is that the optical emitter is a circumbinary disk; however, we disfavour this scenario because it would require a very small donor star. A more plausible scenario is that the disk is large but only the inner annuli are involved in the X-ray outburst. We propose that the recurrent outbursts are caused by an accretion-rate oscillation driven by wind instability in the inner disk. We argue that the system has a long-term-average accretion rate of a few percent Eddington, just below the upper limit of the low/hard state; a wind-driven oscillation can trigger transitions to the high/soft state, with a recurrence period of ~1 year (much longer than the binary period, which we estimate as ~10 days). The oscillation that dominated the system in the last decade is now damped such that the accretion rate no longer reaches the level required to trigger a transition. Finally, we highlight similarities between disk winds in HLX-1 and in the Galactic black hole V404 Cyg.

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.

The X-ray source HLX-1 near the spiral galaxy ESO 243-49 is currently the best intermediate-mass black hole candidate. It has a peak bolometric luminosity of $10^{42}$ erg s$^{-1}$, which implies a mass inflow rate of $\sim10^{-4}$ MSun yr$^{-1}$, but the origin of this mass is unknown. It has been proposed that there is a star on an eccentric orbit around the black hole which transfers mass at pericentre. To investigate the orbital evolution of this system, we perform stellar evolution simulations using mesa and SPH simulations of a stellar orbit around an intermediate-mass black hole using fi. We run and couple these simulations using the amuse framework. We find that mass is lost through both the first and second Lagrange points and that there is a delay of up to 10 days between the pericentre passage and the peak mass loss event. The orbital evolution timescales we find in our simulations are larger than what is predicted by analytical models, but these models fall within the errors of our results. Despite the fast orbital evolution, we are unable to reproduce the observed change in outburst period. We conclude that the change in the stellar orbit with the system parameters investigated here is unable to account for all observed features of HLX-1.

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.