Trajectory surface hopping calculations have been carried out for collisions of Ar ϩ ϩH 2 and Ar ϩ ϩHD on three low-lying potential energy surfaces projected from the original six in the Kuntz and Roach diatomics in molecules surface for this system. The location and probability of hops between surfaces were determined using the new algorithm developed by Parlant and Gislason. In addition to the reactive channel and total charge transfer to H 2 ϩ and HD ϩ , dissociative channels to, for example, Ar ϩ ϩHϩH, and ArϩH ϩ ϩH have been studied. Particular attention was paid to the dissociative charge transfer isotope effect for the processes Ar ϩ ϩHD→ArϩH ϩ ϩD, or ArϩHϩD ϩ ; near threshold the D ϩ product is favored over H ϩ which we attribute to preferential dissociation of excited ArD ϩ products. This is the first theoretical study of these dissociation processes.