Our goal is to create a population pharmacokinetic (PK) model and identify the best loading dose (LD) of intravenous valproic acid for hospitalized Thai patients. Data from patients who received intravenous valproic acid and underwent measurement of serum valproic acid concentrations during hospitalization were collected retrospectively. A nonlinear mixed-effects modeling approach was conducted to estimate the PK parameters of valproic acid. Covariates affecting the PK parameters of valproic acid were examined and ranked based on their impact on the model's performance. Monte Carlo simulations of 1000 patients were conducted to estimate the optimal LD of valproic acid. A total of 120 hospitalized patients (51.7% male) with 167 valproic acid concentrations were included in the study. A linear one-compartment model with constant residual error was the best base model. An age-covariate model was the best predictor of valproic acid clearance (CL). The typical values of CL and volume of distribution for valproic acid were 0.77 L/h and 14.56 L, respectively. The LD of 1000-1200 mg intravenous was identified as the pragmatic option as an empirical regimen for hospitalized Thai patients. The recommended time to initiate maintenance dose (MD) is 4-8 h following the LD. The population PK model and optimal LD of valproic acid in hospitalized Thai patients has been established, and it may be advisable to initiate the MD at a later time for the elderly.