Задание базиса в Gaussian 03W

[Wan]

У меня стоит задача, рассчитать молекулу C1Fe3H6OH. выбрал метод B3LYP. Для атомов О и Н был выбран валентно-расщепленный базис 6-31G(d,p) для остовных электронов атома Fe псевдопотенциал LanL2MB.
Как задать такой расчет в гауссиане?? Интересует задание самого метода расчета. структуру задаю в виде z-матрицы

[Wolf1]

Как задать, чтобы одни атомы расчитывались в одном базисе, а другие в другом, сам хотел бы узнать.
А вот, чтобы задать расчет для всех атомов нужно вписать вот эти строки (сам работаю на G98)
# RB3LYP/LANL2MB (в G98 )
или
# RB3LYP/6-31G(d,p)

[Wolf1]

Немного поправлюсь:
# RB3LYP/LANL2MB Opt Pop=(Minimal,) Test (в G98 )
или
# RB3LYP/6-31G(d,p) Opt Pop=(Minimal,) Test

[Raiden]

Вот пример гауссиановского файла:

%chk=chkfile.chk
# rb3lyp/Gen pseudo=read opt optcyc=99 freq=noraman vshift=200 iop(5/13=1)

Комментарий

<заряд/мультиплетность>
<координаты>

Pb 0
Lanl2DZ
****
O 0
6-31G*
****
C 0
6-31G*
****
B 0
6-31G*
****
H 0
6-31G*
****

Pb 0
Lanl2DZ