Variational reduced density matrix method in the doubly-occupied configuration interaction space using four-particle N-representability conditions: Application to the XXZ model of quantum magnetism

This work deals with the variational determination of the two-particle reduced density matrix (2-RDM) and the energy corresponding to the ground state of N-particle systems within the doubly occupied configuration interaction (DOCI) space. Here, we impose for the first time up to four-particle N-representability constraint conditions in the variational determination of the 2-RDM matrix elements using the standard semidefinite programming algorithms. The energies and 2-RDMs obtained from this treatment and the corresponding computational costs are compared with those arisen from previously reported less restrictive variational methods [D. R. Alcoba et al., J. Chem. Phys. 149, 194105 (2018)] as well as with the exact DOCI values. We apply the different approximations to the one-dimensional XXZ model of quantum magnetism, which has a rich phase diagram with one critical phase and constitutes a stringent test for the method. The numerical results show the usefulness of our treatment to achieve a high degree of accuracy.Fil: Rubio Garcia, Manuel Alejandro. Consejo Superior de Investigaciones Científicas; EspañaFil: Dukelsky, J.. Consejo Superior de Investigaciones Científicas; EspañaFil: Alcoba, Diego Ricardo. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Capuzzi, Pablo. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Oña, Ofelia Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas; ArgentinaFil: Ríos, Elías Daniel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas. Universidad Nacional de La Plata. Facultad de Ciencias Exactas. Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas; ArgentinaFil: Torre, Alicia. Universidad del País Vasco; EspañaFil: Lain, Luis. Universidad del País Vasco; Españ

Variational reduced density matrix method in the doubly-occupied configuration interaction space using four-particle N-representability conditions: Application to the XXZ model of quantum magnetism

11 pags., 5 figs., 1 app.This work deals with the variational determination of the two-particle reduced density matrix (2-RDM) and the energy corresponding to the ground state of N-particle systems within the doubly occupied configuration interaction (DOCI) space. Here, we impose for the first time up to four-particle N-representability constraint conditions in the variational determination of the 2-RDM matrix elements using the standard semidefinite programming algorithms. The energies and 2-RDMs obtained from this treatment and the corresponding computational costs are compared with those arisen from previously reported less restrictive variational methods [D. R. Alcoba et al., J. Chem. Phys. 149, 194105 (2018)] as well as with the exact DOCI values. We apply the different approximations to the one-dimensional XXZ model of quantum magnetism, which has a rich phase diagram with one critical phase and constitutes a stringent test for the method. The numerical results show the usefulness of our treatment to achieve a high degree of accuracy.The authors acknowledge financial support from the Universidad de Buenos Aires (Grant No. 20020150100157BA), the Consejo Nacional de Investigaciones Científicas y Técnicas (Grant Nos. PIP 11220130100377CO, PIP 11220130100311CO, 11220150100442CO, and 2013-1401PCB), the Agencia Nacional de Promoción Científica y Tecnológica, Argentina (Grant No. PICT-201-0381), and the Spanish Ministerio de Ciencia, Innovación y Universidades, and the European Regional Development Fund (FEDER), Grant Nos. FIS2015-63770-P and PGC2018-094180-B-I00. E.R. acknowledges support from the Consejo Nacional de Investigaciones Científicas y Técnicas