SciPost Submission Page
The simplecubic structure of elemental Polonium and its relation to combined charge and orbital order in other elemental chalcogens
by Ana Silva, Jasper van Wezel
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Jasper van Wezel 
Submission information  

Preprint Link:  http://arxiv.org/abs/1712.09533v3 (pdf) 
Date accepted:  20180523 
Date submitted:  20180507 02:00 
Submitted by:  van Wezel, Jasper 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approach:  Theoretical 
Abstract
Polonium is the only element to crystallise into a simple cubic structure at ambient conditions. Moreover, at high temperatures it undergoes a structural phase transition into a less symmetric trigonal configuration. It has long been suspected that the strong spinorbit coupling in Polonium is involved in both peculiarities, but the precise mechanism by which it operates remains controversial. Here, we introduce a single microscopic model capable of capturing the atomic structure of all chalcogen crystals: Selenium, Tellurium, and Polonium. We show that the strong spinorbit coupling in Polonium suppresses the trigonal charge and orbital ordered state known to be the ground state configuration of Selenium and Tellurium, and allows the simple cubic state to prevail instead. We also confirm a recent suggestion based on ab initio calculations that a small increase in the lattice constant may effectively decrease the role of spinorbit coupling, leading to a reemergence of the trigonal orbital ordered state at high temperatures. We conclude that Polonium is a unique element, in which spins, orbitals, electronic charges, and lattice deformations all cooperate and collectively cause the emergence of the only elemental crystal structure with the simplest possible, cubic, lattice.
Author comments upon resubmission
We would like to sincerely thank the referee for their careful evaluation of our submission, and for their supportive comments and suggestions. We address the issues raised by the referee in the order in which they appear.

Strengths: we thank the referee for recognising the value of our work.

Weaknesses:

The referee points out that we make several assumptions. This is indeed true, and it leads to a very simple model for the family of elemental chalcogens. In fact, this should be seen as one of the main strengths of our approach. The experimentally observed structural phase diagram throughout a whole family of elements can be understood within a single, maximally simplified model. This is a clear indication that the model captures the essential physics.

The referee claims that we make no testable predictions. In fact there are at least two clear and testable predictions already in the paper. One is the presence of orbital order in all elemental chalcogens. At the moment, this type of order is hard to detect experimentally. Dedicated STM setups or nonlinear optical experiments however, could in principle detect the predicted orbital order. We point this out in the revised manuscript. The second prediction, is the phase diagram itself. As the referee already hints at in one of their later questions, looking at different isotopes of Po may be a way of tuning the value of spinorbit coupling, which would allow a direct exploration of the predicted phase diagram.
 Report: We thank the referee for their concise summary and positive remarks.

The referee asks about the nesting vector Q. The present manuscript builds on previous work done by the same authors, and published in Phys. Rev. B 97, 045151 (2018), where the possibility of a combined charge and orbitally ordered phase was first introduced. We therefore recognise that some aspects of the present work were only explained implicitly, by reference to our previous work. We have gone through the paper again, and added additional clarifications to address this issue. Concerning specifically the origin of the nesting vector Q: the model assumes as its starting point a simple cubic “parent” lattice, with a 2/3 filled band of porbitals. Considering only the dominant orbital overlap integrals, the simple cubic lattice consists of interwoven but independent onedimensional chains running in all three lattice directions. The resulting electronic structure then contains three pairs of parallel planar Fermi surfaces. This situation is extremely wellnested, and a Peierlstype charge density wave is expected to emerge. In fact, a single nesting vector Q, corresponding to a body diagonal of the cube of intersecting Fermi surfaces, connects any point on the Fermi surface to a point on a parallel Fermi surface sheet. A single nesting vector can therefore gap the entire Fermi surface, and the dominant instability will be towards the formation of charge density waves in each of the three orbital sectors, sharing the same propagation direction Q = (2pi/3a, 2pi/3a, 2pi/3a). We summarise this discussion in the revised manuscript.

The referee asks about the construction of the Hamiltonian. The tightbinding model was formulated in reciprocal space, and because the nesting vector is commensurate with the lattice, there is no problem with the size of the resulting Hamiltonian matrix. We add a remark explaining this in the revised manuscript.

The referee asks about the selfconsistency of the fixed phase relations between the three CDW in our model. The meanﬁeld equations in the present work were first solved selfconsistently without spinorbit coupling, keeping both the amplitudes and phases of the order parameter as free parameters. The reported phase relations for the selfconsistent solutions are the lowest energy solutions found this way. We then assume the phase relations not to change significantly as a function of spinorbit coupling. On the intuitive level, as long as the strength of spinorbit coupling is sufficiently weak compared to the effects of the other interaction terms, the chiral trigonal lattice structure is expected to survive. Since this phase results from a competition between Coulomb and electronphonon interactions, yielding combined charge and orbital order, it is consistent only with the fixed phase relations of eq (4). One may expect that at some critical value of the spinorbit coupling the orbital order breaks down. At that point, rather than slightly modifying the phase relations, the parent simple cubic lattice is expected to emerge as the only selfconsistent solution to the model. That no other stable phases exist was checked at several points in the phase diagram of figure 1. We extend the discussion of this point in the revised manuscript.

a) The referee asks about the linear fit to the temperature dependence of the phonon energy. We agree with the referee that a linear ﬁt can only hold true for a restricted window in temperature. It nevertheless suffices to gain a qualitative understanding of the observed structural phase diagram of Polonium, even if the results may not be quantitatively correct at temperatures far removed from the transition point.
b) The referee asks why we use an Einstein phonon mode in the model. The nesting vector Q in the chalcogens lies far from zero, so that even for acoustic modes the dispersion relation may be approximated to be constant (Einsteinlike) in the momentumspace region of interest. To account for the thermal evolution of the phonon mode, we assume the phonon energy to depend on the lattice constant, which in turn depends on temperature. Equation (5) then, is a firstorder expansion of the phonon energy with varying lattice constant. Combining this with the linear thermal expansion of the lattice discussed above, yields the temperature dependence of the phonon energy. We extend the discussion of this part of the analysis in the revised manuscript.

The referee asks for qualitative tests of the presented model. The first clear experimental prediction coming out of this manuscript is that the trigonal lattice structure of betaPolonium is in fact chiral, and orbital ordered. That is, it should be of precisely the same type as that observed in Selenium and Tellurium. For the latter two elements, the chirality of the lattice can be seen in Xray diffraction as well as optical activity measurements. The orbital order is harder to detect experimentally at the moment. Dedicated STM setups or nonlinear optical experiments, however, could in principle detect it. We point out these predictions in the revised manuscript. A second prediction could be based on the interesting suggestions of the referee. Looking at different isotopes of Po may be a way of tuning the value of spinorbit coupling, which would allow a direct exploration of the predicted phase diagram. Making quantitative predictions in that direction, however, is beyond the scope of the current work.

Requested changes: 1. We now include the mean field Hamiltonian in an appendix, as suggested.
 We rephrased the first line of the manuscript in the suggested manner.
 We break up this sentence as suggested by the referee.
 We adopt the suggested formulation.
List of changes
See reply above
Published as SciPost Phys. 4, 028 (2018)