Conjuntos de funções de base Gaussianos são desenvolvidos para os átomos da primeira fila dos metais de transição Sc - Zn. Esses conjuntos de base foram construídos por meio do método da Coordenada Geradora Hartree-Fock (GCHF - Generator Coordinate Hartre-Fock) baseado em uma expansão polinomial de grau 3 para discretizar as equações Griffin-Wheeler-Hartree-Fock. Neste procedimento, a maneira na qual as equações são discretizadas está baseada em uma malha de pontos flexíveis não igualmente espaçada para cada uma das simetrias orbitais requeridas para descrever os átomos estudados, a diferença do método GCHF original é que a malha é igualmente espaçada para todas as simetrias. Inicialmente, foi gerado um conjunto de base balanceado consistindo de 23s17p13d funções gaussianas primitivas. A partir deste, um conjunto padrão de qualidade 7Z na valência foi construído e, posteriormente, enriquecido com grupos de funções de polarização e gerando assim:pGCHF - 7Z - 2f 1g, pGCHF - 7Z - 3f 2g e pGCHF - 7Z - 3f 2g1h. Energias atômicas Hartree-Fock para os dois estados eletrônicos de menor energia, para os átomos do Sc - Zn, foram calculadas com nossos conjuntos e comparadas com valores da energia numérica. Os resultados apresentaram um erro máximo de 1.02 mHartree, demostrando a capacidade desta expansão polinomial no desenvolvimento de conjuntos de base acurados para átomos dos metais de transição (MT) 3d. Cálculos em nível da Teoria do Funcional da Densidade (DFT - Density Functional Theory), utilizando nove diferentes funcionais, foram realizados com nossos conjuntos de base. A energia eletrônica total e propriedades incluíndo: geometrias otimizadas, cumprimentos de ligação e, frequências vibracionais, foram examinadas para um conjunto de sistemas moleculares (hidretos, dicloretos, dímeros, trímeros e óxidos de metais de transição). Os resultados são comparados com valores teóricos obtidos com conjuntos de base cc-pVnZ (n= Q ou 5) e com valores experimentais, quando disponiveís na literatura. Os resultados mostram que os valores de energia total cc-pV5Z, podem ser atingidos com nossos conjuntos de base com um menor número de funções de polarização. Outros cálculos moleculares dão resultados aproximados com valores experimentais e valores de referência DFT/cc-pV5Z. O ponto mais importante para ser mencionado é que os conjuntos de base da coordenada geradora requerem somente uma pequena fração de tempo computacional para alcançar a convergência, quando comparados com cálculos DFT/cc-pVQZ e DFT/cc-pV5Z. / Gaussian basis set functions have been constructed for the atoms of the first row of the transition metals from Sc to Zn. These basis sets were built by means of the Generator Coordinate Hartree-Fock (GCHF) method based on a polynomial expansion of degree 3 for the discretization of the Grifin-Wheler-Hartree-Fock equations. In this procedure, the equations were discretized through a mesh which is not equally spaced for each one of the orbital symmetries required to describe the atoms studied in this work, differently from the original method, in which the mesh is equally spaced for all symmetries. At first, it was generated a minimal basis set consisting of 23s17p13d primitive functions. Starting with this pattern set was construted a set of 7Z quality in the valence and then, this one was enriched with polarization functions classified as: 2f1g, 3f2g y 3f2g1h, which originated 3 sets of basis functions namely: pGCHF-7Z-2f1g, pGCHF-7Z-3f2g y pGCHF-7Z-3f2g1h. Hartree-Fock energies for the atoms Sc-Zn were calculated with our basis sets, for the two electronic states with the lowest energy and compared with values of numerical energy. The maximum error presented for the results was 1.02mH, showing the ability of this polynomial expansion to create accurate basis sets for atoms of third-rows transition metals. Density Functional Theory (DFT) calculations using a set of nine functionals were realized with our basis sets. The total electronic energy and properties such as: optimized geometries, bond distances and vibrational frequencies were calculated for a set of molecules (hydrides, dichlorides, dimers, trimmers and oxides). The obtained values were compared with theoretical values obtained with calculations using the basis set cc-pVnZ and with some experimental values found in the literature. The outcomes showed that reference values of the total electronic can be reproduced by ours basis sets, although ours have a less degree of polarization. Other molecular calculations yield results very close to experimental values and reference theorical values calculated with DFT/cc-pV5Z. The most important point to be mentioned here is that our generator coordinate basis sets require only a tiny fraction of the computational time when compared to DFT/cc-pV5Z calculations.

No / The design of new functional crystalline materials requires an understanding of the factors that control salt and co-crystal formation. These states often only differ in the location of the proton and are influenced by chemical and crystallographic factors. The interaction between a carboxylic acid and a pyridine is a frequently used supramolecular synthon in crystal engineering which can exist as either a co-crystal (CO2H center dot center dot center dot N) or salt (CO2-center dot center dot center dot HN+). The results of a Cambridge Structure Database search indicate that the nature of the functional groups on the pyridine play a stronger role in selection of the phase than those of the acid. However, the nature of the local hydrogen bonding of the interaction also adjusts the potential for proton transfer. This was demonstrated by ab initio modelling of the energy landscape for binary and ternary co-crystals by inclusion of varying components of the local environment.

Ab initio molecular orbital theory has been successful in predicting the stabilities of many weak complexes; typical of these are the complexes formed between water and various small molecules. To account for the correlation effect, Moller-Plesset perturbation theory truncated at the second order level was employed. In order to account for the hydrogen bonding, the 6-3lG** basis set was used. The geometry optimisations of the complexes were carried out using the Gaussian-92 suite of programs installed on a Hewlett-Packard 720 computer operating under UNIX. The interaction energies of the complexes were subjected to further analysis by applying the Morokuma decomposition scheme. The electrostatic interaction component accounts for over 40% of the total stabilisation energy in all the typical hydrogen bonded complexes. Gas phase enthalpies were computed and compared with the experimental values of similar systems. For the systems studied here, the prediction is that all complexes are stable at 25° C. A second program, Vibra, was used for carrying out a normal coordinate analysis. A third computer program for the graphical representation of molecular and crystallographic models, Schakal-92, was employed to illustrate the predicted equilibrium geometries and the fundamental vibrational modes. The predicted geometries, interaction energies, charge redistributions, vibrational wave numbers, infrared intensities and force constants are listed and compared with those in the literature, where applicable. Correlations between the various predicted properties show some interesting chemistry. / Thesis (Ph.D.)-University of Natal, Durban, 1996.

O método da Coordenada Geradora é uma poderosa ferramenta para gerar funções de base. Sua última versão, chamada de método da Coordenada Geradora polinomial, permite a geração de funções de base mais eficientes e precisas a um baixo custo computacional. Nesta tese, além da geração de funções de base para os átomos do primeiro período da Tabela Periódica, uma estratégia de selecionar expoentes da própria função de base para posteriormente refiná-los, com o intuito de gerar funções de polarização e difusas, é apresentada e testada em cálculos de propriedades eletrônicas de átomos e moléculas. / Ab initio electronic structure calculations for atoms and especially for molecules are mostly carried out within the finite basis set expansion method in the Hartree-Fock theory by Roothaan. The search for ever more efficient basis sets has been a constant quest and here we show a new alternative to develop efficient Gaussian-Type Functions (GTF) basis sets for atomic and molecular calculations by employing the Polynomial Generator Coordinate Hartree-Fock method.

Conjuntos de funções de base Gaussianos são desenvolvidos para os átomos da primeira fila dos metais de transição Sc - Zn. Esses conjuntos de base foram construídos por meio do método da Coordenada Geradora Hartree-Fock (GCHF - Generator Coordinate Hartre-Fock) baseado em uma expansão polinomial de grau 3 para discretizar as equações Griffin-Wheeler-Hartree-Fock. Neste procedimento, a maneira na qual as equações são discretizadas está baseada em uma malha de pontos flexíveis não igualmente espaçada para cada uma das simetrias orbitais requeridas para descrever os átomos estudados, a diferença do método GCHF original é que a malha é igualmente espaçada para todas as simetrias. Inicialmente, foi gerado um conjunto de base balanceado consistindo de 23s17p13d funções gaussianas primitivas. A partir deste, um conjunto padrão de qualidade 7Z na valência foi construído e, posteriormente, enriquecido com grupos de funções de polarização e gerando assim:pGCHF - 7Z - 2f 1g, pGCHF - 7Z - 3f 2g e pGCHF - 7Z - 3f 2g1h. Energias atômicas Hartree-Fock para os dois estados eletrônicos de menor energia, para os átomos do Sc - Zn, foram calculadas com nossos conjuntos e comparadas com valores da energia numérica. Os resultados apresentaram um erro máximo de 1.02 mHartree, demostrando a capacidade desta expansão polinomial no desenvolvimento de conjuntos de base acurados para átomos dos metais de transição (MT) 3d. Cálculos em nível da Teoria do Funcional da Densidade (DFT - Density Functional Theory), utilizando nove diferentes funcionais, foram realizados com nossos conjuntos de base. A energia eletrônica total e propriedades incluíndo: geometrias otimizadas, cumprimentos de ligação e, frequências vibracionais, foram examinadas para um conjunto de sistemas moleculares (hidretos, dicloretos, dímeros, trímeros e óxidos de metais de transição). Os resultados são comparados com valores teóricos obtidos com conjuntos de base cc-pVnZ (n= Q ou 5) e com valores experimentais, quando disponiveís na literatura. Os resultados mostram que os valores de energia total cc-pV5Z, podem ser atingidos com nossos conjuntos de base com um menor número de funções de polarização. Outros cálculos moleculares dão resultados aproximados com valores experimentais e valores de referência DFT/cc-pV5Z. O ponto mais importante para ser mencionado é que os conjuntos de base da coordenada geradora requerem somente uma pequena fração de tempo computacional para alcançar a convergência, quando comparados com cálculos DFT/cc-pVQZ e DFT/cc-pV5Z. / Gaussian basis set functions have been constructed for the atoms of the first row of the transition metals from Sc to Zn. These basis sets were built by means of the Generator Coordinate Hartree-Fock (GCHF) method based on a polynomial expansion of degree 3 for the discretization of the Grifin-Wheler-Hartree-Fock equations. In this procedure, the equations were discretized through a mesh which is not equally spaced for each one of the orbital symmetries required to describe the atoms studied in this work, differently from the original method, in which the mesh is equally spaced for all symmetries. At first, it was generated a minimal basis set consisting of 23s17p13d primitive functions. Starting with this pattern set was construted a set of 7Z quality in the valence and then, this one was enriched with polarization functions classified as: 2f1g, 3f2g y 3f2g1h, which originated 3 sets of basis functions namely: pGCHF-7Z-2f1g, pGCHF-7Z-3f2g y pGCHF-7Z-3f2g1h. Hartree-Fock energies for the atoms Sc-Zn were calculated with our basis sets, for the two electronic states with the lowest energy and compared with values of numerical energy. The maximum error presented for the results was 1.02mH, showing the ability of this polynomial expansion to create accurate basis sets for atoms of third-rows transition metals. Density Functional Theory (DFT) calculations using a set of nine functionals were realized with our basis sets. The total electronic energy and properties such as: optimized geometries, bond distances and vibrational frequencies were calculated for a set of molecules (hydrides, dichlorides, dimers, trimmers and oxides). The obtained values were compared with theoretical values obtained with calculations using the basis set cc-pVnZ and with some experimental values found in the literature. The outcomes showed that reference values of the total electronic can be reproduced by ours basis sets, although ours have a less degree of polarization. Other molecular calculations yield results very close to experimental values and reference theorical values calculated with DFT/cc-pV5Z. The most important point to be mentioned here is that our generator coordinate basis sets require only a tiny fraction of the computational time when compared to DFT/cc-pV5Z calculations.

Computational chemistry lies at the intersection of chemistry, physics, mathematics, and computer science, and can be used to explain the behavior of atoms and molecules, as well as to augment experiment. In this work, computational chemistry methods are used to predict structural and energetic properties of small molecules, i.e. molecules with less than 60 atoms. Different aspects of computational chemistry are examined in this work. The importance of examining the converged orbitals obtained in an electronic structure calculation is explained. The ability to more completely describe the orbital space through the extrapolation of energies obtained at increasing quality of basis set is investigated with the use of the Sapporo-nZP-2012 family of basis set. The correlation consistent Composite Approach (ccCA) is utilized to compute the enthalpies of formation of a set of molecules and the accuracy is compared with the target method, CCSD(T,FC1)/aug-cc-pCV∞Z-DK. Both methodologies are able to produce computed enthalpies of formation that are typically within 1 kcal mol-1 of reliable experiment. This demonstrates that ccCA can be used instead of much more computationally intensive methods (in terms of memory, processors, and time required for a calculation) with the expectation of similar accuracy yet at a reduced computational cost. The enthalpies of formation for systems containing s-block elements have been computed using the multireference variant of ccCA (MR-ccCA), which is designed specifically for systems that require an explicit treatment of nondynamical correlation. Density functional theory (DFT) has been used for the prediction of the structural properties of a set of lanthanide trihalide molecules as well as the reaction energetics for the rearrangement of diphosphine ligands around a triosmium cluster.

O Método da Coordenada Geradora Hartree-Fock Polinomial (pMCG-HF), desenvolvido por R.C. Barbosa e A.B.F. da Silva [1], é uma ferramenta matemática valiosa que permite gerar funções de base (também conhecidas como conjuntos de base). As funções de base geradas por este método têm um bom comportamento e são capazes de calcular valores precisos de propriedades eletrônicas moleculares. Porém, depois de gerar funções de base do hidrogênio até o flúor [2], fez-se necessário a adição de expoentes à função de base, correspondentes a cada átomo, para melhor adaptação à realização dos cálculos moleculares. Estas funções adicionais são o que chamamos de funções de polarização. A adição de funções de polarização, através de otimização computacional, é muito custosa, deste modo o desenvolvimento de uma regra de polarização para se esquivar desta otimização é de grande importância e por isso se transforma na beleza e no objetivo deste trabalho. Portanto, nesta dissertação, estudar-se-á um procedimento para escolher funções de polarização que reduza drasticamente o tempo computacional, no sentido de permitir uma seleção, mais simples, de expoentes da própria função de base primitiva para serem usadas nas funções de polarização p, d, f, g, etc. para a obtenção de propriedades moleculares calculadas através de métodos químico-quânticos / The polynomial generate coordinate method pGCM developed by R.C. Barbosa and A.B.F. da Silva [1] is an remarkble mathematic tool for the generation of basis functions (also known as basis sets). The basis sets generated from this method have a good behavior and are able to produce accurate values for electronic molecular properties. In fact, after generating a basis set [2] we need to add a set of exponent functions in order to better adequate a basis set to perform molecular calculations. These sets of additional functions are called polarizations functions. This work provides a methodology where the polarization functions are obtained from the initial basis set (the primitive set) without optimizing them separately by using optimization algorithms that are, computationally speaking, very costly. This procedure reduces drastically the computational time used to find polarization functions to be used in molecular quantum chemical calculations. Our methodology permits to choose the polarization functions directly from the primitive orbital exponents of each atomic symmetry s, p, d, f etc. in a very simple manner. The finding of polarization functions using our methodology was performed with several quantum chemical methods.

A afinidade eletrônica (AE) é uma importante propriedade de átomos e moléculas, sendo definida como a diferença de energia entre a espécie neutra e seu respectivo íon negativo. Uma vez que a AE é uma fração muito pequena da energia eletrônica total das espécies neutra e aniônica, é necessário que tais energias sejam determinadas com elevado grau de precisão. A receita utilizada para o cálculo teórico acurado da AE atômica e molecular baseia-se na escolha de um conjunto adequado de funções de base juntamente com o emprego de teorias com altos níveis de correlação eletrônica. Durante o cálculo, o mesmo conjunto de base é utilizado para descrever o elemento neutro e seu respectivo ânion. Geralmente, os conjuntos de base para descrever propriedades de ânions possuem seus expoentes otimizados em ambiente neutro, e sua difusibilidade é conferida pela adição de funções difusas para cada valor de momento angular, l. A ideia deste trabalho está no desenvolvimento de conjuntos de base otimizados exclusivamente em ambiente aniônico para cálculos precisos de afinidade eletrônica. Deste modo, foram escolhidos os átomos para serem estudados: B, C, O e F. Os conjuntos de base foram gerados pelo Método da Coordenada Geradora Hartree-Fock, empregando a técnica da Discretização Integral Polinomial para a solução das integrais do problema. Os conjuntos de base obtidos são compostos por (18s13p) primitivas que foram contraídos para [7s6p] via esquema de contração geral proposto por Raffenetti. Os conjuntos contraídos foram polarizados para 4d3f2g e 4d3f2g1h, sendo os expoentes otimizados em ambiente CISD através do método SIMPLEX. Avaliaram-se as funções de base no cálculo de afinidades eletrônicas, tendo seus resultados comparados aos obtidos utilizando as bases aug-cc-pVQZ e aug-cc-pV5Z. A análise dos resultados demonstrou que os conjuntos de base difusos, gerados neste trabalho, reproduzem de maneira satisfatória as afinidades eletrônicas em relação ao valor experimental. Os conjuntos difusos polarizados para 4d3f2g1h apresentaram eficiência superior aos conjuntos aug-cc-pVQZ e, em alguns casos, aos conjuntos aug-cc-pV5Z que são consideravelmente maiores. / The electron affinity (EA) is an important property of atoms and molecules defined as the energy difference between the neutral species and its negative ion. Since the EA is a very small fraction of the total electronic energy of anionic and neutral species, one must determine these energies with high accuracy. The recipe used to calculate accurate atomic and molecular EAs is based on the choice of an adequate basis set and the use of high level of electron correlation calculations. In the computation of EAs, the same basis set is used to describe both neutral and negatively charged species. In general, the basis sets designed to describe anionic properties have their exponents optimized in neutral environment, and its diffuseness is acquired through the addition of diffuse functions for each angular momentum. The main idea of this work is to develop basis sets optimized exclusively in anionic environment that would be applied in accurate calculations of electron affinity. Thus, here follows the chosen atoms to be studied: B, C, O and F. The basis sets were generated by the Generator Coordinate Hartree-Fock Method through the Polynomial Integral Discretization Method. Basis sets were obtained containing (18s13p) primitives that were contracted to [7s6p] via Raffenetti\'s general contraction scheme. The contracted basis sets were polarized to 4d3f2g and 4d3f2g1h, and the exponents of polarization were optimized in a CISD environment through the Simplex algorithm. The basis sets quality was evaluated through the calculation of the electron affinities. The results were compared to those obtained by using the aug-cc-pVQZ and aug-cc-pV5Z basis-sets. The calculation showed that our diffuse basis sets reproduce satisfactorily the electron affinities when compared to the experimental data. The diffuse basis sets polarized to 4d3f2g1h showed to be more efficient than the aug-cc-pVQZ basis sets and in some cases also better than the aug-cc-pV5Z basis sets that are considerably larger.

Orientador: Pedro Antonio Muniz Vazquez / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Quimica / Made available in DSpace on 2018-08-14T17:56:25Z (GMT). No. of bitstreams: 1 Vidal_LucianoNassif_D.pdf: 1273955 bytes, checksum: 79e28497883fdc3e8ebe5907183a15d9 (MD5) Previous issue date: 2009 / Resumo: O presente trabalho aborda o desenvolvimento de metodologias para o cálculo das intensidades absolutas do espalhamento Raman vibracional produzido por moléculas em fase gasosa. Com o objetivo de reduzir a demanda por recursos computacionais nestes cálculos, foram desenvolvidas duas novas famílias de funções de base compactas pela aplicação do método de polarização elétrica de Sadlej às bases para uso com potenciais efetivos de caroço SBKJC e de Stuttgart-Colônia. Utilizando estas novas funções de base, as intensidades Raman podem ser obtidas com a mesma qualidade das bases Sadlej-pVTZ, que são referência no cálculo destas propriedades, porém com um custo computacional sensivelmente menor. Além disso, como estes pseudopotenciais foram modelados para descrever os efeitos relativísticos sobre os elétrons internos, as polarizabilidades e intensidades Raman obtidas no nível Hartree-Fock com estas novas bases concordam, dentro de um erro médio de 6%, com seus respectivos valores relativísticos Dirac-Hartree-Fock/Sadlei-pVTZ com hamiltoniano de Dirac-Coulomb. Também foi desenvolvida uma metodologia para o estudo das intensidades das transições Raman fundamentais, de combinação e sobretom, que inclui as correções para a anarmonicidade cúbica do potencial, introduzidas através de uma transformação de contato. Os resultados obtidos para a molécula de acetileno e seus isotopômeros deuterados mostram que a anarmonicidade mecânica exerce grande influência sobre as intensidades Raman, particularmente das transições de segunda ordem. Excetuando as transições de combinação, em geral, as correções de anarmoniciadade melhoram a concordância dos valores teóricos com os experimentais. Uma terceira parte deste trabalho trata do efeito Raman em condições ressonantes, onde uma expressão para estas intensidades foi derivada, implementada no programa PLACZEK e aplicada no cálculo do espectro Raman da molécula de trans-butadieno nas vizinhanças de sua transição eletrônica 1Bu. Este estudo mostrou que as aproximações utilizadas com maior frequência para simplificar o cálculo desta propriedade afetam significativamente as seções de choque desta molécula, sugerindo que estas aproximações devem ser evitadas em estudos desta natureza. / Abstract: In this work new methodologies for the calculation of absolute vibrational Raman intensities of gaseous systems are presented. In order to reduce the computational requirements in these calculations two families of compact basis functions were generated from the effective core potential valence basis sets SBKJC and Stuttgart-Cologne through the Sadlej's electric polarization procedure. The Raman intensities evaluated with the new bases are close to those obtained with the well successful Sadlej-pVTZ basis but the computational requirements are significatively reduced. Furthermore, since the effective core potentials SBKJC and Stuttgart-Cologne were developed to account for the relativistic effects on the inner electrons, the polarizabilities and Raman intensities evaluated at the Hartree-Fock level with the new bases agree with the relativistic Dirac-Hartree-Fock values, obtained using the Dirac-Coulomb Hamiltonian and the Sadlej-pVTZ set, within the mean error of 6%. In the second part of this work a methodology was developed for the study of fundamental, combination and overtone Raman transitions including a treatment based on the contact transform formalism for the mechanical anharmonicity from the cubic potential energy terms. The results obtained for acetylene and its deutered isotopomers show that anharmonicity effects on the Raman intensities can be very strong, particularly in the second order transitions. With the exception of the combination transitions, in general the corrections for mechanical anharmonicity improve the agreement between ab initio and experimental values. The resonance Raman scattering is the subject of the third part of this work where an expression for the resonance cross section was derived, implemented in the PLACZEK program and applied to the calculation of the resonance Raman spectrum of the trans-butadiene molecule in the region of its. / Doutorado / Físico-Química / Doutor em Ciências

O Método da Coordenada Geradora Hartree-Fock Polinomial (pMCG-HF), desenvolvido por R.C. Barbosa e A.B.F. da Silva [1], é uma ferramenta matemática valiosa que permite gerar funções de base (também conhecidas como conjuntos de base). As funções de base geradas por este método têm um bom comportamento e são capazes de calcular valores precisos de propriedades eletrônicas moleculares. Porém, depois de gerar funções de base do hidrogênio até o flúor [2], fez-se necessário a adição de expoentes à função de base, correspondentes a cada átomo, para melhor adaptação à realização dos cálculos moleculares. Estas funções adicionais são o que chamamos de funções de polarização. A adição de funções de polarização, através de otimização computacional, é muito custosa, deste modo o desenvolvimento de uma regra de polarização para se esquivar desta otimização é de grande importância e por isso se transforma na beleza e no objetivo deste trabalho. Portanto, nesta dissertação, estudar-se-á um procedimento para escolher funções de polarização que reduza drasticamente o tempo computacional, no sentido de permitir uma seleção, mais simples, de expoentes da própria função de base primitiva para serem usadas nas funções de polarização p, d, f, g, etc. para a obtenção de propriedades moleculares calculadas através de métodos químico-quânticos / The polynomial generate coordinate method pGCM developed by R.C. Barbosa and A.B.F. da Silva [1] is an remarkble mathematic tool for the generation of basis functions (also known as basis sets). The basis sets generated from this method have a good behavior and are able to produce accurate values for electronic molecular properties. In fact, after generating a basis set [2] we need to add a set of exponent functions in order to better adequate a basis set to perform molecular calculations. These sets of additional functions are called polarizations functions. This work provides a methodology where the polarization functions are obtained from the initial basis set (the primitive set) without optimizing them separately by using optimization algorithms that are, computationally speaking, very costly. This procedure reduces drastically the computational time used to find polarization functions to be used in molecular quantum chemical calculations. Our methodology permits to choose the polarization functions directly from the primitive orbital exponents of each atomic symmetry s, p, d, f etc. in a very simple manner. The finding of polarization functions using our methodology was performed with several quantum chemical methods.