Related papers: The scattering order problem in Monte Carlo radiat…
The possibility of different interpretations of the stochastic term (or calculi) in the overdamped Langevin equation for the motion of a particle in an inhomogeneous medium is often referred to as the "Ito--Stratonovich dilemma," although…
Monte Carlo simulations are performed to study the in-plane transport of spin-polarized electrons in III-V semiconductor quantum wells. The density matrix description of the spin polarization is incorporated in the simulation algorithm. The…
Monte Carlo simulations are based on the manipulation of random numbers to evaluate probable outcomes, with applicability in a variety of different fields. By assigning probabilities, which can be determined a priori, to various events, it…
We introduce the radiative transfer code Sweep for the cosmological simulation suite Arepo. Sweep is a discrete ordinates method in which the radiative transfer equation is solved under the infinite speed of light, steady state assumption…
To recover the three dimensional (3D) volumetric distribution of matter in an object, images of the object are captured from multiple directions and locations. Using these images tomographic computations extract the distribution. In highly…
Emission features formed through Raman scattering with atomic hydrogen provide unique and crucial information to probe the distribution and kinematics of a thick neutral region illuminated by a strong far UV emission source. We introduce a…
Preliminary low-thrust spacecraft mission design is a global search problem characterized by a complex solution landscape, multiple objectives, and numerous local minima. During this phase, mission parameters are often not yet fully…
Neutron transport along guides is governed by the Liouville theorem and the technology involved has advanced in recent decades. Computer simulations have proven to be useful tools in the design and conception of neutron guide systems in…
We present a new efficient method for Monte Carlo simulations of diffusion-reaction processes. First introduced by us in [Phys. Rev. Lett., 97:230602, 2006], the new algorithm skips the traditional small diffusion hops and propagates the…
The efficient evaluation of high-dimensional integrals is of importance in both theoretical and practical fields of science, such as data science, statistical physics, and machine learning. However, exact computation methods suffer from the…
Multi-MeV flash radiography is often used as the primary diagnostic technique for high energy and density (HED) physics experiments. Primary X-ray which is attenuated by the object offers density information of the object. For a thick metal…
Significance: Changes in the structure of tissue occur in many disease processes, such as the boundaries of cancerous tumors and burn injuries. Spectroscopic and polarimetric alterations of terahertz light caused by Mie scattering patterns…
The validation of a previous work on the optimization of material discrimination in spectral x-ray imaging is reported. Using Monte Carlo simulations based on the BEAMnrc package, material decomposition was performed on the projection…
Multifidelity Monte Carlo methods often rely on a preprocessing phase consisting of standard Monte Carlo sampling to estimate correlation coefficients between models of different fidelity to determine the weights and number of samples for…
Monte Carlo simulations are widely used in many areas including particle accelerators. In this lecture, after a short introduction and reviewing of some statistical backgrounds, we will discuss methods such as direct inversion, rejection…
Synthesizing realistic images involves computing high-dimensional light-transport integrals. In practice, these integrals are numerically estimated via Monte Carlo integration. The error of this estimation manifests itself as conspicuous…
By decomposing the important sampled imaginary time Schr\"odinger evolution operator to fourth order with positive coefficients, we derived a number of distinct fourth order Diffusion Monte Carlo algorithms. These sophisticated algorithms…
The computational cost of a Monte Carlo algorithm can only be meaningfully discussed when taking into account the magnitude of the resulting statistical error. Aiming for a fixed error per particle, we study the scaling behavior of the…
Due to its accuracy and generality, Monte Carlo radiative transfer (MCRT) has emerged as the prevalent method for Ly$\alpha$ radiative transfer in arbitrary geometries. The standard MCRT encounters a significant efficiency barrier in the…
Although Monte Carlo path tracing is a simple and effective algorithm to synthesize photo-realistic images, it is often very slow to converge to noise-free results when involving complex global illumination. One of the most successful…