Simulated historical extreme precipitation is evaluated for Coupled Model Intercomparison Project Phase 6 (CMIP6) models using precipitation indices defined by the Expert Team on Climate Change Detection and Indices (ETCCDI). The indices of 33 Global Circulation Models (GCMs) are evaluated against corresponding indices with observations from the Global Climate Center Precipitation Dataset (GPCC V2020) over five sub-regions across Arid Central Asia (ACA), using the Taylor diagram, interannual variability skill score (IVS) and comprehensive rating index (MR). Moreover, we compare four multi-model ensemble approaches: arithmetic average multi-model ensemble (AMME), median multi-model ensemble (MME), pattern performance-based multi-model ensemble (MM-PERF) and independence weighted mean (IWM). The results show that CMIP6 models have a certain ability to simulate the spatial distribution of extreme precipitation in ACA and the best ability to simulate simple daily intensity (SDII), but it is difficult to capture the spatial bias of consecutive wet days (CWD). Almost all models represent different degrees of wet bias in the southern Xinjiang (SX). Most GCMs are generally able to capture extreme precipitation trends, but to reproduce the performance of interannual variability for heavy precipitation days (R10mm), SDII and CWD need to be improved. The four multi-model ensemble methods can reduce the internal system bias and variability within individual models and outperform individual models in capturing the spatial and temporal variability of extreme precipitation. However, significant uncertainties remain in the simulation of extreme precipitation indices in SX and Tianshan Mountain (TM). Comparatively, IWM simulations of extreme precipitation in the ACA and its sub-regions are more reliable. The results of this study can provide a reference for the application of GCMs in ACA and sub-regions and can also reduce the uncertainty and increase the reliability of future climate change projections through the optimal multi-model ensemble method.