{"id":356,"date":"2023-06-15T11:55:55","date_gmt":"2023-06-15T09:55:55","guid":{"rendered":"https:\/\/perso.ens-lyon.fr\/carine.michel\/?page_id=356"},"modified":"2023-06-15T21:31:26","modified_gmt":"2023-06-15T19:31:26","slug":"summer-school-emie-part-1","status":"publish","type":"page","link":"https:\/\/perso.ens-lyon.fr\/carine.michel\/enseignements\/summer-school-emie-part-1\/","title":{"rendered":"Summer School EMIE &#8211; Part 1"},"content":{"rendered":"<h2><span class=\"category\" style=\"color: #ff6600;\"><strong>\u00a0Proton transfer and Microsolvation<\/strong> <\/span><\/h2>\n<p>Date: June 2023<\/p>\n<p>C. Michel<\/p>\n<p>Our aim is to investigate how a proton transfer can be tuned by the presence of water. We will use the software <tt>Gaussian<\/tt> (<a href=\"http:\/\/gaussian.com\/man\/\">online manual<\/a>). The preparation of input files and the visualisation of results will be done using <tt>Avogadro<\/tt>.<\/p>\n<p>We will study the tautomerization reaction of 2-hydroxypyridine \/ 2-pyridone.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-361\" src=\"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-content\/uploads\/2023\/06\/tautomerization.jpeg\" alt=\"\" width=\"278\" height=\"79\" \/><\/p>\n<p>Our calculations will be performed at the PM6 level for the sake of efficiency. Comparison with experiments and other levels of calculations can be made by comparison with the reference article <a href=\"https:\/\/pubs.acs.org\/doi\/10.1021\/ct800477y\"> J. L. Sonnenberg et al. J. Chem. Theory Comput. 2009, 5, 4, 949-961 <\/a>. The sequence is split in three steps and is expected to take 1h30.<\/p>\n<p>&nbsp;<\/p>\n<ul>\n<li><a href=\"#vac\">1. Tautomerization in vacuum<\/a><\/li>\n<li><a href=\"#PCM\">2. Using a continuum model<\/a><\/li>\n<li><a href=\"#micro\">3. Microsolvation approach<\/a><\/li>\n<\/ul>\n<h3><a name=\"vac\"><\/a><br \/>\n<span class=\"subcategory\"> 1. Tautomerization in vacuum<\/span><\/h3>\n<p>We start by studying the reactant and product.<\/p>\n<ol>\n<li>Build the 2-pyridone (PY) and the 2-hydroxypyridine (HY) using Avogadro.<\/li>\n<li>Optimize those structures at the PM6 level using Gaussian and compute frequencies to get Gibbs free energies.<br \/>\nHere is a typical input file structure:<\/p>\n<pre>#n PM6 Opt Freq\r\n[blanck line]\r\n[Title] Optimization of PY\r\n[blanck line]\r\n[Multiplicity,Charge] 0,1\r\n[xyz formated coordinates]\r\n[blank line]\r\n<\/pre>\n<\/li>\n<li>Check the absence of an imaginary frequency using Avogadro or reading the .log file (<tt>NImag=<\/tt>).<\/li>\n<li>Determine the reaction energy. The Gibbs energy is reported in the <tt> - Thermochemistry - <\/tt> paragraph of the .log file as the <tt> Sum of electronic and thermal Free Energies= <\/tt> in Hartree\/particle.<\/li>\n<\/ol>\n<p>We move to the transition state search. We will first perform a scan along a well-chosen coordinate and consider that the structure at the highest energy is a good approximation of the TS structure, close enough to switch to a eigen-follow method to optimize the TS.<\/p>\n<ol>\n<li>Measure the N-H and O-H distance in both molecules and identify the labels of O, N, and the H to be transferred.<\/li>\n<li>Choose a starting structure and the variable to be used in the scan (typically O&#8230;H or N&#8230;H distance).<\/li>\n<li>Prepare the corresponding input file. You can follow this exemple that corresponds to stirring the N1-H8 distance in PY using 2 steps of 0.15 \u212b. Adapt the number of steps and the step size!\n<pre>#n PM6 Opt=ModRedundant\r\n\r\nH transfer py to HY\r\n\r\n0,1\r\nN  -1.135204 0.239458 0.000000\r\nC  0.000000 1.069179 0.000000\r\nC  1.244256 0.320141 0.000000\r\nC  1.264957 -1.051903 0.000000\r\nC  0.056462 -1.809209 0.000000\r\nC  -1.129901 -1.126884 0.000000\r\nO  -0.123527 2.292845 0.000000\r\nH  -2.017014 0.740588 0.000000\r\nH  2.157927 0.905065 0.000000\r\nH  2.220524 -1.570494 0.000000\r\nH  0.059757 -2.891970 0.000000\r\nH -2.101193 -1.610110 0.000000\r\n\r\n1 8 S 2 0.15\r\n\r\n<\/pre>\n<\/li>\n<li>You can plot the energy in function of the coordinate using this <a href=\"https:\/\/perso.ens-lyon.fr\/carine.michel\/enseignements\/script-extract-pes\/\">script<\/a> which takes only the .log file as an argument and produces a image if the profile as a output, a table of the energies and the coordinates of the optimized structure at each step. If needed, you can refine the scan by adding more points close to the maximum: start from a structure close to this maximum that is provided by the script, and relauch a scan with smaller steps. If needed, change reaction coordinate, switching from N..H to O..H.<\/li>\n<li>Select the structure the closest to the maximum and optimize the transition state using the <tt>TS<\/tt> option in <tt>opt<\/tt>. The <tt>calcfc<\/tt> option allows to compute the hessian prior starting the optimization. The <tt>noeigentest<\/tt> allows to perform the optimisation even if more than one imaginary frequency is found in the starting hessian. You can try with and without those keywords.\n<pre>#n PM6 Opt=(TS,calcfc,noeigentest) freq\r\n<\/pre>\n<\/li>\n<li>Check that the structure has been optimized and shows one and only one imaginary frequency. Visualize the corresponding mode of vibration.<br \/>\nDoes it correspond to the reaction coordinate?<\/li>\n<li>Determine the Gibbs free energy of activation of the tautomerization.<\/li>\n<\/ol>\n<h3><a name=\"PCM\"><\/a><br \/>\n<span class=\"subcategory\"> 2. Using a continuum model<\/span><\/h3>\n<p>This tautomerization is known to be strongly facilitated by the presence of water. Repeat the previous calculation adding a continuum model for the water solvent adding the keyword <tt>SCRF=(Solvent=water)<\/tt>. Compare the reaction energies and activation energies with the ones obtained in vacuum and conclude.<\/p>\n<h3><a name=\"micro\"><\/a><br \/>\n<span class=\"subcategory\"> 3. Microsolvation approach<\/span><\/h3>\n<p>We move now to a micro-solvation approach by adding water molecules. We start with one water molecule that can be H-bonded to O and\/or to N:<\/p>\n<ol>\n<li>Determine the best structure with a water molecule H-bonded to PY in vacuum at the PM6 level of theory. Measure the NH&#8230;Ow and the Hw&#8230;O distances.<\/li>\n<li>Determine the best structure with a water molecule H-bonded to HY in vacuum at the PM6 level of theory. Measure the N&#8230;Hw and the OH&#8230;Ow distances.<\/li>\n<li>How is the reaction energy modified by the inclusion of a H-bonded water molecule?<\/li>\n<li>Determine the best transition state structure for the tautomerization.<\/li>\n<li>How is the reaction barrier modified by the inclusion of a H-bonded water molecule?<\/li>\n<\/ol>\n<p>To ensure that most of the effect is indeed obtained using one water molecule, we need to investigate what happens with more water molecules and upon inclusion of a continuum model.<\/p>\n<ol>\n<li>Repeat with two water molecules.<\/li>\n<li>Repeat with three water molecules.<\/li>\n<li>Repeat including the PCM.<\/li>\n<li>Conclude.<\/li>\n<\/ol>\n<h3>Bonuses<\/h3>\n<ol>\n<li>Let&rsquo;s consider now methanol as a solvent. Is the H-transfer also assisted by this solvent?<\/li>\n<li>Can ammonia also catalyze this tautomerization ? See for instance <span class=\"cit-title\"><i>J. Phys. Chem. A<\/i><\/span> <strong><span class=\"cit-year-info\">2013<\/span><\/strong><span class=\"cit-volume\"><strong>,<\/strong> <em>117<\/em><\/span><span class=\"cit-issue\">, <\/span><span class=\"cit-pageRange\">7523\u20137534<\/span>.<\/li>\n<li>Can you propose an aromatic susbtitution that would favor the HY form?<\/li>\n<li>To move to a fully-explicit description of the water solvent, we could solvate the molecule in a box of liquid water and use ab initio metadynamics. Which collective variables would you suggest to use? Why?<\/li>\n<li>The PY\/HYD can easily dimerize. How could you assess the impact of the solvent on this dimerization?<\/li>\n<li>In this dimer, the proton transfer is not observed in the ground state but can be observed in the exctited state. If interested, see <em>Chemical Physics Letters<\/em><br \/>\n<strong>2019<\/strong>, <em>730<\/em>, 634-637.<\/li>\n<\/ol>\n<h3><a name=\"Optimized structures\"><\/a><br \/>\n<span class=\"subcategory\"> Some optimized structures in .xyz format<\/span><\/h3>\n<pre>12\r\nTS in vacuum\r\n N   -1.14643    0.15373   -0.00000\r\n C   -0.00000    0.98749   -0.00000\r\n C    1.29669    0.45162    0.00000\r\n C    1.39237   -0.93644    0.00000\r\n C    0.24893   -1.76106    0.00000\r\n C   -1.03284   -1.19466    0.00000\r\n O   -0.57315    2.14740   -0.00000\r\n H   -1.78463    1.40726   -0.00000\r\n H    2.16418    1.10316    0.00000\r\n H    2.38019   -1.40446    0.00000\r\n H    0.36258   -2.84078    0.00000\r\n H   -1.94292   -1.80223   -0.00000\r\n<\/pre>\n<pre>15\r\nTS with one water molecule\r\n C   -0.23694    0.68650   -0.01880\r\n C    1.03891    1.32241    0.01873\r\n C    2.16848    0.53222    0.03157\r\n C    2.05929   -0.87898    0.01037\r\n C    0.79986   -1.46218   -0.02365\r\n N   -0.34095   -0.70633   -0.03219\r\n H    3.16022    0.98842    0.05751\r\n H    1.07590    2.40881    0.03077\r\n H    2.95459   -1.49235    0.01835\r\n H    0.66295   -2.54918   -0.04390\r\n O   -1.35507    1.33415   -0.03129\r\n H   -2.38055    0.54579   -0.12789\r\n H   -3.05045   -0.70584    0.88617\r\n O   -2.80644   -0.62790   -0.04509\r\n H   -1.72159   -1.10113   -0.09384\r\n<\/pre>\n<pre>18\r\nTS with two water molecules\r\n N   -0.09023    0.60184   -0.00385\r\n C   -0.30033   -0.77024   -0.02935\r\n C   -1.62996   -1.30705   -0.02201\r\n C   -2.70025   -0.44488    0.00960\r\n C   -2.48379    0.95567    0.03613\r\n C   -1.18535    1.43502    0.03158\r\n O    0.71270   -1.56826   -0.05310\r\n H    1.87034   -1.09008   -0.13952\r\n H   -1.74049   -2.38907   -0.04273\r\n H   -3.72283   -0.82688    0.01480\r\n H   -3.32960    1.63549    0.06187\r\n H   -0.96477    2.50888    0.05793\r\n O    3.05526   -0.73380   -0.01221\r\n H    3.26549   -1.02110    0.88362\r\n H    2.92822    0.42801    0.06722\r\n O    2.29536    1.58299    0.03719\r\n H    1.20178    1.18505    0.05272\r\n H    2.41486    1.89834   -0.85958\r\n<\/pre>\n<pre>21\r\nTS with three water molecules\r\n N   -0.49753   -0.53204    0.22683\r\n C   -0.80998    0.80611    0.06022\r\n C   -2.15069    1.22941   -0.19997\r\n C   -3.14207    0.27662   -0.28479\r\n C   -2.82699   -1.09210   -0.11502\r\n C   -1.51154   -1.45406    0.13540\r\n O    0.13778    1.69352    0.14400\r\n H    1.24819    1.27276    0.36647\r\n H   -2.33953    2.29295   -0.32309\r\n H   -4.17521    0.56753   -0.48290\r\n H   -3.60760   -1.84319   -0.18188\r\n H   -1.21527   -2.50121    0.27370\r\n O    1.91836   -1.43378    0.66260\r\n H    0.87023   -1.03157    0.49043\r\n H    2.29488   -0.76661    1.28205\r\n O    3.05685   -0.42795   -1.26000\r\n H    2.46637   -1.13444   -0.45393\r\n H    3.98072   -0.65299   -1.13009\r\n O    2.49154    0.98336    0.55798\r\n H    2.89246    1.82425    0.73575\r\n H    2.87884    0.57973   -0.57588\r\n<\/pre>\n<hr \/>\n<p>Last update : June 2023<\/p>\n<p>&nbsp;<\/p>\n<p>&nbsp;<\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u00a0Proton transfer and Microsolvation Date: June 2023 C. Michel Our aim is to investigate how a proton transfer can be tuned by the presence of water. We will use the software Gaussian (online manual). The preparation of input files and the visualisation of results will be done using Avogadro. We will study the tautomerization reaction [&hellip;]<\/p>\n","protected":false},"author":2,"featured_media":0,"parent":2,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-356","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/pages\/356","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/comments?post=356"}],"version-history":[{"count":15,"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/pages\/356\/revisions"}],"predecessor-version":[{"id":400,"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/pages\/356\/revisions\/400"}],"up":[{"embeddable":true,"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/pages\/2"}],"wp:attachment":[{"href":"https:\/\/perso.ens-lyon.fr\/carine.michel\/wp-json\/wp\/v2\/media?parent=356"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}