We present Gemini Planet Imager polarized intensity imagery of HD 100453 in Y, J, and K1 bands that reveals an inner gap (9-18 au), an outer disk (18-39 au) with two prominent spiral arms, and two azimuthally localized dark features that are also present in Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) total intensity images. Spectral energy distribution fitting further suggests that the radial gap extends to 1 au. The narrow, wedge-like shape of the dark features appears similar to predictions of shadows cast by an inner disk that is misaligned with respect to the outer disk. Using the Monte Carlo radiative transfer code HOCHUNCK3D, we construct a model of the disk that allows us to determine its physical properties in more detail. From the angular separation of the features, we measure the difference in inclination between the disks (45 degrees) and their major axes, PA = 140 degrees east of north for the outer disk, and 100 degrees for the inner disk. We find an outer-disk inclination of 25 degrees +/- 10 degrees from face-on, in broad agreement with the Wagner et al. measurement of 34 degrees. SPHERE data in J and H bands indicate a reddish disk, which indicates that HD 100453 is evolving into a young debris disk.

The evolution of protoplanetary disks is believed to be driven largely by angular momentum transport resulting from magnetized disk winds and turbulent viscosity. The ionization of the disk that is essential for these processes has been thought to be due to host star coronal X-rays but could also arise from energetic particles produced by coronal flares, or traveling shock waves, and advected by the stellar wind. We have performed test-particle numerical simulations of energetic protons propagating into a realistic T. Tauri stellar wind, including a superposed small-scale magnetostatic turbulence. The isotropic (Kolmogorov power spectrum) turbulent component is synthesized along the individual particle trajectories. We have investigated the energy range [0.1-10] GeV, consistent with expectations from Chandra X-ray observations of large flares on T. Tauri stars and recent indications by the Herschel Space Observatory of a significant contribution of energetic particles to the disk ionization of young stars. In contrast with a previous theoretical study finding a dominance of energetic particles over X-rays in the ionization throughout the disk, we find that the disk ionization is likely dominated by X-rays over much of its area, except within narrow regions where particles are channeled onto the disk by the strongly tangled and turbulent magnetic field. The radial thickness of such regions is 5 stellar radii close to the star and broadens with increasing radial distance. This likely continues out to large distances from the star (10 au or greater), where particles can be copiously advected and diffused by the turbulent wind.

Aims. The aim of this work is to constrain properties of the disk around the archetype FU Orionis object, FU Ori, with as good as similar to 25 au resolution. Methods. We resolved FU Ori at 29-37 GHz using the Karl G. Jansky Very Large Array (JVLA) in the A-array configuration, which provided the highest possible angular resolution to date at this frequency band (similar to 0 ''.07). We also performed complementary JVLA 8-10 GHz observations, Submillimeter Array (SMA) 224 GHz and 272 GHz observations, and compared these with archival Atacama Large Millimeter Array (ALMA) 346 GHz observations to obtain the spectral energy distributions (SEDs). Results. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at frequencies higher than 29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au, which is two orders of magnitude smaller in linear scale than the gaseous disk revealed by the previous Subaru-HiCIAO 1.6 mu m coronagraphic polarization imaging observations. We are struck by the fact that these two spatially compact sources contribute to over 50% of the observed fluxes at 224 GHz, 272 GHz, and 346 GHz. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency), although we cannot rule out that it is due to the time variability of their (sub)millimeter fluxes. Conclusions. The more sophisticated models for SEDs considering the details of the observed spectral indices in the millimeter bands suggest that the >29 GHz emission is contributed by a combination of free-free emission from ionized gas and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks (less than or similar to 0.1 au) has been sublimated, and thus the disks are no longer well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few-au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal instability and the magnetorotational instability when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate.

We have imaged GM Aurigae with the Hubble Space Telescope, detected its disk in scattered light at 1400 and 1650 angstrom, and compared these with observations at 3300 angstrom, 5550 angstrom, 1.1 mu m, and 1.6 mu m. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300 angstrom shows no evidence of the 24 au radius cavity that has been previously observed in submillimeter observations. Comparison with dust grain opacity models indicates that. the surface of the entire disk is populated with submicron grains. We have compiled a. spectral energy distribution from 0.1 mu m to 1 mm. and used it to constrain a model of the star + disk system that includes the submillimeter cavity using the Monte Carlo radiative transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of submicron grains interior to the submillimeter cavity wall. We suggest one explanation for this that. could be due to a planet of mass <9 M-J interior to 24 au. A unique cylindrical structure is detected in the far-UV data from the Advanced Camera for Surveys/ Solar Blind Channel. It is aligned along the system semiminor axis, but does not resemble an accretion-driven jet. The structure is limb. brightened and extends 190 +/- 35 au above the disk midplane. The inner radius of the limb. brightening is 40 +/- 10 au, just beyond the submillimeter cavity wall.

The notion that low- to intermediate-mass young stellar objects (YSOs) gain mass at a constant rate during the early stages of their evolution appears to be challenged by observations of YSOs suffering sudden increases of the rate at which they gain mass from their circumstellar discs. Also, this idea that stars spend most of their lifetime with a low accretion rate and gain most of their final mass during short-lived episodes of high accretion bursts, helps to solve some long-standing problems in stellar evolution. The original classification of eruptive variables divides them in two separate subclasses known as FU Orionis stars (FUors) and EX Lupi stars (EXors). In this classical view FUors are at an early evolutionary stage and are still gaining mass from their parent envelopes, whilst EXors are thought to be older objects only surrounded by an accretion disc. The problem with this classical view is that it excludes younger protostars which have higher accretion rates but are too deeply embedded in circumstellar matter to be observed at optical wavelengths. Optically invisible protostars have been observed to display large variability in the near-infrared. These and some recent discoveries of new eruptive variables, show characteristics that can be attributed to both of the optically-defined subclasses of eruptive variables. The new objects have been proposed to be part of a new class of eruptive variables. However, a more accepted scenario is that in fact the original classes only represent two extremes of the same phenomena. In this sense eruptive variability could be explained as arising from one physical mechanism, i.e. unsteady accretion, where a variation in the parameters of such mechanism can cause the different characteristics observed in the members of this class. With the aim of studying the incidence of episodic accretion among young stellar objects, and to characterize the nature of these eruptive variables we searched for high amplitude variability in two multi-epoch infrared surveys: the UKIDSS Galactic Plane Survey (GPS) and the Vista Variables in the Via Lactea (VVV). In order to further investigate the nature of the selected variable stars, we use photometric information arising from public surveys at near- to farinfrared wavelengths. In addition we have performed spectroscopic and photometric follow-up for a large subset of the samples arising from GPS and VVV. We analyse the widely separated two-epoch K-band photometry in the 5th, 7th and 8th data releases of the UKIDSS Galactic Plane Survey. We find 71 stars with _K > 1 mag, including 2 previously known OH/IR stars and a Nova. Even though the mid-plane is mostly excluded from the dataset, we find the majority (66%) of our sample to be within known star forming regions (SFRs), with two large concentrations in the Serpens OB2 association (11 stars) and the Cygnus-X complex (27 stars). The analysis of the multi-epoch K-band photometry of 2010-2012 data from VVV covering the Galactic disc at |b| < 1◦ yields 816 high amplitude variables, which include known variables of different classes such as high mass X-ray binaries, Novae and eclipsing binaries among others. Remarkably, 65% of the sample are found concentrated towards areas of star formation, similar to the results from GPS. In both surveys, sources in SFRs show spectral energy distributions (SEDs) that support classification as YSOs. This indicates that YSOs dominate the Galactic population of high amplitude infrared variable stars at low luminosities and therefore likely dominate the total high amplitude population. Spectroscopic follow-up allows us to confirm the pre-main sequence nature of several GPS and VVV Objects. Most objects in both samples show spectroscopic signatures that can be attributed to YSOs undergoing high states of accretion, such as veiling of photospheric features and CO emission, or show FUor-like spectra. We also find a large fraction of objects with 2.12 μm H2 emission that can be explained as arising from shock-excited emission caused by molecular outflows. Whether these molecular outflows are related to outbursts events cannot be confirmed from our data. Adding the GPS and VVV spectroscopic results, we find that between 6 and 14 objects are new additions to the FUor class from their close resemblance to the near-infrared spectra of FUors, and at least 23 more objects are new additions to the eruptive variable class. For most of these we are unable to classify them into any of the original definitions for this variable class. In any case, we are adding up to 37 new stars to the eruptive variable class which would double the current number of known objects. We note that most objects are found to be deeply embedded optically invisible stars, thus increasing the number of objects belonging to this subclass by a much larger factor. In general, objects in our samples which are found to be likely eruptive variable stars show a mixture of characteristics that can be attributed to both of the optically-defined classes. This agrees well with the recent discoveries in the literature. Finally, we are able to derive a first rough estimate on the incidence of episodic accretion among class I YSOs in the star-forming complex G305. We find that _ 9% of such objects are in a state of high accretion. This number is in agreement with previous theoretical and observational estimates among class I YSOs.

We report ALMA Cycle 2 observations of 230 GHz (1.3 mm) dust continuum emission, and (CO)-C-12, (CO)-C-13, and (CO)-O-18 J = 2-1 line emission, from the Upper Scorpius transitional disk [PZ99] J160421.7-213028, with an angular resolution of similar to 0''.25 (35 au). Armed with these data and existing H-band scattered light observations, we measure the size and depth of the disk's central cavity, and the sharpness of its outer edge, in three components: sub-mu m-sized "small" dust traced by scattered light, millimeter-sized "big" dust traced by the millimeter continuum, and gas traced by line emission. Both dust populations feature a cavity of radius similar to 70 au that is depleted by factors of at least 1000 relative to the dust density just outside. The millimeter continuum data are well explained by a cavity with a sharp edge. Scattered light observations can be fitted with a cavity in small dust that has either a sharp edge at 60 au, or an edge that transitions smoothly over an annular width of 10 au near 60 au. In gas, the data are consistent with a cavity that is smaller, about 15 au in radius, and whose surface density at 15 au is 10(3 +/- 1) times smaller than the surface density at 70 au; the gas density grades smoothly between these two radii. The CO isotopologue observations rule out a sharp drop in gas surface density at 30 au or a double-drop model, as found by previous modeling. Future observations are needed to assess the nature of these gas and dust cavities (e.g., whether they are opened by multiple as-yet-unseen planets or photoevaporation).

High-contrast imaging instruments such as GPI and SPHERE are discovering gap structures in protoplanetary disks at an ever faster pace. Some of these gaps may be opened by planets forming in the disks. In order to constrain planet formation models using disk observations, it is crucial to find a robust way to quantitatively back out the properties of the gap-opening planets, in particular their masses, from the observed gap properties, such as their depths and widths. Combining 2D and 3D hydrodynamics simulations with 3D radiative transfer simulations, we investigate the morphology of planet-opened gaps in near-infrared scattered-light images. Quantitatively, we obtain correlations that directly link intrinsic gap depths and widths in the gas surface density to observed depths and widths in images of disks at modest inclinations under finite angular resolution. Subsequently, the properties of the surface density gaps enable us to derive the disk scale height at the location of the gap h, and to constrain the quantity M-p(2)/alpha, where Mp is the mass of the gap-opening planet and a characterizes the viscosity in the gap. As examples, we examine the gaps recently imaged by VLT/SPHERE, Gemini/GPI, and Subaru/HiCIAO in HD 97048, TW Hya, HD 169142, LkCa. 15, and RX J1615.3-3255. Scale heights of the disks and possible masses of the gap-opening planets are derived assuming each gap is opened by a single planet. Assuming a = 10(-3), the derived planet masses in all cases are roughly between 0.1 and 1M(J).

The dependence of the mass accretion rate on the stellar properties is a key constraint for star formation and disk evolution studies. Here we present a study of a sample of stars in the Chamaeleon I star-forming region carried out using spectra taken with the ESO VLT/X-shooter spectrograph. The sample is nearly complete down to stellar masses (M-star) similar to 0.1 M-circle dot for the young stars still harboring a disk in this region. We derive the stellar and accretion parameters using a self-consistent method to fit the broadband flux-calibrated medium resolution spectrum. The correlation between accretion luminosity to stellar luminosity, and of mass accretion rate to stellar mass in the logarithmic plane yields slopes of 1.9 +/- 0.1 and 2.3 +/- 0.3, respectively. These slopes and the accretion rates are consistent with previous results in various star-forming regions and with different theoretical frameworks. However, we find that a broken power-law fit, with a steeper slope for stellar luminosity lower than similar to 0.45 L-circle dot and for stellar masses lower than similar to 0.3 M-circle dot is slightly preferred according to different statistical tests, but the single power-law model is not excluded. The steeper relation for lower mass stars can be interpreted as a faster evolution in the past for accretion in disks around these objects, or as different accretion regimes in different stellar mass ranges. Finally, we find two regions on the mass accretion versus stellar mass plane that are empty of objects: one region at high mass accretion rates and low stellar masses, which is related to the steeper dependence of the two parameters we derived. The second region is located just above the observational limits imposed by chromospheric emission, at M-star similar to 0.3-0.4 M-circle dot. These are typical masses where photoevaporation is known to be effective. The mass accretion rates of this region are similar to 10(-10) M-circle dot/yr, which is compatible with the value expected for photoevaporation to rapidly dissipate the inner disk.

HD 141569 A is a pre-main sequence B9.5 Ve star surrounded by a prominent and complex circumstellar disk, likely still in a transition stage from protoplanetary to debris disk phase. Here, we present a new image of the third inner disk component of HD 141569 A made in the L' band (3.8 mu m) during the commissioning of the vector vortex coronagraph that has recently been installed in the near-infrared imager and spectrograph NIRC2 behind the W.M. Keck Observatory Keck II adaptive optics system. We used reference point-spread function subtraction, which reveals the innermost disk component from the inner working distance of similar or equal to 23 au and up to similar or equal to 70 au. The spatial scale of our detection roughly corresponds to the optical and near-infrared scattered light, thermal Q, N, and 8.6 mu m PAH emission reported earlier. We also see an outward progression in dust location from the L' band to the H band (Very Large Telescope/SPHERE image) to the visible (Hubble Space Telescope (HST)/STIS image), which is likely indicative of dust blowout. The warm disk component is nested deep inside the two outer belts imaged by HST-NICMOS in 1999 (at 406 and 245 au, respectively). We fit our new L'-band image and spectral energy distribution of HD 141569 A with the radiative transfer code MCFOST. Our best-fit models favor pure olivine grains and are consistent with the composition of the outer belts. While our image shows a putative very faint point-like clump or source embedded in the inner disk, we did not detect any true companion within the gap between the inner disk and the first outer ring, at a sensitivity of a few Jupiter masses.