The GPRIME approach to finite element modeling

GPRIME, an interactive modeling system, runs on the CDC 6000 computers and the DEC VAX 11/780 minicomputer. This system includes three components: (1) GPRIME, a user friendly geometric language and a processor to translate that language into geometric entities, (2) GGEN, an interactive data generator for 2-D models; and (3) SOLIDGEN, a 3-D solid modeling program. Each component has a computer user interface of an extensive command set. All of these programs make use of a comprehensive B-spline mathematics subroutine library, which can be used for a wide variety of interpolation problems and other geometric calculations. Many other user aids, such as automatic saving of the geometric and finite element data bases and hidden line removal, are available. This interactive finite element modeling capability can produce a complete finite element model, producing an output file of grid and element data

The Atomic to Molecular Transition in Galaxies. I: An Analytic Approximation for Photodissociation Fronts in Finite Clouds

In this series of papers we study the structure of the atomic to molecular transition in the giant atomic-molecular complexes that are the repositories of most molecular gas in galaxies, with the ultimate goal of attaining a better understanding of what determines galaxies' molecular content. Here we derive an approximate analytic solution for the structure of a photodissociation region (PDR) in a cloud of finite size that is bathed in an external dissociating radiation field. Our solution extends previous work, which with few exceptions has been restricted to a one-dimensional treatment of the radiation field. We show that our analytic results compare favorably to exact numerical calculations in the one-dimensional limit. However, our more general geometry provides a more realistic representation than a semi-infinite slab for atomic-molecular complexes exposed to the interstellar radiation field, particularly in environments such as low-metallicity dwarf galaxies where the curvature and finite size of the atomic envelope cannot be neglected. For clouds that are at least 20% molecular we obtain analytic expressions for the molecular fraction in terms of properties of the gas and radiation field that are accurate to tens of percent, while for clouds of lower molecular content we obtain upper limits. As a side benefit, our analysis helps clarify when self-shielding is the dominant process in H_2 formation, and under what circumstances shielding by dust makes a significant contribution.Comment: 19 pages, 11 figures, emulateapj style, accepted to ApJ. Discussion slightly changed from previous version, and some new analytic approximations added. Underlying results unchange

Similarities and contrasts in tectonic and volcanic style and history along the Colorado plateaus-to-basin and range transition zone in Western Arizona: Geologic framework for tertiary extensional tectonics

The overall temporal and spatial relations between middle Tertiary volcanism and tectonism from the Basin and Range province onto the edge of the Colorado Plateaus province suggest that a single magnetic-tectonic episode affected the entire region more or less simultaneously during this period. The episode followed a post-Laramide (late Eocene through Oligocene) period of 25 million years of relative stability. Middle Tertiary volcanism did not migrate gradually eastward in a simple fashion onto the Colorado Plateau. In fact, late Oligocene volcanism appears to be more voluminous near the Aquarius Mountains than throughout the adjacent Basin and Range province westward to the Colorado River. Any model proposed to explain the cause of extension and detachment faulting in the eastern part of the Basin and Range province must consider that the onset of volcanism appears to have been approximately synchronous from the Colorado River region of the Basin and Range across the transition zone and onto the edge of the Colorado Plateaus

Far-Infrared Spectral Energy Distributions and Photometric Redshifts of Dusty Galaxies

We infer the large-scale source parameters of dusty galaxies from their observed spectral energy distributions (SEDs) using the analytic radiative transfer methodology presented in Chakrabarti & McKee (2005). For local ultra-luminous infrared galaxies (ULIRGs), we show that the millimeter to far-infrared (FIR) SEDs can be well fit using the standard dust opacity index of 2 when self-consistent radiative transfer solutions are employed, indicating that the cold dust in local ULIRGs can be described by a single grain model. We develop a method for determining photometric redshifts of ULIRGs and sub-mm galaxies from the millimeter-FIR SED; the resulting value of $1+z$ is typically accurate to about 10%. As such, it is comparable to the accuracy of near-IR photometric redshifts and provides a complementary means of deriving redshifts from far-IR data, such as that from the upcoming $\it{Herschel Space Observatory}$. Since our analytic radiative transfer solution is developed for homogeneous, spherically symmetric, centrally heated, dusty sources, it is relevant for infrared bright galaxies that are primarily powered by compact sources of luminosity that are embedded in a dusty envelope. We discuss how deviations from spherical symmetry may affect the applicability of our solution, and we contrast our self-consistent analytic solution with standard approximations to demonstrate the main differences.Comment: 37 pages, 14 Figures, 3 Tables, submitted to ApJ. submitted to Ap

A search for X-ray emission from rich clusters, extended halos around clusters, and superclusters

The all-sky data base acquired with the HEAO A-2 experiment was searched for X-ray emission on a variety of metagalactic size scales which were either predicted or previously detected. Results in the 0.2-60 keV energy range are presented. The optically richest clusters, including those from which a microwave decrement were observed, appear to be relatively underluminous in X-rays. Observations of Abell 576 show its luminosity to be less than earlier estimates, and moreover less than the luminosity predicted from its microwave decrement, unless the intracluster gas is a factor of approximately 10 hotter than in typical clusters. Near SC0627 there are two X-ray sources, and the identification of the dominant source with SCO627 is probably incorrect. New spectral observations of Abell 401 and 2147, possible superclusters, reveal that they have typical cluster spectra with iron line emission

Radiation-Hydrodynamic Simulations of Collapse and Fragmentation in Massive Protostellar Cores

We simulate the early stages of the evolution of turbulent, virialized, high-mass protostellar cores, with primary attention to how cores fragment, and whether they form a small or large number of protostars. Our simulations use the Orion adaptive mesh refinement code to follow the collapse from ~0.1 pc scales to ~10 AU scales, for durations that cover the main fragmentation phase, using three-dimensional gravito-radiation hydrodynamics. We find that for a wide range of initial conditions radiation feedback from accreting protostars inhibits the formation of fragments, so that the vast majority of the collapsed mass accretes onto one or a few objects. Most of the fragmentation that does occur takes place in massive, self-shielding disks. These are driven to gravitational instability by rapid accretion, producing rapid mass and angular momentum transport that allows most of the gas to accrete onto the central star rather than forming fragments. In contrast, a control run using the same initial conditions but an isothermal equation of state produces much more fragmentation, both in and out of the disk. We conclude that massive cores with observed properties are not likely to fragment into many stars, so that, at least at high masses, the core mass function probably determines the stellar initial mass function. Our results also demonstrate that simulations of massive star forming regions that do not include radiative transfer, and instead rely on a barotropic equation of state or optically thin heating and cooling curves, are likely to produce misleading results.Comment: 23 pages, 18 figures, emulateapj format. Accepted to ApJ. This version has minor typo fixes and small additions, no significant changes. Resolution of images severely degraded to fit within size limit. Download the full paper from http://www.astro.princeton.edu/~krumholz/recent.htm

Measurement of electron screening in muonic lead

Energies of the transitions between high-lying (n≥6) states of muonic lead were accurately determined. The results are interpreted as a ∼2% test of the electron screening. The agreement between experiment and theory is good if it is assumed that the refilling of the electron K shell is fast. The present results furthermore severely restrict possible ionization of the electron L shell

The interface between silicon and a high-k oxide

The ability to follow Moore's Law has been the basis of the tremendous success of the semiconductor industry in the past decades. To date, the greatest challenge for device scaling is the required replacement of silicon dioxide-based gate oxides by high-k oxides in transistors. Around 2010 high-k oxides are required to have an atomically defined interface with silicon without any interfacial SiO2 layer. The first clean interface between silicon and a high-K oxide has been demonstrated by McKee et al. Nevertheless, the interfacial structure is still under debate. Here we report on first-principles calculations of the formation of the interface between silicon and SrTiO3 and its atomic structure. Based on insights into how the chemical environment affects the interface, a way to engineer seemingly intangible electrical properties to meet technological requirements is outlined. The interface structure and its chemistry provide guidance for the selection process of other high-k gate oxides and for controlling their growth. Our study also shows that atomic control of the interfacial structure can dramatically improve the electronic properties of the interface. The interface presented here serves as a model for a variety of other interfaces between high-k oxides and silicon.Comment: 10 pages, 2 figures (one color

The HEAO A-2 survey of Abell clusters and the X-ray luminosity function

The HEAO A-2 all sky data base was surveyed for 2-10KeV X-rau emission from the 225 Abell clusters og galaxies listed in Abell's (1958) catalog which are of distance class four or less, and are within the fraction of the sky surveyed completely by Abell. Thirty-two identifications of clusters with X-ray sources were made, for which 2-10 KeV fluxes and 90% error boxes are presented. Twelve of these identification are new. The X-ray luminosity function was derived for this statistically complete sample and the best exponential fit was found to be f(L) = 20.2 x 10 to the minus 8 power exp (-l(44)/1.9) per cu Mpc 2-10KeV. The relationship between X-ray luminosity and richness was examined and a correlation was found for richness classes 0, 1, and 2. The relationship of X-rau luminosity, Bautz-Morgan type, and Rood-Sastry type was examined. It was found that BM type I's and RS type cD and B have the greatest average luminosity. The contribution of clusters to the X-ray background was calculated from the luminosity function and was found to be 5%, and with 90% certainty, less than 8% in the 2-10 KeV band pass

The beamformer and correlator for the Large European Array for Pulsars

The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tied-array beam from the raw voltages from each of telescopes. We explain the concepts and techniques involved in the process of adding the raw voltages coherently. We further present the software processing pipeline that is specifically designed to deal with data from widely spaced, inhomogeneous radio telescopes and describe the steps involved in preparing, correlating and creating the tied-array beam. This includes polarization calibration, bandpass correction, frequency dependent phase correction, interference mitigation and pulsar gating. A link is provided where the software can be obtained.Comment: 10 pages, 6 figures, accepted for publication in Astronomy and Computin