Transport of ionizing radiation in terrestrial-like exoplanet atmospheres

The propagation of ionizing radiation through model atmospheres of terrestrial-like exoplanets is studied for a large range of column densities and incident photon energies using a Monte Carlo code we have developed to treat Compton scattering and photoabsorption. Incident spectra from parent star flares, supernovae, and gamma-ray bursts are modeled and compared to energetic particles in importance. Large irradiation events with fluences of 10⁶-10⁹ erg cm⁻² at the conventional habitable zone can occur at a rate from many per day (flares from young low-mass parent stars) to ∼100 per Gyr (supernovae and gamma-ray bursts). We find that terrestrial-like exoplanets with atmospheres thinner than about 100 g cm⁻² block nearly all X-rays, but transmit and reprocess a significant fraction of incident γ-rays, producing a characteristic, flat surficial spectrum. Thick atmospheres (?100 g cm⁻²) efficiently block even γ-rays, but nearly all the incident energy is redistributed into diffuse UV and visible aurora-like emission, increasing the effective atmospheric transmission by many orders of magnitude. Depending on the presence of molecular UV absorbers and atmospheric thickness, up to 10% of the incident energy can reach the surface as UV reemission. For the Earth, between 2×10⁻³ and 4×10⁻² of the incident flux reaches the ground in the biologically effective 200-320 nm range, depending on O₂/O₃ shielding. For atmospheres thicker than ∼50 g cm⁻² in the case of pure Rayleigh scattering and ∼100 g cm⁻² in the case of O₂/O₃ absorption, the UV reemission exceeds the surficial transmitted ionizing radiation. We also discuss the effects of angle of incidence and derive a modified two-stream approximation solution for the UV transfer. Finally, we suggest that transient atmospheric ionization layers can be frequently created at altitudes lower than the equilibrium layers that result from steady irradiation and winds from the parent star. We suggest that these events can produce frequent fluctuations in atmospheric ionization levels and surficial UV fluxes on terrestrial-like planets.     

Axisymmetry in protoplanetary nebulae – II. A near-infrared imaging polarimetric survey

In an imaging polarimetry survey of candidate post-AGB stars, scattering envelopes are detected around 20 objects. These envelopes represent the final mass-loss phases at the end of the AGB. In all cases, evidence for axisymmetry in the dust density is seen, suggesting that the presence of an axisymmetric outflow may be a ubiquitous phenomenon of the AGB to post-AGB transition. We use the polarized flux images to classify the objects into detached shell, bipolar and unresolved targets. Modelling based on a simple axisymmetric shell geometry supports the contention that post-AGB objects fall into one of two classes that differ in the amount of dust in the circumstellar environment: the detached shells correspond to stars with an optically thin expanding circumstellar envelope (CSE) whereas the bipolar and unresolved targets have optically thick dust structures, probably in the form of discs, which remain bound to the star, rather than partaking in the expansion of the AGB CSE. It is suggested that this bifurcation in morphology is rooted in the presence or absence of a binary companion, which determines whether or not a disc forms. Because the detached shell objects also appear axisymmetric, an additional mechanism for generating the axisymmetry, such as a magnetically shaped outflow, is needed if they do indeed have single star progenitors.     

Cosmic microwave background polarization due to scattering in clusters

Scattering of the cosmic microwave background (CMB) in clusters of galaxies polarizes the radiation. We explore several polarization components which have their origin in the kinematic quadrupole moments induced by the motion of the scattering electrons, either directed or random. Polarization levels and patterns are determined in a cluster simulated by the hydrodynamical enzo code. We find that polarization signals can be as high as  ∼1 μK  , a level that may be detectable by upcoming CMB experiments.     

Generalized Born scattering of elastic waves in 3-D media

It is well known that when a seismic wave propagates through an elastic medium with gradients in the parameters which describe it (e.g. slowness and density), energy is scattered from the incident wave generating low-frequency partial reflections. Many approximate solutions to the wave equation, e.g. geometrical ray theory (GRT), Maslov theory and Gaussian beams, do not model these signals. The problem of describing partial reflections in 1-D media has been extensively studied in the seismic literature and considerable progress has been made using iterative techniques based on WKBJ, Airy or Langer type ansätze. In this paper we derive a first-order scattering formalism to describe partial reflections in 3-D media. The correction term describing the scattered energy is developed as a volume integral over terms dependent upon the first spatial derivatives (gradients) of the parameters describing the medium and the solution. The relationship we derive could, in principle, be used as the basis for an iterative scheme but the computational expense, particularly for elastic media, will usually prohibit this approach. The result we obtain is closely related to the usual Born approximation, but differs in that the scattering term is not derived from a perturbation to a background model, but rather from the error in an approximate Green's function. We examine analytically the relationship between the results produced by the new formalism and the usual Born approximation for a medium which has no long-wavelength heterogeneities. We show that in such a case the two methods agree approximately as expected, but that in a media with heterogeneities of all wavelengths the new gradient scattering formalism is superior. We establish analytically the connection between the formalism developed here and the iterative approach based on the WKBJ solution which has been used previously in 1-D media. Numerical examples are shown to illustrate the examples discussed.     

Generalized Born inversion of seismic reflection data

Summary . Born inverse methods give accurate and stable results when the source wavelet is impulsive. However, in many practical applications (reflection seismology) an impulsive source cannot be realized and the inversion needs to be generalized to include an arbitrary source function. In this paper, we present a Born solution to the seismic inverse problem which can accommodate an arbitrary source function and give accurate and stable results. It is shown that the form of the generalized inversion algorithm reduces to a Wiener shaping ***filter, which is solved efficiently using a Levinson recursion algorithm. Numerical examples of synthetic and real field data illustrate the validity of our method.     

Radiative Heat Transfer and Hydrostatic Stability in Nocturnal Fog

We have performed a one-dimensional and transient radiative heat transfer analysis in order to investigate interaction between atmospheric radiation and convective instability within a nocturnal fog. The radiation element method using the Ray Emission Model (REM²), which is a generalized numerical method, in conjunction with a line-by-line (LBL) method, is employed to attain high spectral resolution calculations for anisotropically scattering fog. The results show that the convective instability has a strong dependence on radiative properties of the fog. For the condition of a 20-m droplet diameter and liquid water content of 0.1 × 10^–3 kg m^–3;, the temperature profile within the fog becomes S shaped, and a convective instability layer forms in the middle or lower level of the fog. However, for the same water content and a 40-m diameter droplet, no strong convective instability layer forms, whereas for a 10-m diameter droplet a strong convective instability is observed.