• ### Future (of) Synchrotrons for Particle Therapy(1804.08987)

April 24, 2018 physics.med-ph
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.
• ### Kinetic Compressive Sensing(1803.10045)

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 post-reconstruction 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 (OS-MAP-OSL), and kinetic parameters (MAP-LM). 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 post-reconstruction 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 KCS-regularized direct method can produce spatially coherent images and parametric maps, with lower spatial noise and better tissue contrast. A GPU-based open source implementation of the algorithm is provided.
• ### Modeling the Dynamics of Glacial Cycles(1705.07387)

May 21, 2017 math.DS
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. 1955-1963], 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 two-dimensional symmetric version displays many of the essential features of the full model, including equilibrium states, limit cycles, their basic bifurcations, and a Bogdanov-Takens point that serves as an organizing center for the local and global dynamics. Also, symmetry breaking splits the Bogdanov-Takens point into two, with different local dynamics in their neighborhoods.
• ### Radiomics strategies for risk assessment of tumour failure in head-and-neck cancer(1703.08516)

March 24, 2017 cs.CV
Quantitative extraction of high-dimensional 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 pre-treatment FDG-PET 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 head-and-neck cancer. Prediction models combining radiomic and clinical variables were constructed via random forests and imbalance-adjustment 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 Kaplan-Meier 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 chemo-radiation treatments for head-and-neck cancer patients from different risk groups.
• ### Cross and magnetic helicity in the outer heliosphere from Voyager 2 observations(1504.08154)

Aug. 31, 2015 physics.plasm-ph, astro-ph.SR
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 cross-helicity, 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 cross-helicity at smaller scales (over 1-3 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 cross-helicity 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) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (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 X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone 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 two-zone 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 daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone 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 two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone 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.
• Super-Kamiokande (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 upward-going muons produced in the surrounding rock. Compared to the previous result, which used the upward-going muons only, the signal acceptances for light (few-GeV/$c^2$ $\sim$ 200-GeV/$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 atmospheric-neutrino background and the result is interpreted in terms of upper limits on WIMP-nucleon elastic scattering cross sections under different assumptions about the annihilation channel. We set the current best limits on the spin-dependent (SD) WIMP-proton 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 spin-independent (SI) coupling in the few-GeV/$c^2$ mass range.
• ### Snowmass Energy Frontier Simulations using the Open Science Grid (A Snowmass 2013 whitepaper)(1308.0843)

Oct. 2, 2013 hep-ex, cs.DC
Snowmass is a US long-term planning study for the high-energy 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.
• ### A search for the analogue to Cherenkov radiation by high energy neutrinos at superluminal speeds in ICARUS(1110.3763)

March 8, 2012 hep-ex
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_nu-v^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 TPC-LAr 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.
• ### Optical photometry and basic parameters of 10 unstudied open clusters(1001.3001)

Jan. 18, 2010 astro-ph.GA, astro-ph.SR
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(B-V) = 1.0 - 2.3 mag). There is a possibility that some of these difficult-to-study 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.
• ### Structural disjoining potential for grain boundary premelting and grain coalescence from molecular-dynamics simulations(1001.0583)

Jan. 4, 2010 cond-mat.mtrl-sci
We describe a molecular dynamics framework for the direct calculation of the short-ranged structural forces underlying grain-boundary premelting and grain-coalescence 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 embedded-atom 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 grain-boundary width at the melting point. It is demonstrated that both types of behavior can be quantitatively described within a sharp-interface thermodynamic formalism involving a width-dependent 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 grain-boundary premelting.
• ### Structural disjoining potential for grain boundary premelting and grain coalescence from molecular-dynamics simulations(1001.0583)

Jan. 4, 2010 cond-mat.mtrl-sci
We describe a molecular dynamics framework for the direct calculation of the short-ranged structural forces underlying grain-boundary premelting and grain-coalescence 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 embedded-atom 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 grain-boundary width at the melting point. It is demonstrated that both types of behavior can be quantitatively described within a sharp-interface thermodynamic formalism involving a width-dependent 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 grain-boundary premelting.
• ### Intrinsic Optical Transition Energies in Carbon Nanotubes(cond-mat/0505004)

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 70-90 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 band-gap renormalization as the surrounding dielectric constant increases.
• ### Intrinsic Optical Transition Energies in Carbon Nanotubes(cond-mat/0505004)

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 70-90 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 band-gap renormalization as the surrounding dielectric constant increases.
• ### Introduction to statistical physics of media processes: Mediaphysics(physics/0506217)

June 30, 2005 physics.soc-ph, cond-mat.soft
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 (mean-field influential propagation of opinions, synergy effects, etc.) factors, considered as fields. Specifically, the interaction and opinion-influence field can be generalized to incorporate important elements of Ising-spin based sociophysical models and kinetic-equation ones. The distributions were described by a Schrodinger-type equation in terms of Green's functions. The developed approach has been applied to a real mass-media efficiency problem for a large company and generally demonstrated very good results despite low initial correlations of factors and the target variable.
• ### Nanoelectromechanics of Piezoresponse Force Microscopy(cond-mat/0408223)

Aug. 10, 2004 cond-mat.mtrl-sci
To achieve quantitative interpretation of Piezoresponse Force Microscopy (PFM), including resolution limits, tip bias- and strain-induced phenomena and spectroscopy, analytical representations for tip-induced 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 sphere-plane 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.
• ### Quantum dissipation induced noncommutative geometry(quant-ph/0301005)

Jan. 3, 2003 quant-ph
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 Aharanov-Bohm effect.
• ### The superconducting gap of \QTR{em}{in situ} $MgB_{2}$ thin films by microwave surface impedance measurements(cond-mat/0212009)

Dec. 1, 2002 cond-mat.supr-con
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 fully-gapped order parameter. The entire temperature dependence of $R_{s}$ is well described by a Mattis-Bardeen 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(astro-ph/0204213)

April 12, 2002 astro-ph
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 ground-based 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 X-ray through radio wavelengths. GOODS incorporates 3.6-24 micron observations from a SIRTF Legacy Program, four-band ACS imaging from an HST Treasury Program, and extensive new ground-based imaging and spectroscopy. GOODS data products will be made available on a rapid time-scale, 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.
• ### Collective Modes in Strongly Coupled Elecronic Bilayer Liquids(cond-mat/0003340)

March 21, 2000 cond-mat
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}$ frequency-moment-sum-rule 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.
• ### Scaling in a cellular automaton model of earthquake faults(cond-mat/0002459)

Feb. 29, 2000 cond-mat.stat-mech
We present theoretical arguments and simulation data indicating that the scaling of earthquake events in models of faults with long-range 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.
• ### Scaling of the distribution of fluctuations of financial market indices(cond-mat/9905305)

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 13-year period 1984-1996, database (ii) contains 8686 daily records for the 35-year period 1962-1996, and database (iii) contains 852 monthly records for the 71-year period 1926-1996. 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 power-law 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 14-year period 1984-97, and database (v) contains 4649 daily records of the Hang-Seng index for the 18-year period 1980-97. We find estimates of $\alpha$ consistent with those describing the distribution of S&P 500 daily-returns. 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.
• ### Discrete molecular dynamics studies of the folding of a protein-like model(cond-mat/9812291)

Dec. 17, 1998 cond-mat.stat-mech, q-bio
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 recently-proposed 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 inter-residue interactions can serve as a tool to study the thermodynamics and kinetics of protein models.
• ### Identifying the protein folding nucleus using molecular dynamics(cond-mat/9812284)

Dec. 17, 1998 cond-mat.stat-mech, q-bio
Molecular dynamics simulations of folding in an off-lattice protein model reveal a nucleation scenario, in which a few well-defined 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.
• ### Total Widths And Slopes From Complex Regge Trajectories(hep-ph/9804270)

April 10, 1998 hep-ph
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 space-like slopes from their time-like 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 Donnachie-Landshoff total cross section analysis.