ترغب بنشر مسار تعليمي؟ اضغط هنا

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

151   0   0.0 ( 0 )
 نشر من قبل Michael Bonitz
 تاريخ النشر 2014
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

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.



قيم البحث

اقرأ أيضاً

We present extensive new textit{ab intio} path integral Monte Carlo results for the momentum distribution function $n(mathbf{k})$ of the uniform electron gas (UEG) in the warm dense matter (WDM) regime over a broad range of densities and temperatures . This allows us to study the nontrivial exchange--correlation induced increase of low-momentum states around the Fermi temperature, and to investigate its connection to the related lowering of the kinetic energy compared to the ideal Fermi gas. In addition, we investigate the impact of quantum statistics on both $n(mathbf{k})$ and the off-diagonal density matrix in coordinate space, and find that it cannot be neglected even in the strongly coupled electron liquid regime. Our results were derived without any nodal constraints, and thus constitute a benchmark for other methods and approximations.
Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculatio n ultimately requires the ill-defined analytical continuation from the imaginary- to the real-frequency domain. It was recently proposed [Phys. Rev. B 99, 035120 (2019)] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy $Sigma(omega)$ of the doped 32x32 cyclic square-lattice Hubbard model in a non-trivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.
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 simulat ions 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 $mu=175$ MeV, i.e. glueballs, while meson-like bound states are not found.
84 - Mamikon Gulian , Haobo Yang , 2017
Fractional derivatives are nonlocal differential operators of real order that often appear in models of anomalous diffusion and a variety of nonlocal phenomena. Recently, a version of the Schrodinger Equation containing a fractional Laplacian has bee n proposed. In this work, we develop a Fractional Path Integral Monte Carlo algorithm that can be used to study the finite temperature behavior of the time-independent Fractional Schrodinger Equation for a variety of potentials. In so doing, we derive an analytic form for the finite temperature fractional free particle density matrix and demonstrate how it can be sampled to acquire new sets of particle positions. We employ this algorithm to simulate both the free particle and $^{4}$He (Aziz) Hamiltonians. We find that the fractional Laplacian strongly encourages particle delocalization, even in the presence of interactions, suggesting that fractional Hamiltonians may manifest atypical forms of condensation. Our work opens the door to studying fractional Hamiltonians with arbitrarily complex potentials that escape analytical solutions.
Quantum Monte Carlo belongs to the most accurate simulation techniques for quantum many-particle systems. However, for fermions, these simulations are hampered by the sign problem that prohibits simulations in the regime of strong degeneracy. The sit uation changed with the development of configuration path integral Monte Carlo (CPIMC) by Schoof textit{et al.} [T. Schoof textit{et al.}, Contrib. Plasma Phys. textbf{51}, 687 (2011)] that allowed for the first textit{ab initio} simulations for dense quantum plasmas. CPIMC also has a sign problem that occurs when the density is lowered, i.e. in a parameter range that is complementary to traditional QMC formulated in coordinate space. Thus, CPIMC simulations for the warm dense electron gas are limited to small values of the Brueckner parameter -- the ratio of the interparticle distance to the Bohr radius -- $r_s=bar{r}/a_B lesssim 1$. In order to reach the regime of stronger coupling (lower density) with CPIMC, here we investigate additional restrictions on the Monte Carlo procedure. In particular, we introduce two differe
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا