BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael Bonitz

Michael Bonitz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phcond-mat.str-elcond-mat.stat-mechcond-mat.quant-gasphysics.atom-phphysics.comp-phquant-phastro-ph.SRcond-mat.mes-hallphysics.chem-phphysics.hist-ph

Catalog footprint

What is connected

39works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2023arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum fluctuations approach to the nonequilibrium $GW$ approximation

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].

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Doublon production in correlated materials by multiple ion impacts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nonlinear Electronic Density Response in Warm Dense Matter

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Achieving the Ultimate Scaling Limit for Nonequilibrium Green Functions Simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Advanced fluid modelling and PIC/MCC simulations of low-pressure ccrf discharges

Comparative studies of capacitively coupled radio-frequency discharges in helium and argon at pressures between 10 and 80 Pa are presented applying two different fluid modelling approaches as well as two independently developed particle-in-cell/Monte Carlo collision (PIC/MCC) codes. The focus is on the analysis of the range of applicability of a recently proposed fluid model including an improved drift-diffusion approximation for the electron component as well as its comparison with fluid modelling results using the classical drift-diffusion approximation and benchmark results obtained by PIC/MCC simulations. Main features of this time- and space-dependent fluid model are given. It is found that the novel approach shows generally quite good agreement with the macroscopic properties derived by the kinetic simulations and is largely able to characterize qualitatively and quantitatively the discharge behaviour even at conditions when the classical fluid modelling approach fails. Furthermore, the excellent agreement between the two PIC/MCC simulation codes using the velocity Verlet method for the integration of the equations of motion verifies their accuracy and applicability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

{\em Ab initio} Quantum Monte Carlo simulation of the warm dense electron gas in the thermodynamic limit

We perform \emph{ab initio} quantum Monte Carlo (QMC) simulations of the warm dense uniform electron gas in the thermodynamic limit. By combining QMC data with linear response theory we are able to remove finite-size errors from the potential energy over the entire warm dense regime, overcoming the deficiencies of the existing finite-size corrections by Brown \emph{et al.}~[PRL \textbf{110}, 146405 (2013)]. Extensive new QMC results for up to $N=1000$ electrons enable us to compute the potential energy $V$ and the exchange-correlation free energy $F_{xc}$ of the macroscopic electron gas with an unprecedented accuracy of $|ΔV|/|V|, |ΔF_{xc}|/|F|_{xc} \sim 10^{-3}$. A comparison of our new data to the recent parametrization of $F_{xc}$ by Karasiev {\em et al.} [PRL {\bf 112}, 076403 (2014)] reveals significant deviations to the latter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Correlation effects in strong-field ionization of heteronuclear diatomic molecules

We develop a time-dependent theory to investigate electron dynamics and photoionization processes of diatomic molecules interacting with strong laser fields including electron-electron correlation effects. We combine the recently formulated time-dependent generalized-active-space configuration interaction theory [D. Hochstuhl and M. Bonitz, Phys. Rev. A 86, 053424 (2012); S. Bauch, et al., Phys. Rev. A 90, 062508 (2014)] with a prolate spheroidal basis set including localized orbitals and continuum states to describe the bound electrons and the outgoing photoelectron. As an example, we study the strong-field ionization of the two-center four-electron lithium hydride molecule in different intensity regimes. By using single-cycle pulses, two orientations of the asymmetric heteronuclear molecule are investigated: Li-H, with the electrical field pointing from H to Li, and the opposite case of H-Li. The preferred orientation for ionization is determined and we find a transition from H-Li, for low intensity, to Li-H, for high intensity. The influence of electron correlations is studied at different levels of approximation, and we find a significant change in the preferred orientation. For certain intensity regimes, even an interchange of the preferred configuration is observed, relative to the uncorrelated simulations. Further insight is provided by detailed comparisons of photoelectron angular distributions with and without correlation effects taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stopping dynamics of ions passing through correlated honeycomb clusters

A combined nonequilibrium Green functions-Ehrenfest dynamics approach is developed that allows for a time-dependent study of the energy loss of a charged particle penetrating a strongly correlated system at zero and finite temperature. Numerical results are presented for finite inhomogeneous two-dimensional Fermi-Hubbard models, where the many-electron dynamics in the target are treated fully quantum mechanically and the motion of the projectile is treated classically. The simulations are based on the solution of the two-time Dyson (Keldysh-Kadanoff-Baym) equations using the second-order Born, third-order and T-matrix approximations of the self-energy. As application, we consider protons and helium nuclei with a kinetic energy between 1 and 500 keV/u passing through planar fragments of the two-dimensional honeycomb lattice and, in particular, examine the influence of electron-electron correlations on the energy exchange between projectile and electron system. We investigate the time dependence of the projectile's kinetic energy (stopping power), the electron density, the double occupancy and the photoemission spectrum. Finally, we show that, for a suitable choice of the Hubbard model parameters, the results for the stopping power are in fair agreement with ab-initio simulations for particle irradiation of single-layer graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Permutation blocking path integral Monte Carlo approach to the uniform electron gas at finite temperature

The uniform electron gas (UEG) at finite temperature is of high current interest due to its key relevance for many applications including dense plasmas and laser excited solids. In particular, density functional theory heavily relies on accurate thermodynamic data for the UEG. Until recently, the only existing first-principle results had been obtained for $N=33$ electrons with restricted path integral Monte Carlo (RPIMC), for low to moderate density, $r_s = \overline{r}/a_B \gtrsim 1$. This data has been complemented by Configuration path integral Monte Carlo (CPIMC) simulations for $r_s \leq 1$ that substantially deviate from RPIMC towards smaller $r_s$ and low temperature. In this work, we present results from an independent third method---the recently developed permutation blocking path integral Monte Carlo (PB-PIMC) approach [T. Dornheim \textit{et al.}, NJP \textbf{17}, 073017 (2015)] which we extend to the UEG. Interestingly, PB-PIMC allows us to perform simulations over the entire density range down to half the Fermi temperature ($θ=k_BT/E_F=0.5$) and, therefore, to compare our results to both aforementioned methods. While we find excellent agreement with CPIMC, where results are available, we observe deviations from RPIMC that are beyond the statistical errors and increase with density.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Permutation blocking path integral Monte Carlo: A highly efficient approach to the simulation of strongly degenerate non-ideal fermions

Correlated fermions are of high interest in condensed matter (Fermi liquids, Wigner molecules), cold atomic gases and dense plasmas. Here we propose a novel approach to path integral Monte Carlo (PIMC) simulations of strongly degenerate non-ideal fermions at finite temperature by combining a fourth-order factorization of the density matrix with antisymmetric propagators, i.e., determinants, between all imaginary time slices. To efficiently run through the modified configuration space, we introduce a modification of the widely used continuous space worm algorithm, which allows for an efficient sampling at arbitrary system parameters. We demonstrate how the application of determinants achieves an effective blocking of permutations with opposite signs, leading to a significant relieve of the fermion sign problem. To benchmark the capability of our method regarding the simulation of degenerate fermions, we consider multiple electrons in a quantum dot and compare our results with other ab initio techniques, where they are available. The present permutation blocking path integral Monte Carlo approach allows us to obtain accurate results even for $N=20$ electrons at low temperature and arbitrary coupling, where no other ab initio results have been reported, so far.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statically screened ion potential and Bohm potential in a quantum plasma

The effective potential $Φ$ of a classical ion in a weakly correlated quantum plasma in thermodynamic equilibrium at finite temperature is well described by the RPA screened Coulomb potential. Additionally, collision effects can be included via a relaxation time ansatz (Mermin dielectric function). These potentials are used to study the quality of various statically screened potentials that were recently proposed by Shukla and Eliasson (SE) [Phys. Rev. Lett. {\bf 108}, 165007 (2012)], Akbari--Moghanjoughi (AM) [Phys. Plasmas {\bf 22}, 022103 (2015)] and Stanton and Murillo (SM) [Phys. Rev. E {\bf 91}, 033104 (2015)] starting from quantum hydrodynamic theory (QHD). Our analysis reveals that the SE potential is qualitatively different from the full potential, whereas the SM potential (at any temperature) and the AM potential (at zero temperature) are significantly more accurate. This confirms the correctness of the recently derived [Michta {\em et al.}, Contrib. Plasma Phys. {\bf 55}, (2015)] pre-factor $1/9$ in front of the Bohm term of QHD for fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Effect of Correlations on the Heat Transport in a Magnetized Plasma

In a classical ideal plasma, a magnetic field is known to reduce the heat conductivity perpendicular to the field whereas it does not alter the one along the field. Here we show that, in strongly correlated plasmas that are observed at high pressure or/and low temperature, a magnetic field reduces the perpendicular heat transport much less and even {\it enhances} the parallel transport. These surprising observations are explained by the competition of kinetic, potential and collisional contributions to the heat conductivity. Our results are based on first principle molecular dynamics simulations of a one-component plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Controlling strongly correlated dust clusters with lasers

The most attractive feature of dusty plasmas is the possibility to create strong correlations at room temperatures. At the same time, these plasmas allow for a precise diagnostics with single-particle resolution. From such measurements, the structural properties of finite two-dimensional (2D) clusters and three-dimensional (3D) spherical crystals in nearly harmonic traps-Yukawa balls-have been explored in great detail. Their structural properties-the shell compositions and the order within the shells-have been investigated and good agreement to theoretical predictions was found. Open questions on the agenda are the excitation behavior, the structural changes, and phase transitions that occur at elevated temperature. In order to increase the dust temperature in the experiment various techniques have been used. Among them, laser heating appears to have unique capabilities because it affects only the dust particles, leaving the lighter plasma components unchanged. Here we report on recent experimental results where laser heating methods were further improved and applied to finite 2D and 3D clusters. Comparing to simulations, we demonstrate that this indeed allows to increase the temperature in a controlled manner. For the analysis of thermodynamics and phase transitions in these finite systems, we present theoretical and experimental results on the basis of the instantaneous normal modes, pair distribution function and the recently introduced center-two-particle distribution function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hubbard nanoclusters far from equilibrium

The Hubbard model is a prototype for strongly correlated many-particle systems, including electrons in condensed matter and molecules, as well as for fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the nonequilibrium dynamics following a strong non-perturbative excitation only recently came into the focus of experiments and theory. It is of particular interest how the dynamics depend on the coupling strength and on the particle number and whether there exist universal features in the time evolution. Here, we present results for the dynamics of finite Hubbard clusters based on a selfconsistent nonequilibrium Green functions (NEGF) approach invoking the generalized Kadanoff--Baym ansatz (GKBA). We discuss the conserving properties of the GKBA with Hartree--Fock propagators in detail and present a generalized form of the energy conservation criterion of Baym and Kadanoff for NEGF. Furthermore, we demonstrate that the HF-GKBA cures some artifacts of prior two-time NEGF simulations. Besides, this approach substantially speeds up the numerical calculations and thus presents the capability to study comparatively large systems and to extend the analysis to long times allowing for an accurate computation of the excitation spectrum via time propagation. Our data obtained within the second Born approximation compares favorably with exact diagonalization results (available for up to 13 particles) and are expected to have predictive capability for substantially larger systems in the weak coupling limit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Screened Coulomb potential in a flowing magnetized plasma

The electrostatic potential of a moving dust grain in a complex plasma with magnetized ions is computed using linear response theory, thereby extending our previous work for unmagnetized plasmas [P. Ludwig et al., New J. Phys. 14, 053016 (2012)]. In addition to the magnetic field, our approach accounts for a finite ion temperature as well as ion-neutral collisions. Our recently introduced code \texttt{Kielstream} is used for an efficient calculation of the dust potential. Increasing the magnetization of the ions, we find that the shape of the potential crucially depends on the Mach number $M$. In the regime of subsonic ion flow ($M<1$), a strong magnetization gives rise to a potential distribution that is qualitatively different from the unmagnetized limit, while for $M>1$ the magnetic field effectively suppresses the plasma wakefield.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The ion potential in warm dense matter: wake effects due to streaming degenerate electrons

The effective dynamically screened potential of a classical ion in a stationary flowing quantum plasma at finite temperature is investigated. This is a key quantity for thermodynamics and transport of dense plasmas in the warm dense matter regime. This potential has been studied before within hydrodynamic approaches or based on the zero temperature Lindhard dielectric function. Here we extend the kinetic analysis by including the effects of finite temperature and of collisions based on the Mermin dielectric function. The resulting ion potential exhibits an oscillatory structure with attractive minima (wakes) and, thus, strongly deviates from the static Yukawa potential of equilibrium plasmas. This potential is analyzed in detail for high-density plasmas with values of the Brueckner parameter in the range $0.1 \le r_s \le 1$, for a broad range of plasma temperature and electron streaming velocity. It is shown that wake effects become weaker with increasing temperature of the electrons. Finally, we obtain the minimal electron streaming velocity for which attraction between ions occurs. This velocity turns out to be less than the electron Fermi velocity. Our results allow, for the first time, for reliable predictions of the strength of wake effects in nonequilibrium quantum plasmas with fast streaming electrons showing that these effects are crucial for transport under warm dense matter conditions, in particular for laser-matter interaction, electron-ion temperature equilibration and for stopping power.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Heat Transport in Confined Strongly Coupled 2D Dust Clusters

Dusty plasmas are a model system for studying strong correlation. The dust grains' size of a few micro-meters and their characteristic oscillation frequency of a few hertz allows for an investigation of many particle effects on an atomic level. In this article, we model the heat transport through an axially confined 2D dust cluster from the center to the outside. The system behaves particularly interesting since heat is not only conducted within the dust component but also transfered to the neutral gas. Fitting the analytical solution to the obtained radial temperature profiles allows to determine the heat conductivity $\kheat$. The heat conductivity is found to be constant over a wide range of coupling strengths even including the phase transition from solid to liquid here, as it was also found in extended systems by V. Nosenko et al. in 2008 \cite{PhysRevLett.100.025003}

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Magnetoplasmons in rotating dusty plasmas

A rotating dusty plasma apparatus was constructed to provide the possibility of experimental emulation of extremely high magnetic fields by means of the Coriolis force, observable in a co-rotating measurement frame. We present collective excitation spectra for different rotation rates with a magnetic induction equivalent of up to 3200 Tesla. We identify the onset of magnetoplasmon-equivalent mode dispersion in the rotating macroscopic two-dimensional single-layer dusty plasma. The experimental results are supported by molecular dynamics simulations of 2D magnetized Yukawa systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Breathing Mode of Trapped Systems in One and Two Dimensions

We investigate the quantum breathing mode (monopole oscillation) of trapped fermionic particles with Coulomb and dipole interaction in one and two dimensions. This collective oscillation has been shown to reveal detailed information on the many-particle state of interacting trapped systems and is thus a sensitive diagnostics for a variety of finite systems, including cold atomic and molecular gases in traps and optical lattics, electrons in metal clusters and in quantum confined semiconductor structures or nanoplasmas. An improved sum rule formalism allows us to accurately determine the breathing frequencies from the ground state of the system, avoiding complicated time-dependent simulations. In combination with the Hartree-Fock and the Thomas-Fermi approximations this enables us to extend the calculations to large particle numbers $N$ on the order of several million. Tracing the breathing frequency to large $N$ as a function of the coupling parameter of the system reveals a surprising difference of the asymptotic behavior of one-dimensional and two-dimensional harmonically trapped Coulomb systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Hydrodynamics

Quantum plasma physics is a rapidly evolving research field with a very inter-disciplinary scope of potential applications, ranging from nano-scale science in condensed matter to the vast scales of astrophysical objects. The theoretical description of quantum plasmas relies on various approaches, microscopic or macroscopic, some of which have obvious relation to classical plasma models. The appropriate model should, in principle, incorporate the quantum mechanical effects such as diffraction, spin statistics and correlations, operative on the relevant scales. However, first-principle approaches such as quantum Monte Carlo and density functional theory or quantum-statistical methods such as quantum kinetic theory or non-equilibrium Green's functions require substantial theoretical and computational efforts. Therefore, for selected problems, alternative simpler methods have been put forward. In particular, the collective behavior of many-body systems is usually described within a self-consistent scheme of particles and fields on the mean-field level. In classical plasmas, further simplifications are achieved by a transition to hydrodynamic equations. % Similar fluid-type descriptions for quantum plasmas have been proposed and widely used in the recent decade. This chapter is devoted to an overview of the main concepts of quantum hydrodynamics (QHD), thereby critically analyzing its validity range and its main limitations. Furthermore, the results of the linearized QHD in unmagnetized and magnetized plasmas and a few nonlinear solutions are examined with illustrations. The basic concepts and formulation of particle-particle interactions are also reviewed at the end, indicating their possible consequences in quantum many-body problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collective and single-particle excitations in 2D dipolar Bose gases

The Berezinskii-Kosterlitz-Thouless transition in 2D dipolar systems has been studied recently by path integral Monte Carlo (PIMC) simulations [A. Filinov et al., PRL 105, 070401 (2010)]. Here, we complement this analysis and study temperature-coupling strength dependence of the density (particle-hole) and single-particle (SP) excitation spectra both in superfluid and normal phases. The dynamic structure factor, S(q,omega), of the longitudinal excitations is rigorously reconstructed with full information on damping. The SP spectral function, A(q,omega), is worked out from the one-particle Matsubara Green's function. A stochastic optimization method is applied for reconstruction from imaginary times. In the superfluid regime sharp energy resonances are observed both in the density and SP excitations. The involved hybridization of both spectra is discussed. In contrast, in the normal phase, when there is no coupling, the density modes, beyond acoustic phonons, are significantly damped. Our results generalize previous zero temperature analyses based on variational many-body wavefunctions [F. Mazzanti et al., PRL 102, 110405 (2009), D. Hufnagl et al., PRL 107, 065303 (2011)], where the underlying physics of the excitation spectrum and the role of the condensate has not been addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ion-Streaming Induced Order Transition in 3D Dust Clusters

Dust Dynamics Simulations utilizing a dynamical screening approach are performed to study the effect of ion-streaming on the self-organized structures in a three-dimensional spherically confined complex (dusty) plasma. Varying the Mach number M - the ratio of ion drift velocity to the sound velocity, the simulations reproduce the experimentally observed cluster configurations in the two limiting cases: at M=0 strongly correlated crystalline structures consisting of nested spherical shells (Yukawa balls) and, for M\geq1, flow-aligned dust chains, respectively. In addition, our simulations reveal a discontinuous transition between these two limits. It is found that already a moderate ion drift velocity (M\approx0.1) destabilizes the highly ordered Yukawa balls and initiates an abrupt melting transition. The critical value of M is found to be independent of the cluster size.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Wake Structure in Streaming Complex Plasmas

The theoretical description of complex (dusty) plasmas requires multiscale concepts that adequately incorporate the correlated interplay of streaming electrons and ions, neutrals, and dust grains. Knowing the effective dust-dust interaction, the multiscale problem can be effectively reduced to a one-component plasma model of the dust subsystem. The goal of the present publication is a systematic evaluation of the electrostatic potential distribution around a dust grain in the presence of a streaming plasma environment by means of two complementary approaches: (i) a high precision computation of the dynamically screened Coulomb potential from the dynamic dielectric function, and (ii) full 3D particle-in-cell simulations, which self-consistently include dynamical grain charging and non-linear effects. The applicability of these two approaches is addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum Breathing Mode of Interacting Particles in a One-dimensional Harmonic Trap

Extending our previous work, we explore the breathing mode---the [uniform] radial expansion and contraction of a spatially confined system. We study the breathing mode across the transition from the ideal quantum to the classical regime and confirm that it is not independent of the pair interaction strength (coupling parameter). We present the results of time-dependent Hartree-Fock simulations for 2 to 20 fermions with Coulomb interaction and show how the quantum breathing mode depends on the particle number. We validate the accuracy of our results, comparing them to exact Configuration Interaction results for up to 8 particles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Time-dependent restricted active space Configuration Interaction for the photoionization of many-electron atoms

We introduce the time-dependent restricted active space Configuration Interaction method to solve the time-dependent Schrödinger equation for many-electron atoms, and particularly apply it to the treatment of photoionization processes in atoms. The method is presented in a very general formulation and incorporates a wide range of commonly used approximation schemes, like the single-active electron approximation, time-dependent Configuration Interaction with single-excitations, or the time-dependent R-matrix method. We proof the applicability of the method by calculating the photoionization cross sections of Helium and Beryllium, as well as the X-ray--IR pump-probe ionization in Beryllium

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charge Correlations in a Harmonic Trap

A system of N classical Coulomb charges trapped in a harmonic potential displays shell structure and orientational ordering. The local density profile is well understood from theory, simulation, and experiment. Here, pair correlations are considered for this highly inhomogeneous system for both the fluid and ordered states. In the former, it is noted that there is a close relationship to pair correlations in the uniform OCP. For the ordered state, it is shown that the disordered "tiling" is closely related to the ground state Thomson sites for a single sphere.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theoretical description of spherically confined strongly correlated Yukawa plasmas

A theoretical description of the radial density profile for charged particles with Yukawa interaction in a harmonic trap is described. At strong Coulomb coupling shell structure is observed in both computer simulations and experiments. Correlations responsible for such shell structure are described here using a recently developed model based in density functional theory. A wide range of particle number, Coulomb coupling, and screening lengths is considered within the fluid phase. A hypernetted chain approximation shows the formation of shell structure, but fails to give quantitative agreement with Monte Carlo simulation results at strong coupling. Significantly better agreement is obtained within the hypernetted chain structure using a renormalized coupling constant, representing bridge function corrections.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Finite elements and the discrete variable representation in nonequilibrium Green's function calculations. Atomic and molecular models

In this contribution, we discuss the finite-element discrete variable representation (FE-DVR) of the nonequilibrium Green's function and its implications on the description of strongly inhomogeneous quantum systems. In detail, we show that the complementary features of FEs and the DVR allows for a notably more efficient solution of the two-time Schwinger/Keldysh/Kadanoff-Baym equations compared to a general basis approach. Particularly, the use of the FE-DVR leads to an essential speedup in computing the self-energies. As atomic and molecular examples we consider the He atom and the linear version of H$_3^+$ in one spatial dimension. For these closed-shell models we, in Hartree-Fock and second Born approximation, compute the ground-state properties and compare with the exact findings obtained from the solution of the few-particle time-dependent Schrödinger equation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quantum Potential for Diffraction and Exchange Effects

Semi-classical methods of statistical mechanics can incorporate essential quantum effects by using effective quantum potentials. An ideal Fermi gas interacting with an impurity is represented by a classical fluid with effective electron-electron and electron-impurity quantum potentials. The electron-impurity quantum potential is evaluated at weak coupling, leading to a generalization of the Kelbg potential to include both diffraction and degeneracy effects. The electron-electron quantum potential for exchange effects only is the same as that discussed earlier by others.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Linear Response for Confined Particles

The dynamics of fluctuations is considered for electrons near a positive ion or for charges in a confining trap. The stationary nonuniform equilibrium densities are discussed and contrasted. The linear response function for small perturbations of this nonuniform state is calculated from a linear Markov kinetic theory whose generator for the dynamics is exact in the short time limit. The kinetic equation is solved in terms of an effective mean field single particle dynamics determined by the local density and dynamical screening by a dielectric function for the non-uniform system. The autocorrelation function for the total force on the charges is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Path integral Monte Carlo simulation of charged particles in traps

This chapter is devoted to the computation of equilibrium (thermodynamic) properties of quantum systems. In particular, we will be interested in the situation where the interaction between particles is so strong that it cannot be treated as a small perturbation. For weakly coupled systems many efficient theoretical and computational techniques do exist. However, for strongly interacting systems such as nonideal gases or plasmas, strongly correlated electrons and so on, perturbation methods fail and alternative approaches are needed. Among them, an extremely successful one is the Monte Carlo (MC) method which we are going to consider in this chapter.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Wigner function quantum molecular dynamics

Classical molecular dynamics (MD) is a well established and powerful tool in various fields of science, e.g. chemistry, plasma physics, cluster physics and condensed matter physics. Objects of investigation are few-body systems and many-body systems as well. The broadness and level of sophistication of this technique is documented in many monographs and reviews, see for example \cite{Allan,Frenkel,mdhere}. Here we discuss the extension of MD to quantum systems (QMD). There have been many attempts in this direction which differ from one another, depending on the type of system under consideration. One direction of QMD has been developed for condensed matter systems and will not discussed here, e.g. \cite{fermid}. In this chapter we are dealing with unbound electrons as they occur in gases, fluids or plasmas. Here, one strategy is to replace classical point particles by wave packets, e.g. \cite{fermid,KTR94,zwicknagel06} which is quite successful. At the same time, this method struggles with problems related to the dispersion of such a packet and difficulties to properly describe strong electron-ion interaction and bound state formation. We, therefore, avoid such restrictions and consider a completely general alternative approach. We start discussion of quantum dynamics from a general consideration of quantum distribution functions.

0 citations0 reviewsNo code linkedOpen access