Engineering topics

Hartmann, Lee

Publications and source records attributed to Hartmann, Lee.

At least 37 records · Page 2

DOE OSTI2015

THE GOULD’S BELT VERY LARGE ARRAY SURVEY. II. THE SERPENS REGION

We present deep (∼17 μJy) radio continuum observations of the Serpens molecular cloud, the Serpens south cluster, and the W40 region obtained using the Very Large Array in its A configuration. We detect a total of 146 sources, 29 of which are young stellar objects (YSOs), 2 of which are BV stars, and 5 more of which are associated with phenomena related to YSOs. Based on their radio variability and spectral index, we propose that about 16 of the remaining 110 unclassified sources are also YSOs. For approximately 65% of the known YSOs detected here as radio sources, the emission is most likely non-thermal and related to stellar coronal activity. As also recently observed in Ophiuchus, our sample of YSOs with X-ray counterparts lies below the fiducial Güdel and Benz relation. Finally, we analyze the proper motions of nine sources in the W40 region. This allows us to better constrain the membership of the radio sources in the region.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

KINEMATIC AND SPATIAL SUBSTRUCTURE IN NGC 2264

We present an expanded kinematic study of the young cluster NGC 2264 based upon optical radial velocities measured using multi-fiber echelle spectroscopy at the 6.5 m MMT and Magellan telescopes. We report radial velocities for 695 stars, of which approximately 407 stars are confirmed or very likely members. Our results more than double the number of members with radial velocities from Fűrész et al., resulting in a much better defined kinematic relationship between the stellar population and the associated molecular gas. In particular, we find that there is a significant subset of stars that are systematically blueshifted with respect to the molecular ({sup 13}CO) gas. The detection of Lithium absorption and/or infrared excesses in this blueshifted population suggests that at least some of these stars are cluster members; we suggest some speculative scenarios to explain their kinematics. Our results also more clearly define the redshifted population of stars in the northern end of the cluster; we suggest that the stellar and gas kinematics of this region are the result of a bubble driven by the wind from O7 star S Mon. Our results emphasize the complexity of the spatial and kinematic structure of NGC 2264, important for eventually building up a comprehensive picture of cluster formation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

CSI 2264: CHARACTERIZING YOUNG STARS IN NGC 2264 WITH SHORT-DURATION PERIODIC FLUX DIPS IN THEIR LIGHT CURVES

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical CoRoT light curves exhibiting short-duration, shallow periodic flux dips. All of these stars have infrared excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHMs generally less than 1 day, depths almost always less than 15%, and periods (3 < P < 11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected in successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard “disk-locking” models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous Spitzer IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSOs in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

THE GOULD'S BELT VERY LARGE ARRAY SURVEY. IV. THE TAURUS-AURIGA COMPLEX

We present a multi-epoch radio study of the Taurus-Auriga star-forming complex made with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. We detect a total of 610 sources, 59 of which are related to young stellar objects (YSOs) and 18 to field stars. The properties of 56% of the young stars are compatible with non-thermal radio emission. We also show that the radio emission of more evolved YSOs tends to be more non-thermal in origin and, in general, that their radio properties are compatible with those found in other star-forming regions. By comparing our results with previously reported X-ray observations, we notice that YSOs in Taurus-Auriga follow a Güdel-Benz relation with κ = 0.03, as we previously suggested for other regions of star formation. In general, YSOs in Taurus-Auriga and in all the previous studied regions seem to follow this relation with a dispersion of ∼1 dex. Finally, we propose that most of the remaining sources are related with extragalactic objects but provide a list of 46 unidentified radio sources whose radio properties are compatible with a YSO nature.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

CSI 2264: Characterizing Young Stars in NGC 2264 With Short-Duration Periodic Flux Dips in Their Light Curves

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical CoRoT light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian corotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3 < P < 11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard \disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous Spitzer IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel- ow accretion columns arising near the inner disk wall.

CSI 2264↗

DOE OSTI2015

CHARACTERIZING THE YOUNGEST HERSCHEL-DETECTED PROTOSTARS. I. ENVELOPE STRUCTURE REVEALED BY CARMA DUST CONTINUUM OBSERVATIONS

We present Combined Array for Research in Millimeter-wave Astronomy 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample. These objects are characterized by very red 24-70 μm colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from four protostars and one starless core in the fields toward the PBRS; we also report one new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of ∼5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv distance, with more than 50% of the source emission arising from radii <1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with ρ ∝ R {sup –2.5}. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power law. The large amount of mass on scales <1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Accretion outbursts in self-gravitating protoplanetary disks

We improve on our previous treatments of the long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e., PdV work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead zone. The GI-induced spiral waves propagate well inside of the gravitationally unstable region before they trigger outbursts at R ≲ 1 AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local α treatments for GI. In our standard model where zero dead-zone residual viscosity (α{sub rd}) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as 0.01 in terms of the generic α parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite α{sub rd} triggers thermally driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when α{sub rd} = 0. The inclusion of non-zero residual viscosity in the dead zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a 'quiescent' accretion phase with M-dot {sub acc}∼10{sup −8}--10{sup −7} M{sub ⊙} yr{sup −1} over 60% of the time and spend less than 15% of the infall phase in accretion outbursts. While our models indicate that episodic mass accretion during protostellar evolution can qualitatively help explain the low accretion luminosities seen in most low-mass protostars, detailed tests of the mechanism will require model calculations for a range of protostellar masses with some constraint on the initial core angular momentum, which affects the length of time spent in a quasi-steady disk accretion phase.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

A spectroscopic census in young stellar regions: the σ Orionis cluster

We present a spectroscopic survey of the stellar population of the σ Orionis cluster. We have obtained spectral types for 340 stars. Spectroscopic data for spectral typing come from several spectrographs with similar spectroscopic coverage and resolution. More than half of the stars in our sample are members confirmed by the presence of lithium in absorption, strong Hα in emission or weak gravity-sensitive features. In addition, we have obtained high-resolution (R ∼ 34,000) spectra in the Hα region for 169 stars in the region. Radial velocities were calculated from this data set. The radial velocity distribution for members of the cluster is in agreement with previous work. Analysis of the profile of the Hα line and infrared observations reveals two binary systems or fast rotators that mimic the Hα width expected in stars with accretion disks. On the other hand, there are stars with optically thick disks and narrow Hα profiles not expected in stars with accretion disks. This contribution constitutes the largest homogeneous spectroscopic data set of the σ Orionis cluster to date.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The evolution of accretion in young stellar objects: Strong accretors at 3-10 Myr

While the rate of accretion onto T Tauri stars is predicted to decline with age, objects with strong accretion have been detected at ages of up to 10 Myr. We analyze a sample of these old accretors, identified by having a significant U band excess and infrared emission from a circumstellar disk. Objects were selected from the ∼3 Myr σ Ori, 4-6 Myr Orion OB1b, and 7-10 Myr Orion OB1a star forming associations. We use high-resolution spectra from the Magellan Inamori Kyocera Echelle to estimate the veiling of absorption lines and calculate extinction for our T Tauri sample. We also use observations obtained with the Magellan Echellette and, in a few cases, the SWIFT Ultraviolet and Optical Telescope to estimate the excess produced in the accretion shock, which is then fit with accretion shock models to estimate the accretion rate. We find that even objects as old as 10 Myr may have high accretion rates, up to ∼10{sup –8} M{sub ☉} yr{sup –1}. These objects cannot be explained by viscous evolution models, which would deplete the disk in shorter timescales unless the initial disk mass is very high, a situation that is unstable. We show that the infrared spectral energy distribution of one object, CVSO 206, does not reveal evidence of significant dust evolution, which would be expected during the 10 Myr lifetime. We compare this object to predictions from photoevaporation and planet formation models and suggest that neither of these processes have had a strong impact on the disk of CVSO 206.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The Gould's Belt very large array survey. III. The Orion region

We present results from a high-sensitivity (60 μJy), large-scale (2.26 deg{sup 2}) survey obtained with the Karl G. Jansky Very Large Array as part of the Gould's Belt Survey program. We detected 374 and 354 sources at 4.5 and 7.5 GHz, respectively. Of these, 148 are associated with previously known young stellar objects (YSOs). Another 86 sources previously unclassified at either optical or infrared wavelengths exhibit radio properties that are consistent with those of young stars. The overall properties of our sources at radio wavelengths such as their variability and radio to X-ray luminosity relation are consistent with previous results from the Gould's Belt Survey. Our detections provide target lists for follow-up Very Long Baseline Array radio observations to determine their distances as YSOs are located in regions of high nebulosity and extinction, making it difficult to measure optical parallaxes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

CSI 2264: characterizing accretion-burst dominated light curves for young stars in NGC 2264

Based on more than four weeks of continuous high-cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high-quality, multi-wavelength light curves for young stellar objects whose optical variability is dominated by short-duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief—several hours to one day—brightenings at optical and near-infrared wavelengths with amplitudes generally in the range of 5%-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a 30 day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u – g versus g – r color-color diagram with the largest UV excesses. These stars also have large Hα equivalent widths, and either centrally peaked, lumpy Hα emission profiles or profiles with blueshifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Among the stars with the largest UV excesses or largest Hα equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

HOPS 136: An edge-on orion protostar near the end of envelope infall

Edge-on protostars are valuable for understanding the disk and envelope properties of embedded young stellar objects, since the disk, envelope, and envelope cavities are all distinctly visible in resolved images and well constrained in modeling. Comparing Two Micron All Sky Survey, Wide-field Infrared Survey Explorer, Spitzer, Herschel, and APEX photometry and an IRAM limit from 1.2 to 1200 μm, Spitzer spectroscopy from 5 to 40 μm, and high-resolution Hubble imaging at 1.60 and 2.05 μm to radiative transfer modeling, we determine envelope and disk properties for the Class I protostar HOPS 136, an edge-on source in Orion's Lynds 1641 region. The source has a bolometric luminosity of 0.8 L {sub ☉}, a bolometric temperature of 170 K, and a ratio of submillimeter to bolometric luminosity of 0.8%. Via modeling, we find a total luminosity of 4.7 L {sub ☉} (larger than the observed luminosity due to extinction by the disk), an envelope mass of 0.06 M {sub ☉}, and a disk radius and mass of 450 AU and 0.002 M {sub ☉}. The stellar mass is highly uncertain but is estimated to fall between 0.4 and 0.5 M {sub ☉}. To reproduce the flux and wavelength of the near-infrared scattered-light peak in the spectral energy distribution, we require 5.4 × 10{sup –5} M {sub ☉} of gas and dust in each cavity. The disk has a large radius and a mass typical of more evolved T Tauri disks in spite of the significant remaining envelope. HOPS 136 appears to be a key link between the protostellar and optically revealed stages of star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

VLA and CARMA observations of protostars in the Cepheus clouds: Sub-arcsecond proto-binaries formed via disk fragmentation

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of λ = 7.3 mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of ∼0.''06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at λ = 1.3 mm with a best resolution of ∼0.''3 (100 AU) and the CB230 IRS1 observations were taken at λ = 3.4 mm with a best resolution of ∼3'' (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of ∼0.''3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R > 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with ∼20 AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii <20 AU.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THE GOULD's BELT VERY LARGE ARRAY SURVEY. I. THE OPHIUCHUS COMPLEX

We present large-scale (∼2000 arcmin{sup 2}), deep (∼20 μJy), high-resolution (∼1'') radio observations of the Ophiuchus star-forming complex obtained with the Karl G. Jansky Very Large Array at λ = 4 and 6 cm. In total, 189 sources were detected, 56 of them associated with known young stellar sources, and 4 with known extragalactic objects; the other 129 remain unclassified, but most of them are most probably background quasars. The vast majority of the young stars detected at radio wavelengths have spectral types K or M, although we also detect four objects of A/F/B types and two brown dwarf candidates. At least half of these young stars are non-thermal (gyrosynchrotron) sources, with active coronas characterized by high levels of variability, negative spectral indices, and (in some cases) significant circular polarization. As expected, there is a clear tendency for the fraction of non-thermal sources to increase from the younger (Class 0/I or flat spectrum) to the more evolved (Class III or weak line T Tauri) stars. The young stars detected both in X-rays and at radio wavelengths broadly follow a Güdel-Benz relation, but with a different normalization than the most radioactive types of stars. Finally, we detect a ∼70 mJy compact extragalactic source near the center of the Ophiuchus core, which should be used as gain calibrator for any future radio observations of this region.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THE LONG-TERM EVOLUTION OF PHOTOEVAPORATING PROTOPLANETARY DISKS

We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting for {approx}< 20% of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that the initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the M-dot-R{sub wall} plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of the order of 10 AU.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THE DEPENDENCE OF STAR FORMATION EFFICIENCY ON GAS SURFACE DENSITY

Studies by Lada et al. and Heiderman et al. have suggested that star formation mostly occurs above a threshold in gas surface density {Sigma} of {Sigma}{sub c} {approx} 120 M{sub Sun} pc{sup -2} (A{sub K} {approx} 0.8). Heiderman et al. infer a threshold by combining low-mass star-forming regions, which show a steep increase in the star formation rate per unit area {Sigma}{sub SFR} with increasing {Sigma}, and massive cores forming luminous stars which show a linear relation. We argue that these observations do not require a particular density threshold. The steep dependence of {Sigma}{sub SFR}, approaching unity at protostellar core densities, is a natural result of the increasing importance of self-gravity at high densities along with the corresponding decrease in evolutionary timescales. The linear behavior of {Sigma}{sub SFR} versus {Sigma} in massive cores is consistent with probing dense gas in gravitational collapse, forming stars at a characteristic free-fall timescale given by the use of a particular molecular tracer. The low-mass and high-mass regions show different correlations between gas surface density and the area A spanned at that density, with A {approx} {Sigma}{sup -3} for low-mass regions and A {approx} {Sigma}{sup -1} for the massive cores; this difference, along with the use of differing techniques to measure gas surface density and star formation, suggests that connecting the low-mass regions with massive cores is problematic. We show that the approximately linear relationship between dense gas mass and stellar mass used by Lada et al. similarly does not demand a particular threshold for star formation and requires continuing formation of dense gas. Our results are consistent with molecular clouds forming by galactic hydrodynamic flows with subsequent gravitational collapse.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

MODELING THE RESOLVED DISK AROUND THE CLASS 0 PROTOSTAR L1527

We present high-resolution sub/millimeter interferometric imaging of the Class 0 protostar L1527 IRS (IRAS 04368+2557) at {lambda} = 870 {mu}m and 3.4 mm from the Submillimeter Array and Combined Array for Research in Millimeter Astronomy. We detect the signature of an edge-on disk surrounding the protostar with an observed diameter of 180 AU in the sub/millimeter images. The mass of the disk is estimated to be 0.007 M{sub Sun }, assuming optically thin, isothermal dust emission. The millimeter spectral index is observed to be quite shallow at all the spatial scales probed: {alpha} {approx} 2, implying a dust opacity spectral index {beta} {approx} 0. We model the emission from the disk and surrounding envelope using Monte Carlo radiative transfer codes, simultaneously fitting the sub/millimeter visibility amplitudes, sub/millimeter images, resolved L' image, spectral energy distribution, and mid-infrared spectrum. The best-fitting model has a disk radius of R = 125 AU, is highly flared (H{proportional_to}R {sup 1.3}), has a radial density profile {rho}{proportional_to}R {sup -2.5}, and has a mass of 0.0075 M{sub Sun }. The scale height at 100 AU is 48 AU, about a factor of two greater than vertical hydrostatic equilibrium. The resolved millimeter observations indicate that disks may grow rapidly throughout the Class 0 phase. The mass and radius of the young disk around L1527 are comparable to disks around pre-main-sequence stars; however, the disk is considerably more vertically extended, possibly due to a combination of lower protostellar mass, infall onto the disk upper layers, and little settling of {approx}1 {mu}m-sized dust grains.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

A HERSCHEL AND APEX CENSUS OF THE REDDEST SOURCES IN ORION: SEARCHING FOR THE YOUNGEST PROTOSTARS

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument on board the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 {mu}m and 160 {mu}m that are either too faint (m{sub 24} > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 {mu}m band. We find that the 11 reddest protostar candidates with log {lambda}F{sub {lambda}}70/{lambda}F{sub {lambda}}24 > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 {mu}m fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log {lambda}F{sub {lambda}}70/{lambda}F{sub {lambda}}24 > 1.65; we name these sources 'PACS Bright Red sources', or PBRs. Our analysis reveals that the PBR sample is composed of Class 0 like sources characterized by very red spectral energy distributions (SEDs; T{sub bol} < 45 K) and large values of sub-millimeter fluxes (L{sub smm}/L{sub bol} > 0.6%). Modified blackbody fits to the SEDs provide lower limits to the envelope masses of 0.2-2 M{sub Sun} and luminosities of 0.7-10 L{sub Sun }. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates.

79 ASTRONOMY AND ASTROPHYSICS↗