Predicting Organometallic Intermediates in the Surface-Assisted Ullmann Coupling of Chrysene Isomers.

On-surface polymerization of functional organic molecules has been recently recognized as a promising route to persistent low-dimensional structures with tailorable properties. In this contribution, using the coarse-grained Monte Carlo simulation method, we study the initial stage of the Ullmann coupling of doubly halogenated chrysene isomers adsorbed on a catalytically active (111) crystalline surface. To that end, we focus on the formation of labile metal-organic precursor structures preceding the covalent bonding of chrysene monomers. Four monomeric chrysene units with differently distributed halogen substituents were probed in the simulations, and the resulting precursor structures were compared and quantified. Moreover, the effect of (pro)chirality of chrysene tectons on the structure formation was elucidated by running separate simulations in enantiopure and racemic systems. The calculations showed that suitable manipulation of the halogen substitution pattern allows for the creation of diverse precursor architectures, ranging from straight and winded chains to cyclic oligomers with enantiopure, racemic, and nonracemic composition. The obtained findings can be helpful in developing synthetic strategies for covalent polymers with predefined architecture and functionality.

Structural Quantification of the Surface-Confined Metal-Organic Precursors Simulated with the Lattice Monte Carlo Method.

The diversity of surface-confined metal-organic precursor structures, which recently have been observed experimentally, poses a question of how the individual properties of a molecular building block determine those of the resulting superstructure. To answer this question, we use the Monte Carlo simulation technique to model the self-assembly of metal-organic precursors that precede the covalent polymerization of halogenated PAH isomers. For this purpose, a few representative examples of low-dimensional constructs were studied, and their basic structural features were quantified using such descriptors as the orientational order parameter, radial distribution function, and one- and two-dimensional structure factors. The obtained results demonstrated that the morphology of the precursor (and thus the subsequent polymer) could be effectively tuned by a suitable choice of molecular parameters, including size, shape, and intramolecular distribution of halogen substituents. Moreover, our theoretical investigations showed the effect of the main structural features of the precursors on the related indirect characteristics of these constructs. The results reported herein can be helpful in the custom designing and characterization of low-dimensional polymers with adjustable properties.

Theoretical Modeling of the Metal-Organic Precursors of Anthracene-Based Covalent Networks on Surfaces.

Surface-assisted fabrication of molecular network architectures has been a promising route to low-dimensional materials with unique physicochemical properties and functionalities. One versatile way in this field is the Ullmann coupling reaction of halogenated organic monomers on catalytically active metallic surfaces. In this work, using the coarse-grained Monte Carlo simulations, we studied the on-surface self-assembly of metal-organic precursors preceding the covalent Ullman-type linkage of tetrahalogenated anthracene building blocks. To that end, a series of positional isomers was examined and classified with respect to their ability of creation of extended network structures. Our simulations focused on the identification of basic types of self-assembly scenarios distinguishing enantiopure and racemic systems and producing periodic and aperiodic networks. The calculations carried out for selected tectons demonstrated wide possibilities of controlling porosity (e. g. pore size, shape, periodicity, chirality, heterogeneity) of the networks by suitable functionalization of the monomeric unit. The findings reported herein can be helpful in rational designing of 2D polymeric networks with predefined structures and properties.

Designing 2D covalent networks with lattice Monte Carlo simulations: precursor self-assembly.

Organic synthesis reactions in the adsorbed phase have been recently an intensively studied topic in heterogeneous catalysis and material engineering. One of such processes is the Ullmann coupling in which halogenated organic monomers are transformed into covalently bonded polymeric structures. In this work, we use the lattice Monte Carlo simulation method to study the on-surface self-assembly of organometallic precursor architectures comprising tetrasubstituted naphthalene building blocks with differently distributed halogen atoms. In the coarse grained approach adopted herein the molecules and metal atoms were modeled by discrete segments, two connected and one, respectively, placed on a triangular lattice representing a (111) metallic surface. Our simulations focused on the influence of the intramolecular distribution of the substituents on the morphology of the resulting superstructures. Special attention was paid to the molecules that create porous networks characterized by long-range order. Moreover, the structural analysis of the assemblies comprising prochiral building blocks was made by running simulations for the corresponding enantiopure and racemic adsorbed systems. The obtained results demonstrated the possibility of directing the on-surface self-assembly towards networks with controllable pore shape and size. These findings can be helpful in designing covalently bonded 2D superstructures with predefined architecture and functions.

The Macrocycle versus Chain Competition in On-Surface Polymerization: Insights from Reactions of 1,3-Dibromoazulene on Cu(111).

Ring/chain competition in oligomerization reactions represents a long-standing topic of synthetic chemistry and was treated extensively for solution reactions but is not well-understood for the two-dimensional confinement of surface reactions. Here, the kinetic and thermodynamic principles of ring/chain competition in on-surface synthesis are addressed by scanning tunneling microscopy, X-ray photoelectron spectroscopy, and Monte Carlo simulations applied to azulene-based organometallic oligomers on Cu(111). Analysis of experiments and simulations reveals how the ring/chain ratio can be controlled through variation of coverage and temperature. At room temperature, non-equilibrium conditions prevail and kinetic control leads to preferential formation of the entropically favored chains. In contrast, high-temperature equilibrium conditions are associated with thermodynamic control, resulting in increased yields of the energetically favored rings. The optimum conditions for ring formation include the lowest possible temperature within the regime of thermodynamic control and a low coverage. The general implications are discussed and compared to the solution case.

Theoretical Modeling of the Surface-Guided Self-Assembly of Functional Molecules.

Directing the self-assembly of organic building blocks with 2D templates has been a promising method to create molecular superstructures having unique physicochemical properties. In this work the on-surface self-assembly of simple ditopic functional molecules confined inside periodic nanotemplates was modeled by means of the lattice Monte Carlo simulation method. Two types of confinement, that is honeycomb porous networks and parallel grooves of controlled diameter and width were used in the calculations. Additionally, the effect of (pro)chirality of the adsorbing molecules on the outcome of the templated self-assembly was examined. To that end, enantiopure and racemic assemblies were studied and the resulting structures were identified and classified. The obtained findings demonstrated that suitable tuning of the structural parameters of the templates enables directing the self-assembly towards linear and cyclic aggregates with controlled size. Moreover, chiral resolution of the molecular conformers using honeycomb networks with adjusted pore size was found possible. Our theoretical predictions can be helpful in designing structured surfaces to direct self-assembly and polymerization of organic functional building blocks.

Surface-Confined Self-Assembly of Asymmetric Tetratopic Molecular Building Blocks.

Surface-confined self-assembly of functional molecular building blocks has recently been widely used to create low-dimensional, also covalent, superstructures with tailorable geometry and physicochemical properties. In this contribution, using the lattice Monte Carlo simulation method, we demonstrate how the structure-property relation can be established for the 2D self-assembly of a model tetrapod molecule with reduced symmetry. To that end, a rigid functional unit comprising a few interconnected segments arranged in different tetrapod shapes was used and its self-assembly on a triangular lattice representing a (111) crystal surface was simulated. The results of our calculations show strong dependence of the structure formation on the molecular symmetry, in particular on the (pro)chiral nature of the building block. The simulations predicted the formation of unusual ordered racemic networks with unique aperiodic spatial distribution of the surface enantiomers. Molecular symmetry was also found to have significant influence on the enantiopure self-assembly which resulted in the Kagome and brickwall networks and other less ordered extended superstructures with parallelogram pores. The theoretical findings of this contribution can be relevant to designing and on-surface synthesis of molecular superstructures with predefined geometries and functions. In particular, the predicted molecular architectures can stimulate experimental efforts to fabricate and explore new nanostructures, for example graphitic, having the composition and geometry proposed in our study.

On-surface self-assembly of tetratopic molecular building blocks.

Self-assembly of functional molecules on solid substrates has recently attracted special attention as a versatile method for the fabrication of low dimensional nanostructures with tailorable properties. In this contribution, using theoretical modeling, we demonstrate how the architecture of 2D molecular assemblies can be predicted based on the individual properties of elementary building blocks at play. To that end a model star-shaped tetratopic molecule is used and its self-assembly on a (111) surface is simulated using the lattice Monte Carlo method. Several test cases are studied in which the molecule bears terminal arm centers providing interactions with differently encoded directionality. Our theoretical results show that manipulation of the interaction directions can be an effective way to direct the self-assembly towards extended periodic superstructures (2D crystals) as well as to create assemblies characterized by a lower degree of order, including glassy overlayers and quasi one-dimensional molecular connections. The obtained structures are described and classified with respect to their main geometric parameters. A small library of the tetratopic molecules and the corresponding superstructures is provided to categorize the structure-property relationship in the modeled systems. The results of our simulations can be helpful to 2D crystal engineering and surface-confined polymerization techniques as they give hints on how to functionalize tetrapod organic building blocks which would be able to create superstructures with predefined spatial organization and range of order.

Construction of Sierpinski Triangles up to the Fifth Order.

Self-similar fractal structures are of fundamental importance in science, mathematics, and aesthetics. A series of molecular defect-free Sierpinski triangle fractals have been constructed on surfaces recently. However, the highest order of the fractals is only 4 because of the limitation of kinetic growth. Here complete fifth-order Sierpinski triangles with a lateral length of 0.05 µm were successfully prepared in ultrahigh vacuum by a combination of templating and coassembly methods. Fe atoms, 4,4″-dicyano-1,1':3',1″-terphenyl, and 1,3-bis(4-pyridyl)benzene molecules were used to build fractals on the reconstructed Au(100)-(hex) surface. The new strategy may be applied to construct various Sierpinski triangles of higher orders.

Theoretical Modeling of Surface Confined Chiral Nanoporous Networks: Cruciform Molecules as Versatile Building Blocks.

Patterning of solid surfaces with functional organic molecules has been a convenient route to fabricate two-dimensional materials with programmed architecture and activities. One example is the chiral nanoporous networks that can be created via controlled self-assembly of star-shaped molecules under 2D confinement. In this contribution we use computer modeling to predict the formation of molecular networks in adsorbed overlayers comprising cruciform molecular building blocks equipped with discrete interaction centers. To that end, we employ the Monte Carlo simulation method combined with a coarse-grained representation of the adsorbed molecules which are treated as collections of interconnected segments. The interaction centers within the molecules are represented by active segments whose number and distribution are adjusted. Our particular focus is on those distributions that produce prochiral molecules able to occur in adsorbed configurations being mirror images of each other (surface enantiomers). We demonstrate that, depending on size, aspect ratio, and intramolecular distribution of active sites, the surface enantiomers can co-crystallize or segregate into extended homochiral domains with largely diversified nanosized cavities. The insights from our theoretical studies can be helpful in designing 2D chiral porous networks with potential applications in enantioselective adsorption and asymmetric heterogeneous catalysis.

Concentration-dependent supramolecular engineering of hydrogen-bonded nanostructures at surfaces: predicting self-assembly in 2D.

We report a joint computational and experimental study on the concentration-dependent self-assembly of a flat C3-symmetric molecule at surfaces. As a model system we have chosen a rigid molecular module, 1,3,5-tris(pyridine-4-ylethynyl)benzene, which can undergo self-association via hydrogen bonding (H-bonding) to form ordered 2D nanostructures. In particular, the lattice Monte Carlo method, combined with density functional calculations, was employed to explore the spontaneous supramolecular organization of this tripod-shaped molecule under surface confinement. We analyzed the stability of different weak H-bonded patterns and the influence of the concentration of the starting molecule on the 2D supramolecular packing. We found that ordered, densely packed monolayers and 2D porous networks are obtained at high and low concentrations, respectively. A concentration-dependent scanning tunneling microscopy investigation of the molecular self-assembly at a graphite-solution interface revealed supramolecular motifs, which are in perfect agreement with those obtained by simulations. Therefore, our computational approach represents a step forward toward the deterministic prediction of molecular self-assembly at surfaces and interfaces.

Supramolecular assembly of interfacial nanoporous networks with simultaneous expression of metal-organic and organic-bonding motifs.

The formation of 2D surface-confined supramolecular porous networks is scientifically and technologically appealing, notably for hosting guest species and confinement phenomena. In this study, we report a scanning tunneling microscopy (STM) study of the self-assembly of a tripod molecule specifically equipped with pyridyl functional groups to steer a simultaneous expression of lateral pyridyl-pyridyl interactions and Cu-pyridyl coordination bonds. The assembly protocols yield a new class of porous open assemblies, the formation of which is driven by multiple interactions. The tripod forms a purely porous organic network on Ag(111), phase α, in which the presence of the pyridyl groups is crucial for porosity, as confirmed by molecular dynamics and Monte Carlo simulations. Additional deposition of Cu dramatically alters this scenario. For submonolayer coverage, three different porous phases coexist (i.e., ß, γ, and δ). Phases ß and γ are chiral and exhibit a simultaneous expression of lateral pyridyl-pyridyl interactions and twofold Cu-pyridyl linkages, whereas phase δ is just stabilized by twofold Cu-pyridyl bonds. An increase in the lateral molecular coverage results in a rise in molecular pressure, which leads to the formation of a new porous phase (ε), only coexisting with phase α and stabilized by a simultaneous expression of lateral pyridyl-pyridyl interactions and threefold Cu-pyridyl bonds. Our results will open new avenues to create complex porous networks on surfaces by exploiting components specifically designed for molecular recognition through multiple interactions.

Order-Disorder Transition of Two-Dimensional Molecular Networks through a Stoichiometric Design.

Materials with disordered structures may exhibit interesting properties. Metal-organic frameworks (MOFs) are a class of hybrid materials composed of metal nodes and coordinating organic linkers. Recently, there has been growing interest in MOFs with structural disorder and the investigations of amorphous structures on surfaces. Herein, we demonstrate a bottom-up method to construct disordered molecular networks on metal surfaces by selecting two organic molecule linkers with the same symmetry but different sizes for preparing two-component samples with different stoichiometric ratios. The amorphous networks are directly imaged by scanning tunneling microscopy under ultrahigh vacuum with a submolecular resolution, allowing us to quantify its degree of disorder and other structural properties. Furthermore, we resort to molecular dynamics simulations to understand the formation of the amorphous metal-organic networks. The results may advance our understanding of the mechanism of formation of monolayer molecular networks with structural disorders, facilitating the design and exploration of amorphous MOF materials with intriguing properties.

Role of molecular orientational anisotropy in the chiral resolution of enantiomers in adsorbed overlayers.

Separation of chiral molecules using achiral inputs is an interesting alternative to traditional techniques based on the chiral recognition mechanism. In this article we propose a lattice gas Monte Carlo model of two-dimensional chiral segregation induced by breaking of molecular orientational symmetry. Simulations were performed on a square lattice for rigid chain molecules composed of four and five identical segments. Mirror-image flat chain conformations resulting in different enantiomeric pairs were considered for each probe molecule. The enantiomers were assumed to interact via short-ranged segment-segment interaction potential limited to nearest neighbors on the lattice. We considered two qualitatively different situations in which (1) the molecules were allowed to rotate on the surface and adopt any of the four planar orientations and (2) the rotation was blocked, so that only one planar orientation was possible. The results obtained for the racemic overlayers showed clearly that the orientational symmetry breaking can induce spontaneous segregation of the enantiomers into large enantiopure domains. However, this effect was observed only for molecules with sufficiently long linear fragment. In the case of kinked bulky molecules a mixed assembly was formed, demonstrating the role of molecular shape in the orientationally biased segregation of enantiomers in adsorbed films. The insights from this study can be useful in developing strategies for 2D chiral separations in which external directional fields are used.

Computer simulation of chiral nanoporous networks on solid surfaces.

A lattice Monte Carlo (MC) model was proposed with the aim of understanding the factors affecting the chiral self-assembly of tripod-shaped molecules in two dimensions. To that end a system of flat symmetric molecules adsorbed on a triangular lattice was simulated by using the canonical ensemble method. Special attention was paid to the influence of size and composition of the building block on the morphology of the adsorbed overlayer. The obtained results demonstrated a spontaneous self-assembly into extended chiral networks with hexagonal cavities, highlighting the ability of the model to reproduce basic structural features of the corresponding experimental systems. The simulated assemblies were analyzed with respect to their structural and energetic properties resulting in quantitative estimates of the unit cell parameters and mean potential energy of the adsorbed layer. The predictive potential of the model was additionally illustrated by comparison of the obtained superstructures with the recent STM images that have been recorded for different organic tripod-shaped molecules adsorbed at the liquid/pyrolytic graphite interface.

Spontaneous segregation on a hybrid chiral surface.

Segregation of enantiomers in two-dimensional adsorbed layers is a process that is usually controlled by anisotropic directional interactions between adsorbed molecules. In this contribution, we propose a simple theoretical model in which the chiral segregation occurs even though the lateral interactions are neglected. In particular, we consider a solid surface composed of two domains with different patterns of active sites being mirror images of each other. The domains of opposite handedness represent crystal facets of a composite chiral material which are adjoined to form a heterochiral adsorbing surface. To explore equilibrium properties of the system, we use Canonical Ensemble Monte Carlo method for a square lattice. The influence of factors such as energetic properties of the surface and density of the adsorbed layer on the extent of separation is examined. The obtained results indicate that effective two-dimensional separation on the hybrid chiral surface assumed in our model can be achieved only at sufficiently low adsorbate densities. The results also suggest that the segregation on the hybrid surface would be a promising method of enantiodiscrimination for those chiral molecules which do not exhibit strong lateral interactions.

Self-assembly of conformationally flexible molecules under 2D confinement: structural analysis from computer simulations.

Monte Carlo simulations were used to decipher the individual role of tripod molecular conformers in the surface-confined metal-organic self-assembly. Our calculations revealed the extent of structural heterogeneity introduced by the different conformers.

Influence of Relativistic Effects on Assembled Structures of V-Shaped Bispyridine Molecules on M(111) Surfaces Where M = Cu, Ag, Au.

The self-assembly behavior of a V-shaped bispyridine, 1,3-bi(4-pyridyl)benzene (BPyB), was studied by scanning tunneling microscopy on the (111) surfaces of Cu, Ag, and Au. BPyB molecules coordinately bonded with active Cu adatoms on Cu(111) in the form of complete polygonal rings at low coverages. On Ag(111), BPyB molecules aggregated into two-dimensional islands by relatively weak intermolecular hydrogen bonds. The coexistence of hydrogen bonds and coordination interaction was observed on the BPyB-covered Au(111) substrate. Density functional theory calculations of the metal-molecule binding energy and Monte Carlo simulations were performed to help understand the forming mechanism of molecular superstructures on the surfaces. In particular, the comprehensive orbital composition analysis interprets the observed metal-organic complexes and reveals the importance of relativistic effects for the extraordinary activity of gold adatoms. The relativistic effects cause the energy stability of the Au 6s atomic orbital and decrease the energy separation between the Au 6s and 5d orbitals. The enhanced sd hybridization strengthens the N-Au-N bond in BPyB-Au-BPyB complexes.

Theoretical investigations of the chromatographic separation of interacting enantiomers.

Separation of a pair of enantiomers by liquid chromatography is modeled using the equilibrium dispersive (ED) model of chromatography. It is assumed that the chiral stationary phase used for the separation consists of two types of adsorption sites, including chiral selectors linked to the surface and nonselective centers belonging to the achiral matrix. Additionally, intermolecular interactions between adsorbed enantiomers are taken into account. The corresponding equilibrium adsorption isotherms of the enantiomers are derived by means of the mean field approximation (MFA) and used as input data for the ED model. Special attention is paid to the influence of the lateral interactions on the effectiveness of the enantiomer separation. In particular, we examine the effect of the interactions on the shape and relative position of the chromatographic peaks associated with the enantiomers. Furthermore, the influence of the spacer length, which modifies screening of the lateral interactions, on the adsorption process is studied. The obtained results suggest that the lateral interactions combined with the screening effect may cause serious changes in the separation, depending on the nature (attraction or repulsion) and strength of the interactions as well as on the spacer length.

On-surface synthesis of two-dimensional imine polymers with a tunable band gap: a combined STM, DFT and Monte Carlo investigation.

Two-dimensional polymers are of great interest for many potential applications in nanotechnology. The preparation of crystalline 2D polymers with a tunable band gap is critical for their applications in nano-electronics and optoelectronics. In this work, we try to tune the band gap of 2D imine polymers by expanding the conjugation of the backbone of aromatic diamines both laterally and longitudinally. STM characterization reveals that the regularity of the 2D polymers can be affected by the existence of lateral bulky groups. Density functional theory (DFT) simulations discovered a significant narrowing of the band gap of imine 2D polymers upon the expansion of the conjugation of the monomer backbone, which has been confirmed experimentally by UV absorption measurements. Monte Carlo simulations help us to gain further insight into the controlling factors of the formation of regular 2D polymers, which demonstrated that based on the all rigid assumption, the coexistence of different conformations of the imine moiety has a significant effect on the regularity of the imine 2D polymers.