# How do the orbitals in atoms synchronize their phases during the formation of a bond?

Suppose I have two hydrogen atoms, each with one electron in their 1s orbital. Then we can represent the wave functions of these two as say $$|\Psi\rangle_1$$ and $$|\Psi\rangle_2$$ for each of the two atoms respectively. Since we have two identical atoms, $$|\Psi\rangle_2$$ is just $$|\Psi\rangle_1$$ translated to a different point in space with some phase difference, so $$|\Psi\rangle_2 = e^{i\theta}|\Psi\rangle_1'$$ where the $$'$$ denotes that we've shifted the wavefunction in space.

We know that if the two atoms have wave functions in phase with each other ($$\theta = 0$$) then when they come together the orbitals will combine with each other and a $$\sigma$$-bond will be formed. Equally if the wave functions are completely out of phase with each other then the anti-bonding orbital of the combined system would have 2 electrons while the bonding orbital would have 0 and since the anti-bonding orbital has higher energy this transformation is not energetically favourable so no reaction occurs.

I'm curious about what happens when the phase is not equal to either $$0$$ or $$\pi$$, which is a case that happens with probability $$1$$. Clearly the atoms do end up reacting in real life so the phases must get synchronized somehow before the reaction takes place. I'd like to ask how exactly this happens.

• Does this answer your question? Significance of phase of atomic orbitals Commented Jul 10, 2023 at 19:05
• Not really, the answer on that question (or the linked one from that answer) doesn't deal with the formation step. Commented Jul 10, 2023 at 19:32
• If you're really interested in how H2 forms, than you're missing the point. Molecule is initially in excited state and needs to lose surplus energy quickly. Commented Jul 10, 2023 at 19:36
• I'm not interested in how H2 forms, say I have molecule 1 in phase 0 and molecule 2 in phase $\pi/2$. This means that the wavefunction of molecule 2 has an $i$ in front of it so when you add the wavefunctions together you'd get no interference at all (the first one acts in the real dimension, the second one acts in the imaginary dimension). I'm interested in how exactly the phases get synchronized so the reaction can take place. Commented Jul 10, 2023 at 19:39
• In the simplest case when two atoms approach one another their wavefunctions are changed due to mutual interaction, (nuclear-nuclear, electron-electron and nuclear-electron) so phase of isolated and non-interacting wavefunctions no longer matters in the way you suggest. To solve this problem the coordinate system is changed to, say, an elliptical one with each atom at the foci. The integrals of these interactions are then solved as the internuclear distance changes. Commented Jul 11, 2023 at 15:28

The expression $$|\Psi\rangle_2 = e^{i\theta}|\Psi\rangle_1'$$, for $$\theta \neq k \pi$$ ($$k \in \mathbf{Z}$$) is likely to reach to a complex function. This situation is unlikely under an attracting potential because the position tends to be more precise (requires a real wave function) at the expense of the momentum (which requires a complex wave function), even if the exact wavefunction is unknown.
On the other hand if we consider a translation operator $$\hat{T}$$ acting on a $$s$$ orbital of one hydrogen atom of a site $$i$$ transferring the orbital to another site $$j$$ in a chain of $$n > 0$$ hydrogen atoms having $$a$$ as interatomic distance in a $$x$$ axis $$\mathbf{R} = l a \hat{e_x}$$ with $$l=1,2,...,n$$. The case $$n = 2$$ is the hydrogen molecule $$H_2$$. Generally, $$\hat{T}$$ and the Hamiltonian $$\hat{H}$$ do not commute : $$[\hat{H},\hat{T}] \neq 0$$ this means it is "unlikely" that they share the same eigenvectors because of the repulsion $$V_{ee}=+\sum_{ij}\frac{1}{|r_i-r_j|}$$.
We can consider the case $$V_{ee}$$ can be ignored so that $$[\hat{H},\hat{T}] = 0$$, the operator $$\hat{T}$$ gives $$\hat{T}|\Psi\rangle_1 = e^{i\mathbf{k}\mathbf{R}} |\Psi\rangle_1$$, if $$k=\frac{2\pi}{a}$$, therefore $$\hat{T}|\Psi\rangle_1 = e^{i 2\pi l} |\Psi\rangle_1 = |\Psi\rangle_2$$. Here the phase shift is due to the translation operator and can only be a positive integer (1). A similar result compared to the expression above, but there is no anti-bonding.
This result is due to the fact that $$\hat{T}$$ is a generator of bonding orbitals, transferring the wave function from one atom to another. The state described above is similar to a singlet state, where $$V_{ee}$$ is ignored because there is an exchange interaction (attractive) $$V_{x}$$ which lowers significantly the impact of $$V_{ee}$$, for $$n=2$$, $$\Psi_B \propto (|\Psi\rangle_1 + |\Psi\rangle_2)$$. Now for the antibonding is due to the fact that $$V_{x}$$ is only a pairwise interaction, it means the exchange do not act on more than two states, this anti-bonding state has the impact of $$V_{ee}$$ and $$[\hat{H},\hat{T}] \neq 0$$, the state is not an eigenvectors of $$\hat{T}$$, the only possible choice is $$\Psi_A \propto (|\Psi\rangle_1 - |\Psi\rangle_2)$$.