Matthew Connolly

Explicit Hydrogen Molecular Dynamics Simulations of Hexane Deposited onto Graphite at Various Coverages

We present results of molecular dynamics (MD) computer simulations of hexane (C6H14) adlayers physisorbed onto a graphite substrate for coverages in the range 0.5 < or = rho < or = 1 monolayers. The hexane molecules are simulated with explicit hydrogens, and the graphite substrate is modeled as an all-atom structure having six graphene layers. At coverages above about rho congruent with 0.9 the low-temperature herringbone solid loses its orientational order at T(1) = 140 +/- 3 K. At rho = 0.878, the system presents vacancy patches and T(1) decreases to ca. 100 K. As coverage decreases further, the vacancy patches become larger and by rho = 0.614 the solid is a connected network of randomly oriented islands and there is no global herringbone order-disorder transition. In all cases we observe a weak nematic mespohase. The melting temperature for our explicit-hydrogen model is T(2) = 160 +/- 3 K and falls to ca. 145 K by rho = 0.614 (somewhat lower than seen in experiment). The dynamics seen in the fully atomistic model agree well with experiment, as the molecules remain overall flat on the substrate in the solid phase and do not show anomalous tilting behavior at any phase transition observed in earlier simulations in the unified atom (UA) approximation. Energetics and structural parameters also are more reasonable and, collectively, the results from the simulations in this work demonstrate that the explicit-hydrogen model of hexane is substantially more realistic than the UA approximation.

Structural and phase properties of tetracosane (C24H50) monolayers adsorbed on graphite: an explicit hydrogen molecular dynamics study.

We discuss molecular dynamics (MD) computer simulations of a tetracosane (C24H50) monolayer physisorbed onto the basal plane of graphite. The adlayer molecules are simulated with explicit hydrogens, and the graphite substrate is represented as an all-atom structure having six graphene layers. The tetracosane dynamics modeled in the fully atomistic manner agree well with experiment. The low-temperature ordered solid organizes into a rectangularly centered structure that is not commensurate with underlying graphite. Above T=200 K, as the molecules start to lose their translational and orientational order via gauche defect formation a weak smectic mesophase (observed experimentally but never reproduced in united atom (UA) simulations) appears. The phase behavior of the adsorbed layer is critically sensitive to the way the electrostatic interactions are included in the model. If the electrostatic charges are set to zero (as for a UA force field), then the melting temperature increases by approximately 70 K with respect to the experimental value. When the nonbonded 1-4 interaction is not scaled, the melting temperature decreases by approximately 90 K. If the scaling factor is set to 0.5, then melting occurs at T=350 K, in very good agreement with experimental data.