### 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 (β-ET_{2}I_{3}), 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 language | English (US) |
---|---|

Pages (from-to) | 6723-6733 |

Number of pages | 11 |

Journal | Journal of Physical Chemistry B |

Volume | 103 |

Issue number | 32 |

State | Published - Aug 12 1999 |

### Fingerprint

### ASJC Scopus subject areas

- Physical and Theoretical Chemistry

### Cite this

*Journal of Physical Chemistry B*,

*103*(32), 6723-6733.

**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.

Research output: Contribution to journal › Article

*Journal of Physical Chemistry B*, vol. 103, no. 32, pp. 6723-6733.

}

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 -