Publications

Shock-wave experiments on iron preheated to 1,573 K conducted from 14 to 73 GPa, yield new data for sound velocities of the {gamma}- and liquid-phases. Melting was observed in the highest pressure ({approximately} 71 {+-} 2 GPa) experiments at calculated shock temperatures of 2,775 {+-} 160 K. This single crossing of the {gamma}-liquid boundary measured here agrees closely with the {gamma}-iron melting line determined by Boehler [1993], Saxena et al. [1993], and Jephcoat and Besedin [1997]. This {gamma}-iron melting curve is {approximately} 300 C lower than that of Shen et al. [1998b] at 80 GPa.

In this paper we describe a unified, hierarchical computational approach to aging and reliability problems caused by materials changes in the oxide layers of Si-based microelectronic devices. We apply this method to a particular low-dose-rate radiation effects problem

A response surface methodology-based technique is presented for treating discretization error in non-deterministic analysis. The response surface, or metamodel, is estimated from computer experiments which vary both uncertain physical parameters and the fidelity of the computational mesh. The resultant metamodel is then used to propagate the variabilities in the continuous input parameters, while the mesh size is taken to zero, its asymptotic limit. With respect to mesh size, the metamodel is equivalent to Richardson extrapolation, in which solutions on coarser and finer meshes are used to estimate discretization error. The method is demonstrated on a one dimensional prismatic bar, in which uncertainty in the third vibration frequency is estimated by propagating variations in material modulus, density, and bar length. The results demonstrate the efficiency of the method for combining non-deterministic analysis with error estimation to obtain estimates of total simulation uncertainty. The results also show the relative sensitivity of failure estimates to solution bias errors in a reliability analysis, particularly when the physical variability of the system is low.

In this paper we evaluate the design of the hybrid evolutionary algorithms (EAs) that are currently used to perform flexible ligand binding in the Autodock docking software. Hybrid EAs incorporate specialized operators that exploit domain-specific features to accelerate an EA's search. We consider hybrid EAs that use an integrated local search operator to reline individuals within each iteration of the search. We evaluate several factors that impact the efficacy of a hybrid EA, and we propose new hybrid EAs that provide more robust convergence to low-energy docking configurations than the methods currently available in Autodock.

Precision high aspect-ratio micro molds constructed by deep x-ray lithography have been used to batch fabricate accurately shaped bonded rare-earth based permanent magnets with features as small as 5 microns and thicknesses up to 500 microns. Maximum energy products of up to 8 MGOe have been achieved with a 20%/vol. epoxy bonded melt-spun isotropic Nd2Fe14b powder composite. Using individually processed sub- millimeter permanent sections multipole rotors have been assembled. Despite the fact that these permanent magnet structures are small, their magnetic field producing capability remains the same as at any scale. Combining permanent magnet structures with soft magnetic materials and micro-coils makes possible new and more efficient magnetic microdevices.

A test technique has been devised which is suitable for the testing of the bond strength of batch diffusion bonded LIGA or DXRL defined structures. The method uses a torsion tester constructed with the aid of LIGA fabrication and distributed torsion specimens which also make use of the high aspect ratio nature of DXRL based processing. Measurements reveal achieved bond strengths of 130MPa between electroplated nickel with a bond temperature of 450 C at 7 ksi pressure which is a sufficiently low temperature to avoid mechanical strength degradation.

A technique using diffusion bonding (or solid-state welding) has been used to achieve batch fabrication of two- level nickel LIGA structures. Interlayer alignment accuracy of less than 1 micron is achieved using press-fit gauge pins. A mini-scale torsion tester was built to measure the diffusion bond strength of LIGA formed specimens that has shown successful bonding at temperatures of 450"C at 7 ksi pressure with bond strength greater than 100 Mpa. Extensions to this basic process to allow for additional layers and thereby more complex assemblies as well as commensurate packaging are discussed.

Oxidation of layers of high-Al-content III-V materials by water vapor has become the enabling process for high-efficiency vertical cavity surface emitting lasers (VCSELS) and has potential applications for reducing substrate current leakage in GaAs-on-insulator (GOI) MESFETS. Because of the established importance of wet oxidation in optoelectronic devices and its potential applications in electronic devices, it has become increasingly important to understand the mechanism of wet oxidation and how it might be expected to affect both the fabrication and subsequent operation of devices that have been made using this technique. The mechanism of wet oxidation and the consequence of this mechanism for heterostructure design and ultimate device operation are discussed here.

Engineering systems are becoming increasingly complex as state of the art technologies am incorporated into designs. Surety modeling and analysis is an emerging science that permits an engineer to qualitatively and quantitatively predict and assess the completeness and predictability of a design. Surety is a term often used in the Department of Defense (DoD) and Department of Energy (DOE) communities, which refers to the integration of safety, security, reliability and performance aspects of design. Current risk assessment technologies for analyzing complex systems fail to adequately describe the problem, thus making assessment fragmented and non-integrated. To address this problem, we have developed a methodology and extensible software tool set to address model integration and complexity for high consequence systems. The MultiGraph Architecture (MGA) facilitates multi-domain, model-integrated modeling and analyses of complex, high-assurance systems. The MGA modeling environment allows the engineer to customize the modeling environment to match a design paradigm representative of the actual design. Previous modeling tools have a point-defined model space that forms the modeler to work in less than optimal environments. Current approaches for the problem to be bounded and constrained by requirements of the modeling tool and not the actual design problem. In some small cases, this is only maximally adequate MM facilitates the implementation of a surety methodology, which is used to represent high assurance systems with respect to safety and reliability. Formal mathematical models am used to correctly describe design safety and reliability functionality and behavioral fictional and behavioral representations of the design w then analyzed using commercial-off-the-shelf tools.

The authors have used low-energy electron microscopy to investigate the dynamics of the Si(111) 7 x 7 {r_arrow} 1 x 1 phase transition. Because the densities of the two phases differ, the phase transformation is analogous to precipitation in bulk systems: additional material must diffuse to the phase boundaries in order for the transformation to occur. By measuring the size evolution of an ensemble of domains, and comparing the results to simulations, they have identified a new mechanism of precipitate growth. The source of material necessary for the transformation is the random creation of atom/vacancy pairs at the surface. This mechanism contrasts sharply with classical theories of precipitation, in which mass transport kinetics determine the rate of transformation.

Relative concentrations of reactive ions, neutral radicals, resist and substrate etch products have been measured in dielectric etch chemistries using an uncollided beam mass spectrometer / ion extractor from Hiden Analytical. Analysis techniques employed include both electron impact ionization and dissociative ionization of neutral gas, and potential bias extraction of positive ions from the reactor discharge volume. Measurements were made in C{sub 2}F{sub 6} and CHF{sub 3} discharges in an inductively coupled plasma (ICP-GEC) research reactor operating with power densities, pressures, gas compositions and wafer materials typical of those found in etch processing tools. Wafer substrates investigated included blanket silicon wafers and silicon wafers with varying amounts of photo-resist coverage of the surface (20%, 80% and 100%). In C{sub 2}F{sub 6} discharges CF{sub 3}{sup +} was consistently the dominant fluorocarbon ion present, in agreement with published cross sections for dissociative ionization [ 1,2.3,4.5,6]. Smaller concentrations of CF+, CF{sub 2}{sup -}, and C{sub 2}F{sub 5}{sup +}, were also observed, though the dissociative ionization production of C{sub 2}F{sub 5}{sup +} was a factor of five smaller than would be expected from published cross section values. The presence of photo-resist, even in small amounts, was found to produce marked changes in the discharge composition. For example in C{sub 2}F{sub 6} discharges, concentrations of SiF{sub x} etch products relative to concentrations of C{sub x}F{sub y} species were notably diminished and larger concentrations of water vapor were observed when resist was present. In CHF{sub 3} discharges, CF{sub 3}{sup +} and CHF{sub 2}{sup +} were found to be the main species present, along with smaller concentrations of CF{sub 2}{sup +}, CF{sup +}, CHF{sup +}, CH{sup +} and F{sup -}.

We show that the signal and idler beams generated by certain types of unseeded, nanosecond optical parametric oscillators are tilted and angularly dispersed and have anomalously large bandwidths. This effect is demonstrated in both laboratory measurements and a numerical model. We show how the optical cavity design influences the tilts and how they can be eliminated or minimized. We also determine the conditions necessary to injection seed these parametric oscillators.

Galerkin methods are used in separable Hilbert spaces to construct and compute L{sup 2} [0,{pi}] solutions to large classes of differential equations. In this note a Galerkin method is used to construct series solutions of a nonhomogeneous Sturm-Liouville problem defined on [0,{pi}]. The series constructed are shown to converge to a specified du Bois-Reymond function f in L{sup 2} [0,{pi}]. It is then shown that the series solutions can be made to converge uniformly to the specified du Bois-Reymond function when averaged by the Ces{'a}ro-one summability method. Therefore, in the Ces{'a}ro-one sense, every continuous function f on [0,{pi}] is the uniform limit of solutions of nonhomogeneous Sturm-Liouville problems.

This paper describes the microseismic mapping of repeated injections of drill cuttings into two separate formations at a test site near Mounds, OK. Injections were performed in sandstone and shale formations at depths of 830 and 595 m, respectively. Typical injection disposal was simulated using multiple small-volume injections over a three-day period, with long shut-in periods interspersed between the injections. Microseismic monitoring was achieved using a 5-level array of wireline-run, triaxial- accelerometer receivers in a monitor well 76 m from the disposed well. Results of the mapped microseismic locations showed that the disposal domti W= generally aligns with the major horizontal stress with some variations in azimuth and that wide variations in height and length growth occurred with continued injections. These experiments show that the cuttings injection process cm be adequately monitored from a downhole, wireline-run receiver array, thus providing process control and environmental assurance.

LIGA fabricated materials and components exhibit several processing issues affecting their metallurgical and mechanical properties, potentially limiting their usefulness for MEMS applications. For example, LIGA processing by metal electrodeposition is very sensitive to deposition conditions which causes significant processing lot variations of mechanical and metallurgical properties. Furthermore, the process produces a material with a highly textured lenticular rnicrostructural morphology suggesting an anisotropic material response. Understanding and controlling out-of-plane anisotropy is desirable for LIGA components designed for out-of-plane flexures. Previous work by the current authors focused on results from a miniature servo-hydraulic mechanical test frame constructed for characterizing LIGA materials. Those results demonstrated microstructural and mechanical properties dependencies with plating bath current density in LIGA fabricated nickel (LIGA Ni). This presentation builds on that work and fosters a methodology for controlling the properties of LIGA fabricated materials through processing. New results include measurement of mechanical properties of LIGA fabricated copper (LIGA Cu), out-of-plane and localized mechanical property measurements using compression testing and nanoindentation of LIGA Ni and LIGA Cu.

The use of cross-borehole ground penetrating radar (GPR) imaging for determining g the two dimensional (2D) in situ moisture content distribution within the vadose zone is being investigated. The ultimate goal is to use the GPR images as input to a 2D hydrologic inversion scheme for recovering the van Genuchten parameters governing unsaturated ,hydraulic flow. Initial experiments conducted on synthetic data have shown that at least in theory, cross-borehole GPR measurements can provide realistic estimates of the spatial variation in moisture content that are needed for this type of hydrologic inversion scheme. However, the method can not recover exact values of moisture content due to the break down of the empirical expression often employed to convert GPR velocity images to moisture content, and to the smearing nature of the imaging algorithm. To test the applicability of this method in a real world environment cross- borehole GPR measurements were made at a hydrologic/geophysical vadose zone test site in Socorro, New Mexico. Results show that the GPR images compare well with the uncalibrated borehole neutron log data. GPR data acquisition will continue once an infiltration test has started, and the results from these measurements will be employed in a 2D hydrologic inverse scheme.

We are developing a method for analysis of airborne microparticles based on laser ablation of individual molecules in an ion trap mass spectrometer. Airborne particles enter the spectrometer through a differentially-pumped inlet, are detected by light scattered from two CW laser beams, and sampled by a pulsed excimer laser as they pass through the center of the ion trap electrodes. After the laser pulse, the stored ions are separated by conventional ion trap methods. The mass spectra are then analyzed using genetically-trained neural networks (NNs). A number of mass spectra are averaged to obtain training cases which contain a recognizable spectral signature. Averaged spectra for a bacteria and a non-bacteria are shown to the NNs, the response evaluated, and the weights of the connections between neurodes adjusted by a Genetic Algorithm (GA) such that the output from the NN ranges from 0 for non-bacteria to 1 for bacteria. This process is iterated until the population of the GA converges or satisfies predetermined stopping criteria. Using this type of bipolar training we have obtained generalizing NNs able to distinguish five new bacteria from five new non-bacteria, none of which were used in training the NN.

Calculating safety and reliability probabilities with functions of uncertain variables can yield incorrect or misleading results if some precautions are not taken. One important consideration is the application of constrained mathematics for calculating probabilities for functions that contain repeated variables. This paper includes a description of the problem and develops a methodology for obtaining an accurate solution.

The Information-efficient Spectral Imaging Sensor (ISIS) approach to spectral imaging seeks to bridge the gap between tuned multispectral and fixed hyperspectral imaging sensors. By allowing the definition of completely general spectral filter functions, truly optimal measurements can be made for a given task. These optimal measurements significantly improve signal-to-noise ratio (SNR) and speed, minimize data volume and data rate, while preserving classification accuracy. The following paper investigates the application of the ISIS sensing approach in two sample biomedical applications: prostate and colon cancer screening. It is shown that in these applications, two to three optimal measurements are sufficient to capture the majority of classification information for critical sample constituents. In the prostate cancer example, the optimal measurements allow 8% relative improvement in classification accuracy of critical cell constituents over a red, green, blue (RGB) sensor. In the colon cancer example, use of optimal measurements boost the classification accuracy of critical cell constituents by 28% relative to the RGB sensor. In both cases, optimal measurements match the performance achieved by the entire hyperspectral data set. The paper concludes that an ISIS style spectral imager can acquire these optimal spectral images directly, allowing improved classification accuracy over an RGB sensor. Compared to a hyperspectral sensor, the ISIS approach can achieve similar classification accuracy using a significantly lower number of spectral samples, thus minimizing overall sample classification time and cost.

The MELCOR and VICTORIA severe accident analysis codes, which were developed at Sandia National Laboratories for the U. S. Nuclear Regulatory Commission, are designed to estimate fission product releases during nuclear reactor accidents in light water reactors. MELCOR is an integrated plant-assessment code that models the key phenomena in adequate detail for risk-assessment purposes. VICTORIA is a more specialized fission- product code that provides detailed modeling of chemical reactions and aerosol processes under the high-temperature conditions encountered in the reactor coolant system during a severe reactor accident. This paper focuses on recent enhancements and assessments of the two codes in the area of fission product chemistry modeling. Recently, a model for iodine chemistry in aqueous pools in the containment building was incorporated into the MELCOR code. The model calculates dissolution of iodine into the pool and releases of organic and inorganic iodine vapors from the pool into the containment atmosphere. The main purpose of this model is to evaluate the effect of long-term revolatilization of dissolved iodine. Inputs to the model include dose rate in the pool, the amount of chloride-containing polymer, such as Hypalon, and the amount of buffering agents in the containment. Model predictions are compared against the Radioiodine Test Facility (RTF) experiments conduced by Atomic Energy of Canada Limited (AECL), specifically International Standard Problem 41. Improvements to VICTORIA's chemical reactions models were implemented as a result of recommendations from a peer review of VICTORIA that was completed last year. Specifically, an option is now included to model aerosols and deposited fission products as three condensed phases in addition to the original option of a single condensed phase. The three-condensed-phase model results in somewhat higher predicted fission product volatilities than does the single-condensed-phase model. Modeling of U02 thermochemistry was also improved, and results in better prediction of vaporization of uranium from fuel, which can react with released fission products to affect their volatility. This model also improves the prediction of fission product release rates from fuel. Finally, recent comparisons of MELCOR and VICTORIA with International Standard Problem 40 (STORM) data are presented. These comparisons focus on predicted therrnophoretic deposition, which is the dominant deposition mechanism. Sensitivity studies were performed with the codes to examine experimental and modeling uncertainties.

We have performed Monte Carlo calculations to determine the charge accumulation on threading edge dislocations in GaN as a function of the dislocation density and background dopant density. Four possible core structures have been examined, each of which produces defect levels in the gap and may therefore act as electron or hole traps. Our results indicate that charge accumulation, and the resulting electrostatic interactions, can change the relative stabilities of the different core structures. Structures having Ga and N vacancies at the dislocation core are predicted to be stable under nitrogen-rich and gallium-rich growth conditions, respectively. Due to dopant depletion at high dislocation density and the multitude of charge states, the line charge exhibits complex crossover behavior as the dopant and dislocation densities vary.

The linewidth enhancement factor in single quantum-well GRINSCH semiconductor lasers is investigated theoretically and experimentally. For thin wells a small linewidth enhancement factor is obtained which clamps with increasing carrier density, in contrast to the monotonous increase observed for thicker wells. Microscopic many-body calculations reproduce the experimental observations attributing the clamping to a subtle interplay between excitation dependent gain shifts and carrier population distributions.

InGaAsN is a new semiconductor alloy system with the remarkable property that the inclusion of only 2% nitrogen reduces the bandgap by more than 30%. In order to help understand the physical origin of this extreme deviation from the typically observed nearly linear dependence of alloy properties on concentration, we have investigated the pressure dependence of state energies using both experimental and theoretical methods. the excited We report measurements of the low temperature photoluminescence of the material for pressures between ambient and 110 kbar. We also describe a simple, density- functional-theory-based approach to calculating the pressure dependence of low lying excitation energies for low concentration alloys. The theoretically predicted pressure dependence of the bandgap is in excellent agreement with the experimental data. Based on the results of our calculations, we suggest an explanation for the strongly non-linear pressure dependence of the bandgap that, surprisingly, does not involve a nitrogen impurity state.

A low energy electron microscope (LEEM) has been used to investigate the faceting of W(111) as induced by Pt. The atomically rough W(111) surface, when fully covered with a monolayer film of Pt and annealed to temperatures higher than {approximately} 750 K, experiences a significant morphological restructuring: the initially planar surface undergoes a faceting transition and forms three-sided pyramids with {211} faces. The experiments demonstrate the capability of LEEM for imaging both the fully and partially faceted surface. In addition, we have observed the formation of the facets in real time, when Pt is dosed onto the heated surface. We find that the transition from planar surface, to partially faceted surface, and to fully faceted surface proceeds through the nucleation and growth of spatially separated faceted regions.

Leading an audit team goes beyond performance of the duties outlined in any requirement or training course. Anyone can memorize the steps to begin and to complete an audit, but it takes leadership to capitalize on the strengths of each team member and to interact with the auditee. Leadership has been written about and studied for many years. Principles and ideas developed by Covey, Senge, Peters, Blanchard, Hersey, Drucker, Yuki and many, many more but they all come down to some basic issues. There is no magic formula. There are theories and models that when applied work in one situation. Some theories and methodologies work better than others depending on the situation. The presentation today looks at leadership from the perspective of the lead auditor, as he/she has to guide the audit process and deal with many personalities from the audit team to the people being interviewed, Each situation is different, each audit team is different, each audit is unique. The basic principles are applied but it takes understanding leadership to have a successful audit. Applying the Situational Leadership model will enable you to be a good and effective leader and capitalize on the strengths of each team member. It is an invaluable asset to add to your communication toolbox. So let's put our bells on the shelf and put on our learning caps.