Chapter 2 Notes – Computational Methods & Tools

Chemicals & Purity

• 16 key reagents (graphite, strong acids/bases, pesticides, antibiotics, amino acids) with ≥97 % purity from Sigma-Aldrich, Merck, CDH, TCI, Avra.

Quantum-Mechanical Foundation

• Many-electron behaviour governed by time-independent Schrödinger equation: $-\frac{\hbar^2}{2m}\nabla^2\Psi + V\Psi = E\Psi$.

• Exact solutions feasible only for one-electron systems; approximations mandatory for nanoscale modelling.

Electronic Structure Methods

• Semi-empirical (AM1, PM3): parameter-based, fast, limited accuracy.

• Ab initio (no empirical data):

  • Hartree–Fock (HF): single Slater determinant; omits electron correlation.

  • Post-HF improvements:
    • Møller-Plesset perturbation $\text{MP}n$ (commonly $\text{MP}2$).
    • Configuration Interaction (CI): linear combination of excited determinants.
    • Coupled Cluster (CC): exponential ansatz $\Psi = e^{\hat{T}}\Psi_0$.

Density Functional Theory (DFT)

• Uses total electron density $n(\mathbf{r})$; avoids many-electron wavefunction.

• Hohenberg–Kohn theorems: one-to-one mapping between $n(\mathbf{r})$ and external potential; ground-state energy functional $E[n]$.

• Kohn–Sham (KS) formulation:

  • Functional split: $F[n]=T<em>S[n]+E</em>H[n]+E_{XC}[n]$.

  • KS equations: $\left[-\frac{1}{2}\nabla^2 + v<em>{KS}(\mathbf{r})\right]\psi</em>i = \varepsilon<em>i\psi</em>i$, with $v<em>{KS}=v</em>{ext}+v<em>H+v</em>{XC}$.

Exchange–Correlation Functionals

• LDA/LSDA: density-only; e.g., Slater + VWN.

• GGA: density + gradient; popular pairs BLYP, BP86, PW91, PBE.

• Hybrids (HF exchange + DFT): B3LYP, B3P86, mPW1PW91, PBE0.

Basis Sets

• Minimal: STO-nG (e.g., STO-3G).

• Split-valence: 3-21G, 6-31G, 6-311G.

• Polarization: * / (d,p) indicators (e.g., 6-311G(d,p)).

• Diffuse: + / ++ or aug- (e.g., 6-31+G(d), aug-cc-pVDZ).

• Correlation-consistent: cc-pVNZ (N = D,T,Q,5,6).

Gaussian Software

• Gaussian09 with GaussView GUI employed.

• Delivers optimized geometries, energies, vibrational frequencies, transition states via HF, DFT, semi-empirical and post-HF methods.

Solvent Models

• Implicit/continuum: PCM (IEF-PCM), SMx, COSMO.

• Explicit: discrete solvent molecules via MD or MC simulations.

Analysis & Characterization Tools

• Molecular Electrostatic Potential (MEP): maps charge distribution.

• Natural Bond Orbital (NBO): localized lone-pair and bond analysis.

• Intrinsic Reaction Coordinate (IRC): traces minimum-energy path between transition state and minima.

Chemicals & Purity

• 16 key reagents, including high-purity graphite, strong acids (e.g., HCl, H2SO4) and bases (e.g., NaOH, KOH), various pesticides, antibiotics, and amino acids. These are sourced from reputable suppliers like Sigma-Aldrich, Merck, CDH, TCI, and Avra, ensuring a purity of $\ge$97 %. This high purity is crucial for experimental accuracy and reliable calibration in quantum chemical studies, as impurities can significantly alter experimental outcomes and computational models.

Quantum-Mechanical Foundation

• The fundamental behavior of many-electron systems in atoms and molecules is governed by the time-independent Schrödinger equation: $-\frac{\hbar^2}{2m}\nabla^2\Psi + V\Psi = E\Psi$. Here, $\Psi$ represents the many-body wavefunction, which encapsulates all information about the system's quantum state, including the positions and momenta of all electrons and nuclei. Unfortunately, exact analytical solutions to this equation are only feasible for very simple one-electron systems (like the hydrogen atom). For any system with multiple electrons, the complexity of electron-electron interactions makes exact solutions computationally intractable, necessitating the use of various approximation methods for accurate nanoscale modelling of molecular properties and reactions.

Electronic Structure Methods

• Semi-empirical methods (e.g., AM1, PM3): These methods simplify calculations by incorporating empirical parameters, derived from experimental data, to approximate complex integrals. While significantly faster than ab initio methods, their accuracy is limited by the quality and transferability of these parameters, often performing well only for systems similar to those they were parameterized against.

• Ab initio methods (from first principles, no empirical data):

  • Hartree–Fock (HF): This is the foundational ab initio method, which approximates the many-electron wavefunction as a single Slater determinant (an antisymmetrized product of one-electron orbitals). HF considers the average repulsion between electrons but entirely omits instantaneous electron correlation, meaning it doesn't account for the fact that electrons avoid each other due to their instantaneous positions. This omission leads to an overestimation of electron-electron repulsion and an underestimation of binding energies.

  • Post-HF improvements: These methods build upon the HF approximation to systematically incorporate electron correlation, leading to higher accuracy:

    • Møller-Plesset perturbation theory ($\text{MP}n$): Accounts for correlation effects perturbatively, with $\text{MP}2$ being the most common. It adds corrections to the HF energy based on small perturbations from the HF solution, improving accuracy significantly for many systems.

    • Configuration Interaction (CI): Constructs the many-electron wavefunction as a linear combination of the HF ground state and various excited determinants. By including these

Chemicals & Purity

• 16 key reagents, including high-purity graphite, strong acids (e.g., HCl, H2SO4) and bases (e.g., NaOH, KOH), various pesticides, antibiotics, and amino acids. These are sourced from reputable suppliers like Sigma-Aldrich, Merck, CDH, TCI, and Avra, ensuring a purity of $\ge$97 %. This high purity is crucial for experimental accuracy and reliable calibration in quantum chemical studies, as impurities can significantly alter experimental outcomes and computational models.

Quantum-Mechanical Foundation

• The fundamental behavior of many-electron systems in atoms and molecules is governed by the time-independent Schrödinger equation: $-\frac{\hbar^2}{2m}\nabla^2\Psi + V\Psi = E\Psi$. Here, $\Psi$ represents the many-body wavefunction, which encapsulates all information about the system's quantum state, including the positions and momenta of all electrons and nuclei. Unfortunately, exact analytical solutions to this equation are only feasible for very simple one-electron systems (like the hydrogen atom). For any system with multiple electrons, the complexity of electron-electron interactions makes exact solutions computationally intractable, necessitating the use of various approximation methods for accurate nanoscale modelling of molecular properties and reactions.

Electronic Structure Methods

• Semi-empirical methods (e.g., AM1, PM3): These methods simplify calculations by incorporating empirical parameters, derived from experimental data, to approximate complex integrals. While significantly faster than ab initio methods, their accuracy is limited by the quality and transferability of these parameters, often performing well only for systems similar to those they were parameterized against.

• Ab initio methods (from first principles, no empirical data):

  • Hartree–Fock (HF): This is the foundational ab initio method, which approximates the many-electron wavefunction as a single Slater determinant (an antisymmetrized product of one-electron orbitals). HF considers the average repulsion between electrons but entirely omits instantaneous electron correlation, meaning it doesn't account for the fact that electrons avoid each other due to their instantaneous positions. This omission leads to an overestimation of electron-electron repulsion and an underestimation of binding energies.

  • Post-HF improvements: These methods build upon the HF approximation to systematically incorporate electron correlation, leading to higher accuracy:

    • Møller-Plesset perturbation theory ($\text{MP}n$): Accounts for correlation effects perturbatively, with $\text{MP}2$ being the most common. It adds corrections to the HF energy based on small perturbations from the HF solution, improving accuracy significantly for many systems.

    • Configuration Interaction (CI): Constructs the many-electron wavefunction as a linear combination of the HF ground state and various excited determinants. By including these