Entering Gaussian System, Link 0=g03 Initial command: /apps/gaussian/g09_d01/g09/l1.exe "/home/scan-user-1/run/79590/Gau-27924.inp" -scrdir="/home/scan-user-1/run/79590/" Entering Link 1 = /apps/gaussian/g09_d01/g09/l1.exe PID= 27925. Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2013, Gaussian, Inc. All Rights Reserved. This is part of the Gaussian(R) 09 program. It is based on the Gaussian(R) 03 system (copyright 2003, Gaussian, Inc.), the Gaussian(R) 98 system (copyright 1998, Gaussian, Inc.), the Gaussian(R) 94 system (copyright 1995, Gaussian, Inc.), the Gaussian 92(TM) system (copyright 1992, Gaussian, Inc.), the Gaussian 90(TM) system (copyright 1990, Gaussian, Inc.), the Gaussian 88(TM) system (copyright 1988, Gaussian, Inc.), the Gaussian 86(TM) system (copyright 1986, Carnegie Mellon University), and the Gaussian 82(TM) system (copyright 1983, Carnegie Mellon University). Gaussian is a federally registered trademark of Gaussian, Inc. This software contains proprietary and confidential information, including trade secrets, belonging to Gaussian, Inc. This software is provided under written license and may be used, copied, transmitted, or stored only in accord with that written license. The following legend is applicable only to US Government contracts under FAR: RESTRICTED RIGHTS LEGEND Use, reproduction and disclosure by the US Government is subject to restrictions as set forth in subparagraphs (a) and (c) of the Commercial Computer Software - Restricted Rights clause in FAR 52.227-19. Gaussian, Inc. 340 Quinnipiac St., Bldg. 40, Wallingford CT 06492 --------------------------------------------------------------- Warning -- This program may not be used in any manner that competes with the business of Gaussian, Inc. or will provide assistance to any competitor of Gaussian, Inc. The licensee of this program is prohibited from giving any competitor of Gaussian, Inc. access to this program. By using this program, the user acknowledges that Gaussian, Inc. is engaged in the business of creating and licensing software in the field of computational chemistry and represents and warrants to the licensee that it is not a competitor of Gaussian, Inc. and that it will not use this program in any manner prohibited above. --------------------------------------------------------------- Cite this work as: Gaussian 09, Revision D.01, M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, T. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, O. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, and D. J. Fox, Gaussian, Inc., Wallingford CT, 2013. ****************************************** Gaussian 09: ES64L-G09RevD.01 24-Apr-2013 9-Oct-2013 ****************************************** %nprocshared=8 Will use up to 8 processors via shared memory. %mem=13000MB %NoSave %Chk=chk.chk %rwf=/tmp/pbs.5398283.cx1b/rwf -------------------------------------------------------- # freq b3lyp/6-31g(d,p) nosymm pop=nbo geom=connectivity -------------------------------------------------------- 1/10=4,30=1,38=1,57=2/1,3; 2/12=2,15=1,17=6,18=5,40=1/2; 3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=2,74=-5,140=1/1,2,3; 4//1; 5/5=2,38=5,98=1/2; 8/6=4,10=90,11=11/1; 11/6=1,8=1,9=11,15=111,16=1,31=1/1,2,10; 10/6=1,31=1/2; 6/7=2,8=2,9=2,10=2,18=1,28=1,40=1/1,7; 7/8=1,10=1,25=1,30=1/1,2,3,16; 1/10=4,30=1/3; 99//99; -------------- NH3_FREQ_ARRON -------------- Symbolic Z-matrix: Charge = 0 Multiplicity = 1 N -0.11364 0.60656 -0.00001 H 0.28379 -0.33062 0.00003 H 0.28381 1.07511 0.81163 H 0.28376 1.07517 -0.81166 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Initialization pass. Trust Radius=3.00D-01 FncErr=1.00D-07 GrdErr=1.00D-07 Number of steps in this run= 2 maximum allowed number of steps= 2. GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 7 0 -0.113644 0.606557 -0.000005 2 1 0 0.283786 -0.330622 0.000033 3 1 0 0.283805 1.075114 0.811633 4 1 0 0.283757 1.075168 -0.811660 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 1 N 0.000000 2 H 1.017966 0.000000 3 H 1.017972 1.623203 0.000000 4 H 1.017992 1.623296 1.623293 0.000000 Symmetry turned off by external request. Stoichiometry H3N Framework group C1[X(H3N)] Deg. of freedom 6 Full point group C1 NOp 1 Rotational constants (GHZ): 293.7475344 293.7134104 190.3065162 Standard basis: 6-31G(d,p) (6D, 7F) 30 basis functions, 49 primitive gaussians, 30 cartesian basis functions 5 alpha electrons 5 beta electrons nuclear repulsion energy 11.8944651690 Hartrees. NAtoms= 4 NActive= 4 NUniq= 4 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned off. One-electron integrals computed using PRISM. NBasis= 30 RedAO= T EigKep= 2.18D-02 NBF= 30 NBsUse= 30 1.00D-06 EigRej= -1.00D+00 NBFU= 30 ExpMin= 1.61D-01 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=1 IRadAn= 1 AccDes= 0.00D+00 Harris functional with IExCor= 402 and IRadAn= 1 diagonalized for initial guess. HarFok: IExCor= 402 AccDes= 0.00D+00 IRadAn= 1 IDoV= 1 UseB2=F ITyADJ=14 ICtDFT= 3500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0 NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T wScrn= 0.000000 ICntrl= 500 IOpCl= 0 I1Cent= 200000004 NGrid= 0 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Symmetry not used in FoFCou. Keep R1 ints in memory in canonical form, NReq=989971. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Integral accuracy reduced to 1.0D-05 until final iterations. Initial convergence to 1.0D-05 achieved. Increase integral accuracy. SCF Done: E(RB3LYP) = -56.5577685597 A.U. after 10 cycles NFock= 10 Conv=0.37D-09 -V/T= 2.0091 DoSCS=F DFT=T ScalE2(SS,OS)= 1.000000 1.000000 Range of M.O.s used for correlation: 1 30 NBasis= 30 NAE= 5 NBE= 5 NFC= 0 NFV= 0 NROrb= 30 NOA= 5 NOB= 5 NVA= 25 NVB= 25 Symmetrizing basis deriv contribution to polar: IMax=3 JMax=2 DiffMx= 0.00D+00 G2DrvN: will do 5 centers at a time, making 1 passes. Calling FoFCou, ICntrl= 3107 FMM=F I1Cent= 0 AccDes= 0.00D+00. End of G2Drv F.D. properties file 721 does not exist. End of G2Drv F.D. properties file 722 does not exist. End of G2Drv F.D. properties file 788 does not exist. IDoAtm=1111 Differentiating once with respect to electric field. with respect to dipole field. Differentiating once with respect to nuclear coordinates. Keep R1 ints in memory in canonical form, NReq=967932. There are 15 degrees of freedom in the 1st order CPHF. IDoFFX=6 NUNeed= 3. 12 vectors produced by pass 0 Test12= 8.33D-16 6.67D-09 XBig12= 4.73D+00 1.08D+00. AX will form 12 AO Fock derivatives at one time. 12 vectors produced by pass 1 Test12= 8.33D-16 6.67D-09 XBig12= 3.33D-01 2.77D-01. 12 vectors produced by pass 2 Test12= 8.33D-16 6.67D-09 XBig12= 2.61D-03 1.93D-02. 12 vectors produced by pass 3 Test12= 8.33D-16 6.67D-09 XBig12= 1.35D-06 4.92D-04. 11 vectors produced by pass 4 Test12= 8.33D-16 6.67D-09 XBig12= 3.57D-10 1.05D-05. 4 vectors produced by pass 5 Test12= 8.33D-16 6.67D-09 XBig12= 9.33D-14 1.21D-07. InvSVY: IOpt=1 It= 1 EMax= 8.88D-16 Solved reduced A of dimension 63 with 12 vectors. Isotropic polarizability for W= 0.000000 8.57 Bohr**3. End of Minotr F.D. properties file 721 does not exist. End of Minotr F.D. properties file 722 does not exist. End of Minotr F.D. properties file 788 does not exist. ********************************************************************** Population analysis using the SCF density. ********************************************************************** Alpha occ. eigenvalues -- -14.30568 -0.84466 -0.45031 -0.45029 -0.25317 Alpha virt. eigenvalues -- 0.07985 0.16922 0.16923 0.67850 0.67853 Alpha virt. eigenvalues -- 0.71437 0.87555 0.87556 0.88552 1.13371 Alpha virt. eigenvalues -- 1.41877 1.41880 1.83052 2.09376 2.24216 Alpha virt. eigenvalues -- 2.24223 2.34640 2.34643 2.79253 2.95067 Alpha virt. eigenvalues -- 2.95071 3.19853 3.42892 3.42897 3.90462 Condensed to atoms (all electrons): 1 2 3 4 1 N 6.703110 0.337976 0.337975 0.337975 2 H 0.337976 0.487751 -0.032373 -0.032365 3 H 0.337975 -0.032373 0.487754 -0.032366 4 H 0.337975 -0.032365 -0.032366 0.487743 Mulliken charges: 1 1 N -0.717035 2 H 0.239012 3 H 0.239010 4 H 0.239013 Sum of Mulliken charges = 0.00000 Mulliken charges with hydrogens summed into heavy atoms: 1 1 N 0.000000 APT charges: 1 1 N -0.391057 2 H 0.130348 3 H 0.130349 4 H 0.130360 Sum of APT charges = 0.00000 APT charges with hydrogens summed into heavy atoms: 1 1 N 0.000000 Electronic spatial extent (au): = 39.3615 Charge= 0.0000 electrons Dipole moment (field-independent basis, Debye): X= 1.8464 Y= 0.0000 Z= 0.0000 Tot= 1.8464 Quadrupole moment (field-independent basis, Debye-Ang): XX= -8.7019 YY= -6.1592 ZZ= -6.1590 XY= 1.1199 XZ= 0.0001 YZ= -0.0002 Traceless Quadrupole moment (field-independent basis, Debye-Ang): XX= -1.6952 YY= 0.8475 ZZ= 0.8477 XY= 1.1199 XZ= 0.0001 YZ= -0.0002 Octapole moment (field-independent basis, Debye-Ang**2): XXX= 1.4680 YYY= -11.9765 ZZZ= 0.0000 XYY= 1.4943 XXY= -5.2782 XXZ= 0.0001 XZZ= 0.8151 YZZ= -2.9668 YYZ= -0.0002 XYZ= 0.0000 Hexadecapole moment (field-independent basis, Debye-Ang**3): XXXX= -9.6786 YYYY= -25.1778 ZZZZ= -9.7162 XXXY= 0.8904 XXXZ= 0.0000 YYYX= 1.5793 YYYZ= -0.0002 ZZZX= 0.0000 ZZZY= 0.0000 XXYY= -6.4658 XXZZ= -3.2643 YYZZ= -4.5719 XXYZ= 0.0001 YYXZ= 0.0000 ZZXY= 0.8102 N-N= 1.189446516903D+01 E-N=-1.556684602860D+02 KE= 5.604582094468D+01 Exact polarizability: 6.067 0.000 9.827 0.000 0.000 9.827 Approx polarizability: 7.116 0.000 11.922 0.000 -0.001 11.923 ******************************Gaussian NBO Version 3.1****************************** N A T U R A L A T O M I C O R B I T A L A N D N A T U R A L B O N D O R B I T A L A N A L Y S I S ******************************Gaussian NBO Version 3.1****************************** /RESON / : Allow strongly delocalized NBO set Analyzing the SCF density Job title: NH3_FREQ_ARRON Storage needed: 2904 in NPA, 3721 in NBO (1703935888 available) NATURAL POPULATIONS: Natural atomic orbital occupancies NAO Atom No lang Type(AO) Occupancy Energy ---------------------------------------------------------- 1 N 1 S Cor( 1S) 1.99982 -14.16808 2 N 1 S Val( 2S) 1.53303 -0.57736 3 N 1 S Ryd( 3S) 0.00043 1.20835 4 N 1 S Ryd( 4S) 0.00000 3.73005 5 N 1 px Val( 2p) 1.83298 -0.21388 6 N 1 px Ryd( 3p) 0.00520 0.73498 7 N 1 py Val( 2p) 1.37252 -0.16298 8 N 1 py Ryd( 3p) 0.00158 0.77568 9 N 1 pz Val( 2p) 1.37252 -0.16298 10 N 1 pz Ryd( 3p) 0.00158 0.77569 11 N 1 dxy Ryd( 3d) 0.00163 2.29429 12 N 1 dxz Ryd( 3d) 0.00163 2.29431 13 N 1 dyz Ryd( 3d) 0.00016 2.41119 14 N 1 dx2y2 Ryd( 3d) 0.00149 2.16257 15 N 1 dz2 Ryd( 3d) 0.00060 2.32835 16 H 2 S Val( 1S) 0.62250 0.13597 17 H 2 S Ryd( 2S) 0.00093 0.57862 18 H 2 px Ryd( 2p) 0.00066 2.40557 19 H 2 py Ryd( 2p) 0.00053 2.93335 20 H 2 pz Ryd( 2p) 0.00034 2.31981 21 H 3 S Val( 1S) 0.62250 0.13596 22 H 3 S Ryd( 2S) 0.00093 0.57863 23 H 3 px Ryd( 2p) 0.00066 2.40557 24 H 3 py Ryd( 2p) 0.00039 2.47317 25 H 3 pz Ryd( 2p) 0.00048 2.77997 26 H 4 S Val( 1S) 0.62249 0.13594 27 H 4 S Ryd( 2S) 0.00093 0.57863 28 H 4 px Ryd( 2p) 0.00066 2.40553 29 H 4 py Ryd( 2p) 0.00039 2.47317 30 H 4 pz Ryd( 2p) 0.00048 2.77994 Summary of Natural Population Analysis: Natural Population Natural ----------------------------------------------- Atom No Charge Core Valence Rydberg Total ----------------------------------------------------------------------- N 1 -1.12515 1.99982 6.11104 0.01429 8.12515 H 2 0.37505 0.00000 0.62250 0.00246 0.62495 H 3 0.37505 0.00000 0.62250 0.00246 0.62495 H 4 0.37505 0.00000 0.62249 0.00246 0.62495 ======================================================================= * Total * 0.00000 1.99982 7.97852 0.02166 10.00000 Natural Population -------------------------------------------------------- Core 1.99982 ( 99.9908% of 2) Valence 7.97852 ( 99.7316% of 8) Natural Minimal Basis 9.97834 ( 99.7834% of 10) Natural Rydberg Basis 0.02166 ( 0.2166% of 10) -------------------------------------------------------- Atom No Natural Electron Configuration ---------------------------------------------------------------------------- N 1 [core]2S( 1.53)2p( 4.58)3p( 0.01)3d( 0.01) H 2 1S( 0.62) H 3 1S( 0.62) H 4 1S( 0.62) NATURAL BOND ORBITAL ANALYSIS: Occupancies Lewis Structure Low High Occ. ------------------- ----------------- occ occ Cycle Thresh. Lewis Non-Lewis CR BD 3C LP (L) (NL) Dev ============================================================================= 1(1) 1.90 9.99429 0.00571 1 3 0 1 0 0 0.00 ----------------------------------------------------------------------------- Structure accepted: No low occupancy Lewis orbitals -------------------------------------------------------- Core 1.99982 ( 99.991% of 2) Valence Lewis 7.99447 ( 99.931% of 8) ================== ============================ Total Lewis 9.99429 ( 99.943% of 10) ----------------------------------------------------- Valence non-Lewis 0.00000 ( 0.000% of 10) Rydberg non-Lewis 0.00571 ( 0.057% of 10) ================== ============================ Total non-Lewis 0.00571 ( 0.057% of 10) -------------------------------------------------------- (Occupancy) Bond orbital/ Coefficients/ Hybrids --------------------------------------------------------------------------------- 1. (1.99909) BD ( 1) N 1 - H 2 ( 68.83%) 0.8297* N 1 s( 24.87%)p 3.02( 75.05%)d 0.00( 0.09%) -0.0001 -0.4986 -0.0059 0.0000 -0.2910 0.0052 0.8155 0.0277 0.0000 0.0000 0.0281 0.0000 0.0000 0.0032 0.0082 ( 31.17%) 0.5583* H 2 s( 99.91%)p 0.00( 0.09%) -0.9996 0.0000 0.0072 -0.0289 0.0000 2. (1.99909) BD ( 1) N 1 - H 3 ( 68.83%) 0.8297* N 1 s( 24.86%)p 3.02( 75.05%)d 0.00( 0.09%) 0.0001 0.4986 0.0059 0.0000 0.2910 -0.0052 0.4077 0.0138 0.7062 0.0240 0.0140 0.0243 0.0076 0.0033 0.0031 ( 31.17%) 0.5583* H 3 s( 99.91%)p 0.00( 0.09%) 0.9996 0.0000 -0.0072 -0.0145 -0.0250 3. (1.99909) BD ( 1) N 1 - H 4 ( 68.83%) 0.8297* N 1 s( 24.87%)p 3.02( 75.05%)d 0.00( 0.09%) 0.0001 0.4986 0.0059 0.0000 0.2909 -0.0052 0.4077 0.0138 -0.7062 -0.0239 0.0140 -0.0243 -0.0076 0.0033 0.0031 ( 31.17%) 0.5583* H 4 s( 99.91%)p 0.00( 0.09%) 0.9996 0.0000 -0.0072 -0.0145 0.0250 4. (1.99982) CR ( 1) N 1 s(100.00%) 1.0000 -0.0002 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 5. (1.99721) LP ( 1) N 1 s( 25.38%)p 2.94( 74.52%)d 0.00( 0.10%) 0.0001 0.5036 -0.0120 0.0000 -0.8618 0.0505 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0269 0.0155 6. (0.00000) RY*( 1) N 1 s( 99.98%)p 0.00( 0.02%)d 0.00( 0.00%) 7. (0.00000) RY*( 2) N 1 s(100.00%) 8. (0.00000) RY*( 3) N 1 s( 0.03%)p99.99( 99.97%)d 0.01( 0.00%) 9. (0.00000) RY*( 4) N 1 s( 0.00%)p 1.00(100.00%)d 0.00( 0.00%) 10. (0.00000) RY*( 5) N 1 s( 0.00%)p 1.00(100.00%)d 0.00( 0.00%) 11. (0.00000) RY*( 6) N 1 s( 0.00%)p 1.00( 0.12%)d99.99( 99.88%) 12. (0.00000) RY*( 7) N 1 s( 0.00%)p 1.00( 0.12%)d99.99( 99.88%) 13. (0.00000) RY*( 8) N 1 s( 0.00%)p 1.00( 0.01%)d99.99( 99.99%) 14. (0.00000) RY*( 9) N 1 s( 0.01%)p 4.42( 0.06%)d99.99( 99.92%) 15. (0.00000) RY*(10) N 1 s( 0.00%)p 1.00( 0.03%)d99.99( 99.97%) 16. (0.00112) RY*( 1) H 2 s( 72.77%)p 0.37( 27.23%) 0.0038 0.8531 0.5218 -0.0017 -0.0001 17. (0.00045) RY*( 2) H 2 s( 26.59%)p 2.76( 73.41%) -0.0017 0.5157 -0.8435 -0.1501 -0.0007 18. (0.00034) RY*( 3) H 2 s( 0.00%)p 1.00(100.00%) 0.0000 0.0004 -0.0005 -0.0001 1.0000 19. (0.00000) RY*( 4) H 2 s( 0.72%)p99.99( 99.28%) 20. (0.00112) RY*( 1) H 3 s( 72.77%)p 0.37( 27.23%) 0.0038 0.8531 0.5218 0.0009 0.0014 21. (0.00045) RY*( 2) H 3 s( 26.59%)p 2.76( 73.41%) -0.0017 0.5157 -0.8435 0.0756 0.1297 22. (0.00034) RY*( 3) H 3 s( 0.00%)p 1.00(100.00%) 0.0000 0.0004 -0.0005 -0.8660 0.5000 23. (0.00000) RY*( 4) H 3 s( 0.72%)p99.99( 99.28%) 24. (0.00112) RY*( 1) H 4 s( 72.76%)p 0.37( 27.24%) 0.0038 0.8530 0.5219 0.0009 -0.0015 25. (0.00045) RY*( 2) H 4 s( 26.61%)p 2.76( 73.39%) -0.0017 0.5158 -0.8435 0.0750 -0.1300 26. (0.00034) RY*( 3) H 4 s( 0.00%)p 1.00(100.00%) 0.0000 0.0000 0.0000 0.8660 0.5000 27. (0.00000) RY*( 4) H 4 s( 0.72%)p99.99( 99.28%) 28. (0.00000) BD*( 1) N 1 - H 2 ( 31.17%) 0.5583* N 1 s( 24.87%)p 3.02( 75.05%)d 0.00( 0.09%) ( 68.83%) -0.8297* H 2 s( 99.91%)p 0.00( 0.09%) 29. (0.00000) BD*( 1) N 1 - H 3 ( 31.17%) 0.5583* N 1 s( 24.86%)p 3.02( 75.05%)d 0.00( 0.09%) ( 68.83%) -0.8297* H 3 s( 99.91%)p 0.00( 0.09%) 30. (0.00000) BD*( 1) N 1 - H 4 ( 31.17%) 0.5583* N 1 s( 24.87%)p 3.02( 75.05%)d 0.00( 0.09%) ( 68.83%) -0.8297* H 4 s( 99.91%)p 0.00( 0.09%) NHO Directionality and "Bond Bending" (deviations from line of nuclear centers) [Thresholds for printing: angular deviation > 1.0 degree] hybrid p-character > 25.0% orbital occupancy > 0.10e Line of Centers Hybrid 1 Hybrid 2 --------------- ------------------- ------------------ NBO Theta Phi Theta Phi Dev Theta Phi Dev ======================================================================================== 1. BD ( 1) N 1 - H 2 90.0 293.0 90.0 288.7 4.3 -- -- -- 2. BD ( 1) N 1 - H 3 37.1 49.7 34.9 55.9 4.3 -- -- -- 3. BD ( 1) N 1 - H 4 142.9 49.7 145.1 55.9 4.3 -- -- -- 5. LP ( 1) N 1 -- -- 90.0 180.0 -- -- -- -- Second Order Perturbation Theory Analysis of Fock Matrix in NBO Basis Threshold for printing: 0.50 kcal/mol E(2) E(j)-E(i) F(i,j) Donor NBO (i) Acceptor NBO (j) kcal/mol a.u. a.u. =================================================================================================== within unit 1 5. LP ( 1) N 1 / 16. RY*( 1) H 2 1.01 1.43 0.034 5. LP ( 1) N 1 / 17. RY*( 2) H 2 0.67 2.17 0.034 5. LP ( 1) N 1 / 20. RY*( 1) H 3 1.01 1.43 0.034 5. LP ( 1) N 1 / 21. RY*( 2) H 3 0.67 2.17 0.034 5. LP ( 1) N 1 / 24. RY*( 1) H 4 1.01 1.43 0.034 5. LP ( 1) N 1 / 25. RY*( 2) H 4 0.67 2.17 0.034 Natural Bond Orbitals (Summary): Principal Delocalizations NBO Occupancy Energy (geminal,vicinal,remote) ==================================================================================== Molecular unit 1 (H3N) 1. BD ( 1) N 1 - H 2 1.99909 -0.60417 2. BD ( 1) N 1 - H 3 1.99909 -0.60417 3. BD ( 1) N 1 - H 4 1.99909 -0.60416 4. CR ( 1) N 1 1.99982 -14.16768 5. LP ( 1) N 1 1.99721 -0.31756 24(v),16(v),20(v),17(v) 21(v),25(v) 6. RY*( 1) N 1 0.00000 1.20799 7. RY*( 2) N 1 0.00000 3.73005 8. RY*( 3) N 1 0.00000 0.73750 9. RY*( 4) N 1 0.00000 0.77341 10. RY*( 5) N 1 0.00000 0.77342 11. RY*( 6) N 1 0.00000 2.29040 12. RY*( 7) N 1 0.00000 2.29041 13. RY*( 8) N 1 0.00000 2.40939 14. RY*( 9) N 1 0.00000 2.16324 15. RY*( 10) N 1 0.00000 2.32730 16. RY*( 1) H 2 0.00112 1.11323 17. RY*( 2) H 2 0.00045 1.84851 18. RY*( 3) H 2 0.00034 2.31981 19. RY*( 4) H 2 0.00000 2.94718 20. RY*( 1) H 3 0.00112 1.11322 21. RY*( 2) H 3 0.00045 1.84851 22. RY*( 3) H 3 0.00034 2.31980 23. RY*( 4) H 3 0.00000 2.94717 24. RY*( 1) H 4 0.00112 1.11347 25. RY*( 2) H 4 0.00045 1.84827 26. RY*( 3) H 4 0.00034 2.31978 27. RY*( 4) H 4 0.00000 2.94712 28. BD*( 1) N 1 - H 2 0.00000 0.48620 29. BD*( 1) N 1 - H 3 0.00000 0.48619 30. BD*( 1) N 1 - H 4 0.00000 0.48616 ------------------------------- Total Lewis 9.99429 ( 99.9429%) Valence non-Lewis 0.00000 ( 0.0000%) Rydberg non-Lewis 0.00571 ( 0.0571%) ------------------------------- Total unit 1 10.00000 (100.0000%) Charge unit 1 0.00000 Calling FoFJK, ICntrl= 100127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. Full mass-weighted force constant matrix: Low frequencies --- -30.7764 -0.0016 -0.0014 -0.0011 20.3142 28.2484 Low frequencies --- 1089.5557 1694.1237 1694.1868 Diagonal vibrational polarizability: 3.2993697 0.1276763 0.1276583 Harmonic frequencies (cm**-1), IR intensities (KM/Mole), Raman scattering activities (A**4/AMU), depolarization ratios for plane and unpolarized incident light, reduced masses (AMU), force constants (mDyne/A), and normal coordinates: 1 2 3 A A A Frequencies -- 1089.5557 1694.1237 1694.1868 Red. masses -- 1.1800 1.0644 1.0644 Frc consts -- 0.8253 1.8000 1.8001 IR Inten -- 145.4398 13.5557 13.5560 Atom AN X Y Z X Y Z X Y Z 1 7 0.12 0.00 0.00 0.00 -0.02 -0.06 0.00 0.06 -0.02 2 1 -0.53 -0.21 0.00 -0.07 -0.04 0.73 0.25 0.14 0.19 3 1 -0.53 0.11 0.18 0.25 -0.24 -0.02 -0.07 -0.61 0.40 4 1 -0.53 0.11 -0.18 -0.18 0.51 0.18 -0.18 -0.41 -0.36 4 5 6 A A A Frequencies -- 3460.9808 3589.3959 3589.5242 Red. masses -- 1.0272 1.0883 1.0883 Frc consts -- 7.2496 8.2615 8.2621 IR Inten -- 1.0593 0.2700 0.2708 Atom AN X Y Z X Y Z X Y Z 1 7 -0.04 0.00 0.00 0.00 -0.03 0.07 0.00 0.07 0.03 2 1 0.18 -0.55 0.00 -0.13 0.32 0.02 0.28 -0.68 0.01 3 1 0.18 0.27 0.47 -0.18 -0.23 -0.36 -0.25 -0.30 -0.54 4 1 0.18 0.27 -0.48 0.31 0.38 -0.65 -0.02 -0.01 0.06 ------------------- - Thermochemistry - ------------------- Temperature 298.150 Kelvin. Pressure 1.00000 Atm. Atom 1 has atomic number 7 and mass 14.00307 Atom 2 has atomic number 1 and mass 1.00783 Atom 3 has atomic number 1 and mass 1.00783 Atom 4 has atomic number 1 and mass 1.00783 Molecular mass: 17.02655 amu. Principal axes and moments of inertia in atomic units: 1 2 3 Eigenvalues -- 6.14385 6.14457 9.48334 X 0.00005 0.00002 1.00000 Y -0.51365 0.85800 0.00001 Z 0.85800 0.51365 -0.00005 This molecule is an asymmetric top. Rotational symmetry number 1. Rotational temperatures (Kelvin) 14.09764 14.09600 9.13326 Rotational constants (GHZ): 293.74753 293.71341 190.30652 Zero-point vibrational energy 90424.3 (Joules/Mol) 21.61193 (Kcal/Mol) Vibrational temperatures: 1567.63 2437.46 2437.55 4979.57 5164.33 (Kelvin) 5164.52 Zero-point correction= 0.034441 (Hartree/Particle) Thermal correction to Energy= 0.037304 Thermal correction to Enthalpy= 0.038248 Thermal correction to Gibbs Free Energy= 0.015365 Sum of electronic and zero-point Energies= -56.523328 Sum of electronic and thermal Energies= -56.520465 Sum of electronic and thermal Enthalpies= -56.519521 Sum of electronic and thermal Free Energies= -56.542404 E (Thermal) CV S KCal/Mol Cal/Mol-Kelvin Cal/Mol-Kelvin Total 23.408 6.326 48.162 Electronic 0.000 0.000 0.000 Translational 0.889 2.981 34.441 Rotational 0.889 2.981 13.646 Vibrational 21.631 0.364 0.075 Q Log10(Q) Ln(Q) Total Bot 0.856667D-07 -7.067188 -16.272802 Total V=0 0.594908D+09 8.774450 20.203917 Vib (Bot) 0.144835D-15 -15.839126 -36.470936 Vib (V=0) 0.100580D+01 0.002512 0.005783 Electronic 0.100000D+01 0.000000 0.000000 Translational 0.276150D+07 6.441145 14.831285 Rotational 0.214187D+03 2.330793 5.366849 ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 7 0.000007148 0.000001537 -0.000004076 2 1 0.000003163 0.000004032 -0.000009139 3 1 -0.000004158 0.000006306 -0.000008284 4 1 -0.000006153 -0.000011875 0.000021498 ------------------------------------------------------------------- Cartesian Forces: Max 0.000021498 RMS 0.000008881 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Search for a local minimum. Step number 1 out of a maximum of 2 All quantities printed in internal units (Hartrees-Bohrs-Radians) Second derivative matrix not updated -- analytic derivatives used. The second derivative matrix: X1 Y1 Z1 X2 Y2 X1 0.22811 Y1 -0.00002 0.63144 Z1 0.00004 0.00000 0.63146 X2 -0.07603 0.17858 0.00000 0.07581 Y2 0.11892 -0.36057 0.00001 -0.14161 0.39653 Z2 0.00000 0.00001 -0.06041 0.00000 -0.00001 X3 -0.07605 -0.08929 -0.15462 0.00011 0.01135 Y3 -0.05946 -0.13543 -0.12999 -0.01848 -0.01798 Z3 -0.10299 -0.12997 -0.28554 -0.00243 0.00278 X4 -0.07603 -0.08927 0.15459 0.00011 0.01134 Y4 -0.05944 -0.13544 0.12998 -0.01848 -0.01797 Z4 0.10295 0.12997 -0.28551 0.00243 -0.00278 Z2 X3 Y3 Z3 X4 Z2 0.05985 X3 0.01479 0.07584 Y3 -0.03438 0.07080 0.14401 Z3 0.00027 0.12263 0.14578 0.31236 X4 -0.01479 0.00010 0.00714 -0.01721 0.07582 Y4 0.03438 0.00714 0.00940 -0.01858 0.07078 Z4 0.00028 0.01721 0.01858 -0.02709 -0.12259 Y4 Z4 Y4 0.14401 Z4 -0.14577 0.31233 ITU= 0 Eigenvalues --- 0.09784 0.13744 0.13745 0.55420 0.86365 Eigenvalues --- 0.86368 Angle between quadratic step and forces= 40.46 degrees. Linear search not attempted -- first point. TrRot= -0.000002 -0.000004 0.000008 0.000001 0.000003 0.000001 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) X1 -0.21476 0.00001 0.00000 -0.00001 -0.00001 -0.21477 Y1 1.14623 0.00000 0.00000 -0.00001 -0.00001 1.14621 Z1 -0.00001 0.00000 0.00000 0.00002 0.00002 0.00001 X2 0.53628 0.00000 0.00000 0.00000 0.00000 0.53628 Y2 -0.62479 0.00000 0.00000 0.00001 0.00000 -0.62478 Z2 0.00006 -0.00001 0.00000 -0.00007 -0.00007 0.00000 X3 0.53631 0.00000 0.00000 0.00000 0.00000 0.53631 Y3 2.03167 0.00001 0.00000 0.00006 0.00006 2.03173 Z3 1.53376 -0.00001 0.00000 -0.00004 -0.00003 1.53373 X4 0.53622 -0.00001 0.00000 0.00002 0.00001 0.53624 Y4 2.03177 -0.00001 0.00000 -0.00004 -0.00005 2.03173 Z4 -1.53382 0.00002 0.00000 0.00007 0.00008 -1.53374 Item Value Threshold Converged? Maximum Force 0.000021 0.000450 YES RMS Force 0.000009 0.000300 YES Maximum Displacement 0.000077 0.001800 YES RMS Displacement 0.000039 0.001200 YES Predicted change in Energy=-1.603127D-09 Optimization completed. -- Stationary point found. GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad 1\1\GINC-CX1-15-34-2\Freq\RB3LYP\6-31G(d,p)\H3N1\SCAN-USER-1\09-Oct-20 13\0\\# freq b3lyp/6-31g(d,p) nosymm pop=nbo geom=connectivity\\NH3_FR EQ_ARRON\\0,1\N,-0.113644,0.606557,-0.000005\H,0.283786,-0.330622,0.00 0033\H,0.283805,1.075114,0.811633\H,0.283757,1.075168,-0.81166\\Versio n=ES64L-G09RevD.01\HF=-56.5577686\RMSD=3.657e-10\RMSF=8.881e-06\ZeroPo int=0.0344408\Thermal=0.0373037\Dipole=0.7264302,-0.000006,-0.0000002\ DipoleDeriv=-0.5555003,0.0000048,0.0000045,-0.0000116,-0.3088422,-0.00 0005,0.000013,-0.000005,-0.3088294,0.1851579,0.1861313,0.000005,0.0937 758,0.0445709,-0.0000066,-0.0000033,0.0000091,0.1613149,0.1851631,-0.0 930736,-0.1612002,-0.046883,0.1321311,-0.0505415,-0.0812171,-0.0505573 ,0.0737534,0.1851793,-0.0930624,0.1611907,-0.0468812,0.1321402,0.05055 31,0.0812074,0.0505532,0.0737611\Polar=6.067338,0.0000108,9.8266779,0. 0001984,-0.0003454,9.8271686\PG=C01 [X(H3N1)]\NImag=0\\0.22811389,-0.0 0001710,0.63143775,0.00003745,-0.00000208,0.63146234,-0.07603403,0.178 57635,-0.00000269,0.07581230,0.11892066,-0.36057441,0.00001107,-0.1416 1013,0.39652916,-0.00000282,0.00001078,-0.06040532,0.00000460,-0.00000 679,0.05985487,-0.07604981,-0.08928810,-0.15462034,0.00011266,0.011345 55,0.01478512,0.07583504,-0.05946014,-0.13542770,-0.12998736,-0.018483 11,-0.01798130,-0.03437959,0.07079966,0.14401310,-0.10298770,-0.129973 85,-0.28554235,-0.00242704,0.00277507,0.00027419,0.12262520,0.14578243 ,0.31235691,-0.07603005,-0.08927116,0.15458558,0.00010907,0.01134393,- 0.01478690,0.00010212,0.00714360,-0.01721046,0.07581886,-0.05944341,-0 .13543565,0.12997836,-0.01848311,-0.01797346,0.03437560,0.00714289,0.0 0939589,-0.01858365,0.07078363,0.14401321,0.10295306,0.12996515,-0.285 51466,0.00242513,-0.00277935,0.00027626,0.01721003,0.01858452,-0.02708 874,-0.12258822,-0.14577031,0.31232714\\-0.00000715,-0.00000154,0.0000 0408,-0.00000316,-0.00000403,0.00000914,0.00000416,-0.00000631,0.00000 828,0.00000615,0.00001188,-0.00002150\\\@ EXPERIMENTALISTS THINK SILICON IS REALLY FUN TO USE ITS PLACE IN NOVEL COMPOUNDS IS CERTAIN TO AMUSE THEY SIT ALL DAY IN LABORATORIES MAKING ALL THIS SLUDGE "LOADED WITH THE SILICON THEY SAY", TO ME IT LOOKS LIKE FUDGE. FOR HAPPY THOUGH THEY BE WITH CRUD, I'D LIKE TO KNOW A LITTLE ABOUT THE PI BONDS ON THE EDGE AND SIGMAS IN THE MIDDLE. SO LETS DERIVE A WAVEFUNCTION.....6-31G* USE AN OPTIMAL GEOMETRY AND SEE WHERE ELECTRONS ARE. BUT WHAT OF CORRELATION? ASKS THE WIRY LITTLE SKEPTIC. WE'LL THROW IN PERTURBATION AS AN ELECTRON ANTISEPTIC. AND WHEN THE PROGRAM GIVES US ANSWERS IN THEM WE CAN TRUST SINCE NOBODY CAN MAKE THE STUFF, WE HAVE NO CHOICE, WE MUST. SO THEORY GUYS HAVE GOT IT MADE, IN ROOMS FREE OF POLLUTION. INSTEAD OF PROBLEMS WITH THE REFLUX, THEY HAVE ONLY SOLUTIONS. AND WHEN THE FEDS ANNOUNCE THE LIST OF CARCINOGENIC TERRORS, THE THEORISTS SIT SAFELY AT THEIR TERMINALS FIXING ERRORS. IN OTHER WORDS, EXPERIMENTALISTS WILL LIKELY DIE OF CANCER FROM WORKING HARD YET FRUITLESSLY...TILL THEORY GIVES THE ANSWER. -- THOMAS A. HOLME, 1983 Job cpu time: 0 days 0 hours 0 minutes 49.4 seconds. File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1 Normal termination of Gaussian 09 at Wed Oct 9 14:40:04 2013.