Exascale Thinking of Particle Energization Problems

Europe/Stockholm
122:026 (Nordita, Stockholm)

122:026

Nordita, Stockholm

Alexandra Veledina (Nordita), Axel Brandenburg (Nordita and University of Boulder), Joonas Nättilä (University of Turku), Juri Poutanen (Nordita and University of Turku), Åke Nordlund (NBI)
Description

Venue

Nordita, Stockholm, Sweden

Scope

Energization of particles is the cornerstone for the most of bright astrophysical phenomena in the Solar system and beyond. Solar flares, pulsars, jets launched by accreting black hole in active galaxies and microquasars and accretion discs around compact objects shine thanks to the dramatic energy release with associated particle acceleration. Which processes lead to energy output, how particles lose their energy to radiation and what are the observational signatures? These are the open questions to be answered by the modern numerical laboratories.

Particle acceleration in these sources is successfully studied in the kinetic approach, while structure and evolution of the systems can only be probed with fluid approximation. Radiation so far has been treated predominantly using phenomenological description. With the growth of computing capabilities, hybrid approaches start to emerge. Coupling of radiation processes with the PiC simulations as well as embedding the PiC methods into an MHD framework becomes feasible, and the next-generation codes are being developed to perform ab-initio calculation on the global scale, setting new horizons in the studies of high-energy astrophysical phenomena.

The workshop aims to discuss the state-of-the-art approaches, challenges for the next-generation hybrid MHD-PiC and radiation-PiC methods and outline the problems to be answered by the new exascale simulations of Nature in silico. It will bring together experts in high-performance computing, numerical methods, plasma physics, and radiative processes.

Invited Speakers

  • Andrei Beloborodov (Columbia University, USA)
  • Benoit Cerutti (Universite Grenoble Alpes, France)
  • Alexander Chen (Columbia University, USA)
  • Mark Dieckmann (Linkoeping University, Sweden)
  • Kyle Parfrey (UC Berkeley, USA)
  • Lorenzo Sironi (Columbia University, USA)
  • Rami Vainio (University of Turku, Finland)
  • Jean-Luc Vay (Lawrence Berkeley National Lab., USA)
  • Sebastian von Alfthan (Center for Scientific Computing, Espoo, Finland)
  • Marija Vranic (Technical Univ. Lisbon, Portugal)

Format

We will have 3-4 presentations a day with ample time for questions and discussions.

The deadline for registration is August 20, 2017.

    • 09:30 10:00
      Registration + Coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 10:00 11:00
      DISPATCH simulations 1h 122:026

      122:026

      Nordita, Stockholm

      We use the adaptive mesh refinement codes RAMSES and DISPATCH to model protoplanetary disks in realistic star formation environments, covering a range from a 40 pc outer scale to an inner scale of 1/10 of an Earth radius (more than 12 orders of magnitude). The simulations are done in four steps, with the first step following individual star formation in a 40 pc GMC model with 120 AU resolution. In the 2nd step, the neighborhoods of several stars with a final system mass of 1-2 solar masses are followed during the accretion process, with a smallest mesh size of 2 AU, sufficient to follow the development of the large scale structure of the accretion disks and their accretion histories. A selection of these disks are then studied over 100-1000 yr, with mesh sizes down to 0.015 AU. Using the new DISPATCH adaptive mesh refinement code, we then study the dynamics of gas, dust, and ensembles of chondrule size pebbles in the disk, following the accretion of pebbles onto planetary embryos surrounded by proto-atmospheres, with three dimensional radiative energy, and a resolution of 600 km. The unprecedented scale range is made possible by the new DISPATCH code, which is explicitly constructed for the exascale computing era, with the ability to scale to essentially unlimited number of cores with sustained performance per core. This is achieved by using task based scheduling, where tasks do MHD, radiative energy transfer, and evolve the positions of billions of pseudo-particles representing dust or pebbles, using time steps restricted only by local condition.
      Speaker: Prof. Åke Nordlund
      Slides
    • 11:00 12:00
      Techniques for well-behaved force-free regions in relativistic MHD simulations 1h 122:026

      122:026

      Nordita, Stockholm

      Speaker: Dr Kyle Parfrey
      Slides
    • 12:00 13:30
      Lunch 1h 30m Albanova restaurant

      Albanova restaurant

      Nordita, Stockholm

    • 13:30 14:00
      Energisation by reconnection within test-particle approach 30m 122:026

      122:026

      Nordita, Stockholm

      We study the role of turbulence in magnetic reconnection, within the framework of magneto-hydrodynamics, using three-dimensional direct numerical simulations. For small turbulent intensity we find that the reconnection rate obeys Sweet-Parker scaling. For large enough turbulent intensity reconnection rate departs significantly from Sweet-Parker behaviour, becomes almost a constant as a function of the Lundquist number. We further study energisation of test- particles in the same setup. We find that the speed of the energised particles obeys a Maxwellian distribution, whose variance also obeys Sweet-Parker scaling for small turbulent intensity but depends weakly on the Lundquist number for large turbulent intensity. Furthermore, the variance is found to increase with the strength of the reconnecting magnetic field.
      Speaker: Dr Dhrubaditya Mitra
      Slides
    • 14:00 15:00
      Radiation emission, radiation reaction and QED processes for PIC codes 1h 122:026

      122:026

      Nordita, Stockholm

      Particle-in-cell codes have been successfully employed to model particle acceleration in laboratory (e. g. laser-wakefield acceleration) and in space (e.g. in collisionless shocks). With the advent of laser technology, one can experiment with very intense fields, that would otherwise be available only in astrophysical objects (e.g. in pulsars or magnetars). At extreme intensities, new physics becomes relevant for modelling laser-plasma interactions, which brings new challenges for PIC development. For example, radiation can be emitted at high frequencies that are not resolved by the simulation grid. If a fraction of energy this radiation carries is negligible, we can compute the output radiation spectra by post-processing the particle trajectories or by using a real- time diagnostics that does not interfere with the PIC loop itself. However, if this radiation accounts for a large fraction of the particle energy, one needs to correct the particle momentum by including a classical description of radiation reaction (e.g. Landau & Lifshitz equation of motion instead of the Lorentz force). One can expect a correct post-processing account of the emitted radiation only if the particle trajectories themselves are correct and the emissivity calculation includes radiation damping corrections. An even greater computational challenge is modelling a quantum regime of emission - when a particle can emit a single photon that carries a large fraction of its energy. Such a photon is treated as an additional particle species, which can propagate through the simulation box and later decay into an electron- positron pair. The new pairs re-accelerate in the laser field, and they emit new photons. Repeated occurrence of this process can induce a so-called “QED cascade”, that generates an exponentially rising number of particles in the simulation box. Macroparticle merging algorithm is then necessary to keep the simulation load to a manageable level. We have developed a merging scheme that resamples particles in the simulation and preserves the particle distribution function. I will discuss the implementation of the above-mentioned computational developments in OSIRIS and show examples of physical problems where they are essential.
      Speaker: Dr Marija Vranich
      Slides
    • 15:00 15:30
      coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 15:30 16:30
      Discussion on MHD numerics 1h 122:026

      122:026

      Nordita, Stockholm

      Speaker: Prof. Åke Nordlund
    • 17:00 19:00
      Reception at Nordita East 2h 122:026

      122:026

      Nordita, Stockholm

    • 09:30 10:00
      Coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 10:00 11:00
      PIC simulation of the thermal pressure-driven expansion of a blast shell into a magnetized ambient medium 1h 122:026

      122:026

      Nordita, Stockholm

      A large gradient of the thermal pressure in a collision-less plasma triggers the formation of a rarefaction wave. Rarefaction waves can accelerate ions to speeds of the order of hundreds to thousands of km/s in laser-plasma experiments and the collision of these fast ion beams with an ambient plasma triggers the formation of shocks. Forthcoming experimental campaigns will introduce a background magnetic field into the ambient plasma that yields a beta value of the order unity and study the magnetized shocks. I will present results from recent PIC simulation studies that investigated such shocks and show how a fast magnetosonic shock forms in such a plasma, which is trailed by a tangential discontinuity, and how a radially expanding blast shell interacts with a uniform background magnetic field.
      Speaker: Dr Mark Dieckmann
      Slides
    • 11:00 12:00
      Numerical methods in WarpX and PICSAR, new tools toward exascale AMR-PIC modeling of relativistic plasmas 1h 122:026

      122:026

      Nordita, Stockholm

      Speaker: Dr Jean-Luc Vay (Lawrence Berkeley National Laboratory)
      Slides
    • 12:00 13:30
      Lunch 1h 30m AlbaNova restaurant

      AlbaNova restaurant

      Nordita, Stockholm

    • 13:30 14:30
      Reconnection in magnetically-dominated plasmas 1h 122:026

      122:026

      Nordita, Stockholm

      In astrophysical relativistic outflows, the magnetic energy density might exceed even the plasma rest-mass energy density, a regime that is dramatically different from laboratory plasmas. We explore the physics of relativistic magnetic reconnection with fully-kinetic particle-in-cell (PIC) simulations, with an emphasis on large-scale computational domains with open boundary conditions. Our results have implications for magnetically-dominated jets, pulsar winds and accretion disk coronae.
      Speaker: Dr Lorenzo Sironi (Columbia University)
      Slides
    • 14:30 15:30
      Discussion on PiC methods 1h 122:026

      122:026

      Nordita, Stockholm

      Speakers: Dr Jean-Luc Vay (Lawrence Berkeley National Laboratory), Dr Mark Dieckmann
    • 15:30 16:00
      Coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 09:30 15:30
      Boat trip to Birka 6h Departs from Stadshusbron

      Departs from Stadshusbron

      Nordita, Stockholm

    • 09:30 10:00
      Coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 10:00 11:00
      Magnetic energy dissipation 1h 122:026

      122:026

      Nordita, Stockholm

      Speaker: Prof. Andrei Beloborodov (Columbia University)
      Slides
    • 11:00 12:00
      Radiative PIC simulations and their applications to pulsars 1h 122:026

      122:026

      Nordita, Stockholm

      I will present some of the recent efforts to model particle acceleration and emission of energetic radiation in pulsars, using global radiative PIC simulations. These studies show that the equatorial current sheet forming in the pulsar wind is the main site of particle acceleration. Relativistic reconnection dissipates magnetic energy which is then efficiently channeled into energetic particles and high-energy synchrotron radiation. Synthetic lightcurves, spectra and polarization present robust features reminiscent of observed gamma-ray pulsars. I will also discuss how particle acceleration proceeds in the wind far from the star but before it reaches the nebula.
      Speaker: Dr Benoit Cerutti
      Slides
    • 12:00 13:30
      Lunch 1h 30m AlbaNova restaurant

      AlbaNova restaurant

      Nordita, Stockholm

    • 13:30 14:30
      Simulating the Twisted Life of Magnetars 1h 122:026

      122:026

      Nordita, Stockholm

      Speaker: Mr Alexander Chen
      Slides
    • 14:30 15:30
      Discussion on radiative PiC 1h 122:026

      122:026

      Nordita, Stockholm

      Speakers: Prof. Andrei Beloborodov (Columbia University), Dr Benoit Cerutti
    • 15:30 16:00
      Coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 18:00 20:30
      Workshop dinner 2h 30m Prinsen restaurant

      Prinsen restaurant

      Nordita, Stockholm

    • 09:30 10:00
      Coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 10:00 11:00
      Vlasiator - enabling large scale hybrid-Vlasov simulations 1h 122:026

      122:026

      Nordita, Stockholm

      Vlasiator is a 6D hybrid-Vlasov simulation code, that has been developed for simulating Earth's magnetosphere at kinetic scales. Here I will present some of the key techniques used for enabling current 5D simulations on Petascale machines, and also present some of the ongoing development efforts. In particular I will present the semi- Lagrangian solver used for propagating the Vlasov fluid, and describe the sparse representation that enables 5D simulations by reducing the problem size by orders of magnitude. I will also discuss current work on a new heterologous domain decomposition, in which the spatial domain of one mpi task is not identical for the field- and Vlasov solvers. Solving the two problems on separate grids allows load balancing to be performed separately and overcome the previously existing scaling problems. Finally I will also discuss experiences in porting Vlasiator to the Xeon Phi processor (Knights Landing), and share experience in running large scale simulations on a large Xeon Phi based supercomputer (Marconi).
      Speaker: Dr Sebastian von Alfthan (CSC)
      Slides
    • 11:00 12:00
      Collisionless shocks in Vlasov-hybrid simulations 1h 122:026

      122:026

      Nordita, Stockholm

      I will present results from hybrid-Vlasov simulations using Vlasiator, the world's first global magnetospheric hybrid Vlasov code. I will consider both the Earth's bow shock and interplanetary shocks, analysing their similarities and differences in terms of global shock structure and kinetic ion processes at and near the shock fronts.
      Speaker: Prof. Rami Vainio
      Slides
    • 12:00 13:30
      Lunch 1h 30m AlbaNova restaurant

      AlbaNova restaurant

      Nordita, Stockholm

    • 13:30 14:00
      Radiative Vlasov simulations 30m 122:026

      122:026

      Nordita, Stockholm

      Speaker: Mr Joonas Nättilä (University of Turku)
      Slides
    • 14:00 15:00
      Discussion on open questions & Nordita program 1h 122:026

      122:026

      Nordita, Stockholm

      Speaker: Prof. Rami Vainio
    • 15:00 15:30
      Coffee 30m 122:026

      122:026

      Nordita, Stockholm

    • 15:30 15:45
      End of Workshop 15m 122:026

      122:026

      Nordita, Stockholm