(sec:orca_2json.detailed)=
# orca_2json

This utility program supports the exchange of external ORCA data like
geometry, orbitals and basisets but also of internal ORCA data like
1-electron and 2-electron integals with other programs.

(sec:orca_2json.exportData.detailed)=
## Export ORCA data

The program reads information like geometries, basis sets, MOs etc.
stored in the `.gbw` file or other equivalent ones as `.uno`, `.mp2nat`,
`.qro` etc. and calculates integrals to export them in JSON standard
output formats. For density information the `.densities` file must also
be available. The program is called as a standalone via command line.

**Syntax:**
```orca
orca_2json BaseName.gbw -options
or
orca_2json BaseName.mp2nat -options
or
orca_2json BaseName.uno -options
```


The following ASCII and binary JSON-formats are available as command
line options. It is possible to specify more than one format option.

```orca
  -json             Write ASCII JSON file (default) 
  -bson             Write binary JSON file 
  -ubjson           Write universal binary JSON specification file
  -msgpack          Write MessagePack file
```

In addition two more options are availble. The first of them is used to translate a basename.property.txt file to a corresponding one in JSON format (see {ref}`sec:orca_2json.availableinformation.propertyFile.detailed`).

```orca
   -property        Translate a *.txt property file to a *.json one
```

Finally, orca_2json has the ability to create a .gbw file from a json file. For this one needs to use the '-gbw' option (see {ref}`sec:orca_2json.importData.detailed`).

```orca
   -gbw        Create a GBW file from a json one
```

(sec:orca_2json.configurationFiles.detailed)=
## Configuration file

The data stored in the json file can be configured more individually.
Some information like atom information, geometry and charge are always
written in the outputfile (see {ref}`sec:orca_2json.availableinformation.basic.detailed`). Other data can be requested via keywords in a
JSON-formatted configuration file that either exists for every ORCA
output file or for all files in a directory. Without these configuration
files **ALL** available data are stored except for the densities and the
integrals. Because of the huge amount of data these are only available
when explicitly requested in the configuration file.


You can specify a basename-dependend configuration file

```orca
BaseName.json.conf
```

or a global file used for all requests in a directory.

```orca
orca.json.conf
```

Most keywords in the configuration file can be activated or deactivated
with true or false but some keywords like densities or output formats
have more options and require a list of values. If an option is not
specifically selected it is omitted.

**Structure of the configuration file**:

```orca
{
	"keyword": true/false,
	...
	"keyword": true/false,
	"keyword": ["option",...."option"]		
	...
}
```

**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_conffile.conf
:language: orca
```
Using the above configuration file in the working directory, *orca_2json* will not  export the molecular coefficients and the basis set information.

**Example**

Here is an example configuration file with most available keywords where
everything is disabled except for the basis set information and the specified Integrals, all
densities stored in the density file are requested and the output format
should be ascii json and binary json.

```orca
{
	"MOCoefficients": false,
	"Basisset": true,
	"1elPropertyIntegrals": ["dipole","quadrupole","velosity", "printLinMom", "angular_momentum", "higherMoment"],
	"1elPropertyRelIntegrals": ["dipole","quadrupole"],
	"1elIntegrals": ["H","S", "T", "V", "HMO"],
	"1elIntegralsRel": ["H","S","T", "V"],
	"Vaux": false,
	"AuxBasisType": "AuxC"
	"FullTrafo": false,
	"OrbWin": [0,0,0,0,0,0,0,0].
	"2elIntegrals": [ "MO_IJKL", "RI_IAV", "RI_IJKL"],
	"2elNonRedIntegrals": false,
	"2elNonRedRIIntegrals": false,
	"MullikenCharge": false,
	"LoewdinCharge": false,
	"Densities": ["all"],   <--- here you specify the names like "scfp" or "scfr" or all
	"JSONFormats": ["json", "bson"]
}
```

(sec:orca_2json.availableproperties.detailed)=
## Available information

(sec:orca_2json.availableinformation.propertyFile.detailed)=
### Property File

Beside all information that we will see orca_2json can create, it can also translate the property file of ORCA (basename.property.txt) to a JSON file. This option gives access to all properties stored in the property file. (For more information on property file see {ref}`sec:property_file.detailed`).

The **syntax** is:

```orca
orca_2json basename -property
```

*Where*

*basename* is the name of the property file **without** the extension *property.txt*.

**Example**

If we use the following ORCA input (with the name test.inp): 

```{literalinclude} ../../orca_working_input/orca_2json_examples/orca_2json_example_property_file_input.inp
:language: orca
```

ORCA will create and store on disk, a file named *"test.property.txt"*. The start of the file will look like this:

```{literalinclude} ../../orca_working_input/orca_2json_examples/orca_2json_example_property_file_input.property.txt
:language: orca
```

Then running *orca_2json* in the following way:

```orca
orca_2json test -property
```

ORCA will create a new file on disk, named *"test.property.json"*. The start of this file will look like this:

```{literalinclude} ../../orca_working_input/orca_2json_examples/orca_2json_example_property_file_input.property.json
:language: json
```

(sec:orca_2json.availableinformation.basic.detailed)=
### Basic Information
Some basic information will always be written into the JSON-file as

|              |     |                 |
|:-------------|:---:|:----------------|
| per Atom     | \-  | Coords          |
|              | \-  | ElementLabel    |
|              | \-  | ElementNumber   |
|              | \-  | Idx             |
|              | \-  | NuclearCharge   |
| per molecule | \-  | BaseName        |
|              | \-  | Charge          |
|              | \-  | CoordinateUnits |
|              | \-  | HFTyp           |
|              | \-  | Multiplicity    |
|              | \-  | PointGroup      |

For example using the following configuration file where we set everything to false:
```{literalinclude} ../../orca_working_input/orca2json_example_basicInfo.conf
:language: json
```
 will still produce a json file which for the case of a H2 molecule should look like:
 ```{literalinclude} ../../orca_working_input/orca2json_example_basicInfo.json
:language: json
```

(sec:orca_2json.availableinformation.densities.detailed)=
### Densities

*orca_2json* can also export calculated densities in json format.

|                 |     |                                               |
|:----------------|:---:|:----------------------------------------------|
| Densities       | \-  | densities as available in the .densities file |

**Syntax**
*"Densities" : [list of densities]*

Where *list of densities* should be a list of strings with the expected densities.

**NOTE** By default densities, due to their potential size, are not exported to a json file.

**NOTE** An empty bracket syntax (*"Densities" : []*) will cause the program to crash.

**NOTE** There is the string "All" available where the program will export all available densities.

(sec:orca_2json.availableinformation.electronintegrals.detailed)=
### Electron Integrals 

The list of available electron integrals is shown in the next table.

|                         |     |                                               |
|:------------------------|:---:|:----------------------------------------------|
| 1elIntegrals            | \-  | 1-electron integrals                          |
| 1elPropertyIntegrals    | \-  | 1-electron property integrals                 |
| 1elIntegralsRel         | \-  | relativistic 1-electron integrals             |
| 1elPropertyRelIntegrals | \-  | relativistic 1-electron property integrals    |
| 2elIntegrals            | \-  | two-electron integrals                        |
| 2elNonRedIntegrals      | \-  | non-redundant two-electron integrals          |
| 2elNonRedRIIntegrals    | \-  | non-redundant two-electron RI integrals       |

(sec:orca_2json.availableinformation.electronintegrals.1el.detailed)=
#### 1-electron integrals

For 1-electron integrals we use the following notation:
|     |     |                                 |
|:----|:---:|:--------------------------------|
| H   | \-  | one electron matrix             |
| HMO | \-  | one electron matrix in MO basis |
| S   | \-  | overlap matrix                  |
| T   | \-  | kinetic energy matrix           |
| V   | \-  | nuclear attraction matrix       |

**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_1el.json.conf
:language: json
```

will produce a json file were only the H-Matrix and the Overlap matrix are printed (beside the basic information). Please note that for the one electron relativistic integrals there is a separate variable (see {ref}`sec:orca_2json.availableinformation.electronintegrals.1elrel.detailed`)

(sec:orca_2json.availableinformation.electronintegrals.1elproperty.detailed)=
#### 1-electron property integrals
Also available are 1-electron property integrals.

|                      |     |                                               |
|:---------------------|:---:|:----------------------------------------------|
| 1elPropertyIntegrals | \-  | 1-electron property integrals                 |

Currently the following options are valid:
|                   |     |                             |
|:------------------|:---:|:----------------------------|
| angular_momentum  | \-  | Angular momentum integrals  |
| dipole            | \-  | Dipole moment integrals     |
| higherMoment      | \-  | Octupole moment integrals   |
| quadrupole        | \-  | Quadrupole moment integrals |
| velocity          | \-  | Velocity integrals          |

**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_1elprop.json.conf 
:language: json
```

(sec:orca_2json.availableinformation.electronintegrals.1elrel.detailed)=
#### 1-electron relativistic integrals

1-electron relativistic integrals follow the same notation with the corresponding non-relativistic ones (see {ref}`sec:orca_2json.availableinformation.electronintegrals.1el.detailed`).

**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_1elrel.json.conf
:language: json
```
This example will produce and store in the corresponding json file the relativistic H-Matrix and S-Matrix.

(sec:orca_2json.availableinformation.electronintegrals.1elrelprop.detailed)=
#### 1-electron relativistic property integrals
Also available are 1-electron relativistic property integrals similar to the non-relativistic ones but with reduced options.

|                         |     |                                               |
|:------------------------|:---:|:----------------------------------------------|
| 1elPropertyRelIntegrals | \-  | relativistic 1-electron property integrals    |

Currently the following options are valid:
|                   |     |                             |
|:------------------|:---:|:----------------------------|
| dipole            | \-  | Dipole moment integrals     |
| quadrupole        | \-  | Quadrupole moment integrals |


**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_1elproprel.json.conf 
:language: json
```
(sec:orca_2json.availableinformation.electronintegrals.1elorigin.detailed)=
#### Origin setting
The origin of the electric property is per default the Cartesian origin but also the center of mass and the center of nuclear charges can be selected. Additionally an arbitrary position can be given as x,y,z coordinates when ori_el = 3 is chosen.
Currently the following options are valid:
|                   |     |                               |
|:------------------|:---:|:------------------------------|
| ori_el            | \-  | 0 - Cartesian origin          |
|                   | \-  | 1 - center of mass            |
|                   | \-  | 2 - center of nuclear charge  |
|                   | \-  | 3 - arbitrary position        |
| ori_el_xyz        | \-  | position of the origin(x,y,z) |


**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_origin.json.conf 
:language: json
```


(sec:orca_2json.availableinformation.electronintegrals.2el.detailed)=
#### 2-electron integrals

ORCA_2json can produce and write on disk three main categories of 2-electron integrals. 
1. Two-electron integrals in atomic basis ({ref}`sec:orca_2json.availableinformation.electronintegrals.2el.AO.detailed`)
2. Two-electron integrals in molecular basis ({ref}`sec:orca_2json.availableinformation.electronintegrals.2el.MO.detailed`)
3. Two-electron integrals using the resolution of identity approximation (RI). ({ref}`sec:orca_2json.availableinformation.electronintegrals.2el.RI.detailed`)
  
(sec:orca_2json.availableinformation.electronintegrals.2el.AO.detailed)=
##### 2-electron integrals in AO basis

In atomic basis the two-electron integrals can be saved in Coulomb order or in Exchange order. The keywords for the two options are shown in the next table. 

|             |     |                                               |
|:------------|:---:|:----------------------------------------------|
|  AO_PQRS    | \-  |  AO basis integrals in Coulomb order          | 
|  AO_PRQS    | \-  |  AO basis integrals in Exchange order         |


(sec:orca_2json.availableinformation.electronintegrals.2el.MO.detailed)=
##### 2-electron integrals in MO basis

In molecular basis, ORCA follows the accepted notation where by I,J,K and L we specify "internal" orbitals, meaning occupied in the reference wavefunction while by A,B,C and D we specify "external" orbitals, meaning orbitals that are empty in the reference wavefunction.
With P,Q,R and S we specify all possible orbitals, meaning both "internal" and "external".
The available keywords in *"orca_2json"* for the two-electron integrals in the molecular basis are the ones shown in the table below:
|             |     |                                               |
|:------------|:---:|:----------------------------------------------|
|  MO_IJKL    | \-  |  Coulomb 0-external                           |
|  MO_IJKA    | \-  |  Coulomb 1-external                           |
|  MO_IJAB    | \-  |  Coulomb 2-external                           |
|  MO_IABC    | \-  |  Coulomb 3-external                           |
|  MO_ABCD    | \-  |  Coulomb 4-external                           |
|  MO_PQRS    | \-  |  Coulomb ALL integrals                        |
|  MO_IKJL    | \-  |  Exchange 0-external                          |
|  MO_IKJA    | \-  |  Exchange 1-external                          |
|  MO_IAJB    | \-  |  Exchange 2-external                          |
|  MO_IBAC    | \-  |  Exchange 3-external                          |
|  MO_ACBD    | \-  |  Exchange 4-external                          |
|  MO_PRQS    | \-  |  Exchange ALL integrals                       |


**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_2elMO.json.conf 
:language: json
```

(sec:orca_2json.availableinformation.electronintegrals.2el.RI.detailed)=
##### RI 2-electron integrals

Using the Resolution of Identity (RI) one can create the integrals in a more efficient way. There are two main categories of RI integrals: the 3-index integrals, where only only half of the transformation has taken place, and the 4-index integrlas where the integrals are totally transformed in the molecular basis. The notation of the integrals follows the one we just described for the two-electron integrals in the molecular basis.

|             |     |                                               |
|:------------|:---:|:----------------------------------------------|
|  RI_IJV     | \-  |  RI 3-index 0-external                        |
|  RI_IAV     | \-  |  RI 3-index 1-external                        |
|  RI_ABV     | \-  |  RI 3-index 2-external                        |
|  RI_IJKL    | \-  |  RI 4-index Coulomb 0-external                |
|  RI_IJKA    | \-  |  RI 4-index Coulomb 1-external                |
|  RI_IJAB    | \-  |  RI 4-index Coulomb 2-external                |
|  RI_IABC    | \-  |  RI 4-index Coulomb 3-external                |
|  RI_ABCD    | \-  |  RI 4-index Coulomb 4-external                |
|  RI_IKJL    | \-  |  RI 4-index Exchange 0-external               |
|  RI_IKJA    | \-  |  RI 4-index Exchange 1-external               |
|  RI_IAJB    | \-  |  RI 4-index Exchange 2-external               |
|  RI_IBAC    | \-  |  RI 4-index Exchange 3-external               |
|  RI_ACBD    | \-  |  RI 4-index Exchange 4-external               |

In addition to the integrals in case of RI integrals also the used RI Metric is available through the option Vaux:
|             |     |                                               |
|:------------|:---:|:----------------------------------------------|
|  Vaux       | \-  |  true/false                                   |


**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_2elRI.json.conf 
:language: json
```

(sec:orca_2json.availableinformation.electronintegrals.2el.fulltrafo.detailed)=
##### Full Integral Transformation

The full transformation integrals can be selected via the FullTrafo keyword.

|              |     |                     |
|:-------------|:---:|:--------------------|
|  FullTrafo   | \-  | true/false          |


(sec:orca_2json.availableinformation.electronintegrals.2el.nonred.detailed)=
##### More 2-electron integrals

Also the non-redundant 2-electron integrals are available for the RI and the nonRI case. Therefore the options 2elNonRedIntegrals or 2elNonRedRIIntegrals must be specified.

|                       |     |                     |
|:----------------------|:---:|:--------------------|
|  2elNonRedIntegrals   | \-  | true/false          |
|  2elNonRedRIIntegrals | \-  | true/false          |


(sec:orca_2json.availableinformation.electronintegrals.2el.orbwin.detailed)=
##### Orbital Windows, AuxBasisType and Threshold

The orbital window can either be selected automatically by the transformation routine or given by the user via the OrbWin keyword. The internal and external space (i0,i1,a0,a1) is defined via an integer list.

for example:
|              |     |                          |
|:-------------|:---:|:-------------------------|
|  OrbWin      | \-  |  [0,8,9,15]              |
|              | \-  |  [0,12,13,30,0,12,13,30] |


The default AUX basis type is AuxC but can be changed with the keyword AuxBasisType. Please keep in mind that only those basis types used during the ORCA run can be selected.

|               |     |         |
|:--------------|:---:|:--------|
|  AuxBasisType | \-  |  AuxJ   |
|               | \-  |  AuxJK  |
|               | \-  |  AuxC   |

To reduce the number of the integrals the keyword Thresh can be used to decrease the selected integrals to save disk space. This effects ONLY the printing and not the accuracy of the generated integrals. The default print threshold is 1.e-15.

|              |     |                                      |
|:-------------|:---:|:-------------------------------------|
|  Thresh      | \-  |  printout threshold (default 1.e-15) |


**Example**
```{literalinclude} ../../orca_working_input/orca2json_example_integrals_2elRI_orbwin.json.conf 
:language: json
```


(sec:orca_2json.availableinformation.CIS.detailed)=
### TDDFT amplitude data (CIS/RPA) 

The TDDFT amplitudes and root informations can be requested (no triplet information yet).
The available options are:

|              |     |                              |
|:-------------|:---:|:-----------------------------|
| CIS          | \-  | TDDFT amplitudes             |
| CISNRoots    | \-  | informations of all  roots   |
| CISROOT      | \-  | List of roots                |


**Example**
```{literalinclude} ../../orca_working_input/orca_2json_examples/orca2json_example_CIS.json.conf 
:language: json
```
**Example**
```{literalinclude} ../../orca_working_input/orca_2json_examples/orca2json_example_CIS1.json.conf 
:language: json
```

(sec:orca_2json.availableinformation.jsonformat.detailed)=
### JSON Format

The JSON Format that is created can also be defined in the configuration file through the *JSONFormats* variable.

|                 |     |                                               |
|:----------------|:---:|:----------------------------------------------|
| JSONFormats     | \-  | JSON output format                            |

The available JSON formats are:

|         |     |                              |
|:--------|:---:|:-----------------------------|
| json    | \-  | ASCII JSON format            |
| bson    | \-  | binary JSON format           |
| ubjson  | \-  | universal binary json format |
| msgpck  | \-  | MessagePack format           |

**Example**
```{literalinclude} ../../orca_working_input/orca_2json_examples/orca2json_example_jsonFormats.json.conf 
:language: json
```

(sec:orca_2json.availableinformation.mocoefficients.detailed)=
### MO Coefficients 

|                 |     |                                               |
|:----------------|:---:|:----------------------------------------------|
| MOCoefficients  | \-  | molecule orbital information (true\false)     |

**Example**
```{literalinclude} ../../orca_working_input/orca_2json_examples/orca2json_example_MOCoefficients.json.conf 
:language: json
```

(sec:orca_2json.importData.detailed)=
## Import JSON data into ORCA

Some information like geometry, basis set and Molecular orbitals stored in the json
format written by orca_2json can be used to create a new gbw-file to be
specified as an orbital file in the orca input. The program is called as
a standalone via command line.

```orca
orca_2json BaseName.json -gbw
```

(sec:orca_2json.importInformationi.basic.detailed)=
### Basic Information

In order to create a functional gbw-file the following information must
be provided:

|              |     |                   |
|:-------------|:---:|:------------------|
| per Atom     | \-  | Basis             |
|              | \-  | Coords            |
|              | \-  | ElementNumber     |
|              | \-  | NuclearCharge     |
| per molecule | \-  | Charge            |
|              | \-  | CoordinateUnits   |
|              | \-  | HFTyp             |
|              | \-  | Multiplicity      |
|              | \-  | MolecularOrbitals |

**NOTE** 

Please keep in mind that *MolecularOrbitals* is a composite of 2 different components, namely *"EnergyUnit"* and *"MOs"*. Then *"MOs"* contains *"MOCoefficients"*, *"Occupancy"*, *"OrbitalEnergy"*, "OrbitalSymLabel"* and "OrbitalSymmetry"*. 


**Example**

The following file (let's call it *filename.json*) is a json file for H_2 molecule with STO-3G basis set. 

```{literalinclude} ../../orca_working_input/orca2json_example_basicInfo_in.json
:language: json
```

running then the command

```orca
orca_2json filename.json -gbw
```

should create a gbw file that ORCA can read.


### Definition of the real solid harmonic Gaussian orbitals

When integrals over real solid harmonic Gaussian orbitals are issued
into a JSON file, the precise definition of these orbitals becomes
important. ORCA uses its own peculiar conventions for the arrangement
and the phases of the individual components of the orbital shells.

The definitions of the angular parts of all shell components up to
angular momentum $k$ ($\ell=8$) are documented below. These correspond
to the real solid harmonics $S_{\ell,m}(x,y,z)$ that are normalized the
following way:

$$\begin{aligned}
    \int_0^\pi\sin\vartheta\,d\vartheta
    \int_{-\pi}^\pi d\varphi\,r^{-2\ell}S^2_{\ell,m}(x,y,z)=1
  \end{aligned}$$

$R_\ell(r)$ is the common radial part of a shell with angular momentum
$\ell$ that consists of a specific basis set dependent linear
combination of Gaussian primitives. It is normalized independently:

$$\begin{aligned}
    \int_0^\infty r^{2\ell+2}R^2_\ell(r)\,dr=1
  \end{aligned}$$

The factors $r^2$ and $\sin\vartheta$ in these integrals arise from the
volume element in spherical polar coordinates:

$$\begin{aligned}
    x &= r\sin\vartheta\cos\varphi \\
    y &= r\sin\vartheta\sin\varphi \\
    z &= r\cos\vartheta \\
    dx\,dy\,dz &= r^2\sin\vartheta\,dr\,d\vartheta\,d\varphi
  \end{aligned}$$

#### Angular momentum $\boldsymbol{s}$ ($\boldsymbol{\ell=0}$) 

$$\begin{aligned}
    s &= \frac{1}{2\sqrt{\pi} }\,R_s(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{p}$ ($\boldsymbol{\ell=1}$) 

$$\begin{aligned}
    N_p &= \frac{1}{2}\sqrt{\frac{3}{\pi} } \\
    p^{(0) } = p_{0\phantom{+} } &= N_p\,z\,R_p(r) \\
    p^{(1) } = p_{+1} &= N_p\,x\,R_p(r) \\
    p^{(2) } = p_{-1} &= N_p\,y\,R_p(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{d}$ ($\boldsymbol{\ell=2}$) 

$$\begin{aligned}
    N_d &= \frac{1}{2}\sqrt{\frac{15}{\pi} } \\
    d^{(0) } = d_{0\phantom{+} } &= \frac{\sqrt{3} }{6}\,N_d\left(3z^2-r^2\right)R_d(r) \\
    d^{(1) } = d_{+1} &= N_d\,xz\,R_d(r) \\
    d^{(2) } = d_{-1} &= N_d\,yz\,R_d(r) \\
    d^{(3) } = d_{+2} &= N_d\,\frac{x^2-y^2}{2}\,R_d(r) \\
    d^{(4) } = d_{-2} &= N_d\,xy\,R_d(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{f}$ ($\boldsymbol{\ell=3}$) 

$$\begin{aligned}
    N_f &= \frac{1}{2}\sqrt{\frac{105}{\pi} } \\
    f^{(0) } = f_{0\phantom{+} } &= \frac{\sqrt{15} }{30}\,N_f\,z\left(5z^2-3r^2\right)R_f(r) \\
    f^{(1) } = f_{+1} &= \frac{\sqrt{10} }{20}\,N_f\,x\left(5z^2-r^2\right)R_f(r) \\
    f^{(2) } = f_{-1} &= \frac{\sqrt{10} }{20}\,N_f\,y\left(5z^2-r^2\right)R_f(r) \\
    f^{(3) } = f_{+2} &= \frac{1}{2}\,N_f\left(x^2-y^2\right)z\,R_f(r) \\
    f^{(4) } = f_{-2} &= N_f\,xyz\,R_f(r) \\
    f^{(5) } = f_{+3} &= -\frac{\sqrt{6} }{12}\,N_f\,x\left(x^2-3y^2\right)R_f(r) \\
    f^{(6) } = f_{-3} &= -\frac{\sqrt{6} }{12}\,N_f\,y\left(3x^2-y^2\right)R_f(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{g}$ ($\boldsymbol{\ell=4}$) 

$$\begin{aligned}
    N_g &= \frac{3}{2}\sqrt{\frac{35}{\pi} } \\
    g^{(0) } = g_{0\phantom{+} } &= \frac{\sqrt{35} }{280}\,N_g\left(35z^4-30z^2r^2+3r^4\right)R_g(r) \\
    g^{(1) } = g_{+1} &= \frac{\sqrt{14} }{28}\,N_g\,xz\left(7z^2-3r^2\right)R_g(r) \\
    g^{(2) } = g_{-1} &= \frac{\sqrt{14} }{28}\,N_g\,yz\left(7z^2-3r^2\right)R_g(r) \\
    g^{(3) } = g_{+2} &= \frac{\sqrt{7} }{28}\,N_g\left(x^2-y^2\right)\left(7z^2-r^2\right)R_g(r) \\
    g^{(4) } = g_{-2} &= \frac{\sqrt{7} }{14}\,N_g\,xy\left(7z^2-r^2\right)R_g(r) \\
    g^{(5) } = g_{+3} &= -\frac{\sqrt{2} }{4}\,N_g\,x\left(x^2-3y^2\right)z\,R_g(r) \\
    g^{(6) } = g_{-3} &= -\frac{\sqrt{2} }{4}\,N_g\,y\left(3x^2-y^2\right)z\,R_g(r) \\
    g^{(7) } = g_{+4} &= -\frac{1}{8}\,N_g\left(x^4-6x^2y^2+y^4\right)R_g(r) \\
    g^{(8) } = g_{-4} &= -\frac{1}{2}\,N_g\,xy\left(x^2-y^2\right)R_g(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{h}$ ($\boldsymbol{\ell=5}$) 

$$\begin{aligned}
    N_h &= \frac{1}{2}\sqrt{\frac{11}{\pi} } \\
    h^{(0) }  = h_{0\phantom{+} } &= \frac{1}{8}\,N_h\,z\left(63z^4-70z^2r^2+15r^4\right)R_h(r) \\
    h^{(1) }  = h_{+1} &= \frac{\sqrt{15} }{8}\,N_h\,x\left(21z^4-14z^2r^2+r^4\right)R_h(r) \\
    h^{(2) }  = h_{-1} &= \frac{\sqrt{15} }{8}\,N_h\,y\left(21z^4-14z^2r^2+r^4\right)R_h(r) \\
    h^{(3) }  = h_{+2} &= \frac{\sqrt{105} }{4}\,N_h\left(x^2-y^2\right)z\left(3z^2-r^2\right)R_h(r) \\
    h^{(4) }  = h_{-2} &= \frac{\sqrt{105} }{2}\,N_h\,xyz\left(3z^2-r^2\right)R_h(r) \\
    h^{(5) }  = h_{+3} &= -\frac{\sqrt{70} }{16}\,N_h\,x\left(x^2-3y^2\right)\left(9z^2-r^2\right)R_h(r) \\
    h^{(6) }  = h_{-3} &= -\frac{\sqrt{70} }{16}\,N_h\,y\left(3x^2-y^2\right)\left(9z^2-r^2\right)R_h(r) \\
    h^{(7) }  = h_{+4} &= -\frac{3\sqrt{35} }{8}\,N_h\left(x^4-6x^2y^2+y^4\right)z\,R_h(r) \\
    h^{(8) }  = h_{-4} &= -\frac{3\sqrt{35} }{2}\,N_h\,xy\left(x^2-y^2\right)z\,R_h(r) \\
    h^{(9) }  = h_{+5} &= \frac{3\sqrt{14} }{16}\,N_h\,x\left(x^4-10x^2y^2+5y^4\right)R_h(r) \\
    h^{(10) } = h_{-5} &= \frac{3\sqrt{14} }{16}\,N_h\,y\left(5x^4-10x^2y^2+y^4\right)R_h(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{i}$ ($\boldsymbol{\ell=6}$) 

$$\begin{aligned}
    N_i &= \frac{1}{2}\sqrt{\frac{13}{\pi} } \\
    i^{(0) }  = i_{0\phantom{+} } &= \frac{1}{16}\,N_i\left(231z^6-315z^4r^2+105z^2r^4-5r^6\right)R_i(r) \\
    i^{(1) }  = i_{+1} &= \frac{\sqrt{21} }{8}\,N_i\,xz\left(33z^4-30z^2r^2+5r^4\right)R_i(r) \\
    i^{(2) }  = i_{-1} &= \frac{\sqrt{21} }{8}\,N_i\,yz\left(33z^4-30z^2r^2+5r^4\right)R_i(r) \\
    i^{(3) }  = i_{+2} &= \frac{\sqrt{210} }{32}\,N_i\left(x^2-y^2\right)\left(33z^4-18z^2r^2+r^4\right)R_i(r) \\
    i^{(4) }  = i_{-2} &= \frac{\sqrt{210} }{16}\,N_i\,xy\left(33z^4-18z^2r^2+r^4\right)R_i(r) \\
    i^{(5) }  = i_{+3} &= -\frac{\sqrt{210} }{16}\,N_i\,x\left(x^2-3y^2\right)z\left(11z^2-3r^2\right)R_i(r) \\
    i^{(6) }  = i_{-3} &= -\frac{\sqrt{210} }{16}\,N_i\,y\left(3x^2-y^2\right)z\left(11z^2-3r^2\right)R_i(r) \\
    i^{(7) }  = i_{+4} &= -\frac{3\sqrt{7} }{16}\,N_i\left(x^4-6x^2y^2+y^4\right)\left(11z^2-r^2\right)R_i(r) \\
    i^{(8) }  = i_{-4} &= -\frac{3\sqrt{7} }{4}\,N_i\,xy\left(x^2-y^2\right)\left(11z^2-r^2\right)R_i(r) \\
    i^{(9) }  = i_{+5} &= \frac{3\sqrt{154} }{16}\,N_i\,x\left(x^4-10x^2y^2+5y^4\right)z\,R_i(r) \\
    i^{(10) } = i_{-5} &= \frac{3\sqrt{154} }{16}\,N_i\,y\left(5x^4-10x^2y^2+y^4\right)z\,R_i(r) \\
    i^{(11) } = i_{+6} &= \frac{\sqrt{462} }{32}\,N_i\left(x^2-y^2\right)\left(x^4-14x^2y^2+y^4\right)R_i(r) \\
    i^{(12) } = i_{-6} &= \frac{\sqrt{462} }{16}\,N_i\,xy\left(3x^4-10x^2y^2+3y^4\right)R_i(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{j}$ ($\boldsymbol{\ell=7}$) 

$$\begin{aligned}
    N_j &= \frac{1}{2}\sqrt{\frac{15}{\pi} } \\
    j^{(0) }  = j_{0\phantom{+} } &= \frac{1}{16}\,N_j\,z\left(429z^6-693z^4r^2+315z^2r^4-35r^6\right)R_j(r) \\
    j^{(1) }  = j_{+1} &= \frac{\sqrt{7} }{32}\,N_j\,x\left(429z^6-495z^4r^2+135z^2r^4-5r^6\right)R_j(r) \\
    j^{(2) }  = j_{-1} &= \frac{\sqrt{7} }{32}\,N_j\,y\left(429z^6-495z^4r^2+135z^2r^4-5r^6\right)R_j(r) \\
    j^{(3) }  = j_{+2} &= \frac{\sqrt{42} }{32}\,N_j\left(x^2-y^2\right)z\left(143z^4-110z^2r^2+15r^4\right)R_j(r) \\
    j^{(4) }  = j_{-2} &= \frac{\sqrt{42} }{16}\,N_j\,xyz\left(143z^4-110z^2r^2+15r^4\right)R_j(r) \\
    j^{(5) }  = j_{+3} &= -\frac{\sqrt{21} }{32}\,N_j\,x\left(x^2-3y^2\right)\left(143z^4-66z^2r^2+3r^4\right)R_j(r) \\
    j^{(6) }  = j_{-3} &= -\frac{\sqrt{21} }{32}\,N_j\,y\left(3x^2-y^2\right)\left(143z^4-66z^2r^2+3r^4\right)R_j(r) \\
    j^{(7) }  = j_{+4} &= -\frac{\sqrt{231} }{16}\,N_j\left(x^4-6x^2y^2+y^4\right)z\left(13z^2-3r^2\right)R_j(r) \\
    j^{(8) }  = j_{-4} &= -\frac{\sqrt{231} }{4}\,N_j\,xy\left(x^2-y^2\right)z\left(13z^2-3r^2\right)R_j(r) \\
    j^{(9) }  = j_{+5} &= \frac{\sqrt{231} }{32}\,N_j\,x\left(x^4-10x^2y^2+5y^4\right)\left(13z^2-r^2\right)R_j(r) \\
    j^{(10) } = j_{-5} &= \frac{\sqrt{231} }{32}\,N_j\,y\left(5x^4-10x^2y^2+y^4\right)\left(13z^2-r^2\right)R_j(r) \\
    j^{(11) } = j_{+6} &= \frac{\sqrt{6006} }{32}\,N_j\left(x^2-y^2\right)\left(x^4-14x^2y^2+y^4\right)z\,R_j(r) \\
    j^{(12) } = j_{-6} &= \frac{\sqrt{6006} }{16}\,N_j\,xy\left(3x^2-y^2\right)\left(x^2-3y^2\right)z\,R_j(r) \\
    j^{(13) } = j_{+7} &= -\frac{\sqrt{429} }{32}\,N_j\,x\left(x^6-21x^4y^2+35x^2y^4-7y^6\right)R_j(r) \\
    j^{(14) } = j_{-7} &= -\frac{\sqrt{429} }{32}\,N_j\,y\left(7x^6-35x^4y^2+21x^2y^4-y^6\right)R_j(r)
  \end{aligned}$$

#### Angular momentum $\boldsymbol{k}$ ($\boldsymbol{\ell=8}$) 

$$\begin{aligned}
    N_k &= \frac{1}{2}\sqrt{\frac{17}{\pi} } \\
    k^{(0) }  = k_{0\phantom{+} } &= \frac{1}{128}\,N_k\left(6435z^8-12012z^6r^2+6930z^4r^4-1260z^2r^6+35r^8\right)R_k(r) \\
    k^{(1) }  = k_{+1} &= \frac{3}{32}\,N_k\,xz\left(715z^6-1001z^4r^2+385z^2r^4-35r^6\right)R_k(r) \\
    k^{(2) }  = k_{-1} &= \frac{3}{32}\,N_k\,yz\left(715z^6-1001z^4r^2+385z^2r^4-35r^6\right)R_k(r) \\
    k^{(3) }  = k_{+2} &= \frac{3\sqrt{70} }{64}\,N_k\left(x^2-y^2\right)\left(143z^6-143z^4r^2+33z^2r^4-r^6\right)R_k(r) \\
    k^{(4) }  = k_{-2} &= \frac{3\sqrt{70} }{32}\,N_k\,xy\left(143z^6-143z^4r^2+33z^2r^4-r^6\right)R_k(r) \\
    k^{(5) }  = k_{+3} &= -\frac{\sqrt{1155} }{32}\,N_k\,x\left(x^2-3y^2\right)z\left(39z^4-26z^2r^2+3r^4\right)R_k(r) \\
    k^{(6) }  = k_{-3} &= -\frac{\sqrt{1155} }{32}\,N_k\,y\left(3x^2-y^2\right)z\left(39z^4-26z^2r^2+3r^4\right)R_k(r) \\
    k^{(7) }  = k_{+4} &= -\frac{3\sqrt{77} }{64}\,N_k\,\left(x^4-6x^2y^2+y^4\right)\left(65z^4-26z^2r^2+r^4\right)R_k(r) \\
    k^{(8) }  = k_{-4} &= -\frac{3\sqrt{77} }{16}\,N_k\,xy\left(x^2-y^2\right)\left(65z^4-26z^2r^2+r^4\right)R_k(r) \\
    k^{(9) }  = k_{+5} &= \frac{3\sqrt{1001} }{32}\,N_k\,x\left(x^4-10x^2y^2+5y^4\right)z\left(5z^2-r^2\right)R_k(r) \\
    k^{(10) } = k_{-5} &= \frac{3\sqrt{1001} }{32}\,N_k\,y\left(5x^4-10x^2y^2+y^4\right)z\left(5z^2-r^2\right)R_k(r) \\
    k^{(11) } = k_{+6} &= \frac{\sqrt{858} }{64}\,N_k\left(x^2-y^2\right)\left(x^4-14x^2y^2+y^4\right)\left(15z^2-r^2\right)R_k(r) \\
    k^{(12) } = k_{-6} &= \frac{\sqrt{858} }{32}\,N_k\,xy\left(x^2-3y^2\right)\left(3x^2-y^2\right)\left(15z^2-r^2\right)R_k(r) \\
    k^{(13) } = k_{+7} &= -\frac{3\sqrt{715} }{32}\,N_k\,x\left(x^6-21x^4y^2+35x^2y^4-7y^6\right)z\,R_k(r) \\
    k^{(14) } = k_{-7} &= -\frac{3\sqrt{715} }{32}\,N_k\,y\left(7x^6-35x^4y^2+21x^2y^4-y^6\right)z\,R_k(r) \\
    k^{(15) } = k_{+8} &= -\frac{3\sqrt{715} }{128}\,N_k\left(x^8-28x^6y^2+70x^4y^4-28x^2y^6+y^8\right)R_k(r) \\
    k^{(16) } = k_{-8} &= -\frac{3\sqrt{715} }{16}\,N_k\,xy\left(x^2-y^2\right)\left(x^4-6x^2y^2+y^4\right)R_k(r)
  \end{aligned}$$