Optimal Axes of Siberian Snakes for Polarized Proton Acceleration

Accelerating polarized proton beams and storing them for many turns can lead to a loss of polarization when accelerating through energies where a spin rotation frequency is in resonance with orbit oscillation frequencies. First-order resonance effects can be avoided by installing Siberian Snakes in the ring, devices which rotate the spin by 180 degrees around the snake axis while not changing the beam's orbit significantly. For large rings, several Siberian Snakes are required. Here a criterion will be derived that allows to find an optimal choice of the snake axes. Rings with super-period four are analyzed in detail, and the HERA proton ring is used as an example for approximate four-fold symmetry. The proposed arrangement of Siberian Snakes matches their effects so that all spin-orbit coupling integrals vanish at all energies and therefore there is no first-order spin-orbit coupling at all for this choice, which I call snakes matching. It will be shown that in general at least eight Siberian Snakes are needed and that there are exactly four possibilities to arrange their axes. When the betatron phase advance between snakes is chosen suitably, four Siberian Snakes can be sufficient. To show that favorable choice of snakes have been found, polarized protons are tracked for part of HERA-p's acceleration cycle which shows that polarization is preserved best for the here proposed arrangement of Siberian Snakes.Comment: 14 pages, 16 figure

Quasiperiodic spin-orbit motion and spin tunes in storage rings

We present an in-depth analysis of the concept of spin precession frequency for integrable orbital motion in storage rings. Spin motion on the periodic closed orbit of a storage ring can be analyzed in terms of the Floquet theorem for equations of motion with periodic parameters and a spin precession frequency emerges in a Floquet exponent as an additional frequency of the system. To define a spin precession frequency on nonperiodic synchro-betatron orbits we exploit the important concept of quasiperiodicity. This allows a generalization of the Floquet theorem so that a spin precession frequency can be defined in this case too. This frequency appears in a Floquet-like exponent as an additional frequency in the system in analogy with the case of motion on the closed orbit. These circumstances lead naturally to the definition of the uniform precession rate and a definition of spin tune. A spin tune is a uniform precession rate obtained when certain conditions are fulfilled. Having defined spin tune we define spin-orbit resonance on synchro--betatron orbits and examine its consequences. We give conditions for the existence of uniform precession rates and spin tunes (e.g. where small divisors are controlled by applying a Diophantine condition) and illustrate the various aspects of our description with several examples. The formalism also suggests the use of spectral analysis to ``measure'' spin tune during computer simulations of spin motion on synchro-betatron orbits.Comment: 62 pages, 1 figure. A slight extension of the published versio

Two-component radiation model of the sonoluminescing bubble

Based on the experimental data from Weninger, Putterman & Barber, Phys. Rev. (E), 54, R2205 (1996), we offer an alternative interpretation of their experimetal results. A model of sonoluminescing bubble which proposes that the electromagnetic radiation originates from two sources: the isotropic black body or bramsstrahlung emitting core and dipole radiation-emitting shell of accelerated electrons driven by the liquid-bubble interface is outlined.Comment: 5 pages Revtex, submitted to Phys. Rev.

Developing and validating a predictive model for stroke progression

<p><b>Background:</b> Progression is believed to be a common and important complication in acute stroke, and has been associated with increased mortality and morbidity. Reliable identification of predictors of early neurological deterioration could potentially benefit routine clinical care. The aim of this study was to identify predictors of early stroke progression using two independent patient cohorts.</p> <p><b>Methods:</b> Two patient cohorts were used for this study – the first cohort formed the training data set, which included consecutive patients admitted to an urban teaching hospital between 2000 and 2002, and the second cohort formed the test data set, which included patients admitted to the same hospital between 2003 and 2004. A standard definition of stroke progression was used. The first cohort (n = 863) was used to develop the model. Variables that were statistically significant (p < 0.1) on univariate analysis were included in the multivariate model. Logistic regression was the technique employed using backward stepwise regression to drop the least significant variables (p > 0.1) in turn. The second cohort (n = 216) was used to test the performance of the model. The performance of the predictive model was assessed in terms of both calibration and discrimination. Multiple imputation methods were used for dealing with the missing values.</p> <p><b>Results:</b> Variables shown to be significant predictors of stroke progression were conscious level, history of coronary heart disease, presence of hyperosmolarity, CT lesion, living alone on admission, Oxfordshire Community Stroke Project classification, presence of pyrexia and smoking status. The model appears to have reasonable discriminative properties [the median receiver-operating characteristic curve value was 0.72 (range 0.72–0.73)] and to fit well with the observed data, which is indicated by the high goodness-of-fit p value [the median p value from the Hosmer-Lemeshow test was 0.90 (range 0.50–0.92)].</p> <p><b>Conclusion:</b> The predictive model developed in this study contains variables that can be easily collected in practice therefore increasing its usability in clinical practice. Using this analysis approach, the discrimination and calibration of the predictive model appear sufficiently high to provide accurate predictions. This study also offers some discussion around the validation of predictive models for wider use in clinical practice.</p&gt

Accurate and efficient spin integration for particle accelerators

Accurate spin tracking is a valuable tool for understanding spin dynamics in particle accelerators and can help improve the performance of an accelerator. In this paper, we present a detailed discussion of the integrators in the spin tracking code gpuSpinTrack. We have implemented orbital integrators based on drift-kick, bend-kick, and matrix-kick splits. On top of the orbital integrators, we have implemented various integrators for the spin motion. These integrators use quaternions and Romberg quadratures to accelerate both the computation and the convergence of spin rotations. We evaluate their performance and accuracy in quantitative detail for individual elements as well as for the entire RHIC lattice. We exploit the inherently data-parallel nature of spin tracking to accelerate our algorithms on graphics processing units.Comment: 43 pages, 17 figure

Bubble Shape Oscillations and the Onset of Sonoluminescence

An air bubble trapped in water by an oscillating acoustic field undergoes either radial or nonspherical pulsations depending on the strength of the forcing pressure. Two different instability mechanisms (the Rayleigh--Taylor instability and parametric instability) cause deviations from sphericity. Distinguishing these mechanisms allows explanation of many features of recent experiments on sonoluminescence, and suggests methods for finding sonoluminescence in different parameter regimes.Comment: Phys. Rev. Lett., in pres

Three Dimensional Electrical Impedance Tomography

The electrical resistivity of mammalian tissues varies widely and is correlated with physiological function. Electrical impedance tomography (EIT) can be used to probe such variations in vivo, and offers a non-invasive means of imaging the internal conductivity distribution of the human body. But the computational complexity of EIT has severe practical limitations, and previous work has been restricted to considering image reconstruction as an essentially two-dimensional problem. This simplification can limit significantly the imaging capabilities of EIT, as the electric currents used to determine the conductivity variations will not in general be confined to a two-dimensional plane. A few studies have attempted three-dimensional EIT image reconstruction, but have not yet succeeded in generating images of a quality suitable for clinical applications. Here we report the development of a three-dimensional EIT system with greatly improved imaging capabilities, which combines our 64-electrode data-collection apparatus with customized matrix inversion techniques. Our results demonstrate the practical potential of EIT for clinical applications, such as lung or brain imaging and diagnostic screening

The Sound of Sonoluminescence

We consider an air bubble in water under conditions of single bubble sonoluminescence (SBSL) and evaluate the emitted sound field nonperturbatively for subsonic gas-liquid interface motion. Sound emission being the dominant damping mechanism, we also implement the nonperturbative sound damping in the Rayleigh-Plesset equation for the interface motion. We evaluate numerically the sound pulse emitted during bubble collapse and compare the nonperturbative and perturbative results, showing that the usual perturbative description leads to an overestimate of the maximal surface velocity and maximal sound pressure. The radius vs. time relation for a full SBSL cycle remains deceptively unaffected.Comment: 25 pages; LaTex and 6 attached ps figure files. Accepted for publication in Physical Review

Strength of Higher-Order Spin-Orbit Resonances

When polarized particles are accelerated in a synchrotron, the spin precession can be periodically driven by Fourier components of the electromagnetic fields through which the particles travel. This leads to resonant perturbations when the spin-precession frequency is close to a linear combination of the orbital frequencies. When such resonance conditions are crossed, partial depolarization or spin flip can occur. The amount of polarization that survives after resonance crossing is a function of the resonance strength and the crossing speed. This function is commonly called the Froissart-Stora formula. It is very useful for predicting the amount of polarization after an acceleration cycle of a synchrotron or for computing the required speed of the acceleration cycle to maintain a required amount of polarization. However, the resonance strength could in general only be computed for first-order resonances and for synchrotron sidebands. When Siberian Snakes adjust the spin tune to be 1/2, as is required for high energy accelerators, first-order resonances do not appear and higher-order resonances become dominant. Here we will introduce the strength of a higher-order spin-orbit resonance, and also present an efficient method of computing it. Several tracking examples will show that the so computed resonance strength can indeed be used in the Froissart-Stora formula. HERA-p is used for these examples which demonstrate that our results are very relevant for existing accelerators.Comment: 10 pages, 6 figure