CompuServe Thread

#Solvents in conf. analy.

4 messages in this thread
#23801From: John HubbsSep 25, 1992 6:37 PM
Could someone please contrast the effect of merging a large number of solvent molecules with a solute vs. the use of a periodic box when carrying out molecular mechanics or semi-empirical calculations. As near as I can tell, periodic boundary conditions are turned off when using either molecular mechanics or semi-empirical routines within Hyperchem. Thus, what is the function of the periodic box if the calculation is done in vacuuo anyway? The merge option has the distinct advantge of allowing a periodic box type calculation with solvents other than water.
#23868From: Tom Slee [Hypercube]Sep 28, 1992 9:14 AM
John, The periodic boundary conditions are not turned off when doing molecular mechanics calculations within HyperChem. We don't have periodic boundary conditions for semi-empirical methods, but those calculations on a box of any size would take for evern and a day anyway. I think the advantages to the periodic box are: 1. The water molecules are read in from a template file, and have standard charges and geometries that are known and used by many people. 2. The box is equilibrated to a 300K box of water, so the arrangement of molecules is pretty good as a starting point, 3. The periodic boundary conditions stop water molecules from floating away. If you carry out a simulation without the periodic conditions you may see solvent molecules wander away from the rest of the system. Tom
#23915From: John HubbsSep 28, 1992 8:13 PM
Thanks. Your reply seems to verify much of what I suspected. A few semi-random comments in reply follow. When I stated that molecular mechanics turned off the periodic boundary conditions, I really meant that molecular DYNAMICS seem to turn off the periodic boundary conditions. At least that's what I recall from attempting it and remember reading somewhere in one of your three (well written) manuals (which I alas have left in my office). Your state <2. The box is equilibrated to a 300K box of water, so the arrangement of molecules is pretty good as a starting point.> This implies that I could carry out a molecular mechanics calculation in a periodic box until convergence was reached. It should then be possible to allow the geometries (of solute and solvent) to relax using a semi-empirical calculation such as AM-1. I suspect that few molecules would diffuse away under such conditions. You further state <We don't have periodic boundary conditions for semi-empirical methods, but those calculations on a box of any size would take forever and a day anyway> This is a matter of personal opinion (in my personal opinion). Computers are typically turned off at the end of each work day. There is frequently significant availble unused time on most pc's which surround me. The 486 DX2 which I use is however always running even when business interests turn my attention (sometimes for prolonged periods) to matters far removed from chemistry theory. I read somewhere that novices in the field of molecular modelling often make the mistake of calculating too much too quickly. I guess what I'm trying to say is please don't assume that what is slow for some is too slow for everybody. After all, most of us have to go on vacation some time. One of your competitors told me that organic chemists had no need for carring out ab initio calculations on a pc if the calulations would take two weeks to a month (it seemed to escape him that I [>> Continued in next msg]
#23916From: John HubbsSep 28, 1992 8:13 PM
[>> Continued from previous msg] asked the question).