The problems I have been having with NMR seem to be related to the basis sets. I was using the pcSseg basis sets, which seems to work fine for molecules, but was getting huge unphysical chemical shifts for periodic crystals (*****). When I switched to EMSL 6-31G** I am getting reasonable numbers.<div><br></div><div>Ron</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, July 25, 2022 at 4:31:19 PM UTC-4 Ronald Cohen wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">By the way, the regtests for nmr pass, but I do not see any with non-orthogonal lattices.<div>Ron</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Sunday, July 24, 2022 at 9:53:18 PM UTC-4 Ronald Cohen wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I tried another large nonorthogonal system of similar chemistry and find the same thing The shifts are all printed as ****.<div><br></div><div>Ron</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, July 22, 2022 at 3:19:16 PM UTC-4 Ronald Cohen wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thanks! I changed the GAUGE and it now runs! However, the results are no good!<div><br></div><div>Shielding atom at atomic positions. # tensors printed 165 <br> 1O O 2.390434 2.226322 0.031347<br> SIGMA from SOFT J<br> XX = ********** XY = ********** XZ = **********<br> YX = ********** YY = ********** YZ = **********<br> ZX = ********** ZY = ********** ZZ = **********<br> SIGMA from LOCAL J<br> XX = -5315.3694 XY = ********** XZ = **********<br> YX = 98.5820 YY = ********** YZ = **********<br> ZX = 1884.0679 ZY = ********** ZZ = **********<br> SIGMA TOTAL<br> XX = ********** XY = ********** XZ = **********<br> YX = ********** YY = ********** YZ = **********<br> ZX = ********** ZY = ********** ZZ = **********<br> ISOTROPY = ************ ANISOTROPY = ************<br><div><br></div><div><br></div><div>Do you have any suggestions? I simply had changed to</div><div>GAUGE R_AND_STEP_FUNCTION<br></div><div><br></div><div>I let the system restart. Attached is the new out file.</div><div><br></div><div>I also note that it did not try to go to the next step in the trajectory with </div><div>MOTION<br> &MD<br> &REFTRAJ<br> FIRST_SNAPSHOT 1000<br> STRIDE 1000<br> TRAJ_FILE_NAME MD-pos.dcd<br> VARIABLE_VOLUME .TRUE.<br> &END REFTRAJ<br> &END MD<br>&END MOTION<br></div><div><br></div><div>Should that work with LINRES? There are no errors--it is like it didn't read it. </div><div><br></div><div>Also how do I read the checksums in the NMR calculation?</div><div><br></div><div>Thank you!</div><div><br></div><div>Sincerely,</div><div><br></div><div>Ron</div><div><br></div><div><br></div><div><br><br></div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, July 22, 2022 at 10:32:18 AM UTC-4 Marcella Iannuzzi wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Ron Cohen<div><br></div><div>The problem seems to be that for the GAUGE ATOM not-orthorhombic boxes are not accepted.</div><div>Unfortunately there is no earlier stop. I put this in the to-do list</div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, July 22, 2022 at 3:59:03 PM UTC+2 <a rel="nofollow">reco...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I am trying to run a large supercell for nmr with the current development CP2K and after many failures got the linear response to converge. However, it dies from a bound error in qs_linres_current.F:1321 . The last printout is:<div><br><div><div>Inizialization of the NMR environment</div><div><br></div><div> NMR| Shift gapw radius (a.u.) 1.133836E+02</div><div> NMR| Shift factor (ppm) 6.601835E-03</div><div> NMR| Shift factor gapw (ppm) 5.325134E+01</div><div> NMR| Chi factor (SI) 1.972757E+01</div><div> NMR| Conversion Chi (ppm/cgs) 6.022045E-02</div><div> NMR| Conversion Chi to Shift 7.008891E-04</div><div> NMR| Shielding tensor computed for 165 atoms</div><div><br></div><div> *******************************************************************************</div><div> * ___ *</div><div> * / \ *</div><div> * [ABORT] *</div><div> * \___/ CPASSERT failed *</div><div> * | *</div><div> * O/| *</div><div> * /| | *</div><div> * / \ qs_linres_current.F:1321 *</div><div> *******************************************************************************</div><div><br></div><div><br></div><div> ===== Routine Calling Stack ===== </div><div><br></div><div> 5 calculate_jrho_resp</div><div> 4 current_build_current</div><div> 3 linres_calculation_low</div><div> 2 linres_calculation</div><div> 1 CP2K</div></div><div><br></div><div>I attach my inputs and output. Thank you for any ideas or help. Is there something I need to increase, or do I need to choose a different gauge?</div><div><br></div><div>Sincerely,</div><div><br></div><div>Ron Cohen</div><div><br></div></div></blockquote></div></blockquote></div></blockquote></div></blockquote></div></blockquote></div>
<p></p>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d45a8e53-d73b-4886-bb9d-b75f769e14f6n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d45a8e53-d73b-4886-bb9d-b75f769e14f6n%40googlegroups.com</a>.<br />