Properties of the infrared dark cloud G79.2+0.38

The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the ¹²CO and its isotopic molecules ¹³CO and ¹⁸CO. The dense molecular cores defined by the C¹⁸O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H₂) (5 – 12) × 10²² cm⁻² and the mean number density n (3 ± 1) × 10⁴ cm⁻³. Their sizes are 1.7 and 1.2 pc in ¹³CO(1-0) line, 1.2 and 0.6 pc in C¹⁸O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 10² to 2 × 10³ M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p^{−0.34±0.02}. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules ¹³CO and C¹⁸O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X_13/18 between ¹³CO and C¹⁸O with the column density.     

Estimation prospects of the source number density of ultra-high-energy cosmic rays

We discuss the possibility of accurately estimating the source number density of ultra-high-energy cosmic rays (UHECRs) using small-scale anisotropy in their arrival distribution. The arrival distribution has information on their source and source distribution. We calculate the propagation of UHE protons in a structured extragalactic magnetic field (EGMF) and simulate their arrival distribution at the Earth using our previously developed method. The source number density that can best reproduce observational results by Akeno Giant Air Shower Array is estimated at about 10⁻⁵ Mpc⁻³ in a simple source model. Despite having large uncertainties of about one order of magnitude, due to small number of observed events in current status, we find that more detection of UHECRs in the Auger era can sufficiently decrease this so that the source number density can be more robustly estimated. Two hundred event observation above 4 × 10¹⁹ eV in a hemisphere can discriminate between 10⁻⁵ and 10⁻⁶ Mpc⁻³. Number of events to discriminate between 10⁻⁴ and 10⁻⁵ Mpc⁻³ is dependent on EGMF strength. We also discuss the same in another source model in this paper.     

Neutral and ion-exospheres in the solar wind with applications to mercury

A model of planetary neutral and ion-exospheres in the solar wind is formulated for weak or lunar like solar-wind interaction with a planet. The neutral exosphere model allows for density and temperature variations and for rotation at the exobase. The ion-exosphere is produced by ionization of the neutral exosphere in the solar wind and its density distribution is obtained by solving the continuity equation in the drift approximation. Applying to Mercury a surface temperature distribution inferred from infra-red data and a vanishing bound neutral flux at the base, He and He⁺ density distributions are found. When the He atmosphere of Mercury is due entirely to the surface bombardment by solar wind He⁺⁺, the resulting He⁺ density is found to vary from 1.5 × 10⁻¹ to 10⁻³ cm⁻³ over the range 1.5–5 planetocentric radii on the dayside. These densities are found to be detectable by typical solar-wind plasma instruments. The possible effects of cyclotron-resonant scattering by interplanetary magnetic field fluctuations are examined and shown to be negligible. An electromagnetic plasma instability, triggered by the birth of ions in the exosphere, is shown to be important for the thermalization of the energy mode transverse to the interplanetary magnetic field, allowing more ions to be detected by solar-wind ion probes.     

Negative ion composition and sulfuric acid vapour in the upper stratosphere

The nature of negative ions in the altitude region 42–45 km has been investigated by means of a balloon borne mass spectrometer. Apart from the NO₃⁻ and HSO₄⁻ clusters, ions with different cores, which can be identified as CO₃⁻, HCO₃⁻, Cl⁻ and ClO₃⁻ were observed. The spectra have been used to estimate the sulfuric acid number density at 45.2 and 42.3 km altitude.     

Large-scale magnetic fields: galaxy two-point correlation function

We study the effect of large-scale tangled magnetic fields on the galaxy two-point correlation function in the redshift space. We show that (i) the magnetic field effects can be comparable to the gravity-induced clustering for present magnetic field strength   B ₀≃ 5 × 10⁻⁸ G  , (ii) the absence of this signal from the present data gives an upper bound   B ₀≲ 3 × 10⁻⁸ G  and (iii) the future data can probe the magnetic fields of  ≃10⁻⁸ G  . A comparison with other constraints on the present magnetic field shows that they are marginally compatible. However, if the magnetic fields corresponding to   B ₀≃ 10⁻⁸ G  existed at the last scattering surface, they will cause unacceptably large cosmic microwave background radiation anisotropies.     

Emergence of magnetic field due to spin-polarized baryon matter in neutron stars

A model of the ferromagnetic origin of magnetic fields of neutron stars is considered. In this model, the magnetic phase transition occurs inside the core of neutron stars soon after formation. However, owing to the high electrical conductivity the core magnetic field is initially fully screened. We study how this magnetic field emerges for an outside observer. After some time, the induced field that screens the ferromagnetic field decays enough to uncover a detectable fraction of the ferromagnetic field. We calculate the time-scale of decay of the screening field and study how it depends on the size of the ferromagnetic core. We find that the same fractional decay of the screening field occurs earlier for larger cores. We conjecture that weak fields of millisecond pulsars, B ∼10⁸–10⁹ G, could be identified with ferromagnetic fields of unshielded fraction ε ∼10⁻⁴–10⁻³ resulting from the decay of screening fields by a factor 1− ε in ∼10⁸ yr since their birth.     

Three-dimensional simulations of molecular cloud fragmentation regulated by magnetic fields and ambipolar diffusion

We employ the first fully three-dimensional simulation to study the role of magnetic fields and ion–neutral friction in regulating gravitationally driven fragmentation of molecular clouds. The cores in an initially subcritical cloud develop gradually over an ambipolar diffusion time while the cores in an initially supercritical cloud develop in a dynamical time. The infalling speeds on to cores are subsonic in the case of an initially subcritical cloud, while an extended (≳0.1 pc) region of supersonic infall exists in the case of an initially supercritical cloud. These results are consistent with previous two-dimensional simulations. We also found that a snapshot of the relation between density (ρ) and the strength of the magnetic field ( B ) at different spatial points of the cloud coincides with the evolutionary track of an individual core. When the density becomes large, both the relations tend to   B ∝ρ^0.5  .     

The bottom magnetic field and magnetosphere evolution of neutron star in low-mass X-ray binary

The accretion-induced neutron star (NS) magnetic field evolution is studied through considering the accretion flow to drag the field lines aside and dilute the polar-field strength, and as a result the equatorial field strength increases, which is buried inside the crust on account of the accretion-induced global compression of star crust. The main conclusions of model are as follows: (i) the polar field decays with increase in the accreted mass; (ii) the bottom magnetic field strength of about 10⁸ G can occur when the NS magnetosphere radius approaches the star radius, and it depends on the accretion rate as     ; and (iii) the NS magnetosphere radius decreases with accretion until it reaches the star radius, and its evolution is little influenced by the initial field and the accretion rate after accreting  ∼0.01 M_⊙  , which implies that the magnetosphere radii of NSs in low-mass X-ray binaries would be homogeneous if they accreted the comparable masses. As an extension, the physical effects of the possible strong magnetic zone in the X-ray NSs and recycled pulsars are discussed. Moreover, the strong magnetic fields in the binary pulsars PSR 1831−00 and PSR 1718−19 after accreting about  0.5 M_⊙  in the binary-accretion phase,  8.7 × 10¹⁰  and  1.28 × 10¹² G  , respectively, can be explained through considering the incomplete frozen flow in the polar zone. As an expectation of the model, the existence of the low magnetic field  (∼3 × 10⁷ G)  NSs or millisecond pulsars is suggested.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.     

Superbubbles in magnetized interstellar media: blowout or confinement?

The importance of the interstellar magnetic field is studied in relation to the evolution of superbubbles with a three-dimensional (3D) numerical magnetohydrodynamical (MHD) simulation. A superbubble is a large supernova remnant driven by sequential supernova explosions in an OB association. Its evolution is affected by the density stratification in the galactic disc. After the superbubble size reaches 2–3 times the density scaleheight, it expands preferentially in the z -direction, until finally it can punch out a hole in the gas disc (blowout). On the other hand, the magnetic field running parallel to the galactic disc has the effect of preventing it from expanding in the direction perpendicular to the field. The density stratification and the magnetic fields have completely opposite effects on the evolution of the superbubble. We present results of a 3D MHD simulation in which both effects are included. As a result, it is concluded that when the magnetic field has a much larger scaleheight than the density, even for a model in which the bubble would blow out from the disc if the magnetic field were absent, a magnetic field with a strength of 5 μG can confine the bubble in | z |≲300 pc for ≃ 20 Myr (confinement). In a model in which the field strength decreases in the halo as B  ∝ ρ^1/2, the superbubble eventually blows out like a model with B  = 0 even if the magnetic field in the mid-plane is as strong as B  = 5 μG.     

The configuration of the loop-like structure of the solar atmosphere

An integral, governing steady flows in an isolated thin magnetic flux tube in the hydrostatic plane-stratified atmosphere, has been obtained. The integral, that we named as the shape integral, is expressed as (1 − M_A²)B cos θ = const. Here M_A² is the Alfven Mach number, B is the magnetic field strength and θ is the flux tube inclination to the horizontal. The shape integral should hold for most loop models because it represents just the momentum balance laws and has no relation to any energy balance mode. Its application to the isothermal and static cases is discussed and illustrated.     

Magnetic fields in the 3C 129 cluster

We present multifrequency VLA observations of the two radio galaxies 3C 129 and 3C 129.1 embedded in a luminous X-ray cluster. These radio observations reveal a substantial difference in the Faraday Rotation Measures (RMs) toward 3C 129.1 at the cluster centre and 3C 129 at the cluster periphery. After deriving the density profile from available X-ray data, we find that the RM structure of both radio galaxies can be fitted by a tangled cluster magnetic field with strength 6 μG extending at least 3 core radii (450 kpc) from the cluster centre. The magnetic field makes up a small contribution to the total pressure (5 per cent) in the central regions of the cluster. The radio morphology of 3C 129.1 appears disturbed on the southern side, perhaps by the higher pressure environment. In contrast with earlier claims for the presence of a moderately strong cooling flow in the 3C 129 cluster, our analysis of the X-ray data places a limit on the mass deposition rate from any such flow of < 1.2 M_⊙ yr⁻¹.     

Magnetic hydrogen atmosphere models and the neutron star RX J1856.5–3754

RX J1856.5−3754 is one of the brightest nearby isolated neutron stars (INSs), and considerable observational resources have been devoted to it. However, current models are unable to satisfactorily explain the data. We show that our latest models of a thin, magnetic, partially ionized hydrogen atmosphere on top of a condensed surface can fit the entire spectrum, from X-rays to optical, of RX J1856.5−3754, within the uncertainties. In our simplest model, the best-fitting parameters are an interstellar column density   N _H≈ 1 × 10²⁰ cm⁻²  and an emitting area with   R ^∞≈ 17 km  (assuming a distance to RX J1856.5−3754 of 140 pc), temperature   T ^∞≈ 4.3 × 10⁵ K  , gravitational redshift   z _g ∼ 0.22  , atmospheric hydrogen column   y _H≈ 1 g cm⁻²  , and magnetic field   B ≈ (3–4) × 10¹² G  ; the values for the temperature and magnetic field indicate an effective average over the surface. We also calculate a more realistic model, which accounts for magnetic field and temperature variations over the NS surface as well as general relativistic effects, to determine pulsations; we find that there exist viewing geometries that produce pulsations near the currently observed limits. The origin of the thin atmospheres required to fit the data is an important question, and we briefly discuss mechanisms for producing these atmospheres. Our model thus represents the most self-consistent picture to date for explaining all the observations of RX J1856.5−3754.     

Hydrostatic equilibrium in a magnetized, warped Galactic disc

Hydrostatic equilibrium of the multiphase interstellar medium in the solar vicinity is reconsidered, with the regular and turbulent magnetic fields treated separately. The regular magnetic field strength required to support the gas is consistent with independent estimates, provided that energy equipartition is maintained between turbulence and random magnetic fields. Our results indicate that a mid-plane value of B ₀=4 μG for the regular magnetic field near the Sun leads to more attractive models than B ₀=2 μG . The vertical profiles of both the regular and random magnetic fields contain disc and halo components, the parameters of which we have determined. The layer at 1≲| z |≲4 kpc can be overpressured and an outflow at a speed of about 50 km s⁻¹ may occur there, presumably associated with a Galactic fountain flow, if B ₀≃2 μG .
We show that hydrostatic equilibrium in a warped disc must produce asymmetric density distributions in z , in rough agreement with H  i observations in the outer Galaxy. This asymmetry may be a useful diagnostic of the details of the warping mechanism in the Milky Way and other galaxies. We find indications that gas and magnetic field pressures are different above and below the warped midplane in the outer Galaxy, and quantify the difference in terms of turbulent velocity and/or magnetic field strength.     

Fields and filaments in the core of the Centaurus cluster

We present high-resolution images of the Faraday rotation measure (RM) structure of the radio galaxy PKS 1246−410 at the centre of the Centaurus cluster. Comparison with Hα-line and soft X-ray emission reveals a correspondence between the line-emitting gas, the soft X-ray emitting gas, regions with an excess in the RM images and signs of depolarization. Magnetic field strengths of 25 μG, organized on scales of ∼1 kpc and intermixed with gas at a temperature of 5 × 10⁶ K with a density of ∼0.1 cm⁻³, can reproduce the observed RM excess, the depolarization and the observed X-ray surface brightness. This hot gas may be in pressure equilibrium with the optical line-emitting gas, but the magnetic field strength of 25 μG associated with the hot gas provides only 10 per cent of the thermal pressure and is therefore insufficient to account for the stability of the line-emitting filaments.     

Physical and dynamical properties of a low-mass star-forming region

Eleven low-mass cores are found in the Orion Molecular Cloud 2 from VLA observations of the line emission of NH₃ (1,1). They are perhaps clumps prior to gravitatonal collapse with average radius of 0.03 pc and mass of 3.5 M, distributed along the central axes of filaments in the NS direction. We find a velocity gradient of 5 km s⁻¹ pc⁻¹ in the declination direction within a 3′ region. Based on our NH₃(1,1) observations and compared with dust continuum emission at millimilion wavelength as well as in the infrared, we suggest that most of these dense cores are probably protosteller condensations, not yet containing stellar cores, but are self-gravitating systems in thermodynamical equilibrium.     

Hot stars in old stellar populations: a continuing need for intermediate ages

We present numerical investigations into the formation of massive stars from turbulent cores of density structure  ρ∝ r ^−1.5  . The results of five hydrodynamical simulations are described, following the collapse of the core, fragmentation and the formation of small clusters of protostars. We generate two different initial turbulent velocity fields corresponding to power-law spectra   P ∝ k ⁻⁴  and   P ∝ k ^−3.5  , and we apply two different initial core radii. Calculations are included for both completely isothermal collapse, and a non-isothermal equation of state above a critical density  (10⁻¹⁴ g cm⁻³)  . Our calculations reveal the preference of fragmentation over monolithic star formation in turbulent cores. Fragmentation was prevalent in all the isothermal cases. Although disc fragmentation was largely suppressed in the non-isothermal runs due to the small dynamic range between the initial density and the critical density, our results show that some fragmentation still persisted. This is inconsistent with previous suggestions that turbulent cores result in the formation of a single massive star. We conclude that turbulence cannot be measured as an isotropic pressure term.     

Formation of interstellar clouds: Parker instability with phase transitions

We numerically follow the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold H  i interstellar medium in two spatial dimensions. The nonlinear evolution of the system favours modes that allow the magnetic field lines to cross the galactic plane. Cold H  i clouds form with typical masses  ≃10⁵ M_⊙  , mean densities  ≃20 cm⁻³  , mean magnetic-field strengths  ≃4.3 μG  (rms field strengths  ≃6.4 μG  ), mass-to-flux ratios  ≃0.1–0.3  relative to critical, temperatures  ≃50 K  , (two-dimensional) turbulent velocity dispersions  ≃1.6 km s⁻¹  and separations  ≃500 pc  , in agreement with observations. The maximum density and magnetic-field strength are  ≃10³ cm⁻³  and  ≃20 μG  , respectively. Approximately 60 per cent of all H  i mass is in the warm neutral medium. The cold neutral medium is arranged into sheet-like structures both perpendicular and parallel to the galactic plane, but it is also found almost everywhere in the galactic plane, with the density being highest in valleys of the magnetic field lines. 'Cloudlets' also form whose physical properties are in quantitative agreement with those observed for such objects by Heiles. The nonlinear phase of the evolution takes ≲30 Myr, so that, if the instability is triggered by a nonlinear perturbation such as a spiral density shock wave, interstellar clouds can form within a time suggested by observations.     

An Observational Study of the Physical Properties of Star-forming Molecular Cores

We used the 13.7 m millimeter-wave telescope of the Qinghai Station of Purple Mountain Observatory at Delingha to carry out simultaneous mapping observations in the ¹²CO(J = 1 − 0), ¹³CO(J = 1 − 0) and CO¹⁸ (J = 1 − 0) lines towards 11 isolated star-forming molecular cores selected from the sample of Spitzer's c2d program. For all three lines, these sources were all observed to the positions of half the peak intensity of the ¹³CO(J = 1 − 0) line. The volume density, local thermodynamic equilibrium mass M_LTE and virial mass M_VIR were derived for each of the molecular cores. The obtained ratio of M_VIR to M_LTE is 0.85 ± 0.40 for the ¹³CO(J = 1 − 0) cores and 0.77 ± 0.35 for the CO¹⁸ (J = 1 − 0) cores. The density profiles of the cloud cores were also calculated.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.