Acta Physica Sinica
Citation Search Quick Search

ISSN 1000-3290
CN 11-1958/O4
Semimonthly
About
   » About APS
   » Editorial Board
   » SCI IF
   » Staff
   » Contact
Browse APS
   » Accepts
   » In Press
   » Current Issue
   » Past Issues
   » View by Fields
   » Top Downloaded
   » Sci Top Cited
Authors
   » Submit an Article
   » Manuscript Tracking
   » Call for Papers
   » Scope
   » Instruction for Authors
   » Copyright Agreement
   » Templates
   » Author FAQs
   » PACS
Referees
   » Review Policy
   » Referee Login
   » Referee FAQs
   » Editor in Chief Login
   » Office Login
Links
   »
HighLights More»   
Glass transition in binary mixture of colloidal ellipsoids and spheres
Sun Yan-Li, Wang Hua-Guang, Zhang Ze-Xin
Acta Physica Sinica, 2018, 67 (10): 106401
Piezo-electrochemical coupling of AgNbO3 piezoelectric nanomaterials
Hong Yuan-Ting, Ma Jiang-Ping, Wu Zheng, Ying Jing-Shi, You Hui-Lin, Jia Yan-Min
Acta Physica Sinica, 2018, 67 (10): 107702
Molecular dynamics simulations on DNA flexibility: a comparative study of Amber bsc1 and bsc0 force fields
Xiong Kai-Xin, Xi Kun, Bao Lei, Zhang Zhong-Liang, Tan Zhi-Jie
Acta Physica Sinica, 2018, 67 (10): 108701
Current Issue Accepts In Press Earlier Issues Top Downloaded SCI Top Cited
  Acta Physica Sinica--2018, 67 (10)   Published: 20 May 2018
Select | Export to EndNote
CONDENSED MATTER: STRUCTURAL, MECHANICAL, AND THERMAL PROPERTIES

Glass transition in binary mixture of colloidal ellipsoids and spheres Hot!

Sun Yan-Li, Wang Hua-Guang, Zhang Ze-Xin
Acta Physica Sinica. 2018, 67 (10): 106401 doi: 10.7498/aps.67.20180264
Full Text: [PDF 1546 KB] Download:(26)
Show Abstract

The nature of glass and glass transition are considered to be one of the most fundamental research problems in condensed matter physics. Colloidal suspension provides a novel model system for studying glass and glass transition, since the structures and dynamics of a colloidal system can be quantitatively probed by video microscopy. Traditional systems for studying glass transition typically are single-component systems composed of either isotropic or anisotropic colloidal particles. Recently, glass transition of mixture of isotropic and anisotropic colloids has attracted great attention, such as the observation of rotational glass and translational glass, and the establishment of the two-step glass transition. Similarly, computer simulations have also shown that mixture of isotropic and anisotropic colloidal particles could manifest interesting, new glassy behaviors. However, the experimental study of the glass transition in such a colloidal mixture is still rare. In this paper, we experimentally investigate the glass transition of a binary mixture of colloidal ellipsoids and spheres. The colloidal spheres are polystyrene microspheres with a diameter of 1.6 μm, and the ellipsoids are prepared by physically stretching from polystyrene microspheres of 2.5 μm in diameter. The major and minor axes of the as-prepared ellipsoid are 2.0 μm and 1.2 μm, respectively. The mixture is confined between two glass slides to make a quasi-two-dimensional sample. To prevent the mixture from crystallizing, the mixing ratio of ellipsoids and spheres is chosen to be 1/4 in number, which is similar to the mixing ratio used in the classical Kob-Anderson model of binary sphere mixture. We systemically increase the area fraction of colloidal mixture to drive the glass transition. We then employ bright-field video microscopy to record the motion of the particles in the colloidal suspension at a single particle level, and the trajectories of individual particles are obtained by standard particle tracking algorithm. Through the analysis of radial distribution function, Voronoi diagram and local order parameter, we find that the ellipsoids can effectively inhibit the spheres from crystalizing, and the structure of the system remains disordered when increasing the area fraction. For dynamics, mean square displacement and self-intermediate scattering function are calculated. We find that the dynamic process of the system slows down substantially when increasing the area fraction, and the relaxation time of the system increases rapidly and diverges close to the glass transition point predicted by the mode coupling theory. Moreover, we analyze the fast particles that participate in cooperative rearrangement regions (CRRs) in the system, and find that the shapes, sizes and positions of CRRs are closely related to the locations of the ellipsoids in the system.

Abnormal viscosity changes in high-temperature metallic melts

Shang Ji-Xiang, Zhao Yun-Bo, Hu Li-Na
Acta Physica Sinica. 2018, 67 (10): 106402 doi: 10.7498/aps.67.20172721
Full Text: [PDF 1735 KB] Download:(15)
Show Abstract
The viscosity of high-temperature metallic melt, which is an important index for evaluating dynamics of liquid melt, is one of the basic physical properties. It not only influences the mold-filling capacity of melting metal in traditional casting techniques, but also exhibits more distinct influence on the fabrication of advanced material, such as metallic glass. According to the variation tendency of viscosity with temperature in alloy melt, the fragility of superheated melt could be obtained, which has proved to correlate with the ability of alloy to form glass. Besides, the viscosity of alloy well above the liquidus temperature also plays a key role in probing into the characteristic of liquid-liquid phase transition, the fragile-to-strong transition phenomenon, how the potential energy landscape evolves during cooling, etc. It has been generally accepted that the viscosity of metallic melt at high temperatures increases with temperature decreasing and could be fitted by an Arrhenius curve in the whole temperature range. However, recently more and more studies show that the viscosity of metallic melt cannot be fitted by only one Arrhenius curve. Instead, there exists at least one specific temperature below which the viscosity data begins to deviate from the Arrhenius curve at high temperature during cooling. These data could be described by another Arrhenius curve. In order to in depth understand this phenomenon, in this paper we summarize the viscosity data of different metallic melts in the literature. On the basis of introducing the method of detecting high-temperature melt viscosity, we discuss comprehensively the changing tendency of viscosity with temperature and the characteristics of abnormal viscosity changes in pure metal, binary and multivariate alloys well above the liquidus temperature. It is found that the abnormal viscosity changes generally occur in alloys that could form the types of intermetallic compounds. The abnormal viscosity change in metallic melt is accompanied with exothermic or endothermic effect, depending on alloy system, and reflects the existence of liquid-liquid transition well above the liquidus temperature. Besides, such an abnormal change of viscosity influences the ability to form metallic glass liquids. Although the abnormal dynamic change of metallic melt hints the existence of complexity of structural change in liquid during cooling, what is the key factor underlying this phenomenon remains a mystery. By combining the advanced experimental techniques such as high-energy X-ray diffraction and neutron scattering with the computer simulation method, this problem may be understood further. Besides, the relation between viscosity abnormity and the phase diagram is another problem that deserves to be noticed in the future.

Effects of organic cations on performance of halide perovskite solar cell

Zhang Ao, Chen Yun-Lin, Yan Jun, Zhang Chun-Xiu
Acta Physica Sinica. 2018, 67 (10): 106701 doi: 10.7498/aps.67.20180236
Full Text: [PDF 3178 KB] Download:(75)
Show Abstract

The halide perovskite solar cells employing CH3NH3PbX3 (X=Cl-, Br-, I-) and CH3NH3PbI3-xClx as light absorbers each have shown a rapid rise in power conversion efficiency (PCE) from 3.8% to 22.1% in recent years. The excellent photovoltaic performance is attributed to good optical and electrical properties such as appropriate bandgap, large absorption coefficient, high carrier mobility, long carrier lifetime and long carrier diffusion length. However, the physical mechanism of high PCE for halide perovskite solar cells is still unclear. The Gaussian 09 software is utilized to optimize the geometries of isolated CH3NH3+ and CH3NH3 at a B3 LYP/6-311++G(d, p) level, and the Multiwfn software is used to visualize the electrostatic potentials (ESPs) of CH3NH3+ and CH3NH3. Based on the ESPs of CH3NH3+ and CH3NH3, it is found that the CH3NH3+ has a strong electrophilic character, however, the NH3- side and CH3- side of CH3NH3 have weak nucleophilic and electrophilic character, respectively. So the electrostatic characteristics of CH3NH3+ and CH3NH3 are significantly different. The strong electrostatic repulsive interaction between two neighboring CH3NH3+ radicals plays an important role in structural phase transition of CH3NH3PbI3 material. At room temperature, the CH3NH3+ in the inorganic cage is activated and disordered, and has a strong electrophilic character. Due to these characteristics of CH3NH3+, the interfacial electrons at TiO2/CH3NH3PbI3 heterojunction are combined with CH3NH3+ to form CH3NH3 in the inorganic[PbI3]- framework. The CH3NH3 at the heterojunction under the built-in electric field is more easily oriented than CH3NH3+. Two initial geometrical configurations for CH3NH3+:CH3NH3 and CH3NH3:CH3NH3 dimers are optimized by using Gaussian 09 at an MP2/Aug-cc-PVTZ level. On the basis of the electrostatic characteristic of CH3NH3+:CH3NH3 dimer, the interfacial electrons at TiO2/CH3NH3PbI3 heterojunction are easily injected into the CH3NH3PbI3 material, which leads to the strong polarization of CH3NH3PbI3 material at the heterojunction. From the ESP of optimized CH3NH3:CH3NH3 dimer, it is found that the weak electrostatic field of the inorganic framework, parallel to C-N axis, is induced by the CH3NH3 orientational order, which is made for improving the photogenerated electron-hole pair separation and carrier transport. The TiO2/CH3NH3PbI3 heterojunction has more advantage than traditional p-n junction because of no consumption of carrier for CH3NH3PbI3 material in the process of forming built-in electric field. The physical mechanism is the origin of high PCE for CH3NH3PbI3 solar cells. According to the experimental results and first-principle calculations, we can draw an important conclusion that the electrostatic characteristics of organic CH3NH3+ cations in the inorganic[PbI3]- framework result in the high performances of halide perovskite solar cells.

Electronic structure and spin/valley transport properties of monolayer MoS2 under the irradiation of the off-resonant circularly polarized light

Zhang Xin-Cheng, Liao Wen-Hu, Zuo Min
Acta Physica Sinica. 2018, 67 (10): 107101 doi: 10.7498/aps.67.20180213
Full Text: [PDF 436 KB] Download:(22)
Show Abstract
The new-type monolayer semiconductor material molybdenum disulfide (MoS2) is direct band gap semiconductor with a similar geometrical structure to graphene, and as it owns superior physical features such as spin/valley Hall effect, it should be more excellent than graphene from the viewpoint of device design and applications. The manipulation of the spin and valley transport in MoS2-based device has been an interesting subject in both experimental and theoretical researches. Experimentally, the photoninduced quantum spin and valley Hall effects may result in high on-off speed spin and/or valley switching based on MoS2. Theoretically, the off-resonant electromagnetic field induced Floquet effective energy should modulate effectively the electronic structure, spin/valley Hall conductance as well as the spin/valley polarization of the MoS2, through the virtual photon absorption and/or emission processes. Utilizing the low energy effective Hamilton model from the tight-binding approximation and Kubo linear response theorem, we theoretically investigate the electronic structure and spin/valley transport properties of the monolayer MoS2 under the irradiation of the off-resonant circularly polarized light in the present work. The band gaps around the K and K' point of the Brillouin region for monolayer MoS2 proves to increase linearly and decrease firstly and then increase, respectively with the increase of external off-resonant right-circularly polarized light induced effective coupling energy, and decrease firstly and then increase and increase linearly with the increase of left-circularly polarized light induced effective coupling energy, therefore, the interesting transition of semiconducting-semimetallic-semiconducting may be observable in monolayer MoS2. Furthermore, the spin and valley Hall conductance of the monolayer MoS2 for the case without off-resonant circularly polarized light are 0 and 2e2/h, respectively, and they will convert into -2e2/h and 0 when the absolute value of the off-resonant circularly polarized light induced effective coupling energy is in a range of 0.79-0.87 eV. Finally, the spin polarization for monolayer MoS2 increases up to a largest value and changes from positive to negative and/or negative to positive at the vicinity of the effective coupling energy ±0.79 eV of the off-resonant right/left circularly polarized light, while the valley polarization should increase firstly and then decrease with the off-resonant circularly polarized light, and goes up to 100% in the range of 0.79-0.87 eV of the absolute value for effective coupling energy. Therefore, the external off-resonant circularly polarized electromagnetic field should be an effective means in manipulating the electronic structure, spin/valley Hall conductance and spin/valley polarization of the monolayer MoS2, the two-dimensional MoS2 may be tuned into a brand bandgap material with excellent spin/valley and optoelectrical properties.

First-principles study on the optical properties of Fe-doped GaN

Jia Wan-Li, Zhou Miao, Wang Xin-Mei, Ji Wei-Li
Acta Physica Sinica. 2018, 67 (10): 107102 doi: 10.7498/aps.67.20172290
Full Text: [PDF 838 KB] Download:(19)
Show Abstract
Using hybrid density functional theory, we investigate the structural, electronic and optical properties of pristine GaN and Fe-doped GaN with a Fe concentration of 12.5%. Specifically, we first analyze the crystal lattice constant, band structure, and density of states, respectively. Then we predict the dielectric function, absorption coefficient, refractive index, reflectivity, energy-loss spectrum and extinction coefficient. Finally, we analyze the influences of the doping of Fe element on the photoelectric property of Fe doped systems. The calculated lattice constants for perfect GaN are a=b=3.19 Å, c=5.18 Å, which are in good agreement with the experimental values. Furthermore, we find that the doping of Fe element has little effect on the structural properties of GaN. The Band gap of pristine GaN is predicted to be 3.41 eV, very close to the experimental value of 3.39 eV. The band gap of Fe doped GaN (12.5%) significantly decreases to 3.06 eV. By comparing the densities of states of the systems with and without Fe doping, it is found that Fe-3 d state is mainly responsible for the decrease of band gap. The calculated static dielectric constant of perfect GaN is 5.74, and it increases to 6.20 after incorporating the Fe element. The results about the imaginary part of dielectric function show that two equal-strength perfect GaN peaks are observed to be at 6.81 eV and 10.85 eV. The first peak is closely related to the direction transition from the valence band top to the conduction band bottom. Furthermore, it is also observed that a peak is located at 4.04 eV in the low energy, which can be understood as resulting from the electron transition inside the valence band. The optical absorption edge of the intrinsic GaN is 3.25 eV, corresponding to the transition energy. The reason why this energy is smaller than the bandgap is because the electronic band gap equals the sum of optical bandgap and exciton energy. However, the maximum absorption coefficients of these two systems both occur at 13.80 eV in energy. The refractive index for intrinsic system is 2.39, and it increases to 2.48 after doping the Fe element. It is found from the energy-loss spectrum that the maximum energy-loss is at 20.02 eV for a perfect system, while it is at 18.96 eV for a doped system. Additionally, we obtain the reliable reflectivity and excitation coefficient. In conclusion, our calculated results provide a well theoretical basis for the theoretical research on the co-doping of Fe element and other elements. The analyses on the Fe-doped GaN high-voltage photoconductive switch materials and devices also provide a powerful theoretical basis and experimental support in the future research.

First-principles study of hydrogen storage properties of silicene under different Li adsorption components

Sheng Zhe, Dai Xian-Ying, Miao Dong-Ming, Wu Shu-Jing, Zhao Tian-Long, Hao Yue
Acta Physica Sinica. 2018, 67 (10): 107103 doi: 10.7498/aps.67.20172720
Full Text: [PDF 645 KB] Download:(10)
Show Abstract
Alkali metal has predicted to be a promising candidate for decorating silicene surface to obtain the high hydrogen storage capacity, owing to their physical properties of lightweight, lower cohesive energy, and appropriate strength of the interaction with H2 molecules. However, though the high potential in hydrogen storage of alkali metal adatoms-decorated silicene under the fixed adatom adsorption component is well known, the evidence for the hydrogen storage capacity of alkali metal adatoms-decorated silicene under different adatom adsorption components remains largely unexplored, which may be of great significance to make the most advantages of alkali metal adatoms-decorated silicene in hydrogen storage aspects. Herein, according to the first-principles calculation corrected by the van der Waals effect, we take Li-decorated silicene for example and perform the detailed study of the geometry structure, the stability and the hydrogen storage capacity of silicene under different Li adsorption components (LixSi1-x), aiming to maximize the hydrogen storage performance of Li-decorated silicene. The results show that the preferred site of Li changes from the hollow site to the valley site as the Li component increases from 0.11 to 0.50, and binding energy of Li is always greater than the corresponding cohesive energy, showing the high stability of Li-decorated silicene and the feasibility of the method to obtain a higher hydrogen storage capacity by increasing the Li component. The hydrogen storage of silicene under different Li adsorption components is investigated by the sequential addition of H2 molecules nearby Li atoms in a stepwise manner. It can be observed that the hydrogen storage capacity of Li-decorated silicene increases and the average adsorption energy decreases with the increase of the Li component. The corresponding hydrogen storage capacities of Li0.11Si0.89, Li0.20Si0.80, Li0.33Si0.67, Li0.43Si0.57 can reach up to 2.54 wt%, 4.82 wt%, 6.00 wt% and 9.58 wt% with 0.58 eV/H2, 0.47 eV/H2, 0.54 eV/H2 and 0.41 eV/H2 average adsorption energy, respectively. When the Li component increases up to 0.50, Li atoms are saturated with a maximum hydrogen storage capacity of 11.46 wt% and an average adsorption energy of 0.34 eV/H2, which well meet the hydrogen storage standard set by the U.S. Department of Energy and mean that the hydrogen storage can be theoretically improved by increasing the Li adsorption component to a saturated level. Furthermore, we analyze the Mulliken charge population, the charge density difference and the density of states, showing that the charge-induced electrostatic interaction and the orbital hybridization are the key factors for the hydrogen adsorption of Li-decorated silicene. Our results may enhance our fundamental understanding of the hydrogen storage mechanism and explore the applications in areas of hydrogen storage for Li-decorated silicene, which are of great importance for the usage of hydrogen in the future.

Enhanced optical transmission by exciting hybrid states of Tamm and surface plasmon polaritons in single slit with multi-pair groove nanostructure

Qi Yun-Ping, Zhou Pei-Yang, Zhang Xue-Wei, Yan Chun-Man, Wang Xiang-Xian
Acta Physica Sinica. 2018, 67 (10): 107104 doi: 10.7498/aps.67.20180117
Full Text: [PDF 1912 KB] Download:(23)
Show Abstract
In recent years, a metallic single slit nanostructure or slit array structure, due to simple structure and easy-to integration, has been used to construct a light source in the nanostructures based on the surface plasmon polaritons (SPPs). However, the problem of low transmission through an isolated subwavelength single slit nanostructure is still existent. The main reason is that the excitation efficiency of SPPs in the single slit nanostructure is not too high. Therefore, how to effectively enhance the optical transmission has become a research focus. In order to further improve the transmittance of the metallic single slit nanostructure, in this paper, we improve the single slit nanostructure imbedded in the metal silver thin film on a distributed Bragg reflector (DBR) proposed in previous literature. As a result, a novel method of designing a single slit on a DBR is proposed to effectively enhance the optical transmission in a single slit by improving the excitation efficiency of SPPs. Our proposed novel structure is made up of a subwavelength single nano-slit surrounded symmetrically by a pair of grooves on both sides of metal silver film on a distributed Bragg reflector. When the TM polarized light is illuminated from the DBR side of our proposed structure to the DBR-silver slit-grooves nanostructure, the Tamm plasmon polaritons (TPPs) at the interface between the silver film and the DBR and the SPPs in the slit on the entrance side of the silver film are excited at the DBR-silver film interface, and the SPPs in the slit and grooves pair on the exit side of the silver film are excited simultaneously. In our proposed structure, coupling between the TPPs and the SPPs leads to the hybrid state of Tamm and surface plasmon polaritons in the slit and grooves. Finally, taking advantage of constructive interference between SPPs excited by the grooves and exciting hybrid states of TPPs-SPPs in the slit, due to the local field enhancement effect of the TPPs mode and the coupling effect of constructive interference between the pair grooves and the nano-slit, the excitation efficiency of the SPPs can be increased significantly. Furthermore, the quasi Fabry-Pérot resonance effect in the nano-slit is taken into consideration, and the transmittance of our proposed structure is enhanced greatly. In the present paper, the finite element method is used to study the transmission properties of the single nano-slit embedded with paired grooves on the DBR-sliver nanostructure. After a series of parameters are optimized, the maximum transmittance through the single slit in DBR-silver slit-groove nanostructure can increase to 0.22, and this transmittance is expected to be about 22 times the transmittance (0.01) of the light through a single slit in a silver film on the TiO2 substrate (without DBR and grooves), which is higher than the maximum light transsmission 0.166 given in Ref.[23]. The research results of this study have a certain application value in the fields of nano-light source design, photonic integrated circuits and optical signal transmission and so on.

Electronic structures and ferroelectric properties of Ba-doped ZnO

Xu Jia-Nan, Chen Huan-Ming, Pan Feng-Chun, Lin Xue-Ling, Ma Zhi, Chen Zhi-Peng
Acta Physica Sinica. 2018, 67 (10): 107701 doi: 10.7498/aps.67.20172263
Full Text: [PDF 1368 KB] Download:(19)
Show Abstract
Wurtzite ZnO has long been considered to be a promising candidate material for photovoltaic application due to its high power conversion efficiency. More interestingly and very recently, some research results suggested that the ferroelectric property of the photovoltaic material introduced by chemical elements doping can promote its power conversion efficiency significantly. Therefore, in order to understand the effect of Ba doping on the electronic structure and the ferroelectric properties of ZnO and to reveal the potentially optoelectronic properties of Zn1-xBaxO, the energy band structure, the density of states, and the polarizability and the relative dielectric constant of the bulk Ba-doped ZnO supercell system, in which the Zn atoms are partly and uniformly substituted by the Ba atoms, are investigated by using the first-principles method based on the density functional theory and other physical theory. The norm-conserving pseudopotentials and the plane-wave basis set with a cut-off energy of 600 eV are used in the calculation. The generalized gradient approximation refined by Perdew and Zunger (GGA-PBE), the local density approximation (LDA) and the local density approximation added Hubbard energy (LDA+U) are employed for determining the exchange-correlation energy respectively. Brillouin zone is set to be within 4×4×5K point mesh generated by the Monkhorst-Pack scheme. The self-consistent convergence of total energy is at 2.0×10-6 eV/atom. Additionally, in order to obtain a stable and accurate calculation result, the cell structure is optimized prior to calculation. The calculated results suggest that the bulk Ba-doped ZnO semiconductor system is still a semiconductor with a direct wide band gap. The band gap of Zn1-xBaxO increases gradually with Ba atom doping percentage increasing from 12.5% to 87.5%. Consequently, the ferroelectric polarization properties and the dielectric properties of the bulk Ba-doped wurtzite ZnO materials are tailored by doping Ba atoms. It indicates that the polarizability of Zn1-xBaxO system increases with Ba doping atomic percentage increasing, especially, the polarizability reaches to a maximum when the atomic percentage of doping is 75%. Meanwhile, the relative dielectric constant inversely decreases with Ba atomic percentage increasing. This is attributed to the effective contribution of Ba atoms to the density of state at the bottom of the valence band. The diagonalized components of polarizability imply that there are possible micro-domains in the supercell while applying externally electric field to it. And the supercell presents a nearly isotropic polarizability macroscopically due to the strong interaction among the electric dipole moments existing in the different domains.

Piezo-electrochemical coupling of AgNbO3 piezoelectric nanomaterials Hot!

Hong Yuan-Ting, Ma Jiang-Ping, Wu Zheng, Ying Jing-Shi, You Hui-Lin, Jia Yan-Min
Acta Physica Sinica. 2018, 67 (10): 107702 doi: 10.7498/aps.67.20180287
Full Text: [PDF 2034 KB] Download:(37)
Show Abstract

In this work, the AgNbO3 piezoelectric nanomaterials are hydrothermally synthesized, and they have an average particle size of~1 μm, which is obtained from scanning electron microscopy pattern. The AgNbO3 nanomaterial possesses an orthorhombic crystal structure with an mm2 point group symmetry, indicated by the X-ray powder diffraction analysis result. The piezo-electrochemical coupling of AgNbO3 is characterized, and its physical mechanism is discussed. Under an external mechanical vibration, the surfaces of the piezoelectric AgNbO3 nanomaterials will generate a large number of positive and negative electric charges. Due to the existence of spontaneous polarization, these positive and negative electrical carriers are respectively distributed on the top surface and bottom surface of AgNbO3 and can further induce the generation of some strong oxidation middle active species such as hydroxyl radicals in solution on the basis of some special chemical redox reactions, realizing the piezo-electrochemical coupling. Therefore, we can consider the piezo-electrochemical coupling as the product of the piezoelectric effect and the electrochemical redox effect. Utilizing the strong piezo-electrochemical coupling, a practical application in mechano-catalysis is further developed to decompose dye solution under a driven vibration. After experiencing~60 min vibration with AgNbO3 nanomaterial as mechano-catalyst,~70% rhodamine B (~5 mg/L) is decomposed. Prior to the vibration, the rhodamine B solution with the addition of AgNbO3 catalyst is slowly stirred for 30 min to ensure the establishment of the physical adsorptiondesorption equilibrium between catalyst and dye. It is difficult to directly exert a mechanical stress on the micro/nanoparticles. Here, an ultrasonic source with a vibration frequency of~40 kHz is employed to exert a stress to compress and stretch the AgNbO3 particles through utilizing micro-bubble collapse forces during ultrasonic cavitations, which needs the AgNbO3 particle size to be roughly identical with the diameter (~μm) of micro-bubble. Our mechanocatalytic dye decomposition experiment is conducted at room-temperature and in a dark environment to avoid the influence of photocatalysis. The slight increase of temperature of the dye solution in the ultrasonic vibration process has no obvious influence on the dye decomposition efficiency, which has been confirmed from our experiment. Through a technology of fluorescence spectrum trapping, the intermediate active product in the piezo-electrochemical coupling process-the strongly oxidized hydroxyl radicals, is successfully observed. With the increase of vibration time, the number of hydroxyl radicals obviously increases, which proves that the piezo-electrochemical coupling plays a key role in our mechano-catalytic process. After using AgNbO3 catalyst in cyclic decomposition of rhodamine B 5 times, no obvious reduction in the piezo-electrochemical coupling performance occurs. The AgNbO3 nanomaterial possesses an efficient piezo-electrochemical coupling for mechano-catalysis, and it has the advantages of high decomposition efficiency and reusability, and potential applications in vibration decomposing dye.

Polarization-controlled dual-band broadband infrared absorber

Yang Peng, Han Tian-Cheng
Acta Physica Sinica. 2018, 67 (10): 107801 doi: 10.7498/aps.67.20172716
Full Text: [PDF 1063 KB] Download:(10)
Show Abstract
As an important branch of metamaterial-based devices, metamaterial absorber (MA) has aroused great interest and made great progress in the past several years. By manipulating the magnetic resonance and the electric resonance simultaneously, the effective impedance of MA will match the free space impedance, thus resulting in a perfect absorption of incident waves. Due to the advantages of thin thickness, flexible design and tunable property, MA has been extensively studied at various frequencies, e. g. microwave frequency, THz, infrared frequency, and optical frequency. Infrared MA, having important applications in infrared stealth, infrared detection, radiative cooling, and sensors, receives more and more attention, especially for those absorbers based on easy-fabricated one-dimensional grating structure. However, such a grating-based absorber is usually workable in narrow band and effective only for transverse magnetic (TM) wave.
In this paper, a dual-band broadband absorber is proposed based on the easy-fabricated grating structure. The basic unit of the proposed absorber consists of eight gradient subunits, each of which is composed of vertically cascaded two pairs of metal-dielectric bilayers. The as-designed absorber has perfect absorption for both TM and transverse electric (TE) waves. More importantly, the absorption band is different for different polarized wave, which provides more choices and greater flexibility for application. Full-wave simulation shows that the absorption of TM wave is above 90% from 1.68 μm to 2 μm, while the absorption of TE wave is very small (no more than 6%). The absorption of TE wave is above 90% from 3.8 to 3.9 μm, while the absorption of TM wave is very small (no more than 5%). In order to reveal the working principle of the proposed absorber, the electric-field distributions of the whole structure are calculated at different frequency, which demonstrates that the broadband absorption is achieved by exciting multiple resonant coupling. Furthermore, we investigate the performance of the proposed absorber in oblique incidence, and find that the designed absorber can exhibit a good absorption within a broad incident angle ranging from 0 to 60 degrees. It is worth noting that there is an absorption fracture band in the absorption spectrum of TM waves, which is because no resonance occurs in all subunits, resulting in almost no absorption.
In conclusion, we have proposed a dual-band broadband absorber that demonstrates independent absorption of the TM waves and the waves in different bands, which has potential applications in thermal detectors and thermal emitters. The proposed scheme can be extended to microwave, THz, and even visible light band.
INTERDISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY

Modes characteristics analysis of THz waveguides based on three graphene-coated dielectric nanowires

Wei Zhuang-Zhi, Xue Wen-Rui, Peng Yan-Ling, Cheng Xin, Li Chang-Yong
Acta Physica Sinica. 2018, 67 (10): 108101 doi: 10.7498/aps.67.20180036
Full Text: [PDF 714 KB] Download:(22)
Show Abstract
In this paper, the real parts of the effective refractive indexes and the propagating lengths of five low-order modes of the terahertz waveguides based on three graphene-coated dielectric nanowires are analyzed by using the multipole method. The formation of these five lowest order modes can be attributed to the five combinations between the two lowest order modes supported when three nanowires exist alone. Therefore they are named Mode 1, Mode 2, Mode 3, Mode 4, and Mode 5 in sequence. The results show that the mode characteristics of the waveguide can be effectively tuned by changing the operating frequency, the radius of the intermediate nanowire, the gap distance between the nanowires and the Fermi energy of graphene. As the operating frequency increases from 30 THz to 40 THz, the real part of each of the effective refractive indexes increases and the propagation length decreases, and the crossover phenomenon occurs in the process of change. In addition, the real parts of the effective refractive indexes and the propagation lengths of Modes 3 and 4 are basically the same. When the radius of the middle nanowire increases from 25 nm to 75 nm, the real parts of the effective refractive indexes of Modes 1 and 2 increase, and the propagation length of Mode 1 decreases and then increases. Besides the real parts of the effective refractive indexes and the propagation lengths of Modes 3 and 4 are basically not affected by the change of radius, and the values of these two modes are basically the same. For Mode 5, the real part of the effective refractive index and propagation length slowly increase. When the spacing between the nanowires increases from 10 nm to 50 nm, Modes 3 and 4 are basically unaffected by the change of spacing, and the values of these two modes are basically the same. The real parts of the effective refractive indexes of the other modes decrease and the propagation lengths increase and eventually stabilize, and the crossover phenomenon occurs in the process of change. As the Fermi energy of graphene increases from 0.4 eV to 1.2 eV, the real part of the effective refractive index decreases and the propagation length increases. The calculation shows that the result obtained by the multipole method is exactly the same as that obtained by the finite element method. To date, no one has analyzed the terahertz waveguides based on three graphene-coated dielectric nanowires. This work can provide a theoretical basis for the design, fabrication and application of terahertz waveguide based on graphene-coated dielectric nanowires. Such waveguides have potential applications in the field of mode-division multiplexing.

Development of quantum voltage noise source chip for precision measurement of Boltzmann constant

Wang Lan-Ruo, Zhong Yuan, Li Jin-Jin, Qu Ji-Feng, Zhong Qing, Cao Wen-Hui, Wang Xue-Shen, Zhou Zhi-Qiang, Fu Kai, Shi Yong
Acta Physica Sinica. 2018, 67 (10): 108501 doi: 10.7498/aps.67.20172643
Full Text: [PDF 3695 KB] Download:(9)
Show Abstract
The Johnson noise thermometer is used to precisely measure Boltzmann constant by comparing the thermal noise caused by charge movement and the quantized voltage reference noise synthesized by the quantum voltage noise source (QVNS). The QVNS signal is synthesized based on quantized voltage pulses produced by two channels of superconducting Josephson junction arrays, which are designed for cross-correlation electronics. The Nb/NbxSi1-x/Nb Josephson junction is used as a core device of QVNS chip in this work for its non-hysteretic current-voltage (I-V) characteristics and conveniently adjustable barrier parameters.
In this paper, we present the design consideration, fabrication process, and measurement results of the QVNS chip. The QVNS chip contains two Josephson junction arrays, each consists of four 6 μm×12 μm junctions and is embedded in a 50 Ω coplanar waveguide transmission line. The random noise in signals from the two driven channels is eliminated by cross-correlation, and then an accurate quantum noise is obtained. Test chips with different areas of Josephson junctions are also designed on the same mask, aiming at estimating the variation range of Ic. The typical fabrication process for voltage standard chips in our laboratory is used for preparing the QVNS chip.
The sample is measured at 4.2 K. The DC I-V curve shows that the critical current Ic is 6.1 mA. The I-V characteristics of the junctions under 5 GHz microwave radiation are measured. For a series array of four junctions, a 41.44 μV one-stage Shapiro step is observed. Calculation shows that the error between the measurement and theoretical value of 41.36 μV is about 1.9‰, which means that the QVNS chip performs well under microwave radiation and can be used for synthesizing the AC quantum voltage reference noise.
A single-frequency 100 kHz sinusoidal waveform is synthesized by the QVNS chip under pulse driven signal. A spectrum of the synthesized sinusoidal waveform shows a single peak, which means that the digital pulse signal is perfectly filtered by Josephson junction arrays and the synthesized signals possess quantum accuracy. The results indicate that our chip has good dynamic response and works well in synthesizing a single-frequency AC quantum voltage signal. This work can provide core devices for the noise thermometry system and support the precise measurement of Boltzmann constant as well as redefinition of Kelvin in future. As a next step, the design and package will be further improved, and the probe module will be optimized to reduce the measurement uncertainty.

Molecular dynamics simulations on DNA flexibility: a comparative study of Amber bsc1 and bsc0 force fields Hot!

Xiong Kai-Xin, Xi Kun, Bao Lei, Zhang Zhong-Liang, Tan Zhi-Jie
Acta Physica Sinica. 2018, 67 (10): 108701 doi: 10.7498/aps.67.20180326
Full Text: [PDF 763 KB] Download:(25)
Show Abstract

The structural flexibility of DNA plays a key role in many biological processes of DNA, such as protein-DNA interactions, DNA packaging in viruses and nucleosome positioning on genomic DNA. Some experimental techniques have been employed to investigate the structural flexibility of DNA with the combination of elastic models, but these experiments could only provide the macroscopic properties of DNA, and thus, it is still difficult to understand the corresponding microscopic mechanisms. Recently, all-atom molecular dynamics (MD) simulation has emerged as a useful tool to investigate not only the macroscopic properties of DNA, but also the microscopic description of the flexibility of DNA at an atomic level. The most important issue in all-atom MD simulations of DNA is to choose an appropriate force field for simulating DNA. Very recently, a new force field for DNA has been developed based on the last generation force field of Amber bsc0, which was named Amber bsc1. In this work, all-atom MD simulations are employed to study the flexibility of a 30-bp DNA with the force fields of Amber bsc1 and Amber bsc0 in a comparative way. Our aim of the research is to examine the improvement of the new development of force field (Amber bsc1) in the macroscopic and microscopic properties of DNA, in comparison with the corresponding experimental measurements. All the MD simulations are performed with Gromacs 4.6 and lasted with a simulation time of 600 ns. The MD trajectories are analyzed with Curves+ for the last 500 ns, since the system reaches equilibrium approximately after ~100 ns. Our results show that the new force field (Amber bsc1) can lead to the improvements in the macroscopic parameters of DNA flexibility, i.e., stretch modulus S and twist-stretch coupling D become closer to experimental measurements, while bending persistence lengths lp and torsional persistence lengths C from the two force fields (bsc1 and bsc0) are both in good agreement with experimental data. Our microscopic analyses show that the microscopic structure parameters of DNA from the MD simulation with the Amber bsc1 force field are closer to the experimental values than those with the Amber bsc0 force field, except for slide, and the obvious improvements are observed in some microscopic parameters such as twist and inclination. Our further analyses show that the improvements in macroscopic flexibility from the Amber bsc1 force field are tightly related to the microscopic parameters and their fluctuations. This study would be helpful in understanding the performances of Amber bsc1 and bsc0 force fields in the description of DNA flexibility at both macroscopic and microscopic level.

Objective assessment method of image quality based on visual perception of image content

Yao Jun-Cai, Liu Gui-Zhong
Acta Physica Sinica. 2018, 67 (10): 108702 doi: 10.7498/aps.67.20180168
Full Text: [PDF 1345 KB] Download:(42)
Show Abstract
Objective image quality assessment (IQA) plays a very important role in transmission, encoding, and quality of service (QoS) of the image and video data. However, the existing IQA methods often do not consider image content features and their visual perception, so there is a certain gap between the objective IQA sores and the subjective perception. To solve this problem, in the study, we propose an objective IQA method based on the visual perception of image content, which combines the complexity characteristics of image content, and the properties of masking, contrast sensitivity and luminance perception nonlinearity of human visual system (HVS). In the proposed method, the image is first transformed using a nonlinear model of luminance perception to obtain the intensity perception image. Then, the intensity information is summed using the contrast sensitivity values of HVS and the average contrast values of the local image as a weighting factor of the intensity. The summed data information is taken as the content of human perceiving image, and an image perception model is constructed. Finally, the reference images and distorted images are perceived by simulating the HVS with this model. Moreover, the difference in intensity between two perceived images is calculated. Based on the intensity difference and peak signal-to -noise ratio model, an objective IQA model is constructed. Further, the simulation with 47 reference images and 1549 test images in the LIVE, TID2008, and CSIQ databases is conducted. Moreover, the experimental results are compared with those of four typical objective IQA models, namely SSIM, VSNR, FSIM, and PSNRHVS. In addition, we explore the factors that affect the IQA accuracy and a way to improve assessment accuracy by combining HVS characteristics, through analyzing the correlation between IQA results of the proposed model and the subjective mean opinion scores (MOSs) provided in the three image databases from the following two aspects. Namely, (1) all reference images in three image databases are distorted by multiple types, and the distorted images of each reference image are taken as a test sequence. Then, the proposed model is used to evaluate each test sequence to obtain the IQA scores. By analyzing the correlation between the IQA scores of each test sequence and the subjective MOSs and comparing them with the assessment results of SSIM, we explore the influence of the image content complexity on the objective IQA accuracy. (2) The test images which are distorted by each type and many distortion degrees are used as another sequence, and they are evaluated by the proposed IQA model. By analyzing the correlation between the subjective MOSs and the IQA results of each test sequence, and comparing them with assessment results of SSIM, we discuss the influence of image distortion mode on the IQA accuracy. The experimental results show that the coefficient values of Pearson linear correlation and Spearman rank order correlation between the objective IQA scores obtained by the proposed method and the subjective MOSs have been averagely improved by 9.5402% and 3.2852%, respectively, in comparison with IQA results from the SSIM method. Also, they are enhanced more significantly than those fom the PSNRHVS and VSNR methods. In summary, it is shown that the proposed IQA method is an effective and feasible method of objectively assessing the image quality; moreover, it is shown that in the objective assessment of image quality it is very helpful to improve the consistency of subjective and objective assessment of image quality by considering the content perception and complexity analysis of the images.

Concentrating characteristics of Fresnel lens with prism secondary concentrator and optimization of high concentrating photovoltaic module with triple-junction cell

Guo Yin, Shu Bi-Fen, Wang Jing, Yang Qing-Chuan, Jiang Jing-Xiang, Huang Yan, Zhou Zheng-Long
Acta Physica Sinica. 2018, 67 (10): 108801 doi: 10.7498/aps.67.20172778
Full Text: [PDF 2737 KB] Download:(9)
Show Abstract
At present, Fresnel lens is commonly used as a concentrator in high concentrating photovoltaic (HCPV) module, and the triple-junction cell is currently one of the most common multi-junction cells used in it. The triple-junction cell is composed of three p-n junctions in series. Each sub-cell in it absorbs different-wavelength light. The solar cell efficiency of Ⅲ-V multi-junction high concentrating photovoltaic increases up to 46%, which the corresponding module efficiency is quite different from. The output power of the solar cell is related to not only the illumination energy, but also the spectral distribution and the uniformity of the illumination. The loss caused by the non-ideal concentration of the concentrator in the module is as high as 20%. After sunlight enters the lens, the direction of transmission of a monochromatic light is different, because a lens has different refractive index for different-frequency light. So the light disperses when leaving the lens, and thus the colors are arranged in a certain order to form a spectrum. Owing to the dispersion and the differences in refractive index among different spectral bands, the illumination distributions of the three spectral bands are different and non-uniform on the focal plane of lens. The divergence of light will obviously weaken the non-uniformity of the illumination on the solar cell surface. So the divergence angle of the light source has a greater influence on the cell performance because the non-uniformity of illumination has a negative influence on the performance of the cell.
In this paper, according to the establishment of optical model and three-dimensional cell circuit network model under non-uniform illumination, taking Ⅲ-V triple-junction cells for example, we study the concentrating characteristics and photovoltaic characteristics of HCPV module with Fresnel lens concentrator and prism secondary concentrator. The results show that due to the non-parallel incident light and dispersion of the Fresnel lens, the concentrating spots of short-wave light, medium-wave light and long-wave light are divergent and their illuminations are non-uniform, resulting in the spectral response mismatch loss of the three sub-cells in the triple-junction cell, and the photovoltaic performance of the HCPV module also declines. The results show that the secondary optics element is obviously effective in reducing the non-uniformity of the illumination and the temperature which the Fresnel lens creates. However, each waveband of light has a different spot size at the same position, similar to the Fresnel lens without the secondary optics element. So the varieties of cell performance at different positions are similar too. And, by optimizing the focusing characteristics of the three wave bands along the optical axis, the power output of the HCPV module can increase more than 10%. The simulation results are verified experimentally.
GENERAL

Compressed sensing based fast method of solving the electromagnetic scattering problems for threedimensional conductor targets

Ding Ya-Hui, Sun Yu-Fa, Zhu Jin-Yu
Acta Physica Sinica. 2018, 67 (10): 100201 doi: 10.7498/aps.67.20172543
Full Text: [PDF 488 KB] Download:(72)
Show Abstract

The method of moments is one of the most commonly used algorithms for analyzing the electromagnetic scattering problems of conductor targets. However, it is difficult to solve the matrix equation when analyzing the electromagnetic scattering problem of the electric large target. In recent years, the theory of the compressed sensing was introduced into the method of moments to improve the computation efficiency. The random selected impedance matrix is used as a measurement matrix, and the excitation voltage is used as a measurement value when using compressed sensing theory. The recovery algorithm is used to solve the induced current of target. The method can avoid the inverse problem of matrix equation and improve the computational efficiency of the method of moments, but it can be applied only to 2-dimensional or 2.5-dimensional target. The application of compressed sensing needs to know the sparse basis of the current in advance, but the induced current of three-dimensional target which is expressed by an Rao-Wilton-Glisson basis function is not sparse on the commonly used sparse basis, such as discrete cosine transform basis and discrete wavelet basis. To solve this problem, a method of combining compressed sensing with characteristic basis functions is proposed to analyze the electromagnetic scattering problem of three-dimensional conductor target in this paper. The characteristic basis function method is an improved method of moments. The target is divided into several subdomains, the main characteristic basis functions are comprised of current bases arising from the self-interactions within the subdomain, and the secondary characteristic basis functions are obtained from the mutual coupling effects of the rest of the subdomains. Then a reduction matrix is constructed to reduce the order of matrix equation, and the current can be expressed by the characteristic basis function and its weighting coefficient. In the method presented in this paper, the weighting coefficient is considered as a sparse vector to be solved when the characteristic basis function is used as sparse basis. The number of weighting coefficients is less than the number of unknown ones, so it can be obtained from the compressed sensing recovery algorithm. At the same time, the generalized orthogonal matching pursuit algorithm is used as the recovery algorithm to speed up the recovery process. Finally, the proposed method is used to calculate the radar cross sections of a PEC sphere, nine discrete PEC targets and a simple missile model. The numerical results validate the accuracy and efficiency of the method.

Novel dual channel polarization interference imaging system

Wang Tian, Niu Ming-Sheng, Bu Miao-Miao, Han Pei-Gao, Hao Dian-Zhong, Yang Jing-Shun, Song Lian-Ke
Acta Physica Sinica. 2018, 67 (10): 100701 doi: 10.7498/aps.67.20172691
Full Text: [PDF 2367 KB] Download:(38)
Show Abstract

The interference images with fixed spectral resolution can be obtained by using the existing static polarization-difference imaging system because the optical path of the system cannot be changed flexibly. However, for different detection targets, the spectral resolution of the system determined by the optical path difference must be appropriate. To satisfy a variety of application requirements, a novel dual channel polarization-difference interference imaging system (DPDⅡS), based on the lateral shear of the wide-field-of-view Savart polariscope (WSP) and the modulated Savart polariscope (MSP), is presented. The two-dimensional space images of a target and orthogonal interference images can be obtained by adjusting the MSP under different lateral displacements simultaneously. In addition, the remarkable characteristics of the system avoid spilling over rays and optimizing the system optical path effectively. In this paper, by using the Jones matrix, the system structure is demonstrated and the theoretical principle of DPDⅡS is analyzed in detail. The amplitudes of the four beams from the MSP and the interference intensity expressions of the coherent light are derived. Then the splitting characteristics of the Savart polariscope (SP) and WSP are presented. It is concluded that the WSP has better shear ability than SP and the WSP can optimize the optical path effectively compared with Wollaston prism in the DPDⅡS. The change ranges of the optical path difference and lateral displacement produced by the MSP for structure angles α=π/3, π/4, π/6 are analyzed in detail. The reconstructed orthogonal interferograms and the experimental interferograms under 632.8 nm monochromatic light for dMSP=1.00, 1.10, 1.20, 1.30 mm are obtained. A comparison between the experimental interference images and the simulated images proves that the interference fringes with different resolutions can be obtained simultaneously by adjusting the MSP. Meanwhile, the light intensities of the double optical paths are approximately equal and the same optical path difference is generated for the dual channel with the movement of MSP. The experimental results are consistent with the theoretical analyses. The spatial images of parallel and vertical components are detected under 632.8 nm polychromatic light. Then the total intensity image and the polarization-difference image are obtained through data processing. The conclusion that the polarization difference intensity image has a high resolution compared with the polarization intensity image is presented. The study has reference significance and practical value for the dual channel polarization interference imaging system.

ATOMIC AND MOLECULAR PHYSICS

Multiconfiguration time-dependent Hartree-Fock treatment of electron correlation in strong-field ionization of H2 molecules

Zhang Bin, Zhao Jian, Zhao Zeng-Xiu
Acta Physica Sinica. 2018, 67 (10): 103301 doi: 10.7498/aps.67.20172701
Full Text: [PDF 652 KB] Download:(30)
Show Abstract
Electron correlation plays an important role in the multielectron interactions of many physical and chemical processes.The investigation of correlation effects in the non-perturbative electronic dynamics (e.g.non-sequential double ionization) when atoms and molecules are subjected to strong laser fields requires non-perturbative theoretical treatments. The direct numerical integration of the time-dependent Schrödinger equation successfully explains many experimental results,but it is computationally prohibitive for systems with more than two electrons.There is clearly a need for a theory which can treat correlation dynamics self-consistently in strong time-dependent electric fields.In this paper we develop a three-dimensional multiconfiguration time-dependent Hartree-Fock method,which can be applied to the non-perturbative electronic dynamics for diatomic molecules,and it can also investigate the effect of electron correlation in strong-field ionization of H2 molecules.This method adopts the prolate spheroidal coordinates (which can treat the two-center Coulomb potential accurately) and the finite-element method together with discrete-variable representation (which lowers the calculation burden from two-electron integrations).For the temporal propagation,we use the efficient short iterative Lanczos algorithm for the equation which governs the configuration expansion coefficients,while an eight-order Runge-Kutta (RK) method and an Bulirsch-Stoer (BS) extrapolation method,both with adaptive precision controls,are implemented to solve the nonlinear orbital equation.While both methods yield correct results,the BS method displays a better stability in the realtime propagation,while the RK method demands less computation.The alignment-dependent ionization probabilities of H2 molecules in intense extreme ultraviolet pulses are calculated.Comparisons between multi-configuration and single-configuration results show that electron correlation has little effect on the single ionization process,but it plays an important role in double ionization,leading to the decrease in the ionization probability.The double ionization probability from the single-configuration space 1σ is about three times larger that from 4σ1π.The ionization probability increases monotonically when the alignment angle increases from 0° to 90°, yielding a ratio of 2.6(1.5) between 90° and 0° for the double (single) ionization process.This method presents the basis for the future study of electron correlation in strong-field processes.

H/D + Li2 → LiH/LiD + Li reactions studied by quantum time-dependent wave packet approach

Li Wen-Tao, Yu Wen-Tao, Yao Ming-Hai
Acta Physica Sinica. 2018, 67 (10): 103401 doi: 10.7498/aps.67.20180324
Full Text: [PDF 684 KB] Download:(24)
Show Abstract

The isotopic effect is a significant way to further understand the reaction mechanism without greatly changing the system. However, the isotopic effect of the H + Li2 reaction has received little attention in previous theoretical studies. Furthermore, as a deep potential well exists on the reaction path, obtaining convergent result is very time-consuming. So some approximate methods were used in previous theoretical calculations. However the Coriolis coupling effect plays an important role in the reaction, and thus whether these approximate methods are reasonable needs further testing. Based on the potential energy surface (PES) reported by Song et al., the dynamical calculations of H/D + Li2→ LiH/LiD + Li reactions are carried out by time dependent quantum wave packet method with second order split operator in a collision energy range from 0 to 0.4 eV. In order to obtain the convergent results, lots of convergence tests are carried out and because the Coriolis coupling effect plays an important role in the reaction, all the number of projections of total angular momentum J are included in the present calculation. The dynamical properties such as reaction probability, integral cross section, differential cross section are calculated and compared with previous theoretical values. Large discrepancies are found between present results and the values obtained from Gao et al. especially at high collision energies. Owing to the fact that the same PES is applied to the calculation and Gao's results of total angular momentum J=0 accord well with the present values, we suppose that the parameters used in the calculation have little influence on the final results and the main discrepancies are attributed to the number of projections of total angular momentum which are cut off in Gao et al.'s calculation. In order to verify our speculation, the numbers of projections of total angular momentum which are 1, 5, 10, 15, 20, and 25, are considered in the calculation, respectively. The results indicate that the main discrepancy between present values and the results obtained from Gao et al. can be attributed to the number of projections of total angular momentum used in Gao et al.'s calculation that is not convergent, and that the present values are more accurate than previous theoretical studies for all the numbers of projections of total angular momentum which are included in the calculation. Furthermore, when the H atom is substituted by the heavy isotope D atom, the reaction probability and integral cross section become large. However, it does not generate large effect on the reaction mechanism. The forward and backward symmetry differential cross section signals indicate that the complex forming reaction mechanism dominates the reaction.

ELECTROMAGNETISM, OPTICS, ACOUSTICS, HEAT TRANSFER, CLASSICAL MECHANICS, AND FLUID DYNAMICS

Focusing properties of Lucas sieves

Xie Jing, Zhang Jun-Yong, Yue Yang, Zhang Yan-Li
Acta Physica Sinica. 2018, 67 (10): 104201 doi: 10.7498/aps.67.20172260
Full Text: [PDF 883 KB] Download:(26)
Show Abstract

A kind of optical diffractive element named photon sieve, which is essentially Fresnel zone plate in which the transmissive rings are replaced with a large number of randomly distributed isolated pinholes, can be used to focus X-ray and extreme ultraviolet lithography spectrum into spots with sizes smaller than the diameter of the smallest circular pinhole. However, both the traditional photon sieves and Fibonacci sieves have no more than two axial foci. In order to break this limitation, the Lucas sequence is introduced into the design of photon sieves, and thus producing four axial foci. With respect to the previous Fibonacci sequence, Lucas sequence has the same recursion relation as well as the same eigenvalue of golden mean γ=(1 + √5)/2. The only difference between them is the first two initial seeds. Based on Fresnel-Kichhoff diffraction theory, the simulation results show that there exist four focal spots with approximately equal intensity along the optical axis on condition that the hole diameters are set to be 1.16 times the underlying Fresnel zone width. Then in order to verify the validity of our proposed model, a Lucas sieve of diameter 12.11 mm and referred focal length 180 mm is fabricated by photolithography and its focusing properties are precisely measured by the in-line phase-shifting digital holography. In experiment, a quarter wave plate is used to realize two-step phase-shift interferences, and obtain the quad-focal length by auto-focusing algorithm in holography. Meanwhile, the quad-focal spots can also be calculated through the diffraction propagation of reconstructed object wave. Compared with the theoretical values, the measurement results indicate that the maximum deviation of quad-focal lengths is less than 0.9%, and the relative errors of the full width at half maximum of four Airy spots are all less than 5%. The experimental results agree well with the theoretical analysis results. Owing to the advantages of small volume, little weight and easy processing, Lucas sieve has great potential in X-ray microscopy, array imaging for living biological cell and especially in the next generation of synchrotron light sources.

Transmission matrix optimization based on singular value decomposition in strong scattering process

Zhang Xi-Cheng, Fang Long-Jie, Pang Lin
Acta Physica Sinica. 2018, 67 (10): 104202 doi: 10.7498/aps.67.20172688
Full Text: [PDF 3693 KB] Download:(26)
Show Abstract

In the last decade, the scattering medium has been gradually attacking attention from researchers. Among the proposed approaches, the transmission matrix (TM) is considered as an effect way to describe the scattering properties which relate to input optical and output optical fields. However, the acquired transmission matrix and its eigenvalues strongly depend on the experimental conditions, such as the numbers of input channels (limited numerical aperture and illumination area, or the pixel number of the spatial light modulator) and output channels. In other words, the actual transmission matrix of the scattering medium is the acquired transmission matrix with infinite numbers of the input and output channels. We propose an approach to obtaining the actual matrix by evaluating its eigenvalues. First, the matrix is expressed by the singular value decomposition to obtain its inverse matrix. Then first level optimization is to dispose of some extreme singular values to remove the ill-conditioned problem of the matrix, and then, as a second level optimization, the genetic algorithm is to remove the eigenvalues which have the negative contributions to the intensity of the selected focal point. Our experiments show that the gray value of the intensity and the signal-to-noise ratio (SNR) of the focal point after employing the phase pattern are 129 and 7.54, respectively. After the first level optimization, the gray value of the intensity and the SNR rise to 172 and 9.73, respectively. Then, they reach to 192 and 10.29, respectively, after adopting the genetic algorithm. After the second level optimizations, the intensity at the focal point increases 48.8% compared with the case with just the optimized phase pattern from the acquired TM, and the SNR increases by nearly 36.5%. The reason behind the increase of the intensity after the optimizations, we believe, is that the transmission matrix of the scattering medium reaches its actual matrix in certain conditions. The proposed approach opens the way to obtaining the actual transmission matrix by mathematic optimizations without increasing the experimental levels.

Entanglement characteristics of output optical fields in double-cavity optomechanics

Zhang Xiu-Long, Bao Qian-Qian, Yang Ming-Zhu, Tian Xue-Song
Acta Physica Sinica. 2018, 67 (10): 104203 doi: 10.7498/aps.67.20172467
Full Text: [PDF 390 KB] Download:(16)
Show Abstract

Radiation pressure in an optomechanical system can be used to generate various quantum entanglements between the subsystems. Recently, one paid more attention to the study of quantum entanglement in an optomechanical system. Here in this work, we study the properties of output entanglement between two filtered output optical fields by the logarithmic negativity method in a double-cavity optomechanical system. Our calculations show that the decay rate of the mechanical resonator, the bandwidth of filter function, and non-equal-coupling will evidently affect the value of the output entanglement. In particular, under the parameters of equal-coupling and zero filter bandwidth, the output entanglement in the vicinity of resonant frequency (ω=0 in the rotating frame) will decease with mechanical decay rate increasing. But under the parameters of equal-coupling and non-zero filter bandwidth, the output entanglement will be suppressed if the center frequency of output field is in the vicinity of the resonant frequency. However, the output entanglement can be enhanced if we adopt a non-equal-coupling to counteract the suppression effect of the filter bandwidth. Furthermore, we find that there are three peaks in the whole center frequency domain of the output field if we adopt strong non-equal-coupling. This is because the normal mode of Hamiltonian Hint will split into three normal modes in this case. Our results can also be used in other parametrically coupled three-mode bosonic systems and may be applied to realizing the state transfer process and quantum teleportation in an optomechanical system.

Supperssion of higher order modes in gain-guided index-antiguided planar waveguide laser

Kang Da, Luo Bin, Yan Lian-Shan, Pan Wei, Zou Xi-Hua
Acta Physica Sinica. 2018, 67 (10): 104204 doi: 10.7498/aps.67.20180138
Full Text: [PDF 887 KB] Download:(6)
Show Abstract

In order to suppress the higher order modes and improve beam quality in high power waveguide laser, based on gainguided index-antiguided theory, a new symmetric layered waveguide structure is designed, and an interval layer is proposed to be sandwiched between waveguide layer and cladding layer in traditional symmetric GG-IAG waveguide structure. As a result, while reducing the leakage loss of fundamental mode, the threshold gain coefficient differences between fundamental mode and higher order modes will be further increased. When the gain in waveguide layer is between threshold gain coefficient of fundamental mode and that of higher order mode, the fundamental mode will have a greater advantage in mode competition than others, so higher order modes can be suppressed and the laser can obtain a single mode output. In the meantime, the guided-mode principle of this waveguide structure is explained with the theory of wave optics in this paper, the eigen equation of each mode is derived from the wave equation, and the field distributions of fundamental mode and higher order mode are also given. Additionally, in this paper we give the solution process of the threshold gain coefficient of each mode in this waveguide structure. The mode leakage losses of fundamental mode and higher order mode, after adding the interval layer, are numerically calculated, and the parameter optimization process of the interval layer is also given in this paper. In addition, the field distributions of fundamental mode and higher order mode are numerically simulated. The calculation results show that comparing with the traditional symmetric GG-IAG planar waveguide, after adding the interval layer, the loss of fundamental mode can be greatly reduced, while ensuring that the leakage loss of higher order mode reaches a maximum value by reasonably controlling the parameters of interval layer. In this way, we can suppress higher order modes and improve laser efficiency. This paper provides a new idea for improving the beam quality of high power waveguide laser with a large mode area.

Optimized design and epitaxy growth of high speed 850 nm vertical-cavity surface-emitting lasers

Zhou Guang-Zheng, Yao Shun, Yu Hong-Yan, Lü Zhao-Chen, Wang Qing, Zhou Tian-Bao, Li Ying, Lan Tian, Xia Yu, Lang Lu-Guang, Cheng Li-Wen, Dong Guo-Liang, Kang Lian-Hong, Wang Zhi-Yong
Acta Physica Sinica. 2018, 67 (10): 104205 doi: 10.7498/aps.67.20172550
Full Text: [PDF 1670 KB] Download:(58)
Show Abstract

Using transfer matrix method and TFcalc thin film design software,the reflectance spectrum of distributed Bragg reflector (DBR) and vertical cavity surface emitting laser (VCSEL) are simulated.The reflectance spectra from the cavity and surface are compared with each other,thus providing the basis for white light source (WLS) optical reflectance spectrum of the VCSEL epitaxial wafer.When using WLS to characterize VCSEL wafer,it is necessary to combine the simulation results and the shape of optical reflectance spectrum to speculate the reflectance seen from the cavity.The influences of different cap layers on the reflectance of DBRs are discussed theoretically and experimentally.With a 1/4λ GaAs cap layer,the reflectance reaches up to 97.8% seen from the cavity.This design can make the wavelength of the VCSEL etalon picked easily because of avoiding the influence of test noise.
The active region has higher heat accumulation due to the small area and poor thermal conductivity.The characteristics of the gain spectrum of InGaAs/AlGaAs strained quantum well (QW) under different temperatures and the temperature distribution in VCSEL are simulated by Crosslight software.The gain-to-cavity wavelength detuning is used to improve the slope efficiency and the temperature stability.The temperature in active region ranges from 360 K to 370 K.The gain peak wavelength and the Fabry-Perot cavity wavelength are designed in the ranges of 829-832 nm and 845-847 nm,respectively.Epitaxial wafer with top-emitting VCSEL structure grown by metal-organic chemical vapor deposition is characterized.The room temperature photoluminescence peak is at 827.5 nm and the etalon cavity wavelength measured by optical reflectance is 847.7 nm,which are consistent with designed values.
The oxide restricted VCSELs with 7.5 μm oxide aperture are fabricated.The image of the infrared light source CCD shows that the oxide aperture is circular.A passivation layer of 120 nm SiO2 is finally deposited to insulate water vapor.The threshold current is 0.8 mA,and the maximum output power reaches up to 9 mW at 13.5 mA.The optical spectrum at 6.0 mA reveals multiple transverse modes.The center wavelength is 852.3 nm and the root mean square (RMS) spectrum width is 0.6 nm,meeting the high-speed Datacom standards.Shannon theory indicates that the maximum data rate is not only proportional to bandwidth but also related to signal-to-noise ratio (SNR).It is effective to reduce relative intensity noise and enhance the SNR by increasing output power.From the eye diagram of 25 Gbit/s on-off key VCSEL,it is demonstrated that fall time is 38.66 ps,rise time is 41.54 ps,SNR is 5.6,and jitter RMS is 1.57 ps.Clear eye opening is observed from eye diagram of 25 GBaud/s PAM-4 VCSEL,which indicates the qualified 50 Gbit/s high speed performance.

Measurement of molecular absorption spectrum with a laser locked on a high-finesse cavity

Kang Peng, Sun Yu, Wang Jin, Liu An-Wen, Hu Shui-Ming
Acta Physica Sinica. 2018, 67 (10): 104206 doi: 10.7498/aps.67.20172532
Full Text: [PDF 485 KB] Download:(17)
Show Abstract

High-resolution and high-sensitivity molecular spectroscopy is widely used in fundamental molecular physics, atmospheric studies, remote sensing, industrial process monitoring, and medical diagnostics. Accurate determination of the parameters of molecule absorption lines, such as line positions, line strengths, line widths and profiles, is essential to support these studies and applications. For example, in order to retrieve the column density of carbon dioxide with a precision of one part per million (ppm), we need laboratory data of line positions with a uncertainty lower than 0.3 MHz and line intensities with a relative accuracy better than 0.5%.
Here we present precision spectroscopy of molecules using a laser locked with a high-finesse cavity. The cavity made of invar is thermo-stabilized to reduce the drifts of its length and the cavity mode frequencies. The frequency of the probe laser is locked on a longitudinal mode of the cavity by using the Pound-Drever-Hall method. Another beam from the probe laser, which is frequency shifted and on resonance with a nearby longitudinal mode of the cavity, is used for cavity ring-down spectrum (CRDS) measurement. The CRDS absorption spectrum is recorded by stepping the modulation frequency of a fiber electro-optic modulator in increment of the mode spacing of the cavity. Note that the cavity mode frequencies are shifted due to the dispersion introduced by the absorption lines. Prior to the CRDS measurements, the transmittance spectra of the cavity modes are recorded by scanning the probe laser frequencies over the resonance, which allows the determination of the cavity mode frequencies with an accuracy at a Hz level. Therefore, a dispersion spectrum is also obtained using the same setup by measuring the frequency shifts of cavity modes of the samples with and without absorption. The absolute frequency of the probe laser is determined by an optical frequency comb referring to a GPS-disciplined rubidium clock. The long term drift of beat frequency between the optical frequency comb and the probe laser is measured to be about 1.8 MHz per hour, which is consistent with the thermal expansion of the cavity under a temperature drift of 50 mK.
The performance of the spectrometer is demonstrated by measuring the Doppler-broadened spectra of CO2 around 6470.42 cm-1. Precise spectroscopic parameters are derived from both the absorption and dispersion spectra recorded by the same spectrometer. The line position is determined with an accuracy of 0.18 MHz, which is over one order of magnitude better than those given in previous studies and spectral databases.

Frequency locking of fiber laser to 1530.58 nm NH3 sub-Doppler saturation spectrum based on noise-immune cavity-enhanced optical heterodyne molecular spectroscopy technique

Jia Meng-Yuan, Zhao Gang, Zhou Yue-Ting, Liu Jian-Xin, Guo Song-Jie, Wu Yong-Qian, Ma Wei-Guang, Zhang Lei, Dong Lei, Yin Wang-Bao, Xiao Lian-Tuan, Jia Suo-Tang
Acta Physica Sinica. 2018, 67 (10): 104207 doi: 10.7498/aps.67.20172541
Full Text: [PDF 2504 KB] Download:(13)
Show Abstract
Noise-immune cavity-enhanced optical heterodyne molecular spectroscopy (NICE-OHMS) is a powerful tool for trace gas detection, which is based on the combination of frequency modulation spectroscopy (FMS) for reduction of 1/f noise, especially residual intensity noise, and cavity enhanced absorption spectroscopy (CEAS) for prolonging the interaction length between the laser and the targeted gas. Because of the locking of modulation frequency in FMS to the free spectral range (FSR) of the cavity, NICE-OHMS is immune to the frequency-to-amplitude noise, which is a main limitation to CEAS. Moreover, due to the building of high power inside the cavity, NICE-OHMS can easily saturate the molecular absorption thus obtain sub-Doppler spectroscopy, which possess a high resolution and odd symmetry, and thus can act as a frequency discriminator for the locking of the laser frequency to the transition center. In this paper, a fiber laser based NICE-OHMS system is established and the laser frequency is locked to the sub-Doppler absorption line of NH3 by sub-Doppler NICE-OHMS. To avoid the complex design of high-Q-factor bandpass filter at radio frequency, the frequency νpdh, used for Pound-Drever-Hall (PDH) locking, is generated by the beat frequencies νfsr and νdvb, which are used for NICE-OHMS signal and DeVoe-Brewer (DVB) locking, respectively. The performances of PDH and DVB locking are analysed by the frequency distribution deduced from the error signals, which result in frequency deviations of 4.3 kHz and 0.38 kHz, respectively. Then, the CEAS signal and NICE-OHMS signal in the dispersive phase for the measurement of NH3 at 1.53 μm under 70 mTorr are obtained, which show signal-to-noise ratios of 3.3 dB and 45.5 dB, respectively. Due to the high power built in the cavity, the sub-Doppler structure in the NICE-OHMS signal is obtained in the center of the absorption tansition with a satruation degree of 0.22, which is evaluated by the amplitude ratio between sub-Doppler and Doppler-broadened signals. The linewidth (full width at half maximum) of the sub-Doppler signal of 2.05 MHz is obtained, which is calibrated by the time interval between carrier and sideband. The free-running drift of the laser frequency is estimated by the NICE-OHMS signal and results in 50 MHz over 3 h. While, with locking, the relative deviation of the laser frequency is reduced to 16.3 kHz. In order to evaluate the long term stability of the system, the frequency deviation over 3 h is measured. The Allen deviation analysis shows that the white noise is the main noise of the system in the integration time shorter than 10 s. And the frequency stability can reach to 1.6×10-12 in an integration time of 136 s.

Optimized linear wavenumber spectrometer based spectral-domain optical coherence tomography system

Wu Tong, Sun Shuai-Shuai, Wang Xu-Hui, Wang Ji-Ming, He Chong-Jun, Gu Xiao-Rong, Liu You-Wen
Acta Physica Sinica. 2018, 67 (10): 104208 doi: 10.7498/aps.67.20172606
Full Text: [PDF 1622 KB] Download:(9)
Show Abstract
In spectraldomain optical coherence tomography the sample is illuminated by a broadband light source, and the spectrum of the interference light between the light returned from the sample and a reference mirror is detected by a grating spectrometer. Conventionally, the grating spectrometer is comprised of a diffraction grating, a focusing lens, and a line-scan camera. According to the grating equation the diffraction angle from the grating is approximately linearly related to the optical wavelength. Thus the distribution function of the light spectrum at the line-scan camera is nonlinearly dependent on wavenumber. For the high-quality image reconstruction, the numerical resampling of the spectral interference data from wavelength-space to wavenumber-space is commonly required prior to the Fourier Transformation. The nonlinear detection of the spectral interferograms in wavenumber space also degrades the depth-dependent signal sensitivity in conventional linear-wavelength spectrometer based spectraldomain optical coherence tomography. Recently reported spectraldomain optical coherence tomography based on a linearwavenumber spectrometer does not need the resampling or interpolating of the nonlinearwavenumber interference spectral data, which greatly reduces the cost of computation and improves the imaging sensitivity. Various methods based on the different evaluation protocols for optimizing the design of the linear-wavenumber spectrometer have been reported. Here we report an effective optimization method for linear-wavenumber spectrometer used in a high-resolution spectral domain optical coherence tomography system. We take the reciprocal of the fullwidthhalfmaximum of the simulated point spread function as an evaluating criterion to optimize the structure parameters of the linearwavenumber spectrometer, including the refractive index and the vertex angle of the dispersive prism and the rotation angle between the diffraction grating and the dispersive prism. According to the optimization, an F2 equilateral dispersive prism is used to construct the optimized linearwavenumber spectrometer with a rotation angle of 21.8°. We construct an optimized linearwavenumber spectrometer and implement the spectrometer in a developed spectraldomain optical coherence tomography system as a detection unit. We evaluate the performances of the linear-wavenumber spectrometer both theoretically and experimentally. The experimentally measured axial resolution of the spectraldomain optical coherence tomography system based on the linear-wavenumber spectrometer is 8.52 μm, and the sensitivity is measured to be 91 dB with -6 dB sensitivity roll-off within a depth range of 1.2 mm. The experimentally measured sensitivity roll-off curve accords well with the theoretical sensitivity roll-off curve. Utilizing the general parallel computing capability of a GPU card, the highquality spectraldomain optical coherence tomography images of the human finger skin can be reconstructed in real time without any resampling or interpolating process.

Principle analysis of snapshot Mueller matrix imaging polarimeter using birefringent crystal

Cao Qi-Zhi, Yuan Chang-An, Hu Bao-Qing, Ren Wen-Yi, Zhao Yin-Jun, Zhang Jing, Li Jian-Ying, Deng Ting, Mingwu Jin
Acta Physica Sinica. 2018, 67 (10): 104209 doi: 10.7498/aps.67.20172604
Full Text: [PDF 4356 KB] Download:(20)
Show Abstract
Conventional Muller matrix imaging polarimeter (MMIP) with several rotating elements suffers mechanical complexity, vibration noise, heat generation, and other unwanted problems. To overcome those shortcomings, we present a snapshot Muller matrix imaging polarimeter (SMMIP) using a birefringent crystal with high extinction ratio. The snapshot imaging polarimeter allows a single image to be used to measure the polarization of a scene without electronic control units or moving mechanical components. This new polarimeter combines the technique of Muller matrix spectropolarimetry with a snapshot imaging polarimeter through using modified Savart polariscope (MSP-SMMIP). It contains both a generator and an analyzer module. Spatial polarization fringes are localized on a sample by incorporating modified Savart polariscope into a polarization generator module. These fringes modulate the Mueller matrix components of the sample, which are subsequently isolated with modified Savart polariscope in an analyzer module, and the analyzer and the imaging lens combine with 16 beams to create interference, resulting in spatial modulation on the two-dimensional CCD camera. Expressions for interference intensities, optical system analysis, theory of calibration and method of reconstruction are presented. Finally, the numerical simulation is used to demonstrate theoretical analysis and the feasibility of MSP-SMMIP. The layout is very easy to calibrate and the reference target is only a linear polarizer at 22.5°. Moreover, the remarkable advantages of the proposed instrument, compared with conventional Muller matrix imaging polarimeter, are that it is also simple, compact, snapshotted, and static (no moving parts). Therefore we believe that the proposed snapshot imaging polarimeter will be very useful in many applications, such as biomedical imaging and remote sensing.

Characteristics of granular sheet of dense granular jet oblique impact

Wang Yue, Li Wei-Feng, Shi Zhe-Hang, Liu Hai-Feng, Wang Fu-Chen
Acta Physica Sinica. 2018, 67 (10): 104501 doi: 10.7498/aps.67.20172092
Full Text: [PDF 2684 KB] Download:(17)
Show Abstract
Dense granular impinging jets widely exist in natural flow phenomena and industrial processes, such as the rapid heating, cooling or drying, and gasification. It is important to investigate the factors influencing the flow patterns of dense granular impinging jets and reveal the evolution rules of the flow patterns. The dynamic behaviors of the dense granular impinging jets are experimentally studied by a high-speed camera and image processing software of Image J. The effects of the particle diameter, the granular jet velocity (u0) and the solid content of the granular jet (xp) on flow pattern of the granular impinging jet are investigated. Two flow regimes of the dense granular impinging jets, i.e., the liquid-like granular film and the scattering pattern, are identified. The results show that with the increase of the particle diameter and the granular jet velocity, both the solid content of the granular jet and the inter-particle collision frequency decrease, which results in the transition of granular sheet to scattering pattern. With the increase of granular jet velocity, the opening angle of the granular sheet from the side view increases, while the opening angle from the front view increases first and sharply decreases later. The results also show that with the increase of the granular jet velocity, the liquid-like granular film becomes unstable and a non-axisymmetric oscillation appears. And the amplitude and frequency of the liquid-like granular film increase with granular jet velocity increasing, and are significantly affected by particle diameter. The interesting behaviors of the liquid-like surface waves are observed on the granular sheet. The surface waves become remarkable with the increase of the granular jet velocity, and their propagating velocities normalized by the granular jet velocity vary from 0.7 to 0.9. The waves propagating on the granular sheet may emerge, which will reduce the frequencies of the surface waves and increase the surface wavelengths. The results also show that the oscillation frequency of the granular film nearly equals the pulsation frequency of the granular jet. It is indicated that the gas-solid interaction inside the nozzle increases with granular jet velocity increasing, and causes the instability of the granular jet, resulting in the non-axisymmetric oscillation on the granular sheet consequently. The results in this paper present significant knowledge of the dense granular impinging jets and also provide some principles for the applications in dense granular impinging jets in industrial processes.

Elastic behavior of glass-rubber mixed particles system

Zhao Zi-Yuan, Li Yu-Jun, Wang Fu-Shuai, Zhang Qi, Hou Mei-Ying, Li Wen-Hui, Ma Gang
Acta Physica Sinica. 2018, 67 (10): 104502 doi: 10.7498/aps.67.20172772
Full Text: [PDF 731 KB] Download:(21)
Show Abstract
The mixture of scrap rubber particles and sands has been extensively used as geotechnical engineering recycled materials due to its environmental protection performance, light quality and excellent energy dissipation capability. The mechanical properties of the system can be modulated by the mixing ratio between soft and hard components. But the reasons for such a change on a particle scale are not yet clear. In this paper the elastic behaviors of glass-rubber mixed particles are studied by the sound velocity measurement and discrete element simulation. The velocity of compressional wave and the dynamic effective elastic modulus of mixed sample under hydrostatic stress are measured by time-of-flight method. It is found that the wave velocity is almost constant and the modulus decreases slightly with the proportion of rubber particles increasing to 20%. After that the wave velocity and modulus decrease rapidly and the system transforms from rigid-like behavior to soft-like behavior until the proportion of rubber particles reaches to 80%. When the proportion of rubber particles are more than 80%, the compressional wave velocity and the dynamic effective elastic modulus remain stable again. Such experimental results are consistent with discrete element method analyses which provide more in-depth insights into the micromechanics of the mixture. The simulation reveals that at low rubber fraction the main force chain structure is basically composed of glass particles without rubber particles, which accounts for the phenomenon that the velocity of the compressional wave is basically constant. When the glass particles and rubber particles co-construct the main force chain structure, the distribution of the normal contact force is relatively uniform at high rubber fraction. This can be regarded as the glass particles suspending in the rubber particles. An improved effective medium theory is proposed to describe the elastic behavior of the mixed particles system. It is considered that the deformation of the internal particles is relatively uniform for glass dominated mixture which satisfies the isostress hypothesis. A parallel spring model can be used to describe the nonlinear contact model of particles in such materials. On the other hand, rubber dominated mixture approximately satisfies the isostrain hypothesis, which can be described by a series spring model. The outcomes of such models are in agreement with the simulation results for rigid glass dominated mixture and soft rubber dominated mixture. This study is helpful in exploring the mechanisms that are responsible for the macroscale elastic behavior of mixed granular material from the microscopic point of view.
Acta Physica Sinica
Accepts
Note: The papers published below will continue to be available from this page until they are assigned to an issue. To see an article, click its [PDF] link. To review many abstracts, check the boxes to the left of the titles you want, and click the 'Selected articles' button. To see one abstract at a time, click its [Abstract] link.
»

The Effect of Collision Parameter on a Magnetized Electronegative Plasma Sheath Structure

Accept: 2016-10-11
Full Text: PDF (0KB) (107)
Show Abstract
The structure of an electronegative plasma sheath in an oblique magnetic field is investigated. More over, the collisions between positive ions and neutral particles are taken into account. It is assumed that the system consists of hot electrons, hot negative ions and cold positive ions. Also the negative ions and the electrons are assumed to be described by the Boltzmann distributions of their own temperatures, and the accelerated positive ions are treated by means of the continuity and momentum balance equations through the sheath region. In addition, the assumption that the collision cross section has a power law dependence on the positive velocity is introduced. After theoretical derivation, an exact of sheath criterion is obtained. The numerical simulation results include the distributions of the positive ions density for different invariable ion Mach number satisfying Bohm criterion, the comparison of net space charge distributions for variable and invariable ion Mach number. Furthermore, three species of charged particles density, the net space charge and the spatial electric potential in the sheath are studied numerically for different collision parameters under the condition of the fixed ion Mach number. The results show that the ion Mach number has not only the lower limit but also the upper limit. The ion Mach number affects the sheath structure by influencing the distribution of the positive ion density, and different conclusions can be obtained because ion Mach number is adopted as variable or invariable value while discussing the effects of the other variables which can result in the variety of the ion Mach number on the sheath formation. The reason is the actual sheath structure modification brought on by the variation of a parameter can be resolved into two parts. One is the sheath formation change caused directly by the variation of the parameter, the other is the sheath formation change caused by the Bohm criterion modification which the variation of the parameter results in. Therefore, an identical ion Mach number should be adopted when researching the direct effects of a parameter variety on plasma sheath structure. In addition, it is concluded that the collisions between positive ions and neutral particles make positive ions density curve higher and electrons’ lower than the case without collisions. Negative ions density does not alter significantly whether there exists collision or not. Besides there is a peak in the profile of the net space charge while in the presence of ion-neutral collision and the net space charge peak moves toward the sheath edge. The spatial potential increases and the sheath thickness decreases on account of the presence of the collisions between ions and neutral particles.
»

Calculation of Hamilton energy function of dynamical systems by using Helmholtz theorem

null
Accept: 2016-10-11
Full Text: PDF (0KB) (95)
Show Abstract
The Helmholtz theorem confirmed that any vector field could be decomposed of gradient and rotational field. The supply and transmission of energy occur during the propagation of electromagnetic wave accompanied by variation of electromagnetic field, thus the dynamical oscillators and neurons can absorb and release energy in presence of complex electromagnetic condition. Indeed, the energy in nonlinear circuit is often time-varying when the capacitor is in charged or discharged, and occurrence of electromagnetic induction is available. Those nonlinear oscillating circuits can be mapped into dynamical systems by using scale transformation. Based on mean field theory, the energy exchange and transmission between electronic field and magnetic field could be estimated by appropriate nonlinear dynamical equations for oscillating circuits. In this paper, it investigates the calculation of Hamilton energy for a class of dimensionless dynamical systems based on Helmholtz’s theorem. Furthermore, scale transformation could be used to develop dynamical equations from the realistic nonlinear oscillating circuit, so the Hamilton energy function could be approached effectively. These results could be much useful for self-adaptive control of dynamical systems.
»

Ballistic thermal rectification in the three-terminal graphene nanojunction with asymmetric connection angles

null
Accept: 2016-10-11
Full Text: PDF (0KB) (66)
Show Abstract
By using the nonequilibrium Green’s function method, the ballistic thermal rectification in the three-terminal graphene nanojunction is studied. The dynamics of atoms are described by the interatomic fourth-nearest neighbor force-constant model. The nanojunction has a Y-shaped structure, created by a combination of a straight graphene nanoribbon and a leaning branch as the control terminal holding a fixed temperature. No heat flux flows through the control terminal. There exists a temperature bias between the two ends of the graphene nanoribbon served as the left and right terminals, respectively. The primary goal of this paper is to demonstrate that the ballistic thermal rectification can be introduced by the asymmetric structure with different connection angles between terminals. The control terminal has a smaller connection angle with respect to the left terminal than to the right terminal. The forward direction is defined as being from the left terminal to the right terminal. The results demonstrate that, given the same control temperature and absolute temperature bias, the heat flux in the graphene nanoribbon tends to run preferentially along the forward direction. When the difference between the connection angles increases, the rectification ratio rises. Compared to the zigzag graphene nanoribbon, the rectification ratio of the armchair nanoribbon is more sensitive to the direction the control terminal. However, the greatest rectification ratio is found in the zigzag graphene nanoribbon which has a connection angle of 30 degrees with respect to the armchair branch. In addition, the direction of the control terminal can be adjusted to raise more than 50% of the rectification ratio of the graphene thermal recti?er based on the width discrepancy between the left and right terminals. The mechanism of the ballistic thermal recti?cation is also discussed. In the three-terminal graphene nanojunction, a smaller connection angle with respect to the control terminal leads to more phonon scattering. The confirmation of this conclusion comes from a comparison of phonon transmission between different couples of terminals, which shows that, in most of the frequency spectrum, the phonon transmission between the control terminal and the left terminal is smaller than that between the control terminal and the right terminal. Given the same control terminal temperature and temperature bias, the asymmetric connection angles therefore will introduce a higher average temperature of the left and right terminals, and a larger heat flux in the forward process. Moreover, the average temperature difference between in the forward process and in the reverse process is found to be proportional to the temperature bias, and the proportionality coefficient will get bigger if the asymmetry is strengthened.
»

The Propagation Properties of Vortex Beams in a Ring Photonic Crystal Fiber

null
Accept: 2016-10-11
Full Text: PDF (0KB) (19)
Show Abstract
In the last decade, the vortex beams have received lots of attention for their orbital angular momentum. When they are applied to optical fiber communication field, the data channels will increase and information propagation speed will be effectively improved. Recently, researchers have shown the capability of long length stably propagation, nonlinear frequency conversion and mode division multiplexing of vortex modes in a ring fiber. Due to the photonic crystal fiber (PCF) has very flexible design degrees of freedom, it will enable a wide range of propagation properties. In this paper, A SiO2 air-holes ring PCF is proposed for separation and propagation of optical vortex modes. By using COMSOL Multiphysics software, the vortex modes(TE01, HE_21^± and TM01) are simulated and calculated. The differences of the effective refractive index between them are 4.59×〖10〗^(-4) and 3.62×〖10〗^(-4) respectively. One can analyze the propagation properties of vortex beams in the ring PCF by changing the size of first layer air holes’ radius and air hole pitch. When the incident light wavelength of TE01 mode ranges from 1650 nm to 1950 nm, this ring PCF can achieve a total dispersion variation between 44.18 to 45.83 ps?nm^(-1)?km^(-1), which is tend to be flat. When incident light wavelength is 1550 nm, the nonlinear coefficient of TE01 mode vortex light is 1.37 W^(-1)?km^(-1); Due to the long wavelength light is easier to leakage through the cladding than the short wavelength light, the confinement loss increases with the wavelength. When incident light wavelength is 2000 nm, there is still an eight-orders-of-magnitude of the low confinement loss. Theoretically, flat dispersion and low loss vortex beams in this fiber can be beneficial to propagate stably, and the vortex modes lay the foundation for long distance propagation in the optical fiber. In the future, this ring PCF will be used in optical fiber communication field and application in aspects such as continuous spectrum research, which can make it have immense advantage to traditional fibers.
»

Penta-decomposition of instantaneous field in spanwise-rotating turbulent plane Couette flow

null
Accept: 2016-10-11
Full Text: PDF (0KB) (174)
Show Abstract
Spanwise-rotating turbulent plane Couette flow (RPCF) is one of the fundamental prototypes for wall-bounded turbulent flows in the rotating reference frames. In this turbulent problem, there are large-scale roll cells, which are widely studied. In this paper, a penta-decomposition method is proposed to separate the instantaneous velocity and the total kinetic energy into five parts, including a mean part, a streamwise part and a cross-flow part of the secondary flow, and a streamwise part and a cross-flow part of the residual field, aimed to explore the energy balance and transfer among different shares of the turbulent kinetic energy in RPCF at Reynolds number Rew=Uwh/ν=1300 (here, Uw is the half the wall velocity difference, and h is half channel-height) and rotation number Ro=2Ωzh/Uw (Ωz is the constant angular velocity in the spanwise direction) in the range of 0≤Ro≤0.9. The results show that the energy is transferred between streamwise part (cross-flow part) of secondary flows and residual field through the correlation between the vorticity of secondary flows and shear stress of residual field. The rotation term acts as a bridge to transfer the energy between streamwise part and cross-flow part of secondary flows (residual field). Moreover, pressure-strain redistribution term also plays an important role in the energy transfer between streamwise part and cross-flow part in residual field. For the streamwise part of residual field, in certain rotate rates, the energy obtained from the streamwise part of secondary flows is larger than that got from mean flow, implying that the streamwise motions of secondary flows have a significant impact on the streamwise motions of residual field.
»

A fast particle simulation method for calculating the multipactor threshold based on the frequency domain solutions in microwave devices

null
Accept: 2016-10-11
Full Text: PDF (0KB) (98)
Show Abstract
In order to compute the multipactor thresholds of microwave devices with high ef?ciency and precision, a novel fast particle-in-cell (PIC) method is proposed, which takes advantages of the frequency-domain (FD) electromagnetic field solver of CST Microwave Studio (MWS). At the initial stage of multipactor (when there are not many electrons in the devices), the self-consistent field generated by the electrons is much smaller than the applied electromagnetic field. Therefore it can be ignored in calculating the multipactor threshold and this will significantly reduce the computation burden. During simulations of multipactor processes, the FD fields pre-calculated by CST MWS are converted into time-domain (TD) scaling with the square root of the input power. Then the electrons are advanced by Boris algorithm. When the electrons hit the boundaries of the simulation region, where triangular facets from CST are used for discretization, the secondary electrons would be emitted. After series of simulations with variable input powers, the multipactor threshold is determined according to time evolutions of the electron number. As verifications, the multipactor thresholds in a parallel plate and a coaxial transmission line are investigated. Compared with the results of CST Particle Studio (PS), the fast method obtains almost the same thresholds, while the computational efficiency is improved more than 1 order of magnitude. Since the self-consistent field generated by the electrons is ignored in the fast method and it is considered in CST PS, the results validate that the self-consistent field can be ignored in calculating the multipactor threshold. Finally, taking a parallel plate transmission line and a stepped impedance transformer as examples, we studied the effects of the number of initial macro-particles on the calculation precision. When the initial particles are so few that it can hardly reflect the randomness of the multipactor process, it results in a higher calculated value. With the increase of the number of initial macro-particles, the calculated multipactor threshold is lower and more accurate. It is convergent when the number reaches about 2000 for the parallel plate transmission line and 4000 for the stepped impedance transformer, respectively. Taking into account other microwave devices with more complex electromagnetic field distribution, in order to ensure precision, it is recommended to select the number of initial macro-particles 8000. In addition, although CST MWS was used to obtain the electromagnetic fields and boundary information in this paper, of course, other electromagnetic software (such as HFSS) can also be adopted as an alternation.
»

The effect of linear bubble vibration on wave propagation in unsaturated porous media containing air bubbles

null
Accept: 2016-10-11
Full Text: PDF (0KB) (204)
Show Abstract
Biot model is widely applied in geophysics, petroleum engineering, civil engineering and ocean engineering since it has been presented. This leads to a considerable development of the research on the wave propagation in saturated porous medium. However, fully saturated porous medium is rarely found in nature, almost all the rocks or soils contain two kinds of fluid, such as gas and petroleum. So many researches has been done on the wave propagation in unsaturated porous medium by domestic and abroad scholars. It is well known that the presence of a small volume of gas bubbles in a liquid can greatly alter the velocity and attenuation of acoustic waves in the liquid. Evidence is beginning to accumulate that the velocity and attenuation of acoustic waves in a saturated marine sediment can be affected by the presence of gas bubbles in the saturating liquid. To investigate the sound propagation in porous media when the pore water contains a small amount of air bubbles, this paper integrates the volume vibration of bubbles in pore water into the continuity equation of pore-fluid filtration in porous media based on Biot theory, so as to obtain the continuity equation of pore-fluid filtration with bubble volume vibration. On this basis, according to the relationship between the instantaneous radius of bubbles and the background pressure of the medium under the linear vibration of bubbles, as well as the equations of motion of the fluid medium and porous medium, a new displacement vector wave equation of porous media under the influence of bubbles is derived, which establishes the model for the sound velocity dispersion and attenuation prediction under the unsaturated porous media. The presence of air bubbles increases the compressibility of pore fluid, which leads to the decrease in the sound velocity of the bubbly saturated porous media. When the wave frequency equals to the resonance frequency of the bubbles, the bubbles in pore water will produce resonance; the medium will present to be highly dispersive and the velocity can greatly exceed the gas-free velocity, but these have not been measured in field data; and the absorption cross section of the air bubble can reach the maximum, which leads to the maximum attenuation of the porous media. It should be noted that the attenuation coefficient calculated with this model is related to the damping of bubble motion(radiation, thermal and internal friction) and the dissipation of the relative motion between the pore water and porous solid frame. The obtained numerical analysis is consistent with the above conclusions, which indicates that the volume concentration, the bubble size and the excitation frequency of sound field are important parameters affecting the sound wave propagation in the saturated porous media containing few bubbles.
»

Ferroelectric phase transition of perovskite SnTiO3 based on first principles

null
Accept: 2016-10-11
Full Text: PDF (0KB) (78)
Show Abstract
Due to their spontaneous polarization, ferroelectric materials have excellent dielectric, piezoelectric, pyroelectric and other properties, which enable them to be used in many applications, such as capacitors, filters, sensors, detectors, and transducers, among others. In this paper, we employ a first-principles-based effective Hamiltonian method to investigate perovskite SnTiO$_3$, obtaining essential coefficients for the effective Hamiltonian via ab initio computations, which are used in subsequent Monte-Carlo simulations to predict the phase transition temperature of SnTiO$_3$, and different structural phases involved in such phase transition.
»

Nonlocal Symmetries and Interaction Solutions of the (2+1)-dimensional Higher Order Broer-Kaup System

xiangpeng xin Hanze Liu Xi-qiang LIU
Accept: 2016-10-11
Full Text: PDF (0KB) (210)
Show Abstract
The (2+1)-dimensional higher-order Broer-Kaup (HBK) system is studied by nonlocal symmetry method and consistent tanh expansion (CTE) method. In this paper, via the localization of the residual symmetries, the nonlocal symmetries are localized to Lie point symmetries and symmetry groups are also obtained. Many types of soliton solutions and interaction solutions among different nonlinear excitations such as solitons, periodic waves etc. are constructed. In order to study their dynamic behaviors, corresponding images are explicitly given.
»

Development of a intranuclear-cascade code CBIM applicable to the nuclear reaction with incident particle energy above 45MeV

null
Accept: 2016-10-11
Full Text: PDF (0KB) (72)
Show Abstract
The Monte Carlo intra-nuclear cascade program CBIM has been developed for describing nuclear reactions involving protons, neutrons and pions on complex nuclei. In order to describe cascade process, several simplifications have been made in the following: firstly, neither reaction, reflection, refraction, nor ionization will be taken into account before the incident particle enters the target nucleus; secondly, target nucleus is regarded as spherical and the atom number should be greater than 2; thirdly, the knocked nucleon is determined by cross section sampling; last, in the center-of-mass frame, the scattering angle is sampled based on differential cross section distribution.. The basis physics model bases on the above assumptions and Bertini intra-nuclear cascade model; meanwhile, nucleon-nucleon angle differential distributions of INCL in the center-of-mass frame have been introduced to overcome the shortage of Bertini model. The interactions between nucleon and nucleon or between nucleon and pion, for example, elastic scattering, pion production and charge exchange, are simulated in the code. In the particles collision, the nucleon density changes with the target nucleus radius; and the interaction cross sections refer to 22 kinds of experimental cross sections in Bertini model. The intra-nuclear cascades induced by 45MeV~3500MeV neutron, proton or pion below 2500MeV can be simulated by this code. Finally, comparisons with experiment on reaction cross section over the energy range 60~378MeV, and some simulation results by MCNPX, GEANT4 and PHITS over the energy range 65~3000MeV, the CBIM results are in reasonable agreement with them over the broad energy range considered.
»

Omnidirectional photonic bandgap of the one-dimensional plasma photonic crystal based on a novel Fibonacci quasiperiodic structure

Accept: 2016-10-11
Full Text: PDF (0KB) (228)
Show Abstract
Take the binary one-dimensional plasma photonic crystal based on Fibonacci quasiperiodic structure as an object, on the basis of the photonic bandgap characteristics of the structure with different initial sequence and number of period, a novel structure of one-dimensional plasma photonic crystal is proposed in this paper to enlarge the omnidirectional photonic bandgap (OPBG). Compared with previously reported structures in literatures, this structure is simpler in configuration with fewer layers and materials, and its OPBG width is wider. The influence of the parameters of the plasma material, such as the thickness, plasma frequency and collision frequency, on the OPBG characteristics of this structure is systematically discussed and compared with that of the structure in literatures. The research results can provide important theoretical guidance for the design of novel omnidirectional reflectors.
»

Quantum secure direct communication protocol based on the mixture of Bell state particles and single photons

Zheng-Wen CAO
Accept: 2016-10-11
Full Text: PDF (0KB) (26)
Show Abstract
By studying the properties of the mixture of Bell state particles and single photons, the paper designs a quantum code scheme with high coding capacity, and proposes a novel quantum secure direct communication protocol with high transmission efficiency. Alice prepares Bell state particles and single photons, and divides Bell state particles into two sequences $S_A$ and $S_B$. $S_B$ is sent to Bob for the first security check using quantum correlation properties of particles. When the check result shows that the quantum channel is safe, using designed quantum code scheme, Alice encodes her classical message on the mixed quantum state sequence of Bell sequence $S_A$ and single photon sequence $S_S$. Then, some single photons that are used for security check are re-inserted randomly into the encoded sequence, and the order of particles is rearranged to ensure to check Eve's attack. Alice sends the new sequence to Bob. Bob delays and receives it. And then, the quantum channel is conducted security check for the second time. The transmission error rate is calculated, if the error rate is lower than the tolerance threshold, the channel is safe. Bob decodes and reads Alice's message. The first security check is to determine whether quantum channel is safe. The second security check could test whether there are eavesdroppers during information transmission. Safety analysis is done by using quantum information theory to the proposed protocol. The error rate introduced by Eve and the amount of information by Eve are calculated. It is showed that this protocol can effectively resist measurement-resend attack, intercept-resend attack, auxiliary particle attack, denial of service attack and Trojan attack. Among them, auxiliary particle attack is analyzed in details. The transmission efficiency and coding capacity are also analyzed. The transmission efficiency is 2, the quantum bit rate is 1, and the coding capacity is that a quantum state can encode three bits of classical messages. We also compare the proposed protocol to many existing popular protocols in terms of efficiency, e.g., Ping-Pong protocol, Deng,F.G. et al.'s Two-step and One-pad-time quantum secure direct communication protocol, Wang,J. et al.'s quantum secure direct communication protocol based on entanglement swapping and Quan,D.X. et al.'s one-way quantum secure direct communication protocol based on single photon. It is proved that this proposed protocol has higher transmission efficiency. In addition, complex U operation and entanglement swapping are not used, and implementation process is simplified. However, this protocol is devoted to theoretical research of quantum secure direct communication. There are still some difficulties in the practical application. For example, the storage technology of quantum states is not mature at present. It is not easy to prepare and measure Bell state particles and combine them with single photons, and so on. The implementation of this protocol depends on the development of quantum technology in the future.
»

The Relationship between Dielectric Properties and Nanoparticle Dispersion of Nano- SiO2/Epoxy Composites

null
Accept: 2016-10-11
Full Text: PDF (0KB) (38)
Show Abstract
Nano-SiO2 was modified by silane coupling agent and modified nano-SiO2 powder and nano-SiO2 dispersing liquid was obtained. Unmodified and modified nano-SiO2/Epoxy composites made by “mechanical mixing method”, and modified namo-Silica/Epoxy composites made by “bubble mixing method” were prepared, respectively. The content of nano-SiO2 in the composite is 2wt%, 3wt%, 4wt%, 5wt% and 6wt%. Breakdown strength and corona-resistance characteristics of the composites were tested. The results show that, with the increase of nano-SiO2 loading, the breakdown strength and corona-resistance of nano-SiO2/Epoxy composites increase. The maximum breakdown strength of namo-Silica/Epoxy composites was appeared when the nano-Silica content is 5wt%. The SEM images of 5wt% nano-Silica loading composites were analyzed by Software Image J, and the Morisita’s Index method was used to evaluate the dispersion of nano-Silica particles in the matrix quantitatively. The best dispersion was found in the composites made by “bubble mixing method”. The relationship between dielectric properties and nano-particle dispersions of nano-Silica/Epoxy composites was discussed.
»

Combined noise source identification method based on spherical microphone array with random unifrom distribution of elements

null
Accept: 2016-10-11
Full Text: PDF (0KB) (72)
Show Abstract
As the developing of techlology, noise controlling is paied wide attention in recent years. Noise source identification is the key step for noise controlling. Spherical microphone array, which can located the noise source of arbitrary direction in three dimensional space, is widely used for noise source identification in recent years. Conventional methods for noise source localization include spherical near field acoustic holography and spherical focused beamforming. The acoustic quantities are reconstructed by using spherical near field acoustic holography method to realize the noise source identification, while the noise source can also be located by using focused beamforming based on spherical harmonic wave decomposition. However, both these methods have their own limitations while being used in noise source identification. Spherical near field acoustic holography has low resolution in high frequency with far distance from noise source to measurement array for noise source identification, whereas the spherical focused beamforming has low localization resolution in low frequency. Noise source identification is discussed here and a 64-element microphone spherical array with randomly uniform distribution of elements is designed. The combination methods of noise source identification by using spherical near field acoustic holography and mode decomposition focused beamforming are researched. The performance of the proposed combination methods is simulated, and an experiment of noise source identification is carried out based on the designed spherical microphone array to test the validity of proposed method. The dividing frequency point is when selecting noise source identification methods between near field acoustic holography of spherical wave decomposition by using the spherical array designed in this paper. Research results show that high resolution of noise source identification can be obtained by using near field acoustic holography when reconstruction frequency is with a distance from noise source to the center of spherical array, while high resolution of noise source localization can be achieved by using spherical wave decomposition beamforming when signal’s frequency is with a distance from noise source to the center of spherical array. Spherical array with random uniform distribution of elements maintains stable identification ability in all bearing. Spherical near field acoustic holography has high resolution distinguish ability in near field and low frequency, while focused beamforming method has high resolution distinguish ability in far field and high frequency. Therefore the noise source can be efficiently identified by using the proposed combined method of near field holography and focused beamforming with less elements and small aperture spherical microphone array.
»

Optimization design of a Gamma-to-Electron spectrometer for high energy gammas induced by fusion

null
Accept: 2016-10-11
Full Text: PDF (0KB) (203)
Show Abstract
Apart from neutrons, the fusion core produces gamma rays during one fusion reaction. The spectrum of gamma ray can provide very important information for fusion diagnosis. However, due to the gamma energy and yield in one fusion pulse, the gamma spectrometer used should have high detection efficiency and energy resolution. The concept of a Gamma-to-Electron magnetic spectrometer GEMS provides the idea to build up such a spectrometer to meet this requirement. Based on this concept design, four important parts of this facility are investigated. The first part is the gamma-electron converter. The main physics processes include Compton scattering of gamma ray with converter material generating electron, the electron Multiple Coulomb scattering (MCS) inside the converter and the electron attenuation. Affected by the thickness of convector, these processes gives a complex influence on the detection efficiency and angular-energy distribution of the electrons which are emitted from the downstream face of the convector. The Monte Carlo code Geant4 is employed to investigated the functions of Compton scattering, MCS and converter thick on the angular-energy distribution. The second one is the collimation. The collimation is used to select the forward direction election, the performance of cutoff angle of the collimator on the detection efficiency and resolutions, as well as the correlation between electron transportation direction and energy, are also studied using Geant4 code. The third part is the dipole magnetic field. There are several parameters of geometric and magnetic, therefore, a multi-thread parallelized Genetic algorithm is developed to get the best result. Both the irregular geometry (shape) and dipole magnetic field strength are optimized to achieve the best energy resolution and detection efficiency. The obtained magnetic field has intensity less than 100 Gauss, and its performance on gathering elections is also verified by Geant4 code. The last one is the location of electron detectors. The study shows that all the electron detectors should be located according to not a straight line but a quadratic curve. Then the optimized spectrometer is simulated by Geant4 to get the responses of gamma rays with various energies. For the gammas provided by fusion reaction, the simulation shows that when the neutron yield is about 2.5×1015 and 1.2×1016, the energy resolution reaches 0.5 MeV and 0.25 MeV, respectively, provided that different thick Be converters are employed. All in all, this optimized GEMS can be employed to measure the spectrum of gamma rays generated by the fusion reaction.
»

Influnence of Nonspherical Effects on the Secondary Bjerknes Force in a Strong Acoustic Field

null
Accept: 2016-10-11
Full Text: PDF (0KB) (78)
Show Abstract
The secondary Bjerknes force between bubbles in an acoustic field is a well-known acoustic phenomenon. The major theoretically researches of the secondary Bjerknes force were owing to two spherical bubbles. The secondary Bjerknes force between two spherical bubbles which is calculated based on the linear equations is very small and negligible, therefore these theoretically researches did not give a well explanation for the phenomenon, such as “streamer formation” and Multi-bubble sonoluminescence (MBSL). Experiments of sonoluminescence (SL) show that bubbles in a sound field are not entirely spherical bubbles. Nonspherical effects have an important influence on the secondary Bjerknes force when two bubbles come close to each other in a strong acoustic field (>1.0×〖10〗^5 Pa). How does the shape distortion of a nonspherical bubble cause the change of the secondary Bjerknes force between two bubbles, and the secondary Bjerknes force how to affect the oscillation and movement of bubbles are major problems which we wish to solve. The of the secondary Bjerknes force between a nonspherical bubble and a spherical bubble is obtained by considering the shape oscillation of a nonspherical bubble. We numerical simulate the secondary Bjerknes force between a nonspherical bubble and a spherical bubble based on the nonlinear oscillation equations of two bubbles, and compare the secondary Bjerknes force between a nonspherical bubble and a spherical bubble to the secondary Bjerknes force between two spherical bubbles in the same condition. We discuss the influence of nonspherical effects on the secondary Bjerknes force between two bubbles. The results show that when the amplitude of driving pressure is greater than the Blake threshold of a nonspherical bubble and makes the bubble oscillate stably, the secondary Bjerknes force between this nonspherical bubble and a spherical bubble is different to the secondary Bjerknes force between two spherical bubbles in direction and magnitude. The secondary Bjerknes force between a nonspherical bubble and a spherical bubble is much bigger than that of two spherical bubbles. The interactional distance of the secondary Bjerknes force between a nonspherical bubble and a spherical bubble is further than that of two spherical bubbles. The secondary Bjerknes force between a spherical bubble and a nonspherical bubble depends on the radii of two bubbles, distance between two bubbles, shape mode of the nonspherical bubble and the amplitude of driving pressure. Our research is more close to the actual bubbles in liquid. We also prove that big mutual interaction between bubbles is mainly cause for the formation of a stable structure between bubbles. For bubbles, big mutual interaction causes the cavitation become easier. These results are important to explain the phenomenon in an acoustic field, such as “streamer formation” and Multi-bubble sonoluminescence (MBSL).
»

The Principle and Application of Diagonal Reducing Method in the Complex Noise Fields

null
Accept: 2016-10-11
Full Text: PDF (0KB) (208)
Show Abstract
Acoustic environment has low signal-to-noise ratio (SNR); hence, array signal processing is always used for noise reduction and signal enhancement. Because the delay-and-sum beamforming method performs robust, so it is almost widely used, but the array gain is limited by the array aperture. The actual underwater ambient noise is complex, which includes uncorrelated noise and correlated noise. The noise power of each array element is unequal. The noise covariance matrix is not a scaled identity matrix. Consequently, the performance of array signal processing method decreases obviously. Aiming at these two problems, the diagonal reducing method of the covariance matrix in the complex noise fields is proposed. Firstly, a reducing matrix, which is defined as a diagonal matrix with unequal diagonal elements, is subtracted from the covariance matrix so as to reduce the noise, and a new matrix is obtained. Secondly, the delay-and-sum beamforming is done by using the new matrix to obtain the beaming output. The analytic solution and approximate solution of reducing matrix are obtained under the constraint condition that the output SNR attains its maximum. Thirdly, the estimation of the reducing matrix is determined by minimizing the function that is defined as the error between the covariance matrix and the estimated covariance matrix. This minimization problem is accomplished in an iterative method. Fourthly, if the noise is uniform white noise or the nonuniform white noise, this proposed method performs well. While, under the complex noise field the performance of the proposed method may be deteriorated. So the effects of the correlation of the noise field and the input SNR on the estimated error is analyzed. In fact, the weaker the correlation is, or the smaller the input SNR is, the smaller the estimated error is. Lastly, the simulation experiment and the lake trial are implemented. The simulation results show that the diagonal reducing method of the covariance matrix reduces some ambient noise, the noise output power is decreased, the output SNR is increased, and the proposed method improves performance of array signal processing. The experimental results show that the output SNR of the target using the proposed method is increased by about 14 dB. The diagonal reducing method of covariance matrix has definite value to engineering application, and is computationally attractive.
»

Moving target compressive imaging based on improved row scanning measurement matrices

null
Accept: 2016-10-11
Full Text: PDF (0KB) (64)
Show Abstract
Abstract: Moving target imaging (MTI) plays an important role in practical applications. How to capture dynamic images of the targets with high quality is a front-burner issue in the field of MTI. In order to improve the reconstruction quality, a new MTI model based on compressed sensing (CS) is proposed here, applying a sampling protocol of the row-scanning together with a motion measurement matrix constructed by our own. It is proved by the simulation and the experimental results that a relatively higher quality can be achieved through this approach. Furthermore, an evaluation criterion of reconstructed images is introduced to analyze the relationship between the imaging quality and the moving speed of the target. By contrast, the performance of our algorithm is much better than that of traditional CS algorithms under the same moving speed condition. As a result, it suggests that our imaging method may have a great application prospect in the earth observation of unmanned aerial vehicles, video monitoring in the product line and other fields.
»

Spatial Correlation of Underwater Bubble clouds Based on Acoustic Scattering

null
Accept: 2016-10-11
Full Text: PDF (0KB) (210)
Show Abstract
Using effective medium theory to describe acoustic scattering from bubble clouds, one of the underlying assumptions shows that the probability of an individual bubble being located at some position in space is independent of the locations of other bubbles. However, bubbles within naturally occurring clouds are usually influenced by the motion of the fluids which makes they become preferentially concentrated or clustered. According to Weber’s method, it is a useful way to importing spatial correlation function to describe this phenomenon in bubble clouds. The spatial correlation function is contained in acoustic scattering and it is important to notice that the spatial correlation should be dependent of the position and radius of each bubble due to the ‘‘hole correction’’ or the effect of the dynamics of the fluids. Because of these reasons, it is hard to invert the spatial distribution of bubble clouds using spatial correlation function in acoustic scattering. A method is described here in which bubble clouds are separated into many small subareas and the conception called effective spatial correlation function which is the statistic of spatial correlation function used to describe the correlation between each subarea of bubble clouds. Since the effective spatial correlation function is independent of bubbles’ radius and positions, the bubble clouds’ distribution and the trend of clustering can be inverted by using this function. The result of simulation indicates that the effective spatial correlation function can precisely track the position of the clustering center, even the clustering center covered by other bubble clouds can be detected. Using multi-bean sonar measuring the bubbly ship wake generated by a small trial vessel, the method is used to invert the spatial distribution and clustering centers of bubble field in the ship wake. The results show that effective spatial correlation function accurately inverts the distribution and clustering centers of bubbles in ship wake. Furthermore, the method presented in this paper could distinguish the bubble clouds caused by different reasons and detect upper ocean bubble clouds covered by other bubbles generated by wave breaking.
»

Uncertainty Quantification in the Calculation of keff Using Sensitity and Stochastic Sampling method

null
Accept: 2016-10-11
Full Text: PDF (0KB) (37)
Show Abstract
In the neutronics simulation of nuclear reactor, the uncertainties associated to the integral parameters due to the uncertainties in nuclear data are usually quantified using the sensitivity and uncertainty (S/U) analysis method based on the perturbation theory. S/U analysis method is only applicable to the linear model, moreover neutronics code generally can not be directly used in sensitivity analysis. Sampling approach, which evaluating the uncertainties by performing a set of stochastic simulations, is easy to implement and the uncertainties quantified is close to exact. The function of uncertainty quantification based on sampling approach have been added to uncertainty analysis code SURE. Before applying the sampling method to the uncertainty quantification in the simulation of complex problems, it is necessary to carry out a careful verification. The uncertainties of the calculated effective neutron multiplication factor keff for two selected simple critical benchmark experimental model are quantified using SU method and sampling method respectively. The keff uncertainties due to all nuclides and reaction types nuclear data quantified by two methods are in good agreement, and the correctness of the sampling function of SURE code is verified. The keffs distributions from sampling method obey normal distribution, which embodies a linear relation between input nuclear data and output keff in the range of the uncertainty range of nuclear data, and sensitivity analysis method is adaptable to quantify uncertainty of calculated keff.
»

A super-resolution infrared microscopy based on a doughnut pump beam

null
Accept: 2016-10-11
Full Text: PDF (0KB) (208)
Show Abstract
This paper presents an approach to break through the diffraction limitation in infrared microscopies. In this method, instead of Gaussian pump beam, an intensive vortex beam is firstly focused on the sample, leading to saturation absorption of the peripheral molecules in the point spread function (PSF). The vortex beam is followed by a Gaussian beam with the same wavelength, which can only be absorbed by the molecules near the center, resulting in shrunken PSF which means higher resolution. Furthermore, the PSF of a system based on this approach is numerically simulated. With an 100 nJ pulse energy vortex beam and a 0.1 nJ pulse energy probe beam, the theoretical resolution (full width at half maximum, FWHM) is measured to be about 236 nm which is 14 times better than that of the traditional infrared microscopy.
»

Fast Bayesian Blind Restoration for Single Defocus Image with Iterative Joint Bilateral Filters

null
Accept: 2016-10-11
Full Text: PDF (0KB) (57)
Show Abstract
It is significant to realize effective defocus image restoration for acquiring clear image in military and geological examination field. Most of existing algorithms have the problems of large computational cost, ringing and noise sensitivity, hence a novel approach by iterative joint bilateral filtering under Bayesian framework is proposed. Firstly, it utilizes defocus image depth estimation to compute the point spread function in the Bayesian framework. Then a minimum optimization problem is built to represent the blind restoration problem. After inferring the solution procedure of the minimum optimization problem, we find that the joint bilateral filters can be used to search the optimal solution, which not only simplify the searching procedure but also reduce the computational cost. Finally, an iterative joint bilateral filtering was designed to realize the image restoration. That means the original restored image obtained from the bilateral filtering is used to design the guide image for the joint bilateral filters, and the guide image will serve as the input of the optimization problem for acquiring the better optimal result. This procedure was repeated until convergence. The experiment results indicate that this method can yield the ringing, reduce the computational cost and remove the noise. Generally speaking, the average pixel error of 85% images is under 0.03, which has improved 19% comparing with the same error rang of existing algorithms. And 78% shorter than those of compared algorithms. It can be used in the engineering practice of blind restoration for single defocus image.
»

First-principles study on the thermodynamic stabilities and electronic structures of long-period stacking ordered phases in the Mg-Y-Cu alloys

null
Accept: 2016-10-11
Full Text: PDF (0KB) (206)
Show Abstract
A first-principles method based on density functional theory has been used to investigate thermodynamic stability and electronic characteristics of long-period stacking ordered (LPSO) phases 14H and 18R (18R(m),18R(t)) in Mg–Y–Cu alloys. The present calculations are performed using Vienna Ab-initio Simulation Package (VASP) with projector augmented plane wave pseudopotential, and generalized gradient approximation is used to treat with and describe the exchange-correlation interaction. The plane wave cutoff energy is set to 360 eV, the forces on all the atoms is less than 0.02 eV/?. The calculated negative enthalpies of formation show that both 14H and 18R can exist in Mg–Y–Cu system, 14H and 18R are stable with respect to the Mg, Cu and Y elements, the reaction energies indicate that 14H is more stable than 18R. The density of states (DOS) of these phases reveals that the main bonding peaks of 14H is located at energy range between -6.82 eV and 2.09 eV, those of 18R(m) at energy range between -6.82 eV and 2.02 eV, and 18R(t) at energy range between -6.82 eV and 1.98 eV. The Cu 3d orbits, Y 4d orbits, Mg 3s and Mg 2p orbits are broadly distributed in the entire region, while Cu 4s orbits, Y 4s and Y 4p orbits are very weak in whole region. For 14H,18R(m) and 18R(t) phases, the bonding originates mainly from the valence electrons of Mg 3s, Mg2p, Cu 3d and Y 4d orbits. The presence of pseudogap indicates that the bonds in 14H and 18R phases are noticeable covalent. In addition, the charge density on (0 0 0 1) plane of 14H and 18R phases are analyzed, and the results indicate that the Cu-Y bonds exhibits covalent feature in 14H and 18R, the covalent bonding of 14H phase is stronger than that of 18R phase.
»

A broadband low-frequency sound insulation structure based on two-dimensionally inbuilt Helmholtz resonators

null
Accept: 2016-10-11
Full Text: PDF (0KB) (58)
Show Abstract
A man-made acoustic structure with broadband low-frequency sound insulation property is designed based on circularly inbuilt Helmholtz resonators. Beyond this structure, a two-dimensional quiet zone can be created. Being the same as the simulating model, an experimental structure is fabricated. Experiments are carried out to study its sound insulation properties. The experimental results are very coincident with the simulating one, which show that this structure has an excellent sound insulation effect in the frequency band of 680-1050Hz, and the maximum insulation sound pressure level can reach 41dB. Meanwhile, the distribution of the two-dimensional sound field above this structure is measured. The results point out that the range of the insulation area can be changed with the change of the incident frequency. In addition, the sound insulation effect is sensitive to the resonant state of the Helmholtz resonators. This work will be of help for designing new sound protection devices.
»

Total Dose Dependence of Hot Carrier Injection Effect in the NMOS Devices

null
Accept: 2016-10-11
Full Text: PDF (0KB) (212)
Show Abstract
The equipments and devices which were long-time running in space were affected by space radiation effects and hot carrier injection effects at the same time which would reduce their optional times. Normally, the single mechanism test simulation method was used on the ground simulation test but the multi-mechanism effects was affected the space equipments and devices, included total irradiation dose effect, hot carrier injection effect, and so on. The total dose dependence of hot carrier injection (HCI) effect in the 0.35μm NMOS Devices was studied in this paper. Three samples were test with different conditions (sample 1# with TID and HCI test, sample 2# with TID, annealing and HCI test, sample 3# only with HCI test). The results shows that threshold voltage of NMOS devices with 5000s HCI test after 100krad (Si) total dose radiation shift negatively then positively during total dose irradiation test and HCI test,and it was more than the devices without radiation test. But the threshold voltage shift of NMOS devices with 5000s HCI test and 200hours annealing test after TID test was more than the devices without radiation test and lower than the devices without annealing test. That was, the parameters of NMOS device varied faster with the association effects (included total dose irradiation effect and HCI effect) than with single mechanism effect. It was indicated that the hot electrons were trapped by the oxide trap charges induced by irradiation effect and then became recombination centre. And then the oxide trap charges induced by irradiation effect reduced and became to negative electronic. The interface trap charges induced by irradiation effect were reduced and then increased and it was because that the electrons of hole-electron pairs in the Si-SiO2 interface were recombined by oxide traps in the oxide during the forepart of HCI test but then the electrons were trapped by interface traps in the Si-SiO2 interface because the electrons from source area were injected to interface during the HCI test. So the threshold voltage shift was positively due to the negative oxide trap charges and interface trap charges. The association effect was attributed to the reduction of oxide traps induced by recombination with hot electrons and the increase of the interface traps induced by irradiating trapped with hot electrons.
»

Optoelectronic properties of N/B doped graphene

null
Accept: 2016-10-11
Full Text: PDF (0KB) (205)
Show Abstract
Since its discovery in 2004, graphene has attracted great attention because of its unique chemical bonding structure, which has excellent chemical, thermal, mechanical, electrical and optical properties. Due to the zero band gap material, graphene has limited its development in the field of Nano Electronics. Only expanding the band gap of the graphene can promote the application of graphene in Nano Electronics. In this paper, we constructed three models of intrinsic graphene, N-doped graphene and B-doped graphene. The energy band structures, electronic density of states and optical properties of N/B doped graphene with intrinsic graphene and different doping concentrations were studied. The absorption spectra, the reflection spectra, the refractive index, the conductivity and the dielectric function were studied. The study shows that the electronic states near the Fermi level of N/B doped graphene are mainly composed of C-2p and N-2p/B-2p orbitals, and N/B doping can induce the change of the Fermi level and the photoelectric properties of graphene. The conclusion of this paper can provide a theoretical basis for the application of graphene in optoelectronic devices.
»

The study of the dynamic of the slow electrons transmitted through straight glass capillary and tapered glass capillary

null
Accept: 2016-08-18
Full Text: PDF (0KB) (8)
Show Abstract
It was found that the transmission rate of the electrons through insulating capillaries as a function of the time/incident charge is not the same as that for the ions. The question arouse that, by using the electrons, if the negative charge patches can be formed to facilitate the transmission of the followed electrons, substantiating that the so-called guiding effect works also for electrons. This study aims to observe the time evolution of the transmission of electrons through a straight glass tube and a tapered glass capillary. This would reveal the details that how and/or if the negative charge patches can be formed when the electron are being transported through them. In this work, a set of MCP/phosphor two-dimensional detection system based on Labview platform was developed to obtain the time evolution of the angular distribution of the transmitted electrons. The pulsed electron beams through a small hole with the diameter of 0.5 mm was obtained to test our detection system. The time evolution of the angular profile of 1.5 keV electrons transmitted through the glass tube/capillary was observed. The transmitted electrons are observed on the detector for a very short time and disappear for a time and then back again for both the glass tube and tapered glass capillary, leading to an oscillation. The positive charge patches are formed in the insulating glass tube and tapered glass capillary since the secondary electron emission coefficient for the incident energy is larger than 1. It is due to the fast discharge of the deposited charge, leading to an increase of the transmission rate, while the fast blocking of the incident electrons due to the deposited positive charge, leads to a decrease of the transmission rate. The geometrical configuration of the taper glass capillary tends to make the secondary electrons deposited at the exit part to form the negative patches that facilitate transmission of electrons, similar to the guiding of positive charged ions. This suggests that if the stable transmission needs to be reached for the production of the electron micro-beam by using tapered glass capillaries, the steps has to be taken to have the proper grounding and shielding of the glass capillaries and tubes. Our results show a difference for electrons in transmission through the insulating capillary from that of highly charged ions.
»

aaaa

null
Accept: 2016-08-18
Full Text: PDF (0KB) (17)
Show Abstract
bbbbb
»

internuclear-distance-dependent ionization of H$_2^+$ in strong laser fields in a classical perspective

null
Accept: 2016-08-18
Full Text: PDF (0KB) (12)
Show Abstract
The enhanced ionization of H$_2^+$ in strong laser fields is studied by numerically simulating the classical Hamiltonian equation with the fix-nuclei approximation. The classical trajectory of the electron shows the electron gains energy from the laser field by circulating one electron, then passes through the interatomic barrier and move around the other nucleus before ionization. The ionization probability is maximum when the energy difference between the ground state and the the higher value of the interatomic barrier and outatomic Coulomb barrier is minimum. The classical calculation offers a perspective to inspect the intriguing phenomena in quantum systems.
»

Partition and growth of convection patterns in Poiseuille-Rayleigh-Benard flow

null
Accept: 2016-08-18
Full Text: PDF (0KB) (12)
Show Abstract
In this paper, the Simple algorithm is used to numerically simulate the two-dimensional fully hydrodynamic equations. Partition of convection pattern , growth and the effect of horizontal flow on the characteristical parameters of different patterns in Poiseuille-Rayleigh-Benard flow are studied. The result indicated that flow zone is divided into three zones by the upper and lower critical Reynolds numbers , such as traveling wave zone, localized traveling wave zone, horizontal flow zone.and increase with reduced Rayleigh number. In the growth stage of the convection pattern, the growth process of three kinds of patterns with time is different, but the convection rolls grow all from downstream; Variation of characteristic parameters with time is also different, maximum vertical velocity and Nusselt number of traveling wave and localized traveling wave enter into the stable stage of the cycle variation after the exponential growth stage;and of horizontal flow pattern down to a stable constant after slow growth. and of three types of patterns decrease with increasing Reynold number, there are different rules in the different pattern areas. In this paper, formulas on variation ofandwith and formulas on variation ofandwithin different convection patterns are suggested.
More
Acta Physica Sinica
2018 Vol.67      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9
2017 Vol.66      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13      No.14      No.15      No.16      No.17      No.18
     No.19      No.20      No.21      No.22      No.23      No.24
2016 Vol.65      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13      No.14      No.15      No.16      No.17      No.18
     No.19      No.20      No.21      No.22      No.23      No.24
2015 Vol.64      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13      No.14      No.15      No.16      No.17      No.18
     No.19      No.20      No.21      No.22      No.23      No.24
2014 Vol.63      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13      No.14      No.15      No.16      No.17      No.18
     No.19      No.20      No.21      No.22      No.23      No.24
2013 Vol.62      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13      No.14      No.15      No.16      No.17      No.18
     No.19      No.20      No.21      No.22      No.23      No.24
2012 Vol.61      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13      No.14      No.15      No.16      No.17      No.18
     No.19      No.20      No.21      No.22      No.23      No.24
2011 Vol.60      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2010 Vol.59      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2009 Vol.58      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13
2008 Vol.57      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2007 Vol.56      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2006 Vol.55      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2005 Vol.54      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2004 Vol.53      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2003 Vol.52      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2002 Vol.51      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2001 Vol.50      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
2000 Vol.49      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1999 Vol.48      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
     No.13
1998 Vol.47      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1997 Vol.46      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1996 Vol.45      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1995 Vol.44      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1994 Vol.43      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1993 Vol.42      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1992 Vol.41      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1991 Vol.40      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1990 Vol.39      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1989 Vol.38      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1988 Vol.37      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1987 Vol.36      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1986 Vol.35      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1985 Vol.34      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1984 Vol.33      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1983 Vol.32      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1982 Vol.31      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1981 Vol.30      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1980 Vol.29      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1979 Vol.28      No.1      No.2      No.3      No.4      No.5      No.6
1978 Vol.27      No.1      No.2      No.3      No.4      No.5      No.6
1977 Vol.26      No.1      No.2      No.3      No.4      No.5      No.6
1976 Vol.25      No.1      No.2      No.3      No.4      No.5      No.6
1975 Vol.24      No.1      No.2      No.3      No.4      No.5      No.6
1974 Vol.23      No.1      No.2      No.3      No.4      No.5      No.6
1973
1972
1971
1970
1969
1968
1967
1966 Vol.22      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9
1965 Vol.21      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1964 Vol.20      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1963 Vol.19      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1962 Vol.18      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1961 Vol.17      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1960 Vol.16      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8
1959 Vol.15      No.1      No.2      No.3      No.4      No.5      No.6
     No.7      No.8      No.9      No.10      No.11      No.12
1958 Vol.14      No.1      No.2      No.3      No.4      No.5      No.6
1957 Vol.13      No.1      No.2      No.3      No.4      No.5      No.6
1956 Vol.12      No.1      No.2      No.3      No.4      No.5      No.6
1955 Vol.11      No.1      No.2      No.3      No.4      No.5      No.6
1954 Vol.10      No.1      No.2      No.3      No.4
1953 Vol.9      No.1      No.2      No.3      No.4
1952
1951 Vol.8      No.1      No.2      No.3
1950 Vol.7      No.5      No.6
1949 Vol.7      No.4
1948 Vol.7      No.3
1947 Vol.7      No.1      No.2
1946 Vol.6      No.2
1945 Vol.6      No.1
1944 Vol.5      No.1      No.2
1943
1942
1941
1940 Vol.4      No.1
1939 Vol.3      No.2
1938
1937 Vol.3      No.1
1936 Vol.2      No.1      No.2
1935 Vol.1      No.3
1934 Vol.1      No.2
1933 Vol.1      No.1
物理学报
· Numerical simulation of soliton trapping of the supercontinuum in photonic crystal fiber [2012, No.12:124203-124203] (38465)
· Large-eddy simulation and experimental study of deflecting oscillation of planar opposed jets [2013, No.8:84704-084704] (38076)
· Effect of concentration of heavy oxygen vacancy in rutile and anatase (TiO2) on electric conductivity performance studied by simulation and calculation [2013, No.23:237101-237101] (30620)
· Quasiparticle band structure calculation for SiC using self-consistent GW method [2012, No.13:137103-137103] (28890)
· Proximity-effect-induced superconductivity by granular Pb film on the surface of Bi2Te3 topological insulator [2013, No.16:167401-167401] (26945)
Copyright © Acta Physica Sinica
Address: Institute of Physics, Chinese Academy of Sciences, P. O. Box 603,Beijing 100190 China
Tel: 010-82649294,82649829,82649863   E-mail: aps8@iphy.ac.cn