Three-dimensional radiative transfer in an anisotropically scattering, semi-infinite medium: generalized reflection function

Three-dimensional radiative transfer in an anisotropic scattering medium exposed to spatially varying, collimated radiation is studied. The generalized reflection function for a semi-infinite medium with a very general scattering phase function is the focus of this investigation. An integral transform is used to reduce the three-dimensional transport equation to a one-dimensional form, and a modified Ambarzumian's method is applied to formulate a nonlinear integral equation for the generalized reflection function. The integration is over both the polar and azimuthal angles; hence, the integral equation is said to be in the double-integral form. The double-integral, reflection function formulation can handle a variety of anisotropic phase functions and does not require an expansion of the phase function in a Legendre polynomial series. Complicated kernel transformations of previous single-integral studies are eliminated. Single and double scattering approximations are developed. Numerical results are presented for a Rayleigh phase function to illustrate the computational characteristics of the method and are compared to results obtained with the single-integral method. Agreement between the two approaches is excellent; however, as the transform variable increases beyond five the number of quadrature points required for the double-integral method to produce accurate solutions significantly increases. A new interpolation scheme produces accurate results when the transform variable is large.     

Comparisons between model equations for long waves

Summary Considered here are model equations for weakly nonlinear and dispersive long waves, which feature general forms of dispersion and pure power nonlinearity. Two variants of such equations are introduced, one of Korteweg-de Vries type and one of regularized long-wave type. It is proven that solutions of the pure initial-value problem for these two types of model equations are the same, to within the order of accuracy attributable to either, on the long time scale during which nonlinear and dispersive effects may accumulate to make an order-one relative difference to the wave profiles.This research was supported in part by the National Science Foundation. A considerable portion of the project was completed while the first author was resident at the Institute for Mathematics and Its Applications, University of Minnesota.     

Truncated Newton method for sparse unconstrained optimization using automatic differentiation

When solving large complex optimization problems, the user is faced with three major problems. These are (i) the cost in human time in obtaining accurate expressions for the derivatives involved; (ii) the need to store second derivative information; and (iii), of lessening importance, the time taken to solve the problem on the computer. For many problems, a significant part of the latter can be attributed to solving Newton-like equations. In the algorithm described, the equations are solved using a conjugate direction method that only needs the Hessian at the current point when it is multiplied by a trial vector. In this paper, we present a method that finds this product using automatic differentiation while only requiring vector storage. The method takes advantage of any sparsity in the Hessian matrix and computes exact derivatives. It avoids the complexity of symbolic differentiation, the inaccuracy of numerical differentiation, the labor of finding analytic derivatives, and the need for matrix store. When far from a minimum, an accurate solution to the Newton equations is not justified, so an approximate solution is obtained by using a version of Dembo and Steihaug's truncated Newton algorithm (Ref. 1).This paper was presented at the SIAM National Meeting, Boston, Massachusetts, 1986.     

Exact soliton solutions of spinor Bose-Einstein condensates

We study matter-wave solitons in Bose-Einstein condensates of ultracold gaseous atoms with spin degrees of freedom and present a class of exact solutions based on the inverse scattering method. The one-soliton solutions are classified with respect to the spin states. We analyze collisional effects between solitons in the same or different spin state(s), which reveals a very interesting possibility: we can manipulate the spin dynamics by controlling the parameters of colliding solitons.     

Acoustic field in a closed layered shell with resonant systems

The acoustic field inside a shell excited by a spatially inhomogeneous harmonic pressure field is studied. The shell is assumed to have a finite length, a set of orthogonal stiffening ribs, two ends bounding the acoustic volume, and a sound-insulating structure, which includes layers of sound-insulating material, resonant elements, and an interior panel. The shell is considered to be orthotropic with boundary conditions corresponding to a free support. For the acoustic field in the closed volume, analytical expressions are derived with allowance for the elastoacoustic interaction of the shell with the sound-insulating layers and with the medium both inside the shell (with arbitrary impedance values at the ends) and outside it. These expressions are used to investigate the effect of different types of resonant systems on the sound field inside the shell.     

3D-TOFSIMS characterization of black spots in polymer light emitting diodes

The occurrence and formation of black spots areas in PolyLED devices has been studied by time-of-flight SIMS (TOFSIMS). The composition, shape and position of the black spots is visualised by three-dimensional (3D)-TOFSIMS depth-profiling. It has been established that the formation of non-emissive spots is due to the growth of aluminium oxide clusters at the AlBa/polymer interface. Electron injection in the black spots is lost by the resulting local increase of the resistivity of the cathode.     

Specific features of photoisomerization of provitamin D 3 in a nematic liquid crystal

Photoisomerization of provitamin D ₃ (7-dehydrocholesterol) in a nematic liquid crystal (ZLI-1695, Merck) is investigated in detail by UV absorption spectroscopy. It is found that dissolution of chiral molecules of provitamin D ₃ induces the cholesteric phase in a nematic. The spectral kinetics of photoisomerization in this phase changes significantly from that in an ethanol solution. A sharp nonmonotonic dependence of the increase in accumulation of trans isomer tachysterol in a liquid crystal matrix with a decrease in the induced cholesteric pitch from 2200 to 25 μm is revealed.     

Effect of hydrogen on the photoelectronic properties of GaAs/InGaAs quantum-well heterostructures with an island palladium layer on the surface

The effect of hydrogen on the photoluminescence and planar conductivity of GaAs/InGaAs quantum-well heterostructures with an island Pd layer at the anodically oxidized surface was studied. Unlike continuous deposited Pd layers, island layers do not cause the formation of defects in the GaAs surface region and yet the Pd layer maintains high catalytic activity with respect to hydrogen. It is found that the thermal treatment of such a structure in a hydrogen atmosphere causes atomic-hydrogen passivation of the defects in quantum wells. Studies of the characteristics of planar photoresistors with an island Pd layer acting as hydrogen sensors show that their hydrogen detectivity is approximately two orders of magnitude higher than that of diode structures with continuous Pd layers.     

Rehabilitation of the Gauss-Jordan algorithm

Summary In this paper a Gauss-Jordan algorithm with column interchanges is presented and analysed. We show that, in contrast with Gaussian elimination, the Gauss-Jordan algorithm has essentially differing properties when using column interchanges instead of row interchanges for improving the numerical stability. For solutions obtained by Gauss-Jordan with column interchanges, a more satisfactory bound for the residual norm can be given. The analysis gives theoretical evidence that the algorithm yields numerical solutions as good as those obtained by Gaussian elimination and that, in most practical situations, the residuals are equally small. This is confirmed by numerical experiments. Moreover, timing experiments on a Cyber 205 vector computer show that the algorithm presented has good vectorisation properties.