source: ThirdParty/mpqc_open/src/bin/mpqc/validate/ref/opt_b2h6scf631gsd2hopt.in

% Emacs should use -*- KeyVal -*- mode
% this file was automatically generated
% label: optimization test series
% molecule specification
molecule<Molecule>: (
  symmetry = D2H
  unit = angstrom
  { atoms geometry } = {
     H     [     1.000000000000     0.000000000000     1.460000000000 ]
     H     [    -1.000000000000     0.000000000000     1.460000000000 ]
     B     [     0.000000000000     0.000000000000     0.900000000000 ]
     H     [     0.000000000000     0.940000000000     0.000000000000 ]
     H     [     0.000000000000    -0.940000000000     0.000000000000 ]
     B     [     0.000000000000     0.000000000000    -0.900000000000 ]
     H     [     1.000000000000     0.000000000000    -1.460000000000 ]
     H     [    -1.000000000000     0.000000000000    -1.460000000000 ]
  }
)
% basis set specification
basis<GaussianBasisSet>: (
  name = "6-31G*"
  molecule = $:molecule
)
mpqc: (
  checkpoint = no
  savestate = no
  restart = no
  % molecular coordinates for optimization
  coor<SymmMolecularCoor>: (
    molecule = $:molecule
    generator<IntCoorGen>: (
      molecule = $:molecule
    )
  )
  do_energy = yes
  do_gradient = no
  % method for computing the molecule's energy
  mole<CLHF>: (
    molecule = $:molecule
    basis = $:basis
    coor = $..:coor
    memory = 32000000
    total_charge = 0
    multiplicity = 1
    print_npa = yes
    guess_wavefunction<CLHF>: (
      molecule = $:molecule
      total_charge = 0
      multiplicity = 1
      basis<GaussianBasisSet>: (
        molecule = $:molecule
        name = "STO-3G"
      )
      memory = 32000000
    )
  )
  optimize = yes
  % optimizer object for the molecular geometry
  opt<QNewtonOpt>: (
    max_iterations = 20
    function = $..:mole
    update<BFGSUpdate>: ()
    convergence<MolEnergyConvergence>: (
      cartesian = yes
      energy = $..:..:mole
    )
  )
)