Difference between revisions of "AutoMeKin"

From AutoMeKin
Jump to navigation Jump to search
(News)
(News)
 
(9 intermediate revisions by the same user not shown)
Line 3: Line 3:
 
[[File:Automekin.png|center]]
 
[[File:Automekin.png|center]]
 
AutoMeKin (formerly ''tsscds'') has been designed to discover reaction mechanisms in an automated fashion. Transition states are located using MD simulations and Graph Theory algorithms. Monte Carlo simulations afford kinetic results. The only input is a starting structure in XYZ format. The method is described in these two publications: [http://onlinelibrary.wiley.com/doi/10.1002/jcc.23790/abstract 1]
 
AutoMeKin (formerly ''tsscds'') has been designed to discover reaction mechanisms in an automated fashion. Transition states are located using MD simulations and Graph Theory algorithms. Monte Carlo simulations afford kinetic results. The only input is a starting structure in XYZ format. The method is described in these two publications: [http://onlinelibrary.wiley.com/doi/10.1002/jcc.23790/abstract 1]
[http://pubs.rsc.org/en/content/articlelanding/2015/cp/c5cp02175h#!divAbstract 2]. At present [http://openmopac.net/ MOPAC2016], [https://www.entos.ai/qcore/documentation/ Entos Qcore] and Gaussian 09 (G09) are interfaced with AutoMeKin. The program has been tested on the following Linux distros: CentOS 7, Red Hat Enterprise Linux and Ubuntu 20.04 LTS.
+
[http://pubs.rsc.org/en/content/articlelanding/2015/cp/c5cp02175h#!divAbstract 2]. At present [http://openmopac.net/ MOPAC2016], [https://www.entos.ai/qcore/documentation/ Entos Qcore] and Gaussian (G09/G16) are interfaced with AutoMeKin. The program has been tested on the following Linux distros: CentOS 7, Red Hat Enterprise Linux and Ubuntu 20.04 LTS.
  
 
To give you a flavor of the capabilities of the program you can try our [http://rxnkin.usc.es/amk/ web interface]
 
To give you a flavor of the capabilities of the program you can try our [http://rxnkin.usc.es/amk/ web interface]
Line 75: Line 75:
 
<ol start="1">
 
<ol start="1">
 
<li>[https://onlinelibrary.wiley.com/doi/full/10.1002/jcc.26734 <span style="font-size:100%">E Martínez-Nuñez, G.L. Barnes, D.R. Glowacki, S. Kopec, D. Pelaez, A. Rodriguez, R. Rodriguez-Fernandez, R.J. Shannon, J.J.P. Stewart, P.G. Tahoces, S.A. Vazquez, J. Comput. Chem. 2021, 42, 2036</span>]</li>
 
<li>[https://onlinelibrary.wiley.com/doi/full/10.1002/jcc.26734 <span style="font-size:100%">E Martínez-Nuñez, G.L. Barnes, D.R. Glowacki, S. Kopec, D. Pelaez, A. Rodriguez, R. Rodriguez-Fernandez, R.J. Shannon, J.J.P. Stewart, P.G. Tahoces, S.A. Vazquez, J. Comput. Chem. 2021, 42, 2036</span>]</li>
<li>[http://onlinelibrary.wiley.com/doi/10.1002/jcc.23790/abstract <span style="font-size:100%">E. Martínez-Núñez J. Comput. Chem. 2015, 36, 222</span>]</li>
+
<li>[http://onlinelibrary.wiley.com/doi/10.1002/jcc.23790/abstract <span style="font-size:100%">E. Martínez-Núñez, J. Comput. Chem. 2015, 36, 222</span>]</li>
<li>[http://pubs.rsc.org/en/content/articlelanding/2015/cp/c5cp02175h#!divAbstract <span style="font-size:100%">E. Martínez-Núñez Phys. Chem. Chem. Phys. 2015,17, 14912</span>]</li>
+
<li>[http://pubs.rsc.org/en/content/articlelanding/2015/cp/c5cp02175h#!divAbstract <span style="font-size:100%">E. Martínez-Núñez, Phys. Chem. Chem. Phys. 2015,17, 14912</span>]</li>
 
<li>[http://OpenMOPAC.net <span style="font-size:100%">J. J. P. Stewart, MOPAC2016, Stewart Computational Chemistry: Colorado Springs, CO, USA, HTTP://OpenMOPAC.net, 2016.</span>]</li>
 
<li>[http://OpenMOPAC.net <span style="font-size:100%">J. J. P. Stewart, MOPAC2016, Stewart Computational Chemistry: Colorado Springs, CO, USA, HTTP://OpenMOPAC.net, 2016.</span>]</li>
 
</ol>
 
</ol>
Line 95: Line 95:
  
 
==News==
 
==News==
The Python library [https://github.com/dgarayr/amk_tools '''amk-tools'''] can be employed to read, process and visualize reaction networks generated by AutoMeKin and is available in the latest singularity container, available through the [[Installation instructions#Auto installer| autoinstaller]].
+
AutoMeKin is now interfaced with Gaussian 16 (from rev. 1131)
 
 
html lang="en">
 
 
 
  <head>
 
   
 
      <meta charset="utf-8">
 
      <title>Reaction network visualization</title>
 
     
 
     
 
       
 
         
 
       
 
       
 
         
 
        <script type="text/javascript" src="https://cdn.bokeh.org/bokeh/release/bokeh-2.3.3.min.js" integrity="sha384-dM3QQsP+wXdHg42wTqW85BjZQdLNNIXqlPw/BgKoExPmTG7ZLML4EGqLMfqHT6ON" crossorigin="anonymous"></script>
 
        <script type="text/javascript" src="https://cdn.bokeh.org/bokeh/release/bokeh-widgets-2.3.3.min.js" integrity="sha384-3QTqdz9LyAm2i0sG5XTePsHec3UHWwVsrOL68SYRoAXsafvfAyqtQ+h440+qIBhS" crossorigin="anonymous"></script>
 
        <script type="text/javascript">
 
            Bokeh.set_log_level("info");
 
        </script>
 
        <script type="text/javascript">
 
            (function(root, factory) {
 
                    factory(root["Bokeh"]);
 
                })(this, function(Bokeh) {
 
                  var define;
 
                  return (function outer(modules, entry) {
 
                  if (Bokeh != null) {
 
                    return Bokeh.register_plugin(modules, entry);
 
                  } else {
 
                    throw new Error("Cannot find Bokeh. You have to load it prior to loading plugins.");
 
                  }
 
                })
 
                ({
 
                  "custom/main": function(require, module, exports) {
 
                    var models = {
 
                      "JSMol": require("custom/jsmol_bokeh_extension.jsmol.js_mol").JSMol
 
                    };
 
                    require("base").register_models(models);
 
                    module.exports = models;
 
                  },
 
                  "custom/jsmol_bokeh_extension.jsmol.js_mol": function(require, module, exports) {
 
                "use strict";
 
                var __createBinding = (this && this.__createBinding) || (Object.create ? (function(o, m, k, k2) {
 
                    if (k2 === undefined) k2 = k;
 
                    Object.defineProperty(o, k2, { enumerable: true, get: function() { return m[k]; } });
 
                }) : (function(o, m, k, k2) {
 
                    if (k2 === undefined) k2 = k;
 
                    o[k2] = m[k];
 
                }));
 
                var __setModuleDefault = (this && this.__setModuleDefault) || (Object.create ? (function(o, v) {
 
                    Object.defineProperty(o, "default", { enumerable: true, value: v });
 
                }) : function(o, v) {
 
                    o["default"] = v;
 
                });
 
                var __importStar = (this && this.__importStar) || function (mod) {
 
                    if (mod && mod.__esModule) return mod;
 
                    var result = {};
 
                    if (mod != null) for (var k in mod) if (k !== "default" && Object.prototype.hasOwnProperty.call(mod, k)) __createBinding(result, mod, k);
 
                    __setModuleDefault(result, mod);
 
                    return result;
 
                };
 
                Object.defineProperty(exports, "__esModule", { value: true });
 
                exports.JSMol = exports.JSMolView = void 0;
 
                // These "require" lines are similar to python "import" statements
 
                const layout_dom_1 = require("models/layouts/layout_dom");
 
                const layout_1 = require("core/layout");
 
                const p = __importStar(require("core/properties"));
 
                // This defines some default options for JSmol
 
                // See https://gist.github.com/jhjensen2/4701339 for more details.
 
                const INFO = {
 
                    height: "100%",
 
                    width: "100%",
 
                    serverURL: "https://chemapps.stolaf.edu/jmol/jsmol/php/jsmol.php",
 
                    use: "HTML5",
 
                    j2sPath: "https://chemapps.stolaf.edu/jmol/jsmol/j2s",
 
                    script: "background black;load https://dev-www.materialscloud.org/cofs/api/v2/cifs/febd2d02-5690-4a07-9013-505c9a06bc5b/content/download",
 
                };
 
                // To create custom model extensions that will render on to the HTML canvas
 
                // or into the DOM, we must create a View subclass for the model.
 
                //
 
                // In this case we will subclass from the existing BokehJS ``LayoutDOMView``
 
                class JSMolView extends layout_dom_1.LayoutDOMView {
 
                    initialize() {
 
                        super.initialize();
 
                        var url = this.model.js_url;
 
                        if (!url) {
 
                            url = "https://chemapps.stolaf.edu/jmol/jsmol/JSmol.min.js";
 
                        }
 
                        const script = document.createElement('script');
 
                        script.src = url;
 
                        script.async = false;
 
                        script.onload = () => this._init();
 
                        //script.onreadystatechange = (script.onload = () => this._init())
 
                        //script.onreadystatechange = script.onload
 
                        document.querySelector("head").appendChild(script);
 
                    }
 
                    _init() {
 
                        // Create a new Jmol applet using the JSmol.js API. This assumes JSmol.js
 
                        // has already been loaded.
 
                        // Models will be able to specify and load external scripts automatically.
 
                        //
 
                        // BokehJS Views create <div> elements by default, accessible as this.el.
 
                        // if no info dict provided, use reasonable default
 
                        if (!this.model.info) {
 
                            this.model.info = INFO;
 
                        }
 
                        // disable usage tracker - this conflicts with cross-site-scripting policies when served over https
 
                        delete Jmol._tracker;
 
                        // returns html + assigns applet object to "jmolApplet0" variable
 
                        var html = Jmol.getAppletHtml("jmolApplet0", this.model.info);
 
                        this._applet = jmolApplet0;
 
                        this.el.innerHTML = html;
 
                        // avoid creating a "deferred" applet
 
                        // https://sourceforge.net/p/jsmol/discussion/general/thread/48083aa7///10a6/bc1c
 
                        this._applet._cover(false);
 
                        // Set a listener so that when the Bokeh script input changes it is executed
 
                        this.connect(this.model.script_source.properties.data.change, () => {
 
                            console.log("Script source changed");
 
                            Jmol.script(this._applet, this.model.script_source.get_column('script')[0]);
 
                        });
 
                    }
 
                    get child_models() {
 
                        return [];
 
                    }
 
                    _update_layout() {
 
                        this.layout = new layout_1.LayoutItem();
 
                        this.layout.set_sizing(this.box_sizing());
 
                    }
 
                }
 
                exports.JSMolView = JSMolView;
 
                JSMolView.__name__ = "JSMolView";
 
                class JSMol extends layout_dom_1.LayoutDOM {
 
                    constructor(attrs) {
 
                        super(attrs);
 
                    }
 
                    static initClass() {
 
                        // The ``type`` class attribute should generally match exactly the name
 
                        // of the corresponding Python class.
 
                        this.prototype.type = "JSMol";
 
                        // This is usually boilerplate. In some cases there may not be a view.
 
                        this.prototype.default_view = JSMolView;
 
                        // The define block adds corresponding "properties" to the JS model. These
 
                        // should basically line up 1-1 with the Python model class. Most property
 
                        // types have counterparts, e.g. ``bokeh.core.properties.String`` will be
 
                        // ``p.String`` in the JS implementation. Where the JS type system is not yet
 
                        // as rich, you can use ``p.Any`` as a "wildcard" property type.
 
                        this.define({
 
                            script_source: [p.Instance],
 
                            info: [p.Any],
 
                            js_url: [p.String],
 
                        });
 
                    }
 
                }
 
                exports.JSMol = JSMol;
 
                JSMol.__name__ = "JSMol";
 
                JSMol.initClass();
 
                //# sourceMappingURL=jsmol.py:JSMol.js.map
 
                }
 
                }, "custom/main");
 
                ;
 
                });
 
 
 
        </script>
 
       
 
     
 
     
 
<style>
 
.bk-root .bk-btn-default {
 
font-size: 1.2vh;
 
}
 
.bk-root .bk-input {
 
font-size: 1.2vh;
 
padding-bottom: 5px;
 
padding-top: 5px;
 
}
 
.bk-root .bk {
 
font-size: 1.2vh;
 
}
 
.bk-root .bk-clearfix{
 
padding-bottom: 0.8vh;
 
}
 
</style>
 
 
   
 
  </head>
 
 
 
 
 
  <body>
 
   
 
     
 
       
 
         
 
         
 
           
 
              <div class="bk-root" id="3c780241-42fe-478b-8b7d-bb4e3e9dfbe4" data-root-id="1159"></div>
 
           
 
         
 
       
 
     
 
     
 
        <script type="application/json" id="1703">
 
          {"f436dc99-c491-4996-b7ca-3e1293247474":{"defs":[],"roots":{"references":[{"attributes":{"formatter":{"id":"1250"},"major_label_policy":{"id":"1249"},"ticker":{"id":"1171"},"visible":false},"id":"1170","type":"LinearAxis"},{"attributes":{"axis_label":"E (kcal/mol)","axis_label_text_font":"Arial","axis_label_text_font_size":"1.5vh","axis_label_text_font_style":"bold","formatter":{"id":"1253"},"major_label_policy":{"id":"1252"},"major_label_text_font_size":"1.2vh","ticker":{"id":"1175"}},"id":"1174","type":"LinearAxis"},{"attributes":{},"id":"1282","type":"UnionRenderers"},{"attributes":{"data":{"elab":["-8.3","58.2","-37.6"],"lab":["MIN1","TS10","PR2"],"x":{"__ndarray__":"AAAAAAAA4D8AAAAAAAAEQAAAAAAAABJA","dtype":"float64","order":"little","shape":[3]},"y":[-8.323520999999985,58.154156,-37.60977]},"selected":{"id":"1275"},"selection_policy":{"id":"1276"}},"id":"1216","type":"ColumnDataSource"},{"attributes":{},"id":"1171","type":"BasicTicker"},{"attributes":{"axis":{"id":"1170"},"ticker":null},"id":"1173","type":"Grid"},{"attributes":{"text":"Reaction network visualization","text_font_size":"1.5vh"},"id":"1005","type":"Title"},{"attributes":{},"id":"1256","type":"Selection"},{"attributes":{"axis":{"id":"1174"},"dimension":1,"ticker":null},"id":"1177","type":"Grid"},{"attributes":{},"id":"1257","type":"UnionRenderers"},{"attributes":{},"id":"1269","type":"Selection"},{"attributes":{},"id":"1175","type":"BasicTicker"},{"attributes":{},"id":"1270","type":"UnionRenderers"},{"attributes":{},"id":"1258","type":"Selection"},{"attributes":{},"id":"1179","type":"WheelZoomTool"},{"attributes":{},"id":"1178","type":"PanTool"},{"attributes":{},"id":"1259","type":"UnionRenderers"},{"attributes":{"overlay":{"id":"1183"}},"id":"1180","type":"BoxZoomTool"},{"attributes":{},"id":"1181","type":"ResetTool"},{"attributes":{},"id":"1182","type":"SaveTool"},{"attributes":{"args":{"prof":{"id":"1161"},"thrbox":{"id":"1137"}},"code":"\n\t\t// graph - bokeh graph for the network view, use to fetch node and edge renderers\n\t\t// prof - bokeh figure for the profiles, thrbox - textbox for threshold energy\n\t\tvar thr = parseFloat(thrbox.value)\n\t\tvar range = {min: thr, max: thr}\n\t\tfunction overThreshold(element){\n\t\t\treturn element &gt; thr\n\t\t}\n\t\tvar Nlines = prof.renderers.length/2\n\t\tfor (let i = 0 ; i &lt; Nlines ; i++){\n\t\t\tconst index = 2*i\n\t\t\tvar current_data = prof.renderers[index].data_source\n\t\t\tvar energies = current_data.data[\"y\"]\n\t\t\t//hide if anything is above the threshold\n\t\t\tif (energies.some(overThreshold)){\n\t\t\t\tprof.renderers[index].visible = false\n\t\t\t\tprof.renderers[index+1].visible = false\n\t\t\t\t// labels are also duplicated: first tag, then energy.\n\t\t\t\tprof.center[index+2].visible = false\n\t\t\t\tprof.center[index+3].visible = false\n\t\t\t} \n\t\t}\n\t\t"},"id":"1228","type":"CustomJS"},{"attributes":{"args":{"prof":{"id":"1161"},"state":[false]},"code":"\n\t// show ENERGY labels for the profile\n\t// prof - bokeh profile for figures, state - counter for the inner state\n\tstate[0] = !state[0]\n\tvar Nlines = prof.renderers.length/2\n\tfor (let i = 0 ; i &lt; Nlines ; i++){\n\t\t// show or hide for ODD ii values (energy label) AND only for visible renderers\n\t\tvar jj = 2*i + 1\n\t\tif (prof.renderers[jj-1].visible == true){\n\t\t\tprof.center[jj+2].visible = state[0]\n\t\t} \n\t}\n\t"},"id":"1229","type":"CustomJS"},{"attributes":{"align":"center","height":100,"text":"(Click an element)"},"id":"1133","type":"Div"},{"attributes":{},"id":"1271","type":"Selection"},{"attributes":{},"id":"1272","type":"UnionRenderers"},{"attributes":{"callback":{"id":"1229"},"description":"Show energy labels","icon":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAAgAAAAICAIAAABLbSncAAAAKklEQVR4nGP8//8/AzbAwsDAcPLkSTRRc3NzFjgLTY4JqzlQozBNo65RAAfvDWESs4jlAAAAAElFTkSuQmCC"},"id":"1230","type":"CustomAction"},{"attributes":{},"id":"1273","type":"Selection"},{"attributes":{"align":"center","disabled":true,"icon":null,"js_event_callbacks":{"button_click":[{"id":"1228"}]},"label":"Energy filter","max_width":120},"id":"1136","type":"Button"},{"attributes":{},"id":"1252","type":"AllLabels"},{"attributes":{},"id":"1253","type":"BasicTickFormatter"},{"attributes":{},"id":"1274","type":"UnionRenderers"},{"attributes":{},"id":"1260","type":"Selection"},{"attributes":{"js_property_callbacks":{"change:indices":[{"id":"1138"}]}},"id":"1147","type":"Selection"},{"attributes":{"align":"center","disabled":true,"icon":null,"js_event_callbacks":{"menu_item_click":[{"id":"1142"}]},"label":"Normal modes","max_width":400,"menu":[["No vibr. loaded","None"],null]},"id":"1131","type":"Dropdown"},{"attributes":{"align":"center","disabled":true,"icon":null,"js_event_callbacks":{"button_click":[{"id":"1227"}]},"label":"Molec. filter","max_width":120},"id":"1134","type":"Button"},{"attributes":{},"id":"1261","type":"UnionRenderers"},{"attributes":{},"id":"1275","type":"Selection"},{"attributes":{"align":"center","width":300},"id":"1132","type":"Spacer"},{"attributes":{"name":"ENERGYLABELS","source":{"id":"1219"},"text":{"field":"elab"},"text_font_size":{"value":"1.4vh"},"visible":false,"x":{"field":"x"},"x_offset":{"value":-1},"y":{"field":"y"},"y_offset":{"value":-20}},"id":"1220","type":"LabelSet"},{"attributes":{},"id":"1276","type":"UnionRenderers"},{"attributes":{"args":{"graph":{"id":"1033"},"source":{"id":"1125"},"textField":{"id":"1133"}},"code":"\n\t\t// from graph, we fetch nrend - node renderer, erend - edgerenderer\n\t\t//source - source object for JSMol\n\t\tvar nrend = graph.node_renderer.data_source\n\t\tvar erend = graph.edge_renderer.data_source\n\t\tif (cb_obj.indices.length) {\n\t\t\tvar ninds = nrend.selected.indices\n\t\t\tvar einds = erend.selected.indices\n\t\t\tif (ninds.length) {\n\t\t\t\tvar rend = nrend\n\t\t\t} else {\n\t\t\t\tvar rend = erend\n\t\t\t}\n\t\t\tvar ndx = rend.selected.indices[0]\n\t\t\tvar model = rend.data[\"model\"][ndx]\n\t\t\tif (model == null){\n\t\t\t\t// clear view if no model is available\n\t\t\t\tmodel = \"backbone only ; backbone off\"\n\t\t\t}\n\t\t\tvar molname = rend.data[\"name\"][ndx]\n\t\t\ttextField.text = \"&lt;font size=+1&gt;&lt;b&gt;\" + molname + \"&lt;/b&gt;&lt;/font&gt;\"\n\t\t\tsource.data[\"script\"] = [model]\n\t\t\tsource.properties.data.change.emit()\n\t\t}\n\t\t"},"id":"1138","type":"CustomJS"},{"attributes":{"data":{"elab":["-8.3","58.2","-37.6"],"lab":["MIN1","TS10","PR2"],"x":{"__ndarray__":"AAAAAAAA4D8AAAAAAAAEQAAAAAAAABJA","dtype":"float64","order":"little","shape":[3]},"y":[-8.323520999999985,58.154156,-37.60977]},"selected":{"id":"1277"},"selection_policy":{"id":"1278"}},"id":"1219","type":"ColumnDataSource"},{"attributes":{"align":"center","value":"58.15","width":240},"id":"1137","type":"TextInput"},{"attributes":{"end":1.2,"start":-1.2},"id":"1003","type":"Range1d"},{"attributes":{},"id":"1277","type":"Selection"},{"attributes":{"align":"center","js_property_callbacks":{"change:active":[{"id":"1226"}]},"labels":["Show profile"],"max_height":100,"max_width":120},"id":"1135","type":"CheckboxGroup"},{"attributes":{"below":[{"id":"1013"}],"center":[{"id":"1016"},{"id":"1020"},{"id":"1082"},{"id":"1085"}],"left":[{"id":"1017"}],"renderers":[{"id":"1033"}],"title":{"id":"1005"},"toolbar":{"id":"1027"},"toolbar_location":"above","width":800,"x_range":{"id":"1002"},"x_scale":{"id":"1009"},"y_range":{"id":"1003"},"y_scale":{"id":"1011"}},"id":"1004","subtype":"Figure","type":"Plot"},{"attributes":{},"id":"1278","type":"UnionRenderers"},{"attributes":{},"id":"1022","type":"WheelZoomTool"},{"attributes":{"data":{"script":[]},"selected":{"id":"1260"},"selection_policy":{"id":"1261"}},"id":"1125","type":"ColumnDataSource"},{"attributes":{},"id":"1021","type":"PanTool"},{"attributes":{"args":{"graph":{"id":"1033"},"menu":{"id":"1131"},"source":{"id":"1125"},"textField":{"id":"1133"}},"code":"\n\t\t// from graph, we fetch nrend - node renderer, erend - edgerenderer, \n\t\t// source - source object for JSMol, menu - dropdown menu to change vibrations\n\t\tvar nrend = graph.node_renderer.data_source\n\t\tvar erend = graph.edge_renderer.data_source\n\t\t\n\t\tvar ninds = nrend.selected.indices\n\t\tvar einds = erend.selected.indices\n\t\tif (ninds.length) {\n\t\t\tvar rend = nrend\n\t\t} else {\n\t\t\tvar rend = erend\n\t\t}\n\t\tvar ndx = rend.selected.indices[0]\n\t\tvar vibr_models = rend.data[\"vibr_models\"][ndx]\n\t\tif (vibr_models.length){\n\t\t\tsource.data[\"current_model\"] = [vibr_models]\n\t\t\tsource.data[\"script\"] = [vibr_models + \"; vibration on ; vibration scale 0.25\"]\n\t\t\tsource.data[\"current_mol\"] = rend.data[\"name\"][ndx]\n\t\t\t// and also populate the normal mode dropdown: get the frequencies and assign them to menu\n\t\t\t// we must pass indices as strings for it to work, add 1 to match model numbering\n\t\t\tmenu.disabled = false\n\t\t\tvar freqlist = rend.data[\"frequencies\"][ndx].split(\";\")\n\t\t\tvar menulist = freqlist.map(function (frq,ii) {return [frq+\" cm-1\",(ii+1).toString()]})\n\t\t\tmenu.menu = menulist\n\t\t\tvar nw_label = \"Normal modes (\" + source.data[\"current_mol\"] + \")\"\n\t\t\tmenu.label = nw_label\n\t\t\t// modify text element using the BASE frequency\n\t\t\tvar molname = source.data[\"current_mol\"]\n\t\t\tvar frqval = freqlist[0].toString().trim() + \" cm-1\"\n\t\t\ttextField.text = \"&lt;font size=+1&gt;&lt;b&gt;\" + molname + \" (\" + frqval + \")\" + \"&lt;/b&gt;&lt;/font&gt;\"\n\t\t}\n\t\n\t\tsource.properties.data.change.emit()\n\t\t"},"id":"1139","type":"CustomJS"},{"attributes":{"js_property_callbacks":{"change:indices":[{"id":"1138"}]}},"id":"1152","type":"Selection"},{"attributes":{"args":{"button":{"id":"1128"},"graph":{"id":"1033"},"source":{"id":"1125"}},"code":"\n\t\t// from graph, we fetch nrend - node renderer, erend - edgerenderer\n\t\t// source - source object for JSMol, button - button object\n\t\tvar nrend = graph.node_renderer.data_source\n\t\tvar erend = graph.edge_renderer.data_source\n\n\t\tvar ninds = nrend.selected.indices\n\t\tvar einds = erend.selected.indices\n\t\tif (ninds.length) {\n\t\t\tvar rend = nrend\n\t\t} else {\n\t\t\tvar rend = erend\n\t\t}\n\t\tvar ndx = rend.selected.indices[0]\n\t\tvar geo = rend.data[\"geometry\"][ndx]\n\t\tvar name = rend.data[\"name\"][ndx]\n\t\tvar orig_label = button.label\n\t\t// use a function to apply setTimeout and recover original label\n\t\tfunction recover_label(button,orig_label){\n\t\t\tbutton.label = orig_label\n\t\t}\t\n\n\t\tif (geo.length){\n\t\t\tnavigator.clipboard.writeText([geo])\n\t\t\tbutton.label = \"Copied \" + name + \"!\"\n\t\t\tsetTimeout(recover_label,2000,button,orig_label)\n\t\t}\n\t\t"},"id":"1140","type":"CustomJS"},{"attributes":{"data_source":{"id":"1039"},"glyph":{"id":"1038"},"hover_glyph":{"id":"1071"},"muted_glyph":null,"selection_glyph":{"id":"1076"},"view":{"id":"1041"}},"id":"1040","type":"GlyphRenderer"},{"attributes":{"args":{"fig":{"id":"1004"},"graph":{"id":"1033"},"text_input":{"id":"1129"}},"code":"\n\t\t// source - source object for JSMol\n\t\t// pass graph and fetch node and edge renderers\n\t\t// from fig, we modify x_range and y_range. Default plot starts from -1.2 to 1.2,\n\t\tvar nrend = graph.node_renderer.data_source\n\t\tvar erend = graph.edge_renderer.data_source\n\t\tvar layout = graph.layout_provider.graph_layout\n\t\t// fetch the query in the data sources, choosing the appropiate renderer depending on the query\n\t\tvar mol_query = text_input.value\n\t\tif (mol_query.includes(\"TS\")) {\n\t\t\tvar renderer = erend\n\t\t\tvar other_renderer = nrend\n\t\t} else {\n\t\t\tvar renderer = nrend\n\t\t\tvar other_renderer = erend\n\t\t}\n\t\tvar pool_names = renderer.data[\"name\"]\n\t\tvar ndx = pool_names.indexOf(mol_query)\n\t\t// locate positions of the node or of the nodes defining an edge\n\t\tif (mol_query.includes(\"TS\")) {\n\t\t\tvar n1 = renderer.data[\"start\"][ndx]\n\t\t\tvar n2 = renderer.data[\"end\"][ndx]\n\t\t\tvar pos1 = layout[n1]\n\t\t\tvar pos2 = layout[n2]\n\t\t\tvar positions = new Array(2)\n\t\t\tpositions[0] = 0.5*(pos1[0]+pos2[0])\n\t\t\tpositions[1] = 0.5*(pos1[1]+pos2[1])\n\t\t} else {\n\t\t\tvar positions = layout[mol_query]\n\t\t}\n\t\t// returns -1 if element is not present\n\t\tif (ndx &gt;= 0) {\n\t\t\t// clearing other sel. avoids problems for model loading sometimes\n\t\t\tother_renderer.selected.indices = []\n\t\t\trenderer.selected.indices = [ndx]\n\t\t\tfig.x_range.start = positions[0] - 0.5\n\t\t\tfig.x_range.end = positions[0] + 0.5\n\t\t\tfig.y_range.start = positions[1] - 0.5\n\t\t\tfig.y_range.end = positions[1] + 0.5\n\t\t}\n\t\t"},"id":"1141","type":"CustomJS"},{"attributes":{"data":{"form":[null,null,null,null,null,null,null,null,"H2 + CO2","H2 + CO2"],"lab":["MIN1","MIN1","TS9","TS9","MIN3","MIN3","TS8","TS8","PR2","PR2"],"x":{"__ndarray__":"AAAAAAAAAAAAAAAAAADwPwAAAAAAAABAAAAAAAAACEAAAAAAAAAQQAAAAAAAABRAAAAAAAAAGEAAAAAAAAAcQAAAAAAAACBAAAAAAAAAIkA=","dtype":"float64","order":"little","shape":[10]},"y":[-8.323520999999985,-8.323520999999985,53.170985,53.170985,5.312669999999997,5.312669999999997,43.967468999999994,43.967468999999994,-37.60977,-37.60977]},"selected":{"id":"1263"},"selection_policy":{"id":"1264"}},"id":"1190","type":"ColumnDataSource"},{"attributes":{"args":{"menu":{"id":"1131"},"source":{"id":"1125"},"textField":{"id":"1133"}},"code":"\n\t\t// source - source object for JSMol ; menu - dropdown menu object\n\t\tsource.data[\"script\"] = [source.data[\"current_model\"] + \"; vibration on ; vibration scale 0.25 ; model \" + this.item]\n\t\tvar ndx = parseInt(this.item) - 1\n\t\tvar molname = source.data[\"current_mol\"]\n\t\tvar frqval = menu.menu[ndx][0]\n\t\tsource.data[\"current_vibration\"] = frqval\n\t\ttextField.text = \"&lt;font size=+1&gt;&lt;b&gt;\" + molname + \" (\" + frqval.trim() + \")\" + \"&lt;/b&gt;&lt;/font&gt;\"\n\t\tsource.properties.data.change.emit()\n\t\t"},"id":"1142","type":"CustomJS"},{"attributes":{},"id":"1263","type":"Selection"},{"attributes":{"data_source":{"id":"1212"},"glyph":{"id":"1213"},"hover_glyph":null,"muted_glyph":null,"view":{"id":"1215"}},"id":"1214","type":"GlyphRenderer"},{"attributes":{"axis":{"id":"1013"},"grid_line_color":null,"ticker":null},"id":"1016","type":"Grid"},{"attributes":{"source":{"id":"1212"}},"id":"1215","type":"CDSView"},{"attributes":{},"id":"1264","type":"UnionRenderers"},{"attributes":{"name":"THELABELS","source":{"id":"1216"},"text":{"field":"lab"},"text_font_size":{"value":"1.4vh"},"x":{"field":"x"},"x_offset":{"value":-1},"y":{"field":"y"},"y_offset":{"value":2}},"id":"1217","type":"LabelSet"},{"attributes":{"formatter":{"id":"1247"},"major_label_policy":{"id":"1246"},"ticker":{"id":"1018"},"visible":false},"id":"1017","type":"LinearAxis"},{"attributes":{"fill_alpha":{"value":0.9},"fill_color":{"value":"#ff7f0e"},"height":{"value":2},"line_width":{"value":0},"name":"RECTANGLE","width":{"value":1},"x":{"field":"x"},"y":{"field":"y"}},"id":"1213","type":"Rect"},{"attributes":{},"id":"1011","type":"LinearScale"},{"attributes":{"children":[{"id":"1129"},{"id":"1130"}],"height":100},"id":"1154","type":"Row"},{"attributes":{"children":[{"id":"1132"},{"id":"1133"}],"height":100},"id":"1156","type":"Row"},{"attributes":{"axis":{"id":"1017"},"dimension":1,"grid_line_color":null,"ticker":null},"id":"1020","type":"Grid"},{"attributes":{},"id":"1018","type":"BasicTicker"},{"attributes":{"children":[{"id":"1126"}]},"id":"1157","type":"Row"},{"attributes":{"overlay":{"id":"1026"}},"id":"1023","type":"BoxZoomTool"},{"attributes":{"edge_renderer":{"id":"1040"},"inspection_policy":{"id":"1087"},"layout_provider":{"id":"1046"},"node_renderer":{"id":"1036"},"selection_policy":{"id":"1102"}},"id":"1033","type":"GraphRenderer"},{"attributes":{},"id":"1024","type":"ResetTool"},{"attributes":{"children":[{"id":"1127"},{"id":"1131"},{"id":"1128"}],"height":100},"id":"1153","type":"Row"},{"attributes":{},"id":"1025","type":"SaveTool"},{"attributes":{"children":[{"id":"1153"},{"id":"1004"},{"id":"1154"}],"width":800},"id":"1155","type":"Column"},{"attributes":{"high":5.312669999999997,"low":-37.60977,"palette":["#440154","#440255","#440357","#450558","#45065A","#45085B","#46095C","#460B5E","#460C5F","#460E61","#470F62","#471163","#471265","#471466","#471567","#471669","#47186A","#48196B","#481A6C","#481C6E","#481D6F","#481E70","#482071","#482172","#482273","#482374","#472575","#472676","#472777","#472878","#472A79","#472B7A","#472C7B","#462D7C","#462F7C","#46307D","#46317E","#45327F","#45347F","#453580","#453681","#443781","#443982","#433A83","#433B83","#433C84","#423D84","#423E85","#424085","#414186","#414286","#404387","#404487","#3F4587","#3F4788","#3E4888","#3E4989","#3D4A89","#3D4B89","#3D4C89","#3C4D8A","#3C4E8A","#3B508A","#3B518A","#3A528B","#3A538B","#39548B","#39558B","#38568B","#38578C","#37588C","#37598C","#365A8C","#365B8C","#355C8C","#355D8C","#345E8D","#345F8D","#33608D","#33618D","#32628D","#32638D","#31648D","#31658D","#31668D","#30678D","#30688D","#2F698D","#2F6A8D","#2E6B8E","#2E6C8E","#2E6D8E","#2D6E8E","#2D6F8E","#2C708E","#2C718E","#2C728E","#2B738E","#2B748E","#2A758E","#2A768E","#2A778E","#29788E","#29798E","#287A8E","#287A8E","#287B8E","#277C8E","#277D8E","#277E8E","#267F8E","#26808E","#26818E","#25828E","#25838D","#24848D","#24858D","#24868D","#23878D","#23888D","#23898D","#22898D","#228A8D","#228B8D","#218C8D","#218D8C","#218E8C","#208F8C","#20908C","#20918C","#1F928C","#1F938B","#1F948B","#1F958B","#1F968B","#1E978A","#1E988A","#1E998A","#1E998A","#1E9A89","#1E9B89","#1E9C89","#1E9D88","#1E9E88","#1E9F88","#1EA087","#1FA187","#1FA286","#1FA386","#20A485","#20A585","#21A685","#21A784","#22A784","#23A883","#23A982","#24AA82","#25AB81","#26AC81","#27AD80","#28AE7F","#29AF7F","#2AB07E","#2BB17D","#2CB17D","#2EB27C","#2FB37B","#30B47A","#32B57A","#33B679","#35B778","#36B877","#38B976","#39B976","#3BBA75","#3DBB74","#3EBC73","#40BD72","#42BE71","#44BE70","#45BF6F","#47C06E","#49C16D","#4BC26C","#4DC26B","#4FC369","#51C468","#53C567","#55C666","#57C665","#59C764","#5BC862","#5EC961","#60C960","#62CA5F","#64CB5D","#67CC5C","#69CC5B","#6BCD59","#6DCE58","#70CE56","#72CF55","#74D054","#77D052","#79D151","#7CD24F","#7ED24E","#81D34C","#83D34B","#86D449","#88D547","#8BD546","#8DD644","#90D643","#92D741","#95D73F","#97D83E","#9AD83C","#9DD93A","#9FD938","#A2DA37","#A5DA35","#A7DB33","#AADB32","#ADDC30","#AFDC2E","#B2DD2C","#B5DD2B","#B7DD29","#BADE27","#BDDE26","#BFDF24","#C2DF22","#C5DF21","#C7E01F","#CAE01E","#CDE01D","#CFE11C","#D2E11B","#D4E11A","#D7E219","#DAE218","#DCE218","#DFE318","#E1E318","#E4E318","#E7E419","#E9E419","#ECE41A","#EEE51B","#F1E51C","#F3E51E","#F6E61F","#F8E621","#FAE622","#FDE724"]},"id":"1051","type":"LinearColorMapper"},{"attributes":{"children":[{"id":"1156"},{"id":"1157"},{"id":"1161"},{"id":"1160"}],"width":600},"id":"1158","type":"Column"},{"attributes":{},"id":"1265","type":"Selection"},{"attributes":{"source":{"id":"1039"}},"id":"1041","type":"CDSView"},{"attributes":{},"id":"1266","type":"UnionRenderers"},{"attributes":{"graph_layout":{"MIN1":[-0.50000005039689,0.8660253835545058],"MIN3":[0.9999999999999999,4.7706792907989705e-08],"PR2":[-0.49999994960310984,-0.8660254312612986]}},"id":"1046","type":"StaticLayoutProvider"},{"attributes":{"line_alpha":0.1,"line_color":"#ff7f0e","x":{"field":"x"},"y":{"field":"y"}},"id":"1209","type":"Line"},{"attributes":{},"id":"1280","type":"UnionRenderers"},{"attributes":{"children":[{"id":"1135"},{"id":"1134"},{"id":"1137"},{"id":"1136"}],"height":100,"height_policy":"max"},"id":"1160","type":"Row"},{"attributes":{},"id":"1234","type":"Title"},{"attributes":{"source":{"id":"1035"}},"id":"1037","type":"CDSView"},{"attributes":{"data_source":{"id":"1035"},"glyph":{"id":"1052"},"hover_glyph":{"id":"1062"},"muted_glyph":null,"selection_glyph":{"id":"1057"},"view":{"id":"1037"}},"id":"1036","type":"GlyphRenderer"},{"attributes":{"line_color":"#ff7f0e","x":{"field":"x"},"y":{"field":"y"}},"id":"1208","type":"Line"},{"attributes":{"data":{"barrless":[false,false,false],"end":["PR2","MIN1","PR2"],"energy":[43.967468999999994,53.170985,58.154156],"fname":["TS8_ts8_batch9.rxyz","TS9_ts10_batch3.rxyz","TS10_ts3_batch4.rxyz"],"frequencies":["-1949.3;777.7;861.4;895.4;1136.4;1577.9;1608.9;2040.0;2711.8","-1969.9;354.7;496.2;641.8;1064.4;1342.0;1552.5;2462.8;2757.4","-2928.3;451.4;587.2;904.5;997.8;1345.8;1946.9;2103.1;2635.0"],"geometry":["C -0.3115 -0.469 0\nO 0.935099 -0.469 0\nO -0.794898 0.680099 0\nH 1.1494 0.850298 0\nH 0.3382 1.3891 0","C -0.2036 -0.3962 0\nO 1.0655 -0.3962 0\nO -0.911899 0.694699 0\nH 0.4048 -1.5715 0\nH -0.4169 1.5559 0.0001","C 0.008 -0.221 0\nO 1.1831 -0.221 0\nO -1.0411 0.5156 0\nH -1.5835 -0.698202 0\nH -0.764698 -1.3442 0"],"model":["data 'model X'|5|TS8|C -0.3115 -0.469 0|O 0.935099 -0.469 0|O -0.794898 0.680099 0|H 1.1494 0.850298 0|H 0.3382 1.3891 0|end 'model X';show data","data 'model X'|5|TS9|C -0.2036 -0.3962 0|O 1.0655 -0.3962 0|O -0.911899 0.694699 0|H 0.4048 -1.5715 0|H -0.4169 1.5559 0.0001|end 'model X';show data","data 'model X'|5|TS10|C 0.008 -0.221 0|O 1.1831 -0.221 0|O -1.0411 0.5156 0|H -1.5835 -0.698202 0|H -0.764698 -1.3442 0|end 'model X';show data"],"name":["TS8","TS9","TS10"],"start":["MIN3","MIN3","MIN1"],"vibr_displace":[[[["0.033400","0.050400","0.000000"],["-0.025100","-0.035800","0.000000"],["-0.023300","-0.037000","0.000000"],["-0.035800","0.424500","0.000000"],["0.405200","0.131600","0.000000"]],[["0.000000","0.000000","0.112000"],["0.000000","0.000000","-0.065900"],["0.000000","0.000000","-0.065900"],["0.000000","0.000000","0.378100"],["0.000000","0.000000","0.378100"]],[["-0.025600","0.017000","0.000000"],["0.008700","0.042500","0.000000"],["-0.035800","-0.024500","0.000000"],["0.367800","-0.243800","0.000000"],["0.367400","-0.244500","0.000000"]],[["-0.093900","-0.141400","0.000000"],["-0.126600","0.146500","0.000000"],["0.184200","-0.059900","0.000000"],["0.001700","0.221500","0.000000"],["0.203500","0.087500","0.000000"]],[["0.000000","0.000000","0.000000"],["0.000000","0.000000","0.014300"],["0.000000","0.000000","-0.014300"],["0.000000","0.000000","-0.485700"],["0.000000","0.000000","0.485700"]],[["0.111700","0.168200","0.000000"],["-0.122700","0.008100","0.000000"],["0.055000","-0.109900","0.000000"],["-0.005000","-0.276100","0.000000"],["-0.252500","-0.111600","0.000000"]],[["0.234700","-0.155900","0.000000"],["-0.131500","0.009500","0.000000"],["-0.059700","0.117600","0.000000"],["0.216900","0.067700","0.000000"],["0.021700","-0.226200","0.000000"]],[["-0.030700","0.020400","0.000000"],["0.009000","-0.011600","0.000000"],["0.014100","-0.003800","0.000000"],["0.257700","0.389400","0.000000"],["-0.259000","-0.388500","0.000000"]],[["-0.003500","-0.005200","0.000000"],["-0.008900","-0.012700","0.000000"],["-0.008300","-0.013100","0.000000"],["0.480900","0.020500","0.000000"],["-0.167500","0.451300","0.000000"]]],[[["-0.070900","0.030200","0.000000"],["-0.001600","-0.040900","0.000000"],["0.008000","0.011300","0.000000"],["0.799402","0.076100","0.000000"],["-0.056200","0.033000","0.000000"]],[["0.000000","0.000000","-0.131700"],["0.000000","0.000000","0.002600"],["0.000000","0.000000","0.044600"],["0.000000","0.000000","0.759301"],["0.000000","0.000000","0.061900"]],[["-0.081900","-0.140500","0.000000"],["-0.081100","0.124900","0.000000"],["0.108200","-0.007300","0.000000"],["0.100200","-0.017300","0.000000"],["0.445100","-0.175100","0.000000"]],[["0.000000","0.000000","0.066300"],["0.000000","0.000000","-0.034800"],["0.000000","0.000000","-0.067400"],["0.000000","0.000000","0.145900"],["0.000000","0.000000","0.685703"]],[["0.013100","-0.028000","0.000000"],["0.073900","-0.022300","0.000000"],["-0.127500","0.060700","0.000000"],["0.072800","-0.007000","0.000000"],["0.622302","-0.267800","0.000000"]],[["-0.034300","0.194300","0.000000"],["-0.027500","-0.009400","0.000000"],["0.033700","-0.133600","0.000000"],["-0.086800","0.165100","0.000000"],["0.396500","-0.211400","0.000000"]],[["0.248100","0.024700","0.000000"],["-0.187700","-0.009800","0.000000"],["-0.024300","0.012700","0.000000"],["0.276000","-0.265000","0.000000"],["0.132800","-0.075100","0.000000"]],[["0.050400","-0.050500","0.000000"],["-0.031500","-0.024000","0.000000"],["-0.003700","0.009100","0.000000"],["-0.027600","0.855600","0.000000"],["-0.014900","-0.017800","0.000000"]],[["0.004000","-0.015900","0.000000"],["-0.003200","-0.000800","0.000000"],["-0.022500","-0.041400","0.000000"],["0.011900","0.013300","0.000000"],["0.349700","0.845800","0.000000"]]],[[["0.008500","-0.056500","0.000000"],["0.006600","0.008600","0.000000"],["-0.028300","-0.007000","0.000000"],["-0.094400","0.540899","0.000000"],["0.337700","0.105600","0.000000"]],[["0.000000","0.000000","0.207300"],["0.000000","0.000000","-0.083000"],["0.000000","0.000000","-0.029900"],["0.000000","0.000000","-0.641601"],["0.000000","0.000000","-0.038200"]],[["0.090300","0.258300","0.000000"],["0.095500","-0.161800","0.000000"],["-0.183600","-0.040100","0.000000"],["0.082800","0.008000","0.000000"],["0.239400","0.119200","0.000000"]],[["-0.032700","-0.018100","0.000000"],["-0.027100","-0.015300","0.000000"],["0.000800","0.049400","0.000000"],["0.331000","-0.213700","0.000000"],["0.475300","-0.111200","0.000000"]],[["0.000000","0.000000","-0.075400"],["0.000000","0.000000","0.017000"],["0.000000","0.000000","0.032200"],["0.000000","0.000000","-0.378900"],["0.000000","0.000000","0.496400"]],[["0.009300","-0.173300","0.000000"],["0.167300","0.014200","0.000000"],["-0.159100","0.134600","0.000000"],["-0.004100","-0.189100","0.000000"],["-0.237400","-0.108600","0.000000"]],[["-0.100200","0.006300","0.000000"],["0.058300","-0.002100","0.000000"],["0.025200","-0.010600","0.000000"],["-0.340900","-0.223100","0.000000"],["0.208100","0.349100","0.000000"]],[["0.153800","-0.046900","0.000000"],["-0.083400","0.000800","0.000000"],["-0.025600","0.026900","0.000000"],["-0.202800","-0.258100","0.000000"],["0.099600","0.377500","0.000000"]],[["-0.026100","-0.019100","0.000000"],["0.006900","0.001600","0.000000"],["-0.005000","-0.020200","0.000000"],["0.444600","0.060300","0.000000"],["-0.163700","0.463100","0.000000"]]]],"vibr_models":["load data 'model FRQS'\n5\nfrq = -1949.30 cm-1 \nC -0.3115 -0.469 0 0 0.033400 0.050400 0.000000\nO 0.935099 -0.469 0 0 -0.025100 -0.035800 0.000000\nO -0.794898 0.680099 0 0 -0.023300 -0.037000 0.000000\nH 1.1494 0.850298 0 0 -0.035800 0.424500 0.000000\nH 0.3382 1.3891 0 0 0.405200 0.131600 0.000000\n5\nfrq = 777.70 cm-1 \nC -0.3115 -0.469 0 0 0.000000 0.000000 0.112000\nO 0.935099 -0.469 0 0 0.000000 0.000000 -0.065900\nO -0.794898 0.680099 0 0 0.000000 0.000000 -0.065900\nH 1.1494 0.850298 0 0 0.000000 0.000000 0.378100\nH 0.3382 1.3891 0 0 0.000000 0.000000 0.378100\n5\nfrq = 861.40 cm-1 \nC -0.3115 -0.469 0 0 -0.025600 0.017000 0.000000\nO 0.935099 -0.469 0 0 0.008700 0.042500 0.000000\nO -0.794898 0.680099 0 0 -0.035800 -0.024500 0.000000\nH 1.1494 0.850298 0 0 0.367800 -0.243800 0.000000\nH 0.3382 1.3891 0 0 0.367400 -0.244500 0.000000\n5\nfrq = 895.40 cm-1 \nC -0.3115 -0.469 0 0 -0.093900 -0.141400 0.000000\nO 0.935099 -0.469 0 0 -0.126600 0.146500 0.000000\nO -0.794898 0.680099 0 0 0.184200 -0.059900 0.000000\nH 1.1494 0.850298 0 0 0.001700 0.221500 0.000000\nH 0.3382 1.3891 0 0 0.203500 0.087500 0.000000\n5\nfrq = 1136.40 cm-1 \nC -0.3115 -0.469 0 0 0.000000 0.000000 0.000000\nO 0.935099 -0.469 0 0 0.000000 0.000000 0.014300\nO -0.794898 0.680099 0 0 0.000000 0.000000 -0.014300\nH 1.1494 0.850298 0 0 0.000000 0.000000 -0.485700\nH 0.3382 1.3891 0 0 0.000000 0.000000 0.485700\n5\nfrq = 1577.90 cm-1 \nC -0.3115 -0.469 0 0 0.111700 0.168200 0.000000\nO 0.935099 -0.469 0 0 -0.122700 0.008100 0.000000\nO -0.794898 0.680099 0 0 0.055000 -0.109900 0.000000\nH 1.1494 0.850298 0 0 -0.005000 -0.276100 0.000000\nH 0.3382 1.3891 0 0 -0.252500 -0.111600 0.000000\n5\nfrq = 1608.90 cm-1 \nC -0.3115 -0.469 0 0 0.234700 -0.155900 0.000000\nO 0.935099 -0.469 0 0 -0.131500 0.009500 0.000000\nO -0.794898 0.680099 0 0 -0.059700 0.117600 0.000000\nH 1.1494 0.850298 0 0 0.216900 0.067700 0.000000\nH 0.3382 1.3891 0 0 0.021700 -0.226200 0.000000\n5\nfrq = 2040.00 cm-1 \nC -0.3115 -0.469 0 0 -0.030700 0.020400 0.000000\nO 0.935099 -0.469 0 0 0.009000 -0.011600 0.000000\nO -0.794898 0.680099 0 0 0.014100 -0.003800 0.000000\nH 1.1494 0.850298 0 0 0.257700 0.389400 0.000000\nH 0.3382 1.3891 0 0 -0.259000 -0.388500 0.000000\n5\nfrq = 2711.80 cm-1 \nC -0.3115 -0.469 0 0 -0.003500 -0.005200 0.000000\nO 0.935099 -0.469 0 0 -0.008900 -0.012700 0.000000\nO -0.794898 0.680099 0 0 -0.008300 -0.013100 0.000000\nH 1.1494 0.850298 0 0 0.480900 0.020500 0.000000\nH 0.3382 1.3891 0 0 -0.167500 0.451300 0.000000\nend 'model FRQS'","load data 'model FRQS'\n5\nfrq = -1969.90 cm-1 \nC -0.2036 -0.3962 0 0 -0.070900 0.030200 0.000000\nO 1.0655 -0.3962 0 0 -0.001600 -0.040900 0.000000\nO -0.911899 0.694699 0 0 0.008000 0.011300 0.000000\nH 0.4048 -1.5715 0 0 0.799402 0.076100 0.000000\nH -0.4169 1.5559 0.0001 0 -0.056200 0.033000 0.000000\n5\nfrq = 354.70 cm-1 \nC -0.2036 -0.3962 0 0 0.000000 0.000000 -0.131700\nO 1.0655 -0.3962 0 0 0.000000 0.000000 0.002600\nO -0.911899 0.694699 0 0 0.000000 0.000000 0.044600\nH 0.4048 -1.5715 0 0 0.000000 0.000000 0.759301\nH -0.4169 1.5559 0.0001 0 0.000000 0.000000 0.061900\n5\nfrq = 496.20 cm-1 \nC -0.2036 -0.3962 0 0 -0.081900 -0.140500 0.000000\nO 1.0655 -0.3962 0 0 -0.081100 0.124900 0.000000\nO -0.911899 0.694699 0 0 0.108200 -0.007300 0.000000\nH 0.4048 -1.5715 0 0 0.100200 -0.017300 0.000000\nH -0.4169 1.5559 0.0001 0 0.445100 -0.175100 0.000000\n5\nfrq = 641.80 cm-1 \nC -0.2036 -0.3962 0 0 0.000000 0.000000 0.066300\nO 1.0655 -0.3962 0 0 0.000000 0.000000 -0.034800\nO -0.911899 0.694699 0 0 0.000000 0.000000 -0.067400\nH 0.4048 -1.5715 0 0 0.000000 0.000000 0.145900\nH -0.4169 1.5559 0.0001 0 0.000000 0.000000 0.685703\n5\nfrq = 1064.40 cm-1 \nC -0.2036 -0.3962 0 0 0.013100 -0.028000 0.000000\nO 1.0655 -0.3962 0 0 0.073900 -0.022300 0.000000\nO -0.911899 0.694699 0 0 -0.127500 0.060700 0.000000\nH 0.4048 -1.5715 0 0 0.072800 -0.007000 0.000000\nH -0.4169 1.5559 0.0001 0 0.622302 -0.267800 0.000000\n5\nfrq = 1342.00 cm-1 \nC -0.2036 -0.3962 0 0 -0.034300 0.194300 0.000000\nO 1.0655 -0.3962 0 0 -0.027500 -0.009400 0.000000\nO -0.911899 0.694699 0 0 0.033700 -0.133600 0.000000\nH 0.4048 -1.5715 0 0 -0.086800 0.165100 0.000000\nH -0.4169 1.5559 0.0001 0 0.396500 -0.211400 0.000000\n5\nfrq = 1552.50 cm-1 \nC -0.2036 -0.3962 0 0 0.248100 0.024700 0.000000\nO 1.0655 -0.3962 0 0 -0.187700 -0.009800 0.000000\nO -0.911899 0.694699 0 0 -0.024300 0.012700 0.000000\nH 0.4048 -1.5715 0 0 0.276000 -0.265000 0.000000\nH -0.4169 1.5559 0.0001 0 0.132800 -0.075100 0.000000\n5\nfrq = 2462.80 cm-1 \nC -0.2036 -0.3962 0 0 0.050400 -0.050500 0.000000\nO 1.0655 -0.3962 0 0 -0.031500 -0.024000 0.000000\nO -0.911899 0.694699 0 0 -0.003700 0.009100 0.000000\nH 0.4048 -1.5715 0 0 -0.027600 0.855600 0.000000\nH -0.4169 1.5559 0.0001 0 -0.014900 -0.017800 0.000000\n5\nfrq = 2757.40 cm-1 \nC -0.2036 -0.3962 0 0 0.004000 -0.015900 0.000000\nO 1.0655 -0.3962 0 0 -0.003200 -0.000800 0.000000\nO -0.911899 0.694699 0 0 -0.022500 -0.041400 0.000000\nH 0.4048 -1.5715 0 0 0.011900 0.013300 0.000000\nH -0.4169 1.5559 0.0001 0 0.349700 0.845800 0.000000\nend 'model FRQS'","load data 'model FRQS'\n5\nfrq = -2928.30 cm-1 \nC 0.008 -0.221 0 0 0.008500 -0.056500 0.000000\nO 1.1831 -0.221 0 0 0.006600 0.008600 0.000000\nO -1.0411 0.5156 0 0 -0.028300 -0.007000 0.000000\nH -1.5835 -0.698202 0 0 -0.094400 0.540899 0.000000\nH -0.764698 -1.3442 0 0 0.337700 0.105600 0.000000\n5\nfrq = 451.40 cm-1 \nC 0.008 -0.221 0 0 0.000000 0.000000 0.207300\nO 1.1831 -0.221 0 0 0.000000 0.000000 -0.083000\nO -1.0411 0.5156 0 0 0.000000 0.000000 -0.029900\nH -1.5835 -0.698202 0 0 0.000000 0.000000 -0.641601\nH -0.764698 -1.3442 0 0 0.000000 0.000000 -0.038200\n5\nfrq = 587.20 cm-1 \nC 0.008 -0.221 0 0 0.090300 0.258300 0.000000\nO 1.1831 -0.221 0 0 0.095500 -0.161800 0.000000\nO -1.0411 0.5156 0 0 -0.183600 -0.040100 0.000000\nH -1.5835 -0.698202 0 0 0.082800 0.008000 0.000000\nH -0.764698 -1.3442 0 0 0.239400 0.119200 0.000000\n5\nfrq = 904.50 cm-1 \nC 0.008 -0.221 0 0 -0.032700 -0.018100 0.000000\nO 1.1831 -0.221 0 0 -0.027100 -0.015300 0.000000\nO -1.0411 0.5156 0 0 0.000800 0.049400 0.000000\nH -1.5835 -0.698202 0 0 0.331000 -0.213700 0.000000\nH -0.764698 -1.3442 0 0 0.475300 -0.111200 0.000000\n5\nfrq = 997.80 cm-1 \nC 0.008 -0.221 0 0 0.000000 0.000000 -0.075400\nO 1.1831 -0.221 0 0 0.000000 0.000000 0.017000\nO -1.0411 0.5156 0 0 0.000000 0.000000 0.032200\nH -1.5835 -0.698202 0 0 0.000000 0.000000 -0.378900\nH -0.764698 -1.3442 0 0 0.000000 0.000000 0.496400\n5\nfrq = 1345.80 cm-1 \nC 0.008 -0.221 0 0 0.009300 -0.173300 0.000000\nO 1.1831 -0.221 0 0 0.167300 0.014200 0.000000\nO -1.0411 0.5156 0 0 -0.159100 0.134600 0.000000\nH -1.5835 -0.698202 0 0 -0.004100 -0.189100 0.000000\nH -0.764698 -1.3442 0 0 -0.237400 -0.108600 0.000000\n5\nfrq = 1946.90 cm-1 \nC 0.008 -0.221 0 0 -0.100200 0.006300 0.000000\nO 1.1831 -0.221 0 0 0.058300 -0.002100 0.000000\nO -1.0411 0.5156 0 0 0.025200 -0.010600 0.000000\nH -1.5835 -0.698202 0 0 -0.340900 -0.223100 0.000000\nH -0.764698 -1.3442 0 0 0.208100 0.349100 0.000000\n5\nfrq = 2103.10 cm-1 \nC 0.008 -0.221 0 0 0.153800 -0.046900 0.000000\nO 1.1831 -0.221 0 0 -0.083400 0.000800 0.000000\nO -1.0411 0.5156 0 0 -0.025600 0.026900 0.000000\nH -1.5835 -0.698202 0 0 -0.202800 -0.258100 0.000000\nH -0.764698 -1.3442 0 0 0.099600 0.377500 0.000000\n5\nfrq = 2635.00 cm-1 \nC 0.008 -0.221 0 0 -0.026100 -0.019100 0.000000\nO 1.1831 -0.221 0 0 0.006900 0.001600 0.000000\nO -1.0411 0.5156 0 0 -0.005000 -0.020200 0.000000\nH -1.5835 -0.698202 0 0 0.444600 0.060300 0.000000\nH -0.764698 -1.3442 0 0 -0.163700 0.463100 0.000000\nend 'model FRQS'"],"zpe":[16.596,15.256,15.685]},"selected":{"id":"1152"},"selection_policy":{"id":"1280"}},"id":"1039","type":"ColumnDataSource"},{"attributes":{},"id":"1038","type":"MultiLine"},{"attributes":{"bottom_units":"screen","fill_alpha":0.5,"fill_color":"lightgrey","left_units":"screen","level":"overlay","line_alpha":1.0,"line_color":"black","line_dash":[4,4],"line_width":2,"right_units":"screen","syncable":false,"top_units":"screen"},"id":"1026","type":"BoxAnnotation"},{"attributes":{"source":{"id":"1191"}},"id":"1211","type":"CDSView"},{"attributes":{"args":{"jsmol":{"id":"1126"},"layout":{"id":"1159"},"mol_view":[false],"prof":{"id":"1161"}},"code":"\n\t\t//layout - full layout, jsmol - jsmol app (where the profile window will appear), prof - profile\n\t\t//mol_view - boolean for current view\n\t\tvar current_fig = layout.children[1][0].children\n\t\tvar name_elem = current_fig[0]\n\t\tvar jsmol_plot = current_fig[1]\n\t\tvar prof_plot = current_fig[2]\n\t\tvar control_elem = current_fig[3]\n\n\t\t// allow to enable and disable corresponding profile control buttons\n\t\tvar bfilt1 = control_elem.children[1] \n\t\tvar bfilt2 = control_elem.children[3]\n\t\tif (mol_view[0] == true){\n\t\t\tjsmol.visible = true\n\t\t\tprof.visible = false\n\t\t\tbfilt1.disabled = true\n\t\t\tbfilt2.disabled = true\n\t\t\tmol_view[0] = false\n\t\t} else {\n\t\t\tjsmol.visible = false\n\t\t\tprof.visible = true\n\t\t\tbfilt1.disabled = false\n\t\t\tbfilt2.disabled = false\n\t\t\tmol_view[0] = true\n\t\t}\n\n\t\t"},"id":"1226","type":"CustomJS"},{"attributes":{"data_source":{"id":"1191"},"glyph":{"id":"1208"},"hover_glyph":null,"muted_glyph":null,"nonselection_glyph":{"id":"1209"},"view":{"id":"1211"}},"id":"1210","type":"GlyphRenderer"},{"attributes":{},"id":"1267","type":"Selection"},{"attributes":{"args":{"hover":{"id":"1222"}},"code":"\n\t// all renderers are checked at once: when the triggered one is caught (by changes in cb_data.index.line_indices)\n\t// access its data_source directly\n\tif (cb_data.index.line_indices.length &gt; 0) {\n\t\tvar ndx = cb_data.index.line_indices[0]\n\t\tvar rend_obj = cb_data.renderer.data_source\n\t\tvar formula_list = rend_obj.data[\"form\"]\n\n\t// fallback for cases without any fragmentation, where formula is not defined\n\t    if (formula_list){\n\t\t\tvar formula = formula_list[ndx]\n\t\t} else {\n\t\t\tvar formula = null\n\t\t}\n\t    if (!formula){\n\t\t\thover.tooltips = [[\"tag\",\"@lab\"],[\"E\",\"@y{%.2f}\"]]\n \t\t} else {\n\t\t\thover.tooltips = [[\"tag\",\"@lab\"],[\"E\",\"@y{%.2f}\"],[\"formula\",\"@form\"]]\n\t\t}\n\t}\n\t"},"id":"1223","type":"CustomJS"},{"attributes":{},"id":"1268","type":"UnionRenderers"},{"attributes":{"args":{"prof":{"id":"1161"}},"code":"\n\t\t// prof - bokeh figure for the profiles\n\t\t// make everything visible\n\t\tvar Nlines = prof.renderers.length/2\n\t\tfor (let i = 0 ; i &lt; Nlines ; i++){\n\t\t\tconst index = 2*i\n\t\t\tprof.renderers[index].visible = true\n\t\t\tprof.renderers[index+1].visible = true\n\t\t\t// tags are EVEN entries and energies are ODD\n\t\t\tprof.center[index+2].visible = true\n\t\t\tprof.center[index+3].visible = false\n\t\t\tconsole.log(prof.center[index+3].text,prof.center[index+3].visible)\n\t\t}\n\t\t"},"id":"1225","type":"CustomJS"},{"attributes":{"data":{"elab":["-8.3","58.2","-37.6"],"lab":["MIN1","TS10","PR2"],"x":{"__ndarray__":"AAAAAAAA4D8AAAAAAAAEQAAAAAAAABJA","dtype":"float64","order":"little","shape":[3]},"y":[-8.323520999999985,58.154156,-37.60977]},"selected":{"id":"1269"},"selection_policy":{"id":"1270"}},"id":"1212","type":"ColumnDataSource"},{"attributes":{"callback":{"id":"1223"},"description":"Profile hover","formatters":{"@y":"printf"},"line_policy":"interp","renderers":[{"id":"1195"},{"id":"1210"}],"tooltips":[["tag","@lab"],["E","@y{%.2f}"],["formula","@form"]]},"id":"1222","type":"HoverTool"},{"attributes":{},"id":"1014","type":"BasicTicker"},{"attributes":{"data":{"elab":["-8.3","53.2","5.3","44.0","-37.6"],"lab":["MIN1","TS9","MIN3","TS8","PR2"],"x":{"__ndarray__":"AAAAAAAA4D8AAAAAAAAEQAAAAAAAABJAAAAAAAAAGkAAAAAAAAAhQA==","dtype":"float64","order":"little","shape":[5]},"y":[-8.323520999999985,53.170985,5.312669999999997,43.967468999999994,-37.60977]},"selected":{"id":"1273"},"selection_policy":{"id":"1274"}},"id":"1204","type":"ColumnDataSource"},{"attributes":{"align":"center","css_classes":["xtest"],"icon":null,"js_event_callbacks":{"button_click":[{"id":"1140"}]},"label":"To clipboard","max_width":200},"id":"1128","type":"Button"},{"attributes":{"line_alpha":0.1,"line_color":"#1f77b4","x":{"field":"x"},"y":{"field":"y"}},"id":"1194","type":"Line"},{"attributes":{"active_multi":null,"tools":[{"id":"1021"},{"id":"1022"},{"id":"1023"},{"id":"1024"},{"id":"1025"},{"id":"1093"},{"id":"1099"},{"id":"1112"},{"id":"1115"}]},"id":"1027","type":"Toolbar"},{"attributes":{"name":"ENERGYLABELS","source":{"id":"1204"},"text":{"field":"elab"},"text_font_size":{"value":"1.4vh"},"visible":false,"x":{"field":"x"},"x_offset":{"value":-1},"y":{"field":"y"},"y_offset":{"value":-20}},"id":"1205","type":"LabelSet"},{"attributes":{"data_source":{"id":"1197"},"glyph":{"id":"1198"},"hover_glyph":null,"muted_glyph":null,"view":{"id":"1200"}},"id":"1199","type":"GlyphRenderer"},{"attributes":{},"id":"1166","type":"LinearScale"},{"attributes":{},"id":"1235","type":"AllLabels"},{"attributes":{},"id":"1236","type":"BasicTickFormatter"},{"attributes":{"fill_alpha":{"value":0.9},"fill_color":{"value":"#1f77b4"},"height":{"value":2},"line_width":{"value":0},"name":"RECTANGLE","width":{"value":1},"x":{"field":"x"},"y":{"field":"y"}},"id":"1198","type":"Rect"},{"attributes":{"source":{"id":"1197"}},"id":"1200","type":"CDSView"},{"attributes":{"fill_alpha":{"value":0.8},"fill_color":{"field":"energy","transform":{"id":"1051"}},"line_width":{"value":2},"size":{"value":30}},"id":"1057","type":"Circle"},{"attributes":{"data":{"elab":["-8.3","53.2","5.3","44.0","-37.6"],"lab":["MIN1","TS9","MIN3","TS8","PR2"],"x":{"__ndarray__":"AAAAAAAA4D8AAAAAAAAEQAAAAAAAABJAAAAAAAAAGkAAAAAAAAAhQA==","dtype":"float64","order":"little","shape":[5]},"y":[-8.323520999999985,53.170985,5.312669999999997,43.967468999999994,-37.60977]},"selected":{"id":"1265"},"selection_policy":{"id":"1266"}},"id":"1197","type":"ColumnDataSource"},{"attributes":{"fill_alpha":{"value":0.5},"fill_color":{"field":"energy","transform":{"id":"1051"}},"size":{"value":30}},"id":"1052","type":"Circle"},{"attributes":{"data":{"elab":["-8.3","53.2","5.3","44.0","-37.6"],"lab":["MIN1","TS9","MIN3","TS8","PR2"],"x":{"__ndarray__":"AAAAAAAA4D8AAAAAAAAEQAAAAAAAABJAAAAAAAAAGkAAAAAAAAAhQA==","dtype":"float64","order":"little","shape":[5]},"y":[-8.323520999999985,53.170985,5.312669999999997,43.967468999999994,-37.60977]},"selected":{"id":"1271"},"selection_policy":{"id":"1272"}},"id":"1201","type":"ColumnDataSource"},{"attributes":{"source":{"id":"1190"}},"id":"1196","type":"CDSView"},{"attributes":{"data":{"energy":[5.312669999999997,-8.323520999999985,-37.60977],"fname":["MIN3_minf_ts5_batch5.rxyz","MIN1_minf_ts2_torsion.rxyz","PR2_minr_ts3_batch4"],"formula":[null,null,"H2 + CO2"],"frequencies":["468.3;527.1;594.2;820.9;850.1;1243.7;1341.7;2737.5;2762.5","551.3;593.0;956.9;1062.5;1231.9;1308.0;1840.5;2664.1;2787.8","-29.4;150.5;425.8;429.4;517.9;521.2;1375.1;2344.5;5508.7"],"geometry":["C -0.2802 -0.463 0\nO 1.0252 -0.463 0\nO -0.885499 0.693698 0\nH 1.6316 0.2794 0.0005\nH -0.5108 1.5757 0.0003","C -0.1526 -0.426 -0.0001\nO 1.0475 -0.426 -0.0001\nO -0.859601 0.727402 -0.0001\nH -0.8517 -1.2589 -0.000700001\nH -0.3126 1.5507 0.003","C 0.079120000 0.098620000 0.000000000\nO 1.174170000 -0.307940000 0.000000000\nO -1.014880000 0.508180000 -0.000010000\nH -2.088380000 -2.029700000 0.000090000\nH -1.383530000 -2.324960000 0.000080000"],"index":["MIN3","MIN1","PR2"],"model":["data 'model X'|5|MIN3|C -0.2802 -0.463 0|O 1.0252 -0.463 0|O -0.885499 0.693698 0|H 1.6316 0.2794 0.0005|H -0.5108 1.5757 0.0003|end 'model X';show data","data 'model X'|5|MIN1|C -0.1526 -0.426 -0.0001|O 1.0475 -0.426 -0.0001|O -0.859601 0.727402 -0.0001|H -0.8517 -1.2589 -0.000700001|H -0.3126 1.5507 0.003|end 'model X';show data","data 'model X'|5|PR2|C 0.079120000 0.098620000 0.000000000|O 1.174170000 -0.307940000 0.000000000|O -1.014880000 0.508180000 -0.000010000|H -2.088380000 -2.029700000 0.000090000|H -1.383530000 -2.324960000 0.000080000|end 'model X';show data"],"name":["MIN3","MIN1","PR2"],"vibr_displace":[[[["-0.059600","-0.098500","0.000000"],["-0.045800","0.071800","0.000000"],["0.084900","-0.007300","0.000000"],["-0.251300","0.253900","0.000000"],["0.341400","-0.105100","0.000000"]],[["0.000000","0.000000","0.000000"],["0.000000","0.000000","0.033000"],["0.000000","0.000000","-0.033000"],["0.000000","0.000100","-0.467000"],["0.000000","-0.000100","0.467000"]],[["0.000000","0.000000","0.092800"],["0.000000","0.000000","-0.059600"],["0.000000","0.000000","-0.059700"],["0.000000","-0.000100","0.394000"],["0.000000","-0.000100","0.393900"]],[["0.025100","-0.015200","0.000000"],["-0.020800","0.039400","0.000000"],["-0.044500","0.000200","0.000000"],["0.320900","-0.302200","0.000000"],["0.416500","-0.144300","0.000000"]],[["0.006900","0.011400","0.000000"],["0.080300","-0.063100","0.000000"],["-0.093200","0.041900","0.000000"],["-0.258800","0.293300","0.000000"],["0.379800","-0.093400","0.000000"]],[["0.112900","0.186400","0.000000"],["-0.116200","-0.027000","0.000000"],["0.029900","-0.115500","0.000000"],["-0.219100","0.160600","0.000100"],["0.243900","-0.119700","0.000100"]],[["0.270100","-0.163500","0.000000"],["-0.132800","-0.015900","0.000000"],["-0.047500","0.125100","0.000000"],["-0.181800","0.101800","0.000000"],["-0.174500","0.113900","0.000000"]],[["0.005800","-0.003500","0.000000"],["0.016700","0.018600","0.000000"],["-0.008700","-0.023500","0.000000"],["-0.335200","-0.331700","-0.000100"],["0.138600","0.450700","0.000100"]],[["-0.007000","-0.011500","0.000000"],["-0.017800","-0.016600","0.000000"],["-0.006500","-0.023500","0.000000"],["0.339400","0.323600","0.000100"],["0.129400","0.450700","0.000100"]]],[[["-0.079500","-0.129800","-0.000100"],["-0.076600","0.134200","0.000100"],["0.113200","-0.016400","0.000000"],["-0.026200","-0.146600","0.000100"],["0.391900","-0.176900","-0.000100"]],[["0.000000","0.000000","0.071100"],["0.000000","0.000000","-0.032800"],["0.000000","0.000000","-0.072500"],["0.000000","-0.000100","0.107100"],["-0.000200","-0.000500","0.716601"]],[["0.000000","0.000100","-0.140800"],["0.000000","0.000000","0.036800"],["0.000000","0.000000","0.018200"],["-0.000200","-0.000400","0.684300"],["0.000000","-0.000100","0.120000"]],[["0.005700","-0.019500","0.000000"],["0.029900","-0.021800","0.000000"],["-0.078800","0.064200","0.000000"],["0.313800","-0.258400","-0.000100"],["0.394400","-0.182400","0.000000"]],[["-0.041000","0.135200","0.000000"],["-0.062400","-0.017700","0.000000"],["0.065800","-0.065900","0.000000"],["0.511100","-0.332500","-0.000200"],["-0.076800","0.048700","0.000100"]],[["-0.011500","0.105900","0.000000"],["0.010700","-0.001400","0.000000"],["-0.014800","-0.077400","-0.000100"],["-0.207600","0.244600","0.000000"],["0.409900","-0.255800","0.000000"]],[["0.319100","-0.016300","0.000000"],["-0.235000","-0.001400","0.000000"],["-0.020800","0.006400","0.000000"],["0.129800","0.221800","0.000000"],["0.127600","-0.107400","0.000100"]],[["-0.052500","-0.046300","0.000000"],["0.006400","-0.001800","0.000000"],["0.002500","-0.000100","0.000000"],["0.526900","0.652999","0.000600"],["-0.043000","-0.069500","-0.000100"]],[["-0.000300","-0.016800","-0.000100"],["-0.002300","0.000000","0.000100"],["-0.024200","-0.039800","0.000000"],["0.040500","0.061800","-0.000400"],["0.383500","0.770202","0.000900"]]],[]],"vibr_models":["load data 'model FRQS'\n5\nfrq = 468.30 cm-1 \nC -0.2802 -0.463 0 0 -0.059600 -0.098500 0.000000\nO 1.0252 -0.463 0 0 -0.045800 0.071800 0.000000\nO -0.885499 0.693698 0 0 0.084900 -0.007300 0.000000\nH 1.6316 0.2794 0.0005 0 -0.251300 0.253900 0.000000\nH -0.5108 1.5757 0.0003 0 0.341400 -0.105100 0.000000\n5\nfrq = 527.10 cm-1 \nC -0.2802 -0.463 0 0 0.000000 0.000000 0.000000\nO 1.0252 -0.463 0 0 0.000000 0.000000 0.033000\nO -0.885499 0.693698 0 0 0.000000 0.000000 -0.033000\nH 1.6316 0.2794 0.0005 0 0.000000 0.000100 -0.467000\nH -0.5108 1.5757 0.0003 0 0.000000 -0.000100 0.467000\n5\nfrq = 594.20 cm-1 \nC -0.2802 -0.463 0 0 0.000000 0.000000 0.092800\nO 1.0252 -0.463 0 0 0.000000 0.000000 -0.059600\nO -0.885499 0.693698 0 0 0.000000 0.000000 -0.059700\nH 1.6316 0.2794 0.0005 0 0.000000 -0.000100 0.394000\nH -0.5108 1.5757 0.0003 0 0.000000 -0.000100 0.393900\n5\nfrq = 820.90 cm-1 \nC -0.2802 -0.463 0 0 0.025100 -0.015200 0.000000\nO 1.0252 -0.463 0 0 -0.020800 0.039400 0.000000\nO -0.885499 0.693698 0 0 -0.044500 0.000200 0.000000\nH 1.6316 0.2794 0.0005 0 0.320900 -0.302200 0.000000\nH -0.5108 1.5757 0.0003 0 0.416500 -0.144300 0.000000\n5\nfrq = 850.10 cm-1 \nC -0.2802 -0.463 0 0 0.006900 0.011400 0.000000\nO 1.0252 -0.463 0 0 0.080300 -0.063100 0.000000\nO -0.885499 0.693698 0 0 -0.093200 0.041900 0.000000\nH 1.6316 0.2794 0.0005 0 -0.258800 0.293300 0.000000\nH -0.5108 1.5757 0.0003 0 0.379800 -0.093400 0.000000\n5\nfrq = 1243.70 cm-1 \nC -0.2802 -0.463 0 0 0.112900 0.186400 0.000000\nO 1.0252 -0.463 0 0 -0.116200 -0.027000 0.000000\nO -0.885499 0.693698 0 0 0.029900 -0.115500 0.000000\nH 1.6316 0.2794 0.0005 0 -0.219100 0.160600 0.000100\nH -0.5108 1.5757 0.0003 0 0.243900 -0.119700 0.000100\n5\nfrq = 1341.70 cm-1 \nC -0.2802 -0.463 0 0 0.270100 -0.163500 0.000000\nO 1.0252 -0.463 0 0 -0.132800 -0.015900 0.000000\nO -0.885499 0.693698 0 0 -0.047500 0.125100 0.000000\nH 1.6316 0.2794 0.0005 0 -0.181800 0.101800 0.000000\nH -0.5108 1.5757 0.0003 0 -0.174500 0.113900 0.000000\n5\nfrq = 2737.50 cm-1 \nC -0.2802 -0.463 0 0 0.005800 -0.003500 0.000000\nO 1.0252 -0.463 0 0 0.016700 0.018600 0.000000\nO -0.885499 0.693698 0 0 -0.008700 -0.023500 0.000000\nH 1.6316 0.2794 0.0005 0 -0.335200 -0.331700 -0.000100\nH -0.5108 1.5757 0.0003 0 0.138600 0.450700 0.000100\n5\nfrq = 2762.50 cm-1 \nC -0.2802 -0.463 0 0 -0.007000 -0.011500 0.000000\nO 1.0252 -0.463 0 0 -0.017800 -0.016600 0.000000\nO -0.885499 0.693698 0 0 -0.006500 -0.023500 0.000000\nH 1.6316 0.2794 0.0005 0 0.339400 0.323600 0.000100\nH -0.5108 1.5757 0.0003 0 0.129400 0.450700 0.000100\nend 'model FRQS'","load data 'model FRQS'\n5\nfrq = 551.30 cm-1 \nC -0.1526 -0.426 -0.0001 0 -0.079500 -0.129800 -0.000100\nO 1.0475 -0.426 -0.0001 0 -0.076600 0.134200 0.000100\nO -0.859601 0.727402 -0.0001 0 0.113200 -0.016400 0.000000\nH -0.8517 -1.2589 -0.000700001 0 -0.026200 -0.146600 0.000100\nH -0.3126 1.5507 0.003 0 0.391900 -0.176900 -0.000100\n5\nfrq = 593.00 cm-1 \nC -0.1526 -0.426 -0.0001 0 0.000000 0.000000 0.071100\nO 1.0475 -0.426 -0.0001 0 0.000000 0.000000 -0.032800\nO -0.859601 0.727402 -0.0001 0 0.000000 0.000000 -0.072500\nH -0.8517 -1.2589 -0.000700001 0 0.000000 -0.000100 0.107100\nH -0.3126 1.5507 0.003 0 -0.000200 -0.000500 0.716601\n5\nfrq = 956.90 cm-1 \nC -0.1526 -0.426 -0.0001 0 0.000000 0.000100 -0.140800\nO 1.0475 -0.426 -0.0001 0 0.000000 0.000000 0.036800\nO -0.859601 0.727402 -0.0001 0 0.000000 0.000000 0.018200\nH -0.8517 -1.2589 -0.000700001 0 -0.000200 -0.000400 0.684300\nH -0.3126 1.5507 0.003 0 0.000000 -0.000100 0.120000\n5\nfrq = 1062.50 cm-1 \nC -0.1526 -0.426 -0.0001 0 0.005700 -0.019500 0.000000\nO 1.0475 -0.426 -0.0001 0 0.029900 -0.021800 0.000000\nO -0.859601 0.727402 -0.0001 0 -0.078800 0.064200 0.000000\nH -0.8517 -1.2589 -0.000700001 0 0.313800 -0.258400 -0.000100\nH -0.3126 1.5507 0.003 0 0.394400 -0.182400 0.000000\n5\nfrq = 1231.90 cm-1 \nC -0.1526 -0.426 -0.0001 0 -0.041000 0.135200 0.000000\nO 1.0475 -0.426 -0.0001 0 -0.062400 -0.017700 0.000000\nO -0.859601 0.727402 -0.0001 0 0.065800 -0.065900 0.000000\nH -0.8517 -1.2589 -0.000700001 0 0.511100 -0.332500 -0.000200\nH -0.3126 1.5507 0.003 0 -0.076800 0.048700 0.000100\n5\nfrq = 1308.00 cm-1 \nC -0.1526 -0.426 -0.0001 0 -0.011500 0.105900 0.000000\nO 1.0475 -0.426 -0.0001 0 0.010700 -0.001400 0.000000\nO -0.859601 0.727402 -0.0001 0 -0.014800 -0.077400 -0.000100\nH -0.8517 -1.2589 -0.000700001 0 -0.207600 0.244600 0.000000\nH -0.3126 1.5507 0.003 0 0.409900 -0.255800 0.000000\n5\nfrq = 1840.50 cm-1 \nC -0.1526 -0.426 -0.0001 0 0.319100 -0.016300 0.000000\nO 1.0475 -0.426 -0.0001 0 -0.235000 -0.001400 0.000000\nO -0.859601 0.727402 -0.0001 0 -0.020800 0.006400 0.000000\nH -0.8517 -1.2589 -0.000700001 0 0.129800 0.221800 0.000000\nH -0.3126 1.5507 0.003 0 0.127600 -0.107400 0.000100\n5\nfrq = 2664.10 cm-1 \nC -0.1526 -0.426 -0.0001 0 -0.052500 -0.046300 0.000000\nO 1.0475 -0.426 -0.0001 0 0.006400 -0.001800 0.000000\nO -0.859601 0.727402 -0.0001 0 0.002500 -0.000100 0.000000\nH -0.8517 -1.2589 -0.000700001 0 0.526900 0.652999 0.000600\nH -0.3126 1.5507 0.003 0 -0.043000 -0.069500 -0.000100\n5\nfrq = 2787.80 cm-1 \nC -0.1526 -0.426 -0.0001 0 -0.000300 -0.016800 -0.000100\nO 1.0475 -0.426 -0.0001 0 -0.002300 0.000000 0.000100\nO -0.859601 0.727402 -0.0001 0 -0.024200 -0.039800 0.000000\nH -0.8517 -1.2589 -0.000700001 0 0.040500 0.061800 -0.000400\nH -0.3126 1.5507 0.003 0 0.383500 0.770202 0.000900\nend 'model FRQS'",null],"zpe":[16.22,18.579,16.116]},"selected":{"id":"1147"},"selection_policy":{"id":"1282"}},"id":"1035","type":"ColumnDataSource"},{"attributes":{"fill_alpha":{"value":0.7},"fill_color":{"field":"energy","transform":{"id":"1051"}},"line_color":{"field":"energy","transform":{"id":"1051"}},"line_width":{"value":2},"size":{"value":30}},"id":"1062","type":"Circle"},{"attributes":{"line_width":{"value":2}},"id":"1071","type":"MultiLine"},{"attributes":{"name":"THELABELS","source":{"id":"1201"},"text":{"field":"lab"},"text_font_size":{"value":"1.4vh"},"x":{"field":"x"},"x_offset":{"value":-1},"y":{"field":"y"},"y_offset":{"value":2}},"id":"1202","type":"LabelSet"},{"attributes":{"formatter":{"id":"1236"},"major_label_policy":{"id":"1235"},"ticker":{"id":"1014"},"visible":false},"id":"1013","type":"LinearAxis"},{"attributes":{"line_color":"#1f77b4","x":{"field":"x"},"y":{"field":"y"}},"id":"1193","type":"Line"},{"attributes":{"data_source":{"id":"1190"},"glyph":{"id":"1193"},"hover_glyph":null,"muted_glyph":null,"nonselection_glyph":{"id":"1194"},"view":{"id":"1196"}},"id":"1195","type":"GlyphRenderer"},{"attributes":{"data":{"nnames":["MIN3","MIN1","PR2"],"xn":[0.9999999999999999,-0.50000005039689,-0.49999994960310984],"yn":[4.7706792907989705e-08,0.8660253835545058,-0.8660254312612986]},"selected":{"id":"1256"},"selection_policy":{"id":"1257"}},"id":"1081","type":"ColumnDataSource"},{"attributes":{"line_color":{"value":"orange"},"line_width":{"value":3}},"id":"1076","type":"MultiLine"},{"attributes":{"source":{"id":"1081"},"text":{"field":"nnames"},"text_font_size":{"value":"1.5vh"},"text_font_style":{"value":"bold"},"x":{"field":"xn"},"y":{"field":"yn"},"y_offset":{"value":5}},"id":"1082","type":"LabelSet"},{"attributes":{"data":{"enames":["TS8","TS9","TS10"],"etuples":[["MIN3","PR2"],["MIN3","MIN1"],["MIN1","PR2"]],"xe":[0.250000025198445,0.24999997480155495,-0.4999999999999999],"ye":[-0.4330126917772528,0.43301271563064936,-2.385339642474449e-08]},"selected":{"id":"1258"},"selection_policy":{"id":"1259"}},"id":"1084","type":"ColumnDataSource"},{"attributes":{},"id":"1246","type":"AllLabels"},{"attributes":{},"id":"1087","type":"EdgesAndLinkedNodes"},{"attributes":{},"id":"1247","type":"BasicTickFormatter"},{"attributes":{"source":{"id":"1084"},"text":{"field":"enames"},"text_color":{"value":"red"},"text_font_size":{"value":"1.5vh"},"text_font_style":{"value":"bold"},"x":{"field":"xe"},"y":{"field":"ye"}},"id":"1085","type":"LabelSet"},{"attributes":{"active_multi":null,"tools":[{"id":"1178"},{"id":"1179"},{"id":"1180"},{"id":"1181"},{"id":"1182"},{"id":"1222"},{"id":"1230"}]},"id":"1184","type":"Toolbar"},{"attributes":{},"id":"1102","type":"EdgesAndLinkedNodes"},{"attributes":{"args":{"graph":{"id":"1033"},"hover":{"id":"1093"}},"code":"\n\tvar nrend = graph.node_renderer.data_source\n\tif (cb_data.index.indices.length &gt; 0) {\n\t\tvar ndx = cb_data.index.indices[0]\n\t\tvar formula_list = nrend.data[\"formula\"]\n\t\t// fallback for cases without any fragmentation, where formula is not defined\n\t    if (formula_list){\n\t\t\tvar formula = formula_list[ndx]\n\t\t} else {\n\t\t\tvar formula = null\n\t\t}\n\t    if (!formula){\n\t\t\thover.tooltips = [[\"tag\",\"@name\"],[\"E\",\"@energy{%.2f}\"]]\n \t\t} else {\n\t\t\thover.tooltips = [[\"tag\",\"@name\"],[\"E\",\"@energy{%.2f}\"],[\"formula\",\"@formula\"]]\n\t\t}\n\t}\n\t"},"id":"1094","type":"CustomJS"},{"attributes":{},"id":"1009","type":"LinearScale"},{"attributes":{"callback":{"id":"1094"},"description":"Node hover","formatters":{"@energy":"printf"},"renderers":[{"id":"1036"}]},"id":"1093","type":"HoverTool"},{"attributes":{"callback":null,"renderers":[{"id":"1036"},{"id":"1040"}]},"id":"1112","type":"TapTool"},{"attributes":{"args":{"counter":[1],"figure":{"id":"1004"}},"code":"\n\t\t// hide all labels from a reaction network\n\t\t// figure - bokeh FIGURE for the network view, counter - inner counter for state, use list for mutability\n\t\t// edit the 3rd and 4th elements of the figure.center array, containing labelsets\n\t\t// 0 - all labels on, 1 - TS labels off, 2 - all labels off\n\t\tvar ct = counter[0]\n\t\tswitch (ct) {\n\t\t\tcase 0:\n\t\t\t\tfigure.center[2].visible = true\n\t\t\t\tfigure.center[3].visible = true\n\t\t\t\tbreak\n\t\t\tcase 1:\n\t\t\t\tfigure.center[2].visible = true\n\t\t\t\tfigure.center[3].visible = false\n\t\t\t\tbreak\n\t\t\tcase 2:\n\t\t\t\tfigure.center[2].visible = false\n\t\t\t\tfigure.center[3].visible = false\n\t\t\t\tbreak\n\t\t}\n\t\tct = ct + 1\n\t\tif (ct &gt;= 3) {\n\t\t\tct = 0\n\t\t}\n\t\tcounter[0] = ct\n\t\t"},"id":"1114","type":"CustomJS"},{"attributes":{"callback":null,"description":"Edge hover","formatters":{"@energy":"printf"},"line_policy":"interp","renderers":[{"id":"1040"}],"tooltips":[["tag","@name"],["E","@energy{%.2f}"]]},"id":"1099","type":"HoverTool"},{"attributes":{"callback":{"id":"1114"},"description":"Hide labels","icon":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAAgAAAAICAIAAABLbSncAAAAKUlEQVR4nGM8ceIEAzbAwsDAYG5ujiZ68uRJJqzKGRgY6CHBAnEDpgQA24sHx1E8XckAAAAASUVORK5CYII="},"id":"1115","type":"CustomAction"},{"attributes":{"data":{"form":[null,null,null,null,"H2 + CO2","H2 + CO2"],"lab":["MIN1","MIN1","TS10","TS10","PR2","PR2"],"x":{"__ndarray__":"AAAAAAAAAAAAAAAAAADwPwAAAAAAAABAAAAAAAAACEAAAAAAAAAQQAAAAAAAABRA","dtype":"float64","order":"little","shape":[6]},"y":[-8.323520999999985,-8.323520999999985,58.154156,58.154156,-37.60977,-37.60977]},"selected":{"id":"1267"},"selection_policy":{"id":"1268"}},"id":"1191","type":"ColumnDataSource"},{"attributes":{"bottom_units":"screen","fill_alpha":0.5,"fill_color":"lightgrey","left_units":"screen","level":"overlay","line_alpha":1.0,"line_color":"black","line_dash":[4,4],"line_width":2,"right_units":"screen","syncable":false,"top_units":"screen"},"id":"1183","type":"BoxAnnotation"},{"attributes":{"height":600,"info":{"height":"100%","j2sPath":"https://chemapps.stolaf.edu/jmol/jsmol/j2s","script":"background white;","serverURL":"https://chemapps.stolaf.edu/jmol/jsmol/php/jsmol.php","use":"HTML5","width":"100%"},"script_source":{"id":"1125"},"width":600},"id":"1126","type":"JSMol"},{"attributes":{"align":"center","max_width":200,"value":"MIN3"},"id":"1129","type":"TextInput"},{"attributes":{"align":"center","icon":null,"js_event_callbacks":{"button_click":[{"id":"1141"}]},"label":"Locate molecule","max_width":400},"id":"1130","type":"Button"},{"attributes":{"align":"center","css_classes":["xspecial"],"icon":null,"js_event_callbacks":{"button_click":[{"id":"1139"}]},"label":"Load vibrations","max_width":200},"id":"1127","type":"Button"},{"attributes":{"end":1.2,"start":-1.2},"id":"1002","type":"Range1d"},{"attributes":{"args":{"graph":{"id":"1033"},"prof":{"id":"1161"}},"code":"\n\t\t// graph - bokeh graph for the network view, use to fetch node and edge renderers\n\t\t// prof - bokeh figure for the profiles\n\t\tvar nrend = graph.node_renderer.data_source\n\t\tvar erend = graph.edge_renderer.data_source\n\t\tvar ninds = nrend.selected.indices\n\t\tvar einds = erend.selected.indices\n\t\tif (ninds.length) {\n\t\t\tvar rend = nrend\n\t\t} else if (einds.length) {\n\t\t\tvar rend = erend\n\t\t} else {\n\t\t\tvar rend = null\n\t\t\t// make everything visible\n\t\t\tfor (let [index,line] of prof.renderers.entries()){\n\t\t\t\tprof.renderers[index].visible = true\n\t\t\t\t// only EVEN entries for center\n\t\t\t\tif (index % 2 == 0){\n\t\t\t\t\tprof.center[index+2].visible = true\n\t\t\t\t}\n\t\t\t} \n\t\t\treturn true\n\t\t}\n\t\tvar ndx = rend.selected.indices[0]\n\t\tvar sel_spc = rend.data[\"name\"][ndx]\n\t\t// labels are stored in \"center\" property of the figure\n\t\t// and we have TWO renderers per entry: skeleton and rectangle\n\t\tvar Nlines = prof.renderers.length/2\n\t\tfor (let i = 0 ; i &lt; Nlines ; i++){\n\t\t\tconst index = 2*i\n\t\t\tvar current_data = prof.renderers[index].data_source\n\t\t\tvar labels = current_data.data[\"lab\"]\n\t\t\t// use ! to negate inclusion\n\t\t\tif (!labels.includes(sel_spc)){\n\t\t\t\tprof.renderers[index].visible = false\n\t\t\t\tprof.renderers[index+1].visible = false\n\t\t\t\t// labels are also duplicated: first tag, then energy.\n\t\t\t\tprof.center[index+2].visible = false\n\t\t\t\tprof.center[index+3].visible = false\n\t\t\t} \n\t\tprof.properties.renderers.change.emit()\n\t\t}\n\t\t"},"id":"1227","type":"CustomJS"},{"attributes":{},"id":"1249","type":"AllLabels"},{"attributes":{},"id":"1162","type":"DataRange1d"},{"attributes":{},"id":"1250","type":"BasicTickFormatter"},{"attributes":{"children":[[{"id":"1155"},0,0,1,1],[{"id":"1158"},0,1,1,1]]},"id":"1159","type":"GridBox"},{"attributes":{"below":[{"id":"1170"}],"center":[{"id":"1173"},{"id":"1177"},{"id":"1202"},{"id":"1205"},{"id":"1217"},{"id":"1220"}],"js_event_callbacks":{"reset":[{"id":"1225"}]},"left":[{"id":"1174"}],"min_border_left":60,"name":"PROFILE","renderers":[{"id":"1195"},{"id":"1199"},{"id":"1210"},{"id":"1214"}],"title":{"id":"1234"},"toolbar":{"id":"1184"},"visible":false,"x_range":{"id":"1162"},"x_scale":{"id":"1166"},"y_range":{"id":"1164"},"y_scale":{"id":"1168"}},"id":"1161","subtype":"Figure","type":"Plot"},{"attributes":{},"id":"1168","type":"LinearScale"},{"attributes":{},"id":"1164","type":"DataRange1d"}],"root_ids":["1159"]},"title":"Bokeh Application","version":"2.3.3"}}
 
        </script>
 
        <script type="text/javascript">
 
          (function() {
 
            var fn = function() {
 
              Bokeh.safely(function() {
 
                (function(root) {
 
                  function embed_document(root) {
 
                   
 
                  var docs_json = document.getElementById('1703').textContent;
 
                  var render_items = [{"docid":"f436dc99-c491-4996-b7ca-3e1293247474","root_ids":["1159"],"roots":{"1159":"3c780241-42fe-478b-8b7d-bb4e3e9dfbe4"}}];
 
                  root.Bokeh.embed.embed_items(docs_json, render_items);
 
               
 
                  }
 
                  if (root.Bokeh !== undefined) {
 
                    embed_document(root);
 
                  } else {
 
                    var attempts = 0;
 
                    var timer = setInterval(function(root) {
 
                      if (root.Bokeh !== undefined) {
 
                        clearInterval(timer);
 
                        embed_document(root);
 
                      } else {
 
                        attempts++;
 
                        if (attempts > 100) {
 
                          clearInterval(timer);
 
                          console.log("Bokeh: ERROR: Unable to run BokehJS code because BokehJS library is missing");
 
                        }
 
                      }
 
                    }, 10, root)
 
                  }
 
                })(window);
 
              });
 
            };
 
            if (document.readyState != "loading") fn();
 
            else document.addEventListener("DOMContentLoaded", fn);
 
          })();
 
        </script>
 
   
 
  </body>
 
 
 
</html>
 
 
 
 
 
 
 
  
 +
The Python library [https://github.com/dgarayr/amk_tools '''amk-tools'''] can be employed to read, process and visualize reaction networks generated by AutoMeKin and is available in the latest singularity container (and through our web interface). Here is the paper describing the new tool: [https://pubs.acs.org/doi/10.1021/acsphyschemau.1c00051 <span style="font-size:100%"> Garay-Ruiz et al. ACS Phys. Chem. Au 2022]
  
 
AutoMekin got the [https://geqc.rseq.org/emilio-martinez-nunez-and-mireia-via-nadal-recieve-the-geqc-rseq-2021-methodology-development-prize/ GEQC-RSEQ 2021 Methodology Development Prize].
 
AutoMekin got the [https://geqc.rseq.org/emilio-martinez-nunez-and-mireia-via-nadal-recieve-the-geqc-rseq-2021-methodology-development-prize/ GEQC-RSEQ 2021 Methodology Development Prize].

Latest revision as of 13:15, 11 March 2022

Introduction

Automekin.png

AutoMeKin (formerly tsscds) has been designed to discover reaction mechanisms in an automated fashion. Transition states are located using MD simulations and Graph Theory algorithms. Monte Carlo simulations afford kinetic results. The only input is a starting structure in XYZ format. The method is described in these two publications: 1 2. At present MOPAC2016, Entos Qcore and Gaussian (G09/G16) are interfaced with AutoMeKin. The program has been tested on the following Linux distros: CentOS 7, Red Hat Enterprise Linux and Ubuntu 20.04 LTS.

To give you a flavor of the capabilities of the program you can try our web interface


Web.png

Developers team

Emilio Martinez-Nunez, George L. Barnes, Carles Bo, Diego Garay-Ruiz, David R. Glowacki, Sabine Kopec, Daniel Pelaez-Ruiz, Aurelio Rodriguez, Roberto Rodriguez-Fernandez, Robin J. Shannon, James J. P. Stewart, Pablo G. Tahoces and Saulo A. Vazquez

Address
Departamento de Química Física, Facultade de Química, Avda. das Ciencias s/n, 15782 Santiago de Compostela, SPAIN
email me

License

MIT License

Copyright (C) 2021 AutoMeKin

Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

Install the code

To install the code follow the Installation instructions

Tutorial

Download tutorial

Program execution

Unless you donwloaded the singularity container (in that case skip this step), to start using any of the scripts of the program, load the amk/2021 module:

module load amk/2021

To run the low-level calculations use:

nohup llcalcs.sh molecule.dat ntasks niter runningtasks >llcalcs.log 2>&1 &

where:
molecule is the name of your molecule
ntasks is the number of tasks
niter is the number of iterations
runningtasks is the number of simultaneous tasks

To run the high-level calculations use:

nohup hlcalcs.sh molecule.dat runningtasks >hlcalcs.log 2>&1 &

For more details, follow the instructions given in the tutorial

Reference

If you use AutoMeKin, please cite the following publications:

  1. E Martínez-Nuñez, G.L. Barnes, D.R. Glowacki, S. Kopec, D. Pelaez, A. Rodriguez, R. Rodriguez-Fernandez, R.J. Shannon, J.J.P. Stewart, P.G. Tahoces, S.A. Vazquez, J. Comput. Chem. 2021, 42, 2036
  2. E. Martínez-Núñez, J. Comput. Chem. 2015, 36, 222
  3. E. Martínez-Núñez, Phys. Chem. Chem. Phys. 2015,17, 14912
  4. J. J. P. Stewart, MOPAC2016, Stewart Computational Chemistry: Colorado Springs, CO, USA, HTTP://OpenMOPAC.net, 2016.

If you use the older version tsscds, please cite, instead of ref 1 above, the following one:

A. Rodriguez, R. Rodriguez-Fernandez, S.A. Vazquez, G.L. Barnes, J.J.P. Stewart, E Martínez-Nuñez, J. Comput. Chem. 2018, 39, 1922

Works that employ AutoMeKin

The following works employ AutoMeKin: Works

Changelog

Consult the Latest changes

Web interface

AutoMeKin can be used through our web interface.

News

AutoMeKin is now interfaced with Gaussian 16 (from rev. 1131)

The Python library amk-tools can be employed to read, process and visualize reaction networks generated by AutoMeKin and is available in the latest singularity container (and through our web interface). Here is the paper describing the new tool: Garay-Ruiz et al. ACS Phys. Chem. Au 2022

AutoMekin got the GEQC-RSEQ 2021 Methodology Development Prize.

AutoMeKin now includes barrierless channels. However, these channels are not considered in the kinetics.

AutoMeKin has been interfaced with Entos Qcore.

AutoMeKin has been interfaced with BXDE to enhance its efficiency ( R. A. Jara-Toro et al. ChemSystemsChem doi: 10.1002/syst.201900024).

The method has also been recently generalized in a collaboration with Dani Pelaez and co-workers to study van der Waals structures ( S. Kopec et al. Int. J. Quantum Chem. 2019, 119, e26008 ) and also to generate sum-of-products PESs for quantum dynamics ( R. Panades-Barrueta et al. Frontiers in Chemistry 2019, 7, 576).

Return to Contents

Return to Main_Page