Thermal Dust Emission Research Articles

Dusty Gas Accretion onto Massive Black Holes and Infrared Diagnosis of the Eddington Ratio

Abstract Evidence for dust around supermassive black holes (SMBHs) in the early universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black holes (IMBHs) of in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the accretion of dusty gas onto IMBHs proceeds gently with small fluctuations of the accretion rate, whereas that of pristine gas causes more violent periodic bursts. At dust-to-gas mass ratios similar to the solar neighborhood, the time-averaged luminosity becomes smaller than that for primordial gas by one order of magnitude and the time-averaged Eddington ratio ranges from to in clouds with initial gas densities of . Our calculations show that the effect of dust opacity alone is secondary compared to the radiation pressure on dust in regulating the BH growth. We also derive spectral energy distributions at IR bands by calculating dust thermal emission and show that the flux ratio between and is closely related to the Eddington ratio. Thermal emission from hot dust near the BH dominates only during the phase of high accretion, producing higher flux density at . Therefore, we suggest that a combination of mid-IR observations by the James Webb Space Telescope and far-IR observations by ALMA or Spitzer can be used to estimate the Eddington ratio of massive BHs. We also extend our simple modeling to SMBHs of and show that ALMA can detect SMBHs of at .

ALMA Observations of Circumnuclear Disks in Early-type Galaxies: 12CO(2−1) and Continuum Properties

Abstract We present results from an Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 2 program to map CO(2−1) emission in nearby early-type galaxies (ETGs) that host circumnuclear gas disks. We obtained ∼0.″3 resolution Band 6 observations of seven ETGs selected on the basis of dust disks in Hubble Space Telescope images. We detect CO emission in five at high signal-to-noise ratio with the remaining two only faintly detected. All CO emission is coincident with the dust and is in dynamically cold rotation. Four ETGs show evidence of rapid central rotation; these are prime candidates for higher-resolution ALMA observations to measure the black hole masses. In this paper, we focus on the molecular gas and continuum properties. Total gas masses and H2 column densities for our five CO-bright galaxies are on average ∼108 M ☉ and cm−2 over the ∼kpc-scale disks, and analysis suggests that these disks are stabilized against gravitational fragmentation. The continuum emission of all seven galaxies is dominated by a central unresolved source, and in five we also detect a spatially extended component. The ∼230 GHz nuclear continua are modeled as power laws ranging from to within the observed frequency band. The extended continuum profiles of the two radio-bright (and CO-faint) galaxies are roughly aligned with their radio jet and suggest resolved synchrotron jets. The extended continua of the CO-bright disks are coincident with optically thick dust absorption and have spectral slopes that are consistent with thermal dust emission.

Observing Interstellar and Intergalactic Magnetic Fields

Observational results of interstellar and intergalactic magnetic fields are reviewed, including the fields in supernova remnants and loops, interstellar filaments and clouds, Hii regions and bubbles, the Milky Way and nearby galaxies, galaxy clusters, and the cosmic web. A variety of approaches are used to investigate these fields. The orientations of magnetic fields in interstellar filaments and molecular clouds are traced by polarized thermal dust emission and starlight polarization. The field strengths and directions along the line of sight in dense clouds and cores are measured by Zeeman splitting of emission or absorption lines. The large-scale magnetic fields in the Milky Way have been best probed by Faraday rotation measures of a large number of pulsars and extragalactic radio sources. The coherent Galactic magnetic fields are found to follow the spiral arms and have their direction reversals in arms and interarm regions in the disk. The azimuthal fields in the halo reverse their directions below and above the Galactic plane. The orientations of organized magnetic fields in nearby galaxies have been observed through polarized synchrotron emission. Magnetic fields in the intracluster medium have been indicated by diffuse radio halos, polarized radio relics, and Faraday rotations of embedded radio galaxies and background sources. Sparse evidence for very weak magnetic fields in the cosmic web is the detection of the faint radio bridge between the Coma cluster and A1367. Future observations should aim at the 3D tomography of the large-scale coherent magnetic fields in our Galaxy and nearby galaxies, a better description of intracluster field properties, and firm detections of intergalactic magnetic fields in the cosmic web.

Fourier-space combination ofPlanckandHerschelimages

Herschel has revolutionized our ability to measure column densities (N$_{\rm H}$) and temperatures (T) of molecular clouds thanks to its far infrared multiwavelength coverage. However, the lack of a well defined background intensity level in the Herschel data limits the accuracy of the N$_{\rm H}$ and T maps. We provide a method that corrects the missing Herschel background intensity levels using the Planck model for foreground Galactic thermal dust emission. We present a Fourier method that combines the publicly available Planck model on large angular scales with the Herschel images on smaller angular scales. We apply our method to two regions spanning a range of Galactic environments: Perseus and the Galactic plane region around $l = 11\deg$ (HiGal--11). We post-process the combined dust continuum emission images to generate column density and temperature maps. We compare these to previously adopted constant--offset corrections. We find significant differences ($\gtrsim$20\%) over significant ($\sim$15\%) areas of the maps, at low column densities ($N_{\rm H}\lesssim10^{22}$\,cm$^{-2}$) and relatively high temperatures ($T\gtrsim20$\,K). We also apply our method to synthetic observations of a simulated molecular cloud to validate our method. Our method successfully corrects the Herschel images, including both the constant--offset intensity level and the scale-dependent background variations measured by Planck. Our method improves the previous constant--offset corrections, which did not account for variations in the background emission levels.

The relation between the column density structures and the magnetic field orientation in the Vela C molecular complex

We statistically evaluated the relative orientation between gas column density structures, inferred from Herschel submillimetre observations, and the magnetic field projected on the plane of sky, inferred from polarized thermal emission of Galactic dust observed by the Balloon-borne Large-Aperture Submillimetre Telescope for Polarimetry (BLASTPol) at 250, 350, and 500 μm, towards the Vela C molecular complex. First, we find very good agreement between the polarization orientations in the three wavelength-bands, suggesting that, at the considered common angular resolution of 3.́0 that corresponds to a physical scale of approximately 0.61 pc, the inferred magnetic field orientation is not significantly affected by temperature or dust grain alignment effects. Second, we find that the relative orientation between gas column density structures and the magnetic field changes progressively with increasing gas column density, from mostly parallel or having no preferred orientation at low column densities to mostly perpendicular at the highest column densities. This observation is in agreement with previous studies by the Planck collaboration towards more nearby molecular clouds. Finally, we find a correspondencebetween (a) the trends in relative orientation between the column density structures and the projected magnetic field; and (b) the shape of the column density probability distribution functions (PDFs). In the sub-regions of Vela C dominated by one clear filamentary structure, or “ridges”, where the high-column density tails of the PDFs are flatter, we find a sharp transition from preferentially parallel or having no preferred relative orientation at low column densities to preferentially perpendicular at highest column densities. In the sub-regions of Vela C dominated by several filamentary structures with multiple orientations, or “nests”, where the maximum values of the column density are smaller than in the ridge-like sub-regions and the high-column density tails of the PDFs are steeper, such a transition is also present, but it is clearly less sharp than in the ridge-like sub-regions. Both of these results suggest that the magnetic field is dynamically important for the formation of density structures in this region.

The Coldest Place in the Universe: Probing the Ultra-cold Outflow and Dusty Disk in the Boomerang Nebula

Abstract Our Cycle 0 ALMA observations confirmed that the Boomerang Nebula is the coldest known object in the universe, with a massive high-speed outflow that has cooled significantly below the cosmic background temperature. Our new CO 1–0 data reveal heretofore unseen distant regions of this ultra-cold outflow, out to ≳120,000 au. We find that in the ultra-cold outflow, the mass-loss rate ( ) increases with radius, similar to its expansion velocity (V)—taking , we find . The mass in the ultra-cold outflow is M ⊙, and the Boomerang’s main-sequence progenitor mass is M ⊙. Our high angular resolution ( ) CO J = 3–2 map shows the inner bipolar nebula’s precise, highly collimated shape, and a dense central waist of size (FWHM) ∼1740 au × 275 au. The molecular gas and the dust as seen in scattered light via optical Hubble Space Telescope imaging show a detailed correspondence. The waist shows a compact core in thermal dust emission at 0.87–3.3 mm, which harbors M ⊙ of very large (∼millimeter-to-centimeter sized), cold ( K) grains. The central waist (assuming its outer regions to be expanding) and fast bipolar outflow have expansion ages of and : the “jet-lag” (i.e., torus age minus the fast-outflow age) in the Boomerang supports models in which the primary star interacts directly with a binary companion. We argue that this interaction resulted in a common-envelope configuration, while the Boomerang’s primary was an RGB or early-AGB star, with the companion finally merging into the primary’s core, and ejecting the primary’s envelope that now forms the ultra-cold outflow.

Towards understanding the Planck thermal dust models

Understanding the properties of dust emission in the microwave domain is an important premise for the next generation of cosmic microwave background (CMB) experiments, devoted to the measurement of the primordial $B$-modes of polarization. In this paper, we compare three solutions to thermal dust emission by the Planck Collaboration [A. Abergel, P. A. R. Ade et al. (Planck Collaboration), Astron. Astrophys. 571, A11 (2014).,R. Adam, P. A. R. Ade et al. (Planck Collaboration), Astron. Astrophys. 594, A10 (2016)., N. Aghanim, M. Ashdown et al. (Planck Collaboration), Astron. Astrophys. 596, A109 (2016).] to point out significant differences between their respective parameters (the spectral index $\ensuremath{\beta}$, the optical depth $\ensuremath{\tau}$, and the dust temperature ${T}_{d}$). These differences originate from, e.g., the priors on the parameters or the contribution of the cosmic infrared background (CIB). In addition to investigating the angular distributions and statistical properties of each of the $\ensuremath{\beta}$, $\ensuremath{\tau}$, and ${T}_{d}$-maps for the whole sky, we also compute cross-correlations among the maps, specifically the $\ensuremath{\beta}\ensuremath{-}{T}_{d}$ and $\ensuremath{\tau}\ensuremath{-}{T}_{d}$ correlations. All power spectra differ noticeably from each other, which we claim is partly due to the influence of the CIB. Peculiar behavior in the cross-correlations at dust temperatures $\ensuremath{\gtrsim}21\text{ }\text{ }\mathrm{K}$ supports this claim; the precise differences depend on the particular solutions considered. Finally, by the example of two zones on the sky (the BICEP2 zone and a region around the North Celestial Pole), we show that not only the properties of dust are different in these regions on the sky, but moreover the dust emission products do not agree. Furthermore, it is illustrated that the use of average values for dust parameters in one zone will not necessarily be applicable to another zone. In this context, we therefore recommend pixel-based approaches for future analyses, with less stringent constraints in form of priors, despite its higher computational expenditure, and an inclusion of a CIB treatment, which finally allows for a direction dependent removal of dust foregrounds. The central statement of this brief analysis is that while all available solutions are in rough agreement at $\ensuremath{\sim}5%--20%$, further progress must be made to match the goals of planned $B$-mode experiments.

The Python Sky Model: software for simulating the Galactic microwave sky

We present a numerical code to simulate maps of Galactic emission in intensity and polarization at microwave frequencies, aiding in the design of Cosmic Microwave Background experiments. This Python code builds on existing efforts to simulate the sky by providing an easy-to-use interface and is based on publicly available data from the WMAP and Planck satellite missions. We simulate synchrotron, thermal dust, free-free, and anomalous microwave emission over the whole sky, in addition to the Cosmic Microwave Background, and include a set of alternative prescriptions for the frequency dependence of each component that are consistent with current data. We also present a prescription for adding small-scale realizations of these components at resolutions greater than current all-sky measurements. The code is available at https://github.com/bthorne93/PySM_public.

The optical performance of the PILOT instrument from ground end-to-end tests

The Polarized Instrument for Long-wavelength Observation of the Tenuous interstellar medium (PILOT) is a balloon-borne astronomy experiment designed to study the linear polarization of thermal dust emission in two photometric bands centred at wavelengths 240 μm (1.2 THz) and 550 μm (545 GHz), with an angular resolution of a few arcminutes. Several end-to-end tests of the instrument were performed on the ground between 2012 and 2014, in order to prepare for the first scientific flight of the experiment that took place in September 2015 from Timmins, Ontario, Canada. This paper presents the results of those tests, focussing on an evaluation of the instrument’s optical performance. We quantify image quality across the extent of the focal plane, and describe the tests that we conducted to determine the focal plane geometry, the optimal focus position, and sources of internal straylight. We present estimates of the detector response, obtained using an internal calibration source, and estimates of the background intensity and background polarization.

3D radiative transfer of intrinsically polarized dust emission based on aligned aspherical grains

(Sub-)Millimeter observations of the polarized emission of aligned aspherical dust grains enable us to study the magnetic fields within protoplanetary disk. However, the interpretation of these observations is complex. One must consider the various effects that alter the measured polarized signal, such as the shape of dust grains, the efficiency of grain alignment, the magnetic field properties, and the projection of the signal along the line of sight. We aim at analyzing observations of the polarized dust emission by disentangling the effects on the polarization signal in the context of 3D radiative transfer simulations. For this purpose, we developed a code capable of simulating dust grain alignment of aspherical grains and intrinsical polarization of thermal dust emission. We find that the influence of thermal polarization and dust grain alignment on the polarized emission displayed as spatially resolved polarization map or as spectral energy distribution trace disk properties which are not traced in total (unpolarized) emission such as the magnetic field topology. The radiative transfer simulations presented in this work enable the 3D analysis of intrinsically polarized dust emission - observed with, e.g., ALMA - which is essential to constrain magnetic field properties.

Apparent Disk-mass Reduction and Planetisimal Formation in GravitationallyUnstable Disks in Class 0/I Young Stellar Objects

Abstract We investigate the dust structure of gravitationally unstable disks undergoing mass accretion from the envelope, envisioning its application to Class 0/I young stellar objects (YSOs). We find that the dust disk quickly settles into a steady state and that, compared to a disk with interstellar medium (ISM) dust-to-gas mass ratio and micron-sized dust, the dust mass in the steady state decreases by a factor of 1/2 to 1/3, and the dust thermal emission decreases by a factor of 1/3 to 1/5. The latter decrease is caused by dust depletion and opacity decrease owing to dust growth. Our results suggest that the masses of gravitationally unstable disks in Class 0/I YSOs are underestimated by a factor of 1/3 to 1/5 when calculated from the dust thermal emission assuming an ISM dust-to-gas mass ratio and micron-sized dust opacity, and that a larger fraction of disks in Class 0/I YSOs is gravitationally unstable than was previously believed. We also investigate the orbital radius within which planetesimals form via coagulation of porous dust aggregates and show that becomes ∼20 au for a gravitationally unstable disk around a solar mass star. Because increases as the gas surface density increases and a gravitationally unstable disk has maximum gas surface density, is the theoretical maximum radius for planetesimal formation. We suggest that planetesimal formation in the Class 0/I phase is preferable to that in the Class II phase because a large amount of dust is supplied by envelope-to-disk accretion.

Low temperature MIR to submillimeter mass absorption coefficient of interstellar dust analogues

Context.The submillimeter spectral domain has been extensively explored by theHerschelandPlancksatellites and is now reachable from the ground with ALMA. A wealth of data, revealing cold dust thermal emission, is available for astronomical environments ranging from interstellar clouds, cold clumps, circumstellar envelops, and protoplanetary disks. The interpretation of these observations relies on the understanding and modeling of cold dust emission and on the knowledge of the dust optical properties.Aims.The aim of this work is to provide astronomers with a set of spectroscopic data of realistic interstellar dust analogues that can be used to interpret the observations. It pursues the experimental effort aimed at characterizing the spectroscopic properties of interstellar dust analogues at low temperature in the mid-infrared (MIR) to millimeter spectral domain. Compared to previous studies, it extends the range of studied dust analogues in terms of composition and of structure of the material.Methods.Glassy silicates of mean composition (1−x)MgO –xSiO2withx= 0.35 (close to forsterite, Mg2SiO4), 0.50 (close to enstatite, MgSiO3) and 0.40 (close to Mg1.5SiO3.5or MgSiO3:Mg2SiO4= 50:50) were synthesized. The mass absorption coefficient (MAC) of the samples was measured in the spectral domain 30–1000μm for grain temperature in the range 300–10 K and at room temperature in the 5–40μm domain.Results.We find that the MAC of all samples varies with the grains temperature and that its spectral shape cannot be approximated by a single power law inλ−β. In the FIR/submm, and above 30 K, the MAC value at a given wavelength increases with the temperature as thermally activated absorption processes appear. The studied materials exhibit different and complex behaviors at long wavelengths (λ≥ 200 to 700μm depending on the samples). These behaviors are attributed to the amorphous nature of dust and to the amount and nature of the defects within this amorphous structure. We do not observe MAC variations in the 10–30 K range. Above 20μm, the measured MAC are much higher than the MAC calculated from interstellar silicate dust models indicating that the analogues measured in this study are more emissive than the silicates in cosmic dust models.Conclusions.The underestimated value of the MAC deduced from cosmic dust models in the FIR/submm has important astrophysical implications because masses are overestimated by the models. Moreover, constraints on elemental abundance of heavy elements in cosmic dust models are relaxed.

Planck intermediate results

The characterization of the Galactic foregrounds has been shown to be the main obstacle in thechallenging quest to detect primordial B-modes in the polarized microwave sky. We make use of the Planck-HFI 2015 data release at high frequencies to place new constraints on the properties of the polarized thermal dust emission at high Galactic latitudes. Here, we specifically study the spatial variability of the dust polarized spectral energy distribution (SED), and its potential impact on the determination of the tensor-to-scalar ratio, r. We use the correlation ratio of the CBBℓ angular power spectra between the 217 and 353 GHz channels as a tracer of these potential variations, computed on different high Galactic latitude regions, ranging from 80% to 20% of the sky. The new insight from Planck data is a departure of the correlation ratio from unity that cannot be attributed to a spurious decorrelation due to the cosmic microwave background, instrumental noise, or instrumental systematics. The effect is marginally detected on each region, but the statistical combination of all the regions gives more than 99% confidence for this variation in polarized dust properties. In addition, we show that the decorrelation increases when there is a decrease in the mean column density of the region of the sky being considered, and we propose a simple power-law empirical model for this dependence, which matches what is seen in the Planck data. We explore the effect that this measured decorrelation has on simulations of the BICEP2-Keck Array/Planck analysis and show that the 2015 constraints from these data still allow a decorrelation between the dust at 150 and 353 GHz that is compatible with our measured value. Finally, using simplified models, we show that either spatial variation of the dust SED or of the dust polarization angle are able to produce decorrelations between 217 and 353 GHz data similar to the values we observe in the data.

Long-lived Dust Asymmetries at Dead Zone Edges in Protoplanetary Disks

Abstract A number of transition disks exhibit significant azimuthal asymmetries in thermal dust emission. One possible origin for these asymmetries is dust trapping in vortices formed at the edges of dead zones. We carry out high-resolution, two-dimensional hydrodynamic simulations of this scenario, including the effects of dust feedback. We find that, although feedback weakens the vortices and slows down the process of dust accumulation, the dust distribution in the disk can nonetheless remain asymmetric for many thousands of orbits. We show that even after 104 orbits, or 2.5 Myr when scaled to the parameters of Oph IRS 48 (a significant fraction of its age), the dust is not dispersed into an axisymmetric ring, in contrast to the case of a vortex formed by a planet. This is because accumulation of mass at the dead zone edge constantly replenishes the vortex, preventing it from being fully destroyed. We produce synthetic dust emission images using our simulation results. We find that multiple small clumps of dust may be distributed azimuthally. These clumps, if not resolved from one another, appear as a single large feature. A defining characteristic of a disk with a dead zone edge is that an asymmetric feature is accompanied by a ring of dust located about twice as far from the central star.

KINEMATIC STRUCTURE OF MOLECULAR GAS AROUND HIGH-MASS YSO, PAPILLON NEBULA, IN N159 EAST IN THE LARGE MAGELLANIC CLOUD: A NEW PERSPECTIVE WITH ALMA

ABSTRACT We present the ALMA Band 3 and Band 6 results of 12CO(2-1), 13CO(2-1), H30α recombination line, free–free emission around 98 GHz, and the dust thermal emission around 230 GHz toward the N159 East Giant Molecular Cloud (N159E) in the Large Magellanic Cloud (LMC). LMC is the nearest active high-mass star-forming face-on galaxy at a distance of 50 kpc and is the best target for studing high-mass star formation. ALMA observations show that N159E is the complex of filamentary clouds with the width and length of ∼1 pc and several parsecs. The total molecular mass is 0.92 × 105 M ⊙ from the 13CO(2-1) intensity. N159E harbors the well-known Papillon Nebula, a compact high-excitation H ii region. We found that a YSO associated with the Papillon Nebula has the mass of 35 M ⊙ and is located at the intersection of three filamentary clouds. It indicates that the formation of the high-mass YSO was induced by the collision of filamentary clouds. Fukui et al. reported a similar kinematic structure toward two YSOs in the N159 West region, which are the other YSOs that have the mass of ≳35 M ⊙. This suggests that the collision of filamentary clouds is a primary mechanism of high-mass star formation. We found a small molecular hole around the YSO in Papillon Nebula with a sub-parsec scale. It is filled by free–free and H30α emission. The temperature of the molecular gas around the hole reaches ∼80 K. It indicates that this YSO has just started the distruction of parental molecular cloud.

Radio Monitoring of Protoplanetary Discs

Protoplanetary disc systems observed at radio wavelengths often show excess emission above that expected from a simple extrapolation of thermal dust emission observed at short millimetre wavelengths. Monitoring the emission at radio wavelengths can be used to help disentangle the physical mechanisms responsible for this excess, including free-free emission from a wind or jet, and chromospheric emission associated with stellar activity. We present new results from a radio monitoring survey conducted with Australia Telescope Compact Array over the course of several years with observation intervals spanning days, months and years, where the flux variability of 11 T Tauri stars in the Chamaeleon and Lupus star forming regions was measured at 7 and 15 mm and 3 and 6 cm. Results show that for most sources are variable to some degree at 7 mm, indicating the presence of emission mechanisms other than thermal dust in some sources. Additionally, evidence of grain growth to cm-sized pebbles was found for some sources that also have signs of variable flux at 7 mm. We conclude that multiple processes contributing to the emission are common in T Tauri stars at 7 mm and beyond, and that a detection at a single epoch at radio wavelengths should not be used to determine all processes contributing to the emission.

Radio and infrared study of the star-forming region IRAS 20286+4105

A multi-wavelength investigation of the star forming complex IRAS 20286+4105, located in the Cygnus-X region, is presented here. Near-infrared K-band data is used to revisit the cluster / stellar group identified in previous studies. The radio continuum observations, at 610 and 1280 MHz show the presence of a HII region possibly powered by a star of spectral type B0 - B0.5. The cometary morphology of the ionized region is explained by invoking the bow-shock model where the likely association with a nearby supernova remnant is also explored. A compact radio knot with non-thermal spectral index is detected towards the centre of the cloud. Mid-infrared data from the Spitzer Legacy Survey of the Cygnus-X region show the presence of six Class I YSOs inside the cloud. Thermal dust emission in this complex is modelled using Herschel far-infrared data to generate dust temperature and column density maps. Herschel images also show the presence of two clumps in this region, the masses of which are estimated to be {\sim} 175 M{\sun} and 30 M{\sun}. The mass-radius relation and the surface density of the clumps do not qualify them as massive star forming sites. An overall picture of a runaway star ionizing the cloud and a triggered population of intermediate-mass, Class I sources located toward the cloud centre emerges from this multiwavelength study. Variation in the dust emissivity spectral index is shown to exist in this region and is seen to have an inverse relation with the dust temperature.

ALMA BAND 8 CONTINUUM EMISSION FROM ORION SOURCE I

ABSTRACT We have measured continuum flux densities of a high-mass protostar candidate, a radio source I in the Orion KL region (Orion Source I) using the Atacama Large Millimeter/Submillimeter Array (ALMA) at band 8 with an angular resolution of 0.″1. The continuum emission at 430, 460, and 490 GHz associated with Source I shows an elongated structure along the northwest–southeast direction perpendicular to the so-called low-velocity bipolar outflow. The deconvolved size of the continuum source, 90 au × 20 au, is consistent with those reported previously at other millimeter/submillimeter wavelengths. The flux density can be well fitted to the optically thick blackbody spectral energy distribution, and the brightness temperature is evaluated to be 700–800 K. It is much lower than that in the case of proton–electron or H− free–free radiations. Our data are consistent with the latest ALMA results by Plambeck & Wright, in which the continuum emission was proposed to arise from the edge-on circumstellar disk via thermal dust emission, unless the continuum source consists of an unresolved structure with a smaller beam filling factor.

QUANTIFYING THE INTERSTELLAR MEDIUM AND COSMIC RAYS IN THE MBM 53, 54, AND 55 MOLECULAR CLOUDS AND THE PEGASUS LOOP USING FERMI-LAT GAMMA-RAY OBSERVATIONS

ABSTRACT A study of the interstellar medium (ISM) and cosmic rays (CRs) using Fermi Large Area Telescope (LAT) data, in a region encompassing the nearby molecular clouds MBM 53, 54, and 55 and a far-infrared loop-like structure in Pegasus, is reported. By comparing the Planck dust thermal emission model with Fermi-LAT γ-ray data, it was found that neither the dust radiance (R) nor the dust opacity at 353 GHz (τ 353) was proportional to the total gas column density N(Htot) primarily because N(Htot)/R and N(Htot)/τ 353 depend on the dust temperature (T d). The N(Htot) distribution was evaluated using γ-ray data by assuming the regions of high T d to be dominated by optically thin atomic hydrogen ( ) and by employing an empirical linear relation of N(Htot)/R to T d. It was determined that the mass of the gas not traced by the 21 cm or 2.6 mm surveys is ∼25% of the mass of in the optically thin case and is larger than the mass of the molecular gas traced by carbon monoxide by a factor of up to 5. The measured γ-ray emissivity spectrum is consistent with a model based on CR spectra measured at the Earth and a nuclear enhancement factor of ≤1.5. It is, however, lower than local emissivities reported by previous Fermi-LAT studies employing different analysis methods and assumptions on ISM properties by 15%–20% in energies below a few GeV, even if we take account of the statistical and systematic uncertainties. The origin of the discrepancy is also discussed.

THE MAGNETIC FIELD OF L1544. I. NEAR-INFRARED POLARIMETRY AND THE NON-UNIFORM ENVELOPE

ABSTRACT The magnetic field (B-field) of the starless dark cloud L1544 has been studied using near-infrared (NIR) background starlight polarimetry (BSP) and archival data in order to characterize the properties of the plane-of-sky B-field. NIR linear polarization measurements of over 1700 stars were obtained in the H band and 201 of these were also measured in the K band. The NIR BSP properties are correlated with reddening, as traced using the Rayleigh–Jeans color excess (H–M) method, and with thermal dust emission from the L1544 cloud and envelope seen in Herschel maps. The NIR polarization position angles change at the location of the cloud and exhibit their lowest dispersion there, offering strong evidence that NIR polarization traces the plane-of-sky B-field of L1544. In this paper, the uniformity of the plane-of-sky B-field in the envelope region of L1544 is quantitatively assessed. This allows evaluation of the approach of assuming uniform field geometry when measuring relative mass-to-flux ratios in the cloud envelope and core based on averaging of the radio Zeeman observations in the envelope, as done by Crutcher et al. In L1544, the NIR BSP shows the envelope B-field to be significantly non-uniform and likely not suitable for averaging Zeeman properties without treating intrinsic variations. Deeper analyses of the NIR BSP and related data sets, including estimates of the B-field strength and testing how it varies with position and gas density, are the subjects of later papers in this series.