In the standard picture for low-mass star formation, a dense molecular cloud undergoes gravitational collapse to form a protostellar system consisting of a new central star, a circumstellar disk, and a surrounding envelope of remaining material. The mass distribution of the system evolves as matter accretes from the large-scale envelope through the disk and onto the protostar. While this general picture is supported by simulations and indirect observational measurements, the specific timescales related to disk growth and envelope dissipation remain poorly constrained. In this paper we conduct a rigorous test of a method introduced by J{\o}rgensen et al. (2009) to obtain mass measurements of disks and envelopes around embedded protostars with observations that do not resolve the disk (resolution of $\sim$1000\,AU). Using unresolved data from the recent Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, we derive disk and envelope mass estimates for $59$ protostellar systems in the Perseus molecular cloud. We compare our results to independent disk mass measurements from the VLA Nascent Disk and Multiplicity (VANDAM) survey and find a strong linear correlation, suggesting that accurate disk masses can be measured from unresolved observations. Then, leveraging the size of the MASSES sample, we find no significant trend in protostellar mass distribution as a function of age, as approximated from bolometric temperatures. These results may indicate that the disk mass of a protostar is set near the onset of the Class 0 protostellar stage and remains roughly constant throughout the Class I protostellar stage.

(Abridged) Astronomers recently started discovering exoplanets around binary systems. Therefore, understanding the formation and evolution of circumbinary disks is crucial for a complete scenario of planet formation. The aim of this paper is to present the detection of a circumbinary disk around Oph-IRS67 and analyse its structure. We present high-angular-resolution (0.4", 60 AU) observations of C17O, H13CO+ , C34S, SO2, C2H and c-C3H2 molecular transitions with ALMA at 0.8 mm. The spectrally and spatially resolved maps reveal the kinematics of the circumbinary disk as well as its chemistry. Molecular abundances are estimated using RADEX. The continuum emission reveals the presence of a circumbinary disk around the two sources. This disk has a diameter of ~620 AU and is well traced by C17O and H13CO+ emission. C2H and c-C3H2 trace a higher-density region which is spatially offset from the sources (~430 AU). Finally, SO2 shows compact emission around one of the sources, Oph-IRS67 B. The molecular transitions which trace the circumbinary disk are consistent with a Keplerian profile on disk scales (< 200 AU) and an infalling profile for envelope scales (> 200 AU). The Keplerian fit leads to a mass of 2.2 Msun. Inferred CO abundances w.r.t. H2 are comparable to the canonical ISM value of 2.7e-4. This study proves the first detection of the circumbinary disk associated with Oph-IRS67. The disk is chemically differentiated from the nearby high-density region. The lack of methanol emission suggests the extended disk dominates the mass budget in the inner- most regions of the protostellar envelope, generating a flat density profile where less material is exposed to high temperatures. Thus, complex organic molecules would be associated with lower column densities. Finally, Oph-IRS67 is a promising candidate for the detection of both circumstellar disks with higher-angular-resolution observations.

Using zoom-simulations carried out with the adaptive mesh-refinement code RAMSES with a dynamic range of up to $2^{27} \approx 1.34 \times 10^8$ we investigate the accretion profiles around six stars embedded in different environments inside a (40 pc)$^3$ giant molecular cloud, the role of mass infall and disc instabilities on the accretion profile, and thus on the luminosity of the forming protostar. Our results show that the environment in which the protostar is embedded determines the overall accretion profile of the protostar. Infall on to the circumstellar disc may trigger gravitational disc instabilities in the disc at distances of around ~10 to ~50 au leading to rapid transport of angular momentum and strong accretion bursts. These bursts typically last for about ~10 to a ~100 yr, consistent with typical orbital times at the location of the instability, and enhance the luminosity of the protostar. Calculations with the stellar evolution code mesa show that the accretion bursts induce significant changes in the protostellar proper- ties, such as the stellar temperature and radius. We apply the obtained protostellar properties to produce synthetic observables with RADMC3D and predict that accretion bursts lead to ob- servable enhancements around 20 to 200 $\mu$m in the spectral energy distribution of Class 0 type young stellar objects.

We compute the star formation rate (SFR) in molecular clouds (MCs) that originate {\it ab initio} in a new, higher-resolution simulation of supernova-driven turbulence. Because of the large number of well-resolved clouds with self-consistent boundary and initial conditions, we obtain a large range of cloud physical parameters with realistic statistical distributions, an unprecedented sample of star-forming regions to test SFR models and to interpret observational surveys. We confirm the dependence of the SFR per free-fall time, $SFR_{\rm ff}$, on the virial parameter, $\alpha_{\rm vir}$, found in previous simulations, and compare a revised version of our turbulent fragmentation model with the numerical results. The dependences on Mach number, ${\cal M}$, gas to magnetic pressure ratio, $\beta$, and compressive to solenoidal power ratio, $\chi$ at fixed $\alpha_{\rm vir}$ are not well constrained, because of random scatter due to time and cloud-to-cloud variations in $SFR_{\rm ff}$. We find that $SFR_{\rm ff}$ in MCs can take any value in the range $0 \le SFR_{\rm ff} \lesssim 0.2$, and its probability distribution peaks at a value $SFR_{\rm ff}\approx 0.025$, consistent with observations. The values of $SFR_{\rm ff}$ and the scatter in the $SFR_{\rm ff}$--$\alpha_{\rm vir}$ relation are consistent with recent measurements in nearby MCs and in clouds near the Galactic center. Although not explicitly modeled by the theory, the scatter is consistent with the physical assumptions of our revised model and may also result in part from a lack of statistical equilibrium of the turbulence, due to the transient nature of MCs.

Context: Understanding how accretion proceeds is a key question of star formation, with important implications for both the physical and chemical evolution of young stellar objects. In particular, very little is known about the accretion variability in the earliest stages of star formation. Aims: To characterise protostellar accretion histories towards individual sources by utilising sublimation and freeze-out chemistry of CO. Methods: A sample of 24 embedded protostars are observed with the Submillimeter Array (SMA) in context of the large program "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). The size of the C$^{18}$O emitting region, where CO has sublimated into the gas-phase, is measured towards each source and compared to the expected size of the region given the current luminosity. The SMA observations also include 1.3 mm continuum data, which are used to investigate whether a link can be established between accretion bursts and massive circumstellar disks. Results: Depending on the adopted sublimation temperature of the CO ice, between 20% and 50% of the sources in the sample show extended C$^{18}$O emission indicating that the gas was warm enough in the past that CO sublimated and is currently in the process of refreezing; something which we attribute to a recent accretion burst. Given the fraction of sources with extended C$^{18}$O emission, we estimate an average interval between bursts of 20000-50000 yr, which is consistent with previous estimates. No clear link can be established between the presence of circumstellar disks and accretion bursts, however the three closest known binaries in the sample (projected separations <20 AU) all show evidence of a past accretion burst, indicating that close binary interactions may also play a role in inducing accretion variability.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

How protostars accrete mass is one of the fundamental problems of star formation. High column densities and complex kinematical structures make direct observations challenging and they only provide a snapshot. Chemical tracers provide an interesting alternative to characterise the infall histories of protostars. Previous observations of H13CO+ towards the low-mass protostar IRAS15398-3359 showed a depression in the abundance. This is a sign of destruction of HCO+ by an enhanced presence of gaseous water in an extended region, possibly related to a recent burst in the accretion. Direct observations of water vapour can determine the exact extent of the emission and confirm the hypothesis that HCO+ is indeed a good tracer of the water snow-line. IRAS15398 was observed using ALMA at 0.5" resolution. Maps of HDO(101-000) and H218O(414-321) were taken simultaneously with observations of the CS(8-7) and N2H+(5-4) lines and continuum at 0.65 and 0.75 mm. The maps were interpreted using dust radiative transfer calculations of the protostellar infalling envelope with an outflow cavity. HDO is clearly detected and extended over the scales of the H13CO+ depression, although it is displaced by ~500 AU in the direction of the outflow. H218O is tentatively detected towards the red-shifted outflow lobe, but otherwise it is absent from the mapped region, which suggests that temperatures are low. Based on the temperature structure obtained from dust radiative transfer models, we conclude that the water was most likely released from the grains in an extended hour-glass configuration during a recent accretion burst. HDO is only detected in the region closest to the protostar, at distances of up to 500 AU. These signatures can only be explained if the luminosity has recently been increased by orders of magnitudes. Additionally, the densities in the outflow cones must be sufficiently low.

Context: Understanding how protostars accrete their mass is a central question of star formation. One aspect of this is trying to understand whether the time evolution of accretion rates in deeply embedded objects is best characterised by a smooth decline from early to late stages or by intermittent bursts of high accretion. Aims: We create synthetic observations of deeply embedded protostars in a large numerical simulation of a molecular cloud, which are compared directly to real observations. The goal is to compare episodic accretion events in the simulation to observations and to test the methodology used for analysing the observations. Methods: Simple freeze-out and sublimation chemistry is added to the simulation, and synthetic C$^{18}$O line cubes are created for a large number of simulated protostars. The spatial extent of C$^{18}$O is measured for the simulated protostars and compared directly to a sample of 16 deeply embedded protostars observed with the Submillimeter Array. If CO is distributed over a larger area than predicted based on the protostellar luminosity, it may indicate that the luminosity has been higher in the past and that CO is still in the process of refreezing. Results: Approximately 1% of the protostars in the simulation show extended C$^{18}$O emission, as opposed to approximately 50% in the observations, indicating that the magnitude and frequency of episodic accretion events in the simulation is too low relative to observations. The protostellar accretion rates in the simulation are primarily modulated by infall from the larger scales of the molecular cloud, and do not include any disk physics. The discrepancy between simulation and observations is taken as support for the necessity of disks, even in deeply embedded objects, to produce episodic accretion events of sufficient frequency and amplitude.

(abridged) Context: Both observations and simulations of embedded protostars have progressed rapidly in recent years. Bringing them together is an important step in advancing our knowledge about the earliest phases of star formation. Aims: To compare synthetic continuum images and SEDs, calculated from large-scale numerical simulations, to observational studies, thereby aiding in both the interpretation of the observations and in testing the fidelity of the simulations. Methods: The radiative transfer code RADMC-3D is used to create synthetic continuum images and SEDs of protostellar systems in a large numerical simulation of a molecular cloud. More than 13000 unique radiative transfer models are produced of a variety of different protostellar systems. Results: Over the course of 0.76 Myr the simulation forms more than 500 protostars, primarily within two sub-clusters. Synthetic SEDs are used to the calculate evolutionary tracers Tbol and Lsmm/Lbol. It is shown that, while the observed distributions of the tracers are well matched by the simulation, they generally do a poor job of tracking the protostellar ages. Disks form early in the simulation, with 40 % of the Class 0 protostars being encircled by one. The flux emission from the simulated disks is found to be, on average, a factor of 6 too low relative to real observations. The distribution of protostellar luminosities spans more than three order of magnitudes, similar to the observed distribution. Cores and protostars are found to be closely associated with one another, with the distance distribution between them being in excellent agreement with observations. Conclusions: The analysis and statistical comparison of synthetic observations to real ones is established as a powerful tool in the interpretation of observational results.

We have monitored 20 Sun-like stars in the Kepler field-of-view for excess flux with the FIES spectrograph on the Nordic Optical Telescope since the launch of Kepler spacecraft in 2009. These 20 stars were selected based on their asteroseismic properties to sample the parameter space (effective temperature, surface gravity, activity level etc.) around the Sun. Though the ultimate goal is to improve stellar dynamo models, we focus the present paper on the combination of space-based and ground-based observations can be used to test the age-rotation-activity relations. In this paper we describe the considerations behind the selection of these 20 Sun-like stars and present an initial asteroseismic analysis, which includes stellar age estimates. We also describe the observations from the Nordic Optical Telescope and present mean values of measured excess fluxes. These measurements are combined with estimates of the rotation periods obtained from a simple analysis of the modulation in photometric observations from Kepler caused by starspots, and asteroseismic determinations of stellar ages, to test relations between between age, rotation and activity.

We have discovered a new rapidly oscillating Ap star among the Kepler Mission target stars, KIC 10195926. This star shows two pulsation modes with periods that are amongst the longest known for roAp stars at 17.1 min and 18.1 min, indicating that the star is near the terminal age main sequence. The principal pulsation mode is an oblique dipole mode that shows a rotationally split frequency septuplet that provides information on the geometry of the mode. The secondary mode also appears to be a dipole mode with a rotationally split triplet, but we are able to show within the improved oblique pulsator model that these two modes cannot have the same axis of pulsation. This is the first time for any pulsating star that evidence has been found for separate pulsation axes for different modes. The two modes are separated in frequency by 55 microHz, which we model as the large separation. The star is an alpha^2 CVn spotted magnetic variable that shows a complex rotational light variation with a period of Prot = 5.68459 d. For the first time for any spotted magnetic star of the upper main sequence, we find clear evidence of light variation with a period of twice the rotation period; i.e. a subharmonic frequency of $\nu_{\rm rot}/2$. We propose that this and other subharmonics are the first observed manifestation of torsional modes in an roAp star. From high resolution spectra we determine Teff = 7400 K, log g = 3.6 and v sin i = 21 km/s. We have found a magnetic pulsation model with fundamental parameters close to these values that reproduces the rotational variations of the two obliquely pulsating modes with different pulsation axes. The star shows overabundances of the rare earth elements, but these are not as extreme as most other roAp stars. The spectrum is variable with rotation, indicating surface abundance patches.

We present results from the final six months of a survey to search for pulsations in white dwarfs and hot subdwarf stars with the Kepler spacecraft. Spectroscopic observations are used to separate the objects into accurate classes, and we explore the physical parameters of the subdwarf B (sdB) stars and white dwarfs in the sample. From the Kepler photometry and our spectroscopic data, we find that the sample contains 5 new pulsators of the V1093 Her type, one AM CVn type cataclysmic variable, and a number of other binary systems. This completes the survey for compact pulsators with Kepler. No V361 Hya type of short-period pulsating sdB stars were found in this half, leaving us with a total of one single multiperiodic V361 Hya and 13 V1093 Her pulsators for the full survey. Except for the sdB pulsators, no other clearly pulsating hot subdwarfs or white dwarfs were found, although a few low-amplitude candidates still remain. The most interesting targets discovered in this survey will be observed throughout the remainder of the Kepler Mission, providing the most long-term photometric datasets ever made on such compact, evolved stars. Asteroseismic investigations of these datasets will be invaluable in revealing the interior structure of these stars, and will boost our understanding of their evolutionary history.