We study the relationship between the mid-infrared and sub-mm variability of deeply embedded protostars using the multi-epoch data from the Wide Infrared Survey Explorer ($WISE$/NEOWISE) and the ongoing James Clerk Maxwell Telescope (JCMT) transient survey. Our search for signs of stochastic (random) and/or secular (roughly monotonic in time) variability in a sample of 59 young stellar objects (YSOs) revealed that 35 are variable in at least one of the two surveys. This variability is dominated by secular changes. Of those objects with secular variability, 14 objects ($22\%$ of the sample) show correlated secular variability over mid-IR and sub-mm wavelengths. Variable accretion is the likely mechanism responsible for this type of variability. Fluxes of YSOs that vary in both wavelengths follow a relation of $\log_{10} F_{4.6}(t)=\eta \log_{10} F_{850}(t)$ between the mid-IR and sub-mm, with $\eta=5.53\pm0.29$. This relationship arises from the fact that sub-mm fluxes respond to the dust temperature in the larger envelope whereas the mid-IR emissivity is more directly proportional to the accretion luminosity. The exact scaling relation, however, depends on the structure of the envelope, the importance of viscous heating in the disc, and dust opacity laws.

In the episodic accretion scenario, a large fraction of the protostellar mass accretes during repeated and large bursts of accretion. Since outbursts on protostars are typically identified at specific wavelengths, interpreting these outbursts requires converting this change in flux to a change in total luminosity. The Class I young stellar object EC 53 in the Serpens Main cloud has undergone repeated increases in brightness at 850 $\mu$m that are likely caused by bursts of accretion. In this study, we perform two- and three-dimensional continuum radiative transfer modeling to quantify the internal luminosity rise in EC 53 that corresponds to the factor of $\sim$1.5 enhancement in flux at 850 $\mu$m. We model the spectral energy distribution and radial intensity profile in both the quiescent and outburst phases. The internal luminosity in the outburst phase is $\sim 3.3$ times brighter than the luminosity in the quiescent phase. The radial intensity profile analysis demonstrates that the detected sub-mm flux variation of EC 53 comes from the heated envelope by the accretion burst. We also find that the role of external heating of the EC 53 envelope by the interstellar radiation field is insignificant.

Young stars exhibit variability due to changes in the gas accretion rate onto them, an effect that should be quite significant in the early stages of their formation. As protostars are embedded within their natal cloud, this variability may only be inferred through long wavelength observations. We perform radiative transfer simulations of young stellar objects (YSOs) formed in hydrodynamical simulations, varying the structure and luminosity properties in order to estimate the long-wavelength, sub-mm and mm, variations of their flux. We find that the flux increase due to an outburst event depends on the protostellar structure and is more prominent at sub-mm wavelengths than at mm wavelengths; e.g. a factor of 40 increase in the luminosity of the young protostar leads to a flux increase of a factor of 10 at 250 micron but only a factor of 2.5 at 1.3 mm. We find that the interstellar radiation field dilutes the flux increase but that this effect may be avoided if resolution permits the monitoring of the inner regions of a YSO, where the heating is primarily due to protostellar radiation. We also confirm that the bolometric temperature and luminosity of outbursting protostars may result in an incorrect classification of their evolutionary stage.

Accretion onto protostars may occur in sharp bursts. Accretion bursts during the embedded phase of young protostars are probably most intense, but can only be inferred indirectly through long-wavelength observations. We perform radiative transfer calculations for young stellar objects (YSOs) formed in hydrodynamic simulations to predict the long wavelength, sub-mm and mm, flux responses to episodic accretion events, taking into account heating from the young protostar and from the interstellar radiation field. We find that the flux increase due to episodic accretion events is more prominent at sub-mm wavelengths than at mm wavelengths; e.g. a factor of ~570 increase in the luminosity of the young protostar leads to a flux increase of a factor of 47 at 250 micron but only a factor of 10 at 1.3 mm. Heating from the interstellar radiation field may reduce further the flux increase observed at longer wavelengths. We find that during FU Ori-type outbursts the bolometric temperature and luminosity may incorrectly classify a source as a more evolved YSO, due to a larger fraction of the radiation of the object being emitted at shorter wavelengths

Complex organic molecules (COMs), which are the seeds of prebiotic material and precursors of amino acids and sugars, form in the icy mantles of circumstellar dust grains but cannot be detected remotely unless they are heated and released to the gas phase. Around solar-mass stars, water and COMs only sublimate in the inner few au of the disk, making them extremely difficult to spatially resolve and study. Sudden increases in the luminosity of the central star will quickly expand the sublimation front (so-called snow line) to larger radii, as seen previously in the FU Ori outburst of the young star V883 Ori. In this paper, we take advantage of the rapid increase in disk temperature of V883 Ori to detect and analyze five different COMs, methanol, acetone, acetonitrile, acetaldehyde, and methyl formate, in spatially-resolved submillimeter observations. The COMs abundances in V883 Ori is in reasonable agreement with cometary values. This result suggests that outbursting young stars can provide a special opportunity to study the ice composition of material directly related to planet formation.

We present our first results of the survey for high redshift quasars at $5 \lesssim {\rm z} \lesssim 5.7$. The search for quasars in this redshift range has been known to be challenging due to limitations of filter sets used in previous studies. We conducted a quasar survey for two specific redshift ranges, 4.60 $\leq$ z $\leq$ 5.40 and 5.50 $\leq$ z $\leq$ 6.05, using multi-wavelength data that include observations using custom-designed filters, $is$ and $iz$. Using these filters and a new selection technique, we were able to reduce the fraction of interlopers. Through optical spectroscopy, we confirmed seven quasars at 4.7 $\leq$ z $\leq$ 5.4 with $-27.4 < M_{1450} < -26.4$ which were discovered independently by another group recently. We estimated black hole masses and Eddington ratios of four of these quasars from optical and near-infrared spectra, and found that these quasars are undergoing nearly Eddington-limited accretion which is consistent with the rapid growth of supermassive black holes in luminous quasars at z $\sim$ 5.

We carried out photometric observations for HBC 722 in SDSS r, i and z bands from 2011 April to 2013 May with a Camera for Quasars in Early uNiverse attached to the 2.1m Otto Struve telescope at McDonald Observatory. The post-outburst phenomena were classified into five phases according to not only brightness but also color variations, which might be caused by physical changes in the emitting regions of optical and near-infrared bands. A series of spectral energy distribution (SED) is presented to support color variations and track the time evolution of SED in optical/near-infrared bands after the outburst. Given two years of data, possible periodicities of r, i and z bands were checked. We found out three families of signals around ~6, ~10 and ~1 days in three bands, which is broadly consistent with Green et al. (2013). We also examined short term variability (intra-day and day scales) to search for evidences of flickering by using the micro-variability method. We found clear signs of day scale variability and weak indications of intra-day scale fluctuations, which implies that the flickering event occurs in HBC 722 after outburst.

We present the detection of day-timescale periodic variability in the r-band lightcurve of newly outbursting FU Orionis-type object HBC 722, taken from > 42 nights of observation with the CQUEAN instrument on the McDonald Observatory 2.1m telescope. The optical/near-IR lightcurve of HBC 722 shows a complex array of periodic variability, clustering around 5.8 day (0.044 mag amplitude) and 1.28 day (0.016 mag amplitude) periods, after removal of overall baseline variation. We attribute the unusual number of comparable strength signals to a phenomenon related to the temporary increase in accretion rate associated with FUors. We consider semi-random "flickering", magnetic braking/field compression and rotational asymmetries in the disk instability region as potential sources of variability. Assuming the 5.8 day period is due to stellar rotation and the 1.28 day period is indicative of Keplerian rotation at the inner radius of the accretion disk (at 2 R(star)), we derive a B-field strength of 2.2-2.7 kG, slightly larger than typical T Tauri stars. If instead the 5.8 day signal is from a disk asymmetry, the instability region has an outer radius of 5.4 R(star), consistent with models of FUor disks. Further exploration of the time domain in this complicated source and related objects will be key to understanding accretion processes.