What is the meaning of commuting Hamiltonians? 
I have two quantum mechanical Hamiltonians such that \begin{equation}
 [\hat{H}_1,\hat{H}_2] = 0, \end{equation} where $\hat{H}_1$ and
  $\hat{H}_2$ act on the same set of states. What is there to infer
  physically about these two Hamiltonians? Are there further
  mathematical subtleties that were not brought out in "What is the
  Physical Meaning of Commutation of Two Operators?" for the case of
  two Hamiltonians?

From looking around I have these properties:


*

*Treating them as observables I can measure them simultaneously.

*They share the same set of eigenstates and thus any state can be expanded as a sum of these.

*They can both be simultaneously diagonalised.


Are there further properties or subtleties to this relationship?
EDIT : Edited after comments pointing out that there is little to be said if they individually act on different subsystems, apart from the fact that they share eigenstates when acting on both systems together.
 A: Comments to the question (v2): 


*

*It seems that the question does not explain how a 'Hamiltonian' $H$ differs from a self-adjoint operator $A$ (presumably bounded from below). This would make OP's question a duplicate of the linked Phys.SE post. 

*Perhaps a 'Hamiltonian' $H$ is also supposed to generate 'time'-evolution for some distinguished parameter $t$, which may or may not be actual time? Then consider a universe with two 'time' directions $t_1$ and $t_2$, cf. e.g. this Phys.SE post. The two commuting Hamiltonians $[H_1,H_2]=0$ means that one gets the same result if one first 'time'-evolve $e^{iH_1t_1/\hbar}$ wrt. $t_1$ and then 'time'-evolve $e^{iH_2t_2/\hbar}$ wrt. $t_2$, as one would get if one does it the other way around. In other words, $t_1$ and $t_2$ represent commuting flows, and it makes sense to specify a state with two 'time'-coordinates $(t_1,t_2)$.
A: At least a partial answer to your question is that commuting hamiltonians help you to solve the physical system described by one of them: in particular, if your system has $N$ degrees of freedom and you have $N$ commuting hamiltonians, there is good hope that you can trivialize the problem and solve it exactly. In classical mechanics, this is known as Liouville integrability (where there the commutativity is associated to a Poisson bracket). In quantum mechanics, the notion is not completely well-defined, although searching for the term quantum integrability will provide you with ample reading material. Due to the equation of motion in the Heisenberg picture
$$
\frac{dA}{dt} = [A,H],
$$
where $A$ is an operator, we see that if $H_1$ is the hamiltonian governing time evolution, $H_2$ is conserved in time. But since many known Hilbert spaces are infinite dimensional, using the above in practice is harder. 
Note that for quantum mechanical models such as spin chains ($N$ fixed particles interacting via spin degrees of freedom), the space of states is finite dimensional and all of the above does help. 
