Uniform distributions in nonuniform systems: Wall potentials generating constant density profiles in classical density functional theory

Resumo Técnico: Distribuições Uniformes em Sistemas Não Uniformes

Enunciado do Problema
Este trabalho aborda o problema inverso dentro da Teoria do Funcional de Densidade (DFT) Clássica para fluidos inhomogêneos. Enquanto as aplicações padrão de DFT resolvem o "problema direto" — determinando o perfil de densidade de equilíbrio $\rho(r)$ para um potencial externo prescrito $V(r)$ — este estudo busca o inverso: construir o potencial de parede externa específico $V(r)$ necessário para gerar um perfil de densidade de equilíbrio perfeitamente plano e constante $\rho(r) = \rho_b$ na região do fluido adjacente a um substrato.

A motivação decorre da observação de que substratos realistas tipicamente induzem oscilações de empacotamento, camadas de depleção ou adsorção. Os autores visam identificar campos de parede de "cancelamento de estrutura" que suprimem explicitamente essas inhomogeneidades. Este problema está ligado às teorias de interface de molhabilidade (wetting) e secagem (drying), onde perfis de densidade planos surgem sob condições especiais de correspondência entre as interações parede-fluido e fluido-fluido. O estudo investiga como o campo compensador necessário depende da geometria do substrato (planar, esférica e cilíndrica) e das interações do fluido.

Metodologia
Os autores empregam a Teoria de Medida Fundamental (FMT) de Rosenfeld para descrever o sistema de referência de esferas rígidas e um tratamento de campo médio para interações atrativas. Dois modelos de um único componente são analisados:

1. Um fluido de esferas rígidas puramente repulsivo.
2. Um fluido com repulsão de esferas rígidas mais uma cauda atrativa do tipo Lennard-Jones truncada.

O estudo considera três geometrias fundamentais:

• Planar: Uma parede plana em $z=R$.
• Esférica: Uma esfera sólida de raio $\varrho$.
• Cilíndrica: Um cilindro sólido de raio $\varrho$ infinitamente longo.

A metodologia central envolve a substituição do perfil de densidade constante prescrito na equação de Euler-Lagrange derivada do funcional do potencial grand canonical. O potencial externo é recuperado via:
$$\beta V(\xi) = \beta\mu + c^{(1)}(\xi) - \ln[\Lambda^3 \rho(\xi)]$$
onde $c^{(1)}$ é a função de correlação direta de um corpo. Os autores derivam expressões explícitas para as densidades ponderadas e $c^{(1)}$ para cada geometria.

• Tratamento Analítico: Para geometrias plana e esférica, as integrais de convolução da FMT são reduzidas a formas unidimensionais que permitem soluções analíticas fechadas.
• Tratamento Numérico: Para a geometria cilíndrica, a redução dos núcleos da FMT leva a integrais elípticas completas, impedindo uma forma fechada simples; portanto, este caso é avaliado via quadratura numérica.

Resultos Principais

1. Geometria Planar: O problema inverso produz expressões analíticas compactas e explícitas para o potencial de parede. Para esferas rígidas, o potencial é puramente repulsivo e de curto alcance, com uma amplitude que aumenta com a fração de empacotamento $\eta$. Quando as interações atrativas são incluídas, o potencial desenvolve um componente atrativo, refletindo a competição entre empacotamento e coesão. O valor de contato do potencial satisfaz a regra de soma exata $\beta p_{hs}/\rho_b = 1 + \beta V_{hs}^{pl}(R)$.
2. Geometria Esférica: Existem soluções analíticas, mas elas são substancialmente mais complexas, dependendo tanto de $\eta$ quanto do raio do substrato $\varrho$. Os resultados mostram que a curvatura tem um efeito fraco no potencial de esferas rígidas puramente repulsivo, mas um efeito mais pronunciado em fluidos com interações atrativas, onde o poço atrativo torna-se mais raso à medida que o raio da esfera diminui.
3. Geometria Cilíndrica: Os potenciais de parede foram determinados numericamente. Os resultados demonstram que os potenciais cilíndricos situam-se intermediariamente entre os limites planar e esférico para o mesmo raio, refletindo a curvatura intermediária do cilindro.
4. Validação: Cálculos independentes de DFT utilizando os potenciais derivados confirmam que os perfis de densidade planos prescritos são recuperados dentro da precisão numérica (desvios relativos da ordem de $10^{-4}$), validando as construções analíticas e numéricas.

Significância e Alegações
O artigo afirma fornecer uma realização microscópica explícita de campos de parede de cancelamento de estrutura, indo além de inversões puramente formais ou numéricas. As principais contribuições são:

• Construção Analítica Explícita: A derivação de expressões de forma fechada para potenciais de parede em geometrias plana e esférica que cancelam a estrutura de interface.
• Dependência Geométrica: Uma demonstração clara de como a curvatura do substrato modifica o campo compensador, destacando particularmente a dependência de curvatura mais forte em sistemas com interações atrativas em comparação com sistemas puramente repulsivos.
• Material de Referência: O trabalho fornece uma coleção compacta de fórmulas para densidades ponderadas e funções de correlação direta de um corpo nas três geometrias fundamentais, servindo como expressões de referência úteis para implementações de DFT.
• Contexto Teórico: O estudo estabelece que as condições de perfil plano, frequentemente tratadas como condições de correspondência assintóticas ou fenomenológicas em teorias de molhabilidade, podem ser realizadas explicitamente dentro de um funcional de densidade não local microscópico.

Os autores observam que o alcance do potencial compensador é finito, vinculado ao alcance dos núcleos da FMT e à interação atrativa truncada. Eles reconhecem que, além da teoria de campo médio, flutuações críticas introduziriam contribuições de longo alcance, que estão fora do escopo deste tratamento de alcance finito. O trabalho sugere extensões futuras para fluidos multicomponentes, expansões sistemáticas de curvatura e a aplicação deste ponto de vista inverso em contextos dinâmicos e reológicos.