skip to main content


Title: OMC-1 dust polarization in ALMA Band 7: diagnosing grain alignment mechanisms in the vicinity of Orion Source I
ABSTRACT We present ALMA Band 7 polarization observations of the OMC-1 region of the Orion molecular cloud. We find that the polarization pattern observed in the region is likely to have been significantly altered by the radiation field of the >104 L⊙ high-mass protostar Orion Source I. In the protostar’s optically thick disc, polarization is likely to arise from dust self-scattering. In material to the south of Source I – previously identified as a region of ‘anomalous’ polarization emission – we observe a polarization geometry concentric around Source I. We demonstrate that Source I’s extreme luminosity may be sufficient to make the radiative precession time-scale shorter than the Larmor time-scale for moderately large grains ($\gt 0.005\!-\!0.1\, \mu$m), causing them to precess around the radiation anisotropy vector (k-RATs) rather than the magnetic field direction (B-RATs). This requires relatively unobscured emission from Source I, supporting the hypothesis that emission in this region arises from the cavity wall of the Source I outflow. This is one of the first times that evidence for k-RAT alignment has been found outside of a protostellar disc or AGB star envelope. Alternatively, the grains may remain aligned by B-RATs and trace gas infall on to the Main Ridge. Elsewhere, we largely find the magnetic field geometry to be radial around the BN/KL explosion centre, consistent with previous observations. However, in the Main Ridge, the magnetic field geometry appears to remain consistent with the larger-scale magnetic field, perhaps indicative of the ability of the dense Ridge to resist disruption by the BN/KL explosion.  more » « less
Award ID(s):
1715867
NSF-PAR ID:
10388653
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; « less
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
503
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3414 to 3433
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Resonant tunneling diodes (RTDs) have come full-circle in the past 10 years after their demonstration in the early 1990s as the fastest room-temperature semiconductor oscillator, displaying experimental results up to 712 GHz and fmax values exceeding 1.0 THz [1]. Now the RTD is once again the preeminent electronic oscillator above 1.0 THz and is being implemented as a coherent source [2] and a self-oscillating mixer [3], amongst other applications. This paper concerns RTD electroluminescence – an effect that has been studied very little in the past 30+ years of RTD development, and not at room temperature. We present experiments and modeling of an n-type In0.53Ga0.47As/AlAs double-barrier RTD operating as a cross-gap light emitter at ~300K. The MBE-growth stack is shown in Fig. 1(a). A 15-μm-diam-mesa device was defined by standard planar processing including a top annular ohmic contact with a 5-μm-diam pinhole in the center to couple out enough of the internal emission for accurate free-space power measurements [4]. The emission spectra have the behavior displayed in Fig. 1(b), parameterized by bias voltage (VB). The long wavelength emission edge is at  = 1684 nm - close to the In0.53Ga0.47As bandgap energy of Ug ≈ 0.75 eV at 300 K. The spectral peaks for VB = 2.8 and 3.0 V both occur around  = 1550 nm (h = 0.75 eV), so blue-shifted relative to the peak of the “ideal”, bulk InGaAs emission spectrum shown in Fig. 1(b) [5]. These results are consistent with the model displayed in Fig. 1(c), whereby the broad emission peak is attributed to the radiative recombination between electrons accumulated on the emitter side, and holes generated on the emitter side by interband tunneling with current density Jinter. The blue-shifted main peak is attributed to the quantum-size effect on the emitter side, which creates a radiative recombination rate RN,2 comparable to the band-edge cross-gap rate RN,1. Further support for this model is provided by the shorter wavelength and weaker emission peak shown in Fig. 1(b) around = 1148 nm. Our quantum mechanical calculations attribute this to radiative recombination RR,3 in the RTD quantum well between the electron ground-state level E1,e, and the hole level E1,h. To further test the model and estimate quantum efficiencies, we conducted optical power measurements using a large-area Ge photodiode located ≈3 mm away from the RTD pinhole, and having spectral response between 800 and 1800 nm with a peak responsivity of ≈0.85 A/W at  =1550 nm. Simultaneous I-V and L-V plots were obtained and are plotted in Fig. 2(a) with positive bias on the top contact (emitter on the bottom). The I-V curve displays a pronounced NDR region having a current peak-to-valley current ratio of 10.7 (typical for In0.53Ga0.47As RTDs). The external quantum efficiency (EQE) was calculated from EQE = e∙IP/(∙IE∙h) where IP is the photodiode dc current and IE the RTD current. The plot of EQE is shown in Fig. 2(b) where we see a very rapid rise with VB, but a maximum value (at VB= 3.0 V) of only ≈2×10-5. To extract the internal quantum efficiency (IQE), we use the expression EQE= c ∙i ∙r ≡ c∙IQE where ci, and r are the optical-coupling, electrical-injection, and radiative recombination efficiencies, respectively [6]. Our separate optical calculations yield c≈3.4×10-4 (limited primarily by the small pinhole) from which we obtain the curve of IQE plotted in Fig. 2(b) (right-hand scale). The maximum value of IQE (again at VB = 3.0 V) is 6.0%. From the implicit definition of IQE in terms of i and r given above, and the fact that the recombination efficiency in In0.53Ga0.47As is likely limited by Auger scattering, this result for IQE suggests that i might be significantly high. To estimate i, we have used the experimental total current of Fig. 2(a), the Kane two-band model of interband tunneling [7] computed in conjunction with a solution to Poisson’s equation across the entire structure, and a rate-equation model of Auger recombination on the emitter side [6] assuming a free-electron density of 2×1018 cm3. We focus on the high-bias regime above VB = 2.5 V of Fig. 2(a) where most of the interband tunneling should occur in the depletion region on the collector side [Jinter,2 in Fig. 1(c)]. And because of the high-quality of the InGaAs/AlAs heterostructure (very few traps or deep levels), most of the holes should reach the emitter side by some combination of drift, diffusion, and tunneling through the valence-band double barriers (Type-I offset) between InGaAs and AlAs. The computed interband current density Jinter is shown in Fig. 3(a) along with the total current density Jtot. At the maximum Jinter (at VB=3.0 V) of 7.4×102 A/cm2, we get i = Jinter/Jtot = 0.18, which is surprisingly high considering there is no p-type doping in the device. When combined with the Auger-limited r of 0.41 and c ≈ 3.4×10-4, we find a model value of IQE = 7.4% in good agreement with experiment. This leads to the model values for EQE plotted in Fig. 2(b) - also in good agreement with experiment. Finally, we address the high Jinter and consider a possible universal nature of the light-emission mechanism. Fig. 3(b) shows the tunneling probability T according to the Kane two-band model in the three materials, In0.53Ga0.47As, GaAs, and GaN, following our observation of a similar electroluminescence mechanism in GaN/AlN RTDs (due to strong polarization field of wurtzite structures) [8]. The expression is Tinter = (2/9)∙exp[(-2 ∙Ug 2 ∙me)/(2h∙P∙E)], where Ug is the bandgap energy, P is the valence-to-conduction-band momentum matrix element, and E is the electric field. Values for the highest calculated internal E fields for the InGaAs and GaN are also shown, indicating that Tinter in those structures approaches values of ~10-5. As shown, a GaAs RTD would require an internal field of ~6×105 V/cm, which is rarely realized in standard GaAs RTDs, perhaps explaining why there have been few if any reports of room-temperature electroluminescence in the GaAs devices. [1] E.R. Brown,et al., Appl. Phys. Lett., vol. 58, 2291, 1991. [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [2] M. Feiginov et al., Appl. Phys. Lett., 99, 233506, 2011. [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [3] Y. Nishida et al., Nature Sci. Reports, 9, 18125, 2019. [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [4] P. Fakhimi, et al., 2019 DRC Conference Digest. [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). 
    more » « less
  2. Context. Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims. We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods. We present new observations of the dust polarized emission at 850 μ m toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with the James Clerk Maxwell Telescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity ( PI ), the polarization fraction ( PF ), and the plane-of-the-sky B-field angle ( χ B_POS ) toward the whole region, along the 10 pc-long ridge and along the sub-filaments connected to the ridge and the hubs. We derived the power spectra of the intensity and χ B POS along the ridge crest and compared them with the results obtained from simulated filaments. Results. The observations span ~3 orders of magnitude in Stokes I and PI and ~2 orders of magnitude in PF (from ~0.2 to ~ 20%). A large scatter in PI and PF is observed for a given value of I . Our analyses show a complex B-field structure when observed over the whole region (~ 10 pc); however, at smaller scales (~1 pc), χ B POS varies coherently along the crests of the filament network. The observed power spectrum of χ B POS can be well represented with a power law function with a slope of − 1.33 ± 0.23, which is ~20% shallower than that of I . We find that this result is compatible with the properties of simulated filaments and may indicate the physical processes at play in the formation and evolution of star-forming filaments. Along the sub-filaments, χ B POS rotates frombeing mostly perpendicular or randomly oriented with respect to the crests to mostly parallel as the sub-filaments merge with the ridge and hubs. This variation of the B-field structure along the sub-filaments may be tracing local velocity flows of infalling matter in the ridge and hubs. Our analysis also suggests a variation in the energy balance along the crests of these sub-filaments, from magnetically critical or supercritical at their far ends to magnetically subcritical near the ridge and hubs. We also detect an increase in PF toward the high-column density ( N H 2 ≳ 10 23  cm −2 ) star cluster-forming hubs. These latter large PF values may be explained by the increase in grain alignment efficiency due to stellar radiation from the newborn stars, combined with an ordered B-field structure. Conclusions. These observational results reveal for the first time the characteristics of the small-scale (down to ~ 0.1 pc) B-field structure of a 10 pc-long hub-filament system. Our analyses show variations in the polarization properties along the sub-filaments that may be tracing the evolution of their physical properties during their interaction with the ridge and hubs. We also detect an impact of feedback from young high-mass stars on the local B-field structure and the polarization properties, which could put constraints on possible models for dust grain alignment and provide important hints as to the interplay between the star formation activity and interstellar B-fields. 
    more » « less
  3. Abstract Located in the Large Magellanic Cloud and mostly irradiated by the massive star cluster R136, 30 Doradus is an ideal target to test the leading theory of grain alignment and rotational disruption by RAdiative Torques (RATs). Here, we use publicly available polarized thermal dust emission observations of 30 Doradus at 89, 154, and 214 μ m using SOFIA/HAWC+. We analyze the variation of the dust polarization degree ( p ) with the total emission intensity ( I ), the dust temperature ( T d ), and the gas column density ( N H ) constructed from Herschel data. The 30 Doradus complex is divided into two main regions relative to R136, namely North and South. In the North, we find that the polarization degree first decreases and then increases before decreasing again when the dust temperature increases toward the irradiating cluster R136. The first depolarization likely arises from the decrease in grain alignment efficiency toward the dense medium due to the attenuation of the interstellar radiation field and the increase in the gas density. The second trend (the increase of p with T d ) is consistent with the RAT alignment theory. The final trend (the decrease of p with T d ) is consistent with the RAT alignment theory only when the grain rotational disruption by RATs is taken into account. In the South, we find that the polarization degree is nearly independent of the dust temperature, while the grain alignment efficiency is higher around the peak of the gas column density and decreases toward the radiation source. The latter feature is also consistent with the prediction of rotational disruption by RATs. 
    more » « less
  4. Abstract

    The Orion Kleinmann-Low nebula (Orion KL) is notoriously complex and exhibits a range of physical and chemical components. We conducted high-angular-resolution (subarcsecond) observations of13CH3OHν= 0 (∼0.″3 and ∼0.″7) and CH3CNν8= 1 (∼0.″2 and ∼0.″9) line emission with the Atacama Large Millimeter/submillimeter Array (ALMA) to investigate Orion KL’s structure on small spatial scales (≤350 au). Gas kinematics, excitation temperatures, and column densities were derived from the molecular emission via a pixel-by-pixel spectral line fitting of the image cubes, enabling us to examine the small-scale variation of these parameters. Subregions of the Hot Core have a higher excitation temperature in a 0.″2 beam than in a 0.″9 beam, indicative of possible internal sources of heating. Furthermore, the velocity field includes a bipolar ∼7–8 km s−1feature with a southeast–northwest orientation against the surrounding ∼4–5 km s−1velocity field, which may be due to an outflow. We also find evidence of a possible source of internal heating toward the Northwest Clump, since the excitation temperature there is higher in a smaller beam versus a larger beam. Finally, the region southwest of the Hot Core (Hot Core-SW) presents itself as a particularly heterogeneous region bridging the Hot Core and Compact Ridge. Additional studies to identify the (hidden) sources of luminosity and heating within Orion KL are necessary to better understand the nebula and its chemistry.

     
    more » « less
  5. ABSTRACT Polarized continuum emission from aligned grains in discs around young stellar objects can be used to probe the magnetic field, radiation anisotropy, or drift between dust and gas, depending on whether the non-spherical grains are aligned magnetically, radiatively, or mechanically. We show that it can also be used to probe another key disc property – the temperature gradient – along sightlines that are optically thick, independent of the grain alignment mechanism. We first illustrate the technique analytically using a simple 1D slab model, which yields an approximate formula that relates the polarization fraction to the temperature gradient with respect to the optical depth τ at the τ = 1 surface. The formula is then validated using models of stellar irradiated discs with and without accretion heating. The promises and challenges of the technique are illustrated with a number of Class 0 and I discs with ALMA dust polarization data, including NGC 1333 IRAS4A1, IRAS 16293B, BHB 07-11, L1527, HH 212, and HH 111. We find, in particular, that the sightlines passing through the near-side of a highly inclined disc trace different temperature gradient directions than those through the far-side, which can lead to a polarization orientation on the near-side that is orthogonal to that on the far-side, and that the HH 111 disc may be such a case. Our technique for probing the disc temperature gradient through dust polarization can complement other methods, particularly those using molecular lines. 
    more » « less