syllabus - iitg.ac.in
TRANSCRIPT
![Page 1: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/1.jpg)
Syllabus:
Text Books:
1. Nano-Engineering in Science and Technology, M. Reith.
2. Computer Simulation of Liquids, M. P. Allen and D. J. Tildesley.
3. Understanding Molecular Simulation, D. Frenkel and B. Smith. 4. Molecular
Dynamics Simulation, J. M. Haile.
![Page 2: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/2.jpg)
Molecular DynamicsPreliminaries
![Page 3: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/3.jpg)
A Brief History of MD
�1957, 59 Alder-Wainwright, introduction of basic MD
- Of hard sphere particles.
�1964 Rahman, study of liquid Ar (NVE-MD)
- First quantitative study reported.
�1980, 81 Andersen, Parrinello-Rahman
- Constant Pressure (NPT) MD.- Constant Pressure (NPT) MD.
�1984, 86 Nose, Hoover
- Constant Temperature (NVT) MD.
�1985 Car-Parrinello, ab initio MD
- Based on Density Functional Theory (1960s)
� Hohenbrg-Sham
� Kohn-Sham
![Page 4: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/4.jpg)
A movie
System:
Na3Zr2Si2PO12
Na-SuperIonic Conductor
Na+ ions diffuse in the
[Zr2Si2PO12 ] 3- matrix[Zr2Si2PO12 ] 3- matrixconducting electricity.
-A promising material for
better batteries!
![Page 5: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/5.jpg)
![Page 6: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/6.jpg)
Ion Channel’s In NASICON’s
Na3Zr2Si2PO12
Supriya Roy and PKP to be published.
![Page 7: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/7.jpg)
Ion Channel Recognized
System:
Na3Zr2Si2PO12
Na-SuperIonic Conductor
Na+ ions diffuse in the
[Zr2Si2PO12 ] 3- matrix[Zr2Si2PO12 ] 3- matrixconducting electricity.
-A promising material for
better batteries!
![Page 8: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/8.jpg)
…Structural Properties: Site Occupancies
PRL 97, 166401 (2006)
![Page 9: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/9.jpg)
Advanced Simulation Techniques: Simulating Chemical Reaction Dissociation of H2CO3 in gas-phase
H2CO3=H2O + CO2
Indian Institute of Technology Guwahati
![Page 10: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/10.jpg)
Advanced Simulation Techniques: Understanding local structure and reactivity; Carbonates in Solution
![Page 11: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/11.jpg)
Simulation of Life!!
Indian Institute of Technology Guwahati
Charge distribution around the virus
![Page 12: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/12.jpg)
∑≠
=ij
iji FF
••
•2
3
F12
F23
F13
F1 = F12 + F13
F2 = F21 + F23
Molecular Dynamics –The Basic Idea
1
iii mFa /=
•3F13
F3 = F31 + F32
F2 = F21 + F23
Newton’s IInd Law:
ijijUF −∇=
![Page 13: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/13.jpg)
••
•1
2
3
Progress in time…
x(t) → x(t+∆t)
y(t) → y(t+∆t)
z(t) → z(t+∆t)• ••
3
1
2
•3
• •3Taylor Expansion: too crude to
use it as suchas such!!
∆t ~ 1-5 fs (10-15sec)
Advance positions & velocities of each atom:
![Page 14: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/14.jpg)
A good Integrator
Verlet Scheme:
Newton’s equations are time reversible,
Summing the two equations,Summing the two equations,
Now we have advanced our atoms to time t+∆t ! !
Velocity of the atoms:
![Page 15: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/15.jpg)
• ••
3
1
2
• ••
3
1
2F12
F23
F13
…Atoms move forward in time!
1. Calculate 2. Update3. Update •3•3F13
t+Δt t+2Δt
Continue this procedure for several lakhs of Steps. Or as much as you can afford!
2nd MD Step 3rd MD Step
The main O/P of MD is the trajectory.
∆t ~ 1-5 fs (10-15sec)
![Page 16: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/16.jpg)
The missing ingredient… Forces?
Gravitational Potential: r
mGmU 21−
=too weak,
Neglect it!!
Force is the gradient of potential:
The predominant inter-atomic forces are Coulombic in origin!
r
qqU 21
04
1πε
=
However, this pure monopole interaction need not be present!!
![Page 17: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/17.jpg)
Interatomic forces for simple systems
1. Lennard-Jones Potential:
Gives an accurate description of inert gasesinert gases
(Ar, Xe, Kr etc.)
(non-bonded interactions)
Instantaneous dipoles
2. Born-Mayer (Tosi-Fumi) Potential:
(Ar, Xe, Kr etc.)
Faithful in describing purepure ionic solids (NaCl, KCl, NaBr etc.)
![Page 18: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/18.jpg)
The Lennard-Jones Potential
for Ar :
(kB)
Pauli’s repulsion
ijijUF −∇=
Å
Should be known apriori.
“dispersion” interaction
r(nm)
))(612
(48
6
14
12
ji
xij
xxrr
F −−=σσ
ε
iii
xij x
r
r
U
x
UF
∂∂
∂∂
−=∂∂
−=
![Page 19: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/19.jpg)
Length and Times of MD simulation
Typical experiment sample (bulk) contains ~ 1023 atoms!
Typical MD simulations (on a single CPU)
a) Can include 1000 – 10,000 atoms (~20-40 Å in size)!
b) run length ~ 1 –10 nsns (10-9 seconds)!
Consequence of system size:Consequence of system size:
Larger fraction of atoms are on the surface, r
dr
mr
mdrr
N
Ns3
/34
/4
3
2
==ρπ
ρπ
78
10~10
)3(3~.)( −
Α
ΑExpt
N
Ns45.0~
20)3(
3~)(Α
ΑMD
N
Ns
Surface atoms have different environment than bulk atoms!
Bulk system Nano-cluster
![Page 20: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/20.jpg)
The Simulation Cell
Insert the atoms in a perfectly porous box – simulation super-cell.
The length of the box is determined as, L3 = M/Dexp = N*m/Dexp
Dexp = Expt. density; m = At. mass; N= No. of atoms;
If crystal structure/unit cell parameters are unknown (eg., liquids)
For Bulk systems!
L Assign positions and velocities(=0) for each atom.
![Page 21: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/21.jpg)
Periodic Boundary Condition
In 3-D the simulation simulationsimulation--supersuper--cellcell is surrounded by 26 image image cellscells!!
Y
Image coordinates:
Construct Periodic Images:
For Bulk systems!
3D – torus!!
~20 AX
L~20 A
n1, n2, n3 ∈∈∈∈∈∈∈∈ -1,0, 1,
x’ = x + n1 L
y’ = y + n2 L
z’ = z + n3 L
Image coordinates:
Now, there are no surface atoms!!
![Page 22: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/22.jpg)
Minimum Image Convention
Rc
Interactions between atoms separated by a chosen cut-off distance (Rc) or larger (ie, rij > Rc) are neglected.
Rc is chosen such that U(Rc) ~ 0
~20 A
~20 A
A large enough system (ie, bigger simsim..--cellcell) is chosen such that Rc ≤≤≤≤ L/2.
Thus particle ii interact either
with particle j j or one of its images, but not bothbut not both!
![Page 23: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/23.jpg)
•••• •
•••••
••
••
••
•
•••
t = 0
As Time Progress…
iiiiii
jj
•
•jj’’
L••
••
•t = 0
t = 10 pico-sec
How to find the distance to the image of jj that is nearest to ii ??
jj
))(612
(48
6
14
12
jiijij
xij
xxrr
F −−=σσ
εForce between i i & jj ::
•
L
L••
![Page 24: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/24.jpg)
The three lines of code…
dx = x(j) – x(i)
dy = y(j) – y(i)
dz = z(j) – z(i)
dx = dx – L*ANINT(dx/L)
dy = dy – L*ANINT(dy/L)
Define,
ANINT(1.51) = 2
ANINT(3.49) = 3
ANINT(-3.51) = -4
ANINT(-11.39) = -11.0
FORTRAN 77/90/95
dy = dy – L*ANINT(dy/L)
dz = dz – L*ANINT(dz/L)ANINT(-11.39) = -11.0
rij_2 = (dx**2 + dy**2 + dz**2)
))(612
(48
6
14
12
jiijij
xij
xxrr
F −−=σσ
ε
rij_8 = rij_2**4 rij_14 = rij_2**7
![Page 25: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/25.jpg)
•
…demonstration
ii•(2,8)
(-1,6)
dx = dx – L*AINT(dx/L)
dy = dy – L*AINT(dy/L)
• (9,6)
•jj
•jj’’
L=10 A
L=10 A
(0,0)(19,-4)
xi-xj=(2-19)-10*ANINT{(2-19)/10}= -17 - (-10*2)=3
•j”j”
yi-yj=(8+4)-10*ANINT{12/10}= 12 - (10*1)=2
(9,6)
• •
x’ = x + n1 L
y’ = y + n2 L
![Page 26: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/26.jpg)
The structure of a simple MD code
![Page 27: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/27.jpg)
Remarks on Statistical Ensemble
There is no energy coming in Or going out of our system of atoms: micro-canonical (NVE) ensemble.
Thus the total energy (E) and total linear momentum of the system should be conserved – through out our simulation!
∑=N
vmK21 (fluctuates)∑
==
ii
vi
mK
1
221
∑=
=N
ii
UU12
1
UKE += remain const.
(fluctuates)
(fluctuates)
Transient phase Equilibrated state
![Page 28: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/28.jpg)
Remarks on Energy
• Start from the expt. crystal structure if available.
• Else? Start from good guess! (like, in bio-systems, polymers, liquids)
-good check on your code!
-time integrator!
-on time step (∆t) used!
6/ ~ 10 /E
ps psE
−∆• Energy conservation,
• Else? Start from good guess! (like, in bio-systems, polymers, liquids)
• And, perform an energy minimization!
(Routines available in standard packages.
Or, do an MD with constant velocity scaling.)• Reaching a well equilibrated structure (from a random structure)
can be very very costly!
• Fluctuation of U(t) about a mean helps to identify equilibrated system.
![Page 29: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/29.jpg)
Calculating Temperature
Equipartition theorem:
Average temperature : 1
1( )
M
m
mT T tM =
∑< >=
∑=N
iii
vmtTNk2
21
)(23
InstantaneousInstantaneous temperature: ∑=N
iii
vmNk
tT2
31
)(
1mM =
M – no. of MD steps performed
• ••
1
2
3
Even if we start with vi =0, the system picks up non-zero v (hence some TT) as time progress!
This some T(!T(!) need not be what we want!
So how do we control TT?
![Page 30: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/30.jpg)
Controlling the Temperature ?
Velocity rescaling:
∑=N
iii
vmNk
T2
31Actual temp. at some instant.
TTTTTrr
∆+><∆−If T is out side the fluctuation window around Tr:
vvT
rT
r
2
1
=
This instantly bring the T = Tr !!
Then scalescale all velocities:
This phase of the simulation is should not be used for averaging!
However to sustain the temp. around Tr we will need
to do this procedure several times at intervals.
![Page 31: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/31.jpg)
Calculating thermodynamic quantities
Average Pressure,
Average Energy, ∑=
>=<M
i
iUM
U1
1
Heat Capacity, Cv:
![Page 32: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/32.jpg)
Reduced Units
• Length (m) Å
• Mass (kg) a m u
• Time (s) Energy eV
• Charge (C) e-
• Temperature (K) K• Temperature (K) K
• Mole (mol) mol
• Luminous Intensity (candela)
![Page 33: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/33.jpg)
• Assignment a)
Title: Generation of FCC lattice
Make 4x4x4 unit cells of the FCC lattice of Argon. Unit cell
HandHandHandHandHandHandHandHand--------OnOnOnOnOnOnOnOn--------Session: 1 (Today! 2Session: 1 (Today! 2Session: 1 (Today! 2Session: 1 (Today! 2Session: 1 (Today! 2Session: 1 (Today! 2Session: 1 (Today! 2Session: 1 (Today! 2--------4 pm) 4 pm) 4 pm) 4 pm) 4 pm) 4 pm) 4 pm) 4 pm)
Name your codes as "loginid_xy.c".
Name the output files of the code as "loginid_xy.out".
Have your name, name, name, name, rollnorollnorollnorollno, date of the session, and title of the assignment , date of the session, and title of the assignment , date of the session, and title of the assignment , date of the session, and title of the assignment
as the first four lines of your code.
You may bring with you programming books.
Make 4x4x4 unit cells of the FCC lattice of Argon. Unit cell parameters 5.26 A (angstrom).
• Assignment b)
Title: Calculation of total potential energy
Argon atoms interact through Lennard-Jones potential with epsilon = 119.8 kB (kB - Boltzmann constant). Convert epsilon to eV before using it in the code (electron Volts). sigma = 3.4 angstrom.
Read the positions of the Ar atoms from the output of a) and calculate the total potential energy of the system (without any periodic boundary conditions or cut of distances imposed).
![Page 34: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/34.jpg)
Lab session -I!Generation of coordinates of FCC lattice (4x4x4 unit cells)
program fcc
implicit none
integer::i,j,k, n =4 ! No of unit cells
real::a=5.26
open(unit=1,file="fcc.dat")
do i = 0, n - 1
do j = 0, n - 1 do j = 0, n - 1
do k = 0, n - 1
write(1,*) i*a + 0.0, j*a + 0.0, k*a + 0.0
write(1,*) i*a + a/2, j*a + a/2, k*a + 0.0
write(1,*) i*a + 0.0, j*a + a/2, k*a + a/2
write(1,*) i*a + a/2, j*a + 0.0, k*a + a/2
end do
end do
end do
end program fcc
![Page 35: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/35.jpg)
allocate(x(no_atom),y(no_atom),z(no_atom))
sigma6 = sigma**6; sigma12 = sigma6**2; eps4 = 4* eps;
pot_en =0.0
do i=1,no_atom-1
do j=i+1,no_atom
dx = x(i) - x(j) ; dy = y(i) - y(j); dz= z(i)-z(j);
!Potential Energy of fcc lattice (4x4x4 unit cells )
r2=dx**2+dy**2+dz**2
r6= r2*r2*r2
r12= r6 * r6
fact= eps4*(sigma12/r12 - sigma6/r6)
pot_en = pot_en + fact
end do
end do
print*,"potential energy=",pot_en
![Page 36: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/36.jpg)
Integrating Equations of Motion
Verlet Scheme:
Newton’s equations are time reversible,
Summing the two equations,Summing the two equations,
Now we have advanced our atoms to time t+∆t ! !
Velocity of the atoms:
![Page 37: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/37.jpg)
Controlling the Temperature
Velocity scaling:
∑=N
iii
vmNk
T2
31Actual temp. at some instant.
TTTTTrr
∆+><∆−If T is out side the fluctuation window around Tr:
vvT
rT
r
2
1
=Then scalescale all velocities:
Though Good and simple algorithm,
Verlet Algorithm does not involve velocities explicitly!
- Not convenient in implementing velocity scaling!
![Page 38: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/38.jpg)
30-1-12
![Page 39: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/39.jpg)
“Velocity Verlet” Scheme
Calculate the positions at
Advance velocities by half-time step,
As good as the Verlet scheme in terms of numerical quality – energy conservation, storage (3N –words) etc., but has the advantage of having velocities explicitly in the equations.
Complete the velocity-move,
Calculate Forces
![Page 40: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/40.jpg)
Verlet scheme
Velocity-Verlet scheme
CALL FORCE
![Page 41: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/41.jpg)
Structural Characterization
Radial Distribution Function (RDF)/ g(r)
Rdf of fcc-solidsolid
![Page 42: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/42.jpg)
Structural Characterization
Running Coordination Number C(r)
C(r)
Rdf of fcc-solidsolid
![Page 43: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/43.jpg)
Calculation of RDF and RCN
••
••
••
•
••• ii
•
•jj’’
L••
••
•jj
• L••dx = x(j) – x(i)dy = y(j) – y(i)dz = z(j) – z(i)
dx = dx – L*ANINT(dx/L)dy = dy – L*ANINT(dy/L)dz = dz – L*ANINT(dz/L)
r2=dx**2+dy**2+dz**2
The same PBC as that used for MD should bebe used in g(r) calculation!
![Page 44: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/44.jpg)
rij
k = INT(rij/binsize)
rdf (k) = rdf (k) + 1
binsize ∼ 0.02 - 0.1
![Page 45: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/45.jpg)
Complete the assignments of Lab Session-I.
(Demonstrate running the code and the Pot. Energy obtained)
Make a new copy of your Pot. Energy Program &
Name it “loginid_2a.c” (retain your older version!!)
1. Read the fcc (4x4x4) positions (x, y, z -already generated)
Hands-On-Session –II (23/1/12) Time 2-4pm
1. Read the fcc (4x4x4) positions (x, y, z -already generated)
2. Impose PBC & Minimum Image Convention
3. Calculate Pot. Energy U for different cutoff distances (Rc);
(Find U for Rc=4.0 to 10 Å in steps of 0.5 Å )
Bring Your Notebooks/C/Fortran- books;
Do not use internet!
![Page 46: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/46.jpg)
FCC lattice + PBC
U
1. Steps in U coincide with the g(r)! 2. U converges to as r increases.
g(r)
a/√2a
a√3/2
![Page 47: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/47.jpg)
![Page 48: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/48.jpg)
Initialization of Velocities (Random)
1 322 2
.N
i iiK E m v N k T∑= < > =
∑=N
iii
vmNk
T2
31
Uniform random No (0-1)Uniform random No (0-1)
![Page 49: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/49.jpg)
30-1-12
![Page 50: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/50.jpg)
Structural Characterization
Radial Distribution Function (RDF)/ g(r)
Rdf of fcc-solidsolid
![Page 51: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/51.jpg)
“Tail” Correction to Potential Energy!
N(
assuming
NN
N
N
![Page 52: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/52.jpg)
!cmvx,cmvy,cmvz=0; sumvsq=0; rand(3)
do i=1,nCall Random_Number(rand)vx(i) = (rand( 1) -0.5);vy(i) = (rand( 2) -0.5);vz(i) = (rand( 3) -0.5);
Initialization of Velocities (Random)
cmvx = cmvx + vx(i); cmvy = cmvy + vy(i); !For Centre of mass velocitycmvz = cmvz + vz(i); sumvsq = sumvsq + vx(i)**2 + vy(i)**2 + vz(i)**2enddo
cmvx = cmvx/n; cmvy = cmvy/n; cmvz = cmvz/n;
! cmvx~0 cmvy ~ 0 cmvz ~0
![Page 53: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/53.jpg)
do i=1,nvx(i) = vx(i) –cmvx;vy(i) = vy(i) –cmvy;vz(i) = vz(i) –cmvz;
enddo
Then scalescale all velocities: vv rT 2
1
=
3 1 22 2
N
i iiN kT m v∑= 2
3
Nm
N k iiT v∑=
Then scalescale all velocities: vvT
r
r
=
do i=1,nvx(i) = vx(i)*fact;vy(i) = vy(i)*fact;vz(i) = vz(i)*fact;
enddo
fact = SQRT(3*n*kB* Tr / (m*sumvsq))
![Page 54: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/54.jpg)
![Page 55: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/55.jpg)
Dynamical Properties: Diffusion Coefficient
Fick’s Law:
Continuity Eq.:
Diffusion Eq.:
Einstein’s relation: Nr. MSD
TfkDNq B/2=σNernst-Einstein’s relation:
![Page 56: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/56.jpg)
06 Feb 2012
![Page 57: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/57.jpg)
Dynamical Properties: Diffusion Coefficient
Fick’s Law: -
- concentration
D – self diffusivity
(diff. coefficient)
Continuity Eq.:
Diffusion Eq.:
(diff. coefficient)
![Page 58: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/58.jpg)
Prob. Interpretation of c(r,t)of c(r,t)
(Einstein’s relation)
![Page 59: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/59.jpg)
Frenkel & Smith
![Page 60: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/60.jpg)
Einstein’s relation: Nr. M.S.D
(M. S. D – Mean Squared Displacement)
TfkDNq B/2=σNernst-Einstein’s relation:
(M. S. D – Mean Squared Displacement)
![Page 61: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/61.jpg)
CaF2
MS
D
…Diffusion Coefficient: CaF2
F-
Ca2+
time
TfkDNq B/2=σ
PRB, 43, 3180,1991.
![Page 62: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/62.jpg)
In 1-Dimension
- Velocity auto correlation fn.
![Page 63: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/63.jpg)
…Dynamical Properties: Vibrational Spectrum
Velocity Autocorrelation Spectrum, Water @300 K
Power Spectrum,
Helps interpreting IR spectrum.
![Page 64: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/64.jpg)
Demonstrate the lab assignments of sessions-I, II, III.
Make a new copy of your Force Calculation Program &
Name it “loginid_4a.c” (retain your older version!!)
1. Read the fcc (4x4x4) positions (x, y, z -already generated) Keep (0, 0, 0) particle displaced at (-0.2, 0.3, 0.5).
2. Impose PBC & Minimum Image Convention (Rcut = 10A)
Hands-On-Session –IV (13/2/12) Time 2-4pm
2. Impose PBC & Minimum Image Convention (Rcut = 10A)
3. Use the forces to integrate the equations of motion using“velocity Verlet” algorithm. Use time step, dt = 1 fs.
mass- that of argon. (Workout the unit system!!!)
Perform 100 MD steps.
Bring Your Notebooks/C/Fortran- books;
Do not use internet!
![Page 65: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/65.jpg)
Reduced Units
• Length (m) Å
• Mass (kg) a m u
• Time (s) / Energy eV
• Charge (C) e-
• Temperature (K) K• Temperature (K) K
• Mole (mol) mol
• Luminous Intensity (candela)
time step = 1 fs = ???? (In reduced units)
![Page 66: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/66.jpg)
Mid-Sem: 22nd Feb 2012 (10-12pm)
Bring Your Calculators.
Questions are based on NVE-MD of bulk Systems.
Total marks: 30 (6 x 5 marks)
![Page 67: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/67.jpg)
Interatomic forces for simple systems
1. Lennard-Jones Potential:
Gives an accurate description of inert gasesinert gases
(Ar, Xe, Kr etc.)
(non-bonded interactions)
Instantaneous dipoles
2. Born-Mayer (Tosi-Fumi) Potential:
(Ar, Xe, Kr etc.)
Faithful in describing purepure ionic solids (NaCl, KCl, NaBr etc.)
![Page 68: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/68.jpg)
Issues with Coulombic Forces
(Ewald Summation technique)
Minimum Image Convention
This sum is only Conditionally Convergent!
![Page 69: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/69.jpg)
“Tail” Correction to Potential Energy!
N(
assuming
NN
N
N
![Page 70: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/70.jpg)
![Page 71: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/71.jpg)
Remarks on Interatomic Forces
• Development of good force fields (FF) can be a tough task!
• FF assume that electronic clouds around the nucleus of atoms is intact irrespective of the environment around the atom!
FF’s are developed by empirical methods or ab-initio calculations.
This can be a poor assumption for highly polarizable atoms/ions!
Solution?
Develop a shell model of atoms/ions!
Or DFT-based ab-initio (Car-Parrinello) MD calculations !
![Page 72: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/72.jpg)
Comments on Classical MD
Extensively employed to understand Physical processes at atomic resolution
Phase Transitions,
Diffusion and transport properties,
Local structural and short-time relaxation of
� crystalline and amorphous solids
Very powerful in studying a variety of physical phenomena and under several external conditions (T & P).
� crystalline and amorphous solids
� liquids
�solid-fluid interfaces
�nano-clusters
And, serves a very useful bridge between experiment and theory!
Not useful in the study of electronic properties!
Not powerful enough to describe chemical reactions!
![Page 73: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/73.jpg)
Complete the lab assignments of sessions-IV.
(Perform 1000 MD steps; with the particle positions, Potential
Energy and Total Energy written at intervals of 10 MD steps;
and demonstrate energy conservation.)
Use the positions to calculate the Radial Distribution Function (averaged over all particles and over 100 MD steps.)
Hands-On-Session –V (05/03/12) Time 2-4pm
(averaged over all particles and over 100 MD steps.)
Compare the same with the RDF of the static-FCC lattice.
Bring Your Notebooks/C/Fortran- books;
Do not use internet!
![Page 74: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/74.jpg)
Translational & Rotational corrections (from last class)
3 –translational and 3- overall-rotational DOF are to be subtracted. be subtracted.
![Page 75: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/75.jpg)
Ar clusters: Structure
IcosahedronAr13 Cluster
Typically, clusters stabilize in different structures
than their bulk counterparts.
![Page 76: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/76.jpg)
Ar cluster Magic Numbers
Ar13 Ar55 Ar147 Ar309
Stability depends on completion of icosahedral shells.
![Page 77: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/77.jpg)
Magic numbers of Na clusters
l < n
![Page 78: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/78.jpg)
Stability of Na clusters?
Chemical Hardness, I- Ionization potential; A-electron affinity
![Page 79: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/79.jpg)
Experimental evidence: Mass spectrometric dataFor Na cluster
![Page 80: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/80.jpg)
Ar clusters: Structure
IcosahedronAr13 Cluster
Typically, clusters stabilize in different structures
than their bulk counterparts.
![Page 81: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/81.jpg)
Ar cluster Magic Numbers
Ar13 Ar55 Ar147 Ar309
Stability depends on completion of icosahedral shells.
![Page 82: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/82.jpg)
Experimental evidence for magic clusters: Mass spectrometric data for Inert gas clusters
![Page 83: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/83.jpg)
Solid – Liquid transition of Ar clusters
RMS bond length fluctuation,
Ar
(Lindemann Criterion)
Ar13
![Page 84: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/84.jpg)
Solid – Liquid transition of Ar clusters
Heat Capacity Vs Temperature
![Page 85: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/85.jpg)
Solid-Liquid coexistence (25-36 K)
Liquid-like
solid-like
![Page 86: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/86.jpg)
Solid – Liquid transition of Ar clusters
Temp. Vs Tot. Energy Curves!
Tm
Tf
![Page 87: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/87.jpg)
Bimodal distribution of U
![Page 88: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/88.jpg)
Structure of (KCl)32
Born-mayer potential
![Page 89: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/89.jpg)
![Page 90: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/90.jpg)
![Page 91: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/91.jpg)
Quenching &
Steepest Descend - methods
These are all stablestable (solid) structures corresponding to “local minima” of the multi-dimensional Pot. Energy surface.
![Page 92: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/92.jpg)
RDFs of Quenched Structures
1st excited state
ground-state structure
![Page 93: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/93.jpg)
“non-wetted” structure
![Page 94: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/94.jpg)
Amorphous/glassy phase
![Page 95: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/95.jpg)
Statistics of quenched structures crystalline
non-wetted
amorphous
![Page 96: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/96.jpg)
RMS bond length fluctuation,
![Page 97: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/97.jpg)
![Page 98: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/98.jpg)
Velocity Autocorrelation Spectrum,
Power Spectrum,Power Spectrum,
![Page 99: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/99.jpg)
![Page 100: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/100.jpg)
Solid-Liquid Transition: Ar147 vs (KCl)32
T= 38 K“Surface melting”
Inner core remains solid.
“wetting” phenomenon !“wetting” phenomenon !(in Ar147
)
![Page 101: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/101.jpg)
Ar cluster Magic Numbers
Ar13 Ar55 Ar147 Ar309
Stability depends on completion of icosahedral shells.
![Page 102: Syllabus - iitg.ac.in](https://reader030.vdocument.in/reader030/viewer/2022012505/61805be8d15fb6431b0b2b9b/html5/thumbnails/102.jpg)
Solid-Liquid Transition in (KCl)32
�No evidence for a “dynamical” solid-liquid co-existence (like in Ar13) is observed. But still a possibility!
� A crystal-like form co-existing simultaneously with a disordered structure is observed near the transition.
� However this transition is “non-wetting” in nature.((the melted portion does not cover the solid portion uniformly).