Investigation of Projectile Impact Behaviors of Graphene Aerogel Using Molecular Dynamics Simulations

Xinyu Zhang ,Wenjie Xia ,Yang Wang ,Liang Wang,⋆ and Xiaofeng Liu

1Department of Engineering Mechanics,Shanghai Jiao Tong University,Shanghai,200240,China

2Department of Aerospace Engineering,Iowa State University,Ames,IA,50011,USA

3School of Materials Science and Engineering,University of Science&Technology Beijing,Beijing,100083,China

4Zernike Institute for Advanced Materials,University of Groningen,Groningen,9747,The Netherlands

ABSTRACT Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics (MD) simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency (EAE) metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.

KEYWORDS Graphene aerogel;molecular dynamics simulation;impact response;energy absorption

1 Introduction

Graphene is a two-dimensional material with numerous outstanding properties,such as a large surface area[1],excellent thermal[2],and electrical conductivity[3],and a high elastic modulus[4].These characteristics show great potential in engineering applications.In order to enhance material properties,extensive investigations have converted two-dimensional graphene materials into threedimensional structures,creating novel materials through delicate designs.Notably,Gao et al.[5]pioneered the use of freeze-drying techniques in 2013 to fabricate graphene aerogel (GA).This unprecedentedly lightweight solid material has garnered considerable attention from the scientific community.For example,compared to other lightweight 3D graphene networks[6,7],GA possesses significant potential in nanoparticle capture owing to its much lower density.However,observing GA’s mechanical response to high-velocity impacts by nanoparticles presents considerable difficulties,thereby imposing challenges in conducting comprehensive investigations into GA’s projectile impact behaviors.

The synthesis of GA,along with its chemical and conductive characteristics,has been extensively explored in previous studies [8–11].Results show that three-dimensional GA can fully exploit graphene’s inherent properties compared to two-dimensional materials due toπ-πstacking interactions within graphene layers.For example,in energy storage applications,GA has been widely used in electrode fabrication due to properties such as high surface area,high porosity,high specific capacitance,and high cyclic stability [12–14].Furthermore,GA exhibits considerable adsorption capabilities,making it suitable for adsorbing spilled oil and other organic pollutants in environmental conservation applications [15].Although there has been much research on GA,its mechanical properties have rarely been studied in detail[16].

Molecular dynamics (MD) simulation provides a novel approach for in-depth investigations of the mechanical properties and behaviors of two-dimensional materials.Recently,there has been an increase in research focused on using MD simulations to analyze GA.Qin et al.[17] employed MD simulations to study the 3D GA system,emphasizing the investigation of its mechanical characteristics,including tensile strength,compression strength,buckling modulus,and Young’s modulus.The findings demonstrate that GA exhibits superior mechanical properties compared to most polymeric cellular materials.Patil et al.[18] conducted a study on the fracture behaviors of GAs subjected to tensile loading with pre-existing cracks.The investigation focused on the impact of aspect ratio on fracture strength and toughness.The findings revealed a decrease in fracture strength with increasing aspect ratio,while the fracture toughness remained constant across all densities.Moreover,Patil et al.[19]also employed the reverse-ballistic method to investigate GA’s shock wave propagation behaviors.The research findings show a direct proportionality between the shock and particle velocities,while the slope of the Shock-Hugoniot curve exhibits a discernible correlation with GA density.Specifically,as GA density increases,the slope of the Shock-Hugoniot curve also increases.Despite these noteworthy contributions,the mechanical performance of GA under projectile impact loadings has not been reported.

In this work,MD simulation is performed to investigate the mechanical performance of GA subjected to nano-projectile impacts.The remainder of this article is organized such that Section 2 presents the construction of the GA model with various densities and the MD simulation details.Section 3 briefly introduces the penetration equations for low-density materials.Section 4 thoroughly investigates the mechanical performances of GA under projectile impacts,revealing the energy absorption mechanism and capacity of GAs.Finally,Section 5 concludes the paper with some remarks.

2 Method and Materials

In this study,molecular dynamics simulations were performed using the large-scale atomic/molecular massively parallel simulator(LAMMPS)[20].The OVITO software[21]was utilized to visualize the simulation results and obtain snapshots during the simulations.

2.1 Interatomic Potentials

In molecular dynamics simulations,selecting a reasonable interatomic potential function that accurately replicates the system’s mechanical performance is paramount.This work employed the adaptive intermolecular reactive empirical bond order (AIREBO) potential [22,23] to describe the interatomic interactions between carbon atoms of graphene sheets.This potential is a well-established force field extensively utilized in MD simulations.It has been validated as reliable for simulating the mechanical performances of graphene materials,as observed in many experimental studies[24,25].

In the modeling process,the total potential energy of the system is given by:

The first term,EC-C,is the interatomic potential of carbons,which is defined as follows:

whereEijREBOis theREBOterm inAIREBOpotential,which denotes the short-range interactions that govern the potential energy ofC-Ccovalent bonds;EijLJis theLJterm that characterizes the intermolecular interactions within the system;EkijlTORTIONis theTORTIONterm which represents theTORSIONterm that relies on the dihedral angles.

Furthermore,this study adopted the pairwise interaction potentialEC-inclusionto describe the interactions between carbon atoms and inclusion particles and the pairwise interaction potentialEinclusion-inclusionto describe the interactions between two inclusion particles.The 12-6 Lennard-Jones potential was adopted forEC-inclusionandEinclusion-inclusion:

whereεijrepresents the depth of the potential well,σijis the zero potential distance;rijdenotes the interparticle distance.Note that the GA system with different densities can be obtained by changing the values ofσijinEinclusion-inclusion.In this study,parametersεij=0.65 eV andσij=10–14 Å were adopted.This results in GAs with a density range of 275–546 kg/m3.Moreover,the cutoff distance for theLJterm inAIREBOpotential was modified in this study.Specifically,the cutoff distance isrd1=1.7 Å in the modeling process andrd2=10.2 Å in the impact simulations.

2.2 Modeling Process

The all-atom GA models with different densities were constructed to investigate the mechanical response under projectile impact loadings.The GA models consist of 750 randomly distributed graphene sheets and 750 inclusion particles within the simulation box,as shown in Fig.1a.The initial density of the model is 3.9 kg/m3,closely resembling air density.Inclusions were introduced to control the porous microstructure,mimicking the effects of water clusters in the laboratory synthesis of GA using the freeze-casting method[26].Periodic boundary conditions were applied in all three directions of the simulation box during the construction of the MD model.

According to the GA modeling procedures reported in previous work[17–19],the model was first equilibrated under the NPT ensemble(constant number of atoms N,pressure P,and temperature T),with a constant temperature of 300 K and pressure of 1 atm.Thereafter,a program consisting of multiple pretreatment cycles was implemented to achieve a stable GA model,as illustrated in Fig.1d.Each cycle comprised four steps:(i)first increasing the pressure of the system from 1 to 1000 atm at a constant temperature of 300 K under the NPT ensemble;(ii)then increasing the system’s temperature from 300 to 2000 K;(iii)maintaining a constant temperature of 2000 K;(iv)reducing the temperature back to 300 K.The last three steps were carried out under the NVT ensemble (constant number of atoms N,volume V,and temperature T).Stable GA models with different densities were successfully obtained by repeating this cycle eight times,with each step lasting for a simulation time of 25 ps.

Figure 1:Schematic representation of the all-atom model creation process:(a)initial model in which water inclusions and graphene sheets are randomly distributed;(b) equilibrated graphene-inclusion structure after eight cycles of NPT-NVT;(c)GA model with a density of 546 kg/m3 after removing the water inclusions;(d)Pressure p,temperature T as a function of simulation time during the modeling process

2.3 Simulation Details

In this study,nanoscale diamond projectiles were used to impact the GA with different densities.The 12–6 Lennard-Jones potential was selected to represent the interaction between the projectile and the graphene sheets,with the parametersεij=0.035 eV andσij=3.46 Å[27].The simulations were conducted in the NVE ensemble (constant number of atoms N,volume V,and energy E),and the projectile was initially positioned 2 nm away from the upper surface of the GA model.The model dimensions are given in Table 1.The impact region was defined as a cylindrical area with a radius ofR=L/2 Å,while the boundary region remained fixed during impact,as shown in Fig.2.The diamond projectile was used with a diameter ofD=2.6 nm and an initial velocity varying from 1–7 km/s.

Table 1:Description of the MD model

Figure 2:(a)Schematic representation of the initial state of the MD model;(b)Top view of the GA model

3 Penetration Equation of Low-Density Material

When the projectile penetrates a low-density target material,it overcomes both static and dynamic resistance.According to works [28,29],the rigid body penetration process can be described by the following equation:

wherempis the mass of the projectile;V(t)is the velocity of the projectile at timetduring penetration;ρtis the density of GA;Γis the drag coefficient;Sis the cross-sectional area of the projectile that remains constant.Besides,the scaling factorkis introduced to quantify the relationship between the crushing strengthσcand the GA density,namely:

Eq.(4)can be equivalently expressed by dividing its sides by the cross-sectional areaSas follows:

whererrepresents the radius of the projectile;ρpdenotes the density of the projectile.By applying the chain rule of differentiation and subsequently integrating both sides of Eq.(6),the following expression can be obtained:

After a simple algebraic operation,it finally yields:

4 Results and Discussion

4.1 Mechanical Response under Projectile Impact

Fig.3 depicts the Mises stress distributions of the penetration process with three different GA densities,namelyρ=275,436,and 546 kg/m3.The initial impact velocity of the projectile is 10 km/s.Due to the intrinsic porous structure of GA,the impact of a projectile affects only a localized region of the GA structure,as observed in Fig.3.Notably,the impact-induced penetration of the projectile leads to the collapse of the pore structure,and the graphene sheet undergoes stretching,subsequently leading to the rupture of covalent bonds.This contributes to GA’s energy absorption efficiency and its resistance to projectile impact compared to other porous materials,such as silica aerogel and aluminum foam.Moreover,the observed penetration patterns exhibit substantial differences for the three distinct densities under the same impact velocity,thus demonstrating the pivotal influence of GA density on the impact response.

For GA with low density,the large pores allow for greater motion flexibility of graphene sheets during projectile impact,which in turn leads to significant irreversible deformation of GA.In contrast,high-density GA exhibits smaller pore sizes,leading to reduced flexibility and rotational movement of graphene sheets.

The variations of different energy components,namelyEREBO,ELJ,andETORSION,are illustrated in Fig.4.It is observed that theEREBOpotential term,characterizing the alteration of theC-Ccovalent bond energy,experiences the most substantial fluctuation during the projectile impact process.In contrast,the contribution of the Lennard-Jones potential is minimal.Further assessment of theLJpotential’s impact is achieved by comparing the residual projectile velocity curves with and withoutLJinteratomic interactions,as shown in Fig.5.Notably,theLJpotential energy has a minimal influence on the residual projectile velocity,suggesting that the interfacial interactions between GA sheets have negligible effects in halting the projectile impact.

The normalized penetration depthh/2r vs.V0curves are depicted in Fig.6,wherein the analytical curves are fitted using Eq.(8).It is observed that the normalized penetration depth diminishes with an increase in GA density at a constant impact velocity.This observation signifies that GAs with higher density exhibit superior projectile capture capabilities compared to those with lower density.Furthermore,the simulation data presented in Fig.6 enable the determination of the crushing strength of GA with different densities.Herein,a fitting relationship between the crushing strengthσcand the normalized GA densityρ/ρs(ρsis the density of graphene)is obtained,as demonstrated in Fig.7.It is seen that the strength increases nonlinearly with an increase in GA density.

Figure 3:Schematic representation of the GA models to highlight the mises stress of each atom during impact simulation at the impact velocity of 10 km/s

Figure 4:(a) Penetration depth vs.time curves for different GA densities;(b–d) the variations of ΔEREBO,ΔELJ,and ΔETORSION vs.time for each GA with different densities (The initial velocity of the projectile is 10 km/s)

Figure 5:(Continued)

Figure 6:The normalized penetration depth vs.impact velocity curves

Figure 7:Crushing strength vs.the normalized density curve

4.2 Effect of Temperature on Projectile Impact

For many applications,such as in aerospace engineering,structural materials may encounter extensive temperature variations.As a novel material with potential applications in aerospace engineering,the influence of temperature on projectile impact is a key concern.In order to address this issue,the impact simulation is performed at varying temperatures employing two different GA densities,namely 436 and 546 kg/m3.The initial velocity of the projectile was set at 6 km/s,aligning with observations from laboratory experiments.

The research conducted by Burchell et al.[30] focused on examining the particle-capturing capabilities of silica aerogel at temperatures ranging from 175 to 763 K.The findings indicate that the penetration depth of silica aerogel remains unaffected when the temperature is below 400K,as illustrated in Fig.8b.However,beyond this critical temperature,there is a significant decline in penetration depth.When the temperature reaches 763 K,the penetration depth reduces by about 30% compared to GA.In contrast,results show that ambient temperature has minimal impact on the penetration depth of GA at different densities,as illustrated in Fig.8a.This demonstrates GA’s thermal stability and mechanical reliability in terms of impact loadings across a wide range of temperatures.

Figure 8:Normalized penetration depth vs.temperature for:(a)GA;(b)silica aerogel

4.3 Energy Absorption for Projectile Impact

4.3.1 Kinetic Energy of GA

Fig.9 shows the variations in the total kinetic energy of GA with different densities.During the penetration process,the carbon-carbon covalent bonds undergo elongation and rupture.Observations indicate that GA with higher density exhibits a pronounced advantage in terms of kinetic energy absorption,attributable to the greater number of covalent bonds involved in breaking within the impacted region.Notably,all three density curves exhibit non-smooth behavior,characterized by multiple peaks occurring at different intervals.Prior to the peak,the carbon-carbon covalent bonds within the impacted region mainly undergo extension.Thus,the kinetic energy of the projectile converts to the strain energy of GA.As the penetration continues,numerous covalent bonds break,leading to the conversion of strain energy into the kinetic energy of GA.Since GA is a highly porous structure,this process repeats until full penetration,giving rise to the formation of multiple peaks in the kinetic energy curves.

Figure 9:Variation of total kinetic energy of GA vs.normalized penetration depth for different densities.(The initial velocity of the projectile is 12 km/s)

4.3.2 Energy Absorption Density

Previous studies have demonstrated the potential of carbon-based nanomaterials as highly promising candidates for structural energy absorption owing to their excellent mechanical properties[31].In this work,the energy absorption density(EAD)parameter was studied to gain a better understanding of GA’s energy absorption performance.EAD represents the kinetic energy absorbed per unit volume of GA,defined as follows:

whereΔEkrepresents the kinetic energy dissipation of the projectile.In this section,the projectiles are launched at a velocity ranging from 12 to 20 km/s,resulting in the initial kinetic energy of the projectiles from 2.31×10-15to 6.42×10-15J.

As depicted in Fig.10,the EAD value exhibits a positive correlation with both the impact velocity and GA density.During complete penetration,irreversible deformation occurs within GA,accompanied by the conversion of energy from the projectile to GA.At higher impact velocities,the fracture of carbon covalent bonds in the impact region becomes more pronounced,resulting in a higher capacity for kinetic energy absorption and consequently leading to a larger EAD value compared to that observed under lower impact velocities.Additionally,the EAD value increases with increasing GA density,which is explainable by the fact that GA with higher density possesses a greater density of covalent bonds in the impact region,facilitating more covalent bond breakages and thus absorbing more kinetic energy from the projectile.

4.3.3 Energy Absorption Efficiency

In this subsection,the concept of energy absorption efficiency(EAE)was introduced to further evaluate the energy absorption capacity of GA during projectile impact.The EAE parameterηis defined as follows:

whereE0is the initial kinetic energy of the projectile.According to the energy conservation principle,the dissipated kinetic energy of the projectile is transformed into the kinetic energy(ΔEkinetic)and the potential energy(ΔEpotential)of GA as follows:

The energy absorption efficiency is illustrated in Fig.11 with varying initial velocities of the projectile.It is observed thatηexperiences a decrease as the impact velocity increases in Region I.However,once the projectile velocity surpasses 24 km/s,the value ofηremains relatively unchanged in Region II.This phenomenon is attributable to the fact that as the initial impact velocity rises,bothE0andΔEkconcurrently increase.Nevertheless,the increase inE0surpasses that ofΔEk.In accordance with Eq.(10),this results in a decline in EAE,as depicted in Region I.When the initial velocity of the projectile surpasses the critical threshold of 24 km/s,the higher impact velocity leads to a larger impact region.Consequently,a greater number of covalent bonds break compared to those at lower impact velocities,thereby inducing a significant increase inΔEkand leading to a relatively stable value ofηin Region II instead of its anticipated decline in Region I.

Figure 10:A fitting surface of EAD as a function of GA density and impact velocity (GA density varies from 275 to 546 kg/m3,and impact velocity varies from 12 to 15 km/s)

Figure 11:Energy absorption efficiency under different impact velocities(from 12 to 40 km/s)for GA density 546 kg/m3

5 Conclusions

In this study,MD simulation is used to investigate the mechanical performance of graphene aerogel subjected to nano-projectile impacts.A full-atomic GA model with different densities is first developed,followed by MD simulations to thoroughly investigate the mechanical responses and energy absorption mechanisms under impact loadings.The results show that the porous structure and density are key factors determining the mechanical response of GA.Specifically,the impactinduced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in the graphene sheet.Therefore,theREBOterm in theAIREBOpotential contributes to the main part of the energy absorption during all the impact simulations,while the contribution of theLJandTORSIONterms is minimal.

Moreover,the influence of temperature on GA has been investigated and determined to be insignificant,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Besides,the EAD and EAE parameters are adopted to assess the energy absorption capacity of GA during projectile penetration.Results show that at higher initial projectile velocities,a larger number of covalent bonds involved in the impact process contribute to energy absorption,which results in a relatively stable and higher energy absorption efficiency.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.

Acknowledgement:The authors wish to express their appreciation to the reviewers for their helpful suggestions,which greatly improved the presentation of this paper.

Funding Statement:This work was supported by the National Natural Science Foundation of China(No.12102256).

Author Contributions:Xinyu Zhang:Writing—original draft,Software,Visualization,Data curation.Wenjie Xia: Methodology,Writing—review &editing.Yang Wang: Methodology,Writing—review&editing.Liang Wang: Conceptualization,Methodology,Formal analysis,Supervision,Writing—review&editing,Project administration,Funding acquisition.Xiaofeng Liu:Methodology,Writing—review&editing.

Availability of Data and Materials:All materials and data used in this review are readily accessible to interested readers.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.