Polyethylene C2H4 (PE) is a flexible plastic existing in an amorphous state and widely used in consumer products. In particular, PE plays an important role in material science, which is considered as a bright candidate for new industrial materials. When PE is combined with wood, it creates an environmentally Wood Plastic Composite material (WPC)   such as reducing energy consumption in production, lightweight, and sound insulation   ; reducing pollution and greenhouse effect; enhancing biodegradability   . Also, WPC is widely used in industries such as automotive, construction, et al.  . The flexible plastic resins are made of Polyethylene (PE) , Polypropylen , Polyvinyl Clorua (PVC)  and polystyrene (PS) . The Wood components include wood pulp, cotton, wheat straw, bagasse, and rice husk     . PE plastic is used in practice including both High-Density Polyethylene (HDPE) and low-density WPC . PE, as a binder and wood pulp as an additive; HDPE includes material hydrophilic wood pulp and PE hydrophobic wood pulp  . Therefore, PE is considered as an important additive component in industrial applications. PE is divided into many types, mainly based on density, monomer, flexibility , and copolymer   . In particular, with Low-Density Polyethylene (LDPE) which is the most commonly used and commercially produced at high pressure (P), P = 2400 Bar, temperature (T) from T = 363 K to T = 383 K .
To study this material, it is good to use the experimental, theoretical, and simulation methods. As for the experimental method, using pressing methods and multi-layer extrusion method with low cost has high durability, ductility  and increasing stress   as Jatin et al. , Kurz et al. , Pouriayevali et al.  for that plastic deformation of PE is an isotropic function with pressure putting on (the deformation (ε) in the experiment which is always less than ε < 0.12); Epee et al.  suggested that the deformation of the polymer with ε < 800 s−1; Argon et al.  suggested that plastic deformation is due to the twisted bonding pairs along the polymer chain; Roberson suggested that the shear stress is caused by changing angles and movements of molecules ; Eyring et al.  suggested that the plastic deformation is caused by the shear stress, structural changing, and binding energy  - . With the simulation method, Deng et al.  suggested that plastic deformation is caused by the local structure. Besides, Maeda, Takeuchi , Srolovitz et al.  successfully used the molecular statistical method (MS) to study the deformation of three-dimensional metal glasses; Theodorou and Suter   successfully simulated the material in the glass polymer material and studied the deformation. After that, MS. Mott et al. , Hutnick et al.  used this method to study the plastic deformation of Polypropylene and Polycarbonate. While Mott et al. suggested that the plastic deformation is the result of changing the displacement in atoms or molecular branches, Hutnick et al. suggested that the movement of atoms does not depend on chemical reactions. Similarly, Brown, Clarke , Mckechnie, and Clarke  successfully investigated the effect of temperature on the bond in glass polyethylene and mechanics of materials. It is said that increasing the plastic deformation and the length of the initial link is necessary. Although there have been many studies on the plastic deformation of PE in the static state or the dynamic state   -  such as Baltsas et al. , Haefele et al. , using simulation method with Low-Density Polyethylene (LDPE) at high pressure; Asteasuain et al.  use Graphical Optimization Tool (gOPT) of general PROcess Modelling System (gPROMS) simulator program to optimize LDPE; Bezzo et al.  successfully used Fluent, gPROMS for liquids by calculating molecular dynamics. Recently, Clarke  has successfully performed uniaxial deformation of amorphous Polymers with different deformation levels at low temperatures, which is obtained in qualitative form. Capaldi et al.  suggested that the compressive deformation of Polymer is the same as the change of angle, the angle shift along the chain; Li et al.  successfully performed Single-axis plastic deformation of Polyethylene (PE) amorphous by Monte Carlo method (MC) with FCC structure and bonding length (r), r = 1.53 Å. The result shows that there is a dependence on temperature (T), and the heating rate and Uniaxial Tension of the material ; Ospina et al.  using the MC simulation result of the initial plastic deformation stage of polyethylene gives the consistent result with the experimental result, and it does not reduce the plastic deformation   ; factors that alter the structure of PE   as concentrations of impurities playing an important role in compounds  .
The result shows that the phase transition of Polyethylene depends on the temperature    as the liquefaction process which is just below the room temperature T = 300 K, limited by the movement of atoms at transition layer between the crystal area and the amorphous region. This shift occurs very weakly in HDPE with temperature from T = 123 K to T = 173 K and it is linked with the movement of CH2 groups attached to C2H4 , the change in the molecular shape of the polymer depends on the relationship between temperature and pressure ; the phase transition depends on the heating rate, the total energy of the system during the deformation process . The results show that the glass transition temperature (Tg) depends on the movement of atoms with valuable in approx from Tg = 133 K to Tm = 408 K , and this is performed on experimental measurements .
Besides, many authors have successfully studied the plastic deformation of PE by the z-axis stretching method. The results show the influence of the chain length, the number of chains, the strain rate, and the temperature. This depends on the stress-strain     which is ended at the source of plastic deformation. Polyethylene is a problem that has not been explained in detail . To solve the problem, we focus on studying the effect of atomic number, temperature, annealing time on the structure, and the plastic deformation of PE.
2. Method of Calculation
Initially, randomly sow atomic number (N), N = 2000 atoms, 4000 atoms, 6000 atoms, 8000 atoms, 10,000 atoms Polyethylene (C2H4 or PE) into the cube by the Molecular Dynamics (MD) method     Dreading pair interaction (1), cyclic boundary conditions       through the total energy of the system (Etot): Etot = Ebond + Eangle + Edihedral + Enon-bonding (1).
With: Ebond is the bond energy, Eangle is the bond angle energy, the Edihedral is dihedral energy, Enon-bonding is van der Waals energy in the Lennard-Jones interaction, Etot is the total energy of the system, Kb = 350 kcal/mol, Kθ = 60 kcal/mol∙rad2 is the stiffness coefficient, the bond angle coefficient, r0 = 1.53 Å is the bond length, θ0 = 1.911 rad (109.5˚) is the link angle, C0 = 1.736, C1 = −4.490, C2 = 0.776, C3 = 6.99 (kcal/mol) are the coefficients, σ = 4.01 Å is the energy at 0 eV, ε = 0.112 kcal/mol is the dielectric constant, rc = 10Å is the radius interrupt.
After obtaining, all PE samples for running 106 steps molecular dynamics (MD) simulation recovery statistics at temperature (T), T = 500 K; 106 steps NPT (atomic number, pressure, and constant temperature) MD simulation at T = 500 K. After obtaining PE samples at T = 500 K, the samples were lowered from T = 500 K to T = 100 K. Particularly with N = 10000 atoms at T = 500 K, the temperature is lowered to T = 120 K, 100 K, 80 K, 60 K, 40 K. When increasing t, from t = 0 ps to t = 50 ps, 100 ps, 150 ps, 200 ps at T = 100 K. The temperature was set intentionally, which is to study PE material in the crystalline state with temperature the below glass transition temperature (Tg), Tg = 250 K , size (l) of PE material lying in the range from l = 3.73 nm to l = 6.36 nm, the time for each simulation step is Δt = 0.1 fs. The program code used to run the PE sample that is the open-source program code LAMMPS   with the heating rate, which is obeyed following the Nosé-Hoover rule   , van der Waal link     . To study the plastic deformation of PE under the effect of the force in the direction of the z-axis with the constant heating rate at pressure (P), P = 0 GPa by defining quantities in the direction of the z-axis single tension from the NPT equations of the PE shift process . The stress components are calculated from energy contributions related to the length of the link, bonding angle, bipolar angle, and non-bonding interaction, which are combined with the Common Neighbor Analysis (CNA) method    by Ovito software.
3. Results and Discussion
3.1. Effect of Atomic Number
3.1.1. The Phase Transition of Polyethylene
The result of the phase transition Polyethylene C2H4 (PE) is shown in Figure 1.
The results show that after running the Recovery Statistics (RS), RS = 106 steps of moving the Molecular Dynamics (MD) simulation at the temperature (T), T = 500 K, then the pressure of the system (P) decreases, with PE 2000 atoms lead to P decreases from P = 4464 Bar to P = −144 Bar. At a steady state, when N increasing from N = 2000 atoms to N = 10,000 atoms leads to P increases from P = −144 Bar to P = 292 Bar (Figure 1(a)). Similarly, with the process of running 106 steps NPT (with atoms number, pressures, and temperatures is constant) leads to P changes on the range from P = −558 Bar to P = 706 Bar, with PE 10,000 atoms the P has the smallest change (Figure 1(b)). The process of temperature reduction from T = 500 K to T = 100 K shows that the P of PE 10,000 atoms has not changed significantly (Figure 1(c)). The energy values of PE at T = 100 K has value corresponding change. The total energy pair (Epair) decreases from Epair = −3014 eV to Epair = −15,161 eV (Figure 1(d)), the total electron bond energy (Ebond) increases from Ebond = 199 eV to Ebond = 994 eV (Figure 1(e)), the total angle bond energy (Eangle) increases from Eangle = 223eV to Eangle = 1158 eV (Figure 1(f)), the total dihedral energy (Edihed) increases from Edihed = 719 eV to Edihed = 3544 eV (Figure 1(g)) and Etot decreases from Etot = −1873 eV to Etot = −9465 eV (Figure 1(h)). Basing on the given result, after 106 steps moving of recovering statistics, then PE reached equilibrium; 106 steps moving NPT at temperature (T), T = 500 K lead to PE existed in a liquid state. When the temperature decreases from T = 500 K to T = 100 K leading to PE changes from the liquid state to a new crystalline state, running stably 106 steps NPT at T = 100 K obtained PE in the new crystalline state corresponds to the shape of PE (Figure 1(i)). As a result, when increasing atoms number (N) leads to the size (l) increases, but Etot decreases.
3.1.2. Plastic Deformation Process of Polyethylene
The obtained result shows that after the process of the temperature reduction from T = 500 K down to T = 100 K with 106 step MD simulation, the
(a) (b) (c) (d) (e) (f) (g) (h) (i)
Figure 1. The phase transition of polyethylene after 106 steps recovery statistics at the temperature at 500 K (a), 106 steps of NPT at the temperature 500 K (b), 106 steps to the lower temperature from 500 K down to 100 K (c), the energy values (d)-(h), and the shape of the PE (i) with atomic number, the number of different displacement steps.
Table 1. The process of transforming the number of structural units after the number of different MD shift steps.
(a) (b) (c) (d)
Figure 2. The structural unit number shapes of Polyethylene as FCC structure (a), HCP structure (b), BCC structure (c), Amor structure (d).
polyethylene has a cube shape (Figure 2(a)-(d)) with the number of structural units is 4 FCC, 56 HCP, 2 BCC, 1938 Amor (Table 1). When increasing moving the number of steps from 0 to 106 of the MD simulation leading to PE changing from the cube shape to the rectangular box shape (Figures 3(a)-(f)), and the of number structural units of FCC, HCP, BCC, Amor changes corresponding (Table 1). When increasing moving the number of steps of the MD leading to the number of structural units unchanged. Particularly, the shape has a huge change raising the question after compressing the PE time of the z-axis, the PE’s shape has changed leading to changing the plastic deformation characteristics of the PE. To confirm that studying the plastic deformation characteristics with different moving the number of steps, the results shown in Figure 4, Table 2.
The results show that the number of structural units with the atomic number (N), N = 2000 atoms is 12 FCC, 60 HCP, 6 BCC, and the characteristics of the plastic deformation include the bonding energy between atoms, Ebond which is shown by the red color line, Eangle is angular bond energy shown by the blue color line, the Edihed is energy dihedral shown by the blue color line, Etotal shown by the purple color line, (Figure 4(a), Table 2). When increasing N from N = 2000 atoms to N = 4000 atoms, 6000 atoms, 8000 atoms and 10,000 atoms leads to the total energy pair (Epair) decreases from Epair = −2756 eV to Epair = −13,866 eV (Figure 4(b)), energy bond (Ebond) increases from Ebond = 163 eV to Ebond = 812 eV (Figure 4(c)), angular bond energy Eangle increases from Eangle = 256 eV to Eangle = 1264 eV (Figure 4(d)), the energy dihedral increases from Edihed = 751 eV to Edihed = 3899 eV (Figure 4(e)), and Etotal decreases from Etotal = −1586 eV to Etotal = −7891 eV (Figure 4(f)), free volume change the number of atoms (Figure 4(g)), the plastic deformation process of PE (Figure 4(h)), the PE shape (Figure 4(i)) after deformation at T = 100 K, and the number of structural units of FCC
(a) (b) (c) (d) (e) (f)
Figure 3. The PE structure shape of PE with the different number of MD shift steps: In the initial state (a), after 2 × 105 steps MD (b), 4 × 105 steps MD (c), 6 × 105 steps MD (d), 8 × 105 steps MD (e), 106 steps MD (f).
(a) (b) (c) (d) (e) (f) (g) (h) (i)
Figure 4. The plastic deformation process as the energy E (a), the total energy pair Epair (b), the total electron bond energy Ebond (c), the total energy angle Eangle (d), the total dihedral energy Edihed (e), the total energy Etot (f), free volume change the number of atoms (g), stress (h), the PE shape (i) after deformation at T = 100 K with atomic number and different moving the number of steps.
Table 2. The process changes the number of structural units with the different number of MD moving steps.
increases from 12 FCC to 43 FCC, HCP increases from 60 HCP to 265 HCP, BCC increases from 6 BCC to 28 BCC, finishes increasing the number of structural units of FCC, BCC, which is not significant, only HCP increases rapidly (Table 2). Consequently, when increasing the atoms number and the number of shift steps of MD leads to a rapid increase in shape and the number of structural units of FCC, and the plastic deformation energy of PE also varies greatly. The result of the PE modeling process is consistent with the result     and coincide with the previous authors who used the Monte Carlo (MC) method    with a PE transition temperature is 103 K  . Additionally, when increasing the number of the atoms leads to the angle among the atoms is a constant corresponding to 109.5˚, the length of the link r increases from r = 1.529 Å to r = 1.558 Å, the obtained result is correctly consistent with the result of bonding angle of 109.27˚, and the link length r = 1.529 Å  , the size of the PE material increase from l = 3.73 nm to l = 6.63 nm, the total energy Etotal decreases from Etotal = −1586 eV to Etotal = −7891 eV. This means that when increasing the number the atoms leads to increase the link length, Etotal decreases and changes the plastic deformation characteristics of PE.
3.1.3. The Effect of Annealing Time
Similarly, the effect of annealing time of PE material with N = 10,000 atoms at T = 100 K, the results are shown in Figures 5-7.
The results show that with t = 50 ps, the shape (Figure 6(a)), the number of structural units is 54 FCC, 280 HCP, 35 BCC (Figure 7(a)), the plastic deformation characteristics of PE (Figure 5(a)). When increasing the annealing time from t = 50 ps to t = 100 ps, 150 ps, 200 ps lead to FCC, HCP, BCC decreases and Amor increases accordingly: with t = 50 ps, there are 54 FCC, 280 HCP, 35 BCC, 9581 Amor; t = 100 ps has 43 FCC, 265 HCP, 28 BCC, 9564 Amor; t = 150 ps has 51 FCC, 268 HCP, 23 BCC, 9508 Amor; t = 200 ps has 42 FCC, 278 HCP, 19 BCC, 9461 Amor (Figures 7(b)-(d)), shape (Figures 6(b)-(d)), and the plastic deformation characteristics (Figures 5(b)-(d)). Also, Etot increased (Figure 5(e)), free volume change the number of atoms (Figure 5(f)). This showed that after the annealing time, then FCC, HCP, BCC decreased, Amor, increased after finishing which did not make a significant change in the structure.
(a) (b) (c) (d) (e) (f)
Figure 5. The plastic deformation process of PE at temperature T = 100 K with t = 50 ps (a), 100 ps (b), 150 ps (c), 200 ps (d), total energy of system (e), free volume change number atoms (f).
(a) (b) (c) (d)
Figure 6. The shape of PE material after different annealing time as t = 50 ps (a), 100 ps (b), 150 ps (c), 200 ps (d).
(a) (b) (c) (d)
Figure 7. The number of structural units of FCC, HCP, BCC, Amor after the different annealing time: t = 50 ps (a), 100 ps (b), 150 ps (c), 200 ps (d).
3.2. Influence of Temperature
The result shows that at the temperature (T), T = 40 K in which the number of structural units 47 FCC, 307 HCP, 18 BCC (Figure 9(a)) corresponding to the plastic deformation characteristics of PE (Figure 8). When increasing the temperature from T = 40 K to T = 60 K, 80 K, 100 K, 120 K, then the number of structural units of FCC, HCP, BCC decreases, and Amor increases accordingly with T = 40 K, has 47 FCC, 307 HCP, 18 BCC, 9628 Amor; T = 60 K, which has 53 FCC, 267 HCP, 33 BCC, 9647 Amor; T = 80 K has 56 FCC, 295 HCP, 30 BCC, 9619 Amor; T = 100 K has 43 FCC, 265 HCP, 28 BCC, 9664 Amor; T = 120 K has 39 FCC, 303 HCP, 18 BCC, 9640 Amor (Figures 9(b)-(e)), the shape and plastic deformation characteristics of PE (Figures 8(a)-(f)). In addition, the Etot increases (Figure 5(e)), which shows that when increasing temperature leads to the plastic deformation increases, and the structure is unchangeable. The obtained result shows that when increasing N leads to l increasing, Etot decreases, and the number of structural units increases. When increasing the time of elongation and temperature, l increases, Etot increases, and the number of structural units of the system decrease. This result is the basis for future empirical research.
After studying the effect of the atomic number (N), temperature (T), and annealing time (t) on the structure and the plastic deformation of the polyethylene, the result shows that when increasing N, from N = 2000 atoms to N = 10,000 atoms, the moving number of MD leads to the shape and the number of structure units FCC increase. The angle among the atoms is a constant corresponding to 109.50. The length of the link r increases from r = 1.529 Å to r = 1.558 Å. The plastic deformation energy of PE has an enormous change as the bonding angle of 109.270. The length of the link r = 1.529 Å, and the size of the PE material increases from l = 3.73 nm to l = 6.63 nm. The total energy (Etotal) decreases from Etotal = −1586 eV to Etotal = −7891 eV, and it coincides with the previous authors who used the Monte Carlo (MC) method with a PE transition temperature is 103 K. When N increasing leads to the length of the link increases, Etotal decreases, FCC, HCP, BCC, Amor increases, and changes the plastic deformation characteristics of PE with an increase in T and the t. The obtained results are very significant for future experimental research as studying the effect of the number of structural units on the structure, plastic deformation, conductivity, magnetism of PE materials.
We thank computer room, Faculty of Physics of Hanoi National University of Education has created all favorable conditions, help us throughout the process of calculating, simulating and finishing the content of the article.
(a) (b) (c) (d) (e) (f)
Figure 8. The plastic deformation process of PE 10,000 atoms with the total energy pair Epair (a), the total electron bond energy Ebond (b), the total energy angle Eangle (c), the total dihedral energy Edihed (d), the total energy Etot (e), the free volume change the number of atoms (f) at the temperature T = 100 K with different temperatures.
(a) (b) (c) (d) (e)
Figure 9. The number of structural units of FCC, HCP, BCC, Amor after different temperature as T = 40 K (a), 60 K (b), 80 K (c), 100 K (d), 120 K (e).
 Sewda, K. and Maiti, S.N. (2013) Dynamic Mechanical Properties of High Density Polyethylene and Teak Wood Flour Composites. Polymer Bulletin, 70, 2657-2674.
 Joshi, S.V., Drzal, L.T., Mohanty, A.K. and Arora, S. (2004) Are Natural Fiber Composites Environmentally Superior to Glass Fiber Reinforced Composites? Composite Part A, 35, 371-376.
 Orhan, Y. and Büyükgüngör, H. (2000) Enhancement of Biodegradability of Disposable Polyethylene in Controlled Biological Soil. International Biodeterioration & Biodegradation, 45, 49-55.
 Maiti, S.N. and Hassan, M.R. (1989) Melt Rheological Properties of Polypropylene-Wood Flour Composites. Journal of Applied Polymer Science, 37, 2019-2032.
 Li, Q. and Matuana, L.M. (2003) Surface of Cellulosic Materials Modified with Functionalized Polyethylene Coupling Agents. Journal of Applied Polymer Science, 88, 278-286.
 Rizvi, G., Matuana, L.M. and Park, C.B. (2000) Foaming of PS/Wood Fiber Composites Using Moisture as a Blowing Agent. Polymer Engineering & Science, 40, 2124-2132.
 Patil, Y.P., Gajre, B., Dusane, D., Chavan, S. and Mishra, S. (2000) Effect of Maleic Anhydride Treatment on Steam and Water Absorption of Wood Polymer Composites Prepared from Wheat Straw, Cane Bagasse, and Teak Wood Sawdust Using Novolac as Matrix. Journal of Applied Polymer Science, 77, 2963-2967.
 Rowell, R.M., Sanadi, A.R., Caulfield, D.F., Jacobson, E., Leo, A.L., Carvalho, F.X. and Frollini, E. (1997) Utilization of Natural Fibers in Plastic Composites: Problems and Opportunities. Forest, 2, 23-51.
 Zhang, M.Q., Rong, M.Z. and Lu, X. (2005) Fully Biodegradable Natural Fiber Composites from Renewable Resources: All-Plant Fiber Composite. Composites Science and Technology, 65, 2514-2525.
 Saheb, D.N. and Jog, J.P. (1999) Natural Fiber Polymer Composites: A Review. Polymers for Advanced Technologies, 18, 351-363.
 Balasuriya, P.W., Ye, L. and Mai, Y.W. (2001) Mechanical Properties of Wood Flake-Polyethylene Composites. Part I: Effects of Processing Methods and Matrix Melt Flow Behaviour. Compos A, 32, 619-629.
 Stark, N.M. and Matuana, L.M. (2003) Ultraviolet Weathering of Photostabilized Wood-Flour Filled High-Density Polyethylene Composites. Journal of Applied Polymer Science, 90, 2609-2617.
 Stark, N.M. and Matuana, L.M. (2004) Surface Chemistry Changes of Weathered HDPE/Wood-Flour Composites Studied by XPS and FTIR Spectroscopy. Polymer Degradation and Stability, 86, 1-9.
 Kim, D.-M. and Iedema, P.D. (2004) Molecular Weight Distribution in Low-Density Polyethylene Polymerization; Impact of Scission Mechanisms in the Case of a Tubular Reactor. Chemical Engineering Science, 59, 2039-2052.
 Pladis, P. and Kiparissides, C. (1998) A Comprehensive Model for the Calculation of Molecular Weight + Long-Chain Branching Distribution in Free-Radical Polymerizations. Chemical Engineering Science, 53, 3315-3333.
 Cervantes, A.M., Tonelli, S., Brandolin, A., Bandoni, J.A. and Biegler, L.T. (2002) Large-Scale Dynamic Optimization for Grade Transitions in a Low Density Polyethylene Plant. Computers & Chemical Engineering, 26, 227-237.
 Sun, L.Y., Gibson, R.F., Gordaninejad, F. and Suhr, J. (2009) Energy Absorption Capability of Nanocomposites: A Review. Composites Science and Technology, 69, 2392-2409.
 McShane, G.J., Stewart, C., Aronson, M.T., Wadley, H.N.G., Fleck, N.A. and Deshpande, V.S. (2008) Dynamic Rupture of Polymer-Metal Bilayer Plates. International Journal of Solids and Structures, 45, 4407-4426.
 Amini, M.R., Simon, J. and Nemat-Nasser, S. (2010) Numerical Modeling of Effect of Polyurea on Response of Steel Plates to Impulsive Loads in Direct Pressure-Pulse Experiments. Mechanics of Materials, 42, 615-627.
 Jatin, Sudarkodi, V. and Basu, S. (2014) Investigations into the Origins of Plastic Flow and Strain Hardening in Amorphous Glassy Polymers. International Journal of Plasticity, 56, 139-155.
 Kurtz, S.M., Rimnac, C.M., Santner, T.J. and Bartel, D.L. (1996) Exponential Model for the Tensile True Stress-Strain Behavior of As-Irradiated and Oxidatively Degraded Ultra-High Molecular Weight Polyethylene. Journal of Orthopaedic Research, 14, 755-761.
 Pouriayevali, H., Arabnejad, S., Guo, Y.B. and Shim, V.P.W. (2013) A Constitutive Description of the Rate-Sensitive Response of Semi-Crystalline Polymers. International Journal of Impact Engineering, 62, 35-47.
 Epee, A.F., Lauro, F., Bennani, B. and Bourel, B. (2011) Constitutive Model for a Semi-Crystalline Polymer under Dynamic Loading. International Journal of Solids and Structures, 48, 1590-1599.
 Baschnagel, J., Binder, K., Doruker, P., Gusev, A.A., Hahn, O., Kremer, K., et al. (2000) Bridging the Gap between Atomistic and Coarse-Grained Models of Polymers: Status and Perspectives. In: Viscoelasticity, Atomistic Models, Statistical Chemistry, Vol. 152, Springer, Berlin, 41-156.
 Yoon, D.Y., Smith, G.D. and Matsuda, T. (1993) A Comparison of a United Atom and an Explicit Atom Model in Simulations of Polymethylene. The Journal of Chemical Physics, 98, 10037-10043.
 Fukunaga, H., Takimoto, J. and Doi, M. (2002) A Coarse-Graining Procedure for Flexible Polymer Chains with Bonded and Nonbonded Interactions. The Journal of Chemical Physics, 116, 8183-8190.
 Müller-Plathe, F. (2002) Coarse-Graining in Polymer Simulation: From the Atomistic to the Mesoscopic Scale and Back. ChemPhysChem, 3, 754-769.
 Uhlherr, A. and Theodorou, D.N. (1998) Hierarchical Simulation Approach to Structure and Dynamics of Polymers. Current Opinion in Solid State & Materials Science, 3, 544-551.
 Deng, D., Argon, A.S. and Yip, S. (1989) Simulation of Plastic Deformation in a Two-Dimensional Atomic Glass by Molecular Dynamics IV. Philosophical Transactions of the Royal Society of London, 329, 613.
 Hutnik, M., Argon, A.S. and Suter, U. (1993) Simulation of Elastic and Plastic Response in the Glassy Polycarbonate of 4,4’-Isopropylidenediphenol. Macromolecules, 26, 1097.
 Depa, P.K. and Maranas, J.K. (2005) Speed Up of Dynamic Observables in Coarse-Grained Molecular-Dynamics Simulations of Unentangled Polymers. The Journal of Chemical Physics, 123, Article ID: 094901.
 Depa, P.K. and Maranas, J.K. (2007) Dynamic Evolution in Coarse-Grained Molecular Dynamics Simulations of Polyethylene Melts. The Journal of Chemical Physics, 126, Article ID: 054903.
 Lavine, M.S., Waheed, N. and Rutledge, G.C. (2003) Molecular Dynamics Simulation of Orientation and Crystallization of Polyethylene during Uniaxial Extension. Polymer, 44, 1771-1779.
 Pearson, D.S., Ver Strate, G., Von Meerwall, E. and Schilling, F.C. (1987) Viscosity and Self-Diffusion Coefficient of Linear Polyethylene. Macromolecules, 20, 1133-1141.
 Yang, L., Srolovitz, D.J. and Yee, A.F. (1999) Molecular Dynamics Study of Isobaric and Isochoric Glass Transitions in a Model Amorphous Polymer. The Journal of Chemical Physics, 110, 7058-7069.
 Takeuchi, H. and Roe, R.J. (1991) Molecular Dynamics Simulation of Local Chain Motion in Bulk Amorphous Polymers. II. Dynamics at Glass Transition. The Journal of Chemical Physics, 94, 7458-7465.
 Baltsas, A., Papadopoulos, E. and Kiparissides, C. (1998) Application and Validation of the Pseudo-Kinetic Rate Constant Method to High Pressure LDPE Tubular Reactors. Computers & Chemical Engineering, 22, S95-S102.
 Haefele, M., Kienle, A., Boll, M. and Schmidt, C.-U. (2006) Modeling and Analysis of a Plant for the Production of Low Density Polyethylene. Computers & Chemical Engineering, 31, 51-65.
 Asteasuain, M., Tonelli, S.M., Brandolin, A. and Bandoni, J.A. (2001) Dynamic Simulation and Optimisation of Tubular Polymerisation Reactors in gPROMS. Computers & Chemical Engineering, 25, 509-515.
 Bezzo, F., Macchietto, S. and Pantelides, C.C. (2000) A General Framework for the Integration of Computational Fluid Dynamics and Process Simulation. Computers & Chemical Engineering, 24, 653-658.
 Li, J., Mulder, T., Vorselaars, B., Lyulin Alexey, V. and Michels, M.A. (2006) Monte Carlo Simulation of Uniaxial Tension of an Amorphous Polyethylene-Like Polymer Glass. Macromolecules, 39, 7774-7782.
 Li, J., Mulder, T., Vorselaars, B., Lyulin, A.V. and Michels, M.A.J. (2006) Monte Carlo Simulation of Uniaxial Tension of an Amorphous Polyethylene-Like Polymer Glass. Macromolecules, 39, 7774-7782.
 Khonakdar, H.A., Morshedian, J., Wagenknecht, U. and Jafari, S.H. (2003) An Investigation of Chemical Crosslinking Effect on Properties of High Density Polyethylene. Polymer, 44, 4301-4309.
 Goertzen, W.K. and Kessler, M.R. (2007) Dynamic Mechanical Analysis of Carbon/Epoxy Composites for Structural Pipeline Repair. Composites Part B: Engineering, 38, 1-9.
 Lin, Z., Lin, Z., Zhen, S., Kui, N., Ami, M.I., Ayw, Y.W. and Xei, H.K. (1995) Electrical insulating materials. International Symposium on Electrical Insulating Materials, Tokyo, 17-20 September 1995, 173-176.
 Bartczak, Z., Argon, A.S. and Cohen, R.E. (1992) Deformation Mechanisms and Plastic Resistance in Single-Crystal-Textured High-Density Polyethylene. Macromolecules, 25, 5036-5053
 Hossain, D., Tschopp, M.A., Ward, D.K., Bouvard, J.L., Wang, P. and Horstemeyer, M.F. (2010) Molecular Dynamics Simulations of Deformation Mechanisms of Amorphous Polyethylene. Polymer, 51, 6071-6083.
 Tschopp, M.A., Ward, D.K., Bouvard, J.L. and Horstemeyer, M.F. (2011) Atomic Scale Deformation Mechanisms of Amorphous Polyethylene under Tensile Loading. TMS Conference Proceedings, Vol. 2, 1-6.
 Bouvard, J.L., Ward, D., Hossain, D., Marin, E.B. and Horstemeyer, M.F. (2009) Review of Hierarchical Multiscale Modeling to Describe the Mechanical Behavior of Amorphous Polymers. JEMT, Special Issue: Predictive Science & Technology in Mechanics & Materials, 131, Article ID: 041206.
 Gee, R.H. and Boyd, R.H. (1998) The Role of the Torsional Potential in Relaxation Dynamics: A Molecular Dynamics Study of Polyethylene. Computational and Theoretical Polymer Science, 8, 93-98.
 Mayo, S.L., Olafson, B.D. and Goddard III, W.A. (1990) DREIDING: A Generic Force Field for Molecular Simulations. The Journal of Physical Chemistry, 94, 8897-8909.
 Capaldi, F.M., Boyce, M.C. and Rutledge, G.C. (2001) Enhanced Mobility Accompanies the Active Deformation of a Glassy Amorphous Polymer. Physical Review Letters, 89, Article ID: 175505.
 Tschopp, M.A., Bouvard, J.L., Ward, D.K., Bammann, D.J. and Horstemeyer, M.F. (2013) Influence of Ensemble Boundary Conditions (Thermostat and Barostat) on the Deformation of Amorphous Polyethylene by Molecular Dynamics.
 Shepherd, J.E., McDowell, D.L. and Jacob, K.I. (2006) Modeling Morphology Evolution and Mechanical Behavior during Thermo-Mechanical Processing of Semi- Crystalline Polymers. Journal of the Mechanics and Physics of Solids, 54, 467-489.
 Li, Q. and Matuana, L.M. (2003) Effectiveness of Maleated and Acrylic Acid-Functionalized Polyolefin Coupling Agents for HDPE-Wood-Flour Composites. Journal of Thermoplastic Composite Materials, 16, 551-564.
 Lu, J.Z., Wu, Q. and Negulescu, H. (2000) The Influence of Maleation on Polymer Adsorption and Fixation, Wood Surface Wettability, and Interfacial Bonding Strength in Wood-PVC Composites. Wood and Fiber Science, 34, 434-459.
 Wang, Y.E.H., Yeh, F.C., Lai, S.M., Chan, H.C. and Shen, H.F. (2003) Effectiveness of Functionalized Polyolefins as Compatibilizers for Polyethylene/Wood Flour Composites. Polymer Engineering & Science, 43, 933-945.
 Balasuriya, P.W., Ye, L., Mai, Y.W. and Wu, J. (2002) Mechanical Properties of Wood Flake-Polyethylene Composites II. Interface Modification. Journal of Applied Polymer Science, 83, 2505-2521.
 Tsuzuki, H., Branicio, P.S. and Rino, J.P. (2007) Structural Characterization of Deformed Crystals by Analysis of Common Atomic Neighborhood. Computer Physics Communications, 177, 518-523.
 Sankaranarayanan, S.K.R.S., Bhethanabotla, V.R. and Joseph, B. (2005) Molecular Dynamics Simulation Study of the Melting of Pd-Pt Nanoclusters. Physical Review B, 71, Article ID: 195415.
 Honeycutt, J.D. and Andersen, H.C. (1987) Molecular Dynamics Study of Melting and Freezing of Small Lennard-Jones Clusters. The Journal of Physical Chemistry, 91, 4950-4963.