Abrupt acceleration of the syringe of an autoinjector (AI) upon rod-plunger impact may induce undesired severe cavitation events and impose extraneous stresses upon the device, leading to device failure. Cavitation results from a rapid and significant pressure drop in a liquid, leading to the formation and growth of small vapor-filled cavities. Upon collapse, these cavities generate an intense shock wave that may lead to protein aggregation and device container damage and shatter. Since the maximum acceleration of the syringe depends upon the operating conditions of the AI, the severity of cavitation will likewise depend on the operating conditions of the AI. Likewise, injection time and ensuring proper needle displacement before drug release also depend on operating conditions, making optimization of the autoinjector a multiobjective optimization problem.Therefore, in this study, optimization of an autoinjector to limit cavitation severity is pursued via an experimentally validated computational model for cavitation in spring-driven autoinjectors. Our goal is to locate AI design configurations that balance maximizing device performance and patient comfort and minimizing the risks of device damage and severe cavitation upon actuation.Relevant parameters of interest are the drive spring force, air gap height, solution viscosity, friction between the rod and spring, frictional force on the plunger, rates of change of frictional force on the plunger, elasticity of plunger, viscosity of the plunger, and initial displacement between the plunger and the driving rod. The kinematics of the syringe barrel, needle displacement (travel distance) at the start of drug delivery, and injection time are gathered using an experimentally validated autoinjector kinematics model. At the same time, cavitation bubble dynamics are resolved using an experimentally validated cavitation model that takes the temporal displacement of the syringe and temporal air gap pressure as inputs.We use our experimentally validated models to explore the parameter space and understand the driving factors of our desired outcomes. Subsequently, we pose the design problem as a multi-objective optimization problem and develop a deep neural network surrogate model supplemented with iterative learning to speed up optimization. A variance-based sensitivity analysis was performed to determine the sensitivity and influence of design parameters on the outcomes, and the main contributors to the outcomes of interest were isolated. Using a multi-objective optimization framework, we located 300 + successful candidates and evaluated them through uncertainty analysis to identify three promising candidates that meet all criteria for drug viscosities of interest. Finally, we show that this methodology can be used to conduct hypothesis testing, leading to novel design configurations.