MPS sequential state prep algorithm

docs/source/references.bib

@@ -611,6 +611,40 @@ @article{Suzuki1992
     year={1992},
     doi={10.1016/0375-9601(92)90335-J}
 }
+@article{Rupprecht2026,
+    title={Faster matrix product state preparation by exploiting symmetry-induced block-sparsity},
+    author={Felix Rupprecht and Sabine W{"o}lk},
+    journal={arXiv},
+    year={2026},
+    doi={10.48550/arXiv.2605.28489},
+    url={https://arxiv.org/abs/2605.28489}
+}
+@software{Rupprecht2026Zenodo,
+    title={Code and Assets for: Faster matrix product state preparation by exploiting symmetry-induced block-sparsity},
+    author={Felix Rupprecht and Sabine W{"o}lk},
+    year={2026},
+    publisher={Zenodo},
+    doi={10.5281/zenodo.20393500}
+}
+@article{Berry2025,
+    title={Rapid Initial-State Preparation for the Quantum Simulation of Strongly Correlated Molecules},
+    author={Dominic W. Berry and Yu Tong and Tanuj Khattar and Alec White and Tae In Kim and Guang Hao Low and Sergio Boixo and Zhiyan Ding and Lin Lin and Seunghoon Lee and Garnet Kin-Lic Chan and Ryan Babbush and Nicholas C. Rubin},
+    journal={PRX Quantum},
+    volume={6},
+    number={2},
+    pages={020327},
+    year={2025},
+    doi={10.1103/PRXQuantum.6.020327}
+}
+@article{Clements2017,
+    title={An Optimal Design for Universal Multiport Interferometers},
+    author={William R. Clements and Peter C. Humphreys and Benjamin J. Metcalf and W. Steven Kolthammer and Ian A. Walmsley},
+    year={2017},
+    eprint={1603.08788},
+    archivePrefix={arXiv},
+    primaryClass={physics.optics},
+    url={https://arxiv.org/abs/1603.08788},
+}
 @article{Verstraete2005,
     title={Mapping local {Hamiltonians} of fermions onto local {Hamiltonians} of spins},
     author={F. Verstraete and J. I. Cirac},

python/src/qdk_chemistry/algorithms/registry.py

@@ -794,7 +794,11 @@ def _register_python_algorithms():
     from qdk_chemistry.algorithms.propagator import MagnusPropagator  # noqa: PLC0415
     from qdk_chemistry.algorithms.qubit_hamiltonian_solver import DenseMatrixSolver, SparseMatrixSolver  # noqa: PLC0415
     from qdk_chemistry.algorithms.qubit_mapper import QdkQubitMapper  # noqa: PLC0415
-    from qdk_chemistry.algorithms.state_preparation import SparseIsometryGF2XStatePreparation  # noqa: PLC0415
+    from qdk_chemistry.algorithms.state_preparation import (  # noqa: PLC0415
+        MPSSequentialStatePreparation,
+        MPSSparseStatePreparation,
+        SparseIsometryGF2XStatePreparation,
+    )
     from qdk_chemistry.algorithms.term_grouper import (  # noqa: PLC0415
         FullCommutingTermGrouper,
         IdentityTermGrouper,
@@ -804,6 +808,8 @@ def _register_python_algorithms():
 
     register(lambda: QdkEnergyEstimator())
     register(lambda: SparseIsometryGF2XStatePreparation())
+    register(lambda: MPSSequentialStatePreparation())
+    register(lambda: MPSSparseStatePreparation())
     register(lambda: DenseMatrixSolver())
     register(lambda: SparseMatrixSolver())
     register(lambda: QdkQubitMapper())

python/src/qdk_chemistry/algorithms/state_preparation/__init__.py

@@ -9,6 +9,12 @@
 # --------------------------------------------------------------------------------------------
 
 from qdk_chemistry.algorithms.state_preparation.identity import identity_state_prep
+from qdk_chemistry.algorithms.state_preparation.mps_sequential import (
+    MPSSequentialStatePreparation,
+)
+from qdk_chemistry.algorithms.state_preparation.mps_sparse import (
+    MPSSparseStatePreparation,
+)
 from qdk_chemistry.algorithms.state_preparation.sparse_isometry import (
     SparseIsometryGF2XStatePreparation,
 )
@@ -19,6 +25,8 @@
 )
 
 __all__ = [
+    "MPSSequentialStatePreparation",
+    "MPSSparseStatePreparation",
     "SparseIsometryGF2XStatePreparation",
     "StatePreparationFactory",
     "StatePreparationSettings",