-
Parity non-conservation in a condensed matter system
Authors:
Veit Elser
Abstract:
The nuclear spin of a He$^3$ quasiparticle dissolved in superfluid He$^4$ sees an apparent magnetic field proportional to the Fermi coupling constant, the superfluid condensate density, and the electron current at the He$^3$ nucleus. Whereas the direction of the current must be parallel to the quasiparticle momentum, calculating its magnitude presents an interesting theoretical challenge because i…
▽ More
The nuclear spin of a He$^3$ quasiparticle dissolved in superfluid He$^4$ sees an apparent magnetic field proportional to the Fermi coupling constant, the superfluid condensate density, and the electron current at the He$^3$ nucleus. Whereas the direction of the current must be parallel to the quasiparticle momentum, calculating its magnitude presents an interesting theoretical challenge because it vanishes in the Born-Oppenheimer approximation. We find the effect is too small to be observed and present our results in the hope others will be inspired to look for similar effects in other systems.
△ Less
Submitted 17 February, 2020;
originally announced February 2020.
-
Benchmark problems for phase retrieval
Authors:
Veit Elser,
Ti-Yen Lan,
Tamir Bendory
Abstract:
In recent years, the mathematical and algorithmic aspects of the phase retrieval problem have received considerable attention. Many papers in this area mention crystallography as a principal application. In crystallography, the signal to be recovered is periodic and comprised of atomic distributions arranged homogeneously in the unit cell of the crystal. The crystallographic problem is both the le…
▽ More
In recent years, the mathematical and algorithmic aspects of the phase retrieval problem have received considerable attention. Many papers in this area mention crystallography as a principal application. In crystallography, the signal to be recovered is periodic and comprised of atomic distributions arranged homogeneously in the unit cell of the crystal. The crystallographic problem is both the leading application and one of the hardest forms of phase retrieval. We have constructed a graded set of benchmark problems for evaluating algorithms that perform this type of phase retrieval. The data, publicly available online, is provided in an easily interpretable format. We also propose a simple and unambiguous success/failure criterion based on the actual needs in crystallography. Baseline runtimes were obtained with an iterative algorithm that is similar but more transparent than those used in crystallography. Empirically, the runtimes grow exponentially with respect to a new hardness parameter: the sparsity of the signal autocorrelation. We also review the algorithms used by the leading software packages. This set of benchmark problems, we hope, will encourage the development of new algorithms for the phase retrieval problem in general, and crystallography in particular.
△ Less
Submitted 14 June, 2018; v1 submitted 1 June, 2017;
originally announced June 2017.
-
Breaking the Crowther Limit: Combining Depth-Sectioning and Tilt Tomography for High-Resolution, Wide-Field 3D Reconstructions
Authors:
Robert Hovden,
Peter Ercius,
Yi Jiang,
Deli Wang,
Yingchao Yu,
Hector D. Abruna,
Veit Elser,
David A. Muller
Abstract:
To date, high-resolution (< 1 nm) imaging of extended objects in three-dimensions (3D) has not been possible. A restriction known as the Crowther criterion forces a tradeoff between object size and resolution for 3D reconstructions by tomography. Further, the sub-Angstrom resolution of aberration-corrected electron microscopes is accompanied by a greatly diminished depth of field, causing regions…
▽ More
To date, high-resolution (< 1 nm) imaging of extended objects in three-dimensions (3D) has not been possible. A restriction known as the Crowther criterion forces a tradeoff between object size and resolution for 3D reconstructions by tomography. Further, the sub-Angstrom resolution of aberration-corrected electron microscopes is accompanied by a greatly diminished depth of field, causing regions of larger specimens (> 6 nm) to appear blurred or missing. Here we demonstrate a three-dimensional imaging method that overcomes both these limits by combining through-focal depth sectioning and traditional tilt-series tomography to reconstruct extended objects, with high-resolution, in all three dimensions. The large convergence angle in aberration corrected instruments now becomes a benefit and not a hindrance to higher quality reconstructions. A through-focal reconstruction over a 390 nm 3D carbon support containing over one hundred dealloyed and nanoporous PtCu catalyst particles revealed with sub-nanometer detail the extensive and connected interior pore structure that is created by the dealloying instability.
△ Less
Submitted 31 January, 2014;
originally announced February 2014.
-
Real-Space x-ray tomographic reconstruction of randomly oriented objects with sparse data frames
Authors:
Kartik Ayyer,
Hugh T. Philipp,
Mark W. Tate,
Veit Elser,
Sol M. Gruner
Abstract:
Schemes for X-ray imaging single protein molecules using new x-ray sources, like x-ray free electron lasers (XFELs), require processing many frames of data that are obtained by taking temporally short snapshots of identical molecules, each with a random and unknown orientation. Due to the small size of the molecules and short exposure times, average signal levels of much less than 1 photon/pixel/f…
▽ More
Schemes for X-ray imaging single protein molecules using new x-ray sources, like x-ray free electron lasers (XFELs), require processing many frames of data that are obtained by taking temporally short snapshots of identical molecules, each with a random and unknown orientation. Due to the small size of the molecules and short exposure times, average signal levels of much less than 1 photon/pixel/frame are expected, much too low to be processed using standard methods. One approach to process the data is to use statistical methods developed in the EMC algorithm (Loh & Elser, Phys. Rev. E, 2009) which processes the data set as a whole. In this paper we apply this method to a real-space tomographic reconstruction using sparse frames of data (below $10^{-2}$ photons/pixel/frame) obtained by performing x-ray transmission measurements of a low-contrast, randomly-oriented object. This extends the work by Philipp et al. (Optics Express, 2012) to three dimensions and is one step closer to the single molecule reconstruction problem.
△ Less
Submitted 7 November, 2013;
originally announced November 2013.
-
Direct Phasing of Nanocrystal Diffraction
Authors:
Veit Elser
Abstract:
Recent experiments at free-electron laser x-ray sources have been able to resolve the intensity distributions about Bragg peaks in nanocrystals of large biomolecules. Information derived from small shifts in the peak positions augment the Bragg samples of the particle intensity with samples of its gradients. Working on the assumption that the nanocrystal is entirely generated by lattice translatio…
▽ More
Recent experiments at free-electron laser x-ray sources have been able to resolve the intensity distributions about Bragg peaks in nanocrystals of large biomolecules. Information derived from small shifts in the peak positions augment the Bragg samples of the particle intensity with samples of its gradients. Working on the assumption that the nanocrystal is entirely generated by lattice translations of a particle, we develop an algorithm that reconstructs the particle from intensities and intensity gradients. Unlike traditional direct phasing methods that require very high resolution data in order to exploit sparsity of the electron density, our method imposes no constraints on the contrast other than positivity and works well at low resolution. We demonstrate successful reconstructions with simulated P1 lysozyme nanocrystal data down to a signal-to-noise ratio of 2 in the intensity gradients.
△ Less
Submitted 4 June, 2013;
originally announced June 2013.
-
An Improved Three-Weight Message-Passing Algorithm
Authors:
Nate Derbinsky,
José Bento,
Veit Elser,
Jonathan S. Yedidia
Abstract:
We describe how the powerful "Divide and Concur" algorithm for constraint satisfaction can be derived as a special case of a message-passing version of the Alternating Direction Method of Multipliers (ADMM) algorithm for convex optimization, and introduce an improved message-passing algorithm based on ADMM/DC by introducing three distinct weights for messages, with "certain" and "no opinion" weigh…
▽ More
We describe how the powerful "Divide and Concur" algorithm for constraint satisfaction can be derived as a special case of a message-passing version of the Alternating Direction Method of Multipliers (ADMM) algorithm for convex optimization, and introduce an improved message-passing algorithm based on ADMM/DC by introducing three distinct weights for messages, with "certain" and "no opinion" weights, as well as the standard weight used in ADMM/DC. The "certain" messages allow our improved algorithm to implement constraint propagation as a special case, while the "no opinion" messages speed convergence for some problems by making the algorithm focus only on active constraints. We describe how our three-weight version of ADMM/DC can give greatly improved performance for non-convex problems such as circle packing and solving large Sudoku puzzles, while retaining the exact performance of ADMM for convex problems. We also describe the advantages of our algorithm compared to other message-passing algorithms based upon belief propagation.
△ Less
Submitted 8 May, 2013;
originally announced May 2013.
-
Solving Structure with Sparse, Randomly-Oriented X-ray Data
Authors:
Hugh T. Philipp,
Kartik Ayyer,
Mark W. Tate,
Veit Elser,
Sol M. Gruner
Abstract:
Single-particle imaging experiments of biomolecules at x-ray free-electron lasers (XFELs) require processing of hundreds of thousands (or more) of images that contain very few x-rays. Each low-flux image of the diffraction pattern is produced by a single, randomly oriented particle, such as a protein. We demonstrate the feasibility of collecting data at these extremes, averaging only 2.5 photons p…
▽ More
Single-particle imaging experiments of biomolecules at x-ray free-electron lasers (XFELs) require processing of hundreds of thousands (or more) of images that contain very few x-rays. Each low-flux image of the diffraction pattern is produced by a single, randomly oriented particle, such as a protein. We demonstrate the feasibility of collecting data at these extremes, averaging only 2.5 photons per frame, where it seems doubtful there could be information about the state of rotation, let alone the image contrast. This is accomplished with an expectation maximization algorithm that processes the low-flux data in aggregate, and without any prior knowledge of the object or its orientation. The versatility of the method promises, more generally, to redefine what measurement scenarios can provide useful signal in the high-noise regime.
△ Less
Submitted 15 March, 2012;
originally announced March 2012.
-
Three dimensional structure from intensity correlations
Authors:
Veit Elser
Abstract:
We develop the analysis of x-ray intensity correlations from dilute ensembles of identical particles in a number of ways. First, we show that the 3D particle structure can be determined if the particles can be aligned with respect to a single axis having a known angle with respect to the incident beam. Second, we clarify the phase problem in this setting and introduce a data reduction scheme that…
▽ More
We develop the analysis of x-ray intensity correlations from dilute ensembles of identical particles in a number of ways. First, we show that the 3D particle structure can be determined if the particles can be aligned with respect to a single axis having a known angle with respect to the incident beam. Second, we clarify the phase problem in this setting and introduce a data reduction scheme that assesses the integrity of the data even before the particle reconstruction is attempted. Finally, we describe an algorithm that reconstructs intensity and particle density simultaneously, thereby making maximal use of the available constraints.
△ Less
Submitted 20 July, 2011;
originally announced July 2011.
-
Uniqueness transition in noisy phase retrieval
Authors:
Veit Elser,
Stefan Eisebitt
Abstract:
Previous criteria for the feasibility of reconstructing phase information from intensity measurements, both in x-ray crystallography and more recently in coherent x-ray imaging, have been based on the Maxwell constraint counting principle. We propose a new criterion, based on Shannon's mutual information, that is better suited for noisy data or contrast that has strong priors not well modeled by c…
▽ More
Previous criteria for the feasibility of reconstructing phase information from intensity measurements, both in x-ray crystallography and more recently in coherent x-ray imaging, have been based on the Maxwell constraint counting principle. We propose a new criterion, based on Shannon's mutual information, that is better suited for noisy data or contrast that has strong priors not well modeled by continuous variables. A natural application is magnetic domain imaging, where the criterion for uniqueness in the reconstruction takes the form that the number of photons, per pixel of contrast in the image, exceeds a certain minimum. Detailed studies of a simple model show that the uniqueness transition is of the type exhibited by spin glasses.
△ Less
Submitted 23 September, 2010;
originally announced September 2010.
-
Recovering magnetization distributions from their noisy diffraction data
Authors:
Ne-Te Duane Loh,
Stefan Eisebitt,
Samuel Flewett,
Veit Elser
Abstract:
We study, using simulated experiments inspired by thin film magnetic domain patterns, the feasibility of phase retrieval in X-ray diffractive imaging in the presence of intrinsic charge scattering given only photon-shot-noise limited diffraction data. We detail a reconstruction algorithm to recover the sample's magnetization distribution under such conditions, and compare its performance with that…
▽ More
We study, using simulated experiments inspired by thin film magnetic domain patterns, the feasibility of phase retrieval in X-ray diffractive imaging in the presence of intrinsic charge scattering given only photon-shot-noise limited diffraction data. We detail a reconstruction algorithm to recover the sample's magnetization distribution under such conditions, and compare its performance with that of Fourier transform holography. Concerning the design of future experiments, we also chart out the reconstruction limits of diffractive imaging when photon- shot-noise and the intensity of charge scattering noise are independently varied. This work is directly relevant to the time-resolved imaging of magnetic dynamics using coherent and ultrafast radiation from X-ray free electron lasers and also to broader classes of diffractive imaging experiments which suffer noisy data, missing data or both.
△ Less
Submitted 6 August, 2010;
originally announced August 2010.
-
Strategies for processing diffraction data from randomly oriented particles
Authors:
Veit Elser
Abstract:
This note compares the single-shot and intensity cross-correlation proposals for x-ray imaging of randomly oriented particles and shows very directly that the latter will usually not be feasible even when the former is.
This note compares the single-shot and intensity cross-correlation proposals for x-ray imaging of randomly oriented particles and shows very directly that the latter will usually not be feasible even when the former is.
△ Less
Submitted 21 July, 2010;
originally announced July 2010.
-
Cryptotomography: reconstructing 3D Fourier intensities from randomly oriented single-shot diffraction patterns
Authors:
N. D. Loh,
M. Bogan,
V. Elser,
A. Barty,
S. Boutet,
S. Bajt,
J. Hajdu,
T. Ekeberg,
F. R. N. C. Maia,
J. Schulz,
M. M. Seibert,
B. Iwan,
N. Timneanu,
S. Marchesini,
I. Schlichting,
R. L. Shoeman,
L. Lomb,
M. Frank,
M. Liang,
H. N. Chapman
Abstract:
We reconstructed the 3D Fourier intensity distribution of mono-disperse prolate nano-particles using single-shot 2D coherent diffraction patterns collected at DESY's FLASH facility when a bright, coherent, ultrafast X-ray pulse intercepted individual particles of random, unmeasured orientations. This first experimental demonstration of cryptotomography extended the Expansion-Maximization-Compres…
▽ More
We reconstructed the 3D Fourier intensity distribution of mono-disperse prolate nano-particles using single-shot 2D coherent diffraction patterns collected at DESY's FLASH facility when a bright, coherent, ultrafast X-ray pulse intercepted individual particles of random, unmeasured orientations. This first experimental demonstration of cryptotomography extended the Expansion-Maximization-Compression (EMC) framework to accommodate unmeasured fluctuations in photon fluence and loss of data due to saturation or background scatter. This work is an important step towards realizing single-shot diffraction imaging of single biomolecules.
△ Less
Submitted 5 April, 2010; v1 submitted 3 March, 2010;
originally announced March 2010.
-
A reconstruction algorithm for single-particle diffraction imaging experiments
Authors:
Duane Ne-Te Loh,
Veit Elser
Abstract:
We introduce the EMC algorithm for reconstructing a particle's 3D diffraction intensity from very many photon shot-noise limited 2D measurements, when the particle orientation in each measurement is unknown. The algorithm combines a maximization step (M) of the intensity's likelihood function, with expansion (E) and compression (C) steps that map the 3D intensity model to a redundant tomographic…
▽ More
We introduce the EMC algorithm for reconstructing a particle's 3D diffraction intensity from very many photon shot-noise limited 2D measurements, when the particle orientation in each measurement is unknown. The algorithm combines a maximization step (M) of the intensity's likelihood function, with expansion (E) and compression (C) steps that map the 3D intensity model to a redundant tomographic representation and back again. After a few iterations of the EMC update rule, the reconstructed intensity is given to the difference-map algorithm for reconstruction of the particle contrast. We demonstrate reconstructions with simulated data and investigate the effects of particle complexity, number of measurements, and the number of photons per measurement. The relatively transparent scaling behavior of our algorithm provides a first estimate of the data processing resources required for future single-particle imaging experiments.
△ Less
Submitted 14 July, 2009; v1 submitted 16 April, 2009;
originally announced April 2009.
-
Divide and concur: A general approach to constraint satisfaction
Authors:
Simon Gravel,
Veit Elser
Abstract:
Many difficult computational problems involve the simultaneous satisfaction of multiple constraints which are individually easy to satisfy. Such problems occur in diffractive imaging, protein folding, constrained optimization (e.g., spin glasses), and satisfiability testing. We present a simple geometric framework to express and solve such problems and apply it to two benchmarks. In the first ap…
▽ More
Many difficult computational problems involve the simultaneous satisfaction of multiple constraints which are individually easy to satisfy. Such problems occur in diffractive imaging, protein folding, constrained optimization (e.g., spin glasses), and satisfiability testing. We present a simple geometric framework to express and solve such problems and apply it to two benchmarks. In the first application (3SAT, a boolean satisfaction problem), the resulting method exhibits similar performance scaling as a leading context-specific algorithm (walksat). In the second application (sphere packing), the method allowed us to find improved solutions to some old and well-studied optimization problems. Based upon its simplicity and observed efficiency, we argue that this framework provides a competitive alternative to stochastic methods such as simulated annealing.
△ Less
Submitted 31 December, 2007;
originally announced January 2008.
-
Noise limits in the assembly of diffraction data
Authors:
Veit Elser
Abstract:
We obtain an information theoretic criterion for the feasibility of assembling diffraction signals from noisy tomographs when the positions of the tomographs within the signal are unknown. For shot-noise limited data, the minimum number of detected photons per tomograph for successful assembly is much smaller than previously believed necessary, growing only logarithmically with the number of res…
▽ More
We obtain an information theoretic criterion for the feasibility of assembling diffraction signals from noisy tomographs when the positions of the tomographs within the signal are unknown. For shot-noise limited data, the minimum number of detected photons per tomograph for successful assembly is much smaller than previously believed necessary, growing only logarithmically with the number of resolution elements of the diffracting object. We also demonstrate assembly up to the information theoretic limit with a constraint-based algorithm.
△ Less
Submitted 24 September, 2007;
originally announced September 2007.
-
The sloppy model universality class and the Vandermonde matrix
Authors:
Joshua J. Waterfall,
Fergal P. Casey,
Ryan N. Gutenkunst,
Kevin S. Brown,
Christopher R. Myers,
Piet W. Brouwer,
Veit Elser,
James P. Sethna
Abstract:
In a variety of contexts, physicists study complex, nonlinear models with many unknown or tunable parameters to explain experimental data. We explain why such systems so often are sloppy; the system behavior depends only on a few `stiff' combinations of the parameters and is unchanged as other `sloppy' parameter combinations vary by orders of magnitude. We contrast examples of sloppy models (fro…
▽ More
In a variety of contexts, physicists study complex, nonlinear models with many unknown or tunable parameters to explain experimental data. We explain why such systems so often are sloppy; the system behavior depends only on a few `stiff' combinations of the parameters and is unchanged as other `sloppy' parameter combinations vary by orders of magnitude. We contrast examples of sloppy models (from systems biology, variational quantum Monte Carlo, and common data fitting) with systems which are not sloppy (multidimensional linear regression, random matrix ensembles). We observe that the eigenvalue spectra for the sensitivity of sloppy models have a striking, characteristic form, with a density of logarithms of eigenvalues which is roughly constant over a large range. We suggest that the common features of sloppy models indicate that they may belong to a common universality class. In particular, we motivate focusing on a Vandermonde ensemble of multiparameter nonlinear models and show in one limit that they exhibit the universal features of sloppy models.
△ Less
Submitted 6 October, 2006; v1 submitted 15 May, 2006;
originally announced May 2006.
-
Reconstruction of an object from its symmetry-averaged diffraction pattern
Authors:
Veit Elser
Abstract:
By suitably generalizing the Fourier constraint projection in the difference map phasing algorithm, an object can be reconstructed from its diffraction pattern even when the latter has been incoherently averaged over a discrete group of symmetries. This resolves an ambiguity in the recent proposal for aligning molecules by means of their anisotropic dielectric interaction with an intense light f…
▽ More
By suitably generalizing the Fourier constraint projection in the difference map phasing algorithm, an object can be reconstructed from its diffraction pattern even when the latter has been incoherently averaged over a discrete group of symmetries. This resolves an ambiguity in the recent proposal for aligning molecules by means of their anisotropic dielectric interaction with an intense light field. The algorithm is demonstrated with simulated data in two and three dimensions.
△ Less
Submitted 24 May, 2005;
originally announced May 2005.
-
Dynamics of immersed molecules in superfluids
Authors:
Michael J. Quist,
Veit Elser
Abstract:
The dynamics of a molecule immersed in a superfluid medium are considered. Results are derived using a classical hydrodynamic approach followed by canonical quantization. The classical model, a rigid body immersed in incompressible fluid, permits a thorough analysis; its effective Hamiltonian generalizes the usual rigid-rotor Hamiltonian. In contrast to the free rigid rotor, the immersed body is…
▽ More
The dynamics of a molecule immersed in a superfluid medium are considered. Results are derived using a classical hydrodynamic approach followed by canonical quantization. The classical model, a rigid body immersed in incompressible fluid, permits a thorough analysis; its effective Hamiltonian generalizes the usual rigid-rotor Hamiltonian. In contrast to the free rigid rotor, the immersed body is shown to have chaotic dynamics. Quantization of the classical model leads to new and experimentally verifiable features. It is shown, for instance, that chiral molecules can behave as "quantum propellers": the rotational-translational coupling induced by the superfluid leads to a nonzero linear momentum in the ground state. Hydrogen peroxide is a strong candidate for experimental detection of this effect. The signature is a characteristic splitting of rotational absorption lines. The 1_{01} --> 1_{10} line in hydrogen peroxide, for example, is predicted to split into three lines separated by as much as 0.01 cm^{-1}, which is about the experimental linewidth.
△ Less
Submitted 19 November, 2002; v1 submitted 4 December, 2001;
originally announced December 2001.