My question is on the use of the concept of indistinguishable particles (in quantum mechanics) in a very general context and in particular in statistical mechanics. I have made clear some of my opinions on this subject in a recent paper but when thinking even more about it, I cannot seem to comprehend how come this indistinguishability feature is so universally used and how it is justified.
I am aware of the questions asked here and there and the answers therein but I think they do not exactly address my question.
I am surprised by the specific and intrinsic role that seems to be played by the permutation operator (say $\hat{\Pi}_{12}$ for swapping particles 1 and 2) when compared to other symmetry operators.
In particular, in most cases that come to my mind (rotation, translation, time inversion, parity, etc...), the symmetries of the system are contained in the Lagrangian or Hamiltonian of the system. In simple quantum mechanics for instance, if a system has a certain symmetry, then its hamiltonian commutes with the corresponding representation of this symmetry group. It is then standard that, when we want to have a complete description of the state, we seek a Complete Set of Commuting Observables that will induce a particular structure to the quantum state.
To take the example of two particles considered in many cases, if the hamiltonian commutes with the permutation operator, it is then natural to seek a solution to the Schrodinger equation which is an eigenvector of $\hat{\Pi}_{12}$ with either eigenvalue $+1$ or $-1$.
Now, I perfectly understand that the reasoning above is not enough to get to the symmetrisation postulate (as it has been well answered in the other questions cited above) but my point is that, while the tendency seems usually to go from the Hamiltonian (which I conceptually equate to the modelling of the problem) to the symmetries and the structure of the quantum state, I am surprised that in the case of "identical" particles, the fermion/boson reasoning always precedes the hamiltonian (at the exception of the Standard Model of Particle Physics). In other words it is essentially overlooked (I think) that the apparent indistinguishability of a bunch of particles could be induced by the (maybe simplified or at least effective) modelling of the problem.
My concern is more about general particles such as composite particles like atoms, molecules etc... where it is sometimes (if not often) assumed that they are indistinguishable in the quantum mechanical sense while I would think that they are almost never actually indistinguishable.
If I take a gas of molecules for instance, I seriously wonder what it means to consider them as indistinguishable in the quantum sense since they can have different conformations or electronic states at any finite temperature (spectroscopists know that very well).
My question is then the following:
For an assembly of conceptually identical composite systems (i.e. identical in their composition), is it enough to provide an effective hamiltonian that commutes with the permutation operator to conclude that the so-described particles are indistinguishable in the quantum mechanical sense? If not, what is the rationale to justify the use of fermionic and bosonic statistics for composite systems?
EDIT: I will try to be a bit more specific
Let us consider for now a 2-particle system with the effective hamiltonian
$\hat{H} \equiv \frac{\hat{P}_1^2}{2m} + \frac{\hat{P}_2^2}{2m} + V(\hat{X}_1) + V(\hat{X}_2)$
where the external field $V$ acts in the same on the two particles at a given location.
As far as I am concerned, this hamiltonian commutes with the operator $\hat{\Pi}_{12}$ that swaps labels $1 \leftrightarrow 2$.
Let us now imagine that the particles under study have another observable feature $\hat{F}$ (color, size, spin projection etc...) so that $[\hat{H}, \hat{F}] = 0$ and that the two particles we are looking at may possibly have a different eigenvalue for the feature $\hat{F}$.
Let us now denote $\{|\phi_n \rangle \}_{n=1,2..}$ the eigenstates of the individual hamiltonian for each particle and $\{|f_{\alpha} \rangle \}_{\alpha=1,2..}$ the eigenvectors of the feature $\hat{F}$.
Are these ingredients enough to recquire the complete state to be symmetric or anti-symmetric? e.g.
$|\Psi \rangle = \frac{1}{\sqrt{2}}\left(|\phi_{n_1},f_{\alpha_1}; \phi_{n_2},f_{\alpha_2} \rangle \pm |\phi_{n_2},f_{\alpha_2}; \phi_{n_1},f_{\alpha_1} \rangle \right)$
I would suppose so and Hector's answer seems to indicate that it is the case as well for two particles and that it seems "natural" to extend the chosen symmetric or anti-symmetric property to the N-particle case (to be eventually tested against experiments).