Tensor Networks for Ab Initio Quantum Systems
a
E.M. Stoudenmire Sep 2019 - Jülich This talk:
The many-electron problem
Tensor networks
Intro to quantum chemistry
Chemistry approaches for tensor networks The Many-Electron Problem Behavior of electrons in matter: • continuum problem • three dimensional • strong interaction (repulsion) between electrons
Credit: MARK GARLICK/SCIENCE PHOTO LIBRARY/Getty Images Can simplify various ways: • Born-Oppenheimer approximation (classical nuclei) • ignore most relativistic effects
Figure credit: iqmol.org Then problem simplifies to @ i = Hˆ
AAAC9nicdVFNaxsxEJU3bZO6X0567EV0KfRkdk2hvTQEcsnRhToJWIuZlWdtYa12kbRpjbJ/pfTQkmt+R6/tpf+m2o2hTuwMCB5v3tOTZtJSCmOj6G8n2Hnw8NHu3uPuk6fPnr/o7R+cmqLSHEe8kIU+T8GgFApHVliJ56VGyFOJZ+niuOmfXaA2olCf7bLEJIeZEpngYD016R0KyjIN3LEStBUgaf0f2pqyBVrHhkbU9CNlc7DuZJ2c9MKoH7VFN0G8AiFZ1XCy3/nOpgWvclSWSzBmHEelTVyTyCXWXVYZLIEvYIZjDxXkaBLXfrSmbzwzpVmh/VGWtuy6w0FuzDJPvTIHOzd3ew25rTeubPYhcUKVlUXFb4Kyyk+goM3U6FRo5FYuPQCuhX8r5XPwc7N+tl2m8Asv8hzU1DWzqcdx4hgqU2lsstxlGDMNauY/WN9Wpxo21Ey20jC+3KJurx9sNVDvCAf0niS4mN2X5I1rLr/T+O4GN8HpoB9H/fjTu/DoeLXdPfKKvCZvSUzekyNyQoZkRDj5QX6R3+RP8DX4FvwMrm6kQWfleUluVXD9D8dh9Ys= @t| i | i with Hamiltonian
1 2 1 Hˆ = ˆ† + v(r) ˆ + u(r, r0) ˆ† ˆ† ˆ ˆ r r r r0 r0 r
AAAD7nichVJdb9MwFHUWPkb52AaPvFhEiE2DqgkgeEGatJc9Doluk+psunGd1KrjRLZTVHn5GzyB4JHfwSv8Af4NTlep6doJS5FO7jnn3pscJ6Xg2vR6f70N/9btO3c373XuP3j4aGt75/GJLipFWZ8WolBnCWgmuGR9w41gZ6VikCeCnSbjw4Y/nTCleSE/mWnJ4hwyyVNOwbjSxY73lozA2KMaf8AkVUBtWNuoxoRLc0FyMKMktcq9NypSal6fkyFkGVNLrGCpGbzqEAmJgPMI7+PJ7oLfI4pnIxO3urTd+6uT7YJeoBc1rlpdX7aYvf8tuMrapcatzZaI9fX6YjvodXuzg1dBOAcBmp9j96N/kGFBq5xJQwVoPQh7pYktKMOpYHWHVJqVQMeQsYGDEnKmYzsLuMbPXWWI00K5Rxo8q7YdFnKtp3nilM2O+jrXFNdxg8qk72PLZVkZJunVoLQS2BS4uS14yBWjRkwdAKq42xXTEbiwjLtTLm72mRZ5DnJoyZiZehDGljCpK8WaWfYyCIkCmbkPrJfViYIVNREzaRBerlHP2kdrDdg5ggjfMAkm2U2TnLHlcpmG1xNcBSdRN3zdjT6+CQ4O5+luoqfoGdpFIXqHDtAROkZ9RL2v3i/vt/fHL/0v/jf/+5V0w5t7nqCl4//8B3GvVeQ= 2 r 2 Zr Zrr0 ⇥ ⇤ v(r) AAACw3icdVFda9swFFW8rzb7aLs+9kXMDLoXYzsZSd8KZbDHFpa2EJtyrVynIpLsSXJGcP079tr9rP2byW4G7ZZcEBzOPUfnSjcrBTc2DH/3vGfPX7x8tbPbf/3m7bu9/YP3l6aoNMMJK0ShrzMwKLjCieVW4HWpEWQm8CpbnLX9qyVqwwv1za5KTCXMFc85A+uodHmcSLC3WV7r5tPNvh8Go8E4jsc0DMLPJyejgQNDh8IhjYKwK5+s6/zmoHefzApWSVSWCTBmGoWlTWvQljOBTT+pDJbAFjDHqYMKJJq07qZu6EfHzGheaHeUpR372FGDNGYlM6dsZzT/9lpyU29a2Xyc1lyVlUXFHoLySlBb0PYL6IxrZFasHACmuZuVslvQwKz7qH6i8AcrpAQ1q5MF2mYapXWCylQa26z6zo8SDWruHtg8VWca/lMnopP60d0GdXd9vNFAncOP6ZYkWM63JTnjI5fb6d/F0e3gMg6iQRBfDP3Ts/V2d8gR+UCOSURG5JR8JedkQhj5Tn6Se/LL++ItPO3ZB6nXW3sOyZPymj/YoOE7 = one-body potential
u(r, r0) = two-body interaction AAAC0XicdVHLattAFB2rj6TuI06zKt0MFaUpFCHJATu7QDZdprROApYwV+MrZ/BoJGZGCUYRLd32J7Lptv2e/E1GikuT1r4wcDj3nLmvpBBcG9+/7jgPHj56vLH5pPv02fMXW73tl8c6LxXDEctFrk4T0Ci4xJHhRuBpoRCyROBJMj9s8ifnqDTP5RezKDDOYCZ5yhkYS016r8rdKANzlqSVqj/8he/eT3qu7+3vh8NBSH1vz6LB0AK/Hwz6Qxp4fhsuWcbRZLtzFU1zVmYoDROg9TjwCxNXoAxnAutuVGosgM1hhmMLJWSo46qdoaZvLTOlaa7sk4a27F1HBZnWiyyxyqZJ/W+uIVflxqVJh3HFZVEalOy2UFoKanLaLIROuUJmxMICYIrbXik7AwXM2LV1I4kXLM8ykNMqmqOpx0FcRSh1qbCpVV26QaRAzuyA9X11ouA/dSRaqRtcrlC334crDdQ63JCuqQTns3WVrPGOy970z+HoenAcekHfCz/tuQeHy+tuktfkDdklARmQA/KRHJERYeQr+Ul+kd/OZ2fhfHO+30qdztKzQ+6F8+MGMozmew==
The "electronic structure problem" Form of the one-body potential
Z = nucleus v(r)= a r Ra AAAC93icdVFNbxMxEHWWQkv4SuHIxWKFVA6NdgMSvSAq9cKxVKStiFerWcebWrG9K9sbFDn7W+AC4srv4AoH/k2921S0TTqSpeeZ9/zGM1kpuLFR9K8T3Nm4e29z6373wcNHj5/0tp8em6LSlA1pIQp9moFhgis2tNwKdlpqBjIT7CSbHjT1kxnThhfqk52XLJEwUTznFKxPpb33sx0iwZ5ludP1K/wO7xJTyRQwyTVQ9zmF2i3+M/Auvrwc1Sks6rQXRv2oDbwK4iUI0TIO0+3ONzIuaCWZslSAMaM4Km3iQFtOBau7pDKsBDqFCRt5qEAyk7j2pzV+6TNjnBfaH2Vxm72qcCCNmcvMM5s2zc1ak1xXG1U230scV2VlmaIXRnklsC1wMzY85ppRK+YeANXc94rpGfgRWT/cLlHsCy2kBDV2ZMpsPYoTR5gylWaNl1uEMdGgJv6D9XV2pmGFTURLDePFGnb7/GCtAHtFOMC3OMFscpuTF15R+Z3GNze4Co4H/fh1f/DxTbh/sNzuFnqOXqAdFKO3aB99QIdoiCj6jn6jP+hvMA++Bj+CnxfUoLPUPEPXIvh1Dj3v9b0= a X | |
ZAAACuXicdVHBThsxEHWWUtIALbRHLlZXSJyi3QWpSFyQuHAMKgmo2VU060w2JrZ3ZXtB0ZJP4Fqu/Fb/ps4mByBhJEtPb97zG3vSQnBjg+Bfw9v4tPl5q/mltb2z+/Xb3v73nslLzbDLcpHr2xQMCq6wa7kVeFtoBJkKvEknF/P+zT1qw3N1bacFJhIyxUecgXXU7z+DaLDnB+2gLroKwiXwybI6g/3GSzzMWSlRWSbAmH4YFDapQFvOBM5acWmwADaBDPsOKpBokqqedUYPHTOko1y7oyyt2deOCqQxU5k6pQQ7Nu97c3Jdr1/a0WlScVWUFhVbBI1KQW1O5w+nQ66RWTF1AJjmblbKxqCBWfc9rVjhA8ulBDWs4gnaWT9MqhiVKTXOs6pHP4w1qMw9cPZWnWpYUceilvrh4xp1fX201kCdw4/oB0lwn32U5IyvXG6n4fsNroJe1A6P29HViX9+sdxukxyQn+SIhOQXOSeXpEO6hJGMPJG/5Nk788Abe3cLqddYen6QN+WZ/wdw3G4= 2 ZAAACuXicdVFNT+MwEHUD7EJZlq8jF4sIaU9V3EUCiQsSF46soIBoomriTlNvbSeyHVZV6E/gCtf9W/wb3NAD0DKSpac37/mNPWkhhXVR9NIIlpZXvn1fXWuu/9j4ubm1vXNt89Jw7PBc5uY2BYtSaOw44STeFgZBpRJv0tHZtH9zj8aKXF+5cYGJgkyLgeDgPHV512O9rTBqRXXRecBmICSzuuhtN/7H/ZyXCrXjEqztsqhwSQXGCS5x0oxLiwXwEWTY9VCDQptU9awTeuCZPh3kxh/taM2+d1SgrB2r1CsVuKH93JuSi3rd0g2Ok0roonSo+VvQoJTU5XT6cNoXBrmTYw+AG+FnpXwIBrjz39OMNf7juVKg+1U8QjfpsqSKUdvS4DSreghZbEBn/oGTj+rUwJw6lrU0ZA8L1PX17YUG6h1hm36RBPfZV0ne+M7ld8o+b3AeXLdb7Her/ecwPD2bbXeV7JF98oswckROyTm5IB3CSUYeyRN5Dk4CCIbB3zdp0Jh5dsmHCuwrBR3cbQ== 1
RAAACwnicdVHBThsxEHW2paWhtNAee7G6QuopWgckekSihx5p1QBSvIpmndlgYnsX20sVbfY3ei2/1b+pd8kBSBjJ0tOb9/zGnqxU0vkk+deLXrzcevV6+01/5+3uu/d7+x/OXVFZgSNRqMJeZuBQSYMjL73Cy9Ii6EzhRTY/bfsXt2idLMwvvygx1TAzMpcCfKA41+Cvsrz+2UzYZC9OBklXdB2wFYjJqs4m+707Pi1EpdF4ocC5MUtKn9ZgvRQKmz6vHJYg5jDDcYAGNLq07oZu6EFgpjQvbDjG04596KhBO7fQWVC2Q7qnvZbc1BtXPv+a1tKUlUcj7oPySlFf0PYH6FRaFF4tAgBhZZiViiuwIHz4pz43+FsUWoOZ1nyOvhmztOZoXGWxzaqXMeMWzCw8sHmsziysqbnqpDFbblB31w83GmhwxEP6TBLczp5LCsYHrrBT9nSD6+B8OGCHg+GPo/jkdLXdbfKJfCZfCCPH5IR8J2dkRAQpyR/yl9xF36Lr6CZy99Kot/J8JI8qWv4HJY3ghQ== 1
RAAACwnicdVFNT9wwEPWGftClH0CPvViNkHpaJSkSPSLRQ4+AuoC0jlYT72Rx13aC7VCtsvkbXOFv8W9wwh6AXUay9PTmPb+xJyulsC6K7nvBxpu3795vfuhvffz0+cv2zu6ZLSrDccgLWZiLDCxKoXHohJN4URoElUk8z2ZHbf/8Go0Vhf7r5iWmCqZa5IKD8xRjCtxlltenzTgZb4fRIOqKroJ4CUKyrOPxTu+OTQpeKdSOS7B2FEelS2swTnCJTZ9VFkvgM5jiyEMNCm1ad0M3dM8zE5oXxh/taMc+ddSgrJ2rzCvbIe3LXkuu640ql/9Ka6HLyqHmj0F5JakraPsDdCIMcifnHgA3ws9K+SUY4M7/U59p/M8LpUBPajZD14zitGaobWWwzaoXYcwM6Kl/YPNcnRlYUTPZScN4sUbdXZ+sNVDvCBP6ShJcT19L8sYnLr/T+OUGV8FZMoh/DpKT/fDwaLndTfKNfCc/SEwOyCH5Q47JkHBSkhtyS+6C38G/4Cqwj9Kgt/R8Jc8qWDwAJ+Dghg== 2
ZAAACuXicdVFNT+MwEHUD7EIXlq8jF4sIaU9VUlZapL0gcdkjCNoimqiauNPU1HYi22FVhf4ErnDlb/FvcEMPpS0jWXp6857f2JPkghsbBG81b21949v3za36j+2dn7t7+wdtkxWaYYtlItO3CRgUXGHLcivwNtcIMhHYSUYX037nAbXhmbqx4xxjCaniA87AOur6rnfa2/ODRlAVXQbhDPhkVpe9/dpr1M9YIVFZJsCYbhjkNi5BW84ETupRYTAHNoIUuw4qkGjispp1Qk8c06eDTLujLK3YeUcJ0pixTJxSgh2axd6UXNXrFnZwFpdc5YVFxT6CBoWgNqPTh9M+18isGDsATHM3K2VD0MCs+556pPA/y6QE1S+jEdpJN4zLCJUpNE6zykc/jDSo1D1w8lmdaFhSR6KS+uHjCnV1fXOlgTqH36RfJMFD+lWSM8653E7DxQ0ug3azEZ42mle//fOL2XY3yRE5Jr9ISP6Qc/KPXJIWYSQlT+SZvHh/PfCG3v2H1KvNPIfkU3nmHQnD3G8= 3
RAAACwnicdVFNTxsxEHWWz6aUQjn2YnWF1FO0G5DKEYkeOFLUAFK8imad2eDG9i62lyra7N/gWv4W/6beJQcgYSRLT2/e8xt70kIK66LoqROsrW9sbm1/6H7c+bT7eW//y5XNS8NxwHOZm5sULEqhceCEk3hTGASVSrxOp2dN//oejRW5/u1mBSYKJlpkgoPzFGMK3G2aVZf16Gi0F0a9qC26DOIFCMmiLkb7nUc2znmpUDsuwdphHBUuqcA4wSXWXVZaLIBPYYJDDzUotEnVDl3TQ8+MaZYbf7SjLfvSUYGydqZSr2yGtG97DbmqNyxddpJUQhelQ82fg7JSUpfT5gfoWBjkTs48AG6En5XyWzDAnf+nLtP4l+dKgR5XbIquHsZJxVDb0mCTVc3DmBnQE//A+rU6NbCkZrKVhvF8hbq9vr/SQL0j7NN3kuB+8l6SN75w+Z3Gbze4DK76vfio1/91HJ6eLba7Tb6Sb+Q7ickPckrOyQUZEE4K8kD+kcfgZ/AnuAvsszToLDwH5FUF8/8qM+CH 3
Attraction to classical point nuclei, atomic number Z Form of the two-body potential
1 = nucleus u(r, r0)= r r AAAC9nicdVHLbhMxFHWGVwmvFJZsLEaIIkE0DkhlA6rUDcsikbZSPIo8zp3Uiu0Z2Z6WyJlfQSxAbPkOtrDhb/BMIzVt0itZOjr3HN9XVkphXZL860Q3bt66fWfrbvfe/QcPH/W2Hx/aojIchryQhTnOmAUpNAydcBKOSwNMZRKOstl+kz86BWNFoT+7eQmpYlMtcsGZC9S496HaoYq5kyz3pn51AV+8xO8xzQ3jntR+cZF4vaJZ1ONenPSTNvA6IEsQo2UcjLc73+ik4JUC7bhk1o5IUrrUM+MEl1B3aWWhZHzGpjAKUDMFNvXtoDV+HpgJzgsTnna4ZVcdnilr5yoLyqZLezXXkJtyo8rl71IvdFk50Py8UF5J7ArcbA1PhAHu5DwAxo0IvWJ+wsJ2XNhtl2o444VSTE88nYGrRyT1FLStDDS1/CIm1DA9DQPWl9WZYWtqKltpTBYb1O33g40GHBzxAF9TiZ1Or6sUjCuucFNy9YLr4HDQJ2/6g09v47395XW30FP0DO0ggnbRHvqIDtAQcfQd/UZ/0N/oS/Q1+hH9PJdGnaXnCboU0a//eGP13Q== | 0| = electron
rAAAC0XicdVHBbtQwEPWmFMq2wBZOiIvVCIkLqzggwbFSLxyL6LaVNtFq4p1srbWdyHaKVmkE4tqf4MIVvoe/qZPuYdvdjmTp+c17nhlPVkphXRT97wVbj7YfP9l52t/de/b8xWD/5aktKsNxxAtZmPMMLEqhceSEk3heGgSVSTzL5kdt/uwSjRWFPnGLElMFMy1ywcF5ajJ4nShwF1lem2bC3q9c4skgjIZRF3QdsCUIyTKOJ/u938m04JVC7bgEa8csKl1ag3GCS2z6SWWxBD6HGY491KDQpnU3Q0PfemZK88L4ox3t2FVHDcrahcq8sm3S3s+15KbcuHL557QWuqwcan5bKK8kdQVtP4ROhUHu5MID4Eb4Xim/AAPc+W/rJxq/80Ip0NM6maNrxiytE9S2MtjWqq9ClhjQMz9gc1edGVhTJ7KThuxqg7p7Pt5ooN4RxvSBSnA5e6iSN664/E7Z/Q2ug9N4yD4M468fw8Oj5XZ3yBtyQN4RRj6RQ/KFHJMR4eQH+UP+kn/Bt2AR/Ax+3UqD3tLzityJ4PoG0e7mTg== r 1 2
Coulomb repulsion of electrons Accurate ground state energy of electronic structure problem extremely useful
Example: energy of two atoms, distance R
E R
equilibrium bond length
R binding energy
vibrational frequency Most techniques require discrete and finite system To achieve this can: • integrate out "core" electrons (pseudopotentials) • treat high energy states with approximations such as perturbation theory • project electron motion to certain orbitals
Credit: EVA PAVARINI Regardless of discretization, system becomes a "lattice" with four states per "site"
0, , , Four states of a site are AAAC7nicdVHLbhMxFHWGVwmvFJZsLEZILFA0jipgWakblkUibaV4FN1xblIrfoxsT6NoOr/BCgRLvoMt/AB/gzMN0DbplSwdn3uOj+1blEr6kGW/O8mt23fu3tu5333w8NHjJ73dp0feVk7gUFhl3UkBHpU0OAwyKDwpHYIuFB4X84NV//gMnZfWfAzLEnMNMyOnUkCI1Lj3htc0e80pr0pwzi4ipHxiF+bfrvu39Z/mzbiXZv2sLboJ2BqkZF2H493Ot+gXlUYThALvRywrQ16DC1IobGKKxxLEHGY4itCARp/X7QMb+jIyEzq1Li4TaMtedtSgvV/qIio1hFN/vbcit/VGVZi+y2tpyiqgERdB00rRYOnqt+hEOhRBLSMA4WS8KxWn4ECE+KddbnAhrNZgJjWfY2hGLK85Gl85XGXV5ynjDswsPrC5qi4cbKi5aqUpO9+ibo8fbDXQ6EgH9IYkOJvdlBSNl1xxpuz6BDfB0aDPsj77sJfuH6ynu0OekxfkFWHkLdkn78khGRJBPpMf5Cf5lZTJp+RL8vVCmnTWnmfkSiXf/wB7hvFG { " # "#} Regardless of discretization, electronic structure Hamiltonian takes following discrete form
† † † H = tijcˆi cˆj + VVijklcˆi cˆ cˆk 0 cˆl AAACvnicdVFNb9NAEN0YCiVQaOHIZVULiVPkDZ+9BXHhWCSSRoqtaLwZp0v2w9pdl0au/wRXOPC3+m+6cQNqSzLSSqM3782b2clLKZxPkstOdO/+zoOHu4+6j5/sPX22f/B85ExlOQ65kcaOc3AohcahF17iuLQIKpd4ki8+r+onZ2idMPqbX5aYKZhrUQgOPkDj0bQW3xeyme7HSe9dwo7eM5r0kjboP4StkZis43h60PmTzgyvFGrPJTg3YUnpsxqsF1xi000rhyXwBcxxElINCl1WtwM39FVAZrQwNjztaYveVNSgnFuqPDAV+FN3t7YCN9UmlS8+ZrXQZeVR82ujopLUG7rans6ERe7lMiTArQizUn4KFrgPf9RNNf7gRinQszpdoG8mLKtT1K6yuPKqL2KWWtDzsGBzm51b+I+dypYas4sN7LZ9f6OABkXcp1uc4Gy+zSkIb6jCTf8ejm5PRv0ee9Prf30bDz6tr7tLXpJD8pow8oEMyBdyTIaEE0l+kl/kdzSIikhF5poaddaaF+RWROdX0iPfNg== ijkl j AAACvHicdVHLbtswEKTVV+q+kvTYC1GhQE+GqBhxcgkS9NJjCtROAEswVtTaYUxSAkk5MBT9Q67tqb/VvynluEDS2gsQGMzOcJbcrJTCuij63QmePH32/MXOy+6r12/evtvd2x/ZojIch7yQhbnMwKIUGodOOImXpUFQmcSLbP6l7V8s0FhR6O9uWWKqYKbFVHBwnhq5SS2um8luGPWOjw7j/iGNelE0YDFrQTzoH/Qp80xbIVnX+WSv8yvJC14p1I5LsHbMotKlNRgnuMSmm1QWS+BzmOHYQw0KbVqvxm3oJ8/kdFoYf7SjK/ahowZl7VJlXqnAXdl/ey25qTeu3PQorYUuK4ea3wdNK0ldQdu301wY5E4uPQBuhJ+V8iswwJ3/oW6i8YYXSoHO62SOrhmztE5Q28pgm1XfhiwxoGf+gc1jdWbgP3UiV9KQ3W5Qr66PNxqod4Qx3ZIEi9m2JG984PI7/bs4uh2M4h476MXf+uHp2Xq7O+QD+Ug+E0YG5JR8JedkSDi5JnfkB/kZnAR5MA/UvTTorD3vyaMKFn8ABhzefw== ij 0
AAADaXichVHdbtMwFHYafkb5a+ECBDcWEQIJqWoK0nYzaQgJdjkk2k2qS3XiuqkX24lsZ6jy8i68Ao/DM/ASuG02urUTR7Ly6fs5x85JCsGN7XZ/B43w1u07d3fuNe8/ePjocav9ZGDyUlPWp7nI9UkChgmuWN9yK9hJoRnIRLDjJPu00I/PmDY8V9/svGAjCaniU07Bemrc+nWI9zExpRw7flphu/qQGVhHq+9kAmnKtCeJ4amES2XsTi+Yd5fxTFR4cAH+2+KfUrd6s9Y926REnR+3om6nuyy8CeIaRKiuo3E7+EkmOS0lU5YKMGYYdws7cqAtp4JVTVIaVgDNIGVDDxVIZkZu+W8r/NozEzzNtT/K4iW7nnAgjZnLxDsl2Jm5ri3IbdqwtNO9keOqKC1TdDVoWgpsc7xYFJ5wzagVcw+Aau7viukMNFDr19kkiv2guZSgJo5kzFbDeOQIU6bUbDHLnUcx0aBS/8DqqjvRsOEmYmmN4vMt7mX73tYA9omoh2+YBGfpTZN8cC3ldxpf3+AmGPQ68ftO7+uH6OBjvd0d9BK9Qm9RjHbRATpER6iPaPAs2A8+B18af8J2+Dx8sbI2gjrzFF2pMPoL9/8fzg== ij X Xijkl