The Physicochemical Origins of Coincident Epitaxy in Molecular Overlayers: Lattice Modeling vs Potential Energy Calculations

Julie A. Last, Daniel E. Hooks, Andrew C. Hillier, Michael Ward

Research output: Contribution to journalArticle

Abstract

The physicochemical basis for epitaxial stabilization of coincident molecular overlayers is illustrated by comparison of optimum overlayer-substrate configurations calculated with potential energy (PE) methods and a simple geometric lattice misfit modeling algorithm (EpiCalc) that neglects molecular orientations and intermolecular potentials. Using β-bis(ethylenedithio)tetrathiafulvalene triiodide (β-ET2I3), perylenetetracarboxylic diimide (PTCDI), and perylenetetracarboxylic dianhydride (PTCDA) overlayers on a graphite substrate as examples, both methods predict identical optimum azimuthal overlayer orientations for each overlayer that also agree with experimental observations. PE calculations for three hypothetical PTCDA overlayers, with identical lattice parameters but different molecular orientations, predict coincidence at the same azimuthal orientation for all overlayers. Identical results are achieved for PE calculations performed with this lattice when it is occupied by naphthalenetetracarboxylic dianhydride (NTCDA) molecules or argon atoms. These results demonstrate that the epitaxial orientation of coincident overlayers is governed more by geometric lattice matching than specific molecule-substrate interactions and that unambiguous determination of the optimum azimuthal orientation relies on establishing the phase relationship between several overlayer supercells and the substrate. In the case of PE methods, calculations with large overlayer sizes are computationally prohibitive and the energy differences between alternative orientations typically are smaller than the confidence limits of the method. In contrast, the calculation time required by EpiCalc is independent of overlayer size, providing a significant advantage over PE methods with respect to computational speed while enabling unambiguous assignment of the optimum coincident configuration. The reliability of EpiCalc in predicting observed epitaxial overlayer orientations indicates that geometric lattice misfit modeling can be used to screen efficiently for the most favorable epitaxial configuration (overlayer lattice parameters and azimuthal angle), which can then be used in subsequent PE calculations that allow for other degrees of freedom such as molecular orientation.

Original languageEnglish (US)
Pages (from-to)6723-6733
Number of pages11
JournalJournal of Physical Chemistry B
Volume103
Issue number32
StatePublished - Aug 12 1999

Fingerprint

Potential energy
Epitaxial growth
Crystal lattices
epitaxy
potential energy
Molecular orientation
energy methods
Substrates
Lattice constants
Molecules
Graphite
Argon
lattice parameters
configurations
Crystal orientation
confidence limits
Stabilization
Atoms
molecules
graphite

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

The Physicochemical Origins of Coincident Epitaxy in Molecular Overlayers : Lattice Modeling vs Potential Energy Calculations. / Last, Julie A.; Hooks, Daniel E.; Hillier, Andrew C.; Ward, Michael.

In: Journal of Physical Chemistry B, Vol. 103, No. 32, 12.08.1999, p. 6723-6733.

Research output: Contribution to journalArticle

@article{221bd25dd93f473e92ba5a4ea9f06e15,
title = "The Physicochemical Origins of Coincident Epitaxy in Molecular Overlayers: Lattice Modeling vs Potential Energy Calculations",
abstract = "The physicochemical basis for epitaxial stabilization of coincident molecular overlayers is illustrated by comparison of optimum overlayer-substrate configurations calculated with potential energy (PE) methods and a simple geometric lattice misfit modeling algorithm (EpiCalc) that neglects molecular orientations and intermolecular potentials. Using β-bis(ethylenedithio)tetrathiafulvalene triiodide (β-ET2I3), perylenetetracarboxylic diimide (PTCDI), and perylenetetracarboxylic dianhydride (PTCDA) overlayers on a graphite substrate as examples, both methods predict identical optimum azimuthal overlayer orientations for each overlayer that also agree with experimental observations. PE calculations for three hypothetical PTCDA overlayers, with identical lattice parameters but different molecular orientations, predict coincidence at the same azimuthal orientation for all overlayers. Identical results are achieved for PE calculations performed with this lattice when it is occupied by naphthalenetetracarboxylic dianhydride (NTCDA) molecules or argon atoms. These results demonstrate that the epitaxial orientation of coincident overlayers is governed more by geometric lattice matching than specific molecule-substrate interactions and that unambiguous determination of the optimum azimuthal orientation relies on establishing the phase relationship between several overlayer supercells and the substrate. In the case of PE methods, calculations with large overlayer sizes are computationally prohibitive and the energy differences between alternative orientations typically are smaller than the confidence limits of the method. In contrast, the calculation time required by EpiCalc is independent of overlayer size, providing a significant advantage over PE methods with respect to computational speed while enabling unambiguous assignment of the optimum coincident configuration. The reliability of EpiCalc in predicting observed epitaxial overlayer orientations indicates that geometric lattice misfit modeling can be used to screen efficiently for the most favorable epitaxial configuration (overlayer lattice parameters and azimuthal angle), which can then be used in subsequent PE calculations that allow for other degrees of freedom such as molecular orientation.",
author = "Last, {Julie A.} and Hooks, {Daniel E.} and Hillier, {Andrew C.} and Michael Ward",
year = "1999",
month = "8",
day = "12",
language = "English (US)",
volume = "103",
pages = "6723--6733",
journal = "Journal of Physical Chemistry B Materials",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "32",

}

TY - JOUR

T1 - The Physicochemical Origins of Coincident Epitaxy in Molecular Overlayers

T2 - Lattice Modeling vs Potential Energy Calculations

AU - Last, Julie A.

AU - Hooks, Daniel E.

AU - Hillier, Andrew C.

AU - Ward, Michael

PY - 1999/8/12

Y1 - 1999/8/12

N2 - The physicochemical basis for epitaxial stabilization of coincident molecular overlayers is illustrated by comparison of optimum overlayer-substrate configurations calculated with potential energy (PE) methods and a simple geometric lattice misfit modeling algorithm (EpiCalc) that neglects molecular orientations and intermolecular potentials. Using β-bis(ethylenedithio)tetrathiafulvalene triiodide (β-ET2I3), perylenetetracarboxylic diimide (PTCDI), and perylenetetracarboxylic dianhydride (PTCDA) overlayers on a graphite substrate as examples, both methods predict identical optimum azimuthal overlayer orientations for each overlayer that also agree with experimental observations. PE calculations for three hypothetical PTCDA overlayers, with identical lattice parameters but different molecular orientations, predict coincidence at the same azimuthal orientation for all overlayers. Identical results are achieved for PE calculations performed with this lattice when it is occupied by naphthalenetetracarboxylic dianhydride (NTCDA) molecules or argon atoms. These results demonstrate that the epitaxial orientation of coincident overlayers is governed more by geometric lattice matching than specific molecule-substrate interactions and that unambiguous determination of the optimum azimuthal orientation relies on establishing the phase relationship between several overlayer supercells and the substrate. In the case of PE methods, calculations with large overlayer sizes are computationally prohibitive and the energy differences between alternative orientations typically are smaller than the confidence limits of the method. In contrast, the calculation time required by EpiCalc is independent of overlayer size, providing a significant advantage over PE methods with respect to computational speed while enabling unambiguous assignment of the optimum coincident configuration. The reliability of EpiCalc in predicting observed epitaxial overlayer orientations indicates that geometric lattice misfit modeling can be used to screen efficiently for the most favorable epitaxial configuration (overlayer lattice parameters and azimuthal angle), which can then be used in subsequent PE calculations that allow for other degrees of freedom such as molecular orientation.

AB - The physicochemical basis for epitaxial stabilization of coincident molecular overlayers is illustrated by comparison of optimum overlayer-substrate configurations calculated with potential energy (PE) methods and a simple geometric lattice misfit modeling algorithm (EpiCalc) that neglects molecular orientations and intermolecular potentials. Using β-bis(ethylenedithio)tetrathiafulvalene triiodide (β-ET2I3), perylenetetracarboxylic diimide (PTCDI), and perylenetetracarboxylic dianhydride (PTCDA) overlayers on a graphite substrate as examples, both methods predict identical optimum azimuthal overlayer orientations for each overlayer that also agree with experimental observations. PE calculations for three hypothetical PTCDA overlayers, with identical lattice parameters but different molecular orientations, predict coincidence at the same azimuthal orientation for all overlayers. Identical results are achieved for PE calculations performed with this lattice when it is occupied by naphthalenetetracarboxylic dianhydride (NTCDA) molecules or argon atoms. These results demonstrate that the epitaxial orientation of coincident overlayers is governed more by geometric lattice matching than specific molecule-substrate interactions and that unambiguous determination of the optimum azimuthal orientation relies on establishing the phase relationship between several overlayer supercells and the substrate. In the case of PE methods, calculations with large overlayer sizes are computationally prohibitive and the energy differences between alternative orientations typically are smaller than the confidence limits of the method. In contrast, the calculation time required by EpiCalc is independent of overlayer size, providing a significant advantage over PE methods with respect to computational speed while enabling unambiguous assignment of the optimum coincident configuration. The reliability of EpiCalc in predicting observed epitaxial overlayer orientations indicates that geometric lattice misfit modeling can be used to screen efficiently for the most favorable epitaxial configuration (overlayer lattice parameters and azimuthal angle), which can then be used in subsequent PE calculations that allow for other degrees of freedom such as molecular orientation.

UR - http://www.scopus.com/inward/record.url?scp=0000260857&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0000260857&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0000260857

VL - 103

SP - 6723

EP - 6733

JO - Journal of Physical Chemistry B Materials

JF - Journal of Physical Chemistry B Materials

SN - 1520-6106

IS - 32

ER -