Search

Article

x

Vol. 74, No. 13 (2025)

2025-07-05
SPECIAL TOPIC—Order tuning in disordered alloys
Research progress of ordered regulation engineering for developing high-frequency amorphous-based soft magnetic materials
DING Huaping, LIU Lichen, SHAO Liliang, ZHOU Jing, ZUO Dingrong, KE Haibo, WANG Weihua
2025, 74 (13): 136101. doi: 10.7498/aps.74.20250585
Abstract +
The rapid advancement of modern electronics, telecommunications, and artificial intelligence has driven an urgent demand for high-performance soft magnetic materials, particularly those compatible with third-generation semiconductors. These semiconductors, characterized by wide bandgaps, high breakdown fields, and superior thermal conductivity, enable power devices to operate at higher frequencies (> 1 MHz) and power densities. However, traditional soft magnetic materials, such as silicon steels and ferrites, face inherent trade-offs between critical properties: saturation magnetization (Bs) versus coercivity (Hc), permeability versus core loss, and mechanical strength versus magnetic “softness”. These limitations hinder their applications in emerging high-frequency high-efficiency scenarios. Amorphous soft magnetic materials, with their unique hierarchical ordered structures ranging from atomic scale to nano scale, offer a revolutionary platform to overcome these trade-offs. These materials exhibit rich physical properties governed by short-range order (SRO, < 0.5 nm), medium-range order (MRO, 0.5–2.0 nm), and amorphous-nanocrystalline dual-phase architectures. The concept of order modulation strategically tailoring the intrinsic characteristics (e.g., cluster density, topological configuration) and spatial arrangements of these ordered structures has emerged as a transformative approach to decoupling conflicting material properties. This review systematically examines the following key aspects: 1) Historical evolution of soft magnetic materials The development of soft magnetic materials has kept pace with advances in power electronics technology, from early silicon steels and ferrites to modern amorphous and nanocrystalline alloys. The advent of Fe-based amorphous alloys and finemet-type nanocrystalline alloys marks milestones in achieving high Bs (>1.6 T), ultra-low Hc (< 1 A/m), and reduces core losses at high frequencies. However, performance bottlenecks still exist near theoretical limits, and require innovative strategies. 2) Theoretical foundations of order modulation Order parameter theory: Landau’s phase transition theory and synergetics elucidate how magnetic order parameters govern macroscopic properties. In amorphous alloys, magnetic interactions are dominated by SRO clusters and their MRO arrangements.Magnetism-structure relationships: advanced techniques, such as atomic electron tomography (AET) and synchrotron pair distribution function (PDF) analysis, reveal that SRO/MRO structures directly influence exchange coupling, magnetic anisotropy, and domain wall dynamics. For instance, Fe-M (M = Si, B) clusters with dense packing enhance Bs, while MRO homogenization reduces Hc. 3) Advances in order-modulated amorphous soft magnetic materials Atomic-scale modulation: elemental doping (e.g., Co, Mo, Cu) and energy-field treatments (e.g., magnetic annealing and ultrasonic vibration) optimize local atomic configurations. For example, ultrasonic processing of Fe78Si9B13 ribbons induces stress relaxation, forming 2–3 nm Fe-M clusters that increase Bs to 183.2 emu/g while maintaining Hc at 4.2 A/m.Nanoscale dual-phase design: controlled crystallization of α-Fe(Si) nanocrystals (<15 nm) within an amorphous matrix creates exchange-coupled nanocomposites. Co-Mo co-doping in FeSiBCuNb alloys refines grain size to 11.8 nm, achieving a permeability of 65000 H/m at 100 kHz–44% higher than traditional finemet alloys.Interface engineering in soft magnetic composites (SMCs): core-shell architectures (e.g., FeSiB@FeB nanoparticles) with stress-buffering interfaces reduce eddy current losses while preserving permeability. Cold sintering of vortex-domain FeSiAl powders enables GHz-range operation with stable permeability (μi = 13 at 1 GHz). 4) Future directions and challenges Machine learning-driven design: integrating high-throughput simulations with AI models (e.g., XGBoost, random forests) accelerates the discovery of optimal compositions and order parameters. Recent work predicts Bs by using Fe content, mixing enthalpy, and electronegativity differences, guiding the synthesis of (Fe82Co18)85.5Ni1.5B9P3C1 alloys with Bs = 1.92 T.The new magnetic topology structure: the magnetic vortex structures and skyrmion-like configurations in ultrafine powders show the prospect of ultra-high-frequency applications (>100 MHz).Low-stress manufacturing: innovations such as ultrasonic rheoforming reduce compaction pressures by 99% (to 6.2 MPa), alleviating residual stress and improving SMC performance.In situ characterization: neutron scattering and grating-based imaging techniques enable real-time observation of domain dynamics under operational conditions (e.g., stress, magnetic fields).In conclusion, ordered modulation represents a paradigm shift in the design of soft magnetic material, linking atomic-scale interactions with macroscopic performance. By using multi-scale ordered structures and advanced manufacturing technologies, the next-generation amorphous-based materials are expected to revolutionize high-frequency power electronics, electric vehicles, and AI-driven systems. However, challenges in scalable production, cost-effective processing, and standardized evaluation must be addressed to accelerate industrial adoption.
Coupling mechanism between high-temperature rheological behavior and dynamic relaxation in metallic glasses
LIANG Shuyi, ZHANG Langting, ZHU Hangchen, XING Guanghui, QIAO Jichao
2025, 74 (13): 136401. doi: 10.7498/aps.74.20250392
Abstract +
This study aims to establish the intrinsic link between the high-temperature rheological behavior and kinetic relaxation characteristics of La-based metallic glasses. By conducting dynamic mechanical analysis and high-temperature tensile strain-rate jump experiments on three La-based metallic glasses with significant β relaxation, and combining the findings within the free volume theory framework, their high-temperature rheological properties are investigated systematically. The results show that the steady-state flow stress and activation volume evolution trend are consistent within the normalized temperature range. The average activation energy for high-temperature rheology aligns with the activation energy range of α relaxation, confirming the strong association between rheological behavior and α relaxation. The activation energy for β relaxation shows an opposite trend, indicating that it may precede α relaxation. A dynamic competition between defect annihilation and generation governs the rheological behavior, and kinetic parameters reveal the temperature and strain-rate sensitivity of metallic glasses. This study lays a theoretical foundation for optimizing the high-temperature mechanical properties of La-based metallic glasses and also provides new insights into understanding the coupling relationship between multi-scale relaxation behavior and rheological mechanisms in metallic glasses.
Optimizing microstructure and mechanical properties of CoCrFeNi high-entropy alloy microfibers by electric current treatment
BO Le, GAO Xiaoyu, NING Zhiliang, WANG Li, SUN Jianfei, ZHANG Zhenjiang, HUANG Yongjiang
2025, 74 (13): 138102. doi: 10.7498/aps.74.20250518
Abstract +
High-entropy alloy (HEA) microfibers exhibit promising prospects in microscale high-tech applications due to their exceptional mechanical properties and stability. However, the strength-plasticity tradeoff largely hinders their further industrial applications. Heat treatment can optimize the mechanical properties of HEA microfibers. However, the traditional heat treatment (CHT) faces challenges in accurately adjusting the microstructures in a short period of time, while also being prone to grain coarsening, which can affect performance. In this study, an electric current treatment (ECT) technique is used to finely modulate the properties of cold-drawn CoCrFeNi high-entropy alloy microfibers on a microscale (~70 μm in diameter), the effects of thermal and athermal effects during ECT on microstructure and mechanical properties are systematically investigated through electron back scatter diffraction, transmission electron microscopy, and synchrotron radiation. A model of recrystallization, nucleation and growth of HEA microfibers is established. Compared with CHT, the synergistic effects of electron wind force and Joule heating during ECT significantly accelerate recrystallization kinetics, yielding finer and more homogeneous grains with a great decrease in dislocation density, and finally lead to better mechanical properties. The ECT-processed HEA microfibers achieve a yield strength in a range from 400 to 2033 MPa and a tensile elongation reaching 53%, which are much higher than those of CHT samples. These results demonstrate that the ECT is effective for optimizing the microstructure and properties of HEA microfibers, and can also provide both a theoretical foundation and technical guidance for fabricating high-performance metallic microfibers.
DATA PAPERS
First-principles study of electronic structure, elastic properties and hardness of Cr-doped CuZr2
WANG Kun, XU Heyan, ZHENG Xiong, ZHANG Haifeng
2025, 74 (13): 137101. doi: 10.7498/aps.74.20250264
Abstract +
In recent years, the design and development of new high-performance alloys based on first principles have received extensive attention. However, there are few reports on the structural design and thermodynamic properties of Cu-Zr alloys at nanoscale. In this work, based on the crystal structure characteristics of CuZr2, 12 kinds of Cr-doped CuZr2 structures are designed and optimized by the method of Cr atom doping through the first-principle calculation based on the density functional theory, and 6 kinds of mechanically and dynamically stable doped structure models are found. By calculating the electronic structure, elastic properties and hardness of the CuZr2 and its dynamically stable Cr-doped structure, it is found that the studied objects have all energy bands that pass through the Fermi energy level and are metallic. The main contributors to the metallic properties of the CuZr2 are the p and d orbital electrons of Zr, while the main contributors to the metallic properties of the 6 dynamically stable Cr-doped CuZr2 structures are the p and d orbital electrons of Cr and Zr. Meanwhile, CuZr2 has symmetrically distributed spin electrons, which do not show magnetism externally. However, the doping of Cr atoms increases the elemental species of the matrix. In addition to the difference of spin electrons brought by the d-orbital electrons of Cr atoms, the doped Cr atoms destroy the symmetrical distribution of electrons with different spin directions in the p- and d-orbitals of Zr atoms in the matrix, so that the designed 6 dynamically stable Cr-doped CuZr2 structures exhibit ferromagnetic properties with magnetic moments ranging from 0.303μB to 5.243μB. In addition, it is found that Cr atoms can improve the mechanical properties of CuZr2. When the Cr atom is used to replace the Zr atom in the matrix, the elastic modulus and hardness of the material can be improved, and when the Cr atom is used to replace the Cu atom in the matrix, the machining properties of the material can be improved due to the reduction of hardness. The datasets presented in this work, including the band structure, density of states, and phonon dispersion frequency, are available from https://www.doi.org/10.57760/sciencedb.j00213.00122.
RAPID COMMUNICATION
Fano resonance in high-order harmonics of hydrogen atoms driven by intense laser pulse
CHEN Suqi, HE Feng
2025, 74 (13): 133202. doi: 10.7498/aps.74.20250617
Abstract +
We numerically solved the time-dependent Schrödinger equation (TDSE) for a hydrogen atom interacting with intense near-infrared laser fields to investigate the mechanism of below-threshold high-harmonic generation (HHG). The primary focus was on understanding the spectral features, particularly resonant structures, arising in the fifth harmonic region under specific driving conditions. Our simulations utilized a laser wavelength of 608 nm. At this wavelength, hydrogen atoms can resonantly absorb five photons, promoting electrons from the ground stateWe numerically solved the time-dependent Schrödinger equation (TDSE) for a hydrogen atom interacting with intense near-infrared laser fields to investigate the mechanism of below-threshold high-harmonic generation (HHG). The primary focus was on understanding the spectral features, particularly resonant structures, arising in the fifth harmonic region under specific driving conditions. Our simulations utilized a laser wavelength of 608 nm. At this wavelength, hydrogen atoms can resonantly absorb five photons, promoting electrons from the ground state $\left|1{\mathrm{s}}\right\rangle$ to the excited state $\left|2{\mathrm{p}}\right\rangle$. Concurrently, the atom can absorb additional photons leading to ionization. Crucially, due to the AC Stark shift induced by the intense laser field (laser dressing), some laser-dressed continuum states $\left|{\mathrm{c}}\right\rangle$ become energetically degenerate with the laser-dressed $\left|2{\mathrm{p}}\right\rangle$ state. High-harmonic radiation at the fifth harmonic frequency can then be emitted via two distinct quantum paths: 1) Bound-bound recombination: Direct recombination from the laser-dressed $\left|2{\mathrm{p}}\right\rangle$ state back to the ground state $\left|1{\mathrm{s}}\right\rangle$. 2) Continuum-bound recombination: recombination from the laser-dressed continuum states $\left|{\mathrm{c}}\right\rangle$ (reached via ionization) back to $\left|1{\mathrm{s}}\right\rangle$. Both pathways emit photons of identical energy corresponding to the fifth harmonic. Our important finding is that there is significant quantum interference between these two recombination channels. This interference is manifested in the spectrum as an asymmetric Fano lineshape of the fifth harmonic intensity profile. Furthermore, we demonstrate that the shape of this Fano resonance exhibits strong and controllable dependence on the intensity of the driving laser field. This study provides clear evidence that Fano quantum interference, typically associated with multi-electron correlations or autoionizing states in complex systems, can emerge in the fundamental single-electron hydrogen atom system under the condition of intense laser field. The interference arises directly from the coherent superposition of the bound-bound and continuum-bound recombination pathways caused by laser-induced degeneracy. Importantly, by adjusting the laser intensity the spectral profile of the Fano resonance can be actively manipulated, providing a novel method for coherently controlling the harmonic emission in simple atomic systems.
REVIEW
Research progress of single event effect and reinforcement technology of SiC power metal-oxide-semiconductor field-effect transistors
QIU Yiwu, LIANG Di, YIN Yanan, DONG Lei, WANG Tao, ZHOU Xinjie
2025, 74 (13): 136103. doi: 10.7498/aps.74.20250273
Abstract +
In extreme radiation environments, such as space nuclear reactor systems, deep-space probe power modules, and launch vehicle propulsion systems, high-voltage and high-power devices demonstrate significant practical value. Silicon carbide (SiC) metal-oxide-semiconductor field-effect transistors (MOSFETs) possess advantages including high breakdown voltage, thermal stability, and low on-state resistance, enabling further improvements in aerospace power supply efficiency. Therefore, research on radiation effects and radiation-hardening techniques for SiC power MOSFETs has rapidly emerged as a critical focus in the industry. Firstly, this paper reviews the developmental evolution of SiC power MOSFETs, analyzes the necessity of transitioning from planar gate to trench-gate architectures, and provides future prospects for advanced SiC power MOSFET technologies. Secondly, it systematically compiles current research achievements in single event burnout (SEB) and single event gate rupture (SEGR) caused by heavy ion irradiation in SiC power MOSFETs. Finally, based on a mechanistic analysis of radiation-induced single event damage in SiC power MOSFETs, this study summarizes recent progress of radiation-hardening technologies, aiming to provide valuable ideas for understanding radiation induced failure mechanisms and enhancing the radiation tolerance of SiC power MOSFETs.
Research progress of rare earth magnetic materials based on machine learning
LIU Dan, LI Yuan, SUN Ruoxuan, QI Xing, SHEN Baogen
2025, 74 (13): 137102. doi: 10.7498/aps.74.20250431
Abstract +
Rare-earth elements possess unique atomic structures characterized by multiple unpaired 4f orbital electrons in inner shells, high atomic magnetic moments, and strong spin-orbit coupling. These attributes endow them with rich electronic energy levels, enabling them to form compounds with different valence states and coordination environments. Consequently, rare-earth materials typically exhibit excellent magnetic properties and complex magnetic domain structures, making them critical for the development of high-tech industries. The intricate magnetic configurations, different types of magnetic coupling, and direct/indirect magnetic exchange interactions in these materials not only facilitate the development of novel functional devices but also pose significant challenges to fundamental research. With the rapid advancement of data mining techniques, the emergence of big data and artificial intelligence provides researchers with a new method to efficiently analyze vast experimental and computational datasets, thereby accelerating the exploration and development of rare-earth magnetic materials. This work focuses on rare-earth permanent magnetic materials, rare-earth magnetocaloric materials, and rare-earth magnetostrictive materials, detailing the application progress of data mining techniques in property prediction, composition and process optimization, and microstructural analysis. This work also delves into the current challenges and future trends, aiming to provide a theoretical foundation for deepening the integration of data mining technologies with rare-earth magnetic material research.
GENERAL
Machine learning model predicted thermodynamic stability of rare earth compounds
QIN Chenglong, ZHAO Liang, JIANG Gang
2025, 74 (13): 130201. doi: 10.7498/aps.74.20250362
Abstract +
This study aims to predict the thermodynamic stability of rare-earth compounds by using machine learning (ML) models, providing crucial data support for designing advanced materials and facilitating the discovery of new rare-earth compounds.In terms of methods, this study is based on a dataset consisting of 280569 compounds. The formation energies of these compounds are calculated by density functional theory (DFT). A system consisting of 145 feature descriptors is constructed, covering stoichiometric properties, statistical properties of elements, electronic structure properties, and properties of ionic compounds, comprehensively describing the characteristics of rare-earth compounds. Two ML models, i.e. random forest (RF) and neural network (NN), are selected to perform classification and regression tasks respectively. The 5-fold cross-validation is used to improve the reliability of the models. The min-max scaling technique is used for preprocessing data, and an ensemble learning architecture is constructed to address the limitations of single model.In the classification task, the RF and NN algorithms perform remarkably well. With 5-fold cross-validation, the accuracy reaches approximately 0.97, and the F1 score is around 0.98, enabling the precise classification of compounds into stable or unstable categories. In the regression task, the mean absolute errors (MAEs) of the formation energy predictions by the RF and NN models are 0.055 eV/atom and 0.071 eV/atom, respectively. This indicates that the model predictions are highly accurate and can replace complete DFT calculations to a certain extent. In the predictive analysis of system outside the test set, six representative components are selected from the material project database, covering binary, ternary, and quaternary systems. The prediction errors of all compositions are controlled within 0.5 eV/atom, with an error percentage of lower than 25%, indicating that the model has strong ability of extrapolation and prediction. When predicting the binary phase diagrams of rare-earth compounds La-Al and Ce-H by using the trained models, the convex hull phase diagrams constructed through the ensemble learning architecture, which combines the prediction results of the RF and NN models, are highly consistent with those constructed from the open quantum materials database. The models successfully capture several metastable phases that are not present in multiple databases. Moreover, the convex hull distances of the predicted phases are mostly less than 0.1 eV/atom, with the maximum not exceeding 0.2 eV/atom.In conclusion, this study successfully uses ML models to predict the thermodynamic stability of rare-earth compounds. The constructed models demonstrate strong capabilities in classification and regression tasks. The ensemble learning architecture effectively improves the model performance, providing a promising tool for discovering materials in the field of rare-earth science, contributing to the research and development of new rare-earth compounds, and designing advanced materials.
Quantum phase transition and exceptional points in dissipative optical cavity with nonlinear atom-photon interaction
LIU Ni, LUO Yunqing, LIANG Jiuqing
2025, 74 (13): 130301. doi: 10.7498/aps.74.20250075
Abstract +
Dicke model, as an important many-body model in quantum optics, describes the interaction between multiple identical two-level atoms and a quantized electromagnetic field. This spin-boson model shows collective phenomena in light-matter interaction systems and can undergo a second-order quantum phase transition from a normal phase to a superradiant phase when the coupling strength between the two-level atoms and the optical field exceeds a critical value. Dicke model embodies unique many-body quantum theories. And it has been widely studied and obtained many significant research results in quantum information, quantum process, and other quantum systems. Meanwhile, Dicke model also has wide applications in quantum optics and condensed matter physics.The extended Dicke model, describing the interaction of a Bose-Einstein condensate in an optical cavity, provides a remarkable platform for studying extraordinary quantum phase transitions in theory and experiment. Based on the recent experiment on non-Hermitian coupling between two long-lived atomic spin waves in an optical cavity, in this work we use spin-coherent-state variational method and present the macroscopic quantum-state energy of the non-Hermitian Dicke model.The spin coherent state variational method has an advantage in the theoretical research of macroscopic quantum states, especially in the normal and the inverted pseudospin states. In the variational method, optical coherent states and atomic extremum spin coherent states are used as the trial wave functions. A Hermitian transformation operator is proposed to diagonalize the non-Hermitian Hamiltonian, which is different from the ordinary quantum mechanics where the transformation operator must be unitary. Herein, the energy function is not necessarily real in the entire coupling region. Beyond an exceptional point, the spectrum becomes complex and introducing biorthogonal sets of atomic extremum states is necessary to evaluate the average quantities.The normal phase (for the zero average photon number) possesses real energy and atomic population. The non-Hermitian interaction destroys the superradiant phase (for the stable nonzero average photon number) and leads to the absence of quantum phase transition. However, the introduced atom-photon interaction, which is induced by the pump field experimentally, can change the situation, dramatically. The pump field can balance the loss by the non-Hermitian atom-photon interaction to achieve the superradiant phase.An interesting double exceptional point are observed in the energy functional. There is the real spectrum below the first exceptional point and beyond the second exceptional point, while there is a complex spectrum between these two exceptional points. The superradiant phase appears only beyond a critical value, which is related to the nonlinear interaction and the pump laser. A new and inverted quantum phase transition from the superradiant phase to the normal phase, is observed by modulating the atom-field coupling strength. The superradiant phase of the population inversion state appears for a negative effective frequency and a large atom-photon interaction. The influence of the dissipative coupling may be observed in cold atom experiment in an optical cavity. All the parameters adopted in this work are the actual experimental parameters.
Hyperentanglement W state concentration for polarization-time-bin photon systems with linear optics
GUO Pengliang, XI Shun, GAO Chengyan
2025, 74 (13): 130302. doi: 10.7498/aps.74.20241642
Abstract +
In recent years, quantum communication technology has developed rapidly, and quantum communication schemes based on hyperentangled states have attracted widespread attention due to their efficiency and security. However, in practical communication, maximally hyperentangled states are highly susceptible to environmental noise, which causes them to degrade into non-maximally hyperentangled states. This degradation significantly reduces the fidelity of the quantum information and communication efficiency. In this article, we propose an efficient entanglement concentration scheme to restore degraded polarization-time hyperentangled W states, thereby enhancing the reliability and transmission distance of multiparty quantum communication. The protocol employs the parameter-splitting approach, where the receiver performs local operations on received non-maximally hyperentangled photons by using linear optical elements, achieving hyperentanglement concentration through detector responses and post-selection. This method eliminates the need for auxiliary photons, thereby reducing the use of quantum resources and maintaining operational simplicity. Moreover, the scheme can be extended to N-photon hyperentangled W states. The theoretical calculations demonstrate that the success probability of the protocol is determined by the minimal parameter of the hyperentangled state, exhibiting a monotonic increase as this parameter grows. Under ideal conditions, the maximum success probability approaches unity and the success probability improves with the number of entangled photons increasing. When considering the efficiency of practical optical components, the maximal success probabilities for hyperentangled W states with N = 3, 4, and 5 are found to be 0.856, 0.791, and 0.732, respectively. Consequently, the proposed scheme efficiently concentrates the degraded polarization-time hyperentangled W state into the maximally hyperentangled state. This work is of significant importance for long-distance information transmission and provides theoretical references for implementing long-distance multi-party quantum communication.
Design and dynamic mechanism of M-shaped low-frequency isolation structure imitating limb configuration of arthropods
LAN Chunbo, JIA Jie, WANG Yang, WANG Shuo, ZHANG Lu
2025, 74 (13): 130501. doi: 10.7498/aps.74.20241203
Abstract +
Arthropods, including spiders and mantises, can maintain their body stability on shaking surfaces, such as spiderwebs or leaves. This impressive stability can be attributed to the specific geometric shape of their limbs, which exhibit an M-shaped structure. Inspired by this geometry, this work proposes an arthropod-limb-inspired M-shaped structure for low-frequency vibration isolation. First, the design method of the M-shaped quasi-zero-stiffness (QZS) structure is presented. A static analysis of potential energy, restoring force, and equivalent stiffness is conducted, showing that the M-shaped structure enables a horizontal linear spring to generate nonlinear stiffness in the vertical direction. More importantly, this nonlinear stiffness effectively compensates for the negative stiffness in large-displacement responses, thereby achieving a wider quasi-zero-stiffness region than the conventional three-spring-based QZS structure. Subsequently, the harmonic balance method is employed to derive approximate analytical solutions for the M-shaped QZS structure, which are well validated through numerical simulation. A comparison between the proposed M-shaped QZS structure and the conventional three-spring-based QZS structure is performed. Results show that the M-shaped QZS structure is advantageous for reducing both the cut-in isolation frequency and the resonance frequency. In particular, under large excitation or small damping conditions, the performance improvement of the M-shaped QZS structure in terms of reducing the resonance frequency and maximum response becomes more pronounced. The underlying mechanism behind this feature is primarily attributed to the expanded QZS region induced by the M-shaped structure. Finally, since the M-shaped structures vary among different arthropods, the effect of the geometry of M-shaped structures on low-frequency vibration performance is investigated. Interestingly, a trade-off between vibration isolation performance and loading mass is observed. As the M-shaped structure becomes flatter and the QZS region expands, the cut-in isolation frequency, resonance frequency/peak, and loading mass all decrease. This occurs because a flatter M-shaped structure leads to a reduction in the equivalent stiffness generated by the horizontal stiffness. Therefore, as the loading mass capacity decreases, the low-frequency vibration isolation performance is enhanced. This novel finding provides a reasonable explanation for why most arthropods possess many pairs of limbs, allowing the loading mass to be distributed while achieving excellent low-frequency vibration isolation.
NUCLEAR PHYSICS
Measurement of 85Sr(n, γ) cross sections of unstable nuclei in HI-13 tandem accelerators by surrogate-reaction method
WANG Hanyu, QIU Yijia, LIN Chengjian, WU Xiaoguang, HAN Yinlu, WU Hongyi, FENG Jing, ZHENG Yun, YANG Lei, LI Congbo, LUO Tianpeng, CHANG Chang, SUN Qi, ZHU Deyu, ZHAO Yixuan, HUANG Dahu, LI Tianxiao, ZHENG Min, ZHAO Zihao, ZHU Yiwei, ZHAO Kunling, SUN Pengfei, SONG Jinxing, GUO Mingwei, REN Sixi, ZHENG Xiaohai
2025, 74 (13): 132501. doi: 10.7498/aps.74.20250214
Abstract +
Neutron capture cross sections, as important parameters for describing the probability of neutron-nucleus reactions, play a key role in multiple scientific fields. In astrophysics, neutron capture cross section data are essential elements for understanding stellar nucleosynthesis processes. In particular, in extreme environments such as supernova explosions and neutron star mergers, accurate neutron capture cross sections can reveal the secrets of heavy element formation. In the field of national security, neutron capture cross sections are crucial for the design of nuclear weapons and the security of nuclear materials. By accurately grasping the neutron capture characteristics of different nuclides, the nuclear reaction process can be optimized to ensure strategic security. In addition, in the simulation of nuclear power generation, neutron capture cross section data are the basis of reactor design and operational analysis. Through in-depth research on and precise measurements of neutron capture cross sections, the safety and efficiency of nuclear reactors can be improved, thus promoting the sustainable development of nuclear energy. At present, there is little research on the neutron capture cross sections of nuclides with half-lives of only a few years or even shorter, mainly due to the complexity of measurement techniques and the instability of the nuclides themselves. The neutron capture cross section data of these nuclides are crucial for reactor design, nuclear medicine applications, and nuclear waste treatment. Further research requires the development of more advanced detection techniques and theoretical models to accurately measure and predict their neutron capture behavior.The surrogate-reaction method, as an effective measurement means, plays an important role in studying reaction cross sections of short-lived nuclides. Its basic idea is to indirectly obtain the reaction cross section information of short-lived nuclides by measuring the specific particles emitted by stable nuclides. Specifically, when stable nuclides are bombarded by high-energy particles, nuclear reactions will occur and specific particles will be released. By accurately measuring the energies, angles, and numbers of these particles, the cross sections of short-lived nuclides in the corresponding reaction can be inferred. This method can not only overcome the technical difficulties in directly measuring short-lived nuclides, but also improve the accuracy and reliability of the measurement results, which provides important support for nuclear physics research. In addition, the surrogate-reaction method also shows broad application prospects in the fields of nuclear technology application and nuclear data assessment.The experiment is carried out on the Beijing HI-13 tandem accelerator at the China Institute of Atomic Energy. 89Y is bombarded with 22 MeV protons, and the 85Sr(n, γ) cross section is measured through the (p, αγ) reaction. The telescope array composed of silicon strip detectors can effectively identify the reaction products. By precisely measuring parameters such as the energies and angles of particles, the array can distinguish different nuclides, thus determining the outgoing particles. Combined with the γ-ray energy spectrum analysis of the HPGe detector, the (n, γ) reaction cross section data of 85Sr under the Weisskopf-Ewing (W-E) approximation are extracted. Due to the mismatch of the Jπ population between the existing alternative reactions and direct reactions, it is necessary to compensate for this mismatch and then correct the results. In order to obtain relatively reliable results, the Jπ population calculated by TALYS is used to revise the experimental data of the (n, γ) cross section.These results indicate that the cross section of 85Sr(n, γ) varies with neutron energy in a specific energy range, which is consistent with the trend of the existing international evaluation library data. This validates the effectiveness of cross section measurement as an alternative reaction method, thereby providing an important experimental basis for further exploring the nuclear reaction mechanism and nuclear data application. This method has reference significance for the cross section measurement of other nuclides.
Influence of high-momentum distribution within nucleus on production of bremsstrahlung photons
LIU Xun, FAN Yan, GUO Wenjun
2025, 74 (13): 132502. doi: 10.7498/aps.74.20250239
Abstract +
In this work, semi-classical quantum molecular dynamics is used to investigate the influence of high momentum distribution on nuclear reaction systems by using photons produced by nucleon bremsstrahlung as indicators. The research examines the relationship between this influence and the incident energy, collision parameters, and the differences in isotopic spin cross-sections. Under the condition of a 20% high-momentum distribution and ensuring the conservation of nuclear energy, a system different from traditional configurations is constructed by sampling neutrons and protons using the Monte Carlo method, with the selected nucleons exhibiting characteristics of high-momentum nucleons. The influence of high-momentum distribution within the nucleus on bremsstrahlung photons is analyzed through the collision results of heavy ions in nuclear systems spanning from light to heavy species. The results indicate that, at an incident energy value of 50 MeV/u, the collision system studied in this work exhibits higher photon density in the high-energy region of high momentum distribution system than traditional system for nuclear systems ranging from light-mass system (18O+18O), lower-mass system (46Ca+46Ca), and medium-mass system (86Kr+86Sr) to heavy-mass system (124Sn+124Sn), than those, while there is no significant difference in photon density in the low-energy region. At a collision parameter b = 0 fm, the energy shift phenomenon of photons produced by collisions becomes more pronounced with the increase of incident energy, peaking at E = 150 MeV/u. This energy transfer phenomenon induced by high momentum distribution typically exists within the collision parameter range from b = 0 fm to b = 6 fm. When considering isotopic spin cross-sections, high momentum distribution can affect the collision probability of the system. Therefore, high momentum distribution has a significant influence on nuclear reaction system, closely related to incident energy and isotopic spin cross-section.
ATOMIC AND MOLECULAR PHYSICS
Non-adiabtic population transfer in SU(3) Lie algebra based three-level system without the rotating wave approximation
DOU Yuxin, LI Wei
2025, 74 (13): 133201. doi: 10.7498/aps.74.20241504
Abstract +
Coherent population transfer in quantum systems is of fundamental importance in many fields such as atomic and molecular collision dynamics, and information processing for qubit systems. Stimulated Raman nonadiabatic passage technique, when implemented in an externally driven three-level system, provides an efficient approach for realizing accelerated population transfer while maintaining robust quantum coherence, with no need of rotating wave approximation. However, previous protocols employ multiple pulses and imply that Rabi frequencies have a few oscillations during dynamical evolution. In this work, under the condition of two-photon resonance, we use a gauge transformation method to inversely design a Λ-configuration three-level system that can be solved exactly. By using SU(3) transformation, we establish the relationship between Schrödinger representation and gauge representation, where the effective Hamiltonian is an Abelian operator. Subsequently, we construct the desired Hamiltonian and further investigate its dynamic behavior. The result shows that by imposing appropriate boundary conditions on the control parameters, high-fidelity population transfer can be achieved in ideal evolution. In addition, for the practical case with pulse truncation and intermediate state decay, the fidelities of specific models can reach about 99.996% and 99.983%, respectively. Compared with other existing nonadiabatic quantum control schemes, the present scheme has the distinctive advantages. We achieve the required quantum control by applying only a few sets of Stokes and pump pulses without introducing any additional microwave field. This method does not exhibit Rabi oscillations in the dynamic process, nor does it produce singularities in the pulse itself.
ELECTROMAGNETISM, OPTICS, ACOUSTICS, HEAT TRANSFER, CLASSICAL MECHANICS, AND FLUID DYNAMICS
Ultrahigh-resolution mass sensing based on bound states in continuum domain in double-cavity optomechanical system
DONG Yaoyong, WU Yi, ZHENG Xuejun, WANG Denglong, ZHAO Peng
2025, 74 (13): 134201. doi: 10.7498/aps.74.20250063
Abstract +
Resolution is one of the key indicators in the cavity optomechanical mass sensing. The bound states in the continuum (BIC) enable extremely narrow linewidths, which have great potential for enhancing the resolution of cavity optomechanical mass sensors. In order to enhance the resolution of cavity optomechanical mass sensing, we propose a simple double-cavity optomechanical system under the blue-detuning condition to realize the BIC singularity, and present an ultrahigh-resolution mass sensing scheme based on BIC in this paper. By solving the linearized Heisenberg-Langevin equations, the expressions for the susceptibility and transmission rate of the system are derived. Based on the system’s susceptibility, we study the absorption characteristics of the probe field under the blue-detuning condition. The absorption spectrum of the system exhibits three peaks, among which the central narrow peak exhibits optical gain characteristics, collectively forming a phenomenon analogous to double optomechanically induced transparency. Then, analysis of the dressed-state energy-level structure reveals that the formation of the central narrow peak stems from quantum interference effects in a double-Λ-type dark-state resonance. The linewidth evolution of the quasi-BIC central narrow peak is investigated by analyzing the dependence of the real part and imaginary part of the corresponding eigenvalue on the optomechanical coupling strength. It can be found that the imaginary part of the eigenvalue for the central narrow peak becomes zero when the optomechanical cooperativity coefficient equals the double-cavity cooperativity coefficient plus one, enabling the realization of BIC. The linewidth of the central peak is ultrasmall under this BIC condition, and the shift of the transmission peak in the transmission spectrum is linearly related to the adsorbed mass. Based on these characteristics, the system under the BIC condition can achieve mass sensing with an ultrahigh resolution, with a resolution of approximately 1 ag. Meanwhile, the linewidth of the transmission peak can be suppressed below 1 Hz, which is superior to the traditional optomechanical mass sensing schemes based on four-wave mixing, photonic molecules, and plasmon polaritons. Systematic investigation of eigenvalue variations and the corresponding sensitivity enhancement factors under mechanical resonator frequency shift reveals that the real part and the imaginary part of the eigenvalue associated with the central peak exhibit negligible variations under such perturbations. This indicates that the mass sensing scheme based on BIC in the double-cavity optomechanical system can maintain ultrahigh resolution and precise mass measurement under mechanical resonator frequency shift. Our scheme provides an approach for realizing the BIC singularity in optomechanical systems, and presents a new route to improving the resolution of mass sensors based on cavity optomechanical systems.
Numerical inversion of PT-symmetric potential functions for (1+1)-dimensional nonlinear Schrödinger equations
WANG Yang, XU Yinghong, ZHAO Yedan, ZHANG Lipu
2025, 74 (13): 134203. doi: 10.7498/aps.74.20250129
Abstract +
The inverse problem of reconstructing the PT-symmetric potential in a class of (1+1)-dimensional nonlinear Schrödinger equation is investigated in this work. The governing equation is given as follows:      $ \text{i}u_t(x,t) + u_{xx}(x,t) + \alpha\left| u(x,t) \right|^2 u(x,t) + \beta\left| u(x,t) \right|^4 u(x,t) + V_{\rm PT}(x) u(x,t) = 0, $where u(x, t) denotes the wave function in dimensionless coordinates, and the PT-symmetric potential VPT(x) = V(x)+iW(x) consists of a real part V(x) and an imaginary part iW(x), satisfying the symmetry conditions V(x) = V(–x) and W(x) = –W(x).In this inverse problem, partial boundary values of the wave function are known, while the potential $ V_{\rm PT}(x) $ is the unknown to be reconstructed. To address this challenge, we construct a three-level finite difference scheme for the corresponding forward problem, discretizing both the wave function and the potential. This approach leads to a nonlinear system of equations that links the known wave data to the unknown potential values. To simplify the computation, we separate the real and imaginary parts of this system and reformulate it as a real-valued nonlinear system of equations.For the numerical solution, we employ an inexact Newton method to iteratively solve the resulting nonlinear system. In each iteration, the Jacobian matrix is approximated numerically. To ensure that the reconstructed potential strictly satisfies the PT-symmetry, a parity correction mechanism is introduced at the end of the iteration process.We conduct numerical experiments under both noise-free (exact data) and noisy (inexact data) conditions. The results indicate that in both cases, the proposed method converges within a limited number of iterations and keeps the reconstruction error within the order of 10–3. These findings verify the effectiveness and robustness of this method in solving inverse problems of PT-symmetric potentials, and provide a new idea and powerful method for related numerical applications.
Robustness anylysis of periodic oscillations in continuous-time crystals
LI Hongxia, LI Tingmei, WANG Zengpu, CHEN Yuhui, ZHANG Xiangdong
2025, 74 (13): 134204. doi: 10.7498/aps.74.20250036
Abstract +
Continuous time crystals represent a novel state in many-body systems that can self-organize into timeperiodic oscillations without external periodic driving. Recent experiments have achieved such systems in dissipative solid-state materials, where persistent temporal order is autonomously sustained. A decisive characteristic of time crystals is their robustness, meaning that despite various disturbances, including fluctuations in internal parameters and external noise, they can still maintain rhythmic behavior, which has scientific value and echnological application potential. Although previous studies have shown that specific experimental parameters have robustness, thare is a lack of a systematic framework for quantifying and predicting their resilience to disturbances, and the underlying physics of this robustnessis still not fully understood. The key unresolved problems include how nonlinear interactions and feedback mechanisms contribute to stability, and what the critical thresholds are for parameter variations beyond which temporal order collapses.This work addresses these gaps by systematically analyzing how internal parameters and external influences affect the oscillation period and overall stability. Internally, the dynamics are determined by dipole-dipole interactions and atomic transition strengths, which define the temporal symmetry breaking that occurs in the system. Externally, the response of the system is controlled by the strength of the optical driving field and the energy dissipation rate. A key finding is the determination of an intrinsic feedback mechanism for a dynamically stabile time crystal. This mechanism plays a role in restoring force, correcting deviations caused by minor disturbances, and maintaining the coherence of oscillatory phase.Moreover, the system displays nonlinear dynamical behavior, characterized by two different states: one is stable oscillation continuing under moderate disturbance, and the other is stronger disturbance causing dynamical phase transition, resulting in switching between disordered or dynamically unstable state and stable state. These results provide a comprehensive understanding of the various behaviors observed in continuous time crystals and lay an important theoretical foundation for utilizing their unique properties in advanced applications such as quantum information processing and precision metrology.
X-ray tomography based experimental study on packing structure of “hexapod” concave particles under external tapping
LUO Rudan, ZENG Zhikun, GE Zhuan, JIANG Yonglun, WANG Yujie
2025, 74 (13): 134501. doi: 10.7498/aps.74.20250231
Abstract +
Granular materials are ubiquitous in nature and industrial production. Investigating the structure of packing is crucial for understanding the physical properties of granular materials. Owing to their symmetry and simple geometry, spherical particles have long served as an ideal model for studying granular packing, yielding numerous research outcomes.In recent years, the influence of particle shape on packing structures has drawn considerable attention. Non-spherical particles, characterized by complex shapes, tend to interlock and form stable structures. Their significant geometric cohesion notably affects the stability and porosity of granular packing.To investigate the structural evolution and compaction mechanisms of three-dimensional concave particles (hexapod-shaped) under external tapping, focusing on the role of geometric cohesion in enhancing mechanical stability, we employ hexapod-shaped particles that are composed of three mutually orthogonal spherocylinders in this study. The granular system subjected to consecutive tapping can reach a stationary state. In the densifying process of the system, packing structures with different volume fractions will be formed. Meanwhile, by combining with X-ray tomography, we can obtain the microstructure.The findings reveal that the volume fraction of “hexapod” particle packing is significantly lower than that of hard-sphere systems. The compaction curves of “hexapod” particles across varying tapping intensities are accurately described by the Kohlrausch-Williams-Watt (KWW) law, which is consistent with hard-sphere system, suggesting a relaxation process governed by heterogeneous modes. Furthermore, both the volume fraction of the steady-state granular packing and the average contact number exhibit an inverse relationship with tapping intensity, increasing as the intensity decreases. A detailed statistical analysis of contact points indicates that the compaction process of “hexapod” particles is predominantly influenced by two factors: the augmentation in the number of neighboring contacting particles and the modification of contact forms. These factors collectively enhance the degree of interlocking among hexapods within the system. Specifically, the compaction process is primarily propelled by the escalation in neighboring contacts and the refinement of contact types, particularly the increase in cylinder-cylinder (cc) contact. This rise in cc contact significantly enhances mechanical stability through strengthening geometric interlocking.This study reveals the structural evolution characteristics of non-spherical particles in the compaction process and provide important experimental support for understanding the unique mechanical and dynamic properties of concave particle packing. This research not only enriches the experimental data of granular packing structures but also offers a new perspective for exploring the universal laws of packing for particles of different shapes. This study is to lay a more solid foundation for the theoretical research and industrial applications of granular materials, thereby promoting technological progress and innovation in related fields.
Structures of binary hard-sphere packing under tapping
LIU Zelin, YUAN Houfei, ZENG Zhikun, JIANG Yonglun, GE Zhuan, WANG Yujie
2025, 74 (13): 134502. doi: 10.7498/aps.74.20250232
Abstract +
The packing behavior and mechanical properties of granular materials play a critical role in various engineering applications, including materials handling, construction, and energy storage. Although significant progress has been made in understanding the packing of monodisperse spheres, real-world granular systems often exhibit polydispersity, where particles of different sizes coexist. Binary systems, where the particle size ratio is adjustable, serve as a simplified model to study the structural and dynamical properties of granular materials. However, most theoretical studies on binary systems focus on idealized frictionless models, neglecting the coupled effects of friction and preparation history, and experimental data for three-dimensional systems remain limited. This study seeks to address these gaps by investigating the packing behavior of binary hard spheres under tapping through using advanced experimental techniques such as X-ray computed tomography (CT) and tap-driven compaction. The effects of particle size ratio and tap intensity on the packing fraction and local structure of binary granular systems are investigated systematically. The experimental results show that the steady-state packing fraction decreases as tap intensity increases, exhibiting similar behavior at different composition ratios. Additionally, the compaction dynamics are quantified using the Kohlrausch-Williams-Watts (KWW) relaxation function, revealing that the relaxation time decays exponentially with tap intensity increasing, independent of the composition ratio. Voronoi cell analysis demonstrates that the local volume distribution of each component in a bidisperse system composed of big particles and small particles is highly similar to that in a monodisperse system. Notably, as tap intensity decreases, the system density increases, and volume fluctuation decreases, reflecting the trends observed in monodisperse packings. Furthermore, the study highlights the influence of friction on the packing structure. For binary systems, big particles, with rougher surfaces, pack more loosely than smaller particles, and the coordination number increases with the proportion of smaller particles increasing. This suggests that frictional interactions between particles play a significant role in determining the packing density and structural stability of granular materials. The average coordination number and the steady-state packing fraction are found to be weakly dependent on each other, with friction and tap intensity (or effective temperature) being the primary factors affecting the system’s structural characteristics. These findings provide a comprehensive experimental framework for understanding the packing behavior of binary granular systems, with important implications for material design in industrial applications. This study contributes to developing a more complete statistical mechanical theory for granular materials through combining both frictional effects and the influence of preparation history. Future research may extend these findings to more complex particle size distributions and explore the relationship between structural property and mechanical property.
Numerical investigation of droplet generation process in symmetric and asymmetric cross-junction microchannels
LI Xiang, LI Yan, LI Yingyan, DONG Zhiqiang, ZHUANG Xiaoru, ZHONG Zhigang, YU Peng
2025, 74 (13): 134701. doi: 10.7498/aps.74.20250474
Abstract +
Droplet microfluidics technology possesses significant potential applications in chemical analysis, biological detection, and material preparation. Passive droplet generation method can rapidly achieve droplet formation by using the geometric characteristics of microchannels and shear flow. As a typical structure, the influences of fluid parameters and symmetry differences in cross microchannels on the droplet generation process have not been fully studied. Therefore, this paper uses the lattice Boltzmann method to numerically simulate droplet generation in symmetric and asymmetric cross microchannels, thereby systematically analyzing the action mechanisms of capillary number, viscosity ratio, and microchannel symmetry. First, this study verifies the computational reliability of the numerical model through two classic cases, i.e. the droplet deformation under planar shear flow and stationary droplets on ideal solid surfaces. Then, this work focuses on studying the three flow stages in symmetric cross microchannels, i.e. interface immersion stage, shear-induced breakup stage, and the droplet migration and coalescence stage, and analyzes the synergistic mechanism of capillary number and viscosity ratio. In the symmetric cross microchannel structure, the capillary number is the main factor determining the droplet size in the cross microchannel. With the increase of the capillary number, the surface tension gradually weakens, causing the liquid bridge at the droplet neck to break more easily and generate droplets. In contrast, the effect of the viscosity ratio on the droplet size is relatively small, but it can suppress the generation of sub-droplets and improve the uniformity of droplets by adjusting the viscous resistance of the continuous phase. On this basis, this study further quantifies the influence of microchannel symmetry on the droplet generation process in cross microchannels. In the asymmetric cross microchannel structure, the microchannel deviation breaks the flow symmetry and weakens the cooperative shearing effect of the oil-phase fluid on the immersion structure of the water-phase fluid. When the microchannel deviates within the centerline range of the water-phase microchannel, the droplet size increases significantly with the increase of the microchannel deviation. This is mainly because the oil-phase fluid on the side far from the deviation first squeezes the immersion structure of the water-phase fluid, and then the oil-phase fluid near the deviation side exerts a secondary squeeze on the immersion structure, causing the neck liquid bridge of the immersion structure to continuously elongate and the shear position to shift along the microchannel deviation direction. When the microchannel deviation exceeds the centerline position of the water-phase microchannel, the interface fracture of the water-phase immersion structure mainly relies on the double squeezing effect of the oil-phase fluid and the surface tension of water-phase fluid, and the droplet size tends to be stable. The relevant research results lay a theoretical foundation for microchannel design and fluid parameter regulation in droplet microfluidics and thus further promote the application and development of droplet microfluidic technology.
PHYSICS OF GASES, PLASMAS, AND ELECTRIC DISCHARGES
Particle-in-cell/Monte Carlo collision simulations of dielectric barrier discharge packed with mixed dielectrics
ZHANG Lulu, LI Tianxiang, PANG Xuexia, GE Yuqi, LIU Xiaoqian, RAN Junxia, LI Qing, LI Xuechen
2025, 74 (13): 135201. doi: 10.7498/aps.74.20250208
Abstract +
Packed bed dielectric barrier discharge (PB-DBD) is extremely popular in plasma catalysis applications, which can significantly improve the selectivity and energy efficiency of the catalytic processes. In order to achieve some complex chemical reactions, it is necessary to mix different materials in practical applications. In this work, by using the two-dimensional particle-in-cell/Monte Carlo collision (PIC/MCC) method, the discharge evolution in PB-DBD packed with two mixed dielectrics is numerically simulated to reveal the discharge characteristics. Due to the polarization of dielectric columns, the enhancement of electric field induces streamers at the bottom of the dielectric columns with high electrical permittivity (εr). The streamers propagate downward in the voids between the dielectric columns with low εr, which finally converts into volume discharges. Then, a new streamer forms near the upper dielectric plate and propagates downward along the void of the dielectric columns with high εr. Moreover, electron density between the columns with high εr is lower than that between the dielectric columns with low εr. In addition, the numbers of e, $ {\text{N}}_{2}^{+} $, $ {\text{O}}_{2}^{+} $ and $ {\text{O}}_{2}^{-} $ present different profiles versus time. All of e, $ {\text{N}}_{2}^{+} $ and $ {\text{O}}_{2}^{+} $ increase in number before 0.8 ns. After 0.8 ns, the number of electrons decreases with time, while the numbers of $ {\text{N}}_{2}^{+} $ and $ {\text{O}}_{2}^{+} $ keep almost constant. In the whole process, the number of $ {\text{O}}_{2}^{-} $ keeps increasing with time increasing. The reason for the different temporal profiles can be analyzed as follows. The sum of electrons deposited on the dielectric and those lost in attachment reaction is greater than the number of electrons generated by ionization reaction, resulting in the declining trend of electrons. Comparatively, the deposition of $ {\text{N}}_{2}^{+} $ and $ {\text{O}}_{2}^{+} $ on the dielectric almost balances with their generation, leading to the constant numbers of $ {\text{N}}_{2}^{+} $ and $ {\text{O}}_{2}^{+} $. In addition, the variation of averaged electron density ($ {\bar{n}}_{{\mathrm{e}}} $) and averaged electron temperature ($ {\bar{T}}_{{\mathrm{e}}} $) in the voids between the dielectric columns are also analyzed under different experimental parameters. Simulation results indicate that both of them decrease with pressure increasing or voltage amplitude falling. Moreover, they increase with dielectric column radius enlarging. In addition, $ {\bar{n}}_{{\mathrm{e}}} $ increases and then decreases with the increase of N2 content in the working gas, while $ {\bar{T}}_{{\mathrm{e}}} $ monotonically increases. The variations of $ {\bar{n}}_{{\mathrm{e}}} $ and $ {\bar{T}}_{{\mathrm{e}}} $ in the voids can be explained as follows. With the increase of pressure, the increase of collision frequency and the decrease of average free path lead to less energy obtained per unit time by electrons from the electric field, resulting in the decreasing of $ {\bar{T}}_{{\mathrm{e}}} $. Moreover, the first Townsend ionization coefficient decreases with the reduction in $ {\bar{T}}_{{\mathrm{e}}} $, resulting in less electrons produced per unit time. Hence, both $ {\bar{n}}_{{\mathrm{e}}} $ and $ {\bar{T}}_{{\mathrm{e}}} $ decrease with pressure increasing. Additionally, $ {\bar{T}}_{{\mathrm{e}}} $ is mainly determined by electric field strength. Therefore, the rising voltage amplitude results in the increase of and $ {\bar{T}}_{{\mathrm{e}}} $. Based on the same reason for pressure, $ {\bar{n}}_{{\mathrm{e}}} $ also increases with the augment of voltage amplitude. Consequently, both $ {\bar{n}}_{{\mathrm{e}}} $ and $ {\bar{T}}_{{\mathrm{e}}} $ increase with voltage amplitude increasing. In addition, the surface area of dielectric columns increases with dielectric column radius enlarging. Therefore, more polarized charges are induced on the inner surface of the dielectric column, inducing a stronger electric field outside. Accordingly, the enlarging of dielectric column radius leads $ {\bar{n}}_{{\mathrm{e}}} $ and $ {\bar{T}}_{{\mathrm{e}}} $ to increase. Moreover, the variation of $ {\bar{n}}_{{\mathrm{e}}} $ with N2 content is analyzed from the ionization rate, and that of $ {\bar{T}}_{{\mathrm{e}}} $ is obtained by analyzing the ionization thresholds of N2 and O2.
Effect of magnetic field on contrast of Kα emission generated by hot electrons
XIAO Yangyang, WANG Xiaofang
2025, 74 (13): 135202. doi: 10.7498/aps.74.20250144
Abstract +
The interaction of a high-intensity laser with a solid target generates a large number of hot electrons. When these hot electrons are transported in the target material, X-rays, including Kα line and bremsstrahlung emissions are produced. The contrast of Kα line emission, i.e. the intensity of Kα line relative to the intensity of bremsstrahlung continua around the Kα line, depends on the anisotropy of the bremsstrahlung emission and is related to the energy and transportation of the hot electrons. In the past, some researchers used axial or annular magnetic fields to collimate hot electrons, but whether these magnetic fields can enhance the contrast of Kα emission has not been studied. In the present work, the effect of an axially uniform magnetic field or an annular magnetic field with a Gaussian distribution on the contrast of Cu Kα emission is investigated by Monte Carlo simulations. The simulation results and analysis show that the axially uniform magnetic field cannot strengthen the anisotropy of bremsstrahlung emission, so it cannot enhance the contrast of Kα emission efficiently. For the annular magnetic field with a Gaussian distribution, when an electron beam with a Boltzmann energy distribution is incident, due to the weak anisotropy of bremsstrahlung emission by low-energy electrons in the electron beam, the increase of Kα emission contrast is small. When an electron beam with a Boltzmann energy distribution, in which the low-energy part is cut off, or a mono-energetic electron beam is incident, the annular magnetic field with a Gaussian distribution significantly enhances the contrast of Kα emission in the back direction of the electron beam incidence. For an incident electron beam with an energy value in a range of 200–1000 keV, an annular magnetic field with a Gaussian distribution and a peak value of approximately 100 T is optimal for enhancing the contrast of Kα emission. Considering the existing experiments on generating annular magnetic fields and non-Boltzmann energy distribution hot electrons, it is possible to generate higher contrast Kα emissions with the enhancement of magnetic field in future experiments.
Numerical investigation of the secondary electron effect in capacitively coupled plasmas driven by ultra-low frequency/radio frequency sources
SHI Hanxu, LI Xinyang, ZHANG Yuru, WANG Younian
2025, 74 (13): 135203. doi: 10.7498/aps.74.20250341
Abstract +
In recent years, capacitively coupled plasmas driven by ultra-low frequency source have received increasing attention, because they are beneficial to generating ions with high energy and small scattering angle, which aligns well with the current trend in high aspect ratio etching. Since the sheath becomes thicker when an ultra-low frequency source is applied, the secondary electron emission becomes significant. Indeed, these energetic secondary electrons can enhance the ionization process and even affect the discharge mode. In this work, a two-dimensional fluid model is employed to study the influence of secondary electrons on the dual-frequency capacitively coupled plasma under different ultra-low frequency voltages, secondary electron emission coefficients and inter-electrode gaps. The high frequency is fixed at 13.6 MHz, and the ultra-low frequency is fixed at 400 kHz. First, by using the ion energy dependent secondary electron emission coefficient, it is shown that the electron density first decreases and then increases with ultra-low frequency voltage rising. This is because, on the one hand, the higher ultra-low frequency voltage leads to thicker sheath, and therefore, the effective discharge volume is compressed. On the other hand, secondary electrons emitted from electrodes can obtain more energy, thus enhancing the ionization process. By comparing with the results obtained with a fixed secondary electron emission coefficient, it is found that in the low voltage range, the evolution of the electron density is similar to that with a fixed coefficient of 0.1. While, in the high voltage range, the growth of the electron density is even more pronounced than that with a fixed coefficient of 0.2, indicating that the enhancement of the secondary electron effect by ultra-low frequency voltage is non-linear. Finally, the influence of discharge gap on the plasma properties is also discussed. It is shown that with the inter-electrode gap increasing from 2 to 4 cm, the maximum ionization rate becomes lower, but the electron density rises significantly, and the plasma radial uniformity is improved. When inter-electrode gap is large, secondary electrons can completely collide with neutral species, so their influence on the electron density at high ultra-low frequency voltage is more significant. The results obtained in this study contribute to understanding the influence of ultra-low frequency source on the secondary electron effect, and provide some guidance for optimizing plasma processing.
Study of circular cross-section plasmas in HL-2A tokamak: MHD equilibrium, stability and operational $ \boldsymbol{\beta } $ limit
SHEN Yong, DONG Jiaqi, SHI Zhongbing, HE Hongda, ZHAO Kaijun, PENG Xiaodong, QU Hongpeng, LI Jia, SUN Aiping
2025, 74 (13): 135204. doi: 10.7498/aps.74.20250450
Abstract +
Circular cross-section plasma is the most basic form of tokamak plasma and the fundamental configuration for magnetic confinement fusion experiments. Based on the HL-2A limiter discharge experiments, the magnetohydrodynamic (MHD) equilibrium and MHD instability of circular cross-section tokamak plasmas are investigated in this work. The results show that when $ {q}_{0}=0.95 $, the internal kink mode of $ m/n=1/1 $ is always unstable. The increase in plasma $ \beta $ (the ratio of thermal pressure to magnetic pressure) can lead to the appearance of external kink modes. The combination of axial safety factor $ {q}_{0} $ and edge safety factor $ {q}_{{\mathrm{a}}} $ determines the equilibrium configuration of the plasma and also affects the MHD stability of the equilibrium, but its growth rate is also related to the size of $ \beta $. Under the condition of $ {q}_{{\mathrm{a}}} > 2 $ and $ {q}_{0} $ slightly greater than 1, the internal kink mode and surface kink mode can be easily stabilized. However the plasma becomes unstable again and the instability intensity increases as $ {q}_{0} $ continues to increase when $ {q}_{0} $ exceeds $ 1 $. As the poloidal specific pressure ($ {\beta }_{{\mathrm{p}}} $) increases, the MHD instability develops, the equilibrium configuration of MHD elongates laterally, and the Shafranov displacement increases, which in turn has the effect on suppressing instability. Calculations have shown that the maximum $ \beta $ value imposed by the ideal MHD mode in a plasma with free boundary in tokamak experiments is proportional to the normalized current $ {I}_{{\mathrm{N}}} $ ($ {I}_{{\mathrm{N}}}={I}_{{\mathrm{p}}}\left({\mathrm{M}}{\mathrm{A}}\right)/a\left({\mathrm{m}}\right){B}_{0}\left({\mathrm{T}}\right) $), and the maximum specific pressure $ \beta \left({\mathrm{m}}{\mathrm{a}}{\mathrm{x}}\right) $ is calibrated to be $ ~2.01{I}_{{\mathrm{N}}},{\mathrm{ }}{\mathrm{i}}. {\mathrm{e}}. $ $ \beta \left({\mathrm{m}}{\mathrm{a}}{\mathrm{x}}\right)\sim 2.01{I}_{{\mathrm{N}}} $. The operational $ \beta $ limit of HL-2A circular cross-section plasma is approximately $ {\beta }_{{\mathrm{N}}}^{{\mathrm{c}}}\approx 2.0 $. Too high a value of $ {q}_{0} $ is not conducive to MHD stability and leads the $ \beta $ limit value to decrease. When $ {q}_{0}=1.3 $, we obtain a maximum value of $ {\beta }_{{\mathrm{N}}} $ of approximately $ 1.8 $. Finally, based on the existing circular cross-section plasma, some key factors affecting the operational $ \beta $ and the relationship between the achievable high $ \beta $ limit and the calculated ideal $ \beta $ limit are discussed.
Simulation of effect of neoclassical convection on tungsten impurity core accumulation in tokamak
LUO Yuchen, SANG Chaofeng, WANG Yilin, WU Yihan, ZHOU Qingrui, LI Jiaxian, XUE Miao, XUE Lei, ZHENG Guoyao, DU Hailong, WANG Dezhen
2025, 74 (13): 135206. doi: 10.7498/aps.74.20250384
Abstract +
Controlling of tungsten (W) impurity core accumulation is of great significance for the steady-state operation of tokamaks. This work mainly investigates the effect of neoclassical transport on the core accumulation of W impurities by using STRAHL code. The study focuses on the HL-3 device, which will use tungsten divertor and conduct research under argon gas injection discharge conditions. In the simulation, the edge and core background plasma parameters are obtained by SOLPS-ITER and OMFIT simulations, respectively. The distribution of tungsten impurities in the boundary region is simulated using the IMPEDGE code. The edge anomalous transport coefficient in STRAHL is adjusted accordingly, and the simulation results are compared with those from the IMPEDGE to ensure consistency in impurity distribution between the core and edge. In the core region, a numerical scan is performed to adjust the simulation results so that the energy radiation matches the setting values, thereby determining the specific turbulence convection velocity. By setting the coefficients for both the core region and the boundary region, a complete distribution of W impurities from boundary to the core is obtained. To account for the neoclassical transport effects, the neoclassical transport coefficients are calculated using the subroutine NEOART and applied to the impurity transport simulation, and the simulation region is set from $ \rho =0 $ to 0.9. On this basis, the transport of W impurities with and without neoclassical convection is simulated. The simulation results show that without neoclassical convection, anomalous transport dominates the impurity transport, which is inward and enhances impurity accumulation in the core, and the core impurity density reaches $ 1.1\times {10}^{16}\;{{\mathrm{m}}}^{-3} $. After introducing neoclassical convection which is outward, it can offset the inward anomalous convection and significantly reduces the W impurity density in the core, thereby significantly reducing the core tungsten impurity density to $ 4.0\times {10}^{15}\;{{\mathrm{m}}}^{-3} $. In addition, the neoclassical convection in the region of $ \rho$ = 0.72–0.90 plays a more important role in reducing the core impurity density. Further analysis of the components of neoclassical convection shows that the Pfirsche-Schlüter (PS) component dominates the neoclassical convection term, which is mainly driven by the ion temperature gradient term. Therefore, experimentally, plasma heating can be used to enhance the temperature gradient and suppress impurity core accumulation.