共查询到20条相似文献，搜索用时 140 毫秒
1.
The uniform electron gas (UEG) is one of the key models for the understanding of warm dense matter—an exotic, highly compressed state of matter between solid and plasma phases. The difficulty in modelling the UEG arises from the need to simultaneously account for Coulomb correlations, quantum effects, and exchange effects, as well as finite temperature. The most accurate results so far were obtained from quantum Monte Carlo (QMC) simulations with a variety of representations. However, QMC for electrons is hampered by the fermion sign problem. Here, we present results from a novel fermionic propagator path integral Monte Carlo in the restricted grand canonical ensemble. The ab initio simulation results for the spin-resolved pair distribution functions and static structure factor are reported for two isotherms (T in the units of the Fermi temperature). Furthermore, we combine the results from the linear response theory in the Singwi-Tosi-Land-Sjölander scheme with the QMC data to remove finite-size errors in the interaction energy. We present a new corrected parametrization for the interaction energy and the exchange–correlation free energy in the thermodynamic limit, and benchmark our results against the restricted path integral Monte Carlo by Brown et al. [Phys. Rev. Lett. 110 , 146405 (2013)] and configuration path integral Monte Carlo/permutation-blocking path integral Monte Carlo by Dornheim et al. [Phys. Rev. Lett. 117 , 115701 (2016)]. 相似文献
2.
3.
4.
5.
6.
We analyse the simulation of strongly degenerate electrons at finite temperature using the recently introduced permutation blocking path integral Monte Carlo (PB‐PIMC) method [T. Dornheim et al., New J. Phys. 17 , 073017 (2015)]. As a representative example, we consider electrons in a harmonic confinement and carry out simulations for up to P = 2000 so‐called imaginary‐time propagators – an important convergence parameter within the PIMC formalism. This allows us to study the P‐dependence of different observables of the configuration space in the Monte Carlo simulations and of the fermion sign problem. We find a surprisingly persisting effect of the permutation blocking for large P, which is explained by comparing different length scales. Finally, we touch upon the uniform electron gas in the warm dense matter regime. 相似文献
7.
Francis J. Alexander Gregory L. Eyink Juan M. Restrepo 《Journal of statistical physics》2005,119(5-6):1331-1345
Asbstract By casting stochastic optimal estimation of time series in path integral form, one can apply analytical and computational techniques of equilibrium statistical mechanics. In particular, one can use standard or accelerated Monte Carlo methods for smoothing, filtering and/or prediction. Here we demonstrate the applicability and efficiency of generalized (nonlocal) hybrid Monte Carlo and multigrid methods applied to optimal estimation, specifically smoothing. We test these methods on a stochastic diffusion dynamics in a bistable potential. This particular problem has been chosen to illustrate the speedup due to the nonlocal sampling technique, and because there is an available optimal solution which can be used to validate the solution via the hybrid Monte Carlo strategy. In addition to showing that the nonlocal hybrid Monte Carlo is statistically accurate, we demonstrate a significant speedup compared with other strategies, thus making it a practical alternative to smoothing/filtering and data assimilation on problems with state vectors of fairly large dimensions, as well as a large total number of time steps. 相似文献
8.
9.
10.
Naeem A. Tahir Alexander Shutov Antonio R. Piriz Paul Neumayer Igor V. Lomonosov Vincent Bagnoud Sofia A. Piriz 《等离子体物理论文集》2019,59(4-5)
This paper presents detailed 2D hydrodynamic simulations of implosion of a multi‐layered cylindrical target that is driven by an intense uranium beam. The target is comprised of a thick, high‐Z, high‐ρ cylindrical shell that encloses a sample material (Fe in the present case). Two options have been used for the focal spot geometry: an annular form and a circular form. The purpose of this work is to show that an intense heavy‐ion beam can induce the extreme physical conditions in the sample material similar to those that exist in the planetary cores. In this study, we use parameters of the beam that will be generated at the Facility for Antiprotons and Ion Research (FAIR), Darmstadt, in a few years' time. Production of these high‐energy‐density (HED) samples will allow us to study planetary physics in the laboratory. It is to be noted that planetary physics research is an important part of the FAIR HED physics program. A dedicated experiment named LAboratory PLAnetary Sciences (LAPLAS) has been proposed for this purpose. These simulations show that in such experiments an Fe sample can be imploded to the Earth's core conditions and to those in more massive rocky planets called Super‐Earths. Similarly, implosion of hydrogen and water samples will generate the core conditions of solar and extrasolar hydrogen‐rich gas giants and water‐rich icy planets, respectively. The LAPLAS experiments will thus provide very valuable information on the equation of state and transport properties of matter under extreme physical conditions, which will help scientists understand the structure and evolution of the planets in our solar system as well as of the extrasolar planets. 相似文献
11.
A novel path integral Monte Carlo (PIMC) approach for correlated many‐particle systems with arbitrary pair interaction in continuous space at low temperatures is presented. It is based on a representation of the N ‐particle density operator in a basis of (anti‐)symmetrized N ‐particle states (configurations of occupation numbers). The path integral is transformed into a sum over trajectories with the same topology and, finally, the limit of M → ∞, where M is the number of high‐temperature factors, is analytically performed. This yields exact expressions for the thermodynamic quantities and allows to perform efficient simulations for fermions at low temperature and weak to moderate coupling. Our method is expected to be applicable to dense quantum plasmas in the regime of strong degeneracy where conventional PIMC fails due to the fermion sign problem (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim) 相似文献
12.
We investigate a pump‐probe X‐ray Thomson scattering (XRTS) experiment that might be carried out at a free electron laser facility to study warm‐to‐hot states of dense matter. Ultrashort and intense X‐ray pulses with different energies (1,560–1,830 eV) heat a 1 µm thick Al target isochorically and create homogeneous and uncompressed warm‐to‐hot states of dense matter. A second pulse with variable delay probes this heated state via XRTS. The X‐ray laser–target interaction is modelled within radiation‐hydrodynamic simulations applying the HELIOS‐CR code. The HELIOS‐CR results qualitatively agree with Monte‐Carlo simulations, where the laser pulse absorption is simulated based on a uniform random sequence of events. The electron feature in the simultaneously observed X‐ray scattering spectrum is a function of the degree of ionization and the target temperature. Therefore, the temporal evolution of the plasmon peak measures the ionization dynamics on ultra‐short time scales. The XRTS spectrum is calculated based on the Chihara formula utilizing the Born‐Mermin approximation for the free electron dynamic structure factor. The proposed experiment will reveal important details of the ionization dynamics on ultra‐short time scales as well as of the relaxation on ps time scales. 相似文献
13.
Paul Hamann Jan Vorberger Tobias Dornheim Zhandos A. Moldabekov Michael Bonitz 《等离子体物理论文集》2020,60(10):e202000147
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, and 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 ab initio results for S(q, ω) of the homogeneous warm dense electron gas were obtained from path integral Monte Carlo simulations (Dornheim et al., Phys. Rev. Lett., 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 zeroes of the dielectric function can solve this problem for WDM conditions. 相似文献
14.
15.
温稠密物质的物性参数在惯性约束聚变能源、Z箍缩等高能量密度物理领域的实验结果分析和物理过程数值模拟等方面有着重要的应用价值.本文应用部分电离等离子体模型,在理想自由能的基础上考虑了库仑相互作用、排斥体积作用和极化作用等非理想特性,开展了温稠密等离子体物态方程和电离平衡的研究.计算了温稠密铝等离子体的压强等物态方程数据和在密度为1.0×10^-4-3.0 g/cm^3,温度为1.0×10^4-3.0×10^4 K范围内的粒子组分.计算结果显示,铝等离子体的平均电离度在临界密度区域内随着密度的增加而突然增大.根据非理想Saha方程中有效电离能这一关键参数,分析了铝等离子体平均电离度在临界密度区域内随密度迅速增大的现象. 相似文献
16.
ZHAO Xing-Wen CHENG Xin-Lu 《理论物理通讯》2008,50(10):987-989
Path integral Monte Carlo （PIMC） simulations are a powerful computational method to study interacting quantum systems at finite temperatures. In this work, PIMC has been applied to study the finite size effect of the simulated systems of ^4He. We determine the energy as a function of temperature at saturated-vapor-pressure （SVP） conditions in the temperature range of T ∈ [1.0 K,4.0 K], and the equation of state （EOS） in the grmmd state For systems consisted of 32, 64 and 128 ^4He atoms, respectively, We find that the energy at SVP is influenced significantly by the size of the simulated system in the temperature range of T ∈ [2.1 K, 3.0 K] and the larger the system is, the better results are obtained in comparison with the experimental values; while the EOS appeared to be unrelated to it. 相似文献
17.
A new wave packet molecular dynamics–density functional theory (WPMD‐DFT) method is proposed for atomistic simulations of non‐ideal plasma and warm dense matter. The method is based on the WPMD approach, where the electronic exchange and correlation effects are treated using an additional energy term taken from DFT. This term is calculated by integration over the mesh values of the wave packet density. The local density approximation is implemented so far. WPMD‐DFT is meant as a replacement for the anti‐symmetrized WPMD (AWPMD) method which is more time consuming and lacks electron correlation. In this paper, we compare the results obtained by WPMD‐DFT, WPMD, AWPMD, classical molecular dynamics, and path integral Monte Carlo methods for the internal energy of the hydrogen plasma in the temperature range 10–50 kK and electron number density from 1020 to 1024 cm?3. We also demonstrate the ability to handle the simultaneous dynamics of electrons and ions by calculating the electron–ion temperature relaxation. The scalability of the WPMD‐DFT method with the number of electrons is shown for implementations in central processing unit and graphical processing unit. 相似文献
18.
Naeem A. Tahir Paul Neumayer Alexander Shutov Antonio R. Piriz Igor V. Lomonosov Vincent Bagnoud Sofia A. Piriz Claude Deutsch 《等离子体物理论文集》2019,59(4-5)
This paper describes an experiment design based on numerical simulations to measure the equation‐of‐state properties of high‐energy‐density (HED) matter using intense particle beams. The simulations are performed using a 2D hydrodynamic computer code, BIG2, while the beam parameters are considered to match the Facility for Antiprotons and Ion Research beam. This study has shown that in such experiments one can generate different phases of HED lead. Similar calculations are planned for other materials. 相似文献
19.
Path integral Monte Carlo (PIMC) simulations are a powerful computational method to study interacting quantum systems at finite
temperatures. In this work, PIMC has been applied to study the finite size
effect of the simulated systems of
4He. We determine the energy as a function of temperature at saturated-vapor-pressure (SVP) conditions in the temperature range of T∈[1.0 K, 4.0 K], and the equation of state (EOS) in the ground state for systems consisted of 32, 64 and 128
4He atoms, respectively. We find that the energy at SVP is influenced significantly by the size of the simulated system in the temperature range of T∈[2.1 K, 3.0 K] and the larger the system is, the better results are obtained in comparison with the experimental values; while the EOS appeared to be unrelated to it. 相似文献
20.
利用路径积分蒙特卡罗(PIMC)模拟方法研究了由256个4He原子组成的系统在不同条件下的对关联函数g(r).在标准大气压下(SVP),分别在1.38K,1.67K,2.50K和4.24K等四个温度时对系统不同势能的情况的对关联函数进行了计算;另外还讨论了系统的温度和密度对对关联函数的影响.结果发现,系统的势能变化及密度变化都或大或小影响对关联函数,而温度对其影响不是很大.虽然PIMC方法存在有限尺度效应,但是由128和256个4He原子组成系统的对关联函数存在极小的差异,因此可以得到:随着模拟系统尺度的增大,所模拟的有限尺度的系统已逐渐接近于真实系统.所以所得的256原子系统的性质可近似说明液态系统性质. 相似文献