![\"\"](\"http:\/\/pages.vassar.edu\/magnes\/files\/2019\/05\/eq1.png\")
<\/a><\/p>\nWe can approximate the double partial derivative of the wave function as:<\/p>\n
Now that we have an equation for\u00a0the double partial derivative, we can enter this back into the\u00a0Schr\u00f6dinger equation and rearrange the terms in the following way:<\/p>\n We now have an equation that allows us to compute any wavefunction given an energy and potential.<\/p>\n Matching Method<\/strong><\/p>\n Now that we can calculate a wavefunction we must find the eigenenergies associated with a system in order to determine the eigenstate. There are multiple ways to do this, one such is the matching method which is particularly useful for asymmetric potential systems. The goal of the matching function is to generate a smooth and continuous wave function by combining two wavefunctions which are generated from the left and right corners of the potential array.<\/p>\n Figure (1) – Matching method inside an infinite square well.<\/p>\n The\u00a0 inputs\u00a0 of a matching method function are an initial energy guess, the potential array, the spatial step size, and an indexing point which we will use to check the slope of each wavefunction.\u00a0 In order to calculate each wave function, there are unique initial conditions for even and odd parity solutions:<\/p>\n Using equation 3 and the initial conditions, we can calculate the first wavefunction using our initial energy guess. Then, we repeat this process by first flipping the potential array and generating a second wavefunction with the same initial conditions. We then flip this second wavefunction so that we have effectively generated one array from left to right and a second array from right to left. We then normalize each wavefunction at the indexed point so that they intersect.<\/p>\n Now we may calculate the slope of each wavefunction at the index. We then call a secondary matching function which keeps track of the orientation of the slopes. If each slope’s relative orientation remains the same, then our ‘change in energy’ value, which we will add afterwords, remains the same. If each slope’s relative orientation changes, then our ‘change in energy’ value is halved and negated. In either case, we then add a ‘change in energy’ value to the original energy guess, recalculate the left and right wavefunctions and repeat the process.<\/p>\n This function allows us to find eigenstates and hone in on the eigenenergies of a system. Once we reach a certain limit, the ‘change in energy’ value will be so small that we can stop the program and merge the two functions at the index point.<\/p>\n There are many useful applications of this method, one such is the quantum harmonic oscillator. The quantum harmonic oscillator is a useful analog for many arbitrary potentials which have stable equilibrium points.<\/p>\n Figure (2) – Top shows the matching method running on the first 6 energy levels of a quantum harmonic oscillator. Bottom shows the end result.<\/p>\n Figure (3) – The probability distribution of the 100th energy level as a result of the matching method.<\/p>\n Figure (4) – Another useful potential to study is the Lennard Jones potential. This potential is used to model the van der Waals forces between inert atoms.<\/p>\n Monte Carlo Method<\/strong><\/p>\n Another useful computational method is the Monte Carlo method. Rather than guessing at the energy as in the matching method, in this method we will be guessing at the wavefunction and iterating on this guess. The idea of a Monte Carlo method is that we randomly change a system and only keep the changes that bring us closer to a solution we would expect. In the case of quantum mechanics, this will be used to randomly change the wavefunction and only keep the changes that lower the energy of the system.<\/p>\n In order to implement this approach, we must first understand how to calculate energy using equation 4 below:<\/p>\n In order to preform these integrals in discrete space, we can use several different methods, but for this process I chose to use the trapezoidal integration method:<\/p>\n Applying this method to the quotient and denominator of equation 4 leaves us with the following equations:<\/p>\n Now equation 4 becomes:<\/p>\n But now we have energy on both sides of the equation. In order to resolve this, we need to substitute in\u00a0the time independent Schr\u00f6dinger equation (Eq. 1) which gives us<\/p>\n Finally, in order to get rid of the double partial derivatives, we substitute in equation 2 and find:<\/p>\n This is the final solution which will allow us to calculate the energy in any system given a wavefunction and potential array.<\/p>\n Now that we can calculate energy, we can utilize the Monte Carlo method. This method is actually quite simple; we create a function that takes in an arbitrary guess of the wavefunction and calculate the energy of this guess. Then we randomly pick a point of the wavefunction and change it by a random amount. Then we recalculate the energy of the system and if the energy has decreased, we can keep that random change, otherwise discard the change and revert to the previous wavefunction. We repeat this process over and over again until we reach a point in which the energy does not change after a certain number of attempts. At this point we can assume the we reached the ground state of our potential system.<\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a><\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a> <\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n<\/a><\/p>\n