First-principles calculations of the epsilon phase of solid oxygen

The crystal, electronic and magnetic structures of solid oxygen in the epsilon phase have been investigated using the strongly constrained appropriately normed (SCAN) + rVV10 method and the generalized gradient approximation (GGA) + vdW-D + U method. The spin-polarized SCAN + rVV10 method with an 8-atom primitive unit cell provides lattice parameters consistent with the experimental results over the entire pressure range, including the epsilon-zeta structural phase transition at high pressure, but does not provide accurate values of the intermolecular distances d(1) and d(2) at low pressure. The agreement between the intermolecular distances and the experimental values is greatly improved when a 16-atom conventional unit cell is used. Therefore, the SCAN + rVV10 method with a 16-atom unit cell can be considered the most suitable model for the epsilon phase of solid oxygen. The spin-polarized SCAN + rVV10 model predicts a magnetic phase at low pressure. Since the lattice parameters of the predicted magnetic structure are consistent with the experimental lattice parameters measured at room temperature, our results may suggest that the epsilon phase is magnetic even at room temperature. The GGA + vdW-D + U (with an ad hoc value of U(eff) = 2 eV at low pressure instead of the first-principles value of U(lr)(eff) ~ 9 eV) and hybrid functional methods provide similar results to the SCAN + rVV10 method; however, they do not provide reasonable values for the intermolecular distances.