140 C.2. CONFORATEC Full F matrixdo i=1,ndo j=1,iffcc(j,i)=ffcc(i,j)enddoenddoC copy copy copynumb=3*natl=0do i=1,numb/3do i2=1,3l=l+1k=0do j=1,numb/3do jj=1,3k=k+1ffcc(k,l)=ffcc(k,l)/(sqrt(amas(ll,i))*sqrt(amas(ll,j)))C write(6,*) l,i,j,ffcc(k,l),amas(ll,i),amas(ll,j)enddoenddoenddoenddoCwrite(6,*) ((ffcc(i,j),j=1,i),i=1,3*nat)k=0do i=1,numbdo j=1,ik=k+1a(k)=ffcc(i,j)*2240.75D0C write(6,*) a(k)enddoenddoC Diagonalize the FC matrixCALL HQRII(A,NUMB,NUMB,E,FFCC)C write(6,*) (e(i),i=1,3*nat)C write(6,*)C write(6,*)C 2 veces pitpi=6.283185307D0CONS=1.D3*2.0455D0/(2.99792458D0*TPI)do i=1,numbsige(i)=signo(e(i))e(i)=sqrt(abs(e(i)))*consenddoC do i=1,numbC write(6,*) e(i)C enddoC write(6,*)CC nat > 2 normal formatC nat = 2 extra formatif (nat.gt.2) ndumb=6if (nat.eq.2.or.lin.eq.0) ndumb=5Cif (vlts.eqv..true.) dtempor=e(1)do j=3*nat,ndumb+1,-1
C. Programas 141w(ll,j-ndumb)=e(j)sigw(j-ndumb)=sige(j)do i=1,natdLx(1,i,j-ndumb)=FFcc(3*i-2,j)dLx(2,i,j-ndumb)=FFcc(3*i-1,j)dLx(3,i,j-ndumb)=FFcc(3*i,j)enddoenddoCif (vlts.eqv..true.) thenC w(ll,ndumb+1)=-dtemporw(ll,1)=-dtemporendifC write(6,*) ’w’,(w(ll,j),j=1,9)CC l1=1C l2=3C33 do i=1,natC if(i.eq.1) thenCC write(6,5) l1,l2CC write(6,10) (sigw(j)*w(ll,j),j=l1,l2)CC write(6,*)C endifCC write(6,20) ((dLx(ic,i,j),ic=1,3),j=l1,l2)C enddoC l1=l2+1C l2=l2+3C if(l1.le.n) thenC if(l2.le.n) thenC goto 33C else if(l2.gt.n) thenC l2=nC goto 33C elseC continueC endifC endif5 format(//,1x,’w(’,i3,’) to w(’,i3,’):’,/)10 format(2x,3(10x,f8.1,5x))20 format(4x,3(2x,3(1x,f6.3)))C copy copy copyCreturnstopendC======================================================================SUBROUTINE HQRII(A,N,M,E,V)IMPLICIT DOUBLE PRECISION (A-H,O-Z)inclu<strong>de</strong> ’data.h’REAL*8 A((3*natoms)**2), E(3*natoms), V(3*natoms,3*natoms)C REAL*8 A((300)**2), E(300), V(300,300)*************************************************************** HQRII IS A DIAGONALISATION ROUTINE, WRITTEN BY YOSHITAKA BEPPU OF
- Page 5:
VANTONIO FEZNÁNDEZ RAMOS Y SAULO A
- Page 9 and 10:
Índice general1. Introducción y o
- Page 12 and 13:
XIIÍNDICE DE FIGURAS3.4. Confórme
- Page 15 and 16:
Capítulo 1Introducción y objetivo
- Page 17 and 18:
1. Introducción y objetivos 3su ef
- Page 19 and 20:
1. Introducción y objetivos 5Tabla
- Page 21 and 22:
1. Introducción y objetivos 7algun
- Page 23 and 24:
1. Introducción y objetivos 91.3.
- Page 25 and 26:
1. Introducción y objetivos 11Tabl
- Page 27 and 28:
1. Introducción y objetivos 13la p
- Page 29 and 30:
1. Introducción y objetivos 15C 8
- Page 31 and 32:
1. Introducción y objetivos 17Tabl
- Page 33 and 34:
1. Introducción y objetivos 19C 8
- Page 35 and 36:
1. Introducción y objetivos 21Figu
- Page 37 and 38:
1. Introducción y objetivos 23medi
- Page 39 and 40:
Capítulo 2Métodos teóricosPara a
- Page 41 and 42:
2. Métodos teóricos 27Sin embargo
- Page 43 and 44:
2. Métodos teóricos 29mínimo. Mi
- Page 45 and 46:
2. Métodos teóricos 31podrán omi
- Page 47 and 48:
2. Métodos teóricos 33Para poder
- Page 49 and 50:
2. Métodos teóricos 35esférica d
- Page 51 and 52:
2. Métodos teóricos 37recibe el n
- Page 53 and 54:
2. Métodos teóricos 39Las integra
- Page 55 and 56:
2. Métodos teóricos 41electrón l
- Page 57 and 58:
2. Métodos teóricos 43son de áto
- Page 59 and 60:
2. Métodos teóricos 45e INDO. Sin
- Page 61 and 62:
2. Métodos teóricos 47centros y c
- Page 63 and 64:
2. Métodos teóricos 492.3.3. Mét
- Page 65 and 66:
2. Métodos teóricos 51cierto tama
- Page 67 and 68:
2. Métodos teóricos 53Teorema var
- Page 69 and 70:
2. Métodos teóricos 55el funciona
- Page 71 and 72:
2. Métodos teóricos 57• Por úl
- Page 73 and 74:
2. Métodos teóricos 59este nivel
- Page 75 and 76:
2. Métodos teóricos 61transición
- Page 77 and 78:
2. Métodos teóricos 63generalizad
- Page 79 and 80:
2. Métodos teóricos 65expresión
- Page 81 and 82:
2. Métodos teóricos 67(o 2) rotac
- Page 83 and 84:
2. Métodos teóricos 69partición
- Page 85 and 86:
2. Métodos teóricos 71en un incre
- Page 87:
2. Métodos teóricos 73La segunda
- Page 90 and 91:
76 3.1. Superficie de energía pote
- Page 92 and 93:
78 3.1. Superficie de energía pote
- Page 94 and 95:
80 3.1. Superficie de energía pote
- Page 96 and 97:
82 3.1. Superficie de energía pote
- Page 98 and 99:
Tabla 3.2: Energías en kcal/mol y
- Page 100 and 101:
86 3.1. Superficie de energía pote
- Page 102 and 103:
88 3.1. Superficie de energía pote
- Page 104 and 105: Figura 3.9: Estructuras conformacio
- Page 106 and 107: Tabla 3.3: Energías y principales
- Page 108 and 109: 94 3.2. Cálculos de dinámica dire
- Page 110 and 111: 96 3.2. Cálculos de dinámica dire
- Page 112 and 113: 98 3.2. Cálculos de dinámica dire
- Page 114 and 115: 100 3.2. Cálculos de dinámica dir
- Page 117 and 118: Apéndice AMétodo simplex descende
- Page 119: A. Método simplex descendente mult
- Page 122 and 123: 108el que se encuentra en favor de
- Page 124 and 125: 110 C.1. MORATE2end programCsubrout
- Page 126 and 127: 112 C.1. MORATE2if (i.eq.icurt.or.i
- Page 128 and 129: 114 C.2. CONFORATE& ,wts(nconfs,nfr
- Page 130 and 131: 116 C.2. CONFORATEC TSif (nat_ts.gt
- Page 132 and 133: 118 C.2. CONFORATECCCCif (iop6.eq.1
- Page 134 and 135: 120 C.2. CONFORATEakp_vtst(n,i)=ct(
- Page 136 and 137: 122 C.2. CONFORATECCwrite(6,30) tem
- Page 138 and 139: 124 C.2. CONFORATEC Compute traslat
- Page 140 and 141: 126 C.2. CONFORATEC Read input data
- Page 142 and 143: 128 C.2. CONFORATECCCCenddoelsewrit
- Page 144 and 145: 130 C.2. CONFORATEstop1011 write(6,
- Page 146 and 147: 132 C.2. CONFORATElin=0C write(6,*)
- Page 148 and 149: 134 C.2. CONFORATE& ,numsta,numstb,
- Page 150 and 151: 136 C.2. CONFORATEwrite(6,33) (dmir
- Page 152 and 153: 138 C.2. CONFORATEendC Compute the
- Page 156 and 157: 142 C.2. CONFORATE* NAGOYA UNIVERSI
- Page 158 and 159: 144 C.2. CONFORATE* * * * * * * * *
- Page 160 and 161: 146 C.2. CONFORATEDO 400 J=1,N400 V
- Page 162 and 163: 148 C.2. CONFORATE298. 63.0843 # Co
- Page 165 and 166: Bibliografía[1] C. Funk, “On the
- Page 167 and 168: BIBLIOGRAFÍA 153[24] A. Verloop, A
- Page 169 and 170: BIBLIOGRAFÍA 155[47] D. Hartree,
- Page 171: BIBLIOGRAFÍA 157G. Liu, A. Liashen