Homework4

Homework4 - University of Virginia, MSE 4270/6270:...

Info iconThis preview shows pages 1–3. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: University of Virginia, MSE 4270/6270: Introduction to Atomistic Simulations, Leonid Zhigilei Homework #4 1. Review computer code that implements a simple Ising model, mse627-mc.f90, (the code was written recently and not thoroughly tested, if you find any bugs or get ideas for improvement please send them to the class mailing list). 2. Perform Metropolis Monte Carlo simulation (MC moves = switch particle types) to find equilibrium structure (and compositional ordering) for binary alloys with two sets of parameters (E AA = E BB = -0.05 eV, E AB = -0.01 eV) and (E AA = E BB = -0.05 eV, E AB = -0.25 eV) and several temperatures. Choose values of temperature that would allow you to discuss the temperature dependence of the results of the simulations. 3. Plot snapshots from your simulations to illustrate the structural changes during the simulations and the final equilibrium structures observed for different input parameters used (task #2). ( Snapshots can be made using the files written to sub-directory data you should create this directory before running the code ). Make plots showing changes of energy, number of AA, AB, and BB bonds during the simulations. ( Plots can be made using the data written to file EnBnd.out ). Discuss your results based on thermodynamics of binary solid solutions briefly outlined in the next two pages. 4. Adapt the code to model a crystal with vacancies. Perform simulations for a system with energy of vacancy formation of 1 eV. Calculate the equilibrium vacancy concentrations for several (2-3) values of temperature. Compare your results to the theoretical equation derived/discussed below. Discuss agreement/disagreement between the numerical results and the theoretical equation. Objective: Getting experience with Metropolis Monte Carlo simulations, using Ising model to study compositional ordering and segregation in binary alloys. University of Virginia, MSE 4270/6270: Introduction to Atomistic Simulations, Leonid Zhigilei Lets consider mixing of two components A and B:...
View Full Document

This note was uploaded on 03/09/2012 for the course MSE 4270 taught by Professor Zhigilei during the Fall '11 term at UVA.

Page1 / 7

Homework4 - University of Virginia, MSE 4270/6270:...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online