Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “APPL”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

Overexpression of a repressor MdMYB15L negatively regulates anthocyanin and cold tolerance in red-fleshed callus

Highlights: • MdMYB15L inhibit the expression of MdCBF2 by binding its promoter. • MdMYB15L interact with MdbHLH33 via its bHLH binding sequence. • MdMYB15L weaken MdbHLH33-induced anthocyanin accumulation by interacting with MdbHLH33. The cold-induced metabolic pathway and anthocyanin biosynthesis play important roles in plant growth. In this study, we identified a bHLH binding motif in the MdMYB15L protein using protein sequence analyses. Yeast two-hybrid and pull-down assays showed that MdMYB15L could interact with MdbHLH33. Overexpressing MdMYB15L in red-fleshed callus inhibited the expression of MdCBF2 and resulted in reduced cold tolerance but did not affect anthocyanin levels. Chip-PCR and EMSA analysis showed that MdMYB15L could bind the type II cis-acting element found in the MdCBF2 promoter. Overexpressing MdMYB15L in red-fleshed callus overexpressing MdbHLH33 also reduced cold tolerance and reduced MdbHLH33-induced anthocyanin biosynthesis. Knocking out the bHLH binding sequence of MdMYB15L (LBSMdMYB15L) prevented LBSMdMYB15L from interacting with MdbHLH33. Overexpressing LBSMdMYB15L in red-fleshed callus overexpressing MdbHLH33 also reduced cold tolerance and reduced MdbHLH33-induced anthocyanin biosynthesis. Together, these results suggested that an apple repressor MdMYB15L might play a key role in the cold signaling and anthocyanin metabolic pathways.

60 APPLIED LIFE SCIENCES↗

Scalable DPG multigrid solver for Helmholtz problems: A study on convergence

This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from $\mathscr{O}$(10 7 ) degrees of freedom (DOFs) to $\mathscr{O}$(10 9 ) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. Finally, the paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast.

97 MATHEMATICS AND COMPUTING↗

Coupled uptake and conversion of C 12 H 26 and NO on Pd/SSZ-13: Experiments and modeling

The uptake of NO in the presence of C 12 H 26 (dodecane) and H 2 O over a Pd/SSZ-13 washcoated Passive NOx Adsorber (PNA) monolith is reported. When a co-feed containing C 12 H 26 and NO is supplied to an unsaturated (with C 12 H 26 ) sample, the NO uptake is unaffected but during the subsequent temperature ramp the release of trapped NO is delayed from 175°C to over 220°C. The release delay is beneficial for PNA performance as the primary NO x aftertreatment technology, Selective Catalytic Reduction (SCR), is not operated below 200°C. However, pre-saturation of C 12 H 26 followed by the same NO and C 12 H 26 co-feed results in a decrease in the NO uptake compared to the NO-only feed. We conjecture that C 12 H 26 pre-adsorbed on the exterior surface of the sample blocks NO access to the pores, decreasing the number of available sites for NO uptake. Oxidation of C 12 H 26 leads to the generation of partial oxidation product CO at lower temperatures (<250°C) and deep oxidation product CO 2 at higher temperatures. Carbon monoxide binds strongly to Pd sites with NO and can delay NO release. A Pd/SSZ-13 washcoated monolith model developed in an earlier study [Ambast et al., Appl. Catal. B. Environmental (2021)] is upgraded to include C 12 H 26 storage, release, and conversion. Through a systematic combination of judicious experiments, model tuning, and validation, we provide evidence for the underlying NO uptake and release mechanisms in the presence of C 12 H 26 .

42 ENGINEERING↗

Self-consistent clustering analysis: an efficient multiscale scheme for inelastic heterogeneous materials

A predictive computational theory is shown for modeling complex, hierarchical materials ranging from metal alloys to polymer nanocomposites. The theory can capture complex mechanisms such as plasticity and failure that span across multiple length scales. This general multiscale material modeling theory relies on sound principles of mathematics and mechanics, and a cutting-edge reduced order modeling method named self-consistent clustering analysis (SCA) [Zeliang Liu, M.A. Bessa, Wing Kam Liu, “Self-consistent clustering analysis: An efficient multi-scale scheme for inelastic heterogeneous materials,” Comput. Methods Appl. Mech. Engrg. 306 (2016) 319–341]. SCA reduces by several orders of magnitude the computational cost of micromechanical and concurrent multiscale simulations, while retaining the microstructure information. This remarkable increase in efficiency is achieved with a data-driven clustering method. Computationally expensive operations are performed in the so-called offline stage, where degrees of freedom (DOFs) are agglomerated into clusters. The interaction tensor of these clusters is computed. In the online or predictive stage, the Lippmann-Schwinger integral equation is solved cluster-wise using a self-consistent scheme to ensure solution accuracy and avoid path dependence. To construct a concurrent multiscale model, this scheme is applied at each material point in a macroscale structure, replacing a conventional constitutive model with the average response computed from the microscale model using just the SCA online stage. A regularized damage theory is incorporated in the microscale that avoids the mesh and RVE size dependence that commonly plagues microscale damage calculations. The SCA method is illustrated with two cases: a carbon fiber reinforced polymer (CFRP) structure with the concurrent multiscale model and an application to fatigue prediction for additively manufactured metals. For the CFRP problem, a speed up estimated to be about 43,000 is achieved by using the SCA method, as opposed to FE2, enabling the solution of an otherwise computationally intractable problem. The second example uses a crystal plasticity constitutive law and computes the fatigue potency of extrinsic microscale features such as voids. This shows that local stress and strain are capture sufficiently well by SCA. This model has been incorporated in a process-structure-properties prediction framework for process design in additive manufacturing.

Liu, Z.↗

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗

An Efficient Algorithm for Automatic Structure Optimization in X-ray Standing-Wave Experiments

We present that x-ray standing-wave photoemission experiments involving multilayered samples are emerging as unique probes of the buried interfaces that are ubiquitous in current device and materials research. Such data require for their analysis a structure optimization process comparing experiment to theory that is not straightforward. In this work, we present a new computer program for optimizing the analysis of standing-wave data, called SWOPT, that automates this trial-and-error optimization process. The program includes an algorithm that has been developed for computationally expensive problems: so-called black-box simulation optimizations. It also includes a more efficient version of the Yang X-ray Optics Program (YXRO) [Yang, S.-H., Gray, A.X., Kaiser, A.M., Mun, B.S., Sell, B.C., Kortright, J.B., Fadley, C.S., J. Appl. Phys. 113, 1 (2013)] which is about an order of magnitude faster than the original version. Human interaction is not required during optimization. We tested our optimization algorithm on real and hypothetical problems and show that it finds better solutions significantly faster than a random search approach. The total optimization time ranges, depending on the sample structure, from minutes to a few hours on a modern laptop computer, and can be up to 100x faster than a corresponding manual optimization. Finally, these speeds make the SWOPT program a valuable tool for realtime analyses of data during synchrotron experiments.

36 MATERIALS SCIENCE↗

Energy efficiency in industrial drying: A hybrid ultrasonic system with a novel dynamic optimization framework

Drying processes are among the most energy-consuming operations in industrial and manufacturing settings, demanding strategic selection, design, and control for enhanced efficiency. Advancing drying technologies is critical for improving sustainability, lowering energy use, reducing carbon emissions, and minimizing waste. This study explores two innovative strategies aimed at transforming drying processes into sustainable, low-carbon systems by reducing energy consumption, minimizing waste, and maintaining a strong emphasis on preserving product quality. The first strategy showcases a sub-pilot scale hybrid ultrasonic-convective dryer for agrifood products. This technology, powered by electricity (process electrification), integrates non-thermal ultrasonic dehydration with convective heating and is presented as a sustainable and energy-efficient solution that enhances eco-friendly practices. The second strategy involves introducing and implementing a novel, multiobjective, mixed integer dynamic optimization technique to determine the optimal time-dependent process parameter values for the drying operation. This optimization technique yields operating conditions that are piecewise constant in time aiming to maximize the energy efficiency of the hybrid ultrasonic-convective dryer while ensuring strict adherence to product quality constraints. By adopting the hybrid ultrasonic-convective dryer, a notable 35% improvement in energy efficiency was achieved compared to conventional hot-air drying systems for drying apple slices. The proposed optimization framework further enhanced energy efficiency by nearly 14% over the most efficient process on the identical testbed, under static operating conditions. The reported enhancements have been experimentally validated. Regarding drying time (thereby improving production yield), the developed hybrid ultrasonic-convective dryer demonstrates as much as a 41% reduction in total processing time, which is further optimized by an additional 10% using our proposed optimization framework. The research outcomes have profound implications for the design and operation of drying systems, encompassing crucial aspects such as process electrification, cost-effectiveness, energy savings, time efficiency, product yield, product quality, and process automation.

Dynamic optimization↗

A nominally second-order cell-centered Lagrangian scheme for simulating elastic–plastic flows on two-dimensional unstructured grids

In this paper, we describe a cell-centered Lagrangian scheme devoted to the numerical simulation of solid dynamics on two-dimensional unstructured grids in planar geometry. This numerical method, utilizes the classical elastic-perfectly plastic material model initially proposed by Wilkins [M.L. Wilkins, Calculation of elastic–plastic flow, Meth. Comput. Phys. (1964)]. In this model, the Cauchy stress tensor is decomposed into the sum of its deviatoric part and the thermodynamic pressure which is defined by means of an equation of state. Regarding the deviatoric stress, its time evolution is governed by a classical constitutive law for isotropic material. The plasticity model employs the von Mises yield criterion and is implemented by means of the radial return algorithm. The numerical scheme relies on a finite volume cell-centered method wherein numerical fluxes are expressed in terms of sub-cell force. The generic form of the sub-cell force is obtained by requiring the scheme to satisfy a semi-discrete dissipation inequality. Sub-cell force and nodal velocity to move the grid are computed consistently with cell volume variation by means of a node-centered solver, which results from total energy conservation. The nominally second-order extension is achieved by developing a two-dimensional extension in the Lagrangian framework of the Generalized Riemann Problem methodology, introduced by Ben-Artzi and Falcovitz [M. Ben-Artzi, J. Falcovitz, Generalized Riemann Problems in Computational Fluid Dynamics, Cambridge Monogr. Appl. Comput. Math. (2003)]. Finally, the robustness and the accuracy of the numerical scheme are assessed through the computation of several test cases.

97 MATHEMATICS AND COMPUTING↗

Theory of magnetic enhancement in strontium hexaferrite through Zn-Sn pair substitution

We study the site occupancy and magnetic properties of Zn-Sn substituted M-type Sr-hexaferrite SrFe12-x(Zn0.5Sn0.5)(x)O-19 with x=1 using first-principles total-energy calculations. We find that in the lowest-energy configuration Zn2+ and Sn4+ ions preferentially occupy the 4f(1) and 4f(2) sites, respectively, in contrast to the model previously suggested by Ghasemi et al. [J. Appl. Phys, 107, 09A734 (2010)], where Zn2+ and Sn4+ ions occupy the 2b and 4f(2) sites. Density-functional theory calculations show that our model has a lower total energy by more than 0.2 eV per unit cell compared to Ghasemi's model. More importantly, the latter does not show an increase in saturation magnetization (M-s) compared to the pure M-type Sr-hexaferrite, in disagreement with the experiment. On the other hand, our model correctly predicts a rapid increase in M-s as well as a decrease in magnetic anisotropy compared to the pure M-type Sr-hexaferrite, consistent with experimental measurements. (c) 2013 Elsevier B.V. All rights reserved.

Liyanage, LSI↗

A three dimensional field formulation, and isogeometric solutions to point and line defects using Toupin's theory of gradient elasticity at finite strains

In this study, we present a field formulation for defects that draws from the classical representation of the cores as force dipoles. We write these dipoles as singular distributions. Exploiting the key insight that the variational setting is the only appropriate one for the theory of distributions, we arrive at universally applicable weak forms for defects in nonlinear elasticity. Remarkably, the standard, Galerkin finite element method yields numerical solutions for the elastic fields of defects that, when parameterized suitably, match very well with classical, linearized elasticity solutions. The true potential of our approach, however, lies in its easy extension to generate solutions to elastic fields of defects in the regime of nonlinear elasticity, and even more notably for Toupin's theory of gradient elasticity at finite strains (Toupin Arch. Ration. Mech. Anal., 11 (1962) 385). In computing these solutions we adopt recent numerical work on an isogeometric analytic framework that enabled the first three-dimensional solutions to general boundary value problems of Toupin's theory (Rudraraju et al. Comput. Methods Appl. Mech. Eng., 278 (2014) 705). We first present exhaustive solutions to point defects, edge and screw dislocations, and a study on the energetics of interacting dislocations. Then, to demonstrate the generality and potential of our treatment, we apply it to other complex dislocation configurations, including loops and low-angle grain boundaries.

36 MATERIALS SCIENCE↗

Third-order elastic constants of diamond determined from experimental data

The pressure derivatives of the second-order elastic constants (SOECs) of diamond were determined by analyzing previous sound velocity measurements under hydrostatic stress [McSkimin and Andreatch, J. Appl. Phys. 43, 294 (1972)]. Furthermore, our analysis corrects an error in the previously reported results.We present a complete and corrected set of third-order elastic constants (TOECs) using the corrected pressure derivatives, together with published data for the nonlinear elastic response of shock compressed diamond [Lang and Gupta, Phys. Rev. Lett. 106, 125502 (2011)] and it differs significantly from TOECs published previously.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward reducing uncertainty quantification costs in DEM models of particulate flow: Testing simple, sensitivity-based, forward uncertainty propagation techniques

The performance of two conceptually-simple uncertainty quantification techniques are tested against the rigorous nested-loop sampling technique of Roy and Oberkampf (Comput Methods Appl Mech Eng, 200: 2131–2144, 2011) (herein called full-sampling) using two very small-scale DEM-based models of particulate flow (one gas-solid flow and one granular flow). The first simplified forward uncertainty propagation technique, reduced-sampling, uses a sensitivity analysis to eliminate uncertain inputs that have little impact on the model output prior to nested-loop sampling. The second technique, boundary-sampling, uses a sensitivity analysis to inform the selection of two bounding cases for each key model output. In conclusion, the uncertainties in the model outputs obtained via the reduced- and boundary-sampling methods agree well with those from full-sampling for both the gas-solid and granular flow models while yielding computational savings of 65–75% (reduced sampling) and 94–97% (boundary sampling).

42 ENGINEERING↗

Part 2: Forensic attribution profiling of Russian VX in food using liquid chromatography-mass spectrometry

This work is part two of a three-part series in this issue of a Sweden-United States collaborative effort towards the understanding of the chemical attribution signatures of Russian VX (VR) in synthesized samples and complex food matrices. We describe the sourcing of VR present in food based on chemical analysis of attribution signatures by liquid chromatography-tandem mass spectrometry (LC-MS/MS) combined with multivariate data analysis. Analytical data was acquired from seven different foods spiked with VR batches that were synthesized via six different routes in two separate laboratories. The synthesis products were spiked at a lethal dose into seven food matrices: water, orange juice, apple purée, baby food, pea purée, liquid eggs and hot dog. After acetonitrile sample extraction, the samples were analyzed by LC-MS/MS operated in MRM mode. A multivariate statistical calibration model was built on the chemical attribution profiles from 118 VR spiked food samples. Using the model, an external test-set of the six synthesis routes employed for VR production was correctly identified with no observable major impact of the food matrices to the classification. The overall performance of the statistical models was found to be exceptional (94%) for the test set samples retrospectively classified to their synthesis routes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An efficient method for unfolding kinetic pressure driven VISAR data

Velocity Interferometer System for Any Reflector (VISAR) [Barker and Hollenbach, J. Appl. Phys. 43, 4669 (1972)] is a well-known diagnostic that is employed on many shock physics and pulsed-power experiments. With the VISAR diagnostic, the velocity on the surface of any metal flyer can be found. For most experiments employing VISAR, either a kinetic pressure [Grady, Mech. Mater. 29, 181 (1998)] or a magnetic pressure [Lemke et al., Intl J. Impact Eng. 38, 480 (2011)] drives the motion of the flyer. Moreover, reliable prediction of the time-dependent pressure is often a critical component to understanding the physics of these experiments. Although VISAR can provide a precise measurement of a flyer’s surface velocity, the real challenge of this diagnostic implementation is using this velocity to unfold the time-dependent pressure. The purpose of this study is to elucidate a new method for quickly and reliably unfolding VISAR data.

42 ENGINEERING↗

Recent laser upgrades at Sandia’s Z-backlighter facility in order to accommodate new requirements for magnetized liner inertial fusion on the Z-machine

The Z-backlighter laser facility primarily consists of two high energy, high-power laser systems. Z-Beamlet laser (ZBL) (Rambo et al., Appl. Opt. 44, 2421 (2005)) is a multi-kJ-class, nanosecond laser operating at 1054 nm which is frequency doubled to 527 nm in order to provide x-ray backlighting of high energy density events on the Z-machine. Z-Petawatt (ZPW) (Schwarz et al., J. Phys.: Conf. Ser. 112, 032020 (2008)) is a petawatt-class system operating at 1054 nm delivering up to 500 J in 500 fs for backlighting and various short-pulse laser experiments (see also Figure 10 for a facility overview). With the development of the magnetized liner inertial fusion (MagLIF) concept on the Z-machine, the primary backlighting missions of ZBL and ZPW have been adjusted accordingly. As a result, we have focused our recent efforts on increasing the output energy of ZBL from 2 to 4 kJ at 527 nm by modifying the fiber front end to now include extra bandwidth (for stimulated Brillouin scattering suppression). The MagLIF concept requires a well-defined/behaved beam for interaction with the pressurized fuel. Hence we have made great efforts to implement an adaptive optics system on ZBL and have explored the use of phase plates. We are also exploring concepts to use ZPW as a backlighter for ZBL driven MagLIF experiments. Alternatively, ZPW could be used as an additional fusion fuel pre-heater or as a temporally flexible high energy pre-pulse. All of these concepts require the ability to operate the ZPW in a nanosecond long-pulse mode, in which the beam can co-propagate with ZBL. Finally, some of the proposed modifications are complete and most of them are well on their way.

47 OTHER INSTRUMENTATION↗

Differential formulation of the viscous history force on a particle for efficient and accurate computation

It is well known that the computation of the Basset-like history force is very demanding in terms of CPU and memory requirements, since it requires the evaluation of a history integral. We use the recent rational theory of Beylkin & Monzón ( Appl. Comput. Harmon. Anal. , vol. 19, 2005, pp. 17–48) to approximate the history kernel in the form of exponential sums to reformulate the viscous history force in a differential form. This theory allows us to approximate the history kernel in terms of exponential sums to any desired order of accuracy. This removes the need for long-time storage of the acceleration histories of the particle and the fluid. The proposed differential form approximation is applied to compute the history force on a spherical particle in a synthetic turbulent flow and a wall-bounded turbulent channel flow. Particles of various diameters are considered, and results obtained using the present technique are in reasonable agreement with those achieved using the full history integral.

Mechanics↗

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase Imaging: A Compressive Sensing Approach

Since Wolfgang Pauli posed the question in 1933, whether the probability densities |Ψ(r)|² (real-space image) and |Ψ(q)|² (reciprocal space image) uniquely determine the wave function Ψ(r) [1], the so called Pauli Problem sparked numerous methods in all fields of microscopy [2, 3]. Reconstructing the complete wave function Ψ(r) = a(r)e-iφ(r) with the amplitude a(r) and the phase φ(r) from the recorded intensity enables the possibility to directly study the electric and magnetic properties of the sample through the phase. In transmission electron microscopy (TEM), electron holography is by far the most established method for phase reconstruction [4]. Requiring a high stability of the microscope, next to the installation of a biprism in the TEM, holography cannot be applied to any microscope straightforwardly. Recently, a phase retrieval approach was proposed using conventional TEM electron diffractive imaging (EDI). Using the SAD aperture as reciprocal-space constraint, a localized sample structure can be reconstructed from its diffraction pattern and a real-space image using the hybrid input-output algorithm [5]. We present an alternative approach using compressive phase-retrieval [6]. Our approach does not require a real-space image. Instead, random complimentary pairs of checkerboard masks are cut into a 200 nm Pt foil covering a conventional TEM aperture (cf. Figure 1). Used as SAD aperture, subsequently diffraction patterns are recorded from the same sample area. Hereby every mask blocks different parts of gold particles on a carbon support (cf. Figure 2). The compressive sensing problem has the following formulation. First, we note that the complex-valued reciprocal-space wave-function is the Fourier transform of the (also complex-valued) real-space wave-function, Ψ(q) = F[Ψ(r)], and subsequently the diffraction pattern image is given by |Ψ(q)|2 = |F[Ψ(r)]|2. We want to find Ψ(r) given a few differently coded diffraction pattern measurements yn = |F[HnΨ(r)]|2, where the matrices Hn encode the mask structure of the aperture. This is a nonlinear inverse problem, but has been shown to be solvable even in the underdetermined case [6]. Since each diffraction pattern yn contains diffraction information from selected regions of the same sample, the differences in each pattern contain local phase information, which can be combined to form a full estimate of the real-space wave-function[7]. References: [1] W. Pauli in “Die allgemeinen Prinzipien der Wellenmechanik“, ed. H Geiger and W Scheel, (Julius Springer, Berlin). [2] A. Tonomura, Rev. Mod. Phys. 59 (1987), p. 639. [3] J. Miao et al, Nature 400 (1999), p. 342. [4] H. Lichte et al, Annu. Rev. Mater. Res. 37 (2007), p. 539. [5] J. Yamasaki et al, Appl. Phys. Lett. 101 (2012), 234105. [6] P Schniter and S Rangan. Signal Proc., IEEE Trans. on. 64(4), (2015), pp. 1043. [7] Supported by the Chemical Imaging, Signature Discovery, and Analytics in Motion initiatives at PNNL. PNNL is operated by Battelle Memorial Inst. for the US DOE; contract DE-AC05-76RL01830.

Schneider, Sebastian↗