Fine-sorting One-dimensional Particle-In-Cell Algorithm with Monte-Carlo Collisions on a Graphics Processing Unit

Particle-in-cell (PIC) simulations with Monte-Carlo collisions are used in plasma science to explore a variety of kinetic effects. One major problem is the long run-time of such simulations. Even on modern computer systems, PIC codes take a considerable amount of time for convergence. Most of the computations can be massively parallelized, since particles behave independently of each other within one time step. Current graphics processing units (GPUs) offer an attractive means for execution of the parallelized code. In this contribution we show a one-dimensional PIC code running on Nvidia GPUs using the CUDA environment. A distinctive feature of the code is that size of the cells that the code uses to sort the particles with respect to their coordinates is comparable to size of the grid cells used for discretization of the electric field. Hence, we call the corresponding algorithm "fine-sorting". Implementation details and optimization of the code are discussed and the speed-up compared to classical CPU approaches is computed

Simulations of electromagnetic effects in high frequency capacitively coupled discharges using the Darwin approximation

The Darwin approximation is investigated for its possible use in simulation of electromagnetic effects in large size, high frequency capacitively coupled discharges. The approximation is utilized within the framework of two different fluid models which are applied to typical cases showing pronounced standing wave and skin effects. With the first model it is demonstrated that Darwin approximation is valid for treatment of such effects in the range of parameters under consideration. The second approach, a reduced nonlinear Darwin approximation-based model, shows that the electromagnetic phenomena persist in a more realistic setting. The Darwin approximation offers a simple and efficient way of carrying out electromagnetic simulations as it removes the Courant condition plaguing explicit electromagnetic algorithms and can be implemented as a straightforward modification of electrostatic algorithms. The algorithm described here avoids iterative schemes needed for the divergence cleaning and represents a fast and efficient solver, which can be used in fluid and kinetic models for self-consistent description of technical plasmas exhibiting certain electromagnetic activity

Electron dynamics in planar radio frequency magnetron plasmas: I. The mechanism of Hall heating and the {\mu}-mode

The electron dynamics and the mechanisms of power absorption in radio-frequency (RF) driven, magnetically enhanced capacitively coupled plasmas (MECCPs) at low pressure are investigated. The device in focus is a geometrically asymmetric cylindrical magnetron with a radially nonuniform magnetic field in axial direction and an electric field in radial direction. The dynamics is studied analytically using the cold plasma model and a single-particle formalism, and numerically with the inhouse energy and charge conserving particle-in-cell/Monte Carlo collisions code ECCOPIC1S-M. It is found that the dynamics differs significantly from that of an unmagnetized reference discharge. In the magnetized region in front of the powered electrode, an enhanced electric field arises during sheath expansion and a reversed electric field during sheath collapse. Both fields are needed to ensure discharge sustaining electron transport against the confining effect of the magnetic field. The corresponding azimuthal ExB-drift can accelerate electrons into the inelastic energy range which gives rise to a new mechanism of RF power dissipation. It is related to the Hall current and is different in nature from Ohmic heating, as which it has been classified in previous literature. The new heating is expected to be dominant in many magnetized capacitively coupled discharges. It is proposed to term it the "{\mu}-mode" to separate it from other heating modes

Ion energy distribution functions behind the sheaths of magnetized and non magnetized radio frequency discharges

The effect of a magnetic field on the characteristics of capacitively coupled radio frequency discharges is investigated and found to be substantial. A one-dimensional particle-in-cell simulation shows that geometrically symmetric discharges can be asymmetrized by applying a spatially inhomogeneous magnetic field. This effect is similar to the recently discovered electrical asymmetry effect. Both effects act independently, they can work in the same direction or compensate each other. Also the ion energy distribution functions at the electrodes are strongly affected by the magnetic field, although only indirectly. The field influences not the dynamics of the sheath itself but rather its operating conditions, i.e., the ion flux through it and voltage drop across it. To support this interpretation, the particle-in-cell results are compared with the outcome of the recently proposed ensemble-in-spacetime algorithm. Although that scheme resolves only the sheath and neglects magnetization, it is able to reproduce the ion energy distribution functions with very good accuracy, regardless of whether the discharge is magnetized or not

Standalone vertex finding in the ATLAS muon spectrometer

A dedicated reconstruction algorithm to find decay vertices in the ATLAS muon spectrometer is presented. The algorithm searches the region just upstream of or inside the muon spectrometer volume for multi-particle vertices that originate from the decay of particles with long decay paths. The performance of the algorithm is evaluated using both a sample of simulated Higgs boson events, in which the Higgs boson decays to long-lived neutral particles that in turn decay to bbar b final states, and pp collision data at √s = 7 TeV collected with the ATLAS detector at the LHC during 2011

Measurements of Higgs boson production and couplings in diboson final states with the ATLAS detector at the LHC

Measurements are presented of production properties and couplings of the recently discovered Higgs boson using the decays into boson pairs, H →γ γ, H → Z Z∗ →4l and H →W W∗ →lνlν. The results are based on the complete pp collision data sample recorded by the ATLAS experiment at the CERN Large Hadron Collider at centre-of-mass energies of √s = 7 TeV and √s = 8 TeV, corresponding to an integrated luminosity of about 25 fb−1. Evidence for Higgs boson production through vector-boson fusion is reported. Results of combined fits probing Higgs boson couplings to fermions and bosons, as well as anomalous contributions to loop-induced production and decay modes, are presented. All measurements are consistent with expectations for the Standard Model Higgs boson

Measurement of the polarisation of W bosons produced with large transverse momentum in pp collisions at sqrt(s) = 7 TeV with the ATLAS experiment

This paper describes an analysis of the angular distribution of W->enu and W->munu decays, using data from pp collisions at sqrt(s) = 7 TeV recorded with the ATLAS detector at the LHC in 2010, corresponding to an integrated luminosity of about 35 pb^-1. Using the decay lepton transverse momentum and the missing transverse energy, the W decay angular distribution projected onto the transverse plane is obtained and analysed in terms of helicity fractions f0, fL and fR over two ranges of W transverse momentum (ptw): 35 < ptw < 50 GeV and ptw > 50 GeV. Good agreement is found with theoretical predictions. For ptw > 50 GeV, the values of f0 and fL-fR, averaged over charge and lepton flavour, are measured to be : f0 = 0.127 +/- 0.030 +/- 0.108 and fL-fR = 0.252 +/- 0.017 +/- 0.030, where the first uncertainties are statistical, and the second include all systematic effects.Comment: 19 pages plus author list (34 pages total), 9 figures, 11 tables, revised author list, matches European Journal of Physics C versio

Observation of a new chi_b state in radiative transitions to Upsilon(1S) and Upsilon(2S) at ATLAS

The chi_b(nP) quarkonium states are produced in proton-proton collisions at the Large Hadron Collider (LHC) at sqrt(s) = 7 TeV and recorded by the ATLAS detector. Using a data sample corresponding to an integrated luminosity of 4.4 fb^-1, these states are reconstructed through their radiative decays to Upsilon(1S,2S) with Upsilon->mu+mu-. In addition to the mass peaks corresponding to the decay modes chi_b(1P,2P)->Upsilon(1S)gamma, a new structure centered at a mass of 10.530+/-0.005 (stat.)+/-0.009 (syst.) GeV is also observed, in both the Upsilon(1S)gamma and Upsilon(2S)gamma decay modes. This is interpreted as the chi_b(3P) system.Comment: 5 pages plus author list (18 pages total), 2 figures, 1 table, corrected author list, matches final version in Physical Review Letter