Início/Lições/Lição 4

O primeiro princípio da termodinâmica

Trabalho, calor, energia interna e o primeiro princípio.

Energia internaPrimeiro princípioTrabalho de pressãoTransferência de calorFunções de estadoDiferencial exataTransformações quase-estáticasTransformação adiabáticaExtensividade

A lição anterior introduziu os sistemas termodinâmicos, os estados de equilíbrio e as funções de estado. Vamos agora responder a uma questão central: como incorporar a energia do sistema nesta representação macroscópica?

Em certos cursos esta questão é tratada depressa demais: o primeiro princípio é apresentado como uma simples conservação da energia e, logo a seguir, surge a fórmula do balanço energético sem justificações adicionais. Veremos essas fórmulas nas secções 3 e 4.3, e o leitor apressado pode com efeito dirigir-se diretamente a elas (ver também a síntese rápida na Secção 7).

Mas a construção que conduz a essas fórmulas é, na realidade, subtil. Em primeiro lugar, o primeiro princípio afirma mais do que a conservação da energia, que (nos dias de hoje) parece evidente do ponto de vista microscópico. Vamos mostrar que este princípio preenche na realidade o vazio lógico entre a existência de uma energia UmicroU_{\mathrm{micro}} definida mecanicamente num espaço com cerca de 6N6N coordenadas e a de uma função de estado energia UthU_{\mathrm{th}} que depende apenas de algumas variáveis macroscópicas.

Mostraremos em seguida que a aditividade da energia é uma hipótese suplementar, necessária para passar da conservação da energia de um sistema isolado a um balanço de energia na fronteira de um sistema fechado. É esta última etapa que permite precisar enfim o estatuto epistémico e lógico do calor: este é definido como toda a transferência de energia que resta nesse balanço, uma vez contabilizado o trabalho macroscópico realizado sobre um sistema fechado: Q  =def  ΔUWQ \equiv \Delta U - W.

Nesta lição consideraremos primeiro os sistemas fechados: nenhuma matéria atravessa a sua fronteira. A razão é simples: num sistema aberto, a matéria transporta ela própria energia através da fronteira e complica o balanço energético. O caso dos sistemas abertos será abordado mais adiante no curso.

1. A energia interna microscópica

Consideremos um sistema clássico, não quântico e não relativista num referencial de inércia1. Consideremos NN partículas pontuais, de massas mim_i, posições ri\vec r_i e velocidades vi\vec v_i. Em cada instante, o seu estado microscópico fica completamente determinado por

Note 1 : Esses casos exigem um tratamento adaptado, fora do âmbito deste curso.
(r1,...,rN,p1,...,pN)Γ,compi=mivi,\boxed{ (\vec r_1,...,\vec r_N, \vec p_1,...,\vec p_N) \in\Gamma, \, \mathrm{com} \,\,\, \vec p_i=m_i\vec v_i, }

onde Γ\Gamma designa o espaço de fases clássico do sistema, de dimensão dimΓ=6N\dim\Gamma=6N neste modelo simples. Supomos que todas as forças internas, por exemplo as forças de atração ou de repulsão entre moléculas, são conservativas. Podemos então associar-lhes uma energia potencial total Epint(r1,...,rN),comFiint=riEpintE_p^{\mathrm{int}}(\vec r_1,...,\vec r_N), \, \mathrm{com} \,\,\, \vec F_i^{\,\mathrm{int}} =-\nabla_{\vec r_i}E_p^{\mathrm{int}}. O sistema pode ainda estar sujeito a forças conservativas exteriores, por exemplo o seu peso. Associa-se-lhes do mesmo modo uma energia potencial exterior Epext(r1,...,rN)E_p^{\mathrm{ext}}(\vec r_1,...,\vec r_N).

O sistema não está necessariamente em repouso no referencial escolhido. Em cada instante podem definir-se o seu centro de massa e a velocidade V\vec V deste, eventualmente dependente do tempo. Introduzindo então a massa total do sistema e as velocidades relativas ao centro de massa,

M=i=1Nmi,V=1Mi=1Nmivi,ui=viV,M=\sum_{i=1}^{N} m_i, \qquad \vec V=\frac{1}{M}\sum_{i=1}^{N} m_i\vec v_i, \qquad \vec u_i=\vec v_i-\vec V,

a energia cinética total pode escrever-se como

Ec=12i=1Nmivi2=12i=1Nmi(V+ui)2.E_c=\frac12 \sum_{i=1}^{N}m_i \vec {v_i}^2 = \frac12 \sum_{i=1}^{N}m_i (\vec V + \vec u_i)^2.

Ao desenvolver o quadrado, verifica-se que os termos cruzados se anulam exatamente, uma vez que imiui=0\sum_i m_i\vec u_i=\vec 0. Este resultado é conhecido em mecânica clássica como teorema de König:

Ec=12MV2+12i=1Nmiui2+Vi=1NmiuiimiviMV=0E_c =\frac12MV^2 +\frac12 \sum_{i=1}^{N}m_i \vec{u_i}^2 +\vec V\mathbin{\cdot} \underbrace{\sum_{i=1}^{N} m_i\vec u_i}_{\sum_i m_i\vec v_i-M\vec V=\vec0}

Vem então:

Ec=Ecmacro+Ecmicro\boxed{E_c = E_c^{\mathrm{macro}} + E_c^{\mathrm{micro}}}

(1)

A energia cinética total decompõe-se, pois, sempre2 na energia cinética macroscópica do centro de massa e na dos movimentos relativos residuais, os «desordenados», que constituem aquilo a que se chama agitação térmica.

Note 2 : A translação do centro de massa não esgota necessariamente a energia cinética macroscópica. O sistema pode possuir ainda um movimento ordenado de rotação global. Nesse caso pode escrever-se a decomposição vi=V+ω×(riR)+ui,\vec v_i = \vec V +\vec\omega\times(\vec r_i-\vec R) +\vec u_i', onde R\vec R é a posição do centro de massa e onde ui\vec u_i' designa agora o movimento residual após subtração da translação e da rotação globais. Escolhendo ω\vec\omega de tal modo que esse movimento residual tenha momento angular total nulo relativamente ao centro de massa, obtém-se Ec=12MV2+12ωICMω+12i=1Nmiui2,E_c = \frac12 MV^2 + \frac12\vec\omega\cdot \mathbf I_{\mathrm{CM}}\vec\omega + \frac12\sum_{i=1}^{N}m_i {u_i'}^2, onde ICM\mathbf I_{\mathrm{CM}} é o tensor de inércia relativamente ao centro de massa. A separação entre energia cinética macroscópica e energia cinética microscópica permanece, portanto, possível.

Reagrupando os diferentes termos, a energia mecânica total escreve-se finalmente

Etot=Ecmacro+Epext+[Ecmicro+Epint(r1,...,rN)+Eintautres]energia interna.E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +\underbrace{\left[ E_c^{\mathrm{micro}} +E_p^{\mathrm{int}}(\vec r_1,...,\vec r_N) +E_{\mathrm{int}}^{\mathrm{autres}} \right]}_{\text{energia interna}}.

Observação: o termo EintautresE_{\mathrm{int}}^{\mathrm{autres}} reúne as contribuições dos graus de liberdade internos que este modelo pontual mínimo não descreve, por exemplo os graus de liberdade de rotação ou de vibração das moléculas. Se estes forem modelados explicitamente, o espaço de fases Γ\Gamma deve ser alargado em conformidade.

Definição 1 (Energia interna microscópica)
Designaremos por Umicro=Ecmicro+Epint+EintautresU_{\mathrm{micro}}= E_c^{\mathrm{micro}} +E_p^{\mathrm{int}} +E_{\mathrm{int}}^{\mathrm{autres}} a energia interna microscópica. É uma função definida no espaço de fases, Umicro:ΓR,U_{\mathrm{micro}}:\Gamma\longrightarrow\mathbb R,

que depende, portanto, em geral, de um número considerável de variáveis, e que satisfaz:

Etot=Ecmacro+Epext+Umicro.E_{\mathrm{tot}} =E_c^{\mathrm{macro}}+E_p^{\mathrm{ext}}+U_{\mathrm{micro}}.

(2)

No que se segue consideraremos sistemas macroscopicamente em repouso, desprovidos de rotação global e cuja energia potencial exterior não varia. A menos de uma escolha de origem, poderemos então escrever

Etot=Umicro.E_{\mathrm{tot}}=U_{\mathrm{micro}}.

Vamos agora introduzir a função de estado energia termodinâmica, provisoriamente denotada UthU_{\mathrm{th}}, e compreender por que razão não deve ser confundida a priori com UmicroU_{\mathrm{micro}}.

2. A energia interna termodinâmica

Tomemos como sistema uma panela de água colocada sobre uma placa de aquecimento. Mantemo-nos no quadro simples descrito acima: a panela está macroscopicamente imóvel e a sua energia potencial exterior não varia. Podemos portanto escrever Etot=UmicroE_{\mathrm{tot}}=U_{\mathrm{micro}} a menos de uma constante.

Quando se aquece a água, não aparece qualquer energia cinética macroscópica e a panela não muda de posição. Admitindo a conservação estrita da energia, sob todas as suas formas, a energia que lhe é fornecida encontra-se necessariamente na sua energia interna microscópica:

ΔEtot=ΔUmicro.\Delta E_{\mathrm{tot}}=\Delta U_{\mathrm{micro}}.

A termodinâmica procura descrever exatamente a mesma variação. Gostaríamos, pois, de construir uma grandeza UthU_{\mathrm{th}} tal que

ΔUth=ΔUmicro.\Delta U_{\mathrm{th}}=\Delta U_{\mathrm{micro}}.

A solução parece evidente. Não bastará identificar a energia interna termodinâmica com a energia interna microscópica?

Essa identificação seria, na realidade, matematicamente incorreta, uma vez que as duas funções não têm o mesmo domínio de definição. UmicroU_{\mathrm{micro}} é uma função definida no espaço de fases e depende de pelo menos 6N6N variáveis, ao passo que UthU_{\mathrm{th}} depende apenas de um pequeno número dd de variáveis independentes que descrevem o estado de equilíbrio X=(x1,...,xd)EX=(x^1,...,x^d)\in\mathcal E (recorde-se que os xix^i são variáveis de estado macroscópicas, como a temperatura, a pressão ou o volume).

As duas funções devem, contudo, estar relacionadas de modo que Uth(X)U_{\mathrm{th}}(X) represente, à escala macroscópica, a energia interna correspondente ao estado de equilíbrio XX. É preciso explicar como funciona esta redução dimensional. A questão é profunda e ultrapassa o âmbito deste curso de introdução. É necessária uma descrição microscópica se quisermos deduzir esta passagem entre as duas escalas. A subsecção 6.2, fora do programa, apresenta as grandes linhas do modo como a física estatística resolve este problema.

Mas de um ponto de vista estritamente termodinâmico não há justificação possível desta identificação microscópica, e a existência de UthU_{\mathrm{th}} como função de apenas algumas variáveis macroscópicas tem de ser postulada. Será esse o papel do primeiro princípio da termodinâmica, que portanto não afirma apenas a conservação da energia, mas também, e sobretudo, que essa energia pode ser representada, no equilíbrio, por uma função de estado definida no espaço macroscópico E\mathcal E.

Uma vez escolhidas as dd variáveis de estado independentes xix^i, por exemplo (TT, VV, NN) para um gás ideal, a função U(T,V,N)U(T,V,N) define um campo escalar no espaço dos estados; no espaço aumentado das coordenadas U,T,V,NU,T,V,N, o seu gráfico é uma hipersuperfície, dita superfície de equilíbrio, cujo papel no resto deste curso será crucial, cf. Figura 1.

A hipersuperfície energia definida sobre o espaço dos estados de equilíbrio.
Figura 1. A hipersuperfície energia definida sobre o espaço dos estados de equilíbrio.

Notemos por fim um corolário imediato que ilustra bem a diferença entre UthU_{\mathrm{th}} e UmicroU_{\mathrm{micro}}: na termodinâmica de equilíbrio aqui desenvolvida, a energia interna termodinâmica não está, em geral, definida fora dos estados de equilíbrio. Durante uma transformação brusca de um gás, por exemplo, é possível, e mesmo frequente, que já não se lhe possa atribuir uma temperatura ou uma pressão únicas. Nesse caso também não se lhe pode associar diretamente a energia interna termodinâmica UthU_{\mathrm{th}}, uma vez que esta só está definida no espaço dos estados de equilíbrio E\mathcal E. O problema não é que a energia deixe de existir, pois a sua energia microscópica permanece perfeitamente definida em cada instante, mas simplesmente que o sistema já não é representado por um ponto XEX\in\mathcal E.

Dito isto, a partir da secção seguinte a função mecânica UmicroU_{\mathrm{micro}} deixará de intervir diretamente: escreveremos simplesmente UU para a energia interna termodinâmica UthU_{\mathrm{th}}, subentendendo-se que está definida apenas no espaço dos estados de equilíbrio.

Observação 1 (Extensões próximas do equilíbrio)
Assinalemos que, para sistemas suficientemente próximos do equilíbrio, a descrição termodinâmica pode ser alargada aumentando o espaço dos estados. Por exemplo, no caso já encontrado da barra metálica situada entre dois reservatórios térmicos, pode introduzir-se um campo local estacionário de temperatura T(x)T(x), em vez de uma temperatura única, e definir uma densidade local de energia interna. A energia total torna-se então um funcional desses campos. Não desenvolveremos aqui esta termodinâmica dita próxima do equilíbrio.

3. O primeiro princípio da termodinâmica

A secção anterior preparou amplamente o terreno; podemos agora enunciar:

Primeiro princípio da termodinâmica para um sistema fechado
Para todo o sistema termodinâmico fechado existe uma função de estado U:ER,U:\mathcal E\longrightarrow\mathbb R,

dita energia interna, determinada apenas a menos de uma constante aditiva, pelo menos duas vezes diferenciável, tal que, para todo o estado de equilíbrio XEX\in\mathcal E,

Etot=Ecmacro+Epext+U(X).E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +U(X).
(3)

O primeiro princípio afirma igualmente a conservação da energia: para um sistema isolado,

ΔEtot=0.\Delta E_{\mathrm{tot}}=0.

A regularidade da função UU é também uma hipótese não trivial. Na prática poderá considerar-se UU como C\mathcal C^\infty. Como UU só está determinada a menos de uma constante, por causa das energias potenciais envolvidas na igualdade 3, apenas as suas variações têm significado físico. Para simplificar, limitar-nos-emos no que se segue às transformações em que as energias cinética e potencial macroscópicas não variam. Teremos então ΔEtot=ΔU\Delta E_{\mathrm{tot}}=\Delta U.

Consideremos agora um sistema fechado SS, que pode trocar energia, mas não matéria, com o seu ambiente. O sistema global

S=S+ExtS'=S+\mathrm{Ext}

pode sempre ser incluído num sistema isolado. O primeiro princípio impõe então

ΔES=0.\Delta E_{S'}=0.

Sob a hipótese suplementar de que a energia é aditiva, isto é, de que a sua energia de interação é desprezável, tem-se ES=ES+EExtE_{S'}=E_S+E_{\mathrm{Ext}}, de modo que

ΔES=ΔEExt.\Delta E_S=-\Delta E_{\mathrm{Ext}}.

A hipótese de aditividade permite assim traduzir a conservação da energia sob a forma de um balanço entre subsistemas: toda a energia ganha pelo sistema é perdida pelo seu ambiente, e reciprocamente. Resta agora precisar como pode essa energia ser transferida através da fronteira do sistema. Distinguiremos dois modos de transferência: o trabalho, denotado WW, que possui uma definição independente herdada da mecânica, e o calor, denotado QQ.

O balanço energético tomará então a forma

ΔES=ΔU=Q+W\boxed{ \Delta E_S = \Delta U = Q+W }
(4)

onde WW designa o trabalho recebido pelo sistema e QQ o calor por ele recebido.

Note-se que a hipótese de aditividade da energia não é trivial. A um nível elementar poderá simplesmente admitir-se para os sistemas termodinâmicos habituais. Não é, contudo, sempre válida. As interações de longo alcance, em particular, podem pôr em causa simultaneamente a aditividade e a extensividade da energia, ver a Secção 6.4.

Resta agora definir estes dois modos de transferência. Antes disso, porém, é imprescindível registar uma convenção de sinais, bastante padrão, que se utiliza em todo o resto deste curso.

Para lembrar (A convenção de sinais do banqueiro)
Em todo este livro, uma transferência de energia é contada positivamente quando é recebida pelo sistema e negativamente quando é fornecida por ele. Assim, Q>0, W>0Q>0,\ W>0

correspondem a energia recebida pelo sistema, ao passo que

Q<0, W<0Q<0,\ W<0

correspondem a energia cedida ao exterior.

4. Trabalho, calor e balanço energético

4.1. O trabalho

Uma das vantagens do trabalho é possuir uma definição independente da termodinâmica, herdada da mecânica. Quando uma força exterior Fext\vec F_{\mathrm{ext}} atua sobre um ponto da fronteira do sistema e esse ponto se desloca de drd\vec r, o trabalho elementar fornecido por essa força, e portanto recebido pelo sistema, vale

δW=Fextdr.\delta W = \vec F_{\mathrm{ext}}\cdot d\vec r.

Para um sistema extenso, há que somar as contribuições de todas as forças exteriores que realizam trabalho sobre a sua fronteira.

O caso que mais nos servirá na prática é o de um fluido contido num recinto com uma parede móvel, o que conduz a avaliar o trabalho das forças de pressão. Consideremos um gás encerrado num cilindro por um pistão plano de secção S\mathcal S, cf. Figura 2. O eixo xx está orientado para o exterior do gás. Quando o pistão se desloca de uma quantidade dxdx, o volume varia de

dV=Sdx.dV=\mathcal S\,dx.

Denotemos por PextP_{\mathrm{ext}} a pressão exterior exercida sobre a fronteira do gás. A força recebida pelo gás está dirigida para o interior e vale

Fext=PextSex.\vec F_{\mathrm{ext}} = -P_{\mathrm{ext}}\mathcal S\,\vec e_x.

O trabalho elementar recebido pelo gás é, pois,

δWpression=Fext(dxex)=PextSdx,\delta W_{\mathrm{pression}} = \vec F_{\mathrm{ext}}\cdot(dx\,\vec e_x) = -P_{\mathrm{ext}}\mathcal S\,dx,

ou seja

δWpression=PextdV\boxed{ \delta W_{\mathrm{pression}}=-P_{\mathrm{ext}} dV }
(5)

A regra do banqueiro permite verificar imediatamente o sinal. Numa compressão, dV<0dV<0, logo δW>0\delta W>0: o gás recebe trabalho. Numa expansão, dV>0dV>0, logo δW<0\delta W<0: o gás fornece trabalho ao seu ambiente.

Note-se que nesta fórmula se trata efetivamente da pressão exterior ao sistema, que em geral não tem razão para ser igual à pressão do próprio gás. Isto é particularmente importante de ter presente no caso das transformações bruscas, fora do equilíbrio, em que a pressão interna nem sequer está, em geral, definida.

Trabalho recebido por um gás aquando do deslocamento de um pistão de secção S. Numa expansão, dV>0 enquanto a força exercida pela pressão exterior se opõe ao deslocamento: o trabalho recebido é, portanto, negativo.
Figura 2. Trabalho recebido por um gás aquando do deslocamento de um pistão de secção S\mathcal S. Numa expansão, dV>0dV>0 enquanto a força exercida pela pressão exterior se opõe ao deslocamento: o trabalho recebido é, portanto, negativo.

Integrando, obtém-se a fórmula seguinte para uma transformação finita ABA \to B ao longo da qual a pressão exterior é conhecida:

Wpression[AB]=ABPextdV\boxed{ W_{\mathrm{pression}}[A \to B] = -\int_{A}^B P_{\mathrm{ext}} dV }
(6)

Este integral mostra que, para calcular o trabalho total, é preciso conhecer o valor de PextP_{\mathrm{ext}} ao longo de toda a transformação. Isto tem uma consequência maior: o trabalho não pode em geral ser determinado apenas a partir dos estados inicial e final AA e BB: a sua expressão depende do processo seguido.

Três casos particulares merecem ser assinalados aqui.

  1. Se PextP_{\mathrm{ext}} for constante, então Wpression[AB]=PextΔVW_{\mathrm{pression}}[A\to B] = -P_{\mathrm{ext}}\Delta V.
  2. Ao longo de uma transformação isocórica, dV=0dV=0, e o trabalho das forças de pressão é nulo: δWpression=0\delta W_{\mathrm{pression}}=0 e Wpression[AB]=0W_{\mathrm{pression}}[A\to B]=0.
  3. Numa expansão no vazio, Pext=0P_{\mathrm{ext}}=0, pelo que o trabalho das forças de pressão é igualmente nulo.

Observação 2 (Outras formas de trabalho)
O trabalho das forças de pressão é apenas um exemplo entre outros. Para um fio que se estica, se LL designar o seu comprimento e FextF_{\mathrm{ext}} a força de tração recebida, δW=FextdL.\delta W=F_{\mathrm{ext}}\,dL.

Para uma interface cuja área AA se aumenta, a tensão superficial γ\gamma conduz, nas condições habituais, a um trabalho da forma

δW=γdA.\delta W=\gamma\,dA.

Também se pode transferir energia sem deslocamento mecânico visível. Por exemplo, o deslocamento de uma carga dqdq através de uma diferença de potencial pode conduzir a um trabalho elétrico da forma

δWeˊlec=Vextdq,\delta W_{\mathrm{élec}}=V_{\mathrm{ext}}\,dq,

onde VextV_{\mathrm{ext}} é o potencial elétrico exterior, muitas vezes denotado, neste contexto, Φext\Phi_{\mathrm{ext}} para não o confundir com o volume. Também aqui o sinal de dqdq é escolhido em conformidade com a convenção do banqueiro.

A subsecção 6.3 precisará a forma geral dos trabalhos exteriores.

4.2. O calor

O trabalho não pode ser o único modo de transferência de energia, pois aquecer um gás contido num recipiente rígido aumenta a sua temperatura e a sua energia sem que qualquer trabalho exterior seja realizado. Houve, portanto, energia a atravessar a fronteira do sistema sem ter sido transferida sob a forma de trabalho. Este segundo modo de transferência de energia chama-se calor, ou transferência de calor.

Na construção aqui adotada (ver a subsecção 6.1 para uma construção alternativa equivalente), para uma transformação física ABA \to B que liga dois estados de equilíbrio AA e BB, uma vez identificados todos os trabalhos recebidos pelo sistema, o calor recebido é definido por

Q[AB]  =def  U(B)U(A)W[AB].\boxed{ Q[A \to B] \equiv U(B)-U(A)-W[A \to B]. }
(7)

de modo que se reencontra a relação anunciada acima,

ΔU=Q+W\boxed{ \Delta U=Q+W }

Uma vez que o trabalho depende do caminho da transformação ABA \to B, e que a variação de energia ΔU=U(B)U(A)\Delta U = U(B)-U(A) dele não depende, necessariamente a quantidade de calor trocada dele depende também.

Na linguagem termodinâmica, diz-se que a energia interna UU é uma função de estado (por postulado), ao passo que o trabalho WW e o calor QQ não o são. Recorde-se o que isto significa: o facto de WW e QQ não poderem ser funções do estado do sistema num dado instante significa que o sistema não possui «uma quantidade de trabalho ou de calor», mas antes que WW e QQ são transferências de energia na sua fronteira durante uma transformação qualquer ABA \to B, exatamente o que se pretendia obter, à luz das diversas experiências descritas na lição 2.

Notemos que numa transformação cíclica AAA \to A o sistema regressa ao seu estado inicial, pelo que

ΔUcycle=0.\Delta U_{\mathrm{cycle}}=0.

O primeiro princípio impõe então

Qcycle+Wcycle=0.\boxed{Q_{\mathrm{cycle}}+W_{\mathrm{cycle}}=0.}

Uma máquina cíclica não pode, portanto, fornecer trabalho indefinidamente sem receber uma quantidade igual de energia do seu ambiente. É isto, historicamente, o que se chama a impossibilidade do moto perpétuo de primeira espécie.

4.3. Forma diferencial do primeiro princípio

A relação

ΔU=Q+W\Delta U=Q+W

liga dois estados de equilíbrio AA e BB. É sempre verdadeira, no sentido em que não supõe que todos os estados intermédios possam eles próprios ser representados por pontos do espaço dos estados de equilíbrio E\mathcal E. Se a transformação for brusca, o sistema pode abandonar temporariamente a superfície de equilíbrio definida sobre E\mathcal E, ainda que U(A)U(A) e U(B)U(B) permaneçam perfeitamente definidas.

Um caso especial deve agora ser detalhado. Suponhamos que a transformação é quase-estática. Em cada instante o sistema está então, por definição, num estado de equilíbrio. Tal transformação pode ser representada por um caminho

γE\gamma\subset\mathcal E

constituído por estados de equilíbrio infinitamente vizinhos, que não abandona a superfície de equilíbrio. Uma vez que UU é uma função de estado diferenciável em E\mathcal E, a sua variação entre dois estados infinitamente vizinhos é uma diferencial exata, denotada dUdU.

Denotamos então por δQ\delta Q e δW\delta W as transferências elementares de calor e de trabalho correspondentes. Neste caso o primeiro princípio toma a forma

dU=δQ+δW\boxed{ dU=\delta Q+\delta W }
(8)

A diferença de notação entre dUd U, por um lado, e δW\delta W ou δQ\delta Q, por outro, não é meramente cosmética. Indica que a primeira é uma diferencial exata e que as outras duas o não são, o que constitui a tradução matemática do facto de a variação de energia não depender do caminho seguido, ao passo que QQ e WW dele dependem. A secção seguinte recorda ao leitor principiante a matemática necessária para compreender bem este ponto crucial.

No caso particular em que o único trabalho é o das forças de pressão e em que se pode identificar PextP_{\mathrm{ext}} com a pressão PP do sistema, esta relação torna-se

dU=δQPdV.dU=\delta Q-P\,dV.
Para lembrar (Duas formulações do primeiro princípio)
A sua forma integrada ΔU=Q+W\Delta U=Q+W

liga dois estados de equilíbrio e continua utilizável mesmo quando o processo intermédio não é quase-estático. A sua forma diferencial

dU=δQ+δWdU=\delta Q+\delta W

só é válida para uma transformação quase-estática, que possa ser representada por um caminho no espaço dos estados de equilíbrio.

5. Formas diferenciais exatas e inexatas

Esta secção reúne o mínimo de matemática necessário para dar um sentido preciso à distinção entre dUdU, por um lado, e δQ\delta Q e δW\delta W, por outro. O leitor que já a conheça pode saltá-la.

5.1. A diferencial de uma função

Seja ff uma função das variáveis x1,...,xdx^1,...,x^d, suposta derivável. A sua diferencial é a expressão

df=i=1dfxidxi,df=\sum_{i=1}^{d}\frac{\partial f}{\partial x^i}\,dx^i,

que mede a variação de ff quando se passa do ponto XX ao ponto vizinho X+dXX+dX. O ponto importante é o seguinte: se seguirmos um caminho γ\gamma que vai de um ponto AA a um ponto BB, obtemos, somando todas essas variações elementares,

γdf=f(B)f(A).\int_\gamma df=f(B)-f(A).

Em particular, ao longo de um caminho fechado tem-se:

df=0.\oint df=0.

5.2. As formas diferenciais

Consideremos agora uma expressão do mesmo tipo, que não deve de modo algum ser confundida com a diferencial de uma função. Sejam Ai(x1,...,xd)A_i(x^1,...,x^d) funções quaisquer. Definimos então ω\omega por:

ω=i=1dAidxi.\omega=\sum_{i=1}^{d}A_i\,dx^i.

Um objeto deste tipo chama-se uma forma diferencial. Dadas as expressões das funções AiA_i, sabe-se integrá-la sem dificuldade ao longo de um caminho γ\gamma. Mas nada garante que exista uma função ff da qual ω\omega seja a diferencial, isto é, tal que Ai=f/xiA_i=\partial f/\partial x^i para todo o ii.

  • Se tal função ff existir, a forma diz-se exata, e escreve-se ω=df\omega = d f. O integral γω=γdf=f(B)f(A)\int_\gamma\omega = \int_\gamma df = f(B)-f(A) depende então apenas das extremidades.
  • Caso contrário, a forma diz-se inexata, e escreve-se então δω\delta\omega em vez de ω\omega, para lembrar que não é a diferencial de nenhuma função. (Os matemáticos escreveriam simplesmente a forma como ω\omega.)

5.3. O critério de Schwarz

Um critério simples para determinar se uma forma diferencial é exata é o critério de Schwarz. Coloquemo-nos, para simplificar, em duas variáveis (a generalização é imediata), com

ω=A(x,y)dx+B(x,y)dy.\omega=A(x,y)\,dx+B(x,y)\,dy.

Se ω\omega fosse exata, teríamos A=f/xA=\partial f/\partial x e B=f/yB=\partial f/\partial y, logo

Ay=2fyx=2fxy=Bx,\frac{\partial A}{\partial y} =\frac{\partial^2 f}{\partial y\,\partial x} =\frac{\partial^2 f}{\partial x\,\partial y} =\frac{\partial B}{\partial x},

uma vez que a ordem das derivações é indiferente para uma função duas vezes continuamente derivável. Dispomos assim de um teste cómodo:

AyBx    ω na˜eˊ exata.\boxed{ \frac{\partial A}{\partial y}\neq\frac{\partial B}{\partial x} \;\Longrightarrow\; \omega\ \text{não é exata.} }

A menos de subtilezas que não ocorrerão em termodinâmica, o recíproco é verdadeiro.

5.4. Exemplo do trabalho das forças de pressão

Quando as transferências elementares podem ser expressas em função das variáveis de estado ao longo do caminho quase-estático, δQ\delta Q e δW\delta W tornam-se formas diferenciais sobre E\mathcal E, e pode verificar-se se são exatas ou não.

A título de exemplo, tomemos nn moles de gás ideal que sofrem uma transformação quase-estática, com Pext=PP_{\mathrm{ext}} = P. Temos então δW=PdV\delta W=-P\,dV. Mas como P=nRT/VP=nRT/V, o trabalho elementar escreve-se

δW=PdV=nRTVdV=0×dTnRTVdV,ou sejaA=0,B=nRTV,\delta W=-P dV = -\frac{nRT}{V}\,dV = 0\times dT-\frac{nRT}{V}\,dV, \qquad\text{ou seja}\qquad A=0, \quad B=-\frac{nRT}{V},

identificando x=Tx=T e y=Vy=V. Trata-se efetivamente de uma forma diferencial definida no espaço dos estados. O critério de Schwarz dá

AV=0,BT=nRV0.\frac{\partial A}{\partial V}=0, \qquad \frac{\partial B}{\partial T}=-\frac{nR}{V}\neq0.

As duas derivadas cruzadas não são iguais: o trabalho elementar não é, efetivamente, uma diferencial exata. Não existe, portanto, neste caso, nenhuma função de estado W(T,V)W(T,V) da qual δW\delta W fosse a variação, e o trabalho recebido entre dois estados depende de facto do caminho seguido. O mesmo raciocínio vale para δQ=dUδW\delta Q=dU-\delta W: como dUdU é exata e δW\delta W não o é, a sua diferença também não pode ser exata.

Para lembrar (porquê δ\delta e não dd)
dUdU é uma diferencial exata: a sua variação depende apenas dos estados inicial e final, e o seu integral sobre um ciclo é nulo. δQ\delta Q e δW\delta W são formas inexatas: dependem do caminho seguido, e não existe nem uma função QQ nem uma função WW do estado do sistema tais que δW=dW\delta W = d W e δQ=dQ\delta Q = dQ.

6. Para saber mais

Esta secção reúne as considerações mais avançadas anunciadas anteriormente. Pode ser omitida numa primeira leitura.

6.1. Duas construções possíveis do primeiro princípio

Escolhemos aqui a seguinte ordem lógica: a existência da função de estado UU é postulada pelo primeiro princípio, o trabalho é definido independentemente pela mecânica, e o calor é em seguida definido pelo balanço

Q=ΔUW.Q=\Delta U-W.

Daí se retira um facto notável. Se a transformação for adiabática, então ΔU=W\Delta U = W, o que obriga a que o trabalho adiabático WadW_{\mathrm{ad}} seja, por sua vez, sempre independente do caminho seguido.

Isto aponta para outra construção possível da energia interna: operacional, portanto independente de qualquer consideração microscópica, e, neste caso, mais próxima do percurso histórico. Começa-se por caracterizar as transformações adiabáticas sem introduzir previamente a grandeza QQ, a fim de evitar qualquer circularidade de raciocínio.

A ideia é engenhosa: não se define a transferência em si, mas o dispositivo. Uma parede diz-se adiabática quando o estado do sistema que ela encerra só pode ser modificado deslocando as coordenadas mecânicas exteriores, isto é, o pistão, um agitador ou uma corrente elétrica. O teste é então direto: mantêm-se essas coordenadas fixas e modifica-se arbitrariamente o exterior, mergulhando por exemplo o recinto num banho de gelo ou aproximando-o de uma chama. Se nenhuma variável de estado do sistema variar, a parede é adiabática.

Esta caracterização faz intervir apenas estados e deslocamentos mecânicos, nunca uma transferência de energia: precede, portanto, qualquer noção de calor. Na prática aproxima-se dela por um bom isolamento, ou operando depressa face ao tempo de relaxação térmica.

O trabalho recebido, por seu lado, mede-se de forma puramente mecânica ou elétrica. É este todo o conteúdo empírico das experiências de Joule, realizadas entre 1843 e 1850 e retomadas desde então com precisão crescente: uma massa mm que cai de uma altura hh fornece mghmgh, sem que intervenha qualquer transferência de calor. Todos os dispositivos que ele imaginou e realizou num recinto adiabático no sentido acima (moinho de pás, resistência de aquecimento ou compressão de um gás) conduziram-no à mesma constatação: um mesmo trabalho fornecido produz uma mesma mudança de estado, o que constituía para ele «o equivalente mecânico do calor», tal como explicámos na lição 2. Joule foi assim o primeiro a mostrar (uma indicação experimental do facto) que o trabalho adiabático não depende do caminho seguido.

Foram então Carathéodory, em 1909 [2], e sobretudo Born, em 1921, que propuseram erigir este facto experimental em postulado. Admitindo-o, torna-se fácil construir uma função de estado UU pondo

U(B)U(A)=Wad(AB),U(B)-U(A)=W_{\mathrm{ad}}(A\to B),

e definir em seguida o calor por diferença para as transformações gerais, tal como fizemos.

Esta construção é retomada em numerosas obras, nomeadamente as de Pippard [3] e de Callen [4]. É equivalente à nossa: postular UU e deduzir depois a independência do trabalho adiabático relativamente ao caminho equivale a postular essa independência e construir em seguida UU. Encontra-se em Rosenberg [5] o relato detalhado desta evolução conceptual, de Joule a Carathéodory e Born.

6.2. Micro e macroestados

Se não se optar por seguir a via operacional descrita acima, a secção 2 trouxe à luz uma questão fundamental: como pode uma energia definida num espaço de 6N6N coordenadas resumir-se a uma função de apenas algumas variáveis macroscópicas? A física estatística constrói esta passagem entre as duas escalas. Eis as suas grandes linhas.

Em física estatística, chama-se macroestado a um ponto XX do espaço dos estados de equilíbrio E\mathcal E. A um mesmo macroestado corresponde em geral um número gigantesco de microestados γΓ\gamma\in\Gamma compatíveis com os mesmos constrangimentos macroscópicos.

Demonstra-se então que as energias desses microestados se distribuem de forma extremamente concentrada em torno do seu valor médio quando NN\to\infty, limite dito termodinâmico. Por outras palavras, as flutuações relativas da energia tendem para zero:

σ(Umicro)Umicro0.\frac{\sigma(U_{\mathrm{micro}})} {\langle U_{\mathrm{micro}}\rangle}\longrightarrow0.

Pode então pôr-se Uth(X)=UmicroXU_{\mathrm{th}}(X)= \langle U_{\mathrm{micro}} \rangle_X, sendo a média tomada sobre todos os microestados compatíveis com o macroestado. É a concentração acima evocada que dá o seu sentido a esta definição. O sistema ocupa em cada instante um único microestado, e não a média: é apenas porque quase todos têm a mesma energia que se pode falar da energia do macroestado, e que duas preparações idênticas dão a mesma medida.

Para se representar as coisas, pode imaginar-se que a um macroestado corresponde um conjunto muito grande de microestados microscopicamente diferentes mas macroscopicamente indistinguíveis. Esta imagem assemelha-se à de uma classe de equivalência no espaço de fases. A construção precisa será, contudo, ligeiramente diferente e fará intervir uma distribuição de probabilidade no espaço de fases.

A figura seguinte resume este procedimento.

Um mesmo macroestado X é compatível com um número imenso de microestados. À nossa escala, as suas energias estão extremamente concentradas em torno de um mesmo valor macroscópico, o que permite definir U_ th(X).
Figura 3. Um mesmo macroestado XX é compatível com um número imenso de microestados. À nossa escala, as suas energias estão extremamente concentradas em torno de um mesmo valor macroscópico, o que permite definir Uth(X)U_{\mathrm{th}}(X).

6.3. Trabalhos generalizados

As diferentes formas de trabalho encontradas acima possuem uma estrutura comum. Vimos as formas

δW=PextdV\delta W=-P_{\mathrm{ext}}\,dV

para o trabalho de pressão,

δW=FextdL\delta W=F_{\mathrm{ext}}\,dL

para o estiramento de um fio, ou ainda

δW=γextdA\delta W=\gamma_{\mathrm{ext}}\,dA

para um trabalho de superfície. Nota-se que, de cada vez, a grandeza diferenciada é extensiva ao passo que o seu prefator é intensivo. Generalizando, escreveremos qualquer termo de trabalho mecânico na forma:

δW=iνiextdxi,\boxed{ \delta W=\sum_i \nu_i^{\mathrm{ext}}\,dx^i, }

onde νiext\nu_i^{\mathrm{ext}} designa a força generalizada exterior (em geral intensiva) conjugada da coordenada xix^i (em geral extensiva), com o sinal correspondente à nossa convenção.

No caso em que as forças generalizadas exteriores se identificam com as forças termodinâmicas correspondentes do sistema (exemplo: Pext=PP_{\mathrm{ext}}=P), o primeiro princípio pode então escrever-se na forma

dU=δQ+iνidxi.\boxed{ dU=\delta Q+\sum_i\nu_i\,dx^i.}

6.4. A extensividade da energia interna

Apresentámos na lição anterior a energia interna como uma grandeza extensiva: multiplicar o tamanho de um sistema homogéneo por um fator λ\lambda, a variáveis intensivas fixas, multiplica a sua energia por λ\lambda,

UλU.U\longrightarrow\lambda U.

Esta propriedade não é postulada pelo primeiro princípio. É uma hipótese suplementar que faremos muitas vezes, mas que nem sempre se verifica. Espera-se, com efeito, que a existência de forças de longo alcance entre os constituintes do sistema (nomeadamente gravíticas) arruíne a sua extensividade.

Precisemos isto. Consideremos NN constituintes distribuídos a densidade constante num espaço de dimensão DD, e suponhamos que a sua energia potencial de interação se comporta como

Ep(r)1rα.E_p(r)\sim\frac{1}{r^\alpha}.

A densidade ρ=N/V\rho=N/V fixa, o tamanho linear típico do sistema cresce como

LN1/D.L\sim N^{1/D}.

Avaliemos então a energia de interação total em três etapas.

Contar os vizinhos.

Fixemos um constituinte e perguntemo-nos quantos outros há a uma distância compreendida entre rr e r+drr+dr. É a densidade multiplicada pelo volume da coroa correspondente:

dn(r)=ρSDrD1dr,dn(r)=\rho\,S_D\,r^{D-1}\,dr,

onde SDS_D designa a área da esfera unitária em dimensão DD, ou seja 4π4\pi em três dimensões.

Somar sobre as distâncias.

Cada um desses vizinhos contribui com Ep(r)rαE_p(r)\sim r^{-\alpha}. A energia de interação de um único constituinte com todos os outros vale, pois,

uρaLrD11rαdr=ρaLrD1αdr.u\sim\rho\int_a^L r^{D-1}\,\frac{1}{r^{\alpha}}\,dr =\rho\int_a^L r^{D-1-\alpha}\,dr .

O limite inferior aa é a distância mínima de aproximação, abaixo da qual a lei em 1/rα1/r^\alpha deixa de ser válida e que impede o integral de divergir quando r0r\to0. O limite superior LL é o tamanho do sistema: não há vizinhos para além dele.

Somar sobre os constituintes.

Multiplicamos pelo número de constituintes (dividindo por dois para não contar cada par duas vezes):

UintN2uNρaLrD1αdr,U_{\mathrm{int}}\sim\frac{N}{2}\,u \sim N\rho \int_a^L r^{D-1-\alpha}\,dr,

onde se absorveu o fator geométrico SD/2S_D/2 na ordem de grandeza. Olhemos então para o valor do integral:

aLrD1αdr=[rDαDα]aL(αD).\int_a^L r^{D-1-\alpha}\,dr =\left[\frac{r^{D-\alpha}}{D-\alpha}\right]_a^L \qquad(\alpha\neq D).

Se α>D\alpha>D, o expoente é negativo: o integral converge quando LL\to\infty e é dominado pelo seu limite inferior, sendo portanto da ordem da constante aDα/(αD)a^{D-\alpha}/(\alpha-D). Se α<D\alpha<D, é pelo contrário dominado pelo seu limite superior e é da ordem de LDα=N1α/DL^{D-\alpha}=N^{1-\alpha/D}. Por fim, no caso α=D\alpha=D, a primitiva é um logaritmo e o integral vale ln(L/a)=1DlnN\ln(L/a)=\tfrac1D\ln N. Substituindo, obtém-se assim

Uintcste×{N,α>D,NlnN,α=D,N2α/D,α<D.\begin{aligned} U_{\mathrm{int}} \sim \text{cste} \times \begin{cases} N, & \alpha>D,\\[1mm] N\ln N, & \alpha=D,\\[1mm] N^{\,2-\alpha/D}, & \alpha<D. \end{cases} \end{aligned}

Fisicamente, a interpretação é clara. Quando α>D\alpha>D, temos forças de curto alcance: o decrescimento de EpE_p prevalece sobre o crescimento do número de vizinhos. Cada constituinte sente apenas a sua vizinhança imediata, de modo que a sua energia não depende do tamanho do sistema, e o total é proporcional a NN, o que é extensivo.

Quando αD\alpha\leq D, dá-se o inverso: os vizinhos distantes, muito mais numerosos, prevalecem, cada constituinte sente o sistema inteiro, e a sua própria energia cresce com o tamanho deste. O total cresce então mais depressa do que NN, e a extensividade perde-se. Interações suficientemente curtas conduzem, portanto, naturalmente a uma energia extensiva, ao passo que interações de longo alcance podem destruir essa propriedade.

Na gravitação newtoniana tridimensional, por exemplo,

Ep(r)1r,D=3,α=1,E_p(r)\sim-\frac1r, \qquad D=3, \qquad \alpha=1,

e a estimativa anterior dá, a densidade fixa,

UgravN5/3.U_{\mathrm{grav}}\sim-N^{5/3}.

O sinal negativo, que o argumento de ordem de grandeza não fornece por si só, traduz o carácter atrativo da gravitação.

A extensividade da energia será massivamente utilizada no resto desta obra. A gravitação é, evidentemente, universal, mas é uma força muito fraca: se estudarmos dois volumes de gás postos em contacto, podemos desprezá-la na prática. Por outro lado, as forças internas do gás, do tipo van der Waals, decrescem muito rapidamente (como 1/r61/r^6), de modo que a extensividade é, nesse caso, quase exata.

Em contrapartida, para descrever os chamados sistemas autogravitantes (estrelas, galáxias, ...), a gravitação é evidentemente o ingrediente essencial, e toda a análise termodinâmica deve ser refeita desde o início, uma vez que se perde necessariamente a extensividade da energia. Por isso, a termodinâmica desses sistemas é quase uma ciência à parte e mostra comportamentos inesperados: por exemplo, uma estrela que irradia energia aquece em vez de arrefecer. Voltaremos a isto na parte avançada da obra.

7. Síntese

Resumamos o que esta lição estabeleceu e que servirá em tudo o que se segue.

  • O primeiro princípio postula a existência de uma função de estado U:ERU:\mathcal E\longrightarrow\mathbb R, a energia interna, e afirma a conservação da energia de um sistema isolado.
  • No resto da obra supô-la-emos também extensiva, e suficientemente contínua e diferenciável em todas as suas variáveis.
  • Completado pela aditividade da energia, o primeiro princípio permite formular a conservação como um balanço energético entre um sistema fechado e o seu ambiente. Sendo o trabalho definido independentemente pela mecânica, o calor é então definido como a transferência restante, o que conduz à útil fórmula ΔU=Q+W\Delta U=Q+W.
  • Para uma transformação quase-estática, este balanço de energia toma a forma diferencial dU=δQ+δWdU=\delta Q+\delta W, na qual apenas dUdU é uma diferencial exata.

A análise dos sistemas abertos será discutida noutro ponto do curso.

8. Referências

Sobre a evolução conceptual do primeiro princípio, de Joule a Carathéodory e Born, ver Rosenberg [5]. Para uma apresentação clássica, ver Pippard [3] e Callen [4].

  1. J. P. Joule, “On the Mechanical Equivalent of Heat,” Phil. Trans. R. Soc. Lond. 140, 61—82 (1850)
  2. C. Carathéodory, “Untersuchungen über die Grundlagen der Thermodynamik,” Math. Ann. 67, 355—386 (1909)
  3. A. B. Pippard, Elements of Classical Thermodynamics for Advanced Students of Physics, Cambridge University Press (1957)
  4. H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2nd ed., Wiley (1985)
  5. R. M. Rosenberg, “From Joule to Caratheodory and Born: A Conceptual Evolution of the First Law of Thermodynamics,” J. Chem. Educ. 87, 691—693 (2010)
  6. A. M. Steane, “First Law, internal energy,” chap. 7 in Thermodynamics: A Complete Undergraduate Course, Oxford University Press (2017)
  7. E. A. Gislason and N. C. Craig, “Cementing the foundations of thermodynamics: Comparison of system-based and surroundings-based definitions of work and heat,” J. Chem. Thermodynamics 37, 954—966 (2005)