From owner-chemistry@ccl.net Fri Apr 7 02:43:00 2017 From: "MIELCZAREK Detlev Conrad detlev-conrad.mielczarek()ifpen.fr" To: CCL Subject: CCL:G: Thermodynamic Data & Solvation - Calculation Questions: Message-Id: <-52729-170407024208-30578-7NHomkiqazGj/mW55YMGAQ%x%server.ccl.net> X-Original-From: MIELCZAREK Detlev Conrad Content-Language: en-US Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset="utf-8" Date: Fri, 7 Apr 2017 06:41:58 +0000 MIME-Version: 1.0 Sent to CCL by: MIELCZAREK Detlev Conrad [detlev-conrad.mielczarek * ifpen.fr] Dear Uwe, many thanks for the explanation and the reference/paper. Dear Eric, many thanks for your very comprehensive response - with a detailed explanation how and why the terms are chosen - possibly the best I have seen yet. Hopefully I won't have any further questions after sitting down to properly read your response and the linked paper. Once again, many thanks for your two comprehensive responses. Kind regards Detlev Detlev Conrad Mielczarek Scientific Visitor/Post Doctorant IFP Energies nouvelles France detlev-conrad.mielczarek - - ifpen.fr -----Original Message----- > From: owner-chemistry+detlev-conrad.mielczarek==ifpen.fr|ccl.net [mailto:owner-chemistry+detlev-conrad.mielczarek==ifpen.fr|ccl.net] On Behalf Of Eric Hermes erichermes**gmail.com Sent: 06 April 2017 18:57 To: MIELCZAREK Detlev Conrad Subject: CCL:G: Thermodynamic Data & Solvation - Calculation Questions: Sent to CCL by: Eric Hermes [erichermes!^!gmail.com] Dr. Mielczarek, As always, the free energy of a species is defined with respect to some choice of reference state, such that for a reaction of the form: (1) A + B -> C with a standard state free energy of reaction ΔG°, one can write the equilibrium constant as: (2) Keq = exp[-ΔG°/(kB T)] = a_C / (a_A a_C) Where a is the *activity* of a species. In the case of an ideal gas, the activity is simply the density divided by the *standard state* density: (3) a = ρ/ρ° The numerical value of Keq then depends on choice of reference state, but the equilibrium densities do not. Eq. 2 satisfies this requirement. One can show that the free energy of a species at an *arbitrary* reference state is given by: (4) G = G° + kB T ln[ρ/ρ°] This is *general*, but you can show that for an ideal gas the expression derives from the translational entropy. Now, when we are discussing free energy of solvation, the reaction becomes: (5) A_g -> A_soln Which of course brings with it changes to both the enthalpy and the entropy. These issues are *orthogonal* to the reference state issue, though. The standard state for species in the gas phase is usually taken to be 1 bar (occasionally 1 atm is chosen instead), but the standard state for species in *solution* is typically chosen to be 1 M, which is very different! With these choices of reference state, the free energy of solvation ΔG°solv can be used to calculate the equilibrium constant like this (using P for pressure to visually distinguish from the density ρ): (6) Keq = exp[-ΔG°solv/(kB T)] = (ρ_A/(1 M)) / (P_A / (1 bar)) As before, the numerical value of Keq depends on the choice of reference state, but the equilibrium density and pressure do not. So, if your goal is to determine the equilibrium properties of a solvation process, then you do not need to be concerned about the choice of reference state -- you need to be *aware* of it in order to calculate the properties directly, but whichever choice you make you will get the same answer. On the other hand, if you wish to compare your free energies of solvation to experimental values, you should most definitely ensure that the value you are calculating uses a reference state of 1 M for solution-phase species and 1 bar for gas-phase species. In the case of SMD, by default it will use the *same* reference state for the gas-phase and solution-phase species. The actual value of the reference state is irrelevant, as ΔGsolv has the same value for any choice of reference state so long as the reactants and products have the *same* reference state (if you are not convinced, stare at eqs 2, 3, and 4 in the context of eq 5 until you believe me :) ). This means you can arbitrarily say what reference state ΔGsolv is at, say 1 bar, for both reactants and products. Then, if you want to calculate ΔG°solv (the *experimental* standard state value), you just need to use eq 4: (7) ΔGsolv = Gsoln - G°gas (we choice a reference state of 1 bar) (8) Gsoln = G°soln + kB T ln[(1 bar / (kB T))/1 M] (ideal gas law) (9) ΔG°solv = G°soln - G°gas = ΔGsolv - kB T ln[(1 bar / (kB T))/1 M] Note that you would have gotten the *exact same result* if you had arbitrarily chosen a reference state of 1 M instead for ΔGsolv (if you don't believe me, try it yourself -- start with ΔGsolv = G°soln - Ggas). --- Now, as to your point about the entropy of solution phase species -- I'm not sure the paper you are linking is making the claim you are saying it does. It is an *unarguable fact* that (at the typical choice of standard states) a species in solution has less entropy than the same species in the gas phase. This *must* be the case because the molecules in an ideal gas are non-interacting, uncorrelated, and undergoing ballistic motion, whereas in solution molecules are *caged* by the solvent and undergoing diffusive motion, which means they have significantly less freedom of motion. The paper you link talks about *vibrational* motion, which is going to be *significantly* less perturbed by the presence of solvent. The *vast* majority of entropy for gas-phase species comes from translational motion, not vibrational motion. --- Finally, I'm not sure SMD is well-suited to give a *breakdown* of enthalpic and entropic terms to the free energy of solvation. The procedure you discuss will get *some* of that breakdown, but ultimately several contributions to ΔGsolv are all entangled with one another.  What I mean by that is if you simply perform two *single point* calculations on the same species, one with and one without SMD correction, the difference in the potential energies between those two calculations will include a mixture of enthalpic and entropic effects.  Doing the full thermodynamic calculations (i.e. doing a frequency calculation and reading the thermodynamic data printed in the Gaussian output file) will give you a *full* estimate for ΔGsolv. But, in both calculations (with and without SMD) the calculated enthalpic and entropic contributions are arising from the same set of approximations -- harmonic oscillator for vibration, rigid rotor for rotation, and *ideal gas* for translation. This works because the *potential energy* difference between the systems is parameterized such that the free energy difference calculated by Gaussian is a good approximatino of the total free energy of solvation. Since the enthalpy and entropy calculated by Gaussian for both calculations are using the same approximations, they will have very similar (but not identical!) values. All of those complicated factors such as the loss of translational entropy due to solvation are baked into the *electronic energy* that is calculated by the SMD method. The most *robust* way of separating enthalpic and entropic contributions would be to calculate the full value of the free energy at a variety of different temperatures, then numerically evaluating dG/dT to find the entropy. However, SMD as implemented in Gaussian is not temperature-dependent, so you cannot actually do this. There are other continuum solvation models developed by Cramer and Truhlar which do have temperature-dependence, which would allow you to do this, but as far as I can tell they are not included in Gaussian. --- I hope this answers all of the questions that you had. Please let me know if anything I said was unclear or you want additional assistance. Eric Hermes On Thu, 2017-04-06 at 06:24 +0000, MIELCZAREK Detlev Conrad detlev- conrad.mielczarek-.-ifpen.fr wrote: > Sent to CCL by: MIELCZAREK Detlev Conrad [detlev- > conrad.mielczarek]^[ifpen.fr] Dear CCL, a question on thermodynamic > data & solvation from me, maybe you can help me. > > So, the basic problem for me is, that I am calculating/want to > calculate thermodynamic data (Hf, S - hence also dG) in solvation, > using implicit solvation models, SMD with a COSMO cavity to be > specific. For my application, these should be accurate enough. (So no > molecular dynamics simulations etc.) > > Solvation models are normally parametrised for dGsolv - so this value > can be extracted from the quantum chemistry calculation as the > difference in the calculated Gibbs Free Enthalpy. > Hf can calculated easily in the gas phase, and a re-optimisation of > the structure with solvation should capture the majority of the impact > of solvation on the enthalpy. (Which is dominated by molecular > structure.) > (I guess there is the case of stabilisation and complexes, such as are > reported for water. However this is currently beyond the scope of my > work.) > > The topic of solvation has been discussed previously on the CCL here: > http://www.ccl.net/chemistry/resources/messages/2011/12/01.001-dir/ > http://www.ccl.net/chemistry/resources/messages/2011/10/06.005-dir/ > http://www.ccl.net/chemistry/resources/messages/2014/05/01.004-dir/ > And there is the book "Essentials of Computational Chemistry Theories > and Models" from Professor Cramer with a section on phase change (the > source of my confusion). > > Specifically, the discussion concerning the energy change related to > the state conversion causes me some grief. > > On the one hand, the CCL responses read as if this should be applied > in the case of any phase change, but then others suggest this is > applicable only if the process is a second order reaction and thus the > total number of mols changes? - The latter view seems to agree with > the book... > > So if I have compound A in both the gas and liquid phase (from a > quantum chemistry calculation), do I need to account for the phase > change/change of state or not? Or is it something that can be included > in the parametrisation of the solvation model/the quantum chemistry > code already? > > Just to add more confusion to the topic: I have trialled a commercial > product which gives the Gibbs Enthalpy of Solvation in kcal/mol for > mol/L concentrations and using a very low end/fast functional, it > gives values similar to when a correction term is added... on the > other hand, where available, the calculated values without correction > agree with the published values in the SMD paper: http://pubs.acs.org > /doi/abs/10.1021/jp810292n (Supplementary Data) > > In addition, a regular computational chemistry calculation sees very > little (virtually no) difference in the entropy between the gaseous > and solvated phase. This would agree with the CCL-linked paper here h > ttp://pubs.acs.org/doi/abs/10.1021/jp205508z . But this would clash > with the common expectation that entropy in the liquid phase is > reduced... > > Hence, I would highly appreciate if someone knowledgeable in the field > of solvation could guide me onto the correct track. > > Detlev Conrad Mielczarek > Scientific Visitor/Post Doctorant > IFP Energies nouvelles > France > > www.ifpenergiesnouvelles.fr > > __________________________ > Avant d'imprimer, pensez à l'environnement ! Please consider the > environment before printing ! > Ce message et toutes ses pièces jointes sont confidentiels et établis > à l'intention exclusive de ses destinataires. Toute utilisation non > conforme à sa destination, toute diffusion ou toute publication, > totale ou partielle, est interdite, sauf autorisation expresse. IFP > Energies nouvelles décline toute responsabilité au titre de ce > message. This message and any attachments are confidential and > intended solely for the addressees. Any unauthorised use or > dissemination is prohibited. IFP Energies nouvelles should not be > liable for this message. > __________________________ > > > > -= This is automatically added to each message by the mailing script > =- > To recover the email address of the author of the message, please > change>       http://www.ccl.net/cgi-bin/ccl/send_ccl_message>        > change>       http://www.ccl.net/chemistry/sub_unsub.shtml> Conferences: http://server.ccl.net/chemistry/announcements/conference > s/>       http://www.ccl.net/spammers.txthttp://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt__________________________ Avant d'imprimer, pensez à l'environnement ! Please consider the environment before printing ! Ce message et toutes ses pièces jointes sont confidentiels et établis à l'intention exclusive de ses destinataires. Toute utilisation non conforme à sa destination, toute diffusion ou toute publication, totale ou partielle, est interdite, sauf autorisation expresse. IFP Energies nouvelles décline toute responsabilité au titre de ce message. This message and any attachments are confidential and intended solely for the addressees. Any unauthorised use or dissemination is prohibited. IFP Energies nouvelles should not be liable for this message. __________________________ From owner-chemistry@ccl.net Fri Apr 7 03:18:00 2017 From: "adon cumi adonmage|-|gmail.com" To: CCL Subject: CCL:G: Difficulty extracting pdb information form output file Message-Id: <-52730-170406022724-28590-Jup9EXX9O0E+06YVKZi5OQ###server.ccl.net> X-Original-From: adon cumi Content-Type: multipart/alternative; boundary=001a113b02cab7faf0054c799961 Date: Thu, 6 Apr 2017 08:27:16 +0200 MIME-Version: 1.0 Sent to CCL by: adon cumi [adonmage]-[gmail.com] --001a113b02cab7faf0054c799961 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: quoted-printable Hi to all!!! You are triyng to get geometrical data in pdb format. Download the software MERCURY 3.8 from the CCDC database. Save your optimize gaussian g=C3=A9om= =C3=A9try in .mol or .sdf load this .mol file with MERCURY to visualize the structure. Now save as .pdb file! It work perfectly . It can r=C3=A9solve your extensi= on file probl=C3=A8m. Use your .out or .log gaussian file Best regards! Le 5 avr. 2017 22:26, "Close, David M. CLOSED..mail.etsu.edu" < owner-chemistry[A]ccl.net> a =C3=A9crit : > > Sent to CCL by: "Close, David M." [CLOSED:+:mail.etsu.edu] > Vidya: > Not sure what you mean by pdb information? The Protein Data Bank lists > the structure of proteins determined by x-ray crystallography. They list > the xyz coordinates of these structures. Surely the Gaussian out also > lists the xyz coordinates in the output file. Are you trying to convert > one set of data to another? If so there is the old utility called BABEL > that will do just that. > Regards, Dave Close > > -----Original Message----- > > From: owner-chemistry+closed=3D=3Detsu.edu]-[ccl.net [mailto: > owner-chemistry+closed=3D=3Detsu.edu]-[ccl.net] On Behalf Of Vidya > Kaipanchery vidyak.clt-.-gmail.com > Sent: Wednesday, April 05, 2017 8:56 AM > To: Close, David M. > Subject: CCL: Difficulty extracting pdb information form output file > > > Sent to CCL by: "Vidya Kaipanchery" [vidyak.clt**gmail.com] I have done > a ONIOM calculation in G09 and got my optimized structure. It does not > contain any of the pdb (protein data base) information which I require fo= r > further calculation(frequency analysis). I tried extracting it form chk > file using "newzmat" command and it failed. I also tried to convert the c= hk > to fchk and it didnt work. Can you please help me how to extract my > optimized geometry with pdb information from my chk file or even log file= ? > > Kind regards, > Vidyahttp://www.ccl.net/cgi-bin/ccl/send_ccl_messagehttp:/ > /www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.net/spammers.txt > > > -=3D This is automatically added to each message by the mailing script = =3D-> > > --001a113b02cab7faf0054c799961 Content-Type: text/html; charset=UTF-8 Content-Transfer-Encoding: quoted-printable

Hi to all!!!
You are triyng to get geometrical data in pdb format. Download the software= MERCURY 3.8 from the CCDC database. Save your optimize gaussian g=C3=A9om= =C3=A9try in .mol or .sdf load this .mol file with MERCURY to visualize the= structure. Now save as .pdb file! It work perfectly . It can r=C3=A9solve = your extension file probl=C3=A8m. Use your .out or .log gaussian file
Best regards!

Le 5 avr. 2017 22:26, "Close, David M. CLOS= ED..mail.etsu.edu" <owner-chemistry[A]ccl.net> a =C3=A9cr= it=C2=A0:

Sent to CCL by: "Close, David M." [CLOSED:+:mail.etsu.edu]
Vidya:
=C2=A0 Not sure what you mean by pdb information?=C2=A0 The Protein Data Ba= nk lists the structure of proteins determined by x-ray crystallography.=C2= =A0 They list the xyz coordinates of these structures.=C2=A0 Surely the Gau= ssian out also lists the xyz coordinates in the output file.=C2=A0 Are you = trying to convert one set of data to another?=C2=A0 If so there is the old = utility called BABEL that will do just that.
=C2=A0 Regards, Dave Close

-----Original Message-----
> From: owner-chemistry+closed=3D=3Detsu.edu]-[ccl.net [mailto:owner-chemistry+closed=3D=3Detsu.edu]-[ccl.net] On B= ehalf Of Vidya Kaipanchery vidyak.clt-.-gmail.com
Sent: Wednesday, April 05, 2017 8:56 AM
To: Close, David M. <CLOSED]-[mail.etsu.edu>
Subject: CCL: Difficulty extracting pdb information form output file


Sent to CCL by: "Vidya=C2=A0 Kaipanchery" [vidyak.clt**gmail.com] I ha= ve done a ONIOM calculation in G09 and got my optimized structure. It does = not contain any of the pdb (protein data base) information which I require = for further calculation(frequency analysis). I tried extracting it form chk= file using "newzmat" command and it failed. I also tried to conv= ert the chk to fchk and it didnt work. Can you please help me how to extrac= t my optimized geometry with pdb information from my chk file or even log f= ile?

Kind regards,
Vidyahttp://www.ccl.net/cgi-bin/ccl/send_ccl_mes= sagehttp://www.ccl.net/chemistry/sub_unsub.shtmlhttp://www.ccl.ne= t/spammers.txt


-=3D This is automatically added to each message by the mailing script =3D-=
E-mail to subscribers: CHEMISTRY[A]ccl.n= et or use:
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/cgi-bin/= ccl/send_ccl_message

E-mail to administrators: CHEM= ISTRY-REQUEST[A]ccl.net or use
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/cgi-bin/= ccl/send_ccl_message

Subscribe/Unsubscribe:
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/chemistry/s= ub_unsub.shtml

Before posting, check wait time at: http://www.ccl.net

Job: http://www.ccl.net/jobs
Conferences: http://server.ccl.net/ch= emistry/announcements/conferences/

Search Messages: http://www.ccl.net/chemistry/searchccl/index.shtml
=C2=A0 =C2=A0 =C2=A0 http://www.ccl.net/spammers.txt

RTFI: http://www.ccl.net/chemistry/aboutcc= l/instructions/


--001a113b02cab7faf0054c799961-- From owner-chemistry@ccl.net Fri Apr 7 03:53:00 2017 From: "adon cumi adonmage- -gmail.com" To: CCL Subject: CCL:G: Thermodynamic Data & Solvation - Calculation Questions: Message-Id: <-52731-170407031229-7542-JJdb2hnNVHiEfexvG++vlQ ~ server.ccl.net> X-Original-From: adon cumi Content-Transfer-Encoding: 8bit Content-Type: text/plain; charset=UTF-8 Date: Fri, 7 Apr 2017 09:12:22 +0200 MIME-Version: 1.0 Sent to CCL by: adon cumi [adonmage .. gmail.com] HI TO ALL! The most *robust* way of separating enthalpic and entropic contributions would be to calculate the full value of the free energy at a variety of different temperatures, then numerically evaluating dG/dT to find the entropy The JAGUAR program from SCHRODINGER Inc. It provide a tool via jaguar-->single point energy-->vibrational frequencies. in the tab below you should be able to increment the temperature to do your "dG/dT" ! BEST REGARDS, 2017-04-06 17:08 UTC+02:00, Eric Hermes erichermes**gmail.com : > > Sent to CCL by: Eric Hermes [erichermes!^!gmail.com] > Dr. Mielczarek, > > As always, the free energy of a species is defined with respect to some > choice of reference state, such that for a reaction of the form: > > (1) A + B -> C > > with a standard state free energy of reaction ΔG°, one can write the > equilibrium constant as: > > (2) Keq = exp[-ΔG°/(kB T)] = a_C / (a_A a_C) > > Where a is the *activity* of a species. In the case of an ideal gas, > the activity is simply the density divided by the *standard state* > density: > > (3) a = ρ/ρ° > > The numerical value of Keq then depends on choice of reference state, > but the equilibrium densities do not. Eq. 2 satisfies this requirement. > One can show that the free energy of a species at an *arbitrary* > reference state is given by: > > (4) G = G° + kB T ln[ρ/ρ°] > > This is *general*, but you can show that for an ideal gas the > expression derives from the translational entropy. > > Now, when we are discussing free energy of solvation, the reaction > becomes: > > (5) A_g -> A_soln > > Which of course brings with it changes to both the enthalpy and the > entropy. These issues are *orthogonal* to the reference state issue, > though. The standard state for species in the gas phase is usually > taken to be 1 bar (occasionally 1 atm is chosen instead), but the > standard state for species in *solution* is typically chosen to be 1 M, > which is very different! With these choices of reference state, the > free energy of solvation ΔG°solv can be used to calculate the > equilibrium constant like this (using P for pressure to visually > distinguish from the density ρ): > > (6) Keq = exp[-ΔG°solv/(kB T)] = (ρ_A/(1 M)) / (P_A / (1 bar)) > > As before, the numerical value of Keq depends on the choice of > reference state, but the equilibrium density and pressure do not. > > So, if your goal is to determine the equilibrium properties of a > solvation process, then you do not need to be concerned about the > choice of reference state -- you need to be *aware* of it in order to > calculate the properties directly, but whichever choice you make you > will get the same answer. On the other hand, if you wish to compare > your free energies of solvation to experimental values, you should most > definitely ensure that the value you are calculating uses a reference > state of 1 M for solution-phase species and 1 bar for gas-phase > species. > > In the case of SMD, by default it will use the *same* reference state > for the gas-phase and solution-phase species. The actual value of the > reference state is irrelevant, as ΔGsolv has the same value for any > choice of reference state so long as the reactants and products have > the *same* reference state (if you are not convinced, stare at eqs 2, > 3, and 4 in the context of eq 5 until you believe me :) ). > > This means you can arbitrarily say what reference state ΔGsolv is at, > say 1 bar, for both reactants and products. Then, if you want to > calculate ΔG°solv (the *experimental* standard state value), you just > need to use eq 4: > > (7) ΔGsolv = Gsoln - G°gas (we choice a reference state of 1 bar) > (8) Gsoln = G°soln + kB T ln[(1 bar / (kB T))/1 M] (ideal gas law) > (9) ΔG°solv = G°soln - G°gas = ΔGsolv - kB T ln[(1 bar / (kB T))/1 M] > > Note that you would have gotten the *exact same result* if you had > arbitrarily chosen a reference state of 1 M instead for ΔGsolv (if you > don't believe me, try it yourself -- start with ΔGsolv = G°soln - > Ggas). > > --- > > Now, as to your point about the entropy of solution phase species -- > I'm not sure the paper you are linking is making the claim you are > saying it does. It is an *unarguable fact* that (at the typical choice > of standard states) a species in solution has less entropy than the > same species in the gas phase. This *must* be the case because the > molecules in an ideal gas are non-interacting, uncorrelated, and > undergoing ballistic motion, whereas in solution molecules are *caged* > by the solvent and undergoing diffusive motion, which means they have > significantly less freedom of motion. > > The paper you link talks about *vibrational* motion, which is going to > be *significantly* less perturbed by the presence of solvent. The > *vast* majority of entropy for gas-phase species comes from > translational motion, not vibrational motion. > > --- > > Finally, I'm not sure SMD is well-suited to give a *breakdown* of > enthalpic and entropic terms to the free energy of solvation. The > procedure you discuss will get *some* of that breakdown, but ultimately > several contributions to ΔGsolv are all entangled with one another. > > What I mean by that is if you simply perform two *single point* > calculations on the same species, one with and one without SMD > correction, the difference in the potential energies between those two > calculations will include a mixture of enthalpic and entropic effects. > > Doing the full thermodynamic calculations (i.e. doing a frequency > calculation and reading the thermodynamic data printed in the Gaussian > output file) will give you a *full* estimate for ΔGsolv. But, in both > calculations (with and without SMD) the calculated enthalpic and > entropic contributions are arising from the same set of approximations > -- harmonic oscillator for vibration, rigid rotor for rotation, and > *ideal gas* for translation. > > This works because the *potential energy* difference between the > systems is parameterized such that the free energy difference > calculated by Gaussian is a good approximatino of the total free energy > of solvation. Since the enthalpy and entropy calculated by Gaussian for > both calculations are using the same approximations, they will have > very similar (but not identical!) values. All of those complicated > factors such as the loss of translational entropy due to solvation are > baked into the *electronic energy* that is calculated by the SMD > method. > > The most *robust* way of separating enthalpic and entropic > contributions would be to calculate the full value of the free energy > at a variety of different temperatures, then numerically evaluating > dG/dT to find the entropy. However, SMD as implemented in Gaussian is > not temperature-dependent, so you cannot actually do this. There are > other continuum solvation models developed by Cramer and Truhlar which > do have temperature-dependence, which would allow you to do this, but > as far as I can tell they are not included in Gaussian. > > --- > > I hope this answers all of the questions that you had. Please let me > know if anything I said was unclear or you want additional assistance. > > Eric Hermes > > On Thu, 2017-04-06 at 06:24 +0000, MIELCZAREK Detlev Conrad detlev- > conrad.mielczarek-.-ifpen.fr wrote: >> Sent to CCL by: MIELCZAREK Detlev Conrad [detlev- >> conrad.mielczarek]^[ifpen.fr] >> Dear CCL, a question on thermodynamic data & solvation from me, maybe >> you can help me. >> >> So, the basic problem for me is, that I am calculating/want to >> calculate thermodynamic data (Hf, S - hence also dG) in solvation, >> using implicit solvation models, SMD with a COSMO cavity to be >> specific. For my application, these should be accurate enough. (So no >> molecular dynamics simulations etc.) >> >> Solvation models are normally parametrised for dGsolv - so this value >> can be extracted from the quantum chemistry calculation as the >> difference in the calculated Gibbs Free Enthalpy. >> Hf can calculated easily in the gas phase, and a re-optimisation of >> the structure with solvation should capture the majority of the >> impact of solvation on the enthalpy. (Which is dominated by molecular >> structure.) >> (I guess there is the case of stabilisation and complexes, such as >> are reported for water. However this is currently beyond the scope of >> my work.) >> >> The topic of solvation has been discussed previously on the CCL here: >> http://www.ccl.net/chemistry/resources/messages/2011/12/01.001-dir/ >> http://www.ccl.net/chemistry/resources/messages/2011/10/06.005-dir/ >> http://www.ccl.net/chemistry/resources/messages/2014/05/01.004-dir/ >> And there is the book "Essentials of Computational Chemistry Theories >> and Models" from Professor Cramer with a section on phase change (the >> source of my confusion). >> >> Specifically, the discussion concerning the energy change related to >> the state conversion causes me some grief. >> >> On the one hand, the CCL responses read as if this should be applied >> in the case of any phase change, but then others suggest this is >> applicable only if the process is a second order reaction and thus >> the total number of mols changes? - The latter view seems to agree >> with the book... >> >> So if I have compound A in both the gas and liquid phase (from a >> quantum chemistry calculation), do I need to account for the phase >> change/change of state or not? Or is it something that can be >> included in the parametrisation of the solvation model/the quantum >> chemistry code already? >> >> Just to add more confusion to the topic: I have trialled a commercial >> product which gives the Gibbs Enthalpy of Solvation in kcal/mol for >> mol/L concentrations and using a very low end/fast functional, it >> gives values similar to when a correction term is added... on the >> other hand, where available, the calculated values without correction >> agree with the published values in the SMD paper: http://pubs.acs.org >> /doi/abs/10.1021/jp810292n (Supplementary Data) >> >> In addition, a regular computational chemistry calculation sees very >> little (virtually no) difference in the entropy between the gaseous >> and solvated phase. This would agree with the CCL-linked paper here h >> ttp://pubs.acs.org/doi/abs/10.1021/jp205508z . But this would clash >> with the common expectation that entropy in the liquid phase is >> reduced... >> >> Hence, I would highly appreciate if someone knowledgeable in the >> field of solvation could guide me onto the correct track. >> >> Detlev Conrad Mielczarek >> Scientific Visitor/Post Doctorant >> IFP Energies nouvelles >> France >> >> www.ifpenergiesnouvelles.fr >> >> __________________________ >> Avant d'imprimer, pensez à l'environnement ! Please consider the >> environment before printing ! >> Ce message et toutes ses pièces jointes sont confidentiels et établis >> à l'intention exclusive de ses destinataires. Toute utilisation non >> conforme à sa destination, toute diffusion ou toute publication, >> totale ou partielle, est interdite, sauf autorisation expresse. IFP >> Energies nouvelles décline toute responsabilité au titre de ce >> message. This message and any attachments are confidential and >> intended solely for the addressees. Any unauthorised use or >> dissemination is prohibited. IFP Energies nouvelles should not be >> liable for this message. >> __________________________ >> >> >> >> -= This is automatically added to each message by the mailing script >> =- >> To recover the email address of the author of the message, please >> change>       http://www.ccl.net/cgi-bin/ccl/send_ccl_message> >>       http://www.ccl.net/cgi-bin/ccl/send_ccl_message>>       http://www.ccl.net/chemistry/sub_unsub.shtml>> Conferences: http://server.ccl.net/chemistry/announcements/conference >> s/>       http://www.ccl.net/spammers.txt> > > From owner-chemistry@ccl.net Fri Apr 7 16:45:01 2017 From: "AXEL LONTSI adonmage .. gmail.com" To: CCL Subject: CCL:G: service for GAUSSIAN CALCULATION Message-Id: <-52732-170407164350-21303-z8IrXOTAUBlwVEMiamJnNQ(!)server.ccl.net> X-Original-From: "AXEL LONTSI" Date: Fri, 7 Apr 2017 16:43:49 -0400 Sent to CCL by: "AXEL LONTSI" [adonmage-x-gmail.com] HI TO ALL!!! it is been a while that I have start a opt freq simulation of an organic molecule in my laptop using GAUSSIAN 09. but it take too much time!!! I WANT TO KNOW IF THERE ARE SITES WEB OR SERVICES WHICH PROVIDE SERVERS FOR COMPUTATIONAL CHEMISTRY WITH GAUSSIAN (OR JAGUAR) SO THAT I CAN START COMPUTATION ON THEIR SERVERS AND PAY FOR THAT. THIS SOLUTION CAN HELP ME SPEED UP MY RESEARCH!!!! THANK YOU FOR ANY ANSWER!!! AXEL LONTSI TIWA MINISTRY OF RESEARCH AND INNOVATION , CAMEROON