Equação de estado de Murnaghan
A equação de estado de Murnaghan é a relação entre o volume de um corpo e a pressão a que ele está submetido. Esta é uma das muitas equações de estado que têm sido usadas em ciências da Terra e física do choque para modelar o comportamento da matéria sob condições de alta pressão. Deve seu nome a Francis D. Murnaghan[1] que propôs em 1944 refletir o comportamento do material sob uma faixa de pressão tão ampla quanto possível para refletir um fato experimentalmente estabelecido: quanto mais um sólido é comprimido, mais difícil é comprimi-lo ainda mais.
A equação de Murnaghan é derivada, sob certas suposições, a partir das equações da mecânica de meios contínuos. Envolve dois parâmetros ajustáveis: o módulo de incompressibilidade K0 e sua primeira derivada em relação à pressão, K′0, ambos medidos à pressão ambiente. Em geral, esses coeficientes são determinados por uma regressão em valores obtidos experimentalmente de volume V como uma função da pressão P. Esses dados experimentais podem ser obtidos por difração de raios X ou por testes de choque. A regressão também pode ser realizada nos valores da energia em função do volume obtido a partir de ab-initio e cálculos de dinâmica molecular.
A equação de estado de Murnaghan é tipicamente expressa como:
Se a redução do volume sob compressão for baixa, i.e., para V/V0 maior que aproximadamente 90%, a equação de Murnaghan pode modelar dados experimentais com precisão satisfatória. Além disso, ao contrário de muitas equações de estado propostas, fornece uma expressão explícita do volume em função da pressão V(P). Mas o seu alcance de validade é limitado e a interpretação física inadequada. No entanto, esta equação de estado continua a ser amplamente utilizada em modelos de explosivos sólidos. Das equações de estado mais elaboradas, a mais utilizada na física da Terra é a equação de estado de Birch–Murnaghan. Na física de choque de metais e ligas, outra equação de estado amplamente utilizada é a equação de estado de Mie–Grüneisen.
Plano de fundo
[editar | editar código-fonte]O estudo da estrutura interna da Terra através do conhecimento das propriedades mecânicas dos constituintes das camadas internas do planeta envolve condições extremas; a pressão pode ser contada em centenas de gigapascal e as temperaturas em milhares de graus. O estudo das propriedades da matéria nessas condições pode ser feito experimentalmente através de dispositivos como célula de bigorna de diamante para pressões estáticas, ou submetendo o material a ondas de choque. Deu também origem a trabalhos teóricos para determinar a equação de estado, ou seja, as relações entre os diferentes parâmetros que definem neste caso o estado da matéria: o volume (ou densidade), a temperatura e a pressão.
Existem duas abordagens:
- as equações de estado derivadas de potenciais interatômicos, ou possivelmente cálculos ab initio;
- derivado das relações gerais das equações de estado da mecânica e da termodinâmica. A equação de Murnaghan pertence a esta segunda categoria.
Dezenas de equações foram propostas por vários autores.[2] Estas são relações empíricas, a qualidade e a relevância dependem da utilização que delas se faz e podem ser julgadas por diferentes critérios: o número de parâmetros independentes envolvidos, o significado físico que pode ser atribuído a esses parâmetros, a qualidade dos dados experimentais , e a consistência dos pressupostos teóricos que fundamentam sua capacidade de extrapolar o comportamento dos sólidos em alta compressão.[3]
Expressões para a equação de estado
[editar | editar código-fonte]Geralmente, a temperatura constante, o módulo de volume é definido por: A maneira mais fácil de obter uma equação de estado relacionando P e V é assumir que K é constante, ou seja, independente da pressão e da deformação do sólido, então simplesmente encontramos a lei de Hooke. Neste caso, o volume diminui exponencialmente com a pressão. Este não é um resultado satisfatório porque está experimentalmente estabelecido que à medida que um sólido é comprimido, torna-se mais difícil de comprimir. Para ir mais longe, devemos levar em conta as variações das propriedades elásticas do sólido com a compressão.
A suposição de Murnaghan é assumir que o módulo de volume é uma função linear da pressão:[1] A equação de Murnaghan é o resultado da integração da equação diferencial: Também podemos expressar o volume dependendo da pressão:
Esta apresentação simplificada é, no entanto, criticada por Poirier como carente de rigor.[4] A mesma relação pode ser demonstrada de forma diferente do fato de que a incompressibilidade do produto do módulo e do coeficiente de expansão térmica não depende da pressão para um determinado material.[5] Esta equação de estado também é um caso geral da relação politropo mais antiga [6] que também tem uma relação de força constante.
Em algumas circunstâncias, especialmente em relação a cálculos ab initio, , a expressão da energia em função do volume será preferida,[7] a qual pode ser obtido integrando a equação acima de acordo com a relação P = −dE/dV. Pode ser escrita para K′0 diferente de 3,
Derivação da equação de estado de Murnaghan: | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Um sólido tem um certo volume de equilíbrio , e a energia aumenta quadraticamente à medida que o volume aumenta ou diminui uma pequena quantidade desse valor. A dependência plausível mais simples da energia em relação ao volume seria um sólido harmônico, com
O próximo modelo razoável mais simples seria com um módulo volumétrico constante Integrando-se resulta Uma equação de estado mais sofisticada foi derivada por Francis D. Murnaghan da Universidade Johns Hopkins em 1944[1]. Para começar, consideramos a pressão
Experimentalmente, a derivada da pressão do módulo volumétrico
descobriu-se que muda pouco com a pressão. Se tomarmos como sendo uma constante, então
onde é o valor de quando Podemos equiparar isso com (2) e reorganizar como
Integrando isso resulta em
ou equivalentemente
Substituindo (6) em quando então resulta na equação de estado para energia.
Muitas substâncias têm um valor bastante constante de aproximadamente 3,5. |
Vantagens e limitações
[editar | editar código-fonte]Apesar da sua simplicidade, a equação de Murnaghan é capaz de reproduzir os dados experimentais para uma gama de pressões que pode ser bastante grande, da ordem de K0/2.[8] Também permanece satisfatório, uma vez que o razão V/V0 permanece acima de 90%.[9] Nesta faixa, a equação de Murnaghan tem uma vantagem em relação a outras equações de estado se quisermos expressar o volume em função da pressão.[10]
No entanto, outras equações podem fornecer melhores resultados e vários estudos teóricos e experimentais mostram que a equação de Murnaghan é insatisfatória para muitos problemas. Assim, na medida em que a proporção V/V0 torna-se muito baixa, a teoria prevê que K′ tende a 5/3, que é o limite de Thomas–Fermi.[10][11][nota 1] No entanto, na equação de Murnaghan, K′ é constante e definido para seu valor inicial. Em particular, o valor K′0 = 5/3 torna-se inconsistente com a teoria em algumas situações. Na verdade, quando extrapolado, o comportamento previsto pela equação de Murnaghan torna-se rapidamente improvável.[10]
Independentemente deste argumento teórico, a experiência mostra claramente que K′ diminui com a pressão, ou em outras palavras, que a segunda derivada do módulo de incompressibilidade K″ é estritamente negativo. Uma teoria de segunda ordem baseada no mesmo princípio (ver secção seguinte) pode explicar esta observação, mas esta abordagem ainda é insatisfatória. Na verdade, isso leva a um módulo de volume negativo no limite onde a pressão tende ao infinito. Na verdade, esta é uma contradição inevitável, qualquer que seja a expansão polinomial escolhida, porque sempre haverá um termo dominante que diverge para o infinito.[3]
Estas importantes limitações levaram ao abandono da equação de Murnaghan, o que W. Holzapfel chama "uma forma matemática útil sem qualquer justificativa física".[12] Na prática, a análise dos dados de compressão é feita utilizando equações de estado mais sofisticadas. A mais comumente usado na comunidade científica é a equação de Birch–Murnaghan, segunda ou terceira ordem na qualidade dos dados coletados.[13]
Finalmente, uma limitação muito geral deste tipo de equação de estado é a sua incapacidade de levar em conta as transições de fase induzidas pela pressão e temperatura de fusão, mas também múltiplas transições sólido-sólido que podem causar mudanças abruptas na densidade e no módulo aparente baseadas na pressão.[3]
Exemplos
[editar | editar código-fonte]Na prática, a equação de Murnaghan é utilizada para realizar uma regressão em um conjunto de dados, onde se obtém os valores dos coeficientes K0 e K′0. Obtidos estes coeficientes, e conhecendo o valor do volume para as condições ambientais, somos, em princípio, capazes de calcular o volume, a densidade e o módulo aparente para qualquer pressão.
O conjunto de dados é principalmente uma série de medições de volume para diferentes valores de pressão aplicada, obtidos principalmente por difração de raios X. Também é possível trabalhar com dados teóricos, calculando a energia para diferentes valores de volume por métodos ab initio, e então regredindo esses resultados. Isto dá um valor teórico do módulo de elasticidade que pode ser comparado com resultados experimentais.
A tabela seguinte lista alguns resultados de diferentes materiais, com o único propósito de ilustrar algumas análises numéricas que têm sido feitas através da equação de Murnaghan, sem prejuízo da qualidade dos modelos obtidos. Dadas as críticas feitas na seção anterior sobre o significado físico da equação de Murnaghan, estes resultados devem ser considerados com cautela.
Material | (GPa) | |
---|---|---|
NaF[5] | 46,5 | 5,28 |
NaCl[5] | 24,0 | 5,39 |
NaBr[5] | 19,9 | 5,46 |
NaI[5] | 15,1 | 5,59 |
MgO[8] | 156 | 4,7 |
Calcita (CaCO3)[14] | 75,27 | 4,63 |
Magnesita (MgCO3)[15] | 124,73 | 3,08 |
Carbeto de silício (3C-SiC)[16] | 248 | 4,0 |
Extensões e generalizações
[editar | editar código-fonte]Para melhorar os modelos ou evitar críticas descritas acima, foram propostas várias generalizações da equação de Murnaghan. Geralmente consistem em abandonar uma suposição simplificadora e adicionar outro parâmetro ajustável. Isto pode melhorar as qualidades de refinamento, mas também levar a expressões complicadas. A questão do significado físico destes parâmetros adicionais também é levantada.
Uma estratégia possível é incluir um termo adicional P2 no desenvolvimento prévio,[17][18] exigindo que . Resolver esta equação diferencial dá a equação do Murnaghan de segunda ordem: onde . Encontrado naturalmente na equação de primeira ordem tomando-se . Desenvolvimentos para uma ordem superior a 2 são possíveis em princípio,[19] mas ao custo de adicionar um parâmetro ajustável para cada termo.
Outras generalizações podem ser citadas:
- Kumari e Dass propuseram uma generalização abandonando a condição K = 0 mas assumindo o resultado K / K′ independente da pressão;[20]
- Kumar propôs uma generalização levando em consideração a dependência do parâmetro de Anderson em função do volume. Posteriormente, foi demonstrado que esta equação generalizada não era nova, mas sim redutível à equação de Tait.[5][21]
Notas
- ↑ A teoria de Thomas-Fermi considera um sólido fortemente comprimido como um gás de elétrons degenerado (gás de Fermi) com um termo adicional de triagem para levar em conta a presença de núcleos atômicos.
Referências
- ↑ a b Murnaghan, F.D. (1944). «The Compressibility of Media under Extreme Pressures». Proceedings of the National Academy of Sciences of the United States of America. 30 (9): 244–247. Bibcode:1944PNAS...30..244M. PMC 1078704. PMID 16588651. doi:10.1073/pnas.30.9.244
- ↑ Wedepohl, P.T. (1972). «Comparison of a simple two-parameter equation of state with the Murnaghan equation». Solid State Communications. 10 (10): 947–951. Bibcode:1972SSCom..10..947W. doi:10.1016/0038-1098(72)90228-1
- ↑ a b c Stacey, F.D.; Brennan, B.J.; Irvine, R.D. (1981). «Finite strain theories and comparison with seismological data». Surveys in Geophysics. 4 (3): 189–232. Bibcode:1981GeoSu...4..189S. doi:10.1007/bf01449185[ligação inativa]
- ↑ Poirier (2002), p. 65.
- ↑ a b c d e f Kumar, M. (1995). «High pressure equation of state for solids». Physica B: Condensed Matter. 212 (4): 391–394. Bibcode:1995PhyB..212..391K. doi:10.1016/0921-4526(95)00361-C
- ↑ Weppner, S. P., McKelvey, J. P., Thielen, K. D. and Zielinski, A. K., "A variable polytrope index applied to planet and material models", "Monthly Notices of the Royal Astronomical Society", Vol. 452, No. 2 (Sept. 2015), pages 1375–1393, Oxford University Press also found at the arXiv
- ↑ Silvi (1997), p. 122.
- ↑ a b Anderson, O.L. (1995). Equations of state of solids for geophysics and ceramic science. [S.l.]: Oxford University Press. 179 páginas. ISBN 9780195345278
- ↑ Angel, R.J. «Some practical aspects of studying equations of state and structural phase transitions at high pressure». High-Pressure Crystallography. Col: NATO Science Series. Dordrecht: Springer. pp. 21–36. ISBN 978-1-4020-1954-8
- ↑ a b c Holzapfel, W.B. (1996). «Physics of solids under strong compression». Reports on Progress in Physics. 59 (1): 29–90. Bibcode:1996RPPh...59...29H. doi:10.1088/0034-4885/59/1/002
- ↑ Tal, Yoram; Levy, Mel (Fevereiro de 1981). «Expectation values of atoms and ions: The Thomas-Fermi limit». Phys. Rev. A. 23 (2): 408--415. doi:10.1103/PhysRevA.23.408. Consultado em American Physical Society Verifique data em:
|acessodata=
(ajuda) - ↑ Holzapfel, W.B. (2001). «Equations of state for solids under strong compression». Zeitschrift für Kristallographie. 216 (9): 473–488. Bibcode:2001ZK....216..473H. doi:10.1524/zkri.216.9.473.20346
- ↑ Boldyreva, E.; Dera, P.; Ballaran, T. Boffa. «Equations of state and their applications in geosciences». High-Pressure Crystallography: From Fundamental Phenomena to Technological Applications. [S.l.]: Springer. pp. 135–145
- ↑ Silvi,1997. p. 123.
- ↑ Silvi, 1997.
- ↑ Strössner, K.; Cardona, M.; Choyke, W. J. (1987). «High pressure X-ray investigations on 3C-SiC». Solid State Communications. 63 (2): 113–114. Bibcode:1987SSCom..63..113S. doi:10.1016/0038-1098(87)91176-8
- ↑ MacDonald, J.R.; Powell, D.R. (1971). «Discrimination Between Equations of State». Journal of Research of the National Bureau of Standards Section A. 75 (5). 441 páginas. doi:10.6028/jres.075A.035
- ↑ MacDonald, 1969, p. 320
- ↑ Fuchizaki, Kazuhiro (2006). «Murnaghan equation of state revisited». Journal of the Physical Society of Japan. 75 (3). 034601 páginas. Bibcode:2006JPSJ...75c4601F. doi:10.1143/jpsj.75.034601
- ↑ Kumari, M.; Dass, N. (1990). «An equation of state applied to sodium chloride and caesium chloride at high pressures and high temperatures». Journal of Physics: Condensed Matter. 2 (14): 3219–3229. Bibcode:1990JPCM....2.3219K. doi:10.1088/0953-8984/2/14/006
- ↑ Shanker, J.; Singh, B.; Kushwah, S.S. (1997). «On the high-pressure equation of state for solids». Physica B: Condensed Matter. 229 (3–4): 419–420. Bibcode:1997PhyB..229..419S. doi:10.1016/S0921-4526(96)00528-5
Bibliografia
[editar | editar código-fonte]- Poirier, J.P. (2002). Introduction to the physics of the Earth's interior. [S.l.]: Cambridge University Press. ISBN 9780521663922
- Silvi, B.; P., d'Arco (1997). Modelling of Minerals and Silicated Materials. [S.l.]: Kluwer Academic Publishers. ISBN 9780792343332
- MacDonald, J.R. (1969). «Review of Some Experimental and Analytical Equations of State». Reviews of Modern Physics. 41 (2): 316–349. Bibcode:1969RvMP...41..316M. doi:10.1103/revmodphys.41.316