Abstract
Focussed Ion Beam (FIB) milling is a mainstay of nano-scale machining. By manipulating a tightly focussed beam of energetic ions, often gallium (Ga+), FIB can sculpt nanostructures via localised sputtering. This ability to cut solid matter on the nano-scale revolutionised sample preparation across the life, earth and materials sciences. Despite its widespread usage, detailed understanding of the FIB-induced structural damage, intrinsic to the technique, remains elusive. Here we examine the defects caused by FIB in initially pristine objects. Using Bragg Coherent X-ray Diffraction Imaging (BCDI), we are able to spatially-resolve the full lattice strain tensor in FIB-milled gold nano-crystals. We find that every use of FIB causes large lattice distortions. Even very low ion doses, typical of FIB imaging and previously thought negligible, have a dramatic effect. Our results are consistent with a damage microstructure dominated by vacancies, highlighting the importance of free-surfaces in determining which defects are retained. At larger ion fluences, used during FIB-milling, we observe an extended dislocation network that causes stresses far beyond the bulk tensile strength of gold. These observations provide new fundamental insight into the nature of the damage created and the defects that lead to a surprisingly inhomogeneous morphology.
Similar content being viewed by others
Introduction
The ability of FIB to shape materials at the nano-scale has made it central to microchip prototyping1, 3D material analysis2,3, targeted electron microscopy sample extraction4,5 and the nanotechnology behind size-dependent material properties6,7, to name but a few examples. It is tempting to assume FIB milling to be atomically perfect, only removing surface atoms. This is unfortunately not the case: sufficiently energetic incident ions will displace target atoms from their equilibrium lattice positions, causing collision cascades and structural damage8. Predicting this ion-implantation damage and its effect on material properties is not straightforward9,10. The effects of ion-implantation generally remain poorly understood. Yet they have important consequences. For example, FIB-milled nano-structures have been used extensively to investigate the size-dependence of material properties, leading to the âsmaller is strongerâ paradigm6,7. However, several studies suggest that FIB-induced defects may themselves be a major contributor to the observed scale-dependence of material strength11,12.
The damage produced by FIB-milling ranges from amorphization13 to the generation of lattice defects12 and the formation of intermetallic phases14. To examine its effect on material properties, detailed measurements of the lattice strains that govern defect interactions are essential. Previously, FIB-induced strains were inferred by considering the deflection of FIB-milled cantilevers15. However, such coarse measurements cannot capture the rich detail of heterogeneous defect distributions that determine material behaviour.
To investigate in detail the complex changes brought about by FIB milling we concentrate on the non-destructive, three-dimensional nano-scale measurement of FIB-milling-induced lattice strains in initially pristine objects. This is distinct from numerous previous studies where FIB milling was used to relieve pre-existing residual strains16,17,18. Rather our experiments quantify the strains directly caused by the defects introduced by FIB-milling. Gold was chosen as a model system since near-perfect nano-crystals can be reliably grown19. Our experiments use non-destructive Bragg Coherent X-ray Diffraction Imaging (BCDI)20, where a 3D coherent X-ray diffraction pattern (CXDP), i.e. an oversampled 3D reciprocal space map, is collected from a coherently illuminated single crystal (Fig. 1(A))21. The CXDP corresponds to the intensity of the 3D Fourier transform of the electron density in the crystal. By recovering the phase of the CXDP, the real space, complex-valued, electron density can be reconstructed20. Its amplitude provides information about electron density, Ï(r), i.e. the shape of the crystal. Its phase, Ï(r), is linked to displacements, u(r), of atoms from their ideal lattice positions by Ï(r)â=âq.u(r), where q is the scattering vector of the CXDP. By combining at least three CXDPs with linearly independent q vectors, u(r) can be recovered22, and the full lattice strain tensor, ε(r), determined by differentiating u(r). Thus BCDI allows the 3D nano-scale measurement of both crystal morphology and the full lattice strain tensor.
Results
Lattice Distortions due to FIB Imaging
First we consider the effect of FIB-imaging. Before X-ray measurements nano-crystal A was exposed to a Ga+ dose just sufficient to image the sample (30âkeV, 50âpA, 4.2âÃâ104 ions/μm2). The crystal morphology (Fig. 1(B)), reconstructed with a spatial resolution of ~38ânm based on the (1-11) CXDP, is in excellent agreement with scanning electron micrographs (Fig. 1(A) and SI Fig. S1). The 3D field of lattice displacements in the [1-11] crystal direction (Fig. 1(B,C)) shows large displacements near the top, implanted, surface of the crystal. By comparison, the lattice displacements measured in an unimplanted, as-grown crystal, reconstructed with ~12ânm spatial resolution, are small (Fig. 1(D,E)). The unimplanted crystal only shows slight increases in lattice displacement at its vertices due to surface energy effects23. Displacements in the vicinity of the crystal-substrate interface are small, indicating negligible substrate-induced strains in the crystals in contrast to other, previously studied material systems21. This demonstrates that the large displacements in crystal A are caused by the Ga+ bombardment.
To further explore these FIB-induced lattice distortions, CXDPs from five reflections were used to reconstruct the full 3D-resolved lattice strain tensor, ε(r), inside crystal A (Fig. 2). The six independent components of ε(r) are shown on virtual xy and yz sections through crystal A (Fig. 2). εyy(r) is large and negative within ~30ânm of the implanted top surface, indicating a lattice contraction due to Ga+ implantation. The εxy(r) (Fig. 2(C)) and εyz(r) (Fig. 2(D)) shear components show more subtle strain features.
These strains can be understood by direct comparison with numerical calculations. Using the measured 3D morphology, an anisotropically elastic24 finite element (FE) model of crystal A was constructed (Fig. 2(E), SI Fig. S2). Simulations using the Stopping Range of Ions in Matter (SRIM) code25 predict a ~20ânm thick damage layer. The calculated implantation profiles are shown in SI Fig. S3. Accordingly a constant volumetric Eigenstrain, εv, was imposed within a 20ânm thick surface layer in the FE model. εvâ=ââ3.15âÃâ10â3 provides a good match to the experimentally measured lattice displacements. There is striking agreement between calculated strains (Fig. 2(F,G)) and measured strains (Fig. 2(C,D)). Not only are the εyy(r), εxy(r) and εyz(r) components well matched, but finer features in the other strain components also agree. This highlights that even very low Ga+ fluences lead to substantial lattice distortions, and demonstrates the unique capability of BCDI for detailed, 3D-resolved nano-scale strain analysis. The magnitude of the FIB-induced eigenstrain is similar to previous observations in silicon inferred from the distortion of FIB-milled cantilevers15.
The volumetric lattice strain due to defects is given by , where and respectively are number density and relaxation volume of defect type A10. Using density functional theory (DFT) calculations we found that for a gold monovacancy Ωr(V)â=ââ0.38, whilst all self-interstitial atom (SIA) configurations (100 dumbbell, octahedral site, 110 crowdion and 110 dumbbell, SI Fig. 4) have Ωr(SIA)â=â2.0. Hence collision damage in the bulk, involving equal numbers of SIAs and vacancies, will cause a lattice swelling. The observed lattice contraction is surprising, particularly since the relaxation volume of substitution gallium in gold is small and positive (SI Table S1). Thus the lattice contraction we measure must indicate an excess of vacancies. This can be explained by considering the proximity of the crystal surface: vacancies and vacancy clusters with high migration energy (Emââ¥â0.71âeV)26,27,28 are retained, while SIAs, which are mobile even at temperatures of a few K28, escape to the free surface. εv and Ωr(V) allow a lower bound estimate of ~230 retained vacancies/(Ga+), while SRIM calculations provide an upper bound of ~400 vacancies/(Ga+). Based on the implanted Ga+ dose and εv, these bounds correspond to defect concentrations of ~4.9âÃâ1026âvacancies/m3 (lower bound) and ~8.5âÃâ1026âvacancies/m3â+â6.79âÃâ1025 SIAs/m3 (upper bound) respectively. Thus our measurements allow quantitative insight into the nature of the damage formed, even at very low ion fluences.
Higher Dose FIB Milling
At higher Ga doses a distinctly different behaviour is observed. Nano-crystals B and C were exposed to fluences of 1.3âÃâ107 ions/μm2 and 1.5âÃâ108 ions/μm2 respectively, causing the removal of ~3ânm and ~40ânm thick surface layers by sputtering, as predicted by SRIM. Lattice displacements and strains in both crystals were reconstructed using six crystal reflections (SI Fig. S5 and Fig. 3 respectively). The spatial resolutions of these reconstructions are ~45ânm and ~47ânm respectively. Even for these highly damaged crystals agreement of the reconstructed morphology and SEM micrographs is excellent (SI Fig. S1).
The lattice displacement magnitude in crystal C (Fig. 3(A)) shows abrupt variations, in contrast to the gradual changes in crystal A (Fig. 2(A)). The εyy(r) strain (Fig. 3(C)) is no longer uniform and negative in the implanted layer, but contains compressive and tensile regions. Similar variations are present in the other strain components.
The nature of the underlying crystallographic defects can be explored by considering the amplitudes and phases of the complex electron density reconstructed from different crystal reflections. Fig. 3(D) shows the amplitudes and phases of the (200), (020) and (002) reflections for the area marked by a black rectangle in Fig. 3(C). The (020) and (002) reflection phases both show a phase jump of ~4.2 radians, while the (200) reflection phase shows no discontinuity. The structure of the phase jumps suggests dislocations as the underlying defects29,30,31. The magnitude of the phase jump, ÎÏhkl, anticipated due to a dislocation with Burgers vector b observed in a given hkl reflection is ÎÏhklâ=âb.qhkl. The scattering vectors associated with the (200), (020) and (002) reflections respectively are: q200â=â(2Ï/a) [200], q020â=â(2Ï/a) [020] and q002â=â(2Ï/a) [002], where a is the lattice parameter. This suggests that the defect in Fig. 3(D) is a dislocation with Burgers vector (a/3)[01-1]. Such a so-called stair-rod dislocation can be formed through the interaction of two Shockley partial dislocations32. For example in the present case the energetically favourable reaction (a/6)[21-1]â+â(a/6)[â21-1]âââ(a/3) [01-1] would produce a sessile dislocation with the observed Burgers vector.
The amplitude maps associated with reflections where phase jumps are observed show a local reduction in intensity at the defect position. Both (020) and (002) reflections show this reduction (white arrows in Fig. 3(D)), while no such feature is observed in the (200) reflection. This agrees with BCDI observations of âpipes of missing intensityâ at dislocation cores30. Indeed throughout the crystal, several further defects consistent with (a/3) <110> stair-rod dislocations32 can be identified (Supplementary Note 1.1, SI Fig. S6).
The ordering of larger defects in crystal C can be visualized by computing the von Mises stress33. Figure 3(E) shows von Mises stresses >500âMPa, greatly exceeding the yield strength of bulk gold (55â200âMPa)34. The arrangement of defects in lines is unexpected and differs from TEM observations of uniformly distributed (<100ânm) dislocation loops in FIB-milled copper12. The fact that we only observe sessile stair-rod dislocations is surprising. It suggests substantial evolution of the damage microstructure after the initial collision cascade with mobile dislocations escaping to the free surface and only sessile dislocations remaining.
It is interesting to consider the effect of these defects on the average strains induced by FIB milling. Figure 4 shows profiles of εxx, εyy, and εzz for crystals A and C plotted as a function of depth from the implanted surface. εyy in the implanted layer of crystal A is approximately twice as large as in crystal C, but has much less variation. This suggests that at higher Ga+ fluences larger defects, as well as the clustering of point defects10, act to relieve implantation-induced strains by localizing lattice distortion.
Discussion
Our findings show that every use of FIB to image or shape material causes large lattice distortions. Fundamental insight into the underlying damage mechanisms can be gained by combining coherent X-ray measurements with numerical calculations. Surprisingly, FIB-induced lattice strains are not confined to the ion-damaged layer, but can extend far into the material bulk, as visible in crystal C. This is further highlighted by measurements of crystal D into which a central hole was FIB-machined, and which exhibits large strains even far from the ion-damaged surfaces (Fig. 5). These extensive strains may explain the dramatic changes in mechanical properties caused by FIB milling11,35,36.
Our observations emphasise the need to actively consider the defects produced during FIB-based nano-fabrication. For example FIB-assisted deposition of a protective layer (usually Pt, W or C) is often used to prevent sample surface damage4. Our results suggest that even with such a protective layer substantial strains and damage should be expected. The reason is Ga+ damage formed during the initial stages of layer deposition before the layer thickness exceeds the Ga+ range. In the present gold crystals, FIB deposition of a Pt layer37,38,39 would cause defects similar to those in crystals A and B. E-beam-assisted deposition during the initial stages of layer growth could avoid this problem.
To realise the full potential of FIB, new strategies for controlling and minimising FIB damage must be developed. Current approaches include using a final low current milling step40,41,42, low energy ion milling43,44 or flash polishing45 in order to âcleanâ FIB-damaged surfaces. The techniques presented here enable 3D measurements of the complex strain fields caused by FIB-milling and would as such allow a detailed assessment of the effectiveness of these approaches. Conversely the ability of FIB to introduce large lattice strains with high spatial specificity presents an exciting opportunity for modifying material behaviour through strain engineering at the nano-scale.
Methods
Sample Manufacture
Gold crystals were prepared by dewetting a 20ânm thick gold layer, thermally evaporated onto a silicon substrate with a 2ânm titanium adhesion layer. The resulting crystals range from â100ânm to a few μm in size and show facets corresponding to {111} and {100} crystal planes (Fig. 1(A)). No FIB-milling was carried out in the vicinity of the unimplanted reference crystal. FIB-milling of crystals A, B and C was carried out at normal incidence, using a 30âkeV, 50âpA gallium ion beam and fluences of 4.2âÃâ104âions/μm2, 1.3âÃâ107âions/μm2 and 1.5âÃâ108âions/μm2 respectively. Crystal D was exposed to a fluence of 4.2âÃâ104âions/μm2 and a central, nominally 200ânm diameter, region to a fluence of 2.5âÃâ109âions/μm2. To allow reliable measurement of multiple reflections from crystals A, B, C and D, FIB was used to remove any other gold crystals within a 20âμm radius. Scanning electron micrographs of crystals A, B, C and D are shown in the SI Fig. S1. X-ray diffraction measurements were carried out 16 to 20 days after sample manufacture.
Ion Implantation Calculations
Ion implantation calculations used the âmonolayer collision - surface sputteringâ model in the Stopping Range of Ions in Matter (SRIM) code25. For the gold target a displacement energy of 44âeV, binding energy of 3âeV and surface energy of 3.8âeV were used46. Gallium ions were implanted at normal incidence with an energy of 30âkeV, gathering statistics over 105âions. Each ion was estimated to cause on average ~430 target displacements, of which ~30 were replacement collisions. The calculated sputtering rate was ~15.5 gold atoms per gallium ion. For crystal A the amount of material removed by sputtering was negligible. For crystals B and C, an estimated layer of thickness ~3ânm and ~40ânm respectively was removed. Custom MATLAB scripts were used to capture the receding surface effect due to sputtering. The calculated displacement damage and gallium concentration profiles for crystals A, B, C and D, plotted as a function of depth, are shown in the SI Fig. S3.
Experimental Measurements
Synchrotron X-ray micro-beam Laue diffraction at beamline 34-ID-E at the Advanced Photon Source (APS), Argonne National Lab, USA was used to determine the lattice orientations of gold crystals. This served to pre-align crystals for coherent X-ray diffraction measurements at beamline 34-ID-C at the APS. Measurements on the unimplanted reference crystal used an X-ray energy of 9.25âkeV, while diffraction patterns from crystals A, B, C and D were collected at 10.2âkeV. The X-ray beam was focussed to a size of 1.4âÃâ2.1âμm2 (hâÃâv) using KB mirrors. Placing the sample in the KB back-focal plane, within the central maximum of the focus, provides the planar wave front required for BCDI. Diffraction patterns were recorded on a Medipix2 area detector with a 256âÃâ256 pixel matrix and a pixel size of 55âÃâ55âμm2. For crystals A, B, C and D the detector was positioned 1.85âm from the sample and 3D coherent X-ray diffraction patterns (CXDP) were recorded by rotating the crystal through an angular range of 0.4° and recording an image every 0.0025° with 1âs exposure time. For the unimplanted reference crystal a sample-to-detector distance of 0.635âm was used and CXDPs were recorded by rotating through an angular range of 1.5° in 0.01° steps with 0.5âs exposure time. The sample to detector distances were chosen by starting at the distance required for the measurement of an oversampled diffraction pattern and then moving the detector further back until the diffraction pattern filled the detector matrix. To optimize the signal to noise of the CXDPs, multiple repeated scans of each reflection were performed. Repeated scans were then aligned to maximize their cross-correlation coefficient, and scans with a cross-correlation coefficient greater than 0.99 were summed to produce the CXDP for a specific reflection. For each crystal CXDPs from the following reflections were collected (the number of repeat scans that were averaged is noted in [] brackets): unimplanted reference: {111} [30]; crystal A: (1-11) [18], (11-1) [24], (200) [23], (020) [26], (002) [27]; crystal B: (-111) [16], (1-11) [9], (11-1) [14], (200) [11], (020) [16], (002) [17]; crystal C: (-111) [28], (1-11) [14], (11-1) [27], (200) [26], (020) [22], (002) [18]; and crystal D: (-111) [12], (1-11) [16], (11-1) [17], (200) [9], (020) [12], (002) [14]. Unfortunately the (-111) reflection of crystal A was physically inaccessible. Examples of the coherent diffraction patterns recorded from {111} reflections of all crystals are shown in SI Fig. S7.
Phase Retrieval
The phase retrieval algorithm used to recover the real-space complex electron density is adapted from published work30. Each 3D CXDP pattern was treated independently, using a guided phase retrieval approach with 20 random starts and 5 generations. For each generation 330 phase retrieval iterations were performed using Error Reduction and Hybrid-Input-Output algorithms. Trials using larger numbers of iterations showed no significant further evolution of the solution. Partial coherence effects were accounted for19, and the normalised mutual coherence functions, recovered for all reflections, are consistent with an almost fully coherent illumination. After the fifth generation a sharpness metric was used to select the three best estimates, which were then averaged to return the reconstructed complex electron density. Finally all reconstructions were transformed into an orthogonal laboratory reference frame with isotropic real-space pixel spacing. Agreement between the reconstructed crystal morphologies and scanning electron micrographs is excellent (SI Fig. S1). The normalised cross correlation coefficients, found when comparing the sample shape recovered from BCDI with SEM images, are 0.97, 0.98, 0.97 and 0.97 respectively for crystals A, B, C and D, when considering sample shape projected onto the plane of the substrate. Spatial resolution of the reconstructions was determined by taking the derivative of line profiles of the crystal-air-interface and fitting these with a Gaussian. For each reconstruction six profiles (2 in each spatial direction) were measured and the mean resolution value recorded.
3D Reconstruction of Lattice Displacements, Strains and Stresses
To recover the 3D lattice displacement field, u(r), of a given crystal, any phase wraps in the complex electron densities reconstructed from multiple crystal reflections were unwrapped using the algorithm developed by Cusack et al.47. Next all reconstructions were rotated into the same sample coordinate frame. The phase of the electron density reconstructed from a particular hkl peak, Ïhkl(r), is linked to the scattering vector qhkl and lattice displacement u(r) by Ïhkl(r)â=âqhkl.u(r). Thus each reconstruction provides a projection of u(r) along the corresponding qhkl. If 3 reflections with linearly independent qhkl are measured, u(r) can be reconstructed. Here 5 (crystal A) or 6 (crystals B, C and D) reflections with non-collinear q vectors were measured from each crystal. Thus the system of equations is over determined, and a least squares fit was used to calculate u(r). The symmetric Cauchy strain tensor, ε(r), is found by differentiating u(r). The strain uncertainty of our measurements, estimated from line profiles of ε(r) extracted from crystal A (Fig. 4), is ~10â4. Stresses were computed from ε(r) using anisotropic elastic constants for gold24.
Finite Element Calculations
Finite element simulations were performed in Abaqus 6.14, using the experimentally determined crystal morphology as a template for generation of the finite element mesh. Custom Matlab and Python scripts developed for this purpose are available upon request. A global seed size of 10ânm was used, based on mesh dependency studies that showed negligible improvements for finer mesh sizes. The resulting model for crystal A is shown in SI Fig. S2. Material properties were captured using anisotropic linear elastic constants for gold24. A uniform volumetric lattice strain, εv, was imposed within a 20ânm thick layer at the top face of crystal A to represent the effect of ion-implantation damage. εv =ââ3.15âÃâ10â3 provides a good match to the experimentally measured lattice displacement fields in crystal A. Displacements on the bottom surface of the crystal were fixed to capture the substrate effect.
Density Functional Theory Calculations
Ab initio density functional theory (DFT) calculations were performed of a mono-vacancy and of four different self-interstitial defect configurations in fcc gold (100 dumbbell, octahedral site, 110 crowdion and 110 dumbbell). Calculations were carried out in the Vienna ab initio simulation package (VASP)48,49,50,51 using the revised-TPSS exchange functional52,53, and included spin-orbit coupling. Spin-orbit coupling accounts for the band splitting and shape modifications of the 5d bands54,55,56. A plane wave energy cutoff of 450âeV was used and the outermost s- and d-electrons were treated as valence electrons. All samples were relaxed to a stress free condition with residual forces smaller than 0.01âeV/Ã . Formation energies and relaxation volumes were calculated by comparing the energies and volumes of a supercell containing each defect type with those of a perfect crystal supercell of similar size and using the same k-point mesh. Visualizations of the supercells used for these calculations are shown in SI Fig. S4.
The lattice constant and elastic constants were calculated using a 4 atom cubic unit cell. The equilibrium lattice constant is 4.075âà in good agreement with experiments57. The elastic constants were calculated using a finite differences scheme. We obtained c11â=â210.55âGPa, c12â=â168.11âGPa and c44â=â49.96âGPa, which compare well to experimental values at 0âK58. The elastic constants were used to correct both formation energy and relaxation volume of isolated defects according to the method suggested by Varvenne et al.59, which considers the elastic interactions of defects, and image forces due to the finite supercell size and periodic boundary conditions. The formation energies and relaxation volumes are listed in SI Table S1. The relaxation volume of a substitution gallium atom in gold was also calculated.
Data Availability
Diffraction data and selected computer codes used for data analysis and simulations can be obtained from the authors by contacting felix.hofmann@eng.ox.ac.uk.
Additional Information
How to cite this article: Hofmann, F. et al. 3D lattice distortions and defect structures in ion-implanted nano-crystals. Sci. Rep. 7, 45993; doi: 10.1038/srep45993 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
Brousseau, E. B., Dimov, S. S. & Pham, D. T. Some recent advances in multi-material micro- and nano-manufacturing. The International Journal of Advanced Manufacturing Technology 47, 161â180, doi: 10.1007/s00170-009-2214-5 (2009).
Inkson, B. J., Mulvihill, M. & Möbus, G. 3D determination of grain shape in a FeAl-based nanocomposite by 3D FIB tomography. Scripta Materialia 45, 753â758, doi: 10.1016/S1359-6462(01)01090-9 (2001).
Lasagni, F. et al. Three-dimensional characterization of âas-castâ and solution-treated AlSi12(Sr) alloys by high-resolution FIB tomography. Acta Materialia 55, 3875â3882, doi: 10.1016/j.actamat.2007.03.004 (2007).
Mayer, J., Giannuzzi, L. A., Kamino, T. & Michael, J. TEM Sample Preparation and FIB-Induced Damage. MRS Bulletin 32, 400â407, doi: 10.1557/mrs2007.63 (2007).
Giannuzzi, L. A. & Stevie, F. A. A review of focused ion beam milling techniques for TEM specimen preparation. Micron 30, 197â204, doi: 10.1016/S0968-4328(99)00005-0 (1999).
Uchic, M. D., Dimiduk, D. M., Florando, J. N. & Nix, W. D. Sample Dimensions Influence Strength and Crystal Plasticity. Science 305, 986â989, doi: 10.1126/science.1098993 (2004).
Greer, J. R. & De Hosson, J. T. M. Plasticity in small-sized metallic systems: Intrinsic versus extrinsic size effect. Progress in Materials Science 56, 654â724, doi: 10.1016/j.pmatsci.2011.01.005 (2011).
Robinson, M. T. & Torrens, I. M. Computer simulation of atomic-displacement cascades in solids in the binary-collision approximation. Physical Review B 9, 5008â5024 (1974).
Yi, X. et al. Direct observation of size scaling and elastic interaction between nano-scale defects in collision cascades. Europhysics Letters 110, 36001 (2015).
Hofmann, F. et al. Lattice swelling and modulus change in a helium-implanted tungsten alloy: X-ray micro-diffraction, surface acoustic wave measurements, and multiscale modelling. Acta Materialia 89, 352â363, doi: 10.1016/j.actamat.2015.01.055 (2015).
Shim, S., Bei, H., Miller, M. K., Pharr, G. M. & George, E. P. Effects of focused ion beam milling on the compressive behavior of directionally solidified micropillars and the nanoindentation response of an electropolished surface. Acta Materialia 57, 503â510, doi: 10.1016/j.actamat.2008.09.033 (2009).
Kiener, D., Motz, C., Rester, M., Jenko, M. & Dehm, G. FIB damage of Cu and possible consequences for miniaturized mechanical tests. Materials Science and Engineering: A 459, 262â272, doi: 10.1016/j.msea.2007.01.046 (2007).
Giannuzzi, L. A., Geurts, R. & Ringnalda, J. 2 keV Ga+ FIB Milling for Reducing Amorphous Damage in Silicon. Microscopy and Microanalysis 11, 828â829 (2005).
Yu, J., Liu, J., Zhang, J. & Wu, J. TEM investigation of FIB induced damages in preparation of metal material TEM specimens by FIB. Materials Letters 60, 206â209, doi: 10.1016/j.matlet.2005.08.018 (2006).
Salvati, E., Sui, T., Lunt, A. J. G. & Korsunsky, A. M. The effect of eigenstrain induced by ion beam damage on the apparent strain relief in FIB-DIC residual stress evaluation. Materials & Design 92, 649â658, doi: 10.1016/j.matdes.2015.12.015 (2016).
Sabaté, N. et al. FIB-based technique for stress characterization on thin films for reliability purposes. Microelectronic Engineering 84, 1783â1787, doi: 10.1016/j.mee.2007.01.272 (2007).
Korsunsky, A. M., Sebastiani, M. & Bemporad, E. Residual stress evaluation at the micrometer scale: Analysis of thin coatings by FIB milling and digital image correlation. Surface and Coatings Technology 205, 2393â2403, doi: 10.1016/j.surfcoat.2010.09.033 (2010).
Barabash, R. I. et al. Mapping strain gradients in the FIB-structured InGaN/GaN multilayered films with 3D X-ray microbeam. Materials Science and Engineering: A 528, 52â57, doi: 10.1016/j.msea.2010.04.045 (2010).
Clark, J. N., Huang, X., Harder, R. & Robinson, I. K. High-resolution three-dimensional partially coherent diffraction imaging. Nature Communications 3, 993 (2012).
Robinson, I. & Harder, R. Coherent X-ray diffraction imaging of strain at the nanoscale. Nature Materials 8, 291â298 (2009).
Pfeifer, M. A., Williams, G. J., Vartanyants, I. A., Harder, R. & Robinson, I. K. Three-dimensional mapping of a deformation field inside a nanocrystal. Nature 442, 63â66 (2006).
Newton, M. C., Leake, S. J., Harder, R. & Robinson, I. K. Three-dimensional imaging of strain in a single ZnO nanorod. Nat Mater 9, 120â124 (2010).
Robinson, I. Nanoparticle Structure by Coherent X-ray Diffraction. Journal of the Physical Society of Japan 82, 021012, doi: 10.7566/JPSJ.82.021012 (2012).
Hiki, Y. & Granato, A. V. Anharmonicity in Noble Metals; Higher Order Elastic Constants. Physical Review 144, 411â419 (1966).
Ziegler, J. F., Ziegler, M. D. & Biersack, J. P. SRIMâThe stopping and range of ions in matter (2010). Nuclear Instruments and Methods in Physics Research, Section B: Beam Interactions with Materials and Atoms 268, 1818â1823 (2010).
Sonnenberg, K. & Dedek, U. Migration energy of single vacancies in gold. Radiation Effects 61, 175â178, doi: 10.1080/00337578208229930 (1982).
Ehrhart, P., Carstanjen, H. D., Fattah, A. M. & Roberto, J. B. Diffuseâscattering study of vacancies in quenched gold. Philosophical Magazine A 40, 843â858, doi: 10.1080/01418617908234878 (1979).
Atomic Defects in Metals Au (Springer-Verlag Berlin Heidelberg, 1991).
Takahashi, Y. et al. Bragg x-ray ptychography of a silicon crystal: Visualization of the dislocation strain field and the production of a vortex beam. Physical Review B 87, 121201 (2013).
Clark, J. N. et al. Three-dimensional imaging of dislocation propagation during crystal growth and dissolution. Nature Materials 14, 780â784, doi: 10.1038/nmat4320 (2015).
Ulvestad, A. et al. Topological defect dynamics in operando battery nanoparticles. Science 348, 1344â1347, doi: 10.1126/science.aaa1313 (2015).
Hirth, J. P. & Lothe, J. Theory of Dislocations. 2 edn (Wiley, 1982).
Hill, R. The mathematical theory of plasticity(Clarendon Press, 1950).
Espinosa, H. D., Prorok, B. C. & Peng, B. Plasticity size effects in free-standing submicron polycrystalline FCC films subjected to pure tension. Journal of the Mechanics and Physics of Solids 52, 667â689, doi: 10.1016/j.jmps.2003.07.001 (2004).
Bei, H., Shim, S., Miller, M. K., Pharr, G. M. & George, E. P. Effects of focused ion beam milling on the nanomechanical behavior of a molybdenum-alloy single crystal. Applied Physics Letters 91, 111915â111913 (2007).
Lee, S.-W., Mordehai, D., Rabkin, E. & Nix, W. D. Effects of focused-ion-beam irradiation and prestraining on the mechanical properties of FCC Au microparticles on a sapphire substrate. Journal of Materials Research 26, 1653â1661, doi: 10.1557/jmr.2011.221 (2011).
Tao, T., Wilkinson, W. & Melngailis, J. Focused ion beam induced deposition of platinum for repair processes. Journal of Vacuum Science & Technology B 9, 162â164, doi: 10.1116/1.585279 (1991).
Puretz, J. & Swanson, L. W. Focused ion beam deposition of Pt containing films. Journal of Vacuum Science & Technology B 10, 2695â2698, doi: 10.1116/1.586028 (1992).
Telari, K. A. et al. Characterization of platinum films deposited by focused ion beam-assisted chemical vapor deposition. Journal of Vacuum Science & Technology B 20, 590â595, doi: 10.1116/1.1458958 (2002).
Gong, J., Benjamin Britton, T., Cuddihy, M. A., Dunne, F. P. E. & Wilkinson, A. J. ãaã Prismatic, ãaã basal, and ãc+aã slip strengths of commercially pure Zr by micro-cantilever tests. Acta Materialia 96, 249â257, doi: 10.1016/j.actamat.2015.06.020 (2015).
Liu, L., Chen, Z., Liu, C., Wu, Y. & An, B. Micro-mechanical and fracture characteristics of Cu6Sn5 and Cu3Sn intermetallic compounds under micro-cantilever bending. Intermetallics 76, 10â17, doi: 10.1016/j.intermet.2016.06.004 (2016).
Žagar, G., Pejchal, V., Mueller, M. G., Michelet, L. & Mortensen, A. Fracture toughness measurement in fused quartz using triangular chevron-notched micro-cantilevers. Scripta Materialia 112, 132â135, doi: 10.1016/j.scriptamat.2015.09.032 (2016).
McCaffrey, J. P., Phaneuf, M. W. & Madsen, L. D. Surface damage formation during ion-beam thinning of samples for transmission electron microscopy. Ultramicroscopy 87, 97â104, doi: 10.1016/S0304-3991(00)00096-6 (2001).
Thompson, K. et al. In situ site-specific specimen preparation for atom probe tomography. Ultramicroscopy 107, 131â139, doi: 10.1016/j.ultramic.2006.06.008 (2007).
Dawson, K. & Tatlock, G. J. Preparation of micro-foils for TEM/STEM analysis from metallic powders. Micron 74, 54â58, doi: 10.1016/j.micron.2015.04.006 (2015).
Jung, P. Average atomic-displacement energies of cubic metals. Physical Review B 23, 664â670 (1981).
Cusack, R. & Papadakis, N. New Robust 3-D Phase Unwrapping Algorithms: Application to Magnetic Field Mapping and Undistorting Echoplanar Images. NeuroImage 16, 754â764, doi: 10.1006/nimg.2002.1092 (2002).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Physical Review B 54, 11169â11186 (1996).
Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal amorphous-semiconductor transition in germanium. Physical Review B 49, 14251â14269 (1994).
Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Physical Review B 47, 558â561 (1993).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science 6, 15â50, doi: 10.1016/0927-0256(96)00008-0 (1996).
Perdew, J. P., Ruzsinszky, A., Csonka, G. I., Constantin, L. A. & Sun, J. Workhorse Semilocal Density Functional for Condensed Matter Physics and Quantum Chemistry. Physical Review Letters 103, 026403 (2009).
Sun, J. et al. Self-consistent meta-generalized gradient approximation within the projector-augmented-wave method. Physical Review B 84, 035117 (2011).
Corso, A. D. & Conte, A. M. Spin-orbit coupling with ultrasoft pseudopotentials: Application to Au and Pt. Physical Review B 71, 115106 (2005).
Christensen, N. E. & Seraphin, B. O. Relativistic Band Calculation and the Optical Properties of Gold. Physical Review B 4, 3321â3344 (1971).
Rangel, T. et al. Band structure of gold from many-body perturbation theory. Physical Review B 86, 125125 (2012).
Gray, D. E. American Institute of Physics handbook(McGraw-Hill 1963).
Neighbours, J. R. & Alers, G. A. Elastic Constants of Silver and Gold. Physical Review 111, 707â712 (1958).
Varvenne, C., Bruneval, F. & Marinica, M.-C. & Clouet, E. Point defect modeling in materials: Coupling ab initio and elasticity approaches. Physical Review B 88, 134102 (2013).
Acknowledgements
F.H. acknowledges funding from the John Fell fund (122/643) and the Royal Society (RG130308). E.T. acknowledges funding through an EPSRC Early Career Fellowship (EP/N007239/1). J.N.C. acknowledges financial support from the Volkswagen Foundation. Diffraction experiments used the Advanced Photon Source, a US Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357. This work was supported by the Australian Research Council Centre of Excellence in Advanced Molecular Imaging (CE140100011) www.imagingcoe.org. We also acknowledge funding from the United Kingdom Engineering and Physical Sciences Research Council via programme grants EP/H018921/1 and EP/I022562/1. This work has been partly carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014-2018 under grant agreement No 633053 and from the RCUK Energy Programme (grant number EP/I501045). The views and opinions expressed herein do not necessarily reflect those of the European Commission. Work at Brookhaven National Laboratory was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract DE-SC00112704.
Author information
Authors and Affiliations
Contributions
F.H. designed the project. J.N.C., I.K.R. and C.E.B. made the samples. F.H., R.J.H., N.W.P. and W.L. carried out the experiments. F.H. carried out experimental analysis with input from R.J.H. and J.N.C. E.T. carried out the finite element modelling. P.W.M. carried out D.F.T. calculations. F.H. and E.T. wrote the paper. All authors edited the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Hofmann, F., Tarleton, E., Harder, R. et al. 3D lattice distortions and defect structures in ion-implanted nano-crystals. Sci Rep 7, 45993 (2017). https://doi.org/10.1038/srep45993
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep45993
This article is cited by
-
Atomic resolution coherent x-ray imaging with physics-based phase retrieval
npj Computational Materials (2024)
-
Three-dimensional strain imaging of irradiated chromium using multi-reflection Bragg coherent diffraction
npj Materials Degradation (2022)
-
AutoPhaseNN: unsupervised physics-aware deep learning of 3D nanoscale Bragg coherent diffraction imaging
npj Computational Materials (2022)
-
Revealing nano-scale lattice distortions in implanted material with 3D Bragg ptychography
Nature Communications (2021)
-
Three-dimensional strain dynamics govern the hysteresis in heterogeneous catalysis
Nature Communications (2020)