The balancing domain decomposition methods by constraints are extended to solving both nonsymmetric, positive definite and symmetric, indefinite linear systems. In both cases, certain nonstandard primal constraints are included in the coarse problems of BDDC algorithms to accelerate the convergence. Under the assumption that the subdomain size is small enough, a convergence rate estimate for the GMRES iteration is established that the rate is independent of the number of subdomains and depends only slightly on the subdomain problem size. Numerical experiments for several two-dimensional examples illustrate the fast convergence of the proposed algorithms.

This paper assesses the accuracy of the simplified frame cavity conduction/convection and radiation models presented in ISO 15099 and used in software for rating and labeling window products. Temperatures and U-factors for typical horizontal window frames with internal cavities are compared; results from Computational Fluid Dynamics (CFD) simulations with detailed radiation modeling are used as a reference. Four different frames were studied. Two were made of polyvinyl chloride (PVC) and two of aluminum. For each frame, six different simulations were performed, two with a CFD code and four with a building-component thermal-simulation tool using the Finite Element Method (FEM). This FEM tool addresses convection using correlations from ISO 15099; it addressed radiation with either correlations from ISO 15099 or with a detailed, view-factor-based radiation model. Calculations were performed using the CFD code with and without fluid flow in the window frame cavities; the calculations without fluid flow were performed to verify that the CFD code and the building-component thermal-simulation tool produced consistent results. With the FEM-code, the practice of subdividing small frame cavities was examined, in some cases not subdividing, in some cases subdividing cavities with interconnections smaller than five millimeters (mm) (ISO 15099) and in some cases subdividing cavities …

Chemical Cleaning is currently in progress for Tanks 5 and 6 at the Savannah River Site. The Chemical Cleaning process is being utilized to remove the residual waste heel remaining after completion of Mechanical Sludge Removal. This work is required to prepare the tanks for closure. Tanks 5 and 6 are 1950s vintage carbon steel waste tanks that do not meet current containment standards. These tanks are 22.9 meters (75 feet) in diameter, 7.5 meters (24.5 feet) in height, and have a capacity of 2.84E+6 liters (750,000 gallons). Chemical Cleaning adds 8 wt % oxalic acid to the carbon steel tank to dissolve the remaining sludge heel. The resulting acidic waste solution is transferred to Tank 7 where it is pH adjusted to minimize corrosion of the carbon steel tank. The Chemical Cleaning flowsheet includes multiple strikes of acid in each tank. Acid is delivered by tanker truck and is added to the tanks through a hose assembly connected to a pipe penetration through the tank top. The flowsheet also includes spray washing with acid and water. This paper includes an overview of the configuration required for Chemical Cleaning, the planned flowsheet, and an overview of technical concerns associated with …

We present new results of our program to systematically search for strongly lensed galaxies in the Sloan Digital Sky Survey (SDSS) imaging data. In this study six strong lens systems are presented which we have confirmed with followup spectroscopy and imaging using the 3.5m telescope at the Apache Point Observatory. Preliminary mass models indicate that the lenses are group-scale systems with velocity dispersions ranging from 466?878 km s{sup -1} at z = 0.17-0.45 which are strongly lensing source galaxies at z = 0.4-1.4. Galaxy groups are a relatively new mass scale just beginning to be probed with strong lensing. Our sample of lenses roughly doubles the confirmed number of group-scale lenses in the SDSS and complements ongoing strong lens searches in other imaging surveys such as the CFHTLS (Cabanac et al. 2007). As our arcs were discovered in the SDSS imaging data they are all bright (r {approx_equal} 22), making them ideally suited for detailed follow-up studies.

In order to reach the high peak current required for an x-ray free electron laser, two separate magnetic dipole chicanes are used in the Linac Coherent Light Source (LCLS) accelerator to compress the electron bunch length in stages. In these bunch compressors, coherent synchrotron radiation (CSR) can be emitted either by a short electron bunch or by any longitudinal density modulation that may be on the bunch. In this paper, we report detailed measurements of the CSR-induced energy loss and transverse emittance growth in these compressors. Good agreement is found between the experimental results and multi-particle tracking studies. We also describe direct observations of CSR at optical wavelengths and compare with analytical models based on beam microbunching.

Reynolds lubrication approximation is used extensively to study flows between moving machine parts, in narrow channels, and in thin films. The solution of Reynolds equation may be thought of as the zeroth order term in an expansion of the solution of the Stokes equations in powers of the aspect ratio {var_epsilon} of the domain. In this paper, we show how to compute the terms in this expansion to arbitrary order on a two-dimensional, x-periodic domain and derive rigorous, a-priori error bounds for the difference between the exact solution and the truncated expansion solution. Unlike previous studies of this sort, the constants in our error bounds are either independent of the function h(x) describing the geometry, or depend on h and its derivatives in an explicit, intuitive way. Specifically, if the expansion is truncated at order 2k, the error is O({var_epsilon}{sup 2k+2}) and h enters into the error bound only through its first and third inverse moments {integral}{sub 0}{sup 1} h(x){sup -m} dx, m = 1,3 and via the max norms {parallel} 1/{ell}! h{sup {ell}-1}{partial_derivative}{sub x}{sup {ell}}h{parallel}{sub {infinity}}, 1 {le} {ell} {le} 2k + 2. We validate our estimates by comparing with finite element solutions and present numerical evidence that suggests …

This paper discusses the application of existing workflow management systems to a real world science application (LQCD). Typical workflows and execution environment used in production are described. Requirements for the LQCD production system are discussed. The workflow management systems Askalon and Swift were tested by implementing the LQCD workflows and evaluated against the requirements. We report our findings and future work.

We explore the physics of the general CP-conserving MSSM with Minimal Flavor Violation, the pMSSM. The 19 soft SUSY breaking parameters are chosen so to satisfy all existing experimental and theoretical constraints assuming that the WIMP is the lightest neutralino. We scan this parameter space twice using both flat and log priors and compare the results which yield similar conclusions. Constraints from both LEP and the Tevatron play an important role in obtaining our final model samples. Implications for future TeV-scale e{sup +}e{sup -} linear colliders (LC) are discussed.

None

Head-tail modes are described here for the space charge tune shift significantly exceeding the synchrotron tune. General equation for the modes is derived. Spatial shapes of the modes, their frequencies, and coherent growth rates are explored. The Landau damping rates are also found. Suppression of the transverse mode coupling instability by the space charge is explained.

Recent measurements of the branching fraction for D{sub s} {yields} {ell}v disagree with the Standard Model expectation, which relies on calculations of f{sub D{sub s}} from lattice QCD. This paper uses recent preliminary measurements from CLEO and a new preliminary lattice-QCD result from this conference to update the significance of the discrepancy. The 'f{sub D{sub s}} puzzle' stands now at 3.5{sigma}, with {sigma} predominantly from the statistical uncertainty of the experiments. New physics scenarios that could solve the puzzle would also lead to non-Standard amplitudes mediating the semileptonic decays D {yields} K{ell}v. This paper shows where the new amplitudes enter the differential rate and outlines where lattice QCD calculations are needed to confront recent and forthcoming measurements.

The Department of Energy Savannah River Site vitrifies nuclear waste incident to defense programs through its Defense Waste Processing Facility (DWPF). The piping in the DWPF seal pot jumper configuration must withstand the stresses during an unlikely but potential deflagration event, and maintain its safety function for a 20-year service life. Carbon fiber-reinforced epoxy composites (CFR) were proposed for protection and reinforcement of piping during such an event. The proposed CFR materials have been ASME-approved (Section XI, Code Case N-589-1) for post-construction maintenance and is DOT-compliant per 49CFR 192 and 195. The proposed carbon fiber/epoxy composite reinforcement system was originally developed for pipeline rehabilitation and post-construction maintenance in petrochemical, refineries, DOT applications and other industries. The effects of ionizing radiation on polymers and organic materials have been studied for many years. The majority of available data are based on traditional exposures to gamma irradiation at high dose rates ({approx}10,000 Gy/hr) allowing high total dose within reasonable test periods and general comparison of different materials exposed at such conditions. However, studies in recent years have shown that degradation of many polymers are sensitive to dose rate, with more severe degradation often observed at similar or even lower total doses when exposed …

Chemical modification of nanocrystal surfaces is fundamentally important to their assembly, their implementation in biology and medicine, and greatly impacts their electrical and optical properties. However, it remains a major challenge owing to a lack of analytical tools to directly determine nanoparticle surface structure. Early nuclear magnetic resonance (NMR) and X-ray photoelectron spectroscopy (XPS) studies of CdSe nanocrystals prepared in tri-n-octylphosphine oxide (1) and tri-n-octylphosphine (2), suggested these coordinating solvents are datively bound to the particle surface. However, assigning the broad NMR resonances of surface-bound ligands is complicated by significant concentrations of phosphorus-containing impurities in commercial sources of 1, and XPS provides only limited information about the nature of the phosphorus containing molecules in the sample. More recent reports have shown the surface ligands of CdSe nanocrystals prepared in technical grade 1, and in the presence of alkylphosphonic acids, include phosphonic and phosphinic acids. These studies do not, however, distinguish whether these ligands are bound datively, as neutral, L-type ligands, or by X-type interaction of an anionic phosphonate/phosphinate moiety with a surface Cd{sup 2+} ion. Answering this question would help clarify why ligand exchange with such particles does not proceed generally as expected based on a L-type ligand model. By …

We study the QCD evolution for the twist-three quark-gluon correlation functions associated with the transverse momentum odd quark distributions. Different from that for the leading twist quark distributions, these evolution equations involve more general twist-three functions beyond the correlation functions themselves. They provide important information on nucleon structure, and can be studied in the semi-inclusive hadron production in deep inelastic scattering and Drell-Yan lepton pair production in pp scattering process.

We derive an index theorem for the Dirac operator in the background of various topological excitations on an R{sup 3} x S{sup 1} geometry. The index theorem provides more refined data than the APS index for an instanton on R{sup 4} and reproduces it in decompactification limit. In the R{sup 3} limit, it reduces to the Callias index theorem. The index is expressed in terms of topological charge and the {eta}-invariant associated with the boundary Dirac operator. Neither topological charge nor {eta}-invariant is typically an integer, however, the non-integer parts cancel to give an integer-valued index. Our derivation is based on axial current non-conservation--an exact operator identity valid on any four-manifold--and on the existence of a center symmetric, or approximately center symmetric, boundary holonomy (Wilson line). We expect the index theorem to usefully apply to many physical systems of interest, such as low temperature (large S{sup 1}, confined) phases of gauge theories, center stabilized Yang-Mills theories with vector-like or chiral matter (at S{sup 1} of any size), and supersymmetric gauge theories with supersymmetry-preserving boundary conditions (also at any S{sup 1}). In QCD-like and chiral gauge theories, the index theorem should shed light into the nature of topological excitations responsible for …

We classify all bifurcations from traveling waves to non-trivial time-periodic solutions of the Benjamin-Ono equation that are predicted by linearization. We use a spectrally accurate numerical continuation method to study several paths of non-trivial solutions beyond the realm of linear theory. These paths are found to either re-connect with a different traveling wave or to blow up. In the latter case, as the bifurcation parameter approaches a critical value, the amplitude of the initial condition grows without bound and the period approaches zero. We propose a conjecture that gives the mapping from one bifurcation to its counterpart on the other side of the path of non-trivial solutions. By experimentation with data fitting, we identify the form of the exact solutions on the path connecting two traveling waves, which represents the Fourier coefficients of the solution as power sums of a finite number of particle positions whose elementary symmetric functions execute simple orbits in the complex plane (circles or epicycles). We then solve a system of algebraic equations to express the unknown constants in the new representation in terms of the mean, a spatial phase, a temporal phase, four integers (enumerating the bifurcation at each end of the path) and one …

The local coupling of two photons to the fundamental quark currents of a hadron gives an energy-independent contribution to the Compton amplitude proportional to the charge squared of the struck quark, a contribution which has no analog in hadron scattering reactions. We show that this local contribution has a real phase and is universal, giving the same contribution for real or virtual Compton scattering for any photon virtuality and skewness at fixed momentum transfer squared t. The t-dependence of this J = 0 fixed Regge pole is parameterized by a yet unmeasured even charge-conjugation form factor of the target nucleon. The t = 0 limit gives an important constraint on the dependence of the nucleon mass on the quark mass through the Weisberger relation. We discuss how this 1=x form factor can be extracted from high energy deeply virtual Compton scattering and examine predictions given by models of the H generalized parton distribution.

We study the evolutionary Prisoner's Dilemma on two social networks substrates obtained from actual relational data. We find very different cooperation levels on each of them that cannot be easily understood in terms of global statistical properties of both networks. We claim that the result can be understood at the mesoscopic scale, by studying the community structure of the networks. We explain the dependence of the cooperation level on the temptation parameter in terms of the internal structure of the communities and their interconnections. We then test our results on community-structured, specifically designed artificial networks, finding a good agreement with the observations in both real substrates. Our results support the conclusion that studies of evolutionary games on model networks and their interpretation in terms of global properties may not be sufficient to study specific, real social systems. Further, the study allows us to define new quantitative parameters that summarize the mesoscopic structure of any network. In addition, the community perspective may be helpful to interpret the origin and behavior of existing networks as well as to design structures that show resilient cooperative behavior.

Unstable austenitic stainless steels undergo a strain-induced martensite transformation. The effect of hydrogen on this transformation is not well understood. Some researchers believe that hydrogen makes the transformation to martensite more difficult because hydrogen is an austenite stabilizer. Others believe that hydrogen has little or no effect at all on the transformation and claim that the transformation is simply a function of strain and temperature. Still other researchers believe that hydrogen should increase the ability of the metal to transform due to hydrogen-enhanced dislocation mobility and slip planarity. While the role of hydrogen on the martensite transformation is still debated, it has been experimentally verified that this transformation does occur in hydrogen-charged materials. What is the effect of strain-induced martensite on hydrogen embrittlement? Martensite near crack-tips or other highly strained regions could provide much higher hydrogen diffusivity and allow for quicker hydrogen concentration. Martensite may be more intrinsically brittle than austenite and has been shown to be severely embrittled by hydrogen. However, it does not appear to be a necessary condition for embrittlement since Type 21-6-9 stainless steel is more stable than Type 304L stainless steel but susceptible to hydrogen embrittlement. In this study, the effect of hydrogen on strain-induced …

The Silicon Detector proposed for the International Linear Collider (ILC) requires electronic read-out that can be tightly coupled to the silicon detectors envisioned for the tracker and the electromagnetic calorimeter. The KPiX is a 1024-channel read-out chip that bump-bonds to the detector and communicates through a few digital signals, power, and detector bias. The KPiX front-end is a low-noise dual-range charge-amplifier with a dynamic range of 17 bit, achieved by autonomous switching of the feedback capacitor. The device takes advantage of the ILC duty cycle of 1 ms trains at 5 Hz rate by lowering the supply current after the data acquisition cycle for an average power consumption of <20 {micro}W/channel. During the 1 ms train, up to four events exceeding a programmable threshold can be stored, the amplitude as a voltage on a capacitor for subsequent digitization, the event time in digital format. The chip can be configured for other than ILC applications.

The chemistry of Br species associated with sea salt ice and aerosols has been implicated in the episodes of ozone depletion reported at Arctic sunrise. However, Br{sup -} is only a minor component in sea salt, which has a Br{sup -}/Cl{sup -} molar ratio of {approx}0.0015. Sea salt is a complex mixture of many different species, with NaCl as the primary component. In recent years experimental and theoretical studies have reported enhancement of the large, more polarizable halide ion at the liquid/vapor interface of corresponding aqueous alkali halide solutions. The proposed enhancement is likely to influence the availability of sea salt Br{sup -} for heterogeneous reactions such as those involved in the ozone depletion episodes. We report here ambient pressure x-ray photoelectron spectroscopy studies and molecular dynamics simulations showing direct evidence of Br{sup -} enhancement at the interface of an aqueous NaCl solution doped with bromide. The experiments were carried out on samples with Br{sup -}/Cl{sup -} ratios in the range 0.1% to 10%, the latter being also the ratio for which simulations were carried out. This is the first direct measurement of interfacial enhancement of Br{sup -} in a multi-component solution with particular relevance to sea salt chemistry.

The correspondence between theories in anti-de Sitter space and conformal field theories in physical space-time leads to an analytic, semiclassical model for strongly-coupled QCD. Light-front holography allows hadronic amplitudes in the AdS fifth dimension to be mapped to frame-independent light-front wavefunctions of hadrons in physical space-time, thus providing a relativistic description of hadrons at the amplitude level. We identify the AdS coordinate z with an invariant light-front coordinate {zeta} which separates the dynamics of quark and gluon binding from the kinematics of constituent spin and internal orbital angular momentum. The result is a single-variable light-front Schroedinger equation for QCD which determines the eigenspectrum and the light-front wavefunctions of hadrons for general spin and orbital angular momentum. The mapping of electromagnetic and gravitational form factors in AdS space to their corresponding expressions in light-front theory confirms this correspondence. Some novel features of QCD are discussed, including the consequences of confinement for quark and gluon condensates and the behavior of the QCD coupling in the infrared. The distinction between static structure functions such as the probability distributions computed from the square of the light-front wavefunctions versus dynamical structure functions which include the effects of rescattering is emphasized. A new method for computing …

Time-resolved soft X-ray transmission microscopy is applied to image the current-induced resonant dynamics of the magnetic vortex core realized in a micron-sized Permalloy disk. The high spatial resolution better than 25 nm enables us to observe the resonant motion of the vortex core. The result also provides the spin polarization of the current to be 0.67 {+-} 0.16 for Permalloy by fitting the experimental results with an analytical model in the framework of the spin-transfer torque.

Charged lepton-flavor violating processes are extremely rare in the Standard Model, but they are predicted to occur in several beyond-the-Standard Model theories, including Supersymmetry or models with leptoquarks or compositeness. We present a search for such processes in a sample of 117 x 10{sup 6} {Upsilon}(3S) decays recorded with the BABAR detector. We place upper limits on the branching fractions BF({Upsilon}(3S) {yields} e{sup {+-}}{tau}{sup {-+}}) < 5.0 x 10{sup -6} and BF({Upsilon}(3S) {yields} {mu}{sup {+-}}{tau}{sup {-+}}) < 4.1 x 10{sup -6} at 90% confidence level. These results are used to place lower limits on the mass scale of beyond-the-Standard Model physics contributing to lepton-flavor violating decays of the {Upsilon}(3S).