Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Symmetry in Core-Level Spectroscopy #4

Open
ajm143 opened this issue Sep 12, 2019 · 5 comments
Open

Symmetry in Core-Level Spectroscopy #4

ajm143 opened this issue Sep 12, 2019 · 5 comments
Assignees
Labels
bug Something isn't working
Milestone

Comments

@ajm143
Copy link
Contributor

ajm143 commented Sep 12, 2019

No description provided.

@ajm143 ajm143 added the bug Something isn't working label Sep 12, 2019
@ajm143 ajm143 added this to the v2.0 milestone Sep 12, 2019
@Rebecca-Nicholls
Copy link
Contributor

If you do a core-level spectroscopy calculation using symmetry in the castep calculation, you get the wrong answer. To get the correct answer, turn off symmetry in the castep calculation.

@Rebecca-Nicholls
Copy link
Contributor

Rebecca-Nicholls commented Dec 16, 2021

image
image
Calculations done with castep 20.11
Problem exists in latest version of OptaDOS
Problem also with castep 20.11 and elnes (rather than spectral) calculation
Cu2O.zip

@ajm143
Copy link
Contributor Author

ajm143 commented Feb 10, 2022

@ajm143 to make a temporary io_err statement. if you have more than one symmetry element -- stop.

I think adding this to line 104 of core.f90 would do it:

if(num_sym > 1) call io_error('Error: Core loss not currently able to deal with symmetry. Please &
  & re-run CASTEP without symmetry.')

@ajm143
Copy link
Contributor Author

ajm143 commented Dec 5, 2022

N.B. The current test-suite has a set of core examples with symmetry enabled -- so if you make this change the test-suite worn't work.

@ajm143
Copy link
Contributor Author

ajm143 commented Dec 7, 2022

Commit 581c39c now has this temporary io_err statement and the tests have been updated accordingly

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

3 participants