Contents Research Articles, 42 Article(s)
Spatial resolution of novel liquid scintillating capillary array
Mei Zhang, Kui-Nian Li, Yang Li, Liang Sheng, and Yan-Hong Zhang

Scintillating array image plates are allowed high resolution through a thicker detector which increases quantum efficiency without degrading the imaging resolution substantially. Due to limitations imposed by process capability, scintillator fiber array with pixel diameter less than 0.2 mm is hardly manufactured to improve performance. Therefore, a liquid scintillator capillary array with 0.1 mm pixel is developed to improve the detection efficiency and spatial resolution of image plate for low intensity radiation imaging. Its performances are studied and tested by simulation and experiment, and are compared with those of scintillating fiber array. Especially in order to gain high fidelity representation of modulation transfer function of the array image plate, a method of simulating and measuring the slanted knife edge response and an iterative algorithm are introduced. For 14 MeV neutron and 1.25 MeV gamma, the slanted knife edge responses of these array image plates with pixel dimensions in a range from 0.1 mm to 0.5 mm are respectively simulated by MCNPx program and the modulation transfer function (MTF) are obtained. The simulation results show that compared with scintillating fiber array, the liquid scintillator capillaries array has an obvious merit in spatial resolution because of greater stopping power for secondary charged particle in the capillary quartz glass wall with 0.02 mm in thickness. Its ultimate resolution can reach to 1.8 lp/mm for 14 MeV neutron by simulation. At the 4000 Ci 60Co facility, a 5-cm-thick tungsten bar, one side of which has a curvature of 0.1 radian to minimize the misalignment effect, is made a knife edge. The MTF of the scintillating fiber array with 0.3 mm and 0.5 mm pixel and newly developed liquid scintillator capillary array is measured through this tungsten knife edge. Experimental measurement results have also verified that the liquid scintillator capillary array performs well in spatial resolution and luminescent uniformity for 1.25 MeV gamma. The ultimate spatial resolution, 0.9 lp/mm is gained, and those of other scintillating fiber arrays are all less than 0.5 lp/mm. Moreover, experimental test validates the simulating method and simulated results, although the measured value is slight less than the simulated value because of the effect of dimension of 60Co source.

Acta Physica Sinica
Jan. 01, 1900, Vol. 69 Issue 6 062801-1 (2020)
Design of a femtosecond electron diffractometer with adjustable gaps
Duan Luo, Dan-Dan Hui, Wen-Long Wen, Li-Li Li, Li-Wei Xin, Zi-Yuan Zhong, Chao Ji, Ping Chen, Kai He, Xing Wang, and Jin-Shou Tian

One of the grand challenges in ultrafast science is real-time visualization of the microscopic structural evolution on atomic time and length scales. A promising pump-probe technique using a femtosecond laser pulse to initiate the ultrafast dynamics and another ultrashort electron pulse to probe the resulting changes has been developed and widely used to study ultrafast structural dynamics in chemical reactions, phase transitions, charge density waves, and even biological functions. In the past three decades, a number of different ultrafast electron guns have been developed to generate ultashort electron sources, mainly including hybrid electron gun with radio-frequency (RF) cavities for compressing the pulse broadening, relativistic electron gun for suppressing the coulomb interaction, single-electron pulses without space charge effect and compact direct current (DC) electron gun for minimizing the electron propagation distance. At present, these developments with different final electron energy and available total charge have improved the time response of ultrafast electron diffraction (UED) setups to a new frontier approaching to 100 fs regime. Although enormous efforts have been made, the superior capabilities and potentials of ultrafast electron diffraction (UED) are still hindered by space-charge induced pulse broadening. Besides, the penetration depth of electrons increases with the electron energy, while the scattering probability of electrons has the opposite consequence. Thus, in addition to the temporal resolution enhancement, it is also important that the electron energy should be tunable in a wide range to meet the requirements for samples with different thickness. Here in this work, we design a novel ultra-compact electron gun which combines a well-designed cathode profile, thereby providing a uniform field and a movable anode configuration to achieve a temporal resolution on the order of 100 fs over an accelerating voltage range from 10 kV to 125 kV. By optimizing the design of the high-voltage electrode profile, the field enhancement factor on the axis and along the cathode surface are both less than ~4% at different cathode-anode spacings, and thus the maximum on-axis field strength of ~10 MV/m is achieved under various accelerating voltages. This effectively suppresses the space charge broadening effect of the electron pulse. Furthermore, the anode aperture is designed as a stepped hole in which the dense sample grid can be placed, and the sample under study is directly supported by the grid and located at the anode, which reduces the cathode-to-sample distance, thus minimizing the electron pulse broadening from the cathode to sample. Moreover, the defocusing effect caused by the anode hole on the electron beam can be effectively reduced, therefore improving the lateral focusing performance of the electron beam.

Acta Physica Sinica
Jan. 01, 1900, Vol. 69 Issue 5 052901-1 (2020)
Theoretical study of angle-cut collimator based design in high-energy proton radiography
Feng Chen, Hai-Bo Xu, Na Zheng, Qing-Gang Jia, Ruo-Gu She, and Xing-E Li

The angle-cut collimator plays an important role in high-energy proton radiography. By using the collimator, the image contrast can be improved, and the material diagnosis and density reconstruction can be realized through secondary imaging. As all these techniques depend on the flux value, it is of great significance to reduce the error of the detected flux value. The ideal collimator is a much thin surface, but thick enough to block protons outside the collimation region. It is designed by stretching the aperture of the collimation plane. The shape is cylindrical, and it will increase the error of the flux value with the angle truncation. The initial bunch is defined and the phase diagram of the bunch within the angle-cut is ideal in the theoretical model. The equation of designing the collimator is given by theoretical analysis. It is given by the transfer matrix, the radius of the object and the angle-cut. The pore structure is oval-shaped by calculating and simulating. The proton imaging system of 1.6 GeV is established by Geant4 program, and the detector is ideal. The round copper plate and the concentric spheres are chosen as objects respectively. The parameters of the designed collimator is given by this method. The ideal collimator, tensile collimator and designed collimator are used in simulation, the radius of object is 5 cm and the angle-cut is 2 mrad and 3.5 mrad. The results show that when using the ideal and the designed angle-cut collimator, the flux distributions are in good agreement, while when using the tensile collimator, the result is quite different from that obtained by using the ideal collimator. Therefore, the collimator designed by this method can effectively reduce the error of the detected flux value.

Acta Physica Sinica
Jan. 01, 1900, Vol. 69 Issue 3 032901-1 (2020)
Experimental study on neutron single event effects of commercial SRAMs based on CSNS
Xun Wang, Feng-Qi Zhang, Wei Chen, Xiao-Qiang Guo, Li-Li Ding, and Yin-Hong Luo

The experiment of neutron single event effect was carried out at China Spallation Neutron Source (CSNS) back-n on 13 kinds of commercial SRAM. The single event upset (SEU) cross section of each device was obtained, and multiple cell upsets (MCU) were extracted from the SEUs using a statistical method without layout information. The influences of the test pattern, feature size and device layout on the SEU cross section and MCU were studied. The results show that the test pattern has little influence on the SEU cross section of the devices, but has a great influence on the MCU ratio of some devices. The feature size has influence both on the SEU cross section and the MCU ratio of the devices. The influence on SEU cross section is not definite. The influence on the MCU ratio is definite. Both the ratio and the maximum size of the MCUs increase with the decrease of the feature size. The difference of layout has great influence both on the SEU cross section and the MCU ratio of the device. In addition, compared with the results of plateau irradiation, the ratio of MCU in CSNS back-n is less than that of plateau irradiation. There are two reasons for this difference. One is that the energy spectrum of CSNS back-n is softer than that of the atmospheric neutron. The other is the neutron beam at CSNS back-n is perpendicular to the device under test. Therefore, evaluating the atmospheric neutron SEE using CSNS back-n line may underestimate the MCU ratio of the device under test. The experimental data, analytical methods and results obtained in this paper are valuable for the researchers to carry out the atmospheric neutron SEE test and the evaluation of devices on atmospheric neutron SEE.

Acta Physica Sinica
Aug. 05, 2020, Vol. 69 Issue 16 162901-1 (2020)
Study of nuclear charge radius
Ying-Yu Cao, and Jian-You Guo

Based on the existing experimental data of nuclear radius, the previous formula of nuclear charge radius is verified and discussed. Comparing the formula of the single-parameter nuclear charge radius, it is proved that the formula of $Z^{1/3}$ law is better than the formula of $A^{1/3}$ law. We refitted the two-parameter formula and the three-parameter formula that have been proposed and confirmed that the two-parameter and three-parameter formula fit better than the single-parameter formula. It is shown that show that the deformation plays a key role in the nuclear charge radius. The electric quadrupole moment is an important physical quantity representing the properties of the nucleus. Its appearance indicates the deviation from spherical symmetry and also reflects the size of the nuclear deformation. The electric quadrupole moment is also one of the basic observations to understand the distribution of matter within the nucleus, to examine the nuclear model, and to observe nucleon-nuclear interactions. Taking into account the relationship between the nuclear quadrupole moment and the deformation, the electric quadrupole moment factor is added to the original three-parameter formula to obtain a new formula for the nuclear charge radius. Fitting the four-parameter formula, it is found that the theoretical value of the nuclear charge radius is in good agreement with the experimental value, the root-mean-square deviation is 0.0397 fm. Considering the relationship between the total spin and the electric quadrupole moment, the intrinsic electric quadrupole moment is obtained and brought into the formula for fitting, and the root-mean-square deviation further decreases,the root-mean-square deviation is 0.0372 fm. Finally, considering the universality of odd-even staggering, we add the $\delta$ term that can reflect the odd and even oscillation phenomenon, and the root-mean-square deviation obtained by the formula is 0.369 fm, which better reflects the relationship between the deformation and the nuclear charge radius. Compared with the formulas already proposed, the new formula can better reflect the variation trend of nuclear deformation, shell effect, odd-even staggering, etc., and the calculation accuracy is also improved, which can provide a useful reference for future experiments.

Acta Physica Sinica
Aug. 05, 2020, Vol. 69 Issue 16 162101-1 (2020)
Monte Carlo simulation technology based on characteristic γ-ray spectrum analysis
Li Deng, Rui Li, Xin Wang, and Yuan-Guang Fu

Monte Carlo method is an ideal way to simulate criticality, shielding and nuclear detection. JMCT is a multipurpose 3D Mont Carlo (MC) neutron-photon-electron and coupled neutron /photon /electron transport code which is developed by IAPCM. The program is developed based on the combinatorial geometry parallel infrastructure JCOGIN and has the most functions of general Monte Carlo particle transport code, including the various variance reduction techniques. In addition, some new algorithms, such as Doppler broadening on-the-fly (OTF), uniform tally density (UTD), consistent adjoint driven importance sampling (CADIS), fast criticality search of boron concentration (FCSBC), the domain decomposition (DD), the two-level parallel computation of MPI and OpenMP, etc. have been developed, where the number of geometry zones, materials, tallies, depletion zones, memories and period of random number are big enough to simulate various extremely complicated problems. Also the JMCT is hybrid the discrete ordinate SN program JSNT to generate source biasing factors and weight window parameters for deep-penetration shielding problems. The input is based on the CAD modeling, and the result is a visualized output. The JMCT can provide technology support for radiation shielding design, reactor physics and criticality safe analysis. Especially, the JMCT is coupled depletion and thermal-hydraulic code for simulating the reactor feedback effect, including depletion, thermal feedback. In recent years, new function of γ-ray spectrum analysis has been developed.In this paper, the working principles of timing measure are introduced. The advanced calibration count is developed for distinguishing between inelastic γ-ray and capture γ-ray based on time bin tally. On the other hand, when neutron collides with nuclide, the secondary photon is labeled into the primary line photon and primary continuous photon, where energy of primary line photon does not change with the incident neutron energy, such as carbon spectral-line at 4.43 MeV and oxygen spectral-line at 6.13 MeV. The element components of detected object can be determined by the primary line photon. On the other hand, expect value estimator (EVE) is used to produce the secondary photons. The advantage of EVE does not leak any event even with a small probability which is important for detecting the hide exploder. However the shortage of the EVE results in producing a great number of photons with small weight. If all of these small weight photons are simulated one by one, a great amount of computation time and memory will be consumed. For avoiding this case, a new algorithm is design by coupling EVE and DE (direct estimator). The all of secondary photons from EVE only make the direct tally take a little computing time, then end the photon history and return to the DE production photon model (one photon production at most). Final, the total tally is a summation of EVE direct tally and DE scattering tally. The use of new algorithm to realize the analysis of γ-ray spectrum will increase only a little computing time. The numerical tests are done by using own Monte Carlo code JMCT. The correctness and validity of the algorithm are shown preliminarily.

Acta Physica Sinica
Jan. 01, 1900, Vol. 69 Issue 11 112801-1 (2020)
Preparation and properties for X-ray scintillation screen based on ZnO:In nanorod arrays
Qian-Li Li, Ya-Hua Hu, Yi-Fan Ma, Zhi-Xiang Sun, Min Wang, Xiao-Lin Liu, Jing-Tai Zhao, and Zhi-Jun Zhang

X-ray scintillation screens as the core component of X-ray imaging detectors have widespread applications in the medical imaging, security inspection, high energy physics, radiochemistry, and so on. For a long time, the development of X-ray scintillation screen mainly focuses on improving the light yield in order to enhance its detection efficiency. However, a novel tendency has recently emerged for ultrafast time performance of the X-ray imaging detector. The indium doping zinc oxide (ZnO:In) with high radiation hardness, higher light yield(>10000 photons/MeV) and subnanosecond decay time is a promising scintillation material for ultrafast detections. In order to satisfy the requirements of X-ray scintillation screens with ultrafast and high-spatial-resolution in the existing and upcoming high energy physics experiments, the ZnO:In nanorod arrays have been prepared on a 100-nm-thick ZnO-seeded substrate by hydrothermal reaction method and then treated by hydrogen plasma in present work. The results of SEM demonstrate the average diameter and length of the ZnO:In nanorods are about 0.5 and 12 μm, respectively. The XRD shows the ZnO:In nanorods are highly aligned perpendicular to the substrate alongc-axis direction. The X-ray excited luminescence spectra show that two luminescence bands are observed, i.e. an ultraviolet emission peak located at about 395 nm and a visible emission band at 450–750 nm. It is particularly important to point out that hydrogen plasma treatment can enhance the ultraviolet emission of ZnO:In nanorod arrays and suppress its visible emission. The reason is attributed to the formation of shallow donors through hydrogen entering the ZnO and the combination of VO and Oi. In addition, the fluorescence decay times of the ultraviolet and visible emissions for the ZnO:In nanorod arrays are subnanosecond and nanosecond, respectively, satisfying the demand of the fast X-ray imaging. The spatial resolution of ZnO:In nanorod arrays has been characterized in X-ray imaging beamline at the Shanghai Synchrotron Radiation Facility. Under excitation of the X-ray beam with the energy of 20 keV, a system spatial resolution of 1.5 μm could be achieved by using an 12 μm thickness ZnO:In nanorod arrays as the scintillation screen, which is exceeded the highest level had ever been reported on ZnO:In nanorod arrays scintillation screen. In conclusion, this present work shows that it is a feasible solution for X-ray detection and imaging with high temporal and spatial resolution by using ZnO:In nanorod arrays as the X-ray scintillation screen.

Acta Physica Sinica
Jan. 01, 1900, Vol. 69 Issue 10 102902-1 (2020)
Investigation of n-α scattering by combining complex momentum representation and Green’s function
Xiao-Wei Wang, and Jian-You Guo

Nuclear scattering is a very important physical phenomenon in which the resonance state plays an important role. In order to study the two-body system n-α scattering, Green’s function is introduced under the complex momentum representation, so the complex momentum representation-Green’s function approach is established. This method is used to study the elastic scattering of n-α system. By extracting the resonances, it is found that the contributions of resonances in continuum level density, phase shift, and cross section are more important. In the case without introducing any non-physical parameters, it is very helpful to understand the resonant states and the non-resonance continuum states by analyzing the data of scattering states. In this work, we mainly study the p-wave scattering with the orbital angular momentum l = 1, where P1/2 is a wide resonance state and P3/2 is narrow resonance state. The study shows that the sharp resonance peak of p-wave scattering gives rather broad distribution to the scattering phase shift and the cross section of the n-α system. By comparison, we can see that the theoretical calculation results and experimental data are in good consistence.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 9 092101-1 (2019)
Influence of mass sample thickness on quantitative analysis of neutron activation and its modification
Yuan-Ming Sun, Xu Xu, Wan-Yue Tang, Yi Chang, Jing-Bin Lu, Long Zhao, and Yu-Min Liu

With traditional neutron activation analysis, the increase of sample quality leads to some problems in both irradiation process and measurement process. These problems include the neutron flux gradient in the range of the sample, the decrease of the thermal neutron flux rate around the sample and the influence of self-shielding in the sample in the irradiation process, In the process of measurement, the self-attenuation of γ-ray in the sample and the geometric effect of the sample lead to the effect of the detector on the measurement of characteristic γ-ray emissivity due to the difference in the detection efficiency of each point of the sample. So the neutron activation analysis of mass sample needs to make some additional modifications. By using the neutron activation analysis technique, the content of 24Mg and 28Si in a large amount of flour can be detected, and the content of talc powder in the flour can be given, so that the quality of flour can be monitored. The flour mainly contains C, H, O, N, Ca, and F element, but the main chemical constituent of talc powder is Mg3[Si4010](OH)2. Therefore, the measured content of Si and Mg element can be used to judge whether the flour contains talcum powder and to determine its exact content. When the content of 24Mg and 28Si in flour are measured by the neutron activation analysis, the variation of neutron flux and energy with thickness in the measured sample and the effect of γ-ray self-absorption will have great influence on the measurement results. The relationship between the neutron flux and energy and the thickness of the sample is simulated by MCNP5 (Monte Carlo N-particle transport code system 5), and the neutron fluxes at different thickness of the sample are measured by a 3He proportional counter tube. The results show that the simulation results of MCNP5 are basically consistent with the experimental results. Using the simulation by MCNP5 and the measurements by a sodium iodide detector, the relationship between γ-ray self-absorption effect and sample thickness is studied, and the sample thickness is determined to be that is adopted as an optimal experimental condition. Based on the simulated data, the function relationship between the counting of 1.779 MeV γ-ray and the thickness of the sample is obtained as follows: A = 1401 + 3815x – 720x2 + 64x3 – 2.8x4 + 0.05x5. The curve trend of the experimental results is basically the same as that of the simulation results.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 8 082801-1 (2019)
Exploration of the deep-ultraviolet nonlinear optical materials in the derivatives of KBe2BO3F2
Min-Qiang Gai, Ying Wang, and Shi-Lie Pan

The use of nonlinear optical crystal materials to extend the limited range of laser sources to the deep-ultraviolet (deep-UV, λ < 200 nm) regions by various frequency conversion techniques, has become an attractive field for generating deep-UV light. However, the lack of nonlinear optics (NLO) crystal materials capable of frequency conversion in the deep-UV light range, limits the development and application of deep-UV all-solid-state lasers. Therefore, scientists all over the world are actively exploring the new generation of deep-UV NLO crystal materials. At present, only the KBe 2BO3F2 (KBBF) crystal is capable of generating deep-UV light through the direct sixth harmonic generation of the Nd:YAG laser. The infinite ∞[Be2BO3F2] single layers, as the brilliant building blocks in the crystal structures of KBBF family, provide a relatively large second harmonic generation coefficient (d11 = 0.47 pm/V) and a sufficient birefringence (Δn = 0.07@1064 nm). However, the KBBF crystals have insurmountable intrinsic defects, such as the usage of high toxic beryllium oxide, and the serious layer growth habit, which greatly restrict its commercialization process. Since the layered structure of the KBBF crystal is still one of the most brilliant structures for generating deep-UV laser, an effective strategy is to change the interlayer connection mode and develop new NLO materials based on KBBF with less layering growth habit. In this paper, by reviewing the development history of borate deep-UV NLO crystals and the derivatives of KBBF, the relationship between layered structure and optical properties of different interlaminar connections of crystal materials is systematically analyzed. We discuss the main contradictions and solutions of the development of deep-UV NLO crystal materials which are similar to the KBBF structure. In order to provide a reference for the innovative exploration of new materials in the future, several design strategies are also proposed.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024208-1 (2019)
Confinement effect of carbon nanotubes on the chain mobility of conjugated polymer poly(9,9-dioctylfluorenyl-2,7-diyl)
Ling-Dong Li, An-Na Ye, Sheng-Lin Zhou, Xiao-Hua Zhang, and Zhao-Hui Yang

The conjugated polymer polyflourene has been well studied for its strong blue light emission ability and high quantum efficiency behavior. It has wide applications for light emitting diodes, sensors as well as photo-detectors. Therein the \begin{document}$ \beta $\end{document} conformation of PFO crystals is more attractive due to its longer conjugation length, higher carrier mobility and better luminous efficiency. Therefore it is great essential to control the formation and stability of \begin{document}$ \beta $\end{document} conformation of PFO crystals to develop new kind of photo-electronic devices. As is known, polymeric materials confined in a nanometer-sized space often exhibit unique properties compared with their bulk state, such as abnormal chain mobility, molecular assembly and phase transition behavior. These factors are of great significance to develop new kind of material and applications. Generally the confined condition includes quantum dot (zero-dimensional, 0D), nanowire or nanotube (1D), ultrathin film (2D) and nanoparticle (3D). In this paper, we design a unique 1D nanoconfined environment based on vertically aligned carbon nanotube (CNT) array structure. An ultra-high CNT density is achieved through a solvent-induced contraction process. The adjacent narrow carbon nanotube gap thus forms a quasi-1 confined nano-space with the tunable size ranging from 5 to 50 nm. Then we infiltrate the conjugated polymer poly(9,9-dioctylfluorene-2,7-diyl) (PFO) into those nano-gaps of carbon nanotube arrays through a solvent evaporation method to obtain the PFO infilled CNT array composite film. It is found that the chain mobility of PFO molecules in such a 1D nano-confined space of carbon nanotubes is significantly suppressed compared with the scenario of the spin-coated PFO film. The transition speed between different crystal forms of PFO declines greatly, which meanwhile improves the thermal stability of the \begin{document}$ \beta $\end{document} conformation of PFO crystal. Additionally, the aligned carbon nanotubes have great effects on the orientation and distribution of PFO chains. The PFO crystals are confirmed to grow preferentially along the longitudinal direction of CNT array, which is potential to grow PFO crystals with high quality and excellent performance. Thus, such a PFO/CNT array composite film can have great potential to prepare PFO photovoltaic devices with excellent luminescent properties and high stability in the future.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 026402-1 (2019)
Optimal design of large-sized sandwich transducer based on two-dimensional phononic crystal
Sha Wang, and Shu-Yu Lin

Sandwich transducers are extremely versatile, but when the lateral dimension is too large, the displacement of the radiating surface is uneven due to the coupling vibration. Due to the unique vibrational band gap characteristics of phononic crystal, vibrations in the bandgap range can be prohibited from propagating for infinite periodic structure or suppressed for finite periodic structure, which makes it widely used in the field of vibration suppression. In this paper, a two-dimensional phononic crystal structure is formed by processing periodically aligned grooves in the front cover of a large-sized sandwich transducer. Since the periodic grooves are formed in the radial direction, the radial waves cannot propagate, and thus the lateral vibration is suppressed. Subsequently, the finite element method is used to simulate the vibration transmission characteristic, resonance frequency and emission voltage response of a large-sized sandwich transducer based on two-dimensional phononic crystal. The effects of slot height and slot width on its bandgap, resonance and anti-resonant frequency, bandwidth, and displacement profile of the radiating surface are discussed. The results show that the phonon crystal structure can be optimized by using a large-sized sandwich transducer. The large-sized sandwich transducer based on two-dimensional phononic crystal has a lateral band gap. When the operating frequency of the large-sized sandwich transducer is within the band gap range, the two-dimensional phononic crystal structure can effectively suppress the lateral vibration, and the uniformity of the displacement distribution of the radiating surface of the transducer is improved. In addition, when the slot width is constant, the bandwidth of the large-sized sandwich transducer based on the two-dimensional phononic crystal increases as the slot height increases. Similarly, when the slot height is constant, the bandwidth of the large-sized sandwich transducer based on the two-dimensional phononic crystal increases as the slot width increases. The two-dimensional phononic crystal structure is processed in the front cover of the large-sized sandwich transducer, and the working frequency band of the large-sized sandwich transducer is effectively broadened.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024303-1 (2019)
Influence of plasma shielding effect on ground state and excited state energies of Ar16+
Kun Ma, Zhan-Bin Chen, and Shi-Zhong Huang

A systematical knowledge of the atomic properties in plasma is of great interest for various research areas, such as the explanation of the X-ray radiation from universe, plasma diagnostics, extreme ultraviolet (EUV) and X-ray sources and so on. Among these researches, the detailed information about how the plasma influences the atomic energy level and transition spectrum are crucial for understanding the X-ray emission mechanism and the state of plasma. An analytic calculation method of treating the non-relativistic energy and its relativistic corrections for the multi-electron atoms embedded in weakly coupled plasma is developed based on the Rayleigh-Ritz variation method. The systematical investigations are performed for the ground state 1s21S, single excited states 1sns 1,2S (n = 2 5), 1snp 1,3P (n = 2 5) and double excited state 2s2p1P of Ar16+ ion in weak coupled plasma. The analytic formulas for calculating the non-relativistic energy and its relativistic correction energy are derived, which include mass correction, one and two-body Darwin correction, spin-spin contact interaction and orbit-orbit interaction. All the angular integration spin sums involved in the problem are worked out explicitly by using the irreducible theory. The influence of plasma on non-relativistic energy and relativistic correction energy are discussed. The results show that the mass correction and the one-body Darwin correction are the main ones among the terms of relativistic correction, and are three orders of magnitude greater than the other relativistic terms. The plasma shielding effect mainly affects the non-relativistic energy, and has little effect on the relativistic correction. At the same time, it has a more significant selectivity for the electronic configuration. Further research shows that the influence of plasma on the energy of the outer shell electron is greater than that of the inner shell electron. With the increase of the plasma shielding parameters, the outer shell electron extends outward, and the higher the excited state, the greater the degree of extension is. This work should be useful for astrophysical applications where such a plasma environment exists.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 023102-1 (2019)
Molecular simulation of interaction between charged nanoparticles and phase-separated biomembranes containning charged lipids
Yi-Ran Liang, and Qing Liang

Nanoparticles have been widely used in many fields such as nanomedicine and cell imaging. Understanding the microscopic mechanism of the interaction between nanoparticles and biomembranes is very vital for the synthesis and applications of nanoparticles. In this paper, using coarse-grained molecular dynamics simulation, we study the interaction between nanoparticles coated with fully or partially charged ligands and phase-separated biomembranes containing charged lipids. The results show that the final positions or states of nanoparticles on/in the biomembranes can be readily modulated by varying the grafting density, ratio, and type of charged ligands as well as the type of charged lipids. For the nanoparticle with a highly hydrophilic surface, the nanoparticle prefers to be adsorbed on the surface of the biomembrane. In this case, the electrostatic interaction determines that the nanoparticle is adsorbed on the surface of liquid-ordered domain or the surface of liquid-disordered domain. For the nanoparticle with a (partially) hydrophobic surface, the nanoparticle tends to penetrate into the lipid bilayer from the liquid-disordered domain. In this case, the hydrophobicity of the nanoparticle plays a crucial role in the penetrating of the nanoparticle. The hydrophilicity or hydrophobicity of the nanoparticle is affected by the ratio between the charged and neutral ligands, the grafting density of the charged ligands, and the ionic concentration in the system. Furthermore, the microscopic mechanism of the interaction between charged nanoparticles and charged biomembranes is revealed by using the potential of mean force between nanoparticles and lipid domains. The potential of mean force shows that none of the (partially) charged nanoparticles can spontaneously penetrate into the liquid-ordered domain due to a high free energy barrier but they can spontaneously penetrate into the liquid-disordered domain with a certain probability. However, due to the limitation of the simulation time and the number of sampling of the simulations, only some of the partially hydrophobic nanoparticles which are not initially adsorbed onto the surface of liquid-ordered domain are found to finally penetrate into the liquid-disordered domain in this work. This work yields some theoretical insights into the application of nanoparticles in nanomedicine, cell imaging, etc.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 028701-1 (2019)
Ultrafast spin dynamics in double-magnetic-center endohedral fullerene Y2C2@C82-C2(1)
Rui Huang, Chun Li, Wei Jin, Lefkidis Georgios, and Hübner Wolfgang

Spin switching and spin transfer are essential prerequisites for designing the spin-logic devices based on endohedral fullerenes. In this paper by combining the theoreticalΛ-process model with a self-designed genetic algorithm, we are able to theoretically observe spin-switching and spin-transfer scenarios on the subpicosecond time scale in the endohedral fullerene Y2C2@C82-C2(1) from first principles. The results show that the geometry of the optimized enclosed Y2C2 cluster is consistent with the experimental data. There exists a certain repulsive force between the external C82-C2(1) cage and the encaged cluster. However, the whole system still maintains its integral cage structure due to the excellent stability of the fullerene. In the Y2C2@C82-C2(1) system, it is found that the spin density is highly localized on the two Y atoms and only minimally distributed on the carbon cage. By analyzing the spin-density distribution and the evolution of the spin expectation values as influenced by the laser pulses, it is found that global spin switching can be achieved on the two Y atoms, while spin transfer between the two Y atoms actually results from the redistribution of the spin density among the two magnetic centers and the carbon cage under the action of the optimized laser pulses. The achieved spin-switching scenario completes within about 1000 fs and its fidelity reaches 97.8%, while the obtained spin-transfer process completes within 200 fs and its fidelity reaches 95.1%. The electron absorption spectra of the system verify that optical transitions are possible between the main intermediate states and the initial and final states involved in the spin-switching and spin-transfer scenarios. Therefore, by analyzing the electron absorption spectra corresponding to the initial and final states, the energy of the laser pulses adopted for the studied spin-dynamics process can be predicted, and the spin transferability can be evaluated. In addition, it is found that the smaller the detuning between the required energy difference and the applied laser pulse energy is, the greater the probability for spin switching/transfer scenarios becomes. The present results reveal the mechanisms of the laser-induced ultrafast spin dynamics in Y2C2@C82-C2(1) and can provide a theoretical basis for designing the spin-logic devices on realistic endohedral fullerenes.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 023101-1 (2019)
Electrons leakage from upper laser level to high energy levels in active regions of terahertz quantum cascade lasers
Jin-Feng Li, Ting Wan, Teng-Fei Wang, Wen-Hui Zhou, Jie Xin, and Chang-Shui Chen

Terahertz quantum cascade laser is a semiconductor laser that effectively obtains terahertz waves. It uses the semiconductor heterojunction to have a quantum cascade effect under an applied voltage, and then the phonon assists the electron resonance from the upper stage to the next stage, so that a single electron injected externally can emit multiple photons. However, some electrons will deviate from the transport path during transportation and these electrons are called leakage electrons. Electron leakage comes from three ways. The first way is the scattering of electrons from the upper laser level through the long longitudinal phonon to the low energy level; the second way is the scattering of electrons from the lower laser level to the high energy bound level and the continuous level; and the third way is the scattering of electrons from the upper laser level to high energy bound levels and continuous levels. These leakage electrons directly reduce the number of population inversions in the laser system, making the laser output power limited. At present, most of researchers explain the electron leakage through indirect measurements, and there are few studies in which the electron leakage is analyzed by establishing theoretical models. In this paper, the electron leakage model in THz QCL is established by using thermodynamic statistical theory and laser output characteristic theory. The degree of electron leakage is measured by output power. The influence of lattice temperature and quantum well barrier height on electron leakage are studied. It is found that when the lattice temperature rises and the electrons in the upper laser state leak to higher energy levels, the number of electrons leaking to the adjacent bound state and the continuous state increases, and the number of electrons leaking to the next near-bound level is relatively small. In the case of electron leakage, the utilization of electrons becomes lowered, and the laser output power is also lowered. The study also shows that an appropriate increase in the height of the quantum barrier can suppress the leakage of electrons. Using the established theoretical model to optimize the quantum well barrier height of the previously reported laser system, an 8 mW terahertz quantum cascade laser (THz QCL) laser output at 210 K is obtained. Compared with the reported experimental results, the temperature and output power are improved. These results provide a theoretical basis for studying the electron leakage temperature characteristics of THz QCL and also optimally designing the THz QCL active region structure.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 021101-1 (2019)
Size effect on boundary condition at solid-liquid interface in microchannel
Long-Yan Zhang, Jin-Liang Xu, and Jun-Peng Lei

The heat transfer in microchannel has attracted considerable attention due to many important applications in biology, chemistry, physics and engineering. When the fluid size shrinks to nanoscale, the energy transport of micro-system is significantly different from the conventional case. It is of great significance to study the size effect on heat transfer in a micro-system. However, there is a large size gap between existing molecular dynamics simulation and experimental measurement, in which the size effect on solid-liquid interfacial thermal resistance is rarely involved. Non-equilibrium molecular dynamics simulation is performed to investigate the heat transfer through the solid-liquid interface. Simple Lennard-Jones (LJ) fluid is simulated as the ultra-thin liquid film in a non-equilibrium simulation system. The liquid film is confined in a nanochannel composed of two solid surfaces. The potential function between solid and liquid atom is represented by a modified LJ function to control the solid-liquid interfaces of different surface wettabilities. We examine the size effect on temperature jump and thermal resistance at the solid-liquid interface. The fluid number density and temperature distribution in the perpendicular direction of solid wall are evaluated. It is found that the liquid atoms near wall are arranged as a solid-like structure. Particularly in the small channel, liquid atoms confined in the channel are affected by two solid walls. However, with the increase of channel height, the liquid atoms in the middle channel move freely, leading to the decrease of the size effect. The simulation results show that the dependence of thermal resistance on microchannel height exhibits two regimes: (i) monotonically increasing dependence for the small channel and (ii) keeping constant thermal resistance for the large channel. These two distinct trends can be explained by phonon vibrational density of states (VDOS) of solid wall and liquid. For the small channel, a stronger confinement of liquid leads to a weaker mismatch in VDOS of solid wall and liquid, thus resulting in a smaller thermal resistance. Whereas, for the large channel, the vibrational coupling between the solid and the liquid atom remains unchanged and the size effect is negligible. The size thresholds of the two regimes of the thermal resistance are both sensitive to the liquid-solid interaction strength, which decreases with solid-liquid interaction increasing. Furthermore, with the increase of the microchannel height, the temperature jump at the solid-liquid interface monotonically decreases and eventually approaches to the non-jump temperature boundary on a macroscopic scale. These findings may help to understand the mechanism of temperature boundary conditions on a microscopic scale and a macroscopic scale and provide a theoretical support for manufacturing new nano-devices.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 020201-1 (2019)
Effects of seed crystal size on growth of gem-diamond single crystal
Yu-Kun Qin, Hong-Yu Xiao, Li-Na Liu, Rui-Rui Sun, Qiu-Bo Hu, Zhi-Gang Bao, Yong-Sheng Zhang, Shang-Sheng Li, and Xiao-Peng Jia

In the paper, under 5.6 GPa and 1250 1450 ℃, the Ib-ype diamond single crystals chosen as the seed crystals with different sizes, are synthesized in a cubic anvil at high pressure and high temperature. High-purity Fe-Ni-Co solvents are chosen as the catalysts. High-purity graphite powder (99.99%, purity) is selected as the carbon source. Hexahedral abrasive grade high-quality diamonds of 0.8 mm, 1.5 mm or 2.2 mm in diameter are chosen as seed crystals. The effects of seed crystal size on the growth of gem-diamond single crystal are studied in detail. Firstly, the influence of the change of seed size on the cracking of diamond single crystal is investigated. The crystal growth law of increasing the probability of cracking crystal with larger seed crystal is obtained. It can be attributed to the following two points. i) The residual cross section at the separation of the main crystal from the larger seed crystal is too large, thus reducing the overall compressive strength of the crystal. ii) The growth rate of the diamond crystal synthesized by larger seed crystal is too fast, which leads to the increase of impurities and defects and the decrease of compressive strength of the crystal. The decrease of crystal compressive strength leads to cracks in diamond crystals during cooling and depressurizing. Secondly, in the growth time of 25 hours, the relationships between the growth time and the limit growth rate of the diamond single crystals synthesized by choosing three sizes of seed crystals are investigated. The results show that the high-quality single crystal synthesis efficiency can be improved and the synthesis period can be shortened by selecting large seed crystals. This is because the size of the seed crystal becomes larger at each stage of crystal growth, resulting in the enhancement of the ability of diamond single crystal to receive carbon, so that high-quality diamond single crystals can be grown at a faster growth rate. Thirdly, with the help of scanning electron microscope or optical microscope, we calibrate the surface morphologies of diamond single crystals grown with different-size seed crystals. Using the seed crystals of 0.8 mm, 1.5 mm or 2.2 mm in diameter, high-quality diamond single crystals with smooth surfaces can be synthesized. However, with the increase of seed crystal in size, the surface flatness of the grown crystals tends to decrease and the possibility with which surface defects occur and string inclusions increase. The growth rate of high-quality diamond single crystals grown with larger seed crystals must be strictly controlled. Finally, the N impurity content values of diamond single crystals grown with different seed crystals in size are characterized by Fourier transform infrared measurement. The results show that the N impurity content of the crystal increases with the diamond growing rapidly by selecting larger seed crystal.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 020701-1 (2019)
A method of efficiently generating arbitrary vector beams
Shu-Xia Qi, Sheng Liu, Peng Li, Lei Han, Hua-Chao Cheng, Dong-Jing Wu, and Jian-Lin Zhao

Vector beams have been used in scientific and engineering researches due to their unique focusing properties. In recent years, many methods of generating the vector beams have been proposed, among which the spatial light modulator (SLM) is widely used based on the superposition principle with using orthogonally polarized beams. However, the energy waste is generally associated with these superposition methods. How to efficiently generate vector beams is still a hot topic. Recently, we proposed an efficient method to generate tunable vector beams by using two triangular common-path interferometers (TCPIs) as the beam splitting and combining system. However, due to the complex structure of the TCPI, the system is difficult to adjust and unstable. In addition, the optical system brings about a long optical path, and the vector beams consisting of non-eigen modes will be distorted obviously with a long distance propagation. In this paper, an improved method is proposed. We replace the TCPIs with a pair of beam displacers, which act as a beam splitter and combiner, respectively. In this setup, we can arbitrarily manipulate the polarization states and phase distributions of vector beams in real time by managing the phase diagrams load on the SLM. The whole optical system does not involve any diffractive optical elements, and has a higher conversion efficiency. The improved optical system is compact and stable, and makes the adjustment of coaxiality easier. The light energy utilization depends mainly on the reflectivity of SLM. The efficiency of generating vector beams is increased to 58% by using an SLM with a reflectivity value of 79%. Several typical vector beams with phases and tunable amplitude, including cylindrical vector beams, fractional vector beams, and vector beams with double singularities, double-mode, radially variant polarization distribution, and azimuthally and radially variant polarization distribution, are generated and verified well experimentally. This method is also expected to create high-power vector beams and play an important role in laser processing and light trapping.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024201-1 (2019)
Influence of surface duct on the vertical spatial characteristics of wind-generated noise in deep ocean
Guang-Yu Jiang, Chao Sun, Lei Xie, and Xiong-Hou Liu

Vertical spatial characteristics of the wind-generated noise include the noise vertical directionality and the noise vertical coherence, which seriously affect the performance of sonar devices and play an important role in ocean environment parameters inversion. In this paper, we investigate the influence of surface duct on the noise vertical spatial characteristics at the depths below the duct. The Kuperman-Ingenito (K/I) model is employed to describe the distribution of the noise sources, and Pekeris-branch-cut-based normal modes are used to represent the Green's functions between the noise sources and the receivers. Both the noise vertical directionality and the noise vertical coherence are expressed as a function of the normal modes, so that we can investigate the physical reason for the variance of the noise vertical spatial characteristics by analyzing the variance of the normal modes. The numerical simulations on the noise vertical spatial characteristics show that the influence of surface duct above the critical depth is different from that below the critical depth. Above the critical depth, there exists a peak in the noise vertical directionality at the edge of the horizontal notch close to the bottom side. In the presence of surface duct, this peak significantly rises up, and the noise vertical coherence deviates from that in the absence of surface duct and tends to be perfect positive coherence and perfect negative coherence periodically as the vertical distance between the two receivers increases. On the other hand, below the critical depth, the noise power comes from the horizontal direction becomes stronger and the noise vertical coherence tends to be perfect positive coherence in the presence of surface duct as compared with the case in the absence of surface duct. Moreover, the influence will become severer as the thickness of the surface duct increases, while almost keep unchanged when the sound speed gradient in the surface duct varies. The modal analysis indicates that the noise source excites more refracted normal modes with stronger modal intensity in the presence of surface duct, and the excited refracted normal modes become more and stronger if the surface duct is thicker. As the result, the increase of the refracted mode number and the enhancement of their modal intensity cause the vertical spatial characteristics of noise to change.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024302-1 (2019)
Recognition of denatured biological tissue based on variational mode decomposition and multi-scale permutation entropy
Bei Liu, Wei-Peng Hu, Xiao Zou, Ya-Jun Ding, and Sheng-You Qian

It is an important practical problem to accurately recognize whether biological tissue is denatured during high intensity focused ultrasound (HIFU) treatment. Ultrasonic scattering echo signals are related to some physical properties of biological tissues. According to the characteristics of ultrasonic scattering echo signals, the recognition of denatured biological tissues is studied based on the variational mode decomposition (VMD) and multi-scale permutation entropy (MPE) in this paper. The ultrasonic echo signals are decomposed into various modal components by the VMD. The noise components and the useful components are separated according to the power spectrum information entropy of various modal components. The separated useful signals are reconstructed and the MPE are extracted. Furthermore, Gustafson-Kessel (GK) fuzzy clustering analysis is employed to obtain the standard clustering center, and the recognition of denatured biological tissues is carried out by Euclid approach degree and principle of proximity. The proposed method is applied to ultrasonic scattering echo signal during HIFU treatment. In order to determine the parameters of MPE algorithm for ultrasonic scattering echo signals, the embedding dimension of the MPE is discussed, and the scale factor of the MPE algorithm is optimized by genetic algorithm. When the delay time and the embedding dimension are 2 and 7 respectively, the MPE values decrease with scale factor increasing. Assuming that the scale factor is 12 from optimization results, the 293 ultrasonic scattering echo signals from normal tissues and denatured tissues are analyzed by the MPE. It is found that the MPE values of the denatured tissues are higher than those of the normal tissues. The MPE can be used to distinguish normal tissues and denatured tissues. Comparing with the recognition methods of the EMD-MPE-GK fuzzy clustering method and the VMD-WE-GK fuzzy clustering, the proposed method has good clustering performance and separability. Its partition coefficient (PC) is close to 1 and the Xie-Beni (XB) index is smaller. There are fewer feature points in the overlap region between MPE features of denatured tissues and normal tissues. The recognition results of denatured biological tissues in this experimental environment show that the recognition rate based on this method is higher, reaching up to 93.81%.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 028702-1 (2019)
Passivation optimization and performance improvement of planar a-Si:H/c-Si heterojunction cells in perovskite/silicon tandem solar cells
Jun-Fan Chen, Hui-Zhi Ren, Fu-Hua Hou, Zhong-Xin Zhou, Qian-Shang Ren, De-Kun Zhang, Chang-Chun Wei, Xiao-Dan Zhang, Guo-Fu Hou, and Ying Zhao

Recently, the monolithic spin-coating perovskite/planar silicon heterojunction tandem solar cells with high performance have attracted attention mainly due to simple fabrication, low preparation cost and high efficiency, especially compared with fully textured multi-junction perovskite/silicon tandem device. As is well known, the excellent passivation of a-Si:H/c-Si interface is the key to achieving a high-efficiency planar silicon heterojunction solar cell, which further improves the performance of the corresponding tandem cell. Therefore, we investigate the elements affecting a-Si:H/c-Si interface passivation, including the c-Si surface treatment technique, a-Si:H passivation layer and P-type emitter layer and so on. In these experiments, we adjust the immersed time of diluent hydrofluoric acid and pre-deposited hydrogen plasma with different gas mixture flows. Also, the suitable deposition parameters of intrinsic a-Si:H passivation layer are regulated by varying hydrogen dilution and time, and variously slight silane content is embedded into i-a-Si:H /P-type (I/P) emitter interface by hydrogen-rich plasma treating which is for acquiring optimal experimental processing conditions to promote the chemical passivation. In addition, the p-a-Si:H and p-nc-Si:H are comparatively studied as buffer layers to further improve the I/P interface passivation by varying the hydrogen dilution in the gas mixture during deposition. It can be found that p-nc-Si:H buffer layer with high conductivity and wide bandgap can not only reduce the defect density at the I/P interface, but also increase the conductivity of P-type emitter, which further improves the field passivation effect. By the above- mentioned optimization, the highest minority carrier lifetime and implied open-circuit voltage (iVoc) of the structure of P-type emitter/a-Si:H(i)/c-Si/a-Si:H(i)/N-type layer (inip) sample can respectively reach 2855 μs and 709 mV, which demonstrates authentically outstanding passivation performance. An efficiency of 18.76% can be obtained for the planar a-Si/c-Si heterojunction solar cell with a Voc of 681.5 mV, which is 34.3 mV higher than that of the reference device. Regarding the optimized planar a-Si:H/c-Si heterojunction solar cell as the bottom cell, we also obtain an efficiency of 21.24% for perovskite/silicon heterojunction tandem solar cell with an open-circuit voltage of 1780 mV, which proves that the above strategies are very effective for improving the passivation optimization and performance of bottom cell in the tandem device.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 028101-1 (2019)
Non-collinear magnetism and electronic transport of boron or nitrogen doped zigzag graphene nanoribbon
Jin-Tao Liang, Xiao-Hong Yan, Ying Zhang, and Yang Xiao

Zigzag graphene nanoribbon (ZGNR) is important for novel carbon-based spintronic applications. Currently, most of ZGNR spintronic studies focus on the collinear magnetism where the up-spin and down-spin are separated clearly. But in some cases, e.g. doping and adsorption, the magnetization profile can be modulated and thus noncollinear magnetism can occur. In order to shed light on possible noncollinear magnetism in ZGNR, we study non-collinear magnetism and electronic transport of boron or nitrogen-doped zigzag graphene nanoribbon based on noncollinear density functional theory and non-equilibrium Green's function method. For pristine ZGNR, our results show that the ZGNR presents helical magnetization distribution due to noncollinear magnetization in left and right lead. As the ZGNR is doped with boron and nitrogen atoms, the ZGNR shows a characteristic two-zone feature in the magnetization distribution. Near the dopant site, the magnetic moment of carbon atom is small. However, the magnetic moments of carbon atoms in the left (right) region of dopant are close to those of the left (right) lead. Such a feature provides the possibility of constructing domain walls with various widths on the edge of ZGNR. Moreover, the transmission at the Fermi level (E = 0 eV) decreases with the increase of relative angle between magnetizations of left and right lead, indicating that the spin-flip scattering dominates the electronic transport. However, at E = ±0.65 eV, there is a transmission dip with low transmission, which implies that the dopant induces the strong backscattering. To understand the origin of this dip, we calculate the density of states (DOS) and project the DOS onto each atom of doped ZGNR. The projected DOS shows a large and broad peak at E = 0.65 eV for N-doped ZGNR but at E = +0.65 eV for B-doped ZGNR. The consistency between the position of dip in transmission and the position of peak in DOS indicates that the transmission dip mentioned above is attributed to strong backscattering from the dopant-induced bound state. Our theoretical results are expected to be useful for understanding the noncollinear magnetism and spin scattering in the doped ZGNR-based devices. Also, our work provides a considerable insight into the design of ZGNR-based nanoelectronic devices, such as the transistor based on spin transfer torque effect.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 027101-1 (2019)
Density functional study of metal lithium atom adsorption on antimonene
Xiao-Wei Luan, Jian-Ping Sun, Fan-Song Wang, Hui-Lan Wei, and Yi-Fan Hu

Since the discovery of graphene, researchers have been being increasingly attracted by the emerging of a bunch of two-dimensional (2D) materials, such as BN, MoS2 and black phosphorene. These materials possess outstanding physical and chemical properties, which could find great potential applications in nanoelectronics, energy conversion or storage, photocatalysts, etc. Recently, a theoretically predicted pucker layered material consisting of Sb atoms, antimonene, has been prepared, and is attracting the attention in the field of lithium ion batteries. In this paper, based on first-principle density functional theory, the adsorption characteristics of Li atoms on antimony are studied, including the most stable adsorption configuration, the adsorption density and the diffusion path of Li atom on antimonene. The results show that the most stable adsorption configuration of Li atom is in the valley site, i.e. the center of the three Sb atoms in the top layer and one Sb in the bottom layer. The adsorption energy is 1.69 eV and the adsorption distance is 2.81 . The band structure shows that antimony is an indirect band gap semiconductor with a band gap of 1.08 eV. After the absorption of Li atom, the Fermi level enters into the conduction band, which shows an electronic property similar to metal. The analysis of density of states shows that the p-electronic state of Sb atom and the p and s electronic state of Li atom possess distinct resonance peaks, showing hybrid bonding characteristics. With the increase of the number of Li atoms adsorbed, the lattice structure and electronic structure of antimonene change greatly. The nudged elastic band calculation shows that the diffusion barrier of Li atom on antimony surface is 0.07 eV, and a smaller barrier height is beneficial to the rapid charge-discharge process. To sum up, antimony has a good potential as an anode material for lithium ion batteries.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 026802-1 (2019)
Characteristics of gain in Ne-like Ar 69.8 nm laser pumped by capillary discharge based on double-pass amplification
Tao Liu, Yong-Peng Zhao, Huai-Yu Cui, and Xiao-Lin Liu

In this paper, a double-pass amplification experiment of a Ne-like Ar C line 69.8 nm laser is established. The 45-cmlong capillary is used as the discharge load to obtain a double-pass amplification output of a Ne-like Ar C line 69.8 nm laser. Under the same initial experimental conditions that the initial pressure is 15.4 Pa and the main pulse current amplitude is 13.5 kA, the laser pulse intensity and the full width at half maximum (FWHM) of the laser pulse of the single-pass amplification output and the double-pass amplification output are measured by a vacuum X-ray diode (XRD) behind a vacuum ultraviolet (VUV) monochromator (Acton VSN-515) which is used to disperse the extreme ultraviolet (EUV) emission. And then the laser beam divergence of single-pass amplification output and double-pass amplification output are also measured by a space-resolving flat-field EUV spectrograph combined with an EUV CCD (Andor Newton DO920P-BN). The amplitude of the double-pass amplification laser output is 9 times larger than that of single-pass amplification output, and the FWHM of the double-pass amplification laser pulse is nearly 2.4 ns. While the laser beam divergence angle of the double-pass amplification output is 6.6 times wider than that of single-pass amplification output. By comparing the single-pass amplification and double-pass amplification output experimental results, the gain duration of the gain medium in the double-pass amplification and the radial distribution characteristics of the gain medium are analyzed by using the calculation formula of the double-pass amplification laser intensity. The gain duration is more than 4 ns, during this time the gain coefficient decreases at 1.6 ns. And the gain coefficient is the smallest at 2.8 ns, meanwhile the intensity of the single-pass amplification laser is maximum, and the gain medium is in the gain saturation state. So this result indicates that the minimum gain coefficient at this moment is due to the gain saturation effect. Using a similar calculation method to analyze the spatial distribution of gain coefficients, the gain on the plasma axis is larger than that off the plasma axis. These results lay a foundation for the subsequent establishment of resonant cavity and the multi-pass amplification experiment of capillary discharge Ne-like Ar laser.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 025201-1 (2019)
Signal and pointing accuracy of ultraviolet laser in micro-pattern gaseous detector
Hai-Yun Wang, Hui-Rong Qi, Ling Liu, Zhi-Yang Yuan, Yu-Lian Zhang, Zhi-Wen Wen, Jian Zhang, Yuan-Bo Chen, and Qun Ouyang

In the study of the gas detectors, it is an important calibration method to use the ultraviolet (UV) laser with two-photon ionization mechanism for producing ionized signal. In the last decades, micro pattern gas detector, especially gaseous electron multiplier and micromesh gaseous detector, has been widely used in high energy experiments. These kinds of gaseous detectors have the advantages of higher ion backflow suppression ability, smaller E × B effect and good radiation resistance under the relatively higher count rate environment. To obtain a higher spatial resolution with a UV laser calibration system in gaseous electron multiplier detector, two critical technical issues remain to be resolved: the measurability of the laser signal and the accuracy of the laser beam position. In this paper, the studies in simulation and experiment are conducted to discuss these two critical questions. In the simulation section, the simulation results provide an estimation of signal in the gaseous electron multiplier detector with UV laser of 266 nm wavelength in the mixture working gases of Ar/CO2 (70/30), and give an evaluation of the laser pointing accuracy and the possible relative error of the electron drift velocity. In the experiment section, a UV laser calibration prototype is designed and developed. A pulsed laser of 266 nm wavelength is used as a signal source, which has a Gaussian-like cross section with a frequency of 10 Hz. The experimental results indicate that the signal of the UV laser in a triple gaseous electron multiplier detector reaches 400 mV for a readout strip width of 6 mm, a gain of detector of 5000, and a gain of amplifier of 10 mV/fC. For the calibration laser, the angle accuracy is discussed and tested. The angle uncertainty of the laser can be kept under 5′, and the accuracy of the drift velocity can reach 6.4 × 10 4 with a shift of 0.33 mm in the z direction when the laser beam transmits a distance of 400 mm in the gas chamber. All of these results show that the laser beam specific parameters are the main reference for designing the prototype detector. According to the optimal parameters, a gaseous prototype detector will be tested in the next study.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 022901-1 (2019)
Coherent state phase estimation based on digital nonlinear phase-locked loop
Han Xu, Shu-Xin Chen, Hao Wu, Kun Chen, and Lei Hong

The navigation mechanism of obtaining phase parameters based on quantum theory can break through the limitation of classical physical limit to navigation accuracy. In order to achieve the accurate estimation of the coherent state phase, it is usually assumed that the local oscillator phase must be orthogonal to the coherent state phase in the method of quantum homodyne detection. However, the coherent state phase is unknown and the hypothesis cannot be guaranteed to be correct in practice. In this paper we design a nonlinear phase-locked loop to solve the problem. Firstly, in order to obtain the Wigner distribution for a coherent state, we start with the Wigner distribution for the vacuum state and analyze the noise characteristics according to Wigner distribution of coherent state, then the output of homodyne detection is derived. Secondly, in order to avoid introducing errors in theory, caused by linearization and cope with the limiting requirement between local oscillator phase and coherent state phase in the phase tracking of coherent state, we design an orthogonal simplex cubature Kalman filter (OSCKF) algorithm to achieve the function of the nonlinear phase-locked loop. The algorithm converges by updating the state of the local oscillator phase multiple times, and then, the accurate coherent phase is obtained. Finally, according to the design of the phase-locked loop, we observe the data of homodyne detection and then verify the correctness of the OSCKF algorithm. The simulation results show that the OSCKF can converge to the real phase after observing 200 sampled data and the accuracy is higher than extended Kalmn filter (EKF) and cubature Kalman filter (CKF), and the real phase can be obtained under different local oscillator phase. In conclusion, the nonlinear phase-locked loop based on OSCKF algorithm breaks the limitation of traditional way in which the initial local oscillator phase is required to be orthogonal to the coherent state phase, and effectively avoid the linearization error and improve the anti-nonlinear ability. It is very significant in theory and application .

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024204-1 (2019)
X-ray three-dimensional imaging based on light field imaging technology
Jun-Cheng Qi, Bin Liu, Rong-Chang Chen, Zheng-De Xia, and Ti-Qiao Xiao

X-ray three-dimensional (3D) imaging technology is a research hotspot in the field of X-ray imaging. However, for some special imaging targets, the imaging mode of the traditional computer tomography (CT) circular trajectory is prone to lack of projection information, and thus affects the quality of CT reconstruction images, which limites the application of CT imaging. Light field imaging technology, in which a microlens array is inserted between the sensor and main lens in a traditional camera, achieves four-dimensional (4D) light field data with sensor during imaging including both the two-dimensional (2D) directional information of the radiance propagation and 2D spatial distribution information of object radiation. Through computer calculation imaging, 3D imaging such as digital refocusing, slice in the depth direction, stereo imaging, and depth estimation is realized. This article focuses on the 3D X-ray imaging based on the theory of light field imaging in visible light. Based on the model of parallel X-ray of synchrotron radiation source, the data of the X-ray light field with many projection views are acquired by rotating the image sample. Then, the light passing through any voxel in the imaging target is acquired by a geometric projection method, and based on integral imaging theory of light field imaging, the gray value of the slice in depth dimension is reconstructed and the depth information of reconstructed target is acquired. The reconstruction results show that this method can be used to reconstruct the internal slices at any depth in any viewing direction of the imaging target. In the optical imaging, the scene beyond the depth of field is blurred, making the scene more prominent and the imaging effect better. However, for the X-ray imaging, the imaging mode that is completely transmissive, and the light passing through the foreground carry the information about the background. In the refocusing process, the object at the refocusing depth is focused, and other background information is defocused. Excessive background information overwhelms the real useful information, and makes the slice, especially the edge of the image, blurred. Consequently more severe background noise is introduced due to the defocusing phenomenon in the optical refocusing process. Referring to the reconstruction method of the X-ray 3D imaging and light field imaging, the S-L filter is applied to the original data in the article. After filtering the original data, the X-ray "light field refocusing" is processed. The reconstruction results shown that the method can effectively eliminate reconstruction artifacts and improve image reconstruction quality in the reconstruction depth slice. And in this paper, the light field data are collected by rotating the sample with low time resolution. For the fast imaging, according to the digital refocusing theory of the light field imaging, the array X-ray source and detector can be used. After being calibrated, the system can realize the 3D reconstruction of the light field of the target field with high time resolution. This research has not only the theoretical significance in algorithm, but also great application value in the rapid detection of more complicated targets such as industry and medical treatment.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024202-1 (2019)
Phase delay error analysis of wave plate of division-of-amplitude full Stokes simultaneous polarization imaging system
Yu-Long Yin, Xiao-Bing Sun, Mao-Xin Song, Wei Chen, and Fei-Nan Chen

The division-of-amplitude full Stokes simultaneous polarization imaging system has prominent merits, such as real time, high spatial resolution, high precision, etc. The development of the division-of-amplitude full Stokes simultaneous polarization imaging system has a high application value. The division-of-amplitude full Stokes simultaneous polarization imaging system uses polarization beam splitters, a half wave plate (HWP) and a quarter wave plate (QWP) to modulate the incident Stokes vector into four intensity images. Using the four intensity images, the incident Stokes vector can be analyzed. In the system, the phase delay errors of the HWP and the QWP have a direct influence on the measurement accuracy of the incident Stokes vector. A Stokes vector measurement error equation containing the phase delay errors of the HWP and the QWP is established. When there are the phase delay errors of the HWP and the QWP in the system, the Stokes vector measurement errors of the unpolarized light, 0° liner polarized light, 90° liner polarized light, 45° liner polarized light, 135° liner polarized light, right circularly polarized light and left circularly polarized light are analyzed. A method of solving the Stokes vector measurement error of incident light with any polarization state is given. When the Stokes vectors with different degrees of polarization (DOPs) are used as the incident light, the simulation results show that both the Stokes vector measurement error and the DOP measurement error increase with the DOP of incident light increasing. Therefore, we select the polarization measurement accuracy to evaluate the system when the DOP of incident light equals 1. To ensure that the polarization measurement accuracy of the system is within 2%, the phase delay error of the HWP should be within ±1.6° and the phase delay error of the QWP should be within ±0.5°. The analysis results of the phase delay errors of the HWP and the QWP are of great significance for improving the polarization measurement accuracy of the division-of-amplitude full Stokes simultaneous polarization imaging system, and also provide important theoretical guidance in designing and developing the system.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024203-1 (2019)
Scanning broadband optical parametric chirped pulse amplification based on optical beam deflection
Rong Ye, Zhe-Qiang Zhong, and Xian-Yun Wu

One of the goals pursued in laser pulse is to achieve a laser with a shorter duration and higher intensity. In the past two decades, the laser pulse duration has been shortened by more than 7 orders of magnitude due to the development of Q-switched, Mode-locked and pulse compression technology. The peak power of laser pulse has been increased to PW, even EW and ZW from initial MW with the development of pulse amplification technology, whose focused intensity can reach to 1023 W/cm2. Thus, it provides unprecedented extreme conditions, and speeds up the laser applications in ultrafast nonlinear optics, strong field physics, fast ignition of laser nuclear fusion, optic communication, etc. The optical parametric chirped pulse amplification (OPCPA) is one of the important technologies in ultra-short laser pulse field. It is of great significance to increase the gain bandwidth for improving the conversion efficiency of OPCPA and achieving broadband optical parametric amplification. Combining the optical beam deflection and non-collinear OPCPA, a novel scanning broadband OPCPA model is proposed based on the optical beam deflection. The basic principle of increasing the gain bandwidth for the scanning broadband OPCPA is analyzed theoretically, which ensures the phase matching of each frequency component of signal by optical beam deflecting to change the non-collinear angle constantly. Namely, the non-collinear angles of incident frequency components of signal are different from each other, which, however, makes the whole phase matching of signal, i.e. momentum conservation in optics. The optical parametric amplification of signal pulse with 800 nm central wavelength and almost 100 nm bandwidth is simulated numerically by the proposed scanning broadband OPCPA. The results show that the bandwidth after being amplified is almost the same as before and there is no spectral narrowing, and the scanning broadband OPCPA increases the gain bandwidth and conversion efficiency greatly compared with the amplification with a constant given non-collinear angle, which leads to broadband optical parametric amplification. Finally, it is necessary to make sure that the on-load voltage to the KTN crystal matches with the frequency of signal pulse in time and reduces the unfavorable voltage deviation and time-delay for the maximizing gain bandwidth and conversion efficiency and ensuring the phase matching of each signal frequency component. The results of this paper not only provide an approach to increasing the gain bandwidth of OPCPA, but also supply some theoretical references and the basis for the experimental work of OPCPA in ultra-short laser pulse system.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024205-1 (2019)
Molecular dynamics study of release mechanism of stress at Si/Ge interface on a nanoscale
Xian Chen, Jing Zhang, and Zhao-Huan Tang

In this paper, the stress distribution of Si/Ge interface and the mechanism of the point defect buffer layer in the stress release process of Si/Ge structure on a nanoscale are studied by the molecular dynamics method. The results show that in a one-dimensional interface model, the stress relaxation at the Si/Ge interface is closely related to the size of Ge of the simulation, and there is an approximately linear relationship between the decrease rate of interface stress and the size of Ge. The vacancy defect is introduced into the Ge film near the Si/Ge interface, and a point defect rich buffer layer forms, reducing the stress at the Si/Ge interface significantly. As the defect density in the buffer layer increases, the interface stress decreases step by step. Moreover, in the paper examined also is the mechanism by which the buffer layer affects the interface stress. The relationship between the stress variation caused by the point defect in the pure Ge structure and the defect density is compared and analyzed. The introduction of the defect buffer layer and the increase of the defect density can accelerate the release of the Si/Ge interface stress. The internal mechanism of the stress reduction would be that the defect introduces the compressive stress, which can offset the tensile stress generated by the lattice mismatch of the Si and Ge structure. Then the Si/Ge interface stress is reduced. Besides, as the defect density increases, the compressive stress introduced by the defect increases and the interface stress decreases. The examination shows that the introduction of the defect buffer layer into the growth of the pure Ge film on silicon can reduce the probability of the dislocation defect by structuring the buffer layer to reduce interface stress. And this method is indirectly confirmed by preliminary study of silicon-based Ge film growth with low-temperature buffer layer method.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 026801-1 (2019)
Electrical conductivity and infrared ray photoconductivity for lattice distorted SmNiO3 perovskite oxide film
Hai-Yang Hu, Ji-Kun Chen, Fei Shao, Yong Wu, Kang-Kang Meng, Zhi-Peng Li, Jun Miao, Xiao-Guang Xu, Jia-Ou Wang, and Yong Jiang

The metal-to-insulator transitions achieved in rare-earth nickelate (RNiO3) receive considerable attentions owning to their potential applications in areas such as temperature sensors, non-volatile memory devices, electronic switches, etc. In contrast to conventional semiconductors, the RNiO3 is a typical electron correlation system, in which the electronic band structure is dominant by the Coulomb energy relating to the d-band and its hybridized orbitals. It was previously pointed out that lattice distortion can largely influence the electronic band structures and further significantly affect the electronic transportation properties, such as the resistivity and metal-to-insulator transition properties. Apart from directly measuring the transportation performance, the variations in the origin of carrier conduction and orbital transitions relating to the strain distortion of RNiO3 can also be reflected via their optical properties. In this work, we investigate the optical properties of samarium nickel (SmNiO3) thin films when lattice distortions are induced by interfacial strains. To introduce the interfacial strain, the SmNiO3 thin films are epitaxially grown on the strontium titanate (SrTiO3) and lanthanum aluminate (LaAlO3) single crystal substrates by using the pulsed laser deposition. A bi-axial tensile distortion happens when the SmNiO3 thin films are grown on SrTiO3 due to the smaller lattice constant of SmNiO3 than that of SrTiO3, while the one grown on LaAlO3 is strain-relaxed. We measure the infrared radiation (IR) transmission spectra of the SmNiO3 thin films grown on various substrates. The obtained IR transmission spectra are fitted by a Drude-Lorentz model and further converted into the curves of photoconductivity versus IR frequency. Comparing the difference in photoconductance between low frequency and high frequency reflects the two different origins of the conduction, which are related to intraband transition and band-to-band transition, respectively. The smaller photoconductance is observed for SmNiO3/SrTiO3 than for SmNiO3/LaAlO3 at low frequency, and this is expected to be caused by the suppression of free carriers as reported previously for tensile distorted SmNiO3. The consistence is obtained when further measuring the electronic transportation such as temperature-dependent electrical resistivity, as a higher resistivity is observed for SmNiO3/SrTiO3 than for SmNiO3/LaAlO3. The combination of the investigation of electrical transport with that of infrared transmission indicates that the tensile distortion in structure stabilizes the insulating phase to eliminate a pronounced metal-to-insulator transition and elevates the transition temperature. This is related to the respective twisting of the NiO6 octahedron when tensile distortion regulates the valance state of the transition metal and further opens the band gap, which is further confirmed by results of the X-ray absorption spectrum.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 026701-1 (2019)
Acoustic radiation from a cylinder in shallow water by finite element-parabolic equation method
Zhi-Wen Qian, De-Jiang Shang, Qi-Hang Sun, Yuan-An He, and Jing-Sheng Zhai

It can be a difficult problem to precisely predict the acoustic field radiated from a finite elastic structure in shallow water channel because of its strong coupling with up-down boundaries and the fluid medium, whose acoustic field cannot be calculated directly by existing methods, such as Ray theory, normal mode theory and other different methods, which are adaptable to sound fields from idealized point sources in waveguide. So, there is no reliable research method of predicting the acoustic radiation of elastic structure in shallow water at present. Based on the finite element method (FEM) coupled with the parabolic equation (PE), the theoretical model for structure acoustic radiation in shallow water at low frequency is established in this paper. This model mainly consists of three sections. First, obtaining the near-field vibro-acoustic characteristics of the elastic structure in shallow water by the multi-physics coupling model established by FEM, whose FEM model includes the up-down boundaries and the completely absorbent sound boundaries in the horizontal direction. Second, getting the acoustic information in the depth, which is set as the acoustic input condition i.e. starting field for the PE. Third, the acoustic information in the far-field quickly calculated by the PE and the finite difference method (FDM). The accuracy, efficiency and fast convergence of FEM-PE method are validated by numerical simulation and theoretical analysis through using a monopole source and structural source in the Pekeris waveguide, respectively. The vibro-acoustic characteristics of elastic cylinder influenced by upper and lower fluid boundaries of the Pekeris waveguide are calculated and analyzed. The cylindrical shell material is steel, and it is 1 m in radius and 10 m in length. The shallow water channel is a Pekeris waveguide with 30 m in depth, at the upper boundary, i.e., the free surface, the lower boundary is the semi-infinite liquid boundary. The analyzed frequencies range from 50 Hz to 200 Hz. The study shows that when the cylindrical shell approaches to the sea surface or bottom, the coupled frequency is higher or lower respectively than that of the shell immersed in the free field. When the diving depth reaches a certain distance range, the coupled frequency tends to be the same as that in free field. The acoustic field radiated from an elastic shell in Pekeris waveguide is similar to that from a point source at low frequency, but there exists a significant difference in high frequency between them, so the structural source can be equivalent to a point source conditionally. The sound radiation attenuation of the structure happens in sequence according to the near-field acoustic shadow zone, the spherical wave attenuation zone, the region between spherical wave and the cylindrical wave attenuation zone, and the cylindrical wave attenuation zone.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024301-1 (2019)
Two-dimensional silicon annular photonic crystals for realizing polarization-independent unidirectional transmission
Dan Liu, and Sen Hu

Optical diode is a device that can realize unidirectional transmission of light. Its function is similar to that of an electronic diode. It has important applications in the field of optoelectronic integration and all-optical communications. Unidirectional wave transmission requires either time-reversal or spatial inversion symmetry breaking. The magneto-optical effect and optical nonlinearity are usually utilized to break the time-reversal symmetry and obtain the unidirectional transmission. However, these schemes all need high light intensity or magnetic field strength to be realized, and limit the usage. Therefore, spatial inversion symmetry breaking is highly desirable because of totally linear materials under low intensities. Quit a lot of researchers have designed optical diodes based on the photonic crystals and achieved unidirectional transmission for TE-like or TM-like light. The early design realized light unidirectional transmission by PC structures for only one polarization state (TE-like or TM-like incident light). It limits the application for the high integration and reconfigurable optical interconnection. The structure which can achieve unidirectional transmission for both TE and TM polarizations needs to be designed. The annular PCs have been verified to realize polarization-independent phenomena, such as beam splitting, self collimation and waveguide. In this paper, an annular PC is proposed. The plane wave expansion method is used to calculate band structures. The results show that it exhibits a significant directional band gap for both TE and TM mode. Then, the triangular annular PC is constructed, and its transmission spectra and field distributions are calculated by the finite-different time-domain method. It is found that the structure can realize the polarization-independent unidirectional transmission, but the forward transmissivity is too low (about 20%). Moreover, another smaller size annular PC is further introduced to form annular PC heterojunction, which effectively improves the polarization-independent unidirectional transmission performance and the forward transmissivity has doubled. Through the adjustment of the interface structure, the forward transmissivity is further increased. The optimized annular PC heterostructure can realize polarization-independent unidirectional transmission, and the forward transmissivity reaches 44%. The heterostructure can be used to fabricate polarization-independent optical diode, and may have potential applications in complex all-optical integrated circuits.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024206-1 (2019)
Ultrafast azimuthal beam smoothing scheme based on vortex beam
Bo-Yu Tian, Zhe-Qiang Zhong, Zhan Sui, Bin Zhang, and Xiao Yuan

The illumination uniformity of laser beams in inertial confinement fusion (ICF) facility is a key factor, which plays a crucial role in suppressing the laser plasma instabilities. However, the prevailing beam smoothing techniques cannot meet all the requirements for improving the irradiance uniformity of laser beams and mitigating the laser plasma instabilities, which are determined by the high-frequency spatial modulations and the fine-scale speckles of the focal spots. An ultrafast azimuthal beam smoothing scheme based on vortex beams is proposed in this paper. In this scheme, two of the four beams in a laser quad are transformed from super-Gaussian (SG) beams into vortex beams by inserting two spiral phase plates with opposite topological charges into the beam path, whereas the other two SG beams remain unchanged. By controlling the polarization and the center wavelength of each beam, the SG beam and the transformed vortex beam in the quad are coherently superposed on the target plane, so are the remaining two beams. Owing to the difference in central wavelength and the existence of the topological charges, two focal spots rotating in a period of a few picoseconds are generated in the target plane, which can redistribute the speckles quickly in temporal domain and thus improve the irradiance uniformity of the laser quad. By establishing the physical model of the azimuthal smoothing scheme, the smoothing characteristics including the rotation period, the illumination uniformity and the fractional-power-above-intensity of the focal spots are analyzed in detail. In order to improve the smoothing characteristics significantly, the novel smoothing scheme is further combined with another ultrafast smoothing scheme, i.e. radial smoothing scheme. The influence of the key parameters of the combined smoothing scheme on the illumination uniformity and on the smoothing velocity are discussed. Results indicate that the azimuthal smoothing scheme can achieve the ultrafast smooth of the laser quad in the azimuthal direction and the best illumination uniformity within a few picoseconds as well. Though the degree of improvement in the irradiance uniformity of the azimuthal smoothing scheme is lower than that of the radial smoothing, the combination of these two schemes can improve the uniformity effectively and rapidly. The novel smoothing scheme provides a potential smoothing approach for the high-power laser facilities.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 024207-1 (2019)
Phase transition of S4 model on a family of diamond lattice
Xun-Chang Yin, Wan-Fang Liu, Ye-Wan Ma, Xiang-Mu Kong, Jun Wen, and Li-Hua Zhang

The fractal is a kind of geometric figure with self-similar character. Phase transition and critical phenomenon of spin model on fractal lattice have been widely studied and many interesting results have been obtained. The \begin{document}${S^4}$\end{document} model regarded as an extension of the Ising model, can take a continuous spin value. Research of the \begin{document}${S^4}$\end{document} model can give a better understanding of the phase transition in the real ferromagnetic system in nature. In previous work, the phase transition of the \begin{document}${S^4}$\end{document} model on the translation symmetry lattice has been studied with the momentum space renormalization group technique. It is found that the number of the fixed points is related to the space dimensionality. In this paper, we generate a family of diamond hierarchical lattices. The lattice is a typical inhomogenous fractal with self-similar character, whose fractal dimensionality and the order of ramification are \begin{document}${d_{\rm{f}}} = {\rm{1}} + \ln m/\ln {\rm{3}}$\end{document} and \begin{document}$R = \infty $\end{document}, respectively. In order to discuss the phase transition of the \begin{document}${S^4}$\end{document} model on the lattice, we assume that the Gaussian distribution constant \begin{document}${b_i}$\end{document} and the fourth-order interaction parameter \begin{document}${u_i}$\end{document} depend on the coordination number \begin{document}${q_i}$\end{document}<graphic xlink:href="20181315e-wen-revis

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 026401-1 (2019)
Dependence of performance character of photon-atom entanglement source on retrieval efficiency
Sheng-Zhi Wang, Ya-Fei Wen, Chang-Rui Zhang, Deng-Xin Wang, Zhong-Xiao Xu, Shu-Jing Li, and Hai Wang

The photon-atom interface is a basic component of quantum repeater, quantum network, and linear optical quantum computing. Different approaches have been tested in the last decade to develop quantum interface, such as quantum dots, single atoms and ions, color centers and cold atomic ensemble. In the cold atomic ensemble, a normal way to produce photon-atom interface is the Duan-Lukin-Cirac-Zoller (DLCZ) protocol. Used in the DLCZ protocol is an atomic ensemble that can emit single photons while creating a single atomic excitation, which is stored in the ensemble. The atomic excitation can be converted into a photon due to the collective interference. The influences of the retrieval efficiency on the atom-photon entanglement source have been studied in various experiments. But no one has studied the retrieval efficiency threshold of entanglement generation. In our experiment we study the retrieval efficiency dependence on read power and OD. Setting the power of the repump light beam to be 12.2 mW, 5.0 mW, 2.0 mW, 0.5 mW and 0.3 mW, OD of the cold atom ensemble is measured to be 20, 17, 10, 2, and 1, respectively. As we expected, the retrieval efficiency increases with increasing OD value and read power, the curve shows that the retrieval efficiency increases sharply with increasing the OD value and read power, then after a while slowly increases with increasing the OD values and read power. Then we measure the Bell parameter with increasing the retrieval efficiency by increasing the read power. It shows that the Bell parameter sharply increases for retrieval efficiency values ranging from 0 to 3%, but changes very small for retrieval efficiency values ranging from 3% to 18.3%. The maximum Bell parameter is 2.6. We further analysis the result, finding that the Bell parameter can be expressed as \begin{document}$S = \dfrac{{{S_{{\rm{MAX}}}}r}}{{(1 + 2\chi )r + 2B}}$\end{document}. Fitting parameters to the curve are \begin{document}$\chi$\end{document}= 1%, B = 0.073%. To avoid of multi-excitation the write power kept low that \begin{document}$\chi$\end{document} at 1% level. Then we can find out from the function that the signal-to-noise ratio is bigger than 6∶1 the Bell parameter will reach 2. The theoretical analysis and experimental results fit very well. So the further reason that alter the Bell parameter is the signal-to-noise ratio. We should decrease the noise while increasing the retrieval efficiency. This paper will help with rise the quality of entanglement generation through photon-atom interface.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 020301-1 (2019)
Magnetism and electronic properties of LiFeAs superconducting thin filma under two-dimensional strains effect
Xin Wang, Hua Li, Zheng-Chao Dong, and Chong-Gui Zhong

The magnetism, band properties and electronic density of states of LiFeAs superconducting thin film with two-dimensional strain are investigated by using the first principles calculations based on density functional theory, and the influences of different strains on the characteristics of superconducting films are analyzed in detail. The results show that the magnetic ground configuration is the striped antiferromagnetic state of nostrained LiFeAs thin film, and the ground structure of this system is unchanged in the range of applied 1% 6% compressive and tensile strain. The density of states near the Fermi level is mainly from the contribution of Fe-3d orbital and a few As-4p electrons. The electron spin exchange coupling between Fe ions is realized by As ions. Furthermore, unlike the case of the nostrain and the tensile strain, with increasing the compressive strain, the localized antiparallel electron spin magnetic moments of Fe ion decrease, the density of states at the Fermi surface improves, and the itinerant electron magnetism of Fe ions increases, which all greatly suppress the antiferromagnetic properties of thin film and enhance the superconducting phase transition temperature. The superconductivity of LiFeAs thin film originates from the Cooper pairs of electrons between the hole-type and electronic-type bands near the Fermi surface through the antiferromagnetic superexchange coupling effect. Instead, the LiFeAs thin film with the tensile strain presents completely opposite properties, that is to say, the decrease of the electronic density of states in the Fermi level brings about the weakening of the metal properties and the increasing of the antiferromagnetic exchange coupling. Particularly, the band structure of hole-type near the Fermi surface disappears, and the occurrence of Cooper pairs of electrons becomes significantly reduced, resulting in the suppressed superconducting phase transition when the LiFeAs thin film is subjected to tensile strain. In addition, the change of antiferromagnetic exchange coupling and magnetic moments of Fe ions are also explained according to the variation of electronic density of states of the Fe-3d energy levels during the distortion of FeAs tetrahedrons due to compressive strain. In brief, our researches provide an effective way to improve the superconducting properties of LiFeAs thin film and may promote the relevant practical applications of iron-based superconductors in the future.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 027401-1 (2019)
Study of pedestrian flow on stairs with a cellular transmission model
Hui Jin, and Ren-Yong Guo

The aim of this study is to address the following issues: 1) revealing the typical behaviors and properties of pedestrian movement when going upstairs and downstairs; 2) constructing a pedestrian evacuation model to formulate the walking process of pedestrians in stair area; 3) verifying that the cell transmission model widely used in the two-dimensional walking space can also be applied to the three-dimensional staircase area. Firstly, an observation experiment is carried out to gain the pedestrian movement data in the process of going upstairs and downstairs. By collating the data, the relation between density and flow in the unidirectional process of going upstairs or going downstairs, and in the bi-directional process of going upstairs and downstairs, are drawn respectively. Then, by analyzing the fundamental diagrams, several characteristics of pedestrian movement in stair area are revealed. Based on these characteristics, an extended cell transmission model is proposed. In this model, a potential correction coefficient is introduced to change the route choice of pedestrians by using the influence of different directional pedestrians on the potential; a flow modification coefficient is introduced to describe the effect of physical parameters on the maximum flow at the boundary between two neighboring cells; and an offset coefficient is introduced to correct movement rules and strengthen the influence of preferential direction on pedestrian route choice. Further, simulations relied on the proposed model are conducted. By comparing the simulation results with the experimental data, the model is calibrated. Then the calibrated model is employed to formulate the pedestrian movement in stair area, and the sensitivity of the potential correction parameter is also discussed. The simulation results indicate that the proposed model can successfully reproduce the movement of pedestrians on stair. Moreover, the route-choice behaviors of pedestrians can be directed by varying the values of the potential correction coefficient, which can present important information about optimizing the evacuation process of pedestrians on stair, thereby reducing the risk of an accident, such as congesting and treading.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 2 020501-1 (2019)
Tuning of cyclotron resonant frequency and magnetic field of medical cyclotron
Yong-Cun Wu, Xing-Lin Yang, Jin-Shui Shi, Liang-Chao Zhao, and Xiao-Zhong He

The high-frequency resonant cavity is affected by factors such as beam load, gravity and heat loss caused by cavity deformation during the actual operation of the medical cyclotron. The resonant frequency will shift to a certain extent, resulting in the high-frequency operation frequency varying with the resonant frequency of the resonator cavity. In order to meet the requirements for isochronous acceleration, the magnetic field strength should also be changed correspondingly when the high-frequency operation frequency changes, that is, the magnitude of the magnet current needs changing accordingly, so that the particle cyclotron frequency matches the high-frequency resonant frequency to overcome the sliding phase. Firstly, the static magnetic field model is established by finite element simulation software to simulate the average magnetic field of cyclotron under different magnet currents. Then the relationship between the magnetic field and the resonant frequency is theoretically analyzed. Finally, the relationship between the magnet current and the resonant frequency is obtained when the magnet current varies in a small interval. According to the optimal magnet current corresponding to different resonance frequencies, the automatic frequency tracking of magnet current is completed. In the case of ensuring the maximum carbon film beam, the optimal magnet current corresponding to different resonance frequencies is obtained, which makes the theory validated. According to the relationship, the magnet current is automatically adjusted, which overcomes the slip phase and ensures the stable output of the Faraday beam. The method enables the magnet current to be quickly and accurately find and track the cavity frequency, overcomes the slip phase caused by the frequency offset, and completes the stable output of the beam.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 12 122901-1 (2019)
New strategy for global tallying in Monte Carlo criticality calculation
Dan-Hua Shangguan, Zhi-Cheng Ji, Li Deng, Rui Li, Gang Li, and Yuan-Guang Fu

Traditionally, the Monte Carlo criticality calculation must set a maximum inactive step by experience to ensure that a fission source distribution has converged. The tallying process can only be invoked after this maximum inactive step to avoid the system error caused by the non-converged fission source distribution. In the same way, the uniform fission site algorithm for increasing the whole efficiency of global tallying should also be invoked after the fission source distribution has converged fully. The calculation must reach a maximum iteration step, then, this process can be stopped and the tallies can be printed. This old strategy has two defects. Firstly, the appointed maximum inactive step can only be set by experience, which will be insufficient in some cases; secondly, some iteration steps can be wasted because the precision of tallies has been enough and no one knows it. So, a new strategy is proposed in this article to overcome these defects. Based on an on-the-fly diagnostic method for the convergence of Shannon entropy sequence corresponding to the fission source distribution of each iteration step, the uniform fission site algorithm will be invoked after the iteration step whose serial number is the maximum of the first active step and the first converged step diagnosed by the above-mentioned rule. This rule will be helpful in ensuring that the uniform fission site algorithm can use enough accurate data to bias the secondary fission neutron number, thus avoiding the system error to some degree. Then, a global precision index will be calculated at each fixed step to judge whether the precision standard is reached. If so, the whole calculation is stopped. This process will be repeated until the pre-set maximum step number is reached. In this way, superfluous calculations can be skipped but the calculation precision can be guaranteed. In a word, this new strategy is beneficial to increasing the efficiency of global tallying in the Monte Carlo criticality calculation when appropriate parameters are adopted. This conclusion can be proved by the numerical result from the C5G7 benchmark model.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 12 122801-1 (2019)
Time of flight technology based on multi-gap resistive plate chamber
Yi Wang, Qiu-Nan Zhang, Dong Han, and Yuan-Jing Li

Particle identification is very important in nuclear and particle physics experiments. Time of flight system (TOF) plays an important role in particle identification such as the separation of pion, kaon and proton. Multi-gap resistive plate chamber (MRPC) is a new kind of avalanche gas detector and it has excellent time resolution power. The intrinsic time resolution of narrow gap MRPC is less than 10 ps. So the MRPC technology TOF system is widely used in modern physics experiments for particle identification. With the increase of accelerator energy and luminosity, the TOF system is required to indentify definite particles precisely under high rate environment. The MRPC technology TOF system can be defined as three generations according to the timing and rate requirement. The first-generation TOF is based on the float glass MRPC and its time resolution is around 80 ps, but the rate is relatively low (typically lower than 100 Hz/cm2). The typical systems are TOF of RHIC-STAR, LHC-ALICE and BES III endcap. For the second-generation TOF, its time resolution has the same order as that for the first generation, but the rate capability is much higher. Its rate capability can reach 30 kHz/cm2. The typical experiment with this high rate TOF is FAIR-CBM. The biggest challenge is in the third-generation TOF. For example, the momentum upper limit of $ {\rm{K}}/{\text{π}}$ separation is around 7 GeV/c for JLab-SoLID TOF system under high particle rate as high as 20 kHz/cm2, and the time requirement is around 20 ps. The readout electronics of first two generations is based on time over threshold method, and pulse shape sampling technology will be used in the third-generation TOF. In the same time, the machine learning technology LSTM network is also used to analyze the time performance. As a very successful sample, MRPC barrel TOF has been used in RHIC-STAR for more than ten years and many important physics results have been obtained. A prominent result is the observation of antimatter helium-4 nucleus. This discovery proves the existence of antimatter in the early universe. In this paper, we will describe the evolution of MRPC TOF technology and key technology of each generation of TOFs including MRPC detector and related electronics. The industrial and medical usage of MRPC are also introduced in the work finally.

Acta Physica Sinica
Jan. 01, 1900, Vol. 68 Issue 10 102901-1 (2019)
Please enter the answer below before you can view the full text.