Course Material
Online tutorial:
- Introduction to Jupyter notebooks: ccec_notebook_basics.ipynb
- Notebook with solutions: ccec_notebook_basics_instructor.ipynb
- Slides: CCEC_introduction_to_notebooks.pdf
Monday:
- Slides: CCEC_welcome.pdf
- Basis sets: 1_basis-sets_orbitals.ipynb
- Hartree–Fock: 2_hartree-fock.ipynb
- Density Functional Theory (DFT): 3_dft.ipynb
- Environment, QM/MM: 4_qmmm_cpcm.ipynb
- Visualization: 5_visualization.ipynb
Tuesday:
- Session1_Thermochemistry_ReactionEnergies.ipynb
- Session2_Thermochemistry_Conformers.ipynb
- Session3_Electrochemistry.ipynb
- Session4_MolecularProperties.ipynb
Wednesday:
- Optimization: 1_opt.ipynb
- Vibrational analysis: 2_vib.ipynb
- Experimental IR spectrum: ethane_ir.csv
- Transition state search: 3_ts.ipynb
Thursday
Code snippet:
import math
chlorine = vlx.Molecule.read_smiles("[Cl]")
chlorine.get_masses()
chlorine.get_masses()
kB = 1.380649e-23
NA = 6.022141e23
T = 298.15
h = 6.626070e-34
P = 1e5
M = chlorine.get_masses()[0]*0.001
mass = M/NA # in kilos
NA = 6.022141e23
T = 298.15
h = 6.626070e-34
P = 1e5
M = chlorine.get_masses()[0]*0.001
mass = M/NA # in kilos
Htrans = 2.5*kB*NA*T
print(Htrans, "in J/mol")
print(Htrans, "in J/mol")
# Sackur-Tetrode equation
Strans = kB * NA * (math.log((2 * math.pi * mass * kB * T / (h**2))**1.5 * (kB * T / P)) + 2.5)
print(Strans, "in J/(K*mol)")
print(Strans, "in J/(K*mol)")