Collective two-color photoassociation of a freely-interacting 87 Rb Bose-Einstein condensate is theoretically examined, focusing on stimulated Raman adiabatic passage (STIRAP) from an atomic to a stable molecular condensate. In particular, Drummond et al. [Phys. Rev. A 65, 063619 (2002)] have predicted that particle-particle interactions can limit the efficiency of collective atom-molecule STIRAP, and that optimizing the laser parameters can partially overcome this limitation. We suggest that the molecular conversion efficiency can be further improved by treating the initial condensate density as an optimization parameter.PACS numbers: 03.75. Nt,03.65.Ge,05.30.Jp,32.80.Wr Stimulated Raman adiabatic passage (STIRAP) from an atomic to a molecular condensate [1] in the presence of particle-particle interactions was recently investigated [2,3]. In particular, Drummond et al.[3] predict that, for a 87 Rb Bose-Einstein condensate (BEC) of typical density (ρ ∼ 10 14 cm 3 ), the STIRAP conversion efficiency can be limited in practice by two-photon dephasing caused by particle interactions, and that said limitation can be partially overcome by optimizing the laser parameters. The purpose of this Comment is to suggest that the role of particle-particle interactions can be further downplayed, and the conversion efficiency improved, by treating the initial condensate density as an additional optimization parameter.The mean-field equations for collective two-color photoassociation of a freely-interacting gas can be writtenwhere the complex amplitudes a, b, and g represent the respective atomic, excited-molecular, and stablemolecular condensates. The laser-matter interactions that drive the respective atom-molecule and moleculemolecule transitions areand, where χ 0 includes the effects of Bose-enhancement [1], i.e., χ 0 ∝ √ ρ . The two-photon (intermediate) detuning is ∆ (δ), and the mean-field shift due to atom-atom (atommolecule, molecule-molecule) interactions is Λ aa = ρλ aa = 4π ρa aa /m (Λ ag = ρλ gg = 3π ρa ag /m, Λ gg = ρλ gg = 2π ρa gg /m), where m is the atom mass and a aa (a ag , a gg ) is the atom-atom (atom-molecule, moleculemolecule) scattering length. Spontaneous decay is included with the rate γ s , which is generally large enough to justify neglect of any mean-field shifts for the excitedmolecular state.Explicit numbers for 87 Rb are [3] γ s = 7.4 × 10 7 s −1 , χ 0 = 2.1 × 10 6 ρ/ρ 0 s −1 , ρ 0 = 4.3 × 10 14 cm −3 , λ aa = 4.96 × 10 −11 cm 3 /s, λ ag = −6.44 × 10 −11 cm 3 /s; although unknown, the stable-molecule mean-field shift λ gg = 2.48 × 10 −11 cm 3 /s is estimated by assuming equal atom-atom and molecule-molecule scattering lengths.The idea of optimizing the density of the initial atomic Bose-Einstein condensate is drawn from investigations into forming a molecular condensate via Feshbachresonant interactions[4], i.e., magnetoassociation, where collision-induced vibrational relaxation of the molecules is included as a complex particle-particle scattering length, and where a moderate density helps to...