x

## Special topic

DOI: 10.7498/aps.70.010101

超导及其应用专题编者按DOI: 10.7498/aps.70.010101 　1911 年荷兰科学家Heike Kamerlingh Onnes 首次在金属汞中发现超导现象以来, 超导作为人类发现的第一个宏观量子现象已经有百余年的研究历史. 在这百余年的时间里, 人们对传统的低温超导材料的认识及应用已经取得了巨大的成就, 尤其是关于其超导机理的BCS 理论的建立极大地推动了凝聚态物理的发展. 在铜氧化物高温超导体发现后的近三十余年里, 源于对其机理的研究开辟了基础物理新的领域, 也为超导体的应用带来了新的技术. 然而, 非常规高温超导机理的研究和高临界参数的新超导体的探索仍面临许多挑战.

20 世纪80 年代, 铜氧化合物高温超导体的发现为超导研究开辟了一个全新的领域, 在全世界范围内掀起了高温超导研究的热潮. 以赵忠贤院士为代表的中国科学家在铜氧化物高温超导体研究领域做出了重要贡献, 为我国在国际超导界赢得了一席之地. 因铜氧化物高温超导体的发现, 1987 年“三月会议”(March Meeting) 盛况空前, 堪比于音乐界的“摇滚乐狂欢节” (“WoodstockFestival”). 在这次会议上, 赵忠贤先生作为会议特邀的五个科学家之一代表我国介绍了工作,确立了我国在高温超导研究领域中的重要国际地位. 当时, 以赵忠贤先生为代表的老一辈科学家在高温超导研究中取得的成绩激励了国内一大批年轻人(也包括笔者和本专题中许多文章的撰写人) 从事超导相关的研究, 这批“年轻人”现在也大多成为了我国在超导研究领域的中坚力量. 多年来, 国家对基础研究持续投入, 超导研究的基础条件有了长足的进步. 同时, 培养了一批在超导研究领域中具有国际影响力的学者, 我国的超导研究也逐渐走到了世界的前列. 2008 年, 赵忠贤先生和经历铜氧化物高温超导研究培养出来的中坚力量齐心协力, 在铁基高温超导体的研究中再次取得重大突破, 引领了高温超导研究的第二次“热潮”, 实现了我国在高温超导研究领域的全面赶超.

在经历了两次高温超导研究的“热潮”之后, 超导研究领域又面临着新的思考, 未来超导研究将如何发展, 这是摆在超导研究人员面前一个必须回答的问题. 这使我想起了早在1986 年铜氧化物高温超导体发现之前, 《物理》杂志1977 年刊登了赵忠贤先生题为“探索高临界温度超导体”的文章(见: 赵忠贤 1977 《物理》6(4): 211—215). 也许正是这种前瞻性的思考和积累, 才有我国在1986 年和2008 年两次非常规高温超导体研究中取得的成绩和地位. 我们应该有前瞻性的思维和布局, 超导作为一种宏观量子现象, 非常规高温超导的物理机理还没有完全解决, 这不但需要理论物理学家的努力, 也需要实验物理学家和材料科学家的齐心协力. 在未来的若干年里, 这一领域仍然是超导研究的重心之一, 高温超导体新材料的发现是极可能实现重大突破的方向. 在应用方面, 超导材料已经被广泛地应用于我们的生活、科研和生产等许多方面(如: 医院核磁共振成像、大科学装置和实验室的各种超导磁体等), 尤其是高温超导体已开始得到应用, 如超导电磁感应加热应用于铝锭的加工和电网限流器等. 超导材料的广泛应用将会极大改善我们人类的生活品质, 大力开展超导应用领域的研究将是未来超导研究的另一个重要方向. 此外, 超导体还具备许多特殊的“新性质”,例如拓扑超导被视为实现未来量子计算的重要方案之一. 超导体的这些“新性质”也将在未来有重要的发展, 这些“新性质”的潜在应用可提供相关的超导解决方案, 从而实现一些核心技术方面的突破,对其开展研究是未来超导研究的又一重要方向.

为进一步促进我国在超导研究领域的发展, 作为我国中文物理类学术期刊中影响面最广、影响力最大的刊物, 《物理学报》组织出版了“超导及其应用”专题. 本专题针对上述超导研究方向, 邀请了相关领域国内有代表性的学者进行专题评论, 很好地梳理了相关方向的近期进展, 内容紧贴当前的研究前沿, 为我国超导领域的研究人员提供了一个全面的参考资料, 非常及时和必要. 最后, 希望以此专题为契机, 激发国内超导领域的同行对超导未来的发展开展广泛的讨论和思考, 促进我国超导研究事业全面发展.

## EDITOR'S SUGGESTION

2021, 70 (1): 010101. doi: 10.7498/aps.70.010101
Abstract +

## EDITOR'S SUGGESTION

2021, 70 (1): 017101. doi: 10.7498/aps.70.20202122
Abstract +
Based on the common properties exhibited in both cuprates and iron-based high temperature superconductors, we have recently proposed the “gene” concept for unconventional high temperature superconductors: those d-orbitals of transition metal elements with the strongest in-plane bonding to anion p-orbitals must be isolated near Fermi energy. Here we summarized recent progress in this research direction and discussed several electronic environments that meet the “gene” condition. We also discussed the challenge and the possibility in finding new unconventional high temperature superconductors.

## EDITOR'S SUGGESTION

2021, 70 (1): 017408. doi: 10.7498/aps.70.20202180
Abstract +
High-Tc cuprates, iron-based superconductors, heavy-fermion superconductors and κ-type layered organic superconductors share some common features − the proximity of the superconducting state to the magnetic ordered state and the non-s-wave superconducting pairing function. It is generally believed that the Cooper pairings in these unconventional superconductors are mediated by spin fluctuations. In this paper, we present a brief overview on the spin dynamics and unconventional pairing, focusing on high-Tc cuprates and iron-based superconductors. In particular, we will overview the properties of the neutron spin resonance and its possible origin, the pairing mechanism in Hubbard model within the weak-coupling framework and its application to the aforesaid unconventional superconductors. We point out that the interplay between magnetism and superconductivity is still an area of active research.

## EDITOR'S SUGGESTION

2021, 70 (1): 017406. doi: 10.7498/aps.70.20201913
Abstract +
Superconductivity represents a magic macroscopic quantum phenomenon. There have been two major categories of superconductors: the conventional superconductors represented by metals or alloys; and the unconventional superconductors represented by cuprates and iron-based high-temperature superconductors. While the superconductivity mechanism of the conventional superconductors is successfully addressed by the BCS theory of superconductivity, no consensus has been reached in understanding the high temperature superconductivity mechanism for more than 30 years, which has become one of the most prominent issues in condensed matter physics. Revealing the microscopic electronic structure of unconventional superconductors is the prerequisite and foundation in understanding their superconductivity. Angle resolved photoelectron spectroscopy (ARPES) plays an important role in the study of unconventional superconductors because it can directly measure the electronic structure of materials. In this paper, our recent progress in the ARPES study of electronic structure and superconductivity mechanism of high temperature cuprate superconductors and iron-based superconductors is reviewed. It mainly includes the electronic structure of the parent compound, the non-Fermi liquid behavior in the normal state, the band and gap structure of the superconducting state, and the many-body interactions both in the normal and superconducting states. These results will provide important information in understanding the superconductivity mechanism of Cu-based and Fe-based superconductors.

## EDITOR'S SUGGESTION

2021, 70 (1): 017403. doi: 10.7498/aps.70.20202102
Abstract +
Cuprate and iron-based superconductors are known as the only two types of high-Tc superconductors. The mechanism of high-Tc superconductivity is the most challenging issue in the field. Building accurate high-dimensional phase diagram and exploring key parameters that determine Tc, would be essential to the comprehension of high-Tc mechanism. The electronic phase diagrams of cuprate superconductors show complexity and diversity, for the strong coupling and interplay among lattice, orbital, charge and spin degrees of freedom. It is tough to construct a high-dimensional holographic phase diagram and obtain quantitative laws by traditional research methods. Fortunately, the high-throughput synthesis and fast screening techniques enable to probe the phase diagram via line-by-line or map scanning modes, and thereby are expected to obtain high-dimensional phase diagram and key superconducting parameters in a much efficient way. In this article, electronic phase diagrams of cuprate superconductors that are obtained mainly by electrical transport measurements, are briefly summarized in the view of cation substitutions, oxygen variation in the parent compounds, electric double-layer gating (electrostatic/electrochemical manipulation) and magnetic field. We introduce the preparation methods for combinatorial film based on the developed pulsed laser deposition and oxide molecular beam epitaxy techniques, as well as corresponding scale-span high-throughput measurement techniques. These high-throughput techniques have been successfully applied in the research of interface superconductivity, quantum phase transition, and so on. The novel high-throughput superconductivity research mode will play an indispensable role in the construction of the high-dimensional holographic phase diagram, the comprehension of high-Tc mechanism, and practical applications of superconductors.

## EDITOR'S SUGGESTION

2021, 70 (1): 017404. doi: 10.7498/aps.70.20201836
Abstract +
$T_{\rm{c}}$ cuprates. In the hole-doped iron-based superconductors, the Hall coefficient changes its sign in low temperatures, and meanwhile the resistivity shows a broad hump in the same temperature range. Such a behavior is proposed as a crossover from incoherent to coherent transport. The Seebeck coefficients of iron-based superconductors also show remarkable differences from the cuprates. In iron-based superconductors, the absolute value of Seebeck coefficients in the normal state becomes the largest at the optimally doping point with highest $T_{\rm{c}}$, which is probably related to the strong inter-band scattering. The Nernst effect in the normal state of iron-based superconductors indicates that superconducting phase fluctuations is not obvious above $T_{\rm{c}}$, which is also significantly different from the cuprates. These unusual thermoelectric properties observed in iron-based superconductors have not been observed in the nickel-based pnictide superconductors with the analogous structure, i.e., LaNiAsO, and the nickel-based superconductors behave more like a usual metal. All these results above illustrate that these unusual transport properties of iron-based superconductors are inherently associated with their high temperature superconductivity, and these factors should be taken into account in the theory on its superconducting mechanism.">There are a variety of order states in iron-based pnictides, such as electronic nematic phase, spin density wave, and so on, which leads to plenty of novel physical phenomena. The measurements of transport properties can provide extremely useful information for understanding of the low-energy excitations of iron-based superconductors. Due to the multi-band electronic structure in iron-based pnictides, the temperature dependence of resistivity and Hall coefficient varies with different systems, however, there are no evidence for the pseudo-gap opening in the normal state which is a common feature in underdoped high-$T_{\rm{c}}$ cuprates. In the hole-doped iron-based superconductors, the Hall coefficient changes its sign in low temperatures, and meanwhile the resistivity shows a broad hump in the same temperature range. Such a behavior is proposed as a crossover from incoherent to coherent transport. The Seebeck coefficients of iron-based superconductors also show remarkable differences from the cuprates. In iron-based superconductors, the absolute value of Seebeck coefficients in the normal state becomes the largest at the optimally doping point with highest $T_{\rm{c}}$, which is probably related to the strong inter-band scattering. The Nernst effect in the normal state of iron-based superconductors indicates that superconducting phase fluctuations is not obvious above $T_{\rm{c}}$, which is also significantly different from the cuprates. These unusual thermoelectric properties observed in iron-based superconductors have not been observed in the nickel-based pnictide superconductors with the analogous structure, i.e., LaNiAsO, and the nickel-based superconductors behave more like a usual metal. All these results above illustrate that these unusual transport properties of iron-based superconductors are inherently associated with their high temperature superconductivity, and these factors should be taken into account in the theory on its superconducting mechanism.

## EDITOR'S SUGGESTION

2021, 70 (1): 017401. doi: 10.7498/aps.70.20201673
Abstract +
As a novel quantum state in condensed matter physics, Majorana zero mode has become a popular research topic at present because of its potential value in topological quantum computing. Theory predicts that Majorana zero mode appears in the vortex core of the topological superconductor as a unique bound state. However, due to various factors such as the existence of conventional low energy bound states or impurity states, it is difficult to identify the Majorana zero mode and to put it into the specific applications. Nowadays, it is still urgent to find a suitable topological superconducting system and identify the clean Majorana zero mode in experiment. In this paper, we study the vortex states of electron-doped iron-selenium-based superconductors (Li, Fe)OHFeSe and single-layer FeSe/SrTiO3 with extremely high energy resolution STM. There exists a robust and clean Majorana zero mode in the free vortex core of (Li, Fe)OHFeSe, which has the quantized conductance. As for single-layer FeSe/SrTiO3 film, it has only conventional Caroli-de Gennes-Matricon (CdGM) bound states without zero energy mode. These experimental results provide a suitable platform for further studying the physical properties of Majorana zero mode, and also shed light on the source of topological superconductivity in iron-based superconductors.

## EDITOR'S SUGGESTION

2021, 70 (1): 017402. doi: 10.7498/aps.70.20201418
Abstract +
Heavy fermion superconductors belong to a special class of strongly correlated systems and unconventional superconductors. The emergence of superconductivity in these materials is closely associated with the presence of quantum critical fluctuations. Heavy fermion superconductors of different structures often exhibit distinct competing orders and superconducting phase diagrams, implying sensitive dependence of their electronic structures and pairing mechanism on the crystal symmetry. Here we give a brief introduction on recent theoretical and experimental progress in several different material families. We develop a new phenomenological framework of superconductivity combining the Eliashberg theory, a phenomenological form of quantum critical fluctuations, and strongly correlated band structure calculations for real materials. Our theory provides a unified way for systematic understanding of various heavy fermion superconductors.

## EDITOR'S SUGGESTION

2021, 70 (1): 017407. doi: 10.7498/aps.70.20202189
Abstract +
In recent years, hydrogen-rich compounds under extremely high pressure have become the hot target materials for high-temperature superconductors. At present, two landmark progresses have been made in this field. Covalent H3S hydrogen-rich superconductors (Tc = 200 K) and ionic hydrogen-rich superconductors with hydrogen-cage structure, such as LaH10 (Tc = 260 K, –13 ℃), YH6 and YH9, have been successively synthesized, setting a new record of superconducting temperature. These studies have given rise to the hope of discovering room-temperature superconductors in hydrogen-rich compounds under high pressure. This paper focuses on the progress of hydrogen-rich superconductors with high critical temperature under high pressure, discusses the physical mechanism of high-temperature superconductivity in hydrogen-rich compounds, provide an outlook on the possibility of discovering room-temperature superconductors in hydrogen-rich compounds in the future, and offer the candidate system for high superconductivity in multiple hydrogen-rich compounds.

## EDITOR'S SUGGESTION

2021, 70 (1): 017405. doi: 10.7498/aps.70.20201881
Abstract +
$\psi = {\psi _{\rm{0}}}{e^{i\varphi }}$, and the phase is uniform over the space. When applying an external field but still below a certain threshold, a screening current will be established at the surface, which prohibits the entering of magnetic field, that is so-called Meissner effect. When the external field is larger than this threshold, the magnetic flux will penetrate into the sample, forming the interface of superconducting and normal state regions. According to the sign of this interface energy, we can categorize superconductors into type-I (positive interface energy) and type-II (negative interface energy). Most superconductors found so far are type-II in nature. Due to the negative interface energy in type-II superconductors, the penetrated magnetic flux will separate into the smallest bundle, namely the quantum flux line, with a quantized flux ${\varPhi _0} = h/2e$ (h is the Planck constant and e is the charge of an electron). There are weak repulsive interactions among these vortices, thus usually they will form a lattice, called mixed state. When applying a current, a Lorentz force will exert on the flux lines (vortices) and will make them to move, this will induce energy dissipation and the appreciable feature of zero resistance of a superconductor will be lost. By introducing some defects, impurities or dislocations into the system, it is possible to pin down these vortices and restore the state of zero resistance. The study concerning vortex pinning and dynamics is very important, which helps not only the understanding of fundamental physics, but also to the high power application of type-II superconductors. This paper gives a brief introduction to the vortex dynamics of type-II superconductors.">Superconductivity is achieved through macroscopic phase coherence; the charge carriers are Cooper pairs. In absence of an external magnetic field and applied current, the behavior of these Cooper pairs can be described by a single wave function $\psi = {\psi _{\rm{0}}}{e^{i\varphi }}$, and the phase is uniform over the space. When applying an external field but still below a certain threshold, a screening current will be established at the surface, which prohibits the entering of magnetic field, that is so-called Meissner effect. When the external field is larger than this threshold, the magnetic flux will penetrate into the sample, forming the interface of superconducting and normal state regions. According to the sign of this interface energy, we can categorize superconductors into type-I (positive interface energy) and type-II (negative interface energy). Most superconductors found so far are type-II in nature. Due to the negative interface energy in type-II superconductors, the penetrated magnetic flux will separate into the smallest bundle, namely the quantum flux line, with a quantized flux ${\varPhi _0} = h/2e$ (h is the Planck constant and e is the charge of an electron). There are weak repulsive interactions among these vortices, thus usually they will form a lattice, called mixed state. When applying a current, a Lorentz force will exert on the flux lines (vortices) and will make them to move, this will induce energy dissipation and the appreciable feature of zero resistance of a superconductor will be lost. By introducing some defects, impurities or dislocations into the system, it is possible to pin down these vortices and restore the state of zero resistance. The study concerning vortex pinning and dynamics is very important, which helps not only the understanding of fundamental physics, but also to the high power application of type-II superconductors. This paper gives a brief introduction to the vortex dynamics of type-II superconductors.

## EDITOR'S SUGGESTION

2021, 70 (1): 018401. doi: 10.7498/aps.70.20202042
Abstract +
This paper presents a brief review of the development trend of superconducting magnets in large scale applications towards high magnetic fields, depending on and pushing the Nb3Sn wire technics' continuous improvement. The focus is on analysis of the technology challenges of 14 T whole-body superconducting magnets. Using the Bonze Nb3Sn wires and on the base of a combination design of Nb3Sn and NbTi coils, an electromagnetic conception design of a 14 T whole-body MRI magnet is presented, and the thermal stability and quench protection are analyzed by simulations. The critical issues on stress, joints as well as shimming of 14 T whole-body superconducting magnets are also discussed. According to the results, this paper believes: 1) Nb3Sn wires are of the first important issue for 14 T whole-body superconducting magnets—the Bonze Nb3Sn wire is of the best choice but the performance specifications of the current products need to be improved further to match the requirements; 2) quench protection of 14 T whole-body superconducting magnets is one of the most complicated technics that covers design of the copper to superconductor (Cu/SC) ratio, coordination of the operating current and coil inductances, subdivisions of passive protection circuits and quench triggering control of active protection, as well as the stray field limitation during the transient process.

## EDITOR'S SUGGESTION

2021, 70 (1): 018502. doi: 10.7498/aps.70.20202131
Abstract +
Superconductivity is a macroscopic quantum phenomenon. Flux quantization and the Josephson effect are two physical phenomena which can best reflect the macroscopic quantum properties. Superconducting quantum interference device (SQUID) is one type of superconducting devices which uses these two characteristics. SQUID devices are widely used in the sensitive detection of magnetic signals. This paper briefly introduces the background and recent developments of low temperature superconductor and high temperature superconductor SQUID devices.

## EDITOR'S SUGGESTION

2021, 70 (1): 018501. doi: 10.7498/aps.70.20202121
Abstract +
It has been nearly 110 years since the discovery of superconductors, and more than 30 years since the discovery of high temperature superconductors (HTS). Great progress has been made in the application of superconducting electronics in the last two decades. HTS microwave devices have shown much higher perfomance than the traditional ones and have found their ways to the industry applications in mobile communication, radar, and special communication applications. Owing to the ultrahigh sensitivity to magnetic fields and currents, superconducting quantum interference devices (SQUIDs) have been used as the irresplacible sensors in geological surveying, magnetic resonanc imaging, biomagnetic imaging, and other areas. The sensitivity of superconducting radiation detectors such as superconducting SIS mixer, superconducting hot electron bolometer, superconducting transition edge sensor, superconducting nanowire single photon detector, and superconducting microwave kinetic inductance detector are near the quantum limitation. They are now key technology in geophysics, astrophysics, quantum information science, biomedicine, and so on. Superconducting Josephson parametric amplifier has become a key element for superconducting quantum computing. Superconducting integrated circuit has been included in the international roadmap for devices and systems, and shows that having the potential to become one of the mainstreams for post-Moore information processing technology. In metrology, superconducting Josephson effect and Josephson junction array devices have been widely used in the redefinition of quantum voltage reference and basic units of the International system of Units. Superconducting electronics plays an important role in the current quantum information technology boom, which in turn promotes the development of superconducting electronics. This review will brief introduce the research and application of superconducting electronics in China in recent years.

## EDITOR'S SUGGESTION

2021, 70 (1): 010201. doi: 10.7498/aps.70.20201371
Abstract +
Background: The coronavirus disease 2019 (COVID-19) has raged more than 10 months and it has become a major public health concern. It is necessary to account for the intrinsic mechanisms and reveal the transmission pattern. Method: We collect detailed information of 944 COVID-19 cases in Guangdong province from January 23rd to February 16th. According to the age-structured characteristics, the population is divided into four groups such as child group (0–5 years old), adolescent group (6–19 years old), young and middle-aged group (20–64 years old), elderly group (65 and over years old). Coupling with different age-structured contact patterns, we establish a discrete age-structured COVID-19 model, obtain the basic reproduction number and final size. By Markov Chain Monte Carlo numerical method (MCMC), we identify the model parameters, fit the cumulative cases, calculate eradiation time of disease, infection peak and the peak arrival time, etc. Results: We found that the most infected people are the young and middle-aged individuals; Compared with household quarantine measure, the peak value of hospitalizations among young and middle-aged group in community mode will increase of 41%, and the peak will delay two weeks. By analyzing the proportions of the final sizes associated age groups, it is found that the elderly have a higher susceptibility, while the adolescents have a lower susceptibility. Under the household quarantine measure, if infected individuals have been confirmed in time of half a day, the peak size of hospitalizations will be further reduced, and the peak hospitalization will advance one week. The model reveals social contact patterns for impacting on COVID-19 transmission, and evaluates the effectiveness of household quarantine.
2021, 70 (1): 010202. doi: 10.7498/aps.70.20200825
Abstract +
It is a novel and interesting idea to inversely design the scattering structure with the desired scattering field intensity distribution in a given target area as the known information. The inverse design method proposed in this paper does not need to be optimized, and the spatial distribution and dielectric constant distribution of the micro-scatterer array can be quickly analytically calculated according to the desired scattering field intensity in the target area. First, based on the spatial Fourier transform and angular spectrum transformation, the plane wave sources required in all directions are inversely obtained from the electric field intensity distribution required in the target area. Then, based on the theory of induced source, a method of irradiating the array of all-dielectric micro-scatterers with incident electromagnetic field to generate the required plane wave source is proposed. The scattering fields generated by these micro-scatterers will be superimposed on the target area to achieve the desired scattering field strength intensity. Finally, according to the proposed inverse design theory model, a specific three-dimensional (3D) design is carried out. In the 3D example, we study the scattering field intensity distribution of the point-focused shape of the square surface target area, and show an all-dielectric micro-sphere distribution design. Its spatial distribution and permittivity distribution are both obtained through the rapid analytical calculation of the desired scattered field intensity shape in the target area. Finally, based on the principle of linear superposition, we quickly and easily generate the complex shapes of “I”, “T”, and “X” in the target area. The satisfactory results of full-wave simulation show that the proposed inverse design method is effective and feasible.
2021, 70 (1): 010301. doi: 10.7498/aps.70.20201199
Abstract +
The time evolution of multipartite quantum coherence is studied in a three-body spin system with an asymmetric interaction. The l1 norm measurement is used to estimate the degree of quantum coherence in the spin system. The decoherence of all components of quantum coherence in the three-body spin system is analyzed by the exact diagnolization and numerical method based on quantum master equation. The environmental noise induced by the spontaneous decay can be simulated by the quantum amplitude damping model. It is found that the time evolution of quantum coherence component is closely related to the quantum property of the whole initial state. When the initial state is a separable pure one, the asymmetric interactions can conduce to the enhancement of the degree of multipartite quantum coherence in a short time interval. Under the influence of quantum noise, the degree of quantum coherence decreases gradually. We pay much attention to the spatial distribution of the degree of quantum coherence in a many-body system. The additivity relationship of bipartite component and tripartite coherence can exist if the initial state is chosen to be a Werner-like state. This kind of the coherence additivity between all bipartite components and global coherence can be extended to an arbitrary N-body Werner-like state. But this additivity relationship depends on the l1 norm coherence measurement. Owing to the asymmetric interaction and noise, the degree of tripartite quantum coherence is more than the sum of all degrees of bipartite quantum coherence. The difference between the degree of tripartite coherence and the sum of all degrees of bipartite coherence is increased in a short time interval. The environmental noise can also suppress the difference in the coherence degree. The degree of the nearest neighboring bipartite coherence decreases more quickly than those of other bipartite coherences. The asymmetric interaction gives rise to the improvement in the degree of bipartite coherence and tripartite coherence. The coherence of the next-nearest neighboring two systems can be robust against the environmental noise. These results are helpful in preparing the multipartite quantum resources. We can utilize the system of coupled micro-cavities to realize the quantum spin system with controllable asymmetric interaction. In this way, the global coherence and bipartite coherence can be manipulated effectively by the quantum electromagnetic technology.
2021, 70 (1): 010501. doi: 10.7498/aps.70.20200899
Abstract +
Chaotic time series forecasting has been widely used in various domains, and the accurate predicting of the chaotic time series plays a critical role in many public events. Recently, various deep learning algorithms have been used to forecast chaotic time series and achieved good prediction performance. In order to improve the prediction accuracy of chaotic time series, a prediction model (Att-CNN-LSTM) is proposed based on hybrid neural network and attention mechanism. In this paper, the convolutional neural network (CNN) and long short-term memory (LSTM) are used to form a hybrid neural network. In addition, a attention model with softmax activation function is designed to extract the key features. Firstly, phase space reconstruction and data normalization are performed on a chaotic time series, then convolutional neural network (CNN) is used to extract the spatial features of the reconstructed phase space, then the features extracted by CNN are combined with the original chaotic time series, and in the long short-term memory network (LSTM) the combined vector is used to extract the temporal features. And then attention mechanism captures the key spatial-temporal features of chaotic time series. Finally, the prediction results are computed by using spatial-temporal features. To verify the prediction performance of the proposed hybrid model, it is used to predict the Logistic, Lorenz and sunspot chaotic time series. Four kinds of error criteria and model running times are used to evaluate the performance of predictive model. The proposed model is compared with hybrid CNN-LSTM model, the single CNN and LSTM network model and least squares support vector machine(LSSVM), and the experimental results show that the proposed hybrid model has a higher prediction accuracy.
2021, 70 (1): 010701. doi: 10.7498/aps.70.20201228
Abstract +
2021, 70 (1): 010702. doi: 10.7498/aps.70.20201066
Abstract +
NO3 radical is the most important oxidant in atmospheric chemistry at night, and it controls the oxidation and removal of various trace gas components in the atmosphere. The understanding of the chemical process of NO3 radical is of great significance for studying the atmospheric pollution processes such as haze. The NO3 radical has a low concentration and strong activity, so it is relatively difficult to measure accurately. We report here in this paper an instrument for unambiguously measuring NO3 based on broadband cavity enhanced absorption spectroscopy (BBCEAS). To achieve the robust performance and system stability under diverse conditions, this BBCEAS instrument has been developed, with efficient sampling, and resistance against vibration and temperature change improved, and the BBCEAS instrument also has low-power consumption. The 660-nm-wavelemngth light-emitting diode (LED) is used as a light source of the BBCEAS system. The sampling gas path with low loss and suitable for domestic high-particle environment is designed. Through the LED light source test, the optimal working current and temperature can be obtained to achieve the acquisition of NO3 absorption spectrum with high signal-to-noise ratio. Considering the fact that the water vapor absorption is an important interference factor for the measurement of NO3 radical by BBCEAS, the daytime atmospheric measurement spectrum is used as a background spectrum, and participates in spectral fitting of NO3 to reduce the effect of water vapor. The mirror reflectivity and effective cavity length are calibrated, and the Allan variance analysis is also carried out. The reflectance of the mirror can reach about 0.99993 at 662 nm (NO3 absorption peak), and the corresponding theoretical effective optical path can reach more than 7 km, which can meet the measurement requirements of atmospheric NO3 radicals. The detection limit (1σ) of 0.75 pptv for NO3 is achieved with an acquisition time of 10 s and a total measurement error of about 16%. The atmospheric NO3 radical observation is carried out in Hefei. During the observation period, the highest NO3 concentration is 23.4 pptv, demonstrating the promising potential applications in in-situ, sensitive, accurate and fast simultaneous measurements of NO3 in the future by using the developed broadband cavity enhanced absorption spectroscopy.
###### NUCLEAR PHYSICS
2021, 70 (1): 012901. doi: 10.7498/aps.70.20200830
Abstract +
High-energy colliders play an indispensable role in particle physics and high-energy physics. Beam screen is one of the key parts in the high-energy collider. It is used to transfer the heat generated by the beam in the pipeline to a cooling system, and absorb the residual gas to the cold bore through the pumping holes on the wall of the beam screen to ensure the vacuum stability at the same time. However, in the process of transferring thermal load, the deformation caused by temperature change will affect the structural stability of the beam screen. How to reduce the deformation as much as possible while ensuring the good heat transfer performance of the beam screen is one of the key issues in optimizing the structural design of the beam screen. In this paper, the heat transfer performance and mechanical property of the beam screen model are simulated and optimized based on the ANSYS simulation results to ensure the normal and stable operation of the beam in the super proton-proton collider. For the inner surface of the outer screen of the beam screen, the method of reducing the thickness of the copper coating is used to reduce the Lorentz force generated during operation. The calculation results from the relevant theoretical models show that when the thickness of the copper coating varies from 0 to 100 μm, the copper coating with a thickness of 75 μm can reduce the maximum deformation of the outer screen of the beam screen by 70.9%, while the maximum temperature of the beam screen can be increased by 1.1%. For the inner screen of the beam screen, a design scheme in which supporting ribs are arranged at intervals is used to reinforce the structure and improve the overall structural stability of the beam screen. The calculation results show that the maximum deformation of the inner screen of the beam screen can be reduced by 86.8% and the maximum temperature of the beam screen is reduced by 7.69%, compared with the case without supporting fins, when the interval between two adjacent supporting fins is 1 pumping hole. The research results provide important theoretical reference for the design of beam screen, which is the key component of the vacuum system of the new-generation high energy particle accelerator.
###### ELECTROMAGNETISM, OPTICS, ACOUSTICS, HEAT TRANSFER, CLASSICAL MECHANICS, AND FLUID DYNAMICS
2021, 70 (1): 014101. doi: 10.7498/aps.70.20201104
Abstract +
Microstrip phased array has aroused interest of many researchers because of its beam agility. However, a big problem for typical microstrip array is that its main beam can only scan from about –50° to 50°, with a gain loss of 4-5 dB. Meanwhile, the relatively narrow operating bandwidth of microstrip antenna is also a problem in application. These flaws have dramatically limited its applications and spawned many studies on phased array with wide-angle scanning capability. Several methods have been proposed to broaden the scanning coverage of phased array, such as utilizing pattern-reconfigurable antenna as an element of array, taking wide-beam antenna as the element of array, and adopt metasurface as the top cladding of array. However, most of existing researches mainly focus on achieving wide-angle scanning performance within a relatively narrow bandwidth. A phased array that possesses wide-angle scanning capability at both main planes within a relatively wide bandwidth is highly desirable. In this paper, a wide-beam magnetoelectric (ME) dipole antenna is proposed. It consists of an ME dipole antenna in the form of microstrip patch and a pair of magnetic dipoles. Metallic through holes integrated with patches and ground are utilized to form magnetic currents. Extra magnetic dipoles are added to broaden the 3-dB beam-width. The simulated results reveal that the 3-dB beam-width of the proposed antenna is greater than 107° in the E-plane (9 GHz–12 GHz) and 178° in the H-plane (7 GHz–12 GHz) respectively. The impedance bandwidth of the proposed antenna is 53.26% from 7.3 GHz to 12.6 GHz (VWSR < 2). Based on the proposed antenna element, two linear phased arrays are fabricated and measured. To test the wide-angle scanning capability of the arrays, each antenna element is simply fed with alternating currents with identical amplitude and linearly increasing phases. The measured results reveal that the wide-angle scanning capability of H-plane array and E-plane array can be obtained from 9 GHz to 12 GHz. The scanning beam of the H-plane array can cover the range from -90° to 90°. The scanning beam of the E-plane array can cover the range from –70° to 70°. The impedance bandwidth of the central antenna is 27.03% for the H-plane array from 9.6 GHz to 12.6 GHz (active VWSR < 2.5) and 36.36% for the E-plane array from 9 GHz to 13 GHz (active VWSR < 2) respectively. Hence, the proposed method can be used as a reference for designing a wide-beam antenna and wide-angle scanning phased array and the designed phased arrays can be applied to X-band radar systems.
2021, 70 (1): 014102. doi: 10.7498/aps.70.20201089
Abstract +
The reflection and transmission of plane electromagnetic waves on monolayer graphene are studied theoretically in this paper. From an electromagnetic point of view, monolayer graphene is described as an “infinitely thin” graphene sheet characterized by a surface conductivity, and based on a microscopic quantum dynamical approach, the graphene sheet becomes anisotropic in the presence of both an electrostatic and a magnetic bias. In this work, starting from boundary conditions and phase-matching conditions, the propagation matrix for the analysis of the interaction between an electromagnetic field and thin graphene sheet which is biased electrostatically and magnetostatically, and then characterized by an anisotropic conductivity, is derived. Furthermore, the analytical solutions of co- and cross-polarization reflective and transmittance coefficients through an anisotropic graphene planar surface are obtained from the proposal matrix above, which couples the fundamental transverse electric (TE) polarization and transverse magnetic (TM) polarization and includes the possible effects of electrostatic and/or magnetostatic bias. In conclusion, the cross-polarization reflective coefficient of TE wave and that of TM wave are equal, and their cross-polarization transmittance coefficients have opposite phase. Finally, a new propagation matrix for stratified medium containing anisotropic graphene interfaces is deduced by embedding the matrix across graphene sheet mentioned above into the traditional propagation matrix for isotropic stratified medium. The proposed new matrix can be used to investigate the propagation properties of plane wave in a complex structure of layered medium and anisotropic conductivity interfaces (including graphene sheet) analytically and quickly, and represents a very simple tool for the relevant analysis and design.
2021, 70 (1): 014201. doi: 10.7498/aps.70.20201017
Abstract +
$T_{\rm TiO_2} = 20$ nm are chosen as the guideline for designing the dye laser, which generates the resonance wavelength of 820 nm the same as the absorption wavelength of dye molecules. Additionally, the laser characteristics excited by pump light with the wavelength of 820 nm are studied. The continuous laser output is obtained. The energy threshold is about 2.5 mJ/cm2, and the linewidth is about 0.3 nm. The proposed structure can effectively regulate the surface local electric field and enhance the interaction between pump light and gain medium. It can not only be used in lasers, but also provide reference for designing other light-emitting devices.">The enhancement of lasing emission intensity of dye laser is particularly important and urgently required due to a broad range of optical and electrical applications. The guided-mode resonance (GMR) effect occurs in a periodic waveguide structure where an incident wave is coupled to a leaky waveguide mode, and yields a resonance peak. The resonance wavelength can be easily controlled by adjusting the period of the grating, thickness of the waveguide layer, and refractive index of the covering materials. By using band edge states, one may be able to excite optical resonances extended over the entire structure surface, thereby achieving field enhancement over a large area. In this study, mesoporous silica with low refractive index is introduced between the grating layer and the substrate layer of the GMR structure to significantly enhance the contact between local electric field and gain medium. For comparison, another structure using SU-8 with high refractive index as the spacing layer is also proposed. It is clearly observed that the maximum of the electric field intensity is highly localized inside the SU-8 waveguide grating layer. However, it is shifted upward to the gain medium layer in the mesoporous silica structure due to the reverse symmetry waveguide structure design. Therefore, the interaction between laser dye and electric field is increased to further enhance the lasing emission.Besides the refractive index, the waveguide layer, other structural parameters such as thickness of each layer and grating period also affect the electric field distribution in the GMR structure. Based on the finite-difference time-domain method, the structural parameters are analyzed and optimized. According to the simulation results, the structure parameters TWG = 3.5 μm, Λ = 700 nm, and $T_{\rm TiO_2} = 20$ nm are chosen as the guideline for designing the dye laser, which generates the resonance wavelength of 820 nm the same as the absorption wavelength of dye molecules. Additionally, the laser characteristics excited by pump light with the wavelength of 820 nm are studied. The continuous laser output is obtained. The energy threshold is about 2.5 mJ/cm2, and the linewidth is about 0.3 nm. The proposed structure can effectively regulate the surface local electric field and enhance the interaction between pump light and gain medium. It can not only be used in lasers, but also provide reference for designing other light-emitting devices.
2021, 70 (1): 014202. doi: 10.7498/aps.70.20200721
Abstract +
An ultra-compact 1×2 demultiplexer based on directional coupler (DC) waveguide is proposed to separate the 1310 nm wavelength from 1550 nm wavelength, in which a new Si3N4/SiNx/Si3N4 sandwiched structure is used to realize polarization insensitivity. Firstly, the new sandwiched structure is designed to be polarization-independent. The coupling lengths of two orthogonal polarization modes at the same wavelength versus the gap between two parallel SiNx waveguides g1 are calculated with several groups of structure parameters of the demultiplexer. The result shows that the coupling lengths for the two orthogonal polarization modes at the same wavelength can be identical by choosing the proper g1. Then, how to realize the function of wavelength separation is studied. When one wavelength propagates at even multiple of coupling length and the other wavelength propagates at odd multiple of coupling length, and vice versa, the two working wavelengths will output from different output ports, thereby the two wavelengths are successfully separated. Under the premise of satisfying such conditions, a comparison of size and performance among the devices with different groups of structure parameters is given to find the best one. The demultiplexer based on Si3N4/SiO2 platform has a compact structure, easy integration and good tolerance. Three-dimensional(3D) finite-difference time-domain method is used for simulation, and the results show that the length of the DC waveguide is only 23 μm; the insertion loss and crosstalk are as low as 0.1 dB and–26.23 dB respectively; a broad 3-dB bandwidth of 200 nm is achieved. To demonstrate the transmission characteristics of the demultiplexer, the evolution of the excited fundamental mode in the demultiplexer is also given. The novel demultiplexer is polarization-independent and can work at 1310 nm and 1550 nm wavelengths simultaneously. It has a potential application value in future integrated optical circuits.
2021, 70 (1): 014301. doi: 10.7498/aps.70.20200973
Abstract +
The combination of superparamagnetic iron oxide nanoparticles (SPIOs) with ultrasonic contrast agent (UCA) microbubble is called magnetic microbubble (MMB) and has been used to produce multimodal contrast agents to enhance medical ultrasound and magnetic resonance imaging. The nanoparticles are either covalently linked to the shell or physically entrapped into the shell. Considering the effect of the volume fraction of SPIOs on the shell density and viscosity, a nonlinear dynamic equation of magnetic microbubbles (MMBs) with multilayer membrane structure is constructed based on the basic theory of bubble dynamics. The influences of the driving sound pressure and frequency, particle volume fraction, shell thickness and surface tension on the acoustic-dynamics behavior of microbubbles are numerically analyzed. The results show that when the volume fraction of magnetic particles is small and α ≤ 0.1, the acoustic properties of magnetic microbubbles are similar to those of ordinary UCA microbubbles. The acoustic response of the microbubble depends on its initial size and driving pressure. The critical sound pressure of microbubble vibration instability is lowest when the driving sound field frequency is twice the magnetic microbubble resonance frequency f0 (f = 2f0). The presence of magnetic particles inhibits the bubbles from expanding and contracting, but the inhibition effect is very limited. The surface tension parameter K of the outer film material and thickness of the shell also affect the vibration of the microbubble. When K and film thickness are 0.2–0.4 N/m and 50–150 nm respectively, it is observed that the bubble has an unstable vibration response region.
###### PHYSICS OF GASES, PLASMAS, AND ELECTRIC DISCHARGES
2021, 70 (1): 015201. doi: 10.7498/aps.70.20200794
Abstract +
$q$. The influence of parameter $q$ on the criterion is discussed in this paper. The lower limit of the ion Mach number changes with the value of parameter $q$. The lower limit of the ion Mach number increases for $q < 1$. And the lower limit of the ion Mach number decreases for $q>1$. With the increase of $q$, the number of electrons with lower speed increases, ions need less kinetic energy to enter into the sheath and thus enter into the sheath more easily. Through numerical simulation, it is found that compared with the structure of the plasma magnetized sheath with Maxwell distribution ($q=1$), the structure of the plasma magnetized sheath with super-extensive distribution ($q < 1$) and that with sub-extensive ($q>1$) are different, including the distribution of the space potential, the ion density, the electron density, and the space charge density. When $q < 1$, the space potential, the electron density and the ion density fall more slowly, and the peak of the space charge density curve is closer to the wall. When $q>1$, the space potential and the ion electron density fall faster, especially the electron density drops to zero faster, and the peak of the space charge density curve is far away from the wall. The simulation results show that the non-extensive parameter $q$ has a significant influence on the structure of collisional plasma magnetized sheath. The influence of the collision on the magnetized plasma sheath with non-extensive distribution is similar to that with the Maxwell distribution. These conclusions may be useful in solving the problems of plasma boundary.">Many previous researches on the plasma sheath were based on the fact that the electrons satisfy the classical Maxwell velocity distribution function, while the particles in the plasma have long-range electromagnetic interactions. It is more appropriate to use the non-extensive distribution proposed by Tsallis to describe the electrons. In this paper, a collisional magnetized plasma sheath model with non-extensive distribution of electrons is established. Bohm criterion is derived theoretically. With the ion drift motion in the plasma pre-sheath region taken into consideration, the ion Mach number is only related to the angle of the magnetic field, the collision parameters, the electric field at the sheath edge, and non-extensive parameter $q$. The influence of parameter $q$ on the criterion is discussed in this paper. The lower limit of the ion Mach number changes with the value of parameter $q$. The lower limit of the ion Mach number increases for $q < 1$. And the lower limit of the ion Mach number decreases for $q>1$. With the increase of $q$, the number of electrons with lower speed increases, ions need less kinetic energy to enter into the sheath and thus enter into the sheath more easily. Through numerical simulation, it is found that compared with the structure of the plasma magnetized sheath with Maxwell distribution ($q=1$), the structure of the plasma magnetized sheath with super-extensive distribution ($q < 1$) and that with sub-extensive ($q>1$) are different, including the distribution of the space potential, the ion density, the electron density, and the space charge density. When $q < 1$, the space potential, the electron density and the ion density fall more slowly, and the peak of the space charge density curve is closer to the wall. When $q>1$, the space potential and the ion electron density fall faster, especially the electron density drops to zero faster, and the peak of the space charge density curve is far away from the wall. The simulation results show that the non-extensive parameter $q$ has a significant influence on the structure of collisional plasma magnetized sheath. The influence of the collision on the magnetized plasma sheath with non-extensive distribution is similar to that with the Maxwell distribution. These conclusions may be useful in solving the problems of plasma boundary.
2021, 70 (1): 015202. doi: 10.7498/aps.70.20201206
Abstract +
Compared with the two-electrode gas spark switch, the three-electrode gas spark switch has high controllability, low working voltage and small jitter, so the three-electrode gas spark switch is widely used in pulse power technology. The discharge of gas spark switch is high pressure gas discharge, which is characterized by high electron collision frequency (1012 Hz), small mean free path (10–6 m), short breakdown time (10–9 s), and complex physical process (including the secondary electron emission, the generation of seed electrons, the space charge effect and various collision processes between electrons and nitrogen molecules, etc). At present, it is difficult to quantitatively describe the breakdown process of the three-electrode gas switch, and the detailed theoretical research is lacking. Therefore, the breakdown mechanism of atmospheric pressure nitrogen spark switch, including two-electrode and three-electrode, is studied theoretically and numerically in this paper. The purpose of this study is to compare the simulation results of the two different gas spark switches, and obtain the characteristics of stream breakdown in different gas spark switches. Firstly, the numerical simulation and theoretical analysis of two-electrode gas spark switch are carried out. According to theoretical and numerical calculation, it can be found that for the plate-plate two-electrode switch, the stream breakdown cannot be generated under low voltage (less than 6.3 kV), while under high voltage (more than 6.3 kV), first the anode-directed streamer is formed, and then the cathode-directed streamer is created. In addition, the simulation results show that the plasma generated by the trigger can effectively reduce the breakdown voltage. Finally, the three-electrode gas spark switch is studied theoretically and numerically. It can be seen that in the breakdown process of the three-electrode gas spark switch, the breakdown first occurs between the trigger and the insulator, and then this plasma channel expands to the anode and cathode, finally forming the arc channel between the anode and the cathode. Under the calculation conditions in this paper, if the cathode-trigger and the anode-trigger are required to break down simultaneously, the applied voltage between the cathode-trigger should be greater than 1.18 kV, while the applied voltage between the anode-trigger should be greater than 3 KV. When the field emission of the trigger is considered, the breakdown threshold can be significantly reduced.
###### CONDENSED MATTER: STRUCTURAL, MECHANICAL, AND THERMAL PROPERTIES
2021, 70 (1): 016801. doi: 10.7498/aps.70.20201008
Abstract +
###### CONDENSED MATTER: ELECTRONIC STRUCTURE, ELECTRICAL, MAGNETIC, AND OPTICAL PROPERTIES
2021, 70 (1): 017701. doi: 10.7498/aps.70.20201352
Abstract +
2021, 70 (1): 018701. doi: 10.7498/aps.70.20200986
Abstract +
A algorithm of obtaining absolute dose at each tissue depth only by the mathematical calculation of formula is reported. The algorithm is based on integrating the energy flux of the electron generated by X-ray in the range of irradiation field, and the energy spectrum of ray and the secondary scattered ray are considered in this process. In this algorithm, the water phantom in the irradiation field is divided into several thin layers, and the energy flux of the electrons generated by interaction between the ray and thin layer reaching the calculation point is calculated. Finally, the absolute dose of the calculation point can be obtained by accumulating the energy flux contribution of all thin layers. For the X-ray with continuous energy spectrum, the expected mass attenuation coefficient is calculated for obtaining the photon flux at each depth in this process. The absolute dose calculated by this algorithm is verified by Monte Carlo simulation, and the difference between the algorithm and simulation is compensated for by a dose function about multiple scattering photons, and the function shows fast descent and then slow ascent. It is found that the ratio of the dose caused by backscatter to the surface dose, and the relationship among forward scatter, backward scatter and primary ray, and the relationship between the dose and the depth of the secondary scattered rays show a trend of first rising and then declining, and the depth of the peak value deviates from the position of the thin layer. Three-dimensional energy spectra of the secondary photon and the secondary electron are also compared with each other, and the spectrum is a function of particle flux about particle energy and particle direction. From the perspective of Compton effect, the physical meanings of different positions in the three-dimensional energy spectrum of the two particles are explained. It is found that the difference between algorithm percentage depth dose and simulation percentage depth dose is similar to the difference between small irradiation field percentage depth dose and big irradiation field percentage depth dose from simulation, and it is verified that the difference between algorithm and simulation comes from the increase of scattered rays. Finally, the algorithm is applied to the dose calculation of non-uniform phantom, which can accurately reflect the dose distribution characteristics and have less error.
2021, 70 (1): 018702. doi: 10.7498/aps.70.20200904
Abstract +
Epilepsy is an extensive nervous system disease nowadays. Electroencephalogram (EEG) can capture the abnormal discharge of nerves in the brain duration of seizure and provide a non-invasive way to identify epileptogenic sites in the brain. In order to distinguish between focal epilepsy EEG signal and non-focal epilepsy EEG signal, in this paper we propose an automated epileptic EEG detection method based on the elastic variational mode decomposition (EVMD). The proposed EVMD algorithm is a method of analyzing the signals and also a processing method in time-frequency domain, in which the elastic net regression is used to reconstruct a constrained variational model in variational mode decomposition (VMD). Used in the VMD algorithm is the Tikhonov regularization that is also statistically called ridge regression as a solution of recovering the unknown signal and assessing the bandwidth of a mode, namely the variational equation constructed by VMD only has L2 norm. However, the ridge regression cannot select variables when the equation has multiple variables. Another regression method, called lasso regression, only has L1 norm and can select a more accurate model from multiple variables, but it has worse performance when variables have group effect or co-linearity. The elastic net regression has advantages of ridge regression and lasso regression, in other word, the variational equation constructed by EVMD has both L1 regularization item and L2 regularization item, so in this paper we propose the EVMD by elastic net regression. In addition, in this paper the EVMD is used to distinguish between focal epilepsy EEG signal and non-focal epilepsy EEG signal. Firstly, the original EEG signals are divided into several sub-signals where the test signals are divided into sub-signals with shorter durations by time series and a reasonable time overlap is kept between successive sub-signals. After that each sub-signal is decomposed into intrinsic mode functions by using the EVMD. Furthermore, the refined composite multiscale dispersion entropy (RCMDE) as feature is extracted from each intrinsic mode function where a Student’s t-test is used to assess the statistical differences between RCMDEs extracted from focal and non-focal EEG signals respectively. Finally, the support vector machine (SVM) is used to classify their features. For an epilepsy EEG signalspublic data set, the final experimental results show that the performance indices of accuracy, sensitivity, and specificity can reach 92.54%, 93.22% and 91.86% respectively.

## EDITOR'S SUGGESTION

2021, 70 (1): 018801. doi: 10.7498/aps.70.20201194
Abstract +
Ag2ZnSnSe4 is an n-type semiconductor with a suitable bandgap of 1.4 eV. In the present study, a graphene/Ag2ZnSnSe4 induced p-n junction thin film solar cell is proposed and the physical mechanism and performance influencing factors of the solar cell are simulated and analyzed by using the wxAMPS software. The simulation results show that when a high work function graphene contacts an n-type Ag2ZnSnSe4 semiconductor, the energy band of the Ag2ZnSnSe4 absorber layer bends upward, meanwhile a p-type Ag2ZnSnSe4 inversion layer is induced on the surface of n-type Ag2ZnSnSe4, therefore the p-type Ag2ZnSnSe4 and n-type Ag2ZnSnSe4 form an induced p-n homojunction. It is found that the work function of graphene and back contact significantly influence the photogenerated carrier separation, transportation and collection. The graphene work function should be 5.5 eV and the work function of back contact should not be greater than 4.4 eV, which is beneficial to the improving of the device performance. The doping concentration of Ag2ZnSnSe4 absorber mainly affects the short-circuit current of the device, however, the defect density of Ag2ZnSnSe4 absorber affects the whole device performance. When the work function of graphene and back contact are 5.5 eV and 3.8 eV, the doping concentration and defect density of Ag2ZnSnSe4 absorber are 1016 cm–3 and 1014 cm–33, respectively, the conversion efficiency of the graphene/Ag2ZnSnSe4 induced p-n junction thin film solar cell can reach 23.42%. These simulation results provide the idea and physical explanation for designing a novel type of solar cell with high efficiency and low cost.
2021, 70 (1): 018901. doi: 10.7498/aps.70.20201065
Abstract +
$p=0.004$, $p=0.006$, $p=0.008$ and $p=0.01$. The results show that the node hyper degree distribution of this hyper network model complies to the Poisson distribution $p(k)\approx \dfrac{{{\left\langle \lambda \right\rangle }^{k}}}{k!}{{e}^{-\left\langle \lambda \right\rangle }}$, which conforms with the characteristics of random networks and is consistent with the theoretical derivation. Further, in order to more accurately and effectively describe the multiple heterogeneous relationship in real life, in this paper we construct three different kinds of double-layer hyper network models with node hyper degree distribution with bimodal peak characteristics. The three kinds respectively are ER-ER, BA-BA and BA-ER, where ER represents the ER random hyper network, and BA denotes the scale-free hyper network, and the layers are connected by a random manner. The analytical expressions of node hyper degree distribution of the three kinds of double-layer hyper network models are obtained by theoretical analysis, and the average node hyper degrees of the three double-layer hyper networks are closely related to the inter-layer hyper edge probability. As the inter-layer hyper edge probability increases, the average node hyper degree increases. The results of simulation experiments show that the node hyper degree distributions of three kinds of double-layer hyper network models proposed in this paper possess the characteristics of bimodal peaks. The ER random hyper network model and the double-layer hyper network model proposed in this paper provide the theories for further studying the hyper network entropy, hyper network dynamics, hyper network representation learning, hyper network link prediction, and traffic hyper network optimization of such hyper networks in the future, and also it has certain reference significance for studying the evolution of multilayer hyper networks.">With the rapid development of social economy, the relationship between social members and groups has shown more complex and diverse characteristics. As a network depicting complex relation and multi-layer, hyper network has been widely used in different fields. Random network that obeys Poisson distribution is one of the pioneering models studying complex networks. In the existing hyper network researches, the hyper network based on ER random graph is still a blank. In this paper, we first propose an ER random hyper network model which is based on the hypergraph structure and it adopts the ER random graph theory. Furthermore, using this model, the node hyper degree distribution of this hyper network model is analyzed theoretically, and the node hyper degree distribution is simulated under different hyper edge probabilities: $p=0.004$, $p=0.006$, $p=0.008$ and $p=0.01$. The results show that the node hyper degree distribution of this hyper network model complies to the Poisson distribution $p(k)\approx \dfrac{{{\left\langle \lambda \right\rangle }^{k}}}{k!}{{e}^{-\left\langle \lambda \right\rangle }}$, which conforms with the characteristics of random networks and is consistent with the theoretical derivation. Further, in order to more accurately and effectively describe the multiple heterogeneous relationship in real life, in this paper we construct three different kinds of double-layer hyper network models with node hyper degree distribution with bimodal peak characteristics. The three kinds respectively are ER-ER, BA-BA and BA-ER, where ER represents the ER random hyper network, and BA denotes the scale-free hyper network, and the layers are connected by a random manner. The analytical expressions of node hyper degree distribution of the three kinds of double-layer hyper network models are obtained by theoretical analysis, and the average node hyper degrees of the three double-layer hyper networks are closely related to the inter-layer hyper edge probability. As the inter-layer hyper edge probability increases, the average node hyper degree increases. The results of simulation experiments show that the node hyper degree distributions of three kinds of double-layer hyper network models proposed in this paper possess the characteristics of bimodal peaks. The ER random hyper network model and the double-layer hyper network model proposed in this paper provide the theories for further studying the hyper network entropy, hyper network dynamics, hyper network representation learning, hyper network link prediction, and traffic hyper network optimization of such hyper networks in the future, and also it has certain reference significance for studying the evolution of multilayer hyper networks.