Crystallization of carbon-oxygen mixtures in white dwarf stars

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

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.

Actinide Crystallization and Fission Reactions in Cooling White Dwarf Stars

The first solids that form as a cooling white dwarf (WD) starts to crystallize are expected to be greatly enriched in actinides. This is because the melting points of WD matter scale as Z^{5/3} and actinides have the largest charge Z. We estimate that the solids may be so enriched in actinides that they could support a fission chain reaction. This reaction could ignite carbon burning and lead to the explosion of an isolated WD in a thermonuclear supernova (SN Ia). Our mechanism could potentially explain SN Ia with sub-Chandrasekhar ejecta masses and short delay times.

Chemical element transport in stellar evolution models

Stellar evolution computations provide the foundation of several methods applied to study the evolutionary properties of stars and stellar populations, both Galactic and extragalactic. The accuracy of the results obtained with these techniques is linked to the accuracy of the stellar models, and in this context the correct treatment of the transport of chemical elements is crucial. Unfortunately, in many respects calculations of the evolution of the chemical abundance profiles in stars are still affected by sometimes sizable uncertainties. Here, we review the various mechanisms of element transport included in the current generation of stellar evolution calculations, how they are implemented, the free parameters and uncertainties involved, the impact on the models and the observational constraints.

Direct molecular dynamics simulation of liquid-solid phase equilibria for a three-component plasma

The neutron-rich isotope ²²Ne may be a significant impurity in carbon and oxygen white dwarfs and could impact how the stars freeze. We perform molecular dynamics simulations to determine the influence of ²²Ne in carbon-oxygen-neon systems on liquid-solid phase equilibria. Both liquid and solid phases are present simultaneously in our simulation volumes. We identify liquid, solid, and interface regions in our simulations using a bond angle metric. In general we find good agreement for the composition of liquid and solid phases between our MD simulations and the semianalytic model of Medin and Cumming. The trace presence of a third component, neon, does not appear to strongly impact the chemical separation found previously for two-component carbon and oxygen systems. This suggests that small amounts of ²²Ne may not qualitatively change how the material in white dwarf stars freezes. However, we do find systematically lower melting temperatures (higher Γ) in our MD simulations compared to the semianalytic model. This difference seems to grow with impurity parameter Q_{imp} and suggests a problem with simple corrections to the linear mixing rule for the free energy of multicomponent solid mixtures that is used in the semianalytic model.

Direct molecular dynamics simulation of liquid-solid phase equilibria for two-component plasmas

We determine the liquid-solid phase diagram for carbon-oxygen and oxygen-selenium plasma mixtures using two-phase molecular dynamics simulations. We identify liquid, solid, and interface regions using a bond angle metric. To study finite-size effects, we perform 27,648- and 55,296-ion simulations. To help monitor nonequilibrium effects, we calculate diffusion constants D(i). For the carbon-oxygen system we find that D(O) for oxygen ions in the solid is much smaller than D(C) for carbon ions and that both diffusion constants are 80 or more times smaller than diffusion constants in the liquid phase. There is excellent agreement between our carbon-oxygen phase diagram and that predicted by Medin and Cumming. This suggests that errors from finite-size and nonequilibrium effects are small and that the carbon-oxygen phase diagram is now accurately known. The oxygen-selenium system is a simple two-component model for more complex rapid proton capture nucleosynthesis ash compositions for an accreting neutron star. Diffusion of oxygen, in a predominantly selenium crystal, is remarkably fast, comparable to diffusion in the liquid phase. We find a somewhat lower melting temperature for the oxygen-selenium system than that predicted by Medin and Cumming. This is probably because of electron screening effects.

Diffusion in Coulomb crystals

Diffusion in Coulomb crystals can be important for the structure of neutron star crusts. We determine diffusion constants D from molecular dynamics simulations. We find that D for Coulomb crystals with relatively soft-core 1/r interactions may be larger than D for Lennard-Jones or other solids with harder-core interactions. Diffusion, for simulations of nearly perfect body-centered-cubic lattices, involves the exchange of ions in ringlike configurations. Here ions "hop" in unison without the formation of long lived vacancies. Diffusion, for imperfect crystals, involves the motion of defects. Finally, we find that diffusion, for an amorphous system rapidly quenched from Coulomb parameter Γ=175 to Coulomb parameters up to Γ=1750, is fast enough that the system starts to crystalize during long simulation runs. These results strongly suggest that Coulomb solids in cold white dwarf stars, and the crust of neutron stars, will be crystalline and not amorphous.

Diffusion of neon in white dwarf stars

Sedimentation of the neutron rich isotope 22Ne may be an important source of gravitational energy during the cooling of white dwarf stars. This depends on the diffusion constant for 22Ne in strongly coupled plasma mixtures. We calculate self-diffusion constants D(i) from molecular dynamics simulations of carbon, oxygen, and neon mixtures. We find that D(i) in a mixture does not differ greatly from earlier one component plasma results. For strong coupling (coulomb parameter Γ> few), D(i) has a modest dependence on the charge Z(i) of the ion species, D(i)∝Z(i)(-2/3). However, D(i) depends more strongly on Z(i) for weak coupling (smaller Γ). We conclude that the self-diffusion constant D(Ne) for 22Ne in carbon, oxygen, and neon plasma mixtures is accurately known so that uncertainties in D(Ne) should be unimportant for simulations of white dwarf cooling.