We predict a quantum phase transition in fcc Ca under hydrostatic pressure. Using density functional theory, we find at pressures below 80 kbar, the topology of the electron charge density is characterized by nearest neighbor atoms connected through bifurcated bond paths and deep minima in the octahedral holes. At pressures above 80 kbar, the atoms bond through non-nuclear maxima that form in the octahedral holes. This topological change in the charge density softens the C elastic modulus of fcc Ca, while C44 remains unchanged. We propose an order parameter based on applying Morse theory to the charge density, and we show that near the critical point it follows the expected mean-field scaling law with reduced pressure.
PACS numbers:The theory and characterization of solid-solid isostructural phase transitions has been an area of intense experimental and theoretical research since they were described in 1949 [1]. What sets these transitions apart from others is that while crystalline properties change-sometimes discontinuously-crystallographic structure does not, making these transitions purely electronic in nature. Both first order and continuous transitions are known to exist. A well known example of the former is the γ → α transition in fcc Ce [1], while a magneto-volume instability in fcc Fe is an example of the later [2].In the case of Fe, the traditional local order parameter, µ, gave conflicting results computationally, with the transition appearing as first [3] or second [2] order. It was not until the transition was interpreted as a topological change of the elctronic charge density, ρ( r), that the underlying nature of the phase change became apparent [4]. While this approach is promising, no general correlations have been drawn between isostructural phase transitions and ρ( r). We use this approach to predict an isostructural phase transition in calcium under hydrostatic pressure. We then apply Morse theory to the charge density and propose an order parameter for the transition. We show that this order parameter scales in the appropriate way with reduced pressure near the quantum critical point.One of the advantages of casting theories of isostructural phase transition in terms of ρ( r) is that it is a quantum observable. Though ρ( r) is often calculated by way of first principles methods (e.g. density functional theory (DFT)), the total charge density can also be measured via X-ray diffraction techniques, and the spin-polarized charge density can be determined using spin-polarized neutron diffraction. As it is also known that all ground state properties depend on the charge density [5], it seems appropriate to seek relationships between the structure * Electronic address: trjones@mines.edu of ρ( r), changes to that structure, and corresponding changes to properties. We note that previous work has proposed that DFT could be used as a general framework for analyzing quantum phase transitions in electronic systems [6].The electron charge density is a 3 dimensional scalar field. Morse theory tells us the ...