In this work geometric optimization and potential energy surface (PES) scan for the TiO2 molecule were carried out to find out possible stable structures by using quantum chemical calculations. The ground state (S-0 1 (1)A(1)) has the only equilibrium with a symmetric bent structure, whereas the linear conformer is unstable. The first singlet excited state (S-1 1 B-1(2)), on the other hand, possesses two potential minima where one is bent around the Franck-Condon region and the other is linear. The second singlet excited state (S-2 1 (1)A(2)) has only one minimum which is linear, and at which geometry the first two excited states degenerate. Other low-lying singlet and triplet excited states have been investigated accordingly. Among the functionals applied in the present DFT calculations, B3LYP gave results most close to available experimental data, while some recently developed ones including HSE06, omega B97X-D and CAM-B3LYP were not satisfactory for this small system. Ab initio methods such as CASSCF, CASPT2 as well as CCSD-related calculations have been applied, and the latter two showed good results similar to DFT/B3LYP.