Abstract: Among all types of discovered crystals, those formed by organic molecules show the greatest diversity, which results from the intrinsic complexity of the organic molecules and the weak interactions between them. Even for a given compound, different crystal structures can exist. This feature is referred to as polymorphism in the modern crystallographic context and those different crystal forms are called polymorphs. In reality, the crystallization of organic molecules is often performed at the surface of a substrate, giving rise to heterogeneous crystallization. Except for the well-known catalyzing effects, the existence of substrates brings more possibilities to the polymorphic behaviors of organic molecules, promoting the formation of new polymorphs that are only stable in the vicinity of the substrates. For this reason, these new polymorphic forms are often described as substrate-induced polymorphs (SIPs). It is of great importance to understand the formation of SIPs for organic molecules as it has been reported that SIPs can show superior properties with respect to their bulk form counterparts. Up to now, most studies focus on the identifying and characterizing the presence of SIPs, which relies mainly on X-ray diffraction techniques. However, a detailed explanation about the origin of SIPs is still missing. In this work, we have combined several powerful experimental characterization techniques, including X-ray diffraction, transmission electron microscopy (TEM) and scanning tunneling microscopy (STM) in order to reach an integrated view over the formation of SIPs. These experimental studies are strongly supported by computational chemistry simulations, such as density functional theory and molecular dynamics. A big advantage of using atomistic simulations is that it enables the possibility to predict a priori the crystal structures of SIPs and to establish a posteriori the general rules for the formation of SIPs. In practice, this thesis employs state-of-art atomistic simulation approaches in order to bridge substrate-induced polymorphism with a conceptually-connected research area: the self-assembly of molecular networks (SAMNs), also called 2D crystallization. Unlike SIPs, which extend at least several molecular layers, SAMNs are composed of a single layer of molecules with ordered packing. Our simulations have enabled a more comprehensive understanding about the role of substrate during the formation of SIPs and we elucidate how the positional and orientational order of molecules propagates from the substrate to the upper 2D and even 3D crystal layers. In this way, a fundamental understanding of the substrate-induced crystallization is gained by connecting 2D and 3D crystallization using substrate-induced approaches.