An iterative algorithm to find a closed form of solution for Hammerstein nonlinear integral equation constructed by the concept of cosm-rs

In this article, we study Hammerstein nonlinear integral equation introduced by the concept of conductor-like screening model for real solvent. We use a combination of modified homotopy perturbation and Adomian decomposition method to solve the above nonlinear integral equation. The approach is based on preparing a closed form of solution to recognize the Hammerstein integral equation for the determination of the chemical potential of a surface segment as a function of screening charge density. This solution can be more effective in analyzing the future concept on the above chemical problem and decrease computational operations. Finally, we prove convergence of the proposed method and compare our results to some other works to show validity and high accuracy.


Introduction
To introduce of the chemical concept of cosm-rs, we used from [2,6,7,16] references, which the chemical potential of a surface segment with SCD in an ensemble described by normalized distribution function p s ( ) is given by where R is the gas constant, s ( ) is a measure for the affinity of the system S to a surface of polarity at temperature T, and In this paper, we focus on (1) and convert it to a suitable form of integral equation in the following steps, assuming f ( ) = exp − s ( ) RT and Ω( , � ) = exp ( −a eff e( , � ) RT ) , and we conclude that, let us getting , � ∈ [a, b] ⊂ ℜ , k( , � ) = p s ( � )Ω( , � ) and ( � , f ( � )) = (f ( � )) −1 , and then we can give a Hammerstein nonlinear integral equation as follows: Existence of solution of Eq. (4) and some of nonlinear integral equations, you can see [3,4,12]. Also in [9,13,14] introduced some numerical methods to solve Hammerstein nonlinear integral equation. But we decide to solve the above (2)

3
problem by an effective iterative algorithm and comparing it with some other works.

Introducing an iterative algorithm
Homotopy perturbation method is a semianalytic method to solve nonlinear problems which introduced by [5,8] that it is an important concept of topology and perturbations theory. Some modifications and improvement of the above method can be seen in [10,11]. Now, we consider the general operator form of nonlinear problem, where A is a general differential operator, B is a boundary operator, h is a known analytic function, and Γ is the boundary of the domain Ω . We divide the general operator A to two nonlinear operators as N 1 and N 2 . Also, the function h is converted to two functions such as h 1 and h 2 , and thus we have Similar to [10,11], the modified homotopy perturbation can be introduced as, where is an approximation of f and variation p in interval [0, 1] also p = 0 to p = 1 concludes that N 1 ( ) = h 1 ( ) to A( ) = 0 . In fact, we can find a solution of the (5) for p = 1 . In the general case, we introduce a form of solution by series, Similar to (5), we consider the nonlinear integral equation (4) to this form, by using modified homotopy perturbation (7) and Eq. (10), we can get N 1 , N 2 operators in the following form, because function h( ) is equal to zero and then h 2 ( ) = −h 1 ( ) , and replacing (11) operators into (7), we can write, substituting (8) into (12) leads to, To get rid of the nonlinearity of the function , we use a linear combination of Adomian polynomials in this form, Replacing (14) into (13) and sorting out to the power of p, we obtain an iterative algorithm as follows:

Algorithm
where Adomian polynomials are given by (see [15]), To find solution of (10), suitable choice of the 0 (s) is important and since h( ) = 0 so, we set 0 (s) = 0.

Convergence of the proposed method
According to the above explanations, since existence and uniqueness of solution of Eq. (4) is verified, we prove convergence of the series created by algorithm (15) in the proposed method.
Theorem 1 Series (8), which terms of this series are created by algorithm (15), converges to solution of Hammerstein nonlinear integral equation (4) in Banach space (C[a, b], ||f || ∞ ) , if the following conditions are held,

a t i s f y i n L i p s c h i t z c o nd i t i o n r e s p e c t t o s e c o n d c o m p o n e n t ,
Proof According to series (8), we define sequence as, then, with the help of jth step of algorithm (15), we can get Now, we show that the above sequence is a Cuachy sequence in Banach space C [a, b]. To this end ∀n, m ∈ ℕ;n > m and using (i) we have, using (14) and (ii)-(iii),  (4), on the other hand, since (4) has unique solution, then v * ( ) = f ( ). ◻

Application of proposed method and comparison with some other works
In this section, we solve two examples to show the applications of the above method.
Example 1 Consider the following cosmo-rs integral equation for a particular case of the energy expression, namely the electrostatic misfit energy [9], where and So, we consider cosmo-rs integral equation as Hammerstein integral equation as follows: assuming we can rewrite the Hammerstein integral equation in the following form, According to the "Introducing an iterative algorithm" section , we can define the operators N 1 and N 2 as, (24) N 1 (f ( )) = f ( ), because the function h( ) = 0 in problem (20), by choosing h 1 ( ) therefore h 2 ( ) = 0 − h 1 ( ) and also by replacing (26) and (8) with equalling the coefficients of the power p to zero, we get the following iterative algorithm, Algorithm 1 An appropriate choice for 0 ( ) can be 0 ( ) = 1 , thus A 0 ( � ) = 1 0 ( � ) = 1 . To calculate the second step of the algorithm, we do, as follows, by simplifying exp ( + � ) 2 , and we apply the Maclaurin expansion for exp −2 � , thus the factor exp − 2 can come out from integration, We solve the rest of this example by Mathematica 11.2, where for computing of the above integrals, and we used cubic spline and normal series to approximate of and (sin(5 � + 2.5)) 2 respectively. Thus we can give an approximation of solution of the electrostatic misfit energy (20-23) by two first terms of series (8)  (32) f (t) = g(t) + ∫ For solving the example, we consider the following equation According to the definition of modified homotopy perturbation (7), we introduce N 1 and N 2 operators and h(t) to this form, by substituting (33) and (8) into (7) and using Adomian polynomials, we can write where Adomian polynomials are, by sorting and putting the power of p equal to zero, we obtain the following algorithm, Algorithm 2 ,  Table 1, Case 2 we consider to h(t) as the initial guess of a solution, and so in Algorithm 2 for two steps, we obtain 0 (t) = 1, (35)  Therefore, in this case, we have provided exact solution in a closed form as follows, Although [1] and [9] have high accuracy, the closed form of the solution is not provided. Also we compare our result with [1,9] and [14] in Table 2.

Conclusion
In this paper, we solved the famous chemistry problem that led to the Hammerstein nonlinear integral equation by a modified homotopy perturbation method. We obtained a closed form of solution by converting the nonlinear problem into more simple linear problems and we produced an iterative algorithm with a highly accurate solution. Moreover, we presented another example to compare our method with some other techniques. As it can be seen in Table 2, our method led to the exact solution in this case. This method can be useful to solve other nonlinear integral equations especially those involved chemical and physical problems as well.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. (36) Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Bernstein method [14] Sinc method [9]