Monte Carlo studies on neutron interactions in radiobiological experiments.
Ontology highlight
ABSTRACT:
SUBMITTER: Shahmohammadi Beni M
PROVIDER: S-EPMC5509315 | biostudies-other | 2017
REPOSITORIES: biostudies-other
Similar Datasets
Project description:The purpose of this study is to show how Monte Carlo analysis of meta-analytic estimators can be used to select estimators for specific research situations. Our analysis conducts 1620 individual experiments, where each experiment is defined by a unique combination of sample size, effect size, effect size heterogeneity, publication selection mechanism, and other research characteristics. We compare 11 estimators commonly used in medicine, psychology, and the social sciences. These are evaluated on the basis of bias, mean squared error (MSE), and coverage rates. For our experimental design, we reproduce simulation environments from four recent studies. We demonstrate that relative estimator performance differs across performance measures. Estimator performance is a complex interaction of performance indicator and aspects of the application. An estimator that may be especially good with respect to MSE may perform relatively poorly with respect to coverage rates. We also show that the size of the meta-analyst's sample and effect heterogeneity are important determinants of relative estimator performance. We use these results to demonstrate how these observable characteristics can guide the meta-analyst to choose the most appropriate estimator for their research circumstances.
| S-EPMC8074967 | biostudies-literature
Project description:In this study, we present a validated Geant4 Monte Carlo simulation model of the Dingo thermal neutron imaging beamline at the Australian Centre for Neutron Scattering. The model, constructed using CAD drawings of the entire beam transport path and shielding structures, is designed to precisely predict the in-beam neutron field at the position at the sample irradiation stage. The model's performance was assessed by comparing simulation results to various experimental measurements, including planar thermal neutron distribution obtained in-beam using gold foil activation and [Formula: see text]B[Formula: see text]C-coated microdosimeters and the out-of-beam neutron spectra measured with Bonner spheres. The simulation results demonstrated that the predicted neutron fluence at the field's centre is within 8.1% and 2.1% of the gold foil and [Formula: see text]B[Formula: see text]C-coated microdosimeter measurements, respectively. The logarithms of the ratios of average simulated to experimental fluences in the thermal (E[Formula: see text] 0.414 eV), epithermal (0.414 eV < E[Formula: see text] 11.7 keV) and fast (E[Formula: see text] 11.7 keV) spectral regions were approximately - 0.03 to + 0.1, - 0.2 to + 0.15, and - 0.4 to + 0.2, respectively. Furthermore, the predicted thermal, epithermal and fast neutron components in-beam at the sample stage position constituted approximately 18%, 64% and 18% of the total neutron fluence.
| S-EPMC10575880 | biostudies-literature
Project description:Increasingly, microbeams and microcrystals are being used for macromolecular crystallography (MX) experiments at synchrotrons. However, radiation damage remains a major concern since it is a fundamental limiting factor affecting the success of macromolecular structure determination. The rate of radiation damage at cryotemperatures is known to be proportional to the absorbed dose, so to optimize experimental outcomes, accurate dose calculations are required which take into account the physics of the interactions of the crystal constituents. The program RADDOSE-3D estimates the dose absorbed by samples during MX data collection at synchrotron sources, allowing direct comparison of radiation damage between experiments carried out with different samples and beam parameters. This has aided the study of MX radiation damage and enabled prediction of approximately when it will manifest in diffraction patterns so it can potentially be avoided. However, the probability of photoelectron escape from the sample and entry from the surrounding material has not previously been included in RADDOSE-3D, leading to potentially inaccurate does estimates for experiments using microbeams or microcrystals. We present an extension to RADDOSE-3D which performs Monte Carlo simulations of a rotating crystal during MX data collection, taking into account the redistribution of photoelectrons produced both in the sample and the material surrounding the crystal. As well as providing more accurate dose estimates, the Monte Carlo simulations highlight the importance of the size and composition of the surrounding material on the dose and thus the rate of radiation damage to the sample. Minimizing irradiation of the surrounding material or removing it almost completely will be key to extending the lifetime of microcrystals and enhancing the potential benefits of using higher incident X-ray energies.
| S-EPMC7737758 | biostudies-literature
Project description:Carbon-13 (13C) analysis is a commonly used method for estimating reaction rates in biochemical networks. The choice of carbon labeling pattern is an important consideration when designing these experiments. We present a novel Monte Carlo algorithm for finding the optimal substrate input label for a particular experimental objective (flux or flux ratio). Unlike previous work, this method does not require assumption of the flux distribution beforehand.Using a large E. coli isotopomer model, different commercially available substrate labeling patterns were tested computationally for their ability to determine reaction fluxes. The choice of optimal labeled substrate was found to be dependent upon the desired experimental objective. Many commercially available labels are predicted to be outperformed by complex labeling patterns. Based on Monte Carlo Sampling, the dimensionality of experimental data was found to be considerably less than anticipated, suggesting that effectiveness of 13C experiments for determining reaction fluxes across a large-scale metabolic network is less than previously believed.While 13C analysis is a useful tool in systems biology, high redundancy in measurements limits the information that can be obtained from each experiment. It is however possible to compute potential limitations before an experiment is run and predict whether, and to what degree, the rate of each reaction can be resolved.
| S-EPMC3323462 | biostudies-literature
Project description:The PGNAA system for the cement measurement was simulated based on Monte Carlo method. The sizes of the moderator and reflector for the 14 MeV DT neutron generator were optimized for fast and thermal neutron outputs. The DT neutron generator was simulated at the pulse mode, and the gamma-ray detector was set as LaBr3(Ce) scintillator. The characteristic peaks of the major elements (Ca, Si, Al, Fe) can be identified from the gamma-ray spectra which induced at the different time intervals of the neutron radiation. For the different thicknesses of the cement sample the ratios of the gamma-ray peaks were observed, and the result showed that when the thickness was between 20 to 30 cm, the ratios became stable. With the ratios, we can calculate the iron modulus, silica modulus and lime saturation factor.
| S-EPMC8202931 | biostudies-literature
Project description:The development of neutron imaging from a qualitative inspection tool towards a quantitative technique in materials science has increased the requirements for accuracy significantly. Quantifying the thickness or the density of polycrystalline samples with high accuracy using neutron imaging has two main problems: (i) the scattering from the sample creates artefacts on the image and (ii) there is a lack of specific reference attenuation coefficients. This work presents experimental and simulation results to explain and approach these problems. Firstly, a series of neutron radiography and tomography experiments of iron, copper and vanadium are performed and serve as a reference. These materials were selected because they attenuate neutrons mainly through coherent (Fe and Cu) and incoherent (V) scattering. Secondly, an ad hoc Monte Carlo model was developed, based on beamline, sample and detector parameters, in order to simulate experiments, understand the physics involved and interpret the experimental data. The model, developed in the McStas framework, uses a priori information about the sample geometry and crystalline structure, as well as beamline settings, such as spectrum, geometry and detector type. The validity of the simulations is then verified with experimental results for the two problems that motivated this work: (i) the scattering distribution in transmission imaging and (ii) the calculated attenuation coefficients.
| S-EPMC5884388 | biostudies-other
Project description:Biomacromolecule activity is usually related to its ability to keep a specific structure. However, in solution, many parameters (pH, ionic strength) and external compounds (polyelectrolytes, nanoparticles) can modify biomacromolecule structure as well as acid/base properties, thus resulting in a loss of activity and denaturation. In this paper, the impact of neutral and charged nanoparticles (NPs) is investigated by Monte Carlo simulations on polypeptide (PP) chains with primary structure based on bovine serum albumin. The influence of pH, salt valency, and NP surface charge density is systematically studied. It is found that the PP is extended at extreme pH, when no complex formation is observed, and folded at physiological pH. PP adsorption around oppositely-charged NPs strongly limits chain structural changes and modifies its acid/base properties. At physiological pH, the complex formation occurs only with positively-charged NPs. The presence of salts, in particular those with trivalent cations, introduces additional electrostatic interactions, resulting in a mitigation of the impact of negative NPs. Thus, the corona structure is less dense with locally-desorbed segments. On the contrary, very limited impact of salt cation valency is observed when NPs are positive, due to the absence of competitive effects between multivalent cations and NP.
| S-EPMC6432162 | biostudies-literature
Project description:Background: There is currently no reliable or validated tool to delineate and quantify functional lung volumes with ventilation/perfusion (V/Q) SPECT/CT. The main challenges encountered include the physiological non-uniformity of lung function, such as the anterior-to-posterior gradient on perfusion images, and the lack of ground truth to assess the accuracy of delineation algorithms. In that respect, Monte-Carlo simulations would be an interesting tool. Thus, the aim of this study was to develop a realistic model of dual-isotope lung V/Q SPECT-CT Monte-Carlo simulations, integrating the anterior to posterior gradient on perfusion. Methods: Acquisitions and simulations parameters were set in accordance to nuclear medicine guidelines for V/Q lung SPECT-CT. Projections were acquired and simulated, then the reconstructions [with and without attenuation correction (AC)] were compared. A model was built from a patient's CT scan. To model the anterior to posterior gradient, the lungs were divided into sixteen coronal planes, where a rising radioactivity concentration was set. To assess the realism of simulations, they were compared to a normal co-registered normal cases database in terms of pixelwize Z-score map. Results: For ventilation images, mean (SD) Zscores on Zscore maps were -0.2 (0.7) and -0.2 (0.7) for AC and noAC images, respectively. For perfusion images, mean (SD) Zscores were -0.2 (0.6) and -0.1 (0.6) for AC and noAC images, respectively. Conclusion: We developed a model for dual isotopes lung V/Q SPECT-CT, integrating the anterior-to-posterior gradient on perfusion images. This model could be used to build a catalog of clinical scenarios, in order to test delineation methods of functional lung volumes.
| S-EPMC7461964 | biostudies-literature
Project description:The use of Monte-Carlo (MC) p$$ p $$ -values when testing the significance of a large number of hypotheses is now commonplace. In large-scale hypothesis testing, we will typically encounter at least some p$$ p $$ -values near the threshold of significance, which require a larger number of MC replicates than p$$ p $$ -values that are far from the threshold. As a result, some incorrect conclusions can be reached due to MC error alone; for hypotheses near the threshold, even a very large number (eg, 106$$ 1{0}^6 $$ ) of MC replicates may not be enough to guarantee conclusions reached using MC p$$ p $$ -values. Gandy and Hahn (GH)6-8 have developed the only method that directly addresses this problem. They defined a Monte-Carlo error rate (MCER) to be the probability that any decisions on accepting or rejecting a hypothesis based on MC p$$ p $$ -values are different from decisions based on ideal p$$ p $$ -values; their method then makes decisions by controlling the MCER. Unfortunately, the GH method is frequently very conservative, often making no rejections at all and leaving a large number of hypotheses "undecided". In this article, we propose MERIT, a method for large-scale MC hypothesis testing that also controls the MCER but is more statistically efficient than the GH method. Through extensive simulation studies, we demonstrate that MERIT controls the MCER while making more decisions that agree with the ideal p$$ p $$ -values than GH does. We also illustrate our method by an analysis of gene expression data from a prostate cancer study.
| S-EPMC10909586 | biostudies-literature
Project description:The MCPep server (http://bental.tau.ac.il/MCPep/) is designed for non-experts wishing to perform Monte Carlo (MC) simulations of helical peptides in association with lipid membranes. MCPep is a web implementation of a previously developed MC simulation model. The model has been tested on a variety of peptides and protein fragments. The simulations successfully reproduced available empirical data and provided new molecular insights, such as the preferred locations of peptides in the membrane and the contribution of individual amino acids to membrane association. MCPep simulates the peptide in the aqueous phase and membrane environments, both described implicitly. In the former, the peptide is subjected solely to internal conformational changes, and in the latter, each MC cycle includes additional external rigid body rotational and translational motions to allow the peptide to change its location in the membrane. The server can explore the interaction of helical peptides of any amino-acid composition with membranes of various lipid compositions. Given the peptide's sequence or structure and the natural width and surface charge of the membrane, MCPep reports the main determinants of peptide-membrane interactions, e.g. average location and orientation in the membrane, free energy of membrane association and the peptide's helical content. Snapshots of example simulations are also provided.
| S-EPMC3394254 | biostudies-literature