Search

Article

x

Vol. 74, No. 8 (2025)

2025-04-20
SPECIAL TOPIC—Order tuning in disordered alloys
Kinetic simulation of phase diagram and phase transitions in NiCoCr multi-principal element alloy at high temperature and high pressure
XIONG Haozhi, WANG Yunjiang
2025, 74 (8): 086101. doi: 10.7498/aps.74.20250097
Abstract +
Understanding the phase stability and transformation kinetics of multi-principal element alloys (MPEAs) under extreme conditions is critical for optimizing their performance under extreme conditions such as high-temperature and high-pressure environment. In this work the high pressure-temperature (p-T) phase diagram and solid-liquid transition mechanism of an equiatomic NiCoCr alloy are investigated based on embedded atom method (EAM) potential, through advanced molecular dynamics (MD) simulation combined with enhanced sampling techniques. In order to overcome the timescale limitations of traditional MD in capturing phase transitions as rare events, a hybrid approach integrating well-tempered metadynamics (WTMetaD) and the on-the-fly probability-enhanced sampling with expanded ensembles is used in this work. Collective variables such as enthalpy per atom SH, and two-body entropy SS are used to explore the polymorphic states of the NiCoCr alloy. The crystallinity senv, potential energy U, and volume V are utilized to drive phase transitions, and sampling configurations are performed in the range of 1550–1750 K and 0–10 GPa by using multithermal-multibaric-multiumbrella simulation.Several key results about liquid-solid phase transition in NiCoCr alloy are obtained as follows.1) Phase diagram prediction. NiCoCr alloy exhibits a stable body-centered cubic (BCC) phase under high-pressure condition (e.g. 10 GPa) at elevated temperatures (up to 1750 K), rather than a face-centered cubic stable (FCC) phase at room temperature and ambient pressure. The solid-liquid coexistence line shifts upward with the increase of pressure, raising the melting temperature from ~1400 K (ambient pressure) to about 1750 K (over 10 GPa).2) Free energy landscape. The free energy curves corresponding to different thermodynamic conditions are obtained using reweighting techniques and block averaging methods, which reveal that the increase of pressure and decrease of temperature can reduce the free-energy difference ΔGL→BCC, while simultaneously increasing $ G_{ {\mathrm{BCC}} \to{\mathrm{L}}}^* $ required for melting. The combined effects of these changes enhance the stability of the BCC phase in NiCoCr under high-temperature and high-pressure condition.3) Activation parameters and kinetic mechanism. For the activation parameters of solid-liquid dynamic mechanics, $ S_{{\mathrm{L}} \to {\mathrm{BCC}}}^* $ of NiCoCr alloy decreases with the increase of temperature and the decrease of pressure ( from (–4.32 ± 0.16) J·mol–1·K–1 at 1550 K to (–6.71 ± 0.48) J·mol–1·K–1 at 1750 K, 0 GPa ), and |$ V_{{\mathrm{L}} \to {\mathrm{BCC}}}^* $| increases with temperature increasing and pressure decreasing ( from (–88.21 ± 2.57) Å3 at 0 GPa to (–26.09 ± 6.35) Å3 at 10 GPa, 1600 K). At constant temperature, increasing pressure lowers $S^* $ sensitivity to temperature change, whereas higher temperatures amplify pressure’s role in reducing |$ V_{{\mathrm{L}} \to {\mathrm{BCC}}}^* $|, the change of pressure has no significant effect on $ V_{{\mathrm{BCC}}\to {\mathrm{L}}}^* $.These results demonstrate that the synergistic effects of pressure and temperature on $S^* $ and $V^* $ dictate the phase stability and transformation kinetics of NiCoCr alloys under extreme conditions. The predicted p-T phase diagram and quantitative activation parameters provide critical ideas for designing MPEAs with tailored microstructures for high-pressure applications. Limitations of the EAM potential in describing magnetic interactions and non-equilibrium states are discussed, and the necessity of of future validation through first-principles calculations and high-pressure experiments is emphasized.
Dynamic mechanical properties and deformation mechanism of (NiCoV)95W5 medium entropy alloy
LU Shenghan, CHEN Songyang, CUI Guangpeng, ZHOU Dan, CAI Weijin, SONG Min, WANG Zhangwei
2025, 74 (8): 086103. doi: 10.7498/aps.74.20250141
Abstract +
Medium-entropy alloys (MEAs), renowned for their outstanding strength and ductility, possess great potential for high strain-rate applications. This study focuses on a NiCoV-based MEA system, and proposes a novel alloy design strategy to fabricate the (NiCoV)95W5 alloy by introducing 5% (atomic percent) high-melting-point tungsten through vacuum arc melting coupled with thermomechanical processing. Split Hopkinson pressure bar (SHPB) experiments are conducted to elucidate the dynamic response mechanism and deformation behavior under high strain rates (2000-6000 s–1). The results show that due to severe lattice distortion, the enhanced phonon drag effect at elevated strain rates results in a substantial increase in yield strength from 720 MPa (10–3 s–1) to 1887 MPa (6000 s–1), an increase of 162%, accompanied by a relatively high strain-rate sensitivity (m = 0.42). Microscopic analysis reveals the multi-scale cooperative deformation mechanism of the alloy system under high strain rate. When the strain rate is 2000 s–1, the alloy exhibits a low dislocation density dominated by dislocation planar slip. As the strain rate increases to 4000 s–1, the increased flow stress and deformation promote the proliferation and entanglement of a large number of dislocations into high-density dislocation cells. The accumulation of dislocation stress leads to the coordinated deformation of precipitates and releases stress concentration at the phase interface. When the strain rate further increases to 6000–1, severe plastic deformation will lead to the formation of nanotwins within the matrix, which is the main strain hardening. This study elucidates the dynamic response mechanism of NiCoV MEA mediated by tungsten doping, providing a guidance for designing novel structural materials with excellent dynamic mechanical responses.
SPECIAL TOPIC—Ordering control of amorphous alloys
Computational simulation of atomic strain in body-centered cubic multi-principal element alloys
SONG Qianqian, ZHANG Bozhao, DING Jun
2025, 74 (8): 086102. doi: 10.7498/aps.74.20250128
Abstract +
Multi-principal element alloys (MPEAs), also known as high-entropy alloys (HEAs), are novel materials that have received significant attention due to their exceptional mechanical properties, thermal stability, and resistance to wear and corrosion. These alloys are typically composed of multiple principal elements in near-equal atomic proportions, forming solid solution phases such as face-centered cubic (FCC) or body-centered cubic (BCC) structures. Despite the promising applications, a more in-depth understanding of the atomic-level behavior, particularly, lattice distortion and atomic strain, is essential to better design and optimize these materials in extreme environments. This study focuses on systematically investigating the atomic-scale lattice distortion characteristics and their influence on atomic strain in three representative BCC-based MPEAs: TaWNbMo, TiZrNb, and CoFeNiTi. We utilize molecular dynamics (MD) simulations to explore the local atomic strain distributions in these alloys at various temperatures. Von Mises strain and volumetric strain are employed as key descriptors to quantify the atomic strain, providing a clear representation of how lattice distortion on an atomic scale influences the overall strain behavior. The study specifically addresses the effects of atomic radius differences, chemical short-range ordering, and temperature on the strain characteristics of the alloys. The results obtained indicate that an increase in lattice distortion corresponds to a broader distribution of von Mises strain and volumetric strain, with strain values significantly amplified. More precisely, alloys with larger atomic radius differences exhibit greater volumetric strain, reflecting the influence of atomic size disparity on strain distribution. Furthermore, the formation of chemical short-range order (CSRO) significantly mitigates lattice distortion and atomic strain. This finding highlights the importance of short-range atomic ordering in enhancing the stability of the alloy structures, thus potentially improving their mechanical properties. Temperature effects are also investigated, revealing that elevated temperature induces more intense atomic vibration, which in turn increases the atomic strain. The findings underscore the complex interplay between atomic-scale phenomena and macroscopic mechanical properties, offering new insights into the microscopic mechanical behavior of high-entropy alloys. This study contributes to a better understanding of the underlying mechanisms driving atomic strain and lattice distortion in MPEAs. The results provide valuable theoretical insights that can guide the design of high-performance alloys tailored for high-temperature and extreme environments. By addressing the key factors influencing atomic strain, such as atomic radius, chemical ordering, and temperature, this work lays the foundation for future research aimed at enhancing the mechanical performance of MPEAs in various industrial applications.
GENERAL
Influence of B2S3 additive on [111]-oriented diamond crystal synthesized under high pressure condition
WANG Shuai, KANG Ruwei, LI Yong, XIAO Hongyu, WANG Ying, RAN Maowu, MA Hongan
2025, 74 (8): 080701. doi: 10.7498/aps.74.20250028
Abstract +
Diamond is a kind of extremely functional material, which is widely used in the fields of industry, science and technology, military defense, medical and health, jewelry, and others. However, its application in the semiconductor field is still limited, because its electrical transport performance has not yet met the requirements of semiconductor devices. In order to improve the electrical transport performance of diamond as much as possible, the synthesis of diamond single crystal is studied by adding B2S3 to the synthesis system using temperature gradient growth (TGG) method at a pressure of 6.5 GPa in this work. The growth rate of the synthesized diamond crystal decreases from 2.19 mg/h to 1.26 mg/h, indicating that the growth rate of diamond is dependent not only on the growth driving force, but also on the impurity element in the synthetic cavity. Additionally, with the increase of additive dosage, the color of the synthesized diamond crystal changes from yellow to baby blue . Raman measurement results indicate that the obtained diamond appears as a single sp3 hybrid phase without the sp3 hybrid graphite phase. However, the corresponding Raman characteristic peak of the as-grown diamond crystal is located at about 1331 cm–1 and tends to move towards low wave number. According to Fourier Transform Infrared Spectrometer (FTIR) measurement results, the absorption peaks at 1130 cm–1 and 1344 cm–1 are attributed to nitrogen defects. It is found that the nitrogen defect concentration of the synthesized diamond crystal decreases gradually from about 300×10–6 to 60×10–6. Furthermore, the electrical transport performance of the synthesized diamond is characterized by Hall effect measurement. Diamond has insulating properties due to the absence of any additives in the synthetic cavity. However, the results indicate that when B2S3 is introduced into the synthetic system as additive, there is almost no difference in carrier Hall mobility, but the difference in carrier concentration is as high as two orders of magnitude. Furthermore, the resistivity of the synthesized [111]-oriented diamond crystal decreases to 45.4 Ω·cm, due to the addition of B2S3 to the synthesis system. However, it is worth noting that the resistivity of the diamond crystal synthesized with 0.002 g B2S3 and Ti/Cu additives in the synthesis system drops sharply to 0.43 Ω·cm. Therefore, the nitrogen defects in diamond will have an important effect on its conductivity. It provides an important experimental basis for applying diamond to semiconductor field.
Design of in-orbit calibration scheme for scale factor and center-of-mass deviation of inertial sensor of Taiji program for space gravitational wave detection
LIU Chang, WEI Xiaotong, ZHANG Haoyue, DENG Qiong, LIANG Bo, QIANG Li’E, XU Peng, QI Keqi, WANG Shaoxin
2025, 74 (8): 080401. doi: 10.7498/aps.74.20241376
Abstract +
The Taiji program-frequency gravitational waves. The mission's success hinges on the precise operation of its core payloads, particularly the inertial sensors, which are responsible for measuring the residual acceleration noise of the test masses. The duration of a space-based gravitational wave detection mission is 3 to 5 years. During this period, the shift in the satellite’s center of mass due to propellant consumption and other factors, as well as the drift in the scale factors caused by electronic component aging, will gradually degrade the accuracy of inertial sensor data. Therefore, it is necessary to regularly perform in-orbit calibration of inertial sensor parameters.In this work, we develop a calibration scheme, which actively applies controlled satellite oscillations and is tailored according to the installation layout of the inertial sensors in the Taiji program and the noise models. For the calibration of scale factors, high-precision star sensors are used to measure the satellite attitude signal, which is then combined with the driving voltage data from inertial sensors. By using the linear relationship between these signals, the scale factors are estimated using an extended Kalman Filter. For the calibration of center of mass (CoM) offsets, the calibrated scale factors are utilized, along with the driving voltage data from the front-end electronics of inertial sensors, to derive the test mass’s angular acceleration, linear acceleration, and angular velocity. These parameters are then used to complete the CoM offset calibration according to the dynamic equation.The feasibility of the proposed calibration scheme is validated through a simulation experiment. The results demonstrate that the scale factors of the three axes can be calibrated to relative accuracies of 33×10–6 and 173×10–6, respectively, meeting the requirement within 300×10–6. The CoM deviation are calibrated with accuracies of $ {\delta }_{{\boldsymbol{r}}_{1}}= $[15 μm, 31 μm, 34 μm], $ {\delta }_{{\boldsymbol{r}}_{2}}= $[5 μm, 15 μm, 13 μm], satisfying the 75 μm threshold. These results confirm that the proposed scheme can effectively maintain the inertial sensors’ performance within the required accuracy range.All in all, the calibration scheme developed in this study is crucial for maintaining the high performance of inertial sensors in the Taiji program. By achieving the precise calibration of the scale factors and deviation of center of mass within the required accuracy ranges, the scheme ensures the reliability of inertial sensor data, thereby significantly enhancing the sensitivity of space-based gravitational wave detection, which paves the way for groundbreaking discoveries in astrophysics and cosmology.
Design of quantum approximation optimization algorithm for the maximum cut problem
WANG Yunjiang, XI Huiming, XIAO Zhuoyan, WANG Zengbin, SHI Sha
2025, 74 (8): 080301. doi: 10.7498/aps.74.20241223
Abstract +
The max-cut problem (MCP) is a classic problem in the field of combinatorial optimization and has important applications in various fields, including statistical physics and image processing. However, except for some special cases, the MCP still encounters a non-deterministic polynomial complete problem (an NP-complete problem), and there is currently no known efficient classical algorithm that can solve it in polynomial time. The quantum approximate optimization algorithm (QAOA), as a pivotal algorithm in the noisy intermediate-scale quantum (NISQ) computing era, has shown significant potential for solving the MCP. However, due to the lack of quantum error correction, the reliability of computations in NISQ systems sharply declines as the circuit depth of the algorithm increases. Therefore, designing an efficient, shallow-depth, and low-complexity QAOA for the MCP is a critical challenge in demonstrating the advantages of quantum computing in the NISQ era.In this paper, according to the standard QAOA algorithm, we introduce Pauli Y rotation gates into the target Hamiltonian circuit for the MCP. By enhancing the flexibility of quantum trial functions and improving the efficiency of Hilbert space exploration within a single iteration, we significantly improve the performance of QAOA on the MCP.We conduct extensive numerical simulations using the MindSpore quantum platform, and compare the proposed RY-layer-assisted QAOA with standard QAOA and its existing variants, including MA-QAOA and QAOA+. The experiments are performed on various graph types, including complete graphs, 3-regular graphs, 4-regular graphs, and random graphs with edge probabilities between 0.3 and 0.5. Our results show that the RY-layer-assisted QAOA achieves higher approximation ratios in all graph types, particularly in regular and random graphs, where traditional QAOA variants are difficult to implement. Moreover, the proposed method exhibits strong robustness as the graph size increases, and can maintain high performance even for larger graphs. Importantly, the RY-layer-assisted QAOA requires fewer CNOT gates and has a lower circuit depth than the standard QAOA and its variants, making it more suitable for NISQ devices with limited coherence times and high error rates.In conclusion, the RY-layer-assisted QAOA provides a promising approach for solving MCP in the NISQ era. By improving the approximation ratio while reducing circuit complexity, this method demonstrates significant potential for practical quantum computing applications, thus paving the way for developing more efficient and reliable quantum optimization algorithms.
ELECTROMAGNETISM, OPTICS, ACOUSTICS, HEAT TRANSFER, CLASSICAL MECHANICS, AND FLUID DYNAMICS
Single-photon scattering under control of artificial gauge field
WANG Runting, WANG Xudong, MEI Feng, XIAO Liantuan, JIA Suotang
2025, 74 (8): 084205. doi: 10.7498/aps.74.20250021
Abstract +
The mechanism of controlling single-photon scattering in a hybrid system consisting of superconducting qubits coupled to a Su-Schrieffer-Heeger (SSH) topological photonic lattice is investigated under the influence of an artificial gauge field. This research is driven by the growing interest in the intersection between quantum optics and condensed matter physics, particularly in the field of topological quantum optics, where the robustness of photon transport against defects and impurities can be used for quantum information processing. To achieve this, a theoretical model, which incorporates the phase of the artificial gauge field into the coupling between superconducting qubits and the SSH photonic lattice, is developed in this work. The analytical expressions for the reflection and transmission amplitudes of single photons are derived by using the probability-amplitude method. The results show that the artificial gauge field can effectively control single photon scattering in both the upper energy band and the lower energy band of the SSH lattice, thereby enabling total transmission in the upper band and total reflection in the lower band. This band-dependent scattering behavior exhibits a high degree of symmetry with respect to the lattice momentum and energy bands. Importantly, the reflection coefficient can be made independent of the lattice coupling strength and dependent solely on the topological properties of the lattice. This finding suggests a robust method of detecting topological invariants in photonic lattices. Furthermore, our analysis is extended to various coupling configurations between superconducting qubits and the photonic lattice, highlighting the versatility of the artificial gauge field in manipulating photon transport. These findings not only provide new insights into the control of photon transport in topological photonic lattices, but also open the door to the development of novel quantum optical devices and robust quantum information processing platforms.
A method of monitoring image shift based on global fitting Doppler heterodyne interferometer
WEN Zhenqing, LI Juan, HAO Xiongbo, CHANG Chenguang, LI Hongbo, ZUO Haopu, FU Di, FENG Yutao
2025, 74 (8): 084201. doi: 10.7498/aps.74.20250027
Abstract +
Accurate atmospheric wind field measurements are critical for understanding global climate dynamics and facilitating space exploration. Doppler asymmetric spatial heterodyne interferometer (DASH) is used to measure atmospheric wind speed through detecting the phase changes in interferograms induced by Doppler shifts of airglow emission lines. However, environmental temperature fluctuations and mechanical vibrations often cause imaging plane to shift, thereby introducing phase deviations, and degrading the measurement accuracy. In this study, a novel method of monitoring global fitting-based imaging shift is proposed. By etching periodic notches on the diffraction grating surface, the method models and fits the notch patterns formed on the detector plane to achieve precise imaging shift detection and correction. The optimization of notch signal modeling significantly reduces the number of fitting parameters, thus improving computational efficiency and detection precision. Through extensive simulations, the influences of signal-to-noise ratio (SNR) and model parameter variation on detection accuracy are analyzed. The results indicate that when the SNR exceeds 11, the detection uncertainty is still below 6.5 nm. Sensitivity analysis reveals that the detection error stays within acceptable limits when the variations of notch number and notch width are controlled within 40% and 0.7%, respectively, while the influence of edge smoothness parameter of notch pattern is negligible. To validate the performance of the method, the thermal stability is tested by using a near-infrared DASH prototype. The experimental results demonstrate a strong correlation between interferogram phase shifts, imaging plane shifts, and environmental temperature variations. After applying the proposed correction method, local phase fluctuations in the interferogram are significantly reduced, thus the phase stability is improved. Further, artificially applied imaging shifts are accurately detected with errors consistently below 9.96 nm, thereby confirming the reliability and precision of this method. All in all, the proposed method effectively detects and corrects the imaging plane shifts caused by temperature variations, enhancing interferogram phase stability and ensuring high-precision wind speed measurements. This method provides a robust and computationally efficient solution for reducing imaging shifts in DASH systems, and has great potential applications in atmospheric wind field measurement and space-based observation.
Convolutional network single-pixel imaging with fusion attention mechanism
WANG Xiang, ZHOU Yishen, ZHANG Xuange, CHEN Xihao
2025, 74 (8): 084202. doi: 10.7498/aps.74.20250010
Abstract +
This paper presents a novel convolutional neural network-based single-pixel imaging method that integrates a physics-driven fusion attention mechanism. By incorporating a module that combines both channel attention mechanism and spatial attention mechanism into a randomly initialized convolutional network, the method utilizes the physical model constraints of single-pixel imaging to achieve high-quality image reconstruction. Specifically, the spatial and channel attention mechanism are combined into a single module and introduced into various layers of a multi-scale U-net convolutional network. In the spatial attention mechanism, we extract the attention weight features of each spatial region of the pooled feature map by using convolution. In the channel attention mechanism, we pool the three-dimensional feature map into a single-channel signal and input it into a two-layer fully connected network to obtain the attention weight information for each channel. This approach not only uses the critical weighting information provided by the attention mechanism in the three-dimensional data cube but also fully integrates the powerful feature extraction capabilities of the U-net network across different spatial frequencies. This innovative method can effectively capture image details, suppress background noise, and improve image reconstruction quality. During the experimental phase, we employ the optical path of single-pixel imaging to acquire bucket signals for two target images, "snowflake" and "basket". By inputting any noisy image into a randomly initialized neural network with attention mechanism, and using the mean square error between simulated bucket signal and actual bucket signal, we physically constrain the convergence of the network. Ultimately, we achieve a reconstructed image that adheres to the physical model. The experimental results demonstrate that under low sampling rate conditions, the scheme of integrating the attention mechanism can not only intuitively reconstruct image details better, but also demonstrate significant advantages in quantitative evaluation metrics such as peak signal-to-noise ratio (PSNR) and structural similarity (SSIM), confirming its effectiveness and potential application in single-pixel imaging.
Fast generation algorithm of high-quality holograms based on attention convolutional neural network
SUN Kangsheng, HAN Chao, QIN Haifeng, GU Tao, LI Wei, YU Cheng
2025, 74 (8): 084203. doi: 10.7498/aps.74.20241713
Abstract +
In recent years, with the significant improvement of computer performance, deep learning technology has shown an explosive development trend and has been widely used in various fields. In this context, the computer-generated hologram (CGH) generation algorithm based on deep learning provides a new method for displaying the real-time high-quality holograms. The convolutional neural network is a most typical network structure in deep learning algorithms, which can automatically extract key local features from an image and construct more complex global features through operations such as convolution, pooling and full connectivity. Convolutional neural networks have been widely used in the field of holographic displays due to their powerful feature extraction and generalization abilities. Compared with the traditional iterative algorithm, the CGH algorithm based on deep learning has a significantly improved computing speed, but its image quality still needs further improving. In this paper, an attention convolutional neural network based on the angular spectrum diffraction model is proposed to improve the quality as well as the speed of generating holograms. The whole network consists of real-valued and complex-valued convolutional neural networks: the real-valued network is used for phase prediction, while the complex-valued network is used to predict the complex amplitude of the SLM surface, and the phase of the complex amplitude obtained after prediction is used for holographic coding and numerical reconstruction. An attention mechanism is embedded in the down sampling stage of the phase prediction network to improve the feature extraction capability of the whole algorithm, thus improving the quality of the generated phase-only holograms. An accurate diffraction model of the angular spectrum method is embedded in the whole network to avoid labeling the large-scale datasets, and unsupervised learning is used to train the network. The proposed algorithm can generate high-quality 2K (resolution ratio of 11920×1072) holograms within 0.015 s. The average peak signal-to-noise ratio of the reconstruction images reaches up to 32.12 dB and the average structural similarity index measure of the generated holograms can achieve a value as high as 0.934. Numerical simulations and optical experiments verify the feasibility and effectiveness of the proposed attentional convolutional neural network algorithm based on the diffraction model of angular spectrum method, which provides a powerful help for applying the deep learning theory and algorithm to the field of real-time holographic display.
Enhanced sensing of 3.4 GHz microwave in multi-level Rydberg atomic system
XUE Jingjing, LI Ruonan, HU Xuesong, SUN Peisheng, ZHOU HaiTao, ZHANG Junxiang
2025, 74 (8): 084204. doi: 10.7498/aps.74.20250081
Abstract +
The Rydberg-based microwave detection is an all-optical technology that uses the strong coherent interaction between Rydberg atoms and microwave field. Different from the traditional microwave meter, the Rydberg atomic sensing is a new-type microwave detector that transforms the microwave spectrum into a coherent optical spectrum, and arouses increasingly the interests due to its high sensibility. For this kind of sensor, the coherence effect induced by coupling atoms with microwave plays a key role, and the decoherence may reduce the sensitivity. A multi-level Rydberg atomic scheme with optimized quantum coherence, which enhances both the bandwidth and the sensitivity for 4 GHz microwave sensing, is demonstrated experimentally in this work. The enhanced quantum coherence of Rydberg electromagnetically induced transparency (EIT) and microwave induced Autler-Townes (AT) splitting in EIT windows are shown using optical pumping at D1 line. The enhanced sensitivity at 3.4 GHz with 0.3 GHz bandwidth can be realized, based on the enhanced EIT-AT spectrum. The experimental results show that in the stepped Rydberg EIT system, the spectral width of EIT and microwave field EIT-AT can be narrowed by optical pumping (OP), so the sensitivity of microwave electric field measurement can be improved. After optimizing the EIT amplitude and adding single-frequency microwaves, the sensitivity of the microwave electric field measurement observed by the AT splitting interval is improved by 1.3 times. This work provides a reference for utilizing atomic microwave detection.
Sum sideband effect in hybrid optomechanical system with two-level atom ensemble
LIAO Qinghong, TANG Zhian, AO Jiawen
2025, 74 (8): 084206. doi: 10.7498/aps.74.20241432
Abstract +
Cavity optomechanics, as a cross-discipline between nanophotonics and quantum mechanics, provides a unique platform for investigating optomechanical coupling between photons in microcavities and phonons from mechanical modes. It has a wide range of potential applications in quantum physics, and now it has become a hot topic. A theoretical scheme to enhance the sum sideband generation (SSG) via a two-level atom ensemble is proposed. The effect of the atomic ensemble’s detuning frequency on the efficiency of the SSG is considered by introducing a two-level atom medium. The results indicate that the efficiency of the generating sideband can be significantly enhanced under either red or blue detuning of the atoms, with greater dependence and more pronounced enhancement under the red detuning. In addition, we also consider the effect of pump power, which can effectively enhance the intensity of the output signal by selecting the appropriate pump power. More interestingly, the sensitivity of SSG to atomic detuning also indicates that the precise control of the atomic detuning frequency can achieve the fine-tuning of the SSG process. Furthermore, the cavity-atom coupling strength and atom decay rate are discussed for the transmission characteristics of the sum sideband signals. It is found that the efficiency of SSG can be effectively adjusted by the cavity-atom coupling strength and atom decay rate. The results show that the efficiency of SSG can be significantly improved by optimizing system parameters. The method of enhancing SSG may have potential application prospects in measuring high-precision weak forces and on-chip manipulation of light propagation.
CONDENSED MATTER: ELECTRONIC STRUCTURE, ELECTRICAL, MAGNETIC, AND OPTICAL PROPERTIES
A new magnetic phase transition above Morin temperature in DyFeO3
SU Haobin, ZHENG Shiyun, WANG Ning, ZHU Guofeng, JU Xuewei, HUANG Feng, CAO Yiming, WANG Xiangfeng
2025, 74 (8): 087801. doi: 10.7498/aps.74.20250005
Abstract +
Rare-earth orthoferrites (RFeO3) have received significant attention due to their intricate magnetic interactions and potential applications in ultrafast spintronic devices. Among them, DyFeO3 exhibits rich magnetic phase transitions driven by the interplay between Fe3+ and Dy3+ sublattices. Previous studies mainly focused on temperature-induced spin reorientation near the Morin temperature (TM~50 K), but there has been limited exploration of magnetic phase behavior under external fields above TM. This work aims to systematically investigate the temperature- and magnetic-field-dependent magneto-dynamic properties of a-cut DyFeO3 single crystals, with an emphasis on identifying novel phase transitions and elucidating the underlying mechanisms involving Fe3+-Dy3+ anisotropic exchange interactions. High-quality a-cut DyFeO3 single crystals are grown using the optical floating zone method and characterized by X-ray diffraction (XRD) and Laue diffraction. Time-domain terahertz spectroscopy (THz-TDS) coupled with a superconducting magnet (0–7 T, 1.6–300 K) is employed to probe the ferromagnetic resonance (FM) and antiferromagnetic resonance (AFMR) modes. By analyzing the frequency trends in the spectra, the response of internal magnetic moments to external stimuli can be inferred. In the zero magnetic field experiment, it is found that the temperature induced spin reorientation (Γ4Γ1) occurs at Morin temperature (~50 K) with temperature decreasing. A broadband electromagnetic absorption (0.45–0.9 THz) occurs below 4 K, which is attributed to electromagnons activated by broken inversion symmetry in the Dy3+ antiferromagnetic state. Above the Morin temperature, the absorption spectra of the sample are measured at constant temperatures (70, 77, 90, 100 K) and magnetic fields ranging from 0 to 7 T. The experimental results show that with the increase of magnetic field, a new magnetic phase transition occurs (Γ 4 Γ 24 Γ 2 Γ 24 Γ 2 ), and the critical magnetic field of the phase transition varies with temperature. The phase transitions arise from the competition between external magnetic fields and internal effective fields generated by anisotropic Fe3+-Dy3+ exchange. These findings contribute to the further understanding of the magnetoelectric effects in RFeO3 systems and provide a roadmap for using field-tunable phase transitions to design spin-based devices .
Broadband nonreciprocal transmission tuned by pump-induced magnon modes
CHEN Zhijian, ZHAO Kaixin, WANG Chenxiao, WEI Chunke, YAO Bimu
2025, 74 (8): 087501. doi: 10.7498/aps.74.20241666
Abstract +
Nonreciprocal electromagnetic wave transmission is essential for wireless communication, quantum computing, and radar systems, traditionally relying on breaking time-reversal symmetry through static magnetic fields or structural modifications, which face limitations in tunability and integration. Recent advancements in cavity magnonics, particularly the use of bound states in the continuum (BIC) and pump-induced magnon mode (PIM), have enhanced the nonreciprocal isolation and dynamic control of magnon dynamics. In this study, a novel method to achieve broadband-tunable microwave nonreciprocal isolation is presented by introducing multiple modulated pump signals, thereby extending traditional single-mode magnon-based nonreciprocal transmission to multi-channel and broadband regimes. The core method involves exciting multiple PIMs in a cavity magnonics system and strongly coupling them with BIC to generate hybrid modes with pronounced nonreciprocal characteristics. The experimental setup is comprised of a 1-millimeter-diameter yttrium iron garnet (YIG) sphere positioned at the node of a microwave resonator (central frequency: 2.92 GHz), with pump signals injected through a microwave patch antenna. By dynamically tuning the frequency, power, and number of pump signals, the precise control over the number of nonreciprocal isolation channels and their spectral positions is realized. Notably, the continuous tuning of the nonreciprocal bandwidth is achieved by increasing the number of pump signals from 2 to 5, expanding the isolation bandwidth from 6 MHz to 14 MHz. Furthermore, by tailoring the spectral distribution of pump signals, the system realizes flexible switching between bandpass and band-stop isolation states. Importantly, this method eliminates the need of static magnetic field adjustments or structural reconfiguration, relying solely on coherent microwave-photon interactions to modulate PIM-BIC coupling. Experimental results highlight two key physical outcomes: 1) Extending conventional single-mode magnonic nonreciprocal transmission to multi-channel and broadband-tunable regimes; 2) achieving microwave nonreciprocal control without the need of static magnetic field adjustments or structural reconfiguration. These advances establish a robust platform for designing reconfigurable multi-channel isolators and circulators, which can be directly applied to microwave communication systems, quantum information processing, and radar technologies.
Aging characteristics and mechanism of high-dielectric aromatic polyureas considering trap evolutions
FENG Yang, ZHANG Shuo, ZHOU Bin, LIU Peiyan, YANG Xinru, LI Shengtao
2025, 74 (8): 087701. doi: 10.7498/aps.74.20241579
Abstract +
Dielectric capacitors are essential energy storage devices with high power density. The dielectric films of capacitors will age at working temperatures, a leading their performance to degrade. Polyurea (PU) is a potential working dielectric for capacitors with high energy density and low dielectric loss. However, the aging characteristics and underlying mechanism of PU have not been discussed. Considering the operating temperature of commercial dielectric capacitors, the aging characteristics of PU are investigate by being exposed to 80 °C for different durations. Compared with dielectric constant, breakdown strength changes significantly with aging time which can be used as a characteristic parameter to evaluate the aging degree of PU. Combining experimental method and simulation method, the correlation between molecular structure, trap properties and breakdown strength during thermo-oxidative aging is studied and established. The results show that the thermal-oxidative aging of PU can be divided into three stages. In the early stage of aging, the bridging effect of oxygen promotes the order arrangement of molecular chains. In this stage, the molecular chain spacing decreases, but the H-bonding interaction between adjacent urea groups is enhanced slightly as shown in Fig. (a). As a result, the dielectric constant decreases, while the breakdown strength is almost unchanged. In the middle stage of aging, ether bonds break and inducing the formation of biphenyl structures and resulting in a disordered structure as illustrated in Fig. (b). The enhanced mobility effect increases the dielectric constant. Meanwhile, the biphenyl structures deepen the trap depth, resulting in the decrease of carrier mobility and the increase of the breakdown strength. In the late stage of aging, oxygen promotes the decomposition of urea groups, thereby reducing the number of urea groups that form deep traps. At the same time, the main chain undergoes breakage, releasing small molecules such as CO2 and H2O, which is shown in Fig. (c). These factors collectively lead to a significant reduction in the breakdown strength of PU. In addition, the variations of dielectric constant, breakdown strength and energy density in the three stages are summarized in Fig. (d).
A novel LDMOS with dual-drift region and dual-conduction path with ultra-low specific on-resistance
DUAN Baoxing, REN Yuhao, TANG Chunping, YANG Yintang
2025, 74 (8): 087301. doi: 10.7498/aps.74.20241554
Abstract +
In order to improve the contradictory between specific on-resistance (Ron,sp) and breakdown voltage (BV) of lateral double-diffused metal oxide semiconductor (LDMOS) and enhance the turn-off characteristic, this paper proposes a novel LDMOS device with dual-drift regions and dual-conduction paths, which achieves an ultra-low Ron,sp. The key feature of the proposed device is the introduction of a dual-drift region structure with alternating P-type and N-type regions, combined with planar and trench gates to control the P-type and N-type drift regions, respectively. This configuration enables the formation of two independent electron conduction paths within the drift region. When a positive voltage is applied to the planar gate, a voltage difference is generated between the surface of the P-type drift region and the body of device’s drift. Therefore, under the influence of the voltage difference, the electrons are pulled to the surface of the P-type drift region to invert and form a high-density electron inversion layer that connects the channel and the N+ drain, significantly increasing the electron density during conduction and reducing the Ron,sp. The introduction of the trench gate provides an additional electron disappearance path, which shortens the device's turn-off time (toff). Furthermore, the introduction of the P-type drift region facilitates the recombination of electrons with holes within the P-type drift region, accelerating the electron disappearance process and further reducing the device’s toff. Furthermore, the proposed device exhibits a more uniform electric field distribution and higher voltage capability is due to the P+N-N+P+ structure adopted in the PolySi-top layer. During the off-state, both the P+N- junctions and the N+P+ junctions generate electric field peaks at the interfaces. These peaks modulate the electric field distribution across the surface of the drift region. Simulation results indicate that at the BV with a level of 200V, the proposed LDMOS exhibits an Ron,sp of 3.43 mΩ·cm² and a toff of 9 ns. Compared with conventional LDMOS devices, the proposed LDMOS possesses a 90% reduction in Ron,sp and an 11.6% decrease in toff. The proposed device not only achieves an excellent trade-off between Ron,sp and BV but also shortens the toff, demonstrating that the device achieves superior performance.
REVIEW
Advances in single crystal growth methods for novel unconventional superconductor UTe2
XUE Ziwei, YUAN Dengpeng, TAN Shiyong
2025, 74 (8): 087401. doi: 10.7498/aps.74.20241778
Abstract +
Heavy fermion compound UTe2, as a recently discovered unconventional superconductor, has received significant attention due to its potential spin-triplet superconducting pairing, high-field re-entrant superconducting phases, and unique quantum critical characteristics. However, experimental results of this system show significant changes and discrepancies, primarily due to difference in sample quality. The key unresolved issues include whether the system exhibits multi-component superconducting order parameters, whether time-reversal symmetry is spontaneously broken, and whether multiple field-induced superconducting phases share a common origin. These unsolved issues hinder an in-depth understanding of the intrinsic superconducting pairing mechanism in the UTe2 system.This paper reviews recent advances in single-crystal growth methods for UTe2, including chemical vapor transport (CVT), Te-flux, molten salt flux (MSF), and molten salt flux liquid transport (MSFLT). We systematically analyze how growth conditions influence superconductivity and crystal quality. Although the CVT method was initially employed in UTe2 studies, the samples grown by this method exhibit poor quality and significant compositional inhomogeneity, even in individual samples. Consequently, the CVT method has been progressively supplanted by the recently developed MSF method. In contrast, the MSF method and MSFLT method yield high-quality UTe2 single crystals with Tc achieving a value as high as 2.1 K and residual resistivity ratio (RRR) reaching up to 1000; however, the sample sizes are smaller than those grown by the CVT and Te-flux methods. Notably, MSF-grown samples occasionally contain magnetic impurities such as U7Te12, so careful screening is required in the sample collection process. The MSFLT combines the advantages of CVT and MSF methods to grow high-quality UTe2 single crystals while producing larger sample sizes than MSF. Our research findings highlight the importance of optimizing growth parameters such as Te/U ratio, temperature gradient, and cooling rate. For instance, lower growth temperature and precise control of the Te/U ratio can significantly enhance Tc and sample quality. Several controversies have been identified regarding high-quality MSF and MSFLT samples, including clarifying the single-component nature of the superconducting order parameter and confirming the absence of time-reversal symmetry breaking in optimized samples.This review underscores the pivotal role of advanced single-crystal growth techniques in advancing the study of UTe2. Future research should focus on utilizing these high-quality UTe2 samples grown by MSF and MSFLT methods to accurately determine superconducting order parameters, elucidate mechanisms behind high-field re-entrant superconducting phases, and explore topological properties, such as potential Majorana fermions. These efforts will deepen our understanding of unconventional superconductivity, spin fluctuations, and quantum critical phenomena in the UTe2 system.
SPECIAL TOPIC—Precision spectroscopy of few-electron atoms and molecules
Applications of B-spline method in precise calculation of structure of few-electron atoms
ZHANG Yonghui, SHI Tingyun, TANG Liyan
2025, 74 (8): 083101. doi: 10.7498/aps.74.20241728
Abstract +
The precise spectra of few-electron atoms plays a pivotal role in advancing fundamental physics, including the verification of quantum electrodynamics (QED) theory, the determination of the fine-structure constants, and the exploration of nuclear properties. With the rapid development of precision measurement techniques, the demand for atomic structure data has evolved from simply confirming existence to pursuing unprecedented accuracy. To meet the growing needs for precision spectroscopy experiments, we develop a series of high-precision theoretical methods based on B-spline basis sets, such as the non-relativistic configuration interaction (B-NRCI) method, the correlated B-spline basis functions (C-BSBFs) method, and the relativistic configuration interaction (B-RCI) method. These methods use the unique properties of B-spline functions, such as locality, completeness, and numerical stability, to accurately solve the Schrödinger and Dirac equations for few-electron atoms. Our methods yield significant results, particularly for helium and helium-like ions. Using these methods, we obtain accurate energies, polarizabilities, tune-out wavelengths, and magic wavelengths. Specifically, we achieve high-precision measurements of the energy spectra of helium, providing vital theoretical support for conducting related experimental researches. Additionally, we make high-precision theoretical predictions of tune-out wavelengths, paving the way for new tests of QED theory. Furthermore, we propose effective theoretical schemes to suppress Stark shifts, thereby facilitating high-precision spectroscopy experiments of helium. The B-spline-basis methods reviewed in this paper prove exceptionally effective in high-precision calculations for few-electron atoms. These methods not only provide crucial theoretical support for precision spectroscopy experiments but also pave the new way for testing QED. Their ability to handle large-scale configuration interactions and incorporate relativistic and QED corrections makes them versatile tools for advancing atomic physics research. In the future, the high-precision theoretical methods based on B-spline basis sets are expected to be extended to cutting-edge fields, such as quantum state manipulation, determination of nuclear structure properties, formation of ultracold molecules, and exploration of new physics, thus continuously promoting the progress of precision measurement physics.
SPECIAL TOPIC—Correlated electron materials and scattering spectroscopy
Research progress of resonant X-ray scattering of charge order in cuprate superconductors
CHAN Ying, YAN Yujie, WU Yuetong, WANG Qisi
2025, 74 (8): 087402. doi: 10.7498/aps.74.20241402
Abstract +
Unconventional superconductivity often competes or coexists with a variety of complex material states. In cuprate superconductors, there exist states including spin order, charge order, the pseudogap state, and the strange metal phase. A comprehensive understanding of their relationship is fundamental to establishing the mechanism of high-temperature superconductivity. Spin dynamics in cuprates has been extensively investigated using inelastic neutron scattering, but charge correlations remain far less understood. The latest development of resonant X-ray scattering (RXS) has been able to detect charge correlations with unprecedented sensitivity. A series of RXS studies have revealed that there universally exist the charge correlations in cuprate materials, which covers a wide range of the phase diagram. Resonant inelastic X-ray scattering (RIXS) experiments further show the dynamical behaviors of charge order. These findings highlight the important influence of charge correlations on the properties of cuprates. In this paper, we review the latest research progress in the charge order in cuprates by using RXS, with a particular emphasis on RIXS experiments. Our focus is placed on recent works on dynamical charge correlations at high temperatures as well as uniaxial strain tuning of charge order. We discuss topics including the underlying interactions, microscopic structure and symmetries, and the possible influence of charge order on both the superconducting and normal states.
Laser trapping and manipulation of micro/nano-objects on polymer substrates
YIN Yue, DOU Lin, SHEN Tianci, LIU Jiatong, GU Fuxing
2025, 74 (8): 088703. doi: 10.7498/aps.74.20241654
Abstract +
Polymer substrates break through the limitations of rigid planar substrates in spatial deformation scenarios and can be combined with photolithography to fabricate complex, three-dimensional irregular polymer structures. Photothermal-shock tweezer is a laser trapping technique based on the photothermal shock effect. Photothermal-shock tweezer uses pulsed laser induced transient photothermal shock to generate micro-newton-scale thermomechanical strain gradient force, enabling the trapping and manipulation of micro/nano-objects at solid interfaces. Integrating this technique with polymer substrates can meet the demands of new application scenarios. In this work, commonly employed polymethyl methacrylate (PMMA) and negative photoresist (SU-8) are used as polymer substrates, on which SiO2 nanofilms are prepared using the sol-gel method. This method effectively mitigates thermal damage caused by photothermal shock effects, enabling laser trapping and manipulation of micro/nano-objects.The SiO2 nanofilms, characterized by low thermal conductivity, effectively inhibit heat transfer. The nanofilm fabrication technique utilized in this study enables the synthesizing of large-area SiO2 nanofilms with large-area coverage, low surface roughness (Rq ~ 320 pm) and uniform thickness, making them broadly applicable to flexible polymer substrates and irregular structures. Direct contact between the polymer layer and micro/nano-objects during manipulating the photothermal shock tweezer can induce irreversible substrate degradation due to transient photothermal shock effects. Experimental results demonstrate that depositing an SiO2 nanofilm thicker than 110 nm on the polymer substrate can significantly enhance thermal insulation and protection, effectively mitigating laser-induced damage under typical optical manipulation conditions.Additionally, by analyzing the temperature field distribution of the gold nanosheet, PMMA substrate, and SiO2 nanofilm during a single photothermal shock trapping of a gold nanosheet, it is found that the SiO2 nanofilm can reduce the PMMA surface temperature by at least 111 ℃ and delay the time for PMMA to reach its peak temperature by 13.2 ns compared with the the gold nanosheet. The experimental results expand the environmental medium for laser trapping of objects, providing new possibilities for applications in micro/nano-manipulation, micro/nanorobotics, and micro/nano-optoelectronic devices.
High-temperature oxidation and crystallization mechanism of Fe73.5Si13.5B9Cu1Nb3 amorphous alloy
ZHANG Xiang, SONG Yingjie, LIU Haishun, XIONG Xiang, YANG Weiming, HAN Chenkang
2025, 74 (8): 088101. doi: 10.7498/aps.74.20250112
Abstract +
Fe-based amorphous alloys are widely used in electronic devices such as high-frequency transformers and choke cores due to their low coercivity, low loss, and high saturation magnetic induction intensity. However, these alloys have a relatively low crystallization temperature and are prone to oxidation, which limits their applications in high-temperature environments. The addition of copper and niobium elements can suppress the growth of crystal nuclei and improve thermal stability. However, the influences on the alloy's high-temperature oxidation resistance and structural evolution are still unclear. In this work, static air oxidation is used to investigate the microstructure evolution of Fe73.5Si13.5B9Cu1Nb3 amorphous alloy after high-temperature oxidation and its influence on magnetic properties. Besides, long-time oxidation, say, 3000 hours or longer at 500 ℃, is generally hard to perform in the laboratory. Thus, the Van’t Hoff’s rule is used to evaluate outcomes under the condition of the long-time and relatively low-temperature oxidation through using rapid high-temperature oxidation. Based on Van’t Hoff’s rule, the oxidation at 650 ℃ for 5 min will show similar or more severe oxidation effects on the microstructure of Fe73.5Si13.5B9Cu1Nb3 alloy after oxidation at 500 ℃ for 2730 h. The microstructure evolution reveals that silicon and niobium in this alloy will quickly diffuse toward the sample surface during oxidation at 650 ℃, and these two elements will form a dense layer to impede oxygen diffusion. Meanwhile, an α-Fe(Si) phase, mainly composed of iron elements, will be generated in the alloy, with its grain size slowly increasing in the oxidation process. Thermodynamic analysis indicates that the segregation of silicon and niobium can preserve the thermodynamic stability of the alloy system during oxidation and suppress the formation of intermetallic compounds during crystallization. The magnetic hysteresis loop results show that the coercivity of Fe73.5Si13.5B9Cu1Nb3 alloy after 5-min oxidation at 650 ℃ will stay at approximately 0.3 Oe, suggesting that the Fe73.5Si13.5B9Cu1Nb3 alloy may be a candidate for operating at 500 ℃ for more than 2700 h. Subsequently, its coercivity gradually increases to 61 Oe as the oxidation time rises to 0.5 h, while its saturation magnetic induction intensity remains unchanged (~140 emu/g).
An algorithm of reconstructing phaseless radiation source based on singular value decomposition regularization and fast iterative shrinkage-thresholding algorithm
DENG Dianjun, LI Yan
2025, 74 (8): 088701. doi: 10.7498/aps.74.20250062
Abstract +
An algorithm of reconstructing phaseless radiation source based on singular value decomposition (SVD) regularization and fast iterative shrinkage-thresholding algorithm (FISTA) is proposed in this work, aiming at efficiently identifying electromagnetic interference (EMI) sources in integrated circuits (ICs). The method acquires electromagnetic field data through near-field scanning and reconstructs an equivalent dipole array on the surface of the radiation source by using the source reconstruction method (SRM). In the reconstruction process, the SVD regularization term enhances the algorithm's stability and noise resistance, while the FISTA accelerates the convergence speed.In order to validate the effectiveness of the proposed method, dipole array reconstruction is first performed using near-field data at a height of 5 mm for a patch antenna simulation model, followed by analyzing the magnetic field data at a 10 mm validation plane. At the 35th iteration, the total relative error of the reconstruction is 1.21%. The influence of the regularization parameter α on the result is then investigated, and it is found that when α = 0.05 the error is minimized. The method is also tested under different Gaussian white noise conditions, and the relative error is kept below 5%, which demonstrates strong robustness.Finally, the experiments on chips are conducted to verify the method. The proposed method converges stably within 35 iterations, with a relative error of 2.3% in the reconstruction results. The proposed method reduces the total iteration time to 61.7% of the single-layer phaseless interpolation algorithm, while achieving a 52% lower relative error than the double-layer phasless iteration algorithm. The experimental results show that the proposed method can reconstruct phaseless radiation source efficiently and accurately, and has good noise robustness, which is suitable for EMI analysis in ICs.
Effects of Al content on stability, electronic structure, and Li-ion diffusion properties of Li1+xAlxTi2–x(PO4)3 surface
LI Mei, YAN Yi, LAN Wenxin, SUN Baozhen, WU Musheng, XU Bo, OUYANG Chuying
2025, 74 (8): 088201. doi: 10.7498/aps.74.20250016
Abstract +
NASICON-type Li1+xAlxTi2–x(PO4)3 (LATP), as a promising solid-state electrolyte for lithium-ion batteries, has received significant attention due to its simple preparation, low material cost, and good stability in water and air, but the formation of lithium dendrite greatly limits the applications. To elucidate the source of formation of lithium dendrite, in this study, the effects of Al content on the stability, electronic and Li+ mobility properties of the LATP surface with three Al doping concentrations (2AlTi, 4AlTi, 6AlTi) are investigated by combining first-principles calculations and molecular dynamics simulations. The LiTi2(PO4)3 (LTP) surface is also considered for comparison. The results indicate that the (012) surface terminated with Li atoms is the most stable facet. Further, the surface energy of LATP(012) decreases from 0.68 J/m2 to 0.43 J/m2 with the increase of Al content, suggesting that Al doping can effectively improve the stability of the LATP(012) surface. Electronic structure analysis reveals that the surface of LTP(012) retains the semiconductor properties consistent with the bulk phase, whereas the LATP(012) surface exhibits metallicity, which provides an electron pathway for forming the metallic Li . Consequently, the metallic characteristic of the LATP(012) surface is a reason for its lithium dendrite growth. For the Li+ transport properties, two different migration modes: vacancy migration and interstitial migration, are included. When Li+ migrates within the outermost surface, the migration barrier via vacancy is 1.67/1.69 eV for the LTP/LATP (012) surface, while the migration barrier via interstitial is 1.16 eV for LTP(012) and decreases from 1.31 to 0.87 eV with the increase of Al content for LATP(012). Obviously, within the outermost surface, Al doping can reduce the migration barrier of Li+. When Al doping concentration is 6AlTi, the migration barrier is lowest (0.87 eV). Nevertheless, the lowest migration barrier (0.87 eV) for Li+ on the LATP surface is significantly higher than its bulk minimum value of 0.34 eV. When Li+ migrates from the subsurface layer to the outermost surface, the migration barrier is 2.76 eV for LTP(012) and 2.05 eV, 3.20 eV, and 3.06 eV for LATP(012) with 2AlTi, 4AlTi, and 6AlTi content, respectively. All these migration barriers are greater than 2.00 eV, which prevents Li+ from migrating from the subsurface layer to the outermost surface for both LTP and LATP surfaces. Hence, the slow Li+ migration represents another important factor contributing to lithium dendrite growth on the LATP surface. Fortunately, increasing the Al doping concentration can reduce the migration barrier of Li+ and thus enhance its diffusion performance on the LATP surface. Molecular dynamics simulations further reveal that the diffusion behavior of Li+ on the LATP surface is influenced by a combination of factors, including Al content, Li+ occupancy, and ambient temperature. In particular, LATP(012)/6AlTi, LATP(012)/4AlTi, and LATP(012)/2AlTi possess their highest Li+ diffusion coefficients at 900 K, 1100 K, and 1300 K, respectively. Besides, Li+ near the Al doping site is easier to diffuse on the LATP(012) surface. Thus, our study indicates that by changing Al content, Li+ occupation positions, and the temperature, the Li+ diffusion performance of LATP(012) can be effectively modified, thereby suppressing the formation of lithium dendrites on the LATP(012) surface.
Capture of Rayleigh particles by focused double-ring perfect vortex beams
SHI Pengbo, YANG Chen, YIN Xiaojin, WANG Jie
2025, 74 (8): 088702. doi: 10.7498/aps.74.20241768
Abstract +
The double-ring perfect vortex beam (DR-PVB) is generated through the superposition of two concentric perfect vortex beams (PVBs). In this work, firstly, the intensity and phase distribution of the DR-PVB in the source plane are studied. Secondly, based on the Huygens-Fresnel principle and the Collins formula, the intensity distribution of the DR-PVB after being focused by an ABCD optical system that includes a focusing lens is obtained. The results indicate that the intensity distribution of the focused beam is consistent with the interference pattern of two Bessel Gaussian beams. Furthermore, the number of spots in the focused intensity distribution is a multiple of the absolute value of the difference in topological charges between two PVBs. On the other hand, the overall size of the light beam can be adjusted by changing the focal length of the lens. Thirdly, the optical radiation force, exerted by the focused DR-PVB, on Rayleigh particles with different refractive indices, silica and bubbles, are analyzed, respectively. The results show that the focused DR-PVB can capture both high and low refractive index particles in water. In addition, by comparing the focused DR-PVBs under different radius combinations, it found that the light intensity distribution can be changed with the beam radius, which leads the position and quantity of the captured particles to change. This result provides a new idea for adjusting the capture of particles in future experiments. Finally, the gradient forces, scattering, and Brownian forces acting on the particles in the x, y, and z directions are analyzed, respectively. Based on our analysis, the condition for stable particle capture, where the gradient force must overcome the effects of Brownian motion and scattering forces, is established. Therefore, the theoretical size range of particles that can be captured by the focused DR-PVB is determined. Compared with other beams, such as Airy beams and Bessel beams, the focused DR-PVB can be modulated by changing the topological charges of the two PVBs, making it possible to capture multiple particles. These results have potential applications in optical manipulation.
NUCLEAR PHYSICS
Outliers identification of experimental (γ, n) reaction cross section via variational autoencoder
XIE Jinchen, TAO Xi, XU Ruirui, TIAN Yuan, XING Kang, GE Zhigang, NIU Yifei
2025, 74 (8): 082501. doi: 10.7498/aps.74.20241775
Abstract +
The (γ, n) cross-section is important in nuclear engineering transport calculations. The measurements of the (γ, n) reaction for some isotopes show significant discrepancies among different laboratories. Since the analysis of experimental data is the primary task in the evaluation of nuclear data, identifying the measured outlier data is crucial for improving the quality of nuclear data. Therefore, a variational autoencoder (VAE) method is used in this work to analyze experimental measurements of (γ, n) cross sections for nuclear mass ranging from 29 to 207, in order to provide more reliable experimental information for evaluating nuclear data.According to the proton number Z and nuclear mass A, we design a variational autoencoder network for outlier identification in the measurement of (γ, n). The silhouette coefficient method and K-means algorithm are used to cluster the latent variables of VAE. Subsequently, the experimental data with and without the outliers are compared with those from the IAEA-2019-PD to assess the VAE in its application to the evaluation of photoneutron measurements.The results demonstrate that the VAE can effectively identify outliers in the measurements of (γ, n). After excluding outliers, the (γ, n) cross-section for $^{54}\text{Fe}$, $^{63}\text{Cu}$, $^{181}\text{Ta}$, $^{206}\text{Pb}$, and $^{207}\text{Pb}$ showed higher consistency with the IAEA-2019-PD evaluation results. However, $^{29}\text{Si}$ and $^{141}\text{Pr}$ still deviate from the IAEA-2019-PD evaluation results, therefore requiring more analyses of the measurements themselves in future.The VAE method can effectively identify outliers and extract the latent structures in experimental data of (γ, n) reaction. It provides more reliable experimental information for evaluating nuclear data and validating the potential application of this method in nuclear data research. However, the generalizability of VAE method still needs further developing, especially in addressing the issues of uneven energy distribution for various measurements.
CONDENSED MATTER: STRUCTURAL, MECHANICAL, AND THERMAL PROPERTIES
Effect of ZrSiO4 on far-infrared radiation characteristics of Pb3O4/Mg/PTFE infrared decoy
WANG Bing, CHEN Zongsheng, SHI Jiaming, XU Hexiu
2025, 74 (8): 084401. doi: 10.7498/aps.74.20241048
Abstract +
Traditional composite infrared decoy, magnesium/teflon (Mg/PTFE), has been widely used in countering infrared guided weapons since its advent. However, with the development of infrared guidance technology, its drawbacks such as insufficient far-infrared radiation and high combustion temperature emerge, making it difficult to counter novel infrared guided weapons. To address this issue, a strategy of utilizing zirconium silicate (ZrSiO4) as an additive is proposed to improve the infrared radiation of infrared decoy. Therein, seven formulations with different ratios of ZrSiO4 are designed based on the basic formula of trilead tetraoxide/ magnesium/teflon (Pb3O4/Mg/PTFE) mixed powder. And the effect of ZrSiO4 serving as an additive on the performance of Pb3O4/Mg/PTFE infrared decoy is analyzed through experiments. First, initial experiments are conducted on the thermal decomposition characteristics of the basic formula (ZrSiO4 addition ratio is 0%) and its variant counterpart with 12% ZrSiO4. Subsequently, the combustion behaviors of the compacted formulation samples are examined using an infrared thermal imager operating within the 7.5–14 μm range, subsequently, the combution time, combution temperature, burning rate, radiation area, radiance, and radiation intensity of individual samples are computed. These results show that incorporation of ZrSiO4 reduces the intensity of the primary exothermic peak during the reactions with a mixed infrared decoy agent, yielding suboptimal thermal efficiencies. Furthermore, the combustion durations of the samples progressively increase with ZrSiO4 addition increasing, accompanied by consistent reductions in their combustion temperatures. Specifically, the sample reaction time peaks at 3.73 s at a ZrSiO4 addition ratio of 18%, while the combution temperature drops to a minimum value of 765.46 ℃. Moreover, the far-infrared radiance and radiation intensity demonstrate an initial-increase-then-decrease trend with ZrSiO4 addition increasing, thereby achieving the maximum values of 2461 W/(m2·sr) and 142 W/sr, respectively at a ZrSiO4 addition ratio of 6%. Furthermore, the far-infrared radiance and radiation intensity of the base formulation are enhanced when ZrSiO4 addition ratios are kept within 18% and 9% respectively. Based on the comprehensive analysis of the experimental data and considering the requirements for the infrared decoy in practical applications, a formulation with a ZrSiO4 addition ratio of 6% is adopted as an improved formulation for the Pb3O4/Mg/PTFE infrared decoy.
Thin film growth and electronic transmission of picene molecules on Cd(0001) surface
YU Tao, FU Shirun, LU Hongping, WANG Junzhong
2025, 74 (8): 086801. doi: 10.7498/aps.74.20250076
Abstract +
The thin films of picene on the Cd(0001) surface are investigated by scanning tunneling microscopy (STM) in this work. Compared with conventional noble metal substrates such as Au, Ag, and Cu, Cd(0001) has low electronegativity and small work function , which can effectively weaken the molecule-substrate interactions, thereby promoting the intermolecular van der Waals attraction. The experiments are conducted by ultrahigh vacuum low-temperature STM combined with a molecular beam epitaxy system. The crystalline Cd(0001) films are grown on a surface of Si(111)-7×7 by depositing 15-20 monolayers of Cd atoms, followed by being annealed. The picene molecules are deposited on the Cd(0001) surface at 100–120 K, where monolayer (ML) is defined as the critical coverage before second-layer nucleation. All STM measurements are acquired in constant-current mode.It is observed that, in the submonolayer regime, the picene molecules occupy the entire substrate surface and form disordered two-dimensional molecular gas, indicating the existence of electrostatic repulsive interaction among picene molecules. With the coverage increasing, the first layer of molecules undergoes the disorder-order transition, forming the parallel array of molecular stripes of flat-lying molecules. The high-resolution STM images show that the building block of molecular stripes is a picene dimer with the opposite dipole moments. More importantly, under a specific bias voltage, the first layer of molecular stripes exhibits electronic transmission: not only can the underlying Cd substrate atoms be observed, but also the standing waves of scattered electrons can be observed near the defects. When the coverage exceeds 1.0 ML, the second picene layer also forms the stripe array composed of picene dimers of a flat-lying and a side-on molecules, similar to the (110) plane in picene crystals. The above results show that the electrons from the quantum-well states of Cd (0001) thin film have very strong penetration ability, and their vertical tunneling length reaches to the distance of two molecular layers.
ATOMIC AND MOLECULAR PHYSICS
Spectral characteristics of Floquet-electromagnetically induced transparency dressed by radio frequency field in a direct current electric field
DUAN Haonan, JI Zhonghua, LIU Weixin, SU Dianqiang, LI Jingkuan, ZHAO Yanting
2025, 74 (8): 083201. doi: 10.7498/aps.74.20250052
Abstract +
A Rydberg atom is a special type of atom characterized by a high principal quantum number. Electric field sensors based on Rydberg atoms have received widespread attention due to their high polarizability. However, there is currently little research on the use of Rydberg atoms for direct current (DC) or low-frequency electric fields, mainly due to the shielding effect of atomic vapor cells in low-frequency electric fields, which makes accurate measurement of the electric fields extremely challenging.In this paper, we construct a Rydberg ladder configuration by using probe laser at 852 nm and coupling laser at 510 nm in a room-temperature cesium vapor cell with integrated electrode plates, thereby enabling the realizing of a Floquet-EIT (electromagnetically induced transparency) spectrum dressed by a radio frequency (RF) field in the presence of a DC electric field. We further study the influence of DC electric field on spectral characteristic. Experimentally, it is observed that when only the RF electric field is applied, the EIT spectrum displays solely even-order sidebands. Furthermore, when both the RF field and the DC electric field are simultaneously present, the first-order sideband signal of the Floquet-EIT is observed. As the intensity of the DC electric field increases, the amplitude of the first-order sideband gradually increases. However, increasing the DC electric field to a sufficient magnitude induces sideband interference, which leads the sideband amplitudes to decrease. Furthermore, increasing the RF frequency can alleviate the interference effects induced by the DC electric field on the first-order sideband. Finally, comparing the relative standard deviation of the sideband amplitudes of the Floquet-EIT spectra with the frequency shifts of the DC-Stark spectra under weak DC electric fields, we find that the measurement accuracy of the former is significantly superior to that of the latter.This work makes use of a Cs atomic vapor cell with an integrated electrode to avoid shielding effects. By observing Floquet-EIT spectra, the response of the spectra to DC electric fields is investigated. This experiment provides novel insights into the quantum sensing measurements of DC and low-frequency electric fields.
THE PHYSICS OF ELEMENTARY PARTICLES AND FIELDS
Zeno’s paradox and black hole information loss problem
Ge Xianhui
2025, 74 (8): 081101. doi: 10.7498/aps.74.20241751
Abstract +
The black hole information paradox—a fundamental conflict between quantum unitarity and semi-classical gravity—remains unresolved within the framework of quantum gravity. This work builds a conceptual bridge between this modern dilemma and Zeno’s ancient paradox of motion, emphasizing that they all rely on limit processes to reconcile infinite divisibility with finite physical outcomes. Although previous researches have established information conservation in black hole evaporation by copying wormholes, this study transcends mere unitarity verification to answer a deeper question: How do we quantify the statistical evolution and inherent irreversibility of Hawking radiation? Methods : i) Gravitational path integral formalism : We compute Rényi entropies $ S_n = \dfrac{1}{1-n} \ln \text{Tr}(\rho^n) $ by analytic continuation $( n \to 1 )$, incorporating both disconnected geometries (Hawking saddles) and connected replica wormhole saddles. ii) JT gravity model : A two-dimensional Jackiw-Teitelboim (JT) gravity coupled to a conformal bath is analyzed, with End-of-the-World (EOW) branes modeling early Hawking radiation states. iii) Modular thermodynamics : Modular entropy $ S_m $ and entanglement capacity $ C_n $ are defined as:           $ S_m = \log\text{Tr}(\rho^n) - n\partial_n\log\text{Tr}(\rho^n), \quad C_n = n^2\partial_n^2\log\text{Tr}(\rho^n),$mirroring thermodynamic entropy and heat capacity. Key results : i) Page curve restoration : The von Neumann entropy transitions from linear growth (Hawking phase) to zero (wormhole-dominated phase) at the Page time: $ t_{\text{Page}} \propto {S_{\text{BH}}}/({n\kappa c}), \quad S_{\text{BH}} = {A}/({4G_N}), $indicating faster entropy reduction for higher n. ii) Replica wormhole mechanism : The gravitational path integral reveals a topological transition:           $ \text{Tr}(\rho^n) \approx \dfrac{kZ_1^n + k^n Z_n}{(kZ_1)^n} \quad \Rightarrow \quad \text{Dominance shift at } k \sim {\mathrm{e}}^{S_0}$.Quantum non-cloning constraints: A reversibility paradox arises: while n-copies -copies of a known state allow full reconstruction, the non-cloning theorem forbids replication of unknown states. This necessitates the a priori existence of multiple replicas to ensure consistency. iii) Connection to Zeno’s paradox : The resolution of Zeno’s paradox through Newtonian limits $( v = \lim\limits_{\Delta t \to 0} \Delta x / \Delta t $) parall) parallels the replica trick’s role in black hole physics: iv) Infinite-to-finite mapping : Both employ limit operations to compress divergent processes (infinite steps or entropy growth) into finite observables. Discrete vs. continuous : Challenge classical spacetime continuity, hinting at holographic or discrete quantum spacetime. Irreversibility and generalized second law : The relative entropy $ S(\rho \| \sigma) $, defined as:               $ S(\rho \| \sigma) = \text{Tr}(\rho \log \rho) - \text{Tr}(\rho \log \sigma),$monotonically increases during the Hawking-to-wormhole transition, establishing statistical irreversibility in modular space. This generalized second law persists despite microscopic unitarity, analogous to thermodynamic irreversibility in closed systems with the formula: $ \left(S_n-S_{\rho}\right)-n\left(\left\langle H\right\rangle_n-\left\langle H\right\rangle_{\rho}\right)\geqslant 0 $. Conclusion :Our results address the black hole information paradox and the irreversible problem through three pillars: i) Unitarity : Replica wormholes enforce unitarity via non-perturbative spacetime topology changes. ii) Modular thermodynamics : A framework linking entanglement entropy, relative entropy, and irreversibility. iii) Replica priority principle : Multiple radiation replicas exist a priori, circumventing quantum cloning constraints.This work transcends mere information conservation, advancing into the second layer of black hole physics: Characterizing the irreversible evolution of Hawking radiation in modular space.
PHYSICS OF GASES, PLASMAS, AND ELECTRIC DISCHARGES
Diagnosis of number density of negative deuterium ions excited by electron cyclotron resonance by optical emission spectroscopy
CHEN Yue, ZHU Xiaodong
2025, 74 (8): 085201. doi: 10.7498/aps.74.20241605
Abstract +
The electron cyclotron resonance (ECR) plasma is characterized by low working pressure and high dissociation rate, and it has important applications in the deuterium negative ion $ {{\mathrm{D}}}^{-} $ source technology. In this paper, the Yacora collisional-radiative model is applied to the emission spectrum diagnosis of $ {{\mathrm{D}}}^{-} $ in ECR deuterium plasma. The $ {{\mathrm{D}}}^{-} $ density is estimated by using the $ {I}_{{{\mathrm{D}}}_{{\mathrm{\alpha }}}}/{I}_{{{\mathrm{D}}}_{{\mathrm{\beta }}}} $ ratio and the relative intensity of other deuterium molecular lines, thereby avoiding complex calibration procedure of absolute intensity. The spatial structure of $ {{\mathrm{D}}}^{-} $ is studied by the multichannel emission spectrum measured in the source region and diffusion region.The experiments are conducted on a 2.45-GHz ECR plasma source at a deuterium gas pressure of 1 Pa and microwave power of 660 W. The Balmer series of atomic deuterium ($ {{\mathrm{D}}}_{{\mathrm{\alpha }}} $, $ {{\mathrm{D}}}_{{\mathrm{\beta }}} $, $ {{\mathrm{D}}}_{{\mathrm{\gamma }}} $, $ {{\mathrm{D}}}_{{\mathrm{\delta }}} $) and the Fulcher band Q-branches of molecular deuterium are measured in the source region and expanding region of the ECR plasma. It is found that the intensity of $ {{\mathrm{D}}}_{{\mathrm{\alpha }}} $ in the source region is much higher than that of $ {{\mathrm{D}}}_{{\mathrm{\beta }}} $, specifically, the $ {I}_{{{\mathrm{D}}}_{{\mathrm{\alpha }}}}/{I}_{{{\mathrm{D}}}_{{\mathrm{\beta }}}} $ ratio reaches as high as 23, indicating a selective enhancement of Balmer lines due to the mutual neutralization process of $ {{\mathrm{D}}}^{-} $. Furthermore, $ {{\mathrm{D}}}^{-} $ density in the source region is estimated to be about $ 3.6\times {10}^{15}\;{{\mathrm{m}}}^{-3} $, and the $ {{\mathrm{D}}}^{-} $ density in the expanding region decreases significantly. In the ECR plasma source region, the plasma-wall interaction is strong due to the small volume of the cavity. The recombination desorption process produces more vibrationally excited molecules, thereby further enhancing the dissociation attachment reaction, which is beneficial to the generation of deuterium negative ions. On the other hand, the axial electric field within the ECR plasma inhibits the axial transport of $ {{\mathrm{D}}}^{-} $, suggesting that the production and loss of $ {{\mathrm{D}}}^{-} $ are both localized. These characteristics of the ECR plasma source contribute to the formation of a large gradient of $ {{\mathrm{D}}}^{-} $ density between the source region and the expanding region.
Baidu
map