Holding In-plane lattice constants for heterostructures

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Locked
Message
Author
rightsuman
Newbie
Newbie
Posts: 5
Joined: Mon Nov 18, 2019 8:25 pm

Holding In-plane lattice constants for heterostructures

#1 Post by rightsuman » Thu May 26, 2022 4:48 am

If we have an interface, say LaAlO3/SrTiO3, often, it is mentioned that we hold the in-plane lattice constants and allow the growth direction to relax. Does this imply selective dynamics with ISIF 2 or something else?

martin.schlipf
Global Moderator
Global Moderator
Posts: 455
Joined: Fri Nov 08, 2019 7:18 am

Re: Holding In-plane lattice constants for heterostructures

#2 Post by martin.schlipf » Thu May 26, 2022 5:20 am

If you have a vacuum region then that is exactly what you need to do.
If it is a repeated layer structure you need to optimize the z direction as well. Then you could setup multiple calculations with different heights, relax the positions and the lowest energy structure is the most stable one.

User avatar
chengcheng_xiao1
Newbie
Newbie
Posts: 26
Joined: Sun Nov 17, 2019 6:23 pm
Location: London, UK
Contact:

Re: Holding In-plane lattice constants for heterostructures

#3 Post by chengcheng_xiao1 » Thu May 26, 2022 12:51 pm

If I didn't get it wrong and you are trying to relax only one lattice vector while keeping the other two fixed, maybe give this a try? https://github.com/Chengcheng-Xiao/VASP_OPT_AXIS

rightsuman
Newbie
Newbie
Posts: 5
Joined: Mon Nov 18, 2019 8:25 pm

Re: Holding In-plane lattice constants for heterostructures

#4 Post by rightsuman » Mon Jun 06, 2022 6:34 am

I was trying to benchmark a SrTiO3-LaAlO3 system of 9 (SrO + TiO2) and 5 LaAlO3 (LaO + AlO2) layers, but I am not getting the required DOS at the interface whereas the surface DOS seems to be very in agreement with what is available in the literature. I have given a 15 A vacuum to prevent interactions and performed with as well as without dipole corrections. Is there something I am missing while making the structure or performing normal SCF calculations and thereafter DOS. I am attaching the files for scf here, if anyone can suggest anything.

Note: I have fixed the two layers of SrTiO3 in one of them to give a bulk like effect, which is also mentioned in some of the literatures.

rightsuman
Newbie
Newbie
Posts: 5
Joined: Mon Nov 18, 2019 8:25 pm

Re: Holding In-plane lattice constants for heterostructures

#5 Post by rightsuman » Mon Jun 06, 2022 7:08 am

I have attached it here, missed it in the previous reply.
You do not have the required permissions to view the files attached to this post.

martin.schlipf
Global Moderator
Global Moderator
Posts: 455
Joined: Fri Nov 08, 2019 7:18 am

Re: Holding In-plane lattice constants for heterostructures

#6 Post by martin.schlipf » Tue Jun 07, 2022 10:37 am

Sorry, but your input is not specific enough. You mention that the DOS looks wrong, but you do not provide the input/output of that calculation. Also could you clarify better what you see and what you expect, please? These calculations with very long cell are very sensitive, because there may be local minima that you can run into.
You may also want to plot the potential (LVTOT) and make sure that it looks like what you expect. In particular, make sure that the potential has a kink in the vacuum region.

Locked