The Anisotropic Interface Continuum Solvation Model and the Finite-Element Anisotropic Poisson Solver
Fuente:
arXiv
Guardado en:
| Autores principales: | , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2025
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866909728256294912 |
|---|---|
| author | Chai, Ziwei Luber, Sandra |
| author_facet | Chai, Ziwei Luber, Sandra |
| contents | We propose an anisotropic interfacial continuum solvation (AICS) model to simulate the distinct in-plane and out-of-plane dielectric constants of liquids near solid-liquid interfaces and their spatial variations along the surface normal direction. In low-electron-density regions, each dielectric function in the diagonal components of a dielectric tensor varies monotonically with distance from the solid surface along the surface normal; in high-electron-density regions near the surface, each dielectric function adopts the electron-density-based formulation proposed by Andreussi et al. (J. Chem. Phys. 136, 064102 (2012)) The resulting dielectric tensor is continuously differentiable with respect to both electron density and spatial coordinates. We derived analytical expressions for electrostatic contributions to the KS potential and forces, and implemented AICS, including these analytical derivatives, into CP2K. To solve the anisotropic Poisson equations, we developed a parallel finite-element anisotropic Poisson solver (FEAPS) based on the FEniCSx platform and its interface with CP2K. Analytical forces were validated against finite-difference calculations, while electrostatic potentials computed under vacuum and isotropic solvent conditions using AICS and FEAPS were benchmarked against standard vacuum DFT and SCCS results, respectively. In the anisotropic solvent environment characterized by the enhanced in-plane and reduced out-of-plane dielectric functions near the Ag(111) surface, we calculated the resulting work functions and electrostatic potentials, and optimized the adsorption geometry for OH. Compared to the isotropic case, we observed more pronounced work function shifts and spatially modulated electrostatic profiles across different charge states. Our results also showed that OH tilted more towards the plane parallel to the surface under the anisotropic dielectric conditions. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_05577 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | The Anisotropic Interface Continuum Solvation Model and the Finite-Element Anisotropic Poisson Solver Chai, Ziwei Luber, Sandra Chemical Physics Materials Science We propose an anisotropic interfacial continuum solvation (AICS) model to simulate the distinct in-plane and out-of-plane dielectric constants of liquids near solid-liquid interfaces and their spatial variations along the surface normal direction. In low-electron-density regions, each dielectric function in the diagonal components of a dielectric tensor varies monotonically with distance from the solid surface along the surface normal; in high-electron-density regions near the surface, each dielectric function adopts the electron-density-based formulation proposed by Andreussi et al. (J. Chem. Phys. 136, 064102 (2012)) The resulting dielectric tensor is continuously differentiable with respect to both electron density and spatial coordinates. We derived analytical expressions for electrostatic contributions to the KS potential and forces, and implemented AICS, including these analytical derivatives, into CP2K. To solve the anisotropic Poisson equations, we developed a parallel finite-element anisotropic Poisson solver (FEAPS) based on the FEniCSx platform and its interface with CP2K. Analytical forces were validated against finite-difference calculations, while electrostatic potentials computed under vacuum and isotropic solvent conditions using AICS and FEAPS were benchmarked against standard vacuum DFT and SCCS results, respectively. In the anisotropic solvent environment characterized by the enhanced in-plane and reduced out-of-plane dielectric functions near the Ag(111) surface, we calculated the resulting work functions and electrostatic potentials, and optimized the adsorption geometry for OH. Compared to the isotropic case, we observed more pronounced work function shifts and spatially modulated electrostatic profiles across different charge states. Our results also showed that OH tilted more towards the plane parallel to the surface under the anisotropic dielectric conditions. |
| title | The Anisotropic Interface Continuum Solvation Model and the Finite-Element Anisotropic Poisson Solver |
| topic | Chemical Physics Materials Science |
| url | https://arxiv.org/abs/2508.05577 |