No Arabic abstract
We study the disc planet interactions of low-mass protoplanets embedded in a circumstellar disc. We extend the standard theory of planet migration from the usual locally isothermal assumption to include non-barotropic effects, focusing on the validity of linear theory. We compared solutions of the linear equations with results from non-linear hydrodynamic simulations, where in both cases we adopted a background entropy gradient and solved the energy equation. We show that the migration behavior of embedded planets depends critically on the background radial entropy gradient in the disc. The presence of such a gradient not only changes the corotation torque on the planet, but also always guarantees a departure from linear behavior, which gives a singular density response at corotation, in the absence of thermal or viscous diffusion. A negative entropy gradient tends to give rise to positive, non-linear corotation torques apparently produced as material executes horseshoe turns at approximately constant entropy. These torques have no counterpart in linear theory, but can be strong enough to push the planet outwards until saturation starts to occur after a horseshoe libration period. Increased thermal diffusion acts to reduce these non-linear torques, but, at the same time, it can help to prevent their saturation. In combination with a small kinematic viscosity that is able to maintain a smooth density profile the positive torque could be sustained.
Discs of gas and dust surrounding young stars are the birthplace of planets. However, direct detection of protoplanets forming within discs has proved elusive to date. We present the detection of a large, localized deviation from Keplerian velocity in the protoplanetary disc surrounding the young star HD163296. The observed velocity pattern is consistent with the dynamical effect of a two Jupiter-mass planet orbiting at a radius $approx$ 260au from the star.
We present a new method for embedding a stellar disc in a cosmological dark matter halo and provide a worked example from a {Lambda}CDM zoom-in simulation. The disc is inserted into the halo at a redshift z = 3 as a zero-mass rigid body. Its mass and size are then increased adiabatically while its position, velocity, and orientation are determined from rigid-body dynamics. At z = 1, the rigid disc is replaced by an N-body disc whose particles sample a three-integral distribution function (DF). The simulation then proceeds to z = 0 with live disc and halo particles. By comparison, other methods assume one or more of the following: the centre of the rigid disc during the growth phase is pinned to the minimum of the halo potential, the orientation of the rigid disc is fixed, or the live N-body disc is constructed from a two rather than three-integral DF. In general, the presence of a disc makes the halo rounder, more centrally concentrated, and smoother, especially in the innermost regions. We find that methods in which the disc is pinned to the minimum of the halo potential tend to overestimate the amount of adiabatic contraction. Additionally, the effect of the disc on the subhalo distribution appears to be rather insensitive to the disc insertion method. The live disc in our simulation develops a bar that is consistent with the bars seen in late-type spiral galaxies. In addition, particles from the disc are launched or kicked up to high galactic latitudes.
We study mass outflows driven from accretion discs by radiation pressure due to spectral lines. To investigate non-axisymmetric effects, we use the Athena++ code and develop a new module to account for radiation pressure driving. In 2D, our new simulations are consistent with previous 2D axisymmetric solutions by Proga et al. who used the Zeus 2D code. Specifically, we find that the disc winds are time dependent, characterized by a dense stream confined to $sim 45^{circ}$ relative to the disc midplane and bounded on the polar side by a less dense, fast stream. Introducing a vertical, $phi$-dependent, subsonic velocity perturbation in the disc midplane does not change the overall character of the solution but global outflow properties such as the mass, momentum and kinetic energy fluxes are altered by up to 100%. Non-axisymmetric density structures develop and persist mainly at the base of the wind. They are relatively small, and their densities can be a few times higher that the azimuthal average. The structure of the non-axisymmetric and axisymmetric solutions differ also in other ways. Perhaps most importantly from the observational point of view are the differences in the so called clumping factors, that serve as a proxy for emissivity due to two body processes. In particular, the spatially averaged clumping factor over the entire fast stream, while it is of a comparable value in both solutions, it varies about 10 times faster in the non-axisymmetric case.
Transition discs are prime targets to look for protoplanets and study planet-disc interactions. We present VLT/SINFONI observations of PDS~70, a transition disc with a recently claimed embedded protoplanet. We take advantage of the angular and spectral diversity present in our data for an optimal PSF modeling and subtraction using principal component analysis (PCA). We report the redetection of PDS 70 b, both the front and far side of the outer disc edge, and the detection of several extended features in the annular gap. We compare spectral differential imaging applied before (PCA-SADI), and after (PCA-ASDI) angular differential imaging. Our tests suggest that PCA-SADI better recovers extended features, while PCA-ASDI is more sensitive to point sources. We adapted the negative fake companion (NEGFC) technique to infer the astrometry of the companion, and derived $r = 193.5 pm 4.9 mathrm{mas}$ and PA = 158.7deg $pm$ 3.0deg. We used both NEGFC and ANDROMEDA to infer the $K$-band spectro-photometry of the protoplanet, and found results consistent with recent VLT/SPHERE observations, except for their 2018/02 epoch measurement in the $K2$ filter. Finally, we derived an upper limit of $dot{M_b} < 1.26 times 10^{-7} big[ frac{5 M_{rm Jup}}{M_b} big] big[ frac{R_b}{R_{rm Jup}}big] M_{rm Jup} $ yr$^{-1}$ for the accretion rate of the companion based on an adaptation of PCA-SADI/PCA-ASDI around the Br$gamma$ line.
We utilize zoom-in cosmological simulations to study the nature of violent disc instability (VDI) in clumpy galaxies at high redshift, $z=1$--$5$. Our simulated galaxies are not in the ideal state assumed in Toomre instability, of linear fluctuations in an isolated, uniform, rotating disk. There, instability is characterised by a $Q$ parameter below unity, and lower when the disk is thick. Instead, the high-redshift discs are highly perturbed. Over long periods they consist of non-linear perturbations, compact massive clumps and extended structures, with new clumps forming in inter-clump regions. This is while the galaxy is subject to frequent external perturbances. We compute the local, two-component $Q$ parameter for gas and stars, smoothed on a $sim1~{rm kpc}$ scale to capture clumps of $10^{8-9}~{rm M}_odot$. The $Q<1$ regions are confined to collapsed clumps due to the high surface density there, while the inter-clump regions show $Q$ significantly higher than unity. Tracing the clumps back to their relatively smooth Lagrangian patches, we find that $Q$ prior to clump formation typically ranges from unity to a few. This is unlike the expectations from standard Toomre instability. We discuss possible mechanisms for high-$Q$ clump formation, e.g. rapid turbulence decay leading to small clumps that grow by mergers, non-axisymmetric instability, or clump formation induced by non-linear perturbations in the disk. Alternatively, the high-$Q$ non-linear VDI may be stimulated by the external perturbations such as mergers and counter-rotating streams. The high $Q$ may represent excessive compressive modes of turbulence, possibly induced by tidal interactions.