Computing equilibrium measures with power law kernels
We introduce a method to numerically compute equilibrium measures for problems with attractive-repulsive power law kernels of the form $K(x-y) = \frac{|x-y|^\alpha}{\alpha}-\frac{|x-y|^\beta}{\beta}$ using recursively generated banded and approximately banded operators acting on expansions in ultras...
| Main Authors: | , , |
|---|---|
| Format: | Journal article |
| Language: | English |
| Published: |
American Mathematical Society
2022
|
| _version_ | 1797107612821487616 |
|---|---|
| author | Gutleb, TS Carrillo de la Plata, JA Olver, S |
| author_facet | Gutleb, TS Carrillo de la Plata, JA Olver, S |
| author_sort | Gutleb, TS |
| collection | OXFORD |
| description | We introduce a method to numerically compute equilibrium measures for problems with attractive-repulsive power law kernels of the form $K(x-y) = \frac{|x-y|^\alpha}{\alpha}-\frac{|x-y|^\beta}{\beta}$ using recursively generated banded and approximately banded operators acting on expansions in ultraspherical polynomial bases. The proposed method reduces what is naively a difficult to approach optimization problem over a measure space to a straightforward optimization problem over one or two variables fixing the support of the equilibrium measure. The structure and rapid convergence properties of the obtained operators results in high computational efficiency in the individual optimization steps. We discuss stability and convergence of the method under a Tikhonov regularization and use an implementation to showcase comparisons with analytically known solutions as well as discrete particle simulations. Finally, we numerically explore open questions with respect to existence and uniqueness of equilibrium measures as well as gap forming behaviour in parameter ranges of interest for power law kernels, where the support of the equilibrium measure splits into two intervals. |
| first_indexed | 2024-03-07T07:18:28Z |
| format | Journal article |
| id | oxford-uuid:7fc4ddcd-2581-42d5-b5c4-69fbc917e28a |
| institution | University of Oxford |
| language | English |
| last_indexed | 2024-03-07T07:18:28Z |
| publishDate | 2022 |
| publisher | American Mathematical Society |
| record_format | dspace |
| spelling | oxford-uuid:7fc4ddcd-2581-42d5-b5c4-69fbc917e28a2022-09-01T14:21:50ZComputing equilibrium measures with power law kernels Journal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:7fc4ddcd-2581-42d5-b5c4-69fbc917e28aEnglishSymplectic ElementsAmerican Mathematical Society2022Gutleb, TSCarrillo de la Plata, JAOlver, SWe introduce a method to numerically compute equilibrium measures for problems with attractive-repulsive power law kernels of the form $K(x-y) = \frac{|x-y|^\alpha}{\alpha}-\frac{|x-y|^\beta}{\beta}$ using recursively generated banded and approximately banded operators acting on expansions in ultraspherical polynomial bases. The proposed method reduces what is naively a difficult to approach optimization problem over a measure space to a straightforward optimization problem over one or two variables fixing the support of the equilibrium measure. The structure and rapid convergence properties of the obtained operators results in high computational efficiency in the individual optimization steps. We discuss stability and convergence of the method under a Tikhonov regularization and use an implementation to showcase comparisons with analytically known solutions as well as discrete particle simulations. Finally, we numerically explore open questions with respect to existence and uniqueness of equilibrium measures as well as gap forming behaviour in parameter ranges of interest for power law kernels, where the support of the equilibrium measure splits into two intervals. |
| spellingShingle | Gutleb, TS Carrillo de la Plata, JA Olver, S Computing equilibrium measures with power law kernels |
| title | Computing equilibrium measures with power law kernels |
| title_full | Computing equilibrium measures with power law kernels |
| title_fullStr | Computing equilibrium measures with power law kernels |
| title_full_unstemmed | Computing equilibrium measures with power law kernels |
| title_short | Computing equilibrium measures with power law kernels |
| title_sort | computing equilibrium measures with power law kernels |
| work_keys_str_mv | AT gutlebts computingequilibriummeasureswithpowerlawkernels AT carrillodelaplataja computingequilibriummeasureswithpowerlawkernels AT olvers computingequilibriummeasureswithpowerlawkernels |