where r is the distance between two atoms, $\epsilon$ is the depth of the potential well, and $\sigma$ is the distance at which the inter-particle potential is zero.
where r is the distance between two atoms, $\epsilon$ is the depth of the potential well, and $\sigma$ is the distance at which the inter-particle potential is zero.
1. Sketch U(r) as a function of interatomic distance and mark the regions of repulsive and attractive forces acting between atoms.
1. Sketch U(r) as a function of interatomic distance and mark the regions of repulsive and attractive forces acting between atoms.
2. Find the distance, $r_0$ (bond length) at which the potential energy is minimal and find the value of the potential energy at this distance (binding energy of the molecule).
3. Expand U(r) in a Taylor's series around $r_0$ up to second order. By considering harmonic potential around the minimum ($r_0$), find an expression for the spring constant, $K$, in terms of $\epsilon$ and $\sigma$.
4. In the harmonic potential approximation of Lennard-Jones potential of argon, find the ground state energy of the molecule. What is the energy required to break the molecule apart?
...
...
@@ -252,7 +252,7 @@ For this exercise, use Matlab (or any other programming environment you are comf
??? hint
In Matlab use the function [`diag`](https://nl.mathworks.com/help/matlab/ref/diag.html).
In Matlab use the function [`diag`](https://nl.mathworks.com/help/matlab/ref/diag.html).
2. Using numerical diagonalization ([`eig`](https://nl.mathworks.com/help/matlab/ref/eig.html)), compute the eigenfrequencies of this atomic chain. Plot a histogram of these eigenfrequencies.
3. Make the masses of every even atom different from the masses of every odd atom. Compute the eigenfrequencies of this atomic chain and plot a histogram.