We describe a complete set of algorithms for ab initio molecular simulations based on numerically tabulated atom-centered orbitals (NAOs) to capture a wide range of molecular and materials properties from quantum-mechanical first principles. The full algorithmic framework described here is embodied in the Fritz Haber Institute "ab initio molecular simulations" (FHI-aims) computer program package. Its comprehensive description should be relevant to any other first-principles implementation based on NAOs. The focus here is on density-functional theory (DFT) in the local and semilocal (generalized gradient) approximations, but an extension to hybrid functionals, Hartree–Fock theory, and MP2/GW electron self-energies for total energies and excited states is possible within the same underlying algorithms. An all-electron/full-potential treatment that is both computationally efficient and accurate is achieved for periodic and cluster geometries on equal footing, including relaxation and ab initio molecular dynamics. We demonstrate the construction of transferable, hierarchical basis sets, allowing the calculation to range from qualitative tight-binding like accuracy to meV-level total energy convergence with the basis set. Since all basis functions are strictly localized, the otherwise computationally dominant grid-based operations scale as O(N) with system size N. Together with a scalar-relativistic treatment, the basis sets provide access to all elements from light to heavy. Both low-communication parallelization of all real-space grid based algorithms and a ScaLapack-based, customized handling of the linear algebra for all matrix operations are possible, guaranteeing efficient scaling (CPU time and memory) up to massively parallel computer systems with thousands of CPUs

We computationally study the reaction mechanisms of halogen-based covalent self-assembly, a major route for synthesizing molecular nanostructures and nanographenes on surfaces. Focusing on biphenyl as a small model system, we describe the dehalogenation, recombination, and diffusion processes. The kinetics of the different processes are also investigated, in particular how diffusion and coupling barriers affect recombination rates. Trends across the periodic table are derived from three commonly used close-packed ( 111) surfaces (Cu, Ag, and Au) and two halogens (Br and I). We show that the halogen atoms can poison the surface, thus hindering long-range ordering of the self-assembled structures. Finally, we present core-level shifts of the relevant carbon and halogen atoms, to provide reference data for reliably detecting selfassembly without the need for atomic-resolution scanning tunneling microscopy.

Although the local environment of a molecule can play an important role in its chemistry, rarely has it been examined experimentally at the level of individual molecules. Here we report the precise control of intramolecular hydrogen-transfer (tautomerization) reactions in single molecules using scanning tunnelling microscopy. By placing, with atomic precision, a copper adatom close to a porphycene molecule, we found that the tautomerization rates could be tuned up and down in a controlled fashion, surprisingly also at rather large separations. Furthermore, we extended our study to molecular assemblies in which even the arrangement of the pyrrolic hydrogen atoms in the neighbouring molecule influences the tautomerization reaction in a given porphycene, with positive and negative cooperativity effects. Our results highlight the importance of controlling the environment of molecules with atomic precision and demonstrate the potential to regulate processes that occur in a single molecule.

The adsorption of neutral (poly)-aromatic, antiaromatic, and more generally π-conjugated systems on graphene is studied as a prototypical case of π−π stacking. To account for dispersive interactions, we compare the recent van der Waals density functional (vdw-DF) with three semiempirical corrections to density functional theory and two empirical force fields. The adsorption energies of the molecules binding to graphene predicted by the vdw-DF were found to be in excellent agreement with temperature desorption experiments reported in literature, whereas the results of the remaining functionals and force fields only preserve the correct trends. The comparison of the dispersive versus electrostatic contributions to the total binding energies in the aromatic and antiaromatic systems suggests that π−π interactions can be regarded as being prevalently dispersive in nature at large separations, whereas close to the equilibrium bonding distance, it is a complex interplay between dispersive and electrostatic Coulombic interactions. Moreover our results surprisingly indicate that the magnitude of π−π interactions normalized both per number of total atoms and carbon atoms increases significantly with the relative number of hydrogen atoms in the studied systems.

Individual pentacene and naphthalocyanine molecules adsorbed on a bilayer of NaCl grown on Cu(111) were investigated by means of scanning tunneling microscopy using CO-functionalized tips. The images of the frontier molecular orbitals show an increased lateral resolution compared with those of the bare tip and reflect the modulus squared of the lateral gradient of the wave functions. The contrast is explained by tunneling through the p-wave orbitals of the CO molecule. Comparison with calculations using a Tersoff-Hamann approach, including s- and p-wave tip states, demonstrates the significant contribution of p-wave tip states.