We investigate the Liouvillian integrability of Hamiltonian systems describing a universe filled with a scalar field (possibly complex). The tool used is the differential Galois group approach, as introduced by Morales-Ruiz and Ramis. The main result is that the generic systems with minimal coupling are non-integrable, although there still exist some values of parameters for which integrability remains undecided; the conformally coupled systems are only integrable in four known cases. We also draw a connection with the chaos present in such cosmological models, and the issues of the integrability restricted to the real domain.fields [23], as has chaotic dynamics [48].The first of our results is that minimally coupled fields are not integrable in the generic case. There are however special families of the system's parameters which leave the question open. We give the appropriate conditions in the concluding section.There are several physical reasons to study more than just minimally coupled fields. Early works on chaotic inflation found the coupling constant ξ small or negative [30] but some argue [26] that the paradigm of inflation should be generalised to the case with non-zero coupling constant which should not be fine-tuned close to zero, and WMAP observations seem to indicate non-negligible ξ.The coupling could be generated by quantum corrections [10,29], or from the renormalisation of the Klein-Gordon equation as described in [16]. The coupling constant should be fixed by particle physics of the matter composing the scalar field, for example the way ξ = 1/6 was found in the large N approximation to the Nambu-Jona-Lasimo model in [34]. Non-minimally coupled fields are also interesting in the context of description of the dark energy for which the ratio between the pressure and the energy density is less than −1. Such matter is called a phantom matter, and cannot be achieved by standard scalar fields [27].Conformally coupled fields were subject to more rigorous integrability analysis, as opposed to minimally coupled ones, thanks to the natural form of their Hamiltonian. As will be shown in the next section, the kinetic part is of natural form, albeit indefinite, and the potential is polynomial (in the case of real fields).Chaos has been studied in such fields by means of Lyapunov exponents, perturbative approach, breaking up of the KAM tori [17,11]. Also the Painlevé property [35] was employed as an indicator of the system's integrability.Ziglin proved that the system given by (20) is not meromorphically integrable when Λ = λ = 0 and k = 1 [67]. His methods were also used by Yoshida to homogeneous potentials which is the case for the system when k = 0 [61,62,63,64]. Later, Yoshida's results were sharpened by Morales-Ruiz and Ramis [49], and used by the present authors in [44] to obtain countable families of possibly integrable cases with some restrictions on λ and Λ. Also recently, more conditions for integrability have been given in [12], although only for a non-zero spatial curvature k and a generic value of...