Intrinsically disordered proteins (IDPs) are essential components for the formation of membraneless organelles, which play key functional and regulatory roles within biological systems. These complex assemblies form and dissolve spontaneously over time via liquid-liquid phase separation of IDPs. Mutations in their amino acid sequence can alter their phase behavior, which has been linked to the emergence of severe diseases. We study the conformation and phase behavior of a low-complexity domain of heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1), using coarse-grained implicit solvent molecular dynamics simulations. We systematically analyze how these properties are affected by the number of aromatic residues within the examined sequences. We find a significant compaction of the chains and an increase in the critical temperature with increasing number of aromatic residues. The local persistence length is determined in single-chain simulations, revealing strong sequence-specific variations along the chain contour. Comparing single-chain and condensed state simulations, we find much more collapsed polymer conformations in the dilute systems, even at temperatures near the estimated theta temperature of the solution. These observations strongly support the hypothesis that aromatic residues play a dominant role for condensation, which is further corroborated by a detailed analysis of the intermolecular contacts, and conversely that important properties of condensates are captured in coarse-grained simulations. Interestingly, we observe density inhomogeneities within the condensates near criticality, which are driven by electrostatic interactions. Finally, we find that the relatively small fraction of hydrophobic residues in the IDPs results in interfacial tensions which are significantly lower compared to typical combinations of immiscible simple liquids.