An Iterative Fragment Scheme for the ACKS2 Electronic Polarization Model: Application to Molecular Dimers and Chains.

The treatment of electrostatic interactions is a key ingredient in the force field-based simulation of condensed phase systems. Most approaches used fixed, site-specific point charges. Yet, it is now clear that many applications of force fields (FFs) demand more sophisticated treatments, prompting the implementation of charge equilibration methods in polarizable FFs to allow the redistribution of charge within the system. One approach allowing both, charge redistribution and site-specific polarization, while at the same time solving methodological shortcomings of earlier methods, is the first-principles-derived atom-condensed Kohn-Sham density functional theory method approximated to the second order (ACKS2). In this work, we present two fragment approaches to ACKS2, termed f-ACKS2 and a self-consistent version, scf-ACKS2, that treat condensed phase systems as a collection of electronically polarizable molecular fragments. The fragmentation approach to ACKS2 not only leads to a more transferable and less system-specific collection of electronic response parameters but also opens up the method to large condensed phase systems. We validate the accuracies of f-ACKS2 and scf-ACKS2 by comparing polarization energies and induced dipole moments for a number of charged hydrocarbon dimers against DFT reference calculations. Finally, we also apply both fragmented ACKS2 variants to calculate the polarization energy for electron-hole pair separation along a chain of anthracene molecules and find excellent agreement with reference DFT calculations.

Toward First-Principles-Level Polarization Energies in Force Fields: A Gaussian Basis for the Atom-Condensed Kohn-Sham Method.

The last 20 years of force field development have shown that even well parametrized classical models need to at least approximate the dielectric response of molecular systems-based, e.g., on atomic polarizabilities-in order to correctly render their structural and dynamic properties. Yet, despite great advances most approaches tend to be based on ad hoc assumptions and often insufficiently capture the dielectric response of the system to external perturbations, such as, e.g., charge carriers in semiconducting materials. A possible remedy was recently introduced with the atom-condensed Kohn-Sham density-functional theory approximated to second order (ACKS2), which is fully derived from first principles. Unfortunately, specifically its reliance on first-principles derived parameters so far precluded the widespread adoption of ACKS2. Opening up ACKS2 for general use, we here present a reformulation of the method in terms of Gaussian basis functions, which allows us to determine many of the ACKS2 parameters analytically. Two sets of parameters depending on exchange-correlation interactions are still calculated numerically, but we show that they could be straightforwardly parametrized owing to the smoothness of the new basis. Our approach exhibits three crucial benefits for future applications in force fields: i) efficiency, ii) accuracy, and iii) transferability. We numerically validate our Gaussian augmented ACKS2 model for a set of small hydrocarbons which shows a very good agreement with density-functional theory reference calculations. To further demonstrate the method's accuracy and transferability for realistic systems, we calculate polarization responses and energies of anthracene and tetracene, two major building blocks in organic semiconductors.

Energy Dissipation during Diffusion at Metal Surfaces: Disentangling the Role of Phonons versus Electron-Hole Pairs.

Helium spin echo experiments combined with ab initio based Langevin molecular dynamics simulations are used to quantify the adsorbate-substrate coupling during the thermal diffusion of Na atoms on Cu(111). An analysis of trajectories within the local density friction approximation allows the contribution from electron-hole pair excitations to be separated from the total energy dissipation. Despite the minimal electronic friction coefficient of Na and the relatively small mass mismatch to Cu promoting efficient phononic dissipation, about (20±5)% of the total energy loss is attributable to electronic friction. The results suggest a significant role of electronic nonadiabaticity in the rapid thermalization generally relied upon in adiabatic diffusion theories.