Source author record

Ian W. Stephens

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

18 published item(s)

preprint2022arXiv

The Magnetic Field in the Milky Way Filamentary Bone G47

Star formation primarily occurs in filaments where magnetic fields are expected to be dynamically important. The largest and densest filaments trace spiral structure within galaxies. Over a dozen of these dense ($\sim$10$^4$\,cm$^{-3}$) and long ($>$10\,pc) filaments have been found within the Milky Way, and they are often referred to as "bones." Until now, none of these bones have had their magnetic field resolved and mapped in their entirety. We introduce the SOFIA legacy project FIELDMAPS which has begun mapping $\sim$10 of these Milky Way bones using the HAWC+ instrument at 214\,$μ$m and 18$\farcs$2 resolution. Here we present a first result from this survey on the $\sim$60\,pc long bone G47. Contrary to some studies of dense filaments in the Galactic plane, we find that the magnetic field is often not perpendicular to the spine (i.e., the center-line of the bone). Fields tend to be perpendicular in the densest areas of active star formation and more parallel or random in other areas. The average field is neither parallel or perpendicular to the Galactic plane nor the bone. The magnetic field strengths along the spine typically vary from $\sim$20 to $\sim$100\,$μ$G. Magnetic fields tend to be strong enough to suppress collapse along much of the bone, but for areas that are most active in star formation, the fields are notably less able to resist gravitational collapse.

preprint2022arXiv

The Twisted Magnetic Field of the Protobinary L483

We present H-band (1.65 $μ$m) and SOFIA HAWC+ 154 $μ$m polarization observations of the low-mass core L483. Our H-band observations reveal a magnetic field that is overwhelmingly in the E-W direction, which is approximately parallel to the bipolar outflow that is observed in scattered IR light and in single-dish $^{12}$CO observations. From our 154 $μ$m data, we infer a $\sim$ 45$^{\circ}$ twist in the magnetic field within the inner 5" (1000 au) of L483. We compare these new observations with published single-dish 350 $μ$m polarimetry and find that the 10,000 au scale H-band data match the smaller scale 350 $μ$m data, indicating that the collapse of L483 is magnetically regulated on these larger scales. We also present high-resolution 1.3 mm ALMA data of L483 which reveals it is a close binary star with a separation of 34 au. The plane of the binary of L483 is observed to be approximately parallel to the twisted field in the inner 1000 au. Comparing this result to the $\sim$ 1000 au protostellar envelope, we find that the envelope is roughly perpendicular to the 1000 au HAWC+ field. Using the data presented, we speculate that L483 initially formed as a wide binary and the companion star migrated to its current position, causing an extreme shift in angular momentum thereby producing the twisted magnetic field morphology observed. More observations are needed to further test this scenario.

preprint2021arXiv

Evolution and Kinematics of Protostellar Envelopes in the Perseus Molecular Cloud

We present a comprehensive analysis on the evolution of envelopes surrounding protostellar systems in the Perseus molecular cloud using data from the MASSES survey. We focus our attention to the C$^{18}$O(2--1) spectral line, and we characterize the shape, size, and orientation of 54 envelopes and measure their fluxes, velocity gradients, and line widths. To look for evolutionary trends, we compare these parameters to the bolometric temperature Tbol, a tracer of protostellar age. We find evidence that the angular difference between the elongation angle of the C$^{18}$O envelope and the outflow axis direction generally becomes increasingly perpendicular with increasing Tbol, suggesting the envelope evolution is directly affected by the outflow evolution. We show that this angular difference changes at Tbol = $53 \pm 20$ K, which includes the conventional delineation between Class 0 and I protostars of 70K. We compare the C$^{18}$O envelopes with larger gaseous structures in other molecular clouds and show that the velocity gradient increases with decreasing radius ($|\mathcal{G}| \sim R^{-0.72 \pm 0.06}$). From the velocity gradients we show that the specific angular momentum follows a power law fit $J/M \propto R^{1.83 \pm 0.05}$ for scales from 1pc down to $\sim$500 au, and we cannot rule out a possible flattening out at radii smaller than $\sim$1000 au.

preprint2021arXiv

Magnetic Fields in Massive Star-Forming Regions (MagMaR) II. Tomography Through Dust and Molecular Line Polarization in NGC 6334I(N)

Here, we report ALMA detections of polarized emission from dust, CS($J=5 \rightarrow 4$), and C$^{33}$S($J=5 \rightarrow 4$) toward the high-mass star-forming region NGC6334I(N). A clear ``hourglass'' magnetic field morphology was inferred from the polarized dust emission which is also directly seen from the polarized CS emission across velocity, where the polarization appears to be parallel to the field. By considering previous findings, the field retains a pinched shape which can be traced to clump length-scales from the envelope scales traced by ALMA, suggesting that the field is dynamically important across multiple length-scales in this region. The CS total intensity emission is found to be optically thick ($τ_{\mathrm{CS}} = 32 \pm 12$) while the C$^{33}$S emission appears to be optically thin ($τ_{\mathrm{C^{33}S}} = 0.1 \pm 0.01$). This suggests that sources of anisotropy other than large velocity gradients, i.e. anisotropies in the radiation field are required to explain the polarized emission from CS seen by ALMA. By using four variants of the Davis-Chandrasekhar-Fermi technique and the angle dispersion function methods (ADF), we obtain an average of estimates for the magnetic field strength onto the plane of the sky of $\left< \mathrm{B}_{\mathrm{pos}} \right> = 16$ mG from the dust and $\left< \mathrm{B}_{\mathrm{pos}} \right> \sim 2$ mG from the CS emission, where each emission traces different molecular hydrogen number densities. This effectively enables a tomographic view of the magnetic field within a single ALMA observation.

preprint2020arXiv

Magnetic Field Structure in Spheroidal Star-Forming Clouds. II. Estimating Field Structure from Observed Maps

This paper presents models to estimate the structure of density and magnetic field strength in spheroidal condensations, from maps of their column density and their polarization of magnetically aligned dust grains. The density model is obtained by fitting a column density map with an embedded p = 2 Plummer spheroid of any aspect ratio and inclination. The magnetic properties are based on the density model, on the Davis-Chandrasekhar-Fermi (DCF) model of Alfvénic fluctuations, and on the Spheroid Flux Freezing (SFF) model of mass and flux conservation in Paper I. The field strength model has the resolution of the column density map, which is finer than the resolution of the DCF estimate of field strength. The models are applied to ALMA observations of the envelope of the protostar BHR71 IRS1. Column density fits give the density model, from (2.0 +- 0.4) x 10^5 cm^-3 to (7 +- 1) x 10^7 cm^-3 . The density model predicts the field directions map, which fits the polarization map best within 1100 au, with standard deviation of angle differences 17°. In this region the DCF mean field strength is 0.7 +- 0.2 mG and the envelope mass is supercritical, with ratio of mass to magnetic critical mass 1.5 +- 0.4. The SFF field strength profile scales with the DCF field strength, from 60 x 10μG to 4+-1 mG. The spatial resolution of the SFF field strength estimate is finer than the DCF resolution by a factor ~7, and the peak SFF field strength exceeds the DCF field strength by a factor ~5.

preprint2019arXiv

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $μ$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

preprint2019arXiv

The Interaction Between the Supernova Remnant W41 and the Filamentary Infrared Dark Cloud G23.33-0.30

G23.33-0.30 is a 600 $M_{\odot}$ infrared dark molecular filament that exhibits large NH$_3$ velocity dispersions ($σ\sim 8 \ \rm{km \ s^{-1}}$) and bright, narrow NH$_3$(3,3) line emission. We have probed G23.33-0.30 at the $<0.1$ pc scale and confirmed that the narrow NH$_3$(3,3) line is emitted by four rare NH$_3$(3,3) masers, which are excited by a large-scale shock impacting the filament. G23.33-0.30 also displays a velocity gradient along its length, a velocity discontinuity across its width, shock-tracing SiO(5-4) emission extended throughout the filament, broad turbulent line widths in NH$_3$(1,1) through (6,6), CS(5-4), and SiO(5-4), as well as an increased NH$_3$ rotational temperature ($T_{\rm{rot}}$) and velocity dispersion ($σ$) associated with the shocked, blueshifted component. The correlations among $T_{\rm{rot}}$, $σ$, and $V_{\rm{LSR}}$ implies that the shock is accelerating, heating, and adding turbulent energy to the filament gas. Given G23.33-0.30's location within the giant molecular cloud G23.0-0.4, we speculate that the shock and NH$_3$(3,3) masers originated from the supernova remnant W41, which exhibits additional evidence of an interaction with G23.0-0.4. We have also detected the 1.3 mm dust continuum emission from at least three embedded molecular cores associated with G23.33-0.30. Although the cores have moderate gas masses ($M = 7-10$ M$_{\odot}$), their large virial parameters ($α=4-9$) suggest that they will not collapse to form stars. The turbulent line widths of the cores may indicate negative feedback due to the SNR shock.

preprint2016arXiv

Disk Polarization From Both Emission and Scattering of Magnetically Aligned Grains: The Case of NGC 1333 IRAS4A1

Dust polarization in millimeter (and centimeter) has been mapped in disks around an increasing number of young stellar objects. It is usually thought to come from emission by magnetically aligned (non-spherical) grains, but can also be produced by dust scattering. We present a semi-analytic theory of disk polarization that includes both the direction emission and scattering, with an emphasis on their relative importance and how they are affected by the disk inclination. For face-on disks, both emission and scattering tend to produce polarization in the radial direction, making them difficult to distinguish, although the scattering-induced polarization can switch to the azimuthal direction if the incident radiation is beamed strongly enough in the radial direction in the disk plane. Disk inclination affects the polarizations from emission and scattering differently, especially on the major axis where, in the edge-on limit, the former vanishes while the latter reaches a polarization fraction as large as 1/3. The polarizations from the two competing mechanisms tend to cancel each other on the major axis, producing two low polarization "holes" (one on each side of the center) under certain conditions. We find tantalizing evidence for at least one such "hole" in NGC1333 IRAS4A1, whose polarization observed at 8 mm on the 100 AU scale is indicative of a pattern dominated by scattering close to the center and by direction emission in the outer region. If true, it would imply not only that a magnetic field exists on the disk scale, but that it is strong enough to align large, possibly mm-sized, grains.

preprint2016arXiv

Linking Dense Gas from the Milky Way to External Galaxies

In a survey of 65 galaxies, Gao & Solomon (2004) found a tight linear relation between the infrared luminosity ($L_{\rm{IR}}$, a proxy for the star formation rate) and the HCN(1-0) luminosity ($L_{\rm{HCN}}$). Wu et al. (2005, 2010) found that this relation extends from these galaxies to the much less luminous Galactic molecular high-mass star-forming clumps ($\sim$1 pc scales), and posited that there exists a characteristic ratio $L_{\rm{IR}}$ /$L_{\rm{HCN}}$ for high-mass star-forming clumps. The Gao-Solomon relation for galaxies could then be explained as a summation of large numbers of high-mass star-forming clumps, resulting in the same $L_{\rm{IR}}$ /$L_{\rm{HCN}}$ ratio for galaxies. We test this explanation and other possible origins of the Gao-Solomon relation using high-density tracers (including HCN(1-0), N$_2$H$^+$(1-0), HCO$^+$(1-0), HNC(1-0), HC$_3$N(10-9), and C$_2$H(1-0)) for $\sim$300 Galactic clumps from the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey. The MALT90 data show that the Gao-Solomon relation in galaxies cannot be satisfactorily explained by the blending of large numbers of high-mass clumps in the telescope beam. Not only do the clumps have a large scatter in the $L_{\rm{IR}}$/$L_{\rm{HCN}}$ ratio, but also far too many high-mass clumps are required to account for the Galactic IR and HCN luminosities. We suggest that the scatter in the $L_{\rm{IR}}$/$L_{\rm{HCN}}$ ratio converges to the scatter of the Gao-Solomon relation at some size-scale $\gtrsim$1 kpc. We suggest that the Gao-Solomon relation could instead result from of a universal large-scale star formation efficiency, initial mass function, core mass function, and clump mass function.

preprint2015arXiv

Interferometric Observations of High-Mass Star-Forming Clumps with Unusual N2H+/HCO+ Line Ratios

The Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey has detected high-mass star-forming clumps with anomalous N$_2$H$^+$/HCO$^+$(1-0) integrated intensity ratios that are either unusually high ("N$_2$H$^+$ rich") or unusually low ("N$_2$H$^+$ poor"). With 3 mm observations from the Australia Telescope Compact Array (ATCA), we imaged two N$_2$H$^+$ rich clumps, G333.234-00.061 and G345.144-00.216, and two N$_2$H$^+$ poor clumps, G351.409+00.567 and G353.229+00.672. In these clumps, the N$_2$H$^+$ rich anomalies arise from extreme self-absorption of the HCO$^+$ line. G333.234-00.061 contains two of the most massive protostellar cores known with diameters of less than 0.1 pc, separated by a projected distance of only 0.12 pc. Unexpectedly, the higher mass core appears to be at an earlier evolutionary stage than the lower mass core, which may suggest that two different epochs of high-mass star formation can occur in close proximity. Through careful analysis of the ATCA observations and MALT90 clumps (including the G333, NGC 6334, and NGC 6357 star formation regions), we find that N$_2$H$^+$ poor anomalies arise at clump-scales and are caused by lower relative abundances of N$_2$H$^+$ due to the distinct chemistry of H II regions or photodissociation regions.

preprint2015arXiv

Kinematics of the Envelope and Two Bipolar Jets in the Class 0 Protostellar System L1157

A massive envelope and a strong bipolar outflow are the two main structures characterizing the youngest protostellar systems. In order to understand the physical properties of a bipolar outflow and the relationship with those of the envelope, we obtained a mosaic map covering the whole bipolar outflow of the youngest protostellar system L1157 with about $5"$ angular resolution in CO J=2-1 using the Combined Array for Research in Millimeter-wave Astronomy. By utilizing these observations of the whole bipolar outflow, we estimate its physical properties and show that they are consistent with multiple jets. We also constrain a preferred precession direction. In addition, we observed the central envelope structure with $2"$ resolution in the $λ=1.3$ and 3 mm continua and various molecular lines: C$^{17}$O, C$^{18}$O, $^{13}$CO, CS, CN, N$_2$H$^+$, CH$_3$OH, H$_2$O, SO, and SO$_2$. All the CO isotopes and CS, CN, and N$_2$H$^+$ have been detected and imaged. We marginally detected the features that can be interpreted as a rotating inner envelope in C$^{17}$O and C$^{18}$O and as an infalling outer envelope in N$_2$H$^+$. We also estimated the envelope and central protostellar masses and found that the dust opacity spectral index changes with radius.

preprint2014arXiv

Spatially Resolved Magnetic Field Structure in the Disk of a T Tauri Star

Magnetic fields in accretion disks play a dominant role during the star formation process but have hitherto been observationally poorly constrained. Field strengths have been inferred on T Tauri stars themselves and possibly in the innermost part of the accretion disk, but the strength and morphology of the field in the bulk of the disk have not been observed. Unresolved measurements of polarized emission (arising from elongated dust grains aligned perpendicular to the field) imply average fields aligned with the disks. Theoretically, the fields are expected to be largely toroidal, poloidal, or a mixture of the two, which imply different mechanisms for transporting angular momentum in the disks of actively accreting young stars such as HL Tau. Here we report resolved measurements of the polarized 1.25 mm continuum emission from HL Tau's disk. The magnetic field on a scale of 80 AU is coincident with the major axis (~210 AU diameter) of the disk. From this we conclude that the magnetic field inside the disk at this scale cannot be dominated by a vertical component, though a purely toroidal field does not fit the data well either. The unexpected morphology suggests that the magnetic field's role for the accretion of a T Tauri star is more complex than the current theoretical understanding.

preprint2014arXiv

Spitzer Observations of Dust Emission from HII Regions in the Large Magellanic Cloud

Massive stars can alter physical conditions and properties of their ambient interstellar dust grains via radiative heating and shocks. The HII regions in the Large Magellanic Cloud (LMC) offer ideal sites to study the stellar energy feedback effects on dust because stars can be resolved, and the galaxy's nearly face-on orientation allows us to unambiguously associate HII regions with their ionizing massive stars. The Spitzer Space Telescope survey of the LMC provides multi-wavelength (3.6 to 160 $μ$m) photometric data of all HII regions. To investigate the evolution of dust properties around massive stars, we have analyzed spatially-resolved IR dust emission from two classical HII regions (N63 and N180) and two simple superbubbles (N70 and N144) in the LMC. We produce photometric spectral energy distributions (SEDs) of numerous small subregions for each region based on its stellar distributions and nebular morphologies. We use DustEM dust emission model fits to characterize the dust properties. Color-color diagrams and model fits are compared with the radiation field (estimated from photometric and spectroscopic surveys). Strong radial variations of SEDs can be seen throughout the regions, reflecting the available radiative heating. Emission from very small grains drastically increases at locations where the radiation field is the highest, while polycyclic aromatic hydrocarbons (PAHs) appear to be destroyed. PAH emission is the strongest in the presence of molecular clouds, provided that the radiation field is low.

preprint2014arXiv

TADPOL: A 1.3 mm Survey of Dust Polarization in Star-forming Cores and Regions

We present λ1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.

preprint2014arXiv

The Magnetic Field in the Class 0 Protostellar Disk of L1527

We present subarcsecond (~0.35") resolved observations of the 1.3 mm dust polarization from the edge-on circumstellar disk around the Class 0 protostar L1527. The inferred magnetic field is consistent with a dominantly toroidal morphology; there is no significantly detected vertical poloidal component to which observations of an edge-on disk are most sensitive. This suggests that angular momentum transport in Class 0 protostars (when large amounts of material are fed down to the disk from the envelope and accreted onto the protostar) is driven mainly by magnetorotational instability rather than magnetocentrifugal winds at 50 AU scales. In addition, with the data to date there is an early, tentative trend that R>30 AU disks have so far been found in Class 0 systems with average magnetic fields on the 1000 AU scale strongly misaligned with the rotation axis. The absence of such a disk in the aligned case could be due to efficient magnetic braking that disrupts disk formation. If this is the case, this implies that candidate Class 0 disk systems could be identified by the average magnetic field direction at ~1000 AU spatial scales.

preprint2013arXiv

Misalignment of Magnetic Fields and Outflows in Protostellar Cores

We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.

preprint2013arXiv

The Magnetic Field Morphology of the Class 0 Protostar L1157-mm

We present the first detection of polarization around the Class 0 low-mass protostar L1157-mm at two different wavelengths. We show polarimetric maps at large scales (10" resolution at 350 um) from the SHARC-II Polarimeter and at smaller scales (1.2"-4.5" at 1.3 mm) from the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The observations are consistent with each other and show inferred magnetic field lines aligned with the outflow. The CARMA observations suggest a full hourglass magnetic field morphology centered about the core; this is only the second well-defined hourglass detected around a low-mass protostar to date. We apply two different methods to CARMA polarimetric observations to estimate the plane-of-sky magnetic field magnitude, finding values of 1.4 and 3.4 mG.

preprint2010arXiv

The Galactic Magnetic Field's Effect in Star-Forming Region

We investigate the effect of the Milky Way's magnetic field in star forming regions using archived 350 micron polarization data on 52 Galactic star formation regions from the Hertz polarimeter module. The polarization angles and percentages for individual telescope beams were combined in order to produce a large-scale average for each source and for complexes of sources. In more than 80% of the sources, we find a meaningful mean magnetic field direction, implying the existence of an ordered magnetic field component at the scale of these sources. The average polarization angles were analyzed with respect to the Galactic coordinates in order to test for correlations between polarization percentage, polarization angle, intensity, and Galactic location. No correlation was found, which suggests that the magnetic field in dense molecular clouds is decoupled from the large-scale Galactic magnetic field. Finally, we show that the magnetic field directions in the complexes are consistent with a random distribution on the sky.