While accreting through a circumstellar disk, young stellar objects are observed to undergo sudden and powerful accretion events known as FUor or EXor outbursts. Although such episodic accretion is considered to be an integral part of the star formation process, the triggers and mechanisms behind them remain uncertain. We conducted global numerical hydrodynamics simulations of protoplanetary disk formation and evolution in the thin-disk limit, assuming both magnetically layered and fully magnetorotational instability (MRI)-active disk structure. In this paper, we characterize the nature of the outbursts occurring in these simulations. The instability in the dead zone of a typical layered disk results in "MRI outbursts". We explore their progression and their dependence on the layered disk parameters as well as cloud core mass. The simulations of fully MRI-active disks showed an instability analogous to the classical thermal instability. This instability manifested at two temperatures--above approximately 1400 K and 3500 K--due to the steep dependence of Rosseland opacity on the temperature. The origin of these thermally unstable regions is related to the bump in opacity resulting from molecular absorption by water vapor and may be viewed as a novel mechanism behind some of the shorter duration accretion events. Although we demonstrated local thermal instability in the disk, more investigations are needed to confirm that a large-scale global instability will ensue. We conclude that the magnetic structure of a disk, its composition, as well as the stellar mass, can significantly affect the nature of episodic accretion in young stellar objects.

It is believed that large-scale horseshoe-like brightness asymmetries found in dozens of transitional protoplanetary discs are caused by anticyclonic vortices. These vortices can play a key role in planet formation, as mm-sized dust -- the building blocks of planets -- can be accumulated inside them. Anticyclonic vortices are formed by the Rossby wave instability, which can be excited at the gap edges opened by a giant planet or at sharp viscosity transitions of accretionally inactive regions. It is known that vortices are prone to stretching and subsequent dissolution due to disc self-gravity for canonical disc masses in the isothermal approximation. To improve the hydrodynamic model of protoplanetary discs, we include the disc thermodynamics in our model. In this paper, we present our results on the evolution of the vortices formed at the outer edge of an accretionally inactive region (dead zone) assuming an ideal equation of state and taking $PdV$ work, disc cooling in the $\beta$-approximation, and disc self-gravity into account. Thermodynamics affects the offset and the mode number (referring to the number of small vortices at the early phase) of the RWI excitation, as well as the strength, shape, and lifetime of the large-scale vortex formed through merging of the initial small vortices. We found that the inclusion of gas thermodynamics results in stronger, however decreased lifetime vortices. Our results suggest that a hypothetical vortex-aided planet formation scenario favours effectively cooling discs.

Global numerical simulations of protoplanetary disk formation and evolution were conducted in thin-disk limit, where the model included magnetically layered disk structure, a self-consistent treatment for the infall from cloud core as well as the smallest possible inner computational boundary. We compared the evolution of a layered disk with a fully magnetically active disk. We also studied how the evolution depends on the parameters of the layered disk model - the MRI triggering temperature and active layer thickness - as well as the mass of the prestellar cloud core. With the canonical values of parameters a dead zone formed within the inner $\approx$ 15 au region of the magnetically layered disk. The dead zone was not a uniform structure and long-lived, axisymmetric, gaseous rings ubiquitously formed within this region due to the action of viscous torques. The rings showed a remarkable contrast in the disk environment as compared to a fully magnetically active disk and were characterized by high surface density and low effective viscosity. Multiple gaseous rings could form simultaneously in the dead zone region which were highly dynamical and showed complex, time-dependent behavior such as inward migration, vortices, gravitational instability and large-scale spiral waves. An increase in MRI triggering temperature had only marginal effects, while changes in active layer thickness as well as the initial cloud core mass had significant effects on the structure and evolution of the inner disk. Dust with large fragmentation barrier could be trapped in the rings, which may play a key role in planet formation.

Three of the seven rocky planets (e, f, and g) in TRAPPIST-1 system orbit in the habitable zone of the host star. Therefore, water can be in liquid state at their surface being essential for life. Recent studies suggest that these planets formed beyond the snow line in a water-rich region. The initial water reservoir can be lost during the planet formation due to the stellar activity of the infant low-mass star. However, a potential subsequent water delivery event, like the late heavy bombardment (LHB) in the Solar System, can replenish planetary water reservoirs. To study this water delivery process, we set up a simple model in which an additional 5 M_Earth--50 M_Earth planet is embedded in a water-rich asteroid belt beyond the snow line of TRAPPIST-1. Asteroids perturbed out from the chaotic zone of the putative planet can enter into the inner system and accreted by the known planets. Our main finding is that the larger is the orbital distance of planet, the higher is the amount of water delivered to the planet by an LHB-like event.

The particle dynamics of the recently observed ring around dwarf planet Haumea is numerically investigated. The point mass gravitational force, a second degree and order gravity field, and the solar radiation pressure as the main perturbations are considered. The quasi-stationary state of the ring varies for different micron-sized grains and depends also on the spin-orbit resonances between the rotation rate of the main body and the orbital period of the dust particles. The simulations confirm the variable radial width of the ring observed during the transit ingress and egress. Results show that the micron sized grains, initially on circular orbits, become eccentric and form an apse-aligned ring at the observed radial distance near to the 3:1 spin-orbit resonance. It is also demonstrated that this coincidence is only apparent and independent of 3:1 resonance.

Context: Quantifying the gas content inside the dust gaps of transition disks is important to establish their origin. Aims: We seek to constrain the surface density of warm gas in the disk of HD 139614, a Herbig Ae star with a transition disk exhibiting a dust gap from 2.3 to 6 AU. Methods: We have obtained ESO/VLT CRIRES high-resolution spectra of CO ro-vibrational emission. We derived constraints on the disk's structure by modeling the line-profiles, the spectroastrometric signal, and the rotational diagrams using flat Keplerian disk models. Results: We detected v=1-0 12CO, 2-1 12CO, 1-0 13CO, 1-0 C18O, and 1-0 C17O ro-vibrational lines. 12CO v=1-0 lines have an average width of 14 km/s, Tgas of 450 K and an emitting region from 1 to 15 AU. 13CO and C18O lines are on average 70 and 100 K colder, 1 and 4 km/s narrower, and are dominated by emission at R>6 AU. The 12CO v=1-0 line-profile indicates that if there is a gap in the gas it must be narrower than 2 AU. We find that a drop in the gas surface density (delta_gas) at R<5-6 AU is required to simultaneously reproduce the line-profiles and rotational diagrams of the three CO isotopologs. Delta_gas can range from 10^-2 to 10^-4 depending on the gas-to-dust ratio of the outer disk. We find that at 1<R<6 AU the gas surface density profile is flat or increases with radius. We derive a gas column density at 1<R<6 AU of NH=3x10^19 - 10^21 cm^-2. We find a 5sigma upper limit on NCO at R<1 AU of 5x10^15 cm^-2 (NH<5x10^19 cm^-2). Conclusions: The dust gap in the disk of HD 139614 has gas. The gas surface density in the disk at R<6 AU is significantly lower than the surface density expected from HD 139614's accretion rate assuming a viscous alpha-disk model. The gas density drop, the non-negative density gradient of the gas inside 6 AU, and the absence of a wide (>2 AU) gas gap suggest the presence of an embedded <2 MJ planet at around 4 AU.

The circular restricted three body problem is investigated in the context of accretion and scattering processes. In our model a large number of identical non-interacting mass-less planetesimals are considered in planar case orbiting a star-planet system. This description allows us to investigate in dynamical systems approach the gravitational scattering and possible captures of the particles by the forming planetary embryo. Although the problem serves a large variety of complex motion, the results can be easily interpreted because of the low dimensionality of the phase space. We show that initial conditions define isolated regions of the disk, where accretion or escape of the planetesimals occur, these have, in fact, a fractal structure. The fractal geometry of these "basins" implies that the dynamics is very complex. Based on the calculated escape rates and escape times, it is also demonstrated that the planetary accretion rate is exponential for short times and follows a power-law for longer integration. A new numerical calculation of the maximum mass that a planet can reach (described by the expression of the isolation mass) is also derived.

Detecting companions to Cepheids is difficult. In most cases the bright pulsator overshines the fainter secondary. Since Cepheids play a key role in determining the cosmic distance scale it is crucial to find binaries among the calibrating stars of the period-luminosity relation. Here we present an ongoing observing project of searching for faint and close companions of Galactic Cepheids using lucky imaging technique.

In this paper we briefly summarize the design concepts of the Fly's Eye Camera System, a proposed high resolution all-sky monitoring device which intends to perform high cadence time domain astronomy in multiple optical passbands while still accomplish a high \'etendue. Fundings have already been accepted by the Hungarian Academy of Sciences in order to design and build a Fly's Eye device unit. Beyond the technical details and the actual scientific goals, this paper also discusses the possibilities and yields of a network operation involving $\sim10$ sites distributed geographically in a nearly homogeneous manner. Currently, we expect to finalize the mount assembly -- that performs the sidereal tracking during the exposures -- until the end of 2012 and to have a working prototype with a reduced number of individual cameras sometimes in the spring or summer of 2013.

The star CoRoT102781750 reveals a puzzle, showing a very complex and altering variation in different `CoRoT colours'. We established without doubt that more than a single star was situated within the CoRoT mask. Using a search for periodicity as a tool, our aim is to disentangle the composite light curve and identify the type of sources behind the variability. Both flux and magnitude light curves were used. Conversion was applied after a jump- and trend-filtering algorithm. We applied different types of period-finding techniques including MuFrAn and Period04. The amplitude and phase peculiarities obtained from the independent analysis of CoRoT r, g, and b colours and ground-based follow-up photometric observations ruled out the possibility of either a background monoperiodic or a Blazhko type RR Lyrae star being in the mask. The main target, an active star, shows at least two spotted areas that reveal a $P_rot = 8.8$ hours $(f_0 = 2.735$ c d$^{-1})$ mean rotation period. The evolution of the active regions helped to derive a period change of $dP/dt = 1.6\cdot 10^{-6}$ (18 s over the run) and a differential rotation of $\alpha = \Delta\Omega/\Omega = 0.0074$. The $0\fm 015$ linear decrease and a local $0\fm 005$ increase in the dominant period's amplitude are interpreted as a decay of the old spotted region and an appearance of a new one, respectively. A star that is detected only in the CoRoT b domain shows a $f_1 = 7.172$ c d$^{-1}$ pulsation connected to a $14\fd 83$ periodicity via an equidistant triplet structure. The best explanation for our observation is a $\beta$ Cep star with a corotating dust disk.

In this paper we briefly summarize the design concepts of the Fly's Eye Camera System, a proposed high resolution all-sky monitoring device which intends to perform high cadence time domain astronomy in multiple optical passbands while still accomplish a high \'etendue. Fundings have already been accepted by the Hungarian Academy of Sciences in order to design and build a Fly's Eye device unit. Beyond the technical details and the actual scientific goals, this paper also discusses the possibilities and yields of a network operation involving ~10 sites distributed geographically in a nearly homogeneous manner. Currently, we expect to finalize the mount assembly -- that performs the sidereal tracking during the exposures -- until the end of 2012 and to have a working prototype with a reduced number of individual cameras sometimes in the spring or summer of 2013.

HD 15082 (WASP-33) is the hottest and fastest rotating star known to harbor a transiting extrasolar planet (WASP-33b). The lack of high precision radial velocity (RV) data stresses the need for precise light curve analysis and gathering further RV data. By using available photometric and RV data, we perform a blend analysis, compute more accurate system parameters, confine the planetary mass and attempt to cast light on the observed transit anomalies. We combine the original HATNet observations and various followup data to jointly analyze the signal content and extract the transit component and use our RV data to aid the global parameter determination. The blend analysis of the combination of multicolor light curves yields the first independent confirmation of the planetary nature of WASP-33b. We clearly identify three frequency components in the 15-21 1/day regime with amplitudes 7-5 mmag. These frequencies correspond to the delta Scuti-type pulsation of the host star. None of these pulsation frequencies or their low-order linear combinations are in close resonance with the orbital frequency. We show that these pulsation components explain some but not all of the observed transit anomalies. The grand-averaged transit light curve shows that there is a ~1.5 mmag brightening shortly after the planet passes the mid-transit phase. Although the duration and amplitude of this brightening varies, it is visible even through the direct inspections of the individual transit events (some 40-60% of the followup light curves show this phenomenon). We suggest that the most likely explanation of this feature is the presence of a well-populated spot belt which is highly inclined to the orbital plane. This geometry is consistent with the inference from the spectroscopic anomalies. Finally, we constrain the planetary mass to M_p=3.27+/-0.73 M_J by using our RV data collected by the TRES spectrograph.

EX Lup is the prototype of the EXor class of young eruptive stars: objects showing repetitive brightenings due to increased accretion from the circumstellar disk to the star. In this paper, we report on medium-resolution near-infrared spectroscopy of EX\,Lup taken during its extreme outburst in 2008, as well as numerical modeling with the aim of determining the physical conditions around the star. We detect emission lines from atomic hydrogen, helium, and metals, as well as first overtone bandhead emission from carbon monoxide. Our results indicate that the emission lines are originating from gas located in a dust-free region within ~ 0.2 AU of the star. The profile of the CO bandhead indicates that the CO gas has a temperature of 2500 K, and is located in the inner edge of the disk or in the outer parts of funnel flows. The atomic metals are probably co-located with the CO. Some metallic lines are fluorescently excited, suggesting direct exposure to ultraviolet photons. The Brackett series indicates emission from hot (10000 K) and optically thin gas. The hydrogen lines display a strong spectro-astrometric signal, suggesting that the hydrogen emission is probably not coming from an equatorial boundary layer; a funnel flow or disk wind origin is more likely. This picture is broadly consistent with the standard magnetospheric accretion model usually assumed for normally accreting T Tauri stars. Our results also set constraints on the eruption mechanism, supporting a model where material piles up around the corotation radius and episodically falls onto the star.

We conducted a long-term optical photometric and spectroscopic monitoring of the strongly variable, accreting young sun-like star [KOS94] HA11, associated with the dark cloud Lynds 1340, that exhibited large amplitude (5-6 magnitudes in the I_C band) brightness variations on 2-3 years timescales, flat spectral energy distribution (SED), and extremely strong (300 < EW/Angstrom < 900) H alpha emission. In this Letter we describe the basic properties of the star, derived from our observations between 1999 and 2011, and put into context the observed phenomena. The observed variations in the emission spectra, near-infrared colors, and SED suggest that [KOS94] HA11 (spectral type: K7--M0) is an eruptive young star, possibly similar in nature to V1647 Ori: its large-scale photometric variations are governed by variable accretion rate, associated with variations in the inner disk structure. The star recently has undergone strong and rapid brightness variations, thus its further observations may offer a rare opportunity for studying structural and chemical rearrangements of the inner disk, induced by variable central luminosity.

In this Letter we present observations of recent HAT-P-13b transits. The combined analysis of published and newly obtained transit epochs shows evidence for significant transit timing variations since the last publicly available ephemerides. Variation of transit timings result in a sudden switch of transit times. The detected full range of TTV spans ~0.015 days, which is significantly more than the known TTV events exhibited by hot Jupiters. If we have detected a periodic process, its period should be at least ~3 years because there are no signs of variations in the previous observations. This argument makes unlikely that the measured TTV is due to perturbations by HAT-P-13c.

We report monitoring observations of the T Tauri star EX Lupi during its outburst in 2008 in the CO fundamental band at 4.6-5.0 um. The observations were carried out at the VLT and the Subaru Telescope at six epochs from April to August 2008, covering the plateau of the outburst and the fading phase to a quiescent state. The line flux of CO emission declines with the visual brightness of the star and the continuum flux at 5 um, but composed of two subcomponents that decay with different rates. The narrow line emission (50 km s-1 in FWHM) is near the systemic velocity of EX Lupi. These emission lines appear exclusively in v=1-0. The line widths translate to a characteristic orbiting radius of 0.4 AU. The broad line component (FWZI ~ 150 km s-1) is highly excited upto v<=6. The line flux of the component decreases faster than the narrow line emission. Simple modeling of the line profiles implies that the broad-line emitting gas is orbiting around the star at 0.04-0.4 AU. The excitation state, the decay speed of the line flux, and the line profile, indicate that the broad-line emission component is physically distinct from the narrow-line emission component, and more tightly related to the outburst event.

A possible transit of HAT-P-13c has been predicted to occur on 2010 April 28. Here we report on the results of a multi-site campaign that has been organised to detect the event. CCD photometric observations have been carried out at five observatories in five countries. We reached 30% time coverage in a 5 days interval centered on the suspected transit of HAT-P-13c. Two transits of HAT-P-13b were also observed. No transit of HAT-P-13c has been detected while the campaign was on. By a numerical experiment with 10^5 model systems we conclude that HAT-P-13c is not a transiting exoplanet with a significance level from 65% to 72%, depending on the planet parameters and the prior assumptions. We present two times of transit of HAT-P-13b ocurring at BJD 2455141.5522 +- 0.0010 and BJD 2455249.4508 +- 0.0020. The TTV of HAT-P-13b is consistent with zero within 0.001 days. The refined orbital period of HAT-P-13b is 2.916293 +- 0.000010 days.

We study the perturbations of a relatively close third star on a tidally distorted eccentric eclipsing binary. We consider both the observational consequences of the variations of the orbital elements and the interactions of the stellar rotation with the orbital revolution in the presence of dissipation. We concentrate mainly on the effect of a hypothetical third companion on both the real, and the observed apsidal motion period. We investigate how the observed period derived mainly from some variants of the O-C relates to the real apsidal motion period. We carried out both analytical and numerical investigations and give the time variations of the orbital elements of the binary both in the dynamical and the observational reference frames. We give the direct analytical form of an eclipsing O-C affected simultaneously by the mutual tidal forces and the gravitational interactions with a tertiary. We also integrated numerically simultaneously the orbital and rotational equations for the possible hierarchical triple stellar system AS Camelopardalis. We find that there is a significant domain of the possible hierarchical triple system configurations, where both the dynamical and the observational effects tend to measure longer apsidal advance rate than is expected theoretically. This happens when the mutual inclination of the close and the wide orbits is large, and the orbital plane of the tertiary almost coincides with the plane of the sky. We also obtain new numerical results on the interaction of the orbital evolution and stellar rotation in such triplets. The most important fact is that resonances might occur as the stellar rotational rate varies during the dissipation-driven synchronization process...