Context. Low-energy cosmic rays are the dominant source of ionization for molecular cloud cores. The ionization fraction, in turn, controls the coupling of the magnetic field to the gas and hence the ...dynamical evolution of the cores. Aims. The purpose of this work is to compute the attenuation of the cosmic-ray flux rate in a cloud core taking into account magnetic focusing, magnetic mirroring, and all relevant energy loss processes. Methods. We adopt a standard cloud model characterized by a mass-to-flux ratio supercritical by a factor of ~2 to describe the density and magnetic field distribution of a low-mass starless core, and we follow the propagation of cosmic rays through the core along flux tubes enclosing different amount of mass. We then extend our analysis to cores with different mass-to-flux ratios. Results. We find that mirroring always dominates over focusing, implying a reduction of the cosmic-ray ionization rate by a factor of ~2−3 over most of a solar-mass core with respect to the value in the intercloud medium outside the core. For flux tubes enclosing larger masses the reduction factor is smaller, since the field becomes increasingly uniform at larger radii and lower densities. We also find that the cosmic-ray ionization rate is further reduced in clouds with stronger magnetic field, e.g. by a factor ~4 for a marginally critical cloud. Conclusions. The magnetic field threading molecular cloud cores affects the penetration of low-energy cosmic rays and reduces the ionization rate by a factor 3 − 4 depending on the position inside the core and the magnetization of the core.
Cosmic-ray ionization of molecular clouds Padovani, M.; Galli, D.; Glassgold, A. E.
Astronomy and astrophysics (Berlin),
07/2009, Letnik:
501, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Context. Low-energy cosmic rays are a fundamental source of ionization for molecular clouds, influencing their chemical, thermal, and dynamical evolution. Aims. The purpose of this work is to explore ...the possibility that a low-energy component of cosmic rays, not directly measurable from the Earth, can account for the discrepancy between the ionization rate measured in diffuse and dense interstellar clouds. Methods. We collected the most recent experimental and theoretical data on the cross sections for the production of H$_2^+$ and He+ by electron and proton impact and discuss the available constraints on the cosmic-ray fluxes in the local interstellar medium. Starting from different extrapolations at low energies of the demodulated cosmic-ray proton and electron spectra, we computed the propagated spectra in molecular clouds in the continuous slowing-down approximation taking all the relevant energy loss processes into account. Results. The theoretical value of the cosmic-ray ionization rate as a function of the column density of traversed matter agrees with the observational data only if the flux of either cosmic-ray electrons or of protons increases at low energies. The most successful models are characterized by a significant (or even dominant) contribution of the electron component to the ionization rate, in agreement with previous suggestions. However, the large spread of cosmic-ray ionization rates inferred from chemical models of molecular cloud cores remains to be explained. Conclusions. Available data combined with simple propagation models support the existence of a low-energy component (below ~100 MeV) of cosmic-ray electrons or protons responsible for the ionization of molecular cloud cores and dense protostellar envelopes.
Context. Cosmic rays play an important role in dense molecular cores, affecting their thermal and dynamical evolution and initiating the chemistry. Several studies have shown that the formation of ...protostellar discs in collapsing clouds is severely hampered by the braking torque exerted by the entrained magnetic field on the infalling gas, as long as the field remains frozen to the gas. Aims. In this paper we examine the possibility that the concentration and twisting of the field lines in the inner region of collapse can produce a significant reduction of the ionisation fraction. Methods. To check whether the cosmic-ray ionisation rate can fall below the critical value required to maintain good coupling, we first study the propagation of cosmic rays in a model of a static magnetised cloud varying the relative strength of the toroidal/poloidal components and the mass-to-flux ratio. We then follow the path of cosmic rays using realistic magnetic field configurations generated by numerical simulations of a rotating collapsing core with different initial conditions. Results. We find that an increment of the toroidal component of the magnetic field, or, in general, a more twisted configuration of the field lines, results in a decrease in the cosmic-ray flux. This is mainly due to the magnetic mirroring effect that is stronger where larger variations in the field direction are present. In particular, we find a decrease of the cosmic-ray ionisation rate below 10-18 s-1 in the central 300–400 AU, where density is higher than about 109 cm-3. This very low value of the ionisation rate is attained in the cases of intermediate and low magnetisation (mass-to-flux ratio λ = 5 and 17, respectively) and for toroidal fields larger than about 40% of the total field. Conclusions. Magnetic field effects can significantly reduce the ionisation fraction in collapsing clouds. We provide a handy fitting formula to compute approximately the attenuation of the cosmic-ray ionisation rate in a molecular cloud as a function of the density and the magnetic configuration.
Context. Polarized continuum emission at millimeter-to-submillimeter wavelengths is usually attributed to thermal emission from dust grains aligned through radiative torques with the magnetic field. ...However, recent theoretical work has shown that under specific conditions polarization may arise from self-scattering of thermal emission and by radiation fields from a nearby stellar object. Aims. We use multi-frequency polarization observations of a circumbinary disk to investigate how the polarization properties change at distinct frequency bands. Our goal is to discern the main mechanism responsible for the polarization through comparison between our observations and model predictions for each of the proposed mechanisms. Methods. We used the Atacama Large Millimeter/submillimeter Array to perform full polarization observations at 97.5 GHz (Band 3), 233 GHz (Band 6) and 343.5 GHz (Band 7). The ALMA data have a mean spatial resolution of 28 AU. The target is the Class I object BHB07-11, which is the youngest object in the Barnard 59 protocluster. Complementary Karl G. Jansky Very Large Array observations at 34.5 GHz were also performed and revealed a binary system at centimetric continuum emission within the disk. Results. We detect an extended and structured polarization pattern that is remarkably consistent between the three bands. The distribution of polarized intensity resembles a horseshoe shape with polarization angles following this morphology. From the spectral index between Bands 3 and 7, we derived a dust opacity index β ~ 1 consistent with maximum grain sizes larger than expected to produce self-scattering polarization in each band. The polarization morphology and the polarization levels do not match predictions from self-scattering. On the other hand, marginal correspondence is seen between our maps and predictions from a radiation field model assuming the brightest binary component as main radiation source. Previous molecular line data from BHB07-11 indicates disk rotation. We used the DustPol module of the ARTIST radiative transfer tool to produce synthetic polarization maps from a rotating magnetized disk model assuming combined poloidal and toroidal magnetic field components. The magnetic field vectors (i.e., the polarization vectors rotated by 90°) are better represented by a model with poloidal magnetic field strength about three times the toroidal one. Conclusions. The similarity of our polarization patterns among the three bands provides a strong evidence against self-scattering and radiation fields. On the other hand, our data are reasonably well reproduced by a model of disk with toroidal magnetic field components slightly smaller than poloidal ones. The residual is likely to be due to the internal twisting of the magnetic field due to the binary system dynamics, which is not considered in our model.
ABSTRACT Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These ...reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.
Context. Submillimeter Array (SMA) 870 μm polarization observations of the hot molecular core G31.41+0.31 revealed one of the clearest examples up to date of an hourglass-shaped magnetic field ...morphology in a high-mass star-forming region. Aims. To better establish the role that the magnetic field plays in the collapse of G31.41+0.31, we carried out Atacama Large Millimeter/ submillimeter Array (ALMA) observations of the polarized dust continuum emission at 1.3 mm with an angular resolution four times higher than that of the previous (sub)millimeter observations to achieve an unprecedented image of the magnetic field morphology. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6). The resulting synthesized beam is 0′′.28×0′′.20 which, at the distance of the source, corresponds to a spatial resolution of ~875 au. Results. The observations resolve the structure of the magnetic field in G31.41+0.31 and allow us to study the field in detail. The polarized emission in the Main core of G31.41+0.41is successfully fit with a semi-analytical magnetostatic model of a toroid supported by magnetic fields. The best fit model suggests that the magnetic field is well represented by a poloidal field with a possible contribution of a toroidal component of ~10% of the poloidal component, oriented southeast to northwest at approximately −44° and with an inclination of approximately −45°. The magnetic field is oriented perpendicular to the northeast to southwest velocity gradient detected in this core on scales from 103 to 104 au. This supports the hypothesis that the velocity gradient is due to rotation of the core and suggests that such a rotation has little effect on the magnetic field. The strength of the magnetic field estimated in the central region of the core with the Davis–Chandrasekhar-Fermi method is ~8–13 mG and implies that the mass-to-flux ratio in this region is slightly supercritical. Conclusions. The magnetic field in G31.41+0.31 maintains an hourglass-shaped morphology down to scales of <1000 au. Despite the magnetic field being important in G31.41+0.31, it is not enough to prevent fragmentation and collapse of the core, as demonstrated by the presence of at least four sources embedded in the center of the core.
Context. The existence of disks around high-mass stars has yet to be established on a solid ground, as only few reliable candidates are known to date. The disk rotating about the ~104 L⊙ protostar ...IRAS 20126+4104 is probably the most convincing of these. Aims. We would like to resolve the disk structure in IRAS 20126+4104 and, if possible, investigate the relationship between the disk and the associated jet emitted along the rotation axis. Methods. We performed observations at 1.4 mm with the IRAM Plateau de Bure interferometer attaining an angular resolution of ~0.̋4 (~660 AU). We imaged the methyl cyanide J = 12 → 11 ground state and vibrationally excited transitions as well as the CH313CN isotopologue, which had proved to be disk tracers. Results. Our findings confirm the existence of a disk rotating about a ~7–10 M⊙ star in IRAS 20126+4104, with rotation velocity increasing at small radii. The dramatic improvement in sensitivity and spectral and angular resolution with respect to previous observations allows us to establish that higher excitation transitions are emitted closer to the protostar than the ground state lines, which demonstrates that the gas temperature is increasing towards the centre. We also find that the material is asymmetrically distributed in the disk and speculate on the possible origin of such a distribution. Finally, we demonstrate that the jet emitted along the disk axis is co-rotating with the disk. Conclusions. We present iron-clad evidence of the existence of a disk undergoing rotation around a B-type protostar, with rotation velocity increasing towards the centre. We also demonstrate that the disk is not axially symmetric. These results prove that B-type stars may form through disk-mediated accretion as their low-mass siblings do, but also show that the disk structure may be significantly perturbed by tidal interactions with (unseen) companions, even in a relatively poor cluster such as that associated with IRAS 20126+4104.
Context.
Cosmic rays (CRs), which are energetic particles mainly composed of protons and electrons, play an important role in the chemistry and dynamics of the interstellar medium. In dense ...environments, they represent the main ionising agent, hence driving the rich chemistry of molecular ions. Furthermore, they determine the ionisation fraction, which regulates the degree of coupling between the gas and the interstellar magnetic fields, and the heating of the gas. Estimates of the CR ionisation rate of molecular hydrogen (
ζ
2
) span several orders of magnitude, depending on the targeted sources and on the method used.
Aims.
Recent theoretical models have characterised the CR attenuation with increasing density. We aim to test these models for the attenuation of CRs in the low-mass pre-stellar core L1544.
Methods.
We used a state-of-the-art gas-grain chemical model, which accepts the CR ionisation rate profile as input, to predict the abundance profiles of four ions: N
2
H
+
, N
2
D
+
, HC
18
O
+
, and DCO
+
. Non-local thermodynamic equilibrium radiative transfer simulations were run to produce synthetic spectra based on the derived abundances. These were compared with observations obtained with the Institut de Radioastronomie Millimétrique 30 m telescope.
Results.
Our results indicate that a model with high
ζ
2
(>10
−16
s
−1
) is excluded by the observations. Also the model with the standard
ζ
2
= 1.3 × 10
−17
s
−1
produces a worse agreement with respect to the attenuation model based on Voyager observations, which is characterised by an average ⟨
ζ
2
⟩ = 3 × 10
−17
s
−1
at the column densities typical of L1544. The single-dish data, however, are not sensitive to the attenuation of the CR profile, which changes only by a factor of two in the range of column densities spanned by the core model (
N
= 2−50 × 10
21
cm
−2
). Interferometric observations at higher spatial resolution, combined with observations of transitions with lower critical density – hence tracing the low-density envelope – are needed to observe a decrease in the CR ionisation rate with density.
Abstract
Nanostructured Au films fabricated by the assembling of nanoparticles produced in the gas phase have shown properties suitable for neuromorphic computing applications: they are characterized ...by a non-linear and non-local electrical behavior, featuring switches of the electric resistance whose activation is typically triggered by an applied voltage over a certain threshold. These systems can be considered as complex networks of metallic nanojunctions where thermal effects at the nanoscale cause the continuous rearrangement of regions with low and high electrical resistance. In order to gain a deeper understanding of the electrical properties of this nano granular system, we developed a model based on a large three dimensional regular resistor network with non-linear conduction mechanisms and stochastic updates of conductances. Remarkably, by increasing enough the number of nodes in the network, the features experimentally observed in the electrical conduction properties of nanostructured gold films are qualitatively reproduced in the dynamical behavior of the system. In the activated non-linear conduction regime, our model reproduces also the growing trend, as a function of the subsystem size, of quantities like Mutual and Integrated Information, which have been extracted from the experimental resistance series data via an information theoretic analysis. This indicates that nanostructured Au films (and our model) possess a certain degree of activated interconnection among different areas which, in principle, could be exploited for neuromorphic computing applications.
As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41 +0.31 with ALMA at 1.4 mm with an angular resolution of ~0.′′22 (~1700 ...au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE–SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up toward the center. Spectral lines towardtoward the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating toward the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.