Atomistic Investigation of Material Deformation Behavior of Polystyrene in Nanoimprint Lithography

: This research investigates deformation behavior of polystyrene (PS) as a thermoplastic resist material for the thermal nanoimprint lithography (T-NIL) process. Molecular dynamics modeling was conducted on a PS substrate with dimensions 58 × 65 × 61 Å that was imprinted with a rigid, spherical indenter. The e ﬀ ect of indenter size, force, and imprinting duration were evaluated in terms of indentation depth, penetration depth, recovery depth, and recovery percentage of the polymer. The results show that the largest indenter, regardless of force, has the most signiﬁcant impact on deformation behavior. The 40 Å indenter with a 1 µ N of force caused the surface molecules to descend to the lowest point compared to the other indenters. An increase in indenter size resulted in higher penetration depth, recovery depth, and recovery percentage. Higher durations of imprint cycle (400 fs) resulted in plastic deformation of the PS material with minimal recovery (4 Å). The results of this research lay the foundation for explaining the e ﬀ ect of several T-NIL process parameters on virgin PS thermoplastic resist material. nanoimprint lithography


Introduction
Nanoimprint lithography (NIL) is a versatile high precision method to produce nanoscale features by mechanical deformation of an ultra-thin film polymer resist [1]. NIL was developed as a high throughput fabrication technique by Stephen Chou and his research team [2,3]. NIL has been applied to several applications including biological [4], photonic, and semiconductor electronics [5].
Nanoimprint lithography is broadly classified as thermal nanoimprint lithography (T-NIL) and ultra-violet lithography (UV-NIL). In the UV-NIL process, ultra-violet light is used to cure the polymer resist material on a substrate to fabricate nanoscale features using a transparent patterning mold. In this research, our group employs the T-NIL method, wherein a thermoplastic resist is applied onto a rigid substrate and imprinted by a rigid mold [6]. Typically, thermoplastic polymers transform to a glassy phase and aid deformation when heated up to their glass transition temperature (T g ). Different types of mold materials and geometries are implemented depending on the features to be imprinted [7]. Molds are fabricated using focused ion beam (FIB) technology from silicon and silicon dioxide materials. The NIL cycle begins with heating the ultra-thin film thermoplastic polymer to its glass transition temperature (T g ) on a substrate and inserting the mold to replicate its pattern. After complete insertion of the mold, it is held within the polymer for a desired duration followed by cooling to ambient temperature. Further, the mold is retracted from the polymer to retain its negative replica within the ultra-thin film resist material. Materials Studio implements a smart minimizer that combines steepest descent and conjugant gradient (cg) methods. The conjugant gradient is a general method to minimize a function f(x), where f can be any function of points in an N-dimensional space. The conjugant gradient searches the optimal design along the conjugate direction defined by: where β k = ‖∇f(x k )‖ ‖∇f(x k−1 )‖ represents the scale factor, n k is search direction, x k is initial design, f(x k ) is functional value, and ∇f(x k ) is gradient vector. The steepest descent is found using an algorithm that continues in the direction, which will minimize the value of a function, given the current point.
where ( +1) is a point that takes steps of certain length, , at direction ∆ ( ) from the previous point, ( ) . In addition to minimization parameters, specifications about the force field must be identified. In this model, consistent valence force field (CVFF) force field was chosen. CVFF is applicable to many polymer systems based on [36][37][38][39]. CVFF provides for intermolecular interactions that occur within a polymeric MD simulation. The CVFF parameters for PS are shown in Table 1. CVFF is parameterized by the following: Bond stretching energy where ( ) is separation between the atoms constituting the i-th covalent bond, 0 ( ) is equilibrium separation, ( ) represents the angular stiffness of i-th covalent bond, and is the number of covalent bonds.
Angle bending potential where ( ) is the angle between the i-th pair of covalent bonds, with 0 ( ) being equilibrium.
( ) represents the angular stiffness of i-th covalent bond, and is the number of covalent bonds. Torsional rotational potential is given as: where ∅ ( ) is the torsion angle between the planes formed by the first-second and second-third adjacent bond, with ∅ 0 ( ) being equilibrium. Materials Studio implements a smart minimizer that combines steepest descent and conjugant gradient (cg) methods. The conjugant gradient is a general method to minimize a function f(x), where f can be any function of points in an N-dimensional space. The conjugant gradient searches the optimal design along the conjugate direction defined by: where β k = ∇f(x k ) ∇f(x k−1 ) represents the scale factor, n k is search direction, x k is initial design, f x k is functional value, and ∇f x k is gradient vector. The steepest descent is found using an algorithm that continues in the direction, which will minimize the value of a function, given the current point.
where x (k+1) is a point that takes steps of certain length, t k , at direction ∆x (k) from the previous point, x (k) . In addition to minimization parameters, specifications about the force field must be identified. In this model, consistent valence force field (CVFF) force field was chosen. CVFF is applicable to many polymer systems based on [36][37][38][39]. CVFF provides for intermolecular interactions that occur within a polymeric MD simulation. The CVFF parameters for PS are shown in Table 1. CVFF is parameterized by the following: Bond stretching energy where r (i) is separation between the atoms constituting the i-th covalent bond, r (i) 0 is equilibrium separation, K (i) B represents the angular stiffness of i-th covalent bond, and n B is the number of covalent bonds.
Angle bending potential where θ (i) is the angle between the i-th pair of covalent bonds, with θ represents the angular stiffness of i-th covalent bond, and n A is the number of covalent bonds. Torsional rotational potential is given as: where ∅ (i) is the torsion angle between the planes formed by the first-second and second-third adjacent bond, with ∅ (i) T represents the torsional stiffness of the i-th covalent bond, d (i) and n (i) being integers parameterized for specific atomic systems.
Van der Waals pair potential is defined as: where r ij is the distance between the two atoms, and N is the total number of atoms in the system. A ij = A i A j , B ij = B i B j , and A and B are parameters for systems represented by i and j.
Electrostatic pair potential is defined as: where r is the distance between two atoms i and j. q i and q j are the partial charges on the atoms i and j, respectively. 0 is the permittivity of free space. N is the number of atoms in the atomistic systems. The full potential is given by summing all components:

Conversion to LAMMPS Readable Format
Materials Studio outputs data in the .msi format, or alternately, .car and .mdf. The large-scale atomistic massively parallel simulator (LAMMPS) was used as a source code to model the NIL systems [40]. LAMMPS has a built-in package called "msi2lmp" that converts .car and .mdf files to .data files. The molecular configuration used to conduct dynamics for simulations was stored in these files. In order to conduct MD simulation of NIL processes, specific input files must be written to produce the computational representation of all physical phenomena applied to the system. The scripts govern the different phases that represent the MD-NIL process.
Modifying the indenter size and force causes noticeable changes in indentation depth (I D ), penetration depth (P D ), recovery depth (R D ), and recovery percentage (R P ). I D is defined as the amount that a group of surface atoms moves from their initial position at the beginning of the simulation. Only the z dimension of motion was considered for consistency purposes [15,41]. P D differs slightly from I D because penetration is instantaneous while indentation occurs over the entire course of the simulation. Indentation is associated with the molecules displaced, while penetration is identified by the mold. R D is the distance that a group of molecules travels from their lowest point of penetration (L P ) to their final resting positions. Recovery percentage (R P ) is the percentage amount that a group of molecules travels from their position at the indenter's L P to their final resting positions.
R P is defined as follows: where P D and R D are penetration depth and recovery depth, respectively. Figure 2a shows the schematic of the nanoimprint lithography process for polystyrene. Herein the spherical indenter causes plastic deformation of the PS, which results in a permanent set of indenters in the substrate. Figure 2b depicts the comparative cross-section area when comparing cylindrical versus spherical indenters. The hoop stress of both the cylindrical and spherical indenters would exert a force circumferentially (perpendicular both to the axis and radius of the object) in both directions on the polystyrene chains. However, conventional NIL molds with a square profile or linear grating patterns will exert force perpendicular to the mold loading direction. Thus, the stress distribution pattern and resultant polymer deformation behavior are closely dependent on the NIL mold geometry. In this research, we investigate the deformation behavior of the spherical NIL mold based on variations in indenter size and applied force. RP is defined as follows: where PD and RD are penetration depth and recovery depth, respectively. Figure 2a shows the schematic of the nanoimprint lithography process for polystyrene. Herein the spherical indenter causes plastic deformation of the PS, which results in a permanent set of indenters in the substrate. Figure 2b depicts the comparative cross-section area when comparing cylindrical versus spherical indenters. The hoop stress of both the cylindrical and spherical indenters would exert a force circumferentially (perpendicular both to the axis and radius of the object) in both directions on the polystyrene chains. However, conventional NIL molds with a square profile or linear grating patterns will exert force perpendicular to the mold loading direction. Thus, the stress distribution pattern and resultant polymer deformation behavior are closely dependent on the NIL mold geometry. In this research, we investigate the deformation behavior of the spherical NIL mold based on variations in indenter size and applied force. (a)

Importance of Glass Transition Temperature
Choosing a suitable glass transition temperature (T g ) for the polymers is important to produce high quality imprints. T g is identified by plotting the specific volume of the polymer as a function of temperature. When the slope changes, thermal transitions occur [42]. The literature [42] states that the T g of PS can range from 323-373 K. Figure 3 shows that the onset of thermal transitions appeared at about 360 K and continued through about 380 K. This preliminary simulation was essential to be able to determine the T g value for simulations.

Importance of Glass Transition Temperature
Choosing a suitable glass transition temperature (Tg) for the polymers is important to produce high quality imprints. Tg is identified by plotting the specific volume of the polymer as a function of temperature. When the slope changes, thermal transitions occur [42]. The literature [42] states that the Tg of PS can range from 323-373 K. Figure 3 shows that the onset of thermal transitions appeared at about 360 K and continued through about 380 K. This preliminary simulation was essential to be able to determine the Tg value for simulations.

Molecular Dynamics Simulation
In this research, spherical indenters were used to deform the PS resist material. The indenters were treated as a rigid body with no adhesion interaction with the PS resist. Figure 2 presents an overview of the NIL process. The indenter was initiated at a fixed distance away from the resist (Figure 2a). After the temperature was ramped up to the Tg of PS [42], the indenter was lowered into

Molecular Dynamics Simulation
In this research, spherical indenters were used to deform the PS resist material. The indenters were treated as a rigid body with no adhesion interaction with the PS resist. Figure 2 presents an overview of the NIL process. The indenter was initiated at a fixed distance away from the resist (Figure 2a). After the temperature was ramped up to the T g of PS [42], the indenter was lowered into the resist material. The indenter force, size, and initial position are parameters that govern movement of the indenter. The indenter force was varied to evaluate its effect on the deformation behavior of the polymer. Once the indenter had penetrated the resist, it continued its descent until it was held in place for a fixed duration of time. Before the indenter was removed from the resist, the temperature of the entire system was cooled down to room temperature (~300 K) [43].
The constant-energy, constant-volume ensemble (NVE) boundary conditions were applied to simulations with a time step of 0.005 fs. The units assigned to the model were LAMMPS "real." The type of indenter used was modeled after a spherical mold. Figure 3 shows the actual temperature distribution throughout the NIL process. The simulation had an initial stabilization phase of 50,000 time steps. Temperature was ramped up to room temperature (298 K) over 15,000 time steps. Following this initial temperature rise, temperature was ramped up above the T g of PS (373 K) for 50,000 time steps. In the next phase, temperature was kept constant at 373 K for 50,000 time steps. The indenter was lowered into the resist material while holding temperature constant at 373 K for 60,500 time steps. When the indenter reached the lowest point in the resist material, the temperature was cooled from 373 to 298 K for 50,000 time steps. After the system attained room temperature, the mold was retracted from the PS resist. Finally, after mold removal, the system was allowed to stabilize for 80,000 time steps.
The OVITO visualization program was used to analyze molecular interactions [44]. Important physical characteristics such as displacement and position of the PS molecules were used to make inferences about indentation, penetration, and recovery behavior of the resist. The results were post-processed using Microsoft Excel and MATLAB. OVITO was used to track the displacements of the polymer chains. These output metrics were exported to MATLAB, and a custom code was developed to track the molecules directly beneath the lowest point of the spherical mold and surface molecules on the top of the resist material [45].

Results
Molecular dynamics model results were analyzed for different indenter sizes and applied forces. Deformation behavior of PS substrate atoms for different phases of the NIL process revealed distinct patterns depending on the abovementioned process variables. Figure 4 shows a 0.9 nm thick cross-sectional view of PS resist at different phases of the nanoimprint lithography process. During the onset of indentation, the surface atoms of the resist experience deformations for different indenter sizes. The lowest point displays the lowest position of the indenter, consistent with the graphic in Figure 4b. Figure 2 demonstrates how surface atoms were affected by the size of the indenter for different applied forces. The final frame shows the resting position of the surface atoms at the end of the simulation. The 40 Å indenter with a 1 µN force displayed the deepest indentation and lowest recovery of the PMMA resist as compared to other indenters (Figure 5c). Similar trends were observed for the 30 Å indenter, wherein an increase in indenter force resulted in higher plastic deformation and consequently lower recovery of PMMA. The 40 and 30 Å indenters both show consistent deformation within the layers of the PS chains from the surface to subsurface regions. However, the 20 Å indenter had a localized indentation similar to a dimple-shaped asperity with minimal plastic deformation (Figure 4a). This is evident from the localized deformation of the PS chains on the surface layers as shown in Figure 4a for the 20 Å indenter. Higher forces had a deeper imprint for 30 and 40 Å indenters. An increase in the indenter size and force resulted in higher compaction of the PS chains, as seen in Figure 4. Indentation resulted in elongation of the polystyrene chains to accommodate the indenter without breakage of the chains. Thus, differences in indenter sizes and applied forces have a profound impact on the nanoimprinted profiles. chains on the surface layers as shown in Figure 4a for the 20 Å indenter. Higher forces had a deeper imprint for 30 and 40 Å indenters. An increase in the indenter size and force resulted in higher compaction of the PS chains, as seen in Figure 4. Indentation resulted in elongation of the polystyrene chains to accommodate the indenter without breakage of the chains. Thus, differences in indenter sizes and applied forces have a profound impact on the nanoimprinted profiles.

Effect of Indenter Size and Applied Force on Indentation Depth
Indentation depth (ID) is determined by the distance that a group of surface atoms moves from their original position. In this case, ID was graphed and analyzed over the span of the entire simulation. An increase in force for each size of indenter resulted in faster descent into the material. This was particularly observed for the 40 Å indenter (Figure 5c), whereas a lower force of 100 nN was adequate to cause deformation of the resist (Figure 5d). Higher indentation forces did not reveal significant differences in deformation behavior for the 20 Å indenter (Figure 5a). Figure 5d shows a "trough-shaped" profile for lower force (100 nN) that transitions to a deeper profile for higher forces (Figure 5e,f). Thus, higher forces and large diameter indenters had deeper indentations. The 40 Å indenter with 1 µ N of force caused the surface molecules to descend to the lowest point as compared to the other indenters. The surface molecules descended to the lowest position and recovered to a plateau region within the void where permanent material deformation occurred. Varying indenter size and keeping the force constant had a significant influence on the indentation profile versus

Effect on Penetration Depth
Penetration depth (PD) is defined by the average distance between the lowest point of indentation and the atoms of the top row of the resist. Penetration of each indenter size within PS can be correlated to the hoop stresses generated by external applied forces. The semi hemispherical region of the mold deforms the PS chains as they reach the Tg, thereby setting plastic deformation of the PS. Figure 6 shows that deeper penetration occurs in the substrate atoms as the force increases. For a small indenter diameter (20 Å ), less penetration occurred as compared to larger indenters. However, smaller indenters can sink in the PS due to higher stress concentration based on the lower cross-

Effect of Indenter Size and Applied Force on Indentation Depth
Indentation depth (I D ) is determined by the distance that a group of surface atoms moves from their original position. In this case, I D was graphed and analyzed over the span of the entire simulation. An increase in force for each size of indenter resulted in faster descent into the material. This was particularly observed for the 40 Å indenter (Figure 5c), whereas a lower force of 100 nN was adequate to cause deformation of the resist (Figure 5d). Higher indentation forces did not reveal significant differences in deformation behavior for the 20 Å indenter (Figure 5a). Figure 5d shows a "trough-shaped" profile for lower force (100 nN) that transitions to a deeper profile for higher forces (Figure 5e,f). Thus, higher forces and large diameter indenters had deeper indentations. The 40 Å indenter with 1 µN of force caused the surface molecules to descend to the lowest point as compared to the other indenters. The surface molecules descended to the lowest position and recovered to a plateau region within the void where permanent material deformation occurred. Varying indenter size and keeping the force constant had a significant influence on the indentation profile versus keeping the indenter size constant and varying the force.

Effect on Penetration Depth
Penetration depth (P D ) is defined by the average distance between the lowest point of indentation and the atoms of the top row of the resist. Penetration of each indenter size within PS can be correlated to the hoop stresses generated by external applied forces. The semi hemispherical region of the mold deforms the PS chains as they reach the T g , thereby setting plastic deformation of the PS. Figure 6 shows that deeper penetration occurs in the substrate atoms as the force increases. For a small indenter diameter (20 Å), less penetration occurred as compared to larger indenters. However, smaller indenters can sink in the PS due to higher stress concentration based on the lower cross-section region and subsequent higher hoop stress. Higher indenter forces (500 nN and 1 µN) showed marginal differences in penetration behavior. Thus, a threshold indentation force is required to form a permanent indentation in the polymer. Larger force and larger size indenters equaled more penetration and vice versa.

Recovery Behavior
Recovery depth (RD) is the distance a group of atoms traverses from their lowest point of indentation to their final resting position. Figure 7 shows that as larger external forces are applied to the indenter, recovery depth increases. Similarly, larger indenter sizes had higher recovery depths, irrespective of the indenter force. Large force and small size lead to much less RD than indenters with small force and large size. The highest recovery depth was observed for small force and large size indenters and vice versa. This is because lower forces enable elastic recovery of the polymeric chains of PS with limited plastic deformation. In addition, larger indenter sizes displace more PS chains causing a permanent set. Polystyrene displays a localized yielding behavior called crazing in the region of maximum deformation. Localized regions of yielding display an increase in volume (dilatation) through formation of micro-cracks, which are bridged by polymeric chains. Crazing and stress whitening are the typical deformation mechanism for polystyrene. This phenomenon occurs in regions of high dilatational stress around the periphery of the spherical indenter. The abovedescribed results provide guidelines to design NIL molds and also estimate deformation behavior of PS.

Recovery Behavior
Recovery depth (R D ) is the distance a group of atoms traverses from their lowest point of indentation to their final resting position. Figure 7 shows that as larger external forces are applied to the indenter, recovery depth increases. Similarly, larger indenter sizes had higher recovery depths, irrespective of the indenter force. Large force and small size lead to much less R D than indenters with small force and large size. The highest recovery depth was observed for small force and large size indenters and vice versa. This is because lower forces enable elastic recovery of the polymeric chains of PS with limited plastic deformation. In addition, larger indenter sizes displace more PS chains causing a permanent set. Polystyrene displays a localized yielding behavior called crazing in the region of maximum deformation. Localized regions of yielding display an increase in volume (dilatation) through formation of micro-cracks, which are bridged by polymeric chains. Crazing and stress whitening are the typical deformation mechanism for polystyrene. This phenomenon occurs in regions of high dilatational stress around the periphery of the spherical indenter. The above-described results provide guidelines to design NIL molds and also estimate deformation behavior of PS.

Recovery Percentage
Recovery percentage is an internally defined element based on previously determined parameters. As the indentation force increased, recovery percentage increased. Recovery percentage (RP) is the percentage amount that a group of molecules travels from their position at the indenter's lowest point of penetration (LP) to their final resting positions. Figure 8 shows that as the indentation force increased, recovery percentage increased. Similarly, an increase in indenter size resulted in higher recovery percentage. The lowest recovery percentage was observed for small size and large force indentation. In general, smaller indenter sizes result in the lowest RD and RP. The resist experienced similar amounts of penetration and recovery irrespective of applied force. Thus, for all the applied forces an increase in indenter size resulted in higher penetration depth, recovery depth, and recovery percentage. The 20 Å indenter penetrated the least amount into the resist and had the lowest amount of recovery. The 40 Å indenter had around 70% recovery percentage, especially for higher forces. All three indenter sizes showed similar recovery and penetration behavior regardless of the applied force. However, the shear banding deformation phenomena in polystyrene is highly temperature and rate dependent.

Recovery Percentage
Recovery percentage is an internally defined element based on previously determined parameters. As the indentation force increased, recovery percentage increased. Recovery percentage (R P ) is the percentage amount that a group of molecules travels from their position at the indenter's lowest point of penetration (L P ) to their final resting positions. Figure 8 shows that as the indentation force increased, recovery percentage increased. Similarly, an increase in indenter size resulted in higher recovery percentage. The lowest recovery percentage was observed for small size and large force indentation. In general, smaller indenter sizes result in the lowest R D and R P . The resist experienced similar amounts of penetration and recovery irrespective of applied force. Thus, for all the applied forces an increase in indenter size resulted in higher penetration depth, recovery depth, and recovery percentage. The 20 Å indenter penetrated the least amount into the resist and had the lowest amount of recovery. The 40 Å indenter had around 70% recovery percentage, especially for higher forces. All three indenter sizes showed similar recovery and penetration behavior regardless of the applied force. However, the shear banding deformation phenomena in polystyrene is highly temperature and rate dependent.

Recovery Percentage
Recovery percentage is an internally defined element based on previously determined parameters. As the indentation force increased, recovery percentage increased. Recovery percentage (RP) is the percentage amount that a group of molecules travels from their position at the indenter's lowest point of penetration (LP) to their final resting positions. Figure 8 shows that as the indentation force increased, recovery percentage increased. Similarly, an increase in indenter size resulted in higher recovery percentage. The lowest recovery percentage was observed for small size and large force indentation. In general, smaller indenter sizes result in the lowest RD and RP. The resist experienced similar amounts of penetration and recovery irrespective of applied force. Thus, for all the applied forces an increase in indenter size resulted in higher penetration depth, recovery depth, and recovery percentage. The 20 Å indenter penetrated the least amount into the resist and had the lowest amount of recovery. The 40 Å indenter had around 70% recovery percentage, especially for higher forces. All three indenter sizes showed similar recovery and penetration behavior regardless of the applied force. However, the shear banding deformation phenomena in polystyrene is highly temperature and rate dependent.

Effect of Indenter Size and Indentation Hold Time
The indenter size had a significant impact on material deformation of the PS resist. As the size of the indenter increased, surface molecules experienced permanent plastic deformation. The effect of indentation time on deformation behavior was evaluated for the 30 Å indenter with an applied force of 500 nN. Figure 9 shows that an increase in indentation hold time results in permanent plastic deformation of the substrate. The original indentation hold time for all simulations was 175 fs. Herein, the polymer chains are able to recover to higher positions due to the elastic component. However, an increase in the indentation hold time to 400 fs resulted in lower recovery of the polymeric chains due to permanent plastic deformation. These findings have implications for determining the indentation hold time, which in turn can affect imprint quality for different sizes and applied forces to the indenters. Our results are consistent with findings in literature [31], wherein polystyrene displays shear banding and yield behavior in compression mode. The spherical indenter causes shear banding, which involves localized orientation of the polystyrene. Thus, an optimal hold time exists for each mold geometry, size, and material type.

Effect of Indenter Size and Indentation Hold Time
The indenter size had a significant impact on material deformation of the PS resist. As the size of the indenter increased, surface molecules experienced permanent plastic deformation. The effect of indentation time on deformation behavior was evaluated for the 30 Å indenter with an applied force of 500 nN. Figure 9 shows that an increase in indentation hold time results in permanent plastic deformation of the substrate. The original indentation hold time for all simulations was 175 fs. Herein, the polymer chains are able to recover to higher positions due to the elastic component. However, an increase in the indentation hold time to 400 fs resulted in lower recovery of the polymeric chains due to permanent plastic deformation. These findings have implications for determining the indentation hold time, which in turn can affect imprint quality for different sizes and applied forces to the indenters. Our results are consistent with findings in literature [31], wherein polystyrene displays shear banding and yield behavior in compression mode. The spherical indenter causes shear banding, which involves localized orientation of the polystyrene. Thus, an optimal hold time exists for each mold geometry, size, and material type.

Conclusions
This research explains the deformation behavior of polystyrene (PS) in nanoimprint lithography based on variations in process parameters using molecular dynamics modeling. Viscoelastic behavior of PS was modeled using consistent valence force field. Different indenter sizes, forces, and imprinting durations were varied to evaluate their effect on nanoimprinting of PS. The molecular mechanism of nanoimprinting was characterized with respect to indentation depth, penetration depth, recovery depth, and recovery percentage of the polymer. For larger indenter sizes, the polymer material experienced a deeper indentation profile, more penetration, and high recovery rates. The 40 Å indenter with a 1 µ N of force caused surface molecules to descend to the lowest point compared to the other indenters. Thus, varying indenter force and size provides insight into deformation behavior and thereby potential application of PS. An increase in indentation hold time beyond 250 fs resulted in higher plastic deformation of the PS. The molecular model developed in this research will be expanded to include nanoimprint lithography of other polymer composites.
Author Contributions: J.O. contributed to this original research in areas of conceptualization, methodology, software, formal analysis, investigation, data curation, and visualization. S.D. contributed in the areas of

Conclusions
This research explains the deformation behavior of polystyrene (PS) in nanoimprint lithography based on variations in process parameters using molecular dynamics modeling. Viscoelastic behavior of PS was modeled using consistent valence force field. Different indenter sizes, forces, and imprinting durations were varied to evaluate their effect on nanoimprinting of PS. The molecular mechanism of nanoimprinting was characterized with respect to indentation depth, penetration depth, recovery depth, and recovery percentage of the polymer. For larger indenter sizes, the polymer material experienced a deeper indentation profile, more penetration, and high recovery rates. The 40 Å indenter with a 1 µN of force caused surface molecules to descend to the lowest point compared to the other indenters. Thus, varying indenter force and size provides insight into deformation behavior and thereby potential application of PS. An increase in indentation hold time beyond 250 fs resulted in higher plastic deformation of the PS. The molecular model developed in this research will be expanded to include nanoimprint lithography of other polymer composites.