Search NASA⌕ Search

SEARCH · Search NASA

Results for “interstellar clouds”

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 19 records

CF + excitation in the interstellar medium

The detection of CF + in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF + abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF + in collisions with para- and ortho-H 2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF + emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF + as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Elementary Composite Binary and Grain Alignment Locked in Dust Growth

Abstract Planets are known to grow out of a star-encircling disk of the gas and dust inherited from an interstellar cloud; their formation is thought to begin with coagulation of submicron dust grains into aggregates, the first foundational stage of planet formation. However, with nanoscale and submicron solids unobservable directly in the interstellar medium (ISM) and protoplanetary disks, how dust grains grow is unclear, as are the morphology and structure of interstellar grains and the whereabouts and form of “missing iron.” Here we show an elementary composite binary in 3D sub-10 nm detail—and the alignments of its two subunits and nanoinclusions and a population of elongated composite grains locked in a primitive cosmic dust particle—noninvasively uncovered with phase-contrast X-ray nanotomography. The binary comprises a pair of oblate, quasi-spheroidal grains whose alignment and shape meet the astrophysical constraints on polarizing interstellar grains. Each member of the pair contains a high-density core of octahedral nanocrystals whose twin relationship is consistent with the magnetite’s diagnostic property at low temperatures, with a mantle exhibiting nanoscale heterogeneities, rounded edges, and pitted surfaces. This elongated binary evidently formed from an axially aligned collision of the two similar composite grains whose core–mantle structure and density gradients are consistent with interstellar processes and astronomical evidence for differential depletion. Our findings suggest that the ISM is threaded with dust grains containing preferentially oriented iron-rich magnetic nanocrystals that hold answers to astronomical problems from dust evolution, grain alignment, and the structure of magnetic fields to planetesimal growth.

Hu, Z. W. (ORCID:0000000278648830)↗

Polycyclic aromatic hydrocarbons in samples of Ryugu formed in the interstellar medium

Polycyclic aromatic hydrocarbons (PAHs) contain ≲20% of the carbon in the interstellar medium. They are potentially produced in circumstellar environments (at temperatures ≳1000 kelvin), by (~10 kelvin) interstellar clouds, or by processing of carbon-rich dust grains. We report isotopic properties of PAHs extracted from samples of the asteroid Ryugu and the meteorite Murchison. The doubly- 13 C substituted compositions (Δ2× 13 C values) of the PAHs naphthalene, fluoranthene, and pyrene are 9 to 51‰ higher than values expected for a stochastic distribution of isotopes. The Δ2× 13 C values are higher than expected if the PAHs formed in a circumstellar environment, but consistent with formation in the interstellar medium. By contrast, the PAHs phenanthrene and anthracene in Ryugu samples have Δ2× 13 C values consistent with formation by higher-temperature reactions.

Science & Technology - Other Topics↗

Whence the Interstellar Magnetic Field Shaping the Heliosphere?

Measurements of starlight polarized by aligned interstellar dust grains are used to probe the relation between the orientation of the ambient interstellar magnetic field (ISMF) and the ISMF traced by the ribbons of energetic neutral atoms discovered by the Interstellar Boundary Explorer spacecraft. We utilize polarization data, many acquired specifically for this study, to trace the configuration of the ISMF within 40 pc. A statistical analysis yields a best-fit ISMF orientation, B magpol , aligned with Galactic coordinates ℓ = 42°, b = 49°. Further analysis shows the ISMF is more orderly for "downfield" stars located over 90° from Bmagpol. The data subset of downfield stars yields an orientation for the nearby ISMF at ecliptic coordinates λ, β ≈ 219° ± 15°, 43° ± 9° (Galactic coordinates l, b ≈ 40°, 56°, ±17°). This best-fit ISMF orientation from polarization data is close to the field direction obtained from ribbon models. This agreement suggests that the ISMF shaping the heliosphere belongs to an extended ordered magnetic field. Extended filamentary structures are found throughout the sky. A previously discovered filament traversing the heliosphere nose region, "Filament A," extends over 300° of the sky, and crosses the upwind direction of interstellar dust flowing into the heliosphere. Filament A overlaps the locations of the Voyager kilohertz emissions, three quasar intraday variables, cosmic microwave background (CMB) components, and the inflow direction of interstellar grains sampled by Ulysses and Galileo. These features are likely located in the upstream outer heliosheath where ISMF drapes over the heliosphere, suggesting Filament A coincides with a dusty magnetized plasma. A filament 55° long is aligned with a possible shock interface between local interstellar clouds. A dark spot in the CMB is seen within 5° of the filament and within 10° of the downfield ISMF direction. Two large magnetic arcs are centered on the directions of the heliotail. The overlap between CMB components and the aligned dust grains forming Filament A indicates the configuration of dust entrained in the ISMF interacting with the heliosphere provides a measurable foreground to the CMB.

79 ASTRONOMY AND ASTROPHYSICS↗

Collisional excitation of interstellar PN by H 2 : New interaction potential and scattering calculations

Here, rotational excitation of interstellar PN molecules induced by collisions with H 2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN–H 2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H 2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm -1 . The dissociation energies were found to be 40.19 cm -1 and 75.05 cm -1 for complexes of PN with ortho-H 2 and para-H 2 , respectively. Integral cross sections for rotational excitation in PN–H 2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H 2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic survey of faint planetary-nebula nuclei – II. The subdwarf O central star of Fr 2-30

Fr 2-30 = PN? G126.8−15.5 is a faint emission nebula, hosting a 14th-mag central star that we identify here for the first time. Deep Hα and [O III ] images reveal a roughly elliptical nebula with dimensions of at least 22 arcmin × 14 arcmin, fading into a surrounding network of even fainter emission. Optical spectrograms of the central star show it to have a subdwarf O spectral type, with a Gaia parallax distance of 890 pc. A model-atmosphere analysis gives parameters of T eff = 60 000 K, log g = 6.0, and a low helium content of n He /n H = 0.0017. The location of the central star in the log g–T eff plane is inconsistent with a post-asymptotic-giant-branch evolutionary status. Two alternatives are that it is a helium-burning post-extreme-horizontal-branch object, or a hydrogen-burning post-red-giant-branch star. In either case, the evolutionary ages are so long that a detectable planetary nebula (PN) should not be present. We find evidence for a variable radial velocity (RV), suggesting that the star is a close binary. However, there are no photometric variations, and the spectral-energy distribution rules out a companion earlier than M2 V. The RVs of the star and surrounding nebula are discordant, and the nebula lacks typical PN morphology. We suggest that Fr 2-30 is a ‘PN mimic’ – the result of a chance encounter between the hot sdO star and an interstellar cloud. However, we note the puzzling fact that there are several nuclei of genuine PNe that are known to be in evolutionary states similar to that of the Fr 2-30 central star.

79 ASTRONOMY AND ASTROPHYSICS↗

Near L-edge Single and Multiple Photoionization of Doubly Charged Iron Ions

Using the photon-ion merged-beams technique at a synchrotron light source, we have measured relative cross sections for single and up to five-fold photoionization of Fe 2+ ions in the energy range of 690-920 eV. This range contains thresholds and resonances associated with ionization and excitation of 2p and 2s electrons. Calculations were performed to simulate the total absorption spectra. The theoretical results show very good agreement with the experimental data, if overall energy shifts of up to 2.5 eV are applied to the calculated resonance positions and assumptions are made about the initial experimental population of the various levels of the Fe 2+ ([Ar]3d 6 ) ground configuration. Furthermore, we performed extensive calculations of the Auger cascades that result when an electron is removed from the 2p subshell of Fe 2+ . These computations lead to a better agreement with the measured product-charge-state distributions as compared to earlier work. We conclude that the L-shell absorption features of low-charged iron ions are useful for identifying gas-phase iron in the interstellar medium and for discriminating against the various forms of condensed-phase iron bound to composite interstellar dust grains.

79 ASTRONOMY AND ASTROPHYSICS↗

Structural organization of space polymers

Extraterrestrial polymers of glycine with iron have been characterized by mass spectrometry to have a core mass of 1494 Da with dominant rod-like variants at mass to charge ratios of 1567 and 1639 [McGeoch et al., “Meteoritic proteins with glycine, iron and lithium,” arXiv:2102.10700 (2021)]. Several principal macro-structural morphologies are observed in solvent extracts from a chondritic Vigarano class alteration type 3 meteoritic material. The first is an extended sheet of linked (three-legged) triskelia containing the 1494 Da core entity that encloses gas bubbles in the solvent. A second is of fiber-like crystals found here, via x-ray diffraction, to be multiple-walled nanotubes made from a square lattice of the 1494 Da polymer. A third is a dispersion of floating phantom-like short tubes of up to 100 μm length with characteristic angled bends that suggest the influence of a specific underlying protein structure. Here, it is proposed that the angled tubes are the observable result of a space-filling superpolymerization of 1638 Da polymer subunits guided by the tetragonal symmetry of linking silicon bonds. Distorted hexagonal sheets are linked by perpendicular subunits in a three-dimensional hexagonal diamond structure to fill the largest possible volume. This extended very low-density structure is conjectured to have dominated in a process of chemical selection because it captured a maximum amount of molecular raw material in the ultra-low density of molecular clouds or of the proto-solar nebula. This could have led ultimately to the accretion of the earliest planetary bodies.

59 BASIC BIOLOGICAL SCIENCES↗

Synthetic Absorption Lines from Simulations of Multiphase Gas in Galactic Winds

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high-resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T$ ~ 10 4 K supernova-driven wind and a cool $T$ ~ 10 7 K cloud of interstellar material, we create mock observables, including the optical depth ($\tau$) and covering fraction ($C$ $f$ ) of six commonly observed ions (Si II , C II , Si IV , C IV , N V , and O VI ) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth and overestimate the "true" value of $\tau$ for ions with low optical depth relative to the simulated data. In conclusion, the empirically estimated covering fractions match our direct calculations very well for the low-ionization ions; for the high-ionization ions, the empirical covering fractions underestimate the directly computed values by up to ~40%.

79 ASTRONOMY AND ASTROPHYSICS↗

Conditional H I Mass Functions and the H I -to-halo Mass Relation in the Local Universe

We present a new H I mass estimator that relates log 10 (M HI /M * ) to a linear combination of four galaxy properties: stellar surface mass density, color index u - r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H I mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H I detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H I masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H I mass function (HIMF) of local galaxies, as well as the conditional H I mass function in galaxy groups and the H I –halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 10 9 M ⊙ , but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H I mass in a group increases monotonically with halo mass, but for central galaxies, the H I mass shows weak dependence on halo mass when M h ≳ 10 12 M ⊙ . The observed H I –halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high- and low-density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining dark matter-proton scattering from molecular cloud ionization

Optically dense clouds in the interstellar medium composed predominantly of molecular hydrogen, known as molecular clouds, are sensitive to energy injection in the form of photon absorption, cosmic-ray scattering, and dark matter (DM) scattering. The ionization rates in dense molecular clouds are heavily constrained by observations of abundances of various molecular tracers. Recent studies have set constraints on the DM-electron scattering cross section using measurements of ionization rates in dense molecular clouds. Here we calculate the analogous bounds on the DM-proton cross section using the molecular Migdal effect, recently adapted from the neutron scattering literature to the DM context. These bounds may be the strongest limits on a strongly coupled DM subfraction, and represent the first application of the Migdal effect to astrophysical systems. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Physical and dynamical characterization of hyperbolic comet C/2017 U7 (PANSTARRS)

Here we present here a dynamical and observational study of the comet C/2017 U7 (PANSTARRS). This comet was discovered in 2017 and found to have a hyperbolic orbit. Our dynamical analysis shows that the object has probably originated in the Oort cloud, however an interstellar origin cannot be discarded. The observations were obtained in 2018 and 2019 using the Goodman High Throughput Spectrograph (GHTS) at the SOAR telescope. We obtained visible spectra covering the wavelength range of 0.5 - 0.9 μm and also images in the SDSS filters system. Both the low-resolution reflectance spectrum and the reflectance spectra derived from the SDSS filters show an atypical band at ~ 0.595 μm. We conducted a comparative study of the colors and reflectance spectra of different small body populations (e.g., comets, Centaurs, and trans-Neptunian objects or TNOs) from the literature and concluded that the spectra and the colors of this comet are atypical, showing only some overlap with those of some known members of the TNOs and Centaurs, within the large uncertainties of the measurements of those populations. It is found that the feature and overall spectral shape can be reproduced by laboratory spectra of kerite, a template for aliphatic-rich hydrocarbons that has been previously identified in NIR cometary spectra absorptions. It is tentatively proposed that the unusual spectral shape is the result of a particle size distribution of dust grains in the coma or on the surface that has arisen due to a low grain ejection velocity from the surface and large nucleus size.

79 ASTRONOMY AND ASTROPHYSICS↗

Directed Gas-Phase Formation of The Propargyl Family of Resonance-Stabilized Radicals in The Reactions of Ground-State Carbon Atoms (C; 3 P j ) with Butene Isomers (C 4 H 8 ): Dimethylpropargyl and Ethylpropargyl

The propargyl radical (C 3 H 3 ) is the simplest resonance-stabilized free radical (RSFR), but how does stepwise methyl substitution in the alkene reactant affect its dynamics of their formation? We report a crossed molecular beam study of the reactions of atomic carbon (C, 3 P j ) with four butene isomers (C 4 H 8 ) under single collision conditions at a collision energy of 28 ± 2 kJ mol −1 . Barrierless addition of atomic carbon to the alkene C=C bond triggers ring opening to triplet substituted allenes—a de facto insertion mechanism—followed by unimolecular decomposition via atomic hydrogen (H), methyl (CH 3 ) or ethyl (C 2 H 5 ) loss, yielding a family of propargyl‑type RSFRs. RRKM calculations reveal that the branching ratios are highly sensitive to the alkene structure. While the methyl loss channel, affording 1‑methylpropargyl, dominates for 2 butenes (80–90%), the predicted hydrogen atom loss channel (≈10%) leading to 1,3 dimethylpropargyl is identified in the experiment by comparison with theoretical energetics. For isobutene, a near‑equal competition is seen, with the reaction producing 3‑methylpropargyl (≈50%) and 1,1‑dimethylpropargyl (≈40%), along with 2‑vinylallyl (≈5%), whose formation is supported by the experimental data. Most notably, the reaction with 1 butene uniquely favors an enthalpically driven hydrogen shift, eventually producing 1 vinylallyl (≈38%), which is assigned based on the excellent agreement between the measured and calculated reaction exothermicity. Rapid entropically favored fragmentation channels yield ≈40% of propargyl type species (propargyl, 1- and 3-ethylpropargyls), slightly outcompeting the allyl type product. Furthermore, these results establish a systematic progression from C 2 H 4 to C 4 H 8 , where the increasing alkyl substitution unlocks new fragmentation channels, providing a versatile gas phase route to alkylated RSFRs—key intermediates in the growth of methylated and ethylated PAHs and aliphatic chains in combustion and cold interstellar environments (molecular clouds).

Alkyls↗

CO and [C ii] line emission of molecular clouds: the impact of stellar feedback and non-equilibrium chemistry

We analyse synthetic 12 CO, 13 CO, and [C ii] emission maps of molecular cloud (MC) simulations from the SILCC-Zoom project. We present radiation, magnetohydrodynamic zoom-in simulations of individual clouds, both with and without radiative stellar feedback, forming in a turbulent multiphase interstellar medium following on-the-fly the evolution of e.g. H 2 , CO, and C + . We introduce a novel post-processing routine based on cloudy which accounts for higher ionization states of carbon due to stellar radiation in H ii regions. Synthetic emission maps of [C ii] in and around feedback bubbles show that the bubbles are largely devoid of [C ii], as recently found in observations, which we attribute to the further ionization of C+ into C 2+ . For both 12 CO and 13 CO, the cloud-averaged luminosity ratio, $L_\rm {CO}/L_\rm {[C\, \small {II}]}$, can neither be used as a reliable measure of the H 2 mass fraction nor of the evolutionary stage of the clouds. We note a relation between the $I_\rm {CO}/I_\rm {[C\, \small {II}]}$ intensity ratio and the H 2 mass fraction for individual pixels of our synthetic maps. The scatter, however, is too large to reliably infer the H 2 mass fraction. Finally, the assumption of chemical equilibrium overestimates H2 and CO masses by up to 150 and 50 per cent, respectively, and $L_\rm {CO}$ by up to 60 per cent. The masses of H and C + would be underestimated by 65 and 30 per cent, respectively, and $L_\rm {[C\, \small {II}]}$ by up to 35 per cent. Hence, the assumption of chemical equilibrium in MC simulations introduces intrinsic errors of a factor of 2 in chemical abundances, luminosities, and luminosity ratios.

79 ASTRONOMY AND ASTROPHYSICS↗

Applying the Velocity Gradient Technique in NGC 1333: Comparison with Dust Polarization Observations

Magnetic fields (B-fields) are ubiquitous in the interstellar medium (ISM), and they play an essential role in the formation of molecular clouds and subsequent star formation. However, B-fields in interstellar environments remain challenging to measure, and their properties typically need to be inferred from dust polarization observations over multiple physical scales. In this work, we seek to use a recently proposed approach called the velocity gradient technique (VGT) to study B-fields in star-forming regions and compare the results with dust polarization observations in different wavelengths. The VGT is based on the anisotropic properties of eddies in magnetized turbulence to derive B-field properties in the ISM. We investigate that this technique is synergistic with dust polarimetry when applied to a turbulent diffused medium for the purpose of measuring its magnetization. Specifically, we use the VGT on molecular line data toward the NGC 1333 star-forming region ( 12 CO, 13 CO, C 18 O, and N 2 H + ), and we compare the derived B-field properties with those inferred from 214 and 850 μm dust polarization observations of the region using Stratospheric Observatory for Infrared Astronomy/High-Resolution Airborne Wide-band Camera Plus and James Clerk Maxwell Telescope/POL-2, respectively. We estimate both the inclination angle and the 3D Alfvénic Mach number M A from the molecular line gradients. Crucially, testing this technique on gravitationally bound, dynamic, and turbulent regions, and comparing the results with those obtained from polarization observations at different wavelengths, such as the plane-of-sky field orientation, is an important test on the applicability of the VGT in various density regimes of the ISM. We in general do not find a close correlation between the velocity gradient inferred orientations and the dust inferred magnetic field orientations.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmic rays in molecular clouds probed by H 2 rovibrational lines

Context. Low-energy cosmic rays (<1 TeV) play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H 2 . Their characterisation is therefore important both for the interpretation of observations and for the development of theoretical models. However, the methods used so far for estimating the cosmic-ray ionisation rate in molecular clouds have several limitations due to uncertainties in the adopted chemical networks. Aims. Here, we refine and extend a previously proposed method to estimate the cosmic-ray ionisation rate in molecular clouds by observing rovibrational transitions of H 2 at near-infrared wavelengths, which are mainly excited by secondary cosmic-ray electrons. Methods. Combining models of interstellar cosmic-ray propagation and attenuation in molecular clouds with the rigorous calculation of the expected secondary electron spectrum and updated electron-H 2 excitation cross sections, we derive the intensity of the four H 2 rovibrational transitions observable in cold dense gas: (1—0)O(2), (1—0)Q(2), (1—0)S(0), and (1—0)O(4). Results. The proposed method allows the estimation of the cosmic-ray ionisation rate for a given observed line intensity and H 2 column density. We are also able to deduce the shape of the low-energy cosmic-ray proton spectrum impinging upon the molecular cloud. In addition, we present a look-up plot and a web-based application that can be used to constrain the low-energy spectral slope of the interstellar cosmic-ray proton spectrum. We finally comment on the capability of the James Webb Space Telescope to detect these near-infrared H 2 lines, which will make it possible to derive, for the first time, spatial variation in the cosmic-ray ionisation rate in dense gas. Besides the implications for the interpretation of the chemical-dynamic evolution of a molecular cloud, it will finally be possible to test competing models of cosmic-ray propagation and attenuation in the interstellar medium, as well as compare cosmic-ray spectra in different Galactic regions

79 ASTRONOMY AND ASTROPHYSICS↗

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗