JBCS

INTRODUÇÃO

Os cátions metálicos são parte essencial dos fertilizantes, pois eles são responsáveis por várias funções nas plantas, que vão desde a participação na síntese de proteínas, até sua presença em enzimas que participam de processos da fotossíntese.^1-2 Fertilizantes ou adubos são compostos químicos inorgânicos ou orgânicos utilizados para propiciar os nutrientes essenciais ao desenvolvimento das plantas, e são aplicados na agricultura para um maior rendimento da produção.³ No entanto, atualmente, estudos^4-5 demonstram que para a maioria dos fertilizantes usados, parte dos seus nutrientes são perdidos no ambiente, podendo ocasionar uma série de contaminações ambientais. Por exemplo, os fertilizantes que possuem em sua composição nitratos e fosfatos podem perder parte desses compostos em processos de lixiviação, propiciando o processo de eutrofização em águas de mananciais e lençóis freáticos. Ainda, os fertilizantes quando não aplicados de forma correta podem proporcionar desequilíbrios no ecossistema e aumento do número das fontes de toxicidade aguda e carcinogênicas para os seres humanos.^4-5 Desta forma, o uso adequado dos fertilizantes é uma atividade que deve ser levada em consideração, visando à sustentabilidade econômica, ambiental e social do agronegócio mundial.

Nos últimos anos a tecnologia de liberação controlada/modificada tem ganhado destaque, pois ela emergiu como uma importante alternativa capaz de resolver alguns dos problemas que acompanham o uso de substâncias ativas de interesse em alguns setores,⁶ como por exemplo o agrícola (pesticidas e herbicidas),^7-11 alimentício,^12-13 o farmacêutico,^14-17 e, ineditamente, poderá ser aplicada para os cátions metálicos essenciais às plantas visando o aumento da eficiência dos fertilizantes, ou seja, aumento da disponibilidade destes no solo, minimização dos impactos ambientais e garantia de um desenvolvimento mundial considerado sustentável. Esta tecnologia propicia a liberação dos princípios ativos no decorrer de um determinado tempo através de processos físico-químicos, tais como dessorção e/ou erosão, sendo capaz de permitir a alteração das propriedades físico-químicas dos ativos para a liberação controlada/adequada destes,^6,18 possibilitando a atuação de forma contínua e prolongada dos cátions metálicos.

Existem vários parâmetros utilizados para a avaliação da eficiência destes sistemas. Um deles é o cálculo da taxa de encapsulamento dos nutrientes na matriz polimérica, como também a concentração de nutriente que este sistema libera no decorrer de um determinado tempo.^6-18 Assim, são requeridas técnicas analíticas para a quantificação da concentração dos cátions metálicos nesses processos.

A determinação de cátions metálicos, nas mais diferentes matrizes, é tradicionalmente realizada utilizando-se técnicas espectrométricas, em função de sua exatidão, precisão e relativa simplicidade de manuseio. No entanto, amostras que são analisadas por tais técnicas requerem um preparo adequado para possibilitar a detecção e eliminação de possíveis interferentes.¹⁹ Portanto, metodologias para a abertura de amostras dos sistemas de liberação controlada são necessárias e essenciais para a quantificação correta de metais e o não comprometimento dos equipamentos e resultados obtidos.

Atualmente existem vários métodos de preparo para diferentes tipos de amostras. Especialmente para matrizes poliméricas são relatados alguns na literatura.^20-22 Porém, aberturas de amostras de sistemas de liberação controlada poliméricos são inexistentes.

Técnicas de preparo de amostras para análise por absorção atômica devem apresentar simplicidade, rapidez, utilizar pequenos volumes de ácidos, permitir a dissolução de grande número de amostras e, finalmente, produzir resultados precisos e exatos.¹⁹ O sistema tradicional usado é a digestão de amostras com ácido em chapa aquecedora e agitação mecânica, pois é relativamente simples e barato. No entanto, existem muitas variáveis envolvidas neste procedimento que podem interferir diretamente no ótimo do processo, tais como: tipo e volume de ácido empregado, temperatura, tempo de digestão, entre outras. Portanto, estudos que demonstrem a condição ótima para a abertura das partículas poliméricas são fundamentais.

Com isso, torna-se muito útil e mesmo necessário o uso de experimentos estatisticamente planejados de forma multivariada para a investigação do real efeito das possíveis variáveis, assim como das interações entre elas, sobre o sistema acima referenciado. Quando se trabalha com um sistema complexo em que muitas variáveis são capazes de influenciar, porém não se sabe em qual magnitude, pode-se fazer uso como primeira etapa do processo de investigação/otimização de um planejamento de triagem. Fatores significativos selecionados após execução de um planejamento adequado de triagem podem ser empregados em uma metodologia de superfície de resposta (RSM), para a otimização do experimento.^23-25

Uma vez que na literatura é inexistente o relato de condições ótimas para abertura de amostras de partículas poliméricas de processos de liberação controlada carreadoras de cátions metálicos, este trabalho teve como objetivo desenvolver uma metodologia, de fácil aplicação e baixo gasto de reagentes, para a digestão destas matrizes e posterior quantificação de íons metálicos por espectroscopia de absorção atômica por chama (FAAS), com o auxílio de métodos experimentais de otimização multivariada.

PARTE EXPERIMENTAL

Reagentes, soluções e amostras

Todas as soluções foram preparadas com reagentes de grau analítico e os volumes finais ajustados com água deionizada do sistema Milli-Q^® (Millipore Corporation). Todas as vidrarias foram previamente lavadas em banhos contendo ácido nítrico (HNO₃) (Synth) 10% v/v por 24 h. As soluções de referência de carbono total e inorgânico foram preparadas pela dissolução de 2,1250 g de hidrogenoftalato de potássio (C₈H₅KO₄) (Sigma Aldrich) em 1 L de água deionizada para a primeira, e dissolução de 3,4970 g de hidrogenocarbonato de sódio (NaHCO₃) (Sigma Aldrich) em 1 L de água deionizada para a segunda. Já as soluções dos cátions metálicos de cobre, ferro, magnésio, manganês, potássio e zinco foram preparadas com diluições seriais de padrões de 999 ± 4 mg L^-1 (Sigma Aldrich). Todas as soluções foram armazenadas e mantidas refrigeradas a 8 ºC para preservação.

O preparo dos sistemas poliméricos foi realizado segundo o método de nanoprecipitação, que consiste na precipitação ou deposição interfacial de um polímero pré-formado na interface de uma emulsão óleo/água (O/A), promovendo a formação dos sistemas carreadores.²⁶ Inicialmente, pesou-se aproximadamente 5,0000 g de álcool polivinílico (PVA) Sigma Aldrich e procedeu-se a sua solubilização em água deionizada a 70 ºC sob agitação magnética com uso de agitador magnético Hot Lab II - Nalgon, para o preparo de uma solução de PVA 5% m/v. Posteriormente, realizou-se diluições das soluções de referência dos cátions metálicos. Alíquotas destas soluções foram adicionadas a 100 mL da solução de PVA.

Para o preparo da solução orgânica foram pesados aproximadamente 600,0 mg de poli(ε-caprolactona) - Sigma Aldrich - e 300,0 mg de um polímero obtido de rejeito industrial (70% m/m) - Verti Ecotecnologias UFMG, que foram solubilizados, sob agitação mecânica e temperatura de 70 ºC, em uma mistura de água-clorofórmio (Sigma Aldrich) na proporção 1:3. Em seguida, essa solução foi adicionada à solução de PVA com os cátions metálicos, onde foi mantida sob agitação e temperatura de 70 ºC por 4 h. A emulsão foi mantida em capela por 24 h para a evaporação do solvente orgânico. No fim deste processo, completou-se o volume para 200 mL com água deionizada, onde metade desta suspensão foi filtrada com papel filtro quantitativo JP 40 faixa branca. Os materiais retidos no papel de filtro do procedimento citado foram armazenados e mantidos em dessecador. Após completamente seca, certa quantidade de amostra foi coletada e posta sob uma fita de carbono. As amostras foram levadas a um metalizador Q150R S Aname para a deposição de ouro (100 a 200 nm) sob vácuo. Em seguida, a análise da morfologia das amostras por meio das partículas metalizadas foi feita utilizando um microscópio eletrônico de varredura (MEV), Jeol 1200.

Para a digestão das amostras foi utilizado um agitador magnético Hot Lab II NalgonAs, HNO₃ concentrado e termômetro químico (-10 +110 : 1 ºC - 260 mm, Incoterm) que foi preso ao sistema por um suporte universal conectado à chapa. Posteriormente, as soluções de cada ensaio dos planejamentos experimentais utilizados, conforme as designações descritas mais adiante, foram filtradas em sistema de filtração acoplado a uma bomba a vácuo NOF-650 - Max Pump e membrana de filtração 0,45 µm - Merck, para a posterior análise da concentração de carbono orgânico total (TC) e carbono inorgânico (IC).

Instrumentação

Para a avaliação da digestão das amostras poliméricas utilizou-se um analisador de carbono orgânico total TOC-L - Shimadzu, controlado pelo software TOC Control-V Shimadzu^®, com gás de arraste oxigênio (O₂) com grau de pureza analítica (> 99,5%) (White Martins) e os dados foram analisados através do software TOC Control-V Shimadzu. Foram analisados os teores de carbono orgânico, inorgânico e total residuais nos ensaios gerados pelos planejamentos experimentais multivariados empregados no desenvolvimento do método, apresentados a seguir.

Otimização multivariada

Triagem

A definição dos níveis nos quais as variáveis foram estudadas na triagem se deu a partir de testes preliminares em laboratório e avaliação da literatura disponível para as amostras mais similares às deste estudo, a fim de verificar as faixas operacionais de trabalho.

Após definidas as variáveis e os respectivos níveis a serem estudados na triagem, realizou-se um planejamento fatorial completo 2⁴ com triplicata no ponto central, onde se identificou as variáveis que realmente possuem efeito significativo sobre o sistema de digestão das matrizes poliméricas, para um estudo mais detalhado na metodologia de superfície de resposta. Assim, as variáveis independentes investigadas foram volume de amostra, volume de ácido nítrico, temperatura e tempo de digestão. A Tabela 1 apresenta as variáveis e os respectivos níveis avaliados nesse processo.

Utilizou-se planilhas eletrônicas²³ para gerar a codificação dos experimentos a serem realizados, levando-se em consideração as combinações entre os níveis propostos, e para o tratamento estatístico dos resultados.

Para a análise de TOC uma curva analítica padrão foi construída. Esta foi preparada com a diluição serial de soluções padrões de carbono total (1000 mg L^-1) e carbono inorgânico (1000 mg L^-1). Tal diluição foi realizada pelo próprio equipamento, o analisador de carbono orgânico total. Todos os ensaios realizados no planejamento de triagem foram aferidos para 100 mL. Em seguida todas as soluções foram diluídas 100 vezes e analisadas no analisador de carbono. O analisador de carbono orgânico total possibilita a resposta de três valores de concentração: carbono total, carbono orgânico total e carbono inorgânico, onde o TOC é diferença entre TC e IC.

Com isso foram selecionadas as condições que geravam os menores valores de sinal analítico em termos de concentração de carbono, tendo sido possível identificar as variáveis significativas e os níveis que foram estudados na metodologia de superfície de resposta.

Metodologia de superfície de resposta

O planejamento composto central (CCD) foi utilizado neste estudo para a etapa de otimização do processo de digestão das micropartículas. As variáveis e os níveis avaliados estão descritos na Tabela 2.

Utilizou-se planilhas eletrônicas²³ para a geração dos ensaios do planejamento composto central e tratamento dos dados. Todos os ensaios foram analisados pelo analisador de carbono orgânico total conforme descrito na etapa de triagem. Com este estudo foi possível definir a condição ótima de digestão das matrizes poliméricas de sistemas microestruturados carreadores de cátions metálicos para liberação controlada, como se segue.

RESULTADOS E DISCUSSÃO

Matriz polimérica para liberação controlada de cátions metálicos

As análises de MEV dos materiais obtidos a partir da nanoprecipitação possibilitaram observar a formação de micropartículas esféricas e rugosas com o diâmetro médio de 3,8 µm, classificando-as como micropartículas e sem a presença de formação de agregados. Para posteriores determinações dos teores dos elementos metálicos encapsulados é necessário o estabelecimento de um sistema de digestão eficiente das micropartículas, o que será apresentado nos itens a seguir.

Otimização do sistema

Os resultados da triagem mostraram que dentre todas as variáveis estudadas, somente o volume de amostra teve efeito significativo sobre o sistema de avaliação da concentração de TOC (p = 0,0080) e de TC (p = 0,0037), diminuindo o teor de carbono das amostras com o menor nível estudado (10,00 mL) - efeito negativo. Para o IC, nenhuma variável foi apontada como significativa. Os valores de concentração obtidos estão apresentados na Tabela 3.

Analisando a influência do V_am, pôde-se observar que esta variável foi significativa para as concentrações de TC e TOC, com as melhores respostas em 10,00 mL. Devido à matriz polimérica possuir caráter hidrofóbico (PLC), a sua dissolução em meio ácido pode se tornar onerosa pelo aparecimento de resíduos sólidos sobre a superfície do digerido, ou ainda o material pode ser só parcialmente digerido, formando uma escória sobre as paredes das vidrarias utilizadas¹⁹. Assim, uma maior quantidade de agente oxidante (ácido nítrico neste trabalho) torna-se necessário para a dissolução dos polímeros. Dessa forma, pode-se entender que a diminuição do volume de amostra observada na triagem está diretamente relacionada com concentração de carbono total no processo de digestão desta.

No entanto, ao avaliar-se o gráfico de probabilidade normal para os valores de carbono total (Figura 1) observou-se que os efeitos das variáveis 3 e 4 (temperatura e tempo de digestão, respectivamente) mostraram uma possível tendência de não distribuição normal em relação às demais variáveis e suas interações. Observou-se um comportamento similar para os valores de carbono orgânico total com a adição de um efeito possivelmente não distribuído normalmente pela interação das variáveis 1, 2 e 4 (volume de amostra, volume de ácido e temperatura, respectivamente), como ilustrado na Figura 2. O mesmo não aconteceu para a concentração de carbono inorgânico, já que neste caso as variáveis mostraram-se também não significativas no gráfico de probabilidade normal.

Figura 1. Gráfico de probabilidade normal para a concentração de carbono total em relação às variáveis estudadas na digestão de micropartículas poliméricas

Figura 2. Gráfico de probabilidade normal para a concentração de carbono orgânico total em relação às variáveis estudadas na digestão de micropartículas poliméricas

Em regra, preparos de amostras que empregam a dissolução de matrizes com carga orgânica utilizam maiores volumes de ácidos concentrados. O ácido nítrico possui caráter oxidante e, por este fato, é comumente utilizado em digestões via úmida. Esta característica permite a quebra das ligações entre carbono e hidrogênio na maioria dos casos, o que origina uma alta eficiência da degradação da matéria orgânica¹⁹. Portanto, torna-se necessário aumentar o tempo de reação entre o ácido e a matriz; explicação para o fato das interações entre as variáveis T e t para TOC e V_am, V_Ac e t para TC terem mostrado tendência em serem significativas ao sistema, uma vez que têm-se conhecimento de que maiores valores de ácido e de temperatura permitem maior eficiência da digestão polimérica.

Com os resultados obtidos foi possível obter uma modelagem matemática significativa para o TC, com p = 0,0144 para o modelo e p = 0,2670 para a falta de ajuste. Já para o TOC, obteve-se p = 0,0156 para o modelo e p = 0,5330 para a falta de ajuste. Os coeficientes de determinação (R²) para TC e TOC foram 0,9762 e 0,9752, respectivamente. O mesmo não foi observado para as respostas de carbono inorgânico, o que pode ser justificado pela insignificância dos valores dessa resposta frente às concentrações de TC e TOC. E, em virtude disto, adotou-se a concentração de carbono total para o estudo da superfície de resposta.

Com a finalidade de aprofundar o estudo da influência das variáveis e otimizar o sistema de digestão, isto é, obter o menor teor de carbono possível nas amostras após a digestão ácida, realizou-se o planejamento composto central esférico. Este consiste em um conjunto de técnicas estatísticas e matemáticas para o desenvolvimento, melhoria e otimização de processos químicos, físicos, biológicos, entre outros.^23-24 Esta técnica avalia mais detalhadamente as variáveis significativas obtidas da triagem, por apresentar um número maior de níveis.

O planejamento composto central é um tipo de planejamento para a construção de superfícies de resposta que permite a obtenção de um modelo quadrático. Para o sistema de digestão aqui estudado foi estabelecido um modelo correlacionando a concentração de carbono total com as variáveis volume de amostra, volume de ácido, tempo e temperatura de digestão.

Os ensaios para a análise de superfície de resposta foram realizados aleatoriamente, assim como na triagem. Todas as soluções resultantes de cada ensaio foram analisadas através do analisador de ToC como descrito no item Otimização Multivariada - Triagem, e os dados obtidos estão apresentados na Tabela 4.

Os resultados da metodologia de superfície de resposta para a concentração de carbono total (TC) confirmaram que somente o V_am (variável 1) é significativo ao sistema (p = 0,0000), conforme visualizado na triagem, onde constatou-se o decréscimo no teor de carbono das amostras com o menor nível estudado (5,00 mL) - efeito negativo. Os erros associados a cada coeficiente foram avaliados pelo teste t (α = 0,05), e por meio dos resultados da análise da concentração de carbono total foi possível construir um modelo quadrático por regressão linear múltipla (Equação 1). Observou-se também que o sistema possui comportamento não-linear, pois a interação quadrática da variável 2 (V_Ac) foi significativa ao sistema (p = 0,0194). Na Tabela 5 está apresentada a ANOVA para o modelo obtido.

onde TC = concentração de carbono total, V_am = volume de amostra e V_Ac = volume de ácido.

Observa-se pela ANOVA (Tabela 5) que o modelamento matemático proposto pelo planejamento composto central foi uma regressão significativa, ao nível de significância de 0,05, e que não houve falta de ajuste. O gráfico de valores observados versus valores estimados (Figura 3) apresentou-se linear com y = 0,9163x + 0,1583 e R² = 0,9163.

Figura 3. Gráfico de valores observados versus valores estimados para a concentração de carbono total pela metodologia de superfície de resposta

A Figura 4 (gráfico de resíduos) mostra que os resíduos (diferença entre os valores observados e os valores estimados) estão aleatoriamente distribuídos, confirmando que a falta de ajuste não é significativa para o modelo pois, como esperado, a maior parte dos desvios está relacionada ao erro puro, que é o erro experimental.

Figura 4. Gráfico de resíduos para a concentração de carbono total na metodologia de superfície de resposta

A abordagem de otimização multivariada utilizada proporcionou maior média de degradação polimérica, visto que a concentração inicial de carbono total de 4,6154 mg L^-1 foi reduzida a 0,5540 mg L^-1, ou seja, ocorreu a decomposição de 87,80% da matriz carbonácea (Figura 5); a taxa média de degradação previamente aos planejamentos experimentais não ultrapassava 12,00%.

Figura 5. Variação da degradação polimérica entre ensaios sem e com o uso de planejamentos experimentais multivariados. TC = carbono total; TOC = carbono orgânico total; IC = carbono inorgânico e CCD = planejamento experimental composto central

Obteve-se como ótimo de digestão as condições do ensaio 17 do planejamento CCD esférico (Tabela 6). Esta foi a condição que propiciou 87,80% de degradação da matriz polimérica, conforme anteriormente citado.

Os planejamentos experimentais multivariados possibilitaram estabelecer um procedimento que utiliza baixo uso de reagente (ácido) e é relativamente rápido; características necessárias para um processo de preparo de amostras ideal¹⁹, tendo sido possível obter um método de digestão de amostras adequado, simples e de baixo custo para o preparo de micropartículas poliméricas carreadores de cátions metálicos para análise por FAAS.

CONCLUSÃO

Os sistemas de liberação micropolímeros carreadores de cátions metálicos essenciais paras as plantas aqui estudados são inéditos e, portanto, estudos relacionados a eles são inexistentes. Assim, é de suma importância a obtenção de procedimentos com condições ótimas para a avaliação da eficiência desta matriz. A metodologia proposta por este trabalho mostrou-se eficaz para a degradação de matrizes poliméricas de sistemas de liberação controlada carreadoras de cátions metálicos, utilizando um sistema de digestão convencional e simples, por meio do uso de ácido e chapa aquecedora e agitadora. Os planejamentos experimentais multivariados permitiram obter uma condição ótima, onde proporcionou-se a quase completa degradação da matriz orgânica. Tal degradação se mostrou eficiente para garantir a qualidade das análises e o não comprometimento do equipamento, por evitar a formação de depósitos de carbono. Além disso, conseguiu-se alcançar menores valores de volume de reagentes concentrados e de temperatura de digestão para sistemas que ainda não são descritos na literatura.

AGRADECIMENTOS

À Rede Mineira de Química (RQ-MG) (FAPEMIG CEX - RED-00010-14), ao Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Projeto - 479273/2013-8), à Fundação de Amparo à Pesquisa de Minas Gerais (FAPEMIG) e à Pró-Reitoria de Pesquisa e Pós-Graduação (PROPP) da Universidade Federal de Ouro Preto (UFOP) pelo apoio financeiro e bolsas concedidas. Em especial, os autores também expressam seus agradecimentos ao Laboratório de Química Tecnológica e Ambiental (LQTA) / UFOP, sob a coordenação do professor Dr. Sérgio Francisco de Aquino.

REFERÊNCIAS

1. Epstein, E.; Bloom, A. J.; Nutrição mineral de plantas: princípios e perspectivas, 2ª ed., Londrina: Editora Planta, 2006.

2. Fontes, P. C. R.; Hortic. Bras. 2014, 32, 247.

3. Brasil. Decreto nº 4.954 de 14 de janeiro de 2004. Dispõe sobre a inspeção e fiscalização da produção e do comércio de fertilizantes, corretivos, inoculantes ou biofertilizantes destinados à agricultura, e dá outras providências. Ministério da Agricultura, Pecuária e Abastecimento, Rio de Janeiro, 2004. 35 p.

4. Lezcano, M.; Al-Soufi W.; Novo M.; Rodríguez-Núñez E.; Tato J. V.; J. Agric. Food Chem. 2002, 50, 108.

5. Manolikar, M. K.; Sawant, M. R.; Chemosphere, 2003, 51, 811.

6. Grillo, R.; Dissertação de Mestrado, Universidade Estadual de Campinas, Brasil, 2011.

7. Grillo, R.; Pereira A. E.; Nishisaka C. S.; de Lima R.; Oehlke K.; Greiner R.; Fraceto L. F.; J. Hazard. Mater. 2014, 278, 163.

8. Grillo, R.; Pereira dos Santos, N. Z.; Maruyama, C. R.; Rosa, A. H.; de Lima, R.; Fraceto, L. F. J. Hazard. Mater. 2012, 231, 1.

9. Memarizadeh, N.; Ghadamyari, M.; Adeli, M.; Talebi, K.; Ecotoxicol. Environ. Saf. 2014, 107, 77.

10. Sun, C.; Shu, K.; Wang, W.; Ye, Z.; Liu, T.; Gao, Y.; Zheng, H.; He, G.; Yin, Y.; Int. J. Pharm. (Amsterdam, Neth.) 2014, 463, 108.

11. Guo, Y.; Yanga, Q.; Yana, W.; Lia, B.; Qiana, K.; Lia, T.; Xia, W.; Hea, L.; Int. J. Environ. Anal. Chem. 2014, 94, 1001.

12. de Assis, L. M.; Zavareze, E. R.; Prentice-Hernandez, C.; de Souza-Soares, L. A.; Braz. J. Food Technol. 2012, 15, 99.

13. Chopra, M.; Kaur, P.; Bernela, M.; Thankur, R.; Food Control 2014, 37, 158.

14. Natarajan, J. V.; Nugraha, C.; Ng, X. W.; Venkatraman, S.; J. Controlled Release 2014, 193, 122.

15. Li, X.; Yu, Y.; Ji, Q.; Qiu, L.; Nanomedicine-Nanotechnology Biology and Medicine 2015, 11, 175.

16. Pinheiro, A. C.; Bourbon, A. I.; Cerqueira, M. A.; Maricato, É.; Nunes, C.; Coimbra, M. A.; Vicente, A. A.; Carbohydr. Polym. 2015, 115, 1.

17. Cai, K.; He, X.; Song, Z.; Yin, Q.; Zhang, Y.; Uckun, F. M.; Jiang, C.; Cheng, J.; J. Am. Chem. Soc. 2015, 137, 3458.

18. Lobo, F. A.; de Aguirre, C. L.; Silva, M. S.; Grillo, R.; de Melo, N. F. S.; de Oliveira, L. K.; de Morais, L. C.; Campos, V.; Rosa, A. H.; Fraceto, L. F.; Polym. Bull. 2011, 67, 479.

19. Cadore, S.; Matoso, E.; Santos, M. C.; Quim. Nova 2008, 31, 1533.

20. Pereira, J. S. F.; Knorr, C. L.; Pereira, L. S. F.; Moraes, D. P.; Paniz, J. N. G.; Flores, E. M. M.; Knappc, G.; J. Anal. At. Spectrom. 2011, 26, 1849.

21. Sakurai, H.; Noro, J.; Kawase, A.; Fujinami, M.; Oguma, K.; Anal. Sci. 2006, 22, 225.

22. Cho, H. J.; Mysung, S. W.; Bull. Korean Chem. Soc. 2011, 32, 489.

23. Teófilo, R. F.; Ferreira, M. M. C.; Quim. Nova 2006, 29, 338.

24. Neto, B. D. B.; Scarminio, I. S.; Bruns, R. E.; Planejamento e otimização de experimentos. São Paulo: Editora da Unicamp, 1995.

25. Lima, D. R. S.; Baeta, B. E. L.; da Silva, G. A.; Silva, S. Q.; de Aquino, S. F.; Quim. Nova 2014, 37, 827.

26. Souto, E. B.; Severino, P.; Santana, M. H. A.; Polim.: Cienc. Tecnol. 2012, 22, 101.