Geant4 Version Research Articles

Monte Carlo simulations of microdosimetry and radiolytic species production at long time post proton irradiation using GATE and Geant4-DNA.

Radiobiological effectiveness of radiation in cancer treatment can be studied at different scales (molecular till organ scale) and different time post irradiation. The production of free radicals and reactive oxygen species during water radiolysis is particularly relevant to understand the fundamental mechanisms playing a role in observed biological outcomes. The development and validation of Monte Carlo tools integrating the simulation of physical, physico-chemical and chemical stages after radiation is very important to maintain with experiments. Therefore, in this study, we propose to validate a new Geant4-DNA chemistry module through the simulation of water radiolysis and Fricke dosimetry experiments on a proton preclinical beam line. In this study, we used the GATE Monte Carlo simulation platform (version 9.3) to simulate a 67.5 MeV proton beam produced with the ARRONAX isochronous cyclotron (IBA Cyclone 70XP) at conventional dose rate (0.2Gy/s) to simulate the irradiation of ultra-pure liquid water samples and Fricke dosimeter. We compared the depth dose profile with measurements performed with a plane parallel Advanced PTW 34045 Markus ionization chamber. Then, a new Geant4-DNA chemistry application proposed from Geant4 version 11.2 has been used to assess the evolution of , , , , , , and reactive species along time until 1-h post-irradiation. In particular, the effect of oxygen and pH has been investigated through comparisons with experimental measurements of radiolytic yields for and Fe3+. GATE simulations reproduced, within 4%, the depth dose profile in liquid water. With Geant4-DNA, we were able to reproduce experimental radiolytic yields 1-h post-irradiation in aerated and deaerated conditions, showing the impact of small changes in oxygen concentrations on species evolution along time. For the Fricke dosimeter, simulated G(Fe3+) is 15.97±0.2 molecules/100eV which is 11% higher than the measured value (14.4±04 molecules/100 eV). These results aim to be consolidated by new comparisons involving other radiolytic species, such as or to further study the mechanisms underlying the FLASH effect observed at ultra-high dose rates (UHDR).

Range Different of Mono-Energetic 290 MeV/u Carbon Ion Across Geant4 Version with Variation of I-values

The mean ionization potential (I-value) is a primary determinant of the position of the heavy ion Bragg peak. To minimize their impact on beam range errors and quantify their uncertainties, the currently used I-values in Geant4 material database are revisited. The study aims at comparing set of I-values in different Geant4 versions and cross validation with PHITS Monte Carlo (MC) code. The Bragg curves of mono-energetic 290 MeV/u carbon ion beams in a Polymethyl Methacrylate (PMMA) phantom were simulated using Geant4 versions 10.6.2 and 11.2.1. Similar beam energies were replicated using the PHITS code. The Bragg curves showed good agreement between the two versions of Geant4 MC code for an I-value of 74 eV, corresponding to G4_PLEXIGLASS in the material database. When the I-value was lowered to 65 eV and 48 eV, the Bragg curves from Geant4 version 11.2.1 shifted to shallower depths. This research provides insights for evaluating the Geant4 physics model.

Review of chemical models and applications in Geant4-DNA: Report from the ESA BioRad III Project.

A chemistry module has been implemented in Geant4-DNA since Geant4 version 10.1 to simulate the radiolysis of water after irradiation. It has been used in a number of applications, including the calculation of G-values and early DNA damage, allowing the comparison with experimental data. Since the first version, numerous modifications have been made to the module to improve the computational efficiency and extend the simulation to homogeneous kinetics in bulk solution. With these new developments, new applications have been proposed and released as Geant4 examples, showing how to use chemical processes and models. This work reviews the models implemented and application developments for modeling water radiolysis in Geant4-DNA as reported in the ESA BioRad III Project.

A quantitative assessment of Geant4 for predicting the yield and distribution of positron-emitting fragments in ion beam therapy

Objective. To compare the accuracy with which different hadronic inelastic physics models across ten Geant4 Monte Carlo simulation toolkit versions can predict positron-emitting fragments produced along the beam path during carbon and oxygen ion therapy. Approach. Phantoms of polyethylene, gelatin, or poly(methyl methacrylate) were irradiated with monoenergetic carbon and oxygen ion beams. Post-irradiation, 4D PET images were acquired and parent 11C, 10C and 15O radionuclides contributions in each voxel were determined from the extracted time activity curves. Next, the experimental configurations were simulated in Geant4 Monte Carlo versions 10.0 to 11.1, with three different fragmentation models—binary ion cascade (BIC), quantum molecular dynamics (QMD) and the Liege intranuclear cascade (INCL++) - 30 model-version combinations. Total positron annihilation and parent isotope production yields predicted by each simulation were compared between simulations and experiments using normalised mean squared error and Pearson cross-correlation coefficient. Finally, we compared the depth of the maximum positron annihilation yield and the distal point at which the positron yield decreases to 50% of peak between each model and the experimental results. Main results. Performance varied considerably across versions and models, with no one version/model combination providing the best prediction of all positron-emitting fragments in all evaluated target materials and irradiation conditions. BIC in Geant4 10.2 provided the best overall agreement with experimental results in the largest number of test cases. QMD consistently provided the best estimates of both the depth of peak positron yield (10.4 and 10.6) and the distal 50%-of-peak point (10.2), while BIC also performed well and INCL generally performed the worst across most Geant4 versions. Significance. The best predictions of the spatial distribution of positron annihilations and positron-emitting fragment production along the beam path during carbon and oxygen ion therapy was obtained using Geant4 10.2.p03 with BIC or QMD. These version/model combinations are recommended for future heavy ion therapy research.

Fano cavity test and investigation of the response of the Roos chamber irradiated by proton beams in perpendicular magnetic fields up to 1 T

Objective. The aim of this work is to investigate the response of the Roos chamber (type 34001) irradiated by clinical proton beams in magnetic fields. Approach. At first, a Fano test was implemented in Monte Carlo software package GATE version 9.2 (based on Geant4 version 11.0.2) using a cylindrical slab geometry in a magnetic field up to 1 T. In accordance to an experimental setup (Fuchs et al ), the magnetic field correction factors kQB⃗ of the Roos chamber were determined at different energies up to 252 MeV and magnetic field strengths up to 1 T, by separately simulating the ratios of chamber signals MQ/MQB⃗, without and with magnetic field, and the dose-conversion factors Dw,QB⃗/Dw,Q in a small cylinder of water, with and without magnetic field. Additionally, detailed simulations were carried out to understand the observed magnetic field dependence. Main results. The Fano test was passed with deviations smaller than 0.25% between 0 and 1 T. The ratios of the chamber signals show both energy and magnetic field dependence. The maximum deviation of the dose-conversion factors from unity of 0.22% was observed at the lowest investigated proton energy of 97.4 MeV and B⃗ = 1 T. The resulting kQB⃗ factors increase initially with the applied magnetic field and decrease again after reaching a maximum at around 0.5 T; except for the lowest 97.4 MeV beam that show no observable magnetic field dependence. The deviation from unity of the factors is also larger for higher proton energies, where the maximum lies at 1.0035(5), 1.0054(7) and 1.0069(7) for initial energies of E0 = 152, 223.4 and 252 MeV, respectively. Significance. Detailed Monte Carlo studies showed that the observed effect can be mainly attributed to the differences in the transport of electrons produced both outside and inside of the air cavity in the presence of a magnetic field.

The influence of different versions of FLUKA and GEANT4 on the calculation of response functions of ionization chambers in clinical proton beams

Objective. To investigate the influence of different versions of the Monte Carlo codes geant4 and fluka on the calculation of overall response functions f Q of air-filled ionization chambers in clinical proton beams. Approach. f Q factors were calculated for six plane-parallel and four cylindrical ionization chambers with geant4 and fluka. These factors were compared to already published values that were derived using older versions of these codes. Main results. Differences in f Q factors calculated with different versions of the same Monte Carlo code can be up to ∼1%. Especially for geant4, the updated version leads to a more pronounced dependence of f Q on proton energy and to smaller f Q factors for high energies. Significance. Different versions of the same Monte Carlo code can lead to differences in the calculation of f Q factors of up to ∼1% without changing the simulation setup, transport parameters, ionization chamber geometry modeling, or employed physics lists. These findings support the statement that the dominant contributor to the overall uncertainty of Monte Carlo calculated f Q factors are type-B uncertainties.

Evaluation of the neutron spectrum at flight altitude with Geant4 using different parameterizations.

Cosmic-ray-induced particles at flight altitude present an increasing concern because of the development of aircraft with high maximum cruising altitude (higher than 40000ft) and because of the miniaturisation of and increase in the amount of embedded electronic systems. There are different codes used to estimate the radiation field through the atmosphere, using different approaches and Monte Carlo methods. This work aims to evaluate the influence on the neutron spectra at ground and flight altitude of the parameterizations employed to describe spallation reactions and the density profile of the atmosphere. The results show that different versions of Geant4 and the data driven models used significantly change the neutron fluence rate and the ambient dose equivalent rate. The results show that the different atmospheric density profiles considered do not significantly change the neutron spectra. In this work, we also present comparisons with onboard and ground level measurements.

Development of a stand-alone precalculated Monte Carlo code to calculate the dose by alpha and beta emitters from the Ra-224 decay chain.

Recent developments in alpha and beta emitting radionuclide therapy highlight the importance of developing efficient methods for patient-specific dosimetry. Traditional tabulated methods such as Medical Internal Radiation Dose (MIRD) estimate the dose at the organ level while more recent numerical methods based on Monte Carlo (MC) simulations are able to calculate dose at the voxel level. A precalculated MC (PMC) approach was developed in this work as an alternative to time-consuming fully simulated MC. Once the spatial distribution of alpha and beta emitters is determined using imaging and/or numerical methods, the PMC code can be used to achieve an accurate voxelized 3D distribution of the deposited energy without relying on full MC calculations. To implement the PMC method to calculate energy deposited by alpha and beta particles emitted from the Ra-224 decay chain. The GEANT4 (version 10.7) MC toolkit was used to generate databases of precalculated tracks to be integrated in the PMC code as well as to benchmark its output. In this regard, energy spectra of alpha and beta particles emitted by the Ra-224 decay chain were generated using GAMOS (version 6.2.0) and imported into GEANT4 macro files. Either alpha or beta emitting sources were defined at the center of a homogeneous phantom filled with various materials such as soft tissue, bone, and lung where particles were emitted either mono-directionally (for database generation) or isotropically (for benchmarking). Two heterogeneous phantoms were used to demonstrate PMC code compatibility with boundary crossing events. Each precalculated database was generated step-by-step by storing particle track information from GEANT4 simulations followed by its integration in a PMC code developed in MATLAB. For a user-defined number of histories, one of the tracks in a given database was selected randomly and rotated randomly to reflect an isotropic emission. Afterward, deposited energy was divided between voxels based on step length in each voxel using a ray-tracing approach. The radial distribution of deposited energy was benchmarked against fully simulated MC calculations using GEANT4. The effect of the GEANT4 parameter StepMax on the accuracy and speed of the code was also investigated. In the case of alpha decay, primary alpha particles show the highest contribution (>99%) in deposited energy compared to their secondary particles. In most cases, protons act as the main secondary particles in the deposition of energy. However, for a lung phantom, using a range cutoff parameter of 10µm on primary alpha particles yields a higher contribution of secondary electrons than protons. Differences between deposited energy calculated by PMC and fully simulated MC are within 2% for all alpha and beta emitters in homogeneous and heterogeneous phantoms. Additionally, statistical uncertainties are less than 1% for voxels with doses higher than 5% of the maximum dose. Moreover, optimization of the parameter StepMax is necessary to achieve the best tradeoff between code accuracy and speed. The PMC code shows good performance for dose calculations deposited by alpha and beta emitters. As a stand-alone algorithm, it is suitable to be integrated into clinical treatment planning systems.

Simulation of DNA damage using Geant4‐DNA: an overview of the “molecularDNA” example application

AbstractPurposeThe scientific community shows great interest in the study of DNA damage induction, DNA damage repair, and the biological effects on cells and cellular systems after exposure to ionizing radiation. Several in silico methods have been proposed so far to study these mechanisms using Monte Carlo simulations. This study outlines a Geant4‐DNA example application, named “molecularDNA”, publicly released in the 11.1 version of Geant4 (December 2022).MethodsIt was developed for novice Geant4 users and requires only a basic understanding of scripting languages to get started. The example includes two different DNA‐scale geometries of biological targets, namely “cylinders” and “human cell”. This public version is based on a previous prototype and includes new features, such as: the adoption of a new approach for the modeling of the chemical stage, the use of the standard DNA damage format to describe radiation‐induced DNA damage, and upgraded computational tools to estimate DNA damage response.ResultsSimulation data in terms of single‐strand break and double‐strand break yields were produced using each of the available geometries. The results were compared with the literature, to validate the example, producing less than 5% difference in all cases. Conclusion: “molecularDNA” is a prototype tool that can be applied in a wide variety of radiobiology studies, providing the scientific community with an open‐access base for DNA damage quantification calculations. New DNA and cell geometries for the “molecularDNA” example will be included in future versions of Geant4‐DNA.

Multiple Coulomb scattering of muons in lithium hydride

Multiple Coulomb Scattering (MCS) is a well known phenomenon occurring when charged particles traverse materials. Measurements of muons traversing low $Z$ materials made in the MuScat experiment showed that theoretical models and simulation codes, such as GEANT4 (v7.0), over-estimated the scattering. The Muon Ionization Cooling Experiment (MICE) measured the cooling of a muon beam traversing a liquid hydrogen or lithium hydride (LiH) energy absorber as part of a programme to develop muon accelerator facilities, such as a Neutrino Factory or a Muon Collider. The energy loss and MCS that occur in the absorber material are competing effects that alter the performance of the cooling channel. Therefore measurements of MCS are required in order to validate the simulations used to predict the cooling performance in future accelerator facilities. We report measurements made in the MICE apparatus of MCS using a LiH absorber and muons within the momentum range 160 to 245 MeV/c. The measured RMS scattering width is about 9% smaller than that predicted by the approximate formula proposed by the Particle Data Group. Data at 172, 200 and 240 MeV/c are compared to the GEANT4 (v9.6) default scattering model. These measurements show agreement with this more recent GEANT4 (v9.6) version over the range of incident muon momenta.

Experimental validation of proton physics models of Geant4 for calculating stopping power ratio

In this work, we conducted experiments to validate the proton physics models of Geant4 (version 10.6). The stopping power ratios (SPRs) of 11 inserts, such as acrylic, delrin, high density polyethylene, and polytetrafluoroethylene, etc, were measured using a superconducting synchrocyclotron that produces a scattering proton beam. The SPRs of the inserts were also calculated based on Geant4 simulation with six physics lists, i.e. QGSP_ FTFP_ BERT, QGSP_BIC_HP, QGSP_BIC, QGSP_FTFP_BERT, QSGP_BERT, and QBBC. The calculated SPRs were compared to the experimental SPRs, and relative per cent error was used to quantify the accuracy of the simulated SPRs of inserts. The comparison showed that the five physics lists generally agree well with the experimental SPRs with a relative difference of less than 1%. The lowest overall percentage error was observed for QGSP_FTFP_BERT and the highest overall percentage error was observed for QGSP_BIC_HP. The 0.1 mm range cut value consistently led to higher percentage error for all physics lists except for QGSP_BIC_HP and QBBC. Based on the validation, we recommend QGSP_BERT_HP physics list for proton dose calculation.

Experimental and Geant4 Simulation Study of MURA Mask for Scintimammography

Gamma rays have been extensively investigated for breast imaging using collimators; however, the coded-aperture technique needs to be investigated more. In this paper, we propose an experimental study and Geant4 simulations of MURA mask breast imaging. First, we compare the experimental data against the simulation results carried out using Geant4 (version 10.4) and accreditation phantom. Second, we virtually extend our work by changing the tumor-to-background (TBR) and lesion location parameters. We used 99mTc as a radioactive source. Good agreement has been seen for the benchmark stage, especially in terms of tumor localization. Moreover, the calculated full width at half maximum (FWHM) and contrast for decoded images (having average values of 8 and 3.5 for TBR between 2 and 10) permitted us to conclude that we can accurately localize small lesions up to lower TBR values by following the decoding procedure of deducing the image of a “blank phantom” (phantom with TBR = 1) every time within a matlab-based program. Hence, this work can be considered a continuously added value to previous investigations for scintimammography imaging.

Simulation of Gamma-Ray Transmission Buildup Factors for Stratified Spherical Layers.

Deterministic particle transport codes usually take into account scattered photons with correct attenuation laws and application of buildup factor to incident beam. Transmission buildup factors for adipose, bone, muscle, and skin human tissues, as well as for various combinations of these media for point isotropic photon source with energies of .15, 1.5 and 15 MeV, for different thickness of layers, were carried out using Geant4 (version 10.5) simulation toolkit. Also, we performed the analysis of existing multilayered shield fitting models (Lin and Jiang, Kalos, Burke and Beck) of buildup factor and the proposition of a new model. We found that the model combining those of Burke and Beck, for low atomic number (Z) followed by high Z materials and Kalos 1 for high Z followed by low Z materials, accurately reproduces simulation results with approximated deviation of 3 ± 3%, 2 ± 2%, and 3 ± 2% for 2, 3, and 4 layers, respectively. Since buildup factors are the key parameter for point kernel calculations, a correct study can be of great interest to the large community of radiation physicists, in general, and to medical imaging and radiotreatment physicists, especially.

Geant4 Simulation of Gamma-Ray Transmission Buildup Factors for Human Tissues.

Quantification of scattered photons in addition to unscattered primary particles, under realistic exposure scenario, is best dealt with a parameter called “Buildup factor”. The aim of this work is to simulate the transmission buildup factor (BUF) of gamma-ray in the energy range .15–15 MeV for 20 human tissues and organs using the Geant4 (version 10.5) Monte Carlo simulation followed by a geometrical progression (GP) parameterization procedure. Firstly, we verified the accuracy of Geant4 ability to predict the effective transmitted dose according to published data. Also, a comparison of simulated BUF for different geometrical configurations was carried out for some tissues and source energies. Then, we checked out the linear dependency of the K parameter (BUF is function of K) function of mean free path (mfp). Finally, we developed a fitting procedure according to GP method for BUF corresponding to 20 tissues and organs and different mfp (from 1 to 8) for energy range .15–15 MeV. We found a good agreement with previous published data. Proper comprehension of BUF for tissues leads to carefully controlling the energy absorption in the human body. Consequently, provided BUF could be of great interest for estimating safe dose levels in medical imaging and radiation therapy.

Atomic Compton scattering effect on direct dark matter detection

Atomic Compton scattering effect significantly contributes to low-energy electronic recoils below its k-shell energy for the direct dark matter detection. Searches on ADM models, dark photon models, leptophilic dark matter models as well as the conventional WIMPs for background understandings are vitally required to clarify the effect. We employed the relativistic impulse approximation (RIA) together with the ab initio Multi-Configuration Dirac-Fock (MCDF) theory to obtain the atomic Compton scattering for Germanium (Ge) Silicon (Si) and Xenon (Xe) atoms. Comparisons on low momentum transfer regions with our calculations for Ge and Si are achieved. In addition, millicharged dark matter particles estimated by RIA in the atomic ionization for Ge and Xe have been evaluated. A factor-of-two discrepancy on the incoherent-scattering factor (a.k.a. scattering function) near 100 eV/c momentum transfer with the Ge system between our calculation and the latest version of Geant4 (10.07.02) simulation data is observed. Plans on the experimental verification and the perspectives of the atomic Compton scattering effect for the direct detections will be discussed.

Prompt-gamma emission in GEANT4 revisited and confronted with experiment.

The field of online monitoring of the beam range is one of the most researched topics in proton therapy over the last decade. The development of detectors that can be used for beam range verification under clinical conditions is a challenging task. One promising possible solution are modalities that record prompt-gamma radiation produced by the interactions of the proton beam with the target tissue. A good understanding of the energy spectra of the prompt gammas and the yields in certain energy regions is crucial for a successful design of a prompt-gamma detector. Monte-Carlo simulations are an important tool in development and testing of detector concepts, thus the proper modelling of the prompt-gamma emission in those simulations are of vital importance. In this paper, we confront a number of GEANT4 simulations of prompt-gamma emission, performed with different versions of the package and different physics lists, with experimental data obtained from a phantom irradiation with proton beams of four different energies in the range 70-230 MeV. The comparison is made on different levels: features of the prompt-gamma energy spectrum, gamma emission depth profiles for discrete transitions and the width of the distal fall-off in those profiles. The best agreement between the measurements and the simulations is found for the GEANT4 version 10.4.2 and the reference physics list QGSP_BIC_HP. Modifications to prompt-gamma emission modelling in higher versions of the software increase the discrepancy between the simulation results and the experimental data.

Monte Carlo study of patient and medical staff radiation exposures during interventional cardiology.

The aim of this study is to assess the radiation exposure of the patient and the medical staff during interventional cardiology procedures. Realistic exposure scenarios were developed using the adult reference anthropomorphic phantoms adopted by the International Commission on Radiological Protection (ICRP110Male and ICRP110Female), and the radiation transport code Geant4 (version 10.3). The calculated equivalent and effective doses were normalised by the simulated Kerma-Area Product (KAP), resulting in two conversion coefficients HT/KAP and E/KAP. To properly evaluate the risk of exposure, several dose-dependent parameters have been investigated, namely: radiological parameters (tube kilovoltage peak (kVp), type of projection, field size (FOV)), and operator positions. Four projections (AP,PA,LAO25° and RAO25°) were simulated for three X-ray energy spectra (80,100 and 120 kVp) with four different values of FOV (15×15cm2,20×20cm2,25×25cm2 and 30×30cm2). The results showed that the conversion coefficients values increase with increasing tube voltage as well as the FOV size. Recommended projection during the interventional cardiology procedures, whenever possible, should be the PA projection rather than AP projection. The most critical projection for the patient and the main operator is the RAO25° projection and the LAO25° projection respectively. The comparison of our results with the literature data showed good agreement allowing their use in the dosimetric characterization of interventional cardiology procedures.

Validation of Physics Models of Geant4 Versions 10.4.p03, 10.6.p02 and 10.7.p01 using Data from the CMS Experiment

CMS tuned its simulation program and chose a specific physics model of Geant4 by comparing the simulation results with dedicated test beam experiments. Test beam data provide measurements of energy response of the calorimeter as well as resolution for well identified charged hadrons over a large energy region. CMS continues to validate the physics models using the test beam data as well as collision data from the Large Hadron Collider. Isolated charged particles are measured simultaneously in the tracker as well as in the calorimeters. These events are selected using dedicated triggers and are used to measure the response in the calorimeter. Different versions of Geant4 (10.2.p02, 10.4.p03, 10.6.p02) have been used by CMS for its Monte Carlo production and a new version (10.7.p01) is now chosen for future productions. A suitable physics list (collection of physics models) is chosen by optimizing performance against accuracy. A detailed comparison between data and Geant4 predictions is presented in this paper.

CMS Full Simulation for Run 3

We report the status of the CMS full simulation for Run 3. During the long shutdown of the LHC a significant update has been introduced to the CMS code for simulation. The CMS geometry description is reviewed. Several important modifications were needed. CMS detector description software is migrated to the DD4Hep community developed tool. We will report on our experience obtained during the process of this migration. Geant4 10.7 is the CMS choice for Run 3 simulation productions. We will discuss arguments for this choice, the strategy of adaptation of a new Geant4 version, and will report on the physics performance of the CMS simulation. A special Geant4 Physics List configuration FTFP_BERT_EMM will be described, which provides a compromise between simulation accuracy and CPU performance. A significant fraction of time for simulation of CMS events is spent on tracking of charged particles in a magnetic field. In the CMS simulation a dynamic choice of Geant4 parameters for tracking in field is implemented. A new method is introduced into simulation of electromagnetic components of hadronic showers in the electromagnetic calorimeter of CMS. For low-energy electrons and positrons a parametrization of GFlash type is applied. Results of tests of this method will be discussed. In summary, we expect about 25% speedup of the CMS simulation production for Run 3 compared to the Run 2 simulations.

Pulse shape simulations for organic scintillation detectors using Geant4

The accurate simulation of the temporal pulse shapes from organic scintillation detectors capable of pulse shape discrimination (PSD) presents the opportunity to assess the pulse shape discrimination of these detectors prior to fabrication. The aim of this study is the simulation of the temporal pulse shapes from EJ-276, a PSD-capable plastic scintillator developed by Eljen Technologies. PSD plastic scintillators are increasingly replacing organic liquid scintillators for the detection of neutrons in the presence of mixed radiation fields for nuclear security applications. Plastics are inexpensive, robust and can be fabricated in a variety of shapes and sizes. They offer a solid-state alternative to liquid scintillators which can be difficult to transport due to the risk of leakage. However, the PSD performance of plastic scintillators has been observed to decrease due to various factors which combine to influence the overall shape of the pulse. The Monte Carlo toolkit Geant4 has been used to simulate the temporal pulse shapes from an EJ-276 plastic scintillator coupled to a photomultiplier tube (PMT). All three decay time components of EJ-276 have been modelled, utilising new methods available in the latest version of Geant4, for two different scintillator geometries. The simulated n/γ pulse shapes reproduce the features important for PSD. Future work will involve integrating the temporal response of the PMT with existing pulse shape simulations. Simulated data will then be compared with experimental measurements.