Source author record

P. R. Levashov

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

19 published item(s)

preprint2022arXiv

Derivation of the Kelbg potential/functional

The density matrix for a system of particles interacting via the Coulomb potential is obtained in the high--temperature limit following almost entirely the original work by Kelbg. For this purpose the Blöch equation is solved in the first order of perturbation theory. We tried to explain all the transformations in the derivation in order to simplify the understanding of this non-trivial theory. The solution of Kelbg is widely used in path integral simulations of Coulomb systems.

preprint2022arXiv

Momentum distribution functions and pair correlation functions of unpolarized uniform electron gas in warm dense matter regime

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.

preprint2022arXiv

One--Component Plasma of a Million Particles via angular--averaged Ewald potential: A Monte Carlo study

In this work, we derive a correct expression for the one--component plasma (OCP) energy via the angular--averaged Ewald potential (AAEP). Unlike E.~Yakub and C.~Ronchi (J. Low Temp. Phys. 139, 633 (2005)), who had tried to obtain the same energy expression from a two--component plasma model, we used the original Ewald potential for an OCP. A constant in the AAEP was determined using the cluster expansion in the limit of weak coupling. The potential has a simple form suitable for effective numerical simulations. To demonstrate the advantages of the AAEP, we performed a number of Monte--Carlo simulations for an OCP with up to a million particles in a wide range of the coupling parameter. Our computations turned out at least two orders of magnitude more effective than those with a traditional Ewald potential. A unified approach is offered for the determination of the thermodynamic limit in the whole investigated range. Our results are in good agreement with both theoretical data for a weakly coupled OCP and previous numerical simulations. We hope that the AAEP will be useful in path integral Monte Carlo simulations of the uniform electron gas.

preprint2022arXiv

Systematic derivation of angular--averaged Ewald potential

In this work we provide a step by step derivation of an angular--averaged Ewald potential suitable for numerical simulations of disordered Coulomb systems. The potential was first introduced by E.\,Yakub and C.\,Ronchi without a clear derivation. Two methods are used to find the coefficients of the series expansion of the potential: based on the Euler--Maclaurin and Poisson summation formulas. The expressions for each coefficient is represented as a finite series containing derivatives of Jacobi theta functions. We also demonstrate the formal equivalence of the Poisson and Euler--Maclaurin summation formulas in the three-dimensional case. The effectiveness of the angular--averaged Ewald potential is shown by the example of calculating the Madelung constant for a number of crystal lattices.

preprint2021arXiv

Ab initio inspection of thermophysical experiments for zirconium near melting

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.

preprint2020arXiv

Effect of the spin-orbit interaction on thermodynamic properties of liquid uranium

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.

preprint2016arXiv

Ab initio calculation of thermodynamic, transport, and optical properties of CH$_2$ plastics

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.

preprint2014arXiv

Transport and optical properties of warm dense aluminum in the two-temperature regime: Ab initio calculation and semiempirical approximation

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.

preprint2013arXiv

\textit{Ab initio} calculation of transport and optical properties of aluminum: influence of simulation parameters

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.

preprint2013arXiv

Ab initio simulation of complex dielectric function for dense aluminum plasma

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.

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.

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.

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.

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.

preprint2010arXiv

Suppression of ablation in femtosecond double pulse experiments

We report the physical reasons of a curious decrease in the crater depth observed for long delays in experiments with femtosecond double pulses. Detailed hydrodynamic modeling demonstrates that the ablation mechanism is dumped when the delay between the pulses exceeds the electron-ion relaxation time. In this case, the interaction of the second laser pulse with the expanding target material leads to the formation of the second shock wave suppressing the rarefaction wave created by the first pulse. The evidence of this effect follows from the pressure and density profiles obtained at different delays after the first laser pulse.

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.

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.

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.