Solving Dirac equations on a 3D lattice with inverse Hamiltonian and spectral methods

A new method to solve the Dirac equation on a 3D lattice is proposed, in which the variational collapse problem is avoided by the inverse Hamiltonian method and the fermion doubling problem is avoided by performing spatial derivatives in momentum space with the help of the discrete Fourier transform, i.e., the spectral method. This method is demonstrated in solving the Dirac equation for a given spherical potential in 3D lattice space. In comparison with the results obtained by the shooting method, the differences in single particle energy are smaller than $10^{-4}$~MeV, and the densities are almost identical, which demonstrates the high accuracy of the present method. The results obtained by applying this method without any modification to solve the Dirac equations for an axial deformed, non-axial deformed, and octupole deformed potential are provided and discussed.Comment: 18 pages, 6 figure

Covariant description of shape evolution and shape coexistence in neutron-rich nuclei at N\approx60

The shape evolution and shape coexistence phenomena in neutron-rich nuclei at $N\approx60$, including Kr, Sr, Zr, and Mo isotopes, are studied in the covariant density functional theory (DFT) with the new parameter set PC-PK1. Pairing correlations are treated using the BCS approximation with a separable pairing force. Sharp rising in the charge radii of Sr and Zr isotopes at N=60 is observed and shown to be related to the rapid changing in nuclear shapes. The shape evolution is moderate in neighboring Kr and Mo isotopes. Similar as the results of previous Hartree-Fock-Bogogliubov (HFB) calculations with the Gogny force, triaxiality is observed in Mo isotopes and shown to be essential to reproduce quantitatively the corresponding charge radii. In addition, the coexistence of prolate and oblate shapes is found in both $^{98}$Sr and $^{100}$Zr. The observed oblate and prolate minima are related to the low single-particle energy level density around the Fermi surfaces of neutron and proton respectively. Furthermore, the 5-dimensional (5D) collective Hamiltonian determined by the calculations of the PC-PK1 energy functional is solved for $^{98}$Sr and $^{100}$Zr. The resultant excitation energy of $0^+_2$ state and E0 transition strength $\rho^2(E0;0^+_2\rightarrow0^+_1)$ are in rather good agreement with the data. It is found that the lower barrier height separating the two competing minima along the $\gamma$ deformation in $^{100}$Zr gives rise to the larger $\rho^2(E0;0^+_2\rightarrow0^+_1)$ than that in $^{98}$Sr.Comment: 1 table, 11 figures, 23 page

An Efficient Method for GPS Multipath Mitigation Using the Teager-Kaiser-Operator-based MEDLL

An efficient method for GPS multipath mitigation is proposed. The motivation for this proposed method is to integrate the Teager-Kaiser Operator (TKO) with the Multipath Estimating Delay Lock Loop (MEDLL) module to mitigate the GPS multipath efficiently. The general implementation process of the proposed method is that we first utilize the TKO to operate on the received signal’s Auto-Correlation Function (ACF) to get an initial estimate of the multipaths. Then we transfer the initial estimated results to the MEDLL module for a further estimation. Finally, with a few iterations which are less than those of the original MEDLL algorithm, we can get a more accurate estimate of the Line-Of-Sight (LOS) signal, and thus the goal of the GPS multipath mitigation is achieved. The simulation results show that compared to the original MEDLL algorithm, the proposed method can reduce the computation load and the hardware and/or software consumption of the MEDLL module, meanwhile, without decreasing the algorithm accuracy

Direct Instantaneous Torque and Axial Force Control Method for Linear-Rotary Switched Reluctance Motor with Two Radial Windings

IEEEDual-winding linear-rotary switched reluctance motors (LRSRMs) suffer from large torque ripple and severe coupling between the torque winding and the axial-force winding. To address these issues, this paper proposes a direct instantaneous torque and direct axial force control (DITC&DAFC) method to suppress torque ripple and reduce the impact of the coupling between two sets of windings. The DITC&DAFC method divides inductance-rising zone and inductance-falling zone according to the inductance characteristics, and uses hysteresis control to directly control the motor\u27s instantaneous torque in different intervals. Therefore, the generation of negative torque is reduced, which effectively suppresses torque ripples. Meanwhile, the method obviates the calculation of current and flux linkage, thereby alleviating the demands on the controller. In addition, the mechanical structure and operating mechanism of 6/4 pole LRSRM with two radial windings are described in detail. The feasibility of the proposed control method is verified through simulation and experimental results

Low-lying states in 30^{30}Mg: a beyond relativistic mean-field investigation

The recently developed model of three-dimensional angular momentum projection plus generator coordinate method on top of triaxial relativistic mean-field states has been applied to study the low-lying states of $^{30}$Mg. The effects of triaxiality on the low-energy spectra and E0 and E2 transitions are examined.Comment: 6 pages, 3 figures, 1 table, talk presented at the 17th nuclear physics conference "Marie and Pierre Curie" Kazimierz Dolny, 22-26th September 2010, Polan

Mean-field embedding of the dual fermion approach for correlated electron systems

To reduce the rapidly growing computational cost of the dual fermion lattice calculation with increasing system size, we introduce two embedding schemes. One is the real fermion embedding, and the other is the dual fermion embedding. Our numerical tests show that the real fermion and dual fermion embedding approaches converge to essentially the same result. The application on the Anderson disorder and Hubbard models shows that these embedding algorithms converge more quickly with system size as compared to the conventional dual fermion method, for the calculation of both single-particle and two-particle quantities.Comment: 10 pages, 10 figure

Dual Fermion Method for Disordered Electronic Systems

While the coherent potential approximation (CPA) is the prevalent method for the study of disordered electronic systems, it fails to capture non-local correlations and Anderson localization. To incorporate such effects, we extend the dual fermion approach to disordered non-interacting systems using the replica method. Results for single- and two- particle quantities show good agreement with cluster extensions of the CPA; moreover, weak localization is captured. As a natural extension of the CPA, our method presents an alternative to the existing cluster theories. It can be used in various applications, including the study of disordered interacting systems, or for the description of non-local effects in electronic structure calculations.Comment: 5 pages, 4 figure