Equilibrium structure of a dense trimesic acid monolayer on a homogeneous solid surface: from atomistic simulation to
Sergey S Akimenko1,2, Vitaly A Gorbunov1,2, Eugene A Ustinov1
1Ioffe Institute, 26 Polytechnicheskaya, St. Petersburg, 194021, Russian Federation. vitaly_gorbunov@mail.ru.
Abstract:
A general methodology for determining the thermodynamic characteristics of rigid organic crystals on the atomistic level is presented. The proposed approach is based on a combination of grid interpolation of the precalculated intermolecular potential and kinetic Monte Carlo simulation of the gas-crystal system with an explicit interphase. The two-phase system is stabilized in a wide range of external parameters with an imposed external potential and damping field. The damping field reduces the intermolecular potential at the edges of the crystals and turns it off in the gas phase. To determine the thermodynamic characteristics of a crystal the conditions of equality of chemical potentials in coexisting phases are used. The intermolecular pairwise potential can be calculated on the atomistic or quantum level. In the kinetic Monte Carlo simulations, a grid interpolation of the precalculated potential is performed on each iteration of the algorithm. We have applied the approach to the thermodynamic analysis of a dense monolayer of trimesic acid on a homogeneous surface. The calculated free energy and entropy for the dense "superflower" and filled chicken-wire phases obey the Gibbs-Duhem equation, which confirms the thermodynamic consistency of our approach. Using the proposed approach, we have revealed that the dense "superflower" phase becomes metastable at zero pressure and 470-500 K. Under these conditions, the filled chicken-wire structure with partially released hexagonal cages is thermodynamically favourable. The proposed approach is a potentially universal tool for the thermodynamic analysis of crystals formed by "rigid" organic molecules of any complexity on the atomistic level.
More Related Videos
07:31Author Spotlight: Advancing Cell Membrane Biophysics - Exploring Interactions and Challenges Through Experimental and Computational Approaches
Published on: September 1, 2023
10:52Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics
Published on: April 12, 2019
Related Concept Videos
02:20The Small x Assumption
02:46Crystal Field Theory - Tetrahedral and Square Planar Complexes
Crystal field theory (CFT) is applicable to molecules in geometries other than octahedral. In octahedral complexes, the lobes of the dx2−y2 and dz2 orbitals point directly at the ligands. For tetrahedral complexes, the d orbitals remain in place, but with only four ligands located between the axes. None of the orbitals points directly at the tetrahedral ligands. However, the dx2−y2 and dz2 orbitals (along the Cartesian axes) overlap with the ligands less than the dxy,...
01:09Theories of Dissolution: The Danckwerts' Model and Interfacial Barrier Model
02:42Ionic Crystal Structures
Most monatomic ions behave as charged spheres, and their attraction for ions of opposite charge is the same in every direction. Consequently, stable structures for ionic compounds result (1) when ions of one charge are surrounded by as many ions as possible of the opposite...
03:45Formation of Complex Ions
03:05Acid Strength and Molecular Structure
In the absence of any leveling effect, the acid strength of binary compounds of hydrogen with nonmetals (A) increases as the H-A bond strength decreases down a group in the periodic table. For group 17, the order of increasing acidity is HF < HCl < HBr < HI. Likewise, for group 16, the order of increasing acid strength is H2O < H2S < H2Se < H2Te. Across a row in the periodic table, the acid strength of binary hydrogen compounds increases with...