An official website of the United States government
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 is 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 periodic boundary conditions correction scheme for nonperiodic 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 homogeneous 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.
more »
« less
This content will become publicly available on October 28, 2026
Coupling all-electron full-potential density functional theory with grid-based continuum embeddings
Recent advances in continuum embedding models have enabled the incorporation of solvent and electrolyte effects into density functional theory (DFT) simulations of material surfaces, significantly benefiting electrochemistry, catalysis, and other applications. To extend the simulation of diverse systems and properties, the implementation of continuum embedding models into the Environ library adopts a modular programming paradigm, offering a flexible interface for communication with various DFT programs. The speed and scalability of the current implementation rely on a smooth definition of the key physical properties of the atomistic system, in particular, of its electronic density. This has hindered the coupling of Environ with all-electron simulation packages, as the sharp electron density peaks near atomic nuclei are difficult to represent on regular grids. In this work, we introduce a novel smoothing scheme that transforms atom-centered electron densities into a regular grid representation while preserving the accuracy of electrostatic calculations. This approach enables a minimal and generic interface, facilitating seamless interoperability between Environ and all-electron DFT programs. We demonstrate this development through the coupling of Environ with the FHI-aims package and present benchmark simulations that validate the proposed method.
more »
« less
- PAR ID:
- 10656253
- Publisher / Repository:
- AIP
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 163
- Issue:
- 16
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Both experimental results and theoretical models suggest the decisive role of the filler–matrix interfaces on the dielectric, piezoelectric, pyroelectric, and electrocaloric properties of ferroelectric polymer nanocomposites. However, there remains a lack of direct structural evidence to support the so‐called interfacial effect in dielectric nanocomposites. Here, a chemical mapping of the interfacial coupling between the nanofiller and the polymer matrix in ferroelectric polymer nanocomposites by combining atomic force microscopy–infrared spectroscopy (AFM–IR) with first‐principles calculations and phase‐field simulations is provided. The addition of ceramic fillers into a ferroelectric polymer leads to augmentation of the local conformational disorder in the vicinity of the interface, resulting in the local stabilization of the all‐transconformation (i.e., the polar β phase). The formation of highly polar and inhomogeneous interfacial regions, which is further enhanced with a decrease of the filler size, has been identified experimentally and verified by phase‐field simulations and density functional theory (DFT) calculations. This work offers unprecedented structural insights into the configurational disorder‐induced interfacial effect and will enable rational design and molecular engineering of the filler–matrix interfaces of electroactive polymer nanocomposites to boost their collective properties.more » « less
-
Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin–orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm–Dancoff approximation to TDDFT, spin–flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine–Ag 20 complex), we emphasize that caution is needed to properly account for system–environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.more » « less
-
The hydrated electron is of interest to both theorists and experimentalists as a paradigm solution-phase quantum system. Although the bulk of the theoretical work studying the hydrated electron is based on mixed quantum/classical (MQC) methods, recent advances in computer power have allowed several attempts to study this object using ab initio methods. The difficulty with employing ab initio methods for this system is that even with relatively inexpensive quantum chemistry methods such as density functional theory (DFT), such calculations are still limited to at most a few tens of water molecules and only a few picoseconds duration, leaving open the question as to whether the calculations are converged with respect to either system size or dynamical fluctuations. Moreover, the ab initio simulations of the hydrated electron that have been published to date have provided only limited analysis. Most works calculate the electron’s vertical detachment energy, which can be compared to experiment, and occasionally the electronic absorption spectrum is also computed. Structural features, such as pair distribution functions, are rare in the literature, with the majority of the structural analysis being simple statements that the electron resides in a cavity, which are often based only on a small number of simulation snapshots. Importantly, there has been no ab initio work examining the temperature-dependent behavior of the hydrated electron, which has not been satisfactorily explained by MQC simulations. In this work, we attempt to remedy this situation by running DFT-based ab initio simulations of the hydrated electron as a function of both box size and temperature. We show that the calculated properties of the hydrated electron are not converged even with simulation sizes up to 128 water molecules and durations of several tens of picoseconds. The simulations show significant changes in the water coordination and solvation structure with box size. Our temperature-dependent simulations predict a red-shift of the absorption spectrum (computed using TD-DFT with an optimally tuned range-separated hybrid functional) with increasing temperature, but the magnitude of the predicted red-shift is larger than that observed experimentally, and the absolute position of the calculated spectra are off by over half an eV. The spectral red-shift at high temperatures is accompanied by both a partial loss of structure of the electron’s central cavity and an increased radius of gyration that pushes electron density onto and beyond the first solvation shell. Overall, although ab initio simulations can provide some insights into the temperature-dependent behavior of the hydrated electron, the simulation sizes and level of quantum chemistry theory that are currently accessible are inadequate for correctly describing the experimental properties of this fascinating object.more » « less
-
Density functional theory (DFT) prediction of nuclear magnetic resonance (NMR) chemical shifts complements NMR experiments. Predicting chemical shifts accurately with DFT requires many different modeling decisions. Intended for novice modelers and nonexperts, this article discusses the considerations one should take in selecting a density functional, van der Waals dispersion correction, and basis set. It examines different strategies for handling systems in complex environments such as liquids, biomolecules, and crystals. Strategies include the use of cluster models, electrostatic embedding, continuum representations, periodic boundary conditions, and fragment‐based approaches. Finally, approaches for referencing the predicted absolute chemical shieldings for comparison against experimental chemical shifts are discussed.more » « less
