Using atomic scale time-dependent density functional calculations we confirm that both diffusioncontrolled and diffusionless crystallization modes exist in simple 2D systems. We provide theoretical evidence that a faceted to nonfaceted transition is coupled to these crystallization modes, and faceting is governed by the local supersaturation at the fluid-crystalline interface. We also show that competing modes of crystallization have a major influence on mesopattern formation. Irregularly branched and porous structures are emerging at the crossover of the crystallization modes. The proposed branching mechanism differs essentially from dendritic fingering driven by diffusive instability. [6,7]. The amazing complexity of faceted patterns has been generating considerable interest for decades. Various crystal shapes; compact symmetric, dendritic, ramified fractal and needlelike morphologies have been investigated experimentally [1][2][3][4][5]. While robust theoretical interpretations of faceted equilibrium shapes exist, the dynamical aspects of faceted pattern formation are less understood. The main difficulty in establishing a theoretical description lies in the multiscale nature of faceted crystal growth, the mesoscale behavior is influenced by the local arrangement of the particles. Russel et al. [8] have identified some mesoscale characteristics of crystal growth kinetics: (i) a diffusion-controlled (or ''slow'') growth mode that may lead to diffusive instability and dendritic branching in monodisperse colloids [8,9], and (ii) a diffusionless (or ''fast'') steady growth mode, in which crystalline ordering takes place without significant density change. The existence of these modes is verified experimentally by investigating crystallization kinetics in hard-sphere suspensions [10,11]. Besides, in some organic glasses discontinuous enhancement of crystal growth is observed [12][13][14], which is associated with diffusionless crystallization [15]. Although various theories have been developed to explain the existence of diffusion-controlled and diffusionless rapid modes [14][15][16], none of them predicts a facetednonfaceted morphological transition often associated with different crystallization modes, when interpreting experiments [3]. Growth mode (i) and mode (ii) have also been observed in simple mean-field models [17], although these theories address anisotropic crystal growth in a phenomenological way. In contrast, molecular theories based on the classical density functional technique (DFT) predict not only the anisotropic bulk crystal properties, such as the crystal structure and elastic constants, but also anisotropic surface properties (e.g., surface tension [18]). A simple dynamical DFT [19], the phase-field crystal (PFC) model [20] is able to address crystallization kinetics up to the mesoscale [21]. While the phase-field crystal model represents the average of the microscopic states over a coarse graining time, it does describe crystal defects. Grain boundaries and dislocations can be observed directl...