Adsorption studies of C6H6 on Cu (111), Ag (111), and Au (111) within dispersion corrected density functional theory

Probing the thermal resistance of solid-liquid interfaces in nanofluids with molecular dynamics

The significance of interfacial thermal resistance in the thermal conductivity of nanofluids is not well understood, in part because of the absence of measurements of this quantity. Here, we study the interfacial thermal resistance for metal-oil nanofluids of interest as heat transfer fluids for concentrating solar power, using density functional theory and molecular dynamics simulations. Insights on the role of chemical interactions in determining the interfacial thermal resistance are revealed. The results presented here showcase a general picture in which the stronger the chemical interactions between species at the interface, the lower the associated interfacial thermal resistance. The implications toward nanofluid design are discussed. We show that, for this important family of metal-oil nanofluids, the interfacial thermal resistance values are low enough so that it is possible to afford a reduction in particle size, minimizing stability and rheological issues while still offering enhancement in the effective thermal conductivity with respect to the base fluid.

Dispersion-Corrected DFT-D4 Study of the Adsorption of Halobenzenes and 1,3,5-Trihalobenzenes on the Cu(111) Surface─Effect of Sigma Hole Properties

Halogen bonds, characterized by directionality, tunability, hydrophobicity, and variable sizes, are ideal noncovalent interactions to design and control the formation of self-assembled nanostructures. The specific self-assembly cases formed by the halogen-bonding interaction have been well studied by scanning tunneling microscopy (STM) experiments and density functional theory (DFT) calculations. However, there is a lack of systematic theoretical adsorption studies on halogenated molecules. In this work, the adsorption of halobenzenes and 1,3,5-trihalobenzenes on the Cu(111) surface was examined by dispersion-corrected DFT methods. The adsorption geometries, noncovalent molecule-surface interactions, electronic densities, and electrostatic potential maps were examined for their most stable adsorption sites using the DFT-D4 method. Our calculations revealed that the iodo compounds favor a different adsorption geometry from aryl chlorides and bromides. Down the halogen group (Cl to I), the adsorption energy increases and the distance between the halogen atom and Cu surface decreases, which indicates stronger molecule-surface interactions. This is supported by the changes in the density of states upon adsorption. Noncovalent interaction analysis was also employed to further understand the nature and relative strength of the molecule-surface interactions. Electrostatic potential maps revealed that the positive character of the halogen sigma hole becomes stronger upon adsorption. Thus, surface adsorption of the halogenated molecule will enhance the formation of intermolecular halogen bonds. The present theoretical findings are expected to contribute toward a more comprehensive understanding of halogen bonding on the Cu(111) surface.

A first principle evaluation of the adsorption mechanism and stability of volatile organic compounds into NaY zeolite

Removal of volatile organic compounds (VOCs) from indoor or outdoor environments is an urgent challenge for the protection of human populations. Inorganic sorbents such as zeolites are a promising solution to tackle this issue. Using dispersion corrected periodic DFT calculations, we have studied the interaction between sodium-exchanged faujasite zeolite and a large set of VOCs including aromatics, oxygenates and chlorinated compounds. The computed interaction energies range from about −25 (methane) to −130 kJ/mol (styrene). Methane is by far the less interacting specie with the NaY zeolite. All other VOCs present interaction energies higher in absolute value than 69 kJ/mol. Most of them show a similar adsorption strength, between −70 and −100 kJ/mol. While the electrostatic interactions are important in the case of oxygenates and acrylonitrile, van der Waals interactions predominate in hydrocarbons and chlorides. By monitoring the variation of molecular bond lengths of the different VOCs before and after adsorption, we have then evaluated the tendency of adsorbate to react and form by-products, since a significant stretching would evidently lead to the activation of the bond. While hydrocarbons, tetrachloroethylene and acrylonitrile seem to be not activated upon adsorption, all oxygenates and 1,1,2-trichloroethane could possibly react once adsorbed.

First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications

Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.

Surface Adsorption from the Exchange-Hole Dipole Moment Dispersion Model

The accurate calculation of intermolecular interaction energies with density functional theory requires methods that include a treatment of long-range, nonlocal dispersion correlation. In this work, we explore the ability of the exchange-hole dipole moment (XDM) dispersion correction to model molecular surface adsorption. Adsorption energies are calculated for six small aromatic molecules (benzene, furan, pyridine, thiophene, thiophenol, and benzenediamine) and the four DNA nucleobases (adenine, thymine, guanine, and cytosine) on the (111) surfaces of the three coinage metals (copper, silver, and gold). For benzene, where the experimental reference data is most precise, the mean absolute error in the computed absorption energies is 0.04 eV. For the other aromatic molecules, the computed binding energies are found to be within 0.09 eV of the available reference data, on average, which is well below the expected experimental uncertainties for temperature-programmed desorption measurements. Unlike other dispersion-corrected functionals, adequate performance does not require changes to the canonical XDM implementation, and the good performance of XDM is explained in terms of the behavior of the exchange hole. Additionally, the base functional employed (B86bPBE) is also optimal for molecular studies, making B86bPBE-XDM an excellent candidate for studying chemistry on material surfaces. Finally, the noncovalent interaction (NCI) plot technique is shown to detect adsorption effects in real space on the order of tenths of an eV.

Understanding and Designing the Gold-Bio Interface: Insights from Simulations

Gold nanoparticles (AuNPs) are an integral part of many exciting and novel biomedical applications, sparking the urgent need for a thorough understanding of the physicochemical interactions occurring between these inorganic materials, their functional layers, and the biological species they interact with. Computational approaches are instrumental in providing the necessary molecular insight into the structural and dynamic behavior of the Au-bio interface with spatial and temporal resolutions not yet achievable in the laboratory, and are able to facilitate a rational approach to AuNP design for specific applications. A perspective of the current successes and challenges associated with the multiscale computational treatment of Au-bio interfacial systems, from electronic structure calculations to force field methods, is provided to illustrate the links between different approaches and their relationship to experiment and applications.

Understanding molecular self-assembly of a diol compound by considering competitive interactions

With a combination of scanning tunneling microscopy and density functional theory, effects on molecular self-assembly involving two distinct chemical groups were investigated.

Influence of steric hindrance on the molecular packing and the anchoring of quinonoid zwitterions on gold surfaces

The N-substituent on quinonoid zwitterions influences the molecules packing and impacts their anchoring on gold surfaces.

Reversible potential-induced switching of alkyl chain aggregation in octyl-triazatriangulenium adlayers on Au(111)

In situ scanning tunneling microscopy and cyclic voltammetry studies of self-assembled octyl-triazatriangulenium monolayers on Au(111) electrode surfaces in 0.1 M HClO4 reveal a complex surface phase behavior, involving two fast, highly reversible transitions between different ordered adlayer phases: With decreasing potential, the preadsorbed (√19 × √19)R23.4° adlayer first is converted into a (7√3 × 7√3) and then into a (2√3 × 2√3)R30° phase, corresponding to a stepwise increase in the local packing density of the molecules. The (7√3 × 7√3) → (2√3 × 2√3)R30° transition is accompanied by a reorientation of the peripheral octyl chains from a more planar to a close-packed vertical arrangement. This reversible potential-induced switching between a homogeneous adlayer of small vertical extension and a Au surface partially covered by islands of a compact hydrocarbon layer is attributed to changes in the adsorbate charge state and associated changes in the intermolecular interactions.

Biomolecular adsorption at aqueous silver interfaces: first-principles calculations, polarizable force-field simulations, and comparisons with gold

The molecular simulation of biomolecules adsorbed at noble metal interfaces can assist in the development of bionanotechnology applications. In line with advances in polarizable force fields for adsorption at aqueous gold interfaces, there is scope for developing a similar force field for silver. One way to accomplish this is via the generation of in vacuo adsorption energies calculated using first-principles approaches for a wide range of different but biologically relevant small molecules, including water. Here, we present such first-principles data for a comprehensive range of bio-organic molecules obtained from plane-wave density functional theory calculations using the vdW-DF functional. As reported previously for the gold force field, GolP-CHARMM (Wright, L. B.; Rodger, P. M.; Corni, S.; Walsh, T. R. GolP-CHARMM: first-principles based force-fields for the interaction of proteins with Au(111) and Au(100). J. Chem. Theory Comput. 2013, 9, 1616-1630), we have used these data to construct a a new force field, AgP-CHARMM, suitable for the simulation of biomolecules at the aqueous Ag(111) and Ag(100) interfaces. This force field is derived to be consistent with GolP-CHARMM such that adsorption on Ag and Au can be compared on an equal footing. Our force fields are used to evaluate the water overlayer stability on both silver and gold, finding good agreement with known behaviors. We also calculate and compare the structuring (spatial and orientational) of liquid water adsorbed at both silver and gold. Finally, we report the adsorption free energy of a range of amino acids at both the Au(111) and Ag(111) aqueous interfaces, calculated using metadynamics. Stronger adsorption on gold was noted in most cases, with the exception being the carboxylate group present in aspartic acid. Our findings also indicate differences in the binding free energy profile between silver and gold for some amino acids, notably for His and Arg. Our analysis suggests that the relatively stronger structuring of the first water layer on silver, relative to gold, could give rise to these differences.

Adsorption of diferrocenylacetylene on Au(111) studied by scanning tunneling microscopy

Scanning tunneling microscopy images of diferrocenylacetylene (DFA) coadsorbed with benzene on Au(111) show individual and close-packed DFA molecules, either adsorbed alongside benzene or on top of a benzene monolayer. Images acquired over a range of positive and negative tip-sample bias voltages show a shift in contrast, with the acetylene linker appearing brighter than the ferrocenes at positive sample bias (where unoccupied states primarily contribute) and the reverse contrast at negative bias. Density functional theory was used to calculate the electronic structure of the gas-phase DFA molecule, and simulated images produced through two-dimensional projections of these calculations approximate the experimental images. The symmetry of both experimental and calculated molecular features for DFA rules out a cis adsorption geometry, and comparison of experiment to simulation indicates torsion around the inter-ferrocene axis between 90° and 180° (trans); the cyclopentadienyl rings are thus angled with respect to the surface.