Home A DFT study on the molecular properties of synthetic ester under the electric field
Article Open Access

A DFT study on the molecular properties of synthetic ester under the electric field

  • Yachao Wang , Xiaoran Lin EMAIL logo , Mei Wang and Jifang Wang
Published/Copyright: November 6, 2021

Abstract

Synthetic ester can replace the mineral oil traditionally used in transformers to avoid the environmental problems caused by oil leakage. However, the fast discharge phenomenon in a high electric field in transformers using synthetic ester seems to indicate its insulation property is inferior to that of mineral oil. In this paper, typical molecular models of synthetic ester, including F2, F4, F6, F8, and F10, are constructed. We studied the effect of electric fields on the molecular properties of the five molecules by density functional theory and time-dependent density functional theory. According to the electric field intensity required for discharge initiation and propagation in insulating oil, the electric field intensity applied in this study varied from 108 to 109 V/m. The results showed that the molecular bond lengths are obviously dependent on the electric field. The ionization potential (IP) of the F8 and F10 molecules decreases sharply under electric field intensities of 3.1 × 109 and 4.0 × 109 V/m. It can be inferred that the IP reduction of the long carbon chain molecules, such as F8 and F10, is the reason for the formation of fast discharge in the case of synthesis ester. Calculations for excited states show that the introduction of an electric field makes the electron transition more active. The results obtained by this work improve our understanding of the discharge mechanism in synthetic ester dielectrics and provide theoretical support for improvement in the performance of synthetic ester insulating oil.

1 Introduction

Synthetic ester is obtained by chemical synthesis and refining and is widely used as lubricating fluid in automobiles, aircraft, and other machinery [1]. In recent years, ester materials have been used as transformer oil for improved environmental performance [2,3,4,5]. Synthetic ester is degradable and can replace the traditional mineral oil in transformers to avoid the environmental problems caused by oil leakage [6,7]. Compared with traditional transformer oil, the higher ignition point of synthetic ester increases the security of power grids [8,9]. At present, synthetic ester is used in distribution transformers [10]. However, the application of synthetic ester in large power transformer needs further work. Discharge tests show that the breakdown voltage of synthetic ester is almost the same as that of traditional transformer oil under the condition of short oil gap [11]; however, under the condition of long oil gap, the breakdown characteristic of synthetic ester is obviously worse than with traditional mineral oil [12,13]. The fast discharge phenomenon is quick to appear in the case of synthetic ester in a high electric field [14]. Obviously, the molecular characteristics of synthetic ester in an electric field must be different from that of traditional mineral transformer oil. Basic physical research at the molecular level in this regard is still necessary.

The molecular structure of synthetic ester results from the esterification of a pentaerythritol molecule with four fatty acid molecules [14,15]. Figure 1 shows the molecular structure of synthetic ester. To improve the stability of materials, the fatty acid of synthetic ester used for insulating oil is usually saturated fatty acid. Functional group COOR can be the same or different. In this study, five fatty acids, including acetic acid (CH3COOH), butyric acid (CH3(CH2)2COOH), hexanoic acid (CH3(CH2)4COOH), octanoic acid (CH3(CH2)6COOH), and decanoic acid (CH3(CH2)8COOH), were selected to construct the synthetic ester molecule. If each synthetic ester contains only one kind of fatty acid, we can obtain five typical synthetic ester molecules, labeled as F2, F4, F6, F8, and F10. These five typical synthetic ester molecules were the objects of this research study.

Figure 1 
               Molecular structure of synthetic ester.
Figure 1

Molecular structure of synthetic ester.

Density functional theory (DFT) can be used to study the molecular properties of insulating materials [16,17,18]. The effect of C═C double bond on the ionization potential (IP) and electron affinity of unsaturated triglycerides was studied in ref. [16]; In ref. [17], the electronic structures and IPs of triolein and aromatic were compared, and the differences in discharge characteristics between natural ester and mineral oil were explained. The effect of electric fields on the IP and excitation energy of aromatic and alkane molecules was studied in ref. [18]. At present, few studies have been reported on the molecular characteristics of synthetic ester under electric field.

In this study, we carry out a DFT study on the molecular properties of synthetic ester under the electric field. The molecules include F2, F4, F6, F8, and F10. The molecular properties refer to bond length, infrared spectrum, total energy, dipole moment, molecular orbital, IP and the excited state. Existing experimental results show that the electric field intensity required is 108 V/m for discharge initiation and 109 V/m for fast discharge [19,20,21]. So, the electric field intensity applied in this study varies from 108 to 109 V/m. The results obtained in the work are helpful to improve the understanding of discharge mechanism in synthetic ester dielectrics and provide theoretical support for the performance improvement of synthetic ester insulating oil.

2 Theoretical method

The Hamiltonian H of molecule and electric field system can be expressed as [22]

(1) H = H 0 + H int ,

where H 0 is the Hamiltonian without the electric field and H int is the Hamiltonian of the interaction between external electric field and molecule. H int can be expressed as

(2) H int = μ F ,

where μ is the dipole moment and F is the electric field.

Time-dependent density functional theory (TD-DFT) is also used to study the first six excited states. The expression for excitation energy E ex can be written as

(3) E ex ( F ) = E ex ( 0 ) Δ μ F 1 2 Δ α F 2 ,

where Δ μ and ∆α are the variation of the dipole moment and the polarizability, respectively [23,24]. All the tests are carried out by the Gaussian 09 package [25].

Figure 2 shows the model of the F2 molecule. The red, blue, and gray atoms in the figure represent oxygen, hydrogen, and carbon, respectively, which are also marked by numbers. The F4, F6, F8, and F10 molecules have 2, 4, 6, and 8 more carbon atoms in each branch chain than the F2 molecule. The z-axis runs through the two atoms 1C and 2C. The molecular models of F2, F4, F6, F8, and F10 are optimized under the condition of without the electric field. Then, the electric fields of 0.0006, 0.0012, 0.0024, 0.0048, 0.0060, and 0.0078 a.u. are applied in the z-axis direction respectively. The electric field intensity of 1 a.u. stands for 5.1 × 1011 V/m. Consequently, the applied electric field are also 3.1 × 108, 6.2 × 108, 1.2 × 109, 2.5 × 109, 3.1 × 109, and 4.0 × 109 V/m, correspondingly. The molecular models are optimized again when the electric field is applied. Bond lengths, infrared spectra, total energy, dipole moment, and orbital energy are studied. Finally, the first six excited states are calculated by using the optimized molecular structure.

Figure 2 
               Model of F2 molecule with Cartesian axis.
Figure 2

Model of F2 molecule with Cartesian axis.

In terms of specific methods and basis set selection, this paper compares the results of SE/MP6, HF/6-31 G*, and B3LYP/6-31 + G* for the vibrational frequencies of main functional groups of molecules without electric field. 1.06, 0.90, and 0.96 are the scale factors [26]. The results are listed in Table 1. The experimental data for the functional groups are also listed in the table [27,28]. It can be found that the calculation results of B3LYP/6-31 + G* are closer to the experimental values than those of SE/MP6 and HF/6-31 G* scheme. The B3LYP/631 + G* method is selected for all other calculations in this paper.

Table 1

Comparison of the three methods

Molecule Method [26,27,28] Wave number (cm−1)
C–O–C C═O C–H
SE/MP6/1.06 1320.5 1936.6 2784.0
F2 HF/631 + G*/0.90 1250.6 1809.9 2945.7
B3LYP/631 + G*/0.96 1194.0 1751.2 2956.6
SE/MP6/1.06 1190.2 1931.8 2770.5
F4 HF/631 + G*/0.90 1187.6 1802.9 2873.7
B3LYP/631 + G*/0.96 1125.6 1744.7 2985.8
SE/MP6/1.06 1186.8 1932.4 2768.0
F6 HF/631 + G*/0.90 1181.6 1802.7 2868.2
B3LYP/631 + G*/0.96 1120.5 1744.5 2974.4
SE/MP6/1.06 1293.7 1932.4 2767.9
F8 HF/631 + G*/0.90 1179.2 1802.7 2926.1
B3LYP/631 + G*/0.96 1118.6 1744.5 2973.2
SE/MP6/1.06 1189.1 1932.4 2767.9
F10 HF/631 + G*/0.90 1178.2 1802.7 2868.1
B3LYP/631 + G*/0.96 1117.6 1744.5 2901.2
Experimental value 1,150–1,330 1,735–1,750 2,845–2,975

3 Results and discussion

3.1 Bond length and infrared spectra

Bond length is the basic structural parameter of molecules, and the change of bond length plays a key role in the influence of molecular properties. Synthetic ester molecules contain a large number of chemical bonds. Limited to space, the bond lengths of some main chemical bonds of the F2 molecule under different electric fields are given in Table 2. It was found that the bond lengths are obviously dependent on the electric field. When the electric field intensity varies from 0 to 4.0 × 109 V/m, 2C–14O increases from 1.4380 to 1.4437 Å,and the bond length is extended by 0.0057 Å; 18C═19O increases from 1.2083 to 1.2148 Å,and the bond length is extended by 0.0065 Å. From Figure 2 and Table 2, it can be seen that the change of the same type of chemical bond under the electric field is not the same. That is because the angle between the chemical bond and the electric field direction is different. Generally, the variation of C–O with electric field intensity is larger than that of C═O and that of C–C.

Table 2

Bond length of tributyrin versus applied electric field

E (108 V/m) 1C–2C (Å) 2C–14O (Å) 11C–16O (Å) 18C═19O (Å) 18C–38C (Å) 20C═21O (Å) 24C–17O (Å) 24C═25O (Å)
0 1.5410 1.4380 1.4380 1.2083 1.5130 1.2083 1.3670 1.2083
3.1 1.5406 1.4385 1.4381 1.2087 1.5136 1.2085 1.3682 1.2077
6.2 1.5402 1.4391 1.4383 1.2091 1.5141 1.2087 1.3695 1.2070
12 1.5393 1.4404 1.4386 1.2099 1.5152 1.2091 1.3720 1.2057
25 1.5414 1.4405 1.4414 1.2121 1.5166 1.2111 1.3773 1.2031
31 1.5410 1.4417 1.4419 1.2131 1.5172 1.2117 1.3800 1.2018
40 1.5410 1.4437 1.4426 1.2148 1.5181 1.2128 1.3842 1.1999

The reason for the change of bond length is related to the electron charge transfer under the electric field. As we know, the atom is stable at the equilibrium distance under the internal stress of gravity and repulsion. The equilibrium distance is the bond length. When the electric field is applied, the amount of electronic charge on the atom’s transfers, the internal stress and equilibrium distance between atoms change, and the bond length changes accordingly.

Figure 3 shows the infrared spectrum of the F2 molecule under the electric field intensities of 0, 3.1 × 108, 1.2 × 109, and 3.1 × 109 V/m. In the absence of the electric field, the infrared spectrum of the F2 molecule is shown in Figure 2(a). The wavenumbers 1194.0 and 1194.3/cm are mainly contributed by the stretching vibration of C–O–C. The wave number 1751.1/cm is mainly contributed by the stretching vibration of C═O. The wavenumber range of 2938.6–3051.1/cm is mainly contributed by the stretching vibration of C–H.

Figure 3 
                  Infrared spectra of the F2 molecule versus different electric field intensities, including 0 V/m (a), 3.1 × 108 V/m (b), 1.2 × 109 V/m (c) and 3.1 × 109 V/m (d).
Figure 3

Infrared spectra of the F2 molecule versus different electric field intensities, including 0 V/m (a), 3.1 × 108 V/m (b), 1.2 × 109 V/m (c) and 3.1 × 109 V/m (d).

When the electric field is applied, the infrared spectrum of the F2 molecule changes obviously, as shown in Figure 2(b–d). When the electric field changes from 0 to 3.1 × 109 V/m, the vibration frequency of C–O undergoes a redshift, and the wave number peak decreases from 1194.0–1194.3/cm to 1165.2–1190.4/cm. This is due to the increase of C–O bond length with the increasing of electric field intensity. The bond length of 2C–14O, 11C–16O, and 24C–17O is extended by 0.0037, 0.0039, and 0.0130 Å, respectively. With the increasing of the electric field intensity, the wave number peaks of C═O and C–H are broadened and the peak value decreases. In the absence of the electric field, C═O vibration has a main peak of 1751.1/cm; when an electric field intensity of 3.1 × 109 V/m is applied, the wave number peak of C═O vibration ranges from 1727.4 to 1780.3/cm. This is due to the changes of C═O bond length with the increasing of electric field intensity, which is similar to the reason mentioned above.

3.2 Total energy and dipole moment

Figure 4 shows the dipole moments of the five synthetic ester molecules under the different electric fields. It can be found that the dipole moments of these five molecules are all zero in the absence of the electric field, and they are all nonpolar molecules. When the electric field is applied, the positive and negative charge centers of the molecule no longer coincide, and the dipole moment of the molecule increases with the increasing of the electric field intensity. It should be noted that the dipole moment of molecule increases with the increasing of carbon chain length under the same electric field. For example, under the electric field intensity of 4.0 × 109 V/m, the dipole moment of the F10 molecule is more than two times that of the F2 molecule.

Figure 4 
                  Dipole moment versus applied electric field.
Figure 4

Dipole moment versus applied electric field.

With the increasing of the carbon chain length, the volume of triglyceride molecule increases correspondingly. The distance between the positive and negative charge centers of the molecule increases correspondingly under the electric field, resulting in the increasing of dipole moment. In liquid dielectrics, the larger polarity of molecules will bring larger dielectric constant and dielectric loss, and reduce the insulation performance of dielectrics [29,30].

Figure 5 shows the total energy of the five synthetic ester molecules under the different electric fields. With the increasing of electric field intensity, the total energy decreases, which is contrary to the trend of dipole moment with the electric field. As mentioned in Section 2, according to formulas (1) and (2), the total molecular energy is negatively correlated with its dipole moment, so the total molecular energy will be the trend shown in Figure 5.

Figure 5 
                  Total energy versus applied electric field.
Figure 5

Total energy versus applied electric field.

3.3 Orbital energy

The molecular frontier orbitals, including the highest occupied molecular orbital (HOMO) and the low unoccupied molecular orbital (LUMO), have an important influence on the molecular properties [31]. Figure 6 shows the frontier orbital energies of the five synthetic ester molecules under the different electric fields. The LUMO energies of these five molecules decrease with the increasing of the electric field intensity. The reduction of LUMO energy indicates the enhancement of the ability of molecules to obtain electrons. The HOMO energies of these five molecules increase with the increasing of the electric field intensity. when the same electric field intensity is applied, the HOMO energy increases with the increasing of the carbon chain length. The enlargement of HOMO energy indicates that the molecule is easier to lose the orbital electrons.

Figure 6 
                  Frontier molecular orbital energy versus applied electric field.
Figure 6

Frontier molecular orbital energy versus applied electric field.

The LUMO energy minus the HOMO energy is the energy gap. The results are shown in Figure 7. The energy gap reflects the activity of molecules to a certain extent [32]. It can be seen from Figure 7 that the energy gaps of these five molecules decrease with the increasing of the electric field intensity. When the same electric field intensity is applied, with the increasing of carbon chain length, the molecular energy gaps also decrease.

Figure 7 
                  Energy gap versus applied electric field.
Figure 7

Energy gap versus applied electric field.

Figure 8 shows the frontier orbital electron cloud of the F2 molecule under the electric field. The frontier orbital composition of the F2 molecule is given in Table 3, which is calculated by using Multiwfn 3.3.8 software package [33]. It can be seen from Figure 1, Figure 7, and Table 3 that the molecular configuration is obviously distorted under the electric field. In the absence of the electric field, the HOMO of the F2 molecule is mainly contributed by 19O(10.7%), 21O(13.2%), 23O(13.1%), and 25O(10.9%); the LUMO is mainly contributed by 1C(60.8%).

Figure 8 
                  Electron cloud of the HOMO and LUMO of the F2 molecule versus applied electric field.
Figure 8

Electron cloud of the HOMO and LUMO of the F2 molecule versus applied electric field.

Table 3

Composition of the frontier molecular orbitals of the F2 molecule versus applied electric field

E (108 V/m) Composition of the frontier orbitals (>2%)
0 HOMO 1C:4.4, 2C:4.5, 5C:5.7, 8C:4.4, 11C:5.6, 14O:2.3, 15O:2.9, 16O:2.9, 17O:2.4, 19O:10.7, 21O:13.2, 23O:13.1, 25O:10.9, 26C:3.0, 34C:2.5, 38C:2.4
LUMO 1C:60.8, 8C:2.1, 11C:2.1, 18C:5.8, 22C:5.7, 24C:5.7
3.1 HOMO 1C:4.6, 2C:4.4, 8C:14.9, 14O:2.2, 17O:7.8, 19O:9.8, 25O:34.5, 34C:7.7, 38C:2.4
LUMO 1C:47.7, 2C:2.0, 5C:3.4, 8C:1.7, 11C:3.7, 18C:14.7, 24C:9.8, 38C:5.1
6.2 HOMO 1C:4.8, 2C:2.2, 8C:17.4, 17O:9.3, 19O:4.1, 24C:5.0, 25O:41.8, 34C:8.9
LUMO 1C:23.4, 2C:8.2, 5C:4.0, 8C:2.3, 11C:4.6, 18C:22.9, 24C:10.0, 26C:2.3, 30C:2.3, 38C:12.8
12 HOMO 1C:4.6, 8C:16.9, 17O:10.3, 24C:5.6, 25O:45.5, 34C:9.6
LUMO 2C:26.6, 8C:3.8, 11C:2.3, 18C:26.1, 24C:4.8, 26C:2.7, 30C:3.0, 38C:22.8
25 HOMO 1C:6.3, 8C:8.1, 17O:12.2, 24C:6.4, 25O:51.5, 34C:10.9
LUMO 1C:4.4, 2C:39.8, 5C:24.4, 18C:12.7, 38C:12.3
31 HOMO 1C:6.2, 8C:6.3, 17O:12.6, 24C:6.5, 25O:52.4, 34C:11.0
LUMO 1C:4.3, 2C:38.2, 5C:23.3, 18C:13.9, 38C:14.7
40 HOMO 1C:5.9, 8C:4.3, 17O:13.2, 24C:6.6, 25O:53.1, 34C:11.2
LUMO 1C:4.1, 2C:33.3, 5C:21.5, 11C:2.3, 18C:15.9, 38C:18.8

When the electric field is applied, with the increasing of the electric field intensity, the electron cloud of HOMO moves in the opposite direction of the electric field, and the electron cloud of LUMO moves in the positive direction of the electric field. When an electric field intensity of 4.0 × 109 V/m is applied, the HOMO is mainly contributed by 25O(53.1%), 17O(13.2%), and 34C(11.2%); the LUMO is mainly contributed by 2C(33.3%), 5C(21.5%), 38C(18.8%), and 18C(15.9%). The movement of the frontier orbit to the far end of the molecule under the electric field makes it easier for the molecule to gain or lose electrons.

The discharge process in liquid dielectrics is very complicated. Under the strong electric field, the neutral molecules of dielectric ionize to form positive ions and electrons. The positive ions and electrons move to the cathode and anode respectively under the action of electric field force. In this process, electrons can be trapped by neutral molecules to form negative ions. During the whole discharge process, there occurs recombination of positive and negative ions, recombination of positive ions and electrons, and strong acousto-optic phenomena concomitantly. The ionization of neutral molecules runs through the whole discharge process and plays an important role in the development of discharge. Molecular IP can represent the difficulty of molecular ionization. Koopman’s theorem holds that the molecular IP is equal to the opposite number of HOMO energy [34]. According to this theorem, the IP of molecules can be estimated. Koopman’s theorem assumes that when a molecule loses an electron, its molecular structure and orbital energy level remain unchanged, so the IP obtained from this theorem may not be very accurate. However, in liquid dielectrics composed of various molecular components, the discharge process is usually determined by molecules with low IP, so it is still of practical significance to compare the relative value of intermolecular IP [35].

Figure 9 compares the IPs of the five molecules under the different electric fields. The molecular IP decreases with the increasing of carbon chain length. This is because with the increasing of the carbon chain length, the molecular volume becomes larger, the molecule’s attraction to the outermost electrons is weakened, and the molecule is more likely to lose the electrons. It is worth noting that the IPs of the F8 and F10 molecules decrease sharply under the electric field intensities of 3.1 × 109 and 4.0 × 109 V/m. The IP of the F10 molecule is 5.26 eV under the electric field intensity of 4.0 × 109 V/m, which is 2.50 eV lower than that without the electric field, and the decrease rate is more than 30%. The relationship between the molecular IP and the electric field intensity can be expressed as IP ( E ) = IP ( 0 ) γ E , where IP(0) is the molecular IP under the electric field intensity of 0 V/m and γ is a constant [36]. The IP of molecules decreases and the charge production increases under the strong external electric field, which leads to the rapid development of discharge in liquid dielectrics [35,37]. From this point of view, it can be inferred that the long carbon chain molecules, such as the F8 and F10 molecules, are the reason for the formation of fast discharge in the synthesis ester.

Figure 9 
                  IP of different molecules under electric field.
Figure 9

IP of different molecules under electric field.

Tunnel ionization may occur under an intense field, accompanied by molecular dissociation [38,39]. Therefore, the applied electric field intensity has an upper limit. When it exceeds this upper limit, molecular dissociation may occur. The tunnel ionization rate of the molecule can be written as

(4) ω = 4 ω 0 E ( 2 IP ) 5 2 exp 2 3 ( 2 IP ) 3 2 E ,

where ω 0 is the atomic unit of frequency (ω 0 = 4 × 1016/s), E is the electric field intensity, and IP is the ionization potential. The t = 1 / ω is the average tunnel ionization time of the molecule. The vibration period can be written as T = 1 / ( c k ) , where c is velocity of light and k is the wavenumber. With the increasing of electric field density, the value of t will decrease. When t < T, the tunneling effect is obvious and accompanied by molecular decomposition. The electric field density corresponding to t = T is the upper limit of electric field density.

In this study, the maximum applied electric field density is 4.0 × 109 V/m (0.0078 a.u.). The IP of the F10 molecule is 5.26 eV (0.1933 a.u.). According to the formula mentioned above, the value of t is 4.39 × 10−10 s. The wavenumber corresponding to the slowest vibration mode is about 198.1/cm. So, the value of T is 1.68 × 10−13  s. Obviously, t > T. This shows that the electric field in our work is not large enough to cause molecular dissociation. In fact, when the electric field density reaches 6.7 × 109 V/m (0.0132 a.u.), the value of t decreases rapidly to 1.66 × 10−13 s, which is t < T. Molecular dissociation may occur under the electric field.

3.4 Excited state

The excited states of molecules are calculated by time-dependent density functional method. Table 4 shows the excitation energy E ex, transition wavelength λ, and oscillator strength f of the F2, F6, and F10 molecules under the typical electric field intensity. The first six excited states are studied in this work. It is found that the excitation energy of each molecule in each excited state decreases with the increasing of electric field intensity, and the transition wavelength increases accordingly. In the absence of electric field, the oscillator strengths of the three molecules from the ground state to the first two singlet excited states are all zero, which belongs to the forbidden transition. When the field intensity is greater than 3.1 × 109 V/m, the oscillator strength is not zero, and the electronic transition can occur from the ground state to the first and second excited states. Hence, the introduction of electric field makes the electron transition more active. In the discharge test, stronger electric field will be accompanied by stronger light effect. It should be noted that the reduction of excitation energy is much smaller than that of IP under the same electric field. In ref. [36], the IP and excitation energy of molecules, including p-diaminobenzene, 2-propanol, and cyclohexane, is compared under electric field, and the conclusion is similar with ours.

Table 4

Excited state of synthesis ester molecules versus applied electric field

E (108 V/m) m = 1 m = 2 m = 3 m = 4 m = 5 m = 6
0 E ex (eV) 5.6571 5.6578 5.6603 5.6605 6.7224 6.7225
λ (nm) 219.17 219.14 219.04 219.03 184.43 184.43
f 0 0 0.0014 0.0014 0.0106 0.0104
3.1 E ex (eV) 5.6542 5.6572 5.6577 5.663 6.7183 6.722
F2 λ (nm) 219.28 219.16 219.14 218.94 184.55 184.44
f 0 0.0014 0.0003 0.0011 0.0105 0.0106
31 E ex (eV) 5.5967 5.6075 5.6175 5.6923 6.374 6.4994
λ (nm) 221.53 221.1 220.71 217.81 194.51 190.76
f 0.0005 0.001 0.0008 0.0008 0.0012 0.0059
0 E ex (eV) 5.6971 5.6978 5.7007 5.7007 6.7458 6.746
λ (nm) 217.63 217.6 217.49 217.49 183.79 183.79
f 0 0 0.0005 0.0005 0.0093 0.0092
3.1 E ex (eV) 5.6947 5.6977 5.6981 5.7019 6.7426 6.7465
F6 λ (nm) 217.72 217.6 217.59 217.44 183.88 183.78
f 0 0.0001 0.0004 0.0004 0.0095 0.009
31 E ex (eV) 5.6294 5.6498 5.6598 5.7231 6.3722 6.5159
λ (nm) 220.24 219.45 219.06 5.7231 194.57 190.28
f 0.0002 0.0004 0.0003 0.0002 0.0014 0.0054
0 E ex (eV) 5.6959 5.6966 5.6996 5.6996 6.7405 6.7408
λ (nm) 217.67 217.65 217.53 217.53 183.94 183.93
f 0 0 0.0004 0.0004 0.0104 0.0103
3.1 E ex (eV) 5.6931 5.6965 5.6973 5.7017 6.7357 6.7398
F10 λ (nm) 217.78 217.65 217.62 217.45 184.07 183.96
f 0 0.0004 0.0001 0.0004 0.0103 0.0103
31 E ex (eV) 5.6319 5.6487 5.6611 5.7185 6.3833 6.5126
λ (nm) 220.14 219.49 219.01 216.81 194.23 190.37
f 0.0002 0.0004 0.0003 0.0002 0.0017 0.0053

4 Conclusion

In this study, DFT and TD-DFT are used to calculate the molecular properties of the typical synthetic ester molecules under the electric field. The results show that the molecular structure changes under the electric field, and the molecular bond lengths are dependent on the electric field. The dipole moment increases with the increasing of electric field intensity and the carbon chain length. The IPs of the F8 and the F10 molecules decrease sharply under the electric field intensities of 3.1 × 109 and 4.0 × 109 V/m, which is the reason for the formation of fast discharge in the synthesis ester. The calculation results for excited states show that the introduction of electric field makes the electron transition more active. The excitation energy decreases with the increasing of electric field intensity. However, the reduction of excitation energy is much smaller than that of IP under the same electric field. The results obtained in the work are helpful to improve the understanding of discharge mechanism in synthetic ester dielectrics and provide theoretical support for the performance improvement of synthetic ester insulating oil.

  1. Funding information: The authors acknowledge the Science and Technology Research Project of Higher Education of Hebei Province (Grant No. QN2019069) and the Research Foundation of Hebei University of Economics and Business (Grant Nos. 2019YB11). The resources on the computation are provided by National Supercomputing Center in Shenzhen, China.

  2. Author contributions: All authors have accepted responsibility for the entire content of this article and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

References

[1] Waara P, Norrby T, Prakash B. Tribochemical wear of rail steels lubricated with synthetic ester-based model lubricants. Tribol Lett. 2004;17(3):561–8.10.1023/B:TRIL.0000044505.42373.0eSearch in Google Scholar

[2] Pompili M, Mazzetti C, Bartnikas R. Comparative PD pulse burst characteristics of transformer type natural and synthetic ester fluids and mineral oils. IEEE Trans Dielectr Electr Insul. 2009;16(6):1511–8.10.1109/TDEI.2009.5361569Search in Google Scholar

[3] Loiselle L, Rao UM, Fofana I. Influence of aging on oil degradation and gassing tendency for mineral oil and synthetic ester under low energy discharge electrical faults. Energies. 2020;13(3):595.10.3390/en13030595Search in Google Scholar

[4] Wang K, Wang F, Lou Z, Han Q, Li J. Relationship between the electrical characteristics of molecules and fast streamers in ester insulation oil. Int J Mol Sci. 2020;21(3):974.10.3390/ijms21030974Search in Google Scholar PubMed PubMed Central

[5] Ye W, Hao J, Zhu M, Li J, Liao R. Low temperature breakdown behavior analysis of natural esters from the perspective of molecular structure. J Mol Liq. 2020;318(15):114032.10.1016/j.molliq.2020.114032Search in Google Scholar

[6] Novick NJ, Mehta PG, Mcgoldrick PB. Assessment of the biodegradability of mineral oil and synthetic ester base stocks, using CO2 ultimate biodegradability tests and CEC-L-33-T-82. J Synth Lubr. 2010;13(1):19–30.10.1002/jsl.3000130103Search in Google Scholar

[7] Beran E. Experience with evaluating biodegradability of lubricating base oils. Tribol Int. 2008;41(12):1212–8.10.1016/j.triboint.2008.03.003Search in Google Scholar

[8] Rozga P, Stanek M. Characteristics of streamers developing at inception voltage in small gaps of natural ester, synthetic ester and mineral oil under lightning impulse. IET Sci Meas Technol. 2016;10(1):50–7.10.1049/iet-smt.2014.0313Search in Google Scholar

[9] Rozga P. The influence of concentrated heat flux on dielectric properties of synthetic and natural esters. In 2012 IEEE International Power Modulator and High Voltage Conference; 2012. p. 378–81.10.1109/IPMHVC.2012.6518759Search in Google Scholar

[10] Przybylek P, Moranda H, Moscicka H, Szczesniak D. Application of synthetic ester for drying distribution transformer insulation—the influence of cellulose thickness on drying efficiency. Energies. 2019;12(1):3874.10.3390/en12203874Search in Google Scholar

[11] Rozga P. Using the three-parameter Weibull distribution in assessment of threshold strength of pressboard impregnated by different liquid dielectrics. IET Sci Meas Technol. 2016;10(6):665–70.10.1049/iet-smt.2016.0061Search in Google Scholar

[12] Dombek G, Gielniak J. Fire safety and electrical properties of mixtures of synthetic ester/mineral oil and synthetic ester/natural ester. IEEE Trans Dielectr Electr Insul. 2018;25(5):1846–52.10.1109/TDEI.2018.007223Search in Google Scholar

[13] Reffas A, Moulai H, Beroual A. Comparison of dielectric properties of olive oil, mineral oil, and other natural and synthetic ester liquids under AC and lightning impulse stresses. IEEE Trans Dielectr Electr Insul. 2018;25(1):1822–30.10.1109/TDEI.2018.007218Search in Google Scholar

[14] Liu Q, Wang ZD. Streamer characteristic and breakdown in synthetic and natural ester transformer liquids under standard lightning impulse voltage. IEEE Trans Dielectr Electr Insul. 2011;18(1):285–94.10.1109/TDEI.2011.5704520Search in Google Scholar

[15] Rozga P. Studies on behavior of dielectric synthetic ester under the influence of concentrated heat flux. IEEE Trans Dielectr Electr Insul. 2016;23(2):908–14.10.1109/TDEI.2015.005135Search in Google Scholar

[16] Li J, Wang Y, Wang F, Liang S, Lin X, Chen X, Zhou J. A study on ionization potential and electron trap of vegetable insulating oil related to streamer inception and propagation. Phys Lett A. 2017;381(44):3732–8.10.1016/j.physleta.2017.09.037Search in Google Scholar

[17] Wang Y, Wang F, Li J, Liang S, Zhou J. Electronic properties of typical molecules and the discharge mechanism of vegetable and mineral insulating oils. Energies. 2018;11(3):523.10.3390/en11030523Search in Google Scholar

[18] Davari N, Astrand PO, Ingebrigtsen S, Unge M. Excitation energies and ionization potentials at high electric fields for molecules relevant for electrically insulating liquids. J Appl Phys. 2013;113(14):143707.10.1063/1.4800118Search in Google Scholar

[19] Tobazcon R. Prebreakdown phenomena in dielectric liquids. IEEE Trans Dielectr Electr Insul. 1994;1(6):1132–47.10.1109/94.368650Search in Google Scholar

[20] Beroual A, Zahn M, Badent A, Kist K, Torshin Y. Propagation and structure of streamers in liquid dielectrics. IEEE Electr Insul Mag. 1998;14(2):6–17.10.1109/57.662781Search in Google Scholar

[21] Rozga P. Streamer propagation in small gaps of synthetic ester and mineral oil under lightning impulse. IEEE Trans Dielectr Electr Insul. 2015;22(5):2754–62.10.1109/TDEI.2015.005101Search in Google Scholar

[22] Buckingham AD. Direct method of measuring molecular quadrupole moments. J Chem Phys. 1959;30(6):1580–5.10.1063/1.1730242Search in Google Scholar

[23] Gisbergen S, Snijders JG, Baerends EJ. Implementation of time-dependent density functional response equations. Comput Phys Commun. 1999;118(2–3):119–38.10.1016/S0010-4655(99)00187-3Search in Google Scholar

[24] Grozema FC, Telesca R, Jonkman HT, Sieb L, Snijders JG. Excited state polarizabilities of conjugated molecules calculated using time dependent density functional theory. J Chem Phys. 2001;115(21):10014–21.10.1063/1.1415085Search in Google Scholar

[25] Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, et al. Wallingford, CT, USA: Gaussian 09. Gaussian Inc; 2013.Search in Google Scholar

[26] The National Institute of Standards and Technology. Computational chemistry comparison and benchmark data base; 2020. http://cccbdb.nist.gov/vibscalejust.asp.Search in Google Scholar

[27] The National Institute of Standards and Technology. Standard reference database 69. NIST Chemistry WebBook; 2018. https://webbook.nist.gov.Search in Google Scholar

[28] Skoog DA, Holler FJ, Crouch SR. Principles of Instrumental Analysis. Boston, MA, USA: Cengage Learning; 2017.Search in Google Scholar

[29] Oommen TV, Claiborne CC, Walsh EJ. Introduction of a new fully biodegradable dielectric fluid. In IEEE Annual Textile, Fiber and Film Industry Technical Conference; 1998. p. 249–52.10.1109/TEXCON.1998.679223Search in Google Scholar

[30] Li J, Zhang Z, Zou P, Grzybowski S, Zahn M. Preparation of a vegetable oil-based nanofluid and investigation of its breakdown and dielectric properties. IEEE Electr Insul Mag. 2012;28(5):43–50.10.1109/MEI.2012.6268441Search in Google Scholar

[31] Greenwood HH. A molecular orbital theory of reactivity in aromatic hydrocarbons. J Phys Chem B. 1952;20(10):1653.10.1063/1.1700243Search in Google Scholar

[32] Li J, Liu X. Density functional theory study of MgnNi2 (n = 1-6) clusters. Chin Phys B. 2012;21(3):033101.10.1088/1674-1056/21/3/033101Search in Google Scholar

[33] Lu T, Chen F. Multiwfn: a multifunctional wavefunction analyzer. J Comput Chem. 2012;33:580–92.10.1002/jcc.22885Search in Google Scholar PubMed

[34] Smalo HS, Astrand PO, Ingebrigtsen S. Calculation of ionization potentials and electron affinities for molecules relevant for streamer initiation and propagation. IEEE Trans Dielectr Electr Insul. 2010;17(3):733–41.10.1109/ICDL.2008.4622536Search in Google Scholar

[35] Jadidian J, Zahn M, Lavesson N, Widlund O, Borg K. Effects of impulse voltage polarity, peak amplitude, and rise time on streamers initiated from a needle electrode in transformer oil. IEEE Trans Plasma Sci. 2012;40(3):909–18.10.1109/TPS.2011.2181961Search in Google Scholar

[36] Smalo HS, Hestad Ø, Ingebrigtsen S, Åstrand P-O. Field dependence on the molecular ionization potential and excitation energies compared to conductivity models for insulation materials at high electric fields. J Appl Phys. 2011;109:073306.10.1063/1.3562139Search in Google Scholar

[37] Wang Y, Wang F, Li J, Huang Z, Liang S, Zhou J. Molecular structure and electronic properties of triolein molecule under an external electric field related to streamer initiation and propagation. Energies. 2017;10:510.10.3390/en10040510Search in Google Scholar

[38] Zuo T, Bandrauk AD. Charge-resonance-enhanced ionization of diatomic molecular ions by intense lasers. Phys Rev A. 1995;52(4):R2511–4.10.1103/PhysRevA.52.R2511Search in Google Scholar

[39] Xiang Q, Liu Y, Zhang X, Duan Y, Xiang M. Spectrum and physical properties of C70 under the external electric field. J Clust Sci. 2020;31(1):951–60.10.1007/s10876-019-01700-xSearch in Google Scholar

Received: 2021-07-02
Revised: 2021-09-29
Accepted: 2021-10-15
Published Online: 2021-11-06

© 2021 Yachao Wang et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Circular Rydberg states of helium atoms or helium-like ions in a high-frequency laser field
  3. Closed-form solutions and conservation laws of a generalized Hirota–Satsuma coupled KdV system of fluid mechanics
  4. W-Chirped optical solitons and modulation instability analysis of Chen–Lee–Liu equation in optical monomode fibres
  5. The problem of a hydrogen atom in a cavity: Oscillator representation solution versus analytic solution
  6. An analytical model for the Maxwell radiation field in an axially symmetric galaxy
  7. Utilization of updated version of heat flux model for the radiative flow of a non-Newtonian material under Joule heating: OHAM application
  8. Verification of the accommodative responses in viewing an on-axis analog reflection hologram
  9. Irreversibility as thermodynamic time
  10. A self-adaptive prescription dose optimization algorithm for radiotherapy
  11. Algebraic computational methods for solving three nonlinear vital models fractional in mathematical physics
  12. The diffusion mechanism of the application of intelligent manufacturing in SMEs model based on cellular automata
  13. Numerical analysis of free convection from a spinning cone with variable wall temperature and pressure work effect using MD-BSQLM
  14. Numerical simulation of hydrodynamic oscillation of side-by-side double-floating-system with a narrow gap in waves
  15. Closed-form solutions for the Schrödinger wave equation with non-solvable potentials: A perturbation approach
  16. Study of dynamic pressure on the packer for deep-water perforation
  17. Ultrafast dephasing in hydrogen-bonded pyridine–water mixtures
  18. Crystallization law of karst water in tunnel drainage system based on DBL theory
  19. Position-dependent finite symmetric mass harmonic like oscillator: Classical and quantum mechanical study
  20. Application of Fibonacci heap to fast marching method
  21. An analytical investigation of the mixed convective Casson fluid flow past a yawed cylinder with heat transfer analysis
  22. Considering the effect of optical attenuation on photon-enhanced thermionic emission converter of the practical structure
  23. Fractal calculation method of friction parameters: Surface morphology and load of galvanized sheet
  24. Charge identification of fragments with the emulsion spectrometer of the FOOT experiment
  25. Quantization of fractional harmonic oscillator using creation and annihilation operators
  26. Scaling law for velocity of domino toppling motion in curved paths
  27. Frequency synchronization detection method based on adaptive frequency standard tracking
  28. Application of common reflection surface (CRS) to velocity variation with azimuth (VVAz) inversion of the relatively narrow azimuth 3D seismic land data
  29. Study on the adaptability of binary flooding in a certain oil field
  30. CompVision: An open-source five-compartmental software for biokinetic simulations
  31. An electrically switchable wideband metamaterial absorber based on graphene at P band
  32. Effect of annealing temperature on the interface state density of n-ZnO nanorod/p-Si heterojunction diodes
  33. A facile fabrication of superhydrophobic and superoleophilic adsorption material 5A zeolite for oil–water separation with potential use in floating oil
  34. Shannon entropy for Feinberg–Horodecki equation and thermal properties of improved Wei potential model
  35. Hopf bifurcation analysis for liquid-filled Gyrostat chaotic system and design of a novel technique to control slosh in spacecrafts
  36. Optical properties of two-dimensional two-electron quantum dot in parabolic confinement
  37. Optical solitons via the collective variable method for the classical and perturbed Chen–Lee–Liu equations
  38. Stratified heat transfer of magneto-tangent hyperbolic bio-nanofluid flow with gyrotactic microorganisms: Keller-Box solution technique
  39. Analysis of the structure and properties of triangular composite light-screen targets
  40. Magnetic charged particles of optical spherical antiferromagnetic model with fractional system
  41. Study on acoustic radiation response characteristics of sound barriers
  42. The tribological properties of single-layer hybrid PTFE/Nomex fabric/phenolic resin composites underwater
  43. Research on maintenance spare parts requirement prediction based on LSTM recurrent neural network
  44. Quantum computing simulation of the hydrogen molecular ground-state energies with limited resources
  45. A DFT study on the molecular properties of synthetic ester under the electric field
  46. Construction of abundant novel analytical solutions of the space–time fractional nonlinear generalized equal width model via Riemann–Liouville derivative with application of mathematical methods
  47. Some common and dynamic properties of logarithmic Pareto distribution with applications
  48. Soliton structures in optical fiber communications with Kundu–Mukherjee–Naskar model
  49. Fractional modeling of COVID-19 epidemic model with harmonic mean type incidence rate
  50. Liquid metal-based metamaterial with high-temperature sensitivity: Design and computational study
  51. Biosynthesis and characterization of Saudi propolis-mediated silver nanoparticles and their biological properties
  52. New trigonometric B-spline approximation for numerical investigation of the regularized long-wave equation
  53. Modal characteristics of harmonic gear transmission flexspline based on orthogonal design method
  54. Revisiting the Reynolds-averaged Navier–Stokes equations
  55. Time-periodic pulse electroosmotic flow of Jeffreys fluids through a microannulus
  56. Exact wave solutions of the nonlinear Rosenau equation using an analytical method
  57. Computational examination of Jeffrey nanofluid through a stretchable surface employing Tiwari and Das model
  58. Numerical analysis of a single-mode microring resonator on a YAG-on-insulator
  59. Review Articles
  60. Double-layer coating using MHD flow of third-grade fluid with Hall current and heat source/sink
  61. Analysis of aeromagnetic filtering techniques in locating the primary target in sedimentary terrain: A review
  62. Rapid Communications
  63. Nonlinear fitting of multi-compartmental data using Hooke and Jeeves direct search method
  64. Effect of buried depth on thermal performance of a vertical U-tube underground heat exchanger
  65. Knocking characteristics of a high pressure direct injection natural gas engine operating in stratified combustion mode
  66. What dominates heat transfer performance of a double-pipe heat exchanger
  67. Special Issue on Future challenges of advanced computational modeling on nonlinear physical phenomena - Part II
  68. Lump, lump-one stripe, multiwave and breather solutions for the Hunter–Saxton equation
  69. New quantum integral inequalities for some new classes of generalized ψ-convex functions and their scope in physical systems
  70. Computational fluid dynamic simulations and heat transfer characteristic comparisons of various arc-baffled channels
  71. Gaussian radial basis functions method for linear and nonlinear convection–diffusion models in physical phenomena
  72. Investigation of interactional phenomena and multi wave solutions of the quantum hydrodynamic Zakharov–Kuznetsov model
  73. On the optical solutions to nonlinear Schrödinger equation with second-order spatiotemporal dispersion
  74. Analysis of couple stress fluid flow with variable viscosity using two homotopy-based methods
  75. Quantum estimates in two variable forms for Simpson-type inequalities considering generalized Ψ-convex functions with applications
  76. Series solution to fractional contact problem using Caputo’s derivative
  77. Solitary wave solutions of the ionic currents along microtubule dynamical equations via analytical mathematical method
  78. Thermo-viscoelastic orthotropic constraint cylindrical cavity with variable thermal properties heated by laser pulse via the MGT thermoelasticity model
  79. Theoretical and experimental clues to a flux of Doppler transformation energies during processes with energy conservation
  80. On solitons: Propagation of shallow water waves for the fifth-order KdV hierarchy integrable equation
  81. Special Issue on Transport phenomena and thermal analysis in micro/nano-scale structure surfaces - Part II
  82. Numerical study on heat transfer and flow characteristics of nanofluids in a circular tube with trapezoid ribs
  83. Experimental and numerical study of heat transfer and flow characteristics with different placement of the multi-deck display cabinet in supermarket
  84. Thermal-hydraulic performance prediction of two new heat exchangers using RBF based on different DOE
  85. Diesel engine waste heat recovery system comprehensive optimization based on system and heat exchanger simulation
  86. Load forecasting of refrigerated display cabinet based on CEEMD–IPSO–LSTM combined model
  87. Investigation on subcooled flow boiling heat transfer characteristics in ICE-like conditions
  88. Research on materials of solar selective absorption coating based on the first principle
  89. Experimental study on enhancement characteristics of steam/nitrogen condensation inside horizontal multi-start helical channels
  90. Special Issue on Novel Numerical and Analytical Techniques for Fractional Nonlinear Schrodinger Type - Part I
  91. Numerical exploration of thin film flow of MHD pseudo-plastic fluid in fractional space: Utilization of fractional calculus approach
  92. A Haar wavelet-based scheme for finding the control parameter in nonlinear inverse heat conduction equation
  93. Stable novel and accurate solitary wave solutions of an integrable equation: Qiao model
  94. Novel soliton solutions to the Atangana–Baleanu fractional system of equations for the ISALWs
  95. On the oscillation of nonlinear delay differential equations and their applications
  96. Abundant stable novel solutions of fractional-order epidemic model along with saturated treatment and disease transmission
  97. Fully Legendre spectral collocation technique for stochastic heat equations
  98. Special Issue on 5th International Conference on Mechanics, Mathematics and Applied Physics (2021)
  99. Residual service life of erbium-modified AM50 magnesium alloy under corrosion and stress environment
  100. Special Issue on Advanced Topics on the Modelling and Assessment of Complicated Physical Phenomena - Part I
  101. Diverse wave propagation in shallow water waves with the Kadomtsev–Petviashvili–Benjamin–Bona–Mahony and Benney–Luke integrable models
  102. Intensification of thermal stratification on dissipative chemically heating fluid with cross-diffusion and magnetic field over a wedge
Downloaded on 16.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2021-0077/html
Scroll to top button