Publications

A method to calculate the one-body Green's function for ground states of correlated electron materials is formulated by extending the variational Monte Carlo method. We benchmark against the exact diagonalization (ED) for the one- and two-dimensional Hubbard models of 16 site lattices, which proves high accuracy of the method. The application of the method to larger-sized Hubbard model on the square lattice correctly reproduces the Mott insulating behavior at half filling and gap structures of d-wave superconducting state of the hole doped Hubbard model in the ground state optimized by enforcing the charge uniformity, evidencing a wide applicability to strongly correlated electron systems. From the obtained d-wave superconducting gap of the charge uniform state, we find that the gap amplitude at the antinodal point is several times larger than the experimental value, when we employ a realistic parameter as a model of the cuprate superconductors. The effective attractive interaction of carriers in the d-wave superconducting state inferred for an optimized state of the Hubbard model is as large as the order of the nearest-neighbor transfer, which is far beyond the former expectation in the cuprates. We discuss the nature of the superconducting state of the Hubbard model in terms of the overestimate of the gap and the attractive interaction in comparison to the cuprates.

Predicting the properties of a quantum-mechanical system of many interacting particles is a major goal of modern science and an outstanding challenge. We consider a compact and versatile model which captures the essential features of a broad class of systems made of particles obeying Bose–Einstein statistics and which allows one to systematically dial up the effect of quantum entanglement in the presence of particle interaction by tuning a single parameter. We are able to obtain exact numerical results for the phase diagrams of these systems. Possible directions for experimental realization of the predictions are discussed.We present a comprehensive theoretical study of the phase diagram of a system of many Bose particles interacting with a two-body central potential of the so-called Lennard-Jones form. First-principles path-integral computations are carried out, providing essentially exact numerical results on the thermodynamic properties. The theoretical model used here provides a realistic and remarkably general framework for describing simple Bose systems ranging from crystals to normal fluids to superfluids and gases. The interplay between particle interactions on the one hand and quantum indistinguishability and delocalization on the other hand is characterized by a single quantumness parameter, which can be tuned to engineer and explore different regimes. Taking advantage of the rare combination of the versatility of the many-body Hamiltonian and the possibility for exact computations, we systematically investigate the phases of the systems as a function of pressure (P) and temperature (T), as well as the quantumness parameter. We show how the topology of the phase diagram evolves from the known case of 4He, as the system is made more (and less) quantum, and compare our predictions with available results from mean-field theory. Possible realization and observation of the phases and physical regimes predicted here are discussed in various experimental systems, including hypothetical muonic matter.There are no data underlying this work.

Quantum state transformations that are robust to experimental imperfections are important for applications in quantum information science and quantum sensing. Counterdiabatic (CD) approaches, which use knowledge of the underlying system Hamiltonian to actively correct for diabatic effects, are powerful tools for achieving simultaneously fast and stable state transformations. Protocols for CD driving have thus far been limited in their experimental implementation to discrete systems with just two or three levels, as well as bulk systems with scaling symmetries. Here, we extend the tool of CD control to a discrete synthetic lattice system composed of as many as nine sites. Although this system has a vanishing gap and thus no adiabatic support in the thermodynamic limit, we show that CD approaches can still give a substantial, several order-of-magnitude, improvement in fidelity over naive, fast adiabatic protocols.

A stochastic procedure is developed which allows one to express Pontryagin's maximum principle for dissipative quantum system solely in terms of stochastic wave functions. Time-optimal controls can be efficiently computed without computing the density matrix. Specifically, the proper dynamical update rules are presented for the stochastic costate variables introduced by Pontryagin's maximum principle and restrictions on the form of the terminal cost function are discussed. The proposed procedure is confirmed by comparing the results to those obtained from optimal control on Lindbladian dynamics. Numerically, the proposed formalism becomes time and memory efficient for large systems, and it can be generalized to describe non-Markovian dynamics.

Work function-mediated charge transfer in graphene/α-RuCl

We present a technique that enables the evaluation of perturbative expansions based on one-loop-renormalized vertices up to large expansion orders. Specifically, we show how to compute large-order corrections to the random phase approximation in either the particle-hole or particle-particle channels. The algorithm's efficiency is achieved by the summation over contributions of all symmetrized Feynman diagram topologies using determinants, and by integrating out analytically the two-body long-range interactions in order to yield an effective zero-range interaction. Notably, the exponential scaling of the algorithm as a function of perturbation order leads to a polynomial scaling of the approximation error with computational time for a convergent series. To assess the performance of our approach, we apply it to the non-perturbative regime of the square-lattice fermionic Hubbard model away from half-filling and report, as compared to the bare interaction expansion algorithm, significant improvements of the Monte Carlo variance as well as the convergence properties of the resulting perturbative series.

We study disorder induced topological phase transitions in magnetically doped (Bi, Sb)

In this work we introduce the Dual Boson Diagrammatic Monte Carlo technique for strongly interacting electronic systems. This method combines the strength of dynamical mean-filed theory for non-perturbative description of local correlations with the systematic account of non-local corrections in the Dual Boson theory by the diagrammatic Monte Carlo approach. It allows us to get a numerically exact solution of the dual boson theory at the two-particle local vertex level for the extended Hubbard model. We show that it can be efficiently applied to description of single particle observables in a wide range of interaction strengths. We compare our exact results for the self-energy with the ladder Dual Boson approach and determine a physical regime, where description of collective electronic effects requires more accurate consideration beyond the ladder approximation. Additionally, we find that the order-by-order analysis of the perturbative diagrammatic series for the single-particle Green's function allows to estimate the transition point to the charge density wave phase.

We present the theoretical foundations and the implementation details of a density-functional approach for coupled photons, electrons, and effective nuclei in non-relativistic quantum electrodynamics. Starting point of the formalism is a generalization of the Pauli-Fierz field theory for which we establish a one-to-one correspondence between external fields and internal variables. Based on this correspondence, we introduce a Kohn-Sham construction which provides a computationally feasible approach for ab-initio light-matter interactions. In the mean-field limit for the effective nuclei the formalism reduces to coupled Ehrenfest-Maxwell-Pauli-Kohn-Sham equations. We present an implementation of the approach in the real-space real-time code Octopus. For the implementation we use the Riemann-Silberstein formulation of classical electrodynamics and rewrite Maxwell's equations in Schrödinger form. This allows us to use existing time-evolution algorithms developed for quantum-mechanical systems also for Maxwell's equations. Overall, our approach is ideally suited for applications in nano-optics, nano-plasmonics, (photo) electrocatalysis, light-matter coupling in 2D materials, cases where laser pulses carry orbital angular momentum, or light-tailored chemical reactions in optical cavities to name but a few.

Materials with strong electronic correlations host remarkable -- and technologically relevant -- phenomena such as magnetism, superconductivity and metal-insulator transitions. Harnessing and controlling these effects is a major challenge, on which key advances are being made through lattice and strain engineering in thin films and heterostructures, leveraging the complex interplay between electronic and structural degrees of freedom. Here we show that the electronic structure of LaNiO3 can be tuned by means of lattice engineering. We use different substrates to induce compressive and tensile biaxial epitaxial strain in LaNiO3 thin films. Our measurements reveal systematic changes of the optical spectrum as a function of strain and, notably, an increase of the low-frequency free carrier weight as tensile strain is applied. Using density functional theory (DFT) calculations, we show that this apparently counter-intuitive effect is due to a change of orientation of the oxygen octahedra.The calculations also reveal drastic changes of the electronic structure under strain, associated with a Fermi surface Lifshitz transition. We provide an online applet to explore these effects. The experimental value of integrated spectral weight below 2 eV is significantly (up to a factor of 3) smaller than the DFT results, indicating a transfer of spectral weight from the infrared to energies above 2 eV. The suppression of the free carrier weight and the transfer of spectral weight to high energies together indicate a correlation-induced band narrowing and free carrier mass enhancement due to electronic correlations. Our findings provide a promising avenue for the tuning and control of quantum materials employing lattice engineering.