Publications

Interface quantum materials have yielded a plethora of previously unknown phenomena, including unconventional superconductivity, topological phases, and possible Majorana fermions. Typically, such states are detected at the interface between two insulating constituents by electrical transport, but whether either material is conducting, transport techniques become insensitive to interfacial properties. To overcome these limitations, we use angle-resolved photoemission spectroscopy and molecular beam epitaxy to reveal the electronic structure, charge transfer, doping profile, and carrier effective masses in a layer-by-layer fashion for the interface between the Dirac nodal-line semimetal SrIrO3 and the correlated metallic Weyl ferromagnet SrRuO3. We find that electrons are transferred from the SrIrO3 to SrRuO3, with an estimated screening length of λ = 3.2 ± 0.1 Å. In addition, we find that metallicity is preserved even down to a single SrIrO3 layer, where the dimensionality-driven metal-insulator transition typically observed in SrIrO3 is avoided because of strong hybridization of the Ir and Ru t2g states. Tomographic spectroscopy reveals how the properties of topological materials can be engineered at interfaces.

This article summarizes recent work on the many-body (beyond density functional theory) electronic structure of layered rare-earth nickelates, both in the context of the materials themselves and in comparison to the high-temperature superconducting (high-T

The symmetry of the superconducting order parameters in multi-orbital systems involves many quantum numbers. Pairing mediated by electronic correlations being retarded, the frequency structure of superconducting order parameters bears important information. Here we generalize the frequency-dependent theory of superconductivity mediated by spin and charge fluctuations to systems with spin-orbit coupling. This formulation is applied to strontium ruthenate with a normal state obtained using density functional theory. Taking advantage of pseudospin and inversion symmetries, the inter-pseudospin sector of the normal state Eliashberg equation is mapped to a pseudospin-diagonal one. We find ubiquitous entanglement of spin and orbital quantum numbers, along with notable mixing between even- and odd-frequency correlations. We present the phase diagrams for leading and subleading symmetries in the pseudospin-orbital basis of Sr

In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules. We show that the HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant error. The ability of HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems which require a larger active space than possible with complete active space self-consistent field (CASSCF). Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT). We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional. We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet-triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface.

We present a stable and systematically improvable quantum Monte Carlo (QMC) approach to calculating excited-state energies, which we implement using our fast randomized iteration method for the full configuration interaction problem (FCI-FRI). Unlike previous excited-state quantum Monte Carlo methods, our approach, which is an asymmetric variant of subspace iteration, avoids the use of dot products of random vectors and instead relies upon trial vectors to maintain orthogonality and estimate eigenvalues. By leveraging recent advances, we apply our method to calculate ground- and excited-state energies of strongly correlated molecular systems in large active spaces, including the carbon dimer with 8 electrons in 108 orbitals (8e,108o), an oxo-Mn(salen) transition metal complex (28e,28o), ozone (18e,87o), and butadiene (22e,82o). In the majority of these test cases, our approach yields total excited-state energies that agree with those from state-of-the-art methods -- including heat-bath CI, the density matrix renormalization group approach, and FCIQMC -- to within sub-milliHartree accuracy. In all cases, estimated excitation energies agree to within about 0.1 eV.

Auxiliary field Quantum Monte Carlo simulations of interacting fermions require sampling over a Hubbard-Stratonovich field h introduced to decouple the interactions. The weight for a given configuration involves the products of the determinant of matrices M

We perform a data-driven dimensionality reduction of the scale-dependent 4-point vertex function characterizing the functional Renormalization Group (fRG) flow for the widely studied two-dimensional t - t' Hubbard model on the square lattice. We demonstrate that a deep learning architecture based on a Neural Ordinary Differential Equation solver in a low-dimensional latent space efficiently learns the fRG dynamics that delineates the various magnetic and d-wave superconducting regimes of the Hubbard model. We further present a Dynamic Mode Decomposition analysis that confirms that a small number of modes are indeed sufficient to capture the fRG dynamics. Our work demonstrates the possibility of using artificial intelligence to extract compact representations of the 4-point vertex functions for correlated electrons, a goal of utmost importance for the success of cutting-edge quantum field theoretical methods for tackling the many-electron problem.

The optoelectronic and transport properties of two-dimensional transition metal dichalcogenide semiconductors (2D TMDs) are highly susceptible to external perturbation, enabling precise tailoring of material function through post-synthetic modifications. Here we show that nanoscale inhomogeneities known as nanobubbles can be used for both strain and, less invasively, dielectric tuning of exciton transport in bilayer tungsten disulfide (WSe2). We use ultrasensitive spatiotemporally resolved optical scattering microscopy to directly image exciton transport, revealing that dielectric nanobubbles are surprisingly efficient at funneling and trapping excitons at room temperature, even though the energies of the bright excitons are negligibly affected. Our observations suggest that exciton funneling in dielectric inhomogeneities is driven by momentum-indirect (dark) excitons whose energies are more sensitive to dielectric perturbations than bright excitons. These results reveal a new pathway to control exciton transport in 2D semiconductors with exceptional spatial and energetic precision using dielectric engineering of dark state energetic landscapes.

The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is O(N

Competing and intertwined orders including inhomogeneous patterns of spin and charge are observed in many correlated electron materials, such as high-temperature superconductors. Introducing a new development of the constrained-path auxiliary-field quantum Monte Carlo (AFQMC) method, we study the interplay between thermal and quantum fluctuations in the two-dimensional Hubbard model. We obtain an accurate and systematic characterization of the evolution of the spin and charge correlations as a function of temperature T and how it connects to the ground state, at three representative doping levels δ= 1/5, 1/8, and 1/10. We find increasing short-range commensurate antiferromagnetic correlations as T is lowered. As the correlation length grows sufficiently large, a modulated spin-density-wave (SDW) appears. At δ= 1/5, the SDW saturates and remains short-ranged as T →0. In contrast, at δ= 1/8 and 1/10 this evolves into a ground-state stripe phase. We study the relation between spin and charge orders and find that formation of charge order appears to be driven by that of the spin order. We identify a finite-temperature phase transition below which charge ordering sets in and discuss the implications of our results for the nature of this transition.