Cyclopentadienyl complexes of barium have great utility in materials science and engineering, in particular, as precursors in the atomic layer deposition processes, which are required to be fluidic as well as thermally stable and volatile. Here, we investigated the liquid-like properties of cyclopentadienyl barium complexes including (Me5C5)2Ba, (tBu3C5H2)2Ba, (iPr4C5H)2Ba, (iPr5C5)2Ba, and [(SiMe3)3C5H2]2Ba, using molecular dynamics simulations of nanoscale droplets. The compounds were modeled using a recently developed generic force field, GFN-FF. Nanoscale droplets with about 5.0 nm diameters were formed by aggregating 96 molecules of each compound. Simulation results reveal that substituting methyl groups of (Me5C5)2Ba with other alkyl and silyl moieties has a non-negligible effect on the intra- and intermolecular structure and dynamics. In particular, in contrast to more flexible (Me5C5)2Ba, the substitution with five iso-propyl groups to form (iPr5C5)2Ba adds rigidity to the complex with restricted orientational fluctuations for two cyclopentadienyl ligands and arranges molecules parallel to each other with greater probability. In addition, comparison between (tBu3C5H2)2Ba, with three tert-butyl groups, and its silyl analogue, [(SiMe3)3C5H2]2Ba, reveals that intermolecular interactions between the molecules with silyl groups are softer than those with tert-butyl groups and result in broader radial distribution functions, whereas the dynamic properties are similar for both compounds. This work suggests that molecular dynamics simulations contribute to molecular-level understanding of the effect of chemical substitution in organometallic compounds on the intra- and intermolecular properties of molecular liquids.