Monte carlo simulation Definition and 25 Threads

Hello everyone, I would like to start a technical discussion about the TMESH tally in MCNP versions 6.2 and 6.3. From my studies, TMESH is similar to +F6, tied to a mesh, and its results are given in terms of MeV/cm³. I wanted to know if there is any way (a multiplier or conversion factor) to...

I'm working on this Monte Carlo simulation problem for radiation propagation. So far I think I have accomplished creating a histogram of simulated values (iii). Now I need to plot it against the exact PDF. I need to find the normalized PDF, but can't obtain the normalization constant, and so on...

Hello everyone, I am facing difficulties while trying to calculate the dose in the tibia due to brachytherapy in an MCNP simulation. We are working with the radionuclide Ho-166, and therefore, we need to account for both photon and electron contributions to the deposited dose. Initially, I...

Dear experts , when I run the ICRP145_Human Phantoms example following the instructions in the user manual, I get this error response [ VirtualBox:~/Geant4/share/Geant4/examples/advanced/ICRP145_HumanPhantoms/build$./ICRP145phantoms ./Internal -i 9500 -m example.in -o example.out Usage...

Dear experts, I would like to do a double differential cross section calculation for iron at 90 degrees and at 22 MeV in Fluka. What changes should I make in the input and output files of the application and what transformations should I make in order to see the double differential cross section...

here is my attempt to implement using python import numpy as np import matplotlib.pyplot as plt def initialize_spins(L): """Initialize a random spin configuration with unit magnitudes.""" spins = np.random.normal(size=(L, L, L, 3)) magnitudes = np.linalg.norm(spins, axis=-1...

Hi everyone, I've been trying to analyze PTRAC output file from MCNP6 here we can see the location , cell, particle, time, and so on... My question is, I have trouble finding the unit of time listed in PTRAC, (ex, 0.30113E-02), which is hard to find in MCNP manual My intuition is that the...

As far as I understand: Pruning means deleting something. Enrich means to enhance/increase weight. Now my question is, in the case of the PERM algorithm, Are we deleting some polymers chains? Or, are we deleting some beads in all polymer chains? Are we increasing the weights of some...

I am absolutely new to Polymer simulation. I am trying to understand the simulation by analyzing source code written by others. Can anyone tell me what are the differences between the following three source code in terms of their objectives? Monte-Carlo-simulation-of-polymers...

Hi everyone, I want to generate 8 random variables (in reality to form 4 complex numbers) such that the sum of the 8 variables squared is equal to unity. The aim of generating such numbers is to perform a quantum simulation of 4 qubits (thus the 8 parameters). I've been trying to use...

Homework Statement Lo,Im stuck on how to retrieve the specific heat capacity from an MC simulation, with the metropolis algorithm. I want my graph to look something like this: https://i.stack.imgur.com/NXeXs.png Homework Equations C_v = ((<E^2>-<E>^2)/T^2 The Attempt at a Solution My code is...

I am an undergraduate working in my university's high energy physics lab. I'm still new to particle physics, so I have to learn/understand a bunch of things on the fly. We are currently searching for a new state theorized to occur in a specific decay process. I know we have to work with data...

Trying to accomplish a monte carlo simulation on the condensed state of 4He, yet I am in my sophomore year and know only a bit of quantum statistical physics. Is there any documentations recommended for beginners to the algorithm applied to 4He? I've found some but they are not friendly to...

Hi there! I'm new here and have a pretty basic and maybe stupid question: What exactly does one mean when talking about "truth information"? I'm working with a simulation of a particle detector and I've been trying to make sense out of all the data that is stored in the ouput files. Words like...

Hey, When performing Molecular dynamic or Monte Carlo simulation (in NPT or NVT ensemble), I'm wondering whether there is any difference between the average energy of the system and the energy of the average structure. If there is a difference, how munch should it be? and why? Does the the...

Some of the social sciences suffer from "physics envy". This malady causes educators to inject an unnecssary amount of mathematics into the curriculum as a way of gaining scientific letgitimacy. Sadly for most undergrads, the math actually gets in the way. I wrote a paper in which I describe the...

Hello all, I have been working on an NPT monte carlo simulation. I would like to know how I can measure the instantaneous pressure of the system at each monte carlo step?

Hello, 1. Does anybody know which book that gives a good explanation about monte carlo simulation? I've read many tutorials, it mentions about its error, accuracy, variance. But, many of them don't actually show, how to perform monte carlo simulation. Questions 2. Is it actually the same...

The Monte Carlo simulation is a very important tool in particle physics specially to tune in the preciseness of the real time experiment. In particle physics, I had an opportunity to work on the data analysis of neutrino flux produced from the g4numi (the Neutrino beam from main injector, in...

EDIT: I'm not in my sharpest moment. I just found a bunch of posts that discuss this. I'll read them and update this post if I find an answer. I'm working on a hard-sphere MC simulation (for a class). I want to compute the radial distribution function g(r). To put you in context, as to my...

I'm currently working on a project for a code that does umbrella sampling of a 2-D Ising model [size LxL of a magnet (analyzing up or down spins)]. The next step is to take my code to analyze a nucleation region and its growth by varying the temperature above critical. Before I even attempt to...

Hi all! I started learning about Monte Carlo Simulation. However, one thing that I don't quite get is that why for generating any random variable, we have to first generate a Uniform RV? What is the reason behind that? Thanks!

I am trying to model a simple birth and death process in Mathematica using the Gillespie Algorithm. I am using 1 DNA molecule that is transcribed to mRNA with rate k1, \mbox{DNA} \longrightarrow \mbox{DNA + RNA} and the transcribed RNA are subject to degradation with rate k2...

Hi all... Im trying to model a simple lattice-gas ( ising model ) system , using MC metropolis algorithm. I am having trouble with the code. I am trying to do it in matlab. My lattice gas does not show the characteristics it should. I used my old Ising Model code, which I modified to do...

Hello, my colleague has a problem. He has a Monte Carlo simulation written in C. But his program gives different outputs in case he compile it using g++ or icc. He has nothing changed in his source. He is using standard random number generator drand() with initialization srand((unsigned) time...