EDN: GVVEKJ УДК 537.9
Effect of Permittivity on the Ionic Boundary Layer upon Protonation of Lithium Niobate
Vitaly A. Demin* Maxim I. Petukhovt Roman S. Ponomarev*
Perm State University Perm, Russian Federation
Mariana K. Kuneva§
Institute of Solid State Physics BAS Sofia, Bulgaria
Received 20.06.2023, received in revised form 06.07.2023, accepted 07.08.2023 Abstract. The paper presents the results of a numerical study of the behavior of proton exchange reaction products after they enter benzoic acid, which interacts with a lithium niobate crystal. The theoretical model based on the continuous media approximation considers the convective transfer and diffusion of positively charged lithium ions and negatively charged benzoate ions. Additionally, the possibility of their recombination with the subsequent formation of lithium benzoate molecules is taken into account. The results obtained correspond to stationary, uniform along the interface, profiles of the ion concentration and electric potential. In the process of establishment, benzoate ions form an ionic boundary layer. An analysis is made of the influence of the dielectric constant of the carrier liquid on its characteristics.
Keywords: proton exchange, boundary layer, permittivity, numerical simulation.
Citation: V.A. Demin, M.I. Petukhov, R.S. Ponomarev, M.K. Kuneva, Effect of Permittivity on the Ionic Boundary Layer upon Protonation of Lithium Niobate, J. Sib. Fed. Univ. Math. Phys., 2023, 16(5), 611-619. EDN: GVVEKJ.
Introduction
Technologies for production planar or channel waveguides [1] are associated with the implementation of the process of replacing positively charged lithium ions belonging to the lithium niobate crystal lattice with hydrogen ions. In many cases, the latter are formed as a result of thermal dissociation of benzoic acid molecules, which are brought into the contact with crystal at characteristic temperatures of ~ 500 K [2,3]. Other acids, such as palmitic [4], pyrophos-phoric [5,6], etc., can also be a source of protons. The duration of proton exchange can reach several hours. During this time, a near-surface region with an increased hydrogen concentration forms inside the substrate. The value of the refractive index in this region differs markedly from
*demin@psu.ru https://orcid.org/0000-0002-6095-1380
tgeniusmaxp@yandex.ru https://orcid.org/0000-0001-8771-6893
* rsponomarev@gmail.com https://orcid.org/0000-0001-9729-628X § m_kuneva@yahoo.com https://orcid.org/0000-0001-6088-407X © Siberian Federal University. All rights reserved
that which takes place in the bulk of the crystal. The thickness of the protonation region is - 10~6 m [7].
Despite the fact that this technology has been known for several decades, during which a large number of both experimental [1] and theoretical [8,9] works have been accumulated, there are still many questions related to the physical and chemical nature of the processes taking place during the mentioned reaction of substitution. So, at the moment, not all regularities of the formation of many crystalline phases in the protonated region of lithium niobate have been studied. It is known that the phase diagram is affected by both the addition of lithium benzoate to benzoic acid [10], and the temperature at which proton exchange occurs, as well as its duration [11]. In addition, the subject of research at the moment is the preliminary plasma treatment of the crystal surface and its effect on the final characteristics of the waveguide [12].
From a fundamental point of view, it seems important to study not only the processes that take place inside the crystal, but also those that occur in the acid. According to the results obtained during the experiments, the properties of the acid and its impurities have a very strong influence on the final result. For example, dilution of benzoic acid with lithium benzoate reduces the intensity of proton exchange. Already at 5% impurity concentration, the substitution reaction does not occur [10]. The addition of water to benzoic acid, which intensifies the exchange process, has a strong effect on the final waveguide [13].
In [14] a mathematical model that allows one to describe the behavior of ions, which are the products of the proton exchange reaction, in benzoic acid is presented. It was shown that despite the smallness of the region in which the reaction under consideration takes place, the processes occurring in it can be described using the equations of continuum mechanics. This model included the diffusion mechanism of impurity propagation and the recombination mechanism of ion interaction. In the course of solving the problem formulated in [14], the possibility of the existence of stationary ion concentration profiles was discovered. It should be noted that benzoate ions, which have a much lower diffusion coefficient than lithium ions, form the boundary layer. The disadvantages of this work may be that it did not consider the dynamics of the formation of concentration profiles, and also used the model value of the permittivity of the carrier liquid. In [14] it was equal to unit, which, obviously, does not correspond to the realistic value of the permittivity of benzoic acid, which is known to be a polar dielectric. This makes it necessary to refine the previous calculations in terms of quantitative determination of the parameters of the boundary layer.
1. Formulation of the problem
Let us consider a homogeneous surface of a lithium niobate crystal, on which a photolithographic mask is partially deposited (Fig. 1). This configuration is used to create channel waveguides. The resulting cavity, shown in Fig. 1, has dimensions h = 10~6 m, l = 5h [15] and is filled with benzoic acid heated to a temperature of T = 500 K. This temperature favors the conditions for thermal dissociation of molecules acids into benzoate ions C7H5O^ and hydrogen ions H+. The latter tend to penetrate deep into the crystal with further replacement of lithium ions Li+. Thus, from the side of the interface there is a flow of positively charged lithium ions into the acid. After some time, they should meet with benzoate ions and form a lithium benzoate C7H5O2Li molecule.
The estimations in [14] make it possible to say that the equations of concentration convection written in the electroconvective approximation allow to describe the process of ion transport and
Fig. 1. Statement of the problem. System of coordinate
their further recombination using the following system of equations [16,17]:
+ (vV)v^ = —Vp + nAv + e(n+ — n-)E,
divv = 0, div(e0eE ) = e (n+ — n-), E = —V dn±
dt
+ (vV)n± = D±An± ^ V(k±n±E) — kRn+r
(1) (2) (3)
where v, E, y, p, n± are the fields of velocity, electric field strength and electric potential, pressure, concentration of positively and negatively charged ions. Parameters n, e, £, £o are the dynamic viscosity of benzoic acid, electron charge, electrical constant and permittivity. In turn, D±, k±, and kR set the diffusion coefficients of lithium and benzoate ions, their mobility and the recombination coefficient, the value of which was estimated according to [17]:
kR
e(k+ + k-)
£0£
(4)
The concentration fluctuations that may occur in such hydrodynamic system are not able to cause gravitational convection [18], therefore, there is no corresponding term in equation (1). However, unlike the statement presented in [14], the system of equations (1)-(3) will be solved without neglecting the convective terms, as well as the time derivative, which will allow us to study the dynamics of the formation of concentration and electric field profiles.
The walls of the photolithographic mask and the crystal surface were considered to be solid; therefore, the no-slip condition was imposed on them. A constant mass flow of ions of both signs J = 1018 m~2s_1 [14] was maintained at the lower boundary. The remaining boundaries for ions were impenetrable. The electric potential was measured relative to the lower limit.
0, l :
v = 0,
dn± dx
0
' dx
0,
y = 0 :
v = 0,
dn± J
-W = —D±, * = °-
The upper bound of the region in question was free, so
h :
dvx dy
0,
0,
dn± dy
0,
d7T =0-
dy
(5)
(6)
(7)
The fact that the same flux density J is used in boundary conditions (6) for both negatively charged and positively charged ions is due to the overall electroneutrality of the system, since
x=
y
v
y
one hydrogen ion replaces one lithium ion. In other words, the fulfillment of condition (6) is associated with the equality
l fh A fh
n+dxdy = / n-dxdy. (8)
0 J 0 J 0 J 0
The boundary value problem, which includes the system of equations (1)-(3) in combination with boundary conditions (5)-(7), was solved in dimensionless variables. h,v/h,h2/v, Jh/D+, eJh3/e0eD+, pv2/h2, respectively. In these units, v = n/p is the coefficient of kinematic viscosity.
After the exclusion of the electric field strength and procedure of adimensionalization , the equations and boundary conditions have the form:
d v
— + (vV)v = -Vp + Av + A(n+ - n-)Vy, (9)
divv = 0, A^ = (n+ - n-), (10)
ddn± + (vV)n± = s~An± t B±V(n±Vp) - Cn+n-, (11)
x = 0, l : v = 0, =0, ^ =0, (12)
dx dx
dn+ dn_ D+
y=0: v = 0, dn+ =-1, =-D±, *=0, (13)
U dvx n n dn± n n C\A\
y =h : dy = °, v = °, = 0, dy = 0• (14)
The dimensionless control parameters are defined as follows:
A _ e2h6.J2 _ k±eJh3 C _ kRJh3 S (15)
e0epv 2D+ £q£VD+ VD+ ' D±
The formulated problem admits a two-dimensional solution in terms of the stream function 0 and the vorticity 0:
d0 d0 ^ dvx dvy (16)
x By' y dX dy dx
After eliminating velocity and pressure using (16), equations (10)-(13) take the final form:
d0 + d0d0 30 30 _ d20 + d20 + A(3A dp 3A dp\ (17)
dt dy dx dx dy dx2 dy2 \ dy dx dx dy J1
d2p d2p d20 d20
JTi + _ n- - n+ _ A' JT + XT _ 0 (18)
dx2 dy2 dx2 dy2
dn± d0 dn± d0 dn± 1 i d2n± d2 n± \
dt dy dx dx dy Sc± \ dx2 dy2 )
+B I dn± dp + dn± dp + / d2p + d2p\ \ ^
I dx dx dy dy \ dx2 dy2 J I + '
2. Solution technique
The boundary value problem (17)-(19), (12)-(14) was solved numerically by the two-field [19] method implemented using an explicit scheme in the Fortran programming language. One-sided differences were used in the approximation of derivatives with respect to coordinates at the boundaries and derivatives with respect to time. Spatial derivatives in the internal nodes of the computational domain were approximated by central differences. To solve the Poisson equations (18), the method of successive upper relaxation was used. The method of establishment was used in the calculations. The time step was chosen from considerations of the stability of the numerical procedure. The stream function at all boundaries was equal to zero. The vorticity value at the upper boundary was also equal to zero. The vorticity at solid boundaries was determined using the Thom condition. To resolve the boundary layer, the computational domain was discretized using a uniform grid consisting of 40 nodes along the x axis and 100 nodes along the y axis.
The choice of control parameter values was made in accordance with [14]: h = 2 • 10_6 m, k+ = 1.5-10~7 m2/s-V, k_ =2-10-8 m2/s-V, D+ = 10_8 m2/s, D_ = 10_9 m2/s, v = 5 • 10_7 m2/s, p = 103 kg/m 3, J = 1018 s-1m_2.
Separately, it is necessary to comment on the choice of values of the permittivity. At the moment, there are no data in the literature on direct measurements of this value for pure benzoic acid, even under normal conditions, when it is in the solid state. Moreover, according to the available data, the permittivity of benzoic acid varies greatly when impurities are added to the acid [20], which only complicates the search for specific values, given the wide variety of ways to create waveguides [1]. To estimate e, it only remains to find an analogue of a substance whose molecules are identical in their properties to those of benzoic acid and whose permittivity would be measured at the required temperatures corresponding to the proton exchange regimes. Thus, considering the values of dipole moments d of molecules, one can compare benzoic acid (d = 1.78 D [21]) with methanol (d = 1.70 D [22]). The permittivity of the latter was measured over a very wide range of temperatures and pressures in [23]. According to these results, e at temperatures — 500 K can be taken equal to 10. It should be noted that this value was obtained at a pressure — 200 atm, which is very different from the proton exchange conditions. However, according to [23], the dielectric constant of methanol does not strongly depend on pressure, so we can stop at this estimation. When solving the problem, several values of the permittivity were used to potentially consider different ways of protonation. Based on these values, the recombination coefficient kR was also recalculated.
3. Results of numerical simulation
The performed calculations demonstrate a very fast damping of all hydrodynamic perturbations and reaching a stationary regime, regardless of the value of the permittivity (Fig. 2). In dimensional units, the system needs — 10_4 s to settle in a stationary state. It should be noted that these results confirm the assumptions made in [14] regarding the absence of convective mass transfer in the stationary case, since all fields after reaching this regime are uniform along the interface. In other words, this problem can be solved with good accuracy in a one-dimensional approximation. In the future, to visualize the results, only profiles along the y axis for the center of the computational domain will be presented.
An analysis of the concentration profiles shows that negatively charged benzoate ions form a boundary layer, the thickness of which increases with permittivity growth (Fig. 3). This is mainly due to the fact that the recombination coefficient depends on e. An increase in permeability
4E-005
0.3 -i
b
a
2E-005
3E-005
1E-005
0
. i o +
0.3 t 0.4 0
0
0.1
0.2
10
20
30
40 / 50
Fig. 2. Dependencies of stream function maximum value (a) and concentration in the point with coordinates^ = 2.5, y = 0.2 (b) on demensionless time: 1 — e =1, 2 — e =10. On the figure b solid lines — benzoate ions concentration, dashed lines — lithium ions concentration
reduces kR, which makes it more likely that benzoate ions are located at large distances from the crystal. As the results show, on average, the thickness of this layer lies in the range from ~ 10-7 m to ~ 10~6 m. On the other hand, lithium ions, which have higher mobility, have the ability to penetrate deep into benzoic acid at much greater distances. In fact, they occupy the entire computational domain and are held by the electric field that arises according to (11) due to the difference in ion concentrations. The intensity of this field is directed towards the substrate, i.e. the ionic boundary layer is slightly stretched. At the same time, in the concentration profiles of lithium ions, one can notice a bend in the region where benzoate ions are no longer present. This bending is also associated with the action of the electric field, which, on the contrary, tends to "press" the positively charged ions to the substrate.
The electric field also undergoes changes when the permittivity is taken into account (Fig. 4). As can be seen from Fig. 3 an increase in e leads to an increase in the concentration difference. As a result, the potential difference at realistic permeability values can be increased by an order of magnitude compared to the model value e = 1: from ~ 10~4 V to ~ 10~3 V.
Conclusion
The calculation results allow us to state that taking into account realistic values of the permittivity increases the thickness of the ionic boundary layer formed by benzoate ions during proton exchange. If in the model case e = 1 the layer thickness was ~ 2 • 10~7 m, then as the permeability value increases to e = 10, negatively charged ions can be detected at a distance of ~ 4 • 10-7 m from the crystal. Thus, the thickness of the boundary layer becomes comparable with the thickness of the protonated region inside the substrate. In addition, it can be concluded that the use of acids with a high dielectric permittivity makes it possible to influence the "outer" boundary layer to a greater extent during protonation by generating an acid flow along the substrate by swinging the workpiece, which is often used in production.
0 0.2 0.4 0.6 0.8 y 1
Fig. 4. Electric potential profiles in the point x = 2.5 at the time t = 50. 1 — e = 1, 2 — e = 10,
3 — e = 20, 4 — e = 30
The work was supported by the government of the Perm Region, project no. C-26/848.
References
[1] M.Kuneva, Optical waveguides obtained via proton exchange technology in LiNbO3 and LiTaO3 - a short review, International Journal of Scientific Research in Science and Technology, 2(2016), 40-50.
[2] J.L.Jackel, C.E.Rice, J.J.Veselka, Proton exchange for high-index waveguides in LiNbO3, Appl. Phys. Lett., 41(1982), 607-608.
[3] J.L.Jackel, Proton exchange: past, present, and future, Proc. SPIE, 1583(1991), 54-63.
[4] M.De Micheli, J.Botineau, S.Neveu, P.Sibillot, D.B.Ostrowsky, Independent control of index and profiles in proton-exchanged lithium niobate guides, Optics Lett., 8(1983), 114-115.
[5] E.Y.B.Pun, T.C.Kong, P.S.Chung, H.P.Chan, Index profile of proton-exchanged waveguides in LiNbO3 using pyrophosphoric acid, Electr. Lett., 26(1990), 81-82.
[6] N.Goto, G.L.Yip, Characterization of proton-exchange and annealed LiNbO3 waveguides with pyrophosphoric acid, Appl. Optics, 28(1989), 60-65.
[7] I.V.Petukhov, V.I.Kichigin, A.P.Skachkov, S.S.Mushinsky, D.I.Shevtsov, A.BVolyntsev, Microindentation of proton exchange layers on X cut of lithium niobate crystals, Materials Chemistry and Physics, 135(2012), 493-496.
[8] S.T.Vohra, A.R.Mickelson, S.E.Asher, Diffusion characteristics and waveguiding properties of proton-exchanged and annealed LiNbO3 channel waveguides, J. Appl. Phys., 66(1989), 5161-5174.
[9] V.A.Demin, M.I.Petukhov, R.S.Ponomarev, A.V.Topova, On a role of anisotropy and nonlinear diffusive effects during the construction of waveguides in the lithium niobate, Bulletin ofPerm University. Physics, 1(2021), 49-58. DOI: 10.17072/1994-3598-2021-1-49-58
[10] V.I.Kichigin, I.V.Petukhov, S.S.Mushinsky, V.A.Oborin, A.M.Minkin, L.N.Malinina, D.I.Shevtsov, A.B.Volyntsev, Structure and properties of proton exchange waveguides on Z cut of lithium niobate crystal fabricated in molten benzoic acid with the addition of lithium benzoate, International Conference and Seminar of Young Specialists on Mi-cro/Nanotechnologies and Electron Devices, 2012, 238-241.
[11] Yu.N.Korkishko, V.A.Fedorov, Structural phase diagram of HœLii_œNbO3 waveguides: the correlation between optical and structural properties, IEEE J. Sel. Top. Quantum Electron, 2(1992), 187-196.
[12] I.V.Petukhov, V.I.Kichigin, S.S.Mushinskii, D.I.Sidorov, O.R.Semenova, The influence of plasma treatment of lithium niobate crystal surface on the proton exchange process in molten benzoic acid, Vestnik Permskogo universiteta. Seriya "Khimiya", 9(2019), 371-379 (in Russian).
[13] S.S.Mushinsky, A.M.Minkin, I.V.Petukhov, V.I.Kichigin, D.I.Shevtsov, L.N.Malinina, A.B.Volyntsev, M.M.Neradovskiy, V.Ya.Shur, Water Effect on Proton Exchange of X-cut Lithium Niobate in the Melt of Benzoic Acid, Ferroelectrics, 476(2015), 84-93.
DOI: 10.1080/00150193.2015.998530
[14] V.A.Demin, M.I.Petukhov, R.S.Ponomarev, An ionic boundary layer near the lithium niobate surface in the proton exchange process, Surface Engineering and Applied Electrochemistry, 59(2023), 321-328. DOI:10.3103/S1068375523030055
15] S.S.Mushinsky, V.I.Kichigin, I.V.Petukhov, etc., Structure and properties of proton exchanged layers in +Z cut and -Z cut lithium niobate, Ferroelectrics, 443(2013), 20-34. DOI: 10.1080/00150193.2013.773856
16] L.D.Landau, E.M.Lifschitz, Fluid Mechanics, Course of Theoretical Physics, Vol. 6, Oxford, Pergamon press, 2006.
17] F.Pontiga, A.Castellanos, Physical mechanisms of instability in a liquid layer subjected to an electric field and a thermal gradient, Phys. Fluids, 6(1994), 1684.
18] G.Z.Gershuni, E.M.Zhukhovitskii, Convective stability of incompressible fluids, Jerusalem, Keter Publishing House, 1976.
19] E.L.Tarunin, Vychislitel'nyj eksperiment v zadachah svobodnoj konvekcii, Izdatel'stvo Irkutskogo universiteta, 1990.
20] A.Yurquina, M.E.Manzur, P.Brito, R.Manzo, M.A.A.Molina, Solubility and dielectric properties of benzoic acid in a binary solvent: water-ethylene glycol, Journal of Molecular Liquids, 108(2003), 119-133.
21] F.P.Parungo, J.P.Lodge, Molecular structure and ice nucleation of some organics, Journal of Atmospheric Sciences, 22(1965), 309-313.
22] M.Jorge, J.R.B.Gomes, M.C.Barrera, The dipole moment of alcohols in the liquid phase and in solution, Journal of Molecular Liquids, 356(2022), 119033. DOI: 10.1016/j.molliq.2022.119033
23] E.U.Franck, R.Deul, Dielectric behaviour of methanol and related polar fluids at high pressures and temperatures, Faraday Discussions of the Chemical Society, 66(1978), 191-198.
Влияние диэлектрической проницаемости на ионный пограничный слой при протонировании ниобата лития
Виталий А. Демин Максим И. Петухов Роман С. Пономарев
Пермский государственный университет Пермь, Российская Федерация
Мариана К. Кунева
Институт физики твердого тела БАН София, Болгария
Аннотация. В работе представлены результаты численного исследования поведения продуктов реакции, попадающих в бензойную кислоту, которая взаимодействует с кристаллом ниобата лития при протонном обмене. В теоретической модели, построенной в приближении сплошной среды, рассматривается конвективный перенос и диффузия положительно заряженных ионов лития и отрицательно заряженных бензоат-ионов. Дополнительно учитывается возможность их рекомбинации с последующим формированием молекул бензоата лития. В процессе установления бензоат-ионы образуют ионный пограничный слой. В этом пределе полученные результаты отвечают стационарным, однородным вдоль межфазной поверхности профилям концентрации ионов и электрического потенциала.
Ключевые слова: протонный обмен, пограничный слой, диэлектрическая проницаемость, численное моделирование.