Universal Diffusion in Coulomb Crystals

M. E. Caplan,D. Yaacoub

2024-09-12

Abstract:Diffusion coefficients for crystallized Coulomb plasmas are essential microphysics input for modeling white dwarf cores and neutron star crusts but are poorly understood. In this work we present a model for diffusion in Coulomb crystals. We show that melting and diffusion follow the same universal scaling such that diffusion is independent of screening. Our simulations show, contrary to prevailing wisdom, that the formation of vacancies is not suppressed by the large pressure. Rather, vacancy formation and hole diffusion is the dominant mode of self diffusion in Coulomb crystals.

Plasma Physics,High Energy Astrophysical Phenomena

What problem does this paper attempt to address?

The paper attempts to address the issue of diffusion coefficients in Coulomb crystals. Specifically, the authors focus on the diffusion process in Coulomb crystals under extremely high-pressure conditions, such as those found in the cores of white dwarfs and the crusts of neutron stars. Although diffusion coefficients under these conditions are crucial for astrophysical models, the current understanding in this field remains very limited.

### Main Issues:
1. **Universality of Diffusion Coefficients**: The authors aim to establish a diffusion model for Coulomb crystals with different screening lengths, demonstrating the universality of the relationship between diffusion coefficients and melting temperature.
2. **Mechanism of Defect Formation**: The traditional view holds that vacancy formation is suppressed under high pressure, so diffusion mainly occurs through exchanges of neighboring atoms. However, the authors' simulations reveal that the spontaneous formation of vacancies and their subsequent migration within the lattice is the primary mode of diffusion.
3. **Many-body Correlation in Diffusion Process**: The authors emphasize that diffusion in solids is a highly correlated many-body process, significantly different from the diffusion mechanism in liquids.

### Background and Motivation:
- **Strongly Coupled Coulomb Plasma**: This type of plasma exists under various extreme conditions, including dusty plasmas, inertial confinement experiments, the cores of white dwarfs, and the crusts of neutron stars.
- **Microscopic Transport Properties**: These properties are crucial for high-fidelity computational modeling, but due to the extreme temperatures and densities, direct experimental measurements are challenging and often rely on numerical simulations.
- **Importance of Diffusion Coefficients**: Diffusion coefficients are important for modeling direct transport processes of matter (such as chemical separation in white dwarfs and neutron stars) and magnetar elastic phenomena (such as crustal fractures and creep).

### Methods and Results:
- **Molecular Dynamics Simulations**: The authors used large-scale N-body simulations to calculate the diffusion coefficients of Coulomb crystals and demonstrated the detailed mechanisms of the diffusion process.
- **Universal Scaling Law**: The authors proposed a dimensionless diffusion coefficient \( D^* \) and validated its universal relationship with the melting temperature through simulations.
- **Vacancy Diffusion**: Simulation results show that the formation and migration of vacancies are the primary modes of diffusion, rather than the traditional exchange of neighboring atoms.

### Conclusions:
- **Universality**: Regardless of the screening length, the temperature relative to the melting temperature is the only key parameter determining the diffusion coefficient.
- **Vacancy Dominance**: Even under extremely high pressure, the formation and migration of vacancies remain the primary diffusion mechanism.
- **Future Work**: Further research is needed on the equilibrium density of vacancies and the activation energies for their formation and migration to improve the model's accuracy.

Through these studies, the authors provide an important theoretical foundation and experimental basis for understanding and modeling the diffusion process in Coulomb crystals in astrophysical contexts.