Attractive interactions between molecules lead to formation of the liquid phase at sufficiently low temperatures. In the absence of external fields or containers, liquids assume the shape of spherical drops, which minimize their surface area. In systems with 3-D periodic boundary conditions (PBCs), which are routinely employed in molecular simulations, alternate configurations can be stable depending on the extent of the system. The stability of a variety of structures under 3-D PBCs originates from the underlying variation of free energy with density as shown elegantly by MacDowell and co-workers (MacDowell, L. G.; Shen, V. K.; Errington, J. R. J. Chem. Phys. 2006, 125, 3.). Here we present analysis of extensive Monte Carlo and molecular dynamics simulations of Lennard-Jones and water fluids to calculate free energy and explore the phase diagram that governs formation of different liquid assemblies. We also study metastability of different shapes and their interconversions by systematically initializing simulations in various configurations. Further, We present results on the rupture of thin liquid films on solid substrates, with focus on the evolution of liquid structure and the rupture mechanism. Our estimates of important capillary wavelengths from simulations are in good agreement with theoretical predictions of Vrij and Overbeek (Vrij, A.; Overbeek, J. T. J. Am. Chem. Soc. 1968, 90, 3074−3078.). Collectively, our work significantly extends the previous simulation studies of interfacial systems, and especially of thin-film structure, stability, and rupture processes in molecular simulations.