My Coding >
Numerical simulations >
Atoms and Molecules >
Tools for molecular modelling >
2D random distribution with density function by Rejection Sampling techniques
2D random distribution with density function by Rejection Sampling techniquesIn the previous example, I've explained how to generate random distribution for Ndots in the onedimensional case in the presence of the density function. This was a very simple algorithm, but unfortunately, it was not suitable for multidimensional space. In this example, I will show how to solve this problem in a 2D case. And also, this 2D case can be extended to any dimension. Explanation with coding I've shown in the video.
Rejection Sampling techniquesWhat is the idea of rejection sampling? It is very easy. We generate a random dot in the given space and then we calculate the probability, or density in this spot, according to our function. Probability functionThe probability, or density function should be scaled to the same value as a later random threshold. To make everything easier, it is better to scale everything to 1.
In this example, my function calculates random distribution (or density) with maximal value = 1. Random samplingIn the next step, we generate a random dot in our range in the multidimensional space.
And now, we need to make a random test for this dot. We need to generate random values in the same range and compare it with calculated probability. If the calculated probability is bigger than this random value, then we will accept this sample. Otherwise, we will reject it.
This is an idea of this method. Now you can watch the full video to see how to write the full code and display the results. Full code for Rejection Sampling techniqueThis is a full code for the Rejection Sampling technique with some examples of visual representation
As you can see, we've calculated the density distribution for electrons in 2D space. In the next step, I will show you how to apply this code for the 3D case and how to visualise results for an electron cloud in the hydrogen atom. I think for a full understanding of this code, you also need to look into my lecture about meshgrid.

Last 10 artitles
9 popular artitles


© 2020 MyCoding.uk My blog about coding and further learning. This blog was writen with pure Perl and frontend output was performed with TemplateToolkit. 