Julia sets for the super-Newton method,Cauchy's method,and Halley's method

We study numerically and dynamically three cubically convergent iterative root-finding algorithms, namely Cauchy's method, the super-Newton method, and Halley's method. Using the concept of a universal Julia set (motivated by the results of McMullen), we establish that these algorithms converge when applied to any quadratic with distinct roots. We give examples showing the existence of attracting periodic orbits not associated to a root for the super-Newton method and Halley's method applied to cubic polynomials. We include computer plots showing the dynamic structure for each algorithm applied to a variety of polynomials. (c) 2001 American Institute of Physics.     

基于单曲面迭代计算的曲面透镜设计方法

为在平面上实现所需的照明,满足照明系统的要求,需要对光学系统进行二次设计。传统的设计方法往往用规则曲面指定透镜的一个表面,并以此为基础来设计另个一表面,这样会导致其表面形状复杂,不方便加工或装配。对此,提出了一种基于单曲面迭代计算的透镜曲面设计方法,即适当设计指定光线从透镜的一个表面到另一个表面的光路分布,再以此分别迭代计算透镜的各个表面。该设计方法采用两个自由曲面配合成曲面透镜,较以往的设计方法更加灵活、方便,设计得到的曲面镜头较以往更加简单光滑,便于加工和装配。结合一陈列灯镜头设计,给出了基于单曲面迭代计算的曲面透镜设计计算结果。     

SQUID法和Campbell法测量超导材料的研究

本文对比研究了超导材料磁测量中的SQUID法和Campbell法; 并用高压PIT法制备的超导材料MgB₂作为测量样品,用两种方法,测量了超导样品的临界电流密度, 分别得到了样品的J_c-B关系曲线; SQUID法测量样品的外磁场可以达到6 T, 此时材料已经处于失超状态,此方法测得的结果是样品各个小区域结果的平均值, SQUID还可以用来进一步标度材料的钉扎力行为,研究材料磁特性. Campbell法测量只能测量到外磁场强度为0.4 T,外磁场的交流部分的频率可以达到800 Hz, 用这种测量方法得到的是整块样品的电流,由于测量计及材料内部微观结构缺陷等影响电流传输因素, 所测结果小于直流磁化法,但更切近材料实际电流,能用来深入研究材料内部结构差别对材料电性能的影响.     

Relation between heat kernel method and scattering spectral method

Taking into account the available accelerator and astrophysical constraints, the mass of the lightest neutral Higgs boson h in the minimal supersymmetric extension of the Standard Model with universal soft supersymmetry-breaking masses (CMSSM) has been estimated to lie between 114 and ??130?GeV. Recent data from ATLAS and CMS hint that m _h ??125?GeV, though m _h ??119?GeV may still be a possibility. Here we study the consequences for the parameters of the CMSSM and direct dark matter detection if the Higgs hint is confirmed, focusing on the strips in the (m _1/2,m ₀) planes for different tan?? and A ₀ where the relic density of the lightest neutralino ?? falls within the range of the cosmological cold dark matter density allowed by WMAP and other experiments. We find that if m _h ??125?GeV focus-point strips would be disfavoured, as would the low-tan?? ${\tilde{\tau}}$ ?C?? and ${\tilde{t}}_{1} $ ?C?? coannihilation strips, whereas the ${\tilde{\tau}}$ ?C?? coannihilation strip at large tan?? and A ₀>0 would be favoured, together with its extension to a funnel where rapid annihilation via direct-channel H/A poles dominates. On the other hand, if m _h ??119?GeV more options would be open. We give parameterisations of WMAP strips with large tan?? and fixed A ₀/m ₀>0 that include portions compatible with m _h =125?GeV, and present predictions for spin-independent elastic dark matter scattering along these strips. These are generally low for models compatible with m _h =125?GeV, whereas the XENON100 experiment already excludes some portions of strips where m _h is smaller.     

“0.618”法、Fibonacci法、抛物线法搜寻复摆周期极值点

在实验室条件下,用"0.618"法、Fibonacci法和抛物线法来搜寻钟摆复摆周期极值点位置.通过实验发现,Fibonacci法在搜索精度上比"0.618"法高,但在搜寻速度上不占优势,而抛物线法在搜寻速度和精度上比"0.618"法和Fibonacci法都好.     

Hybrid encoding method by assembling the magic-matrix scrambling method and the binary encoding method in image hiding

The magic-matrix scrambling method and the binary encoding method are combined to form a hybrid encoding method for hiding digital covert images. For this hybrid encoding method, a covert image is encoded into a host image to form an overt image. First, the magic-matrix scrambling method is used to rearrange all the pixels of the covert image by using a specified magic matrix modified from a magic square to form a scrambled matrix. Then, all the pixels of the scrambled matrix are denoted by binary data. Finally, the binary data are encoded into the host image to form the overt image. The pixels of the overt image contain nine groups of codes used for decoding the covert image, i.e. identification codes, covert-image dimension codes, scrambling-time codes, magic-square dimension codes, corner codes, shifting codes, arrangement codes, graylevel codes, and information codes. The overt image and the host image look almost the same for eyes. Furthermore, the covert image can be decoded directly from the overt image without using the host image. The most important feature is that the decoded covert image is identical to the original covert image, i.e. there is no distortion in the decoding work.     

A new method for orbit determination: Unit vector method

In this paper, a method of orbit determination is presented according to the principle of unit vector method (UVM). The model and arithmetic are improved and it not only suits initial orbit determination with short arc data, it also suits orbit improvement with data longer. It is also suitable for orbit of any eccentricity and any inclination. It omits most partial derivatives of all the elements which must be calculated in classical differential orbit improvement (DOI), so, it is more efficient than DOI, and the accuracy of orbit determination and convergence of algorithm are also improved appreciably.     

Quantum Noether method

We present a general method for constructing perturbative quantum field theories with global symmetries. We start from a free non-interacting quantum field theory with given global symmetries and we determine all perturbative quantum deformations assuming the construction is not obstructed by anomalies. The method is established within the causal Bogoliubov-Shirkov-Epstein-Glaser approach to perturbative quantum field theory (which leads directly to a finite perturbative series and does not rely on an intermediate regularization). Our construction can be regarded as a direct implementation of Noether's method at the quantum level. We illustrate the method by constructing the pure Yang-Mills theory (where the relevant global symmetry is BRST symmetry), and the N = 1 supersymmetric model of Wess and Zumino. The whole construction is done before the so-called adiabatic limit is taken. Thus, all considerations regarding symmetry, unitarity and anomalies are well defined even for massless theories.     

Basic function method

Science China Physics, Mechanics & Astronomy - A new numerical method-basic function method is proposed. This method can directly discrete differential operators on unstructured grids. By using...     

Divergence-free WKB method

A new semiclassical approach to the linear and nonlinear one-dimensional Schr?dinger equation is presented. For both equations our zeroth-order solutions include nonperturbative quantum corrections to the WKB solution and the Thomas-Fermi solution, thereby allowing us to make uniformly converging perturbative expansions of the wave functions. Our method leads to a new quantization condition that yields exact eigenenergies for the harmonic-oscillator and Morse potentials.     

Complementary method to locate atomic coordinates by combined searching method of structure-sensitive indexes based on bond valence method

Bond valence method illustrates the relation between valence and length of a particular bond type. This theory has been used to predict structure information, but the effect is very limited. In this paper, two indexes, i.e., global instability index(GII) and bond strain index(BSI), are adopted as a judgment of a search-match program for prediction. The results show that with GII and BSI combined as judgment, the predicted atom positions are very close to real ones. The mechanism and validity of this searching program are also discussed. The GII BSI distribution contour map reveals that the predicted function is a reflection of exponential feature of bond valence formula. This combined searching method may be integrated with other structure-determination method, and may be helpful in refining and testifying light atom positions.     

Numerical-experimental hybrid method for stress separation in digital gradient sensing method

A numerical-experimental hybrid method for the stress separation in the digital gradient sensing (DGS) method is proposed in this study. In the proposed hybrid method, boundary conditions for a local finite element model, that is, nodal force along boundaries are inversely determined from experimental values obtained by the digital gradient sensing method. The hybrid method follows two stages. In stage 1, the DGS method measures the Cartesian stress gradient components directly and, subsequently, the sum in Cartesian stresses at all interesting points on the surface; stress sum are used to compute the unknown boundary conditions for the local model. In stage 2, the individual stress components are calculated by the direct finite element method using the computed boundary conditions from stage 1. The effectiveness is demonstrated by applying the proposed method to a stress concentration problem involving concentrated load acting on an edge of a large planar sheet. The individual stress components thus determined are summed and compared with analytical stress sum, confirming the effectiveness and accuracy of the proposed technique.     

Robust wrapping-free phase retrieval method based on weighted least squares method

Phase unwrapping is one of the most challenging processes in many profilometry techniques. To sidestep the phase unwrapping process, Perciante et al. (2015) proposed a wrapping-free method to retrieve the phase based on the direct integration of the spatial derivatives of the fringe patterns. However, this method is only applicable to objects with phase continuity, so it may fail to handle fringe patterns containing complicated singularities such as noise, shadows, shears and surface discontinuities. In the light of this problems, a robust wrapping-free phase retrieval method is proposed that is based on the combined use of Perciante's method and the weighted least squares method. Two partial derivatives of the desired phase are obtained from the fringe patterns, while the carrier is eliminated using the direct phase difference method. The phase singularities are determined using a derivative variance correlation map (DVCM), and the weighting coefficient is obtained from the binary mask of the reverse DVCM. Simulations and experiments are conducted to prove the validity of the proposed method. The results are analyzed and compared with those of Perciante's method demonstrating that in addition to maintaining the advantage of sidestepping the phase unwrapping process, the proposed method is available for measuring objects with some types of singularities sources.     

Using the idea of the cluster-variation method in the crystal-quasiparticle method

A generalization of the method of theoretically investigating the correlation of melt properties with its structure is proposed. A relationship between the new approximation and the cluster-variation method is established. Specific examples are analyzed.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 9, pp. 107–111, September, 1981.     

A new interface tracking method: The polygonal area mapping method

We present a new method, the polygonal area mapping (PAM) method, for tracking a non-diffusive, immiscible material interface between two materials in two-dimensional incompressible flows. This method represents material areas explicitly as piecewise polygons, traces characteristic points on polygon boundaries along pathlines and calculates new material areas inside interface cells via polygon-clippings in a discrete manner. The new method has very little spatial numerical diffusion and tracks the interface singularities naturally and accurately. In addition to high accuracy, the PAM method can be directly used on either a structured rectangular mesh or an unstructured mesh without any modifications. The mass conservation is enforced by heuristic algorithms adjusting the volume of material polygons. The results from a set of widely used benchmark tests show that the PAM method is superior to existing volume-of-fluid (VOF) methods.     

数字散斑相关方法亚像素求解的一种混合方法

亚像素求解是数字散斑相关方法(DSCM)中最重要的技术。亚像素求解算法的精度和效率直接影响DSCM的精度和计算速度。目前广泛应用的两种亚像素求解方法在处理高质量散斑图时都存在一定的缺陷:基于灰度插值的散斑场亚像素恢复方法,其计算精度较高,但计算消耗太大;基于相关系数分布的拟合方法,其计算消耗小,但计算精度较低。在分析两种亚像素求解方法各自优缺点的基础上,提出了一种结合两种算法优势的混合算法。混合算法在保证亚像素求解精度的基础上,可以大大降低整个算法的计算消耗。     

An alternative method of using combined-space method: high index surfaces

This paper presents an alternative method for calculating I(E)-spectra for a high index crystal surface using the SATLEED (symmetrized automated tensor-LEED) program. Two test structures, Cu{2 1 1} and Cu{2 1 1}-c(2 × 2)-Cs, have been analyzed. The SATLEED results are compared with results from earlier published CHANGE analyses for the same systems. Relaxations for the first three interlayer spacings of Cu{2 1 1} are −0.09 ± 0.01, −0.05 ± 0.02 and +0.05 ± 0.02 Å, respectively. The corresponding values for the Cs adsorption system are −0.09 ± 0.01, −0.04 ± 0.01 and +0.06 ± 0.01 Å. The adsorption site and effective radius of Cs determined by SATLEED are slightly different from the earlier CHANGE study. Computational and theoretical issues related to the use of these two programs are also discussed.