Next Article in Journal
Fabrication of Al2O3 Nano-Structure Functional Film on a Cellulose Insulation Polymer Surface and Its Space Charge Suppression Effect
Previous Article in Journal
Influence of Polymer Relaxation Time on the Electrospinning Process: Numerical Investigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interfacial Hydrogen Bonds and Their Influence Mechanism on Increasing the Thermal Stability of Nano-SiO2-Modified Meta-Aramid Fibres

1
College of Engineering and Technology, Southwest University, Chongqing 400715, China
2
School of Electronics and Computer Science, University of Southampton, Southampton SO171BJ, UK
3
Laboratory of Power Transmission Equipment & System Security and New Technology, Chongqing University, Chongqing 400044, China
*
Author to whom correspondence should be addressed.
Polymers 2017, 9(10), 504; https://doi.org/10.3390/polym9100504
Submission received: 19 September 2017 / Revised: 9 October 2017 / Accepted: 9 October 2017 / Published: 12 October 2017

Abstract

:
For further analysis of the effect of nano-doping on the properties of high polymers and research into the mechanism behind modified interfacial hydrogen bonds, a study on the formation probability of nano-SiO2/meta-aramid fibre interfacial hydrogen bonds and the strengthening mechanism behind interfacial hydrogen bonds on the thermal stability of meta-aramid fibres using molecular dynamics is performed in this paper. First, the pure meta-aramid fibre and nano-SiO2/meta-aramid fibre mixed models with nanoparticle radiuses of 3, 5, 7 and 9 Å (1 Å = 10−1 nm) are built, and then the optimization process and dynamics simulation of the models are conducted. The dynamics simulation results indicate that the number of hydrogen bonds increase due to the doping by nano-SiO2 and that the number of interfacial hydrogen bonds increases with the nanoparticle radius. By analysing the hydrogen bond formation probability of all the atom pairs in the mixed model with pair correlation functions (PCFs), it can be observed that the hydrogen bond formation probability between the oxygen atom and hydrogen atom on the nanoparticle surface is the greatest. An effective way to increase the number of interfacial hydrogen bonds in nano-SiO2 and meta-aramid fibres is to increase the number of hydrogen atoms on the nano-silica surface and oxygen atoms in the meta-aramid fibre. By using the radial distribution function (RDF), the conclusion can be further drawn that the hydrogen bond formation probability is at a maximum when the atomic distance is 2.7–2.8 Å; therefore, increasing the number of atoms within this range can significantly increase the formation probability of hydrogen bonds. According to the results of chain movement, the existence of interfacial hydrogen bonds effectively limits the free movement of the molecular chains of meta-aramid fibres and enhances the thermal stability of meta-aramid fibres. The existence of interfacial hydrogen bonds is one of the important reasons for formation of the stable interface structure between nanoparticles and meta-aramid fibres. In addition, a nanoparticle with a small radius improves the interfacial hydrogen bond energy density and interfacial interaction energy density, enhancing the stability of the mixed model interface.

1. Introduction

Through the doping of nanoparticles for modified nano materials, the raw materials are provided with many new characteristics, and many properties of the raw materials are improved [1,2]. With the rapid development of nanotechnology, the application field of nano materials has been further expanded [3]. At present, there is a trend to incorporate nanoparticles in high polymer insulating materials to improve the properties of the materials, especially for the improvement of the voltage level in power systems, and this process requires further improvements in the insulating property, thermal stability, mechanical property, and so on, of the electrical insulating material [4]. As an important electrical insulating material, the aramid insulating paper will gradually be widely used [5,6,7,8], and on the study of the aspect of property improvement, modification by doping with nano-silica will gradually become a key research direction widely assumed by many scholars [9,10,11,12].
At present, the focus of the study of modified nanomaterials has gradually been transferred from the improvement of the macro-properties to the exploration of the microscopic mechanism, and, for example, in the documents of [13,14,15,16], and so on, the mechanism behind the improvement to the macro-properties by the nanoparticles is explained through an exploration of the microscopic mechanism. With the exploration of the microcosms, the microscopic mechanism behind the effect of nanomaterials on polymer properties is gradually revealed further. For example, in the research field of biotoxicology of nanomaterials, study of the interaction between nanomaterials and the interface of the biological system [17] is mainly focused on the dynamic process of physical and chemical interactions between the surfaces of nanomaterials, and the biological composition (such as protein, membrane, organelle, DNA, etc.) is included in the “nano-biology” interface; these processes include the exchange reaction kinetics and thermodynamic levels between the nanomaterials and various biomolecules in the biological system [18]. It is widely believed that the DNA double helix owes its stability to the hydrogen bonds holding the two strands together. However, in the dissociated state, the individual strands can form even stronger hydrogen bonds with the solvent water molecules [19]. Klaus Lackner’s group studied hydrogen-bond effects on the ions on the surface of polymer [20] and the hydrogen-bond can affect the interactions of different anions on polymers [21]. Qibin Li’s group found that the hydrogen bonds play an important role in the surface of hydrate [22,23]. The number of average hydrogen bonds per water molecule is little different in the mixture region and the hydrate region, close to the value of the saturation number of water molecules. The radial distribution function (RDF) of H2O particles and F3 order parameter indicate that the water is a crystal structure in the whole system [22]. The phase transition in the water/graphene-oxide system needs more time to break the network of hydrogen bonds among water molecules, thereby reducing the evaporation rate [23]. It can be seen that the hydrogen bond exists widely, and it has certain effects on the structure and properties of materials. Therefore, the hydrogen bond is an important non-negligible factor affecting the interface interactions [24].
The hydrogen bond is a weak force that widely exists between molecules or within the molecules and is similar to the electrostatic interaction. The hydrogen bond that forms between two molecules is called the intermolecular hydrogen bond. The hydrogen bond that forms between atoms in the same molecule is called the intramolecular hydrogen bond. A close relation [13] exists between the hydrogen bonding network of polymers and the melting point, boiling point, structure, cohesive energy density, and so on, of polymers, and a large effect [25] of hydrogen bonds on the mechanical properties and anti-aging properties of cellulose, and so on, has been considered. The intramolecular hydrogen bond in the body plays a pivotal role, for example, a close relation [26] exists between the helical structure of protein, activity of some biochemical molecules, and so on, and hydrogen bonding. In the preparation of macromolecular polymers, the purpose of changing the performance is achieved by mixing various polymers and changing the blend compositions. However, the number of discovered miscible polymer blends is small, and most polymers are not miscible with other polymers; their interface is clear, the bonding between compositions is not strong, and it can be seen in a study that the compatibility of immiscible copolymers can be strengthened by forming intermolecular hydrogen bonds to improve the properties of the polymer blend [27,28]. Furthermore, because a large number of hydroxide radicals exist on the surface of nano-SiO2, a coacervate with a dimensional network structure forms through hydrogen bond interactions among the nanoparticles. To decrease such an agglomeration phenomenon, the surface modification [3] of nanoparticles is carried out, and one of the purposes is to decrease the possibility of the occurrence of the agglomeration phenomenon by decreasing the number of hydrogen bonds at the nanoparticle interface. Xi Chen’s group studied different number of hydrogen bonds can affect the energy barriers of chemical reactions [29]. Hossein Eslami’s group found that: Addition of water molecules to PA-6,6 + water system, confined in wide pores, accumulates water molecules into clusters. In such wide pores, the amide groups are not as accessible to water molecules, therefore hydrogen bonding between the water polymer cannot stabilize the mixture [30]. Studies in the literature [31] show that, The number of hydrogen bonds per donor for free polymer chains decreases with increasing the grafting density, with a faster decrease for systems containing bigger nanoparticles compared to those with smaller nanoparticles. The maximum probability occurs at a distance about 0.15 nm from the surface. Further reports in the literature [32,33] show that: hydrogen-bond formation depends on the layering effect and on the geometrical restrictions and reveals an oscillatory behavior like the solvation force oscillations. The hydrogen bonding in PA-6,6 is responsible for the formation of a 3-D network, giving rise to the special mechanical properties of this polymer. Although the abovementioned results suggest the importance of interfacial hydrogen bonds and that the number of interfacial hydrogen bonds is an important factor affecting interaction, a further study on the important factors affecting the number of interatomic hydrogen bonds on the atomic scale has not been carried out, namely, a study on the formation probability of interatomic hydrogen bonding.
Additionally, the performance of nanocomposites is closely associated with the particle size and morphology of the materials, and the many properties of the composite, such as optical property, electricity property, magnetic property, and so on, can only be viewed at the macroscopic level when the particle size and shape of the materials are uniform [34,35]; thus, control of nanoparticle size has become a popular topic of research. For example, studies of the effects of nanoparticles with different sizes on the properties and the mechanisms behind the modified nanomaterials have been performed [36,37,38,39,40]. Although the studies carried out by many scholars at home and abroad show that nanoscale is one of the most important parameters for characterization of nanomaterials and that the effects of nanoparticles with different sizes on material properties and in the aspect of biological effects are different [41,42,43,44,45,46], further exploration is still required on the microcosmic mechanism behind the effect of the size of the nanoparticle on the properties of the modified materials.
On that account, this paper describes the doping modification of meta-aramid fibres by nano-silica, and studies the hydrogen bond interactions in the SiO2/meta-aramid fibre materials at the atomic scale. First, the pure meta-aramid fibre and mixed SiO2/meta-aramid fibre (with nanoparticle radiuses of 3, 5, 7 and 9 Å) models are built with the Materials Studio (MS) software, and then the optimization processing and dynamics simulation of the model are carried out. In this paper, the effect of the nanoparticle radius on the number of different types of hydrogen bonds is analysed by the statistics of the number of different types of hydrogen bonds in all the models. The formation probability of the hydrogen bonds among different atomic and SiO2/meta-aramid fibre interfaces is reflected through the pair correlation function and radial distribution function The effect of the nano-silica particle size on the thermal stability of the meta-aramid fibres is represented by using the mean square displacement (MSD) of meta-aramid fibres in different models. The effect of the interfacial hydrogen bond interactions on the overall interaction between the nanoparticle and meta-aramid fibre is described through changes in the interaction energy and hydrogen bond energy of the interface.

2. Formation, Model Building and Parameter Setting of Hydrogen Bonds

2.1. Formation of Hydrogen Bonds

The hydrogen bond generally refers to the bond formed through the certain bonding force generated between the H atom combined with a high-electronegativity element in a molecule and a high-electronegativity atom in another molecule. That is, D–H…A by formula, where D refers to the donor atom (donor), A refers to the acceptor atom (acceptor), and both of them represent atoms with great electronegativity and a small radius, such as F, O, N, and so on.
In D–H…A, D–H refers to the strong polar covalent bond, and because the electronegativity of D is great and the electronic attraction is strong, the hydrogen atom becomes a “naked” electron nearly without an electron cloud and with some positive charge. The radius of D is especially small, the electric field intensity of D is great, and no inner electron exists for D. The other atom, A (that is, the atom with great electronegativity, small radius and lone pair electrons), with some negative charge, is allowed to sufficiently reach D; therefore, a strong electrostatic interaction is generated, and the hydrogen bond is formed. Hydrogen bonding can be in the range of 2–8 kcal, and compared with the intermolecular action force (i.e., van der Waals force), the bonding is slightly strong; however, compared with the covalent bond and ionic bond, the bonding is weak. The hydrogen bond energy in the range of 25–40 kJ/mol belongs to a medium-strength hydrogen bond, and it is generally acknowledged that a bond is of weak hydrogen bond strength when the bond energy is less than 25 kJ/mol, whereas the bond is of a strong hydrogen bond when the bond energy is greater than 40 kJ/mol.
In classical molecular dynamics simulation, the definition of the hydrogen bond is vague [47,48]; therefore, two methods are usually adopted: energy criterion and the geometry rule. The definition of the hydrogen bond by the energy criterion is when the interaction energy between molecular pairs is greater than a certain value, and this interaction between the molecular pairs is deemed a hydrogen bond interaction [49]. The definition of hydrogen bond by the geometry rule is the relative position between two molecules [48,50,51], as shown in Figure 1; the maximum distance (rAH) between the acceptor atom and hydrogen atom is defined as 3 Å (angstrom, 1 Å = 10−1 nm), the minimum value of the included angle (β) is 90°, and the donor atom and acceptor atom are described using the oxygen atom as an example. In this paper, the definition of the hydrogen bond is made by the geometry rule [13].

2.2. Model Building

First, the meta-aramid fibre single-chain model is built, and the polymerization degree (DP) is 10, as shown in Figure 2. The left half of the figure indicates the repeating unit structure in the meta-aramid chain and a single meta-aramid molecular model, and the right half indicates the bond angle among atoms N, H and O and bond length of the atoms H and O.
By using the Visualizer module in MS, an amorphous model can be built through the Amorphous Cell (AC) module, for example, Model A in Figure 3 is a pure meta-aramid fibre model. A periodic boundary condition is established. Then, the silica unit cell is imported to build the nano-silica cluster model. Because unsaturated residual bonds and hydroxide radicals with different bonds and status [52,53] exist on the surface of nano-SiO2, during the building of the nano-SiO2 cluster model, the position supplement is conducted with the hydrogen atom for the unsaturated chemical bond existing on the surface because of shearing. The mixed model of the SiO2/meta-aramid fibre is built with the Packing tool in the AC module. For the model building, the periodic boundary condition is established, the density is 1.2 g/cc, the chain number (N aramid) of the meta-aramid fibres is 4, the number of nano-SiO2 is 1, and the radiuses of nanoparticles are 3, 5, 7 and 9 Å. The figure for the building of the mixed models of nano-SiO2/meta-aramid fibre (Model B, Model C, Model D and Model E) is shown in Figure 3.

2.3. Simulation Parameter Setting

During the geometric optimization, annealing treatment and molecular dynamics simulation of the model, the COMPASS (condensed-phase optimized molecular potentials for atomistic simulation studies) force field applicable to the organic and inorganic molecules [44] is adopted as the force field. After the model is built, the energy is great and the model is unstable; therefore, a geometric optimization is carried out for all the models, and then an annealing treatment is conducted with the number of annealing cycles being 10, an initial temperature of 300 K, and the middle temperature of each cycle set as 900 K. After the complete annealing treatment, the model energy reaches a reasonable value, and the dynamics simulation is subsequently conducted. A 300 ps simulation for the model after relaxation is first carried out with the Discover module in the NVT (constant volume and temperature) ensemble, the simulation time step is 1 fs, the cut-off radius is set as 9.5 Å, and the spline width is set as 1 Å; subsequently, the NPT (constant pressure and temperature) ensemble is selected, the Andersen method is adopted as the temperature control method, the pressure intensity is set as the standard atmospheric pressure, and the information of the primary dynamics is collected with the Berendsen pressure controlling method every 50 steps.

3. Simulation Results and Discussion

3.1. Change in Number of Hydrogen Bonds in all the Models

In the mixed model of the nano-SiO2/meta-aramid fibre, because the –NH radical and oxygen atom exist in the meta-aramid insulating paper fibre, the –NH…O and –NH…N hydrogen bonds are easily formed in the molecules and between the molecules of the meta-aramid fibre. –OH…O hydrogen bonds are easily formed on the surface of nano-silica. –OH…N, –OH…O and –NH…O hydrogen bonds are easily formed at the meta-aramid insulating paper fibre and silica nanoparticle interface. For a simplified description and instruction, the surficial hydrogen atom on the nano-silica, oxygen atom on the nano-silica, hydrogen atom in the meta-aramid fibre chain–NH radical, nitrogen atom in the meta-aramid fibre chain and oxygen atom in the fibre chain of the meta-aramid insulating paper are respectively marked as atoms ①, ②, ③, ④ and ⑤. Therefore, only ①…②, ①…④, ①…⑤, ②…③, ③…④ and ③…⑤ hydrogen bonds are formed in the model. Therein, ①…② hydrogen bonds are the hydrogen bonds on the nanoparticle; ③…④ and ③…⑤ hydrogen bonds are the hydrogen bonds in the meta-aramid fibre; and ①…④, ①…⑤ and ②…③ hydrogen bonds are the interfacial hydrogen bonds. The hydrogen bonds formed between the nano-silica and meta-aramid fibre are shown in Figure 4.
The number of hydrogen bonds in the model can be expressed as shown in the relational expression (1), where Ntotal means the total number of hydrogen bonds in the model, Npolymer means the number of meta-aramid fibre intramolecular and intermolecular hydrogen bonds, Nnanosilica means the total number of hydrogen bonds in nano-SiO2, and Ninterface means the number of the hydrogen bonds between the nanoparticle and meta-aramid fibre molecule. The statistics of different types of hydrogen bonds in all the models after the dynamics simulation are shown in Figure 5.
N total = N polymer + N nanosilica + N interface
It can be seen from Figure 5a that Ntotal of all the models doped with nanoparticles is greater than that of Model A, and this can indicate that doping with nanoparticles can beneficially increase the number of hydrogen bonds. With the increase in the radius of the nanoparticle, except for the total amount (Npolymer) of intramolecular or intermolecular hydrogen bonds of meta-aramid fibres decreasing, the total number (Ntotal) of the system hydrogen bonds, total number (Nnanosilica) of the hydrogen bonds on the nano-silica particles and total number (Ninterface) of the nanoparticle and meta-aramid insulating fibre interfacial hydrogen bonds increase, and the rate of increase of these values follows the trend of the radius. The reasons for this change are as follows. On the one hand, in the model built in this paper, the number of surface hydroxide radicals increases sharply with the increase of nanoparticle radius, and this leads to an increase in the formation probability of hydrogen bonds on the nanoparticle and at the interface. On the other hand, in the dynamics process, because the distance and bond angle between atoms change due to the torsion and deformation of the molecular chains of the meta-aramid fibres, the fibre is readily found close to the nanoparticle, and the formation probability of hydrogen bonds is further increased, therefore, the Nnanosilica and Ninterface values are increased. However, a decreasing trend is presented for the Npolymer value; this is mainly due to the fact that the greater the nanoparticle radius, the greater the volume of the mixed model, and the fewer contact opportunities for the atoms between the meta-aramid fibre chains, leading to the decrease of the formation probability of intramolecular or intermolecular hydrogen bonds in meta-aramid fibres. From the trend in the figure, when the nanoparticle radius is greater than 9 Å (Model E), the contribution of the Nnanosilica value to Ntotal value is the greatest; this is mainly due to the proportion of the nanoparticles gradually increasing, and the number of surface atoms and number of hydrogen bonds formed between the hydrogen and oxygen atoms on the nanoparticles increase with the increase of radius.
The change in the number of the different hydrogen bonds formed with the nanoparticle radius is shown in Figure 5b. In this figure, with the increase in radius, the increasing trend of nano-SiO2 intramolecular ①…② hydrogen bonds is evident, and when the radius is 7 Å, the nano-SiO2 bonds become the main contributor to the total number of hydrogen bonds. In the meta-aramid insulating paper fibre, a clear downtrend is presented for the number of ③…⑤ hydrogen bonds, while the number of ③…④ hydrogen bonds basically remains unchanged, and this indicates that a strong effect is generated on the ③…⑤ hydrogen bonds in the meta-aramid fibre from the dynamics motion limit. In the hydrogen bonds formed at the interface, the number of ①…⑤ hydrogen bonds keeps rising below the radius is 7 Å, and the number of hydroxide radicals on the nanoparticles increases sharply after the radius is 7 Å. The number of ①…⑤ hydrogen bonds shows a downward trend after the radius is 7 Å; the main reason for this is that the number of oxygen atoms contained in the meta-aramid fibres of the model remains the same and the motion of the meta-aramid fibre molecules are limited, which cause a decrease in the number of ①…⑤ hydrogen bonds instead of an increase when the radius is larger than 7 Å. Furthermore, in the hydrogen bonds formed at the interface, with the increase of radius, ②…③ hydrogen bonds and ①…④ hydrogen bonds basically show an upward trend, and ②…③ hydrogen bonds become the main contributor of the interfacial hydrogen bonds when the radius is greater than 7 Å. The number of ②…③ hydrogen bonds is always greater than the number of ①…④ hydrogen bonds, which indicates the bonding probability of the hydrogen atom bonded with the nitrogen atom in the meta-aramid fibre and the oxygen atom on the nano-silica molecule is greater than the bonding probability of the hydroxide radical on the nano-silica molecule and the N atom in the meta-aramid fibre.

3.2. Analysis of Formation Probability of Interatomic Hydrogen Bonds

3.2.1. Pair Correlation Function

The PCF can be represented with the symbol g(r), and the distribution of the local space particles becomes an effective method for studying the interaction between the material structure and the particle [54]. This means that the conditional probability density that a particle β is found in the range of the radius from r to r + dr, where the particle α is centred.
In this paper, the interaction condition between atoms and molecules is further obtained through the pair correlation function between these pairs of atoms. Here, Model C is taken as an example. The pair correlation functions for all the atoms are shown in Figure 6.
In Figure 6, the horizontal coordinate represents the distance between atoms, and the vertical ordinate represents the probability density [55]. The action range of hydrogen bonds is generally within 2.6–3.1 Å [54].
From Figure 6, when r = 0.95 Å, the first peak appears, similarly, in Figure 6d, when r = 1.05 Å and g(r) = 127.23, for the distance between atoms ③ and ④, a peak value appears. Because hydroxide radical exists in the nano-silica molecule and –NH bond exists in the meta-aramid fibre, the bond length of the oxygen atom and hydrogen atom in the hydroxide radical and the bond length of the nitrogen atom and hydrogen atom in the –NH bond are roughly equal to the distance (r) at the peak value. After a comparison of all the figures, it can be shown that the g(r) values are different with increasing radius, and this indicates that the occurrence probability of these atom pairs at the distance is different.
A further study is carried out through the hydrogen bond interactions of the model in the curve shown in Figure 6, where the statistics for certain hydrogen bond interactions is obtained by removing the curve with a radius of 2.6–3.1 Å. The calculation results are given by Formula (2), and if the formula is replaced with the indefinite integral of gij(r), the integral will be the same as in KB (Kirkwood-Buff) theory [56]. The larger S is, the greater the hydrogen bond formation probability is, and P indicates the hydrogen bond formation probability. Because the minimum value of g(r) is 0, the minimum value of S is ‒0.5.
S = 2.6 3.1 [ g ( r ) 1 ] d r
The calculation results are as shown in Figure 7.
It can be seen in Figure 7 that decrease of S is great with increasing radius of the nano-silica, and S and S maintains a downtrend that is relatively gentle. For S in Model D, there is a slight increase, but the overall trend is downward; this indicates that the hydrogen bond formation probability of this atom pair decreases with the increase in the radius of the nanoparticle. S exhibits a downtrend at first, with the minimum value occurring when the radius is equal to 7 Å, then an uptrend occurs, and the change of S is opposite to that of S. According to the changing relationships of the number of hydrogen bonds in the models, it can be seen that the number of hydrogen bonds is not positively related to the hydrogen bond formation probability of atom pairs, and the main reason is that the factors affecting the number of hydrogen bonds include the formation probability of a hydrogen bond and the number of atoms forming these hydrogen bonds. Although the formation probability of a hydrogen bond between atoms ① and ② decreases with the increase in radius, the number of hydroxide radicals and oxygen atoms on the nano-silica molecule surface increase sharply with the radius, and so this finally leads to an increase in the number of O–H…O (nanosilica) hydrogen bonds, with other atom pairs being similar.
P can be sorted in accordance with the value of S, taking Model C as an example, at 90 °C, the formation probability sequence of hydrogen bonds between the atom pairs in Model C is P > P > P > P > P > P.

3.2.2. Radial Distribution Function

There are three concepts of the radial distribution function, namely radial distribution functions of the wave function, the electron cloud and probability [57], and there are two definitions shown in Formulas (3) and (4) [58]. Formula (3) is adopted in this paper to express the radial distribution function [59], and the meaning of the formula is a statistical average of the atomic radial distribution in space.
G ( r ) = r 2 R 2 ( r )
G ( r ) = 4 π r 2 R 2 ( r )
The radial distribution function shows that if the particle α is deemed as the centre, the relationship of the probability of seeking a β particle in the range of the radius from r to r + dr with PCF is:
g ( r ) = G ( r ) 4 π r 2 ρ
where
ρ = N β V
ρ means the average density of particles, Nβ means the total number of β particles, and V means the space volume. The PCF is a “normalization” process result of the RDF.
In the RDF curve, the horizontal coordinate at the peak indicates the distance most likely to occur among atoms, the half width of the peak mainly shows the distribution of atomic spacing, and the peak area indicates the number of atoms distributed in the spacing, that is, the coordination number [60]. The distance range of hydrogen bonds is generally within 2.6–3.1 Å [54,55]; therefore, only total RDFs (atoms ①, ②, ③, ④ and ⑤) of all the models when r is within 2.6–3.1 Å are given, as shown in Figure 8.
A peak value of the model appears when r is within 2.75–2.80 Å, indicating that a hydrogen bond is most likely to appear at this position, and the occurrence probability of the hydrogen bonds of all the models in this position is different, with the sequence of Model E > Model B > Model A > Model C > Model D. When r is within 3.0–3.1 Å, the formation possibility of hydrogen bonds is small.
By analysing the formation probability of hydrogen bonds in atom pairs in Model C through PCF, the hydrogen bond formation probability sequence is obtained: P > P > P > P > P > P. The results are consistent with the results drawn by using the number of hydrogen bonds. From the radial distribution function, the hydrogen bond formation probability in the case of the radiuses of 2.75–2.80 Å in all the models is the greatest, the possibility to form hydrogen bonds within 3.0–3.1 Å is the smallest, and to increase the number of hydrogen bonds, the number of atoms within a radius in the range of 2.75–2.8 Å should be increased as far as possible. In the interfacial hydrogen bond, the formation probability of ①…⑤ hydrogen bonds is greatest, that is, an effective way to increase the number of interfacial hydrogen bonds in nano-SiO2 and meta-aramid fibres is to increase the number of hydrogen atoms on the nano-silica surface and oxygen atoms in the meta-aramid fibre.

3.3. Thermal Stability of Meta-Aramid Fibres

From the MSD, the chain motion of the polymer can be obtained, and the greater the slope in the time curve, the fiercer the chain movement of the polymer and the weaker the thermal stability of the polymer. The relationship between the MSD and time [61] is represented by Formula (7); in the formula, r i ( t ) represents the coordinates of the i-th atom at the moment of t in the system, and r i ( 0 ) represents the initial coordinates of the atom.
M S D = | r i ( t ) r i ( 0 ) | 2
The MSD of meta-aramid fibres in different models was shown in Figure 9. It indicates that the addition of nanoparticles reduces the strength of the meta-aramid fibres chain motion. For doping with nano-silica, a hydrogen bond is formed between the nanoparticle and the meta-aramid fibre, with binding for both caused by the interfacial hydrogen bond. The free movement of both is limited, thus, the chain movement of meta-aramid fibres is weakened. Therefore, compared with the unmodified model, the chain motion intensity of the modified model is clearly weakened.
Moreover, the inhibition degree of nano-silica with different radiuses to the meta-aramid chain motion varies. In a certain case, this is mainly due to the contact area of the nanoparticle and meta-aramid increasing with the radius of the nano-silica particle; therefore, the number of the interfacial hydrogen bonds increases, and the binding effect on the meta-aramid fibre strengthens. Therefore, the larger the radius is, the weaker the chain movement of meta-aramid fibres is.
As the dynamics of hydrogen bonds are very fast and their length scales are very short, a more detailed analysis is needed to decide either the hydrogen bonds are responsible for the differences in MSDs as shown or there might be processed occurring at larger time and length scales (such as Coulombic interactions) which are mainly responsible for this. It might be that the effect of hydrogen bonds at the interface does not extend to such large time and length scales to affect the whole chain translation. The same phenomenon is addressed in the literature in [62], in charged systems capable of hydrogen bond formation. Therefore, the binding of interfacial hydrogen bonds to the free movement of meta-aramid fibres is limited and some more detailed analysis will be presented in a future paper.

3.4. Effects on Interaction Energy

The binding energy of the molecular polymer and nanoparticle surface greatly affects the properties of the high polymer and nanoparticle composites, and the impact degree often depends on the structural characteristics and microchemical properties of the high polymer and nanoparticle interface [63]. The intensity of the intermolecular interaction can be represented with intermolecular interaction energy, Eint. If the value is less than 0, the molecular structure is steady, and the smaller the value is, the steadier the structural stability is. The interaction energy between the nanoparticle and meta-aramid is calculated with Formula (8).
E int = E total ( E SiO 2 + E ploymer )
where Etotal indicates the total potential energy of the meta-aramid and silica composite system, ESiO2 indicates the total potential energy of the nano-silica particle after the removal of the meta-aramid fibre, and Epolymer indicates the total potential energy of the meta-aramid fibre molecule chain after the removal of silica. Etotal consists of the EValence energy value and ENon-bond energy value, and the ENon-bond energy value consists of EHydrogen bond energy value, Evan der Waals energy value, EElectrostatic energy value, E3-Body energy value and ERestraint energy value. The energies of different models are shown in Table 1; therein, the EH-interface indicates the bond energy of hydrogen bonds formed at the nanoparticle and meta-aramid interface.
As shown in Table 1, with the increase in the nanoparticle radius, Eint increases in the negative direction, the amplitude increases, and the trend of EH-interface rises. This shows that EH-interface provides a positive contribution to Eint, and it indicates that the existence of interfacial hydrogen bonds is beneficial to the bonding of the two substances; that is, the interfacial hydrogen bond promotes the nano-silica particle to disperse in the meta-aramid fibres, to form a new type of insulating material. This shows that the existence of interfacial hydrogen bonds is one of the important reasons for the formation of the stable interface structure between nanoparticles and meta-aramid fibres. In spite of this, the hydrogen bond energy occupies a small proportion of the interaction energy; therefore, the effect of hydrogen bonds to the combination of both substances is still limited, and the effect of the degree of interfacial hydrogen bonds in different models on the interaction between the two substances is different. To analyse such a difference, the interfacial hydrogen bond energy density is utilized.
The interfacial hydrogen bond energy density here is defined as the ratio of the interfacial hydrogen bond energy and the surface area of the nanoparticle, which is expressed as ρH-interface. As given by Formula (9), where EH-interface indicates the interfacial hydrogen bond energy, and Snanoparticle indicates the surface area of the nanoparticle. Because the hydrogen bond energy is a negative value, the smaller the ρH-interface is, the stronger the effect of hydrogen bonds on the interaction between the two substances is. In the calculation, the circular constant is taken as 3.14, and the final results retain four decimal places as per the rounding-off method, as shown in Figure 10a.
ρ H interface = E H interface S nanoparticle
With the increase in radius, the ρH-interface gradually increases but the increase in amplitude is gradually reduced, which shows that with the increase of the nanoparticle radius, the effect of hydrogen bonds on the interaction between the two substances is weakened. Therefore, to improve the hydrogen bond interaction, a surface modification of the nanoparticle can be performed [64,65].
A similar function can be adopted to define the interfacial interaction energy density, ρint, as shown in Formula (10). For the different models, the change in ρH-interface and ρint are shown in Figure 10b. The two functions have similar laws of change, and the results show that with the increase in the nanoparticle radius, the interaction per unit surface area weakens. Therefore, to improve ρint and ρH-interface, when nano-silica is selected for the modification of meta-aramid fibres, nano-silica materials with small radiuses should be selected to the greatest extent.
ρ int = E int S nanoparticle
In summary, the existence of interfacial hydrogen bonds effectively restricts the free movement of meta-aramid fibres, reduces the chain movement of meta-aramid fibres, and increases the thermal stability of meta-aramid fibres. The interfacial hydrogen bond plays a critical role in the promotion of the bonding of the nano-silica particle and meta-aramid fibre. The interfacial hydrogen bond energy is an important constituent of the interaction energy, that is, the existence of interfacial hydrogen bonds promotes the interaction between the two substances, forming a relatively stable interface structure. Nanomaterials with small radiuses improve the interfacial hydrogen bond energy density and interface interaction energy density, enhancing the stability of the mixed model interface.

4. Conclusions

In this paper, the formation of hydrogen bonds between the nano-silica with different radiuses and the meta-aramid fibre interface is studied through molecular simulation. The hydrogen bond formation probability between atoms is analysed through the number of hydrogen bonds, pair correlation function and radial distribution function. The effect of the hydrogen bond on chain movement of meta-aramid fibres is analysed through the mean square displacement. The effect of the hydrogen bond on the interaction at the nano-silica and meta-aramid fibre interface is analysed by means of the interaction energy and hydrogen bond energy. The conclusions are as follows:
  • The total number of the hydrogen bonds of the mixed model increases with doping by nano-silica. With the increase of the nanoparticle radius, the number of hydrogen bonds on the nanoparticle surface and the mixed model interface of nano-SiO2/meta-aramid fibre increases; thus, the increasing trend of the –OH…O hydrogen bonds on the surface of nano-SiO2 is most evident.
  • For all of the hydrogen bond types of the SiO2/meta-aramid fibre mixed models, the hydrogen bond formation probability between the oxygen atom and hydrogen atom on the nanoparticle surface is the greatest. One of the effective ways to increase the hydrogen bonds in a SiO2/meta-aramid fibre is to increase the number of hydrogen atoms on the nano-silica surface and oxygen atoms in the meta-aramid fibre. Furthermore, to increase the formation probability of hydrogen bonds in all the models, the atomic distance should be enlarged to 2.7–2.8 Å as far as possible.
  • The existence of interfacial hydrogen bonds can restrict the free movement of meta-aramid fibres, lower the chain movement of meta-aramid fibres, and improve the thermal stability of fibres. Additionally, the existence of interfacial hydrogen bonds is one of the important reasons for the formation of the stable interface structure of nanoparticles and meta-aramid fibres. However, at larger time and length scales, some other influence factors (such as Coulombic interactions) which can act on the MSDs should be taken into consideration.
In the modification process of the meta-aramid fibre, although the increase of the nanoparticle radius results in an increase in the number of interfacial hydrogen bonds of the mixed model, generally, a nanoparticle with a small radius can preferably promote the interfacial hydrogen bond energy density and interface interaction energy density of the mixed model and then promote stability of the mixed model interface, facilitating a better nanoparticle doping effect.

Acknowledgments

The authors wish to thank the National Key R&D Program of China (Grant Nos. 2017YFB0902700, 2017YBF0902702) for the financial support.

Author Contributions

Chao Tang and Xu Li conceived and designed the simulations; Zhiwei Li and Jian Hao performed the simulations; all authors analyzed the data; Chao Tang and Xu Li wrote the paper, while other authors offered their modification suggestions for the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, H.; Chen, M.; Rice, P.M.; Wang, S.X.; White, R.L.; Sun, S.H. Dumbbell-like bifunctional Au–Fe3O4 nanoparticles. Nano Lett. 2005, 5, 379–382. [Google Scholar] [CrossRef] [PubMed]
  2. Salgueiriño-Maceira, V.; Correa-Duarte, M.A.; Spasova, M.; Liz-Marzán, L.M.; Farle, M. Composite silica spheres with magnetic and luminescent functionalities. Adv. Funct. Mater. 2006, 16, 509–514. [Google Scholar] [CrossRef]
  3. Zhang, D.H.; Zhang, H.; Zhang, Z.; Chen, Y.F. Industry applications of nanotechnology in high performance insulation composites. Sci. Sin. Chim. 2013, 43, 725–743. [Google Scholar] [CrossRef]
  4. Bai, G.; Liao, R.J.; Liu, N.; Liu, H.B.; Yang, L.J.; Shakeel, A. Influence of nano-AlN modification on the dielectric properties of meta-aramid paper. High Volt. Eng. 2015, 41, 461–467. [Google Scholar] [CrossRef]
  5. Chen, L. The development of M-aramid paper and its application. Chain Pulp. Ind. 2016, 37, 28–29. [Google Scholar] [CrossRef]
  6. Qin, M.; Kong, H.; Zhang, K.; Teng, C.; Yu, M.; Liao, Y. Simple synthesis of hydroxyl and ethylene functionalized aromatic polyamides as sizing agents to improve adhesion properties of aramid fiber/vinyl epoxy composites. Polymers 2017, 9, 143. [Google Scholar] [CrossRef]
  7. Schuh, C.; Schuh, K.; Lechmann, M.C.; Garnier, L.; Kraft, A. Shape-memory properties of segmented polymers containing aramid hard segments and polycaprolactone soft segments. Polymers 2010, 2, 71–85. [Google Scholar] [CrossRef]
  8. Seyhan, E.C.; Goksu, C.; Uzunhasanoglu, A.; Ilki, A. Seismic behavior of substandard RC columns retrofitted with embedded aramid fiber reinforced polymer (AFRP) reinforcement. Polymers 2015, 7, 2535–2557. [Google Scholar] [CrossRef]
  9. Lu, Z.Q.; Wei, N.; Ding, W.; Zhang, M.J. Improving the properties of polyimide fiber-based materials by nanosilica. J. Shaanxi Univ. Sci. Technol. 2017, 35, 6–10. [Google Scholar] [CrossRef]
  10. Liao, R.J.; Liu, H.B.; Bai, G.; Hao, J.; Li, X.; Lin, Y.D. Space charge characteristics and dielectric properties of nano-SiO2/aramid paper composite. Trans. CES 2016, 31, 40–48. [Google Scholar]
  11. Li, H.Y.; Wu, Q.L.; Zhou, D.G. Preparation and property analysis of cellulose nano-fibril and nano-silicon dioxide composites. Trans. CSAE 2015, 31, 299–303. [Google Scholar] [CrossRef]
  12. Ardalan, R.B.; Jamshidi, N.; Arabameri, H.; Joshaghani, A.; Mehrinejad, M.; Sharafi, P. Enhancing the permeability and abrasion resistance of concrete using colloidal nano-SiO2 oxide and spraying nanosilicon practices. Constr. Build. Mater. 2017, 146, 128–135. [Google Scholar] [CrossRef]
  13. Zhang, S.; Tang, C.; Chen, G.; Zhou, Q.; Lv, C.; Li, X. The influence and mechanism of nano Al2O3 to the thermal stability of cellulose insulation paper. Sci. Sin. Technol. 2015, 45, 1167–1179. [Google Scholar] [CrossRef]
  14. Yang, C.; Li, J.S.; Li, Y.R. Density-functional study of the adsorption of ZnO on α-Al2O3(0001) surface. Chem. J. Chin. Univ. 2004, 25, 1897–1900. [Google Scholar]
  15. Lewis, T.J. Interface are the dominaant feature of dielectrics at the nanometric level. IEEE Trans. Dielectr. Electr. Insul. 2004, 11, 739–753. [Google Scholar] [CrossRef]
  16. Cheong, W.C.D.; Zhang, L.C. Molecular dynamics simulation of phase transformations in silicon monocrystals due to nano-indentation. Nanotechnology 2000, 11, 173–180. [Google Scholar] [CrossRef]
  17. Nel, A.E.; Mädler, L.; Velegol, D.; Xia, T.; Hoek, E.M.V.; Somasundaran, P.; Klaessig, F.; Castranova, V.; Thompson, M. Understanding biophysicochemical interactions at the nano-bio interface. Nat. Mater. 2009, 8, 543–557. [Google Scholar] [CrossRef] [PubMed]
  18. Chang, X.L.; Zu, Y.; Zhao, Y.L. Size and structure effects in the nanotoxic response of nanomaterials. Chin. Sci. Bull. 2011, 56, 108–118. [Google Scholar] [CrossRef]
  19. Israelachvili, J.; Wennerström, H. Role of hydration and water structure in biological and colloidal interactions. Nature 1996, 379, 219–225. [Google Scholar] [CrossRef] [PubMed]
  20. Shi, X.; Xiao, H.; Lackner, K.S.; Chen, X. Capture CO2 from ambient air using nanoconfined ion hydration. Angew. Chem. 2016, 128, 4094–4097. [Google Scholar] [CrossRef]
  21. Shi, X.; Xiao, H.; Chen, X.; Lackner, K.S. The Effect of Moisture on the Hydrolysis of Basic Salts. Chem. Eur. J. 2016, 22, 18326–18330. [Google Scholar] [CrossRef] [PubMed]
  22. Li, Q.; Tang, Q.; Peng, T.; Zhang, X.; Liu, C.; Shi, X. Molecular characteristics of H2O in hydrate/ice/liquid water mixture. Int. J. Mod. Phys. B 2015, 29, 1550185. [Google Scholar] [CrossRef]
  23. Li, Q.; Xiao, Y.; Shi, X.; Song, S. Rapid evaporation of water on graphene/graphene-oxide: A molecular dynamics study. Nanomaterials 2017, 7, 265. [Google Scholar] [CrossRef] [PubMed]
  24. Li, X.; Tang, C.; Wang, Q.; Li, X.P.; Hao, J. Molecular simulation research on the micro effect mechanism of interfacial properties of nano SiO2/meta-aramid fiber. Int. J. Heat Technol. 2017, 35, 123–129. [Google Scholar] [CrossRef]
  25. Heinera, A.P.; Kuuttib, L.; Telemana, O. Comparison of the interface between water and four surfaces of native crystalline cellulose by molecular dynamics simulations. Carbohydr. Res. 1998, 306, 205–220. [Google Scholar] [CrossRef]
  26. Yin, K.L.; Zou, D.H.; Yang, B.; Zhang, X.H.; Xia, Q.; Xu, D.J. Investigation of H-bonding for the related force fields in materials studio software. Comput. Appl. Chem. 2006, 23, 1335–1340. [Google Scholar] [CrossRef]
  27. Li, F.; Dong, J.Q.; Shen, Q. The function of hydrogen bond in polymer blends I. Characters and affecting factors of hydrogen bond. Chin. Polym. Bull. 2009, 7, 45–52. [Google Scholar] [CrossRef]
  28. Li, F.; Dong, J.Q.; Shen, Q. The function of hydrogen bond in polymer blends II. Characters and affecting factors of hydrogen bond. Chin. Polym. Bull. 2009, 8, 54–65. [Google Scholar] [CrossRef]
  29. Xiao, H.; Shi, X.; Zhang, Y.; Liao, X.; Hao, F.; Lackner, K.S.; Chen, X. The catalytic effect of water in basic hydrolysis of CO32− in hydrated clusters and its implication in the humidity driven CO2 air capture. Phys. Chem. Chem. Phys. 2017. [Google Scholar] [CrossRef] [PubMed]
  30. Eslami, H.; Müller-Plathe, F. Molecular dynamics simulation of water influence on local structure of nanoconfined polyamide-6,6. J. Phys. Chem. B 2011, 115, 9720–9731. [Google Scholar] [CrossRef] [PubMed]
  31. Eslami, H.; Rahimi, M.; Müller-Plathe, F. Molecular dynamics simulation of a silica nanoparticle in oligomeric poly(methyl methacrylate): A model system for studying the interphase thickness in a polymer–nanocomposite via different properties. Macromolecules 2013, 46, 8680–8692. [Google Scholar] [CrossRef]
  32. Eslami, H.; Müller-Plathe, F. Structure and mobility of nanoconfined polyamide-6,6 oligomers: Application of a molecular dynamics technique with constant temperature, Surface Area, and parallel pressure. J. Phys. Chem. B 2009, 113, 5568–5581. [Google Scholar] [CrossRef] [PubMed]
  33. Eslami, H.; Müller-Plathe, F. How thick is the interphase in an ultrathin polymer film? Coarse-grained molecular dynamics simulations of polyamide-6,6 on graphene. J. Phys. Chem. C 2013, 117, 5249–5257. [Google Scholar] [CrossRef]
  34. Wang, J.Z.; Zhao, Y.; Zhang, C.B. Study of mono-dispersed spherical silica in sol-gel process with orthogonal design. China Powder Sci. Technol. 2003, 9, 4–8. [Google Scholar] [CrossRef]
  35. Bao, M.R.; Zhu, G.R.; Wang, M.; Gao, C.J. Progress on preparation of mono-dispersed spherical nano-silica. Mater. Rev. 2011, 25, 135–139. [Google Scholar]
  36. Li, Z.; Okamoto, K.; Ohki, Y.; Tanaka, T. The role of nano and micro particles on partial discharge and breakdown strength in epoxy composites. IEEE Trans. Dielectr. Electr. Insul. 2011, 18, 675–681. [Google Scholar] [CrossRef]
  37. Roman, K.; Thomas, A.; Peter, H.F.M.; Johan, J.S. Anomalous behaviour of the dielectric spectroscopy response of nanocomposites. IEEE Trans. Dielectr. Electr. Insul. 2012, 19, 107–117. [Google Scholar] [CrossRef]
  38. Yang, K.D.; Lee, C.W.; Jang, J.H.; Ha, T.R.; Nam, K.T. Rise of nano effects in electrode during electrocatalytic CO2 conversion. Nanotechnology 2017, 28, 352001. [Google Scholar] [CrossRef] [PubMed]
  39. Zha, J.W.; Zhu, Y.H.; Li, W.K.; Bai, J.B.; Dang, Z.M. Low dielectric permittivity and high thermal conductivity silicone rubber composites with micro-nano-sized particles. Appl. Phy. Lett. 2012, 101. [Google Scholar] [CrossRef]
  40. Zha, J.W.; Li, W.K.; Zhu, Y.H.; Dang, Z.M.; Chen, G. Effect of Micro-Si3N4-nano-Al2O3 Cofilled Particles on Thermal Conductivity, Dielectric and Mechanical Properties of Sillicone Rubber Composites. IEEE Trans. Dielectr. Electr. Insul. 2014, 21, 1989–1996. [Google Scholar] [CrossRef]
  41. Xu, N.Q.; Gu, J.X.; Luo, K.; Zheng, S.B. Preparation of nanon silicon dioxide of large particle size. Chem. Ind. Eng. Prog. (China) 2004, 23, 747–750. [Google Scholar] [CrossRef]
  42. Qamar, R.; Mohtashim, L.; Elke, D.; Heidemarie, P.; Ludwig, J.; Dieter, G.W.; Dietmar, S. Evidence that ultrafine titanium dioxide induces micronuclei and apoptosis in syrian hamster embryo fibroblasts. Environ. Health Perspect. 2012, 110, 797–800. [Google Scholar]
  43. Park, M.V.D.Z.; Annema, W.; Salvati, A.; Lesniak, A.; Elsaesser, A.; Barnes, C.; McKerr, G.; Howard, C.V.; Lynch, I.; Dawson, K.A.; et al. In vitro developmental toxicity test detects inhibition of stem cell differentiation by silica nanoparticles. Toxicol. Appl. Pharmacol. 2009, 240, 108–116. [Google Scholar] [CrossRef] [PubMed]
  44. Li, A.S.; Huang, Y.P.; Liu, J.W.; Lan, M.B.; Gao, F.; Xu, H.L.; Chen, M.C.; Xu, Y.C. Oxidative lesions induced by medical SiO2 nanoparticles on normal human lung cells. Chin. J. Clin. Pharmacol. Ther. 2009, 14, 1115–1120. [Google Scholar]
  45. Yu, K.O.; Grabinski, C.M.; Schrand, A.M.; Murdock, R.C.; Wang, W.; Gu, B.H.; Schlager, J.J.; Hussain, S.M. Toxicity of amorphous silica nanoparticles in mouse keratinocytes. J. Nanopart. Res. 2009, 11, 15–24. [Google Scholar] [CrossRef]
  46. Napierska, D.; Thomassen, L.C.J.; Rabolli, V.; Lison, D.; Gonzalez, L.; Kirsch-Volders, M.; Martens, J.A.; Hoet, P.H. Size-dependent cytotoxicity of monodisperse silica nanoparticles in human endothelial cells. Small 2009, 5, 846–853. [Google Scholar] [CrossRef] [PubMed]
  47. Chelli, R.; Procacci, P.; Cardini, G.; Califano, S. Glycerol condensed phases Part II. A molecular dynamics study of the conformational structure and hydrogen bonding. Phys. Chem. Chem. Phys. 1999, 1, 879–885. [Google Scholar] [CrossRef]
  48. Chelli, R.; Procacci, P.; Cardini, G.; Valle, R.G.D.; Califano, S. Glycerol condensed phases Part I. A molecular dynamics study. Phys. Chem. Chem. Phys. 1999, 1, 871–877. [Google Scholar] [CrossRef]
  49. Chen, C.; Li, W.Z. Molecular dynamics simulation of hydrogen bonding characteristics in aqueous glycerol solutions. Acta Phys. Chim. Sin. 2009, 25, 507–512. [Google Scholar]
  50. Callam, C.S.; Singer, S.J.; Lowary, T.L.; Hadad, C.M. Computational analysis of the potential energy surfaces of glycerol in the gas and aqueous phases: Effects of level of theory, basis set, and solvation on strongly intramolecularly hydrogen-bonded systems. J. Am. Chem. Soc. 2001, 123, 11743–11754. [Google Scholar] [CrossRef] [PubMed]
  51. Guàrdia, E.; Martí, J.; Padró, J.A.; Saiz, L.; Komolkin, A.V. Dynamics in hydrogen bonded liquids: Water and alcohols. J. Mol. Liq. 2002, 96–97, 3–17. [Google Scholar] [CrossRef]
  52. Chen, H.T.; Wang, J.F.; Wu, X.L.; Chen, X.H. Technological conditions and properties of nano-SiO2 modified epoxy resin. Insul. Mater. 2008, 41, 40–43. [Google Scholar] [CrossRef]
  53. Li, C.Y.; Qiu, D.J.; Xie, G.X.; Xiao, X.D.; Ni, X.X. Study on nano-silica toughening modification for epoxy resin. Mater. Prot. 2007, 41, 21–23. [Google Scholar] [CrossRef]
  54. Liu, D.M. Molecular Dynamics Simulation Study on the Structure and Properties of HMX and PETN and Their Composite Materials. Master’s Thesis, Nanjing University of Science & Technology, Nanjing, China, 2014. [Google Scholar]
  55. Zhang, Y.; Gu, J.H.; Shang, J.X.; Ma, Y. Fundamentals of Computational Materials, 1st ed.; Beihang University Press: Beijing, China, 2007; pp. 121–122. ISBN 978-7-81077-788-9. [Google Scholar]
  56. Wernersson, E.; Stenqvist, B.; Lund, M. The mechanism of cellulose solubilization by urea studied by molecular simulaton. Cellulose 2015, 22, 991–1001. [Google Scholar] [CrossRef]
  57. Liu, X.; Luo, D.F. The concept and discussion of the radial distribution function. J. Liaocheng Teach. Univ. (Nat. Sci.) 2002, 15, 48–49. [Google Scholar]
  58. Li, B.R.; Li, D. Two definitions of radial distribution function of atomic orbits and their correlations. Chem. Bull. 1982, 4, 62–63. [Google Scholar] [CrossRef]
  59. Li, B.R. Discrimination of common misconceptions about atomic orbits and electronic cloud patterns. Chem. Bull. 1984, 3, 46–50. [Google Scholar] [CrossRef]
  60. Hu, J. Molecular Simulation Studies on Aging Mechanism in Crazing Zone of the Insulation Paper in Power Transformers. Master Thesis, Chongqing University, Chongqing, China, 2009. [Google Scholar]
  61. Hofmann, D.; Fritz, L.; Ulbrich, J.; Paul, D. Molecular simulation of small molecule diffusion and solution in dense amorphous polysiloxanes and polyimides. Comput. Theor. Polym. Sci. 2000, 10, 419–436. [Google Scholar] [CrossRef]
  62. Afandak, A.; Eslami, H. Ion-Pairing and Electrical Conductivity in the Ionic Liquid 1-n-Butyl-3-methylimidazolium Methylsulfate [Bmim][MeSO4]: Molecular Dynamics Simulation Study. J. Phys. Chem. B 2017, 121, 7699–7708. [Google Scholar] [CrossRef] [PubMed]
  63. Tong, H.Y.; Shi, F.G.; Zhao, B.; Wang, S.Q.; Brongo, M.; Vasudev, P.K. Interfacial microstructure and reaction at the spin-coated fluorinated polyimide/Al interface: Surface enhanced X-ray diffraction and TEM studies. Appl. Phys. A 1997, 65, 287–290. [Google Scholar] [CrossRef]
  64. Wang, B. Surface Modification of Inorganic Nanoparticles by Macromolecules. Master’s Thesis, Tsinghua University, Beijing, China, 2005. [Google Scholar]
  65. Jesionowski, T.; Krysztafkiewicz, A. Influence of silane coupling agents on surface properties of precipitated silicas. Appl. Surf. Sci. 2001, 172, 18–32. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram for the definition of a hydrogen bond.
Figure 1. Schematic diagram for the definition of a hydrogen bond.
Polymers 09 00504 g001
Figure 2. Meta-aramid fibre molecular model, polymerization degree (DP) = 10.
Figure 2. Meta-aramid fibre molecular model, polymerization degree (DP) = 10.
Polymers 09 00504 g002
Figure 3. Model building process and built models; Model A is the pure meta-aramid fibre model, and Model B, Model C, Model D and Model E are the mixed models of the SiO2/meta-aramid fibre.
Figure 3. Model building process and built models; Model A is the pure meta-aramid fibre model, and Model B, Model C, Model D and Model E are the mixed models of the SiO2/meta-aramid fibre.
Polymers 09 00504 g003
Figure 4. Hydrogen bonds formed between the nano-silica and meta-aramid fibre, ① indicates the hydrogen atom on the nano-silica surface, ② indicates the oxygen atom on the nano-silica surface, ③ indicates the hydrogen atom in the –NH radical in the meta-aramid fibre molecule, ④ indicates the nitrogen atom in the –NH radical in the meta-aramid fibre molecule, and ⑤ indicates the oxygen atom in the meta-aramid fibre molecule.
Figure 4. Hydrogen bonds formed between the nano-silica and meta-aramid fibre, ① indicates the hydrogen atom on the nano-silica surface, ② indicates the oxygen atom on the nano-silica surface, ③ indicates the hydrogen atom in the –NH radical in the meta-aramid fibre molecule, ④ indicates the nitrogen atom in the –NH radical in the meta-aramid fibre molecule, and ⑤ indicates the oxygen atom in the meta-aramid fibre molecule.
Polymers 09 00504 g004
Figure 5. Number of hydrogen bonds, (a) indicates the number of different types of hydrogen bonds in all the models; and (b) indicates the number of hydrogen bonds between different atoms at the interface.
Figure 5. Number of hydrogen bonds, (a) indicates the number of different types of hydrogen bonds in all the models; and (b) indicates the number of hydrogen bonds between different atoms at the interface.
Polymers 09 00504 g005
Figure 6. Pair correlation function for different types of hydrogen bond atoms of Model C at 90 °C: (a) ①…② hydrogen bonds; (b) ①…④ hydrogen bonds; (c) ①…⑤ and ②…③ hydrogen bonds; (d) ③…④ and ③…⑤ hydrogen bonds.
Figure 6. Pair correlation function for different types of hydrogen bond atoms of Model C at 90 °C: (a) ①…② hydrogen bonds; (b) ①…④ hydrogen bonds; (c) ①…⑤ and ②…③ hydrogen bonds; (d) ③…④ and ③…⑤ hydrogen bonds.
Polymers 09 00504 g006
Figure 7. Relationship between all the atom pairs (S) and the radius of nano-silica at 90 °C.
Figure 7. Relationship between all the atom pairs (S) and the radius of nano-silica at 90 °C.
Polymers 09 00504 g007
Figure 8. Total radial distribution functions (RDFs) of all models when r is within 2.6–3.1 Å.
Figure 8. Total radial distribution functions (RDFs) of all models when r is within 2.6–3.1 Å.
Polymers 09 00504 g008
Figure 9. Mean square displacement of all the models.
Figure 9. Mean square displacement of all the models.
Polymers 09 00504 g009
Figure 10. Change in bond energy density with the nanoparticle radius: (a) the change in hydrogen bond energy density with the radius of the nanoparticle and (b) the change in the interaction energy density with the nanoparticle radius.
Figure 10. Change in bond energy density with the nanoparticle radius: (a) the change in hydrogen bond energy density with the radius of the nanoparticle and (b) the change in the interaction energy density with the nanoparticle radius.
Polymers 09 00504 g010
Table 1. Energies of different models (kcal/mol) (1 cal = 4.1868 J).
Table 1. Energies of different models (kcal/mol) (1 cal = 4.1868 J).
ModelEintESiO2EPolymerEtotalEH-interface
B−432.3091−636.4831−1171.2478−2240.0400−53.4254
C−620.5648−856.8470−1476.5172−2953.9290−99.2254
D−1116.4885−5631.9281−1316.9705−8065.3871−147.7210
E−1742.1470−13,262.3598−1506.8031−16,511.3099−214.4546

Share and Cite

MDPI and ACS Style

Tang, C.; Li, X.; Li, Z.; Hao, J. Interfacial Hydrogen Bonds and Their Influence Mechanism on Increasing the Thermal Stability of Nano-SiO2-Modified Meta-Aramid Fibres. Polymers 2017, 9, 504. https://doi.org/10.3390/polym9100504

AMA Style

Tang C, Li X, Li Z, Hao J. Interfacial Hydrogen Bonds and Their Influence Mechanism on Increasing the Thermal Stability of Nano-SiO2-Modified Meta-Aramid Fibres. Polymers. 2017; 9(10):504. https://doi.org/10.3390/polym9100504

Chicago/Turabian Style

Tang, Chao, Xu Li, Zhiwei Li, and Jian Hao. 2017. "Interfacial Hydrogen Bonds and Their Influence Mechanism on Increasing the Thermal Stability of Nano-SiO2-Modified Meta-Aramid Fibres" Polymers 9, no. 10: 504. https://doi.org/10.3390/polym9100504

APA Style

Tang, C., Li, X., Li, Z., & Hao, J. (2017). Interfacial Hydrogen Bonds and Their Influence Mechanism on Increasing the Thermal Stability of Nano-SiO2-Modified Meta-Aramid Fibres. Polymers, 9(10), 504. https://doi.org/10.3390/polym9100504

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop