-
Multiomics-based Outcome Prediction for the Treatment of Brain Metastases with Personalized Ultra-fractionated Stereotactic Adaptive Radiotherapy (PULSAR)
Authors:
Haozhao Zhang,
Michael Dohopolski,
Strahinja Stojadinovic,
Luiza Giuliani,
Soummitra Anand,
Heejung Kim,
Arnold Pompos,
Andrew Godley,
Steve Jiang,
Tu Dan,
Zabi Wardak,
Robert Timmerman,
Hao Peng
Abstract:
Purpose: We aimed to develop a data-driven multiomics approach integrating radiomics, dosiomics, and delta features to predict treatment response at an earlier stage (intra-treatment) for brain metastases (BMs) patients treated with PULSAR. Methods: We conducted a retrospective study of 39 patients with 69 BMs treated with PULSAR. Radiomics, dosiomics, and delta features were extracted from pretre…
▽ More
Purpose: We aimed to develop a data-driven multiomics approach integrating radiomics, dosiomics, and delta features to predict treatment response at an earlier stage (intra-treatment) for brain metastases (BMs) patients treated with PULSAR. Methods: We conducted a retrospective study of 39 patients with 69 BMs treated with PULSAR. Radiomics, dosiomics, and delta features were extracted from pretreatment and intra-treatment MRI scans and dose distributions. Six individual models and an ensemble feature selection (EFS) model were constructed using SVM and evaluated via stratified 5-fold cross-validation. The classification task distinguished lesions with >20% volume reduction at follow-up. We assessed performance metrics including sensitivity, specificity, accuracy, precision, F1 score, and AUC. Various feature extraction and ensemble selection scenarios were explored to enhance model robustness and reduce overfitting. Results: The EFS model, integrating features from pre-treatment radiomics and dosiomics, intra-treatment radiomics, and delta-radiomics, outperformed six individual models. It achieved an AUC of 0.979, accuracy of 0.917, and F1 score of 0.821. Among the top 9 features, six were derived from post-wavelet transformation, and three from original images. The discrete wavelet transform decomposes volumetric images into multi-resolution components, offering a more comprehensive characterization of underlying structures. Conclusion: Our study demonstrated the feasibility of employing a data-driven multiomics approach to predict tumor volume changes in BMs patients with PULSAR. The EFS model demonstrates enhanced performance compared with six individual models, emphasizing the importance of integrating both pretreatment and intra-treatment data. This application holds promise for optimizing BMs management, potentially mitigating risks associated with under- or over-treatment.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Simulating a pulsed power-driven plasma with ideal MHD
Authors:
A. Beresnyak,
A. L. Velikovich,
J. L. Giuliani,
S. L. Jackson,
J. T. Engelbrecht,
A. S. Richardson,
A. Dasgupta
Abstract:
We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a r…
▽ More
We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a recipe for coupling pulsed power generators to the MHD domain by approximating the perfectly insulating vacuum as a light, perfectly conducting, inviscid MHD fluid and discuss the applicability of this counter-intuitive technique. This, much more affordable ideal MHD representation, is particularly useful in situations where a plasma exhibits interesting three-dimensional phenomena, either due to the design of the experiment or due to developing instabilities. We verified that this coupling recipe works by modeling an exactly solvable flux compression generator as well as a self-similar Noh-like solution and demonstrated convergence to the theoretical solution. We also showed examples of simulating complex three-dimensional pulsed power devices with this technique. We release our code implementation to the public.
△ Less
Submitted 6 May, 2022;
originally announced May 2022.
-
Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma
Authors:
Andrey Beresnyak,
Alexander L. Velikovich,
John L. Giuliani,
Arati Dasgupta
Abstract:
The Naval Research Laboratory "Mag Noh problem", described in this paper, is a self-similar magnetized implosion flow, which contains a fast MHD outward propagating shock of constant velocity. We generalize the classic Noh (1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five-parametric, each solution having its own self-simila…
▽ More
The Naval Research Laboratory "Mag Noh problem", described in this paper, is a self-similar magnetized implosion flow, which contains a fast MHD outward propagating shock of constant velocity. We generalize the classic Noh (1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five-parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field, and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$. Our self-similar solutions are indeed realized when we solve the initial value MHD problem with finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence, and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field, and a shock wave.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
Stability of stagnation via an expanding accretion shock wave
Authors:
A. L. Velikovich,
M. Murakami,
B. D. Taylor,
J. L. Giuliani,
S. T. Zalesak,
Y. Iwamoto
Abstract:
Stagnation of a cold plasma streaming to the center or axis of symmetry via an expanding accretion shock wave is ubiquitous in inertial confinement fusion (ICF) and high-energy-density plasma physics, the examples ranging from plasma flows in x-ray-generating Z pinches [Y. Maron et al., Phys. Rev. Lett. 111, 035001 (2013)] to the experiments in support of the recently suggested concept of impact i…
▽ More
Stagnation of a cold plasma streaming to the center or axis of symmetry via an expanding accretion shock wave is ubiquitous in inertial confinement fusion (ICF) and high-energy-density plasma physics, the examples ranging from plasma flows in x-ray-generating Z pinches [Y. Maron et al., Phys. Rev. Lett. 111, 035001 (2013)] to the experiments in support of the recently suggested concept of impact ignition in ICF [H. Azechi et al., Phys. Rev. Lett. 102, 235002 (2009); M. Murakami et al., Nucl. Fusion 54, 054007 (2014)]. Some experimental evidence indicates that stagnation via an expanding shock wave is stable, but its stability has never been studied theoretically. We present such analysis for the stagnation that does not involve a rarefaction wave behind the expanding shock front and is described by the classic ideal-gas Noh solution in spherical and cylindrical geometry. In either case the stagnated flow has been demonstrated to be stable, initial perturbations exhibiting a power-law, oscillatory or monotonic, decay with time for all the eigenmodes. This conclusion has been supported by our simulations done both on a Cartesian grid and on a curvilinear grid in spherical coordinates. Dispersion equation determining the eigenvalues of the problem and explicit formulas for the eigenfunction profiles corresponding to these eigenvalues are presented, making it possible to use the theory for hydro code verification in two and three dimensions.
△ Less
Submitted 16 March, 2016;
originally announced March 2016.