I'm reading the chapter about the renormalization group in Yeoman's book "Statistical mechanics of phase transitions" and I'm puzzled about how the author relates the scaling of the RG with the critical exponents. We have some RG map on the Hamiltonian $H\rightarrow R(H)$. We suppose that we are close to the fixed point $H^* $, so
$$H'=R(H^*+\delta H)=H^* + A(H^*)\delta H$$
where $A$ is a matrix and $\delta H$ is seen as a vector with the coupling constants as components. This matrix can be diagonalized and we can write
$$ A(H^*)\delta H= A(H^*)\sum_k\mu_k \Phi_k=\sum_k\lambda_k\mu_k \Phi_k\tag {$\star$}$$
where $\Phi_k$ are functions of the lattice and $\lambda_k$ are the eigenvalues of $A$. It's easy to argue that they must have the form
$$ \lambda_k=b^{y_k}$$
where $b$ is the scaling factor of the map. No problem until here. If $y_k>0$ we call it relevant, otherwise irrelevant.
Then the author says that for the Ising model the relevant couplings are the temperature and the magnetic field, and here I'm already confused. First, because $A$ depends on what renormalization scheme we choose, so how can we predict what the eigenvectors are without saying more about $R$? Second, the Ising Hamiltonian (if we absorb temperature in it) is
$$ H_I=-\beta J \sum_{\langle i,j\rangle}s_is_j-\beta h\sum_i s_i$$
I don't see how $\beta$ and $h$ could ever appear as in $(\star)$ as two linearly superposed terms. We would need something like
$$ H_I=H^*+\beta\Phi_\beta+h\Phi_h$$
but I don't see how this can be true because $\beta h$ appears in the Hamiltonian, it seems to me that $\beta h$ must be treated as a single coupling, and you can't take $\beta$ and $h$ separately. In short, my first question is
How can one treat temperature as a coupling constant if it appears in all the coupling constants?
There were originally two parts of this question, following advice of Adam I posted the second part in a separate question