**Principal Investigator:**

Carlo Pierleoni

**Affiliation:**

Universita' dell'Aquila (Italy)

**Local Project ID:**

QMCSim

**HPC Platform used:**

Hermit of HLRS

**Date published:**

**Quantum Mechanics is the fundamental theory on which our present understanding of the physical world is based. During the last century the “Founding Fathers” have developed the theory and established its fundamental character, and now we are facing the problem of using the laws to predict the behavior of systems of many particles in a wide range of physical conditions, from the ultra-cold states of matter to the extreme conditions of temperature and pressure found in the interiors of planets. This is a phenomenal challenge that can only be attacked by numerical methods.**

The current description for extended systems of nuclei and electrons is based on a self-consistent mean field theory which is successful for many applications but has some well known limitations. In the specific case of hydrogen at the interesting physical conditions, related to planetary science, to high pressure experiments and to energy applications, this theory is not accurate enough and new fundamental developments are required.

The main goals of a project under the leadership of Carlo Pierleoni of the University of L’Aquila, Italy, were: 1) to elucidate important aspects of the hydrogen phase diagram related to the pressure-induced molecular dissociation and metallization. In particular a) the occurrence of two coexisting liquid phases (one molecular and insulating, the other monatomic and metallic) and the related reentrant melting line of the molecular solid; b) the possible existence of a liquid ground state of metallic hydrogen, a new state of matter. 2) to improve the treatment of electronic correlation by developing algorithms based on Quantum Monte Carlo (QMC) methods. QMC is much more accurate than the mean-field theory but requires at least an order of magnitude more computer resources. However it is highly parallelizable and scales extremely well on many thousands of core, therefore is ideal for Tier-0 machines like Hermit.

The main results achieved are: 1) a new prediction of the liquid-liquid transition line with Coupled Electron-Ion Monte Carlo (a QMC-based method) including a proper treatment of nuclear quantum effects. These results are to date the most accurate predictions of this phase transition. Our results will be a reference for future experiments and also a benchmark for better approximations within the mean-field framework. 2) Benchmark studies of the mean-field method with several approximation which has been the subject of three publications. 3) Evidence for a stable crystalline low temperature phase of metallic hydrogen which rules out the ground state liquid.

The study will be very useful for the large community of scientists that uses the approximate theory to better understand its limitations. More in general the new QMC algorithms will constitute the method to overcome the mean-field treatment of electronic correlations in the future.

The project was made possible through the Partnership for Advanced Computing in Europe, PRACE, which allocated 24 million core hours of computing time on GCS supercomputer Hermit of HLRS Stuttgart to the project.

The team was international with a nearly equal number of participants from Italy and from the USA, and one additional member from France. The same team also benefitted from other computer time allocations within the INCITE-USA program.

Carlo Pierleoni

Dipartimento di Scienze Fisiche e Chimiche

Universita' dell'Aquila

Via Vetoio, 10 I-67100 L'AQUILA/Italy

e-mail: carlo.pierleoni@aquila.infn.it

[Datchi2000] Datchi, F., P. Loubeyre, and R. LeToullec, 2000, Phys. Rev. B 61, 6535.

[Gregorianz2013] Gregoryanz, E., A.F. Goncharov, K. Matsuishi, H.-k. Mao, and R. J. Hemley, 2003, Phys. Rev. Lett. 90, 175701.

[Deemyad2008] Deemyad, S., and I. F. Silvera, 2008, Phys. Rev. Lett. 100, 155701.

[Eremetes2009] Eremets, M., and I. Trojan, 2009, JETP Lett. 89, 174.

[Subramanian2011] Subramanian, N., A. Goncharov, V. Struzhkin, M. Somayazulu, and R. Hemley, 2011, Proc. Natl. Acad. Sci. U.S.A. 108, 6014.

[Morales2010] Morales, M. A., C. Pierleoni, E. Schwegler, and D. M.Ceperley, 2010, Proc. Natl. Acad. Sci. U.S.A. 107, 12799.

[Bonev2004] Bonev, S., E. Schwegler, T. Ogitsu, and G. Galli, 2004, Nature 431, 669.

[Mazin1997] Mazin, I. I., and R. E. Cohen, 1995, Phys. Rev. B 52, R8597.

[Eremetes2011] Eremets, M. I., and I. A. Troyan, 2011, Nature Mater. 10, 927.

[Attaccalite2008] Attaccalite, C., and S. Sorella, 2008, Phys. Rev. Lett. 100, 114501.

[Howie2012] Howie, R. T., C. L. Guillaume, T. Scheler, A. F. Goncharov, and E.Gregoryanz, 2012, Phys. Rev. Lett. 108, 125501.

[McMahon2014] J.M. McMahon, M.A. Morales, C. Pierleoni and D.M. Ceperley, “Monatomic hydrogen: a ground-state fluid or a metallic solid?”, submitted.

[Chen2013] J. Chen, X.-Z. Li, Q. Zhang, M. I. J. Probert, C. J. Pickard, R. J. Needs, A. Michaelides, and E. Wang, Nature Communications 4, 2064 (2013).

[Morales2013] M. A. Morales, J. M. McMahon, C. Pierleoni, and D. M. Ceperley, Phys. Rev. Lett. 110, 065702 (2013).