QED calculation tutorials

QED-HF calculation

A QED-Hartree-Fock (HF) calculation is performed by specifying a QED-HF (qed-hf) wave function in the method section of the input file. To start a ground state QED-HF we specify

- method
  qed-hf

- do
  ground state

A minimal working example for a QED-HF single-point calculation on water with wavevector \(\vec{k}=(0,0,1)\), frequency \(\omega=0.5\) au and coupling \(\lambda=0.05\) au:

- system
  charge: 0

- do
  ground state

- method
  qed-hf

- boson
  interaction type: photon
  modes: 2
  frequency: {0.5, 0.5}
  wavevector: {0.0, 0.0, 1.0}
  coupling: {0.05, 0.05}

- geometry
  basis: cc-pVDZ
  H          0.86681        0.60144        0.00000
  H         -0.86681        0.60144        0.00000
  O          0.00000       -0.07579        0.00000

Save this as h2o.inp and invoke the launch script.

path/to/eT_launch.py h2o.inp

After the calculation finished you should find h2o.out and h2o.timing.out in your working directory. If the calculation exited successfully (look for eT terminated successfully! at the bottom of the file), a summary of the calculation should be printed at the end.

SC-QED-HF calculation

The wave function has a coherent state parameter per each orbital in the basis that diagonalizes the dipole operator. This way electron photon correlation is naturally included in the methodology. A Strong Coupling QED-Hartree-Fock (SC-QED-HF) calculation is performed by specifying a SC-QED-HF (sc-qed-hf) wave function in the method section of the input file. To start a ground state SC-QED-HF we specify

- method
  sc-qed-hf

- do
  ground state

A minimal working example for a SC-QED-HF single-point calculation on water with polarization \(\vec{\epsilon}=(0,0,1)\), frequency \(\omega=0.5\) au and coupling \(\lambda=0.05\) au:

- system
  charge: 0

- do
  ground state

- method
  sc-qed-hf

- boson
  interaction type: photon
  modes: 1
  frequency: {0.5}
  polarization: {0.0, 0.0, 1.0}
  coupling: {0.05}

- geometry
  basis: cc-pVDZ
  H          0.86681        0.60144        0.00000
  H         -0.86681        0.60144        0.00000
  O          0.00000       -0.07579        0.00000

Save this as h2o.inp and invoke the launch script.

path/to/eT_launch.py h2o.inp

After the calculation finished you should find h2o.out and h2o.timing.out in your working directory. If the calculation exited successfully (look for eT terminated successfully! at the bottom of the file), a summary of the calculation should be printed at the end.

QED-FCI calculation

A QED Full Configuration Interaction (FCI) calculation is performed by specifying a qed-hf/hf/sc-qed-hf (for closed-shell systems) or a rohf/uhf wave function (for open-shell systems) and the QED-FCI wave function (qed-fci) in the method section of the input file. E.g., for a QED-HF calculation we write

- method
  qed-hf
  qed-fci

To perform a single-point QED-FCI calculation it is necessary to specify ground state do section:

- do
  ground state

As in standard FCI the calculations require the solver ci section in which it is necessary to specify the number of states you want to compute in your calculation, for instance:

- solver fci
  states: 5

A minimal working example for a QED-FCI single-point calculation on water with polarization \(\vec{\epsilon}=(0,0,1)\), frequency \(\omega=0.5\) au and coupling \(\lambda=0.05\) au:

- do
  ground state

- memory
  available: 16

- method
  hf
  fci

- boson
  interaction type:         photon
  modes:             1
  boson states:     {1}
  coupling:         {0.05}
  frequency:        {0.5}
  polarization:     {0, 0, 1}

- solver fci
  states: 5

- geometry
  basis: sto-3g
  H          0.86681        0.60144        0.00000
  H         -0.86681        0.60144        0.00000
  O          0.00000       -0.07579        0.00000

Note

For open-shell calculations you might want to specify the multiplicity in the system section. The default multiplicity is 1.

Save this as h2o.inp and invoke the launch script.

path/to/eT_launch.py h2o.inp

After the calculation finished you should find h2o.out and h2o.timing.out in your working directory. If the calculation exited successfully (look for eT terminated successfully! at the bottom of the file), a summary of the calculation should be printed at the end.

QED-FCI calculation

It is also possible to perform QED-FCI calculations only in an active space of canonical molecular orbitals, leading to a Complete Active Space QED Configuration Interaction (QED-CASCI) calculation.

To perform a QED-CASCI calculation the keyword qed-fci needs to be exchanged with qed-casci in the method section. Similartly to QED-FCI, the QED-CASCI calculations can be computed using the QED-HF orbitals or the SC-QED-HF orbitals. This is specified by the qed-hf or sc-qed-hf wave function in the method section of the input file. E.g.,

- method
  sc-qed-hf
  qed-casci

Additionally, it is necessary to specify the number of active electrons and active canonical orbitals included in the active space. This is done using the canonical keyword in the active space section. For a CASCI calculation with 2 active electrons in 4 active canonical orbitals we write

- active space
   canonical: {2, 4}

Warning

For a charged system, energy invariance under translation is not guaranteed if the QED-HF orbitals are employed. The use of the SC-QED-HF wave function is thus advised.

Relevant input sections

Solver SCF

HF mean value

CI solver section

boson

Active space section