We point out that neutrino events observed at Kamiokande andIMB from SN1987A disfavor the neutrino oscillation parameters preferredby the LSND experiment. For Delta m2>0 (the light side), theelectron neutrinos from the neutronization burst would be lost, while thefirst event at Kamiokande is quite likely to be due to an electronneutrino. For Delta m2<0 (the dark side), the average energy of thedominantly bar nu e events is already lower than the theoreticalexpectations, which would get aggravated by a complete conversion frombar nu mu to bar nu e. If taken seriously, the LSND data are disfavoredindependent of the existence of a sterile neutrino. A possible remedy isCPT violation, which allows different mass spectra for neutrinos andanti-neutrinos and hence can accommodate atmospheric, solar and LSND datawithout a sterile neutrino. If this is the case, Mini-BooNE must run inbar nu rather than the planned nu mode to test the LSND signal. Wespeculate on a possible origin of CPT violation.

The energetic and kinetic properties of surfaces play a critical role in defining the microstructural changes that occur during sintering and high-temperature use of ceramics. Characterization of surface diffusion in ceramics is particularly difficult, and significant variations in reported values of surface diffusivities arise even in well-studied systems. Effects of impurities, surface energy anisotropy, and the onset of surface attachment limited kinetics (SALK) are believed to contribute to this variability. An overview of the use of Rayleigh instabilities as a means of characterizing surface diffusivities is presented. The development of models of morphological evolution that account for effects of surface energy anisotropy is reviewed, and the potential interplay between impurities and surface energy anisotropy is addressed. The status of experimental studies of Rayleigh instabilities in sapphire utilizing lithographically introduced pore channels of controlled geometry and crystallography is summarized. Results of model studies indicate that impurities can significantly influence both the spatial and temporal characteristics of Rayleigh instabilities; this is attributed at least in part to impurity effects on the surface energy anisotropy. Related model experiments indicate that the onset of SALK may also contribute significantly to apparent variations in surface diffusion coefficients.

I present a critical account of the final-state interaction (FSI) in two-body B decays from viewpoint of the hadron picture. I emphasize that the phase and the magnitude of decay amplitude are related to each other by a dispersion relation. In a model phase of FSI motivated by experiment, I illustrate how much the magnitude of amplitude can deviate from its factorization value by the FSI.

Following the May 14, 1997 chemical explosion at Hanford's Plutonium Reclamation Facility, the Department of Energy Richland Operations Office and its prime contractor, Fluor Hanford, Inc., completed an extensive assessment to identify and address chemical and radiological safety vulnerabilities at all facilities under the Project Hanford Management Contract. This was a challenging undertaking because of the immense size of the problem, unique technical issues, and competing priorities. This paper focuses on the assessment process, including the criteria and methodology for data collection, evaluation, and risk-based scoring. It does not provide details on the facility-specific results and corrective actions, but discusses the approach taken to address the identified vulnerabilities.

We examine the constraints that satellite-acquired Type Ia and IIP supernova apparent magnitude versus redshift data will place on cosmological model parameters in models with and without a constant or time-variable cosmological constant lambda. High-quality data which could be acquired in the near future will result in tight constraints on these parameters. For example, if all other parameters of a spatially-flat model with a constant lambda are known, the supernova data should constrain the non-relativistic matter density parameter omega to better than 1 (2, 0.5) at 1 sigma with neutral (worst case, best case) assumptions about data quality.

Monitoring of geologic sequestration projects will require the measurement of many different parameters and processes at many different locations at the surface and in the subsurface. The greatest need for technology development is for monitoring of processes in the subsurface in the region between wells. The approach to fitting this need is to build upon decades of experience in use of geophysics in the oil and gas industry. These methods can be optimized for CO2 monitoring, and customized and extended in order to meet the need for cost-effective methods applicable to saline disposal sites, coal bed methane sites, as well as oil and gas reservoir sequestration sites. The strategy for development of cost-effective methods follows a three step iterative process of sensitivity analysis using numerical and experimental techniques, field testing at a range of scale in different formations, and analysis and integration of complimentary types of data.

Reliably processing, imaging, and interpreting seismic data from areas with complicated structures, such as sub-salt, requires a thorough understanding of elastic as well as acoustic wave propagation. Elastic numerical modeling is an essential tool to develop that understanding. While 2-D elastic modeling is in common use, 3-D elastic modeling has been too computationally intensive to be used routinely. Recent advances in computing hardware, including commodity-based hardware, have substantially reduced computing costs. These advances are making 3-D elastic numerical modeling more feasible. A series of example 3-D elastic calculations were performed using a complicated structure, the SEG/EAGE salt structure. The synthetic traces show that the effects of shear wave propagation can be important for imaging and interpretation of images, and also for AVO and other applications that rely on trace amplitudes. Additional calculations are needed to better identify and understand the complex wave propagation effects produced in complicated structures, such as the SEG/EAGE salt structure.

Elastic scattering of low (10-50 eV) kinetic energy electrons from free diatomic molecules is studied using a single-center expansion of the full molecular potential. Dynamic exchange and polarization are included in a local form. The calculated elastic differential scattering cross-sections (DCS) for electron impact on CO and N2 are in good agreement with available experimental data. The importance of using the full molecular potential instead of a two-center potential approach is pointed out. These corrections are small for energies above 50 eV, but they become increasingly important at lower energies. When discussing the angular distributions of elastically-scattered electrons from oriented molecules (like surface adsorbates), we show that these corrections are particularly significant. The results have implications for other electron scattering problems such as those encountered in low-energy photoelectron diffraction from both core and valence levels.

The Galerkin Finite Element Method was used to predict a natural convection flow in an enclosed cavity. The problem considered was a differentially heated, tall (8:1), rectangular cavity with a Rayleigh number of 3.4 x 10{sup 5} and Prandtl number of 0.71. The incompressible Navier-Stokes equations were solved using a Boussinesq approximation for the buoyancy force. The algorithm was developed for efficient use on massively parallel computer systems. Emphasis was on time-accurate simulations. It was found that the average temperature and velocity values can be captured with a relatively coarse grid, while the oscillation amplitude and period appear to be grid sensitive and require a refined computation.

The monitoring of nuclear explosions on a global basis requires accurate event locations. As an example, under the Comprehensive Test Ban Treaty, the size of an on-site inspection search area is 1,000 square kilometers or approximately 17 km accuracy assuming a circular area. This level of accuracy is a significant challenge for small events that are recorded using a sparse regional network. In such cases, the travel-time of seismic energy is strongly affected by crustal and upper mantle heterogeneity and large biases can result. This can lead to large systematic errors in location and, more importantly, to invalid error bounds associated with location estimates. Corrections can be developed and integrated to correct for these biases. These path corrections take the form of both three-dimensional model corrections along with three-dimensional empirically based travel time corrections. LLNL is currently working to integrate a diverse set of three-dimensional velocity model and empirical based travel-time products into one consistent and validated calibration set. To perform this task, we have developed a hybrid approach that uses three-dimensional model corrections for a region and then uses reference events when available to improve the path correction. This Bayesian kriging approach uses the best apriori three-dimensional velocity model …

We carried out a set of experiments on the Omega laser facility at Rochester with Laser MegaJoule (LMJ) like indirect drive irradiation. We studied the irradiation non-uniformity with the foam ball radiography technique and the implosion symmetry with (D{sub 2} + Argon) filled capsules core emission. Cylindrical ''Nova scale 1'' thin wall hohlraums were used. Forty of the Omega beams, arranged in three cones on each side of the hohlraum (5, 5, and lo), were used to create the X-ray drive. Eight additional beams were used on a Ti source to radiograph the foam balls. The shaped laser pulse was about 3 ns duration. The radiation drive was measured on each shot. The images were recorded with a 5 prn resolution Gated X-ray Imager coupled to a CCD camera.

The Yucca Mountain Site Characterization Project is conducting a drift scale heater test, known as the Drift Scale Test (DST), in an alcove of the Exploratory Studies Facility at Yucca Mountain, Nevada. The DST is a large-scale, long-term thermal test designed to investigate coupled thermal-mechanical-hydrological-chemical behavior in a fractured, welded tuff rock mass. The general layout of the DST is shown in Figure 1a, along with the locations of several of the boreholes being used to monitor deformation during the test. Electric heaters are being used to heat a planar region of rock that is approximately 50 m long and 27 m wide for 4 years, followed by 4 years of cooling. Both in-drift and ''wing'' heaters are being used to heat the rock. The heating portion of the DST was started in December, 1997, and the target drift wall temperature of 200 C was reached in summer 2000. A drift-scale distinct element model (DSDE) is being used to analyze the geomechanical response of the rock mass forming the DST. The distinct element method was chosen to permit explicit modeling of fracture deformations. Shear deformations and normal mode opening of fractures are expected to increase fracture permeability and thereby alter …

We report the development of techniques to diagnose plasmas produced by X-ray photoionization of thin foils placed near the Z-pinch on the Sandia Z Machine. The development of 100+ TW X-ray sources enables access to novel plasma regimes, such as the photoionization equilibrium. To diagnose these plasmas one must simultaneously characterize both the foil and the driving pinch. The desired photoionized plasma equilibrium is only reached transiently for a 2-ns window, placing stringent requirements on diagnostic synchronization. We have adapted existing Sandia diagnostics and fielded an additional gated 3-crystal Johann spectrometer with dual lines of sight to meet these requirements. We present sample data from experiments in which 1 cm, 180 eV tungsten pinches photoionized foils composed of 200{angstrom} Fe and 300{angstrom} NaF co-mixed and sandwiched between 1000{angstrom} layers of Lexan (CHO), and discuss the application of this work to benchmarking astrophysical models.

The Rayleigh-Taylor (RT) instability, which occurs when a lower-density fluid accelerates a higher-density layer, is common in nature. At an ablation front a sharp reduction in the growth rate of the instability at short wave-lengths can occur, in marked contrast to the classical case where growth rates are highest at the shortest wavelengths. Theoretical and numerical investigations of the ablative RT instability are numerous and differ considerably on the level of stabilization expected. We present here the results of a series of laser experiments designed to probe the roll-over and cutoff region of the ablation-front RT dispersion curve in indirect drive. Aluminum foils with imposed sinusoidal perturbations ranging in wavelength from 10 to 70 pm were ablatively accelerated with a radiation drive generated in a gold cylindrical hohlraum. A strong shock wave compresses the package followed by an {approx}2 ns period of roughly constant acceleration and the experiment is diagnosed via face-on radiography. Perturbations with wavelengths {ge} 20 {micro}m experienced substantial growth during the acceleration phase while shorter wavelengths showed a sharp drop off in overall growth. These experimental results compared favorably to calculations with a 2-D radiation-hydrodynamics code, however, the growth is significantly affected by the rippled shock launched …

In the frame of a CEA/US DOE collaboration, radiation driven spherically convergent experiments were performed on the Nova laser in order m measure the Rayleigh-Taylor growth at the ablation front. Numerical simulations using the 2D Lagrangian code FCI2 have correctly reproduced experiments in moderate convergent geometry. [C. Cherfils et al., PRL 83, 5507 (1999)]. Experiments have addressed convergence ratios up to 4 by considering larger capsules, larger hohlraum and longer laser pulses [S.G. Glendinning et al., to be published in Physics of Plasmas]. Numerical analysis of these high convergence implosions is presented, and the effect of convergence on the Rayleigh-Taylor growth is investigated.

Coscheduling is essential for obtaining good performance in a time-shared symmetric multiprocessor (SMP) cluster environment. However, the most common technique, gang scheduling, has limitations such as poor scalability and vulnerability to faults mainly due to explicit synchronization between its components. A decentralized approach called dynamic coscheduling (DCS) has been shown to be effective for network of workstations (NOW), but this technique is not suitable for the workloads on a very large SMP-cluster with thousands of processors. Furthermore, its implementation can be prohibitively expensive for such a large-scale machine. In this paper, we propose a novel coscheduling technique based on the DCS approach which can achieve coscheduling on very large SMP-clusters in a scalable, efficient, and cost-effective way. In the proposed technique, each local scheduler achieves coscheduling based upon message traffic between the components of parallel jobs. Message trapping is carried out at the user-level, eliminating the need for unsupported hardware or device-level programming. A sending process attaches its status to outgoing messages so local schedulers on remote nodes can make more intelligent scheduling decisions. Once scheduled, processes are guaranteed some minimum period of time to execute. This provides an opportunity to synchronize the parallel job's components across all nodes and …

Heavy ion beam transport through the containment chamber plays a crucial role in all heavy ion fusion (HIF) scenarios. Here, several parameters are used to characterize the operating space for HIF beams; transport modes are assessed in relation to evolving target/accelerator requirements; results of recent relevant experiments and simulations of HIF transport are summarized; and relevant instabilities are reviewed. All transport options still exist, including (1) vacuum ballistic transport, (2) neutralized ballistic transport, and (3) channel-like transport. Presently, the European HIF program favors vacuum ballistic transport, while the US HIF program favors neutralized ballistic transport with channel-like transport as an alternate approach. Further transport research is needed to clearly guide selection of the most attractive, integrated HIF system.

Numerical modeling of airflow and pollutant dispersion around buildings is a challenging task due to the geometrical variations of buildings and the extremely complex flow created by such surface-mounted obstacles. The airflow around buildings inevitably involves impingement and separation regions, a multiple vortex system with building wakes, and jetting effects in street canyons. The interference from adjacent buildings further complicates the flow and dispersion patterns. Thus accurate simulations of such flow and pollutant transport require not only appropriate physics submodels but also accurate numerics and significant computing resources. We have developed an efficient, high resolution CFD model for such purposes, with a primary goal to support incident response and preparedness in emergency response planning, vulnerability analysis, and the development of mitigation techniques.

The dynamics of particle flows in the DIII-D tokamak for two divertor configurations is considered. Fuel and intrinsic carbon impurity flows are analyzed using experimental data and 2D fluid plasma simulations. The flows in puff and pump experiments done in an open and a closed divertor geometry are described. It is shown that the flow of fuel particles is sensitive to divertor geometry. The pumping efficiency of the DIII-D cryopumps is a factor of 2 higher in a closed geometry than an open. The core refueling rate of an open divertor is a factor of 2 higher than that of a closed divertor. In contrast, the flow of impurity carbon particles is insensitive to divertor geometry. Both the core carbon content and the fraction of the carbon source which penetrates to the core is unchanged between an open and closed divertor. In addition, the core impurity content is found to be insensitive to the amplitude of gas puffing in the simulations.

We propose an alternative orthonormalization method that computes the orthonormal basis from the right singular vectors of a matrix. Its advantage are: (a) all operations are matrix-matrix multiplications and thus cache-efficient, (b) only one synchronization point is required in parallel implementations, (c) could be more stable than Gram-Schmidt. In addition, we consider the problem of incremental orthonormalization where a block of vectors is orthonormalized against a previously orthonormal set of vectors and among itself. We solve this problem by alternating iteratively between a phase of Gram-Schmidt and a phase of the new method. We provide error analysis and use it to derive bounds on how accurately the two successive orthonormalization phases should be performed to minimize total work performed. Our experiments confirm the favorable numerical behavior of the new method and its effectiveness on modern parallel computers.

None

An experimental and calculational study has been performed to understand and prevent inadvertent activation of the criticality alarm system (CAS) from fuel-handling operations at the Irradiated Fuel Storage Facility. In conjunction with the study, the CAS neutron detectors were tested to verify the design specifications for gamma rejection capability and zero response limit. A minimum physical restrictive boundary around the CAS location was established based on a gamma ray dose rate limit of 10 rad/hr. The canister loaded with spent nuclear fuel must be moved in the area outside the exclusion zone so as not to trigger a false alarm from the CAS detectors.

Amorphous carbon films with an s p{sup 3} content up to 25% and a negligible amount of hydrogen have been grown by evaporation of graphite and concurrent Ar{sup +} ion bombardment. The s p{sup 3} content is maximized for Ar{sup +} energies between 200 and 300 eV following a subplantation mechanism. Higher ion energies deteriorate the film due to sputtering and heating processes. The hardness of the films increases in the optimal assisting range from 8 to 18 GPa, and is explained by the crosslinking of graphitic planes through s p {sup 3} connecting site.

The structural properties, energetics, and dynamics of Ca{sup 2+} and Mn{sup 2+} substituents in KTaO{sub 3} are investigated from first principles. It is found that Ca substitutes for both K and Ta ions. Oxygen vacancies bind to isolated Ca ions residing at Ta-sites, causing off-center Ca displacement and forming large dipoles. There is also evidence that oppositely charged defects may cluster together. The calculations predict that the activation energy for dipole reorientation via oxygen vacancy hopping within the first neighbor shell of Ta-substituting Ca or Mn exceeds 2 eV. On the other hand, Mn{sup 2+} substituting at the K-site displaces off center along the (100) direction, also forming a dipole. This dipole can reorient via Mn hopping motion with an activation energy of {approximately} 0.18 eV, in reasonable agreement with experiments. The authors argue that, in general, metal ion hopping at the A-site, not oxygen vacancy hopping, is responsible for the small activation energies found in experiments.