Efficient Algorithms for Electrostatic Interactions Including Dielectric Contrasts

A simple efficient algorithm for molecular simulations of constant potential electrodes

Increasingly, society requires high power, high energy storage devices for applications ranging from electric vehicles to buffers on the electric grid. Supercapacitors are a promising contribution to meeting these demands, though there still remain unsolved practical problems. Molecular dynamics simulations can shed light on the relevant molecular level processes in electric double layer capacitors, but these simulations are computationally very demanding. Our focus here is on the algorithmic complexity of the constant potential method (CPM), which uses dedicated electrostatics solvers to maintain a fixed potential difference between two conducting electrodes. We show how any standard electrostatics solver-capable of calculating the energies and forces on all atoms-can be used to implement CPM with a minimum of coding. As an example, we compare our generalized implementation of CPM, based on invocations of the particle-particle-particle-mesh routine of the Large-scale Atomic/Molecular Massively Parallel Simulator, with a traditional implementation based on a dedicated re-implementation of Ewald summation. Both methods yield comparable results on four test systems, with the former achieving a substantial gain in speed and improved scalability. The step from dedicated electrostatic solvers to generic routines is made possible by noting that CPM's traditional narrow Gaussian point-spread of atomic charges on the electrodes effectively endows point-like atoms with chemical hardness, i.e., an intra-atomic energy quadratic in the charge.

Theory and Simulations of Ionic Liquids in Nanoconfinement

Room-temperature ionic liquids (RTILs) have exciting properties such as nonvolatility, large electrochemical windows, and remarkable variety, drawing much interest in energy storage, gating, electrocatalysis, tunable lubrication, and other applications. Confined RTILs appear in various situations, for instance, in pores of nanostructured electrodes of supercapacitors and batteries, as such electrodes increase the contact area with RTILs and enhance the total capacitance and stored energy, between crossed cylinders in surface force balance experiments, between a tip and a sample in atomic force microscopy, and between sliding surfaces in tribology experiments, where RTILs act as lubricants. The properties and functioning of RTILs in confinement, especially nanoconfinement, result in fascinating structural and dynamic phenomena, including layering, overscreening and crowding, nanoscale capillary freezing, quantized and electrotunable friction, and superionic state. This review offers a comprehensive analysis of the fundamental physical phenomena controlling the properties of such systems and the current state-of-the-art theoretical and simulation approaches developed for their description. We discuss these approaches sequentially by increasing atomistic complexity, paying particular attention to new physical phenomena emerging in nanoscale confinement. This review covers theoretical models, most of which are based on mapping the problems on pertinent statistical mechanics models with exact analytical solutions, allowing systematic analysis and new physical insights to develop more easily. We also describe a classical density functional theory, which offers a reliable and computationally inexpensive tool to account for some microscopic details and correlations that simplified models often fail to consider. Molecular simulations play a vital role in studying confined ionic liquids, enabling deep microscopic insights otherwise unavailable to researchers. We describe the basics of various simulation approaches and discuss their challenges and applicability to specific problems, focusing on RTIL structure in cylindrical and slit confinement and how it relates to friction and capacitive and dynamic properties of confined ions.

Frequency and field-dependent response of confined electrolytes from Brownian dynamics simulations

Using Brownian dynamics simulations, we investigate the effects of confinement, adsorption on surfaces, and ion-ion interactions on the response of confined electrolyte solutions to oscillating electric fields in the direction perpendicular to the confining walls. Nonequilibrium simulations allows to characterize the transitions between linear and nonlinear regimes when varying the magnitude and frequency of the applied field, but the linear response, characterized by the frequency-dependent conductivity, is more efficiently predicted from the equilibrium current fluctuations. To that end, we (rederive and) use the Green-Kubo relation appropriate for overdamped dynamics, which differs from the standard one for Newtonian or underdamped Langevin dynamics. This expression highlights the contributions of the underlying Brownian fluctuations and of the interactions of the particles between them and with external potentials. Although already known in the literature, this relation has rarely been used to date, beyond the static limit to determine the effective diffusion coefficient or the DC conductivity. The frequency-dependent conductivity always decays from a bulk-like behavior at high frequency to a vanishing conductivity at low frequency due to the confinement of the charge carriers by the walls. We discuss the characteristic features of the crossover between the two regimes, most importantly how the crossover frequency depends on the confining distance and the salt concentration, and the fact that adsorption on the walls may lead to significant changes both at high and low frequencies. Conversely, our results illustrate the possibility to obtain information on diffusion between walls, charge relaxation, and adsorption by analyzing the frequency-dependent conductivity.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

Electronic screening using a virtual Thomas-Fermi fluid for predicting wetting and phase transitions of ionic liquids at metal surfaces

Of relevance to energy storage, electrochemistry and catalysis, ionic and dipolar liquids display unexpected behaviours-especially in confinement. Beyond adsorption, over-screening and crowding effects, experiments have highlighted novel phenomena, such as unconventional screening and the impact of the electronic nature-metallic versus insulating-of the confining surface. Such behaviours, which challenge existing frameworks, highlight the need for tools to fully embrace the properties of confined liquids. Here we introduce a novel approach that involves electronic screening while capturing molecular aspects of interfacial fluids. Although available strategies consider perfect metal or insulator surfaces, we build on the Thomas-Fermi formalism to develop an effective approach that deals with any imperfect metal between these asymptotes. Our approach describes electrostatic interactions within the metal through a 'virtual' Thomas-Fermi fluid of charged particles, whose Debye length sets the screening length λ. We show that this method captures the electrostatic interaction decay and electrochemical behaviour on varying λ. By applying this strategy to an ionic liquid, we unveil a wetting transition on switching from insulating to metallic conditions.

Microscopic origin of the effect of substrate metallicity on interfacial free energies

We investigate the effect of the metallic character of solid substrates on solid-liquid interfacial thermodynamics using molecular simulations. Building on the recent development of a semiclassical Thomas-Fermi model to tune the metallicity in classical molecular dynamics simulations, we introduce a thermodynamic integration framework to compute the evolution of the interfacial free energy as a function of the Thomas-Fermi screening length. We validate this approach against analytical results for empty capacitors and by comparing the predictions in the presence of an electrolyte with values determined from the contact angle of droplets on the surface. The general expression derived in this work highlights the role of the charge distribution within the metal. We further propose a simple model to interpret the evolution of the interfacial free energy with voltage and Thomas-Fermi length, which allows us to identify the charge correlations within the metal as the microscopic origin of the evolution of the interfacial free energy with the metallic character of the substrate. This methodology opens the door to the molecular-scale study of the effect of the metallic character of the substrate on confinement-induced transitions in ionic systems, as reported in recent atomic force microscopy and surface force apparatus experiments.

Conformation and Ionization Behavior of Charge-Regulating Polyelectrolyte Brushes in a Poor Solvent

Understanding the structural response of weak polyelectrolyte brushes upon external stimuli is crucial for their applications ranging from modifying surface properties to the development of smart and intelligent materials. In this work, coarse-grained molecular dynamics simulations were carried out to investigate the conformation and ionization behavior of charge-regulating polyelectrolyte brushes under poor solvent conditions, using an implicit solvent model. The results show that, while the thickness of a sparse polyelectrolyte brush shows a similar behavior to that of a single chain, namely, a monotonic change as a function of solvent quality (modeled by an effective segment-segment attraction strength parameter) and solution pH, a dense polyelectrolyte brush exhibits more complex behavior. An unexpected reexpansion is observed when the effective segment-segment attraction strength is further increased, especially in the case of a high pH. In the latter case, strong attraction in polymer segments promotes the formation of large, interchain, cylindrical aggregates, leading to an increase in brush thickness.

Pimples reduce and dimples enhance flat dielectric surface image repulsion

In solid-liquid, or liquid-liquid, interfaces with dielectric contrast, charged particles interact with the induced polarization charge of the interface. These interactions contribute to an effective self-energy of the bulk ions and mediate ion-ion interactions. For flat interfaces, the self-energy and the mediated interactions are neatly constructed by the image charge method. For other geometries, explicit results are scarce and the problem must be treated via approximations or direct computation. The case of interfaces with roughness is of great practical importance. This article provides analytical results, valid to first-order in perturbation theory, for the self-energy of particles near rough substrates. Explicit formulas are provided for the case of a sinusoidal deformation of a flat surface. Generic deformations can be treated by superposition. In addition to results for the self-energy, the surface polarization charge is presented as a quadrature. The interaction between an ion and the deformed surface is modified by the change in relative distance as well as by the local curvature of the surface. Solid walls, with a lower dielectric constant than the liquid, repel all ions. We show that the repulsion is reduced by local convexity and enhanced by concavity; dimples are more repulsive than pimples.

A fast spectral method for electrostatics in doubly periodic slit channels

We develop a fast method for computing the electrostatic energy and forces for a collection of charges in doubly periodic slabs with jumps in the dielectric permittivity at the slab boundaries. Our method achieves spectral accuracy by using Ewald splitting to replace the original Poisson equation for nearly singular sources with a smooth far-field Poisson equation, combined with a localized near-field correction. Unlike existing spectral Ewald methods, which make use of the Fourier transform in the aperiodic direction, we recast the problem as a two-point boundary value problem in the aperiodic direction for each transverse Fourier mode for which exact analytic boundary conditions are available. We solve each of these boundary value problems using a fast, well-conditioned Chebyshev method. In the presence of dielectric jumps, combining Ewald splitting with the classical method of images results in smoothed charge distributions, which overlap the dielectric boundaries themselves. We show how to preserve the spectral accuracy in this case through the use of a harmonic correction, which involves solving a simple Laplace equation with smooth boundary data. We implement our method on graphical processing units and combine our doubly periodic Poisson solver with Brownian dynamics to study the equilibrium structure of double layers in binary electrolytes confined by dielectric boundaries. Consistent with prior studies, we find strong charge depletion near the interfaces due to repulsive interactions with image charges, which points to the need for incorporating polarization effects in understanding confined electrolytes, both theoretically and computationally.

Molecular Simulation of Electrode-Solution Interfaces

Many key industrial processes, from electricity production, conversion, and storage to electrocatalysis or electrochemistry in general, rely on physical mechanisms occurring at the interface between a metallic electrode and an electrolyte solution, summarized by the concept of an electric double layer, with the accumulation/depletion of electrons on the metal side and of ions on the liquid side. While electrostatic interactions play an essential role in the structure, thermodynamics, dynamics, and reactivity of electrode-electrolyte interfaces, these properties also crucially depend on the nature of the ions and solvent, as well as that of the metal itself. Such interfaces pose many challenges for modeling because they are a place where quantum chemistry meets statistical physics. In the present review, we explore the recent advances in the description and understanding of electrode-electrolyte interfaces with classical molecular simulations, with a focus on planar interfaces and solvent-based liquids, from pure solvent to water-in-salt electrolytes.

How to speed up ion transport in nanopores

Electrolyte-filled subnanometre pores exhibit exciting physics and play an increasingly important role in science and technology. In supercapacitors, for instance, ultranarrow pores provide excellent capacitive characteristics. However, ions experience difficulties in entering and leaving such pores, which slows down charging and discharging processes. In an earlier work we showed for a simple model that a slow voltage sweep charges ultranarrow pores quicker than an abrupt voltage step. A slowly applied voltage avoids ionic clogging and co-ion trapping—a problem known to occur when the applied potential is varied too quickly—causing sluggish dynamics. Herein, we verify this finding experimentally. Guided by theoretical considerations, we also develop a non-linear voltage sweep and demonstrate, with molecular dynamics simulations, that it can charge a nanopore even faster than the corresponding optimized linear sweep. For discharging we find, with simulations and in experiments, that if we reverse the applied potential and then sweep it to zero, the pores lose their charge much quicker than they do for a short-circuited discharge over their internal resistance. Our findings open up opportunities to greatly accelerate charging and discharging of subnanometre pores without compromising the capacitive characteristics, improving their importance for energy storage, capacitive deionization, and electrochemical heat harvesting.

Narrowing pores filled with an electrolyte usually slows down their charge-discharge dynamics. Here the authors demonstrate through molecular dynamics simulations and experiments with novolac-derived carbon electrodes how non-linear voltage sweeps can accelerate charging and discharging of subnanometer pores.

Translocation through a narrow pore under a pulling force

We employ a three-dimensional molecular dynamics to simulate a driven polymer translocation through a nanopore by applying an external force, for four pore diameters and two external forces. To see the polymer and pore interaction effects on translocation time, we studied nine interaction energies. Moreover, to better understand the simulation results, we investigate polymer center of mass, shape factor and the monomer spatial distribution through the translocation process. Our results reveal that increasing the polymer-pore interaction energy is accompanied by an increase in the translocation time and decrease in the process rate. Furthermore, for pores with greater diameter, the translocation becomes faster. The shape analysis of the polymer indicates that the polymer shape is highly sensitive to the interaction energy. In great interactions, the monomers come close to the pore from both sides. As a result, the translocation becomes fast at first and slows down at last. Overall, it can be concluded that the external force does not play a major role in the shape and distribution of translocated monomers. However, the interaction energy between monomer and nanopore has a major effect especially on the distribution of translocated monomers on the trans side.

Molecular Electron Density Distribution and X-Ray Diffraction Patterns of Smectic A Liquid Crystals - A Simulation Study

X-ray diffraction (XRD) is one of the most important methods to assess the long-range translational order in smectic A (SmA) liquid crystals. Nevertheless, the knowledge about the influence of the molecular electron density distribution (MEDD) on the XRD pattern is rather limited because it is not possible to vary the orientational order, the translational order and the MEDD independently in an experiment. We here present a systematic simulation study in which we examine this effect and show that the MEDD indeed has a major impact on the general appearance of the XRD pattern. More specifically, we find that the smectic layer peaks and the intensity ratios thereof strongly depend on the width of the MEDD. The classic approach by Leadbetter et al. to determine the smectic translational order parameter ∑ from XRD intensities works if the MEDD is quite narrow. In all other cases the influence of the MEDD has to be taken into account.

Feeling Your Neighbors across the Walls: How Interpore Ionic Interactions Affect Capacitive Energy Storage

Progress in low-dimensional carbon materials has intensified research on supercapacitors with nanostructured/nanoporous electrodes. The theoretical and simulation work so far has focused on charging single nanopores or nanoporous networks and the effects due to ionic interactions inside the pores, while the effect of interpore ion-ion correlations has received less attention. Herein, we study how the interactions between the ions in the neighboring pores across the pore walls affect capacitive energy storage. We develop a simple lattice model for the ions in a stack of parallel-aligned nanotubes, solve it by using the perturbation and "semi-mean-field" theories, and test the results by Monte Carlo simulations. We demonstrate that the interpore ionic interactions can have a profound effect on charge storage; in particular, such interactions can enhance or diminish the stored energy density, depending on the sign of like-charge interactions. We also find that charging can proceed either continuously or via a phase transition. Our results call for more detailed investigations of the properties of carbon pore walls and suggest that tuning their electrostatic response may be promising for the rational design of an optimal supercapacitor.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Structure and stimuli-responsiveness of all-DNA dendrimers: theory and experiment

We present a comprehensive theoretical and experimental study of the solution phase properties of a DNA-based family of nanoparticles - dendrimer-like DNA molecules (DL-DNA). These charged DNA dendrimers are novel macromolecular aggregates, which hold high promise in targeted self-assembly of soft matter systems in the bulk and at interfaces. To describe the behaviour of this family of dendrimers (with generations ranging from G1 to G7), we use a theoretical model in which base-pairs of a single DL-DNA molecule are modeled by charged monomers, whose interactions are chosen to mimic the equilibrium properties of DNA correctly. Experimental results on the sizes and conformations of DL-DNA are based on static and dynamic light scattering; and molecular dynamics simulations are employed to model the equilibrium properties of DL-DNA, which compare favorably to the findings from experiments while at the same time providing a host of additional information and insight into the molecular structure of the nanostructures. We also examine the salt-responsiveness of these macromolecules, finding that despite the strong screening of electrostatic interactions brought about by the added salt, the macromolecules shrink only slightly, their size robustness stemming from the high bending rigidity of the DNA-segments. The study of these charged dendrimer systems is an important field of research in the area of soft matter due to their potential role for various interdisciplinary applications, ranging from molecular cages and carriers for drug delivery in a living organism to the development of dendrimer- and dendron-based ultra-thin films in the area of nanotechnology. These findings are essential to determine if DL-DNA is a viable candidate for the experimental realization of cluster crystals in the bulk, a novel form of solid with multiple site occupancy.

Charge Me Slowly, I Am in a Hurry: Optimizing Charge-Discharge Cycles in Nanoporous Supercapacitors

Nanoporous supercapacitors attract much attention as green energy storage devices with remarkable cyclability and high power and energy densities. However, their use in high-frequency applications is limited by relatively slow charging processes, while accelerating charging without compromising the energy storage still remains a challenging task. Here, we study in detail the charging and discharging behavior of nanoporous supercapacitors with narrow pores, which provide exceptionally high capacitances and stored energy densities. We scrutinize the dynamic modes of charging, revealing, in particular, a transient formation of crowded and dilute ionic-liquid phases inside the pores, which leads to co-ion trapping and correspondingly slow charging. We show how trapping can be circumvented by applying a slow voltage sweep, and we demonstrate that it can accelerate the overall charging process considerably if the sweep rate is chosen appropriately. While one might be tempted to apply a similar strategy to discharging, we find that the best discharge rates are obtained when the voltage is switched off in a step-like fashion, whereby the optimal charge and discharge times differ a few-fold. We unveil the scaling laws for such optimal quantities, which allow one to predict quantitatively the charging behavior for realistically long pores. On the basis of our findings, we propose an optimal charge-discharge cycle and elaborate on optimization strategies.

Lattice model of ionic liquid confined by metal electrodes

We study, using Monte Carlo simulations, the density profiles and differential capacitance of ionic liquids confined by metal electrodes. To compute the electrostatic energy, we use the recently developed approach based on periodic Green's functions. The method also allows us to easily calculate the induced charge on the electrodes permitting an efficient implementation of simulations in a constant electrostatic potential ensemble. To speed up the simulations further, we model the ionic liquid as a lattice Coulomb gas and precalculate the interaction potential between the ions. We show that the lattice model captures the transition between camel-shaped and bell-shaped capacitance curves-the latter characteristic of ionic liquids (strong coupling limit) and the former of electrolytes (weak coupling). We observe the appearance of a second peak in the differential capacitance at ≈0.5 V for 2:1 ionic liquids, as the packing fraction is increased. Finally, we show that ionic size asymmetry decreases substantially the capacitance maximum, when all other parameters are kept fixed.

Charged nanorods at heterogeneously charged surfaces

We study the spatial and orientational distribution of charged nanorods (rodlike counterions) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The nanorods are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude.

The effect of finite pore length on ion structure and charging

Nanoporous supercapacitors play an important role in modern energy storage systems, and their modeling is essential to predict and optimize the charging behaviour. Two classes of models have been developed that consist of finite and infinitely long pores. Here, we show that although both types of models predict qualitatively consistent results, there are important differences emerging due to the finite pore length. In particular, we find that the ion density inside a finite pore is not constant, but increases linearly from the pore entrance to the pore end, where the ions form a strongly layered structure. This hinders a direct quantitative comparison between the two models. In addition, we show that although the ion density between the electrodes changes appreciably with the applied potential, this change has a minor effect on charging. Our simulations also reveal a complex charging behaviour, which is adsorption-driven at high voltages, but it is dominated either by co-ion desorption or by adsorption of both types of ions at low voltages, depending on the ion concentration.

Simulations of Coulomb systems confined by polarizable surfaces using periodic Green functions

We present an efficient approach for simulating Coulomb systems confined by planar polarizable surfaces. The method is based on the solution of the Poisson equation using periodic Green functions. It is shown that the electrostatic energy arising from the surface polarization can be decoupled from the energy due to the direct Coulomb interaction between the ions. This allows us to combine an efficient Ewald summation method, or any other fast method for summing over the replicas, with the polarization contribution calculated using Green function techniques. We apply the method to calculate density profiles of ions confined between the charged dielectric and metal surfaces.

Computer simulations of single particles in external electric fields

Applying electric fields is an attractive way to control and manipulate single particles or molecules, e.g., in lab-on-a-chip devices. However, the response of nanosize objects in electrolyte solution to external fields is far from trivial. It is the result of a variety of dynamical processes taking place in the ion cloud surrounding charged particles and in the bulk electrolyte, and it is governed by an intricate interplay of electrostatic and hydrodynamic interactions. Already systems composed of one single particle in electrolyte solution exhibit a complex dynamical behaviour. In this review, we discuss recent coarse-grained simulations that have been performed to obtain a molecular-level understanding of the dynamic and dielectric response of single particles and single macromolecules to external electric fields. We address both the response of charged particles to constant fields (DC fields), which can be characterized by an electrophoretic mobility, and the dielectric response of both uncharged and charged particles to alternating fields (AC fields), which is described by a complex polarizability. Furthermore, we give a brief survey of simulation algorithms and highlight some recent developments.

Convexity and Stiffness in Energy Functions for Electrostatic Simulations

We study the properties of convex functionals which have been proposed for the simulation of charged molecular systems within the Poisson-Boltzmann approximation. We consider the extent to which the functionals reproduce the true fluctuations of electrolytes and thus the one-loop correction to mean field theory-including the Debye-Hückel correction to the free energy of ionic solutions. We also compare the functionals for use in numerical optimization of a mean field model of a charged polymer and show that different functionals have very different stiffnesses leading to substantial differences in accuracy and speed.

Computing the Coulomb interaction in inhomogeneous dielectric media via a local electrostatics lattice algorithm

The local approach to computing electrostatic interactions proposed by Maggs and adapted by Rottler and Pasichnyk for molecular-dynamics simulations is extended to situations where the dielectric background medium is inhomogeneous. We furthermore correct a problem of the original algorithm related to the correct treatment of the global dipole moment, provide an error estimate for the accuracy of the algorithm, and suggest a different form of the treatment of the self-energy problem. Our implementation is highly scalable on many cores, and we have validated and compared its performance against theoretical predictions and simulation data obtained by other algorithmic approaches.