Showing papers on "Radiative transfer published in 2014"

Improved Reflection Models of Black Hole Accretion Disks: Treating the Angular Distribution of X-Rays

Abstract: X-ray reflection models are used to constrain the properties of the accretion disk, such as the degree of ionization of the gas and the elemental abundances. In combination with general relativistic ray tracing codes, additional parameters like the spin of the black hole and the inclination to the system can be determined. However, current reflection models used for such studies only provide angle-averaged solutions for the flux reflected at the surface of the disk. Moreover, the emission angle of the photons changes over the disk due to relativistic light bending. To overcome this simplification, we have constructed an angle-dependent reflection model with the XILLVER code and self-consistently connected it with the relativistic blurring code RELLINE. The new model, relxill, calculates the proper emission angle of the radiation at each point on the accretion disk and then takes the corresponding reflection spectrum into account. We show that the reflected spectra from illuminated disks follow a limb-brightening law highly dependent on the ionization of disk and yet different from the commonly assumed form Iln (1 + 1/μ). A detailed comparison with the angle-averaged model is carried out in order to determine the bias in the parameters obtained by fitting a typical relativistic reflection spectrum. These simulations reveal that although the spin and inclination are mildly affected, the Fe abundance can be overestimated by up to a factor of two when derived from angle-averaged models. The fit of the new model to the Suzaku observation of the Seyfert galaxy Ark 120 clearly shows a significant improvement in the constraint of the physical parameters, in particular by enhancing the accuracy in the inclination angle and the spin determinations.

Land Surface Temperature Retrieval from Landsat 8 TIRS—Comparison between Radiative Transfer Equation-Based Method, Split Window Algorithm and Single Channel Method

Abstract: Accurate inversion of land surface geo/biophysical variables from remote sensing data for earth observation applications is an essential and challenging topic for the global change research. Land surface temperature (LST) is one of the key parameters in the physics of earth surface processes from local to global scales. The importance of LST is being increasingly recognized and there is a strong interest in developing methodologies to measure LST from the space. Landsat 8 Thermal Infrared Sensor (TIRS) is the newest thermal infrared sensor for the Landsat project, providing two adjacent thermal bands, which has a great benefit for the LST inversion. In this paper, we compared three different approaches for LST inversion from TIRS, including the radiative transfer equation-based method, the split-window algorithm and the single channel method. Four selected energy balance monitoring sites from the Surface Radiation Budget Network (SURFRAD) were used for validation, combining with the MODIS 8 day emissivity product. For the investigated sites and scenes, results show that the LST inverted from the radiative transfer equation-based method using band 10 has the highest accuracy with RMSE lower than 1 K, while the SW algorithm has moderate accuracy and the SC method has the lowest accuracy.

Planck intermediate results. XIX. An overview of the polarized thermal emission from Galactic dust

Abstract: This paper presents the large-scale polarized sky as seen by Planck HFI at 353 GHz, which is the most sensitive Planck channel for dust polarization. We construct and analyse large-scale maps of dust polarization fraction and polarization direction, while taking account of noise bias and possible systematic effects. We find that the maximum observed dust polarization fraction is high (pmax > 18%), in particular in some of the intermediate dust column density (AV < 1mag) regions. There is a systematic decrease in the dust polarization fraction with increasing dust column density, and we interpret the features of this correlation in light of both radiative grain alignment predictions and fluctuations in the magnetic field orientation. We also characterize the spatial structure of the polarization angle using the angle dispersion function and find that, in nearby fields at intermediate latitudes, the polarization angle is ordered over extended areas that are separated by filamentary structures, which appear as interfaces where the magnetic field sky projection rotates abruptly without apparent variations in the dust column density. The polarization fraction is found to be anti-correlated with the dispersion of the polarization angle, implying that the variations are likely due to fluctuations in the 3D magnetic field orientation along the line of sight sampling the diffuse interstellar medium.We also compare the dust emission with the polarized synchrotron emission measured with the Planck LFI, with low-frequency radio data, and with Faraday rotation measurements of extragalactic sources. The two polarized components are globally similar in structure along the plane and notably in the Fan and North Polar Spur regions. A detailed comparison of these three tracers shows, however, that dust and cosmic rays generally sample different parts of the line of sight and confirms that much of the variation observed in the Planck data is due to the 3D structure of the magnetic field.

A Global Three Dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks

Abstract: We study super-Eddington accretion flows onto black holes using a global three-dimensional radiation magneto-hydrodynamical simulation. We solve the time-dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ∼220 L {sub Edd}/c {sup 2} and forms a radiation-driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ∼20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ∼10 L {sub Edd}. This yields a radiative efficiency ∼4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed publishedmore » results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.« less

Three-dimensional general relativistic radiation magnetohydrodynamical simulation of super-Eddington accretion, using a new code HARMRAD with M1 closure

Abstract: Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant application using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20\%$, the large-scale electromagnetic jet efficiency is of order $10\%$, and the total radiative efficiency that reaches large distances remains low at only order $1\%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs.\ disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $z\sim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.

Seismic diagnostics for transport of angular momentum in stars 1. Rotational splittings from the PMS to the RGB

Abstract: Context. Rotational splittings are currently measured for several main sequence stars and a large number of red giants with the space mission Kepler. This will provide stringent constraints on rotation profiles. Aims. Our aim is to obtain seismic constraints on the internal tran sport and surface loss of angular momentum of oscillating solar-like stars. To this end, we study the evolution of rotational spli ttings from the pre-main sequence to the red-giant branch for stochastically excited oscillation modes. Methods. We modified the evolutionary code CESAM2K to take rotational ly induced transport in radiative zones into account. Linear rotational splittings were computed for a sequence of 1.3M⊙ models. Rotation profiles were derived from our evolutionar y models and eigenfunctions from linear adiabatic oscillation calc ulations. Results. We find that transport by meridional circulation and shear tu rbulence yields far too high a core rotation rate for red-gia nt models compared with recent seismic observations. We discuss several uncertainties in the physical description of sta rs that could have an impact on the rotation profiles. For instance, we find t hat the Goldreich-Schubert-Fricke instability does not extract enough angular momentum from the core to account for the discrepancy. In contrast, an increase of the horizontal turbulent visc osity by 2 orders of magnitude is able to significantly decrease the cen tral rotation rate on the red-giant branch. Conclusions. Our results indicate that it is possible that the prescripti on for the horizontal turbulent viscosity largely underest imates its actual value or else a mechanism not included in current stellar models of low mass stars is needed to slow down the rotation in the radiative core of red-giant stars.

Numerical simulations of super-critical black hole accretion flows in general relativity

Abstract: A new general relativistic radiation magnetohydrodynamical code KORAL, is described, which employs the M1 scheme to close the radiation moment equations. The code has been successfully verified against a number of tests. Axisymmetric simulations of super-critical magnetized accretion on a non-rotating black hole (a=0.0) and a spinning black hole (a=0.9) are presented. The accretion rates in the two models are \dot M = 100-200 \dot M_Edd. These first general relativistic simulations of super-critical black hole accretion are potentially relevant to tidal disruption events and hyper-accreting supermassive black holes in the early universe. Both simulated models are optically and geometrically thick, and have funnels through which energy escapes in the form of relativistic gas, Poynting flux and radiative flux. The jet is significantly more powerful in the a=0.9 run. The net energy outflow rate in the two runs correspond to efficiencies of 5% (a=0) and 33% (a=0.9), as measured with respect to the mass accretion rate at the black hole. These efficiencies agree well with those measured in previous simulations of non-radiative geometrically thick disks. Furthermore, in the a=0.9 run, the outflow power appears to originate in the spinning black hole, suggesting that the associated physics is again similar in non-radiative and super-critical accretion flows. While the two simulations are efficient in terms of total energy outflow, both runs are radiatively inefficient. Their luminosities are only \sim 1-10 L_Edd, which corresponds to a radiative efficiency \sim 0.1%. Interestingly, most of the radiative luminosity emerges through the funnels, which subtend a very small solid angle. Therefore, measured in terms of a local radiative flux, the emitted radiation is highly super-Eddington.

Lyman Alpha Emitting Galaxies as a Probe of Reionization

Abstract: The Epoch of Reionization (EoR) represents a milestone in the evolution of our Universe. Star-forming galaxies that existed during the EoR likely emitted a significant fraction (~5-40%) of their bolometric luminosity as Lyman Alpha (Lya) line emission. However, neutral intergalactic gas that existed during the EoR was opaque to Lya emission that escaped from galaxies during this epoch, which makes it difficult to observe. The neutral intergalactic medium (IGM) may thus reveal itself by suppressing the Lya flux from background galaxies. Interestingly, a `sudden' reduction in the observed Lya flux has now been observed in galaxies at z >6. This review contains a detailed summary of Lya radiative processes: I describe (i) the main Lya emission processes, including collisional-excitation & recombination (and derive the famous factor `0.68'), and (ii) basic radiative transfer concepts, including e.g. partially coherent scattering, frequency diffusion, resonant versus wing scattering, optically thick versus 'extremely' optically thick (static/outflowing/collapsing) media, and multiphase media. Following this review, I derive expressions for the Gunn-Peterson optical depth of the IGM during (inhomogeneous) reionization and post-reionization. I then describe why current observations appear to require a very rapid evolution of volume-averaged neutral fraction of hydrogen in the context of realistic inhomogeneous reionization models, and discuss uncertainties in this interpretation. Finally, I describe how existing & futures surveys and instruments can help reduce these uncertainties, and allow us to fully exploit Lya emitting galaxies as a probe of the EoR.

Radiative transfer through terrestrial atmosphere and ocean: Software package SCIATRAN

Abstract: SCIATRAN is a comprehensive software package for the modeling of radiative transfer processes in the terrestrial atmosphere and ocean in the spectral range from the ultraviolet to the thermal infrared ( 0.18 – 40 μ m ) including multiple scattering processes, polarization, thermal emission and ocean–atmosphere coupling. The software is capable of modeling spectral and angular distributions of the intensity or the Stokes vector of the transmitted, scattered, reflected, and emitted radiation assuming either a plane-parallel or a spherical atmosphere. Simulations are done either in the scalar or in the vector mode (i.e. accounting for the polarization) for observations by space-, air-, ship- and balloon-borne, ground-based, and underwater instruments in various viewing geometries (nadir, off-nadir, limb, occultation, zenith-sky, off-axis). All significant radiative transfer processes are accounted for. These are, e.g. the Rayleigh scattering, scattering by aerosol and cloud particles, absorption by gaseous components, and bidirectional reflection by an underlying surface including Fresnel reflection from a flat or roughened ocean surface. The software package contains several radiative transfer solvers including finite difference and discrete-ordinate techniques, an extensive database, and a specific module for solving inverse problems. In contrast to many other radiative transfer codes, SCIATRAN incorporates an efficient approach to calculate the so-called Jacobians, i.e. derivatives of the intensity with respect to various atmospheric and surface parameters. In this paper we discuss numerical methods used in SCIATRAN to solve the scalar and vector radiative transfer equation, describe databases of atmospheric, oceanic, and surface parameters incorporated in SCIATRAN, and demonstrate how to solve some selected radiative transfer problems using the SCIATRAN package. During the last decades, a lot of studies have been published demonstrating that SCIATRAN is a valuable tool for a wide range of remote sensing applications. Here, we present some selected comparisons of SCIATRAN simulations to published benchmark results, independent radiative transfer models, and various measurements from satellite, ground-based, and ship instruments. Methods for solving inverse problems related to remote sensing of the Earth's atmosphere using the SCIATRAN software are outside the scope of this study and will be discussed in a follow-up paper. The SCIATRAN software package along with a detailed User's Guide is freely available for non-commercial use via the webpage of the Institute of Environmental Physics (IUP), University of Bremen: http://www.iup.physik.uni-bremen.de/sciatran .

Lyα emitting galaxies as a probe of reionisation

Abstract: The Epoch of Reionization (EoR) represents a milestone in the evolution of our Universe. Star-forming galaxies that existed during the EoR likely emitted a significant fraction ( ~ 5 − 40%) of their bolometric luminosity as Lyα line emission. However, neutral intergalactic gas that existed during the EoR was opaque to Lyα emission that escaped from galaxies during this epoch, which makes it difficult to observe. The neutral intergalactic medium (IGM) may thus reveal itself by suppressing the Lyα flux from background galaxies. Interestingly, a ‘sudden’ reduction in the observed Lyα flux has now been observed in galaxies at z > 6. This review contains a detailed summary of Lyα radiative processes: I describe (i) the main Lyα emission processes, including collisional-excitation & recombination (and derive the origin of the famous factor ‘0.68’), and (ii) basic radiative transfer concepts, including e.g. partially coherent scattering, frequency diffusion, resonant versus wing scattering, optically thick versus ‘extremely’ optically thick (static/outflowing/collapsing) media, and multiphase media. Following this review, I derive expressions for the Gunn-Peterson optical depth of the IGM during (inhomogeneous) reionisation and post-reionisation. I then describe why current observations appear to require a very rapid evolution of volume-averaged neutral fraction of hydrogen in the context of realistic inhomogeneous reionisation models, and discuss uncertainties in this interpretation. Finally, I describe how existing & futures surveys and instruments can help reduce these uncertainties, and allow us to fully exploit Lyα emitting galaxies as a probe of the EoR.

Anisotropic failure of Fourier theory in time-domain thermoreflectance experiments

Abstract: The applicability of Fourier's law to heat transfer problems relies on the assumption that heat carriers have mean free paths smaller than important length scales of the temperature profile. This assumption is not generally valid in nanoscale thermal transport problems where spacing between boundaries is small (<1 μm), and temperature gradients vary rapidly in space. Here we study the limits to Fourier theory for analysing three-dimensional heat transfer problems in systems with an interface. We characterize the relationship between the failure of Fourier theory, phonon mean free paths, important length scales of the temperature profile and interfacial-phonon scattering by time-domain thermoreflectance experiments on Si, Si0.99Ge0.01, boron-doped Si and MgO crystals. The failure of Fourier theory causes anisotropic thermal transport. In situations where Fourier theory fails, a simple radiative boundary condition on the heat diffusion equation cannot adequately describe interfacial thermal transport.

Control of radiative processes using tunable plasmonic nanopatch antennas.

Abstract: The radiative processes associated with fluorophores and other radiating systems can be profoundly modified by their interaction with nanoplasmonic structures. Extreme electromagnetic environments can be created in plasmonic nanostructures or nanocavities, such as within the nanoscale gap region between two plasmonic nanoparticles, where the illuminating optical fields and the density of radiating modes are dramatically enhanced relative to vacuum. Unraveling the various mechanisms present in such coupled systems, and their impact on spontaneous emission and other radiative phenomena, however, requires a suitably reliable and precise means of tuning the plasmon resonance of the nanostructure while simultaneously preserving the electromagnetic characteristics of the enhancement region. Here, we achieve this control using a plasmonic platform consisting of colloidally synthesized nanocubes electromagnetically coupled to a metallic film. Each nanocube resembles a nanoscale patch antenna (or nanopatch) whose ...

Global climate modeling of the Martian water cycle with improved microphysics and radiatively active water ice clouds

Abstract: Water ice clouds play a key role in the radiative transfer of the Martian atmosphere, impacting its thermal structure, its circulation, and, in turn, the water cycle. Recent studies including the radiative effects of clouds in global climate models (GCMs) have found that the corresponding feedbacks amplify the model defaults. In particular, it prevents models with simple microphysics from reproducing even the basic characteristics of the water cycle. Within that context, we propose a new implementation of the water cycle in GCMs, including a detailed cloud microphysics taking into account nucleation on dust particles, ice particle growth, and scavenging of dust particles due to the condensation of ice. We implement these new methods in the Laboratoire de Meteorologie Dynamique GCM and find satisfying agreement with the Thermal Emission Spectrometer observations of both water vapor and cloud opacities, with a significant improvement when compared to GCMs taking into account radiative effects of water ice clouds without this implementation. However, a lack of water vapor in the tropics after Ls = 180° is persistent in simulations compared to observations, as a consequence of aphelion cloud radiative effects strengthening the Hadley cell. Our improvements also allow us to explore questions raised by recent observations of the Martian atmosphere. Supersaturation above the hygropause is predicted in line with Spectroscopy for Investigation of Characteristics of the Atmosphere of Mars observations. The model also suggests for the first time that the scavenging of dust by water ice clouds alone fails to fully account for the detached dust layers observed by the Mars Climate Sounder.

HI-to-H2 Transitions and H I Column Densities in Galaxy Star-Forming Regions

Abstract: We present new analytic theory and radiative transfer computations for the atomic to molecular (HI-to-H2) transitions, and the build-up of atomic-hydrogen (HI) gas columns, in optically thick interstellar clouds, irradiated by far-ultraviolet photodissociating radiation fields. We derive analytic expressions for the total HI column densities for (1D) planar slabs, for beamed or isotropic radiation fields, from the weak- to strong-field limits, for gradual or sharp atomic to molecular transitions, and for arbitrary metallicity. Our expressions may be used to evaluate the HI column densities as functions of the radiation field intensity and the H2-dust-limited dissociation flux, the hydrogen gas density, and the metallicity-dependent H2 formation rate-coefficient and far-UV dust-grain absorption cross-section. We make the distinction between "HI-dust" and "H2-dust" opacity, and we present computations for the "universal H2-dust-limited effective dissociation bandwidth". We validate our analytic formulae with Meudon PDR code computations for the HI-to-H2 density profiles, and total HI column densities. We show that our general 1D formulae predict HI columns and H2 mass fractions that are essentially identical to those found in more complicated (and approximate) spherical (shell/core) models. We apply our theory to compute H2 mass fractions and star-formation thresholds for individual clouds in self-regulated galaxy disks, for a wide range of metallicities. Our formulae for the HI columns and H2 mass fractions may be incorporated into hydrodynamics simulations for galaxy evolution.

Contrasting the direct radiative effect and direct radiative forcing of aerosols

Abstract: . The direct radiative effect (DRE) of aerosols, which is the instantaneous radiative impact of all atmospheric particles on the Earth's energy balance, is sometimes confused with the direct radiative forcing (DRF), which is the change in DRE from pre-industrial to present-day (not including climate feedbacks). In this study we couple a global chemical transport model (GEOS-Chem) with a radiative transfer model (RRTMG) to contrast these concepts. We estimate a global mean all-sky aerosol DRF of −0.36 Wm−2 and a DRE of −1.83 Wm−2 for 2010. Therefore, natural sources of aerosol (here including fire) affect the global energy balance over four times more than do present-day anthropogenic aerosols. If global anthropogenic emissions of aerosols and their precursors continue to decline as projected in recent scenarios due to effective pollution emission controls, the DRF will shrink (−0.22 Wm−2 for 2100). Secondary metrics, like DRE, that quantify temporal changes in both natural and anthropogenic aerosol burdens are therefore needed to quantify the total effect of aerosols on climate.

Observational Signatures of Planets in Protoplanetary Disks I: Gaps Opened by Single and Multiple Young Planets in Disks

Abstract: It has been suggested that the gaps and cavities recently discovered in transitional disks are opened by planets. To explore this scenario, we combine two-dimensional two fluid (gas + particle) hydrodynamical calculations with three-dimensional Monte Carlo Radiative Transfer simulations, and study the observational signatures of gaps opened by one or several planets, making qualitative comparisons with observations. We find that a single planet as small as 0.2 MJ can produce a deep gap at millimeter (mm) wavelengths and almost no features at near-infrared (NIR) wavelengths, while multiple planets can open up a few *10 AU wide common gap at both wavelengths. Both the contrast ratio of the gaps and the wavelength dependence of the gap sizes are broadly consistent with data. We also confirm previous results that NIR gap sizes may be smaller than mm gap sizes due to dust-gas coupling and radiative transfer effects. When viewed at a moderate inclination angle, a physically circular on-centered gap could appear to be off-centered from the star due to shadowing. Planet-induced spiral arms are more apparent at NIR than at mm wavelengths. Overall, our results suggest that the planet-opening-gap scenario is a promising way to explain the origin of the transitional disks. Finally, inspired by the recent ALMA release of the image of the HL Tau disk, we show that multiple narrow gaps, well separated by bright rings, can be opened by 0.2 MJ planets soon after their formation in a relatively massive disk.

Electromagnetic Scattering by Particles and Particle Groups: An Introduction

Abstract: Preface Acknowledgements 1. Introduction 2. The macroscopic Maxwell equations and monochromatic fields 3. Fundamental homogeneous-medium solutions of the macroscopic Maxwell equations 4. Basic theory of frequency-domain electromagnetic scattering by a fixed finite object 5. Far-field scattering 6. The Foldy equations 7. The Stokes parameters 8. Poynting-Stokes tensor 9. Polychromatic electromagnetic fields 10. Polychromatic scattering by fixed and randomly changing objects 11. Measurement of electromagnetic energy flow 12. Measurement of the Stokes parameters 13. Description of far-field scattering in terms of actual optical observables 14. Electromagnetic scattering by a small random group of sparsely distributed particles 15. Statistically isotropic and mirror-symmetric random particles 16. Numerical computations and laboratory measurements of electromagnetic scattering 17. Far-field observables: qualitative and quantitative traits 18. Electromagnetic scattering by discrete random media: far field 19. Near-field scattering by a sparse discrete random medium: microphysical radiative transfer theory 20. Radiative transfer in plane-parallel particulate media 21. Weak localization 22. Epilogue Appendix A. Dyads and dyadics Appendix B. Free-space dyadic Green's function Appendix C. Euler rotation angles Appendix D. Spherical-wave expansion of a plane wave in the far zone Appendix E. Integration quadrature formulas Appendix F. Wigner d-functions Appendix G. Stationary phase evolution of a double integral Appendix H. Hints and answers to selected problems Appendix I. List of acronyms References Index.

Radiative bistability and thermal memory.

Abstract: We predict the existence of a thermal bistability in many-body systems out of thermal equilibrium which exchange heat by thermal radiation using insulator-metal transition materials. We propose a writing-reading procedure and demonstrate the possibility to exploit the thermal bistability to make a volatile thermal memory. We show that this thermal memory can be used to store heat and thermal information (via an encoding temperature) for arbitrary long times. The radiative thermal bistability could find broad applications in the domains of thermal management, information processing, and energy storage.

Accuracy tests of radiation schemes used in hot Jupiter global circulation models

Abstract: The treatment of radiation transport in global circulation models (GCMs) is crucial for correctly describing Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate the rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods’ applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high-temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10% for both heating rates and fluxes using ~10 k -coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~100%, and should thus be used with caution.

Experimental investigation of radiative thermal rectifier using vanadium dioxide

Abstract: Vanadium dioxide (VO2) exhibits a phase-change behavior from the insulating state to the metallic state around 340 K. By using this effect, we experimentally demonstrate a radiative thermal rectifier in the far-field regime with a thin film VO2 deposited on the silicon wafer. A rectification contrast ratio as large as two is accurately obtained by utilizing a one-dimensional steady-state heat flux measurement system. We develop a theoretical model of the thermal rectifier with optical responses of the materials retrieved from the measured mid-infrared reflection spectra, which is cross-checked with experimentally measured heat flux. Furthermore, we tune the operating temperatures by doping the VO2 film with tungsten (W). These results open up prospects in the fields of thermal management and thermal information processing.

Cosmic Reionization on Computers. I. Design and Calibration of Simulations

Abstract: Cosmic Reionization On Computers is a long-term program of numerical simulations of cosmic reionization. Its goal is to model fully self-consistently (albeit not necessarily from the first principles) all relevant physics, from radiative transfer to gas dynamics and star formation, in simulation volumes of up to 100 comoving Mpc, and with spatial resolution approaching 100 pc in physical units. In this method paper, we describe our numerical method, the design of simulations, and the calibration of numerical parameters. Using several sets (ensembles) of simulations in 20 h –1 Mpc and 40 h –1 Mpc boxes with spatial resolution reaching 125 pc at z = 6, we are able to match the observed galaxy UV luminosity functions at all redshifts between 6 and 10, as well as obtain reasonable agreement with the observational measurements of the Gunn-Peterson optical depth at z < 6.

16 yr of RXTE Monitoring of Five Anomalous X-Ray Pulsars

Abstract: We present a summary of the long-term evolution of various properties of the five non-transient anomalous X-ray pulsars (AXPs) 1E 1841–045, RXS J170849.0–400910, 1E 2259+586, 4U 0142+61, and 1E 1048.1–5937, regularly monitored with RXTE from 1996 to 2012. We focus on three properties of these sources: the evolution of the timing, pulsed flux, and pulse profile. We report several new timing anomalies and radiative events, including a putative anti-glitch seen in 1E 2259+586 in 2009, and a second epoch of very large spin-down rate fluctuations in 1E 1048.1–5937 following a large flux outburst. We compile the properties of the 11 glitches and 4 glitch candidates observed from these 5 AXPs between 1996 and 2012. Overall, these monitoring observations reveal several apparent patterns in the behavior of this sample of AXPs: large radiative changes in AXPs (including long-lived flux enhancements, short bursts, and pulse profile changes) are rare, occurring typically only every few years per source; large radiative changes are almost always accompanied by some form of timing anomaly, usually a spin-up glitch; only 20%-30% of timing anomalies are accompanied by any form of radiative change. We find that AXP radiative behavior at the times of radiatively loud glitches is sufficiently similar to suggest common physical origins. The similarity in glitch properties when comparing radiatively loud and radiatively silent glitches in AXPs suggests a common physical origin in the stellar interior. Finally, the overall similarity of AXP and radio pulsar glitches suggests a common physical origin for both phenomena.

Radiative-convective instability

Abstract: [1] Radiative-moist-convective equilibrium (RCE) is a simple paradigm for the statistical equilibrium the earth's climate would exhibit in the absence of lateral energy transport. It has generally been assumed that for a given solar forcing and long-lived greenhouse gas concentration, such a state would be unique, but recent work suggests that more than one stable equilibrium may be possible. Here we show that above a critical specified sea surface temperature, the ordinary RCE state becomes linearly unstable to large-scale overturning circulations. The instability migrates the RCE state toward one of the two stable equilibria first found by Raymond and Zeng (2000). It occurs when the clear-sky infrared opacity of the lower troposphere becomes so large, owing to high water vapor concentration, that variations of the radiative cooling of the lower troposphere are governed principally by variations in upper tropospheric water vapor. We show that the instability represents a subcritical bifurcation of the ordinary RCE state, leading to either a dry state with large-scale descent, or to a moist state with mean ascent; these states may be accessed by finite amplitude perturbations to ordinary RCE in the subcritical state, or spontaneously in the supercritical state. As first suggested by Raymond (2000) and Sobel et al. (2007), the latter corresponds to the phenomenon of self-aggregation of moist convection, taking the form of cloud clusters or tropical cyclones. We argue that the nonrobustness of self-aggregation in cloud system resolving models may be an artifact of running such models close to the critical temperature for instability.

The nebular spectra of SN 2012aw and constraints on stellar nucleosynthesis from oxygen emission lines

Abstract: We present nebular-phase optical and near-infrared spectroscopy of the Type IIP supernova SN 2012aw combined with non-local thermodynamic equilibrium radiative transfer calculations applied to ejecta from stellar evolution/explosion models Our spectral synthesis models generally show good agreement with the ejecta from a M-ZAMS = 15 M-circle dot progenitor star The emission lines of oxygen, sodium, and magnesium are all consistent with the nucleosynthesis in a progenitor in the 14-18 M-circle dot range We also demonstrate how the evolution of the oxygen cooling lines of [O i] lambda 5577, [O i] lambda 6300, and [O i] lambda 6364 can be used to constrain the mass of oxygen in the non-molecularly cooled ashes to 20 M-circle dot progenitor

A non-grey analytical model for irradiated atmospheres - I. Derivation

Abstract: Context. Semi-grey atmospheric models (with one opacity for the visible and one opacity for the infrared) are useful for understanding the global structure of irradiated atmospheres, their dynamics, and the interior structure and evolution of planets, brown dwarfs, and stars. When compared to direct numerical radiative transfer calculations for irradiated exoplanets, however, these models systematically overestimate the temperatures at low optical depths, independently of the opacity parameters.Aims. We investigate why semi-grey models fail at low optical depths and provide a more accurate approximation to the atmospheric structure by accounting for the variable opacity in the infrared.Methods. Using the Eddington approximation, we derive an analytical model to account for lines and/or bands in the infrared. Four parameters (instead of two for the semi-grey models) are used: a visible opacity (κ v ), two infrared opacities, (κ 1 and κ 2 ), and β (the fraction of the energy in the beam with opacities κ 1 ). We consider that the atmosphere receives an incident irradiation in the visible with an effective temperature T irr and at an angle μ ∗ , and that it is heated from below with an effective temperature T int .Results. Our non-grey, irradiated line model is found to provide a range of temperatures that is consistent with that obtained by numerical calculations. We find that if the stellar flux is absorbed at optical depth larger than τ lim = (κ R /κ 1 κ 2 )(κ R κ P /3)1/2 , it is mainly transported by the channel of lowest opacity whereas if it is absorbed at τ ≳ τ lim it is mainly transported by the channel of highest opacity, independently of the spectral width of those channels. For low values of β (expected when lines are dominant), we find that the non-grey effects significantly cool the upper atmosphere. However, for β ≳ 1/2 (appropriate in the presence of bands with a wavelength-dependence smaller than or comparable to the width of the Planck function), we find that the temperature structure is affected down to infrared optical depths unity and deeper as a result of the so-called blanketing effect.Conclusions. The expressions that we derive can be used to provide a proper functional form for algorithms that invert the atmospheric properties from spectral information. Because a full atmospheric structure can be calculated directly, these expressions should be useful for simulations of the dynamics of these atmospheres and of the thermal evolution of the planets. Finally, they should be used to test full radiative transfer models and to improve their convergence.

Seasonal and Elevational Variations of Black Carbon and Dust in Snow and Ice in the Solu-Khumbu, Nepal and Estimated Radiative Forcings

Abstract: Abstract. Black carbon (BC) and dust deposited on snow and glacier surfaces can reduce the surface albedo, accelerate snow and ice melt, and trigger albedo feedback. Assessing BC and dust concentrations in snow and ice in the Himalaya is of interest because this region borders large BC and dust sources, and seasonal snow and glacier ice in this region are an important source of water resources. Snow and ice samples were collected from crevasse profiles and snow pits at elevations between 5400 and 6400 m a.s.l. from Mera glacier located in the Solu-Khumbu region of Nepal during spring and fall 2009, providing the first observational data of BC concentrations in snow and ice from the southern slope of the Himalaya. The samples were measured for Fe concentrations (used as a dust proxy) via ICP-MS, total impurity content gravimetrically, and BC concentrations using a Single Particle Soot Photometer (SP2). Measured BC concentrations underestimate actual BC concentrations due to changes to the sample during storage and loss of BC particles in the ultrasonic nebulizer; thus, we correct for the underestimated BC mass. BC and Fe concentrations are substantially higher at elevations

Vertical shear instability in accretion disc models with radiation transport

Abstract: The origin of turbulence in accretion discs is still not fully understood. While the magneto-rotational instability is considered to operate in sufficiently ionized discs, its role in the poorly ionized protoplanetary disc is questionable. Recently, the vertical shear instability (VSI) has been suggested as a possible alternative. Our goal is to study the characteristics of this instability and the efficiency of angular momentum transport, in extended discs, under the influence of radiative transport and irradiation from the central star. We use multi-dimensional hydrodynamic simulations to model a larger section of an accretion disc. First we study inviscid and weakly viscous discs using a fixed radial temperature profile in two and three spatial dimensions. The simulations are then extended to include radiative transport and irradiation from the central star. In agreement with previous studies we find for the isothermal disc a sustained unstable state with a weak positive angular momentum transport of the order of $\alpha \approx 10^{-4}$. Under the inclusion of radiative transport the disc cools off and the turbulence terminates. For discs irradiated from the central star we find again a persistent instability with a similar $\alpha$ value as for the isothermal case. We find that the VSI can indeed generate sustained turbulence in discs albeit at a relatively low level with $\alpha$ about few times $10^{-4}$

Reconstructing the density and temperature structure of prestellar cores from Herschel data: A case study for B68 and L1689B

Abstract: Utilizing multiwavelength dust emission maps acquired with Herschel, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using an inverse-Abel transform-based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight, which were previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in the core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 ± 0.5 K and 9.8 ± 0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M⊙ and 11 M⊙, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, κλ = 0.1 × (λ/300 μm)-2 cm2 g-1 agrees to within 50% with the dust extinction constraints.