-
Proposal of quantum repeater architecture based on Rydberg atom quantum processors
Authors:
Yan-Lei Zhang,
Qing-Xuan Jie,
Ming Li,
Shu-Hao Wu,
Zhu-Bo Wang,
Xu-Bo Zou,
Peng-Fei Zhang,
Gang Li,
Tiancai Zhang,
Guang-Can Guo,
Chang-Ling Zou
Abstract:
Realizing large-scale quantum networks requires the generation of high-fidelity quantum entanglement states between remote quantum nodes, a key resource for quantum communication, distributed computation and sensing applications. However, entanglement distribution between quantum network nodes is hindered by optical transmission loss and local operation errors. Here, we propose a novel quantum rep…
▽ More
Realizing large-scale quantum networks requires the generation of high-fidelity quantum entanglement states between remote quantum nodes, a key resource for quantum communication, distributed computation and sensing applications. However, entanglement distribution between quantum network nodes is hindered by optical transmission loss and local operation errors. Here, we propose a novel quantum repeater architecture that synergistically integrates Rydberg atom quantum processors with optical cavities to overcome these challenges. Our scheme leverages cavity-mediated interactions for efficient remote entanglement generation, followed by Rydberg interaction-based entanglement purification and swapping. Numerical simulations, incorporating realistic experimental parameters, demonstrate the generation of Bell states with 99\% fidelity at rates of 1.1\,kHz between two nodes in local-area network (distance $0.1\,\mathrm{km}$), and can be extend to metropolitan-area ($25\,\mathrm{km}$) or intercity ($\mathrm{250\,\mathrm{km}}$, with the assitance of frequency converters) network with a rate of 0.1\,kHz. This scalable approach opens up near-term opportunities for exploring quantum network applications and investigating the advantages of distributed quantum information processing.
△ Less
Submitted 16 October, 2024;
originally announced October 2024.
-
Two Rapid Power Iterative DOA Estimators for UAV Emitter Using Massive/Ultra-massive Receive Array
Authors:
Yiwen Chen,
Feng Shu,
Qijuan Jie,
Xichao Zhan,
Xuehui Wang,
Zhongwen Sun,
Shihao Yan,
Wenlong Cai,
Peng Zhang,
Peng Chen
Abstract:
To provide rapid direction finding (DF) for unmanned aerial vehicle (UAV) emitter in future wireless networks, a low-complexity direction of arrival (DOA) estimation architecture for massive multiple input multiple output (MIMO) receiver arrays is constructed. In this paper, we propose two strategies to address the extremely high complexity caused by eigenvalue decomposition of the received signal…
▽ More
To provide rapid direction finding (DF) for unmanned aerial vehicle (UAV) emitter in future wireless networks, a low-complexity direction of arrival (DOA) estimation architecture for massive multiple input multiple output (MIMO) receiver arrays is constructed. In this paper, we propose two strategies to address the extremely high complexity caused by eigenvalue decomposition of the received signal covariance matrix. Firstly, a rapid power-iterative rotational invariance (RPI-RI) method is proposed, which adopts the signal subspace generated by power iteration to gets the final direction estimation through rotational invariance between subarrays. RPI-RI makes a significant complexity reduction at the cost of a substantial performance loss. In order to further reduce the complexity and provide a good directional measurement result, a rapid power-iterative Polynomial rooting (RPI-PR) method is proposed, which utilizes the noise subspace combined with polynomial solution method to get the optimal direction estimation. In addition, the influence of initial vector selection on convergence in the power iteration is analyzed, especially when the initial vector is orthogonal to the incident wave. Simulation results show that the two proposed methods outperform the conventional DOA estimation methods in terms of computational complexity. In particular, the RPIPR method achieves more than two orders of magnitude lower complexity than conventional methods and achieves performance close to CRLB. Moreover, it is verified that the initial vector and the relative error have a significant impact on the performance of the computational complexity.
△ Less
Submitted 23 April, 2023; v1 submitted 6 May, 2022;
originally announced May 2022.
-
Two Low-complexity DOA Estimators for Massive/Ultra-massive MIMO Receive Array
Authors:
Yiwen Chen,
Xichao Zhan,
Feng Shu,
Qijuan Jie,
Xin Cheng,
Zhihong Zhuang,
Jiangzhou Wang
Abstract:
Eigen-decomposition-based direction finding methods of using large-scale/ultra-large-scale fully-digital receive antenna arrays lead to a high or ultra-high complexity. To address the complexity dilemma, in this paper, three low-complexity estimators are proposed: partitioned subarray auto-correlation combining (PSAC), partitioned subarray cross-correlation combining (PSCC) and power iteration max…
▽ More
Eigen-decomposition-based direction finding methods of using large-scale/ultra-large-scale fully-digital receive antenna arrays lead to a high or ultra-high complexity. To address the complexity dilemma, in this paper, three low-complexity estimators are proposed: partitioned subarray auto-correlation combining (PSAC), partitioned subarray cross-correlation combining (PSCC) and power iteration max correlation successive convex approximation (PI-Max-CSCA). Compared with the conventional no-partitioned direction finding method like root multiple signal classification (Root-MUSIC), in the PSAC method, the total set of antennas are equally partitioned into subsets of antennas, called subarrays, each subarray performs independent DOA estimation, and all DOA estimates are coherently combined to give the final estimation. For a better performance, the cross-correlation among sub-arrays is further exploited in the PSCC method to achieve the near-Cramer-Rao lower bound (CRLB) performance with the help of auto-correlation. To further reduce the complexity, in the PI-Max-CSCA method, using a fraction of all subarrays to make an initial coarse direction measurement (ICDM), the power iterative method is adopted to compute the more precise steering vector (SV) by exploiting the total array, and a more accurate DOA value is found using ICDM and SV through the maximum correlation method solved by successive convex approximation. Simulation results show that as the number of antennas goes to large-scale, the proposed three methods can achieve a dramatic complexity reduction over conventional Root-MUISC. Particularly, the PSCC and PI-Max-CSCA can reach the CRLB while the PSAC shows a substantial performance loss.
△ Less
Submitted 10 August, 2022; v1 submitted 20 April, 2022;
originally announced April 2022.
-
Machine-learning-aided Massive Hybrid Analog and Digital MIMO DOA Estimation for Future Wireless Networks
Authors:
Feng Shu,
Yiwen Chen,
Xichao Zhan,
Wenlong Cai,
Mengxing Huang,
Qijuan Jie,
Yifang Li,
Baihua Shi,
Jiangzhou Wang,
Xiaohu You
Abstract:
Due to a high spatial angle resolution and low circuit cost of massive hybrid analog and digital (HAD) multiple-input multiple-output (MIMO), it is viewed as a valuable green communication technology for future wireless networks. Combining a massive HAD-MIMO with direction of arrival (DOA) will provide a high-precision even ultra-high-precision DOA measurement performance approaching the fully-dig…
▽ More
Due to a high spatial angle resolution and low circuit cost of massive hybrid analog and digital (HAD) multiple-input multiple-output (MIMO), it is viewed as a valuable green communication technology for future wireless networks. Combining a massive HAD-MIMO with direction of arrival (DOA) will provide a high-precision even ultra-high-precision DOA measurement performance approaching the fully-digital (FD) MIMO. However, phase ambiguity is a challenge issue for a massive HAD-MIMO DOA estimation. In this paper, we review three aspects: detection, estimation, and Cramer-Rao lower bound (CRLB) with low-resolution ADCs at receiver. First, a multi-layer-neural-network (MLNN) detector is proposed to infer the existence of passive emitters. Then, a two-layer HAD (TLHAD) MIMO structure is proposed to eliminate phase ambiguity using only one-snapshot. Simulation results show that the proposed MLNN detector is much better than both the existing generalized likelihood ratio test (GRLT) and the ratio of maximum eigen-value (Max-EV) to minimum eigen-value (R-MaxEV-MinEV) in terms of detection probability. Additionally, the proposed TLHAD structure can achieve the corresponding CRLB using single snapshot.
△ Less
Submitted 5 August, 2023; v1 submitted 12 January, 2022;
originally announced January 2022.
-
DOA Estimation Using Massive Receive MIMO: Basic Principle and Key Techniques
Authors:
Jiangzhou Wang,
Baihua Shi,
Feng Shu,
Qi Zhang,
Di Wu,
Qijuan Jie,
Zhihong Zhuang,
Siling Feng,
Yijin Zhang
Abstract:
As massive multiple-input multiple-output (MIMO) becomes popular, direction of arrival (DOA) measurement has been made a real renaissance due to the high-resolution achieved. Thus, there is no doubt about DOA estimation using massive MIMO. The purpose of this paper is to describe its basic principles and key techniques, to present the performance analysis, and to appreciate its engineering applica…
▽ More
As massive multiple-input multiple-output (MIMO) becomes popular, direction of arrival (DOA) measurement has been made a real renaissance due to the high-resolution achieved. Thus, there is no doubt about DOA estimation using massive MIMO. The purpose of this paper is to describe its basic principles and key techniques, to present the performance analysis, and to appreciate its engineering applications. It is anticipated that there are still many challenges in DOA estimation using massive receive MIMO, such as high circuit cost, high energy consumption and high complexity of the algorithm implementation. New researches and breakthroughs are illustrated to deal with those problems. Then, a new architecture, hybrid analog and digital (HAD) massive receive MIMO with low-resolution ADCs, is presented to strike a good balance among circuit cost, complexity and performance. Then, a novel three-dimensional (3D) angle of arrival (AOA) localization method based on geometrical center is proposed to compute the position of a passive emitter using single base station equipped with an ultra-massive MIMO system. And, it can achieve the Cramer-Rao low bound (CRLB). Here, the performance loss is also analyzed to quantify the minimum number of bits. DOA estimation will play a key role in lots of applications, such as directional modulation, beamforming tracking and alignment for 6G.
△ Less
Submitted 15 July, 2023; v1 submitted 31 August, 2021;
originally announced September 2021.
-
High-performance Passive Eigen-model-based Detectors of Single Emitter Using Massive MIMO Receivers
Authors:
Qijuan Jie,
Xichao Zhan,
Feng Shu,
Yaohui Ding,
Baihua Shi,
Yifan Li,
Jiangzhou Wang
Abstract:
For a passive direction of arrival (DoA) measurement system using massive multiple input multiple output (MIMO), it is mandatory to infer whether the emitter exists or not before performing DOA estimation operation. Inspired by the detection idea from radio detection and ranging (radar), three high-performance detectors are proposed to infer the existence of single passive emitter from the eigen-s…
▽ More
For a passive direction of arrival (DoA) measurement system using massive multiple input multiple output (MIMO), it is mandatory to infer whether the emitter exists or not before performing DOA estimation operation. Inspired by the detection idea from radio detection and ranging (radar), three high-performance detectors are proposed to infer the existence of single passive emitter from the eigen-space of sample covariance matrix of receive signal vector. The test statistic (TS) of the first method is defined as the ratio of maximum eigen-value (Max-EV) to minimum eigen-value (R-MaxEV-MinEV) while that of the second one is defined as the ratio of Max-EV to noise variance (R-MaxEV-NV). The TS of the third method is the mean of maximum eigen-value (EV) and minimum EV(M-MaxEV-MinEV). Their closed-form expressions are presented and the corresponding detection performance is given. Simulation results show that the proposed M-MaxEV-MinEV and R-MaxEV-NV methods can approximately achieve the same detection performance that is better than the traditional generalized likelihood ratio test method with false alarm probability being less than 0.3.
△ Less
Submitted 3 August, 2021;
originally announced August 2021.
-
DOA Estimation for Hybrid Massive MIMO Systems using Mixed-ADCs: Performance Loss and Energy Efficiency
Authors:
Baihua Shi,
Qi Zhang,
Rongen Dong,
Qijuan Jie,
Shihao Yan,
Feng Shu,
Jiangzhou Wang
Abstract:
Due to the power consumption and high circuit cost in antenna arrays, the practical application of massive multiple-input multiple-output (MIMO) in the sixth generation (6G) and future wireless networks is still challenging. Employing low-resolution analog-to-digital converters (ADCs) and hybrid analog and digital (HAD) structure is two low-cost choice with acceptable performance loss.In this pape…
▽ More
Due to the power consumption and high circuit cost in antenna arrays, the practical application of massive multiple-input multiple-output (MIMO) in the sixth generation (6G) and future wireless networks is still challenging. Employing low-resolution analog-to-digital converters (ADCs) and hybrid analog and digital (HAD) structure is two low-cost choice with acceptable performance loss.In this paper, the combination of the mixed-ADC architecture and HAD structure employed at receiver is proposed for direction of arrival (DOA) estimation, which will be applied to the beamforming tracking and alignment in 6G. By adopting the additive quantization noise model, the exact closed-form expression of the Cramér-Rao lower bound (CRLB) for the HAD architecture with mixed-ADCs is derived. Moreover, the closed-form expression of the performance loss factor is derived as a benchmark. In addition, to take power consumption into account, energy efficiency is also investigated in our paper. The numerical results reveal that the HAD structure with mixed-ADCs can significantly reduce the power consumption and hardware cost. Furthermore, that architecture is able to achieve a better trade-off between the performance loss and the power consumption. Finally, adopting 2-4 bits of resolution may be a good choice in practical massive MIMO systems.
△ Less
Submitted 19 May, 2023; v1 submitted 21 July, 2021;
originally announced July 2021.
-
Higgsino Asymmetry and Direct-Detection Constraints of Light Dark Matter in the NMSSM with Non-Universal Higgs Masses
Authors:
Kun Wang,
Jingya Zhu,
Quanlin Jie
Abstract:
In this work, we study the direct-detection constraints of light dark matter in the next-to minimal supersymmetric standard model (NMSSM) with non-universal Higgs masses (NUHM), especially the correlation between higgsino asymmetry and spin-dependent cross section. Finally, we get the following conclusions: (i) The spin-dependent cross section is proportional to the square of higgsino asymmetry in…
▽ More
In this work, we study the direct-detection constraints of light dark matter in the next-to minimal supersymmetric standard model (NMSSM) with non-universal Higgs masses (NUHM), especially the correlation between higgsino asymmetry and spin-dependent cross section. Finally, we get the following conclusions: (i) The spin-dependent cross section is proportional to the square of higgsino asymmetry in dark matter $\tildeχ^0_1$ in the NMSSM-NUHM. (ii) For highly singlino-dominated dark matter $\tildeχ^0_1$, the relic density can be sufficient, but the higgsino asymmetry and spin-dependent cross section is always small. (iii) With a sizeable higgsino component in the light dark matter, the higgsino asymmetry and spin-dependent cross section can be sizeable, but dark matter relic density is always small, thus it can escape the direct detections. (iv) Light dark matter in the $h_2$- and $Z$-funnel annihilation channels with sufficient relic density can be covered by future LUX-ZEPLIN (LZ) 7-ton in spin-dependent detections.
△ Less
Submitted 25 November, 2020;
originally announced November 2020.
-
High Entropy Alloys Mined From Phase Diagrams
Authors:
Qi Jie,
Andrew Cheung,
S. Joseph Poon
Abstract:
High entropy alloys (HEA) show promise as a new type of high-performance structural material. Their vast degrees of freedom provide for extensive opportunities to design alloys with tailored properties. However, the compositional complexities of HEAs present great challenges for alloy design. Current approaches have shown limited reliability in accounting for the compositional regions of single so…
▽ More
High entropy alloys (HEA) show promise as a new type of high-performance structural material. Their vast degrees of freedom provide for extensive opportunities to design alloys with tailored properties. However, the compositional complexities of HEAs present great challenges for alloy design. Current approaches have shown limited reliability in accounting for the compositional regions of single solid solution and composite phases. We present a phenomenological method, analyzing binary phase diagrams to predict HEA phase formation on the hypothesis that the structural stability of HEAs is encoded within. Accordingly, we introduce a small number of phase-diagram inspired parameters and employ machine learning to partition the formation region of 500+ reported HEA compositions. The model achieved a single phase HEA prediction rate >80 %. To validate our method, we demonstrated the capability of this method in predicting HEA solid solution phases with and without intermetallics in 30 randomly selected complex compositions, with a success rate of 77 %. The presented efficient search approach with high predictive capability can be exploited to complement computation-intense methods in providing a powerful platform for the design of high entropy alloys.
△ Less
Submitted 18 April, 2019;
originally announced April 2019.
-
Initial state radiation correction and its effect to data-taking scheme for $σ^{\mathrm{B}}(e^+e^-\to ZH)$ measurement
Authors:
Shuang Han,
Gang Li,
Xiang Zhou,
Quan-Lin Jie
Abstract:
The measurement of Born cross section of $e^+e^-\to ZH$ process is one of the major goals of the future Circular Electron Positron Collider, which may reach a precision of 0.5%. Such unprecedented precision must be guaranteed by both theoretical and experimental sides, such as the calculations of high order corrections, the knowledge of the $σ^{\mathrm{B}}(e^+e^-\to ZH)$ line shape. The uncertaint…
▽ More
The measurement of Born cross section of $e^+e^-\to ZH$ process is one of the major goals of the future Circular Electron Positron Collider, which may reach a precision of 0.5%. Such unprecedented precision must be guaranteed by both theoretical and experimental sides, such as the calculations of high order corrections, the knowledge of the $σ^{\mathrm{B}}(e^+e^-\to ZH)$ line shape. The uncertainty of radiative correction factor due to lack of knowledge of the $σ^{\mathrm{B}}(e^+e^-\to ZH)$ line shape is evaluated in this work. Therefore, the dedicated data-taking schemes are proposed in order to precisely calculate the ISR correction factor.
△ Less
Submitted 23 April, 2019; v1 submitted 25 February, 2019;
originally announced February 2019.
-
The semi-constrained NMSSM in light of muon g-2, LHC, and dark matter constraints
Authors:
Kun Wang,
Fei Wang,
Jingya Zhu,
Quanlin Jie
Abstract:
The semi-constrained NMSSM (scNMSSM) extends the MSSM by a singlet field, and requires unification of the soft SUSY breaking terms in the squark and slepton sectors, while it allows that in the Higgs sector to be different. We try to interpret the muon g-2 in the scNMSSM, under the constraints of 125 GeV Higgs data, B physics, searches for low and high mass resonances, searches for SUSY particles…
▽ More
The semi-constrained NMSSM (scNMSSM) extends the MSSM by a singlet field, and requires unification of the soft SUSY breaking terms in the squark and slepton sectors, while it allows that in the Higgs sector to be different. We try to interpret the muon g-2 in the scNMSSM, under the constraints of 125 GeV Higgs data, B physics, searches for low and high mass resonances, searches for SUSY particles at the LHC, dark matter relic density by WMAP/Planck, and direct searches for dark matter by LUX, XENON1T, and PandaX-II. We find that under the above constraints, the scNMSSM can still (i) satisfy muon g-2 at 1$σ$ level, with a light muon sneutrino and light chargino; (ii) predict a highly-singlet-dominated 95~GeV Higgs, with a diphoton rate as hinted at by CMS data, because of a light higgsino-like chargino and moderate $λ$; (iii) get low fine tuning from the GUT scale with small $μ_{\rm eff},\, M_0,\, M_{1/2},\, {\rm and}\, A_0$, with a lighter stop mass which can be as low as about 500 GeV, which can be further checked in future studies with search results from the 13~TeV LHC; (iv) have the lightest neutralino be singlino-dominated or higgsino-dominated, while the bino and wino are heavier because of high gluino bounds at the LHC and universal gaugino conditions at the GUT scale; (v) satisfy all the above constraints, although it is not easy for the lightest neutralino, as the only dark matter candidate, to get enough relic density. Several ways to increase relic density are discussed.
△ Less
Submitted 24 November, 2018; v1 submitted 11 November, 2018;
originally announced November 2018.
-
Search for $h_{c}\rightarrowπ^+π^-J/ψ$ via $ψ(3686)\rightarrow π^0π^+π^-J/ψ$
Authors:
BESIII Collaboration,
M. Ablikim,
M. N. Achasov,
S. Ahmed,
X. C. Ai,
M. Albrecht,
D. J. Ambrose,
A. Amoroso,
F. F. An,
Q. An,
J. Z. Bai,
O. Bakina,
R. Baldini Ferroli,
Y. Ban,
D. W. Bennett,
J. V. Bennett,
N. Berger,
M. Bertani,
D. Bettoni,
J. M. Bian,
F. Bianchi,
E. Boger,
I. Boyko,
R. A. Briere,
H. Cai
, et al. (393 additional authors not shown)
Abstract:
Using a data sample of $448.1\times10^6$ $ψ(3686)$ events collected with the BESIII detector operating at the BEPCII, we perform search for the hadronic transition $h_c\rightarrowπ^+π^-J/ψ$ via $ψ(3686)\rightarrowπ^0h_c$. No signals of the transition are observed, and the upper limit on the product branching fraction $\mathcal{B}(ψ(3686)\rightarrowπ^0h_c)\mathcal{B}(h_c\rightarrowπ^+π^-J/ψ)$ at th…
▽ More
Using a data sample of $448.1\times10^6$ $ψ(3686)$ events collected with the BESIII detector operating at the BEPCII, we perform search for the hadronic transition $h_c\rightarrowπ^+π^-J/ψ$ via $ψ(3686)\rightarrowπ^0h_c$. No signals of the transition are observed, and the upper limit on the product branching fraction $\mathcal{B}(ψ(3686)\rightarrowπ^0h_c)\mathcal{B}(h_c\rightarrowπ^+π^-J/ψ)$ at the 90% confidence level is determined to be $2.0\times10^{-6}$. This is the most stringent upper limit to date.
△ Less
Submitted 16 March, 2018; v1 submitted 10 January, 2018;
originally announced January 2018.
-
Preliminary study of light yield dependence on LAB liquid scintillator composition
Authors:
Xing-Chen Ye,
Bo-Xiang Yu,
Xiang Zhou,
Li Zhao,
Ya-Yun Ding,
Quan-Lin Jie,
Shun-Li Niu,
Meng-Chao Liu,
Xue-Feng Ding,
Xuan Zhang,
Li Zhou,
Jian Fang,
Hai-Tao Chen,
Wei Hu,
Jia-Qing Yan,
Hang Zhao,
Dao-Jin Hong
Abstract:
Liquid scintillator (LS) will be adopted as the detector material in JUNO (Jiangmen Underground Neutrino Observatory). The energy resolution requirement of JUNO is 3%, which has never previously been reached. To achieve this energy resolution, the light yield of liquid scintillator is an important factor. PPO (the fluor) and bis-MSB (the wavelength shifter) are the two main materials dissolved in…
▽ More
Liquid scintillator (LS) will be adopted as the detector material in JUNO (Jiangmen Underground Neutrino Observatory). The energy resolution requirement of JUNO is 3%, which has never previously been reached. To achieve this energy resolution, the light yield of liquid scintillator is an important factor. PPO (the fluor) and bis-MSB (the wavelength shifter) are the two main materials dissolved in LAB. To study the influence of these two materials on the transmission of scintillation photons in LS, 25 and 12 cm-long quartz vessels were used in a light yield experiment. LS samples with different concentration of PPO and bis-MSB were tested. At these lengths, the light yield growth is not obvious when the concentration of PPO is higher than 4 g/L. The influence from bis-MSB becomes insignificant when its concentration is higher than 8 mg/L. This result could provide some useful suggestions for the JUNO LS.
△ Less
Submitted 31 May, 2015;
originally announced June 2015.
-
Cluster Density Matrix Embedding Theory for Quantum Spin Systems
Authors:
Zhuo Fan,
Quanlin Jie
Abstract:
We applied cluster density matrix embedding theory, with some modifications, to a spin lattice system. The reduced density matrix of the impurity cluster is embedded in the bath states, which are a set of block-product states. The ground state of the impurity model is formulated using a variational wave function. We tested this theory in a two-dimensional (2-D) spin-1/2 J1-J2 model for a square la…
▽ More
We applied cluster density matrix embedding theory, with some modifications, to a spin lattice system. The reduced density matrix of the impurity cluster is embedded in the bath states, which are a set of block-product states. The ground state of the impurity model is formulated using a variational wave function. We tested this theory in a two-dimensional (2-D) spin-1/2 J1-J2 model for a square lattice. The ground-state energy (GSE) and the location of the phase boundaries agree well with the most accurate previous results obtained using the quantum Monte Carlo and coupled cluster methods. Moreover, this cluster density matrix embedding theory is cost-effective and convenient for calculating the von Neumann entropy, which is related to the quantum phase transition.
△ Less
Submitted 21 April, 2015;
originally announced April 2015.
-
Parameterized excitation operators for coupled cluster method
Authors:
Quanlin Jie
Abstract:
We present a coupled cluster method (CCM) with optimized excitation operators. The efficiency comes from a parameterized form of excitation operators. The parameters are found by variational optimization procedure. The resultant number of excitation operators is much smaller than that of the conventional CCM theory. This property makes it possible to apply the method in systems of solid state phys…
▽ More
We present a coupled cluster method (CCM) with optimized excitation operators. The efficiency comes from a parameterized form of excitation operators. The parameters are found by variational optimization procedure. The resultant number of excitation operators is much smaller than that of the conventional CCM theory. This property makes it possible to apply the method in systems of solid state physics. Starting from Hartree-Fock state as the reference state, i.e., the Fermi sea, we search for particle-hole excitation operators such that the wave function of configuration interaction in terms of these excitation operators spans a good approximation to the ground state. The Match-pursuit algorithm is capable of doing the search of the excitation operators. The resultant operators are our excitation operators for the CCM wave function. We test the method by two dimensional fermionic Hubbard model on a square lattice.
△ Less
Submitted 8 December, 2014;
originally announced December 2014.
-
Problem Complexity Research from Energy Perspective
Authors:
Feng Pan,
Heng-liang Zhang,
Qi Jie
Abstract:
Computational complexity is a particularly important objective. The idea of Landauer principle was extended through mapping three classic problems (sorting,ordered searching and max of N unordered numbers) into Maxwell demon thought experiment in this paper. The problems'complexity is defined on the entropy basis and the minimum energy required to solve them are rigorous deduced from the perspecti…
▽ More
Computational complexity is a particularly important objective. The idea of Landauer principle was extended through mapping three classic problems (sorting,ordered searching and max of N unordered numbers) into Maxwell demon thought experiment in this paper. The problems'complexity is defined on the entropy basis and the minimum energy required to solve them are rigorous deduced from the perspective of energy (entropy) and the second law of thermodynamics. Then the theoretical energy consumed by real program and basic operators of classical computer are both analyzed, the time complexity lower bounds of three problems'all possible algorithms are derived in this way. The lower bound is also deduced for the two n*n matrix multiplication problem. In the end, the reason why reversible computation is impossible and the possibility of super-linear energy consumption capacity which may be the power behind quantum computation are discussed, a conjecture is proposed which may prove NP!=P. The study will bring fresh and profound understanding of computation complexity.
△ Less
Submitted 16 September, 2013;
originally announced September 2013.
-
Electronic thermoelectric power factor and metal-insulator transition in FeSb2
Authors:
Qing Jie,
Rongwei Hu,
Emil Bozin,
A. Llobet,
I. Zaliznyak,
C. Petrovic,
Q. Li
Abstract:
We show that synthesis-induced Metal -Insulator transition (MIT) for electronic transport along the orthorombic c axis of FeSb$_{2}$ single crystals has greatly enhanced electrical conductivity while keeping the thermopower at a relatively high level. By this means, the thermoelectric power factor is enhanced to a new record high S$^{2}$$σ$ $\sim$ 8000 $μ$WK$^{-2}$cm$^{-1}$ at 28 K. We find that t…
▽ More
We show that synthesis-induced Metal -Insulator transition (MIT) for electronic transport along the orthorombic c axis of FeSb$_{2}$ single crystals has greatly enhanced electrical conductivity while keeping the thermopower at a relatively high level. By this means, the thermoelectric power factor is enhanced to a new record high S$^{2}$$σ$ $\sim$ 8000 $μ$WK$^{-2}$cm$^{-1}$ at 28 K. We find that the large thermopower in FeSb$_{2}$ can be rationalized within the correlated electron model with two bands having large quasiparaticle disparity, whereas MIT is induced by subtle structural differences. The results in this work testify that correlated electrons can produce extreme power factor values.
△ Less
Submitted 11 October, 2012;
originally announced October 2012.
-
Thermal Stability of Thermoelectric Materials via In Situ Resistivity Measurements
Authors:
K. C. Lukas,
W. S. Liu,
Q. Jie,
Z. F. Ren,
C. P. Opeil
Abstract:
An experimental setup for determining the electrical resistivity of several types of thermoelectric materials over the temperature range 20 < T < 550 C is described in detail. One resistivity measurement during temperature cycling is also explained for Cu0.01Bi2Te2.7Se0.3 while a second measurement is made on Yb0.35Co4Sb12 as a function of time at 400 C. Both measurements confirm that the material…
▽ More
An experimental setup for determining the electrical resistivity of several types of thermoelectric materials over the temperature range 20 < T < 550 C is described in detail. One resistivity measurement during temperature cycling is also explained for Cu0.01Bi2Te2.7Se0.3 while a second measurement is made on Yb0.35Co4Sb12 as a function of time at 400 C. Both measurements confirm that the materials are thermally stable for the temperature range and time period measured. Measurements made during temperature cycling show an irreversible decrease in the electrical resistivity of Cu0.01Bi2Te2.7Se0.3 when the measuring temperature exceeds the pressing temperature. Several other possible uses of such a system include but are not limited to studying the effects of annealing and/or oxidation as a function of both temperature and time.
△ Less
Submitted 11 June, 2012;
originally announced June 2012.
-
Transport Properties of Stripe-Ordered High Tc Cuprates
Authors:
Qing Jie,
Su Jung Han,
Ivo Dimitrov,
J. M. Tranquada,
Qiang Li
Abstract:
Transport measurements provide important characterizations of the nature of stripe order in the cuprates. Initial studies of systems such as La(1.6-x)Nd(0.4)Sr(x)CuO(4) demonstrated the strong anisotropy between in-plane and c-axis resistivities, but also suggested that stripe order results in a tendency towards insulating behavior within the planes at low temperature. More recent work on La(2-x)B…
▽ More
Transport measurements provide important characterizations of the nature of stripe order in the cuprates. Initial studies of systems such as La(1.6-x)Nd(0.4)Sr(x)CuO(4) demonstrated the strong anisotropy between in-plane and c-axis resistivities, but also suggested that stripe order results in a tendency towards insulating behavior within the planes at low temperature. More recent work on La(2-x)Ba(x)CuO(4) with x=1/8 has revealed the occurrence of quasi-two-dimensional superconductivity that onsets with spin-stripe order. The suppression of three-dimensional superconductivity indicates a frustration of the interlayer Josephson coupling, motivating a proposal that superconductivity and stripe order are intertwined in a pair-density-wave state. Complementary characterizations of the low-energy states near the Fermi level are provided by measurements of the Hall and Nernst effects, each revealing intriguing signatures of stripe correlations and ordering. We review and discuss this work.
△ Less
Submitted 24 May, 2012;
originally announced May 2012.
-
Probing the connections between superconductivity, stripe order, and structure in La1.905Ba0.095Cu1-yZnyO4
Authors:
Jinsheng Wen,
Zhijun Xu,
Guangyong Xu,
Qing Jie,
M. Huecker,
A. Zheludev,
Wei Tian,
B. L. Winn,
J. L. Zarestky,
D. K. Singh,
Tao Hong,
Qiang Li,
Genda Gu,
J. M. Tranquada
Abstract:
The superconducting system La2-xBaxCuO4 is known to show a minimum in the transition temperature, Tc, at x = 1/8 where maximal stripe order is pinned by the anisotropy within the CuO2 planes that occurs in the low-temperature-tetragonal (LTT) crystal structure. For x = 0.095, where Tc reaches its maximum value of 32 K, there is a roughly coincident structural transition to a phase that is very clo…
▽ More
The superconducting system La2-xBaxCuO4 is known to show a minimum in the transition temperature, Tc, at x = 1/8 where maximal stripe order is pinned by the anisotropy within the CuO2 planes that occurs in the low-temperature-tetragonal (LTT) crystal structure. For x = 0.095, where Tc reaches its maximum value of 32 K, there is a roughly coincident structural transition to a phase that is very close to LTT. Here we present a neutron scattering study of the structural transition, and demonstrate how features of it correlate with anomalies in the magnetic susceptibility, electrical resistivity, thermal conductivity, and thermoelectric power. We also present measurements on a crystal with 1% Zn substituted for Cu, which reduces Tc to 17 K, enhances the spin stripe order, but has much less effect on the structural transition. We make the case that the structural transition correlates with a reduction of the Josephson coupling between the CuO2 layers, which interrupts the growth of the superconducting order. We also discuss evidence for two-dimensional superconducting fluctuations in the normal state, analyze the effective magnetic moment per Zn impurity, and consider the significance of the anomalous thermopower often reported in the stripe-ordered phase.
△ Less
Submitted 12 April, 2012; v1 submitted 22 November, 2011;
originally announced November 2011.
-
Iron-chalcogenide FeSe$_{0.5}$Te$_{0.5}$ coated superconducting tapes for high field applications
Authors:
Weidong Si,
Juan Zhou,
Qing Jie,
Ivo Dimitrov,
V. Solovyov,
P. D. Johnson,
J. Jaroszynski,
V. Matias,
C. Sheehan,
Qiang Li
Abstract:
The high upper critical field characteristic of the recently discovered iron-based superconducting chalcogenides opens the possibility of developing a new type of non-oxide high-field superconducting wires. In this work, we utilize a buffered metal template on which we grow a textured FeSe$_{0.5}$Te$_{0.5}$ layer, an approach developed originally for high temperature superconducting coated conduct…
▽ More
The high upper critical field characteristic of the recently discovered iron-based superconducting chalcogenides opens the possibility of developing a new type of non-oxide high-field superconducting wires. In this work, we utilize a buffered metal template on which we grow a textured FeSe$_{0.5}$Te$_{0.5}$ layer, an approach developed originally for high temperature superconducting coated conductors. These tapes carry high critical current densities (>1$\times10^{4}$A/cm$^{2}$) at about 4.2K under magnetic field as high as 25 T, which are nearly isotropic to the field direction. This demonstrates a very promising future for iron chalcogenides for high field applications at liquid helium temperatures. Flux pinning force analysis indicates a point defect pinning mechanism, creating prospects for a straightforward approach to conductor optimization.
△ Less
Submitted 13 June, 2011;
originally announced June 2011.
-
Nanostructures and Defects in Non-equilibrium Synthesized Filled Skutterudite CeFe4Sb12
Authors:
Juan Zhou,
Qing Jie,
Lijun Wu,
Ivo Dimitrov,
Xun Shi,
Qiang Li
Abstract:
We studied nanoprecipitates and defects in p-type filled skutterudite CeFe4Sb12 prepared by non-equilibrium melt-spinning plus spark plasma sintering method using transmission electron microscopy. Nanoprecipitates with mostly spherical shapes and different sizes (from several nm to several tens of nm) have been observed. The most typically observed nanoprecipitates are shown to be Sb-rich. Superla…
▽ More
We studied nanoprecipitates and defects in p-type filled skutterudite CeFe4Sb12 prepared by non-equilibrium melt-spinning plus spark plasma sintering method using transmission electron microscopy. Nanoprecipitates with mostly spherical shapes and different sizes (from several nm to several tens of nm) have been observed. The most typically observed nanoprecipitates are shown to be Sb-rich. Superlattices with a periodicity of about 3.576 nm were induced by the ordering of excessive Sb atoms along the c direction. These nanoprecipitates usually share coherent interfaces with the surrounding matrix and induce anisotropic and strong strain fields in the surrounding matrix. Nanoprecipitates with compositions close to CeSb2 are much larger in size (~ 30 nm) and have orthorhombic structures. Various defects were typically observed on the interfaces between these nanoprecipitates and the matrix. The strain fields induced by these nanoprecipitates are less distinct, possibly because part of the strains has been released by the formation of defects.
△ Less
Submitted 11 December, 2010;
originally announced December 2010.
-
Uniaxial linear resistivity of superconducting La(1.905)Ba(0.095)CuO(4) induced by an external magnetic field
Authors:
Jinsheng Wen,
Qing Jie,
Qiang Li,
M. Huecker,
M. v. Zimmermann,
Su Jung Han,
Zhijun Xu,
D. K. Singh,
R. M. Konik,
Liyuan Zhang,
Genda Gu,
J. M. Tranquada
Abstract:
We present an experimental study of the anisotropic resistivity of superconducting La(2-x)Ba(x)CuO(4) with x=0.095 and transition temperature Tc=32 K. In a magnetic field perpendicular to the CuO(2) layers, H(perp), we observe that the resistivity perpendicular to the layers, ρ(perp), becomes finite at a temperature consistent with previous studies on very similar materials; however, the onset of…
▽ More
We present an experimental study of the anisotropic resistivity of superconducting La(2-x)Ba(x)CuO(4) with x=0.095 and transition temperature Tc=32 K. In a magnetic field perpendicular to the CuO(2) layers, H(perp), we observe that the resistivity perpendicular to the layers, ρ(perp), becomes finite at a temperature consistent with previous studies on very similar materials; however, the onset of finite parallel resistivity, ρ(par), occurs at a much higher temperature. This behavior contradicts conventional theory, which predicts that ρ(perp) and ρ(par) should become finite at the same temperature. Voltage vs. current measurements near the threshold of voltage detectability indicate linear behavior perpendicular to the layers, becoming nonlinear at higher currents, while the behavior is nonlinear from the onset parallel to the layers. These results, in the presence of moderate H(perp), appear consistent with superconducting order parallel to the layers with voltage fluctuations between the layers due to thermal noise. In search of uncommon effects that might help to explain this behavior, we have performed diffraction measurements that provide evidence for H(perp)-induced charge and spin stripe order. The field-induced decoupling of superconducting layers is similar to the decoupled phase observed previously in La(2-x)Ba(x)CuO(4) with x=1/8 in zero field.
△ Less
Submitted 28 March, 2012; v1 submitted 31 August, 2010;
originally announced September 2010.
-
Strong impact of grain boundaries on the thermoelectric properties of non-equilibrium synthesized p-type Ce1.05Fe4Sb12.04 filled skutterudites with nanostructure
Authors:
Qing Jie,
Juan Zhou,
Xun Shi,
Ivo K. Dimitrov,
Qiang Li
Abstract:
p-type Ce1.05Fe4Sb12.04 filled skutterudites with much improved thermoelectric properties have been synthesized by rapidly converting nearly amorphous ribbons into crystalline pellets under pressure. It is found that this process greatly suppresses grain growth and second phase formation/segregation, and hence results in the samples consisting of nano-sized grains with strongly-coupled grain bound…
▽ More
p-type Ce1.05Fe4Sb12.04 filled skutterudites with much improved thermoelectric properties have been synthesized by rapidly converting nearly amorphous ribbons into crystalline pellets under pressure. It is found that this process greatly suppresses grain growth and second phase formation/segregation, and hence results in the samples consisting of nano-sized grains with strongly-coupled grain boundaries, as observed by transmission electron microscopy. The room temperature carrier mobility in these samples is significantly higher (nearly double) than those in the samples of the same starting composition made by the conventional solid-state reaction. Nanostructure reduces the lattice thermal conductivity, while cleaner grain boundaries permit higher electron conduction.
△ Less
Submitted 29 June, 2010;
originally announced June 2010.
-
Disappearance of static magnetic order and evolution of spin fluctuations in Fe$_{1+δ}$Se$_{x}$Te$_{1-x}$
Authors:
Zhijun Xu,
Jinsheng Wen,
Guangyong Xu,
Qing Jie,
Zhiwei Lin,
Qiang Li,
Songxue Chi,
D. K. Singh,
Genda Gu,
J. M. Tranquada
Abstract:
We report neutron scattering studies on static magnetic orders and spin excitations in the Fe-based chalcogenide system Fe$_{1+δ}$Se$_{x}$Te$_{1-x}$ with different Fe and Se compositions. Short-range static magnetic order with the "bicollinear" spin configuration is found in all non-superconducting samples, with strong low-energy magnetic excitations near the $(0.5,0)$ in-plane wave-vector (using…
▽ More
We report neutron scattering studies on static magnetic orders and spin excitations in the Fe-based chalcogenide system Fe$_{1+δ}$Se$_{x}$Te$_{1-x}$ with different Fe and Se compositions. Short-range static magnetic order with the "bicollinear" spin configuration is found in all non-superconducting samples, with strong low-energy magnetic excitations near the $(0.5,0)$ in-plane wave-vector (using the two-Fe unit cell) for Se doping up to 45%. When the static order disappears and bulk superconductivity emerges, the spectral weight of the magnetic excitations shifts to the region of reciprocal space near the in-plane wave-vector $(0.5,0.5)$, corresponding to the "collinear" spin configuration. Our results suggest that spin fluctuations associated with the collinear magnetic structure appear to be universal in all Fe-based superconductors, and there is a strong correlation between superconductivity and the character of the magnetic order/fluctuations in
△ Less
Submitted 26 May, 2010;
originally announced May 2010.
-
Einstein Modes in the Phonon Density of States of the Single-Filled Skutterudite Yb$_{0.2}$Co$_{4}$Sb$_{12}$
Authors:
I. K. Dimitrov,
M. E. Manley,
S. M. Shapiro,
J. Yang,
W. Zhang,
L. D. Chen,
Q. Jie,
G. Ehlers,
A. Podlesnyak,
J. Camacho,
Qiang Li
Abstract:
Measurements of the phonon density of states by inelastic neutron \emph{time-of-flight} scattering and specific heat measurements along with first principles calculations, provide compelling evidence for the existence of an Einstein oscillator ("\emph{rattler}") at $ω_{E1} \approx$ 5.0 meV in the filled skutterudite Yb$_{0.2}$Co$_{4}$Sb$_{12}$. Multiple dispersionless modes in the measured density…
▽ More
Measurements of the phonon density of states by inelastic neutron \emph{time-of-flight} scattering and specific heat measurements along with first principles calculations, provide compelling evidence for the existence of an Einstein oscillator ("\emph{rattler}") at $ω_{E1} \approx$ 5.0 meV in the filled skutterudite Yb$_{0.2}$Co$_{4}$Sb$_{12}$. Multiple dispersionless modes in the measured density of states of Yb$_{0.2}$Co$_{4}$Sb$_{12}$ at intermediate transfer energies (14 meV $\leq$ \emph{$ω$} $\leq$ 20 meV) are exhibited in both the experimental and theoretical \emph{density-of-states} of the Yb-filled specimen. A peak at 12.4 meV is shown to coincide with a second Einstein mode at \emph{$ω_{E2} \approx$} 12.8 meV obtained from heat capacity data. The local modes at intermediate transfer energies are attributed to altered properties of the host CoSb$_{3}$ cage as a result of Yb-filling. It is suggested that these modes are owed to a complementary mechanism for the scattering of heat-carrying phonons in addition to the mode observed at $ω_{E1} \, \approx$ 5.0 meV. Our observations offer a plausible explanation for the significantly-higher \textit{dimensionless figures of merit} of filled skutterudites, compared to their parent compounds.
△ Less
Submitted 4 October, 2010; v1 submitted 20 May, 2010;
originally announced May 2010.
-
Local Structure of Thermoelectric Ca3Co4O9
Authors:
T. A. Tyson,
Z. Chen,
Q. Jie,
Q. Li,
J. J. Tu
Abstract:
We have combined temperature dependent local structural measurements with first principles density functional calculations to develop a three dimensional local structure model of the misfit system [Ca2CoO3][CoO2]1.61 (referred to as Ca3Co4O9) which has a rock salt structure stacked incommensurately on a hexagonal CoO2 lattice. The local structural measurements reveal a low coordination of Co(2)-…
▽ More
We have combined temperature dependent local structural measurements with first principles density functional calculations to develop a three dimensional local structure model of the misfit system [Ca2CoO3][CoO2]1.61 (referred to as Ca3Co4O9) which has a rock salt structure stacked incommensurately on a hexagonal CoO2 lattice. The local structural measurements reveal a low coordination of Co(2)-O bonds in the rock salt layer with large static structural disorder. The temperature dependence of the Co(1)-Co(1) bond correlations in the CoO2 layer are found to be normal above ~75K and with a very small static disorder component. An anomalous enhancement in the Co(1)-Co(1) correlations occurs at the onset of long-range magnetic order. Density functional computations suggest that the reduction of the coordination of Co(2) is due to the formation of chains of Co(2)Ox in the a-b plane linked to the Ca-O layers by c-axis Co(2)-O bonds. The reduced dimensionality introduced by the chain-like structure in the rock salt layer and high atomic order in the C
△ Less
Submitted 4 November, 2008;
originally announced November 2008.
-
Evidence for unusual superconducting correlations coexisting with stripe order in La(1.875)Ba(0.125)CuO(4)
Authors:
J. M. Tranquada,
G. D. Gu,
M. Huecker,
Q. Jie,
H. -J. Kang,
R. Klingeler,
Q. Li,
N. Tristan,
J. S. Wen,
G. Y. Xu,
Z. J. Xu,
J. Zhou,
M. v. Zimmermann
Abstract:
We present new x-ray and neutron scattering measurements of stripe order in La(1.875)Ba(0.125)CuO(4), along with low-field susceptibility, thermal conductivity, and specific heat data. We compare these with previously reported results for resistivity and thermopower. Temperature-dependent features indicating transitions (or crossovers) are correlated among the various experimental quantities. Ta…
▽ More
We present new x-ray and neutron scattering measurements of stripe order in La(1.875)Ba(0.125)CuO(4), along with low-field susceptibility, thermal conductivity, and specific heat data. We compare these with previously reported results for resistivity and thermopower. Temperature-dependent features indicating transitions (or crossovers) are correlated among the various experimental quantities. Taking into account recent spectroscopic studies, we argue that the most likely interpretation of the complete collection of results is that an unusual form of two-dimensional superconducting correlations appears together with the onset of spin-stripe order. Recent theoretical proposals for a sinusoidally-modulated superconducting state compatible with stripe order provide an intriguing explanation of our results and motivate further experimental tests. We also discuss evidence for one-dimensional pairing correlations that appear together with the charge order. With regard to the overall phenomenology, we consider the degree to which similar behavior may have been observed in other cuprates, and describe possible connections to various puzzling phenomena in cuprate superconductors.
△ Less
Submitted 3 November, 2008; v1 submitted 3 September, 2008;
originally announced September 2008.
-
Approximating the ground state of fermion system by multiple determinant states: matching pursuit approach
Authors:
Quanlin Jie
Abstract:
We present a simple and stable numerical method to approximate the ground state of a quantum many-body system by multiple determinant states. This method searches these determinant states one by one according to the matching pursuit algorithm. The first determinant state is identical to that of the Hartree-Fock theory. Calculations for two-dimensional Hubbard model serve as a demonstration.
We present a simple and stable numerical method to approximate the ground state of a quantum many-body system by multiple determinant states. This method searches these determinant states one by one according to the matching pursuit algorithm. The first determinant state is identical to that of the Hartree-Fock theory. Calculations for two-dimensional Hubbard model serve as a demonstration.
△ Less
Submitted 13 February, 2008;
originally announced February 2008.
-
Quantum-classical correspondence in entanglement production: Entropy and classical tori
Authors:
Shi-hui Zhang,
Quan-lin Jie
Abstract:
We analyze the connections between entanglement dynamics and classical trajectories in a semiclassi-cal regime for two systems: A pair of coupled oscillators and the Jaynes-Cummings model. We find that entanglement production depends on classical invariant tori and such phenomenon is closely related to the power spectra of classical trajectories. Classical power spectrum with a larger number of…
▽ More
We analyze the connections between entanglement dynamics and classical trajectories in a semiclassi-cal regime for two systems: A pair of coupled oscillators and the Jaynes-Cummings model. We find that entanglement production depends on classical invariant tori and such phenomenon is closely related to the power spectra of classical trajectories. Classical power spectrum with a larger number of frequency com-ponents corresponds to larger entanglement. We introduce a frequency entropy to describe the classical frequency distribution, and find that there is good correspondence between the classical frequency entro-pies and the maximum von Neumann entropies.
△ Less
Submitted 19 January, 2008;
originally announced January 2008.
-
Decoherence and localization in tunneling process under influence of one external degree of freedom
Authors:
Quanlin Jie,
Bambi Hu,
Guangjong Dong
Abstract:
We investigate numerically the tunneling effect under influence of another particle in a double well system. Such influence from only one degree of freedom makes decoherence and quantum-classical transition, i.e., suppression of the tunneling effect. The decoherence happens even for cases that the influence is from a particle of very small mass, and it has virtually no effect in the correspondin…
▽ More
We investigate numerically the tunneling effect under influence of another particle in a double well system. Such influence from only one degree of freedom makes decoherence and quantum-classical transition, i.e., suppression of the tunneling effect. The decoherence happens even for cases that the influence is from a particle of very small mass, and it has virtually no effect in the corresponding classical dynamics. There are cases similar to dynamical localization that the suppressed tunneling rate is several times smaller than the classical counterpart. This result is relevant for understanding quantitatively the dynamical process of decoherence and quantum to classical transition.
△ Less
Submitted 4 January, 2006;
originally announced January 2006.
-
An alternative implementation of the Lanczos algorithm for wave function propagation
Authors:
Quanlin Jie,
Dunhuan Liu
Abstract:
We reformulate the Lanczos algorithm for quantum wave function propagation in terms of variational principle. By including some basis states of previous time steps into the variational subspace, the resultant accuracy increases by several orders. Numerical errors of the alternative method accumulate much slower than that of the original Lanczos method. There is almost no extra numeric cost for t…
▽ More
We reformulate the Lanczos algorithm for quantum wave function propagation in terms of variational principle. By including some basis states of previous time steps into the variational subspace, the resultant accuracy increases by several orders. Numerical errors of the alternative method accumulate much slower than that of the original Lanczos method. There is almost no extra numeric cost for the gaining of the accuracy, i.e., the accuracy increase needs no extra operations of the Hamiltonian acting on state vectors, which are the major numeric cost for wave function propagation. A wave packet moving in a 2-dimensional Hénon-Heiles model serves as an illustration. This method is suitable for small time step propagation of quantum wave functions in large scale time dependent calculations where the operations of the Hamiltonian acting on state vectors are expensive.
△ Less
Submitted 4 January, 2006;
originally announced January 2006.
-
Modification to the pre-factor of the semiclassical propagator
Authors:
Quanlin Jie,
Bambi Hu,
Baowen Li
Abstract:
We modify the pre-factor of the semiclassical propagator to improve its efficiency in practical implementations. The new pre-factor represents the smooth portion of an orbit's contribution, and leads to fast convergence in numerical calculations. As an illustration of the accuracy and efficiency of the resultant propagator, we numerically calculate overlaps between quantum and semiclassical wave…
▽ More
We modify the pre-factor of the semiclassical propagator to improve its efficiency in practical implementations. The new pre-factor represents the smooth portion of an orbit's contribution, and leads to fast convergence in numerical calculations. As an illustration of the accuracy and efficiency of the resultant propagator, we numerically calculate overlaps between quantum and semiclassical wave functions, as well as low-lying spectrum density in a 10-dimensional system contains unstable classical orbits. This sheds light on applying semiclassical propagator to high dimensional systems.
△ Less
Submitted 8 March, 2005;
originally announced March 2005.
-
Modified conjugated gradient method for diagonalising large matrices
Authors:
Quanlin Jie,
Dunhuan Liu
Abstract:
We present an iterative method to diagonalise large matrices. The basic idea is the same as the conjugated gradient (CG) method, i.e, minimizing the Rayleigh quotient via its gradient and avoiding reintroduce errors to the directions of previous gradients. Each iteration step is to find lowest eigenvector of the matrix in a subspace spanned by the current trial vector and the corresponding gradi…
▽ More
We present an iterative method to diagonalise large matrices. The basic idea is the same as the conjugated gradient (CG) method, i.e, minimizing the Rayleigh quotient via its gradient and avoiding reintroduce errors to the directions of previous gradients. Each iteration step is to find lowest eigenvector of the matrix in a subspace spanned by the current trial vector and the corresponding gradient of the Rayleigh quotient, as well as some previous trial vectors. The gradient, together with the previous trail vectors, play a similar role of the conjugated gradient of the original CG algorithm. Our numeric tests indicate that this method converges significantly faster than the original CG method. And the computational cost of one iteration step is about the same as the original CG method. It is suitably for first principle calculations.
△ Less
Submitted 20 September, 2003;
originally announced September 2003.