Determination of muscle shape deformations of the tibialis anterior during dynamic contractions using 3D ultrasound.

Purpose: In this paper, we introduce a novel method for determining 3D deformations of the human tibialis anterior (TA) muscle during dynamic movements using 3D ultrasound. Materials and Methods: An existing automated 3D ultrasound system is used for data acquisition, which consists of three moveable axes, along which the probe can move. While the subjects perform continuous plantar- and dorsiflexion movements in two different controlled velocities, the ultrasound probe sweeps cyclically from the ankle to the knee along the anterior shin. The ankle joint angle can be determined using reflective motion capture markers. Since we considered the movement direction of the foot, i.e., active or passive TA, four conditions occur: slow active, slow passive, fast active, fast passive. By employing an algorithm which defines ankle joint angle intervals, i.e., intervals of range of motion (ROM), 3D images of the volumes during movement can be reconstructed. Results: We found constant muscle volumes between different muscle lengths, i.e., ROM intervals. The results show an increase in mean cross-sectional area (CSA) for TA muscle shortening. Furthermore, a shift in maximum CSA towards the proximal side of the muscle could be observed for muscle shortening. We found significantly different maximum CSA values between the fast active and all other conditions, which might be caused by higher muscle activation due to the faster velocity. Conclusion: In summary, we present a method for determining muscle volume deformation during dynamic contraction using ultrasound, which will enable future empirical studies and 3D computational models of skeletal muscles.

Rheological Analysis and Evaluation of Measurement Techniques for Curing Poly(Methyl Methacrylate) Bone Cement in Vertebroplasty.

Vertebroplasty is a minimally invasive surgical procedure used to treat vertebral fractures, which conventionally involves injecting poly(methyl methacrylate) (PMMA) bone cement into the fractured vertebra. A common risk associated with vertebroplasty is cement leaking out of the vertebra during the injection, which may occur due to a lack of understanding of the complex flow behavior. Therefore, experiments to quantify the cement's flow properties are necessary for understanding and proper handling of the bone cement. In this study, we aimed to characterize the behavior of PMMA bone cement in its curing stages to obtain parameters that govern the flow behavior during injection. We used rotational and oscillatory rheometry for our measurements, as well as a custom-made injector setup that replicated a typical vertebroplasty setting. Our results showed that the complex viscoelastic behavior of bone cement is significantly affected by deformations and temperature. We found that the results from rotational tests, often used for characterizing the bone cement, are susceptible to measurement artifacts caused by wall slip and "ridge"-like formations in the test sample. We also found the Cox-Merz rule to be conditionally valid, which affects the use of oscillatory tests to obtain the shear-thinning characteristics of bone cement. Our findings identify important differences in the measured flow behavior of PMMA bone cement when assessed by different rheological methods, an understanding that is crucial for its risk-free usage in downstream medical applications.

Determining a musculoskeletal system's pre-stretched state using continuum-mechanical forward modelling and joint range optimization.

The subject-specific range of motion (RoM) of a musculoskeletal joint system is balanced by pre-tension levels of individual muscles, which affects their contraction capability. Such an inherent pre-tension or pre-stretch of muscles is not measureable with in vivo experiments. Using a 3D continuum mechanical forward simulation approach for motion analysis of the musculoskeletal system of the forearm with 3 flexor and 2 extensor muscles, we developed an optimization process to determine the muscle fibre pre-stretches for an initial arm position, which is given human dataset. We used RoM values of a healthy person to balance the motion in extension and flexion. The performed sensitivity study shows that the fibre pre-stretches of the m. brachialis, m. biceps brachii and m. triceps brachii with 91 % dominate the objective flexion ratio, while m. brachiradialis and m. anconeus amount 7.8 % and 1.2 % . Within the multi-dimensional space of the surrogate model, 3D sub-spaces of primary variables, namely the dominant muscles and the global objective, flexion ratio, exhibit a path of optimal solutions. Within this optimal path, the muscle fibre pre-stretch of two flexors demonstrate a negative correlation, while, in contrast, the primary extensor, m. triceps brachii correlates positively to each of the flexors. Comparing the global optimum with four other designs along the optimal path, we saw large deviations, e.g., up to 15 ∘ in motion and up to 40% in muscle force. This underlines the importance of accurate determination of fibre pre-stretch in muscles, especially, their role in pathological muscular disorders and surgical applications such as free muscle or tendon transfer.

3D ultrasound-based determination of skeletal muscle fascicle orientations.

Architectural parameters of skeletal muscle such as pennation angle provide valuable information on muscle function, since they can be related to the muscle force generating capacity, fiber packing, and contraction velocity. In this paper, we introduce a 3D ultrasound-based workflow for determining 3D fascicle orientations of skeletal muscles. We used a custom-designed automated motor driven 3D ultrasound scanning system for obtaining 3D ultrasound images. From these, we applied a custom-developed multiscale-vessel enhancement filter-based fascicle detection algorithm and determined muscle volume and pennation angle. We conducted trials on a phantom and on the human tibialis anterior (TA) muscle of 10 healthy subjects in plantarflexion (157 ± 7 ∘ ), neutral position (109 ± 7 ∘ , corresponding to neutral standing), and one resting position in between (145 ± 6 ∘ ). The results of the phantom trials showed a high accuracy with a mean absolute error of 0.92 ± 0.59 ∘ . TA pennation angles were significantly different between all positions for the deep muscle compartment; for the superficial compartment, angles are significantly increased for neutral position compared to plantarflexion and resting position. Pennation angles were also significantly different between superficial and deep compartment. The results of constant muscle volumes across the 3 ankle joint angles indicate the suitability of the method for capturing 3D muscle geometry. Absolute pennation angles in our study were slightly lower than recent literature. Decreased pennation angles during plantarflexion are consistent with previous studies. The presented method demonstrates the possibility of determining 3D fascicle orientations of the TA muscle in vivo.

A System for Reproducible 3D Ultrasound Measurements of Skeletal Muscles.

In 3D freehand ultrasound imaging, operator dependent variations in applied forces and movements can lead to errors in the reconstructed images. In this paper, we introduce an automated 3D ultrasound system, which enables acquisitions with controlled movement trajectories by using motors, which electrically move the probe. Due to integrated encoders there is no need of position sensors. An included force control mechanism ensures a constant contact force to the skin. We conducted 8 trials with the automated 3D ultrasound system on 2 different phantoms with 3 force settings and 10 trials on a human tibialis anterior muscle with 2 force settings. For comparison, we also conducted 8 freehand 3D ultrasound scans from 2 operators (4 force settings) on one phantom and 10 with one operator on the tibialis anterior muscle. Both freehand and automated trials showed small errors in volume and length computations of the reconstructions, however the freehand trials showed larger standard deviations. We also computed the thickness of the phantom and the tibialis anterior muscle. We found significant differences in force settings for the operators and higher coefficients of variation for the freehand trials. Overall, the automated 3D ultrasound system shows a high accuracy in reconstruction. Due to the smaller coefficients of variation, the automated 3D ultrasound system enables more reproducible ultrasound examinations than the freehand scanning. Therefore, the automated 3D ultrasound system is a reliable tool for 3D investigations of skeletal muscle.

Postinhibitory excitation in motoneurons can be facilitated by hyperpolarization-activated inward currents: A simulation study.

Postinhibitory excitation is a transient overshoot of a neuron's baseline firing rate following an inhibitory stimulus and can be observed in vivo in human motoneurons. However, the biophysical origin of this phenomenon is still unknown and both reflex pathways and intrinsic motoneuron properties have been proposed. We hypothesized that postinhibitory excitation in motoneurons can be facilitated by hyperpolarization-activated inward currents (h-currents). Using an electrical circuit model, we investigated how h-currents can modulate the postinhibitory response of motoneurons. Further, we analyzed the spike trains of human motor units from the tibialis anterior muscle during reciprocal inhibition. The simulations revealed that the activation of h-currents by an inhibitory postsynaptic potential can cause a short-term increase in a motoneuron's firing probability. This result suggests that the neuron can be excited by an inhibitory stimulus. In detail, the modulation of the firing probability depends on the time delay between the inhibitory stimulus and the previous action potential. Further, the postinhibitory excitation's strength correlates with the inhibitory stimulus's amplitude and is negatively correlated with the baseline firing rate as well as the level of input noise. Hallmarks of h-current activity, as identified from the modeling study, were found in 50% of the human motor units that showed postinhibitory excitation. This study suggests that h-currents can facilitate postinhibitory excitation and act as a modulatory system to increase motoneuron excitability after a strong inhibition.

Mapping the role of oral cavity physiological factors into the viscoelastic model of denture adhesives for numerical implementation.

Physiological parameters of the oral cavity have a profound impact on any restorative solutions designed for edentulous patients including denture adhesives. This study aims to mathematically quantify the influence of three such variables, namely: the temperature, pH, and the swelling of such adhesives under the influence of saliva on its mechanical behavior. The mathematical quantification is further aimed to implement a material model for such adhesives which considers the impact of such physiological factors. The denture adhesive is experimentally investigated by means of rheological steady state frequency sweep tests to obtain the relaxation spectrum of the material. The relaxation behavior is measured for a wide range of oral cavity temperatures and pH. Also, the adhesive is hydrated and upon swelling to different levels again tested to understand the impact of swelling on the mechanical behavior. The experimentally measured continuous relaxation spectrum is modeled as a viscoelastic material using a discrete set of points based on the Prony series discretization technique. The relaxation spectrums for various temperatures are compared and the possibility of a time-temperature superposition is explored for the model. Similarly, the measured values of Storage and loss modulus are investigated to understand the role of pH and swelling. The results in this study clearly indicated a horizontal shift in the relaxation behavior with increase in temperature. And hence, the time-temperature shift factor was calculated for the adhesive. The relaxation spectrum also showed a strong correlation with swelling of the adhesive and the pH. The influence of these two parameters were captured into the model based on the relaxation time parameter in the Prony series approach. Based on this study the impact of these parameters could be appreciated on the performance and mechanical behavior of denture adhesives and implemented into a Prony series based viscoelastic material model which can be used with numerical simulations.

High-density magnetomyography is superior to high-density surface electromyography for motor unit decomposition: a simulation study.

Objective.Studying motor units is essential for understanding motor control, the detection of neuromuscular disorders and the control of human-machine interfaces. Individual motor unit firings are currently identifiedin vivoby decomposing electromyographic (EMG) signals. Due to our body's properties and anatomy, individual motor units can only be separated to a limited extent with surface EMG. Unlike electrical signals, magnetic fields do not interact with human tissues. This physical property and the emerging technology of quantum sensors make magnetomyography (MMG) a highly promising methodology. However, the full potential of MMG to study neuromuscular physiology has not yet been explored.Approach.In this work, we performin silicotrials that combine a biophysical model of EMG and MMG with state-of-the-art algorithms for the decomposition of motor units. This allows the prediction of an upper-bound for the motor unit decomposition accuracy.Main results.It is shown that non-invasive high-density MMG data is superior over comparable high-density surface EMG data for the robust identification of the discharge patterns of individual motor units. Decomposing MMG instead of EMG increased the number of identifiable motor units by 76%. Notably, MMG exhibits a less pronounced bias to detect superficial motor units.Significance.The presented simulations provide insights into methods to study the neuromuscular system non-invasively andin vivothat would not be easily feasible by other means. Hence, this study provides guidance for the development of novel biomedical technologies.

Linking cortex and contraction-Integrating models along the corticomuscular pathway.

Computational models of the neuromusculoskeletal system provide a deterministic approach to investigate input-output relationships in the human motor system. Neuromusculoskeletal models are typically used to estimate muscle activations and forces that are consistent with observed motion under healthy and pathological conditions. However, many movement pathologies originate in the brain, including stroke, cerebral palsy, and Parkinson's disease, while most neuromusculoskeletal models deal exclusively with the peripheral nervous system and do not incorporate models of the motor cortex, cerebellum, or spinal cord. An integrated understanding of motor control is necessary to reveal underlying neural-input and motor-output relationships. To facilitate the development of integrated corticomuscular motor pathway models, we provide an overview of the neuromusculoskeletal modelling landscape with a focus on integrating computational models of the motor cortex, spinal cord circuitry, α-motoneurons and skeletal muscle in regard to their role in generating voluntary muscle contraction. Further, we highlight the challenges and opportunities associated with an integrated corticomuscular pathway model, such as challenges in defining neuron connectivities, modelling standardisation, and opportunities in applying models to study emergent behaviour. Integrated corticomuscular pathway models have applications in brain-machine-interaction, education, and our understanding of neurological disease.

A continuum mechanical porous media model for vertebroplasty: Numerical simulations and experimental validation.

The outcome of vertebroplasty is hard to predict due to its dependence on complex factors like bone cement and marrow rheologies. Cement leakage could occur if the procedure is done incorrectly, potentially causing adverse complications. A reliable simulation could predict the patient-specific outcome preoperatively and avoid the risk of cement leakage. Therefore, the aim of this work was to introduce a computationally feasible and experimentally validated model for simulating vertebroplasty. The developed model is a multiphase continuum-mechanical macro-scale model based on the Theory of Porous Media. The related governing equations were discretized using a combined finite element-finite volume approach by the so-called Box discretization. Three different rheological upscaling methods were used to compare and determine the most suitable approach for this application. For validation, a benchmark experiment was set up and simulated using the model. The influence of bone marrow and parameters like permeability, porosity, etc., was investigated to study the effect of varying conditions on vertebroplasty. The presented model could realistically simulate the injection of bone cement in porous materials when used with the correct rheological upscaling models, of which the semi-analytical averaging of the viscosity gave the best results. The marrow viscosity is identified as the crucial reference to categorize bone cements as 'high- 'or 'low-' viscosity in the context of vertebroplasty. It is confirmed that a cement with higher viscosity than the marrow ensures stable development of the injection and a proper cement interdigitation inside the vertebra.

Modelling motor units in 3D: influence on muscle contraction and joint force via a proof of concept simulation.

Functional heterogeneity is a skeletal muscle's ability to generate diverse force vectors through localised motor unit (MU) recruitment. Existing 3D macroscopic continuum-mechanical finite element (FE) muscle models neglect MU anatomy and recruit muscle volume simultaneously, making them unsuitable for studying functional heterogeneity. Here, we develop a method to incorporate MU anatomy and information in 3D models. Virtual fibres in the muscle are grouped into MUs via a novel "virtual innervation" technique, which can control the units' size, shape, position, and overlap. The discrete MU anatomy is then mapped to the FE mesh via statistical averaging, resulting in a volumetric MU distribution. Mesh dependency is investigated using a 2D idealised model and revealed that the amount of MU overlap is inversely proportional to mesh dependency. Simultaneous recruitment of a MU's volume implies that action potentials (AP) propagate instantaneously. A 3D idealised model is used to verify this assumption, revealing that neglecting AP propagation results in a slightly less-steady force, advanced in time by approximately 20 ms, at the tendons. Lastly, the method is applied to a 3D, anatomically realistic model of the masticatory system to demonstrate the functional heterogeneity of masseter muscles in producing bite force. We found that the MU anatomy significantly affected bite force direction compared to bite force magnitude. MU position was much more efficacious in bringing about bite force changes than MU overlap. These results highlight the relevance of MU anatomy to muscle function and joint force, particularly for muscles with complex neuromuscular architecture.

Multi-scale mechanobiological model for skeletal muscle hypertrophy.

Skeletal muscle adaptation is correlated to training exercise by triggering different signaling pathways that target many functions; in particular, the IGF1-AKT pathway controls protein synthesis and degradation. These two functions regulate the adaptation in size and strength of muscles. Computational models for muscle adaptation have focused on: the biochemical description of signaling pathways or the mechanical description of muscle function at organ scale; however, an interrelation between these two models should be considered to understand how an adaptation in muscle size affects the protein synthesis rate. In this research, a dynamical model for the IGF1-AKT signaling pathway is linked to a continuum-mechanical model describing the active and passive mechanical response of a muscle; this model is used to study the impact of the adaptive muscle geometry on the protein synthesis at the fiber scale. This new computational model links the signaling pathway to the mechanical response by introducing a growth tensor, and links the mechanical response to the signaling pathway through the evolution of the protein synthesis rate. The predicted increase in cross sectional area (CSA) due to an 8 weeks training protocol excellently agreed with experimental data. Further, our results show that muscle growth rate decreases, if the correlation between protein synthesis and CSA is negative. The outcome of this study suggests that multi-scale models coupling continuum mechanical properties and molecular functions may improve muscular therapies and training protocols.

A biophysically guided constitutive law of the musculotendon-complex: modelling and numerical implementation in Abaqus.

BACKGROUND AND OBJECTIVE: Many biomedical, clinical, and industrial applications may benefit from musculoskeletal simulations. Three-dimensional macroscopic muscle models (3D models) can more accurately represent muscle architecture than their 1D (line-segment) counterparts. Nevertheless, 3D models remain underutilised in academic, clinical, and commercial environments. Among the reasons for this is a lack of modelling and simulation standardisation, verification, and validation. Here, we strive towards a solution by providing an open-access, characterised, constitutive relation (CR) for 3D musculotendon models. METHODS: The musculotendon complex is modelled following the state-of-the-art active stress approach and is treated as hyperelastic, transversely isotropic, and nearly incompressible. Furthermore, force-length and -velocity relationships are incorporated, and muscle activation is derived from motor-unit information. The CR was implemented within the commercial finite-element software package Abaqus as a user-subroutine. A masticatory system model with left and right masseters was used to demonstrate active and passive movement. RESULTS: The CR was characterised by various experimental data sets and was able to capture a wide variety of passive and active behaviours. Furthermore, the masticatory simulations revealed that joint movement was sensitive to the muscle's in-fibre passive response. CONCLUSIONS: This user-material provides a "plug and play" template for 3D neuro-musculoskeletal finite element modelling. We hope that this reduces modelling effort, fosters exchange, and contributes to the standardisation of such models.

Dynamic 3D Ultrasound Imaging of the Tibialis Anterior Muscle.

Skeletal muscle volume has been mainly investigated under static conditions, i.e. isometric contractions. The aim of our study is to use ultrasound imaging to determine muscle deformation during movement. We used a custom-designed scanning rig to obtain 3D ultrasound images of a subject moving the foot from plantarflexion to dorsiflexion at constant velocity. Using motion capture, we computed the respective angle of the ankle for each frame and collected them in bins based on the measured angle (rounded on the next normal number). For each degree, we used Stradwin for the 3D reconstruction of the respective volume. We found increasing cross-sectional areas for increasing dorsiflexion angles. The proposed method is a promising approach for determining muscle volume during movement. Future studies aim at collecting more data to compute muscle volume and length during contraction and compare the results to isometric measurements.

Modelling intra-muscular contraction dynamics using in silico to in vivo domain translation.

BACKGROUND: Advances in sports medicine, rehabilitation applications and diagnostics of neuromuscular disorders are based on the analysis of skeletal muscle contractions. Recently, medical imaging techniques have transformed the study of muscle contractions, by allowing identification of individual motor units' activity, within the whole studied muscle. However, appropriate image-based simulation models, which would assist the continued development of these new imaging methods are missing. This is mainly due to a lack of models that describe the complex interaction between tissues within a muscle and its surroundings, e.g., muscle fibres, fascia, vasculature, bone, skin, and subcutaneous fat. Herein, we propose a new approach to overcome this limitation. METHODS: In this work, we propose to use deep learning to model the authentic intra-muscular skeletal muscle contraction pattern using domain-to-domain translation between in silico (simulated) and in vivo (experimental) image sequences of skeletal muscle contraction dynamics. For this purpose, the 3D cycle generative adversarial network (cycleGAN) models were evaluated on several hyperparameter settings and modifications. The results show that there were large differences between the spatial features of in silico and in vivo data, and that a model could be trained to generate authentic spatio-temporal features similar to those obtained from in vivo experimental data. In addition, we used difference maps between input and output of the trained model generator to study the translated characteristics of in vivo data. RESULTS: This work provides a model to generate authentic intra-muscular skeletal muscle contraction dynamics that could be used to gain further and much needed physiological and pathological insights and assess and overcome limitations within the newly developed research field of neuromuscular imaging.

Investigating the spatial resolution of EMG and MMG based on a systemic multi-scale model.

While electromyography (EMG) and magnetomyography (MMG) are both methods to measure the electrical activity of skeletal muscles, no systematic comparison between both signals exists. Within this work, we propose a novel in silico model for EMG and MMG and test the hypothesis that MMG surpasses EMG in terms of spatial selectivity, i.e. the ability to distinguish spatially shifted sources. The results show that MMG provides a slightly better spatial selectivity than EMG when recorded directly on the muscle surface. However, there is a remarkable difference in spatial selectivity for non-invasive surface measurements. The spatial selectivity of the MMG components aligned with the muscle fibres and normal to the body surface outperforms the spatial selectivity of surface EMG. Particularly, for the MMG's normal-to-the-surface component the influence of subcutaneous fat is minimal. Further, for the first time, we analyse the contribution of different structural components, i.e. muscle fibres from different motor units and the extracellular space, to the measurable biomagnetic field. Notably, the simulations show that for the normal-to-the-surface MMG component, the contribution from volume currents in the extracellular space and in surrounding inactive tissues, is negligible. Further, our model predicts a surprisingly high contribution of the passive muscle fibres to the observable magnetic field.

Numerical study of the stress state on the oral mucosa and abutment tooth upon insertion of partial dentures in the mandible.

The introduction of a removable partial denture onto the dental arch significantly influences the mechanical stress characteristics of both the jawbone and oral mucosa. The aim of this study was to analyze the stress state caused by biting forces upon insertion of partial dentures into the assembly, and to understand the influence of the resulting contact pressure on its retention behavior. For this purpose, a numerical model of a removable partial denture is proposed based on 3D models developed using computer tomography data of the jawbone and the removable partial denture. The denture system rests on the oral mucosa surface and three abutment teeth. The application of bite forces on the denture generated a stick condition on the loaded regions of the denture-oral mucosa interface, which indicates positive retention of the denture onto the oral mucosa surface. Slip and negative retention were observed in the regions of the contact space that were not directly loaded. The contact pressures observed in the regions of the oral mucosa in contact with the denture were below the clinical pressure pain threshold value for soft tissue, which potentially lowers the risk of pain being experienced by denture users. Further, the variation of the retention behavior and contact pressures across different regions of the denture assembly was observed. Thus, there is a need for adhesives or restraining mechanisms for the denture system in order to avoid bending and deformation of sections of the denture as a consequence of the applied bite force.

Muscle-Specific High-Density Electromyography Arrays for Hand Gesture Classification.

OBJECTIVE: Dexterous hand motion is critical for object manipulation. Electrophysiological studies of the hand are key to understanding its underlying mechanisms. High-density electromyography (HD-EMG) provides spatio-temporal information about the underlying electrical activity of muscles, which can be used in neurophysiological research, rehabilitation and control applications. However, existing EMG electrodes platforms are not muscle-specific, which makes the assessment of intrinsic hand muscles difficult. METHODS: Muscle-specific flexible HD-EMG electrode arrays were developed to capture intrinsic hand muscle myoelectric activity during manipulation tasks. The arrays consist of 60 individual electrodes targeting 10 intrinsic hand muscles. Myoelectric activity was displayed as spatio-temporal amplitude maps to visualize muscle activation. Time-domain and temporal-spatial HD-EMG features were extracted to train cubic support vector machine machine-learning classifiers to classify the intended user motion. RESULTS: Experimental data was collected from 5 subjects performing a range of 10 common hand motions. Spatio-temporal EMG maps showed distinct activation areas correlated to the muscles recruited during each movement. The thenar muscle fiber conduction velocity (CV) was estimated to be at 4.7±0.3 m/s for all subjects. Hand motions were successfully classified and average accuracy for all subjects was directly related to spatial resolution based on the number of channels used as inputs; ranging from 74±4% when using only 5 channels and up to 92±2% when using 41 channels. Temporal-spatial features were shown to provide increased motion-specific accuracy when similar muscles were recruited for different gestures. CONCLUSIONS: Muscle-specific electrodes were capable of accurately recording HD-EMG signals from intrinsic hand muscles and accurately predicting motion. SIGNIFICANCE: The muscle-specific electrode arrays could improve electrophysiological research studies using EMG decomposition techniques to assess motor unit activity and in applications involving the analysis of dexterous hand motions.

Optimizing NV magnetometry for Magnetoneurography and Magnetomyography applications.

Magnetometers based on color centers in diamond are setting new frontiers for sensing capabilities due to their combined extraordinary performances in sensitivity, bandwidth, dynamic range, and spatial resolution, with stable operability in a wide range of conditions ranging from room to low temperatures. This has allowed for its wide range of applications, from biology and chemical studies to industrial applications. Among the many, sensing of bio-magnetic fields from muscular and neurophysiology has been one of the most attractive applications for NV magnetometry due to its compact and proximal sensing capability. Although SQUID magnetometers and optically pumped magnetometers (OPM) have made huge progress in Magnetomyography (MMG) and Magnetoneurography (MNG), exploring the same with NV magnetometry is scant at best. Given the room temperature operability and gradiometric applications of the NV magnetometer, it could be highly sensitive in the pT / Hz -range even without magnetic shielding, bringing it close to industrial applications. The presented work here elaborates on the performance metrics of these magnetometers to the state-of-the-art techniques by analyzing the sensitivity, dynamic range, and bandwidth, and discusses the potential benefits of using NV magnetometers for MMG and MNG applications.