Laboratory 4 ORCA Avogadro and Data analysis

Pagina 1

DATA ANALYSIS In order to plot only 2 columns we can use the indexing, but if you want ONLY IMPORT the second and third column of your text file, you use the loadtxt construct within numPy, and using the construct usecols. Moving to carrion data, here is how to plot the 0 column as x axis and 1 column (the second one) as y axis Let’s move to the OPTIMIZE module. The callable function we enter in the curve_fit function is created with “model”. NOTE THAT THE X VARIABLE, IN OUR CASE THE T, MUST COME FIRST IN THE DEFINITION OF THE MODEL, OTHERWISE THE OPTIMIZATION WILL BE DONE WRONGLY. In such a way, the function knows that the parameters a and k are the ones to be optimized. Calculate basic statistics (see Notebook) ORCA We will perform some electronic calculations on Carbon Monoxide. In order to invoke orca, you do it with the Anaconda Poweshell, not double clicking on it ( otherwise you can use the prompt). Remember to set the correct directory, or type the full path, or adding the the environmental variables the path. (you will notice that removing the path, orca.exe will not be recognized once you launch it from the console). Now, in the temp directory (but any directory is fine) you just create a file. mkdir temp  cd temp  code CO-rfh.inp It will open this file in VS studio. Now , we’re going to copy the code lines. SAVE it, and now the .inp file has been successfully created! # A simple single point RHF calculation ! RHF def2-svp LargePrint # Carbon Monoxide at the experimental geometry * xyz 0 1 C 0 0 0 O 0 0 1.128 * Code Let’s look at the code. The first line starts with * ! RHF def2-svp LargePrint The fist line (in red) tells the program what type of calculation to do. In our case RHF stands for restricted Harthree-Fock Method. def2-svp stands for split valence…