Properties of B4 C in the shocked state for pressures up to 1.5 TPa

Density functional theory calculations using the quasiharmonic approximation have been used to calculate the solid Hugoniot of two polytypes of boron carbide up to 100 GPa. Under the assumption that segregation into the elemental phases occurs around the pressure that the B11Cp(CBC) polytype becomes thermodynamically unstable with respect to boron and carbon, two discontinuities in the Hugoniot, one at 50 GPa and one at 90-100 GPa, are predicted. The former is a result of phase segregation, and the latter a phase transition within boron. First-principles molecular dynamics (FPMD) simulations were employed to calculate the liquid Hugoniot of B4C up to 1.5 TPa, and the results are compared to recent experiments carried out at the Omega Laser Facility up to 700 GPa [D. E. Fratanduono, Phys. Rev. B 94, 184107 (2016)1098-012110.1103/PhysRevB.94.184107]. A generally good agreement between theory and experiment was observed. Analysis of the FPMD simulations provides evidence for an amorphous, but covalently bound, fluid below 438 GPa, and an atomistic fluid at higher pressures and temperatures.