Monte Carlo Solution of High Electric Field Hole Transport Processes in Avalanche Amorphous Selenium.
Ontology highlight
ABSTRACT:
SUBMITTER: Mukherjee A
PROVIDER: S-EPMC7905821 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In this work, we present a robust and powerful method for the verification, with arbitrary accuracy, of Monte Carlo codes for simulating random walks in complex media. Such random walks are typical of photon propagation in turbid media, scattering of particles, i.e., neutrons in a nuclear reactor or animal/humans' migration. Among the numerous applications, Monte Carlo method is also considered a gold standard for numerically "solving" the scalar radiative transport equation even in complex geometries and distributions of the optical properties. In this work, we apply the verification method to a Monte Carlo code which is a forward problem solver extensively used for typical applications in the field of tissue optics. The method is based on the well-known law of average path length invariance when the entrance of the entities/particles in a medium obeys to a simple cosine law, i.e., Lambertian entrance, and annihilation of particles inside the medium is absent. By using this law we achieve two important points: (1) the invariance of the average path length guarantees that the expected value is known regardless of the complexity of the medium; (2) the accuracy of a Monte Carlo code can be assessed by simple statistical tests. We will show that we can reach an arbitrary accuracy of the estimated average pathlength as the number of simulated trajectories increases. The method can be applied in complete generality versus the scattering and geometrical properties of the medium, as well as in presence of refractive index mismatches in the optical case. In particular, this verification method is reliable to detect inaccuracies in the treatment of boundaries of finite media. The results presented in this paper, obtained by a standard computer machine, show a verification of our Monte Carlo code up to the sixth decimal digit. We discuss how this method can provide a fundamental tool for the verification of Monte Carlo codes in the geometry of interest, without resorting to simpler geometries and uniform distribution of the scattering properties.
| S-EPMC8484597 | biostudies-literature
Project description:Recent technical and strategical developments have increased the survival chances for avalanche victims. Still hundreds of people, primarily recreationists, get caught and buried by snow avalanches every year. About 100 die each year in the European Alps-and many more worldwide. Refining concepts for avalanche rescue means to optimize the procedures such that the survival chances are maximized in order to save the greatest possible number of lives. Avalanche rescue includes several parameters related to terrain, natural hazards, the people affected by the event, the rescuers, and the applied search and rescue equipment. The numerous parameters and their complex interaction make it unrealistic for a rescuer to take, in the urgency of the situation, the best possible decisions without clearly structured, easily applicable decision support systems. In order to analyse which measures lead to the best possible survival outcome in the complex environment of an avalanche accident, we present a numerical approach, namely a Monte Carlo simulation. We demonstrate the application of Monte Carlo simulations for two typical, yet tricky questions in avalanche rescue: (1) calculating how deep one should probe in the first passage of a probe line depending on search area, and (2) determining for how long resuscitation should be performed on a specific patient while others are still buried. In both cases, we demonstrate that optimized strategies can be calculated with the Monte Carlo method, provided that the necessary input data are available. Our Monte Carlo simulations also suggest that with a strict focus on the "greatest good for the greatest number", today's rescue strategies can be further optimized in the best interest of patients involved in an avalanche accident.
| S-EPMC5414947 | biostudies-literature
Project description:A practical detector structure is proposed to achieve stable avalanche multiplication gain in direct-conversion amorphous selenium radiation detectors.The detector structure is referred to as a field-shaping multi-well avalanche detector. Stable avalanche multiplication gain is achieved by eliminating field hot spots using high-density avalanche wells with insulated walls and field-shaping inside each well.The authors demonstrate the impact of high-density insulated wells and field-shaping to eliminate the formation of both field hot spots in the avalanche region and high fields at the metal-semiconductor interface. Results show a semi-Gaussian field distribution inside each well using the field-shaping electrodes, and the electric field at the metal-semiconductor interface can be one order-of-magnitude lower than the peak value where avalanche occurs.This is the first attempt to design a practical direct-conversion amorphous selenium detector with avalanche gain.
| S-EPMC3543372 | biostudies-other
Project description:MotivationThe binding between a peptide and a major histocompatibility complex (MHC) is one of the most important processes for the induction of an adaptive immune response. Many algorithms have been developed to predict peptide/MHC (pMHC) binding. However, no approach has yet been able to give structural insight into how peptides detach from the MHC.ResultsIn this study, we used a combination of coarse graining, hierarchical natural move Monte Carlo and stochastic conformational optimization to explore the detachment processes of 32 different peptides from HLA-A*02:01. We performed 100 independent repeats of each stochastic simulation and found that the presence of experimentally known anchor amino acids affects the detachment trajectories of our peptides. Comparison with experimental binding affinity data indicates the reliability of our approach (area under the receiver operating characteristic curve 0.85). We also compared to a 1000 ns molecular dynamics simulation of a non-binding peptide (AAAKTPVIV) and HLA-A*02:01. Even in this simulation, the longest published for pMHC, the peptide does not fully detach. Our approach is orders of magnitude faster and as such allows us to explore pMHC detachment processes in a way not possible with all-atom molecular dynamics simulations.Availability and implementationThe source code is freely available for download at http://www.cs.ox.ac.uk/mosaics/.Contactbernhard.knapp@stats.ox.ac.ukSupplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC4708099 | biostudies-literature
Project description:SignificancePerturbation and differential Monte Carlo (pMC/dMC) methods, used in conjunction with nonlinear optimization methods, have been successfully applied to solve inverse problems in diffuse optics. Application of pMC to systems over a large range of optical properties requires optimal "placement" of baseline conventional Monte Carlo (cMC) simulations to minimize the pMC variance. The inability to predict the growth in pMC solution uncertainty with perturbation size limits the application of pMC, especially for multispectral datasets where the variation of optical properties can be substantial.AimWe aim to predict the variation of pMC variance with perturbation size without explicit computation of perturbed photon weights. Our proposed method can be used to determine the range of optical properties over which pMC predictions provide sufficient accuracy. This method can be used to specify the optical properties for the reference cMC simulations that pMC utilizes to provide accurate predictions over a desired optical property range.ApproachWe utilize a conventional error propagation methodology to calculate changes in pMC relative error for Monte Carlo simulations. We demonstrate this methodology for spatially resolved diffuse reflectance measurements with ±20% scattering perturbations. We examine the performance of our method for reference simulations spanning a broad range of optical properties relevant for diffuse optical imaging of biological tissues. Our predictions are computed using the variance, covariance, and skewness of the photon weight, path length, and collision distributions generated by the reference simulation.ResultsWe find that our methodology performs best when used in conjunction with reference cMC simulations that utilize Russian Roulette (RR) method. Specifically, we demonstrate that for a proximal detector placed immediately adjacent to the source, we can estimate the pMC relative error within 5% of the true value for scattering perturbations in the range of [-15%,+20%]. For a distal detector placed at ∼3 transport mean free paths relative to the source, our method provides relative error estimates within 20% for scattering perturbations in the range of [-8%,+15%]. Moreover, reference simulations performed at lower (μs'/μa) values showed better performance for both proximal and distal detectors.ConclusionsThese findings indicate that reference simulations utilizing continuous absorption weighting (CAW) with the Russian Roulette method and executed using optical properties with a low (μs'/μa) ratio spanning the desired range of μs values, are highly advantageous for the deployment of pMC to obtain radiative transport estimates over a wide range of optical properties.
| S-EPMC10245552 | 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:We implement a Monte Carlo heuristic algorithm to model drug release from a solid dosage form. We show that with Monte Carlo simulations it is possible to identify and explain the causes of the unsatisfactory predictive power of current drug release models. It is well known that the power-law, the exponential models, as well as those derived from or inspired by them accurately reproduce only the first 60% of the release curve of a drug from a dosage form. In this study, by using Monte Carlo simulation approaches, we show that these models fit quite accurately almost the entire release profile when the release kinetics is not governed by the coexistence of different physico-chemical mechanisms. We show that the accuracy of the traditional models are comparable with those of Monte Carlo heuristics when these heuristics approximate and oversimply the phenomenology of drug release. This observation suggests to develop and use novel Monte Carlo simulation heuristics able to describe the complexity of the release kinetics, and consequently to generate data more similar to those observed in real experiments. Implementing Monte Carlo simulation heuristics of the drug release phenomenology may be much straightforward and efficient than hypothesizing and implementing from scratch complex mathematical models of the physical processes involved in drug release. Identifying and understanding through simulation heuristics what processes of this phenomenology reproduce the observed data and then formalize them in mathematics may allow avoiding time-consuming, trial-error based regression procedures. Three bullet points, highlighting the customization of the procedure. •An efficient heuristics based on Monte Carlo methods for simulating drug release from solid dosage form encodes is presented. It specifies the model of the physical process in a simple but accurate way in the formula of the Monte Carlo Micro Step (MCS) time interval.•Given the experimentally observed curve of drug release, we point out how Monte Carlo heuristics can be integrated in an evolutionary algorithmic approach to infer the mode of MCS best fitting the observed data, and thus the observed release kinetics.•The software implementing the method is written in R language, the free most used language in the bioinformaticians community.
| S-EPMC5959741 | biostudies-literature
Project description:Nonlinear state-space models are powerful tools to describe dynamical structures in complex time series. In a streaming setting where data are processed one sample at a time, simultaneous inference of the state and its nonlinear dynamics has posed significant challenges in practice. We develop a novel online learning framework, leveraging variational inference and sequential Monte Carlo, which enables flexible and accurate Bayesian joint filtering. Our method provides an approximation of the filtering posterior which can be made arbitrarily close to the true filtering distribution for a wide class of dynamics models and observation models. Specifically, the proposed framework can efficiently approximate a posterior over the dynamics using sparse Gaussian processes, allowing for an interpretable model of the latent dynamics. Constant time complexity per sample makes our approach amenable to online learning scenarios and suitable for real-time applications.
| S-EPMC10082974 | biostudies-literature
Project description:The data presented in this article describes a multi-criteria decision problem, where 13 criteria and 14 alternatives have been taken into account, consisting in the selection of an electric vehicle. The data set contains: (1) the parameters of the electric vehicles concerned included in the alternative performance model, (2) the weights of the criteria for assessing the vehicles, preference functions and thresholds constituting the preference model, (3) the overall performances and rankings of the alternatives (electric vehicles concerned). The data on vehicle parameters were collected from reports, catalogues and websites of car manufacturers and then processed into a decision table. In turn, data constituting various random preference models were generated using the Monte Carlo method. The overall performances and ranks of the alternatives were obtained using the MCDA (multi-criteria decision aid) method called NEAT F-PROMETHEE (New Easy Approach To Fuzzy Preference Ranking Organization METHod for Enrichment Evaluation), based on the performance model (decision table) and individual preference models. By linking vehicle parameters, preference models and vehicle rankings, the data allow, among other things, determining the impact of the preference model (weights of criteria, preference functions, thresholds) on the obtained vehicle rankings. The data also allow determining the probability of individual vehicles taking a specific position in the ranking on the basis of vehicle parameters, and regardless of the preferences of decision makers. Therefore, the data presented are valuable for practitioners and theorists dealing with electric vehicles and management, and in particular decision support. In the context of decision support, this data is also valuable to consumers considering the purchase of an electric vehicle, electric vehicle manufacturers, and dealers because it indicates the vehicles with the greatest market potential and user acceptance. This fact was confirmed by the research article entitled "Multi-criteria approach to stochastic and fuzzy uncertainty in the selection of electric vehicles with high social acceptance" [1] linked to this data article.
| S-EPMC8144355 | biostudies-literature
Project description:Monte Carlo is famous for accepting model extensions and model refinements up to infinite dimension. However, this powerful incremental design is based on a premise which has severely limited its application so far: a state-variable can only be recursively defined as a function of underlying state-variables if this function is linear. Here we show that this premise can be alleviated by projecting nonlinearities onto a polynomial basis and increasing the configuration space dimension. Considering phytoplankton growth in light-limited environments, radiative transfer in planetary atmospheres, electromagnetic scattering by particles, and concentrated solar power plant production, we prove the real-world usability of this advance in four test cases which were previously regarded as impracticable using Monte Carlo approaches. We also illustrate an outstanding feature of our method when applied to acute problems with interacting particles: handling rare events is now straightforward. Overall, our extension preserves the features that made the method popular: addressing nonlinearities does not compromise on model refinement or system complexity, and convergence rates remain independent of dimension.
| S-EPMC6125312 | biostudies-literature