If to particles with wave functions $\psi_a$ and $\psi_b$ are „combined“ their total wave function is given by:
$$\psi(r_1,r_2)= A[\psi_a(r_1)\psi_b(r_2) \pm \psi_b(r_1)\psi_a(r_2)]$$
(+ for bosons and - for fermions, $A$ is a constant, here I‘m ignoring spin.)
But if we combine the wave functions of two atomic orbitals $\psi_c$ and $\psi_d$ it is done by a simple Linear Combination of Atomic Orbitals (LCAO):
$$\psi = C_1\psi_c + C_2\psi_d$$
Where the $C_1$, $C_2$ are constants.
These two processes of generating the total wave function out of existing ones feel to me as if they should be carried out in the same way. And I understand that we do it this way in the first case since we can’t tell the particles apart and the total wave function therefore has either to be symmetric or antisymmetric. So what is the justification that we can simply add them in the second LCAO case?