Publications

In-situ filling through hydrolysis and condensation of silicon alkoxides dissolved into polymers has been utilized to generate nanocomposites in which the filler phase can be intimately associated with the polymer on relatively small length scales. One problem of the method has been achieving useful fill volumes without bulk phase separation of the growing inorganic component from the polymer. In this paper, we describe the preparation of a new class of nanocomposite materials in which the inorganic filler phase is pre-assembled before copolymerization with an organic monomer. Maleimide monomers, prepared from alkoxysilylpropyl amines and maleic anhydride, were protected against side reactions by forming the oxonorbornene Diels-Alder adduct with furan. The monomers were then reacted under sol-gel conditions to form oligomers or polymers making up the filler phase. The material was activated by thermal deprotection of the maleimide and reacted with organic monomers or polymers to form the filled nanocomposite.

Progress towards the development of such algorithms as been reported for waveguide analysis'-3and vertical-cavity laser simulation. In all these cases, the higher accuracy order was obtained for a single spatial dimension. More recently, this concept was extended to differencing of the Helmholtz Equation on a 2-D grid, with uniform regions treated to 4th order and dielectric interfaces to 3'd order5. No attempt was made to treat corners properly. In this talk I will describe the extension of this concept to allow differencing of the Helmholtz Equation on a 2-D grid to 6* order in uniform regions and 5* order at dielectric interfaces. In addition, the first known derivation of a finite difference equation for a dielectric comer that allows correct satisfaction of all boundary conditions will be presented. This equation is only accurate to first order, but as will be shown, results in simulations that are third-order-accurate. In contrast to a previous approach3 that utilized a generalized Douglas scheme to increase the accuracy order of the difference second derivative, the present method invokes the Helmholtz Equation itself to convert derivatives of high order in a single direction into mixed

The computing power available to scientists and engineers has increased dramatically in the past decade, due in part to progress in making massively parallel computing practical and available. The expectation for these machines has been great. The reality is that progress has been slower than expected. Nevertheless, massively parallel computing is beginning to realize its potential for enabling significant breakthroughs in science and engineering. This paper provides a perspective on the state of the field, colored by the authors' experiences using large-scale parallel machines at Sandia National Laboratories. We address trends in hardware, system software and algorithms, and we also offer our view of the forces shaping the parallel computing industry.

To physically investigate permeability upscaling, over 13,000 permeability values were measured with four different sample supports (i.e., sample volumes) on a block of Berea Sandstone. At each sample support, spatially exhaustive permeability datasets were measured, subject to consistent flow geometry and boundary conditions, with a specially adapted minipermeameter test system. Here, we present and analyze a subset of the data consisting of 2304 permeability values collected from a single block face oriented normal to stratification. Results reveal a number of distinct and consistent trends (i.e., upscaling) relating changes in key summary statistics to an increasing sample support. Examples include the sample mean and semivariogram range that increase with increasing sample support and the sample variance that decreases. To help interpret the measured mean upscaling, we compared it to theoretical models that are only available for somewhat different flow geometries. The comparison suggests that the nonuniform flow imposed by the minipermeameter coupled with permeability anisotropy at the scale of the local support (i.e., smallest sample support for which data is available) are the primary controls on the measured upscaling. This work demonstrates, experimentally, that it is not always appropriate to treat the local-support permeability as an intrinsic feature of the porous medium, that is, independent of its conditions of measurement.

Isotopic labeling experiments have revealed correlations between hydrogen reactions, Ga desorption, and ammonia decomposition in GaN CVD. Low energy electron diffraction (LEED) and temperature programmed desorption (TPD) were used to demonstrate that hydrogen atoms are available on the surface for reaction after exposing GaN(001) to deuterium at elevated temperatures. Hydrogen reactions also lowered the temperature for Ga desorption significantly. Ammonia did not decompose on the surface before hydrogen exposure. However, after hydrogen reactions altered the surface, N15H3 did undergo both reversible and irreversible decomposition. This also resulted in the desorption of N2 of mixed isotopes below the onset of GaN sublimation. This suggests that the driving force of the high nitrogen-nitrogen bond strength (226 kcal/mol) can lead to the removal of nitrogen from the substrate when the surface is nitrogen rich. Overall, these findings indicate that hydrogen can influence GaN CVD significantly, being a common factor in the reactivity of the surface, the desorption of Ga, and the decomposition of ammonia.

Advanced rechargeable lithium-ion batteries are presently being developed and commercialized worldwide for use in consumer electronics, military and space applications. The motivation behind these efforts involves, among other things, a favorable combination of energy and power density. For some of the applications the power sources may need to perform at a reasonable rate at subambient temperatures. Given the nature of the lithium-ion cell chemistry the low temperature performance of the cells may not be very good. At Sandia National Laboratories, we have used different electrochemical techniques such as impedance and charge/discharge at ambient and subambient temperatures to probe the various electrochemical processes that are occurring in Li-ion cells. The purpose of this study is to identify the component that reduces the cell performance at subambient temperatures. We carried out 3-electrode impedance measurements on the cells which allowed us to measure the anode and cathode impedances separately. Our impedance data suggests that while the variation in the electrolyte resistance between room temperature and -20°C is negligible, the cathode electrolyte interfacial resistance increases substantially in the same temperature span. We believe that the slow interfacial charge transfer kinetics at the cathode electrolyte may be responsible for the increase in cell impedance and poor cell performance. © Copyright 1999 Society of Automotive Engineers, Inc.

This paper addresses key issues for the cost-effective use of COTS (Commercially available Off The Shelf) microelectronics in radiation environments that enable circuit or system designers to manage risks and ensure mission success. We review several factors and tradeoffs affecting the successful application of COTS parts including (1) hardness assurance and qualification issues, (2) system hardening techniques, and (3) life-cycle costs. The paper also describes several experimental studies that address trends in total-dose, transient, and single-event radiation hardness as COTS technology scales to smaller feature sizes. As an example, the level at which dose-rate upset occurs in Samsung SRAMs increases from 1.4x10s rad(Si)/s for a 256K SRAM to 7.7xl09rad(Si)/s for a 4M SRAM, indicating unintentional hardening improvements in the design or process of a commercial technology. Additional experiments were performed to quantify variations in radiation hardness for COTS parts. In one study, only small (10-15%) variations were found in the dose-rate upset and latchup thresholds for Samsung 4M SRAMs from three different date codes. In another study, irradiations of 4M SRAMs from Samsung, Hitachi, and Toshiba indicate large differences in total-dose radiation hardness. The paper attempts to carefully define terms and clear up misunderstandings about the definitions of "COTS'1 and "radiation-hardened (RH)" technology. © 1999 IEEE.

The rutile TiO2(110) (1 × 1) surface is considered the prototypical 'well-defined' system in the surface science of metal oxides. Its popularity results partly from two experimental advantages: (i) bulk-reduced single crystals do not exhibit charging, and (ii) stoichiometric surfaces, as judged by electron spectroscopies, can be prepared reproducibly by sputtering and annealing in oxygen. We present results that show that this commonly applied preparation procedure may result in a surface structure that is by far more complex than generally anticipated. Flat, (1 × 1)-terminated surfaces are obtained by sputtering and annealing in ultrahigh vacuum. When re-annealed in oxygen at moderate temperatures (470-660 K), irregular networks of partially connected, pseudohexagonal rosettes (6.5 × 6 Å wide), one-unit cell wide strands, and small (≈tens of Å) (1 × 1) islands appear. This new surface phase is formed through reaction of oxygen gas with interstitial Ti from the reduced bulk. Because it consists of an incomplete, kinetically limited (1 × 1) layer, this phenomenon has been termed 'restructuring'. We report a combined experimental and theoretical study that systematically explores this restructuring process. The influence of several parameters (annealing time, temperature, pressure, sample history, gas) on the surface morphology is investigated using STM. The surface coverage of the added phase as well as the kinetics of the restructuring process are quantified by LEIS and SSIMS measurements in combination with annealing in 18O-enriched gas. Atomic models of the essential structural elements are presented and are shown to be stable with first-principles density functional calculations. The effect of oxygen-induced restructuring on surface chemistry and its importance for TiO2 and other bulk-reduced oxide materials is briefly discussed.

We describe the microfabrication of an extremely compact optical system as a key element in an integrated capillary channel electrochromatograph with fluorescence detection. The optical system consists of a vertical cavity surface-emitting laser (VCSEL), two high performance microlenses, and a commercial photodetector. The microlenses are multilevel diffractive optics patterned by electron beam lithography and etched by reactive ion etching in fused silica. The design uses substrate-mode propagation within the fused silica substrate. Two generations of optical subsystems are described. The first generation design has a 6 mm optical length and is integrated directly onto the capillary channel-containing substrate. The second generation design separates the optical system onto its own substrate module and the optical path length is further compressed to 3.5 mm. The first generation design has been tested using direct fluorescence detection with a 750 nm VCSEL pumping a 10-4 M solution of CY-7 dye. The observed signal-to-noise ratio of better than 100:1 demonstrates that the background signal from scattered pump light is low despite the compact size of the optical system and is adequate for system sensitivity requirements. © 1999 American Vacuum Society.

A computational procedure for extracting substructure-by-substructure flexibility properties from global modal parameters is presented. The present procedure consists of two key features: an element-based direct flexibility method, which uniquely determines the global flexibility without resorting to case-dependent redundancy selections, and the projection of kinematically inadmissible modes that are contained in the iterated substructural matrices. The direct flexibility method is used as the basis of an inverse problem, whose goal is to determine substructural flexibilities given the global flexibility, geometrically determined substructural rigid-body modes, and the local-to-global assembly operators. The resulting procedure, given accurate global flexibility, extracts the exact element-by-element substructural flexibilities for determinate structures. For indeterminate structures, the accuracy depends on the iteration tolerance limits. The procedure is illustrated using both simple and complex numerical examples and appears to be effective for structural applications such as damage localization and finite element model reconciliation.

Evidence for capillary waves at a liquid-vapor interface are presented from extensive molecular dynamics simulations of a system containing up to 1.24 million Lennard-Jones particles. Careful measurements show that the total interfacial width depends logarithmically on L∥, the length of the simulation cell parallel to the interface, as predicted theoretically. The strength of the divergence of the interfacial width on L∥ depends inversely on the surface tension γ. This allows us to measure γ two ways since γ can also be obtained from the difference in the pressure parallel and perpendicular to the interface. These two independent measures of γ agree provided that the interfacial order parameter profile is fit to an error function and not a hyperbolic tangent, as often assumed. We explore why these two common fitting functions give different results for γ.

Previous photopumping studies of coupled vertical cavity laser structures have demonstrated three mode coupling (two photonic and one excitonic), dual wavelength emission, and short pulse generation. This paper reports on electrically injected coupled resonator vertical-cavity laser diodes, including two novel modulation approaches.

Consider the 2-matching problem defined on the complete graph, with edge costs which satisfy the triangle inequality. We prove that the value of a minimum cost 2-matching is bounded above by 4/3 times the value of its linear programming relaxation, the fractional 2-matching problem. This lends credibility to a long-standing conjecture that the optimal value for the traveling salesman problem is bounded above by 4/3 times the value of its linear programming relaxation, the subtour elimination problem. © Springer-Verlag 1999.

Diamond films were grown by chemical deposition of hydrocarbon species in a vapor composed predominantly of hydrogen. The rate constants of the surface reactions and the concentrations of the gas-phase species were used as input to a variable time step Monte Carlo algorithm, which simulates the evolution of the diamond growth surface by tracking the occupancies of surface sites. The results of the combined tight binding and density functional theory quantum mechanical calculations are presented, suggesting that the etching of isolated, monomolecular moieties occurred at an appreciable rate, while etching from larger carbon islands was nor favorable.

We study the low-temperature in-plane magnetoresistance of tunnel-coupled quasi-one-dimensional quantum wires. The wires are defined by two pairs of mutually aligned split gates on opposite sides of a ⩽1-μm-thick (Formula presented) double-quantum-well heterostructure, allowing independent control of the width of each quantum well. In the ballistic regime, when both wires are defined and the field is perpendicular to the current, a large resistance peak at ∼6 T is observed with a strong gate voltage dependence. The data are consistent with a counting model whereby the number of subbands crossing the Fermi level changes with field due to the formation of an anticrossing in each pair of one-dimensional subbands. © 1999 The American Physical Society.

We describe a very-low power consumption circuit for processing the pulses from a semiconductor radiation detector. The circuit was designed for use with a cadmium zinc telluride (CZT) detector for unattended monitoring of stored nuclear materials. The device is intended to be battery powered and operate at low duty-cycles over a long period of time. This system will provide adequate performance for medium resolution gamma-ray pulse-height spectroscopy applications. The circuit incorporates the functions of a charge sensitive preamplifier, shaping amplifier, and peak sample and hold circuit. An application specific integrated circuit (ASIC) version of the design has been designed, built and tested. With the exception of the input field effect transistor (FET), the circuit is constructed using bipolar components. In this paper the design philosophy and measured performance characteristics of the circuit are described.

We report the stability of TlBa2CaCu2O7 and Tl2Ba2CaCu2O8 on LaAlO3(100) epitaxial thin films, under a variety of conditions. All films are stable in acetone and methanol and with repeated thermal cycling to cryogenic temperatures. Moisture, especially vapor, degrades film quality rapidly. These materials are stable to high temperatures in either N2 or O2 ambients. While total degradation, resulting from Tl depletion, occurs at the same temperatures for both phases, 600 °C in N2 and 700 °C in O2, the onset of degradation occurs at somewhat lower temperatures for TlBa2CaCu2O7 than for Tl2Ba2CaCu2O8.

Tailoring of porous materials involves not only chemical synthetic techniques for tailoring microscopic properties such as pore size, pore shape, pore connectivity, and pore surface reactivity, but also materials processing techniques for tailoring the meso- and the macroscopic properties of bulk materials in the form of fibers, thin films and monoliths. These issues are addressed in the context of five specific classes of porous materials: oxide molecular sieves, porous coordination solids, porous carbons, sol-gel derived oxides, and porous heteropolyanion salts. Reviews of these specific areas are preceded by a presentation of background material and review of current theoretical approaches to adsorption phenomena. A concluding section outlines current research needs and opportunities.

A new multidimensional high lateral resolution ion beam analysis technique, ion-electron emission microscopy (IEEM) is described. Using MeV energy ions, IEEM is shown to be capable of ion beam induced charge collection (IBICC) measurements in semiconductors. IEEM should also be capable of microscopically and multidimensionally mapping the surface and bulk composition of solids. As such, IEEM has nearly identical capabilities as traditional nuclear microprobe analysis, with the advantage that the ion beam does not have to be focused. The technique is based on determining the position where an individual ion enters the surface of the sample by projection secondary electron emission microscopy. The x-y origination point of a secondary electron, and hence the impact coordinates of the corresponding incident ion, is recorded with a position sensitive detector connected to a standard photoemission electron microscope (PEEM). These signals are then used to establish coincidence with IBICC, atomic, or nuclear reaction induced ion beam analysis signals simultaneously caused by the incident ion. © 1999 Elsevier Science B.V. All rights reserved.

A series of controlled impact experiments has been performed to determine the shock loading and release behavior of two types of concrete, differentiated by aggregate size, but with average densities varying by less than 2 percent. Hugoniot stress and subsequent release data was collected over a range of approximately 3 to 25 GPa using a plate reverberation technique in combination with velocity interferometry. The results of the current data are compared to those obtained in previous studies on concrete with a different aggregate size but similar density. Results indicate that the average loading and release behavior are comparable for the three types of concrete discussed in this paper. Residual strain is also indicated from these measurements. © 1999 Elsevier Science Ltd. All rights reserved.

Project 'Rolling Thunder' is a dish/Stirling demonstration project at Ft. Huachuca in southeastern Arizona. As part of the project, Sandia decided to retrofit a SOLO 161 Stirling engine on the CPG-460 at Ft. Huachuca. Although the SOLO 161 PCU has operated nearly flawlessly and the CPG-460 has been, for the most part, a solid and reliable component, integration of the SOLO PCU with the CPG-460 into a functional system required significant attention.

Resistive bolometry is an accurate, robust, spectrally broadband technique for measuring absolute x-ray fluence and flux. Bolometry is an independent technique for x-ray measurements that is based on a different set of physical properties than other diagnostics such as x-ray diodes, photoconducting detectors, and P-I-N diodes. Bolometers use the temperature-driven change in element resistivity to determine the total deposited energy. The calibration of such a device is based on fundamental material properties and its physical dimensions. We describe the use of nickel and gold bolometers to measure x rays generated by high-power z pinches on Sandia's Saturn and Z accelerators. The Sandia bolometer design described herein has a pulse response of ∼1 ns. We describe in detail the fabrication, fielding, and data analysis issues leading to highly accurate x-ray measurements. The fundamental accuracy of resistive bolometry will be discussed. © 1999 American Institute of Physics.

We have developed a comprehensive diagnostic package for observing z-pinch radiation along the pinch axis on the Z accelerator. The instrumentation, fielded on the axial package, are x-ray diagnostics requiring direct lines of sight to the target. The diagnostics require vacuum access to the center of the accelerator. The environment is a hostile one, where we must deal with an intense, energetic photon flux (>100 keV) debris (e.g., bullets or shrapnel), and mechanical shock in order for the diagnostics to survive. In addition, practical constraints require the package be refurbished and utilized on a once a day shot schedule. In spite of this harsh environment, we have successfully fielded the diagnostic package with a high survivability of the data and the instruments. In this article, we describe the environment and issues related to the reentrant diagnostic package's implementation and maintenance. © 1999 American Institute of Physics.

Lattice defects are introduced into the structure to suppress the motion of magnetic vortices and enhance the critical current density in high temperature superconductors. Point defects are not very effective pinning sites for the cuprate superconductors; however, extended defects, such as linear tracks, have been shown to be strong pinning sites. We study the superconducting cuprate Tl-2212 (the numbers designate Tl-Ba-Ca-Cu stoichiometry). Large enhancements of vortex pinning potential were observed in Tl-2212 after high-intermediate energy heavy-ion irradiations where non-continuous extended defects were induced at dE/dx of 9 to 15.2 keV/nm (60 MeV Au, 60 MeV Cu, and 30 MeV Au) and continuous linear defects were induced at 19.5 keV/nm (88 MeV Au). Our research addresses the question of pinning in highly anisotropic materials like Tl-2212 where the vortices are `pancakes' rather than `rods' and suitable defect structures may be discontinuous extended damage domains. The defect microstructure and the effectiveness of the pinning potential in Tl-2212 after irradiation by intermediate energy Au at lower dE/dx of 5-15 keV/nm, where recoils are more significant, is studied using high resolution transmission electron microscopy digital imaging and a SQUID magnetometer. The nature of the ion irradiation damage at these intermediate dE/dx will be correlated to the average vortex pinning potential and the TRIMRC calculations for recoils.

Sol-gel processing of materials is plagued by shrinkage during polymerization of the alkoxide monomers and processing (aging and drying) of the resulting gels. We have developed a new class of hybrid organic-inorganic materials based on the solventless ring-opening polymerization (ROP) of monomers bearing the 2,2,5,5-tetramethyl-2,5-disilaoxacyclopentyl group, which permits us to drastically reduce shrinkage in sol-gel processed materials. Because the monomers are polymerized through a chain growth mechanism catalyzed by base rather than the step growth mechanism normally used in sol-gel systems, hydrolysis and condensation products are entirely eliminated. Furthermore, since water is not required for hydrolysis, an alcohol solvent is not necessary. Monomers with two disilaoxacyclopentyl groups, separated by a rigid phenylene group or a more flexible alkylene group, were prepared through disilylation of the corresponding diacetylenes, followed by ring closure and hydrogenation. Anionic polymerization of these materials, either neat or with 2,2,5,5-tetramethyl-2,5-disila-1-oxacyclopentane as a copolymer, affords thermally stable transparent gels with no visible shrinkage. These materials provide an easy route to the introduction of sol-gel type materials in encapsulation of microelectronics, which we have successfully demonstrated.

Thermally stimulated current (TSC) techniques provide information about oxide-trap charge densities and energy distributions in MOS (metal-oxide-semiconductor) capacitors exposed to ionizing radiation or high-field stress that is difficult or impossible to obtain via standard capacitance-voltage or current-voltage techniques. The precision and reproducibility of measurements through repeated irradiation/TSC cycles on a single capacitor is demonstrated with a radiation-hardened oxide, and small sample-to-sample variations are observed. A small increase in E′δ center density may occur in some non-radiation-hardened oxides during repeated irradiation/TSC measurement cycles. The importance of choosing an appropriate bias to obtain accurate measurements of trapped charge densities and energy distributions is emphasized. A 10 nm deposited oxide with no subsequent annealing above 400 °C shows a different trapped-hole energy distribution than thermally grown oxides, but a similar distribution to thermal oxides is found for deposited oxides annealed at higher temperatures. Charge neutralization during switched-bias irradiation is found to occur both because of hole-electron annihilation and increased electron trapping in the near-interfacial SiO2. Limitations in applying TSC to oxides thinner than approximately 5 nm are discussed.

The grown-in tensile strain, due to a lattice mismatch between AlGaN and GaN, is responsible for the observed cracking that seriously limits the feasibility of nitride-based ultraviolet (UV) emitters. We report in-situ monitoring of strain/stress during MOCVD of AlGaN based on a wafer-curvature measurement technique. The strain/stress measurement confirms the presence of tensile strain during growth of AlGaN pseudomorphically on a thick GaN layer. Further growth leads to the onset of stress relief through crack generation. We find that the growth of AlGaN directly on low-temperature (LT) GaN or AlN buffer layers results in a reduced and possibly controllable strain.

A comparison of KOH, NaOH and AZ400K solutions for UV photo-assisted etching of undoped and n+ GaN is discussed. The etching is diffusion-limited (Ea < 6kCal·mol-1) under all conditions and is significantly faster with bias applied to the sample during light exposure. No etching of InN was observed, due to the very high n-type background doping (> 1020cm-3) in the material.

Stiff polyelectrolytes are found to spontaneously form oriented bundles. Conditions under which bundling occurs are found. Molecular dynamics simulations show that divalent counterions are necessary, and the chains must be sufficiently long and stiff. No aggregation occurs for monovalent counterions. For flexible or short chains, aggregation occurs but bundle formation does not. Because of dynamical constraints, the systems tend to order into a network of connected bundles, not a single bundle. © 1999 The American Physical Society.

A critical issue in the long-term reliability of solder connections used in electronic packages is joint failure during thermal cycling. At present, solder is assumed to be a homogeneous single-phase metal in most finite element analyses to predict solder joint fatigue failures. However, in the last decade, several metallurgical studies have shown that solder microstructure may have a role in early solder joint failures. Investigators have observed that solder microstructure coarsens in local bands during aging and during thermal cycle fatigue. In a failed solder joint, the fatigue cracks are found in these bands of coarse grains. It is speculated that the grain coarsening increases local strains within the microstructure, thereby increasing the likelihood for a crack to initiate. The objective of this study is to model and simulate the effect of grain coarsening on local stresses and strains. During solidification of eutectic Pb/Sn solder, two types of microstructures form: lamellar and equiaxed. In this study, the author has developed a computer code to generate both types of microstructures of varying grain coarseness. This code is incorporated into the finite element code that analyzes the local stresses and strains within the computer-generated microstructure. The FE code, specifically developed for this study, uses an algorithm involving the sparse matrix and iterative solver. This code on a typical single-processor machine will allow the analyst to use over 1 million degrees of freedom.

The minimum expended energy for fracture is the free energy required to form two new surfaces. For intergranular fracture, the minimum surface formation energy is complicated by the rough fracture surface, with area greater than the specimen cross-section. We utilize network optimization algorithms (max-flow/min-cut) to determine the minimum surface formation energies and surfaces for intergranular fracture in 3D polycrystals. For equiaxed grains and uniform boundary strength, the minimum energy fracture area is independent of grain size and is 45% larger than the specimen cross-section, and intergranular fracture will occur when surface energy is less than 1.6 times the grain boundary energy. The 3D fracture area is larger than projected from 2D systems. In systems with microcracked boundaries, the fracture surface deviates to preferentially include microcracked boundaries, creating interlocking grain configurations. Two-dimensional percolation of microcracks occurs at about 80% microcracked boundaries.

An integrated microsystem to detect traces of chemical agents (pChemLab™) is being developed at Sandia for counter-terrorism and nonproliferation applications. This microsystem has two modes of operation: liquid and gas phase detection. For the gas phase detection. we are integrating these critical components: a preconcentrator for sample collection, a gas chromatographic (GC) separator, a chemically selective flexural plate wave (FPW) anay mass detector, and a latching valve onto a single chip. By fabricating these components onto a single integrated system (pChemLab™on a chip), the advantages of reduced dead volume, lower power consumption, and smaller physical size can be realized. In this paper, the development of a latching valve will be presented. The key design parameters for this latching valve are: a volumetric flow rate of 1 mL/min, a maximum hold-off pressure of 40 kPa (6 psi), a relatively low power, and a fast response time. These requirements have led to the design of a magnetically actuated latching relay diaphragm valve. Magnetic actuation is chosen because it can achieve sufficient force to effectively seal against back pressure and its power consumption is relatively low. The actuation time is rapid, and valve can latch in either an open or closed state. A corrugated parylene membrane is used to separate the working fluid from internal components of the valve. Corrugations in the parylene ensure that the diaphragm presents minimum resistance to the actuator for a relativley large deflection. Two different designs and their performance of the magnetic actuation have been evaluated. The first uses a linear magnetic drive mechanism, and the second uses a relay mechanism. Preliminary results of the valve performance indicates that the required driving voltage is about 10 volts, the measured flow rate is about 50 mL/min, and it can hold off pressure of about 5 psi (34 kPa). Latest modifications of the design show excellent performance improvements.

In the ℓ-phylogeny problem, one wishes to construct an evolutionary tree for a set of species represented by characters, in which each state of each character induces no more than ℓ connected components. We consider the fixed-topology version of this problem for fixed-topologies of arbitrary degree. This version of the problem is known to be NP-comPlete for ℓ ≥ 3 even for degree-3 trees in which no state labels more than ℓ+ 1 leaves (and therefore there is a trivial ℓ + 1 phytogeny). We give a 2-approximation algorithm for all ℓ ≥ 3 for arbitrary input topologies and we give an optimal approximation algorithm that constructs a 4-phylogeny when a 3-phylogeny exists. Dynamic programming techniques, which are typically used in fixed-topology problems, cannot be applied to ℓ-phylogeny problems. Our 2-approximation algorithm is the first application of linear programming to approximation algorithms for phylogeny problems. We extend our results to a related problem in which characters are polymorphic.

Evaluation of battery and other energy-storage technologies for stationary uses is progressing rapidly toward application-specific testing. This testing uses computer-based data acquisition and control equipment, active electronic loads and power supplies, and customized software, to enable sophisticated test regimes which simulate actual use conditions. These simulated-use tests provide more accurate performance and life evaluations than simple constant resistance or current testing regimes. Several organizations are cooperating to develop simulated-use tests for utility-scale storage systems, especially battery energy-storage systems (BESSs). Some of the tests use stepped constant-power charge and discharge regimes to simulate conditions created by electric utility applications such as frequency regulation (FR) and spinning reserve (SR). Other test profiles under development simulate conditions for the energy-storage component of remote-area power supplies (RAPSs) which include renewable and/or fossil-fuelled generators. Various RAPS applications have unique sets of service conditions that require specialized test profiles. Almost all RAPS tests and many tests that represent other stationary applications need, however, to simulate significant time periods that storage devices operate at low-to-medium states-of-charge without full recharge. Consideration of these and similar issues in simulated-use test regimes is necessary to predict effectively the responses of the various types of batteries in specific stationary applications. This paper describes existing and evolving stationary applications for energy-storage technologies and test regimes which are designed to simulate them. The paper also discusses efforts to develop international testing standards.

SolarPACES (Solar Power and Chemical Energy Systems) is the International Energy Agency's solar thermal working group. To date, research and development activities sponsored by the group have helped reduce the cost of solar thermal systems to one-fifth that of the early pilot plants. This report presents the collective position of the SolarPACES community on solar thermal electricity-generating technology. Topics discussed include the current status of the technology and likely near-term improvements, the needs of target markets, and important technical and financial issues that must be resolved for success in near-term global markets.

Some time ago we presented evidence that, under nonhydrostatic loading, the F{sub R1} {r_arrow} A{sub O} polymorphic phase transformation in unpoled PZT 95/5-2Nb ceramic began when the maximum compressive stress equaled the hydrostatic pressure at which the transformation otherwise took place. More recently, we showed that this simple stress criterion did not apply to nonhydrostatically compressed, poled ceramic. However, unpoled ceramic is isotropic, whereas poled ceramic has a preferred crystallographic orientation and is mechanically anisotropic. If we further assume that the transformation depends not only on the magnitude of the compressive stress, but also its orientation relative to some feature(s) of PZT 95/5-2Nb's crystallography, then these disparate results can be qualitatively resolved. In this report, we first summarize the existing results for unpoled and poled ceramic. Using our orientation-dependent hypothesis and these results, we derive simple arithmetic expressions that accurately describe our previously-observed effects of nonhydrostatic stress on the transformation of unpoled ceramic. We then go on to test new predictions based on the orientation-dependent model. It has long been known that the transformation can be triggered in uniaxial compression: the model specifically requires a steadily increasing axial stress to drive the transformation of a randomly-oriented polycrystal to completion. We show that when the stress is held constant during uniaxial compression experiments, the transformation stops, supporting our hypothesis. We close with a discussion of implications of our model, and ways to test it using poled ceramic.

A computer code has been developed to solve the linear Boltzmann transport equation on an unstructured mesh of triangles, from a Pro/E model. An arbitriwy arrangement of distinct material regions is allowed. Energy dependence is handled by solving over an arbitrary number of discrete energy groups. Angular de- pendence is treated by Legendre-polynomial expansion of the particle cross sections and a discrete ordinates treatment of the particle fluence. The resulting linear system is solved in parallel with a preconditioned conjugate-gradients method. The solution method is unique, in that the space-angle dependence is solved si- multaneously, eliminating the need for the usual inner iterations. Electron cross sections are obtained from a Goudsrnit-Saunderson modifed version of the CEPXS code. A one-dimensional version of the code has also been develop@ for testing and development purposes.

Small, reliable chemical sensors are needed for a wide range of applications, such as weapon state-of-health monitoring, nonproliferation activities, and manufacturing emission monitoring. Significant improvements in present surface acoustic wave sensors could be achieved by developing a flexural plate-wave (FPW) architecture, in which acoustic waves are excited in a thin sensor membrane. Further enhancement of device performance could be realized by integrating a piezoelectric thin film on top of the membrane. These new FPW-piezoelectric thin film devices would improve sensitivity, reduce size, enhance ruggedness and reduce the operating frequency so that the FPW devices would be compatible with standard digital microelectronics. Development of these piezoelectric thin film // FPW devices requires integration of (1) acoustic sensor technology, (2) silicon rnicromachining techniques to fabricate thin membranes, and (3) piezoelectric thin films. Two piezoelectric thin film technologies were emphasized in this study: Pb(Zr,Ti)O{sub 3} (PZT) and AlN. PZT thin films were of sufficient quality such that the first high frequency SAW measurements on PZT thin films were measured during the course of this study. Further, reasonable ferroelectric properties were obtained from PZT films deposited on Si surface micromachined FPW device membranes. Fundamental understanding of the effect of nanodimension interfacial layers on AlN thin film domain configurations and piezoelectric response was developed. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the US Department of Energy under contract DE-AC04-94AL85000.

As part of the Laboratory-Directed Research and Development (LDRD) Program at Sandia National Laboratories, an investigation into the existence of enhanced vapor-phase diffusion (EVD) in porous media has been conducted. A thorough literature review was initially performed across multiple disciplines (soil science and engineering), and based on this review, the existence of EVD was found to be questionable. As a result, modeling and experiments were initiated to investigate the existence of EVD. In this LDRD, the first mechanistic model of EVD was developed which demonstrated the mechanisms responsible for EVD. The first direct measurements of EVD have also been conducted at multiple scales. Measurements have been made at the pore scale, in a two- dimensional network as represented by a fracture aperture, and in a porous medium. Significant enhancement of vapor-phase transport relative to Fickian diffusion was measured in all cases. The modeling and experimental results provide additional mechanisms for EVD beyond those presented by the generally accepted model of Philip and deVries (1957), which required a thermal gradient for EVD to exist. Modeling and experimental results show significant enhancement under isothermal conditions. Application of EVD to vapor transport in the near-surface vadose zone show a significant variation between no enhancement, the model of Philip and deVries, and the present results. Based on this information, the model of Philip and deVries may need to be modified, and additional studies are recommended.

This report documents a collection of papers on a family of uniform strain tetrahedral finite elements and their connection to different element types. Also included in the report are two papers which address the general problem of connecting dissimilar meshes in two and three dimensions. Much of the work presented here was motivated by the development of the tetrahedral element described in the report "A Suitable Low-Order, Eight-Node Tetrahedral Finite Element For Solids," by S. W. Key {ital et al.}, SAND98-0756, March 1998. Two basic issues addressed by the papers are: (1) the performance of alternative tetrahedral elements with uniform strain and enhanced uniform strain formulations, and (2) the proper connection of tetrahedral and other element types when two meshes are "tied" together to represent a single continuous domain.

Optical channeling or refractive guiding processes involving the nonlinear interaction of intense femtosecond optical pulses with matter in the self-focussing regime has created exciting opportunities for next-generation laser plasma-based x-ray sources and directed energy applications. This fundamentally new form of extended paraxial electromagnetic propagation in nonlinear dispersive media such as underdense plasma is attributed to the interplay between normal optical diffraction and intensity-dependent nonlinear focussing and refraction contributions in the dielectric response. Superposition of these mechanisms on the intrinsic index profile acts to confine the propagating energy in a dynamic self-guiding longitudinal waveguide structure which is stable for power transmission and robust compression. The laser-driven channels are hypothesized to support a degree of solitonic transport behavior, simultaneously stable in the space and time domains (group velocity dispersion balances self-phase modulation), and are believed to be self-compensating for diffraction and dispersion over many Rayleigh lengths in contrast with the defining characteristics of conventional diffractive imaging and beamforming. By combining concentrated power deposition with well-ordered spatial localization, this phenomena will also create new possibilities for production and regulation of physical interactions, including electron beams, enhanced material coupling, and self-modulated plasma wakefields, over extended gain distances with unprecedented energy densities. Harmonious combination of short-pulse x-ray production with plasma channeling resulting from a relativistic charge displacement nonlinearity mechanism in the terawatt regime (10{sup 18} W/cm{sup 2}) has been shown to generate high-field conditions conducive to efficient multi-kilovolt x-ray amplification and peak spectral brightness. Channeled optical propagation with intense short-pulse lasers is expected to impact several critical mission areas at Sandia including x-ray backlighting of pinch implosions, nondestructive radiographic imaging of aging weapons components, high-power electromagnetic pulse generation, particle acceleration, and remote sensing.

This effort studied the integration of innovative methods of key management crypto synchronization, and key agility while scaling encryption speed. Viability of these methods for encryption of ATM cell payloads at the SONET OC- 192 data rate (10 Gb/s), and for operation at OC-48 rates (2.5 Gb/s) was shown. An SNL-Developed pipelined DES design was adapted for the encryption of ATM cells. A proof-of-principle prototype circuit board containing 11 Electronically Programmable Logic Devices (each holding the equivalent of 100,000 gates) was designed, built, and used to prototype a high speed encryptor.

This study evaluated multiple, long-term environmental oil-contamination risk scenarios that could result from the potential leakage of UP to 1.5 million barrels of crude oil entombed in the Weeks Island SPR mine following site decommissioning and abandonment, and up to 100 years thereafter. This risk assessment also provides continuity with similar risk evaluations performed earlier and documented in the 1995 DOE Environmental Assessment for Decommissioning the Strategic Petroleum Reserve Weeks Island Facility (EA). This current study was requested by the DOE to help them determine if their previous Finding of No Significant Impact (FONSI), in the EA, is still valid or needs to be rescinded. Based on the calculated environmental risk results (in terms of clean-up and remediation expenses) presented in this risk assessment, including the calculated average likelihoods of oil release and potential oil-leakage volumes, none of the evaluated risk events would appear to satisfy the definition of significant environmental impact in National Environmental Policy Act (NEPA) terminology. The DOE may combine these current results with their earlier evaluations and interpretations in the 1995 EA in order to assess whether the existing FONSI is still accurate, acceptable, and valid. However, from a risk evaluation standpoint, the assessment of impacts appears to be the same whether only 10,000 to 30,000 barrels of crude oil (as considered in the 1995 EA), or up to 1.5 million barrels of oil (as considered herein) are abandoned in the Weeks Island SPR facility.

The theoretical background for the finite element computer program, MPSalsa Version 1.5, is presented in detail. MPSalsa is designed to solve laminar or turbulent low Mach number, two- or three-dimensional incompressible and variable density reacting fluid flows on massively parallel computers, using a Petrov-Galerkin finite element formulation. The code has the capability to solve coupled fluid flow (with auxiliary turbulence equations), heat transport, multicomponent species transport, and finite-rate chemical reactions, and to solve coupled multiple Poisson or advection-diffusion-reaction equations. The program employs the CHEMKIN library to provide a rigorous treatment of multicomponent ideal gas kinetics and transport. Chemical reactions occurring in the gas phase and on surfaces are treated by calls to CHEMKIN and SURFACE CHEMK3N, respectively. The code employs unstructured meshes, using the EXODUS II finite element database suite of programs for its input and output files. MPSalsa solves both transient and steady flows by using fully implicit time integration, an inexact Newton method and iterative solvers based on preconditioned Krylov methods as implemented in the Aztec. solver library.

The U.S. Department of Energy (DOE) has been developing a nuclear waste disposal facility, the Waste Isolation Pilot Plant (WIPP), located approximately 42 km east of Carlsbad, New Mexico. The WIPP is designed to demonstrate the safe disposal of transuranic wastes produced by the defense nuclear-weapons program. Performance assessment analyses (U.S. DOE, 1996) indicate that human intrusion by inadvertent and intermittent drilling for resources provide the only credible mechanisms for significant releases of radionuclides horn the disposal system. These releases may occur by five mechanisms: (1) cuttings, (2) cavings, (3) spallings, (4) direct brine releases, and (5) long-term brine releases. The first four mechanisms could result in immediate release of contaminant to the accessible environment. For the last mechanism, migration pathways through the permeable layers of rock above the Salado are important, and major emphasis is placed on the Culebra Member of the Rustler Formation because this is the most transmissive geologic layer in the disposal system. For reasons of initial quantity, half-life, and specific radioactivity, certain isotopes of Th, U, Am, and Pu would dominate calculated releases from the WIPP. In order to help quanti~ parameters for the calculated releases, radionuclide transport experiments have been carried out using five intact-core columns obtained from the Culebra dolomite member of the Rustler Formation within the Waste Isolation Pilot Plant (WIPP) site in southeastern New Mexico. This report deals primarily with results of mathematical analyses related to the retardation of %J%, 24%, and 24'Am in two of these cores (B-Core - VPX26-11A and C-Core - VPX28-6C). All B-Core transport experiments were done using Culebra-simukmt brine relevant to the core recovery location (the WIPP air-intake shaft - AIS). Most experiments with C-Core were done with AIS brine with some admixture of a brine composition (ERDA-6) that simulated deeper formation brines. No significant changes in transport behavior were observed for changes in brine. Hydraulic characteristics (i.e., apparent porosity and apparent dispersion coefficient) for the cores were obtained via experiments using conservative tracer `Na. Elution experiments carried out over periods of a few days with tracers `*U and %Np indicated that these tracers were weakly retarded as indicated by delayed elution of these species. Elution experiments with tracers `%, 24'Pu, and 24'Ani were performed, but no elution of any of these species was observed in any flow experiment to date, including experiments of up to two years duration. However, B-Core was subjected to tomographic analysis from which a retardation factor can be inferred for%. Moreover, the fact of non- elution for 24*Pu and 24'Am after more than two years brine flow through C-Core can be coupled with the minimum detectable activity for each of these species to compute minimum retardation factors in C-Core. The retardation factors for all three species can then be coupled with the apparent hydraulic characteristics to estimate an apparent minimum solutionhock distribution coefficient, &, for each actinide. The specific radionuclide isotopes used in these experiments were chosen to facilitate analysis. Even though these isotopes are not necessarily the same as those that are most important to WIPP performance, they are isotopes of the same elements, and . their chemical and transport properties are therefore identical to those of isotopes in the WIPP inventory. The retardation factors and & values deduced from experimental results strongly support the contention that sorption in the Culebra provides an effective barrier to release of Th, Pu, and Am during the regulatory period.

This paper presents the Local Area Network Distributed Supervisory Control and Programming Environment (LandScape) commands set that provides a Generic Device Subsystem Application Programmers Interface (API). These commands are implemented using the Common Object Request Broker Architecture (CORBA) specification with Orbix from Iona Technologies.

The purpose of this project was to create a reliable, 3D sensing and visualization system for unattended monitoring. The system provides benefits for several of Sandia's initiatives including nonproliferation, treaty verification, national security and critical infrastructure surety. The robust qualities of the system make it suitable for both interior and exterior monitoring applications. The 3D sensing system combines two existing sensor technologies in a new way to continuously maintain accurate 3D models of both static and dynamic components of monitored areas (e.g., portions of buildings, roads, and secured perimeters in addition to real-time estimates of the shape, location, and motion of humans and moving objects). A key strength of this system is the ability to monitor simultaneous activities on a continuous basis, such as several humans working independently within a controlled workspace, while also detecting unauthorized entry into the workspace. Data from the sensing system is used to identi~ activities or conditions that can signi~ potential surety (safety, security, and reliability) threats. The system could alert a security operator of potential threats or could be used to cue other detection, inspection or warning systems. An interactive, Web-based, 3D visualization capability was also developed using the Virtual Reality Modeling Language (VRML). The intex%ace allows remote, interactive inspection of a monitored area (via the Internet or Satellite Links) using a 3D computer model of the area that is rendered from actual sensor data.

This report provides an overview of the results of the Vital the Nicaraguan Water Resources Management Initiative, Issues process as implemented for a collaborative effort between the Nicaraguan Ministry of Environment and Natural Resources and Sandia National Laboratories. This initiative is being developed to assist in the development of an efficient and sustainable water resources management system for Nicamgua. The Vital Issues process was used to provide information for developing a project that will develop and implement an advanced information system for managing Nicaragua's water resources. Three Vital Issues panel meetings were convened to 1) develop a mission statement and evaluation criteria for identifying and ranking the issues vital to water resources management in Nicaragua 2) define and rank the vital issues; and 3) identify a preliminary list of information needed to address the vital issues. The selection of panelists from the four basic institutional perspectives- government, industiy, academe, and citizens' groups (through nongovernmental organizations (NGOs))-ensured a high level of stakeholder representation on the panels. The already existing need for a water resource management information system has been magnified in the aftemnath of Hurricane Mitch. This information system would be beneficial for an early warning system in emergencies, and the modeling and simulation capabilities of the system would allow for advanced planning. Additionally, the outreach program will provide education to help Nicaraguan improve their water hygiene practices.

Materials compatibility studies of aged, engineered materials and hardware are critical to understanding and predicting component reliability, particularly for systems with extended stockpile life requirements. Nondestructive testing capabilities for component reliability would significantly enhance lifetime predictions. For example, if the detection of crack propagation through a solder joint can be demonstrated, this technique could be used to develop baseline information to statistically determine solder joint lifelengths. This report will investigate high frequency signal response techniques for nondestructively evaluating the electrical behavior of thick film hybrid transmission lines.

Parameters have been identified that can be modeled stochastically using PORFLOW and Latin Hypercube Sampling (LHS). These parameters include hydrologic and transport properties in the vadose and saturated zones, as well as source-term parameters and infiltration rates. A number of resources were used to define the parameter distributions, primarily those provided in the Retrieval Performance Evaluation Report (Jacobs, 1998). A linear rank regression was performed on the vadose-zone hydrologic parameters given in Khaleel and Freeman (1995) to determine if correlations existed between pairs of parameters. No strong correlations were found among the vadose-zone hydrologic parameters, and it was recommended that these parameters be sampled independently until future data or analyses reveal a strong correlation or functional relationship between parameters. Other distributions for source-term parameters, infiltration rates, and saturated-zone parameters that are required to stochastically analyze the performance of the AX Tank Farm using LHS/PORFLOW were adapted from distributions and values reported in Jacobs (1998) and other literature sources. Discussions pertaining to the geologic conceptualization, vadose-zone modeling, and saturated-zone modeling of the AX Tank Farm are also presented.