Calculating the spectroscopic properties of complex conjugated organic molecules in their relaxed state is far from simple. An additional complexity arises for flexible molecules in solution, where the rotational energy barriers are low enough so that nonminimum conformations may become dynamically populated. These metastable conformations quickly relax during the minimization procedures preliminary to density functional theory calculations, and so accounting for their contribution to the experimentally observed properties is problematic. We describe a strategy for stabilizing these nonminimum conformations in silico, allowing their properties to be calculated. Diadinoxanthin and alloxanthin present atypical vibrational properties in solution, indicating the presence of several conformations. Performing energy calculations in vacuo and polarizable continuum model calculations in different solvents, we found three different conformations with values for the δ dihedral angle of the end ring ca. 0, 180, and 90° with respect to the plane of the conjugated chain. The latter conformation, a nonglobal minimum, is not stable during the minimization necessary for modeling its spectroscopic properties. To circumvent this classical problem, we used a Car-Parinello MD supermolecular approach, in which diadinoxanthin was solvated by water molecules so that metastable conformations were stabilized by hydrogen-bonding interactions. We progressively removed the number of solvating waters to find the minimum required for this stabilization. This strategy represents the first modeling of a carotenoid in a distorted conformation and provides an accurate interpretation of the experimental data.