19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

24 THE CLOSED SHELL CCSD PROGRAM 183<br />

24.2 Quadratic configuration interaction, QCI<br />

QCI or QCISD performs quadratic configuration interaction, QCISD. Using the QCI(T) or<br />

QCISD(T) commands, the contributions of connected triples are also computed by perturbation<br />

theory. Normally, no further input is needed if the QCI card follows the corresponding HF-SCF.<br />

Otherwise, occupancies and orbitals can be specified as in the CI program. For modifying DIIS<br />

directives, see section 24.5. For first-order QCISD and QCIDF(T) properties see section .<br />

For avoiding error exits in case of no convergence, see CCSD(T).<br />

24.3 Brueckner coupled-cluster calculations, BCCD<br />

BCCD,[SAVE=record],[PRINT],[TYPE=,type]<br />

BCCD performs a Brueckner coupled-cluster calculation and computes Brueckner orbitals. With<br />

these orbitals, the amplitudes of the singles vanish at convergence. Using the BCCD(T) command,<br />

the contributions of connected triples are also computed by perturbation theory. Normally,<br />

no further input is needed if the BCCD card follows the corresponding HF-SCF. Otherwise,<br />

occupancies and orbitals can be specified as in the CI program. BRUECKNER parameters<br />

can be modified using the BRUECKNER directive.<br />

The Brueckner orbitals and approximate density matrix can be saved on a <strong>MOLPRO</strong> dump record<br />

using the SAVE option. The orbitals are printed if the PRINT option is given. TYPE can be used<br />

to specify the type of the approximate density to be computed:<br />

TYPE=REF<br />

TYPE=TOT<br />

TYPE=ALL<br />

Compute and store density of reference determinant only (default).<br />

This corresponds to the BOX (Brueckner orbital expectation value)<br />

method of Chem. Phys. Lett. 315, 248 (1999).<br />

Compute and store density with contribution of pair amplitudes (linear<br />

terms). Normally, this does not seem to lead to an improvement.<br />

Compute and store both densities<br />

Note: The expectation variables are stored in variables as usual. In the case that both densities<br />

are made, the variables contain two values, the first corresponding to REF and the second to<br />

TOT (e.g., DMZ(1) and DMZ(2)). If TYPE=REF or TYPE=TOT is give, only the corresponding<br />

values are stored.<br />

For avoiding error exits in case of no convergence, see CCSD(T).<br />

24.3.1 The BRUECKNER directive<br />

BRUECKNER,orbbrk,ibrstr,ibrueck,brsfak;<br />

This directive allows the modification of options for Brueckner calculations. Normally, none of<br />

the options has to be specified, and the BCCD command can be used to perform a Brueckner<br />

CCD calculation.<br />

orbbrk:<br />

ibrstr:<br />

ibrueck:<br />

brsfak:<br />

if nonzero, the Brueckner orbitals are saved on this record.<br />

First iteration in which orbitals are modified (default=3).<br />

Iteration increment between orbital updates (default=1).<br />

Scaling factor for singles in orbital updates (default=1).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!