It should be noted that the s,p,d... etc. orbitals you refer to are solutions to the one-electron Schroedinger eqn. only. For atoms and molecules with more electrons, multi-electron systems, no exact closed solutions are known. Nonetheless, multi-electron systems can be approximated using linear combinations of hydrogenic orbitals (LCAO). Moreover, many properties of atomic and molecular systems can be explained, if not exactly predicted, using either AO's or LCAO's. Some of the strongest evidence of the fundamental role of atomic orbitals may be found in VSEPR theory. This qualitative theory predicts the structure of molecules of all kinds by considering repulsions between electron pairs occupying valence shell orbitals and hybrids. (This theory is normally described in any modern first year college chemisty text.) So, the existence of atomic orbitals in multi-electron systems has proven quite useful for a wide range of predictions, even if their existence hasn't been mathematically proved.
Regards,
Jim Ritchie