A Charge−Dipole Model for the Static Polarizability of Nanostructures Including Aliphatic, Olephinic, and Aromatic Systems

Demonstrating the Connection between the Nonvalence Correlation-Bound Anions of Polyaromatic Hydrocarbons and the Image Potential States of Graphene Using a One-Electron Model Hamiltonian

The ground and excited state nonvalence correlation-bound (NVCB) anion states of the C 6 n 2 H 6 n hexagonal polycylic aromatic hydrocarbons and of hexagonal C 6 n 2 graphene nanoflakes are characterized using a one-electron model Hamiltonian which incorporates atomic electrostatic moments up to the quadrupole, coupled inducible charges and dipoles, and atom-centered Gaussians to describe the short-range repulsive interactions. Extrapolation of the calculated electron binding energies of the lowest energy symmetric and antisymmetric (with respect to the molecular plane) NVCB anions of both the polycylic aromatic hydrocarbons and the carbon nanoflakes to the n → ∞ limit yields binding energies that are in good agreement with those of the most stable symmetric and antisymmetric image potential states of freestanding graphene as determined from two-photon photoemission spectroscopy (2PPE) experiments.

Application of a Fluctuating Charge Polarization Model to Large Polyaromatic Hydrocarbons and Graphene Nanoflakes

We present a polarization model incorporating coupled fluctuating charges and point inducible dipoles that is able to accurately describe the dipole polarizabilities of small hydrocarbons and, for sufficiently large graphene nanoflakes, reproduce the classical image potential of an infinite conducting sheet. When our fluctuating charge model is applied to the hexagonal carbon nanoflake C60000 we attain excellent agreement with the image potential and induced charge distribution of a conducting sheet. With the inclusion of inducible dipole terms, the model predicts an image plane of zim = 1.3334 a0, which falls in line with prior estimates for graphene. We consider the case of two charges placed on opposite sides of C60000 and find that the fluctuating charge model reproduces classical electrostatics once again. By testing opposing and similar signs of the external charges, we conclude that an atomically thin molecule or extended system does not fully screen their interaction.

Data-driven tailoring of molecular dipole polarizability and frontier orbital energies in chemical compound space

Understanding correlations - or lack thereof - between molecular properties is crucial for enabling fast and accurate molecular design strategies. In this contribution, we explore the relation between two key quantities describing the electronic structure and chemical properties of molecular systems: the energy gap between the frontier orbitals and the dipole polarizability. Based on the recently introduced QM7-X dataset, augmented with accurate molecular polarizability calculations as well as analysis of functional group compositions, we show that polarizability and HOMO-LUMO gap are uncorrelated when considering sufficiently extended subsets of the chemical compound space. The relation between these two properties is further analyzed on specific examples of molecules with similar composition as well as homooligomers. Remarkably, the freedom brought by the lack of correlation between molecular polarizability and HOMO-LUMO gap enables the design of novel materials, as we demonstrate on the example of organic photodetector candidates.

Computation of Flexoelectric Coefficients of a MoS2 monolayer with a Model of Self-consistently Distributed Effective Charges and Dipoles

Flexoelectricity is an electromechanical coupling phenomenon, that can generate noticeable electric polarization in dielectric materials for nanoscale strain gradients. It is gaining an increasing attention because of its potential applications, and the fact that experimental results were initially an order of magnitude higher than initial theoretical predictions. This stimulated intense experimental and theoretical researches to investigate flexoelectric coefficients in dielectric materials such as two-dimensional materials. In this work, we concentrate on the calculation of the flexoelectric coefficients of 2D-MoS2 thanks to a model using self-consistently determined charges and dipoles on the atoms. More specifically, we study the importance of two contributions which were neglected/omitted in previous papers using this model, namely the charge term in the total polarization and the conservation of electric charge through a Lagrange multiplier. Our calculations demonstrate that the results for flexoelectric coefficient computed with this improved definition of polarization agree better with experimental measurements, provided consistent definitions for signs are used. Additionally, we show how two physical contributions with opposite signs compete to give net values of flexoelectric coefficients that can be either positive of negative depending on their relative importance, and give net values for the case of MoS2.

Effect of carbon doping on electromechanical response of boron nitride nanosheets

The electromechanical response of hexagonal-boron nitride nanosheets (h-BNSs) was studied via molecular dynamics simulations (MDS) with a three-body Tersoff potential force field using a charge-dipole (C-D) potential model. Carbon (C)-doped h-BNSs with triangular, trapezoidal and circular pores were considered. The elastic and piezoelectric coefficients of h-BNSs under tension and shear loading conditions were determined. The induced polarization in h-BNSs was found to depend on the local arrangement of C atoms around B and N atoms, and the polarization increases if C atoms are surrounded by N atoms. This is attributed to the generation of higher dipole moments due to C-N bonds compared with C-B bonds. At ∼5.5% C-doping concentration, the axial piezoelectric coefficient of doped h-BNSs with triangular and trapezoidal pores increased by 18.5% and 3.5%, respectively, while it reduced by 22.5% in the case of circular pores compared to pristine h-BNS. The shear piezoelectric coefficient of C-doped h-BNSs with triangular and trapezoidal pores increased by 20.5% and 1%, respectively, while it reduced by 7% in case of circular pores. Young's moduli of C-doped h-BNSs with triangular, trapezoidal and circular pores increased by 9%, 7.5% and 5.5%, respectively, due to the C-C bonds being stronger than all other bonds. The respective improvements in shear moduli are 8.5%, 5% and 5%. The elastic and piezoelectric properties of armchair h-BNSs were found to be higher than zigzag h-BNSs. The results also reveal that the piezoelectric coefficient increases as doping increases; it reaches its maximum value around 0.41 C m^-2 at 12.6% C-doping concentration and then starts decreasing. The present work shows that we can engineer the electromechanical response of h-BNSs via novel pathways such as different types and size of pores as well as C-doping concentration to suit a particular nanoelectromechanical systems (NEMS) application.

New scaling relations to compute atom-in-material polarizabilities and dispersion coefficients: part 1. Theory and accuracy

A new method was developed to compute atom-in-material polarizabilities and dispersion coefficients for diverse material types.

Model potential study of non-valence correlation-bound anions of (C60)n clusters: the role of electric field-induced charge transfer

Orbitals associated with the non-valence correlation-bound anions of the C₆₀ dimer and linear trimer from calculations allowing charge transfer.

An atomistic model for the charge distribution in layered MoS2

We present an atomistic model for predicting the distribution of doping electric charges in layered molybdenum disulfide (MoS₂). This model mimics the charge around each ion as a net Gaussian-spatially distributed charge plus an induced dipole, and is able to predict the distribution of doping charges in layered MoS₂ in a self-consistent scheme. The profiles of doping charges in monolayer MoS₂ flakes computed by this charge-dipole model are in good agreement with those obtained by density-functional-theory calculations. Using this model, we quantitatively predict the charge enhancement in MoS₂ monolayer nanoribbons, with which strong ionic charge-localization effects are shown.

Molecular adsorbates as probes of the local properties of doped graphene

Graphene-based sensors are among the most promising of graphene's applications. The ability to signal the presence of molecular species adsorbed on this atomically thin substrate has been explored from electric measurements to light scattering. Here we show that the adsorbed molecules can be used to sense graphene properties. The interaction of porphyrin molecules with nitrogen-doped graphene has been investigated using scanning tunneling microscopy and ab initio calculations. Molecular manipulation was used to reveal the surface below the adsorbed molecules allowing to achieve an atomic-scale measure of the interaction of molecules with doped graphene. The adsorbate's frontier electronic states are downshifted in energy as the molecule approaches the doping site, with largest effect when the molecule sits over the nitrogen dopant. Theoretical calculations showed that, due to graphene's high polarizability, the adsorption of porphyrin induces a charge rearrangement on the substrate similar to the image charges on a metal. This charge polarization is enhanced around nitrogen site, leading to an increased interaction of molecules with their image charges on graphene. Consequently, the molecular states are stabilized and shift to lower energies. These findings reveal the local variation of polarizability induced by nitrogen dopant opening new routes towards the electronic tuning of graphene.

Local electric field factors by a combined charge-transfer and point–dipole interaction model

A model for the local electric field as a linear response to a frequency-dependent external electric field is presented based on a combined charge-transfer and point–dipole interaction force-field model.

Nonvalence correlation-bound anion states of spherical fullerenes

We present a one-electron model Hamiltonian for characterizing nonvalence correlation-bound anion states of fullerene molecules. These states are the finite system analogs of image potential states of metallic surfaces. The model potential accounts for both atomic and charge-flow polarization and is used to characterize the nonvalence correlation-bound anion states of the C60, (C60)2, C240, and C60@C240 fullerene systems. Although C60 is found to have a single (s-type) nonvalence correlation-bound anion state, the larger fullerenes are demonstrated to have multiple nonvalence correlation-bound anion states.

A periodic charge-dipole electrostatic model. II. A kinetic-exchange-correlation correction

We extend the periodic charge-dipole electrostatic model, see I. V. Bodrenko, M. Sierka, E. Fabiano, and F. Della Sala, J. Chem. Phys. 137, 134702 (2012), to include a kinetic-exchange-correlation (KXC) correction. The KXC correction is approximated by means of an extended-Hückel-type formula, it is exact in the infinite jellium model and it is also computationally efficient as it requires only the computation of overlap integrals. Tests on the linear response of silver slabs to an external electrostatic perturbation show that the KXC correction yields a very accurate description of induced dipole and of the whole induced charge density profile. We also show that the KXC parameters are quite transferable and related to the atomic polarizability.

A periodic charge-dipole electrostatic model: parametrization for silver slabs

We present an extension of the charge-dipole model for the description of periodic systems. This periodic charge-dipole electrostatic model (PCDEM) allows one to describe the linear response of periodic structures in terms of charge- and dipole-type gaussian basis functions. The long-range electrostatic interaction is efficiently described by means of the continuous fast multipole method. As a first application, the PCDEM method is applied to describe the polarizability of silver slabs. We find that for a correct description of the polarizability of the slabs both charges and dipoles are required. However a continuum set of parametrizations, i.e., different values of the width of charge- and dipole-type gaussians, leads to an equivalent and accurate description of the slabs polarizability but a completely unphysical description of induced charge-density inside the slab. We introduced the integral squared density measure which allows one to obtain a unique parametrization which accurately describes both the polarizability and the induced density profile inside the slab. Finally the limits of the electrostatic approximations are also pointed out.

Accurate and efficient method for many-body van der Waals interactions

An efficient method is developed for the microscopic description of the frequency-dependent polarizability of finite-gap molecules and solids. This is achieved by combining the Tkatchenko-Scheffler van der Waals (vdW) method [Phys. Rev. Lett. 102, 073005 (2009)] with the self-consistent screening equation of classical electrodynamics. This leads to a seamless description of polarization and depolarization for the polarizability tensor of molecules and solids. The screened long-range many-body vdW energy is obtained from the solution of the Schrödinger equation for a system of coupled oscillators. We show that the screening and the many-body vdW energy play a significant role even for rather small molecules, becoming crucial for an accurate treatment of conformational energies for biomolecules and binding of molecular crystals. The computational cost of the developed theory is negligible compared to the underlying electronic structure calculation.

Molecular modeling and dynamics studies with explicit inclusion of electronic polarizability. Theory and applications

A current emphasis in empirical force fields is on the development of potential functions that explicitly treat electronic polarizability. In the present article, the commonly used methodologies for modelling electronic polarization are presented along with an overview of selected application studies. Models presented include induced point-dipoles, classical Drude oscillators, and fluctuating charge methods. The theoretical background of each method is followed by an introduction to extended Langrangian integrators required for computationally tractable molecular dynamics simulations using polarizable force fields. The remainder of the review focuses on application studies using these methods. Emphasis is placed on water models, for which numerous examples exist, with a more thorough discussion presented on the recently published models associated with the Drude-based CHARMM and the AMOEBA force fields. The utility of polarizable models for the study of ion solvation is then presented followed by an overview of studies of small molecules (e.g. CCl(4), alkanes, etc) and macromolecule (proteins, nucleic acids and lipid bilayers) application studies. The review is written with the goal of providing a general overview of the current status of the field and to facilitate future application and developments.

Enhanced polarizability of aromatic molecules placed in the vicinity of silver clusters

We use a charge-dipole interaction model to study the polarizability of aromatic molecules that are placed between two silver clusters. In particular we examine the enhancement in polarizability induced by the clusters at plasmon-like resonant frequencies of the cluster-molecule-cluster system. The model used for these simulations relies on representation of the atoms by both a net electric charge and a dipole. By relating the time variation of the atomic charges to the currents that flow through the bonds of the structures considered, a least-action principle can be formulated that enables the atomic charges and dipoles to be determined. We consider benzene, naphthalene and anthracene for this study, comparing the polarizability of these aromatic molecules when placed in the middle between two Ag(120) clusters, with their polarizability as isolated molecules. We find that the polarizability of these molecules is enhanced by the clusters, and this increases the electromagnetic coupling between the two clusters. This results in significant red-shifting (by up to 0.8 eV) of the lowest energy optical transition in the cluster-molecule-cluster system compared to plasmon-like excitation in the cluster-cluster system. The resulting resonant polarizability enhancement leads to an electromagnetic enhancement in surface-enhanced Raman scattering of over 10(6).

A charge-dipole interaction model for the frequency-dependent polarizability of silver clusters

We present a charge-dipole interaction model for the calculation of the frequency-dependent polarizability of silver clusters. The model relies on the representation of silver atoms by both a net electric charge and a dipole. Time variations of the atomic charges are related to the currents that flow through the bonds of the structures considered and the atomic charges and dipoles are eventually determined from the application of a least-action principle. After a generalization that enables the bonds of the bulk and surface atoms to have specific resistances, the model is parameterized on data obtained by the time-dependent density functional theory for tetrahedral Ag(20), Ag(84) and Ag(120) clusters. We then study the polarization properties of dimers of silver clusters. We compare in particular the polarizability of the dimers with that of the isolated clusters, for a range of gap distances and frequencies. We also consider the field enhancements one can achieve with these systems. The results are in good agreement with reference data and enable an extension of these data to a wider range of situations. They show that significant field enhancements are achieved at frequencies associated with resonant polarization along the axis of the dimer.