Next Article in Journal
Investigation on the Electrochemical Corrosion Behavior of TP2 Copper and Influence of BTA in Organic Acid Environment
Previous Article in Journal
Application of Transformation Treatment to Commercial Low-Grade Electrical Steels under Different Processing Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Nanoscale Deformation Fields from Phase Field Crystal Simulations

Division of Solid Mechanics, Lund University, P.O. Box 118, SE-221 00 Lund, Sweden
*
Author to whom correspondence should be addressed.
Metals 2022, 12(10), 1630; https://doi.org/10.3390/met12101630
Submission received: 13 September 2022 / Revised: 26 September 2022 / Accepted: 26 September 2022 / Published: 28 September 2022
(This article belongs to the Section Computation and Simulation on Metals)

Abstract

:
Different methods for evaluation of displacement and strain fields based on phase field crystal (PFC) simulations are shown. Methods originally devised for molecular dynamics (MD) simulations or analysis of high-resolution microscopy images are adapted to a PFC setting, providing access to displacement and strain fields for systems of discrete atoms, such as in MD, as well as to continuous deformation fields. The latter being achieved by geometrical phase analysis. As part of the study, the application of prescribed non-affine deformations in a 3D structural PFC (XPFC) setting is demonstrated as well as an efficient numerical scheme for evaluation of PFC phase diagrams, such as, for example, those required to stabilize solid/liquid coexistence. The present study provides an expanded toolbox for using PFC simulations as a versatile numerical method in the analysis of material behavior at the atomic scale.

1. Introduction

The Phase Field Crystal (PFC) method, introduced in [1], is a non-phenomenological continuum approach to simulating microstructure evolution in materials at the atomic scale. PFC models a time-averaged atomic number density field over long, diffusive, time scales while retaining atomic-level spatial resolution. This provides access to modeling regimes that fall beyond the reach of conventional techniques, such as Density Functional Theory (DFT) and Molecular Dynamics (MD), on the atomic scale, and phase fields or level sets on the mesoscale. The incorporation of atomic length scales means that dislocation mechanisms, as well as heterogeneities such as phase or grain boundaries in polycrystal aggregates, are automatically captured. PFC was introduced as an elaboration of standard phase field formulations, but it has later been shown that PFC can be derived directly from DFT, as discussed in [2,3]. PFC modeling has been employed for investigation of a wide range of phenomena, many of which involving solid–liquid transitions, including nucleation and growth, melting and glass formation [4,5]. However, solid-state phenomena have also been targeted, such as dislocation mechanics [6,7], grain boundary properties [8,9], grain boundary and particle interaction [10], precipitation [11] and atomic-scale fracture [12]. Various applications of PFC modeling are reviewed in, for example, [13,14].
In most PFC-based studies, lattice deformation and corresponding strains are not quantified, or done so by relatively crude measures. One example is the average local strain used in [15], which is a scalar quantity that, for example, does not distinguish shear deformation from a state of tension or compression. In the present study, different methods for evaluation of deformation and strain fields from PFC simulations are discussed. The investigated methods were originally introduced for use in either MD simulations or for analysis of experimental images obtained by High-Resolution Transmission Electron Microscopy (HRTEM). In this work, the methods are adapted to a PFC setting, taking advantage of PFC providing access to both a continuous density field as well as the discrete atom/particle positions through interpolation of the density maxima. Further, when PFC is employed for analysis of structures under mechanical load, the applied deformation is usually assumed to be homogeneous, for example in [16]. As part of the present study, it is demonstrated how displacement and strain field fields under non-homogeneous, or non-affine, deformation can be evaluated. In addition, control of the PFC model requires establishing a phase diagram to be able to stabilize different phases and to maintain phase coexistence. A numerical scheme for achieving this in a flexible and efficient manner is provided in the appendix. This study concerns analysis of simulation results provided by the wide range of PFC models based on [1], as well as those based on the structural PFC (XPFC) method introduced in [17]. An alternative approach, not covered here, is the amplitude expansion PFC (APFC) method, as discussed in [18], which offers other means for strain evaluation. While PFC and XPFC mainly differ in terms of how the direct pair correlation between particles or atoms is handled, APFC constitutes a more significant step in a different direction.
This paper is structured in such way that the PFC model is discussed first in Section 2. Various methods for evaluation of displacement and strain fields are discussed in Section 3. Different example applications of the methods are demonstrated in Section 4, along with a discussion of the results. Some concluding remarks close the paper in Section 5.

2. Phase Field Crystal Model

The starting point for the PFC model formulation is the definition of a time-averaged and normalized particle density field
n ( r ) = ρ ( r ) ρ l ρ l
where r = r ( x , y , z ) is a position vector based on the spatial coordinates ( x , y , z ) in a domain Ω , ρ is the particle number density and ρ l a reference density in the liquid state. The equilibrium density of the solid phase can be approximated by the Fourier series
n ( r ) = n 0 + j A j e 2 π a i k j · r
with n 0 denoting a reference density, A j is the amplitude of mode j, assumed to be identical for all lattice planes of family j, a is the lattice parameter and k j are reciprocal lattice vectors.
The Helmholtz free energy for the system under consideration—scaled by ρ l k B T , with k B and T denoting the Boltzmann constant and the absolute temperature, respectively—can be written as a functional of the particle density n according to
F [ n ] = F i d [ n ] + F e x c [ n ] + F e x t [ n ]
where F i d is related to the ideal gas behavior of non-interacting particles, F e x c is the excess free energy due to particle interaction and F e x t is a free energy component related to external forces. Based on the definition of the number density field in Equation (1), the Helmholtz free energy functional in Equation (3) describes the difference in energy with respect to a reference state at which ρ = ρ 0 . The first two components on the right-hand side of Equation (3) are quite standard in PFC modeling and are taken, in non-dimensional form, as
F i d [ n ] = Ω d r n ( r ) 2 2 n ( r ) 3 6 + n ( r ) 4 12
and
F e x c [ n ] = 1 2 Ω d r n ( r ) C 2 n ( r )
where C 2 ( | r r | ; σ ) is a two-point correlation function, describing the interaction between two particles located at positions r and r , respectively. The correlation function also depends on an effective temperature, denoted by σ . The convolution operation C 2 n , appearing in Equation (5), is defined by
C 2 n = Ω d r C 2 ( | r r | ; σ ) n ( r ) .
Following [17], the correlation function can be approximated as a sum of Gaussian peaks in Fourier space by
C ˜ 2 ( k ) = max j exp σ 2 k j 2 2 ρ j β j exp ( k k j ) 2 2 α j 2
for all reciprocal lattice vectors k , of length k = | k | . A particular lattice is constructed from a subset of j density waves, with wave numbers k j . In the present case, wave numbers k 1 = 3 k 0 and k 2 = 2 k 0 , corresponding to the [ 111 ] and [ 200 ] planes, are used for FCC and wave numbers k 1 = 2 k 0 and k 2 = 2 k 0 , corresponding to the [ 110 ] and [ 200 ] planes, are used for BCC. With k 0 = 2 π , a unit lattice parameter is defined by a = 2 π / k 0 .
Appearing in Equation (7), ρ j is the atomic density and β j denotes the number of planes. For FCC, these parameters are ( ρ 1 , ρ 2 , β 1 , β 2 ) = ( 4 / 3 , 2 , 8 , 6 ) and for BCC they can be found as ( ρ 1 , ρ 2 , β 1 , β 2 ) = ( 2 / 2 , 1 , 12 , 6 ) . The widths of the Gaussian peaks in Equation (7) are set by α j and can be used to tune the elastic properties of the crystal structure to comply with a particular material behavior. For simplicity, α j = 1 is used here.
The last term on the right-hand side of Equation (3), related to the external free energy, is included to permit prescribed deformation of the PFC system. Following [15], the external free energy is incorporated in the free energy functional as a penalty function, taken as
F e x t [ n ] = Ω d r ξ ( r ) n ( r ) n e x t ( r + Δ r ) 2 .
The field n e x t ( r + Δ r ) is based on an approximation of the density field, as shown in Equations (A3) and (A4). The external field n e x t corresponds to the XPFC density field n ( r ) , but is subject to a prescribed change Δ r = ( Δ x , Δ y , Δ z ) in the spatial coordinates whereby r : = r + Δ r . Appearing in Equation (8), ξ ( r ) is a modulation function which is zero outside the region in which the prescribed external field n e x t is applied. Following [19], ξ ( r ) is here taken as a normalized Gaussian of width σ ξ , centered at r , according to
ξ ( r ) = 1 2 π σ ξ 2 exp ( r r ) 2 2 σ ξ 2 .
In order to accommodate both the diffusion time-scale governing the underlying XPFC dynamics, as well as the much faster relaxation of the elastic strains imposed by Equation (8), the modified PFC (MPFC) model proposed in [15,20] is adopted. The evolution of the normalized particle density field n ( r ) is thus described by the damped wave equation
γ 2 n t 2 + β n t = α 2 2 δ F [ n ] δ n
where standard diffusion-based PFC dynamics are retrieved for the parameter choices ( α , β , γ ) = ( 1 , 1 , 0 ) . PFC dynamics on the form described by Equation (10) is also adopted in [19,21]. The evolution of the XPFC system, as described by Equation (10), is implemented using the semi-implicit spectral scheme discussed in [10].

3. Evaluation of Nanoscale PFC Deformation and Strain Fields

Deformation and strain is not available directly from PFC (or XPFC) data, but must be evaluated based on the density field n. Tracing the evolution of n by Equation (10) gives access to a continuous density field and, as discussed in [8], interpolation of the density maxima can be employed to also provide the locations of the maxima, corresponding to the discrete “atom” positions. The latter makes it possible to adapt procedures introduced for evaluation of deformation and strain in MD into a PFC setting. As an initial example related to strain evaluation in PFC, it can be noted that the discrete atom positions are used in [15] to find a local average strain by considering the β = 1 N nearest neighbors of each atom α . Knowledge of the initial distance a α β between atoms α and β , in terms of the lattice parameter, the deformed distance d a α β between the same atoms can be used to find the average local strain for each atom α according to
ε ¯ α = 1 N β = 1 N d a α β a α β
This is a simple approach to obtain a strain measure that conveys a local average strain magnitude and that gives a hint towards whether atom α is under an average state of compression ( ε ¯ α < 0 ) or tension ( ε ¯ α > 0 ). However, as Equation (11) is a scalar measure no information is provided on the strain directions nor is any distinction made between states of shear and tension/compression, i.e., the individual strain components are not provided.
In order to establish a more rigorous method for PFC strain analysis, the discrete atom positions are used in conjunction with a method devised for MD in Section 3.1 to evaluate a local deformation gradient and associated strain tensors. As an alternative approach, the continuous PFC density field is used in Section 3.2 to demonstrate how procedures originally devised for analysis of high-resolution microscopy images can be used to evaluate continuous, rather than discrete, displacement and strain fields.

3.1. Discrete Deformation Gradient and Related Strain Measures

The development of a discrete or local—in the sense of relating to an individual atom—deformation gradient can be motivated based on continuum mechanics. In continuum mechanics, a continuous body Ω 0 R 3 can be defined in the material reference configuration, comprising particles positioned at X Ω 0 . At a certain time, the nonlinear deformation map ϕ ( X ) maps the reference configuration Ω 0 onto the current, spatial, configuration Ω = ϕ ( Ω 0 ) R 3 . Accordingly, at a certain time a particle at X Ω 0 can be identified in the current configuration at x = ϕ ( X ) . The continuum mechanics deformation gradient F is defined as linear mapping between tangent spaces, transforming tangent vectors in Ω 0 to tangent vectors in Ω , by
F = d x d X .
The local volume change can be found from the determinant J of the deformation gradient, providing J = det ( F ) V / V 0 where V 0 and V is the local volume in the reference and in the current configuration, respectively. Different strain measures can be defined based on the deformation gradient. The Lagrangian (Green) strain tensor is defined entirely in terms of reference coordinates and is found as
E = 1 2 F T F I
where I is the second-order identity tensor and ( · ) T denotes a transpose. Alternatively, the Eulerian (Almansi) strain tensor can be defined in the current configuration by
e = 1 2 I F F T 1 .
In the small strain limit, the Lagrangian and Eulerian strain tensors provided by Equations (13) and (14) will be approximately equal to each other and can be represented by the small strain tensor
ε = 1 2 u T + u 1 2 F T + F I
with the displacement field being defined by u = x X and where the last identity is based on the relation F = I + u . It can also be recalled that the volumetric dilatation is provided by the first strain tensor invariant
I 1 = tr ε
where tr · denotes the trace of a tensorial quantity. Further, based on the first equality in Equation (15), it is noted that the displacement gradient can be stated as u = ε + ω with the skew-symmetric rotation tensor being defined by
ω = 1 2 u u T .
Turning to the present atom-scale representation of a material body, no continuous displacement field is available from which the deformation gradient can be evaluated. Instead, a discrete version of Equation (12) is sought. Here, the method proposed in [22,23] for evaluation of a discrete atom-scale deformation gradient in molecular dynamics simulations is considered. For this purpose, it is noted that an atom with index α will be positioned at X α in the reference configuration and at x α in the current configuration. The deformation in the local neighborhood of atom α is determined by the positions of its neighbors β , located at X β and x β , respectively. The relative positions between the atom α and its neighbors β are found as
Δ X α β = X β X α Ω 0 Δ x α β = x β x α Ω .
A linear mapping between these relative positions is defined by the deformation gradient F α of atom α , appearing as
F α = Δ x α β Δ X α β .
In contrast to the continuum deformation gradient F , defined in Equation (12), the discrete analogue F α in Equation (19) is not a single unique mapping between the reference configuration Ω 0 and the current configuration Ω . Rather, evaluation of Equation (19) for atom α and each of its β neighbors results in a system of linear equations. Following [22,23], a discrete local deformation gradient F ^ α is determined to minimize the mapping error between atom α and an individual neighbor β in a least-squares sense. The squared mapping error is defined as
η α = β = 1 N Δ x α β F ^ α Δ X α β T Δ x α β F ^ α Δ X α β
where the summation is performed over all N nearest neighbors β of atom α . Minimization of η α with respect to the components of F ^ , i.e., η α / F ^ = 0 , yields a linear system of equations that can be written as
F ^ α D α = A α F ^ α = A α D α 1
where the matrices
D α = β Δ X α β Δ X α β T
and
A α = β Δ x α β Δ X α β T
were introduced. The discrete deformation gradient F ^ , available from Equation (21), represents a local averaged deformation gradient, which can be used with Equation (13) and Equation (14) to evaluate the corresponding discrete strain tensors E ^ , e ^ and ε ^ , respectively. However, application of Equation (21) requires access to at least three neighbors β in order to determine the nine components of F ^ . In addition, these neighbors need to span R 3 to render D α invertible. This implies that the β neighbors can be neither co-linear nor co-planar with atom α .
The evaluation of F ^ by Equation (21), which was proposed in in [22,23] for strain evaluation in MD simulations, is here adopted for PFC simulations. To accomplish this, discrete atom positions are identified by performing a quadratic interpolation of the density maxima in the n ( r ) -field to provide the positions of the N nearest neighbors of each atom (8 in BCC and 12 in FCC) in the current configuration, described by Δ X α β . In addition, the ideal lattice locations of the N nearest neighbors are used to provide Δ x α β in the reference configuration. The local, atomic-scale, deformation gradient F ^ is then evaluated for each atom (density maxima) using Equations (21)–(23).

3.2. Evaluation of Continuous Displacement and Strain Fields by Geometrical Phase Analysis

Geometrical phase analysis (GPA) was introduced in [24,25] as a method to evaluate displacement and strain fields from high resolution transmission electron microscopy (HRTEM) images of crystal lattices. Central to the method is evaluation of the phase components of the intensity peaks in the HRTEM image in Fourier space. The same approach is adopted here by recognizing the correspondence between the HRTEM intensity peaks and the maxima in the PFC density field. This permits evaluation of the PFC displacement and strain fields without prior interpolation of density peaks and use of discrete atom positions. Instead, the continuous PFC density field can be used directly.
GPA is based on formulating the intensity of a 2D crystal lattice image as the Fourier series
I ( r ) = g A g exp i P g + 2 π i g · r
where the periodicity of the intensities is defined by the vector g in reciprocal lattice space. Further, A g and P g denote the amplitude and phase of the periodicity, respectively. With access to the intensity field in Fourier space, single g vectors can be filtered out by masking the corresponding Bragg peaks. In the present implementation, a simple Gaussian mask is employed, taken as
M ( k ) = exp k g 2 2 σ M 2
with the width of the Gaussian being defined by the parameter σ M . In the next step, the masked intensity image is brought back to real space by an inverse Fourier transform to provide the complex image H g , which can be used to evaluate
P g ( r ) = phase H g r 2 π g · r .
This is the real-space phase difference between the masked image and the vector g , used to define the mask in Equation (25) and thereby choosing the reference lattice relative to which displacements and strains will be evaluated. If needed, the vector g can be refined by evaluating the deviation Δ g from the reference lattice g . This is achieved by fitting the gradient the gradient of P g in a region of the original image in which the strain is assumed to be uniform, using
Δ g = 1 2 π P g .
To evaluate the full 2D displacement and strain fields, the steps involving Equations (25)–(27) can be repeated for two non-parallel vectors g 1 and g 2 , corresponding to the real space vectors a 1 and a 2 . Denoting the components of these vectors by subscripts x and y, the following matrices can be defined
G = g 1 x g 2 x g 1 y g 2 y and A = a 1 x a 2 x a 1 y a 2 y
that satisfy G T A = I , as shown in [25]. From the phase differences P g 1 and P g 2 corresponding to g 1 and g 2 , the deviation from the local reciprocal reference lattice can be found using Equations (27) and (28), providing
Δ G = 1 2 π P g 1 x P g 2 x P g 1 y P g 2 y .
This makes it possible to formulate a local 2D deformation gradient F ˜ ( r ) , akin to Equation (12), between the reference lattice A and the deformed lattice A ˜ as
A ˜ r = F ˜ r A .
Inverting both sides of Equation (30) and using the relation G T = A 1 provides
G ˜ T r = G T F ˜ 1 r F ˜ r = A G ˜ T r 1 .
Finally, substituting Equation (29) for G ˜ yields
F ˜ r = I + A Δ G T r 1 .
Equation (32) permits evaluation of the local deformation gradient F ˜ based on the phase gradients. Subsequently, also the corresponding (2D) Green and Almansi strain tensors E ˜ and e ˜ can be evaluated, based on GPA, using Equations (13) and (14). It is emphasized that the strain evaluated by GPA constitutes the lattice strain relative to a reference lattice, defined by the reciprocal lattice vector g . In the limit of small deformations, Equation (32) can be approximated by
F ˜ r = I A Δ G T r
which can be recast into
F ˜ r = I + u r .
The latter formulation provides the gradient of the displacement field as
u r = A Δ G T r .
Based on Equation (35), the small strain tensor ε ˜ and rotation ω ˜ , corresponding to Equation (15) and Equation (17), are now available, based on GPA.
The phase difference maps P g 1 and P g 2 , relating the deformed lattice to a chosen reference lattice, provide the displacement field components in the direction of g 1 and g 2 , respectively, by
P g 1 ( r ) = 2 π g 1 · u ( r ) P g 2 ( r ) = 2 π g 2 · u ( r ) .
Under the requirement that g 1 and g 2 are not co-linear, the two equations in Equation (36) can be written on matrix form and inverted to obtain displacement field by
u ( r ) = 1 2 π G T 1 P g 1 ( r ) P g 2 ( r ) 1 2 π A P g 1 ( r ) P g 2 ( r )
where the identity G T A = I was used in the last step. Using Equation (37), the displacement field is also available from the PFC-based GPA.

4. Examples

To illustrate how the methods for deformation analysis, discussed in Section 3, perform in a PFC setting, some simulation examples are provided here. Two 3D simulation models are considered, as shown in Figure 1. Model 1, seen in Figure 1a, is a single crystal specimen with two notches, placed symmetrically at the top and bottom edges. The specimen is subject to tensile loading along the x-axis by an external force field defined according to Equation (8). The force field is applied by the Gaussian modulation function in Equation (9) being non-zero in the hatched regions in Figure 1a. The coordinates r in Equation (9) are thus set to define lines along the y-axis through the center of each of the two hatched regions in Figure 1a. Periodicity of the simulation domain, as required by the spectral solution method, is maintained by defining liquid regions along the left and right domain boundaries. The same liquid state is defined in the cut-out regions in the notches. This approach to prescribing deformations in a PFC setting follows the method introduced in [20]. From the same study, the parameter values ( α , β , γ ) = ( 15 , 0.9 , 1 ) are adopted, cf. Equation (10). Working with a PFC formulation on non-dimensional form, a unit lattice parameter a = 1 is used and model dimensions are in units of a. Following [20], the grid spacing is adjusted to be as close as possible to π / 8 in all coordinate directions and a non-dimensional time increment of Δ t = 1 × 10 3 is used. To model the notches and to incorporate the liquid regions along the left and right domain edges, the PFC system must be in a state of stable solid/liquid coexistence. This requires identification of the appropriate reference density n 0 and amplitudes A j in Equation (2), as well as the corresponding value of the effective temperature σ in Equation (7). Phase diagrams can be used to achieve this and an efficient numerical scheme for establishing the PFC phase diagrams, while avoiding common approximations on the amplitudes A j , is discussed in Appendix A. For the BCC structure considered in Model 1, the density field approximation in Equation (A3) is used at the effective temperature σ = 0.05 to find the amplitudes [ A 1 , A 2 ] = [ 0.11594 , 0.053652 ] as well as the reference densities in the liquid and solid states [ n l , n s ] = [ 0.17069 , 0.11550 ] . These parameter settings place the PFC system in a stable state of solid/liquid coexistence, cf. Figure A1a.
Application of the deformation by MPFC is attractive as it allows simulation of a non-affine, or heterogeneous, deformation field. This is a more relevant approach compared to alternative methods used in other PFC studies in which a homogeneous, or affine, deformation is used when considering boundary value problems. Homogeneous deformation in a PFC setting is for example used in [16] based on a constant volume assumption and involve a continuous modification of the computational grid.
Model 2, illustrated in Figure 1b, is the second simulation set-up under consideration. The simulation model comprises two crystals, defined by orientations R 1 and R 2 , which are separated by two symmetrical tilt grain boundaries with tilt angle θ . The grain boundaries are initiated as narrow liquid strips and over the course of the simulation the two crystals rapidly grow into the liquid strips to form the minimum energy grain boundary configuration. The steps involved in setting up this type of PFC simulation domain, for example in terms of ensuring periodicity, are detailed in [8]. For Model 2, no external force field is applied as deformation is caused by the presence of the grain boundaries. Consequently, the F e x t component is omitted in the free energy functional in Equation (3) and the parameters ( α , β , γ ) = ( 1 , 1 , 0 ) are used in Equation (10). The time increment and the grid spacing are both set identical to Model 1.

4.1. Discrete Deformation Fields under Non-Affine Deformation

The notches are defined by the parameters r c = 10 and l c = 10 , cf. Figure 1a, and the domain size is set to L x , L y , L z = 121 , 81 , 3 . The external force field is applied over the hatched regions in Figure 1a, of width 10 and the liquid regions along the left and right domain boundaries also have a width of 10. A tensile deformation rate of 5 × 10 2 , without units in the present non-dimensional formulation, along the x-axis is used. As shown in Figure 1a, the single crystal BCC lattice is oriented with the [ 100 ] , [ 010 ] and [ 001 ] axes parallel to the x, y and z axes, respectively. The model comprises approximately 50,000 atoms (density peaks) which are localized by performing a quadratic interpolation in the computational grid at regular intervals. The intervals are chosen such that the magnitude of the movement Δ r α of an individual atom α between two interpolations is much less than the lattice parameter. The atom positions in the previously interpolated state are stored and each atom can be uniquely traced throughout the process by finding its incremental displacement relative to the position of its nearest neighbor among the previously stored position data set. The total accumulated displacement is also stored, along each coordinate direction, to provide the full displacement history for each atom by performing the update u α : = b m u α + Δ r α . This is in line with the procedure outlined in [26], where the small strain tensor in Equation (15) is evaluated based on the atom displacement field in post-processing of PFC data. In the present study, a more general approach is taken as the discrete deformation gradient defined in Equation (19) is considered instead. For illustration purposes, however, the x-component u x of the displacement field obtained from the present PFC simulations using Model 1, cf. Figure 1a, is shown in Figure 2 and the y-component u y in Figure 3. The sequence of images in these figures shows how deformation is localized at the notch tips and the subsequent growth of a cleavage crack between the [ 100 ] planes across the vertical direction of the simulation domain.
The same sequence of simulation steps as in Figure 2 is also shown in Figure 4, but in terms of the average strain ε ¯ obtained using Equation (11). A gradually increasing concentration of a tensile strain state at the notch tips is evident.
Similar to the results in Figure 4, the x-component E ^ x x of the Green–Lagrange strain tensor E ^ , obtained by using the discrete deformation gradient F ^ from Equation (21) in Equation (13), is shown in Figure 5. However, while no additional information can be drawn from the average strain ε ¯ , the individual components of the strain tensor adds further possibilities to describe the strain state.
Figure 6 shows the strain E ^ y y in the vertical y-direction, indicating a narrow zone of compressive strain in the vertical direction at the notch tips. This is in accord with the macroscale strain field in notched specimens, observed using digital image correlation in [27]. As a further example of the additional information residing in the full tensorial strain description, compared to the average strain ε ¯ , Figure 7 shows the shear strain component E ^ x y which also traces the deformation localization at the notch tip throughout the simulated fracture process.

4.2. Continuous Deformation Fields Based on Geometrical Phase Analysis

Model 2, shown in Figure 1b, is used to illustrate the use of Geometrical Phase Analysis (GPA) on the PFC data. Two BCC crystals are defined with the [ 100 ] , [ 010 ] and [ 001 ] axes initially aligned with the x, y and z axes. The crystal orientations R 1 and R 2 are then defined by rotating the crystals through an angle of θ / 2 in opposite directions around the z-axis. In this study, a tilt angle of θ = 12 . 68 is used. This will result in low-angle tilt grain boundaries, cf. Figure 1b, comprising clearly distinguishable arrays of grain boundary dislocations. The domain size is taken as L x , L y , L z = 61.3 , 66.5 , 9.00 , with the fractional numbers being due to the domain periodicity requirement.
The different steps in the GPA analysis are illustrated in Figure 8. The starting point is the PFC density field, shown in Figure 8a, where it is noted that the density maxima correspond to the intensity peaks seen in HRTEM images of crystal lattices. This observation conveniently makes GPA directly applicable to the analysis of PFC density fields. A Fourier transform (FFT) of the density field provides the intensity map in Figure 8b and the two peaks used for the subsequent analysis are indicated by red circles, showing the reciprocal lattice vectors g 1 = [ 100 ] and g 2 = [ 010 ] . A mask is placed, in turn, at g 1 and g 2 using Equation (25). The mask is multiplied by the intensity data in Fourier space (a convolution in real space) and the result is then brought back to real space by a reverse FFT. Figure 8c shows the real part of the resulting complex H g 1 map after masking by g 1 . The same procedure is repeated using g 2 instead, providing H g 2 . Using Equation (26), the phase difference maps P g 1 and P g 2 related to g 1 and g 2 , respectively, are evaluated. As an example, Figure 8d shows the phase difference map P g 1 , normalized to the range [ π , π ] . A phase difference of 2 π , indicated by an abrupt discontinuous shift from black to white in the figure, corresponds to the presence of additional lattice planes. Using Equation (37), the 2D displacement field can be evaluated based on the phase difference maps P g 1 and P g 2 . Figure 8e,f show the displacement components along the x and y axis, respectively. It is emphasized, however, that in contrast to the discrete atom displacements shown in Figure 2 and Figure 3, the GPA-based displacements in Figure 8e–f are measured relative to the chosen reference lattice, defined by g 1 and g 2 .
Based on Equation (32), a deformation gradient F ˜ can be obtained by GPA and corresponding strain tensors can be evaluated by using Equations (13)–(15). It is emphasized that just as for the displacements also the strains obtained from GPA will be relative to the chosen reference lattice. Figure 9 shows some strain and rotation components obtained from GPA applied to the PFC data provided by Model 2. Characteristic alternating tension/compression regions at the individual dislocations are clearly distinguishable in the E ˜ x x field shown in Figure 9a. In addition, Figure 9b shows the concurrent E ˜ y y strain component, considerably smaller in magnitude than E ˜ x x , and Figure 9c shows the shear strain component E ˜ x y . The PFC-based strain fields evaluated by GPA are in qualitative agreement with the analytical elastic strain fields expected for edge dislocations [28]. The rotation component ω ˜ x y , cf. Equation (17), is shown in Figure 9d and agrees with the crystal rotations described by ± θ / 2 . The rotation magnitude in Figure 9d is evaluated as tan ( ω ˜ x y ) and is shown in degrees.
The GPA method, introduced in [25] and described in Section 3.2, provides the displacement and strain fields over 2D surfaces. Unfortunately, the extension to 3D is not straightforward as the selection of Bragg peaks for masking of the intensity field becomes quite involved in the 3D analogue to Figure 8b. A multitude of density peaks are present due to the symmetries of the lattices. In addition, the GPA results are highly sensitive to the selection of the reciprocal lattice vectors g and manual refinement using Equation (27) is usually required for the individual vectors. Moreover, the size σ M of the mask in Equation (25) can in general not be held constant, but must be manually adjusted for each g -vector. Taken together, these factors defy an autonomous algorithmic implementation of GPA for the general 3D case. As a compromise, however, different planar 2D sections through the 3D PFC data can be analyzed by GPA to provide the strain components on these planes. Figure 10 gives some examples of how strain components can be evaluated over orthogonal planes through a point in the 3D PFC data. This procedure is, due to the manual steps involved, mainly useful to access the full strain tensor in singular points or small regions of particular interest.

5. Conclusions

The present study demonstrates how numerical methods for the evaluation of deformation fields, originally intended for molecular dynamics simulations or for analysis of high-resolution microscopy images, can be adapted to a PFC setting. Two main strategies are identified and demonstrated: either evaluation of discrete deformation fields based on the interpolated positions of individual PFC density maxima (identified as atoms) or evaluation of continuous deformation fields, based on a spectral analysis of the full density field using Geometrical Phase Analysis (GPA). It is shown that both the discrete and the continuous field approach can be used to post-process PFC data to obtain access to the full displacement and strain fields. Working with discrete points, 3D data sets can be conveniently handled and this method is also attractive for large-scale simulations where it becomes unfeasible to store the full PFC density field in comparison to just storing data at the discrete density maxima. Alternatively, working with the continuous data using GPA, more highly resolved deformation fields can be evaluated in 2D. GPA-based strain analysis becomes restricted in 3D as it is not straightforward to perform the required phase analysis in three dimensions. A compromise is demonstrated in the present work, however, evaluating the full deformation field by its individual components on orthogonal 2D sections through the 3D PFC density field. Extension of 3D PFC strain analysis based on spectral methods will be a target for future work. In addition to the different approaches to evaluation of PFC deformation fields, the present study also demonstrates application of prescribing non-affine deformations in a 3D structural PFC (XPFC) setting (to the author’s knowledge only 2D implementations exist in the literature) and an efficient and general numerical scheme for evaluation of PFC phase diagrams is provided in Appendix A. Taken together, the findings of this study expand the toolbox available for using PFC as a versatile method for numerical simulations of material behavior at the atomic scale.

Author Contributions

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

Funding

This research was funded by the Swedish Research Council under grant number 2019-03945 and by the strategic innovation program LIGHTer provided by the Swedish Innovation Agency.

Data Availability Statement

The data presented in this study are contained within the article.

Acknowledgments

This research was supported by a hardware grant from NVIDIA and utilized NVIDIA Quadro RTX8000 GPUs.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Identification of PFC Equilibrium States

In order to have a system with coexistence between solid and liquid phases, phase diagrams are established to identify the conditions for such phase coexistence. The equilibrium density of the solid phase can be approximated by Equation (1) to yield
n s ( r ) = n 0 + j A j e 2 π a i k j · r .
The equilibrium density in the liquid state is found by setting A j = 0 in Equation (A1) which provides
n l ( r ) = n 0 .
Note that subscripts s and l are introduced on the left-hand sides of Equation (A1) and Equation (A2) to distinguish components in the solid and liquid phase, respectively. It is also noted that the reference density n 0 in general differ between the liquid and solid phases. In the present case, explicit formats of Equation (A1) for BCC and FCC lattices are provided by
n s , b c c ( r ) = n 0 + 4 A 1 cos ( k 0 x ) cos ( k 0 y ) + cos ( k 0 x ) cos ( k 0 z ) + cos ( k 0 y ) cos ( k 0 z ) + 2 A 2 cos ( 2 k 0 x ) + cos ( 2 k 0 y ) + cos ( 2 k 0 z )
and
n s , f c c ( r ) = n 0 + 8 A 1 cos ( k 0 x ) cos ( k 0 y ) cos ( k 0 z ) + 2 A 2 cos ( 2 k 0 x ) + cos ( 2 k 0 y ) + cos ( 2 k 0 z )
where k 0 = 2 π / a which means spatial dimensions are taken in units of the lattice parameter a. It is emphasized that the method outlined here, with BCC and FCC as example structures, is identical and applicable for also for other crystal structures, for example involving SC and HCP lattices.
As a next step, either of the approximations in Equation (A1) or Equation (A2) are inserted into the free energy functional F [ n ] = F i d [ n ] + F e x c [ n ] in Equations (5) and (6), neglecting the external free energy F e x t . The resulting free energy functional is integrated over a 3D unit cell, placing the origin at one of the cell’s lattice points, and subsequently divided by the unit cell volume. This procedure yields one free energy density f ¯ s ( A j , n 0 , σ ) for the solid state and another free energy density f l ( n 0 , σ ) for the liquid state.
As noted in [17,29], the integration over a unit cell provides the excess free energy as
F e x c = 1 2 β 1 A 1 2 C ˜ 2 ( k 1 ) 1 2 β 2 A 2 2 C ˜ 2 ( k 2 )
where it can be assumed that the correlation kernels satisfy C ˜ 2 ( k 1 ) = C ˜ 2 ( k 2 ) = 1 , as noted in [17,29]. To eliminate the dependence on the amplitudes A j , f ¯ s is minimized with respect to A j by evaluating
f ¯ s A j = 0 A j
which provides A j = A j ( n 0 , σ ) and hence f s = f ¯ s ( A j ( n 0 , σ ) , n 0 , σ ) , allowing construction of a phase diagram in the n 0 σ plane. Under isothermal conditions, the phases are in equilibrium when they have uniform chemical potential μ , requiring
f s n 0 | n 0 = n 0 , s = f l n 0 | n 0 = n 0 , l = μ
with n 0 , s and n 0 , l being the average density in the solid and liquid phase, respectively. Coexistence also requires a uniform pressure in the phases, as the grand potentials of the solid and liquid phases should be in equilibrium, whereby the standard common tangent construction provides
f l ( n 0 , l , σ ) f s ( n 0 , s , σ ) n 0 , l n 0 , s = μ .
In previous studies, such as [17], phase diagram construction is simplified by assuming that the mode amplitudes satisfy either A 1 = A 2 or that A 2 is a function of A 1 . To avoid any a priori assumptions on the amplitudes A j , Equations (A6)–(A8) are combined into a system of nonlinear residual equations to provide
R s = f l ( n 0 , l , σ ) f s ( n 0 , s , σ ) n 0 , l n 0 , s f s n 0 | n 0 = n 0 , s = 0 R l = f l ( n 0 , l , σ ) f s ( n 0 , s , σ ) n 0 , l n 0 , s f l n 0 | n 0 = n 0 , l = 0 R A 1 = f ¯ s A 1 = 0 R A 2 = f ¯ s A 2 = 0 .
The residuals defined in Equation (A9) are collected in a vector 𝓡 = R s , R l , R A 1 , R A 2 and the iteration variables are collected in a vector Y = n s , n l , A 1 , A 2 . For a given value of the effective temperature σ , the iteration variables can now be identified by adopting the Newton-Raphson scheme
Y : = Y J 1 𝓡
where the components of the jacobian J = 𝓡 / Y are found by differentiating the residuals comprising 𝓡 , as defined in Equation (A9). The scheme in Equation (A10) can be conveniently put inside a loop over σ to produce phase diagrams for a chosen range of σ . Figure A1 shows the results for BCC- and FCC-structures, respectively. It can be noted from the figures that the relations between the amplitudes A 1 and A 2 are neither constant, as proposed in [17], nor are the quotients A 2 / A 1 confined to small variations over the inspected range of σ -values, as seen in Figure A1b.
Figure A1. (a) Phase diagram for liquid and BCC phases. (b) Density field amplitudes for BCC, Equation (A3). (c) Phase diagram for liquid and FCC phases. (d) Density field amplitudes for FCC, Equation (A4). Note that the horizontal axes in figures a and c show both n l and n s as the lines separating liquid/coexistence and solid/coexistence are plotted in the same figure.
Figure A1. (a) Phase diagram for liquid and BCC phases. (b) Density field amplitudes for BCC, Equation (A3). (c) Phase diagram for liquid and FCC phases. (d) Density field amplitudes for FCC, Equation (A4). Note that the horizontal axes in figures a and c show both n l and n s as the lines separating liquid/coexistence and solid/coexistence are plotted in the same figure.
Metals 12 01630 g0a1

References

  1. Elder, K.; Katakowski, M.; Haataja, M.; Grant, M. Modeling Elasticity in Crystal Growth. Phys. Rev. Lett. 2002, 88, 245701. [Google Scholar] [CrossRef] [PubMed]
  2. van Teeffelen, S.; Backofen, R.; Voigt, A.; Löwen, H. Derivation of the phase-field-crystal model for colloidal solidification. Phys. Rev. E 2009, 79, 051404. [Google Scholar] [CrossRef] [PubMed]
  3. Huang, Z.F.; Elder, K.; Provatas, N. Phase-field-crystal dynamics for binary systems: Derivation from dynamical density functional theory, amplitude equation formalism, and applications to alloy heterostructures. Phys. Rev. E 2010, 82, 021605. [Google Scholar] [CrossRef] [PubMed]
  4. Asadi, E.; Zaeem, M.A. Quantitative phase-field crystal modeling of solid-liquid interfaces for FCC metals. Comput. Mater. Sci. 2017, 127, 236–243. [Google Scholar] [CrossRef]
  5. Jreidini, P.; Provatas, N. Phase-field crystal model for materials with anomalous expansion during solidification and its application to the cavitation of supercooled water droplets. Phys. Rev. Mater. 2022, 6, 053404. [Google Scholar] [CrossRef]
  6. Gao, Y.; Huang, L.; Deng, Q.; Zhou, W.; Luo, Z.; Lin, K. Phase field crystal simulation of dislocation configuration evolution in dynamic recovery in two dimensions. Acta Mater. 2016, 117, 238–251. [Google Scholar] [CrossRef]
  7. Li, H.; Wang, X.; Zhang, H.; Tian, X.; Hou, H.; Zhao, Y. Arrangement and Decomposition of Grain Boundary Dislocations: Two-Mode Phase-Field Crystal Simulation. Front. Mater. 2022, 9, 875519. [Google Scholar] [CrossRef]
  8. Blixt, K.H.; Hallberg, H. Evaluation of grain boundary energy, structure and stiffness from phase field crystal simulations. Model. Simul. Mater. Sci. Eng. 2022, 30, 014002. [Google Scholar] [CrossRef]
  9. Blixt, K.H.; Hallberg, H. Grain boundary stiffness based on phase field crystal simulations. Mater. Lett. 2022, 318, 132178. [Google Scholar] [CrossRef]
  10. Blixt, K.H.; Hallberg, H. Grain boundary and particle interaction: Enveloping and pass-through mechanisms studied by 3D phase field crystal simulations. Mater. Des. 2022, 220, 110845. [Google Scholar] [CrossRef]
  11. Shuai, X.; Mao, H.; Tang, S.; Kong, Y.; Du, Y. Growth modes of grain boundary precipitate in aluminum alloys under different lattice misfits. J. Mater. Sci. 2022, 57, 2744–2757. [Google Scholar] [CrossRef]
  12. Gao, Y.; Deng, Q.; Huang, L.; Ye, L.; Wen, Z.; Luo, Z. Atomistic modeling for mechanism of crack cleavage extension on nano-scale. Comput. Mater. Sci. 2017, 130, 64–75. [Google Scholar] [CrossRef]
  13. Emmerich, H.; Löwen, H.; Wittkowski, R.; Gruhn, T.; Tóth, G.; Tegze, G.; Gránásy, L. Phase-field-crystal models for condensed matter dynamics on atomic length and diffusive time scales: An overview. Adv. Phys. 2012, 61, 665–743. [Google Scholar] [CrossRef]
  14. Asadi, E.; Zaeem, M.A. A Review of Quantitative Phase-Field Crystal Modeling of Solid-Liquid Structures. JOM 2015, 67, 186–201. [Google Scholar] [CrossRef]
  15. Stefanovic, P.; Haataja, M.; Provatas, N. Phase field crystal study of deformation and plasticity in nanocrystalline materials. Phys. Rev. E 2009, 80, 046107. [Google Scholar] [CrossRef] [PubMed]
  16. Hirouchi, T.; Takaki, T.; Tomita, Y. Development of numerical scheme for phase field crystal deformation simulation. Comput. Mater. Sci. 2009, 44, 1192–1197. [Google Scholar] [CrossRef]
  17. Greenwood, M.; Ofori-Opoku, N.; Rottler, J.; Provatas, N. Modeling structural transformations in binary alloys with phase field crystals. Phys. Rev. B 2011, 84, 064104. [Google Scholar] [CrossRef]
  18. Salvalaglio, M.; Elder, K. Coarse-grained modeling of crystals by the amplitude expansion of the phase-field crystal model: An overview. Model. Simul. Mater. Sci. Eng. 2022, 30, 053001. [Google Scholar] [CrossRef]
  19. Trautt, Z.; Adland, A.; Karma, A.; Mishin, Y. Coupled motion of asymmetrical tilt grain boundaries: Molecular dynamics and phase field crystal simulations. Acta Mater. 2012, 60, 6528–6546. [Google Scholar] [CrossRef]
  20. Stefanovic, P.; Haataja, M.; Provatas, N. Phase-Field Crystals with Elastic Interactions. Phys. Rev. Lett. 2006, 96, 225504. [Google Scholar] [CrossRef] [Green Version]
  21. Zhou, W.; Wang, J.; Wang, Z.; Huang, Z.F. Mechanical relaxation and fracture of phase field crystals. Phys. Rev. E 2019, 99, 013302. [Google Scholar] [CrossRef] [PubMed]
  22. Gullett, P.; Horstemeyer, M.; Baskes, M.; Fang, H. A deformation gradient tensor and strain tensors for atomistic simulations. Model. Simul. Mater. Sci. Eng. 2008, 16, 015001. [Google Scholar] [CrossRef]
  23. Zimmerman, J.; Bammann, D.; Gao, H. Deformation gradients for continuum mechanical analysis of atomistic simulations. Int. J. Solids Struct. 2009, 46, 238–253. [Google Scholar] [CrossRef]
  24. Hÿtch, M. Geometric Phase analysis of High Resolution electron Microscope Images. Scanning Microsc. 1997, 11, 53–66. [Google Scholar]
  25. Hÿtch, M.; Snoeck, E.; Kilaas, R. Quantitative measurement of displacement and strain fields from HREM micrographs. Ultramicroscopy 1998, 74, 131–146. [Google Scholar] [CrossRef]
  26. Bernal, F.; Backofen, R.; Voigt, A. Elastic interactions in phase-field crystal models: Numerics and postprocessing. Int. J. Mater. Res. 2010, 101, 467–472. [Google Scholar] [CrossRef]
  27. Rickhey, F.; Hong, S. Stress Triaxiality in Anisotropic Metal Sheets - Definition and Experimental Acquisition for Numerical Damage Prediction. Materials 2022, 15, 3738. [Google Scholar] [CrossRef] [PubMed]
  28. Weertman, J.; Weertman, J. Elementary Dislocation Theory; Oxford University Press: Oxford, UK, 1992. [Google Scholar]
  29. Greenwood, M.; Rottler, J.; Provatas, N. Phase-field-crystal methodology for modeling of structural transformations. Phys. Rev. E 2011, 83, 031601. [Google Scholar] [CrossRef]
Figure 1. Schematic illustration of the two 3D PFC simulation models. (a) Model 1: Notched sample subject to prescribed deformation. The deformation is applied in the hatched regions by an external force field described using Equations (8) and (9). Liquid regions are defined along the left and right domain edges to permit use of periodic boundary conditions, as required by the spectral solution scheme. (b) Model 2: Domain with two crystals having orientations R 1 and R 2 , separated by symmetrical tilt grain boundaries with tilt angle θ .
Figure 1. Schematic illustration of the two 3D PFC simulation models. (a) Model 1: Notched sample subject to prescribed deformation. The deformation is applied in the hatched regions by an external force field described using Equations (8) and (9). Liquid regions are defined along the left and right domain edges to permit use of periodic boundary conditions, as required by the spectral solution scheme. (b) Model 2: Domain with two crystals having orientations R 1 and R 2 , separated by symmetrical tilt grain boundaries with tilt angle θ .
Metals 12 01630 g001
Figure 2. Discrete atom displacement field. The images show the displacement component u x along the x-axis at (a) the initial state, (b) an intermediate state and (c) the final state.
Figure 2. Discrete atom displacement field. The images show the displacement component u x along the x-axis at (a) the initial state, (b) an intermediate state and (c) the final state.
Metals 12 01630 g002
Figure 3. Discrete atom displacement field. The images show the displacement component u y along the y-axis at (a) the initial state, (b) an intermediate state and (c) the final state.
Figure 3. Discrete atom displacement field. The images show the displacement component u y along the y-axis at (a) the initial state, (b) an intermediate state and (c) the final state.
Metals 12 01630 g003
Figure 4. Average discrete strain field ε ¯ , evaluated by Equation (11), at (a) the initial state, (b) an intermediate state and (c) the final state.
Figure 4. Average discrete strain field ε ¯ , evaluated by Equation (11), at (a) the initial state, (b) an intermediate state and (c) the final state.
Metals 12 01630 g004
Figure 5. Discrete Green strain component E ^ x x , evaluated by Equation (13) based on the discrete deformation gradient in Equation (21), at (a) the initial state, (b) an intermediate state and (c) the final state.
Figure 5. Discrete Green strain component E ^ x x , evaluated by Equation (13) based on the discrete deformation gradient in Equation (21), at (a) the initial state, (b) an intermediate state and (c) the final state.
Metals 12 01630 g005
Figure 6. Discrete Green strain component E ^ y y , evaluated by Equation (13) based on the discrete deformation gradient in Equation (21), at (a) the initial state, (b) an intermediate state and (c) the final state.
Figure 6. Discrete Green strain component E ^ y y , evaluated by Equation (13) based on the discrete deformation gradient in Equation (21), at (a) the initial state, (b) an intermediate state and (c) the final state.
Metals 12 01630 g006
Figure 7. Discrete Green strain component E ^ x y , evaluated by Equation (13) based on the discrete deformation gradient in Equation (21), at (a) the initial state, (b) an intermediate state and (c) the final state.
Figure 7. Discrete Green strain component E ^ x y , evaluated by Equation (13) based on the discrete deformation gradient in Equation (21), at (a) the initial state, (b) an intermediate state and (c) the final state.
Metals 12 01630 g007
Figure 8. Steps in strain field evaluation using GPA based on PFC density field data: (a) Continuous PFC density field; (b) Intensity of the Fourier transform of the density field in a. Darker spots indicate high intensity peaks; (c) Result after masking the intensity at g 1 in Fourier space and then transforming it back to real space. The real part of the resulting complex map is shown; (d) Phase difference map evaluated according to Equation (26); (e,f) Displacement field components in the x- and y-directions, respectively, evaluated by Equation (37).
Figure 8. Steps in strain field evaluation using GPA based on PFC density field data: (a) Continuous PFC density field; (b) Intensity of the Fourier transform of the density field in a. Darker spots indicate high intensity peaks; (c) Result after masking the intensity at g 1 in Fourier space and then transforming it back to real space. The real part of the resulting complex map is shown; (d) Phase difference map evaluated according to Equation (26); (e,f) Displacement field components in the x- and y-directions, respectively, evaluated by Equation (37).
Metals 12 01630 g008
Figure 9. Examples of strain and rotation components based on the discrete deformation gradient F ˜ in Equation (32), evaluated by GPA. Strains and rotations are found by Equation (13) and Equation (17), respectively. (a) Strain component E ˜ x x ; (b) Strain component E ˜ y y ; (c) Strain component E ˜ x y ; (d) Rotation component ω ˜ x y (rotation around the z-axis, measured in degrees [ ] ). Note the different scale of the color legend in figure (b), compared to those in figures (a,c).
Figure 9. Examples of strain and rotation components based on the discrete deformation gradient F ˜ in Equation (32), evaluated by GPA. Strains and rotations are found by Equation (13) and Equation (17), respectively. (a) Strain component E ˜ x x ; (b) Strain component E ˜ y y ; (c) Strain component E ˜ x y ; (d) Rotation component ω ˜ x y (rotation around the z-axis, measured in degrees [ ] ). Note the different scale of the color legend in figure (b), compared to those in figures (a,c).
Metals 12 01630 g009
Figure 10. (a) PFC density field n ( r ) , shown on three orthogonal planes through the 3D simulation domain. Figures (bd) show examples of strain components obtained by GPA applied to the PFC results. The strain components are shown on the same planes as indicated in figure (a); (b) Strain component E ˜ x x ; (c) Strain component E ˜ y y ; (d) Strain component E ˜ z z . Note the different scale of the color legend in figure (c), compared to those in figures (b,d).
Figure 10. (a) PFC density field n ( r ) , shown on three orthogonal planes through the 3D simulation domain. Figures (bd) show examples of strain components obtained by GPA applied to the PFC results. The strain components are shown on the same planes as indicated in figure (a); (b) Strain component E ˜ x x ; (c) Strain component E ˜ y y ; (d) Strain component E ˜ z z . Note the different scale of the color legend in figure (c), compared to those in figures (b,d).
Metals 12 01630 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hallberg, H.; Hult Blixt, K. Evaluation of Nanoscale Deformation Fields from Phase Field Crystal Simulations. Metals 2022, 12, 1630. https://doi.org/10.3390/met12101630

AMA Style

Hallberg H, Hult Blixt K. Evaluation of Nanoscale Deformation Fields from Phase Field Crystal Simulations. Metals. 2022; 12(10):1630. https://doi.org/10.3390/met12101630

Chicago/Turabian Style

Hallberg, Håkan, and Kevin Hult Blixt. 2022. "Evaluation of Nanoscale Deformation Fields from Phase Field Crystal Simulations" Metals 12, no. 10: 1630. https://doi.org/10.3390/met12101630

APA Style

Hallberg, H., & Hult Blixt, K. (2022). Evaluation of Nanoscale Deformation Fields from Phase Field Crystal Simulations. Metals, 12(10), 1630. https://doi.org/10.3390/met12101630

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