Liquid−liquid phase separation (LLPS) of transactive response DNAbinding protein of 43 kDa (TDP-43), which exerts multiple functions in the splicing, trafficking, and stabilization of RNA, mediates the formation of membraneless condensates with crucial physiological roles, while its aberrant LLPS is linked to multiple neurodegenerative diseases. However, due to the heterogeneous and dynamic nature of LLPS, major gaps remain in understanding the precise intermolecular interactions driving LLPS and how specific mutations alter LLPS dynamics. Here, we investigated the molecular mechanisms underlying the LLPS of the TDP-43 lowcomplexity domain (LCD) by simulating the dimerization process using all-atom discrete molecular dynamics with microsecond-long simulations. Our results showed that the TDP-43 LCD was intrinsically disordered, with helical structures consistent with prior nuclear magnetic resonance studies. Phase separation propensity was assessed by simulating the dimerization of the TDP-43 LCD and four mutants, showing that A321G, W334G, and M337V inhibited self-association, while G335D promoted it, fully consistent with experimental reports. During the dimerization process, two peptides experienced both elastic and nonelastic collisions, and the self-associated dimer featured both high-and low-contact states. These results suggested that the dimerization process of the TDP-43 LCD was accordingly dynamic and heterogeneous. Additionally, we identified crucial regions containing hydrophobic clusters and aromatic residues in the Nterminus, central region, and C-terminus that were essential for the self-association of the TDP-43 LCD. These residues with high binding affinities can act as stickers to form peptide networks in LLPS. Together, our simulation provides a comprehensive picture of the intermolecular interactions driving the phase separation of the TDP-43 LCD, offering insights into both physiological functions and pathological mechanisms.