CyberSec Research

Skin tension plays a pivotal role in clinical settings, it affects scarring, wound healing and skin necrosis. Despite its importance, there is no widely accepted method for assessing in vivo skin tension or its natural pre-stretch. This study aims to utilise modern machine learning (ML) methods to develop a model that uses non-invasive measurements of surface wave speed to predict clinically useful skin properties such as stress and natural pre-stretch. A large dataset consisting of simulated wave propagation experiments was created using a simplified two-dimensional finite element (FE) model. Using this dataset, a sensitivity analysis was performed, highlighting the effect of the material parameters and material model on the Rayleigh and supersonic shear wave speeds. Then, a Gaussian process regression model was trained to solve the ill-posed inverse problem of predicting stress and pre-stretch of skin using measurements of surface wave speed. This model had good predictive performance (R2 = 0.9570) and it was possible to interpolate simplified parametric equations to calculate the stress and pre-stretch. To demonstrate that wave speed measurements could be obtained cheaply and easily, a simple experiment was devised to obtain wave speed measurements from synthetic skin at different values of pre-stretch. These experimental wave speeds agree well with the FE simulations and a model trained solely on the FE data provided accurate predictions of synthetic skin stiffness. Both the simulated and experimental results provide further evidence that elastic wave measurements coupled with ML models are a viable non-invasive method to determine in vivo skin tension.

We investigate subset-based optimization methods for positron emission tomography (PET) image reconstruction incorporating a regularizing prior. PET reconstruction methods that use a prior, such as the relative difference prior (RDP), are of particular relevance, as they are widely used in clinical practice and have been shown to outperform conventional early-stopped and post-smoothed ordered subsets expectation maximization (OSEM). Our study evaluates these methods on both simulated data and real brain PET scans from the 2024 PET Rapid Image Reconstruction Challenge (PETRIC), where the main objective was to achieve RDP-regularized reconstructions as fast as possible, making it an ideal benchmark. Our key finding is that incorporating the effect of the prior into the preconditioner is crucial for ensuring fast and stable convergence. In extensive simulation experiments, we compare several stochastic algorithms -- including Stochastic Gradient Descent (SGD), Stochastic Averaged Gradient Amelior\'e (SAGA), and Stochastic Variance Reduced Gradient (SVRG) -- under various algorithmic design choices and evaluate their performance for varying count levels and regularization strengths. The results show that SVRG and SAGA outperformed SGD, with SVRG demonstrating a slight overall advantage. The insights gained from these simulations directly contributed to the design of our submitted algorithms, which formed the basis of the winning contribution to the PETRIC 2024 challenge.

Infrared thermography has gained interest as a tool for non-contact measurement of blood circulation and skin blood flow due to cardiac activity. Partiularly, blood vessels on the surface, such as on the back of the hand, are suited for visualization. However, standardized methodologies have not yet been established for areas such as the face and neck, where many blood vessels are lie deeper beneath the surface, and external stimulation for measurement could be harmful. Here we propose Synchro-Thermography for stable monitoring of facial temperature changes associated with heart rate variability. We conducted experiments with eight subjects and measured minute temperature changes with an amplitude of about \SI{10}{mK} on the forehead and chin. The proposed method improves the temperature resolution by a factor of 2 or more, and can stably measure skin temperature changes caused by blood flow. This skin temperature change could be applied to physiological sensing such as blood flow changes due to injury or disease, or as an indicator of stress.

Ultra-low-field (ULF) MRI is emerging as an alternative modality to high-field (HF) MRI due to its lower cost, minimal siting requirements, portability, and enhanced accessibility factors that enable large-scale deployment. Although ULF-MRI exhibits lower signal-to-noise ratio (SNR), advanced imaging and data-driven denoising methods enabled by high-performance computing have made contrasts like diffusion-weighted imaging (DWI) feasible at ULF. This study investigates the potential and limitations of ULF tractography, using data acquired on a 0.064 T commercially available mobile point-of-care MRI scanner. The results demonstrate that most major white matter bundles can be successfully retrieved in healthy adult brains within clinically tolerable scan times. This study also examines the recovery of diffusion tensor imaging (DTI)-derived scalar maps, including fractional anisotropy and mean diffusivity. Strong correspondence is observed between scalar maps obtained with ULF-MRI and those acquired at high field strengths. Furthermore, fibre orientation distribution functions reconstructed from ULF data show good agreement with high-field references, supporting the feasibility of using ULF-MRI for reliable tractography. These findings open new opportunities to use ULF-MRI in studies of brain health, development, and disease progression particularly in populations traditionally underserved due to geographic or economic constraints. The results show that robust assessments of white matter microstructure can be achieved with ULF-MRI, effectively democratising microstructural MRI and extending advanced imaging capabilities to a broader range of research and clinical settings where resources are typically limited.

Purpose: Intensity-modulated proton therapy (IMPT) offers precise tumor coverage while sparing organs at risk (OARs) in head and neck (H&N) cancer. However, its sensitivity to anatomical changes requires frequent adaptation through online adaptive radiation therapy (oART), which depends on fast, accurate dose calculation via Monte Carlo (MC) simulations. Reducing particle count accelerates MC but degrades accuracy. To address this, denoising low-statistics MC dose maps is proposed to enable fast, high-quality dose generation. Methods: We developed a diffusion transformer-based denoising framework. IMPT plans and 3D CT images from 80 H&N patients were used to generate noisy and high-statistics dose maps using MCsquare (1 min and 10 min per plan, respectively). Data were standardized into uniform chunks with zero-padding, normalized, and transformed into quasi-Gaussian distributions. Testing was done on 10 H&N, 10 lung, 10 breast, and 10 prostate cancer cases, preprocessed identically. The model was trained with noisy dose maps and CT images as input and high-statistics dose maps as ground truth, using a combined loss of mean square error (MSE), residual loss, and regional MAE (focusing on top/bottom 10% dose voxels). Performance was assessed via MAE, 3D Gamma passing rate, and DVH indices. Results: The model achieved MAEs of 0.195 (H&N), 0.120 (lung), 0.172 (breast), and 0.376 Gy[RBE] (prostate). 3D Gamma passing rates exceeded 92% (3%/2mm) across all sites. DVH indices for clinical target volumes (CTVs) and OARs closely matched the ground truth. Conclusion: A diffusion transformer-based denoising framework was developed and, though trained only on H&N data, generalizes well across multiple disease sites.

Among the genetic algorithms generally used for optimization problems in the recent decades, quantum-inspired variants are known for fast and high-fitness convergence and small resource requirement. Here the application to the patient scheduling problem in proton therapy is reported. Quantum chromosomes are tailored to possess the superposed data of patient IDs and gantry statuses. Selection and repair strategies are also elaborated for reliable convergence to a clinically feasible schedule although the employed model is not complex. Clear advantage in population size is shown over the classical counterpart in our numerical results for both a medium-size test case and a large-size practical problem instance. It is, however, observed that program run time is rather long for the large-size practical case, which is due to the limitation of classical emulation and demands the forthcoming true quantum computation. Our results also revalidate the stability of the conventional classical genetic algorithm.

Background: Mechanical Thrombectomy (MT) is a widely accepted first-line treatment for Acute Ischemic Stroke (AIS) and it has been studied using in vitro and in silico models. Thrombectomy outcomes have been performed for patient-specific cases using in silico models. However, until now, in vivo friction coefficients for stent-vessel, stent-clot, and clot-vessel interactions are unknown, but in vitro experiments have been attempted with significant standard deviations. These interactions and friction coefficients have been considered an important aspect of thrombectomy success. Objectives: In the current study, we explored the influence of variation in friction forces for stent-vessel, stent-clot, and clot-vessel interactions using virtual mechanical thrombectomy (VMT). We have performed three simulations for each interaction and varied friction coefficients around the standard deviation observed in the past in vitro studies. Results: (i) clot-vessel friction: higher friction leads to clot fragmentation and VMT failure. (ii) stent-clot friction: it is susceptible to VMT outcomes, with lower values showing the slippage of the clot while higher values lead to fragmentation. (iii) stent-vessel friction: higher friction shows compression of the stent in curved vessels and dislodgment of clot from stent retriever (SR) due to its compression, which leads to VMT failure. (iv) retrieval speed (RS): higher RS (>30 mm/s) leads to significant stent compression and unrealistic behavior of the SR. Conclusions: Analysis of results proposes the necessity for calculating accurate friction factor values and their implementation into in silico models, due to their sensitivity towards thrombectomy outcomes. Such in silico models mimic in vivo thrombectomy more closely and can be used in mechanical thrombectomy planning, management, and decision-making.

Recent work has shown improved lesion detectability and flexibility to reconstruction hyperparameters (e.g. scanner geometry or dose level) when PET images are reconstructed by leveraging pre-trained diffusion models. Such methods train a diffusion model (without sinogram data) on high-quality, but still noisy, PET images. In this work, we propose a simple method for generating subject-specific PET images from a dataset of multi-subject PET-MR scans, synthesizing "pseudo-PET" images by transforming between different patients' anatomy using image registration. The images we synthesize retain information from the subject's MR scan, leading to higher resolution and the retention of anatomical features compared to the original set of PET images. With simulated and real [$^{18}$F]FDG datasets, we show that pre-training a personalized diffusion model with subject-specific "pseudo-PET" images improves reconstruction accuracy with low-count data. In particular, the method shows promise in combining information from a guidance MR scan without overly imposing anatomical features, demonstrating an improved trade-off between reconstructing PET-unique image features versus features present in both PET and MR. We believe this approach for generating and utilizing synthetic data has further applications to medical imaging tasks, particularly because patient-specific PET images can be generated without resorting to generative deep learning or large training datasets.

Time-resolved CT is an advanced measurement technique that has been widely used to observe dynamic objects, including periodically varying structures such as hearts, lungs, or hearing structures. To reconstruct these objects from CT projections, a common approach is to divide the projections into several collections based on their motion phases and perform reconstruction within each collection, assuming they originate from a static object. This describes the gating-based method, which is the standard approach for time-periodic reconstruction. However, the gating-based reconstruction algorithm only utilizes a limited subset of projections within each collection and ignores the correlation between different collections, leading to inefficient use of the radiation dose. To address this issue, we propose two analytical reconstruction pipelines in this paper, and validate them with experimental data captured using tomographic synchrotron microscopy. We demonstrate that our approaches significantly reduce random noise in the reconstructed images without blurring the sharp features of the observed objects. Equivalently, our methods can achieve the same reconstruction quality as gating-based methods but with a lower radiation dose. Our code is available at github.com/PeriodRecon.

In order to obtain accurate contact parameters for the discrete element simulation of salt particles used in animal husbandry, the principle of particle contact scaling and dimensional analysis were used for particle scaling. Firstly, the Plackett Burman experiment was used to screen the parameters that significantly affect the angle of repose: salt salt rolling friction coefficient, salt salt recovery coefficient, and salt steel rolling friction coefficient. Considering the influence of other parameters, a combination of bench and simulation experiments was used to calibrate the contact parameters between salt particles and steel plates used in animal husbandry in EDEM. Finally, through the stacking test, steepest climbing test, and orthogonal rotation combination test, the salt salt rolling friction coefficient was obtained to be 0.23, the salt salt recovery coefficient was 0.544, and the salt steel rolling friction coefficient was 0.368, which were verified through bench tests. The experimental results show that the relative error between the actual value of the stacking angle and the simulation results is 0.6%. The results indicate that the calibrated contact parameters can be used for discrete element simulation of salt particles for animal husbandry, providing reference for the design of quantitative feeding screws and silos.

X-ray imaging, traditionally relying on attenuation contrast, struggles to differentiate materials with similar attenuation coefficients like soft tissues. X-ray phase contrast imaging (XPCI) and dark-field (DF) imaging provide enhanced contrast by detecting phase shifts and ultra-small-angle X-ray scattering (USAXS). However, they typically require complex and costly setups, along with multiple exposures to retrieve various contrast features. In this study, we introduce a novel single-mask X-ray imaging system design that simultaneously captures attenuation, differential phase contrast (DPC), and dark-field images in a single exposure. Most importantly, our proposed system design requires just a single mask alignment with relatively low-resolution detectors. Using our novel light transport models derived for these specific system designs, we show intuitive understanding of contrast formation and retrieval method of different contrast features. Our approach eliminates the need for highly coherent X-ray sources, ultra-high-resolution detectors, spectral detectors or intricate gratings. We propose three variations of the single-mask setup, each optimized for different contrast types, offering flexibility and efficiency in a variety of applications. The versatility of this single-mask approach along with the use of befitting light transport models holds promise for broader use in clinical diagnostics and industrial inspection, making advanced X-ray imaging more accessible and cost-effective.

Craniosynostosis is a medical condition that affects the growth of babies' heads, caused by an early fusion of cranial sutures. In recent decades, surgical treatments for craniosynostosis have significantly improved, leading to reduced invasiveness, faster recovery, and less blood loss. At Great Ormond Street Hospital (GOSH), the main surgical treatment for patients diagnosed with sagittal craniosynostosis (SC) is spring assisted cranioplasty (SAC). This procedure involves a 15x15 mm2 osteotomy, where two springs are inserted to induce distraction. Despite the numerous advantages of this surgical technique for patients, the outcome remains unpredictable due to the lack of efficient preoperative planning tools. The surgeon's experience and the baby's age are currently relied upon to determine the osteotomy location and spring selection. Previous tools for predicting the surgical outcome of SC relied on finite element modeling (FEM), which involved computed tomography (CT) imaging and required engineering expertise and lengthy calculations. The main goal of this research is to develop a real-time prediction tool for the surgical outcome of patients, eliminating the need for CT scans to minimise radiation exposure during preoperative planning. The proposed methodology involves creating personalised synthetic skulls based on three-dimensional (3D) photographs, incorporating population average values of suture location, skull thickness, and soft tissue properties. A machine learning (ML) surrogate model is employed to achieve the desired surgical outcome. The resulting multi-output support vector regressor model achieves a R2 metric of 0.95 and MSE and MAE below 0.13. Furthermore, in the future, this model could not only simulate various surgical scenarios but also provide optimal parameters for achieving a maximum cranial index (CI).

Precise radiation delivery is critical for effective radiotherapy, and gold nanoparticles (AuNPs) have emerged as promising tools to enhance local dose deposition while sparing the surrounding healthy tissue. In this study, the PENELOPE Monte Carlo code was used to investigate the dosimetry of AuNPs under different conditions and models. The Dose Enhancement Ratio (DER) was studied in water and breast tissue with spherical shapes and in agreement with previously published results. To further analyse the physical interactions of the particles around the AuNP, a Phase Space File (PSF) in a volume around the AuNPs was created. This showed that larger AuNPs lead to increased doses, as expected, yielding DER values exceeding 100 times. Finally, results reveal that in the volume surrounding the AuNP, 80% of emitted electrons originate from photoelectric absorption, leading to Auger electron emission cascades which were analysed in detail. It was also possible to establish a direct relation between number of secondaries and the particle volumes. The Local Effect Model (LEM) was used to determine survival curves in AuNPs of different sizes at different gold concentrations. The last part of this work consisted in analysing a distribution of AuNPs within a flattened cell typical of clonogenic assays where a log-normal distribution of dose was observed. This led to the development of a new, mechanistic, Local Effect Model which, if further validated, can have further applications in-vitro and in-silico.

Psoriasis is a long-term inflammatory skin disease that remains difficult to treat. In this study, we developed a new topical treatment by combining metal oxide nanoparticles: cerium oxide (CeO2), zinc oxide (ZnO), and silver (Ag), with natural plant extracts in a gel made from fish collagen and agar. The nanoparticles were characterized using UV-Vis spectroscopy, dynamic light scattering (DLS), Fourier-transform infrared spectroscopy (FTIR), and scanning electron microscopy (SEM), showing good stability and a uniform particle size distribution (ZnO averaged 66 nm). To enhance therapeutic potential, the gel was enriched with plant-derived antioxidants from bitter melon, ginger, and neem. This formulation was tested on an animal model of psoriasis. The treated group exhibited faster wound healing and reduced inflammation compared to both placebo and untreated groups, with statistically significant results (p < 0.01 to p < 0.001) observed from Day 3, becoming more pronounced by Day 14. These results indicate that the combination of nanoparticles with plant-based components in a topical gel may provide a promising new approach to psoriasis treatment. Further studies are recommended to evaluate long-term safety and therapeutic effectiveness.

We present a stabilised finite element method for modelling proton transport in tissue, incorporating both inelastic energy loss and elastic angular scattering. A key innovation is a positivity-preserving formulation that guarantees non-negative fluence and dose, even on coarse meshes. This enables reliable computation of clinically relevant quantities for treatment planning. We derive a priori error estimates demonstrating optimal convergence rates and validate the method through numerical benchmarks. The proposed framework provides a robust, accurate and efficient tool for advancing proton beam therapy.

Noninvasive imaging deep into the adult brain at submillimeter and millisecond scales remains a challenge in medical imaging. Here, we report a helmet based ultrasound brain imager built from a customized helmet, a scanned ultrasound array, and three dimensional printing for real time imaging of human brain anatomical and functional information. Through its application to post hemicraniectomy patients in a sitting position, we achieved volumetric brain tissue structural, vascular, and blood flow images at centimeter scale depths with submillimeter and millisecond spatiotemporal resolutions. We also demonstrated the system capability to track cerebral blood flow over repeated imaging sessions, including during motion prone conditions. Our brain imager circumvents the skull and bridges the gap between high resolution human brain imaging and wearable convenience. This imager may serve as a platform for further investigations into human brain dynamics in post hemicraniectomy patients and offer insights into the brain that could surpass those obtained from non human primate studies.

Magnetic hyperthermia treatment (MHT) utilizes heat generated from magnetic nanoparticles (MNPs) under an alternating magnetic field (AMF) for therapeutic applications. Gadolinium silicide (Gd5Si4) has emerged as a promising MHT candidate due to its self-regulating heating properties and potential biocompatibility. However, the impact of high-dose X-ray irradiation on its magnetic behavior remains uncertain. This study examines Gd5Si4 nanoparticles exposed to 36 and 72 kGy X-ray irradiation at a high-dose rate (120 Gy/min). While X-ray diffraction, scanning electron microscopy, and energy dispersive spectroscopy confirm no structural or compositional changes, transmission electron microscopy reveals localized lattice distortions, along with observable changes in magnetic properties, as evidenced in magnetization vs. temperature and hysteresis measurements. Despite this, magnetocaloric properties and specific loss power (SLP) remain unaffected. Our findings confirm the stability of Gd5Si4 under high-dose X-ray irradiation, supporting its potential for radiotherapy (RT) and magnetocaloric cooling in deep-space applications.

For nonlinear multispectral computed tomography (CT), accurate and fast image reconstruction is challenging when the scanning geometries under different X-ray energy spectra are inconsistent or mismatched. Motivated by this, we propose an accurate and fast algorithm named AFIRE to address such problem in the case of mildly full scan. We discover that the derivative operator (gradient) of the involved nonlinear mapping at some special points, for example, at zero, can be represented as a composition (block multiplication) of a diagonal operator (matrix) composed of X-ray transforms (projection matrices) and a very small-scale matrix. Based on the insights, the AFIRE is proposed respectively from the continuous, discrete and actual-use perspectives by leveraging the simplified Newton method. Under proper conditions, we establish the convergence theory of the proposed algorithm. Furthermore, numerical experiments are also carried out to verify that the proposed algorithm can accurately and effectively reconstruct the basis images in completely geometric-inconsistency dual-energy CT with noiseless and noisy projection data. Particularly, the proposed algorithm significantly outperforms some state-of-the-art methods in terms of accuracy and efficiency. Finally, the flexibility and extensibility of the proposed algorithm are also demonstrated.