We present classical density functional theory calculations of the free energy landscape for fluids below their triple point as a function of density and crystallinity. We find that for both a model globular protein and for a simple atomic fluid modeled with a Lennard-Jones interaction, it is free-energetically easier to crystallize by passing through a metastable dense fluid in accord with the Ostwald rule of stages but in contrast to the alternative of ordering and densifying at once as assumed in the classical picture of crystallization.PACS numbers: 82.60. Nh,87.15.Nn,05.20.Jj Crystallization is an intricate process of fundamental importance in many areas of physics, chemistry and engineering. The classical picture of crystallization from supersaturated solutions goes back to Gibbs and consists of the spontaneous formation of crystalline clusters which then either grow or shrink depending on the relative importance of the free energy gain due to the lower bulk free energy of the crystal cluster and the free energy penalty due to the surface tension between the two phases. In this picture, the local density is the only order parameter: the crystalline cluster is (in general) denser than the fluid. In recent years, this picture has been called into question by simulation, theory and experiment for the particular and important case of the crystallization of globular proteins. ten Wolde and Frenkel (hereafter tWF) showed by means of simulation that the free energy landscape of protein crystal clusters as a function of the number of atoms in the cluster and the "crystallinity " favored paths leading from no clusters to clusters with low order to ordered clusters over paths moving from no clusters directly to ordered clusters [1]. This picture was confirmed by Talanquer and Oxtoby [2] and Shiryayev and Gunton[3] who showed using a parameterized van der Waals-type model of globular proteins that surface wetting did indeed lower the free energy of crystal clusters. More recently, the simple picture has also been challenged by novel experimental investigations. Vekilov and co-workers have shown that, prior to crystallization, protein solutions harbor metastable droplets of dense fluid and they have suggested that these droplets are necessary precursors of crystallization [4,5,6,7]. The picture that emerges is one of the formation of metastable droplets of dense fluid which then subsequently crystallize. In this paper, we show by means of classical density functional theory calculations that there is an intrinsic free-energy advantage in first densifying into a metastable densefluid state and then crystallizing rather than following the classical path which goes directly from gas to crystal. Furthermore, our calculations suggest that a similar advantage exists for fluids of small molecules, modeled here via the Lennard-Jones (LJ) interaction, thus indicating that this mechanism may underlie most crystallization processes.The starting point for our analysis is classical density functional theory (DFT) which is based...