Motivated by the recent claim of hot superconductivity with critical temperatures up to 550 K in La + x hydrides (arXiv:2006.03004), we investigate the high-pressure phase diagram of possible compounds that may have formed in the experiment, using fi
rst-principles calculations for evolutionary crystal structure prediction and superconductivity. Starting from the hypothesis that the observed T$_c$ may be realized by successive heating upon a pre-formed LaH$_{10}$ phase, we examine plausible ternaries of lanthanum, hydrogen and other elements present in the diamond anvil cell: boron, nitrogen, carbon, platinum, gallium, gold. We find that only boron forms superhydride-like structures that can host high-Tc superconductivity, but the predicted Tc are incompatibe with the experimental reports. Our results indicate that, while the claims of hot superconductivity should be reconsidered, it is very likely that unkwown H-rich ternary or multinary phases containing lanthanum, boron and hydrogen may have formed under the reported experimental conditions, and that these may exhibit superconducting properties comparable, or even superior, to those of currently known hydrides.
In the last five years a large number of new high-temperature superconductors have been predicted and experimentally discovered among hydrogen-rich crystals, at pressures which are way too high to meet any practical application. In this work, we repo
rt the computational prediction of a hydride superconductor, LaBH$_{8}$, with a T$_{c}$ of 126 K at a pressure of 50 GPa, thermodynamically stable above 100 GPa, and dynamically stable down to 40 GPa, an unprecedentedly low pressure for high-T$_{c}$ hydrides. LaBH$_{8}$ can be seen as a ternary sodalite-like hydride, in which a metallic hydrogen sublattice is stabilized by the chemical pressure exerted by the guest elements. The combination of two elements with different atomic sizes in LaBH$_{8}$ realizes a more efficient packing of atoms than in binary sodalite hydrides. A suitable choice of elements may be exploited to further reduce the stabilization pressure to ambient conditions.
Motivated by the recent discovery of near-room temperature superconductivity in high-pressure superhydrides, we investigate from first-principles the high-pressure superconducting phase diagram of the ternary Ca-B-H system, using ab-initio evolutiona
ry crystal structure prediction, and Density Functional Perturbation Theory. We find that below 100 GPa all stable and weakly metastable phases are insulating. This pressure marks the appearance of several new chemically-forbidden phases on the hull of stability, and the first onset of metalization in CaBH$_5$. Metallization is then gradually achieved at higher pressure at different compositions. Among the metallic phases stable in the Megabar regime, we predict two high-$T_c$ superconducting phases with CaBH$_6$ and Ca$_2$B$_2$H$_{13}$ compositions, with critical temperatures of 119 and 89 K at 300 GPa, respectively, surviving to lower pressures. Ternary hydrides will most likely play a major role in superconductivity research in the coming years; our study suggests that, in order to reduce the pressure for the onset of metallicity and superconductivity, further explorations of ternary hydrides should focus on elements less electronegative than boron.
Bayesian parametric analytic continuation (BPAC) is proposed for the analytic continuation of noisy imaginary-time Greens function data as, e.g., obtained by continuous-time quantum Monte Carlo simulations (CTQMC). Within BPAC, the spectral function
is inferred from a suitable set of parametrized basis functions. Bayesian model comparison then allows to assess the reliability of different parametrizations. The required evidence integrals of such a model comparison are determined by nested sampling. Compared to the maximum entropy method (MEM), routinely used for the analytic continuation of CTQMC data, the presented approach allows to infer whether the data support specific structures of the spectral function. We demonstrate the capability of BPAC in terms of CTQMC data for an Anderson impurity model (AIM) that shows a generalized Kondo scenario and compare the BPAC reconstruction to the MEM as well as to the spectral function obtained from the real-time fork tensor product state impurity solver where no analytic continuation is required. Furthermore, we present a combination of MEM and BPAC and its application to an AIM arising from the ab initio treatment of SrVO$_3$.
We present a general scheme to map correlated nonequilibrium quantum impurity problems onto an auxiliary open quantum system of small size. The infinite fermionic reservoirs of the original system are thereby replaced by a small number $N_B$ of nonin
teracting auxiliary bath sites whose dynamics is described by a Lindblad equation. Due to the presence of the intermediate bath sites, the overall dynamics acting on the impurity site is non-Markovian. With the help of an optimization scheme for the auxiliary Lindblad parameters, an accurate mapping is achieved, which becomes exponentially exact upon increasing $N_B$. The basic idea for this scheme was presented previously in the context of nonequilibrium dynamical mean field theory. In successive works on improved manybody solution strategies for the auxiliary Lindblad equation, such as Lanczos exact diagonalization or matrix product states, we applied the approach to study the nonequilibrium Kondo regime. In the present paper, we address in detail the mapping procedure itself, rather than the many-body solution. In particular, we investigate the effects of the geometry of the auxiliary system on the accuracy of the mapping for given $N_B$. Specifically, we present a detailed convergence study for five different geometries which, besides being of practical utility, reveals important insights into the underlying mechanisms of the mapping. For setups with onsite or nearest-neighbor Lindblad parameters we find that a representation adopting two separate bath chains is by far more accurate with respect to other choices based on a single chain or a commonly used star geometry. A significant improvement is obtained by allowing for long-ranged and complex Lindblad parameters. These results can be of great value when studying Lindblad-type approaches to correlated systems.
A numerical approach is presented that allows to compute nonequilibrium steady state properties of strongly correlated quantum many-body systems. The method is imbedded in the Keldysh Greens function formalism and is based upon the idea of the variat
ional cluster approach as far as the treatment of strong correlations is concerned. It appears that the variational aspect is crucial as it allows for a suitable optimization of a reference system to the nonequilibrium target state. The approach is neither perturbative in the many-body interaction nor in the field, that drives the system out of equilibrium, and it allows to study strong perturbations and nonlinear responses of systems in which also the correlated region is spatially extended. We apply the presented approach to non-linear transport across a strongly correlated quantum wire described by the fermionic Hubbard model. We illustrate how the method bridges to cluster dynamical mean-field theory upon coupling two baths containing and increasing number of uncorrelated sites.
Among the various numerical techniques to study the physics of strongly correlated quantum many-body systems, the self-energy functional approach (SFA) has become increasingly important. In its previous form, however, SFA is not applicable to Bose-Ei
nstein condensation or superfluidity. In this paper we show how to overcome this shortcoming. To this end we identify an appropriate quantity, which we term $D$, that represents the correlation correction of the condensate order parameter, as it does the self-energy for the Greens function. An appropriate functional is derived, which is stationary at the exact physical realizations of $D$ and of the self-energy. Its derivation is based on a functional-integral representation of the grand potential followed by an appropriate sequence of Legendre transformations. The approach is not perturbative and therefore applicable to a wide range of models with local interactions. We show that the variational cluster approach based on the extended self-energy functional is equivalent to the pseudoparticle approach introduced in Phys. Rev. B, 83, 134507 (2011). We present results for the superfluid density in the two-dimensional Bose-Hubbard model, which show a remarkable agreement with those of Quantum-Monte-Carlo calculations.
