Monte carlo simulation Definition and 30 Threads

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...

Introducing the spacetime spherical symmetric lattice, I use the following notifications in my program. i - index enumerating the nodes along t-coordinate, j - along the r-coordinate, k - along the theta-coordinate, l - along the phi-coordinate. N_t - the number of nodes along t-coordinate. N_r...

Hi guys. I'm using the Monte Carlo method to simulate a spin lattice. If I have a square lattice, L x L, I can plot the phase transition temperature by the inverse of the lattice length (1/L) to find the phase transition temperature in the thermodynamic limit (extrapolating the curve for 1/L =...

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 [/B] I'm trying to write a program for caclulating Green's function using Monte Carlo method (Metropolis algorithm) in scalar field theory with a potential λφ4 in 4D. I'm writing it in python. N_t, N_x, N_y, N_z - total number of lattice sites in each directions. Field...

Homework Statement The number of busy lines in a trunk group (Erlang system) is given by a truncated Poisson distribution. I am asked to generate values from this distribution by applying the Metropolis-Hastings algorithm. Homework Equations The distribution is given in the attached picture...

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...

Homework Statement Calculate 5-day 1% Value at Risk of a portfolio using Monte Carlo simulation. Homework EquationsThe Attempt at a Solution I've found an article explaining how to perform Monte Carlo simulation on portfolio returns therefore calculating the Value At Risk (pdf attachment)...

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!

Hi there, I have a question about probability, conditional expectations, copula functions and their use in a Monte carlo simulation. I'd appreciate any help or comments that you can offer. I'll describe briefly what I am trying to simulate and I'll ask the specific question at the end...

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...

Let \bold{X} be a discrete random variable whose set of possible values is \bold{x}_j, \ j \geq 1 . Let the probability mass function of \bold{X} be given by P \{\bold{X} = \bold{x}_j \}, \ j \geq 1 , and suppose we are interested in calculating \theta = E[h(\bold{X})] =...

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...