Source author record

Jie Meng

Jie Meng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

65works
21topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

65 published item(s)

preprint2026arXiv

Strain Engineering of Intrinsic Anomalous Hall and Nernst Effects in Altermagnetic MnTe at Realistic Doping Levels

Hexagonal MnTe has emerged as a prototypical g-wave altermagnet, hosting time-reversal symmetry breaking in momentum space despite a vanishing net magnetization. While this symmetry breaking theoretically allows for an intrinsic anomalous Hall effect, experimentally observed signals have remained weak. In this work, we investigate the origin of this suppression and demonstrate a strategy to amplify anomalous transport responses within the experimentally accessible doping regime. Using a $\bm{k}\cdot\bm{p}$ effective model, we reveal that near the valence band maximum, which corresponds to the energy window relevant for typical hole doping ($\sim10^{19}cm^{-3}$), the intrinsic Hall effect is suppressed due to a symmetry-enforced cancellation of opposing Berry curvature contributions. We propose that breaking the crystalline symmetry via volume-conserving biaxial strain lifts this cancellation, resulting in a significant enhancement of the anomalous Hall conductivity by orders of magnitude. This strain-induced Fermi surface distortion also amplifies the anomalous Nernst effect. Furthermore, the analysis of the spin texture confirms that these strain-enabled anomalous transport signatures emerge while preserving the zero net magnetization.

preprint2022arXiv

Accurate relativistic chiral nucleon-nucleon interaction up to NNLO

We construct a relativistic chiral nucleon-nucleon interaction up to the next-to-next-to-leading order in covariant baryon chiral perturbation theory. We show that a good description of the $np$ phase shifts up to $T_\mathrm{lab}=200$ MeV and even higher can be achieved with a $\tildeχ^2/\mathrm{d.o.f.}$ less than 1. Both the next-to-leading order results and the next-to-next-to-leading order results describe the phase shifts equally well up to $T_\mathrm{lab}=200$ MeV, but for higher energies, the latter behaves better, showing satisfactory convergence. The relativistic chiral potential provides the most essential inputs for relativistic ab initio studies of nuclear structure and reactions, which has been in need for almost two decades.

preprint2022arXiv

Asymmetric Nuclear Matter and Neutron Star Properties in Relativistic ab initio Theory in the Full Dirac Space

The long-standing controversy about the isospin dependence of the effective Dirac mass in ab initio calculations of asymmetric nuclear matter is clarified by solving the relativistic Brueckner-Hartree-Fock equations in the full Dirac space. The symmetry energy and its slope parameter at the saturation density are $E_{\text{sym}}(ρ_0)=33.1$ MeV and $L=65.2$ MeV, in agreement with empirical and experimental values. Further applications predict the neutron star radius $R_{1.4M_\odot}\approx 12$ km and the maximum mass of a neutron star $M_{\text{max}}\leq 2.4M_\odot$.

preprint2022arXiv

Finite amplitude method on the deformed relativistic Hartree-Bogoliubov theory in continuum: The isoscalar giant monopole resonance in exotic nuclei

Finite amplitude method based on the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc-FAM) is developed and applied to study isoscalar giant monopole resonance in exotic nuclei. Validation of the numerical implementation is examined for $^{208}\textrm{Pb}$. The isoscalar giant monopole resonances for even-even calcium isotopes from $^{40}\textrm{Ca}$ to the last bound neutron-rich nucleus $^{80}\textrm{Ca}$ are calculated, and a good agreement with the available experimental centroid energies is obtained for $^{40-48}\textrm{Ca}$. For the exotic calcium isotopes, e.g., $^{68}\textrm{Ca}$ and $^{80}\textrm{Ca}$, the DRHBc-FAM calculated results are closer to the energy weighted sum rule than the calculations on the harmonic oscillator basis, which highlights the advantages of DRHBc-FAM in describing giant resonances for exotic nuclei. In order to explore the soft monopole mode in the exotic nuclei, the giant monopole resonance for the deformed exotic nucleus $^{200}\textrm{Nd}$ is investigated, where the prolate shape and the oblate shape coexist. A soft monopole mode near 6.0 MeV is found in the prolate case, and another one near 4.5 MeV is found in the oblate case. The transition density of the soft monopole mode shows in phase or out-of-phase vibrations near the surface region, which is generated by quadrupole vibrations.

preprint2022arXiv

Massively scalable stencil algorithm

Stencil computations lie at the heart of many scientific and industrial applications. Unfortunately, stencil algorithms perform poorly on machines with cache based memory hierarchy, due to low re-use of memory accesses. This work shows that for stencil computation a novel algorithm that leverages a localized communication strategy effectively exploits the Cerebras WSE-2, which has no cache hierarchy. This study focuses on a 25-point stencil finite-difference method for the 3D wave equation, a kernel frequently used in earth modeling as numerical simulation. In essence, the algorithm trades memory accesses for data communication and takes advantage of the fast communication fabric provided by the architecture. The algorithm -- historically memory bound -- becomes compute bound. This allows the implementation to achieve near perfect weak scaling, reaching up to 503 TFLOPs on WSE-2, a figure that only full clusters can eventually yield.

preprint2022arXiv

Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum: I. even-even nuclei

Ground-state properties of even-even nuclei with $8\le Z\le120$ from the proton drip line to the neutron drip line have been investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the density functional PC-PK1. With the effects of deformation and continuum included simultaneously, 2583 even-even nuclei are predicted to be bound. The calculated binding energies, two-nucleon separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. The rms deviation from the 637 mass data is 1.518 MeV, providing one of the best microscopic descriptions for nuclear masses. The drip lines obtained from DRHBc calculations are compared with other calculations, including the spherical relativistic continuum Hartree-Bogoliubov (RCHB) and triaxial relativistic Hartree-Bogoliubov (TRHB) calculations with PC-PK1. The deformation and continuum effects on the limits of the nuclear landscape are discussed. Possible peninsulas consisting of bound nuclei beyond the two-neutron drip line are predicted. The systematics of the two-nucleon separation energies, two-nucleon gaps, rms radii, quadrupole deformations, potential energy curves, neutron densities, neutron mean-field potentials, and pairing energies in the DRHBc calculations are also discussed. In addition, the $α$ decay energies extracted are in good agreement with available data.

preprint2022arXiv

Triaxiality induced monopole-quadrupole-hexadecupole coupling in the isoscalar giant resonances for $^{86}\textrm{Ge}$

The isoscalar giant resonances for $^{86}\textrm{Ge}$ are studied by the quasiparticle finite amplitude method based on the covariant density functional theory. In addition to the well-known monopole-quadrupole coupling that splits the isoscalar giant monopole resonance in axially deformed nuclei, a monopole-quadrupole-hexadecupole coupling is identified in the triaxially deformed nucleus $^{86}\textrm{Ge}$, leading to the emergence of a distinct resonance peak at the low energy side of the isoscalar monopole strength function. The transition density of the triaxiality induced resonance peak shows a strong interplay among monopole, quadrupole, and hexadecupole vibrations. The resonance peak responses to monopole, quadrupole, and hexadecupole perturbations simultaneously, which could be regarded as a fingerprint of the triaxiality in $^{86}\textrm{Ge}$.

preprint2021arXiv

Geometric means of quasi-Toeplitz matrices

We study means of geometric type of quasi-Toeplitz matrices, that are semi-infinite matrices $A=(a_{i,j})_{i,j=1,2,\ldots}$ of the form $A=T(a)+E$, where $E$ represents a compact operator, and $T(a)$ is a semi-infinite Toeplitz matrix associated with the function $a$, with Fourier series $\sum_{\ell=-\infty}^{\infty} a_\ell e^{\mathfrak i \ell t}$, in the sense that $(T(a))_{i,j}=a_{j-i}$. If $a$ is \rv\ and essentially bounded, then these matrices represent bounded self-adjoint operators on $\ell^2$. We consider the case where $a$ is a continuous function, where quasi-Toeplitz matrices coincide with a classical Toeplitz algebra, and the case where $a$ is in the Wiener algebra, that is, has absolutely convergent Fourier series. We prove that if $a_1,\ldots,a_p$ are continuous and positive functions, or are in the Wiener algebra with some further conditions, then means of geometric type, such as the ALM, the NBMP and the Karcher mean of quasi-Toeplitz positive definite matrices associated with $a_1,\ldots,a_p$, are quasi-Toeplitz matrices associated with the geometric mean $(a_1\cdots a_p)^{1/p}$, which differ only by the compact correction. We show by numerical tests that these operator means can be practically approximated.

preprint2020arXiv

Accelerating High-Order Stencils on GPUs

Stencil computations are widely used in HPC applications. Today, many HPC platforms use GPUs as accelerators. As a result, understanding how to perform stencil computations fast on GPUs is important. While implementation strategies for low-order stencils on GPUs have been well-studied in the literature, not all of proposed enhancements work well for high-order stencils, such as those used for seismic modeling. Furthermore, coping with boundary conditions often requires different computational logic, which complicates efficient exploitation of the thread-level parallelism on GPUs. In this paper, we study high-order stencils and their unique characteristics on GPUs. We manually crafted a collection of implementations of a 25-point seismic modeling stencil in CUDA and related boundary conditions. We evaluate their code shapes, memory hierarchy usage, data-fetching patterns, and other performance attributes. We conducted an empirical evaluation of these stencils using several mature and emerging tools and discuss our quantitative findings. Among our implementations, we achieve twice the performance of a proprietary code developed in C and mapped to GPUs using OpenACC. Additionally, several of our implementations have excellent performance portability.

preprint2020arXiv

CasGCN: Predicting future cascade growth based on information diffusion graph

Sudden bursts of information cascades can lead to unexpected consequences such as extreme opinions, changes in fashion trends, and uncontrollable spread of rumors. It has become an important problem on how to effectively predict a cascade' size in the future, especially for large-scale cascades on social media platforms such as Twitter and Weibo. However, existing methods are insufficient in dealing with this challenging prediction problem. Conventional methods heavily rely on either hand crafted features or unrealistic assumptions. End-to-end deep learning models, such as recurrent neural networks, are not suitable to work with graphical inputs directly and cannot handle structural information that is embedded in the cascade graphs. In this paper, we propose a novel deep learning architecture for cascade growth prediction, called CasGCN, which employs the graph convolutional network to extract structural features from a graphical input, followed by the application of the attention mechanism on both the extracted features and the temporal information before conducting cascade size prediction. We conduct experiments on two real-world cascade growth prediction scenarios (i.e., retweet popularity on Sina Weibo and academic paper citations on DBLP), with the experimental results showing that CasGCN enjoys a superior performance over several baseline methods, particularly when the cascades are of large scale.

preprint2020arXiv

Deformed relativistic Hartree-Bogoliubov theory in continuum with point coupling functional: examples of even-even Nd isotopes

The aim of this work is to develop the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) theory based on the point-coupling density functionals and extend it to provide a unified description for all even-even nuclei in the nuclear chart by overcoming all possible challenges. The nuclear superfluidity is considered via Bogoliubov transformation. Densities and potentials are expanded in terms of Legendre polynomials to include the axial deformation degrees of freedom. Sophisticated relativistic Hartree-Bogoliubov equations in coordinate space are solved in the DiracWoods-Saxon basis to consider the continuum effects. Numerical checks are performed from light nuclei to heavy nuclei. The techniques to construct the DRHBc mass table for even-even nuclei are explored. The DRHBc theory is extended to study heavier nuclei beyond magnesium isotopes. Taking Nd isotopes as examples, the experimental binding energies, two-neutron separation energies, quadrupole deformations, and charge radii are reproduced rather well. The deformation and continuum play essential roles in the description of nuclear masses and prediction of drip-line nuclei. By examining the single-particle levels in the canonical basis and their contributions to the total density, the thickness of the neutron skin, the particles number in continuum, and the Coulomb barrier, the exotic structures including the neutron skin and the proton radioactivity are predicted.

preprint2020arXiv

Minimod: A Finite Difference solver for Seismic Modeling

This article introduces a benchmark application for seismic modeling using finite difference method, which is namedMiniMod, a mini application for seismic modeling. The purpose is to provide a benchmark suite that is, on one hand easy to build and adapt to the state of the art in programming models and changing high performance hardware landscape. On the other hand, the intention is to have a proxy application to actual production geophysical exploration workloads for Oil & Gas exploration, and other geosciences applications based on the wave equation. From top to bottom, we describe the design concepts, algorithms, code structure of the application, and present the benchmark results on different current computer architectures.

preprint2019arXiv

Solving quadratic matrix equations arising in random walks in the quarter plane

Quadratic matrix equations of the kind $A_1X^2+A_0X+A_{-1}=X$ are encountered in the analysis of Quasi--Birth-Death stochastic processes where the solution of interest is the minimal nonnegative solution $G$. In many queueing models, described by random walks in the quarter plane, the coefficients $A_1,A_0,A_{-1}$ are infinite tridiagonal matrices with an almost Toeplitz structure. Here, we analyze some fixed point iterations, including Newton's iteration, for the computation of $G$ and introduce effective algorithms and acceleration strategies which fully exploit the Toeplitz structure of the matrix coefficients and of the current approximation. Moreover, we provide a structured perturbation analysis for the solution $G$. The results of some numerical experiments which demonstrate the effectiveness of our approach are reported.

preprint2019arXiv

Symmetry energy at supra-saturation densities via the Gravitational Waves from GW170817

Motivated by the historical detection of gravitational waves from GW170817, the neutron star and the neutron drop, i.e., a certain number of neutrons confined in an external field, are systematically investigated by ab initio calculations as well as the nonrelativistic and relativistic state-of-art density functional theories. Strong correlations are found among the neutron star tidal deformability, the neutron star radius, the root-mean-square radii of neutron drops, and the symmetry energies of nuclear matter at supra-saturation densities. From these correlations and the upper limit on the tidal deformability extracted from GW170817, the neutron star radii, the neutron drop radii, and the symmetry energy at twice saturation density are respectively constrained as $R_{1.4M_{\odot}}\leqslant 12.94$ km, $R_{\rm nd} \leqslant 2.36$ fm, and $E_{\mathrm{sym}}(2ρ_0) \leqslant 53.2$ MeV.

preprint2016arXiv

Nuclear chiral and magnetic rotation in covariant density functional theory

Excitations of chiral rotation observed in triaxial nuclei and magnetic and/or antimagnetic rotations seen in near-spherical nuclei have attracted a lot of attention. Unlike conventional rotation in well-deformed or superdeformed nuclei, here the rotational axis is not necessary coinciding with any principal axis of the nuclear density distribution. Thus, tilted axis cranking is mandatory to describe these excitations self-consistently in the framework of covariant density functional theory (CDFT). We will briefly introduce the formalism of tilted axis cranking CDFT and its application for magnetic and antimagnetic rotation phenomena. Configuration-fixed CDFT and its predictions for nuclear chiral configurations and for favorable triaxial deformation parameters are also presented, and the discoveries of the multiple chiral doublets (M\c{hi}D) in 133Ce and 103Rh are discussed.

preprint2016arXiv

Octet baryon masses and sigma terms in covariant baryon chiral perturbation theory

We report on a recent study of the ground-state octet baryon masses and sigma terms in covariant baryon chiral perturbation theory with the extended-on-mass-shell scheme up to next-to-next-to-next-to-leading order. To take into account lattice QCD artifacts, the finite-volume corrections and finite lattice spacing discretization effects are carefully examined. We performed a simultaneous fit of all the $n_f = 2+1$ lattice octet baryon masses and found that the various lattice simulations are consistent with each other. Although the finite lattice spacing discretization effects up to $\mathcal{O}(a^2)$ can be safely ignored, but the finite volume corrections cannot even for configurations with $M_ϕL>4$. As an application, we predicted the octet baryon sigma terms using the Feynman-Hellmann theorem. In particular, the pion- and strangeness-nucleon sigma terms are found to be $σ_{πN} = 55(1)(4)$ MeV and $σ_{sN} = 27(27)(4)$ MeV, respectively.

preprint2015arXiv

Correlations between neutrons and protons near Fermi surface and $Q_α$ of super-heavy nuclei

The shell corrections and shell gaps in nuclei are systematically studied with the latest Weizsäcker-Skyrme (WS4) mass model. We find that most of asymmetric nuclei with (sub)-shell closures locate along the shell stability line (SSL), $N=1.37Z+13.5$, which might be due to a strong correlation between neutrons and protons near Fermi surface. The double magicity of nuclei $^{46}$Si and $^{78}$Ni is predicted according to the corresponding shell gaps, shell corrections and nuclear deformations. The unmeasured super-heavy nuclei $^{296}$118 and $^{298}$120, with relatively large shell gaps and shell corrections, also locate along the SSL, whereas the traditional magic nucleus $^{298}$Fl evidently deviates from the line. The $α$-decay energies of super-heavy nuclei with $Z=113-126$ are simultaneously investigated by using the WS4 model together with the radial basis function corrections. For super-heavy nuclei with large shell corrections, the smallest $α$-decay energy for elements $Z=116$, 117 and 118 in their isotope chains locates at $N=178$ rather than $184$.

preprint2015arXiv

Halos in medium-heavy and heavy nuclei with covariant density functional theory in continuum

The covariant density functional theory with a few number of parameters has been widely used to describe the ground-state and excited-state properties for the nuclei all over the nuclear chart. In order to describe exotic properties of unstable nuclei, the contribution of the continuum and its coupling with bound states should be treated properly. In this Topical Review, the development of the covariant density functional theory in continuum will be introduced, including the relativistic continuum Hartree-Bogoliubov theory, the relativistic Hartree-Fock-Bogoliubov theory in continuum, and the deformed relativistic Hartree-Bogoliubov theory in continuum. Then the descriptions and predictions of the neutron halo phenomena in both spherical and deformed nuclei will be reviewed. The diffuseness of the nuclear potentials, nuclear shapes and density distributions, and the impact of the pairing correlations on nuclear size will be discussed.

preprint2015arXiv

Scalar strangeness content of the nucleon and baryon sigma terms

The scalar strangeness content of the nucleon, characterized by the so-called strangeness-nucleon sigma term, is of fundamental importance in understanding its sea-quark flavor structure. We report a determination of the octet baryon sigma terms via the Feynman-Hellmann theorem by analyzing the latest high-statistics $n_f=2+1$ lattice QCD simulations with covariant baryon chiral perturbation theory up to next-to-next-to-next-to-leading order. In particular, we predict $σ_{πN}=55(1)(4)$ MeV and $σ_{sN}=27(27)(4)$ MeV, while the first error is statistical and the second systematic due to different lattice scales. The predicted $σ_{sN}$ is consistent with the latest LQCD results and the results based on the next-to-next-to-leading order chiral perturbation theory. Several key factors in determining the sigma terms are systematically taken into account and clarified for the first time, including the effects of lattice scale setting, systematic uncertainties originating from chiral expansion truncations, and constraint of strong-interaction isospin breaking effects.

preprint2014arXiv

Baryon chiral perturbation theory with Wilson fermions up to $\mathcal{O}(a^2)$ and discretization effects of latest $n_f=2+1$ LQCD octet baryon masses

We construct the chiral Lagrangians relevant in studies of the ground-state octet baryon masses up to $\mathcal{O}(a^2)$ by taking into account discretization effects and calculate the masses up to $\mathcal{O}(p^4)$ in the extended-on-mass-shell scheme. As an application, we study the latest $n_f=2+1$ LQCD data on the ground-state octet baryon masses from the PACS-CS, QCDSF-UKQCD, HSC, and NPLQCD Collaborations. It is shown that the discretization effects for the studied LQCD simulations are at the order of one to two percent for lattice spacings up to $0.15$ fm and the pion mass up to 500 MeV.

preprint2014arXiv

Decuplet baryon masses in covariant baryon chiral perturbation theory

We present an analysis of the lowest-lying decuplet baryon masses in the covariant baryon chiral perturbation theory with the extended-on-mass-shell scheme up to next-to-next-to-next-to-leading order. In order to determine the $14$ low-energy constants, we perform a simultaneous fit of the $n_f=2+1$ lattice QCD data from the PACS-CS, QCDSF-UKQCD, and HSC Collaborations, taking finite-volume corrections into account self-consistently. We show that up to next-to-next-to-next-to-leading order one can achieve a good description of the lattice QCD and experimental data. Surprisingly, we note that the current lattice decuplet baryon masses can be fitted rather well by the next-to-leading order baryon chiral perturbation theory, which, however, misses the experimental data a little bit. Furthermore, we predict the pion- and strangeness-sigma terms of the decuplet baryons by use of the Feynman-Hellmann theorem.

preprint2014arXiv

Hidden pseudospin and spin symmetries and their origins in atomic nuclei

Symmetry plays a fundamental role in physics. The quasi-degeneracy between single-particle orbitals $(n, l, j = l + 1/2)$ and $(n-1, l + 2, j = l + 3/2)$ indicates a hidden symmetry in atomic nuclei, the so-called pseudospin symmetry (PSS). Since the introduction of the concept of PSS in atomic nuclei, there have been comprehensive efforts to understand its origin. Both splittings of spin doublets and pseudospin doublets play critical roles in the evolution of magic numbers in exotic nuclei discovered by modern spectroscopic studies with radioactive ion beam facilities. Since the PSS was recognized as a relativistic symmetry in 1990s, many special features, including the spin symmetry (SS) for anti-nucleon, and many new concepts have been introduced. In the present Review, we focus on the recent progress on the PSS and SS in various systems and potentials, including extensions of the PSS study from stable to exotic nuclei, from non-confining to confining potentials, from local to non-local potentials, from central to tensor potentials, from bound to resonant states, from nucleon to anti-nucleon spectra, from nucleon to hyperon spectra, and from spherical to deformed nuclei. Open issues in this field are also discussed in detail, including the perturbative nature, the supersymmetric representation with similarity renormalization group, and the puzzle of intruder states.

preprint2014arXiv

Surface diffuseness correction in global mass formula

By taking into account the surface diffuseness correction for unstable nuclei, the accuracy of the macroscopic-microscopic mass formula is further improved. The rms deviation with respect to essentially all the available mass data falls to 298 keV, crossing the 0.3 MeV accuracy threshold for the first time within the mean-field framework. Considering the surface effect of the symmetry potential which plays an important role in the evolution of the "neutron skin" toward the "neutron halo" of nuclei approaching the neutron drip line, we obtain an optimal value of the symmetry energy coefficient J=30.16 MeV. With an accuracy of 258 keV for all the available neutron separation energies and of 237 keV for the alpha-decay Q-values of super-heavy nuclei, the proposed mass formula is particularly important not only for the reliable description of the r-process of nucleosynthesis but also for the study of the synthesis of super-heavy nuclei.

preprint2014arXiv

Test of the $h_1(1830)$ made of $K^{*}\bar{K}^*$ with the $η_c\rightarrow ϕK^{*}\bar{K}^{*}$ decay

We present a new reaction, complementary to the $J/ψ\rightarrow ηK^{*0}\bar{K}^{*0}$ from which an $h_1$ resonance with mass around $1830$ MeV was reported from a BESIII experiment. The new reaction is $η_c\rightarrow ϕK^*\bar{K}^*$, or $η_c(2S)\rightarrow ϕK^*\bar{K}^*$. Using the information from the analysis of $J/ψ\rightarrow ηK^{*0}\bar{K}^{*0}$, we find that the $K^*\bar{K}^*$ invariant mass distribution for those two $η_c$ decays exhibits a clear peak around $1830$ MeV perfectly distinguishable from what one obtains with pure phase space. We suggest the implementation of these reactions to assert the existence of this elusive resonance which, by its nature as a vector-vector molecule with $0^-(1^{+-})$ quantum numbers, only couples to the $K^*\bar{K}^*$ channel.

preprint2014arXiv

Treating Coulomb exchange contributions in relativistic mean field calculations: why and how

The energy density functional (EDF) method is very widely used in nuclear physics, and among the various existing functionals those based on the relativistic Hartree (RH) approximation are very popular because the exchange contributions (Fock terms) are numerically rather onerous to calculate. Although it is possible to somehow 'mock up' the effects of meson-induced exchange terms by adjusting the meson-nucleon couplings, the lack of Coulomb exchange contributions hampers the accuracy of predictions. In this note, we show that the Coulomb exchange effects can be easily included with a good accuracy in a perturbative approach. Therefore, it would be desirable for future relativistic EDF models to incorporate Coulomb exchange effects, at least to some order of perturbation.

preprint2013arXiv

Extending the nuclear chart by continuum: from oxygen to titanium

Nuclear masses ranging from O to Ti isotopes are systematically investigated with relativistic continuum Hartree-Bogoliubov (RCHB) theory, which can provide a proper treatment of pairing correlations in the presence of the continuum. From O to Ti isotopes, there are 402 nuclei predicted to be bound by the density functional PC-PK1. For the 234 nuclei with mass measured, the root mean square (rms) deviation is 2.23 MeV. It is found that the proton drip-lines predicted with various mass models are roughly the same and basically agree with the observation. The neutron drip-lines predicted, however, are quite different. Due to the continuum couplings, the neutron drip-line nuclei predicted are extended further neutron-rich than other mass models. By comparison with finite-range droplet model (FRDM), the neutron drip-line nucleus predicted by RCHB theory has respectively 2(O), 10(Ne), 10(Na), 6(Mg), 8(Al), 6(Si), 8(P), 6(S), 14(K), 10(Ca), 10(Sc), and 12(Ti) more neutrons.

preprint2013arXiv

Feasibility of the finite amplitude method in covariant density functional theory

Self-consistent relativistic random-phase approximation (RPA) in the radial coordinate representation is established by using the finite amplitude method (FAM). Taking the isoscalar giant monopole resonance in spherical nuclei as example, the feasibility of the FAM for the covariant density functionals is demonstrated, and the newly developed methods are verified by the conventional RPA calculations. In the present relativistic RPA calculations, the effects of the Dirac sea can be automatically taken into account in the coordinate-space representation. The rearrangement terms due to the density-dependent couplings can be implicitly calculated without extra computational costs in both iterative and matrix FAM schemes.

preprint2013arXiv

Finite-amplitude method: An extension to the covariant density functionals

The finite-amplitude method (FAM) is one of the most promising methods for optimizing the computational performance of the random-phase approximation (RPA) calculations in deformed nuclei. In this report, we will mainly focus on our recent progress in the self-consistent relativistic RPA established by using the FAM. It is found that the effects of Dirac sea can be taken into account implicitly in the coordinate-space representation and the rearrangement terms due to the density-dependent couplings can be treated without extra computational costs.

preprint2013arXiv

Lowest-lying octet baryon masses in covariant baryon chiral perturbation theory

We report on a systematic study of the ground-state octet baryon masses in the covariant baryon chiral perturbation theory with the extended-on-mass-shell renormalization scheme up to next-to-next-to-next-to-leading order, taking into account the contributions of the virtual decuplet baryons. A reasonable description of the lattice results is achieved by fitting simultaneously all the publicly available $n_f = 2+1$ lattice QCD data. It confirms that the various lattice simulations are consistent with each other. We stress that a self-consistent treatment of finite-volume corrections is important to obtain a $χ^2/\rm{d.o.f.}$ about 1.

preprint2013arXiv

Microscopic and self-consistent description for neutron halo in deformed nuclei

A deformed relativistic Hartree-Bogoliubov theory in continuum has been developed for the study of neutron halos in deformed nuclei and the halo phenomenon in deformed weakly bound nuclei is investigated. Magnesium and neon isotopes are studied and some results are presented for the deformed neutron-rich and weakly bound nuclei 44Mg and 36Ne. The core of the former nucleus is prolate, but the halo has a slightly oblate shape. This indicates a decoupling of the halo orbitals from the deformation of the core. The generic conditions for the existence of halos in deformed nuclei and for the occurrence of this decoupling effect are discussed.

preprint2013arXiv

Nuclear superfluidity for antimagnetic rotation in $^{105}$Cd and $^{106}$Cd

The effect of nuclear superfluidity on antimagnetic rotation bands in $^{105}$Cd and $^{106}$Cd are investigated by the cranked shell model with the pairing correlations and the blocking effects treated by a particle-number conserving method. The experimental moments of inertia and the reduced $B(E2)$ transition values are excellently reproduced. The nuclear superfluidity is essential to reproduce the experimental moments of inertia. The two-shears-like mechanism for the antimagnetic rotation is investigated by examining the shears angle, i.e., the closing of the two proton hole angular momenta, and its sensitive dependence on the nuclear superfluidity is revealed.

preprint2013arXiv

Pairing strength on the nuclear size in relativistic continuum Hartree-Bogoliubov theory

The influence of pairing correlations on the nuclear size and in particular on the formation of nuclear halos is studied in the framework of relativistic continuum Hartree-Bogoliubov (RCHB) theory. It turns out that the contributions from the weakly-bound orbits with low orbital angular momenta $l$ play an important role. As an example, we investigate the neutron-rich Mg isotopes as a function of the pairing strength in situations, where the neuron Fermi surface are below, between and above two weakly-bound $2p$ levels. We find that the size of the pairing correlations has a two-fold influence on the density distribution of the neutrons and therefore on the total nuclear size. First it can change the root mean square (rms) radius of the individual weakly-bound orbits and second it can enhance the occupation probabilities of these orbits in the nuclear system. On one side increasing pairing correlations reduce the rms radii of the orbits with small orbital angular momenta close to the continuum limit (pairing anti-halo effect), on the other side they also can lead to an enhanced occupation of low-$l$ orbits above the Fermi surface producing in this way a strong increase of the total radius of the nuclear system. As a consequence, a nuclear halo can form even in cases, where non of the individual low-$l$ orbits is very close to the continuum. This leads to the fact that compared with well-bound nuclei, the impact of the pairing strength on the nuclear size is more pronounced in weakly-bound nuclei than in well-bound systems.

preprint2013arXiv

Pseudospin symmetry in supersymmetric quantum mechanics: II. Spin-orbit effects

Following a previous paper [Haozhao Liang \textit{et al.}, Phys. Rev. C \textbf{87}, 014334 (2013)], we discuss the spin-orbit effects on the pseudospin symmetry (PSS) within the framework of supersymmetric quantum mechanics. By using the perturbation theory, we demonstrate that the perturbative nature of PSS maintains when a substantial spin-orbit potential is included. With the explicit PSS-breaking potential, the spin-orbit effects on the pseudospin-orbit splittings are investigated in a quantitative way.

preprint2013arXiv

Pseudospin symmetry in supersymmetric quantum mechanics: Schrödinger equations

The origin of pseudospin symmetry (PSS) and its breaking mechanism are explored by combining supersymmetry (SUSY) quantum mechanics, perturbation theory, and the similarity renormalization group (SRG) method. The Schrödinger equation is taken as an example, corresponding to the lowest-order approximation in transforming a Dirac equation into a diagonal form by using the SRG. It is shown that while the spin-symmetry-conserving term appears in the single-particle Hamiltonian $H$, the PSS-conserving term appears naturally in its SUSY partner Hamiltonian $\tilde{H}$. The eigenstates of Hamiltonians $H$ and $\tilde{H}$ are exactly one-to-one identical except for the so-called intruder states. In such a way, the origin of PSS deeply hidden in $H$ can be traced in its SUSY partner Hamiltonian $\tilde{H}$. The perturbative nature of PSS in the present potential without spin-orbit term is demonstrated by the perturbation calculations, and the PSS-breaking term can be regarded as a very small perturbation on the exact PSS limits. A general tendency that the pseudospin-orbit splittings become smaller with increasing single-particle energies can also be interpreted in an explicit way.

preprint2013arXiv

Pseudospin symmetry: Recent progress with supersymmetric quantum mechanics

It is an interesting and open problem to trace the origin of the pseudospin symmetry in nuclear single-particle spectra and its symmetry breaking mechanism in actual nuclei. In this report, we mainly focus on our recent progress on this topic by combining the similarity renormalization group technique, supersymmetric quantum mechanics, and perturbation theory. We found that it is a promising direction to understand the pseudospin symmetry in a quantitative way.

preprint2013arXiv

Rotational properties of the superheavy nucleus 256Rf and its neighboring even-even nuclei in particle-number conserving cranked shell model

The ground state band was recently observed in the superheavy nucleus 256Rf. We study the rotational properties of 256Rf and its neighboring even-even nuclei by using a cranked shell model (CSM) with the pairing correlations treated by a particle-number conserving (PNC) method in which the blocking effects are taken into account exactly. The kinematic and dynamic moments of inertia of the ground state bands in these nuclei are well reproduced by the theory. The spin of the lowest observed state in 256Rf is determined by comparing the experimental kinematic moments of inertia with the PNC-CSM calculations and agrees with previous spin assignment. The effects of the high order deformation varepsilon6 on the angular momentum alignments and dynamic moments of inertia in these nuclei are discussed.

preprint2013arXiv

The Slater approximation for Coulomb exchange effects in nuclear covariant density functional theory

The relativistic local density approximation (LDA) for the Coulomb exchange functional in nuclear systems is presented. This approximation is composed of the well-known Slater approximation in the non-relativistic scheme and the corrections due to the relativistic effects. Its validity in finite nuclei is examined by comparing with the exact treatment of the Coulomb exchange term in the relativistic Hartree-Fock-Bogoliubov theory. The relativistic effects are found to be important and the exact Coulomb exchange energies can be reproduced by the relativistic LDA within 5% demonstrated by the semi-magic Ca, Ni, Zr, Sn, and Pb isotopes from proton drip line to neutron drip line.

preprint2013arXiv

Virtual decuplet effects on octet baryon masses in covariant baryon chiral perturbation theory

We extend a previous analysis of the lowest-lying octet baryon masses in covariant baryon chiral perturbation theory (ChPT) by explicitly taking into account the contribution of the virtual decuplet baryons. Up to next-to-next-to-next-to-leading order (N$^3$LO), the effects of these heavier degrees of freedom are systematically studied. Their effects on the light-quark mass dependence of the octet baryon masses are shown to be relatively small and can be absorbed by the available low-energy constants up to N$^3$LO. Nevertheless, a better description of the finite-volume corrections of the lattice QCD data can be achieved, particularly those with small $M_ϕL$ ($<4$), which is demonstrated by a careful study of the NPLQCD and QCDSF-UKQCD small-volume data. Finally, we show that the predicted pion- and strangeness-baryon sigma terms are only slightly changed by the inclusion of the virtual decuplet baryons.

preprint2012arXiv

BCS-BEC crossover in nuclear matter with the relativistic Hartree-Bogoliubov theory

Based on the relativistic Hartree-Bogoliubov theory, the influence of the pairing interaction strength on the di-neutron correlations and the crossover from superfluidity of neutron Cooper pairs in the $^{1}S_{0}$ channel to Bose-Einstein condensation of di-neutron pairs is systematically investigated in the nuclear matter. The bare nucleon-nucleon interaction Bonn-B is taken in the particle-particle channel with an effective factor to simulate the medium effects and take into account the possible ambiguity of pairing force, and the effective interaction PK1 is used in the particle-hole channel. If the effective factor is larger than 1.10, a di-neutron BEC state appears in the low-density limit, and if it is smaller than 0.85, the neutron Cooper pairs are found totally in the weak coupling BCS region. The reference values of several characteristic quantities which characterize the BCS-BEC crossover are obtained respectively from the dimensionless parameter $1/(k_{\rm Fn}a)$ with $a$ the scattering length and $k_{\rm{Fn}}$ the neutron Fermi momentum, the zero-momentum transfer density correlation function D(0) and the effective chemical potential $ν_{\rm n}$.

preprint2012arXiv

Deformed relativistic Hartree Bogoliubov theory in continuum

A deformed relativistic Hartree Bogoliubov (RHB) theory in continuum is developed aiming at a proper description of exotic nuclei, particularly those with a large spatial extension. In order to give an adequate consideration of both the contribution of the continuum and the large spatial distribution in exotic nuclei, the deformed RHB equations are solved in a Woods-Saxon (WS) basis in which the radial wave functions have a proper asymptotic behavior at large distance from the nuclear center. This is crucial for the proper description of a possible halo. The formalism of deformed RHB theory in continuum is presented. A stable nucleus, 20Mg and a weakly-bound nucleus, 42Mg, are taken as examples to present numerical details and to carry out necessary numerical checks. In addition, the ground state properties of even-even magnesium isotopes are investigated. The generic conditions of the formation of a halo in weakly bound deformed systems and the shape of the halo in deformed nuclei are discussed. We show that the existence and the deformation of a possible neutron halo depend essentially on the quantum numbers of the main components of the single particle orbitals in the vicinity of the Fermi surface.

preprint2012arXiv

Density-dependent deformed relativistic Hartree-Bogoliubov theory in continuum

The deformed relativistic Hartree-Bogoliubov theory in continuum with the density-dependent meson-nucleon couplings is developed. The formulism is briefly presented with the emphasis on handling the density-dependent couplings, meson fields, and potentials in axially deformed system with partial wave method. Taking the neutron-rich nucleus $^{38}$Mg as an example, the newly developed code is verified by the spherical relativistic continuum Hartree-Bogoliubov calculations, where only the spherical components of the densities are considered. When the deformation is included self-consistently, it is shown that the spherical components of density-dependent coupling strengths are dominant, while the contributions from low-order deformed components are not negligible.

preprint2012arXiv

Fine structure of charge-exchange spin-dipole excitations in $^{16}$O

The charge-exchange spin-dipole (SD) excitations for both $(p,n)$ and $(n,p)$ channels in $^{16}$O are investigated in the fully self-consistent random phase approximation based on the covariant density functional theory. The fine structure of SD excitations in the most up-to-date $^{16}$O($\vec p, \vec n$)$^{16}$F experiment is excellently reproduced without any readjustment in the functional. The characteristics of SD excitations are understood with the delicate balance between the $σ$- and $ω$-meson fields via the exchange terms. The fine structure of SD excitations for $^{16}$O($n,p$)$^{16}$N channel is predicted for future experiments.

preprint2012arXiv

Halos in a deformed Relativistic Hartree-Bogoliubov theory in continuum

In this contribution we present some recent results about neutron halos in deformed nuclei. A deformed relativistic Hartree-Bogoliubov theory in continuum has been developed and the halo phenomenon in deformed weakly bound nuclei is investigated. These weakly bound quantum systems present interesting examples for the study of the interdependence between the deformation of the core and the particles in the halo. Magnesium and neon isotopes are studied and detailed results are presented for the deformed neutron-rich and weakly bound nuclei 42Mg. The core of this nucleus is prolate, but the halo has a slightly oblate shape. This indicates a decoupling of the halo orbitals from the deformation of the core. The generic conditions for the existence of halos in deformed nuclei and for the occurrence of this decoupling effect are discussed.

preprint2012arXiv

Localized form of Fock terms in nuclear covariant density functional theory

In most of the successful versions of covariant density functional theory in nuclei, the Fock terms are not included explicitly, which leads to local functionals and forms the basis of their widespread applicability at present. However, it has serious consequences for the description of Gamow-Teller resonances (GTR) and spin-dipole resonances (SDR) which can only be cured by adding further phenomenological parameters. Relativistic Hartree-Fock models do not suffer from these problems. They can successfully describe the GTR and SDR as well as the isovector part of the Dirac effective mass without any additional parameters. However, they are non-local and require considerable numerical efforts. By the zero-range reduction and the Fierz transformation, a new method is proposed to take into account the Fock terms in local functionals, which retains the simplicity of conventional models and provides proper descriptions of the spin-isospin channels and the Dirac masses.

preprint2012arXiv

Odd systems in deformed relativistic Hartree Bogoliubov theory in continuum

In order to describe the exotic nuclear structure in unstable odd-$A$ or odd-odd nuclei, the deformed relativistic Hartree Bogoliubov theory in continuum has been extended to incorporate the blocking effect due to the odd nucleon. For a microscopic and self-consistent description of pairing correlations, continuum, deformation, blocking effects, and the extended spatial density distribution in exotic nuclei, the deformed relativistic Hartree Bogoliubov equations are solved in a Woods-Saxon basis in which the radial wave functions have a proper asymptotic behavior at large $r$. The formalism and numerical details are provided. The code is checked by comparing the results with those of spherical relativistic continuum Hartree Bogoliubov theory in the nucleus $^{19}$O. The prolate deformed nucleus $^{15}$C is studied by examining the neutron levels and density distributions.

preprint2012arXiv

The Gamow-Teller response within Skyrme random-phase approximation plus particle-vibration coupling

Although many random-phase approximation (RPA) calculations of the Gamow-Teller (GT) response exist, this is not the case for calculations going beyond the mean-field approximation. We apply a consistent model that includes the coupling of the GT resonance to low-lying vibrations, to nuclei of the $fp$ shell. Among other motivations, our goal is to see if the particle-vibration coupling can redistribute the low-lying GT$^+$ strength that is relevant for electron-capture processes in core-collapse supernova. We conclude that the lowering and fragmentation of that strength are consistent with the experimental findings and validate our model. However, the particle-vibration coupling cannot account for the quenching of the total value of the low-lying strength.

preprint2011arXiv

Lambda and Anti-Lambda Hypernuclei in Relativistic Mean-field Theory

Several aspects about $Λ$-hypernuclei in the relativistic mean field theory, including the effective $Λ$-nucleon coupling strengths based on the successful effective nucleon-nucleon interaction PK1, hypernuclear magnetic moment and $\barΛ$-hypernuclei, have been presented. The effect of tensor coupling in $Λ$-hypernuclei and the impurity effect of $\barΛ$ to nuclear structure have been discussed in detail.

preprint2011arXiv

Neutron halo in deformed nuclei from a relativistic Hartree-Bogoliubov model in a Woods-Saxon basis

Halo phenomenon in deformed nuclei is studied by using a fully self-consistent deformed relativistic Hartree-Bogoliubov model in a spherical Woods-Saxon basis with the proper asymptotic behavior at large distance from the nuclear center. Taking a deformed neutron-rich and weakly bound nucleus $^{44}$Mg as an example and by examining contributions of the halo, deformation effects, and large spatial extensions, we show a decoupling of the halo orbitals from the deformation of the core.

preprint2011arXiv

Perturbative interpretation of relativistic symmetries in nuclei

Perturbation theory is used systematically to investigate the symmetries of the Dirac Hamiltonian and their breaking in atomic nuclei. Using the perturbation corrections to the single-particle energies and wave functions, the link between the single-particle states in realistic nuclei and their counterparts in the symmetry limits is discussed. It is shown that the limit of S-V=const and relativistic harmonic oscillator (RHO) potentials can be connected to the actual Dirac Hamiltonian by the perturbation method, while the limit of S+V=const cannot, where S and V are the scalar and vector potentials, respectively. This indicates that the realistic system can be treated as a perturbation of spin-symmetric Hamiltonians, and the energy splitting of the pseudospin doublets can be regarded as a result of small perturbation around the Hamiltonian with RHO potentials, where the pseudospin doublets are quasidegenerate.

preprint2011arXiv

Stellar electron-capture rates calculated with the finite-temperature relativistic random-phase approximation

We introduce a self-consistent microscopic theoretical framework for modelling the process of electron capture on nuclei in stellar environment, based on relativistic energy density functionals. The finite-temperature relativistic mean-field model is used to calculate the single-nucleon basis and the occupation factors in a target nucleus, and $J^π = 0^{\pm}$, $1^{\pm}$, $2^{\pm}$ charge-exchange transitions are described by the self-consistent finite-temperature relativistic random-phase approximation. Cross sections and rates are calculated for electron capture on 54,56Fe and 76,78Ge in stellar environment, and results compared with predictions of similar and complementary model calculations.

preprint2010arXiv

Neutron halo in deformed nuclei

Halo phenomena in deformed nuclei are investigated within a deformed relativistic Hartree Bogoliubov (DRHB) theory. These weakly bound quantum systems present interesting examples for the study of the interdependence between the deformation of the core and the particles in the halo. Contributions of the halo, deformation effects, and large spatial extensions of these systems are described in a fully self-consistent way by the DRHB equations in a spherical Woods-Saxon basis with the proper asymptotic behavior at large distance from the nuclear center. Magnesium and neon isotopes are studied and detailed results are presented for the deformed neutron-rich and weakly bound nucleus $^{44}$Mg. The core of this nucleus is prolate, but the halo has a slightly oblate shape. This indicates a decoupling of the halo orbitals from the deformation of the core. The generic conditions for the occurence of this decoupling effects are discussed.

preprint2010arXiv

Nuclear halo structure and pseudo-spin symmetry

Nuclear halo structure and restoration of relativistic symmetry are studied within the framework of the relativistic Hartree-Fock-Bogoliubov (RHFB) theory. Giant halos as well as ordinary ones are found in Cerium isotopes close to the neutron drip line. Bridged by T=0 {channel}, the restoration of pseudo-spin symmetry (PSS) plays an essential role in stabilizing the neutron halo structures. The Fock terms, especially the $ρ$-tensor couplings, not only play significant role in the PSS restoration but also present substantial contributions to the T=0 {channel}, from which is well demonstrated the necessity of Fock terms.

preprint2010arXiv

Open Problems in Understanding the Nuclear Chirality

Open problems in the interpretation of the observed pair of near degenerate $ΔI = 1$ bands with the same parity as the chiral doublet bands are discussed. The ambiguities for the existing fingerprints of the chirality in atomic nuclei and problems in existing theory are discussed, including the description of quantum tunneling in the mean field approximation as well as the deformation, core polarization and configuration of particle rotor model (PRM). Future developments of the theoretical approach are prospected.

preprint2010arXiv

Persistent contribution of unbound quasiparticles to the pair correlation in continuum Skyrme-Hartree-Fock-Bogoliubov approach

The neutron pair correlation in nuclei near the neutron drip-line is investigated using the selfconsistent continuum Skyrme-Hartree-Fock-Bogoliubov theory formulated with the coordinate-space Green's function technique. Numerical analysis is performed for even-even N = 86 isotones in the Mo-Sn region, where the 3p3/2 and 3p1/2 orbits lying near the Fermi energy are either weakly bound or unbound. The quasiparticle states originating from the l = 1 orbits form resonances with large widths, which are due to the low barrier height and the strong continuum coupling caused by the pair potential. Analyzing in detail the pairing properties and roles of the quasiparticle resonances, we found that the l = 1 broad quasiparticle resonances persist to feel the pair potential and contribute to the pair correlation even when their widths are comparable with the resonance energy.

preprint2010arXiv

Spin symmetry in Dirac negative energy spectrum in density-dependent relativistic Hartree-Fock theory

The spin symmetry in the Dirac negative energy spectrum and its origin are investigated for the first time within the density-dependent relativistic Hartree-Fock (DDRHF) theory. Taking the nucleus $^{16}$O as an example, the spin symmetry in the negative energy spectrum is found to be a good approximation and the dominant components of the Dirac wave functions for the spin doublets are nearly identical. In comparison with the relativistic Hartree approximation where the origin of spin symmetry lies in the equality of the scalar and vector potentials, in DDRHF the cancellation between the Hartree and Fock terms is responsible for the better spin symmetry properties and determines the subtle spin-orbit splitting. These conclusions hold even in the case when significant deviations from the G-parity values of the meson-antinucleon couplings occur.

preprint2010arXiv

Spin-orbit and orbit-orbit strengths for radioactive neutron-rich doubly magic nucleus $^{132}$Sn in relativistic mean field theory

Relativistic mean field (RMF) theory is applied to investigate the properties of the radioactive neutron-rich doubly magic nucleus $^{132}$Sn and the corresponding isotopes and isotones. The two-neutron and two-proton separation energies are well reproduced by the RMF theory. In particular, the RMF results agree with the experimental single-particle spectrum in $^{132}$Sn as well as the Nilsson spin-orbit parameter $C$ and orbit-orbit parameter $D$ thus extracted, but remarkably differ from the traditional Nilsson parameters. Furthermore, the present results provide a guideline for the isospin dependence of the Nilsson parameters.

preprint2010arXiv

Tensor coupling effects on spin symmetry in anti-Lambda spectrum of hypernuclei

The effects of $\barΛ\barΛω$-tensor coupling on the spin symmetry of $\barΛ$ spectra in $\barΛ$-nucleus systems have been studied with the relativistic mean-field theory. Taking $^{12}$C+$\barΛ$ as an example, it is found that the tensor coupling enlarges the spin-orbit splittings of $\barΛ$ by an order of magnitude although its effects on the wave functions of $\barΛ$ are negligible. Similar conclusions has been observed in $\barΛ$-nucleus of different mass regions, including $^{16}$O+$\barΛ$, $^{40}$Ca+$\barΛ$ and $^{208}$Pb+$\barΛ$. It indicates that the spin symmetry in anti-lambda-nucleus systems is still good irrespective of the tensor coupling.

preprint2009arXiv

Avoid the Tsunami of the Dirac sea in the Imaginary Time Step method

The discrete single-particle spectra in both the Fermi and Dirac sea have been calculated by the imaginary time step (ITS) method for the Schrödinger-like equation after avoiding the "tsunami" of the Dirac sea, i.e., the diving behavior of the single-particle level into the Dirac sea in the direct application of the ITS method for the Dirac equation. It is found that by the transform from the Dirac equation to the Schrödinger-like equation, the single-particle spectra, which extend from the positive to the negative infinity, can be separately obtained by the ITS evolution in either the Fermi sea or the Dirac sea. Identical results with those in the conventional shooting method have been obtained via the ITS evolution for the equivalent Schrödinger-like equation, which demonstrates the feasibility, practicality and reliability of the present algorithm and dispels the doubts on the ITS method in the relativistic system.

preprint2009arXiv

Influence of nuclear physics inputs and astrophysical conditions on Th/U chronometer

The productions of thorium and uranium are key ingredients in $r$-process nucleo-cosmochronology. With the combination of improved nuclear and stellar data, we have made detailed investigations on the $r$-process abundance pattern in the very metal-poor halo stars based on the classical $r$-process approach. It is found that the results are almost independent of specified simulations to observed abundances. The influence from nuclear mass uncertainties on Th/U chronometer can approach 2 Gyr. Moreover, the ages of the metal-poor stars HE 1523-0901, CS 31082-001, and BD +17$^\circ$3248 are determined as $11.8\pm 3.7$, $13.5\pm 2.9$, and $10.9 \pm 2.9$ Gyr, respectively. The results can serve as an independent check for age estimate of the universe.

preprint2009arXiv

Neutron star properties in density-dependent relativistic Hartree-Fock theory

With the equations of state provided by the newly developed density dependent relativistic Hartree-Fock (DDRHF) theory for hadronic matter, the properties of the static and $β$-equilibrium neutron stars without hyperons are studied for the first time, and compared to the predictions of the relativistic mean field (RMF) models and recent observational data. The influences of Fock terms on properties of asymmetric nuclear matter at high densities are discussed in details. Because of the significant contributions from the $σ$- and $ω$-exchange terms to the symmetry energy, large proton fractions in neutron stars are predicted by the DDRHF calculations, which strongly affect the cooling process of the star. The critical mass about 1.45 $M_\odot$, close to the limit 1.5 $M_\odot$ determined by the modern soft X-ray data analysis, is obtained by DDRHF with the effective interactions PKO2 and PKO3 for the occurrence of direct Urca process in neutron stars. The maximum masses of neutron stars given by the DDRHF calculations lie between 2.45 M$_\odot$ and 2.49 M$_\odot$, which are in reasonable agreement with high pulsar mass $2.08 \pm 0.19 M_\odot$ from PSR B1516+02B. It is also found that the mass-radius relations of neutron stars determined by DDRHF are consistent with the observational data from thermal radiation measurement in the isolated neutron star RX J1856, QPOs frequency limits in LMXBs 4U 0614+09 and 4U 1636-536, and redshift determined in LMXBs EXO 0748-676.

preprint2009arXiv

Pairing Properties of Symmetric Nuclear Matter in Relativistic Mean Field Theory

The properties of pairing correlations in symmetric nuclear matter are studied in the relativistic mean field (RMF) theory with the effective interaction PK1. Considering well-known problem that the pairing gap at Fermi surface calculated with RMF effective interactions are three times larger than that with Gogny force, an effective factor in the particle-particle channel is introduced. For the RMF calculation with PK1, an effective factor 0.76 give a maximum pairing gap 3.2 MeV at Fermi momentum 0.9 fm$^{-1}$, which are consistent with the result with Gogny force.

preprint2009arXiv

Relativistic description of BCS-BEC crossover in nuclear matter

We study theoretically the di-neutron spatial correlations and the crossover from superfluidity of neutron Cooper pairs in the $^1S_0$ pairing channel to Bose-Einstein condensation (BEC) of di-neutron pairs for both symmetric and neutron matter in the microscopic relativistic pairing theory. We take the bare nucleon-nucleon interaction Bonn-B in the particle-particle channel and the effective interaction PK1 of the relativistic mean-field approach in the particle-hole channel. It is found that the spatial structure of neutron Cooper pair wave function evolves continuously from BCS-type to BEC-type as density decreases. We see a strong concentration of the probability density revealed for the neutron pairs in the fairly small relative distance around $1.5 {\rm fm}$ and the neutron Fermi momentum $k_{Fn}\in[0.6,1.0] {\rm fm^{-1}}$. However, from the effective chemical potential and the quasiparticle excitation spectrum, there is no evidence for the appearance of a true BEC state of neutron pairs at any density. The most BEC-like state may appear at $k_{Fn}\thicksim0.2 {\rm fm^{-1}}$ by examining the density correlation function. From the coherence length and the probability distribution of neutron Cooper pairs as well as the ratio between the neutron pairing gap and the kinetic energy at the Fermi surface, some features of the BCS-BEC crossover are seen in the density regions, $0.05 {\rm fm^{-1}}<k_{Fn}<0.7 (0.75) {\rm fm^{-1}}$, for the symmetric nuclear (pure neutron) matter.

preprint2009arXiv

Relativistic Hartree-Fock-Bogoliubov theory with Density Dependent Meson-Nucleon Couplings

Relativistic Hartree-Fock-Bogoliubov (RHFB) theory with density-dependent meson-nucleon couplings is presented. The integro-differential RHFB equations are solved by expanding the different components of the quasi-particle spinors in the complete set of eigen-solutions of the Dirac equations with Woods-Saxon potentials. Using the finite-range Gogny force D1S as an effective interaction in the pairing channel, systematic RHFB calculations are performed for Sn isotopes and N=82 isotones. It is demonstrated that an appropriate description of both mean field and pairing effects can be obtained within RHFB theory with finite range Gogny pairing forces. Better systematics are also found in the regions from the stable to the neutron-rich side with the inclusion of Fock terms, especially in the presence of $ρ$-tensor couplings.

preprint2009arXiv

Solving the Dirac equation with nonlocal potential by Imaginary Time Step method

The Imaginary Time Step (ITS) method is applied to solve the Dirac equation with the nonlocal potential in coordinate space by the ITS evolution for the corresponding Schrödinger-like equation for the upper component. It is demonstrated that the ITS evolution can be equivalently performed for the Schrödinger-like equation with or without localization. The latter algorithm is recommended in the application for the reason of simplicity and efficiency. The feasibility and reliability of this algorithm are also illustrated by taking the nucleus $^{16}$O as an example, where the same results as the shooting method for the Dirac equation with localized effective potentials are obtained.

preprint2009arXiv

Spin Symmetry for Anti-Lambda Spectrum in atomic nucleus

The spin symmetry of anti-Lambda spectrum in nucleus $^{16}$O has been studied in the relativistic mean field theory. The spin-orbit splittings of spin doublets are found to be around 0.03-0.07 MeV and the dominant components of the Dirac spinor for the anti-Lambda spin doublets are found to be near identical. It indicates that there is an even better spin symmetry in the anti-Lambda spectrum than that in the anti-nucleon spectrum.