unrestricted.cpp File Reference

#include <electronsystems/gaussian/gaussiancore.h>
#include <electronsystems/gaussian/gaussiansystem.h>
#include <solvers/unrestrictedhartreefocksolver.h>
#include <iostream>
#include <unittest++/UnitTest++.h>

Include dependency graph for unrestricted.cpp:

Go to the source code of this file.

Functions
	SUITE (Unrestricted)

Function Documentation

SUITE

(

Unrestricted

)

Definition at line 12 of file unrestricted.cpp.

                    {
    TEST(HydrogenMolecule) {
        vector<GaussianCore> cores;
        cores.push_back(GaussianCore({0,0,0}, "atom_1_basis_3-21G.tm"));
        cores.push_back(GaussianCore({1.4,0.0,0}, "atom_1_basis_3-21G.tm"));
        GaussianSystem system;
        for(const GaussianCore &core : cores) {
            system.addCore(core);
        }
        UnrestrictedHartreeFockSolver solver(&system);
        solver.setConvergenceTreshold(1e-12);
        solver.setNIterationsMax(1e3);
        solver.setDensityMixFactor(0.5);
        solver.solve();
        CHECK_CLOSE(-1.122933363617109, solver.energy(), 1e-6);
    }

    TEST(HydrogenMoleculeLongDistance) {
        vector<GaussianCore> cores;
        cores.push_back(GaussianCore({0,0,0}, "atom_1_basis_4-31G.tm"));
        cores.push_back(GaussianCore({100.0,0,0}, "atom_1_basis_4-31G.tm"));
        GaussianSystem system;
        for(const GaussianCore &core : cores) {
            system.addCore(core);
        }
        UnrestrictedHartreeFockSolver solver(&system);
        solver.setConvergenceTreshold(1e-9);
        solver.setNIterationsMax(1e3);
        solver.setDensityMixFactor(0.5);
        solver.solve();
        CHECK_CLOSE(-0.9964658189712963, solver.energy(), 1e-6);
    }
    TEST(Water) {
        vector<GaussianCore> cores;
        cores.push_back(GaussianCore({0,0,0}, "atom_8_basis_4-31G.tm"));
        cores.push_back(GaussianCore({-1.43,1.108,0}, "atom_1_basis_4-31G.tm"));
        cores.push_back(GaussianCore({1.43,1.108,0}, "atom_1_basis_4-31G.tm"));
        GaussianSystem system;
        for(const GaussianCore &core : cores) {
            system.addCore(core);
        }
        UnrestrictedHartreeFockSolver solver(&system);
        solver.setConvergenceTreshold(1e-12);
        solver.setNIterationsMax(1e3);
        solver.setDensityMixFactor(0.5);
        solver.solve();
        CHECK_CLOSE(-75.90736859918989, solver.energy(), 1e-6);
    }
}