Parallel Implementation of the PHOENIX Generalized Stellar Atmosphere Program. II: Wavelength Parallelization

We describe an important addition to the parallel implementation of our generalized NLTE stellar atmosphere and radiative transfer computer program PHOENIX. In a previous paper in this series we described data and task parallel algorithms we have developed for radiative transfer, spectral line opacity, and NLTE opacity and rate calculations. These algorithms divided the work spatially or by spectral lines, that is distributing the radial zones, individual spectral lines, or characteristic rays among different processors and employ, in addition task parallelism for logically independent functions (such as atomic and molecular line opacities). For finite, monotonic velocity fields, the radiative transfer equation is an initial value problem in wavelength, and hence each wavelength point depends upon the previous one. However, for sophisticated NLTE models of both static and moving atmospheres needed to accurately describe, e.g., novae and supernovae, the number of wavelength points is very large (200,000--300,000) and hence parallelization over wavelength can lead both to considerable speedup in calculation time and the ability to make use of the aggregate memory available on massively parallel supercomputers. Here, we describe an implementation of a pipelined design for the wavelength parallelization of PHOENIX, where the necessary data from the processor working on a previous wavelength point is sent to the processor working on the succeeding wavelength point as soon as it is known. Our implementation uses a MIMD design based on a relatively small number of standard MPI library calls and is fully portable between serial and parallel computers.Comment: AAS-TeX, 15 pages, full text with figures available at ftp://calvin.physast.uga.edu/pub/preprints/Wavelength-Parallel.ps.gz ApJ, in pres

Kinematic and morphological modeling of the bipolar nebula Sa2-237

We present [OIII]500.7nm and Halpha+[NII] images and long-slit, high resolution echelle spectra in the same spectral regions of Sa2--237, a possible bipolar planetary nebula. The image shows a bipolar nebula of about 34" extent, with a narrow waist, and showing strong point symmetry about the central object, indicating it's likely binary nature. The long slit spectra were taken over the long axis of the nebula, and show a distinct ``eight'' shaped pattern in the velocity--space plot, and a maximum projected outflow velocity of V=106km/s, both typical of expanding bipolar planetary nebulae. By model fitting the shape and spectrum of the nebula simultaneously, we derive the inclination of the long axis to be 70 degrees, and the maximum space velocity of expansion to be 308 km/s. Due to asymmetries in the velocities we adopt a new value for the system's heliocentric radial velocity of -30km/s. We use the IRAS and 21cm radio fluxes, the energy distribution, and the projected size of Sa2-237 to estimate it's distance to be 2.1+-0.37kpc. At this distance Sa2-237 has a luminosity of 340 Lsun, a size of 0.37pc, and -- assuming constant expansion velocity -- a nebular age of 624 years. The above radial velocity and distance place Sa2--237 in the disk of the Galaxy at z=255pc, albeit with somewhat peculiar kinematics.Comment: 10pp, 4 fig

Nature of the Earth's Interior

n/

Kaluza-Klein electrically charged black branes in M-theory

We present a class of Kaluza-Klein electrically charged black p-brane solutions of ten-dimensional, type IIA superstring theory. Uplifting to eleven dimensions these solutions are studied in the context of M-theory. They can be interpreted either as a p+1 extended object trapped around the eleventh dimension along which momentum is flowing or as a boost of the following backgrounds: the Schwarzschild black (p+1)-brane or the product of the (10-p)-dimensional Euclidean Schwarzschild manifold with the (p+1)-dimensional Minkowski spacetime.Comment: 16 pages, uses latex and epsf macro, figures include

Relativistic J-matrix method

The relativistic version of the J-matrix method for a scattering problem on the potential vanishing faster than the Coulomb one is formulated. As in the non-relativistic case it leads to a finite algebraic eigenvalue problem. The derived expression for the tangent of phase shift is simply related to the non-relativistic case formula and gives the latter as a limit case. It is due to the fact that the used basis set satisfies the ``kinetic balance condition''.Comment: 21 pages, RevTeX, accepted for publication in Phys. Rev.

Abelian Duality

We show that on three-dimensional Riemannian manifolds without boundaries and with trivial first real de Rham cohomology group (and in no other dimensions) scalar field theory and Maxwell theory are equivalent: the ratio of the partition functions is given by the Ray-Singer torsion of the manifold. On the level of interaction with external currents, the equivalence persists provided there is a fixed relation between the charges and the currents.Comment: 11 pages, LaTeX, no figures, a reference added, submitted to Phys. Rev.

3-D Photoionization Structure and Distances of Planetary Nebulae II. Menzel 1

We present the results of a spatio-kinematic study of the planetary nebula Menzel 1 using spectro-photometric mapping and a 3-D photoionization code. We create several 2-D emission line images from our long-slit spectra, and use these to derive the line fluxes for 15 lines, the Halpha/Hbeta extinction map, and the [SII] line ratio density map of the nebula. We use our photoionization code constrained by these data to derive the three-dimensional nebular structure and ionizing star parameters of Menzel 1 by simultaneously fitting the integrated line intensities, the density map, and the observed morphologies in several lines, as well as the velocity structure. Using theoretical evolutionary tracks of intermediate and low mass stars, we derive a mass for the central star of 0.63+-0.05 Msolar. We also derive a distance of 1050+_150 pc to Menzel 1.Comment: To be published in ApJ of 10th February 2005. 12 figure

Effect of Poisson ratio on cellular structure formation

Mechanically active cells in soft media act as force dipoles. The resulting elastic interactions are long-ranged and favor the formation of strings. We show analytically that due to screening, the effective interaction between strings decays exponentially, with a decay length determined only by geometry. Both for disordered and ordered arrangements of cells, we predict novel phase transitions from paraelastic to ferroelastic and anti-ferroelastic phases as a function of Poisson ratio.Comment: 4 pages, Revtex, 4 Postscript figures include