Automated Defect and Correlation Length Analysis of Block Copolymer Thin Film Nanopatterns

Line patterns produced by lamellae- and cylinder-forming block copolymer (BCP) thin films are of widespread interest for their potential to enable nanoscale patterning over large areas. In order for such patterning methods to effectively integrate with current technologies, the resulting patterns need to have low defect densities, and be produced in a short timescale. To understand whether a given polymer or annealing method might potentially meet such challenges, it is necessary to examine the evolution of defects. Unfortunately, few tools are readily available to researchers, particularly those engaged in the synthesis and design of new polymeric systems with the potential for patterning, to measure defects in such line patterns. To this end, we present an image analysis tool, which we have developed and made available, to measure the characteristics of such patterns in an automated fashion. Additionally we apply the tool to six cylinder-forming polystyrene-block-poly(2-vinylpyridine) polymers thermally annealed to explore the relationship between the size of each polymer and measured characteristics including line period, line-width, defect density, line-edge roughness (LER), line-width roughness (LWR), and correlation length. Finally, we explore the line-edge roughness, line-width roughness, defect density, and correlation length as a function of the image area sampled to determine each in a more rigorous fashion.

section_content: The ability of block copolymers (BCPs) to self-assemble into periodic structures, with periods ranging from 5 nm to well over 100 nm, has prompted investigation into their potential applications for nanopatterning of integrated circuits, [1] [2] [3] [4] bit-patterned storage media, [4] [5] [6] [7] optical devices, [8, 9] tissue interfacing, [10] [11] [12] and others. Lamellar or cylindrical domains of block copolymers can be used to create linear structures, [13, 14] both large [8, 15] and small, [16] when confined in one dimension as thin films on substrates with appropriate wetting characteristics. [17] [18] [19] Such patterns can be used as lithographic masks through etching or as scaffolds to create other nanostructured surfaces and materials. [20, 21] For application in semiconductor fabrication, the International Technology Roadmap for Semiconductors (ITRS) has, in its Directed Self-Assembly Critical Assessment (where the term directed self-assembly is represented by the acronym DSA), identified challenges in 15 metrics, including: [22] Feature sizes of under 10 nm, the ability to "add, exclude or trim individual DSA. . . features with simple lithography", 22 a low degree of line edge roughness (LER, 3σ) < 0. 6 nm, defect density less than 0. 01 cm -2, and an annealing time of less than one minute. In addition, surfaces require appropriate wetting characteristics and surface energies in order to enable the process of self-assembly in the desired orientation with respect to the surface plane. These metrics have been correctly identified as challenges as they are daunting goals, but they represent very clear, quantified metrics that need to be attained. Lacking, however, is a unified method of accurately determining each parameter 'in the field', with actual samples of surfaces patterned via block copolymer self-assembly (vide infra). 

Much of the work to optimize BCP DSA has been carried out with a narrow range of polymers, namely polystyrene-block-poly(methyl methacrylate) (PS-b-PMMA), [1, 23, 24] polystyrene-block-polydimethylsiloxane (PS-b-PDMS), [2, [25] [26] [27] and polystyrene-block-poly (2-vinylpyridine) (PS-b-P2VP). [10, 14, 28, 29] Each of these polymers possesses favorable characteristics for nanopatterning, but many other block copolymer systems still remain to be designed, synthesized, and investigated, as the exploration of the space of possible systems, including structural classes and chemical motifs (monomers) is nowhere near complete [16, 17, [30] [31] [32] [33] [34] Dimensions of polymer-space available for exploration include triblock, comb, or other architectures [32] [33] [34] [35] and topologies; [31] alternate chemical moieties such as silicon-containing polymers other than PDMS, [36, 37] and oxygen-rich groups such as oligosaccharides and poly(lactic acid); [16] or tailoring polydispersity to modify morphological stability and domain sizes. [38] [39] [40] A consequence is that there remains much to be explored synthetically in order to optimize pattern formation, etch selectivity (or resistivity), polymer reactivity, surface energies, Flory-Huggins parameters, [41] LER, and annealing conditions. In particular, there is a persistent analytical barrier that synthetic chemists must overcome in order to readily determine whether their polymeric creations may be applicable to novel DSA applications: they require access to a toolbox capable of analyzing critical features such as the defect density, correlation lengths, and LER of their patterns in order to determine whether their block copolymers have promise. The dearth of accessible tools remains a significant obstacle for the area of directed self-assembly. 

Defects themselves also warrant a more in-depth investigation, which can only be achieved by studying defects "in the wild", in the actual nanopatterns as they progress through various stages of annealing. While simulations can find matches to thin film defect structures, [42] [43] [44] automated analyses of defects in block copolymer thin films in an experimental setting allow access to statistical data about the frequency and distribution of various defects. Statistical data is generally inaccessible via modeling due to computational limits for defects beyond the simplest examples. [42] Furthermore, for cylindrical block copolymer domains, defects do not always have liquid crystal analogues, [43] rendering past defect-detection methods, originally developed for patterns formed in liquid crystal thin films, [45] inappropriate. Hence identification of structures beyond simple counting of disclinations and dislocations [46] would be advantageous. 

Computerized analyses that are widely used to study images of BCP patterns include to determination of periodicity using azimuthally averaged fast Fourier transform (FFT) images, [47] image filtering and regional analysis of domain orientation using FFTs [48] and defect density measurements. Less commonly measured is LER. [49] Rarely, however, are more than one or two methods packaged together in a published work, which leaves unanswered questions since there is typically a trade-off between factors such as defect density, orientation, nanostructure spacing, line-width, line-width roughness (LWR), LER, and correlation length (i. e., grain size). Initially, in the course of investigating the process of microwave annealing of block copolymer thin films, we developed an in-house algorithm to quantify defects in block copolymer thin films, utilizing particle analysis and skeletonization to identify defect features. [28, 50, 51] Later, while analyzing density doubled cylindrical line patterns, a separate process for measuring the LER was created, which was limited to analyzing nearly-straight segments of nanowire structures. [29] To remedy the lack of a readily available and straightforward analytical tool, we developed an accessible and free-to-download application for analyzing the defects in BCP thin films using a combined particle and skeleton based analysis of the pattern, called Automated Defect Analysis of Block Copolymers, or ADAblock for short (links to the tool provided in S1 Instructions). The application was constructed using the ImageJ platform, a free, open-source, Javabased image analysis program, which provides a full and easy-to-understand output. [52, 53] The tool identifies not only the type(s) of defects found in a sample, but also quantifies the density of defects over a range of length scales, accompanied by additional information regarding LER and LWR, as well as an alternate means of accessing the correlation length. In this work, we screened the ADAblock application against a range of nanopatterns prepared via block copolymer self-assembly and show the effects of polymer molecular weight on the defect densities of self-assembled BCP films. Additionally, we demonstrate how ADAblock can simultaneously track LER, defects, and correlation lengths. To our knowledge, no previous work analyzing 2D block copolymer line patterns has brought together data on defects, LER, LWR, and correlation lengths into one application or analysis. We believe that this omission is likely due, in part, to the lack of readily available, widely applicable, easy-to-use tools for analysis, and on occasion, may be a result of selection of the 'makes-it-look-best metric', rather than a complete description of pattern quality over larger areas of the sample. In this paper, we show how such data can ideally be combined to better describe line patterns derived from BCP assembly in thin films. Images can be deceiving, and we hope that ADAblock will assist researchers in avoiding pitfalls resulting from performing incomplete defect density analyses.
section_name: Introduction

section_content: Examples of four different nanopatterns derived from BCP self-assembled templates are shown in Fig 1. The patterns have been converted into easily visible platinum lines through a welldescribed platinization of three different PS-b-P2VP BCPs; [14, 28, 29] the Pt nanolines are derived from the P2VP blocks. Although these scanning electron micrographs (SEMs) are similar in appearance, each is subtly different, and thus the question to be posed is how to distinguish one pattern from another and to determine which is more defective. As shown in Fig 1 and Table 1, the pattern in 1A has 30% more defect pairs than the patterns in Fig 1C or 1D. Moreover, the correlation length of 1B is shorter than any of the others (in part due to the shorter period). Additionally, in terms of the line edge roughness (LER), they all appear at first glance to be quite smooth, but the measured roughness of these lines, as summarized in Table 1, would put them out of contention for ITRS targets. With respect to LER, the values of ~4 nm are significantly larger than the maximum 0. 6 nm suggested, [22] but the feature size here is also ~2x larger than the 10 nm features sought by the ITRS; similarly, the defect density is ~10,000 times higher than ITRS goals. However the present samples lack any features to guide alignment, as is the case with graphoepitaxy, which assists in significantly lowering the observed defect density. [28, 44, 54]
section_name: Results and Discussion

section_content: The analysis is briefly outlined in Fig 2, breaking down the ADAblock sequence into eight broad stages. Details of each stage are provided in detail, vide infra. The first stage is simply the representation of the original SEM image, which must be smoothed to reduce noise, while retaining all features of interest. Next, the smoothed image is thresholded, in order to produce a binary image from which data like area, perimeter, and shape can be determined. Next, the period and line-widths are calculated, followed by particle analysis to determine the shapes of the binary objects and to classify line and dot features. The line features identified are then isolated and converted into a skeleton, from which the connectivity can be determined. This resulting structure is groomed and then analyzed for defects. Lastly, the data is recorded and confirmation images are produced for user inspection. All stages noted in the text correspond to the stages represented diagrammatically in Fig 2. Units are in nm, except for defect density, which is given as defect pairs per μm 2. Line-edge roughness (LER) is given as three times the standard deviation (3σ) in the edge position, relative to the center of the line; line-width roughness (LWR) is three times the standard deviation (3σ) in the width of the line. 

doi:10. 1371/journal. pone. 0133088. t001

Stages 1 and 2: Input of original image and smoothing. In order to extract information regarding defect density and LER (above, Table 1 ), a number of factors must be taken into account. Starting with the image itself, basic parameters must be adapted to (1) the image resolution [for instance, determination of how many nanometres are represented by each pixel (nm/px)]; (2) the contrast of the image, which can vary considerably image-to-image and instrument-to-instrument; (3) image noise, which creates artifacts not inherent to the actual structure under investigation; and (4) the period and line-width of the block copolymer. Given the nature of block copolymer patterns typically observed, certain presumptions about the structures observed within the images can be made. To begin, predominant structures within a given image are primarily limited to dots, lines, and meshes. Classification into the basic families of structures in turn constrains certain shape characteristics for the features. Additionally, the period can be defined for a relatively narrow range (e. g. 20 to 40 nm for the samples described here, but modifiable for a given system), as block copolymer samples for a given image data set can be manually selected to those having similar period values. As preliminary background data, the period of a pattern can be obtained via azimuthal averaging of the image's fast Fourier transform before application of ADAblock. [47] The first item, the image resolution, is frequently embedded within the image's metadata and hence can be called by the program or input by the operator. The preponderance of our BCP pattern images were obtained using a Hitachi S4800 scanning electron microscope (SEM), which provided information in a legend at the bottom of the image; the consistency of this feature also provided a means for automated extraction of resolution parameters. 

Combined with its high resolution and high throughput, SEM can be the ideal imaging tool for BCPs, although it does have some drawbacks: For our work, smoothing was necessary due to random noise, charging effects, and edge effects. In the case of SEM images, edges can possess enhanced brightness, [14, 28, 29, 50] and white noise results in speckling of the image with bright and dark pixels. Without some smoothing, such salt-and-pepper noise can result in unwanted extra features. SEM owes much of its brilliance to edge effects, which result in objects protruding from the surface (such as Pt nanowires on Si) appearing much brighter than surrounding substrate. [14, 28, 29, 50] Typically, smoothing images involves trial-and-error, but linking the smoothing to the period of the pattern and the image resolution gives consistent results: Gaussian and/or median filtering are automatically applied with filter radii calculated in proportion to the period of the pattern to avoid under-and over-smoothing. Median filtering typically is best, as it can preserve and even enhance the structure of the line pattern, as shown in S1 Fig

Stage 3 : Threshold to binary. In order to analyze the pattern, the two phases, each corresponding to one of the blocks, must be clearly identified and separated by thresholding. This assigns each feature in the image to one of the two phases, referred to herein as positive (i. e. bright) and negative (dark). Contrast enhancement and thresholding typically requires manual intervention as well. For SEM and atomic force microscopy (AFM) images, which are typically used for block copolymer thin films and patterns, as well as helium-ion microscopy images, a bimodal histogram is either typical or attainable given the nature of the pattern. Such a bimodal histogram can occur either globally (i. e., over the whole image) or locally (over smaller subregions); a suitable thresholding filter can be applied on either scale. Several "auto-local" thresholding plugins are available for ImageJ; analysis of our images typically works best utilizing an auto-local threshold which applies Otsu's clustering method [55] locally across the image, however, other thresholds implemented in ImageJ are available options. [56] When the surface is not uniformly covered by features (e. g. featureless regions), however, automated thresholding can result in additional artifacts, hence subsequent steps are taken to remove noise and incorrectly phased features. 

Stage 4: Initial Line-Width Analysis. The first pieces of data that must be determined are the dimensions. While period can be readily and automatically measured from azimuthally averaged FFT patterns, [47] line-widths and spacings cannot be derived directly from the image of the BCP nanopattern. Profile plots can make for easy manual measurement of these features when patterns are regular and aligned, but that is not always the case. Knowledge of the dimensions is useful, even necessary, in contexts where the pattern is poorly ordered. Particle analysis can measure the area and perimeter of each particle accurately. While Feret measurements (See S2 Fig) work for simple particles, the tortuous nature of BCP "lines" calls for more nuanced measurement. Imagine a spaghetti noodle shape confined in 2D; there exists a relationship between the perimeter of the noodle's edge and the area covered by the noodle. Using these easily measurable geometric quantities-particle area and perimeter-the width of lines can be calculated, straightness and degree of branching notwithstanding. Provided there are enough lines available, the particle area plotted as a function of perimeter is linear as shown in Fig 3; the slope of the plot is half of the width of the line. 

For lines without junctions and only uniform tips, perimeter, P, can be broken into

where P t is the perimeter of each tip region (see Fig 3H ) and L is the length of the main portion of the line. The area can be calculated similarly

where w is the width and A t is the area of the line's tip region. Area as a function of perimeter can be calculated by substituting

into the area equation, giving

As previously stated, the slope is 0. 5w. For patterns where junctions rather than terminal points are predominant, P t becomes zero, and the intercept is positive (due to an additional area term derived from the junction). In practice the contribution of the ends (thus intercepts) is negligible relative to the segment lengths (L >> w) for images with junctions. With exclusively semicircular terminal points or triangular junctions, one would expect intercepts of -0. 785w 2 and 0. 289w 2 respectively. (See SI for calculation. )

Intermediate combinations can be avoided by temporarily excluding junctions and breaking down the binary pattern into smaller, junctionless particles, as shown in Fig 3B. By excluding the junctions, along with excessively small particles and sections of particles on the edge, a better fit can be obtained, providing a better estimate of the line-width. (See S1 Appendix for a schematic depiction of how junction exclusion achieves this. ) The value of the intercept in Fig 3A is given as approximately -200 nm 2, which is reasonably close to the predicted value of -250 nm 2, detailed in S1 Appendix. 

Repeating the process for the negative phase, separately, gives a measure of the spacing between lines. Summing the two measurements to approximate the period has, in most cases, been found to come within 5% of the period measured by FFT, usually slightly greater. It appears that this discrepancy may be due, in part, to the particle perimeters being larger than non-discrete analogues, and also due to the approximations made herein. Alternatively, knowledge of the line-width and the period would give the line spacing by difference. Stage 5: Particle analysis. In the course of annealing spin-coated BCP thin films, the pattern may evolve from a dot pattern, or similarly disconnected collection of features, into an array of lines, with numerous defect-rich intermediate states. The binary image can be analyzed using ImageJ's built-in particle analysis routine to determine characteristics of particle size and shape descriptors such as circularity and Feret measurements, mean pixel values, and relationships to the image boundary for each feature. Particle analysis is done separately for the positive and negative phases in order to access all the features. Particle analysis data is then used to separate dots (or other objects), which cannot be accurately treated as lines, and identifies them as a specific type of defect. It can also provide information on the evolution of particles in the course of the annealing process (e. g. increases in the average size or length of lines). Moreover, the creation of a binary pattern further enables distinction between noise and misclassified particles. 

Stage 6: Skeletonization. The most effective means to analyze the topology of lines and meshes is via analysing the connectivity of the pattern by creating a skeleton of it. Skeletonization reduces lines or meshes to binary objects which maintains the connectivity of the original by "thinning" the pattern to create a simplified, single-pixel-wide version of the shape, suitable for pixel-by-pixel analysis. Skeletonization algorithms and skeleton analysis has been widely used in other fields to study the topology of structures, from text recognition algorithms in computer science to numerous subfields in biology, including bone analysis ("bonej"), [57] and studying the structure of neurons, as well as for the recognition of typographic characters. In all cases, skeletonization is used to simplify collections of interconnected shapes and objects into networks to study their properties. Although other algorithms do exist, the default technique is implemented in ImageJ: [58] the skeletonize function in ImageJ uses a lookup table to progressively thin the structure based on each pixel's 3x3 neighbourhood, leaving a 1-pixel wide topological skeleton. [58] At least two other groups have applied skeletonization as a means to interpret BCP thin film patterns. [59, 60] Rehse and coworkers utilized skeletonization of one phase of the polymer pattern to study frame-to-frame correlations between junctions as a measure of BCP dynamics; [59] their work followed that of Scherdel [60] and Vigild [61] who used 3D interpenetrating skeletons to describe gyroidal phases. 

The skeletonization process itself is quite straightforward, as shown in In order to actually access the defects in striped BCP patterns, it is necessary to investigate each phase separately, thus requiring parallel particle analysis and skeletonization of each phase. It is worth noting that for BCPs, certain kinds of defects will prefer one phase to the other, resulting in a surplus of terminal points or junctions in either phase. This tends to limit the frequency of spatially paired defects. (a) Locating defects. The ideal, defect-free line pattern derived from lamellar or cylindrical domains of block copolymers, consists of perfectly parallel straight lines extending across the entire substrate without interruption, as by breaks or junctions in the lines. It is with respect to this ideal that topological defects are defined. The analogy between block copolymers and liquid crystals (nematic and lyotropic phases in particular) inspired previous defect analyses [63] utilizing winding numbers to identify and measure topological defects. While this does work in principle, and many previous analyses have utilized it and other defect-detection methods, [46, 59] these methods typically are published without full working details or code. Kléman suggested in 1983 that defects in two-dimensional line patterns could be simplified to junctions, terminal points, and dots, shown in Fig 6, rather than the more conventional approach using winding numbers to determine the type of defect. [64] While such methods correctly describe the type of defect, [46] high levels of defects and variability of the patterns, including variations in line-width (or LER) can make it difficult to correctly identify and quantify defects. One author noted an order of magnitude change in the density of defects for one particular image depending on the amount of Gaussian filtering applied; the filtering parameters ultimately selected appeared to be arbitrary. [45] Furthermore, images depicting disclinations and dislocations suggest imprecision in the identification of closely associated dislocations. [46] Moreover, while such methods provide the magnitude of each defect, they do not provide information about the connectivity of defects or orientation of surrounding features. 

Many frequently encountered defect structures are not isolated dislocations (junction-terminal point pairs immediately adjacent) or disclinations (either terminal points or junctions), but are part of more complex defect structures. These structures can be broken down into component dots, junctions, and terminal points, the elementary components of defects shown in Fig 6. Dots can be determined best using the particle analysis data, so from the skeleton analysis we locate and characterize the junctions and terminal points. 

At this point in the analysis, the skeleton is a binary object where lines are represented by a series of 2-connected pixels in any of the 8 directions; terminal points are singly-connected pixels; and junctions occur where more than 2 pixels are neighbouring a given pixel. Dots do not, however, always reduce to single-pixel objects and hence they are treated separately. Furthermore, junctions exist in numerous possible configurations, often with multiple (3+)-connected pixels per junction, hence there will not be a one-to-one correspondence between junction pixels and either the number of junctions or junction types (see Fig 7). Identification of defects is done in a manner, which is, in essence, analogous to playing the minesweeper-type [65] games: by counting the number of skeleton pixels adjacent to any given pixel that is part of the skeleton, thus providing a connectivity value for that pixel, as shown in Fig 7. 

At its simplest, any connection or disconnection that breaks the 2-connected topology of the skeleton, resulting in a new local topology (or connectedness), is a defect with a corresponding value. For junction points (JP), with each additional branch beyond two (which, on its own, would constitute a line without topological defect), the defect increases in magnitude by ½:

where B is the number of branches. Typically junctions come with only 3 branches, but 4-way, or even 5-way, intersections can be found, on occasion, between clusters of dots or other complex features. A 4-way intersection would be n jp = -1, which can be imagined as being derived from two adjacent junctions, each with n jp = -½, with a common line-segment, where the intervening line segment's length decreases to zero; the same approach can be generalized for any number of additional line segments (as shown in S3 Fig). 

Terminal points (TP) possess only one configuration, hence their value is assigned:

Dots can be considered as a line with two terminal points, collapsed to a single point (as depicted in S3 Fig), hence their value is twice that of a terminal point:

Other more complex structures, such as spirals (containing terminal point dislocations) can be counted via their component structures in this regard. Large, solid spots in the bright phase or large regions without any pattern (i. e., large spots in the dark phase) possibly formed due wetting (or other causes) may exist, these regions may be treated as dots with radiating arms, however, we found it was more effective to separate the core of the dot prior to skeletonization. The result treats the dot as a kind of enlarged junction, with defects existing only at the periphery. 

For all defects to be counted, skeletons for both phases must be generated and connectivity analyzed separately. This raises an important point for defect analysis: that defects exist in a particular phase. The phase dependency of defects has not been sufficiently explored, however hints are seen in the literature, as it directly controls the topology of a system. [66] The dual phase analysis brings about an addition rule for description of the system: on average, the sum of all defects in the pattern should be zero. Alternatively, this can be stated that every defect is "paired", hence for every junction there is a terminal point; Moreover, for every dot, there will be two junctions: consider this an analogue of a unit cell (see S7 Fig), as two defect pairs are produced if spontaneously generated and two are required to cancel out through annihilation. [46] Typically there exists a small imbalance between the two measures, which others have observed as well. [45] The pairing of defects does not imply however that the number of defects in the two phases of the block copolymer will be equal. 

Skeletons also provide a description of the connectivity of defects, which merits further exploration. Defects can also be associated with particles using this method, but perhaps the greatest benefit is derived from the ability to search for and positively identify particular clusters of defects. One such example is an H-junction, which results from a break in the line or a bridging of two adjacent lines, shown in S4 Fig These junctions are supposedly not the result of a defect in the actual thin film structure, but result from (a) incomplete metallization or other means of pattern transfer, (b) image noise, or (c) the smoothing-thresholding process. Hence it may be prudent to recognize them and count them separately or to "correct" such errors in the binary image itself. 

(b) Grooming the Skeleton. Grooming the skeleton consists of trimming away short branches, which may result as artifacts from small "bumps" on the edge of a line. With the dimensions determined in Stage 5, we can create a metric to selectively prune away any branches resulting from variations in line-width or simply from sharp points or edge effects that can influence the skeletonization algorithm. It may be a point for philosophical debate what constitutes a branch, justifying a junction and terminal point, but objectivity can be introduced by basing the grooming procedure on the measured LER. For this purpose, any end point separated from a junction by less than 1. 5 × line-width (for a given phase) is considered roughness, rather than an additional defect pair, and is hence pruned, as shown by the example in Fig 8. 

Because any image represents a finite sample of a larger structure, defects at image edges must be carefully treated. Depending on the resolution of the image and the domain size of the block copolymer, these can for smaller images, represent a significant fraction of defects; additionally, in otherwise low-defect patterns, features cut off at the edge may appear as additional defects. In particular, three rules must be applied:

1. Any "dot" (or sufficiently small object without junctions) touching 2 edges is not a defect. 

(See S4 Fig) 2. Lines that run roughly parallel to the edge, touching at all times, are not defects. 

3. Lines that terminate at the edge of an image are not defects, as it is not a true terminal point. 

The third rule requires some manipulation of skeleton points & component terminal points near edges, as ImageJ's native skeletonize algorithm can produce limited edge artifacts. What these rules do not address is particulate matter; other analysis methods tend to default to manual identification or require equipment unavailable to most researchers. [54] (c) Line-Edge Roughness and Line-Width Roughness. One of the chief questions posed for BCP lithography is whether lines can be produced with sufficient uniformity and with smooth edges. LER measures the variation in the position of the edge of a line, which can have different frequency components, leading to undulation of the edge and variation in the width of the line, or LWR (Fig 9 ). The variation in position is measured as the standard deviation in the position of the edge, and LER is reported as 3σ. Such variations are deleterious for circuit elements: For transistor gate features with widths < 85 nm, line roughness causes significant variations in the off-current, as well as affecting threshold voltages. [67, 68] For nanometre-scale interconnects, line roughness increases both resistance and capacitance, [69, 70] resulting in degraded transistor performance. 

The line roughness of block copolymer nanostructures has been considered theoretically and has been shown to depend on χN, [71] [72] [73] and polymer polydispersity; [73] results have suggested that the Flory-Huggins χ parameter may need to be increased by a factor of 3 to 4, relative to that of PS-b-PMMA, [72] in order to decrease LER sufficiently to accommodate ITRS targets. [74] It has been specifically noted that there are few reports on the topic of LER/ LWR in the literature; [49] typically, the actual position of the edge is measured relative to the ideal or average edge position for straight or aligned lithographic patterns. In order to achieve the same measurements for block copolymers, films aligned via graphoepitaxy would typically be required in order to have linear lines representing ideal edges. However, we [29, 75] and others [76] have taken the approach of measuring LER for unaligned patterns. One may measure edge positions relative to the centre of the line, rather than with respect to a linear ideal edge position; the standard deviation in the edge position will be the same either way. As lines get narrower, however, the influence of pixel position can begin to slightly increase the measured LER, up to 0. 5 nm in our previous work using high resolution (ca. 100,000x) BCP patterns. We mitigate this, in part, by smoothing both the centre line of the skeleton and the outer edge, while constraining the positions of the edge points. Edge-to-skeleton distances are determined for all points on the smoothed line edge, matching with the nearest points (shown in Fig 10A ) on the smoothed skeleton line which satisfy:

As derived from the dot product of the vector on the edge-to-skeleton distance and the orthogonal vector (1, slope) of the skeleton at that point, an interpolated point on the skeleton can be obtained (shown in Fig 10B). 

Line-width measurements can be made in conjunction with edge-to-skeleton measurements by finding a line segment on the opposing edge, which is intersected by the vector made between the edge point and skeleton point of the previous step (shown in Fig 10C ). The solution exists at a point on the line segment formed by the vector between the edge (x edge, y edge ) and the skeleton (x skel, y skel ) is scaled by a factor, a, and on the line segment formed by the vector between two consecutive points on the transverse edge (x trans1, y trans1 ) & (x trans2, y trans2 ), scaled by a factor, b (shown in Fig 10D ). Provided that the two vectors are not parallel, the equations [77] for the scalars, a and b, are:

An intersection is considered valid when 1 < a < 4, indicating that the side opposite would have a width ranging from 0 to 3 times the width of the first side. The limit, a < 4, prevents identification of points on parallel segments, as with a hairpin, from being identified as valid; typically the period is on the order of 2 times the width of a given line, hence 4 times the half-width of a line. In practice, the values of a are in the range 1. 5 < a < 2. 5, as can be seen typified in Fig 11 via the histograms. The second limit for valid points is that 0 b ! 1, which ensures that the point of intersection is within the line segment formed by the two consecutive edge points. 

In order to obtain reasonable measurements of LER and LWR, the blocky structures of binary lines and skeletons need to be smoothed. The smoothing process, which we have utilized here, involves 4 stages: A. Centring of the skeleton points by adding 0. 5 px to each x and y coordinate. This accounts for the slight truncation from the skeletonization process and makes the edge-to-skeleton distances more equidistant on each side. 

B. Shifting all edge points to the midpoints between consecutive points. This averaging reduces roughness introduced by the shape of individual pixels. 

C. Smoothing the skeleton by iteratively averaging the positions of points, while limiting the displacement to within 0. 25 pixels. This provides a smooth, continuous, reasonably centred skeleton line. 

D. Smoothing the edges likewise provides a smooth edge while maintaining the shape and deviations in width, from which roughness can be measured. 

Fig 11 shows the data for a single line as it is modified by each of these four smoothing processes (A, B, C, D). By the fourth stage (D), the data shows considerably less noise. In particular, the histograms of edge widths, depicting which edge is further from the skeleton, for each point on the edge, begins to approach a normal distribution, as one would expect for a line with random variations in width. Visually, the line becomes sufficiently smooth that pixels are no longer apparent, while variations in width are in keeping with the original image, and the sequential widths and edge positions measured do not have large point-to-point changes in displacement. While the skeletonization algorithm is largely effective in finding the centre line, it is imperfect. In particular for lines with pixelated widths less than 7 pixels, the centre will tend to be skewed preferentially depending on the orientation of the line. However this does not affect LWR measurements and smoothing does help to limit the impact on LER. 

In order for BCPs to be relevant in industrial manufacturing, they must achieve a low frequency LWR (3σ) of 1. 1 nm on features 16 nm wide; in order to "significantly exceed" conventional lithography, the patterns would need to be better than 0. 6 nm LWR on features 9 nm wide. [74] Presently our best measured samples have a LER (3σ) of ~2 to 3 nm, [29] however, no aspect of the process has, as of yet, been explored with respect to minimizing LER or LWR. To avoid the local effects of junctions and to increase the speed of the calculation, the lines are modified, as shown in Fig 3B, to render all lines junction-free. Additionally, points where lines contact image edges are selectively modified, erasing large contacts, to prevent any effects of the image edge. 

(d) Correlation Lengths & Order Parameters. Correlation lengths (or orientational persistence lengths) are typically calculated for large images, often with low resolution (pixels/nm), by subdividing the area into overlapping squares, for which azimuthal angles are derived from two-dimensional FFTs of each region. [78] Lack of clarity for such images sometimes necessitates filtering in order to avoid disordered regions. In this work we implemented an alternative means of determining the 2D correlation function using the skeletonized lines. Skeletons are groomed to remove junctions and loops are broken to provide isolated lines. Orientation along the skeletonized lines can be calculated using a rolling average of each line's tangent to provide smoothly varying angles along the lines. In a typical image, there can be over 20000 points in the lines; calculating the correlation length using every point is feasible, however for expediency, the set of points can be downsampled or randomly sampled to a smaller set of ~4000 points, which provides faster calculation with minimal trade-off in terms of accuracy. From the set of orientation angles, φ(r), the correlation function, C(r-r'), can be calculated. 

Advantages of this method include ease of applicability to higher resolution, smaller-area images and images with disordered regions where, due to defects, line segments are particularly short, and φ(r) might not be determinable via FFT. This is demonstrated in Fig 11. 

The correlation function is fit using an exponential function,

where κ is the correlation length, a characteristic measure of the degree of ordering in the film, which describes the average distance over which orientational order is preserved. The correlation length should be proportionate to the grain size, as illustrated by the circles in Fig 11, which are approximately keeping in proportion with the domains visible in the orientationallycolour-mapped pattern image. However the circles are unquestionably smaller than the observed domains. One disadvantage of this method of determining κ via skeletonization is that one observes a periodic variation (corresponding to the periodicity of the pattern) in the correlation function, as shown in Fig 12. This periodic variation is a result of features separated by non-integer line spacings tending toward greater disorder than points separated by integer spacings. This appears to be due to influence by neighbouring defects. The large undulation in the curve can be partially compensated by using both the positive phase and negative phase skeletons (thus reducing the period and amplitude of the variation, however exclusion of non-line areas may be necessary), by binning measurements, as is typically done in FFT-based methods, [79, 80] or by smoothing, as we apply in the algorithm. Herman's orientational parameter, S, [81] gives a measure of how uniformly oriented the lines within an image frame are. It can also be readily calculated using the set of orientational data:

The reference angle can be set as the average orientation for the whole image, thus giving the best orientation parameter for a disordered image. Because it is widely used, we implemented this calculation into our code, however, Herman's orientation parameter tends to be less useful than the correlation length, as it can be significantly influenced by the size of the area sampled. That is to say one can typically choose a sample area small enough to give S 2D ffi 1 (perfect net order) or an area large enough to give S 2D ffi 0 (no net order). The code may, however, be adapted to set an angle where a particular direction is induced via processes such as directional annealing [81] or graphoepitaxy; in such cases, S 2D = -0. 5 is a possibility for samples where the line orientation is orthogonal to the desired orientation. [2] Finally, this skeleton-based approach facilitates generation of pseudo-coloured orientation maps, as in Fig 12, which also avoid grain-edge averaging problems exhibited with other methods. [82] Such images may assist researchers in qualitatively grasping the orientational ordering in their system. Such visual checks, can provide researchers with an accessible means of confirming numeric results, as it allows for a qualitative, direct measure of grain size on the image. 

Stage 8: Output and confirmation images. Finally, as a result of these considerations, we seek to provide self-assurance and quality control by creating confirmation images, wherein features described numerically are mapped onto real images to provide visual feedback of the accuracy of the measurement, as shown in Fig 13, which shows the defects found alongside the associated SEM images. This step is ultimately the means to determine whether the defects identified are (1) a true representation of the pattern and ( 2 ) are in the correct location. Such images of pattern orientation, line roughness, defects, and thresholding provide visual confirmation that all stages of the analysis proceeded correctly. Specifically, one can check simultaneously whether the thresholding, connectivity, grooming, and defect identification have all functioned as expected. 

Such visual feedback also lets researchers, particularly those presently involved in synthetic work, to tangibly grasp the important aspects of the pattern quality. By encoding the information spatially with colours and shapes rather than relying purely on the abstraction of defect densities and correlation lengths, ADAblock's visual feedback can function as a guiding indicator for selection of optimum structures and conditions. The data output, both numerical and visual, make it possible to engage in exploratory data analysis [83] to discover new trends, motifs, and outliers in the data available, as demonstrated in
section_name: Outline of the analysis

section_content: In order to demonstrate the utility and versatility of this application, two different scenarios and questions are posed. First, what is the effect of image resolution, and the area sampled, on the measured defect densities, LER, and other parameters for patterns derived from self-assembled BCP thin films? Secondly, what can we learn from investigating the data provided by these samples, by examining the relationships between different features, to identify features that warrant further investigation-and what does this suggest about the resulting properties of a self-assembled BCP film?

Effect of Resolution and Sampling Area. When measuring defect densities, correlation lengths, LER, and LWR, the area sampled and the resolution can potentially affect the measured results. Ideally, for any measurement, the effect of sample size must be analyzed and understood in order to obtain reliable results. To develop a general sense of how this and different polymer sizes are affected in the analysis, we annealed 5 different polymer types, each with approximately ideal thicknesses, for 20 minutes at 200°C and imaged the resulting metallized patterns at different magnifications. 

The effect of resolution in the LWR measurements in Fig 15A appears to be minimal, although there is a slight downward trend with increasing resolution (smaller image area) for the two smallest polymers, where the LWR (1σ) values decrease from 0. 18 to 0. 14. The increase in LWR is primarily observed for those samples with the smallest period, which would likely be on account of pixelation of the lines, as suggested by S8 Fig. A confounding effect may also result from the decreased length of line sampled for images of higher resolution. LER data, on the other hand, shows a more consistent trend of decreasing LER with increasing resolution in Fig 15B. LER is likely more affected by pixelation due to the inability of the skeletonization process to precisely locate the line center, in particular when line-widths are a small, even number of pixels. In contrast, the line-width is not strongly constrained by the determination of the line centre. The magnitude of the decrease (-0. 03 to -0. 05 pixels) here is still small, given that image area changes by a factor of up to 100. 

Sampling effects can be observed in the measurement of defect pair density at various resolutions in affected, due to having a larger grain size. The same effect can be observed for correlation length measurements in Fig 15D, although this affects all of the polymers. In order for the correlation length measurement to be meaningful, the measured value should be significantly shorter than the dimensions of the image. The decrease in the average measured correlation length as a function of the image area suggests that one may be able to estimate the true value based on the size of the image. The plot of 23. 6k-b-10. 4k is particularly telling because it shows the effect of sampling within a single grain or few grains (at low resolution) and the sudden decrease once more grains become involved. The limitation of large grains may be partly avoided by using automated data collection, combined image stitching, which has been demonstrated to be effective for imaging large areas with electron microscopy, [84, 85] however as ordering approaches perfection, grain sizes become infinite, [86] and the correlation function will approach unity. 

Feature Relationships. In order to derive lessons from the data, we undertake a form of exploratory data analysis to chart the relationships of different parameters observed. In particular, whether parameters such as LER and LWR are independent of the feature size, and how line-widths, polymers, and periods have a simple relationship. 

Taking all of the data (across resolutions) for each polymer, we note that as a proportion of the line-widths, the standard deviations in the edge position (LER, 1σ) and line-width (LWR, 1σ) stay constant, about 10% and 16% respectively, indicating that the LER and LWR scale with the line-width dimension of the polymer, as shown in Fig 16A and 16B. The set point may be a property of a given BCP's Flory-Huggins parameter, indicating a higher χ required. However we must caution that other factors, such as the processing, metallization, plasma treatment, and lack of alignment are convoluted with the roughness inherent to the polymer, preventing a direct conclusion. However this method should enable comparison between polymer templates and patterns translated from the BCP via etching or other means. The values observed here would however all exceed LWR targets set by the ITRS for LWR (3σ) of less than 6%: 1. 1 nm for patterns with 18 nm feature size; or <0. 6 nm for patterns with 10 nm feature size. [22] For aligned patterns, solvent annealed with water as a co-solvent, we have observed significantly better LER and LWR values. [29] We hypothesize that it may be the result of the water selectively partitioning inside of the P2VP block during annealing, resulting in a higher effective χ, leading to a smoother interface than we attain here with thermal annealing. 

Line-width in Fig 16C and 16D shows the expected relationship of being proportionate to the period, although there does appear to be a greater spread in the width of lines than in the FFT-measured periods. This is likely an effect of thresholding, which needs to be done relative to each image. It may be possible for a specific polymer or a series of images to constrain the threshold, as a fraction of area, in order to obtain a narrower distribution of line-widths. 

Limitations of the code. As with any programmed analysis, there are drawbacks and tradeoffs made in analysis to optimize for speed or accuracy. The approximations we implemented are one reason that necessitates a full sharing of the code. ImageJ's macro language is interpreted, hence it is slower in processing compared to plugins or other compiled programs. It is, however, easily edited and modified, which enables adaptation where modification may be required. The code was written so that it can be operated in a batch mode to process a folder of images, meaning that a series of images can be processed overnight, or while attending to other tasks. It should be cautioned that in the present state, as ADAblock continues to be developed, the code may produce a reproducible error for ~4% of images at present. Further refinement should reduce this error rate, but at present may limit a series from being completed. With manual intervention, however, the image can be skipped, or the settings modified, and the queue re-continued. 

Typically an image with dimensions of 1280 x 896 pixels (the default of our SEM, for example) requires ~7 minutes to process when run on the standard personal computers that we used for testing. Higher pixel-resolution (e. g. 2560 x 1792) images require more time to process, roughly in proportion to the number of pixels. Given the automated nature of the program, it's possible to run a queue of images overnight, rendering the increased processing time irrelevant. In addition to images showing the locations of defects, the code saves several check images to act as references to help determine whether any errors have taken place or other undesirable operations. Consequently, ~16 MB is recorded to the disk for each image processed, as presently conFigd, although non-graphic data only accounts for less than 300 kB. (See S1 Instructions for a list of files output by the program. )
section_name: Application of ADAblock

section_content: We have developed a facile, automated, and reliable analysis for striped patterns derived from the self-assembly of BCP thin films, that integrates both conventional and newly developed techniques. This analysis is done in order to quantify defects and their types using a skeletonization-based method; to measure line-edge roughness; and to calculate Herman's order parameter and the correlation length in a novel fashion, based upon the skeletonized structure. Moreover, the skeletonized structure provides information about the connectivity of patterns. We expect that this will be of use to others carrying out annealing studies and preliminary characterizations of novel self-assembling polymeric materials. Finally, for 5 block copolymers of similar composition, we have found the metallized patterns to have LER and LWR in roughly constant proportion to the line-width. 

Ultimately, no one measurement provides a "complete description" of pattern quality; typically they are complementary. Hence this work represents an attempt to broaden the scope of analysis and to make tools which may not be readily accessible to all. Additionally having shared protocols, or at least protocols derived from a common origin, we might be able to standardize a broad toolset, providing consistent analysis via fully shared code. [87, 88] We hope this aids comparisons between polymers, between papers, and between scientists seeking to understand the characteristics of block copolymers, and in addressing the numerous critical issues associated with block copolymer lithography. [89]
section_name: Conclusions

section_name: Materials and Methods

section_content: 100 mm diameter, single-side polished silicon wafers were diced into squares with dimensions 1. 0 cm x 1. 0 cm. Prior to cleaning, substrates were scribed, on the unpolished side, with a diamond-tip, to mark the identity of each substrate as part of a set of 10. The samples were then immersed in methanol and sonicated for 15 minutes in glass beakers. Next, after rinsing each substrate square in a series of beakers filled with 18. 2 MOÁcm water, the substrates were placed polished-side-up in PTFE beakers, and immersed in 6. 0 mL of concentrated H 2 SO 4, to which was added 2. 0 mL of 30% H 2 O 2, before placing the beaker to stand in an 80°C hot water bath for 20 minutes. The piranha solution was then decanted to a glass flask to cool prior to neutralization. 

Following several rinses with water, the substrates were immersed in aqueous 1% NH 4 OH solution for 5 minutes to remove any surface sulfonate groups, prior to a final decant and replacement of the solution with 18 MOÁcm water. Typically samples were stored immersed in water with the top sealed with paraffin wax.
section_name: Substrate Preparation

section_content: Immediately prior to spin coating, each wafer was dried under a nitrogen stream. Once dry, the sample was analyzed using fixed-angle, single-wavelength ellipsometry (632. 8 nm) to determine the thickness of the thermal oxide at the center; typically 2 nm. Spin coating was carried out under argon or nitrogen gas. Each substrate's polished side was evenly coated with 10 μL of 10-15 g/L BCP solution; any bubbles were manually removed; then the substrates were spun for up to 15 s, between 3000 rpm and 4000 rpm, with an initial acceleration of 1500 rpm/s. Following this, the film was reanalyzed by ellipsometry, prior to quartering the sample and annealing.
section_name: Solutions & Spin Coating

section_content: Thermal annealing was carried out in ambient atmosphere on a hotplate covered with a thin aluminum sheet. Temperature was monitored directly at the wafer using an OSENSA fiberoptic fluorescence-based temperature probe. For the thickness measurements and for the comparison of the 5 polymers, the substrates were annealed for 20 minutes at 200°C.
section_name: Annealing

section_content: A solution of 20 mM Na 2 PtCl 4 in 0. 9 M HCl (aq) was used for metallizing PS-b-P2VP samples. Samples were submerged for at least 2-3 hours prior to removal and rinsing with 18. 2 MOÁcm water.
section_name: Metallization

section_content: Following metallization, sample sets were placed together in a plasma chamber, and the chamber was evacuated to < 200 mTorr to remove contaminant gases or adsorbates. Finally, O 2 gas was leaked into the chamber to a pressure of ~750 mTorr. The RF coils were then energized and a faint lavender-blue O 2 plasma was maintained for ~60 s (depending on the film thickness) to etch the organic materials from the substrate. Finally, samples were imaged using a Hitachi S-4800 scanning electron microscope, sampling regions near the centre of each substrate.
section_name: Plasma processing

section_content: For image analysis, ImageJ, [52] version 1. 49 and above, was used. It is freely available at http:// imagej. nih. gov/ij/. The code for performing the analyses is available on our institutional repository; updated versions will be available on GitHub. Python scripts used in preparing the data shown here are also available to assist with processing and plotting output from multiple runs. They are available under an MIT license, allowing users to freely copy, redistribute, and modify the code.
section_name: Computation

section_content: Raw, unedited SEM images, in TIFF format, analyzed for defects, along with summary files of the analysis of each, have been deposited in the University of Alberta's Education and Research Archive and are available via: https://hdl. handle. net/10402/era. 41438. Code for use with ImageJ for analysis of the images, as well as Python scripts for collecting the data across images, will be made available on GitHub under the MIT License at: https://doi. org/10. 5281/zenodo. 19644, and online, with updated versions at: https://github. com/ MurphysLab/ADAblock

section_content: A dot can be thought of as being a junctionless line which has been reduced in length to its width. The +½ defect value associated with each of the two terminal points can be viewed as combining to give the +1 value of the dot. This can also be seen when the dots are paired with-½ defects. The 3-branch junction has a value of-½. Although 3-branch junctions are the most common, (C) 4-branch junctions (-1) and (C) even 5-branch junctions (-1½) can, on rare occasion, be observed. For each additional branch, the value decreases by ½; this can be viewed as being equivalent to sliding an additional branch from a 3-branch junction to increase the junction by 1 branch.