| 查看: 2615 | 回复: 4 | |||||
| 当前主题已经存档。 | |||||
czyzsu专家顾问 (文坛精英)
|
[交流]
【转载】高斯中溶剂模型没有相应溶剂的解决办法[zz]
|
||||
|
CCL:G: Dioxane in PCM with G03 -------------------------------------------------------------------------------- From: Madi Subject: CCL:G: Dioxane in PCM with G03 Date: Thu, 14 May 2009 17:17:29 -0400 -------------------------------------------------------------------------------- Sent to CCL by: Madi [madi7sk{:}gmail.com] Hi, I was trying to specify DMF using the PCM model in gaussian and was having trouble getting it to work. I wrote to gaussian and their reply was very helpful. I am attaching it here for your reference. After you are done specifying all the parameters, make sure you deliberately leave a blank line. Hope this helps madi ****************************************************************** Regarding the specification of the DMF solvent, even though only the solvents listed in the Gaussian documentation of the SCRF keyword are available through the "Solvent=" option, a different solvent may be explicitly defined in the PCM input. However, please, note that it is not possible to input all parameters that define an arbitrary solvent. Specifically, one cannot define all of the terms necessary to calculate the non-electrostatic contributions. For this reason, there are two recommended alternatives: 1) Check in the list of solvents provided for the SCRF keyword whether there is a solvent similar enough to the one you would like to use in your calculation and provide this name in the "Solvent" option of the SCRF keyword. Then, add also the "Read" option to provide the modifications of "EPS", "RSOLV", "Density", and "VMOL". This way, you would be introducing the parameters for your desired solvent, for the cases of parameters that are available for user input, and the program would be using the closest possible parameters for those other parameters that are not available for user input. (Note: if no solvent is specified, then the program defaults to water). So, probably, your best choice here would be to select a solvent with similar properties to DMF (such as Acetonitrile, DMSO, Acetone, Nitromethane) and then give the corresponding "EPS", "RSOLV", "Density", and "VMOL" parameters for DMF. 2) Alternatively, you can use the options "NoDis NoRep NoCav" in order to prevent the use of the parameters for the internally stored solvent (water, by default), so the non-electrostatic contributions to the solvation energy will not be computed. Note also that this is true if the default "UA0" radii are used, the electrostatic terms are then correctly computed this way. If this is not the case, and the UAHF or UAKS radii are used, then there is a dependency of the radius with the type of solvent so it would not be correct to use a modified water as a model for a solvent with properties different than water. Thus, if UAHF or UAKS radii are used, then one has to still do what it was explained in the previous paragraph, 1), (that is to use the "Solvent=" option to select, from the list of predefined solvents, a solvent similar to the one you are interested in). For the calculation of only the electrostatic part, you can define the solvent by defining the dielectric constant using the "EPS" option, the solvent radius using the "RSOLV" option, and the dielectric constant at infinite frequency using the "EPSINF". EPSINF is the dielectric constant at infinite frequency. It is an experimentally measured value. The value of EPSINF will not affect the results of ground-state energy calculation of DeltaG(solv), so you do not need to worry about getting a non-standard value. EPSINF is only necessary when performing excited state PCM calculations. Other than EPS, RSOLV usually makes the biggest difference in the result. As for choosing a value of RSOLV, our developers recommended a couple of references they use for this parameter: * O. Sinanoglu in "Molecular Interactions" Z. Ratajczack and Orville-Thomas eds, Wiley 1982 vol. 3 chap. 6 * M.H. Abraham and A. Nasehzadeh J. Chem. Soc. Faraday Trans. 1 1981, 77, 340. If the desired RSOLV parameter is not found there, the molar volume may give some hints for RSOLV values for new solvents, even though attention must be paid to the chemical nature of the solvent, because strong interactions modify, to some extent, the definition of the hard sphere analogue of the solvent molecular volume. Thus, although a simple approximation, the Stearn-Eyring formula (Stearn, Eyring J. Chem. Phys. 1937, 5, 113) can be used to derive a value for RSOLV, given the molar volume of the solvent. Finally, about other parameters, which only affect the non-electrostatic terms, the DENSITY is a "number density", i.e. number of particles per volume, and the units are Angstroms^–3. Taking the density of the solvent in "g * cm^–3", divide by the molar weight of the solvent and multiply by Avogadro's number to transform the grams into number of particles, and, finally, divide by 10^24 to transform cm^–3 into Angstrom^–3. Although the density and the molar volume are obviously related, these two quantities are used by the program to compute two different energy contributions and, unfortunately, they were introduced in the code as two independent parameters. Thus, the user needs to provide both the density ("DENSITY" PCM keyword) and the molar volume ("VMOL" PCM keyword) in order to set the desired modifications, so the two of them are consistent with each other. The units of the density are Angstrom^–3 while the units of the molar volume are cm^3. So, from the value of the density of the solvent in "g * cm^–3", take the inverse of the density and then divide by the molar weight of the solvent in order to get the molar volume in cm^3. Enter this parameter in your input file with the "VMOL" PCM keyword. The following input will use the standard parameters for the solvent given in "Solvent= additional PCM options (the values given below are not representative of any particular solvent and they are only given as an example of how to use these keywords): % #P SCRF=(PCM,Solvent= EPS=6.0 RSOLV=1.6 DENSITY=0.045 VMOL=19.05 [zz] http://www.ccl.net/cgi-bin/ccl/message-new?2009+05+14+009 |
» 猜你喜欢
研究生招生
已经有0人回复
河海大学 (211.双一流高校)钙钛矿光电实验室丁勇团队欢迎有志青年!
已经有3人回复
物理化学论文润色/翻译怎么收费?
已经有240人回复
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
26博士申请
已经有0人回复
一志愿南师大0703化学 275求调剂
已经有6人回复
招生材料与化工、环境科学调剂研究生
已经有5人回复
A 区过线就可申!重庆双非公办高校 0856调剂,接受跨专业!纳米材料 / 润滑热门方向
已经有0人回复
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
中国电机工程学报的投稿系统
已经有0人回复
wxjbuilder
木虫 (知名作家)
- 应助: 1 (幼儿园)
- 贵宾: 0.05
- 金币: 2986.4
- 散金: 1761
- 红花: 2
- 沙发: 1
- 帖子: 5392
- 在线: 273.2小时
- 虫号: 442356
- 注册: 2007-10-27
- 性别: GG
- 专业: 理论和计算化学

2楼2010-01-21 14:15:44
3楼2010-01-22 09:30:38
quantum999
金虫 (著名写手)
- 应助: 15 (小学生)
- 金币: 2146.3
- 散金: 39
- 帖子: 1426
- 在线: 198.2小时
- 虫号: 41283
- 注册: 2004-03-05
- 专业: 凝聚态物性 II :电子结构
4楼2010-01-27 15:59:47
5楼2010-03-23 21:43:21















回复此楼