The Self-Assembly of Nano-Objects Code: Applications to supramolecular organic monolayers adsorbed on metal surfaces

The Self-Assembly of Nano-Objects (SANO) code we implemented demonstrates the ability to predict the molecular self-assembly of different structural motifs by tuning the molecular building blocks as well as the metallic substrate. It consists in a two-dimensional Grand Canonical Monte-Carlo (GCMC) approach developed to perform atomistic simulations of thousands of large organic molecules self-assembling on metal surfaces. Computing adsorption isotherms at room temperature and spanning over the characteristic sub-micrometric scales, we confront the robustness of the approach with three different well-known systems: ZnPcCl8 on Ag(111), CuPcF16 on Au(111) and PTBC on Ag(111). We retrieve respectively their square, oblique and hexagonal supramolecular tilling. The code incorporates generalized force fields to describe the molecular interactions, which provides transferability and versatility to many organic building blocks and metal surfaces.

💡 Research Summary

The paper introduces SANO (Self‑Assembly of Nano‑Objects), a computational framework designed to predict the self‑assembly of large organic molecules on metal surfaces with atomistic detail. The core of the methodology is a two‑dimensional Grand Canonical Monte‑Carlo (GCMC) algorithm that directly controls the chemical potential, allowing the simulation to reproduce adsorption isotherms at room temperature and to explore a wide range of surface coverages. Molecules are treated as rigid or semi‑rigid bodies with translational and rotational degrees of freedom, while the metal substrate is represented by a fixed (111) lattice of interaction sites. Inter‑molecular and molecule‑substrate forces are described by a generalized force field that combines Lennard‑Jones dispersion, Coulombic electrostatics, and specific terms for π‑π stacking and hydrogen bonding. The force‑field parameters are calibrated against experimental data and density‑functional theory calculations, providing a transferable description that can be applied to many different organic building blocks and metal surfaces without re‑parameterization.

Implementation details focus on scalability: the code employs spatial decomposition and neighbor‑list updates to enable parallel execution, allowing simulations of several thousand molecules—far beyond the capability of traditional Monte‑Carlo approaches that are limited to a few hundred. Each Monte‑Carlo cycle consists of insertion, deletion, translation, and rotation attempts, accepted according to the Metropolis criterion. By running tens to hundreds of millions of cycles, the system reaches equilibrium and yields statistically reliable structural information.

The authors validate SANO on three benchmark systems that are well documented in the literature: (i) octachloro‑zinc phthalocyanine (ZnPcCl8) on Ag(111), which experimentally forms a square lattice; (ii) hexadecafluoro‑copper phthalocyanine (CuPcF16) on Au(111), known to produce an oblique (tilted square) lattice; and (iii) PTBC (1,3,5‑triphenyl‑benzene‑carboxylate) on Ag(111), which assembles into a hexagonal tiling. For each case the simulated adsorption isotherms match experimental measurements, and the predicted supramolecular tilings reproduce the correct symmetry, lattice constants, and domain orientations. Temperature variations in the simulations capture the experimentally observed phase transitions, such as domain growth, defect formation, and lattice distortion.

The analysis highlights the physical mechanisms governing each system. In ZnPcCl8/Ag(111) the dominant π‑π interactions between planar macrocycles drive square ordering, while the weak electronic coupling to the silver surface stabilizes the lattice without inducing significant charge transfer. In CuPcF16/Au(111) the fluorine substituents modify the electrostatic landscape, leading to a slight shear of the square lattice into an oblique geometry; the gold surface contributes a modest charge‑induced dipole that further tilts the arrangement. PTBC/Ag(111) demonstrates that molecular shape alone can dictate a hexagonal network, with the silver substrate acting mainly as a template that minimizes strain and defect formation.

The authors discuss both strengths and limitations of the approach. Strengths include (1) the ability to handle thousands of large molecules, (2) direct control of experimental variables (temperature, chemical potential), and (3) a transferable force‑field framework that reduces the need for system‑specific re‑parameterization. Limitations arise from the rigid‑body approximation (neglecting intramolecular vibrations), the simplified representation of the metal surface as a static lattice (which omits explicit electronic structure effects), and the reliance on classical force fields that may not capture subtle many‑body dispersion at very high coverages. The paper proposes future extensions such as integrating density‑functional‑theory‑derived potentials, employing machine‑learning corrections to the force field, and developing hybrid quantum‑mechanical/molecular‑mechanical (QM/MM) schemes to treat the metal‑adsorbate interface more accurately.

In conclusion, SANO provides a robust, versatile, and computationally efficient tool for predicting and rationalizing the self‑assembly of organic monolayers on metal substrates. Its capacity to reproduce experimentally observed square, oblique, and hexagonal supramolecular tilings demonstrates its potential to guide the design of functional nanostructured interfaces in organic electronics, catalysis, and sensor technologies.