BCS wave function in Neutron stars I've heard mentioned in various classes that neutron stars, like superconductors, are described by BCS theory. I know that in superconductors a key element in forming cooper pairs is a net attractive force between the electrons which would normally repel one another. That attractive force is accounted for via lattice vibrations (phonons) created and "absorbed" by electrons. 
So my question is: what provides the attractive force between neutrons? just gravity? If it is true that neutron stars follow BCS theory, by what means was someone able to verify that?
 A: I'm afraid the actual situation is much more complicated than you've been told. For one thing, the superconductivity does not occur between neutrons, but between quarks themselves. The topic of high density QCD is a very cool interplay of condensed matter and high energy physics, and a very nice review is available by Frank Wilczek. However, that article does need some background in QCD and superconductivity simultaneously to appreciate.
A shortened version might go something like this:


*

*Inspiration: free fermions are incredibly unstable to superconductivity, in that any attractive interaction will cause it (in fact, there's an old theorem by Pierls (?) that almost all interactions (even repulsive ones) will cause superconductivity if you cool far enough). In QCD, quarks naturally attract already! So at sufficiently high density and low temperature, we can imagine that QCD can cause a strong attractive instability to a Fermi gas of quarks.

*Complications: 6 flavours, chiralities, masses of quarks are different, etc.

*Simplification by complication: realise that the normal state of QCD (i.e. 3- and 2-quark combinations) is just that: one possible state. Other phases of QCD exist, and we can study the phase boundaries and so forth even if we can't compute things exactly (universality saves the day!). We find that at really high densities, quarks pair up to give a background diquark condensate, through which single quarks move, and through the Anderson-Higgs mechanism gains a large mass by eating some Goldstone modes. All gluons become gapped (again, Anderson-Higgs), apart from one which mixes with the photon. The symmetries of this solution is actually the same as that of normal matter --- replace baryons with quarks (+ their diquark condensed background) and mesons with diquarks; this suggests that they are really the same phase in theory. In practise, getting from one to another requires some other phases in the middle, which are more complicated and arise due to the quark masses and the number of flavours, etc.

