How a Single Point Energy Calculation Works¶
A Single Point Energy calculation determines the electronic energy of a molecule while keeping the positions of all nuclei fixed. Unlike a geometry optimization, the atomic coordinates are not modified during the calculation. Instead, Gaussian solves for the most stable distribution of electrons corresponding to the supplied molecular geometry.
The calculation begins with the molecular structure provided in the input file. Using the selected computational method and basis set, Gaussian determines the electronic wavefunction and calculates the total electronic energy of the molecule.
The entire process is performed without moving any atoms, making the Single Point Energy calculation one of the fastest and most fundamental calculations available in Gaussian.
Overall Workflow¶
The sequence of operations performed during a Single Point Energy calculation can be summarized as
Input Molecular Geometry
│
▼
Read Atomic Coordinates
│
▼
Construct Basis Functions
│
▼
Generate Initial Electron Density
│
▼
Perform SCF Iterations
│
▼
Check Convergence
│
▼
Compute Molecular Properties
│
▼
Print Final Electronic Energy
Each of these steps is explained below.
Step 1: Reading the Molecular Geometry¶
The calculation begins by reading the Cartesian coordinates supplied in the input file.
For example,
These coordinates define the positions of all nuclei in three-dimensional space.
During a Single Point Energy calculation, these coordinates remain unchanged from the beginning to the end of the calculation.
Step 2: Constructing the Basis Functions¶
Instead of solving the electronic Schrödinger equation exactly, Gaussian represents each molecular orbital as a combination of mathematical functions known as basis functions.
The chosen basis set determines
- how accurately the electron density is described,
- the computational cost,
- and the overall quality of the calculation.
For this tutorial, the 6-31G basis set is used.
Step 3: Initial Electron Density Guess¶
Before solving for the electronic structure, Gaussian must generate an initial estimate of where the electrons are likely to be located.
This initial guess is constructed automatically using built-in procedures based on the atomic composition of the molecule.
Although the first estimate is only approximate, it serves as the starting point for the iterative solution.
Step 4: Self-Consistent Field (SCF) Iterations¶
The heart of the calculation is the Self-Consistent Field (SCF) procedure.
Starting from the initial electron density, Gaussian repeatedly performs the following cycle:
- Construct the electronic Hamiltonian (or Kohn–Sham matrix for DFT).
- Solve for the molecular orbitals.
- Generate a new electron density.
- Compare the new density with the previous one.
- Repeat the process until the changes become sufficiently small.
Each repetition is called an SCF iteration.
As the iterations proceed, the electronic energy gradually approaches a stable value.
Step 5: Convergence of the Electronic Structure¶
Gaussian continuously monitors the change in electronic energy and electron density between successive SCF iterations.
Once these changes fall below predefined convergence thresholds, the calculation is considered converged.
At this point, the electronic wavefunction is regarded as self-consistent, meaning that the electron density used to construct the Hamiltonian is the same as the density obtained after solving it.
Step 6: Calculating Molecular Properties¶
After SCF convergence has been achieved, Gaussian evaluates additional molecular properties using the converged wavefunction.
Depending on the requested keywords, these properties may include
- Total electronic energy
- Molecular orbital energies
- Dipole moment
- Electron density
- Mulliken charges
- Electrostatic potential
- Population analysis
For a basic Single Point Energy calculation, the primary quantity of interest is the total electronic energy.
Step 7: Writing the Output File¶
Finally, Gaussian writes all computed information to the output file.
The output contains
- Input summary
- SCF convergence history
- Final electronic energy
- Molecular properties
- Timing information
- Job completion status
If the calculation finishes successfully, the output ends with
Why Are the Atoms Not Moving?¶
One of the most common questions asked by beginners is why the molecular geometry remains unchanged.
The answer is simple:
A Single Point Energy calculation is designed only to evaluate the electronic energy of a given molecular structure.
Since no optimization keyword (such as Opt) is present in the input file, Gaussian performs no force calculations to update the atomic positions. The nuclei remain fixed while only the electrons are allowed to relax.
Relationship to Other Gaussian Calculations¶
Many other Gaussian calculations are built upon repeated Single Point Energy calculations.
For example,
Single Point Energy
│
├────────► Electronic Energy
│
├────────► Geometry Optimization
│ │
│ └── Many Single Point Calculations
│
├────────► Frequency Analysis
│
├────────► IRC Calculation
│
├────────► Coordinate Scan
│
└────────► Transition-State Search
A geometry optimization, for example, repeatedly performs Single Point Energy calculations at different molecular geometries until the lowest-energy structure is found.
For this reason, understanding the Single Point Energy calculation provides the foundation for understanding nearly every other calculation available in Gaussian.
Summary¶
A Single Point Energy calculation evaluates the electronic energy of a molecule for a fixed molecular geometry. Gaussian begins by reading the atomic coordinates, constructs the basis functions, generates an initial electron density, and then performs Self-Consistent Field iterations until the electronic structure converges. Once convergence is achieved, the program computes the requested molecular properties and writes the results to the output file.
Although conceptually simple, the Single Point Energy calculation is the building block for almost all electronic structure calculations performed in Gaussian.