Coherent states arise naturally in many quantum systems besides the free electromagnetic field. Many of them are mentioned in the following review by WM Zhang, which is actually dedicated to the same question asked here.
I should emphasize that the examples given in this answer refer to the usual Heisenberg-Weyl coherent states as well as generalized coherent states including squeezed states coherent states on Lie groups etc. Please see the following book by Perelomov.
There are many examples of coherent states besides the electromagnetic field; may be the most known examples are:
- Ground states of Bose-Einstein condensates can be approximated by coherent states of the condensate fields.
- The BCS ground state of superconductor can be approximated by a coherent state of the Cooper pairs.
- Ground states of a superfluid can be approximated by a coherent state.
- Ground states of certain phases of ferromagnets and anti- ferromagnets can be approximated by coherent states, please see the following lecture note by Christopher Henley .
The above examples correspond to large systems with spontaneously broken symmetry. This fact stems from the general rule that in spontaneously broken theories, the condensate state can be approximated by a coherent state, which can be heuristically explained as follows, suppose that the number of condensate atoms is $N$. The action of the annihilation operator per unit volume on the condensate state is:
$$\frac{a}{\sqrt{V}}|N\rangle = \sqrt{\frac{N}{V}}|N-1\rangle$$
When $N\rightarrow\infty$, the states on both sides of the equation have the same limit and we obtain:
$$ \lim_{N, V, \rightarrow \infty} \frac{a}{\sqrt{V}}|\mathrm{Condensate}\rangle = \sqrt{\rho}|\mathrm{Condensate}\rangle$$
Where $\rho$ is the condensate density. This shows that the condensate state is a coherent state.
In systems of interacting of atoms with radiation fields, the atoms can be approximated by two state systems. An atomic gas in this approximation, has an SU(2) (angular momentum) symmetry of the two state systems. In this case also when this symmetry is broken, the atomic gas will reside in a coherent state of the angular momentum. This is a generalized coherent state of Perelomov- Gilmore type and still the above rule applies.
Coherent states are the solution of the Hartree-Fock approximations in atomic and molecular physics. The slater determinant is actually a coherent state over a Grassmann manifold. Please see the following review by Chiumiento, and Melgaard.
In quantum field theories with massless particles, it was found that using coherent state bases of asymptotic states instead of the particle number basis leads to a faster cancelation of the infrared divergences, please, see for example the following work by More and Misra and the following thesis by Catana. This is because the asymptotic charged fermion states are surrounded by a coherent state cloud of soft photons.
It is also worthwhile also to mention that electrons restricted to the lowest Landau level are described by coherent states. Please see the following lecture note by Christopher Henley.
Although not requested in the question, I have to mention that coherent states have many applications in mathematics. Their importance lies in the fact that they generate coherent state representations of groups, please see for example the following work by Berceanu and Schlichenmaier. Their importance lies in the fact that their construction philosophy can be extended to other mathematical objects such as quantum groups, groupoids which have not found extensive application in physics yet.
Coherent states (in systems other than the free electromagnetic field) can also be manufactured, rather than naturally obtained. This is actually a problem in quantum control and one of the principles that can be applied is the following:
In contrast orthonormal bases in a Hilbert space, which are parametrized by a discrete variable, coherent states are parametrized by a phase space. For example a harmonic oscillator coherent state $|z\rangle$ is parametrized by $z\in \mathbb{C}$ the phase space of a single degree of freedom particle. There is a family of Hamiltonians whose action on the coherent states can be implemented as a classical action on the phase space variable:
$$ e^{i\hat{H}t} |z\rangle = |\phi_t \cdot z\rangle$$
Where $\phi_t \cdot z$ is the Hamiltonian $H$ flow of the point $z$, namely the solution of the Hamilton equation of motion. The relation between the quantum Hamiltonian $\hat{H}$ and its classical symbol $H$ is given be certain quantization rules, in our case they can be chosen to be the Wick or anti-Wick ordering.
In the case of the harmonic oscillator coherent states, the distinguished family of Hamiltonians are the linear Hamiltonians:
$$\hat{H} = \lambda(t) a + \bar{\lambda(t)} a^{\dagger}$$
Thus these Hamiltonians make the system evolve within the coherent state manifold. In order to bring the system to a desired coherent state, we must start from an easily prepared coherent state (in the Harmonic oscillator the ground state) and act on it by this type of Hamiltonian. The system will evolve to a general coherent state controlled by our choice of the Hamiltonian.
For example, the Harmonic oscillator whose total Hamiltonian
$$ H_T = \omega a^{\dagger} a + \lambda(t) a + \bar{\lambda(t)} a^{\dagger}$$
evolves to a coherent state $|\Lambda(t)\rangle$, where
$$\Lambda(t) = \int_{0}^{t} e^{i \omega\tau} \lambda(\tau) d\tau$$
This method is used for example to generate phonon coherent states.
There is another principle allowing the manufacture of coherent states: There are cases where coherent states appear as steady states of open system quantum dynamics, approximated by Lindblad equation. For example those implemented in cooling processes in which Bose-Einstein condensates are prepared. See for example the following work by Albert, Bradlyn, Fraas, and Jiang