We assess the status of big-bang nucleosynthesis (BBN) in light of the final Planck data release and other recent developments, and in anticipation of future measurements. Planck data fix the cosmic baryon density to 0.9% precision, and determine the helium abundance and effective number of neutrinos with precision approaching that of astronomical and BBN determinations respectively. In addition, new high-redshift measurements give D/H to better precision than theoretical predictions, and new Li/H data reconfirm the lithium problem. We present new ${}^{7}{\rm Be}(n,p){}^{7}{\rm Li}$ rates using new neutron capture measurements; we have also examined the effect of proposed changes in the $d(p,\gamma){}^{3}{\rm He}$ rates. Using these results we perform a series of likelihood analyses. We assess BBN/CMB consistency, with attention to how our results depend on the choice of Planck data, as well as how the results depend on the choice of non-BBN, non-Planck data sets. Most importantly the lithium problem remains, and indeed is more acute given the very tight D/H observational constraints; new neutron capture data reveals systematics that somewhat increases uncertainty and thus slightly reduces but does not essentially change the problem. We confirm that $d(p,\gamma){}^{3}{\rm He}$ theoretical rates brings D/H out of agreement and slightly increases 7Li; new experimental data are needed at BBN energies. Setting the lithium problem aside, we find the effective number of neutrino species at BBN is $N_\nu = 2.86 \pm 0.15$. Future CMB Stage-4 measurements promise substantial improvements in BBN parameters: helium abundance determinations will be competitive with the best astronomical determinations, and $N_{\rm eff}$ will approach sensitivities capable of detecting the effects of Standard Model neutrino heating of the primordial plasma. (Abridged)

This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy.

LOX is a lunar-orbiting astrophysics mission that will probe the cosmos at MeV energies. It is guided by open questions regarding thermonuclear, or Type-Ia, supernovae (SNeIa) and will characterize these inherently radioactive objects by enabling a systematic survey of SNeIa at gamma-ray energies for the first time. Astronomical investigations from lunar orbit afford new opportunities to advance our understanding of the cosmos. The foundation of LOX is an observational approach well suited to the all-sky monitoring demands of supernova investigations and time-domain astronomy. Its inherently wide field-of-view and continuous all-sky monitoring provides an innovative way of addressing decadal survey questions at MeV energies (0.1-10 MeV). The LOX approach achieves high sensitivity with a simple, high-heritage instrument design that eliminates the need for complex, position-sensitive detectors, kinematic event reconstruction, masks, or other insensitive detector mass, while also mitigating technology development, implementation complexity, and their associated costs. LOX can be realized within existing programs, like Explorer.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

We present the first spectral features obtained from Cassini's Composite Infrared Spectrometer (CIRS) for any icy moon. The spectral region covered by CIRS focal planes (FP) 3 and 4 is rich in emissivity features, but previous studies at these wavelengths have been limited by low signal to noise ratios (S/Rs) for individual spectra. Our approach is to average CIRS FP3 spectra to increase the S/R and use emissivity spectra to constrain the composition of the dark material on Iapetus. We find an emissivity feature at ~855 cm-1 and a possible doublet at 660 and 690 cm-1 that do not correspond to any known instrument artifacts. We attribute the 855 cm-1 feature to fine-grained silicates, similar to those found in dust on Mars and in meteorites, which are nearly featureless at shorter wavelengths. Silicates on the dark terrains of Saturn's icy moons have been suspected for decades, but there have been no definitive detections until now. Serpentines reported in the literature at ambient temperature and pressure have features near 855 and 660 cm-1. However, peaks can shift depending on temperature and pressure, so measurements at Iapetus-like conditions are necessary for more positive feature identifications. As a first investigation, we measured muscovite at 125K in a vacuum and found that this spectrum does match the emissivity feature near 855 cm-1 and the location of the doublet. Further measurements are needed to robustly identify a specific silicate, which would provide clues regarding the origin and implications of the dark material.

The Spitzer Space Telescope mapped the Perseus molecular cloud complex with IRAC and MIPS as part of the c2d Spitzer Legacy project. This paper combines the observations from both instruments giving an overview of low-mass star formation across Perseus from 3.6 to 70 micron. We provide an updated list of young stellar objects with new classifications and source fluxes from previous works, identifying 369 YSOs in Perseus with the Spitzer dataset. By synthesizing the IRAC and MIPS maps of Perseus and building on the work of previous papers in this series (Jorgensen et al. 2006, Rebull et al. 2007), we present a current census of star formation across the cloud and within smaller regions. 67% of the YSOs are associated with the young clusters NGC 1333 and IC 348. The majority of the star formation activity in Perseus occurs in the regions around the clusters, to the eastern and western ends of the cloud complex. The middle of the cloud is nearly empty of YSOs despite containing regions of high visual extinction. The western half of Perseus contains three-quarters of the total number of embedded YSOs (Class 0+I and Flat SED sources) in the cloud and nearly as many embedded YSOs as Class II and III sources. Class II and III greatly outnumber Class 0+I objects in eastern Perseus and IC 348. These results are consistent with previous age estimates for the clusters. Across the cloud, 56% of YSOs and 91% of the Class 0+I and Flat sources are in areas where Av > 5 mag, indicating a possible extinction threshold for star formation.

Solar filaments exhibit a range of eruptive-like dynamic activity, ranging from the full or partial eruption of the filament mass and surrounding magnetic structure as a coronal mass ejection (CME), to a fully confined or 'failed' eruption. On 2011 June 7, a dramatic partial eruption of a filament was observed by multiple instruments on SDO and STEREO. One of the interesting aspects of this event is the response of the solar atmosphere as non-escaping material falls inward under the influence of gravity. The impact sites show clear evidence of brightening in the observed EUV wavelengths due to energy release. Two plausible physical mechanisms explaining the brightening are considered: heating of the plasma due to the kinetic energy of impacting material compressing the plasma, or reconnection between the magnetic field of low-lying loops and the field carried by the impacting material. By analyzing the emission of the brightenings in several SDO/AIA wavelengths, and comparing the kinetic energy of the impacting material (7.6 x 10^26 - 5.8 x 10^27 ergs) to the radiative energy (1.9 x 10^25 - 2.5 x 10^26 ergs) we find the dominant mechanism of energy release involved in the observed brightening is plasma compression.

Solar flare X-ray emission results from rapidly increasing temperatures and emission measures in flaring active region loops. To date, observations from the X-Ray Sensor (XRS) onboard the Geostationary Operational Environmental Satellite (GOES) have been used to derive these properties, but have been limited by a number of factors, including the lack of a consistent background subtraction method capable of being automatically applied to large numbers of flares. In this paper, we describe an automated temperature and emission measure-based background subtraction method (TEBBS), which builds on the methods of Bornmann (1990). Our algorithm ensures that the derived temperature is always greater than the instrumental limit and the pre-flare background temperature, and that the temperature and emission measure are increasing during the flare rise phase. Additionally, TEBBS utilizes the improved estimates of GOES temperatures and emission measures from White et al. (2005). TEBBS was successfully applied to over 50,000 solar flares occurring over nearly three solar cycles (1980-2007), and used to create an extensive catalog of the solar flare thermal properties. We confirm that the peak emission measure and total radiative losses scale with background subtracted GOES X-ray flux as power-laws, while the peak temperature scales logarithmically. As expected, the peak emission measure shows an increasing trend with peak temperature, although the total radiative losses do not. While these results are comparable to previous studies, we find that flares of a given GOES class have lower peak temperatures and higher peak emission measures than previously reported. The resulting TEBBS database of thermal flare plasma properties is publicly available on Solar Monitor (www.solarmonitor.org/TEBBS/) and will be available on Heliophysics Integrated Observatory (www.helio-vo.eu).

Sunspot numbers form a comprehensive, long-duration proxy of solar activity and have been used numerous times to empirically investigate the properties of the solar cycle. A number of correlations have been discovered over the 24 cycles for which observational records are available. Here we carry out a sophisticated statistical analysis of the sunspot record that reaffirms these correlations, and sets up an empirical predictive framework for future cycles. An advantage of our approach is that it allows for rigorous assessment of both the statistical significance of various cycle features and the uncertainty associated with predictions. We summarize the data into three sequential relations that estimate the amplitude, duration, and time of rise to maximum for any cycle, given the values from the previous cycle. We find that there is no indication of a persistence in predictive power beyond one cycle, and conclude that the dynamo does not retain memory beyond one cycle. Based on sunspot records up to October 2011, we obtain, for Cycle 24, an estimated maximum smoothed monthly sunspot number of 97 +- 15, to occur in January--February 2014 +- 6 months.

Coronal mass ejections (CMEs) are large-scale eruptions of plasma and magnetic feld that can produce adverse space weather at Earth and other locations in the Heliosphere. Due to the intrinsic multiscale nature of features in coronagraph images, wavelet and multiscale image processing techniques are well suited to enhancing the visibility of CMEs and supressing noise. However, wavelets are better suited to identifying point-like features, such as noise or background stars, than to enhancing the visibility of the curved form of a typical CME front. Higher order multiscale techniques, such as ridgelets and curvelets, were therefore explored to characterise the morphology (width, curvature) and kinematics (position, velocity, acceleration) of CMEs. Curvelets in particular were found to be well suited to characterising CME properties in a self-consistent manner. Curvelets are thus likely to be of benefit to autonomous monitoring of CME properties for space weather applications.

Recently observed oscillations in the solar atmosphere have been interpreted and modeled as magnetohydrodynamic wave modes. This has allowed the estimation of parameters that are otherwise hard to derive, such as the coronal magnetic-field strength. This work crucially relies on the initial detection of the oscillations, which is commonly done manually. The volume of Solar Dynamics Observatory (SDO) data will make manual detection inefficient for detecting all of the oscillating regions. An algorithm is presented which automates the detection of areas of the solar atmosphere that support spatially extended oscillations. The algorithm identifies areas in the solar atmosphere whose oscillation content is described by a single, dominant oscillation within a user-defined frequency range. The method is based on Bayesian spectral analysis of time-series and image filtering. A Bayesian approach sidesteps the need for an a-priori noise estimate to calculate rejection criteria for the observed signal, and it also provides estimates of oscillation frequency, amplitude and noise, and the error in all these quantities, in a self-consistent way. The algorithm also introduces the notion of quality measures to those regions for which a positive detection is claimed, allowing simple post-detection discrimination by the user. The algorithm is demonstrated on two Transition Region and Coronal Explorer (TRACE) datasets, and comments regarding its suitability for oscillation detection in SDO are made.

Observational evidence is presented for the merging of a downward-propagating plasmoid with a looptop kernel during an occulted limb event on 2007 January 25. RHESSI lightcurves in the 9-18 keV energy range, as well as that of the 245 MHz channel of the Learmonth Solar Observatory, show enhanced nonthermal emission in the corona at the time of the merging suggesting that additional particle acceleration took place. This was attributed to a secondary episode of reconnection in the current sheet that formed between the two merging sources. RHESSI images were used to establish a mean downward velocity of the plasmoid of 12 km/s. Complementary observations from the SECCHI suite of instruments onboard STEREO-Behind showed that this process occurred during the acceleration phase of the associated CME. From wavelet-enhanced EUVI, images evidence of inflowing magnetic field lines prior to the CME eruption is also presented. The derived inflow velocity was found to be 1.5 km/s. This combination of observations supports a recent numerical simulation of plasmoid formation, propagation and subsequent particle acceleration due to the tearing mode instability during current sheet formation.

A long-standing problem in low-mass star formation is the "luminosity problem," whereby protostars are underluminous compared to the accretion luminosity expected both from theoretical collapse calculations and arguments based on the minimum accretion rate necessary to form a star within the embedded phase duration. Motivated by this luminosity problem, we present a set of evolutionary models describing the collapse of low-mass, dense cores into protostars, using the Young & Evans (2005) model as our starting point. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse in two dimensions. From the resulting spectral energy distributions, we calculate standard observational signatures to directly compare to observations. We incorporate several modifications and additions to the original Young & Evans model in an effort to better match observations with model predictions. We find that scattering, 2-D geometry, mass-loss, and outflow cavities all affect the model predictions, as expected, but none resolve the luminosity problem. A cycle of episodic mass accretion, however, can resolve this problem and bring the model predictions into better agreement with observations. Standard assumptions about the interplay between mass accretion and mass loss in our model give star formation efficiencies consistent with recent observations that compare the core mass function (CMF) and stellar initial mass function (IMF). The combination of outflow cavities and episodic mass accretion reduce the connection between observational Class and physical Stage to the point where neither of the two common observational signatures (bolometric temperature and ratio of bolometric to submillimeter luminosity) can be considered reliable indicators of physical Stage.

Observations of Lynds Dark Nebula 1221 from the Spitzer Space Telescope are presented. These data show three candidate protostars towards L1221, only two of which were previously known. The infrared observations also show signatures of outflowing material, an interpretation which is also supported by radio observations with the Very Large Array. In addition, molecular line maps from the Five College Radio Astronomy Observatory are shown. One-dimensional dust continuum modelling of two of these protostars, IRS1 and IRS3, is described. These models show two distinctly different protostars forming in very similar environments. IRS1 shows a higher luminosity and larger inner radius of the envelope than IRS3. The disparity could be caused by a difference in age or mass, orientation of outflow cavities, or the impact of a binary in the IRS1 core.

The diffuse morphology and transient nature of coronal mass ejections (CMEs) make them difficult to identify and track using traditional image processing techniques. We apply multiscale methods to enhance the visibility of the faint CME front. This enables an ellipse characterisation to objectively study the changing morphology and kinematics of a sample of events imaged by the Large Angle Spectrometric Coronagraph (LASCO) onboard the Solar and Heliospheric Observatory (SOHO) and the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) onboard the Solar Terrestrial Relations Observatory (STEREO). The accuracy of these methods allows us to test the CMEs for non-constant acceleration and expansion. We exploit the multiscale nature of CMEs to extract structure with a multiscale decomposition, akin to a Canny edge detector. Spatio-temporal filtering highlights the CME front as it propagates in time. We apply an ellipse parameterisation of the front to extract the kinematics (height, velocity, acceleration) and changing morphology (width, orientation). The kinematic evolution of the CMEs discussed in this paper have been shown to differ from existing catalogues. These catalogues are based upon running-difference techniques that can lead to over-estimating CME heights. Our resulting kinematic curves are not well-fitted with the constant acceleration model. It is shown that some events have high acceleration below $\sim$5 R$_{\sun}$. Furthermore, we find that the CME angular widths measured by these catalogues are over-estimated, and indeed for some events our analysis shows non-constant CME expansion across the plane-of-sky.

A new method for automated detection of polar coronal holes is presented. This method, called perimeter tracing, uses a series of 171, 195, and 304 \AA\ full disk images from the Extreme ultraviolet Imaging Telescope (EIT) on SOHO over solar cycle 23 to measure the perimeter of polar coronal holes as they appear on the limbs. Perimeter tracing minimizes line-of-sight obscurations caused by the emitting plasma of the various wavelengths by taking measurements at the solar limb. Perimeter tracing also allows for the polar rotation period to emerge organically from the data as 33 days. We have called this the Harvey rotation rate and count Harvey rotations starting 4 January 1900. From the measured perimeter, we are then able to fit a curve to the data and derive an area within the line of best fit. We observe the area of the northern polar hole area in 1996, at the beginning of solar cycle 23, to be about 4.2% of the total solar surface area and about 3.6% in 2007. The area of the southern polar hole is observed to be about 4.0% in 1996 and about 3.4% in 2007. Thus, both the north and south polar hole areas are no more than 15% smaller now than they were at the beginning of cycle 23. This compares to the polar magnetic field measured to be about 40% less now than it was a cycle ago.

Since Baade's photographic study of M32 in the mid 1940s, it has been accepted as an established fact that M32 is a compact dwarf satellite of M31. The purpose of this paper is to report on the findings of our investigation into the nature of the existing evidence. We find that the case for M32 being a satellite of M31 rests upon Hubble Space Telescope (HST) based stellar population studies which have resolved red-giant branch (RGB) and red clump stars in M32 as well as other nearby galaxies. Taken in isolation, this recent evidence could be considered to be conclusive in favour of the existing view. However, the conventional scenario does not explain M32's anomalously high central velocity dispersion for a dwarf galaxy (several times that of either NGC 147, NGC 185 or NGC 205) or existing planetary nebula observations (which suggest that M32 is more than twice as distant as M31) and also requires an elaborate physical explanation for M32's inferred compactness. Conversely, we find that the case for M32 being a normal galaxy, of the order of three times as distant as M31, is supported by: (1) a central velocity dispersion typical of intermediate galaxies, (2) the published planetary nebula observations, and (3) known scaling relationships for normal early-type galaxies. However, this novel scenario cannot account for the high apparent luminosities of the RGB stars resolved in the M32 direction by HST observations. We are therefore left with two apparently irreconcilable scenarios, only one of which can be correct, but both of which suffer from potentially fatal evidence to the contrary. This suggests that current understanding of some relevant fields is still very far from adequate.

Two different multi-resolution analyses are used to decompose the structure of active region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average and standard deviation of the magnetic flux gradient for alpha, beta, beta-gamma and beta-gamma-delta active regions increase in the order listed, and that the order is maintained over all length-scales. This study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length-scales in the active region, and not just those length-scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length-scales. It is also shown that the average gradient content of active regions that have large flares (GOES class 'M' and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length-scales.

Coronal Mass Ejections (CMEs) are challenging objects to detect using automated techniques, due to their high velocity and diffuse, irregular morphology. A necessary step to automating the detection process is to first remove the subjectivity introduced by the observer used in the current, standard, CME detection and tracking method. Here we describe and demonstrate a multiscale edge detection technique that addresses this step and could serve as one part of an automated CME detection system. This method provides a way to objectively define a CME front with associated error estimates. These fronts can then be used to extract CME morphology and kinematics. We apply this technique to a CME observed on 18 April 2000 by the Large Angle Solar COronagraph experiment (LASCO) C2/C3 and a CME observed on 21 April 2002 by LASCO C2/C3 and the Transition Region and Coronal Explorer (TRACE). For the two examples in this work, the heights determined by the standard manual method are larger than those determined with the multiscale method by approximately 10% using LASCO data and approximately 20% using TRACE data.

Observations by the Cores to Disk Legacy Team with the Spitzer Space Telescope have identified a low luminosity, mid-infrared source within the dense core, Lynds 1014, which was previously thought to harbor no internal source. Followup near-infrared and submillimeter interferometric observations have confirmed the protostellar nature of this source by detecting scattered light from an outflow cavity and a weak molecular outflow. In this paper, we report the detection of cm continuum emission with the VLA. The emission is characterized by a quiescent, unresolved 90 uJy 6 cm source within 0.2" of the Spitzer source. The spectral index of the quiescent component is $\alpha = 0.37\pm 0.34$ between 6 cm and 3.6 cm. A factor of two increase in 6 cm emission was detected during one epoch and circular polarization was marginally detected at the $5\sigma$ level with Stokes {V/I} $= 48 \pm 16$% . We have searched for 22 GHz H2O maser emission toward L1014-IRS, but no masers were detected during 7 epochs of observations between June 2004 and December 2006. L1014-IRS appears to be a low-mass, accreting protostar which exhibits cm emission from a thermal jet or a wind, with a variable non-thermal emission component. The quiescent cm radio emission is noticeably above the correlation of 3.6 cm and 6 cm luminosity versus bolometric luminosity, indicating more radio emission than expected. We characterize the cm continuum emission in terms of observations of other low-mass protostars, including updated correlations of centimeter continuum emission with bolometric luminosity and outflow force, and discuss the implications of recent larger distance estimates on the physical attributes of the protostar and dense molecular core.

We present the results of a submillimeter survey of 53 low-mass dense cores with the Submillimeter High Angular Resolution Camera II (SHARC-II). The survey is a follow-up project to the Spitzer Legacy Program ``From Molecular Cores to Planet-Forming Disks'', with the purpose being to create a complete data set of nearby low-mass dense cores from the infrared to the millimeter. We present maps of 52 cores at 350 microns and three cores at 450 microns, two of which were observed at both wavelengths. Of these 52 cores, 41 were detected by SHARC-II: 32 contained one submillimeter source while 9 contained multiple sources. For each submillimeter source detected, we report various source properties including source position, fluxes in various apertures, size, aspect ratio, and position angle. For the 12 cores that were not detected we present upper limits. The sources detected by SHARC-II have, on average, smaller sizes at the 2sigma contours than those derived from longer-wavelength bolometer observations. We conclude that this is not caused by a failure to integrate long enough to detect the full extent of the core; instead it arises primarily from the fact that the observations presented in this survey are insensitive to smoothly varying extended emission. We find that SHARC-II observations of low-mass cores are much better suited to distinguishing between starless and protostellar cores than observations at longer wavelengths. Very Low Luminosity Objects, a new class of objects being discovered by the Spitzer Space Telescope in cores previously classified as starless, look very similar at 350 microns to other cores with more luminous protostars.

Infrared images of the dark cloud core B59 were obtained with the Spitzer Space Telescope as part of the "Cores to Disks" Legacy Science project. Photometry from 3.6-70 microns indicates at least 20 candidate low-mass young stars near the core, more than doubling the previously known population. Out of this group, 13 are located within about 0.1 pc in projection of the molecular gas peak, where a new embedded source is detected. Spectral energy distributions span the range from small excesses above photospheric levels to rising in the mid-infrared. One other embedded object, probably associated with the millimeter source B59-MMS1, with a bolometric luminosity L(bol) roughly 2 L(sun), has extended structure at 3.6 and 4.5 microns, possibly tracing the edges of an outflow cavity. The measured extinction through the central part of the core is A(V) greater than of order 45 mag. The B59 core is producing young stars with a high efficiency.

We present submillimeter observations of dark clouds that are part of the Spitzer Legacy Program, From Molecular Cores to Planet-Forming Disks (c2d). We used the Submillimetre Common User's Bolometer Array to map the regions observed by Spitzer by the c2d program to create a census of dense molecular cores including data from the infrared to the submillimeter. In this paper, we present the basic data from these observations: maps, fluxes, and source attributes. We also show data for an object just outside the Perseus cloud that was serendipitously observed in our program. We propose that this object is a newly discovered, evolved protostar.

We present Spitzer Space Telescope observations of the "evolved starless core" L1521F which reveal the presence of a very low luminosity object (L < 0.07 Lsun). The object, L1521F-IRS, is directly detected at mid-infrared wavelengths (>5 micron) but only in scattered light at shorter infrared wavelengths, showing a bipolar nebula oriented east-west which is probably tracing an outflow cavity. The nebula strongly suggests that L1521F-IRS is embedded in the L1521F core. Thus L1521F-IRS is similar to the recently discovered L1014-IRS and the previously known IRAM 04191 in its substellar luminosity and dense core environment. However these objects differ significantly in their core density, core chemistry, and outflow properties, and some may be destined to be brown dwarfs rather than stars.

We report the first detections of the Class 0 protostellar source IRAM 04191+1522 at wavelengths shortward of 60 microns with the Spitzer Space Telescope. We see extended emission in the Spitzer images that suggests the presence of an outflow cavity in the circumstellar envelope. We combine the Spitzer observations with existing data to form a complete dataset ranging from 3.6 to 1300 microns and use these data to construct radiative transfer models of the source. We conclude that the internal luminosity of IRAM 04191+1522, defined to be the sum of the luminosity from the internal sources (a star and a disk), is L_int = 0.08 +/- 0.04 L_sun, placing it among the lowest luminosity protostars known. Though it was discovered before the launch of the Spitzer Space Telescope, IRAM 04191+1522 falls within a new class of Very Low Luminosity Objects being discovered by Spitzer. Unlike the two other well-studied objects in this class, which are associated either with weak, compact outflows or no outflows at all, IRAM 04191+1522 has a well-defined molecular outflow with properties consistent with those expected based on relations derived from higher luminosity (L_int > 1 L_sun) protostars. We discuss the difficulties in understanding IRAM 04191+1522 in the context of the standard model of star formation, and suggest a possible explanation for the very low luminosity of this source.

We present evidence that the 8 $\mu$m (dust) and 24 $\mu$m luminosities of star-forming galaxies are both strongly correlated with their 1.4 GHz and H$\alpha$ luminosities over a range in luminosity of two-to-three orders of magnitude. At the bright end, the correlations are found to be essentially linear over a luminosity range of about two orders of magnitude (corresponding to star-formation rates of several-tenths to several tens of solar masses per year). However, at the faint end, there appears to be a slope change for dwarf galaxies, possibly due to the lower dust-to-gas ratios and lower metallicities of the dwarfs. The correlations suggest that PAH features and mid-IR continuum emissions are good measures of the star formation rates of galaxies, and we present calibrations of star-formation rates based on existing radio and H$\alpha$ relations. Our findings are based on a sample of star-forming galaxies selected from the main field of the Spitzer First Look Survey with the aid of spectroscopic data from the Sloan Digital Sky Survey and VLA 1.4 GHz data.

Using the Submillimeter Array we report the discovery of a compact low mass bipolar molecular outflow from L1014-IRS and confirm its association with the L1014 dense core at 200 pc. Consequently, L1014-IRS is the lowest luminosity (L \~0.09 Lsun) and perhaps the lowest mass source known to be driving a bipolar molecular outflow, which is one of the smallest known in size (~500 AU), mass (< 10^{-4} Msun), and energetics (e.g., force < 10^{-7} Msun km/s/yr).

We present an evolutionary picture of a forming star. We assume a singular, isothermal sphere as the initial state of the core that undergoes collapse as described by \citet{shu77}. We include the evolution of a first hydrostatic core at early times and allow a disk to grow as predicted by \citet{adams86}. We use a 1-dimensional radiative transfer code to calculate the spectral energy distribution for the evolving protostar from the beginning of collapse to the point when all envelope material has accreted onto the star+disk system. Then, we calculate various observational signatures ($T_{bol}$, $L_{bol}/L_{smm}$, and infrared colors) as a function of time. As defined by the bolometric temperature criterion, the Class 0 stage should be very short, while the Class I stage persists for much of the protostar's early life. We present physical distinctions among the classes of forming stars and calculate the observational signatures for these classes. Finally, we present models of infrared color-magnitude diagrams, as observed by the Spitzer Space Telescope, that should be strong discriminators in determining the stage of evolution for a protostar.

Spitzer Space Telescope observations of a point-like source, L1014-IRS, close to the dust peak of the low-mass dense core L1014 have questioned its starless nature. The presence of an object with colors expected for an embedded protostar makes L1014-IRS the lowest luminosity isolated protostar known, and an ideal target with which to test star formation theories at the low mass end. In order to study its molecular content and to search for the presence of a molecular outflow, we mapped L1014 in at least one transition of 12CO, N2H+, HCO+, CS and of their isotopologues 13CO, C18O, C17O, N2D+ and H13CO+, using the FCRAO, the IRAM 30 meter and the CSO. The data show physical and chemical properties in L1014 typical of the less evolved starless cores: i.e. H2 central density of a few 10^5 molecules cm^-3, estimated mass of ~2M_sun, CO integrated depletion factor less than 10, N(N2H+)~6*10^12 cm^-2, N(N2D+)/N(N2H+) equal to 10% and relatively broad N2H+(1--0) lines (0.35 km/s). Infall signatures and significant velocity shifts between optically thick and optically thin tracers are not observed in the line profiles. No classical signatures of molecular outflow are found in the 12CO and 13CO observations. In particular, no high velocity wings are found, and no well-defined blue-red lobes of 12CO emission are seen in the channel maps. If sensitive, higher resolution observations confirm the absence of an outflow on a smaller scale than probed by our observations, L1014-IRS would be the only protostellar object known to be formed without driving an outflow.

We present observations of L1014, a dense core in the Cygnus region previously thought to be starless, but data from the Spitzer Space Telescope shows the presence of an embedded source. We propose a model for this source that includes a cold core, heated by the interstellar radiation field, and a low-luminosity internal source. The low luminosity of the internal source suggests a substellar object. If L1014 is representative, other "starless" cores may turn out to harbor central sources.

We present results from the observations and modeling of seventeen Class I cores with the Submillimetre Common Users Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). By modeling the transfer of radiation through the envelope for nine cores, we find, for a power law distribution n(r)=n_f(r/r_f)^-p, the average and standard deviation p=1.6 +/- 0.4 and a median of p=1.8. However, the inclusion of a disk or other point-like component can cause the derived p to be shallower by as much as 0.5. In addition, we test the Shu collapse model for our sources and discuss the application of simpler analyses that derive a density power law distribution directly from the slope of the intensity radial profile. The total mass of the envelope in our sample has a range from 0.04 to 5.0 M_sun, but these masses disagree with the virial masses derived from molecular line observations. Finally, we discuss the nature of these sources in light of various evolutionary indicators and find that T_bol and L_obs/L_smm are often inconsistent in distinguishing Class 0 from Class I cores.

Studying the physical environments of low mass and high mass cores using dust continuum emission provides important observational constraints on theoretical models of star formation. The motivation and procedure for modeling dust continuum emission is reviewed and the results of recent surveys towards low mass and high mass star forming regions are compared.

We present the results of a multi-faceted study of the high-mass star-forming region S235. With the CalTech Submillimeter Observatory (CSO), we have observed and/or mapped over fifty molecular transitions and 350 micron dust emission in an attempt to learn more about the chemistry, kinematics, and mass distribution within this rich region. We have found two dust cores that are not resolved in our molecular line observations and suggest that one of these is a less evolved star-forming core.

We present multi-frequency radio observations from the afterglow of GRB 980519 beginning 7.2 hours after the gamma-ray burst and ending 63 days later. The fast decline in the optical and X-ray light curves for this burst has been interpreted either as afterglow emission originating from a collimated outflow -- a jet -- or the result of a blast wave propagating into a medium whose density is shaped by the wind of an evolved massive star. These two models predict divergent behavior for the radio afterglow, and therefore, radio observations are capable, in principle, of discriminating between the two. We show that a wind model describes the subsequent evolution of the radio afterglow rather well. However, we see strong modulation of the light curve, which we interpret as diffractive scintillation. These variations prevent us from decisively rejecting the jet model.

GRB 990123 was the first burst from which simultaneous optical, X-ray and gamma-ray emission was detected; its afterglow has been followed by an extensive set of radio, optical and X-ray observations. We have studied the gamma-ray burst itself as observed by the CGRO detectors. We find that gamma-ray fluxes are not correlated with the simultaneous optical observations, and the gamma-ray spectra cannot be extrapolated simply to the optical fluxes. The burst is well fit by the standard four-parameter GRB function, with the exception that excess emission compared to this function is observed below ~15 keV during some time intervals. The burst is characterized by the typical hard-to-soft and hardness-intensity correlation spectral evolution patterns. The energy of the peak of the nu f_nu spectrum, E_p, reaches an unusually high value during the first intensity spike, 1470 +/- 110 keV, and then falls to \~300 keV during the tail of the burst. The high-energy spectrum above ~MeV is consistent with a power law with a photon index of about -3. By fluence, GRB 990123 is brighter than all but 0.4% of the GRBs observed with BATSE, clearly placing it on the -3/2 power-law portion of the intensity distribution. However, the redshift measured for the afterglow is inconsistent with the Euclidean interpretation of the -3/2 power-law. Using the redshift value of >= 1.61 and assuming isotropic emission, the gamma-ray fluence exceeds 10E54 ergs.

In a recent paper, Binggeli & Jerjen (1998) question the value of the extra- galactic distance indicators presented by Young & Currie (1994 & 1995) and state that they have refuted `the claim that the Virgo dEs [dwarf-elliptical galaxies]...are distributed in a prolate structure stretching from 8 to 20 Mpc distance (Young & Currie 1995).' even though no such claim was ever made. In this paper, we examine Binggeli & Jerjen's claims that intrinsic scatter rather than spatial depth must be the main cause of the large scatters observed in the relevant scaling relationships for Virgo galaxies. We investigate the accuracy of Binggeli & Jerjen's photometric parameters and find that while their profile curvature and scalelength measurements are probably useful, their total magnitude and central surface-brightness measurements are not useful for the purpose of investigating scaling laws because they suffer from serious systematic and random errors. We also investigate Binggeli & Jerjen's criti- cisms of our (1995) analysis. We demonstrate that their test for strong mutual dependence between distance estimates based on the two different scaling laws is invalid because of its prior assumption of negligible cluster depth. We further demonstrate that the [relative] distance estimates on which their kinematical arguments are based cannot be meaningful, not only because of the seriousness of the photometric errors, but also because they are undermined by the prior assumption that depth effects can again be neglected.

The `Virgo Photometry Catalogue' (VPC) is the first independently calibrated general catalogue of galaxies to cover the Virgo cluster since the `Catalog of Galaxies and Clusters of Galaxies' of Zwicky et al. (1961,1963). It contains 1180 galaxies (including background objects) within a 23 square-degree region centred on the cluster's core. Photographic surface photometry is presented for 1067 galaxies in the U band, for 1020 galaxies in the B_J band and for 1020 galaxies in the R_C band. All total magnitudes and total colours are extrapolated according to a new system; denoted `t' to distinguish it from the `T' system already in use. This paper outlines: the scope of the VPC, the new extrapolation system, some recent findings based on the catalogue and further work in progress.

We present a new, but simple, procedure for estimating the total magnitudes of galaxies. This procedure involves the out-focusing of digital galaxy images numerically, the fitting of the resulting surface-brightness profiles with a single generalised profile model and the extrapolation of the fitted profiles to infinite radial distances. This new system, which we denote `t', differs fundamentally from the `T' system (of the `Reference Catalog of Bright Galaxies' series) in that: (1) it enables a galaxy's luminosity profile to be extrapolated without the need for any prior morphological classification, and (2) it is applicable to images of widely different spatial resolutions (including unresolved ones) because it takes into account systematic effects due to differential image resolution. It also differs fundamentally from the Kron system in that: (1) it can be derived directly from surface photometry without the need to go back to the plate scans or CCD frames (unless the surface photometry is of high resolution and/or the galaxies being measured are very bright), and (2) it can cope with merged images (provided they are separable by image-segmentation software). Through worked examples, we demonstrate the stability of `t'-system total magnitudes with respect to morphological type, the seeing conditions at the time of observation, degree of smoothing and limiting isophote. We also compare and contrast the new system with both the `T' system and the Kron system, and investigate the advantages and limitations of each of the three systems.