Block-localized Wave Function Research Articles

Solvent Effect on Dative and Ionic Bond Strengths: A Unified Theory from Potential Analysis and Valence-Bond Computations.

Solvent molecules interact with a solute through various intermolecular forces. Here we employed a potential energy surface (PES) analysis to interpret the solvent-induced variations in the strengths of dative (Me3NBH3) and ionic (LiCl) bonds, which possess both ionic and covalent (neutral) characteristics. The change of a bond is driven by the gradient (force) of the solvent-solute interaction energy with respect to the focused bond length. Positive force shortens the bond length and increases the bond force constant, leading to a blue-shift of the bond stretching vibrational frequency upon solvation. Conversely, negative force elongates the bond, resulting in a reduced bond force constant and red-shift of the stretching vibrational frequency. The different responses of Me3NBH3 and LiCl to solvation are studied with valence bond (VB) theory, as Me3NBH3 and LiCl are dominated by the neutral covalent VB structure and the ionic VB structure, respectively. The dipole moment of an ionic VB structure increases along the increasing bond distance, while the dipole moment of a neutral covalent VB structure increases with the decreasing bond distance. The roles of the dominating VB structures are further examined by the geometry optimizations and frequency calculations with the block-localized wavefunction (BLW) method.

Conformational Preference of Lithium Polysulfide Clusters Li2Sx (x = 4-8) in Lithium-Sulfur Batteries.

Structures are of fundamental importance for diverse studies of lithium polysulfide clusters, which govern the performance of lithium-sulfur batteries. The ring-like geometries were regarded as the most stable structures, but their physical origin remains elusive. In this work, we systematically explored the minimal structures of Li2Sx (x = 4-8) clusters to uncover the driving force for their conformational preferences. All low-lying isomers were generated by performing global searches using the ABCluster program, and the ionic nature of the Li···S interactions was evidenced with the energy decomposition analysis based on the block-localized wave function (BLW-ED) approach and further confirmed with the quantum theory of atoms in molecule (QTAIM). By analysis of the contributions of various energy components to the relative stability with the references of the lowest-lying isomers, the controlling factor for isomer preferences was found to be the polarization interaction. Notably, although the electrostatic interaction dominates the binding energies, it contributes favorably to the relative stabilities of most isomers. The Li+···Li+ distance is identified as the key geometrical parameter that correlates with the strength of the polarization of the Sx2- fragment imposed by the Li+ cations. Further BLW-ED analyses reveal that the cooperativity of the Li+ cations primarily determines the relative strength of the polarization.

Metastability of π-π stacking between the closed-shell ions of like charges

Planar cations or anions can form stacks in crystals or solutions, where the surrounding or environment plays a decisive role as demonstrated in previous studies. However, it remains unclear whether these counterintuitive interactions possess any inherent stability or are thoroughly repulsive if the constraint of environment is removed. In this work, we explored the inherent stability of π-π stacking between closed-shell ions of like charges with prototypes derived from experimental studies. The inherent metastability was identified by the characteristic local minima and the transition states preventing their dissociation and verified by ab initio molecular dynamics (AIMD) simulations. The nature of involved interactions was deciphered with the energy decomposition approach based on the block-localized wavefunction method (BLW-ED). Like the conventional neutral π-π stacking interactions, electron correlation is the most attractive energy component. But it is overturned by the Coulombic repulsion between net charges for all modes of dimerization, resulting in the overall repulsive inter-cation or anion interactions. Contributions from van der Waals interactions were also observed in the reduced density gradient analysis. The origin of the metastability was elucidated by examining the contributions of individual physical factors to the well-depths. The inherent metastability originates from the electron correlation, which dramatically increases due to the enhanced overlap between ions from a transition state to its corresponding minimum.

Planar Four-Membered Diboron Actinide Compound with Double Möbius Aromaticity.

The Möbius rule predicts that a planar four-membered metallacycle can be aromatic with four mobile electrons, but such a simple ring has escaped recognition because it usually favors Hückel anti-aromaticity. Here, we report that a quasi-square four-membered actinide compound (Pa2B2) is doubly Möbius aromatic. Chemical bonding analyses reveal that this diboron protactinium molecule has four delocalized π electrons in addition to four delocalized σ electrons, satisfying the 4n Möbius rule for both σ and π components. Energetically, the block-localized wavefunction method, which is the simplest variant of ab initio valence bond theory, shows that the delocalization energy for the π and σ electrons reaches up to 65.0 and 72.3 kcal/mol, respectively, while the extra cyclic resonance energy (ECRE) amounts to 45 kcal/mol. The large positive ECRE values strongly confirm the unprecedented double Möbius aromaticity in Pa2B2. We anticipate that this new type of aromatic molecule can enrich the concept of Möbius aromaticity and open a new avenue for actinide compounds.

Selectivity Rule of Cryptands for Anions: Molecular Rigidity and Bonding Site.

Cryptands utilize inside CH or NH groups as hydrogen bond (H-bond) donors to capture anions such as halides. In this work, the nature and selectivity of confined hydrogen bonds inside cryptands were computationally analyzed with the energy decomposition scheme based on the block-localized wavefunction method (BLW-ED), aiming at an elucidation of governing factors in the binding between cryptands and anions. It was revealed that the intrinsic strengths of inward hydrogen bonds are dominated by the electrostatic attraction, while the anion preferences (selectivity) of inner CH and NH hydrogen bonds are governed by the Pauli exchange repulsion and electrostatic interaction, respectively. Typical conformers of cages are classified into two groups, including the C3(h) -symmetrical conformers, in which all halide anions are located near the centroids of cages, and the "semi-open" conformers, which exhibit shifted bonding sites for different halide anions. Accordingly, the difference in governing factors of selectivity is attributed to either the rigidity of cages or the binding site of anions for these two groups. In details, the C3 conformers of NH cryptands can be enlarged more remarkably than the C3(h) -symmetrical conformers of CH cryptands as the size of anion (ionic radius) increases, resulting in the relaxation of the Pauli repulsion and a dramatic reduction in electrostatic attraction, which eventually rules the selectivity of NH cryptands for halide anions. By contrary, the CH cryptands are more rigid and cannot effectively reduce the Pauli repulsion, which subsequently governs the anion preference. Unlike C3 conformers whose rigidity determines the selectivity, semi-open conformers exhibit different binding sites for different anions. From F- to I- , the bonding site shifts toward the outside end of the pocket inside the semi-open NH cryptand, leading to the significant reduction of the electrostatic interaction that dominates the anion preference. Differently, binding sites are much less affected by the size of anion inside the semi-open CH cryptand, in which the Pauli exchange repulsion remains the key factor for the selectivity of inner hydrogen bonds.

On the nature of inter-anion coinage bonds.

To explore the binding energy profiles and elucidate the bonding nature in counter-intuitive anion⋯anion coinage bonds (CiBs), thirty-one complexes were constructed, and the inter-anion CiBs were studied theoretically. The metastability was evidenced by the characteristic potential wells in six cases, demonstrating that anions [Au(CN)4]-, [Ag(CN)2]- and [AuO]- are appropriate building blocks for CiBs. The kinetic stability was further supported by ab initio molecular dynamics (AIMD) simulations and the analyses based on the local vibrational mode and quantum theory of atoms in molecules (QTAIM) methods. The anion⋯anion CiBs in the dimers of [AuCl4]- and [Au(CN)4]- previously observed in condensed phases were confirmed to be thoroughly repulsive under vacuum, but turned attractive in the crystal environment which was simulated using the solvation model based on density (SMD). However, the intrinsic strength of the inter-anion bonding is barely variated by the environment, as it is the combination of the inter-anion interaction and the environment effect that stabilizes the anion pairs. The block-localized wavefunction (BLW) method and its corresponding energy decomposition (BLW-ED) approach were further employed aiming at a chemically meaningful explanation for these counterintuitive phenomena. By inspecting the profiles of energy components, we identified the vital distinction between inter-anion CiBs and conventional non-covalent interactions lying in the electrostatic interaction, which variates nonmonotonically in the inter-anion complexes. The electrostatic interaction also dominates the depth of potential wells, which is commonly used to evaluate the kinetic stability, while Pauli exchange repulsion is the most repulsive factor preventing the formation of anion adducts. The importance of the Pauli exchange repulsion was further highlighted by comparing cases with and without metastability, in which the absence of a potential well is solely caused by the enhancement of the Pauli exchange repulsion.

The roles of charge transfer and polarization in non-covalent interactions: a perspective from ab initio valence bond methods.

Noncovalent interactions are ubiquitous and have been well recognized in chemistry, biology and material science. Yet, there are still recurring controversies over their natures, due to the wide range of noncovalent interaction terms. In this Essay, we employed the Valence Bond (VB) methods to address two types of interactions which recently have drawn intensive attention, i.e., the halogen bonding and the CH‧‧‧HC dihydrogen bonding. The VB methods have the advantage of interpreting molecular structures and properties in the term of electron-localized Lewis (resonance) states (structures), which thereby shed specific light on the alteration of the bonding patterns. Due to the electron localization nature of Lewis states, it is possible to define individually and measure both polarization and charge transfer effects which have different physical origins. We demonstrated that both the ab initio VB method and the block-localized wavefunction (BLW) method can provide consistent pictures for halogen bonding systems, where strong Lewis bases NH3, H2O and NMe3 partake as the halogen bond acceptors, and the halogen bond donors include dihalogen molecules and XNO2 (X = Cl, Br, I). Based on the structural, spectral, and energetic changes, we confirm the remarkable roles of charge transfer in these halogen bonding complexes. Although the weak C-H∙∙∙H-C interactions in alkane dimers and graphene sheets are thought to involve dispersion only, we show that this term embeds delicate yet important charge transfer, bond reorganization and polarization interactions.

On the Bonding Nature in the Crystalline Tri-Thorium Cluster: Core-Shell Syngenetic σ-Aromaticity.

A unique thorium‐thorium bond was observed in the crystalline tri‐thorium cluster [{Th(η8‐C8H8)(μ3‐Cl)2}3{K(THF)2}2]∞, though the claim of σ‐aromaticity for Th3 bond has been questioned. Herein, a new type of core–shell syngenetic bonding model is proposed to describe the stability of this tri‐thorium cluster. The model involves a 3c–2e bond in the Th3 core and a multicentered (ThCl2)3 charge‐shift bond with 12 electrons scattering along the outer shell. To differentiate the strengths of the 3c–2e bond and the charge‐shift bond, the block‐localized wavefunction (BLW) method which falls into the ab initio valence bond (VB) theory is employed to construct a strictly core/shell localized state and its contributing covalent resonance structure for the Th3 core bond. By comparing with the σ‐aromatic H3 + and nonaromatic Li3 +, the computed resonance energies and extra cyclic resonance energies confirm that this Th3 core bond is truly delocalized and σ‐aromatic.

On the Bonding Nature in the Crystalline Tri‐Thorium Cluster: Core‐Shell Syngenetic σ‐Aromaticity

AbstractA unique thorium‐thorium bond was observed in the crystalline tri‐thorium cluster [{Th(η8‐C8H8)(μ3‐Cl)2}3{K(THF)2}2]∞, though the claim of σ‐aromaticity for Th3 bond has been questioned. Herein, a new type of core–shell syngenetic bonding model is proposed to describe the stability of this tri‐thorium cluster. The model involves a 3c–2e bond in the Th3 core and a multicentered (ThCl2)3 charge‐shift bond with 12 electrons scattering along the outer shell. To differentiate the strengths of the 3c–2e bond and the charge‐shift bond, the block‐localized wavefunction (BLW) method which falls into the ab initio valence bond (VB) theory is employed to construct a strictly core/shell localized state and its contributing covalent resonance structure for the Th3 core bond. By comparing with the σ‐aromatic H3+ and nonaromatic Li3+, the computed resonance energies and extra cyclic resonance energies confirm that this Th3 core bond is truly delocalized and σ‐aromatic.

The strength and selectivity of perfluorinated nano-hoops and buckybowls for anion binding and the nature of anion-π interactions.

Perfluorinated cycloparaphenylenes (F-[n]CPP, n= 5-8), boron nitride nanohoop (F-[5]BNNH), and buckybowls (F-BBs) were proposed as anion receptors via anion-π interactions with halide anions (Cl- , Br- and I- ), and remarkable binding strengths up to -294.8kJ/mol were computationally verified. The energy decomposition approach based on the block-localized wavefunction method, which combines the computational efficiency of molecular orbital theory and the chemical intuition of ab initio valence bond theory, was applied to the above anion-π complexes, in order to elucidate the nature and selectivity of these interactions. The overall attraction is mainly governed by the frozen energy component, in which the electrostatic interaction is included. Remarkable binding strengths with F-[n]CPPs can be attributed to the accumulated anion-π interactions between the anion and each conjugated ring on the hoop, while for F-BBs, additional stability results from the curved frameworks, which distribute electron densities unequally on π-faces. Interestingly, the strongest host was proved to be the F-[5]BNNH, which exhibits the most significant anisotropy of the electrostatic potential surface due to the difference in the electronegativities of nitrogen and boron. The selectivity of each host for anions was explored and the importance of the often-overlooked Pauli exchange repulsion was illustrated. Chloride anion turns out to be the most favorable anion for all receptors, due to the smallest ionic radius and the weakest destabilizing Pauli exchange repulsion.

Inter-anion chalcogen bonds: Are they anti-electrostatic in nature?

Inter-anion hydrogen and halogen bonds have emerged as counterintuitive linkers and inspired us to expand the range of this unconventional bonding pattern. Here, the inter-anion chalcogen bond (IAChB) was proposed and theoretically analyzed in a series of complexes formed by negatively charged bidentate chalcogen bond donors with chloride anions. The kinetic stability of IAChB was evidenced by the minima on binding energy profiles and further supported by ab initio molecular dynamic simulations. The block-localized wave function (BLW) method and its subsequent energy decomposition (BLW-ED) approach were employed to elucidate the physical origin of IAChB. While all other energy components vary monotonically as anions get together, the electrostatic interaction behaves exceptionally as it experiences a Coulombic repulsion barrier. Before reaching the barrier, the electrostatic repulsion increases with the shortening Ch⋯Cl- distance as expected from classical electrostatics. However, after passing the barrier, the electrostatic repulsion decreases with the Ch⋯Cl- distance shortening and subsequently turns into the most favorable trend among all energy terms at short ranges, representing a dominating force for the kinetic stability of inter-anions. For comparison, all energy components exhibit the same trends and vary monotonically in the conventional counterparts where donors are neutral. By comparing inter-anions and their conventional counterparts, we found that only the electrostatic energy term is affected by the extra negative charge. Remarkably, the distinctive (nonmonotonic) electrostatic energy profiles were reproduced using quantum mechanical-based atomic multipoles, suggesting that the crucial electrostatic interaction in IAChB can be rationalized within the classical electrostatic theory just like conventional non-covalent interactions.

Classical Electrostatics Remains the Driving Force for Interanion Hydrogen and Halogen Bonding.

Interanion hydrogen bonding (IAHB) and halogen bonding (IAXB) have emerged as a counterintuitive linker in a range of fascinating applications. Despite the overall repulsive (positive) binding energy, anions are trapped in a local minimum with its corresponding transition state (TS) preventing dissociation. In other words, the adduct of anions is metastable. Seemingly, the electrostatic paradigm and force field description of hydrogen/halogen bonding (HB/XB) are challenged, because of the preconceived Coulombic repulsion. Aiming at an insightful understanding of these interanion phenomena, we employed the energy decomposition approach based on the block-localized wavefunction method (BLW-ED) to investigate a series of exemplary interanion complexes. As expected, the key distinction from the conventional HB/XB lies in the electrostatic interaction, which is not increasingly repulsive as anions gradually approach to each other. Rather, there is a Coulombic barrier at a certain point. After this point, the electrostatic repulsion diminishes with the decreasing distance between anions. Differently, other energy components vary monotonically just like in conventional cases. The nonmonotonic characteristic of the electrostatic interaction in interanion complexes was reproduced using the multipole expansion in AMOEBA polarizable force field in which the state-specified atomic multipoles were adopted. This suggests that the nonmonotonicity can be well interpreted by classical electrostatic theory and there is no conceptual difference between conventional HB/XB and IAHB/IAXB. The stability of IAHB/IAXB depends on the competition between the local attractive HB/XB and the global Coulombic repulsion of net charges, though there is cooperativity between these two contrasting forces. This concise model was supported by the attractive IAHB/IAXB in modified molecular capsules, which exhibit strong quadruple HB/XBs and a considerable distance between charged substituents.

Hydrogen and Halogen Bonding in Homogeneous External Electric Fields: Modulating the Bond Strengths.

Recent years have witnessed various fascinating phenomena arising from the interactions of noncovalent bonds with homogeneous external electric fields (EEFs). Here we performed a computational study to interpret the sensitivity of intrinsic bond strengths to EEFs in terms of steric effect and orbital interactions. The block-localized wavefunction (BLW) method, which combines the advantages of both ab initio valence bond (VB) theory and molecular orbital (MO) theory, and the subsequent energy decomposition (BLW-ED) approach were adopted. The sensitivity was monitored and analyzed using the induced energy term, which is the variation in each energy component along the EEF strength. Systems with single or multiple hydrogen (H) or halogen (X) bond(s) were also examined. It was found that the X-bond strength change to EEFs mainly stems from the covalency change, while generally the steric effect rules the response of H-bonds to EEFs. Furthermore, X-bonds are more sensitive to EEFs, with the key difference between H- and X-bonds lying in the charge transfer interaction. Since phenylboronic acid has been experimentally used as a smart linker in EEFs, switchable sensitivity was scrutinized with the example of the phenylboronic acid dimer, which exhibits two conformations with either antiparallel or parallel H-bonds, thereby, opposite or consistent responses to EEFs. Among the studied systems, the quadruple X-bonds in molecular capsules exhibit remarkable sensitivity, with its interaction energy increased by -95.2 kJ mol-1 at the EEF strength 0.005 a.u.

Anti-Electrostatic Main Group Metal–Metal Bonds That Activate CO2

There has been growing interest in the CO2 capture and reduction by transition-metal-free catalysts. Here we performed a proof-of-concept study using an ab initio valence bond method called the block-localized wave function (BLW) method. The integrated BLW and density function theory (DFT) computations demonstrated that heterobimetallic Ae+/Al(I) (Ae represents alkaline earth metals Mg and Ca) Lewis acid/base combinations without transition metals can facilely capture and activate CO2. There are two remarkable findings in this study. The first concerns the ionic nature of the metal-metal bonds. The experimentally synthesized low valent aluminum compound with a bidentate β-diketiminate (BDI) ligand, or (BDI)Al(I) in brief, is a Lewis base due to the lone pair on the aluminum cation though overall Al(I) is positively charged. Al(I) can form ionic metal-metal bonds with the alkaline earth metals of the positively charged Lewis acids (BDI)Ae+. This type of ionic metal-metal bonds is counterintuitive and antielectrostatic as both metals carry positive charges. The second finding is the CO2 activation mechanism. (BDI)Al(I) can effectively bind and activate CO2 by transferring one electron to CO2, and the resulting complex can be best expressed as [(BDI)Al(I)]+[CO2]-. The participation of (BDI)Ae+ further enhances the capture and activation of CO2 by (BDI)Al(I).

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF.

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to a relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time-reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break the symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time-reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

Role of Charge Transfer in Halogen Bonding.

Halogen bonding has received intensive attention recently for its applications in the construction of supramolecular assemblies and crystal engineering and its implications and potentials in chemical and biological processes and rational drug design. Peculiarly, in intermolecular interactions, halogen atoms are known as electron-donating groups carrying partial negative charges in molecules due to its high electronegativity, but they can counterintuitively act as Lewis acids and bind with Lewis bases in the form of a halogen bond. The unsettling issue regarding the nature of the halogen bonding is whether the electrostatics or charge transfer interaction dominates. The recently proposed σ-hole concept nicely reinforces the role of electrostatic attraction. Also, good correlations between the halogen bonding strength and the interaction energy from the simple point-charge model have been found. This leads to the claim that there is no need to invoke the charge transfer concept in the halogen bond. But there is alternative evidence supporting the importance of charge transfer interaction. Here, we visited a series of prominent halogen bonded complexes of the types Y3C-X···Z (X = Br, I; Y = F, Cl, Br; Z = F-, Cl-, Br-, I-, NMe3) with the block-localized wave function (BLW) method at the M06-2X-D3/6-311+G(d,p) (def2-SVP for iodine) level of theory. As the simplest variant of ab initio valence bond (VB) theory, the BLW method is unique in the strict localization of electrons within interacting moieties, allowing for quantitative evaluation of the charge transfer effect on geometries, spectral properties, and energetics in halogen bonding complexes. By comparing the halogen bonding complexes with and without the charge transfer interaction, we proved that the charge transfer interaction significantly shortens the X···Z bonding distance and stretches the C-X bonds. But the shortening of the halogen bonding results in the less favorable steric effect, which is composed of Pauli repulsion, electrostatics, and electron correlation. There are approximate linear correlations between the charge transfer effect and binding energy and between bonding distance and binding energy. These correlations may lead to the illusion that the charge transfer interaction is unimportant or irrelevant, but further analyses showed that the inclusion of charge transfer is critical for the proper description of the halogen bonding, as considering only electrostatics and polarization leads to only about 45-60% of the binding strengths and much elongated bonding distances.

Generalization of Block-Localized Wave Function for Constrained Optimization of Excited Determinants.

The block-localized wave function method is useful to provide insights on chemical bonding and intermolecular interactions through energy decomposition analysis. The method relies on block localization of molecular orbitals (MOs) by constraining the orbitals to basis functions within given blocks. Here, a generalized block-localized orbital (GBLO) method is described to allow both physically localized and delocalized MOs to be constrained in orbital-block definitions. Consequently, GBLO optimization can be conveniently tailored by imposing specific constraints. The GBLO method is illustrated by three examples: (1) constrained polarization response orbitals through dipole and quadrupole perturbation in a water dimer complex, (2) the ground and first excited-state potential energy curves of ethene about its C-C bond rotation, and (3) excitation energies of double electron excited states. Multistate density functional theory is used to determine the energies of the adiabatic ground and excited states using a minimal active space (MAS) comprising specifically charge-constrained and excited determinant configurations that are variationally optimized by the GBLO method. We find that the GBLO expansion that includes delocalized MOs in configurational blocks significantly reduces computational errors in comparison with physical block localization, and the computed ground- and excited-state energies are in good accordance with experiments and results obtained from multireference configuration interaction calculations.

Resonance-Assisted but Antielectrostatic Intramolecular Au···H-O Hydrogen Bonding in Gold(I) Complexes: A Computational Verification.

The peculiar Au···H hydrogen bonding has garnered significant interests, but its existence and nature remain an open question in gold chemistry. In this paper, we established and authenticated the first intramolecular Au···H-O hydrogen bonding in gold(I) complexes. Our computational results based on an ab initio valence bond (VB) method, namely the block-localized wave function (BLW) method, clarified that this Au···H-O hydrogen bonding is a new type of resonance-assisted hydrogen bond (RAHB). In this RAHB, π conjugation in the hydrogen bond donor has the capability to modulate the Au···H-O hydrogen bonding. In contrast to conventional RAHBs, however, the Au···H-O hydrogen bonding is antielectrostatic and is dominated by the charge transfer interaction including cooperative σ electron donation from the d orbital of the metal center to the antibonding orbital of the hydroxyl group and π conjugation from the hydroxyl group to adjacent π deficient or strong π electron withdrawing groups. This novel theoretical perspective not only confirms the existence and reveals the nature of intramolecular Au···H hydrogen bonding but also provides a promising strategy to rationally design strong gold hydrogen bonds.

External Electric Fields Interrupt the Concerted Cope Rearrangement of Semibullvalene.

The topic of this paper is whether the mechanism of the degenerate Cope rearrangement of semibullvalene can be affected by the presence of electrostatic fields. Herein, we report that the shape of the energy surface, as demonstrated by an "interrupted" (stepwise) mechanism, is altered in the presence of a copper cation, Cu+. Natural bond-orbital and block-localized wave-function energy decomposition analyses suggest that orbital and electrostatic interactions play a major role in altering the shape of the energy surface. Applying additional external electric fields (EEFs) induces a significant change to the energy surface with Cu+ present but negligible effects in the absence of Cu+. These findings are consistent with recent studies that demonstrate that EEFs more readily stabilize/destabilize systems with larger, more polarizable, dipole moments.

Water adlayers on noble metal surfaces: Insights from energy decomposition analysis

Water molecules adsorbed on noble metal surfaces are of fundamental interest in surface science, in heterogeneous catalysis, and as a model for the metal/water interface. Herein, we analyze 28 water structures adsorbed on five noble metal surfaces (Cu, Ag, Au, Pd, and Pt) via density functional theory and energy decomposition analysis based on the block localized wave function technique. Structures, ranging from monomers to ice adlayers, reveal that the charge transfer from water to the surface is nearly independent from the charge transfer between the water molecules, while the polarization energies are cooperative. Dense water-water networks with small surface dipoles, such as the 39×39 unit cell [experimentally observed on Pt(111)], are favored compared to the highly ordered and popular Hup and Hdown phases. The second main result of our study is that the many-body interactions, which stabilize the water assemblies on the metal surfaces, are dominated by the polarization energies, with the charge transfer scaling with the polarization energies. Hence, if an empirical model could be found that reproduces the polarization energies, the charge transfer could be predicted as well, opening exciting perspectives for force field development.