Phase-only filtering based on sliced orthogonal nonlinear generalized decomposition

The sliced orthogonal nonlinear generalized (SONG) decomposition and correlation have been demonstrated to be powerful tools in digital image processing and promising for nonlinear optical information processing. In this paper, we propose an optical phase-only filtering system based on SONG decomposition (PBS), in which the phase-only filtering of the target and the input scene binary slices are performed separately by pairs and then added together. We numerically show that the PBS has extremely high and sharp output correlation peak compared with other optical correlators. Furthermore, such a SONG decomposition based phase-only filtering naturally inherits the nature of SONG decomposition, which has strong robustness to additive Gaussian noise and substitutive noise, and also the high light efficiency of phase-only filtering. We demonstrate that the PBS may serve as an optimized optical correlation scheme, which is promising in nonlinear optical pattern recognition.     

一种分数阶非线性系统稳定性判定定理的问题及分析

分数阶非线性系统稳定性理论的研究对于分数阶混沌系统同步控制的应用具有重要价值,将分数阶非线性系统稳定性判断转化为相应整数阶非线性系统稳定性判断的探讨很有意义.通过实例表明了:对于时变系数矩阵,如果整数阶系统稳定,其对应的阶次小于1的分数阶系统也稳定的判定定理是错误的,并分析了问题产生的原因.     

Two model equations with a second degree logarithmic nonlinearity and their Gaussian solutions

In the paper, we try to study the mechanism of the existence of Gaussian waves in high degree logarithmic nonlinear wave motions. We first construct two model equations which include the high order dispersion and a second degree logarithmic nonlinearity. And then we prove that the Gaussian waves can exist for high degree logarithmic nonlinear wave equations if the balance between the dispersion and logarithmic nonlinearity is kept. Our mathematical tool is the logarithmic trial equation method.     

Order-by-disorder in classical oscillator systems

We consider classical nonlinear oscillators on hexagonal lattices. When the coupling between the elements is repulsive, we observe coexisting states, each one with its own basin of attraction. These states differ by their degree of synchronization and by patterns of phase-locked motion. When disorder is introduced into the system by additive or multiplicative Gaussian noise, we observe a non-monotonic dependence of the degree of order in the system as a function of the noise intensity: intervals of noise intensity with low synchronization between the oscillators alternate with intervals where more oscillators are synchronized. In the latter case, noise induces a higher degree of order in the sense of a larger number of nearly coinciding phases. This order-by-disorder effect is reminiscent to the analogous phenomenon known from spin systems. Surprisingly, this non-monotonic evolution of the degree of order is found not only for a single interval of intermediate noise strength, but repeatedly as a function of increasing noise intensity. We observe noise-driven migration of oscillator phases in a rough potential landscape.     

Effects of harmonic phase on nonlinear surface acoustic waves in the (111) plane of cubic crystals

Spectral evolution equations are used to perform numerical studies of nonlinear surface acoustic waves in the (111) plane of several nonpiezoelectric cubic crystals. Nonlinearity matrix elements which describe the coupling of harmonic interactions are used to characterize velocity waveform distortion. In contrast to isotropic solids and the (001) plane of cubic crystals, the nonlinearity matrix elements usually cannot be written in a real-valued form. As a result, the harmonic components are not necessarily in phase, and dramatic variations in waveforms and propagation curves can be observed. Simulations are performed for initially monofrequency surface waves. In some directions the waveforms distort in a manner similar to nonlinear Rayleigh waves, while in other directions the velocity waveforms distort asymmetrically and the formation of shocks and cusped peaks is less distinct. In some cases, oscillations occur near the shocks and peaks because of phase differences between harmonics. A mathematical transformation based on the phase of the matrix elements is shown to provide a reasonable approximation of asymmetric waveform distortion in cases where the matrix elements have similar phase.     

Hartree-fock equations for atoms with two open shells

We have formulated the Hartree-Fock equations for multielectron systems with two open shells (the Huzinaga method) in terms of the density matrix in the LCAO approximation. In order to solve the Hartree-Fock equations, in the algebraic approximation we have obtained expressions for the derivatives of the energy with respect to the density matrix elements and the nonlinear atomic orbital parameters (the orbital exponents). We discuss the question of calculating the open shell parameters (the vector coupling coefficients) in the configurations s¹p^N and s¹d^N within the Huzinaga method. We have calculated the energy for a series of atoms with two open shells in these configurations. Using rather narrow basis sets of Slater-type atomic orbitals), we have obtained energy values close to the results of a numerical solution of the Hartree-Fock equations with sufficiently high accuracy of the virial ratio. __________ Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 74, No. 2, pp. 145–152, March–April, 2007     

Semi-automatic sparse preconditioners for high-order finite element methods on non-uniform meshes

High-order finite elements often have a higher accuracy per degree of freedom than the classical low-order finite elements. However, in the context of implicit time-stepping methods, high-order finite elements present challenges to the construction of efficient simulations due to the high cost of inverting the denser finite element matrix. There are many cases where simulations are limited by the memory required to store the matrix and/or the algorithmic components of the linear solver. We are particularly interested in preconditioned Krylov methods for linear systems generated by discretization of elliptic partial differential equations with high-order finite elements. Using a preconditioner like Algebraic Multigrid can be costly in terms of memory due to the need to store matrix information at the various levels. We present a novel method for defining a preconditioner for systems generated by high-order finite elements that is based on a much sparser system than the original high-order finite element system. We investigate the performance for non-uniform meshes on a cube and a cubed sphere mesh, showing that the sparser preconditioner is more efficient and uses significantly less memory. Finally, we explore new methods to construct the sparse preconditioner and examine their effectiveness for non-uniform meshes. We compare results to a direct use of Algebraic Multigrid as a preconditioner and to a two-level additive Schwarz method.     

How random are matrix elements of the nuclear shell model Hamiltonian?

In this paper we study the general behavior of matrix elements of the nuclear shell model Hamiltonian.We find that nonzero off-diagonal elements exhibit a regular pattern,if one sorts the diagonal matrix elements from smaller to larger values.The correlation between eigenvalues and diagonal matrix elements for the shell model Hamiltonian is more remarkable than that for random matrices with the same distribution unless the dimension is small.     

Pathological tremors: Deterministic chaos or nonlinear stochastic oscillators?

Pathological tremors exhibit a nonlinear oscillation that is not strictly periodic. We investigate whether the deviation from periodicity is due to nonlinear deterministic chaotic dynamics or due to nonlinear stochastic dynamics. To do so, we apply various methods from linear and nonlinear time series analysis to tremor time series. The results of the different methods suggest that the considered types of pathological tremors represent nonlinear stochastic second order processes. Finally, we evaluate whether two earlier proposed features capturing nonlinear effects in the time series allow for a discrimination between two pathological forms of tremor for a much larger sample of time series than previously investigated. (c) 2000 American Institute of Physics.     

Nonlinear moment method for the isotropic Boltzmann equation and invariance of the collision integral

A new approach is proposed for the development of a nonlinear moment method of solving the Boltzmann equation. This approach is based on the principle of invariance of the collision integral with respect to the choice of basis functions. Sonine polynomials with a Maxwellian weighting function are taken as these basis functions for the velocity-isotropic Boltzmann equation. It is shown that for arbitrary interaction cross sections the matrix elements corresponding to the moments of the nonlinear collision integral are not independent but are coupled by simple recurrence formulas by means of which all the nonlinear matrix elements are expressed in terms of linear ones. As a result, a highly efficient numerical scheme is constructed for calculating the nonlinear matrix elements. The proposed approach opens up prospects for calculating relaxation processes at high velocities and also for solving more complex kinetic problems. Zh. Tekh. Fiz. 69, 22–29 (June 1999)     

An investigation into the linear and nonlinear correlation of two music walk sequences

In order to investigate the features of a multiple-part musical score which enhance its appeal to the listener’s ear, this study performs a robust analysis of the correlation between two musical sequences. In the proposed approach, a series of notes are extracted from seven well-known classical pieces of music and are converted into one-variable “music walks”. The linear correlation between pairs of music walks is assessed using the conventional linear correlation coefficient, while the nonlinear correlation is examined using the mutual information concept. The results show that even though two musical time walks may exhibit virtually no linear correlation, they invariably have a certain degree of nonlinear correlation. In other words, to truly understand the correlation between two musical sequences, it is necessary to consider not only the linear correlation between them, but also the nonlinear correlation. In addition, it is shown that the normalized mutual information coefficient between musical sequences has a relatively low value and varies significantly over the course of the musical score. Thus, it can be inferred that the appeal of a musical score stems at least in part from significant variations in both the melody and the rhythm of the constituent parts such that the overall score has a rich and unpredictable property.     

Control of higher order antibunching

Based on general physical arguments we derive a new general criterion for higher order antibunching. We show that the degree of antibunching increases with the order of antibunched photon states. We use the criterion for the case of the interaction of an intense electromagnetic field with a (k-1)-th order nonlinear medium. For a third order nonlinear and nonabsorbing medium with inversion symmetry, we show that the degree of antibunching can be tuned by changing the phase of the external field. PACS 42.50.Ar; 42.50.ct; 42.50.Dv     

用于心脏电活动成像的空间滤波器输出噪声抑制方法

利用人体表面测量的心脏磁场数据无创成像心脏电活动,需要解决的关键问题是提高其重建分布电流源偶极矩强度的分辨率.本文在最小方差波束成形(MVB)方法的基础上,提出了一种可抑制空间滤波器输出噪声功率增益(SONG)的波束成形方法,目的是通过构造一种新的滤波器权矩阵,约束空间滤波器的输出噪声功率增益,提高重建分布电流源偶极矩强度的分辨率,即分布电流源空间谱估计的源分辨能力,从而增强心脏电活动磁成像的分辨率.文中给出了电流源重建的理论分析和仿真结果;比较了该方法与MVB方法的差别;并给出了两个健康人36通道心脏磁场数据的电活动成像.结果表明,SONG方法分辨电流源的能力较强,能够观察到心脏电磁场信号R峰时刻健康人的心室电活动较强,心房电活动较弱等特征.     

Complete deterministic linear optics Bell state analysis

We show how hyperentanglement allows us to deterministically distinguish between all four polarization Bell states of two photons. In this proof-of-principle experiment, we employ the intrinsic time-energy correlation of photon pairs generated with high temporal definition in addition to the polarization entanglement obtained from parametric down-conversion. For the identification, no nonlinear optical elements or auxiliary photons are needed. The new possibilities this complete Bell measurement offers are demonstrated by realizing an optimal dense coding protocol.     

Isotropic matrix elements of the collision integral for the Boltzmann equation

We have proposed an algorithm for constructing matrix elements of the collision integral for the nonlinear Boltzmann equation isotropic in velocities. These matrix elements have been used to start the recurrent procedure for calculating matrix elements of the velocity-nonisotropic collision integral described in our previous publication. In addition, isotropic matrix elements are of independent interest for calculating isotropic relaxation in a number of physical kinetics problems. It has been shown that the coefficients of expansion of isotropic matrix elements in Ω integrals are connected by the recurrent relations that make it possible to construct the procedure of their sequential determination.     

A phase-synchronization and random-matrix based approach to multichannel time-series analysis with application to epilepsy

We present a general method to analyze multichannel time series that are becoming increasingly common in many areas of science and engineering. Of particular interest is the degree of synchrony among various channels, motivated by the recognition that characterization of synchrony in a system consisting of many interacting components can provide insights into its fundamental dynamics. Often such a system is complex, high-dimensional, nonlinear, nonstationary, and noisy, rendering unlikely complete synchronization in which the dynamical variables from individual components approach each other asymptotically. Nonetheless, a weaker type of synchrony that lasts for a finite amount of time, namely, phase synchronization, can be expected. Our idea is to calculate the average phase-synchronization times from all available pairs of channels and then to construct a matrix. Due to nonlinearity and stochasticity, the matrix is effectively random. Moreover, since the diagonal elements of the matrix can be arbitrarily large, the matrix can be singular. To overcome this difficulty, we develop a random-matrix based criterion for proper choosing of the diagonal matrix elements. Monitoring of the eigenvalues and the determinant provides a powerful way to assess changes in synchrony. The method is tested using a prototype nonstationary noisy dynamical system, electroencephalogram (scalp) data from absence seizures for which enhanced cortico-thalamic synchrony is presumed, and electrocorticogram (intracranial) data from subjects having partial seizures with secondary generalization for which enhanced local synchrony is similarly presumed.     

Frequency dependence of phase-synchronization time in nonlinear dynamical systems

It has been found recently that the averaged phase-synchronization time between the input and the output signals of a nonlinear dynamical system can exhibit an extremely high sensitivity to variations in the noise level. In real-world signal-processing applications, sensitivity to frequency variations may be of considerable interest. Here we investigate the dependence of the averaged phase-synchronization time on frequency of the input signal. Our finding is that, for typical nonlinear oscillator systems, there can be a frequency regime where the time exhibits significant sensitivity to frequency variations. We obtain an analytic formula to quantify the frequency dependence, provide numerical support, and present experimental evidence from a simple nonlinear circuit system.     

Multiple parton scattering in nuclei: twist-four nuclear matrix elements and off-forward parton distributions

Multiple parton scatterings inside a large nucleus generally involve higher-twist nuclear parton matrix elements. The gluon bremsstrahlung induced by multiple scattering depends not only on direct parton matrix elements but also on momentum-crossed ones, due to the Landau–Pomeranchuk–Migdal interference effect. We show that both types of twist-four nuclear parton matrix elements can be factorized approximately into the product of twist-two nucleon matrix elements in the limit of extremely large nuclei, A→∞, as assumed in previous studies. Due to the correlative nature of the twist-four matrix elements under consideration, it is actually the off-forward parton distributions that appear naturally in this decomposition, rather than the ordinary diagonal distributions probed in deeply-inelastic scattering. However, we argue that the difference between these two distribution classes is small in certain kinematic regimes. In these regions, the twist-four nuclear parton matrix elements are evaluated numerically and compared to the factorized form for different nuclear sizes within a schematic model of the two-nucleon correlation function. The nuclear size dependence is found to be A^4/3 in the limit of large A, as expected. We find that the factorization is reasonably good when the momentum fraction carried by the gluon field is moderate. The deviation can be more than a factor of 2, however, for small gluon momentum fractions, where the gluon distribution is very large.     

Asymptotics of collision integral matrix elements in the isotropic case

The method of nonlinear moments, when used to solve the Boltzmann equation, necessitates the calculation of collision integral matrix elements. The matrix elements are hard to calculate numerically, especially at large indices. The asymptotics of the matrix elements are constructed. In terms of the model of pseudopower particle interaction, a formula free of summation is derived. This makes it possible to find the asymptotic behavior of linear and nonlinear elements when two indices are large. For an arbitrary interaction cross section, asymptotic expansions of linear and nonlinear matrix elements in one index are obtained. For Maxwellian molecules, asymptotic formulas are derived for three large indices.     

Statistical analysis and model determination for thermal desorption spectra: Nitrogen on tungsten

Thermal desorption spectra of nitrogen chemisorbed on tungsten have been examined in order to discriminate between two simple alternative kinetic models for the desorption process. Both these models consider the β₂ desorption as a second order process, while for the β₁ peak they assume first order and second order kinetics, respectively. The maximum likelihood method has been used for fitting the spectra and the resulting kinetic parameters have been discussed through the analysis of the correlation matrix. The strict correlation which exists between activation energies and frequency factors prevents effective discrimination and should diminish confidence in too simple and rigid models.