Symmetry of XANES L-edge for Si defect in Graphene

Problems running VASP: crashes, internal errors, "wrong" results.

Moderators: Global Moderator, Moderator

Message
Author
aohara
Newbie
Posts: 10
Joined: Sat Apr 16, 2016 5:54 am

Symmetry of XANES L-edge for Si defect in Graphene

Dear VASP team and community,

I am currently working to use the new XANES functionality in VASP 6.2.1 to calculate core-loss spectra for a project. I am technically interested in ELNES, but they should essentially be the same.

First, however, I am trying to reproduce some of the results in the following paper:
W. Zhou et al, "Direct Determination of the Chemical Bonding of Individual Impurities in Graphene," Phys Rev Lett 109, 206803 (2012).
In particular, I am attempting to calculate the spectra for the defect in which an Si atom fills a carbon divacancy (Si-4C), this would be the L2/3 edge of silicon and the composite theoretical spectra is shown in Fig. 3(b) of the paper. In the SI, the breakdown into the in-plane and out-of-plane contribution from the dielectric function are shown in S2.

Now, I've used a slightly larger than this paper. However, when I calculate the L2/3 edge for the silicon, I obtain a spectra that is not anisotropic. That is the core dielectric tensor has the same value for the xx, yy, and zz at each energy. In the attached plot the two lines are exactly on top of each other - the plot was re-broadened from the raw output, but this is not the issue. I am not entirely sure what the issue here is as I expect the spectra to not be isotropic between the in-plane and out-of-plane components. The peaks obtained via VASP match approximately to those obtained before but their intensities are quite different which suggest to me that some sort of averaging effect is being introduced.
Graphene_Si-4C_L_edge_calculated.png
If instead I calculate the carbon K-edge for pristine graphene (also re-broadened in postprocessing), I obtain an anisotropic tensor as one would expect. I can then average them based on the experimental beam parameters.
Pristine_Graphene_C_Kedge_calculated.png
It is not clear to me if I am missing a setting in the INCAR or if something is happening inside the code though. Any help would be much appreciated. I have included my input files that were used to perform the calculation as well as some images. Please let me know if additional information is needed to help.

Andrew
You do not have the required permissions to view the files attached to this post.

henrique_miranda
Global Moderator
Posts: 255
Joined: Mon Nov 04, 2019 12:41 pm
Contact:

Re: Symmetry of XANES L-edge for Si defect in Graphene

Dear Andrew,

Thanks for the detailed report.
Indeed something seems to be wrong. We will look into it
I have one request and a suggestion:

1. Do you think you could try to reproduce this issue on a system with fewer atoms?
Same defect but a smaller graphene supercell.
This would help us a lot to track down the problem.

2. Could you try to run the calculation with ISYM=-1 and see if the result changes?
(This could be done on the system with fewer atoms)

aohara
Newbie
Posts: 10
Joined: Sat Apr 16, 2016 5:54 am

Re: Symmetry of XANES L-edge for Si defect in Graphene

Henrique,

Thank you for looking into this to help resolve the issue. I apologize for the delay in providing the additional calculation due to the American holiday.

I re-ran the calculation using a 5x5 graphene supercell with the same defect of a silicon atom sitting inside a carbon divacancy (but unrelaxed - the only changes would be light in-plane relaxations). For me, this calculation can run on a single node (48 cores) in under 5 minutes. I also set ISYM = -1 as requested. The spectra below is the result (unlike before, no additional broadening was employed - only the setting in VASP):