Showing 1–50 of 217 results for author: Bonitz, M

First principles simulations of dense hydrogen

Authors: Michael Bonitz, Jan Vorberger, Mandy Bethkenhagen, Maximilian Böhme, David Ceperley, Alexey Filinov, Thomas Gawne, Frank Graziani, Gianluca Gregori, Paul Hamann, Stephanie Hansen, Markus Holzmann, S. X. Hu, Hanno Kählert, Valentin Karasiev, Uwe Kleinschmidt, Linda Kordts, Christopher Makait, Burkhard Militzer, Zhandos Moldabekov, Carlo Pierleoni, Martin Preising, Kushal Ramakrishna, Ronald Redmer, Sebastian Schwalbe , et al. (2 additional authors not shown)

Abstract: Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extre… ▽ More Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data. rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange-correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principle PIMC and DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Yuri Lvovich Klimontovich, his theory of fluctuations and its impact on the kinetic theory

Authors: Michael Bonitz, Anatoly Zagorodny

Abstract: Yuri L'vovich Klimontovich (28.09.1924--26.10.2002) was an outstanding theoretical physicist who made major contributions to kinetic theory. On the occasion of his 100th birthday we recall his main scientific achievements. Yuri L'vovich Klimontovich (28.09.1924--26.10.2002) was an outstanding theoretical physicist who made major contributions to kinetic theory. On the occasion of his 100th birthday we recall his main scientific achievements. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Classical and Quantum Theory of Fluctuations for Many-Particle Systems out of Equilibrium

Authors: Erik Schroedter, Michael Bonitz

Abstract: Correlated classical and quantum many-particle systems out of equilibrium are of high interest in many fields, including dense plasmas, correlated solids, and ultracold atoms. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. While for classical systems, in principle, exact simulations are possible via molecular dynam… ▽ More Correlated classical and quantum many-particle systems out of equilibrium are of high interest in many fields, including dense plasmas, correlated solids, and ultracold atoms. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. While for classical systems, in principle, exact simulations are possible via molecular dynamics, this is not the case for quantum systems. Alternatively, one can use many-particle approaches such as hydrodynamics, kinetic theory or nonequilibrium Green functions (NEGF). However, NEGF exhibit a very unfavorable cubic scaling of the CPU time with the number of time steps. An alternative is the G1--G2 scheme [N. Schlünzen et al., Phys. Rev. Lett. \textbf{124}, 076601 (2020)] which allows for NEGF simulations with time linear scaling, however, at the cost of large memory consumption. The reason is the need to store the two-particle correlation function. This problem can be overcome for a number of approximations by reformulating the kinetic equations in terms of fluctuations -- an approach that was developed, for classical systems, by Yu.L. Klimontovich [JETP \textbf{33}, 982 (1957)]. Here we present an overview of his ideas and extend them to quantum systems. In particular, we demonstrate that this quantum fluctuations approach can reproduce the nonequilibrium $GW$ approximation [E. Schroedter \textit{et al.}, Cond. Matt. Phys. \textbf{25}, 23401 (2022)] promising high accuracy at low computational cost which arises from an effective semiclassical stochastic sampling procedure. We also demonstrate how to extend the approach to the two-time exchange-correlation functions and the density response properties. [E. Schroedter \textit{et al.}, Phys. Rev. B \textbf{108}, 205109 (2023)]. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Two-Time Quantum Fluctuations Approach and its Relation to the Bethe--Salpeter Equation

Authors: Erik Schroedter, Michael Bonitz

Abstract: Correlated quantum many-particle systems out of equilibrium are of high interest in many fields, including correlated solids, ultracold atoms or dense plasmas. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. We have recently presented a quantum fluctuations approach which is equivalent to the nonequilibrium $GW$ app… ▽ More Correlated quantum many-particle systems out of equilibrium are of high interest in many fields, including correlated solids, ultracold atoms or dense plasmas. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. We have recently presented a quantum fluctuations approach which is equivalent to the nonequilibrium $GW$ approximation [E. Schroedter \textit{et al.}, Cond. Matt. Phys. \textbf{25}, 23401 (2022)] that promises high accuracy at low computational cost. In a second publication [E. Schroedter \textit{et al.}, Phys. Rev. B \textbf{108}, 205109 (2023)], this approach was extended to the two-time exchange-correlation functions and the density response properties. Here, we analyze the properties of this approach in more detail. We demonstrate that the method is equivalent to the Bethe--Salpeter equation for the two-time exchange-correlation function when the generalized Kadanoff-Baym ansatz with Hartree-Fock propagators is applied. △ Less

Submitted 7 February, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Accelerating Nonequilibrium Green functions simulations: the G1-G2 scheme and beyond

Authors: Michael Bonitz, Jan-Philip Joost, Christopher Makait, Erik Schroedter, Tim Karsberger, Karsten Balzer

Abstract: The theory of Nonequilibrium Green functions (NEGF) has seen a rapid development over the recent three decades. Applications include diverse correlated many-body systems in and out of equilibrium. Very good agreement with experiments and available exact theoretical results could be demonstrated if the proper selfenergy approximations were used. However, full two-time NEGF simulations are computati… ▽ More The theory of Nonequilibrium Green functions (NEGF) has seen a rapid development over the recent three decades. Applications include diverse correlated many-body systems in and out of equilibrium. Very good agreement with experiments and available exact theoretical results could be demonstrated if the proper selfenergy approximations were used. However, full two-time NEGF simulations are computationally costly, as they suffer from a cubic scaling of the computation time with the simulation duration. Recently we have introduced the G1-G2 scheme that exactly reformulates the Kadanoff-Baym ansatz with Hartree-Fock propagators (HF-GKBA) into time-local equations, allowing for a dramatic reduction of the scaling to time-linear scaling [Schluenzen et al., Phys. Rev. Lett. \textbf{124}, 076601 (2020)]. Remarkably, this scaling is achieved quickly, and also for high-level selfenergies, including the nonequilibrium $GW$ and $T$-matrix approximations [Joost et al., Phys. Rev. B \textbf{101}, 245101 (2020)]. Even the dynamically screened ladder approximation is now feasible [Joost et al., Phys. Rev. B \textbf{105}, 165155 (2022)], and also applications to electron-boson systems were demonstrated. Here we present an overview on recent results that were achieved with the G1--G2 scheme. We discuss problems and open questions and present further ideas how to overcome the current limitations of the scheme.We illustrate the G1--G2 scheme by presenting applying it to the excitation dynamics of Hubbard clusters, to optical excitation of graphene, and to charge transfer during stopping of ions by correlated materials. △ Less

Submitted 25 February, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

The equation of state of partially ionized hydrogen and deuterium plasma revisited

Authors: A. V. Filinov, M. Bonitz

Abstract: We present novel first-principle fermionic path integral Monte Carlo (PIMC) simulation results for a dense partially ionized hydrogen (deuterium) plasma, for temperatures in the range $15,000$K $\leq T \leq 400,000$K and densities $7 \cdot 10^{-7}$g/cm$^{3}\leq ρ_H \leq 0.085$ g/cm$^{3}$ ($1.4 \cdot 10^{-6}$g/cm$^{3}\leq ρ_D \leq 0.17$ g/cm$^{3}$), corresponding to $100\geq r_s\geq 2$, where… ▽ More We present novel first-principle fermionic path integral Monte Carlo (PIMC) simulation results for a dense partially ionized hydrogen (deuterium) plasma, for temperatures in the range $15,000$K $\leq T \leq 400,000$K and densities $7 \cdot 10^{-7}$g/cm$^{3}\leq ρ_H \leq 0.085$ g/cm$^{3}$ ($1.4 \cdot 10^{-6}$g/cm$^{3}\leq ρ_D \leq 0.17$ g/cm$^{3}$), corresponding to $100\geq r_s\geq 2$, where $r_s=\bar r/a_B$ is the ratio of the mean interparticle distance to the Bohr radius. These simulations are based on the fermionic propagator PIMC (FP-PIMC) approach in the grand canonical ensemble [A. Filinov \textit{et al.}, Contrib. Plasma Phys. \textbf{61}, e202100112 (2021)] and fully account for correlation and quantum degeneracy and spin effects. For the application to hydrogen and deuterium, we develop a combination of the fourth-order factorization and the pair product ansatz for the density matrix. Moreover, we avoid the fixed node approximation that may lead to uncontrolled errors in restricted PIMC (RPIMC). Our results allow us to critically re-evaluate the accuracy of the RPIMC simulations for hydrogen by Hu \textit{et al.} [Phys. Rev. B \textbf{84}, 224109 (2011)] and of various chemical models. The deviations are generally found to be small, but for the lowest temperature, $T=15,640$~K they reach several percent. We present detailed tables with our first principles results for the pressure and energy isotherms. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Journal ref: Phys. Rev. E 108, 055212 (2023)

Quantum Fluctuations Approach to the Nonequilibrium $GW$-Approximation II: Density Correlations and Dynamic Structure Factor

Authors: Erik Schroedter, Björn Jakob Wurst, Jan-Philip Joost, Michael Bonitz

Abstract: The quantum dynamics of correlated fermionic or bosonic many-body systems following external excitation can be successfully studied using nonequilibrium Green functions (NEGF) or reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle selfenergy or decoupling of the BBGKY-hierarchy, respectively. These approximations are based on Feynman's diagram app… ▽ More The quantum dynamics of correlated fermionic or bosonic many-body systems following external excitation can be successfully studied using nonequilibrium Green functions (NEGF) or reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle selfenergy or decoupling of the BBGKY-hierarchy, respectively. These approximations are based on Feynman's diagram approaches or on cluster expansions into single-particle and correlation operators. In a recent paper [E. Schroedter, J.-P. Joost, and M. Bonitz, Cond. Matt. Phys. \textbf{25}, 23401 (2022)] we have presented a different approach where, instead of equations of motion for the many-particle NEGF (or density operators), equations for the correlation functions of fluctuations are analyzed. In particular, we derived the stochastic GW and polarization approximations that are closely related to the nonequilibrium GW approximation. Here, we extend this approach to the computation of two-time observables depending on the specific ordering of the underlying operators. In particular, we apply this extension to the calculation of the density correlation function and dynamic structure factor of correlated Hubbard clusters in and out of equilbrium. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Prediction of a roton-type feature in warm dense hydrogen

Authors: Paul Hamann, Linda Kordts, Alexey Filinov, Michael Bonitz, Tobias Dornheim, Jan Vorberger

Abstract: In a recent Letter [T. Dornheim \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 255001 (2018)], it was predicted on the basis of \textit{ab initio} quantum Monte Carlo simulations that, in a uniform electron gas, the peak $ω_0$ of the dynamic structure factor $S(q,ω)$ exhibits an unusual non-monotonic wave number dependence, where $dω_0/dq < 0$, at intermediate $q$, under strong coupling condition… ▽ More In a recent Letter [T. Dornheim \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 255001 (2018)], it was predicted on the basis of \textit{ab initio} quantum Monte Carlo simulations that, in a uniform electron gas, the peak $ω_0$ of the dynamic structure factor $S(q,ω)$ exhibits an unusual non-monotonic wave number dependence, where $dω_0/dq < 0$, at intermediate $q$, under strong coupling conditions. This effect was subsequently explained by the pair alignment of electrons %at an intermediate range of wave numbers [T. Dornheim \textit{et al.}, Comm. Phys. \textbf{5}, 304 (2022)]. Here we predict that this non-monotonic dispersion resembling the roton-type behavior known from superfluids should be observable in a dense, partially ionized hydrogen plasma. Based on a combination of path integral Monte Carlo simulations and linear response results for the density response function, we present the approximate range of densities, temperatures and wave numbers and make predictions for possible experimental observations. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Journal ref: Phys. Rev. Research 5, 033039 (2023)

Non-equilibrium correlation dynamics in the one-dimensional Fermi-Hubbard model: A testbed for the two-particle reduced density matrix theory

Authors: Stefan Donsa, Fabian Lackner, Joachim Burgdörfer, Michael Bonitz, Benedikt Kloss, Angel Rubio, Iva Březinová

Abstract: We explore the non-equilibrium dynamics of a one-dimensional Fermi-Hubbard system as a sensitive testbed for the capabilities of the time-dependent two-particle reduced density matrix (TD2RDM) theory to accurately describe time-dependent correlated systems. We follow the time evolution of the out-of-equilibrium finite-size Fermi-Hubbard model initialized by a quench over extended periods of time.… ▽ More We explore the non-equilibrium dynamics of a one-dimensional Fermi-Hubbard system as a sensitive testbed for the capabilities of the time-dependent two-particle reduced density matrix (TD2RDM) theory to accurately describe time-dependent correlated systems. We follow the time evolution of the out-of-equilibrium finite-size Fermi-Hubbard model initialized by a quench over extended periods of time. By comparison with exact calculations for small systems and with matrix product state (MPS) calculations for larger systems but limited to short times, we demonstrate that the TD2RDM theory can accurately account for the non-equilibrium dynamics in the regime from weak to moderately strong inter-particle correlations. We find that the quality of the approximate reconstruction of the three-particle cumulant (or correlation) required for the closure of the equations of motion for the reduced density matrix is key to the accuracy of the numerical TD2RDM results. We identify the size of the dynamically induced three-particle correlations and the amplitude of cross correlations between the two- and three-particle cumulants as critical parameters that control the accuracy of the TD2RDM theory when current state-of-the art reconstruction functionals are employed. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 18 pages, 13 figures

Journal ref: Phys. Rev. Research 5, 033022 (2023)

Time-dependent charged particle stopping in quantum plasmas: testing the G1-G2 scheme for quasi-one-dimensional systems

Authors: C. Makait, F. Borges Fajardo, M. Bonitz

Abstract: Warm dense matter--an exotic, highly compressed state on the boarder between solid and plasma phases is of high current interest, in particular for compact astrophysical objects, high pressure laboratory systems, and inertial confinement fusion. For many applications the interaction of quantum plasmas with energetic particles is crucial. Moreover, often the system is driven far out of equilibrium.… ▽ More Warm dense matter--an exotic, highly compressed state on the boarder between solid and plasma phases is of high current interest, in particular for compact astrophysical objects, high pressure laboratory systems, and inertial confinement fusion. For many applications the interaction of quantum plasmas with energetic particles is crucial. Moreover, often the system is driven far out of equilibrium. In that case, there is high interest in time-dependent simulations to understand the physics, in particular, during thermalization. Recently a novel many-particle technique, the G1--G2 scheme was presented [N. Schlünzen et al., Phys. Rev. Lett. \textbf{124}, 076601 (2020)] which allows for first-principle simulations of the time evolution of interacting quantum systems. Here we apply this scheme to a spatially uniform dense quantum plasma (jellium) and explore its performance. To this end the G1-G2 scheme is transformed into momentum representation, and first results are presented for a quasi-one-dimensional model system. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Dynamic structure factor and excitation spectrum of the one-component plasma: the case of weak to moderate magnetization

Authors: Hanno Kählert, Michael Bonitz

Abstract: Magnetized plasmas are well known to exhibit a rich spectrum of collective modes. Here, we focus on the density modes in dense or cold plasmas, where strong coupling effects alter the mode spectrum known from traditional weakly coupled plasmas. In particular, we study the dynamic structure factor (DSF) of the magnetized one-component plasma with molecular dynamics simulations. Extending our previo… ▽ More Magnetized plasmas are well known to exhibit a rich spectrum of collective modes. Here, we focus on the density modes in dense or cold plasmas, where strong coupling effects alter the mode spectrum known from traditional weakly coupled plasmas. In particular, we study the dynamic structure factor (DSF) of the magnetized one-component plasma with molecular dynamics simulations. Extending our previous results [H.~Kählert and M.~Bonitz, Phys. Rev. Research \textbf{2022}, 4, 013197], it is shown that Bernstein modes can be observed in the weakly magnetized regime, where they are found below the upper hybrid frequency, provided the coupling strength is sufficiently low. We investigate the DSF for a variety of different wave numbers and plasma parameters and show that even small magnetization can give rise to a strong zero-frequency mode perpendicular to the magnetic field and change the dispersion as well as the damping of the upper hybrid mode. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Electronic Density Response of Warm Dense Matter

Authors: Tobias Dornheim, Zhandos A. Moldabekov, Kushal Ramakrishna, Panagiotis Tolias, Andrew D. Baczewski, Dominik Kraus, Thomas R. Preston, David A. Chapman, Maximilian P. Böhme, Tilo Döppner, Frank Graziani, Michael Bonitz, Attila Cangi, Jan Vorberger

Abstract: Matter at extreme temperatures and pressures -- commonly known as warm dense matter (WDM) in the literature -- is ubiquitous throughout our Universe and occurs in a number of astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion, and is realized in the laboratory using different tech… ▽ More Matter at extreme temperatures and pressures -- commonly known as warm dense matter (WDM) in the literature -- is ubiquitous throughout our Universe and occurs in a number of astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion, and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are routinely probed in x-ray Thomson scattering (XRTS) experiments, and, in addition, are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear-response theory as well as nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques including ab initio path integral Monte Carlo (PIMC) simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrarily complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments. △ Less

Submitted 19 December, 2022; v1 submitted 16 December, 2022; originally announced December 2022.

Accelerating Nonequilibrium Green functions simulations with embedding selfenergies

Authors: Karsten Balzer, Niclas Schlünzen, Hannes Ohldag, Jan-Philip Joost, Michael Bonitz

Abstract: Real-time nonequilibrium Green functions (NEGF) have been very successful to simulate the dynamics of correlated many-particle systems far from equilibrium. However, NEGF simulations are computationally expensive since the effort scales cubically with the simulation duration. Recently we have introduced the G1--G2 scheme that allows for a dramatic reduction to time-linear scaling [Schlünzen, Phys.… ▽ More Real-time nonequilibrium Green functions (NEGF) have been very successful to simulate the dynamics of correlated many-particle systems far from equilibrium. However, NEGF simulations are computationally expensive since the effort scales cubically with the simulation duration. Recently we have introduced the G1--G2 scheme that allows for a dramatic reduction to time-linear scaling [Schlünzen, Phys. Rev. Lett. 124, 076601 (2020); Joost et al., Phys. Rev. B 101, 245101 (2020)]. Here we tackle another problem: the rapid growth of the computational effort with the system size. In many situations where the system of interest is coupled to a bath, to electric contacts or similar macroscopic systems for which a microscopic resolution of the electronic properties is not necessary, efficient simplifications are possible. This is achieved by the introduction of an embedding selfenergy -- a concept that has been successful in standard NEGF simulations. Here, we demonstrate how the embedding concept can be introduced into the G1--G2 scheme, allowing us to drastically accelerate NEGF embedding simulations. The approach is compatible with all advanced selfenergies that can be represented by the G1--G2 scheme [as described in Joost et al., Phys. Rev. B 105, 165155 (2022)] and retains the memory-less structure of the equations and their time linear scaling. As a numerical illustration we investigate the charge transfer between a Hubbard nanocluster and an additional site which is of relevance for the neutralization of ions in matter. △ Less

Submitted 7 April, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Quantum fluctuations approach to the nonequilibrium $GW$ approximation

Authors: Erik Schroedter, Jan-Philip Joost, Michael Bonitz

Abstract: The quantum dynamics of fermionic or bosonic many-body systems following external excitation can be successfully studied using two-time nonequilibrium Green's functions (NEGF) or single-time reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle self-energy or decoupling of the BBGKY hierarchy. These approximations are based on Feynman's diagram appr… ▽ More The quantum dynamics of fermionic or bosonic many-body systems following external excitation can be successfully studied using two-time nonequilibrium Green's functions (NEGF) or single-time reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle self-energy or decoupling of the BBGKY hierarchy. These approximations are based on Feynman's diagram approaches or on cluster expansions into single-particle and correlation operators. Here, we develop a different approach where, instead of equations of motion for the many-particle NEGF (or density operators), single-time equations for the correlation functions of fluctuations are analyzed. We present a derivation of the first two equations of the alternative hierarchy of fluctuations and discuss possible decoupling approximations. In particular, we derive the polarization approximation (PA) which is shown to be equivalent to the single-time version [following by applying the generalized Kadanoff-Baym ansatz (GKBA)] of the nonequilibrium $GW$ approximation with exchange effects of NEGF theory, for weak coupling. The main advantage of the quantum fluctuations approach is that the standard ensemble average can be replaced by a semiclassical average over different initial realizations, as was demonstrated before by Lacroix and co-workers [see e.g. D. Lacroix et al., Phys. Rev. B, 2014, 90, 125112]. Here, we introduce the stochastic $GW$ (SGW) approximation and the stochastic polarization approximation (SPA) which are demonstrated to be equivalent to the single-time $GW$ approximation without and with exchange, respectively, in the weak coupling limit. Our numerical tests confirm that our approach has the same favorable linear scaling with the computation time as the recently developed G1-G2 scheme [Schluenzen et al., Phys. Rev. Lett., 2020, 124, 076601]. △ Less

Submitted 9 July, 2022; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: 36 pages, 14 figures

Journal ref: Condensed Matter Physics, 2022, vol. 25, No. 2, 23401

Electronic pair alignment and roton feature in the warm dense electron gas

Authors: Tobias Dornheim, Zhandos Moldabekov, Jan Vorberger, Hanno Kählert, Michael Bonitz

Abstract: The study of matter under extreme densities and temperatures as they occur e.g. in astrophysical objects and nuclear fusion applications has emerged as one of the most active frontiers in physics, material science, and related disciplines. In this context, a key quantity is given by the dynamic structure factor $S(\mathbf{q},ω)$, which is probed in scattering experiments -- the most widely used me… ▽ More The study of matter under extreme densities and temperatures as they occur e.g. in astrophysical objects and nuclear fusion applications has emerged as one of the most active frontiers in physics, material science, and related disciplines. In this context, a key quantity is given by the dynamic structure factor $S(\mathbf{q},ω)$, which is probed in scattering experiments -- the most widely used method of diagnostics at these extreme conditions. In addition to its crucial importance for the study of warm dense matter, the modelling of such dynamic properties of correlated quantum many-body systems constitutes one of the most fundamental theoretical challenges of our time. Here we report a hitherto unexplained \emph{roton feature} in $S(\mathbf{q},ω)$ of the warm dense electron gas, and introduce a microscopic explanation in terms of a new \emph{electronic pair alignment} model. This new paradigm will be highly important for the understanding of warm dense matter, and has a direct impact on the interpretation of scattering experiments. Moreover, we expect our results to give unprecedented insights into the dynamics of a number of correlated quantum many-body systems such as ultracold helium, dipolar supersolids, and bilayer heterostructures. △ Less

Submitted 23 March, 2022; originally announced March 2022.

The dynamically screened ladder approximation: Simultaneous treatment of strong electronic correlations and dynamical screening out of equilibrium

Authors: Jan-Philip Joost, Niclas Schlünzen, Hannes Ohldag, Michael Bonitz, Fabian Lackner, Iva Březinová

Abstract: Dynamical screening is a key property of charged many-particle systems. Its theoretical description is based on the $GW$ approximation that is extensively applied for ground-state and equilibrium situations but also for systems driven out of equilibrium. The main limitation of the $GW$ approximation is the neglect of strong electronic correlation effects that are important in many materials as wel… ▽ More Dynamical screening is a key property of charged many-particle systems. Its theoretical description is based on the $GW$ approximation that is extensively applied for ground-state and equilibrium situations but also for systems driven out of equilibrium. The main limitation of the $GW$ approximation is the neglect of strong electronic correlation effects that are important in many materials as well as in dense plasmas. Here we derive the dynamically screened ladder (DSL) approximation that selfconsistently includes, in addition to the $GW$ diagrams, also particle--particle and particle--hole $T$-matrix diagrams. The derivation is based on reduced-density-operator theory and the result is equivalent to the recently presented G1--G2 scheme [Schlünzen \textit{et al.}, Phys. Rev. Lett. \textbf{124}, 076601 (2020); Joost \textit{et al.}, Phys. Rev. B \textbf{101}, 245101 (2020)]. We perform extensive time-dependent DSL simulations for finite Hubbard clusters and present tests against exact results that confirm excellent accuracy as well as total energy conservation of the approximation. At strong coupling and for long simulation durations, instabilities are observed. These problems are solved by enforcing contraction consistency and applying a purification approach. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Journal ref: Phys. Rev. B 105 (16), 165155 (2022)

Dynamic structure factor of the magnetized one-component plasma: crossover from weak to strong coupling

Authors: Hanno Kählert, Michael Bonitz

Abstract: Plasmas in strong magnetic fields have been mainly studied in two distinct limiting cases--that of weak and strong nonideality with very different physical properties. While the former is well described by the familiar theory of Braginskii, the latter regime is closer to the behavior of a Coulomb liquid. Here we study in detail the transition between both regimes. We focus on the evolution of the… ▽ More Plasmas in strong magnetic fields have been mainly studied in two distinct limiting cases--that of weak and strong nonideality with very different physical properties. While the former is well described by the familiar theory of Braginskii, the latter regime is closer to the behavior of a Coulomb liquid. Here we study in detail the transition between both regimes. We focus on the evolution of the dynamic structure factor of the magnetized one-component plasma from weak to strong coupling, which is studied with first-principle molecular dynamics simulations. The simulations show the vanishing of Bernstein modes and the emergence of higher harmonics of the upper hybrid mode across the magnetic field, a redistribution of spectral power between the two main collective modes under oblique angles, and a suppression of plasmon damping along the magnetic field. Comparison with results from various models, including the random phase approximation, a Mermin-type dielectric function, and the Quasi-Localized Charge Approximation show that none of the theories is capable of reproducing the crossover that occurs when the coupling parameter is on the order of unity. The findings are relevant to the scattering spectra, stopping power, and transport coefficients of correlated magnetized plasmas. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Nonlinear electronic density response of the warm dense electron gas: multiple perturbations and mode coupling

Authors: Tobias Dornheim, Jan Vorberger, Zhandos Moldabekov, Michael Bonitz

Abstract: We present extensive new ab initio path integral Monte Carlo (PIMC) results for an electron gas at warm dense matter conditions that is subject to multiple harmonic perturbations. In addition to the previously investigated nonlinear effects at the original wave number [Dornheim \emph{et al.}, PRL \textbf{125}, 085001 (2020)] and the excitation of higher harmonics [Dornheim \emph{et al.}, PRR \text… ▽ More We present extensive new ab initio path integral Monte Carlo (PIMC) results for an electron gas at warm dense matter conditions that is subject to multiple harmonic perturbations. In addition to the previously investigated nonlinear effects at the original wave number [Dornheim \emph{et al.}, PRL \textbf{125}, 085001 (2020)] and the excitation of higher harmonics [Dornheim \emph{et al.}, PRR \textbf{3}, 033231 (2021)], the presence of multiple external potentials leads to mode-coupling effects, which constitute the dominant nonlinear effect and lead to a substantially more complicated density response compared to linear response theory. One possibility to estimate mode-coupling effects from a PIMC simulation of the unperturbed system is given in terms of generalized imaginary-time correlation functions that have been recently introduced by Dornheim \emph{et al.}~[JCP \textbf{155}, 054110 (2021)]. In addition, we extend our previous analytical theory of the nonlinear density response of the electron gas in terms of the static local field correction [Dornheim \emph{et al.}, PRL \textbf{125}, 235001 (2020)], which allows for a highly accurate description of the PIMC results with negligible computational cost. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Doublon production in correlated materials by multiple ion impacts

Authors: Lotte Borkowski, Niclas Schlünzen, Jan Philip Joost, Franziska Reiser, Michael Bonitz

Abstract: In a recent Letter [Balzer \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 267602 (2018)] it was demonstrated that ions impacting a correlated graphene cluster can excite strongly nonequilibrium states. In particular, this can lead to an enhanced population of bound pairs of electrons with opposite spin -- doublons -- where the doublon number can be increased via multiple ion impacts. These predic… ▽ More In a recent Letter [Balzer \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 267602 (2018)] it was demonstrated that ions impacting a correlated graphene cluster can excite strongly nonequilibrium states. In particular, this can lead to an enhanced population of bound pairs of electrons with opposite spin -- doublons -- where the doublon number can be increased via multiple ion impacts. These predictions were made based on nonequilibrium Green functions (NEGF) simulations allowing for a time-dependent non-perturbative study of the energy loss of charged particles penetrating a strongly correlated system. Here we extend these simulations to larger clusters and longer simulation times, utilizing the recently developed G1--G2 scheme [Schünzen \textit{et al.}, Phys. Rev. Lett. \textbf{124}, 076601 (2020)] which allows for a dramatic speedup of NEGF simulations. Furthermore, we investigate the dependence of the energy and doublon number on the time interval between ion impacts and on the impact point. △ Less

Submitted 13 October, 2021; originally announced October 2021.

Shock Physics in Warm Dense Matter--a quantum hydrodynamics perspective

Authors: F. Graziani, Z. Moldabekov, B. Olson, M. Bonitz

Abstract: Warm dense matter (WDM)--an exotic, highly compressed state of matter between solid and plasma phases is of high current interest, in particular for astrophysics and inertial confinement fusion. For the latter, in particular the propagation of compression shocks is crucial. The main unknown in the shock propagation in WDM is the behavior of the electrons since they are governed by correlations, qu… ▽ More Warm dense matter (WDM)--an exotic, highly compressed state of matter between solid and plasma phases is of high current interest, in particular for astrophysics and inertial confinement fusion. For the latter, in particular the propagation of compression shocks is crucial. The main unknown in the shock propagation in WDM is the behavior of the electrons since they are governed by correlations, quantum and spin effects that need to be accounted for simultaneously. Here we describe the shock dynamics of the warm dense electron gas using a quantum hydrodynamic model. From the numerical hydrodynamic simulations we observe that the quantum Bohm pressure induces shear force which weakens the formation and strength of the shock. This is confirmed by the theoretical analysis of the early stage of the shock formation. Our theoretical and numerical analysis allows us to identify characteristic dimensionless shock propagation parameters at which the effect of the Bohm force is important. △ Less

Submitted 23 November, 2021; v1 submitted 19 September, 2021; originally announced September 2021.

Journal ref: Contributions to Plasma Physics 2022, e202100170

Density Response of the Warm Dense Electron Gas beyond Linear Response Theory: Excitation of Harmonics

Authors: Tobias Dornheim, Maximilian Böhme, Zhandos A. Moldabekov, Jan Vorberger, Michael Bonitz

Abstract: In a recent Letter, Dornheim et al. [PRL 125, 085001 (2020)] have investigated the nonlinear density response of the uniform electron gas in the warm dense matter regime. More specifically, they have studied the cubic response function at the first harmonic, which cannot be neglected in many situations of experimental relevance. In this work, we go one step further and study the full spectrum of e… ▽ More In a recent Letter, Dornheim et al. [PRL 125, 085001 (2020)] have investigated the nonlinear density response of the uniform electron gas in the warm dense matter regime. More specifically, they have studied the cubic response function at the first harmonic, which cannot be neglected in many situations of experimental relevance. In this work, we go one step further and study the full spectrum of excitations at the higher harmonics of the original perturbation based on extensive new ab initio path integral Monte Carlo (PIMC) simulations. We find that the dominant contribution to the density response beyond linear response theory is given by the quadratic response function at the second harmonic in the moderately nonlinear regime. Furthermore, we show that the nonlinear density response is highly sensitive to exchange-correlation effects, which makes it a potentially valuable new tool of diagnostics. To this end, we present a new theoretical description of the nonlinear electronic density response based on the recent effective static approximation to the local field correction [PRL 125, 235001 (2020)], which accurately reproduces our PIMC data with negligible computational cost. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. Research 3, 033231 (2021)

Towards a Quantum Fluid Theory of Correlated Many-Fermion Systems from First Principles

Authors: Zh. A. Moldabekov, T. Dornheim, G. Gregori, F. Graziani, M. Bonitz, A. Cangi

Abstract: Correlated many-fermion systems emerge in a broad range of phenomena in warm dense matter, plasmonics, and ultracold atoms. Quantum hydrodynamics (QHD) complements common first-principles methods for many-fermion systems and enables simulations at larger length and longer time scales. While the quantum Bohm potential is central to QHD, we illustrate its failure for strong perturbations. We extend… ▽ More Correlated many-fermion systems emerge in a broad range of phenomena in warm dense matter, plasmonics, and ultracold atoms. Quantum hydrodynamics (QHD) complements common first-principles methods for many-fermion systems and enables simulations at larger length and longer time scales. While the quantum Bohm potential is central to QHD, we illustrate its failure for strong perturbations. We extend QHD to this regime by utilizing the many-fermion quantum Bohm potential. This opens up the path to more accurate simulations in strongly perturbed warm dense matter, inhomogeneous quantum plasmas, and on nano-structure surfaces at scales unattainable with first-principles algorithms. The many-fermion quantum Bohm potential might also have important astrophysical applications in developing conformal-invariant cosmologies. △ Less

Submitted 22 April, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Journal ref: SciPost Phys. 12, 062 (2022)

Neutralization dynamics of slow highly charged ions passing through graphene nanoflakes--an embedding self-energy approach

Authors: Karsten Balzer, Michael Bonitz

Abstract: We study the time-dependent neutralization of a slow highly charged ion that penetrates a hexagonal hollow-centred graphene nanoflake. To compute the ultrafast charge transfer dynamics, we apply an effective Hubbard nanocluster model and use the method of nonequilibrium Green functions (NEGF) in conjunction with an embedding self-energy scheme which allows one to follow the temporal changes of the… ▽ More We study the time-dependent neutralization of a slow highly charged ion that penetrates a hexagonal hollow-centred graphene nanoflake. To compute the ultrafast charge transfer dynamics, we apply an effective Hubbard nanocluster model and use the method of nonequilibrium Green functions (NEGF) in conjunction with an embedding self-energy scheme which allows one to follow the temporal changes of the number of electrons in the nanoflake. We perform extensive simulations of the charge transfer dynamics for a broad range of ion charge states and impact velocities. The results are used to put forward a simple semi-analytical model of the neutralization dynamics that is in very good agreement with transmission experiments, in which highly charged xenon ions pass through sheets of single-layer graphene. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Journal ref: Contrib. Plasma Phys. 62(2), e202100041 (2021)

Momentum distribution function and short-range correlations of the warm dense electron gas -- ab initio quantum Monte Carlo results

Authors: Kai Hunger, Tim Schoof, Tobias Dornheim, Michael Bonitz, Alexey Filinov

Abstract: In a classical plasma the momentum distribution, $n(k)$, decays exponentially, for large $k$, and the same is observed for an ideal Fermi gas. However, when quantum and correlation effects are relevant simultaneously, an algebraic decay, $n_\infty(k)\sim k^{-8}$ has been predicted. This is of relevance for cross sections and threshold processes in dense plasmas that depend on the number of energet… ▽ More In a classical plasma the momentum distribution, $n(k)$, decays exponentially, for large $k$, and the same is observed for an ideal Fermi gas. However, when quantum and correlation effects are relevant simultaneously, an algebraic decay, $n_\infty(k)\sim k^{-8}$ has been predicted. This is of relevance for cross sections and threshold processes in dense plasmas that depend on the number of energetic particles. Here we present extensive \textit{ab initio} results for the momentum distribution of the nonideal uniform electron gas at warm dense matter conditions. Our results are based on first principle fermionic path integral Monte Carlo (CPIMC) simulations and clearly confirm the $k^{-8}$ asymptotic. This asymptotic behavior is directly linked to short-range correlations which are analyzed via the on-top pair distribution function (on-top PDF), i.e. the PDF of electrons with opposite spin. We present extensive results for the density and temperature dependence of the on-top PDF and for the momentum distribution in the entire momentum range. △ Less

Submitted 8 March, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Journal ref: Phys. Rev. E 103, 053204 (2021)

Screening of a test charge in a free-electron gas at warm dense matter and dense non-ideal plasma conditions

Authors: Zhandos Moldabekov, Tobias Dornheim, Michael Bonitz

Abstract: The screening of a test charge by partially degenerate non-ideal free electrons at conditions related to warm dense matter and dense plasmas is investigated using linear response theory and the local field correction based on ab inito Quantum Monte-Carlo simulations data. The analysis of the obtained results is performed by comparing to the random phase approximation and the Singwi-Tosi-Land-Sjöla… ▽ More The screening of a test charge by partially degenerate non-ideal free electrons at conditions related to warm dense matter and dense plasmas is investigated using linear response theory and the local field correction based on ab inito Quantum Monte-Carlo simulations data. The analysis of the obtained results is performed by comparing to the random phase approximation and the Singwi-Tosi-Land-Sjölander approximation. The applicability of the long-wavelength approximation for the description of screening is investigated. The impact of electronic exchange-correlations effects on structural properties and the applicability of the screened potential from linear response theory for the simulation of the dynamics of ions are discussed. △ Less

Submitted 19 September, 2020; originally announced September 2020.

Journal ref: Contributions to Plasma Physics 61, e202000176 (2021)

Finite temperature density functional theory investigation to the nonequilibrium transient warm dense state created by laser excitation

Authors: Hengyu Zhang, Shen Zhang, Dongdong Kang, Jiayu Dai, M. Bonitz

Abstract: We present a finite-temperature density functional theory investigation of the nonequilibrium transient electronic structure of warm dense Li, Al, Cu, and Au created by laser excitation. Photons excite electrons either from the inner shell orbitals or from the valence bands according to the photon energy, and give rise to isochoric heating of the sample. Localized states related to the 3d orbital… ▽ More We present a finite-temperature density functional theory investigation of the nonequilibrium transient electronic structure of warm dense Li, Al, Cu, and Au created by laser excitation. Photons excite electrons either from the inner shell orbitals or from the valence bands according to the photon energy, and give rise to isochoric heating of the sample. Localized states related to the 3d orbital are observed for Cu when the hole lies in the inner shell 3s orbital. The electrical conductivity for these materials at nonequilibrium states is calculated using the Kubo-Greenwood formula. The change of the electrical conductivity, compared to the equilibrium state, is different for the case of holes in inner shell orbitals or the valence band. This is attributed to the competition of two factors: the shift of the orbital energies due to reduced screening of core electrons, and the increase of chemical potential due to the excitation of electrons. The finite temperature effect of both the electrons and the ions on the electrical conductivity is discussed in detail. This work is helpful to better understand the physics of laser excitation experiments of warm dense matter. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 111 pages, 9 figures

Journal ref: Phys. Rev. E 103: 013210 (2021)

Ab initio results for the plasmon dispersion and damping of the warm dense electron gas

Authors: Paul Hamann, Jan Vorberger, Tobias Dornheim, Zhandos Moldabekov, Michael Bonitz

Abstract: Warm dense matter (WDM) is an exotic state on the border between condensed matter and dense plasmas. Important occurrences of WDM include dense astrophysical objects, matter in the core of our Earth, as well as matter produced in strong compression experiments. As of late, x-ray Thomson scattering has become an advanced tool to diagnose WDM. The interpretation of the data requires model input for… ▽ More Warm dense matter (WDM) is an exotic state on the border between condensed matter and dense plasmas. Important occurrences of WDM include dense astrophysical objects, matter in the core of our Earth, as well as matter produced in strong compression experiments. As of late, x-ray Thomson scattering has become an advanced tool to diagnose WDM. The interpretation of the data requires model input for the dynamic structure factor $S(q,ω)$ and the plasmon dispersion $ω(q)$. Recently the first \textit{ab initio} results for $S(q,ω)$ of the homogeneous warm dense electron gas were obtained from path integral Monte Carlo simulations, [Dornheim \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 255001 (2018)]. Here, we analyse the effects of correlations and finite temperature on the dynamic dielectric function and the plasmon dispersion. Our results for the plasmon dispersion and damping differ significantly from the random phase approximation and from earlier models of the correlated electron gas. Moreover, we show when commonly used weak damping approximations break down and how the method of complex zeros of the dielectric function can solve this problem for WDM conditions. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Dynamic properties of the warm dense electron gas: an ab initio path integral Monte Carlo approach

Authors: Paul Hamann, Tobias Dornheim, Jan Vorberger, Zhandos A. Moldabekov, Michael Bonitz

Abstract: There is growing interest in warm dense matter (WDM) -- an exotic state on the border between condensed matter and plasmas. Due to the simultaneous importance of quantum and correlation effects WDM is complicated to treat theoretically. A key role has been played by \textit{ab initio} path integral Monte Carlo (PIMC) simulations, and recently extensive results for thermodynamic quantities have bee… ▽ More There is growing interest in warm dense matter (WDM) -- an exotic state on the border between condensed matter and plasmas. Due to the simultaneous importance of quantum and correlation effects WDM is complicated to treat theoretically. A key role has been played by \textit{ab initio} path integral Monte Carlo (PIMC) simulations, and recently extensive results for thermodynamic quantities have been obtained. The first extension of PIMC simulations to the dynamic structure factor of the uniform electron gas were reported by Dornheim \textit{et al.} [Phys. Rev. Lett. \textbf{121}, 255001 (2018)]. This was based on an accurate reconstruction of the dynamic local field correction. Here we extend this concept to other dynamical quantities of the warm dense electron gas including the dynamic susceptibility, the dielectric function and the conductivity. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. B 102, 125150 (2020)

Restricted configuration path integral Monte Carlo

Authors: A. Yilmaz, K. Hunger, T. Dornheim, S. Groth, M. Bonitz

Abstract: Quantum Monte Carlo belongs to the most accurate simulation techniques for quantum many-particle systems. However, for fermions, these simulations are hampered by the sign problem that prohibits simulations in the regime of strong degeneracy. The situation changed with the development of configuration path integral Monte Carlo (CPIMC) by Schoof \textit{et al.} [T. Schoof \textit{et al.}, Contrib.… ▽ More Quantum Monte Carlo belongs to the most accurate simulation techniques for quantum many-particle systems. However, for fermions, these simulations are hampered by the sign problem that prohibits simulations in the regime of strong degeneracy. The situation changed with the development of configuration path integral Monte Carlo (CPIMC) by Schoof \textit{et al.} [T. Schoof \textit{et al.}, Contrib. Plasma Phys. \textbf{51}, 687 (2011)] that allowed for the first \textit{ab initio} simulations for dense quantum plasmas. CPIMC also has a sign problem that occurs when the density is lowered, i.e. in a parameter range that is complementary to traditional QMC formulated in coordinate space. Thus, CPIMC simulations for the warm dense electron gas are limited to small values of the Brueckner parameter -- the ratio of the interparticle distance to the Bohr radius -- $r_s=\bar{r}/a_B \lesssim 1$. In order to reach the regime of stronger coupling (lower density) with CPIMC, here we investigate additional restrictions on the Monte Carlo procedure. In particular, we introduce two different versions of "restricted CPIMC" where certain sign changing Monte Carlo updates are being omitted. Interestingly, one of the methods (RCPIMC) has no sign problem at all, but it is less accurate than RCPIMC+ which neglects only a smaller class of the Monte Carlo steps. Here we report extensive simulations for the ferromagnetic uniform electron gas with which we investigate the properties and accuracy of RCPIMC and RCPIMC+. Further, we establish the parameter range in the density-temperature plane where these simulations are both feasible and accurate. The conclusion is that RCPIMC and RCPIMC+ work best at temperatures in the range of $Θ\sim 0.1\dots 0.5$ allowing to reach density parameters up to $r_s \sim 3\dots 5$, thereby partially filling a gap left open by existing \textit{ab initio} QMC methods. △ Less

Submitted 24 July, 2020; originally announced July 2020.

Löwdin's symmetry dilemma within Green functions theory for the one-dimensional Hubbard model

Authors: J. -P. Joost, N. Schlünzen, S. Hese, M. Bonitz, C. Verdozzi, P. Schmitteckert, M. Hopjan

Abstract: The energy gap of correlated Hubbard clusters is well studied for one-dimensional systems using analytical methods and density-matrix-renormalization-group (DMRG) simulations. Beyond 1D, however, exact results are available only for small systems by quantum Monte Carlo. For this reason and, due to the problems of DMRG in simulating 2D and 3D systems, alternative methods such as Green functions com… ▽ More The energy gap of correlated Hubbard clusters is well studied for one-dimensional systems using analytical methods and density-matrix-renormalization-group (DMRG) simulations. Beyond 1D, however, exact results are available only for small systems by quantum Monte Carlo. For this reason and, due to the problems of DMRG in simulating 2D and 3D systems, alternative methods such as Green functions combined with many-body approximations (GFMBA), that do not have this restriction, are highly important. However, it has remained open whether the approximate character of GFMBA simulations prevents the computation of the Hubbard gap. Here we present new GFMBA results that demonstrate that GFMBA simulations are capable of producing reliable data for the gap which agrees well with the DMRG benchmarks in 1D. An interesting observation is that the accuracy of the gap can be significantly increased when the simulations give up certain symmetry restriction of the exact system, such as spin symmetry and spatial homogeneity. This is seen as manifestation and generalization of the "symmetry dilemma" introduced by Löwdin for Hartree--Fock wave function calculations. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Journal ref: Contributions to Plasma Physics 61, e202000220 (2021)

Nonlinear Electronic Density Response in Warm Dense Matter

Authors: Tobias Dornheim, Jan Vorberger, Michael Bonitz

Abstract: Warm dense matter (WDM)---an extreme state with high temperatures and densities that occurs e.g. in astrophysical objects---constitutes one of the most active fields in plasma physics and materials science. These conditions can be realized in the lab by shock compression or laser excitation, and the most accurate experimental diagnostics is achieved with lasers and free electron lasers which is th… ▽ More Warm dense matter (WDM)---an extreme state with high temperatures and densities that occurs e.g. in astrophysical objects---constitutes one of the most active fields in plasma physics and materials science. These conditions can be realized in the lab by shock compression or laser excitation, and the most accurate experimental diagnostics is achieved with lasers and free electron lasers which is theoretically modeled using linear response theory. Here, we present first \textit{ab initio} path integral Monte Carlo results for the nonlinear density response of correlated electrons in WDM and show that for many situations of experimental relevance nonlinear effects cannot be neglected. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. Lett. 125, 085001 (2020)

The G1--G2 Scheme: Dramatic Acceleration of Nonequilibrium Green Functions Simulations Within the Hartree--Fock-GKBA

Authors: J. -P. Joost, N. Schlünzen, M. Bonitz

Abstract: The time evolution in quantum many-body systems after external excitations is attracting high interest in many fields. The theoretical modeling of these processes is challenging, and the only rigorous quantum-dynamics approach that can treat correlated fermions in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their… ▽ More The time evolution in quantum many-body systems after external excitations is attracting high interest in many fields. The theoretical modeling of these processes is challenging, and the only rigorous quantum-dynamics approach that can treat correlated fermions in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their $T^3$-scaling with the simulation duration $T$. Recently, $T^2$-scaling was achieved with the generalized Kadanoff--Baym ansatz (GKBA), for the second-order Born (SOA) selfenergy, which has substantially extended the scope of NEGF simulations. In a recent Letter [Schlünzen \textit{et al.}, Phys. Rev. Lett. \textbf{124}, 076601 (2020)] we demonstrated that GKBA-NEGF simulations can be efficiently mapped onto coupled time-local equations for the single-particle and two-particle Green functions on the time diagonal, hence the method has been called G1--G2 scheme. This allows one to perform the same simulations with order $T^1$-scaling, both for SOA and $GW$ selfenergies giving rise to a dramatic speedup. Here we present more details on the G1--G2 scheme, including derivations of the basic equations including results for a general basis, for Hubbard systems and for jellium. Also, we demonstrate how to incorporate initial correlations into the G1--G2 scheme. Further, the derivations are extended to a broader class of selfenergies, including the $T$ matrix in the particle--particle and particle--hole channels, and the dynamically screened-ladder approximation. Finally, we demonstrate that, for all selfenergies, the CPU time scaling of the G1--G2 scheme with the basis dimension, $N_b$, can be improved compared to our first report: the overhead compared to the original GKBA, is not more than an additional factor $N_b$. △ Less

Submitted 28 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. B 101, 245101 (2020)

Energy loss and friction characteristics of electrons at warm dense matter and non-ideal dense plasma conditions

Authors: Zh. A. Moldabekov, T. Dornheim, M. Bonitz, T. S. Ramazanov

Abstract: We investigate the energy loss characteristics of warm dense matter (WDM) and dense plasmas concentrating on the influence of electronic correlations. The basis for our analysis is a recently developed ab initio Quantum Monte-Carlo (QMC) based machine-learning representation of the static local field correction (LFC) [Dornheim et al., J. Chem. Phys. 151, 194104 (2019)], which provides an accurate… ▽ More We investigate the energy loss characteristics of warm dense matter (WDM) and dense plasmas concentrating on the influence of electronic correlations. The basis for our analysis is a recently developed ab initio Quantum Monte-Carlo (QMC) based machine-learning representation of the static local field correction (LFC) [Dornheim et al., J. Chem. Phys. 151, 194104 (2019)], which provides an accurate description of the dynamical density response function of the electron gas at the considered parameters. We focus on the polarization-induced stopping power due to free electrons, the friction function, and the straggling rate. In addition, we compute the friction coefficient which constitutes a key quantity for the adequate Langevin dynamics simulation of ions. Considering typical experimental WDM parameters with partially degenerate electrons, we find that the friction coefficient is of the order of $γ/ω_{pi}=0.01$, where $ω_{pi}$ is the ionic plasma frequency. This analysis is performed by comparing QMC based data to results from the random phase approximation (RPA), the Mermin dielectric function, and the Singwi-Tosi-Land-Sjölander (STLS) approximation. It is revealed that the widely used relaxation time approximation (Mermin dielectric function) has severe limitations regarding the description of the energy loss properties of correlated partially degenerate electrons. Moreover, by comparing QMC based data with the results obtained using STLS, we find that energy loss properties are not sensitive to the inaccuracy of the static LFC at large wave numbers $k/k_{F}>2$ (with $k_F$ being the usual Fermi wave number), but that a correct description of the static LFC at $k/k_{F}\lesssim 1.5$ is important. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. E 101, 053203 (2020)

Ab initio simulation of warm dense matter

Authors: M. Bonitz, T. Dornheim, Zh. A. Moldabekov, S. Zhang, P. Hamann, H. Kählert, A. Filinov, K. Ramakrishna, J. Vorberger

Abstract: Warm dense matter (WDM) -- an exotic state of highly compressed matter -- has attracted high interest in recent years in astrophysics and for dense laboratory systems. At the same time, this state is extremely difficult to treat theoretically. This is due to the simultaneous appearance of quantum degeneracy, Coulomb correlations and thermal effects, as well as the overlap of plasma and condensed p… ▽ More Warm dense matter (WDM) -- an exotic state of highly compressed matter -- has attracted high interest in recent years in astrophysics and for dense laboratory systems. At the same time, this state is extremely difficult to treat theoretically. This is due to the simultaneous appearance of quantum degeneracy, Coulomb correlations and thermal effects, as well as the overlap of plasma and condensed phases. Recent breakthroughs are due to the successful application of density functional theory (DFT) methods which, however, often lack the necessary accuracy and predictive capability for WDM applications. The situation has changed with the availability of the first \textit{ab initio} data for the exchange-correlation free energy of the warm dense uniform electron gas (UEG) that were obtained by quantum Monte Carlo (QMC) simulations, for recent reviews, see Dornheim \textit{et al.}, Phys. Plasmas \textbf{24}, 056303 (2017) and Phys. Rep. \textbf{744}, 1-86 (2018). In the present article we review recent further progress in QMC simulations of the warm dense UEG: namely, \textit{ab initio} results for the static local field correction $G(q)$ and for the dynamic structure factor $S(q,ω)$. These data are of key relevance for the comparison with x-ray scattering experiments at free electron laser facilities and for the improvement of theoretical models. In the second part of this paper we discuss simulations of WDM out of equilibrium. The theoretical approaches include Born-Oppenheimer molecular dynamics, quantum kinetic theory, time-dependent DFT and hydrodynamics. Here we analyze strengths and limitations of these methods and argue that progress in WDM simulations will require a suitable combination of all methods. A particular role might be played by quantum hydrodynamics, and we concentrate on problems, recent progress, and possible improvements of this method. △ Less

Submitted 28 February, 2020; v1 submitted 20 December, 2019; originally announced December 2019.

Correlated Topological States in Graphene Nanoribbon Heterostructures

Authors: Jan-Philip Joost, Antti-Pekka Jauho, Michael Bonitz

Abstract: Finite graphene nanoribbon (GNR) heterostructures host intriguing topological in-gap states (Rizzo, D. J. et al.~\textit{Nature} \textbf{2018}, \textit{560}, 204]). These states may be localized either at the bulk edges, or at the ends of the structure. Here we show that correlation effects (not included in previous density functional simulations) play a key role in these systems: they result in i… ▽ More Finite graphene nanoribbon (GNR) heterostructures host intriguing topological in-gap states (Rizzo, D. J. et al.~\textit{Nature} \textbf{2018}, \textit{560}, 204]). These states may be localized either at the bulk edges, or at the ends of the structure. Here we show that correlation effects (not included in previous density functional simulations) play a key role in these systems: they result in increased magnetic moments at the ribbon edges accompanied by a significant energy renormalization of the topological end states -- even in the presence of a metallic substrate. Our computed results are in excellent agreement with the experiments. Furthermore, we discover a striking, novel mechanism that causes an energy splitting of the non-zero-energy topological end states for a weakly screened system. We predict that similar effects should be observable in other GNR heterostructures as well. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Journal ref: Nano Letters 19, 9045-9050 (2019)

Achieving the Ultimate Scaling Limit for Nonequilibrium Green Functions Simulations

Authors: Niclas Schluenzen, Jan-Philip Joost, Michael Bonitz

Abstract: The dynamics of strongly correlated fermions following an external excitation reveals extremely rich collective quantum effects. Examples are fermionic atoms in optical lattices, electrons in correlated materials, and dense quantum plasmas. Presently, the only quantum-dynamics approach that rigorously describes these processes in two and three dimensions is nonequilibrium Green functions (NEGF). H… ▽ More The dynamics of strongly correlated fermions following an external excitation reveals extremely rich collective quantum effects. Examples are fermionic atoms in optical lattices, electrons in correlated materials, and dense quantum plasmas. Presently, the only quantum-dynamics approach that rigorously describes these processes in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their $T^3$ scaling with the simulation duration $T$. Recently, $T^2$ scaling was achieved with the generalized Kadanoff--Baym ansatz (GKBA) which has substantially extended the scope of NEGF simulations. Here we present a novel approach to GKBA-NEGF simulations that is of order $T$, and demonstrate its remarkable capabilities. △ Less

Submitted 26 November, 2019; v1 submitted 25 September, 2019; originally announced September 2019.

Journal ref: Phys. Rev. Lett. 124, 076601 (2020)

The Static Local Field Correction of the Warm Dense Electron Gas: An ab Initio Path Integral Monte Carlo Study and Machine Learning Representation

Authors: Tobias Dornheim, Jan Vorberger, Simon Groth, Nico Hoffmann, Zhandos Moldabekov, Michael Bonitz

Abstract: The study of matter at extreme densities and temperatures as they occur in astrophysical objects and state-of-the art experiments with high-intensity lasers is of high current interest for many applications. While no overarching theory for this regime exists, accurate data for the density response of correlated electrons to an external perturbation are of paramount importance. In this context, the… ▽ More The study of matter at extreme densities and temperatures as they occur in astrophysical objects and state-of-the art experiments with high-intensity lasers is of high current interest for many applications. While no overarching theory for this regime exists, accurate data for the density response of correlated electrons to an external perturbation are of paramount importance. In this context, the key quantity is given by the local field correction (LFC), which provides a wave-vector resolved description of exchange-correlation effects. In this work, we present extensive new path integral Monte Carlo (PIMC) results for the static LFC of the uniform electron gas, which are subsequently used to train a fully connected deep neural network. This allows us to present a continuous representation of the LFC with respect to wave-vector, density, and temperature covering the entire warm dense matter regime. Both the PIMC data and neural-net results are available online. Moreover, we expect the presented combination of ab initio calculations with machine-learning methods to be a promising strategy for many applications. △ Less

Submitted 19 July, 2019; originally announced July 2019.

Dynamical structure factor of strongly coupled ions in a dense quantum plasma

Authors: Zh. A. Moldabekov, H. Kählert, T. Dornheim, S. Groth, M. Bonitz, T. S. Ramazanov

Abstract: The dynamical structure factor (DSF) of strongly coupled ions in dense plasmas with partially and strongly degenerate electrons is investigated. The main focus is on the impact of electronic correlations (non-ideality) on the ionic DSF. The latter is computed by carrying out molecular dynamics (MD) simulations with a screened ion-ion interaction potential. The electronic screening is taken into ac… ▽ More The dynamical structure factor (DSF) of strongly coupled ions in dense plasmas with partially and strongly degenerate electrons is investigated. The main focus is on the impact of electronic correlations (non-ideality) on the ionic DSF. The latter is computed by carrying out molecular dynamics (MD) simulations with a screened ion-ion interaction potential. The electronic screening is taken into account by invoking the Singwi-Tosi-Land-Sjölander approximation, and compared to the MD simulation data obtained considering the electronic screening in the random phase approximation and using the Yukawa potential. %This allows us to gain insight into the impact of the electronic non-ideality on the ionic DSF. We find that electronic correlations lead to lower values of the ion-acoustic mode frequencies and to an extension of the applicability limit with respect to the wave-number of a hydrodynamic description. Moreover, we show that even in the limit of weak electronic coupling, electronic correlations have a non-negligible impact on the ionic longitudinal sound speed. Additionally, the applicability of the Yukawa potential with an adjustable screening parameter is discussed, which will be of interest, e.g., for the interpretation of experimental results for the ionic DSF of dense plasmas. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. E 99, 053203 (2019)

Ultrafast Dynamics of Strongly Correlated Fermions -- Nonequilibrium Green Functions and Selfenergy Approximations

Authors: N. Schluenzen, S. Hermanns, M. Scharnke, M. Bonitz

Abstract: This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). NEGF, presently, are the only \textit{ab-initio} quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the seco… ▽ More This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). NEGF, presently, are the only \textit{ab-initio} quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the second-order Born approximation (SOA). While they correctly capture the qualitative trends of the relaxation towards equilibrium, the reliability and accuracy of these NEGF simulations has remained open, for a long time. Here we report on recent tests of NEGF simulations for finite lattice systems against exact-diagonalization and density-matrix-renormalization-group benchmark data. The results confirm the high accuracy and predictive capability of NEGF simulations---provided selfenergies are used that go beyond the SOA and adequately include strong correlation and dynamical-screening effects. We present a selfcontained introduction to the theory of NEGF and give an overview on recent numerical applications to compute the ultrafast relaxation dynamics of correlated fermions. In the second part we give a detailed introduction to selfenergies beyond the SOA. Important examples are the third-order approximation, the \GWAx, the \TMA and the fluctuating-exchange approximation. We give a comprehensive summary of the explicit selfenergy expressions for a variety of systems of practical relevance, starting from the most general expressions and the Feynman diagrams, and including also the important cases of diagonal basis sets, the Hubbard model and the differences occuring for bosons and fermions. With these details, and information on the computational effort and scaling with the basis size and propagation duration, an easy use of these approximations in numerical applications is made possible. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Journal ref: J. Phys. Cond. Matt. 32 (10) 2020

Path Integral Monte Carlo Simulation of Degenerate Electrons: Permutation-Cycle Properties

Authors: Tobias Dornheim, Simon Groth, Alexei Filinov, Michael Bonitz

Abstract: Being motivated by the surge of fermionic quantum Monte Carlo simulations at finite temperature, we present a detailed analysis of the permutation-cycle properties of path integral Monte Carlo (PIMC) simulations of degenerate electrons. Particular emphasis is put onto the uniform electron gas in the warm dense matter regime. We carry out PIMC simulations of up to $N=100$ electrons and investigate… ▽ More Being motivated by the surge of fermionic quantum Monte Carlo simulations at finite temperature, we present a detailed analysis of the permutation-cycle properties of path integral Monte Carlo (PIMC) simulations of degenerate electrons. Particular emphasis is put onto the uniform electron gas in the warm dense matter regime. We carry out PIMC simulations of up to $N=100$ electrons and investigate exchange-cycle frequencies, which are found not to follow any simple exponential law even in the case of ideal fermions due to the finite size of the simulation box. Moreover, we introduce a permutation-cycle correlation function, which allows us to analyse the joint probability to simultaneously find cycles of different lengths within a single configuration. Again, we find that finite-size effects predominate the observed behaviour. Finally, we briefly consider an inhomogeneous system, namely electrons in a $2D$ harmonic trap. We expect our results to be of interest for the further development of fermionic PIMC methods, in particular to alleviate the notorious fermion sign problem. △ Less

Submitted 17 February, 2019; originally announced February 2019.

A critique of "quantum dusty plasmas"

Authors: Zh. A. Moldabekov, M. Bonitz, T. S. Ramazanov

Abstract: During the recent decade a steadily increasing number of papers devoted to so-called "quantum dusty plasmas" (QDP) appeared. These systems are combining the properties of dusty (i.e. micrometer to nanometer-sized particles containing) plasmas with quantum effects %(of electrons) appearing at high-density and/or low-temperature. Many exciting properties of QDP were predicted including nontrivial co… ▽ More During the recent decade a steadily increasing number of papers devoted to so-called "quantum dusty plasmas" (QDP) appeared. These systems are combining the properties of dusty (i.e. micrometer to nanometer-sized particles containing) plasmas with quantum effects %(of electrons) appearing at high-density and/or low-temperature. Many exciting properties of QDP were predicted including nontrivial collective oscillations and linear and nonlinear waves. It was predicted that the results have relevance for dense astrophysical plasmas such as white dwarf stars, neutron stars, magnetars and laboratory plasmas in micro- and nano-devices. Here we present a careful examination of the plasma parameters for the expected applications and compare them to the conditions under which a dust particle can exist in a quantum plasma. The conclusion is that quantum dusty plasmas do not exist in reality. △ Less

Submitted 30 January, 2019; originally announced January 2019.

In memoriam Leonid V. Keldysh

Authors: Michael Bonitz, Antti-Pekka Jauho, Michael Sadovski, Sergei Tikhodeev

Abstract: Leonid Keldysh -- one of the most influential theoretical physicists of the 20th century -- passed away in November 2016. Keldysh is best known for the diagrammatic formulation of real-time (nonequilibrium) Green functions theory and for the theory of strong field ionization of atoms. Both theories profoundly changed large areas of theoretical physics and stimulated important experiments. Both the… ▽ More Leonid Keldysh -- one of the most influential theoretical physicists of the 20th century -- passed away in November 2016. Keldysh is best known for the diagrammatic formulation of real-time (nonequilibrium) Green functions theory and for the theory of strong field ionization of atoms. Both theories profoundly changed large areas of theoretical physics and stimulated important experiments. Both these discoveries emerged almost simultaneously -- like Einstein, also Keldysh had his \textit{annus mirabilis} -- the year 1964. But the list of his theoretical developments is much broader and is briefly reviewed here. △ Less

Submitted 4 January, 2019; originally announced January 2019.

Journal ref: phys. stat. sol. (b) 257, issue 7, article 1800600 (2019)

Self-diffusion in two-dimensional quasi-magnetized rotating dusty plasmas

Authors: Peter Hartmann, Jorge C. Reyes, Evdokiya G. Kostadinova, Lorin S. Matthews, Truell W. Hyde, Ranna U. Masheyeva, Karlygash N. Dzhumagulova, Tlekkabul S. Ramazanov, Torben Ott, Hanno Kählert, Michael Bonitz, Ihor Korolov, Zoltán Donkó

Abstract: The self-diffusion phenomenon in a two-dimensional dusty plasma at extremely strong (effective) magnetic fields is studied experimentally and by means of molecular dynamics simulations. In the experiment the high magnetic field is introduced by rotating the particle cloud and observing the particle trajectories in a co-rotating frame, which allows reaching effective magnetic fields up to 3000 Tesl… ▽ More The self-diffusion phenomenon in a two-dimensional dusty plasma at extremely strong (effective) magnetic fields is studied experimentally and by means of molecular dynamics simulations. In the experiment the high magnetic field is introduced by rotating the particle cloud and observing the particle trajectories in a co-rotating frame, which allows reaching effective magnetic fields up to 3000 Tesla. The experimental results confirm the predictions of the simulations: (i) super-diffusive behavior is found at intermediate time-scales and (ii) the dependence of the self-diffusion coefficient on the magnetic field is well reproduced. △ Less

Submitted 26 December, 2018; originally announced December 2018.

Comments: accepted by Physical Review E

\textit{Ab Initio} Path Integral Monte Carlo Results for the Dynamic Structure Factor of Correlated Electrons: From the Electron Liquid to Warm Dense Matter

Authors: T. Dornheim, S. Groth, J. Vorberger, M. Bonitz

Abstract: The accurate description of electrons at extreme density and temperature is of paramount importance for, e.g., the understanding of astrophysical objects and inertial confinement fusion. In this context, the dynamic structure factor $S(\mathbf{q},ω)$ constitutes a key quantity as it is directly measured in X-ray Thomson (XRTS) scattering experiments and governs transport properties like the dynami… ▽ More The accurate description of electrons at extreme density and temperature is of paramount importance for, e.g., the understanding of astrophysical objects and inertial confinement fusion. In this context, the dynamic structure factor $S(\mathbf{q},ω)$ constitutes a key quantity as it is directly measured in X-ray Thomson (XRTS) scattering experiments and governs transport properties like the dynamic conductivity. In this work, we present the first \textit{ab initio} results for $S(\mathbf{q},ω)$ by carrying out extensive path integral Monte Carlo simulations and developing a new method for the required analytic continuation, which is based on the stochastic sampling of the dynamic local field correction $G(\mathbf{q},ω)$. In addition, we find that the so-called static approximation constitutes a promising opportunity to obtain high-quality data for $S(\mathbf{q},ω)$ over substantial parts of the warm dense matter regime. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. Lett. 121, 255001 (2018)

Towards an integrated modeling of the plasma-solid interface

Authors: M Bonitz, A Filinov, J W Abraham, K Balzer, H Kählert, E Pehlke, FX Bronold, M Pamperin, M M Becker, D Loffhagen, H Fehske

Abstract: Solids facing a plasma are a common situation in many astrophysical systems and laboratory setups. Moreover, many plasma technology applications rely on the control of the plasma-surface interaction. However, presently often a fundamental understanding of them is missing, so most technological applications are being developed via trial and error. In the majority of plasma simulations surface proce… ▽ More Solids facing a plasma are a common situation in many astrophysical systems and laboratory setups. Moreover, many plasma technology applications rely on the control of the plasma-surface interaction. However, presently often a fundamental understanding of them is missing, so most technological applications are being developed via trial and error. In the majority of plasma simulations surface processes are either neglected or treated via phenomenological parameters such as sticking coefficients, sputter rates or secondary electron emission coefficients. However, those parameters are known only in some cases and with very limited accuracy. Similarly, while surface physics simulations have often studied the impact of single ions or neutrals, so far, the influence of a plasma medium and correlations between successive impacts have not been taken into account. Such an approach cannot have predictive power. In this paper we discuss in some detail the physical processes a the plasma-solid interface which brings us to the necessity of coupled plasma-solid simulations. We briefly summarize relevant theoretical methods from solid state and surface physics that are suitable to contribute to such an approach and identify four methods. The first are mesoscopic simulations such as kinetic Monte Carlo (KMC) and molecular dynamics (MD) that are able to treat complex processes on large scales but neglect electronic effects. The second are quantum kinetic methods based on the quantum Boltzmann equation that give access to a more accurate treatment of surface processes using simplifying models for the solid. The third approach are ab initio simulations of surface process that are based on density functional theory (DFT) and time-dependent DFT. The fourths are nonequilibrium Green functions that able to treat correlation effects in the material and at the interface. △ Less

Submitted 7 September, 2018; originally announced September 2018.

Journal ref: Frontiers of Chemical Sciences and Engineering 13(2), 201-237 (2019)

Ion Impact Induced Ultrafast Electron Dynamics in Correlated Materials and Finite Graphene Clusters

Authors: M. Bonitz, K. Balzer, N. Schlünzen, M. Rasmussen, J. -P. Joost

Abstract: Strongly correlated systems of fermions have an interesting phase diagram arising from the Hubbard gap. Excitation across the gap leads to the formation of doubly occupied lattice sites (doublons). This state offers interesting electronic and optical properties. Moreover, when the system is driven out of equilibrium interesting collective dynamics may arise that are related to the spatial propagat… ▽ More Strongly correlated systems of fermions have an interesting phase diagram arising from the Hubbard gap. Excitation across the gap leads to the formation of doubly occupied lattice sites (doublons). This state offers interesting electronic and optical properties. Moreover, when the system is driven out of equilibrium interesting collective dynamics may arise that are related to the spatial propagation of doublons. Here, a novel mechanism that was recently proposed by us [Balzer \textit{et al.}, submitted for publication] is verified by exact diagonalization and nonequilibrium Green functions (NEGF) simulations---fermionic doublon creation by the impact of energetic ions. We report the formation of a nonequilibrium steady state with homogeneous doublon distribution. A physically intuitive picture is given in terms of an analytical model for a two-site system where the doublon formation is explained in terms of a two-fold passage of an avoided crossing (Landau-Zener picture). The effect should be particularly important for strongly correlated finite systems, such as graphene nanoribbons, and directly observable with fermionic atoms in optical lattices. We demonstrate that doublon formation and propagation in correlated lattice systems can be accurately simulated with NEGF. In addition to two-time results we present single-time results within the generalized Kadanoff-Baym ansatz (GKBA) with Hartree-Fock propagators (HF-GKBA), and we present systematic improvements that use correlated propagators (correlated GKBA). △ Less

Submitted 19 September, 2018; v1 submitted 23 August, 2018; originally announced August 2018.

Journal ref: phys. stat. sol. (b) 257, issue 7, article 1800490 (2019)

Structural characteristics of strongly coupled ions in a dense quantum plasma

Authors: Zh. A. Moldabekov, S. Groth, T. Dornheim, H. Kaehlert, M. Bonitz, T. S. Ramazanov

Abstract: The structural properties of strongly coupled ions in dense plasmas with moderately to strongly degenerate electrons are investigated in the framework of the one-component plasma model of ions interacting through a screened pair interaction potential. Special focus is put on the description of the electronic screening in the Singwi-Tosi-Land-Sjoelander (STLS) approximation. Different crosschecks a… ▽ More The structural properties of strongly coupled ions in dense plasmas with moderately to strongly degenerate electrons are investigated in the framework of the one-component plasma model of ions interacting through a screened pair interaction potential. Special focus is put on the description of the electronic screening in the Singwi-Tosi-Land-Sjoelander (STLS) approximation. Different crosschecks and analyses using ion potentials obtained from ground-state quantum Monte Carlo data, the random phase approximation (RPA), and existing analytical models are presented for the computation of the structural properties, such as the pair distribution and the static structure factor, of strongly coupled ions. The results are highly sensitive to the features of the screened pair interaction potential. This effect is particularly visible in the static structure factor. The applicability range of the screened potential computed from STLS is identified in terms of density and temperature of the electrons. It is demonstrated that at r_s > 1, where rs is the ratio of the mean inter-electronic distance to the Bohr radius, electronic correlations beyond RPA have a non-negligible effect on the structural properties. Additionally, the applicability of the hypernetted chain approximation for the calculation of the structural properties using the screened pair interaction potential is analyzed employing the effective coupling parameter approach. △ Less

Submitted 29 June, 2018; v1 submitted 30 April, 2018; originally announced April 2018.

Journal ref: Phys. Rev. E 98, 023207 (2018)

Extending \textit{ab initio} plasma-surface simulations to experimentally relevant scales

Authors: M. Bonitz, A. Filinov, J. W. Abraham, D. Loffhagen

Abstract: The physical processes at the interface of a low-temperature plasma and a solid are extremely complex. They involve a huge number of elementary processes in the plasma, in the solid as well as charge, momentum and energy transfer across the interface. In the majority of plasma simulations these surface processes are either neglected or treated via phenomenological parameters. However, those parame… ▽ More The physical processes at the interface of a low-temperature plasma and a solid are extremely complex. They involve a huge number of elementary processes in the plasma, in the solid as well as charge, momentum and energy transfer across the interface. In the majority of plasma simulations these surface processes are either neglected or treated via phenomenological parameters. However, those parameters are known only in some cases, so such an approach is very inaccurate and does not have predictive capability. Therefore, improvements are highly needed. In this paper we briefly summarize relevant theoretical methods from solid state and surface physics that are able to contribute to an improved simulation of plasma-surface interaction in the near future. Full \textit{ab initio} quantum simulations are feasible only for extremely short times and/or small system sizes. A substantial simplification is achieved when electronic quantum effects are not treated explicitly. Then one arrives at semi-classical molecular dynamics (MD) simulations for the heavy particles that have become the main workhorse in surface science simulations. Using microscopically founded potentials and force fields as an input, these MD simulations approach the quality of \textit{ab initio} simulations, in many cases. However, despite their simplified nature, these simulations require a time step that is of the order or below one femtosecond making it prohibitive to reach experimentally relevant scales of minutes. To bridge this gap in length and time scales without compromising the first principles character of the simulations, many physical and computational strategies have been put forward in surface science. This paper presents a brief overview on different methods and their underlying physical ideas, and we compare their strengths and weaknesses. △ Less

Submitted 28 May, 2018; v1 submitted 23 February, 2018; originally announced February 2018.

Journal ref: Plasma Sources Science Technol. 27 No.6, 064005 (2018)

Microscopic modeling of gas-surface scattering. I. A combined molecular dynamics-rate equation approach

Authors: A Filinov, M Bonitz, D Loffhagen

Abstract: A combination of first principle molecular dynamics (MD) simulations with a rate equation model (MD-RE approach) is presented to study the trapping and the scattering of rare gas atoms from metal surfaces. The temporal evolution of the atom fractions that are either adsorbed or scattered into the continuum is investigated in detail. We demonstrate that for this description one has to consider trap… ▽ More A combination of first principle molecular dynamics (MD) simulations with a rate equation model (MD-RE approach) is presented to study the trapping and the scattering of rare gas atoms from metal surfaces. The temporal evolution of the atom fractions that are either adsorbed or scattered into the continuum is investigated in detail. We demonstrate that for this description one has to consider trapped, quasi-trapped and scattering states, and present an energetic definition of these states. The rate equations contain the transition probabilities between the states. We demonstrate how these rate equations can be derived from kinetic theory. Moreover, we present a rigorous way to determine the transition probabilities from a microscopic analysis of the particle trajectories generated by MD simulations. Once the system reaches quasi-equilibrium, the rates converge to stationary values, and the subsequent thermal adsorption/desorption dynamics is completely described by the rate equations without the need to perform further time-consuming MD simulations. As a proof of concept of our approach, MD simulations for argon atoms interacting with a platinum (111) surface are presented. A detailed deterministic trajectory analysis is performed, and the transition rates are constructed. The dependence of the rates on the incidence conditions and the lattice temperature is analyzed. Based on this example, we analyze the time scale of the gas-surface system to approach the quasi-stationary state. The MD-rate equation model has great relevance for the plasma-surface modeling as it makes an extension of accurate simulations to long, experimentally relevant time scales possible. Its application to the computation of atomic sticking probabilities is given in the second part (paper II). △ Less

Submitted 28 May, 2018; v1 submitted 6 February, 2018; originally announced February 2018.

Journal ref: Plasma Sources Science Technol. 27 No.6, 064003 (2018)

Microscopic modeling of gas-surface scattering. II. Application to argon atom adsorption on a platinum (111) surface

Authors: A Filinov, M Bonitz, D Loffhagen

Abstract: A new combination of first principle molecular dynamics (MD) simulations with a rate equation model presented in the preceding paper (paper I) is applied to analyze in detail the scattering of argon atoms from a platinum (111) surface. The combined model is based on a classification of all atom trajectories according to their energies into trapped, quasi-trapped and scattering states. The number o… ▽ More A new combination of first principle molecular dynamics (MD) simulations with a rate equation model presented in the preceding paper (paper I) is applied to analyze in detail the scattering of argon atoms from a platinum (111) surface. The combined model is based on a classification of all atom trajectories according to their energies into trapped, quasi-trapped and scattering states. The number of particles in each of the three classes obeys coupled rate equations. The coefficients in the rate equations are the transition probabilities between these states which are obtained from MD simulations. While these rates are generally time-dependent, after a characteristic time scale $t_E$ of several tens of picoseconds they become stationary allowing for a rather simple analysis. Here, we investigate this time scale by analyzing in detail the temporal evolution of the energy distribution functions of the adsorbate atoms. We separately study the energy loss distribution function of the atoms and the distribution function of in-plane and perpendicular energy components. Further, we compute the sticking probability of argon atoms as a function of incident energy, angle and lattice temperature. Our model is important for plasma-surface modeling as it allows to extend accurate simulations to longer time scales. △ Less

Submitted 13 February, 2018; v1 submitted 6 February, 2018; originally announced February 2018.

Journal ref: Plasma Sources Science Technol. 27 No.6, 064002 (2018)