Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars (Class 0/I). The Herschel-PACS spectral maps covering $\sim55-210$ $\mu$m with a field-of-view of $\sim$50'' are used to quantify the gas excitation conditions and spatial extent using rotational transitions of H$_{2}$O, high-$J$ CO, and OH, as well as [O I] and [C II]. We confirm that a warm ($\sim$300 K) CO reservoir is ubiquitous and that a hotter component ($760\pm170$ K) is frequently detected around protostars. The line emission is extended beyond $\sim$1000 AU spatial scales in 40/90 objects, typically in molecular tracers in Class 0 and atomic tracers in Class I objects. High-velocity emission ($\gtrsim90$ km s$^{-1}$) is detected in only 10 sources in the [O I] line, suggesting that the bulk of [O I] arises from gas that is moving slower than typical jets. Line flux ratios show an excellent agreement with models of $C$-shocks illuminated by UV photons for pre-shock densities of $\sim$$10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-IR molecular and atomic lines are a unique diagnostic of feedback from UV emission and shocks in envelopes of deeply embedded protostars.

We analyze results from the first eighteen months of monthly sub-mm monitoring of eight star-forming regions in the JCMT Transient Survey. In our search for stochastic variability in 1643 bright peaks, only the previously identified source, EC53, shows behavior well above the expected measurement uncertainty. Another four sources, two disks and two protostars, show moderately-enhanced standard deviations in brightness, as expected for stochastic variables. For the two protostars, this apparent variability is the result of single epochs that are much brighter than the mean. In our search for secular brightness variations that are linear in time, we measure the fractional brightness change per year for 150 bright peaks, fifty of which are protostellar. The ensemble distribution of slopes is well fit by a normal distribution with sigma ~ 0.023. Most sources are not rapidly brightening or fading in the sub-mm. Comparison against time-randomized realizations shows that the width of the distribution is dominated by the uncertainty in the individual brightness measurements of the sources. A toy model for secular variability reveals that an underlying Gaussian distribution of linear fractional brightness change sigma = 0.005 would be unobservable in the present sample, whereas an underlying distribution with sigma = 0.02 is ruled out. Five protostellar sources, 10% of the protostellar sample, are found to have robust secular measures deviating from a constant flux. The sensitivity to secular brightness variations will improve significantly with a larger time sample, with a factor of two improvement expected by the conclusion of our 36-month survey.

Hydrogen-poor superluminous supernovae (SLSNe-I) have been predominantly found in low-metallicity, star-forming dwarf galaxies. Here we identify Gaia17biu/SN 2017egm as an SLSN-I occurring in a "normal" spiral galaxy (NGC 3191) in terms of stellar mass (several times 10^10 M_sun) and metallicity (roughly Solar). At redshift z=0.031, Gaia17biu is also the lowest redshift SLSN-I to date, and the absence of a larger population of SLSNe-I in dwarf galaxies of similar redshift suggests that metallicity is likely less important to the production of SLSNe-I than previously believed. With the smallest distance and highest apparent brightness for an SLSN-I, we are able to study Gaia17biu in unprecedented detail. Its pre-peak near-ultraviolet to optical color is similar to that of Gaia16apd and among the bluest observed for an SLSN-I while its peak luminosity (M_g = -21 mag) is substantially lower than Gaia16apd. Thanks to the high signal-to-noise ratios of our spectra, we identify several new spectroscopic features that may help to probe the properties of these enigmatic explosions. We detect polarization at the ~0.5% level that is not strongly dependent on wavelength, suggesting a modest, global departure from spherical symmetry. In addition, we put the tightest upper limit yet on the radio luminosity of an SLSN-I with <5.4x10^26 erg/s/Hz (at 10 GHz), which is almost a factor of 40 better than previous upper limits and one of the few measured at an early stage in the evolution of an SLSN-I. This limit largely rules out an association of this SLSNe-I with known populations of gamma-ray burst (GRB) like central engines.

The variability of young stellar objects is mostly driven by star-disk interactions. In long-term photometric monitoring of the accreting T Tauri star GI Tau, we detect extinction events with typical depths of $\Delta V \sim 2.5$ mag that last for days-to-months and often appear to occur stochastically. In 2014 - 2015, extinctions that repeated with a quasi-period of 21 days over several months is the first empirical evidence of slow warps predicted from MHD simulations to form at a few stellar radii away from the central star. The reddening is consistent with $R_V=3.85\pm0.5$ and, along with an absence of diffuse interstellar bands, indicates that some dust processing has occurred in the disk. The 2015 -- 2016 multi-band lightcurve includes variations in spot coverage, extinction, and accretion, each of which results in different traces in color-magnitude diagrams. This lightcurve is initially dominated by a month-long extinction event and return to the unocculted brightness. The subsequent light-curve then features spot modulation with a 7.03 day period, punctuated by brief, randomly-spaced extinction events. The accretion rate measured from $U$-band photometry ranges from $1.3\times10^{-8}$ to $1.1\times10^{-10}$ M$_\odot$ yr$^{-1}$ (excluding the highest and lowest 5% of high- and low- accretion rate outliers), with an average of $4.7 \times 10^{-9}$ M$_\odot$ yr$^{-1}$. A total of 50% of the mass is accreted during bursts of $>12.8\times10^{-9}$ M$_\odot$ yr${^{-1}}$, which indicates limitations on analyses of disk evolution using single-epoch accretion rates.

Investigating variability at the earliest stages of low-mass star formation is fundamental in understanding how a protostar assembles mass. While many simulations of protostellar disks predict non-steady accretion onto protostars, deeper investigation requires robust observational constraints on the frequency and amplitude of variability events characterised across the observable SED. In this study, we develop methods to robustly analyse repeated observations of an area of the sky for submillimetre variability in order to determine constraints on the magnitude and frequency of deeply embedded protostars. We compare \mbox{850 $\mu$m} JCMT Transient Survey data with archival JCMT Gould Belt Survey data to investigate variability over 2-4 year timescales. Out of 175 bright, independent emission sources identified in the overlapping fields, we find 7 variable candidates, 5 of which we classify as \textit{Strong} and the remaining 2 as \textit{Extended} to indicate the latter are associated with larger-scale structure. For the \textit{Strong} variable candidates, we find an average fractional peak brightness change per year of |4.0|\% yr$^{-1}$ with a standard deviation of $2.7\%\mathrm{\:yr}^{-1}$. In total, 7\% of the protostars associated with \mbox{850 $\mu$m} emission in our sample show signs of variability. Four of the five \textit{Strong} sources are associated with a known protostar. The remaining source is a good follow-up target for an object that is anticipated to contain an enshrouded, deeply embedded protostar. In addition, we estimate the \mbox{850 $\mu$m} periodicity of the submillimetre variable source, EC 53, to be \mbox{567 $\pm$ 32 days} based on the archival Gould Belt Survey data.

During the protostellar phase of stellar evolution, accretion onto the star is expected to be variable, but this suspected variability has been difficult to detect because protostars are deeply embedded. In this paper, we describe a sub-mm luminosity burst of the Class I protostar EC 53 in Serpens Main, the first variable found during our dedicated JCMT/SCUBA-2 monitoring program of eight nearby star-forming regions. EC 53 remained quiescent for the first 6 months of our survey, from February to August 2016. The sub-mm emission began to brighten in September 2016, reached a peak brightness of $1.5$ times the faint state, and has been decaying slowly since February 2017. The change in sub-mm brightness is interpreted as dust heating in the envelope, generated by a luminosity increase of the protostar of a factor of $\ge 4$. The 850~$\mu$m lightcurve resembles the historical $K$-band lightcurve, which varies by a factor of $\sim 6$ with a 543 period and is interpreted as accretion variability excited by interactions between the accretion disk and a close binary system. The predictable detections of accretion variability observed at both near-infrared and sub-mm wavelengths make the system a unique test-bed, enabling us to capture the moment of the accretion burst and to study the consequences of the outburst on the protostellar disk and envelope.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

We investigate the relation between disk mass and mass accretion rate to constrain the mechanism of angular momentum transport in protoplanetary disks. Dust mass and mass accretion rate in Chamaeleon I are correlated with a slope close to linear, similar to the one recently identified in Lupus. We investigate the effect of stellar mass and find that the intrinsic scatter around the best-fit Mdust-Mstar and Macc-Mstar relations is uncorrelated. Disks with a constant alpha viscosity can fit the observed relations between dust mass, mass accretion rate, and stellar mass, but over-predict the strength of the correlation between disk mass and mass accretion rate when using standard initial conditions. We find two possible solutions. 1) The observed scatter in Mdust and Macc is not primoridal, but arises from additional physical processes or uncertainties in estimating the disk gas mass. Most likely grain growth and radial drift affect the observable dust mass, while variability on large time scales affects the mass accretion rates. 2) The observed scatter is primordial, but disks have not evolved substantially at the age of Lupus and Chamaeleon I due to a low viscosity or a large initial disk radius. More accurate estimates of the disk mass and gas disk sizes in a large sample of protoplanetary disks, either through direct observations of the gas or spatially resolved multi-wavelength observations of the dust with ALMA, are needed to discriminate between both scenarios or to constrain alternative angular momentum transport mechanisms such as MHD disk winds.

TWA 3A is the most recent addition to a small group of young binary systems that both actively accrete from a circumbinary disk and have spectroscopic orbital solutions. As such, it provides a unique opportunity to test binary accretion theory in a well-constrained setting. To examine TWA 3A's time-variable accretion behavior, we have conducted a two-year, optical photometric monitoring campaign, obtaining dense orbital phase coverage (~20 observations per orbit) for ~15 orbital periods. From U-band measurements we derive the time-dependent binary mass accretion rate, finding bursts of accretion near each periastron passage. On average, these enhanced accretion events evolve over orbital phases 0.85 to 1.05, reaching their peak at periastron. The specific accretion rate increases above the quiescent value by a factor of ~4 on average but the peak can be as high as an order of magnitude in a given orbit. The phase dependence and amplitude of TWA 3A accretion is in good agreement with numerical simulations of binary accretion with similar orbital parameters. In these simulations, periastron accretion bursts are fueled by periodic streams of material from the circumbinary disk that are driven by the binary orbit. We find that TWA 3A's average accretion behavior is remarkably similar to DQ Tau, another T Tauri binary with similar orbital parameters, but with significantly less variability from orbit to orbit. This is only the second clear case of orbital-phase-dependent accretion in a T Tauri binary.

Past estimates for the age of the Upper Sco Association are typically 11-13 Myr for intermediate-mass stars and 4-5 Myr for low-mass stars. In this study, we simulate populations of young stars to investigate whether this apparent dependence of estimated age on spectral type may be explained by the star formation history of the association. Solar and intermediate mass stars begin their pre-main sequence evolution on the Hayashi track, with fully convective interiors and cool photospheres. Intermediate mass stars quickly heat up and transition onto the radiative Henyey track. As a consequence, for clusters in which star formation occurs on a similar timescale as the transition from a convective to a radiative interior, discrepancies in ages will arise when ages are calculated as a function of temperature instead of mass. Simple simulations of a cluster with constant star formation over several Myr may explain about half of the difference in inferred ages versus photospheric temperature; speculative constructions that consist of a constant star formation followed by a large supernova-driven burst could fully explain the differences, including those between F and G stars where evolutionary tracks may be more accurate. The age spreads of low-mass stars predicted from these prescriptions for star formation are consistent with the observed luminosity spread of Upper Sco. The conclusion that a lengthy star formation history will yield a temperature dependence in ages is expected from the basic physics of pre-main sequence evolution and is qualitatively robust to the large uncertainties in pre-main sequence evolutionary models.

The mass of a protoplanetary disk limits the formation and future growth of any planet. Masses of protoplanetary disks are usually calculated from measurements of the dust continuum emission by assuming an interstellar gas-to-dust ratio. To investigate the utility of CO as an alternate probe of disk mass, we use ALMA to survey $^{13}$CO and C$^{18}$O J = $3-2$ line emission from a sample of 93 protoplanetary disks around stars and brown dwarfs with masses from 0.03 -- 2 M$_{\odot}$ in the nearby Chamaeleon I star-forming region. We detect $^{13}$CO emission from 17 sources and C$^{18}$O from only one source. Gas masses for disks are then estimated by comparing the CO line luminosities to results from published disk models that include CO freeze-out and isotope-selective photodissociation. Under the assumption of a typical ISM CO-to-H$_2$ ratios of $10^{-4}$, the resulting gas masses are implausibly low, with an average gas mass of $\sim$ 0.05 M$_{Jup}$ as inferred from the average flux of stacked $^{13}$CO lines. The low gas masses and gas-to-dust ratios for Cha I disks are both consistent with similar results from disks in the Lupus star-forming region. The faint CO line emission may instead be explained if disks have much higher gas masses, but freeze-out of CO or complex C-bearing molecules is underestimated in disk models. The conversion of CO flux to CO gas mass also suffers from uncertainties in disk structures, which could affect gas temperatures. CO emission lines will only be a good tracer of the disk mass when models for C and CO depletion are confirmed to be accurate.

The Taurus-Auriga association and its associated molecular cloud are a benchmark population for studies of star and planet formation. The census of Taurus-Auriga has been assembled over seven decades and has inherited the biases, incompleteness, and systematic uncertainties of the input studies. The notably unusual shape of the inferred IMF and the existence of several isolated disk-bearing stars suggest that additional (likely disk-free) members might remain to be discovered. We therefore have begun a global reassessment of the membership of Taurus-Auriga that exploits new data and better definitions of youth and kinematic membership. As a first step, we reconsider the membership of all disk-free candidate members from the literature with spectral type $\ge$F0, $3^h50^m<\alpha<5^h40^m$, and $14^{\circ}<\delta<34^{\circ}$. We combine data from the literature with Keck/HIRES and UH88/SNIFS spectra to test the membership of these candidates using HR diagram positions, proper motions, RVs, H$\alpha$, lithium, and surface gravity. We find 218 confirmed or likely Taurus members, 160 confirmed or likely interlopers, and only 18 that still lack sufficient evidence to draw firm conclusions. A significant fraction of these stars (81/218=37%) are not included in the most recent canonical member lists. Intriguingly, there are few additional members in the immediate vicinity of the molecular clouds, preserving the IMFs that have been deemed anomalous in past work. Many of the likely Taurus members are distributed broadly across the search area. When combined with known disk hosts, our updated census reveals two regimes: a high-density population with a high disk fraction (indicative of youth) that broadly traces the molecular clouds, and a low-density population with low disk fraction (hence likely older) that most likely represents previous generations of star formation.

Ages and masses of young stars are often estimated by comparing their luminosities and effective temperatures to pre-main sequence stellar evolution tracks, but magnetic fields and starspots complicate both the observations and evolution. To understand their influence, we study the heavily-spotted weak-lined T-Tauri star LkCa 4 by searching for spectral signatures of radiation originating from the starspot or starspot groups. We introduce a new methodology for constraining both the starspot filling factor and the spot temperature by fitting two-temperature stellar atmosphere models constructed from Phoenix synthetic spectra to a high-resolution near-IR IGRINS spectrum. Clearly discernable spectral features arise from both a hot photospheric component $T_{\mathrm{hot}} \sim4100$ K and to a cool component $T_{\mathrm{cool}} \sim2700-3000$ K, which covers $\sim80\%$ of the visible surface. This mix of hot and cool emission is supported by analyses of the spectral energy distribution, rotational modulation of colors and of TiO band strengths, and features in low-resolution optical/near-IR spectroscopy. Although the revised effective temperature and luminosity make LkCa 4 appear much younger and lower mass than previous estimates from unspotted stellar evolution models, appropriate estimates will require the production and adoption of spotted evolutionary models. Biases from starspots likely afflict most fully convective young stars and contribute to uncertainties in ages and age spreads of open clusters. In some spectral regions starspots act as a featureless veiling continuum owing to high rotational broadening and heavy line-blanketing in cool star spectra. Some evidence is also found for an anti-correlation between the velocities of the warm and cool components.

Binary star-formation theory predicts that close binaries (a<100 AU) will experience periodic pulsed accretion events as streams of material form at the inner edge of a circumbinary disk, cross a dynamically cleared gap, and feed circumstellar disks or accrete directly onto the stars. The archetype for the pulsed-accretion theory is the eccentric, short-period, classical T Tauri binary DQ Tau. Low-cadence (~daily) broadband photometry has shown brightening events near most periastron passages, just as numerical simulations would predict for an eccentric binary. Magnetic reconnection events (flares) during the collision of stellar magnetospheres near periastron could, however, produce the same periodic, broadband behavior when observed at a one-day cadence. To reveal the dominate physical mechanism seen in DQ Tau's low-cadence observations, we have obtained continuous, moderate-cadence, multi-band photometry over 10 orbital periods, supplemented with 27 nights of minute-cadence photometry centered on 4 separate periastron passages. While both accretion and stellar flares are present, the dominant timescale and morphology of brightening events are characteristic of accretion. On average, the mass accretion rate increases by a factor of 5 near periastron, in good agreement with recent models. Large variability is observed in the morphology and amplitude of accretion events from orbit-to-orbit. We argue this is due to the absence of stable circumstellar disks around each star, compounded by inhomogeneities at the inner edge of the circumbinary disk and within the accretion streams themselves. Quasi-periodic apastron accretion events are also observed, which are not predicted by binary accretion theory.

The evolution of HII regions/supershells can trigger a new generation of stars/clusters at their peripheries, with environmental conditions that may affect the initial mass function, disk evolution and star formation efficiency. In this paper we study the stellar content and star formation processes in the young cluster Stock 8, which itself is thought to be formed during the expansion of a supershell. We present deep optical photometry along with JHK and 3.6, 4.5 {\mu}m photometry from UKIDSS and Spitzer-IRAC. We use multi-color criteria to identify the candidate young stellar objects in the region. Using evolutionary models, we obtain a median log(age) of ~6.5 (~3.0 Myr) with an observed age spread of ~0.25 dex for the cluster. Monte Carlo simulations of the population of Stock 8, based on estimates for the photometric uncertainty, differential reddening, binarity, and variability, indicate that these uncertainties introduce an age spread of ~0.15 dex. The intrinsic age spread in the cluster is ~0.2 dex. The fraction of young stellar objects surrounded by disk is ~35%. The K-band luminosity function of Stock 8 is similar to that of the Trapezium cluster. The IMF of Stock 8 has a Salpeter- like slope at >0.5 Msun and the IMF flattens and peaks at ~0.4 Msun, below which declines into the substellar regime. Although Stock 8 is surrounded by several massive stars, there seems to be no severe environmental effect in the form of IMF due to the proximity of massive stars around the cluster.

We present ALMA observations of the GQ Lup system, a young Sun-like star with a substellar mass companion in a wide-separation orbit. These observations of 870 $\mu$m continuum and CO J=3-2 line emission with beam size $\sim0.3''$ ($\sim45$ AU) resolve the disk of dust and gas surrounding the primary star, GQ Lup A, and provide deep limits on any circumplanetary disk surrounding the companion, GQ Lup b. The circumprimary dust disk is compact with a FWHM of $59\pm12$ AU, while the gas has a larger extent with a characteristic radius of $46.5\pm1.8$ AU. By forward-modeling the velocity field of the circumprimary disk based on the CO emission, we constrain the mass of GQ Lup A to be $M_* = (1.03\pm0.05)*(d/156\text{ pc})$ $M_\odot$, where $d$ is a known distance, and determine that we view the disk at an inclination angle of $60.5^\circ\pm0.5^\circ$ and a position angle of $346^\circ \pm1^\circ$. The $3\sigma$ upper limit on the 870 $\mu$m flux density of any circumplanetary disk associated with GQ Lup b of $<0.15$ mJy implies an upper limit on the dust disk mass of $<0.04$ $M_\oplus$ for standard assumptions about optically thin emission. We discuss proposed mechanisms for the formation of wide-separation substellar companions given the non-detection of circumplanetary disks around GQ Lup b and other similar systems.

Forbidden neon emission from jets of low-mass young stars can be used to probe the underlying high-energy processes in these systems. We analyze spectra of the jet of DG Tau obtained with the Very Large Telescope/X-Shooter spectrograph in 2010. [Ne III] $\lambda$3869 is clearly detected in the innermost 3" microjet and the outer knot located at $\sim$6".5. The velocity structure of the inner microjet can be decomposed into the low-velocity component (LVC) at $\sim -70$ km/s and the high-velocity component (HVC) at $\sim -180$ km/s. Based on the observed [Ne III] flux and its spatial extent, we suggest the origins of the [Ne III] emission regions and their relation with known X-ray sources along the jet. The flares from the hard X-ray source close to the star may be the main ionization source of the innermost microjet. The fainter soft X-ray source at 0".2 from the star may provide sufficient heating to help to sustain the ionization fraction against the recombination in the flow. The outer knot may be reionized by shocks faster than 100 km/s such that [Ne III] emission reappears and that the soft X-ray emission at 5".5 is produced. Velocity decomposition of the archival Hubble Space Telescope spectra obtained in 1999 shows that the HVC had been faster, with a velocity centroid of $\sim -260$ km/s. Such a decrease in velocity may potentially be explained by the expansion of the stellar magnetosphere, changing the truncation radius and thus the launching speed of the jet. The energy released by magnetic reconnections during relaxation of the transition can heat the gas up to several tens of megakelvin and provide the explanation for on-source keV X-ray flares that ionize the neon microjet.

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

One of the key questions in the field of star formation is the role of stellar feedback on subsequent star formation process. The W3 giant molecular cloud complex at the western border of the W4 super bubble is thought to be influenced by the stellar winds of the massive stars in W4. AFGL333 is a ~10^4 Msun cloud within W3. This paper presents a study of the star formation activity within AFGL333 using deep JHKs photometry obtained from the NOAO Extremely Wide-Field Infrared Imager combined with Spitzer-IRAC-MIPS photometry. Based on the infrared excess, we identify 812 candidate young stellar objects in the complex, of which 99 are classified as Class I and 713 are classified as Class II sources. The stellar density analysis of young stellar objects reveals three major stellar aggregates within AFGL333, named here AFGL333-main, AFGL333-NW1 and AFGL333-NW2. The disk fraction within AFGL333 is estimated to be ~50-60%. We use the extinction map made from the H-Ks colors of the background stars to understand the cloud structure and to estimate the cloud mass. The CO-derived extinction map corroborates the cloud structure and mass estimates from NIR color method. From the stellar mass and cloud mass associated with AFGL333, we infer that the region is currently forming stars with an efficiency of ~4.5% and at a rate of ~2 - 3 Msun Myr-1pc-2. In general, the star formation activity within AFGL333 is comparable to that of nearby low mass star-forming regions. We do not find any strong evidence to suggest that the stellar feedback from the massive stars of nearby W4 super bubble has affected the global star formation properties of the AFGL333 region.

We have constructed a series of non-rotating quasi-hydrostatic evolutionary models for the M2 Iab supergiant Betelgeuse ($\alpha~Orionis$). Our models are constrained by multiple observed values for the temperature, luminosity, surface composition and mass loss for this star, along with the parallax distance and high resolution imagery that determines its radius. We have then applied our best-fit models to analyze the observed variations in surface luminosity and the size of detected surface bright spots as the result of up-flowing convective material from regions of high temperature in the surface convective zone. We also attempt to explain the intermittently observed periodic variability in a simple radial linear adiabatic pulsation model. Based upon the best fit to all observed data, we suggest a best progenitor mass estimate of $ 20 ^{+5}_{-3} M_\odot$ and a current age from the start of the zero-age main sequence of $8.0 - 8.5$ Myr based upon the observed ejected mass while on the giant branch.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

Absolute ages of young stars are important for many issues in pre-main sequence stellar and circumstellar evolution but are long recognized as difficult to derive and calibrate. In this paper, we use literature spectral types and photometry to construct empirical isochrones in HR diagrams for low-mass stars and brown dwarfs in the eta Cha, epsilon Cha, and TW Hya Associations and the beta Pic and Tuc-Hor Moving Groups. A successful theory of pre-main sequence evolution should match the shapes of the stellar loci for these groups of young stars. However, when comparing the combined empirical isochrones to isochrones predicted from evolutionary models, discrepancies lead to a spectral type (mass) dependence in stellar age estimates. Improved prescriptions for convection and boundary conditions in the latest models of pre-main sequence models lead to a significantly improved correspondence between empirical and model isochrones, with small offsets at low temperatures that may be explained by observational uncertainties or by model limitations. Independent of model predictions, linear fits to combined stellar loci of these regions provide a simple empirical method to order clusters by luminosity with a reduced dependence on spectral type. Age estimates calculated from various sets of modern models that reproduce Li depletion boundary ages of the beta Pic Moving Group also imply a ~4 Myr age for the low mass members of the Upper Sco OB Association, which is younger than the 11 Myr age that has been recently estimated for intermediate mass members.

We present optical and near-infrared adaptive optics (AO) imaging and spectroscopy of 13 ultracool (>M6) companions to late-type stars (K7-M4.5), most of which have recently been identified as candidate members of nearby young moving groups (YMGs; 8-120 Myr) in the literature. The inferred masses of the companions (~10-100 Mjup) are highly sensitive to the ages of the primary stars so we critically examine the kinematic and spectroscopic properties of each system to distinguish bona fide YMG members from old field interlopers. 2MASS J02155892-0929121 C is a new M7 substellar companion (40-60 Mjup) with clear spectroscopic signs of low gravity and hence youth. The primary, possibly a member of the ~40 Myr Tuc-Hor moving group, is visually resolved into three components, making it a young low-mass quadruple system in a compact (<100 AU) configuration. In addition, Li 1 $\lambda$6708 absorption in the intermediate-gravity M7.5 companion 2MASS J15594729+4403595 B provides unambiguous evidence that it is young (<200 Myr) and resides below the hydrogen burning limit. Three new close-separation (<1") companions (2MASS J06475229-2523304 B, PYC J11519+0731 B, and GJ 4378 Ab) orbit stars previously reported as candidate YMG members, but instead are likely old (>1 Gyr) tidally-locked spectroscopic binaries without convincing kinematic associations with any known moving group. The high rate of false positives in the form of old active stars with YMG-like kinematics underscores the importance of radial velocity and parallax measurements to validate candidate young stars identified via proper motion and activity selection alone. Finally, we spectroscopically confirm the cool temperature and substellar nature of HD 23514 B, a recently discovered M8 benchmark brown dwarf orbiting the dustiest-known member of the Pleiades. [Abridged]

We conducted Herschel/PACS observations of five very low-mass stars or brown dwarfs located in the TW Hya association with the goal of characterizing the properties of disks in the low stellar mass regime. We detected all five targets at $70\,\mu{\rm{m}}$ and $100\,\mu{\rm{m}}$ and three targets at $160\,\mu{\rm{m}}$. Our observations, combined with previous photometry from 2MASS, WISE, and SCUBA-2, enabled us to construct SEDs with extended wavelength coverage. Using sophisticated radiative transfer models, we analyzed the observed SEDs of the five detected objects with a hybrid fitting strategy that combines the model grids and the simulated annealing algorithm and evaluated the constraints on the disk properties via the Bayesian inference method. The modelling suggests that disks around low-mass stars and brown dwarfs are generally flatter than their higher mass counterparts, but the range of disk mass extends to well below the value found in T Tauri stars, and the disk scale heights are comparable in both groups. The inferred disk properties (i.e., disk mass, flaring, and scale height) in the low stellar mass regime are consistent with previous findings from large samples of brown dwarfs and very low-mass stars. We discuss the dependence of disk properties on their host stellar parameters and find a significant correlation between the Herschel far-IR fluxes and the stellar effective temperatures, probably indicating that the scaling between the stellar and disk masses (i.e., $M_{\rm{disk}} \propto M_{\star}$) observed mainly in low-mass stars may extend down to the brown dwarf regime.

The 9 Myr old TW Hya Association (TWA) is the nearest group (typical distances of $\sim$50 pc) of pre-main-sequence (PMS) stars with ages less than 10 Myr and contains stars with both actively accreting disks and debris disks. We have studied the coronal X-ray emission from a group of low mass TWA common proper motion binaries using the {\it{Chandra}} and {\it{Swift}} satellites. Our aim is to understand better their coronal properties and how high energy photons affect the conditions around young stars and their role in photo-exciting atoms, molecules and dust grains in circumstellar disks and lower density circumstellar gas. Once planet formation is underway, this emission influences protoplanetary evolution and the atmospheric conditions of the newly-formed planets. The X-ray properties for 7 individual stars (TWA 13A, TWA 13B, TWA 9A, TWA 9B, TWA 8A, TWA 8B, and TWA 7) and 2 combined binary systems (TWA 3AB and TWA 2AB) have been measured. All the stars with sufficient signal require two-component fits to their CCD-resolution X-ray spectra, typically with a dominant hot (~2 kev (25 MK)) component and a cooler component at ~0.4 keV (4 MK). The brighter sources all show significant X-ray variability (at a level of 50-100\% of quiescence) over the course of 5-15 ksec observations due to flares. We present the X-ray properties for each of the stars and find that the coronal emission is in the super-saturated rotational domain.