Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

A Bootstrap Test for the Probability of Ruin in the Compound Poisson Risk Process

In this article we propose a bootstrap test for the probability of ruin in the compound Poisson risk process. We adopt the P-value approach, which leads to a more complete assessment of the underlying risk than the probability of ruin alone. We provide second-order accurate P-values for this testing problem and consider both parametric and nonparametric estimators of the individual claim amount distribution. Simulation studies show that the suggested bootstrap P-values are very accurate and outperform their analogues based on the asymptotic normal approximation.

Asphyxial suicide by inhalation of helium inside a plastic bag

Asphyxial suicide by placing a plastic bag over the head, especially in combination with inhalation of gases, is a rarely described method of committing suicide. This article reports a case of suicidal asphyxiation by inhaling the inert gas helium inside a plastic bag. A 64-year-old man probably followed the instructions described in an article about committing suicide written by a medical practitioner from Zürich. This form of suicide is recommended by right-to-die groups and in the internet as a certain, fast, and painless suicide method. Additionally, it leaves only seldom externally visible marks or pathomorphological findings on the body. If the plastic bag and other auxiliary means are removed by another person, the forensic death investigation of cause and manner of death may be very difficult. Therefore, the death scene investigation and the inquiry ordered in the environment of the deceased are very important.

Accurate computations of the structures and binding energies of the imidazole ... benzene and pyrrole ... benzene complexes

Using explicitly-correlated coupled-cluster theory with single and double excitations, the intermolecular distances and interaction energies of the T-shaped imidazole⋯⋯benzene and pyrrole⋯⋯benzene complexes have been computed in a large augmented correlation-consistent quadruple-zeta basis set, adding also corrections for connected triple excitations and remaining basis-set-superposition errors. The results of these computations are used to assess other methods such as Møller–Plesset perturbation theory (MP2), spin-component-scaled MP2 theory, dispersion-weighted MP2 theory, interference-corrected explicitly-correlated MP2 theory, dispersion-corrected double-hybrid density-functional theory (DFT), DFT-based symmetry-adapted perturbation theory, the random-phase approximation, explicitly-correlated ring-coupled-cluster-doubles theory, and double-hybrid DFT with a correlation energy computed in the random-phase approximation.

Line focusing for soft x-ray laser-plasma lasing

A computational study of line-focus generation was done using a self-written ray-tracing code and compared to experimental data. Two line-focusing geometries were compared, i.e., either exploiting the sagittal astigmatism of a tilted spherical mirror or using the spherical aberration of an off-axis- illuminated spherical mirror. Line focusing by means of astigmatism or spherical aberration showed identical results as expected for the equivalence of the two frames of reference. The variation of the incidence angle on the target affects the line-focus length, which affects the amplification length such that as long as the irradiance is above the amplification threshold, it is advantageous to have a longer line focus. The amplification threshold is physically dependent on operating parameters and plasma-column conditions and in the present study addresses four possible cases.

Comparison of six radiographic projections to assess femoral head/neck asphericity

Early radiographic detection of femoroacetabular impingement might prevent initiation and progression of osteoarthritis. The structural abnormality in femoral-induced femoroacetabular impingement (cam type) is frequently asphericity at the anterosuperior head/neck contour. To determine which of six radiographic projections (anteroposterior, Dunn, Dunn/45 degrees flexion, cross-table/15 degrees internal rotation, cross-table/neutral rotation, and cross-table/15 degrees external rotation) best identifies femoral head/neck asphericity, we studied 21 desiccated femurs; 11 with an aspherical femoral head/neck contour and 10 with a spherical femoral head/neck contour. To radiographically quantify femoral head asphericity, we measured the angle where the femoral head/neck leaves sphericity (angle alpha). The aspherical femoral head/neck contours had a greater maximum angle alpha (70 degrees ) compared with the spherical head/neck contours (50 degrees ). The angle alpha varied depending on the radiographic projection: it was greatest in the Dunn view with 45 degrees hip flexion (71 degrees +/- 10 degrees ) and least in the cross-table view in 15 degrees external rotation (51 degrees +/- 7 degrees ). Diagnosis of a pathologic femoral head/neck contour depends on the radiologic projection. The Dunn view in 45 degrees or 90 degrees flexion or a cross-table projection in internal rotation best show femoral head/neck asphericity, whereas anteroposterior or externally rotated cross-table views are likely to miss asphericity. Level of Evidence: Prognostic study, level II-1 (retrospective study).

Frequency-domain source localization shows state-dependent diazepam effects in 47-channel EEG

The topic of this study was to evaluate state-dependent effects of diazepam on the frequency characteristics of 47-channel spontaneous EEG maps. A novel method, the FFT-Dipole-Approximation (Lehmann and Michel, 1990), was used to study effects on the strength and the topography of the maps in the different frequency bands. Map topography was characterized by the 3-dimensional location of the equivalent dipole source and map strength was defined as the spatial standard deviation (the Global Field Power) of the maps of each frequency point. The Global Field Power can be considered as a measure of the amount of energy produced by the system, while the source location gives an estimate of the center of gravity of all sources in the brain that were active at a certain frequency. State-dependency was studied by evaluating the drug effects before and after a continuous performance task of 25 min duration. Clear interactions between drug (diazepam vs. placebo) and time after drug intake (before and after the task) were found, especially in the inferior-superior location of the dipole sources. It supports the hypothesis that diazepam, like other drugs, has different effects on brain functions depending on the momentary functional state of the brain. In addition to the drug effects, clearly different source locations and Global Field Power were found for the different frequency bands, replicating earlier reports (Michel et al., 1992).

Smell and taste of chewing gum affect frequency domain EEG source localizations

We investigated brain electric field signatures of subjective feelings after chewing regular gum or gum base without flavor. 19-channel eyes-closed EEG from 20 healthy males before and after 5 minutes of chewing the two gum types in random sequence was source modeled in the frequency domain using the FFT-Dipole-Approximation. 3-dimensional brain locations and strengths (Global Field Power, GFP) of the equivalent sources of five frequency bands were computed as changes from pre-chewing baseline. Gum types differed (ANOVA) in pre-post changes of source locations for the alpha-2 band (to anterior and right after regular gum, opposite after gum base) and beta-2 band (to anterior and inferior after regular gum, opposite after gum base), and of GFP for delta-theta, alpha-2 and beta-1 (regular gum: increase, gum base: decrease). Subjective feeling changed to more positive values after regular gum than gum base (ANOVA).—Thus, chewing gum with and without taste-smell activates different brain neuronal populations.

Thermal 2-loop master spectral function at finite momentum

When considering NLO corrections to thermal particle production in the “relativistic” regime, in which the invariant mass squared of the produced particle is K2 ~ (πT)2, then the production rate can be expressed as a sum of a few universal “master” spectral functions. Taking the most complicated 2-loop master as an example, a general strategy for obtaining a convergent 2-dimensional integral representation is suggested. The analysis applies both to bosonic and fermionic statistics, and shows that for this master the non-relativistic approximation is only accurate for K2 ~(8πT)2, whereas the zero-momentum approximation works surprisingly well. Once the simpler masters have been similarly resolved, NLO results for quantities such as the right-handed neutrino production rate from a Standard Model plasma or the dilepton production rate from a QCD plasma can be assembled for K2 ~ (πT)2.

From transient to steady state deformation and grain size: A thermodynamic approach using elasto-visco-plastic numerical modeling

Numerical simulation experiments give insight into the evolving energy partitioning during high-strain torsion experiments of calcite. Our numerical experiments are designed to derive a generic macroscopic grain size sensitive flow law capable of describing the full evolution from the transient regime to steady state. The transient regime is crucial for understanding the importance of micro structural processes that may lead to strain localization phenomena in deforming materials. This is particularly important in geological and geodynamic applications where the phenomenon of strain localization happens outside the time frame that can be observed under controlled laboratory conditions. Ourmethod is based on an extension of the paleowattmeter approach to the transient regime. We add an empirical hardening law using the Ramberg-Osgood approximation and assess the experiments by an evolution test function of stored over dissipated energy (lambda factor). Parameter studies of, strain hardening, dislocation creep parameter, strain rates, temperature, and lambda factor as well asmesh sensitivity are presented to explore the sensitivity of the newly derived transient/steady state flow law. Our analysis can be seen as one of the first steps in a hybrid computational-laboratory-field modeling workflow. The analysis could be improved through independent verifications by thermographic analysis in physical laboratory experiments to independently assess lambda factor evolution under laboratory conditions.

Convolution roots and differentiability of isotropic positive definite functions on spheres

We prove that any isotropic positive definite function on the sphere can be written as the spherical self-convolution of an isotropic real-valued function. It is known that isotropic positive definite functions on d-dimensional Euclidean space admit a continuous derivative of order [(d − 1)/2]. We show that the same holds true for isotropic positive definite functions on spheres and prove that this result is optimal for all odd dimensions.

Saddlepoint approximations to sensitivities of tail probabilities of random sums and comparisons with Monte Carlo estimators

This article proposes computing sensitivities of upper tail probabilities of random sums by the saddlepoint approximation. The considered sensitivity is the derivative of the upper tail probability with respect to the parameter of the summation index distribution. Random sums with Poisson or Geometric distributed summation indices and Gamma or Weibull distributed summands are considered. The score method with importance sampling is considered as an alternative approximation. Numerical studies show that the saddlepoint approximation and the method of score with importance sampling are very accurate. But the saddlepoint approximation is substantially faster than the score method with importance sampling. Thus, the suggested saddlepoint approximation can be conveniently used in various scientific problems.

A new, fast and flexible radiative transfer method for Venus general circulation models

We present a new radiation scheme for the Oxford Planetary Unified Model System for Venus, suitable for the solar and thermal bands. This new and fast radiative parameterization uses a different approach in the two main radiative wavelength bands: solar radiation (0.1-5.5 mu m) and thermal radiation (1.7-260 mu m). The solar radiation calculation is based on the delta-Eddington approximation (two-stream-type) with an adding layer method. For the thermal radiation case, a code based on an absorptivity/emissivity formulation is used. The new radiative transfer formulation implemented is intended to be computationally light, to allow its incorporation in 3D global circulation models, but still allowing for the calculation of the effect of atmospheric conditions on radiative fluxes. This will allow us to investigate the dynamical-radiative-microphysical feedbacks. The model flexibility can be also used to explore the uncertainties in the Venus atmosphere such as the optical properties in the deep atmosphere or cloud amount. The results of radiative cooling and heating rates and the global-mean radiative-convective equilibrium temperature profiles for different atmospheric conditions are presented and discussed. This new scheme works in an atmospheric column and can be easily implemented in 3D Venus global circulation models. (C) 2014 Elsevier Ltd. All rights reserved.

Self-adjoint extensions for confined electrons: From a particle in a spherical cavity to the hydrogen atom in a sphere and on a cone

In a recent study of the self-adjoint extensions of the Hamiltonian of a particle confined to a finite region of space, in which we generalized the Heisenberg uncertainty relation to a finite volume, we encountered bound states localized at the wall of the cavity. In this paper, we study this situation in detail both for a free particle and for a hydrogen atom centered in a spherical cavity. For appropriate values of the self-adjoint extension parameter, the bound states localized at the wall resonate with the standard hydrogen bound states. We also examine the accidental symmetry generated by the Runge–Lenz vector, which is explicitly broken in a spherical cavity with general Robin boundary conditions. However, for specific radii of the confining sphere, a remnant of the accidental symmetry persists. The same is true for an electron moving on the surface of a finite circular cone, bound to its tip by a 1/r1/r potential.