The mass, rate, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Advanced Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previous work showed that stellar mass black holes (sBHs) in an AGN disk migrate to regions of the disk, known as migration traps, where positive and negative gas torques cancel out, leading to frequent BBH formation. Here we build on that work by simulating the evolution of additional sBHs that enter the inner disk by either migration or inclination reduction. We also examine whether the BBHs formed in our models have retrograde or prograde orbits around their centers of mass with respect to the disk, determining the orientation, relative to the disk, of the spin of the merged BBHs. We find that orbiters entering the inner disk form BBHs, that are commonly asymmetric in mass, with sBHs on resonant orbits near the migration trap. When these black holes reach 80 solar masses, they form BBHs with sBHs in the migration trap, which over 10 Myr reach about 1000 solar masses. We also find that 68% of the BBHs in our simulation orbit in the retrograde direction, which implies our merger channel could easily produce the low values of the dimensionless aligned spin, chi_eff, observed in a majority of BBH mergers detected thus far.

The merger rate of stellar-mass black hole binaries (sBHBs) inferred by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) suggests the need for an efficient source of sBHB formation. Active galactic nucleus (AGN) disks are a promising location for the formation of these sBHBs, as well as binaries of other compact objects, because of powerful torques exerted by the gas disk. These gas torques cause orbiting compact objects to migrate towards regions in the disk where inward and outward torques cancel, known as migration traps. We simulate the migration of stellar mass black holes in an example of a model AGN disk, using an augmented N-body code that includes analytic approximations to migration torques, stochastic gravitational forces exerted by turbulent density fluctuations in the disk, and inclination and eccentricity dampening produced by passages through the gas disk, in addition to the standard gravitational forces between objects. We find that sBHBs form rapidly in our model disk as stellar-mass black holes migrate towards the migration trap. These sBHBs are likely to subsequently merge on short time-scales. The process continues, leading to the build-up of a population of over-massive stellar-mass black holes. The formation of sBHBs in AGN disks could contribute significantly to the sBHB merger rate inferred by LIGO.

Protoplanetary disks dissipate rapidly after the central star forms, on time-scales comparable to those inferred for planet formation. In order to allow the formation of planets, disks must survive the dispersive effects of UV and X-ray photoevaporation for at least a few Myr. Viscous accretion depletes significant amounts of the mass in gas and solids, while photoevaporative flows driven by internal and external irradiation remove most of the gas. A reasonably large fraction of the mass in solids and some gas get incorporated into planets. Here, we review our current understanding of disk evolution and dispersal, and discuss how these might affect planet formation. We also discuss existing observational constraints on dispersal mechanisms and future directions.

In this paper the migration of a 10 Earth-mass planetary core is investigated at the outer boundary of the dead zone of a protoplanetary disc by means of 2D hydrodynamic simulations done with the graphics processor unit version of the FARGO code. In the dead zone, the effective viscosity is greatly reduced due to the disc self-shielding against stellar UV radiation, X-rays from the stellar magnetosphere and interstellar cosmic rays. As a consequence, mass accumulation occurs near the outer dead zone edge, which is assumed to trap planetary cores enhancing the efficiency of the core-accretion scenario to form giant planets. Contrary to the perfect trapping of planetary cores in 1D models, our 2D numerical simulations show that the trapping effect is greatly dependent on the width of the region where viscosity reduction is taking place. Planet trapping happens exclusively if the viscosity reduction is sharp enough to allow the development of large-scale vortices due to the Rossby wave instability. The trapping is only temporarily, and its duration is inversely proportional to the width of the viscosity transition. However, if the Rossby wave instability is not excited, a ring-like axisymmetric density jump forms, which cannot trap the 10 Earth-mass planetary cores. We revealed that the stellar torque exerted on the planet plays an important role in the migration history as the barycentre of the system significantly shifts away from the star due to highly non-axisymmetric density distribution of the disc. Our results still support the idea of planet formation at density/pressure maximum, since the migration of cores is considerably slowed down enabling them further growth and runaway gas accretion in the vicinity of an overdense region.

Submillimetre images of transition discs are expected to reflect the distribution of the optically thin dust. Former observation of three transition discs LkHa330, SR21N, and HD1353444B at submillimetre wavelengths revealed images which cannot be modelled by a simple axisymmetric disc. We show that a large-scale anticyclonic vortex that develops where the viscosity has a large gradient (e.g., at the edge of the disc dead zone), might be accountable for these large-scale asymmetries. We modelled the long-term evolution of vortices being triggered by the Rossby wave instability. We found that a horseshoe-shaped (azimuthal wavenumber m=1) large-scale vortex forms by coalescing of smaller vortices within 5x10^4 yr, and can survive on the disc life-time (~5x10^6 yr), depending on the magnitude of global viscosity and the thickness of the viscosity gradient. The two-dimensional grid-based global disc simulations with local isothermal approximation and compressible-gas model have been done by the GPU version of hydrodynamic code FARGO (GFARGO). To calculate the dust continuum image at submillimetre wavelengths, we combined our hydrodynamical results with a 3D radiative transfer code. By the striking similarities of the calculated and observed submillimetre images, we suggest that the three transition discs can be modelled by a disc possessing a large-scale vortex formed near the disc dead zone edge. Since the larger dust grains (larger than mm in size) are collected in these vortices, the non-axisymmetric submillimetre images of the above transition discs might be interpreted as active planet and planetesimal forming regions situated far (> 50 AU) from the central stars.

Star formation occurs via fragmentation of molecular clouds, which means that the majority of stars born are a members of binaries. There is growing evidence that planets might form in circumprimary disks of medium-separation binaries. The tidal forces caused by the secondary generally act to distort the originally circular disk to an eccentric one. To infer the disk eccentricity from high-res NIR spectroscopy, we calculate the fundamental band emission lines of the CO molecule emerging from the atmosphere of the disk. We model circumprimary disk evolution under the gravitational perturbation of the orbiting secondary using a 2D grid-based hydrodynamical code, assuming alpha-type viscosity. The hydrodynamical results are combined with our spectral code based on the double-layer disk model to calculate the CO molecular line profiles. We find that the orbital velocity distribution of the gas parcels differs significantly from the circular Keplerian fashion, thus the line profiles are asymmetric in shape. The magnitude of asymmetry is insensitive to the binary mass ratio, the magnitude of viscosity, and the disk mass. In contrast, the disk eccentricity, thus the level of the line profile asymmetry, is influenced significantly by the binary eccentricity and the disk geometrical thickness. We demonstrate that the disk eccentricity profile in the planet-forming region can be determined by fitting the high-resolution CO line profile asymmetry using a simple 2D spectral model that accounts for the velocity distortions caused by the disk eccentricity. Thus, with our novel approach the disk eccentricity can be inferred with high-resolution near-IR spectroscopy prior to the era of high angular resolution optical or radio direct-imaging. By determining the disk eccentricity in medium-separation young binaries, we might be able to constrain the planet formation theories.

One of the longstanding unsolved problems of planet formation is how solid bodies of a few decimeters in size can "stick" to form large planetesimals. This is known as the "meter size barrier". In recent years it has become increasingly clear that some form of "particle trapping" must have played a role in overcoming the meter size barrier. Particles can be trapped in long-lived local pressure maxima, such as those in anticyclonic vortices, zonal flows or those believed to occur near ice lines or at dead zone boundaries. Such pressure traps are the ideal sites for the formation of planetesimals and small planetary embryos. Moreover, they likely produce large quantities of such bodies in a small region, making it likely that subsequent N-body evolution may lead to even larger planetary embryos. The goal of this Letter is to show that this indeed happens, and to study how efficient it is. In particular, we wish to find out if rocky/icy bodies as large as 10 M_earth can form within 1 Myr, since such bodies are the precursors of gas giant planets in the core accretion scenario.

Theories of planet formation predict the birth of giant planets in the inner, dense, and gas-rich regions of the circumstellar disks around young stars. These are the regions from which strong CO emission is expected. Observations have so far been unable to confirm the presence of planets caught in formation. We have developed a novel method to detect a giant planet still embedded in a circumstellar disk by the distortions of the CO molecular line profiles emerging from the protoplanetary disk's surface. The method is based on the fact that a giant planet significantly perturbs the gas velocity flow in addition to distorting the disk surface density. We have calculated the emerging molecular line profiles by combining hydrodynamical models with semianalytic radiative transfer calculations. Our results have shown that a giant Jupiter-like planet can be detected using contemporary or future high-resolution near-IR spectrographs such as VLT/CRIRES or ELT/METIS. We have also studied the effects of binarity on disk perturbations. The most interesting results have been found for eccentric circumprimary disks in mid-separation binaries, for which the disk eccentricity - detectable from the asymmetric line profiles - arises from the gravitational effects of the companion star. Our detailed simulations shed new light on how to constrain the disk kinematical state as well as its eccentricity profile. Recent findings by independent groups have shown that core-accretion is severely affected by disk eccentricity, hence detection of an eccentric protoplanetary disk in a young binary system would further constrain planet formation theories.

We examine the first phase of the core accretion model, namely the dust growth/fragmentation in binary systems. In our model, a gas and dust disk is present around the primary star and is perturbed by the secondary. We study the effects of a secondary with/without eccentricity on the dust population to determine what sizes the aggregates can reach and how that compares to the dust population in disks around single stars. We find that the secondary star has two effects on the dust population. 1.) The disk is truncated due to the presence of the secondary star and the maximum mass of the particles is decreased in the lowered gas densities. This effect is dominant in the outer disk. 2.) The perturbation of the secondary pumps up the eccentricity of the gas disk, which in turn increases the relative velocity between the dust and the gas. Therefore the maximum particle sizes are further decreased. The second effect of the secondary influences the entire disk. Coagulation is efficiently reduced even at the very inner parts of the disk. The average mass of the particles is reduced by four orders of magnitude (as a consequence, the stopping time is reduced by one order of magnitude) in disks around binary systems compared to dust in disks around single stars.

In this paper we intend to provide an indirect method to detect Jovian planets by studying near infrared emission spectra originating in the protoplanetary disks around T Tauri stars. Our idea is to investigate whether a massive planet could induce any observable effect on the spectral lines emerging in the disks atmosphere. As a tracer molecule we propose CO, which is excited in the ro-vibrational fundamental band in the disk atmosphere to a distance of ~2-3 AU (depending on the stellar mass) where terrestrial planets are thought to form. The synthetic molecular spectral line profiles were calculated by an own developed semi-analytical double layer disk model. 2D gas dynamics were incorporated in the calculation of synthetic spectral lines. We demonstrate that a massive planet embedded in a protoplanetary disk strongly influences the originally circular Keplerian gas dynamics. The perturbed motion of the gas can be detected by comparing the CO line profiles in emission, which is emerge from planet-bearing to those of planet-free disk models. The planet signal has two major characteristics: a permanent line profile asymmetry, and short timescale variability correlated with the orbital phase of the giant planet. We have found that the strength of the asymmetry depends on the physical parameters of the star-planet-disk system, such as the disk inclination angle, the planetary and stellar masses, the orbital distance, and the size of the disk inner cavity. The permanent line profile asymmetry is caused by a disk in an eccentric state in the gap opened by the giant planet. However, the variable component is a consequence of the local dynamical perturbation by the orbiting giant planet. We show that a forming giant planet, still embedded in the protoplanetary disk, can be detected using contemporary or future high-resolution near-IR spectrographs like VLT/CRIRES and ELT/METIS.

Among multi-planet planetary systems there are a large fraction of resonant systems. Studying the dynamics and formation of these systems can provide valuable informations on processes taking place in protoplanetary disks where the planets are thought have been formed. The recently discovered resonant system HD 60532 is the only confirmed case, in which the central star hosts a pair of giant planets in 3:1 mean motion resonance. We intend to provide a physical scenario for the formation of HD 60532, which is consistent with the orbital solutions derived from the radial velocity measurements. Observations indicate that the system is in an antisymmetric configuration, while previous theoretical investigations indicate an asymmetric equilibrium state. The paper aims at answering this discrepancy as well. We performed two-dimensional hydrodynamical simulations of thin disks with an embedded pair of massive planets. Additionally, migration and resonant capture are studied by gravitational N-body simulations that apply properly parametrized non-conservative forces. Our simulations suggest that the capture into the 3:1 mean motion resonance takes place only for higher planetary masses, thus favouring orbital solutions having relatively smaller inclination i=20 degrees. The system formed by numerical simulations qualitatively show the same behaviour as HD 60532. We also find that the presence of an inner disk (between the inner planet and the star) plays a very important role in determining the final configurations of resonant planetary systems. Its damping effect on the inner planet's eccentricity is responsible for the observed antisymmetric state of HD 60532.

The size distribution of the stability region around the Lagrangian point L4 is investigated in the elliptic restricted three-body problem as the function of the mass parameter and the orbital eccentricity of the primaries. It is shown that there are minimum zones in the size distribution of the stability regions, and these zones are connected with secondary resonances between the frequencies of librational motions around L4. The results can be applied to hypothetical Trojan planets for predicting values of the mass parameter and the eccentricity for which such objects can be expected or their existence is less probable.

Based on radial velocity measurements it has been found recently that the two giant planets detected around the star HD 73526 are in 2:1 resonance. However, as our numerical integration shows, the derived orbital data for this system result in chaotic behavior of the giant planets, which is uncommon among the resonant extrasolar planetary systems. We intend to present regular (non-chaotic) orbital solutions for the giant planets in the system HD 73526 and offer formation scenarios based on combining planetary migration and sudden perturbative effects such as planet-planet scattering or rapid dispersal of the protoplanetary disk. A comparison with the already studied resonant system HD 128311, exhibiting similar behavior, is also done. The new sets of orbital solutions have been derived by the Systemic Console (www.oklo.org). The stability of these solutions has been investigated by the Relative Lyapunov indicator, while the migration and scattering effects are studied by gravitational N-body simulations applying non-conservative forces as well. Additionally, hydrodynamic simulations of embedded planets in protoplanetary disks are performed to follow the capture into resonance. For the system HD 73526 we demonstrate that the observational radial velocity data are consistent with a coplanar planetary system engaged in a stable 2:1 resonance exhibiting apsidal corotation. We have shown that, similarly to the system HD 128311, the present dynamical state of HD 73526 could be the result of a mixed evolutionary process melting together planetary migration and a perturbative event.

The probability of the detection of Earth-like exoplanets may increase in the near future after the launch of the space missions using the transit photometry as observation method. By using this technique only the semi-major axis of the detected planet can be determined, and there will be no information on the upper limit of its orbital eccentricity. However, the orbital eccentricity is a very important parameter, not only from a dynamical point of view, since it gives also information on the climate and the habitability of the Earth-like planets. In this paper a possible procedure is suggested for confining the eccentricity of an exoplanet discovered by transit photometry if an already known giant planet orbits also in the system.

A significant number of the known multiple exoplanetary systems are containing a pair of giant planets engaged in a low order mean motion resonance. Such a resonant condition protects the dynamics of these planets resulting in very stable orbits. According to recent studies the capture into a resonance is the result of a planetary migration process induced by the interaction of the planets with a protoplanetary disk. If the migration is slow enough (adiabatic) next to a mean motion resonance, the two planets will also be in apsidal corotation. The recently refined orbital parameters of the system HD 128311 suggest that the two giant planets are in a 2:1 mean motion resonance, however without exhibiting apsidal corotation. Thus the evolution of this system can not be described by an adiabatic migration process alone. We present possible evolution scenarios of this system combining migration processes and sudden perturbations. We model migration scenarios through numerical integration of the gravitational N-body problem with additional non-conservative forces. Planet-planet scattering has been investigated by N-body simulations. We show that the present dynamical state of the system HD128311 may be explained by such evolutionary processes.