Modeling the surface charge evolution of spherical nanoparticles by considering dielectric discontinuity effects at the solid/electrolyte solution interface

It is well known that the electrostatic repulsions between charges on neighboring sites decrease the effective charge at the surface of a charged nanoparticle (NP). However, the situation is more complex close to a dielectric discontinuity, since charged sites are interacting not only with their neighbors but also with their own image charges and the image charges of all neighbors. Titrating site positions, solution ionic concentration, dielectric discontinuity effects, and surface charge variations with pH are investigated here using a grand canonical Monte Carlo method. A Tanford and Kirkwood approach is used to calculate the interaction potentials between the discrete charged sites. Homogeneous, heterogeneous, and patch site distributions are considered to reproduce the various titrating site distributions at the solid/solution interface of spherical NPs. By considering Coulomb, salt, and image charges effects, results show that for different ionic concentrations, modifications of the dielectric constant of NPs having homogeneous and heterogeneous site distributions have little effect on their charging process. Thus, the reaction field, due to the presence of image charges, fully counterbalances the Coulomb interactions. This is not the case for patch distributions, where Coulomb interactions are not completely counterbalanced by the reaction field. Application of the present model to pyrogenic silica is also performed and comparison is made with published experimental data of titration curves at various ionic concentrations.     

Extension of the fast multipole method for the rectangular cells with an anisotropic partition tree structure

The fast multipole method (FMM) is an order N method for the numerically rigorous calculation of the electrostatic interactions among point charges in a system of interest. The FMM is utilized for massively parallelized software for molecular dynamics (MD) calculations. However, an inconvenient limitation is imposed on the implementation of the FMM: In three-dimensional case, a cubic MD unit cell is hierarchically divided by the octree partitioning under isotropic periodic boundary conditions along three axes. Here, we extended the FMM algorithm adaptive to a rectangular MD unit cell with different periodicity along the axes by applying an anisotropic hierarchical partitioning. The algorithm was implemented into the parallelized general-purpose MD calculation software designed for a system with uniform distribution of point charges in the unit cell. The partition tree can be a mixture of binary and ternary branches, the branches being chosen arbitrarily with respect to the coordinate axes at any levels. Errors in the calculated electrostatic interactions are discussed in detail for a selected partition tree structure. The extension enables us to execute MD calculations under more general conditions for the shape of the unit cell, partition tree, and boundary conditions, keeping the accuracy of the calculated electrostatic interactions as high as that with the conventional FMM. An extension of the present FMM algorithm to other prime number branches, such as 5 and 7, is straightforward.     

Effects of surface site distribution and dielectric discontinuity on the charging behavior of nanoparticles: a grand canonical Monte Carlo study

The surface site distribution and the dielectric discontinuity effects on the charging process of a spherical nanoparticle (NP) have been investigated. It is well known that electrostatic repulsion between charges on neighbouring sites tends to decrease the effective charge of a NP. The situation is more complicated close to a dielectric breakdown, since here a charged site is not only interacting with its neighbours but also with its own image charge and the image charges of all its neighbours. Coexistence of opposite charges, titration sites positions, and pH dependence are systematically studied using a grand canonical Monte Carlo method. A Tanford and Kirkwood approach has been applied to describe the interaction potentials between explicit discrete ampholytic charging sites. Homogeneous, heterogeneous and patch site distributions were considered to reproduce the titration site distribution at the solid/solution interface of natural NPs. Results show that the charging process is controlled by the balance between Coulomb interactions and the reaction field through the solid-liquid interface. They also show that the site distribution plays a crucial role in the charging process. In patch distributions, charges accumulate at the perimeter of each patch due to finite size effects. When homogeneous and heterogeneous distributions are compared, three different charging regimes are obtained. In homogeneous and heterogeneous (with quite low polydispersity indexes) distributions, the effects of the NP dielectric constant on Coulomb interactions are counterbalanced by the reaction field and in this case, the dielectric breakdown has no significant effect on the charging process. This is not the case in patch distributions, where the dielectric breakdown plays a crucial role in the charging process.     

The electrostatic surface term: (I) periodic systems

The authors propose a new approach to understand the electrostatic surface contributions to the interactions of large but finite periodic distributions of charges. They present a simple method to derive and interpret the surface contribution to any electrostatic field produced by a periodic distribution of charges. They discuss the physical and mathematical interpretations of this term. They present several examples and physical details associated with the calculation of the surface term. Finally, they provide a simple derivation of the surface contribution to the virial. This term does not disappear even if tinfoil boundary conditions are applied.     

Conformational energies of trans-1,4-difluoro- and trans-1-chloro-4-fluorocyclohexane. The role of electrostatic interactions.

The conformational energies of the title compounds were extrapolated from solution ¹⁹F NMR measurements to vapour phase dielectric medium. Vapour phase values were satisfactorily reproduced by MM2 computations including electrostatic interactions between charges obtained by the PEOE method.     

Ionic electrets: electrostatic charging of surfaces by transferring mobile ions upon contact

This paper describes the fabrication and characterization of ionic electrets-materials that bear a long-lived electrostatic charge because of an imbalance between the number of cationic and anionic charges in the material. Crosslinked polystyrene microspheres that contain covalently bound ions and mobile counterions transfer some of their mobile ions in air, in the absence of bulk liquid, to another material upon contact. According to the ion-transfer model of contact electrification, this selective transfer of mobile ions yields microspheres that have a net electrostatic charge. A tool that operates on the principle of electrostatic induction measures the charge on individual microspheres (50-450 microm in diameter). Microspheres with a variety of covalently bound ionic functional groups (tetraalkylammonium, alkyltriphenylphosphonium, alkylsulfonate, and arylsulfonate) acquire charges consistent with this ion-transfer mechanism. The charge on a microsphere is proportional to its surface area (ca. 1 elementary charge per 2000 nm2) and close to the theoretical limit imposed by the dielectric breakdown of air. The charge density in an atmosphere of SF6 is more than twice that in an atmosphere of N2. These observations suggest that the charge density of these ionic electret microspheres is limited by the dielectric breakdown of the surrounding gas. Functionalizing the surfaces of glass or silicon with covalently bound ions and mobile counterions generates ionic electrets from these inorganic substrates. Soft lithography can pattern charge on a planar silicon surface (with oxide) and on the surface of 250-mum glass microspheres.     

Fragment-based quantum mechanical methods for periodic systems with Ewald summation and mean image charge convention for long-range electrostatic interactions

We describe an Ewald-summation method to incorporate long-range electrostatic interactions into fragment-based electronic structure methods for periodic systems. The present method is an extension of the particle-mesh Ewald technique for combined quantum mechanical and molecular mechanical (QM/MM) calculations, and it has been implemented into the explicit polarization (X-Pol) potential to illustrate the computational details. As in the QM/MM-Ewald method, the X-Pol-Ewald approach is a linear-scaling electrostatic method, in which the short-range electrostatic interactions are determined explicitly in real space and the long-range Ewald pair potential is incorporated into the Fock matrix as a correction. To avoid the time-consuming Fock matrix update during the self-consistent field procedure, a mean image charge (MIC) approximation is introduced, in which the running average with a user-chosen correlation time is used to represent the long-range electrostatic correction as an average effect. Test simulations on liquid water show that the present X-Pol-Ewald method takes about 25% more CPU time than the usual X-Pol method using spherical cutoff, whereas the use of the MIC approximation reduces the extra costs for long-range electrostatic interactions by 15%. The present X-Pol-Ewald method provides a general procedure for incorporating long-range electrostatic effects into fragment-based electronic structure methods for treating biomolecular and condensed-phase systems under periodic boundary conditions.     

A self-organized 2-dimensional bifunctional structure formed by supramolecular design

The production of a novel two-dimensional bimolecular surface structure using weak noncovalent interactions is demonstrated and observed by scanning tunneling microscopy. This work follows the three-dimensional (3D) ideas of crystal engineering and applies the concepts of supramolecular synthons to molecular systems constrained to 2D by physisorption on a conducting surface. We demonstrate a well-ordered planar structure that self-assembles through the influence of fluorine-phenyl interactions. This study provides a concrete example of the "bottom up" construction of nanostructures and of the systematic design of self-organized layers. To our knowledge, this is the first in a new class of fully 2D materials based both upon weak intermolecular interactions and upon image charges and weak interactions associated with adsorption on metal surfaces.     

Calculation of the microscopic and macroscopic linear and nonlinear optical properties of liquid acetonitrile. II. Local fields and linear and nonlinear susceptibilities in quadrupolar approximation

A discrete model based on the multipolar expansion including terms up to hexadecapoles was employed to describe the electrostatic interactions in liquid acetonitrile. Liquid structures obtained form molecular dynamics simulations with different classical, nonpolarizable potentials were used to analyze the electrostatic interactions. The computed average local field was employed for the determination of the environmental effects on the linear and nonlinear electrical molecular properties. Dipole-dipole interactions yield the dominant contribution to the local field, whereas higher multipolar contributions are small but not negligible. Using the effective in-phase properties, macroscopic linear and nonlinear susceptibilities of the liquid were computed. Depending on the partial charges describing the Coulomb interactions of the force field employed, either the linear properties (refractive index and dielectric constant) were reproduced in good agreement with experiment or the nonlinear properties [third-harmonic generation (THG) and electric field induced second-harmonic (EFISH) generation] and the bulk density but never both sets of properties together. It is concluded that the partial charges of the force fields investigated are not suitable for reliable dielectric properties. New methods are probably necessary for the determination of partial charges, which should take into account the collective and long-range nature of electrostatic interactions more precisely.     

甲硫氨酸-脑啡肽的活性位点

用ＡＭ１方法中计算静电势的ＰＭＥＰ子程序了甲硫氨酸－脑啡肽的活性位点。通过计算得到了整个分子的三维空间静电势分布和由静电势导出的各原子的电荷分布。进一步分析,确定模型分子活性位点为酪氨酸残在的叔氨原子和苯酚基、苯丙氨酸残基的苯基及部分氧原子。     

Ligand-Mediated Regioselective Rhodium-Catalyzed Benzotriazole–Allene Coupling: Mechanistic Exploration and Quantum Chemical Analysis

The ligand-controlled rhodium-catalyzed regioselective coupling of 1,2,3-benzotriazoles and allenes was investigated by DFT calculations. Because allylation can occur at either the N1 or N2 position of the 1,2,3-benzotriazole, the complete Gibbs free energy profiles for both pathways were computed. A kinetic preference emerged for the experimentally observed N1 allylation with the JoSPOphos ligand, whereas N2 allylation was favored with DPEphos. Analysis of the regiodetermining oxidative addition step by using the activation strain model in conjunction with a matching energy decomposition analysis has revealed that the unprecedented N2 reaction regioselectivity is dictated by the strength of the electrostatic interactions between the 1,2,3-benzotriazole and the rhodium catalyst. The nature of the electrostatic interaction was rationalized by analysis of the electrostatic potential maps and Hirshfeld charges: a stabilizing electrostatic interaction was found between the key atoms involved in the oxidative addition for the N2 pathway, analogous interactions are weaker in the N1 case.     

用第一性原理方法获取周期体系中原子的部分电荷

利用量子化学软件包Crystal计算了立方周期性边界条件下液态水体系的静电势(ESP)和静电场(EF).在此基础上,提出了一种由第一性原理方法获取周期体系中原子的部分电荷的快捷方法.该方法把由周期性边界条件引入的平均静电势准mean作为一个拟合参数,通过对第一性原理静电势与Ewald加和法静电势的最小二乘法拟合而实现.值得说明的是,比较静电势与静电场拟合方法,前者的相对拟合误差仅为2%-3%,比后者小一个数量级.考察了四种电荷限制条件下,静电势、静电场拟合的水分子原子部分电荷及偶极矩的分布情况.     

Long-range Lennard-Jones and electrostatic interactions in interfaces: application of the isotropic periodic sum method

Molecular dynamics (MD) simulations of heptane/vapor, hexadecane/vapor, water/vapor, hexadecane/water, and dipalmitoylphosphatidylcholine (DPPC) bilayers and monolayers are analyzed to determine the accuracy of treating long-range interactions in interfaces with the isotropic periodic sum (IPS) method. The method and cutoff (rc) dependences of surface tensions, density profiles, water dipole orientation, and electrostatic potential profiles are used as metrics. The water/vapor, heptane/vapor, and hexadecane/vapor interfaces are accurately and efficiently calculated with 2D IPS (rc=10 A). It is demonstrated that 3D IPS is not practical for any of the interfacial systems studied. However, the hybrid method PME/IPS [Particle Mesh Ewald for electrostatics and 3D IPS for Lennard-Jones (LJ) interactions] provides an efficient way to include both types of long-range forces in simulations of large liquid/vacuum and all liquid/liquid interfaces, including lipid monolayers and bilayers. A previously published pressure-based long-range LJ correction yields results similar to those of PME/IPS for liquid/liquid interfaces. The contributions to surface tension of LJ terms arising from interactions beyond 10 A range from 13 dyn/cm for the hexadecane/vapor interface to approximately 3 dyn/cm for hexadecane/water and DPPC bilayers and monolayers. Surface tensions of alkane/vapor, hexadecane/water, and DPPC monolayers based on the CHARMM lipid force fields agree very well with experiment, whereas surface tensions of the TIP3P and TIP4P-Ew water models underestimate experiment by 16 and 11 dyn/cm, respectively. Dipole potential drops (DeltaPsi) are less sensitive to long-range LJ interactions than surface tensions. However, DeltaPsi for the DPPC bilayer (845+/-3 mV proceeding from water to lipid) and water (547+/-2 mV for TIP4P-Ew and 521+/-3 mV for TIP3P) overestimate experiment by factors of 3 and 5, respectively, and represent expected deficiencies in nonpolarizable force fields.     

MMM1D: a method for calculating electrostatic interactions in one-dimensional periodic geometries

We present a new method to accurately calculate the electrostatic energy and forces on charges in a system with periodic boundary conditions in one of three spatial dimensions. We transform the Coulomb sum via a convergence factor into a series of fast decaying functions similar to the Lekner method. Rigorous error bounds for the energies and the forces are derived and numerically verified. The method has a computational complexity of O(N(2)), but is faster and easier to use than previously reported methods.     

Bipolar Tribocharging Signal During Friction Force Fluctuations at Metal–Insulator Interfaces

Friction and triboelectrification of materials show a strong correlation during sliding contacts. Friction force fluctuations are always accompanied by two tribocharging events at metal–insulator [e.g., polytetrafluoroethylene (PTFE)] interfaces: injection of charged species from the metal into PTFE followed by the flow of charges from PTFE to the metal surface. Adhesion maps that were obtained by atomic force microscopy (AFM) show that the region of contact increases the pull‐off force from 10 to 150 nN, reflecting on a resilient electrostatic adhesion between PTFE and the metallic surface. The reported results suggest that friction and triboelectrification have a common origin that must be associated with the occurrence of strong electrostatic interactions at the interface.     

A universal approach to fabricate ordered colloidal crystals arrays based on electrostatic self-assembly

We present a novel and simple method to fabricate two-dimensional (2D) poly(styrene sulfate) (PSS, negatively charged) colloidal crystals on a positively charged substrate. Our strategy contains two separate steps: one is the three-dimensional (3D) assembly of PSS particles in ethanol, and the other is electrostatic adsorption in water. First, 3D assembly in ethanol phase eliminates electrostatic attractions between colloids and the substrate. As a result, high-quality colloidal crystals are easily generated, for electrostatic attractions are unfavorable for the movement of colloidal particles during convective self-assembly. Subsequently, top layers of colloidal spheres are washed away in the water phase, whereas well-packed PSS colloids that are in contact with the substrate are tightly linked due to electrostatic interactions, resulting in the formation of ordered arrays of 2D colloidal spheres. Cycling these processes leads to the layer-by-layer assembly of 3D colloidal crystals with controllable layers. In addition, this strategy can be extended to the fabrication of patterned 2D colloidal crystals on patterned polyelectrolyte surfaces, not only on planar substrates but also on nonplanar substrates. This straightforward method may open up new possibilities for practical use of colloidal crystals of excellent quality, various patterns, and controllable fashions.     

Simulating van der Waals interactions in water/hydrocarbon-based complex fluids

In systems composed of water and hydrocarbons, Van der Waals interactions are dominated by the nonretarded, classical (Keesom) part of the Lifshitz interaction; the interaction is screened by salt and extends over mesoscopic distances of the order of the size of the (micellar) constituents of complex fluids. We show that these interactions are included intrinsically in a recently introduced local Monte Carlo algorithm for simulating electrostatic interactions between charges in the presence of nonhomogeneous dielectric media.