Studying antibody dynamics following re-exposure to infection and/or vaccination is crucial for a better understanding of fundamental immunological processes, vaccine development, and health policy research. We adopted a nonlinear mixed modeling approach based on ordinary differential equations (ODE) to characterize varicella-zoster virus specific antibody dynamics during and after clinical herpes zoster. Our ODEs models convert underlying immunological processes into mathematical formulations, allowing for testable data analysis. In order to cope with inter- and intra-individual variability, mixed models include population-averaged parameters (fixed effects) and individual-specific parameters (random effects). We explored the use of various ODE-based nonlinear mixed models to describe longitudinally collected markers of immunological response in 61 herpes zoster patients. Starting from a general formulation of such models, we study different plausible processes underlying observed antibody titer concentrations over time, including various individual-specific parameters. Among the converged models, the best fitting and most parsimonious model implies that once Varicella-zoster virus (VZV) reactivation is clinically apparent (i.e., Herpes-zoster (HZ) can be diagnosed), short-living and long-living antibody secreting cells (SASC and LASC, respectively) will not expand anymore. Additionally, we investigated the relationship between age and viral load on SASC using a covariate model to gain a deeper understanding of the population's characteristics. The results of this study provide crucial and unique insights that can aid in improving our understanding of VZV antibody dynamics and in making more accurate projections regarding the potential impact of vaccines.