Progress toward a better BOCS: Systematic coarse-graining with local density potentials

Resumo Técnico: Progresso em direção a um melhor BOCS: Granularização sistemática com potenciais de densidade local

Definição do Problema
Modelos de granularização fina (CG) de baixo para cima (bottom-up) oferecem eficiência computacional para simular escalas de comprimento e tempo inacessíveis aos modelos atômicos detalhados (AA). Embora o princípio variacional de ajuste de força (force-matching - FM) permita a parametrização de potenciais CG para aproximar de forma ideal o potencial CG exato, as implementações padrão frequentemente dependem de potenciais de par convencionais. Esses potenciais de par têm dificuldade em capturar as propriedades termodinâmicas e estruturais de sistemas desiguais, como interfaces líquido-vapor, e frequentemente carecem de transferibilidade entre ambientes de bulk e de interface. Além disso, potenciais de volume, embora eficazes para reproduzir equações de estado de pressão-densidade no bulk, não podem modelar interfaces. Há uma necessidade de um arcabouço sistemático que incorpore potenciais dependentes da densidade local (LD) e do gradiente ao quadrado (SG) dessa densidade para melhorar a fidelidade e a transferibilidade de modelos CG para materiais moles complexos.

Metodologia
Os autores descrevem a versão 5.0 do pacote Bottom-up Open-source Coarse-graining Software (BOCS) e uma extensão correspondente ao pacote PKG-BOCS para LAMMPS.

• Arcabouço Teórico: O software parametriza potenciais CG usando o princípio de ajuste de força, que minimiza a diferença média dos quadrados das forças entre o modelo AA e o modelo CG sobre o conjunto mapeado. O potencial não ligado total é expandido para incluir três componentes:

  1. Potenciais de par centrais convencionais ($U_{pair}$).
  2. Potenciais de densidade local ($U_{LD}$), que dependem da densidade local instantânea $\rho_I$ ao redor de cada partícula $I$.
  3. Potenciais de gradiente ao quadrado ($U_{SG}$), que dependem do quadrado do gradiente da densidade local, $|\nabla_I \rho_I|^2$.
     A densidade local é definida como uma soma de contribuições de partículas vizinhas ponderadas por uma função $w(r)$, para a qual cinco diferentes formas funcionais (ex: DPD, Lucy, Shell, smooth, esfera) são implementadas.

• Implementação de Software:

  • BOCS v5.0: Esta versão introduz ferramentas (ldcalc) para calcular densidades locais e seus gradientes a partir de trajetórias AA. Ela emprega um arcabouço generalizado de Yvon-Born-Green (g-YBG) para realizar o ajuste de força sem informações explícitas de força para termos ligados, de par e de LD. O software utiliza funções de base baseadas em splines para representar funções de força e coeficientes, permitindo formas funcionais flexíveis. Ele resolve as equações normais resultantes via decomposição de valores singulares para determinar os parâmetros ótimos.
  • PKG-BOCS (LAMMPS): Um novo estilo de par (ldd) foi implementado no LAMMPS para simular esses potenciais. Esta implementação suporta topologias moleculares complexas, múltiplos tipos de partículas e várias funções de ponderação. Ela contabiliza adequadamente as contribuições intramoleculares e exclusões.

• Validação e Aplicação: Os autores validaram a implementação através de três testes de referência:

  1. Um sistema de 400 partículas com potenciais LD lineares.
  2. Um sistema de 2 partículas com potenciais LD quadráticos e SG lineares.
  3. Um sistema complexo de 125 moléculas ramificadas com 4 tipos de sítios, envolvendo interações ligadas, de par, de LD e de SG.
     Adicionalmente, o software foi aplicado para parametrizar modelos de água CG de 1 sítio baseados no modelo AA SPC/E, comparando modelos parametrizados a partir de simulações de bulk, simulações de slab e simulações de bulk aumentadas com potenciais de LD.

Principais Resultados

• Precisão de Implementação: Em todos os testes de referência, o BOCS v5.0 recuperou quantitativamente os potenciais e forças de "verdade fundamental" (ground truth) conhecidos usados nas simulações LAMMPS. Os cálculos de ajuste de força produziram erros desprezíveis ($\chi^2 \approx 0$), demonstrando que o software implementa corretamente o princípio variacional para potenciais de LD e SG.
• Percepção Física: O estudo revelou que as interações de par efetivas geradas por potenciais de LD e SG são altamente sensíveis à escolha da função de ponderação $w(r)$. Por exemplo, diferentes funções de ponderação podem levar a potenciais de par e perfis de força significativamente diferentes, mesmo quando a forma funcional de LD ou SG subjacente é idêntica.
• Desempenho do Modelo de Água:
  • Um modelo padrão apenas de par, parametrizado a partir de água em bulk (PAIR-BULK), falhou em estabilizar a fase líquida a 1 bar, resultando em vaporização.
  • Um modelo apenas de par, parametrizado a partir de uma simulação de slab (PAIR-SLAB), estabilizou o líquido, mas superestimou significativamente a densidade e a compressibilidade e reproduziu mal a largura da interface líquido-vapor.
  • O modelo LD-BULK (par + potenciais de LD parametrizados a partir de bulk) melhorou significativamente o desempenho. Ele reproduziu a equação de estado do SPC/E com alta precisão (subestimando a densidade em apenas 0,1%) e reproduziu razoavelmente a largura da interface líquido-vapor e o perfil de densidade, apesar de ter sido parametrizado unicamente a partir de uma simulação de bulk.
  • A análise mostrou que o potencial de LD modula as forças de par com base no ambiente local, efetivamente aprofundando as interações atrativas em regiões de baixa densidade (imitando o comportamento de interface) enquanto permanece repulsivo em regiões de alta densidade.

Significância e Alegações
Os autores afirmam que o BOCS versão 5.0 e a extensão PKG-BOCS fornecem uma ferramenta robusta e precisa para parametrizar sistematicamente modelos CG com potenciais de densidade local e de gradiente ao quadrado. Estes potenciais oferecem uma melhoria significativa na fidelidade estrutural, propriedades termodinâmicas e transferibilidade de modelos CG de baixo para cima, particularmente para sistemas desiguais como interfaces líquido-vapor. O trabalho demonstra que os potenciais de LD podem capturar efetivamente interações de muitos corpos cooperativas e melhorar a descrição das flutuações de densidade sem exigir termos de força de muitos corpos explícitos na simulação, uma vez que a dependência ambiental é capturada através da densidade local. Os autores posicionam este software como um complemento útil aos pacotes existentes para o desenvolvimento de modelos CG de materiais moles complexos.