In this paper, we present a 3D virtual phantom design software, which was developed based on object-oriented programming methodology and dedicated to medical physics research. This software was named Magical Phantom (MPhantom), which is composed of 3D visual builder module and virtual CT scanner. The users can conveniently construct any complex 3D phantom, and then export the phantom as DICOM 3.0 CT images. MPhantom is a user-friendly and powerful software for 3D phantom configuration, and has passed the real scene's application test. MPhantom will accelerate the Monte Carlo simulation for dose calculation in radiation therapy and X ray imaging reconstruction algorithm research.
In the present work, Monte Carlo simulations were employed to study the characteristics of the dose distribution of high energy electron beam in the presence of uniform transverse magnetic field. The simulations carried out the transport processes of the 30 MeV electron beam in the homogeneous water phantom with different magnetic field. It was found that the dose distribution of the 30 MeV electron beam had changed significantly because of the magnetic field. The result showed that the range of the electron beam was decreased obviously and it formed a very high dose peak at the end of the range, and the ratio of maximum dose to the dose of the surface was greatly increased. The results of this study demonstrated that we could change the depth dose distribution of electron beam which is analogous to the heavy ion by modulating the energy of the electron and magnetic field. It means that using magnetic fields in conjunction with electron radiation therapy has great application prospect, but it also has brought new challenges for the research of dose algorithm.
The deoxyribonucleic acid (DNA) molecule damage simulations with an atom level geometric model use the traversal algorithm that has the disadvantages of quite time-consuming, slow convergence and high-performance computer requirement. Therefore, this work presents a density-based spatial clustering of applications with noise (DBSCAN) clustering algorithm based on the spatial distributions of energy depositions and hydroxyl radicals (·OH). The algorithm with probability and statistics can quickly get the DNA strand break yields and help to study the variation pattern of the clustered DNA damage. Firstly, we simulated the transportation of protons and secondary particles through the nucleus, as well as the ionization and excitation of water molecules by using Geant4-DNA that is the Monte Carlo simulation toolkit for radiobiology, and got the distributions of energy depositions and hydroxyl radicals. Then we used the damage probability functions to get the spatial distribution dataset of DNA damage points in a simplified geometric model. The DBSCAN clustering algorithm based on damage points density was used to determine the single-strand break (SSB) yield and double-strand break (DSB) yield. Finally, we analyzed the DNA strand break yield variation trend with particle linear energy transfer (LET) and summarized the variation pattern of damage clusters. The simulation results show that the new algorithm has a faster simulation speed than the traversal algorithm and a good precision result. The simulation results have consistency when compared to other experiments and simulations. This work achieves more precise information on clustered DNA damage induced by proton radiation at the molecular level with high speed, so that it provides an essential and powerful research method for the study of radiation biological damage mechanism.
The Monte Carlo N-Particle (MCNP) is often used to calculate the radiation dose during computed tomography (CT) scans. However, the physical calculation process of the model is complicated, the input file structure of the program is complex, and the three-dimensional (3D) display of the geometric model is not supported, so that the researchers cannot establish an accurate CT radiation system model, which affects the accuracy of the dose calculation results. Aiming at these two problems, this study designed a software that visualized CT modeling and automatically generated input files. In terms of model calculation, the theoretical basis was based on the integration of CT modeling improvement schemes of major researchers. For 3D model visualization, LabVIEW was used as the new development platform, constructive solid geometry (CSG) was used as the algorithm principle, and the introduction of editing of MCNP input files was used to visualize CT geometry modeling. Compared with a CT model established by a recent study, the root mean square error between the results simulated by this visual CT modeling software and the actual measurement was smaller. In conclusion, the proposed CT visualization modeling software can not only help researchers to obtain an accurate CT radiation system model, but also provide a new research idea for the geometric modeling visualization method of MCNP.