Fate of density functional theory in high-pressure solid hydrogen

This paper investigates some of the successes and failures of density functional theory in the study of high-pressure solid hydrogen at low temperature. We calculate the phase diagram, metallization pressure, phonon spectrum, and proton zero-point energy using three popular exchange-correlation functionals: the local density approximation (LDA), the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation, and the semi-local Becke-Lee-Yang-Parr (BLYP) functional. We focus on the solid molecular P$6_3$/m, C2/c, Cmca-12, and Cmca structures in the pressure range from $100<P<500$ GPa over which phases I, II and III are observed experimentally. At the static level of theory, in which proton zero-point energy is ignored, the LDA, PBE and BLYP functionals give very different structural transition and metallization pressures, with the BLYP phase diagram in better agreement with experiment. Nevertheless, all three functionals provide qualitatively the same information about the band gaps of the four structures and the phase transitions between them. Going beyond the static level, we find that the frequencies of the vibron modes observed above 3000 cm$^{-1}$ depend strongly on the choice of exchange-correlation functional, although the low-frequency part of the phonon spectrum is little affected. The largest and smallest values of the proton zero-point energy, obtained using the BLYP and LDA functionals, respectively, differ by more than 10 meV/proton. Including the proton zero-point energy calculated from the phonon spectrum within the harmonic approximation improves the agreement of the BLYP and PBE phase diagrams with experiment. Taken as a whole, our results demonstrate the inadequacy of mean-field-like density functional calculations of solid molecular hydrogen in phases I, II and III and emphasize the need for more sophisticated methods.Comment: Accepted for publicatio

Beyond "the Relationship between the Individual and Society": broadening and deepening relational thinking in group analysis

The question of ‘the relationship between the individual and society’ has troubled group analysis since its inception. This paper offers a reading of Foulkes that highlights the emergent, yet evanescent, psychosocial ontology in his writings, and argues for the development of a truly psychosocial group analysis, which moves beyond the individual/society dualism. It argues for a shift towards a language of relationality, and proposes new theoretical resources for such a move from relational sociology, relational psychoanalysis and the ‘matrixial thinking’ of Bracha Ettinger which would broaden and deepen group analytic understandings of relationality

Quantum Monte Carlo Analysis of Exchange and Correlation in the Strongly Inhomogeneous Electron Gas

We use variational quantum Monte Carlo to calculate the density-functional exchange-correlation hole n_{xc}, the exchange-correlation energy density e_{xc}, and the total exchange-correlation energy E_{xc}, of several electron gas systems in which strong density inhomogeneities are induced by a cosine-wave potential. We compare our results with the local density approximation and the generalized gradient approximation. It is found that the nonlocal contributions to e_{xc} contain an energetically significant component, the magnitude, shape, and sign of which are controlled by the Laplacian of the electron density.Comment: 4 pages, 3 figure

Quantum Monte Carlo Study of High Pressure Solid Molecular Hydrogen

We use the diffusion quantum Monte Carlo (DMC) method to calculate the ground state phase diagram of solid molecular hydrogen and examine the stability of the most important insulating phases relative to metallic crystalline molecular hydrogen. We develop a new method to account for finite-size errors by combining the use of twist-averaged boundary conditions with corrections obtained using the Kwee-Zhang-Krakauer (KZK) functional in density functional theory. To study band-gap closure and find the metallization pressure, we perform accurate quasi-particle many-body calculations using the $GW$ method. In the static approximation, our DMC simulations indicate a transition from the insulating Cmca-12 structure to the metallic Cmca structure at around 375 GPa. The $GW$ band gap of Cmca-12 closes at roughly the same pressure. In the dynamic DMC phase diagram, which includes the effects of zero-point energy, the Cmca-12 structure remains stable up to 430 GPa, well above the pressure at which the $GW$ band gap closes. Our results predict that the semimetallic state observed experimentally at around 360 GPa [Phys. Rev. Lett. {\bf 108}, 146402 (2012)] may correspond to the Cmca-12 structure near the pressure at which the band gap closes. The dynamic DMC phase diagram indicates that the hexagonal close packed $P6_3/m$ structure, which has the largest band gap of the insulating structures considered, is stable up to 220 GPa. This is consistent with recent X-ray data taken at pressures up to 183 GPa [Phys. Rev. B {\bf 82}, 060101(R) (2010)], which also reported a hexagonal close packed arrangement of hydrogen molecules

Dissociation of high-pressure solid molecular hydrogen: Quantum Monte Carlo and anharmonic vibrational study

A theoretical study is reported of the molecular-to-atomic transition in solid hydrogen at high pressure. We use the diffusion quantum Monte Carlo method to calculate the static lattice energies of the competing phases and a density-functional-theory-based vibrational self-consistent field method to calculate anharmonic vibrational properties. We find a small but significant contribution to the vibrational energy from anharmonicity. A transition from the molecular Cmca-12 direct to the atomic I4_1/amd phase is found at 374 GPa. The vibrational contribution lowers the transition pressure by 91 GPa. The dissociation pressure is not very sensitive to the isotopic composition. Our results suggest that quantum melting occurs at finite temperature.Comment: Accepted for publication by Phys. Rev. Let

Hydrodynamic modelling of accretion flows

In the proceedings of this, and of several recent close binary conferences, there have been several contributions describing smoothed particle hydrodynamics simulations of accretion disks. It is apposite therefore to review the numerical scheme itself with emphasis on its advantages for disk modelling, and the methods used for modelling viscous processes.Comment: 3 pages, to appear in proceedings of IAU Colloquium 194: Compact binaries in the galaxy and beyon

The effect of quantization on the FCIQMC sign problem

The sign problem in Full Configuration Interaction Quantum Monte Carlo (FCIQMC) without annihilation can be understood as an instability of the psi-particle population to the ground state of the matrix obtained by making all off-diagonal elements of the Hamiltonian negative. Such a matrix, and hence the sign problem, is basis dependent. In this paper we discuss the properties of a physically important basis choice: first versus second quantization. For a given choice of single-particle orbitals, we identify the conditions under which the fermion sign problem in the second quantized basis of antisymmetric Slater determinants is identical to the sign problem in the first quantized basis of unsymmetrized Hartree products. We also show that, when the two differ, the fermion sign problem is always less severe in the second quantized basis. This supports the idea that FCIQMC, even in the absence of annihilation, improves the sign problem relative to first quantized methods. Finally, we point out some theoretically interesting classes of Hamiltonians where first and second quantized sign problems differ, and others where they do not.Comment: 4 pages w/ 2 page appendix, 2 figures, 1 tabl