Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PROTOSTARS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

On the Origin of C4H and CH3OH in Protostellar Envelopes

The formation pathways of different types of organic molecules in protostellar envelopes and other regions of star formation are subjects of intense current interest. We present here observations of C4H and CH3OH, tracing two distinct groups of interstellar organic molecules, toward 16 protostars in the Ophiuchus and Corona Australis molecular clouds. Together with observations in the literature, we present C4H and CH3OH data from single-dish observations of 40 embedded protostars. We nd no correlation between the C4H and CH3OH column densities in this large sample. Based on this lack of correlation, a difference in line proles between C4H and CH3OH, and previous interferometric observations of similar sources, we propose that the emission from these two molecules is spatially separated, with the CH3OH tracing gas that has been transiently heated to high (approx. 70-100 K) temperatures and the C4H tracing the cooler large-scale envelope where CH4 molecules have been liberated from ices. These results provide insight in the differentiation between hot corino and warm carbon-chain chemistry in embedded protostars.

Ophiuchus star-forming region↗

Mass assembly of stellar systems and their evolution with the SMA (MASSES). multiplicity and the physical environment in L1448N

We present continuum and molecular line observations at 230 and 345 GHz from the Submillimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project “Mass Assembly of Stellar Systems and their Evolution with the SMA.” Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C{sup 18}O (J=2→1), {sup 13}CO (J=2→1), and HCO{sup +} (J=4→3), and have offsets with N{sub 2}D{sup +} (J=3→2) peaks. High angular resolution data from a complementary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to <100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C{sup 18}O emission are found to be perpendicular to within 20° of the outflow directions as revealed by {sup 12}CO (J=2→1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

79 ASTRONOMY AND ASTROPHYSICS↗

DEEPLY EMBEDDED PROTOSTELLAR POPULATION IN THE 20 km s{sup −1} CLOUD OF THE CENTRAL MOLECULAR ZONE

We report the discovery of a population of deeply embedded protostellar candidates in the 20 km s{sup −1} cloud, one of the massive molecular clouds in the Central Molecular Zone (CMZ) of the Milky Way, using interferometric submillimeter continuum and H{sub 2}O maser observations. The submillimeter continuum emission shows five 1 pc scale clumps, each of which further fragments into several 0.1 pc scale cores. We identify 17 dense cores, among which 12 are gravitationally bound. Among the 18 H{sub 2}O masers detected, 13 coincide with the cores and probably trace outflows emanating from the protostars. There are also 5 gravitationally bound dense cores without H{sub 2}O maser detection. In total, the 13 masers and 5 cores may represent 18 protostars with spectral types later than B1 or potentially growing more massive stars at earlier evolutionary stages, given the non-detection in the centimeter radio continuum. In combination with previous studies of CH{sub 3}OH masers, we conclude that the star formation in this cloud is at an early evolutionary phase, before the presence of any significant ionizing or heating sources. Our findings indicate that star formation in this cloud may be triggered by a tidal compression as it approaches pericenter, similar to the case of G0.253+0.016 but with a higher star formation rate, and demonstrate that high angular resolution, high-sensitivity maser, and submillimeter observations are promising techniques to unveil deeply embedded star formation in the CMZ.

79 ASTRONOMY AND ASTROPHYSICS↗

The effects of magnetic fields and protostellar feedback on low-mass cluster formation

In this paper, we present a large suite of simulations of the formation of low-mass star clusters. Our simulations include an extensive set of physical processes – magnetohydrodynamics, radiative transfer, and protostellar outflows – and span a wide range of virial parameters and magnetic field strengths. Comparing the outcomes of our simulations to observations, we find that simulations remaining close to virial balance throughout their history produce star formation efficiencies and initial mass function (IMF) peaks that are stable in time and in reasonable agreement with observations. Our results indicate that small-scale dissipation effects near the protostellar surface provide a feedback loop for stabilizing the star formation efficiency. This is true regardless of whether the balance is maintained by input of energy from large-scale forcing or by strong magnetic fields that inhibit collapse. In contrast, simulations that leave virial balance and undergo runaway collapse form stars too efficiently and produce an IMF that becomes increasingly top heavy with time. In all cases, we find that the competition between magnetic flux advection towards the protostar and outward advection due to magnetic interchange instabilities, and the competition between turbulent amplification and reconnection close to newly formed protostars renders the local magnetic field structure insensitive to the strength of the large-scale field, ensuring that radiation is always more important than magnetic support in setting the fragmentation scale and thus the IMF peak mass. The statistics of multiple stellar systems are similarly insensitive to variations in the initial conditions and generally agree with observations within the range of statistical uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗

THE SPITZER SPACE TELESCOPE SURVEY OF THE ORION A AND B MOLECULAR CLOUDS. II. THE SPATIAL DISTRIBUTION AND DEMOGRAPHICS OF DUSTY YOUNG STELLAR OBJECTS

We analyze the spatial distribution of dusty young stellar objects (YSOs) identified in the Spitzer Survey of the Orion Molecular clouds, augmenting these data with Chandra X-ray observations to correct for incompleteness in dense clustered regions. We also devise a scheme to correct for spatially varying incompleteness when X-ray data are not available. The local surface densities of the YSOs range from 1 pc{sup −2} to over 10,000 pc{sup −2}, with protostars tending to be in higher density regions. This range of densities is similar to other surveyed molecular clouds with clusters, but broader than clouds without clusters. By identifying clusters and groups as continuous regions with surface densities ≥10 pc{sup −2}, we find that 59% of the YSOs are in the largest cluster, the Orion Nebula Cluster (ONC), while 13% of the YSOs are found in a distributed population. A lower fraction of protostars in the distributed population is evidence that it is somewhat older than the groups and clusters. An examination of the structural properties of the clusters and groups shows that the peak surface densities of the clusters increase approximately linearly with the number of members. Furthermore, all clusters with more than 70 members exhibit asymmetric and/or highly elongated structures. The ONC becomes azimuthally symmetric in the inner 0.1 pc, suggesting that the cluster is only ∼2 Myr in age. We find that the star formation efficiency (SFE) of the Orion B cloud is unusually low, and that the SFEs of individual groups and clusters are an order of magnitude higher than those of the clouds. Finally, we discuss the relationship between the young low mass stars in the Orion clouds and the Orion OB 1 association, and we determine upper limits to the fraction of disks that may be affected by UV radiation from OB stars or dynamical interactions in dense, clustered regions.

79 ASTRONOMY AND ASTROPHYSICS↗

THE JCMT GOULD BELT SURVEY: A FIRST LOOK AT DENSE CORES IN ORION B

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 × 10{sup 23} cm{sup −2}, most of the mass is found within dense cores, while at lower cloud column densities, below 1 × 10{sup 23} cm{sup −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.

79 ASTRONOMY AND ASTROPHYSICS↗

Outflow-confined H ii regions. I. First signposts of massive star formation

We present an evolutionary sequence of models of the photoionized disk-wind outflow around forming massive stars based on the Core Accretion model. The outflow is expected to be the first structure to be ionized by the protostar and can confine the expansion of the H ii region, especially in lateral directions in the plane of the accretion disk. The ionizing luminosity increases as Kelvin–Helmholz contraction proceeds, and the H ii region is formed when the stellar mass reaches ∼10–20M{sub ⊙} depending on the initial cloud core properties. Although some part of the outer disk surface remains neutral due to shielding by the inner disk and the disk wind, almost the whole of the outflow is ionized in 10{sup 3}–10{sup 4} yr after initial H ii region formation. Having calculated the extent and temperature structure of the H ii region within the immediate protostellar environment, we then make predictions for the strength of its free–free continuum and recombination line emission. The free–free radio emission from the ionized outflow has a flux density of ∼(20–200) × (ν/10GHz){sup p} mJy for a source at a distance of 1 kpc with a spectral index p ≃ 0.4–0.7, and the apparent size is typically ∼500 AU at 10 GHz. The H40α line profile has a width of about 100 km s{sup −1}. These properties of our model are consistent with observed radio winds and jets around forming massive protostars.

79 ASTRONOMY AND ASTROPHYSICS↗

FORMATION AND RECONDENSATION OF COMPLEX ORGANIC MOLECULES DURING PROTOSTELLAR LUMINOSITY OUTBURSTS

During the formation of stars, the accretion of surrounding material toward the central object is thought to undergo strong luminosity outbursts followed by long periods of relative quiescence, even at the early stages of star formation when the protostar is still embedded in a large envelope. We investigated the gas-phase formation and recondensation of the complex organic molecules (COMs) di-methyl ether and methyl formate, induced by sudden ice evaporation processes occurring during luminosity outbursts of different amplitudes in protostellar envelopes. For this purpose, we updated a gas-phase chemical network forming COMs in which ammonia plays a key role. The model calculations presented here demonstrate that ion–molecule reactions alone could account for the observed presence of di-methyl ether and methyl formate in a large fraction of protostellar cores without recourse to grain-surface chemistry, although they depend on uncertain ice abundances and gas-phase reaction branching ratios. In spite of the short outburst timescales of about 100 years, abundance ratios of the considered species higher than 10% with respect to methanol are predicted during outbursts due to their low binding energies relative to water and methanol which delay their recondensation during cooling. Although the current luminosity of most embedded protostars would be too low to produce complex organics in the hot-core regions that are observable with current sub-millimetric interferometers, previous luminosity outburst events would induce the formation of COMs in extended regions of protostellar envelopes with sizes increasing by up to one order of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗

A WISE CENSUS OF YOUNG STELLAR OBJECTS IN CANIS MAJOR

With the Wide-field Infrared Survey Explorer (WISE), we searched for young stellar objects (YSOs) in a 100 deg{sup 2} region centered on the lightly studied Canis Major star-forming region. Applying stringent magnitude cuts to exclude the majority of extragalactic contaminants, we find 144 Class I candidates and 335 Class II candidates. The sensitivity to Class II candidates is limited by their faintness at the distance to Canis Major (assumed as 1000 pc). More than half the candidates (53%) are found in 16 groups of more than four members, including four groups with more than 25 members each. The ratio of Class II to Class I objects, N {sub II}/ N {sub I}, varies from 0.4 to 8.3 in just the largest four groups. We compare our results to those obtainable with combined Two Micron All Sky Survey and post-cryogenic Spitzer Space Telescope data; the latter approach recovers missing Class II sources. Via a comparison to protostars characterized with the Herschel Space Observatory , we propose new WISE color criteria for flat-spectrum and Class 0 protostars, finding 80 and 7 of these, respectively. The distribution of YSOs in CMa OB1 is consistent with supernova-induced star formation, although the diverse N {sub II}/ N {sub I} ratios are unexpected if this parameter traces age and the YSOs are due to the same supernova. Less massive clouds feature larger N {sub II}/ N {sub I} ratios, suggesting that initial conditions play a role in determining this quantity.

79 ASTRONOMY AND ASTROPHYSICS↗

ALMA Cycle 1 observations of the HH 46/47 molecular outflow: structure, entrainment, and core impact

We present Atacama Large Millimeter/sub-millimeter Array Cycle 1 observations of the HH 46/47 molecular outflow using combined 12 m array and Atacama Compact Array observations. The improved angular resolution and sensitivity of our multi-line maps reveal structures that help us study the entrainment process in much more detail and allow us to obtain more precise estimates of outflow properties than in previous observations. We use {sup 13}CO (1–0) and C{sup 18}O (1–0) emission to correct for the {sup 12}CO (1–0) optical depth to accurately estimate the outflow mass, momentum, and kinetic energy. This correction increases the estimates of the mass, momentum, and kinetic energy by factors of about 9, 5, and 2, respectively, with respect to estimates assuming optically thin emission. The new {sup 13}CO and C{sup 18}O data also allow us to trace denser and slower outflow material than that traced by the {sup 12}CO maps, and they reveal an outflow cavity wall at very low velocities (as low as 0.2 km s{sup −1} with respect to the core’s central velocity). Adding the slower material traced only by {sup 13}CO and C{sup 18}O, there is another factor of three increase in the mass estimate and 50% increase in the momentum estimate. The estimated outflow properties indicate that the outflow is capable of dispersing the parent core within the typical lifetime of the embedded phase of a low-mass protostar and that it is responsible for a core-to-star efficiency of 1/4 to 1/3. We find that the outflow cavity wall is composed of multiple shells associated with a series of jet bow-shock events. Within about 3000 au of the protostar the {sup 13}CO and C{sup 18}O emission trace a circumstellar envelope with both rotation and infall motions, which we compare with a simple analytic model. The CS (2–1) emission reveals tentative evidence of a slowly moving rotating outflow, which we suggest is entrained not only poloidally but also toroidally by a disk wind that is launched from relatively large radii from the source.

79 ASTRONOMY AND ASTROPHYSICS↗

An Estimation of the Star Formation Rate in the Perseus Complex

We present the results of our investigation of the star-forming potential in the Perseus star-forming complex. We build on previous starless core, protostellar core, and young stellar object (YSO) catalogs from Spitzer (3.6–70 μm), Herschel (70–500 μm), and SCUBA (850 μm) observations in the literature. We place the cores and YSOs within seven star-forming clumps based on column densities greater than 5×10{sup 21} cm{sup −2}. We calculate the mean density and free-fall time for 69 starless cores as ∼5.55 ×10{sup −19} g cm{sup −3} and ∼0.1 Myr, respectively, and we estimate the star formation rate for the near future as ∼150 M {sub ⊙} Myr{sup −1}. According to Bonnor–Ebert stability analysis, we find that majority of starless cores in Perseus are unstable. Broadly, these cores can collapse to form the next generation of stars. We found a relation between starless cores and YSOs, where the numbers of young protostars (Class 0 + Class I) are similar to the numbers of starless cores. This similarity, which shows a one-to-one relation, suggests that these starless cores may form the next generation of stars with approximately the same formation rate as the current generation, as identified by the Class 0 and Class I protostars. It follows that if such a relation between starless cores and any YSO stage exists, the SFR values of these two populations must be nearly constant. In brief, we propose that this one-to-one relation is an important factor in better understanding the star formation process within a cloud.

79 ASTRONOMY AND ASTROPHYSICS↗

A Search for Water Maser Emission from Brown Dwarfs and Low-luminosity Young Stellar Objects

We present a survey for water maser emission toward a sample of 44 low-luminosity young objects, comprising (proto-)brown dwarfs, first hydrostatic cores (FHCs), and other young stellar objects (YSOs) with bolometric luminosities lower than 0.4 L {sub ⊙}. Water maser emission is a good tracer of energetic processes, such as mass-loss and/or accretion, and is a useful tool to study these processes with very high angular resolution. This type of emission has been confirmed in objects with L {sub bol} ≳ 1 L {sub ⊙}. Objects with lower luminosities also undergo mass-loss and accretion, and thus, are prospective sites of maser emission. Our sensitive single-dish observations provided a single detection when pointing toward the FHC L1448 IRS 2E. However, follow-up interferometric observations showed water maser emission associated with the nearby YSO L1448 IRS 2 (a Class 0 protostar of L {sub bol} ≃ 3.6–5.3 L {sub ⊙}) and did not find any emission toward L1448 IRS 2E. The upper limits for water maser emission determined by our observations are one order of magnitude lower than expected from the correlation between water maser luminosities and bolometric luminosities found for YSOs. This suggests that this correlation does not hold at the lower end of the (sub)stellar mass spectrum. Possible reasons are that the slope of this correlation is steeper at L {sub bol} ≤ 1 L {sub ⊙} or that there is an absolute luminosity threshold below which water maser emission cannot be produced. Alternatively, if the correlation still stands at low luminosity, the detection rates of masers would be significantly lower than the values obtained in higher-luminosity Class 0 protostars.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovery of 22 GHz Water Masers in the Serpens South Region

Using the Karl G. Jansky Very Large Array (VLA), we have conducted a survey for 22 GHz, 6{sub 1,6}–5{sub 2,3} H{sub 2}O masers toward the Serpens South region. The masers were also observed with the Very Long Baseline Array following the VLA detections. We detect for the first time H{sub 2}O masers in the Serpens South region that are found to be associated to three Class 0–Class I objects, including the two brightest protostars in the Serpens South cluster, known as CARMA-6 and CARMA-7. We also detect H{sub 2}O masers associated to a source with no outflow or jet features. We suggest that this source is most probably a background asymptotic giant branch star projected in the direction of Serpens South. The spatial distribution of the emission spots suggest that the masers in the three Class 0–Class I objects emerge very close to the protostars and are likely excited in shocks driven by the interaction between a protostellar jet and the circumstellar material. Based on the comparison of the distributions of bolometric luminosity of sources hosting 22 GHz H{sub 2}O masers and 162 young stellar objects covered by our observations, we identify a limit of L {sub Bol} ≈ 10L {sub ⊙} for a source to host water masers. However, the maser emission shows strong variability in both intensity and velocity spread, and therefore masers associated to lower-luminosity sources may have been missed by our observations. We also report 11 new sources with radio continuum emission at 22 GHz.

47 OTHER INSTRUMENTATION↗

Cold CO Gas in the Envelopes of FU Orionis-type Young Eruptive Stars

FU Orionis-type objects (FUors) are young stellar objects experiencing large optical outbursts due to highly enhanced accretion from the circumstellar disk onto the star. FUors are often surrounded by massive envelopes, which play a significant role in the outburst mechanism. Conversely, the subsequent eruptions might gradually clear up the obscuring envelope material and drive the protostar on its way to become a disk-only T Tauri star. Here we present an APEX {sup 12}CO and {sup 13}CO survey of eight southern and equatorial FUors. We measure the mass of the gaseous material surrounding our targets, locate the source of the CO emission, and derive physical parameters for the envelopes and outflows, where detected. Our results support the evolutionary scenario where FUors represent a transition phase from envelope-surrounded protostars to classical T Tauri stars.

79 ASTRONOMY AND ASTROPHYSICS↗

L483: Warm Carbon-chain Chemistry Source Harboring Hot Corino Activity

The Class 0 protostar, L483, has been observed in various molecular lines in the 1.2 mm band at a subarcsecond resolution with ALMA. An infalling–rotating envelope is traced by the CS line, while a very compact component with a broad velocity width is observed for the CS, SO, HNCO, NH{sub 2}CHO, and HCOOCH{sub 3} lines. Although this source is regarded as the warm carbon-chain chemistry (WCCC) candidate source at a 1000 au scale, complex organic molecules characteristic of hot corinos such as NH{sub 2}CHO and HCOOCH{sub 3} are detected in the vicinity of the protostar. Thus, both hot corino chemistry and WCCC are seen in L483. Although such a mixed chemical character source has been recognized as an intermediate source in previous single-dish observations, we here report the first spatially resolved detection. A kinematic structure of the infalling–rotating envelope is roughly explained by a simple ballistic model with a protostellar mass of 0.1–0.2 M{sub ⊙} and a radius of the centrifugal barrier (half of the centrifugal radius) of 30–200 au, assuming an inclination angle of 80° (0° for face-on). The broad-line emission observed in the above molecules most likely comes from the disk component inside the centrifugal barrier. Thus, a drastic chemical change is seen around the centrifugal barrier.

79 ASTRONOMY AND ASTROPHYSICS↗

The JCMT Transient Survey: Data Reduction and Calibration Methods

Though there has been a significant amount of work investigating the early stages of low-mass star formation in recent years, the evolution of the mass assembly rate onto the central protostar remains largely unconstrained. Examining in depth the variation in this rate is critical to understanding the physics of star formation. Instabilities in the outer and inner circumstellar disk can lead to episodic outbursts. Observing these brightness variations at infrared or submillimeter wavelengths constrains the current accretion models. The JCMT Transient Survey is a three-year project dedicated to studying the continuum variability of deeply embedded protostars in eight nearby star-forming regions at a one-month cadence. We use the SCUBA-2 instrument to simultaneously observe these regions at wavelengths of 450 and 850 μ m. In this paper, we present the data reduction techniques, image alignment procedures, and relative flux calibration methods for 850 μ m data. We compare the properties and locations of bright, compact emission sources fitted with Gaussians over time. Doing so, we achieve a spatial alignment of better than 1″ between the repeated observations and an uncertainty of 2%–3% in the relative peak brightness of significant, localized emission. This combination of imaging performance is unprecedented in ground-based, single-dish submillimeter observations. Finally, we identify a few sources that show possible and confirmed brightness variations. These sources will be closely monitored and presented in further detail in additional studies throughout the duration of the survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Zoom-in Simulations of Protoplanetary Disks Starting from GMC Scales

We investigate the formation of protoplanetary disks around nine solar-mass stars formed in the context of a (40 pc){sup 3} Giant Molecular Cloud model, using ramses adaptive mesh refinement simulations extending over a scale range of about 4 million, from an outer scale of 40 pc down to cell sizes of 2 au. Our most important result is that the accretion process is heterogeneous in multiple ways: in time, in space, and among protostars of otherwise similar mass. Accretion is heterogeneous in time, in the sense that accretion rates vary during the evolution, with generally decreasing profiles, whose slopes vary over a wide range, and where accretion can increase again if a protostar enters a region with increased density and low speed. Accretion is heterogeneous in space, because of the mass distribution, with mass approaching the accreting star–disk system in filaments and sheets. Finally, accretion is heterogeneous among stars, since the detailed conditions and dynamics in the neighborhood of each star can vary widely. We also investigate the sensitivity of disk formation to physical conditions and test their robustness by varying numerical parameters. We find that disk formation is robust even when choosing the least favorable sink particle parameters, and that turbulence cascading from larger scales is a decisive factor in disk formation. We also investigate the transport of angular momentum, finding that the net inward mechanical transport is compensated for mainly by an outward-directed magnetic transport, with a contribution from gravitational torques usually subordinate to the magnetic transport.

79 ASTRONOMY AND ASTROPHYSICS↗

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗