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

M. Bonitz

M. 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-elphysics.atom-phcond-mat.stat-mechcond-mat.quant-gasnucl-thhep-phcond-mat.mes-hallcond-mat.mtrl-sciphysics.comp-phquant-phastro-ph.HEcond-mat.dis-nnhep-latphysics.class-phphysics.flu-dynphysics.gen-ph

Catalog footprint

What is connected

83works
17topics
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

83 published item(s)

preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ab initio simulation of warm dense matter

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ab Initio Quantum Monte Carlo Simulations of the Uniform Electron Gas without Fixed Nodes II: Unpolarized Case

In a recent publication [S. Groth \textit{et al.}, PRB (2016)], we have shown that the combination of two novel complementary quantum Monte Carlo approaches, namely configuration path integral Monte Carlo (CPIMC) [T. Schoof \textit{et al.}, PRL \textbf{115}, 130402 (2015)] and permutation blocking path integral Monte Carlo (PB-PIMC) [T. Dornheim \textit{et al.}, NJP \textbf{17}, 073017 (2015)], allows for the accurate computation of thermodynamic properties of the spin-polarized uniform electron gas (UEG) over a wide range of temperatures and densities without the fixed-node approximation. In the present work, we extend this concept to the unpolarized case, which requires non-trivial enhancements that we describe in detail. We compare our new simulation results with recent restricted path integral Monte Carlo data [E. Brown \textit{et al}., PRL \textbf{110}, 146405 (2013)] for different energy contributions and pair distribution functions and find, for the exchange correlation energy, overall better agreement than for the spin-polarized case, while the separate kinetic and potential contributions substantially deviate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Molecular dynamics simulation of gold cluster growth during sputter deposition

We present a molecular dynamics simulation scheme that we apply to study the time evolution of the self-organized growth process of metal cluster assemblies formed by sputter-deposited gold atoms on a planar surface. The simulation model incorporates the characteristics of the plasma-assisted deposition process and allows for an investigation over a wide range of deposition parameters. It is used to obtain data for the cluster properties which can directly be compared to recently published experimental data for gold on polystyrene (M. Schwartzkopf \textit{et al}., ACS Appl. Mater. Interfaces \textbf{7}, 13547 (2015)). While good agreement is found between the two, the simulations additionally provide valuable time-dependent real-space data of the surface morphology some of whose details are hidden in the reciprocal-space scattering images that were used for the experimental analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab Initio Quantum Monte Carlo Simulations of the Uniform Electron Gas without Fixed Nodes

The uniform electron gas (UEG) at finite temperature is of key relevance for many applications in the warm dense matter regime, e.g. dense plasmas and laser excited solids. Also, the quality of density functional theory calculations crucially relies on the availability of accurate data for the exchange-correlation energy. Recently, new benchmark results for the N = 33 spin-polarized electrons at high density, r_s = r/a_B <= 4 and low temperature, have been obtained with the configuration path integral Monte Carlo (CPIMC) method [T. Schoof et al., Phys. Rev. Lett. 115, 130402 (2015)]. To achieve these results, the original CPIMC algorithm [T. Schoof et al., Contrib. Plasma Phys. 51, 687 (2011)] had to be further optimized to cope with the fermion sign problem (FSP). It is the purpose of this paper to give detailed information on the manifestation of the FSP in CPIMC simulations of the UEG and to demonstrate how it can be turned into a controllable convergence problem. In addition, we present new thermodynamic results for higher temperatures. Finally, to overcome the limitations of CPIMC towards strong coupling, we invoke an independent method|the recently developed permutation blocking path integral Monte Carlo approach [T. Dornheim et al., accepted for publication in J. Chem Phys., arXiv:1508.03221]. The combination of both approaches is able to yield ab initio data for the UEG over the entire density range, above a temperature of about one half of the Fermi temperature. Comparison with restricted path integral Monte Carlo data [E. W. Brown et al., Phys. Rev. Lett. 110, 146405 (2013)] allows us to quantify the systematic error arising from the free particle nodes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab initio thermodynamic results for the degenerate electron gas at finite temperature

The uniform electron gas (UEG) at finite temperature is of key relevance for many applications in dense plasmas, warm dense matter, laser excited solids and much more. Accurate thermodynamic data for the UEG are an essential ingredient for many-body theories, in particular, density functional theory. Recently, first-principle restricted path integral Monte Carlo results became available which, however, due to the fermion sign problem, had to be restricted to moderate degeneracy, i.e. low to moderate densities with $r_s={\bar r}/a_B \gtrsim 1$. Here we present novel first-principle configuration PIMC results for electrons for $r_s \leq 1$. We also present quantum statistical data within the $e^4$-approximation that are in good agreement with the simulations at small to moderate $r_s$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab initio transport results for strongly correlated fermions

Quantum transport of strongly correlated fermions is of central interest in condensed matter physics. Here, we present first-principle nonequilibrium Green functions results using $T$-matrix selfenergies for finite Hubbard clusters of dimension $1,2,3$. We compute the expansion dynamics following a potential quench and predict its dependence on the interaction strength and particle number. We discover a universal scaling, allowing an extrapolation to infinite-size systems, which shows excellent agreement with recent cold atom diffusion experiments [Schneider et al., Nat. Phys. 8, 213 (2012)].

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Influence of the nature of confinement on the melting of Wigner molecules in quantum dots

We analyze the quantum melting of two-dimensional Wigner molecules (WM) in confined geometries with distinct symmetries and compare it with corresponding thermal melting. Our findings unfold complementary mechanisms that drive the quantum and thermal crossovers in a WM and show that the symmetry of the confinement plays no significant role in determining the quantum crossover scale $n_X$. This is because the zero-point motion screens the boundary effects within short distances. The phase diagram as a function of thermal and quantum fluctuations determined from independent criteria is unique, and shows "melting" from the WM to both the classical and quantum "liquids." An intriguing signature of weakening liquidity with increasing temperature, $T$, is found in the extreme quantum regime. The crossover is associated with production of defects. However, these defects appear to play distinct roles in driving the quantum and thermal "melting." Our study will help comprehending melting in a variety of experimental traps - from quantum dots to complex plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum hydrodynamics for plasmas -- a Thomas-Fermi theory perspective

The idea to describe quantum systems within a hydrodynamic framework (quantum hydrodynamics, QHD) goes back to Madelung and Bohm. While such a description is formally exact for a single particle, more recently the concept has been applied to many-particle systems by Manfredi and Haas [Phys. Rev. B {\bf 64}, 075316 (2001)] and received high popularity in parts of the quantum plasma community. Thereby, often the applicability limits of these equations are ignored, giving rise to unphysical predictions. Here we demonstrate that modified QHD equations for plasmas can be derived from Thomas-Fermi theory including gradient corrections. This puts QHD on firm grounds. At the same time this derivation yields a different prefactor, $γ=(D-2/3D)$, in front of the quantum (Bohm) potential which depends on the system dimensionality $D$. Our approach allows one to identify the limitations of QHD and to outline systematic improvements.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolving structural transitions in spherical dust clusters

Finite systems in confining potentials are known to undergo structural transitions similar to phase transitions. However, these systems are inhomogeneous, and their "melting" point may depend on the position in the trap and vary with the particle number. Focusing on three-dimensional Coulomb systems in a harmonic trap a rich physics is revealed: in addition to radial melting we demonstrate the existence of intrashell disordering and inter-shell angular melting. Our analysis takes advantage of a novel melting criterion that is based on the spatial two and three-particle distribution functions and the associated reduced entropy which can be directly measured in complex plasma experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coupling Strength in Coulomb and Yukawa One-Component Plasmas

In a non-ideal classical Coulomb one-component plasma (OCP) all thermodynamic properties are known to depend only on a single parameter -- the coupling parameter $Γ$. In contrast, if the pair interaction is screened by background charges (Yukawa OCP) the thermodynamic state depends, in addition, on the range of the interaction via the screening parameter $κ$. How to determine in this case an effective coupling parameter has been a matter of intensive debate. Here we propose a consistent approach for defining and measuring the coupling strength in Coulomb and Yukawa OCPs based on a fundamental structural quantity, the radial pair distribution function (RPDF). The RPDF is often accessible in experiments by direct observation or indirectly through the static structure factor. Alternatively, it is directly computed in theoretical models or simulations. Our approach is based on the observation that the build-up of correlation from a weakly coupled system proceeds in two steps: First, a monotonically increasing volume around each particle becomes devoid of other particles (correlation hole), and second (upon further increase of the coupling), a shell structure emerges around each particle giving rise to growing peaks of the RPDF. Using molecular dynamics simulation, we present a systematic study for the dependence of these features of the RPDF on $Γ$ and $κ$ and derive a simple expression for the effective coupling parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fermionic path integral Monte Carlo results for the uniform electron gas at finite temperature

The uniform electron gas (UEG) at finite temperature has recently attracted substantial interest due to the epxerimental progress in the field of warm dense matter. To explain the experimental data accurate theoretical models for high density plasmas are needed which crucially depend on the quality of the thermodynamic properties of the quantum degenerate correlated electrons. Recent fixed node path integral Monte Carlo (RPIMC) data are the most accurate for the UEG at finite temperature, but they become questionable at high degeneracy when the Brueckner parameter $r_s$ becomes smaller than $1$. Here we present new improved direct fermionic PIMC simulations that are exptected to be more accurate than RPIMC at high densities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Introduction to Configuration Path Integral Monte Carlo

In low-temperature high-density plasmas quantum effects of the electrons are becoming increasingly important. This requires the development of new theoretical and computational tools. Quantum Monte Carlo methods are among the most successful approaches to first-principle simulations of many-body quantum systems. In this chapter we present a recently developed method---the configuration path integral Monte Carlo (CPIMC) method for moderately coupled, highly degenerate fermions at finite temperatures. It is based on the second quantization representation of the $N$-particle density operator in a basis of (anti-)symmetrized $N$-particle states (configurations of occupation numbers) and allows to tread arbitrary pair interactions in a continuous space. We give a detailed description of the method and discuss the application to electrons or, more generally, Coulomb-interacting fermions. As a test case we consider a few quantum particles in a one-dimensional harmonic trap. Depending on the coupling parameter (ratio of the interaction energy to kinetic energy), the method strongly reduces the sign problem as compared to direct path integral Monte Carlo (DPIMC) simulations in the regime of strong degeneracy which is of particular importance for dense matter in laser plasmas or compact stars. In order to provide a self-contained introduction, the chapter includes a short introduction to Metropolis Monte Carlo methods and the second quantization of quantum mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Superfluidity of strongly correlated bosons in two- and three-dimensional traps

We analyze the superfluid phase transition of harmonically confined bosons with long-range interaction in both two and three dimensions in a broad parameter range from weak to strong coupling. We observe that the onset of superfluidity occurs in $3D$ at significantly lower temperatures compared to $2D$. This is demonstrated to be a quantum degeneracy effect. In addition, the spatial distribution of superfluidity across the shells of the clusters is investigated. It is found that superfluidity is substantially reduced in the outer layers due to increased correlation effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thermodynamics of the quark-gluon plasma at finite chemical potential: color path integral Monte Carlo results

Based on the constituent quasiparticle model of the quark-gluon plasma (QGP), color quantum path-integral Monte-Carlo (PIMC) calculations of the thermodynamic properties of the QGP are performed. We extend our previous zero chemical potential simulations to the QGP at finite baryon chemical potential. The results indicate that color PIMC can be applied not only above the QCD critical temperature $T_c$ but also below $T_c$. Besides reproducing the lattice equation of state our approach yields also valuable additional insight into the internal structure of the QGP, via the pair distribution functions of the various quasiparticles. In particular, the pair distribution function of gluons reflects the existence of gluon-gluon bound states at low temperatures and $μ=175$ MeV, i.e. glueballs, while meson-like bound states are not found.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Toward a Nonequilibrium Green functions approach to diffusion in strongly coupled finite quantum systems

Transport properties of strongly correlated quantum systems are of central interest in condensed matter, ultracold atoms and in dense plasmas. There, the proper treatment of strong correlations poses a great challenge to theory. Here, we apply a Nonequilibrium Green Functions approach using a lattice model as a basic system. This allow us to treat a finite spatially inhomogeneous system with an arbitrary nonequilibrium initial state. Placing all particles initially to one side of the system allows for a nonequilibrium study of diffusion. Strong correlation effects are incorporated via T-matrix selfenergies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Towards ab initio thermodynamics of the electron gas at strong degeneracy

Recently a number of theoretical studies of the uniform electron gas (UEG) at finite temperature have appeared that are of relevance for dense plasmas, warm dense matter and laser excited solids and thermodynamic density functional theory simulations. In particular, restricted path integral Monte Carlo (RPIMC) results became available which, however, due to the Fermion sign problem, are confined to moderate quantum degeneracy, i.e. low to moderate densities. We have recently developed an alternative approach---configuration PIMC [T. Schoof {\em et al.}, Contrib. Plasma Phys. {\bf 51}, 687 (2011)] that allows one to study the so far not accessible high degeneracy regime. Here we present the first step towards UEG simulations using CPIMC by studying implementation and performance of the method for the model case of $N=4$ particles. We also provide benchmark data for the total energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ultrafast dynamics of finite Hubbard clusters - a stochastic mean-field approach

Finite lattice models are a prototype for strongly correlated quantum systems and capture essential properties of condensed matter systems. With the dramatic progress in ultracold atoms in optical lattices, finite fermionic Hubbard systems have become directly accessible in experiments, including their ultrafast dynamics far from equilibrium. Here, we present a theoretical approach that is able to treat these dynamics in any dimension and fully includes inhomogeneity effects. The method consists in stochastic sampling of mean-field trajectories and is found to be more accurate and efficient than current nonequilibrium Green functions approaches. This is demonstrated for Hubbard clusters with up to 512 particles in one, two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Color path-integral Monte-Carlo simulations of quark-gluon plasma: Thermodynamic and transport properties

Based on the quasiparticle model of the quark-gluon plasma (QGP), a color quantum path-integral Monte-Carlo (PIMC) method for calculation of thermodynamic properties and -- closely related to the latter -- a Wigner dynamics method for calculation of transport properties of the QGP are formulated. The QGP partition function is presented in the form of a color path integral with a new relativistic measure instead of the Gaussian one traditionally used in the Feynman-Wiener path integral. It is shown that the PIMC method is able to reproduce the lattice QCD equation of state at zero baryon chemical potential at realistic model parameters (i.e. quasiparticle masses and coupling constant) and also yields valuable insight into the internal structure of the QGP. Our results indicate that the QGP reveals quantum liquid-like (rather than gas-like) properties up to the highest considered temperature of 525 MeV. The pair distribution functions clearly reflect the existence of gluon-gluon bound states, i.e. glueballs, at temperatures just above the phase transition, while meson-like $q\bar{q}$ bound states are not found. The calculated self-diffusion coefficient agrees well with some estimates of the heavy-quark diffusion constant available from recent lattice data and also with an analysis of heavy-quark quenching in experiments on ultrarelativistic heavy ion collisions, however, appreciably exceeds other estimates. The lattice and heavy-quark-quenching results on the heavy-quark diffusion are still rather diverse. The obtained results for the shear viscosity are in the range of those deduced from an analysis of the experimental elliptic flow in ultrarelativistic heavy ions collisions, i.e. in terms the viscosity-to-entropy ratio, $1/4π< η/S < 2.5/4π$, in the temperature range from 170 to 440 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comment on "Discussion on `Novel attractive force between ions in quantum plasmas -- failure of simulations based on a density functional approach"

In a recent article [P.K. Shukla, B. Eliasson and M. Akbari-Moghanjoughi, Physica Scripta {\bf 87}, 018202 (2013)] the authors criticized our analysis of the screened proton potential in dense hydrogen that was based on {\em ab initio} density functional theory (DFT) simulations [M. Bonitz, E. Pehlke, and T. Schoof, Phys. Rev. E {\bf 87}, 037102 (2013)]. In particular, they attributed the absence of the Shukla-Eliasson attractive force between protons in the DFT simulations to a failure of DFT. Here we discuss in detail their arguments and show that their conclusions are incorrect.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of Hubbard nano-clusters following strong excitation

The Hubbard model is a prototype for strongly correlated electrons in condensed matter, for molecules and fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the excitation and relaxation dynamics following a perturbation of the system are only poorly explored. Here, we present results for the dynamics of electrons following nonlinear strong excitation that are based on a nonequilibrium Green functions approach. We focus on small systems---"Hubbard nano-clusters"---that contain just a few particles where, in addition to the correlation effects, finite size effects and spatial inhomegeneity can be studied systematically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of two-dimensional complex plasmas in a magnetic field

We consider a two-dimensional complex plasma layer containing charged dust particles in a perpendicular magnetic field. Computer simulations of both one-component and binary systems are used to explore the equilibrium particle dynamics in the fluid state. The mobility is found to scale with the inverse of the magnetic field strength (Bohm diffusion) for strong fields. For bidisperse mixtures, the magnetic field dependence of the long-time mobility depends on the particle species providing an external control of their mobility ratio. For large magnetic fields, even a two-dimensional model porous matrix can be realized composed by the almost immobilized high-charge particles which act as obstacles for the mobile low-charge particles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Formation of magnetic nanocolumns during vapor phase deposition of a metal-polymer nanocomposite: experiments and kinetic Monte Carlo simulations

Metal-polymer nanocomposites have been investigated extensively during the last years due to their interesting functional applications. They are often produced by vapor phase deposition which generally leads to the self-organized formation of spherical metallic nanoparticles in the organic matrix, while nanocolumns are only obtained under very specific conditions. Experiments\cite{Grev+06} have shown that co-evaporation of the metallic and organic components in a simple single-step process can give rise to the formation of ultrahigh-density Fe-Ni-Co nanocolumnar structures embedded in a fluoropolymer matrix. Here we present a kinetic Monte Carlo approach which is based on an new model involving the depression of the melting point on the nanoscale and a critical nanoparticle size required for solidification. In addition we present new experimental results down to a deposition temperature of \cel{-70} and also report the magnetic properties. The simulations provide a detailed understanding of the transition from spherical cluster growth to formation of elongated structures and are in quantitative agreement with the experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On "Novel attractive forces" between ions in quantum plasmas -- failure of linearized quantum hydrodynamics

In a recent letter [P.K. Shukla and B. Eliasson, Phys. Rev. Lett. 108, 165007 (2012)] the discovery of a new attractive force between protons in a hydrogen plasma was reported that would be responsible for the formation of molecules and of a proton lattice. Here we show, based on ab initio density functional theory calculations, that these predictions are wrong and caused by using linearized quantum hydrodynamics beyond the limits of its applicability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reply to comment of Shukla et al. on 'On "Novel attractive forces" between ions in quantum plasmas - failure of linearized quantum hydrodynamics'

In an earlier paper we demonstrated that the "novel attractive force" between protons in dense hydrogen reported by Shukla and Eliasson based on linearized quantum hydrodynamics (LQHD) is wrong. As benchmark results we used state of the art density functional theory (DFT) simulations. In their response to that paper Shukla et al. (arXiv:1112.5556) claimed that the disagreement is due to a failure of DFT. Here we show - by comparing the properties and limitations of DFT and QHD - that their statements have no scientific basis and that DFT is a suitable benchmark for more approximate approaches such as QHD and LQHD.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Color path-integral Monte Carlo simulations of quark-gluon plasma

Thermodynamic properties of a strongly coupled quark-gluon plasma (QGP) of constituent quasiparticles are studied by a color path-integral Monte Carlo simulations (CPIMC). For our simulations we have presented QGP partition function in the form of color path integral with new relativistic measure instead of Gaussian one used in Feynman and Wiener path integrals. For integration over color variable we have also developed procedure of sampling color variables according to the group SU(3) Haar measure. It is shown that this method is able to reproduce the available quantum lattice chromodynamics (QCD) data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electronic double-excitations in quantum wells: solving the two-time Kadanoff-Baym equations

For a quantum many-body system, the direct population of states of double-excitation character is a clear indication that correlations importantly contribute to its nonequilibrium properties. We analyze such correlation-induced transitions by propagating the nonequilibrium Green's functions in real-time within the second Born approximation. As crucial benchmarks, we compute the absorption spectrum of few electrons confined in quantum wells of different width. Our results include the full two-time solution of the Kadanoff-Baym equations as well as of their time-diagonal limit and are compared to Hartree-Fock and exact diagonalization data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Few-particle quantum dynamics---comparing Nonequilibrium Green's functions with the generalized Kadanoff-Baym ansatz to density operator theory

The nonequilibrium description of quantum systems requires, for more than two or three particles, the use of a reduced description to be numerically tractable. Two possible approaches are based on either reduced density matrices or nonequilibrium Green's functions (NEGF). Both concepts are formulated in terms of hierarchies of coupled equations---the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy for the reduced density operators and the Martin-Schwinger-hierarchy (MS) for the Green's functions, respectively. In both cases, similar approximations are introduced to decouple the hierarchy, yet still many questions regarding the correspondence of both approaches remain open.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetized strongly coupled plasmas and how to realize them in a dusty plasma setup

Strongly coupled plasmas in which the interaction energy exceeds the kinetic energy play an important role in many astrophysical and laboratory systems including compact stars, laser plasmas and dusty plasmas. They exhibit many unusual collective properties, such as liquid or crystalline behaviour, peculiar oscillation spectra and transport properties. Recently, strongly coupled plasmas were studied in the presence of a strong magnetic field by computer simulations, and strong modifications of their transport properties and oscillation spectra were observed. While strong magnetization is common in stellar systems it is practically impossible to achieve in complex plasmas due to the large mass of the dust particles. Here we discuss a recently demonstrated approach to achieve very strong "magnetization" by a rotation of the neutral gas, and we present new results for macroscopic two-dimensional systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetizing a complex plasma without a magnetic field

We propose and demonstrate a concept that mimics the magnetization of the heavy dust particles in a complex plasma while leaving the properties of the light species practically unaffected. It makes use of the frictional coupling between a complex plasma and the neutral gas, which allows to transfer angular momentum from a rotating gas column to a well-controlled rotation of the dust cloud. This induces a Coriolis force that acts exactly as the Lorentz force in a magnetic field. Experimental normal mode measurements for a small dust cluster with four particles show excellent agreement with theoretical predictions for a magnetized plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-equilibrium Green's function approach to inhomogeneous quantum many-body systems using the Generalized Kadanoff Baym Ansatz

In non-equilibrium Green's function calculations the use of the Generalized Kadanoff-Baym Ansatz (GKBA) allows for a simple approximate reconstruction of the two-time Green's function from its time-diagonal value. With this a drastic reduction of the computational needs is achieved in time-dependent calculations, making longer time propagation possible and more complex systems accessible. This paper gives credit to the GKBA that was introduced 25 years ago. After a detailed derivation of the GKBA, we recall its application to homogeneous systems and show how to extend it to strongly correlated, inhomogeneous systems. As a proof of concept, we present results for a 2-electron quantum well, where the correct treatment of the correlated electron dynamics is crucial for the correct description of the equilibrium and dynamic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonequilibrium Green's function approach to the pair distribution function of quantum many-body systems out of equilibrium

The pair distribution function (PDF) is a key quantity for the analysis of correlation effects of a quantum system both in equilibrium and far from equilibrium. We derive an expression for the PDF in terms of the single-particle Green's functions---the solutions of the Keldysh/Kadanoff-Baym equations in the two-time plane---for a one- or two-component system. The result includes initial correlations and generalizes previous density matrix expressions from single-time quantum kinetic theory. Explicit expressions for the PDF are obtained in second Born approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase diagram of bilayer electron-hole plasmas

We investigate exciton bound-state formation and crystallization effects in two-dimensional electron-hole bilayers. Performing unbiased path integral Monte Carlo simulations all quantum and Coulomb correlation effects are treated on first principles. We analyze diverse pair distribution functions in dependence on the layer separation, particle density and hole-to-electron mass ratio and derive a schematic phase diagram for the neutral mass-asymmetric bilayer system. Our simulations reveal a great variety of possible phases namely an exciton gas, an exciton crystal, an electron-hole liquid and a hole crystal embedded in an electron gas.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The generalized Kadanoff-Baym ansatz. Computing nonlinear response properties of finite systems

For a minimal Hubbard-type system at different interaction strengths U, we investigate the density-response for an excitation beyond the linear regime using the generalized Kadanoff-Baym ansatz (GKBA) and the second Born (2B) approximation. We find strong correlation features in the response spectra and establish the connection to an involved double excitation process. By comparing approximate and exact Green's function results, we also observe an anomalous U-dependence of the energy of this double excitation in 2B+GKBA. This is in accordance with earlier findings [K. Balzer et al., EPL 98, 67002 (2012)] on double excitations in quantum wells.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Theoretical description of field-assisted post-collision interaction in Auger decay of atoms

In a recent publication [B. Schuette, S. Bauch et al., accepted for publication in Phys. Rev. Lett. 2012] it was demonstrated both, experimentally and theoretically, that Auger electrons are subject to an energetic chirp if a three-body-interaction dubbed post-collision interaction (PCI) is involved. Here, we extend previous theoretical work and give a detailed analysis of field-assisted PCI based on numerical solutions of the time-dependent Schrödinger equation, extensive Monte-Carlo averaged molecular dynamics simulations and analytical theory. The dependence on various streaking and excitation conditions is investigated, and we discuss, how these findings may help to improve XUV pulse characterization as well as understanding of ultrafast atomic processes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Crystallization of an exciton superfluid

Indirect excitons -- pairs of electrons and holes spatially separated in semiconductor bilayers or quantum wells -- are known to undergo Bose-Einstein condensation and to form a quantum fluid. Here we show that this superfluid may crystallize upon compression. However, further compression results in quantum melting back to a superfluid. This unusual behavior is explained by the effective interaction potential between indirect excitons which strongly deviates from a dipole potential at small distances due to many-particle and quantum effects. Based on first principle path integral Monte Carlo simulations, we compute the complete phase diagram of this system and predict the relevant parameters necessary to experimentally observe exciton crystallization in semiconductor quantum wells.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diffusion and growth of metal clusters in nanocomposites: a Kinetic Monte Carlo study

Nobel metals that are deposited on a polymer surface exhibit surface diffusion and diffusion into the bulk. At the same time the metal atoms tend to form clusters because their cohesive energy is about two orders of magnitude higher than the cohesive energy of polymers. To selfconsistently simulate these coupled processes, we present in this paper a Kinetic Monte Carlo approach. Using a simple model with diffusion coefficients taken as input parameters allows us to perform a systematic study of the behavior of a large ensemble of metal atoms on a polymer surface eventually leading to polymer nanocomposites. Special emphasis is placed on the cluster growth, cluster size distribution and the penetration of clusters into the substrate. We also study the influence of surface defects and analyze how the properties of the resulting material can be controlled by variation of the deposition rate.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diffusion in a strongly coupled magnetized plasma

A first-principle study of diffusion in a strongly coupled one-component plasma in a magnetic field B is presented. As in a weakly coupled plasma, the diffusion coefficient perpendicular to the field exhibits a Bohm-like 1/B-behavior in the strong-field limit, but its overall scaling is substantially different. The diffusion coefficient parallel to the field is strongly affected by the field as well and also approaches a 1/B scaling, in striking contrast to earlier predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Higher Harmonics Generation in Strongly Coupled Magnetized Two-Dimensional Yukawa Liquids

The excitation spectra of two-dimensional strongly coupled (liquid-like) Yukawa systems under the influence of a magnetic field perpendicular to the plane are found to sustain additional high-frequency modes at multiples of the magnetoplasmon. These modes are reminiscent of the well-known Bernstein modes but show a number of important differences due to the strong coupling of the particles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Oscillation Spectrum of a Magnetized Strongly Coupled One-Component-Plasma

A first-principle study of the collective oscillation spectrum of a strongly correlated one-component plasma in a strong magnetic field is presented. The spectrum consists of six fundamental modes which are found to be in good agreement with results from the Quasi-Localized Charge approximation. At high frequencies additional modes are observed that include Bernstein-type oscillations and their transverse counterparts which are of importance for the high-frequency optical and transport properties of these plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum simulations of thermodynamic properties of strongly coupled quark-gluon plasma

A strongly coupled quark-gluon plasma (QGP) of heavy constituent quasi-particles is studied by a path-integral Monte-Carlo method. This approach is a quantum generalization of the model developed by Gelman, Shuryak and Zahed. It is shown that this method is able to reproduce the QCD lattice equation of state and also yields valuable insight into the internal structure of the QGP. The results indicate that the QGP reveals liquid-like rather than gas-like properties. At temperatures just above the critical one it was found that bound quark-antiquark states still survive. These states are bound by effective string-like forces and turns out to be colorless. At the temperature as large as twice the critical one no bound states are observed. Quantum effects turned out to be of prime importance in these simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Berezinskii-Kosterlitz-Thouless transition in two-dimensional dipole systems

The superfluid to normal fluid transition of dipolar bosons in two dimensions is studied throughout the whole density range using path integral Monte Carlo simulations and summarized in the phase diagram. While at low densities, we find good agreement with the universal results depending only on the scattering length $a_s$, at moderate and high densities, the transition temperature is strongly affected by interactions and the elementary excitation spectrum. The results are expected to be of relevance to dipolar atomic and molecular systems and indirect excitons in quantum wells.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Collective excitations of a spherically confined Yukawa plasma

The results of a recent fluid theory for the multipole modes of a Yukawa plasma in a spherical confinement [H. Kählert and M. Bonitz, Phys. Rev. E \textbf{82}, 036407 (2010)] are compared with molecular dynamics simulations and the exact $N$-particle eigenmodes in the crystalline phase. Simulations confirm the existence of high order modes found in cold fluid theory. We investigate the influence of screening, coupling and friction on the mode spectra in detail. Good agreement between theory and simulation is found for weak to moderate screening and low order modes. The relations between the breathing mode in the fluid theory, simulation and the crystal eigenmode are investigated in further detail.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effective coupling parameter for 2D Yukawa liquids and non-invasive measurement of plasma parameters

We derive an effective coupling parameter for two-dimensional Yukawa systems based on the height of the first maximum of the pair distribution function. Two variants -- one valid in the high-coupling range, the other for arbitrary couplings of the liquid -- are derived. Comparison to previous approaches to Yukawa coupling parameters shows that the present concept is more general and more accurate. Using, in addition, dynamical information contained in the velocity autocorrelation function, we outline a reference data method that can be employed as a non-invasive measurement scheme of the plasma parameters -- the coupling strength and the screening length. This approach requires only input from a time-series of configuration snapshots and particle velocities with no recourse to additional information about the system. Our results should be directly applicable as a simple, yet reliable diagnostic method for a variety of experiments, including dusty plasmas, colloidal suspensions and ions in traps, and can be employed to facilitate comparisons between experiments, theory and simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electronic correlations in double ionization of atoms in pump-probe experiments

The ionization dynamics of a two-electron atom in an attosecond XUV-infrared pump-probe experiment is simulated by solving the time-dependent two-electron Schrödinger equation. A dramatic change of the double ionization (DI) yield with variation of the pump-probe delay is reported and the governing role of electron-electron correlations is shown. The results allow for a direct control of the DI yield and of the relative strength of double and single ionization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fluid Modes of a Spherically Confined Yukawa Plasma

The normal modes of a three-dimensional Yukawa plasma in an isotropic, harmonic confinement are investigated by solving the linearized cold fluid equations. The eigenmodes are found analytically and expressed in terms of hypergeometric functions. It is found that the mode frequencies solely depend on the dimensionless plasma parameter $ξ=κR$, where $R$ is the plasma radius and $κ$ the inverse screening length. The eigenfrequencies increase monotonically with $ξ$ and saturate in the limit $ξ\to\infty$. Compared with the results in the Coulomb limit~[D. H. E. Dubin, Phys. Rev. Lett. \textbf{66}, 2076 (1991)], we find a new class of modes characterized by the number $n$ which determines the number of radial nodes in the perturbed potential. These modes originate from the degenerate bulk modes of the Coulomb system. Analytical formulas for the eigenfrequencies are derived for limiting cases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Higher harmonics of the magnetoplasmon in strongly coupled Coulomb and Yukawa systems

The generation of higher harmonics of the magnetoplasmon frequency which has recently been reported in strongly coupled two-dimensional Yukawa systems is investigated in detail and, in addition, extended to two-dimensional Coulomb systems. We observe higher harmonics over a much larger frequency range than before and compare the theoretical prediction with the simulations. The influence of the coupling, structure, and thermal energy on the excitation of these modes is examined in detail. We also report on the effect of friction on the mode spectra to make predictions about the experimental observability of this new effect.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nonlinear magnetoplasmons in strongly coupled Yukawa plasmas

The existence of plasma oscillations at multiples of the magnetoplasmon frequency in a strongly coupled two-dimensional magnetized Yukawa plasma is reported, based on extensive molecular dynamics simulations. These modes are the analogues of Bernstein modes which are renormalized by strong interparticle correlations. Their properties are theoretically explained by a dielectric function incorporating the combined effect of a magnetic field, strong correlations and finite temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum simulations of strongly coupled quark-gluon plasma

A strongly coupled quark-gluon plasma (QGP) of heavy constituent quasiparticles is studied by a path-integral Monte-Carlo method, which improves the corresponding classical simulations by extending them to the quantum regime. It is shown that this method is able to reproduce the lattice equation of state and also yields valuable insight into the internal structure of the QGP. The results indicate that the QGP reveals liquid-like rather than gas-like properties. At temperatures just above the critical one it was found that bound quark-antiquark states still survive. These states are bound by effective string-like forces. Quantum effects turned out to be of prime importance in these simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Time-dependent second Born calculations for model atoms and molecules in strong laser fields

Using the finite-element discrete variable representation of the nonequilibrium Green's function (NEGF) we extend previous work [K.~Balzer et al., Phys. Rev. A \textbf{81}, 022510 (2010)] to nonequilibrium situations and compute---from the two-time Schwinger-Keldysh-Kadanoff-Baym equations---the response of the helium atom and the heteronuclear molecule lithium hydride to laser fields in the uv and xuv regime. In particular, by comparing the one-electron density and the dipole moment to time-dependent Hartree-Fock results on one hand and the full solution of the time-dependent Schrödinger equation on the other hand, we demonstrate that the time-dependent second Born approximation carries valuable information about electron-electron correlation effects. Also, we outline an efficient distributed memory concept which enables a parallel and well scalable algorithm for computing the NEGF in the two-time domain.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Two-photon ionization of Helium studied with the multiconfigurational time-dependent Hartree-Fock method

The multiconfigurational time-dependent Hartree-Fock method (MCTDHF) is applied for simulations of the two-photon ionization of Helium. We present results for the single- and double ionization from the groundstate for photon energies in the non-sequential regime, and compare them to direct solutions of the Schrödinger equation using the time-dependent (full) Configuration Interaction method (TDCI). We find that the single-ionization is accurately reproduced by MCTDHF, whereas the double ionization results correctly capture the main trends of TDCI.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Efficient grid-based method in nonequilibrium Green's function calculations. Application to model atoms and molecules

We propose and apply the finite-element discrete variable representation to express the nonequilibrium Green's function for strongly inhomogeneous quantum systems. This method is highly favorable against a general basis approach with regard to numerical complexity, memory resources, and computation time. Its flexibility also allows for an accurate representation of spatially extended hamiltonians, and thus opens the way towards a direct solution of the two-time Schwinger/Keldysh/Kadanoff-Baym equations on spatial grids, including e.g. the description of highly excited states in atoms. As first benchmarks, we compute and characterize, in Hartree-Fock and second Born approximation, the ground states of the He atom, the H$_2$ molecule and the LiH molecule in one spatial dimension. Thereby, the ground-state/binding energies, densities and bond-lengths are compared with the direct solution of the time-dependent Schrödinger equation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Equation of state of strongly coupled quark--gluon plasma -- Path integral Monte Carlo results

A strongly coupled plasma of quark and gluon quasiparticles at temperatures from $ 1.1 T_c$ to $3 T_c$ is studied by path integral Monte Carlo simulations. This method extends previous classical nonrelativistic simulations based on a color Coulomb interaction to the quantum regime. We present the equation of state and find good agreement with lattice results. Further, pair distribution functions and color correlation functions are computed indicating strong correlations and liquid-like behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fast electron generation by Coulomb scattering on spatially correlated ions in a strong laser field

Electrons colliding with spatially fixed ions in strong laser fields are investigated by solving the time-dependent Schrödinger equation. Considering first simple one-dimensional model systems, the mechanisms and energy spectra of fast electrons are analyzed, starting from collisions on a single ion. By using these electrons as projectiles for a second and third collision, the maximum possible energy obtained can be significantly increased. We then generalize the analysis to 2D systems where additional angular degrees of freedom lead to a drastic loss of efficiency. This problem can be overcome by introducing external confinements, which allow to focus the electrons and increase the intensity of high-energy electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

How spherical plasma crystals form

The correlation buildup and the formation dynamics of the shell structure in a spherically confined one-component plasma are studied. Using Langevin dynamics simulations the relaxation processes and characteristic time scales and their dependence on the pair interaction and dissipation in the plasma are investigated. While in systems with Coulomb interaction (e.g. trapped ions) in a harmonic confinement shell formation starts at the plasma edge and proceeds inward, this trend is significantly weakened for dusty plasmas with Yukawa interaction. With a suitable change of the confinement conditions the crystallization scenario can be externally controlled.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Interacting Bosons beyond the Gross-Pitaevskii Mean-Field

Systems consisting of cold interacting bosons show interesting collective phenomena such as Bose-Einstein condensation or superfluidity and are currently studied in condensed matter and atomic physics. Of particular interest are nonideal bosons which exhibit strong correlation and spatial localization effects. Here we analyse the ground-state of a two-dimensional Bose system with a Hartree-Fock type approximation that was first introduced by Romanovsky et al. [Phys. Rev. Lett. {\bf 93}, 230405 (2004)]. We apply this method to a one dimensional system of charged bosons and analyze the behaviour at strong coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Is diffusion anomalous in two-dimensional Yukawa liquids?

There have recently been many predictions of "superdiffusion" in two-dimensional strongly coupled Yukawa systems, both by computer simulations and in dusty plasma experiments, with substantially varying diffusion exponents. Here we show that the results crucially depend on the strength of dissipation and the time instant of the measurement. For sufficiently large friction even subdiffusion is possible. However, there are strong indications that, in the long-time limit, anomalous diffusion vanishes for dissipative as well as for frictionless systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nonequilibrium Green function approach to photoionization processes in atoms

We present a quantum kinetic approach for the time-resolved description of many-body effects in photoionization processes in atoms. The method is based on the non-equilibrium Green functions formalism and solves the Keldysh/Kadanoff-Baym equations in second Born approximation. An approximation scheme is introduced and discussed, which provides a complete single-particle description of the continuum, while the atom is treated fully correlated.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quantum breathing mode of trapped bosons and fermions at arbitrary coupling

An analysis of the quantum breathing behavior of few-particle Coulomb systems in one- and two-dimensional harmonic traps is presented. We report the existence of \emph{two independent breathing modes} and present exact numerical results for two particles at any coupling strength which smoothly connect the two known limits of an \emph {ideal} quantum and a strongly coupled \emph {classical} system. Substantial differences in the breathing frequency of two-dimensional fermions and bosons are observed which may be used as a sensitive experimental tool to probe confined interacting quantum systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quantum kinetic approach to time-resolved photoionization of atoms

Theoretical approaches to the photoionization of few-electron atoms are discussed. These include nonequilibrium Greens functions and wave function based approaches. In particular, the Multiconfiguration Time-Dependent Hartree-Fock method is discussed and applied to a model one-dimensional atom with four electrons. We compute ground state energies and the time-dependent photoionization by the field a strong laser pulse with two different frequencies in the ultraviolet.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Shell Structure of Confined Charges at Strong Coupling

A theoretical description of shell structure for charged particles in a harmonic trap is explored at strong coupling conditions of $Γ$ = 50 and 100. The theory is based on an extension of the hypernetted chain approximation to confined systems plus a phenomenological representation of associated bridge functions. Predictions are compared to corresponding Monte Carlo simulations and quantitative agreement for the radial density profile is obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Structure and Phase transitions of Yukawa balls

In this review, an overview of structural properties and phase transitions in finite spherical dusty (complex) plasma crystals -- so-called Yukawa balls -- is given. These novel kinds of Wigner crystals can be directly analyzed experimentally with video cameras. The experiments clearly reveal a shell structure and allow to determine the shell populations, to observe metastable states and transitions between configurations as well as phase transitions. The experimental observations of the static properties are well explained by a rather simple theoretical model which treats the dust particles as being confined by a parabolic potential and interacting via an isotropic Yukawa pair potential. The excitation properties of the Yukawa balls such as normal modes and the dynamic behavior, including the time-dependent formation of the crystal requires, in addition, to include the effect of friction between the dust particles and the neutral gas. Aside from first-principle molecular dynamics and Monte Carlo simulations several analytical approaches are reviewed which include shell models and a continuum theory. A summary of recent results and theory-experiment comparisons is given and questions for future research activities are outlined.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Theoretical Description of Coulomb Balls - Fluid Phase

A theoretical description for the radial density profile of a finite number of identical charged particles confined in a harmonic trap is developed for application over a wide range of Coulomb coupling (or, equivalently, temperatures) and particle numbers. A simple mean field approximation neglecting correlations yields a density profile which is monotonically decreasing with radius for all temperatures, in contrast to molecular dynamics simulations and experiments showing shell structure at lower temperatures. A more complete theoretical description including charge correlations is developed here by an extension of the hypernetted chain approximation, developed for bulk fluids, to the confined charges. The results reproduce all of the qualitative features observed in molecular dynamics simulations and experiments. These predictions are then tested quantitatively by comparison with new benchmark Monte Carlo simulations. Quantitative accuracy of the theory is obtained for the selected conditions by correcting the hypernetted chain approximation with a representation for the associated bridge functions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Towards a theory of attosecond transient recorder

Laser assisted photoemission by a chirped subfemtosecond extreme ultraviolet (XUV) pulse is considered within an exactly solvable quantum-mechanical model. Special emphasis is given to the energy dependence of photoexcitation cross-section. The streaked spectra are analyzed within the classical picture of initial time-momentum distribution r_ini(p,t) of photoelectrons mapped to the final energy scale. The actual time-momentum distribution in the absence of the probe laser field is shown to be a poor choice for r_ini, and a more adequate ansatz is suggested. The semiclassical theory offers a simple practically useful approximation for streaked spectra. Its limitations for sufficiently long chirped XUV pulses are established.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Angular distributions of atomic photoelectrons produced in the UV and XUV regimes

We present angular distributions of photoelectrons of atomic model systems excited by intense linearly polarized laser pulses in the VUV- and XUV-regime. We solve the multi-dimensional time-dependent Schrödinger equation for one particle on large spatial grids and investigate the direction dependence of the ionized electrons for isotropic s-states as well as p-states. Although the ponderomotive potential is small compared to the binding energy of the initially bound electron and the photon energy of the exciting laser field, richly structured photoelectron angular distributions are found which sensitively depend on the laser frequency and intensity as well as on the number of absorbed photons. The occuring shapes are explained in terms of scattering mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Effective interaction potential and superfluid-solid transition of spatially indirect excitons

Using an adiabatic approximation we derive an effective interaction potentially for spatially indirect excitons. Using this potential and path integral Monte Carlo simulations we study exciton crystllization and the quantum melting phase transition in a macroscopic system of 2D excitons. Furthermore, the superfluid fraction is calculated as a function of density and shown to vanish upon crystallization. We show that the commonly used dipole model fails to correctly describe indirect excitons in quantum well structures.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Nonequilibrium Green's functions approach to strongly correlated few-electron quantum dots

The effect of electron-electron scattering on the equilibrium properties of few-electron quantum dots is investigated by means of nonequilibrium Green's functions theory. The ground and equilibrium state is self-consistently computed from the Matsubara (imaginary time) Green's function for the spatially inhomogeneous quantum dot system whose constituent charge carriers are treated as spin-polarized. To include correlations, the Dyson equation is solved, starting from a Hartree-Fock reference state, within a conserving (second order) self-energy approximation where direct and exchange contributions to the electron-electron interaction are included on the same footing. We present results for the zero and finite temperature charge carrier density, the orbital-resolved distribution functions and the self-consistent total energies and spectral functions for isotropic, two-dimensional parabolic confinement as well as for the limit of large anisotropy--quasi-one-dimensional entrapment. For the considered quantum dots with N=2, 3 and 6 electrons, the analysis comprises the crossover from Fermi gas/liquid (at large carrier density) to Wigner molecule or crystal behavior (in the low-density limit).

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Structures of quantum 2D electron-hole plasmas

We investigate structures of 2D quantum electron-hole (e-h) plasmas by the direct path integral Monte Carlo method (PIMC) in a wide range of temperature, density and hole-to-electron mass ratio. Our simulation includes a region of appearance and decay of the bound states (excitons and biexcitons), the Mott transition from the neutral e-h plasma to metallic-like clusters, formation from clusters the hexatic-like liquid and formation of the crystal-like lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Thermodynamic properties and electrical conductivity of strongly correlated plasma media

We study thermodynamic properties and the electrical conductivity of dense hydrogen and deuterium using three methods: classical reactive Monte Carlo (REMC), direct path integral Monte Carlo (PIMC) and a quantum dynamics method in the Wigner representation of quantum mechanics. We report the calculation of the deuterium compression quasi-isentrope in good agreement with experiments. We also solve the Wigner-Liouville equation of dense degenerate hydrogen calculating the initial equilibrium state by the PIMC method. The obtained particle trajectories determine the momentum-momentum correlation functions and the electrical conductivity and are compared with available theories and simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Kadanoff-Baym equations and non-Markovian Boltzmann equation in generalized T-matrix approximation

A recently developed method for incorporating initial binary correlations into the Kadanoff-Baym equations (KBE) is used to derive a generalized T-matrix approximation for the self-energies. It is shown that the T-matrix obtains additional contributions arising from initial correlations. Using these results and taking the time-diagonal limit of the KBE, a generalized quantum kinetic equation in binary collision approximation is derived. This equation is a far-reaching generalization of Boltzmann-type kinetic equations: it selfconsistently includes memory effects (retardation, off-shell T-matrices) as well as many-particle effects (damping, in-medium T-Matrices) and spin-statistics effects (Pauli-blocking).

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Plasmon dispersion of a weakly degenerate nonideal one-component plasma

Classical Molecular Dynamics simulations (MD) for a one-component weakly degenerate plasma are presented. Using an effective quantum pair potential (Kelbg potential), the dynamic structure factor and the dispersion of Langmuir waves are computed. The influence of the coupling strength $Γ$ and degree of degeneracy $ρΛ^3$ on these properties is discussed. The results are compared with predictions of mean-field theories.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Two-dimensional electron crystals in single and double layers

We present results of Monte-Carlo simulations for finite 2D single and bilayer systems. Strong Coulomb correlations lead to arrangement of particles in configurations resembling a crystal lattice. For binary layers, there exists a particularly rich variety of lattice symmetries which depend on the interlayer separation $d$. We demonstrate that in these mesoscopic lattices there exist two fundamental types of ordering: radial and orientational. The dependence of the melting temperature on $d$ is analyzed, and a stabilization of the crystal compared to a single layer is found.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Quantum Kinetic Theory for Laser Plasmas. Dynamical Screening in Strong Fields

A quantum kinetic theory for correlated charged-particle systems in strong time-dependent electromagnetic fields is developed. Our approach is based on a systematic gauge-invariant nonequilibrium Green's functions formulation. Extending our previous analysis (Kremp et al., Phys. Rev. E vol. 60, p. 4725 (1999) we concentrate on the selfconsistent treatment of dynamical screening and electromagnetic fields which is applicable to arbitrary nonequilibrium situations. The resulting kinetic equation generalizes previous results to quantum plasmas with full dynamical screening and includes many-body effects. It is, in particular, applicable to the interaction of dense plasmas with strong electromagnetic fields, including laser fields and x-rays. Furthermore, results for the modification of the plasma screening and the longitudinal field fluctuations due to the electromagnetic field are presented.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Short time dynamics with initial correlations

The short-time dynamics of correlated systems is strongly influenced by initial correlations giving rise to an additional collision integral in the non-Markovian kinetic equation. Exact cancellation of the two integrals is found if the initial state is thermal equilibrium which is an important consistency criterion. Analytical results are given for the time evolution of the correlation energy which are confirmed by comparisons with molecular dynamics simulations (MD).

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Thermodynamic Properties of Strongly Degenerate Interacting Fermi Systems

A numerical approach is presented which allows to calculate the equilibrium properties of Fermi systems which are both, strongly coupled and strongly degenerate. Based on a novel path integral representation of the many-particle density operator, quantum and spin effects are included. As an illustration, results for the pressure and energy of an electron-proton plasma are presented.

0 citations0 reviewsNo code linkedOpen access