Next Article in Journal
Microsatellites as Agents of Adaptive Change: An RNA-Seq-Based Comparative Study of Transcriptomes from Five Helianthus Species
Next Article in Special Issue
Identifying the Locations of Atmospheric Pollution Point Source by Using a Hybrid Particle Swarm Optimization
Previous Article in Journal
Ensemble Empirical Mode Decomposition with Adaptive Noise with Convolution Based Gated Recurrent Neural Network: A New Deep Learning Model for South Asian High Intensity Forecasting
Previous Article in Special Issue
Multinomial Logit Model Building via TreeNet and Association Rules Analysis: An Application via a Thyroid Dataset
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Properties of Discrete Minimal s,logt-Energy Constants and Configurations

by
Nichakan Loesatapornpipit
1,†,‡ and
Nattapong Bosuwan
1,2,*,‡
1
Department of Mathematics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand
2
Centre of Excellence in Mathematics, Commission on Higher Education, Ministry of Education, Si Ayutthaya Road, Bangkok 10400, Thailand
*
Author to whom correspondence should be addressed.
Results of this article constitute part of Nichakan Loesatapornpipit’s senior project under the mentorship of Nattapong Bosuwan at Mahidol University.
These authors contributed equally to this work.
Submission received: 18 April 2021 / Revised: 18 May 2021 / Accepted: 20 May 2021 / Published: 24 May 2021
(This article belongs to the Special Issue Modelling and Simulation of Natural Phenomena of Current Interest)

Abstract

:
We investigated the energy of N points on an infinite compact metric space ( A , d ) of a diameter less than 1 that interact through the potential ( 1 / d s ) ( log 1 / d ) t , where s , t 0 and d is the metric distance. With E log t s ( A , N ) denoting the minimal energy for such N-point configurations, we studied certain continuity and differentiability properties of E log t s ( A , N ) in the variable s . Then, we showed that in the limits, as s and as s s 0 > 0 , N-point configurations that minimize the s , log t -energy tends to an N-point best-packing configuration and an N-point configuration that minimizes the s 0 , log t -energy, respectively. Furthermore, we considered when A are circles in the Euclidean space R 2 . In particular, we proved the minimality of N distinct equally spaced points on circles in R 2 for some certain s and t. The study on circles shows a possibility for the utilization of N points generated through such new potential to uniformly discretize on objects with very high symmetry.

1. Introduction

The general setting of the discrete minimal energy problem is the following. Let ( A , d ) be an infinite compact metric space and K : A × A R { } be a lower semicontinuous kernel. Note that in some contexts, the kernel K ( x , y ) is called a potential. For a fixed set of N points ω N A , we define the K-energy of ω N as follows:
E K ( ω N ) : = x y x , y ω N K ( x , y ) .
The minimal N-point K-energy of the set A is defined by
E K ( A , N ) : = min ω N A # ω N = N E K ( ω N ) ,
where # ω N stands for the cardinality of the set ω N . A minimal N-point K-energy configuration is a configuration ω N K of N points in A that minimizes such energy, namely
E K ( ω N K ) = min ω N A # ω N = N E K ( ω N ) .
It is known that ω N K always exists and in general ω N K may not be unique.
Two important kernels in the theory of minimal energy are Riesz and logarithmic kernels. The (Riesz) s-kernel and log-kernel are defined by
K s ( x , y ) : = 1 d ( x , y ) s , s 0 ,
and
K log ( x , y ) : = log 1 d ( x , y ) ,
for all ( x , y ) A × A , respectively. It is not difficult to check that both kernels are lower semicontinuous on A × A . The s-energy of ω N and the minimal N-point s-energy of the set A are
E s ( ω N ) : = E K s ( ω N ) and E s ( A , N ) : = min ω N A # ω N = N E s ( ω N )
and we denote by ω N s : = ω N K s and call this configuration a minimal N-point s-energy configuration. Similarly, the log-energy of ω N and the N-point log-energy of the set A are
E log ( ω N ) : = E K log ( ω N ) and E log ( A , N ) : = min ω N A # ω N = N E log ( ω N )
and we denote by ω N log : = ω N K log and call this configuration a minimal N-point log-energy configuration.
Let us provide a short survey of these two energy problems.
The study of s-energy constants and configurations has a long history in physics, chemistry, and mathematics. Finding the arrangements of ω N s where the set A is the unit sphere S 2 in the Euclidean space R 3 has been an active area since the beginning of the 19th century. The problem is known as the generalized Thomson problem (see [1] and [2] (Chapter 2.4)). Candidates for ω N s for several numbers of N are available (see, e.g., [3]). However, the solutions (with rigorous proof) are obtainable for a handful of values of N (see, e.g., [4,5,6,7]). For example, when N = 5 , the generalized Thomson problem becomes surprisingly difficult. Schwartz, using computer-aided proof, showed that such ω 5 s on S 2 can be the vertices of the triangular bipyramid or a square-based pyramid (depending on s) in a single monograph of 180 pages [8] (see also a synopsis of his work [7]). For a general compact set A in the Euclidean space R m , the study of the distribution of a minimal N-point s-energy configurations of A as N can be found in [9,10]. In [10], it was shown that when s is any fixed number greater than the Hausdorff dimension of A , the minimal N-point s-energy configurations of A are “good points” to represent the set A . This is because such configurations are asymptotically uniformly distributed over the set A as N (see the precise statement in [10] (Theorems 2.1 and 2.2)). The results in [10] have wide ranging applications in various fields of computational science, such as computer-aided geometric design, finite element tesselations, statistical sampling, etc.
The log-energy problem has been heavily studied when A is a subset of the Euclidean space R 2 (or C ) because it has had a profound influence in approximation theory (see, e.g., [11,12,13,14,15]). For A C , the points in ω N log are commonly known as Fekete points or Chebyshev points which can be used as interpolation and integration nodes (see [16]). The log-energy problem received other special attention when Steven Smale posed Problem #7 in his book chapter entitled “Mathematical problems for the next century” [17]. Problem #7 asks for a construction of an algorithm which on input N 2 outputs a configuration ω N = { x 1 , , x N } of distinct points on S 2 embedded in R 3 such that
E log ( ω N ) E log ( S 2 , N ) c log N
(where c is a constant independent of N and ω N ) and requires that its running time grows at most polynomially in N . This problem arose from his joint work with Shub [18] on complexity theory. In order to answer this question, it is natural to understand the asymptotic expansion of E log ( S 2 , N ) in the variable N (see [19] for conjectures and the progress). The problem concerning the arrangements of ω N log on the unit sphere S 2 in R 3 is posed by Whyte [20] in 1952. Whyte’s problem is also attractive and intractable. We refer to [21] for a glimpse of this problem.
In [2], Borodachov, Hardin, and Saff investigated asymptotic properties of minimal N-point s-energy constants and configurations for fixed N and varying s . Since this is our main focus in this paper, we will state these results below.
The first theorem (Ref. [2] (Theorem 2.7.1 and Theorem 2.7.3)) concerns the continuity and differentiability of the function
f ( s ) : = E s ( A , N ) , s 0 .
In order to state such a theorem, let us define a set
G log s ( A , N ) : = x y x , y ω N K s ( x , y ) K log ( x , y ) : ω N A   and   E s ( ω N ) = E s ( A , N ) ,
for s 0 .
Theorem 1.
Let ( A , d ) be an infinite compact metric space and let N 2 be fixed. Then:
(a) 
the function f ( s ) defined in (2) is continuous on [ 0 , ) ;
(b) 
the function f ( s ) is right differentiable on [ 0 , ) and left differentiable on ( 0 , ) with:
f + ( s ) : = lim r s + f ( r ) f ( s ) r s = inf G s ( A , N ) , s 0 ,
and:
f ( s ) : = lim r s f ( r ) f ( s ) r s = sup G s ( A , N ) , s > 0 .
We will see in Theorems 2 and 3 below that there are certain relations between minimal s-energy problems as s and the best-packing problem defined as follows. The N-point best-packing distance of the set A is defined as
δ N ( A ) : = max { δ ( ω N ) : ω N A } ,
where
δ ( ω N ) : = min 1 i j N d ( x i , x j )
denotes the separation distance of an N-point configuration ω N = { x 1 , , x N } , and N-point best-packing configurations are N-point configurations attaining the maximum in (3). For further details on the best-packing problem, we refer the reader to [2] (Chapter 3).
The following theorem [2] (Corollary 2.7.5 and Proposition 3.1.2) explains the behavior of E s ( A , N ) as s 0 + and s .
Theorem 2.
For N 2 and an infinite compact metric space ( A , d ) ,
lim s 0 + E s ( A , N ) N ( N 1 ) s = E log ( A , N )
and
lim s E s ( A , N ) 1 / s = 1 δ N ( A ) .
Before we state more results, let us define a cluster configuration. Let s 0 [ 0 , ] . We say that
  • An N-point configuration ω N A is a cluster configuration of ω N s as s s 0 + if there is a sequence { s k } k = 1 ( s 0 , ) such that lim k s k = s 0 and lim k ω N s k = ω N in the topology of A N induced by the metric d.
  • An N-point configuration ω N A is a cluster configuration of ω N s as s s 0 if there is a sequence { s k } k = 1 [ 0 , s 0 ) such that lim k s k = s 0 and lim k ω N s k = ω N in the topology of A N induced by the metric d.
  • An N-point configuration ω N A is a cluster configuration of ω N s as s s 0 if there is a sequence { s k } k = 1 [ 0 , ) such that lim k s k = s 0 and lim k ω N s k = ω N in the topology of A N induced by the metric d.
The properties of the cluster configurations of minimal N-point s-energy configurations as s varies (see [2] (Theorem 2.7.1 and Proposition 3.1.2)) are described i Theorem 3.
Theorem 3.
Let ( A , d ) be an infinite compact metric space and for s 0 and N 2 , let ω N s denote a minimal N-point s-energy configuration on A . Then,
(a) 
For s 0 > 0 , any cluster configuration of ω N s as s s 0 is a minimal N-point s 0 -energy configuration;
(b) 
Any cluster configuration of ω N s as s 0 + is a minimal N-point log-energy configuration;
(c) 
Any cluster configuration of ω N s as s is a N-point best-packing configuration.
In this paper, we consider the following s , log t -kernel:
K log t s ( x , y ) = 1 d ( x , y ) s log 1 d ( x , y ) t , s 0 , t 0 .
with a corresponding s , log t -energy of ω N and minimal N-point s , log t -energy of the set A:
E log t s ( ω N ) : = E K log t s ( ω N ) and E log t s ( A , N ) : = min ω N A # ω N = N E log t s ( ω N ) ,
respectively. We set
ω N s , log t : = ω N K log t s ,
and call it a minimal N-point s , log t -energy configuration. Note that the kernel K log t s ( x , y ) is lower semicontinuous on A × A and this s , log t -energy can be viewed as a generalization of both s-energy and log-energy. The kernel in (4) first appeared in the study of the differentiability of the function f ( s ) = E s ( A , N ) in [2] (Theorem 2.7.3). To the authors’ knowledge, no study involving s , log t -energy constants and configurations has appeared in the literature previously.
The first goal of this paper was to prove the analogues of Theorems 1–3 for s , log t -energy constants and configurations. We would like to emphasize that we will limit our interest to the sets A with diam ( A ) < 1 , where
diam ( A ) : = sup x , y A d ( x , y )
denotes the diameter of A . For the case where diam ( A ) 1 , the values of the kernel K log t s ( x , y ) can be 0 or negative and the analysis becomes laborious.
The second goal was to investigate the arrangements of ω N s , log t on circles in R 2 . Using an available tool in (Ref. [2] (Theorem 2.3.1)), we show that ω N s , log t on any circle with a diameter less than 1 are N distinct equally spaced points. The motivation of this study for objects with very high symmetry comes from the study of the limiting distributions of ω N s , s > 0 and ω N log on the m-dimensional sphere S m in the Euclidean space R m + 1 in [2] (Theorem 6.1.7). In [2] (Theorem 6.1.7), it was shown that ω N s , s > 0 and ω N log on S m are asymptotically uniformly distributed with respect to the surface area measured on S m as N (see also [22] for applications of this result). Our study on circles exhibits a possibility for the utilization of ω N s , log t to uniformly discretize m-dimensional spheres in R m + 1 .
The remainder of this article is organized as follows. The main results of this paper are stated in Section 2. The proof of the main results are in Section 3. We keep all auxiliary lemmas in Section 4. Finally, conclusions and future work are discussed in Section 5.

2. Main Results

2.1. Asymptotic Properties of Discrete Minimal s , log t -Energy

The asymptotic behavior of minimal N-point s , log t -energy constants and configurations as s can be explained in the following theorem.
Theorem 4.
Let N 2 and t 0 be fixed. Assume that ( A , d ) is an infinite compact metric space with diam ( A ) < 1 . Then,
lim s E log t s ( A , N ) 1 / s = 1 δ N ( A ) .
Furthermore, every cluster configuration of ω N s , log t as s is an N-point best-packing configuration on A.
For a fixed t 0 , we define
g ( s ) : = E log t s ( A , N ) , s 0 .
The continuity of g ( s ) is stated below.
Theorem 5.
Let N 2 and t 0 be fixed. Assume that ( A , d ) is an infinite compact metric space with diam ( A ) < 1 . Then, the function g ( s ) is continuous on [ 0 , ) .
As a consequence of the continuity of g ( s ) , we analyze a property of cluster configurations of ω N s , log t as s s 0 > 0 in the following theorem.
Theorem 6.
Let N 2 and t 0 be fixed. Assume that ( A , d ) is an infinite compact metric space with diam ( A ) < 1 . Denote by ω N s , log t a minimal N-point s , log t -energy configuration on A. Then, for any s 0 > 0 , any cluster configuration of ω N s , log t , as s s 0 , is a minimal N-point s 0 , log t -energy configuration on A.
For s 0 and t 0 , we set
G log t + 1 s ( A , N ) : = { E log t + 1 s ( ω N ) : ω N A and E log t s ( ω N ) = E log t s ( A , N ) } .
The differentiability properties of g ( s ) are described in Theorems 7 and 8.
Theorem 7.
Let N 2 and t 0 be fixed. Assume that ( A , d ) is an infinite compact metric space with diam ( A ) < 1 . Then, the function g ( s ) is right differentiable on [ 0 , ) and left differentiable on ( 0 , ) with
g + ( s ) : = lim r s + g ( r ) g ( s ) r s = inf G log t + 1 s ( A , N ) , s 0 ,
and
g ( s ) : = lim r s g ( r ) g ( s ) r s = sup G log t + 1 s ( A , N ) , s > 0 .
Observing that
inf G log 0 ( A , N ) = E log ( A , N ) and g ( 0 ) = N ( N 1 ) ,
when s , t = 0 , Theorem 7 simply implies that
lim s 0 + E s ( A , N ) N ( N 1 ) s = E log ( A , N ) .
Theorem 8.
Let N 2 and t 0 be fixed. Assume that ( A , d ) is an infinite compact metric space with diam ( A ) < 1 . Then,
(a) 
The function g ( s ) is differentiable at s = s 0 > 0 if and only if
inf G log t s 0 ( A , N ) = sup G log t s 0 ( A , N ) ;
(b) 
If ω N * is a cluster point of ω N s , log t as s s 0 + 0 , then
E log t + 1 s 0 ( ω N * ) = inf G log t + 1 s 0 ( A , N ) = g + ( s 0 ) ;
(c) 
If ω N * * is a cluster point of ω N s , log t as s s 0 > 0 , then
E log t + 1 s 0 ( ω N * * ) = sup G log t + 1 s 0 ( A , N ) = g ( s 0 ) ;
(d) 
For s 0 > 0 , if there exists a configuration ω N * that is both cluster configurations of ω N s , log t as s s 0 + and s s 0 , then the function g ( s ) is differentiable at s = s 0 with
E log t + 1 s 0 ( ω N * ) = g ( s 0 ) .
In Theorem 8, we provide criteria for the differentiability of g ( s ) . In particular, the part (a) in Theorem 8 implies that if all minimal N-point s 0 , log t -energy configurations on A have the same distribution of distances, then g ( s ) is differentiable at s 0 .

2.2. Minimality of ω N s , log t on Circles

Let d u be the 2-dimensional Euclidean metric of R 2 . For α > 0 , we denote by
S α 1 : = { x R 2 : d u ( 0 , x ) = α }
the circle centered at 0 of radius α . We let L ( x , y ) be the geodesic distance between the points x and y on S α 1 ; that is, the length of the shorter arc of S α 1 connecting the points x and y .
The minimality of N distinct equally spaced points on S α 1 with the Euclidean metric d u or the geodesic distance L for the certain s , log t -energy problems is stated in Propositions 1–3.
Proposition 1.
Let N 2 , s 0 , t 1 , and 0 < α < π 1 . Then, ω N is a minimal N-point s , log t -energy configuration on S α 1 with the geodesic distance L if and only if ω N is a configuration of N distinct equally spaced points on S α 1 .
Proposition 2.
Let N 2 , 0 < α < ( e π ) 1 , and s , t satisfy s > 0 , t 0 or s = 0 , t > 0 . Then, ω N is a minimal N-point s , log t -energy configuration on S α 1 with the geodesic distance L if and only if ω N is a configuration of N distinct equally spaced points on S α 1 .
Proposition 3.
Let N 2 , s 0 , t 1 , and 0 < α < 1 / 2 . Then, ω N is a minimal N-point s , log t -energy configuration on S α 1 with the Euclidean metric d u if and only if ω N is a configuration of N distinct equally spaced points on S α 1 .
Note that our approach works only for the case diam ( S α 1 ) < 1 and the conditions 0 < α < π 1 in Proposition 1 and 0 < α < 1 / 2 in Proposition 3 are required for diam ( S α 1 ) < 1 . The case diam ( S α 1 ) 1 remains open for further investigation.

3. Proofs of Main Results

We keep all proof of the main results in this section. In our proof, we may sometimes refer to lemmas. In order to avoid any interruption, we keep all lemmas in Section 4.
Proof of Theorem 4.
Let t 0 be fixed, s > 0 , ω N s , log t be a minimal N-point s , log t -energy configuration on A, and let ω N be an N-point best-packing configuration on A. Since diam ( A ) < 1 and the points in ω N s , log t are distinct, there is a constant c > 0 such that
0 < δ ( ω N s , log t ) c < 1
where the constant c only depends on the set A . This implies that
log 1 c t log 1 δ ( ω N s , log t ) t .
Then,
1 δ N ( A ) log 1 c t / s 1 δ ( ω N s , log t ) log 1 c t / s 1 δ ( ω N s , log t ) log 1 δ ( ω N s , log t ) t / s
E log t s ( ω N s , log t ) 1 / s = E log t s ( A , N ) 1 / s E log t s ( ω N ) 1 / s 1 δ N ( A ) E log t ( ω N ) 1 / s .
Since
lim s 1 δ N ( A ) log 1 c t / s = 1 δ N ( A )
and
lim s 1 δ N ( A ) E log t ( ω N ) 1 / s = 1 δ N ( A ) ,
it follows that
lim s E log t s ( A , N ) 1 / s = 1 δ N ( A ) .
Let ω N * be a cluster configuration of ω N s , log t as s . This implies that there is a sequence { s k } k = 1 R such that s k and ω N s k , log t ω N * as k . Arguing as in (7), we have
1 δ ( ω N s k , log t ) log 1 c t / s k E log t s k ( ω N s k , log t ) 1 / s k = E log t s k ( A , N ) 1 / s k E log t s k ( ω N ) 1 / s k
1 δ ( ω N ) E log t ( ω N ) 1 / s k .
Taking k , we obtain
δ N ( A ) = δ ( ω N ) δ ( ω N * ) .
This means that ω N * is also an N-point best-packing configuration on A. □
Proof of Theorem 5.
First of all, we show that g ( s ) is continuous on ( 0 , ) . Let s > 0 and let ω N s , log t be a minimal N-point s , log t -energy configuration on A. Using Lemma 4, we obtain for any ω N s , log t ,
lim inf r s g ( r ) g ( s ) r s lim inf r s E log t r ( ω N s , log t ) E log t s ( ω N s , log t ) r s
lim r s E log t + 1 r ( ω N s , log t ) = E log t + 1 s ( ω N s , log t ) sup G log t + 1 s ( A , N ) > 0 ,
and
lim sup r s g ( r ) g ( s ) r s lim sup r s E log t r ( ω N r , log t ) E log t s ( ω N r , log t ) r s lim sup r s E log t + 1 s ( ω N r , log t ) ,
where the second inequality in (8) follows from the arbitrariness of ω N s , log t and the last inequality in (8) follows from Lemma 3.
Let ω N be a fixed configuration of N distinct points of A. Note that 0 < δ ( ω N ) < 1 . For all r ( s / 2 , s ) , we have
1 δ ( ω N r , log t ) s / 2 log 1 δ ( ω N r , log t ) t 1 δ ( ω N r , log t ) r log 1 δ ( ω N r , log t ) t E log t r ( ω N r , log t )
E log t r ( ω N ) 1 δ ( ω N ) r log 1 δ ( ω N ) t N ( N 1 )
1 δ ( ω N ) s log 1 δ ( ω N ) t N ( N 1 ) .
That is,
( δ ( ω N r , log t ) ) s / 2 log 1 δ ( ω N r , log t ) t ( δ ( ω N ) ) s log 1 δ ( ω N ) t ( N ( N 1 ) ) 1 .
This implies that for all r ( s / 2 , s ) ,
δ ( ω N r , log t ) log 1 δ ( ω N r , log t ) 2 t / s ( δ ( ω N ) ) 2 log 1 δ ( ω N ) 2 t / s ( N ( N 1 ) ) 2 / s = : c 1 > 0 .
Since by Lemma 1,
h ( x ) : = x log 1 x β , β > 0 ,
is a strictly increasing function on ( 0 , 1 ) , there exists a constant c 2 > 0 such that for all r ( s / 2 , s ) ,
δ ( ω N r , log t ) c 2 > 0 .
Therefore, E log t + 1 s ( ω N r , log t ) are bounded above where r ( s / 2 , s ) . From this and (9),
lim sup r s g ( r ) g ( s ) r s lim sup r s E log t + 1 s ( ω N r , log t ) < .
Let s 0 . Using Lemma 4, we also obtain for any ω N s , log t ,
lim sup r s + g ( r ) g ( s ) r s lim sup r s + E log t r ( ω N s , log t ) E log t s ( ω N s , log t ) r s
lim r s + E log t + 1 r ( ω N s , log t ) = E log t + 1 s ( ω N s , log t ) inf G log t + 1 s ( A , N ) < ,
and
lim inf r s + g ( r ) g ( s ) r s lim inf r s + E log t r ( ω N r , log t ) E log t s ( ω N r , log t ) r s lim inf r s + E log t + 1 s ( ω N r , log t ) > 0 ,
where the second inequality in (11) follows from the arbitrariness of ω N s , log t and the last inequality in (12) follows from Lemma 3.
The inequalities (8) and (10)–(12) imply that for all s > 0 ,
0 < lim inf r s g ( r ) g ( s ) r s lim sup r s g ( r ) g ( s ) r s <
and for all s 0 ,
0 < lim inf r s + g ( r ) g ( s ) r s lim sup r s + g ( r ) g ( s ) r s < .
The inequalities in (13) and (14) further imply that g ( s ) is continuous for all s > 0 and is right continuous at s = 0 . □
Proof of Theorem 6.
Let s 0 > 0 . In order to show Theorem 6, it suffices to show that any cluster configuration of ω N s , log t as s s 0 + or as s s 0 is a minimal N-point s 0 , log t -energy configuration on A.
Let ω N * be a cluster configuration of ω N s , log t , as s s 0 + . Then, there is a sequence { s k } k = 1 ( s o , ) such that s k s 0 and ω N s k , log t ω N * as k . Let α = diam ( A ) . For any configuration of N distinct points ω N on A, notice that α s E log t s ( ω N ) is an increasing function of s. Applying the continuity of g ( s ) : = E log t s ( A , N ) at s 0 , we have
α s 0 E log t s 0 ( ω N * ) = lim k α s 0 E log t s 0 ( ω N s k , log t ) lim k α s k E log t s k ( ω N s k , log t )
= lim k α s k E log t s k ( A , N ) = α s 0 E log t s 0 ( A , N ) .
This implies that E log t s 0 ( ω N * ) = E log t s 0 ( A , N ) . Hence, ω N * is a minimal N-point s 0 , log t -energy configuration on A.
Let ω N * * be a cluster configuration of ω N s , log t , as s s 0 . Then, there is a sequence { s k } k = 1 [ 0 , s 0 ) such that s k s 0 and ω N s k , log t ω N * * as k . Without loss of generality, we may assume that s 0 / 2 < s k < s 0 for all k. For any configuration of N distinct points ω N of A, observe that δ ( ω N ) s E log t s ( ω N ) is a decreasing function of s. It follows from the continuity of the function g ( s ) that g ( s ) is bounded above by some number M > 1 for all s ( s 0 / 2 , s 0 ) . For every s 0 / 2 < s k < s 0 ,
( δ ( ω N s k , log t ) ) s 0 / 2 log 1 δ ( ω N s k , log t ) t ( δ ( ω N s k , log t ) ) s k log 1 δ ( ω N s k , log t ) t
E log t s k ( ω N s k , log t ) M .
Then,
δ ( ω N s k , log t ) log 1 δ ( ω N s k , log t ) 2 t / s 0 M 2 / s 0 > 0 .
Using Lemma 1, there is a constant c > 0 such that
δ ( ω N s k , log t ) c > 0 for all k N .
Using the continuity of g ( s ) : = E log t s ( A , N ) at s 0 , we have
( δ ( ω N * * ) ) s 0 E log t s 0 ( ω N * * ) = lim k ( δ ( ω N s k , log t ) ) s 0 E log t s 0 ( ω N s k , log t )
lim k ( δ ( ω N s k , log t ) ) s k E log t s k ( ω N s k , log t ) = lim k ( δ ( ω N s k , log t ) ) s k E log t s k ( A , N )
= ( δ ( ω N * * ) ) s 0 E log t s 0 ( A , N ) .
This implies that E log t s 0 ( ω N * * ) = E log t s 0 ( A , N ) . Hence, ω N * * is a minimal N-point s 0 , log t -energy configuration on A. □
Proof of Theorem 7.
Firstly, we show (5). Let s 0 be fixed and { r k } k = 1 ( s , ) be a sequence such that r k s as k and
lim k E log t + 1 s ( ω N r k , log t ) = lim inf r s + E log t + 1 s ( ω N r , log t ) .
Since A N is compact, there exists a subsequence { s } = 1 { r k } k = 1 such that
lim ω N s , log t = ω N *
and ω N * is a minimal N-point s , log t -energy configuration by Theorem 6. By
lim k E log t + 1 s ( ω N r k , log t ) = lim E log t + 1 s ( ω N s , log t ) ,
(11), (12), (15) and (16), we get
lim inf r s + g ( r ) g ( s ) r s lim inf r s + E log t + 1 s ( ω N r , log t ) = lim E log t + 1 s ( ω N s , log t )
= E log t + 1 s ( ω N * ) inf G log t + 1 s ( A , N ) lim sup r s + g ( r ) g ( s ) r s .
Then,
g + ( s ) = inf G log t + 1 s ( A , N ) .
It is easy to check that from Lemma 3, the constant inf G log t + 1 s ( A , N ) in (18) is finite. This verifies (5).
Then, we prove (6). Let s > 0 be fixed and { r k } k = 1 [ 0 , s ) be a sequence such that r k s as k and
lim k E log t + 1 s ( ω N r k , log t ) = lim sup r s E log t + 1 s ( ω N r , log t ) .
Because A N is compact, there exists a subsequence { s } = 1 { r k } k = 1 such that
lim ω N s , log t = ω N * *
and ω N * * is a minimal N-point s , log t -energy configuration by Theorem 6. Then, we get
lim k E log t + 1 s ( ω N r k , log t ) = lim E log t + 1 s ( ω N s , log t ) .
Using (8), (9), (19) and (20), we obtain
lim inf r s g ( r ) g ( s ) r s sup G log t + 1 s ( A , N ) E log t + 1 s ( ω N * * )
= lim E log t + 1 s ( ω N s , log t ) = lim sup r s E log t + 1 s ( ω N r , log t ) lim sup r t g ( r ) g ( s ) r s .
Then,
g ( s ) = sup G log t + 1 s ( A , N ) .
Then, we want to show that sup G log t + 1 s ( A , N ) is finite. Let ω N be a fixed configuration of N distinct points on A and let ω N s , log t be any minimal N-point s , log t configurations. Then,
( δ ( ω N s , log t ) ) s log 1 δ ( ω N s , log t ) t E log t s ( ω N s , log t )
E log t s ( ω N ) ( δ ( ω N ) ) s log 1 δ ( ω N ) t N ( N 1 ) .
That is,
δ ( ω N s , log t ) log 1 δ ( ω N s , log t ) t / s δ ( ω N ) log 1 δ ( ω N ) t / s ( N ( N 1 ) ) 1 / s = : c 1 > 0 .
It follows from Lemma 1 that there is a constant c 2 > 0 such that for any ω N s , log t ,
δ ( ω N s , log t ) c 2 > 0 .
Since by Lemma 2,
p ( x ) : = 1 x s log 1 x t + 1 ,
is a strictly decreasing function on ( 0 , 1 ) , the set G log t + 1 s ( A , N ) is bounded above. This implies that sup G log t + 1 s ( A , N ) in (21) is finite. Hence, (6) is proved. □
Proof of Theorem 8.
 
(a): This is a direct consequence of Theorem 7.
(b): Let s 0 0 and ω N * be a cluster configuration of { ω N s , log t } as s s 0 + . Then, there exists a sequence { s k } k = 1 ( s 0 , ) such that lim k s k = s 0 and lim k ω N s k , log t = ω N * . Then, ω N * is a minimal N-point s 0 , log t -energy configuration by Theorem 6. Using (5) and the similar argument used to show (12), we have
E log t + 1 s 0 ( ω N * ) = lim k E log t + 1 s 0 ( ω N s k , log t ) lim k g ( s k ) g ( s 0 ) s k s 0 = g + ( s 0 ) = inf G log t + 1 s 0 ( A , N ) .
Since inf G log t + 1 s 0 ( A , N ) E log t + 1 s 0 ( ω N * ) ,
E log t + 1 s 0 ( ω N * ) = inf G log t + 1 s 0 ( A , N ) = g + ( s 0 ) .
(c): Let s 0 > 0 and ω N * * be a cluster configuration of { ω N s , log t } as s s 0 . Then, there exists a sequence { s k } k = 1 [ 0 , s 0 ) such that lim k s k = s 0 and lim k ω N s k , log t = ω N * * . Then, ω N * * is a minimal N-point s 0 , log t -energy configuration by Theorem 6. Using (6) and the similar argument used to show (10), we have
E log t + 1 s 0 ( ω N * * ) = lim k E log t + 1 s 0 ( ω N s k , log t ) lim k g ( s k ) g ( s 0 ) s k s 0 = g ( s 0 ) = sup G log t + 1 s 0 ( A , N ) .
Since E log t + 1 s 0 ( ω N * * ) sup G log t + 1 s 0 ( A , N ) ,
E log t + 1 s 0 ( ω N * * ) = sup G log t + 1 s 0 ( A , N ) = g ( s 0 ) .
(d): This is a direct consequence of (b) and (c). □
Proof of Proposition 1.
Let N 2 , s 0 , t 1 , and 0 < α < π 1 . We prove this proposition using Lemma 5. The function k : ( 0 , 1 ) : R in the lemma is
k ( x ) = 1 x s log 1 x t .
By Lemma 2, k ( x ) is strictly decreasing on ( 0 , 1 ) . Since for all x ( 0 , 1 ) ,
k ( x ) = 1 x s + 2 log 1 x 2 + t ( 1 + t ) t + ( t + 2 s t ) log 1 x + s ( 1 + s ) log 2 1 x > 0 ,
k ( x ) is strictly convex on ( 0 , 1 ) . Hence, because the function k ( x ) satisfies all required properties in Lemma 5, all minimal N-point K-energy configurations on S α 1 are configurations of N distinct equally spaced points on S α 1 with respect to the arc length and vice versa. □
Proof of Proposition 2.
Let N 2 , 0 < α < ( e π ) 1 , and s , t satisfy s > 0 , t 0 or s = 0 , t > 0 . We can use the same lines of reasoning as in the proof of Proposition 1 except the function k is considered on ( 0 , 1 / e ) and for all x ( 0 , 1 / e ) ,
k ( x ) = 1 x s + 2 log 1 x 2 + t ( 1 + t ) t + ( t + 2 s t ) log 1 x + s ( 1 + s ) log 2 1 x
1 x s + 2 log 1 x 2 + t t 2 + 2 s t log 1 x + s ( 1 + s ) log 2 1 x + log 1 x 1 t > 0 .
Hence, because the function k ( x ) satisfies all required properties in Lemma 5, Proposition 2 is proved. □
Proof of Proposition 3.
Let N 2 , s 0 , t 1 , and 0 < α < 1 / 2 . Again, we want to use Lemma 5. The function k : ( 0 , π α ] R in the lemma is
k ( x ) = 1 2 α sin ( x / 2 α ) s log 1 2 α sin ( x / 2 α ) t .
Since 2 α sin ( x / 2 α ) is strictly increasing on ( 0 , π α ] and ( 1 / x s ) ( log ( 1 / x ) ) t is strictly decreasing on ( 0 , 1 ) , k ( x ) is strictly decreasing on ( 0 , π α ] . Then, we want to show that k ( x ) is strictly convex on ( 0 , π α ] , meaning that
k ( x ) > 0 for all x ( 0 , π α ) .
To show (23), it suffices to show that q ( x ) > 0 for all x ( 0 , π / 2 ) , where
q ( x ) : = 1 2 α sin x s log 1 2 α sin x t .
Because for all x ( 0 , π / 2 ) ,
q ( x ) = s ( cot 2 x ) ( 2 α sin x ) s log 1 2 α sin x t 1
+ ( t 1 ) ( cot 2 x ) ( 2 α sin x ) s log 1 2 α sin x t 2 s log 1 2 α sin x + t
+ ( csc 2 x + s cot 2 x ) ( 2 α sin x ) s log 1 2 α sin x t 1 s log 1 2 α sin x + t > 0 ,
k ( x ) is strictly convex on ( 0 , π α ] . Hence, the function k ( x ) satisfies all required properties in Lemma 5. This completes the proof. □

4. Appendix: Auxiliary Lemmas

Lemmas 1–3 are very fundamental but highly important. For example, making use of Lemma 3 and the assumption that diam ( A ) < 1 , we can conclude that
E log s ( A , N ) N ( N 1 ) ( diam ( A ) ) s log 1 diam ( A ) t > 0 .
Lemma 1.
Let β 0 and h : ( 0 , 1 ) ( 0 , ) be a function defined by
h ( x ) : = x log 1 x β for all x ( 0 , 1 ) .
Then, h ( x ) is strictly increasing on ( 0 , 1 ) .
Proof of Lemma 1.
 
Because
h ( x ) = β log 1 x ( β + 1 ) + log 1 x β
and ( log ( 1 / x ) ) β > 0 for all x ( 0 , 1 ) and β 0 , h ( x ) > 0 for all x ( 0 , 1 ) . Therefore, h ( x ) is strictly increasing on ( 0 , 1 ) . □
Lemma 2.
Let ( s , t ) [ 0 , ) × [ 0 , ) ] { ( 0 , 0 ) } and p : ( 0 , 1 ) ( 0 , ) be a function defined by
p ( x ) : = 1 x s log 1 x t for all x ( 0 , 1 ) .
Then, p ( x ) is strictly decreasing on ( 0 , 1 ) .
Proof of Lemma 2.
Using Lemma 1, we set β = t / s and
p ( x ) = 1 h ( x ) s = 1 x s log 1 x t
is strictly decreasing on ( 0 , 1 ) . □
Lemma 3.
Let ( A , d ) be an infinite compact metric space with diam ( A ) < 1 and s , t 0 . Then, for all N-point configurations ω N A ,
E log t s ( ω N ) N ( N 1 ) ( diam ( A ) ) s log 1 diam ( A ) t .
Proof of Lemma 3.
The proof relies on the fact that p ( x ) in Lemma 2 is strictly decreasing on ( 0 , 1 ) . □
The following is the main lemma of this paper. It allows us to prove analogues of Theorems 1–3.
Lemma 4.
Let ( A , d ) be an infinite compact metric space with diam ( A ) < 1 and ω N = { x 1 , , x N } be any configuration of N distinct points of A. Then, for any s > r 0 and t 0 ,
E log t + 1 r ( ω N ) E log t s ( ω N ) E log t r ( ω N ) s r E log t + 1 s ( ω N ) .
Proof of Lemma 4.
Let x i , x j ω N where 1 i j N , let s > r 0 , and let t 0 . Then,
1 d ( x i , x j ) r log 1 d ( x i , x j ) 1 d ( x i , x j ) s 1 d ( x i , x j ) r s r 1 d ( x i , x j ) s log 1 d ( x i , x j ) .
Since log 1 d ( x i , x j ) t > 0 ,
1 d ( x i , x j ) r log 1 d ( x i , x j ) t + 1 1 d ( x i , x j ) s log 1 d ( x i , x j ) t 1 d ( x i , x j ) r log 1 d ( x i , x j ) t s r
1 d ( x i , x j ) s log 1 d ( x i , x j ) t + 1 .
It follows that
E log t + 1 r ( ω N ) E log t s ( ω N ) E log t r ( ω N ) s r E log t + 1 s ( ω N ) .
Let Γ be a rectifiable simple closed curve in R m , m 2 , of length | Γ | with a chosen orientation. We recall that L ( x , y ) is the geodesic distance between the points x and y on Γ . With the help of the following lemma [2] (Theorem 2.3.1), we can prove propositions 1–3.
Lemma 5.
Let k : ( 0 , | Γ | / 2 ] R be a strictly convex and decreasing function defined at u = 0 by the (possibly infinite) value lim u 0 + k ( u ) and let K be the kernel on Γ × Γ defined by K ( x , y ) = k ( L ( x , y ) ) . Then, all minimal N-point K-energy configurations on Γ are configurations of N distinct equally spaced points on Γ with respect to the arc length and vice versa.

5. Discussion and Conclusions

We introduce minimal N-point s , log t -energy constants and configurations of an infinite compact metric space ( A , d ) . Such constants and configurations are generated using the kernel (or potential):
K log t s ( x , y ) = 1 d ( x , y ) s log 1 d ( x , y ) t , s 0 , t 0 .
In this paper, we studied the asymptotic properties of minimal N-point s , log t -energy constants and configurations of A with diam ( A ) < 1 where N , t are fixed and s is varying. In Theorem 4, we show that
lim s E log t s ( A , N ) 1 / s = 1 δ N ( A ) .
and minimal N-point s , log t -energy configurations on A tend to an N-point best-packing configuration on A as s . Then, we show that the s , log t -energy
g ( s ) : = E log t s ( A , N )
is continuous and right differentiable on [ 0 , ) and is left differentiable on ( 0 , ) in Theorems 5 and 7. Using the continuity of E log t s ( A , N ) in the variable s, we show in Theorem 6 that for any s 0 > 0 , any cluster configuration of ω N s , log t , as s s 0 , is a minimal N-point s 0 , log t -energy configuration on A.
We want to emphasize that when t = 0 , our proof of Theorems 4–8 can handle the case diam ( A ) 1 . However, when t > 0 , we require that diam ( A ) < 1 . This is because our methods rely on the positivity of the kernel K log t s ( x , y ) and the property that K log t s ( x , y ) decreases as d ( x , y ) increases. These limitations would leave room for future improvement (when t > 0 and diam ( A ) 1 ).
Note that the kernel K log t s ( x , y ) is symmetric, namely K log t s ( x , y ) = K log t s ( y , x ) . When the metric space ( A , d ) has a great symmetry, we observe that such minimal N-point s , log t -energy configurations should be evenly distributed over the set A . The most prominent sets with a great symmetry are the spheres:
S m : = { x R m + 1 : d u ( x , 0 ) = 1 } ,
where d u is the m + 1 -dimensional Euclidean metric. As a motivated result, it is known that for s > 0 , minimal N-point s-energy configurations and minimal N-point log-energy configurations on the metric space ( S m , d u ) are asymptotically uniformly distributed with respect to the surface area measure on S m as N (see [2] (Theorem 6.1.7)). We refer the reader to the review article [22] for a number of applications of uniformly distributed points on the sphere S m . Our investigation in this paper on circles in R 2 serves as a basis example of our observation. In Propositions 1–3, we prove that for certain values of s and t , all minimal N-point s , log t -energy configurations on the circle S α 1 with diam ( S α 1 ) < 1 are the configurations of N distinct equally spaced points S α 1 .
In addition to the problem on the sphere S m , explaining the limiting distribution as the N of minimal N-point s , log t -energy configurations on a compact set in a finite dimensional Euclidean space would be another interesting problem. We refer the reader to Chapters 4 and 8 in [2] or [9,10] for the study of such a problem for the minimal N-point s-energy and log-energy configurations. The study of such limiting distribution problem is important in both theoretical and computational sciences. For example, it shows applications in computer-aided geometric design, finite element tesselations, and statistical sampling.

Author Contributions

Conceptualization, N.B.; formal analysis, N.L. and N.B.; writing—original draft, N.L. and N.B.; supervision, N.B.; writing—review and editing, N.B. Both authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thomson, J.J. On the Structure of the Atom: An Investigation of the Stability and Periods of Oscillation of a number of Corpuscles arranged at equal intervals around the Circumference of a Circle; with Application of the results to the Theory of Atomic Structure. Philos. Mag. 1904, 7, 237–265. [Google Scholar] [CrossRef] [Green Version]
  2. Borodachov, S.V.; Hardin, D.P.; Saff, E.B. Discrete Energy on Rectifiable Sets; Springer Monographs in Mathematics; Springer: New York, NY, USA, 2019. [Google Scholar]
  3. Wales, D.J.; Ulker, S. Structure and dynamics of spherical crystals characterized for the Thomson problem. Phys. Lett. B 2006, 74, 212101. [Google Scholar] [CrossRef] [Green Version]
  4. Föppl, L. Stabile Anordnungen von Elektronen im Atom. J. Reine Angew. Math. 1912, 141, 251–301. [Google Scholar]
  5. Yudin, V.A. The minimum of potential energy of a system of point charges. Discret. Math. Appl. 1992, 4, 112–115. Discret. Math. Appl.1993, 3, 75–81. (In Russian) [Google Scholar] [CrossRef]
  6. Andreev, N.N. An extremal property of the icosahedron. East J. Approx. 1996, 2, 459–462. [Google Scholar]
  7. Schwartz, R.E. Five Point Energy Minimization: A Synopsis. Constr. Approx. 2020, 51, 537–564. [Google Scholar] [CrossRef]
  8. Schwartz, R.E. The phase transition in five point energy minimization, research monograph. arXiv 2016, arXiv:1610.03303. [Google Scholar]
  9. Landkof, N.S. Foundations of Modern Potential Theory; Springer: Berlin, Germany, 1972. [Google Scholar]
  10. Hardin, D.P.; Saff, E.B. Minimal Riesz energy point configurations for rectifiable d-dimensional manifolds. Adv. Math. 2005, 193, 174–204. [Google Scholar] [CrossRef]
  11. Mhaskar, H.N.; Saff, E.B. Where does the sup norm of a weighted polynomial live? Constr. Approx. 1985, 1, 71–91. [Google Scholar] [CrossRef]
  12. Gonchar, A.A.; Rakhmanov, E.A. Equilibrium distributions and the degree of rational approximation of analytic functions. Math. USSR-Sb. 1989, 62, 305–348. [Google Scholar] [CrossRef]
  13. Lubinsky, D.S.; Mhaskar, H.N.; Saff, E.B. Freud’s conjecture for exponential weights. Bull. Amer. Math. Soc. 1986, 15, 217–221. [Google Scholar] [CrossRef] [Green Version]
  14. Totik, V. Weighted polynomial approximation for convex external fields. Constr. Approx. 2000, 16, 261–281. [Google Scholar] [CrossRef]
  15. Saff, E.B.; Totik, V. Logarithmic Potentials with External Fields; Springer: New York, NY, USA, 1997. [Google Scholar]
  16. Trefethen, L.N. Approximation Theory and Approximation Practice; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2013. [Google Scholar]
  17. Smale, S. Mathematical Problems for the Next Century; Mathematics: Frontiers and Perspectives; American Mathematical Society: Providence, RI, USA, 2000. [Google Scholar]
  18. Shub, M.; Smale, S. Complexity of Bezout’s theorem. III. Condition number and packing. J. Complex. 1993, 9, 4–14. [Google Scholar] [CrossRef] [Green Version]
  19. Brauchart, J.S.; Hardin, D.P.; Edward, B.S. The next-order term for optimal Riesz and logarithmic energy asymptotics on the sphere. Contemp. Math 2012, 578, 31–61. [Google Scholar]
  20. Whyte, L.L. Unique arrangements of points on a sphere. Am. Math. Mon. 1952, 59, 606–611. [Google Scholar] [CrossRef]
  21. Dragnev, P.D.; Legg, D.A.; Townsend, D.W. Discrete logarithmic energy on the sphere. Pac. J. Appl. Math. 2002, 207, 345–358. [Google Scholar] [CrossRef]
  22. Brauchart, J.S.; Grabner, P.J. Distributing many points on spheres: Minimal energy and designs. J. Complexity 2015, 31, 293–326. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Loesatapornpipit, N.; Bosuwan, N. Asymptotic Properties of Discrete Minimal s,logt-Energy Constants and Configurations. Symmetry 2021, 13, 932. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13060932

AMA Style

Loesatapornpipit N, Bosuwan N. Asymptotic Properties of Discrete Minimal s,logt-Energy Constants and Configurations. Symmetry. 2021; 13(6):932. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13060932

Chicago/Turabian Style

Loesatapornpipit, Nichakan, and Nattapong Bosuwan. 2021. "Asymptotic Properties of Discrete Minimal s,logt-Energy Constants and Configurations" Symmetry 13, no. 6: 932. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13060932

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop