Next Article in Journal
Kinetic Study of the Diels–Alder Reaction between Maleimide and Furan-Containing Polystyrene Using Infrared Spectroscopy
Previous Article in Journal
Machine Learning Backpropagation Prediction and Analysis of the Thermal Degradation of Poly (Vinyl Alcohol)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Polymers
Volume 16
Issue 3
10.3390/polym16030438
polymers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MD-DFT Calculations on Dissociative Absorption Configurations of FOX-7 on (001)- and (101)-Oriented Crystalline Parylene Protective Membranes

by
Weihui Luo
1,
Liang Bian
1,2,3,4,*,
Faqin Dong
1,
Jianan Nie
1 and
Jingjie Yang
1,3
1
Key Laboratory of Solid Waste Treatment and Resource Recycle, Ministry of Education, South West University of Science and Technology, Mianyang 621010, China
2
Institute of Gem and Material Technology, Hebei GEO University, Shijiazhuang 050000, China
3
State Key Laboratory of Environment-friendly Energy Materials, South West University of Science and Technology, Mianyang 621010, China
4
Tianfu Institute of Research and Innovation, Southwest University of Science and Technology, Chengdu 610299, China
*
Author to whom correspondence should be addressed.
Polymers 2024, 16(3), 438; https://doi.org/10.3390/polym16030438
Submission received: 19 February 2023 / Revised: 29 September 2023 / Accepted: 7 October 2023 / Published: 5 February 2024
(This article belongs to the Section Polymer Physics and Theory)
Browse Figures
Review Reports Versions Notes

Abstract

:
Crystalline poly-para-xylylene (parylene) has the potential for use as a protective membrane to delay the nucleation of explosives by separating the explosives and their decomposition products to decrease the explosive sensitivity. Here, molecular dynamics (MD) and density functional theory (DFT) techniques were used to calculate the dissociative adsorption configurations of 1,1-diamino-2,2-dinitroethylene (FOX-7) on (001)- and (101)-oriented crystalline parylene membranes. Based on the results of the calculations, this work demonstrates that the -NO2–π electrostatic interactions are the dominant passivation mechanism of FOX-7 on these oriented surfaces. FOX-7 can dissociatively adsorb on oriented parylene membranes due to the interactions between the LUMO of the toluene (or methyl) groups on parylene and the HOMO of the -NO2 (or -NH2) groups on FOX-7. The formation of a new intermolecular H-bond with the ONO group leads to FOX-7 decomposition via intramolecular C-NO2 bond fission and nitro-to-nitrite rearrangement. The most likely adsorption configurations are described in terms of the decomposition products, surface active groups of parylene, binding behaviors, and N charge transfer. Importantly, the (001)-oriented parylene AF8 membrane is promising for use as a protective membrane to passivate the high-energy -NO2 bonds during the dissociative adsorption of FOX-7. This study offers a new perspective on the development of protective membranes for explosives.

1. Introduction

The novel high-energy material 1,1-diamino-2,2-dinitroethylene (FOX-7, DADNE) was first synthesized in 1998 and instantly received much attention due to its high thermal stability and low impact and friction sensitivities [1]. FOX-7 has a highly polarized carbon–carbon double bond, and the positive and negative charges on the carbons are stabilized by two amino and two nitro groups [2]. It is a nitroenamine with a “push-pull” alkene structure [3]. Under long-term storage conditions, these structures decompose due to the significant aging effects of temperature [4], pressure [5], and pH [6]. NVT (Canonical Ensemble) and AIMD (Ab Initio Molecular Dynamics) simulations revealed that the thermal decomposition mechanism of FOX-7 involves three different initial decomposition steps, namely, C-NO2 bond fission and intermolecular and intramolecular hydrogen transfer [7,8]. The C-NO2 bond fission mechanism was determined to be the most common route [9]. Because this fission mechanism causes performance issues, particularly for FOX-7 gas and its decomposition products, this explosive requires the use of a protective membrane overcoat that inhibits the slow decomposition [10]. For example, the oxygen atoms of the -NO2 groups in FOX-7 partially or completely dissociate from the molecule to oxidize the matrix surface [11]. The FOX-7 molecule decomposes into -ONO and -OC(NO2)=CN2H4, which have unsaturated bonds with high electrostatic potentials that can attract or repel the surface electrons of protective membranes [12]. Hence, a dissociative adsorption mechanism is possible, depending on the particular surface sites involved [13]. Therefore, the protective membranes required for this compound should have low sensitivity, good flexibility, and mechanical strength, and should be incompatible with organic molecules [14].
Recently, it was suggested that polymers can be coated on the surfaces of energetic materials (such as explosives, propellants, and pyrotechnics) to act as protective membranes [15]. Surface crystalline polymers are well known to decrease explosive sensitivity and infiltration [16]. Moreover, the explosive cannot crystallize on the surface of the protective membrane. These membranes can delay the nucleation of the explosive by inhibiting its volumetric expansion [17] and separating the explosive molecules from their decomposition products, hindering heat conduction [18]. Of the various possible polymer membranes, pinhole-free poly-para-xylylene (parylene) and substituted parylene membranes have high stabilities, high chemical and electrical resistances, and low gas permeability [19,20]. Substituted parylene membranes, such as poly-chloro-para-xylylene (parylene C), poly-tetrafluoro-para-xylylene (parylene F), α,α,α′,α′-poly-tetrafluoro-para-xylylene (parylene AF4), and fluorinated poly-para-xylylene (parylene AF8) membranes, were shown to exhibit improved gas impermeability [21]. The primary roles of these membranes are to reduce gas permeation and prevent interactions with organic compounds [22]. The stabilities of fluoroelastomer-coated cyclotrimethylenetrinitramine (RDX) and cyclotetramethylene tetranitramine (HMX) were previously verified [23]. The polar C-Cl and C-F bonds in parylene cause the formation of polymorphic structures [24], including those with (101) [25] and (001) [26] orientations. The exposed C-Cl and C-F bonds are different on the different surfaces [27], which improves the mitigation effects on the explosive thermolysis pathways [28]. Consequently, in this study, based on previous studies of polymer–explosive systems that examined their interfacial structures, phase transitions, heat conduction barriers, etc., four parylene protective membranes were investigated for use as a barrier for FOX-7 and its decomposition products.
Currently, few experimental data are available to explain the efficiency of energetic compounds coated with parylene membranes. Moreover, the mechanism describing the interactions between FOX-7 (or its decomposition products) and parylene surfaces is unknown. To elucidate the fundamental issues related to the nitro-to-nitrite rearrangement of FOX-7 deposited on (101)- and (001)-oriented parylene membranes, this work presents an atomic-level description of the intermolecular hydrogen bonds (H-bonds) formed between the decomposition products of FOX-7 and parylene surfaces. Specifically, MD and DFT calculations were used to study the corresponding local chemical adsorption processes.

2. Computational Details

The chain units of various parylene (parylene C, parylene F, parylene AF4, and parylene AF8) membranes were constructed using molecular mechanics (MM), and optimized by 200 ps NPT and NVT (NPT 100 ps before NVT 100 ps) molecular dynamics (MD) simulations at 298 K and 101 KPa using the COMPASS force field [21]. In this work, gradually reduced size (GRS) methods were employed to optimize the weak interaction potentials, reducing the error in the van der Waals and electrostatic potentials between chain units, and regular arrangements of (101)- and (001)-oriented parylene membranes were constructed. These membranes were optimized and relaxed via 110 ps NPT annealing and 500 ps NPT and NVT MD simulations, and the T and P control method was set to Nose Q ratio with 0.01 timestep and Berend decay constant with 0.1 timestep. Short-range van der Waals interactions were calculated using an atom-based method. The Ewald+Group summation method was routinely used to evaluate the long-range electrostatic interactions in reasonably small models [29].
The surface electrostatic energies of the crystalline parylene membranes were determined by density functional theory (DFT) calculations within the generalized gradient approximation (GGA) using the Perdew–Burke–Ernzerhof (PBE) functional. Initially, several tests were performed to verify the accuracy of the method for bulk FOX-7 and its decomposition products on (101)- and (001)-oriented parylene membranes. Specifically, both the physical atomic distributions and chemical electron transfer processes were investigated. The atomic distances and total energies of the models of FOX-7 and its decomposition products on the parylene surfaces were tested for convergence [30,31]. The atomic configurations were optimized by conjugate gradient minimization of the total energy via 500 ps NPT and NVT MD simulations [32]. To elucidate the interactions between the dissociatively adsorbed FOX-7 molecules and the polymer surfaces, the electron localization functions (ELFs) were optimized [11]. The corresponding kinetic energies were calculated by a spin-polarized GGA-PBE DFT method [33]. Other areas of the model and the relative boundary area were simulated by MM using periodic boundary conditions [34].

3. Results and Discussion

3.1. Surface Structures of the (001)- and (101)-Oriented Crystalline Parylene Protective Membranes

The main parylene polymers used as protective membranes include parylene C ((C8H7Cl)n), parylene F ((C8H7F)n), parylene AF4 ((C8H4F4)n), and parylene AF8 ((C8F8)n) [21,29], which are shown in Figure 1b. This figure shows that the errors in the bond lengths and angles are less than 0.5%, which is in contrast to the literature data [35]. The high-energy areas are the most likely sites for the condensation of adjacent chain units to construct low-energy parylene membranes (Figure 1c). Thus, the surface electrons in the -C-Cl (or -C-F), benzyl (π), or -C-F2 (or -C-H2) bonds of the parylene chain units tend to overlap. For example, the (101) and (001) orientations of parylene C were constructed by exploiting the coupling interactions between the π–π and -C-Cl-Cl-C- bonds, as shown in Figure 1a. The four optimized crystalline parylene configurations are given in Table 1. The calculated crystalline structures and densities are close to the experimental data with errors of less than 1% [21,35]. The errors in the calculated Tg values, which depend on the number and nature of the Cl and F atoms, are 4.4%, 6.5%, 8.7%, and 9.4%, which is in contrast to the literature data [21,35,36,37,38,39,40,41].
The parylene surfaces are represented by slab models of these membranes with a 1 nm vacuum. Specifically, (101)- and (001)-oriented parylene membranes with two layers containing 1920 and 2880 atoms were used to investigate the adsorption of FOX-7 and its decomposition products. Typically, the electrons in the benzene π-bonds are the closest to the surface, as shown in Figure 2. Although the chlorine (Cl) and fluorine (F) atoms pull electrons towards themselves [22], the ortho-hydrogen (H) atoms in benzene have the highest electrostatic densities, as shown by the highest occupied molecular orbital (HOMO) analysis. They are the source of the high electrostatic energies (25–70 eV) of parylene C (or parylene F), as shown in Figure 3. However, replacing the H atoms on the benzyl groups with fluorine (F) atoms increases the total surface electron densities [21]. This substitution results in the transfer of the active electrons to the benzyl groups and the formation of the lowest unoccupied molecular orbital (LUMO). The -CF2 bonds in parylene AF4 and parylene AF8 have the highest electrostatic energies (350–385 eV). These surface potential energies correspond to positive charges, which facilitate the adsorption of the negatively charged FOX-7 (surface potential energy: −255–2235 eV).

3.2. Adsorption Configurations of FOX-7 on the (001)- and (101)-Oriented Crystalline Parylene Protective Membranes

After determining the potential adsorption sites on the (101)- and (001)-oriented parylene membranes, the possible adsorption configurations of the FOX-7 molecule on these parylene membranes were investigated, as illustrated in Figure 4. These configurations were obtained by defining the distributions of the O and H atoms. The intermolecular H-bond lengths between the -NH2 and -NO2 groups in FOX-7 are approximately 0.2 nm [12,42]. The -NH2 and -NO2 bonds rearrange through interactions with adjacent active functional groups, a process that was tested and verified by experiments and simulations [3]. On the surfaces of the FOX-7 molecules, the electrons are transferred between the O-O and H-O bonds through 2p4 → 2p4 and 1s1 → 2p4 orbital processes [10,11]. The O-O bonds have high electron transfer ratios, that is, the O-O p → p orbitals have high electrostatic energies that facilitate electron transfer between the adjacent active groups [12]. Figure 5a indicates that the -NO2 groups on FOX-7 first bind to the π-bonds of benzene in (101)-oriented parylene C and F or to the benzyl groups in parylene AF4 and AF8. These electrostatic interactions cause the FOX-7 surface electrons to overlap with the parylene membranes, leading to the formation of new short-range intermolecular H-bonds. The adjacent -NH2 groups continue to interact with the π-bonds to provide the missing electrons to the -NO2-benzene unit until the O-H bonds break [13]. The corresponding energy differences (ΔE) between the adsorption models and single systems are the lowest (parylene C: −11 kcal∙mol−1; parylene F: −130 kcal∙mol−1; parylene AF4: −7.5 kcal∙mol−1; parylene AF8: −10 kcal∙mol−1), as shown in Figure 5b. Thus, the FOX-7 molecules lie flat on the (101)-oriented parylene membranes because the HOMO on the -NO2 group of FOX-7 interacts with the LUMO on the -C-Cl group on the parylene benzene to lower the energy of the system, as shown by the HOMO–LUMO results presented in Figure 6. The electrostatic interactions are considered to be the dominant adsorption mechanism for FOX-7 on the (101)-oriented parylene membranes. The -NO2 bonds (FOX-7), which are the high-energy origin of detonation, are passivated at the two basic adsorption sites, that is, at the toluene and methyl groups, due to the electron transfer from the FOX-7 HOMO to the parylene LUMO [2,3]. Simultaneously, the dissociative adsorption of FOX-7 is also possible and involves the breaking of the intramolecular H-bond in FOX-7 (-NO2-H2N-) and the formation of a new intermolecular H-bond, which has the smallest calculated bond length (0.2–0.24 nm) with the ONO group.
Another useful tool for investigating the coupling interactions between the FOX-7 active groups and the (101)-oriented parylene surfaces is the electron localization function (ELF). Isosurfaces of the ELF can be calculated and used to explain bond and charge data. The relative energy differences and bond energies are discussed below and are shown in Figure 7. The results of these calculations indicate that both the -NO2 and -NH2 groups of FOX-7 are bound to toluene groups based on the lowest energy principle. The surface -Cl (or -F) atoms of parylene preferentially capture the surface electrons of the benzene groups, which results in high-energy active electrons that can be transferred to the unsaturated -NH2 electron orbitals on FOX-7 [12]. Then, the paired positive holes bind to the negatively charged -NO2 groups. In contrast to the electron-withdrawing nature of the -benzene-Cl(F) bonds, the -CF2- surface potentials homogenize and push the surface electrons of benzene [23]. These low-energy sites are positively charged and thus preferentially bind to the negatively charged -NO2 groups. The surface potentials at the negatively charged methyl bonds are relatively uniform, and these bonds interact with the positively charged -NH2 bonds. A common tool for providing a semi-quantitative measure of the charge transfer is the Mulliken charge analysis in which the electronic charge is partitioned among the individual atoms, as shown in Figure 8. Although all the N atoms in the FOX-7 molecule have negative Mulliken charges on the parylene surfaces [43], the N atoms in both the -NH2 and -NO2 groups have highly negative charges (−2–0.2 e) on the (101)-oriented parylene membranes.
In contrast to the observed interactions between both the -NH2 and -NO2 groups and the (101)-oriented parylene membranes during adsorption, active methyl groups with positive charges participate in the adsorption process on the (001)-oriented parylene membranes. Based on the lowest energy principle, the exposed paired positive holes on the benzenes can capture the active electrons of the -NO2 bonds. The adsorption configurations of FOX-7 on the (001)-oriented parylene AF4 or AF8 membranes are similar to those on the (101)-oriented membranes. However, it should be noted that the -NO2 surface electrons are transferred to the exposed low-energy orbitals on the toluene groups on the (001)-oriented parylene C and F membranes, as shown by the lowest energy results presented in Figure 7. Due to the coupling interaction between the high- and low-energy areas of adjacent -benzene-Cl(F) bonds, the exposed methyl groups do not have enough active electrons to break the H-bonds in the FOX-7 molecules. This conclusion is supported by the Mulliken charge results, which show that the positive Mulliken charges are mainly located on the N atoms of the -NO2 groups (II2 on parylene C and F: ~0.11 e; II4 on parylene AF4 and AF8: 0.1–0.12 e), as shown in Figure 8. Thus, the four-electron orbitals of the methyl C-H2 2p2-1s1 (or C-F2 2p2-2p5) groups mainly capture three surface electrons from the -NO2 2p3-2p4 orbitals on FOX-7 by forming sp electron orbitals, and the adsorption angles are 45° based on the long-range half-full bonds (complex of seven-electron orbitals and five-electron orbitals (-NH2 2p3-1s1)), as shown in Figure 5c,d. Because the high-energy -NO2 bonds are passivated on the (001)-oriented parylene surfaces, these parylene materials can be used as protective membranes for FOX-7, which is in contrast to the (101)-oriented membranes. Furthermore, the (001)-oriented parylene AF8 is a potential protective membrane for FOX-7 because the binding energies of -NO2-benzene and -NO2-methyl are the lowest at ~0 eV and 2.5 eV, respectively.

3.3. Nitro-to-Nitrite Rearrangement of the Decomposition Products of FOX-7 on the Parylene Membranes

According to Rashkeev et al. [12], FOX-7 molecules decompose via C-NO2 bond fission and nitro-to-nitrite rearrangement, and the oxygen atoms in the -NO2 groups of FOX-7 partially or completely dissociate to produce (NH2)(NO2)C=C(NH2)O and ONO. Under continuous aging by temperature [4], pressure [5], and pH [6], the -N-O- and -C-O- bonds in the (NH2)(NO2)C=C(NH2)O product rearrange in the same manner as those in FOX-7. The oxygen atoms in the -NO2 groups of (a) (NH2)(NO2)C=C(NH2)O partially or completely dissociate to give (b) NH(NH2)C=CO and (c) ONO. To analyze the active adsorption sites on the parylene membranes, the probable configurations of the decomposition products of FOX-7 were determined. The probable configurations of the decomposition products on the active adsorption sites (A: methyl (M) and toluene (T) groups) of the parylene membranes are presented as (a1) (NH2)(NO2)C=C(NH2)O and (a2) A-(NH2)(NO2)C=C(NH2)OH in Figure 9; (b1) NH(NH2)C=CO, (b2) A-NH(NH2)C=CO and (b3) (NH2)2C=CO in Figure 10; and (c1) HON=O and (c2) O=NO-A in Figure 11. Considering the effects of the active adsorption sites (A: methyl and toluene groups) and orientations of the parylene membranes, the most likely dissociative adsorption processes for the decomposition products of FOX-7 were analyzed based on the lowest binding energy principle. The -NH- and -O=N-O- bonds in the decomposition products have unsaturated orbitals and thus high electrostatic potentials that attract the surface electrons of the exposed toluene and methyl groups of the (101)- and (001)-oriented parylene membranes [11], as illustrated in Figure 12. Thus, active electrons are transferred from the highest occupied molecular orbitals (HOMOs) on the -NH- and -O=N-O- groups of the decomposition products to the lowest unoccupied molecular orbitals (LUMOs) on the active adsorption sites on the parylene surfaces [13]. Therefore, the high-energy -NO2 bonds in the decomposition products of FOX-7 are passivated by the toluene and methyl groups, which is consistent with the results for FOX-7.
ELF calculations were used to provide bond information to distinguish between the exposed toluene and methyl groups. The adsorption behaviors involve both -NO2–π and -NH2–π interactions, which result in (NH2)(NO2)C=C(NH2)O adsorption at the toluene and methyl groups on the (101)- and (001)-oriented parylene surfaces by H-bond (0.15~0.2 nm) fission, as shown in Figure 13a. The secondary dissociative products of (NH2)(NO2)C=C(NH2)O (a), including NH(NH2)C=CO (b) and HON=O (c), might also adsorb at the basic toluene and methyl adsorption sites. On the (101)-oriented parylene surfaces, new short-range intermolecular H-bonds (0.21–0.23 nm) involving both the -NH–π and -O=N-O–π (or -NH-methyl and -O=N-O-methyl) hybridized orbitals (b2) form due to electrostatic interactions. This result can be explained by the fact that the unsaturated, hybridized N-H bonds in the -C-NH and -C-NH2 groups in (b) have half-full p4 (2p3-1s1) orbitals that can interact with the half-full p4 (2p2-1s12 or 2p2-2p52) orbitals of the hybridized -C-H2 or -C-F2 bonds, as shown in Figure 13b. Simultaneously, both the -Cl (or -F) p7 orbitals on parylene and the complex p7 orbitals (-N-O- 2p3-2p4 orbitals of the -ONO product) compete to couple with the surface benzene π electrons on the (101)-oriented parylene C and F surfaces, as shown in Figure 13c. The p7 (-N-O- 2p3-2p4) orbital captures an active electron from the hybridized benzene π orbital via -C-F2 (or the perfluorobenzene ring) in parylene AF4 (or AF8), creating a hole to offset the electron loss from the -C-F2 bond (or perfluorobenzene ring). The corresponding -NO–π and -NH–π bond energies (3.46 and 1.42 e, respectively) are lower than the methyl-ON- bond energies (5.48 and 1.83 e).
In contrast to the observed interactions with the toluene π groups, the half-full hybridized methyl (-C-F2) p4 (2p2-2p52) orbitals on the (001)-oriented membranes mainly interact with the half-full p4 (N-H 2p3-1s1) and approximately half-full p3 (O-N=O 2p4-2p3-2p4) orbitals on the NH(NH2)C=CO and ONO products, respectively. These results are supported by the low energy results for both the -NO–π and -NH–π interactions (V), as shown in Figure 13. The half-full p4 and approximately half-full p3 orbitals of the -NH2 and -NO2 bonds, respectively, bind to the active electron–hole sites on the toluene π groups. It should be noted that the positive Mulliken charges (0.09–0.145 e) on the N atoms become negative ((NH2)(NO2)C=C(NH2)O: −1.31–0.25 e; NH(NH2)C=CO: −0.34–0.04 e; ONO: −0.36–0.02 e), as shown in the Supporting Information (Figure S1). The N 2p3 electrons are transferred from the positively charged -NO 2p3-2p4 orbitals to the negatively charged -NH 2p3-1s1 (or neighboring O 2p4) orbitals through the electronic transitions of the intramolecular -N-C-N- (or -O-N=O) bonds. Consequently, the enhanced -NH orbital strength induced by the (001)-oriented parylene surfaces can delay the nucleation of the two neighboring decomposition products of FOX-7 and inhibit the nitro-to-nitrite rearrangement by hindering intermolecular and intramolecular hydrogen transfer [7]. Therefore, the (001)-oriented parylene AF8 membrane, which has the lowest bond energy, is proposed to be a potential protective membrane for inhibiting FOX-7 decomposition pathways.

4. Summary and Conclusions

In conclusion, the dissociative adsorption of FOX-7 on (101)- and (001)-oriented crystalline parylene C, F, AF4, and AF8 membranes was demonstrated. Based on the data provided by this study, it can be concluded that FOX-7 adsorbs on the oriented parylene membranes via -NO2–π electrostatic interactions, which appear to be the dominant passivation mechanism of FOX-7. FOX-7 molecules lie flat on the (101)-oriented parylene membranes due to both -NO2–π and -NH2–π electrostatic interactions, and the -NO2 groups mainly bind to the exposed toluene (or methyl) groups on the (001)-oriented parylene membranes. FOX-7 dissociates into (NH2)(NO2)C=C(NH2)O, NH(NH2)C=CO, and ONO via C-NO2 bond fission and nitro-to-nitrite rearrangement. The high-energy -NO2 bonds of (NH2)(NO2)C=C(NH2)O and ONO are passivated by the exposed toluene and methyl groups. In particular, the positively charged -NO 2p3-2p4 orbitals become negatively charged -NH 2p3-1s1 (or neighboring O 2p4) orbitals through intramolecular -N-C-N- (or -O-N=O) electronic transitions on the (001)-oriented membranes. Based on the ELF and Mulliken charge analysis results, the (001)-oriented parylene AF8 membrane is proposed to be a potential protective membrane for passivating the high-energy -NO2 bonds in FOX-7 by a dissociative adsorption process.
The results of this work are insufficient for analyzing the effects of oriented parylene membranes on FOX-7 nucleation and thermolysis mitigation. Further investigations will be focused on the complexation of two FOX-7 molecules and their decomposition products under long-term storage or extreme conditions, such as widely fluctuating temperatures or high pressures. However, this study provides useful information for determining the dissociative processes of FOX-7 on crystalline parylene surfaces and understanding the interactions between explosive products and protective membranes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/polym16030438/s1. Figure S1: Mulliken charges of total N atoms of FOX-7’s decomposition products on (101) and (001) orientated parylene membranes.

Author Contributions

Conceptualization, L.B. and F.D.; methodology, L.B.; software, L.B., J.N. and W.L.; validation, J.N. and W.L.; formal analysis, J.N. and W.L.; investigation, J.N., W.L. and J.Y.; resources, W.L.; data curation, W.L. and J.N.; writing—original draft preparation, W.L., J.N., L.B. and F.D.; writing—review and editing, L.B. and F.D.; visualization, W.L.; supervision, L.B. and F.D.; project administration, L.B. and F.D.; funding acquisition, L.B. and F.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the 973 Project (2014CB8460003), the One Thousand Talents Scheme of Sichuan Province, the National Natural Science Foundation of China (61504168), Hebei Outstanding Young Scholars, Science and Technology Program of Hebei Province (D2016403064, 16044601Z and 15211121), and the Longshan fund of the Southwest University of Science and Technology (17QR004).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, Liang Bian, upon reasonable request.

Acknowledgments

The authors acknowledge the financial support provided by the 973 Project (2014CB8460003), the One Thousand Talents Scheme of Sichuan Province, the National Natural Science Foundation of China (61504168), Hebei Outstanding Young Scholars, Science and Technology Program of Hebei Province (D2016403064, 16044601Z and 15211121), and the Longshan fund of the Southwest University of Science and Technology (17QR004).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rose, A.; Zhu, Z.; Madigan, C.F.; Swager, T.M.; Bulovi, V. Sensitivity gains in chemosensing by lasing action in organic polymers. Nature 2005, 434, 876–879. [Google Scholar] [CrossRef] [PubMed]
  2. Anniyappan, M.; Talawar, M.B.; Venugopalan, S.B.; Gandhe, R. Synthesis Method for preparation of fine TATB (2–5 μm) and its evaluation in plastic bonded explosive (PBX) formulations. J. Hazard. Mater. 2006, 137, 812–819. [Google Scholar] [CrossRef] [PubMed]
  3. Xu, K.; Song, J.; Zhao, F.; Ma, H.; Gao, H.; Chang, C.; Ren, Y.; Hu, R. Thermal behavior, specific heat capacity and adiabatic time-to-explosion of G(FOX-7). J. Hazard. Mater. 2008, 158, 333–339. [Google Scholar] [CrossRef] [PubMed]
  4. Evers, J.; Klapötke, T.M.; Mayer, P.; Oehlinger, G.; Welch, J. alpha and beta-FOX-7, polymorphs of a high energy density material, studied by X-ray single crystal and powder investigations in the temperature range from 200 to 423 K. Inorg. Chem. 2006, 45, 4996–5007. [Google Scholar] [CrossRef]
  5. Zhao, J.; Liu, H. High-pressure behavior of crystalline FOX-7 by density functional theory calculations. Comput. Mater. Sci. 2008, 42, 689–703. [Google Scholar] [CrossRef]
  6. Vo, T.T.; Parrish, D.A.; Shreeve, J.M. 1,1-Diamino-2,2-dintroethene (FOX-7) in copper and nickel diamine complexes and copper FOX-7. Inorg. Chem. 2012, 51, 1963–1968. [Google Scholar] [CrossRef]
  7. Booth, R.S.; Lam, C.S.; Brynteson, M.D.; Wang, L.; Butler, L.J. Elucidating the Decomposition Mechanism of Energetic Materials with Geminal Dinitro Groups Using 2-Bromo-2-nitropropane Photodissociation. J. Chem. Phys. A 2013, 117, 9531–9547. [Google Scholar] [CrossRef]
  8. Bishop, M.M.; Velisavljevic, N.; Chellappa, R.; Vohra, Y.K. High Pressure–Temperature Phase Diagram of 1,1-Diamino-2,2-dinitroethylene (FOX-7). J. Phys. Chem. A 2015, 119, 9739–9747. [Google Scholar] [CrossRef]
  9. Yuan, X.T.; Yang, H.S.; Wang, L.L. Effects of aestivation on the energy budget of sea cucumber Apostichopus japonicus (Selenka) (Echinodermata: Holothuroidea). Acta Ecol. Sin. 2008, 27, 3155–3161. [Google Scholar] [CrossRef]
  10. Hu, A.; Larade, B.; Rachid, H.A.; Lussier, L.S.; Guo, H. Theoretical Prediction of Heats of Sublimation of Energetic Materials Using Pseudo-Atomic Orbital Density Functional Theory Calculations. Propellants Explos. Pyrotech. 2007, 32, 331–337. [Google Scholar] [CrossRef]
  11. Sorescu, D.C.; Boatz, J.A.; Thompson, D.L. First-principles calculations of the adsorption of nitromethane and 1,1-diamino-2,2-dinitroethylene (FOX-7) molecules on the Al(111) surface. J. Phys. Chem. B 2004, 107, 8953–8964. [Google Scholar] [CrossRef]
  12. Rashkeev, S.N.; Kuklja, M.M.; Zerilli, F.J. Electronic excitations and decomposition of 1,1-diamino-2,2-dinitroethylene. Appl. Phys. Lett. 2003, 82, 1371–1373. [Google Scholar] [CrossRef]
  13. Yan, A.H.; Son, S.F.; Jackson, T.L.; Venugopal, P. Validation of numerical simulations for nano-aluminum composite solid propellants. J. Propuls. Power 2011, 27, 1280–1287. [Google Scholar] [CrossRef]
  14. Sonmez, G.; Shen, C.K.F.; Rubin, Y.; Wudl, F. A Red, Green, and Blue (RGB) Polymeric Electrochromic Device (PECD): The Dawning of the PECD Era. Angew. Chem. J. Ger. Chem. Soc. 2005, 43, 1489–1502. [Google Scholar] [CrossRef]
  15. Tompa, A.S. The ultrastructure and mineralogy of the dart from Philomycus carolinianus (Pulmonata: Gastropoda) with a brief survey of the occurrence of darts in land snails. Veliger 1980, 23, 35–42. [Google Scholar]
  16. Castorina, T.C.; Smetana, A.F. Effect of polymer coating on ammonium nitrate substrate. J. Appl. Polym. Sci. 2010, 18, 1373–1383. [Google Scholar] [CrossRef]
  17. Brown, K.E.; Monroe, D.C.; Brown, G.W. Effect of pellet coatings on PETN porosity and slapper detonator efficacy. AIP Conf. Proc. 2017, 1793, 040011. [Google Scholar] [CrossRef]
  18. Dai, X.; Xu, J.; Wen, Y.; Li, Y.; Huang, F.; Li, M.; Zeng, Q. Delay Mechanism of β→δ Phase Transition of Cyclotetramethylene Tetranitramine in Polymer Bonded Explosive Formulations by Heat Conduction Obstacle. Propellants Explos. Pyrotech. 2016, 41, 637–640. [Google Scholar] [CrossRef]
  19. Smalara, K.; Giełdo, A.; Bobrowski, M.; Rybicki, J.; Czaplewski, C. Theoretical study of polymerization mechanism of p-xylylene based polymers. J. Phys. Chem. A 2010, 114, 4296–4303. [Google Scholar] [CrossRef]
  20. Bae, J.; Suh, M.J.; Kim, B.S. Optimal production of 7,10-dihydroxy-8(E)-hexadecenoic acid from palmitoleic acid by Pseudomonas aeruginosa PR3. New Biotechnol. 2010, 27, 352–357. [Google Scholar] [CrossRef]
  21. Bian, L.; Shu, Y.; Wang, X. Computational investigation on the new high energy density material of aluminum enriched 1,1-diamino-2,2-dinitroethylene. J. Mol. Model. 2013, 19, 131–138. [Google Scholar] [CrossRef]
  22. Kahouli, A.; Sylvestre, A.; Pairis, S. Effect of ClH aromatic substitution on structural and dielectric properties of poly(p-xylylene). Polymer 2012, 53, 3001–3007. [Google Scholar] [CrossRef]
  23. Mattos, E.C.; Moreira, E.D.; Diniz, M.F. Characterization of Polymer-Coated RDX and HMX Particles. Propellants Explos. Pyrotech. 2008, 33, 44–50. [Google Scholar] [CrossRef]
  24. Zhao, L.; Gou, Y.; Li, H. A Network Pharmacology Approach to Determine Active Compounds and Action Mechanisms of Ge-Gen-Qin-Lian Decoction for Treatment of Type 2 Diabetes. J. Evid.-Based Complement. Altern. Med. eCAM 2014, 2014, 495840. [Google Scholar] [CrossRef]
  25. Streltsov, D.R.; Mailyan, K.A.; Gusev, A.V. Structure and optical properties of thin poly(p-xylylene)—Silver nanocomposite films prepared by low-temperature vapor deposition polymerization. Polymer 2015, 71, 60–69. [Google Scholar] [CrossRef]
  26. Frac, I.; Kucinska, M.; Gawrys, P. Ambipolar organic thin film transistors prepared with a one step solution technique. Synth. Met. 2016, 220, 197–201. [Google Scholar] [CrossRef]
  27. Sroka-Bartnicka, A.; Olejniczak, S.; Ciesielski, W. Solid State NMR Study and Density Functional Theory (DFT) Calculations of Structure and Dynamics of Poly(p-xylylenes). J. Phys. Chem. B 2009, 113, 5464–5472. [Google Scholar] [CrossRef] [PubMed]
  28. Yan, Q.; Zeman, S.; Sánchez Jiménez, P.E. The effect of polymer matrices on the thermal hazard properties of RDX-based PBXs by using model-free and combined kinetic analysis. J. Phys. Chem. C 2014, 271, 185–195. [Google Scholar] [CrossRef] [PubMed]
  29. Bian, L.; Shu, Y.; Wang, X. A molecular dynamics study on permeabilityof gases through parylene AF8 membranes. Polym. Adv. Technol. 2012, 23, 1429–1528. [Google Scholar] [CrossRef]
  30. Bian, L.; Shu, Y.; Xu, J. Computational Investigation on Adsorption and Diffusion of 4 Gas Molecules in High Pressure PPX Structure. Integr. Ferroelectr. 2011, 127, 83–90. [Google Scholar] [CrossRef]
  31. Fox, T.G.; Flory, P.J. Second-Order Transition Temperatures and Related Properties of Polystyrene. I. Influence of Molecular Weight. J. Appl. Phys. 1950, 21, 581–591. [Google Scholar] [CrossRef]
  32. Bian, L.; Shu, Y.; Li, H. Sorption and permeation of gaseous molecules in amorphous and crystalline PPX C membranes: Molecular dynamics and grand canonical Monte Carlo simulation studies. Chin. Phys. B 2012, 7, 074208. [Google Scholar] [CrossRef]
  33. Sorescu, D.C.; Boatz, J.A.; Thompson, D.L. Classical and Quantum-Mechanical Studies of Crystalline FOX-7 (1,1-Diamino-2,2-dinitroethylene). J. Phys. Chem. A 2001, 105, 5010–5021. [Google Scholar] [CrossRef]
  34. Taylor, D.E.; Rob, F.; Rice, B.M. A molecular dynamics study of 1,1-diamino-2,2-dinitroethylene (FOX-7) crystal using a symmetry adapted perturbation theory-based intermolecular force field. Phys. Chem. Chem. Phys. 2011, 13, 16629–16636. [Google Scholar] [CrossRef]
  35. Isoda, S.; Tsuji, M.; Ohara, M. Structural analysis of β-form poly(p-xylyene) starting from a high-resolution image. Polymers 1983, 24, 1155–1161. [Google Scholar] [CrossRef]
  36. Chung, Y.; Murmann, B.; Selvarasah, S. Low-voltage and short-channel pentacene field-effect transistors with top-contact geometry using parylene-C shadow masks. Appl. Phys. Lett. 2010, 96, 133306. [Google Scholar] [CrossRef]
  37. Chang, T.Y.; Yadav, V.G.; Sarah, D.L. Cell and protein compatibility of parylene-C surfaces. Langmuir 2007, 23, 11718–11725. [Google Scholar] [CrossRef]
  38. He, X.; Zhang, F.; Wang, R. Preparation of a carbon nanotube/carbon fiber multi-scale reinforcement by grafting multi-walled carbon nanotubes onto the fibers. Carbon 2007, 45, 2559–2563. [Google Scholar] [CrossRef]
  39. Hanyaloglu, B.; Aydinli, A.; Oye, M. Low dielectric constant Parylene-F-like films for intermetal dielectric applications. Appl. Phys. Lett. 1999, 74, 606–608. [Google Scholar] [CrossRef]
  40. Diaham, S.; Bechara, M.; Locatelli, M.L. Influence of crystallization-induced amorphous phase confinement on α- and β-relaxation molecular mobility in parylene F. J. Appl. Phys. 2011, 110, 063703. [Google Scholar] [CrossRef]
  41. Sutcliffe, R.; Lee, W.W.; Gaynor, J.F. Characterization and aluminum metallization of a parylene AF-4 surface. Appl. Surf. Sci. 1998, 126, 43–56. [Google Scholar] [CrossRef]
  42. Sabri, S.S.; Levesque, P.L.; Aguirre, C.M. Graphene field effect transistors with parylene gate dielectric. Appl. Phys. Lett. 2009, 95, 242104. [Google Scholar] [CrossRef]
  43. Zhang, J.X.; Zhang, L.D.; Xu, W. Surface plasmon polaritons: Physics and applications. J. Phys. D Appl. Phys. 2012, 45, 113001. [Google Scholar] [CrossRef]
Polymers 16 00438 g001
Figure 1. Chain polymerization configurations (a), chemical formulas (b), potential energy and crystalline models (c) of parylene membranes. The color code is the following: grey, C; green, Cl; blue, F and white H. Blue area of parylene chains, low energy area; and yellow area, high energy area.
Figure 1. Chain polymerization configurations (a), chemical formulas (b), potential energy and crystalline models (c) of parylene membranes. The color code is the following: grey, C; green, Cl; blue, F and white H. Blue area of parylene chains, low energy area; and yellow area, high energy area.
Polymers 16 00438 g001
Polymers 16 00438 g002
Figure 2. Illustrations of frontier orbitals (Highest Occupied Molecular Orbital: HOMO; Lowest Unoccupied Molecular Orbital: LUMO) of FOX-7, (101)- and (001)-oriented parylene F surfaces. The color code is the following: grey, C; green, Cl; blue, F and white H. Blue area of parylene chains, low energy area; and yellow area, high energy area.
Figure 2. Illustrations of frontier orbitals (Highest Occupied Molecular Orbital: HOMO; Lowest Unoccupied Molecular Orbital: LUMO) of FOX-7, (101)- and (001)-oriented parylene F surfaces. The color code is the following: grey, C; green, Cl; blue, F and white H. Blue area of parylene chains, low energy area; and yellow area, high energy area.
Polymers 16 00438 g002
Polymers 16 00438 g003
Figure 3. Surface potential energies (eV) of parylene membranes and FOX-7.
Figure 3. Surface potential energies (eV) of parylene membranes and FOX-7.
Polymers 16 00438 g003
Polymers 16 00438 g004
Figure 4. (I1–I4) Possible adsorption configurations of FOX-7 molecules on the (101)-oriented parylene membranes. (II1–II4) Possible adsorption configurations of FOX-7 molecules on the (001)-oriented parylene membranes.
Figure 4. (I1–I4) Possible adsorption configurations of FOX-7 molecules on the (101)-oriented parylene membranes. (II1–II4) Possible adsorption configurations of FOX-7 molecules on the (001)-oriented parylene membranes.
Polymers 16 00438 g004
Polymers 16 00438 g005
Figure 5. (a,c) The -NO2 groups on FOX-7 bind to the groups in parylene; (b,d) The corresponding energy differences (ΔE) between the adsorption models and single systems.
Figure 5. (a,c) The -NO2 groups on FOX-7 bind to the groups in parylene; (b,d) The corresponding energy differences (ΔE) between the adsorption models and single systems.
Polymers 16 00438 g005
Polymers 16 00438 g006
Figure 6. Illustrations of frontier orbitals of FOX-7 on (101)- and (001)-oriented parylene F membranes. (I1–I4) Possible adsorption configurations of FOX-7 molecules on the (101)-oriented parylene membranes. (II1–II4) Possible adsorption configurations of FOX-7 molecules on the (001)-oriented parylene membranes.
Figure 6. Illustrations of frontier orbitals of FOX-7 on (101)- and (001)-oriented parylene F membranes. (I1–I4) Possible adsorption configurations of FOX-7 molecules on the (101)-oriented parylene membranes. (II1–II4) Possible adsorption configurations of FOX-7 molecules on the (001)-oriented parylene membranes.
Polymers 16 00438 g006
Polymers 16 00438 g007
Figure 7. Binding energies (eV) of various probable configurations of FOX-7 adsorbed on (101)- and (001)-oriented parylene membranes. In the figure, “M” refers to the methyl group, “T” refers to the toluene group, and “O” and “H” refer to the -NO2 and -NH2 functional groups of FOX-7, respectively.
Figure 7. Binding energies (eV) of various probable configurations of FOX-7 adsorbed on (101)- and (001)-oriented parylene membranes. In the figure, “M” refers to the methyl group, “T” refers to the toluene group, and “O” and “H” refer to the -NO2 and -NH2 functional groups of FOX-7, respectively.
Polymers 16 00438 g007
Polymers 16 00438 g008
Figure 8. Mulliken charges of total N atoms of FOX-7 on (101)- and (001)-oriented parylene membranes.
Figure 8. Mulliken charges of total N atoms of FOX-7 on (101)- and (001)-oriented parylene membranes.
Polymers 16 00438 g008
Polymers 16 00438 g009
Figure 9. The most probable absorption configurations of (NH2)(NO2)C=C(NH2)O and A-(NH2)(NO2)C=C(NH2)OH on (101)- and (001)-oriented parylene membranes. In the figure, “A” refers to active adsorption sites of parylene surfaces. (a1,a1′) (NH2)(NO2)C=C(NH2)O; (a2,a2′) A-(NH2)(NO2)C=C(NH2)OH.
Figure 9. The most probable absorption configurations of (NH2)(NO2)C=C(NH2)O and A-(NH2)(NO2)C=C(NH2)OH on (101)- and (001)-oriented parylene membranes. In the figure, “A” refers to active adsorption sites of parylene surfaces. (a1,a1′) (NH2)(NO2)C=C(NH2)O; (a2,a2′) A-(NH2)(NO2)C=C(NH2)OH.
Polymers 16 00438 g009
Polymers 16 00438 g010
Figure 10. The most probable absorption configurations of NH(NH2)C=CO, A-NH(NH2)C=CO, and (NH2)2C=CO on (101)- and (001)-oriented parylene membranes. (b1,b1′) NH(NH2)C=CO; (b2,b2′) A-NH(NH2)C=CO; (b3, b3′) (NH2)2C=CO.
Figure 10. The most probable absorption configurations of NH(NH2)C=CO, A-NH(NH2)C=CO, and (NH2)2C=CO on (101)- and (001)-oriented parylene membranes. (b1,b1′) NH(NH2)C=CO; (b2,b2′) A-NH(NH2)C=CO; (b3, b3′) (NH2)2C=CO.
Polymers 16 00438 g010
Polymers 16 00438 g011
Figure 11. The most probable absorption configurations of HON=O and O=NO-A (FOX-7) on (101)- and (001)-oriented parylene membranes. (c1,c1′) HON=O; (c2,c2′) O=NO-A.
Figure 11. The most probable absorption configurations of HON=O and O=NO-A (FOX-7) on (101)- and (001)-oriented parylene membranes. (c1,c1′) HON=O; (c2,c2′) O=NO-A.
Polymers 16 00438 g011
Polymers 16 00438 g012
Figure 12. Illustrations of frontier orbitals of FOX-7’s decomposition products on (101)- and (001)-oriented parylene F membranes. (a2,a2′) A-(NH2)(NO2)C=C(NH2)OH; (b1′) NH(NH2)C=CO; (b2,b2′) A-NH(NH2)C=CO; (c1,c1′) HON=O; (c2,c2′) O=NO-A.
Figure 12. Illustrations of frontier orbitals of FOX-7’s decomposition products on (101)- and (001)-oriented parylene F membranes. (a2,a2′) A-(NH2)(NO2)C=C(NH2)OH; (b1′) NH(NH2)C=CO; (b2,b2′) A-NH(NH2)C=CO; (c1,c1′) HON=O; (c2,c2′) O=NO-A.
Polymers 16 00438 g012
Polymers 16 00438 g013
Figure 13. Binding energies (eV) of various probable configurations of FOX-7’s decomposition products adsorbed on (101)- and (001)-oriented parylene membranes. (a–c) reflect the adsorption models of parylene-(NH2)(NO2)C=C(NH2)OH, parylene-NH(NH2)C=CO, and O=NO-parylene, respectively. In the figure, “V” refers to the adsorption behavior of both -NO–π and -NH–π interactions of FOX-7’s decomposition products being bound by parylene toluene groups through intramolecular H-bond breakage, “M” refers to methyl-ON-(or -NH-) adsorption, and “T” refers to toluene-ON-(or -NH-) adsorption.
Figure 13. Binding energies (eV) of various probable configurations of FOX-7’s decomposition products adsorbed on (101)- and (001)-oriented parylene membranes. (a–c) reflect the adsorption models of parylene-(NH2)(NO2)C=C(NH2)OH, parylene-NH(NH2)C=CO, and O=NO-parylene, respectively. In the figure, “V” refers to the adsorption behavior of both -NO–π and -NH–π interactions of FOX-7’s decomposition products being bound by parylene toluene groups through intramolecular H-bond breakage, “M” refers to methyl-ON-(or -NH-) adsorption, and “T” refers to toluene-ON-(or -NH-) adsorption.
Polymers 16 00438 g013
Table 1. Calculated structures, densities, and glass transition temperatures (Tg) of parylene membranes. “a–c” show the lengths of parylene cells.
Table 1. Calculated structures, densities, and glass transition temperatures (Tg) of parylene membranes. “a–c” show the lengths of parylene cells.
Parylene C
[36,37,38]		Parylene F
[39,40]		Parylene AF4
[41]		Parylene AF8
[21]
a (nm)0.950.881.151.17
b (nm)2.152.151.952.13
c (nm)4.214.643.723.85
Density (g·cm−3)1.29 (1.27)1.13 (1.11)1.71 (1.58)2.17 (1.93)
Tg (K)332 (347)312 (289)840 (773)610 (673)
The experimental and computer data are tested by SEM and MD results, as shown in Refs. [21,35,36,37,38,39,40,41].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, W.; Bian, L.; Dong, F.; Nie, J.; Yang, J. MD-DFT Calculations on Dissociative Absorption Configurations of FOX-7 on (001)- and (101)-Oriented Crystalline Parylene Protective Membranes. Polymers 2024, 16, 438. https://doi.org/10.3390/polym16030438

AMA Style

Luo W, Bian L, Dong F, Nie J, Yang J. MD-DFT Calculations on Dissociative Absorption Configurations of FOX-7 on (001)- and (101)-Oriented Crystalline Parylene Protective Membranes. Polymers. 2024; 16(3):438. https://doi.org/10.3390/polym16030438

Chicago/Turabian Style

Luo, Weihui, Liang Bian, Faqin Dong, Jianan Nie, and Jingjie Yang. 2024. "MD-DFT Calculations on Dissociative Absorption Configurations of FOX-7 on (001)- and (101)-Oriented Crystalline Parylene Protective Membranes" Polymers 16, no. 3: 438. https://doi.org/10.3390/polym16030438

APA Style

Luo, W., Bian, L., Dong, F., Nie, J., & Yang, J. (2024). MD-DFT Calculations on Dissociative Absorption Configurations of FOX-7 on (001)- and (101)-Oriented Crystalline Parylene Protective Membranes. Polymers, 16(3), 438. https://doi.org/10.3390/polym16030438

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