Densities of the quaternary system consisting of {1-chlorobutane (1) + butylamine (2) + butylacetate (3) + isobutanol (4)}, related ternary systems of {1-chlorobutane (1) + butylamine (2) + butylacetate (3)}, (1-chlorobutane (1) + butylacetate (3) + isobutanol (4)}, and {butylamine (2) + butylacetate (3) + isobutanol (4)}, and binary systems of {1-chlorobutane (1) + butylacetate (3)} and {butylamine (2) + butylacetate (3)} were measured over the whole range of composition at 298.15 K and ambient pressure. Excess molar volumes, V(m)(E), for the mixtures were derived and correlated as a function of mole fraction by using the Redlich-Kister and the Cibulka equations for binary and ternary mixtures, respectively. From the experimental data, partial molar volumes, V(m,i), and excess partial molar volumes, V(m,i)(E), were also calculated for binary systems. The experimental results of the constituted binary mixtures have been used to test the applicability of the Prigogine-Flory-Paterson (PFP) theory. Values of densities and excess molar volumes, V(m)(E) , of other related binary and ternary systems were reported in our previous works. Some geometrical models which are the Kohler, Muggianu, Toop, Li, and a new model by Zhang as symmetric and asymmetric models have also been extended to estimate the physicochemical properties of quaternary systems based on the data of four subternary and six sub-binary systems. The experimental data were used to evaluate the nature and type of intermolecular interactions in multicomponent systems.