Abstract:Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at "uninteresting" state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble ( NPTΔμ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles. https://doi.org/10.1103/PhysRevE.103.043204 ©2021 American Physical Society

Virialized equation of state for warm and dense stellar plasmas in proto-neutron stars and Supernova matter

Virialized equation of state for warm and dense stellar plasmas in proto-neutron stars and supernova matter

Equation of state of dense plasmas by ab initio simulations: Bridging the gap between quantum molecular dynamics and orbital-free molecular dynamics at high temperature

\emph{Ab initio} Quantum Monte Carlo simulation of the warm dense electron gas

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Finite-temperature expansion of the dense-matter equation of state

Modeling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

Generalised description of Neutron Star matter with nucleonic Relativistic Density Functional

Anomalous thermal relaxation in warm ion plasmas

Equation of state of neutron star matter and its warm extension with an interacting hadron resonance gas

Constraints on the dense matter equation of state from young and cold isolated neutron stars

Plasma screening of nuclear fusion reactions in liquid layers of compact degenerate stars: a first-principle study

Statistical description of complex nuclear phases in supernovae and proto-neutron stars

Far-from-equilibrium bulk-viscous transport coefficients in neutron star mergers

Thermal x-ray studies of neutron stars and the equation of state

Equation of State of Hot Neutron Star Matter using Finite Range Simple Effective Interaction

Impact of Hot Inner Crust on Compact Stars at Finite Temperature

The proto-neutron star inner crust in the liquid phase

Hot and highly magnetized neutron star matter properties with Skyrme interactions

Virial coefficients of the uniform electron gas from path-integral Monte Carlo simulations