Function Of Column Density Research Articles

A halo model for cosmological Lyman-limit systems

ABSTRACT We present an analytical model for cosmological Lyman-limit systems (lls) that successfully reproduces the observed evolution of the mean free path (λeff) of ionizing photons. The evolution of the co-moving mean free path is predominantly a consequence of the changing meta galactic photoionization rate and the increase with cosmic time of the minimum mass below which haloes lose their gas due to photoheating. In the model, Lyman-limit absorption is caused by highly ionized gas in the outskirt of dark matter haloes. We exploit the association with haloes to compute statistical properties of λeff and of the bias, b, of lls. The latter increases from b ∼ 1.5 → 2.6 from redshifts z = 2 → 6. Combined with the rapid increase with redshift of the bias of the haloes that host a quasar, the model predicts a rapid drop in the value of λeff when measured in quasar spectra from z = 5 → 6, whereas the actual value of λeff falls more smoothly. We derive an expression for the effective optical depth due to Lyman limit absorption as a function of wavelength and show that it depends sensitively on the poorly constrained number density of lls as a function of column density. The optical depth drops below unity for all wavelengths below a redshift of ∼2.5 which is therefore the epoch when the Universe first became transparent to ionizing photons.

Precipitating Solar Wind Hydrogen as Observed by the MAVEN Spacecraft: Distribution as a Function of Column Density, Altitude, and Solar Zenith Angle

AbstractThe Mars Atmosphere and Volatile EvolutioN (MAVEN) mission observes precipitating solar wind hydrogen ENAs in the upper Martian ionosphere with various instruments onboard. These precipitating solar wind hydrogen atoms result from the direct interaction between the solar wind and Mars's extensive hydrogen corona, creating a beam of particles traveling at solar wind speed that can penetrate deep into the ionosphere of Mars. Using particle data collected by both the Solar Wind Ion Analyzer (SWIA) and Neutral Gas and Ion Mass Spectrometer (NGIMS) onboard MAVEN, we seek to characterize the global behavior of these particles. By using a combination of energy spectra and neutral gas density data, we are able to determine the behavior of this particle population as a function of altitude, column density, and solar zenith angle. Through our analysis, we see that these precipitating solar wind hydrogen ENAs undergo multiple charge exchange interactions and ultimately scatter significantly at an altitude of 160 km or column density of 5 × 1016 cm−2. This observed column density at which equilibrium is reached is at a value ∼17 times greater than the theoretical value. We also demonstrate that this particle population's energy flux follows a cosine dependence of solar zenith angle, as is predicted by models.

Dust polarization studies on MHD simulations of molecular clouds: comparison of methods for the relative-orientation analysis

Context. The all-sky survey from the Planck space telescope has revealed that thermal emission from Galactic dust is polarized on scales ranging from the whole sky down to the inner regions of molecular clouds. Polarized dust emission can therefore be used as a probe for magnetic fields on different scales. In particular, the analysis of the relative orientation between the density structures and the magnetic field projected on the plane of the sky can provide information on the role of magnetic fields in shaping the structure of molecular clouds where star formation takes place. Aims. The orientation of the magnetic field with respect to the density structures has been investigated using different methods. The goal of this paper is to explicitly compare two of these: the Rolling Hough Transform (RHT) and the gradient technique (GRAD). Methods. We generated synthetic surface brightness maps at 353 GHz (850 μm) via magnetohydrodynamic simulations. We applied RHT and GRAD to two morphologically different regions identified in our maps. Region 1 is dominated by a dense and thick filamentary structure with some branches, while Region 2 includes a thinner filament with denser knots immersed in a more tenuous medium. Both methods derive the relative orientation between the magnetic field and the density structures, to which we applied two statistics, the histogram of relative orientation and the projected Rayleigh statistic, to quantify the variations of the relative orientation as a function of column density. Results. Both methods find areas with significant signal, and these areas are substantially different. In terms of relative orientations, in all our considered cases the predominant orientation of the density structures is perpendicular to the direction of the magnetic field. When the methods are applied to the same selected areas the results are consistent with each other in Region 2 but show some noticeable differences in Region 1. In Region 1, RHT globally finds the relative orientation becoming more perpendicular for increasing column density, while GRAD, applied at the same resolution as RHT, gives the opposite trend. These disparities are caused by the intrinsic differences in the methods and in the structures that they select. Conclusions. Our results indicate that the interpretation of the relative orientation between the magnetic field and density structures should take into account the specificity of the methods used to determine such orientation. The combined use of complementary techniques such as RHT and GRAD provides more complete information, which can be advantageously used to better understand the physical mechanisms operating in magnetized molecular clouds.

Dust and polarization of cold clumps

AbstractThe Planck catalogue of Galactic cold clumps, PGCC, contains sources of ongoing and future star formation. The data show clear variations also in their dust properties.We use Planck polarization measurements to investigate the polarization fraction in PGCC clumps and the relative orientation of filamentary structures and magnetic fields (Alina et al. 2017). The decrease of polarization fraction as a function of column density can be related to the field geometry but also suggest some loss of grain alignment.PGCCs have been studied with ground-based observations (Liu et al. 2018). The first SCUBA-2/POL-2 polarization studies have targeted the infrared dark cloud G35.39-0.33. The magnetic field is found to be mostly perpendicular to the main filament. The plane-of-the-sky field strength is $\sim 50\mu \,{\rm{G}}$ , a noticeable support against gravity. The polarization fraction decreases with increasing column density. This matches predictions of RAT grain alignment models but the relative contribution of the field morphology is hard to quantify (Juvela et al. 2018).We continue to use MHD simulations to study the same phenomena, with synthetic observations of clumps and filaments.

Simulating the effect of high column density absorbers on the one-dimensional Lyman α forest flux power spectrum

We measure the effect of high column density absorbing systems of neutral hydrogen (HI) on the one-dimensional (1D) Lyman-alpha forest flux power spectrum using cosmological hydrodynamical simulations from the Illustris project. High column density absorbers (which we define to be those with HI column densities $N(\mathrm{HI}) > 1.6 \times 10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$) cause broadened absorption lines with characteristic damping wings. These damping wings bias the 1D Lyman-alpha forest flux power spectrum by causing absorption in quasar spectra away from the location of the absorber itself. We investigate the effect of high column density absorbers on the Lyman-alpha forest using hydrodynamical simulations for the first time. We provide templates as a function of column density and redshift, allowing the flexibility to accurately model residual contamination, i.e., if an analysis selectively clips out the largest damping wings. This flexibility will improve cosmological parameter estimation, e.g., allowing more accurate measurement of the shape of the power spectrum, with implications for cosmological models containing massive neutrinos or a running of the spectral index. We provide fitting functions to reproduce these results so that they can be incorporated straightforwardly into a data analysis pipeline.

The observable properties of cool winds from galaxies, AGN, and star clusters – I. Theoretical framework

Winds arising from galaxies, star clusters, and active galactic nuclei are crucial players in star and galaxy formation, but it has proven remarkably difficult to use observations of them to determine physical properties of interest, particularly mass fluxes. Much of the difficulty stems from a lack of a theory that links a physically-realistic model for winds' density, velocity, and covering factors to calculations of light emission and absorption. In this paper we provide such a model. We consider a wind launched from a turbulent region with a range of column densities, derive the differential acceleration of gas as a function of column density, and use this result to compute winds' absorption profiles, emission profiles, and emission intensity maps in both optically thin and optically thick species. The model is sufficiently simple that all required computations can be done analytically up to straightforward numerical integrals, rendering it suitable for the problem of deriving physical parameters by fitting models to observed data. We show that our model produces realistic absorption and emission profiles for some example cases, and argue that the most promising methods of deducing mass fluxes are based on combinations of absorption lines of different optical depths, or on combining absorption with measurements of molecular line emission. In the second paper in this series, we expand on these ideas by introducing a set of observational diagnostics that are significantly more robust that those commonly in use, and that can be used to obtain improved estimates of wind properties.

Galactic cold cores

The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the whole sample and inside individual fields. Based on IRAS and Planck data, the fields are characterised by a median colour temperature of 16.1 K and a median opacity spectral index of beta=1.84. We observe a clear T-beta anti-correlation. In Herschel observations, constrained at lower resolution by Planck data, the variations follow the column density structure and beta(FIR) can rise to ~2.2 in individual clumps. The Planck 217 GHz band shows a systematic excess that is consistent with a general flattening of the dust emission spectrum at millimetre wavelengths. When fitted separately below and above 700 um, the median spectral index values are beta(FIR) ~ 1.91 and beta(mm) ~ 1.66. The spectral index changes as a function of column density and wavelength. Beta variations are partly masked by temperature gradients and the changes in the intrinsic grain properties may be even greater.

Polarized thermal dust emission as seen by Planck : A comparison with MHD simulations and lessons from a toy model

AbstractThe Planck satellite has mapped the polarized microwave sky (from 30 GHz to 353 GHz) with unprecedented sensitivity and angular resolution. This wealth of data yields the first complete map of polarized thermal emission from dust in our own Galaxy, shedding new light on the formation of dense cold structures within which new stars and planetary systems are born, under the combined effects of gravity, turbulence and magnetic fields. We present a statistical analysis of this polarized emission from nearby molecular clouds, with an emphasis on the evolution of the maximum polarization fraction observed as a function of column density, and on the anti-correlation between the polarization fraction and the local dispersion of polarization angles. To interpret this data, numerical simulations of anisotropic MHD turbulence underline the essential role played by the topology of the interstellar magnetic field, in particular its large-scale component. Indeed, the polarization of dust thermal emission at the scales observed by Planck is essentially related to the geometry of the magnetic field. Polarization fractions anti-correlate with column densities, which may be due to a succession of variously polarized structures on the line of sight. They also anti-correlate with the local dispersion of polarization angles. These features are well reproduced by MHD simulations of the diffuse ISM, with comparable correlation coefficients. As an extension to this work published in Planck Intermediate Results XX (A&A, 576, 105, 2015), the statistical properties of the random component of the interstellar magnetic field are explored using a toy model of the turbulent magnetized ISM based on fractional Brownian motion (fBm) fields. A least-squares analysis to retrieve the statistical properties of the interstellar magnetic field from Planck observations is pursued. Application of this method on the toy model shows good promise, and we are currently working towards its application on Planck data.

Polarized thermal dust emission as seen by Planck : A comparison with MHD simulations and lessons from a toy model

AbstractThe Planck satellite (Planck 2015 Results I) has mapped the polarized microwave sky (from 30 GHz to 353 GHz) with unprecedented sensitivity and angular resolution. This wealth of data yields the first complete map of polarized thermal emission from dust in our own Galaxy (Planck Intermediate Results XIX),, shedding new light on the formation of dense cold structures within which new stars and planetary systems are born, under the combined effects of gravity, turbulence and magnetic fields. We present a statistical analysis of this polarized emission from nearby molecular clouds, with an emphasis on the evolution of the maximum polarization fraction observed as a function of column density, and on the anti-correlation between the polarization fraction and the local dispersion of polarization angles. To interpret this data, numerical simulations of anisotropic MHD turbulence (Fromang, Hennebelle, & Teyssier 2006, Hennebelle et al. 2008) underline the essential role played by the topology of the interstellar magnetic field, in particular its large-scale component (Planck Intermediate Results XX). Indeed, the polarization of dust thermal emission at the scales observed by Planck is essentially related to the geometry of the magnetic field. Polarization fractions anti-correlate with column densities, which may be due to a succession of variously polarized structures on the line of sight. They also anti-correlate with the local dispersion of polarization angles. These features are well reproduced by MHD simulations of the diffuse ISM, with comparable correlation coefficients. As an extension to this work published in Planck Intermediate Results XX, the statistical properties of the random component of the interstellar magnetic field are explored using a toy model of the turbulent magnetized ISM based on fractional Brownian motion (fBm) fields. A least-squares analysis to retrieve the statistical properties of the interstellar magnetic field from Planck observations is pursued. Application of this method on the toy model shows good promise, and we are currently working towards its application on Planck data.

Polarized thermal dust emission as seen by Planck : A comparison with MHD simulations and lessons from a toy model

AbstractThe Planck satellite has mapped the polarized microwave sky (from 30 GHz to 353 GHz) with unprecedented sensitivity and angular resolution. This wealth of data yields the first complete map of polarized thermal emission from dust in our own Galaxy, shedding new light on the formation of dense cold structures within which new stars and planetary systems are born, under the combined effects of gravity, turbulence and magnetic fields. We present a statistical analysis of this polarized emission from nearby molecular clouds, with an emphasis on the evolution of the maximum polarization fraction observed as a function of column density, and on the anti-correlation between the polarization fraction and the local dispersion of polarization angles. To interpret this data, numerical simulations of anisotropic MHD turbulence underline the essential role played by the topology of the interstellar magnetic field, in particular its large-scale component. As an extension to this work published in Planck Intermediate Results XX (A&A, 576, 105, 2015), the statistical properties of the random component of the interstellar magnetic field are explored using a toy model based on fractional Brownian motion (fBm) fields.

What probability distribution functions tell us about the processes of star formation

Probability distribution functions of column density (N-PDFs) are used to evaluate the relative importance of gravity, turbulence, magnetic fields, geometry, and radiative feedback governing the cloud's density structure and star-formation activity. These proce- sses influence the N-PDF and thus determine their shape, ranging from a purely lognormal distribution to one with power-law tail(s).

Moving-mesh cosmology: properties of neutral hydrogen in absorption

We examine the distribution of neutral hydrogen in cosmological simulations carried out with the new moving-mesh code AREPO and compare it with the corresponding GADGET simulations based on the smoothed particle hydrodynamics (SPH) technique. The two codes use identical gravity solvers and baryonic physics implementations, but very different methods for solving the Euler equations, allowing us to assess how numerical effects associated with the hydro-solver impact the results of simulations. Here we focus on an analysis of the neutral gas, as detected in quasar absorption lines. We find that the high column density regime probed by Damped Lyman-alpha (DLA) and Lyman Limit Systems (LLS) exhibits significant differences between the codes. GADGET produces spurious artefacts in large halos in the form of gaseous clumps, boosting the LLS cross-section. Furthermore, it forms halos with denser central baryonic cores than AREPO, which leads to a substantially greater DLA cross-section from smaller halos. AREPO thus produces a significantly lower cumulative abundance of DLAs, which is intriguingly in much closer agreement with observations. The column density function, however, is not altered enough to significantly reduce the discrepancy with the observed value. For the low column density gas probed by the Lyman-alpha forest, the codes differ only at the level of a few percent, suggesting that this regime is quite well described by both methods, a fact that is reassuring for the many Lyman-alpha studies carried out with SPH thus far. While the residual differences are smaller than the errors on current Lyman-alpha forest data, we note that this will likely change for future precision experiments.

THE LOW-REDSHIFT INTERGALACTIC MEDIUM AS SEEN IN ARCHIVAL LEGACYHST/STIS ANDFUSEDATA

We present a comprehensive catalog of ultraviolet HST/STIS and FUSE absorbers in the low-redshift IGM at z<0.4. The catalog draws from the extensive literature on IGM absorption, and it reconciles discrepancies among previous catalogs through a critical evaluation of all reported absorption features in light of new HST/COS data. We report on 746 HI absorbers down to a rest-frame equivalent width of 12 milliAngstroms over a maximum redshift path length Deltaz=5.38. We also confirm 111 OVI absorbers, 29 CIV absorbers, and numerous absorption features due to other metal ions. We characterize the distribution of absorber line frequency as a function of column density as a power law, dN/dz \propto N^{-beta}, where beta=2.08+-0.12 for OVI and beta=1.68+-0.03 for HI. Utilizing a more sophisticated accounting technique than past work, the catalog accounts for ~43% of the baryons: 24+-2% in the photoionized Ly-alpha forest and 19+-2% in the WHIM as traced by OVI. We discuss the large systematic effects of various assumed metallicities and ionization states on these calculations, and we implement recent simulation results in our estimates.

Feedback of Active Galactic Nuclei in Seyfert 2 Galaxies

It is well accepted that feedback from active galactic nuclei (AGNs) plays an important role in the coevolution of the supermassive black hole (SMBH) and its host galaxy, but the concrete mechanism of feedback remains unclear. A considerable body of evidence suggests that AGN feedback suppresses star formation in the host galaxy. We assemble a sample of Seyfert 2 galaxies with recent observational data of compact nuclear starbursts and estimate the gas surface density as a function of column density to illuminate the relation between feedback and AGN properties. Although there are some uncertainties, our data still imply the deviation from the star formation law (Kennicutt-Schmidt law). Further, they indicate that: (1) Feedback correlates with the Eddington ratio, rather than with the mass of SMBH, as a result of decreasing star formation efficiency. (2) The SMBH and the torus are probably undergoing coevolution. Conclusions presented here can be refined through future high resolution CO or HCN observations.

The Origin of Diffuse X‐Ray and γ‐Ray Emission from the Galactic Center Region: Cosmic‐Ray Particles

The inner couple of hundred parsecs of our Galaxy are characterized by a significant amount of synchrotron-emitting gas. Many of best studied sources in this region exhibit a mixture of 6.4 keV Fe Kα emission, molecular line emission, and nonthermal radio continuum radiation. The spatial correlation between fluorescent Fe Kα line emission at 6.4 keV and molecular line emission from Galactic center molecular clouds has been explained as reflected X-rays from a past outburst of Sgr A*. Here we present a multiwavelength study of this region and find a correlation between nonthermal radio filaments and X-ray features. This correlation, when combined with distribution of molecular gas, suggests against irradiation model. Instead, we account for this distribution in terms of impact of relativistic particles from local (nonthermal filaments) and extended sources with diffuse neutral gas producing both nonthermal bremsstrahlung X-ray continuum emission and diffuse 6.4 keV line emission. The production rate of Fe Kα photons associated with injection of electrons into a cloud as a function of column density is calculated. The required energy density of low-energy cosmic rays associated with synchrotron-emitting radio filaments or extended features is estimated to be in range between 20 and ~103 eV cm-3 for Sgr C, Sgr B1, Sgr B2, and the 45 and -30 km s-1 clouds. We also generalize this idea to explain cosmic-ray heating of molecular gas, interstellar cosmic-ray ionization, pervasive production of diffuse Kα line, and TeV emission from Galactic center molecular clouds. In particular, we suggest that inverse Compton scattering of submillimeter radiation from dust by relativistic electrons may contribute substantially to large-scale diffuse TeV emission observed toward central regions of Galaxy.

Physical effects on the Lyα forest flux power spectrum: damping wings, ionizing radiation fluctuations and galactic winds

We explore several physical effects on the power spectrum of the Lyα forest transmitted flux. The effects we investigate here are not usually part of hydrodynamic simulations and so need to be estimated separately. The most important effect is that of high column density absorbers with damping wings, which add power on large scales. We compute their effect using the observational constraints on their abundance as a function of column density. Ignoring their effect leads to an underestimation of the slope of the linear theory power spectrum. The second effect we investigate is that of fluctuations in the ionizing radiation field. For this purpose we use a very large high-resolution N-body simulation, which allows us to simulate both the fluctuations in the ionizing radiation and the small-scale Lyα forest within the same simulation. We find an enhancement of power on large scales for quasars and a suppression for galaxies. The strength of the effect rapidly increases with increasing redshift, allowing it to be uniquely identified in cases where it is significant. We develop templates that can be used to search for this effect as a function of quasar lifetime, quasar luminosity function and attenuation length. Finally, we explore the effects of galactic winds using hydrodynamic simulations. We find the wind effects on the Lyα forest power spectrum to be degenerate with parameters related to the temperature of the gas that are already marginalized over in cosmological fits. While more work is needed to conclusively exclude all possible systematic errors, our results suggest that, in the context of data analysis procedures, where parameters of the Lyα forest model are properly marginalized over, the flux power spectrum is a reliable tracer of cosmological information. Ke yw ords: galaxies: high-redshift ‐ intergalactic medium ‐ quasars: absorption lines ‐ cosmology: theory ‐ diffuse radiation ‐ large-scale structure of Universe.

The thermal history of the intergalactic medium

At redshifts z ≳ 2, most of the baryons reside in the smooth intergalactic medium which is responsible for the low column density Lyα forest. This photoheated gas follows a tight temperature-density relation which introduces a cut-off in the distribution of widths of the Lyα absorption lines (b-parameters) as a function of column density. We have measured this cut-off in a sample of nine high-resolution, high signal-to-noise ratio quasar spectra and determined the thermal evolution of the intergalactic medium in the redshift range 2.0–4.5. At a redshift z∼3, the temperature at the mean density shows a peak and the gas becomes nearly isothermal. We interpret this as evidence for the reionization of He ii.

Measuring the equation of state of the intergalactic medium

(Abridged) Numerical simulations indicate that the intergalactic medium (IGM) responsible for the low column density Lyman-alpha forest follows a well defined temperature-density relation. We demonstrate that such an equation of state results in a cutoff in the distribution of line widths (b-parameters) as a function of column density (N) for the low column density absorption lines. This explains the existence of the lower envelope which is clearly seen in scatter plots of the b(N)-distribution in observed QSO spectra. We show that the parameters of the cutoff in the b(N)-distribution are strongly correlated with the parameters of the underlying equation of state. We use simulations to determine these relations, which can then be applied to the observed cutoff in the b(N)-distribution to measure the equation of state of the IGM. Using Monte Carlo simulations of Keck spectra at z=3, we show that determining the slope of the equation of state will be difficult, but that the amplitude can be determined to within ten per cent, even from a single QSO spectrum. Measuring the evolution of the equation of state with redshift will allow us to put tight constraints on the reionization history of the universe.

On the role of nonthermal electrons in EUV and X-ray line emissions from solar flares

The energy and angular distribution of electrons as a function of column densities initially for monoenergetic and monodirectional electron beams and incidence angles of 0‡, 30‡ and 60‡ have been studied by combining small angle scattering using analytical treatment with large angle collisions using Monte Carlo calculations. Using these distributions, X-ray and EUV-line flux have been studied as a function of column density. It is observed that the line flux increases with the increase in column density, becoming significant at intermediate column densities where the electron energies and angular distributions have a non-Maxwellian nature.

On the role of reverse current on the hard X-ray production in solar flares and time-lags between high- and low-energy photons

The evolution of energy and angular distributions of electrons has been studied accounting for the reverse current effect by combining analytically treated small angle multiple scatterings with large angle Monte-Carlo calculations. Reverse current and potential variations as function of column density have been computed. It is found that the reverse current decreases steeply with increase in electron energy. However, it becomes significant for low-energy electrons. By use of these distributions and bremsstrahlung crosssection, the X-ray energy spectrum has been calculated. The nature of the resulting X-ray spectrum integrated over all column depths is similar to the one without reverse current. The time-lag between high-and low-energy photon production has been calculated. It is found that there is a small difference between time-lags as function of observation angles. This fact can be used to test the validity of the beamed thick target model.