Despite the alarming incidence rate and potential for debilitating outcomes of sports-related concussion, the underlying mechanisms of injury remain to be expounded. Since as early as 1950, researchers have aimed to characterize head impact biomechanics through in-lab and in-game investigations. The ever-growing body of literature within this area has supported the inherent connection between head kinematics during impact and injury outcomes. Even so, traditional metrics of peak acceleration, time window, and HIC have outlived their potential. More sophisticated analysis techniques are required to advance the understanding of concussive vs subconcussive impacts. The work presented within this thesis was motivated by the exploration of advanced approaches to 1) experimental theory and design of impact reconstructions and 2) characterization of kinematic profiles for model building. These two areas of investigation resulted in the presentation of refined, systematic approaches to head impact analysis that should begin to replace outdated standards and metrics.

<p>Bridging the accuracy of ab initio (AI) QM/MM with the efficiency of semi-empirical<br> (SE) QM/MM methods has long been a goal in computational chemistry. This dissertation<br> presents four ∆-Machine learning schemes aimed at achieving this objective. Firstly, the in-<br> corporation of negative force observations into the Gaussian process regression (GPR) model,<br> resulting in GPR with derivative observations, demonstrates the remarkable capability to<br> attain high-quality potential energy surfaces, accurate Cartesian force descriptions, and reli-<br> able free energy profiles using a training set of just 80 points. Secondly, the adaptation of the<br> sparse streaming GPR algorithm showcases the potential of memory retention from previous<br> phasespace, enabling energy-only models to converge using simple descriptors while faith-<br> fully reproducing high-quality potential energy surfaces and accurate free energy profiles.<br> Thirdly, the utilization of GPR with atomic environmental vectors as input features proves<br> effective in enhancing both potential energy surface and free energy description. Further-<br> more, incorporating derivative information on solute atoms further improves the accuracy<br> of force predictions on molecular mechanical (MM) atoms, addressing discrepancies arising<br> from QM/MM interaction energies between the target and base levels of theory. Finally, a<br> comprehensive comparison of three distinct GPR schemes, namely GAP, GPR with an aver-<br> age kernel, and GPR with a system-specific sum kernel, is conducted to evaluate the impact<br> of permutational invariance and atomistic learning on the model’s quality. Additionally, this<br> dissertation introduces the adaptation of the GAP method to be compatible with the sparse<br> variational Gaussian processes scheme and the streaming sparse GPR scheme, enhancing<br> their efficiency and applicability. Through these four ∆-Machine learning schemes, this dis-<br> sertation makes significant contributions to the field of computational chemistry, advancing<br> the quest for accurate potential energy surfaces, reliable force descriptions, and informative<br> free energy profiles in QM/MM simulations.<br> </p>