Measurements of the pd → pspdK + K − reaction, where psp is a spectator proton, have been undertaken at the Cooler Synchrotron COSY-Jülich by detecting a fast deuteron in coincidence with a K + K − pair in the ANKE facility. Although the proton beam energy was fixed, the moving target neutron allowed values of the non-resonant quasi-free pn → dK + K − total cross section to be deduced up to an excess energy ǫ ≈ 100 MeV. Evidence is found for the effects of K − d and KK final state interactions. The comparison of these data with those of pp → ppKshows that all the total cross sections are very similar in magnitude. We have recently published measurements of the differential and total cross sections for the pp → ppK + K − reaction at three energies close to threshold [1]. A major challenge in the analysis was the separation of the contribution from the production and decay of the φ meson from that of the non-φ component [2]. One of the striking features of the non-φ results is the strong attraction between the K − and each of the final protons seen in the differential distributions. This also has a major effect on the energy dependence of the total cross section, enhancing it at low energies. Although tantalizing, these results do not, however, resolve the ongoing question as to whether the interaction is sufficiently strong to allow the K − to form a bound state with the two protons [3,4,5].The isospin dependence of φ production has been studied through an investigation of pd → p sp dK. By identifying the final deuteron and kaon pair and measuring their momenta, it was possible to construct the momentum of the recoil proton p sp to show that it was consistent with being a spectator, whose only significant participation in a reaction is through a change in the kinematics. Interpreting the results in this way, it was possible to extract values of the quasi-free pn → dK + K − cross section. Moreover, although the experiment was carried out at one fixed beam energy, the movement of the target neutron enabled data to be obtained over a wide range of excess energy ǫ = √ s − m d − 2m K on an event-by-event basis. Just as in the pp → ppK + K −