Context

An extremely powerful and ubiquitous concept in theoretical physics is that of duality, loosely defined as a non-trivial physical equivalence between distinct theories. The existence of a duality is typically rooted in a common mathematical structure underlying seemingly very different theories. One appeal of dualities is that certain computations that are difficult in a given theory may be much easier in a dual theory. In the context of quantum lattice models, dualities have recently been classified and constructed in terms of tensor networks in one spatial dimension for closed boundary conditions [1,2,3]. This entanglement based approach to dualities has not only allowed for a concrete realisation of the operators implementing the dualities explicitly on the lattice but has also led to an understanding of how quantum phases of matter can be converted into each other.

An integrable quantum system is roughly speaking a system with an extensive number of conserved charges – local operators commuting with the Hamiltonian – to which they owe their exact solvability. Similar to dualities, tensor networks provide a natural language for describing integrable models, where the transfer matrix describing the model can be written as a matrix product operator and the central Yang-Baxter equation can be constructed following the fundamental theorem of matrix product states, a particular class of tensor networks [4].

Goal

Recent attention has been devoted to the fate of integrable models under duality transformations [5]. Both the transfer matrix of the integrable model and the duality operator can be described by a matrix product operator. Moreover, it has recently been proven that the central charges of the integrable model can likewise be encoded as matrix product operators [6]. The approach in [5] applies these duality transformations for integrable models with closed boundary conditions.

The goal of this project is to systematically construct dualities of integrable models on open boundary conditions extending the approach of [5], along with the explicit mapping of their conserved charges as defined in [6]. This project is well suited for a student with a keen interest in mathematical physics and entanglement theory.

References

[1] I. Cirac, D. Perez-Garcia, N. Schuch, F. Verstraete - arxiv: 2011.12127

[2] L. Lootens, C. Delcamp, G. Ortiz, F. Verstraete - arxiv: 2112.09091

[3] L. Lootens, C. Delcamp, F. Verstraete - arxiv: 2211.03777

[4] B. Vancraeynest-De Cuiper, W. Wiesiolek, F. Verstraete - arxiv: 2512.24390

[5] Y. Miao, A. Molnar, N. Jones - arxiv: 2602.17436

[6] P. Fendley, S. Gehrmann, E. Vernier, F. Verstraete - arxiv: 2511.04674

Context

The phase of a quantum wave function is usually a microscopic detail, yet in coherent systems, it can drive macroscopic behavior. The Josephson effect is perhaps the most famous example of this phenomenon, where a phase difference between two coupled fluids results in a measurable supercurrent. While typically associated with solid-state materials, the effect is a general feature of coupled quantum matter. In the context of ultracold atomic gases, this can be realized by coupling two internal hyperfine states of a Bose gas with an externally driven Rabi field. This coupling acts to synchronize the relative phase between the species, a process which reduces to pendulum-like dynamics in the mean-field limit [1].

In one dimension, however, phase coherence is suppressed by strong quantum fluctuations that prevent the formation of a macroscopic condensate. The system instead behaves as a Tomonaga-Luttinger liquid [2], where phase correlations decay algebraically. The low-energy physics of such a system is captured by an effective Sine-Gordon field theory [3] in the presence of weak Rabi coupling. This framework shows that while even an infinitesimal coupling pins the relative phase and opens an excitation gap, the quantum fluctuations strongly oppose this pinning. This leads to a renormalized scaling of the gap that deviates significantly from the prediction of mean-field theory.

Capturing this behavior goes beyond the reach of standard descriptions like the Gross-Pitaevskii Equation (GPE) [4], which treats the system as a classical object and ignores the non-trivial effects of many-body correlations. As an alternative, one may utilize the recently developed multi-component continuous Matrix Product State (cMPS) formalism [5] to resolve these effects. This approach is particularly powerful because the ansatz is parametrized by the bond dimension D, which acts as a direct probe of entanglement. This allows for a systematic study of how many-body correlations modify the Josephson physics as the system moves beyond the D=1 mean-field limit.

References:

[1] S. Raghavan et al., Phys. Rev. A 59, 620 (1999)

[2] M. A. Cazalilla, J. Phys. B 37, S1 (2004)

[3] Thierry Giamarchi. Quantum Physics in One Dimension, volume 121. Clarendon Press, Oxford, 2003.

[4] Lev Pitaevskii and Sandro Stringari. Bose-Einstein condensation and superfluidity, volume 164. Oxford University Press, 2016.

[5] W. Tang, B. Tuybens, and J. Haegeman, arXiv:2512.24998 (2025)

Goal

The goal of this project is to use the cMPS variational framework to characterize the ground-state properties of a Rabi-coupled 1D Bose gas and quantify the departure from mean-field theory directly in the continuum and the thermodynamic limit.

The project will begin with a literature study on one-dimensional quantum fluids. The student will familiarize themselves with the Lieb-Liniger model, Luttinger liquid theory, and the effective Sine-Gordon description of coupled systems. In parallel, the student will become acquainted with the theoretical construction and numerical implementation of the cMPS variational ansatz.

Building on this, the student will use the cMPS toolbox to investigate how quantum fluctuations renormalize the Josephson excitation gap. The focus will be on testing Sine-Gordon predictions through a finite entanglement scaling analysis in the weak coupling regime. This work will also involve quantifying the effect of inter-species interaction in the system.

This project is well suited for a student with a strong interest in theoretical and computatonal many-body physics. No prior knowledge of tensor networks or coupled bosonic systems is required.

Context

Understanding the interplay between interactions and disorder is a central problem in condensed matter physics. While clean one-dimensional bosonic systems are well described by Luttinger liquid theory and exhibit the well-known superfluid–Mott insulator transition, the presence of disorder fundamentally alters their behaviour. Even weak randomness can destabilise long-range coherence and generate new quantum phases that have no analogue in clean systems.

A paradigmatic model capturing this physics is the disordered Bose–Hubbard chain. In addition to the superfluid and Mott insulating phases, disorder induces the Bose glass phase, which is insulating yet gapless and compressible. Unlike conventional phases, the Bose glass is characterised by strong spatial inhomogeneity and the presence of rare regions that locally resemble different phases. As a result, its physical properties are governed by broad distributions of observables and rare fluctuations rather than uniform bulk behaviour.

Renormalisation group [1] approaches suggest that disorder-driven transitions in one dimension may be governed by strong-disorder fixed points with unconventional scaling properties. In these regimes, rare-region (Griffiths) effects play a dominant role, leading to anomalous dynamical behaviour and broad distributions of energy scales. Despite substantial theoretical work, several fundamental questions remain open. In particular, the precise nature of the superfluid–Bose glass transition is still debated, including the universality of its critical behaviour and the role played by rare regions. Understanding how these features manifest in physical observables such as correlation functions, entanglement properties, and dynamical response remains an active area of research.

Numerical simulations are essential for addressing these questions. One-dimensional systems are particularly well suited to tensor network methods [2], which provide efficient representations of weakly entangled quantum states. Matrix product state techniques enable accurate calculations of ground-state properties, correlations, and entanglement for system sizes far beyond exact diagonalisation. Recent algorithmic developments for disordered systems, such as those introduced in [3], further improve the stability and efficiency of tensor network simulations in strongly inhomogeneous settings.

References:

[1] M. Fisher, Phys. Rev. B 40, 546

[2] J. I. Cirac, D. Pérez-García, N. Schuch, F. Verstraete, arXiv:2011.12127

[3] K. Vervoort, W. Tang, N. Bultinck, arXiv:2504.21089

Goal

The goal of this project is to study disordered bosonic chains using tensor network methods, with a focus on the physics of the Bose glass phase and the critical behaviour of disorder-driven transitions.

The project will begin with a literature study introducing the theoretical framework of one-dimensional bosonic systems. The student will familiarise themselves with Luttinger liquid theory, the Bose–Hubbard model, and the superfluid–Mott insulator transition in clean systems. The effects of disorder will then be explored, including the emergence of the Bose glass phase, renormalisation group descriptions of disordered systems, and the role of rare-region physics in low-dimensional quantum matter. In parallel, the student will become acquainted with tensor network techniques for bosonic lattice models.

Building on this background, the student will use the extensive tensor network toolbox for simulations of disordered Bose–Hubbard chains and use them to investigate the properties of the different phases. Particular attention will be given to identifying signatures of superfluid, insulating, and glassy behaviour through correlation functions, compressibility, excitation gaps, and entanglement entropy. The transition between the superfluid and Bose glass phases will be studied in detail, with the aim of probing its scaling behaviour and identifying possible signatures of strong-disorder physics.

This project is well suited for students with a strong interest in theoretical and computational condensed matter physics. It combines concepts from quantum many-body theory, statistical mechanics of disordered systems, and modern numerical methods. No prior knowledge of tensor networks or disordered bosonic systems is required.

Context

One of the most exciting developments in modern quantum physics has been the discovery of topological phases of matter, whose exotic properties are robust against local perturbations and whose classification bridges the fields of topology and quantum physics. Photonic systems have emerged as a powerful and flexible platform for exploring such phases in a highly controllable setting [1]. A particularly fruitful technique in this and many other areas is Floquet engineering, where periodic driving is used to tailor the time-averaged dynamics of a system to realize a Hamiltonian of interest. One compelling application of this idea is the creation of synthetic dimensions: by coupling internal states of the system through a carefully designed driving protocol, one engineers an effective lattice dimension along those internal states. In photonic ring resonators, only a discrete evenly spaced set of frequencies is held by the resonator. Periodic modulation can coherently couple these modes, with each frequency then acting as a site along the synthetic dimension [2]. Extending this construction by coupling a 1D array of resonators adds a real spatial dimension, and the resulting hybrid system supports effective magnetic fields acting on the photons, opening up the possibility of realizing the paradigmatic quantum Hall effect in a photonic setting [2]. Beyond its fundamental interest as a quantum simulation platform, such a system supports unidirectional topological edge states that are inherently robust, making it equally relevant for practical applications in photonics engineering.

Goal

The goal of this project is to build a theoretical and numerical understanding of the quantum Hall strip realized in coupled ring resonators, with a concrete experimental implementation in mind. The student will first work through the theoretical background needed to understand the setup including topological band theory and the bulk-boundary correspondence, Floquet theory, and how synthetic dimensions are constructed in photonic ring resonator systems. With this in hand, the student will numerically model a specific realization of a band structure spectroscopy measurement along the synthetic dimension, after extending the protocol from [3]. A central question is which features of the bulk and edge states show up clearly in such a measurement and what experimental conditions are needed to resolve them. The model is then extended to include practical imperfections such as fabrication disorder, optical loss, and driving errors, to see how robust the results are. This project benefits from close collaboration with the photonics research group who will provide experimental guidance and context throughout. It is well suited for students who are excited by the prospect of connecting abstract ideas from topological quantum matter to a concrete and timely photonic platform.

References:

[1] Ozawa, T. et al., ”Topological photonics”, Rev. Mod. Phys. 91, 015006 (2019).

[2] Yuan, L. et al., ”Synthetic dimension in photonics,” Optica 5, 1396-1405 (2018)

[3] Dutt, A., et al., ”Experimental band structure spectroscopy along a synthetic dimension.” Nat Commun 10, 3122 (2019).

Context

Classifying quantum phases of matters is one of the most significant objectives in condensed matter physics. By understanding the physical properties and their mechanisms, we can design the materials for our purposes. Yet, the macroscopic mechanism is often not obvious from the microscopic models, making it a challenging task for physicists. One strategy is to look at the phase boundaries. There, the systems often exhibit universal properties that reflect the underlying macroscopic models. The renormalisation group (RG) is an idea that captures this: we can classify universality classes of models sharing universal properties at phase boundaries, which are fixed points of RG flows.

RG flows can capture various phase transitions. In particular, fixed points of second-order phase transitions can often be described by a conformal field theory (CFT), which unlocks a powerful toolbox to understand universality classes. First-order phase transitions, on the other hand, do not show the necessary conformal invariance to accurately be captured by a CFT. However, it has been known historically that some weakly-first order phase transitions such as that of the five-state Potts model, exhibit approximate conformal invariance. Recently, such phase transitions has been shown to be still described by an analytic continuation of CFTs, called complex CFTs [2, 3]. These have been shown to be applicable in various condensed matter/statistical physics systems, as well as in high-energy physics.

CFTs are analytically tractable in two dimensions due to conformal symmetry considerations. For this reason, most of the works in the literature studying critical lattice models do so in two dimensions. However, numerics allow us to access various conformal data in other dimensions, notably in 3D. One such development is the fuzzy sphere regularisation, which has proven to detect the fingerprints of CFTs accurately [4], and opens the gateway to realisations of various 3D CFTs. In particular, a possible connection to complex CFTs has been made with the 3-state Potts model in 2+1D [5], where the phase transition is first-order, but shows approximate conformal invariance.

Goal

This project will start off with a literature study into the relevance of complex conformal field theories in condensed matter contexts. The student will first familiarise themselves with real CFT, and how these are related to criticality and renormalisation group. Afterwards, the student is fully equipped to tackle complex CFT. Finally, the student is asked to perform some original work related to [4], which can be approached from various angles, and especially from the numerical point of view. Tensor networks, a most popular theoretical and numerical toolbox that has been developed in our research group to study strongly correlated quantum systems, will come into play on the fuzzy sphere.

The project can be tailored to the wishes and interests of the student in terms of which research topics they want to delve into.

This project is well-suited for students with a keen interest in mathematical, theoretical and computational physics. No prior knowledge on any of the mentioned topics above is required; a large part of this project concerns understanding these connections more deeply.

[1] D. Kaplan, et al., PhysRevD.80.125005

[2] V. Gorbenko, et al., arXiv:1807.11512

[3] V. Vander Linden, et al., arXiv:2507:14732

[4] W. Zhu, et al., PhysRevX.13.021009

[5] S. Yang, et al., arXiv:2501.14320

Context

The search for new functional quantum phases is a central theme in condensed-matter physics, driven by the desire understand how strong electronic interactions and quantum mechanics conspire to produce unexpected collective behaviour. Among the most intriguing of these phases are quantum spin liquids: highly entangled states that evade conventional symmetry breaking even at zero temperature [1]. Unlike ordinary magnets, their electron spins remain in a fluctuating, liquid-like state, stabilised by strong quantum to fluctuations and geometric frustration. The triangular lattice provides a canonical setting for such physics: antiferromagnetic interactions on this geometry cannot all be satisfied simultaneously, giving rise to frustration, competing ground states, and a rich landscape of possible behaviours, including Mott insulating, metallic, and spin-liquid regimes.

Since 1991, a prominent material candidate realising this landscape is κ-(BEDT-TTF)₂Cu₂(CN)₃, an organic charge-transfer salt whose low-energy electronic structure is well described by half-filled bands on an anisotropic triangular lattice of BEDT-TTF dimers (Fig. 1) [2]. At ambient pressure, it exhibits a Mott insulating state widely regarded as a quantum spin-liquid candidate. Upon applying pressure, the bandwidth increases relative to the on-site Coulomb repulsion, effectively tuning the ratio and driving the system across a metal-insulator transition, with superconductivity emerging at intermediate pressures. However, how pressure alters the effective interactions, and whether the corresponding simplified two-dimensional model accurately captures the resulting phases, remains unresolved, yet crucial for guiding the design of new functional materials.

The reason is that ab initio simulations for these types of materials are notoriously difficult, because they exhibit very strong electron correlations. Indeed, typical methods such as density-functional theory (DFT) fail in the presence of strong correlations and more accurate methods become unwieldy very quickly for realistic materials. Recently, the group of Profs. Frank Verstraete and Jutho Haegeman on the one side and the group of Prof. Veronique Van Speybroeck have teamed up for tackling precisely this question. The strategy is to use advanced ab initio methods such as extensions of DFT, GW-methods and random-phase approximations to derive effective models for the electron bands that exhibit strong correlations, resulting in a (extended) Hubbard model describing these bands. Afterwards, variational entanglement-based methods (known as tensor-network methods) can be used to solve this effective model in the regime of strong electron correlations.

Goal

In this thesis, the goal is to get more insight into both the Hubbard model on the triangular lattice and the physics of strongly-correlated materials with a triangular geometry. This means that we will construct effective models for κ-(BEDT-TTF)₂Cu₂(CN)₃ at different hydrostatic pressures. The resulting triangular-lattice Hubbard model (Fig. 1) will be solved using tensor networks to probe the pressure-driven stability of the different phases.

In the literature, there is a well-established tensor-network method for simulating these types of model, which boils down to placing the triangular-lattice model on a cylindrical geometry and using matrix product state (MPS) techniques for solving this effective one-dimensional model. Although quite successful, this method scales exponentially with the cylinder circumference, and it is unclear whether the physics is accurately represented by this method.

In order to have a more systematic approach, the full two-dimensional model can be solved with tensor networks using the formalism of projected entangled-pair states (PEPS). While previous methods have done this by mapping the triangular to a square lattice, which breaks the lattice symmetry, we will develop a methodology where we tackle the model directly on the infinite triangular lattice [3].

During the course of the thesis, the student will develop a performant framework to simulate the physics of the triangular-lattice Hubbard model, which is still under heavy debate in the community. This involves both code to find the ground state of the Hubbard model and code to calculate its physical properties. If this first part is successful, the student can proceed to simulate the physics of the effective models for κ-(BEDT-TTF)₂Cu₂(CN)₃ at different hydrostatic pressures. Furthermore, the methodology as developed by the student will serve as an integral part of the global strategy for the ab initio simulation of triangular-lattice materials.

[1] Anderson, P. W. (1973). Resonating valence bonds: A new kind of insulator? Materials Research Bulletin, 8(2), 153–160. https://doi.org/10.1016/0025-5408(73)90167-0

[2] Pustogow, A. (2022). Thirty-Year Anniversary of κ-(BEDT-TTF)2Cu2(CN)3: Reconciling the Spin Gap in a Spin-Liquid Candidate. Solids, 3(1), 93-110. https://doi.org/10.3390/solids3010007

[3] Naumann, J., Eisert, J., & Schmoll, P. (2026). Variational optimization of projected entangled-pair states on the triangular lattice. Physical Review B, 113(4). https://doi.org/10.1103/g5gm-tzf8

Context

Solving strongly correlated quantum systems is notoriously difficult, but nevertheless very important to unravel the mysteries behind high-temperature superconductvitiy and strange quantum phases like spin liquids. One of the most widely-used models to describe the properties of stronlgy-correlated materials is the Hubbard model. It describes fermions hopping on a lattice, with a hopping term and a repulsive term [1]

\(H = -t \sum_{<i j> \sigma} \left(\hat{c}_{i \sigma}^\dagger \hat{c}_{j \sigma} + h.c.\right) + U \sum_i \hat{n}_i \hat{n}_j\).

However, solving this model remains a significant challenge, leading to the development of numerous computational approaches. Among these, Tensor Networks (TN) have emerged as a powerful tool, explicitly capturing strong correlations without suffering from the sign problem (like Quantum Monte Carlo methods). In particular, the two-dimensional TN variant, Projected Entangled-Pair States (PEPS), is widely used but has been primarily applied to ground-state (0K) calculations and static properties. To get more insight into the Hubbard model, it is crucial to extend these methods to finite-temperature states and their time evolution [2]. While a lot of progress has been made in the former, the dynamics of these quantum states is much less investigated. Understanding what happens to a system when there is a sudden quench or when a perturbation on top of the ground state evolves through time nonetheless contains crucial information about the underlying properties of the system.

To investigate the real-time evolution of a quantum state, the evolution operator has to be calculated, which is given by

Since this is exponentially hard in the system size, approximations are necessary. The traditional approach relies on the Suzuki-Trotter decomposition, but a more recent and accurate alternative is to use cluster expansions (see Figure 2) [3,4]. For both of these methods, the key challenge is to systemetically find accurate representations of this evolution operator at high values of \(t\). The cluster expansion method has up until now mostly been used for thermal states in an imaginary-time evolution approach, whereas here real-time evolution would be needed [5].

Goal

The thesis will explore different methods for simulating real-time evolution of a quantum system. Since the state-of-the-art methods are scarce in (2+1)-dimensional systems, the first step will consist of benchmarking the current methods on a simpler model. Depending on the interests of the student, the next step can be to either apply the methods used for the simple model to the tJ and Hubbard model to investigate the properties of these models and the materials they describe, or focus more on the development and implementation of different schemes that could be more performant.

[1] Qin M. (2021) arxiv:2104.00064

[2] Sinha A. (2022) arxiv:2209.00985

[3] Vanhecke B. (2021) arxiv:2112.01507

[4] Vanhecke B. (2019) arxiv:1912.10512

[5] De Meyer S. (2026) arxiv:2602.22113

Context

High-temperature superconductivity was first discovered in 1986 when a system of barium, lanthanum, copper, and oxide exhibited a critical temperature of 30K—then the highest recorded. However, BCS theory, the original framework for superconductivity, failed to explain this phenomenon. A complete theory of this effect is still missing, more than 40 years after its discovery.

To study and simulate superconductivity, the Hubbard model was introduced as a minimal model. It describes fermions hopping on a (in the case of the cuprates: square) lattice, with a hopping term and a repulsive term [1]

\(H = -t \sum_{<i j> \sigma} \left(\hat{c}_{i \sigma}^\dagger \hat{c}_{j \sigma} + h.c.\right) + U \sum_i \hat{n}_i \hat{n}_j\).

However, solving it remains a significant challenge, leading to the development of numerous computational approaches. Among these, Tensor Networks (TN) have emerged as a powerful tool, explicitly capturing strong correlations without suffering from the sign problem (like Quantum Monte Carlo methods). In particular, the two-dimensional TN variant, Projected Entangled-Pair States (PEPS), is widely used but has been primarily applied to ground-state (0K) calculations and static properties. To get more insight into the Hubbard model, it is crucial to extend these methods to finite-temperature states and their time evolution [2]. This would allow us to accuretly predict the finite-temperature phase diagram of the Hubbard model (Figure 1).

A key computational challenge in the latter is evaluating the exponential of the Hamiltonian. For finite temperatures, this requires computing the density operator:

Since this is exponentially hard in the system size, approximations are necessary. The traditional approach relies on the Suzuki-Trotter decomposition, but a more recent and accurate alternative is to use cluster expansions (see Figure 1) [3,4]. For both of these methods, the key challenge is to systemetically find accurate representations of this thermal state at high values of \(\beta\) (low temperatures).

Goal

This thesis will explore different methods for probing the finite-temperature phase diagram of the Hubbard model. There are multiple directions that could be taken depending on the interests of the student. One possibility is using the existing methodologies to study various parameter regimes of the Hubbard model, looking both at the feasibility and computational cost of obtaining accurate results. Another is working further on the development of these algorithms to further optimize their efficiency and accuracy. In both approaches, the first step will consist of looking at the strong-coupling limit of the Hubbard model, where the system can be modeled by the computationally less challenging tJ-model. Once the results of this model (which are interesting in their own right) are benchmarked, several different aspects of the Hubbard model can be investigated.

[1] Qin M. (2021) arxiv:2104.00064

[2] Sinha A. (2022) arxiv:2209.00985

[3] Vanhecke B. (2021) arxiv:2112.01507

[4] Vanhecke B. (2019) arxiv:1912.10512

[5] Czarnik P. (2018) arxiv:1811.05497

Context

Many interesting problems in and around physics can be formulated as the contraction of an infinite two-dimensional tensor network [1].

• Partition function of 2d classical models and (1+1)D quantum systems.
• Ground state entropy of 1d quantum mechanical systems.
• 2- (or more) point correlation functions.
• Counting problems like the perfect matching problem / dimer covering problem.

However, contracting an infinite tensor network exactly would require infinite computational resources. To address this, researchers have developed various approximate contraction techniques that achieve surprisingly high accuracy.

One common approach uses the ideas of the renormalization group to perform coarse-graining on tensors within a network. These techniques - Tensor Network Renormalisation (TNR) techniques - are well-established and serve as a useful benchmark for this thesis.

A promising alternative approach to contract 2d tensor networks, that borrows from statistical physics, is belief propagation (BP) [2]. It has been shown that BP is equivalent with a mean-field approximation to the problem. To increase its accuracy, researchers have put forward several improvements. Loop series expansions and loop cluster expansions to name a few [3, 4].

These novel methods remain an active field of research with many open questions to be answered and possible improvements to discovered.

Goal

This thesis will have a major numerical component, requiring programming in the Julia programming language (no prior knowledge required). It consists of three main objectives:

1. Implement belief propagation Implement belief propagation for the contraction of infinite 2D tensor networks and validate it against exact analytical results. We will make use of the TensorKit Julia library which makes writing tensor network code a breeze.

2. Systematic Benchmarking Against TNR Methods Perform a systematic comparison of BP-based contraction against established techniques (e.g. TRG, HOTRG) in terms of accuracy, computational cost, and scaling with bond dimension. Identify regimes where BP is competitive, and where it breaks down.

3. BP Improvements Implement several improvements on top of BP put forward by [3] and [4].

4. Exploration Explore novel improvements for BP and/or come up with hybrid TNR/BP methods.

References

[1] Bridgeman, J. (2016) arxiv.org:1603.03039

[2] Alkabetz R. (2020) arXiv:2008.04433

[3] Evenbly, G. (2024) arxiv.org:2409.03108

[4] Gray J. (2025) arXiv:2510.05647

Context

In classical circuit, one can use a series of binary code (i.e., ‘0’ and ‘1’) to represent the information and use logic gate to realize all kinds of operations. While in quantum circuit, the binary code is replaced by the quantum state (i.e., ‘｜0>’ and ‘｜1>’ in Fock space) and the gate operation becomes local unitary evolution. From a product state, the system can generate long-range entanglement after a long enough unitary evolution, while the measurement (local Hermitian projection) could reduce it. Thus, one should expect there is a Measurement-Induced Phase Transitions (MIPT) [1-3]. However, most of the recent quantum circuit research focus on the bosonic states and operations (i.e., hard core bosons) and fermionic quantum circuit which shares the same Hilbert space seems less attractive. The reason can come from the redundancy of Jordan-Wigner string in a matrix product operator (MPO) which is the main tool used. If one is familiar with the Projective Quantum Monte Carlo (PQMC) [5], he/she should find out it is a perfect tool to study fermionic quantum circuit since there is an isomorphic mapping. The decoupled repulsive interaction corresponds to the unitary evolution, and the kinetic term becomes the Hermitian projection. The interesting observations in bosonic quantum circuit like MIPT should also be expected in a fermionic quantum circuit with a potentially different critical behavior.

Goal

The goal of this thesis is to explore the possibility of simulating a fermionic quantum circuit within PQMC framework. A first step could be comparing the one-dimensional critical behavior of MIPT between fermionic and bosonic random quantum circuit. Due to PQMC has no dimensional limitation from entanglement area law, one can also study the critical behavior in higher dimension (e.g., two-dimension) fermion quantum circuit.

References

[1] Fisher, M.P., Khemani, V., Nahum, A. and Vijay, S., arXiv:2207.14280

[2] Skinner B, Ruhman J, Nahum A., arXiv:1808.05953

[3] Li Y, Chen X, Fisher M.P., arXiv:1808.06134

[4] Li Y, Chen X, Fisher M.P., arXiv:1901.08092

[5] F Assaad, H Evertz, Computational many-particle physics