This site uses cookies. By continuing to use this site you agree to our use of cookies. To find out more, see our Privacy and Cookies policy.

Table of contents

Volume 22

Number 6, June 2013

Previous issue Next issue

Topical reviews

Magnetism, magnetic materials and interdisciplinary research

067202

, , and

Electric-field-induced resistance switching (RS) phenomena have been studied for over 60 years in metal/dielectrics/metal structures. In these experiments a wide range of dielectrics have been studied including binary transition metal oxides, perovskite oxides, chalcogenides, carbon- and silicon-based materials, as well as organic materials. RS phenomena can be used to store information and offer an attractive performance, which encompasses fast switching speeds, high scalability, and the desirable compatibility with Si-based complementary metal—oxide—semiconductor fabrication. This is promising for nonvolatile memory technology, i.e., resistance random access memory (RRAM). However, a comprehensive understanding of the underlying mechanism is still lacking. This impedes faster product development as well as accurate assessment of the device performance potential. Generally speaking, RS occurs not in the entire dielectric but only in a small, confined region, which results from the local variation of conductivity in dielectrics. In this review, we focus on the RS in oxides with such an inhomogeneous conductivity. According to the origin of the conductivity inhomogeneity, the RS phenomena and their working mechanism are reviewed by dividing them into two aspects: interface RS, based on the change of contact resistance at metal/oxide interface due to the change of Schottky barrier and interface chemical layer, and bulk RS, realized by the formation, connection, and disconnection of conductive channels in the oxides. Finally the current challenges of RS investigation and the potential improvement of the RS performance for the nonvolatile memories are discussed.

067501

, and

Our recent research on the Mn-based antiperovskite functional materials AXMn3 (A: metal or semiconducting elements; X: C or N) is outlined. Antiperovskite carbides (e.g., AlCMn3) show large magnetocaloric effect comparable to those of typical magnetic refrigerant materials. Enhanced giant magnetoresistance up to 70% at 50 kOe (1 Oe = 79.5775 Am−1) over a wide temperature span was obtained in Ga1−xZnxCMn3 and GaCMn3−xNix. In Cu0.3Sn0.5NMn3.2, negative thermal expansion (NTE) was achieved in a wide temperature region covering room temperature (α = −6.8 ppm/K, 150 K–400 K). Neutron pair distribution function analysis suggests the Cu/Sn-Mn bond fluctuation is the driving force for the NTE in Cu1−xSnxNMn3. In CuN1−xCxMn3 and CuNMn3−yCoy, the temperature coefficient of resistivity (TCR) decreases monotonically from positive to negative as Co or C content increases. TCR is extremely low when the composition approaches the critical points. For example, TCR is ∼ 1.29 ppm/K between 240 K and 320 K in CuN0.95C0.05Mn3, which is one twentieth of that in the typical low-TCR materials (∼ 25 ppm/K). By studying the critical scaling behavior and X deficiency effect, some clues of localized-electron magnetism have been found against the background of electronic itinerant magnetism.

Topological insulator

067201

, , , and

The quantum spin Hall (QSH) effect is considered to be unstable to perturbations violating the time-reversal (TR) symmetry. We review some recent developments in the search of the QSH effect in the absence of the TR symmetry. The possibility to realize a robust QSH effect by artificial removal of the TR symmetry of the edge states is explored. As a useful tool to characterize topological phases without the TR symmetry, the spin-Chern number theory is introduced.

067301

and

Topological insulators as a new type of quantum matter materials are characterized by a full insulating gap in the bulk and gapless edge/surface states protected by the time-reversal symmetry. We propose that the interference patterns caused by the elastic scattering of defects or impurities are dominated by the surface states at the extremal points on the constant energy contour. Within such a formalism, we summarize our recent theoretical investigations on the elastic scattering of topological surface states by various imperfections, including non-magnetic impurities, magnetic impurities, step edges, and various other defects, in comparison with the recent related experiments in typical topological materials such as BiSb alloys, Bi2Te3, and Bi2Se3 crystals.

067302

, , , and

The last several years have witnessed the rapid developments in the study and understanding of topological insulators. In this review, after a brief summary of the history of topological insulators, we focus on the recent progress made in transport experiments on topological insulator films and nanowires. Some quantum phenomena, including the weak antilocalization, the Aharonov—Bohm effect, and the Shubnikov—de Haas oscillations, observed in these nanostructures are described. In addition, the electronic transport evidence of the superconducting proximity effect as well as an anomalous resistance enhancement in topological insulator/superconductor hybrid structures is included.

067304

, , , and

Providing the strong spin-orbital interaction, Bismuth is the key element in the family of three-dimensional topological insulators. At the same time, Bismuth itself also has very unusual behavior, existing from the thinnest unit to bulk crystals. Ultrathin Bi (111) bilayers have been theoretically proposed as a two-dimensional topological insulator. The related experimental realization achieved only recently, by growing Bi (111) ultrathin bilayers on topological insulator Bi2Te3 or Bi2Se3 substrates. In this review, we started from the growth mode of Bi (111) bilayers and reviewed our recent progress in the studies of the electronic structures and the one-dimensional topological edge states using scanning tunneling microscopy/spectroscopy (STM/STS), angle-resolved photoemission spectroscopy (ARPES), and first principles calculations.

067305

, , , , and

Quantum Hall effect (QHE), as a class of quantum phenomena that occur in macroscopic scale, is one of the most important topics in condensed matter physics. It has long been expected that QHE may occur without Landau levels so that neither external magnetic field nor high sample mobility is required for its study and application. Such a QHE free of Landau levels, can appear in topological insulators (TIs) with ferromagnetism as the quantized version of the anomalous Hall effect, i.e., quantum anomalous Hall (QAH) effect. Here we review our recent work on experimental realization of the QAH effect in magnetically doped TIs. With molecular beam epitaxy, we prepare thin films of Cr-doped (Bi,Sb)2Te3 TIs with well-controlled chemical potential and long-range ferromagnetic order that can survive the insulating phase. In such thin films, we eventually observed the quantization of the Hall resistance at h/e2 at zero field, accompanied by a considerable drop in the longitudinal resistance. Under a strong magnetic field, the longitudinal resistance vanishes, whereas the Hall resistance remains at the quantized value. The realization of the QAH effect provides a foundation for many other novel quantum phenomena predicted in TIs, and opens a route to practical applications of quantum Hall physics in low-power-consumption electronics.

068101

, , and

This paper presents an overview of the growth of Bi2Se3, a prototypical three-dimensional topological insulator, by molecular-beam epitaxy on various substrates. Comparison is made between the growth of Bi2Se3 (111) on van der Waals (vdW) and non-vdW substrates, with attention paid to twin suppression and strain. Growth along the [221] direction of Bi2Se3 on InP (001) and GaAs (001) substrates is also discussed.

Papers

General

060201

, and

In this paper, a procedure for constructing discrete models of the high dimensional nonlinear evolution equations is presented. In order to construct the difference model, with the aid of the potential system of the original equation and compatibility condition, the difference equations which preserve all Lie point symmetries can be obtained. As an example, invariant difference models of the (2+1)-dimensional Burgers equation are presented.

060202

, and

We consider the D-dimensional Schrödinger equation under the hyperbolic potential V0(1 − coth(αr)) + V1(1 − coth(αr))2. Using a Pekeris-type approximation, the approximate analytical solutions of the problem are obtained via the supersymmetric quantum mechanics. The behaviors of energy eigenvalues versus dimension are discussed for various quantum numbers. Useful expectation values as well as the oscillator strength are obtained.

060203

, and

In this paper we solve spin-weighted spheroidal wave equations through super-symmetric quantum mechanics with a different expression of the super-potential. We use the shape invariance property to compute the "excited" eigenvalues and eigenfunctions. The results are beneficial to researchers for understanding the properties of the spin-weighted spheroidal wave more deeply, especially its integrability.

060204

, and

This paper addresses the extended (G'/G)-expansion method and applies it to a couple of nonlinear wave equations. These equations are modified the Benjamin—Bona—Mahoney equation and the Boussinesq equation. This extended method reveals several solutions to these equations. Additionally, the singular soliton solutions are revealed, for these two equations, with the aid of the ansatz method.

060205

and

A qualitative analysis method to efficiently solve the shallow wave equations is improved, so that a more complicated nonlinear Schrödinger equation can be considered. By using the detailed study, some quite strange optical solitary waves are obtained in which the bright and dark optical solitary waves are allowed to coexist.

060206

, , and

The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method (FEM) as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline (TPS) radial basis functions.

060207

and

We derive a new method for a coupled nonlinear Schrödinger system by using the square of first-order Fourier spectral differentiation matrix D1 instead of traditional second-order Fourier spectral differentiation matrix D2 to approximate the second derivative. We prove the proposed method preserves the charge and energy conservation laws exactly. In numerical tests, we display the accuracy of numerical solution and the role of the nonlinear coupling parameter in cases of soliton collisions. Numerical experiments also exhibit the excellent performance of the method in preserving the charge and energy conservation laws. These numerical results verify that the proposed method is both a charge-preserving and an energy-preserving algorithm.

060208

, , and

A new car-following model is proposed by considering information from a number of preceding vehicles with inter-vehicle communication. A supernetwork architecture is first described, which has two layers: a traffic network and a communication network. The two networks interact with and depend on each other. The error dynamic system around the steady state of the model is theoretically analyzed and some nonjam criteria are derived. A simple control signal is added to the model to analyze the criteria of suppressing traffic jams. The corresponding numerical simulations confirm the correctness of the theoretical analysis. Compared with previous studies concerning coupled map models, the controlled model proposed in this paper is more reasonable and also more effective in the sense that it takes into account the formation of traffic congestion.

060209

and

In this paper, we analyze the generalized Camassa and Holm (CH) equation by the improved element-free Galerkin (IEFG) method. By employing the improved moving least-square (IMLS) approximation, we derive the formulas for the generalized CH equation with the IEFG method. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Because there are fewer coefficients in the IMLS approximation than in the MLS approximation, and in the IEFG method, fewer nodes are selected in the entire domain than in the conventional EFG method, the IEFG method should result in a higher computing speed. The effectiveness of the IEFG method for the generalized CH equation is investigated by numerical examples in this paper.

060210

, and

A generalized Fisher equation (GFE) relates the time derivative of the average of the intrinsic rate of growth to its variance. The exact mathematical result of the GFE has been widely used in population dynamics and genetics, where it originated. Many researchers have studied the numerical solutions of the GFE, up to now. In this paper, we introduce an element-free Galerkin (EFG) method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics. Compared with other numerical methods, the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. In comparison with the traditional method, numerical solutions show that the new method has higher accuracy and better convergence. Several numerical examples are presented to demonstrate the effectiveness of the method.

060301

, and

By virtue of the entangled state representation we concisely derive some new operator identities with regard to the two-variable Hermite polynomial (TVHP). By them and the technique of integration within an ordered product (IWOP) of operators we further derive new generating function formulas of the TVHP. They are useful in quantum optical theoretical calculations. It is seen from this work that by combining the IWOP technique and quantum mechanical representations one can derive some new integration formulas even without really performing the integration.

060302

and

By virtue of the coherent state representation of the newly introduced Fresnel operator and its group product property, we obtain new decomposition of the Fresnel operator as the product of the quadratic phase operator, the squeezing operator, and the fractional Fourier transformation operator, which in turn sheds light on the matrix optics design of ABCD-systems. The new decomposition for the two-mode Fresnel operator is also obtained by the use of entangled state representation.

060303

, , and

We apply an approximation to the centrifugal term and solve the two-body spinless-Salpeter equation (SSE) with the Yukawa potential via the supersymmetric quantum mechanics (SUSYQM) for arbitrary quantum numbers. Useful figures and tables are also included.

060304

and

In this work, we study superintegrable quantum systems in two-dimensional Euclidean space and on a complex two-sphere with second-order constants of motion. We show that these constants of motion satisfy the deformed oscillator algebra. Then, we easily calculate the energy eigenvalues in an algebraic way by solving of a system of two equations satisfied by its structure function. The results are in agreement to the ones obtained from the solution of the relevant Schrödinger equation.

060305

and

The Dirac equation is solved to obtain its approximate bound states for a spin-1/2 particle in the presence of trigonometric Pöschl—Teller (tPT) potential including a Coulomb-like tensor interaction with arbitrary spin—orbit quantum number κ using an approximation scheme to substitute the centrifugal terms κ(κ ± 1)r−2. In view of spin and pseudo-spin (p-spin) symmetries, the relativistic energy eigenvalues and the corresponding two-component wave functions of a particle moving in the field of attractive and repulsive tPT potentials are obtained using the asymptotic iteration method (AIM). We present numerical results in the absence and presence of tensor coupling A and for various values of spin and p-spin constants and quantum numbers n and κ. The non-relativistic limit is also obtained.

060307

, , and

Quantum steganography that utilizes the quantum mechanical effect to achieve the purpose of information hiding is a popular topic of quantum information. Recently, El Allati et al. proposed a new quantum steganography using the GHZ4 state. Since all of the 8 groups of unitary transformations used in the secret message encoding rule change the GHZ4 state into 6 instead of 8 different quantum states when the global phase is not considered, we point out that a 2-bit instead of a 3-bit secret message can be encoded by one group of the given unitary transformations. To encode a 3-bit secret message by performing a group of unitary transformations on the GHZ4 state, we give another 8 groups of unitary transformations that can change the GHZ4 state into 8 different quantum states. Due to the symmetry of the GHZ4 state, all the possible 16 groups of unitary transformations change the GHZ4 state into 8 different quantum states, so the improved protocol achieves a high efficiency.

060308

and

This paper points out that, due to a flaw in the sender's encoding, the receiver in Gao et al.'s controlled quantum secret direct communication (CQSDC) protocol [Chin. Phys.14 (2005), No. 5, p. 893] can reveal the whole secret message without permission from the controller. An improvement is proposed to avoid this flaw.

060309

, , , , , , , and

Besides serving as promising candidates for realizing quantum computing, superconducting quantum circuits are one of a few macroscopic physical systems in which fundamental quantum phenomena can be directly demonstrated and tested, giving rise to a vast field of intensive research work both theoretically and experimentally. In this paper we report our work on the fabrication of superconducting quantum circuits, starting from its building blocks: Al/AlOx/Al Josephson junctions. By using electron beam lithography patterning and shadow evaporation, we have fabricated aluminum Josephson junctions with a controllable critical current density (jc) and wide range of junction sizes from 0.01 μm2 up to 1 μm2. We have carried out systematical studies on the oxidation process in fabricating Al/AlOx/Al Josephson junctions suitable for superconducting flux qubits. Furthermore, we have also fabricated superconducting quantum circuits such as superconducting flux qubits and charge-flux qubits.

060310

, and

We study the entanglement of dressed atom and its spontaneous emission in a three-level Λ-type closed-loop atomic system in a multi-photon resonance condition and beyond it. It is shown that the von Neumann entropy in such a system is phase-dependent, and it can be controlled by either the intensity or relative phase of applied fields. It is demonstrated that for the special case of the Rabi frequency of applied fields, the system is disentangled. In addition, we take into account the effect of Doppler broadening on the entanglement and it is found that a suitable choice of laser propagation direction allows us to obtain the steady state degree of entanglement (DEM) even in the presence of the Doppler effect.

060501

, , and

We present a study of the physical properties of the vibrational excitation in α-helicoidal macromolecular chains, caused by the interaction with acoustical and optical phonon modes. The influence of the temperature and the basic system parameters on the vibron dressing have been analyzed by employing the simple mean-field approach based on the variational extension of the Lang—Firsov unitary transformation. The applied approach predicts a region in system parameter space where one has an abrupt transition from a partially dressed (light and mobile) to a fully dressed (immobile) vibron state. We found that the boundary of this region depends on system temperature and the type of bond among structural elements in the macromolecular chain.

060502

The dynamical properties of a tumor cell growth system described by the logistic system with coupling between non-Gaussian and Gaussian noise terms are investigated. The effects of the nonextensive index q on the stationary properties and the transient properties are discussed, respectively. The results show that the nonextensive index q can induce the tumor cell numbers to decrease greatly in the case of q > 1. Moreover, the switch from the steady stable state to the extinct state is speeded up as the increases of q, and the tumor cell numbers can be more obviously restrained for a large value of q. The numerical results are found to be in basic agreement with the theoretical predictions.

060503

, , , and

The transient properties of a three-level atomic optical bistable system in the presence of multiplicative and additive noises are investigated. The explicit expressions of the mean first-passage time (MFPT) of the transition from the high intracavity intensity state to the low one are obtained by numerical computations. The impacts of the intensities of the multiplicative noise DM and the additive noise DA, the intensity of correlation between two noises λ, and the intensity of the incident light y on the MFPT are discussed, respectively. Our results show: (i) for the case of no correlation between two noises (λ = 0.0), the increase in DM and DA can lead to an increase in the probability of the transition to the low intracavity intensity state, while the increase in y can lead to a retardation of the transition; and (ii) for the case of correlation between two noises (λ = 0.0), the increase in λ can cause an increase in the probability of the transition, and the increase in DA can cause a retardation of the transition firstly and then an increase in the probability of the transition, i.e., the noise-enhanced stability is observed for the case of correlation between two noises.

060504

, and

In this paper, we propose a new formula of the real-time minimum safety headway based on the relative velocity of consecutive trains and present a dynamic model of high-speed passenger train movements in the rail line based on the proposed formula of the minimum safety headway. Moreover, we provide the control strategies of the high-speed passenger train operations based on the proposed formula of the real-time minimum safety headway and the dynamic model of high-speed passenger train movements. The simulation results demonstrate that the proposed control strategies of the passenger train operations can greatly reduce the delay propagation in the high-speed rail line when a random delay occurs.

060505

, , , , and

The statistical distribution of natural phenomena is of great significance in studying the laws of nature. In order to study the statistical characteristics of a random pulse signal, a random process model is proposed theoretically for better studying of the random law of measured results. Moreover, a simple random pulse signal generation and testing system is designed for studying the counting distributions of three typical objects including particles suspended in the air, standard particles, and background noise. Both normal and lognormal distribution fittings are used for analyzing the experimental results and testified by chi-square distribution fit test and correlation coefficient for comparison. In addition, the statistical laws of three typical objects and the relations between them are discussed in detail. The relation is also the non-integral dimension fractal relation of statistical distributions of different random laser scattering pulse signal groups.

060506

Recent studies have shown that explosive synchronization transitions can be observed in networks of phase oscillators [Gómez-Gardeñes J, Gómez S, Arenas A and Moreno Y 2011 Phys. Rev. Lett.106 128701] and chaotic oscillators [Leyva I, Sevilla-Escoboza R, Buldú J M, Sendiña-Nadal I, Gómez-Gardeñes J, Arenas A, Moreno Y, Gómez S, Jaimes-Reátegui R and Boccaletti S 2012 Phys. Rev. Lett.108 168702]. Here, we study the effect of different chaotic dynamics on the synchronization transitions in small world networks and scale free networks. The continuous transition is discovered for Rössler systems in both of the above complex networks. However, explosive transitions take place for the coupled Lorenz systems, and the main reason is the abrupt change of dynamics before achieving complete synchronization. Our results show that the explosive synchronization transitions are accompanied by the change of system dynamics.

060507

, , , and

An evolutionary network driven by dynamics is studied and applied to the graph coloring problem. From an initial structure, both the topology and the coupling weights evolve according to the dynamics. On the other hand, the dynamics of the network are determined by the topology and the coupling weights, so an interesting structure-dynamics co-evolutionary scheme appears. By providing two evolutionary strategies, a network described by the complement of a graph will evolve into several clusters of nodes according to their dynamics. The nodes in each cluster can be assigned the same color and nodes in different clusters assigned different colors. In this way, a co-evolution phenomenon is applied to the graph coloring problem. The proposed scheme is tested on several benchmark graphs for graph coloring.

060508

, and

The aim of this paper is to study the control of hyperchaotic complex nonlinear systems with unknown parameters using passive control theory. An approach is stated to design the passive controller and estimate the unknown parameters based on the property of the passive system. The feasibility and effectiveness of the proposed approach is demonstrated through its application to the hyperchaotic complex Lü system, as an example. The estimated values of the unknown parameters are calculated. The analytical form of the complex controller is derived and used in the numerical simulation to control the hyperchaotic attractors of this example. Block diagrams of this example using Matlab/Simulink are constructed after and before the control to ensure the validity of the analytical results. Other examples of hyperchaotic complex nonlinear systems can be similarly treated.

060509

, and

By deriving the discrete-time models of a digitally controlled H-bridge inverter system modulated by bipolar sinusoidal pulse width modulation (BSPWM) and unipolar double-frequency sinusoidal pulse width modulation (UDFSPWM) respectively, the performances of the two modulation strategies are analyzed in detail. The circuit parameters, used in this paper, are fixed. When the systems, modulated by BSPWM and UDFSPWM, have the same switching frequency, the stability boundaries of the two systems are the same. However, when the equivalent switching frequencies of the two systems are the same, the BSPWM modulated system is more stable than the UDFSPWM modulated system. In addition, a convenient method of establishing the discrete-time model of piecewise smooth system is presented. Finally, the analytical results are confirmed by circuit simulations and experimental measurements.

060510

The stability of impulsive fractional-order systems is discussed. A new synchronization criterion of fractional-order chaotic systems is proposed based on the stability theory of impulsive fractional-order systems. The synchronization criterion is suitable for the case of the order 0 < q ≤ 1. It is more general than those of the known results. Simulation results are given to show the effectiveness of the proposed synchronization criterion.

060511

, , and

Based on the full velocity difference model, Jiang et al. put forward the speed gradient model through the micro-macro linkage (Jiang R, Wu Q S and Zhu Z J 2001 Chin. Sci. Bull.46 345 and Jiang R, Wu Q S and Zhu Z J 2002 Trans. Res. B36 405). In this paper, the Taylor expansion is adopted to modify the model. The backward travel problem is overcome by our model, which exists in many higher-order continuum models. The neutral stability condition of the model is obtained through the linear stability analysis. Nonlinear analysis shows clearly that the density fluctuation in traffic flow leads to a variety of density waves. Moreover, the Korteweg-de Vries—Burgers (KdV—Burgers) equation is derived to describe the traffic flow near the neutral stability line and the corresponding solution for traffic density wave is derived. The numerical simulation is carried out to investigate the local cluster effects. The results are consistent with the realistic traffic flow and also further verify the results of nonlinear analysis.

060512

, and

Totally asymmetric exclusion processes at constrained m-input n-output junction points under random update are studied by theoretical calculation and computer simulation in this paper. At the junction points, the hopping rate of particles from m-input parallel lattices to n-output parallel lattices is assumed to be equal to r/n (0 < r < 1). The mean-field approach and Monte Carlo simulations show that the phase diagram can be classified into three regions at any value of r. More interestingly, there is a threshold . In the cases of r > rc and r < rc, qualitatively different phases exist in the system. With the increase of the value of m/n, the regions of (LD, LD) and (MC, LD) or (HD, LD) decrease, and the (HD, HD) is the only phase that increases in the region (LD stands for low density, HD stands for high density, and MC for maximal current). Stationary current and density profiles are calculated, showing that they are in good agreement with Monte Carlo simulations.

060513

, and

A minimal system-plus-reservoir model yielding a nonergodic Langevin equation is proposed, which originates from the cubic-spectral density of environmental oscillators and momentum-dependent coupling. This model allows ballistic diffusion and classical localization simultaneously, in which the fluctuation-dissipation relation is still satisfied but the Khinchin theorem is broken. The asymptotical equilibrium for a nonergodic system requires the initial thermal equilibrium, however, when the system starts from nonthermal conditions, it does not approach the equilibration even though a nonlinear potential is used to bound the particle, this can be confirmed by the zeroth law of thermodynamics. In the dynamics of Brownian localization, due to the memory damping function inducing a constant term, our results show that the stationary distribution of the system depends on its initial preparation of coordinate rather than momentum. The coupled oscillator chain with a fixed end boundary acts as a heat bath, which has long been used in studies of collinear atom/solid-surface scattering and lattice vibration, we investigate this problem from the viewpoint of nonergodicity.

Atomic and molecular physics

063101

, , , , , , and

A study of the electronic and structural properties of iron phthalocyanine (FePc) molecules adsorbed on coinage metal surfaces Cu (100) and Cu (110) has been conducted by means of density functional theory calculations. The strength of the molecule-substrate interactions is interpreted in terms of the lateral adsorption geometry and the site specific electronic structure of the molecule. In the case of FePc on a (100)-oriented copper surface, the benzopyrrole leg is found to be oriented at an angle of 9° or 3° from the [01-1] substrate direction. Further, an upward bend in the molecular plane ranging from 7 to 10° is also observed; giving an almost buckled shape to the molecule. However, in the case of FePc on Cu (110), neither a bend nor a sizable rotation is observed. From the knowledge of the principle structural and electronic properties, it is concluded that FePc—Cu (100) interaction is relatively stronger than FePc—Cu (110) interaction, which is further evidenced by the charge transfer, work function changes, changes in the shape of the adsorbed molecular orbitals, and the orbital shifts. Furthermore, density of states analysis shows that the valence band level shift is surface- and site-dependent.

063102

, , , and

A full-dimensional analytical potential energy surface (APES) for the F + CH4 →HF + CH3 reaction is developed based on 7127 ab initio energy points at the unrestricted coupled-cluster with single, double, and perturbative triple excitations. The correlation-consistent polarized triple-split valence basis set is used. The APES is represented with a many-body expansion containing 239 parameters determined by the least square fitting method. The two-body terms of the APES are fitted by potential energy curves with multi-reference configuration interaction, which can describe the diatomic molecules (CH, H2, HF, and CF) accurately. It is found that the APES can reproduce the geometry and vibrational frequencies of the saddle point better than those available in the literature. The rate constants based on the present APES support the experimental results of Moore et al. [Int. J. Chem. Kin.26, 813 (1994)]. The analytical first-order derivation of energy is also provided, making the present APES convenient and efficient for investigating the title reaction with quasiclassical trajectory calculations.

063201

and

Investigations show that X-ray-boosted photoionization (XBP) has the following advantages for in-situ measurements of ultrahigh laser intensity I and field envelope F(t) (time t, pulse duration τL, carrier-envelope-phase Φ): accuracy, dynamic range, and rapidness. The calculated XBP spectra resemble inversely proportional functions of the photoelectron momentum shift. The maximum momentum p and the observable value Q (defined as a double integration of a normalized photoelectron energy spectrum, PES) linearly depend on I1/2 and τL, respectively. Φ and F(t) can be determined from the PES cut-off energy and peak positions. The measurable laser intensity can be up to and over 1018 W/cm2 by using high energy X-rays and highly charged inert gases.

063401

, and

The stereodynamic properties of the F + HO (v, j) reaction are explored by quasi-classical trajectory (QCT) calculations performed on the 1A' and 3A' potential energy surfaces (PESs). Based on the polarization-dependent differential cross sections (PDDCSs) and the angular distributions of the product angular momentum with the reactant at different values of initial v or j, the results show that the product scattering and product polarization have strong links with initial vibrational—rotational numbers of v and j. The significant manifestation of the normal DCSs is that the forward scattering gradually becomes predominant with the initial vibrational excitation increasing, and the scattering angle of the HF product taking place on the 3A' potential energy surface is found to be more sensitive to the initial value of v. The product orientation and alignment are strongly dependent on the initial rovibrational excitation effect. With enhancement in the initial rovibrational excitation effect, there is an overall decrease in the product orientation as well as in the product alignment either perpendicular to the reagent relative velocity vector k or along the direction of the y axis, for which the initial rotational excitation effect is much more noticeable than the vibrational excitation effect. Moreover, the initial rovibrational excitation effect on the product polarization is more pronounced for the 3A' potential energy surface than for the 1A' potential energy surface.

063402

, and

The stereodynamic properties of the reaction C (3P) + NO(X2Π) → CN (X2Σ+) + O (3P) in different rotational states of reactant NO are studied theoretically by using the quasiclassical trajectory method on 2A'' and 2A' potential energy surfaces (PESs) at a collision energy of 0.06 eV. The vector properties in different rotational states on the two surfaces are discussed in detail. The results indicate that the rotational excitation of NO has considerable influence on the stereodynamic property of the reaction occurring on the two surfaces. At the same time, the calculated polarization-dependent differential cross sections (PDDCSs) in different initial rotational states manifest that products are strongly polarized at three scattering angles.

Classical areas of phenomenology

064101

, and

Analytical expressions for the three components of the nonparaxial propagation of a Hermite—Laguerre—Gaussian (HLG) beam in uniaxial crystal orthogonal to the optical axis are derived. The intensity distribution of an HLG beam and its three components propagating in a uniaxial crystal orthogonal to the optical axis are demonstrated by numerical examples. Although the y and z components of an HLG beam in the incident plane are both equal to zero, they emerge upon propagation inside the uniaxial crystal. Moreover, the beam profile of the x component is relatively stable and the beam profiles of the y and z components have the same evolution law. If the ratio of the extraordinary refractive index to the ordinary refractive index is larger than unity, the beam profile of the HLG beam is elongated in the x direction and generally rotates clockwise. Otherwise, the beam profile of the HLG beam is elongated in the y direction and generally rotates anticlockwise. This research is beneficial to the optical trapping and nonlinear optics involved in the rotation of a beam profile.

064201

The scattering of electromagnetic wave by an array of parallel metallic single-walled carbon nanotubes is investigated based on the boundary-value method. Electronic excitations over each nanotube surface are modeled as an infinitesimally thin cylindrical layer of the free-electron gas. The scattering cross section of both transverse magnetic (TM) and transverse electric (TE) uniform plane waves by the system at normal incidences is obtained.

064202

, , , , and

The theoretical and experimental results of tightly focused radially polarized vortex beams are demonstrated. An auto-focus technology is introduced into the measurement system in order to enhance the measurement precision, and the radially polarized vortex beams are generated by a liquid-crystal polarization converter and a vortex phase plate. The focused fields of radially polarized vortex beams with different topological charges at numerical apertures (NAs) of 0.65 and 0.85 are measured respectively, and the results indicate that the total intensity distribution at focus is dependent not only on the NA of the focusing objective lens and polarization pattern of the beam but also on the topological charge l of the beam. Some unique focusing properties of radially polarized vortex beams with fractional topological charges are presented based on numerical calculations. The experimental verification paves the way for some practical applications of radially polarized vortex beams, such as in optical trapping, near-field microscopy, and material processing.

064203

, and

Multispectral time delay and integration charge coupled device (TDICCD) image compression requires a low-complexity encoder because it is usually completed on board where the energy and memory are limited. The Consultative Committee for Space Data Systems (CCSDS) has proposed an image data compression (CCSDS-IDC) algorithm which is so far most widely implemented in hardware. However, it cannot reduce spectral redundancy in multispectral images. In this paper, we propose a low-complexity improved CCSDS-IDC (ICCSDS-IDC)-based distributed source coding (DSC) scheme for multispectral TDICCD image consisting of a few bands. Our scheme is based on an ICCSDS-IDC approach that uses a bit plane extractor to parse the differences in the original image and its wavelet transformed coefficient. The output of bit plane extractor will be encoded by a first order entropy coder. Low-density parity-check-based Slepian—Wolf (SW) coder is adopted to implement the DSC strategy. Experimental results on space multispectral TDICCD images show that the proposed scheme significantly outperforms the CCSDS-IDC-based coder in each band.

064204

, , , and

In this paper we study the atom dissipation effect in a laser cavity. The cavity field mode is described by the Fox—Li quasimode due to the leakiness of the cavity. Our results show that the atom decay rate versus the decay rate of the quasimode is a Lorentz type. Effects of the atom—cavity detuning as well as cavity size are also discussed.

064205

, , , , and

An effective scheme within two displaced bosonic operators with equal positive and negative displacements is extended to study qubit—oscillator systems analytically in a unified way. Many previous analytical treatments, such as generalized rotating-wave approximation (GRWA) [Phys. Rev. Lett.99, 173601 (2007)] and an expansion in the qubit tunneling matrix element in the deep strong coupling regime [Phys. Rev. Lett.105, 263603 (2010)] can be recovered straightforwardly within the present scheme. Moreover, further improving GRWA and the extension to the finite-bias case are implemented easily. The algebraic formulae for the eigensolutions are then derived explicitly and uniquely, which work well in a wide range of the coupling strengths, detunings, and static bias including the recent experimentally accessible parameters. The dynamics of the qubit for an oscillator in the ground state is also studied. At the experimentally accessible coupling regime, GRWA can always work well. When the coupling is enhanced to the intermediate regime, only the improving GRWA can give the correct description, while the result of GRWA shows strong deviations. The previous Van Vleck perturbation theory is not valid to describe the dynamics in the present-day experimentally accessible regime, except for the strongly biased cases.

064206

, and

We theoretically investigate the stationary entanglement of a optomechanical system with an additional Kerr medium in the cavity. There are two kinds of interactions in the system, photon—mirror interaction and photon—photon interaction. The optomechanical entanglement created by the former interaction can be effectively controlled by the latter one. We find that the optomechanical entanglement is suppressed by Kerr interaction due to photon blockage. We also find that the Kerr interaction can create the stationary entanglement and induce the resonance of entanglement in the small detuning regime. These results show that the Kerr interaction is an effective control for the optomechanical system.

064207

, , , and

A curviform surface breaks the symmetrical shape of silicon quantum dots on which some bonds can produce localized electronic states in the bandgap. The calculation results show that the bonding energy and electronic states of silicon quantum dots are different on various curved surfaces, for example, a Si—O—Si bridge bond on curved surface provides localized levels in bandgap and its bonding energy is shallower than that on the facet. The red-shifting of the photoluminescence spectrum on smaller silicon quantum dots can be explained by the curved surface effect. Experiments demonstrate that silicon quantum dots are activated for emission due to the localized levels provided by the curved surface effect.

064208

, , , and

In this paper, the theoretical rate equation model of an in-band pumped gain-switched thulium-doped fiber (TDF) laser is investigated. The analytical formulations of pump energy threshold, peak power extraction efficiency, and pulse extraction efficiency are derived through analyzing the interaction process between the pump pulse and the laser pulse. They are useful for understanding, designing, and optimizing the in-band pumped TDF lasers in a 1.9 μm–2.1 μm wavelength region. The experiment with an all-fiber gain-switched TDF laser pumped by a 1.558-μm pulse amplifier is conducted, and our experimental results show good agreement with theoretical analysis.

064209

, , , , and

The output performance of a 980-nm broad-area vertical-cavity surface-emitting laser (VCSEL) is improved by optimizing the p-electrode diameter in this study. Based on a three-dimensional finite-element method, the current density distribution within the active region of the VCSEL is optimized through the appropriate adjustment of the p-electrode diameter, and uniform current-density distribution is achieved. Then, the effects of this optimization are studied experimentally. The LIV characteristics under different temperatures of the VCSELs with different p-electrode diameters are investigated, and better temperature stability is demonstrated in the VCSEL with an optimized p-electrode diameter. The far-field measurements show that with an injected current of 2 A, the far-field divergence angle of the VCSEL with an optimized p-electrode diameter is 9°, which is much lower than the far-field angle of the VCSEL without this optimization. Also the VCSEL with an optimized p-electrode diameter shows a better near-field distribution.

064210

, , , and

Looking for new light sources, especially short wavelength laser light sources has attracted widespread attention. This paper analytically describes the radiation of a crystalline undulator field by the sine-squared potential. In the classical mechanics and the dipole approximation, the motion equation of a particle is reduced to a generalized pendulum equation with a damping term and a forcing term. The bifurcation behavior of periodic orbits is analyzed by using the Melnikov method and the numerical method, and the stability of the system is discussed. The results show that, in principle, the stability of the system relates to its parameters, and only by adjusting these parameters appropriately can the occurrence of bifurcation be avoided or suppressed.

064211

, , , , and

Broadband normal dispersion pumping supercontinuum (SC) generation in silica photonic crystal fiber (PCF) is investigated in this paper. A 1064-nm picosecond fiber laser is used to pump silica PCF for the SC generation. The length of PCF is optimized for the most efficient stimulated Raman scattering process in the picosecond pump pulse region. The first stimulated Raman Stokes peak is located in the anomalous dispersion regime of the PCF and near the zero dispersion wavelength; thus the SC generation process can benefit from both a normal dispersion pumping scheme and an anomalous dispersion pumping scheme. The 51.7-W SC spanning from about 700 nm to beyond 1700 nm is generated with an all-fiber configuration, and the pump-to-SC conversion efficiency is up to 90%. In order to avoid the output fiber end face damage and increase the stability of the system, an improved output solution for the high power SC is proposed in our experiment. This high-efficiency near-infrared SC source is very suitable for applications in which average output power and spectral power density are firstly desirable.

064212

, , , , , , , , , et al

We develop a picosecond widely tunable laser in a deep-ultraviolet region from 175 nm to 210 nm, generated by two stages of frequency doubling of a 80-MHz mode-locked picosecond Ti:sapphire laser. A β-BaB2O4 walk-off compensation configuration and a KBe2BO3F2 prism-coupled device are adopted for the generation of second harmonic and fourth harmonics, respectively. The highest power is 3.72 mW at 193 nm, and the fluctuation at 2.85 mW in 130 min is less than ±2%.

064213

, , , , , and

In this paper, spectral hole depth dependence on temperature below 10 K in Tm3+:YAG crystal is investigated in detail. A novel model is proposed to analyze the temperature dependence on the spectral hole. By using the proposed model, we theoretically deduce the temperature dependence of spectral hole depth. The results are compared with experimental results and they are in good agreement. According to the theoretic results, the optimum temperature in experiment can be found.

064214

, , and

We report on the generation of self-oscillations from a continuously pumped singly resonant frequency doubler based on a periodically poled potassium titanyl phosphate crystal (PPKTP). The sustained square-wave and staircase curve of self-oscillations are obtained when the incident pump powers are below and above the threshold of subharmonic-pumped parametric oscillation (SPO), respectively. The self-oscillations can be explained by the competition between the phase shifts induced by cascading nonlinearity and thermal effect, and the influence of fundamental nonlinear phase shift by the generation of SPO. The simulation results are in good agreement with the experiment data.

064215

and

We theoretically investigate high-order harmonic generation by employing strong-field approximation (SFA) and present a new approach to the extension of the high-order harmonic cutoff frequency via an exploration of the dependence of high-order harmonic generation on the waveform of laser fields. The dependence is investigated via detailed analysis of the classical trajectories of the ionized electron moving in the continuum in the velocity—position plane. The classical trajectory consists of three sections (Acceleration Away, Deceleration Away, and Acceleration Back), and their relationship with the electron recollision energy is investigated. The analysis of classical trajectories indicates that, besides the final (Acceleration Back) section, the electron recollision energy also relies on the previous two sections. We simultaneously optimize the waveform in all three sections to increase the electron recollision energy, and an extension of the cutoff frequency up to Ip + 20.26Up is presented with a theoretically synthesized waveform of the laser field.

064216

, , , and

A novel tunable comb filter composed of a single-mode/multimode/polarization-maintaining-fiber-based Sagnac fiber loop is proposed and experimentally demonstrated. The filter tunability is achieved by rotating the polarization controller. The spectral shift is dependent on rotation direction and the position of the polarization controller. In addition, the adjustable range achieved by rotating the half-wave-plate polarization controller is twice higher than that of the quarter-wave-plate one.

064217

, , and

The deposition of a Cu seed layer film is investigated by supercritical fluid deposition (SCFD) using H2 as a reducing agent for Bis(2,2,6,6-tetramethyl-3,5- heptanedionato) copper in supercritical CO2 (scCO2). The effects of deposition temperature, precursor, and H2 concentration are investigated to optimize Cu deposition. Continuous metallic Cu films are deposited on Ru substrates at 190 °C when a 0.002 mol/L Cu precursor is introduced with 0.75 mol/L H2. A Cu precursor concentration higher than 0.002 mol/L is found to have negative effects on the surface qualities of Cu films. For a H2 concentration above 0.56 mol/L, the root-mean-square (RMS) roughness of a Cu film decreases as the H2 concentration increases. Finally, a 20-nm thick Cu film with a smooth surface, which is required as a seed layer in advanced interconnects, is successfully deposited at a high H2 concentration (0.75 mol/L).

064401

, , , and

Mesoporous silica SBA-15 consists of uniform hexagonal, unconnected cylindrical channels with diameters that can be tuned within a range of 1.5 nm–30 nm, and is thought to have a special thermal conductivity. The theoretical investigation of the shell thermal conductivity of the mesoporous silica is performed in the relaxation time approximation in this paper and an available one-dimensional heat transfer model is used to predict the effective thermal conductivity (ETC) of the mesoporous silica. The experimental result of the ETC is also presented for comparison. The shell thermal conductivity of the mesoporous silica decreases with mesochannel radius increasing or wall thickness decreasing, but does not strictly decrease with porosity increasing. The thermal radiation possibly plays a primary role in heat transfer at the large porosity scale. The predicted ETC of SBA-15 with only conduction considered is quite low at the large porosity, even lower than the thermal conductivity of the silica aerogels. To realize it, doping carbon or other matters which can strongly absorb infrared light into SBA-15 is a possible way.

064701

, , and

The characteristics of a nitrogen arc using a graphite cathode and a melting anode in a pilot-scale plasma furnace are investigated. The voltage is examined as a function of current and apparent plasma length. The voltage increases non-linearly with the increase of apparent plasma length, with the current fixed. The experimental data so obtained are compared with the predictions of the Bowman model for the electric arc, and with numerical simulations as well. The level of agreement between the experimental data at the melting anode and the numerical predictions confirms the suitability of the proposed the Bowman model. These characteristics are relevant to the engineering design and evaluation of a DC plasma furnace and reactor for the treatment of hazardous fly ash waste.

Fluids, plasmas and electric discharges

065201

, , , , , and

An experiment with thin titanium foils irradiated by two pulses delayed in time is conducted on the Shenguang-II laser facility. A prepulse induces an underdense plasma, 2-ns later a main pulse (λL = 0.35 μm, EL ≈ 120 J, τL ≈ 100 ps) is injected into the underdense plasma and produces strong line emission from the titanium K shell (i.e., Heα at 4.7 keV). Data show that the intensity of 4.7-keV X-ray emission with the prepulse is approximately twice more than without the prepulse, and can be used as a backlighting source satisfying the diagnostic requirements for dense plasma probing. High-quality plasma images are obtained with the backlighting 4.7-keV X-rays in a Rayleigh—Taylor hydrodynamic instability experiment.

Condensed matter: structure, thermal and mechanical properties

066101

, , , , , , and

We report on the Au-assisted vapour—liquid—solid (VLS) growth of GaAs/InxGa1−xAs/GaAs (0.2 ≤ x ≤ 1) axial double-heterostructure nanowires on GaAs (111) B substrates via the metal-organic chemical vapor deposition (MOCVD) technique. The influence of the indium (In) content in an Au particle on the morphology of nanowires is investigated systematically. A short period of pre-introduced In precursor before the growth of InxGa1−xAs segment, coupled with a group III precursor interruption, is conducive to obtaining symmetrical heterointerfaces as well as the desired In/Ga ratio in the InxGa1−xAs section. The nanowire morphology, such as kinking and tapering, are thought to be related to the In composition in the catalyst alloy as well as the VLS growth mechanism.

066201

and

The structural and elastic properties of multiferroic Ca3Mn2O7 with ferroelectric orthorhombic (O-phase) and para-electric tetragonal structures (T-phase) have been studied by first-principles calculations within the generalized gradient approximation (GGA) and the GGA plus Hubbard U approaches (GGA + U). The calculated theoretical structures are in good agreement with the experimental values. The T-phase is found to be antiferromagnetic (AFM) and the AFM O-phase is more stable than the T-phase, which also agree with the experiments. On these bases, the single-crystal elastic constants (Cijs) and elastic properties of polycrystalline aggregates are investigated for the two phases. Our elasticity calculations indicate Ca3Mn2O7 is mechanically stable against volume expansions. The AFM O-phase is found to be a ductile material, while the AFM T-phase shows brittle nature and tends to be elastically isotropic. We also investigate the influence of strong correlation effects on the elastic properties, qualitatively consistent results are obtained in a reasonable range of values of U. Finally, the ionicity is discussed by Bader analysis. Our work provides useful guidance for the experimental elasticity measurements of Ca3Mn2O7, and makes the strain energy calculation in multiferroic Ca3Mn2O7 thin films possible.

066202

, , , , , , and

A new structural phase of MgV2O6 was obtained by a high-pressure, high-temperature (HPHT) synthesis method. The new phase was investigated by the Rietveld analysis of X-ray powder diffraction data, showing space group Pbcn (No. 60) symmetry and a = 13.6113(6) Å (1 Å = 0.1 nm), b = 5.5809(1) Å, c = 4.8566(3) Å, V = 368.93(2) Å3 (Z = 4). High pressure behavior was studied by Raman spectroscopy at room temperature. Under 22.5 GPa, there was no sign of a structural phase transition in the spectra, demonstrating stability of the HPHT phase up to the highest pressure.

066301

, and

This paper studies the force network properties of marginally and deeply jammed packings of frictionless soft particles from the perspective of complex network theory. We generate zero-temperature granular packings at different pressures by minimizing the inter-particle potential energy. The force networks are constructed as nodes representing particles and links representing normal forces between the particles. Deeply jammed solids show remarkably different behavior from marginally jammed solids in their degree distribution, strength distribution, degree correlation, and clustering coefficient. Bimodal and multi-modal distributions emerge when the system enters the deep jamming region. The results also show that small and large particles can show different correlation behavior in this simple system.

066801

and

An Sr/Si(100)-c(2×4) surface is investigated by high-resolution scanning tunneling microscopy (STM) and scanning tunneling spectroscopy (STS). The semiconductor property of this surface is confirmed by STS. The STM images of this surface shows that it is bias-voltage dependent and an atomic resolution image can be obtained at an empty state under a bias voltage of 1.5 V. Furthermore, one-dimensional (1D) diffusion of vacancies can be found in the room-temperature STM images. Sr vacancies diffuse along the valley channels, which are constructed by silicon dimers in the surface. Weak interaction between Sr and silicon dimers, low metal coverage, surface vacancy, and energy of thermal fluctuation at room temperature all contribute to this 1D diffusion.

066802

, , and

In the present paper we give a detailed report on the results of our first-principles investigations of Ar adsorptions at the four high symmetry sites on M (111) (M = Pd, Pt, Cu, and Rh) surfaces. Our studies indicate that the most stable adsorption sites of Ar on Pd (111) and Pt (111) surfaces are found to be the fcc-hollow sites. However, for Ar adsorptions on Cu (111) and Rh (111) surfaces, the most favorable site is the on-top site. The density of states (DOS) is analyzed for Ar adsorption on M (111) surfaces, and it is concluded that the adsorption behavior is dominated by the interaction between 3s, 3p orbits of Ar atoms and the d orbit of the base metal atoms.

066803

, , , , and

Thin films of tungsten (W)-doped thermochromic vanadium dioxide (VO2) were deposited onto soda-lime glass and fused silica by radio frequency magnetron sputtering. The doped VO2 films were characterized by X-ray diffraction, optical transmittance measurement, and near field optical microscopy with Raman spectroscopy. X-ray diffraction patterns show that the (011) peak of W-doped thermochromic VO2 film shifts to a lower diffraction angle with the increase of W concentration. The optical measurements indicated that the transmittance change (ΔT) at wavelength of 2500 nm drops from 65% (ΔT at 35 °C and 80 °C for undoped VO2 film) to 38% (ΔT at 30 °C and 42 °C for the doped VO2 film). At the same time, phase transition temperature drops from 65 °C to room temperature or lower with the increase of W concentration. Near field optical microscopy image shows that the surface of W-doped VO2 film is smooth. Raman results show that the main Raman modes of W-doped VO2 are centered at 614 cm−1, the same as that of undoped VO2, suggesting no Raman mode changes for lightly W-doped VO2 at room temperature, due to no phase transition appearing under this condition.

066804

, , and

Thin films of ternary compounds CuxInyN and CuxTiyN were grown by magnetron sputtering to improve the thermal stability of Cu3N, a material that decomposes below 300 °C, and thus promises many interesting applications in direct-writing. The effect of In or Ti incorporation in altering the structure and physical properties of copper nitride was evaluated by characterizing the film structure, surface morphology, and temperature dependence of electrical resistivity. More Ti than In can be accommodated by copper nitride without completely deteriorating the Cu3N lattice. A small amount of In or Ti can improve the crystallinity, and consequently the surface morphology. While the decomposition temperature is rarely influenced by In, the Ti-doped sample, Cu59.31Ti2.64N38.05, shows an X-ray diffraction pattern dominated by characteristic Cu3N peaks, even after annealing at 500 °C. Both In and Ti reduce the bandgap of the original Cu3N phase, resulting in a smaller electrical resistivity at room temperature. The samples with more Ti content manifest metal-semiconductor transition when cooled from room temperature down to 50 K. These results can be useful in improving the applicability of copper—nitride-based thin films.

Condensed matter: electronic structure, electrical, magnetic and optical properties

067101

, , , , and

The adsorption of hydrogen molecules on titanium-decorated (Ti-decorated) single-layer and bilayer graphenes is studied using density functional theory (DFT) with the relativistic effect. Both the local density approximation (LDA) and the generalized gradient approximation (GGA) are used for obtaining the region of the adsorption energy of H2 molecules on Ti-decorated graphene. We find that a graphene layer with titanium (Ti) atoms adsorbed on both sides can store hydrogen up to 9.51 wt% with average adsorption energy in a range from −0.170 eV to −0.518 eV. Based on the adsorption energy criterion, we find that chemisorption is predominant for H2 molecules when the concentration of H2 molecules absorbed is low while physisorption is predominant when the concentration is high. The computation results for the bilayer graphene decorated with Ti atoms show that the lower carbon layer makes no contribution to hydrogen adsorption.

067102

The electronic structures, optical dielectric functions, elastic properties, and lattice dynamics of Ba2ZnWO6 have been investigated by using the generalized gradient approximation. The density of states and distributions of charge density show that O and Ba tend toward ionic bond, but O and W or Zn display the covalent bond character. The calculated energy band structure shows that Ba2ZnWO6 is a wide indirect band gap semiconductor. The static value 2.28 of the refractive index is attained. The analysis of the elastic properties of Ba2ZnWO6 indicates a rather weak elastic anisotropy. The phonon dispersion is calculated, suggesting no structural instability, which is agreement with the recent low temperature neutron diffraction experiments. The mensurability Cv (phonon heat capacity) as the function of the temperature is also calculated to judge our results for future experiment.

067103

, , and

The inelastic scattering of oppositely charge polarons in polymer heterojunctions is believed to be of fundamental importance for the light-emitting and transport properties of conjugated polymers. Based on the tight-binding SSH model, and by using a nonadiabatic molecular dynamic method, we investigate the effects of interface hopping on inelastic scattering of oppositely charged polarons in a polymer heterojunction. It is found that the scattering processes of the charge and lattice defect depend sensitively on the hopping integrals at the polymer/polymer interface when the interface potential barrier and applied electric field strength are constant. In particular, at an intermediate electric field, when the interface hopping integral of the polymer/polymer heterojunction material is increased beyond a critical value, two polarons can combine to become a lattice deformation in one of the two polymer chains, with the electron and the hole bound together, i.e., a self-trapped polaron—exciton. The yield of excitons then increases to a peak value. These results show that interface hopping is of fundamental importance and facilitates the formation of polaron—excitons.

067104

, , , , , , , , , et al

Using measured capacitance—voltage curves and current—voltage characteristics for the AlGaN/AlN/GaN heterostructure field-effect transistors with different gate lengths and drain-to-source distances, the influence of drain bias on the electron mobility is investigated. It is found that below the knee voltage the longitudinal optical (LO) phonon scattering and interface roughness scattering are dominant for the sample with a large ratio of gate length to drain-to-source distance (here 4/5), and the polarization Coulomb field scattering is dominant for the sample with a small ratio (here 1/5). However, the above polarization Coulomb field scattering is weakened in the sample with a small drain-to-source distance (here 20 μm) compared with the one with a large distance (here 100 μm). This is due to the induced strain in the AlGaN layer caused by the drain bias.

067105

, , , , , and

An iterative method is used to find the values of the Hamiltonian parameters for Yb3+ in a given low-symmetry crystalline site. Samples of Yb3+:RETaO4 (RE = Gd, Y, and Sc) were prepared and their structures were determined. Based on the obtained structural data, their orbital-spin parameters and crystal field parameters were fitted by the superposition model (SM). Using the crystal field parameters obtained by the SM fitting as the initial parameters, the Hamiltonian parameters were fitted iteratively. The calculated and experimental energy levels for Yb3+:RETaO4 are consistent, and the maximal mean-root-square deviation is only 2.84 cm−1, indicating that the method is effective to determine the Hamiltonian parameters of Yb3+ in low-symmetry crystalline sites.

067203

, , , , and

We simulate the current—voltage (IV) characteristics of AlGaN/AlN/GaN heterostructure field-effect transistors (HFETs) with different gate lengths using the quasi-two-dimensional (quasi-2D) model. The calculation results obtained using the modified mobility model are found to accord well with the experimental data. By analyzing the variation of the electron mobility for the two-dimensional electron gas (2DEG) with the electric field in the linear region of the AlGaN/AlN/GaN HFET IV output characteristics, it is found that the polarization Coulomb field scattering still plays an important role in the electron mobility of AlGaN/AlN/GaN HFETs at the higher drain voltage and channel electric field. As drain voltage and channel electric field increase, the 2DEG density reduces and the polarization Coulomb field scattering increases, as a result, the 2DEG electron mobility decreases.

067303

We present a metallic/dielectric multi-wedge model to investigate the coupled edge plasmon modes (CEPMs), where all wedges have a common edge and the wave propagates along the edge direction. A general theoretical method valid to many practical structures is presented. The analytical dispersion relations of CEPMs in these structures are obtained and the CEPM properties are discussed with numerical results and the dispersion relations. For all structures mentioned in this paper, we find that the structures containing an even number of metallic wedges have four CEPMs and those with an odd-number of metallic wedges have two CEPMs. Further, the periodic structures containing any odd number of periods and any even number of periods possess their common CEPMs, respectively.

067306

, , , , and

A novel low specific on-resistance (Ron,sp) silicon-on-insulator (SOI) p-channel lateral double-diffused metal—oxide semiconductor (pLDMOS) compatible with high voltage (HV) n-channel LDMOS (nLDMOS) is proposed. The pLDMOS is built in the N-type SOI layer with a buried P-type layer acting as a current conduction path in the on-state (BP SOI pLDMOS). Its superior compatibility with the HV nLDMOS and low voltage (LV) complementary metal-oxide semiconductor (CMOS) circuitry which are formed on the N-SOI layer can be obtained. In the off-state the P-buried layer built in the N-SOI layer causes multiple depletion and electric field reshaping, leading to an enhanced (reduced) surface field (RESURF) effect. The proposed BP SOI pLDMOS achieves not only an improved breakdown voltage (BV) but also a significantly reduced Ron,sp. The BV of the BP SOI pLDMOS increases to 319 V from 215 V of the conventional SOI pLDMOS at the same half cell pitch of 25 μm, and Ron,sp decreases from 157 mΩ·cm2 to 55 mΩ·cm2. Compared with the PW SOI pLDMOS, the BP SOI pLDMOS also reduces the Ron,sp by 34% with almost the same BV.

067502

and

In this work we study the correlation function of the ground state of a two-dimensional fully frustrated Ising model as well as spin glass. The Pfaffian method is used to calculate free energy and entropy as well as the correlation function. We estimate the exponent of spin correlation function for the fully frustrated model and spin glass. In this paper an overview of the latest results on the spin correlation function is presented.

067503

, , , , , and

Zn0.99Cu0.01O films were studied experimentally and theoretically. The films were prepared by pulsed-laser deposition on Pt(111)/Ti/SiO2/Si substrates under various oxygen pressures to investigate the growth-dependence of the ferromagnetic properties. The structural, magnetic, and optical properties were studied, and it was found that all the samples possess a typical wurtzite structure, and that the films exhibit room-temperature ferromagnetism. The sample deposited at 600 °C and an oxygen pressure of 10 Pa showed a large saturation magnetization of 0.83 μB/Cu. The enhanced ferromagnetism in the (Cu, Li)-codoped ZnO is attributable to the existence of Zn vacancies (VZn), as shown by first-principles calculations. The photoluminescence analysis demonstrated the existence of VZn in both Zn0.99Cu0.01O and (Cu, Li)-codoped ZnO thin films, and this plays an important role in the increase of ferromagnetism, according to the results of first-principles calculations.

067504

, , , and

Interfacial magnetic anisotropy in a Pt/Co1−xFex/Pt multilayer is tuned by doping iron atoms into the cobalt layer. The perpendicular magnetic anisotropy and out-of-plane coercivity are found to decrease with increasing x. For a specific x, the out-of-plane coercivity acquires a maximal value as a function of the thickness of the CoFe layer. At low temperature, the coercivity is enhanced. Small coercivity but reasonably large perpendicular magnetic anisotropy can be obtained by controlling the x and CoFe layer thickness.

067701

, , , , , , , and

Aluminum—oxide films deposited as gate dielectrics on germanium (Ge) by atomic layer deposition were post oxidized in an ozone atmosphere. No additional interfacial layer was detected by the high-resolution cross-sectional transmission electron microscopy and X-ray photoelectron spectroscopy measurements made after the ozone post oxidation (OPO) treatment. Decreases in the equivalent oxide thickness of the OPO-treated Al2O3/Ge MOS capacitors were confirmed. Furthermore, a continuous decrease in the gate leakage current was achieved with increasing OPO treatment time. The results can be attributed to the film quality having been improved by the OPO treatment.

067702

, and

ZrO2 nanocrystallite-based charge trap flash memory capacitors incorporating a (ZrO2)0.6(SiO2)0.4 pseudobinary high-k oxide film as the charge trapping layer were prepared and investigated. The precipitation reaction in the charge trapping layer, forming ZrO2 nanocrystallites during rapid thermal annealing, was investigated by transmission electron microscopy and X-ray diffraction. It was observed that a ZrO2 nanocrystallite-based memory capacitor after post-annealing at 850 °C for 60 s exhibits a maximum memory window of about 6.8 V, good endurance and a low charge loss of ∼25% over a period of 10 years (determined by extrapolating the charge loss curve measured experimentally), even at 85 °C. Such 850 °C-annealed memory capacitors appear to be candidates for future nonvolatile flash memory device applications.

067801

, , , and

We propose an ultrathin wide-band metamaterial absorber (MA) based on a Minkowski (MIK) fractal frequency selective surface and resistive film. This absorber consists of a periodic arrangement of dielectric substrates sandwiched with an MIK fractal loop structure electric resonator and a resistive film. The finite element method is used to simulate and analyze the absorption of the MA. Compared with the MA-backed copper film, the designed MA-backed resistive film exhibits an absorption of 90% at a frequency region of 2 GHz–20 GHz. The power loss density distribution of the MA is further illustrated to explain the mechanism of the proposed MA. Simulated absorptions at different incidence cases indicate that this absorber is polarization-insensitive and wide-angled. Finally, further simulated results indicate that the surface resistance of the resistive film and the dielectric constant of the substrate can affect the absorbing property of the MA. This absorber may be used in many military fields.

067802

and

In this paper, we study a three-resonant metamaterial with the combination of dual-resonant and single-resonant metamaterials. We present a new method to design multi-resonant metamaterial, which has a smaller dimension than general symmetric and asymmetric multi-resonant metamaterials. Theoretical and experimental results show that the structure has three distinct absorption frequencies centering around 0.29 THz, 0.46 THz, and 0.92 THz, and that each of them corresponds to a different resonant mode. Due to the good separation of the different resonances, this design provides a unique and effective method to construct multiband terahertz devices.

067803

, , and

Based on space-charge wave theory, the formulae of the beam—wave coupling coefficient and the beam-loaded conductance are given for the beam—wave interaction in an N-gap Hughes-type coupled cavity chain. The ratio of the non-beam-loaded quality factor of the coupled cavity chain to the beam quality factor is used to determine the stability of the beam—wave interaction. As an example, the stabilities of the beam—wave interaction in a three-gap Hughes-type coupled cavity chain are discussed with the formulae and the CST code for the operations of the 2π, π, and π/2 modes, respectively. The results show that stable operation of the 2π, π, and π/2 modes may all be realized in an extended-interaction klystron with the three-gap Hughes-type coupled cavity chain.

Cross-disciplinary physics and related areas of science and technology

068102

, , , , , and

Different aluminum-doped ZnO (AZO)/metal composite thin films, including AZO/Ag/Al, AZO/Ag/nickel—chromium alloy (NiCr), and AZO/Ag/NiCr/Al, are utilized as the back reflectors of p—i—n amorphous silicon germanium thin film solar cells. NiCr is used as diffusion barrier layer between Ag and Al to prevent mutual diffusion, which increases the short circuit current density of solar cell. NiCr and NiCr/Al layers are used as protective layers of Ag layer against oxidation and sulfurization, the higher efficiency of solar cell is achieved. The experimental results show that the performance of a-SiGe solar cell with AZO/Ag/NiCr/Al back reflector is best. The initial conversion efficiency is achieved to be 8.05%.

068103

, , , and

For both the vibrating and steady supporting surfaces of a scanning disk in a Besocke-style piezoelectric scanner, a theoretical model is given by considering the nonlinear lateral friction at the micro-contact interface between the positioning legs and the supporting surface. Numerical simulations demonstrate that unexpected flexural vibrations can arise from a vibrating ramp, and their frequencies are lower than the eigenfrequencies of the scanner in the linearly elastic regime. The vibrations essentially depend on 1) the vibrational states of the supporting ramp and the steel ball tips on the three piezoelectric positioning legs, and 2) the tribological characteristics of the contacts between the tips and the ramp. The results give an insight into the intrinsic vibrations of the scanners, and are applicable in designing and optimizing piezoelectric scanning systems.

068201

, , and

The quasi-classical trajectory (QCT) is calculated to study the stereodynamics properties of the title reaction H(2S) + NH(X3 Σ) → N(4S) + H2 on the ground state 4A'' potential energy surface (PES) constructed by Zhai and Han [2011 J. Chem. Phys.135 104314]. The calculated QCT reaction probabilities and cross sections are in good agreement with the previous theoretical results. The effects of the collision energy on the kk' distribution and the product polarization of H2 are studied in detail. It is found that the scattering direction of the product is strongly dependent on the collision energy. With the increase in the collision energy, the scattering directions of the products change from backward scattering to forward scattering. The distribution of Pr) is strongly dependent on the collision energy below the lower collision energy (about 11.53 kcal/mol). In addition, the P(φr) distribution dramatically changes as the collision energy increases. The calculated QCT results indicate that the collision energy plays an important role in determining the stereodynamics of the title reaction.

068202

, , and

A new London—Eyring—Polanyi—Sato potential energy surface is employed in this work to study the stereo properties of the O(3P) + CH4 → H + CH3O reaction in its rovibrationally ground state using the quasiclassical trajectory method (QCT). Our calculations are performed at a range of collision energies, Ec = 1.5 eV∼ 3.5 eV, and the excitation function obtained by the QCT method accords well with the experimental data. The product rotational polarization is calculated, and the product shows a strong rotational polarization in the centre-of-mass coordinate system. The orientation of the product rotational angular momenta is sensitive to the increase in collision energy, and the alignment of the product rotational angular momenta shows some of the properties of the heavy heavy—light mass combination reactions. In the isotopic substituted reaction study, when the H atoms in methane are replaced by D atoms, the rotational polarization is obviously reduced. The polarization-dependent differential cross section is also studied by this QCT calculation to provide detailed information about the rotational alignment and orientation of the product.

068203

, , , , and

The quasi-classical trajectory (QCT) method is employed to calculate the stereodynamics of the abstraction reactions H/D + HS/DS based on an accurate potential energy surface [Lü S J, Zhang P Y, Han K L and He G Z 2012 J. Chem. Phys. 136 094308]. The reaction cross sections of the title reaction are computed, and the vector correlations for different collision energies and different initial vibrational states are presented. The influences of the collision energy and reagent vibration on the product polarization are studied, and the product polarizations of the title reactions are found to be distinctly different, which arises from the different mass factors, collision energies, and reagent vibrational states.

068204

, , , and

WO3 thin films were sputtered onto alumina substrates by DC facing-target magnetron sputtering. One sample was rapid-thermal-annealed (RTA) at 600 °C in a gas mixture of N2:O2 = 4:1, and as a comparison, another was conventionally thermal-annealed at 600 °C in air. The morphology of both was investigated by scanning electron microscopy (SEM) and atomic force microscopy (AFM), and the crystallization structure and phase identification were characterized by X-ray diffraction (XRD). The NO2-sensing measurements were taken under LED light at room temperature. The sensitivity of the RTA-treated sample was found to be high, up to nearly 100, whereas the sensitivity of the conventionally thermal-annealed sample was about five under the same conditions. From the much better selectivity and response-recovery characteristics, it can be concluded that compared to conventional thermal annealing, RTA has a greater effect on the NO2-sensing properties of WO3 thin films.

068401

, , and

A flux-controlled memristor characterized by smooth cubic nonlinearity is taken as an example, upon which the voltage—current relationships (VCRs) between two parallel memristive circuits — a parallel memristor and capacitor circuit (the parallel MC circuit), and a parallel memristor and inductor circuit (the parallel ML circuit) — are investigated. The results indicate that the VCR between these two parallel memristive circuits is closely related to the circuit parameters, and the frequency and amplitude of the sinusoidal voltage stimulus. An equivalent circuit model of the memristor is built, upon which the circuit simulations and experimental measurements of both the parallel MC circuit and the parallel ML circuit are performed, and the results verify the theoretical analysis results.

068402

In this paper, the electrical parameters of Au/n-Si (MS) and Au/Si3N4/n-Si (MIS) Schottky diodes are obtained from the forward bias current—voltage (I—V) and capacitance—voltage (C—V) measurements at room temperature. Experimental results show that the rectifying ratios of the MS and MIS diodes at ± 5 V are found to be 1.25 × 103 and 1.27 × 104, respectively. The main electrical parameters of the MS and MIS diodes, such as the zero-bias barrier height (ΦBo) and ideality factor (n), are calculated to be 0.51 eV (IV), 0.53 eV (CV), and 4.43, and 0.65 eV (IV), 0.70 eV (C—V), and 3.44, respectively. In addition, the energy density distribution profile of the interface states (Nss) is obtained from the forward bias I—V, and the series resistance (Rs) values for the two diodes are calculated from Cheung's method and Ohm's law.

068501

, , , and

A novel super-junction lateral double-diffused metal—oxide semiconductor (SJ-LDMOS) with a partial lightly doped P pillar (PD) is proposed. Firstly, the reduction in the partial P pillar charges ensures the charge balance and suppresses the substrate-assisted depletion effect. Secondly, the new electric field peak produced by the P/P junction modulates the surface electric field distribution. Both of these result in a high breakdown voltage (BV). In addition, due to the same conduction paths, the specific on-resistance (Ron,sp) of the PD SJ-LDMOS is approximately identical to the conventional SJ-LDMOS. Simulation results indicate that the average value of the surface lateral electric field of the PD SJ-LDMOS reaches 20 V/μm at a 15 μm drift length, resulting in a BV of 300 V.

068502

, , , , and

In the present paper we study the influences of the bias voltage and the external components on the damage progress of a bipolar transistor induced by high-power microwaves. The mechanism is presented by analyzing the variation in the internal distribution of the temperature in the device. The findings show that the device becomes less vulnerable to damage with an increase in bias voltage. Both the series diode at the base and the relatively low series resistance at the emitter, Re, can obviously prolong the burnout time of the device. However, Re will aid damage to the device when the value is sufficiently high due to the fact that the highest hot spot shifts from the base-emitter junction to the base region. Moreover, the series resistance at the base Rb will weaken the capability of the device to withstand microwave damage.

068503

, , , , and

AlGaN/GaN high-electron-mobility transistors (HEMTs) with Al-doped ZnO (AZO) transparent gate electrodes are fabricated, and Ni/Au/Ni-gated HEMTs are produced in comparison. The AZO-gated HEMTs show good DC characteristics and Schottky rectifying characteristics, and the gate electrodes achieve excellent transparencies. Compared with Ni/Au/Ni-gated HEMTs, AZO-gated HEMTs show a low saturation current, high threshold voltage, high Schottky barrier height, and low gate reverse leakage current. Due to the higher gate resistivity, AZO-gated HEMTs exhibit a current—gain cutoff frequency (fT) of 10 GHz and a power gain cutoff frequency (fmax) of 5 GHz, and lower maximum oscillation frequency than Ni/Au/Ni-gated HEMTs. Moreover, the C—V characteristics are measured and the gate interface characteristics of the AZO-gated devices are investigated by a C—V dual sweep.

068504

In this paper, a vacuum system is employed to compare the emission stabilities of the same ZnO cathode in a sealed field emission (FE) device and under ultrahigh vacuum (UHV) conditions. It is observed that the emission current is more stable under the UHV level than in the device. When all conditions except the ambient gases are kept unchanged, the emission current degradation is mainly caused by the residual gases in the sealed device. The quadrupole mass spectrometer (QMS) equipped on the vacuum system is used to investigate the residual gas components. Based on the obtained QMS data, the following conclusions can be drawn: the residual gases in ZnO-FE devices are H2, CH4, CO, Ar, and CO2. These residual gases can change the work function at the surface through adsorption or ion bombardment, thereby degrading the emission current of the cathode.

068505

, , , , , , , , and

GaN/InGaN superlattice barriers are used in InGaN-based light-emitting diodes (LEDs). The electrostatic field in the quantum wells, electron hole wavefunction overlap, carrier concentration, spontaneous emission spectrum, light-current performance curve, and internal quantum efficiency are numerically investigated using the APSYS simulation software. It is found that the structure with GaN/InGaN superlattice barriers shows improved light output power, and lower current leakage and efficiency droop. According to our numerical simulation and analysis, these improvements in the electrical and optical characteristics are mainly attributed to the alleviation of the electrostatic field in the active region.

068506

, , , and

Compositing soft and hard materials is a promising method to decrease the coercivity of L10 FePt, which is considered to be a suitable material for bit-patterned media. This paper reports the simulation of three models of FeCo/L10 FePt exchange-coupled composite particles for bit patterned media by the OOMMF micromagnetic simulation software: the enclosed model, the side-enclosed model, and the top-covered model. All of them have the same volumes of the soft and hard parts but different shapes. Simulation results show that the switching fields for the three models can be reduced to about 10 kOe (1 Oe = 79.5775 A/m) and the factor gain can be improved to 1.4 when the interface exchange coefficient has a proper value. Compared to the other models, the enclosed model has a wider range of interface exchange coefficient values, in which a low switching field and high gain can be obtained. The dependence of the switching fields on the angle of the applied field shows that none of the three models are easily affected by the stray field of a magnetic head.

068801

, , , , and

The combustion characteristics of styrene—butadiene—styrene (SBS) asphalt are studied by thermogravimetric analysis (TG/DTG) at four different heating rates. According to the saturates/aromatics/resins/asphaltenes (SARA) fractionation method, the combustion process of SBS asphalt can be divided by Gaussian peak fitting into three main stages: oil content release, resin pyrolysis, and asphaltene and char combustion. When the heating rate increases, the mass losses of the oil content and resin pyrolysis increase, and less asphaltenes are formed at a higher temperature. The activation energy values are calculated by the Coats—Redfern method to be in the range 61.6 kJ/mol–142.9 kJ/mol. The Popescu method is used for the kinetic analysis, and the result shows that the three stages of asphalt combustion can be explained by the sphere phase boundary reaction model, the second order chemical reaction model, nucleation, and its subsequent growth model, respectively.

068802

, , , , , , , and

InGaN/GaN p—i—n solar cells, each with an undoped In0.12Ga0.88N absorption layer, are grown on c-plane sapphire substrates by metal—organic chemical vapor deposition. The effects of the thickness and dislocation density of the absorption layer on the collection efficiency of InGaN-based solar cells are analyzed, and the experimental results demonstrate that the thickness of the InGaN layer and the dislocation density significantly affect the performance. An optimized InGaN-based solar cell with a peak external quantum efficiency of 57% at a wavelength of 371 nm is reported. The full width at half maximum of the rocking curve of the (0002) InGaN layer is 180 arcsec.

068901

, and

For further research on the gravity mechanism of the routing protocol in complex networks, we introduce the concept of routing awareness depth, which is represented by ρ. On this basis, we define the calculation formula of the gravity of the transmission route for the packet, and propose a routing strategy based on the gravitational field of the node and the routing awareness depth. In order to characterize the efficiency of the method, we introduce an order parameter, ζ, to measure the throughput of the network by the critical value of phase transition from free flow to congestion, and use the node betweenness centrality, B, to test the transmission efficiency of the network and congestion distribution. We simulate the network transmission performance under different values of the routing awareness depth, ρ. Simulation results show that if the value of the routing awareness depth ρ is too small, then the gravity of the route is composed of the attraction of very few nodes on the route, which cannot improve the capacity of the network effectively. If the value of the routing awareness depth ρ is greater than the network's average distance 〈l〉, then the capacity of the network may be improved greatly and no longer change with the sustainable increment of routing awareness depth ρ, and the routing strategy performance enters into a constant state. Moreover, whatever the value of the routing awareness depth ρ, our algorithm always effectively balances the distribution of the betweenness centrality and realizes equal distribution of the network load.