@article{21340,
  abstract     = {Equilibrium quantum systems are often described by a gas of weakly interacting normal modes. Bringing such systems far from equilibrium, however, can drastically enhance mode-to-mode interactions. Understanding the resulting liquid is a fundamental question for quantum statistical mechanics and a practical question for engineering driven quantum devices. To tackle this question, we probe the non-equilibrium kinetics of one-dimensional plasmons in a long chain of Josephson junctions. We introduce multimode spectroscopy to controllably study the departure from equilibrium, witnessing the evolution from pairwise coupling between plasma modes at weak driving to dramatic, high-order, cascaded couplings at strong driving. Scaling to many-mode drives, we stimulate interactions between hundreds of modes, resulting in near-continuum internal dynamics. Imaging the resulting non-equilibrium plasmon populations, we then resolve the nonlocal redistribution of energy in the response to a weak perturbation—an explicit verification of the emergence of a strongly interacting, non-equilibrium liquid of plasmons.},
  author       = {Bubis, Anton and Vigliotti, Lucia and Serbyn, Maksym and Higginbotham, Andrew P},
  issn         = {2375-2548},
  journal      = {Science Advances},
  number       = {7},
  publisher    = {American Association for the Advancement of Science},
  title        = {{Non-equilibrium plasmon liquid in a Josephson junction chain}},
  doi          = {10.1126/sciadv.ady7222},
  volume       = {12},
  year         = {2026},
}

@article{21480,
  abstract     = {We present and test a protocol to learn the matrix-product operator (MPO) representation of an experimentally prepared quantum state. The protocol takes as input classical shadows corresponding to local randomized measurements, and outputs the tensors of an MPO maximizing a suitably defined fidelity with the experimental state. The tensor optimization is carried out sequentially, similarly to the well-known density matrix renormalization group algorithm. Our approach is provably efficient under certain technical conditions expected to be met in short-range correlated states and in typical noisy experimental settings. Under the same conditions, we also provide an efficient scheme to estimate fidelities between the learned and the experimental states. We experimentally demonstrate our protocol by learning entangled quantum states of up to N = 96 qubits in a superconducting quantum processor. Our method upgrades classical shadows to large-scale quantum computation and simulation experiments.},
  author       = {Votto, Matteo and Ljubotina, Marko and Lancien, Cécilia and Cirac, J. Ignacio and Zoller, Peter and Serbyn, Maksym and Piroli, Lorenzo and Vermersch, Benoît},
  issn         = {1079-7114},
  journal      = {Physical Review Letters},
  number       = {9},
  publisher    = {American Physical Society},
  title        = {{Learning mixed quantum states in large-scale experiments}},
  doi          = {10.1103/rbg2-f61m},
  volume       = {136},
  year         = {2026},
}

@article{21501,
  abstract     = {Kinetically constrained models were originally introduced to capture slow relaxation in glassy systems, where dynamics are hindered by local constraints instead of energy barriers. Their quantum counterparts have recently drawn attention for exhibiting highly degenerate eigenstates at zero energy—known as zero modes—stemming from chiral symmetry. Yet, the structure and implications of these zero modes remain poorly understood. In this work, we focus on the properties of the zero mode subspace in quantum kinetically constrained models with a U(1) particle-conservation symmetry. We use the U(1) East, which lacks inversion symmetry, and the inversion-symmetric U(1) East-West models to illustrate our two main results. First, we observe that the simultaneous presence of constraints and chiral symmetry generally leads to a parametric increase in the number of zero modes due to the fragmentation of the many-body
Hilbert space into disconnected sectors. Second, we generalize the concept of compact localized states from single-particle physics and introduce the notion of collective bound states, a special kind of nonergodic eigenstates that are robust to enlarging the system size. We formulate sufficient criteria for their existence, arguing that the degenerate zero mode subspace plays a central role, and demonstrate bound states in both example models and in a two-dimensional model, the U(1) North-East, and in the pairflip model, a system without particle conservation. Our results motivate a systematic study of bound states and their relation to ergodicity breaking, transport, and other properties of quantum kinetically constrained
models. },
  author       = {Nicolau Jimenez, Eulalia and Ljubotina, Marko and Serbyn, Maksym},
  issn         = {2691-3399},
  journal      = {PRX Quantum},
  publisher    = {American Physical Society},
  title        = {{Fragmentation, zero modes, and collective bound states in constrained models}},
  doi          = {10.1103/sl79-1xgb},
  volume       = {7},
  year         = {2026},
}

@article{19833,
  abstract     = {Eigenstates of quantum many-body systems are often used to define phases of matter in and out of equilibrium; however, experimentally accessing highly excited eigenstates is a challenging task, calling for alternative strategies to dynamically probe nonequilibrium phases. In this work, we characterize the dynamical properties of a disordered spin chain, focusing on the spin-glass regime. Using tensor-network simulations, we observe oscillatory behavior of local expectation values and bipartite entanglement entropy. We explain these oscillations deep in the many-body localized spin-glass regime via a simple theoretical model. From perturbation theory, we predict the timescales up to which our analytical description is valid and confirm it with numerical simulations. Finally, we study the correlation length dynamics, which, after a long-time plateau, resume growing in line with renormalization group (RG) expectations. Our work suggests that RG predictions can be quantitatively tested against numerical simulations and experiments, potentially enabling microscopic descriptions of dynamical phases in large systems.},
  author       = {Brighi, Pietro and Ljubotina, Marko and Serbyn, Maksym},
  issn         = {2469-9969},
  journal      = {Physical Review B},
  number       = {22},
  publisher    = {American Physical Society},
  title        = {{Probing the many-body localized spin-glass phase through quench dynamics}},
  doi          = {10.1103/9fms-ygfz},
  volume       = {111},
  year         = {2025},
}

@article{19852,
  abstract     = {Technology involving hybrid superconductor–semiconductor materials is a promising avenue for engineering quantum devices for information storage, manipulation, and transmission. Proximity-induced superconducting correlations are an essential part of such devices. While the proximity effect in the conduction band of common semiconductors is well understood, its manifestation in confined hole gases, realized for instance in germanium, is an active area of research. Lower-dimensional hole-based systems, particularly in germanium, are emerging as an attractive platform for a variety of solid-state quantum devices, due to their combination of efficient spin and charge control and long coherence times. The recent experimental realization of the proximity effect in germanium thus calls for a theoretical description that is tailored to hole gases. In this work, we propose a simple model to describe proximity-induced superconductivity in two-dimensional hole gases, incorporating both the heavy-hole (HH) and light-hole (LH) bands. We start from the Luttinger–Kohn model, introduce three parameters that characterize hopping across the superconductor–semiconductor interface, and derive explicit intraband and interband effective pairing terms for the HH and LH bands. Unlike previous approaches, our theory provides a quantitative relationship between induced pairings and interface properties. Restricting our general model to an experimentally relevant case where only the HH band crosses the chemical potential, we predict the coexistence of 𝑠-wave and 𝑑-wave singlet pairings, along with triplet-type pairings, and modified Zeeman and Rashba spin–orbit couplings. Our results thus present a starting point for theoretical modeling of quantum devices based on proximitized hole gases, fueling further progress in quantum technology.},
  author       = {Babkin, Serafim and Joecker, Benjamin and Flensberg, Karsten and Serbyn, Maksym and Danon, Jeroen},
  issn         = {2469-9969},
  journal      = {Physical Review B},
  number       = {21},
  publisher    = {American Physical Society},
  title        = {{Superconducting proximity effect in two-dimensional hole gases}},
  doi          = {10.1103/k4jh-pnxy},
  volume       = {111},
  year         = {2025},
}

@article{20646,
  abstract     = {Describing general quantum many-body dynamics is a challenging task due to the exponential growth of the Hilbert space with system size. The time-dependent variational principle (TDVP) provides a powerful tool to tackle this task by projecting quantum evolution onto a classical dynamical system within a variational manifold. In classical systems, periodic orbits play a crucial role in understanding the structure of the phase space and the long-term behavior of the system. However, finding periodic orbits is generally difficult, and their existence and properties in generic TDVP dynamics over matrix product states have remained largely unexplored. In this work, we develop an algorithm to systematically identify and characterize periodic orbits in TDVP dynamics. Applying our method to the periodically kicked Ising model, we uncover both stable and unstable periodic orbits. We characterize the Kolmogorov-Arnold-Moser tori in the vicinity of stable periodic orbits and track the change of the periodic orbits as we modify the Hamiltonian parameters. We observe that periodic orbits exist at any value of the coupling constant of the kicked Ising model between prethermal and fully thermalizing regimes, but their relevance to quantum dynamics and imprint on quantum eigenstates diminishes as the system leaves the prethermal regime. Our results demonstrate that periodic orbits provide valuable insights into the TDVP approximation of quantum many-body evolution and establish a closer connection between quantum and classical chaos.},
  author       = {Petrova, Elena and Ljubotina, Marko and Yalniz, Gökhan and Serbyn, Maksym},
  issn         = {2691-3399},
  journal      = {PRX Quantum},
  number       = {4},
  publisher    = {American Physical Society},
  title        = {{Finding periodic orbits in projected quantum many-body dynamics}},
  doi          = {10.1103/tldp-kvkd},
  volume       = {6},
  year         = {2025},
}

@article{19664,
  abstract     = {Persistent revivals recently observed in Rydberg atom simulators have challenged our understanding of thermalization and attracted much interest to the concept of quantum many-body scars (QMBSs). QMBSs are non-thermal highly excited eigenstates that coexist with typical eigenstates in the spectrum of many-body Hamiltonians, and have since been reported in multiple theoretical models, including the so-called PXP model, approximately realized by Rydberg simulators. At the same time, questions of how common QMBSs are and in what models they are physically realized remain open. In this Letter, we demonstrate that QMBSs exist in a broader family of models that includes and generalizes PXP to longer-range constraints and states with different periodicity. We show that in each model, multiple QMBS families can be found. Each of them relies on a different approximate algebra, leading to oscillatory dynamics in all cases. However, in contrast to the PXP model, their observation requires launching dynamics from weakly entangled initial states rather than from a product state. QMBSs reported here may be experimentally probed using Rydberg atom simulator in the regime of longer-range Rydberg blockades.},
  author       = {Kerschbaumer, Aron and Ljubotina, Marko and Serbyn, Maksym and Desaules, Jean-Yves Marc},
  issn         = {1079-7114},
  journal      = {Physical Review Letters},
  number       = {16},
  publisher    = {American Physical Society},
  title        = {{Quantum many-body scars beyond the PXP model in Rydberg simulators}},
  doi          = {10.1103/PhysRevLett.134.160401},
  volume       = {134},
  year         = {2025},
}

@article{18488,
  abstract     = {The advancement of quantum simulators motivates the development of a theoretical framework to assist with efficient state preparation in quantum many-body systems. Generally, preparing a target entangled state via unitary evolution with time-dependent couplings is a challenging task and very little is known about the existence of solutions and their properties. In this work we develop a constructive approach for preparing matrix product states (MPS) via continuous unitary evolution. We provide an explicit construction of the operator that exactly implements the evolution of a given MPS along a specified direction in its tangent space. This operator can be written as a sum of local terms of finite range, yet it is in general non-Hermitian. Relying on the explicit construction of the non-Hermitian generator of the dynamics, we demonstrate the existence of a Hermitian sequence of operators that implements the desired MPS evolution with an error that decreases exponentially with the operator range. The construction is benchmarked on an explicit periodic trajectory in a translationally invariant MPS manifold. We demonstrate that the Floquet unitary generating the dynamics over one period of the trajectory features an approximate MPS-like eigenstate embedded among a sea of thermalizing eigenstates. These results show that our construction is not only useful for state preparation and control of many-body systems, but also provides a generic route towards Floquet scars—periodically driven models with quasilocal generators of dynamics that have exact MPS eigenstates in their spectrum.},
  author       = {Ljubotina, Marko and Petrova, Elena and Schuch, Norbert and Serbyn, Maksym},
  issn         = {2691-3399},
  journal      = {PRX Quantum},
  number       = {4},
  publisher    = {American Physical Society},
  title        = {{Tangent space generators of matrix product states and exact floquet quantum scars}},
  doi          = {10.1103/prxquantum.5.040311},
  volume       = {5},
  year         = {2024},
}

@article{18616,
  abstract     = {By patterning an ultrathin layered structure with tiny wells, physicists have created and imaged peculiar states known as quantum scars — revealing behaviour that could be used to boost the performance of electronic devices.},
  author       = {Abanin, Dmitry and Serbyn, Maksym},
  issn         = {1476-4687},
  journal      = {Nature},
  number       = {8040},
  pages        = {825--826},
  publisher    = {Springer Nature},
  title        = {{Quantum scars make their mark in graphene}},
  doi          = {10.1038/d41586-024-03649-y},
  volume       = {635},
  year         = {2024},
}

@article{15367,
  abstract     = {Two-dimensional semiconductor-superconductor heterostructures form the foundation of numerous nanoscale physical systems. However, measuring the properties of such heterostructures, and characterizing the semiconductor in-situ is challenging. A recent experimental study by [Phys. Rev. Lett. 128, 107701 (2022)] was able to probe the semiconductor within the heterostructure using microwave measurements of the superfluid density. This work revealed a rapid depletion of superfluid density in semiconductor, caused by the in-plane magnetic field which in presence of spin-orbit coupling creates so-called Bogoliubov Fermi surfaces. The experimental work used a simplified theoretical model that neglected the presence of non-magnetic disorder in the semiconductor, hence describing the data only qualitatively. Motivated by experiments, we introduce a theoretical model describing a disordered semiconductor with strong spin-orbit coupling that is proximitized by a superconductor. Our model provides specific predictions for the density of states and superfluid density. Presence of disorder leads to the emergence of a gapless superconducting phase, that may be viewed as a manifestation of Bogoliubov Fermi surface. When applied to real experimental data, our model showcases excellent quantitative agreement, enabling the extraction of material parameters such as mean free path and mobility, and estimating g-tensor after taking into account the orbital contribution of magnetic field. Our model can be used to probe in-situ parameters of other superconductor-semiconductor heterostructures and can be further extended to give access to transport properties.},
  author       = {Babkin, Serafim and Higginbotham, Andrew P and Serbyn, Maksym},
  issn         = {2542-4653},
  journal      = {SciPost Physics},
  number       = {5},
  publisher    = {SciPost Foundation},
  title        = {{Proximity-induced gapless superconductivity in two-dimensional Rashba semiconductor in magnetic field}},
  doi          = {10.21468/scipostphys.16.5.115},
  volume       = {16},
  year         = {2024},
}

@article{17493,
  abstract     = {Estimating global properties of many-body quantum systems such as entropy or bipartite entanglement is a notoriously difficult task, typically requiring a number of measurements or classical postprocessing resources growing exponentially in the system size. In this work, we address the problem of estimating global entropies and mixed-state entanglement via partial-transposed (PT) moments and show that efficient estimation strategies exist under the assumption that all the spatial correlation lengths are finite. Focusing on one-dimensional systems, we identify a set of approximate factorization conditions (AFCs) on the system density matrix, which allow us to reconstruct entropies and PT moments from information on local subsystems. This identification yields a simple and efficient strategy for entropy and entanglement estimation. Our method could be implemented in different ways, depending on how information on local subsystems is extracted. Focusing on randomized measurements providing a practical and common measurement scheme, we prove that our protocol requires only polynomially many measurements and postprocessing operations, assuming that the state to be measured satisfies the AFCs. We prove that the AFCs hold for finite-depth quantum-circuit states and translation-invariant matrix-product density operators and provide numerical evidence that they are satisfied in more general, physically interesting cases, including thermal states of local Hamiltonians. We argue that our method could be practically useful to detect bipartite mixed-state entanglement for large numbers of qubits available in today’s quantum platforms.},
  author       = {Vermersch, Benoît and Ljubotina, Marko and Cirac, J. Ignacio and Zoller, Peter and Serbyn, Maksym and Piroli, Lorenzo},
  issn         = {2160-3308},
  journal      = {Physical Review X},
  number       = {3},
  publisher    = {American Physical Society},
  title        = {{Many-body entropies and entanglement from polynomially many local measurements}},
  doi          = {10.1103/physrevx.14.031035},
  volume       = {14},
  year         = {2024},
}

@unpublished{17222,
  abstract     = {The quantum approximate optimization algorithm (QAOA) uses a quantum computer
to implement a variational method with $2p$ layers of alternating unitary
operators, optimized by a classical computer to minimize a cost function. While
rigorous performance guarantees exist for the QAOA at small depths $p$, the
behavior at large depths remains less clear, though simulations suggest
exponentially fast convergence for certain problems. In this work, we gain
insights into the deep QAOA using an analytic expansion of the cost function
around transition states. Transition states are constructed in a recursive
manner: from the local minima of the QAOA with $p$ layers we obtain transition
states of the QAOA with $p+1$ layers, which are stationary points characterized
by a unique direction of negative curvature. We construct an analytic estimate
of the negative curvature and the corresponding direction in parameter space at
each transition state. The expansion of the QAOA cost function along the
negative direction to the quartic order gives a lower bound of the QAOA cost
function improvement. We provide physical intuition behind the analytic
expressions for the local curvature and quartic expansion coefficient. Our
numerical study confirms the accuracy of our approximations and reveals that
the obtained bound and the true value of the QAOA cost function gain have a
characteristic exponential decrease with the number of layers $p$, with the
bound decreasing more rapidly. Our study establishes an analytical method for
recursively studying the QAOA that is applicable in the regime of high circuit
depth.},
  author       = {Medina Ramos, Raimel A and Serbyn, Maksym},
  booktitle    = {arXiv},
  title        = {{A recursive lower bound on the energy improvement of the quantum approximate optimization algorithm}},
  doi          = {10.48550/arXiv.2405.10125},
  year         = {2024},
}

@article{13963,
  abstract     = {The many-body localization (MBL) proximity effect is an intriguing phenomenon where a thermal bath localizes due to the interaction with a disordered system. The interplay of thermal and nonergodic behavior in these systems gives rise to a rich phase diagram, whose exploration is an active field of research. In this paper, we study a bosonic Hubbard model featuring two particle species representing the bath and the disordered system. Using state-of-the-art numerical techniques, we investigate the dynamics of the model in different regimes, based on which we obtain a tentative phase diagram as a function of coupling strength and bath size. When the bath is composed of a single particle, we observe clear signatures of a transition from an MBL proximity effect to a delocalized phase. Increasing the bath size, however, its thermalizing effect becomes stronger and eventually the whole system delocalizes in the range of moderate interaction strengths studied. In this regime, we characterize particle transport, revealing diffusive behavior of the originally localized bosons.},
  author       = {Brighi, Pietro and Ljubotina, Marko and Abanin, Dmitry A. and Serbyn, Maksym},
  issn         = {2469-9969},
  journal      = {Physical Review B},
  number       = {5},
  publisher    = {American Physical Society},
  title        = {{Many-body localization proximity effect in a two-species bosonic Hubbard model}},
  doi          = {10.1103/physrevb.108.054201},
  volume       = {108},
  year         = {2023},
}

@article{14320,
  abstract     = {The development of two-dimensional materials has resulted in a diverse range of novel, high-quality compounds with increasing complexity. A key requirement for a comprehensive quantitative theory is the accurate determination of these materials' band structure parameters. However, this task is challenging due to the intricate band structures and the indirect nature of experimental probes. In this work, we introduce a general framework to derive band structure parameters from experimental data using deep neural networks. We applied our method to the penetration field capacitance measurement of trilayer graphene, an effective probe of its density of states. First, we demonstrate that a trained deep network gives accurate predictions for the penetration field capacitance as a function of tight-binding parameters. Next, we use the fast and accurate predictions from the trained network to automatically determine tight-binding parameters directly from experimental data, with extracted parameters being in a good agreement with values in the literature. We conclude by discussing potential applications of our method to other materials and experimental techniques beyond penetration field capacitance.},
  author       = {Henderson, Paul M and Ghazaryan, Areg and Zibrov, Alexander A. and Young, Andrea F. and Serbyn, Maksym},
  issn         = {2469-9969},
  journal      = {Physical Review B},
  number       = {12},
  publisher    = {American Physical Society},
  title        = {{Deep learning extraction of band structure parameters from density of states: A case study on trilayer graphene}},
  doi          = {10.1103/physrevb.108.125411},
  volume       = {108},
  year         = {2023},
}

@article{14334,
  abstract     = {Quantum kinetically constrained models have recently attracted significant attention due to their anomalous dynamics and thermalization. In this work, we introduce a hitherto unexplored family of kinetically constrained models featuring conserved particle number and strong inversion-symmetry breaking due to facilitated hopping. We demonstrate that these models provide a generic example of so-called quantum Hilbert space fragmentation, that is manifested in disconnected sectors in the Hilbert space that are not apparent in the computational basis. Quantum Hilbert space fragmentation leads to an exponential in system size number of eigenstates with exactly zero entanglement entropy across several bipartite cuts. These eigenstates can be probed dynamically using quenches from simple initial product states. In addition, we study the particle spreading under unitary dynamics launched from the domain wall state, and find faster than diffusive dynamics at high particle densities, that crosses over into logarithmically slow relaxation at smaller densities. Using a classically simulable cellular automaton, we reproduce the logarithmic dynamics observed in the quantum case. Our work suggests that particle conserving constrained models with inversion symmetry breaking realize so far unexplored dynamical behavior and invite their further theoretical and experimental studies.},
  author       = {Brighi, Pietro and Ljubotina, Marko and Serbyn, Maksym},
  issn         = {2542-4653},
  journal      = {SciPost Physics},
  keywords     = {General Physics and Astronomy},
  number       = {3},
  publisher    = {SciPost Foundation},
  title        = {{Hilbert space fragmentation and slow dynamics in particle-conserving quantum East models}},
  doi          = {10.21468/scipostphys.15.3.093},
  volume       = {15},
  year         = {2023},
}

@article{12790,
  abstract     = {Motivated by the recent discoveries of superconductivity in bilayer and trilayer graphene, we theoretically investigate superconductivity and other interaction-driven phases in multilayer graphene stacks. To this end, we study the density of states of multilayer graphene with up to four layers at the single-particle band structure level in the presence of a transverse electric field. Among the considered structures, tetralayer graphene with rhombohedral (ABCA) stacking reaches the highest density of states. We study the phases that can arise in ABCA graphene by tuning the carrier density and transverse electric field. For a broad region of the tuning parameters, the presence of strong Coulomb repulsion leads to a spontaneous spin and valley symmetry breaking via Stoner transitions. Using a model that incorporates the spontaneous spin and valley polarization, we explore the Kohn-Luttinger mechanism for superconductivity driven by repulsive Coulomb interactions. We find that the strongest superconducting instability is in the p-wave channel, and occurs in proximity to the onset of Stoner transitions. Interestingly, we find a range of densities and transverse electric fields where superconductivity develops out of a strongly corrugated, singly connected Fermi surface in each valley, leading to a topologically nontrivial chiral p+ip superconducting state with an even number of copropagating chiral Majorana edge modes. Our work establishes ABCA-stacked tetralayer graphene as a promising platform for observing strongly correlated physics and topological superconductivity.},
  author       = {Ghazaryan, Areg and Holder, Tobias and Berg, Erez and Serbyn, Maksym},
  issn         = {2469-9969},
  journal      = {Physical Review B},
  number       = {10},
  publisher    = {American Physical Society},
  title        = {{Multilayer graphenes as a platform for interaction-driven physics and topological superconductivity}},
  doi          = {10.1103/PhysRevB.107.104502},
  volume       = {107},
  year         = {2023},
}

@article{12839,
  abstract     = {Universal nonequilibrium properties of isolated quantum systems are typically probed by studying transport of conserved quantities, such as charge or spin, while transport of energy has received considerably less attention. Here, we study infinite-temperature energy transport in the kinetically constrained PXP model describing Rydberg atom quantum simulators. Our state-of-the-art numerical simulations, including exact diagonalization and time-evolving block decimation methods, reveal the existence of two distinct transport regimes. At moderate times, the energy-energy correlation function displays periodic oscillations due to families of eigenstates forming different su(2) representations hidden within the spectrum. These families of eigenstates generalize the quantum many-body scarred states found in previous works and leave an imprint on the infinite-temperature energy transport. At later times, we observe a long-lived superdiffusive transport regime that we attribute to the proximity of a nearby integrable point. While generic strong deformations of the PXP model indeed restore diffusive transport, adding a strong chemical potential intriguingly gives rise to a well-converged superdiffusive exponent z≈3/2. Our results suggest constrained models to be potential hosts of novel transport regimes and call for developing an analytic understanding of their energy transport.},
  author       = {Ljubotina, Marko and Desaules, Jean Yves and Serbyn, Maksym and Papić, Zlatko},
  issn         = {2160-3308},
  journal      = {Physical Review X},
  number       = {1},
  publisher    = {American Physical Society},
  title        = {{Superdiffusive energy transport in kinetically constrained models}},
  doi          = {10.1103/PhysRevX.13.011033},
  volume       = {13},
  year         = {2023},
}

@article{13125,
  abstract     = {The quantum approximate optimization algorithm (QAOA) is a variational quantum algorithm, where a quantum computer implements a variational ansatz consisting of p layers of alternating unitary operators and a classical computer is used to optimize the variational parameters. For a random initialization, the optimization typically leads to local minima with poor performance, motivating the search for initialization strategies of QAOA variational parameters. Although numerous heuristic initializations exist, an analytical understanding and performance guarantees for large p remain evasive.We introduce a greedy initialization of QAOA which guarantees improving performance with an increasing number of layers. Our main result is an analytic construction of 2p + 1 transition states—saddle points with a unique negative curvature direction—for QAOA with p + 1 layers that use the local minimum of QAOA with p layers. Transition states connect to new local minima, which are guaranteed to lower the energy compared to the minimum found for p layers. We use the GREEDY procedure to navigate the exponentially increasing with p number of local minima resulting from the recursive application of our analytic construction. The performance of the GREEDY procedure matches available initialization strategies while providing a guarantee for the minimal energy to decrease with an increasing number of layers p. },
  author       = {Sack, Stefan and Medina Ramos, Raimel A and Kueng, Richard and Serbyn, Maksym},
  issn         = {2469-9934},
  journal      = {Physical Review A},
  number       = {6},
  publisher    = {American Physical Society},
  title        = {{Recursive greedy initialization of the quantum approximate optimization algorithm with guaranteed improvement}},
  doi          = {10.1103/physreva.107.062404},
  volume       = {107},
  year         = {2023},
}

@article{12259,
  abstract     = {Theoretical foundations of chaos have been predominantly laid out for finite-dimensional dynamical systems, such as the three-body problem in classical mechanics and the Lorenz model in dissipative systems. In contrast, many real-world chaotic phenomena, e.g., weather, arise in systems with many (formally infinite) degrees of freedom, which limits direct quantitative analysis of such systems using chaos theory. In the present work, we demonstrate that the hydrodynamic pilot-wave systems offer a bridge between low- and high-dimensional chaotic phenomena by allowing for a systematic study of how the former connects to the latter. Specifically, we present experimental results, which show the formation of low-dimensional chaotic attractors upon destabilization of regular dynamics and a final transition to high-dimensional chaos via the merging of distinct chaotic regions through a crisis bifurcation. Moreover, we show that the post-crisis dynamics of the system can be rationalized as consecutive scatterings from the nonattracting chaotic sets with lifetimes following exponential distributions. },
  author       = {Choueiri, George H and Suri, Balachandra and Merrin, Jack and Serbyn, Maksym and Hof, Björn and Budanur, Nazmi B},
  issn         = {1089-7682},
  journal      = {Chaos: An Interdisciplinary Journal of Nonlinear Science},
  keywords     = {Applied Mathematics, General Physics and Astronomy, Mathematical Physics, Statistical and Nonlinear Physics},
  number       = {9},
  publisher    = {AIP Publishing},
  title        = {{Crises and chaotic scattering in hydrodynamic pilot-wave experiments}},
  doi          = {10.1063/5.0102904},
  volume       = {32},
  year         = {2022},
}

@article{12276,
  abstract     = {Ongoing development of quantum simulators allows for a progressively finer degree of control of quantum many-body systems. This motivates the development of efficient approaches to facilitate the control of such systems and enable the preparation of nontrivial quantum states. Here we formulate an approach to control quantum systems based on matrix product states (MPSs). We compare counterdiabatic and leakage minimization approaches to the so-called local steering problem that consists in finding the best value of the control parameters for generating a unitary evolution of the specific MPS in a given direction. In order to benchmark the different approaches, we apply them to the generalization of the PXP model known to exhibit coherent quantum dynamics due to quantum many-body scars. We find that the leakage-based approach generally outperforms the counterdiabatic framework and use it to construct a Floquet model with quantum scars. We perform the first steps towards global trajectory optimization and demonstrate entanglement steering capabilities in the generalized PXP model. Finally, we apply our leakage minimization approach to construct quantum scars in the periodically driven nonintegrable Ising model.},
  author       = {Ljubotina, Marko and Roos, Barbara and Abanin, Dmitry A. and Serbyn, Maksym},
  issn         = {2691-3399},
  journal      = {PRX Quantum},
  keywords     = {General Medicine},
  number       = {3},
  publisher    = {American Physical Society},
  title        = {{Optimal steering of matrix product states and quantum many-body scars}},
  doi          = {10.1103/prxquantum.3.030343},
  volume       = {3},
  year         = {2022},
}

