@article{18923, abstract = {Combinatorial optimization is a challenging problem applicable in a wide range of fields from logistics to finance. Recently, quantum computing has been used to attempt to solve these problems using a range of algorithms, including parameterized quantum circuits, adiabatic protocols, and quantum annealing. These solutions typically have several challenges: 1) there is little to no performance gain over classical methods; 2) not all constraints and objectives may be efficiently encoded in the quantum ansatz; and 3) the solution domain of the objective function may not be the same as the bit strings of measurement outcomes. This work presents “nonnative hybrid algorithms”: a framework to overcome these challenges by integrating quantum and classical resources with a hybrid approach. By designing nonnative quantum variational anosatzes that inherit some but not all problem structure, measurement outcomes from the quantum computer can act as a resource to be used by classical routines to indirectly compute optimal solutions, partially overcoming the challenges of contemporary quantum optimization approaches. These methods are demonstrated using a publicly available neutral-atom quantum computer on two simple problems of Max k-Cut and maximum independent set. We find improvements in solution quality when comparing the hybrid algorithm to its “no quantum” version, a demonstration of a “comparative advantage.”}, author = {Wurtz, Jonathan and Sack, Stefan and Wang, Sheng-Tao}, issn = {2689-1808}, journal = {IEEE Transactions on Quantum Engineering}, pages = {1--14}, publisher = {Institute of Electrical and Electronics Engineers }, title = {{Solving nonnative combinatorial optimization problems using hybrid quantum–classical algorithms}}, doi = {10.1109/tqe.2024.3443660}, volume = {5}, year = {2024}, } @article{15002, abstract = {The lattice Schwinger model, the discrete version of QED in 1 + 1 dimensions, is a well-studied test bench for lattice gauge theories. Here, we study the fractal properties of this model. We reveal the self-similarity of the ground state, which allows us to develop a recurrent procedure for finding the ground-state wave functions and predicting ground-state energies. We present the results of recurrently calculating ground-state wave functions using the fractal Ansatz and automized software package for fractal image processing. In certain parameter regimes, just a few terms are enough for our recurrent procedure to predict ground-state energies close to the exact ones for several hundreds of sites. Our findings pave the way to understanding the complexity of calculating many-body wave functions in terms of their fractal properties as well as finding new links between condensed matter and high-energy lattice models.}, author = {Petrova, Elena and Tiunov, Egor S. and Bañuls, Mari Carmen and Fedorov, Aleksey K.}, issn = {1079-7114}, journal = {Physical Review Letters}, number = {5}, publisher = {American Physical Society}, title = {{Fractal states of the Schwinger model}}, doi = {10.1103/PhysRevLett.132.050401}, volume = {132}, year = {2024}, } @article{15122, abstract = {Quantum computers are increasing in size and quality but are still very noisy. Error mitigation extends the size of the quantum circuits that noisy devices can meaningfully execute. However, state-of-the-art error mitigation methods are hard to implement and the limited qubit connectivity in superconducting qubit devices restricts most applications to the hardware's native topology. Here we show a quantum approximate optimization algorithm (QAOA) on nonplanar random regular graphs with up to 40 nodes enabled by a machine learning-based error mitigation. We use a swap network with careful decision-variable-to-qubit mapping and a feed-forward neural network to optimize a depth-two QAOA on up to 40 qubits. We observe a meaningful parameter optimization for the largest graph which requires running quantum circuits with 958 two-qubit gates. Our paper emphasizes the need to mitigate samples, and not only expectation values, in quantum approximate optimization. These results are a step towards executing quantum approximate optimization at a scale that is not classically simulable. Reaching such system sizes is key to properly understanding the true potential of heuristic algorithms like QAOA.}, author = {Sack, Stefan and Egger, Daniel J.}, issn = {2643-1564}, journal = {Physical Review Research}, number = {1}, publisher = {American Physical Society}, title = {{Large-scale quantum approximate optimization on nonplanar graphs with machine learning noise mitigation}}, doi = {10.1103/PhysRevResearch.6.013223}, volume = {6}, 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}, } @misc{17471, abstract = {Mechanisms for suppressing thermalization in disorder-free many-body systems, such as Hilbert space fragmentation and quantum many-body scars, have recently attracted much interest in foundations of quantum statistical physics and potential quantum information processing applications. However, their sensitivity to realistic effects such as finite temperature remains largely unexplored. Here, we have utilized IBM's Kolkata quantum processor to demonstrate an unexpected robustness of quantum many-body scars at finite temperatures when the system is prepared in a thermal Gibbs ensemble. We identify such robustness in the PXP model, which describes quantum many-body scars in experimental systems of Rydberg atom arrays and ultracold atoms in tilted Bose--Hubbard optical lattices. By contrast, other theoretical models which host exact quantum many-body scars are found to lack such robustness, and their scarring properties quickly decay with temperature. Our study sheds light on the important differences between scarred models in terms of their algebraic structures, which impacts their resilience to finite temperature.}, author = {Desaules, Jean-Yves Marc}, keywords = {quantum many-body scars, non-equilibrium physics, non-Hermitian physics}, publisher = {Institute of Science and Technology Austria}, title = {{Data for "Enhanced many-body quantum scars from the non-Hermitian Fock skin effect"}}, doi = {10.15479/AT:ISTA:17471}, 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}, } @phdthesis{17208, abstract = {Can current quantum computers provide a speedup over their classical counterparts for some kinds of problems? In this thesis, with a focus on ground state search/preparation, we address some of the challenges that both quantum annealing and variational quantum algorithms suffer from, hindering any possible practical speedup in comparison to the best classical counterparts. In the first part of the thesis, we study the performance of quantum annealing for solving a particular combinatorial optimization problem called 3-XOR satisfability (3-XORSAT). The classical problem is mapped into a ground state search of a 3-local classical Hamiltonian $H_C$. We consider how modifying the initial problem, by adding more interaction terms to the corresponding Hamiltonian, leads to the emergence of a first-order phase transition during the annealing process. This phenomenon causes the total annealing duration, $T$, required to prepare the ground state of $H_C$ with a high probability to increase exponentially with the size of the problem. Our findings indicate that with the growing complexity of problem instances, the likelihood of encountering first-order phase transitions also increases, making quantum annealing an impractical solution for these types of combinatorial optimization problems. In the second part, we focus on the problem of barren plateaus in generic variational quantum algorithms. Barren plateaus correspond to flat regions in the parameter space where the gradient of the cost function is zero in expectation, and with the variance decaying exponentially with the system size, thus obstructing an efficient parameter optimization. We propose an algorithm to circumvent Barren Plateaus by monitoring the entanglement entropy of k-local reduced density matrices, alongside a method for estimating entanglement entropy via classical shadow tomography. We illustrate the approach with the paradigmatic example of the variational quantum eigensolver, and show that our algorithm effectively avoids barren plateaus in the initialization as well as during the optimization stage. Lastly, in the last two Chapters of this thesis, we focus on the quantum approximate optimization algorithm (QAOA), originally introduced as an algorithm for solving generic combinatorial optimization problems in near-term quantum devices. Specifically, we focus on how to develop rigorous initialization strategies with guarantee improvement. Our motivation for this study lies in that for random initialization, the optimization typically leads to local minima with poor performance. Our main result corresponds to the analytical construction of index-1 saddle points or transition states, stationary points with a single direction of descent, as a tool for systematically exploring the QAOA optimization landscape. This leads us to propose a novel greedy parameter initialization strategy that guarantees for the energy to decrease with an increasing number of circuit layers. Furthermore, with precise estimates for the negative Hessian eigenvalue and its eigenvector, we establish a lower bound for energy improvement following a QAOA iteration.}, author = {Medina Ramos, Raimel A}, issn = {2663-337X}, keywords = {Quantum computing, Variational Quantum Algorithms, Optimization}, pages = {133}, publisher = {Institute of Science and Technology Austria}, title = {{Exploring the optimization landscape of variational quantum algorithms}}, doi = {10.15479/at:ista:17208}, 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{14406, abstract = {Recently, a concept of generalized multifractality, which characterizes fluctuations and correlations of critical eigenstates, was introduced and explored for all 10 symmetry classes of disordered systems. Here, by using the nonlinear sigma-model ( NL σ M ) field theory, we extend the theory of generalized multifractality to boundaries of systems at criticality. Our numerical simulations on two-dimensional systems of symmetry classes A, C, and AII fully confirm the analytical predictions of pure-scaling observables and Weyl symmetry relations between critical exponents of surface generalized multifractality. This demonstrates the validity of the NL σ M for the description of Anderson-localization critical phenomena, not only in the bulk but also on the boundary. The critical exponents strongly violate generalized parabolicity, in analogy with earlier results for the bulk, corroborating the conclusion that the considered Anderson-localization critical points are not described by conformal field theories. We further derive relations between generalized surface multifractal spectra and linear combinations of Lyapunov exponents of a strip in quasi-one-dimensional geometry, which hold under the assumption of invariance with respect to a logarithmic conformal map. Our numerics demonstrate that these relations hold with an excellent accuracy. Taken together, our results indicate an intriguing situation: the conformal invariance is broken but holds partially at critical points of Anderson localization.}, author = {Babkin, Serafim and Karcher, Jonas F. and Burmistrov, Igor S. and Mirlin, Alexander D.}, issn = {2469-9969}, journal = {Physical Review B}, number = {10}, publisher = {American Physical Society}, title = {{Generalized surface multifractality in two-dimensional disordered systems}}, doi = {10.1103/PhysRevB.108.104205}, volume = {108}, year = {2023}, } @article{14690, abstract = {Generalized multifractality characterizes system size dependence of pure scaling local observables at Anderson transitions in all 10 symmetry classes of disordered systems. Recently, the concept of generalized multifractality has been extended to boundaries of critical disordered noninteracting systems. Here we study the generalized boundary multifractality in the presence of electron-electron interaction, focusing on the spin quantum Hall symmetry class (class C). Employing the two-loop renormalization group analysis within the Finkel'stein nonlinear sigma model, we compute the anomalous dimensions of the pure scaling operators located at the boundary of the system. We find that generalized boundary multifractal exponents are twice larger than their bulk counterparts. Exact symmetry relations between generalized boundary multifractal exponents in the case of noninteracting systems are explicitly broken by the interaction.}, author = {Babkin, Serafim and Burmistrov, I}, issn = {2469-9969}, journal = {Physical Review B}, number = {20}, publisher = {American Physical Society}, title = {{Boundary multifractality in the spin quantum Hall symmetry class with interaction}}, doi = {10.1103/PhysRevB.108.205429}, volume = {108}, 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{13277, abstract = {Recent experimental advances have inspired the development of theoretical tools to describe the non-equilibrium dynamics of quantum systems. Among them an exact representation of quantum spin systems in terms of classical stochastic processes has been proposed. Here we provide first steps towards the extension of this stochastic approach to bosonic systems by considering the one-dimensional quantum quartic oscillator. We show how to exactly parameterize the time evolution of this prototypical model via the dynamics of a set of classical variables. We interpret these variables as stochastic processes, which allows us to propose a novel way to numerically simulate the time evolution of the system. We benchmark our findings by considering analytically solvable limits and providing alternative derivations of known results.}, author = {Tucci, Gennaro and De Nicola, Stefano and Wald, Sascha and Gambassi, Andrea}, issn = {2666-9366}, journal = {SciPost Physics Core}, keywords = {Statistical and Nonlinear Physics, Atomic and Molecular Physics, and Optics, Nuclear and High Energy Physics, Condensed Matter Physics}, number = {2}, publisher = {SciPost Foundation}, title = {{Stochastic representation of the quantum quartic oscillator}}, doi = {10.21468/scipostphyscore.6.2.029}, volume = {6}, year = {2023}, } @phdthesis{12732, abstract = {Nonergodic systems, whose out-of-equilibrium dynamics fail to thermalize, provide a fascinating research direction both for fundamental reasons and for application in state of the art quantum devices. Going beyond the description of statistical mechanics, ergodicity breaking yields a new paradigm in quantum many-body physics, introducing novel phases of matter with no counterpart at equilibrium. In this Thesis, we address different open questions in the field, focusing on disorder-induced many-body localization (MBL) and on weak ergodicity breaking in kinetically constrained models. In particular, we contribute to the debate about transport in kinetically constrained models, studying the effect of $U(1)$ conservation and inversion-symmetry breaking in a family of quantum East models. Using tensor network techniques, we analyze the dynamics of large MBL systems beyond the limit of exact numerical methods. In this setting, we approach the debated topic of the coexistence of localized and thermal eigenstates separated by energy thresholds known as many-body mobility edges. Inspired by recent experiments, our work further investigates the localization of a small bath induced by the coupling to a large localized chain, the so-called MBL proximity effect. In the first Chapter, we introduce a family of particle-conserving kinetically constrained models, inspired by the quantum East model. The system we study features strong inversion-symmetry breaking, due to the nature of the correlated hopping. We show that these models host so-called quantum Hilbert space fragmentation, consisting of disconnected subsectors in an entangled basis, and further provide an analytical description of this phenomenon. We further probe its effect on dynamics of simple product states, showing revivals in fidelity and local observalbes. The study of dynamics within the largest subsector reveals an anomalous transient superdiffusive behavior crossing over to slow logarithmic dynamics at later times. This work suggests that particle conserving constrained models with inversion-symmetry breaking realize new universality classes of dynamics and invite their further theoretical and experimental studies. Next, we use kinetic constraints and disorder to design a model with many-body mobility edges in particle density. This feature allows to study the dynamics of localized and thermal states in large systems beyond the limitations of previous studies. The time-evolution shows typical signatures of localization at small densities, replaced by thermal behavior at larger densities. Our results provide evidence in favor of the stability of many-body mobility edges, which was recently challenged by a theoretical argument. To support our findings, we probe the mechanism proposed as a cause of delocalization in many-body localized systems with mobility edges suggesting its ineffectiveness in the model studied. In the last Chapter of this Thesis, we address the topic of many-body localization proximity effect. We study a model inspired by recent experiments, featuring Anderson localized coupled to a small bath of free hard-core bosons. The interaction among the two particle species results in non-trivial dynamics, which we probe using tensor network techniques. Our simulations show convincing evidence of many-body localization proximity effect when the bath is composed by a single free particle and interactions are strong. We furthter observe an anomalous entanglement dynamics, which we explain through a phenomenological theory. Finally, we extract highly excited eigenstates of large systems, providing supplementary evidence in favor of our findings.}, author = {Brighi, Pietro}, issn = {2663-337X}, pages = {158}, publisher = {Institute of Science and Technology Austria}, title = {{Ergodicity breaking in disordered and kinetically constrained quantum many-body systems}}, doi = {10.15479/at:ista:12732}, year = {2023}, } @phdthesis{14622, abstract = {This Ph.D. thesis presents a detailed investigation into Variational Quantum Algorithms (VQAs), a promising class of quantum algorithms that are well suited for near-term quantum computation due to their moderate hardware requirements and resilience to noise. Our primary focus lies on two particular types of VQAs: the Quantum Approximate Optimization Algorithm (QAOA), used for solving binary optimization problems, and the Variational Quantum Eigensolver (VQE), utilized for finding ground states of quantum many-body systems. In the first part of the thesis, we examine the issue of effective parameter initialization for the QAOA. The work demonstrates that random initialization of the QAOA often leads to convergence in local minima with sub-optimal performance. To mitigate this issue, we propose an initialization of QAOA parameters based on the Trotterized Quantum Annealing (TQA). We show that TQA initialization leads to the same performance as the best of an exponentially scaling number of random initializations. The second study introduces Transition States (TS), stationary points with a single direction of descent, as a tool for systematically exploring the QAOA optimization landscape. This leads us to propose a novel greedy parameter initialization strategy that guarantees for the energy to decrease with increasing number of circuit layers. In the third section, we extend the QAOA to qudit systems, which are higher-dimensional generalizations of qubits. This chapter provides theoretical insights and practical strategies for leveraging the increased computational power of qudits in the context of quantum optimization algorithms and suggests a quantum circuit for implementing the algorithm on an ion trap quantum computer. Finally, we propose an algorithm to avoid “barren plateaus”, regions in parameter space with vanishing gradients that obstruct efficient parameter optimization. This novel approach relies on defining a notion of weak barren plateaus based on the entropies of local reduced density matrices and showcases how these can be efficiently quantified using shadow tomography. To illustrate the approach we employ the strategy in the VQE and show that it allows to successfully avoid barren plateaus in the initialization and throughout the optimization. Taken together, this thesis greatly enhances our understanding of parameter initialization and optimization in VQAs, expands the scope of QAOA to higher-dimensional quantum systems, and presents a method to address the challenge of barren plateaus using the VQE. These insights are instrumental in advancing the field of near-term quantum computation.}, author = {Sack, Stefan}, issn = {2663-337X}, pages = {142}, publisher = {Institute of Science and Technology Austria}, title = {{Improving variational quantum algorithms : Innovative initialization techniques and extensions to qudit systems}}, doi = {10.15479/at:ista:14622}, 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{12111, abstract = {Quantum impurities exhibit fascinating many-body phenomena when the small interacting impurity changes the physics of a large noninteracting environment. The characterisation of such strongly correlated nonperturbative effects is particularly challenging due to the infinite size of the environment, and the inability of local correlators to capture the buildup of long-ranged entanglement in the system. Here, we harness an entanglement-based observable—the purity of the impurity—as a witness for the formation of strong correlations. We showcase the utility of our scheme by exactly solving the open Kondo box model in the small box limit, and thus describe all-electronic dot-cavity devices. Specifically, we conclusively characterize the metal-to-insulator phase transition in the system and identify how the (conducting) dot-lead Kondo singlet is quenched by an (insulating) intraimpurity singlet formation. Furthermore, we propose an experimentally feasible tomography protocol for the measurement of the purity, which motivates the observation of impurity physics through their entanglement build up.}, author = {Stocker, Lidia and Sack, Stefan and Ferguson, Michael S. and Zilberberg, Oded}, issn = {2643-1564}, journal = {Physical Review Research}, number = {4}, publisher = {American Physical Society}, title = {{Entanglement-based observables for quantum impurities}}, doi = {10.1103/PhysRevResearch.4.043177}, volume = {4}, year = {2022}, }