We aim to describe the population pharmacokinetics (PPK) of tacrolimus in Chinese pediatric patients under 4 years old after liver transplantation and to develop individualized tacrolimus dosing software.A total of 663 blood concentrations from 85 patients aged 4.57 months to 3.97 years were collected in this study. PPK analysis was performed using a nonlinear mixed effects modeling approach with the software, Phoenix. Using C#, an individualized tacrolimus dosing software was created. The software was then used to predict the concentrations of another ten pediatric liver transplantation patients to verify the accuracy of said software. The predictive error (PE) and the absolute predictive error (APE) for each predicted time point were computed.A one-compartment model with first-order elimination best fitted the data. The apparent volume of distribution (V/F) and apparent clearance (CL/F) were 198.65 L and 2.41 L/h. Postoperative days (POD), total bilirubin (TBIL), and the use of voriconazole significantly influenced tacrolimus apparent clearance. The incorporation of an increasing number of actual blood drug concentrations into the prediction resulted in a decrease in both PE (72%, 17%, 7%) and APE (87%, 53%, 26%).A qualified PPK model of tacrolimus was developed in Chinese pediatric patients. The individualized tacrolimus dosing software could be used as a suitable tool for the personalization of tacrolimus dosing for pediatric patients after liver transplantation.© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.