Showing papers on "Radiative transfer published in 2011"

Physics of the Interstellar and Intergalactic Medium

Abstract: This is a comprehensive and richly illustrated textbook on the astrophysics of the interstellar and intergalactic medium--the gas and dust, as well as the electromagnetic radiation, cosmic rays, and magnetic and gravitational fields, present between the stars in a galaxy and also between galaxies themselves. Topics include radiative processes across the electromagnetic spectrum; radiative transfer; ionization; heating and cooling; astrochemistry; interstellar dust; fluid dynamics, including ionization fronts and shock waves; cosmic rays; distribution and evolution of the interstellar medium; and star formation. While it is assumed that the reader has a background in undergraduate-level physics, including some prior exposure to atomic and molecular physics, statistical mechanics, and electromagnetism, the first six chapters of the book include a review of the basic physics that is used in later chapters. This graduate-level textbook includes references for further reading, and serves as an invaluable resource for working astrophysicists. * Essential textbook on the physics of the interstellar and intergalactic medium * Based on a course taught by the author for more than twenty years at Princeton University * Covers radiative processes, fluid dynamics, cosmic rays, astrochemistry, interstellar dust, and more * Discusses the physical state and distribution of the ionized, atomic, and molecular phases of the interstellar medium * Reviews diagnostics using emission and absorption lines * Features color illustrations and detailed reference materials in appendices * Instructor's manual with problems and solutions (available only to teachers)

Finite difference solution of mhd radiative boundary layer flow of a nanofluid past a stretching sheet

Abstract: Unsteady heat and mass flow of a nanofluid past a stretching sheet with thermal radiation in the presence of magnetic field is studied. To obtain non-similar equation, continuity, momentum, energy and concentration equations have been non-dimensionalised by usual transformation. The non-similar solutions are presented here which depends on the Magnetic parameter M respectively . The obtained equations have been solved by explicit finite difference method. The temperature and concentration profiles are discussed for the different values of the above parameters with different time steps.

The stellar atmosphere simulation code Bifrost

Abstract: Context: Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very different physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the different layers interact strongly. These physical regimes are very diverse and it takes a highly efficient massively parallel numerical code to solve the associated equations. Aims: The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods: The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results: Bifrost completes the tests with good results and shows near linear efficiency scaling to thousands of computing cores.

The stellar atmosphere simulation code Bifrost - Code description and validation

Abstract: Context. Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very di erent physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the di erent layers interact strongly. These physical regimes are very diverse and it takes a highly e cient massively parallel numerical code to solve the associated equations. Aims. The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods. The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results. Bifrost completes the tests with good results and shows near linear e ciency scaling to thousands of computing cores.

A scaling theory for the size distribution of emitted dust aerosols suggests climate models underestimate the size of the global dust cycle.

Abstract: Mineral dust aerosols impact Earth’s radiation budget through interactions with clouds, ecosystems, and radiation, which constitutes a substantial uncertainty in understanding past and predicting future climate changes. One of the causes of this large uncertainty is that the size distribution of emitted dust aerosols is poorly understood. The present study shows that regional and global circulation models (GCMs) overestimate the emitted fraction of clay aerosols (< 2 μm diameter) by a factor of ∼2–8 relative to measurements. This discrepancy is resolved by deriving a simple theoretical expression of the emitted dust size distribution that is in excellent agreement with measurements. This expression is based on the physics of the scale-invariant fragmentation of brittle materials, which is shown to be applicable to dust emission. Because clay aerosols produce a strong radiative cooling, the overestimation of the clay fraction causes GCMs to also overestimate the radiative cooling of a given quantity of emitted dust. On local and regional scales, this affects the magnitude and possibly the sign of the dust radiative forcing, with implications for numerical weather forecasting and regional climate predictions in dusty regions. On a global scale, the dust cycle in most GCMs is tuned to match radiative measurements, such that the overestimation of the radiative cooling of a given quantity of emitted dust has likely caused GCMs to underestimate the global dust emission rate. This implies that the deposition flux of dust and its fertilizing effects on ecosystems may be substantially larger than thought.

An XMM–Newton spectral survey of 12 μm selected galaxies – I. X-ray data

Abstract: We present an X-ray spectral analysis of 126 galaxies of the 12 μm galaxy sample. By studying this sample at X-ray wavelengths, we aim to determine the intrinsic power, continuum shape and obscuration level in these sources. We improve upon previous works by the use of superior data in the form of higher signal-to-noise ratio spectra, finer spectral resolution and a broader bandpass from XMM–Newton. We pay particular attention to Compton thick active galactic nucleus (AGN) with the help of new spectral fitting models that we have produced, which are based on Monte Carlo simulations of X-ray radiative transfer, using both a spherical and torus geometry, and taking into account Compton scattering and iron fluorescence. We use this data to show that with a torus geometry, unobscured sightlines can achieve a maximum equivalent width of the Fe Kα line of ∼150 eV, originally shown by Ghisellini et al. In order for this to be exceeded, the line of sight must be obscured with NH > 10 23 cm −2 , as we show for one case, NGC 3690. We also calculate flux suppression factors from the simulated data, the main conclusion from which is that for NH ≥ 10 25 cm −2 , the X-ray flux is suppressed by a factor of at least 10 in all X-ray bands and at all redshifts, revealing the biases present against these extremely heavily obscured systems inherent in all X-ray surveys. Furthermore, we confirm previous results from Murphy & Yaqoob that show that the reflection fraction determined from slab geometries is underestimated with respect to toroidal geometries. For the 12 μm selected galaxies, we investigate the distribution of X-ray power-law indices, finding that the mean (� � �= 1.90 +0.05 −0.07 and σ � = 0.31 +0.05 −0.05 ) is consistent with previous works, and that the distribution = = = =

Electromagnetic transients powered by nuclear decay in the tidal tails of coalescing compact binaries

Abstract: We investigate the possibility that long tidal tails formed during compact object mergers may produce optical transients powered by the decay of freshly synthesized r-process material. Precise modeling of the merger dynamics allows for a realistic determination of the thermodynamic conditions in the ejected debris. We combine hydrodynamic and full nuclear network calculations to determine the resultant r-process abundances and the heating of the material by their decays. The subsequent homologous structure is mapped into a radiative transfer code to synthesize emergent model light curves and determine how their properties (variability and color evolution) depend on the mass ratio and orientation of the merging binary. The radiation emanating from the ejected debris, though less spectacular than a typical supernova, should be observable in transient surveys and we estimate the associated detection rates. We find that it is unlikely that photometry alone will be able to distinguish between different binary mass ratios and the nature of the compact objects, emphasizing the need for spectroscopic follow-up of these events. The case for (or against) compact object mergers as the progenitors of short gamma-ray bursts can be tested if such electromagnetic transients are detected (or not) in coincidence with some bursts, although they may be obscured by on-axis afterglows.

The extratropical upper troposphere and lower stratosphere

Abstract: [1] The extratropical upper troposphere and lower stratosphere (Ex-UTLS) is a transition region between the stratosphere and the troposphere. The Ex-UTLS includes the tropopause, a strong static stability gradient and dynamic barrier to transport. The barrier is reflected in tracer profiles. This region exhibits complex dynamical, radiative, and chemical characteristics that place stringent spatial and temporal requirements on observing and modeling systems. The Ex-UTLS couples the stratosphere to the troposphere through chemical constituent transport (of, e.g., ozone), by dynamically linking the stratospheric circulation with tropospheric wave patterns, and via radiative processes tied to optically thick clouds and clear-sky gradients of radiatively active gases. A comprehensive picture of the Ex-UTLS is presented that brings together different definitions of the tropopause, focusing on observed dynamical and chemical structure and their coupling. This integral view recognizes that thermal gradients and dynamic barriers are necessarily linked, that these barriers inhibit mixing and give rise to specific trace gas distributions, and that there are radiative feedbacks that help maintain this structure. The impacts of 21st century anthropogenic changes to the atmosphere due to ozone recovery and climate change will be felt in the Ex-UTLS, and recent simulations of these effects are summarized and placed in context.

The TW Hya Disk at 870 microns: Comparison of CO and Dust Radial Structures

Abstract: We present high resolution (0.3" = 16 AU), high signal-to-noise ratio Submillimeter Array observations of the 870 microns (345 GHz) continuum and CO J=3--2 line emission from the protoplanetary disk around TW Hya. Using continuum and line radiative transfer calculations, those data and the multiwavelength spectral energy distribution are analyzed together in the context of simple two-dimensional parametric disk structure models. Under the assumptions of a radially invariant dust population and (vertically integrated) gas-to-dust mass ratio, we are unable to simultaneously reproduce the CO and dust observations with model structures that employ either a single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared scattered light) that extends to a much larger radius of at least 215 AU. We discuss some possible explanations for the observed radial distribution of millimeter-sized dust grains and the apparent CO-dust size discrepancy, and suggest that they may be hallmarks of substructure in the dust disk or natural signatures of the growth and radial drift of solids that might be expected for disks around older pre-main sequence stars like TW Hya.

Three-dimensional radiative transfer modeling of AGN dusty tori as a clumpy two-phase medium

Abstract: We investigate the emission of active galactic nuclei (AGN) dusty tori in the infrared domain. Following theoretical predictions coming from hydrodynamical simulations, we model the dusty torus as a 3D two-phase medium with high-density clumps and low-density medium filling the space between the clumps. Spectral energy distributions (SED) and images of the torus at different wavelengths are obtained using 3D Monte Carlo radiative transfer code SKIRT. Our approach of generating clumpy structure allows us to model tori with single clumps, complex structures of merged clumps or interconnected sponge-like structure. A corresponding set of clumps-only models and models with smooth dust distribution is calculated for comparison. We found that dust distribution, optical depth, clump size and their actual arrangement in the innermost region, all have an impact on the shape of near- and mid-infrared SED. The 10 µm silicate feature can be suppressed for some parameters, but models with smooth dust distribution are also able to produce a wide range of the sili cate feature strength. Finally, we find that having the dust distributed in a two-phase medium, m ight offer a natural solution to the lack of emission in the near-infrared, compared to observed data, which affects clumpy models currently available in the literature.

Efficient three-dimensional nlte dust radiative transfer with skirt

Abstract: We present an updated version of SKIRT, a three-dimensional (3D) Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate, and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, two-dimensional, local thermal equilibrium (LTE) model to a fully 3D model that includes non-LTE (NLTE) dust emission and a clumpy structure of the dusty interstellar medium. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures.

Efficient 3D NLTE dust radiative transfer with SKIRT

Abstract: We present an updated version of SKIRT, a 3D Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, 2D, LTE model to a fully 3D model that includes NLTE dust emission and a clumpy structure of the dusty ISM. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures.

A Comparison of Methods for Determining the Molecular Content of Model Galaxies

Abstract: Recent observations indicate that star formation occurs only in the molecular phase of a galaxy's interstellar medium. A realistic treatment of star formation in simulations and analytic models of galaxies therefore requires that one determine where the transition from the atomic to molecular gas occurs. In this paper, we compare two methods for making this determination in cosmological simulations where the internal structures of molecular clouds are unresolved: a complex time-dependent chemistry network coupled to a radiative transfer calculation of the dissociating ultraviolet (UV) radiation field and a simple time-independent analytic approximation. We show that these two methods produce excellent agreement at all metallicities 10–2 of the Milky Way value across a very wide range of UV fields. At lower metallicities the agreement is worse, likely because time-dependent effects become important; however, there are no observational calibrations of molecular gas content at such low metallicities, so it is unclear if either method is accurate. The comparison suggests that, in many but not all applications, the analytic approximation provides a viable and nearly cost-free alternative to full time-dependent chemistry and radiative transfer.

Multiangle implementation of atmospheric correction (MAIAC): 1. Radiative transfer basis and look-up tables

Abstract: [1] This paper describes a radiative transfer basis of the algorithm MAIAC which performs simultaneous retrievals of atmospheric aerosol and bidirectional surface reflectance from the Moderate Resolution Imaging Spectroradiometer (MODIS). The retrievals are based on an accurate semianalytical solution for the top-of-atmosphere reflectance expressed as an explicit function of three parameters of the Ross–Thick Li–Sparse model of surface bidirectional reflectance. This solution depends on certain functions of atmospheric properties and geometry which are precomputed in the look-up table (LUT). This paper further considers correction of the LUT functions for variations of surface pressure/height and of atmospheric water vapor, which is a common task in the operational remote sensing. It introduces a new analytical method for the water vapor correction of the multiple-scattering path radiance. It also summarizes the few basic principles that provide a high efficiency and accuracy of the LUT-based radiative transfer for the aerosol/surface retrievals and optimize the size of LUT. For example, the single-scattering path radiance is calculated analytically for a given surface pressure and atmospheric water vapor. The same is true for the direct surface-reflected radiance, which along with the single-scattering path radiance largely defines the angular dependence of measurements. For these calculations, the aerosol phase functions and kernels of the surface bidirectional reflectance model are precalculated at a high angular resolution. The other radiative transfer functions depend rather smoothly on angles because of multiple scattering and can be calculated at coarser angular resolution to reduce the LUT size. At the same time, this resolution should be high enough to use the nearest neighbor geometry angles to avoid costly three-dimensional interpolation. The pressure correction is implemented via linear interpolation between two LUTs computed for the standard and reduced pressure levels. A linear mixture and a modified linear mixture methods are used to represent different aerosol types in the aerosol/surface retrievals from several base models of the fine and coarse aerosol fractions. In summary, the developed LUT algorithm allows fast high-accuracy simulations of the outgoing radiance with full variability of the atmospheric and surface bidirectional reflectance properties for the aerosol/surface remote sensing.

Dust size distributions in coagulation/fragmentation equilibrium: numerical solutions and analytical fits

Abstract: Context. Grains in circumstellar disks are believed to grow by mutual collisions and subsequent sticking due to surface forces. Results of many fields of research involving circumstellar disks, such as radiative transfer calculations, disk chemistry, magneto-hydrodynamic simulations largely depend on the unknown grain size distribution.Aims. As detailed calculations of grain growth and fragmentation are both numerically challenging and computationally expensive, we aim to find simple recipes and analytical solutions for the grain size distribution in circumstellar disks for a scenario in which grain growth is limited by fragmentation and radial drift can be neglected.Methods. We generalize previous analytical work on self-similar steady-state grain distributions. Numerical simulations are carried out to identify under which conditions the grain size distributions can be understood in terms of a combination of power-law distributions. A physically motivated fitting formula for grain size distributions is derived using our analytical predictions and numerical simulations.Results. We find good agreement between analytical results and numerical solutions of the Smoluchowski equation for simple shapes of the kernel function. The results for more complicated and realistic cases can be fitted with a physically motivated “black box” recipe presented in this paper. Our results show that the shape of the dust distribution is mostly dominated by the gas surface density (not the dust-to-gas ratio), the turbulence strength and the temperature and does not obey an MRN type distribution.

Environmental Dependence of the Kennicutt?Schmidt Relation in Galaxies

Abstract: We present a detailed description of a phenomenological H2 formation model and local star formation prescription based on the density of molecular (rather than total) gas. Such an approach allows us to avoid the arbitrary density and temperature thresholds typically used in star formation recipes in galaxy formation simulations. We present results of the model based on realistic cosmological simulations of high-z galaxy formation for a grid of numerical models with varied dust-to-gas ratios and interstellar far-UV (FUV) fluxes. Our results show that both the atomic-to-molecular transition on small, tens-of-parsec scales and the Kennicutt-Schmidt (K-S) relation on large, kiloparsec scales are sensitive to the dust-to-gas ratio and the FUV flux. The atomic-to-molecular transition as a function of gas density or column density has a large scatter but is rather sharp and shifts to higher densities with decreasing dust-to-gas ratio and/or increasing FUV flux. Consequently, star formation is concentrated to higher gas surface density regions, resulting in steeper slope and lower amplitude of the K-S relation at a given ΣH, in less dusty and/or higher FUV flux environments. We parameterize the dependences observed in our simulations in convenient fitting formulae, which can be used to model the dependence of the K-S relation on the dust-to-gas ratio and FUV flux in semi-analytic models and in cosmological simulations that do not include radiative transfer and H2 formation. Finally, we show that ionized gas can contribute a significant fraction of the total gas surface density in environments typical for high-redshift galaxies.

ARTS, the atmospheric radiative transfer simulator, version 2

Abstract: The second version of the atmospheric radiative transfer simulator, ARTS, is introduced. This is a general software package for long wavelength radiative transfer simulations, with a focus on passive microwave observations. The core part provides a workspace environment, in line with script languages. New for this version is an agenda mechanism that gives a high degree of modularity. The framework is intended to be as general as possible: the polarisation state can be fully described, the model atmosphere can be one- (1D), two- (2D) or three-dimensional (3D), a full description of geoid and surface is possible, observation geometries from the ground, from satellite, and from aeroplane or balloon are handled, and surface reflection can be treated in simple or complex manners. Remote sensing applications are supported by a comprehensive and efficient treatment of sensor characteristics. Jacobians can be calculated for the most important atmospheric variables in non-scattering conditions. Finally, the most prominent feature is the rigorous treatment of scattering that has been implemented in two modules: a discrete ordinate iterative approach mainly used for 1D atmospheres, and a Monte Carlo approach which is the preferred algorithm for 3D atmospheres. ARTS is freely available, and maintained as an open-source project.

Collapse of massive magnetized dense cores using radiation magnetohydrodynamics: early fragmentation inhibition

Abstract: We report the results of radiation-magnetohydrodynamics calculations in the context of high-mass star formation, using for the first time a self-consistent model for photon emission (i.e., via thermal emission and in radiative shocks) and with the high resolution necessary to properly resolve magnetic braking effects and radiative shocks on scales <100 AU. We investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores. We identify a new mechanism that inhibits initial fragmentation of massive dense cores where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. Magnetic braking is transporting angular momentum outward and is lowering the rotational support and is thus increasing the infall velocity. This enhances the radiative feedback owing to the accretion shock on the first core. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation while moderately magnetized massive dense cores are more appropriate forming OB associations or small star clusters.

HYPERION: an open-source parallelized three-dimensional dust continuum radiative transfer code

Abstract: Hyperion is a new three-dimensional dust continuum Monte-Carlo radiative transfer code that is designed to be as generic as possible, allowing radiative transfer to be computed through a variety of three-dimensional grids. The main part of the code is problemindependent, and only requires an arbitrary three-dimensional density structure, dust properties, the position and properties of the illuminating sources, and parameters controlling the running and output of the code. Hyperion is parallelized, and is shown to scale well to thousands of processes. Two common benchmark models for protoplanetary disks were computed, and the results are found to be in excellent agreement with those from other codes. Finally, to demonstrate the capabilities of the code, dust temperatures, SEDs, and synthetic multi-wavelength images were computed for a dynamical simulation of a low-mass star formation region. Hyperion is being actively developed to include new features, and is publicly available �� .

Photodissociation of H2 in protogalaxies: modelling self‐shielding in three‐dimensional simulations

Abstract: The ability of primordial gas to cool in protogalactic haloes exposed to ultraviolet (UV) radiation is critically dependent on the self-shielding of H2. We perform radiative transfer calculations of LW line photons, post-processing outputs from three-dimensional adaptive mesh refinement simulations of haloes with Tvir≳ 104 K at z ∼ 10. We calculate the optically thick photodissociation rate numerically, including the effects of density, temperature and velocity gradients in the gas, as well as line overlap and shielding of H2 by H i, over a large number of sightlines. In low-density regions (n ≲ 104 cm−3) the dissociation rates exceed those obtained using most previous approximations by more than an order of magnitude; the correction is smaller at higher densities. We trace the origin of the deviations primarily to inaccuracies of (i) the most common fitting formula for the suppression of the dissociation rate, from Draine and Bertoldi and (ii) estimates for the effective shielding column density from local properties of the gas. The combined effects of gas temperature and velocity gradients are comparatively less important, typically altering the spherically averaged rate only by a factor of ≲2. We present a simple modification to the Draine & Bertoldi fitting formula for the optically thick rate which improves agreement with our numerical results to within ∼15 per cent, and can be adopted in future simulations. We find that estimates for the effective shielding column can be improved by using the local Sobolev length. Our correction to the H2 self-shielding reduces the critical LW flux to suppress H2 cooling in Tvir≳ 104 K haloes by an order of magnitude; this increases the number of such haloes in which supermassive (M ∼ 105 M⊙) black holes may have formed.

RAMI4PILPS: An intercomparison of formulations for the partitioning of solar radiation in land surface models

Abstract: [1] Remotely sensed, multiannual data sets of shortwave radiative surface fluxes are now available for assimilation into land surface schemes (LSSs) of climate and/or numerical weather prediction models. The RAMI4PILPS suite of virtual experiments assesses the accuracy and consistency of the radiative transfer formulations that provide the magnitudes of absorbed, reflected, and transmitted shortwave radiative fluxes in LSSs. RAMI4PILPS evaluates models under perfectly controlled experimental conditions in order to eliminate uncertainties arising from an incomplete or erroneous knowledge of the structural, spectral and illumination related canopy characteristics typical for model comparison with in situ observations. More specifically, the shortwave radiation is separated into a visible and near-infrared spectral region, and the quality of the simulated radiative fluxes is evaluated by direct comparison with a 3-D Monte Carlo reference model identified during the third phase of the Radiation transfer Model Intercomparison (RAMI) exercise. The RAMI4PILPS setup thus allows to focus in particular on the numerical accuracy of shortwave radiative transfer formulations and to pinpoint to areas where future model improvements should concentrate. The impact of increasing degrees of structural and spectral subgrid variability on the simulated fluxes is documented and the relevance of any thus emerging biases with respect to gross primary production estimates and shortwave radiative forcings due to snow and fire events are investigated.

Photodissociation of H2 in Protogalaxies: Modeling Self-Shielding in 3D Simulations

Abstract: The ability of primordial gas to cool in proto-galactic haloes exposed to Lyman-Werner (LW) radiation is critically dependent on the self-shielding of H_2. We perform radiative transfer calculations of LW line photons, post-processing outputs from three-dimensional adaptive mesh refinement (AMR) simulations of haloes with T_vir > 10^4 K at redshifts around z=10. We calculate the optically thick photodissociation rate numerically, including the effects of density, temperature, and velocity gradients in the gas, as well as line overlap and shielding of H_2 by HI, over a large number of sight-lines. In low-density regions (n 10^4 K haloes by an order of magnitude; this increases the number of such haloes in which supermassive (approx. M=10^5 M_sun) black holes may have formed.

Empirical constraints on turbulence in protoplanetary accretion disks

Abstract: We present arcsecond-scale Submillimeter Array observations of the CO(3-2) line emission from the disks around the young stars HD 163296 and TW Hya at a spectral resolution of 44 m s–1. These observations probe below the ~100 m s–1 turbulent linewidth inferred from lower-resolution observations, and allow us to place constraints on the turbulent linewidth in the disk atmospheres. We reproduce the observed CO(3-2) emission using two physical models of disk structure: (1) a power-law temperature distribution with a tapered density distribution following a simple functional form for an evolving accretion disk, and (2) the radiative transfer models developed by D'Alessio et al. that can reproduce the dust emission probed by the spectral energy distribution. Both types of models yield a low upper limit on the turbulent linewidth (Doppler b-parameter) in the TW Hya system (40 m s–1) and a tentative (3σ) detection of a ~300 m s–1 turbulent linewidth in the upper layers of the HD 163296 disk. These correspond to roughly ≤10% and 40% of the sound speed at size scales commensurate with the resolution of the data. The derived linewidths imply a turbulent viscosity coefficient, α, of order 0.01 and provide observational support for theoretical predictions of subsonic turbulence in protoplanetary accretion disks.

High‐resolution solar spectral irradiance from extreme ultraviolet to far infrared

Abstract: [1] This paper presents new extremely high-resolution solar spectral irradiance (SSI) calculations covering wavelengths from 0.12 nm to 100 micron obtained by the Solar Irradiance Physical Modeling (SRPM) system. Daily solar irradiance spectra were constructed for most of Solar Cycle 23 based on a set of physical models of the solar features and non-LTE calculations of their emitted spectra as function of viewing angle, and solar images specifying the distribution of features on the solar disk. Various observational tests are used to assess the quality of the spectra provided here. The present work emphasizes the effects on the SSI of the upper chromosphere and full-non-LTE radiative transfer calculation of level populations and ionizations that are essential for physically consistent results at UV wavelengths and for deep lines in the visible and IR. This paper also considers the photodissociation continuum opacity of molecular species, e.g., CH and OH, and proposes the consideration of NH photodissociation which can solve the puzzle of the missing near-UV opacity in the spectral range of the near-UV. Finally, this paper is based on physical models of the solar atmosphere and extends the previous lower-layer models into the upper-transition-region and coronal layers that are the dominant source of photons at wavelengths shorter than ∼50 nm (except for the He II 30.4 nm line, mainly formed in the lower-transition-region).

enzo+moray: radiation hydrodynamics adaptive mesh refinement simulations with adaptive ray tracing

Abstract: We describe a photon-conserving radiative transfer algorithm, using a spatially-adaptive ray-tracing scheme, and its parallel implementation into the adaptive mesh refinement cosmological hydrodynamics code ENZO. By coupling the solver with the energy equation and non-equilibrium chemistry network, our radiation hydrodynamics framework can be utilized to study a broad range of astrophysical problems, such as stellar and black hole feedback. Inaccuracies can arise from large time-steps and poor sampling; therefore, we devised an adaptive time-stepping scheme and a fast approximation of the optically-thin radiation field with multiple sources. We test the method with several radiative transfer and radiation hydrodynamics tests that are given in Iliev et al. We further test our method with more dynamical situations, for example, the propagation of an ionization front through a Rayleigh―Taylor instability, time-varying luminosities and collimated radiation. The test suite also includes an expanding H II region in a magnetized medium, utilizing the newly implemented magnetohydrodynamics module in ENZO. This method linearly scales with the number of point sources and number of grid cells. Our implementation is scalable to 512 processors on distributed memory machines and can include the radiation pressure and secondary ionizations from X-ray radiation. It is included in the newest public release of ENZO.

Revisiting the radiative impact of dust on Mars using the LMD Global Climate Model

Abstract: [1] Airborne dust is the main driver of Martian atmospheric temperature, and accurately accounting for its radiative effect in Global Climate Models (GCMs) is essential. This requires the modeling of the dust distribution and radiative properties, and when trying to simulate the true climate variability, the use of the observed dust column opacity to guide the model. A recurrent problem has been the inability of Mars GCMs to predict realistic temperatures while using both the observed dust radiative properties and column opacity. One would have to drive the model with a tuned opacity to reach an agreement with the observations, thereby losing its self-consistency. In this paper, we show that using the most recently derived dust radiative properties in the LMD (Laboratoire de Meteorologie Dynamique) GCM solves this problem, which was mainly due to the underestimation of the dust single scattering albedo in the solar domain. However, an overall warm temperature bias remains above the 1 hPa pressure level. We therefore refine the model by implementing a “semi-interactive” dust transport scheme which is coupled to the radiative transfer calculations. This scheme allows a better representation of the dust layer depth in the model and thereby removes the remaining warm bias. The LMD/GCM is now able to predict accurate temperatures without any tuning of the dust opacity used to guide the model. Remaining discrepancies are discussed, and seem to be primarily due to the neglect of the radiative effect of water-ice clouds, and secondarily to persisting uncertainties in the dust spatial distribution.

Modelling the spectral energy distribution of galaxies - V. The dust and PAH emission SEDs of disk galaxies

Abstract: We present a self-consistent model of the spectral energy distributions (SEDs) of spiral galaxies from the ultraviolet (UV) to the mid-infrared (MIR)/far-infrared (FIR)/submillimeter (submm) based on a full radiative transfer calculation of the propagation of starlight in galaxy disks. This model predicts not only the total integrated energy absorbed in the UV/optical and re-emitted in the infrared/submm, but also the colours of the dust emission based on an explicit calculation of the strength and colour of the UV/optical radiation fields heating the dust, and incorporating a full calculation of the stochastic heating of small dust grains and PAH molecules. The geometry of the translucent components of the model is empirically constrained using the results from the radiation transfer analysis of Xilouris et al. on spirals in the middle range of the Hubble sequence, while the geometry of the optically thick components is constrained from physical considerations with a posteriori checks of the model predictions with observational data. Following the observational constraints, the model has both a distribution of diffuse dust associated with the old and young disk stellar populations as well as a clumpy component arising from dust in the parent molecular clouds in star forming regions. In accordance with the fragmented nature of dense molecular gas in typical star-forming regions, UV light from massive stars is allowed to either freely stream away into the diffuse medium in some fraction of directions or be geometrically blocked and locally absorbed in clumps. These geometrical constraints enable the dust emission to be predicted in terms of a minimum set of free parameters: the central face-on dust opacity in the B -band , a clumpiness factor F for the star-forming regions, the star-formation rate SFR , the normalised luminosity of the old stellar population old and the bulge-to-disk ratio B /D . We show that these parameters are almost orthogonal in their predicted effect on the colours of the dust/PAH emission. In most practical applications B /D will actually not be a free parameter but (together with the angular size θ gal and inclination i of the disk) act as a constraint derived from morphological decomposition of higher resolution optical images. This also extends the range of applicability of the model along the Hubble sequence. We further show that the dependence of the dust emission SED on the colour of the stellar photon field depends primarily on the ratio between the luminosities of the young and old stellar populations (as specified by the parameters SFR and old ) rather than on the detailed colour of the emissions from either of these populations. The model is thereby independent of a priori assumptions of the detailed mathematical form of the dependence of SFR on time, allowing UV/optical SEDs to be dereddened without recourse to population synthesis models. Utilising these findings, we show how the predictive power of radiative transfer calculations can be combined with measurements of θ gal , i and B /D from optical images to self-consistently fit UV/optical-MIR/FIR/submm SEDs observed in large statistical surveys in a fast and flexible way, deriving physical parameters on an object-by-object basis. We also identify a non-parametric test of the fidelity of the model in practical applications through comparison of the model predictions for FIR colour and surface brightness with the corresponding observed quantities. This should be effective in identifying objects such as AGNs or star-forming galaxies with markedly different geometries to those of the calibrators of Xilouris et al. The results of the calculations are made available in the form of a large library of simulated dust emission SEDs spanning the whole parameter space of our model, together with the corresponding library of dust attenuation calculated using the same model.

Comparison of reionization models: radiative transfer simulations and approximate, seminumeric models

Abstract: We compare the predictions of four different algorithms for the distribution of ionized gas during the Epoch of Reionization. These algorithms are all used to run a 100 Mpc h−1 simulation of reionization with the same initial conditions. Two of the algorithms are state-of-the-art ray-tracing radiative transfer codes that use disparate methods to calculate the ionization history. The other two algorithms are fast but more approximate schemes based on iterative applications of a smoothing filter on the underlying source and density fields. We compare these algorithms' resulting ionization and 21-cm fields using several different statistical measures. The two radiative transfer schemes are in excellent agreement judging by the power spectra of both the ionization fields and the 21-cm emission fields (agreeing to better than 10 per cent) and are in good agreement with the analytic schemes (better than 50 per cent) over the range of ionized fractions and wavevectors we compare ( Mpc−1). This agreement suggests that the different approximations involved in the ray-tracing algorithms are sensible and that seminumerical schemes provide a numerically inexpensive, yet fairly accurate, description of the reionization process.

MiRS: An All-Weather 1DVAR Satellite Data Assimilation and Retrieval System

Abstract: A 1-D variational system has been developed to process spaceborne measurements. It is an iterative physical inversion system that finds a consistent geophysical solution to fit all radiometric measurements simultaneously. One of the particularities of the system is its applicability in cloudy and precipitating conditions. Although valid, in principle, for all sensors for which the radiative transfer model applies, it has only been tested for passive microwave sensors to date. The Microwave Integrated Retrieval System (MiRS) inverts the radiative transfer equation by finding radiometrically appropriate profiles of temperature, moisture, liquid cloud, and hydrometeors, as well as the surface emissivity spectrum and skin temperature. The inclusion of the emissivity spectrum in the state vector makes the system applicable globally, with the only differences between land, ocean, sea ice, and snow backgrounds residing in the covariance matrix chosen to spectrally constrain the emissivity. Similarly, the inclusion of the cloud and hydrometeor parameters within the inverted state vector makes the assimilation/inversion of cloudy and rainy radiances possible, and therefore, it provides an all-weather capability to the system. Furthermore, MiRS is highly flexible, and it could be used as a retrieval tool (independent of numerical weather prediction) or as an assimilation system when combined with a forecast field used as a first guess and/or background. In the MiRS, the fundamental products are inverted first and then are interpreted into secondary or derived products such as sea ice concentration, snow water equivalent (based on the retrieved emissivity) rainfall rate, total precipitable water, integrated cloud liquid amount, and ice water path (based on the retrieved atmospheric and hydrometeor products). The MiRS system was implemented operationally at the U.S. National Oceanic and Atmospheric Administration (NOAA) in 2007 for the NOAA-18 satellite. Since then, it has been extended to run for NOAA-19, Metop-A, and DMSP-F16 and F18 SSMI/S. This paper gives an overview of the system and presents brief results of the assessment effort for all fundamental and derived products.