Researcher profile

Giuseppe Lodato

Giuseppe Lodato contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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.

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.

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.

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.

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]$.

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.

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.

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.

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.

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.

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.

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 &#34;kink&#34; 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.

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.

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.

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}$&#39;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.

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.

preprint2019arXiv

Effects of photoevaporation on protoplanetary disc `isochrones&#39;

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.

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.

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.