Student registration form


Education


Supervisor registration form


Monte Carlo vs. Two Component DFT

In the course of working with DFT-based software packages, we were forced to circumvent numerous problems arising from the nature of the spatial-wave mismatch with real physical models. For this reason, it is interesting for us to approach the study of the structural features of materials in a new way, as well as the defects that arise in them as a result of different types of interventions. The study of matter by means of annihilation spectroscopy is a relatively well-modeled process, but the exact comparison with the experiment determines a significant deviation from the current models we use. For this reason, we would like to try to model correlation effects of positron-electron wave functions including long-range correlations by using a quantum Monte Carlo approach.

Tasks

As mentioned above, the task of accurately determining the positron lifetime in a perfect lattice was performed by colleagues K. A. Simula et al., which opens up an opportunity for us to test our capabilities and continue work towards determining the momentum density by the Monte Carlo method. Thus, we formulate the main task of this project: "Theoretical justification of the electron momentum density by calculating positron annihilation by the Monte Carlo method".

Preliminary schedule by topics/tasks

1 To understand and prepare a brief review of the Two-Component DFT.
2 To review the work of K. A. Simula et al. and try to establish contact with their working group in order to understand what steps they have taken.
3 To attempt to add additional code to the current version of the program - an upgrade.

Required skills

Fortran programming, Basic concepts of what TCDFT is, Knowledge of Monte Carlo modeling.

Acquired skills and experience

Acquisition of skills for working with algorithms serving modeling of positron annihilation in condensed matter.

Recommended literature

[1] N. D. Drummond, P. L´opez R´ıos, R. J. Needs, and C. J. Pickard, Quantum Monte Carlo study of a positron in an electron gas, Phys. Rev. Lett. 107, 207402 (2011),
[2] B. Barbiellini and J. Kuriplach, Proposed parameter-free model for interpreting the measured positron annihila tion spectra of materials using a generalized gradient ap proximation, Phys. Rev. Lett. 114, 147401 (2015),
[3] K. A. Simula, J. E. Muff, I. Makkonen, and N. D. Drummond, Quantum Monte Carlo study of positron lifetimes in solids, Phys Rev Lett. 2022 Oct 14;129(16):166403,
[4] R. Q. Hood, P.R.C. Kent, R.J. Needs, and P.R. Briddon, Quantum monte carlo study of the optical and diffusive properties of the vacancy defect in diamond, Phys. Rev. Lett. 91, 076403 (2003),
[5] D. M. Ceperley and B. J. Alder, Ground state of the electron gas by a stochastic method, Phys. Rev. Lett. 45, 566 (1980),
[6] W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Ra jagopal, Quantum Monte Carlo simulations of solids, Rev. Mod. Phys. 73, 33 (2001),
[7] R. J. Needs, M. D. Towler, N. D. Drummond, and P. L´opez R´ıos, Continuum variational and diffusion quan tumMonteCarlocalculations, J. Phys.: Condens. Matter 22, 023201 (2009),
[8] R. J. Needs, M. D. Towler, N. D. Drummond, P. L´opez R´ıos, and J. R. Trail, Variational and diffusion quantum Monte Carlo calculations with the CASINO code, J. Chem. Phys. 152, 154106 (2020),
[9] J. B. Anderson, A random-walk simulation of the Schr¨odinger equation: H+ 3, J. Chem. Phys. 63, 1499 (1975),
[10] R. Jastrow, Many-body problem with strong forces, Phys. Rev. 98, 1479 (1955),
[11] P. L´opez R´ıos, A. Ma, N. D. Drummond, M. D. Towler, and R. J. Needs, Inhomogeneous backflow transforma tions in quantum Monte Carlo calculations, Phys. Rev. E 74, 066701 (2006),
[12] W. Kohn and L. J. Sham, Self-consistent equations in cluding exchange and correlation effects, Phys. Rev. 140, A1133 (1965),
[13] P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ cioni, I. Dabo, et al., Quantum ESPRESSO: a modular and open-source software project for quantum simula tions of materials, J. Phys.: Condens. Matter 21, 395502 (2009)
[14] T. Torsti, T. Eirola, J. Enkovaara, T. Hakala, P. Havu, V. Havu, T. H¨oyn¨al¨anmaa, J. Ignatius, M. Lyly, I. Makkonen, et al., Three real-space discretization tech niques in electronic structure calculations, Phys. Status Solidi (b) 243, 1016 (2006),
[15] I. Makkonen, M. Hakala, and M. J. Puska, Modeling the momentum distributions of annihilating electron positron pairs in solids, Phys. Rev. B 73, 035103 (2006),
[16] J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996),
[17] E. Boroński and R. M. Nieminen, Electron-positron density-functional theory, Phys. Rev. B 34, 3820 (1986),
[18] D. Alf`e and M. J. Gillan, Efficient localized basis set for quantum Monte Carlo calculations on condensed matter, Phys. Rev. B 70, 161101(R) (2004),
[19] N. D. Drummond, M. D. Towler, and R. J. Needs, Jas trow correlation factor for atoms, molecules, and solids, Phys. Rev. B 70, 235119 (2004),
[20] N. D. Drummond and R. J. Needs, Variance minimization scheme for optimizing Jastrow factors, Phys. Rev. B 72, 085124 (2005),
[21] C.J. Umrigar, J. Toulouse, C. Filippi, S. Sorella, and R. G. Hennig, Alleviation of the fermion-sign problem by optimization of many-body wave functions, Phys. Rev. Lett. 98, 110201 (2007).

Full project description

download file