Source author record

Robert Gutermuth

Robert Gutermuth appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2026arXiv

The Structure and Kinematics of Three Class 0 Protostellar Jets from JWST

We present observations of jets within 2000 au of three deeply embedded protostars using 2.9-27 micron observations with JWST. These observations show the morphologies and kinematics of the collimated jets from three protostars, the low-mass Class 0 protostars B335 and HOPS 153, and the intermediate-mass protostar HOPS 370. These jets are traced by shock-ionized fine-structure line emission observed with the JWST NIRSpec and MIRI IFUs. We find that [Fe II] emission traces the full extent of the inner 1000 to 2000 au of the jets, depending on distance to the protostar, while other ions mostly trace isolated shocked knots. The jets show evidence of wiggling motion in the plane of the sky as well as asymmetries between blue and red-shifted lobes. The widths of the jets increase non-monotonically with distance from the central protostar, with opening angles ranging from 2.1 degrees to < 10.1 degrees for the three protostars in the sample. The jets have total velocities ranging from 147 to 184 km/s after correcting for disk inclination. For B335, an 8-month gap between NIRSpec and MIRI MRS observations enabled measurement of the tangential velocity of a shocked knot; in combination with the radial velocity, this shows that the jet has a different inclination than the outflow cavity. We find multiple knots before and during a recent outburst in B335, although the knots were more frequent during the burst. The asymmetries between blue- and red-shifted lobes strongly suggest complex interactions between the circumstellar disks and magnetic fields.

preprint2022arXiv

The Rate, Amplitude and Duration of Outbursts from Class 0 Protostars in Orion

At least half of a protostar's mass is accreted in the Class 0 phase, when the central protostar is deeply embedded in a dense, infalling envelope. We present the first systematic search for outbursts from Class 0 protostars in the Orion clouds. Using photometry from Spitzer/IRAC spanning 2004 to 2017, we detect three outbursts from Class 0 protostars with $\ge 2$ mag changes at 3.6 or 4.5 $μ$m. This is comparable to the magnitude change of a known protostellar FU Ori outburst. Two are newly detected bursts from the protostars HOPS 12 and 124. The number of detections implies that Class 0 protostars burst every 438 yr, with a 95% confidence interval of 161 to 1884 yr. Combining Spitzer and WISE/NEOWISE data spanning 2004-2019, we show that the bursts persist for more than nine years with significant variability during each burst. Finally, we use $19-100$ $μ$m photometry from SOFIA, Spitzer and Herschel to measure the amplitudes of the bursts. Based on the burst interval, a duration of 15 yr, and the range of observed amplitudes, 3-100% of the mass accretion during the Class 0 phase occurs during bursts. In total, we show that bursts from Class 0 protostars are as frequent, or even more frequent, than those from more evolved protostars. This is consistent with bursts being driven by instabilities in disks triggered by rapid mass infall. Furthermore, we find that bursts may be a significant, if not dominant, mode of mass accretion during the Class 0 phase.

preprint2021arXiv

A Census of Protostellar Outflows in Nearby Molecular Clouds

We adopt the deep learning method CASI-3D (Convolutional Approach to Structure Identification-3D) to systemically identify protostellar outflows in 12CO and 13CO observations of the nearby molecular clouds, Ophiuchus, Taurus, Perseus and Orion. The total outflow masses are 267 Msun, 795 Msun, 1305 Msun and 6332 Msun for Ophiuchus, Taurus, Perseus and Orion, respectively. We show the outflow mass in each cloud is linearly proportional to the total number of young stellar objects. The estimated total 3D deprojected outflow energies are 9e45 ergs, 6e46 ergs, 1.2e47 ergs and 6e47 ergs for Ophiuchus, Taurus, Perseus and Orion, respectively. The energy associated with outflows is sufficient to offset turbulent dissipation at the current epoch for all four clouds. All clouds also exhibit a break point in the spatial power spectrum of the outflow prediction map, which likely corresponds to the typical outflow mass and energy injection scale.

preprint2020arXiv

Application of Convolutional Neural Networks to Identify Protostellar Outflows in CO Emission

We adopt the deep learning method CASI-3D (Convolutional Approach to Structure Identification-3D) to identify protostellar outflows in molecular line spectra. We conduct magneto-hydrodynamics simulations that model forming stars that launch protostellar outflows and use these to generate synthetic observations. We apply the 3D radiation transfer code RADMC-3D to model 12CO (J=1-0) line emission from the simulated clouds. We train two CASI-3D models: ME1 is trained to predict only the position of outflows, while MF is trained to predict the fraction of the mass coming from outflows in each voxel. The two models successfully identify all 60 previously visually identified outflows in Perseus. Additionally, CASI-3D finds 20 new high-confidence outflows. All of these have coherent high-velocity structures, and 17 of them have nearby young stellar objects, while the remaining three are outside the Spitzer survey coverage. The mass, momentum and energy of individual outflows in Perseus predicted by model MF is comparable to the previous estimations. This similarity is due to a cancelation in errors: previous calculations missed outflow material with velocities comparable to the cloud velocity, however, they compensate for this by over-estimating the amount of mass at higher velocities that has contamination from non-outflow gas. We show outflows likely driven by older sources have more high-velocity gas compared to those driven by younger sources.

preprint2020arXiv

Application of Convolutional Neural Networks to Identify Stellar Feedback Bubbles in CO Emission

We adopt the deep learning method CASI (Convolutional Approach to Shell Identification) and extend it to 3D (CASI-3D) to identify signatures of stellar feedback in molecular line spectra, such as 13CO. We adopt magneto-hydrodynamics simulations that study the impact of stellar winds in a turbulent molecular cloud as an input to generate synthetic observations. We apply the 3D radiation transfer code radmc-3d to model 13CO (J=1-0) line emission from the simulated clouds. We train two CASI-3d models: ME1 predicts only the position of feedback, while MF predicts the fraction of the mass coming from feedback in each voxel. We adopt 75% of the synthetic observations as the training set and assess the accuracy of the two models with the remaining data. We demonstrate that model ME1 identifies bubbles in simulated data with 95% accuracy, and model MF predicts the bubble mass within 4% of the true value. We use bubbles previously visually identified in Taurus in 13CO to validate the models and show both perform well on the highest confidence bubbles. We apply our two models on the full 98 deg2 FCRAO 13CO survey of the Taurus cloud. Models ME1 and MF predict feedback gas mass of 2894 M and 302 M, respectively. When including a correction factor for missing energy due to the limited velocity range of the 13CO data cube, model ME1 predicts feedback kinetic energies of 4.0*1e46 ergs and 1.5*1e47 ergs with/without subtracting the cloud velocity gradient. Model MF predicts feedback kinetic energy of 9.6*1e45 ergs and 2.8*1e46 ergs with/without subtracting the cloud velocity gradient. Model ME1 predicts bubble locations and properties consistent with previous visually identified bubbles. However, model MF demonstrates that feedback properties computed based on visual identifications are significantly overestimated due to line of sight confusion and contamination from background and foreground gas.

preprint2016arXiv

Spitzer Observations of Long Term Infrared Variability Among Young Stellar Objects in Chamaeleon I

Infrared variability is common among young stellar objects, with surveys finding daily to weekly fluctuations of a few tenths of a magnitude. Space-based observations can produce highly sampled infrared light curves, but are often limited to total baselines of about a month due to the orientation of the spacecraft. Here we present observations of the Chameleon I cluster whose low declination makes it observable by the Spitzer space telescope over a 200 day period. We observe 30 young stellar objects with a daily cadence to better sample variability on timescales of months. We find such variability is common, occurring in ~80% of the detected cluster members. The change in [3.6]-[4.5] color over 200 days for many of the sources falls between that expected for extinction and fluctuations in disk emission. With our high cadence and long baseline we can derive power spectral density curves covering two orders of magnitude in frequency and find significant power at low frequencies, up to the boundaries of our 200 day survey. Such long timescales are difficult to explain with variations driven by the interaction between the disk and stellar magnetic field, which has a dynamical timescale of days to weeks. The most likely explanation is either structural or temperature fluctuations spread throughout the inner ~0.5 au of the disk, suggesting that the intrinsic dust structure is highly dynamic.

preprint2015arXiv

Radio Properties of Young Stellar Objects in the Core of the Serpens South Infrared Dark Cloud

We present deep radio continuum observations of the star-forming core of the Serpens South Infrared Dark Cloud with the Karl G. Jansky Very Large Array (VLA). Observations were conducted in two bands centered at 7.25 GHz (4.14 cm) and 4.75 GHz (6.31 cm) with an rms of 8.5 and 11.1 microJy/beam, respectively. We also use 2MASS, Spitzer and Herschel data to put our radio observations in the context of young stellar populations characterized by near and far infrared observations. Within a 5 arcmin x 5 arcmin region of interest around the central cluster, we detect roughly eighteen radio sources, seven of which we determine are protostellar in nature due to their radio spectral indices and their association with infrared sources. We find evidence for a previously undetected embedded Class 0 protostar and reaffirm Class 0 protostellar classifications determined by previous millimeter wavelength continuum studies. We use our infrared data to derive mid-infrared luminosities for three of our protostellar sources and find relative agreement between the known YSO radio luminosity vs bolometric luminosity correlation. Lastly, we marginally detect an additional six radio sources at the 2-3 sigma level that lie within two arcseconds of infrared YSO candidates, providing motivation for higher sensitivity studies to clarify the nature of these sources and further probe embedded and/or low luminosity YSOs in Serpens South.

preprint2014arXiv

Connecting X-ray and Infrared Variability among Young Stellar Objects: Ruling out potential sources of disk fluctuations

Variability in the infrared emission from disks around pre-main sequence stars over the course of days to weeks appears to be common, but the physical cause of the changes in disk structure are not constrained. Here we present coordinated monitoring of one young cluster with the Spitzer and Chandra space telescopes aimed at studying the physical source of the variability. In fall 2011 we obtained ten epochs of Chandra ACIS photometry over a period of 30 days with a roughly 3 day cadence contemporaneous with 20 epochs of Spitzer [3.6],[4.5] photometry over 40 days with a roughly 2 day cadence of the IC 348 cluster. This cadence allows us to search for week to month long responses of the infrared emission to changes in the high-energy flux. We find no strong evidence for a direct link between the X-ray and infrared variability on these timescales among 39 cluster members with circumstellar disks. There is no significant correlation between the shape of the infrared and X-ray light curves, or between the size of the X-ray and infrared variability. Among the stars with an X-ray flare none showed evidence of a correlated change in the infrared photometry on timescales of days to weeks following the flare. This lack of connection implies that X-ray heating of the planet forming region of the disk is not significant, although we cannot rule out rapid or instantaneous changes in infrared emission.

preprint2014arXiv

CSI 2264: Simultaneous optical and infrared light curves of young disk-bearing stars in NGC 2264 with CoRoT and Spitzer-- evidence for multiple origins of variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$α$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

preprint2014arXiv

Star Cluster Formation and Feedback

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

preprint2013arXiv

Pulsed Accretion in a Variable Protostar

Periodic increases in luminosity arising from variable accretion rates have been predicted for some close pre-main sequence binary stars as they grow from circumbinary disks. The phenomenon is known as "pulsed accretion" and can affect the orbital evolution and mass distribution of young binaries, as well as the potential for planet formation in the circumbinary environment. Accretion variability is a common feature of young stars, with a large range of amplitudes and timescales as measured from multi-epoch observations at optical and infrared wavelengths. Periodic variations consistent with pulsed accretion have been seen in only a few young binaries via optical accretion tracers, albeit intermittently with accretion luminosity variations ranging from 0-50 percent from orbit to orbit. Here we report on a young protostar (age ~10^5 yr) that exhibits periodic variability in which the infrared luminosity increases by a factor of 10 in roughly one week every 25.34 days. We attribute this to pulsed accretion associated with an unseen binary companion. The strength and regularity of this accretion signal is surprising; it may be related to the very young age of the system, which is a factor of ten younger than the other pulsed accretors previously studied.

preprint2012arXiv

Infrared Variability of Evolved Protoplanetary Disks: Evidence for Scale Height Variations in the Inner Disk

We present the results of a multi-wavelength multi-epoch survey of five evolved protoplanetary disks in the IC 348 cluster that show significant infrared variability. Using 3-8micron and 24micron photometry along with 5-40micron spectroscopy from the Spitzer Space Telescope, as well as ground-based 0.8-5micron spectroscopy, optical spectroscopy and near-infrared photometry, covering timescales of days to years, we examine the variability in the disk, stellar and accretion flux. We find substantial variations (10-60%) at all infrared wavelengths on timescales of weeks to months for all of these young stellar objects. This behavior is not unique when compared to other cluster members and is consistent with changes in the structure of the inner disk, most likely scale height fluctuations on a dynamical timescale. Previous observations, along with our near-infrared photometry, indicate that the stellar fluxes are relatively constant; stellar variability does not appear to drive the large changes in the infrared fluxes. Based on our near-infrared spectroscopy of the Pa-beta and Br-gamma lines we find that the accretion rates are variable in most of the evolved disks but the overall rates are probably too small to cause the infrared variability. We discuss other possible physical causes for the variability, including the influence of a companion, magnetic fields threading the disk, and X-ray flares.

preprint2012arXiv

Potential Drivers of Mid-Infrared Variability in Young Stars: testing physical models with multi-epoch near-infrared spectra of YSOs in ρ Oph

Recent studies have identified several young stellar objects (YSOs) which exhibit significant mid-infrared (mid-IR) variability. A wide range of physical mechanisms may be responsible for these variations, including changes in a YSO's accretion rate or in the extinction or emission from the inner disk. We have obtained and analyzed multi-epoch near-infrared (NIR) spectra for five actively accreting YSOs in the $ρ$ Oph star-forming region along with contemporaneous mid-IR light curves obtained as part of the YSOVAR Spitzer/IRAC survey. Four of the five YSOs exhibit mid-IR light curves with modest ($\sim 0.2$--0.4 mag) but statistically significant variations over our 40-day observation window. Measuring the strengths of prominent photospheric absorption lines and accretion sensitive \ion{H}{1} and \ion{He}{1} lines in each NIR spectrum, we derive estimates of each YSO's spectral type, effective temperature ({\Teff}), and $H$ band extinction ($A_H$), and analyze the time evolution of their NIR veiling ($r_H$ and $r_K$) and mass accretion rates ({\Macc}). Defining a YSO's evolutionary stage such that heavily veiled, high accretion rate objects are less evolved than those with lower levels of veiling and ongoing accretion, we infer that GY 314 is the most evolved YSO in our sample, with GY 308 and GY 292 at progressively earlier evolutionary stages. Leveraging our multi-epoch, multi-wavelength dataset, we detect significant variations in mass accretion rates over timescales of days to weeks, but find that extinction levels in these YSOs remain relatively constant. We find no correlation between these YSO mid-IR light curves and time-resolved veiling or mass accretion rates, such that we are unable to link their mid-IR variability with physical processes localized near the inner edge of the circumstellar disk or within regions which are directly responsive to mass accretion...

preprint2012arXiv

The Spitzer Survey of Interstellar Clouds in the Gould Belt. V. Ophiuchus North Observed with IRAC and MIPS

We present Spitzer IRAC (2.1 sq. deg.) and MIPS (6.5 sq. deg.) observations of star formation in the Ophiuchus North molecular clouds. This fragmentary cloud complex lies on the edge of the Sco-Cen OB association, several degrees to the north of the well-known rho Oph star-forming region, at an approximate distance of 130 pc. The Ophiuchus North clouds were mapped as part of the Spitzer Gould Belt project under the working name `Scorpius'. In the regions mapped, selected to encompass all the cloud with visual extinction AV>3, eleven Young Stellar Object (YSO) candidates are identified, eight from IRAC/MIPS colour-based selection and three from 2MASS K/MIPS colours. Adding to one source previously identified in L43 (Chen et al. 2009), this increases the number of YSOcs identified in Oph N to twelve. During the selection process, four colour-based YSO candidates were rejected as probable AGB stars and one as a known galaxy. The sources span the full range of YSOc classifications from Class 0/1 to Class III, and starless cores are also present. Twelve high-extinction (AV>10) cores are identified with a total mass of approx. 100 solar masses. These results confirm that there is little ongoing star formation in this region (instantaneous star formation efficiency <0.34%) and that the bottleneck lies in the formation of dense cores. The influence of the nearby Upper Sco OB association, including the 09V star zeta Oph, is seen in dynamical interactions and raised dust temperatures but has not enhanced levels of star formation in Ophiuchus North.

preprint2011arXiv

Do All Stars in the Solar Neighbourhood Form in Clusters?

We present a global study of low mass, young stellar object (YSO) surface densities (Sigma) in nearby (< 500 pc) star forming regions based on a comprehensive collection of Spitzer Space Telescope surveys. We show that the distribution of YSO surface densities is a smooth distribution, being adequately described by a lognormal function from a few to 10^3 YSOs per pc^2, with a peak at ~22 stars pc^-2. The observed lognormal Sigma is consistent with predictions of hierarchically structured star-formation at scales below 10 pc, arising from the molecular cloud structures. We do not find evidence for multiple discrete modes of star-formation (e.g. clustered and distributed). Comparing the observed Sigma distribution to previous Sigma threshold definitions of clusters show that they are arbitrary. We find that only a low fraction (< 26$) of stars are formed in dense environments where their formation/evolution (along with their circumstellar disks and/or planets) may be affected by the close proximity of their low-mass neighbours.

preprint2011arXiv

The Highly Dynamic Behavior of the Innermost Dust and Gas in the Transition Disk Variable LRLL 31

We describe extensive synoptic multi-wavelength observations of the transition disk LRLL 31 in the young cluster IC 348. We combined four epochs of IRS spectra, nine epochs of MIPS photometry, seven epochs of cold-mission IRAC photometry and 36 epochs of warm mission IRAC photometry along with multi-epoch near-infrared spectra, optical spectra and polarimetry to explore the nature of the rapid variability of this object. We find that the inner disk, as traced by the 2-5micron excess stays at the dust sublimation radius while the strength of the excess changes by a factor of 8 on weekly timescales, and the 3.6 and 4.5micron photometry shows a drop of 0.35 magnitudes in one week followed by a slow 0.5 magnitude increase over the next three weeks. The accretion rate, as measured by PaBeta and BrGamma emission lines, varies by a factor of five with evidence for a correlation between the accretion rate and the infrared excess. While the gas and dust in the inner disk are fluctuating the central star stays relatively static. Our observations allow us to put constraints on the physical mechanism responsible for the variability. The variabile accretion, and wind, are unlikely to be causes of the variability, but both are effects of the same physical process that disturbs the disk. The lack of periodicity in our infrared monitoring indicates that it is unlikely that there is a companion within ~0.4 AU that is perturbing the disk. The most likely explanation is either a companion beyond ~0.4 AU or a dynamic interface between the stellar magnetic field and the disk leading to a variable scale height and/or warping of the inner disk.

preprint2011arXiv

The Structure of the Star-forming Cluster RCW 38

We present a study of the structure of the high mass star-forming region RCW~38 and the spatial distribution of its young stellar population. Spitzer IRAC photometry 3-8um are combined with 2MASS near-IR data to identify young stellar objects by IR-excess emission from their circumstellar material. Chandra X-ray data are used to identify class III pre-main sequence stars lacking circumstellar material. We identify 624 YSOs: 23 class 0/I and 90 flat spectrum protostars, 437 Class II stars, and 74 Class III stars. We also identify 29 (27 new) O star candidates over the IRAC field. Seventy-two stars exhibit IR-variability, including seven class 0/I and 12 flat spectrum YSOs. A further 177 tentative candidates are identified by their location in the IRAC [3.6] vs. [3.6]-[5.8] cmd. We find strong evidence of subclustering in the region. Three subclusters were identified surrounding the central cluster, with massive and variable stars in each subcluster. The central region shows evidence of distinct spatial distributions of the protostars and pre-main sequence stars. A previously detected IR cluster, DB2001_Obj36, has been established as a subcluster of RCW 38. This suggests that star formation in RCW 38 occurs over a more extended area than previously thought. The gas to dust ratio is examined using the X-ray derived hydrogen column density, N_H and the K-band extinction, and found to be consistent with the diffuse ISM, in contrast with Serpens & NGC1333. We posit that the high photoionising flux of massive stars in RCW 38 affects the agglomeration of the dust grains.

preprint2010arXiv

Circumstellar Structure around Evolved Stars in the Cygnus-X Star Formation Region

We present observations of newly discovered 24 micron circumstellar structures detected with the Multiband Imaging Photometer for Spitzer (MIPS) around three evolved stars in the Cygnus-X star forming region. One of the objects, BD+43 3710, has a bipolar nebula, possibly due to an outflow or a torus of material. A second, HBHA 4202-22, a Wolf-Rayet candidate, shows a circular shell of 24 micron emission suggestive of either a limb-brightened shell or disk seen face-on. No diffuse emission was detected around either of these two objects in the Spitzer 3.6-8 micron Infrared Array Camera (IRAC) bands. The third object is the luminous blue variable candidate G79.29+0.46. We resolved the previously known inner ring in all four IRAC bands. The 24 micron emission from the inner ring extends ~1.2 arcmin beyond the shorter wavelength emission, well beyond what can be attributed to the difference in resolutions between MIPS and IRAC. Additionally, we have discovered an outer ring of 24 micron emission, possibly due to an earlier episode of mass loss. For the two shell stars, we present the results of radiative transfer models, constraining the stellar and dust shell parameters. The shells are composed of amorphous carbon grains, plus polycyclic aromatic hydrocarbons in the case of G79.29+0.46. Both G79.29+0.46 and HBHA 4202-22 lie behind the main Cygnus-X cloud. Although G79.29+0.46 may simply be on the far side of the cloud, HBHA 4202-22 is unrelated to the Cygnus-X star formation region.

preprint2008arXiv

The Spitzer c2d Survey of Nearby Dense Cores. V. Discovery of a VeLLO in the "Starless" Dense Core L328

This paper reports the discovery of a Very Low Luminosity Object (VeLLO) in the "starless" dense core L328, using the Spitzer Space Telescope and ground based observations from near-infrared to millimeter wavelengths. The Spitzer 8 micron image indicates that L328 consists of three subcores of which the smallest one may harbor a source, L328-IRS while two other subcores remain starless. L328-IRS is a Class 0 protostar according to its bolometric temperature (44 K) and the high fraction ~72 % of its luminosity emitted at sub-millimeter wavelengths. Its inferred "internal luminosity" (0.04 - 0.06 Lsun) using a radiative transfer model under the most plausible assumption of its distance as 200 pc is much fainter than for a typical protostar, and even fainter than other VeLLOs studied previously. Note, however, that its inferred luminosity may be uncertain by a factor of 2-3 if we consider two extreme values of the distance of L328-IRS (125 or 310 pc). Low angular resolution observations of CO do not show any clear evidence of a molecular outflow activity. But broad line widths toward L328, and Spitzer and near-infrared images showing nebulosity possibly tracing an outflow cavity, strongly suggest the existence of outflow activity. Provided that an envelope of at most ~0.1 Msunis the only mass accretion reservoir for L328-IRS, and the star formation efficiency is close to the canonical value ~30%, L328-IRS has not yet accreted more than 0.05 Msun. At the assumed distance of 200 pc, L328-IRS is destined to be a brown dwarf.