Related papers: Increased Performance of Matsubara space calculati…
We present an alternative representation for the anisotropic Eliashberg equations of superconductivity, whose numerical solution yields an efficiency gain of several orders of magnitude with respect to the conventional representation in…
Experimentally relevant signatures of superconductivity require access to real-frequency quantities, such as the spectral functions, optical response, and transport properties, yet Migdal-Eliashberg calculations are commonly performed on…
We present a compressive sensing approach for the long standing problem of Matsubara summation in many-body perturbation theory. By constructing low-dimensional, almost isometric subspaces of the Hilbert space we obtain optimum imaginary…
Migdal-Eliashberg theory is one of the state-of-the-art methods for describing conventional superconductors from first principles. However, widely used implementations assume a constant density of states around the Fermi level, which…
We solve the anisotropic, full-bandwidth and non-adiabatic Eliashberg equations for phonon-mediated superconductivity by fully including the first vertex correction in the electronic self-energy. The non-adiabatic equations are solved…
This paper introduces the Julia package IsoME, an easy-to-use yet accurate and robust computational tool designed to calculate superconducting properties. Multiple levels of approximation are supported, ranging from the basic…
We present numerical results showing how our recently proposed relativistic three-particle quantization condition can be used in practice. Using the isotropic (generalized $s$-wave) approximation, and keeping only the leading terms in the…
This article begins with a brief introduction to numerical relativity aimed at readers who have a background in applied mathematics but not necessarily in general relativity. I then introduce and summarise my work on the problem of treating…
We introduce a generalized finite difference method for solving a large range of fully nonlinear elliptic partial differential equations in three dimensions. Methods are based on Cartesian grids, augmented by additional points carefully…
In this paper, we introduce a new scheme for the efficient numerical treatment of the electronic Schr\"odinger equation for molecules. It is based on the combination of a many-body expansion, which corresponds to the so-called bond order…
This study investigates numerical methods to solve nonlinear transport problems characterized by various sorption isotherms with a focus on the Freundlich type of isotherms. We describe and compare second order accurate numerical schemes,…
We introduce a modification of the Fast Marching Algorithm, which solves the generalized eikonal equation associated to an arbitrary continuous riemannian metric, on a two or three dimensional domain. The algorithm has a logarithmic…
The effectiveness of the hyperbolic relaxation method for solving the Einstein constraint equations numerically is studied here on a variety of compact orientable three-manifolds. Convergent numerical solutions are found using this method…
In order to reduce the computational cost of the simulation of electromagnetic responses in geophysical settings that involve highly heterogeneous media, we develop a multiscale finite volume method with oversampling for the quasi-static…
In this work numerical methods for solving Einstein's equations are developed and applied to the study of inhomogeneous cosmological models. A two-dimensional computer code is described which implements two advanced numerical methods:…
We present different methods to increase the performance of Hybrid Monte Carlo simulations of the Hubbard model in two-dimensions. Our simulations concentrate on a hexagonal lattice, though can be easily generalized to other lattices. It is…
We investigate the possibility to assist the numerically ill-posed calculation of spectral properties of interacting quantum systems in thermal equilibrium by extending the imaginary-time simulation to a finite Schwinger-Keldysh contour.…
This paper develops a strong computational approach to simulate a three-dimensional nonlinear radiation-conduction model in optically thick media, subject to suitable initial and boundary conditions. The space derivatives are approximated…
Application of integral equation theory to complex fluids is reviewed, with particular emphasis to the effects of polydispersity and anisotropy on their structural and thermodynamic properties. Both analytical and numerical solutions of…
In this paper we present a methodology for increasing the accuracy and accelerating the convergence of numerical methods for solution of Maxwell's equations in the frequency domain by taking into account the be-havior of the electromagnetic…