Abstract:
A method of imaging a specimen using ptychography includes directing a charged-particle beam from a source through an illuminator so as to traverse the specimen and land upon a detector, detecting a flux of radiation emanating from the specimen with the detector, calculating at least one property of a charged-particle wavefront exiting the specimen based on using an output of the detector in combination with applying a mathematical reconstruction technique, wherein the at least one property comprises a phase of the wavefront, and wherein applying the mathematical construction technique comprises directly reconstructing the phase of the wavefront to determine a reconstructed phase of the wavefront. An associated apparatus is also described.

Description:
CROSS REFERENCE TO RELATED APPLICATION 
       [0001]    The present application claims priority from pending European Patent Application No. 16179172.8, filed Jul. 13, 2016, which is incorporated herein by reference. 
       BACKGROUND 
       [0002]    The invention relates to a method of imaging a specimen using ptychography, whereby a charged-particle beam is directed from a source through an illuminator so as to traverse the specimen and land upon a detector, an output of the detector being used in combination with a mathematical reconstruction technique so as to calculate at least one property of a charged-particle wavefront exiting the specimen. 
         [0003]    The invention additionally relates to an apparatus for performing such a method. 
         [0004]    The invention further relates to a charged-particle microscope in which such a method can be enacted and/or in which such an apparatus can be comprised. 
         [0005]    Charged-particle microscopy is a well-known and increasingly important technique for imaging microscopic objects, particularly in the form of electron microscopy. Historically, the basic genus of electron microscope has undergone evolution into a number of well-known apparatus species, such as the Transmission Electron Microscope (TEM), Scanning Electron Microscope (SEM), and Scanning Transmission Electron Microscope (STEM), and also into various sub-species, such as so-called “dual-beam” tools (eg. a FIB-SEM), which additionally employ a “machining” Focused Ion Beam (FIB), allowing supportive activities such as ion-beam milling or Ion-Beam-Induced Deposition (IBID), for example. More specifically:
       In a SEM, irradiation of a specimen by a scanning electron beam precipitates emanation of “auxiliary” radiation from the specimen, in the form of secondary electrons, backscattered electrons, X-rays and cathodoluminescence (infrared, visible and/or ultraviolet photons), for example; one or more components of this emanating radiation is/are then detected and used for image accumulation purposes.   In a TEM, the electron beam used to irradiate the specimen is chosen to be of a high-enough energy to penetrate the specimen (which, to this end, will generally be thinner than in the case of a SEM specimen); the transmitted electrons emanating from the specimen can then be used to create an image. When such a TEM is operated in scanning mode (thus becoming a STEM), the image in question will be accumulated during a scanning motion of the irradiating electron beam.
 
More information on some of the topics elucidated here can, for example, be gleaned from the following Wikipedia links:
       http://en.wikipedia.org/wiki/Electron_microscope   http://en.wikipedia.orgiwiki/Scanning_electron_microscope   http://en.wikipedia.org/wiki/Transmission_electron_microscopy   http://en.wikipedia.org/wiki/Scanning_transmission_electron_microscopy
 
As an alternative to the use of electrons as irradiating beam, charged particle microscopy can also be performed using other species of charged particle. In this respect, the phrase “charged particle” should be broadly interpreted as encompassing electrons, positive ions (e.g. Ga or He ions), negative ions, protons and positrons, for instance. As regards non-electron-based charged particle microscopy, some further information can, for example, be gleaned from references such as the following:
   https://en.wikipedia.org/wiki/Focused_ion_beam   http://en.wikipedia.org/wiki/Scanning_Helium_Ion_Microscope   W. H. Escovitz, T. R. Fox and R. Levi-Setti,  Scanning Transmission Ion Microscope with a Field Ion Source , Proc. Nat. Acad. Sci. USA 72(5), pp 1826-1828 (1975).   http://www.ncbi.nlm.nih.gov/pubmed/22472444
 
It should be noted that, in addition to imaging and performing (localized) surface modification (e.g. milling, etching, deposition, etc.), a charged particle microscope may also have other functionalities, such as performing spectroscopy, examining diffractograms, etc.
   
 
         [0016]    In all cases, a Charged-Particle Microscope (CPM) will comprise at least the following components:
       A radiation source, such as a Schottky electron source or ion gun.   An illuminator, which serves to manipulate a “raw” radiation beam from the source and perform upon it certain operations such as focusing, aberration mitigation, cropping (with an aperture), filtering, etc. It will generally comprise one or more (charged-particle) lenses, and may comprise other types of (particle-) optical component also. If desired, the illuminator can be provided with a deflector system that can be invoked to cause its exit beam to perform a scanning motion across the specimen being investigated.   A specimen holder, on which a specimen under investigation can be held and positioned (e.g. tilted, rotated). If desired, this holder can be moved so as to effect scanning motion of the beam w.r.t. the specimen. In general, such a specimen holder will be connected to a positioning system.   A detector (for detecting radiation emanating from an irradiated specimen), which may be unitary or compound/distributed in nature, and which can take many different forms, depending on the radiation being detected. Examples include photodiodes, CMOS detectors, CCD detectors, photovoltaic cells, X-ray detectors (such as Silicon Drift Detectors and Si(Li) detectors), etc. In general, a CPM may comprise several different types of detector, selections of which can be invoked in different situations.
 
In the case of a transmission-type microscope (such as a (S)TEM, for example), a CPM will specifically comprise:
   An imaging system, which essentially takes charged particles that are transmitted through a specimen (plane) and directs (focuses) them onto analysis apparatus, such as a detection/imaging device, spectroscopic apparatus (such as an EELS device), etc. As with the illuminator referred to above, the imaging system may also perform other functions, such as aberration mitigation, cropping, filtering, etc., and it will generally comprise one or more charged-particle lenses and/or other types of particle-optical components.
 
In what follows, the invention may—by way of example—sometimes be set forth in the specific context of electron microscopy; however, such simplification is intended solely for clarity/illustrative purposes, and should not be interpreted as limiting.
       
 
         [0022]    A method as set forth in the opening paragraph above is, for example, elucidated in an article by M. J. Humphry et al. in Nature Communications, 3:730, DOI 10:1038/ncomms1733, Macmillan Publishers Limited (2011), see:
   http://www.nature.com/ncomms/journal/v3/n3/full/ncomms1733.html
 
This article discusses electron beam ptychography, and its application to conduct a form of “lens-less microscopy” in a SEM. The approach disclosed in the article can be regarded as a modification of related techniques from the field of X-ray imaging, where lens-less techniques are attractive because of the difficulty in producing satisfactory X-ray optical systems. In the article, an illuminator produces a convergent electron beam that is focused onto and through a specimen (in a specimen plane) so as to impinge upon a CCD detector. The convergent electron beam is referred to as a “probe”, and this probe is scanned across the specimen in a series of measurement sessions. In each such session, an associated diffraction pattern is recorded by the CCD, and these various patterns are used as input to an iterative mathematical inverse problem solving technique in an attempt to reconstruct the amplitude/phase of the electron-beam wavefront exiting the specimen (somewhat analogous to a deconvolution technique). This, in turn, provides information on the structure of the employed specimen. For more information on this approach, see, for example, the article by J. M. Rodenburg and H. M. L. Faulkner in Appl. Phys. Lett. 85, pp. 4795-4798 (2004) [see the following link]:
   http://scitation.aip.org/content/aip/journal/apl/85/20/10.1063/1.1823034
 
Some general information on ptychography can be gleaned from the following Wikipedia link:
   https://en.wikipedia.org/wiki/Ptychography
 
And an example of a specific ptychographic method/apparatus is set forth in U.S. Pat. No. 9,202,670 (with the same assignee as the present patent application), which is incorporated herein by reference.
   
 
         [0026]    For good order, it is pointed out that what is effectively being reconstructed in such ptychography techniques is a change in a (planar) wavefront as it traverses the specimen. Although ptychography techniques use a relative narrow beam that only illuminates a localized area of the specimen (and that is moved to a series of different locations on the specimen in a corresponding series of measurement sessions), the reconstruction effectively calculates changes to a (virtual) broad wavefront that traverses the full area of the specimen in one go. This point will be fully grasped by the skilled artisan. 
         [0027]    A problem with current ptychographic techniques is that they are subject to so-called “phase wrapping” (also called “phase vortexing”)—which can be particularly problematic for specimens that are relatively thick and/or comprise material with a relatively high atomic number. This phenomenon has to do with the fact that conventional wavefront reconstruction techniques intrinsically limit the reconstructed wavefront phase to a truncated range [0, 2π] (or, equivalently, [−π, +π]); in reality, however, the actual phase can have a value outside this range, in which case the phase will have to be “wrapped up” in order to fit within the truncated range. When the reconstructed phase is used to produce an image of the specimen, such “phase wrapping” can cause distortion of the image (see, for example,  FIG. 2A ); in order to (attempt to) restore such a distorted image, a “phase unwrapping” algorithm must be applied. 
       SUMMARY 
       [0028]    It is an object of the invention to provide an alternative to known ptychographic techniques. In particular, it is an object of the invention that said alternative should mitigate the abovementioned phase wrapping issue. 
         [0029]    These and other objects are achieved in a method as set forth in the opening paragraph above, which method is characterized in that:
       Said property is a phase of the wavefront;   Said mathematical reconstruction technique directly reconstructs said phase, rather than deriving it indirectly from a reconstructed function of amplitude and phase.
 
In line with comments already made above, the term “phase” as referred to here is intended to encompass a phase change to the wavefront.
       
 
     
    
     
       BRIEF DESCRIPTION OF DRAWINGS 
         [0032]    The invention will now be elucidated in more detail on the basis of exemplary embodiments and the accompanying schematic drawings, in which: 
           [0033]      FIG. 1  is a longitudinal cross-sectional view of a CPM, depicted schematically, in which an embodiment of the current invention can be carried out. 
           [0034]      FIGS. 2A and 2B  show ptychographic images obtained using a prior-art technique ( FIG. 2A ) and an embodiment of the present invention ( FIG. 2B ), for comparison purposes. 
       
    
    
       [0035]    In the Figures, where pertinent, corresponding parts are indicated using corresponding reference symbols. It should be noted that, in general, the Figures are not to scale. 
       DETAILED DESCRIPTION 
       [0036]    In traditional mathematical reconstruction techniques used in ptychography (such as the so-called “Rodenburg scheme” [see above]), one attempts to reconstruct a so-called “transmission function” (T), which is a functional expression of wavefront amplitude (A) and phase (φ) at a (vector) point {right arrow over (r)} in (or just after) the specimen plane, and typically takes a form such as: 
         [0000]        T ({right arrow over ( r )})= A ({right arrow over ( r )}) e   iφ({right arrow over (r)})    
         [0000]    For large φ, a series of expansion of the term e iφ(r)  will include higher powers of φ (higher harmonics), thereby distributing power to high spatial frequencies. Consequently, if one would like to limit the reconstructed entity to a given resolution, it will not be possible to successfully filter the transmission function itself—instead, one will need to be able to directly filter the phase. However, with phase jumps present—due to the aforementioned problem of phase wrapping—this cannot be done. The present invention solves this dilemma by directly reconstructing φ rather than indirectly deriving it from a reconstructed transmission function. Apart from addressing the issues described above, this has the additional advantage that, for a relatively thin specimen, φ is proportional to electrical potential, so that the reconstructed entity has direct physical meaning. 
         [0037]    A more mathematical elucidation of the invention can be given as follows. The invention is applicable to any ptychographic scheme that involves iteration over an object, which, in the case of charged particle microscopy, takes the form of a transmission function T(r). This transmission function is, in general, a complex-valued function of (vector) position r, and can be written as: 
         [0000]        T ( r )= A ( r ) e   iφ(t)   (1)
 
         [0000]    where A(r) is an amplitude and φ(r) is the phase. Both A(r) and φ(r) are real-valued functions and, without loss of generality, one can require A(r)&gt;0.
 
For a pure phase object, the amplitude can be neglected. Note that, for typical samples that are used in (S)TEM imaging, the amplitude is close to one for all positions, as most of the electron beam is transmitted through the typically thin samples. In general, it is possible to absorb the amplitude into the imaginary part of φ(r). Therefore, one can re-write the transmission function as:
 
         [0000]        T ( r )= e   iφ(r)   (2)
 
         [0000]    where φ(r) is now understood to be a complex-valued function representing the object. Any iterative Ptychographic reconstruction scheme has the characteristic that it takes an estimate, T n (r), of the transmission function at the current iteration index, n, and updates it with some correction Δ n (r): 
         [0000]        T   n+1 ( r )= T   n ( r )+α n Δ n ( r )  (3)
 
         [0000]    which includes a parameter α n  that controls the (rate of) convergence of the scheme. The subscript n indicates the iteration index and the scheme is initialized with some appropriate guess for T 0 (r).
 
One can now rewrite (3) using (2) as:
 
         [0000]        e   iφ     n+1     (r)   =e   iφ     n     (r) +α n Δ n ( r ),  (4)
 
         [0000]    where φ n (r) is the estimate for φ(r) at iteration index n.
 
Equation (4) can also be re-written as:
 
         [0000]        e   i(φ     n+1     (r)−φ     n     (r)) =1+ e   −iφ     n     (r) α n Δ n ( r ).  (5)
 
         [0000]    The insight at the heart of the present invention is that, physically, one is interested in φ and not in T. It is possible to add (integer multiples of) 2π to the real part of φ (the phase of the transmission function) and obtain the same transmission function from a mathematical point of view. This can even be done at random for any position r and nothing changes. Therefore, mathematically, the phase is usually limited to a range of extent 2π (e.g. [−π,π] or [0, 2π]). This, however, does not make sense physically, as this real part of φ represents the projected potential inside the sample, which is a continuous function that can have any real value.
 
The present invention iterates over the logarithm of (5), assuming that:
       (a) The function T is nowhere zero (which in Ptychography would correspond to either a measurement where one completely blocked the beam, or measuring no signal at all due to detector limitations), and;   (b) The norm of the correction term α n Δ n  is small (at least smaller than 1), which can be arranged by a proper choice of α n .
 
The result of taking the natural logarithm of (5) and approximating it with the lowest term in its series expansion around 1 is:
       
 
         [0000]        i (φ n+1 ( r )−φ n ( r ))=ln(1+ e   −φ     n     (r) α n Δ n ( r ))≈ e   −iφ     n     (r) α n Δ n ( r )  (6)
 
         [0000]      which can be rearranged as: 
         [0000]      φ n+1 ( r )=φ n ( r )− ie   −iφ     n     (r) α n Δ n ( r ).  (7)
 
         [0000]    One important aspect of the described approach is to replace the iteration step (3) with (7), thereby keeping track of the function φ(r). This ensures continuity of (i.e. prevents phase jumps from building up in) the estimate of φ(r) (assuming that both the correction terms and the initial guess are continuous) and allows spatial filtering of the physical object instead of the mathematical object which is the transmission function.
 
This procedure can be applied to a whole class of Ptychographic reconstruction schemes as will be demonstrated below. This can be done even in the case where, in the formulation of the algorithm, there is no specific reference to an iteration scheme as in (3).
 
         [0040]    A specific embodiment of the invention modifies the so-called Rodenburg scheme, e.g. as originally presented in the journal article by A. M. Maiden &amp; J. M. Rodenburg,  An improved ptychographical phase retrieval algorithm for diffractive imaging , Ultramicroscopy, 109 (2009), pp. 1256-1262. In this scheme, there are two main steps: the first one is iteration over estimates of the probe P n (r) (which is another complex-valued function of position), and the second one is iteration over estimates of the object (transmission function). As a starting point, one uses two parameters, α&gt;2 and β&gt;2, and an initial guess of the probe, P 0 (r), and of the object, T 0 (r). Furthermore, one assumes a set of m measurements M j (k) (taken in the diffraction plane with coordinate k) at probe positions r p   j  for j=0 . . . m. For each probe position, an exit wave, ψ 0 (r,r p   j ), which is consistent with the measurements, is defined. This set of exit waves can, for example, be constructed according to equation (12) below, using the initial guesses for the probe and the object. Other choices are also allowed. 
         [0000]    The iteration step for the probe is given by the following update step: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         P 
                         
                           n 
                           + 
                           1 
                         
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     ← 
                     
                       
                         
                           P 
                           n 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                       - 
                       
                         
                           2 
                           
                             α 
                             n 
                           
                         
                          
                         
                           
                             ∑ 
                             
                               j 
                               = 
                               1 
                             
                             m 
                           
                            
                           
                             
                               
                                 
                                   T 
                                   n 
                                 
                                  
                                 
                                   ( 
                                   
                                     r 
                                     + 
                                     
                                       r 
                                       p 
                                       J 
                                     
                                   
                                   ) 
                                 
                               
                               _ 
                             
                              
                             
                               ( 
                               
                                 
                                   
                                     
                                       T 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       
                                         r 
                                         + 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                   
                                   · 
                                   
                                     
                                       P 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       r 
                                       ) 
                                     
                                   
                                 
                                 - 
                                 
                                   
                                     ψ 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       r 
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   where 
                 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       α 
                       n 
                     
                     = 
                     
                       α 
                        
                       
                         
                            
                           
                             
                               ∑ 
                               
                                 j 
                                 = 
                                 1 
                               
                               m 
                             
                              
                             
                               
                                 
                                   
                                     T 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       r 
                                       + 
                                       
                                         r 
                                         p 
                                         J 
                                       
                                     
                                     ) 
                                   
                                 
                                 _ 
                               
                               · 
                               
                                 
                                   T 
                                   n 
                                 
                                  
                                 
                                   ( 
                                   
                                     r 
                                     + 
                                     
                                       r 
                                       p 
                                       j 
                                     
                                   
                                   ) 
                                 
                               
                             
                           
                            
                         
                         ∞ 
                       
                     
                   
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
           
         
       
     
         [0000]    With ∥ X ∥ ∞  denoting the maximum value of the argument X, and the horizontal bar denoting complex conjugation.
 
For the object (transmission function) the iteration step is:
 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         T 
                         
                           n 
                           + 
                           1 
                         
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     ← 
                     
                       
                         
                           T 
                           n 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                       - 
                       
                         
                           2 
                           
                             β 
                             n 
                           
                         
                          
                         
                           
                             ∑ 
                             
                               j 
                               = 
                               1 
                             
                             m 
                           
                            
                           
                             
                               
                                 
                                   P 
                                   n 
                                 
                                  
                                 
                                   ( 
                                   
                                     r 
                                     - 
                                     
                                       r 
                                       p 
                                       J 
                                     
                                   
                                   ) 
                                 
                               
                               _ 
                             
                              
                             
                               ( 
                               
                                 
                                   
                                     
                                       T 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       r 
                                       ) 
                                     
                                   
                                   · 
                                   
                                     
                                       P 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       
                                         r 
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                   
                                 
                                 - 
                                 
                                   
                                     ψ 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         r 
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   where 
                 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       β 
                       n 
                     
                     = 
                     
                       β 
                        
                       
                         
                           
                              
                             
                               
                                 ∑ 
                                 
                                   j 
                                   = 
                                   1 
                                 
                                 m 
                               
                                
                               
                                 
                                   
                                     
                                       P 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       
                                         r 
                                         - 
                                         
                                           r 
                                           p 
                                           J 
                                         
                                       
                                       ) 
                                     
                                   
                                   _ 
                                 
                                 · 
                                 
                                   
                                     P 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       r 
                                       + 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                             
                              
                           
                           ∞ 
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Using the updated probe and object, one can form a new estimate of the output wave for each probe position that satisfies the set of measurements as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ψ 
                       
                         n 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         r 
                         , 
                         
                           r 
                           p 
                           j 
                         
                       
                       ) 
                     
                   
                   ← 
                   
                     
                       ℱ 
                       
                         - 
                         1 
                       
                     
                      
                     
                       { 
                       
                         
                           
                             
                               
                                 
                                   M 
                                   j 
                                 
                                  
                                 
                                   ( 
                                   k 
                                   ) 
                                 
                               
                             
                             · 
                             ℱ 
                           
                            
                           
                             { 
                             
                               
                                 
                                   P 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                                  
                                 
                                   ( 
                                   
                                     
                                       r 
                                       ′ 
                                     
                                     - 
                                     
                                       r 
                                       p 
                                       j 
                                     
                                   
                                   ) 
                                 
                               
                               · 
                               
                                 
                                   T 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                                  
                                 
                                   ( 
                                   
                                     r 
                                     ′ 
                                   
                                   ) 
                                 
                               
                             
                             } 
                           
                            
                           
                             ( 
                             k 
                             ) 
                           
                         
                         
                            
                           
                             ℱ 
                              
                             
                               { 
                               
                                 
                                   
                                     P 
                                     
                                       n 
                                       + 
                                       1 
                                     
                                   
                                   ( 
                                   
                                     
                                       r 
                                       ′ 
                                     
                                     - 
                                     
                                       r 
                                       p 
                                       j 
                                     
                                   
                                   ) 
                                 
                                 · 
                                 
                                   
                                     T 
                                     
                                       n 
                                       + 
                                       1 
                                     
                                   
                                    
                                   
                                     ( 
                                     
                                       r 
                                       ′ 
                                     
                                     ) 
                                   
                                 
                               
                               } 
                             
                              
                             
                               ( 
                               k 
                               ) 
                             
                           
                            
                         
                       
                       } 
                     
                      
                     
                       
                         ( 
                         r 
                         ) 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   12 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Here,          { } and            −1 { } denote a 2D Fourier transform and its inverse, respectively. The forward transform is over the dummy variable r′ in both cases, and the inverse transform is over the dummy variable k.
 
Iteration proceeds until some stopping criterion is fulfilled.
 
The modification according to the present invention is applied in (10), resulting in:
 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ϕ 
                       
                         n 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       r 
                       ) 
                     
                   
                   ← 
                   
                     
                       
                         ϕ 
                         n 
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     + 
                     
                       ɛ 
                        
                       
                           
                       
                        
                       
                         e 
                         
                           
                             - 
                             i 
                           
                            
                           
                               
                           
                            
                           
                             
                               ϕ 
                               n 
                             
                              
                             
                               ( 
                               r 
                               ) 
                             
                           
                         
                       
                        
                       
                         2 
                         
                           β 
                           n 
                         
                       
                        
                       
                         
                           ∑ 
                           
                             j 
                             = 
                             1 
                           
                           m 
                         
                          
                         
                           
                             
                               
                                 P 
                                 n 
                               
                                
                               
                                 ( 
                                 
                                   r 
                                   - 
                                   
                                     r 
                                     p 
                                     J 
                                   
                                 
                                 ) 
                               
                             
                             _ 
                           
                            
                           
                             ( 
                             
                               
                                 
                                   
                                     T 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     r 
                                     ) 
                                   
                                 
                                 · 
                                 
                                   
                                     P 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       r 
                                       - 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               - 
                               
                                 
                                   ψ 
                                   n 
                                 
                                  
                                 
                                   ( 
                                   
                                     
                                       r 
                                       - 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     , 
                                     
                                       r 
                                       p 
                                       j 
                                     
                                   
                                   ) 
                                 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   13 
                   ) 
                 
               
             
           
         
       
     
         [0000]    with ε a parameter that is used to keep the changes in φ small so that (6) holds. The updated φ is subsequently used to form the updated T which is needed in the rest of the scheme, i.e. T n (r)=e iφ     n     (r) . 
         [0041]    Another embodiment of the present invention modifies the so-called PheBIE scheme, e.g. as originally presented in the journal article by R. Hesse, D. R. Luke, S. Sabach &amp; M. K. Tam,  Proximal Heterogeneous Block Implicit - Explicit Method and Application to Blind Ptychographic Diffraction Imaging , SIAM J. Imaging Sciences 8(1) (2015), pp. 426-457. 
         [0000]    In the PHeBIE scheme, essentially the same steps are needed as in the Rodenburg scheme elucidated above, albeit with some adaptations in the coefficients and with a different ordering of the update steps.
 
As a starting point, one has two real-valued functions, α(r)&gt;1 and β(r)&gt;1, a parameter γ&gt;0, an initial guess of the probe, P 0 (r), and of the object, T 0 (r), and a set of m measurements M j (k) at probe positions r p   j  for j=0 . . . m. For each probe position, an exit wave, ψ 0 (r,r p   j ), that is consistent with the measurements is defined. This set of exit waves can, for example, be constructed according to equation (12) above, using the initial guesses for the probe and the object. Other choices are also allowed.
 
The iteration step over the probe is given by:
 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         P 
                         
                           n 
                           + 
                           1 
                         
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     ← 
                     
                       
                         
                           P 
                           n 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                       - 
                       
                         
                           2 
                           
                             
                               α 
                               n 
                             
                              
                             
                               ( 
                               r 
                               ) 
                             
                           
                         
                          
                         
                           
                             ∑ 
                             
                               j 
                               = 
                               1 
                             
                             m 
                           
                            
                           
                             
                               
                                 
                                   T 
                                   n 
                                 
                                  
                                 
                                   ( 
                                   
                                     r 
                                     + 
                                     
                                       r 
                                       p 
                                       J 
                                     
                                   
                                   ) 
                                 
                               
                               _ 
                             
                              
                             
                               ( 
                               
                                 
                                   
                                     
                                       T 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       
                                         r 
                                         + 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                   
                                   · 
                                   
                                     
                                       P 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       r 
                                       ) 
                                     
                                   
                                 
                                 - 
                                 
                                   
                                     ψ 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       r 
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   14 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   where 
                 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       α 
                       n 
                     
                     = 
                     
                       
                         α 
                         ( 
                         r 
                         ) 
                       
                        
                       
                         
                           ∑ 
                           
                             j 
                             = 
                             1 
                           
                           m 
                         
                          
                         
                           
                             
                               
                                 T 
                                 n 
                               
                                
                               
                                 ( 
                                 
                                   r 
                                   + 
                                   
                                     r 
                                     p 
                                     J 
                                   
                                 
                                 ) 
                               
                             
                             _ 
                           
                           · 
                           
                             
                               
                                 T 
                                 n 
                               
                                
                               
                                 ( 
                                 
                                   r 
                                   + 
                                   
                                     r 
                                     p 
                                     j 
                                   
                                 
                                 ) 
                               
                             
                             . 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   15 
                   ) 
                 
               
             
           
         
       
     
         [0000]    For the object (transmission function) one has: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         T 
                         
                           n 
                           + 
                           1 
                         
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     ← 
                     
                       
                         
                           T 
                           n 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                       - 
                       
                         
                           2 
                           
                             
                               β 
                               n 
                             
                              
                             
                               ( 
                               r 
                               ) 
                             
                           
                         
                          
                         
                           
                             ∑ 
                             
                               j 
                               = 
                               1 
                             
                             m 
                           
                            
                           
                             
                               
                                 
                                   P 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                                  
                                 
                                   ( 
                                   
                                     r 
                                     - 
                                     
                                       r 
                                       p 
                                       J 
                                     
                                   
                                   ) 
                                 
                               
                               _ 
                             
                              
                             
                               ( 
                               
                                 
                                   
                                     
                                       P 
                                       
                                         n 
                                         + 
                                         1 
                                       
                                     
                                      
                                     
                                       ( 
                                       
                                         r 
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                   
                                   · 
                                   
                                     
                                       T 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       r 
                                       ) 
                                     
                                   
                                 
                                 - 
                                 
                                   
                                     ψ 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         r 
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   16 
                   ) 
                 
               
             
           
         
       
     
         [0000]    which now immediately uses the updated probe (14) and where 
         [0000]      β n ( r )=β( r )Σ j=1   m     P   n+1 ( r−r   p   j ) · P   n+1 ( r−r   p   j ).  (17)
 
         [0000]    Note that, compared to the Rodenburg scheme, the coefficients α n  and β n  have now become real-valued functions of position instead of constants. Furthermore, in (16), the updated probe is used instead of the previous estimate of the probe.
 
Using the updated probe and object, a new estimate of the output wave for each probe position is formed using the set of measurements as:
 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ψ 
                       
                         n 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         r 
                         , 
                         
                           r 
                           p 
                           j 
                         
                       
                       ) 
                     
                   
                   ← 
                   
                     
                       ℱ 
                       
                         - 
                         1 
                       
                     
                      
                     
                       { 
                       
                         
                           
                             
                               
                                 
                                   M 
                                   j 
                                 
                                  
                                 
                                   ( 
                                   k 
                                   ) 
                                 
                               
                             
                             · 
                             ℱ 
                           
                            
                           
                             { 
                             
                               
                                 
                                   2 
                                   
                                     2 
                                     + 
                                     γ 
                                   
                                 
                                  
                                 
                                   
                                     
                                       P 
                                       
                                         n 
                                         + 
                                         1 
                                       
                                     
                                      
                                     
                                       ( 
                                       
                                         
                                           r 
                                           ′ 
                                         
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                   
                                   · 
                                   
                                     
                                       T 
                                       
                                         n 
                                         + 
                                         1 
                                       
                                     
                                      
                                     
                                       ( 
                                       
                                         r 
                                         ′ 
                                       
                                       ) 
                                     
                                   
                                 
                               
                               + 
                               
                                 
                                   γ 
                                   
                                     2 
                                     + 
                                     γ 
                                   
                                 
                                  
                                 
                                   
                                     ψ 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         r 
                                         ′ 
                                       
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                             
                             } 
                           
                            
                           
                             ( 
                             k 
                             ) 
                           
                         
                         
                            
                           
                             ℱ 
                              
                             
                               { 
                               
                                 
                                   
                                     2 
                                     
                                       2 
                                       + 
                                       γ 
                                     
                                   
                                    
                                   
                                     
                                       
                                         P 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                       ( 
                                       
                                         
                                           r 
                                           ′ 
                                         
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                     · 
                                     
                                       
                                         T 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                        
                                       
                                         ( 
                                         
                                           r 
                                           ′ 
                                         
                                         ) 
                                       
                                     
                                   
                                 
                                 + 
                                 
                                   
                                     γ 
                                     
                                       2 
                                       + 
                                       γ 
                                     
                                   
                                    
                                   
                                     
                                       ψ 
                                       n 
                                     
                                     ( 
                                     
                                       
                                         r 
                                         ′ 
                                       
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               } 
                             
                              
                             
                               ( 
                               k 
                               ) 
                             
                           
                            
                         
                       
                       } 
                     
                      
                     
                       
                         ( 
                         r 
                         ) 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   18 
                   ) 
                 
               
             
           
         
       
     
         [0000]    The parameter γ is used to influence the convergence of the scheme. Iteration proceeds until a stopping criterion is fulfilled.
 
The modification according to the present invention is applied in (16), resulting in:
 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         ϕ 
                         
                           n 
                           + 
                           1 
                         
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     ← 
                     
                       
                         
                           ϕ 
                           n 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                       + 
                       
                         ɛ 
                          
                         
                             
                         
                          
                         
                           e 
                           
                             
                               - 
                               i 
                             
                              
                             
                                 
                             
                              
                             
                               
                                 ϕ 
                                 n 
                               
                                
                               
                                 ( 
                                 r 
                                 ) 
                               
                             
                           
                         
                          
                         
                           2 
                           
                             
                               β 
                               n 
                             
                              
                             
                               ( 
                               r 
                               ) 
                             
                           
                         
                          
                         
                           
                             ∑ 
                             
                               j 
                               = 
                               1 
                             
                             m 
                           
                            
                           
                             
                               
                                 
                                   P 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                                  
                                 
                                   ( 
                                   
                                     r 
                                     - 
                                     
                                       r 
                                       p 
                                       J 
                                     
                                   
                                   ) 
                                 
                               
                               _ 
                             
                              
                             
                               ( 
                               
                                 
                                   
                                     P 
                                     
                                       n 
                                       + 
                                       1 
                                     
                                   
                                    
                                   
                                     
                                       ( 
                                       
                                         r 
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                     · 
                                     
                                       
                                         T 
                                         n 
                                       
                                        
                                       
                                         ( 
                                         r 
                                         ) 
                                       
                                     
                                   
                                 
                                 - 
                                 
                                   
                                     ψ 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         r 
                                         - 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   19 
                   ) 
                 
               
             
           
         
       
     
         [0000]    with ε a parameter that is used to keep the changes in φ small so that (6) holds. The updated φ is subsequently used to form the updated T, which is needed in the rest of the scheme, i.e. T n (r)=e iφ     n     (r) . 
         [0042]    Yet another embodiment of the present invention modifies the so-called Thibault scheme, e.g. as originally presented in the journal article by P. Thibault, M. Dierolf, O. Bunk, A. Menzel &amp; F. Pfeiffer,  Probe retrieval in ptychographic coherent diffractive imaging,  Ultramicroscopy, 109 (2009), pp. 338-343. The Thibault algorithm uses two nested iteration loops, whereby the outer iteration loop uses the outcome of the inner iteration loop to update the output waves at each probe position. 
         [0000]    As a starting point, one has an initial guess of the probe P 0 (r) and of the object T 0 (r), and a set of m measurements M j (k) at probe positions r p   j  for j=0 . . . m. For each probe position, an exit wave φ 0 (r,r p   j ), which is consistent with the measurements, is defined. This set of exit waves can, for example, be constructed according to equation (12) above, using the initial guesses for the probe and the object. Other choices are also allowed.
 
The inner iteration loop does a number Λ of iterations (indicated with the superscript l) on the probe, initialized with P n   0  (r)=P n (r), and the object, initialized with T n   0 (r)=T n (r), as given by:
 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       P 
                       n 
                       
                         l 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       r 
                       ) 
                     
                   
                   ← 
                   
                     
                       1 
                       
                         
                           α 
                           l 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                     
                      
                     
                       
                         ∑ 
                         
                           j 
                           = 
                           1 
                         
                         m 
                       
                        
                       
                         
                           
                             
                               T 
                               n 
                               l 
                             
                              
                             
                               ( 
                               
                                 r 
                                 + 
                                 
                                   r 
                                   p 
                                   J 
                                 
                               
                               ) 
                             
                           
                           _ 
                         
                         · 
                         
                           
                             ψ 
                             n 
                           
                            
                           
                             ( 
                             
                               r 
                               , 
                               
                                 r 
                                 p 
                                 j 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   20 
                   ) 
                 
               
             
             
               
                 and 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                     
                       T 
                       n 
                       
                         l 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       r 
                       ) 
                     
                   
                   ← 
                   
                     
                       1 
                       
                         
                           β 
                           l 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                     
                      
                     
                       
                         ∑ 
                         
                           j 
                           = 
                           1 
                         
                         m 
                       
                        
                       
                         
                           
                             
                               P 
                               n 
                               l 
                             
                              
                             
                               ( 
                               
                                 r 
                                 - 
                                 
                                   r 
                                   p 
                                   J 
                                 
                               
                               ) 
                             
                           
                           _ 
                         
                         · 
                         
                           
                             ψ 
                             n 
                           
                            
                           
                             ( 
                             
                               
                                 r 
                                 - 
                                 
                                   r 
                                   p 
                                   j 
                                 
                               
                               , 
                               
                                 r 
                                 p 
                                 j 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   21 
                   ) 
                 
               
             
             
               
                 where 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                     
                       α 
                       l 
                     
                      
                     
                       ( 
                       r 
                       ) 
                     
                   
                   = 
                   
                     2 
                      
                     
                       
                         ∑ 
                         
                           j 
                           = 
                           1 
                         
                         m 
                       
                        
                       
                         
                           
                             
                               T 
                               n 
                               l 
                             
                              
                             
                               ( 
                               
                                 r 
                                 + 
                                 
                                   r 
                                   p 
                                   J 
                                 
                               
                               ) 
                             
                           
                           _ 
                         
                         · 
                         
                           
                             T 
                             n 
                             l 
                           
                            
                           
                             ( 
                             
                               r 
                               + 
                               
                                 r 
                                 p 
                                 j 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   22 
                 
               
             
             
               
                 and 
               
               
                 
                     
                 
               
             
             
               
                 
                   
                     
                       β 
                       l 
                     
                      
                     
                       ( 
                       r 
                       ) 
                     
                   
                   = 
                   
                     2 
                      
                     
                       
                         ∑ 
                         
                           j 
                           = 
                           1 
                         
                         m 
                       
                        
                       
                         
                           
                             
                               P 
                               n 
                               l 
                             
                              
                             
                               ( 
                               
                                 r 
                                 - 
                                 
                                   r 
                                   p 
                                   J 
                                 
                               
                               ) 
                             
                           
                           _ 
                         
                         · 
                         
                           
                             
                               P 
                               n 
                               l 
                             
                              
                             
                               ( 
                               
                                 r 
                                 + 
                                 
                                   r 
                                   p 
                                   j 
                                 
                               
                               ) 
                             
                           
                           . 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   23 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Note that (22) and (23) resemble (15) and (17), respectively, with α(r)=β(r)=2.
 
At the end of this set of iterations, the results are an updated probe and object given by P n+1 (r)=P n   Λ+1 (r) and T n+1 (r)=T n   Λ+1 (r), respectively. These are then used in the outer iteration loop to form the following updates for each probe position:
 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ψ 
                       
                         n 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         r 
                         , 
                         
                           r 
                           p 
                           j 
                         
                       
                       ) 
                     
                   
                   ← 
                   
                     
                       
                         ψ 
                         n 
                       
                        
                       
                         ( 
                         
                           r 
                           , 
                           
                             r 
                             p 
                             j 
                           
                         
                         ) 
                       
                     
                     + 
                     
                       
                         ℱ 
                         
                           - 
                           1 
                         
                       
                        
                       
                         { 
                         
                           
                             
                               
                                 
                                   
                                     M 
                                     j 
                                   
                                    
                                   
                                     ( 
                                     k 
                                     ) 
                                   
                                 
                               
                               · 
                               ℱ 
                             
                              
                             
                               { 
                               
                                 
                                   2 
                                    
                                   
                                     
                                       
                                         P 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                        
                                       
                                         ( 
                                         
                                           
                                             r 
                                             ′ 
                                           
                                           - 
                                           
                                             r 
                                             p 
                                             j 
                                           
                                         
                                         ) 
                                       
                                     
                                     · 
                                     
                                       
                                         T 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                        
                                       
                                         ( 
                                         
                                           r 
                                           ′ 
                                         
                                         ) 
                                       
                                     
                                   
                                 
                                 - 
                                 
                                   
                                     ψ 
                                     n 
                                   
                                    
                                   
                                     ( 
                                     
                                       r 
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               } 
                             
                              
                             
                               ( 
                               k 
                               ) 
                             
                           
                           
                              
                             
                               ℱ 
                                
                               
                                 { 
                                 
                                   
                                     2 
                                      
                                     
                                       
                                         
                                           P 
                                           
                                             n 
                                             + 
                                             1 
                                           
                                         
                                         ( 
                                         
                                           
                                             r 
                                             ′ 
                                           
                                           - 
                                           
                                             r 
                                             p 
                                             j 
                                           
                                         
                                         ) 
                                       
                                       · 
                                       
                                         
                                           T 
                                           
                                             n 
                                             + 
                                             1 
                                           
                                         
                                          
                                         
                                           ( 
                                           
                                             r 
                                             ′ 
                                           
                                           ) 
                                         
                                       
                                     
                                   
                                   - 
                                   
                                     
                                       ψ 
                                       n 
                                     
                                     ( 
                                     
                                       r 
                                       , 
                                       
                                         r 
                                         p 
                                         j 
                                       
                                     
                                     ) 
                                   
                                 
                                 } 
                               
                                
                               
                                 ( 
                                 k 
                                 ) 
                               
                             
                              
                           
                         
                         } 
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     - 
                     
                       
                         P 
                         
                           n 
                           + 
                           1 
                         
                       
                       · 
                       
                         
                           
                             T 
                             
                               n 
                               + 
                               1 
                             
                           
                            
                           
                             ( 
                             r 
                             ) 
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   24 
                   ) 
                 
               
             
           
         
       
     
         [0000]    This algorithm can be adapted according to the present invention by first rewriting (21) as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       T 
                       n 
                       
                         l 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       r 
                       ) 
                     
                   
                   ← 
                   
                     
                       
                         T 
                         n 
                         l 
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     - 
                     
                       
                         T 
                         n 
                         l 
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     + 
                     
                       
                         1 
                         
                           
                             α 
                             l 
                           
                            
                           
                             ( 
                             r 
                             ) 
                           
                         
                       
                        
                       
                         
                           ∑ 
                           
                             j 
                             = 
                             1 
                           
                           m 
                         
                          
                         
                           
                             
                               
                                 P 
                                 n 
                                 l 
                               
                                
                               
                                 ( 
                                 
                                   r 
                                   - 
                                   
                                     r 
                                     p 
                                     J 
                                   
                                 
                                 ) 
                               
                             
                             _ 
                           
                           · 
                           
                             
                               
                                 ψ 
                                 n 
                               
                                
                               
                                 ( 
                                 
                                   
                                     r 
                                     - 
                                     
                                       r 
                                       p 
                                       j 
                                     
                                   
                                   , 
                                   
                                     r 
                                     p 
                                     j 
                                   
                                 
                                 ) 
                               
                             
                             . 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   25 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Next, (25) can be rewritten as an iteration over the phase φ according to (6): 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         ϕ 
                         
                           l 
                           + 
                           1 
                         
                       
                        
                       
                         ( 
                         r 
                         ) 
                       
                     
                     ← 
                     
                       
                         
                           ϕ 
                           l 
                         
                          
                         
                           ( 
                           r 
                           ) 
                         
                       
                       - 
                       
                         ɛ 
                          
                         
                             
                         
                          
                         
                           
                             e 
                             
                               
                                 - 
                                 i 
                               
                                
                               
                                   
                               
                                
                               
                                 
                                   ϕ 
                                   l 
                                 
                                  
                                 
                                   ( 
                                   r 
                                   ) 
                                 
                               
                             
                           
                           ( 
                           
                             
                               - 
                               
                                 
                                   T 
                                   n 
                                   l 
                                 
                                  
                                 
                                   ( 
                                   r 
                                   ) 
                                 
                               
                             
                             + 
                             
                               
                                 1 
                                 
                                   
                                     α 
                                     l 
                                   
                                    
                                   
                                     ( 
                                     r 
                                     ) 
                                   
                                 
                               
                                
                               
                                 
                                   ∑ 
                                   
                                     j 
                                     = 
                                     1 
                                   
                                   m 
                                 
                                  
                                 
                                   
                                     
                                       
                                         P 
                                         n 
                                         l 
                                       
                                        
                                       
                                         ( 
                                         
                                           r 
                                           - 
                                           
                                             r 
                                             p 
                                             J 
                                           
                                         
                                         ) 
                                       
                                     
                                     _ 
                                   
                                   · 
                                   
                                     
                                       ψ 
                                       n 
                                     
                                      
                                     
                                       ( 
                                       
                                         
                                           r 
                                           - 
                                           
                                             r 
                                             p 
                                             j 
                                           
                                         
                                         , 
                                         
                                           r 
                                           p 
                                           j 
                                         
                                       
                                       ) 
                                     
                                   
                                 
                               
                             
                           
                           ) 
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   26 
                   ) 
                 
               
             
           
         
       
     
         [0000]    with ε a parameter that is used to keep the changes in φ small so that (6) holds. The updated φ is subsequently used to form the updated T that is needed in the rest of the scheme, i.e. T n   l (r)=e iφ     l     (r) . 
         [0043]    The examples above are non-exhaustive. For more information on (unmodified) prior-art ptychographic reconstruction techniques, reference is made to the journal article by R. Hesse et al.,  Proximal Heterogeneous Block Implicit - Explicit Method and Application to Blind Ptychographic Diffractive Imaging , SIAM J. Imaging Sciences 8(1), 2015, pp. 426-457 [see the following link]:
   https://www.researchgate.net/publication/264624024_Proximal_Heterogeneous_Block_Implicit-Explicit_Method_and_Application_to_Blind_Ptychographic_Diffraction_Imaging   
 
         [0045]    It should be noted that, although charged-particle ptychography can be conveniently carried out in a CPM, one is not limited to use of a CPM in this context, and one can instead choose to construct/use a dedicated ptychography apparatus, if one so wishes. 
       Embodiment 1 
       [0046]      FIG. 1  is a highly schematic depiction of an embodiment of a CPM M according to the current invention, which, in this case, is a (S)TEM (though, in the context of the current invention, it could just as validly be an ion-based or proton microscope, for example, or a SEM operating in transmissive mode (TSEM)). In  FIG. 1 , within a vacuum enclosure E, an electron source  4  (such as a Schottky emitter, for example) produces a beam (B) of electrons that traverse an electron-optical illuminator  6 , serving to direct/focus them onto a chosen part of a specimen S (which may, for example, be (locally) thinned/planarized). This illuminator  6  has an electron-optical axis B′, and will generally comprise a variety of electrostatic/magnetic lenses, (scan) deflector(s) D, correctors (such as stigmators), etc.; typically, it can also comprise a condenser system (in fact, the whole of item  6  is sometimes referred to as “a condenser system”). 
         [0047]    The specimen S is held on a specimen holder H. As here illustrated, part of this holder H (inside enclosure E) is mounted in a cradle A′ that can be positioned/moved in multiple degrees of freedom by a positioning device (stage) A; for example, the cradle A′ may (inter alia) be displaceable in the X, Y and Z directions (see the depicted Cartesian coordinate system), and may be rotated about a longitudinal axis parallel to X. Such movement allows different parts of the specimen S to be irradiated/imaged/inspected by the electron beam traveling along axis B′ (and/or allows scanning motion to be performed as an alternative to beam scanning [using deflector(s) D], and/or allows selected parts of the specimen S to be machined by a (non-depicted) focused ion beam, for example). 
         [0048]    The (focused) electron beam B traveling along axis B′ will interact with the specimen S in such a manner as to cause various types of “stimulated” radiation to emanate from the specimen S, including (for example) secondary electrons, backscattered electrons, X-rays and optical radiation (cathodoluminescence). If desired, one or more of these radiation types can be detected with the aid of sensor  22 , which might be a combined scintillator/photomultiplier or EDX (Energy-Dispersive X-Ray Spectroscopy) module, for instance; in such a case, an image could be constructed using basically the same principle as in a SEM. However, of principal importance in a (S)TEM, one can instead/supplementally study electrons that traverse (pass through) the specimen  5 , emerge (emanate) from it and continue to propagate (substantially, though generally with some deflection/scattering) along axis B′. Such a transmitted electron flux enters an imaging system (combined objective/projection lens)  24 , which will generally comprise a variety of electrostatic/magnetic lenses, deflectors, correctors (such as stigmators), etc. In normal (non-scanning) TEM mode, this imaging system  24  can focus the transmitted electron flux onto a fluorescent screen  26 , which, if desired, can be retracted/withdrawn (as schematically indicated by arrows  26 ′) so as to get it out of the way of axis B′. An image (or diffractogram) of (part of) the specimen S will be formed by imaging system  24  on screen  26 , and this may be viewed through viewing port  28  located in a suitable part of a wall of enclosure E. The retraction mechanism for screen  26  may, for example, be mechanical and/or electrical in nature, and is not depicted here. 
         [0049]    As an alternative to viewing an image on screen  26 , one can instead make use of the fact that the depth of focus of the electron flux emerging from imaging system  24  is generally quite large (e.g. of the order of 1 meter). Consequently, various types of sensing device/analysis apparatus can be used downstream of screen  26 , such as:
       TEM camera  30 . At camera  30 , the electron flux can form a static image (or diffractogram) that can be processed by controller  10  and displayed on a display device (not depicted), such as a flat panel display, for example. When not required, camera  30  can be retracted/withdrawn (as schematically indicated by arrows  30 ′) so as to get it out of the way of axis B′.   STEM detector  32 . An output from detector  32  can be recorded as a function of (X,Y) scanning position of the beam B on the specimen  5 , and an image can be constructed that is a “map” of output from detector  32  as a function of X,Y. Typically, detector  32  will have a much higher acquisition rate (e.g. 106 points per second) than camera  30  (e.g. 102 images per second). In conventional tools, detector  32  can comprise a single pixel with a diameter of e.g. 20 mm, as opposed to the matrix of pixels characteristically present in camera  30 . Once again, when not required, detector  32  can be retracted/withdrawn (as schematically indicated by arrows  32 ′) so as to get it out of the way of axis B′ (although such retraction would not be a necessity in the case of a donut-shaped annular dark field detector  32 , for example; in such a detector, a central hole would allow beam passage when the detector was not in use).   As an alternative to imaging using camera  30  or detector  32 , one can also invoke spectroscopic apparatus  34 , which could be an EELS module, for example.
 
It should be noted that the order/location of items  30 ,  32  and  34  is not strict, and many possible variations are conceivable. For example, spectroscopic apparatus  34  can also be integrated into the imaging system  24 .
       
 
         [0053]    Note that the controller/computer processor  10  is connected to various illustrated components via control lines (buses)  10 ′. This controller  10  can provide a variety of functions, such as synchronizing actions, providing setpoints, processing signals, performing calculations, and displaying messages/information on a display device (not depicted). Needless to say, the (schematically depicted) controller  10  may be (partially) inside or outside the enclosure E, and may have a unitary or composite structure, as desired. The skilled artisan will understand that the interior of the enclosure E does not have to be kept at a strict vacuum; for example, in a so-called “Environmental (S)TEM”, a background atmosphere of a given gas is deliberately introduced/maintained within the enclosure E. The skilled artisan will also understand that, in practice, it may be advantageous to confine the volume of enclosure E so that, where possible, it essentially hugs the axis B′, taking the form of a small tube (e.g. of the order of 1 cm in diameter) through which the employed electron beam passes, but widening out to accommodate structures such as the source  4 , specimen holder H, screen  26 , camera  30 , detector  32 , spectroscopic apparatus  34 , etc. 
         [0054]    In the context of the current invention, the illuminator  6  may, for example, comprise a beam shaping element such as a lens and/or aperture plate/diaphragm, to appropriately shape (focus) the beam B into a relatively narrow “pencil” of charged particles, causing it to irradiate only a relatively small area (footprint) of the specimen S at any given time. Relative motion between the specimen S and the footprint of the beam B (so as to move the footprint onto another area of the specimen) can be produced by one or more of the following methods:
       Using the positioning device A to move the specimen S relative to the beam B;   Using the deflector D to deflect the beam B relative to the specimen S;   Moving the source  4  or/and the abovementioned beam shaping element so as to displace the beam B relative to the specimen S.
 
For each such chosen position of the (footprint of the) beam B relative to the specimen S, one can (for example) use the TEM camera  30  [(position-dependent) intensity detector] to capture a diffractogram (diffraction pattern). Specifically, the controller  10  (or another processor device) can be configured/invoked to:
   Acquire a ptychographic measurement set, by recording an output of (for example) TEM camera  30  for each of a series of different positions of the beam B upon the specimen S (achieved, for example, by sending an appropriate series of setpoints to items A, D, etc.).   Process the recorded outputs of said camera  30 , and use them as input to perform a mathematical reconstruction algorithm according to the invention. The skilled artisan in the field of ptychography will understand that this step typically involves calculating (mean) intensity values of features in the recorded diffractograms: see, for example, the references cited above.   Display the results of said reconstruction algorithm, e.g. in the form of an image on a display device (not shown).       
 
       Embodiment 2 
       [0061]      FIGS. 2A and 2B  show ptychographic images obtained using a prior-art technique ( FIG. 2A ) and an embodiment of the present invention ( FIG. 2B ), for comparison purposes. Both images show a GaN crystalline lattice, imaged in a TEM with the following settings/parameters:
       Specimen thickness: 5 nm.   Input (primary/probing) beam energy: 300 kV.   Field of view: 5 nm.   Defocus: 70 nm from disc of least confusion.
 
 FIG. 2A  shows the results of a standard Rodenburg reconstruction algorithm. Note the artifacts caused by phase wrapping, which cause the relatively heavy Ga atoms (atomic number Z=31) to appear as contorted, irregular, dark “loops” surrounded by bright halos. In contrast,  FIG. 2B  shows the results of a direct phase-reconstructing algorithm according to the present invention. Note in this case that the Ga atoms take the form of uniform dots.
       
 
         [0066]    In view of the many possible embodiments to which the principles of the disclosed technology may be applied, it should be recognized that the illustrated embodiments are only preferred examples of the disclosed technology and should not be taken as limiting the scope of protection. Rather, the scope of protection is defined by the following claims. We therefore claim all that comes within the scope of these claims.