Uncoupling System and Environment Simulation Cells for Fast-Scaling Modeling of Complex Continuum Embeddings

Abstract

Continuum solvation models are becoming increasingly relevant in condensed matter simulations, allowing to characterize materials interfaces in the presence of wet electrified environments at a reduced computational cost with respect to all atomistic simulations. However, some challenges with the implementation of these models in plane-wave simulation packages still persists, especially when the goal is to simulate complex and heterogeneous environments. Among these challenges is the computational cost associated with large heterogeneous environments, which in plane-wave simulations has a direct effect on the basis-set size and, as a result, on the cost of the electronic structure calculation. Moreover, the use of periodic simulation cells are not well-suited for modeling systems embedded in semi-infinite media, which is often the case in continuum solvation models. To address these challenges, we present the implementation of a double-cell formalism, in which the simulation cell used for the continuum environment is uncoupled from the one used for the electronic-structure simulation of the quantum-mechanical system. This allows for a larger simulation cell to be used for the environment, without significantly increasing computational time. In this work, we show how the double-cell formalism can be used as an effective PBC correction scheme for non-periodic and partially periodic systems. The accuracy of the double-cell formalism is tested using representative examples with different dimensionalities, both in vacuum and in a continuum dielectric environment. Fast convergence and good speedups are observed for all the simulation setups, provided the quantum-mechanical simulation cell is chosen to completely fit the electronic density of the system.

References in corpus (20)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Advanced capabilities for materials modelling with Quantum ESPRESSO
• Maximally-localized generalized Wannier functions for composite energy bands
• AFLOW: An automatic framework for high-throughput materials discovery
• Precision and efficiency in solid-state pseudopotential calculations
• Revised self-consistent continuum solvation in electronic-structure calculations
• Grand canonical electronic density-functional theory: algorithms and applications to electrochemistry
• JDFTx: software for joint density-functional theory
• AiiDA 1.0, a scalable computational infrastructure for automated reproducible workflows and data provenance
• Workflows in AiiDA: Engineering a high-throughput, event-based engine for robust and modular computational workflows
• Electrostatics in Periodic Boundary Conditions and Real-space Corrections
• Hybrid-functional calculations with plane-wave basis sets: The effect of the singularity correction on total energies, energy eigenvalues, and defect energy levels
• A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments
• Evaluating continuum solvation models for the electrode-electrolyte interface: challenges and strategies for improvement
• Electrostatics of solvated systems in periodic boundary conditions
• Continuum models of the electrochemical diffuse layer in electronic-structure calculations
• Function-Space Based Solution Scheme for the Size-Modified Poisson-Boltzmann Equation in Full-Potential DFT
• eQE 2.0: Subsystem DFT Beyond GGA Functionals
• Electronic levels and electrical response of periodic molecular structures from plane-wave orbital-dependent calculations
• Quantum-continuum calculation of the surface states and electrical response of silicon in solution