Publications

The Hund-metal route to strong correlations continues to attract large interest in the condensed-matter community. The question arose to what extent it applies to the infinite-layer nickelates and, as a related question, to two-orbital systems in general. Here, we provide a low-energy perspective on this topic through a dynamical mean-field study using the numerical renormalization group (NRG) as a real-frequency impurity solver. We find that the RG flow from high to low energy is a uniquely adequate tool to reveal two-stage Kondo screening (2SKS), a fascinating mechanism for Hund physics. Further, we show that 2SKS takes place in a quarter-filled two-orbital system, but can be easily suppressed by a sufficiently large crystal-field splitting. We apply these insights to LaNiO2 using a recently proposed two-orbital model and show that it is indeed the crystal-field splitting that suppresses multiorbital phenomena in this scenario. Our general findings open the way for further explorations of 2SKS, and we propose a way of reviving low-energy Hund physics in LaNiO2 by counteracting the crystal field.

Conceptually, the Matsubara formalism (MF), using imaginary frequencies, and the Keldysh formalism (KF), formulated in real frequencies, give equivalent results for systems in thermal equilibrium. The MF has less complexity and is thus more convenient than the KF. However, computing dynamical observables in the MF requires the analytic continuation from imaginary to real frequencies. The analytic continuation is well-known for two-point correlation functions (having one frequency argument), but, for multipoint correlators, a straightforward recipe for deducing all Keldysh components from the MF correlator had not been formulated yet. Recently, a representation of MF and KF correlators in terms of formalism-independent partial spectral functions and formalism-specific kernels was introduced by Kugler, Lee, and von Delft [Phys. Rev. X 11, 041006 (2021)]. We use this representation to formally elucidate the connection between both formalisms. We show how a multipoint MF correlator can be analytically continued to recover all partial spectral functions and yield all Keldysh components of its KF counterpart. The procedure is illustrated for various correlators of the Hubbard atom.

Effectively compressing and optimizing tensor networks requires reliable methods for fixing the latent degrees of freedom of the tensors, known as the gauge. Here we introduce a new algorithm for gauging tensor networks using belief propagation, a method that was originally formulated for performing statistical inference on graphical models and has recently found applications in tensor network algorithms. We show that this method is closely related to known tensor network gauging methods. It has the practical advantage, however, that existing belief propagation implementations can be repurposed for tensor network gauging, and that belief propagation is a very simple algorithm based on just tensor contractions so it can be easier to implement, optimize, and generalize. We present numerical evidence and scaling arguments that this algorithm is faster than existing gauging algorithms, demonstrating its usage on structured, unstructured, and infinite tensor networks. Additionally, we apply this method to improve the accuracy of the widely used simple update gate evolution algorithm.

We explore all the chiral limits of perturbed Dirac field theories applicable to C

The phenomena of superconductivity and charge density waves are observed in close vicinity in many strongly correlated materials. Increasing evidence from experiments and numerical simulations suggests both phenomena can also occur in an intertwined manner, where the superconducting order parameter is coupled to the electronic density. Employing density matrix renormalization group simulations, we investigate the nature of such an intertwined state of matter stabilized in the phase diagram of the elementary t-t-U Hubbard model in the strong coupling regime. Remarkably, the condensate of Cooper pairs is shown to be fragmented in the presence of a charge density wave where more than one pairing wave function is macroscopically occupied. Moreover, we provide conclusive evidence that the macroscopic wave functions of the superconducting fragments are well-described by soliton solutions of a Ginzburg-Landau equation in a periodic potential constituted by the charge density wave. In the presence of an orbital magnetic field, the order parameters are gauge invariant, and superconducting vortices are pinned between the stripes. This intertwined Ginzburg-Landau theory is proposed as an effective low-energy description of the stripe fragmented superconductor.

One of the demanding frontiers in ultracold science is identifying laser cooling schemes for complex atoms and molecules, out of their vast spectra of internal states. Motivated by a need to expand the set of available ultracold molecules for applications in fundamental physics, chemistry, astrochemistry, and quantum simulation, we propose and demonstrate an automated graph-based search approach for viable laser cooling schemes. The method is time efficient and the outcomes greatly surpass the results of manual searches used so far. We discover new laser cooling schemes for C

Synthetic dimension platforms offer unique pathways for engineering quantum matter. We compute the phase diagram of a many-body system of ultracold atoms (or polar molecules) with a set of Rydberg states (or rotational states) as a synthetic dimension, where the particles are arranged in real space in optical microtrap arrays and interact via dipole-dipole exchange interaction. Using mean-field theory, we find three ordered phases - two are localized in the synthetic dimension, predicted as zero-temperature ground states in Refs. [Sci. Rep., 8, 1 (2018) and Phys. Rev. A 99, 013624 (2019)], and a delocalized phase. We characterize them by identifying the spontaneously broken discrete symmetries of the Hamiltonian. We also compute the phase diagram as a function of temperature and interaction strength, for both signs of the interaction. For system sizes with more than six synthetic sites and attractive interactions, we find that the thermal phase transitions can be first or second order, which leads to a tri-critical point on the phase boundary. By examining the dependence of the tri-critical point and other special points of the phase boundary on the synthetic dimension size, we shed light on the physics for thermodynamically large synthetic dimension.

The performance of semiconductor devices can be enhanced by understanding and characterizing the bound states of deep defects. However, simulating these states using conventional

In conventional

Using scanning-tunneling-microscopy and theoretical modeling on heterostructures of twisted bilayer graphene and hexagonal Boron-Nitride, we show that the emergent super-moire structures display a rich landscape of moire-crystals and quasicrystals. We reveal a phase-diagram comprised of commensurate moire-crystals embedded in swaths of moire quasicrystals. The 1:1 commensurate crystal, expected to be a Chern insulator, should only exist at one point on the phase-diagram, implying that it ought to be practically undetectable. Surprisingly we find that the commensurate crystals exist over a much wider than predicted range, providing evidence of an unexpected self-alignment mechanism that is explained using an elastic-network model. The remainder of the phase-diagram, where we observe tunable quasicrystals, affords a new platform for exploring the unique electronic-properties of these rarely found in nature structures.