PolCat is a software package designed to generate and constrain models of triaxial, magnetized molecular cloud cores using submillimeter polarization and continuum intensity data. Such models often ...compare well with observational data, but they are mathematically degenerate. As a result, many allowed models are either very elongated or flattened, and far from virial equilibrium. We present a tensor virial analysis of PolCat models with the aim of developing an optional new regularization constraint that can be used on the fly during the PolCat modeling process. This constraint is intended to guide PolCat toward models that are in or close to equilibrium on all three principal axes. While we have found the expected families of spheroidal solutions in tensor virial balance, we have also found a population of triaxial cores that are in three-axis tensor virial equilibrium. We find that models generated using the virial constraint have much more realistic shapes, with very elongated or flattened solutions now absent. Thus, these show that the tensor virial constraint may be useful for restricting PolCat's solution space to more realistic models, thereby also reducing the degeneracy inherent in our approach. We also perform a Monte Carlo analysis to predict the range of projected axis ratios present in the population of equilibrium solutions to suggest how these cores may appear when observed.
In this paper, we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We ...demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N2H+ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialized. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N2H+ is a good tracer of the bound material of prestellar cores, although we find some evidence for N2H+ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C18O and that traced by N2H+, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor–Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from south-west to north-east, indicating sequential star formation across the region.
ABSTRACT The dust emissivity spectral index, β, is a critical parameter for deriving the mass and temperature of star-forming structures and, consequently, their gravitational stability. The β value ...is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present β, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting spectral energy distributions to combined Herschel and JCMT observations in the 160, 250, 350, 500, and 850 m bands. Most of the derived β and dust temperature values fall within the ranges of 1.0-2.7 and 8-20 K, respectively. In Perseus, we find the β distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant localized β variations within individual clumps and find low-β regions correlate with local temperature peaks, hinting at the possible origins of low-β grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.
ABSTRACT The equation of state of cold supra-nuclear-density matter, such as in neutron stars, is an open question in astrophysics. A promising method for constraining the neutron star equation of ...state is modeling pulse profiles of thermonuclear X-ray burst oscillations from hot spots on accreting neutron stars. The pulse profiles, constructed using spherical and oblate neutron star models, are comparable to what would be observed by a next-generation X-ray timing instrument like ASTROSAT, NICER, or a mission similar to LOFT. In this paper, we showcase the use of an evolutionary optimization algorithm to fit pulse profiles to determine the best-fit masses and radii. By fitting synthetic data, we assess how well the optimization algorithm can recover the input parameters. Multiple Poisson realizations of the synthetic pulse profiles, constructed with 1.6 million counts and no background, were fitted with the Ferret algorithm to analyze both statistical and degeneracy-related uncertainty and to explore how the goodness of fit depends on the input parameters. For the regions of parameter space sampled by our tests, the best-determined parameter is the projected velocity of the spot along the observer's line of sight, with an accuracy of ≤3% compared to the true value and with ≤5% statistical uncertainty. The next best determined are the mass and radius; for a neutron star with a spin frequency of 600 Hz, the best-fit mass and radius are accurate to ≤5%, with respective uncertainties of ≤7% and ≤10%. The accuracy and precision depend on the observer inclination and spot colatitude, with values of ∼1% achievable in mass and radius if both the inclination and colatitude are 60°.
We solve the focused transport equation of cosmic rays numerically to investigate non-isotropic models of the pitch-angle scattering coefficient. In previous work, the Fokker-Planck equation was ...solved either analytically by using approximations, or by using a numerical approach together with simple models for the pitch-angle scattering coefficient. It is the purpose of the current article so compute particle distribution functions as well as the parallel diffusion coefficient by solving numerically the focused transport equation for a more realistic Fokker-Planck coefficient of pitch-angle scattering. Our analytical form for the scattering parameter is based on non-linear diffusion theory that takes into account realistic scattering at pitch-angles close to 90°. This general form contains the isotropic form as well as the quasi-linear limit as special cases. We show that the ratio of the diffusion coefficients with and without focusing sensitively depends on the ratio of the turbulent magnetic field and the mean field. The assumed form of the pitch-angle Fokker-Planck coefficient has an influence on the parallel diffusion coefficient. In all considered cases we found a reduction of the ratio of the diffusion coefficients if the ratio of magnetic fields is reduced.
We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud (δ ≤ −5:31:27.5). Employing a two-step structure identification process, we construct ...individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 μm continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only ∼72 per cent of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 μm emission structures (column densities >3.7 × 1021 cm−2). The remaining 28 per cent of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of YSOs such that sources which are more evolved are associated with higher velocities.
We study the equilibrium of pressure truncated, filamentary molecular clouds that are threaded by rather general helical magnetic fields. We first apply the virial theorem to filamentary molecular ...clouds, including the effects of non-thermal motions and the turbulent pressure of the surrounding ISM. When compared with the data, we find that many filamentary clouds have a mass per unit length that is significantly reduced by the effects of external pressure, and that toroidal fields play a significant role in squeezing such clouds. We also develop exact numerical MHD models of filamentary molecular clouds with more general helical field configurations than have previously been considered. We examine the effects of the equation of state by comparing ‘isothermal’ filaments, with constant total (thermal plus turbulent) velocity dispersion, with equilibria constructed using a logatropic equation of state. Our theoretical models involve three parameters: two to describe the mass loading of the toroidal and poloidal fields, and a third that describes the radial concentration of the filament. We thoroughly explore our parameter space to determine which choices of parameters result in models that agree with the available observational constraints. We find that both equations of state result in equilibria that agree with the observational results. Moreover, we find that models with helical fields have more realistic density profiles than either unmagnetized models or those with purely poloidal fields; we find that most isothermal models have density distributions that fall off as r−1.8 to r−2, while logatropes have density profiles that range from r−1 to r−1.8. We find that purely poloidal fields produce filaments with steep radial density gradients that are not allowed by the observations.
ABSTRACT We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. ...We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 m map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 m peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 × 1023 cm−2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 × 1023 cm−2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023/2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.
We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of ...sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that if modelled as thermally supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival 13CO velocity dispersion measurements and find that our cores are typically pressure confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.
We present 450 and 850 m submillimeter continuum observations of the IC 5146 star-forming region taken as part of the James Clerk Maxwell Telescope Gould Belt Survey. We investigate the location of ...bright submillimeter (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimeter clumps are identified using FellWalker, and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5 to 116 with a mean mass of 8 and a median mass of 3.7 . A stability analysis for the clumps suggests that the majority are (thermally) Jeans stable, with . We further compare the locations of known protostars with the observed submillimeter emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimeter peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC 5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.