BACKGROUND. CT is increasingly detecting thyroid nodules. Prior studies indicated a potential role of CT-based radiomics models in characterizing thyroid nodules, although these studies lacked external validation. OBJECTIVE. The purpose of this study was to develop and validate a CT-based radiomics model for the differentiation of benign and malignant thyroid nodules. METHODS. This retrospective study included 378 patients (mean age, 46.3 ± 13.9 [SD] years; 86 men, 292 women) with 408 resected thyroid nodules (145 benign, 263 malignant) from two centers (center 1: 293 nodules, January 2018 to December 2022; center 2: 115 nodules, January 2020 to December 2022) who underwent preoperative multiphase neck CT (noncontrast, arterial, and venous phases). Nodules from center 1 were divided into training (n = 206) and internal validation (n = 87) sets; all nodules from center 2 formed an external validation set. Radiologists assessed nodules for morphologic CT features. Nodules were manually segmented on all phases, and radiomic features were extracted. Conventional (clinical and morphologic CT), noncontrast CT radiomics, arterial phase CT radiomics, venous phase CT radiomics, multiphase CT radiomics, and combined (clinical, morphologic CT, and multiphase CT radiomics) models were established using feature selection methods and evaluated by ROC curve analysis, calibration-curve analysis, and decision-curve analysis. RESULTS. The combined model included patient age, three morphologic features (cystic change, "edge interruption" sign, abnormal cervical lymph nodes), and 28 radiomic features (from all three phases). In the external validation set, the combined model had an AUC of 0.923, and, at an optimal threshold derived in the training set, sensitivity of 84.0%, specificity of 94.1%, and accuracy of 87.0%. In the external validation set, the AUC was significantly higher for the combined model than for the conventional model (0.827), noncontrast CT radiomics model (0.847), arterial phase CT radiomics model (0.826), venous phase CT radiomics model (0.773), and multiphase CT radiomics model (0.824) (all p < .05). In the external validation set, the calibration curves indicated the lowest (i.e., best) Brier score for the combined model; in the decision-curve analysis, the combined model had the highest net benefit for most of the range of threshold probabilities. CONCLUSION. A combined model incorporating clinical, morphologic CT, and multiphase CT radiomics features exhibited robust performance in differentiating benign and malignant thyroid nodules. CLINICAL IMPACT. The combined radiomics model may help guide further management for thyroid nodules detected on CT.