Personal Info
Full Name: Joaquin Barroso-Flores
Born: March 15th, 1978 Mexico City
Education:
- Chemistry Faculty, National Autonomous University of Mexico (UNAM). Chemist (graduated with honors, 2001)
- Chemistry Faculty, National Autonomous University of Mexico (UNAM). PhD in chemistry (2001 – 2005)
Current work: Associate researcher at the Joint Center for Research in Sustainable Chemistry (CCIQS), which is part of the Institute of Chemistry from the National Autonomous University of Mexico (UNAM). This center is located at the city of Toluca, 100 km west from Mexico City. (Link to official site on Instituto de Química)
Current projects: DFT and Molecular Dynamics calculations on calixarene systems as drug delivery agents. Excitonic transference in photosynthesis. Non-canonical DNA. Organic reaction mechanisms driven by unusual quantum effects.
ResearchGate profile: Dr. Joaquín Barroso-Flores
Academia.edu profile: Dr. Joaquín Barroso-Flores
Google Scholar Profile: Dr. Joaquín Barroso-Flores
Open Researcher & Contributor ID (ORCID): 0000-0003-0554-7569
ResearcherID: D-1790-2011
Follow us on Twitter! @joaquinbarroso
Professional Experience:
- Associate Researcher at Instituto de Química, UNAM. Mexico. (May 2010 – present)
- Postdoc researcher at Facultatea de Chemie si Ingenieria Chimica, Universitatea Babes-Bolyai; Cluj-Napoca, Romania. Click on the following link to find out more about the MOLSEN project at Babes-Bolyai. (Oct. 2008 – Jan 2010)
- Associate Research Fellow at Centro de Investigacion en Polimeros (COMEXgroup) (Jan. 2006 – Oct. 2008)
- TA electromagnetism with Prof. Dr. Germund Höjer at the Chemistry Faculty, UNAM (2003 – 2005)
- High School Chemistry teacher at Escuela Preparatoria Logos (Chemistry I & II) (2002 – 2003)
Research interests: Theoretical chemistry; Computational chemistry; Description of molecular properties in weakly interacting systems through assessment of electron density distribution.
Published papers (Not necessarily an exhaustive/updated list):
- Jacinto Sandoval-Lira, Gustavo Mondragón-Solórzano, Leonardo I Lugo-Fuentes, Joaquín Barroso-Flores* “Accurate Estimation of pKb Values for Amino Groups from Surface Electrostatic Potential (VS, min) Calculations: The Isoelectric Points of Amino Acids as a Case Study” J. Chem. Inf. Model. 2020, 60, 3, 1445-1452
- Reynosa-Martínez, A. C.; Navarro-Tovar, G.; Gallegos, W. R.; Rodríguez-Meléndez, H.; Torres-Cadena, R.; Mondragón-Solórzano, G.; Barroso-Flores, J.; Álvarez-Lemus, M. A.; García-Montalvo, V.; López-Honorato, E. “Effect of the degree of oxidation of Graphene oxide on As(III) adsorption” J. Hazard. Mater. 2020, 384, 15, 121440 DOI: 10.1016/j.jhazmat.2019.121440
- Medina-Mercado, I.; Asomoza-Solís, E. O.; Martínez-González, E.; Ugalde-Saldivar, V. M.; Ledesma-Olvera, L. G.; Barquera-Lozada, J. E.; Gómez-Vidales, V.; Barroso-Flores, J.; Frontana-Uribe, B. A.; Porcel, S.* “Ascorbic Acid as an Aryl Radical Inducer in the Gold-Mediated Arylation of Indoles with Aryldiaonium Chlorides” Chem. Eur. J. 2020, 26, 634-642 DOI: 10.1002/chem.201904413
- Mancillas-Salas, S.; Barroso-Flores, J.; Villaurrutia, R.; García-Montalvo, V.; López-Honorato, E. “Production of few-layer graphene by wet media milling using organic solvents and different types of graphite” Ceram. Int. 2020, 46, 2413–2420 DOI: 10.1016/j.ceramint.2019.09.235
- Galindo-Murillo, R.; Barroso-Flores, J* “Hydrophobic unnatural base pairs show a Watson-Crick pairing in micro-second molecular dynamics simulations”. J. Biomol. Struct. Dyn. 2019, DOI: 10.1080/07391102.2019.1671898
- Hernandez-Martinez, Á. R.; Molina, G. A.; Rodríguez-Torres, A.; Ledesma-Mendoza, B.; Del Real, A.; Barroso-Flores, J.; Estevez, M. “Fluorescence decay rate of selected compounds from Eysenhardtia polystachya extracts and their viability as biosensors” Mater. Sci. Eng. C 2019, 104 (July), 109978-91.
- Caballero-García, G.; Mondragón-Solórzano, G.; Torres-Cadena, R.; Díaz-García, M.; Sandoval-Lira, J.; Barroso-Flores, J.* “Calculation of VS,max and its use as a descriptor for the Theoretical Calculation of pKa values for Carboxylic Acids” Molecules, 2019, 24(1) 79
- Salomón-Flores, M. K.; Hernández-Juárez, C. L.; Bazany-Rodríguez, I. J.; Barroso-Flores, J.; Martínez-Otero, D.; López-Arteaga, R.; Valdés-Martínez, J.; Dorazco-González, A.* “Efficient fluorescent chemosensing of iodide based on a cationic meso-tetraarylporphyrin in pure water” Sensors Actuators, B Chem. 2019, 281, 462–470
- Silva-Nigenda, E.; Martínez-Gómez, A.; Cruz-de la Cruz, J.; Barroso-Flores, J.; González-Romero, C.; Fuentes-Benites, A.; Jankowski, C. K.; Cuevas-Yáñez, E.; Díaz-Torres, E.; Corona-Becerril, D.* “Long range1H–19F coupling through multiple bond in thienopyridines, isoquinolines and 2-aza-carbazoles derivatives.” J. Mol. Struct. 2019, 1176, 562–566 IF 2.011
- Cheng, H. F.; d’Aquino, A. I.; Barroso-Flores, J.; Mirkin, C. A.* “A Redox-Switchable, Allosteric Coordination Complex” J. Am. Chem. Soc. 2018, 140, 14590–14594
- A.; Heger, Z.; Milosavljevic, V.; Richtera, L.; Barroso-Flores, J.; Merlos Rodrigo, M. A.; Buchtelova, H.; Podgajny, R.; Hynek, D.; Kopel, P.; Adam, V. “Real-Time Visualization of Cell Membrane Damage Using Gadolinium–Schiff Base Complex-Doped Quantum Dots” ACS Appl.Mater. Interfaces 2018, 10, 35859–35868
- Mondragón-Solórzano, G.; Barroso-Flores, J.* “Spectroscopical UV-Vis implications of an intramolecular η 2 -Mg coordination in bacteriochlorophyll-a from the Fenna-Matthews-Olson complex” Int. J. Quantum Chem. 2018, 118, e25663
- Andrea I. d’Aquino, Ho Fung Cheng, Joaquin Barroso-Flores, Zachary S. Kean, Jose Mendez-Arroyo, C. Michael McGuirk, Chad A. Mirkin*, “An Allosterically Regulated, Four-State Macrocycle” Inorg. Chem., (2018) doi: 10.1021/acs.inorgchem.7b02745
- Joaquín Barroso-Flores* “Evolution of the Fenna–Matthews–Olson Complex and Its Quantum Coherence Features. Which Led the Way?” ACS Central Science (2017) 3 (10), 1061-1062
- Erandi Bernabé-Pablo, Vojtech Jancik, Diego Martínez-Otero, Joaquín Barroso-Flores, Mónica Moya-Cabrera* “Molecular Group 13 Metallaborates Derived from M–O–M Cleavage Promoted by BH3” Inorg. Chem., (2017), 56 (14), 7890–7899 DOI: 10.1021/acs.inorgchem.7b00634
- Ana Cristina García-Álvarez, Erandi Bernabé-Pablo, Joaquín Barroso-Flores, Vojtech Jancik, Diego Martínez-Otero, T. Jesús Morales-Juárez, MónicaMoya-Cabrera* “Multinuclear rare-earth metal complexes supported by chalcogen-based 1,2,3-triazole” Polyhedron (2017) 135, 10-16 DOI: 10.1016/j.poly.2017.06.047
- Rodrigo Galindo-Murillo, Joaquín Barroso-Flores* “Structural and Dynamical Instability of DNA Caused by High Occurrence of d5SICS and dNaM Unnatural Nucleotides” Phys. Chem. Chem. Phys. (2017) 19, 10571-10580 DOI: 10.1039/C7CP01477E
- Caballero-García, Guillermo; Romero-Ortega, Moisés, Barroso-Flores, Joaquín* “Reactivity of electrophilic Chlorine Atoms Due to σ-holes. A Mechanistic Assessment of the Chemical Reduction of the Trichloromethyl Group by Sulfur Nucleophiles” Phys. Chem. Chem. Phys. 2016 18, 27300-27307
- Caballero-García, Guillermo; Romero-Ortega, Moisés; Barroso-Flores, Joaquín* “Exocyclic delocalization inhibits the aromatization of pyridinylidenes. A Theoretical assesment of the electronic density and the associated transition states” Tetrahedron (2016) 72, 29, 4194-4200
- Mondragón Solórzano, Gustavo; López-Honorato, Eddie; Sierra-Álvarez, Reyes; Barroso –Flores, Joaquín* “In Silico Design of Calixarene-Based Arsenic Acid Removal Agents” J. Incl. Phenom. (2016) 85, 1, 169-174
- Seballos-Resendiz, Arturo; Lechuga-Eduardo, Harim; Martínez-Otero, Diego; Barroso-Flores, Joaquín; Romero-Ortega, Moisés “Synthesis and Crystal Structures of Stable 4-Aryl-2-trichloromethyl-1,3-diaza-1,3-butadienes” Synthesis (2016) Online First DOI: 10.1055/s-0035-1561936
- Galindo-Murillo, Rodrigo; Aguilar-Suárez, Luis Enrique; Barroso-Flores, Joaquín* “A Mixed DFT – MD Methodology for the In Silico Development of Drug Releasing Macrocycles. Calix and thia-calix[n]arenes as Carriers for Bosutinib and Sorafenib” J. Comp. Chem. (2016) 37, 10, 940-946
- Lukas Richtera,*, Vojtech Jancik,*, Joaquín Barroso-Flores, Petr Nykel, Jiri Touzin and Jan Taraba “Synthesis and Crystal Structure of the First Selenonyl Bis(carboxylate) SeO2(O2CCH3)2” Eur. J. Inorg. Chem. (2015) 18, 2923-2927
- Bazany-Rodríguez, I. J., Martínez-Otero, D., Barroso-Flores, J, Yatsimirsky, A. K., & Dorazco-González, A. ” J. Comp. Chem. (2015) Accepted.Sensitive water-soluble fluorescent chemosensor for chloride based on a bisquinolinium pyridine-dicarboxamide compound”. Sensors and Actuators B: Chemical, (2015) 221, 1348–1355
- Mendez-Arroyo, J.; Barroso-Flores, J.; Lifschitz, A. M.; Sarjeant, A. a; Stern, C. L.; Mirkin, C. a. A Multi-State, Allosterically-Regulated Molecular Receptor With Switchable Selectivity. J. Am. Chem. Soc. 2014, 136, 10340–10348.
- Zarabadi-poor, P.; Barroso-Flores, J. Theoretical Assessment of the Selective Fluorescence Quenching of 1‑Amino-8-Naphthol-3,6-Disulfonic Acid (H-Acid) Complexes with Zn2+,Cd2+, and Hg2+: A DFT and TD-DFT Study. J. Phys. Chem. A 2014, 12178–12183.
- Galindo-Murillo, R.; Olmedo-Romero, A.; Cruz-Flores, E.; Petrar, P. M.; Kunsagi-Mate, S.; Barroso-Flores, J. Calix[n]arene-Based Drug Carriers: A DFT Study of Their Electronic Interactions with a Chemotherapeutic Agent Used against Leukemia. Comput. Theor. Chem. 2014, 1035, 84–91.
- Rodrigo Galindo-Murillo, María Eugenia Sandoval-Salinas, and Joaquín Barroso-Flores “In Silico Design of Monomolecular Drug Carriers for the Tyrosine Kinase Inhibitor Drug Imatinib Based on Calix- and Thiacalix[n]arene Host Molecules: A DFT and Molecular Dynamics Study” J. Chem. Theory Comp. (2014) DOI: 10.1021/ct4004178
- Pezhman Zarabadi-Poor; Alireza Badiei; Ali Akbar Yousefi; Joaquín Barroso-Flores “SELECTIVE OPTICAL SENSING OF HG(II) IN AQUEOUS MEDIA BY H-ACID/SBA-15: A COMBINED EXPERIMENTAL AND THEORETICAL STUDY” Journal of Chemical Physics C 2013, 117 (18), 9281-9289
- Sandra Hidalgo-Bonilla; Ricardo Peyrot; Vojtech Jancik; Joaquín Barroso-Flores; Marisol Reyes-Lezama; Mónica Moya-Cabrera “MOLECULAR HETEROBIMETALLIC ALUMINOXANES AND ALUMINOXANE SULFIDES CONTAINING GROUP 4 METALS” European Journal of Inorganic Chemistry 2013, 16, 2849-2857
- Azucena Garduño-Alva, M. Carmen Ortega-Alfaro, José G. López-Cortés, Isabel Chávez, Joaquin Barroso-Flores, Rubén A. Toscano, Henri Rudler, and Cecilio Álvarez-Toledano “Synthesis of new γ-lactones from preactivated monosubstituted pyrazines and TMS–ketene acetals” Can. J. Chem. 2012, 90 (5): 469-482
- Joaquín Barroso-Flores, Petronela Maria Petrar; Ioan Silaghi-Dumitrescu, Sandor Kunsagi-Mate “Ab initio calculations of electronic interactions in inclusion comoplexes of calix- and thiacalix[n]arenes and block s cations” J. Incl. Phenom. Macrocycl. Chem. 2012 online first DOI: 10.1007/s10847-012-0144-6
- Jocelyn Alcántara-García, Vojtech Jancik, Joaquín Barroso, Raymundo Cea-Olivares, Mónica Moya-Cabrera, “Coordination Diversity of Aluminum centers Molded by Triazole Based Chalcogen Ligands” Inorganic Chemistry (2009) 48, 13 5545-5558.
- Joaquín Barroso-Flores, J. A. Cogordan “Influence of Intramolecular Sn-Chalcogen interactions on the conformational preferences for three Diorganotin (IV) Xanthates.” Journal of Organometallic Chemistry 691 (2006) 4937-4944.
- Óscar Baldovino-Pantaleón, Joaquín Barroso-Flores, J.A. Cogordan, Simón Hernández-Ortega, Rubén A. Toscano, David Morales-Morales, “Phosphane-free C-C Heck couplings catalyzed by Pd(II) fluorinated aniline complexes of the type trans-[PdCl2(NH2ArF)2]” Journal of Molecular Catalysis A. Chemical 247 (2005) 65-72.
- Joaquín Barroso-Flores, Raymundo Cea-Olivares, Rubén A. Toscano, J. A. Cogordan, “Synthesis of the anisobidentate compound bis(2-amino-cyclopent-1-ene-carbodithioate)diethyltin (IV). Experimental and theoretical study” Journal of Organometallic Chemistry 689 (2004) 2096-2102.
Dr. Joaquin Barroso's Blog 
Apreciado Joaquin: Conozco algunos de los próceres que mencionas en el Bolg. Cea Olivares, e algunosde los vampiritos. Haiduc tal vez. Estuve varias veces en México, una de ellas convidaddo por Jacobo Gomez Lara, ya fallecido, otras por Victor Castaño. bueno..
mi pregunta esla siguiente como buen novato. Estudiando el complexo bis-dietilditiocarbamato de Cu(II) enconbtré transferencia de carga entre:
Donador ACEEPTOR ENERGIA
85 LP(2) S8 LP* Cu1 31.14
esto está muy padre y pudeo visualizar el orbital 85
en Gaussian View y en el Chemcraft. Como visualizao el orbital 83? Si es antiligante deve estar arriba de l LUMO y no abajo del HOMO.
Puede usted ayudarme co la ayuda? Apreciaria enormemente.
Un abrazo extensivo a L. Haiduc.
Calduio Téllez
Apreciado Joaquin: Conozco algunos de los próceres que mencionas en el Blog. Cea Olivares, e algunos de los vampiritos. Haiduc tal vez. Estuve varias veces en México, una de ellas convidado por Jacobo Gomez Lara, ya fallecido, otras por Victor Castaño. bueno..
mi pregunta es la siguiente como buen novato. Estudiando el complexo bis-dietilditiocarbamato de Cu(II) encontré transferencia de carga entre:
Donador ACEPTOR ENERGIA
85 LP(2) S8 LP* Cu1 31.14
esto está muy padre y puedo visualizar el orbital 85
en Gaussian View y en el Chemcraft. Como visualizo el orbital 83? Si es antiligante deve estar arriba de l LUMO y no abajo del HOMO.
Puede usted ayudarme con mi duda? Apreciaria enormemente.
Un abrazo extensivo a L. Haiduc.
(texto corrigido!)
Claudio Téllez
Apreciado Claudio
Me da gusto saber que a través del blog has viajado mentalmente a México nuevamente, por supuesto que pasaré tus saludos al Dr. Haiduc en breve.
En cuanto a tu pregunta: Creo que hay un error de interpretación; según entiendo las lineas que me pones corresponden al análisis de perturbaciones a segundo orden (el cual se obtiene por default o con la opción E2PERT según el programa que uses). En esta sección el número 85 NO hace referencia al orbital, sino a la interacción (la interacción perturbativa número 85 que pasó el umbral -threshold- de 0.5 kCal/mol para ser impreso). Si estás graficando el ORBITAL NATURAL 85 quien sabe que sea lo que estás viendo, pero no corresponde a la interacción que describes. Lo que podrías hacer (asumiendo que etiendo tu intención) es buscar a que NBO corresponde el LP(2) del S8, es decir buscar su número en donde se encuentra la descripción de todos los NBO’s (ve mi post sobre el análisis de resultados NBO), graficarlo y hacer lo mismo para el LP* de Cu1. Molekel y Gaussview pueden presentar superficies simultaneamente (OJO! no es lo mismo que sumar las superficies, esto es completamente distinto). Así graficadas puedes analizar como se están traslapando ambos NBOs. Desconozco el resto de tu archivo pero me imagino que esta interacción es importante en la formación de tu complejo (¿tiolato?),revisa que no tengas interacciones perturbativas con mayores energías que pudieran ser más importantes que esta.
Espero que te sea de ayuda. No dudes en contactarme por este medio o por mail para cualquier otra duda.
Saludos hasta allá
I’m post graduated chemsitry student and I’m trying to learn gaussian, I use gaussian 03 for win. I am studying on TD DFT study a dinitrile compound. After 6 days processing, because of electricity
cutting, my gaussian job is terminated and computer is closed. I want to resumption my job but I dont know how to do it. is it possible? or should I open
latest out file and restart process.
I tried “Geom=check guess=read” functions but I thought I didnt manage
My input file;
%chk=C:/Users/TURKR/Desktop/tddft50trs.chk
%mem=100MW
%nproc=1
#p td=(singlets,nstates=50) b3lyp/6-311+g(d) guess=save geom=connectivity
pop=full gfinput
final part of my out file;
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX: IP= 1 NPass= 3 NMax= 34.
CISAX will form 34 AO SS matrices at one time.
NMat= 34 NSing= 34.
NMat= 34 NSing= 34.
NMat= 32 NSing= 32.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX: IP= 1 NPass= 3 NMax= 34.
CISAX will form 34 AO SS matrices at one time.
NMat= 34 NSing= 34.
NMat= 34 NSing= 34.
NMat= 32 NSing= 32.
Root 1 has converged.
Root 2 has converged.
Root 3 has converged.
Root 4 has converged.
Root 5 has converged.
Root 6 has converged.
Root 7 has converged.
Root 8 has converged.
Root 9 has converged.
Root 10 has converged.
Root 11 has converged.
Root 12 has converged.
Root 13 not converged, maximum delta is 0.001146383849072
Root 14 not converged, maximum delta is 0.001076901825947
Root 15 has converged.
Root 16 has converged.
Root 17 has converged.
Root 18 has converged.
Root 19 has converged.
Root 20 has converged.
Root 21 has converged.
Root 22 has converged.
Root 23 has converged.
Root 24 has converged.
Root 25 has converged.
Root 26 has converged.
Root 27 has converged.
Root 28 not converged, maximum delta is 0.001298203075257
Root 29 has converged.
Root 30 has converged.
Root 31 has converged.
Root 32 not converged, maximum delta is 0.001487599984012
Root 33 not converged, maximum delta is 0.002130187054738
Root 34 not converged, maximum delta is 0.009687228977448
Root 35 not converged, maximum delta is 0.014708674044360
Root 36 has converged.
Root 37 not converged, maximum delta is 0.001056178564194
Root 38 not converged, maximum delta is 0.001307476560090
Root 39 not converged, maximum delta is 0.001281081401660
Root 40 has converged.
Root 41 not converged, maximum delta is 0.001061548241041
Root 42 not converged, maximum delta is 0.001573623655971
Root 43 has converged.
Root 44 has converged.
Root 45 has converged.
Root 46 has converged.
Root 47 not converged, maximum delta is 0.001567486009780
Root 48 not converged, maximum delta is 0.005544365303514
Root 49 not converged, maximum delta is 0.015140942251070
Root 50 not converged, maximum delta is 0.016487405708613
Excitation Energies [eV] at current iteration:
Root 1 : 3.097975009587616 Change is -0.000001105943253
Root 2 : 3.583514852565895 Change is -0.000007448515000
Root 3 : 3.617241820182027 Change is -0.000021813870515
Root 4 : 3.991502586858061 Change is -0.000023174927667
Root 5 : 4.138876536189863 Change is -0.000007145969739
Root 6 : 4.211283026540273 Change is -0.000007302963600
Root 7 : 4.270327709174773 Change is -0.000012531626331
Root 8 : 4.368328284874000 Change is -0.000005937553695
Root 9 : 4.445535867701103 Change is -0.000032779976579
Root 10 : 4.540640063389401 Change is -0.000002225280201
Root 11 : 4.656404960218242 Change is -0.000023209382109
Root 12 : 4.714797549284374 Change is -0.000019602079988
Root 13 : 4.754991510970658 Change is -0.00001
Hi Turker
I’m catching up with my email and my blog comments today and I saw you already asked the CCL community. To restart a job you have to use the keywords geom=check guess=check but also you have to allocate the checkpoint file in the default location. NOW! you have to also type %NoSave below the %chk=filename.chk line, that way it will read it and wont write it! Be careful because if you don’t write that line then you will loose your checkpoint. I think the idea of copying and renaming it is a good one, try it but also use the %NoSave option I write about here.
Have a nice day!
Hello I’m surfing the net and I found you are know everything about gaussian03.
I have some problems and don’t know how can solve my problems.
In my thesis I being involved to find HOMO and LOMO wave function but I found HOMO only in wfn file
Please help to find LOMO wave function
Please
We are carrying out our calculations using G03W and GVW3.09.I have recently purchased G09W. But it is not possible now to submit the job from G03W to G09W.
Can you kindly guide me in this direction?
If I understand correctly you have GaussView installed along G03 and now you changed to G09 and as a consequence you cant link GaussView with G09. It has to do with the way Gaussview is installed so:
1) uninstall everything!!
2) install G09W
3) Install GaussView
I hope this helps
We are carrying out our calculations using G03W and GV3.09. I have recently purchased G09W. But it is not possible now to submit the job from GV3.09 to G09W.
Can you kindly guide me in this direction?
I think you need to reinstall GaussView, or just re-link it. I’m afraid you need to contact Gaussian Inc for this one or maybe your local IT guy can fix you up.
Have a nice day!
Hi Dr. Joaquin
I am trying to optimize metal complexes of phthalocyanine derivatives by using ONIOM method. I completed optimization of Zn-Phthalocyanine and calculated IR and UV. But I couldn’t optimize Cu-Pht and Co-Pht and I haven’t tried Ni-Pht yet. I met an error message in optimization of Cu-Pht. I think that For metal complex, ONIOM method is more complicated than general single layer methods.
Could you help me, please
Best Regards
Command line my out file:
#p opt=(tight,maxcycle=512,gdiis) freq=noraman oniom(b3lyp/6-31g:pm3mm) geom=connectivity int=ultrafine scf=maxcycle=512
Error message:
Warning! Cu atom 209 has 11 valence electrons but only 4 basis functions.This is less than a minimal basis set!
IExCor= 0 DFT=F Ex=HF Corr=None ExCW=0 ScaHFX= 1.000000ScaDFX= 1.000000 1.000000 1.000000 1.000000IRadAn= 5 IRanWt= -1 IRanGd= 0 ICorTp=0NAtoms= 209 NActive= 209 NUniq= 209 SFac= 1.00D+00 NAtFMM= 60 Big=TLeave Link 301 at Thu Nov 25 11:42:58 2010, MaxMem= 104857600 cpu: 0.0
(Enter C:\G03W\l401.exe)Simple Huckel Guess.NBasis= 620 NMin= 625 so simple Huckel guess is impossible.Error termination via Lnk1e in C:\G03W\l401.exe at Thu Nov 25 11:43:02 2010.Job cpu time: 0 days 0 hours 0 minutes 12.0 seconds.File lengths (MBytes): RWF= 52 Int= 0 D2E= 0 Chk= 1 Scr= 1
Hi Turker
Your basis set is just too small, try changing it for 6-311G(d,p) perhaps, or get one with a pseudopotential such as SDD or LANL2DZ. Check Gaussian’s manual for a wider variety of basis sets available for Cu and the other metallic ions you are studying.
Hi Dr. Joaquin
I am appriciated for your advices. You are right, I think semiemprical method isn’t suitable for Cu, Co and Ni atoms. LANL2DZ is good but too slow for the compound, so I have tried ONIOM method (DFT631G:PM3MM). I also tried LANL2DZ:PM3MM and only PM3MM. Results were same. One more interesting point; when I used molecular mechanics instead of PM3MM (for example ONIOM631G:UFF), I didn’t meet an error.
I am an amateur, I am studying and studying 🙂
thanks, happy new years
Dear Dr. Barroso:
I am trying to optimize a structure in which Zn is octahedrally coordinated with two citrate molecules. When I optimize this structure the bonding between Zn and the two citrate molecules dissapears. I know this complex is possible since I got the structure from the authors of a published crystallography paper. This is my input:
%chk=octacomplex.chk
%mem=70MW
%nprocshared=2
# opt rb3lyp/6-311+g(d,p) guess=(local,save) geom=connectivity
Title Card Required
-4 1
C 0.00000000 0.00000000 0.00000000
O 0.00000000 0.00000000 1.44957233
H 0.84063190 0.00000000 1.70043798
C -1.45846423 -0.22289327 -0.45896769
O -2.29631461 -0.61293155 0.41986215
O -1.72466070 -0.03349046 -1.65301136
C 0.88489416 -1.13110147 -0.54914157
H 0.78695471 -1.14345589 -1.51425010
H 1.80867833 -0.90423719 -0.35833379
O -0.35546025 -2.74486386 0.71897856
C 0.55671496 1.33042743 -0.51066400
H 1.47945129 1.40232696 -0.21933614
H 0.56235260 1.30216071 -1.48011487
C -0.17273032 2.58702456 -0.08041309
O -1.40015281 2.55436371 0.08757794
O 0.51874797 3.64592327 0.04721742
C -2.70592252 -3.08341975 4.08248350
O -2.70592252 -3.08341975 2.63291117
H -3.54664302 -3.08346576 2.38205127
C -1.24759139 -2.86049159 4.54141854
O -0.40960791 -2.47048821 3.66262135
O -0.98130632 -3.04984838 5.73545646
C -3.59086117 -1.95223738 4.63158667
H -3.49287722 -1.93996386 5.59673360
H -4.51473395 -2.17914766 4.44078464
C -3.34303653 -0.55032925 4.12272545
O -2.35059538 -0.33852100 3.36347229
O -4.12000079 0.34887973 4.48811035
C -3.26277058 -4.41381229 4.59311486
H -4.18537380 -4.48574672 4.30181964
H -3.26836372 -4.38562647 5.56260412
C -2.53332530 -5.67040942 4.16286395
O -1.30581420 -5.63770256 3.99486717
O -3.22475909 -6.72938904 4.03527183
Zn -1.35300556 -1.54173288 2.04124462
C 0.63698091 -2.53305560 -0.04027460
O 1.42338885 -3.44322386 -0.41004057
1 2 1.0 4 1.0 7 1.0 11 1.0
2 3 1.0 35 1.0
3
4 5 1.0 6 2.0
5 35 1.0
6
7 8 1.0 9 1.0 36 1.0
8
9
10 36 1.0 35 1.0
11 12 1.0 13 1.0 14 1.0
12
13
14 15 2.0 16 2.0
15
16
17 18 1.0 20 1.0 23 1.0 29 1.0
18 19 1.0 35 1.0
19
20 21 2.0 22 2.0
21 35 1.0
22
23 24 1.0 25 1.0 26 1.0
24
25
26 27 1.0 28 2.0
27 35 1.0
28
29 30 1.0 31 1.0 32 1.0
30
31
32 33 2.0 34 2.0
33
34
35
36 37 2.0
37
1 2 1.0 4 1.0 7 1.0 11 1.0
2 3 1.0 35 1.0
3
4 5 1.5 6 2.0
5 35 1.0
6
7 8 1.0 9 1.0 36 1.0
8
9
10 36 1.5 35 1.0
11 12 1.0 13 1.0 14 1.0
12
13
14 15 2.0 16 2.0
15
16
17 18 1.0 20 1.0 23 1.0 29 1.0
18 19 1.0 35 1.0
19
20 21 1.5 22 2.0
21 35 1.0
22
23 24 1.0 25 1.0 26 1.0
24
25
26 27 1.5 28 2.0
27 35 1.0
28
29 30 1.0 31 1.0 32 1.0
30
31
32 33 2.0 34 2.0
33
34
35
36 37 2.0
37
Any help would be greatly appreciated.
Best regards,
Gabriela Arias
Hola Gabriela!
What do you mean it disappears? If you mean you can no longer see it in a visualizer such as GaussView then you have nothing to worry about, sometimes visualizers have different thresholds for painting a bond or not.
Your initial bond lengths (Zn-O) where taken from the crystallographic data, right? How much longer are those bonds after optimization?
Take a look on the bond indexes to assert that the bond is lost or if it is just an artifact from the visualizer.
I need more information to help you out but I hope this sets you on the right path.
Happy holidays!
Dear Dr. Barroso:
Thanks for your quick response. In the optimized structure, Zn is no longer coordinated with the citrate molecules, so it is not a visualization problem. The distance between Zn and the binding sites in the optimized structure is about 0.10 longer than in the crystallographic data.
Happy Christmas,
Gabriela Arias
Hola de nuevo Gabriela!
Wouldn’t this mean your bonding is mainly ionic in nature? I mean, if the distances are so similar after the opt procedure and the bond order is so low then there is probably no covalent bond there to begin with! Try running a bond order calculation (I suggest NBO with BNDIDX keyword for a Wiberg analysis)on the initial and final structures just to make sure. Also employing Generalize Valence Bond theory could throw some light on to why are your bonds disappearing in the gas phase. Thanks for stopping by.
Have a great year!
Hi Dr. Barboso
I’m facing some troubles trying to visualize NBOs. I’ve optimized some molecules using g09 and then I made a single point calculation for each one and added some keywords for a NBO analysis. The problem is: when I try to visualize the NBO orbitals using GV 5.0, only the the usual MOs show up.
The CMO, NRT and other default analysis performed by NBO were in the .log file as they should be. The only problem is visualizing the orbitals.
I used the following keywords:
——-
#T M062X/cc-pvtz sp Pop=(NBORead,savenbos)
…
$NBO NRT CMO plot file=job $END
——-
I appreciate any help you can provide me.
With regards,
Italo
Hi Italo
I assume you already read my post (in this same blog) about how to visualize NBO’s with GaussView. Maybe there is something about the 5.0 version. I don’t know if I have it available but if I do I will run some tests and will let you know.
In the mean time I hope you have a nice day
Dear Italo,
Did you get solution for your error?
Kindly share with me, I am having the same kind of error.
Many thanks,
Bijan
I am trying to visualize the HOMO and LUMO of a compound on MOLDEN of which I ran gaussian calculations on. I read an article online that says that for MOLDEN to find all required information in the output file, I need to run a single point calculation at the desired level of theory including the pop=full and gfinput keywords. I entered the above in my input file and the following is my input:
%chk=Chlorin-1_b3lyp_go
%mem=23000mb
%nprocshared=8
# b3lyp/6-31G(d) sp scf=tight pop=full gfinput maxdisk=50gb
b3lyp geometry optimization using Gaussian
Chlorin-sp-1
The mem value in the run script is 23500mb and maxdisk=50gb. However, when I run this job, an error message comes up which says: ‘galloc: cannot allocate memory.’ I tried increasing the memory on the run script to 24000mb but I keep getting the message ‘Requested node configuration unavailable.’ I do not know what the problem is. Is there a problem with my input file or what else do you think I should do. I am sorry if I did not do a good job expalaining the problem; I am new to this. Thanks a lot for your anticipated response.
Hi Tijesunimi!
Sorry for the delay in my response. The error is very simple, you are trying to allocate A LOT of memory in the %mem part of your link0!
I don’t know what kind of processor are you using but if you already have the optimized structure and you only want to plot the orbitals then run this job (at the SAME level of theory than your opt) but get rid of those options (I mean %mem and maxdisk) let the program use the defaults. Also you can drop the %nprocshared=8 Gaussian will never use them all efficiently anyway. I hope this helps, please try it and PLEASE LET ME KNOW how it went ok?
I see you are working in photosynthesis, interesting because I’ve been thinking on doing that myself 🙂 If your chlorin molecule is too big then perhaps you can run it with %mem=300MB or something. Asking for 23GB just wont work 🙂
Take care and thanks for reading!
Joaquin, primero que todo, mis mas sinceras y grandes felicitaciones por tu blog, y tu esmero por ayudar a lagente a tra ves de el.
Segundo quisiera saber si existe alguna manera usando Gaussian 03 de hacer el calculo de un “ground state” pero “quitando” el lone pair de el nitrogeno, o los dos lone pairs del oxigeno y asi obtener un nuevo “ground state” donde estos lone pairs no son tenidos en cuenta y ver la diferencia entre estructuras….mil gracias por cualquier ayuda que me puedas dar
Hola Henry
Mil perdones por la tardanza de mi respuesta. Muchas gracias por tus palabras acerca del blog, espero que pueda seguir a la altura por mucho tiempo (contestar prontamente podría ser una buena manera de empezar).
Vi que pusiste esta misma pregunta en CCL pero ya no vi que respuestas obtuviste.
Sugiero que le entres al esquema NBODel, este método elimina interacciones particulares (elementos de la Matriz de Fock) y calcula la nueva función de onda con las condiciones dadas.
la forma de usarlo sería poner en el route section pop=NBODel y al final de la molécula
$NBO $END
$DEL
delete n orbitals
i
j
k
$END
donde n es el numero de orbitales a eliminar e ijk son las etiquetas NBO de cada uno. Corre un calculo NBO inicial (pop=NBO) luego localiza los pares libre etiquetados como LP, copia los números y ahora sí corres el cálculo con NBODel. Busca en este mismo blog el post sobre como interpretar NBO para más ayuda.
Espero que te haya servido Henry! Gracias por tus comentarios y por leer este blog!
Cheers!
Send me the details of calculating first order hyperpolarizability from G03output usning freq=raman keyword.
Thanking you
Dear Dr. Barroso,
thank you for your illuminative blog!
Perhaps you have a quick answer to my question:
I managed to display the HOMO of a molecule in GaussView 5.04 – but I can’t do two things:
– change the orbital colors to anything other than red and green;
– export the surface as an image file (preferably a vector image, of course)
Do you know how to do this?
Thanks a million!
David
Hi David!
Thank you very much for your kind words. About your questions: 1) I don’t think you can change the colors, I haven’t tried. Probably you need to do it externally with photoshop or something.
2) you can save it as a picture file (gif, tiff, jpg, etc.) but not as a vector file.
GaussView is not the best solution for this sort of high performance rendering you are trying to achieve. Let me suggest VMD, I believe that could be more up your alley 😉
Thanks for reading!
Hi,
I was actually looking for some hints on the internet about changing the orbital colours in Gaussview and ended up with this blog. Though, I was not able to get the answer here, I had a play around Gaussview and finally found the way to change orbital colours. Go to File—Preferences—Colors—Surface Colors. You should be able to change the colors there. I hope this is useful.
Thanks
Shanthi
Hi Dr. Barroso,
Sorry, I would like to ask a simple question. I ran a basic job as follow from a Gaussian textbook.
#T RHF/6-31G(d) Pop=Full Test
Formaldehyde Single Point
0 1
C 0.0 0.0 0.0
O 0.0 1.22 0.0
H 0.94 -0.54 0.0
H -0.94 -0.54 0.0
But, an error message showed up saying,
galloc: could not allocate memory.: Cannot allocate memory
galloc: could not allocate memory.
I saw one of your previous reply, but I don’t have %mem in the input.
Could you please tell me where has gone wrong?
Thanks a lot,
jovis
Hi Jovis
Include the %mem line at the top of your input file and then set it accordingly to your available computation resources, say
%mem=128MW
#T RHF…
Check this link for more info http://gaussian.com/g_tech/g_ur/m_eff.htm
I hope this helps and thanks for reading!
Dear Dr. Barroso, I have a problems with my input file for G03W
here is my input
%mem=100MB
%chk=dai.chk
#n pop=(NBO,savenbo)density=current rmp2/aug-cc-pvdz scf=tight
NBO-HF at MP2 wfn
0 1
S,0.,0.,0.5909815593
C,0.,0.,-1.0398818668
H,0.,0.9356771308,-1.6155519262
H,0.,-0.9356771308,-1.6155519262
but when I was running programme I got the error with
Sorting of NBOs: 48 30 47 45 46 64 56 29 28 44
Sorting of NBOs: 43 42 27 26 25 41 63 55 40 62
Sorting of NBOs: 54 24 23 22 53 61 60 52 59 51
Sorting of NBOs: 21 39 38 37 36 58 50 20 57 49
Sorting of NBOs: 35 34 33 19 18 17 32 16 31 15
Sorting of NBOs: 14 13 65 67 68 66 12 2 1 11
Sorting of NBOs: 3 4 6 10 9 7 8 5
NBStor is confused about NOcc.
Error termination via Lnk1e in C:\G03W\l607.exe at Mon Jul 18 07:45:50 2011
Please tell me why? And could you help me to fix it?
Many thanks to you!
Regards!
Dai.
Hello Dai
I get this question a lot and I still haven’t found a right solution, I need to work on it a little bit more. This has to do with a warning that NBO sends you if the occupation of each orbital doesn’t follow the aufbau principle, i.e., the energy and the population don’t follow the same order. Try the following:
1) Make sure you are working with a minimum of the PES by performing a vibrational frequency analysis (you probably already did but its better to be sure)
2) Change your basis set. Although NBO is not so sensitive to change in basis set there could be a reorthonormalization problem with the particular one you are using.
Hope this works. Thanks for reading
Hi Dr.!
I guess I’ve just found that my problem depends on my program,, that is G03 doesn’t include NBO5.G so I can’t perform NBO analysis with MP2 wavefunction (keyword: density=current). DO you agree with me?
Regards!
Dai
I think its a good idea, Dai! Calculate the mp2 wavefunction and then perform the nbo on it, it might work. Try it and please let me know how it goes ok?
Have a nice day and keep coming back!
It’s ok but I used G09, not G03. You can’t perform NBO on G03, you only perform NBO analysis from HF-wfn. I’ve just know that.
Thank for you kind, Dr.!
Regards!
p.s: you dont mind if I will have questions for you in the future?
Of course I don’t mind!!! I’m glad I can be of help and I like to interact with other comp.chemists from around the world.
Have a nice day!
Dear Dr. Barroso,
I have performed two calculations using Nosymm keyword with G09. I know that in this way all calculation are performed in the Input Orientation. However in one calculations I obteined a Z-Matrix orientation, while for the other I have an Input Orientation.
I ckecked the output, in both cases the symmetry was turned off. The only difference between the calculations is the size of the system, for the system with 61 atoms I have Z-Matrix Orientation , for smaller system I have Input Orientation (In both cases the command used where the same).
Do you where these differences come from. Is there some way to control it?
Thanks in advance.
Hello Fernanda
First I appologize for the lateness of my response. As far as I know there is no way to control this since the program will re-orientate the molecule to where its best suited even if the symmetry is turned off. the symmetry option is more related to the symmetry of the wavefunction, which is a consequence of the symetry of the molecule; a molecule could be symetric and the wavefunction can be made non-symetric with this option. However orientations are a bit harder to control.
Hope this helps!
dear Dr Barroso.
please help me to know how i can calculate bond order in NBO (G03).
I know how can i calculate wiberg bond index, are these two parameters equal?
Hello Mohamad!
They are definitely not the same thing since their definitions are diferent. If you run a full population analysis with nbo pop=(full,nboread) and then you use bndidx in the line:
$NBO BNDIDX NAO $END
you will find both bond orders (Wiberg bond index and Bond Order in the NAO basis)
I hope this helps
Have a nice day
what should be the input file for the optimization of metal complex of phthalocyanine in g03w rather than using oniom calculations?
Hi, Dr. I research in computational chemistry with use Gaussian 03/gaussview 05 software . so , I need a more tutorial or learning in field .please guide me for useful link, literature, book and journal particular commends in Gaussian program. thanks
Hello Fardin!
Go to http://www.gaussian.com and look for tutorials under the Technical Support page. Also try to find a book called “Exploring Chemistry with Electronic Structure Methods” by Frisch and Frisch. This book is very illuminating about how to work with Gaussian
I hope this helps
hi Dr.
ok and thanks.
Hello Fardin!
Go to http://www.gaussian.com and look for tutorials under the Technical Support page. Also try to find a book called “Exploring Chemistry with Electronic Structure Methods” by Frisch and Frisch. This book is very illuminating about how to work with Gaussian
I hope this helps
Greetings! Thanks for following my blog “Word Play – A Place For Pleasurable Procrastination”! I’ll try my best to keep you entertained. 🙂
Thanks, Chris! I hope my blog seems entertaining to you too 🙂
Sir, I try to optimised a cluster of TiO2 which is designed by PBC. But after 6 minute of running in Gaussian 09 windows version, an error via Lnk1e in C:\G09W\l401.exe is appeared. I m very much new in this field. So, I aspect your kind help…
Did you find a solution to your problem?
Hi Dr. Barroso,
I just stumpled upon your great blog! My question is closely realted to Debajyotis. I am trying to set up a PBC calculation on a measured crystal structure, but I always end up with the same error in the MO guess:
Dimensions of /B2/ exceeded in AtP2B2.
Error termination via Lnk1e in c:\program files\G09W\l401.exe
I am using #p pbepbe/3-21g/auto and I really dont have a clue whats the problem. He has almost 200 atoms in the unit cell, maybe its to large? But he doesnt complain about the memory….
Thanks a lot for your kind assistance!
Hello, I’m interested in the same gaussian error, if you already found a solution please let me know.
Hi, I’m having the same trouble when optimizing a crystal structure using gaussian09. If you have already found an answer, please let me know, thanks!
PBC in Gaussian is a mystery. You should try using CRYSTAL. In this case it means that the number of atoms as defined by the unit cell is exceeded within the k-mesh, that is that two adjacent units of the mesh are actually overlapping. The way to fix it is re-defining the k-mesh. i’m afraid that’s the best I can do.
Hello,
This is my first time posting. I am new to computational chemistry and Gaussian. I’m using Gaussian09.
My general concern is that when something goes wrong – as it seems to every single day – I don’t know where to look for resolution of the various error messages. I usually Google the error and look to CCL.net. But, some errors are so obscure that I don’t get any relevant search results. Even though I am in grad school, I’m the only person in my lab-group who does computational chemistry. One other person does some work with ADF.
My specific concerns are with organometallic complexes – finding reaction paths for now. I have been having a lot of trouble optimizing a lanthanide complex Ce(IV) Ce(III). I finally got an successful SCF calculation using SCF=QC. But, using guess=read, I keep getting an error that says Gaussian doesn’t see the basis set info in my checkpoint file. But, I’m supplying it in my input file.
Hello! Welcome to my blog and definitely welcome to the obscure world of Gaussian error messages. You will be surprised when you begin to solve them on your own in no time but it will take patience and many trials.
Your address suggests you attend the UCLA so you must have a valid Gaussian license, don’t be shy and directly ask their help desk they are very helpful and prompt.
If you use the option guess=read then the program gets the coefficients from the check file but tries to get the corresppnding basis set from it too. Try using the GEN option to override the search for the BS in the chk file along with the previous keyword. Not sure if itll work but is worth the try!
Hope this helps. Have a nice day
Thanks for the suggestion and the welcome! Since I am using pseudo potentials, I used gen and specified my basis sets in the input file. So, it was odd that Gaussian didn’t read them.
I submitted my query to Gaussian tech support.
Are there any other resources you can recommend?
I frequently come across old CCL.net postings about error messages. The Aeleen Frisch book was good for starters, but there are a lot of topics not covered.
Thanks again.
I would suggest to subscribe to the CCL, if you have an email with an unlimited amount of space like yahoo or gmail (practically) you can store all the threads. This way you can create your own database; sometimes when you search their archives you find the same post various times but answers come harder to find. This is what I do, anyway.
Why do you need the guess to be read from the chk file? is your system too big?
Have a nice day and please share the answer you get from Gaussian Inc if you can.
Hi. I am with same problem of the Anja. If the pabukha user found the answer, could post here? I am working with glycine crystal. Thanks.
Hi, I got a response from Gaussian. BTW – I checked the box to alert me when there’s a response but I didn’t get one.
Here’s the response.
(The first response was that there may be a problem with the location or corruption of the chk file I was referring the input to. But I saw that that wasn’t it. Here is the follow up e-mail from Gaussian.)
There are a few things to try that would save a good amount of computational time over simply switching to the QC SCF procedure. If you are running into SCF convergence difficulties, here are a few general guidelines to follow that can be helpful to you:
1) “SCF=NoVarAcc”. This one requests the use of full integral accuracy from the first SCF cycle. The default is to use lower accuracy in the first few cycles in interest of speed and then switch to full accuracy as the SCF approaches convergence. In some difficult cases, the initial lower integral accuracy is not good enough and the SCF procedure might do a lot of useless iterations. This option never hurts, it would make the first few cycles a little slower but the overall time may be lower if it takes fewer cycles to get close to convergence.
2) Increase integral accuracy. The default thresholds are set for an integral accuracy of 10^-10. In some cases when using very diffuse functions, one might need to use tighter thresholds since there will be very small quantities being computed. The integral accuracy is controlled with “Integral=(Acc2E=N)” in G09 or by “IOp(5/87=N)” in G03. The default is N=10. One may try increasing the accuracy by one or two orders of magnitude (i.e. N=11 or N=12).
3) Tighten the grid for numerical integration in DFT (This does not apply to non-DFT jobs such as HF, MP2, etc.). The default (“Fine”) grid for numerical integration is adequate for most purposes but in some cases the small numerical noise might make difficult SCF convergence and the SCF procedure takes too many iterations when it is close to convergence. If the problem is with large oscillations in early SCF cycles, this is likely not going to help. One can tighten the grid with “Integral=(Grid=XXXXX)”, where XXXXX is the value of the integration grid (see the documentation of the “Integral” keyword for more details). In most cases, tightening the grid to the “Ultrafine” grid is sufficient.
4) “SCF=NoIncFock”. Turning off the incremental Fock formation can help in cases with very diffuse functions. This increases the computational cost of the SCF calculation as it requests the recomputation of the full Fock matrix at each SCF cycle, instead of only doing an incremental Fock formation (default). The increase in cost with respect to the default SCF procedure is lower than doing QC SCF (QC SCF is not available for SCRF calculations).
5) “SCF=Fermi”. This one does a combination of things including Fermi broadening, damping and level shifting, deciding dynamically the most suitable strategy. It can be helpful in cases with oscillations in early cycles of the SCF procedure.
6) “SCF=XQC” (Note that QC SCF is not available for SCRF calculations). Turns on the QC SCF procedure after the maximum number of regular SCF cycles have been done without convergence (if the SCF is converged in fewer cycles with the regular SCF procedure, then the QC step is skipped). The maximum number of regular SCF cycles by default is 128, and this can be controlled by “SCF=(MaxConventionalCyc=N)”. The QC procedure is much slower than the regular SCF procedure but it is very robust and almost guaranteed to converge to an SCF solution. It also has higher chance to converge to the nearest SCF solution (which could be an unstable solution) than the regular SCF procedure, so it is strongly recommended to check for instabilities. The “SCF=XQC” option thus provides a fall-back solution when the regular SCF procedure fails to converge.
Hola Dr. Barroso me encantaría realizar mi verano de investigación con usted
Por supuesto! Ponte en contacto conmigo para ponernos de acuerdo. Gracias por tu interés.
Buenas Tardes!!! Dr. Joaquin
A mi también me gustaría trabajar con usted en el verano próximo.
Y estuve revisando los programas de veranos de investigación y el que más me llamo la atención fue el de Delfín y vi que usted estuvo registrado el del año pasado, y quisiera saber si para este año también estará?
Bueno saludos desde León, Gto! y de antemano muchas gracias.
Claro Fernando, espero acordarme a tiempo de darme de alta en el programa delfín de este verano. Muchas gracias por tu interés en nuestra investigación. Te esperamos!
Cordial saludo.
Dr Barroso le escribo porque tengo una duda, tratando de hacer calculos con “counterpoise correction” no se a que distancia debo colocar los reactantes, ya que el valor obtenido de “Counterpoise: corrected energy” disminuye a medida que la distancia entre los reactivos aumenta; si tiene alguna sugerencia le agradezco su puntual colaboracion.
Emiliano,
me parece que lo más pertinente en este momento es que busques la literatura correspondiente al método Counterpoise para calcular el BSSE (Basis Set Superposition Error).
Cuando quieres calcular la energía de una ‘supermolécula’ (más de un fragmento unidos no covalentemente) se produce un artificio matemático por el cual la densidad electrónica de una molécula es transferida a la base de la otra (de ahí el nombre de superposición de funciones de base), con lo que la energía queda sobreestimada. El método Counterpoise corrige este efecto. Ese traslapamiento será mayor a medida que estén más cerca los fragmentos, como ya lo viste. La pregunta no es a que distancia los colocas, sino que es lo que quieres calcular y porque la corrección por counterpoise es importante. Otra manera de decirlo es que el BSSE (calculado por counterpoise) no es una cantidad/propiedad física sino un artificio ocasionado por la manera en que son construidas las funciones de onda.
Espero que te sirva. Gracias por leer el blog!
#Teaching Hello, could you please explain what is the difference between MRMP2 method introduced by Hirao and CASPT2. Thanks in advance!
HI Dr, thank you for this useful blog. I have a question please: I want to evaluate the solvent effect on fluorine molecule, I tried variable solvent but the energy is the same energy of F2 in a gase phase. Can you help me please to find the right key words to run the job on gaussian09.
Can I modify the dielectric constant for a solvent? when I use the word dielectric i obtain a message error saying that there is a syntax error.
thank you very much for your help.
Hello Wiem,
Thank you very much for your kind words about my blog!
What keywords are you using now? Perhaps you are using the Onsager method (which for a non polar molecule as F2 is equal to a gas phase calculation). Try using the PCM or CPCM models with available solvents in gaussian. You can change the dielectric constant but I don’t recommend it; if you want to define your own solvent you need data as the molar volume of it among other things so try to stick to what is already well defined.
I hope this helps!
Hello Dr Barrosso, thank you for your quick response, the problem is that i tried every key word existing in the manual of G09, but i can’t find a notable difference with the gas phase or between solvents! (notice: i did try just few of them) Here are some energies: E(phase gase)= -195.9816245 u.a , using CPCM E(water)=-195.9816792 u.a E(CCl4)= -195.981655u.a E( Ethanol)=-195.9816777 u.a ,, using PCM E(water) = -195.9816786 u.a , using SMD E(water)=-195.9781132 u.a , E(ethanol) = -195.9806675 u.a
Do you think that it is useless to modulate the solvent effect on this molecule?
Thank you very much.
Hello sir,
I’m running gaussian calculations for isomers in different solvents using different methods and different basis sets..for example with B3LYP and MP2 with 2 basis sets each..
when I’m performing Optimization getting error as error in internal coordinate system with some methods..I can not change the geometry bcoz i need to use same geometry for all methods..i tried by giving different keywords likeGuess=AM1, SCF=QC,
OPT=mod redundant .OPT Guess=Mix etc…. but for some it didn’t work with any of these…pls help me regarding this..
One more problem is during freq calculations some times its giving so many errors like could not allocate memory , use 2340000 instead of 11196670 …. i tried by changing mem% and MaxDisk…
Regarding your first question, there is little that can be done if you can’t change the geometry. Have you considered making partial optimizations? That might work.
About your second question, try using rwf partition files. Google “efficiency considerations site:gaussian.com” for a link to a useful tutorial on the matter
Have a nice day!
I like Doctor Joaquin investigate computational chemistry issues with you that would be possible … and good conduct publications
Sure! Send me an email and we’ll talk about it.
Have a nice day!
Hello Doctor
I am working on my thesis and I need to obtain uv-vis spectra of some quinoline carboxaldehyde molecules in gaussian 03. please if you can help me necessary procedure for performing uv-vis because I am totally clueless.
Thank alot
Hello Dr. Joaquin Barroso-Flores
I am new to Gaussian and have access to Gaussian 9.0 and Gaussview. I want to generate potential energy surface for a molecule (using scan option around one torison angle). Could you please let me know if there is any link available providing step wise directions to submit and analyse job ?
Thanks & Regards
Rajni
Hello, Rajni!
There is a fine link right here on this blog, please check it out.
Hope this helps
Hello Dr. Joaquin Barroso-Flores,
I am trying to draw a simple crystal structure of TiO2 using Gaussian03 and GaussView 5. I need to use PBC option, and when i press calculation setup, Gaussian03 starts an error appears :
Dimensions of /B2/ exceeded in AtP2B2.
Error termination via Lnk1e in /home/gaussian/g09/l401.exe at
copy-paste from .LOG file.
Can you please help me?
Tengo la misma duda que Djohra apreciado Doctor.
¿Podría ayudarnos a ambos por favor?
Djohra: I guess i have the keyword to avoid these error, just write guess=indo on the keywords list. It works for me
Dear Professor Joaquin please correct me if i am wrong.
Sir
I am from India.
I have been trying optimize the structure of cyclopropenyl carbene(triplet) using MP2/6-31G(d), I am sending the input file. But it says that there is a #2070 error after the optimization. Can you please check whether my level of theory and basis set is sufficient for this system, or suggest a better alternative.
%mem=6MW
%nproc=1
%chk=CPTRIP2.chk
# opt ump2/6-31g(d) nosymm geom=connectivity
Title Card Required
0 3
C
C 1 1.299101
C 2 1.463395 1 63.649311
H 1 1.079741 2 149.639600 3 -170.548274
H 2 1.079777 1 149.623659 3 170.566591
1 2 2.0 3 1.0 4 1.0
2 3 1.0 5 1.0
3
4
5
Hola doctor Barroso, buenos días. Tengo una pregunta respecto al keyword EPT en Gaussian 09 para calcular los potenciales de ionización. Mi input es el siguiente:
#T b3lyp/6-31+g(3d) EPT
Cuando empiezo a correr el programa, en el output aparece el siguiente error:
Internal consistency failure #1 in ROv08.
Usted podría orientarme para saber cuál es mi error y cómo corregir este problema?
Saludos y buenas tardes.
atte.
Daniel
Hola Daniel,
Intenta usar otra geometría (si no está optimizada la que estás usando). Por otro lado, si estás usando la última revisión, la D.01, intenta correr la anterior si la tienen disponible. A veces este error sale por algún error interno en el código. Repórtalo a Gaussian para que tengan más información y se corriga en la siguiente revisión.
Saludos y espero que te sirva!
please doctor how to encapsule guest molecule manualy with gaussview in cyclodextrrin host; pleasewith details.
Hi Sir;
how can i creat .46file (nbo) for visualize orbitals…by jmol ?
for example my input is:
%chk=D:\hosseini\NBOGly.chk
#P HF/6-31G(d) pop=(nboread,savenbo) geom=connectivity
Title Card Required
0 1
N -1.96650000 0.01218900 -0.00005900
H -2.02138500 0.61821600 0.81086400
C -0.72608900 -0.72384900 0.00017900
C 0.54190200 0.11216500 -0.00000200
H -0.68266900 -1.38671000 0.87122100
H -0.68273900 -1.38736800 -0.87030100
O 0.58818600 1.31491000 0.00004900
O 1.64526900 -0.66519300 -0.00011900
H 2.41145000 -0.07420000 -0.00009500
H -2.02167800 0.61710900 -0.81178500
1 2 1.0 3 1.0 10 1.0
2
3 4 1.0 5 1.0 6 1.0
4 7 2.0 8 1.0
5
6
7
8 9 1.0
9
10
$NBO plot $END
but exist .31…41 files in output except .46 file,
Try including the keyword ARCHIVE=FILENAME after PLOT in the NBO section. Let me know if it helps.
Have a nice day
Hi sir, thanks a lot for your attention and your answer,but .46file could not created
Mehran
Hi sir, Hi freinds
i see how i can creat .46 file (with .46 file we can visualize orbital interaction … by NBO5 or 6 softwear) and i want to share it with you.
for creating .46 we must open .47 file by notepad++ and change the secound line :
$ nbo $end to $nbo plot file=addres $end then save it.
now run gennbo5 and open .47 file , gennbo5 produce .46 file and …
Thanks for the tip!
Hello Dear Dr. Joaquin
How can I calculate NICS(π), NICS(σ), NICS(π)zz and MO- NICS with Gaussian. Actually, I need to calculate NICS for pi and Sigma separately. Please help me for solving this problem.
thank you
Hi, I would like also to get some kind of tutorial. I have been looking at the literature and programs like Aroma (from Prof. Stanger) but I’m not yet capable to easily calculate NICS values by myself looking just at a log/out Gaussian file.
Thanks!
Hello sir,
I’m running Gaussian 09 TD-DFT calculations for truncated part of my selected Covalent organic framework. With .#p opt td=(nstates=10,root=1) rb3lyp/6-31+g(d,p)geom=connectivity scf =(maxcycle=700,conver=6) My current running job is Td optimization but it could not get termination for many days although in ground state it is easily completed. ..
I brought little changes in bond length between COF and sensor molecule and ran the job again but there was no much of differences in all completed cycles while only shifting of double bond is occurred towards phenyl sides attached to ethylene (four Phenyl ring attached to C=C bond) with same order. Please I need your valuable suggestion for job completion.
hOur company offers supreme quality weight loss products. Take a look at our health contributing portal in case you want to improve your health. Our company offers supreme quality non prescription products. Take a look at our health contributing portal in case you want to improve your health. Our company provides herbal weight loss products. Look at our health contributing portal in case you want to feel healthier. Our company provides a wide variety of healthcare products. Look at our health contributing portal in case you want to feel better. Our site offers a wide variety of non prescription drugs. Take a look at our health website in case you want to look healthier with a help generic supplements. Our company provides a wide variety of non prescription products. Visit our health portal in case you want to feel better.
Our company provides supreme quality general health products. Take a look at our health contributing site in case you want to feel better. Our company provides herbal non prescription products. Look at our health contributing site in case you want to look healthier. Our company offers safe non prescription products. Visit our health contributing portal in case you want to feel healthier. Our company provides a wide variety of non prescription products. Take a look at our health site in case you want to strengthen your health with a help of health products. Our company provides herbal healthcare products. Visit our health contributing portal in case you want to look healthier.
Hi,
When I try to perform reaction parth calculations I got following error message.
Severe Error Message # 2070
My input file is as follows: If you could please help me to solve it. Thank you.
%chk=D:\Calculations Reaction parth\New folder\QST2CAL\Reactance_CAL_QST2nEW_LABLE_submitnew.chk
# opt=qst2 hf/6-31g(d) nosymm geom=connectivity
Title Card Required
0 2
C -1.91543797 7.29206211 -0.58482610
C -0.48352897 9.70276911 -0.72666610
C -0.50255697 7.24656511 -0.53358710
C -2.59895397 8.50147111 -0.67975410
C -1.88132197 9.70087111 -0.74237010
C 0.17440103 8.48460211 -0.62890910
H -3.68472497 8.48757611 -0.71605210
H 0.07216703 10.63019611 -0.79189110
H 3.58663903 9.65364611 0.08975490
C 4.25050803 8.87737111 -0.27080910
C 6.01194303 6.91676211 -1.19848010
C 5.60954903 9.14240211 -0.33292510
C 3.73616703 7.61943211 -0.66334810
C 4.64323703 6.64725511 -1.14524510
C 6.50071803 8.15109811 -0.78950810
H 4.28081003 5.71141211 -1.55658610
O -2.50273197 10.91303011 -0.83766710
H -3.46255197 10.77452411 -0.84014110
O -2.67168297 6.15444511 -0.52879710
H -2.10088997 5.38572111 -0.68051510
O 2.36354703 4.98190311 -0.18329210
H 3.26822403 5.21041511 0.08990490
C 0.27529703 6.05694411 -0.35073410
H -0.16867197 5.09742411 -0.10752110
C 1.66128203 6.14284811 -0.38051010
C 2.30716803 7.36856411 -0.58354910
O 1.54292503 8.52581111 -0.61954510
H 6.69755603 6.16411611 -1.58365610
O 6.08295503 10.35869011 0.05707090
H 7.04784903 10.35053911 -0.06468510
O 7.82768003 8.51956011 -0.80032610
H 8.36422003 7.79986611 -1.16482110
C 10.60247067 5.13437726 -0.01990632
C 9.48271481 4.34724963 -0.01993858
C 9.83081630 2.96815434 -0.01996448
C 11.21062623 2.72859474 -0.02000698
C 10.45636119 0.53838262 -0.02002535
C 11.57187534 1.36766920 -0.02004009
H 8.48431551 4.76773557 -0.01985678
H 12.58824334 0.99630883 -0.02002798
S 12.09930637 4.23303385 -0.01989543
S 8.94562388 1.47908566 -0.02003854
C 10.33719745 -0.88790671 -0.02002897
H 9.32217874 -1.27556669 -0.01985610
C 11.31664825 -1.84517802 -0.02017125
C 10.99450926 -3.29331368 -0.02020841
C 12.70920778 -1.53307343 -0.02024885
N 13.83303924 -1.22892949 -0.01939226
O 9.65718228 -3.53794348 -0.01972967
H 9.56940368 -4.50870525 -0.01985347
O 11.82183130 -4.17977476 -0.02063055
H 10.64030505 6.21596625 -0.01996381
1 3 1.5 4 1.5 19 1.0
2 5 1.5 6 1.5 8 1.0
3 6 1.5 23 1.5
4 5 1.5 7 1.0
5 17 1.0
6 27 1.0
7
8
9 10 1.0
10 12 1.5 13 1.5
11 14 1.5 15 1.5 28 1.0
12 15 1.5 29 1.0
13 14 1.5 26 1.0
14 16 1.0
15 31 1.0
16
17 18 1.0
18
19 20 1.0
20
21 22 1.0 25 1.0
22
23 24 1.0 25 1.5
24
25 26 1.5
26 27 1.0
27
28
29 30 1.0
30
31 32 1.0
32
33 34 2.0 41 1.0 52 1.0
34 35 1.5 39 1.0
35 36 1.5 42 1.0
36 38 1.5 41 1.0
37 38 1.5 42 1.0 43 1.5
38 40 1.0
39
40
41
42
43 44 1.0 45 2.0
44
45 46 1.0 47 1.5
46 49 1.0 51 2.0
47 48 3.0
48
49 50 1.0
50
51
52
P03 QST2 CAL
0 2
C 8.19728381 0.73296959 -0.09356053
C 6.48119981 2.94327659 0.18289447
C 6.80406481 0.50594559 0.01661147
C 8.72317181 2.02018559 -0.04628553
C 7.86572081 3.11789159 0.09853147
C 5.98124181 1.65033259 0.13519347
H 9.79802781 2.14929659 -0.13559953
H 5.81796581 3.79408459 0.27875047
H 2.50707081 2.33225959 1.04372647
C 1.92216781 1.51545359 0.63981947
C 0.35540881 -0.55272341 -0.34114253
C 0.53726481 1.62890859 0.67458147
C 2.55861781 0.36339059 0.13349347
C 1.73134381 -0.67794741 -0.37193853
C -0.29467719 0.58831759 0.17971247
H 2.16541781 -1.52896741 -0.88519053
O 8.33119181 4.39724359 0.14739247
H 9.29849881 4.38713859 0.07643847
O 9.08525281 -0.29372241 -0.23519753
H 8.60756881 -1.10752841 -0.45660753
O 4.26134981 -2.12961741 0.18957847
H 3.36714181 -2.06910741 0.56489947
C 6.18223481 -0.78178341 0.06984847
H 6.75112681 -1.70274041 0.14571047
C 4.80141581 -0.87371741 0.12057847
C 3.99941781 0.28271259 0.12239147
O 4.62257781 1.51332259 0.21709947
H -1.51095593 2.51918030 0.10218368
O -1.51095593 3.28074230 0.70078368
H -1.51095593 4.04230430 0.10218368
O 0.05288381 2.78098659 1.21907747
H -0.89856319 2.69942959 1.40267347
C -1.75227519 0.63643459 0.19127947
C -2.61018419 -0.43310641 0.37212047
C -3.97024419 -0.05930641 0.27772147
C -4.18094019 1.30458259 0.03368347
C -6.38774019 0.61330559 0.15723047
C -5.52864019 1.69297259 -0.03373653
H -2.25832219 -1.43168041 0.60117847
H -5.87551319 2.70223059 -0.21290653
S -2.65631719 2.14186559 -0.12513953
S -5.48066419 -0.89615241 0.42704147
C -7.81331519 0.53029659 0.17021747
H -8.22648719 -0.45892341 0.34645247
C -8.74796219 1.51881959 -0.00692053
C -10.20175719 1.23671659 0.04149947
C -8.40091619 2.88128059 -0.24761553
N -8.06655719 3.97995859 -0.44054153
O -10.48002419 -0.07459441 0.27353847
H -11.45270819 -0.13542841 0.28351447
O -11.06940119 2.07133959 -0.10735753
H -0.25471919 -1.33470941 -0.78328453
1 3 1.5 4 1.5 19 1.0
2 5 1.5 6 1.5 8 1.0
3 6 1.5 23 1.5
4 5 1.5 7 1.0
5 17 1.0
6 27 1.0
7
8
9 10 1.0
10 12 1.5 13 1.5
11 14 2.0 15 1.5 52 1.0
12 15 1.5 31 1.0
13 14 1.5 26 1.5
14 16 1.0
15 33 1.0
16
17 18 1.0
18
19 20 1.0
20
21 22 1.0 25 1.0
22
23 24 1.0 25 2.0
24
25 26 1.5
26 27 1.0
27
28 29 1.0
29 30 1.0
30
31 32 1.0
32
33 34 2.0 41 1.0
34 35 1.5 39 1.0
35 36 1.5 42 1.0
36 38 1.5 41 1.0
37 38 1.5 42 1.0 43 1.5
38 40 1.0
39
40
41
42
43 44 1.0 45 2.0
44
45 46 1.0 47 1.5
46 49 1.0 51 2.0
47 48 3.0
48
49 50 1.0
50
51
52
SCF=Tight
Dear Dr. Joaquin Barroso-Flores
I run a single point calculation with Gaussian09, revision A.02. My input file contains fragment-combined wavefunction as initial guess. Part of my input file is as follows:
When I run this input file, I receive an error (Error termination via Lnk1e in C:\G09W\l401.exe
). Could you please kindly guide me in this problem? I appreciate in advance for the time you devote to answer my question.
Best regards,
Mehdi
Chk=new.chk
%mem=700mb
%nprocshared=2
#p rcam-b3lyp/gen pseudo=read nosymm guess=cards
Please check this file to ensure validity
2 1
C 1.90992500 3.16721400 -0.50060600
C 0.54789600 2.74660600 -0.44522800
C 0.19538600 1.36319800 -0.49336500
C 1.25707900 0.46851800 -0.75010000
C 2.59624300 0.86019300 -0.71804900
C 2.93638500 2.21919900 -0.59908800
C -1.20085300 0.98099500 -0.29788200
C -1.66491700 -0.33704700 -0.09830500
C -1.05101200 -2.62529500 -0.15795600
C 0.02904800 -3.53159000 -0.41070300
C 1.40536500 -3.15465800 -0.60944800
C 1.87884200 -1.81648000 -0.80490800
C 3.22699900 -1.49410400 -0.76687600
C 3.60293700 -0.12763500 -0.78042900
C -0.47035900 3.75364400 -0.35159700
C -1.83057600 3.38527900 -0.31726600
C -2.20234800 1.99958600 -0.28581900
C -3.59026600 1.67436500 -0.23011600
C -3.99860200 0.33507800 -0.19172200
C -3.01935600 -0.67180300 -0.13135900
C -2.37847600 -3.02318500 -0.20008200
C -3.39102800 -2.03386800 -0.12896800
C 2.26924200 4.54648900 -0.43784800
C 4.30754300 2.59176200 -0.60208700
C 4.65641100 3.95365500 -0.53107300
C 3.63124900 4.93130100 -0.45735500
C 5.31259000 1.60285900 -0.67171400
C 4.95884500 0.23362700 -0.73382300
C 5.95028900 -0.77580700 -0.73003900
C 5.57480400 -2.14614800 -0.72175100
C 4.21039800 -2.50612800 -0.71503000
C 2.40959500 -4.17137600 -0.63363800
C 3.79869800 -3.84685100 -0.68822600
C -0.32664800 -4.91577500 -0.44955900
C 2.03733700 -5.54933900 -0.59725100
C 0.68804700 -5.91503900 -0.54820300
C -2.71940600 -4.38906800 -0.31263300
C -1.68935800 -5.33710800 -0.39764900
C -4.74500400 -2.40228400 -0.17770400
C -5.09348300 -3.77203100 -0.24082600
C -4.07947200 -4.76513900 -0.30740600
C -5.37102000 -0.03357400 -0.19102100
C -5.74064800 -1.39578600 -0.18152000
C -4.59566600 2.68650100 -0.23268800
C -5.96609200 2.33228500 -0.21534100
C -6.35873000 0.96926400 -0.20312000
C -2.83779700 4.39378800 -0.27848200
C -4.21587000 4.04367300 -0.24931000
C -0.10823800 5.13257200 -0.32917700
C -1.11058900 6.12683400 -0.28094300
C -2.47288900 5.75843800 -0.26504900
C 1.25817400 5.52530700 -0.35923900
C 6.56651700 -3.16191800 -0.71505300
C 7.94916700 -2.76323900 -0.71782100
H 8.70976500 -3.53694200 -0.71794400
C 8.30788100 -1.45534400 -0.71437400
H 9.35728200 -1.17873700 -0.71009300
C 7.32198900 -0.40538000 -0.71254400
C 4.80199800 -4.85975200 -0.69108300
C 6.15412100 -4.49824900 -0.70493000
H 6.90486200 -5.28263700 -0.70200600
C 3.06419000 -6.57130700 -0.63971400
C 4.38550800 -6.22209900 -0.67683600
H 5.14305000 -6.99936600 -0.69842100
C 0.31584100 -7.31593300 -0.56888500
C 1.36468100 -8.31080200 -0.62344100
H 1.07229700 -9.35552200 -0.63703600
C 2.65986900 -7.96004300 -0.64758400
H 3.43860500 -8.71483200 -0.68094400
C -1.00071700 -7.68129000 -0.53153300
H -1.26397200 -8.73440200 -0.55736400
C -2.04620700 -6.71645900 -0.45667500
C -3.39654000 -7.08540400 -0.44291300
H -3.65037900 -8.13986400 -0.49304900
C -4.42507300 -6.14048100 -0.37478400
C -6.46534800 -4.14190100 -0.25970400
C -6.78869400 -5.54402200 -0.31898100
H -7.83467500 -5.83339600 -0.32417400
C -5.82048600 -6.49220000 -0.37134400
H -6.08838700 -7.54257000 -0.41703700
C -7.43580900 -3.13846000 -0.23061400
H -8.48408600 -3.42530100 -0.24374300
C -7.11206700 -1.76557700 -0.20159100
C -8.09085800 -0.73904900 -0.20340900
H -9.14122700 -1.01579900 -0.21063200
C -7.74400100 0.60004400 -0.20093900
C -8.72656900 1.65923800 -0.20155500
H -9.77676100 1.38705800 -0.19781000
C -6.96625200 3.35315000 -0.21314300
C -8.35622700 2.95686200 -0.20528600
H -9.10796200 3.73934700 -0.20404100
C -5.21662200 5.05187500 -0.23009200
C -6.57593200 4.67911500 -0.21595600
H -7.33176300 5.45855900 -0.20706100
C -3.49214700 6.77391200 -0.23072100
C -4.81679700 6.41093400 -0.22111500
H -5.58212500 7.18089900 -0.20122400
C -0.74158600 7.51769500 -0.25329700
C -1.78972800 8.51265900 -0.21459400
H -1.49917200 9.55794300 -0.19419200
C -3.08799600 8.16171000 -0.20688400
H -3.86542700 8.91823200 -0.18039300
C 0.58543200 7.87152200 -0.26483100
H 0.85898700 8.92201700 -0.23810900
C 1.61512600 6.89992900 -0.31782700
C 2.97697100 7.26295300 -0.33523700
H 3.23770400 8.31639700 -0.29758700
C 3.98147100 6.31576100 -0.39879400
C 6.03911000 4.33229200 -0.53555700
C 5.38205700 6.67350800 -0.41085100
H 5.63769300 7.72726600 -0.36594200
C 6.35467200 5.74082400 -0.47333100
H 7.39916700 6.03490200 -0.47877300
C 6.68388400 1.97311700 -0.65352300
C 7.01286700 3.35143200 -0.59266500
H 8.05965900 3.64127600 -0.58675200
C 7.65507000 0.95005100 -0.68438900
H 8.70529500 1.22994300 -0.67268100
O 1.00595500 -0.82777600 -1.05517500
O -0.79162200 -1.34393700 0.14890500
Hg 1.57958200 -0.64661800 3.20888800
C H O 0
6-31G**
****
Hg 0
S 3 1.00
48.013786990 0.58008164912E-02
21.239875095 -0.17328165235
15.876100879 0.36416685034
S 1 1.00
5.5180531820 1.0000000
S 1 1.00
1.5099145479 1.0000000
S 1 1.00
0.70683482159 1.0000000
S 1 1.00
0.16266265576 1.0000000
S 1 1.00
0.56960589112E-01 1.0000000
P 4 1.00
17.500000000 -0.86457187331E-01
15.252594855 0.15327318234
6.4404715169 -0.30274955249
1.8180159920 0.54361285026
P 1 1.00
0.90067981801 1.0000000
P 1 1.00
0.41304540122 1.0000000
P 1 1.00
0.11845702655 1.0000000
P 1 1.00
0.36087619349E-01 1.0000000
D 4 1.00
10.028197701 0.39922353165E-01
7.5920661493 -0.94774230522E-01
1.9144256118 0.31164157463
0.88641552102 0.45567078406
D 1 1.00
0.38154235320 1.0000000
D 1 1.00
0.14880740086 1.0000000
F 1 1.00
1.5521700 1.0000000
F 1 1.00
0.5243100 1.0000000
G 1 1.00
1.2895000 1.0000000
****
HG 0
HG-ECP 3 60
f-ul potential
1
2 3.88579112 30.36499643
s-ul potential
3
2 12.98154870 275.73721174
2 6.49077440 49.08921249
2 3.88579112 -30.36499643
p-ul potential
3
2 10.53809580 241.54007398
2 5.26904790 27.39659081
2 3.88579112 -30.36499643
d-ul potential
3
2 8.10172051 127.86700761
2 4.05086026 16.60831151
2 3.88579112 -30.36499643
5(E16.9)
-1
! Orbital: 1 Occ: 2.000000 from fragment 1
0.100000000E-04 0.500000000E-04-0.100000000E-04-0.300000000E-04 0.000000000E+00
-0.110000000E-02 0.150000000E-03 0.301000000E-02 0.150000000E-03-0.200000000E-04
0.000000000E+00 0.100000000E-04 0.300000000E-04 0.000000000E+00 0.000000000E+00
0.000000000E+00-0.300000000E-04 0.000000000E+00-0.200000000E-04 0.000000000E+00
0.157000000E-02 0.117000000E-02-0.162000000E-02-0.110000000E-03 0.200000000E-04
0.300000000E-04 0.100000000E-04 0.000000000E+00 0.000000000E+00-0.100000000E-04
0.000000000E+00 0.100000000E-03 0.500000000E-04 0.000000000E+00-0.200000000E-04
-0.101000000E-02-0.283000000E-02-0.128000000E-02 0.130000000E-03-0.400000000E-04
-0.800000000E-04 0.000000000E+00 0.100000000E-04 0.000000000E+00-0.100000000E-04
-0.500000000E-04 0.300000000E-04-0.130000000E-03-0.480000000E-03-0.170000000E-03
0.234000000E-02 0.250000000E-02 0.406000000E-02 0.850000000E-03 0.300000000E-04
0.160000000E-03-0.900000000E-04 0.900000000E-04 0.300000000E-04 0.110000000E-03
0.200000000E-04 0.230000000E-03-0.300000000E-04 0.400000000E-04 0.100000000E-04
-0.340000000E-02 0.362000000E-02-0.199000000E-02-0.340000000E-03-0.130000000E-03
Dear Dr Barroso,
I am runnunig SP energy calculations on ruthenocene with external basis set functions (LANL2DZ) and they are doing well up to some point at which the job terminates because of the following error:
IErr= 2 Err= 1.74D-10 NSaved= 2 NIter= 3000 NITot= 3001. DSDGMF failed for ItDFit. Error termination via Lnk1e in C:\G16W\l502.exe
The key words I use are the following:
#p pbepbe/genecp/svpfit gfinput gfprint pop=full pseudo=read
I tried to add SCF(XQC) in the key words, but it did not help, unfortunately, and the severe error 2070 occured again.
I would be very grateful for any help or tips.
Kind regards,
Iwona
Dear Dr Barroso
Can you please tell me how somebody could usually predict the distance between the metal and organic molecule in Gaussian 09?
Thanks