Search results for "Monte carlo method"

2020

Discriminant validity was originally presented as a set of empirical criteria that can be assessed from multitrait-multimethod (MTMM) matrices. Because datasets used by applied researchers rarely lend themselves to MTMM analysis, the need to assess discriminant validity in empirical research has led to the introduction of numerous techniques, some of which have been introduced in an ad hoc manner and without rigorous methodological support. We review various definitions of and techniques for assessing discriminant validity and provide a generalized definition of discriminant validity based on the correlation between two measures after measurement error has been considered. We then review t…

Quantifying sensitivity and uncertainty analysis of a new mathematical model for the evaluation of greenhouse gas emissions from membrane bioreactors

2015

Abstract A new mathematical model able to quantify greenhouse gas (GHG) emissions in terms of carbon dioxide (CO 2 ) and nitrous oxide (N 2 O) for a Membrane Bioreactor (MBR) is presented. The proposed mathematical model is of the Activated Sludge Model (ASM) family and takes into account simultaneously both biological and physical processes (e.g., membrane fouling). An analysis of the key factors and sources of uncertainty influencing GHG emissions is also presented. Specifically, the standardized regression coefficient, the Extended-FAST and a Monte Carlo based method are employed for assessing model factors which influence three performance indicators: effluent quality index, operational…

Particle Group Metropolis Methods for Tracking the Leaf Area Index

2020

Monte Carlo (MC) algorithms are widely used for Bayesian inference in statistics, signal processing, and machine learning. In this work, we introduce an Markov Chain Monte Carlo (MCMC) technique driven by a particle filter. The resulting scheme is a generalization of the so-called Particle Metropolis-Hastings (PMH) method, where a suitable Markov chain of sets of weighted samples is generated. We also introduce a marginal version for the goal of jointly inferring dynamic and static variables. The proposed algorithms outperform the corresponding standard PMH schemes, as shown by numerical experiments.

Dosimetry comparison between TG-43 and Monte Carlo calculations using the Freiburg flap for skin high-dose-rate brachytherapy

2012

Abstract Purpose The purpose of this work was to evaluate whether the delivered dose to the skin surface and at the prescription depth when using a Freiburg flap applicator is in agreement with the one predicted by the treatment planning system (TPS) using the TG-43 dose-calculation formalism. Methods and Materials Monte Carlo (MC) simulations and radiochromic film measurements have been performed to obtain dose distributions with the source located at the center of one of the spheres and between two spheres. Primary and scatter dose contributions were evaluated to understand the role played by the scatter component. A standard treatment plan was generated using MC- and TG-43-based TPS appl…

Correction factors for ionization chamber measurements with the ‘Valencia’ and ‘large field Valencia’ brachytherapy applicators

2018

Treatment of small skin lesions using HDR brachytherapy applicators is a widely used technique. The shielded applicators currently available in clinical practice are based on a tungsten-alloy cup that collimates the source-emitted radiation into a small region, hence protecting nearby tissues. The goal of this manuscript is to evaluate the correction factors required for dose measurements with a plane-parallel ionization chamber typically used in clinical brachytherapy for the 'Valencia' and 'large field Valencia' shielded applicators. Monte Carlo simulations have been performed using the PENELOPE-2014 system to determine the absorbed dose deposited in a water phantom and in the chamber act…

Scientific simulations and optimization of the XGIS instrument on board THESEUS

2020

The XGIS (X and Gamma Imaging Spectrometer) is one of the three instruments onboard the THESEUS mission (ESA M5, currently in Phase-A). Thanks to its wide field of view and good imaging capabilities, it will efficiently detect and localize gamma-ray bursts and other transients in the 2-150 keV sky, and also provide spectroscopy up to 10 MeV. Its current design has been optimized by means of scientific simulations based on a Monte Carlo model of the instrument coupled to a state-of-the-art description of the populations of long and short GRBs extending to high redshifts. We describe the optimization process that led to the current design of the XGIS, based on two identical units with partial…

Reduction of stored-particle background by a magnetic pulse method at the KATRIN experiment

2018

Arenz, M., et al. “Reduction of Stored-Particle Background by a Magnetic Pulse Method at the KATRIN Experiment.” The European Physical Journal C, vol. 78, no. 9, Sept. 2018. © 2018 The Authors

Shell structure and level spacing distribution in metallic clusters

1993

The lattice gas Monte Carlo and tight binding method is used to study the electronic shell structure in large metallic clusters. The average density of states of a large ensemble of deformed clusters shows the same shell structure as the most spherical geometry. The level spacing distribution at the Fermi level is a Wigner distribution.

Static and dynamic glass transitions in the 10-state Potts glass: What can Monte Carlo simulations contribute?

2002

The p-state Potts glass with infinite range Gaussian interactions can be solved exactly in the thermodynamic limit and exhibits an unconventional phase behavior if p >4: A dynamical transition from ergodic to non-ergodic behavior at a temperature T D is followed by a first order transition at T 0 < T D, where a glass order parameter appears discontinuously, although the latent heat is zero. If one assumes that a similar scenario occurs for the structural glass transition as well (though with the singular behavior at T D rounded off), the p-state Potts glass should be a good test case to develop methods to deal with finite size effects for the static as well as the dynamic transition, and to…

Performance potential for simulating spin models on GPU

2012

Graphics processing units (GPUs) are recently being used to an increasing degree for general computational purposes. This development is motivated by their theoretical peak performance, which significantly exceeds that of broadly available CPUs. For practical purposes, however, it is far from clear how much of this theoretical performance can be realized in actual scientific applications. As is discussed here for the case of studying classical spin models of statistical mechanics by Monte Carlo simulations, only an explicit tailoring of the involved algorithms to the specific architecture under consideration allows to harvest the computational power of GPU systems. A number of examples, ran…