4.7 Article

Two-phase simulation of the crystalline silicon melting line at pressures from-1 to 3 GPa

Journal

JOURNAL OF CHEMICAL PHYSICS
Volume 137, Issue 5, Pages -

Publisher

AMER INST PHYSICS
DOI: 10.1063/1.4739085

Keywords

-

Funding

  1. Russian Foundation for Basic Research [10-08-01266, 11-08-01225]
  2. Council of the President of the Russian Federation [NSh-7241.2012.2]
  3. Ministry of Education and Science of the Russian Federation [07.514.12.4002]

Ask authors/readers for more resources

Results of a numerical investigation of crystalline silicon melting line within the range of pressures from -1 to 3 GPa are presented. A two-phase molecular dynamics method is applied to obtain temperature, pressure, and densities of solid and liquid phases on the melting line. Using a special procedure we ensure the strict control of the two-phase equilibrium in the simulation cell. To describe the interaction between the atoms four classic potentials have been chosen: the Stillinger-Weber one and three modified variants of the Tersoff potential. For the Stillinger-Weber and Tersoff potentials in the modification by Kumagai-Izumi-Hara-Sakai a good coincidence with experimental data on crystalline Si melting temperature is obtained within the range of pressure from 0 to 3 GPa. Calculations of the solid and liquid phase densities on the silicon melting line for the Stillinger-Weber potential are also in close agreement with experiments. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4739085]

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Physics, Fluids & Plasmas

Thermodynamic properties of the finite-temperature electron gas by the fermionic path integral Monte Carlo method

V. S. Filinov, P. R. Levashov, A. S. Larkin

Summary: The quantum path integral Monte Carlo approach has been applied to the entropy difference calculations for strongly coupled degenerated uniform electron gas (UEG), providing data on pressure, internal energy, and entropy change, which may be beneficial for density functional theory. Published by AIP Publishing under an exclusive license.

PHYSICS OF PLASMAS (2021)

Article Physics, Multidisciplinary

Monte Carlo simulations of the electron short-range quantum ordering in Coulomb systems and the 'fermionic sign problem'

V Filinov, P. Levashov, A. Larkin

Summary: A new path integral representation of the density matrix in the canonical ensemble at finite temperatures has been developed to address the interference effects of the Coulomb and exchange interactions of electrons, reducing the 'fermionic sign problem' in Monte Carlo simulations. The results obtained for pair distribution functions in plasma and uniform electron gas demonstrate short-range quantum ordering of electrons associated with exchange-correlation excitons in the literature. Comparisons with available literature also show that the short-range ordering does not significantly contribute to integral thermodynamic characteristics.

JOURNAL OF PHYSICS A-MATHEMATICAL AND THEORETICAL (2022)

Article Physics, Fluids & Plasmas

Single-momentum path integral Monte Carlo simulations of uniform electron gas in warm dense matter regime

A. S. Larkin, V. S. Filinov, P. R. Levashov

Summary: This paper presents a single-momentum path integral Monte Carlo method adapted for simulations of the uniform electron gas, based on the combination of Wigner formalism and the path integral approach. By taking into account the exchange interaction between electrons and using the Gram determinants of the exchange matrix, the fermionic sign problem is significantly reduced and completely eliminated in coordinate-dependent variables, making it a valuable tool for studying thermodynamic properties.

PHYSICS OF PLASMAS (2021)

Article Physics, Applied

Ab initio inspection of thermophysical experiments for zirconium near melting

M. A. Paramonov, D. V. Minakov, V. B. Fokin, D. V. Knyazev, G. S. Demyanov, P. R. Levashov

Summary: This study presents quantum molecular dynamics calculations of thermophysical properties of solid and liquid zirconium near the melting point and provides an overview of available experimental data. By comparing the simulation results with experiments, possible reasons for discrepancies are discussed. The study reveals a significant volume change and low enthalpy of fusion on melting, and observes underestimation of electrical resistivity in the simulations while the slope of resistivity temperature dependencies agrees with experiments. The calculations predict almost constant normal spectral emissivity in liquid zirconium.

JOURNAL OF APPLIED PHYSICS (2022)

Article Physics, Multidisciplinary

Systematic derivation of angular-averaged Ewald potential

G. S. Demyanov, P. R. Levashov

Summary: 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's coefficients are found using Euler-Maclaurin and Poisson summation formulas, and are represented as finite series containing derivatives of Jacobi theta functions. The effectiveness of the potential is demonstrated through the calculation of the Madelung constant for various crystal lattices.

JOURNAL OF PHYSICS A-MATHEMATICAL AND THEORETICAL (2022)

Article Mathematics

Momentum Distribution Functions and Pair Correlation Functions of Unpolarized Uniform Electron Gas in Warm Dense Matter Regime

Alexander Larkin, Vladimir Filinov, Pavel Levashov

Summary: This paper investigates the uniform electron gas in a warm dense matter regime and reveals the presence of quantum tails and short-range order under non-ideal conditions.

MATHEMATICS (2022)

Article Physics, Fluids & Plasmas

A wide-range semiclassical self-consistent average atom model

A. S. Poliukhin, S. A. Dyachkov, A. A. Malyugin, P. R. Levashov

Summary: The discovery of material properties at extremes is crucial for the development of high energy density physics. Advanced experimental facilities, theories, and computations are needed for this purpose. The use of density functional theory (DFT) and path-integral Monte Carlo approaches have provided precise models for matter properties under extreme conditions. However, their practical usage is limited due to fundamental and computational constraints, and wide-range thermodynamic and transport models of plasma are required.

PHYSICS OF PLASMAS (2023)

Article Multidisciplinary Sciences

Ab Initio Calculations of Transport and Optical Properties of Dense Zr Plasma Near Melting

Vladimir Fokin, Dmitry Minakov, Pavel Levashov

Summary: The dynamic electrical conductivity of dense Zr plasma near melting is calculated using ab initio molecular dynamics and the Kubo-Greenwood formula. The influence of computational parameters and inner shell electrons on the results is thoroughly investigated. The convergence of the computations and comparison with experimental data are demonstrated.

SYMMETRY-BASEL (2023)

Article Computer Science, Interdisciplinary Applications

New features of the TFmix code: Thermodynamic properties of electrons in mixtures

O. P. Shemyakin, P. R. Levashov, P. A. Krasnova

Summary: In this paper, the authors review the finite-temperature Thomas-Fermi model for a mixture of elements and describe the modifications made for the new version (1.5) of the code. The model is used to calculate thermodynamic properties of electrons in a mixture of elements. The modifications include adding density as an input parameter, changing the calculation of thermal energy, adding more thermodynamic functions in the output, providing the ability to plot isotherms and isochores in the graphical user interface, and making minor corrections.

COMPUTER PHYSICS COMMUNICATIONS (2023)

Article Physics, Multidisciplinary

Density of states of a 2D system of soft-sphere fermions by path integral Monte Carlo simulations

V Filinov, P. Levashov, A. Larkin

Summary: The Wigner formulation of quantum mechanics is utilized to derive a new path integral representation of quantum density of states. A path integral Monte Carlo approach is developed for investigating the DOSs, internal energy, and spin-resolved radial distribution functions of a 2D system of strongly correlated soft-sphere fermions. The study examines and explains the peculiarities of DOSs and internal energy distributions depending on the softness of the potential and particle density, where the DOSs tend to a constant value at high densities, resembling an ideal 2D fermion system.

JOURNAL OF PHYSICS A-MATHEMATICAL AND THEORETICAL (2023)

Article Materials Science, Multidisciplinary

Ab initio calculation of hafnium and zirconium melting curves via the Lindemann criterion

D. V. Minakov, M. A. Paramonov, G. S. Demyanov, V. B. Fokin, P. R. Levashov

Summary: In this study, the melting curves of hafnium and zirconium are obtained using quantum molecular dynamics calculations. The results show that the melting curves of both materials are consistent with experimental data within the error range. Additionally, we provide the first estimate of the melting curve of hafnium under ultra-high pressures.

PHYSICAL REVIEW B (2022)

Article Physics, Applied

Quenching of Liquid Carbon on the Surface of a Diamond Substrate

V. S. Dozhdikov, A. Yu. Basharin, P. R. Levashov

Summary: Amorphous carbon is obtained by quenching liquid carbon for the first time. Laser melting of graphite flakes on a diamond substrate allows liquid carbon to fall on cold diamond for rapid cooling. Diamond acts as a unique quenching medium for extremely high-temperature molten carbon.

HIGH TEMPERATURE (2022)

Article Physics, Fluids & Plasmas

One-component plasma of a million particles via angular-averaged Ewald potential: A Monte Carlo study

G. S. Demyanov, P. R. Levashov

Summary: In this study, a correct expression for the one-component plasma energy is derived using the angular-averaged Ewald potential. The derived potential shows significant advantages in large-scale numerical simulations, offering a unified approach for determining the thermodynamic limit.

PHYSICAL REVIEW E (2022)

Article Physics, Fluids & Plasmas

Continuous Kubo-Greenwood formula: Theory and numerical implementation

G. S. Demyanov, D. Knyazev, P. R. Levashov

Summary: In this paper, the continuous Kubo-Greenwood formula for calculating the dynamic Onsager coefficients and real part of dynamic electrical conductivity is presented. It is formulated as an integral over the whole energy range, unlike the usual formula which involves summation over discrete transitions between electron energy levels. The continuous formula includes smooth functions that can be plotted and analyzed, allowing for the contribution analysis of various parts of the electron spectrum to the obtained properties. The formula has been implemented in the parallel code CUBOGRAM and used to study the influence of technical parameters on simulation results for liquid aluminum.

PHYSICAL REVIEW E (2022)

Article Physics, Fluids & Plasmas

Equilibrium properties of warm dense deuterium calculated by the wave packet molecular dynamics and density functional theory method

Yaroslav Lavrinenko, Pavel R. Levashov, Dmitry Minakov, Igor Morozov, Ilya A. Valuev

Summary: The joint simulation method of WPMD-DFT is utilized to study warm dense deuterium, providing a more accurate representation of quantum effects compared to classical molecular dynamics and WPMD simulations. The method allows for studying nonadiabatic dynamics of electrons and ions in equilibrium and nonequilibrium states, being more accurate and efficient at high densities. Results obtained by WPMD-DFT for shock-compressed deuterium are compared with experimental data and other simulation approaches to validate its accuracy for further application in studying nonequilibrium states and relaxation processes.

PHYSICAL REVIEW E (2021)

No Data Available