Efficient calculation of local dose distribution for response modelling in proton and ion beams

We present an algorithm for fast and accurate computation of the local dose distribution in MeV beams of protons, carbon ions or other heavy-charged particles. It uses compound Poisson-process modelling of track interaction and succesive convolutions for fast computation. It can handle mixed particle fields over a wide range of fluences. Since the local dose distribution is the essential part of several approaches to model detector efficiency or cellular response it has potential use in ion-beam dosimetry and radiotherapy.Comment: 9 pages, 3 figure

W Plus Multiple Jets at the LHC with High Energy Jets

We study the production of a W boson in association with n hard QCD jets (for n>=2), with a particular emphasis on results relevant for the Large Hadron Collider (7 TeV and 8 TeV). We present predictions for this process from High Energy Jets, a framework for all-order resummation of the dominant contributions from wide-angle QCD emissions. We first compare predictions against recent ATLAS data and then shift focus to observables and regions of phase space where effects beyond NLO are expected to be large.Comment: 19 pages, 9 figure

Linear response separation of a solid into atomic constituents: Li, Al, and their evolution under pressure

We present the first realization of the generalized pseudoatom concept introduced by Ball, and adopt the name enatom to minimize confusion. This enatom, which consists of a unique decomposition of the total charge density (or potential) of any solid into a sum of overlapping atomiclike contributions that move rigidly with the nuclei to first order, is calculated using (numerical) linear response methods, and is analyzed for both fcc Li and Al at pressures of 0, 35, and 50 GPa. These two simple fcc metals (Li is fcc and a good superconductor in the 20-40 GPa range) show different physical behaviors under pressure, which reflects the increasing covalency in Li and the lack of it in Al. The nonrigid (deformation) parts of the enatom charge and potential have opposite signs in Li and Al; they become larger under pressure only in Li. These results establish a method of construction of the enatom, whose potential can be used to obtain a real-space understanding of the vibrational properties and electron-phonon interaction in solids.Comment: 13 pages, 9 figures, 1 table, V2: fixed problem with Fig. 7, V3: minor correction

Enumeration of chord diagrams on many intervals and their non-orientable analogs

Two types of connected chord diagrams with chord endpoints lying in a collection of ordered and oriented real segments are considered here: the real segments may contain additional bivalent vertices in one model but not in the other. In the former case, we record in a generating function the number of fatgraph boundary cycles containing a fixed number of bivalent vertices while in the latter, we instead record the number of boundary cycles of each fixed length. Second order, non-linear, algebraic partial differential equations are derived which are satisfied by these generating functions in each case giving efficient enumerative schemes. Moreover, these generating functions provide multi-parameter families of solutions to the KP hierarchy. For each model, there is furthermore a non-orientable analog, and each such model likewise has its own associated differential equation. The enumerative problems we solve are interpreted in terms of certain polygon gluings. As specific applications, we discuss models of several interacting RNA molecules. We also study a matrix integral which computes numbers of chord diagrams in both orientable and non-orientable cases in the model with bivalent vertices, and the large-N limit is computed using techniques of free probability.Comment: 23 pages, 7 figures; revised and extended versio

Negative Specific Heat in a Quasi-2D Generalized Vorticity Model

Negative specific heat is a dramatic phenomenon where processes decrease in temperature when adding energy. It has been observed in gravo-thermal collapse of globular clusters. We now report finding this phenomenon in bundles of nearly parallel, periodic, single-sign generalized vortex filaments in the electron magnetohydrodynamic (EMH) model for the unbounded plane under strong magnetic confinement. We derive the specific heat using a steepest descent method and a mean field property. Our derivations show that as temperature increases, the overall size of the system increases exponentially and the energy drops. The implication of negative specific heat is a runaway reaction, resulting in a collapsing inner core surrounded by an expanding halo of filaments.Comment: 12 pages, 3 figures; updated with revision

Interplay between nanometer-scale strain variations and externally applied strain in graphene

We present a molecular modeling study analyzing nanometer-scale strain variations in graphene as a function of externally applied tensile strain. We consider two different mechanisms that could underlie nanometer-scale strain variations: static perturbations from lattice imperfections of an underlying substrate and thermal fluctuations. For both cases we observe a decrease in the out-of-plane atomic displacements with increasing strain, which is accompanied by an increase in the in-plane displacements. Reflecting the non-linear elastic properties of graphene, both trends together yield a non-monotonic variation of the total displacements with increasing tensile strain. This variation allows to test the role of nanometer-scale strain variations in limiting the carrier mobility of high-quality graphene samples