# Mathematical Modelling for Joining Boron Nitride Graphene with Other BN Nanostructures

Authors: Nawa A. Alshammari
Journal: Advances in Mathematical Physics
 (2020)

---

## Abstract

Boron nitride (BN) nanomaterials such as boron nitride graphenes, boron nitride nanotubes, and boron nitride nanocones are attracting attention among the most promising nanomaterials due to their physical, chemical, and electronic properties when compared to other nanomaterials. BN nanomaterials suggest many exciting potential applications in various fields. Joining between BN nanostructures gives new enhanced structures with outstanding properties and potential applications for design of probes for scanning tunnelling microscopy and other nanoscale devices. This paper uses calculus of variations to model the joining between BN graphene with other BN nanostructures: BNNTs and BNNCs. Furthermore, during the joining between these BN nanostructures, this research examines two models which are depending on the curvature of the join profile. For the first case, Model I refers to when the join profile only includes positive curvature where for the second case, Model II is considered for both positive and negative curvatures. Thus, the purpose of this research is to formulate the basic underlying structure to present simple models based on joining BN graphene to other BN nanostructures.


---

## Body

## 1. Introduction

Materials in nanoscale have attracted research in different scientific fields because of their chemical, physical, and electronic properties. In particular, carbon in nanoscale has different shapes and structures such as carbon nanotubes, carbon nanocones, carbon nanosheets (graphene), and fullerene.Due to similarity to carbon nanostructures, boron nitride nanostructures have received specific attention due to their remarkable properties that lead them to a wide range of applications such as developing many nanoscale devices [1, 2]. Hexagonal BN structure involves boron and nitrogen atoms bound by strong covalent bonds. Moreover, many forms of boron nitride in nanoscale have been reported in the previous literature such as boron nitride nanotubes (BNNTs), boron nitride graphene (BN graphene), boron nitride nanocones (BNNCs), and boron nitride fullerene (BN fullerene) [3].BNNTs were first discovered in 1994 [4] and contain a tubular structure as carbon nanotubes with replacing carbon atoms by boron and nitrogen atoms arranging in a hexagonal lattice. BNNTs have gained significant attention in recent years because of their interesting properties that are not available in other nanomaterials. While they are structurally close to carbon nanotubes, BNNTs have completely different physical properties. For example, BNNTs have mechanical, thermal, electrical, and chemical properties like high tensile stiffness and high thermal conductivity. Also, their electronic properties are independent of their chirality and radii of the tube [4–6].Another boron nitride nanostructure has attracted much attention, that is, boron nitride graphene. Geometrically, it is similar to carbon nanosheet with boron and nitrogen atoms instead of carbon atoms. BN graphene has physical and chemical distinct properties, such as high-temperature stability, intrinsic electrical insulation, and semiconducting and antioxidation ability. As a result, BN graphene has a promising application in various areas; namely, it is used as a two-dimensional filler, nonwetting coatings, and field emitters [7–9].Boron nitride nanocones (BNNCs) are important structures of BN nanostructures. These conical structures have been considered in many studies experimentally and theoretically. Their special structure offers suitability for fabrication of nanodevices such as electron field emitters, sensors, nanoindenters, and probe tips. Carbon nanocones (CNCs) can be formed by curling the fan-shaped graphite sheet with a disclination angle of an integer multiple of 60° where BNNCs can be produced by only rolling up the white graphene with a disclination angle of an integer multiple of 120° [10]. In addition, carbon nanocones have five possible structures depending on the number of pentagons, and the conical structure of boron nitride has been studied using transmission electron microscopy and observed their angles to be 84°, 19.2°, and 38.9° [11].New structures produced from the joining between two nanostructures can enhance the physiochemical and electrochemical performances of the joined materials. For example, carbon nanotubes and carbon nanosheets are useful in nanosensors and nanooscillators, and connections of these carbon nanostructures may also be employed in the same applications. Furthermore, the same applications exist for joining BNNTs and BN graphene [1]. Moreover, the composite structures are useful for the design of probes for scanning tunnelling microscopy, energy storage, and other electronic devices as carriers for drug delivery [2].Research in [12] minimizes the elastic energy (which depends on the axial curvature only) using calculus of variations to determine the joining area between two carbon nanostructures. The chemical issues are ignored in this model such as the positions of the carbon atoms and bonds. Based on the length constraint of the joining between two nanostructures, the curvature of the join curve can be either positive or both positive and negative. In addition, the study in [13] studies different cases for joining between carbon nanostructures using the same method described above, including the joining of two carbon nanotubes, tube with fullerene, tube with cone, and two fullerenes. Furthermore, this model has been extended to study the joining of some boron nitride nanostructures as in [14].In this paper, we use the same model to determine the shape of the joining between BN graphene with BNNTs and BNNCs, where the arc length of the curve and the size of the defect in the BN graphene are specified, and the distance of the BNNTs and BNNCs from the BN graphene is not prescribed and is determined as a part of the solution.Finally, we comment that using Willmore energy which depends on the axial and rotational curvatures to determine the joining between two nanostructures gives rise to similar join profiles of using elastic energy as studied in [15]. Moreover, the mathematical techniques in [16–19] have similar ideas as in this research.In the following section, we state the fundamental equations of the calculus of variations to model the joining region between BN graphene with BNNTs and BNNCs. In Subsection2.1, the curvature is assumed to remain positive in the join region, and we denote this case as Model I. Subsection 2.2 presents Model II which is assumed to have two join regions, one of positive curvature and the second of negative curvature. Results and discussion are given in Section 4. Section 4 provides the summary.


## 2. Model

Here, the basic variational equations of the model that is used to join boron nitride graphene to other boron nitride nanostructures are formulated. In details, variational calculus is used to identify the curve adopted by a line smoothly connecting a BN graphene base to a vertical other boron nitride nanostructure, where the arc length of the joining curve and the defect site at the BN graphene base are specified. Therefore, the distance in they-direction y0 of the joining to the tube or cone is not specified, and it is found as a part of the solution.We position the first nanostructure (BN graphene) in thex,z-plane with a circular defect of radius x0 centred on the origin. Assuming that a second nanostructure (BNNTs or BNNCs) of radius a is located with its axis colinear with the y-axis starting from an unknown positive distance above the x,z-plane that is denoted by y0, the defect and second nanostructure are assumed to be rotationally symmetric about the y-axis; the problem can be reduced to two-dimensional x,y-plane. The total prescribed arc length ℓ is assumed to connect the defects at the points x0,0 and a,y0. We comment that this technique is used first by Cox and Hill in [12] and then by Baowan et al. [2, 13].Using calculus of variations to find the curveyx, with an element of arc length ds, minimizes the energy functional Jy that is given by
(1)Jy=∫0ℓκ2ds+λ∫0ℓds,where κ is the curvature, λ is a Lagrange multiplier corresponding to the fixed length constraint, and ℓ is the length of the join curve. The boundaries of the join region are x0 and x1, where at x=x0, we have s=0, and at x=x1, we have s=ℓ. The curve in two dimensions described as a graph y=yx, we have κ=y¨/1+ẏ23/2, and ds=1+ẏ21/2dx, then equation (1) becomes
(2)Jy=∫ax0y¨21+ẏ25/2dx+λ∫ax01+ẏ21/2dx,where in this paper, dot denotes differentiation with respect to x. We have two different models based on the curvature’s sign, and the join curvature remains positive as shown in Figure 1 and both positive and negative curvatures as in Figure 2. For these models, we impose the continuity boundary conditions at the first nanostructure (BN graphene) where θ0=0,
(3)yx0=0,ẏx0=0.Figure 1
Model I: positive curvature [2].Figure 2
Model II: positive and negative curvatures [2].To determine the boundary conditions at the second nanostructure (BNNTs or BNNCs), we integrate by parts to derive the standard equation(4)δJy=Fẏ−ddxFy¨δy+Fy¨δẏx0x1+∫x0x1Fy−ddxFẏ+d2dx2Fy¨δydx,where subscripts denote partial derivatives, and here, the function F is given by
(5)Fẏ,y¨=y¨21+ẏ25/2+λ1+ẏ21/2.In this case, onlyẏ is prescribed at x=a, as y0 is unknown, and at x=a, we require the natural or alternative boundary condition given by
(6)Fẏ−ddxFy¨x=a=0.For the case of BNNTs,θ1=π/2. The value of ẏ in Model I ranges from 0 at x=x0 to −∞ at x=a. Thus, the boundary condition for this model is ẏa=−∞. For Model II, ẏ ranges from 0 to −∞, where it changes sign and then ranges from ∞ down to some finite positive value before turning to ∞. Therefore, the boundary condition in the case of Model II is ẏa=∞. For both models, we have at x=a, equation (6) along with the boundary condition for Model I:
(7)ẏa=−∞,and for Model II is
(8)ẏa=∞.Noting that, for the case of BNNCs,θ1 is the angle of the nanocone. From equation (4), we have the usual Euler-Lagrange equation for Fx,y,ẏ,y¨, which is given by
(9)Fy−ddxFẏ+d2dx2Fy¨=0.However, since our functionfẏ,y¨ is independent of y, integrating equation (9) gives
(10)Fẏ−ddxFy¨=C1,where C1 is an arbitrary constant. From the natural boundary condition in equation (6), we find that C1=0, so we obtain
(11)Fẏ=ddxFy¨.Using the definition of the full derivative, we have(12)dFdx=Fx+ẏFy+y¨Fẏ+y⃛Fy¨.SinceFx=Fy=0, from equation (11), we obtain
(13)ddxF−y¨Fy¨=0,by integrating with respect to x, we get
(14)F−y¨Fy¨=−α,where α is an arbitrary constant. Now, we substitute equation (5) into equation (14) and obtain
(15)y¨21+ẏ23=λ+α1+ẏ21/2.As a result, we can write the curvatureκ as
(16)κ=±λ+α1+ẏ21/21/2.


### 2.1. Model I: Positive Curvature

The curvature in Model I is positive through the arc lengthℓ as shown in Figures 3 and 4. As a result, we only consider the positive case from equation (16). Using ẏ=tanθ, equation (8) becomes
(17)κ=λ+αcosθ1/2.Figure 3
Joining of BN graphene and BN nanotube.Figure 4
Joining of BN graphene and BN nanocone.Ifκ=y¨/1+ẏ23/2, and using the same substitution for ẏ, we obtain
(18)dydθ=sinθλ+αcosθ1/2.By introducing the constantk=λ+α/2α1/2, and a new parametric variable ϕ which is defined by
(19)cosθ=1−2k2sin2ϕ,then, equation (18) can be written as
(20)dθdϕ=2kcosθ1−k2sin2ϕ1/2,and on introducing another constant β=2/α1/2, we obtain
(21)dydϕ=2βksinϕ.By the integration and the boundary condition at the pointx0,y0, we obtain
(22)yϕ=2βk1−cosϕ.We note that equation (22) can be considered as a parametric equation for y in terms of the parameter ϕ and denote the value of ϕ at the point where θ=−π/2, with ϕ0=sin−11/2k, and for Model I, −ϕ0<ϕ≤0.Similarly, we derive(23)dxdθ=cosθλ+αcosθ1/2,and in terms of the variable ϕ, we have
(24)dxdϕ=β1−2k2sin2ϕ1−k2sin2ϕ1/2=β21−k2sin2ϕ1/2−1−k2sin2ϕ−1/2.Integrating (24) gives
(25)xϕ=x0+β2Eϕ,k−Fϕ,k,where Fϕ,k and Eϕ,k denote the usual Legendre incomplete elliptic integrals of the first and second kinds, respectively. Using equations (22) and (25), we obtain
(26)x0−a=β2Eϕ0,k−Fϕ0,k,y0=2βk1−cosϕ0.From the definition of the arc length, we have(27)ℓ=∫ax01+ẏ21/2dx.Upon substitutingẏ=tanθ, changing the parameter to ϕ as in cosθ=1−2k2sin2ϕ and integrating, we have
(28)ℓ=βFϕ0,k.Now, we define a dimensionless parameterμ=x0−a/ℓ, which can be shown to be
(29)μ=2Eϕ0,kFϕ0,k−1,where μ=x0−a/ℓ and −1<μ<1. For prescribed values of a, x0, and ℓ, equation (29) can be solved numerically to determine the value of k. Then, substituting k into equation (28), the value of β can be determined, and therefore, y0 can be obtained from equation (28).


### 2.2. Model II: Positive and Negative Curvatures

For the first region, the curvature is positive atx0,0 up until the point xc,yc, where the curvature changes sign to be negative until the point a,y0 is reached as in Figures 5 and 6. At xc,yc, κ=0, and by solving equation (17), we have θc=−cos−1−γ/α, noting that from geometrical considerations, we have −π<θc<−ϕ/2. By making the substitution used in equation (19) for ϕ, we have ϕc=−π/2. By substituting ϕc into equations (22) and (25), we can calculate xc and yc,
(30)xc=x0−β2Ek−Kk,yc=2βk,noting that Kk and Ek are the complete elliptic integrals of the first and second kinds, respectively. In the second joining region which is between xc,yc and a,y0, we take the negative sign of equation (16), and following the same method used in Model I, we obtain
(31)xϕ=x0−β22Ek−E−ϕ,k−2Kk−F−ϕ,k,(32)yϕ=2βk1+cosϕ.Figure 5
Joining of BN graphene and BN nanotube.Figure 6
Joining of BN graphene and BN nanocone.As a result, from the boundary conditions at the pointa,y0, we know that ϕ=−ϕ0. By substitution into (32), we find
(33)x0−a=β22Ek−Eϕ0,k−2Kk−Fϕ0,k,(34)y0=2βk1+cosϕ0.In this model, the arc length is determined in two parts, which is given by(35)ℓ=∫θc01+ẏ21/2dx+∫θc−ϕ/21+ẏ21/2dx=β2Kk−Fϕ0,k.As a result, we find a dimensionless parameterμ=x0−a/ℓ for this case to be
(36)μ=22Ek−Eϕ0,k2Kk−Fϕ0,k−1.Also, for prescribed values ofx0,x1, and ℓ, we can solve equation (36) numerically to find the value of k, and by substituting k into equation (35), we can determine the value of β so that y0 can be determined from equation (34) [12]. Here, we note that equation (29) coincides with (36) for the value k=1/2, and the value of μ at this point is
(37)μ0=2E1/2K1/2−1=0.4569465810⋯.


## 2.1. Model I: Positive Curvature

The curvature in Model I is positive through the arc lengthℓ as shown in Figures 3 and 4. As a result, we only consider the positive case from equation (16). Using ẏ=tanθ, equation (8) becomes
(17)κ=λ+αcosθ1/2.Figure 3
Joining of BN graphene and BN nanotube.Figure 4
Joining of BN graphene and BN nanocone.Ifκ=y¨/1+ẏ23/2, and using the same substitution for ẏ, we obtain
(18)dydθ=sinθλ+αcosθ1/2.By introducing the constantk=λ+α/2α1/2, and a new parametric variable ϕ which is defined by
(19)cosθ=1−2k2sin2ϕ,then, equation (18) can be written as
(20)dθdϕ=2kcosθ1−k2sin2ϕ1/2,and on introducing another constant β=2/α1/2, we obtain
(21)dydϕ=2βksinϕ.By the integration and the boundary condition at the pointx0,y0, we obtain
(22)yϕ=2βk1−cosϕ.We note that equation (22) can be considered as a parametric equation for y in terms of the parameter ϕ and denote the value of ϕ at the point where θ=−π/2, with ϕ0=sin−11/2k, and for Model I, −ϕ0<ϕ≤0.Similarly, we derive(23)dxdθ=cosθλ+αcosθ1/2,and in terms of the variable ϕ, we have
(24)dxdϕ=β1−2k2sin2ϕ1−k2sin2ϕ1/2=β21−k2sin2ϕ1/2−1−k2sin2ϕ−1/2.Integrating (24) gives
(25)xϕ=x0+β2Eϕ,k−Fϕ,k,where Fϕ,k and Eϕ,k denote the usual Legendre incomplete elliptic integrals of the first and second kinds, respectively. Using equations (22) and (25), we obtain
(26)x0−a=β2Eϕ0,k−Fϕ0,k,y0=2βk1−cosϕ0.From the definition of the arc length, we have(27)ℓ=∫ax01+ẏ21/2dx.Upon substitutingẏ=tanθ, changing the parameter to ϕ as in cosθ=1−2k2sin2ϕ and integrating, we have
(28)ℓ=βFϕ0,k.Now, we define a dimensionless parameterμ=x0−a/ℓ, which can be shown to be
(29)μ=2Eϕ0,kFϕ0,k−1,where μ=x0−a/ℓ and −1<μ<1. For prescribed values of a, x0, and ℓ, equation (29) can be solved numerically to determine the value of k. Then, substituting k into equation (28), the value of β can be determined, and therefore, y0 can be obtained from equation (28).


## 2.2. Model II: Positive and Negative Curvatures

For the first region, the curvature is positive atx0,0 up until the point xc,yc, where the curvature changes sign to be negative until the point a,y0 is reached as in Figures 5 and 6. At xc,yc, κ=0, and by solving equation (17), we have θc=−cos−1−γ/α, noting that from geometrical considerations, we have −π<θc<−ϕ/2. By making the substitution used in equation (19) for ϕ, we have ϕc=−π/2. By substituting ϕc into equations (22) and (25), we can calculate xc and yc,
(30)xc=x0−β2Ek−Kk,yc=2βk,noting that Kk and Ek are the complete elliptic integrals of the first and second kinds, respectively. In the second joining region which is between xc,yc and a,y0, we take the negative sign of equation (16), and following the same method used in Model I, we obtain
(31)xϕ=x0−β22Ek−E−ϕ,k−2Kk−F−ϕ,k,(32)yϕ=2βk1+cosϕ.Figure 5
Joining of BN graphene and BN nanotube.Figure 6
Joining of BN graphene and BN nanocone.As a result, from the boundary conditions at the pointa,y0, we know that ϕ=−ϕ0. By substitution into (32), we find
(33)x0−a=β22Ek−Eϕ0,k−2Kk−Fϕ0,k,(34)y0=2βk1+cosϕ0.In this model, the arc length is determined in two parts, which is given by(35)ℓ=∫θc01+ẏ21/2dx+∫θc−ϕ/21+ẏ21/2dx=β2Kk−Fϕ0,k.As a result, we find a dimensionless parameterμ=x0−a/ℓ for this case to be
(36)μ=22Ek−Eϕ0,k2Kk−Fϕ0,k−1.Also, for prescribed values ofx0,x1, and ℓ, we can solve equation (36) numerically to find the value of k, and by substituting k into equation (35), we can determine the value of β so that y0 can be determined from equation (34) [12]. Here, we note that equation (29) coincides with (36) for the value k=1/2, and the value of μ at this point is
(37)μ0=2E1/2K1/2−1=0.4569465810⋯.


## 3. Results

In this section, we investigate the numerical solutions for equations (29) and (37) when they are characterized by the nondimensional parameter μ=x0−a/ℓ, subjecting to the constraint −1<μ<1, as shown in [12]. Figure 7 shows that the three distinct regions are evident: the first one is where the parameter B is negative that is corresponding to Model I; the second region is 0<B≤2, which also corresponds to Model I; and the third region is 1<B≤2, which corresponds to Model II.Figure 7
The relation between the characteristic parameterμ and constant B for Model I and Model II [2].Results for the joining between BN graphene and BN nanotube for both models I and II are shown in Figures8 and 9, respectively. Furthermore, Figures 10 and 11 illustrate the join profile for BN graphene and BN nanocone with both models.Figure 8
Plot of the join curvey=yx of Model I for joining BN nanotube with BN graphene where B=−4 and ℓ=2.Figure 9
Plot of the join curvey=yx of Model II for joining BN nanotube with BN graphene where B=1.1 and ℓ=1.Figure 10
Plot of the join curvey=yx of Model I for joining BN nanocone 38.9° with BN graphene where B=−4 and ℓ=2.Figure 11
Plot of the join curvey=yx of Model II for joining BN nanocone 38.9° with graphene sheet where B=1.1 and ℓ=1.


## 4. Conclusion

In this paper, conventional applied mathematical modelling is considered to identify an approximate analytical solution to the shape of a surface joining BN graphene to other BN nanostructures. The resulted improved nanostructures could be useful for the design of probes for scanning tunnelling microscopy and other nanoscale devices. Using calculus of variations, this research minimizes the elastic energy of the joining curve which relates to minimization of the covalent bond energy. In particular, the joining is considered based on two different models that depend on the sign of the curvature: Model I refers to positive curvature only and Model II refers to positive and negative curvatures. The main purpose of this work is to formulate the axially symmetric model to have a reference basis for the comparison of real physical structures. As a result, these models lead to significative approximations to complex structures that are formulated by these certain BN nanostructures, BN graphene with two different BN nanostructures, BNNTs and BNNCs. Finally, other different BN nanostructures can be joined using the same model, and this is an area for future work.



---

*Source: 2394037-2020-08-30.xml*