We study the evolution of planet-induced vortices in radially stratified disks, with initial conditions allowing for radial buoyancy. For this purpose we run global two dimensional hydrodynamical simulations, using the PLUTO code. Planet-induced vortices are a product of the Rossby wave instability (RWI) triggered in the edges of a planetary gap. In this work we assess the influence of radial buoyancy for the development of the vortices. We found that radial buoyancy leads to smoother planetary gaps, which generates weaker vortices. This effect is less pronounced for locally isothermal and quasi-isothermal (very small cooling rate) disks. We observed the formation of two generations of vortices. The first generation of vortices is formed in the outer wall of the planetary gap. The merged primary vortex induces accretion, depleting the mass on its orbit. This process creates a surface density enhancement beyond the primary vortex position. The second generation of vortices arise in this surface density enhancement, indicating that the bump in this region is sufficient to trigger the RWI. The merged secondary vortex is a promising explanation for the location of the vortex in the Oph IRS48 system. Finally, we observed a nonmonotonic behavior for the vortex lifetimes as a function of the thermal relaxation timescale, agreeing with previous studies. The birth times of the secondary vortices also display a nonmonotonic behavior, which is correlated with the growth time of the primary vortex and its induced accretion.

Circumstellar disks are considered to be the birthplace of planets. Specific structures like spiral arms, gaps, and cavities are characteristic indicators of planet-disk interaction. Investigating these structures can provide insights into the growth of protoplanets and the physical properties of the disk. We investigate the feasibility of using molecular lines to trace planet-induced structures in circumstellar disks. Based on 3D hydrodynamic simulations of planet-disk interactions, we perform self-consistent temperature calculations and produce N-LTE molecular line velocity-channel maps and spectra of these disks using our new N-LTE line radiative transfer code Mol3D. Subsequently, we simulate ALMA observations using the CASA simulator. We consider two nearly face-on inclinations, 5 disk masses, 7 disk radii, and 2 different typical pre-main-sequence host stars (T Tauri, Herbig Ae). We calculate up to 141 individual velocity-channel maps for five molecules/isotopoloques in a total of 32 rotational transitions to investigate the frequency dependence of the structures indicated above. We find that the majority of protoplanetary disks in our parameter space could be detected in the molecular lines considered. However, unlike the continuum case, gap detection is not straightforward in lines. For example, gaps are not seen in symmetric rings but are masked by the pattern caused by the global (Keplerian) velocity field. We identify specific regions in the velocity-channel maps that are characteristic of planet-induced structures. Simulations of high angular resolution molecular line observations demonstrate the potential of ALMA to provide complementary information about the planet-disk interaction as compared to continuum observations. In particular, the detection of planet-induced gaps is possible under certain conditions.(abridged)

Close-in exoplanets interact with their host stars gravitationally as well as via their magnetized plasma outflows. The rich dynamics that arises may result in distinct observable features. Our objective is to study and classify the morphology of the different types of interaction that can take place between a giant close-in planet (a Hot Jupiter) and its host star, based on the physical parameters that characterize the system. We perform 3D magnetohydrodynamic numerical simulations to model the star--planet interaction, incorporating a star, a Hot Jupiter, and realistic stellar and planetary outflows. We explore a wide range of parameters and analyze the flow structures and magnetic topologies that develop. Our study suggests the classification of star--planet interactions into four general types, based on the relative magnitudes of three characteristic length scales that quantify the effects of the planetary magnetic field, the planetary outflow, and the stellar gravitational field in the interaction region. We describe the dynamics of these interactions and the flow structures that they give rise to, which include bow shocks, cometary-type tails, and inspiraling accretion streams. We point out the distinguishing features of each of the classified cases and discuss some of their observationally relevant properties. The magnetized interactions of star--planet systems can be categorized, and their general morphologies predicted, based on a set of basic stellar, planetary, and orbital parameters.

We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a combination of numerical simulations in 2D and 3D and a linear perturbation analysis. Steady-state models of such disks require the inclusion of magnetic diffusivity. To make progress using ideal MHD, we focus on simplified field configurations, involving purely vertical ($B_z$) and azimuthal ($B_\varphi$) field components and a combination of the two. For each of the models we calculate the locations of the relevant resonances and of the turning points, which delineate the propagation regions of the MHD waves that transport angular momentum from the planet to the disk. We use both numerical and semianalytic methods to evaluate the cumulative back torque acting on the planet, and explore the effect of spatial gradients in the disk's physical variables on the results. We conclude that, under realistic (3D) circumstances, a large-scale magnetic field can slow down the inward migration that characterizes the underlying unmagnetized disk --- by up to a factor of $\sim 2$ when the magnetic pressure approaches the thermal pressure --- but it cannot reverse it. A previous inference that a pure-$B_\phi$ field whose amplitude decreases fast enough with radius leads to outward migration applies only in 2D. In fact, we find that, in 3D, a pure-$B_\phi$ disk undergoes a rapid transition to turbulence on account of a magnetorotational instability that is triggered by the planet-induced appearance of a weak $B_z$ component.

We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a linear perturbation analysis in the ideal-MHD limit. We focus on wind-driving disks, in which a magnetic torque $\propto B_{0z} \partial B_{0\varphi}/\partial z$ (where $B_{0z}$ and $B_{0\varphi}$ are the equilibrium vertical and azimuthal field components) induces vertical angular momentum transport. We derive the governing differential equation for the disk response and identify its resonances and turning points. For a disk containing a slightly subthermal, pure-$B_{0z}$ field, the total 3D torque is close to its value in the 2D limit but remains lower than the hydrodynamic torque. In contrast with the 2D pure-$B_{0\varphi}$ field model considered by Terquem (2003), inward migration is not reduced in this case when the field amplitude decreases with radius. The presence of a subdominant $B_{0\varphi}$ component whose amplitude increases from zero at $z=0$ has little effect on the torque when acting alone, but in conjunction with a $B_{0z}$ component it gives rise to a strong torque that speeds up the inward migration by a factor $\gtrsim 200$. This factor could, however, be reduced in a real disk by dissipation and magnetic diffusivity effects. Unlike all previously studied disk migration models, in the $B_{0z}+\partial B_{0\varphi}/\partial z$ case the dominant contributions to the torque add with the same sign from the two sides of the planet. We attribute this behavior to a new mode of interaction wherein a planet moves inward by plugging into the disk's underlying angular momentum transport mechanism.

Young giant planets, which are embedded in a circumstellar disk, will significantly perturb the disk density distribution. This effect can potentially be used as an indirect tracer for planets. We investigate the feasibility of observing planet-induced gaps in circumstellar disks in scattered light. We perform 3D hydrodynamical disk simulations combined with subsequent radiative transfer calculations in scattered light for different star, disk, and planet configurations. The results are compared to those of a corresponding study for the (sub)mm thermal re-emission. The feasibility of detecting planet-induced gaps in scattered light is mainly influenced by the optical depth of the disk and therefore by the disk size and mass. Planet-induced gaps are in general only detectable if the photosphere of the disks is sufficiently disturbed. Within the limitations given by the parameter space here considered, we find that gap detection is possible in the case of disks with masses below $\sim 10^{-4\dots-3} \, \rm M_\odot$. Compared to the disk mass that marks the lower Atacama Large (Sub)Millimeter Array (ALMA) detection limit for the thermal radiation re-emitted by the disk, it is possible to detect the same gap both in re-emission and scattered light only in a narrow range of disk masses around $\sim 10^{-4} \, \rm M_\odot$, corresponding to $16\%$ of cases considered in our study.

We have performed three-dimensional magneto-hydrodynamical simulations of stellar accretion disks, using the PLUTO code, and studied the accretion of gas onto a Jupiter-mass planet and the structure of the circumplanetary gas flow after opening a gap in the disk. We compare our results with simulations of laminar, yet viscous disks with different levels of an $\alpha$-type viscosity. In all cases, we find that the accretion flow across the surface of the Hill sphere of the planet is not spherically or azimuthally symmetric, and is predominantly restricted to the mid-plane region of the disk. Even in the turbulent case, we find no significant vertical flow of mass into the Hill sphere. The outer parts of the circumplanetary disk are shown to rotate significantly below Keplerian speed, independent of viscosity, while the circumplanetary disk density (therefore the angular momentum) increases with viscosity. For a simulation of a magnetized turbulent disk, where the global averaged alpha stress is $\alpha_{MHD}=10^{-3}$, we find the accretion rate onto the planet to be $\dot{M}\approx 2\times10^{-6}M_{J}yr^{-1}$ for a gap surface density of $12 g cm^{-2}$. This is about a third of the accretion rate obtained in a laminar viscous simulation with equivalent $\alpha$ parameter.

Planets are supposed to form in circumstellar disks. The gravitational potential of a planet perturbs the disk and leads to characteristic structures, i.e. spiral waves and gaps, in the disk's density profile. We perform a large-scale parameter study of the observability of these planet-induced structures in circumstellar disks with ALMA. On the basis of HD and MHD simulations, we calculated the disk temperature structure and (sub)mm images of these systems. These were used to derive simulated ALMA images. Because appropriate objects are frequent in Taurus, we focused on a distance of 140pc and a declination of 20{\deg}. The explored range of star-disk-planet configurations consists of 6 HD simulations (including magnetic fields and different planet masses), 9 disk sizes, 15 total disk masses, 6 different central stars, and two different grain size distributions. On almost all scales and in particular down to a scale of a few AU, ALMA is able to trace disk structures induced by planet-disk interaction or by the influence of magnetic fields on the wavelength range between 0.4 and 2.0mm. In most cases, the optimum angular resolution is limited by the sensitivity. However, within the range of typical masses of protoplanetary disks (0.1-0.001Msun) the disk mass has a minor impact on the observability. It is possible to resolve disks down to 2.67e-6Msun and trace gaps induced by a planet with M_p/M_s = 0.001 in disks with 2.67e-4Msun with a signal-to-noise ratio greater than three. The central star has a major impact on the observability of gaps, as well as the considered maximum grainsize of the dust in the disk. In general, it is more likely to trace planet-induced gaps in our MHD models, because gaps are wider in the presence of magnetic fields. We also find that zonal flows resulting from MRI create gap-like structures in the disk's re-emission radiation, which are observable with ALMA.

In order to explain grain growth to mm sized particles and their retention in outer regions of protoplanetary disks, as it is observed at sub-mm and mm wavelengths, we investigate if strong inhomogeneities in the gas density profiles can slow down excessive radial drift and can help dust particles to grow. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million years time scales. In addition, we introduced a comparison between the pressure inhomogeneities considered in this work and the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study if there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. We present the favorable conditions to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of the Taurus Auriga or Ophiucus star forming regions.

We performed 3D MHD simulations of planet migration in stratified disks using the Godunov code PLUTO, where the disk is turbulent due to the magnetorotational instability. We study the migration for planets with different planet-star mass ratios $q=M_{p}/M_{s}$. In agreement with previous studies, for the low-mass planet cases ($q=5\times10^{-6}$ and $10^{-5}$), migration is dominated by random fluctuations in the torque. For a Jupiter-mass planet $(q=M_{p}/M_{s}=10^{-3}$ for $M_{s}=1M_{\odot})$, we find a reduction of the magnetic stress inside the orbit of the planet and around the gap region. After an initial stage where the torque on the planet is positive, it reverses and we recover migration rates similar to those found in disks where the turbulent viscosity is modelled by an $\alpha$ viscosity. For the intermediate-mass planets ($q=5\times10^{-5}, 10^{-4}$ and $2\times10^{-4}$) we find a new and so far unexpected behavior. In some cases they experience sustained and systematic outwards migration for the entire duration of the simulation. For this case, the horseshoe region is resolved and torques coming from the corotation region can remain unsaturated due to the stresses in the disk. These stresses are generated directly by the magnetic field. The magnitude of the horseshoe drag can overcome the negative Lindblad contribution when the local surface density profile is flat or increasing outwards, which we see in certain locations in our simulations due to the presence of a zonal flow. The intermediate-mass planet is migrating radially outwards in locations where there is a positive gradient of a pressure bump (zonal flow).