
The field of particle therapy is quickly growing and yet it's more widespread
adoption is limited by size, cost and adaptation to the more conformal
treatment techniques. In order to realize the benefits of this modality the
equipment used to generate and deliver the beam is evolving. The accelerator is
one of the key components and its future is dictated by the ability to
accommodate the clinical requirements. This lecture is intended to provide an
introduction to these requirements and identify how synchrotrons are designed
to deliver the desired beams as well as what limitations exist and expectations
for the future of synchrotrons.

Parametric images provide insight into the spatial distribution of
physiological parameters, but they are often extremely noisy, due to low SNR of
tomographic data. Direct estimation from projections allows accurate noise
modeling, improving the results of postreconstruction fitting. We propose a
method, which we name kinetic compressive sensing (KCS), based on a
hierarchical Bayesian model and on a novel reconstruction algorithm, that
encodes sparsity of kinetic parameters. Parametric maps are reconstructed by
maximizing the joint probability, with an Iterated Conditional Modes (ICM)
approach, alternating the optimization of activity time series (OSMAPOSL),
and kinetic parameters (MAPLM). We evaluated the proposed algorithm on a
simulated dynamic phantom: a bias/variance study confirmed how direct estimates
can improve the quality of parametric maps over a postreconstruction fitting,
and showed how the novel sparsity prior can further reduce their variance,
without affecting bias. Real FDG PET human brain data (Siemens mMR, 40min)
images were also processed. Results enforced how the proposed KCSregularized
direct method can produce spatially coherent images and parametric maps, with
lower spatial noise and better tissue contrast. A GPUbased open source
implementation of the algorithm is provided.

This article is concerned with the dynamics of glacial cycles observed in the
geological record of the Pleistocene Epoch. It focuses on a conceptual model
proposed by Maasch and Saltzman [J. Geophys. Res.,95, D2 (1990), pp.
19551963], which is based on physical arguments and emphasizes the role of
atmospheric CO2 in the generation and persistence of periodic orbits (limit
cycles). The model consists of three ordinary differential equations with four
parameters for the anomalies of the total global ice mass, the atmospheric CO2
concentration, and the volume of the North Atlantic Deep Water (NADW). In this
article, it is shown that a simplified twodimensional symmetric version
displays many of the essential features of the full model, including
equilibrium states, limit cycles, their basic bifurcations, and a
BogdanovTakens point that serves as an organizing center for the local and
global dynamics. Also, symmetry breaking splits the BogdanovTakens point into
two, with different local dynamics in their neighborhoods.

Quantitative extraction of highdimensional mineable data from medical images
is a process known as radiomics. Radiomics is foreseen as an essential
prognostic tool for cancer risk assessment and the quantification of
intratumoural heterogeneity. In this work, 1615 radiomic features (quantifying
tumour image intensity, shape, texture) extracted from pretreatment FDGPET
and CT images of 300 patients from four different cohorts were analyzed for the
risk assessment of locoregional recurrences (LR) and distant metastases (DM) in
headandneck cancer. Prediction models combining radiomic and clinical
variables were constructed via random forests and imbalanceadjustment
strategies using two of the four cohorts. Independent validation of the
prediction and prognostic performance of the models was carried out on the
other two cohorts (LR: AUC = 0.69 and CI = 0.67; DM: AUC = 0.86 and CI = 0.88).
Furthermore, the results obtained via KaplanMeier analysis demonstrated the
potential of radiomics for assessing the risk of specific tumour outcomes using
multiple stratification groups. This could have important clinical impact,
notably by allowing for a better personalization of chemoradiation treatments
for headandneck cancer patients from different risk groups.

Plasma velocity and magnetic field measurements from the Voyager 2 mission
are used to study solar wind turbulence in the slow solar wind at two different
heliocentric distances, 5 and 29 astronomical units, sufficiently far apart to
provide information on the radial evolution of this turbulence. The magnetic
helicity and the crosshelicity, which express the correlation between the
plasma velocity and the magnetic field, are used to characterize the
turbulence. Wave number spectra are computed by means of the Taylor hypothesis
applied to time resolved single point Voyager 2 measurements. The overall
picture we get is complex and difficult to interpret. A substantial decrease of
the crosshelicity at smaller scales (over 13 hours of observation) with
increasing heliocentric distance is observed. At 5 AU the only peak in the
probability density of the normalized residual energy is negative, near 0.5.
At 29 AU the probability density becomes doubly peaked, with a negative peak at
0.5 and a smaller peak at a positive values of about 0.7. A decrease of the
crosshelicity for increasing heliocentric distance is observed, together with
a reduction of the unbalance toward the magnetic energy of the energy of the
fluctuations. For the smaller scales, we found that at 29 AU the normalized
polarization is small and positive on average (about 0.1), it is instead zero
at 5 AU. For the larger scales, the polarization is low and positive at 5 AU
(average around 0.1) while it is negative (around  0.15) at 29 AU.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed
for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV)
gammarays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift,
GASPWEBT, and several optical and radio telescopes. We model the dayscale
SEDs with onezone and twozone synchrotron selfCompton (SSC) models,
investigate the physical parameters, and evaluate whether the observed
broadband SED variability can be associated to variations in the relativistic
particle population.
Flux variability was remarkable in the Xray and VHE bands while it was minor
or not significant in the other bands. The onezone SSC model can describe
reasonably well the SED of each day for the 13 consecutive days. This flaring
activity is also very well described by a twozone SSC model, where one zone is
responsible for the quiescent emission while the other smaller zone, which is
spatially separated from the first one, contributes to the dailyvariable
emission occurring in Xrays and VHE gammarays.
Both the onezone SSC and the twozone SSC models can describe the daily SEDs
via the variation of only four or five model parameters, under the hypothesis
that the variability is associated mostly to the underlying particle
population. This shows that the particle acceleration and cooling mechanism
producing the radiating particles could be the main one responsible for the
broadband SED variations during the flaring episodes in blazars. The twozone
SSC model provides a better agreement to the observed SED at the narrow peaks
of the low and highenergy bumps during the highest activity, although the
reported onezone SSC model could be further improved by the variation of the
parameters related to the emitting region itself ($\delta$, $B$ and $R$), in
addition to the parameters related to the particle population.

SuperKamiokande (SK) can search for weakly interacting massive particles
(WIMPs) by detecting neutrinos produced from WIMP annihilations occurring
inside the Sun. In this analysis, we include neutrino events with interaction
vertices in the detector in addition to upwardgoing muons produced in the
surrounding rock. Compared to the previous result, which used the upwardgoing
muons only, the signal acceptances for light (fewGeV/$c^2$ $\sim$
200GeV/$c^2$) WIMPs are significantly increased. We fit 3903 days of SK data
to search for the contribution of neutrinos from WIMP annihilation in the Sun.
We found no significant excess over expected atmosphericneutrino background
and the result is interpreted in terms of upper limits on WIMPnucleon elastic
scattering cross sections under different assumptions about the annihilation
channel. We set the current best limits on the spindependent (SD) WIMPproton
cross section for WIMP masses below 200 GeV/$c^2$ (at 10 GeV/$c^2$, 1.49$\times
10^{39}$ cm$^2$ for $\chi\chi\rightarrow b \bar{b}$ and 1.31$\times 10^{40}$
cm$^2$ for $\chi\chi\rightarrow\tau^+\tau^$ annihilation channels), also
ruling out some fraction of WIMP candidates with spinindependent (SI) coupling
in the fewGeV/$c^2$ mass range.

Snowmass is a US longterm planning study for the highenergy community by
the American Physical Society's Division of Particles and Fields. For its
simulation studies, opportunistic resources are harnessed using the Open
Science Grid infrastructure. Late binding grid technology, GlideinWMS, was used
for distributed scheduling of the simulation jobs across many sites mainly in
the US. The pilot infrastructure also uses the Parrot mechanism to dynamically
access CvmFS in order to ascertain a homogeneous environment across the nodes.
This report presents the resource usage and the storage model used for
simulating large statistics Standard Model backgrounds needed for Snowmass
Energy Frontier studies.

The OPERA collaboration has claimed evidence of superluminal {\nu}{_\mu}
propagation between CERN and the LNGS. Cohen and Glashow argued that such
neutrinos should lose energy by producing photons and e+e pairs, through Z0
mediated processes analogous to Cherenkov radiation. In terms of the parameter
delta=(v^2_nuv^2_c)/v^2_c, the OPERA result implies delta = 5 x 10^5. For
this value of \delta a very significant deformation of the neutrino energy
spectrum and an abundant production of photons and e+e pairs should be
observed at LNGS. We present an analysis based on the 2010 and part of the 2011
data sets from the ICARUS experiment, located at Gran Sasso National Laboratory
and using the same neutrino beam from CERN. We find that the rates and
deposited energy distributions of neutrino events in ICARUS agree with the
expectations for an unperturbed spectrum of the CERN neutrino beam. Our results
therefore refute a superluminal interpretation of the OPERA result according to
the Cohen and Glashow prediction for a weak current analog to Cherenkov
radiation. In particular no superluminal Cherenkov like e+e pair or gamma
emission event has been directly observed inside the fiducial volume of the
"bubble chamber like" ICARUS TPCLAr detector, setting the much stricter limit
of delta < 2.5 10^8 at the 90% confidence level, comparable with the one due
to the observations from the SN1987A.

We present BVI CCD photometry of 10 northern open clusters, Berkeley 43,
Berkeley 45, Berkeley 47, NGC 6846, Berkeley 49, Berkeley 51, Berkeley 89,
Berkeley 91, Tombaugh 4 and Berkeley 9, and estimate their fundamental
parameters. Eight of the clusters are located in the first galactic quadrant
and 2 are in the second. This is the first optical photometry for 8 clusters.
All of them are embedded in rich galactic fields and have large reddening
towards them (E(BV) = 1.0  2.3 mag). There is a possibility that some of
these difficulttostudy clusters may be asterisms rather than physical
systems, but assuming they are physical clusters, we find that 8 of them are
located beyond 2 kpc, and 6 clusters (60% of the sample) are located well above
or below the Galactic plane. Seven clusters have ages 500 Myr or less and the
other 3 are 1 Gyr or more in age. This sample of clusters has increased the
optical photometry of clusters in the second half of the first galactic
quadrant, beyond 2 kpc, from 10 to 15. NGC 6846 is found to be one of the most
distant clusters in this region of the Galaxy.

We describe a molecular dynamics framework for the direct calculation of the
shortranged structural forces underlying grainboundary premelting and
graincoalescence in solidification. The method is applied in a comparative
study of (i) a Sigma 9 <115> 120 degress twist and (ii) a Sigma 9 <110> {411}
symmetric tilt boundary in a classical embeddedatom model of elemental Ni.
Although both boundaries feature highly disordered structures near the melting
point, the nature of the temperature dependence of the width of the disordered
regions in these boundaries is qualitatively different. The former boundary
displays behavior consistent with a logarithmically diverging premelted layer
thickness as the melting temperature is approached from below, while the latter
displays behavior featuring a finite grainboundary width at the melting point.
It is demonstrated that both types of behavior can be quantitatively described
within a sharpinterface thermodynamic formalism involving a widthdependent
interfacial free energy, referred to as the disjoining potential. The
disjoining potential for boundary (i) is calculated to display a monotonic
exponential dependence on width, while that of boundary (ii) features a weak
attractive minimum. The results of this work are discussed in relation to
recent simulation and theoretical studies of the thermodynamic forces
underlying grainboundary premelting.

We describe a molecular dynamics framework for the direct calculation of the
shortranged structural forces underlying grainboundary premelting and
graincoalescence in solidification. The method is applied in a comparative
study of (i) a Sigma 9 <115> 120 degress twist and (ii) a Sigma 9 <110> {411}
symmetric tilt boundary in a classical embeddedatom model of elemental Ni.
Although both boundaries feature highly disordered structures near the melting
point, the nature of the temperature dependence of the width of the disordered
regions in these boundaries is qualitatively different. The former boundary
displays behavior consistent with a logarithmically diverging premelted layer
thickness as the melting temperature is approached from below, while the latter
displays behavior featuring a finite grainboundary width at the melting point.
It is demonstrated that both types of behavior can be quantitatively described
within a sharpinterface thermodynamic formalism involving a widthdependent
interfacial free energy, referred to as the disjoining potential. The
disjoining potential for boundary (i) is calculated to display a monotonic
exponential dependence on width, while that of boundary (ii) features a weak
attractive minimum. The results of this work are discussed in relation to
recent simulation and theoretical studies of the thermodynamic forces
underlying grainboundary premelting.

Intrinsic optical transition energies for isolated and individual single wall
carbon nanotubes grown over trenches are measured using tunable resonant Raman
scattering. Previously measured E22_S optical transitions from nanotubes in
surfactants are blue shifted 7090 meV with respect to our measurements of
nanotubes in air. This large shift in the exciton energy is attributed to a
larger change of the exciton binding energy than the bandgap renormalization
as the surrounding dielectric constant increases.

Intrinsic optical transition energies for isolated and individual single wall
carbon nanotubes grown over trenches are measured using tunable resonant Raman
scattering. Previously measured E22_S optical transitions from nanotubes in
surfactants are blue shifted 7090 meV with respect to our measurements of
nanotubes in air. This large shift in the exciton energy is attributed to a
larger change of the exciton binding energy than the bandgap renormalization
as the surrounding dielectric constant increases.

Processes of mass communications in complicated social or sociobiological
systems such as marketing, economics, politics, animal populations, etc. as a
subject for the special scientific discipline  "mediaphysics"  are considered
in its relation with sociophysics. A new statistical physics approach to
analyze these phenomena is proposed. A keystone of the approach is an analysis
of population distribution between two or many alternatives: brands, political
affiliations, or opinions. Relative distances between a state of a "person's
mind" and the alternatives are measures of propensity to buy (to affiliate, or
to have a certain opinion). The distribution of population by those relative
distances is time dependent and affected by external (economic, social,
marketing, natural) and internal (meanfield influential propagation of
opinions, synergy effects, etc.) factors, considered as fields. Specifically,
the interaction and opinioninfluence field can be generalized to incorporate
important elements of Isingspin based sociophysical models and
kineticequation ones. The distributions were described by a Schrodingertype
equation in terms of Green's functions. The developed approach has been applied
to a real massmedia efficiency problem for a large company and generally
demonstrated very good results despite low initial correlations of factors and
the target variable.

To achieve quantitative interpretation of Piezoresponse Force Microscopy
(PFM), including resolution limits, tip bias and straininduced phenomena and
spectroscopy, analytical representations for tipinduced electroelastic fields
inside the material are derived for the cases of weak and strong indentation.
In the weak indentation case, electrostatic field distribution is calculated
using image charge model. In the strong indentation case, the solution of the
coupled electroelastic problem for piezoelectric indentation is used to obtain
the electric field and strain distribution in the ferroelectric material. This
establishes a complete continuum mechanics description of the PFM contact
mechanics and imaging mechanism. The electroelastic field distribution allows
signal generation volume in PFM to be determined. These rigorous solutions are
compared with the electrostatic point charge and sphereplane models, and the
applicability limits for asymptotic point charge and point force models are
established. The implications of these results for ferroelectric polarization
switching processes are analyzed.

The quantum statistical dynamics of a position coordinate x coupled to a
reservoir requires theoretically two copies of the position coordinate within
the reduced density matrix description. One coordinate moves forward in time
while the other coordinate moves backward in time. It is shown that quantum
dissipation induces, in the plane of the forward and backward motions, a
noncommutative geometry. The noncommutative geometric plane is a consequence of
a quantum dissipation induced phase interference which is closely analogous to
the AharanovBohm effect.

Precision measurements of the microwave surface resistance $R_{s}$ of in situ
$MgB_{2}$ films directly reveal an exponential behavior of $R_{s}$ at low
temperature indicating a fullygapped order parameter. The entire temperature
dependence of $R_{s}$ is well described by a MattisBardeen formalism but with
a small gap ratio of $\Delta (0)/kT_{c}=0.72$, corresponding to $\Delta
(0)=1.9meV$.

The Great Observatories Origins Deep Survey (GOODS) is designed to gather the
best and deepest multiwavelength data for studying the formation and evolution
of galaxies and active galactic nuclei, the distribution of dark and luminous
matter at high redshift, the cosmological parameters from distant supernovae,
and the extragalactic background light. The program uses the most powerful
space and groundbased telescopes to cover two fields, each 10'x16', centered
on the Hubble Deep Field North and the Chandra Deep Field South, already the
sites of extensive observations from Xray through radio wavelengths. GOODS
incorporates 3.624 micron observations from a SIRTF Legacy Program, fourband
ACS imaging from an HST Treasury Program, and extensive new groundbased
imaging and spectroscopy. GOODS data products will be made available on a rapid
timescale, enabling community research on a wide variety of topics. Here we
describe the project, emphasizing its application for studying the mass
assembly history of galaxies.

We present the first reliable calculation of the collective mode structure of
a strongly coupled electronic bilayer. The calculation is based on a classical
model through the $3^{rd}$ frequencymomentsumrule preserving Quasi Localized
Charge Approximation, using the recently calculated Hypernetted Chain pair
correlation functions. The spectrum shows an energy gap at $k=0$ and the
absence of a previously conjectured dynamical instability.

We present theoretical arguments and simulation data indicating that the
scaling of earthquake events in models of faults with longrange stress
transfer is composed of at least three distinct regions. These regions
correspond to three classes of earthquakes with different underlying physical
mechanisms. In addition to the events that exhibit scaling, there are larger
``breakout'' events that are not on the scaling plot. We discuss the
interpretation of these events as fluctuations in the vicinity of a spinodal
critical point.

We study the distribution of fluctuations over a time scale $\Delta t$ (i.e.,
the returns) of the S&P 500 index by analyzing three distinct databases.
Database (i) contains approximately 1 million records sampled at 1 min
intervals for the 13year period 19841996, database (ii) contains 8686 daily
records for the 35year period 19621996, and database (iii) contains 852
monthly records for the 71year period 19261996. We compute the probability
distributions of returns over a time scale $\Delta t$, where $\Delta t$ varies
approximately over a factor of 10^4  from 1 min up to more than 1 month. We
find that the distributions for $\Delta t \leq$ 4 days (1560 mins) are
consistent with a powerlaw asymptotic behavior, characterized by an exponent
$\alpha \approx 3$, well outside the stable L\'evy regime $0 < \alpha < 2$. To
test the robustness of the S&P result, we perform a parallel analysis on two
other financial market indices. Database (iv) contains 3560 daily records of
the NIKKEI index for the 14year period 198497, and database (v) contains 4649
daily records of the HangSeng index for the 18year period 198097. We find
estimates of $\alpha$ consistent with those describing the distribution of S&P
500 dailyreturns. One possible reason for the scaling of these distributions
is the long persistence of the autocorrelation function of the volatility. For
time scales longer than $(\Delta t)_{\times} \approx 4$ days, our results are
consistent with slow convergence to Gaussian behavior.

Background: Many attempts have been made to resolve in time the folding of
model proteins in computer simulations. Different computational approaches have
emerged. Some of these approaches suffer from the insensitivity to the
geometrical properties of the proteins (lattice models), while others are
computationally heavy (traditional MD).
Results: We use a recentlyproposed approach of Zhou and Karplus to study the
folding of the protein model based on the discrete time molecular dynamics
algorithm. We show that this algorithm resolves with respect to time the
folding  unfolding transition. In addition, we demonstrate the ability to
study the coreof the model protein.
Conclusion: The algorithm along with the model of interresidue interactions
can serve as a tool to study the thermodynamics and kinetics of protein models.

Molecular dynamics simulations of folding in an offlattice protein model
reveal a nucleation scenario, in which a few welldefined contacts are formed
with high probability in the transition state ensemble of conformations. Their
appearance determines folding cooperativity and drives the model protein into
its folded conformation.

Maximally complex Regge trajectories are introduced for which both Re
$\alpha(s)$ and Im $\alpha(s)$ grow as $s^{1\epsilon}$ ($\epsilon$ small and
positive). Our expression reduces to the standard real linear form as the
imaginary part (proportional to $\epsilon$) goes to zero. A scaling formula for
the total widths emerges: $\Gamma_{TOT}/M\to$ constant for large M, in very
good agreement with data for mesons and baryons. The unitarity corrections also
enhance the spacelike slopes from their timelike values, thereby resolving an
old problem with the $\rho$ trajectory in $\pi N$ charge exchange. Finally, the
unitarily enhanced intercept, $\alpha_{\rho}\approx 0.525$, \nolinebreak is in
good accord with the DonnachieLandshoff total cross section analysis.