Multifluid, Magnetohydrodynamic Shock Waves with Grain Dynamics II. Dust and the Critical Speed for C Shocks

This is the second in a series of papers on the effects of dust on multifluid, MHD shock waves in weakly ionized molecular gas. We investigate the influence of dust on the critical shock speed, v_crit, above which C shocks cease to exist. Chernoff showed that v_crit cannot exceed the grain magnetosound speed, v_gms, if dust grains are dynamically well coupled to the magnetic field. We present numerical simulations of steady shocks where the grains may be well- or poorly coupled to the field. We use a time-dependent, multifluid MHD code that models the plasma as a system of interacting fluids: neutral particles, ions, electrons, and various ``dust fluids'' comprised of grains with different sizes and charges. Our simulations include grain inertia and grain charge fluctuations but to highlight the essential physics we assume adiabatic flow, single-size grains, and neglect the effects of chemistry. We show that the existence of a phase speed v_phi does not necessarily mean that C shocks will form for all shock speeds v_s less than v_phi. When the grains are weakly coupled to the field, steady, adiabatic shocks resemble shocks with no dust: the transition to J type flow occurs at v_crit = 2.76 v_nA, where v_nA is the neutral Alfven speed, and steady shocks with v_s > 2.76 v_nA are J shocks with magnetic precursors in the ion-electron fluid. When the grains are strongly coupled to the field, v_crit = min(2.76 v_nA, v_gms). Shocks with v_crit < v_s < v_gms have magnetic precursors in the ion-electron-dust fluid. Shocks with v_s > v_gms have no magnetic precursor in any fluid. We present time-dependent calculations to study the formation of steady multifluid shocks. The dynamics differ qualitatively depending on whether or not the grains and field are well coupled.Comment: 43 pages with 17 figures, aastex, accepted by The Astrophysical Journa

Complex Langevin Equation and the Many-Fermion Problem

We study the utility of a complex Langevin (CL) equation as an alternative for the Monte Carlo (MC) procedure in the evaluation of expectation values occurring in fermionic many-body problems. We find that a CL approach is natural in cases where non-positive definite probability measures occur, and remains accurate even when the corresponding MC calculation develops a severe ``sign problem''. While the convergence of CL averages cannot be guaranteed in principle, we show how convergent results can be obtained in three examples ranging from simple one-dimensional integrals over quantum mechanical models to a schematic shell model path integral.Comment: 19 pages, 10 PS figures embedded in tex

Discovery of interstellar mercapto radicals (SH) with the GREAT instrument on SOFIA

We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unassociated with W49N, that is located along the sight-line. For the foreground material at velocities in the range 37 - 44 km/s with respect to the local standard of rest, we infer a total SH column density ~ 2.6 E+12 cm-2, corresponding to an abundance of ~ 7 E-9 relative to H2, and yielding an SH/H2S abundance ratio ~ 0.13. The observed SH/H2S abundance ratio is much smaller than that predicted by standard models for the production of SH and H2S in turbulent dissipation regions and shocks, and suggests that the endothermic neutral-neutral reaction SH + H2 -> H2S + H must be enhanced along with the ion-neutral reactions believed to produce CH+ and SH+ in diffuse molecular clouds.Comment: Accepted for publication in Astronomy and Astrophysics (SOFIA/GREAT special issue

Low, Milky-Way like, Molecular Gas Excitation of Massive Disk Galaxies at z~1.5

We present evidence for Milky-Way-like, low-excitation molecular gas reservoirs in near-IR selected massive galaxies at z~1.5, based on IRAM Plateau de Bure Interferometer CO[3-2] and NRAO Very Large Array CO[1-0] line observations for two galaxies that had been previously detected in CO[2-1] emission. The CO[3-2] flux of BzK-21000 at z=1.522 is comparable within the errors to its CO[2-1] flux, implying that the CO[3-2] transition is significantly sub-thermally excited. The combined CO[1-0] observations of the two sources result in a detection at the 3 sigma level that is consistent with a higher CO[1-0] luminosity than that of CO[2-1]. Contrary to what is observed in submillimeter galaxies and QSOs, in which the CO transitions are thermally excited up to J>=3, these galaxies have low-excitation molecular gas, similar to that in the Milky Way and local spirals. This is the first time that such conditions have been observed at high redshift. A Large Velocity Gradient analysis suggests that molecular clouds with density and kinetic temperature comparable to local spirals can reproduce our observations. The similarity in the CO excitation properties suggests that a high, Milky-Way-like, CO to H_2 conversion factor could be appropriate for these systems. If such low-excitation properties are representative of ordinary galaxies at high redshift, centimeter telescopes such as the Expanded Very Large Array and the longest wavelength Atacama Large Millimeter Array bands will be the best tools for studying the molecular gas content in these systems through the observations of CO emission lines.Comment: 5 pages, 4 figures. ApJ Letters in pres

Initial Ionization of Compressible Turbulence

We study the effects of the initial conditions of turbulent molecular clouds on the ionization structure in newly formed H_{ii} regions, using three-dimensional, photon-conserving radiative transfer in a pre-computed density field from three-dimensional compressible turbulence. Our results show that the initial density structure of the gas cloud can play an important role in the resulting structure of the H_{ii} region. The propagation of the ionization fronts, the shape of the resulting H_{ii} region, and the total mass ionized depend on the properties of the turbulent density field. Cuts through the ionized regions generally show ``butterfly'' shapes rather than spherical ones, while emission measure maps are more spherical if the turbulence is driven on scales small compared to the size of the H_{ii} region. The ionization structure can be described by an effective clumping factor $\zeta=< n > \cdot /^2$, where $n$ is number density of the gas. The larger the value of $\zeta$, the less mass is ionized, and the more irregular the H_{ii} region shapes. Because we do not follow dynamics, our results apply only to the early stage of ionization when the speed of the ionization fronts remains much larger than the sound speed of the ionized gas, or Alfv\'en speed in magnetized clouds if it is larger, so that the dynamical effects can be negligible.Comment: 9 pages, 10 figures, version with high quality color images can be found in http://research.amnh.org/~yuexing/astro-ph/0407249.pd

Detection of the 13CO(J=6-5) Transition in the Starburst Galaxy NGC 253

We report the detection of 13CO(J=6-5) emission from the nucleus of the starburst galaxy NGC 253 with the redshift (z) and Early Universe Spectrometer (ZEUS), a new submillimeter grating spectrometer. This is the first extragalactic detection of the 13CO(J=6-5) transition, which traces warm, dense molecular gas. We employ a multi-line LVG analysis and find ~ 35% - 60% of the molecular ISM is both warm (T ~ 110 K) and dense (n(H2) ~ 10^4 cm^-3). We analyze the potential heat sources, and conclude that UV and X-ray photons are unlikely to be energetically important. Instead, the molecular gas is most likely heated by an elevated density of cosmic rays or by the decay of supersonic turbulence through shocks. If the cosmic rays and turbulence are created by stellar feedback within the starburst, then our analysis suggests the starburst may be self-limiting.Comment: 4 pages, 2 figures, accepted by ApJ Letter

Epigenome-wide association study reveals decreased average methylation levels years before breast cancer diagnosis

Interest in the potential of DNA methylation in peripheral blood as a biomarker of cancer risk is increasing. We aimed to assess whether epigenome-wide DNA methylation measured in peripheral blood samples obtained before onset of the disease is associated with increased risk of breast cancer. We report on three independent prospective nested case-control studies from the European Prospective Investigation into Cancer and Nutrition (EPIC-Italy; n = 162 matched case-control pairs), the Norwegian Women and Cancer study (NOWAC; n = 168 matched pairs), and the Breakthrough Generations Study (BGS; n = 548 matched pairs). We used the Illumina 450k array to measure methylation in the EPIC and NOWAC cohorts. Whole-genome bisulphite sequencing (WGBS) was performed on the BGS cohort using pooled DNA samples, combined to reach 50× coverage across ~16 million CpG sites in the genome including 450k array CpG sites. Mean β values over all probes were calculated as a measurement for epigenome-wide methylation

Transcriptional implications of intragenic DNA methylation in the oestrogen receptor alpha gene in breast cancer cells and tissues

Background DNA methylation variability regions (MVRs) across the oestrogen receptor alpha (ESR1) gene have been identified in peripheral blood cells from breast cancer patients and healthy individuals. In contrast to promoter methylation, gene body methylation may be important in maintaining active transcription. This study aimed to assess MVRs in ESR1 in breast cancer cell lines, tumour biopsies and exfoliated epithelial cells from expressed breast milk (EBM), to determine their significance for ESR1 transcription. Methods DNA methylation levels in eight MVRs across ESR1 were assessed by pyrosequencing bisulphite-converted DNA from three oestrogen receptor (ER)-positive and three ER-negative breast cancer cell lines. DNA methylation and expression were assessed following treatment with DAC (1 μM), or DMSO (controls). ESR1 methylation levels were also assayed in DNA from 155 invasive ductal carcinoma biopsies provided by the Breast Cancer Campaign Tissue Bank, and validated with DNA methylation profiles from the TCGA breast tumours (n = 356 ER-pos, n = 109 ER-neg). DNA methylation was profiled in exfoliated breast epithelial cells from EBM using the Illumina 450 K (n = 36) and pyrosequencing in a further 53 donor samples. ESR1 mRNA levels were measured by qRT-PCR. Results We show that ER-positive cell lines had unmethylated ESR1 promoter regions and highly methylated intragenic regions (median, 80.45%) while ER-negative cells had methylated promoters and lower intragenic methylation levels (median, 38.62%). DAC treatment increased ESR1 expression in ER-negative cells, but significantly reduced methylation and expression of ESR1 in ER-positive cells. The ESR1 promoter was unmethylated in breast tumour biopsies with high levels of intragenic methylation, independent of ER status. However, ESR1 methylation in the strongly ER-positive EBM DNA samples were very similar to ER-positive tumour cell lines. Conclusion DAC treatment inhibited ESR1 transcription in cells with an unmethylated ESR1 promoter and reduced intragenic DNA methylation. Intragenic methylation levels correlated with ESR1 expression in homogenous cell populations (cell lines and exfoliated primary breast epithelial cells), but not in heterogeneous tumour biopsies, highlighting the significant differences between the in vivo tumour microenvironment and individual homogenous cell types. These findings emphasise the need for care when choosing material for epigenetic research and highlights the presence of aberrant intragenic methylation levels in tumour tissue

The Herschel HIFI water line survey in the low-mass proto-stellar outflow L1448

As part of the WISH (Water In Star-forming regions with Herschel) key project, we report on the observations of several ortho- and para-H2O lines performed with the HIFI instrument towards two bright shock spots (R4 and B2) along the outflow driven by the L1448 low-mass proto-stellar system, located in the Perseus cloud. These data are used to identify the physical conditions giving rise to the H2O emission and infer any dependence with velocity. These observations provide evidence that the observed water lines probe a warm (T_kin~400-600 K) and very dense (n 10^6 - 10^7 cm^-3) gas, not traced by other molecules, such as low-J CO and SiO, but rather traced by mid-IR H2 emission. In particular, H2O shows strong differences with SiO in the excitation conditions and in the line profiles in the two observed shocked positions, pointing to chemical variations across the various velocity regimes and chemical evolution in the different shock spots. Physical and kinematical differences can be seen at the two shocked positions. At the R4 position, two velocity components with different excitation can be distinguished, with the component at higher velocity (R4-HV) being less extended and less dense than the low velocity component (R4-LV). H2O column densities of about 2 10^13 and 4 10^14 cm^-2 have been derived for the R4-LV and the R4-HV components, respectively. The conditions inferred for the B2 position are similar to those of the R4-HV component, with H2O column density in the range 10^14 - 5 10^14 cm^-2, corresponding to H2O/H2 abundances in the range 0.5 - 1 10^-5. The observed line ratios and the derived physical conditions seem to be more consistent with excitation in a low velocity J-type shock with large compression rather than in a stationary C-shock, although none of these stationary models seems able to reproduce all the characteristics of the observed emission.Comment: Accepted for publication in A&

The gas temperature in flaring disks around pre-main sequence stars

A model is presented which calculates the gas temperature and chemistry in the surface layers of flaring circumstellar disks using a code developed for photon-dominated regions. Special attention is given to the influence of dust settling. It is found that the gas temperature exceeds the dust temperature by up to several hundreds of Kelvins in the part of the disk that is optically thin to ultraviolet radiation, indicating that the common assumption that Tgas=Tdust is not valid throughout the disk. In the optically thick part, gas and dust are strongly coupled and the gas temperature equals the dust temperature. Dust settling has little effect on the chemistry in the disk, but increases the amount of hot gas deeper in the disk. The effects of the higher gas temperature on several emission lines arising in the surface layer are examined. The higher gas temperatures increase the intensities of molecular and fine-structure lines by up to an order of magnitude, and can also have an important effect on the line shapes.Comment: 14 pages, 10 figures, accepted for publication in A&