1
|
Ziaja B, Bekx JJ, Masek M, Medvedev N, Lipp V, Saxena V, Stransky M. Application of Boltzmann kinetic equations to model X-ray-created warm dense matter and plasma. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220216. [PMID: 37393933 PMCID: PMC10876064 DOI: 10.1098/rsta.2022.0216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Accepted: 03/28/2023] [Indexed: 07/04/2023]
Abstract
In this review, we describe the application of Boltzmann kinetic equations for modelling warm dense matter and plasma formed after irradiation of solid materials with intense femtosecond X-ray pulses. Classical Boltzmann kinetic equations are derived from the reduced N-particle Liouville equations. They include only single-particle densities of ions and free electrons present in the sample. The first version of the Boltzmann kinetic equation solver was completed in 2006. It could model non-equilibrium evolution of X-ray-irradiated finite-size atomic systems. In 2016, the code was adapted to study plasma created from X-ray-irradiated materials. Additional extension of the code was then also performed, enabling simulations in the hard X-ray irradiation regime. In order to avoid treatment of a very high number of active atomic configurations involved in the excitation and relaxation of X-ray-irradiated materials, an approach called 'predominant excitation and relaxation path' (PERP) was introduced. It limited the number of active atomic configurations by following the sample evolution only along most PERPs. The performance of the Boltzmann code is illustrated in the examples of X-ray-heated solid carbon and gold. Actual model limitations and further model developments are discussed. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.
Collapse
Affiliation(s)
- Beata Ziaja
- Center for Free-Electron Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany
- Institute of Nuclear Physics, Polish Academy of Sciences, Radzikowskiego 152, 31-342 Krakow, Poland
| | - John Jasper Bekx
- Center for Free-Electron Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany
| | - Martin Masek
- Institute of Physics, Czech Academy of Sciences, Na Slovance 2,182 21 Prague 8, Czech Republic
| | - Nikita Medvedev
- Institute of Physics, Czech Academy of Sciences, Na Slovance 2,182 21 Prague 8, Czech Republic
- Institute of Plasma Physics, Czech Academy of Sciences, Za Slovankou 3, 182 00 Prague 8, Czech Republic
| | - Vladimir Lipp
- Center for Free-Electron Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany
- Institute of Nuclear Physics, Polish Academy of Sciences, Radzikowskiego 152, 31-342 Krakow, Poland
| | - Vikrant Saxena
- Center for Free-Electron Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany
- Department of Physics, Indian Institute of Technology Delhi,New Delhi 110016, India
| | - Michal Stransky
- Institute of Physics, Czech Academy of Sciences, Na Slovance 2,182 21 Prague 8, Czech Republic
- European XFEL, Holzkoppel 4, 22869 Schenefeld, Germany
| |
Collapse
|
2
|
Density functional tight binding approach utilized to study X-ray-induced transitions in solid materials. Sci Rep 2022; 12:1551. [PMID: 35091574 PMCID: PMC8799736 DOI: 10.1038/s41598-022-04775-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2021] [Accepted: 12/02/2021] [Indexed: 01/17/2023] Open
Abstract
Intense X-ray pulses from free-electron lasers can trigger ultrafast electronic, structural and magnetic transitions in solid materials, within a material volume which can be precisely shaped through adjustment of X-ray beam parameters. This opens unique prospects for material processing with X rays. However, any fundamental and applicational studies are in need of computational tools, able to predict material response to X-ray radiation. Here we present a dedicated computational approach developed to study X-ray induced transitions in a broad range of solid materials, including those of high chemical complexity. The latter becomes possible due to the implementation of the versatile density functional tight binding code DFTB+ to follow band structure evolution in irradiated materials. The outstanding performance of the implementation is demonstrated with a comparative study of XUV induced graphitization in diamond.
Collapse
|
3
|
Sprenkle T, Dodson A, McKnight Q, Spencer R, Bergeson S, Diaw A, Murillo MS. Ion friction at small values of the Coulomb logarithm. Phys Rev E 2019; 99:053206. [PMID: 31212549 DOI: 10.1103/physreve.99.053206] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2018] [Indexed: 11/07/2022]
Abstract
Transport properties of high-energy-density plasmas are influenced by the ion collision rate. Traditionally, this rate involves the Coulomb logarithm, lnΛ. Typical values of lnΛ are ≈10-20 in kinetic theories where transport properties are dominated by weak-scattering events caused by long-range forces. The validity of these theories breaks down for strongly coupled plasmas, when lnΛ is of order one. We present measurements and simulations of collision data in strongly coupled plasmas when lnΛ is small. Experiments are carried out in the first dual-species ultracold neutral plasma (UNP), using Ca^{+} and Yb^{+} ions. We find strong collisional coupling between the different ion species in the bulk of the plasma. We simulate the plasma using a two-species fluid code that includes Coulomb logarithms derived from either a screened Coulomb potential or a the potential of mean force. We find generally good agreement between the experimental measurements and the simulations. With some improvements, the mixed Ca^{+} and Yb^{+} dual-species UNP will be a promising platform for testing theoretical expressions for lnΛ and collision cross-sections from kinetic theories through measurements of energy relaxation, stopping power, two-stream instabilities, and the evolution of sculpted distribution functions in an idealized environment in which the initial temperatures, densities, and charge states are accurately known.
Collapse
Affiliation(s)
- Tucker Sprenkle
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Adam Dodson
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Quinton McKnight
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Ross Spencer
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Scott Bergeson
- Department of Physics and Astronomy, Brigham Young University, Provo, Utah 84602, USA
| | - Abdourahmane Diaw
- Computational Physics and Methods Group, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Michael S Murillo
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
4
|
Yabuuchi T, Kon A, Inubushi Y, Togahi T, Sueda K, Itoga T, Nakajima K, Habara H, Kodama R, Tomizawa H, Yabashi M. An experimental platform using high-power, high-intensity optical lasers with the hard X-ray free-electron laser at SACLA. JOURNAL OF SYNCHROTRON RADIATION 2019; 26:585-594. [PMID: 30855271 PMCID: PMC6412175 DOI: 10.1107/s1600577519000882] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Accepted: 01/17/2019] [Indexed: 05/25/2023]
Abstract
An experimental platform using X-ray free-electron laser (XFEL) pulses with high-intensity optical laser pulses is open for early users' experiments at the SACLA XFEL facility after completion of the commissioning. The combination of the hard XFEL and the high-intensity laser provides capabilities to open new frontiers of laser-based high-energy-density science. During the commissioning phase, characterization of the XFEL and the laser at the platform has been carried out for the combinative utilization as well as the development of instruments and basic diagnostics for user experiments. An overview of the commissioning and the current capabilities of the experimental platform is presented.
Collapse
Affiliation(s)
| | - Akira Kon
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
- Japan Synchrotoron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
| | - Yuichi Inubushi
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
- Japan Synchrotoron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
| | - Tadashi Togahi
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
- Japan Synchrotoron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
| | - Keiichi Sueda
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
| | - Toshiro Itoga
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
- Japan Synchrotoron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
| | - Kyo Nakajima
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
| | - Hideaki Habara
- Graduate School of Engineering, 2-1 Yamada-oka, Suita, Osaka 565-0871, Japan
| | - Ryosuke Kodama
- Graduate School of Engineering, 2-1 Yamada-oka, Suita, Osaka 565-0871, Japan
- Institute of Laser Engineering, Osaka University, 2-6 Yamada-oka, Suita, Osaka 565-0871, Japan
| | - Hiromitsu Tomizawa
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
- Japan Synchrotoron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
| | - Makina Yabashi
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan
- Japan Synchrotoron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan
| |
Collapse
|
5
|
Abdullah MM, Jurek Z, Son SK, Santra R. Molecular-dynamics approach for studying the nonequilibrium behavior of x-ray-heated solid-density matter. Phys Rev E 2017; 96:023205. [PMID: 28950476 DOI: 10.1103/physreve.96.023205] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2017] [Indexed: 06/07/2023]
Abstract
When matter is exposed to a high-intensity x-ray free-electron-laser pulse, the x rays excite inner-shell electrons leading to the ionization of the electrons through various atomic processes and creating high-energy-density plasma, i.e., warm or hot dense matter. The resulting system consists of atoms in various electronic configurations, thermalizing on subpicosecond to picosecond timescales after photoexcitation. We present a simulation study of x-ray-heated solid-density matter. For this we use XMDYN, a Monte Carlo molecular-dynamics-based code with periodic boundary conditions, which allows one to investigate nonequilibrium dynamics. XMDYN is capable of treating systems containing light and heavy atomic species with full electronic configuration space and three-dimensional spatial inhomogeneity. For the validation of our approach we compare for a model system the electron temperatures and the ion charge-state distribution from XMDYN to results for the thermalized system based on the average-atom model implemented in XATOM, an ab initio x-ray atomic physics toolkit extended to include a plasma environment. Further, we also compare the average charge evolution of diamond with the predictions of a Boltzmann continuum approach. We demonstrate that XMDYN results are in good quantitative agreement with the above-mentioned approaches, suggesting that the current implementation of XMDYN is a viable approach to simulate the dynamics of x-ray-driven nonequilibrium dynamics in solids. To illustrate the potential of XMDYN for treating complex systems, we present calculations on the triiodo benzene derivative 5-amino-2,4,6-triiodoisophthalic acid (I3C), a compound of relevance of biomolecular imaging, consisting of heavy and light atomic species.
Collapse
Affiliation(s)
- Malik Muhammad Abdullah
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
- Department of Physics, University of Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
| | - Zoltan Jurek
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Sang-Kil Son
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Robin Santra
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
- Department of Physics, University of Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
| |
Collapse
|
6
|
Ziaja B, Saxena V, Son SK, Medvedev N, Barbrel B, Woloncewicz B, Stransky M. Kinetic Boltzmann approach adapted for modeling highly ionized matter created by x-ray irradiation of a solid. Phys Rev E 2016; 93:053210. [PMID: 27300998 DOI: 10.1103/physreve.93.053210] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2015] [Indexed: 11/07/2022]
Abstract
We report on the kinetic Boltzmann approach adapted for simulations of highly ionized matter created from a solid by its x-ray irradiation. X rays can excite inner-shell electrons, which leads to the creation of deeply lying core holes. Their relaxation, especially in heavier elements, can take complicated paths, leading to a large number of active configurations. Their number can be so large that solving the set of respective evolution equations becomes computationally inefficient and another modeling approach should be used instead. To circumvent this complexity, the commonly used continuum models employ a superconfiguration scheme. Here, we propose an alternative approach which still uses "true" atomic configurations but limits their number by restricting the sample relaxation to the predominant relaxation paths. We test its reliability, performing respective calculations for a bulk material consisting of light atoms and comparing the results with a full calculation including all relaxation paths. Prospective application for heavy elements is discussed.
Collapse
Affiliation(s)
- Beata Ziaja
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany.,Institute of Nuclear Physics, Polish Academy of Sciences, Radzikowskiego 152, 31-342 Kraków, Poland
| | - Vikrant Saxena
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Sang-Kil Son
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Nikita Medvedev
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Benjamin Barbrel
- Center for Intense Lasers and Applications (CELIA), University of Bordeaux 1, 351 Cours de la Liberation, F-33405 Talence, France
| | - Bianca Woloncewicz
- Institute of Experimental Physics, University of Gdansk, ulica Wita Stwosza 57, 80-952 Gdansk, Poland
| | - Michal Stransky
- Department of Radiation and Chemical Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague, Czech Republic
| |
Collapse
|