Distribuição de peso molecular em poliolefinas: contribuição ao estudo da bimodalidade e da desconvolução

Molecular weight distribution of polyolefins: deconvolution and bimodality

ARTIGO TÉCNICO CIENTÍFICO

Distribuição de peso molecular em poliolefinas: contribuição ao estudo da bimodalidade e da desconvolução

Molecular weight distribution of polyolefins: deconvolution and bimodality

Montserrat Fortuny; Márcio Nele; Príamo A. Melo; José Carlos Pinto

Programa de Engenharia Química, COPPE, UFRJ

^{Endereço para correspondência}

RESUMO

A distribuição de pesos moleculares de resinas poliolefínicas é freqüentemente analisada em termos da existência de bimodalidades e com auxílio de técnicas de desconvolução. As condições necessárias para a obtenção de uma resina poliolefínica bimodal são apresentadas em função das características cinéticas dos sítios ativos presentes no catalisador, sendo que um certo grau de bimodalidade sempre pode ser definido para essas resinas. Distribuições de peso molecular de Schulz-Flory dinâmicas são neste trabalho derivadas para a utilização em problemas de desconvolução, quando a dinâmica da concentração de monômero não pode ser desprezada. É demonstrado que a utilização de distribuições de Schulz-Flory padrão, nestas condições, leva à detecção de sítios ativos inexistentes no sistema catalítico.

Palavras-chave: Desconvolução, modelagem, distribuição de peso molecular, polimerização.

ABSTRACT

The deconvolution of molecular weight distributions (MWD) may be useful to obtain information about the polymerization kinetics and properties of catalytic systems. However, deconvolution techniques are normally based on steady-state assumptions and very little has been reported about the use of non-stationary approaches for deconvolution of MWDs. In spite of that, polymerization reactions are often performed in batch or semi-batch modes. For this reason, dynamic solutions are proposed here for simple kinetic models and are then used for deconvolution of actual MWD data. Deconvolution results obtained with dynamic models are compared to deconvolution results obtained with the standard stationary Flory-Schulz distributions. For coordination polymerizations, results show that dynamic MWD models are able to describe experimental data with fewer catalytic sites, which indicates that the proper interpretation of the reaction dynamics may be of fundamental importance for kinetic characterization. On the other hand, reaction dynamics induced by modification of chain transfer agent concentrations seems to play a minor role on the shape of the MWD in free-radical polymerizations

Keywords: Deconvolution, modeling, molecular weight distribution, polymerization.

Introdução

As propriedades finais dos polímeros são fortemente dependentes dos seus pesos moleculares médios, da largura da distribuição de peso molecular (normalmente expressa através da polidispersão, ou MFR) e também da forma desta distribuição^[1]. O controle da forma da distribuição de peso molecular (DPM) é efetuado na produção de resinas bimodais pela indústria de poliolefinas. Estas resinas apresentam excelente processabilidade, quando comparadas às resinas lineares com as mesmas propriedades mecânicas, devido à presença simultânea de frações de pesos moleculares bastante distintas misturadas em nível molecular. A fração de menor peso molecular é responsável pela boa processabilidade da resina, enquanto a fração de maior peso molecular proporciona boas propriedades mecânicas^[2-4]. As principais técnicas utilizadas para produzir estas resinas incluem a produção em reatores sequenciais, alimentação de mistura de catalisadores e a imposição de variações súbitas das condições de reação^[5].

A distribuição de peso molecular das poliolefinas produzidas por processos catalíticos pode ser descrita como uma soma de distribuições de Schulz-Flory^[6], e um arcabouço teórico baseado nestas distribuições poder ser utilizado para a obtenção de informações cinéticas sobre os catalisadores, a partir da desconvolução das curvas de distribuição de peso molecular^[7].

Neste trabalho, faz-se a expansão do arcabouço teórico necessário para representar curvas de peso molecular descritas por somas de distribuições de Schulz-Flory dinâmicas, obtidas quando a dinâmica de concentração do monômero não pode ser desprezada. Será apresentada uma caracterização matemática da bimodalidade das curvas de distribuição de peso molecular descritas por estas distribuições, que mostra quando uma resina bimodal pode ser obtida, definindo-se também um grau de bimodalidade para as distribuições de peso molecular. As distribuições de Schulz-Flory podem ser generalizadas para distribuições de Schulz-Flory dinâmicas, que podem ser então utilizadas em problemas de desconvolução nas quais a concentração de monômero muda durante a polimerização. Nestes casos, um número de sítios ativos maior do que o real pode ser encontrado se distribuições de Schulz-Flory estacionárias são utilizadas como base para o procedimento de desconvolução.

Este manuscrito apresenta alguns resultados obtidos no campo das distribuições de peso molecular de poliolefinas. Um relato mais amplo e detalhado destes resultados encontra-se disponível na literatura, incluindo outros tipos de polimerização^[8], que não Ziegler-Natta, caracterização de multimodalidade e estudo de outras distribuições dinâmicas^[9].

Resultados e Discussão

Bimodalidade em Poliolefinas

A curva de distribuição de peso molecular de polímeros pode ser interpretada como uma distribuição estatística. Como tal apresenta algumas características que são brevemente discutidas a seguir. Uma distribuição P_i no intervalo discreto [1,µ) pode ser definida como uma função de i tal que:

Uma moda da distribuição P_i é um máximo local tal que:

A curva de distribuição é dita unimodal se possuir apenas uma moda, bimodal, se possuir duas modas e multimodal se possuir mais de duas modas.

A distribuição de Schulz-Flory pode ser definida pela relação de recorrência:

onde q é denominado probabilidade de propagação e está relacionado com as condições de polimerização na forma:

quando o crescimento da cadeia é limitado predominantemente por reações de transferência de cadeia. M é a concentração de monômero, e X_m, m=1...NX é a concentração dos agentes de transferência de cadeia. A Equação (3) pode ser reescrita em termos da espécie P₁ na forma:

Utilizando-se a Equação (1), a Equação (5) pode ser normalizada em termos de q na forma:

que é a forma usual da distribuição de Schulz-Flory em número.

Se a curva de distribuição de peso molecular numérica é formada por NS distribuições de Flory, ela pode ser escrita na forma:

onde a_j são escalares tais que:

e q_j é a probabilidade de propagação da distribuição j. NS é o número de frações de polímero com distribuições de Schulz-Flory, possivelmente produzidos em NS sítios diferentes, que constituem a DPM do polímero final.

A distribuição NS-Flory representada pela Equação (7) é monotônica, em relação a i, com um máximo dentro do intervalo i Î [1,¥) em i =1. Então a distribuição numérica de peso molecular do tipo NS-Flory é sempre unimodal, com moda em i =1. A Figura 1 apresenta um exemplo de curva de distribuição de peso molecular numérico de Flory.

A distribuição mássica de peso molecular de Schulz-Flory pode ser definida como:

A distribuição mássica de Flory tem um máximo localizado no intervalo i Î [1,¥) quando q > 0.5. Se q < 0.5, o máximo desta distribuição é localizado em i =1. Em sistemas de polimerização, q é muito próximo de um, de forma que a curva de distribuição de peso molecular mássica tem a forma de sino (Figura 2), sendo crescente no intervalo [0, q/(1q)] e decrescente no intervalo [q/(1-q), ¥)].

Se a distribuição é formada pela soma de NS distribuições de peso molecular mássicas, então ela pode ser escrita na forma:

onde a_j são escalares com as propriedades definidas pela Equação (8).

A distribuição mássica de Flory para NS =2 (dois sítios) pode ser escrita na forma:

Admitindo que q₂ > q₁, é possível demonstrar que a distribuição mássica de Flory para dois sítios é uma função crescente em i no intervalo [0, q/(1q)] e decrescente no intervalo [q/(1-q), ¥)]. Se a distribuição de Flory para dois sítios for aproximada por uma função contínua, os extremos são localizados em pontos tais que:

enquanto os pontos de inflexão devem satisfazer a:

A Teoria da Bifurcação^[10] demonstra que o número de soluções de um sistema algébrico muda nos pontos nos quais o Jacobiano (matriz de derivadas) é singular. Portanto, a condição para que a distribuição mássica de Flory com dois sítios seja bimodal é que as Equações (13) e (15) sejam satisfeitas simultaneamente. Descartando a solução trivial, a solução deste sistema é dada por:

i está no domínio real se:

A Equação (17) é satisfeita se:

que é a relação entre os valores de q₁ e q₂ que deve ser satisfeita para que a distribuição mássica de Flory com dois sítios seja bimodal. A Figura 3 mostra a faixa de valores que q₁ e q₂ que leva a uma distribuição de peso molecular mássica bimodal.

Comparando as distribuições de peso molecular numérico (Equação (6)) e mássico (Equação (9)), é possível observar que a distribuição mássica é uma convolução de primeira ordem da distribuição numérica; ou seja, é o resultado normalizado da multiplicação da distribuição numérica por i¹. Em muitas situações práticas, convoluções da distribuição numérica são encontradas. Por exemplo, curvas de distribuição de peso molecular obtidas por cromatografia de exclusão (Equação (19)) são apresentadas como uma convolução de segunda ordem da distribuição numérica de Schulz-Flory, na forma.

onde F_i é a concentração cumulativa mássica de espécies poliméricas de tamanho i.

A enésima convolução da distribuição de Schulz-Flory pode ser dada por:

onde

e G_N é o enésimo momento da distribuição numérica de Schulz-Flory, dada pela seguinte relação de recorrência:

A convolução de enésima ordem de uma distribuição de Flory para dois sítios tem a seguinte forma:

com 0 < a < 1 e q₂ > q₁. Analogamente, ao desenvolvimento anterior para a distribuição mássica de Flory com dois sítios, a solução dos sistema de equações que define os pontos de bimodalidade é dada por:

i é real se:

que admite as soluções:

As faixas de valores de q₁ e q₂ que levam a convoluções de ordem N bimodais é dada pela Equação (26). É interessante notar que, à medida que a ordem da convolução aumenta, a faixa de valores de q₂ que leva à bimodalidade aumenta. Quando a ordem da convolução N tende a infinito, a convolução da distribuição é bimodal para todos os valores de q₂, uma vez que.

A Figura 4 mostra as regiões de q₁ e q₂ onde uma distribuição bimodal pode aparecer, para diferentes ordens de convolução. Como algumas distribuições são bimodais até para N=1, e outras somente quando N tende a infinito, é natural definir um grau de bimodalidade (r) para uma distribuição de peso molecular:

onde N_b é a menor ordem de convolução que leva à bimodalidade. Então o grau de bimodalidade é 1, se a distribuição mássica for bimodal, e aproxima-se de zero à medida que ordens superiores de convolução são necessárias para a bimodalidade aparecer.

Desconvolução de Distribuições de Peso Molecular

A desconvolução de curvas experimentais de peso molecular para a obtenção de parâmetros cinéticos pode ser realizada descrevendo-se a curva experimental como uma soma de distribuições de Schulz-Flory^[7]. Uma das principais hipóteses em que se baseia esta análise é a de que as resinas poliméricas são produzidas sob condições estacionárias ou quase-estacionárias. No entanto, diversas polimerizações são conduzidas em condições nas quais a dinâmica é relevante e não pode ser desprezada, como por exemplo nos processos em batelada, e nos sistemas semi-batelada com elevada resistência à transferência de massa. Nestas condições, a curva de peso molecular experimental deveria ser descrita como uma soma de curvas de Schulz-Flory modificadas, que levam em conta explicitamente o comportamento dinâmico do sistema. A derivação destas curvas, denominadas de distribuições de Schulz-Flory dinâmicas, e sua utilização na desconvolução de curvas de pesos moleculares experimentais é apresentada a seguir.

Derivação da Distribuição de Schulz-Flory Dinâmica

As distribuições Schulz-Flory modificadas podem ser derivadas a partir de um esquema cinético simplificado que contém as etapas básicas de um processo linear de polimerização:

Neste esquema estão incluídas as etapas mais importantes para a determinação da distribuição de peso molecular final das cadeias poliméricas. De acordo com este mecanismo e admitindo que o catalisador seja uma mistura de NS sítios ativos com constantes cinéticas diferentes, os balanços de massa do monômero e das cadeias inativas de polímero podem ser representados de acordo com as seguintes equações :

onde R_i^j representa a concentração das cadeias ativas com tamanho i formadas no sítio ativo de tipo j.

O balanço das cadeias ativas pode ser resolvido admitindo-se condições estacionárias ou quase-estacionárias para estas espécies (Equação 34), uma vez que suas concentrações variam desprezivelmente no decorrer da polimerização na maioria dos sistemas de interesse prático. Entretanto, para as cadeias poliméricas mortas, cuja dinâmica não é desprezível, a concentração pode ser obtida como função da conversão do monômero (X). Introduzindo a Equação (34) nas Equações (32) e (33), obtém-se então a Equação (35):

onde R_o^j representa a concentração inicial dos sítios ativos de tipo j.

A Equação (35) pode ser simplificada, admitindo-se que as constantes de desativação dos sítios ativos são similares para todos os sítios,^[7] e então resolvida analiticamente. A expressão para a DPM como uma convolução de segunda ordem, nesse caso, pode ser dada pela Equação (36):

onde Aj, q_j, q_oj e a_Rj são definidos na forma:

e M_o é concentração molar inicial de monômero.

Analisando-se a Equação (36), é possível observar que a DPM é função da conversão do monômero (através de q_j), ao contrário da distribuição de Schulz-Flory, na qual a DPM é função somente das características cinéticas do sítio ativo e das condições de polimerização. É possível observar também que as distribuições dinâmicas podem possuir índices de polidispersão superiores a 2; ou seja, índices de polidispersão maiores que aqueles obtidos para DPM produzidas sob condições estacionárias.

A Figura 5 mostra a utilização da distribuição de Schulz-Flory estacionária [Equação (19)] para a desconvolução de uma distribuição de peso molecular cuja dinâmica de polimerização não é desprezível [Equação (36)]. Apesar da curva dinâmica ter sido gerada a partir de um único sítio, são necessários duas curvas de Schulz-Flory estacionárias (dois sítios) para descrevê-la. Portanto, se a dinâmica não pode ser desprezada, pode-se obter um número de sítios maior do que o verdadeiro, se a dinâmica de polimerização não é desprezível, quando curvas de Schulz-Flory estacionárias são utilizadas na descrição da DPM.

É interessante comparar os valores de q₁ e q₂ que levam à produção de polímeros com distribuições de pesos moleculares mássicas bimodais em condições estacionárias e dinâmicas. A Figura 6 mostra que a faixa de valores de q₁ e q₂ que levam à distribuições unimodais é maior para distribuições obtidas em condições dinâmicas do que para aquelas obtidas em condições estacionárias. Este comportamento pode ser atribuído ao fato que nas distribuições dinâmicas os valores de q₁ e q₂ seguem trajetórias que aproximam os pesos moleculares dos polímeros produzidos por estes sítios, enquanto na polimerização estacionária a diferença na distribuição de peso molecular dos polímeros produzidos por cada sítio é constante durante a polimerização, não ocorrendo esta aproximação.

Desconvolução de curvas de peso molecular experimentais

A distribuição dinâmica [Equação (36)] e a estacionária [Equação (5)] foram utilizadas para a desconvolução de curvas de peso molecular experimentais de amostras de polipropileno produzidas a partir de um processos contínuo de polimerização em lama empregando um catalisador Ziegler-Natta convencional^[7], e em solução, utilizando um catalisador metalocênico^[11].

Os catalisadores metalocênicos destacam-se dos sistemas Ziegler-Natta convencionais por usualmente apresentarem um sítio ativo definido^[12]. Entretanto, se a polimerização não for realizada em condições controladas que minimizem efeitos de transferência de massa, procedimentos de desconvolução utilizando a distribuição de Schulz-Flory estacionária podem levar à detecção de um sítio adicional falso, como sugerido na Figura 5.

A Figura 7 mostra a desconvolução, utilizando distribuições de Flory estacionárias, da distribuição de peso molecular de uma amostra de PP metalocênico na qual foram necessários dois sítios ativos para a descrição apropriada da curva, sugerindo a existência de dois sítios ativos. Entretanto, se a distribuição dinâmica de Flory for utilzada na desconvolução, verifica-se que apenas um sítio é de fato necessário para descrever a distribuição (Figura 8), como esperado pela natureza dos catalisadores metalocênicos. O aparecimento de um sítio ativo adicional falso pode ser explicado pela variação na concentração de monômero no meio reacional, devido ao aumento da viscosidade no decorrer da polimerização. Esse efeito viscoso, muito comum em polimerizações em solução, pode alterar bastante a interpretação final dos resultados obtidos.

Ao contrário dos reatores de bancada, nos quais os efeitos dinâmicos na distribuição de peso molecular devem ser cuidadosamente controlados, reatores industriais de polimerização de olefinas operam em estado estacionário e a distribuição de peso molecular deve ser adequadamente descrita por distribuições de Schulz-Flory estacionárias. As Figuras 9 e 10 mostram a desconvolução da distribuição de peso molecular de uma amostra de polipropileno produzida em um reator de escala industrial. As deconvoluções com distribuições de Schulz-Flory dinâmicas e estacionárias levam a resultados similares. Particularmente, obtém-se o mesmo número de sítios ativos, o que enfatiza que a distribuição dinâmica converge para a estacionária quando a dinâmica de polimerização é desprezível. Este resultado demonstra a robustez da utilização das distribuições dinâmicas para verificação da existência de efeitos significativos da dinâmica de polimerização na distribuição de peso molecular do polímero.

Conclusões

As distribuições de peso molecular das resinas poliolefínicas podem ser descritas por somas de distribuições de Schulz-Flory. Se a distribuição de peso molecular de uma resina for descrita como uma soma de duas distribuições de Schulz-Flory, a distribuição de peso molecular mássica será bimodal se as probabilidades de propagação dos sítios obedecerem a relação ou . Foi demonstrado que existe sempre uma convolução desta distribuição que é bimodal, e o inverso da ordem de convolução cuja bimodalidade (N_b) aparece, foi definido como grau de bimodalidade (r = 1 / N_b).

A distribuição de peso molecular de resinas poliolefínicas preparadas em reações em que ocorrem variações significativas das condições de polimerização não pode ser descrita como uma soma de distribuições de Schulz-Flory estacionárias. Caso contrário, as conclusões obtidas da desconvolução podem ser errôneas. Para estes casos, distribuições Schulz-Flory dinâmicas derivadas neste trabalho devem ser utilizadas.

A distribuição de Flory dinâmica é capaz de detectar corretamente o número de sítios ativos em uma polimerização com variações de concentração de monômero empregando catalisador metalocênico (sítio único), enquanto a distribuição de Schulz-Flory estacionária detecta a presença de um sítio espúrio adicional. No caso da desconvolução da DPM de uma amostra produzida em condições estacionárias, os resultados da distribuição dinâmica convergiram para os da distribuição estacionária, a correta para este caso, demonstrando a robustez da utilização destas distribuições.

Enviado: 15/12/03

Aprovado: 09/06/04

• 1. Hinchliffe, M.; Montague, G.; Willis, M. - AIChE J., 49, p.2609 (2003).
• 2. Thomas, E. L. - "Structure and Properties of Polymers", VCH, Weinheim (1993).
• 3. Bicerano, B. - "Prediction of Polymer Properties", Marcel Dekker, 2nd Ed., New York, (1996).
• 4. Dan, N.; Tirrell, M. - Macromolecules, 26, p.6467, (1993).
• 5. Nele, M.; Pinto, J. C. - J. Appl. Polym. Sci., , 77, p.437 (2000).
• 6. Nele, M.; Pinto, J. C. - Macromol. Theory and Simul. 11, p.293, (2002).
• 7. Matos, V.; Mattos Neto, A.G.; Pinto, J.C. - J. Appl. Polym. Sci., , 79, p.2076 (2001).
• 8. Pinto, J. C. - Macromol. Theory Simul., 10, p.79 (2001).
• 9. Fortuny, M.; Nele, M.; Melo, P. A.; Pinto, J. C. - Macromol. Theory and Simul., submetido (2003).
• 10. Iooss, G.; Joseph, D. D. - "Elementary Stability and Bifurcation Theory", Springer, New York (1980).
• 11. Lovisi, H. - "Desenvolvimento de copolímeros de propileno elastoméricos via catálisadores metalocênicos isoespecíficos", Tese de Doutorado, Universidade Federal do Rio de Janeiro, Brasil (2001).
• 12. Brintzinger, H. H.; Fischer, D.; Mulhaupt, R.; Rieger, B.; Waymouth, R. M. - Angew. Chem. Int. Ed. Engl., 34, p.1143 (1995).

Datas de Publicação

• Publicação nesta coleção
  29 Set 2004
• Data do Fascículo
  Set 2004

Histórico

• Aceito
  09 Jun 2004
• Recebido
  15 Dez 2003