##### Prospective Graduate Students / Postdocs

This faculty member is currently not actively recruiting graduate students or Postdoctoral Fellows, but might consider co-supervision together with another faculty member.

- Graduate School
- »
- Prospective Students »
- Research Supervisors »
- Mona Berciu

This faculty member is currently not actively recruiting graduate students or Postdoctoral Fellows, but might consider co-supervision together with another faculty member.

Professor

condensed matter theory

polarons, bipolarons

strongly correlated systems

I am available and interested in collaborations (e.g. clusters, grants).

I am interested in and conduct interdisciplinary research.

I am interested in working with undergraduate students on research projects.

Dissertations completed in 2010 or later are listed below. Please note that there is a 6-12 month delay to add the latest dissertations.

Electron-phonon coupling in insulators : beyond the Migdal limit (2023)

Motivated by the success of the momentum average (MA) approximation for the single-polaron and bipolaron cases, in this dissertation we pursue the first generalization of this technique to the case of weakly-doped insulators beyond the Migdal limit. In search of a suitable starting point, we investigate the ground states of a variety of model Hamiltonians, inspired by real materials such as rare-earth nickelates, polyacetylene, barium bismuthates, and buckminsterfullerene. Our study of the ground state of the 2D rare-earth nickelate sheet shows that it exhibits cuprate-like orbital polarization unimpeded by weak to moderate electron-phonon coupling. This lends further support to the proposal for cuprate-like superconductivity in nickelate thin films and heterostructures. Investigating polyacetylene leads us to develop a semi-analytical technique that can be used to account for zero-point lattice energy in determining the equilibrium geometry of any crystal. The study of a simplified, one-dimensional model of buckminsterfullerene results in the first generalization of MA to a system with two polaron clouds, and a phonon-driven mechanism for exciton dissociation that could be responsible for photocurrent generation in these materials. Because both the nickelates and polyacetylene turn out to be too complex as starting points for MA generalization, we focus on a simplified, s-p chain model of the barium bismuthates. We ultimately succeed in generalizing MA for this band insulator – the first application of its kind. We compare the results to those from determinant quantum Monte Carlo to establish the accuracy of the MA results (with the MA results obtained at lower computational cost). Finally, we share in the development of the first open-sourced Python package implementing automatic MA equation generation via the GGCE method, which allows MA-style methods to be easily used by non-experts – and is expected to expand in the coming years, providing access to many more MA-style methods in the future.

View record

Peierls coupling in multi-orbital metal oxides (2022)

In this thesis we study the effects of the Peierls electron-phonon coupling in multi-band metal oxides systems. In contrast to the more commonly employed Holstein coupling, which is used in single-band models and is momentum-independent, the momentum-dependent Peierls coupling can explicitly treat coupling to multiple bands. We demonstrate the importance of using the Peierls coupling in modelling complex systems by looking at two examples. First, we investigate single polaron physics on a perovskite lattice inspired by BaBiO₃. We find that with Peierls coupling, the ground state momentum of the polaron jumps between high-symmetry points in the Brillouin zone as the coupling strength is increased. Because such sharp transitions are not possible in the Holstein model, it follows that it is not always feasible to map the more complex Peierls model onto the simpler Holstein model. Then, we add the Peierls coupling to Emery's three-band model for cuprate layers and study its effect on the polaron effective mass. We show that although the hole-phonon coupling strength is moderate to strong, it only causes a negligible increase in the effective mass, indicating that the effective coupling to the magnon-dressed quasiparticle is much reduced by the dressing. We explain the reason for this and describe how to treat the difference between lattice coupling to bare holes versus to correlations-dressed quasiparticles. Our results prove that it is essential to choose the proper coupling when describing multi-orbital systems, in order to deduce their properties appropriately. We propose to investigate more systematically the generic nontrivial behaviour brought by the momentum-dependent Peierls coupling by studying a diatomic 1D chain, to reinforce our conclusion.

View record

Peierls Bipolarons and Localization in Solid-State and Molecular Systems (2019)

In this thesis, I investigate the behavior of particles dressed by quantum field excitations and random interactions.First I consider two-carrier states in the Peierls model describing the modulation of the particle hopping due to lattice distortions. I compute the spectral response using the Momentum Average approximation. Combining accurate numerical techniques and analytical arguments, I provide a complete picture of the Peierls bipolarons. It is found that polarons bind into strongly bound yet light bipolarons in the singlet sector, even at large values of the electron-phonon coupling strength. At finite electron fillings, these bipolarons may condense into a high-Tc superconductor. On the other hand, phonons mediate a repulsive interaction in the triplet sector, or equivalently (in one dimension), between two hard-core particles, in which case the ground-state dimers bound by sufficiently attractive bare interactions exhibit two sharp transitions, one of which is the first known example of a self-trapping transition at the two-carrier level. In both situations, phonons mediate pair-hopping effective interactions between the carriers. I further study some aspects of the excited spectrum for the two hard-core particles, a situation relevant to ultracold quantum simulators. It is found that the repulsive phonon-mediated interaction binds a repulsive bipolaron embedded in the excited spectrum.I then turn to the study of quenched randomness in an ultracold molecular plasma. I argue that the quenched ultracold plasma presents an experimental platform for studying quantum many-body physics of disordered systems in the long-time and finite energy-density limits. I analyze an experiment that quenches a plasma of nitric oxide to an ultracold system of Rydberg molecules, ions and electrons that exhibits a long-lived state of arrested relaxation. The qualitative features of this state fail to conform with classical models. I develop a microscopic quantum description for the arrested phase based on an effective many-body spin Hamiltonian that includes both dipole-dipole and van der Waals interactions. This effective model appears to offer a way to envision the essential quantum disordered non-equilibrium physics of this system.This thesis thus examines the quantum many-body response in interacting systems coupled to bosonic fields or in disordered environments.

View record

Temperature-driven spectral weight transfer in doped magnetic insulators (2017)

In this thesis we study the effects of finite temperature (T ) on the single-electron spectral function of doped magnetic insulators. First, we derive the low-temperature correction to the self-energy of a charge carrier injected with parallel spin into a ferromagnetic background which is modeled with both Heisenberg and Ising Hamiltonians so that differences due to gapless versus gapped magnons can be understood. Beside the expected thermal broadening of the T = 0 quasiparticle peak which becomes a resonance inside a continuum, we find that spectral weight istransferred to regions lying outside this continuum, because the carrier and a thermal magnon can bind into a spin-polaron. This work is valid in dimensions d ≥ 2, because it does not include the role of magnetic domains which are important in 1d. We then consider the role of these magnetic domains in 1d systems, for models where spin-polaron formation is impossible. We present Monte Carlo simulations for the spectral function of three related models of a charge carrier that is injected into an Ising chain. Both ferromagnetic and antiferromagnetic coupling between the Ising spins are considered. The interaction between the carrier and the Ising spins is also of Ising type. In two of the models the charge carrier is hosted by a different band, while in the third model it is hosted by the same band as the Ising spins. We find that thecarrier’s spectral function exhibits a distinctive fine structure due to its temporary entrapping inside small magnetic domains, and use these results to construct an accurate (quasi)analytic approximation for low and medium T . While at T = 0, for ferromagnetic order all three models have identical, low-energy quasiparticles, at finite T the low-energy behavior of the first two models remains equivalent, but that of the third model is controlled by rare events due to thermal fluctuations, which transfer spectral weight below the T = 0 quasiparticle peak, generating a pseudogaplike phenomenology. Taken together, our results show that the temperature evolution of the spectral weight of weakly doped magnetic insulators can be very complex.

View record

Extensions Beyond Standard Models (2016)

In this thesis, we investigated a set of theoretical models frequently used in the field of solid state physics. These models describe coupling of charge carriers to bosonic modes such as phonons or magnons. In particular, the Holstein model describes coupling of charge carriers to dispersionless phonons, whereas the Emery model describes coupling of charge carriers to magnons in hole-doped antiferromagnets.For the Holstein-like models, we studied how extending the model of the coupling beyond terms that are merely linear in the lattice distortion affects the ground state properties. Using appropriate extensions of the momentum average approximation, we could show that even small nonlinearities have a dramatic effect on the resulting quasi-particle's properties. We further investigated a particular type of nonlinear coupling, the double-well coupling model. After studying the properties of a single quasi-particle, we also showed that this system allows the formation of bound states between two charge carriers and a phonon cloud, the so-called bi-polaron. In contrast to the linear variation of the Holstein model, the resulting bi-polaron can be strongly bound yet lightweight. For the Emery model, we consider an experimentally relevant extension. The original model describes a single layer of CuO₂, relevant for the hole-doped cuprate superconductors. We consider recently synthesized layers of CuO, which can be viewed as two intercalated layers of CuO₂. The resulting system is similar to CuO₂, but different in important aspects. We use a variational method similar in spirit to MA but applicable to magnons instead of phonons to obtain the system's dispersion and compare it to that of the original CuO₂ layer. We observe a discrepancy between these dispersions that cannot be accounted for with a single-band model that is commonly used to model the CuO₂ dispersion. However, it has been a long-standing question whether or not this and other single-band models are appropriate for the description of cuprate physics. With our study of CuO, we demonstrated how a careful experimental analysis of this system can resolve that question.

View record

Variational studies of dressed quasiparticles' properties and their interactions with external potentials (2015)

In this thesis, we investigated the spectral properties of polaronic quasiparticles resulting from the coupling of a charge carrier to the bosonic excitations of ordered environments. Holstein polarons, describing an electron locally coupled to dispersionless phonons, and spin polarons in hole-doped antiferromagnets are considered.The Green's function of the Holstein polaron in a lattice with various kinds of impurities is calculated using the momentum average approximation. The main finding is that the scenario where the mere effect of the coupling to phonons is to enhance the electron's effective mass is incomplete as the phonons are found to also renormalize the impurity potential the polaron interacts with. This formalism is applied first to the case of a single impurity and the range of parameters where the polaron's ground state is localized is identified. The lifetime of the polaron due to scattering from weak but extended disorder is then studied and it is shown that the renormalization of the disorder potential leads to deviation of the strong coupling results from Fermi's golden rule's predictions. Motivated by the hole-doped cuprate superconductors, the motion of a hole in a two-dimensional Ising antiferromagnet and its binding to an attractive impurity is studied next, based on a variational scheme that allows for configurations with a certain maximum number of spin flips. The role of the magnetic sublattices in determining the symmetry of the resulting bound states is discussed. Next, a more realistic model describing a hole in a CuO₂ layer which retains the O explicitly, is considered. By neglecting the fluctuations of the Cu spins and using a variational principle similar to that of the previous chapter, a semi-analytic solution for the Green's function of the hole in an infinite 2D lattice is constructed. The resulting quasiparticle dispersion shows the proper ground state and other features observed in experiments. The lack of importance of the background spin fluctuations is justified based on the hole's ability to move on the O sublattice without disturbing the Cu spins. Finally, the model is generalized to gauge the importance of other relevant O orbitals.

View record

Polaron physics beyond the Holstein model (2011)

Many condensed matter problems involve a particle coupled to its environment. The polaron, originally introduced to describe electrons in a polarizable medium, describes a particle coupled to a bosonic field. The Holstein polaron model, although simple, including only optical Einstein phonons and an interaction that couples them to the electron density, captures almost all of the standard polaronic properties. We herein investigate polarons that differ significantly from this behaviour. We study a model with phonon-modulated hopping, and find a radically different behaviour at strong couplings. We report a sharp transition, not a crossover, with a diverging effective mass at the critical coupling. We also look at a model with acoustic phonons, away from the perturbative limit, and again discover unusual polaron properties. Our work relies on the Bold Diagrammatic Monte Carlo (BDMC) method, which samples Feynman diagrammatic expansions efficiently, even those with weak sign problems. Proposed by Prokof'ev and Svistunov, it is extended to lattice polarons for the first time here. We also use the Momentum Average (MA) approximation, an analytical method proposed by Berciu, and find an excellent agreement with the BDMC results. A novel MA approximation able to treat dispersive phonons is also presented, along with a new exact solution for finite systems, inspired by the same formalism.

View record

Quantum Transport Through Open Systems (2011)

To study transport properties, one needs to investigate the system of interest when coupled to biased external baths. This requires solving a master equation for this open quantum system. Obtaining this solution is very challenging, especially for large systems. This limits applications of the theories of open quantum systems, especially insofar as studies of transport in large quantum systems, of interest in condensed matter, is concerned.In this thesis, I propose three efficient methods to solve the Redfield equation --- an example of such a master equation. The first is an open-system Kubo formula, valid in the limit of weak bias. The second is a solution in terms of Green's functions, based on a BBGKY (Bogoliubov--Born--Green--Kirkwood--Yvon)-like hierarchy. In the third, the Redfield equation is mapped to a generalized Fokker-Planck equation using the coherent-state representation. All three methods, but especially the latter two, have much better efficiency than direct methods such as numerical integration of the Redfield equation via the Runge-Kutta method. For a central system with a d-dimensional Hilbert space, the direct methods have complexity of d³, while that of the latter two methods is on the order of order of polynomials of log d. The first method, besides converting the task of solving the Redfield equation to solving the much easier Schrödinger's equation, also provides an even more important conceptual lesson: the standard Kubo formula is not applicable to open systems.Besides these general methodologies, I also investigate transport properties of spin systems using the framework of the Redfield equation and with direct methods. Normal energy and spin transport is found for integrable but interacting systems. This conflicts with the well-known conjecture linking anomalous conductivity to integrability, and it also contradicts the relationship, suggested by some, between gapped systems (Jz > Jxy) and normal spin conductivity. I propose a new conjecture, linking anomalous transport to the existence of a mapping of the problem to one for non-interacting particles.

View record

Study of polaron properties using the momentum average approximation (2009)

In this work we present a highly efficient and accurate analytical approximation for the Green's function of a polaron: the momentum average (MA) approximation. It is obtained by summing all of the self-energy diagrams, but with each diagram averaged over the momenta of its free propagators. The result becomes exact for both zero bandwidth and for zero electron-phonon coupling, and is accurate everywhere in the parameter space. The approximation is first used to investigate the Holstein model. A detailed analysis of its accuracy is provided through diagrammatics and spectral weight sum rules. It is shown that the resulting Green's function satisfies exactly the first six spectral weight sum rules, and all higher order sum rules are satisfied with great accuracy. Comparison with numerical data also confirms this accuracy. We then show how to improve the MA approximation by systematically improving the accuracy of the self-energy diagrams in such a way that they can still all be summed efficiently. This allows us to fix some of the problems of the zeroth-order MA approximation. The quantitative agreement with numerical data is also improved.Next, we generalize the MA approximation to study the properties of models with momentum-dependent electron-phonon coupling, and then show that further improvements can be obtained based on variational considerations, using the 1D breathing-mode Hamiltonian as a specific example. For example, by using this variational MA, we obtain ground state energies within at most 0.3% error of the numerical data.Finally, we study the effects of a nearby surface on the spectral weight of a Holstein polaron. The broken translational symmetry is accounted for without any additional approximations, and the resulting inhomogeneous MA approximation continues to be accurate for all coupling strengths. We show that the surface changes properties significantly, with bulk values being recovered only very far away from it. We find that the electron-phonon coupling gives rise to an additional surface potential which is responsible for binding surface states even when they are not normally expected. These results demonstrate that interpretation in terms of bulk properties of spectroscopic data sensitive only to a few surface layers is not straightforward.

View record

Theses completed in 2010 or later are listed below. Please note that there is a 6-12 month delay to add the latest theses.

Effects of the lattice distortion on the magnetic order in rare-earth nickelates (2019)

The low-temperature magnetic order in the rare-earth nickelates is a subject of vigorous debate in the literature. Recent work emphasized the primary role of the electron-phonon coupling for the metal-insulator transition in the nickelates, and suggested that lattice distortions are the driver of the transition, leading to the observed charge order. However, to our knowledge there has been little work on the impact of lattice distortions on the magnetic order, in particular whether distortions favour some orders over others. In this thesis, we study the magnetic order in the nickelates at zero temperature, and investigate whether the breathing-mode lattice distortions select a preferred ground state. An effective two-band Hubbard model for the nickelates is constructed and coupled to the lattice distortions with an on-site Holstein-like term. The distortions are treated semiclassically. Using the Hartree-Fock approximation, we obtain the magnetic phase diagram, then turn on the coupling to the lattice to observe its impact on the various phases. Our model reproduces the earlier work showing the stronger charge disproportionation and insulating behaviour in the phase space due to increased coupling to the lattice. Furthermore, we find numerous 4-site magnetic orders that are self-consistent within the model, including all of the main suggestions in the literature (states such as ↑↑↓↓,↑→↓← and ⇑0⇓0). However, in this model a magnetic order can only couple to the lattice distortions if there is nonzero charge disproportionation. As a result, we find that coupling to the lattice distortions broadly favours the ⇑0⇓0 order in large sectors of the parameter space. Finally, we considered the impact of longer range hopping on the magnetic order: we find that the shape of the density of states, rather than overall bandwidth, primarily determines the magnetic ground state. A van Hove singularity arises even for small 2nd-nearest hopping amplitudes, which results in robust ferromagnetism across most of the phase diagram in a Stoner-like fashion. On the contrary, even small 4th-nearest amplitudes decrease the Fermi level density of states, resulting in ballooning of the metallic phase despite a barely renormalized bandwidth.

View record

Fabrication and Raman study of twisted stacked few layer black phosphorus (2018)

Two-dimensional black phosphorus (BP) has attracted much interest as a promising semiconductor with a tunable direct bandgap of 0.3-2 eV, depending on its thickness. However, BP-based heterostructures, such as van der Waals junctions, and their interlayer couplings have not been well studied. In this work, we first provided a new“film-based” transfer method, which can be used to transfer BP to specific positions. This is a very convenient transfer method for BP, since commonly used dry transfer and wet transfer methods do not work well for BP. Based on this new method, we successfully fabricated twisted stacked few layer BP with different twist angles (7 to 90 degrees) and thicknesses (10 -30 nm). Raman measurements revealed an “abnormal blue shift” for the twisted few-layer BP. Density functional theory (DFT) calculations by Teren Liu provided an explanation for this phenomenon in terms of changes in charge distributions. Based on the calculated charge distributions, we propose that the interlayer coupling in BP is not just through van der Waals interactions. Other interactions, such as weak valence bond, may exist in this system.

View record

Higgs spectroscopy of superconductors: a new method to identify the superconducting gap symmetry (2018)

In this thesis we study the response of a BCS superconductor to an external ultra-fast terahertz electromagnetic field, which we choose so as to mimic the setup in a pump-probe experiment.We begin by considering an optical experimental setup and demonstrate that in an optical pump-probe experiment, the superconducting amplitude Higgs mode can be excited and measured with ultra-fast terahertz pump pulses. Moreover, for an anistropic d-wave superconductor, there are two Higgs mode, one at the usual 2∆ energy and one with a lower energy. The latter can be used to differentiate the d-wave symmetry from isotropic s-wave, by varying the polarization of the pump relative to the sample. For a linearly polarized pump with a vector potential aligned along a d-wave node we find only a single Higgs mode, while for a direction along an antinode we find two Higgs modes.Next, we consider an angle resolved photoemission spectroscopy (ARPES) experiment and derive a new set of equations of motion, for which we can analyze the two-time nonequilibrium Green’s functions. We show that the Higgs mode can also be studied in an ARPES pump-probe experiment. Moreover, we show how an ARPES pump-probe experiment can be used to differentiate between different momentum-dependent nonequilibrium Higgs modes. Our results suggest that in a d-wave superconductor, the second low-energy Higgs mode is of osculating, B1g character, which corresponds to a symmetry breaking along the d-wave nodal lines. Further study of the role of momentum symmetry breaking promises to provide deeper insight into generating new nonequilibrium states.

View record

Large Electropositive Cations as Surfactants for the Growth of Polar Epitaxial Films (2015)

Using density functional theory (DFT) we demonstrate that the adsorption of large cations such as potassium or cesium facilitates the epitaxial growth of polar LaAlO₃ (LAO) on SrTiO₃ (STO). The low ionization potential of K favors efficient electron transfer to the STO conduction band and results in a 2D electron gas which exactly compensates for the diverging potential with increasing layer thickness. For large cations like K or Cs, DFT total energy considerations show that they remain adsorbed on the LAO surface and do not enter substitutionally into LAO. These results suggest a novel scheme for growing clean LAO/STO interface systems, and polar systems in general, by performing the growth process in the presence of large, low ionization potential alkali metal ions.

View record

Nonlocal Dielectric Response of Water (2012)

Water, possessing a permanent dipolar moment, generates a peculiar potential of mean force between solvated ions, the investigation of which is the goal of this thesis. Other authors have observed overscreening between two charges in water, leading to attractive regimes in otherwise repulsive potentials, and vice versa. In this thesis, atomistic molecular dynamics simulations of Extended Simple Point Charge (SPC/E) water are used to generate a nonlocal dielectric function that matches with other computational and experimental results from literature. Nonlocality means the consideration of the orientation of neighbouring water molecules when finding the displacement field caused by an electric field. Specifically, instead of the displacement field at a given point depending only on the electric field at that point, with nonlocality the displacement field involves an integral of the electric field over all space. To compare to this nonlocal dielectric function derived from simulations of only water, simulations of water with ions are performed. Simulations of two ions in water demonstrate the existence of nonlocality, also in agreement with current scientific understanding. Further, simulations of three ions show that the effective mean force in more complicated systems is not simply a summation of pair-wise forces but instead must account for multi-body effects.

View record