This article is a response to an article by M. Adell et al. [Phy. Rev. Lett. 94, 139701 (2005)] about semiconductor-based spintronics research.

A drift shadow is an area immediately beneath an undergroundvoidthat, in theory, will be relatively drier than the surrounding rockmass. Numerical and analytical models of water flow through unsaturatedrock predict the existence of a drift shadow, but field tests confirmingits existence have yet to be performed. Proving the existence of driftshadows and understanding their hydrologic and transport characteristicscould provide a better understanding of how contaminants move in thesubsurface if released from waste emplacement drifts such as the proposednuclear waste repository at Yucca Mountain, Nevada. We describe the fieldprogram that will be used to investigate the existence of a drift shadowand the corresponding hydrological process at the Hazel-Atlas silica-sandmine located at the Black Diamond Mines Regional Preserve in Antioch,California. The location and configuration of this mine makes it anexcellent site to observe and measure drift shadow characteristics. Themine is located in a porous sandstone unit of the Domengine Formation, anapproximately 230 meter thick series of interbedded Eocene-age shales,coals, and massive-bedded sandstones. The mining method used at the minerequired the development of two parallel drifts, one above the other,driven along the strike of the mined sandstone stratum. Thisconfiguration provides the opportunity to introduce water into the rockmass in the upper drift and …

We present experimental data of electron energy distributions from ultra-intense (>10{sup 19} W/cm{sup 2}) laser-solid interactions using the Rutherford Appleton Laboratory Vulcan petawatt laser. These measurements were made using a CCD-based magnetic spectrometer. We present details on the distinct effective temperatures that were obtained for a wide variety of targets as a function of laser intensity. It is found that as the intensity increases from 10{sup 17} W/cm{sup 2} to 10{sup 19} W/cm{sup 2}, a 0.4 dependence on the laser intensity is found. Between 10{sup 19} W/cm{sup 2} and 10{sup 20} W/cm{sup 2}, a gradual rolling off of temperature with intensity is observed.

Z-Beamlet [1] is a single-beam high-energy Nd:glass laser used for backlighting high energy density (HED) plasma physics experiments at Sandia's Z-accelerator facility. The system currently generates a single backlit image per experiment, and has been employed on approximately 50% of Z-accelerator system shots in recent years. We have designed and are currently building a system that uses Z-Beamlet to generate two distinct backlit images with adjustable time delay ranging from 2 to 20 ns between frames. The new system will double the rate of data collection and allow the temporal evolution of high energy density phenomena to be recorded on a single shot.

Phylogenetic relationships among species of the salamanderfamily Salamandridae are investigated using nearly 3000 nucleotide basesof newly reported mitochondrial DNA sequence data from the mtDNA genicregion spanning the genes tRNALeu-COI. This study uses nearlycomprehensive species-level sampling to provide the first completephylogeny for the Salamandridae. Deep phylogenetic relationships amongthe three most divergent lineages in the family Salamandrina terdigitata,a clade comprising the "True" salamanders, and a clade comprising allnewts except S. terdigitata are difficult to resolve. However, mostrelationships within the latter two lineages are resolved with robustlevels of branch support. The genera Euproctus and Triturus arestatistically shown to be nonmonophyletic, instead each contains adiverse set of lineages positioned within the large newt clade. The genusParamesotriton is also resolve as a nonmonophyletic group, with the newlydescribed species P. laoensis constituting a divergent lineage placed ina sister position to clade containing all Pachytriton species and allremaining Paramesotriton species. Sequence divergences between P.laoensis and other Paramesotriton species are as great as those comparingP. laoensis and species of the genera Cynops and Pachytriton. Analyses oflineage diversification across the Salamandridae indicate that, despiteits exceptional diversity, lineage accumulation appears to have beenconstant across time, indicating that it does not represent a truespecies radiation.

The paramagnetic 1:1 coordination complexes of (C5Me5)2Ybwith a series of 4,4'-disubstituted bipyridines, bipy-X, where X is Me,tert-Bu, OMe, Ph, CO2Me, and CO2Et have been prepared. All of thecomplexes are paramagnetic and the values of the magnetic susceptibilityas a function of temperature show that these values are less thanexpected for the cation, [(C5Me5)2Yb(III)(bipy-X)]+, which have beenisolated as the cation-anion ion-pairs[(C5Me5)2Yb(III)(bipy-X)]+[(C5Me5)2YbI2]f fnfn where X is CO2Et, OMe andMe. The 1H NMR chemical shifts (293 K) for the methine resonances locatedat the 6,6' site in the bipy-X ring show a linear relationship with thevalues of chiT (300 K) for the neutral complexes which illustrates thatthe molecular behavior does not depend upon the phase with one exception,viz., (C5Me5)2Yb(bipy-Me). Single crystals of the 4,4'-dimethylbipyridinecomplex undergo an irreversible, abrupt first order phase change at 228 Kthat shatters the single crystals. The magnetic susceptibility,represented in a delta vs. T plot, on this complex, in polycrystallineform undergoes reversible abrupt changes in the temperature regime 205 -212 K, which is suggested to be due to the way the individual molec ularunits pack in the unit cell. A qualitative model is proposed thataccounts for the sub-normal magnetic moments in theseytterbocene-bipyridine complexes.

The Kinetic Stabilizer (K-S) concept [1] represents a means for stabilizing axisymmetric mirror and tandem-mirror (T-M) magnetic fusion systems against MHD interchange instability modes. Magnetic fusion research has given us examples of axisymmetric mirror confinement devices in which radial transport rates approach the classical ''Spitzer'' level, i.e. situations in which turbulence if present at all, is at too low a level to adversely affect the radial transport [2,3,4]. If such a low-turbulence condition could be achieved in a T-M system it could lead to a fusion power system that would be simpler, smaller, and easier to develop than one based on closed-field confinement, e.g., the tokamak, where the transport is known to be dominated by turbulence. However, since conventional axisymmetric mirror systems suffer from the MHD interchange instability, the key to exploiting this new opportunity is to find a practical way to stabilize this mode. The K-S represents one avenue to achieving this goal. The starting point for the K-S concept is a theoretical analysis by Ryutov [5]. He showed that a MHD-unstable plasma contained in an axisymmetric mirror cell can be MHD-stabilized by the presence of a low-density plasma on the expanding field lines outside the mirrors. If this …

Interactions between surface and ground water are a key component of the hydrologic budget on the watershed scale. Models that honor these interactions are commonly based on the conductance concept that presumes a distinct interface at the land surface, separating the surface from the subsurface domain. These types of models link the subsurface and surface domains via an exchange flux that depends upon the magnitude and direction of the hydraulic gradient across the interface and a proportionality constant (a measure of the hydraulic connectivity). Because experimental evidence of such a distinct interface is often lacking in field systems, there is a need for a more general coupled modeling approach. A more general coupled model is presented that incorporates a new two-dimensional overland flow simulator into the parallel three-dimensional variable saturated subsurface flow code ParFlow. In ParFlow, the overland flow simulator takes the form of an upper boundary condition and is, thus, fully integrated without relying on the conductance concept. Another important advantage of this approach is the efficient parallelism incorporated into ParFlow, which is efficiently exploited by the overland flow simulator. Several verification and simulation examples are presented that focus on the two main processes of runoff production: excess infiltration …

We have synthesized GaN{sub x}As{sub 1-y}P{sub y} alloys (x {approx} 0.3-1% and y = 0-0.4) using nitrogen N ion implantation into GaAsP epilayers followed by pulsed laser melting and rapid thermal annealing techniques. As predicted by the band anticrossing model, the incorporation of N splits the conduction band (E{sub M}) of the GaAs{sub 1-y}P{sub y} substrate, and strong optical transitions from the valence band to the lower (E{sub -}) and upper (E{sub +}) conduction subbands are observed. The relative strengths of the E{sub -} and E{sub +} transition change as the localized N level E{sub N} emerges from the conduction band forming narrow intermediate band for y > 0.3. The results show that GaN{sub x}As{sub 1-x-y}P{sub y} alloys with y > 0.3 is a three band semiconductor alloy with potential applications for high-efficiency intermediate band solar cells.

A new method for representing subgrid-scale cloud structure, in which each model column is decomposed into a set of subcolumns, has been introduced into the Geophysical Fluid Dynamics Laboratory's global climate model AM2. Each subcolumn in the decomposition is homogeneous but the ensemble reproduces the initial profiles of cloud properties including cloud fraction, internal variability (if any) in cloud condensate, and arbitrary overlap assumptions that describe vertical correlations. These subcolumns are used in radiation and diagnostic calculations, and have allowed the introduction of more realistic overlap assumptions. This paper describes the impact of these new methods for representing cloud structure in instantaneous calculations and long-term integrations. Shortwave radiation computed using subcolumns and the random overlap assumption differs in the global annual average by more than 4 W/m{sup 2} from the operational radiation scheme in instantaneous calculations; much of this difference is counteracted by a change in the overlap assumption to one in which overlap varies continuously with the separation distance between layers. Internal variability in cloud condensate, diagnosed from the mean condensate amount and cloud fraction, has about the same effect on radiative fluxes as does the ad hoc tuning accounting for this effect in the operational radiation scheme. Long …

Hiding communication latency is an important optimization for parallel programs. Programmers or compilers achieve this by using non-blocking communication primitives and overlapping communication with computation or other communication operations. Using non-blocking communication raises two issues: performance and programmability. In terms of performance, optimizers need to find a good communication schedule and are sometimes constrained by lack of full application knowledge. In terms of programmability, efficiently managing non-blocking communication can prove cumbersome for complex applications. In this paper we present the design principles of HUNT, a runtime system designed to search and exploit some of the available overlap present at execution time in UPC programs. Using virtual memory support, our runtime implements demand-driven synchronization for data involved in communication operations. It also employs message decomposition and scheduling heuristics to transparently improve the non-blocking behavior of applications. We provide a user level implementation of HUNT on a variety of modern high performance computing systems. Results indicate that our approach is successful in finding some of the overlap available at execution time. While system and application characteristics influence performance, perhaps the determining factor is the time taken by the CPU to execute a signal handler. Demand driven synchronization at execution time eliminates the …

We performed ab initio quantum chemical studies for the development of intra and intermolecular interaction potentials for formic acid for use in molecular dynamics simulations of formic acid molecular crystal. The formic acid structures considered in the ab initio studies include both the cis and trans monomers which are the conformers that have been postulated as part of chains constituting liquid and crystal phases under extreme conditions. Although the cis to trans transformation is not energetically favored, the trans isomer was found as a component of stable gas-phase species. Our decomposition scheme for the interaction energy indicates that the hydrogen bonded complexes are dominated by the Hartree-Fock forces while parallel clusters are stabilized by the electron correlation energy. The calculated three-body and higher interactions are found to be negligible, thus rationalizing the development of an atom-atom pair potential for formic acid based on high-level ab initio calculations of small formic acid clusters. Here we present an atom-atom pair potential that includes both intra- and inter-molecular degrees of freedom for formic acid. The newly developed pair potential is used to examine formic acid in the condensed phase via molecular dynamics simulations. The isothermal compression under hydrostatic pressure obtained from molecular dynamics …

In Biot's theory of poroelasticity, elastic materials contain connected voids or pores and these pores may be filled with fluids under pressure. The fluid pressure then couples to the mechanical effects of stress or strain applied externally to the solid matrix. Eshelby's formula for the response of a single ellipsoidal elastic inclusion in an elastic whole space to a strain imposed at a distant boundary is a very well-known and important result in elasticity. Having a rigorous generalization of Eshelby's results valid for poroelasticity means that the hard part of Eshelby's work (in computing the elliptic integrals needed to evaluate the fourth-rank tensors for inclusions shaped like spheres, oblate and prolate spheroids, needles and disks) can be carried over from elasticity to poroelasticity--and also thermoelasticity--with only relatively minor modifications. Effective medium theories for poroelastic composites such as rocks can then be formulated easily by analogy to well-established methods used for elastic composites. An identity analogous to Eshelby's classic result has been derived [Physical Review Letters 79:1142-1145 (1997)] for use in these more complex and more realistic problems in rock mechanics analysis. Descriptions of the application of this result as the starting point for new methods of estimation are presented, including …

Large scale molecular dynamics simulations of bulk melts of polar (poly(vinylidene fluoride) (pVDF)) polymers are utilized to study chain conformation and ordering prior to crystallization under cooling. While the late stages of polymer crystallization have been studied in great detail, recent theoretical and experimental evidence indicates that there are important phenomena occurring in the early stages of polymer crystallization that are not understood to the same degree. When the polymer melt is quenched from a temperature above the melting temperature to the crystallization temperature, crystallization does not occur instantaneously. This initial interval without crystalline order is characterized as an induction period. It has been thought of as a nucleation period in the classical theories of polymer crystallization, but recent experiments, computer simulations, and theoretical work suggest that the initial period in polymer crystallization is assisted by a spinodal decomposition type mechanism. In this study we have achieved physically realistic length scales to study early stages of polymer ordering, and show that spinodal-assisted ordering prior to crystallization is operative in polar polymers suggesting general applicability of this process.

To enable high-energy petawatt laser operation we have developed the processing methods and tooling that produced both the world's largest multilayer dielectric reflection grating and the world's highest laser damage resistant gratings. We have successfully delivered the first ever 80 cm aperture multilayer dielectric grating to LLNL's Titan Intense Short Pulse Laser Facility. We report on the design, fabrication and characterization of multilayer dielectric diffraction gratings.

Three specific failure criteria for the transversely isotropic fiber composite case will be discussed. All three use the polynomial expansion method. The three criteria are the Tsai-Wu criterion, the Hashin criterion and the Christensen criterion. All three criteria will be given in forms that admit direct and easy comparison, which has not usually been done. The central differences between these three criteria will be discussed, and steps will be taken toward the evaluation of them.

The inductive coupling between various MICE magnet circuits is described. The consequences of this coupling on magnet charging and quenching are discussed. Magnet quench protection is achieved through the use of quench-back. Calculations of the quenching of a magnet due to quench-back resulting from circulating currents induced in the magnet mandrel due to quenching of an adjacent magnet are discussed. This report describes how the MICE magnet channel will react when magnets in that channel are quenched.

As low temperature cryocoolers become more frequently used to cool superconducting magnets, it becomes increasingly apparent that the connection between the cooler and the magnet has an effect on the design and performance of the magnet. In general, the use of small coolers can be considered in two different temperature ranges; (1) from 3.8 to 4.8 K for magnet fabricated with LTS conductor and (2) from 18 to 35 K for magnets fabricated using HTS conductor. In general, both temperature ranges call for the use of a two-stage cooler. The best method for connecting a cooler to the magnet depends on a number of factors. The factors include: (1) whether the cooler must be used to cool down the magnet from room temperature, (2) whether the magnet must have one or more reservoirs of liquid cryogen to keep the magnet cold during a loss of cooling, and (3) constraints on the distance from the cooler cold heads and the magnet and its shield. Two methods for connecting low temperature coolers to superconducting magnets have been studied. The first method uses a cold strap to connect the cold heads directly to the loads. This method is commonly used for cryogen-free magnets. …

We describe the progress on the design of the RF coupling coil (RFCC) module for the international Muon Ionization Cooling Experiment (MICE) at Rutherford Appleton Laboratory (RAL) in the UK. The MICE cooling channel design consists of one SFOFO cell that is similar to that of the US Study-II of a neutrino factory. The MICE RFCC module comprises a superconducting solenoid, mounted around four normal conducting 201.25-MHz RF cavities. Each cavity has a pair of thin curved beryllium windows to close the conventional open beam irises, which allows for independent control of the phase in each cavity and for the RF power to be fed separately. The coil package that surrounds the RF cavities is mounted on a vacuum vessel. The RF vacuum is shared between the cavities and the vacuum vessel around the cavities such that there is no differential pressure on the thin beryllium windows. This paper discusses the design progress of the RFCC module and the fabrication progress of a prototype 201.25-MHz cavity.

The transition from arid glacial to moist early Holocene conditions represented a profound change in northern lowland Neotropical climate. Here we report a detailed record of changes in moisture availability during the latter part of this transition ({approx}11,250 to 7,500 cal yr BP) inferred from sediment cores retrieved in Lake Peten Itza, northern Guatemala. Pollen assemblages demonstrate that a mesic forest had been largely established by {approx}11,250 cal yr BP, but sediment properties indicate that lake level was more than 35 m below modern stage. From 11,250 to 10,350 cal yr BP, during the Preboreal period, lithologic changes in sediments from deep-water cores (>50 m below modern water level) indicate several wet-dry cycles that suggest distinct changes in effective moisture. Four dry events (designated PBE1-4) occurred at 11,200, 10,900, 10,700, and 10,400 cal yr BP and correlate with similar variability observed in the Cariaco Basin titanium record and glacial meltwater pulses into the Gulf of Mexico. After 10,350 cal yr BP, multiple sediment proxies suggest a shift to a more persistently moist early Holocene climate. Comparison of results from Lake Peten Itza with other records from the circum-Caribbean demonstrates a coherent climate response during the entire span of our record. …

Habitat use and food web linkages are critical data for fish conservation and habitat restoration efforts, particularly for threatened salmonids species. Otolith microchemistry has been shown to be a powerful tool for reconstructing fish movement, but over small distances (kilometers), geology-derived differences in otolith chemistry are rare. Here, we demonstrate that otolith {sup 13}C/{sup 12}C ratio (i.e. {delta}{sup 13}C) of anadromous steelhead trout can be used to distinguish residence in small streams from residence in larger streams and rivers. While previous research has shown that water dissolved inorganic carbon {delta}{sup 13}C is the primary source of carbon in otoliths, the downstream change in food {delta}{sup 13}C in this watershed appears to be the primary control on otolith {delta}{sup 13}C. As a result, this method can also be applied to the problem of reconstructing feeding history at a location.

Oceanic uptake and transport of bomb radiocarbon as {sup 14}CO{sub 2} created by atmospheric nuclear weapons testing in the 1950s and 1960s has been a useful diagnostic to determine the carbon transfer between the ocean and atmosphere. In addition, the distribution of radiocarbon in the ocean can be used as a tracer of oceanic circulation. Results obtained from samples collected in the Gulf of Alaska in the summer of 2002 provide a direct comparison with results in the 1970s during GEOSECS and in the early 1990s during WOCE. The open gyre values are 20-40{per_thousand} more negative than those documented in 1991 and 1993 (WOCE) although the general trends as a function of latitude are reproduced. Surface values are still significantly higher than pre-bomb levels ({approx}-105{per_thousand} or lower). In the central gyre, we observe {Delta}{sup 14}C-values that are lower in comparison to GEOSECS (stn 218) and WOCE P16/P17 to a density of {approx}26.8{sigma}t. This observation is consistent with the overall decrease in surface {Delta}{sup 14}C values, and reflects the erosion of the bomb-{sup 14}C transient. We propose that erosion of the bomb-{sup 14}C transient is accomplished by entrainment of low {sup 14}C water via vertical exchange within the Gulf of Alaska …

The authors present molecular dynamics simulations of aggregation kinetics in a colloidal suspension modeled as a highly asymmetric binary mixture. Starting from a configuration with largely uncorrelated colloidal particles the system relaxes by coagulation-fragmentation dynamics to a structured state of low-dimensionality clusters with an exponential size distribution. The results show that short range repulsive interactions alone can give rise to so-called cluster phases. For the present model and probably other, more common colloids, the observed clusters appear to be equilibrium phase fluctuations induced by the entropic inter-colloidal attractions.

None