Effect of the local energy distribution of x-ray beams generated through inverse Compton scattering in dual-energy imaging applications.

X-ray sources based on the inverse Compton interaction between a laser and a relativistic electron beam are emerging as a promising compact alternative to synchrotron for the production of intense monochromatic and tunable radiation. The emission characteristics enable several innovative imaging techniques, including dual-energy K-edge subtraction (KES) imaging. The performance of these techniques is optimal in the case of perfectly monochromatic x-ray beams, and the implementation of KES was proven to be very effective with synchrotron radiation. Nonetheless, the features of inverse Compton scattering (ICS) sources make them good candidates for a more compact implementation of KES techniques. The energy and intensity distribution of the emitted radiation is related to the emission direction, which means different beam qualities in different spatial positions. In fact, as the polar angle increases, the average energy decreases, while the local energy bandwidth increases and the emission intensity decreases. The scope of this work is to describe the impact of the local energy distribution variations on KES imaging performance. By means of analytical simulations, the reconstructed signal, signal-to-noise ratio, and background contamination were evaluated as a function of the position of each detector pixel. The results show that KES imaging is possible with ICS x-ray beams, even if the image quality slightly degrades at the detector borders for a fixed collimation angle and, in general, as the beam divergence increases. Finally, an approach for the optimization of specific imaging tasks is proposed by considering the characteristics of a given source.

A deep learning classifier for digital breast tomosynthesis.

PURPOSE: To develop a computerized detection system for the automatic classification of the presence/absence of mass lesions in digital breast tomosynthesis (DBT) annotated exams, based on a deep convolutional neural network (DCNN). MATERIALS AND METHODS: Three DCNN architectures working at image-level (DBT slice) were compared: two state-of-the-art pre-trained DCNN architectures (AlexNet and VGG19) customized through transfer learning, and one developed from scratch (DBT-DCNN). To evaluate these DCNN-based architectures we analysed their classification performance on two different datasets provided by two hospital radiology departments.DBT slice images were processed following normalization, background correction and data augmentation procedures. The accuracy, sensitivity, and area-under-the-curve (AUC) values were evaluated on both datasets, using receiver operating characteristic curves. A Grad-CAM technique was also implemented providing anindication of the lesion position in the DBT slice. RESULTS: Accuracy, sensitivity and AUC for the investigated DCNN are in-line with the best performance reported in the field. The DBT-DCNN network developed in this work showed an accuracy and a sensitivity of (90% ± 4%) and (96% ± 3%), respectively, with an AUC as good as 0.89 ± 0.04. Ak-fold cross validation test (withk = 4) showed an accuracy of 94.0% ± 0.2%, and a F1-score test provided a value as good as 0.93 ± 0.03. Grad-CAM maps show high activation in correspondence of pixels within the tumour regions. CONCLUSIONS: We developed a deep learning-based framework (DBT-DCNN) to classify DBT images from clinical exams. We investigated also apossible application of the Grad-CAM technique to identify the lesion position.

BriXS, a new X-ray inverse Compton source for medical applications.

MariX is a research infrastructure conceived for multi-disciplinary studies, based on a cutting-edge system of combined electron accelerators at the forefront of the world-wide scenario of X-ray sources. The generation of X-rays over a large photon energy range will be enabled by two unique X-ray sources: a Free Electron Laser and an inverse Compton source, called BriXS (Bright compact X-ray Source). The X-ray beam provided by BriXS is expected to have an average energy tunable in the range 20-180 keV and intensities between 1011 and 1013 photon/s within a relative bandwidth ΔE/E=1-10%. These characteristics, together with a very small source size (~20 µm) and a good transverse coherence, will enable a wide range of applications in the bio-medical field. An additional unique feature of BriXS will be the possibility to make a quick switch of the X-ray energy between two values for dual-energy and K-edge subtraction imaging. In this paper, the expected characteristics of BriXS will be presented, with a particular focus on the features of interest to its possible medical applications.

GEANT4 Monte Carlo simulations for virtual clinical trials in breast X-ray imaging: Proof of concept.

Virtual clinical trials (VCT) are in-silico reproductions of medical examinations, which adopt digital models of patients and simulated devices. They are intended to produce clinically equivalent outcome data avoiding long execution times, ethical issues related to radiation induced risks and huge costs related to real clinical trials with a patient population. In this work, we present a platform for VCT in 2D and 3D X-ray breast imaging. The VCT platform uses Monte Carlo simulations based on the Geant4 toolkit and patient breast models derived from a cohort of high resolution dedicated breast CT (BCT) volume data sets. Projection images of the breast and three-dimensional glandular dose maps are generated for a given breast model, by simulating both 2D full-field digital mammography (DM) and 3D BCT examinations. Uncompressed voxelized breast models were derived from segmented patient images. Compressed versions of the digital breast phantoms for DM were generated using a previously published digital compression algorithm. The Monte Carlo simulation framework has the capability of generating and tracking ~105 photons/s using a server equipped with 16-cores and 3.0 GHz clock speed. The VCT platform will provide a framework for scanner design optimization, comparison between different scanner designs and between different modalities or protocols on computational breast models, without the need for scanning actual patients as in conventional clinical trials.

Dose intercomparison at Italian hadrontherapy centers.

PURPOSE: To perform the first dosimetric intercomparison for proton beams in Italy using ionization chambers, according to the IAEA TRS-398 code of practice. METHODS: Measurement sites included: National Center for Oncological Hadron Therapy (CNAO, Pavia), Center for Proton Therapy (CTP, Trento) and Center for Hadron Therapy and for advanced Nuclear Applications (CATANA, Catania). For comparison we also included a 6 MV photon beam produced at Istituti Clinici Scientifici Maugeri (ICSM, Pavia). For proton beams, both single pseudo-monoenergetic layers (in order to obtain a planned dose of 2 Gy at the reference depth of 2 cm in a water phantom) and Spread-out Bragg peaks (SOBP) have been delivered. Measurements were performed with a PTW Farmer 30010-1 and a PTW Advanced Markus type 34,045 ionization chamber. RESULTS: Data obtained at CATANA, CNAO and CPT in terms of absorbed dose to water depth show good consistency within the experimental uncertainties, with a weighted mean of 1.99 ± 0.01 Gy and a standard error of 0.003 Gy, with reference to a nominal dose of 2 Gy as designed by the treatment planning system. CONCLUSIONS: The results showed a standard deviation of less than 1% for single layer and SOBP beams, for all chambers and a percent deviation less than 1.5% for single layer measurements. The weighted means of the absorbed doses for clinical proton beams (118.19 MeV and 173.61 MeV) are consistent within less than 1%. These results agree within the 1.5% difference considered acceptable for national dose intercomparison.

Rotational radiotherapy of breast cancer with polyenergetic kilovoltage X-ray beams: An experimental and Monte Carlo phantom study.

PURPOSE: We investigated the feasibility of kilovoltage rotational radiotherapy for breast cancer (kV-EBRT) via Monte Carlo simulations and measurements on phantoms. METHODS: We derived the dose distributions for X-ray beams at 150 kV, 300 kVp and 320 kV irradiating breast cylindrical phantoms of 14 cm diameter, mimicking the pendant breast. Simulations were based on the Geant4 toolkit. The point-like X-ray source was rotated either over a full circle or on a limited arc around the phantom. We studied the influence on the surface dose of the distance between the tumor lesion to the skin, of the irradiation protocol (full scan or partial scan) and of the X-ray tube current modulation. RESULTS: Rotational kV-EBRT permitted a periphery-to-center dose ratio from 13% to 9% in homogeneous breast phantoms. Dose distributions in phantoms with off-center simulated lesions, showed a skin-to-tumor dose ratio of 16% and 34% for lesions at 3.25 and 5.25 cm from cylinder axis, respectively. Simulation of the X-ray tube current modulation during the rotation, permits to reach a dose ratio of 20% for the lesion located at 5.25 cm from phantom axis. CONCLUSIONS: We showed the possibility of using low-energy X-ray spectra for kV-EBRT with collimated beams, for obtaining a periphery-to-center dose ratio in the same order of conventional accelerator based megavoltage radiotherapy, when the irradiated area is localized in the center of the breast. For tumors localized near the breast border, we showed that the tube current modulation can be a good solution in order to reduce the skin-to-tumor dose ratio.

Monte Carlo calculation of monoenergetic and polyenergetic DgN coefficients for mean glandular dose estimates in mammography using a homogeneous breast model.

We computed normalized glandular dose (DgN) coefficients for mean glandular dose estimates in contemporary 2D mammography units, taking into account a homogeneous model for the breast which reflects recent literature reports. We developed a Monte Carlo code based on the simulation toolkit GEANT4 ver. 10.00. The breast was modelled as a cylinder with a semi-cylindrical section with a radius of 10 cm, enveloped in a 1.45 mm thick skin layer, as found out in recent reports in the analysis of breast computed tomography clinical scans. The compressed breast thickness was between 3 cm and 8 cm. The DgN coefficients were calculated for monoenergetic x-ray beams between 4.25 keV and 49.25 keV and were fitted with polynomial curves. Polyenergetic DgN coefficients were then computed for spectra obtained for various anode/filter combinations as adopted in routine clinical practice: Mo/Mo 30 µm (25-40 kV), Mo/Rh 25 µm (25-40 kV), Rh/Rh 25 µm (25-40 kV), W/Ag 50 µm (26-34 kV), W/Al 500 µm (26-38 kV), W/Al 700 µm (28-40 kV) and W/Rh 50 µm (24-35 kV). Monoenergetic DgN curve fit coefficients and polyenergetic DgNp coefficients were released for research and clinical work. Polyenergetic DgNp coefficients were 6% higher than those provided in the recent literature, on average. The differences range between -18% and 30%; up to 50% of the computed coefficients differed by less than 10%. The dataset of DgN coefficients are provided as tables for varying glandular fraction by mass and compressed breast thickness. Moreover, a computer code has been developed for generating user specific coefficients DgNp for user defined x-ray spectra up to 49 kV, calculated by spectral weighting from the dataset of monoenergetic DgN coefficients.

Investigation of the refractive index decrement of 3D printing materials for manufacturing breast phantoms for phase contrast imaging.

3D breast modelling for 2D and 3D breast x-ray imaging would benefit from the availability of digital and physical phantoms that reproduce accurately the complexity of the breast anatomy. While a number of groups have produced digital phantoms with increasing level of complexity, physical phantoms reproducing that software approach have been scarcely developed. One possibility is offered by 3D printing technology. This implies the assessment of the energy dependent absorption index ß of 3D printing materials for absorption based imaging, as well as the assessment of the refractive index decrement, δ, of the printing material, for phase contrast imaging studies, at the energies of interest for breast imaging. In this work we set-up a procedure and performed a series of measurements (at 30, 45 and 60 keV, at the European Synchrotron Radiation Facility) for assessing the relative value of δ with respect to that of breast tissues, for twelve 3D printing materials. The method included propagation based phase contrast 2D imaging and retrieval of the estimated phase shift map, using the Paganin's algorithm. Breast glandular, adipose and skin tissues were used as reference materials of known ratio δ/ß. A percentage difference Δδ was introduced to assess the suitability of the printing materials as tissue substitutes. The accuracy of the method (about 4%) was assessed based on the properties of PMMA and Nylon, acting as gold standard. Results show that, for the above photon energies, ABS is a good substitute for adipose tissue, Hybrid as a substitute of the glandular tissue and PET-G for simulating the skin. We plan to realize a breast phantom manufactured by fused deposition modelling (FDM) technology using ABS, Hybrid and PET-G as substitutes of the glandular and skin tissue and a second phantom by stereolithography (SLA) technology with the resins Flex, Tough and Black.

Normalized glandular dose coefficients in mammography, digital breast tomosynthesis and dedicated breast CT.

PURPOSE: To provide mean glandular dose (MGD) estimates via Monte Carlo (MC) simulations as a function of the breast models and scan parameters in mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (BCT). METHODS: The MC code was based on GEANT4 toolkit. The simulated compressed breast was either a cylinder with a semi-circular section or ad hoc shaped for oblique view (MLO). In DBT we studied the influence of breast models and exam parameters on the T-factors (i.e. the conversion factor for the calculation of the MGD in DBT from that for a 0-degree projection), and in BCT we investigated the influence on the MGD estimates of the ion chamber volume used for the air kerma measurements. RESULTS: In mammography, a model representative of a breast undergoing an MLO view exam did not produce substantial differences (0.4%) in MGD estimates, when compared to a conventional cranio-caudal (CC) view breast model. The beam half value layer did not present a significant influence on T-factors in DBT (<0.8%), while the skin model presented significant influence on MGD estimates (up to 3.3% at 30 degrees scan angle), increasing for larger scan angles. We derived a correction factor for taking into account the different ion chamber volume used in MGD estimates in BCT. CONCLUSIONS: A series of MC code modules for MGD estimates in 2D and 3D breast imaging have been developed in order to take into account the most recent advances in breast models.

Homogeneous vs. patient specific breast models for Monte Carlo evaluation of mean glandular dose in mammography.

PURPOSE: To compare, via Monte Carlo simulations, homogeneous and non-homogenous breast models adopted for mean glandular dose (MGD) estimates in mammography vs. patient specific digital breast phantoms. METHODS: We developed a GEANT4 Monte Carlo code simulating four homogenous cylindrical breast models featured as follows: (1) semi-cylindrical section enveloped in a 5-mm adipose layer; (2) semi-elliptical section with a 4-mm thick skin; (3) semi-cylindrical section with a 1.45-mm skin layer; (4) semi-cylindrical section in a 1.45-mm skin layer and 2-mm subcutaneous adipose layer. Twenty patient specific digital breast phantoms produced from a dedicated CT scanner were assumed as reference in the comparison. We simulated two spectra produced from two anode/filter combinations. An additional digital breast phantom was produced via BreastSimulator software. RESULTS: With reference to the results for patient-specific breast phantoms and for W/Al spectra, models #1 and #3 showed higher MGD values by about 1% (ranges [-33%; +28%] and [-31%; +30%], respectively), while for model #4 it was 2% lower (range [-34%; +26%]) and for model #2 -11% (range [-39%; +14%]), on average. On the other hand, for W/Rh spectra, models #1 and #4 showed lower MGD values by 2% and 1%, while for model #2 and #3 it was 14% and 8% lower, respectively (ranges [-43%; +13%] and [-41%; +21%]). The simulation with the digital breast phantom produced with BreastSimulator showed a MGD overestimation of +33%. CONCLUSIONS: The homogeneous breast models led to maximum MGD underestimation and overestimation of 43% and 28%, respectively, when compared to patient specific breast phantoms derived from clinical CT scans.

Evaluation of the BreastSimulator software platform for breast tomography.

The aim of this work was the evaluation of the software BreastSimulator, a breast x-ray imaging simulation software, as a tool for the creation of 3D uncompressed breast digital models and for the simulation and the optimization of computed tomography (CT) scanners dedicated to the breast. Eight 3D digital breast phantoms were created with glandular fractions in the range 10%-35%. The models are characterised by different sizes and modelled realistic anatomical features. X-ray CT projections were simulated for a dedicated cone-beam CT scanner and reconstructed with the FDK algorithm. X-ray projection images were simulated for 5 mono-energetic (27, 32, 35, 43 and 51 keV) and 3 poly-energetic x-ray spectra typically employed in current CT scanners dedicated to the breast (49, 60, or 80 kVp). Clinical CT images acquired from two different clinical breast CT scanners were used for comparison purposes. The quantitative evaluation included calculation of the power-law exponent, ß, from simulated and real breast tomograms, based on the power spectrum fitted with a function of the spatial frequency, f, of the form S(f) = α/f ß . The breast models were validated by comparison against clinical breast CT and published data. We found that the calculated ß coefficients were close to that of clinical CT data from a dedicated breast CT scanner and reported data in the literature. In evaluating the software package BreastSimulator to generate breast models suitable for use with breast CT imaging, we found that the breast phantoms produced with the software tool can reproduce the anatomical structure of real breasts, as evaluated by calculating the ß exponent from the power spectral analysis of simulated images. As such, this research tool might contribute considerably to the further development, testing and optimisation of breast CT imaging techniques.

In-line phase-contrast breast tomosynthesis: a phantom feasibility study at a synchrotron radiation facility.

The major objective is to adopt, apply and test developed in-house algorithms for volumetric breast reconstructions from projection images, obtained in in-line phase-contrast mode. Four angular sets, each consisting of 17 projection images obtained from four physical phantoms, were acquired at beamline ID17, European Synchroton Radiation Facility, Grenoble, France. The tomosynthesis arc was ±32°. The physical phantoms differed in complexity of texture and introduced features of interest. Three of the used phantoms were in-house developed, and made of epoxy resin, polymethyl-methacrylate and paraffin wax, while the fourth phantom was the CIRS BR3D. The projection images had a pixel size of 47 µm × 47 µm. Tomosynthesis images were reconstructed with standard shift-and-add (SAA) and filtered backprojection (FBP) algorithms. It was found that the edge enhancement observed in planar x-ray images is preserved in tomosynthesis images from both phantoms with homogeneous and highly heterogeneous backgrounds. In case of BR3D, it was found that features not visible in the planar case were well outlined in the tomosynthesis slices. In addition, the edge enhancement index calculated for features of interest was found to be much higher in tomosynthesis images reconstructed with FBP than in planar images and tomosynthesis images reconstructed with SAA. The comparison between images reconstructed by the two reconstruction algorithms shows an advantage for the FBP method in terms of better edge enhancement. Phase-contrast breast tomosynthesis realized in in-line mode benefits the detection of suspicious areas in mammography images by adding the edge enhancement effect to the reconstructed slices.

Imaging performance of phase-contrast breast computed tomography with synchrotron radiation and a CdTe photon-counting detector.

PURPOSE: Within the SYRMA-CT collaboration based at the ELETTRA synchrotron radiation (SR) facility the authors investigated the imaging performance of the phase-contrast computed tomography (CT) system dedicated to monochromatic in vivo 3D imaging of the female breast, for breast cancer diagnosis. METHODS: Test objects were imaged at 38keV using monochromatic SR and a high-resolution CdTe photon-counting detector. Signal and noise performance were evaluated using modulation transfer function (MTF) and noise power spectrum. The analysis was performed on the images obtained with the application of a phase retrieval algorithm as well as on those obtained without phase retrieval. The contrast to noise ratio (CNR) and the capability of detecting test microcalcification clusters and soft masses were investigated. RESULTS: For a voxel size of (60µm)(3), images without phase retrieval showed higher spatial resolution (6.7mm(-1) at 10% MTF) than corresponding images with phase retrieval (2.5mm(-1)). Phase retrieval produced a reduction of the noise level and an increase of the CNR by more than one order of magnitude, compared to raw phase-contrast images. Microcalcifications with a diameter down to 130µm could be detected in both types of images. CONCLUSIONS: The investigation on test objects indicates that breast CT with a monochromatic SR source is technically feasible in terms of spatial resolution, image noise and contrast, for in vivo 3D imaging with a dose comparable to that of two-view mammography. Images obtained with the phase retrieval algorithm showed the best performance in the trade-off between spatial resolution and image noise.

Towards breast tomography with synchrotron radiation at Elettra: first images.

The aim of the SYRMA-CT collaboration is to set-up the first clinical trial of phase-contrast breast CT with synchrotron radiation (SR). In order to combine high image quality and low delivered dose a number of innovative elements are merged: a CdTe single photon counting detector, state-of-the-art CT reconstruction and phase retrieval algorithms. To facilitate an accurate exam optimization, a Monte Carlo model was developed for dose calculation using GEANT4. In this study, high isotropic spatial resolution (120 µm)(3) CT scans of objects with dimensions and attenuation similar to a human breast were acquired, delivering mean glandular doses in the range of those delivered in clinical breast CT (5-25 mGy). Due to the spatial coherence of the SR beam and the long distance between sample and detector, the images contain, not only absorption, but also phase information from the samples. The application of a phase-retrieval procedure increases the contrast-to-noise ratio of the tomographic images, while the contrast remains almost constant. After applying the simultaneous algebraic reconstruction technique to low-dose phase-retrieved data sets (about 5 mGy) with a reduced number of projections, the spatial resolution was found to be equal to filtered back projection utilizing a four fold higher dose, while the contrast-to-noise ratio was reduced by 30%. These first results indicate the feasibility of clinical breast CT with SR.

Energy dependent calibration of XR-QA2 radiochromic film with monochromatic and polychromatic x-ray beams.

PURPOSE: This work investigates the energy response and dose-response curve determinations for XR-QA2 radiochromic film dosimetry system used for synchrotron radiation work and for quality assurance in diagnostic radiology, in the range of effective energies 18-46.5 keV. METHODS: Pieces of XR-QA2 films were irradiated, in a plane transverse to the beam axis, with a monochromatic beam of energy in the range 18-40 keV at the ELETTRA synchrotron radiation facility (Trieste, Italy) and with a polychromatic beam from a laboratory x-ray tube operated at 80, 100, and 120 kV. The film calibration curve was expressed as air kerma (measured free-in-air with an ionization chamber) versus the net optical reflectance change (netΔR) derived from the red channel of the RGB scanned film image. Four functional relationships (rational, linear exponential, power, and logarithm) were tested to evaluate the best curve for fitting the calibration data. The adequacy of the various fitting functions was tested by using the uncertainty analysis and by assessing the average of the absolute air kerma error calculated as the difference between calculated and delivered air kerma. The sensitivity of the film was evaluated as the ratio of the change in net reflectance to the corresponding air kerma. RESULTS: The sensitivity of XR-QA2 films increased in the energy range 18-39 keV, with a maximum variation of about 170%, and decreased in the energy range 38-46.5 keV. The present results confirmed and extended previous findings by this and other groups, as regards the dose response of the radiochromic film XR-QA2 to monochromatic and polychromatic x-ray beams, respectively. CONCLUSIONS: The XR-QA2 radiochromic film response showed a strong dependence on beam energy for both monochromatic and polychromatic beams in the range of half value layer values from 0.55 to 6.1 mm Al and corresponding effective energies from 18 to 46.5 keV. In this range, the film response varied by 170%, from a minimum sensitivity of 0.0127 to a maximum sensitivity of 0.0219 at 10 mGy air kerma in air. The more suitable function for air kerma calibration of the XR-QA2 radiochromic film was the power function. A significant batch-to-batch variation, up to 55%, in film response at 120 kV (46.5 keV effective energy) was observed in comparison with published data.

Energy response of GR-200A thermoluminescence dosemeters to 60Co and to monoenergetic synchrotron radiation in the energy range 28-40 keV.

The response of LiF:Mg,Cu,P thermoluminescence dosemeters (type GR-200A) to monoenergetic radiation of energy 28, 35, 38 and 40 keV was evaluated with respect to irradiation with a calibrated (60)Co gamma-ray source. High-precision measurements of the relative air kerma response performed at the SYRMEP beamline of the ELETTRA synchrotron radiation facility (Trieste, Italy) showed a significant deviation of the average response to low-energy X-rays from that to (60)Co, with an over-response from 6 % (at 28 keV) to 22 % (at 40 keV). These data are not consistent with literature data for these dosemeters, where model predictions gave deviation from unity of the relative air kerma response of about 10 %. The authors conclude for the need of additional determinations of the low-energy relative response of GR-200A dosemeters, covering a wider range of monoenergetic energies sampled at a fine energy step, as planned in future experiments by their group at the ELETTRA facility.

Glandular dose in breast computed tomography with synchrotron radiation.

The purpose of this work is to provide an evaluation of the mean glandular dose (MGD) for breast computed tomography (CT) with synchrotron radiation in an axial scanning configuration with a partial or total organ volume irradiation, for the in vivo program of breast CT ongoing at the ELETTRA facility (Trieste, Italy). A Geant4 Monte Carlo code was implemented, simulating the photon irradiation from a synchrotron radiation source in the energetic range from 8 to 50 keV with 1 keV intervals, to evaluate the MGD. The code was validated with literature data, in terms of mammographic normalized glandular dose coefficients (DgN) and with ad hoc experimental data, in terms of computed tomography dose index (CTDI). Simulated cylindrical phantoms of different sizes (diameter at phantom base 8, 10, 12, 14 or 16 cm, axial length 1.5 times the radius) and glandular fraction by weight (0%, 14.3%, 25%, 50%, 75% and 100%) were implemented into the code. The validation of the code shows an excellent agreement both with previously published work and in terms of DgN and CDTI measurements. The implemented simulations show a dependence of the glandular dose estimate on the vertical dimension of the irradiated zone when a partial organ irradiation was implemented. Specific normalized coefficients for calculating the MGD to the whole breast or to the single irradiated slice were reported.

GEANT4 for breast dosimetry: parameters optimization study.

Mean glandular dose (MGD) is the main dosimetric quantity in mammography. MGD evaluation is obtained by multiplying the entrance skin air kerma (ESAK) by normalized glandular dose (DgN) coefficients. While ESAK is an empirical quantity, DgN coefficients can only be estimated with Monte Carlo (MC) methods. Thus, a MC parameters benchmark is needed for effectively evaluating DgN coefficients. GEANT4 is a MC toolkit suitable for medical purposes that offers to the users several computational choices. In this work we investigate the GEANT4 performances testing the main PhysicsLists for medical applications. Four electromagnetic PhysicsLists were implemented: the linear attenuation coefficients were calculated for breast glandularity 0%, 50%, 100% in the energetic range 8-50 keV and DgN coefficients were evaluated. The results were compared with published data. Fit equations for the estimation of the G-factor parameter, introduced by the literature for converting the dose delivered in the heterogeneous medium to that in the glandular tissue, are proposed and the application of this parameter interaction-by-interaction or retrospectively is discussed. G4EmLivermorePhysicsList shows the best agreement for the linear attenuation coefficients both with theoretical values and published data. Moreover, excellent correlation factor (r2>0.99) is found for the DgN coefficients with the literature. The final goal of this study is to identify, for the first time, a benchmark of parameters that could be useful for future breast dosimetry studies with GEANT4.

A software platform for phase contrast x-ray breast imaging research.

PURPOSE: To present and validate a computer-based simulation platform dedicated for phase contrast x-ray breast imaging research. METHODS: The software platform, developed at the Technical University of Varna on the basis of a previously validated x-ray imaging software simulator, comprises modules for object creation and for x-ray image formation. These modules were updated to take into account the refractive index for phase contrast imaging as well as implementation of the Fresnel-Kirchhoff diffraction theory of the propagating x-ray waves. Projection images are generated in an in-line acquisition geometry. To test and validate the platform, several phantoms differing in their complexity were constructed and imaged at 25 keV and 60 keV at the beamline ID17 of the European Synchrotron Radiation Facility. The software platform was used to design computational phantoms that mimic those used in the experimental study and to generate x-ray images in absorption and phase contrast modes. RESULTS: The visual and quantitative results of the validation process showed an overall good correlation between simulated and experimental images and show the potential of this platform for research in phase contrast x-ray imaging of the breast. The application of the platform is demonstrated in a feasibility study for phase contrast images of complex inhomogeneous and anthropomorphic breast phantoms, compared to x-ray images generated in absorption mode. CONCLUSIONS: The improved visibility of mammographic structures suggests further investigation and optimisation of phase contrast x-ray breast imaging, especially when abnormalities are present. The software platform can be exploited also for educational purposes.