Molecular Dynamics Simulations of Molecular Diffusion Equilibrium and Breakdown Mechanism of Oil-Impregnated Pressboard with Water Impurity

The water molecule migration and aggregation behaviors in oil-impregnated pressboard are investigated by molecular dynamics simulations in combination with Monte Carlo molecular simulation technique. The free energy and phase diagram of H2O-dodecylbenzene (DDB) and H2O-cellulose mixtures are calculated by Monte Carlo technique combined with the modified Flory-Huggins model, demonstrating that H2O molecules can hardly dissolved with infinitesimal content in cellulose system at temperature lower than 650 K, based on which the oil/cellulose layered structure with water impurity representing three-phase coexistence in oil-impregnated pressboard are modeled and performed for molecular dynamics. The molecular dynamics of H2O/DDB/cellulose three-phase mixture simulating oil-paper insulating system with H2O impurity indicates that DDB molecules can thermally intrude into the cellulose-water interface so as to separate the water phase and cellulose fiber. The first-principles electronic structure calculations for local region of H2O/DDB interface show that H2O molecules can introduce bound states to trap electrons and acquire negative charges, so that they will obtain sufficient energy from applied electric field to break DDB molecular chain by collision, which are verified by subsequent molecular dynamics simulations of H2O−/DDB interface model. The electric breakdown mechanism under higher than 100 kV/m electric field is presented based on the further first-principles calculations of the produced carbonized fragments being dissolved and diffusing in DDB phase. The resulted broken DDB fragments will introduce impurity band between valence and conduction bands of DDB system, evidently decreasing bandgap as to that of conducting materials in their existence space. The conductance channel of these carbonized DDB fragments will eventually be formed to initiate the avalanche breakdown process by the cycle-feedback of injected charge carriers with carbonized channels.


Introduction
Synthetic liquid insulating material such as synthetic oil is a mixture of many hydrocarbons, the purification and preparation process of which is complex. Because the general insulating oil is easy to burn and has low heat resistance and dielectric constant, it is necessary to study and develop novel kinds of liquid insulating materials with excellent properties [1]. Therefore, the formation mechanism of dielectric properties of liquid insulating materials should be adequately understood. 3 and 0.5 mol % H2O-cellulose mixtures. Employing the rotational isomeric state (RIS) model [13], the simulated cellulose molecules with tacticity and random torsion are constructed from glucose monomers to be connected by glycosidic bonds in 5-10 polymerization degrees, as shown in Figure  1b, with the consideration that the cellulose molecules are in very low polymerization degree in amorphous regions of insulation pressboard where the H2O molecules can possibly dissolved [14].    [15,16]. The molecular dynamics simulations as implemented in Materials Studio 8.0 Forcite code are performed for H 2 O-DDB, H 2 O-cellulose mixtures to obtain equilibrium structures and for H 2 O/DDB/cellulose three-phase system to investigate the H 2 O transport and aggregation at the interface of mineral oil and cellulose in oil-impregnated pressboard. The modeled H 2 O/DDB/cellulose three-phase system for oil-impregnated pressboard insulation are constituted by packing tool of Amorphous Cell code, as shown in Figure 1c, based on the pre-constructed DDB-cellulose double layer structure with layer thickness ratio of approximate 1:1, and interface thickness of 5Å which is larger and shorter than chemical bonding and van der Waals interaction distances. The Amorphous Cell calculations have performed a geometry optimization, but the constructed initial structures may still contain stress points and density inhomogeneities that need to be relaxed dynamically. The procedure from molecule and mixture model to free energy/phase diagram calculation and molecular dynamics simulation is schematically exhibited in Figure 2. respectively; (c) Schematic microstructure of porous intervals between cellulose fibers with H2O molecule cluster (water impurity) absorbed on them in pressboard.
Molecular interaction binding energy, free energy, and phase diagram of H2O-DDB and H2Ocellulose mixtures are calculated by Monte Carlo molecular simulation technique combined with the modified Flory-Huggins model, as implemented in Blends module of Materials Studio 8.0 package, to investigate the compatibility of H2O in DDB and cellulose [15,16]. The molecular dynamics simulations as implemented in Materials Studio 8.0 Forcite code are performed for H2O-DDB, H2Ocellulose mixtures to obtain equilibrium structures and for H2O/DDB/cellulose three-phase system to investigate the H2O transport and aggregation at the interface of mineral oil and cellulose in oil-impregnated pressboard. The modeled H2O/DDB/cellulose three-phase system for oil-impregnated pressboard insulation are constituted by packing tool of Amorphous Cell code, as shown in Figure 1c, based on the pre-constructed DDB-cellulose double layer structure with layer thickness ratio of approximate 1:1, and interface thickness of 5Å which is larger and shorter than chemical bonding and van der Waals interaction distances. The Amorphous Cell calculations have performed a geometry optimization, but the constructed initial structures may still contain stress points and density inhomogeneities that need to be relaxed dynamically. The procedure from molecule and mixture model to free energy/phase diagram calculation and molecular dynamics simulation is schematically exhibited in Figure 2.   The geometry optimization for atomic structure and MD simulations are carried out with the PCFF (polymer consistent forcefield) atomic force field which has been verified to be valid for a majority of polymers and compounds [17]. The PCFF has been developed from an earlier forcefield CFF91 intending to be utilized for polymers, organic and inorganic materials, demonstrated to be qualified for such as melamine resins, polycarbonates, polysaccharides and about 20 inorganic metals, as well as for carbohydrates, lipids and nucleic acids. PCFF can be effectively used to calculate mechanical property, cohesive energy, compressibility, elastic constants and heat capacity. The geometry optimization for atomic structure and MD simulations are carried out with the PCFF (polymer consistent forcefield) atomic force field which has been verified to be valid for a majority of polymers and compounds [17]. The PCFF has been developed from an earlier forcefield CFF91 intending to be utilized for polymers, organic and inorganic materials, demonstrated to be qualified for such as melamine resins, polycarbonates, polysaccharides and about 20 inorganic metals, as well as for carbohydrates, lipids and nucleic acids. PCFF can be effectively used to calculate mechanical property, cohesive energy, compressibility, elastic constants and heat capacity. The detailed scheme and parameter setup adopted for geometry optimization and MD simulations are listed in the Table 1. The time step of dynamics taken to sufficiently obtain numerical stability depends on the maximum frequency of atomic vibration formalized with dynamics equation and the corresponding integration algorithms. It can be reasonably assumed the velocity and acceleration being constant in 1/15-1/20 of vibration period, which complies to the common rule-of-thumb that the maximum frequency vibration should be sampled between 15 and 20 times through one cycle. The highest vibration frequency in the carbon-based system originates from the stretching C-H bonds of~10 fs. Hence the time step should be set as 0.5-0.7 fs. For the amorphous system consisting of large polymer molecules, the MD simulation of NPT ensemble will approach thermodynamic equilibrium after 10 ps; while for the composite nanolayered structure to adequately stabilize molecule diffusion at interface and relax interfacial structure, the total simulation time is required to be set substantially longer. Therefore, the time step and total simulation time are set as 0.5 fs and 100.0 ps respectively to acquire definite thermodynamic ensemble and accurate dynamics calculations.

Free Energy and Phase Diagram
The best known thermodynamics theory of mixing and phase separation in binary systems is the Flory-Huggins model [15,16], in which the free energy for mixing a binary system is represented as: where ∆G denotes the free energy of mixing (per mole), ϕ i symbolizes the volume fraction of component i, n i represents the degree of polymerization of component i, χ is the interaction parameter, R and T identify gas constant and absolute temperature respectively. The first two terms represent the combinatorial entropy, being always negative and thus preferring a mixed state over pure components, while the last term characterizes the free energy from interaction, deprecating a mixed state if the interaction parameter χ is positive. The balance between the two contributions gives rise to various phase diagrams. Phase diagrams are useful in illustrating the compatibility of binary mixtures. A phase diagram generally contains three pieces of information: critical points, binodals and spinodals. At the critical point, both the second derivative and third derivatives of the free energy with respect to composition vanish. The coexistence region is bound by the binodal (the blue line in Figure 4), in which the mixture can lower its free energy by separating into two phases. The coexistence region also gives rise to two compositions where the second derivative of the free energy is zero. The line through those points is called the spinodal. The spinodal separates the coexistence region into two regions: the mixture is metastable only start phase separating after a sufficiently large fluctuation between the binodal and spinodal. The mixture is unstable in the region bounded by the spinodal, and thus any fluctuation will cause the mixture to spontaneously phase separate. According to Equation (1), the increment of polymerization degree generally decreases miscibility, since the energy term (prefer to phase separation) increases with respect to the entropy term (favoring mixing). Therefore, applying for H 2 O molecules solute in amorphous cellulose system, the H 2 O solubility decreases appreciably with polymerization degree of cellulose chains, as the free energy and phase diagram results shown in Figure 3 from our calculations. These results are inconsistent with the verdict concluded by K. Mazeau that no evident variation has been found for physico-chemical properties of cellulose when length of molecular chain varies [14]. Although H 2 O solubility in DDB and cellulose are pretty small, there are substantial metastable regions between spinodal and binodals in phase diagram presenting significantly higher H 2 O molecule concentrations approaching to 5% in Mole fraction, as illustrated in Figure 3. These results predict there may exist about 0.4 wt. % metastable H 2 O molecules in DDB and amorphous cellulose system at 400 K and 650 K respectively, based on which for the subsequent molecular dynamics simulation, we construct unsaturated mixture model with Mole percentage under these values determined from calculated phase diagram.
Polymers 2018, 10, x FOR PEER REVIEW 6 of 12 6 The coexistence region also gives rise to two compositions where the second derivative of the free energy is zero. The line through those points is called the spinodal. The spinodal separates the coexistence region into two regions: the mixture is metastable only start phase separating after a sufficiently large fluctuation between the binodal and spinodal. The mixture is unstable in the region bounded by the spinodal, and thus any fluctuation will cause the mixture to spontaneously phase separate. According to Equation (1), the increment of polymerization degree generally decreases miscibility, since the energy term (prefer to phase separation) increases with respect to the entropy term (favoring mixing). Therefore, applying for H2O molecules solute in amorphous cellulose system, the H2O solubility decreases appreciably with polymerization degree of cellulose chains, as the free energy and phase diagram results shown in Figure 3 from our calculations. These results are inconsistent with the verdict concluded by K. Mazeau that no evident variation has been found for physico-chemical properties of cellulose when length of molecular chain varies [14]. Although H2O solubility in DDB and cellulose are pretty small, there are substantial metastable regions between spinodal and binodals in phase diagram presenting significantly higher H2O molecule concentrations approaching to 5% in Mole fraction, as illustrated in Figure 3. These results predict there may exist about 0.4 wt. % metastable H2O molecules in DDB and amorphous cellulose system at 400 K and 650 K respectively, based on which for the subsequent molecular dynamics simulation, we construct unsaturated mixture model with Mole percentage under these values determined from calculated phase diagram. Base on the experiments, the substantial water constituents can been absorbed into pressboard, while H2O molecules cannot saluted with nearly neglected content in cellulose system at temperature lower than 650 K from our Monte Carlo calculation results as shown in Figure 3c,d. This is due to that the H2O absorbed in the pressboard actually aggregate in the pores formed between cellulose fibers without permeating into the fibers constructed from cellulose molecules, as illustrated from Figure 1c.   Base on the experiments, the substantial water constituents can been absorbed into pressboard, while H 2 O molecules cannot saluted with nearly neglected content in cellulose system at temperature lower than 650 K from our Monte Carlo calculation results as shown in Figure 3c,d. This is due to that the H 2 O absorbed in the pressboard actually aggregate in the pores formed between cellulose fibers without permeating into the fibers constructed from cellulose molecules, as illustrated from Figure 1c.

Diffusion and Phase Equilibrium
Although H 2 O molecule densities in the DDB and cellulose are appreciably different, when the two equilibrium phases of DDB and cellulose are contacted to each other, H 2 O molecules in cellulose with much lower density could diffuse towards DDB side because that it is the interactions O migration at the DDB/cellulose interface is interface free energy which in thermodynamics identify the interaction between DDB and cellulose molecules at interface. It is also definitely noted from the molecular dynamics simulation that the H 2 O molecules cannot dissolve or permeate into the cellulose phase so that the H 2 O molecules have been excluded from the cellulose system in the completed modeling of H 2 O/DDB/cellulose three-phase coexistence.
The H 2 O/DDB/cellulose three-phase mixture system is firstly constructed by multi-molecules Monte Carlo method of packing technique as implemented by Amorphous Cell program, in which the atomic geometry has also been optimized. Then the initially modeled configuration is thermodynamically equilibrated by performing molecular dynamics simulation to represent oil-paper insulating system containing H 2 O impurity at ambient temperature and pressure. In practical transformer insulation, before pressboard being impregnated in DDB to form oil-paper complex dielectrics, the condensed H 2 O cluster have absorbed on cellulose fibers in porous intervals of pressboard even after heating treatment in vacuum drying oven. Therefore, in the practical process of preparing oil (DDB)-impregnated pressboard, the DDB permeate into the micro-pores as interval space between cellulose fibers of pressboard, and clad cellulose fibers together with H 2 O condensed cluster absorbed on them, as initial three-phase model shown in left Figure 4. Because of the relatively better miscibility of DDB-cellulose than that of cellulose-water, under the thermodynamic driving forces from interface free energies (interface tensions) of three coupling interfaces, the cellulose-water interface spontaneously shrinks and finally disappears, which exactly means the DDB molecules thermally diffuse as intruding into the cellulose-water interface and eventually separate the water phase and cellulose fiber, as shown in right Figure 4 of phase equilibrium after thermodynamics simulation.

Diffusion and Phase Equilibrium
Although H2O molecule densities in the DDB and cellulose are appreciably different, when the two equilibrium phases of DDB and cellulose are contacted to each other, H2O molecules in cellulose with much lower density could diffuse towards DDB side because that it is the interactions from the based solvent of DDB or cellulose, not the random collisions between H2O molecules, determining the molecule migration of solute H2O. The fugacity of solute component in solution or mixture characterizes the solute-solvent molecular interaction and thus the pressure deriving from molecule diffusion of solute, therefore the H2O diffusion trend at DDB/cellulose interface lies thermodynamically on the individual H2O fugacities in DDB and cellulose. Another thermodynamic element that may dominate the H2O migration at the DDB/cellulose interface is interface free energy which in thermodynamics identify the interaction between DDB and cellulose molecules at interface. It is also definitely noted from the molecular dynamics simulation that the H2O molecules cannot dissolve or permeate into the cellulose phase so that the H2O molecules have been excluded from the cellulose system in the completed modeling of H2O/DDB/cellulose three-phase coexistence.
The H2O/DDB/cellulose three-phase mixture system is firstly constructed by multi-molecules Monte Carlo method of packing technique as implemented by Amorphous Cell program, in which the atomic geometry has also been optimized. Then the initially modeled configuration is thermodynamically equilibrated by performing molecular dynamics simulation to represent oil-paper insulating system containing H2O impurity at ambient temperature and pressure. In practical transformer insulation, before pressboard being impregnated in DDB to form oil-paper complex dielectrics, the condensed H2O cluster have absorbed on cellulose fibers in porous intervals of pressboard even after heating treatment in vacuum drying oven. Therefore, in the practical process of preparing oil (DDB)-impregnated pressboard, the DDB permeate into the micro-pores as interval space between cellulose fibers of pressboard, and clad cellulose fibers together with H2O condensed cluster absorbed on them, as initial three-phase model shown in left Figure 4. Because of the relatively better miscibility of DDB-cellulose than that of cellulose-water, under the thermodynamic driving forces from interface free energies (interface tensions) of three coupling interfaces, the cellulose-water interface spontaneously shrinks and finally disappears, which exactly means the DDB molecules thermally diffuse as intruding into the cellulose-water interface and eventually separate the water phase and cellulose fiber, as shown in right Figure 4 of phase equilibrium after thermodynamics simulation.

Electric Breakdown Process
In order to determine the charge distribution and the motion of H 2 O molecule under the electric field, the electronic structure and charge distribution of the simplified model of H 2 O/DDB mixture (including six DDB molecules and one water molecule, representing the local region of the H 2 O/DDB two-phase interface) under different temperatures and electric fields are calculated by all-electron numerical orbital first-principles method with PBEsol form of GGA correlation-exchanged functional as implemented in DMol3 code [20]. The localized trap state introduced by H 2 O molecule and its electrical transport behavior in DDB are analyzed based on the energy density of states (DOS), as the calculated results shown in Figure 5. Calculation of Mulliken population charge of H 2 O (sum of Mulliken charges of hydrogen and oxygen atoms) is calculated from self-consistent charge density. The calculated results of DOS indicate that the unoccupied bound state (the wave function is localized near the H 2 O molecule) introduced by H 2 O in the band gap of DDB is 0.5 eV lower than the intrinsic electron trap level of DDB matrix which is about 1.0 eV lower than conduction band minimum, as shown in Figure 5b. This bound state can also be looked as the deeper trap state for the electrons of conduction band, which is energetically preferable for electron occupation and will move together with H 2 O thermal diffusion. Meanwhile, the water molecule and the DDB matrix molecule interact by van der Waals force without electronic bonding. Therefore, it is necessary to leap over the energy barrier between H 2 O and DDB molecules which turns to be lowest in the thermal collision between H 2 O and DDB molecules at H 2 O/DDB interface. Consequently, the electrons in conduction band or intrinsic trap states of DDB phase will trap into the unoccupied bound state of H 2 O molecule so as to make H 2 O being loaded with negative charge after the molecular collision, and more hole will be produced simultaneously. The frequency of collision between H 2 O and DDB due to thermal motion increases with temperature increment. Thus, the probability of producing negatively charged H 2 O ion from electron transfer of trapping transition will be raised when temperature goes up. Population analysis of Mulliken atomic charge confirms that the charge on H 2 O is negative, hence the H 2 O − ion produced by thermal collision and electron transfer will transport in the electric field as net motion in the inverse direction of electric field.

Electric Breakdown Process
In order to determine the charge distribution and the motion of H2O molecule under the electric field, the electronic structure and charge distribution of the simplified model of H2O/DDB mixture (including six DDB molecules and one water molecule, representing the local region of the H2O/DDB two-phase interface) under different temperatures and electric fields are calculated by all-electron numerical orbital first-principles method with PBEsol form of GGA correlation-exchanged functional as implemented in DMol3 code [20]. The localized trap state introduced by H2O molecule and its electrical transport behavior in DDB are analyzed based on the energy density of states (DOS), as the calculated results shown in Figure 5. Calculation of Mulliken population charge of H2O (sum of Mulliken charges of hydrogen and oxygen atoms) is calculated from self-consistent charge density. The calculated results of DOS indicate that the unoccupied bound state (the wave function is localized near the H2O molecule) introduced by H2O in the band gap of DDB is 0.5 eV lower than the intrinsic electron trap level of DDB matrix which is about 1.0 eV lower than conduction band minimum, as shown in Figure 5b. This bound state can also be looked as the deeper trap state for the electrons of conduction band, which is energetically preferable for electron occupation and will move together with H2O thermal diffusion. Meanwhile, the water molecule and the DDB matrix molecule interact by van der Waals force without electronic bonding. Therefore, it is necessary to leap over the energy barrier between H2O and DDB molecules which turns to be lowest in the thermal collision between H2O and DDB molecules at H2O/DDB interface. Consequently, the electrons in conduction band or intrinsic trap states of DDB phase will trap into the unoccupied bound state of H2O molecule so as to make H2O being loaded with negative charge after the molecular collision, and more hole will be produced simultaneously. The frequency of collision between H2O and DDB due to thermal motion increases with temperature increment. Thus, the probability of producing negatively charged H2O ion from electron transfer of trapping transition will be raised when temperature goes up. Population analysis of Mulliken atomic charge confirms that the charge on H2O is negative, hence the H2O − ion produced by thermal collision and electron transfer will transport in the electric field as net motion in the inverse direction of electric field. The molecular dynamics simulation utilizes the semi-empirical force field with assigning immutable charge to individual atoms, which cannot reflect the electron transfer process as calculated by the first principles kinetics. However, the first principles calculation is limited by the computational cost, which cannot be used to calculate the whole multi-phases model. In actual electric breakdown experiments, the electrodes are applied in contact with the DDB oil, and the electrons transferred to H2O molecule actually come from the electrode-injected charge carriers rather than a small number of thermally excited electrons. Therefore, in the molecular dynamics simulation, the assigning first-principles calculated charges to H2O molecules can truly represent the real situation of the charge injection by the external electrodes in experiments. Accordingly, The H2O/DDB two-layer model is constructed to simulate the interface of water/oil two-phase coexistence system, in which H2O molecules on H2O/DDB two-phase interface are assigned by charge equivalent to first principles calculated Mulliken charge (the positive charge of the matrix is negligible), and the molecular dynamics simulations with ReaxFF 6.0 force field are performed at ambient temperature and under high electric field with the direction along two-phase interface, as implemented by General Utility Lattice Program (GULP) of materials studio 8.0 package [21].
The molecular dynamics simulations of H2O − /DDB two-phase system under high electric field indicate that when the electric field is high enough (>100 kV/m), the ambient temperature collisions between H2O with DDB molecules at two-phase interface of oil-water in transformer insulating oil will lead to the breaking of the DDB molecular chain and the formation of CH, CH3, C3H5, C7H6, and other unsaturated carbonic fragments which will diffuse into the DDB liquid phase. Figure 6 exhibits structural result of molecular dynamics simulation for H2O − /DDB two-phase system at 300 K temperature and 200 kV/m electric field, in comparison with the geometrically optimized structure without atomic kinetics and electric field (before thermodynamic equilibrium). According to the first principles electronic structure calculations, the molecular orbital energy levels of these fracture fragments are just located in the DDB band-gap, as shown in Figure 7. Therefore, under the extremely high electric field, H2O − ions can obtain high enough kinetic energy as to break C-C and C-H bonds of DDB molecular chain when they collide with DDB molecules. The resulted broken DDB fragments as described by carbonized DDB constituents will introduce impurity band between valence and conduction bands of DDB matrix, and thus cause bandgap to decrease lower than kT (k is Boltzmann constant, T denotes temperature) as that of conducting materials in their existence The molecular dynamics simulation utilizes the semi-empirical force field with assigning immutable charge to individual atoms, which cannot reflect the electron transfer process as calculated by the first principles kinetics. However, the first principles calculation is limited by the computational cost, which cannot be used to calculate the whole multi-phases model. In actual electric breakdown experiments, the electrodes are applied in contact with the DDB oil, and the electrons transferred to H 2 O molecule actually come from the electrode-injected charge carriers rather than a small number of thermally excited electrons. Therefore, in the molecular dynamics simulation, the assigning first-principles calculated charges to H 2 O molecules can truly represent the real situation of the charge injection by the external electrodes in experiments. Accordingly, The H 2 O/DDB two-layer model is constructed to simulate the interface of water/oil two-phase coexistence system, in which H 2 O molecules on H 2 O/DDB two-phase interface are assigned by charge equivalent to first principles calculated Mulliken charge (the positive charge of the matrix is negligible), and the molecular dynamics simulations with ReaxFF 6.0 force field are performed at ambient temperature and under high electric field with the direction along two-phase interface, as implemented by General Utility Lattice Program (GULP) of materials studio 8.0 package [21].
The molecular dynamics simulations of H 2 O − /DDB two-phase system under high electric field indicate that when the electric field is high enough (>100 kV/m), the ambient temperature collisions between H 2 O with DDB molecules at two-phase interface of oil-water in transformer insulating oil will lead to the breaking of the DDB molecular chain and the formation of CH, CH 3 , C 3 H 5 , C 7 H 6 , and other unsaturated carbonic fragments which will diffuse into the DDB liquid phase. Figure 6 exhibits structural result of molecular dynamics simulation for H 2 O − /DDB two-phase system at 300 K temperature and 200 kV/m electric field, in comparison with the geometrically optimized structure without atomic kinetics and electric field (before thermodynamic equilibrium). According to the first principles electronic structure calculations, the molecular orbital energy levels of these fracture fragments are just located in the DDB band-gap, as shown in Figure 7. Therefore, under the extremely high electric field, H 2 O − ions can obtain high enough kinetic energy as to break C-C and C-H bonds of DDB molecular chain when they collide with DDB molecules. The resulted broken DDB fragments as described by carbonized DDB constituents will introduce impurity band between valence and conduction bands of DDB matrix, and thus cause bandgap to decrease lower than kT (k is Boltzmann constant, T denotes temperature) as that of conducting materials in their existence space. With these carbonized DDB constituents being continuously produced at the water/oil interface under high electric field and diffusing through the DDB oil, the conductance channels constituted by high density of them will been constructed. In the actual electric breakdown experiment, the electrodes will be applied at two terminals of oil-impregnated pressboard in direct contact with DDB oil, and after the conductance channel penetrates the oil, the avalanche breakdown process will be initiated: a large number of carriers are injected from both electrodes and obtain sufficient kinetic energy under the action of high electric field so as to break C-C or C-H bonds of DDB molecules and sharply increase the density and region of carbonized DDB constituents, which as a feedback aggravates the carrier injections by conductance channel expansion to the surrounding area. The formation of "avalanche" similar to the discharge process results in the ultimate damage of DDB dielectrics.
Polymers 2018, 10, x FOR PEER REVIEW 10 of 12 space. With these carbonized DDB constituents being continuously produced at the water/oil interface under high electric field and diffusing through the DDB oil, the conductance channels constituted by high density of them will been constructed. In the actual electric breakdown experiment, the electrodes will be applied at two terminals of oil-impregnated pressboard in direct contact with DDB oil, and after the conductance channel penetrates the oil, the avalanche breakdown process will be initiated: a large number of carriers are injected from both electrodes and obtain sufficient kinetic energy under the action of high electric field so as to break C-C or C-H bonds of DDB molecules and sharply increase the density and region of carbonized DDB constituents, which as a feedback aggravates the carrier injections by conductance channel expansion to the surrounding area. The formation of "avalanche" similar to the discharge process results in the ultimate damage of DDB dielectrics.

Conclusions
Employing molecular dynamics simulations in combination with Monte Carlo technique and first-principles calculations, the H2O molecule diffusion and aggregation at interfaces of H2O/DDB/cellulose three-phase coexistence in oil-impregnated pressboard with water impurity are space. With these carbonized DDB constituents being continuously produced at the water/oil interface under high electric field and diffusing through the DDB oil, the conductance channels constituted by high density of them will been constructed. In the actual electric breakdown experiment, the electrodes will be applied at two terminals of oil-impregnated pressboard in direct contact with DDB oil, and after the conductance channel penetrates the oil, the avalanche breakdown process will be initiated: a large number of carriers are injected from both electrodes and obtain sufficient kinetic energy under the action of high electric field so as to break C-C or C-H bonds of DDB molecules and sharply increase the density and region of carbonized DDB constituents, which as a feedback aggravates the carrier injections by conductance channel expansion to the surrounding area. The formation of "avalanche" similar to the discharge process results in the ultimate damage of DDB dielectrics.  Electronic density of states for mixture system containing DDB molecules and CH, CH3, C3H5, and C7H6 unsaturated broken fragments of DDB molecular chain, with the highest occupied state as reference energy zero, Gauss smearing is set as 0.1 eV.

Conclusions
Employing molecular dynamics simulations in combination with Monte Carlo technique and first-principles calculations, the H2O molecule diffusion and aggregation at interfaces of H2O/DDB/cellulose three-phase coexistence in oil-impregnated pressboard with water impurity are  . Electronic density of states for mixture system containing DDB molecules and CH, CH 3 , C 3 H 5 , and C 7 H 6 unsaturated broken fragments of DDB molecular chain, with the highest occupied state as reference energy zero, Gauss smearing is set as 0.1 eV.

Conclusions
Employing molecular dynamics simulations in combination with Monte Carlo technique and first-principles calculations, the H 2 O molecule diffusion and aggregation at interfaces of H 2 O/DDB/cellulose three-phase coexistence in oil-impregnated pressboard with water impurity are simulated to investigate electric breakdown mechanism of oil-paper composite insulation utilized in high voltage transformers. The phase diagrams of H 2 O-DDB and H 2 O-cellulose mixtures from Monte Carlo molecular simulations indicate that H 2 O molecules cannot dissolved into cellulose system at temperature lower than 650 K, attributed to which the layered three-phase coexistence system of oil/cellulose interface structure with water impurity are modeled with molecular dynamics. The molecular dynamics of H 2 O/DDB/cellulose three-phase dominating the breakdown process of oil-paper insulating system containing H 2 O impurity demonstrate that DDB molecules will diffuse thermally through cellulose-water interface, resulting in the separation of water phase and cellulose fiber. The subsequent first-principles electronic structure calculations for H 2 O/DDB interface region prove that H 2 O molecules can acquire negative charges by trapping electrons into their bound states through molecular collisions and will obtain adequate energy from electric field to fracture DDB molecular chain by collision, which are confirmed by correlated molecular dynamics simulations of H 2 O − /DDB interface model. The electric breakdown mechanism under higher than 100 kV/m electric field is revealed theoretically based on the further first-principles calculation of the carbonized fragments resulted from breaking DDB molecules by high-electric-field-triggered violent H 2 O − collisions. The produced carbonized DDB fractures will introduce assistant band between valence and conduction bands of DDB matrix, which decreases bandgap in their existence space to the value as that of conducting material. It is reasonably predicted that the conductance channel constituted by these continuously condensed and extended carbonized regions will be promptly formed so as to establish the avalanche breakdown process through cycle-feedback of injected charge carriers with carbonized channels under high voltage.