The homogeneous cosmological models with a Liouville scalar field are investigated in classical and quantum context of Wheeler-DeWitt geometrodynamics. In the quantum case of quintessence field with potential unbounded from below and phantom field, the energy density operators are not essentially self-adjoint and self-adjoint extensions contain ambiguities. Therefore the same classical actions correspond to a family of distinct quantum models. For the phantom field the energy spectrum happens to be discrete. The probability conservation and appropriate classical limit can be achieved with a certain restriction of the functional class. The appropriately localized wave packets are studied numerically using the Schrödinger's norm and a conserved Mostafazadeh's norm introduced from techniques of pseudo-Hermitian quantum mechanics. These norms give a similar packet evolution that is confronted with analytical classical solutions.