Understanding the interactions between molecules on surfaces is crucial for advancing technologies in sensing, catalysis, and energy harvesting. In this study we explore the complex surface chemistry resulting from the interaction of Co(II)octaethylporphyrin (CoOEP) and iodine, I2, both in solution and at the phenyloctane/HOPG interface. In pursuit of this goal, we report results from electrochemistry, NMR and UV-Vis spectroscopy, X-ray crystallography, scanning tunneling microscopy (STM), and density functional theory (DFT). Both spectroscopic methods of analysis confirmed that at and above the stoichiometric ratio of one CoOEP to one I2 the reaction product was metal centered CoIII(OEP)I. X-ray crystallography verified that a single iodine is bonded to each cobalt ion in the triclinic, P-1 system. The surface chemistry of CoOEP and I2 is complicated and remarkably dependent on the iodine concentration. STM images of CoOEP and I2 in phenyloctane on highly oriented pyrolytic graphite (HOPG) at low halogen concentrations (1:<2 Co:I ratios) presented random individual Co(OEP)I molecules weakly adsorbed onto a hexagonal (HEX) CoOEP monolayer. Images of 1:2 Co:I ratio solutions, showed phase segregated HEX CoOEP and pseudo-rectangular (REC) Co(OEP)I incorporating one solvent molecule per Co(OEP)I. The REC structure formed in long parallel rows with the number of rows increasing with increasing solution I2. In this case, the presence of CoOEP on the surface was attributed to the spontaneous reduction of Co(OEP)I by the graphite substrate. DFT calculations indicate that the REC Co(OEP)I:PhO form is energetically more stable than the HEX form of Co(OEP)I on HOPG. Experimental STM images and DFT calculated adsorption energies and STM images support our interpretation of the observed structures.
more »
« less
Computational Study of Solvent Incorporation into a Porphyrin Monolayer
Density functional theory (DFT) is used to investigate the conversion from a solvent incorporated pseudo-polymorph into a single component monolayer. Calculations of thermodynamic properties both for the surfaces in contact with gas phase and with solvent are reported. In the case of wetted surfaces, a simple bond-additivity model, first proposed by Campbell and modified here, is used to augment the DFT calculations. The model predicts a dramatic reduction in desorption energies in solvent as compared to gas phase. Eyring’s reaction rate theory is used to predict limiting desorption rates for guest (solvent) molecules from the pockets in the pseudo-polymorph and for cobalt octaethylporphyrin (COEP) molecules in all structures. The pseudo-polymorph studied here is a nearly rectangular lattice (REC) composed of two CoOEP and 2 molecules of either 1,2,4-trichlorobenzene (TCB) or toluene (TOL) supported on 63 atoms of Au(111). At sufficiently high initial concentrations of CoOEP, only a hexagonal unit cell (HEX) with two molecules of CoOEP, supported on 50 atoms of gold is observed. Experimentally, the TCB-REC structure is more stable than the TOL-REC structure existing in solution at initial mM concentrations of CoOEP in TCB as opposed to initial M concentration of CoOEP in toluene. Calculations here show that the HEX structure is the thermodynamically stable structure at all practical concentrations of CoOEP. Once the REC structure forms kinetically at low concentration because of the vast excess of solvent on the surface, it is difficult to convert to the more stable HEX structure. The difference in stability is primarily due to the difference in electronic adsorption energy of the solvents (TOL or TCB) and to the very low desorption rate of CoOEP. The adsorption energy of TCB has two important contributors: the adsorption energy onto Au alone, and the intermolecular interactions between TCB and the CoOEP host lattice. Neither factor can be neglected. We also find that planar adsorption of both TOL and TCB on Au(111) is the energetically preferred orientation when space is available on the surface. Rates of desorption are very sensitive to the solvent free activation energy and to the thermodynamic parameters required to convert the solvent free activation energy to one for the solvated surface. Small changes in the computed energy (of the order of 5%) can lead to one order of magnitude change in rates. Further, the solvation model used does not provide the barrier to adsorption in solution needed to determine values for the desorption activation energy. Thus, the rates computed here for desorption into solvent are limiting values.
more »
« less
- Award ID(s):
- 2306316
- PAR ID:
- 10506700
- Publisher / Repository:
- The American Chemical Society
- Date Published:
- Journal Name:
- The Journal of Physical Chemistry C
- Volume:
- 128
- Issue:
- 4
- ISSN:
- 1932-7447
- Page Range / eLocation ID:
- 1827 to 1839
- Subject(s) / Keyword(s):
- adsorption polymorph STM solvent effects liquid-solid interface equilibrium kinetic control
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
One of the common practices in the literature of molecular desorption is the comparison of theoretically (mostly using DFT) calculated single molecule adsorption energies with experimental desorption energies from studies like temperature programmed desorption (TPD) etc. Comparisons like those do not consider that the experimental desorption energies are obtained via ensemble techniques while theoretical values are calculated at the single molecule level. Theoretical values are generally based upon desorption of a single molecule from a clean surface, or upon desorption of an entire monolayer. On the other hand, coverage dependent molecule–molecule interactions add to and modify molecule–substrate interactions that contribute to the experimentally determined desorption energies. In this work, we explore the suitability of an additive nearest neighbor model for determining general coverage dependent single molecule desorption energies in non-covalent self-assembled monolayers (SAMs). These coverage dependent values serve as essential input to any model attempting to reproduce coverage dependent desorption or for understanding the time dependent desorption from a partially covered surface. This method is tested using a case study of coronene adsorbed on Au(111) and HOPG substrates with periodic DFT calculations. Calculations show that coronene exhibits coverage and substrate dependence in molecular desorption. We found that intermolecular contact energies in the coronene monolayer are not strongly influenced by the HOPG substrate, while coronene desorption on Au(111) exhibits strong cooperativity where the additive model fails.more » « less
-
Adsorption involves molecules colliding at the surface of a solid and losing their incidence energy by traversing a dynamical pathway to equilibrium. The interactions responsible for energy loss generally include both chemical bond formation (chemisorption) and nonbonding interactions (physisorption). In this work, we present experiments that revealed a quantitative energy landscape and the microscopic pathways underlying a molecule’s equilibration with a surface in a prototypical system: CO adsorption on Au(111). Although the minimum energy state was physisorbed, initial capture of the gas-phase molecule, dosed with an energetic molecular beam, was into a metastable chemisorption state. Subsequent thermal decay of the chemisorbed state led molecules to the physisorption minimum. We found, through detailed balance, that thermal adsorption into both binding states was important at all temperatures.more » « less
-
Abstract Fullertubes, that is, fullerenes consisting of a carbon nanotube moiety capped by hemifullerene ends, are emerging carbon nanomaterials whose properties show both fullerene and carbon nanotube (CNT) traits. Albeit it may be expected that their electronic states show a certain resemblance to those of the extended nanotube, such a correlation has not yet been found or described. Here it shows a scanning tunneling microscopy (STM) and spectroscopy (STS) characterization of the adsorption, self‐assembly, and electronic structure of 2D arrays of [5,5]‐C90fullertube molecules on two different noble metal surfaces, Ag(111) and Au(111). The results demonstrate that the shape of the molecular orbitals of the adsorbed fullertubes corresponds closely to those expected for isolated species on the grounds of density functional theory calculations. Moreover, a comparison between the electronic density profiles in the bands of the extended [5,5]‐CNT and in the molecules reveals that some of the frontier orbitals of the fullertube molecules can be described as the result of the quantum confinement imposed by the hemifullerene caps to the delocalized band states in the extended CNT. The results thus provide a conceptual framework for the rational design of custom fullertube molecules and can potentially become a cornerstone in the understanding of these new carbon nanoforms.more » « less
-
H2 activation is fundamental in catalysis. Single-atom catalysts (SACs) can be highly selective hydrogenation catalysts due to their tunable geometric and electronic properties. In this work, H2 activation (adsorption, splitting, and diffusion) on the anatase TiO2-supported SAC has been modeled in detail. The stable configurations of 14 transition metals from 3d to 5d (Fe, Co, Ni, Cu, Zn, Ru, Rh, Pd, Ag, Cd, Os, Ir, Pt, and Au) and Sn have been screened. We compared H and H2 adsorption and H2 heterolytic and homolytic splitting on SA/TiO2. H on the SAC in neutral, hydridic, and proton forms and the preferred H2 dissociation paths are revealed. We found that the metal adatoms strengthen the Brønsted acids via forming the SA-O bonds and promote the H adsorption on Ti sites via forming the Ti3+ sites. The electronic descriptor using the energy level of the frontier d orbital, referenced to vacuum, can predict the single H and H2 dissociative adsorption energies on the metal site. As the SA-Hδ- interaction is stronger than Ti-Hδ-, the activation barriers for heterolytic paths over SA-O sites are lower than over Ti-O sites. H2 adsorption is activated on Au, Ru, Rh, Pd, and Ir in a dihydrogen complex structure with an elongated H-H bond. Homolytic splitting over SA sites is favored thermodynamically and kinetically on Rh, Pd, Os, Ir, and Pt. In contrast, for the remaining SA/TiO2, H-H splitting at the SA-O is kinetically favored compared to the Ti-O sites, but the products are less thermodynamically favored.more » « less
An official website of the United States government

