-
Phase space path integral representation of the dynamic structure factor. Monte Carlo simulation of strongly correlated soft-sphere fermions
Authors:
V. S. Filinov,
P. R. Levashov,
A. S. Larkin
Abstract:
The dynamic structure factor (DSF) is a mathematical function that contains information about inter-particle correlations and their time evolution. Mostly the classical molecular dynamics is used to calculate the DSF of the classical systems. On the contrary this article deals with quantum systems and the quantum dynamic structure factor. The Wigner formulation of quantum mechanics was used to der…
▽ More
The dynamic structure factor (DSF) is a mathematical function that contains information about inter-particle correlations and their time evolution. Mostly the classical molecular dynamics is used to calculate the DSF of the classical systems. On the contrary this article deals with quantum systems and the quantum dynamic structure factor. The Wigner formulation of quantum mechanics was used to derive the path integral representation of the DSF, which is based on the Wiener-Khinchin theorem showing relation of the the power spectrum of a random paths to their correlation function. The $3 {\rm D} $ quantum system of strongly correlated soft-sphere fermions was considered as an interesting physical example.
The developed Wigner path integral Monte Carlo (WPIMC) approach has been developed to calculate the spin--resolved DSFs, the radial distribution functions (RDFs) and other thermodynamic functions in a wide range of density and temperatures. The physical meaning of the peaks arising on the RDFs and DSF have been analyzed and explained by the manifestation of the interference effects of the exchange and interparticle interactions and, as well as, the wave interference between multiple-scattering. This phenomenon in the system of the soft-sphere scatterers may to be the precursor effect of the Anderson localization, which finds its origin in the wave interference between multiple-scattering paths.
△ Less
Submitted 26 June, 2024;
originally announced July 2024.
-
Pressure of Coulomb systems with volume-dependent long-range potentials
Authors:
A. S. Onegin,
G. S. Demyanov,
P. R. Levashov
Abstract:
In this work, we consider the pressure of Coulomb systems, in which particles interact via a volume-dependent potential (in particular, the Ewald potential). We confirm that the expression for virial pressure should be corrected in this case. We show that the corrected virial pressure coincides with the formula obtained by differentiation of free energy if the potential energy is a homogeneous fun…
▽ More
In this work, we consider the pressure of Coulomb systems, in which particles interact via a volume-dependent potential (in particular, the Ewald potential). We confirm that the expression for virial pressure should be corrected in this case. We show that the corrected virial pressure coincides with the formula obtained by differentiation of free energy if the potential energy is a homogeneous function of particle coordinates and a cell length. As a consequence, we find out that the expression for pressure in the recent paper by J. Liang \textit{et al.} [\href{https://doi.org/10.1063/5.0107140}{J. Chem. Phys. \textbf{157}, 144102 (2022)}] is incorrect.
△ Less
Submitted 11 September, 2023;
originally announced September 2023.
-
Exchange--correlation bound states of the triplet soft--sphere fermions by the path integral Monte Carlo simulations
Authors:
V. S. Filinov,
R. A. Syrovatka,
P. R. Levashov
Abstract:
Path integral Monte Carlo simulations in the Wigner approach to quantum mechanics has been applied to calculate momentum and spin--resolved radial distribution functions of the strongly correlated soft--sphere quantum fermions. The obtained spin--resolved radial distribution functions demonstrate arising triplet clusters of fermions, that is the consequence of the interference of exchange and inte…
▽ More
Path integral Monte Carlo simulations in the Wigner approach to quantum mechanics has been applied to calculate momentum and spin--resolved radial distribution functions of the strongly correlated soft--sphere quantum fermions. The obtained spin--resolved radial distribution functions demonstrate arising triplet clusters of fermions, that is the consequence of the interference of exchange and interparticle interactions. The semiclassical analysis in the framework of the Bohr--Sommerfeld quantization condition applied to the potential of the mean force corresponding to the same--spin radial distribution functions allows to detect exchange--correlation bound states in triplet clusters and to estimate corresponding averaged energy levels. The obtained momentum distribution functions demonstrate the narrow sharp separated peaks corresponding to bound states and disturbing the Maxwellian distribution.
△ Less
Submitted 30 May, 2023;
originally announced June 2023.
-
Momentum distribution functions and pair correlation functions of unpolarized uniform electron gas in warm dense matter regime
Authors:
A. S. Larkin,
V. S. Filinov,
P. R. Levashov
Abstract:
In this paper we continued our research of the uniform electron gas, using the single--momentum path integral Monte Carlo method, and studied the momentum distribution functions and the pair correlation functions in the warm dense matter regime. We discovered that the single--particle momentum distribution function deviates from the Fermi distribution and forms so-called "quantum tails" at high mo…
▽ More
In this paper we continued our research of the uniform electron gas, using the single--momentum path integral Monte Carlo method, and studied the momentum distribution functions and the pair correlation functions in the warm dense matter regime. We discovered that the single--particle momentum distribution function deviates from the Fermi distribution and forms so-called "quantum tails" at high momenta, if non-ideality is strong enough in both degenerate and non-degenerate cases. This effect is always followed by the appearance of the short--range order on the pair distribution functions and can be explained via the tunneling through the effective potential wells surrounding the electrons. Also we calculated the average kinetic and potential energies in the wide range of states, expanding our previous results significantly.
△ Less
Submitted 23 March, 2022;
originally announced March 2022.
-
Ab initio inspection of thermophysical experiments for zirconium near melting
Authors:
M. A. Paramonov,
D. V. Minakov,
V. B. Fokin,
D. V. Knyazev,
G. S. Demyanov,
P. R. Levashov
Abstract:
We present quantum molecular dynamics calculations of thermophysical properties of solid and liquid zirconium in the vicinity of melting. An overview of available experimental data is also presented. We focus on the analysis of thermal expansion, molar enthalpy, resistivity and normal spectral emissivity of solid and liquid Zr. Possible reasons of discrepancies between the first-principle simulati…
▽ More
We present quantum molecular dynamics calculations of thermophysical properties of solid and liquid zirconium in the vicinity of melting. An overview of available experimental data is also presented. We focus on the analysis of thermal expansion, molar enthalpy, resistivity and normal spectral emissivity of solid and liquid Zr. Possible reasons of discrepancies between the first-principle simulations and experiments are discussed. Our calculations reveal a significant volume change on melting in agreement with electrostatic levitation experiments. Meanwhile, we confirm a low value of enthalpy of fusion obtained in some pulse-heating experiments. Electrical resistivity of solid and liquid Zr is systematically underestimated in our simulations, however the slope of resistivity temperature dependencies agrees with experiment. Our calculations predict almost constant normal spectral emissivity in liquid Zr.
△ Less
Submitted 27 April, 2021; v1 submitted 23 April, 2021;
originally announced April 2021.
-
Effect of the spin-orbit interaction on thermodynamic properties of liquid uranium
Authors:
D. V. Minakov,
M. A. Paramonov,
P. R. Levashov
Abstract:
We present the first quantum molecular dynamics calculation of zero-pressure isobar of solid and liquid uranium that account for spin-orbit coupling. We demonstrate that inclusion of spin-orbit interaction leads to higher degree of the thermal expansion of uranium, especially in the liquid phase. Full accounting of relativistic effects for valence electrons, particularly spin-orbital splitting of…
▽ More
We present the first quantum molecular dynamics calculation of zero-pressure isobar of solid and liquid uranium that account for spin-orbit coupling. We demonstrate that inclusion of spin-orbit interaction leads to higher degree of the thermal expansion of uranium, especially in the liquid phase. Full accounting of relativistic effects for valence electrons, particularly spin-orbital splitting of the 5f band, is substantial for the reproduction of the experimental density of molten uranium at the melting temperature. Influence of the spin-orbit interaction on the thermodynamic properties at high temperatures and pressures is also analyzed.
△ Less
Submitted 11 May, 2020;
originally announced May 2020.
-
Ab initio calculation of thermodynamic, transport, and optical properties of CH$_2$ plastics
Authors:
D. V. Knyazev,
P. R. Levashov
Abstract:
This work covers an ab initio calculation of thermodynamic, transport, and optical properties of plastics of the effective composition CH$_2$ at density 0.954 g/cm$^3$ in the temperature range from 5 kK up to 100 kK. The calculation is based on the quantum molecular dynamics, density functional theory and the Kubo-Greenwood formula. The temperature dependence of the static electrical conductivity…
▽ More
This work covers an ab initio calculation of thermodynamic, transport, and optical properties of plastics of the effective composition CH$_2$ at density 0.954 g/cm$^3$ in the temperature range from 5 kK up to 100 kK. The calculation is based on the quantum molecular dynamics, density functional theory and the Kubo-Greenwood formula. The temperature dependence of the static electrical conductivity $σ(T)$ has a step-like shape: $σ(T)$ grows rapidly for 5 kK <= $T$ <= 10 kK and is almost constant for 20 kK <= $T$ <= 60 kK. The additional analysis based on the investigation of the electron density of states (DOS) is performed. The rapid growth of $σ(T)$ at 5 kK<= $T$ <= 10 kK is connected with the increase of DOS at the electron energy equal to the chemical potential $ε= μ$. The frequency dependence of the dynamic electrical conductivity $σ_1(ω)$ at 5 kK has the distinct non-Drude shape with the peak at $ω\approx 10$ eV. This behavior of $σ_1(ω)$ was explained by the dip at the electron DOS.
△ Less
Submitted 19 July, 2016;
originally announced July 2016.
-
Implementation of kinetics of phase transitions into hydrocode for simulation of laser ablation
Authors:
Mikhail E. Povarnitsyn,
Pavel R. Levashov,
Konstantin V. Khishchenko
Abstract:
We model an interaction of femtosecond laser pulses (800 nm, 100 fs, 1-100 TW/cm$^2$) with metal targets (Al, Au, Cu, and Ni). A detailed analysis of laser-induced phase transitions, melting wave propagation and material decomposition is performed using a thermodynamically complete two-temperature equation of state with stable and metastable phases. Material evaporation from the surface of the tar…
▽ More
We model an interaction of femtosecond laser pulses (800 nm, 100 fs, 1-100 TW/cm$^2$) with metal targets (Al, Au, Cu, and Ni). A detailed analysis of laser-induced phase transitions, melting wave propagation and material decomposition is performed using a thermodynamically complete two-temperature equation of state with stable and metastable phases. Material evaporation from the surface of the target and fast melting wave propagation into the bulk are observed. On rarefaction, the liquid phase becomes metastable and its lifetime is estimated using the theory of homogeneous nucleation. Mechanical fragmentation of the target material at high strain rates is also possible because of void growth and confluence. In our simulation several ablation mechanisms are observed but the major output of the material is found to originate from the metastable liquid state. It can be decomposed either into a liquid-gas mixture near the critical point, or into droplets at high strain rates and negative pressure. The simulation results correlate with available experiments.
△ Less
Submitted 9 November, 2014;
originally announced November 2014.
-
Transport and optical properties of warm dense aluminum in the two-temperature regime: Ab initio calculation and semiempirical approximation
Authors:
D. V. Knyazev,
P. R. Levashov
Abstract:
This work is devoted to the investigation of transport and optical properties of liquid aluminum in the two-temperature case. At first optical properties, static electrical and thermal conductivities were obtained in the \textit{ab initio} calculation. The \textit{ab initio} calculation is based on the quantum molecular dynamics, density functional theory and the Kubo-Greenwood formula. The semiem…
▽ More
This work is devoted to the investigation of transport and optical properties of liquid aluminum in the two-temperature case. At first optical properties, static electrical and thermal conductivities were obtained in the \textit{ab initio} calculation. The \textit{ab initio} calculation is based on the quantum molecular dynamics, density functional theory and the Kubo-Greenwood formula. The semiempirical approximation was constructed based on the results of the \textit{ab initio} caculation. The approximation yields the dependences $σ_{1_\mathrm{DC}}\propto1/T_i^{0.25}$ and $K\propto T_e/T_i^{0.25}$ for the static electrical conductivity and thermal conductivity, respectively. The approximation is valid for liquid aluminum at $ρ=2.70$~g/cm$^3$, 3~kK~$\leq T_i\leq T_e\leq20$~kK. Our results are well described by the Drude model with the effective relaxation time $τ\propto T_i^{-0.25}$. We have compared our results with a number of other models. They are all reduced in the low-temperature limit to the Drude model with different expressions for the relaxation time $τ$. Our results are not consistent with the models in which $τ\propto T_i^{-1}$ and support the models which use the expressions with the slower decrease of the relaxation time.
△ Less
Submitted 10 September, 2014;
originally announced September 2014.
-
Ab initio simulation of complex dielectric function for dense aluminum plasma
Authors:
M. E. Povarnitsyn,
D. V. Knyazev,
P. R. Levashov
Abstract:
We present calculations of frequency-dependent complex dielectric function of dense aluminum plasma by quantum molecular dynamics method for temperatures up to 20 kK. Analysis shows that the dependencies for real and imaginary parts can be interpolated by the Drude formula with two effective parameters: the mean charge of ions and the effective frequency of collisions. The rise of these parameters…
▽ More
We present calculations of frequency-dependent complex dielectric function of dense aluminum plasma by quantum molecular dynamics method for temperatures up to 20 kK. Analysis shows that the dependencies for real and imaginary parts can be interpolated by the Drude formula with two effective parameters: the mean charge of ions and the effective frequency of collisions. The rise of these parameters with temperature deviates from simple theoretical predictions.
△ Less
Submitted 21 October, 2013;
originally announced October 2013.
-
\textit{Ab initio} calculation of transport and optical properties of aluminum: influence of simulation parameters
Authors:
D. V. Knyazev,
P. R. Levashov
Abstract:
This work is devoted to the \textit{ab initio} calculation of transport and optical properties of aluminum. The calculation is based on the quantum molecular dynamics simulation, density functional theory and the Kubo-Greenwood formula. Mainly the calculations are performed for liquid aluminum at near-normal densities for the temperatures from melting up to 20000 K. The results on dynamic electric…
▽ More
This work is devoted to the \textit{ab initio} calculation of transport and optical properties of aluminum. The calculation is based on the quantum molecular dynamics simulation, density functional theory and the Kubo-Greenwood formula. Mainly the calculations are performed for liquid aluminum at near-normal densities for the temperatures from melting up to 20000 K. The results on dynamic electrical conductivity, static electrical conductivity and thermal conductivity are obtained and compared with available reference and experimental data and the calculations of other authors. The influence of the technical parameters on the results is investigated in detail. The error of static electrical conductivity calculation is estimated to be about 20%; more accurate results require bigger number of atoms.
△ Less
Submitted 15 October, 2013;
originally announced October 2013.
-
Hole crystallization in semiconductors
Authors:
M. Bonitz,
V. S. Filinov,
V. E. Fortov,
P. R. Levashov,
H. Fehske
Abstract:
When electrons in a solid are excited to a higher energy band they leave behind a vacancy (hole) in the original band which behaves like a positively charged particle. Here we predict that holes can spontaneously order into a regular lattice in semiconductors with sufficiently flat valence bands. The critical hole to electron effective mass ratio required for this phase transition is found to be…
▽ More
When electrons in a solid are excited to a higher energy band they leave behind a vacancy (hole) in the original band which behaves like a positively charged particle. Here we predict that holes can spontaneously order into a regular lattice in semiconductors with sufficiently flat valence bands. The critical hole to electron effective mass ratio required for this phase transition is found to be of the order of 80.
△ Less
Submitted 16 December, 2005;
originally announced December 2005.
-
Crystallization in two-component Coulomb systems
Authors:
M. Bonitz,
V. S. Filinov,
V. E. Fortov,
P. R. Levashov,
H. Fehske
Abstract:
The analysis of Coulomb crystallization is extended from one-component to two-component plasmas. Critical parameters for the existence of Coulomb crystals are derived for both classical and quantum crystals. In the latter case, a critical mass ratio of the two charged components is found which is of the order of 80. Thus, holes in semiconductors with sufficiently flat valence bands are predicted…
▽ More
The analysis of Coulomb crystallization is extended from one-component to two-component plasmas. Critical parameters for the existence of Coulomb crystals are derived for both classical and quantum crystals. In the latter case, a critical mass ratio of the two charged components is found which is of the order of 80. Thus, holes in semiconductors with sufficiently flat valence bands are predicted to spontaneously order into a regular lattice. Such hole crystals are intimately related to ion Coulomb crystals in white dwarf and neutron stars as well as to ion crystals produced in the laboratory. A unified phase diagram of two-component Coulomb crystals is presented and is verified by first-principle computer simulations.
△ Less
Submitted 7 December, 2005; v1 submitted 10 July, 2005;
originally announced July 2005.
-
Phase Transition in Strongly Degenerate Hydrogen Plasma
Authors:
V. S. Filinov,
V. E. Fortov,
M. Bonitz,
P. R. Levashov
Abstract:
Direct fermionic path-integral Monte-Carlo simulations of strongly coupled hydrogen are presented. Our results show evidence for the hypothetical plasma phase transition. Its most remarkable manifestation is the appearance of metallic droplets which are predicted to be crucial for the electrical conductivity allowing to explain the rapid increase observed in recent shock compression measurments.
Direct fermionic path-integral Monte-Carlo simulations of strongly coupled hydrogen are presented. Our results show evidence for the hypothetical plasma phase transition. Its most remarkable manifestation is the appearance of metallic droplets which are predicted to be crucial for the electrical conductivity allowing to explain the rapid increase observed in recent shock compression measurments.
△ Less
Submitted 22 August, 2001;
originally announced August 2001.
-
Thermodynamic Properties of Correlated Strongly Degenerate Plasmas
Authors:
V. S. Filinov,
P. R. Levashov,
V. E. Fortov,
M. Bonitz
Abstract:
An efficient numerical approach to equilibrium properties of strongly coupled systems which include a subsystem of fermionic quantum particles and a subsystem of classical particles is presented. It uses an improved path integral representation of the many-particle density operator and allows to describe situations of strong coupling and strong degeneracy, where analytical theories fail. A novel…
▽ More
An efficient numerical approach to equilibrium properties of strongly coupled systems which include a subsystem of fermionic quantum particles and a subsystem of classical particles is presented. It uses an improved path integral representation of the many-particle density operator and allows to describe situations of strong coupling and strong degeneracy, where analytical theories fail. A novel numerical method is developed, which allows to treat degenerate systems with full account of the spin scatistics. Numerical results for thermodynamic properties such as internal energy, pressure and pair correlation functions are presented over a wide range of degeneracy parameter.
△ Less
Submitted 3 December, 1999;
originally announced December 1999.