Hybridisation and the Schrödinger equation I am slightly confused about hybridisation and how it relates to molecular and atomic orbitals, despite having pored through many sources online. I was hoping someone could verify whether my current understanding is correct, in particular regarding what hybridisation actually is/does because I have not read this explicitly, but am assuming it is the case from what I have read so far:


The Schrödinger equation can be solved to give the atomic orbitals (at least, some simplification involving the effective nuclear charge can be used to find the outer (and inner?) atomic orbitals.
In hybridisation, we first consider the shape of a molecule and then we consider each atom seperately and how we could combine the atomic orbitals so that the geometry about each atom is correct. Now what I am particularly unsure about is: Are the hybridised orbitals entirely a conceptual construct, or are they mathematical solutions to the Schrödinger equation? I was thinking that maybe the hybrid atomic orbitals that are created are a linear combination of the atomic orbitals, and thus this linear combination also solved the Schrödinger equation and can exist but would have a different energy to the individual atomic orbitals? I am finding it difficult to think how a linear combination of the atomic orbitals could produce something with a completely different shape though (although I suppose the orbitals are orientated and essentially vectorial, so I think it could work). Perhaps  the hybrid atomic orbital form could also be, say, the product of the atomic orbital forms? Although then I do not think this would solve the Schrödinger equation anymore. But I am not sure. And anyway, all of these ideas assume that hybridisation has a mathematical basis in the Schrödinger equation, which I am not sure about at all!
Then molecular orbital theory takes the orbitals of each atom (I think theoretically it takes every orbital in every atom?) in the molecule and combines them to form a molecular orbital for the whole molecule - no longer simply an overlap between two atoms. This certainly also solves the Schrödinger equation. From what I understand, theoretically the calculation to form MOs should involve only the AOs of each atom, but to simplify the situation sometimes hybrid atomic orbitals are used in MO construction? If so, surely the hybrid atomic orbitals must be a linear combination of the atomic orbitals that solve the Schrödinger equation, so that this can be used in molecular obrital theory which is based in the solution to the Schrödinger equation? Also, where in the calulcation are the coefficients for each atomic orbital arising? Is it from the boundary conditions including things like molecule geometry, or perhaps that we know the energy of a molecule and when we put this particular energy in to the Schrödinger equation it gives us the coefficients?


I apologise for the long post and I realise there seem to be many questions within here, however they are all linked and about the mathematical underpinnings of hybridisation/valence bond theory and molecular orbitals, so I thought they belonged in one post.
Mathematically, atomic/molecular orbitals are 1-electron wavefunctions (hydrogen-like wavefunctions) that are used as a basis with-which the total N-election wavefunction is expanded.  The N-election wavefunction is a determinant (or a linear combination of determinants) built from these 1-electron wavefunctions: an anti-symmetric (Fermi statistics) linear combination of Hartree products (usually a product of MOs expanded in a basis of AOs).

For simplicity's sake, lets assume we are interested only in 1-determinant wavefunctions.  The Hartree-Fock method will (when used where HF is appropriate) result in a variationally-optimal 1-determinant wavefunction that is an energy eigenfunction of the applicable Hamiltonian operator.  Not only is the final, self-consistently optimized HF total N-electron wavefunction an energy eigenfunction, but the individual HF canoncial MOs are also eigenfunctions of the Hamiltonian/Fock operator, with eigenvalues that are the so-called 'orbital energies' (canonical orbitals are defined as the optimzed orbitals which diagonalized the Fock matrix).

However, there is nothing special about these canonical orbitals in the context of the total system energy.  Any linear, norm-preserving transformation of the occupied orbital space will result in an N-electron wavefunction that is still an eigenfunction of the Fock operator, while the individual MOs are no longer eigenfunctions of the Fock operator (for example, localized MOs, LMOs).

So, once we have a HF wavefunction and optimized MOs, we can linearly transform the MOs however we want.  Such a linear transformation might have the effect of producing orbitals that resemble hybridized orbitals.  Or maybe we are interested in generating localized MOs.  The point is, a linear transformation within the occupied orbital space results in a wavefunction that has the exact same electon density and energy eigenvalue as the wavefunction HF originally provides us in terms of Canonical orbitals.

Orbitals transformed to look like hybrid-orbitals are on the exact same footing as any other choice of orbital.