The Modular Neutron Array (MoNA) and 4 Tm Sweeper magnet were used to measure the free neutrons and heavy charged particles from the radioactive ion beam induced 32 Mg + 9 Be reaction. The fragmentation reaction was simulated with the constrained molecular dynamics model (CoMD), which demonstrated that the N/Z of the heavy fragments and free neutron multiplicities were observables sensitive to the density dependence of the symmetry energy at subsaturation densities. Through comparison of these simulations with the experimental data, constraints on the density dependence of the symmetry energy were extracted. The advantage of radioactive ion beams as a probe of the symmetry energy is demonstrated through examination of CoMD calculations for stable and radioactive-beam-induced reactions. PACS number(s): 21.65. Mn, 25.70.Mn, 21.65.Ef Introduction. The desire to extend our understanding of nuclear matter at densities, temperatures, pressures, and neutron-to-proton ratios (N/Z) away from that of ground-state nuclei has become a driving force in the nuclear science community. In particular, the emergence of radioactive ion beam (RIB) facilities has placed an emphasis on exploring nuclear matter along the isospin degree of freedom. The symmetry energy is the critical component which defines how the properties of nuclear matter, or the nuclear equation of state (EoS), change as a function of isospin. The nuclear EoS can be approximated aswhere the energy per nucleon of infinite nuclear matter, E(ρ, δ), is a function of the density (ρ) and isospin concentration (δ) [1][2][3]. The isospin concentration is the difference in the neutron and proton densities, δ = (ρ n − ρ p )/(ρ n + ρ p ) ≈ (N − Z)/A. The first term of the EoS is isopin independent and thus represents the binding energy of symmetric (N = Z) nuclear matter. The second term of the EoS has a strong dependence on the asymmetry of the nuclear matter and has been historically termed the symmetry energy (however, a more appropriate term is the asymmetry energy). The EoS for symmetric nuclear matter is relatively well constrained around the saturation density (ρ 0 = 0.16 fm −3 ) from isoscalar giant monopole and dipole resonances [4,5] and at higher densities (up to ρ/ρ 0 ∼4.5) from heavy-ion collisions [6,7]. Constraining the form of the density dependence of the symmetry energy [E sym (ρ)] is essential for developing a complete description of asymmetric nuclear matter. For example, the properties of neutron stars are strongly correlated to the * kohley@nscl.msu.edu † Present address: TRIUMF, 4004 Wesbrook Mall, Vancouver, British Columbia V6T 2A3, Canada.