Neste trabalho, a reação O(3P) + HBr → OH + Br, que se insere no contexto da química atmosférica, mais especificamente nos ciclos de destruição catalítica do ozônio, foi estudada empregando superfícies de energia potencial ab initio de alta qualidade. Construímos superfícies para os estados eletrônicos 3A\" e 3A\', baseadas na expansão de muitos corpos, com potenciais de dois corpos ajustados por funções \"switched-MLJ\" e potenciais de três corpos interpolados utilizando o método RKHS. Os pontos ab initio utilizados foram obtidos no nível de teoria MRCISD+Q/CBS+SO. A superfície final para o estado eletrônico 3A\" foi corrigida por um fator multiplicativo de modo que reproduzisse cálculos benchmark para a altura de barreira, em que foram considerados efeitos relativísticos (escalares e spin-órbita), correlação caroço-valência e correlação de ordem alta (excitações triplas e quádruplas). As principais características da SEP 3A\" são dois mínimos de van der Walls no canal de entrada e saída e um estado de transição com uma barreira energética de 5,01 kcal mol-1. A SEP 3A\' tem um estado de transição linear em 6,45 kcal mol-1. O valor obtido para a entalpia de reação, a 0 K, foi de -15,7 kcal mol-1, muito próximo do valor experimental de -15,14 kcal mol-1. Calculamos constantes de velocidade, no intervalo de 200 a 1000 K, utilizando a teoria do estado de transição variacional com contribuições de tunelamento multidimensional (ICVT/µOMT) e com uma aproximação de espalhamento reativo quântico (QM/JS). Estes valores estão em ótima concordância com os dados experimentais da literatura, em todo intervalo de temperatura em que estão disponíveis: de 221 a 554 K para reação O + HBr e de 295 a 419 K para reação O + DBr. As constantes de velocidade, a 298 K e em cm3 molécula-1 s-1, obtidas para a reação O + HBr são: 3,62·10-14 (ICVT/µOMT) e 3,35·10-14 (QM/JS), enquanto que o valor experimental é 3,78·10-14. A qualidade destes resultados reforça nossa confiança nos procedimentos e aproximações utilizadas, abrindo caminho para a caracterização, em alto nível, de uma grande variedade de reações em fase gasosa. / In this work, the O(3P) + HBr → OH + Br reaction, which is relevant to atmospheric chemistry, specially for the catalytic ozone depletion, was studied using high-level ab initio potential energy surfaces. We constructed surfaces for the 3A\" and 3A\' electronic states, based on the many-body expansion, with the two-body potentials adjusted by the switched-MLJ function and the three-body potentials interpolated using the RKHS method. The ab initio points were calculated at the MRCISD+Q/CBS+SO level of theory. The main features of the 3A\" are the presence of two van der Waals minima, at the entrance and exit channels, and a transition state with a barrier height of 5.01 kcal mol-1. The 3A\' PES has a linear transition state at 6.45 kcal mol-1. We obtained the enthalpy of reaction, at 0 K, of -15.7 kcal mol-1, in close agreement with the experimental value of -15.14 kcal mol-1. Rate constants, in the temperature range from 200 to 1000 K, were calculated using the variational transition state theory with contributions of multidimensional tunneling (ICVT/µOMT) and also a quantum reactive scattering approach (QM/JS). Their values are in fair agreement with the experimental data in the literature in the whole temperature range available: from 221 to 554 K for the O + HBr reaction, and from 295 to 419 K for the O + DBr reaction. At 298 K, the calculated rate constants, in cm3 molecule-1 s-1, for the O + HBr reaction are 3.62·10-14 (ICVT/µOMT) and 3.35·10-14 (QM/JS); and the experimental value is 3.78·10-14. The quality of these results reinforces our confidence in the procedures and approximations used, leading to the possibility of high-level characterization of a variety of gas phase reactions.

Este trabalho teve como foco principal a caracterização espectroscópica dos haletos de berílio, BeCl, BeBr e BeI. O conhecimento acerca dessas espécies foi ampliado significativamente através do cálculo de parâmetros de difícil caracterização experimental. Como os poucos trabalhos experimentais realizados para esses haletos apresentam resultados inconclusivos para a transição C 2Σ+ - X 2Σ+, através do cálculo acurado de probabilidades de transição foi possível oferecer uma nova interpretação para os dados existentes, apresentando inclusive resultados para as outras transições eletrônicas até o presente desconhecidas. Cálculos utilizando o método de Interação de Configurações Multirreferencial (MRCI) com um conjunto de funções de base consistentes na correlação de qualidade quintupla-zeta foram utilizados para se obter as curvas de energia potencial associadas ao primeiro e segundo canais de dissociação para os estados dubleto e quarteto dessas espécies. Parâmetros espectroscópicos como as constantes vibracionais we, wexe, a constante rotacional Be e a distância de equilíbrio foram determinados para os estados de mais baixa energia. Além disso, foram calculadas a energia de excitação (Te) e a energia de dissociação (De) com e sem a inclusão de efeitos spin-órbita. Os resultados obtidos expandiram significativamente nosso conhecimento sobre os estados eletrônicos dessa espécies. Para os estados já investigados experimentalmente, houve boa coerência entre os parâmetros calculados e experimentais. As diferenças de energia entre os dois canais no limite de dissociação para as moléculas BeCl, BeBr e BeI obtidas neste trabalho foram, respectivamente: 21835, 21889 e 21998 cm-1, em concordância com o respectivo resultado experimental das três espécies (21980 cm-1). A determinação teórica da energia de dissociação D0 foi bastante satisfatória. Obtivemos 92,24; 72,77 e 51,75 kcal·mol-1, respectivamente, para as moléculas BeCl, BeBr e BeI incluindo os efeitos spin-órbita, comparados a 91,78, 71 e 57 kcal·mol-1. Para uma análise mais completa das curvas de energia potencial, elas também foram obtidas considerando-se os efeitos relativísticos. A constante de acoplamento spin-órbita calculada na região de equilíbrio de cada uma das moléculas BeCl, BeBr e BeI foi, respectivamente: 41, 207 e 324 cm-1, em boa concordância com os resultados experimentais: 52,8, 198 e 361,1 cm-1, respectivamente. A separação spin-órbita no limite de dissociação calculada foi de 823, 3446 e 6975 cm-1 (BeCl, BeBr e BeI), também coerentes com resultados experimentais: 882, 3685,24 e 7603,15 cm-1, respectivamente. Foi ainda realizada uma análise detalhada sobre os canais de dissociação das moléculas HBeP e BePH. A construção de um diagrama de energias relativas permitiu obter um conjunto de dados importantes a respeito de todos os possíveis canais de dissociação. Em particular é mostrado a estabilidade relativa e a energia no limite de dissociação. Para a energia no limite de dissociação obtivemos resultados consistentes com os valores experimentais. Nossos cálculos para os três primeiros canais de dissociação HBeP = H + Be resultaram nos seguintes valores: H (2Sg) + Be (X 4Σ+), 0 cm-1; H (2Sg) + Be (a 2π), 426 cm-1 (exp. 502 cm-1); H (2Sg) + Be (b 2Σ-), 1896 cm-1 (exp. 1976 cm-1). Por fim, pela primeira vez na literatura, o calor de formação e de atomização foram calculados para estas espécies. Para os isômeros HBeP e BePH, a energia de atomização a 298,15 K calculada foi de 119,02 kcal·mol-1 e 107,40 kcal·mol-1, respectivamente; para o calor de formação a 298,15 K, obtivemos 86,14 e 97,76 kcal·mol-1 para as espécies HBeP e BePH, respectivamente. / This work has mainly focused on the spectroscopic characterization of the beryllium halides, BeCl, BeBr and BeI. Knowledge about these species was substantially increased through the calculation of parameters that are difficult to access in an experimental characterization. As the few experimental works carried out for these halides presented inconclusive results for the C 2Σ+ - X2Σ+ transition, we carried out accurate calculations of transition probabilities, thus offering a new interpretation about the existing data, including also results for other electronic transitions as yet unknown. Calculations using the method Multireference Configuration Interaction (MRCI) along with correlation-consistent basis set functions of quintuple-zeta quality were used to obtain the potential energy curves associated with the first and second dissociation channels for the doublet and quartet states of these species. Spectroscopic parameters such as vibrational constants We, WeXe, the rotational constant Be and the equilibrium distance were determined for the low-lying states. Furthermore, we also calculated the excitation energy (Te) and the dissociation energy (De) with and without the inclusion of spin-orbit effects. The obtained results significantly expanded our knowledge about the electronic states of this species. For states already investigated experimentally, there was good consistency between the calculated and the experimental parameters. The energy differences between the two channels in the dissociation limit for BeCI, BeBr and Bel molecules obtained in this work were: 21835, 21889 and 21998 cm-1, in agreement with the corresponding experimental results of the three species (21 980 cm-1). The theoretical determination of the dissociation energy D0 was very satisfactory. We obtained 92.24, 72.77, and 51.75 kcal·mol-1, respectively, for the BeCl, BeBr and BeI molecules including spin-orbit effects, compared to 91.78, 71 and 57 kcal . mol-1. For a more complete analysis of the potential energy curves, they also were obtained considering relativistic effects. The calculated spin-orbit coupling constants in the equilibrium region of BeCl, BeBr and BeI molecules were respectively 41, 207 and 324 cm-1, in good agreement with the experimental results: 52.8, 198, and 361.1 cm-1, respectively. The calculated spin-orbit splitting in the dissociation limit was 823, 3446 and 6975 cm-1 (BeCl, BeBr and BeI) also consistent with the experimental results: 882, 3685.24 and 7603.15 cm-1, respectively. Additionally, a detailed analysis of the dissociation channels of HBeP and BePH molecules was performed. The construction of a relative energies diagram allowed us to derive a set of important data for all the possible dissociation channels. In particular, it is shown the relative stability and the energy in the dissociation limit. For the energy in the dissociation limit, we obtained results consistent with the experimental values. Our calculations for the first three dissociation channels HBeP = H + Be resulted in the following values: H (2Sg) + Be (X 4Σ+), 0 cm-1; H (2Sg) + Be (a 2π), 426 cm-1 (exp. 502 cm-1); H (2Sg) + Be (b 2Σ-), 1896 cm-1 (exp. 1976 cm-1). Finally, for the first time in the literature, the heat of formation and the atomization energy were calculated for these species. For the HBeP and BePH isomers, the atomization energy calculated at 298.15 K was 119.02 and 107.40 kcal . mol-1, respectively; for the heat of formation at 298.15 K, we obtained 86.14 and 97.76 kcal· mol-1 for HBeP and BePH species, respectively.

Neste trabalho, a reação O(3P) + HBr → OH + Br, que se insere no contexto da química atmosférica, mais especificamente nos ciclos de destruição catalítica do ozônio, foi estudada empregando superfícies de energia potencial ab initio de alta qualidade. Construímos superfícies para os estados eletrônicos 3A\" e 3A\', baseadas na expansão de muitos corpos, com potenciais de dois corpos ajustados por funções \"switched-MLJ\" e potenciais de três corpos interpolados utilizando o método RKHS. Os pontos ab initio utilizados foram obtidos no nível de teoria MRCISD+Q/CBS+SO. A superfície final para o estado eletrônico 3A\" foi corrigida por um fator multiplicativo de modo que reproduzisse cálculos benchmark para a altura de barreira, em que foram considerados efeitos relativísticos (escalares e spin-órbita), correlação caroço-valência e correlação de ordem alta (excitações triplas e quádruplas). As principais características da SEP 3A\" são dois mínimos de van der Walls no canal de entrada e saída e um estado de transição com uma barreira energética de 5,01 kcal mol-1. A SEP 3A\' tem um estado de transição linear em 6,45 kcal mol-1. O valor obtido para a entalpia de reação, a 0 K, foi de -15,7 kcal mol-1, muito próximo do valor experimental de -15,14 kcal mol-1. Calculamos constantes de velocidade, no intervalo de 200 a 1000 K, utilizando a teoria do estado de transição variacional com contribuições de tunelamento multidimensional (ICVT/µOMT) e com uma aproximação de espalhamento reativo quântico (QM/JS). Estes valores estão em ótima concordância com os dados experimentais da literatura, em todo intervalo de temperatura em que estão disponíveis: de 221 a 554 K para reação O + HBr e de 295 a 419 K para reação O + DBr. As constantes de velocidade, a 298 K e em cm3 molécula-1 s-1, obtidas para a reação O + HBr são: 3,62·10-14 (ICVT/µOMT) e 3,35·10-14 (QM/JS), enquanto que o valor experimental é 3,78·10-14. A qualidade destes resultados reforça nossa confiança nos procedimentos e aproximações utilizadas, abrindo caminho para a caracterização, em alto nível, de uma grande variedade de reações em fase gasosa. / In this work, the O(3P) + HBr → OH + Br reaction, which is relevant to atmospheric chemistry, specially for the catalytic ozone depletion, was studied using high-level ab initio potential energy surfaces. We constructed surfaces for the 3A\" and 3A\' electronic states, based on the many-body expansion, with the two-body potentials adjusted by the switched-MLJ function and the three-body potentials interpolated using the RKHS method. The ab initio points were calculated at the MRCISD+Q/CBS+SO level of theory. The main features of the 3A\" are the presence of two van der Waals minima, at the entrance and exit channels, and a transition state with a barrier height of 5.01 kcal mol-1. The 3A\' PES has a linear transition state at 6.45 kcal mol-1. We obtained the enthalpy of reaction, at 0 K, of -15.7 kcal mol-1, in close agreement with the experimental value of -15.14 kcal mol-1. Rate constants, in the temperature range from 200 to 1000 K, were calculated using the variational transition state theory with contributions of multidimensional tunneling (ICVT/µOMT) and also a quantum reactive scattering approach (QM/JS). Their values are in fair agreement with the experimental data in the literature in the whole temperature range available: from 221 to 554 K for the O + HBr reaction, and from 295 to 419 K for the O + DBr reaction. At 298 K, the calculated rate constants, in cm3 molecule-1 s-1, for the O + HBr reaction are 3.62·10-14 (ICVT/µOMT) and 3.35·10-14 (QM/JS); and the experimental value is 3.78·10-14. The quality of these results reinforces our confidence in the procedures and approximations used, leading to the possibility of high-level characterization of a variety of gas phase reactions.

Simulações de sais de carbonato fundidos pelo método de Dinâmica Molecular (MD) foram efetuadas com o modelo polarizável de cargas flutuantes (FC). O modelo de cargas flutuantes implementa os efeitos de polarização pelo método de Lagrangiano estendido, onde as variáveis extras são as próprias cargas parciais do íon poliatômico. O modelo FC foi parametrizado por meio de cálculos ab inito, aplicado ao ânion carbonato. Cálculos de Química Quântica ab initio foram utilizados para corroborar o modelo proposto para o ânion carbonato. Os sistemas investigados consistem em misturas de carbonatos alcalinos fundidos, Li2CO3/K2CO3, os quais são utilizados como eletrólitos em células a combustível. As simulações MD foram utilizadas para verificar o efeito da polarização dos ânions sobre a estrutura e dinâmica do líquido. Estudamos o efeito da inclusão de polarização sobre a condutividade do eletrólito. / Simulations of molten carbonate salts by the method of Molecular Dynamics (MD) have been performed with the fluctuanting charge (FC) model. The FC model implements the effect of polarization by using method of extended Lagrangian, where the extra variables are the partial charges of the poliatomic ion. The FC model was parametrized by ab inito calculations os a single carbonate anion. Quantum Chemistry calculations bave been used to corroborate the model for the carbonate anion. The investigated systems consist of alkaline carbonate mixtures, Li2CO3/K2CO3, which are used as electrolytes in fuel cells. MD simulations have been used to verify polarization effects on structure and dynamics of the liquid. We study the effect of including anion polarization on the condutivity of the electrolyte.

Neste trabalho foram sintetizados a poli-2-etinilpiridina (P2EP), poli-4etinilpiridina (P4EP), o iodeto de poli(2-N-t-butilpiridiniumilacetileno) (P2EPtBu) e a poli-β-etinilnaftaleno (P&#946:EN), os quais são poliacetilenos substituídos. Estes polímeros, juntamente com o cloreto de poli(2-piridínio-2-piridilacetileno) (P2EPH), foram caracterizados por espectroscopia vibracional no infravermelho e Raman. Estes polímeros apresentaram variação na posição das bandas Raman com a energia da radiação excitante - chamada dispersão Raman ou fotosseletividade - da ordem de 10 cm-1, bem inferior ao apresentado pelo poliacetileno (cerca de 60 cm-1). Este deslocamento foi interpretado utilizando-se dois dos modelos existentes para descrever este fenômeno: o Modelo de Modo de Amplitude (AMM) e o Modelo de Coordenada de Conjugação Efetiva (ECCM), os quais fornecem informações sobre a estrutura polimérica e sobre seus níveis eletrônicos. Utilizando-se o AMM foi possível obter informações sobre os níveis eletrônicos excitados de mesma simetria que o estado eletrônico fundamental. Por outro lado, o ECCM, com a ajuda de cálculos DFT, mostrou diferenças na extensão da conjugação e no grau de dimerização entre o P2EP na forma cis e trans e indicou que este polímero apresentava, predominantemente, a estrutura cis, fato este confirmado pelos espectros no infravermelho. A dopagem com I2 provocou efeitos diferentes na estrutura dos polímeros. Os espectros no infravermelho dos polímeros dopados indicaram que o P2EP e o P2EPH apresentaram aumento na quantidade de segmentos cis enquanto o P2EPtBu apresentou diminuição na quantidade desses segmentos. Os espectros Raman dos polímeros dopados confirmaram os dados dos espectros no infravermelho. Esta diferença foi interpretada como sendo devida à diferença no volume do substituinte, pois grupos volumosos favorecem o isômero trans-cisóide onde a distância entre os substituintes é maior. A dopagem também levou a um aumento na condutividade dos polímeros, porém os valores de condutividade obtidos foram bem inferiores que os apresentados pelo poliacetileno dopado (10-5 a 10-7 contra 102 S cm-1, tipicamente). / In this work poly-2-ethynylpyridine (P2EP), poly-4-ethynylpyridine (P2EP), poly(2-N-t-butylpyridiniumylacetylene) iodide (P2EPtBu) and poly-β-pethynylnaphthalene (PβEN) were synthesized. These polymers, which are polyacetylenes derivatives, and poly(2-pyridinium-2-pyridylacetylene) chloride (P2EPH) were analyzed by vibrational techniques. It was observed that some Raman band positions shifted when excitation energy was changed, however this shift (ca. 10 cm-1) was substantially smaller than the presented by polyacetylene (ca. 60 cm-1). This phenomenon called Raman dispersion or photoselectivity - was examinated using Amplitude Mode Model (AMM) and Effective Conjugation Coordinate Model (ECCM) and gave us important information about polymer structure and its electronic energies levels. Using AMM it was possible to infere information about electronic excited states with the same symmetry of the electronic ground state. ECCM gave us evidences about differences in conjugation extension and dimerization degree between cis and trans structures. Besides this, ECCM was successful in indicating cis structure to P2EP. Infrared spectrum confirmed that this polymer is mainly cis. I2 doping of P2EP, P2EPH and P2EPtBu provoked differents perturbations in these systems. P2EP and P2EPH increased their cis segments amount, as showed by infrared spectra. On the other hand, this doping induced the decrease of cis segments in P2EPtBu. Raman spectroscopy supported these conclusions too. This difference was interpretated as due to the side group volume, which could favour one or another isomerization. Electronic conductivity increased by doping for all these three polymers, however the conductivity is considerably lower than the presented by polyacetylene (from 10-5 to 10-7 S cm-1 for these substituted polyacetylenes against typically 10-2 S cm-1 for polyacetylene).

Neste trabalho foram sintetizados a poli-2-etinilpiridina (P2EP), poli-4etinilpiridina (P4EP), o iodeto de poli(2-N-t-butilpiridiniumilacetileno) (P2EPtBu) e a poli-β-etinilnaftaleno (P&#946:EN), os quais são poliacetilenos substituídos. Estes polímeros, juntamente com o cloreto de poli(2-piridínio-2-piridilacetileno) (P2EPH), foram caracterizados por espectroscopia vibracional no infravermelho e Raman. Estes polímeros apresentaram variação na posição das bandas Raman com a energia da radiação excitante - chamada dispersão Raman ou fotosseletividade - da ordem de 10 cm-1, bem inferior ao apresentado pelo poliacetileno (cerca de 60 cm-1). Este deslocamento foi interpretado utilizando-se dois dos modelos existentes para descrever este fenômeno: o Modelo de Modo de Amplitude (AMM) e o Modelo de Coordenada de Conjugação Efetiva (ECCM), os quais fornecem informações sobre a estrutura polimérica e sobre seus níveis eletrônicos. Utilizando-se o AMM foi possível obter informações sobre os níveis eletrônicos excitados de mesma simetria que o estado eletrônico fundamental. Por outro lado, o ECCM, com a ajuda de cálculos DFT, mostrou diferenças na extensão da conjugação e no grau de dimerização entre o P2EP na forma cis e trans e indicou que este polímero apresentava, predominantemente, a estrutura cis, fato este confirmado pelos espectros no infravermelho. A dopagem com I2 provocou efeitos diferentes na estrutura dos polímeros. Os espectros no infravermelho dos polímeros dopados indicaram que o P2EP e o P2EPH apresentaram aumento na quantidade de segmentos cis enquanto o P2EPtBu apresentou diminuição na quantidade desses segmentos. Os espectros Raman dos polímeros dopados confirmaram os dados dos espectros no infravermelho. Esta diferença foi interpretada como sendo devida à diferença no volume do substituinte, pois grupos volumosos favorecem o isômero trans-cisóide onde a distância entre os substituintes é maior. A dopagem também levou a um aumento na condutividade dos polímeros, porém os valores de condutividade obtidos foram bem inferiores que os apresentados pelo poliacetileno dopado (10-5 a 10-7 contra 102 S cm-1, tipicamente). / In this work poly-2-ethynylpyridine (P2EP), poly-4-ethynylpyridine (P2EP), poly(2-N-t-butylpyridiniumylacetylene) iodide (P2EPtBu) and poly-β-pethynylnaphthalene (PβEN) were synthesized. These polymers, which are polyacetylenes derivatives, and poly(2-pyridinium-2-pyridylacetylene) chloride (P2EPH) were analyzed by vibrational techniques. It was observed that some Raman band positions shifted when excitation energy was changed, however this shift (ca. 10 cm-1) was substantially smaller than the presented by polyacetylene (ca. 60 cm-1). This phenomenon called Raman dispersion or photoselectivity - was examinated using Amplitude Mode Model (AMM) and Effective Conjugation Coordinate Model (ECCM) and gave us important information about polymer structure and its electronic energies levels. Using AMM it was possible to infere information about electronic excited states with the same symmetry of the electronic ground state. ECCM gave us evidences about differences in conjugation extension and dimerization degree between cis and trans structures. Besides this, ECCM was successful in indicating cis structure to P2EP. Infrared spectrum confirmed that this polymer is mainly cis. I2 doping of P2EP, P2EPH and P2EPtBu provoked differents perturbations in these systems. P2EP and P2EPH increased their cis segments amount, as showed by infrared spectra. On the other hand, this doping induced the decrease of cis segments in P2EPtBu. Raman spectroscopy supported these conclusions too. This difference was interpretated as due to the side group volume, which could favour one or another isomerization. Electronic conductivity increased by doping for all these three polymers, however the conductivity is considerably lower than the presented by polyacetylene (from 10-5 to 10-7 S cm-1 for these substituted polyacetylenes against typically 10-2 S cm-1 for polyacetylene).

Simulações de sais de carbonato fundidos pelo método de Dinâmica Molecular (MD) foram efetuadas com o modelo polarizável de cargas flutuantes (FC). O modelo de cargas flutuantes implementa os efeitos de polarização pelo método de Lagrangiano estendido, onde as variáveis extras são as próprias cargas parciais do íon poliatômico. O modelo FC foi parametrizado por meio de cálculos ab inito, aplicado ao ânion carbonato. Cálculos de Química Quântica ab initio foram utilizados para corroborar o modelo proposto para o ânion carbonato. Os sistemas investigados consistem em misturas de carbonatos alcalinos fundidos, Li2CO3/K2CO3, os quais são utilizados como eletrólitos em células a combustível. As simulações MD foram utilizadas para verificar o efeito da polarização dos ânions sobre a estrutura e dinâmica do líquido. Estudamos o efeito da inclusão de polarização sobre a condutividade do eletrólito. / Simulations of molten carbonate salts by the method of Molecular Dynamics (MD) have been performed with the fluctuanting charge (FC) model. The FC model implements the effect of polarization by using method of extended Lagrangian, where the extra variables are the partial charges of the poliatomic ion. The FC model was parametrized by ab inito calculations os a single carbonate anion. Quantum Chemistry calculations bave been used to corroborate the model for the carbonate anion. The investigated systems consist of alkaline carbonate mixtures, Li2CO3/K2CO3, which are used as electrolytes in fuel cells. MD simulations have been used to verify polarization effects on structure and dynamics of the liquid. We study the effect of including anion polarization on the condutivity of the electrolyte.