12.07.2015 Views

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

C. Programas 117call wgt(qa(j,i),temp(i),enra(j))enddoC Atomselse if (nat_ra.eq.1) thenqa(j,i)=1.D0/ansyma(j)if (iop6.eq.1) thenqeleca=qelect(numsta,<strong>de</strong>ga,difena,temp(i))qa(j,i)=qa(j,i)*qelecaendifendifCC Bimolecu<strong>la</strong>r reactionCC Reactant Bif (iop1.eq.2) thendo j=1,n_rbC Moleculeif (nat_rb.gt.1) thenCC Rotational Partition FunctionCC REACT-Bif (linb(j).eq.0) thenC 1: Linear Moleculeqrotb(j,i)=zrotlin(dmirb(j,1),temp(i))nfreq=3*nat_rb-5elseC 2: Non-Linear Moleculeqrotb(j,i)=zrotnolin(dmirb(j,1),dmirb(j,2),dmirb(j,3),& temp(i))nfreq=3*nat_rb-6endifCCCCCVibrational Partition Functionqvibb(j,i)=1.D0do m=1,nfreqqvibbp(j,i,m)=qqvib(wb(j,m),j,temp(i))qvibb(j,i)=qvibb(j,i)*qvibbp(j,i,m)enddoAtomelse if (nat_rb.eq.1) thenqb(j,i)=1.D0qrotb(j,i)=1.D0qvibb(j,i)=1.D0endifTras<strong>la</strong>tional Partition Functionamatot=0.D0ambtot=0.D0do k=1,nat_raamatot=amatot+amasa(j,k)enddodo k=1,nat_rbambtot=ambtot+amasb(j,k)enddoqtrans=ztrans(amatot,ambtot,temp(i))qb(j,i)=qrotb(j,i)*qvibb(j,i)Electronic Partition Function

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

Saved successfully!

Ooh no, something went wrong!