Estudo de solvatocromismo em líquidos moleculares orgânicos via método seqüencial Monte Carlo/Mecânica Quântica

Utilizou-se o método seqüencial Monte Carlo / Mecânica Quântica para obterem-se os desvios de solvatocromismo e os momentos de dipolo dos sistemas de moléculas orgânicas: Uracil em meio aquoso, -Caroteno em Ácido Oléico, Ácido Ricinoléico em metanol e em Etanol e Ácido Oléico em metanol e em Etanol. As otimizações das geometrias e as distribuições de cargas foram obtidas através da Teoria do Funcional Densidade com o funcional B3LYP e os conjuntos de funções de base 6-31G(d) para todas as moléculas exceto para a água e Uracil, as quais, foram utilizadas o conjunto de funções de base 6-311++G(d,p). No tratamento clássico, Monte Carlo, aplicou-se o algoritmo Metropólis através do programa DICE. A separação de configurações estatisticamente relevantes para os cálculos das propriedades médias foi implementada com a utilização da função de auto-correlação calculada para cada sistema. A função de distribuição radial dos líquidos moleculares foi utilizada para a separação da primeira camada de solvatação, a qual, estabelece a principal interação entre soluto-solvente. As configurações relevantes da primeira camada de solvatação de cada sistema foram submetidas a cálculos quânticos a nível semi-empírico com o método ZINDO/S-CI. Os espectros de absorção foram obtidos para os solutos em fase gasosa e para os sistemas de líquidos moleculares comentados. Os momentos de dipolo elétrico dos mesmos também foram obtidos. Todas as bandas dos espectros de absorção dos sistemas tiveram um desvio para o azul, exceto a segunda banda do sistema de Beta-Caroteno em Ácido Oléico que apresentou um desvio para o vermelho. Os resultados encontrados apresentam-se em excelente concordância com os valores experimentais encontrados na literatura. Todos os sistemas tiveram aumento no momento de dipolo elétrico devido às moléculas dos solventes serem moléculas polares. Os sistemas de ácidos graxos em álcoois apresentaram resultados muito semelhantes, ou seja, os ácidos graxos mencionados possuem comportamentos espectroscópicos semelhantes submetidos aos mesmos solventes. As simulações através do método seqüencial Monte Carlo / Mecânica Quântica estudadas demonstraram que a metodologia é eficaz para a obtenção das propriedades espectroscópicas dos líquidos moleculares analisados.

Estimativas da condutividade térmica dos minerais e rochas e influência de parâmetros térmicos e petrofísicos na resistividade aparente da formação

O presente estudo realiza estimativas da condutividade térmica dos principais minerais formadores de rochas, bem como estimativas da condutividade média da fase sólida de cinco litologias básicas (arenitos, calcários, dolomitos, anidritas e litologias argilosas). Alguns modelos térmicos foram comparados entre si, possibilitando a verificação daquele mais apropriado para representar o agregado de minerais e fluidos que compõem as rochas. Os resultados obtidos podem ser aplicados a modelamentos térmicos os mais variados. A metodologia empregada baseia-se em um algoritmo de regressão não-linear denominado de Busca Aleatória Controlada. O comportamento do algoritmo é avaliado para dados sintéticos antes de ser usado em dados reais. O modelo usado na regressão para obter a condutividade térmica dos minerais é o modelo geométrico médio. O método de regressão, usado em cada subconjunto litológico, forneceu os seguintes valores para a condutividade térmica média da fase sólida: arenitos 5,9 ± 1,33 W/mK, calcários 3.1 ± 0.12 W/mK, dolomitos 4.7 ± 0.56 W/mK, anidritas 6.3 ± 0.27 W/mK e para litologias argilosas 3.4 ± 0.48 W/mK. Na sequência, são fornecidas as bases para o estudo da difusão do calor em coordenadas cilíndricas, considerando o efeito de invasão do filtrado da lama na formação, através de uma adaptação da simulação de injeção de poços proveniente das teorias relativas à engenharia de reservatório. Com isto, estimam-se os erros relativos sobre a resistividade aparente assumindo como referência a temperatura original da formação. Nesta etapa do trabalho, faz-se uso do método de diferenças finitas para avaliar a distribuição de temperatura poço-formação. A simulação da invasão é realizada, em coordenadas cilíndricas, através da adaptação da equação de Buckley-Leverett em coordenadas cartesianas. Efeitos como o aparecimento do reboco de lama na parede do poço, gravidade e pressão capilar não são levados em consideração. A partir das distribuições de saturação e temperatura, obtém-se a distribuição radial de resistividade, a qual é convolvida com a resposta radial da ferramenta de indução (transmissor-receptor) resultando na resistividade aparente da formação. Admitindo como referência a temperatura original da formação, são obtidos os erros relativos da resistividade aparente. Através da variação de alguns parâmetros, verifica-se que a porosidade e a saturação original da formação podem ser responsáveis por enormes erros na obtenção da resistividade, principalmente se tais "leituras" forem realizadas logo após a perfuração (MWD). A diferença de temperatura entre poço e formação é a principal causadora de tais erros, indicando que em situações onde esta diferença de temperatura seja grande, perfilagens com ferramentas de indução devam ser realizadas de um a dois dias após a perfuração do poço.

Modelamento da permissividade dielétrica de rochas saturadas de óleo e água e suas aplicações em perfilagem de poço

A ferramenta de propagação eletromagnética (EPT) fornece o tempo de propagação (T_pl) e a atenuação (A) de uma onda eletromagnética que se propaga num meio com perdas. Estas respostas da EPT são funções da permissividade dielétrica do meio. Existem vários modelos e fórmulas de misturas sobre a permissividade dielétrica de rochas reservatório que podem ser utilizados na interpretação da ferramenta de alta frequência. No entanto, as fórmulas de mistura não consideram a distribuição e a geometria do espaço poroso, e estes parâmetros são essenciais para que sejam obtidas respostas dielétricas mais próximas de uma rocha real. Foi selecionado um modelo baseado nos parâmetros descritos acima e este foi aplicado à dados dielétricos disponíveis na literatura. Foi obtida uma boa concordância entre as curvas teóricas e os dados experimentais, comprovando assim que a distribuição e a geometria dos poros têm que ser levadas em conta no desenvolvimento de um modelo realista. Foram conseguidas também funções de distribuição de razão de aspecto de poros, através das quais geramos várias curvas relacionando as respostas da EPT com diversas saturações de óleo/gás. Estas curvas foram aplicadas na análise de perfis. Como o modelo selecionado ajusta-se bem aos dados dielétricos disponíveis na literatura, torna-se atraente aplicá-lo à dados experimentais obtidos em rochas de campos brasileiros produtores de hidrocarbonetos para interpretação da EPT corrida em poços destes campos petrolíferos.

Relationships between dendritic morphology, spatial distribution and firing patterns in rat layer 1 neurons

The cortical layer 1 contains mainly small interneurons, which have traditionally been classified according to their axonal morphology. The dendritic morphology of these cells, however, has received little attention and remains ill defined. Very little is known about how the dendritic morphology and spatial distribution of these cells may relate to functional neuronal properties. We used biocytin labeling and whole cell patch clamp recordings, associated with digital reconstruction and quantitative morphological analysis, to assess correlations between dendritic morphology, spatial distribution and membrane properties of rat layer 1 neurons. A total of 106 cells were recorded, labeled and subjected to morphological analysis. Based on the quantitative patterns of their dendritic arbor, cells were divided into four major morphotypes: horizontal, radial, ascendant, and descendant cells. Descendant cells exhibited a highly distinct spatial distribution in relation to other morphotypes, suggesting that they may have a distinct function in these cortical circuits. A significant difference was also found in the distribution of firing patterns between each morphotype and between the neuronal populations of each sublayer. Passive membrane properties were, however, statistically homogeneous among all subgroups. We speculate that the differences observed in active membrane properties might be related to differences in the synaptic input of specific types of afferent fibers and to differences in the computational roles of each morphotype in layer 1 circuits. Our findings provide new insights into dendritic morphology and neuronal spatial distribution in layer 1 circuits, indicating that variations in these properties may be correlated with distinct physiological functions.

Metodologia de discretização espacial para o método radial de interpolação por pontos (RPIM) aplicada para solução numérica das equações de Maxwell

Neste trabalho, foi desenvolvido e implementado um método de discretização espacial baseado na lei de Coulomb para geração de pontos que possam ser usados em métodos meshless para solução das equações de Maxwell. Tal método aplica a lei de Coulomb para gerar o equilíbrio espacial necessário para gerar alta qualidade de discretização espacial para um domínio de análise. Este método é denominado aqui de CLDM (Coulomb Law Discretization Method ) e é aplicado a problemas bidimensionais. Utiliza-se o método RPIM (Radial Point Interpolation Method) com truncagem por UPML (Uniaxial Perfectlly Matched Layers) para solução das equações de Maxwell no domínio do tempo (modo TMz).

An automatic methodology for obtaining optimum shape factors for the radial point interpolation method

In this letter, a methodology is proposed for automatically (and locally) obtaining the shape factor c for the Gaussian basis functions, for each support domain, in order to increase numerical precision and mainly to avoid matrix inversion impossibilities. The concept of calibration function is introduced, which is used for obtaining c. The methodology developed was applied for a 2-D numerical experiment, which results are compared to analytical solution. This comparison revels that the results associated to the developed methodology are very close to the analytical solution for the entire bandwidth of the excitation pulse. The proposed methodology is called in this work Local Shape Factor Calibration Method (LSFCM).

Aspectos morfológicos comparativos entre neurônios da camada I do córtex visual de duas espécies de roedores: Cavia porcellus e Rattus norvegicus

A camada I tem como característica principal a baixa quantidade de neurônios e uma alta densidade de fibras nervosas. A morfologia dos neurônios da camada I ainda é pouco estudada, tanto que nos estudos que avaliaram a morfologia desses neurônios não se chegou ainda a um consenso sobre as formas e funções desses neurônios. Este estudo avaliou a morfologia dos neurônios da camada I do córtex visual de duas espécies de roedores: Cavia porcellus, popularmente conhecido no Brasil como porquinho-da-índia e Rattus norvegicus, que é o rato e foi utilizada a linhagem Wistar, comumente usado nas pesquisas científicas. O porquinho-da-índia é um modelo animal muito estudado, utilizado em diversos segmentos da ciência. Apesar dessa espécie ser bem estudada, trabalhos na camada I desse animal são relativamente raros, especialmente em relação à morfologia e eletrofisiologia dos neurônios dessa região cortical. Pesquisas em ratos sobre os neurônios da camada I são mais frequentes, tanto em relação a morfologia quanto a eletrofisiologia. Para discriminar as possibilidades de diferenças na morfologia dos neurônios da camada I do córtex visual do porquinho-da-índia e do rato, este estudo classificou esses neurônios de acordo com a trajetória de seus dendritos e analisou as medidas dendríticas utilizando a técnica de injeção intracelular de biocitina. Após a classificação dos neurônios as comparações foram feitas entre os mesmos tipos celulares de cada roedor. Foram utilizados 35 porquinhos-da-índia da variedade Dunkin-Hartley de pêlo curto de ambos os sexos com idades de 4 a 5 dias de vida pós-natal. Quanto aos ratos, foram utilizados 30 ratos da variedade Wistar, de ambos os sexos com idades de 14 a 21 dias de vida pós-natal. Os animais foram anestesiados e tiveram seus encéfalos removidos, hemisférios separados e foram realizados cortes no plano coronal na região occipital onde se localiza a área visual dos roedores. As fatias foram mantidas em líquido cérebro-espinhal artificial e em seguida levadas ao microscópio para injeção de biocitina e posteriormente foram fixadas e tratadas para montagem em lâmina e contracoradas com Nissl para melhor visualização. Os neurônios encontrados foram classificados como: horizontais, ascendentes, descendentes e radias. Foram analisadas as seguintes medidas dendríticas: área do campo receptor, comprimento dendrítico total e médio, área total do corpo celular, número de dendritos, distância da pia-máter e análise da distribuição de Sholl. Dos resultados obtidos os mais notáveis foram o alcance dos ramos dendríticos e o tamanho do corpo celular dos neurônios da camada I do porquinho-da-índia quando comparados aos do rato. Isso sugere que, nessa espécie, um maior número de microcircuitos neurais podem ser estabelecidos, e por conseguinte maior taxa metabólica, justificada pelo maior tamanho do corpo celular.