Clustering in low density nuclear matter has been investigated using the NIMROD multidetector at Texas A&M University. Thermal coalescence modes were employed to extract densities, ρ, and temperatures, T, for evolving systems formed in collisions of 47A MeV (40)Ar+(112)Sn, (124)Sn and (64)Zn+(112)Sn, (124)Sn. The yields of d, t, (3)He, and (4)He have been determined at ρ=0.002 to 0.03 nucleons/fm(3) and T=5 to 11 MeV. The experimentally derived equilibrium constants for α particle production are compared with those predicted by a number of astrophysical equations of state. The data provide important new constraints on the model calculations.