Microscopic crystallographic analysis of dislocations in molecular crystals

Main

Across materials classes, dislocations in crystals have profound effects on properties—from the mechanical properties of metals7,8 to deleterious effects on performance in semiconductors9 as well as the advantageous modulation of charge transport properties in oxides10. These findings have relied heavily on electron microscopy analysis. In characteristically anisotropic molecular crystals, with molecules rather than atoms at the lattice sites, modelling and experiments have outlined how the structures at particular dislocation cores govern reactions in energetic materials11 and pharmaceutical tableting12, as well as producing trap states in optoelectronics5. Yet, there is a critical knowledge gap in analysing these dislocations at the nanoscale.

Dislocations are, in effect, a planar cut through a crystal followed by atomic displacements in a particular direction, the Burgers vector B (Fig. 1a). The singular displacement direction creates a linear defect, and this dislocation line vector u coincides with the location of maximum lattice distortion. Together, these vectors describe the operative slip systems and the character of the dislocation (that is, edge, screw or mixed). For edge dislocations, Bu, whereas for screw dislocations, Bu (Fig. 1a). The classical electron microscopy approach to dislocation analysis involves the combination of diffraction and diffraction contrast imaging to evaluate the behaviour of lattice planes and associated diffraction vectors ghkl in the vicinity of a dislocation. For ghkl ∙ (B × u) = 0, a unique invisibility criterion is established for ghkl ∙ B= 0. By exploring a range of diffraction conditions, the Burgers vector direction can be determined unambiguously.

Fig. 1: Low-dose dislocation analysis by SED.
figure 1

a, Illustration of the distortion of planes due to edge, mixed and screw dislocations, with no distortion for planes corresponding to diffraction vectors ghkl at the invisibility criterion (ghkl ∙ B = 0). b, Schematic of an edge dislocation core in a p-terphenyl crystal with Burgers vector B = [010]. The blue parallelograms represent molecules in a plane below those in orange. A black rectangle or line marks the ab plane of the p-terphenyl lattice in the plane below (blue) the dislocation core (). c, A bend contour appears where the sample curvature brings planes (hkl) to the exact Bragg condition, with displacements in the bend contour due to abrupt changes in the orientation of planes around the dislocation core, marked by the dislocation line u. The electron beam is scanned in (x, y), enabling the reconstruction of bend contours from multiple diffraction vectors ghkl in parallel from two-dimensional diffraction patterns (kx, ky) recorded at each probe position to give a four-dimensional dataset (x, y, kx, ky). The sample bending is exaggerated for visual effect. The axes are shown as in an experiment, with the scan and diffraction cameras not generally aligned, requiring rotation calibration (Methods).

Full size image

Although scanning probe microscopy13, X-ray topography14 and modelling3 have supported dislocation characterization in molecular crystals, these techniques lack the spatial resolution to image individual dislocations at high density. Atomic-resolution imaging of dislocations15, including imaging at electron fluences of <200 e Å−2 in halide perovskites16, and three-dimensional imaging modes17,18 have been demonstrated, but none of the existing approaches support dislocation analysis at ~10 e Å−2 as required for the electron microscopy of many organic molecular crystals6. Moreover, the complex molecular packing of organic crystals in projection (Fig. 1b) generally precludes the lattice imaging of dislocations beyond a few high-symmetry orientations19,20. Early work on real-space crystallography supported the inference of plausible Burgers vector directions21,22,23, but these analyses have only been reported for molecular crystals with relatively high electron-beam stability. Our aim is to develop a low-dose method that is compatible with diverse sample preparations and supports practical application across varied instrumentation and sample environments.

Scanning electron diffraction (SED), a type of four-dimensional scanning transmission electron microscopy (4D-STEM) with a nearly parallel, low-convergence-angle probe24, suggests a route to low-dose electron microscopy for dislocation analysis. SED combines real- and reciprocal-space data acquisition simultaneously and can be operated at sufficiently low fluences to detect mosaicity in peptide crystals25 and molecular packing in organic semiconductors26,27,28. We now establish a methodology for single-exposure SED analysis of in-plane dislocations in thin, electron-transparent organic crystals at fluences as low as 5 e Å−2. We apply this approach to aromatic molecules, long-chain hydrocarbons and hydrogen-bonding crystals as key optoelectronic, wax and pharmaceutical models, respectively.

SED enables the parallel acquisition of many diffraction vectors by recording a two-dimensional diffraction pattern at each probe position in a scan (Fig. 1c). By selecting individual ghkl, we construct virtual dark-field (VDF) images. Bending of the electron-transparent crystals deposited onto electron microscopy grids results in the appearance of bands (denoted as bend contours), marking the positions in the sample that are exactly at the Bragg condition. For planes that are distorted by the dislocation (ghkl ∙ B ≠ 0), a break in the otherwise-continuous bend contour appears at the dislocation line u (refs. 21,22,29), akin to features observed in convergent-beam electron diffraction30. This break corresponds to a displacement of the bend contour along u (Supplementary Note 1). For ghkl at the invisibility criterion, the contour remains continuous. By contrast, planar and point defects would exhibit a change in all ghkl values or, if correlated and symmetry breaking, would introduce additional ghkl.

In SED, all the recorded ghkl values in the diffraction plane of sufficient intensity can be used to construct VDF images for analysis as a function of the azimuthal angle in the diffraction plane φ (Supplementary Fig. 1). A geometric model of the tilt of lattice planes across the dislocation core describes the angle-dependent contour displacement f(φ) for in-plane dislocations (Supplementary Figs. 2–5):

$$fleft(varphi right)=Aarctan left(B{cos }^{2}(varphi -C)right),$$
(1)

where A, B and C are the fitting coefficients related to the effective radius of curvature, Burgers vector and azimuthal orientation of the diffraction pattern on the detector, respectively. The displacement is zero at the invisibility criterion and maximum along the Burgers vector. The crystal orientation given by the diffraction pattern together with the Burgers vector direction (equation (1)) provides the necessary information to identify the operative slip system (Supplementary Note 1), as well as the dislocation character by inspection of the relative directions of B and u.

Figure 2 presents a first demonstration applied to p-terphenyl (C18H14; CCDC 1269381) and anthracene (C14H10; CCDC 1269381). Both crystallize in the P21/a space group with herringbone packing. Although p-terphenyl is a relatively stable benchmark organic crystal22 (with a critical fluence (CF) > 330 e Å−2; Supplementary Fig. 6), anthracene is not (CF ≈ 20 e Å−2; Supplementary Fig. 7). Figure 2b shows an annular dark-field (ADF) STEM image of a p-terphenyl film with a dense dislocation network retrieved from multiple VDF images (Methods). The length of individual dislocations, for example, >2 μm, indicates that they are predominantly in plane (<3° out of plane for film thicknesses of 10–100 nm). The dislocation density is 0.32–3.2 × 1010 cm−2 (Supplementary Table 1). The area-averaged diffraction pattern was indexed to [001], as expected22. Figure 2c presents the recorded bend contour displacements as a function of azimuthal angle φ (relative to the horizontal axis of the diffraction pattern; Supplementary Note 1) in a polar plot. A fit to equation (1) is overlaid. Supplementary Note 2, Supplementary Figs. 8 and 9 and Supplementary Table 2 provide details of the curve-fitting process and parameters. Although sizable displacements are visible for ghkl, for example, g210 (d spacing, 3.28 Å), none of them are visible for g200 or ({{bf{g}}}_{bar{2}00}) (d spacing, 4.05 Å) (Fig. 2c). On the basis of the in-plane assumption, we show that B = [010] (equation (6) in the Supplementary Information), with mixed dislocation character; here the operative slip system is identified as [010](001), again as expected22. Supplementary Fig. 10 shows two further dislocations from the same area, similarly assigned as mixed, B = [010] dislocations.

Fig. 2: Dislocation analysis in organic optoelectronic materials.
figure 2

a, p-Terphenyl molecule and p-terphenyl unit cell. b, ADF-STEM image, dislocation network (green lines tracing the diffraction contrast features of dislocations) and average diffraction pattern extracted from a single SED (4D-STEM) dataset. The diffraction pattern indexes to the [001] zone axis. c, Polar plot of the bend contour displacements at the solid dislocation line highlighted in b for a series of ghkl values as a function of φ. The red trace marks a fit to equation (1), with the measured displacements in blue. Two VDF images corresponding to the displacements marked as (i) and (ii) depicting bend contours at the ghkl ∙ B = 0 and ghkl ∙ B ≠ 0 conditions, respectively. The indices hkl are marked on the individual VDF images. d, p-Terphenyl unit cell oriented to match the sample orientation. The magenta arrow marks the Burgers vector B = [010]. e, Anthracene molecule and unit cell. f, ADF-STEM image and the overlaid dislocation network extracted from the image contrast. The image in the inset shows the average diffraction pattern of the entire field of view indexed to the [101] zone axis. g, Polar plot of the bend contour displacement as a function of φ for the solid green dislocation line highlighted in f. The red trace marks a fit to equation (1) with the measured displacements in blue. h, Anthracene unit cell oriented to match the sample orientation. The magenta arrow marks the Burgers vector B = [010].

Source data

Full size image

Critically, the fitting approach (equation (1)) means that SED-based dislocation analysis is amenable even for a few ghkl values, for example, for a limited area with few bend contours or crystals tilted off high-symmetry zone axis orientations (Supplementary Figs. 11–14). In p-terphenyl, ghkl > 1 Å−1 was recorded, presenting many spacings for identifying contours (Supplementary Fig. 13). Figure 2f presents an ADF-STEM image containing a broad bend contour (top) and corresponding dislocation network in anthracene (confirmed by VDF and diffraction pattern inspection; Supplementary Figs. 15–18). No single VDF showed a zero-displacement bend contour (ghkl ∙ B = 0). Nevertheless, fitting of the polar plot is unambiguous (Fig. 2g) and supports the determination of B = [010] for this mixed-type dislocation (Fig. 2h). Not all analyses are sufficiently constrained for Burgers vector determination, but dislocation character and handedness can be retrieved and combined with the knowledge of molecular packing to produce useful analyses (Supplementary Figs. 19 and 20). The analysis of dislocations in spin-coated anthracene films demonstrates the success of the approach for varied sample preparations (Supplementary Fig. 21).

The versatility of our approach is further demonstrated in crystal plates (third sample preparation; Fig. 3) of the model pharmaceutical compound theophylline (C7H8N4O2; CF ≈ 26–36 e Å−2 for the stable form II (ref. 6)), and a paraffin n-hentriacontane (C31H64; measured CF ≈ 6 e Å−2, similar to other paraffins19). Their defects offer insights into hydration/rehydration pathways on drug storage31 and diffusion pathways in leaf waxes32. Theophylline exhibits a number of polymorphs and hydrates, and we first show an example from theophylline anhydrous form II (Fig. 3a–d). The average diffraction pattern was indexed to [141] (Fig. 3c). The dislocation is short, suggesting that the dislocation line is not in plane, testing our model’s generalization. Still, VDF analysis of the available ghkl values enables fitting to equation (1) (Fig. 3d), and the geometrically recovered B aligns with B = [102] in projection, exhibiting mixed character. Additional screw-type dislocations with B along 〈100〉 were identified in the monohydrate form (form M; Supplementary Fig. 22) and metastable anhydrous form III (Supplementary Figs. 23 and 24). Given the differences in hydrogen bonding in the resulting dislocation cores, these dislocations may influence the dehydration pathway from form M to form II (ref. 33).

Fig. 3: Low-dose dislocation analysis in theophylline and paraffin crystals.
figure 3

a, Theophylline molecule and unit cell of form II. b, ADF-STEM image with a single dislocation line overlaid in green. c, Corresponding electron diffraction pattern, indexed to the [141] zone axis orientation. d, Burgers vector determination by fitting a polar plot of the bend contour displacement as a function of φ to equation (1) for the dislocation line highlighted in b. The measured displacements are shown in blue, with the fit shown in red. A dashed magenta line marks the experimentally determined Burgers vector direction ({bf{B}}=[1bar{16}4]), which aligns in projection to within ~3° with the assigned vector (B = [102]). The unit cell is depicted at the same orientation as given by the diffraction pattern in c as well as along the [010] direction to better visualize the assigned Burgers vector. e, n-Hentriacontane (C31H64) molecule and unit cell. f, ADF-STEM image with dislocation lines marked in green. g, Corresponding electron diffraction pattern, indexed to the [001] zone axis orientation. h, Burgers vector determination by fitting a polar plot of the bend contour displacement as a function of φ to equation (1) for the dislocation line highlighted in f. The measured displacements are shown in blue, with the fit shown in red. The unit cell is depicted at the same orientation as given by the diffraction pattern in g. The magenta arrow follows the long axis of the lobes in the polar plot, giving the Burgers vector direction as B = [010].

Source data

Full size image

Figure 3e depicts the n-hentriacontane unit cell34. The plate crystal in Fig. 3f was oriented along the [001] zone axis (Fig. 3g). Here the dislocation was unambiguously identified with the Burgers vector direction along B = [010]. The dislocation line illustrates a full range of edge, screw and mixed character along its curved path, with the length of the dislocation line suggesting that the dislocation lies along (001), adding to reported screw dislocations3. Supplementary Fig. 25 presents additional pure screw dislocations in wax crystals with equivalent Burgers vectors. Going beyond edge dislocations imaged in helium-cooled paraffins19, in-plane mixed and screw dislocations can now be identified.

VDF images additionally support the analysis of the ‘twisting’ of bend contours at dislocations, a marker of the handedness of the screw component29,30 (Supplementary Note 3 and Supplementary Fig. 19). By combining the total Burgers vector direction (180° ambiguity) and the handedness for screw and mixed dislocations, we identify the handedness and the edge-component direction. Figure 4 demonstrates this approach for p-terphenyl and theophylline form IIIb crystals, depicting networks with aligned screw components, consistent with repulsion between such interacting screw components. Supplementary Figs. 20, 22, 24 and 25 present additional examples from anthracene, theophylline forms M and IIIb and wax crystals, respectively, also pointing to features that may pin the dislocations of opposite handedness and inviting experimentally constrained evaluation of dislocation dynamics at densities only accessible with nanoscale probes. In specific cases, the magnitude of the Burgers vector can also be determined using the established relationship29,30 ghkl ∙ B = n (Supplementary Figs. 22–24).

Fig. 4: Sign analysis of interacting dislocations.
figure 4

a,b, Summed VDF images displaying bend contours crossing dislocation networks in p-terphenyl (a) and theophylline form IIIb (b). c,d, VDF images showing the selected pairs of bend contours at scattering vectors ghkl that are parallel or nearly parallel to the dislocation lines, that is, along the screw component of the Burgers vector, for p-terphenyl (c) and theophylline form IIIb (d). The green dashed lines mark the dislocation cores (lines u), and the red and yellow arrows indicate the scattering vectors ghkl. The grey solid line in d outlines a grain boundary. The white arrows mark the screw components determined using the Cherns and Preston rules30. The magenta arrows mark the full Burgers vectors and the blue arrows mark the edge components, determined by geometric self-consistency with each screw component and Burgers vector direction from the ghkl ∙ Bscrew = 0 analysis (Supplementary Fig. 19).

Source data

Full size image

In developing this approach, we also note some limitations. Measuring displacements depends on the precision in locating the centre of the bend contour. Consequently, samples containing multiple finely spaced dislocations (Supplementary Figs. 26 and 27) present challenges at finite resolution or where broad or tortuous bend contours preclude precise contour-centre localization. Excluding unusable ghkl values is possible for fitting such cases. Further extensions may also be possible by analysing contours crossing several dislocations, though prior knowledge of contributing Burgers vectors may be required, and we cannot rule out deviations from linear, elastic behaviour in such regions. During beam rastering, slight reorientation of the sample can occur, which may result in shifts in bend contours mid-scan. However, these features are readily distinguished by their uniform behaviour for all ghkl values (Supplementary Fig. 28). In more beam-tolerant systems (or where sample cooling supports multiple exposures), controlled tilting strategies may support Burgers vector determination from non-coplanar ghkl. Similarly, where sample preparation can be adjusted, controlling the curvature may support the analysis of additional model parameters. Acceleration of data processing would further enable selective measurements from pristine adjacent areas for targeted contour evaluation.

In summary, we have demonstrated in-plane dislocation analysis in varied organic molecular crystals using low-dose, single-exposure SED measurements. Introducing a geometric model for fitting, although imposing limited assumptions on slowly varying curvature, and component sign analysis enables the determination of Burgers vectors with no presumed knowledge of the crystal structure and from highly limited sampling of ghkl. Observed dislocations exhibit Burgers vector directions that lie predominantly in planes containing weak intermolecular interactions between layers of molecules without cutting across the molecular structure—structures that are distinct from inorganic solids. The generalization of this approach provides a means for the routine analysis of dislocations in molecular crystals and other materials with CFs down to a few e Å−2 (for example, using bend contours in framework materials35,36 and halide perovskites37), driving forward research into which intermolecular interactions at a dislocation core (dictated by the Burgers vector) are active in mechanical, optoelectronic and chemical properties and the development of methods for inhibiting or designing dislocation formation.

Methods

Sample preparation

Thin films of p-terphenyl and anthracene were prepared by solvent evaporation from xylene solutions22. Briefly, p-terphenyl (99.5%, Merck) or anthracene (99.0%, Merck) powder was dissolved in xylene to obtain a clear solution (~0.45–0.50 mg ml−1). Here 0.4 ml of the xylene solution was left to evaporate slowly on the surface of deionized water in a glass dish covered with Al foil to slow the rate of evaporation. Single-crystal films were obtained after 2 days of evaporation followed by half a day of aging. Interference colours were visible indicators of crystal thickness. Films that appeared blue or gold were sufficiently thin (120–200 nm) for electron microscopy38, and these films were selected for the SED experiments. The single-crystal films were transferred onto lacey-carbon-film transmission electron microscopy (TEM) support grids (EM Resolutions) by placing the grid under the film and lifting to catch the film on the TEM grid. Before the transfer process, the TEM grids were treated with an oxygen–argon (25:75) plasma for 10 s using an HPT-100 plasma treatment system (Henniker Plasma). Grids were allowed to dry in ambient air. Spin-coated anthracene films were prepared on SiNx-membrane-window TEM grids (3 × 3 SiNx windows and 30 nm thickness; Norcada) exposed to an argon plasma for a minute. For spin coating, a 5–6 mg ml−1 solution of anthracene in chloroform was used. A SiNx membrane window was fixed on a custom chuck using tape and 30 μl of the casting solution was dropped on the membrane window at 200 r.p.m. for 10 s. The spinning was then increased to 3,000 r.p.m. and the film was spun for a further minute. Finally, the SiNx grid was annealed at 80 °C on a hot plate for 15 min.

Single-crystal films of theophylline were prepared by the solvent evaporation approach39,40. Briefly, theophylline powder was dissolved in nitromethane followed by heating to ~60 °C to obtain a saturated solution. The solution was then drop cast onto a holey-carbon-film TEM support grid (EM Resolutions) and left to evaporate. Highly faulted particles of partially hydrated form II plates were produced by immersion in water. A sample of theophylline monohydrate (form M) was prepared by evaporative crystallization from either ethanol:water or nitromethane:water solutions at 3:2 molar ratio followed by heating to ~55 °C. The solution was allowed to cool slowly to ~35 °C under continuous stirring for a period of 3 h. The solution was then drop cast onto a holey-carbon-film TEM support grid and left to evaporate. Form III (indexed to the form IIIb structure33) was identified on a TEM grid with theophylline form M (recrystallized from ethanol:water) and then immersed in liquid nitrogen. Data on form IIIb were similarly acquired under liquid-nitrogen cooling (<110 K) using a Fischione cryo-transfer holder. All other reported data were acquired under ambient conditions. Single-crystal films of n-hentriacontane and 1-triacontanol were similarly prepared by solvent evaporation41 by drop casting a supersaturated hexane solution (2 μl) onto a continuous-carbon-film TEM support grid (EM Resolutions).

TEM

Electron diffraction pattern series as a function of cumulative electron fluence for p-terphenyl and anthracene were recorded using an FEI Titan3 Themis 300 (X-FEG high-brightness electron source, operated at 300 kV) microscope in a standard parallel-beam TEM configuration (Supplementary Note 4). For p-terphenyl, an electron-beam flux of 2.5 e Å−2 s−1 was used, whereas for anthracene, an electron-beam flux of 0.085 e Å−2 s−1 was used. These values were determined from the current read-out at the viewing screen, calibrated using a Faraday cup. A time series of the selected-area electron diffraction patterns were recorded for each compound across multiple crystals and multiple areas on the TEM grid.

SED

SED data were acquired on a JEOL ARM300CF instrument (electron Physical Sciences Imaging Centre, Diamond Light Source) equipped with a high-resolution pole piece, a cold field emission gun, aberration correctors in both probe-forming and image-forming optics, and a four-chip Merlin/Medipix pixelated electron-counting STEM detector. The instrument was operated at 300 kV. Selected samples were analysed at 200 kV (identified in the corresponding figure captions). The electron optics were configured for nanobeam diffraction by switching off the aberration corrector in the probe-forming optics and adjusting the condenser lens system to produce a convergence semiangle of 0.8 mrad using a 10 μm condenser aperture. At 300 kV, this produces a 3 nm diffraction-limited probe diameter ddiff = 1.22λ/α, where λ is the electron de Broglie wavelength and α is the convergence semiangle. The probe current was measured using a Faraday cup as ~1 pA, and the exposure time at each probe position was set as 1 ms. The electron fluence was approximately 8.8 e Å−2 at 300 kV or 5.4 e Å−2 at 200 kV in a single scan, assuming a disc-like probe with a diameter equal to ddiff. All the SED measurements were conducted over a scan size of 256 × 256 probe positions. Image and diffraction calibration data, including reference calibration of residual elliptical distortion in the diffraction plane, were acquired using a gold diffraction cross-grating with a period of 500 nm (Ted Pella). These cross-gratings contain a small fraction of Pd (AuPd); compared with an evaporated Au sample, we estimate ~1% error in Å−1 per pixel calibrations at the camera lengths used in this work when calibrating with respect to the pure Au crystal structure (with no effect on elliptical distortion calibration). To minimize these effects on our analysis, we used full pattern matching for indexation (see the ‘Data analysis’ section). The relative rotation between the diffraction pattern and raster pattern was calibrated using standard MoO3 crystals (Agar Scientific). Calibration data were acquired under identical conditions as the molecular crystal samples. When accounting for this calibration, the diffraction patterns are rotated, resulting in the rotation of a ‘cross’ arising from a gap in the pixel readouts between the four quadrants of the Merlin/Medipix detector.

Data analysis

SED data were processed, aligned and calibrated using Pyxem-0.11.0 and supporting tools from the HyperSpy package (v. 1.6.5)42 according to previously reported procedures43. The direct beam in each diffraction pattern was aligned and moved to the pattern centre using a cross-correlation function in Pyxem followed by the application of an affine transformation matrix to correct the elliptical distortion and elimination of the rotation offset between the SED raster pattern and diffraction data. In some datasets, the orientation of the cross pattern (arising from the gap in pixel read-out between quadrants of the Merlin/Medipix detector) varies because the beam was positioned on different quadrants across different experiments. The ADF images were formed by integrating the diffraction pattern at each probe position between an inner radius of 0.12 Å−1 (2θ = 2.7 mrad) and an outer radius of 1 Å−1 (2θ = 22.4 mrad) to produce an image dominated by diffraction contrast. Average electron diffraction patterns were obtained by taking the mean intensity from all the diffraction patterns contributing to an image. Diffraction patterns are presented as the square root of the recorded intensity (applied in ImageJ software, v. 2.16.0) for an improved visualization of the low- and high-intensity features.

Diffraction patterns were indexed using CrystalMaker software (v. 10.7.3). The simulated kinematical electron diffraction patterns were generated from the CIF files for p-terphenyl (CCDC 1269381), anthracene (CCDC 1103062), theophylline (CCDC 128707 for form II, CCDC 183780 for form M and ref. 33 for metastable form IIIb). A unit cell for n-hentriacontane (C31H64) was constructed from previously reported odd-number, long-chain hydrocarbons34. The Pnam polyethylene unit cell was used for indexing the diffraction patterns from 1-triacontanol wax crystals, following the established polyethylene-like packing known for long-chain hydrocarbons44. Diffraction patterns at different major orientations for each material were compared with the calibrated experimental diffraction patterns in overlay. In most cases, the samples were indexed to low-index zone axes such as [001] for p-terphenyl, [101] for anthracene and [001] for n-hentriacontane and 1-triacontanol. Mis-tilt from high-symmetry zone axes (Supplementary Fig. 12) was estimated by the inspection of simulated diffraction patterns tilted away from the high-symmetry zone axis to match the experimental data. For theophylline, auto-indexing in CrystalMaker was used to assist in finding the possible zone axes as these samples could not be indexed to low-index zone axes. All the simulated diffraction patterns at the candidate zone axes were inspected and compared with the calibrated experimental diffraction patterns. Following this approach, we indexed the theophylline crystals to the zone axes [141] for form II theophylline, [(bar{2})13] for form M theophylline and [211] for metastable form IIIb.

Automated constructions of multiple VDF images were scripted with Python using functions from Pyxem-0.11.0 and HyperSpy-1.6.5 packages. Diffraction patterns, originally 515 × 515 pixels, were cropped and then binned by a factor of two to produce diffraction patterns with dimensions of 257 × 257 pixels. Bragg discs were detected using a peak-finding function using the difference of Gaussians method. Filter settings were tuned iteratively and assessed by manual inspection to capture all disc-like diffraction features in a randomized sub-sample of diffraction patterns. After peak finding, the position of each Bragg disc in the diffraction pattern was recorded in the coordinates of the diffraction vector (kx, ky). The peak position was further refined using a local centre-of-mass search. VDF images were then constructed from the integrated intensity within a circular aperture centred at the peak position for each probe position in the SED scan. The aperture radius was adjusted to minimize Bragg-disc overlap, typically set in the range of 0.016–0.026 Å−1. The corresponding diameters (0.032–0.052 Å−1) were fixed to be slightly smaller than the diffraction-disc diameter of 2α/λ = 0.08 Å−1 (at 300 kV) to limit the intensity from any adjacent diffraction spots in large unit cells with finely spaced diffraction spots (otherwise advantageous for accessing many bend contours). For a given spot at coordinate (kx, ky), the angle of the associated diffraction vector ghkl for each Bragg disc taken as the inverse tangent of ky/kx. Dislocation lines (obscured in annular integration) were mapped by examining multiple VDF images constructed from the ghkl values at which ghkl ∙ B ≠ 0, that is, where the displacement of the bend contours is recorded. Under such conditions, the dislocation is readily traced between successive bend contour displacements across the field of view.

Bend contour displacements were measured using ImageJ software. The displacement of the bend contours on crossing the dislocation line is determined as the distance, tracing along the dislocation line, between the middle of the two bend contours (Supplementary Fig. 1 (red dashed lines) and Supplementary Note 1). The azimuthal angular variation in bend contour displacements was analysed by nonlinear least squares curve fitting, following equation (1). Curve fitting was implemented using the curve_fit function in the SciPy Python package with further evaluation of confidence intervals using the Lmfit-1.2.1 Python package45 (Supplementary Note 2).

Related Articles

Too big to purge: persistence of deleterious Mutations in Island populations of the European Barn Owl (Tyto alba)

A key aspect of assessing the risk of extinction/extirpation for a particular wild species or population is the status of inbreeding, but the origin of inbreeding and the current mutational load are also two crucial factors to consider when determining survival probability of a population. In this study, we used samples from 502 barn owls from continental and island populations across Europe, with the aim of quantifying and comparing the level of inbreeding between populations with differing demographic histories. In addition to comparing inbreeding status, we determined whether inbreeding is due to non-random mating or high co-ancestry within the population. We show that islands have higher levels of inbreeding than continental populations, and that this is mainly due to small effective population sizes rather than recent consanguineous mating. We assess the probability that a region is autozygous along the genome and show that this probability decreased as the number of genes present in that region increased. Finally, we looked for evidence of reduced selection efficiency and purging in island populations. Among island populations, we found an increase in numbers of both neutral and deleterious minor alleles, possibly as a result of drift and decreased selection efficiency but we found no evidence of purging.

Suppressed ballistic transport of dislocations at strain rates up to 109 s–1 in a stable nanocrystalline alloy

Dislocations are crucial to plastic deformation in crystals. At extreme strain rates, their motion shifts from thermally activated glide to ballistic transport, causing significant drag due to interactions with phonons, which can lead to embrittlement and failure in metals. The concept of dislons, quantized dislocations, has emerged to better understand these types of interactions. Similar to quantum treatment of dislocation-electron interactions, confining dislocations to nanometer scales, especially in nanocrystalline metals, could also yield unique mechanical behaviors different from bulk materials. Here, we present evidence showing that in Cu-3Ta, a thermo-mechanically stable nanocrystalline alloy, the phonon drag effect is entirely suppressed even at ultra-high strain rates (109 s−1). This is due to the stable confinement of dislocations within several-nanometer range, limiting their velocity and interaction with phonons. Our study indicates that in confined environments, the dislocation-phonon drag effect is minimal, potentially improving material performance under extreme conditions.

Improving the fast-charging capability of NbWO-based Li-ion batteries

The discovery of Nb-W-O materials years ago marks the milestone of charging a lithium-ion battery in minutes. Nevertheless, for many applications, charging lithium-ion battery within one minute is urgently demanded, the bottleneck of which largely lies in the lack of fundamental understanding of Li+ storage mechanisms in these materials. Herein, by visualizing Li+ intercalated into representative Nb16W5O55, we find that the fast-charging nature of such material originates from an interesting rate-dependent lattice relaxation process associated with the Jahn-Teller effect. Furthermore, in situ electron microscopy further reveals a directional, [010]-preferred Li+ transport mechanism in Nb16W5O55 crystals being the “bottleneck” toward fast charging that deprives the entry of any desolvated Li+ through the prevailing non-(010) surfaces. Hence, we propose a machine learning-assisted interface engineering strategy to swiftly collect desolvated Li+ and relocate them to (010) surfaces for their fast intercalation. As a result, a capacity of ≈ 116 mAh g−1 (68.5% of the theoretical capacity) at 80 C (45 s) is achieved when coupled with a Li negative electrode.

High strength and ductility in a rare-earth free magnesium alloy processed by rotary swaging and flash annealing

Magnesium (Mg) is the lightest structural metal. It is promising for aerospace applications if its mechanical properties can be improved at a reasonable cost, without using expensive alloying elements or time-consuming processing routes. Here, we develop a rare-earth free Mg-Al-Ca alloy, processed through rotary swaging followed by flash annealing. The resulting alloy exhibits superior strength and ductility, surpassing nearly all reported rare-earth free magnesium alloys. Nanosized Al-Ca precipitates and clusters, formed largely during rotary swaging and flash annealing, significantly strengthen the alloy. Deformation twinning is suppressed, necessitating the formation of more <c + a> dislocations to accommodate the severe plastic strain induced by rotary swaging. These <c + a> dislocations are retained during flash annealing due to the high climb energy barrier and pinning by those nanosized Al-Ca precipitates and clusters. This retention of <c + a> dislocations contributes to high ductility and provides forest hardening, further increasing strength. This study offers a simple strategy for designing and fabricating high-performance magnesium alloys.

Revealing the molecular interplay of coverage, wettability, and capacitive response at the Pt(111)-water solution interface under bias

While electrified interfaces are crucial for electrocatalysis and corrosion, their molecular morphology remains largely unknown. Through highly realistic ab initio molecular dynamics simulations of the Pt(111)-water solution interface in reducing conditions, we reveal a deep interconnection among electrode coverage, wettability, capacitive response, and catalytic activity. We identify computationally the experimentally hypothesised states for adsorbed hydrogen on Pt, HUPD and HOPD, revealing their role in governing interfacial water reorientation and hydrogen evolution. The transition between these two H states with increasing potential, induces a shift from a hydrophobic to a hydrophilic interface and correlates with a change in the primary electrode screening mechanism. This results in a slope change in differential capacitance, marking the onset of the experimentally observed peak around the potential of zero charge. Our work produces crucial insights for advancing electrocatalytic energy conversion, developing deep understanding of electrified interfaces.

Responses

Your email address will not be published. Required fields are marked *