Core-shell model and QM/MM ECOUPL GAUSS are not compatible?

zuohen... at gmail.com zuohen... at gmail.com
Wed Aug 22 12:58:50 UTC 2018


Hi Matt,

Thank you for your work. I am not sure if the fault is caused by my 
incorrect set.

Best wishes,
Zuoheng

在 2018年8月21日星期二 UTC+2下午8:32:12,Matt W写道:
>
> Thanks for the report, I'll try and have a look.
>
> Core-shell QM-MM should be working, but has never really been tested.
>
> Matt
>
> On Tuesday, August 21, 2018 at 4:42:39 PM UTC+1, zuoh... at gmail.com 
> wrote:
>>
>> Thank you for your response.
>> I am using the EWALD_TYPE SPME for MM. It is OK for MM force field. 
>> However, I want to combine QM and MM and I have to specifies the type of 
>> the QM - MM electrostatic coupling (ECOUPL). This is different from 
>> EWALD_TYPE.
>>
>> Best wishes,
>> Zuoheng
>>
>>
>> 在 2018年8月21日星期二 UTC+2下午4:21:35,Sergey Chulkov写道:
>>>
>>> Thank you for sharing your experience.
>>>
>>> Core shell model is currently not implemented for the ewald sums. Try to 
>>> use a different EWALD_TYPE, for example, SPME.
>>>
>>> Best wishes,
>>>
>>> Sergey
>>>
>>> On Tuesday, August 21, 2018 at 1:39:19 PM UTC+1, zuoh... at gmail.com 
>>> wrote:
>>>>
>>>> Dear CP2Kers,
>>>>
>>>> I am using core-shell model in MM. I found a Segmentation fault during 
>>>> QM/MM simulation with ECOUPL GAUSS. 
>>>> I am not sure if it is a bug? or do I have to use ECOUPL NONE? (I am 
>>>> using the CP2K of 6.1 version)
>>>> The attached file is the input file I used to test this fault. This is 
>>>> changed from the test file of CP2K.
>>>>
>>>> The fault is following:
>>>>
>>>> Program received signal SIGSEGV: Segmentation fault - invalid memory 
>>>> reference.
>>>>
>>>> Backtrace for this error:
>>>> #0  0x7FDD41DCCE08
>>>> #1  0x7FDD41DCBF90
>>>> #2  0x7FDD412C44AF
>>>> #3  0x13180AF in __qmmm_gpw_energy_MOD_qmmm_elec_with_gaussian_lr at 
>>>> qmmm_gpw_energy.F:800
>>>> #4  0x131A147 in __qmmm_gpw_energy_MOD_qmmm_elec_with_gaussian_low at 
>>>> qmmm_gpw_energy.F:503
>>>> #5  0x131A645 in __qmmm_gpw_energy_MOD_qmmm_elec_with_gaussian at 
>>>> qmmm_gpw_energy.F:310
>>>> #6  0x131B643 in __qmmm_gpw_energy_MOD_qmmm_el_coupling at 
>>>> qmmm_gpw_energy.F:197
>>>> #7  0x9E3635 in __qmmm_force_MOD_qmmm_calc_energy_force at 
>>>> qmmm_force.F:160
>>>> #8  0x8015AF in __force_env_methods_MOD_force_env_calc_energy_force at 
>>>> force_env_methods.F:251
>>>> #9  0x5CF918 in __md_run_MOD_qs_mol_dyn_low at md_run.F:371
>>>> #10  0x5D095E in __md_run_MOD_qs_mol_dyn at md_run.F:152
>>>> #11  0x4EF51E in __cp2k_runs_MOD_cp2k_run at cp2k_runs.F:356
>>>> #12  0x4F1D6B in __cp2k_runs_MOD_run_input at cp2k_runs.F:1103
>>>> #13  0x4DBB21 in MAIN__ at cp2k.F:281
>>>>
>>>>
>>>> ===================================================================================
>>>> =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
>>>> =   PID 13646 RUNNING AT ubuntu
>>>> =   EXIT CODE: 139
>>>> =   CLEANING UP REMAINING PROCESSES
>>>> =   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
>>>>
>>>> ===================================================================================
>>>> YOUR APPLICATION TERMINATED WITH THE EXIT STRING: Segmentation fault 
>>>> (signal 11)
>>>> This typically refers to a problem with your application.
>>>> Please see the FAQ page for debugging suggestions
>>>>
>>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180822/82d63da6/attachment.htm>


More information about the CP2K-user mailing list