This work develops a novel superstructure considering a reaction−separation system for aromatics production from methanol, which consists of methanol aromatization, water removal, product separation, and light hydrocarbon aromatization units. On the basis of rigorous simulation of these units in Aspen HYSYS, a multiobjective simulation−optimization model is established for synthesis of methanol to aromatics process to simultaneously optimize the net present value and unit ReCiPe score. A matrix-based method is proposed to determine the Pareto-optimal curve and identify the optimal production routes through data interaction between MATLAB and Aspen HYSYS. The most profitable production route generates a net present value 39.0% higher than the lowest one, and the most environmentally sustainable production route shows a unit ReCiPe score 66.42% lower than the worst one. Monte Carlo simulation is used to investigate the effects of uncertain economic parameters on the three obtained Pareto-optimal production routes, which present robustness in the face of fluctuations of raw material and product prices as well as capital investment. Besides, the optimal production routes for variant reaction effluents are also analyzed to show the insight in production route selection.