by Nissanke, Samaya and Vallisneri, Michele and Nelemans, Gijs and Prince, Thomas A.
17 pages, 3 figures, 5 tables, to be submitted to the Astrophysical Journal

Compact Galactic binaries where at least one member is a white dwarf (WD) or neutron star constitute the majority of individually detectable sources for future low-frequency space-based gravitational-wave (GW) observatories; in addition, they form an unresolved continuum, the dominant Galactic foreground at frequencies below a few mHz. A handful of ultra-compact binaries, observed at optical, ultraviolet and X-ray wavelengths, are known verification sources for space-based GW interferometers. Due to the paucity of electromagnetic observations, the majority of studies of Galactic-binary populations so far have been based on population-synthesis simulations. However, recent surveys have reported several new detections of compact binaries including double WDs, providing new constraints for population estimates. In this article, we evaluate the impact of revised local densities of interacting WD binaries on future low-frequency GW observations. Specifically: we consider five scenarios that explain these densities with different assumptions on the formation of interacting systems; we simulate corresponding populations of detached and interacting WD binaries; we estimate the number of individually detectable GW sources and the magnitude of the confusion-noise foreground, in the case of two GW interferometers with armlengths of 1 and 5 Mkm. We confirm earlier estimates of thousands of detached-binary detections, but project only a few ten to a few hundred detections of interacting systems. We also confirm estimates for the confusion-noise foreground (except in one scenario that explains smaller local densities of interacting systems with fewer progenitor detached systems). Last, we provide a general scaling argument that shows that the magnitude of the GW foreground can be derived robustly from the merger rate of Galactic WD binaries, and depends only weakly on the structure of the Galaxy.

by Garcia-Senz, Domingo and Cabezon, Ruben M. and Escartin, Jose Antonio
15 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

In this paper we develop and check a fully conservative SPH scheme based on a tensor formulation which can be applied to simulate astrophysical systems. In the proposed scheme derivatives are calculated from an integral expression which leads to a tensor, rather than vectorial, estimation of gradients and reduces to the standard formulation in the continuum limit. The new formulation improves the interpolation of physical magnitudes, leading to a set of conservative equations which looks similar to those of standard SPH. The resulting scheme was checked using a variety of well known tests, all of them simulated in two dimensions. An application of the proposed tensor method to astrophysics was also discussed by simulating the stability of a sun-like polytrope calculated in three dimensions.

by Rice, W. K. M. and Forgan, D. H. and Armitage, P. J.
9 pages, 11 figures, MNRAS in press

Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.

by Benacquista, Matthew J. and Downing, Jonathan M. B.
88 pages, 13 figures. Submitted update of Living Reviews article

Galactic globular clusters are old, dense star systems typically containing 10\super{4}–10\super{7} stars. As an old population of stars, globular clusters contain many collapsed and degenerate objects. As a dense population of stars, globular clusters are the scene of many interesting close dynamical interactions between stars. These dynamical interactions can alter the evolution of individual stars and can produce tight binary systems containing one or two compact objects. In this review, we discuss theoretical models of globular cluster evolution and binary evolution, techniques for simulating this evolution that leads to relativistic binaries, and current and possible future observational evidence for this population. Our discussion of globular cluster evolution will focus on the processes that boost the production of hard binary systems and the subsequent interaction of these binaries that can alter the properties of both bodies and can lead to exotic objects. Direct {\it N}-body integrations and Fokker–Planck simulations of the evolution of globular clusters that incorporate tidal interactions and lead to predictions of relativistic binary populations are also discussed. We discuss the current observational evidence for cataclysmic variables, millisecond pulsars, and low-mass X-ray binaries as well as possible future detection of relativistic binaries with gravitational radiation.

by Benacquista, Matthew
5 pages, submitted to ApJL

The recently discovered J0651+2844 is a detached, eclipsing white dwarf binary with an orbital period of 765 s. We investigate the prospects for the detection of gravitational radiation from this system and estimate the effect of the tidal deformation of the low-mass component on the period evolution of the system. Because of the high inclination of the system, the amplitude of the gravitational waves at Earth will be as much as a factor of two lower than that from an optimally oriented system. The dominant contribution of tidal corrections to the period evolution comes from the increase in rotational energy of the components as they spin up to remain tied to the orbital period. This contribution results in an advance of the timing of the eclipses by an additional 0.3 s after one year.

by Cordero-Carrión, Isabel and Cerdá-Durán, Pablo and Ibáñez, José María
18 pages, 9 figures, 2 tables

The Fully Constrained Formulation (FCF) of General Relativity is a novel framework introduced as an alternative to the hyperbolic formulations traditionally used in numerical relativity. The FCF equations form a hybrid elliptic-hyperbolic system of equations including explicitly the constraints. We present an implicit-explicit numerical algorithm to solve the hyperbolic part, whereas the elliptic sector shares the form and properties with the well known Conformally Flat Condition (CFC) approximation. We show the stability andconvergence properties of the numerical scheme with numerical simulations of vacuum solutions. We have performed the first numerical evolutions of the coupled system of hydrodynamics and Einstein equations within FCF. As a proof of principle of the viability of the formalism, we present 2D axisymmetric simulations of an oscillating neutron star. In order to simplify the analysis we have neglected the back-reaction of the gravitational waves into the dynamics, which is small (<2 %) for the system considered in this work. We use spherical coordinates grids which are well adapted for simulations of starsand allow for extended grids that marginally reach the wave zone. We have extracted the gravitational wave signature and compared to the Newtonian quadrupole and hexadecapole formulae. Both extraction methods show agreement within the numerical errors and the approximations used (~5 %).

by Rice, W. K. M. and Armitage, P. J. and Mamatsashvili, G. R. and Lodato, G. and Clarke, C. J.
MNRAS, in press

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress – which can be represented as an effective viscous alpha – out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

by Kiuchi, Kenta and Shibata, Masaru and Montero, Pedro J. and Font, José A.
4 pages, 4 figure, to be published in PRL

Black hole (BH)–torus systems are promising candidates for the central engine of gamma-ray bursts (GRBs), and also possible outcomes of the collapse of supermassive stars to supermassive black holes (SMBHs). By three-dimensional general relativistic numerical simulations, we show that an $latex m=1$ nonaxisymmetric instability grows for a wide range of self-gravitating tori orbiting BHs. The resulting nonaxisymmetric structure persists for a timescale much longer than the dynamical one, becoming a strong emitter of large amplitude, quasiperiodic gravitational waves. Our results indicate that both, the central engine of GRBs and newly formed SMBHs, can be strong gravitational wave sources observable by forthcoming ground-based and spacecraft detectors.

by Valsecchi, Francesca and Farr, Will M. and Willems, Bart and Deloye, Christopher . J. and Kalogera, Vassiliki
21 pages, 7 figures, submitted to APJ

Galactic interacting double white dwarfs (DWD) are guaranteed gravitational wave (GW) sources for the GW detector LISA, with more than 10^4 binaries expected to be detected over the mission’s lifetime. While the majority of DWDs are expected to be circular, dynamical interactions in globular clusters can lead to a sub-population of eccentric DWDs detectable by LISA. Here we investigate the potential for constraining the white dwarf (WD) properties through apsidal precession in these binaries. We analyze the tidal, rotational, and general relativistic contributions to apsidal precession by using detailed He WD models, where the evolution of the star’s interior is followed throughout the cooling phase. In agreement with previous studies of zero-temperature WDs, we find that apsidal precession in eccentric DWDs can lead to a detectable shift in the emitted GW signal when binaries with cool (old) components are considered. This shift increases significantly for hot (young) WDs. We find that apsidal motion in hot (cool) DWDs is dominated by tides at orbital frequencies above ~10^(-4) Hz [10^(-3) Hz]. The analysis of apsidal precession in these sources while ignoring the tidal component would lead to an extreme bias in the mass determination, and could lead us to misidentify WDs as neutron stars or black holes. We use the detailed WD models to show that for older, cold WDs, there is a unique relationship that ties the radius and apsidal precession constant to the WD masses, therefore allowing tides to be used as a tool to constrain the source masses.

by Lu, Ting-Ni and Kong, Albert K. H.
Accepted by ApJL

We present results of deep radio continuum observations of two galactic globular clusters 47 Tucanae (47 Tuc) and Omega Centauri (Omega Cen) with Australia Telescope Compact Array (ATCA). No statistically significant evidence for radio emission was found from the central region for the two clusters. However, both clusters show a 2.5-sigma detection near the center that may be confirmed by future deeper radio observations. The 3-sigma upper limits of the radio observations is 20 and 40 microJy/beam for Omega Cen and 47 Tuc, respectively. By using the fundamental plane of accreting black holes which describes the relationship between radio luminosity, X-ray luminosity and black hole mass, we constrain the mass of a possible intermediate-mass black hole (IMBH) in the globualar clusters. We also compare our results with other globular clusters and discuss the existence of IMBHs in globular clusters.

by Marsh, T. R.
11 pages, 2 figures, in press Classical and Quantum Gravity. Summary of talk given at 8th LISA Symposium, Stanford, July 2010

Close pairs of white dwarfs are potential progenitors of Type~Ia supernovae and they are common, with of order 100 — 300 million in the Galaxy. As such they will be significant, probably dominant, sources of the gravitational waves detectable by LISA. In the context of LISA’s goals for fundamental physics, double white dwarfs are a source of noise, but from an astrophysical perspective, they are of considerable interest in their own right. In this paper I discuss our current knowledge of double white dwarfs and their close relatives (and possible descendants) the AM~CVn stars. LISA will add to our knowledge of these systems by providing the following unique constraints: (i) an almost direct measurement of the Galactic merger rate of DWDs from the detection of short period systems and their period evolution, (ii) an accurate and precise normalisation of binary evolution models at the shortest periods, (iii) a determination of the evolutionary pathways to the formation of AM~CVn stars, (iv) measurements of the influence of tidal coupling in white dwarfs and its significance for stabilising mass transfer, and (v) discovery of numerous examples of eclipsing white dwarfs with the potential for optical follow-up to test models of white dwarfs.

by Lodato, Giuseppe and Clarke, Cathie C.
7 pages, 2 figures. Accepted by MNRAS

Stimulated by recent results by Meru and Bate (2010a,b), we revisit the issue of resolution requirements for simulating self-gravitating accretion discs with Smoothed Particle Hydrodynamics (SPH). We show that the results by Meru and Bate (2010a) are consistent with those of Meru and Bate (2010b) if they are both interpreted as driven by resolution effects, therefore implying that the resolution criterion for cooling gaseous discs is a function of the imposed cooling rate. We discuss two possible numerical origins of such dependence, which are both consistent with the limited number of available data. Our results tentatively indicate that convergence for current simulations is being reached for a number of SPH particles approaching 10 millions (for a disc mass of order 10 per cent of the central object mass), which would set the critical cooling time for fragmentation at about $latex 15\Omega^{-1}$, roughly a factor two larger than previously thought. More in general, we discuss the extent to which the large number of recent numerical results are reliable or not. We argue that those results that pertain to the dynamics associated with gravitational instabilities (such as the locality of angular momentum transport, and the relationship between density perturbation and induced stress) are robust, while those pertaining to the thermodynamics of the system (such as the determination of the critical cooling time for fragmentation) can be affected by poor resolution.

by Ferreira, Barbara T.
PhD thesis, 199 pages, 32 figures. For higher quality images, please contact the author

This thesis investigates phenomena occurring in black-hole accretion discs which are likely to induce high-frequency quasi-periodic variability. Two classes of pseudo-relativistic theoretical models are studied. The first is based on the stability of transonic accretion flows and its connection to a disc instability that takes the form of propagating waves (viscous overstability). The second class of models looks at accretion-disc oscillations which are trapped due to the non-monotonic variation of the epicyclic frequency in relativistic flows. In particular, it focuses on inertial waves trapped below the maximum of the epicyclic frequency which are excited in deformed, warped or eccentric, discs. The influence of a transonic background on the propagation of these inertial modes is also investigated.

by Moeckel, Nickolas and Clarke, Cathie J.
8 pages, submitted to MNRAS. Includes revisions per referee’s comments

We investigate the contraction of accreting protoclusters using an extension of n-body techniques that incorporates the accretional growth of stars from the gaseous reservoir in which they are embedded. Following on from Monte Carlo studies by Davis et al., we target our experiments toward populous clusters likely to experience collisions as a result of accretion-driven contraction. We verify that in less extreme star forming environments, similar to Orion, the stellar density is low enough that collisions are unimportant, but that conditions suitable for stellar collisions are much more easily satisfied in large-n clusters, i.e. n ~ 30,000 (we argue, however, that the density of the Arches cluster is insufficient for us to expect stellar collisions to have occurred in the cluster’s prior evolution). We find that the character of the collision process is not such that it is a route toward smoothly filling the top end of the mass spectrum. Instead, runaway growth of one or two extreme objects can occur within less than 1 Myr after accretion is shut off, resulting in a few objects with masses several times the maximum reached by accretion. The rapid formation of these objects is due to not just the post-formation dynamical evolution of the clusters, but an interplay of dynamics and the accretional growth of the stars. We find that accretion-driven cluster shrinkage results in a distribution of gas and stars that offsets the disruptive effect of gas expulsion, and we propose that the process can lead to massive binaries and early mass segregation in star clusters.

by Lousto, Carlos O. and Zlochower, Yosef
4 pages, 4 figures, 3 tables

We perform the first fully nonlinear numerical simulations of black-hole binaries with mass ratios 100:1. Our technique for evolving such extreme mass ratios is based on the moving puncture approach with a new gauge condition and an optimal choice of the mesh refinement (plus large computational resources). We achieve a convergent set of results for simulations starting with a small nonspinning black hole just outside the ISCO that then performs over two orbits before plunging into the 100 times more massive black hole. We compute the gravitational energy and momenta radiated as well as the final remnant parameters and compare these quantities with the corresponding perturbative estimates. The results show a close agreement. We briefly discuss the relevance of this simulations for Advanced LIGO, third-generation ground based detectors, and LISA observations, and self-force computations.

by Colaiuda, A. and Ferrari, V. and Gualtieri, L. and Pons, J. A.
25 pages, 9 figures, submitted to MNRAS

We find numerical solutions of the coupled system of Einstein-Maxwell’s equations with a linear approach, in which the magnetic field acts as a perturbation of a spherical neutron star. In our study, magnetic fields having both poloidal and toroidal components are considered, and higher order multipoles are also included. We evaluate the deformations induced by different field configurations, paying special attention to those for which the star has a prolate shape. We also explore the dependence of the stellar deformation on the particular choice of the equation of state and on the mass of the star. Our results show that, for neutron stars with mass M = 1.4 Msun and surface magnetic fields of the order of 10^15 G, a quadrupole ellipticity of the order of 10^(-6) – 10^(-5) should be expected. Low mass neutron stars are in principle subject to larger deformations (quadrupole ellipticities up to 10^(-3) in the most extreme case). The effect of quadrupolar magnetic fields is comparable to that of dipolar components. A magnetic field permeating the whole star is normally needed to obtain negative quadrupole ellipticities, while fields confined to the crust typically produce positive quadrupole ellipticities.

by Lesur, G. and Ogilvie, G. I.
6 pages, 5 figures, accepted in MNRAS

The mechanism of angular momentum transport in accretion discs has long been debated. Although the magnetorotational instability appears to be a promising process, poorly ionized regions of accretion discs may not undergo this instability. In this letter, we revisit the possibility of transporting angular momentum by turbulent thermal convection. Using high-resolution spectral methods, we show that strongly turbulent convection can drive outward angular momentum transport at a rate that is, under certain conditions, compatible with observations of discs. We find however that the angular momentum transport is always much weaker than the vertical heat transport. These results indicate that convection might be another way to explain global disc evolution, provided that a sufficiently unstable vertical temperature profile can be maintained.

by Zwart, Simon Portegies and McMillan, Steve and Gieles, Mark
Only 88 pages. To be published in ARAA. Final version to be submitted on Friday 12 February

Young massive clusters are dense aggregates of young stars that form the fundamental building blocks of galaxies. Several examples exist in the Milky Way Galaxy and the Local Group, but they are particularly abundant in starburst and interacting galaxies. The few young massive clusters that are close enough to resolve are of prime interest for studying the stellar mass function and the ecological interplay between stellar evolution and stellar dynamics. The distant unresolved clusters may be effectively used to study the star-cluster mass function, and they provide excellent constraints on the formation mechanisms of young cluster populations. Young massive clusters are expected to be the nurseries for many unusual objects, including a wide range of exotic stars and binaries. So far only a few such objects have been found in young massive clusters, although their older cousins, the globular clusters, are unusually rich in stellar exotica. In this review we focus on star clusters younger than $latex \sim100$ Myr, more than a few current crossing times old, and more massive than $latex \sim10^4$ \Msun, irrespective of cluster size or environment. We describe the global properties of the currently known young massive star clusters in the Local Group and beyond, and discuss the state of the art in observations and dynamical modeling of these systems. In order to make this review readable by observers, theorists, and computational astrophysicists, we also review the cross-disciplinary terminology.

by Lattanzi, Massimiliano and Montani, Giovanni
6 pages, 4 figures. To appear in Europhysics Letters

We provide a new analysis of the system of partial differential equations describing the radial and vertical equilibria of the plasma in accretion disks. In particular, we show that the partial differential system can be separated once a definite, oscillatory (or hyperbolic) form for the radial dependence of the relevant physical quantities is assumed. The system is thus reduced to an ordinary differential system in the vertical dimensionless coordinate. The resulting equations can be integrated analytically in the limit of small magnetic pressure. We complete our analysis with a direct numerical integration of the more general case. The main result is that a ring-like density profile (i.e., radial oscillations in the mass density) can appear even in the limit of small magnetic pressure.

by Downing, J. M. B. and Benacquista, M. J. and Giersz, M. and Spurzem, R.
18 pages, 5 Tables, 10 Figures, Submitted to MNRAS

We study the compact binary population in star clusters, focusing on binaries containing neutron stars and black holes, using a self-consistent Monte Carlo treatment of dynamics and full stellar evolution. We find that the black holes experience strong mass segregation and become centrally concentrated. In the core the black holes interact strongly with each other and black hole-black hole binaries are formed very efficiently. The strong interactions, however, also destroy or eject the black hole-black hole binaries. We find no black hole-black hole mergers within our simulations but produce many hard escapers that will merge in the galactic field within a Hubble time. We also find two highly eccentric black hole-black hole binaries that are potential LISA sources, suggesting that star clusters are interesting targets for space-based detectors. We conclude that star clusters must be taken into account when predicting compact binary population statistics.

by Gvaramadze, V. V. and Gualandris, A. and Zwart, S. Portegies
4 pages, 2 figure, to appear in Star Clusters — Basic Galactic Building Blocks throughout Time and Space, Proceed. of the IAU Symp. 266, eds. R. de Grijs and J. Lepine

We performed numerical simulations of dynamical encounters between hard massive binaries and a very massive star (VMS; formed through runaway mergers of ordinary stars in the dense core of a young massive star cluster), in order to explore the hypothesis that this dynamical process could be responsible for the origin of high-velocity (\geq 200-400 km/s) early or late B-type stars. We estimated the typical velocities produced in encounters between very tight massive binaries and VMSs (of mass of \geq 200 Msun) and found that about 3-4 per cent of all encounters produce velocities of \geq 400 km/s, while in about 2 per cent of encounters the escapers attain velocities exceeding the Milky Ways’s escape velocity. We therefore argue that the origin of high-velocity (\geq 200-400 km/s) runaway stars and at least some so-called hypervelocity stars could be associated with dynamical encounters between the tightest massive binaries and VMSs formed in the cores of star clusters. We also simulated dynamical encounters between tight massive binaries and single ordinary 50-100 Msun stars. We found that from 1 to \simeq 4 per cent of these encounters can produce runaway stars with velocities of \geq 300-400 km/s (typical of the bound population of high-velocity halo B-type stars) and occasionally (in less than 1 per cent of encounters) produce hypervelocity (\geq 700 km/s) late B-type escapers.

by Stroeer, Alexander and Nelemans, Gijs
5 pages, 3 figures, accepted for publication in MNRAS Letters

We study the effect of short term variations of the evolution of AM CVn systems on their gravitational wave emissions and in particular LISA observations. We model the systems according to their equilibrium mass-transfer evolution as driven by gravitational wave emission and tidal interaction, and determine their reaction to a sudden perturbation of the system. This is inspired by the suggestion to explain the orbital period evolution of the ultra-compact binary systems V407 Vul and RX-J0806+1527 by non-equilibrium mass transfer. The characteristics of the emitted gravitational wave signal are deduced from a Taylor expansion of a Newtonian quadrupolar emission model, and the changes in signal structure as visible to the LISA mission are determined. We show that short term variations can significantly change the higher order terms in the expansion, and thus lead to spurious (non) detection of frequency derivatives. This may hamper the estimation of the parameters of the system, in particular their masses and distances. However, we find that overall detection is still secured as signals still can be described by general templates. We conclude that a better modelling of the effects of short term variations is needed to prepare the community for astrophysical evaluations of real gravitational wave data of AM CVn systems.

by Beuther, H. and Walsh, A. J. and Longmore, S. N.
21 pages, 32 figures, accepted for ApJS. A high-resolution version can be found at http://www.mpia.de/homes/beuther/papers.html

To better understand the physical properties of accretion disks in high-mass star formation, we present a study of a 12 high-mass accretion disk candidates observed at high spatial resolution with the Australia Telescope Compact Array (ATCA) in the NH3 (4,4) and (5,5) lines. Almost all sources were detected in NH3, directly associated with CH3OH Class II maser emission. From the remaining eleven sources, six show clear signatures of rotation and/or infall motions. These signatures vary from velocity gradients perpendicular to the outflows, to infall signatures in absorption against ultracompact HII regions, to more spherical infall signatures in emission. Although our spatial resolution is ~1000AU, we do not find clear Keplerian signatures in any of the sources. Furthermore, we also do not find flattened structures. In contrast to this, in several of the sources with rotational signatures, the spatial structure is approximately spherical with sizes exceeding 10^4 AU, showing considerable clumpy sub-structure at even smaller scales. This implies that on average typical Keplerian accretion disks — if they exist as expected — should be confined to regions usually smaller than 1000AU. It is likely that these disks are fed by the larger-scale rotating envelope structure we observe here. Furthermore, we do detect 1.25cm continuum emission in most fields of view.