Human skeletal muscle energy metabolism: when a physiological model promotes the search for new technologies

The idea that muscle function and, in fact, the function of all living tissues may be described by physical laws appeared in the second half of the 19th century. During this period, Helmoltz (Über die Erhaltung der Kraft, 1847) showed that the principle of energy conservation may be applied to living systems. In the field of exercise physiology, this idea has been subsequently developed by Rodolfo Margaria's (1901-1983) School and the mathematical formalization of the theory has succeeded in the bioenergetic model. During the last 20years methodical and critical study of the bioenergetic model has been carried out by Paolo Cerretelli, one of the most important heirs of Rodolfo Margaria's School. Original results and technological developments have been generated by his activity and many young scientists have been educated in this approach. The present paper wants to present the modern history of the bioenergetic model and is dedicated to Paolo Cerretelli on the occasion of his seventieth birthda

Probability density function for random photon steps in a binary (isotropic-Poisson) statistical mixture

Monte Carlo (MC) simulations allowing to describe photons propagation in statistical mixtures represent an interest that goes way beyond the domain of optics, and can cover, e.g., nuclear reactor physics, image analysis or life science just to name a few. MC simulations are considered a ``gold standard'' because they give exact solutions (in the statistical sense), however, in the case of statistical mixtures they are enormously time consuming and their implementation is often extremely complex. For this reason, the aim of the present contribution is to propose a new approach that should allow us in the future to simplify the MC approach. This is done through an explanatory example, i.e.; by deriving the `exact' analytical expression for the probability density function of photons' random steps (single step function, SSF) propagating in a medium represented as a binary (isotropic-Poisson) statistical mixture. The use of the SSF reduces the problem to an `equivalent' homogeneous medium behaving exactly as the original binary statistical mixture. This will reduce hundreds time-consuming MC simulations to only one equivalent simple MC simulation. To the best of our knowledge the analytically `exact' SSF for a binary (isotropic-Poisson) statistical mixture has never been derived before

Topological complexity of photons' paths in biological tissues

In the present contribution three means of measuring the geometrical and topological complexity of photons' paths in random media are proposed. This is realized by investigating the behavior of the average crossing number, the mean writhe, and the minimal crossing number of photons' paths generated by Monte Carlo (MC) simulations, for different sets of optical parameters. It is observed that the complexity of the photons' paths increases for increasing light source/detector spacing, and that highly "knotted" paths are formed. Due to the particular rules utilized to generate the MC photons' paths, the present results may have an interest not only for the biomedical optics community, but also from a pure mathematical point of view

Use of near-infrared systems for investigations of hemodynamics in human in vivo bone tissue: a systematic review

This is the author accepted manuscript. The final version is available from Wiley via the DOI in this record.A range of technologies using near infrared (NIR) light have shown promise at providing real time measurements of hemodynamic markers in bone tissue in vivo, an exciting prospect given existing difficulties in measuring hemodynamics in bone tissue. This systematic review aimed to evaluate the evidence for this potential use of NIR systems, establishing their potential as a research tool in this field. Major electronic databases including MEDLINE and EMBASE were searched using pre‐planned search strategies with broad scope for any in vivo use of NIR technologies in human bone tissue. Following identification of studies by title and abstract screening, full text inclusion was determined by double blind assessment using predefined criteria. Full text studies for inclusion were data extracted using a predesigned proforma and quality assessed. Narrative synthesis was appropriate given the wide heterogeneity of included studies. Eighty‐eight full text studies fulfilled the inclusion criteria, 57 addressing laser Doppler flowmetry (56 intra‐operatively), 21 near infrared spectroscopy, and 10 photoplethysmography. The heterogeneity of the methodologies included differing hemodynamic markers, measurement protocols, anatomical locations, and research applications, making meaningful direct comparisons impossible. Further, studies were often limited by small sample sizes with potential selection biases, detection biases, and wide variability in results between participants. Despite promising potential in the use of NIR light to interrogate bone circulation, the application of NIR systems in bone requires rigorous assessment of the reproducibility of potential hemodynamic markers and further validation of these markers against alternative physiologically relevant reference standards.This systematic review was supported by the College of Radiographers Industry Partnership Scheme (CORIPS) Doctoral Fellowship Grant (Applicant 003). The CORIPS are providing financial support but have no input into the design, performance or analysis of this systematic review. WDS, FC and CT would like to acknowledge the NIHR Exeter Clinical Research Facility and the NIHR Collaboration for Leadership in Applied Health Research and Care (CLAHRC) for the South West Peninsula. The views expressed in this publication are those of the author(s) and not necessarily those of the NIHR Exeter Clinical Research Facility, the NHS, the NIHR or the Department of Health in England

There's plenty of light at the bottom: Statistics of photon penetration depth in random media

We propose a comprehensive statistical approach describing the penetration depth of light in random media. The presented theory exploits the concept of probability density function f(z|ρ, t) for the maximum depth reached by the photons that are eventually re-emitted from the surface of the medium at distance ρ and time t. Analytical formulas for f, for the mean maximum depth 〈zmax〉 and for the mean average depth 〈z〉 reached by the detected photons at the surface of a diffusive slab are derived within the framework of the diffusion approximation to the radiative transfer equation, both in the time domain and the continuous wave domain. Validation of the theory by means of comparisons with Monte Carlo simulations is also presented. The results are of interest for many research fields such as biomedical optics, advanced microscopy and disordered photonics

Time-domain Raman analytical forward solvers

A set of time-domain analytical forward solvers for Raman signals detected from homogeneous diffusive media is presented. The time-domain solvers have been developed for two geometries: the parallelepiped and the finite cylinder. The potential presence of a background fluorescence emission, contaminating the Raman signal, has also been taken into account. All the solvers have been obtained as solutions of the time dependent diffusion equation. The validation of the solvers has been performed by means of comparisons with the results of "gold standard" Monte Carlo simulations. These forward solvers provide an accurate tool to explore the information content encoded in the time-resolved Raman measurements

The Influence of Maturity Status on Muscle Architecture in School-Aged Boys

Purpose: To determine the differences in muscle architecture of the lower limb in pre-peak height velocity (PHV), circa-PHV, and post-PHV boys. Method: Muscle architecture variables from both the gastrocnemius medialis (GM) and vastus lateralis (VL) were derived from ultrasonographic images in 126 school-aged boys. One-way analysis of variance using Bonferroni post hoc comparisons was employed to determine between-group differences, and effect sizes were calculated to establish the magnitude of these differences. Results: All muscle architecture variables showed significant small to large increases from pre-PHV to post-PHV, excluding GM fascicle length (d = 0.59–1.39; P  0.57; P < .05); however, only the VL muscle thickness and physiological thickness increased from circa-PHV to post-PHV (d = 0.68; P < .05). The post-PHV group also showed larger GM pennation angles than the circa-PHV group (d = 0.59; P < .05). Conclusion: The combined results showed that maturation is associated with changes in muscle morphology. These data quantify that the maturity-related changes in muscle architecture variables provide a reference to differentiate between training-induced adaptations versus changes associated with normal growth and maturation

Two-step verification method for Monte Carlo codes in biomedical optics applications

Significance: Code verification is an unavoidable step prior to using a Monte Carlo (MC) code. Indeed, in biomedical optics, a widespread verification procedure for MC codes is still missing. Analytical benchmarks that can be easily used for the verification of different MC routines offer an important resource.Aim: We aim to provide a two-step verification procedure for MC codes enabling the two main tasks of an MC simulator: (1) the generation of photons' trajectories and (2) the intersections of trajectories with boundaries separating the regions with different optical properties. The proposed method is purely based on elementary analytical benchmarks, therefore, the correctness of an MC code can be assessed with a one-sample t-test.Approach: The two-step verification is based on the following two analytical benchmarks: (1) the exact analytical formulas for the statistical moments of the spatial coordinates where the scattering events occur in an infinite medium and (2) the exact invariant solutions of the radiative transfer equation for radiance, fluence rate, and mean path length in media subjected to a Lambertian illumination.Results: We carried out a wide set of comparisons between MC results and the two analytical benchmarks for a wide range of optical properties (from non-scattering to highly scattering media, with different types of scattering functions) in an infinite non-absorbing medium (step 1) and in a non-absorbing slab (step 2). The deviations between MC results and exact analytical values are usually within two standard errors (i.e., t-tests not rejected at a 5% level of significance). The comparisons show that the accuracy of the verification increases with the number of simulated trajectories so that, in principle, an arbitrary accuracy can be obtained.Conclusions: Given the simplicity of the verification method proposed, we envision that it can be widely used in the field of biomedical optics.</p

Depth sensitivity of frequency domain optical measurements in diffusive media

The depth sensitivity functions for AC amplitude, phase (PH) and DC intensity signals have been obtained in the frequency domain (where the source amplitude is modulated at radio-frequencies) by making use of analytical solutions of the photon diffusion equation in an infinite slab geometry. Furthermore, solutions for the relative contrast of AC, PH and DC signals when a totally absorbing plane is placed at a fixed depth of the slab have also been obtained. The solutions have been validated by comparisons with gold standard Monte Carlo simulations. The obtained results show that the AC signal, for modulation frequencies < 200 MHz, has a depth sensitivity with similar characteristics to that of the continuous-wave (CW) domain (source modulation frequency of zero). Thus, the depth probed by such a signal can be estimated by using the formula of penetration depth for the CW domain (Sci. Rep. 6, 27057 (2016)). However, the PH signal has a different behavior compared to the CW domain, showing a larger depth sensitivity at shallow depths and a less steep relative contrast as a function of depth. These results mark a clear difference in term of depth sensitivity between AC and PH signals, and highlight the complexity of the estimation of the actual depth probed in tissue spectroscopy.Postprint (published version

Comparison of NIRS, laser Doppler flowmetry, photoplethysmography, and pulse oximetry during vascular occlusion challenges

© 2016 Institute of Physics and Engineering in Medicine. Monitoring changes in blood volume, blood flow, and oxygenation in tissues is of vital importance in fields such as reconstructive surgery and trauma medicine. Near infrared spectroscopy (NIRS), laser Doppler (LDF) flowmetry, photoplethysmography (PPG), and pulse oximetry (PO) contribute to such fields due to their safe and noninvasive nature. However, the techniques have been rarely investigated simultaneously or altogether. The aim of this study was to investigate all the techniques simultaneously on healthy subjects during vascular occlusion challenges. Sensors were attached on the forearm (NIRS and LDF) and fingers (PPG and PO) of 19 healthy volunteers. Different degrees of vascular occlusion were induced by inflating a pressure cuff on the upper arm. The responses of tissue oxygenation index (NIRS), tissue haemoglobin index (NIRS), flux (LDF), perfusion index (PPG), and arterial oxygen saturation (PO) have been recorded and analyzed. Moreover, the optical densities were calculated from slow varying dc PPG, in order to distinguish changes in venous blood volumes. The indexes showed significant changes (p < 0.05) in almost all occlusions, either venous or over-systolic occlusions. However, differentiation between venous and arterial occlusion by LDF may be challenging and the perfusion index (PI) may not be adequate to indicate venous occlusions. Optical densities may be an additional tool to detect venous occlusions by PPG