Bond Breaking in MD

From NWChem

You are viewing a single post from the thread title above
Jump to: navigation, search

Click here for full thread
Clicked A Few Times
Threads 10
Posts 16
I am running a molecular dynamics simulation on a small molecule for which reported force field does not include bond stretching parameters. When I use just these reported Lenard-Johns parameters, bonds break and I do not have a true molecule. When I add excess parameters for bond stretching, Although is not a correct manner, I get a positive value for total energy that resembles my system is not stable. Also I tried to define these bond lengths as z-matrixes in segment file according to some segment files in the test folder, but I got " Fortran runtime error: End of file". Any help is appreciated.

Best Regards
Saeed


Who's here now Members 0 Guests 0 Bots/Crawler 0


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC