exercises:2014_uzh_molsim:h2o_diff
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
exercise:mm_uzh:h2o_diff [2014/05/12 20:26] – talirz | exercises:2014_uzh_molsim:h2o_diff [2020/08/21 10:15] (current) – external edit 127.0.0.1 | ||
---|---|---|---|
Line 2: | Line 2: | ||
When simulating liquids or solids under periodic boundary conditions, we are making two fundamental approximations: | When simulating liquids or solids under periodic boundary conditions, we are making two fundamental approximations: | ||
- | - We simulate an infinite system, thus neglecting the fact that any real-world system | + | - We simulate an infinite system, thus neglecting the fact that any real-world system |
- We impose the condition that the properties of the system under study repeat //exactly// from one simulation cell to the next. The quality of this approximation depends on the system under study and the quantity of interest. | - We impose the condition that the properties of the system under study repeat //exactly// from one simulation cell to the next. The quality of this approximation depends on the system under study and the quantity of interest. | ||
Line 20: | Line 20: | ||
- We want to simulate diffusion at room temperature. Why aren't we using the $NVT$ ensemble? //Hint:// Think about how thermostats work. | - We want to simulate diffusion at room temperature. Why aren't we using the $NVT$ ensemble? //Hint:// Think about how thermostats work. | ||
- Use the provided script '' | - Use the provided script '' | ||
- | - How are temperature fluctuations expected to depend on system size? Use gnuplot' | + | - How are temperature fluctuations expected to depend on system size? Use gnuplot' |
</ | </ | ||
Line 41: | Line 41: | ||
- We have precalculated trajectories for 64, 128 and 256 water molecules (ask your teaching assistant). Use '' | - We have precalculated trajectories for 64, 128 and 256 water molecules (ask your teaching assistant). Use '' | ||
- | - Plot the msd as a function of time on a double logarithmic scale. Can you identify different regimes? Why does the signal become noisy towards long times? | + | - Plot the msd as a function of time on a double logarithmic scale. Can you identify different regimes? Why does the signal become noisy towards long times? |
- Obtain the diffusion constant $D_{pbc}$ by fitting a line through the mean square displacement data in the range $2-10$ ps. | - Obtain the diffusion constant $D_{pbc}$ by fitting a line through the mean square displacement data in the range $2-10$ ps. | ||
- Compare against the values in Table I of the article. //Note:// We are using a slightly different force field, but the values should be of a similar magnitude. If not, check your units! | - Compare against the values in Table I of the article. //Note:// We are using a slightly different force field, but the values should be of a similar magnitude. If not, check your units! | ||
Line 52: | Line 52: | ||
- Plot $D_{PBC}$ as a function of $1/L$, where $L$ is the length of the edge of the simulation box. | - Plot $D_{PBC}$ as a function of $1/L$, where $L$ is the length of the edge of the simulation box. | ||
- Perform a linear fit of this curve to obtain the diffusion constant $D=D_{pbc}(L=\infty)$ | - Perform a linear fit of this curve to obtain the diffusion constant $D=D_{pbc}(L=\infty)$ | ||
- | - Use Eq. 12 in the article to calculate the viscosity. | + | - Use equation (12) in the article to calculate the viscosity |
- Compare the results to the data in the paper. | - Compare the results to the data in the paper. | ||
</ | </ |
exercises/2014_uzh_molsim/h2o_diff.1399926390.txt.gz · Last modified: 2020/08/21 10:14 (external edit)