Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$\mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $\alpha \simeq 10^{-4}$. The value of $\alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $\alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.

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 New Horizons flyby of the cold classical Kuiper Belt object MU69 showed it to be a contact binary. The existence of other contact binaries in the 1-10km range raises the question of how common these bodies are and how they evolved into contact. Here we consider that the lobes of MU69 formed as a binary in the Solar nebula, and calculate its orbital evolution in the presence of gas. We find that the sub-Keplerian wind of the disk brings the drag timescales for 10km bodies under 1 Myr for quadratic-velocity drag. In the Kuiper belt, however, the drag is linear with velocity and the effect of the wind cancels out as the angular momentum gained in half an orbit is lost in the other half: the drag timescales for 10 km bodies remain over 10 Myr. In this situation we find that a combination of drag and Kozai-Lidov oscillations is a promising channel for collapse. We analytically solve the hierarchical three-body problem with drag and implement it into a Kozai plus tidal friction model. The permanent quadrupoles of the lobes make the Kozai oscillations stochastic, and as drag shrinks the semimajor axis it more easily allows the fluctuations to bring the system into contact. Evolution to contact happens very rapidly (within $10^4$ yr) in the pure, double-average quadrupole, Kozai region between $\approx 85-95^\circ$, and within 3 Myr in the drag-assisted region beyond it. The synergy between $J_2$ and gas drag widens the window of contact to 80-100$^\circ$ initial inclination, over a larger range of semimajor axes. As such, the model predicts an initial contact binary fraction of about 10\% for the cold classicals in the Kuiper belt. The speed at contact deviates from the escape velocity only because of the oblateness. For MU69, the contact velocity should be 3.3-4.2 m/s, in line with the observational evidence from the lack of deformation features and estimate of the tensile strength.

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.

Accretion disks around supermassive black holes (SMBH) are promising sites for stellar mass black hole (BH) mergers due to mass segregation and merger acceleration by disk gas torques. Here we show that a GW-kick at BH merger causes ram-pressure stripping of gas within the BH Hill sphere. If $R_{H}\geq H$, the disk height, an off-center UV flare at $a_{\rm BH} \sim 10^{3}r_{g}$ emerges within $t_{\rm UV} \sim \rm{O}(2 {\rm days})(a_{\rm BH}/10^{3}r_{g})(M_{\rm SMBH}/10^{8}M_{\odot})(v_{\rm kick}/10^{2}\rm{km/s})$ post-merger and lasts O$(R_{H}/v_{\rm kick}) \sim \rm{O}(5 t_{\rm UV}$). The flare emerges with luminosity O($10^{42}{\rm erg/s})(t_{\rm UV}/2{\rm days})^{-1}(M_{\rm Hill}/1M_{\odot})(v_{\rm kick}/10^{2}{\rm km/s})^{2}$. AGN optical/UV photometry alters and asymmetric broad emission line profiles can develop after weeks. If $R_{H}<H$, detectability depends on disk optical depth. Follow-up by large optical sky surveys is optimized for small GW error volumes and for LIGO/Virgo triggers $>50M_{\odot}$.

Structures in circumstellar disks such as gaps and rings are often attributed to planets. This connection has been difficult to show unequivocally, as other processes may also produce these features. Particularly, a photoelectric instability (PEI) has been proposed, operating in gas-rich optically thin disks, that generates structures predicted by planet--disk interactions. We examine the question of how to disentangle planetary effects on disk structure from the effects of the PEI. We use the Pencil Code to perform 2D global hydrodynamical models of the dynamics of gas and dust in a thin disk, with and without planetary perturbers. Photoelectric heating is modeled with an equation of state where pressure is proportional to dust surface density. The drag force on grains and its backreaction on the gas are included. Analyzing the situation without PEI, we find that gas--dust interactions alter the shape of the planetary gap from the dust-free case when the local dust-to-gas ratio $\varepsilon$ approaches unity. This result applies also to primordial disks, because dust drifting inwards accumulates at the edge of the planetary gap, and any initial dust-to-gas ratio eventually achieves $\varepsilon=1$ if the dust reservoir is sufficient. We find a result particular to high dust-to-gas ratio disks as well: as dust drifts inwards, the dust front becomes a sharp transition, and the backreaction triggers the Rossby wave instability. When PEI is included, we find that it obscures structures induced by planets unless the planet's mass is sufficiently large to carve a noticeable gap. Specifically, the instability generates arcs and rings of regular spacing: a planet is discernible when it carves a dust gap wider than the wavelength of the PEI.

Two studies utilizing sparse aperture masking (SAM) interferometry and $H_{\rm \alpha}$ differential imaging have reported multiple jovian companions around the young solar-mass star, LkCa 15 (LkCa 15 bcd): the first claimed direct detection of infant, newly-formed planets ("protoplanets"). We present new near-infrared direct imaging/spectroscopy from the SCExAO system coupled with the CHARIS integral field spectrograph and multi-epoch thermal infrared imaging from Keck/NIRC2 of LkCa 15 at high Strehl ratios. These data provide the first direct imaging look at the same wavelengths and in the same locations where previous studies identified the LkCa 15 protoplanets and thus offer the first decisive test of their existence. The data do not reveal these planets. Instead, we resolve extended emission tracing a dust disk with a brightness and location comparable to that claimed for LkCa 15 bcd. Forward-models attributing this signal to orbiting planets are inconsistent with the combined SCExAO/CHARIS and Keck/NIRC2 data. An inner disk provides a more compelling explanation for the SAM detections and perhaps also the claimed $H_{\alpha}$ detection of LkCa 15 b. We conclude that there is currently no clear, direct evidence for multiple protoplanets orbiting LkCa 15, although the system likely contains at least one unseen jovian companion. To identify jovian companions around LkCa 15 from future observations, the inner disk should be detected and its effect modeled, removed, and shown to be distinguishable from planets. Protoplanet candidates identified from similar systems should likewise be clearly distinguished from disk emission through modeling.

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities.

Evidence for exocomets, icy bodies in extrasolar planetary systems, has rapidly increased over the past decade. Volatiles are detected through the gas that exocomets release as they collide and grind down within their natal belts, or as they sublimate once scattered inwards to the regions closest to their host star. Most detections are in young, 10 to a few 100 Myr-old systems that are undergoing the final stages of terrestrial planet formation. This opens the exciting possibility to study exocomets at the epoch of volatile delivery to the inner regions of planetary systems. Detection of molecular and atomic gas in exocometary belts allows us to estimate molecular ice abundances and overall elemental abundances, enabling comparison with the Solar Nebula and Solar System comets. At the same time, observing star-grazing exocomets transiting in front of their star (for planetary systems viewed edge-on) and exozodiacal dust in the systems' innermost regions gives unique dynamical insights into the inward scattering process producing delivery to inner rocky planets. The rapid advances of this budding subfield of exoplanetary science will continue in the short term with the upcoming JWST, WFIRST and PLATO missions. In the longer term, the priority should be to explore the full composition of exocomets, including species crucial for delivery and later prebiotic synthesis. Doing so around an increasingly large population of exoplanetary systems is equally important, to enable comparative studies of young exocomets at the epoch of volatile delivery. We identify the proposed LUVOIR and Origins flagship missions as the most promising for a large-scale exploration of exocometary gas, a crucial component of the chemical heritage of young exo-Earths.

Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength.

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to support liquid water oceans, and thus Earth-like biology, despite differing in many ways from our own planet. In addition to their quantitative abundance, super-Earths are relatively large and are thus more easily detected than true Earth twins. As a result, super-Earths represent a uniquely powerful opportunity to discover and explore a panoply of fascinating and potentially habitable planets in 2020 - 2030 and beyond.

Ground-based telescopes coupled with adaptive optics (AO) have been playing a leading role in exoplanet direct imaging science and technological development for the past two decades and will continue to have an indispensable role for the next decade and beyond. Over the next decade, extreme AO systems on 8-10m telescopes will 1) mitigate risk for WFIRST-CGI by identifying numerous planets the mission can spectrally characterize, 2) validate performance requirements and motivate improvements to atmosphere models needed to unambiguously characterize solar system-analogues from space, and 3) mature novel technological innovations useful for space. Extremely Large Telescopes can deliver the first thermal infrared (10 $\mu m$) images of rocky planets around Sun-like stars and identify biomarkers. These data provide a future NASA direct imaging flagship mission (i.e. HabEx, LUVOIR) with numerous exo-Earth candidates and critical ancillary information to help clarify whether these planets are habitable.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We present simulations of the capabilities of the Next Generation Very Large Array to image at high angular resolution substructures in the dust emission of protoplanetary disks. The main goal of this study is to investigate the kinds of substructures that are expected by state-of-the-art 3D simulations of disks and that an instrument like the ngVLA, with its current design, can detect. The disk simulations adopted in this investigation consist of global 3D radiation-hydrodynamics models with embedded particles, the latter representing dust grains. Our work shows that the ngVLA can detect and spatially resolve, down to sub-astronomical unit scales in disks in nearby star forming regions, the dust continuum emission at 3mm from azimuthal asymmetric structures, as well as from weak rings and gaps produced in these models as a consequence of the vertical shear instability (VSI). This hydrodynamical instability has been proposed to generate turbulence in regions of weak coupling between the disk gas and magnetic field, as well as to form vortices which may be preferred locations of planetesimal formation.

The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates followup VLA observations. Deep multiconfiguration integrations reveal the morphology of Clump 1 and additional cm-wave components which we characterize via comparison with a deconvolution of recent 342GHz data (~1mm). Clump1, which concentrates ~1/3 of the whole disk flux density at ~1cm, is resolved as a narrow arc with a deprojected aspect ratio Chi>5.6, and with half the azimuthal width than at 342 GHz. The spectral trends in the morphology of Clump1 are quantitatively consistent with the Lyra-Lin prescriptions for dust trapping in an anticyclonic vortex, provided with porous grains (f~0.2+-0.2) in a very elongated (Chi~14+-3) and cold (T~23+-2K) vortex. The same prescriptions constrain the turbulence parameter alpha and the gas surface density Sigma_g through log10( alpha x Sigma_g /g/cm2)~-2.3+-0.4, thus requiring values for Sigma_g larger than a factor of a few compared to that reported in the literature from the CO isotopologues, if alpha <~ 1E-3. Such physical conditions imply an appreciably optically thick continuum even at cm-wavelengths (tau(33GHz)~0.2). A secondary and shallower peak at 342GHz is about twice fainter relative to Clump1 at 33GHz. Clump2 appears to be less efficient at trapping large grains.

This review examines recent theoretical developments in our understanding of turbulence in cold, non-magnetically active, planetesimal forming regions of protoplanetary disks which we refer to throughout as "Ohmic zones". We give a brief background introduction to the subject of disk turbulence followed by a terse pedagogical review of the phenomenology of hydrodynamic turbulence. The equations governing the dynamics of cold astrophysical disks are given and basic flow states are described. We discuss the Solberg-Hoiland conditions required for stability, and the three recently identified turbulence generating mechanisms possibly active in protoplanetary disk Ohmic zones, namely, (i) the Vertical Shear Instability, (ii) The Convective Overstability and (iii) the Zombie Vortex Instability. We summarize the properties of these processes, identify their limitations and discuss where and under what conditions these processes are active in protoplanetary disk Ohmic zones.

The formation of planets remains one of the most challenging problems of contemporary astrophysics. Starting with micron-sized dust grains, coagulation models predict growth up to centimeter (pebbles), but growth beyond this size is difficult because of fragmentation and drift. Ways to bypass this problem have focused on inhomogeneities in the flow, be that zonal flows, streaming instability, or vortices. Because vortices are in equilibrium between the Coriolis and the pressure force, the pressureless grains will orbit along a vortex streamline experiencing a drag force. This is a very effective mechanism to concentrate pebbles as also seen in numerical simulations and possibly in ALMA observations. Yet, a high pebble load is dangerous for the vortex, and we showed that in two-dimensional simulations the backreaction eventually leads to vortex disruption. We investigate whether the same happens in three dimensions. We perform 3D simulations with pebbles in a local box finding that, although the pebbles disturb the vortex around the midplane, the column does not get destroyed. This result is important because, based on the previous 2D result suggesting complete disruption, the vortex interpretation of ALMA observations has been called into question. We show instead that the vortex behaves like a Taylor column, and the pebbles as obstacles to the flow. Pebble accumulation in the center of the vortices proceeds to roughly the same concentration as in the control run without backreaction.

WD 1145+017 is currently the only white dwarf known to exhibit periodic transits of planetary debris as well as absorption lines from circumstellar gas. We present the first simultaneous fast optical spectrophotometry and broad-band photometry of the system, obtained with the Gran Telescopio Canarias (GTC) and the Liverpool Telescope (LT), respectively. The observations spanned $5.5$ h, somewhat longer than the $4.5$-h orbital period of the debris. Dividing the GTC spectrophotometry into five wavelength bands reveals no significant colour differences, confirming grey transits in the optical. We argue that absorption by an optically thick structure is a plausible alternative explanation for the achromatic nature of the transits that can allow the presence of small-sized ($\sim\mu$m) particles. The longest ($87$ min) and deepest ($50$ per cent attenuation) transit recorded in our data exhibits a complex structure around minimum light that can be well modelled by multiple overlapping dust clouds. The strongest circumstellar absorption line, Fe II $\lambda$5169, significantly weakens during this transit, with its equivalent width reducing from a mean out-of-transit value of $2$ \AA\ to $1$ \AA\ in-transit, supporting spatial correlation between the circumstellar gas and dust. Finally, we made use of the Gaia Data Release 2 and archival photometry to determine the white dwarf parameters. Adopting a helium-dominated atmosphere containing traces of hydrogen and metals, and a reddening $E(B-V)=0.01$ we find $T_\mathrm{eff}=15\,020 \pm 520$ K, $\log g=8.07\pm0.07$, corresponding to $M_\mathrm{WD}=0.63\pm0.05\ \mbox{$\mathrm{M}_{\odot}$}$ and a cooling age of $224\pm30$ Myr.

Recent developments in non-ideal magnetohydrodynamic simulations of protoplanetary disks suggest that instead of being traditional turbulent (viscous) accretion disks, they have a largely laminar flow with accretion driven by large-scale wind torques. These disks are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields. We have examined the dynamics of the corotation region of a low mass planet embedded in such a disk and the evolution of the associated migration torque. These disks lack strong turbulence and associated turbulent diffusion, and the presence of a magnetic field and radial gas flow presents a situation outside the applicability of previous corotation torque theory. We summarize the analytical analysis of the corotation torque, give details on the numerical methods used, and in particular the relative merits of different numerical schemes for the inviscid problem.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

Black hole mergers detectable with LIGO can occur in active galactic nucleus (AGN) disks. Here we parameterize the merger rates, the mass spectrum and the spin spectrum of black holes (BH) in AGN disks. The predicted merger rate spans $\sim 10^{-4}-10^{4} \rm{Gpc}^{-1} \rm{yr}^{-1}$, so upper limits from LIGO ($<212\rm{Gpc}^{-1}\rm{yr}^{-1}$) already constrain it. The predicted mass spectrum has the form of a broken power-law consisting of a pre-existing BH powerlaw mass spectrum and a harder powerlaw mass spectrum resulting from mergers. The predicted spin spectrum is multi-peaked with the evolution of retrograde spin BH in the gas disk playing a key role. We outline the large uncertainties in each of these LIGO observables for this channel and we discuss ways in which they can be constrained in the future.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.

Planets are born in protostellar disks, which are now observed with enough resolution to address questions about internal gas flows. Candidates for driving the flows include magnetic forces, but ionization state estimates suggest much of the gas mass decouples from magnetic fields. Thus, hydrodynamical instabilities could play a major role. We investigate disk dynamics under conditions typical for a T Tauri system, using global 3D radiation hydrodynamics simulations with embedded particles and a resolution of 70 cells per scale height. Stellar irradiation heating is included with realistic dust opacities. The disk starts in joint radiative balance and hydrostatic equilibrium. The vertical shear instability (VSI) develops into turbulence that persists up to at least 1600 inner orbits (143 outer orbits). Turbulent speeds are a few percent of the local sound speed at the midplane, increasing to 20%, or 100 m/s, in the corona. These are consistent with recent upper limits on turbulent speeds from optically thin and thick molecular line observations of TW Hya and HD 163296. The predominantly vertical motions induced by the VSI efficiently lift particles upwards. Grains 0.1 and 1 mm in size achieve scale heights greater than expected in isotropic turbulence. We conclude that while kinematic constraints from molecular line emission do not directly discriminate between magnetic and nonmagnetic disk models, the small dust scale heights measured in HL Tau and HD 163296 favor turbulent magnetic models, which reach lower ratios of the vertical kinetic energy density to the accretion stress.

Previous theoretical works have shown that in optically thin disks, dust grains are photoelectrically stripped of electrons by starlight, heating nearby gas and possibly creating a dust clumping instability, the photoelectric instability (PeI), that significantly alters global disk structure. In the current work, we use the Pencil Code to perform the first numerical models of the PeI that include stellar radiation pressure on dust grains in order to explore the parameter regime in which the instability operates. In models with gas surface densities greater than $\sim$$10^{-4}~\mathrm{g}~\mathrm{cm}^{-2}$, we see a variety of dust structures, including sharp concentric rings and non-axisymmetric arcs and clumps that represent dust surface density enhancements of factors of $\sim$$5-20$ depending on the run parameters. The gas distributions show various structures as well, including clumps and arcs formed from spiral arms. In models with lower gas surface densities, vortices and smooth spiral arms form in the gas distribution, but the dust is too weakly coupled to the gas to be significantly perturbed. In one high gas surface density model, we include a large, low-order gas viscosity, and, in agreement with previous radiation pressure-free models, find that it observably smooths the structures that form in the gas and dust, suggesting that resolved images of a given disk may be useful for deriving constraints on the effective viscosity of its gas. Broadly, our models show that radiation pressure does not preclude the formation of complex structure from the PeI, but the qualitative manifestation of the PeI depends strongly on the parameters of the system. The PeI may provide an explanation for unusual disk morphologies such as the moving blobs of the AU Mic disk, the asymmetric dust distribution of the 49 Ceti disk, and the rings and arcs found in the disk around HD 141569A.

Recent observations of the protoplanetary disk around the Herbig Be star HD 100546 show two bright features in infrared (H and L' bands) at about 50 AU, with one so far unexplained. We explore the observational signatures of a high mass planet causing shock heating in order to determine if it could be the source of the unexplained infrared feature in HD 100546. More fundamentally, we identify and characterize planetary shocks as an extra, hitherto ignored, source of luminosity in transition disks. The RADMC-3D code is used to perform dust radiative transfer calculations on the hydrodynamical disk models, including volumetric heating. A stronger shock heating rate by a factor 20 would be necessary to qualitatively reproduce the morphology of the second infrared source. Instead, we find that the outer edge of the gap carved by the planet heats up by about 50% relative to the initial reference temperature, which leads to an increase in the scale height. The bulge is illuminated by the central star, producing a lopsided feature in scattered light, as the outer gap edge shows an asymmetry in density and temperature attributable to a secondary spiral arm launched not from the Lindblad resonances but from the 2:1 resonance. We conclude that high-mass planets lead to shocks in disks that may be directly observed, particularly at wavelengths of 10 mm or longer, but that they are more likely to reveal their presence in scattered light by puffing up their outer gap edges and exciting multiple spiral arms.

Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear term, and we find that we do not need an advective gauge as solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited for non-magnetic problems. Benchmarked results of (hydrodynamical) planet-disk interaction and the of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficiently in straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community.

Motivated by models suggesting that the inner planet forming regions of protoplanetary discs are predominantly lacking in viscosity-inducing turbulence, and are possibly threaded by Hall-effect generated large-scale horizontal magnetic fields, we examine the dynamics of the corotation region of a low-mass planet in such an environment. The corotation torque in an inviscid, isothermal, dead zone ought to saturate, with the libration region becoming both symmetrical and of a uniform vortensity, leading to fast inward migration driven by the Lindblad torques alone. However, in such a low viscosity situation, the material on librating streamlines essentially preserves its vortensity. If there is relative radial motion between the disc gas and the planet, the librating streamlines will no longer be symmetrical. Hence, if the gas is torqued by a large scale magnetic field so that it undergoes a net inflow or outflow past the planet, driving evolution of the vortensity and inducing asymmetry of the corotation region, the corotation torque can grow, leading to a positive torque. In this paper we treat this effect by applying a symmetry argument to the previously studied case of a migrating planet in an inviscid disc. Our results show that the corotation torque due to a laminar Hall-induced magnetic field in a dead zone behaves quite differently from that studied previously for a viscous disc. Furthermore, the magnetic field induced corotation torque and the dynamical corotation torque in a low viscosity disc can be regarded as one unified effect.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted decades ago by numerical simulations of disk-planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for the observed signatures. However, such interpretation is not free of problems. The spirals are found to have large pitch angles, and in at least one case (HD 100546) the spiral feature appears effectively unpolarized, which implies thermal emission of the order of 1000K (465$\pm$40K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. In this paper we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5$M_J$ planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient disk gas at that radius. The gas is accelerated vertically away from the midplane by the shocks to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf near the Lindblad resonances. This turbulence, although localized, has high $\alpha$ values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We also find evidence that the disk regions heated up by the planetary shocks eventually becomes superadiabatic, generating convection far from the planet's orbit.

The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a compact emission region in the outer disk indicating a strong concentration of big dust grains. Tracing smaller grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to the north-west of the central star, which matches the VLA clump position. The compactness of the Ka emission is expected in the context of dust trapping, as big grains are trapped more easily than smaller grains in gas pressure maxima. We develop a non-axisymmetric parametric model inspired by a steady state vortex solution with parameters adequately selected to reproduce the observations, including the spectral energy distribution. Finally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially coincides with a spiral-like feature seen in scattered light, while the VLA clump is offset from the scattered light maximum. Moreover, the ALMA map shows a decrement that matches a region devoid of scattered polarised emission. Continuum observations at a different wavelength are necessary to conclude if the VLA-ALMA difference is an opacity or a real dust segregation.

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet-disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1-10 M_jup for a 1 M_sun star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets (M>5 M_jup) produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a 10 M_jup planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet-disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory.

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3<P<11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard "disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous {\em Spitzer} IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

We analyse the concentration of solid particles in vortices created and sustained by radial buoyancy in protoplanetary disks, i.e. baroclinic vortex growth. Besides the gas drag acting on particles we also allow for back-reaction from dust onto the gas. This becomes important when the local dust-to-gas ratio approaches unity. In our 2D, local, shearing sheet simulations we see high concentrations of grains inside the vortices for a broad range of Stokes numbers, ${\rm St}$. An initial dust-to-gas ratio of 1:100 can easily be reversed to 100:1 for ${\rm St}=1$. The increased dust-to-gas ratio triggers the streaming instability, thus counter-intuitively limiting the maximal achievable overdensities. We find that particle trapping inside vortices opens the possibility for gravity-assisted planetesimal formation even for small particles ($\rm{St}=0.01$) and low initial dust-to-gas ratios (1:$10^4$)

Early in the study of viscous accretion disks it was realized that energy transfers from distant sources must be important, not least because the flow at the disk midplane in the bulk of the disk is likely outwards, out of the gravitational potential well. If the source of the viscosity is powered by accretion, such as in the case of the magneto-rotational instability, such distant energy sources must lie in the innermost regions of the disk, where accretion occurs even at the midplane. We argue here that modulations in this energy supply can alter the accretion rate on dynamical, rather than far longer viscous, time scales. This means that both the steady state value of and fluctuations in the inner disk's accretion rate, depending on the details of the inner boundary condition and occurring on the inner disk's rapid evolution time, can affect the outer disk. This is particularly interesting because observations have shown that disk accretion is not steady (e.g.~EX Lupi type objects). We also note that the power supplied to shearing boxes is set by the boxes themselves rather than the physical energy fluxes in a global disk. That is, their saturated magnetic field is not subject to the full set of energy constraints present in an actual disk. Our analysis suggests that large scale radial transport of energy has a critical impact on the evolution and variability of accretion disks.

Rossby wave instability (RWI) at dead zone boundaries may play an important role in planet formation. Viscous hydrodynamics results suggest RWI is excited only when the viscosity changes over a radial distance less than two density scale heights. However in the disks around Solar-mass T Tauri stars, it is not viscosity but magnetic forces that provide the accretion stress beyond about 10 AU, where surface densities are low enough so stellar X-rays and interstellar cosmic rays can penetrate. Here we aim to explore the conditions for RWI in the smooth transition with increasing distance, from resistive and magnetically-dead to conducting and magnetically-active. We perform 3D unstratified MHD simulations with the Pencil Code, using static resistivity profiles. As a result, we find that in MHD, contrary to viscous models, the RWI is triggered even with a gradual change in resistivity extending from 10 to 40 AU (i.e., spanning 15 scale heights for aspect ratio 0.1). This is because magneto-rotational turbulence sets in abruptly when the resistivity reaches a threshold level. At higher resistivities the longest unstable wavelength is quenched, resulting in a sharp decline of the Maxwell stress towards the star. The sharp gradient in the magnetic forces leads to a localized density bump, that is in turn Rossby wave unstable. We conclude that even weak gradients in the resistivity can lead to sharp transitions in the Maxwell stress. The upshot is that the RWI is more easily activated in the outer disk than previously thought. Rossby vortices at the outer dead zone boundary thus could underlie the dust asymmetries seen in the outer reaches of transition disks.

Recently, Klahr & Hubbard (2014) claimed that a hydrodynamical linear overstability exists in protoplanetary disks, powered by buoyancy in the presence of thermal relaxation. We analyse this claim, confirming it through rigorous compressible linear analysis. We model the system numerically, reproducing the linear growth rate for all cases studied. We also study the saturated properties of the overstability in the shearing box, finding that the saturated state produces finite amplitude fluctuations strong enough to trigger the subcritical baroclinic instability. Saturation leads to a fast burst of enstrophy in the box, and a large-scale vortex develops in the course of the next $\approx$100 orbits. The amount of angular momentum transport achieved is of the order of $\alpha \approx 10^{-3}$, as in compressible SBI models. For the first time, a self-sustained 3D vortex is produced from linear amplitude perturbation of a quiescent base state.

If intermediate mass black holes (IMBHs) grow efficiently in gas disks around supermassive black holes, their host active galactic nucleus (AGN) disks should exhibit myriad observational signatures. Gap-opening IMBHs in AGN disks can exhibit spectral features and variability analagous to gapped protoplanetary disks. A gap-opening IMBH in the innermost disk imprints ripples and oscillations on the broad Fe K$\alpha$ line which may be detectable with future X-ray missions. A non-gap-opening IMBH will accrete and produce a soft X-ray excess relative to continuum emission. An IMBH on a retrograde orbit in an AGN disk will not open a gap and will generate soft X-rays from a bow-shock 'headwind'. Accreting IMBH in a large cavity can generate ULX-like X-ray luminosities and LINER-like optical line ratios from local ionized gas. We propose that many LINERs house a weakly accreting MBH binary in a large central disk cavity and will be luminous sources of gravitational waves (GW). IMBHs in galactic nuclei may also be detected via intermittent observational signatures including: UV/X-ray flares due to tidal disruption events, asymmetric X-ray intensity distributions as revealed by AGN transits, quasi-periodic oscillations and underluminous Type Ia supernovae. GW emitted during IMBH inspiral and collisions may be detected with eLISA and LIGO, particularly from LINERs. We summarize observational signatures and compare to current data where possible or suggest future observations.

Debris disks should not be completely gas-free, since there is second generation gas from outgassing of planetesimals and dust grains via sublimation, photodesorption, or collisions, generating a system of dust-to-gas ratio close to unity, where hydrodynamics cannot be ignored. A clumping instability exists in this configuration, that has been hitherto explored only in one-dimensional, incompressible models. We performed 2D numerical compressible models of a disk with comparable amounts of gas and dust to study the growth and development of this instability. Our model solves the momentum equation for the gas and dust, together with energy and continuity equations. We uncover that the backreaction of the drag force from the gas onto the dust shepherds rings, similar to those observed in debris disks and usually attributed to the presence of hypothetical undetected planets. We also uncover that the eccentricity of these rings, usually presented as convincing evidence for the presence of a planet, can actually be simply explained by a standing wave propagating along the ring. The rings support a spectrum of oscillations, with one particular mode representing epicyclic motion. The apparent eccentricity matches the eccentricity in observed systems. This suggests that the planet possibility, though thrilling, is not necessarily required to explain these systems.

"Debris disks" around young stars (analogues of the Kuiper Belt in our Solar System) show a variety of non-trivial structures attributed to planetary perturbations and used to constrain the properties of the planets. However, these analyses have largely ignored the fact that some debris disks are found to contain small quantities of gas, a component that all such disks should contain at some level. Several debris disks have been measured with a dust-to-gas ratio around unity at which the effect of hydrodynamics on the structure of the disk cannot be ignored. Here we report linear and nonlinear modelling that shows that dust-gas interactions can produce some of the key patterns attributed to planets. We find a robust clumping instability that organizes the dust into narrow, eccentric rings, similar to the Fomalhaut debris disk. The conclusion that such disks might contain planets is not necessarily required to explain these systems.

The Atacama Large Millimeter Array (ALMA) has been returning images of transitional disks in which large asymmetries are seen in the distribution of mm-sized dust in the outer disk. The explanation in vogue borrows from the vortex literature by suggesting that these asymmetries are the result of dust trapping in giant vortices, excited via Rossby wave instability (RWI) at planetary gap edges. Due to the drag force, dust trapped in vortices will accumulate in the center, and diffusion is needed to maintain a steady state over the lifetime of the disk. While previous work derived semi-analytical models of the process, in this paper we provide analytical steady-state solutions. Exact solutions exist for certain vortex models. The solution is determined by the vortex rotation profile, the gas scale height, the vortex aspect ratio, and the ratio of dust diffusion to gas-dust friction. In principle, all these quantities can be derived from observations, which would give validation of the model, also giving constrains on the strength of the turbulence inside the vortex core. Based on our solution, we derive quantities such as the gas-dust contrast, the trapped dust mass, and the dust contrast at the same orbital location. We apply our model to the recently imaged Oph IRS 48 system, finding values within the range of the observational uncertainties.

Vortices are the fundamental units of turbulent flow. Understanding their stability properties therefore provides fundamental insights on the nature of turbulence itself. In this contribution I briefly review the phenomenological aspects of the instability of elliptic streamlines, in the hydro (elliptic instability) and hydromagnetic (magneto-elliptic instability) regimes. Vortex survival in disks is a balance between vortex destruction by these mechanisms, and vortex production by others, namely, the Rossby wave instability and the baroclinic instability.

Recent studies have shown that baroclinic vortex amplification is strongly dependent on certain factors, namely, the global entropy gradient, the efficiency of thermal diffusion and/or relaxation as well as numerical resolution. We conduct a comprehensive study of a broad range and combination of various entropy gradients, thermal diffusion and thermal relaxation time-scales via local shearing sheet simulations covering the parameter space relevant for protoplanetary disks. We measure the Reynolds stresses as a function of our control parameters and see that there is angular momentum transport even for entropy gradients as low as $\beta=-{d\ln s}/{d\ln r}={1}/{2}$, which corresponds to values observed in protoplanetary accretion disks. The amplification-rate of the perturbations, $\Gamma$, appears to be proportional to $\beta^2$ and thus proportional to the square of the \BV ($\Gamma \propto \beta^2 \propto N^2$). The saturation level of Reynolds stresses on the other hand seems to be proportional to $\beta^{1/2}$. This highlights the importance of baroclinic effects even for the low entropy gradients expected in protoplanetary disks.

It has been suggested that the transition between magnetorotationally active and dead zones in protoplanetary disks should be prone to the excitation of vortices via Rossby wave instability (RWI). However, the only numerical evidence for this has come from alpha disk models, where the magnetic field evolution is not followed, and the effect of turbulence is parametrized by Laplacian viscosity. We aim to establish the phenomenology of the flow in the transition in 3D resistive-magnetohydrodynamical models. We model the transition by a sharp jump in resistivity, as expected in the inner dead zone boundary, using the Pencil Code to simulate the flow. We find that vortices are readily excited in the dead side of the transition. We measure the mass accretion rate finding similar levels of Reynolds stress at the dead and active zones, at the $\alpha\approx 10^{-2}$ level. The vortex sits in a pressure maximum and does not migrate, surviving until the end of the simulation. A pressure maximum in the active zone also triggers the RWI. The magnetized vortex that results should be disrupted by parasitical magneto-elliptic instabilities, yet it subsists in high resolution. This suggests that either the parasitic modes are still numerically damped, or that the RWI supplies vorticity faster than they can destroy it. We conclude that the resistive transition between the active and dead zones in the inner regions of protoplanetary disks, if sharp enough, can indeed excite vortices via RWI. Our results lend credence to previous works that relied on the alpha-disk approximation, and caution against the use of overly reduced azimuthal coverage on modeling this transition.

Here we propose a mechanism for efficiently growing intermediate mass black holes (IMBH) in disks around supermassive black holes. Stellar mass objects can efficiently agglomerate when facilitated by the gas disk. Stars, compact objects and binaries can migrate, accrete and merge within disks around supermassive black holes. While dynamical heating by cusp stars excites the velocity dispersion of nuclear cluster objects (NCOs) in the disk, gas in the disk damps NCO orbits. If gas damping dominates, NCOs remain in the disk with circularized orbits and large collision cross-sections. IMBH seeds can grow extremely rapidly by collisions with disk NCOs at low relative velocities, allowing for super-Eddington growth rates. Once an IMBH seed has cleared out its feeding zone of disk NCOs, growth of IMBH seeds can become dominated by gas accretion from the AGN disk. However, the IMBH can migrate in the disk and expand its feeding zone, permitting a super-Eddington accretion rate to continue. Growth of IMBH seeds via NCO collisions is enhanced by a pile-up of migrators. We highlight the remarkable parallel between the growth of IMBH in AGN disks with models of giant planet growth in protoplanetary disks. If an IMBH becomes massive enough it can open a gap in the AGN disk. IMBH migration in AGN disks may stall, allowing them to survive the end of the AGN phase and remain in galactic nuclei. Our proposed mechanisms should be more efficient at growing IMBH in AGN disks than the standard model of IMBH growth in stellar clusters. Dynamical heating of disk NCOs by cusp stars is transferred to the gas in a AGN disk helping to maintain the outer disk against gravitational instability. Model predictions, observational constraints and implications are discussed in a companion paper (Paper II).

Recently, there has been a significant level of discussion of the correct treatment of Kelvin-Helmholtz instability in the astrophysical community. This discussion relies largely on how the KHI test is posed and analyzed. We pose a stringent test of the initial growth of the instability. The goal is to provide a rigorous methodology for verifying a code on two dimensional Kelvin-Helmholtz instability. We ran the problem in the Pencil Code, Athena, Enzo, NDSPHMHD, and Phurbas. A strict comparison, judgment, or ranking, between codes is beyond the scope of this work, though this work provides the mathematical framework needed for such a study. Nonetheless, how the test is posed circumvents the issues raised by tests starting from a sharp contact discontinuity yet it still shows the poor performance of Smoothed Particle Hydrodynamics. We then comment on the connection between this behavior to the underlying lack of zeroth-order consistency in Smoothed Particle Hydrodynamics interpolation. We comment on the tendency of some methods, particularly those with very low numerical diffusion, to produce secondary Kelvin-Helmholtz billows on similar tests. Though the lack of a fixed, physical diffusive scale in the Euler equations lies at the root of the issue, we suggest that in some methods an extra diffusion operator should be used to damp the growth of instabilities arising from grid noise. This statement applies particularly to moving-mesh tessellation codes, but also to fixed-grid Godunov schemes.

It has been recently suggested that the magneto-rotational instability (MRI) is a limiting case of the magneto-elliptic instability (MEI). This limit is obtained for horizontal modes in the presence of rotation and an external vertical magnetic field, when the aspect ratio of the elliptic streamlines tends to infinite. In this paper we unveil the link between these previously unconnected mechanisms, explaining both the MEI and the MRI as different manifestations of the same Magneto-Elliptic-Rotational Instability (MERI). The growth rates are found and the influence of the magnetic and rotational effects is explained, in particular the effect of the magnetic field on the range of negative Rossby numbers at which the horizontal instability is excited. Furthermore, we show how the horizontal rotational MEI in the rotating shear flow limit links to the MRI by the use of the local shearing box model, typically used in the study of accretion discs. In such limit the growth rates of the two instability types coincide for any power-type background angular velocity radial profile with negative exponent corresponding to the value of the Rossby number of the rotating shear flow. The MRI requirement for instability is that the background angular velocity profile is a decreasing function of the distance from the centre of the disk which corresponds to the horizontal rotational MEI requirement of negative Rossby numbers. Finally a physical interpretation of the horizontal instability, based on a balance between the strain, the Lorentz force and the Coriolis force is given.

The torques exerted by a locally isothermal disk on an embedded planet lead to rapid inward migration. Recent work has shown that modeling the thermodynamics without the assumption of local isothermality reveals regions where the net torque on an embedded planet is positive, leading to outward migration of the planet. When a region with negative torque lies directly exterior to this, planets in the inner region migrate outwards and planets in the outer region migrate inwards, converging where the torque is zero. We incorporate the torques from an evolving non-isothermal disk into an N-body simulation to examine the behavior of planets or planetary embryos interacting in the convergence zone. We find that mutual interactions do not eject objects from the convergence zone. Small numbers of objects in a laminar disk settle into near resonant orbits that remain stable over the 10 Myr periods that we examine. However, either or both increasing the number of planets or including a correlated, stochastic force to represent turbulence drives orbit crossings and mergers in the convergence zone. These processes can build gas giant cores with masses of order ten Earth masses from sub-Earth mass embryos in 2-3 Myr.

Based on the viscous disk theory, a number of recent studies have suggested there is large scale meridional circulation in protoplanetary disks. Such a flow could account for the presence of crystalline silicates, including calcium- and aluminum-rich inclusions (CAIs), at large distances from the sun. This paper aims at examining whether such large-scale flows exist in turbulent protoplanetary disks. High-resolution global hydrodynamical and magnetohydrodynamical (MHD) numerical simulations of turbulent protoplanetary disks were used to infer the properties of the flow in such disks. By performing hydrodynamic simulations using explicit viscosity, we demonstrate that our numerical setup does not suffer from any numerical artifact. The aforementioned meridional circulation is easily recovered in viscous and laminar disks and is quickly established. In MHD simulations, the magnetorotational instability drives turbulence in the disks. Averaging out the turbulent fluctuations on a long timescale, the results fail to show any large-scale meridional circulation. A detailed analysis of the simulations show that this lack of meridional circulation is due to the turbulent stress tensor having a vertical profile different from the viscous stress tensor. A simple model is provided that successfully accounts for the structure of the flow in the bulk of the disk. In addition to those results, possible deviations from standard vertically averaged alpha disk models are suggested by the simulations and should be the focus of future work. Global MHD numerical simulations of fully ionized and turbulent protoplanetary disks are not consistent with the existence of a large-scale meridional flow. As a consequence, the presence of crystalline silicates at large distance for the central star cannot be accounted for by that process as suggested by recent models based on viscous disk theory.

Galactic nuclei should contain a cluster of stars and compact objects in the vicinity of the central supermassive black hole due to stellar evolution, minor mergers and gravitational dynamical friction. By analogy with protoplanetary migration, nuclear cluster objects (NCOs) can migrate in the accretion disks that power active galactic nuclei by exchanging angular momentum with disk gas. Here we show that an individual NCO undergoing runaway outward migration comparable to Type III protoplanetary migration can generate an accretion rate corresponding to Seyfert AGN or quasar luminosities. Multiple migrating NCOs in an AGN disk can dominate traditional viscous disk accretion and at large disk radii, ensemble NCO migration and accretion could provide sufficient heating to prevent the gravitational instability from consuming disk gas in star formation. The magnitude and energy of the X-ray soft excess observed at ~0.1-1keV in Seyfert AGN could be explained by a small population of ~10^{2}-10^{3} accreting stellar mass black holes or a few ULXs. NCO migration and accretion in AGN disks are therefore extremely important mechanisms to add to realistic models of AGN disks.

Turbulence and angular momentum transport in accretion disks remains a topic of debate. With the realization that dead zones are robust features of protoplanetary disks, the search for hydrodynamical sources of turbulence continues. A possible source is the baroclinic instability (BI), which has been shown to exist in unmagnetized non-barotropic disks. We present shearing box simulations of baroclinicly unstable, magnetized, 3D disks, in order to assess the interplay between the BI and other instabilities, namely the magneto-rotational instability (MRI) and the magneto-elliptical instability. We find that the vortices generated and sustained by the baroclinic instability in the purely hydrodynamical regime do not survive when magnetic fields are included. The MRI by far supersedes the BI in growth rate and strength at saturation. The resulting turbulence is virtually identical to an MRI-only scenario. We measured the intrinsic vorticity profile of the vortex, finding little radial variation in the vortex core. Nevertheless, the core is disrupted by an MHD instability, which we identify with the magneto-elliptic instability. This instability has nearly the same range of unstable wavelengths as the MRI, but has higher growth rates. In fact, we identify the MRI as a limiting case of the magneto-elliptic instability, when the vortex aspect ratio tends to infinity (pure shear flow). We conclude that vortex excitation and self-sustenance by the baroclinic instability in protoplanetary disks is viable only in low ionization, i.e., the dead zone. Our results are thus in accordance with the layered accretion paradigm. A baroclinicly unstable dead zone should be characterized by the presence of large-scale vortices whose cores are elliptically unstable, yet sustained by the baroclinic feedback. As magnetic fields destroy the vortices and the MRI outweighs the BI, the active layers are unmodified.

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.

Outward migration of low-mass planets has recently been shown to be a possibility in non-barotropic disks. We examine the consequences of this result in evolutionary models of protoplanetary disks. Planet migration occurs towards equilibrium radii with zero torque. These radii themselves migrate inwards because of viscous accretion and photoevaporation. We show that as the surface density and temperature fall, the planet orbital migration and disk depletion timescales eventually become comparable, with the precise timing depending on the mass of the planet. When this occurs, the planet decouples from the equilibrium radius. At this time, however, the gas surface density is already too low to drive substantial further migration. A higher mass planet, of 10 Earth masses, can open a gap during the late evolution of the disk, and stops migrating. Low mass planets, with 1 or 0.1 Earth masses, released beyond 1 AU in our models, avoid migrating into the star. Our results provide support for the reduced migration rates adopted in recent planet population synthesis models.

Extrasolar planets are not named and are referred to only by their assigned scientific designation. The reason given by the IAU to not name the planets is that it is considered impractical as planets are expected to be common. I advance some reasons as to why this logic is flawed, and suggest names for the 403 extrasolar planet candidates known as of Oct 2009. The names follow a scheme of association with the constellation that the host star pertains to, and therefore are mostly drawn from Roman-Greek mythology. Other mythologies may also be used given that a suitable association is established.

As accretion in protoplanetary disks is enabled by turbulent viscosity, the border between active and inactive (dead) zones constitutes a location where there is an abrupt change in the accretion flow. The gas accumulation that ensues triggers the Rossby wave instability, that in turn saturates into anticyclonic vortices. It was suggested that the trapping of solids within them leads to a burst of planet formation on very short timescales. We perform two-dimensional global simulations of the dynamics of gas and solids in a non-magnetized thin protoplanetary disk with the Pencil Code. We use multiple particle species of radius 1, 10, 30, and 100 cm, solving for the particles' gravitational interaction by a particle-mesh method. The dead zone is modeled as a region of low viscosity. Adiabatic and locally isothermal equations of state are used. We find that the Rossby wave instability is triggered under a variety of conditions, thus making vortex formation a robust process. Inside the vortices, fast accumulation of solids occurs and the particles collapse into objects of planetary mass in timescales as short as five orbits. Because the drag force is size-dependent, aerodynamical sorting ensues within the vortical motion, and the first bound structures formed are composed primarily of similarly-sized particles. In addition to erosion due to ram pressure, we identify gas tides from the massive vortices as a disrupting agent of formed protoplanetary embryos. We also estimate the collisional velocity history of the particles that compose the most massive embryo by the end of the simulation, finding that the vast majority of them never experienced a collision with another particle at speeds faster than 1 m/s.

A calibration of H-alpha as both a chromospheric diagnostic and an age indicator is presented, complementing the works previously done on this subject (Herbig 1985, Pasquini & Pallavicini 1991. The chromospheric diagnostic was built with a statistically significant sample, covering nine years of observations, and including 175 solar neighborhood stars. Regarding the age indicator, the presence of stars for which very accurate ages are determined, such as those belonging to clusters and kinematic groups, lends confidence to our analysis. We also investigate the possibility that stars of the same age might have gone through different tracks of chromospheric decay, identifying - within the same age range - effects of metallicity and mass. These parameters, however, as well as age, seem to be significant only for dwarf stars, losing their meaning when we analyze stars in the subgiant branch. This result suggests that, in these evolved stars, the emission mechanism cannot be magnetohydrodynamical in nature, in agreement with recent models (Fawzy et al. 2002c, and references therein). The Sun is found to be a typical star in its H-alpha chromospheric flux, for its age, mass and metallicity. As a byproduct of this work, we developed an automatic method to determine temperatures from the wings of H-alpha, which means the suppression of the error inherent to the visual procedure used in the literature.

We present local numerical models of accretion disk turbulence driven by the magnetorotational instability with varying shear rate. The resulting turbulent stresses are compared with predictions of a closure model in which triple correlations are modelled in terms of quadratic correlations. This local model uses five nondimensional parameters to describe the properties of the flow. We attempt to determine these closure parameters for our simulations and find that the model does produce qualitatively correct behaviour. In addition, we present results concerning the shear rate dependency of the magnetic to kinetic energy ratio. We find both the turbulent stress ratio and the total stress to be strongly dependent on the shear rate.

Centimeter and meter sized solid particles in protoplanetary disks are trapped within long lived high pressure regions, creating opportunities for collapse into planetesimals and planetary embryos. We study the accumulations in the stable Lagrangian points of a giant planet, as well as in the Rossby vortices launched at the edges of the gap it carves. We employ the Pencil Code, tracing the solids with a large number of interacting Lagrangian particles, usually 100,000. For particles of 1 cm to 10 cm radii, gravitational collapse occurs in the Lagrangian points in less than 200 orbits. For 5 cm particles, a 2 Earth mass planet is formed. For 10 cm, the final maximum collapsed mass is around 3 Earth masses. The collapse of the 1 cm particles is indirect, following the timescale of depletion of gas from the tadpole orbits. In the edges of the gap vortices are excited, trapping preferentially particles of 30 cm radii. The rocky planet that is formed is as massive as 17 Earth masses, constituting a Super-Earth. By using multiple particle species, we find that gas drag modifies the streamlines in the tadpole region around the classical L4 and L5 points. As a result, particles of different radii have their stable points shifted to different locations. Collapse therefore takes longer and produces planets of lower mass. Three super-Earths are formed in the vortices, the most massive having 4.4 Earth masses. We conclude that a Jupiter mass planet can induce the formation of other planetary embryos in the outer edge of its gas gap. Trojan Earth mass planets are readily formed, and although not existing in the solar system, might be common in the exoplanetary zoo.

In the borders of the dead zones of protoplanetary disks, the inflow of gas produces a local density maximum that triggers the Rossby wave instability. The vortices that form are efficient in trapping solids. We aim to assess the possibility of gravitational collapse of the solids within the Rossby vortices. We perform global simulations of the dynamics of gas and solids in a low mass non-magnetized self-gravitating thin protoplanetary disk with the Pencil code. We use multiple particle species of radius 1, 10, 30, and 100 cm. The dead zone is modeled as a region of low viscosity. The Rossby vortices excited in the edges of the dead zone are very efficient particle traps. Within 5 orbits after their appearance, the solids achieve critical density and undergo gravitational collapse into Mars sized objects. The velocity dispersions are of the order of 10 m/s for newly formed embryos, later lowering to less than 1 m/s by drag force cooling. After 200 orbits, 38 gravitationally bound embryos were formed inside the vortices, half of them being more massive than Mars. The embryos are composed primarily of same-sized particles. We conclude that the presence of a dead zone naturally gives rise to a population of protoplanetary cores in the mass range of 0.1-0.6 Earth masses, on very short timescales.

The alpha Centauri binary system, owing to its duplicity, proximity and brightness, and its components' likeness to the Sun, is a fundamental calibrating object for the theory of stellar structure and evolution and the determination of stellar atmospheric parameters. This role, however, is hindered by a considerable disagreement in the published analyses of its atmospheric parameters and abundances. We report a new spectroscopic analysis of both components of the alpha Centauri binary system and compare published analyses of the system. The analysis is differential with respect to the Sun, based on high-quality spectra, and employed spectroscopic and photometric methods to obtain as many independent Teff determinations as possible. The atmospheric parameters are also checked for consistency against the results of the dynamical analysis and the positions of the components in a theoretical HR diagram. We discuss possible origins of discrepancies, concluding that the presence of NLTE effects is a probable candidate, but we note that there is as yet no consensus on the existence and cause of an offset between the spectroscopic and photometric Teff scales of cool dwarfs. The spectroscopic surface gravities also agree with those derived from directly measured masses and radii. The abundance pattern can be deemed normal in the context of recent data on metal-rich stars. The position of alpha Cen A in an up-to-date theoretical evolutionary diagrams yields a good match of the evolutionary mass and age with those from the dynamical solution and seismology.

We present global 3D MHD simulations of disks of gas and solids, aiming at developing models that can be used to study various scenarios of planet formation and planet-disk interaction in turbulent accretion disks. A second goal is to show that Cartesian codes are comparable to cylindrical and spherical ones in handling the magnetohydrodynamics of the disk simulations, as the disk-in-a-box models presented here develop and sustain MHD turbulence. We investigate the dependence of the magnetorotational instability on disk scale height, finding evidence that the turbulence generated by the magnetorotational instability grows with thermal pressure. The turbulent stresses depend on the thermal pressure obeying a power law of 0.24+/-0.03, compatible with the value of 0.25 found in shearing box calculations. The ratio of stresses decreased with increasing temperature. We also study the dynamics of boulders in the hydromagnetic turbulence. The vertical turbulent diffusion of the embedded boulders is comparable to the turbulent viscosity of the flow. Significant overdensities arise in the solid component as boulders concentrate in high pressure regions.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

Context:Ages derived from the low mass stars still contracting onto the main sequence often differ from ages derived from the high mass ones that have already evolved away from it. Aims:We investigate the general claim of disagreement between these two independent age determinations by presenting UBVRI photometry of the young galactic open clusters NGC 2232, NGC 2516, NGC 2547 and NGC 4755, spanning the age range ~10-150 Myr Methods: We derive reddenings, distances and nuclear ages by fitting ZAMS and isochrones to color-magnitudes and color-color diagrams. To derive contraction ages, we use four different pre-main sequence models, with an empirically calibrated color-temperature relation to match the Pleiades cluster sequence. Results: When using exclusively the V vs. V-I color-magnitude diagram and empirically calibrated isochrones, there is consistency between nuclear and contraction ages for the studied clusters. Although the contraction ages seem systematically underestimated, in none of cases they deviate by more than one standard deviation from the nuclear ages.

With the discovery of the Sagittarius dwarf spheroidal (Ibata et al. 1994), a galaxy caught in the process of merging with the Milky Way, the hunt for other such accretion events has become a very active field of astrophysical research. The identification of a stellar ring-like structure in Monoceros, spanning more than 100 degrees (Newberg et al. 2002), and the detection of an overdensity of stars in the direction of the constellation of Canis Major (CMa, Martin et al. 2004), apparently associated to the ring, has led to the widespread belief that a second galaxy being cannibalised by the Milky Way had been found. In this scenario, the overdensity would be the remaining core of the disrupted galaxy and the ring would be the tidal debris left behind. However, unlike the Sagittarius dwarf, which is well below the Galactic plane and whose orbit, and thus tidal tail, is nearly perpendicular to the plane of the Milky Way, the putative CMa galaxy and ring are nearly co-planar with the Galactic disk. This severely complicates the interpretation of observations. In this letter, we show that our new description of the Milky Way leads to a completely different picture. We argue that the Norma-Cygnus spiral arm defines a distant stellar ring crossing Monoceros and the overdensity is simply a projection effect of looking along the nearby local arm. Our perspective sheds new light on a very poorly known region, the third Galactic quadrant (3GQ), where CMa is located.

The fact that most extrasolar planets found to date are orbiting metal-rich stars lends credence to the core accretion mechanism of gas giant planet formation over its competitor, the disc instability mechanism. However, the core accretion mechanism is not refined to the point of explaining orbital parameters such as their unexpected semi-major axes and eccentricities. We propose a model, which correlates the metallicity of the host star with the original semi-major axis of its most massive planet, prior to migration, considering that the core accretion scenario governs giant gas planet formation. The model predicts that the optimum regions for planetary formation shift inward as stellar metallicity decreases, providing an explanation for the observed absence of long period planets in metal-poor stars. We compare our predictions with the available data on extrasolar planets for stars with masses similar to the mass of the Sun. A fitting procedure produces an estimate of what we define as the Zero Age Planetary Orbit (ZAPO) curve as a function of the metallicity of the star. The model also hints that the lack of planets circling metal-poor stars may be partly caused by an enhanced destruction probability during the migration process, since the planets lie initially closer to the central stars.