Quantum circuit simulation with a local time-dependent variational principle

Quantum circuit simulation with a local time-dependent variational principle | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Quantum circuit simulation with a local time-dependent variational principle Jens Eisert, Aaron Sander, Maximilian Fröhlich, Mazen Ali, Martin Eigel, and 4 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8196476/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted You are reading this latest preprint version Abstract Classical simulations of quantum circuits are vital for assessing potential quantum advantage and benchmarking devices, yet they require sophisticated methods to avoid the exponential growth of resources. Tensor network approaches, in particular matrix product states (MPS) combined with the time-evolving block decimation (TEBD) algorithm, currently dominate large-scale circuit simulations. These methods scale efficiently when entanglement is limited but suffer rapid bond dimension growth with increasing entanglement and handle long-range gates via costly SWAP insertions. Motivated by the success of the time-dependent variational principle (TDVP) in many-body physics, we reinterpret quantum circuits as series of discrete time evolutions, using gate generators to construct an MPS-based circuit simulation via a local TDVP formulation. This addresses TEBD's key limitations by (1) naturally accommodating long-range gates and (2) optimally representing states on the MPS manifold. By diffusing entanglement more globally, the method suppresses local bond growth and reduces memory and runtime costs. We benchmark the approach on five 49-qubit circuits: three Hamiltonian circuits (1D open and periodic Heisenberg, 2D 7x7 Ising) and two algorithmic ones (quantum approximate optimization, hardware-efficient ansatz). Across all cases, our method yields substantial resource reductions over standard tools, establishing a new state-of-the-art for circuit simulation and enabling advances across quantum computing, condensed matter, and beyond. Physical sciences/Physics/Quantum physics/Quantum information Physical sciences/Physics/Quantum physics/Qubits Physical sciences/Mathematics and computing/Computational science Full Text Additional Declarations There is NO Competing Interest. Supplementary Files supplement.pdf Supplementary material Cite Share Download PDF Status: Under Review Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8196476","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":554810257,"identity":"80b2605d-42e6-4708-8fb9-573702814e68","order_by":0,"name":"Jens Eisert","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA40lEQVRIiWNgGAWjYBACxgYQaQAimA8ACQkZErSwsSWAtPCQYB8bD9guwlqY23sPf/hQUCdnPr/n86sbNRY8DOyHj27A67Cec2mSMwwOG8sc491mnXMM6DCetLQbeLXMyDFj5jE4kDiDjXebcQ4bUIsEjxl+LfPfGH/+Y1BXP4ON55lxzj9itMzgMZBmMGBOkGDjYX6c20aMlp4cM8keg8OGM9jSzJhz+yR42Aj5xbD9jPGHH3/q5CWYDz/+nPOtTo6f/fAx/FoaEGw2CTCJTzkIyCOxmT8QUj0KRsEoGAUjEwAAE+NBFcgJGPUAAAAASUVORK5CYII=","orcid":"https://orcid.org/0000-0003-3033-1292","institution":"Freie Universität Berlin","correspondingAuthor":true,"prefix":"","firstName":"Jens","middleName":"","lastName":"Eisert","suffix":""},{"id":554810258,"identity":"47ce2568-13de-47db-a67d-c02188f9ff01","order_by":1,"name":"Aaron Sander","email":"","orcid":"https://orcid.org/0009-0007-9166-6113","institution":"Technical University of Munich","correspondingAuthor":false,"prefix":"","firstName":"Aaron","middleName":"","lastName":"Sander","suffix":""},{"id":554810259,"identity":"d4434b64-f4c4-4b8e-b47d-d7a85af1a473","order_by":2,"name":"Maximilian Fröhlich","email":"","orcid":"https://orcid.org/0009-0007-5276-2858","institution":"WIAS Berlin","correspondingAuthor":false,"prefix":"","firstName":"Maximilian","middleName":"","lastName":"Fröhlich","suffix":""},{"id":554810260,"identity":"3d202654-0ec9-4c85-9d74-057eb4ae17a3","order_by":3,"name":"Mazen Ali","email":"","orcid":"","institution":"Department of Flow and Material Simulation, Fraunhofer ITWM","correspondingAuthor":false,"prefix":"","firstName":"Mazen","middleName":"","lastName":"Ali","suffix":""},{"id":554810261,"identity":"91567357-2893-4489-9cf0-2b973122f4e6","order_by":4,"name":"Martin Eigel","email":"","orcid":"","institution":"WIAS Berlin","correspondingAuthor":false,"prefix":"","firstName":"Martin","middleName":"","lastName":"Eigel","suffix":""},{"id":554810262,"identity":"c0754bc4-6d94-4f29-bee5-317c66361615","order_by":5,"name":"Michael Hintermüller","email":"","orcid":"","institution":"Weierstrass Institute for Applied Analysis and Stochastics","correspondingAuthor":false,"prefix":"","firstName":"Michael","middleName":"","lastName":"Hintermüller","suffix":""},{"id":554810263,"identity":"a4fcd5cb-bbf5-45a1-b1c2-cd96c4deed74","order_by":6,"name":"Christian Mendl","email":"","orcid":"https://orcid.org/0000-0002-6386-0230","institution":"Technical University of Munich","correspondingAuthor":false,"prefix":"","firstName":"Christian","middleName":"","lastName":"Mendl","suffix":""},{"id":554810264,"identity":"f7e74926-611b-44fc-a769-3cabefa50217","order_by":7,"name":"Richard Milbradt","email":"","orcid":"https://orcid.org/0000-0001-8630-9356","institution":"Technical University of Munich","correspondingAuthor":false,"prefix":"","firstName":"Richard","middleName":"","lastName":"Milbradt","suffix":""},{"id":554810265,"identity":"bebc5ac9-ab05-4295-a01d-e07f49254197","order_by":8,"name":"Robert Wille","email":"","orcid":"","institution":"Technical University of Munich (TUM)","correspondingAuthor":false,"prefix":"","firstName":"Robert","middleName":"","lastName":"Wille","suffix":""}],"badges":[],"createdAt":"2025-11-24 19:15:10","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-8196476/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-8196476/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":97498384,"identity":"fedbbc3e-c83e-409f-9fb6-1a81d6d4ae3e","added_by":"auto","created_at":"2025-12-05 05:28:26","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1044572,"visible":true,"origin":"","legend":"","description":"","filename":"main.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8196476/v1/473af19a5bbbb94e564ac85b.pdf"},{"id":97498381,"identity":"982c9bbe-0f5f-40d4-9d49-584c0cac3891","added_by":"auto","created_at":"2025-12-05 05:28:26","extension":"json","order_by":1,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":9811,"visible":true,"origin":"","legend":"","description":"","filename":"NCOMMS2595172.json","url":"https://assets-eu.researchsquare.com/files/rs-8196476/v1/2129509f40f2b75d11d533c9.json"},{"id":97498383,"identity":"29f9d6e1-0c2f-41a7-8f90-800aa3a3390a","added_by":"auto","created_at":"2025-12-05 05:28:26","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":533589,"visible":true,"origin":"","legend":"","description":"","filename":"supplement.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8196476/v1/51c4aad699daaec37eb430c8.pdf"},{"id":97498385,"identity":"c94e51ee-cda7-46b9-808a-bbfef6b0fb91","added_by":"auto","created_at":"2025-12-05 05:28:31","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":387694,"visible":true,"origin":"","legend":"Article File","description":"","filename":"main.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8196476/v1_covered_d9953c8d-f0ad-4478-888b-b30edccdc27d.pdf"},{"id":97498382,"identity":"2f076894-8748-43b5-8645-184930ff49b6","added_by":"auto","created_at":"2025-12-05 05:28:26","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":533589,"visible":true,"origin":"","legend":"Supplementary material","description":"","filename":"supplement.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8196476/v1/ed574ee0238457eb7b2f2717.pdf"}],"financialInterests":"There is \u003cb\u003eNO\u003c/b\u003e Competing Interest.","formattedTitle":"Quantum circuit simulation with a local time-dependent variational principle","fulltext":[],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":true,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":true,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"nature-portfolio","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"","title":"Nature Portfolio","twitterHandle":"","acdcEnabled":false,"dfaEnabled":false,"editorialSystem":"ejp","reportingPortfolio":"","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-8196476/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8196476/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"Classical simulations of quantum circuits are vital for assessing potential quantum advantage and benchmarking devices, yet they require sophisticated methods to avoid the exponential growth of resources. Tensor network approaches, in particular matrix product states (MPS) combined with the time-evolving block decimation (TEBD) algorithm, currently dominate large-scale circuit simulations. These methods scale efficiently when entanglement is limited but suffer rapid bond dimension growth with increasing entanglement and handle long-range gates via costly SWAP insertions. Motivated by the success of the time-dependent variational principle (TDVP) in many-body physics, we reinterpret quantum circuits as series of discrete time evolutions, using gate generators to construct an MPS-based circuit simulation via a local TDVP formulation. This addresses TEBD's key limitations by (1) naturally accommodating long-range gates and (2) optimally representing states on the MPS manifold. By diffusing entanglement more globally, the method suppresses local bond growth and reduces memory and runtime costs. We benchmark the approach on five 49-qubit circuits: three Hamiltonian circuits (1D open and periodic Heisenberg, 2D 7x7 Ising) and two algorithmic ones (quantum approximate optimization, hardware-efficient ansatz). Across all cases, our method yields substantial resource reductions over standard tools, establishing a new state-of-the-art for circuit simulation and enabling advances across quantum computing, condensed matter, and beyond.","manuscriptTitle":"Quantum circuit simulation with a local time-dependent variational principle","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-12-05 05:28:21","doi":"10.21203/rs.3.rs-8196476/v1","editorialEvents":[],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"nature-communications","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"NCOMMS","sideBox":"Learn more about [Nature Communications](http://www.nature.com/ncomms/)","snPcode":"","submissionUrl":"https://mts-ncomms.nature.com/","title":"Nature Communications","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"ejp","reportingPortfolio":"Nature Communications","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"ab7dfc4a-e9c8-4bb4-b28e-d0e0f1fb28d2","owner":[],"postedDate":"December 5th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[{"id":59045653,"name":"Physical sciences/Physics/Quantum physics/Quantum information"},{"id":59045654,"name":"Physical sciences/Physics/Quantum physics/Qubits"},{"id":59045655,"name":"Physical sciences/Mathematics and computing/Computational science"}],"tags":[],"updatedAt":"2026-01-22T10:11:34+00:00","versionOfRecord":[],"versionCreatedAt":"2025-12-05 05:28:21","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-8196476","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8196476","identity":"rs-8196476","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}