IntroductionLong-term simulation of tooth movement is crucial for clear aligner (CA) treatment. This study aimed to investigate the effect of maxillary molar distalization with CA via an automatic staging simulation.MethodA finite-element method (FEM) model of maxillary dentition, periodontal ligaments, attachments, and corresponding CA was established, and a prescribed 2-mm distalization with 0.1 mm each step of the second molar was simulated. The long-term tooth movement under orthodontic force was simulated with an iterative computation method. The morphologic changes of CA during staging were simulated with the thermal expansion method.ResultsTwenty steps of molar distalization were simulated. Significant distal tilting of the second molar was revealed, along with the proclination of anterior teeth, which caused the ‘reversed bow effect’. For the second molar, 4.63°distal tilting at the 20th step was revealed. The intrusion of the incisors and the second molar were 0.43 mm, 0.39 mm, and 0.45 mm, respectively, at step 20. All the anterior teeth showed a proclination of approximately 1.41°–2.01° at the 20th step. The expression rate of the designed distalization of the second molar was relatively low (approximately 68%) compared to the high efficacy of interdental space opening between molars with CA (approximately 89%).ConclusionA novel method of simulating long-term molar distalization with CA with FEM was developed. The FEM results suggested distal tilting of the second molar and the proclination of anterior teeth during the molar distalization.