BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Giuseppe Lodato

Giuseppe Lodato appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.HEastro-ph.GAastro-ph.COgr-qcastro-ph.IM

Catalog footprint

What is connected

51works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2025arXiv

Brightness variability in polar circumbinary disks

In binary systems with a strongly misaligned disk, the central binary stars can travel a significant vertical distance above and below the disk's orbital plane. This can cause large changes in illumination of the disk over the course of the binary orbital period. We use both analytic and radiative transfer models to examine the effect of changes in stellar illumination on the appearance of the disk, particularly in the case of the polar disk HD 98800B. We find that the observed flux from the disk can vary significantly over the binary orbital period, producing a periodically varying lightcurve which peaks twice each binary orbit. The amount of flux variation is strongly influenced by the disk geometry. We suggest that these flux variations produce several observable signatures, and that these observables may provide constraints on different properties of the disk such as its vertical structure, geometry, and cooling rate.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Probing initial distributions of orbital eccentricity and disc misalignment via polar discs

In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the distributions of initial eccentricity and mutual inclination in multiple stellar systems at birth.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accretion rates in hierarchical triple systems with discs

Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical configurations. Despite its astrophysical relevance, differential accretion in hierarchical systems remains to be understood. In this work, we investigate whether the accretion trends expected in binaries are valid for higher order multiples. We performed a set of 3D Smoothed Particle Hydrodynamics simulations of binaries and of hierarchical triples (HTs) embedded in an accretion disc, with the code Phantom. We identify for the first time accretion trends in HTs and their deviations compared to binaries. These deviations, due to the interaction of the small binary with the infalling material from the circum-triple disc, can be described with a semi-analytical prescription. Generally, the smaller binary of a HT accretes more mass than a single star of the same mass as the smaller binary. We found that in a HT, if the small binary is heavier than the third body, the standard differential accretion scenario (whereby the secondary accretes more of the mass) is hampered. Reciprocally, if the small binary is lighter than the third body, the standard differential accretion scenario is enhanced. The peculiar differential accretion mechanism we find in HTs is expected to affect their mass ratio distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

External or internal companion exciting the spiral arms in CQ Tau?

We present new high-contrast images in near-infrared wavelengths ($λ$ = 1.04, 1.24, 1.62, 2.18 and 3.78$μ$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M$_{\rm Jup}$ immediately outside the spiral arms at $\sim$60 au and above 2-3 M$_{\rm Jup}$ beyond 100 au to 5$σ$ confidence assuming hot-start models. We do, however, detect four spiral arms, a double-arc and evidence for shadows in scattered light cast by a misaligned inner disc. Our observations may be explained by an unseen close-in companion on an inclined and eccentric orbit. Such a hypothesis would also account for the disc CO cavity and disturbed kinematics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extreme Mass Ratio Inspirals triggered by Massive Black Hole Binaries: from Relativistic Dynamics to Cosmological Rates

Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass-ratio $q=m/M$ in the range $~10^{-9}-10^{-4}$ and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the central massive black hole, a number of alternative formation channels has been proposed, including binary tidal break-up, migration in accretion disks and secular and chaotic dynamics around a massive black hole binary (MBHB). In this work, we take an extensive closer look at this latter scenario, investigating how EMRIs can be triggered by a MBHBs, formed in the aftermath of galaxy mergers. By employing a suite of relativistic three-body simulations, we evaluate the efficiency of EMRI formation for different parameters of the MBHB, assessing the importance of both secular and chaotic dynamics. By modelling the distribution of compact objects in galaxy nuclei, we estimate the resulting EMRI formation rate, finding that EMRI are produced in a sharp burst, with peak rates that are 10-100 times higher than the standard two-body relaxation channel, lasting for 10$^6$--10$^8$ years. By coupling our results with an estimate of the cosmic MBHB merger rate, we finally forecast that LISA could observe ${\cal O}(10)$ EMRIs per year formed by this channel.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Investigating Protoplanetary Disk Cooling through Kinematics: Analytical GI Wiggle

It is likely that young protostellar disks undergo a self-gravitating phase. Such systems are characterized by the presence of a spiral pattern that can be either in a quasi-steady state or in a nonlinear unstable condition. This spiral wave affects both the gas dynamics and kinematics, resulting in deviations from the Keplerian rotation. Recently, a lot of attention has been devoted to kinematic studies of planet-forming environments, and we are now able to measure even small perturbations of velocity field ($\lesssim 1 \%$ of the Keplerian speed) thanks to high spatial and spectral resolution observations of protostellar disks. In this work, we investigate the kinematic signatures of gravitational instability: we perform an analytical study of the linear response of a self-gravitating disk to a spiral-like perturbation, focusing our attention on the velocity field perturbations. We show that unstable disks have clear kinematic imprints into the gas component across the entire disk extent, due to the GI spiral wave perturbation, resulting in deviations from Keplerian rotation. The shape of these signatures depends on several parameters, but they are significantly affected by the cooling factor: by detecting these features, we can put constraints on protoplanetary disk cooling.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mapping the Planetary Wake in HD 163296 with Kinematics

We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD~163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk major axis confirm that the wake induces predominantly radial motion, as predicted by models of planet-disk interaction. These results provide the first direct confirmation of planet wakes generated by Lindblad resonances. Mapping the wake provides a constraint on the disk aspect ratio, which is required to measure the mass of the planet.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the time evolution of the $M_{\rm d} - M_\star$ and $\dot M - M_\star$ correlations for protoplanetary discs: the viscous timescale increases with stellar mass

Large surveys of star-forming regions have unveiled power-law correlations between the stellar mass and the disc parameters, such as the disc mass $M_{\mathrm{d}} \propto {M_{\star}}^{λ_{\mathrm{m}}}$ and the accretion rate $\dot M \propto {M_{\star}}^{λ_{\mathrm{acc}}}$. The observed slopes appear to be increasing with time, but the reason behind the establishment of these correlations and their subsequent evolution is still uncertain. We conduct a theoretical analysis of the impact of viscous evolution on power-law initial conditions for a population of protoplanetary discs. We find that, for evolved populations, viscous evolution enforces the two correlations to have the same slope, $λ_{\mathrm{m}}$ = $λ_{\mathrm{acc}}$, and that this limit is uniquely determined by the initial slopes $λ_{\mathrm{m}, 0}$ and $λ_{\mathrm{acc}, 0}$. We recover the increasing trend claimed from the observations when the difference in the initial values, $δ_0 = λ_{\mathrm{m}, 0} - λ_{\mathrm{acc}, 0}$, is larger than $1/2$; moreover, we find that this increasing trend is a consequence of a positive correlation between the viscous timescale and the stellar mass. We also present the results of disc population synthesis numerical simulations, that allow us to introduce a spread and analyse the effect of sampling, which show a good agreement with our analytical predictions. Finally, we perform a preliminary comparison of our numerical results with observational data, which allows us to constrain the parameter space of the initial conditions to $λ_{\mathrm{m}, 0} \in [1.2, 2.1]$, $λ_{\mathrm{acc}, 0} \in [0.7, 1.5]$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Physics of Accretion Discs, Winds And Jets in Tidal Disruption Events

Accretion onto black holes is an efficient mechanism in converting the gas mass-energy into energetic outputs as radiation, wind and jet. Tidal disruption events, in which stars are tidally torn apart and then accreted onto supermassive black holes, offer unique opportunities of studying the accretion physics as well as the wind and jet launching physics across different accretion regimes. In this review, we systematically describe and discuss the models that have been developed to study the accretion flows and jets in tidal disruption events. A good knowledge of these physics is not only needed for understanding the emissions of the observed events, but also crucial for probing the general relativistic space-time around black holes and the demographics of supermassive black holes via tidal disruption events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discovery of a low-mass companion embedded in the disk of the young massive star MWC 297 with VLT/SPHERE

We report the discovery of a low-mass stellar companion around the young Herbig Be star MWC 297. We performed multi-epoch high-contrast imaging in the near infrared (NIR) with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. The companion is found at projected separation of 244.7$\pm$13.2 au and a position angle of 176.4$\pm$0.1 deg. The large separation supports formation via gravitational instability. From the spectrum, we estimate a mass of 0.1-0.5 M$_{\odot}$, the range conveying uncertainties in the extinction of the companion and in evolutionary models at young ages. The orbit coincides with a gap in the dust disk inferred from the Spectral Energy Distribution (SED). The young age ($\lesssim$ 1 Myr) and mass ratio with the central star ($\sim 0.01$) makes the companion comparable to PDS 70~b, suggesting a relation between formation scenarios and disk dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dual-Wavelength ALMA Observations of Dust Rings in Protoplanetary Disks

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($α_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($τ\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient Dust Ring Formation in Misaligned Circumbinary Discs

Binary systems exert a gravitational torque on misaligned discs orbiting them, causing differential precession which may produce disc warping and tearing. While this is well understood for gas-only discs, misaligned cirumbinary discs of gas and dust have not been thoroughly investigated. We perform SPH simulations of misaligned gas and dust discs around binaries to investigate the different evolution of these two components. We choose two different disc aspect ratios: A thin case for which the gas disc always breaks, and a thick one where a smooth warp develops throughout the disc. For each case, we run simulations of five different dust species with different degrees of coupling with the gas component, varying in Stokes number from 0.002 (strongly coupled dust) to 1000 (effectively decoupled dust). We report two new phenomena: First, large dust grains in thick discs pile up at the warp location, forming narrow dust rings, due to a difference in precession between the gas and dust components. These pile ups do not form at gas pressure maxima, and hence are different from conventional dust traps. This effect is most evident for St ~ 10 - 100. Second, thin discs tear and break only in the gas, while dust particles with St > 10 form a dense dust trap due to the steep pressure gradient caused by the break in the gas. We find that dust with St < 0.02 closely follow the gas particles, for both thin and thick discs, with radial drift becoming noticeable only for the largest grains in this range.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Future Simulations of Tidal Disruption Events

Tidal disruption events involve numerous physical processes (fluid dynamics, magnetohydrodynamics, radiation transport, self-gravity, general relativistic dynamics) in highly nonlinear ways, and, because TDEs are transients by definition, frequently in non-equilibrium states. For these reasons, numerical solution of the relevant equations can be an essential tool for studying these events. In this chapter, we present a summary of the key problems of the field for which simulations offer the greatest promise and identify the capabilities required to make progress on them. We then discuss what has been---and what cannot be---done with existing numerical methods. We close with an overview of what methods now under development may do to expand our ability to understand these events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HST survey of the Orion Nebula Cluster in the H$_2$O 1.4 $μ$m absorption band: I. A census of substellar and planetary mass objects

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 μ$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Is the gap in the DS Tau disc hiding a planet?

Recent mm-wavelength surveys performed with the Atacama Large Millimeter Array (ALMA) have revealed protoplanetary discs characterized by rings and gaps. A possible explanation for the origin of such rings is the tidal interaction with an unseen planetary companion. The protoplanetary disc around DS Tau shows a wide gap in the ALMA observation at 1.3 mm. We construct a hydrodynamical model for the dust continuum observed by ALMA assuming the observed gap is carved by a planet between one and five Jupiter masses. We fit the shape of the radial intensity profile along the disc major axis varying the planet mass, the dust disc mass, and the evolution time of the system. The best fitting model is obtained for a planet with $M_{\rm p}=3.5\,M_{\rm Jup}$ and a disc with $M_{\rm dust}= 9.6\cdot10^{-5}\,M_{\odot}$. Starting from this result, we also compute the expected signature of the planet in the gas kinematics, as traced by CO emission. We find that such a signature (in the form of a "kink" in the channel maps) could be observed by ALMA with a velocity resolution between $0.2-0.5\,\rm{kms}^{-1}$ and a beam size between 30 and 50 mas.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Orbital and mass constraints of the young binary system IRAS 16293-2422 A

We present 3 mm ALMA continuum and line observations at resolutions of 6.5 au and 13 au respectively, toward the Class 0 system IRAS 16293-2422 A. The continuum observations reveal two compact sources towards IRAS 16293-2422 A, coinciding with compact ionized gas emission previously observed at radio wavelengths (A1 and A2), confirming the long-known radio sources as protostellar. The emission towards A2 is resolved and traces a dust disk with a FWHM size of ~12 au, while the emission towards A1 sets a limit to the FWHM size of the dust disk of ~4 au. We also detect spatially resolved molecular kinematic tracers near the protostellar disks. Several lines of the J=5-4 rotational transition of HNCO, NH2CHO and t-HCOOH are detected, with which we derived individual line-of-sight velocities. Using these together with the CS (J=2-1), we fit Keplerian profiles towards the individual compact sources and derive masses of the central protostars. The kinematic analysis indicates that A1 and A2 are a bound binary system. Using this new context for the previous 30 years of VLA observations, we fit orbital parameters to the relative motion between A1 and A2 and find the combined protostellar mass derived from the orbit is consistent with the masses derived from the gas kinematics. Both estimations indicate masses consistently higher (0.5< M1<M2<2 Msun) than previous estimations using lower resolution observations of the gas kinematics. The ALMA high-resolution data provides a unique insight into the gas kinematics and masses of a young deeply embedded bound binary system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simulations of Tidal Disruption Events

Numerical simulations have historically played a major role in understanding the hydrodynamics of the tidal disruption process. Given the complexity of the geometry of the system, the challenges posed by the problem have indeed stimulated much work on the numerical side. Smoothed Particles Hydrodynamics methods, for example, have seen their very first applications in the context of tidal disruption and still play a major role to this day. Likewise, initial attempts at simulating the evolution of the disrupted star with the so-called affine method have been historically very useful. In this Chapter, we provide an overview of the numerical techniques used in the field and of their limitations, and summarize the work that has been done to simulate numerically the tidal disruption process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The gravitational wave background signal from tidal disruption events

In this paper we derive the gravitational wave stochastic background from tidal disruption events (TDEs). We focus on both the signal emitted by main sequence stars disrupted by super-massive black holes (SMBHs) in galaxy nuclei, and on that from disruptions of white dwarfs by intermediate mass black holes (IMBHs) located in globular clusters. We show that the characteristic strain $h_{\rm c}$'s dependence on frequency is shaped by the pericenter distribution of events within the tidal radius, and under standard assumptions $h_{\rm c} \propto f^{-1/2}$. This is because the TDE signal is a burst of gravitational waves at the orbital frequency of the closest approach. In addition, we compare the background characteristic strains with the sensitivity curves of the upcoming generation of space-based gravitational wave interferometers: the Laser Interferometer Space Antenna (LISA), TianQin, ALIA, the DECI-hertz inteferometer Gravitational wave Observatory (DECIGO) and the Big Bang Observer (BBO). We find that the background produced by main sequence stars might be just detected by BBO in its lowest frequency coverage, but it is too weak for all the other instruments. On the other hand, the background signal from TDEs with white dwarfs will be within reach of ALIA, and especially of DECIGO and BBO, while it is below the LISA and TianQin sensitive curves. This background signal detection will not only provide evidence for the existence of IMBHs up to redshift $z\sim 3$, but it will also inform us on the number of globular clusters per galaxy and on the occupation fraction of IMBHs in these environments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

What causes the fragmentation of debris streams in TDEs?

A tidal disruption event (TDE) occurs when a star passes too close to a supermassive black hole and gets torn apart by its gravitational tidal field. After the disruption, the stellar debris form an expanding gaseous stream. The morphology and evolution of this stream is particularly interesting as it ultimately determines the observational properties of the event itself. In this work we perform 3D hydrodynamical simulations of the TDE of a star modelled as a polytropic sphere of index γ = 5/3, and study the gravitational stability of the resulting gas stream. We provide an analytical solution for the evolution of the stream in the bound, unbound and marginally bound case, that allows us to describe the stream properties and analyse the time-scales of the physical processes involved, applying a formalism developed in star formation context. Our results are that, when fragmentation occurs, it is fueled by the failure of pressure in supporting the gas against its self-gravity. We also show that a stability criterion that includes also the stream gas pressure proves to be far more accurate than one that only considers the black hole tidal forces, giving analytical predictions of the time evolution of the various forces associated to the stream. Our results point out that fragmentation occurs on timescales longer compared with the observational windows of these events and is thus not expected to give rise to significant observational features.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Effects of photoevaporation on protoplanetary disc `isochrones'

Protoplanetary discs are the site of star and planet formation, and their evolution and consequent dispersal deeply affect the formation of planetary systems. In the standard scenario they evolve on timescales ~Myr due to the viscous transport of angular momentum. The analytical self-similar solution for their evolution predicts also specific disc isochrones in the accretion rate - disc mass plane. However, photoevaporation by radiation emitted by the central star is likely to dominate the gas disc dispersal of the innermost region, introducing another (shorter) timescale for this process. In this paper, we include the effect of internal (X and EUV) photoevaporation on the disc evolution, finding numerical solutions for a population of protoplanetary discs. Our models naturally reproduce the expected quick dispersal of the inner region of discs when their accretion rates match the rate of photoevaporative mass loss, in line with previous studies. We find that photoevaporation preferentially removes the lightest discs in the sample. The net result is that, counter-intuitively, photoevaporation increases the average disc mass in the sample, by dispersing the lightest discs. At the same time, photoevaporation also reduces the mass accretion rate by cutting the supply of material from the outer to the inner disc. In a purely viscous framework, this would be interpreted as the result of a longer viscous evolution, leading to an overestimate of the disc age. Our results thus show that photoevaporation is a necessary ingredient to include when interpreting observations of large disc samples with measured mass accretion rates and disc masses. Photoevaporation leaves a characteristic imprint on the shape of the isochrone. Accurate data in the accretion rate - disc mass plane in the low disc mass region therefore give clues on the typical photoevaporation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Long lived dust rings around HD169142

Recent ALMA observations of the protoplanetary disc around HD~169142 reveal a peculiar structure made of concentric dusty rings: a main ring at $\sim$20 au, a triple system of rings at $\sim 55-75$ au in millimetric continuum emission and a perturbed gas surface density from the $^{12}$CO,$^{13}$CO and C$^{18}$O$(J=2-1)$ surface brightness profile. In this Letter, we perform three-dimensional numerical simulations and radiative transfer modeling exploring the possibility that two giant planets interacting with the disc and orbiting in resonant locking can be responsible for the origin of the observed dust inner rings structure. We find that in this configuration the dust structure is actually long lived while the gas mass of the disc is accreted onto the star and the giant planets, emptying the inner region. In addition, we also find that the innermost planet is located at the inner edge of the dust ring, and can accrete mass from the disc, generating a signature in the dust ring shape that can be observed in mm ALMA observations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Type II migration strikes back -- An old paradigm for planet migration in discs

In this paper we analyse giant gap-opening planet migration in protoplanetary discs, focusing on the type II migration regime. According to standard type II theory, planets migrate at the same rate as the gas in the disc, as they are coupled to the disc viscous evolution; however, recent studies questioned this paradigm, suggesting that planets migrate faster than the disc material. We study the problem through 2D long-time simulations of systems consistent with type II regime, using the hydrodynamical grid code FARGO3D. Even though our simulations confirm the presence of an initial phase characterised by fast migration, they also reveal that the migration velocity slows down and eventually reaches the theoretical prediction if we allow the system to evolve for enough time. We find the same tendency to evolve towards the theoretical predictions at later times when we analyse the mass flow through the gap and the torques acting on the planet. This transient is related to the initial conditions of our (and previous) simulations, and is due to the fact that the shape of the gap has to adjust to a new profile, once the planet is set into motion. Secondly, we test whether the type II theory expectation that giant planet migration is driven by viscosity is consistent with our simulation by comparing simulations with the same viscosity and different disc mass (or viceversa). We find a good agreement with the theory, since when the discs are characterised by the same viscosity, the migration properties are the same.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bad prospects for the detection of giant stars' tidal disruption: effect of the ambient medium on bound debris

Most massive galaxies are thought to contain a supermassive black hole in their centre surrounded by a tenuous gas environment, leading to no significant emission. In these quiescent galaxies, tidal disruption events represent a powerful detection method for the central black hole. Following the disruption, the stellar debris evolves into an elongated gas stream, which partly falls back towards the disruption site and accretes onto the black hole producing a luminous flare. Using an analytical treatment, we investigate the interaction between the debris stream and the gas environment of quiescent galaxies. Although we find dynamical effects to be negligible, we demonstrate that Kelvin-Helmholtz instability can lead to the dissolution of the stream into the ambient medium before it reaches the black hole, likely dimming the associated flare. This result is robust against the presence of a typical stellar magnetic field and fast cooling within the stream. Furthermore, we find this effect to be enhanced for disruptions involving more massive black holes and/or giant stars. Consequently, although disruptions of evolved stars have been proposed as a useful probe of black holes with masses $\gtrsim 10^8 \, M_{\odot}$, we argue that the associated flares are likely less luminous than expected.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gas squeezing during the merger of a supermassive black hole binary

We study accretion rates during the gravitational wave-driven merger of a binary supermassive black hole embedded in an accretion disc, formed by gas driven to the centre of the galaxy. We use 3D simulations performed with PHANTOM, a Smoothed Particle Hydrodynamics code. Contrary to previous investigations, we show that there is evidence of a "squeezing phenomenon", caused by the compression of the inner disc gas when the secondary black hole spirals towards the primary. This causes an increase in the accretion rates that always exceed the Eddington rate. We have studied the main features of the phenomenon for a mass ratio $q = 10^{-3}$ between the black holes, including the effects of numerical resolution, the secondary accretion radius and the disc thickness. With our disc model with a low aspect ratio, we show that the mass expelled from the orbit of the secondary is negligible ($< 5\%$ of the initial disc mass), different to the findings of previous 2D simulations with thicker discs. The increase in the accretion rates in the last stages of the merger leads to an increase in luminosity, making it possible to detect an electromagnetic precursor of the gravitational wave signal emitted by the coalescence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Long-term stream evolution in tidal disruption events

A large number of tidal disruption event (TDE) candidates have been observed recently, often differing in their observational features. Two classes appear to stand out: X-ray and optical TDEs, the latter featuring lower effective temperatures and luminosities. These differences can be explained if the radiation detected from the two categories of events originates from different locations. In practice, this location is set by the evolution of the debris stream around the black hole and by the energy dissipation associated with it. In this paper, we build an analytical model for the stream evolution, whose dynamics is determined by both magnetic stresses and shocks. Without magnetic stresses, the stream always circularizes. The ratio of the circularization timescale to the initial stream period is $t_{\rm ev}/t_{\rm min} = 8.3 (M_{\rm h}/10^6 M_{\odot})^{-5/3} β^{-3}$, where $M_{\rm h}$ is the black hole mass and $β$ is the penetration factor. If magnetic stresses are strong, they can lead to the stream ballistic accretion. The boundary between circularization and ballistic accretion corresponds to a critical magnetic stresses efficiency $v_{\rm A}/v_{\rm c} \approx 10^{-1}$, largely independent of $M_{\rm h}$ and $β$. However, the main effect of magnetic stresses is to accelerate the stream evolution by strengthening self-crossing shocks. Ballistic accretion therefore necessarily occurs on the stream dynamical timescale. The shock luminosity associated to energy dissipation is sub-Eddington but decays as $t^{-5/3}$ only for a slow stream evolution. Finally, we find that the stream thickness rapidly increases if the stream is unable to cool completely efficiently. A likely outcome is its fast evolution into a thick torus, or even an envelope completely surrounding the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the origin of horseshoes in transitional discs

We investigate whether the rings, lopsided features and horseshoes observed at millimetre wavelengths in transitional discs can be explained by the dynamics of gas and dust at the edge of the cavity in circumbinary discs. We use 3D dusty smoothed particle hydrodynamics calculations to show that binaries with mass ratio $q \gtrsim 0.04$ drive eccentricity in the central cavity, naturally leading to a crescent-like feature in the gas density, which is accentuated in the mm dust grain population with intensity contrasts in mm-continuum emission of 10 or higher. We perform mock observations to demonstrate that these features closely match those observed by ALMA, suggesting that the origin of rings, dust horseshoes and other non-axisymmetric structures in transition discs can be explained by the presence of massive companions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Suppression of the accretion rate in thin discs around binary black holes

We present three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations investigating the dependence of the accretion rate on the disc thickness around an equal-mass, circular black hole binary system. We find that for thick/hot discs, with $H/R\gtrsim 0.1$, the binary torque does not prevent the gas from penetrating the cavity formed in the disc by the binary (in line with previous investigations). The situation drastically changes for thinner discs, in this case the mass accretion rate is suppressed, such that only a fraction (linearly dependent on $H/R$) of the available gas is able to flow within the cavity and accrete on to the binary. Extrapolating this result to the cold and thin accretion discs expected around supermassive black hole binary systems implies that this kind of systems accretes less material than predicted so far, with consequences not only for the electromagnetic and gravitational waves emissions during the late inspiral phase but also for the recoil speed of the black hole formed after binary coalescence, thus influencing also the evolutionary path both of the binary and of the host galaxy. Our results, being scale-free, are also applicable to equal mass, circular binaries of stellar mass black holes, such as the progenitor of the recently discovered gravitational wave source GW150914.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two mechanisms for dust gap opening in protoplanetary discs

We identify two distinct physical mechanisms for dust gap opening by embedded planets in protoplanetary discs based on the symmetry of the drag-induced motion around the planet: I) A mechanism where low mass planets, that do not disturb the gas, open gaps in dust by tidal torques assisted by drag in the inner disc, but resisted by drag in the outer disc; and II) The usual, drag assisted, mechanism where higher mass planets create pressure maxima in the gas disc which the drag torque then acts to evacuate further in the dust. The first mechanism produces gaps in dust but not gas, while the second produces partial or total gas gaps which are deeper in the dust phase. Dust gaps do not necessarily indicate gas gaps.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

"Spoon-feeding" an AGN

Tidal disruption events (TDEs) occur when a star, passing too close to a massive black hole, is ripped apart by tidal forces. A less dramatic event occurs if the star orbits just outside the tidal radius, resulting in a mild stripping of mass. Thus, if a star orbits a central black hole on one of these bound eccentric orbits, weaker outbursts will occur recurring every orbital period. Thanks to five Swift observations, we observed a recent flare from the close by (92 Mpc) galaxy IC 3599, where a possible TDE was already observed in December 1990 during the Rosat All-Sky Survey. By light curve modeling and spectral fitting, we account for all these events as the non-disruptive tidal stripping of a single star into a 9.5 yr highly eccentric bound orbit. This is the first example of periodic partial tidal disruptions, possibly spoon-feeding the central black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disc formation from tidal disruptions of stars on eccentric orbits by Schwarzschild black holes

The potential of tidal disruption of stars to probe otherwise quiescent supermassive black holes cannot be exploited, if their dynamics is not fully understood. So far, the observational appearance of these events has been derived from analytical extrapolations of the debris dynamical properties just after disruption. By means of hydrodynamical simulations, we investigate the subsequent fallback of the stream of debris towards the black hole for stars already bound to the black hole on eccentric orbits. We demonstrate that the debris circularize due to relativistic apsidal precession which causes the stream to self-cross. The circularization timescale varies between 1 and 10 times the period of the star, being shorter for more eccentric and/or deeper encounters. This self-crossing leads to the formation of shocks that increase the thermal energy of the debris. If this thermal energy is efficiently radiated away, the debris settle in a narrow ring at the circularization radius with shock-induced luminosities of $\sim 10-10^3 \, L_{\rm Edd}$. If instead cooling is impeded, the debris form an extended torus located between the circularization radius and the semi-major axis of the star with heating rates $\sim 1-10^2 \, L_{\rm Edd}$. Extrapolating our results to parabolic orbits, we infer that circularization would occur via the same mechanism in $\sim 1$ period of the most bound debris for deeply penetrating encounters to $\sim 10$ for grazing ones. We also anticipate the same effect of the cooling efficiency on the structure of the disc with associated luminosities of $\sim 1-10 \, L_{\rm Edd}$ and heating rates of $\sim 0.1-1 \, L_{\rm Edd}$. In the latter case of inefficient cooling, we deduce a viscous timescale generally shorter than the circularization timescale. This suggests an accretion rate through the disc tracing the fallback rate, if viscosity starts acting promptly.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dust trapping by spiral arms in gravitationally unstable protostellar discs

In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Estimating the fossil disc mass during supermassive black hole mergers: the importance of torque implementation

In this paper, we revisit the issue of estimating the "fossil" disc mass in the circumprimary disc, during the merger of a supermassive black hole binary. As the binary orbital decay speeds up due to the emission of gravitational waves, the gas in the circumprimary disc might be forced to accrete rapidly and could in principle provide a significant electromagnetic counterpart to the gravitational wave emission. Since the luminosity of such flare is proportional to the gaseous mass in the circumprimary disc, estimating such mass accurately is important. Previous investigations of this issue have produced contradictory results, with some authors estimating super-Eddington flares and large disc mass, while others suggesting that the "fossil" disc mass is very low, even less than a Jupiter mass. Here, we perform simple 1D calculations to show that such very low estimates of the disc mass are an artifact of the specific implementation of the tidal torque in 1D models. In particular, for moderate mass ratios of the binary, the usual formula for the torque used in 1D models significantly overestimates the width of the gap induced by the secondary and this artificially leads to a very small leftover circumprimary disc. Using a modified torque, calibrated to reproduce the correct gap width as estimated by 3D models, leads to fossil disc masses of the order of one solar mass. The rapid accretion of the whole circumprimary disc would produce peak luminosities of the order of 1-20 times the Eddington luminosity. Even if a significant fraction of the gas escapes accretion by flowing out the secondary orbit during the merger (an effect not included in our calculations), we would still predict close to Eddington luminosities that might be easily detected.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On planet formation in HL Tau

We explain the axisymmetric gaps seen in recent long-baseline observations of the HL Tau protoplanetary disc with the Atacama Large Millimetre/Submillimetre Array (ALMA) as being due to the different response of gas and dust to embedded planets in protoplanetary discs. We perform global, three dimensional dusty smoothed particle hydrodynamics calculations of multiple planets embedded in dust/gas discs which successfully reproduce most of the structures seen in the ALMA image. We find a best match to the observations using three embedded planets with masses of 0.2, 0.27 and 0.55 $M_{\rm J}$ in the three main gaps observed by ALMA, though there remain uncertainties in the exact planet masses from the disc model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spin alignment and differential accretion in merging black hole binaries

Interactions between a supermassive black hole binary and the surrounding accretion disc can both assist the binary inspiral and align the black hole spins to the disc angular momentum. While binary migration is due to angular-momentum transfer within the circumbinary disc, the spin-alignment process is driven by the mass accreting on to each black hole. Mass transfer between different disc components thus couples the inspiral and the alignment process together. Mass is expected to leak through the cavity cleared by the binary, and preferentially accretes on to the lighter (secondary) black hole which orbits closer to the disc edge. Low accretion rate on to the heavier (primary) black hole slows the alignment process down. We revisit the problem and develop a semi-analytical model to describe the coupling between gas-driven inspiral and spin alignment, finding that binaries with mass ratio q<~0.2 approach the gravitational-wave driven inspiral in differential misalignment: light secondaries prevent primaries from aligning. Binary black holes with misaligned primaries are ideal candidates for precession effects in the strong-gravity regime and may suffer from moderately large (~1500 km/s) recoil velocities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Episodic Accretion in Young Stars

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How to Detect the Signatures of Self-Gravitating Circumstellar Discs with the Atacama Large Millimetre/sub-millimetre Array

In this paper we present simulated Atacama Large Millimetre/sub-millimetre Array (ALMA) observations of self-gravitating circumstellar discs with different properties in size, mass and inclination, located in four of the most extensively studied and surveyed star-forming regions. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation and representative dust opacities, we have initially constructed maps of the expected emission at sub-mm wavelengths of a large sample of discs with different properties. We have then simulated realistic observations of discs as they may appear with ALMA using the Common Astronomy Software Application ALMA simulator. We find that, with a proper combination of antenna configuration and integration time, the spiral structure characteristic of self-gravitating discs is readily detectable by ALMA over a wide range of wavelengths at distances comparable to TW Hydrae ($\sim 50 \,$pc), Taurus - Auriga and Ophiucus ($\sim 140 \,$pc) star-forming regions. However, for discs located in Orion complex ($\sim 400 \,$pc) only the largest discs in our sample (outer radius of 100 au) show a spatially resolved structure while the smaller ones (outer radius of 25 au) are characterized by a spiral structure that is not conclusively detectable with ALMA.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probing the presence of planets in transition discs' cavities via warps: the case of TW Hya

We are entering the era in which observations of protoplanetary discs properties can indirectly probe the presence of massive planets or low mass stellar companions interacting with the disc. In particular, the detection of warped discs can provide important clues to the properties of the star-disc system. In this paper we show how observations of warped discs can be used to infer the dynamical properties of the systems. We concentrate on circumbinary discs, where the mass of the secondary can be planetary. First, we provide some simple relations that link the amplitude of the warp in the linear regime to the parameters of the system. Secondly, we apply our method to the case of TW Hya, a transition disc for which a warp has been proposed based on spectroscopic observations. Assuming values for the disc and stellar parameters from observations, we conclude that, in order for a warp induced by a planetary companion to be detectable, the planet mass should be large ($M_{\rm p} \approx 10 - 14M_{\rm J}$) and the disc should be viscous ($α\approx 0.15 - 0.25$). We also apply our model to LkCa 15 and T Cha, where a substellar companion has been detected within the central cavity of the transition discs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The puzzling source IGR J17361-4441 in NGC 6388: a possible planetary tidal disruption event

On 2011 August 11, INTEGRAL discovered the hard X-ray source IGR J17361-4441 near the centre of the globular cluster NGC 6388. Follow up observations with Chandra showed the position of the transient was inconsistent with the cluster dynamical centre, and thus not related to its possible intermediate mass black hole. The source showed a peculiar hard spectrum (Gamma \approx 0.8) and no evidence of QPOs, pulsations, type-I bursts, or radio emission. Based on its peak luminosity, IGR J17361-4441 was classified as a very faint X-ray transient, and most likely a low-mass X-ray binary. We re-analysed 200 days of Swift/XRT observations, covering the whole outburst of IGR J17361-4441 and find a t^{-5/3} trend evident in the light curve, and a thermal emission component that does not evolve significantly with time. We investigate whether this source could be a tidal disruption event, and for certain assumptions find an accretion efficiency epsilon \approx 3.5E-04 (M_{Ch}/M) consistent with a massive white dwarf, and a disrupted minor body mass M_{mb}=1.9E+27(M/M_{Ch}) g in the terrestrial-icy planet regime. These numbers yield an inner disc temperature of the order kT_{in} \approx 0.04 keV, consistent with the blackbody temperature of kT_{in} \approx 0.08 keV estimated by spectral fitting. Although the density of white dwarfs and the number of free-floating planets are uncertain, we estimate the rate of planetary tidal disruptions in NGC 6388 to be in the range 3E-06 to 3E-04 yr^{-1}. Averaged over the Milky Way globular clusters, the upper limit value corresponds to 0.05 yr^{-1}, consistent with the observation of a single event by INTEGRAL and Swift.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Protoplanetary Disk Masses from Stars to Brown Dwarfs

We present SCUBA-2 850um observations for 7 very low mass stars (VLMS) and brown dwarfs (BDs): 3 in Taurus, 4 in the TWA, and all classical T Tauri (cTT) analogs. We detect 2 of the 3 Taurus disks, but none of the TWA ones. Our 3sigma limits correspond to a dust mass of 1.2 MEarth in Taurus and a mere 0.2 MEarth in the TWA (3--10x deeper than previous work). We combine our data with other sub-mm/mm surveys of Taurus, rho Oph and the TWA to investigate trends in disk mass and grain growth during the cTT phase. We find : (1) The minimum disk outer radius required to explain the upper envelope of sub-mm/mm fluxes is 100 AU for intermediate-mass stars, solar-types and VLMS, and 20 AU for BDs. (2) While the upper envelope of disk masses increases with Mstar from BDs to VLMS to solar-types, no increase is seen from solar-type to intermediate-mass stars. We propose this is due to enhanced photoevaporation around intermediate masses. (3) Many disks around Taurus and rho Oph intermediate-mass and solar-type stars evince an opacity index beta of 0--1, indicating large grains. Of the only four VLMS/BDs in these regions with multi-wavelength data, three are consistent with large grains, though optically thick disks are not ruled out. (4) For the TWA VLMS (TWA 30A, B), combining our fluxes with Mdot and ages suggests substantial grain growth by 10 Myr. The degree of grain growth in the TWA BDs (2M1207A, SSPM1102) remains largely unknown. (5) A Bayesian analysis shows that mean(log[Mdisk/Mstar]) = -2.4, roughly constant all the way from intermediate-mass stars to VLMS/BDs, and (6) the disk mass in close solar-type Taurus binaries is significantly lower than in singles (by a factor of 10), while that in wide solar-type Taurus binaries is closer to that in singles (lower by a factor of 3). (7) We discuss the implications for planet formation, and for the dependence of Mdot on Mstar.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: Luminous extragalactic transients

We discuss the importance and potential contribution of Athena+ to the science questions related to gamma-ray bursts, tidal disruption events and supernova shock break-out. Athena+ will allow breakthrough observations involving high-resolution X-ray spectroscopic observations of high-z gamma-ray bursts, observations of tidal disruption events tailored to determine the mass and potentially the spin of the black hole responsible for the tidal disruption and observations of X-rays from the supernova shock breakout providing a measure of the radius of the exploding star or of the companion in the case of type Ia supernovae. We briefly discuss survey facilities that will provide triggers to these events envisaged to be operational around 2028.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wave-like warp propagation in circumbinary discs I. Analytic theory and numerical simulations

In this paper we analyse the propagation of warps in protostellar circumbinary discs. We use these systems as a test environment in which to study warp propagation in the bending-wave regime, with the addition of an external torque due to the binary gravitational potential. In particular, we want to test the linear regime, for which an analytic theory has been developed. In order to do so, we first compute analytically the steady state shape of an inviscid disc subject to the binary torques. The steady state tilt is a monotonically increasing function of radius. In the absence of viscosity, the disc does not present any twist. Then, we compare the time-dependent evolution of the warped disc calculated via the known linearised equations both with the analytic solutions and with full 3D numerical simulations, which have been performed with the PHANTOM SPH code using 2 million particles. We find a good agreement both in the tilt and in the phase evolution for small inclinations, even at very low viscosities. Moreover, we have verified that the linearised equations are able to reproduce the diffusive behaviour when α > H/R, where α is the disc viscosity parameter. Finally, we have used the 3D simulations to explore the non-linear regime. We observe a strongly non-linear behaviour, which leads to the breaking of the disc. Then, the inner disc starts precessing with its own precessional frequency. This behaviour has already been observed with numerical simulations in accretion discs around spinning black holes. The evolution of circumstellar accretion discs strongly depends on the warp evolution. Therefore the issue explored in this paper could be of fundamental importance in order to understand the evolution of accretion discs in crowded environments, when the gravitational interaction with other stars is highly likely, and in multiple systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wave-like warp propagation in circumbinary discs II. Application to KH 15D

KH 15D is a protostellar binary system that shows a peculiar light curve. In order to model it, a narrow circumbinary precessing disc has been invoked, but a proper dynamical model has never been developed. In this paper, we analytically address the issue of whether such a disc can rigidly precess around KH 15D, and we relate the precessional period to the main parameters of the system. Then, we simulate the disc's dynamics by using a 1D model developed in a companion paper, such that the warp propagates into the disc as a bending wave, which is expected to be the case for protostellar discs. The validity of such an approach has been confirmed by comparing its results with full 3D SPH simulations on extended discs. In the present case, we use this 1D code to model the propagation of the warp in a narrow disc. If the inner truncation radius of the disc is set by the binary tidal torques at {\sim} 1 AU, we find that the disc should extend out to 6-10 AU (depending on the models), and is therefore wider than previously suggested. Our simulations show that such a disc does reach an almost steady state, and then precesses as a rigid body. The disc displays a very small warp, with a tilt inclination that increases with radius in order to keep the disc in equilibrium against the binary torque. However, for such wider discs, the presence of viscosity leads to a secular decay of the tilt on a timescale of {\approx} 3000 (α/0.05)^(-1) years, where α is the disc viscosity parameter. The presence of a third body (such as a planet), orbiting at roughly 10 AU might simultaneously explain the outer truncation of the disc and the maintenance of the tilt for a prolonged time.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Black hole mergers: do gas discs lead to spin alignment?

In this Letter we revisit arguments suggesting that the Bardeen-Petterson effect can coalign the spins of a central supermassive black hole binary accreting from a circumbinary (or circumnuclear) gas disc. We improve on previous estimates by adding the dependence on system parameters, and noting that the nonlinear nature of warp propagation in a thin viscous disc affects alignment. This reduces the disc's ability to communicate the warp, and can severely reduce the effectiveness of disc-assisted spin alignment. We test our predictions with a Monte Carlo realization of random misalignments and accretion rates and we find that the outcome depends strongly on the spin magnitude. We estimate a generous upper limit to the probability of alignment by making assumptions which favour it throughout. Even with these assumptions, about 40% of black holes with $a \gtrsim 0.5$ do not have time to align with the disc. If the residual misalignment is not small and it is maintained down to the final coalescence phase this can give a powerful recoil velocity to the merged hole. Highly spinning black holes are thus more likely of being subject to strong recoils, the occurrence of which is currently debated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Challenges in the modelling of tidal disruption events lightcurves

In this contribution, I review the recent developments on the modelling of the lightcurve of tidal disruption events. Our understanding has evolved significantly from the earlier seminal results that imply a simple power-law decay of the bolometric light curve as $t^{-5/3}$. We now know that the details of the rise to the peak of the lightcurve is determined mainly by the internal structure of the disrupted star. We also have improved models for the disc thermal emission, showing that in this case the decline of the luminosity with time should be much flatter than the standard $t^{-5/3}$ law, especially in optical and UV wavelengths, while the X-ray lightcurve is generally best suited to track the bolometric one. Finally, we are just starting to explore the interesting general relativistic effects that might arise for such events, for which the tidal radius lies very close to the black hole event horizon. In particular, I describe here some possible evidences for relativistic Lense-Thirring precession from the light curve of the event Swift J1644.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Response of a circumbinary accretion disc to black hole mass loss

We investigate the evolution of the surface density of a circumbinary accretion disc after the mass loss induced by the merger of two supermassive black holes. We first introduce an analytical model, under the assumption of a disc composed of test particles, to derive the surface density evolution of the disc following the mass loss. The model predicts the formation of sharp density peaks in the disc; the model also allows us to compute the typical timescale for the formation of these peaks. To test and validate the model, we run numerical simulations of the process using the Smoothed Particle Hydrodynamics (SPH) code PHANTOM, taking fluid effects into account. We find good agreement in the shape and position of the peaks between the model and the simulations. In a fluid disc, however, the epicyclic oscillations induced by the mass loss can dissipate, and only some of the predicted peaks form in the simulation. To quantify how fast this dissipation proceeds, we introduce an appropriate parameter, and we show that it is effective in explaining the differences between the analytical, collisionless model and a real fluid disc.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fu Ori outbursts and the planet-disc mass exchange

It has been recently proposed that giant protoplanets migrating inward through the disc more rapidly than they contract could be tidally disrupted when they fill their Roche lobes ~ 0.1 AU away from their parent protostars. Here we consider the process of mass and angular momentum exchange between the tidally disrupted planet and the surrounding disc in detail. We find that the planet's adiabatic mass-radius relation and its ability to open a deep gap in the disc determine whether the disruption proceeds as a sudden runaway or a balanced quasi-static process. In the latter case the planet feeds the inner disc through its Lagrangian L1 point like a secondary star in a stellar binary system. As the planet loses mass it gains specific angular momentum and normally migrates in the outward direction until the gap closes. Numerical experiments show that planet disruption outbursts are preceded by long "quiescent" periods during which the disc inward of the planet is empty. The hole in the disc is created when the planet opens a deep gap, letting the inner disc to drain onto the star while keeping the outer one stalled behind the planet. We find that the mass-losing planet embedded in a realistic protoplanetary disc spawns an extremely rich set of variability patterns. In a subset of parameter space, there is a limit cycle behaviour caused by non-linear interaction between the planet mass loss and the disc hydrogen ionisation instability. We suggest that tidal disruptions of young massive planets near their stars may be responsible for the observed variability of young accreting protostars such as FU Ori, EXor and T Tauri stars in general.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multiband lightcurves of tidal disruption events

Unambiguous detection of the tidal disruption of a star would allow an assessment of the presence and masses of supermassive black holes in quiescent galaxies. It would also provide invaluable information on bulge scale stellar processes (such as two-body relaxation) via the rate at which stars are injected into the tidal sphere of influence of the black holes. This rate, in turn, is essential to predict gravitational radiation emission by compact object inspirals. The signature of a tidal disruption event is thought to be a fallback rate for the stellar debris onto the black hole that decreases as $t^{-5/3}$. This mass flux is often assumed to yield a luminous signal that decreases in time at the same rate. In this paper, we calculate the monochromatic lightcurves arising from such an accretion event. Differently from previous studies, we adopt a more realistic description of the fallback rate and of the super-Eddigton accretion physics. We also provide simultaneous lightcurves in optical, UV and X-rays. We show that, after a few months, optical and UV lightcurves scale as $t^{-5/12}$, and are thus substantially flatter than the $t^{-5/3}$ behaviour, which is a prerogative of the bolometric lightcurve, only. At earlier times and for black hole masses $< 10^7 M_{\sun}$, the wind emission dominates: after reaching a peak of $10^{41}-10^{43}$ erg/s at roughly a month, the lightcurve decreases steeply as $\sim t^{-2.6}$, until the disc contribution takes over. The X-ray band, instead, is the best place to detect the $t^{-5/3}$ "smoking gun" behaviour, although it is displayed only for roughly a year, before the emission steepens exponentially.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resolution requirements for Smoothed Particle Hydrodynamics simulations of self-gravitating accretion discs

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 $15Ω^{-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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Nature of Angular Momentum Transport in Radiative Self-Gravitating Protostellar Discs

Semi-analytic models of self-gravitating discs often approximate the angular momentum transport generated by the gravitational instability using the phenomenology of viscosity. This allows the employment of the standard viscous evolution equations, and gives promising results. It is, however, still not clear when such an approximation is appropriate. This paper tests this approximation using high resolution 3D smoothed particle hydrodynamics (SPH) simulations of self-gravitating protostellar discs with radiative transfer. The nature of angular momentum transport associated with the gravitational instability is characterised as a function of both the stellar mass and the disc-to-star mass ratio. The effective viscosity is calculated from the Reynolds and gravitational stresses in the disc. This is then compared to what would be expected if the effective viscosity were determined by assuming local thermodynamic equilibrium or, equivalently, that the local dissipation rate matches the local cooling rate. In general, all the discs considered here settle into a self-regulated state where the heating generated by the gravitational instability is modulated by the local radiative cooling. It is found that low-mass discs can indeed be represented by a local "alpha-parametrisation", provided that the disc aspect ratio is small (H/R < 0.1) which is generally the case when the disc-to-star mass ratio q <0.5. However, this result does not extend to discs with masses approaching that of the central object. These are subject to transient burst events and global wave transport, and the effective viscosity is not well modelled by assuming local thermodynamic equilibrium. In spite of these effects, it is shown that massive (compact) discs can remain stable and not fragment, evolving rapidly to reduce their disc-to-star mass ratios through stellar accretion and radial spreading.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The effects of opacity on gravitational stability in protoplanetary discs

In this paper we consider the effects of opacity regimes on the stability of self-gravitating protoplanetary discs to fragmentation into bound objects. Using a self-consistent 1-D viscous disc model, we show that the ratio of local cooling to dynamical timescales Omega*tcool has a strong dependence on the local temperature. We investigate the effects of temperature-dependent cooling functions on the disc's gravitational stability through controlled numerical experiments using an SPH code. We find that such cooling functions raise the susceptibility of discs to fragmentation through the influence of temperature perturbations - the average value of Omega*tcool has to increase to prevent local variability leading to collapse. We find the effects of temperature dependence to be most significant in the "opacity gap" associated with dust sublimation, where the average value of Omega*tcool at fragmentation is increased by over an order of magnitude. We then use this result to predict where protoplanetary discs will fragment into bound objects, in terms of radius and accretion rate. We find that without temperature dependence, for radii < ~10AU a very large accretion rate ~10^-3 Msun/yr is required for fragmentation, but that this is reduced to 10^-4 Msun/yr with temperature-dependent cooling. We also find that the stability of discs with accretion rates < ~10^-7 Msun/yr at radii > ~50AU is enhanced by a lower background temperature if the disc becomes optically thin.

0 citations0 reviewsNo code linkedOpen access