PSPW QM/MM Lennard-Jones Ion Parameters

From NWChem

Viewed 1984 times, With a total of 3 Posts
Jump to: navigation, search

Just Got Here
Threads 2
Posts 4

I'm a new NWchem user trying to run a decently large QM/MM simulation. My QM region consists of ~150 atoms, and my MM region will be ~500 ethanol molecules. I'd like to apply this procedure to larger systems in the future, so I figured I would use the QM/MM capabilities in the PSPW module. The documentation has a very instructive example (, however I'm a bit confused about how to translate this procedure to my molecule.

More specifically, I've no idea where the Lennard-Jones ion parameters came from, and thus I've no clue where to obtain the relevant values for my molecule. Additionally, it's unclear to me if two carbon atoms in different chemical environments would have different LJ parameters, or if the LJ parameters of atoms in the QM region change as a function of solvent molecules in the MM region.

I'm aware of the more standard QM/MM Gaussian-based module and am working to get that going for my systems, too, but any help with the PSPW module and these ion parameters would be greatly appreciated.


Edited On 1:40:10 PM PDT - Tue, Aug 27th 2013 by Sblau

Just Got Here
Threads 2
Posts 4
Answering my own question - LJ parameters can be found throughout the literature, especially from papers about the OPSL-AA forcefield. If anyone else is struggling to find them, feel free to respond and I can send over some helpful papers.

Clicked A Few Times
Threads 10
Posts 16
Dear All

We are a small group starting to play MD with nwchem. Unfortunately we are new with MD. Hereby I would like to have the form of Lennard-Jones equation from which parameter files are written in nwchem. Also, symbols for definition of TiO2 bulck and surface atoms in .par and .frg files. Also, the form of definition of an organic liquid which is adsorbed on the TiO2: number of organic molecules in the PDB and FRG files. I am searching for several weeks in nwchem forum and other places, but I can not find any help except to the QMMM appendix page which is not useful for me. 

Best Regards

Just Got Here
Threads 1
Posts 4
Dear Sam:
   Have you solved this problem? I also encountered the same problem. Two carbon atoms in different chemical environments in QM region would have different LJ parameters. Can you tell me how to do? Can you tell me your email address? My email address is
Thank you!

Forum >> NWChem's corner >> QMMM

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

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