Author Topic: Question about Boron diffusion in bulk silicon tutorial  (Read 354 times)

0 Members and 1 Guest are viewing this topic.

Offline pshinyeong

  • Regular QuantumATK user
  • **
  • Posts: 6
  • Country: kr
  • Reputation: 0
    • View Profile

I was looking at tutorial titled "Boron diffusion in bulk silicon" (link = "" ) to perform Adaptive Kinetic Monte Carlo (AKMC) simulation. I realized that silicon was initially geometrically optimized but with the two atoms (which is all the atoms there is in the primitive cell) fixed. Does performing a geometry optimization calculation with all the atoms fixed have any meaningful purpose? I didn't understand why the Silicon atoms were fixed here.


+++ This was a bad question, lattice vectors were unfixed. Problem solved.
« Last Edit: August 6, 2023, 07:49 by pshinyeong »

Offline Anders Blom

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 5317
  • Country: dk
  • Reputation: 88
    • View Profile
    • QuantumATK at Synopsys
Re: Question about Boron diffusion in bulk silicon tutorial
« Reply #1 on: August 7, 2023, 19:54 »
This old example is quite simplistic and was written before we had more fancy constraint. The purpose of the setup as written was indeed to optimize the cell size (lattice vectors) without disrupting the space group symmetry. Today you can do that with a space group constrain, but actually a better way to do that for a cubic material is to run a quick equation of state curve and find the lattice constant giving the minimum energy.