Hot-mix asphalt is a mixture of aggregates and asphalt binder in appropriate ratios to produce a high-performing material for asphalt pavements. The aggregate structure, which depends on the gradation, is an important factor in determining the volumetric properties of HMA. The design process to determine the optimal aggregate blend is currently iterative and engineers rely almost exclusively on experience. This approach is time consuming and often results in suboptimal HMA mixtures. This study presents linear programming optimization models and attendant solution procedures that minimize HMA aggregate cost and produce high-quality HMA. The models are validated with real-life examples, and results indicate that the models are useful to replicate HMA mixes during field modifications, reduce the aggregate cost in a mixture, and manage stockpile inventory. The application of optimization models will increase the application of the Bailey method in the United States. © 2011 American Society of Civil Engineers.