This research focuses on the predictive modeling between rocks' dynamic properties and the optimization of neural network models. For this purpose, the rocks' dynamic properties were measured in terms of quality factor (Q), resonance frequency (FR), acoustic impedance (Z), oscillation decay factor (α), and dynamic Poisson’s ratio (v). Rock samples were tested in both longitudinal and torsion modes. Their ratios were taken to reduce data variability and make them dimensionless for analysis. Results showed that with the increase in excitation frequencies, the stiffness of the rocks got increased because of the plastic deformation of pre-existing cracks and then started to decrease due to the development of new microcracks. After the evaluation of the rocks’ dynamic behavior, thevwas estimated by the prediction modeling. Overall, 15 models were developed by using the backpropagation neural network algorithms including feed-forward, cascade-forward, and Elman. Among all models, the feed-forward model with 40 neurons was considered as best one due to its comparatively good performance in the learning and validation phases. The value of the coefficient of determination (R2= 0.797) for the feed-forward model was found higher than the rest of the models. To further improve its quality, the model was optimized using the meta-heuristic algorithm (i.e. particle swarm optimizer). The optimizer ameliorated its R2values from 0.797 to 0.954. The outcomes of this study exhibit the effective utilization of a meta-heuristic algorithm to improve model quality that can be used as a reference to solve several problems regarding data modeling, pattern recognition, data classification, etc.