Simulation of the secondary electrons energy deposition produced by proton beams in PMMA: influence of the target electronic excitation description

We have studied the radial dependence of the energy deposition of the secondary electron generated by swift proton beams incident with energies T = 50 keV–5 MeV on poly(methylmethacrylate) (PMMA). Two different approaches have been used to model the electronic excitation spectrum of PMMA through its energy loss function (ELF), namely the extended-Drude ELF and the Mermin ELF. The singly differential cross section and the total cross section for ionization, as well as the average energy of the generated secondary electrons, show sizeable differences at T ⩽ 0.1 MeV when evaluated with these two ELF models. In order to know the radial distribution around the proton track of the energy deposited by the cascade of secondary electrons, a simulation has been performed that follows the motion of the electrons through the target taking into account both the inelastic interactions (via electronic ionizations and excitations as well as electron-phonon and electron trapping by polaron creation) and the elastic interactions. The radial distribution of the energy deposited by the secondary electrons around the proton track shows notable differences between the simulations performed with the extended-Drude ELF or the Mermin ELF, being the former more spread out (and, therefore, less peaked) than the latter. The highest intensity and sharpness of the deposited energy distributions takes place for proton beams incident with T ~ 0.1–1 MeV. We have also studied the influence in the radial distribution of deposited energy of using a full energy distribution of secondary electrons generated by proton impact or using a single value (namely, the average value of the distribution); our results show that differences between both simulations become important for proton energies larger than ~0.1 MeV. The results presented in this work have potential applications in materials science, as well as hadron therapy (due to the use of PMMA as a tissue phantom) in order to properly consider the generation of electrons by proton beams and their subsequent transport and energy deposition through the target in nanometric scales.

SiD Letter of Intent

Letter of intent describing SiD (Silicon Detector) for consideration by the International Linear Collider IDAG panel. This detector concept is founded on the use of silicon detectors for vertexing, tracking, and electromagnetic calorimetry. The detector has been cost-optimized as a general-purpose detector for a 500 GeV electron-positron linear collider.

SiD Letter of Intent

Letter of intent describing SiD (Silicon Detector) for consideration by the International Linear Collider IDAG panel. This detector concept is founded on the use of silicon detectors for vertexing, tracking, and electromagnetic calorimetry. The detector has been cost-optimized as a general-purpose detector for a 500 GeV electron-positron linear collider.

Ortho-positronium annihilation: steps toward computing the first order radiative corrections

The 0.2% experimental accuracy of the 1968 Beers and Hughes measurement of the annihilation lifetime of ortho-positronium motivates the attempt to compute the first order quantum electrodynamic corrections to this lifetime. The theoretical problems arising in this computation are here studied in detail up to the point of preparing the necessary computer programs and using them to carry out some of the less demanding steps -- but the computation has not yet been completed. Analytic evaluation of the contributing Feynman diagrams is superior to numerical evaluation, and for this process can be carried out with the aid of the Reduce algebra manipulation computer program.

The relation of the positronium decay rate to the electronpositron annihilation-in-flight amplitude is derived in detail, and it is shown that at threshold annihilation-in-flight, Coulomb divergences appear while infrared divergences vanish. The threshold Coulomb divergences in the amplitude cancel against like divergences in the modulating continuum wave function.

Using the lowest order diagrams of electron-positron annihilation into three photons as a test case, various pitfalls of computer algebraic manipulation are discussed along with ways of avoiding them. The computer manipulation of artificial polynomial expressions is preferable to the direct treatment of rational expressions, even though redundant variables may have to be introduced.

Special properties of the contributing Feynman diagrams are discussed, including the need to restore gauge invariance to the sum of the virtual photon-photon scattering box diagrams by means of a finite subtraction.

A systematic approach to the Feynman-Brown method of Decomposition of single loop diagram integrals with spin-related tensor numerators is developed in detail. This approach allows the Feynman-Brown method to be straightforwardly programmed in the Reduce algebra manipulation language.

The fundamental integrals needed in the wake of the application of the Feynman-Brown decomposition are exhibited and the methods which were used to evaluate them -- primarily dis persion techniques are briefly discussed.

Finally, it is pointed out that while the techniques discussed have permitted the computation of a fair number of the simpler integrals and diagrams contributing to the first order correction of the ortho-positronium annihilation rate, further progress with the more complicated diagrams and with the evaluation of traces is heavily contingent on obtaining access to adequate computer time and core capacity.

A measurement of B0-B0 mixing in e+e- annihilation at 29 GeV

This thesis describes a measurement of B⁰- B⁰ mixing in events produced by electron-positron annihilation at a center of mass energy of 29 GeV. The data were taken by the Mark II detector in the PEP storage ring at the Stanford Linear Accelerator Center between 1981 and 1987, and correspond to a total integrated luminosity of 224pb^-1.

We used a new method, based on the kinematics of hadronic events containing two leptons, to provide a measurement of the probability, x, that a hadron, initially containing a b (b) quark decays to a positive (negative) lepton to be X = 0.17^+0.15_-0.08, with 90% confidence level upper and lower limits of 0.38 and 0.06, respectively, including all estimated systematic errors. Because of the good separation of signal and background, this result is relatively insensitive to various systematic effects which have complicated previous measurements.

We interpret this result as evidence for the mixing of neutral B mesons. Based on existing B⁰_d mixing rate measurements, and some assumptions about the fractions of B⁰_d and B⁰_s mesons present in the data, this result favors maximal mixing of B⁰_s mesons, although it cannot rule out zero B⁰_s mixing at the 90% confidence level.

Supernovae explosions induced by pair production instability

Stars with a core mass greater than about 30 M_⊙ become dynamically unstable due to electron-positron pair production when their central temperature reaches 1.5-2.0 x 10^{9 0}K. The collapse and subsequent explosion of stars with core masses of 45, 52, and 60 M_⊙ is calculated. The range of the final velocity of expansion (3,400 – 8,500 km/sec) and of the mass ejected (1 – 40 M_⊙) is comparable to that observed for type II supernovae.

An implicit scheme of hydrodynamic difference equations (stable for large time steps) used for the calculation of the evolution is described.

For fast evolution the turbulence caused by convective instability does not produce the zero entropy gradient and perfect mixing found for slower evolution. A dynamical model of the convection is derived from the equations of motion and then incorporated into the difference equations.

I. Baryon-antibaryon phase transition at high temperature. II. Inclusive virtual photon-hadron reactions in the parton model

Part I

Present experimental data on nucleon-antinucleon scattering allow a study of the possibility of a phase transition in a nucleon-antinucleon gas at high temperature. Estimates can be made of the general behavior of the elastic phase shifts without resorting to theoretical derivation. A phase transition which separates nucleons from antinucleons is found at about 280 MeV in the approximation of the second virial coefficient to the free energy of the gas.

Part II

The parton model is used to derive scaling laws for the hadrons observed in deep inelastic electron-nucleon scattering which lie in the fragmentation region of the virtual photon. Scaling relations are obtained in the Bjorken and Regge regions. It is proposed that the distribution functions become independent of both q² and ν where the Bjorken and Regge regions overlap. The quark density functions are discussed in the limit x→1 for the nucleon octet and the pseudoscalar mesons. Under certain plausible assumptions it is found that only one or two quarks of the six types of quarks and antiquarks have an appreciable density function in the limit x→1. This has implications for the quark fragmentation functions near the large momentum boundary of their fragmentation region. These results are used to propose a method of measuring the proton and neutron quark density functions for all x by making measurements on inclusively produced hadrons in electroproduction only. Implications are also discussed for the hadrons produced in electron-positron annihilation.

Modulated electrostatic modes in pair plasmas: Modulational stability profile and envelope excitations

A pair plasma consisting of two types of ions, possessing equal masses and opposite charges, is considered. The nonlinear propagation of modulated electrostatic wave packets is studied by employing a two-fluid plasma model. Considering propagation parallel to the external magnetic field, two distinct electrostatic modes are obtained, namely a quasiacoustic lower moddfe and a Langmuir-like, as optic-type upper one, in agreement with experimental observations and theoretical predictions. Considering small yet weakly nonlinear deviations from equilibrium, and adopting a multiple-scale technique, the basic set of model equations is reduced to a nonlinear Schrodinger equation for the slowly varying electric field perturbation amplitude. The analysis reveals that the lower (acoustic) mode is stable and may propagate in the form of a dark-type envelope soliton (a void) modulating a carrier wave packet, while the upper linear mode is intrinsically unstable, and may favor the formation of bright-type envelope soliton (pulse) modulated wave packets. These results are relevant to recent observations of electrostatic waves in pair-ion (fullerene) plasmas, and also with respect to electron-positron plasma emission in pulsar magnetospheres. (c) 2006 American Institute of Physics.

Electrostatic mode envelope excitations in e-p-i plasmas - application in warm pair ion plasmas with a small fraction of stationary ions

The nonlinear propagation of amplitude-modulated electrostatic wavepackets in an electron-positron-ion (e-p-i) plasma is considered, by employing a two-fluid plasma model. Considering propagation parallel to the external magnetic field, two distinct electrostatic modes are obtained, namely a quasi-thermal acoustic-like lower mode and a Langmuir-like optic-type upper one. These results equally apply in warm pair ion ( e. g. fullerene) plasmas contaminated by a small fraction of stationary ions ( or dust), in agreement with experimental observations and theoretical predictions in pair plasmas. Considering small yet weakly nonlinear deviations from equilibrium, and adopting a multiple-scales perturbation technique, the basic set of model equations is reduced to a nonlinear Schrodinger (NLS) equation for the slowly varying electric field perturbation amplitude. The analysis reveals that the lower ( acoustic) mode is mostly stable for large wavelengths, and may propagate in the form of a dark-type envelope soliton ( a void) modulating a carrier wavepacket, while the upper linear mode is intrinsically unstable, and thus favours the formation of bright-type envelope soliton ( pulse) modulated wavepackets. The stability ( instability) range for the acoustic ( Langmuir-like optic) mode shifts to larger wavenumbers as the positive-to-negative ion temperature ( density) ratio increases. These results may be of relevance in astrophysical contexts, where e-p-i plasmas are encountered, and may also serve as prediction of the behaviour of doped ( or dust-contaminated) fullerene plasmas, in the laboratory.

Acoustic solitary waves in dusty and/or multi-ion plasmas with cold, adiabatic, and hot constituents

Large nonlinear acoustic waves are discussed in a four-component plasma, made up of two superhot isothermal species, and two species with lower thermal velocities, being, respectively, adiabatic and cold. First a model is considered in which the isothermal species are electrons and ions, while the cooler species are positive and/or negative dust. Using a Sagdeev pseudopotential formalism, large dust-acoustic structures have been studied in a systematic way, to delimit the compositional parameter space in which they can be found, without restrictions on the charges and masses of the dust species and their charge signs. Solitary waves can only occur for nonlinear structure velocities smaller than the adiabatic dust thermal velocity, leading to a novel dust-acoustic-like mode based on the interplay between the two dust species. If the cold and adiabatic dust are oppositely charged, only solitary waves exist, having the polarity of the cold dust, their parameter range being limited by infinite compression of the cold dust. However, when the charges of the cold and adiabatic species have the same sign, solitary structures are limited for increasing Mach numbers successively by infinite cold dust compression, by encountering the adiabatic dust sonic point, and by the occurrence of double layers. The latter have, for smaller Mach numbers, the same polarity as the charged dust, but switch at the high Mach number end to the opposite polarity. Typical Sagdeev pseudopotentials and solitary wave profiles have been presented. Finally, the analysis has nowhere used the assumption that the dust would be much more massive than the ions and hence, one or both dust species can easily be replaced by positive and/or negative ions and the conclusions will apply to that plasma model equally well. This would cover a number of different scenarios, such as, for example, very hot electrons and ions, together with a mix of adiabatic ions and dust (of either polarity) or a very hot electron-positron mix, together with a two-ion mix or together with adiabatic ions and cold dust (both of either charge sign), to name but some of the possible plasma compositions.

Low-frequency electrostatic defect mode in doped pair-ion plasmas

The nonlinear amplitude modulation dynamics of electrostatic oscillations of massive charged defects in a three-component pair plasma is investigated; i.e. doped pair-ion plasmas (anticipating the injection of a massive charged component in the background; e.g. in fullerene experiments). Ton-acoustic oscillations in electron-positron-ion (e-p-i) plasmas are also covered, in the appropriate limit. Linear and nonlinear effects (MI, envelope modes) are discussed. The role of the temperature and density ratio between the pair species is stressed.

Nonlinear Dynamics of Rotating Multi-Component Pair Plasmas and e-p-i Plasmas

The propagation of small amplitude stationary profile nonlinear electrostatic excitations in a pair plasma is investigated, mainly drawing inspiration from experiments on fullerene pair-ion plasmas. Two distinct pair ion species are considered of opposite polarity and same mass, in addition to a massive charged background species, which is assumed to be stationary, given the frequency scale of interest. In the pair-ion context, the third species is thought of as a background defect (e.g. charged dust) component. On the other hand, the model also applies formally to electron-positron-ion (e-p-i) plasmas, if one neglects electron-positron annihilation. A two-fluid plasma model is employed, incorporating both Lorentz and Coriolis forces, thus taking into account the interplay between the gyroscopic (Larmor) frequency ?c and the (intrinsic) plasma rotation frequency O0. By employing a multi-dimensional reductive perturbation technique, a Zakharov-Kuznetsov (ZK) type equation is derived for the evolution of the electric potential perturbation. Assuming an arbitrary direction of propagation, with respect to the magnetic field, we derive the exact form of nonlinear solutions, and study their characteristics. A parametric analysis is carried out, as regards the effect of the dusty plasma composition (background number density), species temperature(s) and the relative strength of rotation to Larmor frequencies. It is shown that the Larmor and mechanical rotation affect the pulse dynamics via a parallel-to-transverse mode coupling diffusion term, which in fact diverges at ?c ? ±2O0. Pulses collapse at this limit, as nonlinearity fails to balance dispersion. The analysis is complemented by investigating critical plasma compositions, in fact near-symmetric (T- ˜ T+) “pure” (n- ˜ n+) pair plasmas, i.e. when the concentration of the 3rd background species is negligible, case in which the (quadratic) nonlinearity vanishes, so one needs to resort to higher order nonlinear theory. A modified ZK equation is derived and analyzed. Our results are of relevance in pair-ion (fullerene) experiments and also potentially in astrophysical environments, e.g. in pulsars.

Slowly fading super-luminous supernovae that are not pair-instability explosions

Super-luminous supernovae that radiate more than 10 44 ergs per second at their peak luminosity have recently been discovered in faint galaxies at redshifts of 0.1-4. Some evolve slowly, resembling models of 'pair-instability' supernovae. Such models involve stars with original masses 140-260 times that of the Sun that now have carbon-oxygen cores of 65-130 solar masses. In these stars, the photons that prevent gravitational collapse are converted to electron-positron pairs, causing rapid contraction and thermonuclear explosions. Many solar masses of 56 Ni are synthesized; this isotope decays to 56 Fe via 56 Co, powering bright light curves. Such massive progenitors are expected to have formed from metal-poor gas in the early Universe. Recently, supernova 2007bi in a galaxy at redshift 0.127 (about 12 billion years after the Big Bang) with a metallicity one-third that of the Sun was observed to look like a fading pair-instability supernova. Here we report observations of two slow-to-fade super-luminous supernovae that show relatively fast rise times and blue colours, which are incompatible with pair-instability models. Their late-time light-curve and spectral similarities to supernova 2007bi call the nature of that event into question. Our early spectra closely resemble typical fast-declining super-luminous supernovae, which are not powered by radioactivity. Modelling our observations with 10-16 solar masses of magnetar-energized ejecta demonstrates the possibility of a common explosion mechanism. The lack of unambiguous nearby pair-instability events suggests that their local rate of occurrence is less than 6 × 10 -6 times that of the core-collapse rate. © 2013 Macmillan Publishers Limited. All rights reserved.

ON THE DOUBLE NATURED SOLUTIONS OF THE TWO-TEMPERATURE EXTERNAL SOFT PHOTON COMPTONIZED ACCRETION DISKS

We have analyzed pair production in the innermost region of a two-temperature external soft photon Comptonized accretion disk. We have shown that, if the viscosity parameter is greater than a critical value alpha(c), the solution to the disk equation is double valued: one, advection dominated, and the other, radiation dominated. When alpha <= alpha(c), the accretion rate has to satisfy (m) over dot(1) <= (m) over dot <= (m) over dot(c) in order to have two steady-state solutions. It is shown that these critical parameters (m) over dot(1), (m) over dot(c) are functions of r, alpha, and theta(e), and alpha(c) is a function of r and theta(e). Depending on the combination of the parameters, the advection-dominated solution may not be physically consistent. It is also shown that the electronic temperature is maximum at the onset of the thermal instability, from which results this inner region. These solutions are stable against perturbations in the electron temperature and in the density of pairs.

Dilepton mass spectra in p plus p collisions at root s=200 GeV

PHENIX has measured the electron-positron pair mass spectrum from 0 to 8 GeV/c(2) in p + p collisions at root s = 200 GeV. The contributions from light meson decays to e(+)e(-) pairs have been determined based on measurements of hadron production cross sections by PHENIX. Within the systematic uncertainty of similar to 20% they account for all e(+)e(-) pairs in the mass region below similar to 1 GeV/c(2). The e(+)e(-) pair yield remaining after subtracting these contributions is dominated by semileptonic decays of charmed hadrons correlated through flavor conservation. Using the spectral shape predicted by PYTHIA, we estimate the charm production cross section to be 544 +/- 39(stat) +/- 142(syst) +/- 200(model) pb. which is consistent with QCD calculations and measurements of single leptons by PHENIX. (C) 2008 Elsevier BV. All rights reserved.