The self-reaction of propargyl (C3H3) radicals has been widely suggested as one of the key routes forming benzene in a variety of aliphatic flames. Currently, in the majority of aromatic models, the C3H3 + C3H3 submechanism often contains one or two C6H6 isomers and a few global reaction steps, which do not adequately represent the actual recombination chemistry. Recent experimental and theoretical studies on the direct propargyl recombination and subsequent C6H6 isomerization have provided sufficient information to revisit and revise the C3H3 + C3H3 reaction submechanism. In the present work, a semidetailed kinetic model consisting of seven isomeric C6H6 species and 14 reaction steps was constructed based on the most recent potential energy surface for this system. The trial model was subjected to systemic optimization by use of a recently developed physically bounded Gauss-Newton (PGN) method against detailed species profiles of direct propargyl recombination and 1,5-hexadiyne (15HD) isomerization obtained front experiments at high temperatures in a shock tube and at low temperatures in a flow reactor, which were all measured at very high pressure (shock tube) or atmospheric (flow reactor) conditions. Predictions of the optimized model were in excellent agreement with all experimental measurements. The optimized C3H3 + C3H3 reaction subset was also tested for flame modeling. Two different aromatic chemistry models that incorporate benzene formation from propargyl radicals as a single step reaction were modified to include the complete submechanism for propargyl recombination. The updated models predict significant percentages of three isomeric species [2-ethynyl-1,3-butadiene (2E13BD), fulvene, and benzene] in premixed fuel-rich acetylene and ethylene flames, reflecting the observed flame structures.