ElemCo.jl <img align="right" src="files/coil.png" height=74> <br/> <span style="font-weight:normal"> elemcoil</span><br/>

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

| Class | Methods | |-----------------|---------------------------------------| | Mean-field | DF-HF, DF-UHF, BO-HF, BO-UHF, DF-MCSCF| | Perturbative | DF-MP2, MP2, UMP2, SOS-LT-DF-MP2 | | Coupled Cluster | CCSD, RCCSD, UCCSD, DCSD, RDCSD, UDCSD, CCSD(T), RCCSD(T), UCCSD(T), ΛCCSD(T), ΛUCCSD(T), FR-CCSD, FR-DCSD, 2D-CCSD, 2D-DCSD, CCSDT, UCCSDT, DC-CCSDT, UDC-CCSDT, FR-CCSDT, FR-DC-CCSDT, SVD-DF-DCSD, SVD-DC-CCSDT | | CI | FCI, CIPHI | | DMRG | ITensors.jl interface | | Excited states | EOM-CCSD, EOM-DCSD, FCI, CIPHI |

Getting started

Requirements: julia (>1.9)

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>'

Start jlm in the terminal and run the following command to install the dependencies,

julia> using Pkg
julia> Pkg.instantiate()

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.

Input generation and orbital visualizer

ElemCo.jl input files can be generated using jlmol, which can also be used to visualize the molecular orbitals. The browser version of jlmol can be found at app.jlmol.com.

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.

Various options can be set using the @set macro. It is advisable to set options locally for each calculation to avoid unintended side effects. For example, to change the convergence threshold of the Hartree-Fock calculation to 1e-8, the following line can be added before the @dfhf macro:

@set scf thr=1e-8

or locally for the @dfhf macro as:

@dfhf begin
  @set scf thr=1e-8
end

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