ElemCo.jl

Julia implementation of various electron-correlation methods (main focus on coupled cluster methods).
Author fkfest
Popularity
8 Stars
Updated Last
2 Months Ago
Started In
February 2023

ElemCo.jl

Julia implementation of various electron-correlation methods (main focus on coupled cluster methods). The integrals are obtained from a FCIDUMP file or calculated using an interface to libcint.

Capabilities

Canonical DF Only closed-shell
RHF ✔️
UHF ✔️
BO-HF ✔️
MCSCF 🔧
MP2 ✔️
CCSD ✔️ ✔️
RCCSD ✔️
UCCSD ✔️
ΛCCSD ✔️ ✔️
ΛUCCSD ✔️
CCSD(T) ✔️ ✔️
UCCSD(T) ✔️
ΛCCSD(T) ✔️ ✔️
ΛUCCSD(T) ✔️
FR-CCSD ✔️
2D-CCSD ✔️
DCSD ✔️ ✔️
RDCSD ✔️
UDCSD ✔️
ΛDCSD ✔️ ✔️
ΛUDCSD ✔️
FR-DCSD ✔️
2D-DCSD ✔️
SVD-DCSD ✔️ ✔️
SVD-DC-CCSDT ✔️ ✔️
DMRG ✔️

Getting started

Requirements: julia (>1.8)

Usage

Installation

To install ElemCo.jl, run the following command in the Julia REPL,

julia> using Pkg
julia> Pkg.add("ElemCo")

For a development version of ElemCo.jl, clone the ElemCo.jl-devel repository and create an alias to set the project to the ElemCo.jl directory,

alias jlm='julia --project=<path_to_ElemCo.jl>'

Now the command jlm can be used to start the calculations,

jlm input.jl

Input file

The input file is a Julia script that contains the calculation details. The script should start with the following lines,

using ElemCo
@print_input

The @print_input macro prints the input file to the standard output. The calculation details are specified using the macros provided by ElemCo.jl.

Macros

The following macros are available in ElemCo.jl (see the documentation for more details and macros),

  • @dfhf - Performs a density-fitted Hartree-Fock calculation.
  • @cc <method> - Performs a coupled cluster calculation.
  • @dfcc <method> - Performs a coupled cluster calculation using density fitting.
  • @set <option> <setting> - Sets the options for the calculation.

etc.

Default scratch directory path on Windows is the first environment variable found in the ordered list TMP, TEMP, USERPROFILE. On all other operating systems TMPDIR, TMP, TEMP, and TEMPDIR. If none of these are found, the path /tmp is used. Default scratch folder name is elemcojlscr.

Variable names fcidump, geometry and basis are reserved for the file name of FCIDUMP, geometry specification and basis sets, respectively.

Example

DCSD calculation using integrals from a FCIDUMP file

The ground state energy can be calculated using the DCSD method with the following script:

using ElemCo
@print_input
fcidump = "../test/H2O.FCIDUMP"
@cc dcsd

DCSD calculation of the water molecule using density-fitted integrals

In order to calculate the ground state energy of the water molecule using the DCSD method, the following script can be used:

using ElemCo
@print_input
geometry="bohr
     O      0.000000000    0.000000000   -0.130186067
     H1     0.000000000    1.489124508    1.033245507
     H2     0.000000000   -1.489124508    1.033245507"

basis = "vdz"
@dfhf
@cc dcsd

The @dfhf macro calculates the density-fitted Hartree-Fock energy and orbitals and then DCSD calculation is performed using density-fitted integrals.

Further example scripts are provided in the examples directory.

Documentation is available at https://elem.co.il.

Electron coil
A poem by Bing

In the heart of an atom lies a tiny core
Where protons and neutrons are tightly bound
But around this nucleus, there's so much more
A cloud of electrons that swirls around

They don't orbit in circles like planets do
But jump and spin in quantum states
They can be here and there and everywhere too
And sometimes they even change their mates

This is the electron coil, the turmoil of the shell
The source of light and heat and power
The force that makes the atoms repel or gel
The spark that ignites the cosmic flower

Used By Packages

No packages found.