How can I derive fusion rules for anyons?

Naively, the collection of data which describes a model of anyons is not its fusion rules, but rather its modular data - that is, a pair of matrices $S$ and $T$ which generate a representation of the (modular) group $SL(2,\mathbb Z)$.

The (diagonal) matrix $T$ encodes the mutual statistics of quasi-particles - that is what happens when you exchange two identical excitations. The entries of this matrix are called twists. Vafa's Theorem says that all of these must be roots of unity. An excitation is anyonic if and only if it has twist not equal to $1$ or $-1$. The $S$ matrix encodes evaluations of colored Hopf links when braiding pairs of different excitations, but that's less important at the moment. It is however important that the $S$ matrix be non-degenerate.

A very important thing is that the Verlinde formula can be used to derive the fusion rules from the $(S,T)$ matrices. The converse is not true - there are plenty of examples (Ising, Fibonacci) where a set of fusion rules may have give rise to more than one set of $S$ and $T$ matrices. Taking everything together through the magic of polynomials, these give rise to mathematical objects called modular tensor categories, and it is really here that the story begins. MTCs are the objects which describe the algebraic properties for a given model of anyons.

Important note: I take issue with the phrase "you can just make them up" because it makes a set of very big problems seem trivial. It also takes as fact something which is only at this point conjecture - that every pair $(S,T)$ is uniquely realized by some MTC.

We can now give a partial answer the parenthetical to the first question: If a set of fusion rules does not give rise to a set of modular data, there does not exist a model (though we've still not said what a model is) with anyonic excitations realizing those fusion rules.

Another note: This statement sweeps several things under the rug like $G$-crossed modular categories and minimal modularizations of ribbon categories. Let's not complicate matters.

The next important question then is: What is a model with anyonic excitations? Beyond their exchange properties we also know some other physical things our model with anyons should do. The big ones are that there is degeneracy in the ground state and that our anyonic excitations should be robust against local perturbations. Add all of the bits and pieces up and you get that the lowest energy Hilbert space should be described by a topological quantum field theory or TQFT.

Usually, the information for a physical model is contained in some Hamiltonian $H$ or Lagrangian $L$, and here the situation will be no different. To say that I have a a model with anyonic excitations will be to say that I have some $H$ or $L$ whose lowest energy Hilbert space is described by a TQFT. A choice of TQFT uniquely determines the MTC (i.e. the anyonic excitations), and given an MTC and a manifold there is a way to construct a TQFT.

Now then, the questions can be rephrased as follows:

1. Given a modular tensor category $\mathcal C$, how does one construct a Hamiltonian $H$ or Lagrangian $L$ realizing a TQFT whose quasi-particle excitations are described by $\mathcal C$?
2. Given a Hamiltonian $H$ or Lagrangian $L$ realizing a TQFT, how does one determine the MTC $\mathcal C$ describing its anyonic excitations?

As far as I know, the answer to both questions is incomplete. I do not know a general recipe for producing a Hamiltonian/Lagrangian which realizes an arbitrary MTC. Conversely, my understanding of the opposite direction is that given $H$/$L$, one must first show that the ground state is degenerate and then go through and construct the appropriate operators describing the anyonic excitations. I don't know much about that, but it seems as if it could be a tricky proposition.

Because of this, what I've included below are classes of examples in which both questions have known answers.

1. For MTCs which can be realized as centers of some unitary fusion categories, one can construct a Levin-Wen Hamiltonian on a trivalent lattice. In fact, Levin-Wen models take as input the UFC and lattice, and there is a formula to follow for constructing the quasi-particle operators which correspond to objects in the center. Levin's thesis has of the details of this.
2. A particularly nice subclass of these are Kitaev's double models. What one does in that case is play the game.
3. To provide a partial example which starts in the other direction, one could start with a manifold $M$, Lie group $G$ and "level" (essentially an integer) $k$ and construct its Chern-Simons action. After some work - which I won't try to describe because I'm not terribly familiar with it - one obtains a TQFT whose quasiparticle excitations are described by the MTC for the category of representations for the affine Lie algebra associated to $G$ with highest integral weight $k$. Both Ising and Fibonacci are examples of this. Ising is the category $SU(2)_4$ and Fibonacci is the even half of the category $SU(2)_3$.