Scalar field ϕCDM models provide an alternative to the standard ΛCDM paradigm, while being physically better motivated. Dynamical scalar field ϕCDM models are divided into two classes: the quintessence (minimally and non-minimally interacting with gravity) and phantom models. These models explain the phenomenology of late-time dark energy. In these models, energy density and pressure are time-dependent functions under the assumption that the scalar field is described by the ideal barotropic fluid model. As a consequence of this, the equation of state parameter of the ϕCDM models is also a time-dependent function. The interaction between dark energy and dark matter, namely their transformation into each other, is considered in the interacting dark energy models. The evolution of the universe from the inflationary epoch to the present dark energy epoch is investigated in quintessential inflation models, in which a single scalar field plays a role of both the inflaton field at the inflationary epoch and of the quintessence scalar field at the present epoch. We start with an overview of the motivation behind these classes of models, the basic mathematical formalism, and the different classes of models. We then present a compilation of recent results of applying different observational probes to constraining ϕCDM model parameters. Over the last two decades, the precision of observational data has increased immensely, leading to ever tighter constraints. A combination of the recent measurements favors the spatially flat ΛCDM model but a large class of ϕCDM models is still not ruled out.