-
Notifications
You must be signed in to change notification settings - Fork 42
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
QCG is having issue with adding a large number of solvent molecules #202
Comments
Hi, |
Hi,
Thank you very much for your help. I've attached the input coordinates.
czw., 1 cze 2023 o 13:06 Christoph Plett ***@***.***>
napisał(a):
… Hi,
This might be indeed a problem of the wall potential that becomes too
small.
If you can send me your input coordinates, please, then I can check on
this.
—
Reply to this email directly, view it on GitHub
<#202 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/BAHUK6AH5JIK2RDDTGPVYLDXJBZTBANCNFSM6AAAAAAYWU5LQA>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
Hi, I'm unsure whether my input coordinates were sent alongside my message. So, I've attached the coordinates below. Thank you. 71 N 4.7999300130 1.5396351138 -0.4639876729 |
Hi, Thank you for sending the results. According to the visualisation of the system, I guess that the NADH molecule was not constrained during the adding water molecules. The final NADH structure surrounded by 200 water molecules is quite far from the starting structure, which I must maintain during the cluster growth. To constrain the NADH molecule, I added '--cinp constrain' to the command line. Could it be a problem for growing the cluster? The contents of the constrain file: I use the release versions of a crest (2.12), xtb (6.6.0), and xtbiff (1.1) downloaded from the github. According to my understanding, I've been using only the crest with the xtbiff algorithm. Is there a way to force using the new aISS algorithm of xtb to grow the cluster? Thank you in advance for your help. |
Hi, |
@cplett actually the continous release should include the aISS usage capability, shouldn't it? |
Somehow, the continuous release didn't update after the last commits. |
According to the description the release has the status of commit #189 of April 3rd, so it should be up to date. You can clearly see the binaries were uploaded on that date. The release tag never changes from Dec 23, 2022. |
Thank you for all your suggestions. I tested the aISS algorithm for the cluster growth from the continuous release, and it has the same problem as the xtbiff. When I freeze the solute (using the .xcontrol file and adding the --fixsolute and --nopreopt flags), I cannot add more than 94 water molecules to the system. Formally, the NADH molecule has a charge of -2. Manipulating the charge, I wanted to check whether the negatively charged system could affect the cluster growth, but it's not. Any advice would be appreciated. |
Hi, |
Hi, I've attached the whole printout of crest. Thank you. |
Hi, |
Hi, Thank you for the suggestions. Increasing the value of the OMP_STACKSIZE variable helped. I tested both the aISS and xtbiff approaches, and everything works. |
This issue had no activity for 6 months. It will be closed in 1 week unless there is some new activity. |
Hello,
I'm using the QCQ extension of CREST to create a water cluster (200 water molecules) around frozen stacked chromophores. The problem appears during the addition of the 93rd water molecule to the system. Below, I've attached the last lines of the output file. I was trying to change various parameters such as --wscal, -mdtime, --mddump, but the problems occurred exactly at the same step. I would like kindly ask to any help, suggestion or advice.
92 -615.535733 -16.42 -1229.06 1.362 -32.902 14.2 21.3 19243.3 normal Cluster energy not found Total Energy of cluster not found. 93 0.000000 ******* ******** 1.357 0.000 14.3 20.2 19459.7 normal
The text was updated successfully, but these errors were encountered: