d_{x^2-y^2} pairing of composite excitations in the 2D Hubbard model
  We report on a strong coupling approach (on-site Coulomb repulsion, U larger than the nearest-neighbour hopping energy |t|) to the Hubbard model. Starting from the Hubbard operators which diagonalize the interaction term, we generate a hierarchy of composite operators from the equations of motion. Using the Hubbard operators as a basis, we are able to compute the associated Green functions including the anomalous Green functions which describe pair formation. We show explicitly that these anomalous Green functions are non-zero in the d_{x^2-y^2} channel; however, the entities that pair up are not single electron-like particles but rather composite excitations (which we call cexons) made out of an electron and a hole on nearest-neighbour sites. Cexons are fermionic in nature as they have spin 1/2 and also have unit charge. Our calculations of the chemical potential reveal that negative compressibility in the 2D Hubbard model and composite excitation pairing are intimately connected, namely, the larger the negative compressibility, the larger the pairing amplitude. Our observation of negative compressibility in the under-doped regime is consistent with phase segregation or stripe formation in the normal state. While pairing ameliorates the negative compresssibility, it does not eliminate it entirely. In addition, we find that the anomalous correlation functions are particle-hole symmetric and exhibit a maximum at a doping level of roughly 10% as measured from half-filling. For U=8|t|, the onset temperature for pair formation is 0.02|t|. 
Tudor D. Stanescu, Ivar Martin, and Philip Phillips