  A strategy is presented for variational orbital optimization in
time-independent calculations of excited electronic states. The approach
involves minimizing the energy while constraining specific degrees of freedom
corresponding to negative curvature on the electronic energy surface, followed
by a fully unconstrained saddle point search. Both steps of this
freeze-and-release strategy are carried out via direct orbital optimization at
similar cost as ground state calculations. The approach is applied in orbital
optimized density functional calculations of intramolecular charge transfer
excited states where the maximum overlap method proves unable to prevent
variational collapse to unphysical, charge delocalized solutions. Calculations
with the local density approximation and generalized gradient approximation
functionals are carried out for a large set of charge transfer excitations in
organic molecules and the results compared to linear-response adiabatic
time-dependent density functional theory (TD-DFT) calculations. The
time-independent approach is fully variational and provides a relaxed excited
state electron density that can be used to quantify the extent of charge
transfer. The TD-DFT calculations are found to generally overestimate the
charge transfer distance, even when including first-order orbital relaxation
effects with the Z-vector method in the so-called ''relaxed'' TD-DFT excited
state density. Furthermore, while the error on the excitation energy relative
to theoretical best estimates is found to increase with the extent of charge
transfer for TD-DFT, no correlation is observed for the orbital optimized
calculations, leading to a significant improvement for the excitations with
strong charge transfer character.
