Compositional data are represented as vector variables with individual vector components ranging between zero and a positive maximum value representing a constant sum constraint, usually unity (or 100 percent). The earth sciences are flooded with spatial distributions of compositional data, such as concentrations of major ion constituents in natural waters (e.g. mole, mass, or volume fractions), mineral percentages, ore grades, or proportions of mutually exclusive categories (e.g. a water-oil-rock system). While geostatistical techniques have become popular in earth science applications since the 1970s, very little attention has been paid to the unique mathematical properties of geostatistical formulations involving compositional variables. The book 'Geostatistical Analysis of Compositional Data' by Vera Pawlowsky-Glahn and Ricardo Olea (Oxford University Press, 2004), unlike any previous book on geostatistics, directly confronts the mathematical difficulties inherent to applying geostatistics to compositional variables. The book righteously justifies itself with prodigious referencing to previous work addressing nonsensical ranges of estimated values and error, spurious correlation, and singular cross-covariance matrices.

In the context of eternal inflation we discuss the fate of Lambda = 0 bubbles when they collide with Lambda< 0 crunching bubbles. When the Lambda = 0 bubble is supersymmetric, it is not completely destroyed by collisions. If the domain wall separating the bubbles has higher tension than the BPS bound, it is expelled from the Lambda = 0 bubble and does not alter its long time behavior. If the domain wall saturates the BPS bound, then it stays inside the Lambda = 0 bubble and removes a finite fraction of future infinity. In this case, the crunch singularity is hidden behind the horizon of a stable hyperbolic black hole.

The evaporation rate of D{sub 2}O has been determined by Raman thermometry of a droplet train (12-15 {micro}m diameter) injected into vacuum ({approx}10{sup -5} torr). The cooling rate measured as a function of time in vacuum was fit to a model that accounts for temperature gradients between the surface and the core of the droplets, yielding an evaporation coefficient ({gamma}{sub e}) of 0.57 {+-} 0.06. This is nearly identical to that found for H{sub 2}O (0.62 {+-} 0.09) using the same experimental method and model, and indicates the existence of a kinetic barrier to evaporation. The application of a recently developed transition state theory (TST) model suggests that the kinetic barrier is due to librational and hindered translational motions at the liquid surface, and that the lack of an isotope effect is due to competing energetic and entropic factors. The implications of these results for cloud and aerosol particles in the atmosphere are discussed.

None

We investigate the climate forcing from and response to projected changes in short-lived species and methane under the A1B scenario from 2000-2050 in the GISS climate model. We present a meta-analysis of new simulations of the full evolution of gas and aerosol species and other existing experiments with variations of the same model. The comparison highlights the importance of several physical processes in determining radiative forcing, especially the effect of climate change on stratosphere-troposphere exchange, heterogeneous sulfate-nitrate-dust chemistry, and changes in methane oxidation and natural emissions. However, the impact of these fairly uncertain physical effects is substantially less than the difference between alternative emission scenarios for all short-lived species. The net global mean annual average direct radiative forcing from the short-lived species is .02 W/m{sup 2} or less in our projections, as substantial positive ozone forcing is largely offset by negative aerosol direct forcing. Since aerosol reductions also lead to a reduced indirect effect, the global mean surface temperature warms by {approx}0.07 C by 2030 and {approx}0.13 C by 2050, adding 19% and 17%, respectively, to the warming induced by long-lived greenhouse gases. Regional direct forcings are large, up to 3.8 W/m{sup 2}. The ensemble-mean climate response shows little regional …

Direct laser heating of molecular N2 to above 1400 K at 120-130 GPa results in the formation of a reddish amorphous phase and a transparent crystalline solid above 2000 K. Raman and x-ray data confirm that the transparent phase is cubic-gauche nitrogen (cg-N), while the reddish color of the amorphous phase might indicate the presence of N=N dish bonds. The quenched amorphous phase is stable down to at least 70GPa, analogous to cg-N, and could be a new non-molecular phase or an extension of the already known {eta}-phase. A chemo-physical phase diagram is presented which emphasizes the difference between pressure- and temperature-induced transitions from molecular to non-molecular solids, as found in other low Z systems.

We study the energy dependence of the transverse momentum (pT) spectra for charged pions, protons and anti-protons for Au+Au collisions at sqrt sNN = 62.4 and 200 GeV. Data are presented at mid-rapidity (lbar y rbar< 0.5) for 0.2< pT< 12 GeV/c. In the intermediate pT region (2< pT< 6 GeV/c), the nuclear modification factor is higher at 62.4 GeV than at 200 GeV, while at higher pT (pT> 7 GeV/c) the modification is similar for both energies. The p/pi+ and pbar/pi- ratios for central collisions at sqrt sNN = 62.4 GeV peak at pT _~;; 2 GeV/c. In the pT range where recombination is expected to dominate, the p/pi+ ratios at 62.4 GeV are larger than at 200 GeV, while the pbar/pi- ratios are smaller. For pT> 2 GeV/c, the pbar/pi- ratios at the two beam energies are independent of pT and centrality indicating that the dependence of the pbar/pi- ratio on pT does not change between 62.4 and 200 GeV. These findings challenge various models incorporating jet quenching and/or constituent quark coalescence.

We develop a new ab initio many-body approach capable of describing simultaneously both bound and scattering states in light nuclei, by combining the resonating-group method with the use of realistic interactions, and a microscopic and consistent description of the nucleon clusters. This approach preserves translational symmetry and Pauli principle. We present phase shifts for neutron scattering on {sup 3}H, {sup 4}He and {sup 10}Be and proton scattering on {sup 3,4}He, using realistic nucleon-nucleon potentials. Our A = 4 scattering results are compared to earlier ab initio calculations. We demonstrate that a proper treatment of the coupling to the n-{sup 10}Be continuum is essential to explain the parity-inverted ground state in {sup 11}Be.

In contrast to failure approaches at the lamina level or the micromechanics level the present work concerns failure characterization at the laminate level. Specifically, attention is given to the ultimate failure characterization for quasi-isotropic laminates. This is in further contrast to the commonly used approaches for initial damage or progressive damage. It is shown that the analytical failure forms decompose into two modes, one for out of plane, delamination type failure and one for in plane, fiber controlled type failure. The work here is mainly given over to the delamination mode of failure. Experimental results are presented for laminates in this mode of failure. These results are then integrated with the analytical forms to give a simple criterion for delamination failure.

The usual approach to deal with noise present in many real-world optimization problems is to take an arbitrary number of samples of the objective function and use the sample average as an estimate of the true objective value. The number of samples is typically chosen arbitrarily and remains constant for the entire optimization process. This paper studies an adaptive sampling technique that varies the number of samples based on the uncertainty of deciding between two individuals. Experiments demonstrate the effect of adaptive sampling on the final solution quality reached by a genetic algorithm and the computational cost required to find the solution. The results suggest that the adaptive technique can effectively eliminate the need to set the sample size a priori, but in many cases it requires high computational costs.

QR and LU factorizations for dense matrices are important linear algebra computations that are widely used in scientific applications. To efficiently perform these computations on modern computers, the factorization algorithms need to be blocked when operating on large matrices to effectively exploit the deep cache hierarchy prevalent in today's computer memory systems. Because both QR (based on Householder transformations) and LU factorization algorithms contain complex loop structures, few compilers can fully automate the blocking of these algorithms. Though linear algebra libraries such as LAPACK provides manually blocked implementations of these algorithms, by automatically generating blocked versions of the computations, more benefit can be gained such as automatic adaptation of different blocking strategies. This paper demonstrates how to apply an aggressive loop transformation technique, dependence hoisting, to produce efficient blockings for both QR and LU with partial pivoting. We present different blocking strategies that can be generated by our optimizer and compare the performance of auto-blocked versions with manually tuned versions in LAPACK, both using reference BLAS, ATLAS BLAS and native BLAS specially tuned for the underlying machine architectures.

Nestled in a valley between the whitecaps of the Pacific and the snowcapped crests of the Sierra Nevada, Lawrence Livermore National Laboratory (LLNL) is home to the nearly complete National Ignition Facility (NIF). The purpose of NIF is to create a miniature star-on demand. An enormous amount of laser light energy (1.8 MJ in a pulse that is 20 ns in duration) will be focused into a small gold cylinder approximately the size of a pencil eraser. Centered in the gold cylinder (or hohlraum) will be a nearly perfect sphere filled with a complex mixture of hydrogen gas isotopes that is similar to the atmosphere of our Sun. During experiments, the laser light will hit the inside of the gold cylinder, heating the metal until it emits X-rays (similar to how your electric stove coil emits visible red light when heated). The X-rays will be used to compress the hydrogen-like gas with such pressure that the gas atoms will combine or 'fuse' together, producing the next heavier element (helium) and releasing energy in the form of energetic particles. 2010 will mark the first credible attempt at this world-changing event: the achievement of fusion energy 'break-even' on Earth using NIF, the …

The feasibility to coat large SNF/HLW containers with a structurally amorphous material (SAM) was demonstrated on sub-scale models fabricated from Type 316L stainless steel. The sub-scale model were coated with SAM 1651 material using kerosene high velocity oxygen fuel (HVOF) torch to thicknesses ranging from 1 mm to 2 mm. The process parameters such as standoff distance, oxygen flow, and kerosene flow, were optimized in order to improve the corrosion properties of the coatings. Testing in an electrochemical cell and long-term exposure to a salt spray environment were used to guide the selection of process parameters.

The science of complexity has moved to center stage within the past few decades. Complex systems range from glasses to the immune system and the brain. Glasses are too simple to possess all aspects of complexity; brains are too complex to expose common concepts and laws of complexity. Proteins, however, are systems where many concepts and laws of complexity can be explored experimentally, theoretically, and computationally. Such studies have elucidated crucial aspects. The energy landscape has emerged as one central concept; it describes the free energy of a system as a function of temperature and the coordinates of all relevant atoms. A second concept is that of fluctuations. Without fluctuations, proteins would be dead and life impossible. A third concept is slaving. Proteins are not isolated systems; they are embedded in cells and membranes. Slaving arises when the fluctuations in the surroundings of a protein dominate many of the motions of the protein proper.

We study the Collins mechanism for the single transverse spin asymmetry in the collinear factorization approach. The correspondent twist-three fragmentation function is identified. We show that the Collins function calculated in this approach is universal.We further examine its contribution to the single transverse spin asymmetry of semi-inclusive hadron production in deep inelastic scattering and demonstrate that the transverse momentum dependent and twist-three collinear approaches are consistent in the intermediate transverse momentum region where both apply.

The objective of this paper is to share the SRS methodology for identifying the reliability requirements and documenting the expected performance of Safety Instrumented Functions (SIFs) used as criticality defenses. Nuclear Criticality SIFs are comprised of sensors, logic solvers, and final control elements, which may be either automatic or manual, to detect a process hazard and respond to prevent a criticality. The Savannah River Site (SRS) has invoked the chemical process industry safety standard (ANSI/ISA 84.00.01) for the design of safety significant instrumented systems. The ISA standard provides a graded approach to design based on the amount of risk reduction that is required of an SIF. SRS is embarking on application of this standard to nuclear criticality defenses, thus integrating criticality safety requirements with verifiable design methodology. Per the DOE G 421.1-1 discussion of the double contingency principle, guidance for a single contingency barrier includes, ''The estimated probability that the control will fail (when called upon for protection) is not greater than 1 in 100 demands''. The application of this standard to nuclear criticality SIFs will provide clear requirements in terms of safety availability and testing to assure that the instrumented criticality system as designed, installed, and maintained will meet …

We perform a maximum likelihood analysis of the cluster abundance measured in the SDSS using the maxBCG cluster finding algorithm. Our analysis is aimed at constraining the power spectrum normalization {sigma}{sub 8}, and assumes flat cosmologies with a scale invariant spectrum, massless neutrinos, and CMB and supernova priors {Omega}{sub m}h{sup 2} = 0.128 {+-} 0.01 and h = 0.72 {+-} 0.05 respectively. Following the method described in the companion paper Rozo et al. (2007), we derive {sigma}{sub 8} = 0.92 {+-} 0.10 (1{sigma}) after marginalizing over all major systematic uncertainties. We place strong lower limits on the normalization, {sigma}{sub 8} > 0.76 (95% CL) (> 0.68 at 99% CL). We also find that our analysis favors relatively low values for the slope of the Halo Occupation Distribution (HOD), {alpha} = 0.83 {+-} 0.06. The uncertainties of these determinations will substantially improve upon completion of an ongoing campaign to estimate dynamical, weak lensing, and X-ray cluster masses in the SDSS maxBCG cluster sample.

For the licensing of the current fleet of commercial nuclear power plants (NPPs), the Nuclear Regulatory Commission (NRC) used two key documents, NUREG-0800 and Regulatory Guide (RG) 1.70. RG 1.70 provided guidance to applicants on the contents needed in their Safety Analysis Reports (SARs) submitted as part of their application to construct or operate an NPP. NUREG-0800, the NRC Standard Review Plan (SRP), provides guidance to the NRR staff reviewers on performing their safety reviews of these applications. As part of the preparation for a new wave of improved NPP designs the NRC is in the process of updating the SRP and is also developing a new RG designated as draft RG or DG-1145, ''Combined License Applications for Nuclear Power Plants (LWR Edition).'' This will eventually become RG 1.206 and will take the place of RG 1.70. This will provide guidance for combined license (COL) applicants, as well as for other 10CFR Part 52 variations that are permitted.

I compute the v/c correction to the gravitational time delayfor light passing by a massive object moving with speed v, and I finddisagreement with previously published results. It is also argued thatthe speed of gravity formula that was recently used in the conjunction ofJupiter and quasar J0842+1845 is frame dependent.

At high temperatures or densities matter formed by strongly interacting elementary particles (hadronic matter) is expected to undergo a transition to a new form of matter--the quark gluon plasma--in which elementary particles (quarks and gluons) are no longer confined inside hadrons but are free to propagate in a thermal medium much larger in extent than the typical size of a hadron. The transition to this new form of matter as well as properties of the plasma phase are studied in large scale numerical calculations based on the theory of strong interactions--Quantum Chromo Dynamics (QCD). Experimentally properties of hot and dense elementary particle matter are studied in relativistic heavy ion collisions such as those currently performed at the relativistic heavy ion collider (RHIC) at BNL. We review here recent results from studies of thermodynamic properties of strongly interacting elementary particle matter performed on Teraflops-Computer. We present results on the QCD equation of state and discuss the status of studies of the phase diagram at non-vanishing baryon number density.

The effects of fast neutron irradiation on SiC and SiC composites have been studied. The materials used were chemical vapor deposition (CVD) SiC and SiC/SiC composites reinforced with either Hi-Nicalon{trademark} Type-S, Hi-Nicalon{trademark} or Sylramic{trademark} fibers fabricated by chemical vapor infiltration. Statistically significant numbers of flexural samples were irradiated up to 4.6 x 10{sup 25} n/m{sup 2} (E>0.1 MeV) at 300, 500 and 800 C in the High Flux Isotope Reactor at Oak Ridge National Laboratory. Dimensions and weights of the flexural bars were measured before and after the neutron irradiation. Mechanical properties were evaluated by four point flexural testing. Volume increase was seen for all bend bars following neutron irradiation. Magnitude of swelling depended on irradiation temperature and material, while it was nearly independent of irradiation fluence over the fluence range studied. Flexural strength of CVD SiC increased following irradiation depending on irradiation temperature. Over the temperature range studied, no significant degradation in mechanical properties was seen for composites fabricated with Hi-Nicalon{trademark} Type-S, while composites reinforced with Hi-Nicalon{trademark} or Sylramic fibers showed significant degradation. The effects of irradiation on the Weibull failure statistics are also presented suggesting a reduction in the Weibull modulus upon irradiation. The cause of this potential …

None

We introduce a framework for describing the halo selection function of optical cluster finders. We treat the problem as being separable into a term that describes the intrinsic galaxy content of a halo (the Halo Occupation Distribution, or HOD) and a term that captures the effects of projection and selection by the particular cluster finding algorithm. Using mock galaxy catalogs tuned to reproduce the luminosity dependent correlation function and the empirical color-density relation measured in the SDSS, we characterize the maxBCG algorithm applied by Koester et al. to the SDSS galaxy catalog. We define and calibrate measures of completeness and purity for this algorithm, and demonstrate successful recovery of the underlying cosmology and HOD when applied to the mock catalogs. We identify principal components--combinations of cosmology and HOD parameters--that are recovered by survey counts as a function of richness, and demonstrate that percent-level accuracies are possible in the first two components, if the selection function can be understood to {approx} 15% accuracy.

In recent years, one of the advances of the World Wide Web is social media and one of the fastest growing aspects of social media is the blogosphere. Blogs make content creation easy and are highly accessible through web pages and syndication. With their growing influence, a need has arisen to be able to monitor the opinions and insight revealed within their content. In this paper we describe a technical approach for analyzing the content of blog data using a visual analytic tool, IN-SPIRE, developed by Pacific Northwest National Laboratory. We highlight the capabilities of this tool that are particularly useful for information gathering from blog data.