van der Waals corrected density functional calculations of the adsorption of benzene on the Cu (111) surface

We investigate the performance of several van der Waals (vdW) functionals at calculating the interactions between benzene and the copper (111) surface, using the local orbital approach in the SIESTA code. We demonstrate the importance of using surface optimized basis sets to calculate properties of pure surfaces, including surface energies and the work function. We quantify the errors created using (3 × 3) supercells to study adsorbate interactions using much larger supercells, and show non‐negligible errors in the binding energies and separation distances. We examine the eight high‐symmetry orientations of benzene on the Cu (111) surface, reporting the binding energies, separation distance, and change in work function. The optimized vdW‐DF(optB88‐vdW) functional provides superior results to the vdW‐DF(revPBE) and vdW‐DF2(rPW86) functionals, and closely matches the experimental and experimentally deduced values. This work demonstrates that local orbital methods using appropriate basis sets combined with a vdW functional can model adsorption between metal surfaces and organic molecules.     

First principles study of gallium‐cleaning for hydrogen‐contaminated α‐Al2O3(0001) surfaces

The use of gallium for cleaning hydrogen‐contaminated Al₂O₃ surfaces is explored by performing first principles density functional calculations of gallium adsorption on a hydrogen‐contaminated Al‐terminated α‐Al₂O₃(0001) surface. Both physisorbed and chemisorbed H‐contaminated α‐Al₂O₃(0001) surfaces with one monolayer (ML) gallium coverage are investigated. The thermodynamics of gallium cleaning are considered for a variety of different asymptotic products, and are found to be favorable in all cases. Physisorbed H atoms have very weak interactions with the Al₂O₃ surface and can be removed easily by the Ga ML. Chemisorbed H atoms form stronger interactions with the surface Al atoms. Bonding energy analysis and departure simulations indicate, however, that chemisorbed H atoms can be effectively removed by the Ga ML. © 2013 Wiley Periodicals, Inc.     

Dichloro[N,N‐dimethyl‐1‐(1‐methyl‐1H‐tetrazol‐5‐yl‐κN4)methanamine‐κN]copper(II)

The title compound, [CuCl₂(C₅H₁₁N₅)], is the first structurally characterized molecular chelate complex involving an α‐­amino­alkyl­tetrazole. There are two complex mol­ecules in the asymmetric unit. The ligand mol­ecules are bidentate. Both Cu atoms reveal rather distorted square‐planar coordinations. The complex mol­ecules are linked together by van der Waals interactions only.     

Removal of 4‐chlorophenol using graphene,graphene oxide,and a‐doped graphene (A = N,B): A computational study

Density functional theory (including van der Waals correction with the PBE‐D functional) is applied to the study of 4‐chlorophenol (4‐CP) adsorption on graphene oxide (GO), A‐doped graphene (A = N, B), and pristine graphene and test their possible application for 4‐CP removal. Results show that on GO adsorption is improved by the hydrogen bond interactions between the adsorbents and 4‐CP, suggesting that functionalized graphene is a preferable alternative than pristine graphene for 4‐CP removal. In addition, the stability of hydrogen bonds is confirmed by molecular dynamics calculations using the PM6 potential. Without hydrogen bonds, A‐doped graphene models show a comparable performance for 4‐CP removal than pristine graphene. Finally, even in a solvent medium, 4‐CP adsorption is strong. © 2013 Wiley Periodicals, Inc.     

On the van der Waals interactions and the stability of polypeptide chains in helical conformations

In this work, we aim to investigate the contribution of van der Waals (vdW) interactions to the stability of polypeptides in helical conformations studying infinitely long chains of alanine and glycine with density functional theory. To account for vdW interactions, we have used the interatomic pairwise dispersion approach proposed by Tkatchenko–Scheffler (TS), the TS approach with self‐consistent screening (SCS) that self‐consistently includes long‐range electrostatic effects (TS + SCS), the D2 and D3 methods of Grimme et al., and the Langreth–Lundqvist procedure that treats nonlocally the correlation part of the approximation to the exchange‐correlation (xc) functional (called DF). First, we have tested the performance of these strategies studying a set of representative hydrogen bonded dimers. Next, we have studied polyalanine and polyglicine in π‐helix, α‐helix, ‐helix, 2₇, and polyproline‐II conformations and in a fully extended structure. We have found that the DF methodology in combination with a modified version for the Becke approximation to the exchange (optB86b), the D2, D3, TS, and TS + SCS strategies in combination with the Perdew–Burke–Ernzerhof approximation to the xc functional, describe fairly well dimer association energies. Furthermore, the DF method and the D2, D3, TS, and TS + SCS strategies predict very similar helical stabilities even though the approximation used in DF for describing the long‐range dispersion interactions is different that the one used in D2/D3 and TS/TS + SCS. We found that the stability doubles for π and α helices if vdW interactions are taken into account. © 2015 Wiley Periodicals, Inc.     

A comparative study of two polymorphs of bis(1‐hydroxy‐2‐methylpropan‐2‐aminium) carbonate

Alkanolamines have been known for their high CO₂ absorption for over 60 years and are used widely in the natural gas industry for reversible CO₂ capture. In an attempt to crystallize a salt of (RS)‐2‐(3‐benzoylphenyl)propionic acid with 2‐amino‐2‐methylpropan‐1‐ol, we obtained instead a polymorph (denoted polymorph II) of bis(1‐hydroxy‐2‐methylpropan‐2‐aminium) carbonate, 2C₄H₁₂NO⁺·CO₃²⁻, (I), suggesting that the amine group of the former compound captured CO₂ from the atmosphere forming the aminium carbonate salt. This new polymorph was characterized by single‐crystal X‐ray diffraction analysis at low temperature (100 K). The salt crystallizes in the monoclinic system (space group C2/c, Z = 4), while a previously reported form of the same salt (denoted polymorph I) crystallizes in the triclinic system (space group P, Z = 2) [Barzagli et al. (2012). ChemSusChem, 5 , 1724–1731]. The asymmetric unit of polymorph II contains one 1‐hydroxy‐2‐methylpropan‐2‐aminium cation and half a carbonate anion, located on a twofold axis, while the asymmetric unit of polymorph I contains two cations and one anion. These polymorphs exhibit similar structural features in their three‐dimensional packing. Indeed, similar layers of an alternating cation–anion–cation neutral structure are observed in their molecular arrangements. Within each layer, carbonate anions and 1‐hydroxy‐2‐methylpropan‐2‐aminium cations form planes bound to each other through N—H…O and O—H…O hydrogen bonds. In both polymorphs, the layers are linked to each other via van der Waals interactions and C—H…O contacts. In polymorph II, a highly directional C—H…O contact (C—H…O = 156°) shows as a hydrogen‐bonding interaction. Periodic theoretical density functional theory (DFT) calculations indicate that both polymorphs present very similar stabilities.     

(S)‐3‐(Ammoniomethyl)‐5‐methylhexanoate (pregabalin)

The title compound, C₈H₁₇NO₂, exists as a zwitterion, adopting a propeller conformation. Molecules self‐assemble to form a hydrogen‐bonded layer parallel to the ab crystallographic plane connected by N⁺—H...O⁻ and C—H...O⁻ hydrogen bonds. These layers are stacked along the c axis and are stabilized by van der Waals interactions.     

First‐principles calculations of oxygen adsorption on the Ti3Al (0001) surface

Adsorption energies and density of states for O atoms adsorption on the Ti₃Al (0001) surface have been calculated using first‐principles calculations based on density functional theory. It is found that the order of O atom adsorption on the Ti₃Al (0001) surface is associated with the adsorption energy as well as the distance of O atoms because of the interaction. The adsorption energy mainly depends on the bond number and bond strength between O and Ti atoms, and the adsorption site with rich‐Ti surface (H_I and HCP_Al) is first priority. The adsorption energy decreases with the increase of the oxygen coverage because of the characteristics of the valence d‐orbitals of transition metals surface. Furthermore, the density of states indicates that the hybridization peak of O and Ti atoms is mainly from the contribution of Ti 3d‐ and O 2p‐orbitals, and the hybridization peak of O and Al atoms from the contribution of Al 2p‐ and O 2p‐orbitals. Copyright © 2016 John Wiley & Sons, Ltd.     

(4‐Aminopyridine‐κN1)(phthalocyaninato‐κ4N)zinc(II) tetrahydrofuran disolvate

The title compound, [Zn(C₃₂H₁₆N₈)(C₅H₆N₂)]·2C₄H₈O, consists of one (phthalocyaninato)zinc (ZnPc) unit, a coordinated 4‐aminopyridine (4‐ap) molecule and two tetrahydrofuran (THF) solvent molecules. The central Zn atom is (4+1)‐coordinated by four isoindole N atoms of the Pc core and by the pyridine N atom of 4‐aminopyridine. The Zn atom is displaced by 0.4464 (8) Å from the isoindole N₄ plane towards the pyridine N atom. The crystal structure is stabilized by intermolecular amine–phthalocyaninate N—H...N hydrogen bonds and π–π interactions between the aggregated Pc rings, which form molecular layers, and by weak van der Waals interactions between the layers. As well as hindering the aggregation of ZnPc molecules by occupying an axial position, the amino group will add new interactions which will favor applications of ZnPc, for example, as a sensitizer of photodynamic therapy.     

1‐(4‐Bromo­phenyl)‐3‐(4‐methyl­phenyl)­prop‐2‐en‐1‐one

In the mol­ecule of the title compound, C₁₆H₁₃BrO, the two benzene rings are rotated in opposite directions with respect to the central C—C=C—C part of the mol­ecule. The phenone O atom deviates from the least‐squares plane of the mol­ecule by 0.300 (3) Å. In the crystal structure, mol­ecules are paired through C—H⋯π interactions. The molecular pairs along [001] are hydrogen bonded through three translation‐related co‐operative hydrogen bonds in the `bay area', forming molecular chains, which are further hydrogen bonded through C—H⋯Br weak interactions, forming (010) molecular layers. In the third direction, there are only weak van der Waals interactions. The co‐operative hydrogen bonds in the `bay area' are discussed briefly.     

van der Waals interactions are critical in Car–Parrinello molecular dynamics simulations of porphyrin–fullerene dyads

The interplay between electrostatic and van der Waals (vdW) interactions in porphyrin‐C₆₀ dyads is still under debate despite its importance in influencing the structural characteristics of such complexes considered for various applications in molecular photovoltaics. In this article, we sample the conformational space of a porphyrin‐C₆₀ dyad using Car–Parrinello molecular dynamics simulations with and without empirical vdW corrections. Long‐range vdW interactions, which are poorly described by the commonly used density functional theory functionals, prove to be essential for a proper dynamics of the dyad moieties. Inclusion of vdW corrections brings porphyrin and C₆₀ close together in an orientation that is in agreement with experimental observations. The structural differences arising from the vdW corrections are shown to be significant for several properties and potentially less important for others. Additionally, our Mulliken population analysis reveals that contrary to the common belief, porphyrin is not the primary electron donating moiety for C₆₀. In the considered dyad, fullerene's affinity for electrons is primarily satisfied by charge transfer from the amide group of the linker. However, we show that in the absence of another suitable bound donor, C₆₀ can withdraw electrons from porphyrin if it is sufficiently close. © 2015 Wiley Periodicals, Inc.     

Formation of self-assembled chains of tetrathiafulvalene on a Cu(100) surface

Formation of self-assembled chains of tetrathiafulvalene (TTF) on the Cu(100) surface has been investigated by scanning tunneling microscopy and density functional theory calculations that include semiempirical van der Waals (vdW) interaction corrections. The calculations show that the chain structures observed in the experiments can only be explained by including the vdW interactions. The molecules are tilted along the chain in order to achieve maximal intermolecular interaction. The chains are metastable on the surface, which is consistent with the experimental observation that they disappear after annealing. The fact that all TTF chains observed in the experiment are short might be possibly explained by the interplay between the stabilizing vdW molecule-molecule interaction and the destabilizing rearrangement of surface atoms due to the strong molecule-substrate interaction.     

(E)‐2‐[(4‐Chlorophenyl)iminomethyl]‐5‐methoxyphenol and (E)‐2‐[(2‐chlorophenyl)iminomethyl]‐5‐methoxyphenol: X‐ray and DFT‐calculated structures

The crystal structures of the title 4‐chlorophenyl, (I), and 2‐chlorophenyl, (II), compounds, both C₁₄H₁₂ClNO₂, have been determined using X‐ray diffraction techniques and the molecular structures have also been optimized at the B3LYP/6‐31 G(d,p) level using density functional theory (DFT). The X‐ray study shows that the title compounds both have strong intramolecular O—H...N hydrogen bonds and that the crystal networks are primarily determined by weak C—H...π and van der Waals interactions. The strong intramolecular O—H...N hydrogen bond is evidence of the preference for the phenol–imine tautomeric form in the solid state. The IR spectra of the compounds were recorded experimentally and also calculated for comparison. The results from both the experiment and theoretical calculations are compared in this study.     

Halogen‐bonded adduct of 1,2‐dibromo‐1,1,2,2‐tetrafluoroethane and 1,4‐diazabicyclo[2.2.2]octane

Halogen bonding is an intermolecular interaction capable of being used to direct extended structures. Typical halogen‐bonding systems involve a noncovalent interaction between a Lewis base, such as an amine, as an acceptor and a halogen atom of a halofluorocarbon as a donor. Vapour‐phase diffusion of 1,4‐diazabicyclo[2.2.2]octane (DABCO) with 1,2‐dibromotetrafluoroethane results in crystals of the 1:1 adduct, C₂Br₂F₄·C₆H₁₂N₂, which crystallizes as an infinite one‐dimensional polymeric structure linked by intermolecular N...Br halogen bonds [2.829 (3) Å], which are 0.57 Å shorter than the sum of the van der Waals radii.     

cis‐[(4R,5R)‐4,5‐Bis­(amino­methyl)‐2,2‐di­methyl‐1,3‐dioxolane‐N,N′](malonato‐O,O′)­platinum(II), an anticancer agent

In the title compound, [Pt(C₃H₂O₄)(C₇H₁₆N₂O₂)], the Pt atom is coordinated to two O and two N atoms in a square‐planar arrangement. The two independent mol­ecules, which have very similar structures, are approximately related by pseudo‐twofold screw‐axis symmetry. The six‐membered chelate ring in the leaving ligand assumes a conformation intermediate between the half‐chair and boat forms. The seven‐membered ring in the carrier ligand assumes a twist‐chair conformation and the oxolane ring assumes an envelope conformation. The crystal packing consists of extensive hydrogen‐bonding networks which form two‐dimensional molecular layers, and there are weak van der Waals interactions between these layers.     

Short but Weak: The Z‐DNA Lone‐Pair⋅⋅⋅π Conundrum Challenges Standard Carbon Van der Waals Radii

Current interest in lone‐pair???π (lp???π) interactions is gaining momentum in biochemistry and (supramolecular) chemistry. However, the physicochemical origin of the exceptionally short (ca. 2.8 Å) oxygen‐to‐nucleobase plane distances observed in prototypical Z‐DNA CpG steps remains unclear. High‐level quantum mechanical calculations, including SAPT2+3 interaction energy decompositions, demonstrate that lp???π contacts do not result from n→π* orbital overlaps but from weak dispersion and electrostatic interactions combined with stereochemical effects imposed by the locally strained structural context. They also suggest that the carbon van der Waals (vdW) radii, originally derived for sp³ carbons, should not be used for smaller sp² carbons attached to electron‐withdrawing groups. Using a more adapted carbon vdW radius results in these lp???π contacts being no longer of the sub‐vdW type. These findings challenge the whole lp???π concept that refers to elusive orbital interactions that fail to explain short interatomic contact distances.     

Adsorption and dissociation of the methanethiol (CH3SH) molecule on the Fe(100) surface

These contributions explore interaction modes between the methanethoil (CH₃SH) molecule and the Fe(100) surface via implementing accurate density functional theory (DFT) calculations with the inclusion of van der Waals corrections. We consider three adsorption sites over the Fe(100) surface, namely, top(T), bridge (B), and hollow (H) sites as potential catalytic active sites for the molecular and dissociative adsorption of the CH₃SH molecule. The molecular adsorption structures are found to occupy either B or T sites with former sites holding higher stability by 0.17 eV. The inclusion of van der Waals corrections refound to slightly alter adsorption energies. For instance, adsorption energies increased by ~ 0.18 and ~ 0.21 eV for B and T structure, respectively, in reference to values obtained by the plain generalized gradient approximation (GGA) functional. A stability ordering of the dissociation products was found to follow the sequence (CH₄, S) > (CH₃, S, H) > (─SCH_3, H) > (─CH₃, SH). The differential charge density distributions were examined to underpin prominent electronic contributing factors. Direct fission of C─S bond in the CH₃SH molecule attains exothermic values in the range 2.0 to 2.1 eV. The most energetically favorable sites for the surface-mediated fission of the thiol's S─H bond correspond to the structure where the ─SCH₃ and H are both situated on hollow sites with an adsorption energy of −2.43 eV. Overall, we found that inclusion of van der Waals functional to change the binding energies more noticeably in case of dissociative adsorption structures. The results presented herein should be instrumental in efforts that aim to design stand-alone Fe desulfurization catalysts.     

Bis(tert‐butyl isocyanide‐κC)[4‐fluoro‐N‐({2‐[N‐(4‐fluorophenyl)carboximidoyl]cyclopenta‐2,4‐dien‐1‐ylidene}methyl)anilinido‐κ2N,N′]copper(I)

The solid‐state structure of the title compound, [Cu(C₁₉H₁₃F₂N₂)(C₅H₉N)₂], shows that the Cu^I centre adopts a distorted tetrahedral coordination geometry, being coordinated by two N atoms of the 6‐aminofulvene‐2‐aldimine (AFA) chelating ligand and by the bridgehead C atoms of the two isocyanide ligands. The cyclopentadienyl and imine components of the AFA ligand are approximately coplanar, with an angle between the planes of 5.00 (3)°. The Cu atom lies 0.6460 (3) Å above the imine plane defined by the N and C atoms of the seven‐membered metallocycle. There is also an uncommon C—H...Cu anagostic interaction, with an intramolecular Cu...H distance of 2.67 Å, which is less than the sum of the van der Waals radii.     

Application of computer modeling to surface organometallic chemistry

Computer modeling studies of several surface organometallic complexes have been performed. These studies demonstrate the utility of the modeling technique as well as providing insight into the interaction between the surface organometallic species and the support surface. Models for hydroxylated oxide supports for silica, γ-alumina and magnesia are described. The surface structures for MORh(η³-C₃H₅)₂ (MAl, Si), Mg/[HFeOs₃(CO)₁₃]⁻, and (μ-H) (μ-OMM) Os₃(CO)₁₀ (MAl,Si) are modeled to determine the preferred arrangement of the organometallic species on the hydroxylated support. In the modeling procedures van der Waals interaction energies and non-bonded contacts as a function of orientation of the cluster with respect to the surface are considered. Relaxation of the cluster and Coulombic interaction energies are also considered where appropriate.