Reconstruction of the fine structure of an acoustic scatterer against the distorting influence of its large-scale inhomogeneities

In the ultrasonic diagnostics of small-size neoplasms of biological tissues at the earliest stage of their development, an efficient way to eliminate the distorting influence of high-contrast or large inhomogeneities of the biological medium is to apply the iterative technique. A simple approach is proposed, which makes it possible with only two iteration steps to achieve an efficient focusing of the tomograph array. At the first step, the unknown distribution of the large-scale inhomogeneities of sound velocity and absorption over the scatterer is reconstructed, where the large-scale inhomogeneities are those whose size exceeds several wavelengths. At the second step, the fine structure of the scatterer is reconstructed against the large-scale background, which can be performed with a high accuracy owing to the evaluation of the background at the first step. The possibility of simultaneous reconstruction of the large-scale and fine structures by the noniterative Grinevich-Novikov algorithm is considered as an alternative. This algorithm reconstructs in an explicit form two-dimensional refractive-absorbing acoustic scatterers of almost arbitrary shape and strength. Taking into account the effects of multiple scattering, this algorithm provides resolution of the fine structure almost as good as that achieved in reconstructing the same structure against an undistorting homogeneous background. The results of numerical simulations of both algorithms are presented.     

Imaging tumor angiogenesis by use of combined near-infrared diffusive light and ultrasound

A novel two-step reconstruction scheme using a combined near-infrared and ultrasound technique and its utility in imaging distributions of optical absorption and hemoglobin concentration of breast lesions are demonstrated. In the first-step image reconstruction, the entire tissue volume is segmented based on initial coregistered ultrasound measurements into lesion and background regions. Reconstruction is performed by use of a finer grid for lesion region and a coarse grid for the background tissue. As a result, the total number of voxels with unknown absorption can be maintained on the same order of total measurements, and the matrix with unknown total absorption distribution is appropriately scaled for inversion. In the second step, image reconstruction is refined by optimization of lesion parameters measured from ultrasound images. It is shown that detailed distributions of wavelength-dependent absorption and hemoglobin concentration of breast carcinoma can be obtained with the new reconstruction scheme.     

I. The properties of hot Ca-like fragments from the 40Ca + 40Ca reaction at 35AMeV

The creation of hot Ca-like fragments was investigated in the ⁴⁰Ca + ⁴⁰Ca reaction at 35 AMeV. Using the AMPHORA 4π detector system, the primary projectile-like fragment was reconstructed and its properties were determined. Both primary and secondary distributions are compared with the predictions of a Monte Carlo code describing a heavy-ion collision as a two-step process. Some of the nucleons which are identified as participants in the first step are transferred in the second step to these final states, which correspond on the average to the maximum value of entropy (thermodynamic probability). The model allows for competition between mean-field effects and nucleon-nucleon interactions in the overlap zone of the interacting nuclei. The analysis presented here suggests a thermalized source picture of the decay of the projectile-like fragment. The validity of the reconstruction procedure for projectile-like fragments is discussed. Received: 12 March 2001 / Accepted: 20 June 2001     

A simple method for NMR photography

The application of a multi-frequency weak pulse to a liquid crystal can excite narrow 1H NMR peaks at the applied frequencies. By using two-level amplitude coding, namely setting the amplitudes of some of the harmonics to zero, this method can be used to store up to 1024 bits of binary information in the liquid crystal molecules. When the information is retrieved in the form of a spectrum, which is plotted as an array of 32 or 16 segments, the stack of spectral segments reproduces a 2D input pattern quite well. This technique is called "NMR photography." The original method was a pseudo-2D technique that applies the 1024-frequency pulse in the first step, and reads the signal row-by-row in the second step. The present improvement involves subtracting two spectra obtained with the same 1024-frequency pulse but with different durations, so that the method becomes a 1D technique, with a tremendous saving of experimental time. Several examples are given to illustrate the results.     

Image processing pipeline for synchrotron‐radiation‐based tomographic microscopy

With synchrotron‐radiation‐based tomographic microscopy, three‐dimensional structures down to the micrometer level can be visualized. Tomographic data sets typically consist of 1000 to 1500 projections of 1024 × 1024 to 2048 × 2048 pixels and are acquired in 5–15 min. A processing pipeline has been developed to handle this large amount of data efficiently and to reconstruct the tomographic volume within a few minutes after the end of a scan. Just a few seconds after the raw data have been acquired, a selection of reconstructed slices is accessible through a web interface for preview and to fine tune the reconstruction parameters. The same interface allows initiation and control of the reconstruction process on the computer cluster. By integrating all programs and tools, required for tomographic reconstruction into the pipeline, the necessary user interaction is reduced to a minimum. The modularity of the pipeline allows functionality for new scan protocols to be added, such as an extended field of view, or new physical signals such as phase‐contrast or dark‐field imaging etc.     

Two-dimensional phased arrays for surgery: Movement of a single focus

Numerical simulation and comparative analysis of acoustic fields generated by two-dimensional phased arrays designed for ultrasonic surgery is conducted. The case of movement of a single focus by an array with the surface shaped as a part of a spherical shell with the curvature radius 120 mm is considered. The influence of the number of elements (varying from 64 to 1024), their diameter (from 2.5 to 10 mm), frequency (from 1 to 2 MHz), and the degree of sparseness of the elements at the array surface on the field characteristics is studied. The calculations are performed for arrays with elements positioned at the surface both regularly (in square, annular, or hexagonal patterns) and randomly. Criteria for the evaluation of the “quality” of intensity distributions in the field generated by an array in the case of movement of a single focus are suggested. Of all arrays studied, the best quality of distributions is obtained for an array containing 256 elements of diameter 5 mm randomly positioned at the array surface. The quality of the intensity distributions for arrays consisting of 255, 256, and 1024 elements positioned regularly (in square, annular, and hexagonal patterns) is inferior to the corresponding quality for arrays with randomly positioned elements. The irregularity in elements’ positioning considerably improves the distribution quality by suppressing the secondary intensity peaks in the field generated by the array; or, alternatively, it provides an opportunity to obtain the same distribution quality with a fraction of the number of elements in the array. The effects of the number and shape of elements, errors in phase setting, frequency modulation of signals, and non-uniform distribution of amplitudes over the array surface on the distribution quality are analyzed.     

Generalized parton distributions from hadronic observables: non-zero skewness

We propose a physically motivated parameterization for the unpolarized generalized parton distributions, H and E, valid at both zero and non-zero values of the skewness variable, ζ. Our approach follows a previous detailed study of the ζ=0 case where H and E were determined using constraints from simultaneous fits of the experimental data on both the nucleon elastic form factors and the deep inelastic structure functions in the non-singlet sector. Additional constraints at ζ≠0 are provided by lattice calculations of the higher moments of generalized parton distributions. We illustrate a method for extracting generalized parton distributions from lattice moments based on a reconstruction using sets of orthogonal polynomials. The inclusion in our fit of data on Deeply Virtual Compton Scattering is also discussed. Our method provides a step towards an extraction of generalized distributions based on a global fit of the available data within the given set of constraints.     

Theoretical studies of fast He atom scattering from a W{100} surface

Theoretical calculations of the scattering of fast neutral He atoms from W{100} are presented which are directly compared to the results from recent experiments of Nielsen and Delchar. The experiments which were performed for He atom energies between 150 and 1000 eV, and for incident polar angles between 0° and 65° as measured from the surface normal, displayed peaks in the polar angle distributions at 72°, 56° and ±10°. The results from classical dynamics calculations are employed here to explain the scattering phenomena that give rise to these peaks in the polar distributions. The calculations indicate that the peak at 72° is primarily due to scattering from the first layer W atoms. The peaks at ± 10° and 56° are unusual in that there are a multitude of different collision paths that result in the He atom being scattered into the same final angle. The peaks at ± 10° result from He atoms scattering mainly from the second, third and fourth layers of W atoms. The He atoms are focused on the outward path into the near normal direction by two first and two second layer W atoms. Subsurface scattering is also responsible for the peak at 56°. In this case the channel of first and second layer W atoms that focuses the outgoing He atoms is oriented at 54.7° with respect to the surface normal. It is proposed that slight variations of the experimental data from the calculated values are due to surface reconstruction of W{100} and that a more thorough analysis could reveal the microscopic nature of this structure.     

Damping of high-lying single-particle modes in heavy nuclei

The recent experimental and theoretical results on the damping of high-lying single-particle modes in heavy nuclei are reviewed. In one-nucleon transfer reactions these states manifest themselves as broad “resonance”-like structures superimposed on a large continuum. The advantages and the limitations of the transfer reaction approach will be presented using the results from neutron and proton pick-up and stripping reactions. The problem raised by the subtraction of the underlying background, the assumptions made to describe the reaction process and the method used to extract the strength distributions are presented. The existing empirical systematics is summarized for nuclei ranging from ⁹⁰Zr to ²⁰⁸Pb.The theoretical approaches used to explain the damping of the high-lying single-particle modes are based on the coupling between collective and single-particle degrees of freedom. In a first step the bare single-particle mode is spread over several doorway collective states due to the interaction with surface vibrations. In a second step the doorway states spread their strengths over many other degrees of freedom. These two steps of the damping mechanism are discussed in detail within the framework of the quasiparticle-phonon nuclear model. A large-scale comparison between the measured and calculated average energies, spreading widths and spectroscopic strengths of the high-lying single-particle (hole) states in heavy nuclei is presented. The systematic features of the damping (energy, angular momentum and isotopic dependence) are discussed. Recent advances of the experimental approaches, such as the γ-decay of the high-lying states or the use of heavy-ion transfer reactions at intermediate energies, are outlined.The detailed study of the damping mechanism of high-lying single-particle modes reveals new features and leads us to a new field in nuclear structure: “the spectroscopy of inner and outer subshells”.     

Rough surface characterization using off-axis digital holographic microscopy compensated with self-hologram rotation

In this paper, an off-axis digital holographic microscopy compensated with self-hologram rotation is presented. The process is implemented via subtracting the unwrapped phase maps of the off-axis parabolic hologram and its rotation 180° to eliminate the tilt induced by the angle between the spherical object wave O and the plane reference wave R. Merit of the proposed method is that it can be done without prior knowledge of physical parameters and hence can reconstruct a parabolic hologram of 1024 × 768 pixels within tens of milliseconds since it doesn't require a digital reference wave. The method is applied to characterize rough gold bumps and the obtained results were compared with those extracted from the conventional reconstruction method. The comparison showed that the proposed method can characterize rough surfaces with excellent contrast and in real-time. Merit of the proposed method is that it can be used for monitoring smaller biological cells and micro-fluidic devices.     

Visualization of three-dimensional temperature distributions in a large-scale furnace via regularized reconstruction from radiative energy images: numerical studies

In this paper, the possibility of visualization of three-dimensional (3-D) temperature distributions in large-scale boiler furnaces from radiative energy images captured by multiple charge-coupled device (CCD) cameras mounted around the furnace is studied numerically. For the calculation of the radiative energy image formation, a fast algorithm proposed by the authors for pinhole imaging is used in this paper, which is based on the Monte Carlo method and combined with a concept of angular factor effective for image formation. This algorithm is applicable for the emitting, absorbing, and isotropic scattering medium. For the inversion of the 3-D temperature distributions which is an ill-posed problem, a modified Tikhonov regularization method is improved, where the finite difference regularizer is defined and can be used in 3-D cases, and the optimal regularization parameter is suggested to be selected by using a post-treatment method. For a 3-D unimodal temperature distribution, the numerical simulation results show that the reconstruction errors for the 3-D temperature distribution can be maintained at levels similar to the measurement error and the visualization quality of the temperature distribution is satisfactory. For a kind of bimodal temperature distribution, the reconstruction errors are higher than those for the unimodal distribution, but the bimodal feature of the temperature distribution can also be reproduced clearly.     

Calibration Invariance of the MaxEnt Distribution in the Maximum Entropy Principle

The maximum entropy principle consists of two steps: The first step is to find the distribution which maximizes entropy under given constraints. The second step is to calculate the corresponding thermodynamic quantities. The second part is determined by Lagrange multipliers’ relation to the measurable physical quantities as temperature or Helmholtz free energy/free entropy. We show that for a given MaxEnt distribution, the whole class of entropies and constraints leads to the same distribution but generally different thermodynamics. Two simple classes of transformations that preserve the MaxEnt distributions are studied: The first case is a transform of the entropy to an arbitrary increasing function of that entropy. The second case is the transform of the energetic constraint to a combination of the normalization and energetic constraints. We derive group transformations of the Lagrange multipliers corresponding to these transformations and determine their connections to thermodynamic quantities. For each case, we provide a simple example of this transformation.     

Azimuth-invariant mueller-matrix differentiation of the optical anisotropy of biological tissues

A Mueller-matrix model is proposed for analysis of the optical anisotropy of protein networks of optically thin nondepolarizing layers of biological tissues with allowance for birefringence and dichroism. The model is used to construct algorithms for reconstruction of coordinate distributions of phase shifts and coefficient of linear dichroism. Objective criteria for differentiation of benign and malignant tissues of female genitals are formulated in the framework of the statistical analysis of such distributions. Approaches of evidence-based medicine are used to determine the working characteristics (sensitivity, specificity, and accuracy) of the Mueller-matrix method for the reconstruction of the parameters of optical anisotropy and show its efficiency in the differentiation of benign and malignant tumors.     

Energy transfers in dynamos with small magnetic Prandtl numbers

We perform numerical simulation of dynamo with magnetic Prandtl number Pm = 0.2 on 1024³ grid, and compute the energy fluxes and the shell-to-shell energy transfers. These computations indicate that the magnetic energy growth takes place mainly due to the energy transfers from large-scale velocity field to large-scale magnetic field and that the magnetic energy flux is forward. The steady-state magnetic energy is much smaller than the kinetic energy, rather than equipartition; this is because the magnetic Reynolds number is near the dynamo transition regime. We also contrast our results with those for dynamo with Pm = 20 and decaying dynamo.     

Local Au coverage as driving force for Au induced faceting of vicinal Si(0 0 1): a LEEM and XPEEM study

Adsorption of Au at 850°C on a regular stepped 4° vicinal Si(0 0 1) surface results in a dramatic change of the step morphology: the surface decomposes into areas which are perfectly flat with a (0 0 1) orientation and (1 1 9) facets. Low energy electron microscopy shows the dynamics of the faceting process in real space while X-ray photoemission electron microscopy (XPEEM) allows a spatially resolved determination of the Au coverage at different stages of the faceting process. At a critical Au coverage of ≈1/3 ML (0 0 1) terraces are formed which extend anisotropically along the step edges of the surface. The steps in between the terraces bunch and form step bands in order to conserve the macroscopic miscut of the sample. Driving force for this morphological transformation is a complex (5×3.2) reconstruction formed on the (0 0 1) terraces. XPEEM shows this phase separation also for the Au coverage: on the (0 0 1) terraces the Au coverage is up to 40% higher compared to the step bands. With further increasing Au coverage the width of the Au rich terraces increases while the step bands become steeper. In a second step Au adsorbs on the step bands transforming them into well defined and smooth (1 1 9) facets.     

The mesoscopic and microscopic structural consequences from decomposition and desorption of ultrathin oxide layers on Si(100) studied by scanning tunneling microscopy

The spatially inhomogeneous decomposition and desorption reaction of oxide layers with coverage 1-0.3 monolayers (ML) from a silicon (100) surface has been studied using scanning tunneling microscopy (STM). After desorption, microscopic changes to the (2 × 1) reconstruction produce two variations on the dimer row reconstruction with decreased surface atom density. A (2 × n) vacancy chain reconstruction and a c(4 × 4) incomplete row reconstruction were observed; a structure for the latter is proposed. Both reconstructions are metastable, reforming the (2 × 1) reconstruction upon heating. At greater length scales during desorption from an initial 1.0 ML coverage, the mesoscopic changes to the surface structure include pitting and roughening, with up to a measured 20 fold increase in the edge density as compared to the clean Si(100) surface.

These structural changes suggest a reaction mechanism involving a substantial rearrangement of the substrate silicon. From an initial 1.0 ML oxygen coverage, using measured void size distributions at total desorption levels of 13% and below — before voids have begun to coalesce — the evolution of void sizes during initial desorption can be followed. A mechanism for desorption is proposed in which silicon atoms must diffuse from adjacent clean surface area to the oxide boundary, producing a reactive complex from which SiO is desorbed. Void growth rates derived from two rate limiting cases for this desorption reaction mechanism can be compared to measured results. We show that the measured void area evolution is consistent with a reaction mechanism where the rate limiting step for monolayer desorption is the promotion of a silicon atom in a lattice site to a mobile monomer within the void.     

Imprint of large-scale flows on turbulence

We investigate the locality of interactions in hydrodynamic turbulence using data from a direct numerical simulation on a grid of 1024(3) points; the flow is forced with the Taylor-Green vortex. An inertial range for the energy is obtained in which the flux is constant and the spectrum follows an approximate Kolmogorov law. Nonlinear triadic interactions are dominated by their nonlocal components, involving widely separated scales. The resulting nonlinear transfer itself is local at each scale but the step in the energy cascade is independent of that scale and directly related to the integral scale of the flow. Interactions with large scales represent 20% of the total energy flux. Possible explanations for the deviation from self-similar models, the link between these findings and intermittency, and their consequences for modeling of turbulent flows are briefly discussed.     

中能Ne~(4+)离子诱导的羰基硫分子三体碎裂动力学分析

利用反应显微成像谱仪开展了56 keV/u的Ne~(4+)离子与羰基硫(OCS)气体的交叉碰撞实验,研究了Ne~(4+)离子诱导的OCS~(3+)的碎裂动力学.通过符合探测三个末态离子,重构了OCS~(3+)离子三体碎裂对应的牛顿图和Dalitz图,并明确区分了直接解离和次序解离两种碎裂过程.重构了OCS~(3+)离子解离过程的动能释放(KER)分布,发现其峰值在25 eV处,同时在18 eV处有肩膀结构的存在,其中25 eV左右的峰来源于直接解离过程,18 eV处的肩膀结构来源于次序解离和非次序解离两种过程.通过分析不同能量和不同电荷态下重离子碰撞实验所得到的KER谱,发现微扰强度不是影响态布居的主要因素.OCS~(3+)次序解离中的第二步KER的峰值在6.2 eV处.结合以往的实验结果,我们得出结论:多电离态的分子发生次序碎裂的根源在于二价离子碎片存在亚稳态,而重构得到的第二步KER可以反映亚稳态离子的电子态信息.     

用改进的Fienup迭代算法进行数字全息重建

数字全息成像中往往包括再现像、共轭像和零级项,再现像的分离是常见而又难以彻底解决的问题。提出一种新的全息图重建算法,即利用相干衍射成像（CDI）中的迭代方法处理数字全息图,实现仅有一个实像的再现结果,从而彻底解决该问题。该方法包括两个步骤：首先通过CCD分别测量样品单独存在、样品和参考光同时存在以及参考光的远场衍射分布;然后通过计算机进行迭代运算。由于干涉光的存在,该算法比传统的CDI算法有更快的收敛速度和重建质量。同时进行了数值模拟验证并对生物切片进行了物像重建。