01 February 2013, Volume 22 Issue 3
GENERAL
 Mei symmetry and conservation laws of discrete nonholonomic dynamical systems with regular and irregular lattices Zhao Gang-Ling (赵纲领), Chen Li-Qun (陈立群), Fu Jing-Li (傅景礼), Hong Fang-Yu (洪方昱) Chin. Phys. B, 2013, 22 (3):  030201.  DOI: 10.1088/1674-1056/22/3/030201 Abstract ( 735 )   PDF (240KB) ( 566 )   In this paper, Noether symmetry and Mei symmetry of discrete nonholonomic dynamical systems with regular and the irregular lattices are investigated. Firstly, the equations of motion of discrete nonholonomic systems are introduced on the regular and rregular lattices. Secondly, for cases of the two lattices, based on the invariance of the Hamiltomian functional under the infinitesimal transformation of time and generalized coordinates, we present the quasi-extremal equation, the discrete analogues of Noether identity, Noether theorems, and the Noether conservation laws of the systems. Thirdly, in cases of the two lattices, we study the Mei symmetry in which we give the discrete analogues of the criterion, the theorem, and the conservative laws of Mei symmetry for the systems. Finally, an example is discussed for applications of the results.
 Melting phenomenon in magneto hydro-dynamics steady flow and heat transfer over a moving surfacein the presence of thermal radiation Reda G. Abdel-Rahman, M. M. Khader, Ahmed M. Megahed Chin. Phys. B, 2013, 22 (3):  030202.  DOI: 10.1088/1674-1056/22/3/030202 Abstract ( 771 )   PDF (270KB) ( 763 )   The Lie group method is applied to present an analysis of the magneto hydro-dynamics (MHD) steady laminar flow and the heat transfer from a warm laminar liquid flow to a melting moving surface in the presence of thermal radiation. By using the Lie group method, we have presented the transformation groups for the problem with a part from the scaling group. The application of this method reduces the partial differential equations (PDEs) with their boundary conditions governing the flow and heat transfer to a system of nonlinear ordinary differential equations (ODEs) with appropriate boundary conditions. The resulting nonlinear system of ODEs is solved numerically using an implicit finite difference method (FDM). The local skin-friction coefficients and the local Nusselt numbers for different physical parameters are presented in a table.
 Two new discrete integrable systems Chen Xiao-Hong (陈晓红), Zhang Hong-Qing (张鸿庆) Chin. Phys. B, 2013, 22 (3):  030203.  DOI: 10.1088/1674-1056/22/3/030203 Abstract ( 654 )   PDF (183KB) ( 449 )   In this paper, we focus on the construction of new (1+1)-dimensional discrete integrable systems according to a subalgebra of loop algebra Â1. By designing two new (1+1)-dimensional discrete spectral problems, two new discrete integrable systems are obtained, namely, a 2-field lattice hierarchy and a 3-field lattice hierarchy. When deriving the two new discrete integrable systems, we find the generalized relativistic Toda lattice hierarchy and the generalized modified Toda lattice hierarchy. Moreover, we also obtain the Hamiltonian structures of the two lattice hierarchies by means of the discrete trace identity.
 Stability and Neimark-Sacker bifurcation analysis of food-limited population model with time delay Jiang Xiao-Wei (姜晓伟), Guan Zhi-Hong (关治洪), Zhang Xian-He (张先鹤), Zhang Ding-Xue (张顶学), Liu Feng (刘峰) Chin. Phys. B, 2013, 22 (3):  030204.  DOI: 10.1088/1674-1056/22/3/030204 Abstract ( 804 )   PDF (387KB) ( 808 )   In this paper, a kind of discrete delay food-limited model obtained by Euler method is investigated, where the discrete delay τ is regarded as a parameter. By analyzing the associated characteristic equation, the linear stability of this model is studied. It is shown that Neimark–Sacker bifurcation occurs when τ crosses some critical values. The explicit formulae which determine the stability, direction, and other properties of bifurcating periodic solution are derived by means of the theory of center manifold and normal form. Finally, numerical simulations are performed to verify the analytical results.
 Asymptotic solution of weak nonlinear model for mid-latitude stationary wind field of two-layer barotropic ocean Lin Wan-Tao (林万涛), Zhang Yu (张宇), Mo Jia-Qi (莫嘉琪) Chin. Phys. B, 2013, 22 (3):  030205.  DOI: 10.1088/1674-1056/22/3/030205 Abstract ( 783 )   PDF (463KB) ( 403 )   The weak nonlinear model of two-layer barotropic ocean with Rayleigh dissipation is built. The analytic asymptotic solution is derived in mid-latitude stationary wind field, and the physical meaning of the corresponding problem is discussed.
 Stability analysis and control synthesis of uncertain Roesser-type discrete-time two-dimensional systems Wang Jia (王佳), Hui Guo-Tao (会国涛), Xie Xiang-Peng (解相朋) Chin. Phys. B, 2013, 22 (3):  030206.  DOI: 10.1088/1674-1056/22/3/030206 Abstract ( 784 )   PDF (1000KB) ( 988 )   We study the stability analysis and control synthesis of uncertain discrete-time two-dimensional (2D) systems. The mathematical model of the discrete-time 2D system is established upon the well-known Roesser model, and the uncertainty phenomenon, which appears typically in practical environments, is modeled by a convex bounded (polytope type) uncertain domain. Then, the stability analysis and control synthesis of uncertain discrete-time 2D systems are developed by applying the Lyapunov stability theory. In the processes of stability analysis and control synthesis, the obtained stability/stabilzaition conditions become less conservative by applying some novel relaxed techniques. Moreover, the obtained results are formulated in the form of linear matrix inequalities, which can be easily solved via the standard numerical software. Finally, numerical examples are given to demonstrate the effectiveness of the obtained results.
 Dynamical analysis of a sexually transmitted disease model on complex networks Yuan Xin-Peng (原新鹏), Xue Ya-Kui (薛亚奎), Liu Mao-Xing (刘茂省) Chin. Phys. B, 2013, 22 (3):  030207.  DOI: 10.1088/1674-1056/22/3/030207 Abstract ( 656 )   PDF (3715KB) ( 665 )   In this paper, a sexually transmitted disease model is proposed on complex networks, where contacts between humans are treated as a scale-free social network. There are three groups in our model, which are dangerous male, non-dangerous male, and female. By mathematical analysis, we obtain the basic reproduction number for the existence of endemic equilibrium and study the effects of various immunization schemes about different groups. Furthermore, numerical simulations are undertook to reach and verify more conclusions.
 A new complex variable meshless method for advection–diffusion problems Wang Jian-Fei (王健菲), Cheng Yu-Min (程玉民) Chin. Phys. B, 2013, 22 (3):  030208.  DOI: 10.1088/1674-1056/22/3/030208 Abstract ( 736 )   PDF (317KB) ( 633 )   In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, an improved complex variable meshless method (ICVMM) for two-dimensional advection–diffusion problems is developed. The equivalent functional of two-dimensional advection–diffusion problems is formed, the variation method is used to obtain the equation system, and the penalty method is employed to impose the essential boundary conditions. The difference method for two-point boundary value problems is used to obtain the discrete equations. Then the corresponding formulas of the ICVMM for advection–diffusion problems are presented. Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper. It is shown that the ICVMM is very effective for advection–diffusion problems, and has good convergent character, accuracy, and computational efficiency.
 Multisymplectic implicit and explicit methods for Klein–Gordon–Schrödinger equations Cai Jia-Xiang (蔡加祥), Yang Bin (杨斌), Liang Hua (梁华) Chin. Phys. B, 2013, 22 (3):  030209.  DOI: 10.1088/1674-1056/22/3/030209 Abstract ( 695 )   PDF (1830KB) ( 505 )   We propose multisymplectic implicit and explicit Fourier pseudospectral methods for the Klein–Gordon–Schrödinger equations. We prove that the implicit method satisfies the charge conservation law exactly. Both methods provide accurate solutions in long-time computations and simulate the soliton collision well. Numerical results show the abilities of the two methods in preserving charge, energy, and momentum conservation laws.
 Modified Burgers equation by local discontinuous Galerkin method Zhang Rong-Pei (张荣培), Yu Xi-Jun (蔚喜军), Zhao Guo-Zhong (赵国忠) Chin. Phys. B, 2013, 22 (3):  030210.  DOI: 10.1088/1674-1056/22/3/030210 Abstract ( 920 )   PDF (279KB) ( 1240 )   In this paper, we present the local discontinuous Galerkin method for solving Burgers’ equation and the modified Burgers’ equation. We describe the algorithm formulation and practical implementation of the local discontinuous Galerkin method in detail. The method is applied to the solution of the one-dimensional viscous Burgers’ equation and two forms of the modified Burgers’ equation. The numerical results indicate that the method is very accurate and efficient.
 Complex Maxwell's equations A. I. Arbab Chin. Phys. B, 2013, 22 (3):  030301.  DOI: 10.1088/1674-1056/22/3/030301 Abstract ( 981 )   PDF (218KB) ( 5165 )   A unified complex model of Maxwell's equations is presented. The wave nature of electromagnetic field vector is related to the temporal and spatial distributions and the circulation of charge and current densities. A new vacuum solution is obtained. A new transformation under which Maxwell's equations are invariant is proposed. This transformation extends the ordinary gauge transformation to include charge-current besides scalar-vector potential. An electric dipole moment is found to be related to magnetic charges. The Dirac's quantization is found to determine an uncertainty relation expressing the indeterminacy of electric and magnetic charges. We generalize Maxwell's equations to include longitudinal waves. A formal analogy between this formulation and Dirac's equation is discussed.
 Exact solutions of Dirac equation with Pöschl–Teller double-ring-shaped Coulomb potential via Nikiforov–Uvarov method E. Maghsoodi, H. Hassanabadi, S. Zarrinkamar Chin. Phys. B, 2013, 22 (3):  030302.  DOI: 10.1088/1674-1056/22/3/030302 Abstract ( 783 )   PDF (196KB) ( 646 )   Exact analytical solutions of the Dirac equation are reported for the Pöschl–Teller double-ring-shaped Coulomb potential. The radial, polar, and azimuthal parts of the Dirac equation are solved by using the Nikiforov–Uvarov method, and exact bound state energy eigenvalues and the corresponding two-component spinor wavefunctions are reported.
 Relativistic symmetry of position-dependent mass particle in Coulomb field including tensor interaction M. Eshghi, M. Hamzavi, S. M. Ikhdair Chin. Phys. B, 2013, 22 (3):  030303.  DOI: 10.1088/1674-1056/22/3/030303 Abstract ( 677 )   PDF (230KB) ( 562 )   The spatially-dependent mass Dirac equation is solved exactly for attractive scalar and repulsive vector Coulomb potentials including a tensor interaction under the spin and pseudospin symmetric limit. Closed forms of the energy eigenvalue equation and wave functions are obtained for arbitrary spin-orbit quantum number κ. Some numerical results are given too. The effect of the tensor interaction on the bound states is presented. It is shown that the tensor interaction removes the degeneracy between two states in the spin doublets. We also investigate the effects of the spatially-dependent mass on the bound states under the conditions of the spin symmetric limit in the absence of tensor interaction.
 Comparison between non-Markovian dynamics with and without rotating wave approximation Tang Ning (唐宁), Xu Tian-Tian (徐甜甜), Zeng Hao-Sheng (曾浩生) Chin. Phys. B, 2013, 22 (3):  030304.  DOI: 10.1088/1674-1056/22/3/030304 Abstract ( 720 )   PDF (621KB) ( 619 )   In the limit of weak coupling between the system and its reservoir, we derive the time-convolutionless (TCL) non-Markovian master equation for a two-level system interacting with a zero-temperature structured environment with no rotating wave approximation (NRWA). By comparing with the dynamics with RWA, we demonstrate the impact of the RWA on the system dynamics, as well as the effects of non-Markovianity on the preservation of atomic coherence, squeezing, and entanglement.
 Efficient three-step entanglement concentration for an arbitrary four-photon cluster state Si Bin (司斌), Su Shi-Lei (苏石磊), Sun Li-Li (孙立莉), Cheng Liu-Yong (程留永), Wang Hong-Fu (王洪福), Zhang Shou (张寿) Chin. Phys. B, 2013, 22 (3):  030305.  DOI: 10.1088/1674-1056/22/3/030305 Abstract ( 706 )   PDF (266KB) ( 1015 )   We propose an entanglement concentration protocol (ECP) to concentrate an arbitrary partially-entangled four-photon cluster state. As a pioneering three-step entanglement concentration scheme, our protocol only needs single-photon resource to assist the concentration in each step, which makes this protocol more economical. With the help of the linear optical elements and weak cross-Kerr nonlinearity, one can obtain a maximally-entangled cluster state via local operations and classical communication. Moreover, the protocol can be iterated to obtain a higher success probability and is feasible under the current experimental conditions.
 Squeezing entangled state of two particles with unequal masses Yang Yang (杨阳), Fan Hong-Yi (范洪义) Chin. Phys. B, 2013, 22 (3):  030306.  DOI: 10.1088/1674-1056/22/3/030306 Abstract ( 748 )   PDF (176KB) ( 415 )   For two unequal-mass particles, we construct the entangled state representation and then derive the corresponding squeezing operator. This squeezing operator has a natural realization in the entangled state representation, which exhibits the intrinsic relation between squeezing and quantum entanglement. This squeezing operator involves both two-mode squeezing and the direct product of two single-mode squeezings. The maximum squeezing occurs when the two particles possess equal mass. When two particles' mass difference becomes large, the component of the two single-mode squeezings becomes dominant.
 New approach for deriving the exact time evolution of density operator for diffusive anharmonic oscillator and its Wigner distribution function Meng Xiang-Guo (孟祥国), Wang Ji-Suo (王继锁), Liang Bao-Long (梁宝龙) Chin. Phys. B, 2013, 22 (3):  030307.  DOI: 10.1088/1674-1056/22/3/030307 Abstract ( 879 )   PDF (239KB) ( 479 )   Using the thermal entangled state representation, we solve the master equation of a diffusive anharmonic oscillator (AHO) to obtain the exact time evolution formula for the density operator in the infinitive operator-sum representation. We present a new evolution formula of the Wigner function (WF) for any initial state of the diffusive AHO by converting the calculation of the WF to an overlap between two pure states in an enlarged Fock space. It is found that this formula brings us much convenience to investigate the WF's evolution of any known initial state. As applications, this formula is used to obtain the evolution of the WF for a coherent state and the evolution of the photon-number distribution of the diffusive AHO.
 Induced modification of geometric phase of a qubit coupled to an XY spin chain by the Dzyaloshinsky–Moriya interaction Zhang Ai-Ping (张爱萍), Li Fu-Li (李福利) Chin. Phys. B, 2013, 22 (3):  030308.  DOI: 10.1088/1674-1056/22/3/030308 Abstract ( 842 )   PDF (331KB) ( 531 )   We consider a qubit symmetrically and transversely coupled to an XY spin chain with Dzyaloshinsky–Moriya (DM) interaction in the presence of a transverse magnetic field. An analytical expression for the geometric phase of the qubit is obtained in the weak coupling limit. We find that the modification of the geometrical phase induced by the spin chain environment is greatly enhanced by the DM interaction in the weak coupling limit around the quantum phase transition point of the spin chain.
 Electromagnetically induced transparency of single Λ-type three-level atom in high-finesse optical cavity Sun Yan-Fen (孙燕芬), Tan Lei (谭磊), Xu Yan (徐岩) Chin. Phys. B, 2013, 22 (3):  030309.  DOI: 10.1088/1674-1056/22/3/030309 Abstract ( 1003 )   PDF (497KB) ( 1344 )   We study the features of the electromagnetically induced transparency (EIT) in a single Λ-type three-level atom placed in a high finesse cavity under the action of a coupling laser and a probe laser. Our calculations show that three transparency windows appear when the pump strength is big enough. It can be explained by the residual pump in the cavity resulting mostly in the energy splitting. Level |3〉is split into four slightly different energy levels. An interference takes place between excitation pathways. Furthermore, it is also shown that the frequencies of the EIT windows can be tuned by changing the coupling field detuning Δ2 and the reflection profile is very sensitive to the cavity field detuning Δc.
RAPID COMMUNICATION
 Average position in quantum walks with a U(2) coin Hot! Li Min (李敏), Zhang Yong-Sheng (张永生), Guo Guang-Can (郭光灿) Chin. Phys. B, 2013, 22 (3):  030310.  DOI: 10.1088/1674-1056/22/3/030310 Abstract ( 753 )   PDF (266KB) ( 389 )   We investigated discrete-time quantum walks with an arbitary unitary coin. Here we discover that the average position〈x〉=max(〈x〉)sin(α+γ), while the initial state is 1/√2(|0L〉+ i|0R〉). We verify the result, and obtain some symmetry properties of quantum walks with a U(2) coin with |0L〉 and |0R〉 as the initial state.
GENERAL
 Nonclassicality and decoherence of coherent superposition operation of photon subtraction and photon addition on squeezed state Xu Li-Juan (徐莉娟), Tan Guo-Bin (谭国斌), Ma Shan-Jun (马善钧), Guo Qin (郭琴) Chin. Phys. B, 2013, 22 (3):  030311.  DOI: 10.1088/1674-1056/22/3/030311 Abstract ( 843 )   PDF (13850KB) ( 439 )   The statistical properties of m-coherent superposition operation (μa+νa+)m on the single-mode squeezed vacuum state (M-SSVS) and its decoherence in a thermal environment have been studied. Converting the M-SSVS to a squeezed Hermite polynomial excitation state, we obtain a compact expression for the normalization factor of M-SSVS, which is the Legendre polynomial of the squeezing parameter. We also derive the explicit expression of Wigner function (WF) of M-SSVS, and find the negative region of WF in phase space. The decoherence effect on this state is then discussed by deriving the time evolution of the WF. Using the negativity of WF, the loss of nonclassicality has been discussed.
 Maximum confidence measurements via probabilistic quantum cloning Zhang Wen-Hai (张文海), Yu Long-Bao (余龙宝), Cao Zhuo-Liang (曹卓良), Ye Liu (叶柳) Chin. Phys. B, 2013, 22 (3):  030312.  DOI: 10.1088/1674-1056/22/3/030312 Abstract ( 748 )   PDF (196KB) ( 509 )   Probabilistic quantum cloning (PQC) cannot copy a set of linearly dependent quantum states. In this paper, we show that if incorrect copies are allowed to be produced, linearly dependent quantum states may also be cloned by the PQC. By exploiting this kind of PQC to clone a special set of three linearly dependent quantum states, we derive the upper bound of the maximum confidence measure of a set. An explicit transformation of the maximum confidence measure is presented.
 Two-qubit and three-qubit controlled gates with cross-Kerr nonlinearity Zhao Rui-Tong (赵瑞通), Guo Qi (郭奇), Cheng Liu-Yong (程留永), Sun Li-Li (孙立莉), Wang Hong-Fu (王洪福), Zhang Shou (张寿) Chin. Phys. B, 2013, 22 (3):  030313.  DOI: 10.1088/1674-1056/22/3/030313 Abstract ( 699 )   PDF (247KB) ( 947 )   Schemes for two-qubit and three-qubit controlled gates based on cross-Kerr nonlinearity are proposed in this paper. The success probability of these gates can be increased by the quantum nondemolition detectors which are used to judge which paths the signal photons pass through. These schemes are nearly deterministic and require no ancilla photon. The advantages of these gates over the existing ones include Less resource consumption and higher success probability, which make our schemes more feasible with current technology.
 Complete Bell-state analysis for single-photon hybrid entangled state Hot! Sheng Yu-Bo (盛宇波), Zhou Lan (周澜), Cheng Wei-Wen (程维文), Gong Long-Yan (巩龙龑), Wang Lei (王磊), Zhao Sheng-Mei (赵生妹) Chin. Phys. B, 2013, 22 (3):  030314.  DOI: 10.1088/1674-1056/22/3/030314 Abstract ( 684 )   PDF (257KB) ( 728 )   We propose a scheme capable of performing complete Bell-state analysis for single-photon hybrid entangled state. Our single-photon state is encoded in both polarization and frequency degrees of freedom. The setup of scheme is composed of polarizing beam splitters, half wave plates, frequency shifters, and independent wavelength division multiplexers, which are feasible with current technology. We also show that, with this setup, we can perform the schemes for complete two-photon Bell-state analysis for polarization degree of freedom. Moreover, it can also be used to perform the teleportation scheme between different degrees of freedom. This setup may allow extensive applications in current quantum communications.
 Three-dimensional Bose–Einstein condensate vortex soliton sunder optical lattice and harmonic confinements Wang Ying (王莹), Zong Feng-De (宗丰德), Li Feng-Bo (李峰波) Chin. Phys. B, 2013, 22 (3):  030315.  DOI: 10.1088/1674-1056/22/3/030315 Abstract ( 844 )   PDF (1183KB) ( 720 )   We predict three-dimensional vortex solitons in a Bose–Einstein condensate under a complex potential which is the combination of a two-dimensional parabolic trap along the transverse radial direction and a one-dimensional optical-lattice potential along the z axis direction. The vortex solitons are built in the form of layer-chain structure made up of several fundamental vortices along the optical-lattice direction, which were not reported before in the three-dimensional Bose–Einstein condensate. By using the combination of the energy density functional method with the direct numerical simulation, we find three-dimensional vortex solitons with topological charge χ=1, χ=2, and χ=3. Moreover, the macroscopic quantum tunneling and the chirp phenomena of the vortex solitons are shown in the evolution. Thereinto, the occurrence of the macroscopic quantum tunneling provides a possibility for the realization of the quantum tunneling in experiment. Specifically, we manipulate the vortex solitons along the optical lattice direction successfully. The stability limits for dragging the vortex solitons from an initial fixed position to a prescribed location are further pursued.
 Singularities of noncompact charged objects M. Sharif, G. Abbas Chin. Phys. B, 2013, 22 (3):  030401.  DOI: 10.1088/1674-1056/22/3/030401 Abstract ( 760 )   PDF (187KB) ( 427 )   We formulate a model of noncompact spherical charged objects in the framework of noncommutative field theory. The Einstein–Maxwell field equations are solved with charged anisotropic fluid. We choose matter and charge densities as functions of two parameters instead of defining these quantities in terms of Gaussian distribution function. It is found that the corresponding densities and the Ricci scalar are singular at origin, whereas the metric is nonsingular, indicating a spacelike singularity. The numerical solution of the horizon equation implies that there are two or one or no horizon(s) depending on the mass. We also evaluate the Hawking temperature, and find that a black hole with two horizons is evaporated to an extremal black hole with one horizon.
 Thermodynamic properties of Reissner–Nordström–de Sitter quintessence black holes Wei Yi-Huan (魏益焕), Ren Jun (任军) Chin. Phys. B, 2013, 22 (3):  030402.  DOI: 10.1088/1674-1056/22/3/030402 Abstract ( 1033 )   PDF (225KB) ( 541 )   Mass, electric charge, and temperature of black holes in the Reissner–Nordström–de Sitter quintessence (RN–dSQ) spacetime are obtained. The heat capacities of the RN–dSQ black hole for certain electric charge and mass are analyzed. The electrostatic energy and the dark energy in the RN–dSQ black hole are also calculated.
 Quantum mechanical version of classical Liouville theorem Xie Chuan-Mei (谢传梅), Fan Hong-Yi (范洪义) Chin. Phys. B, 2013, 22 (3):  030501.  DOI: 10.1088/1674-1056/22/3/030501 Abstract ( 915 )   PDF (192KB) ( 1803 )   In terms of the coherent state evolution in phase space, we present a quantum mechanical version of the classical Liouville theorem. The evolution of coherent state from |z〉to |sz-rz*〉angle corresponds to the motion from a point z(q,p) to another point sz-rz* with |s|2-|r|2=1. The evolution is governed by the so-called Fresnel operator U(s,r) recently proposed in quantum optics theory, which classically corresponds to the matrix optics law and the optical Fresnel transformation and obeys the group product rules. In another word, we can recapitulate the Liouville theorem in the context of quantum mechanics by virtue of coherent state evolution in phase space, which seems to be a combination of quantum statistics and quantum optics.
 Noise destroys the coexisting of periodic orbits of a piecewise linear map Wang Can-Jun (王参军), Yang Ke-Li (杨科利), Qu Shi-Xian (屈世显) Chin. Phys. B, 2013, 22 (3):  030502.  DOI: 10.1088/1674-1056/22/3/030502 Abstract ( 828 )   PDF (406KB) ( 414 )   The effects of the Gaussian white noise and the Gaussian colored noise on the periodic orbits of period-5 (P-5) and period-6 (P-6) in their coexisting domain of a piecewise linear map are investigated numerically. The probability densities of some orbits are calculated. When the noise intensity is D=0.0001, only the orbits of P-5 exist, and the coexisting phenomenon is destroyed. On the other hand, the self-correlation time τ of the colored noise also affects the coexisting phenomenon. When τc<τ<τc', only the orbits of P-5 appear, and the stability of the orbits of P-5 is enhanced. However, when τ>τc' (τc and τc' are critical values), only the orbits of P-6 exist, and the stability of the orbits of P-6 is enhanced greatly. When τ<τc, the orbits of P-5 and P-6 coexist, but the stability of the orbits of P-5 is enhanced and that of P-6 is weakened with τ increasing.
 High frequency forcing on nonlinear systems Yao Cheng-Gui (姚成贵), He Zhi-Wei (何志威), Zhan Meng (占萌) Chin. Phys. B, 2013, 22 (3):  030503.  DOI: 10.1088/1674-1056/22/3/030503 Abstract ( 699 )   PDF (3954KB) ( 961 )   High-frequency signals are pervasive in many science and engineering fields. In this work, the effect of high-frequency driving on general nonlinear systems is investigated, and an effective equation for the slow motion is derived by extending the inertial approximation for the direct separation of fast and slow motions. Based on this theory, a high-frequency force can induce various phase transitions of a system by changing its amplitude and frequency. Numerical simulations on several nonlinear oscillator systems show very good agreements with the theoretic result. These findings may shed light on our understanding of the dynamics of nonlinear systems subject to a periodic force.
 Correlation between detrended fluctuation analysis and the Lempel–Ziv complexity in nonlinear time series analysis Tang You-Fu (唐友福), Liu Shu-Lin (刘树林), Jiang Rui-Hong (姜锐红), Liu Ying-Hui (刘颖慧) Chin. Phys. B, 2013, 22 (3):  030504.  DOI: 10.1088/1674-1056/22/3/030504 Abstract ( 727 )   PDF (724KB) ( 778 )   We focus on the study of the correlation between the detrended fluctuation analysis (DFA) and the Lempel–Ziv complexity (LZC) in nonlinear time series analysis in this paper. Typical dynamical systems including logistic map and Duffing model are investigated. Moreover, the influences of the Gaussian random noise on both DFA and LZC are analyzed. The results show a high correlation between DFA and LZC, which can quantify the non-stationarity and the nonlinearity of the time series, respectively. With the enhancement of the random component, the exponent α and the normalized complexity index C show increasing trends. In addition, C is found to be more sensitive to the fluctuation in the nonlinear time series than α. Finally, the correlation between DFA and LZC is applied to the feature extraction of vibration signals for a reciprocating compressor gas valve, and an effective fault diagnosis result is obtained.
 Hybrid three-dimensional variation and particle filtering for nonlinear systems Leng Hong-Ze (冷洪泽), Song Jun-Qiang (宋君强) Chin. Phys. B, 2013, 22 (3):  030505.  DOI: 10.1088/1674-1056/22/3/030505 Abstract ( 932 )   PDF (415KB) ( 473 )   This work addresses the problem of estimating the states of nonlinear dynamic systems with sparse observations. We present a hybrid three-dimensional variation (3DVar) and particle piltering (PF) method, which combines the advantages of 3DVar and particle-based filters. By minimizing the cost function, this approach will produce a better proposal distribution of the state. Afterwards the stochastic resampling step in standard PF can be avoided through a deterministic scheme. The simulation results show that the performance of the new method is superior to the traditional ensemble Kalman filtering (EnKF) and the standard PF, especially in highly nonlinear systems.
 Transfer function modeling and analysis of the open-loop Buck converter using the fractional calculus Wang Fa-Qiang (王发强), Ma Xi-Kui (马西奎) Chin. Phys. B, 2013, 22 (3):  030506.  DOI: 10.1088/1674-1056/22/3/030506 Abstract ( 1075 )   PDF (457KB) ( 21064 )   Based on the fact that the real inductor and the real capacitor are fractional order in nature and the fractional calculus, the transfer function modeling and analysis of the open-loop Buck converter in continuous conduction mode (CCM) operation are carried out in this paper. The fractional order small signal model and the corresponding equivalent circuit of the open-loop Buck converter in CCM operation are presented. The transfer functions from the input voltage to the output voltage, from the input voltage to the inductor current, from the duty cycle to the output voltage, from the duty cycle to the inductor current, and the output impedance of the open-loop Buck converter in CCM operation are derived, and their bode diagrams and step responses are calculated, respectively. It is found that all the derived fractional order transfer functions of the system are influenced by the fractional orders of the inductor and the capacitor. Finally, the realization of the fractional order inductor and the fractional order capacitor is designed, and the corresponding PSIM circuit simulation results of the open-loop Buck converter in CCM operation are given to confirm the correctness of the derivations and the theoretical analysis.
 Generation of transient chaos from two-dimensional systems via a switching approach Sun Chang-Chun (孙常春), Xu Qi-Cheng (徐启程), Sui Ying (隋英) Chin. Phys. B, 2013, 22 (3):  030507.  DOI: 10.1088/1674-1056/22/3/030507 Abstract ( 842 )   PDF (334KB) ( 566 )   A new approach of generating transient chaos from two-dimensional (2D) continuous autonomous systems within finite time is presented. Based on an absolute-value switching law, the phenomenon of transient chaos takes place by switching between three 2D systems. Basic dynamical behaviors of the systems are investigated. Numerical examples illustrate the validity of the results.
 Adaptive control of bifurcation and chaos in a time-delayed system Li Ning (李宁), Yuan Hui-Qun (袁惠群), Sun Hai-Yi (孙海义), Zhang Qing-Ling (张庆灵) Chin. Phys. B, 2013, 22 (3):  030508.  DOI: 10.1088/1674-1056/22/3/030508 Abstract ( 756 )   PDF (391KB) ( 758 )   In this paper, the stabilization of a continuous time-delayed system is considered. To control the bifurcation and chaos in a time-delayed system, a parameter perturbation control and a hybrid control are proposed. Then, to ensure the asymptotic stability of the system in the presence of unexpected system parameter changes, the adaptive control idea is introduced, i.e., the perturbation control parameter and the hybrid control parameter are automatically tuned according to the adaptation laws, respectively. The adaptation algorithms are constructed based on the Lyapunov–Krasovskii stability theorem. The adaptive parameter perturbation control and the adaptive hybrid control methods improve the corresponding constant control methods. They have the advantages of increased stability, adaptability to the changes of the system parameters, control cost saving, and simplicity. Numerical simulations for a well-known chaotic time-delayed system are performed to demonstrate the feasibility and superiority of the proposed control methods. Besides, comparison of the two adaptive control methods are made in an experimental study.
 Chaotic dynamic behavior analysis and control for a financial risk system Zhang Xiao-Dan (张晓丹), Liu Xiang-Dong (刘祥东), Zheng Yuan (郑媛), Liu Cheng (刘澄) Chin. Phys. B, 2013, 22 (3):  030509.  DOI: 10.1088/1674-1056/22/3/030509 Abstract ( 906 )   PDF (1490KB) ( 1930 )   According to the risk management process of financial market, a financial risk dynamic system is constructed in this paper. Through analyzing the basic dynamic property, we obtain the conditions for stability and bifurcation of the system based on Hopf bifurcation theory of nonlinear dynamical systems. In order to make the system's chaos disappear, we select the feedback gain matrix to design a class of chaotic controller. Numerical simulations are performed to reveal the change process of financial market risk. It is shown that, when the parameter of risk transmission rate changes, the system gradually comes into chaos from the asymptotically stable state through bifurcation. Besides, the controller can control chaos effectively.
 Average consensus of multi-agent systems with communication time delays and noisy links Sun Yong-Zheng (孙永征), Li Wang (李望), Ruan Jiong (阮炯) Chin. Phys. B, 2013, 22 (3):  030510.  DOI: 10.1088/1674-1056/22/3/030510 Abstract ( 722 )   PDF (787KB) ( 874 )   In this paper, we consider the average-consensus problem with communication time delays and noisy links. We analyze two different cases of coupling topologies: fixed and switching topologies. By utilizing the stability theory of the stochastic differential equations, we analytically show that the average consensus could be achieved almost surely with the perturbation of noise and the communication time delays even if the time delay is time-varying. The theoretical results show that the multi-agent systems can tolerate relatively large time delays if the noise is weak, and it can tolerate relatively strong noise if the time delays are low. The simulation results show that systems with strong noise intensities yield slow convergence.
 Novel loop-like solitons for generalized Vakhnenko equation Zhang Min (张旻), Ma Yu-Lan (马玉兰), Li Bang-Qing (李帮庆) Chin. Phys. B, 2013, 22 (3):  030511.  DOI: 10.1088/1674-1056/22/3/030511 Abstract ( 820 )   PDF (1313KB) ( 516 )   A nontraveling wave solution of a generalized Vakhnenko equation arising from the high-frequent wave motion in a relaxing medium is derived via the extended Riccati mapping method. The solution includes an arbitrary function of an independent variable. Based on the solution, two hyperbolic functions are chosen to construct new solitons. Novel single-loop-like and double-loop-like solitons are found for the equation.
 Propagation of electromagnetic soliton in anisotropic biquadratic ferromagnetic medium L. Kavitha, M. Saravanan, D. Gopi Chin. Phys. B, 2013, 22 (3):  030512.  DOI: 10.1088/1674-1056/22/3/030512 Abstract ( 851 )   PDF (1073KB) ( 939 )   The information storage technology based on anisotropic ferromagnets with sufficiently high magneto-optical effects has received much attention in recent years. The magneto-optical recording combines the merits of magnetic and optical techniques. We investigate the magneto-optical effects on a biquadratic ferromagnet and show that the dynamics of the system is governed by a perturbed nonlinear Schrödinger equation. The evolutions of amplitude and velocity of the soliton are found to be time independent, thereby admitting the lossless propagation of electromagnetic soliton in the medium, which may have potential applications in accordance with the soliton based optical communication systems. We also exploit the role of perturbation, which has a significant impact on the propagation of electromagnetic soliton.
REVIEW
 Mathematical structure of the three-dimensional (3D) Ising model Zhang Zhi-Dong (张志东) Chin. Phys. B, 2013, 22 (3):  030513.  DOI: 10.1088/1674-1056/22/3/030513 Abstract ( 1133 )   PDF (394KB) ( 3893 )   An overview of the mathematical structure of the three-dimensional (3D) Ising model is given from the viewpoints of topology, algebra, and geometry. By analyzing the relationships among transfer matrices of the 3D Ising model, Reidemeister moves in the knot theory, Yang-Baxter and tetrahedron equations, the following facts are illustrated for the 3D Ising model. 1) The complexified quaternion basis constructed for the 3D Ising model naturally represents the rotation in a (3+1)-dimensional space-time as a relativistic quantum statistical mechanics model, which is consistent with the 4-fold integrand of the partition function obtained by taking the time average. 2) A unitary transformation with a matrix that is a spin representation in 2n·l·o-space corresponds to a rotation in 2n·l·o-space, which serves to smooth all the crossings in the transfer matrices and contributes the non-trivial topological part of the partition function of the 3D Ising model. 3) A tetrahedron relation would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model, and its existence is guaranteed by the Jordan algebra and the Jordan-von Neumann-Wigner procedures. 4) The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases φx, φy, and φz. The relationship with quantum field and gauge theories and the physical significance of the weight factors are discussed in detail. The conjectured exact solution is compared with numerical results, and the singularities at/near infinite temperature are inspected. The analyticity in β=1/(kBT) of both the hard-core and the Ising models has been proved only for β>0, not for β=0. Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.
ATOMIC AND MOLECULAR PHYSICS
 New atomic data for Kr XXXV useful in fusion plasma Sunny Aggarwal, Jagjit Singh, Man Mohan Chin. Phys. B, 2013, 22 (3):  033201.  DOI: 10.1088/1674-1056/22/3/033201 Abstract ( 879 )   PDF (177KB) ( 585 )   Energy levels and emission line wavelengths of high-Z materials are useful for impurity diagnostics due to its potential application in the next generation fusion devices. For this purpose, we have calculated the fine structural energies of the 67 levels belonging to the 1s2, 1s2l, 1s3l, 1s4l, 1s5l, and 1s6l configurations of Kr XXXV using GRASP (general purpose relativistic atomic structure package) code. Additionally, we have reported the transition probabilities, oscillator strengths, line strengths, and transition wavelengths for some electric dipole (E1) transitions among these levels. We predict new energy levels and radiative rates, which have not been reported experimentally or theoretically, forming the basis for future experimental work.
 Linewidth of electromagnetically induced transparency under motional averaging in coated vapor cell Xu Zhi-Xiang (徐智翔), Qu Wei-Zhi (曲伟智), Gao Ran (高然), Hu Xin-Hua (胡新华), Xiao Yan-Hong (肖艳红) Chin. Phys. B, 2013, 22 (3):  033202.  DOI: 10.1088/1674-1056/22/3/033202 Abstract ( 824 )   PDF (3295KB) ( 1388 )   The linewidth of electromagnetically induced transparency (EIT) in a coated Rb vapor cell was studied under a magnetic field gradient. The nonlinear broadening of the EIT linewidth with the magnetic field gradient was observed. It was found that the motional averaging of the field gradient was more pronounced at higher laser intensities and larger beam sizes. In the same regime, there was a small linewidth decrease with the increasing magnetic field gradient. We have established a Monte–Carlo model, which gave results in good qualitative agreement with our experiment. Physics pictures for the above phenomena were also suggested. These results provide an understanding of the EIT linewidth behavior under the motional averaging, and should be useful for applications in quantum optics and metrology based on coated vapor cells.
 High-intensity molecular harmonic generation without ionization Wang Jun (王俊), Chen Gao (陈高), Guo Fu-Ming (郭福明), Li Su-Yu (李苏宇), Chen Ji-Gen (陈基根), Yang Yu-Jun (杨玉军) Chin. Phys. B, 2013, 22 (3):  033203.  DOI: 10.1088/1674-1056/22/3/033203 Abstract ( 665 )   PDF (612KB) ( 517 )   We theoretically investigate the high-order harmonic generation from H2+ in an infrared laser field. Our numerical simulations show that there exists a highly efficient plateau structure in the molecular harmonic spectrum. Under the action of the infrared laser pulse, the bound electronic wave packet in a potential well has enough time to tunnel through the effective potential barrier, which is formed by the molecular potential and the infrared laser field, and then recombine with the neighboring nucleus emitting a harmonic photon. During the entire dynamic process, because the wave packet is mainly located in the effective potential, the diffusion effect is of no significance, and thus a highly efficient harmonic plateau can be achieved. Specifically, the cut-off frequency of the plateau is linearly scaled with the peak amplitude of the infrared laser electric field, which may open another route to examine the internuclear distance of the molecule. Furthermore, one may detect the molecular bond lengths using the harmonic plateau.
ELECTROMAGNETISM, OPTICS, ACOUSTICS, HEAT TRANSFER, CLASSICAL MECHANICS, AND FLUID DYNAMICS
 Switching the direction of spin accumulation in the spin Hall effect of light by adjusting the optical axis of an uniaxial crystal Tang Ming (唐明), Zhou Xin-Xing (周新星), Xiao Zhi-Cheng (肖志成), Luo Hai-Lu (罗海陆), Wen Shuang-Chun (文双春) Chin. Phys. B, 2013, 22 (3):  034101.  DOI: 10.1088/1674-1056/22/3/034101 Abstract ( 816 )   PDF (744KB) ( 532 )   We theoretically and experimentally investigate a switchable spin Hall effect (SHE) of light in reflection near the Brewster angle at an air–uniaxial crystal interface. We find a large transverse spin splitting near the Brewster angle, whose sign can be altered by rotating the optical axis. As an analogy of the SHE in an electronic system, a switchable spin accumulation in the SHE of light is detected. We are able to switch the direction of the spin accumulation by adjusting the optical axis angle of the uniaxial crystal. These findings may give opportunities for photon spin manipulating and developing a new generation of nano-photonic devices.
 Manipulating surface plasmon waves by transformation optics: Design examples of beam squeezer, bend, and omnidirectional absorber Yu Zhen-Zhong (余振中), Feng Yi-Jun (冯一军), Wang Zheng-Bin (王正斌), Zhao Jun-Ming (赵俊明), Jiang Tian (姜田) Chin. Phys. B, 2013, 22 (3):  034102.  DOI: 10.1088/1674-1056/22/3/034102 Abstract ( 634 )   PDF (2974KB) ( 1132 )   We present several design examples of how to apply the transformation optics and curved space under coordinate transformation to manipulating the surface plasmon waves in a controlled manner. We demonstrate in detail the design procedure of the plasmonic wave squeezer, in-plane bend and omnidirectional absorber. We show that the approximation method of modifying only the dielectric material of a dielectric–metal surface of the plasmonic device could lead to acceptable performance, which facilitates the fabrication of the device. The functionality of the proposed plasmonic device is verified using three-dimensional full-wave electromagnetic simulations. Aiming at practical realization, we also show the design of plasmonic in-plane bend and omnidirectional absorber by an alternative transformation scheme, which results in simple device structure with a tapered isotropic dielectric cladding layer on the top of the metal surface that can be fabricated with the existing nanotechnology.
RAPID COMMUNICATION
 A single diffractive optical element implementing spectrum-splitting and beam-concentration functions simultaneously with high diffraction efficiency Ye Jia-Sheng (叶佳声), Wang Jin-Ze (王进泽), Huang Qing-Li (黄庆礼), Dong Bi-Zhen (董碧珍), Zhang Yan (张岩), Yang Guo-Zhen (杨国桢) Chin. Phys. B, 2013, 22 (3):  034201.  DOI: 10.1088/1674-1056/22/3/034201 Abstract ( 875 )   PDF (581KB) ( 519 )   In this paper, a novel method is proposed and employed to design a single diffractive optical element (DOE) for implementing spectrum-splitting and beam-concentration (SSBC) functions simultaneously. We develop an optimization algorithm, through which the SSBC DOE can be optimized within an arbitrary thickness range according to the limitations of the modern photolithography technology. Theoretical simulation results reveal that the designed SSBC DOE has a high optical focusing efficiency. It is expected that the designed SSBC DOE should have practical applications in high-efficiency solar cell systems.
ELECTROMAGNETISM, OPTICS, ACOUSTICS, HEAT TRANSFER, CLASSICAL MECHANICS, AND FLUID DYNAMICS
 New formulas for normalizing photon-added (-subtracted) two-mode squeezed thermal states Hu Li-Yun (胡利云), Fan Hong-Yi (范洪义), Zhang Zhi-Ming (张智明) Chin. Phys. B, 2013, 22 (3):  034202.  DOI: 10.1088/1674-1056/22/3/034202 Abstract ( 842 )   PDF (236KB) ( 651 )   For the first time, we derive the compact forms of normalization factors for photon-added (-subtracted) two-mode squeezed thermal states by using the P-representation and the integration within an ordered product of operators (IWOP) technique. It is found that these two factors are related to the Jacobi polynomials. In addition, some new relations are presented for the Jacobi polynomials.
 Dynamic thermal modeling and parameter identification for monolithic laser diode module Li Jin-Yi (李金义), Du Zhen-Hui (杜振辉), Ma Yi-Wen (马艺闻), Xu Ke-Xin (徐可欣) Chin. Phys. B, 2013, 22 (3):  034203.  DOI: 10.1088/1674-1056/22/3/034203 Abstract ( 815 )   PDF (311KB) ( 1296 )   We improve the thermal equivalent-circuit model of the laser diode module (LDM) to evaluate its thermal dynamic property and calculate the junction temperature of the laser diode with a high accuracy. The thermal parameters and the transient junction temperature of LDM are modeled and obtained according to the temperature of the thermistor integrated in the module. Our improved thermal model is verified indirectly by monitoring the emission wavelength of the laser diode against gas absorption lines, and several thermal parameters are obtained with the temperature uncertainty of 0.01 K in the thermal dynamic process.
 Pulse collapse and blue-shifted enhanced supercontinuumin photonic crystal fiber Liu Wen-Jun (刘文军), Pang Li-Hui (庞利辉), Lin Xiang (林翔), Gao Ren-Xi (高仁喜), Song Xiao-Wei (宋晓伟) Chin. Phys. B, 2013, 22 (3):  034204.  DOI: 10.1088/1674-1056/22/3/034204 Abstract ( 857 )   PDF (1060KB) ( 547 )   The blue-shifted supercontinuum generation in a photonic crystal fiber pumped by high peak power femtosecond pulses with wavelength located in the anomalous dispersion region is investigated experimentally and numerically. The formation of a blue-shifted enhanced supercontinuum due to the pulse collapse is demonstrated. The process of the pulse collapse is measured by using the grating-eliminated no-nonsense observation of ultrafast incident laser light e-fields technique (GRENOUILLE). Numerical simulations in spectral and temporal domains are conducted. The data from numerical simulations are in good agreement with the experimental results. Our experimental results and numerical simulations show that the pulse collapse is the determining factor in the generation of blue-shifted supercontinuum.
 Impacts of higher-order dispersions and saturable nonlinearities on modulation instability in negative-refractive metamaterials Zhong Xian-Qiong (钟先琼), Cheng Ke (程科), Xiang An-Ping (向安平) Chin. Phys. B, 2013, 22 (3):  034205.  DOI: 10.1088/1674-1056/22/3/034205 Abstract ( 772 )   PDF (971KB) ( 631 )   On the basis of the standard linear stability analysis and Drude electromagnetic model, the impacts of higher-order dispersions and three kinds of typical saturable nonlinearities on modulation instability (MI) have been analyzed and calculated for negative-refractive metamaterials (MMs). Our results show that the MI gain spectra consist of only one spectral region instead of one or two regions in ordinary materials, which may be close to or far from the zero point. Particularly, the spectrum far from the zero point has a high cut-off frequency but a narrow spectral width, which is obviously beneficial to the generation of high-repetition-rate pulse trains. Moreover, MI characteristics here will vary with the normalized angular frequency which can be modified by adjusting the structures of negative-refractive MMs, signifying the controllability of bistable solitons and MI based applications. The effects of saturable nonlinearities are similar to those in ordinary materials.
 Optical manipulation of gold nanoparticles using an optical nanofiber Li Ying (李英), Hu Yan-Jun (胡艳军) Chin. Phys. B, 2013, 22 (3):  034206.  DOI: 10.1088/1674-1056/22/3/034206 Abstract ( 760 )   PDF (996KB) ( 600 )   Gold nanoparticles are gaining increasing attention due to their biological and medical applications. In this letter, we experimentally demonstrate an optical manipulation of 250-nm-diameter gold nanoparticles along an optical nanofiber (550 nm in diameter) injected by an 808-nm laser light. The nanoparticles situated in the evanescent optical field are trapped by optical gradient force and move along the direction of light propagation due to optical scattering force. The velocities reach as high as 132 μm/s at an optical power of 80 mW.
 Infrared-induced temperature distributions of solid D2 ices Bi Peng (毕鹏), Xie Duan (谢端), Lin Wei (林伟), Wang Kai (王凯), Liu Jiang-Ping (刘江平), Tang Yong-Jian (唐永建), Yang Xiang-Dong (杨向东), Lei Hai-Le (雷海乐) Chin. Phys. B, 2013, 22 (3):  034401.  DOI: 10.1088/1674-1056/22/3/034401 Abstract ( 757 )   PDF (1232KB) ( 678 )   A specific-wavelength infrared (IR) light (λ=3140 nm) was irradiated into a solid D2 ice prepared in a cylinder target cell. The temperature in the solid D2 ice oscillated periodically with a high amplitude when irradiated by the IR light. The temperature oscillation has been well explained based on the two-dimensional heat transfer theory plus the IR-irradiation effect. The transmission optical imaging reveals that such a temperature oscillation is favorable to recrystallize the solid D2 ice from multicrystal to quasi single crystal. This suggests an efficient method to layer the solid hydrogen-isotope ice for the inertial-confinement-fusion (ICF) experiments.
 Mean and fluctuating velocity fields of a diamond turbulent jet Xu Min-Yi (徐敏义), Zhang Jian-Peng (张健鹏), Mi Jian-Chun (米建春), Nathan G. J., Kalt P. A. M. Chin. Phys. B, 2013, 22 (3):  034701.  DOI: 10.1088/1674-1056/22/3/034701 Abstract ( 734 )   PDF (1892KB) ( 559 )   The present paper reports the first investigation on a turbulent jet issuing from a diamond orifice (hereafter termed ''diamond jet'') with an aspect ratio of 1.7. Velocity measurements were conducted in the transitional region, and the exit Reynolds number of the jet was 50000. For comparison, a round jet with identical normalized boundary conditions was also measured. It is shown that the diamond jet decays and spreads faster than the round jet does over the measured flow region. The axis-switching phenomenon is observed in the diamond jet. Although both jets display primary coherent structures in the near field, these structures are found to break down more rapidly in the diamond jet, due to the higher three-dimensionality of the flow. Moreover, the streamwise components of the Reynolds normal stress and all the shear stresses reach their maxima around the location of the maximal mean shear while the maxima of the lateral components of the Reynolds normal stresses occur around the centreline of the jet.
PHYSICS OF GASES, PLASMAS, AND ELECTRIC DISCHARGES
 Effective ionization coefficients and electron drift velocities in gasmixtures of CF3I with N2 and CO2 obtained from Boltzmannequation analysis Deng Yun-Kun (邓云坤), Xiao Deng-Ming (肖登明) Chin. Phys. B, 2013, 22 (3):  035101.  DOI: 10.1088/1674-1056/22/3/035101 Abstract ( 662 )   PDF (428KB) ( 1368 )   The electron swarm parameters including the density-normalized effective ionization coefficients (α-η)/N and the electron drift velocities Ve are calculated for the gas mixture of CF3I with N2 and CO2 by solving the Boltzmann equation in the condition of steady-state Townsend (SST) experiment. The overall density-reduced electric field strength is from 100 Td to 1000 Td (1 Td=10-17 V·cm2), while the CF3I content k in the gas mixture can be varied over the range from 0% to 100%. From the variation of (α-η)/N with the CF3I mixture ratio k, the limiting field strength (E/N)lim for each CF3I concentration is derived. It is found that for the mixtures with 70% CF3I, the values of (E/N)lim are essentially the same as that for pure SF6. Additionally, the global warming potential (GWP) and the liquefaction temperature of the gas mixtures are also taken into account to evaluate the possibility of applying in the gas insulation of power equipment.
 Effect on Landau damping rates for non-Maxwellian distribution function consisting of two electron populations M. N. S. Qureshi, S. Sehar, H. A. Shah, J. B. Cao Chin. Phys. B, 2013, 22 (3):  035201.  DOI: 10.1088/1674-1056/22/3/035201 Abstract ( 866 )   PDF (352KB) ( 803 )   In many physical situations where a laser or electron beam passes through a dense plasma, hot low-density electron populations can be generated, resulting in a particle distribution function consisting of a dense cold population and a small hot population. Presence of such low-density electron distributions can alter the wave damping rate. Kinetic model is employed to study the Landau damping of Langmuir waves when a small hot electron population is present in the dense cold electron population with non-Maxwellian distribution functions. Departure of plasma from Maxwellian distributions significantly alters the damping rates as compared to the Maxwellian plasma. Strong damping is found for highly non-Maxwellian distributions as well as plasmas with higher dense and hot electron population. Existence of weak damping is also established when the distribution contains broadened flat tops at the low energies or tends to be Maxwellian. These results may be applied in both experimental and space physics regimes.
 Kadomtsev–Petviashvili equation for dust ion-acoustic solitons in pair-ion plasmas Hafeez Ur-Rehman Chin. Phys. B, 2013, 22 (3):  035202.  DOI: 10.1088/1674-1056/22/3/035202 Abstract ( 838 )   PDF (387KB) ( 1258 )   Using the reductive perturbation method, we have derived the Kadomtsev–Petviashvili (KP) equation to study the nonlinear properties of electrostatic collisionless dust ion-acoustic solitons in the pair-ion (p-i) plasmas. We have chosen the fluid model for the positive ions, the negative ions, and a fraction of static charged (both positively and negatively) dust particles. Numerical solutions of these dust ion-acoustic solitons are plotted and their characteristics are discussed. It is found that only the amplitudes of the electrostatic dust ion-acoustic solitons vary when the dust is introduced in the pair-ion plasma. It is also noticed that the amplitude and the width of these solitons both vary when the thermal energy of the positive or negative ions is varied. It is shown that potential hump structures are formed when the temperature of the negative ions is higher than that of the positive ions, and potential dip structures are observed when the temperature of the positive ions supersedes that of the negative ions. As the pair-ion plasma mimics the electron–positron plasma, thus our results might be helpful in understanding the nonlinear dust ion acoustic solitary waves in super dense astronomical bodies.
 Comparative investigation of resistance and ability to trigger high voltage discharge for single and multiple femtosecond filaments in air Guo Kai-Min (郭凯敏), Hao Zuo-Qiang (郝作强), Lin Jing-Quan (林景全), Sun Chang-Kai (孙长凯), Gao Xun (高勋), Zhao Zhen-Ming (赵振明) Chin. Phys. B, 2013, 22 (3):  035203.  DOI: 10.1088/1674-1056/22/3/035203 Abstract ( 824 )   PDF (274KB) ( 494 )   A comparative investigation of resistance and ability to trigger high voltage (HV) discharge for single filament (SF) and multiple filaments (MF) has been carried out. The experimental results show that the trend of the breakdown threshold of the SF exactly follows that of its resistance, but this is not the case for the MF. The MF's resistance is much smaller than SF's. However, the MF shows a bit higher HV breakdown threshold than the SF. The underlying physics is that the measured resistance of the MF is collectively contributed by every filament in the MF while the HV breakdown threshold is determined by only one single discharging path.
 Effect of pulsed bias on properties of ZrN/TiZrN films deposited by cathodic vacuum arc Zhang Guo-Ping (张国平), Wang Xing-Quan (王兴权), Lü Guo-Hua (吕国华), Zhou Lan (周澜), Huang Jun (黄骏), Chen Wei (陈维), Yang Si-Ze (杨思泽) Chin. Phys. B, 2013, 22 (3):  035204.  DOI: 10.1088/1674-1056/22/3/035204 Abstract ( 682 )   PDF (835KB) ( 800 )   ZrN/TiZrN multilayer are deposited by cathodic vacuum arc method with different substrate bias (from 0 to -800 V), using Ti and Zr plasma flows in residual N2 atmosphere, combined with ion bombardment of sample surfaces. The effect of pulsed bias on structure and properties of films is investigated. Microstructure of the coating is analyzed by X-ray diffraction (XRD), and scanning electron microscopy (SEM). Meanwhile, the nanohardness, Young's modulus, and scratch tests are performed. The experimental results show that the films exhibit a nanoscale multilayer structure consisting of TiZrN and ZrN phases. Solid solutions are formed for component TiZrN films. The dominant preferred orientation of TiZrN films is (111) and (220). At pulsed bias of -200 V, the nanohardness and the adhesion strength of ZrN/TiZrN multilayer reach a maximum of 38 GPa, and 78 N, respectively. The ZrN/TiZrN multilayer demonstrates an enhanced nanohardness compared with binary TiN and ZrN films deposited under equivalent conditions.
CONDENSED MATTER: STRUCTURAL, MECHANICAL, AND THERMAL PROPERTIES
 Band structure characteristics of T-square fractal phononic crystals Liu Xiao-Jian (刘晓健), Fan You-Hua (凡友华) Chin. Phys. B, 2013, 22 (3):  036101.  DOI: 10.1088/1674-1056/22/3/036101 Abstract ( 764 )   PDF (521KB) ( 700 )   The T-square fractal two-dimensional phononic crystal model is presented in this article. A comprehensive study is performed for the Bragg scattering and locally resonant fractal phononic crystal. We find that the band structures of the fractal and non-fractal phononic crystals at the same filling ratio are quite different through using the finite element method. The fractal design has important impact on the band structures of the two-dimensional phononic crystals.
 Optical study of D–D neutron irradiation-induced defects in Co- and Cu-doped ZnO wafers Wang Yun-Bo (王云波), Li Gong-Ping (李公平), Xu Nan-Nan (许楠楠), Pan Xiao-Dong (潘小东) Chin. Phys. B, 2013, 22 (3):  036102.  DOI: 10.1088/1674-1056/22/3/036102 Abstract ( 669 )   PDF (487KB) ( 617 )   Room-temperature photoluminescence and optical transmittance spectroscopy of Co-doped (1×1014,5×1016, and 1×1017 cm-2) and Cu-doped (5×1016 cm-2) ZnO wafers irradiated by D–D neutrons (fluence of 2.9×1010 cm-2) have been investigated. After irradiation, the Co or Cu metal and oxide clusters in doped ZnO wafers are dissolved, and the würtzite structure of ZnO substrate for each sample remains unchanged and keeps in high c-axis preferential orientation. The degree of irradiation-induced crystal disorder reflected from absorption band tail parameter (E0) is far greater for doped ZnO than undoped one. Under the same doping concentration, the Cu-doped ZnO wafer has much higher irradiation-induced disorder than the Co-doped one. Photoluminescence measurements indicate that the introduction rate of both zinc vacancy and zinc interstitial is much higher for the doped ZnO wafer with high doping level than the undoped one. In addition, both crystal lattice distortion and defect complexes are suggested to be formed in doped ZnO wafers. Consequently, the Co- or Cu-doped ZnO wafer (especially with high doping level) exhibits very low radiation hardness compared with the undoped one, and the Cu-doped ZnO wafer is much less radiation-hard than the Co-doped one.
 Radiation damage effects on power VDMOS devices with composite SiO2–Si3N4 films Gao Bo (高博), Liu Gang (刘刚), Wang Li-Xin (王立新), Han Zheng-Sheng (韩郑生), Song Li-Mei (宋李梅), Zhang Yan-Fei (张彦飞), Teng Rui (腾瑞), Wu Hai-Zhou (吴海舟) Chin. Phys. B, 2013, 22 (3):  036103.  DOI: 10.1088/1674-1056/22/3/036103 Abstract ( 799 )   PDF (523KB) ( 634 )   Total dose effects and single event effects on radiation-hardened power vertical double-diffusion metal oxide semiconductor (VDMOS) devices with composite SiO2-–Si3N4 film gate are investigated. The relationships among the important electrical parameters of the samples with different thickness SiO2-–Si3N4 films, such as threshold voltage, breakdown voltage, and on-state resistance in accumulated dose, are discussed. The total dose experiment results show that the breakdown voltage and the on-state resistance barely change with the accumulated dose. However, the relationships between the threshold voltages of the samples and the accumulated dose are more complex, not only positive drift, but also negative drift. At the end of the total dose experiment, we select the group of samples which have the smaller threshold voltage shift to carry out the single event effect studies. We find that the samples with appropriate thickness ratio SiO2-–Si3N4 films have a good radiation-hardening ability. This method may be useful in solving both the SEGR and the total dose problems with the composite SiO2-–Si3N4 films.
 A Dugdale–Barenblatt model for a strip with a semi-infinite crack embedded in decagonal quasicrystals Li Wu (李梧), Xie Ling-Yun (解凌云) Chin. Phys. B, 2013, 22 (3):  036201.  DOI: 10.1088/1674-1056/22/3/036201 Abstract ( 798 )   PDF (253KB) ( 715 )   The present study is to determine the solution of a strip with a semi-infinite crack embedded in decagonal quasicrystals, which transforms a physically and mathematically daunting problem. Then cohesive forces are incorporated into a plastic strip in the elastic body for nonlinear deformation. By superposing the two linear elastic fields, one is evaluated with internal loadings and the other with cohesive forces, the problem is treated in Dugdale–Barenblatt manner. A simple but yet rigorous version of the complex analysis theory is employed here, which involves conformal mapping technique. The analytical approach leads to the establishment of a few equations, which allows the exact calculation of the size of cohesive force zone and the most important physical quantity in crack theory: stress intensity factor. The analytical results of the present study may be used as the basis of fracture theory of decagonal quasicrystals.
 Site occupation of doping La3+ cations and phase transition in Na0.5Bi0.5TiO3–BaTiO3 solid solution Liu Li-Ying (刘立英), Wang Ru-Zhi (王如志), Zhu Man-Kang (朱满康), Hou Yu-Dong (侯育冬) Chin. Phys. B, 2013, 22 (3):  036401.  DOI: 10.1088/1674-1056/22/3/036401 Abstract ( 756 )   PDF (880KB) ( 669 )   Effects of La doping on the ferroelectric properties of 0.92Na0.5Bi0.5TiO3–0.08BaTiO3 (NBT–BT) solid solution have been studied both experimentally and theoretically. The experimental results show that an abnormal ferro-to-antiferroelectric phase transition is induced by La doping in NBT–BT. The first-principles calculations indicate that La3+ cations selectively substitute for the A site in NBT–BT as donors. Furthermore, the computed binding energy reveals that La cations is most likely to substitute for Ba2+ or Na+ to Bi3+ at A site as donors in NBT–BT, as supported by our Raman spectra. The ferro-to-antiferroelectric phase transition of La-doped NBT–BT is believed to originate from the lattice aberrance and redistribution of valence electrons, thus strengthening the bonding of A–O, enhancing the hybridization between the A cation d orbital and O 2p orbital, and resulting in the deflection of the polar direction of NBT–BT lattice.
 Relation between Tolman length and isothermal compressibility for simple liquids Wang Xiao-Song (王小松), Zhu Ru-Zeng (朱如曾) Chin. Phys. B, 2013, 22 (3):  036801.  DOI: 10.1088/1674-1056/22/3/036801 Abstract ( 856 )   PDF (180KB) ( 583 )   Tolman length δ0 of the liquid with a plane surface has attracted an increasing theoretical attention in recent years, while the expression of Tolman length in terms of observable quantities is still not very clear. In 2001, Bartell gave a simple expression of Tolman length δ0 in terms of isothermal compressibility. However, this expression predicts that Tolman length is always negative, which is contrary to the results of molecular dynamics simulations (MDS) for simple liquids. In this paper, this contradiction is analyzed and the reason for the discrepancy in the sign is found. In addition, we introduce a new expression of Tolman length in terms of isothermal compressibility for simple fluids not near the critical points under some week restrictions. The Tolman length of simple liquids calculated by using this formula is consistent with that obtained using MDS regarding the sign.
CONDENSED MATTER: ELECTRONIC STRUCTURE, ELECTRICAL, MAGNETIC, AND OPTICAL PROPERTIES
 Phase transition and thermodynamic properties of BiFeO3 from first-principles calculations Li Qiang (李强), Huang Duo-Hui (黄多辉), Cao Qi-Long (曹启龙), Wang Fan-Hou (王藩侯) Chin. Phys. B, 2013, 22 (3):  037101.  DOI: 10.1088/1674-1056/22/3/037101 Abstract ( 702 )   PDF (432KB) ( 1841 )   The first-principles projector-augmented wave method employing quasi-harmonic Debye model, is applied to investigate the thermodynamic properties and the phase transition between the trigonal R3c structure and the orthorhombic Pnma structure. It is found that at ambient temperature, the phase transition from the trigonal R3c phase to the orthorhombic Pnma phase is a first-order antiferromagnetic–nonmagnetic and insulator–metal transition, and occurs at 10.56 GPa, which is in good agreement with experimental data. With increasing temperature, the transition pressure decreases almost linearly. Moreover, the thermodynamic properties including Grüneisen parameter, heat capacity, entropy, and the dependences of thermal expansion coefficient on temperature and pressure are also obtained.
 Ferromagnetism in ZnO with (Mn,Li) codoping Ma Shi-Jia (马世甲), Lu Peng-Fei (芦鹏飞), Yu Zhong-Yuan (俞重远), Zhao Long (赵龙), Li Qiong-Yao (李琼瑶), Wu Cheng-Jie (武成洁), Ding Lu (丁路) Chin. Phys. B, 2013, 22 (3):  037102.  DOI: 10.1088/1674-1056/22/3/037102 Abstract ( 703 )   PDF (390KB) ( 618 )   First-principles calculations were performed to investigate the magnetic properties of Zn(Mn,Li)O based on the Perdew–Burke–Ernzerhof form of generalized gradient approximation. Antiferromagnetic (AFM) ordering is the ground state in Mn-doped ZnO system without the codopant of Li, while seven different geometrical configurations of Zn(Mn,Li)O prefer stable ferromagnetic (FM) ordering. We found that dopant Li can effectively change the magnetic coupling in the ZnMnO system. The Curie temperature (TC) of FM ordering depends on the geometric configuration, and the highest TC is about 1388 K. The FM stabilization is greatly affected by Mn–Mn distance rather than by the position of dopant Li. We propose that dopant Li mediates FM coupling through a double exchange interaction or an RKKY interaction when Li is located, respectively, near or far from Mn ions.
 Abnormal lattice contraction of plutonium hydrides studied by first-principles calculations Ao Bing-Yun (敖冰云), Shi Peng (史鹏), Guo Yong (郭咏), Gao Tao (高涛) Chin. Phys. B, 2013, 22 (3):  037103.  DOI: 10.1088/1674-1056/22/3/037103 Abstract ( 602 )   PDF (1086KB) ( 729 )   Pu can be loaded with H forming complicated continuous solid solutions and compounds, and causing remarkable electronic and structural changes. Full potential linearized augmented plane wave methods combining with Hubbard parameter U and the spin–orbit effects are employed to investigate the electronic and structural properties of stoichiometric and non-stoichiometric face-centered cubic Pu hydrides (PuHx, x=2, 2.25, 2.5, 2.75, 3). The decreasing trend with increasing x of the calculated lattice parameters is in reasonable agreement with the experimental findings. A comparative analysis of the electronic-structure results for a series of PuHx compositions reveals that the lattice contraction is resulted from the associated effects of the enhanced chemical bonding and the size effects involving the interstitial atoms. We find that the size effects are the driving force for the abnormal lattice contraction.
 Optical absorption in asymmetric double quantum wells driven by two intense terahertz fields Wu Hong-Wei (吴宏伟), Mi Xian-Wu (米贤武) Chin. Phys. B, 2013, 22 (3):  037104.  DOI: 10.1088/1674-1056/22/3/037104 Abstract ( 629 )   PDF (316KB) ( 482 )   Optical absorption is investigated for asymmetric double quantum wells driven by a resonant terahertz field and a varied terahertz field both polarized along the growth direction. Rich nonlinear dynamics of the replica peak and the Autler–Townes splitting of various dressed states are systematically studied in undoped asymmetric double quantum wells by taking account of multiple factors, such as the frequency of the varied terahertz field and the strength of the resonant terahertz field. Each electron subband splits into two dressed states when the resonant terahertz field is applied in the absence of the varied terahertz field, the optical absorption spectrum shows the first order Autler–Townes splitting of the electron subbands. When a varied terahertz field is added into the resonant system, the replica peak and the second order Autler–Townes splitting of the dressed states near the band edge respectively emerge when the varied terahertz field is non-resonant and resonant with these dressed states. When the strength of the resonant terahertz field is increased, the first order Autler–Townes double peaks and the replica peak in the optical absorption spectrum shift with the shifts of the dressed states. The presented results have potential applications in electro-optical devices.
 Analysis of resistive switching behaviors of vanadium oxide thin film Wei Xiao-Ying (韦晓莹), Hu Ming (胡明), Zhang Kai-Liang (张楷亮), Wang Fang (王芳), Zhao Jin-Shi (赵金石), Miao Yin-Ping (苗银萍) Chin. Phys. B, 2013, 22 (3):  037201.  DOI: 10.1088/1674-1056/22/3/037201 Abstract ( 832 )   PDF (689KB) ( 1421 )   We demonstrated the polarization of resistive switching for Cu/VOx/Cu memory cell. Switching behaviors of Cu/VOx/Cu cell were tested by semiconductor device analyzer (Agilent B1500A), and the relative micro-analysis of I–V characteristics of VOx/Cu was characterized by conductive atomic force microscope (CAFM). The I–V test results indicated that both forming and the reversible resistive switching between low resistance state (LRS) and high resistance state (HRS) can be observed under either positive or negative sweep. The CAFM images for LRS and HRS directly exhibited evidences of the formation and rupture of filaments based on positive or negative voltage. Cu/VOx/Cu sandwiched structure exhibits a reversible resistive switching behavior and shows potential applications in the next generation nonvolatile memory field.
 Tandem organic light-emitting diode with molybdenum tri-oxide thin film interconnector layer Lu Fei-Ping (路飞平), Wang Qian (王倩), Zhou Xiang (周翔) Chin. Phys. B, 2013, 22 (3):  037202.  DOI: 10.1088/1674-1056/22/1/037202 Abstract ( 753 )   PDF (316KB) ( 1463 )   A 10-nm thickness molybdenum tri-oxide (MoO3) thin film was used as the interconnector layer in tandem organic light-emitting devices (OLEDs). The tandem OLEDs with two identical emissive units consisting of N,N'-bis(naphthalen-1-yl)-N,N'-bis(phenyl)-benzidine (NPB)/tris(8-hydroxyquinoline) aluminum (Alq3) exhibited current efficiency–current density characteristics superior to the conventional single-unit devices. At 20 mA/cm2, the current efficiency of the tandem OLEDs using the interconnector layers of MoO3 thin film was about 4.0 cd/A, which is about twice of that of the corresponding conventional single-unit device (1.8 cd/A). The tandem OLED showed a higher power efficiency than the conventional single-unit device for luminance over 1200 cd/m2. The experimental results demonstrated that a MoO3 thin film with a proper thickness can be used as an effective interconnector layer in tandem OLEDs. Such an interconnector layer can be easily fabricated by simple thermal evaporation, greatly simplifying the device processing and fabrication processes required by previously reported interconnector layers. A possible explanation was proposed for the carrier generation of the MoO3 interconnector layer.
 Simulation of near-infrared photodiode detectors based on β-FeSi2/4H-SiC heterojunction Pu Hong-Bin (蒲红斌), He Xin (贺欣), Quan Ru-Dai (全汝岱), Cao Lin (曹琳), Chen Zhi-Ming (陈治明) Chin. Phys. B, 2013, 22 (3):  037301.  DOI: 10.1088/1674-1056/22/3/037301 Abstract ( 680 )   PDF (344KB) ( 902 )   In this paper, we propose the near-infrared p-type β-FeSi2/n-type 4H-SiC heterojunction photodetector with semiconducting silicide (β-FeSi2) as the active region for the first time. Optoelectronic characteristics of the photodetector are simulated using a commercial simulator at room temperature. The results show that the photodetector has a good rectifying character and a good response to the near-infrared light. Interface states should be minimized to obtain a lower reverse leakage current. The response spectrum of the β-FeSi2/4H-SiC detector, which consists of a p-type β-FeSi2 absorption layer with a doping concentration of 1×1015 cm-3 and a thickness of 2.5 μm, has a peak of 755 mA/W at 1.42 μm. The illumination of the SiC side obtains a higher responsivity than that of the β-FeSi2 side. The results illustrate that the β-FeSi2/4H-SiC heterojunction can be used as a near-infrared photodetector compatible with near-infrared optically-activated SiC-based power switching devices.
 Effect of microstructure on thermoelectric properties of CSD-grown Bi2Sr2Co2Oy thin films Wang Shu-Fang (王淑芳), Yan Guo-Ying (闫国英), Chen Shan-Shan (陈珊珊), Bai Zi-Long (白子龙), Wang Jiang-Long (王江龙), Yu Wei (于威), Fu Guang-Sheng (傅广生) Chin. Phys. B, 2013, 22 (3):  037302.  DOI: 10.1088/1674-1056/22/3/037302 Abstract ( 918 )   PDF (836KB) ( 695 )   Three Bi2Sr2Co2Oy thin films with different microstructures have been prepared by chemical solution deposition on LaAlO3(001) through varying the annealing temperature. With the decrease in the annealing temperature, both the size and c-axis alignment degree of grains in the film decrease as well, leading to an increase in the film resistivity. In addition, the decrease in the annealing temperature also results in a slight increase in the seebeck coefficient due to the enhanced energy filtering effect of small-grain film. The nanostructured Bi2Sr2Co2Oy film with the average grain size of about 100 nm shows a power factor comparable to that of the films with larger grains. Since the thermal conductivity of the nanostrcutured films can be depressed due to the enhanced phonon scattering by grain boundary, a higher figure of merit is expected in Bi2Sr2Co2Oy thin film with grains in nanometer size.
 Dislocation-mediated creep process in nanocrystalline Cu Mu Jun-Wei (穆君伟), Sun Shi-Cheng (孙世成), Jiang Zhong-Hao (江忠浩), Lian Jian-She (连建设), Jiang Qing (蒋青) Chin. Phys. B, 2013, 22 (3):  037303.  DOI: 10.1088/1674-1056/22/3/037303 Abstract ( 723 )   PDF (3158KB) ( 588 )   Nanocrystalline Cu with average grain sizes ranging from ～24.4 to 131.3 nm were prepared by electric brush-plating technique. Nanoindentation tests were performed within a wide strain rate range, and the creep process of nanocrystalline Cu during holding period and its relationship to dislocation and twin structures were examined. It was demonstrated that creep strain and creep strain rate are considerably significant for smaller grain size and higher loading strain rate, and are far higher than those predicted by the models of Cobble creep and grain boundary sliding. The analysis based on the calculations and experiments reveals that the significant creep deformation arises from the rapid absorption of high density dislocations stored in loading regime. Our experiments imply that stored dislocations during loading are highly unstable and dislocation activity can proceed and lead to significant post-loading plasticity.
 Wave-vector filtering effect of electric–magnetic barrier in HgTe quantum wells Zou Yong-Lian (邹永连), Song Jun-Tao (宋俊涛) Chin. Phys. B, 2013, 22 (3):  037304.  DOI: 10.1088/1674-1056/22/3/037304 Abstract ( 980 )   PDF (304KB) ( 535 )   Because of helicity of electrons in HgTe quantum wells (QWs) with inverted band structure, the electrons cannot be confined by electric barriers since electrons can tunnel the barriers perfectly without backscattering in HgTe QWs. This behavior is similar to Dirac electrons in graphene. In this paper, we propose a scheme to confine carriers in HgTe QWs using an electric–magnetic barrier. We calculate the transmission of carriers in 2-dimensional HgTe QWs and find that the wave-vector filtering effect of local magnetic fields can confine the carriers. The confining effect will have potential application in nanodevices based on HgTe QWs.
 Dynamic compensation temperature in kinetic spin-5/2 Ising model on hexagonal lattice Ümüt Temizer, Ayşegül Özkılıç Chin. Phys. B, 2013, 22 (3):  037501.  DOI: 10.1088/1674-1056/22/3/037501 Abstract ( 747 )   PDF (779KB) ( 465 )   We present a study of the dynamic behavior of a two-sublattice spin-5/2 Ising model with bilinear and crystal-field interactions in the presence of a time-dependent oscillating external magnetic field on alternate layers of a hexagonal lattice by using the Glauber-type stochastic dynamics. The lattice is formed by alternate layers of spins σ=5/2 and S=5/2. We employ the Glauber transition rates to construct the mean-field dynamic equations. First, we investigate the time variations of the average sublattice magnetizations to find the phases in the system and then the thermal behavior of the dynamic sublattice magnetizations to characterize the nature (first- or second-order) of the phase transitions and to obtain the dynamic phase transition (DPT) points. We also study the thermal behavior of the dynamic total magnetization to find the dynamic compensation temperature and to determine the type of the dynamic compensation behavior. We present the dynamic phase diagrams, including the dynamic compensation temperatures, in nine different planes. The phase diagrams contain seven different fundamental phases, thirteen different mixed phases, in which the binary and ternary combination of fundamental phases and the compensation temperature or the L-type behavior strongly depend on the interaction parameters.
 Dynamics of one-dimensional random quantum XY system with Dzyaloshinskii–Moriya interaction Li Yin-Fang (李银芳), Kong Xiang-Mu (孔祥木) Chin. Phys. B, 2013, 22 (3):  037502.  DOI: 10.1088/1674-1056/22/3/037502 Abstract ( 600 )   PDF (363KB) ( 437 )   In this paper, the effects of random variables on the dynamics of the s=1/2 XY model with the Dzyaloshinskii–Moriya interaction are studied. By means of the recurrence relation method in the high-temperature limit, we calculate the spin autocorrelation functions as well as the corresponding spectral densities for the cases that the exchange couplings between spins or external magnetic fields satisfy the double-Gaussian distribution. It is found that when the standard deviation of random exchange coupling δj (or the standard deviation of random external field δB) is small, the dynamics of the system undergoes a crossover from a collective-mode behavior to a central-peak one. However, when δJ (or δB) is large, the crossover vanishes, and the system shows a central-peak behavior or the most disordered one. We also analyze the cases in which the exchange couplings or the external fields satisfy the bimodal and the Gaussian distributions. Our results show that for all the cases considered, the dynamics of the above system is similar to that of the one-dimensional random XY model.
 Magnetic properties and magnetostriction of PrxNd1-xFe1.9 (0≤x≤1.0) alloys at low temperature Hot! Wang Yong (王永), Tang Shao-Long (唐少龙), Li Yu-Long (李玉龙), Xie Ren (谢忍), Du You-Wei (都有为) Chin. Phys. B, 2013, 22 (3):  037503.  DOI: 10.1088/1674-1056/22/3/037503 Abstract ( 785 )   PDF (385KB) ( 595 )   The crystal structure, magnetic and magnetostrictive properties of high-pressure synthesized PrxNd1-xFe1.9 (0≤x≤1.0) alloys were studied. The alloys exhibit single cubic Laves phase with MgCu2-type structure. The initial magnetization curve reveals that Pr0.2Nd0.8Fe1.9 has a minimum magnetocrystalline anisotropy at 5 K. The magnetostriction curve at 5 K shows that Pr0.2Nd0.8Fe1.9 has a very good low-field magnetostrictive property, and the magnetostriction of the PrxNd1-xFe1.9 alloy in high magnetic field is attributable mainly to Pr. The temperature dependence of the magnetostriction (λ||) at the field of 5 kOe shows that the substitution of Nd reduces the K1 remarkably, and the values of λ|| of Pr0.2Nd0.8Fe1.9 and Pr0.8Nd0.2Fe1.9 alloys are nearly five times larger than that of PrFe1.9 alloy below 50 K; the λ|| of Pr0.8Nd0.2Fe1.9 reaches up to 1082 ppm at 100 K, which makes it a potential candidate for application in this temperature range.
 Defect types and room-temperature ferromagnetism in undoped rutile TiO2 single crystals Li Dong-Xiang (李东翔), Qin Xiu-Bo (秦秀波), Zheng Li-Rong (郑黎荣), Li Yu-Xiao (李玉晓), Cao Xing-Zhong (曹兴忠), Li Zhuo-Xin (李卓昕), Yang Jing (杨静), Wang Bao-Yi (王宝义) Chin. Phys. B, 2013, 22 (3):  037504.  DOI: 10.1088/1674-1056/22/3/037504 Abstract ( 897 )   PDF (264KB) ( 947 )   Room-temperature ferromagnetism has been experimentally observed in annealed rutile TiO2 single crystals when magnetic field is applied parallel to the sample plane. By combining X-ray absorption near edge structure spectrum and positron annihilation lifetime spectroscopy, Ti3+–VO defect complexes (or clusters) have been identified in annealed crystals at high vacuum. We elucidate that the unpaired 3d electrons in Ti3+ ions provide the observed room-temperature ferromagnetism. Besides, excess oxygen ions in TiO2 lattice could induce a number of Ti vacancies which increase magnetic moments obviously.
TOPICAL REVIEW—Magnetism, magnetic materials, and interdisciplinary research
 Magnetic entropy change involving martensitic transition in NiMn-based Heusler alloys Hu Feng-Xia(胡凤霞), Shen Bao-Gen (沈保根), Sun Ji-Rong(孙继荣) Chin. Phys. B, 2013, 22 (3):  037505.  DOI: 10.1088/1674-1056/22/3/037505 Abstract ( 1203 )   PDF (1295KB) ( 1238 )   Our recent progress on magnetic entropy change (ΔS) involving martensitic transition in both conventional and metamagnetic NiMn-based Heusler alloys is reviewed. For the conventional alloys, where both martensite and austenite exhibit ferromagnetic (FM) behavior but show differentmagnetic anisotropies, a positive ΔS as large as 4.1 J·kg-1·K-1 under a field change of 0–0.9 T was first observed at martensitic transition temperature TM ~ 197 K. Through adjusting the Ni:Mn:Ga ratio to affect valence electron concentration e/a, TM was successfully tuned to room temperature, and a large negative ΔS was observed in a single crystal. The -ΔS attained 18.0 J·kg-1·K-1 under a field change of 0–5 T. We also focused on the metamagnetic alloys that show mechanisms different from the conventional ones. It was found that post-annealing in suitable conditions or introducing interstitial H atoms can shift the TM across a wide temperature range while retaining the strong metamagnetic behavior, and hence, retaining large magnetocaloric effect (MCE) and magnetoresistance (MR). The melt-spun technique can disorder atoms and make the ribbons display a B2 structure, but the metamagnetic behavior, as well as the MCE, becomes weak due to the enhanced saturated magnetization of martensites. We also studied the effect of Fe/Co co-doping in Ni45(Co1-xFex)5Mn36.6In13.4 metamagnetic alloys. Introduction of Fe atoms can assist the conversion of the Mn–Mn coupling from antiferromagnetic to ferromagnetic, thus maintaining the strong metamagnetic behavior and large MCE and MR. Furthermore, a small thermal hysteresis but significant magnetic hysteresis was observed around TM in Ni51Mn49-xInx metamagnetic systems, which must be related to different nucleation mechanisms of structural transition under different external perturbations.
 Phase transitions and magnetocaloric effects in intermetallic compounds MnFeX (X=P, As, Si, Ge) O. Tegus (特古斯), Bao Li-Hong (包黎红), Song Lin (松林) Chin. Phys. B, 2013, 22 (3):  037506.  DOI: 10.1088/1674-1056/22/3/037506 Abstract ( 1100 )   PDF (688KB) ( 1169 )   Since the discovery of giant magnetocaloric effect in MnFeP1-xAsx compounds, much valuable work has been performed to develop and improve Fe2P-type transition-metal-based magnetic refrigerants. In this article, the recent progress of our studies on fundamental aspects of theoretical considerations and experimental techniques, effects of atomic substitution on the magnetism and magnetocalorics of Fe2P-type intermetallic compounds MnFeX (X=P, As, Ge, Si) is reviewed. Substituting Si (or Ge) for As leads to an As-free new magnetic material MnFeP1-xSi(Ge)x. These new materials show large magnetocaloric effects resembling MnFe(P, As) near room temperature. Some new physical phenomena, such as huge thermal hysteresis and ‘virgin’ effect, were found in new materials. On the basis of Landau theory, a theoretical model was developed for studying the mechanism of phase transition in these materials. Our studies reveal that MnFe(P, Si) compound is a very promising material for room-temperature magnetic refrigeration and thermo-magnetic power generation.
CONDENSED MATTER: ELECTRONIC STRUCTURE, ELECTRICAL, MAGNETIC, AND OPTICAL PROPERTIES
 Carbon materials with quasi-graphene layers: the dielectric, percolation properties and the electronic transport mechanism Lu Ming-Ming (卢明明), Yuan Jie (袁杰), Wen Bo (温博), Liu Jia (刘甲), Cao Wen-Qiang (曹文强), Cao Mao-Sheng (曹茂盛) Chin. Phys. B, 2013, 22 (3):  037701.  DOI: 10.1088/1674-1056/22/3/037701 Abstract ( 810 )   PDF (1126KB) ( 931 )   We investigate the dielectric properties of multi-walled carbon nanotubes (MWCNTs) and graphite filling in SiO2 with the filling concentration of 2–20 wt.% in the frequency range of 102–107 Hz. MWCNTs and graphite have general electrical properties and percolation phenomena owing to their quasi-structure made up of graphene layers. Both permittivity ε and conductivity σ exhibit jumps around the percolation threshold. Variations of dielectric properties of the composites are in agreement with the percolation theory. All the percolation phenomena are determined by hopping and migrating electrons, which are attributed to the special electronic transport mechanism of the fillers in the composites. However, the twin-percolation phenomenon exists when the concentration of MWCNTs is between 5–10 wt.% and 15–20 wt.% in the MWCNTs/SiO2 composites, while in the graphite/SiO2 composites, there is only one percolation phenomenon in the graphite concentration of 10–15 wt.%. The unique twin-percolation phenomenon of MWCNTs/SiO2 is described and attributed to the electronic transfer mechanism, especially the network effect of MWCNTs in the composites. The formation of network plays an essential role in determining the second percolation threshold of MWCNTs/SiO2.
 Evidence of GeO volatilization and its effect on the characteristics of HfO2 grown on Ge substrate Fan Ji-Bin (樊继斌), Liu Hong-Xia (刘红侠), Fei Cheng-Xi (费晨曦), Ma Fei (马飞), Fan Xiao-Jiao (范晓娇), Hao Yue (郝跃) Chin. Phys. B, 2013, 22 (3):  037702.  DOI: 10.1088/1674-1056/22/3/037702 Abstract ( 906 )   PDF (467KB) ( 609 )   HfO2 films are deposited by atomic layer deposition (ALD) using tetrakis ethylmethylamino hafnium (TEMAH) as the hafnium precursor, while O3 or H2O is used as the oxygen precursor. After annealing at 500℃ in nitrogen, the thickness of Ge oxide's interfacial layer decreases, and the presence of GeO is observed at the H2O-based HfO2 interface due to GeO volatilization, while it is not observed for the O3-based HfO2. The difference is attributed to the residue hydroxyl groups or H2O molecules in H2O-based HfO2 hydrolyzing GeO2 and forming GeO, whereas GeO is only formed by the typical reaction mechanism between GeO2 and the Ge substrate for O3-based HfO2 after annealing. The volatilization of GeO deteriorates the characteristics of the high-κ films after annealing, which has effects on the variation of valence band offset and the C–V characteristics of HfO2/Ge after annealing. The results are confirmed by X-ray photoelectron spectroscopy (XPS) and electrical measurements.
 Spectral and laser properties of Yb and Ho co-doped (YLa)2O3 transparent ceramic Huang Dong-Dong (黄栋栋), Yang Qiu-Hong (杨秋红), Wang Yong-Gang (王永刚), Zhang Hao-Jia (张浩佳), Lu Shen-Zhou (陆神洲), Zou Yu-Wan (邹育婉), Wei Zhi-Yi (魏志义) Chin. Phys. B, 2013, 22 (3):  037801.  DOI: 10.1088/1674-1056/22/3/037801 Abstract ( 782 )   PDF (449KB) ( 643 )   Highly transparent Yb, Ho doped (YLa)2O3 ceramic was fabricated by conventional ceramic processing with nanopowders. The absorption and emission spectra of the ceramic were investigated. The energy transfer mechanism between Yb3+ and Ho3+ was also discussed. The strong emission band around 2 μm indicated that the Yb–Ho: (Y0.90La0.10)2O3 transparent ceramic is a promising gain medium for the generation of 2 μm laser emissions. The laser operation of Yb–Ho co-doped (YLa)2O3 ceramic at 2.1 μm is first reported.
 Optical properties of InGaAsBi/GaAs strained quantum wells studied by temperature-dependent photoluminescence Gu Yi (顾溢), Zhang Yong-Gang (张永刚), Song Yu-Xin (宋禹忻), Ye Hong (叶虹), Cao Yuan-Ying (曹远迎), Li Ai-Zhen (李爱珍), Wang Shu-Min (王庶民) Chin. Phys. B, 2013, 22 (3):  037802.  DOI: 10.1088/1674-1056/22/3/037802 Abstract ( 778 )   PDF (305KB) ( 1161 )   The effect of bismuth on the optical properties of InGaAsBi/GaAs quantum well structures is investigated by using the temperature-dependent photoluminescence from 12 K to 450 K. The incorporation of bismuth in the InGaAsBi quantum well is confirmed and found to result in a red shift of photoluminescence wavelength of 27.3 meV at 300 K. The photoluminescence intensity is significantly enhanced by about 50 times at 12 K with respect to that of the InGaAs quantum well due to the surfactant effect of bismuth. The temperature-dependent integrated photoluminescence intensities of the two samples reveal different behaviors related to various non-radiative recombination processes. The incorporation of bismuth also induces alloy non-uniformity in the quantum well, leading to an increased photoluminescence linewidth.
INTERDISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY
 Correlation between lithium storage and diffusion properties and electrochromic characteristics of WO3 thin films Yu Peng-Fei (于鹏飞), Cui Zhong-Hui (崔忠慧), Fan Wu-Gang (范武刚), Guo Xiang-Xin (郭向欣) Chin. Phys. B, 2013, 22 (3):  038101.  DOI: 10.1088/1674-1056/22/3/038101 Abstract ( 671 )   PDF (603KB) ( 885 )   As essential electrochromic (EC) materials are related to energy savings in fenestration technology, tungsten oxide (WO3) films have been intensively studied recently. In order to achieve better understanding of the mechanism of EC properties, and thus facilitate optimization of device performance, clarification of correlation between the cation storage and transfer properties and the coloration performance is needed. In this study, transparent polycrystalline and amorphous WO3 thin films were deposited on SnO2:F-coated glass substrates by the pulsed laser deposition technique. Investigation into optical transmittance in wavelength range of 400–800 nm measured at the current density of 130 μA/cm2 with the applied potential ranging from 3.2 to 2.2 V indicates that the polycrystalline films have a larger optical modulation of ～30% at 600 nm and a larger coloration switch time of 95 s in the whole wavelength range compared with the amorphous films (～24% and 50 s). Meanwhile, under the same conditions, the polycrystalline films show a larger lithium storage capacity corresponding to a Li/W ratio of 0.5, a smaller lithium diffusion coefficient (2×10–12 cm2·–1 for Li/W=0.24) compared with the amorphous ones with the Li/W ratio of 0.29 and the coefficient of ～2.5×10–11 cm2· s–1 as Li/W=0.24. These results demonstrate that the large optical modulation relates to the large lithium storage capacity, and the fast coloration transition associates with the fast lithium diffusion.
 Binary-component micelle and vesicle: Free energy and asymmetric distributions of amphiphiles between monolayers of vesicle Zhang Qi-Yi (张启义), Xiang Xun (向洵) Chin. Phys. B, 2013, 22 (3):  038201.  DOI: 10.1088/1674-1056/22/3/038201 Abstract ( 685 )   PDF (4215KB) ( 405 )   The real-space two-dimensional self-consistent field theory (SCFT) is employed to study the free energies of micelles and vesicles constituted by binary amphiphilic diblock copolymer AB in homopolymer A. With increasing volume fraction of copolymer AB, there are morphological transitions from the circle micelles to oblate circle-like micelles, to compound structure with inverted micelles in the inner center and micelles outer layer, and to vesicles. Special attentions are paid to the role of the copolymer AB in controlling free energies of the micelles and vesicles, by examining the effect of length ratio of A/B with the fixed whole chain length of AB copolymer, the length effect of A or B block with the corresponding fixed length of B or A block, for one component of copolymer, and the effect of different amphiphile compositions for binary-component copolymer system. The quantity η is provided to describe the asymmetric density distribution of amphiphiles between the inner and outer monolayers of vesicles, and to quantify the relative asymmetric extent of the density distribution between two species of copolymers in binary component vesicles.
 Analysis and modeling of resistive switching mechanisms oriented to resistive random-access memory Huang Da (黄达), Wu Jun-Jie (吴俊杰), Tang Yu-Hua (唐玉华) Chin. Phys. B, 2013, 22 (3):  038401.  DOI: 10.1088/1674-1056/22/3/038401 Abstract ( 838 )   PDF (273KB) ( 1205 )   With the progress of the semiconductor industry, the resistive random-access memory (RAM) has drawn increasing attention. The discovery of the memristor brings much attention to this study. Those researches focus on resistive switching characteristics of different materials and the analysis of resistive switching mechanisms. We discuss the resistive switching mechanisms of different materials in this paper and analyze the differences of those mechanisms from the view point of circuitry to establish their respective circuit models. Finally, simulations are present. We give the prospect of different materials used in resistive RAM on account of resistive switching mechanisms, which are applied to explain their resistive switchings.
 Two-dimensional threshold voltage model of nanoscale silicon-on-insulator tunneling field-effect transistor Li Yu-Chen (李妤晨), Zhang He-Ming (张鹤鸣), Zhang Yu-Ming (张玉明), Hu Hui-Yong (胡辉勇), Wang Bin (王斌), Lou Yong-Le (娄永乐), Zhou Chun-Yu (周春宇) Chin. Phys. B, 2013, 22 (3):  038501.  DOI: 10.1088/1674-1056/22/3/038501 Abstract ( 939 )   PDF (417KB) ( 1499 )   The tunneling field-effect transistor (TFET) is a potential candidate for the post-CMOS era. In this paper, a threshold voltage model is developed for this new kind of device. First, two-dimensional (2D) models are used to describe the distributions of potential and electric field in the channel and two depletion regions. Then based on the physical definition of threshold voltage for the nanoscale TFET, the threshold voltage model is developed. The accuracy of the proposed model is verified by comparing the calculated results with the 2D device simulation data. It has been demonstrated that the effects of varying the device parameters can easily be investigated using the model presented in this paper. This threshold voltage model provides a valuable reference to the TFET device design, simulation, and fabrication.
 Ballistic diffusion induced by non-Gaussian noise Qin Li (覃莉), Li Qiang (李强) Chin. Phys. B, 2013, 22 (3):  038701.  DOI: 10.1088/1674-1056/22/3/038701 Abstract ( 727 )   PDF (260KB) ( 578 )   In this letter, we have analyzed diffusive behavior of a Brownian particle subject to both internal Gaussian thermal and external non-Gaussian noise sources. We discuss two time correlation functions C(t) of the non-Gaussian stochastic process, and find that they depend on the parameter q, indicating the departure of the non-Gaussian noise from Gaussian behavior: for q≤1, C(t) is fitted very well by the first-order exponentially decaying curve and approaches zero in the long-time limit, whereas for q>1, C(t) can be approximated by a second-order exponentially decaying function and converges to a non-zero constant. Due to the properties of C(t), the particle exhibits a normal diffusion for q≤1, while for q>1 the non-Gaussian noise induces a ballistic diffusion, i.e., long-time mean square displacement of the free particle reads 〈[x(t)-〈x(t)〉]2〉∞t2.
 Feasibility of similarity coefficient map in improving morphological evaluation of T2* weighted MRI for renal cancer Wang Hao-Yu (王浩宇), Hu Jiani, Xie Yao-Qin (谢耀钦), Chen Jie (陈杰), Yu Amy, Wei Xin-Hua (魏新华), Dai Yong-Ming (戴勇鸣), Li Meng, Bao Shang-Lian (包尚联), E. M. Haacke Chin. Phys. B, 2013, 22 (3):  038702.  DOI: 10.1088/1674-1056/22/3/038702 Abstract ( 676 )   PDF (3349KB) ( 555 )   The purpose of this paper is to investigate the feasibility of similarity coefficient map (SCM) in improving morphological evaluation of T2* weighted (T2*W) magnatic resonance imaging (MRI) for renal cancer. Simulation studies and in vivo 12-echo T2*W experiments for renal cancers were performed for this purpose. The results of the first simulation study suggest that SCM can reveal small structures which are hard to be distinguished from the background tissue in T2*W images and the corresponding T2* map. The capability of improving morphological evaluation is likely due to the improvement in signal to noise ratio (SNR) and carrier to noise ratio (CNR) by SCM technique. Compared with T2*W images, SCM can improve SNR by a factor ranging from 1.87 to 2.47. Compared with T2* maps, SCM can improve SNR by a factor ranging from 3.85 to 33.31. Compared with T2*W images, SCM can improve CNR by a factor raging from 2.09 to 2.43. Compared with T2* maps SCM can improve CNR by a factor raging from 1.94 to 8.14. For a given noise level, the improvements of SNR and CNR depend mainly on the original SNRs and CNRs in T2*W images, respectively. In vivo experiments confirmed the results of the first simulation study. The results of the second simulation study suggest that more echoes are used to generate SCM, and higher SNR and CNR can be achieved in SCM. In conclusion, SCM can provide improved morphological evaluation of T2*W MR images for renal cancer by unveiling fine structures which are ambiguous or invisible in the corresponding T2*W MR images and T2* maps. What is more, in practical application, for a fixed total sampling time, one should increase the number of echoes as much as possible to achieve SCMs with better SNR and CNR.
 Finite-time consensus of heterogeneous multi-agent systems Zhu Ya-Kun (朱亚锟), Guan Xin-Ping (关新平), Luo Xiao-Yuan (罗小元) Chin. Phys. B, 2013, 22 (3):  038901.  DOI: 10.1088/1674-1056/22/3/038901 Abstract ( 638 )   PDF (1119KB) ( 1080 )   We investigate the finite-time consensus problem for the heterogeneous multi-agent systems composed of first-order and second-order agents. A novel continuous nonlinear distributed consensus protocol is constructed, and finite-time consensus criteria are obtained for the heterogeneous multi-agent systems. Compared with the existing results, the stationary and the kinetic consensuses of the heterogeneous multi-agent systems can be achieved in a finite time respectively. Moreover, the leader can be a first-order or a second-order integrator agent. Finally, some simulation examples are employed to verify the efficiency of the theoretical results.
 Measuring causality by directional symbolic mutual information approach Chen Gui (陈贵), Xie Lei (谢磊), Chu Jian (褚健) Chin. Phys. B, 2013, 22 (3):  038902.  DOI: 10.1088/1674-1056/22/3/038902 Abstract ( 961 )   PDF (253KB) ( 493 )   We propose a novel measure to assess causality through the comparison of symbolic mutual information between the future of one random quantity and the past of the other. This provides a new perspective different from the conventional conceptions. Based on this point of view, a new causality index is derived that uses the definition of directional symbolic mutual information. This measure presents properties different from the time delayed mutual information since the symbolization captures the dynamic features of the analyzed time series. In addition to characterizing the direction and the amplitude of the information flow, it can also detect coupling delays. This method has the property of robustness, conceptual simplicity, and fast computational speed.
GEOPHYSICS, ASTRONOMY, AND ASTROPHYSICS