Source author record

Bing-Nan Lu

Bing-Nan Lu 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

28works
9topics
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

28 published item(s)

preprint2026arXiv

Charge-dependent nucleon-nucleon interaction at N$^3$LO in nuclear lattice effective field theory

The nuclear lattice effective field theory (NLEFT) is an efficient tool for solving nuclear many-body problems, which takes high-fidelity lattice chiral interactions as input and computes nuclear low-energy observables via quantum Monte Carlo techniques. In this work, we present the first next-to-next-to-next-to-leading order (N$^3$LO) chiral forces on the lattice with the isospin-breaking effects fully taken into account. We focus on both the charge-independence breaking (CIB) and charge-symmetry breaking (CSB) effects. Specifically, we include the isospin-breaking effect from the mass difference between the charged and neutral pions in the one-pion-exchange potential (OPEP), the Coulomb force for the $pp$ interaction and the contribution of two additional charge-dependent contact operators. We also explicitly incorporate the two-pion-exchange potentials which was mostly neglected in previous NLEFT calculations. With these improvements, we are able to accurately reproduce the $np$ and $pp$ scattering phase shifts up to relative momentum $p \sim 200$ MeV as well as the deuteron properties. The construction of these charge-dependent lattice nuclear forces establishes a solid foundation for future high-precision nuclear ab initio calculations within the NLEFT framework.

preprint2026arXiv

Investigating nuclear beta decay using lattice quantum Monte Carlo approach

We present an \textit{ab initio} calculation of nuclear $β$ decay within the framework of nuclear lattice effective field theory (NLEFT), employing auxiliary-field quantum Monte Carlo methods to solve the nuclear many-body problem. Our approach combines next-to-next-to-leading order two- and three-body chiral interactions with one- and two-body axial current operators, all consistently derived in chiral effective field theory. Low-energy constants are determined exclusively from nucleon-nucleon scattering phase shifts and few-body observables for systems with $A \leq 3$. Using these interactions and transition operators, we perform two-channel Monte Carlo simulations to compute the $β$-decay matrix element for $^6$He, obtaining results in reasonable agreement with experimental measurements. To address the Monte Carlo sign problem, we implement a perturbative expansion around a leading-order Hamiltonian with approximate Wigner-SU(4) symmetry. This systematic approach provides a foundation for extending NLEFT simulations to precision studies of weak processes in medium-mass nuclei.

preprint2026arXiv

Sign-Problem-Free Nuclear Quantum Monte Carlo Simulation

Quantum Monte Carlo (QMC) methods offer exact solutions for quantum many-body systems but face severe limitations in fermionic systems like atomic nuclei due to the sign problem. While sign-problem-free QMC algorithms exist and provide valuable insights across disciplines, they have been restricted to simple models with limited quantitative predictive power. Here we overcome this barrier by developing a novel lattice nuclear force that is rigorously sign-problem-free for even-even nuclei. This interaction achieves a standard deviation of $σ= 2.932$ MeV from experimental binding energies for 76 even-even nuclei ($N,Z \leq 28$), matching state-of-the-art phenomenological mean-field models. Key innovations include the first sign-problem-free implementation of spin-orbit coupling for shell evolutions and an efficient QMC-optimized framework for global parameter fitting. Using this approach, we compute binding energies from $^4$He to $^{132}$Sn with unprecedented one-thousandth level numerical precision, reproduce symmetric nuclear matter saturation, and reveal novel spin-orbit-driven clustering in light nuclei. This work transforms sign-problem-free QMC into a scalable and predictive nuclear structure tool, while establishing a high-fidelity, non-perturbative foundation for \textit{ab initio} calculations of heavy nuclei.

preprint2022arXiv

Accelerated quantum Monte Carlo with mitigated error on noisy quantum computer

Quantum Monte Carlo and quantum simulation are both important tools for understanding quantum many-body systems. As a classical algorithm, quantum Monte Carlo suffers from the sign problem, preventing its application to most fermion systems and real time dynamics. In this paper, we introduce a novel non-variational algorithm using quantum simulation as a subroutine to accelerate quantum Monte Carlo by easing the sign problem. The quantum subroutine can be implemented with shallow circuits and, by incorporating error mitigation, can reduce the Monte Carlo variance by several orders of magnitude even when the circuit noise is significant. As such, the proposed quantum algorithm is applicable to near-term noisy quantum hardware.

preprint2022arXiv

Angular momentum and parity projected multidimensionally constrained relativistic Hartree-Bogoliubov model

The nuclear deformations are of fundamental importance in nuclear physics. Recently we developed a multi-dimensionally constrained relativistic Hartree-Bogoliubov (MDCRHB) model, in which all multipole deformations respecting the $V_4$ symmetry can be considered self-consistently. In this work we extend this model by incorporating the angular momentum projection (AMP) and parity projection (PP) to restore the rotational and parity symmetries broken in the mean-field level. This projected-MDCRHB (p-MDCRHB) model enables us to connect certain nuclear spectra to exotic intrinsic shapes such as triangle or tetrahedron. We present the details of the method and an exemplary calculation for $^{12}$C. We develop a triangular moment constraint to generate the triangular configurations consisting of three $α$ clusters arranged as an equilateral triangle. The resulting $^{12}$C spectra are consistent with that from a triangular rigid rotor for large separations between the $α$ clusters. We also calculate the $B(E2)$ and $B(E3)$ values for low-lying states and find good agreement with the experiments.

preprint2022arXiv

Perturbative quantum Monte Carlo method for nuclear physics

While first order perturbation theory is routinely used in quantum Monte Carlo (QMC) calculations, higher-order terms present significant numerical challenges. We present a new approach for computing perturbative corrections in projection QMC calculations. We demonstrate the method by computing nuclear ground state energies up to second order for a realistic chiral interaction. We calculate the binding energies of several light nuclei up to $^{16}$O by expanding the Hamiltonian around the Wigner SU(4) limit and find good agreement with data. In contrast to the natural ordering of the perturbative series, we find remarkably large second order energy corrections. This occurs because the perturbing interactions break the symmetries of the unperturbed Hamiltonian. Our method is free from the sign problem and can be applied to QMC calculations for many-body systems in nuclear physics, condensed matter physics, ultracold atoms, and quantum chemistry.

preprint2020arXiv

Projected Cooling Algorithm for Quantum Computation

In the current era of noisy quantum devices, there is a need for quantum algorithms that are efficient and robust against noise. Towards this end, we introduce the projected cooling algorithm for quantum computation. The projected cooling algorithm is able to construct the localized ground state of any Hamiltonian with a translationally-invariant kinetic energy and interactions that vanish at large distances. The term "localized" refers to localization in position space. The method can be viewed as the quantum analog of evaporative cooling. We start with an initial state with support over a compact region of a large volume. We then drive the excited quantum states to disperse and measure the remaining portion of the wave function left behind. For the nontrivial examples we consider here, the improvement over other methods is substantial. The only additional resource required is performing the operations in a volume significantly larger than the size of the localized state. These characteristics make the projected cooling algorithm a promising tool for calculations of self-bound systems such as atomic nuclei.

preprint2020arXiv

Superfluid Condensate Fraction and Pairing Wave Function of the Unitary Fermi Gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gases and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this work we present auxiliary-field lattice Monte Carlo simulations using a novel lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test we compute the ground state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{FG}$, we find $E_0/E_{FG}= 0.369(2), 0.372(2)$, using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is $α= 0.43(2)$. We also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F = 1.8(3) \hbar$, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic P-wave and D-wave superfluids.

preprint2016arXiv

Multidimensionally-constrained relativistic mean-field study of spontaneous fission: coupling between shape and pairing degrees of freedom

Studies of fission dynamics, based on nuclear energy density functionals, have shown that the coupling between shape and pairing degrees of freedom has a pronounced effect on the nonperturbative collective inertia and, therefore, on dynamic (least-action) spontaneous fission paths and half-lives. Collective potentials and nonperturbative cranking collective inertia tensors are calculated using the multidimensionally-constrained relativistic mean-field (MDC-RMF) model. Pairing correlations are treated in the BCS approximation using a separable pairing force of finite range. Pairing fluctuations are included as a collective variable using a constraint on particle-number dispersion. Fission paths are determined with the dynamic programming method by minimizing the action in multidimensional collective spaces. The dynamics of spontaneous fission of $^{264}$Fm and $^{250}$Fm are explored. Fission paths, action integrals and corresponding half-lives computed in the three-dimensional collective space of shape and pairing coordinates, using the relativistic functional DD-PC1 and a separable pairing force of finite range, are compared with results obtained without pairing fluctuations. Results for $^{264}$Fm are also discussed in relation with those recently obtained using the HFB model. The inclusion of pairing correlations in the space of collective coordinates favors axially symmetric shapes along the dynamic path of the fissioning system, amplifies pairing as the path traverses the fission barriers, significantly reduces the action integral and shortens the corresponding SF half-life.

preprint2016arXiv

Nuclear binding near a quantum phase transition

How do protons and neutrons bind to form nuclei? This is the central question of ab initio nuclear structure theory. While the answer may seem as simple as the fact that nuclear forces are attractive, the full story is more complex and interesting. In this work we present numerical evidence from ab initio lattice simulations showing that nature is near a quantum phase transition, a zero-temperature transition driven by quantum fluctuations. Using lattice effective field theory, we perform Monte Carlo simulations for systems with up to twenty nucleons. For even and equal numbers of protons and neutrons, we discover a first-order transition at zero temperature from a Bose-condensed gas of alpha particles (4He nuclei) to a nuclear liquid. Whether one has an alpha-particle gas or nuclear liquid is determined by the strength of the alpha-alpha interactions, and we show that the alpha-alpha interactions depend on the strength and locality of the nucleon-nucleon interactions. This insight should be useful in improving calculations of nuclear structure and important astrophysical reactions involving alpha capture on nuclei. Our findings also provide a tool to probe the structure of alpha cluster states such as the Hoyle state responsible for the production of carbon in red giant stars and point to a connection between nuclear states and the universal physics of bosons at large scattering length.

preprint2016arXiv

Precise determination of lattice phase shifts and mixing angles

We introduce a general and accurate method for determining lattice phase shifts and mixing angles, which is applicable to arbitrary, non-cubic lattices. Our method combines angular momentum projection, spherical wall boundaries and an adjustable auxiliary potential. This allows us to construct radial lattice wave functions and to determine phase shifts at arbitrary energies. For coupled partial waves, we use a complex-valued auxiliary potential that breaks time-reversal invariance. We benchmark our method using a system of two spin-1/2 particles interacting through a finite-range potential with a strong tensor component. We are able to extract phase shifts and mixing angles for all angular momenta and energies, with precision greater than that of extant methods. We discuss a wide range of applications from nuclear lattice simulations to optical lattice experiments.

preprint2015arXiv

Breaking and restoration of rotational symmetry for irreducible tensor operators on the lattice

We study the breaking of rotational symmetry on the lattice for irreducible tensor operators and practical methods for suppressing this breaking. We illustrate the features of the general problem using an $α$ cluster model for $^{8}$Be. We focus on the lowest states with non-zero angular momentum and examine the matrix elements of multipole moment operators. We show that the physical reduced matrix element is well reproduced by averaging over all possible orientations of the quantum state, and this is expressed as a sum of matrix elements weighted by the corresponding Clebsch-Gordan coefficients. For our $α$ cluster model we find that the effects of rotational symmetry breaking can be largely eliminated for lattice spacings of $a\leq 1.7$ fm, and we expect similar improvement for actual lattice Monte Carlo calculations.

preprint2015arXiv

Multidimensionally-constrained relativistic Hartree-Bogoliubov study of nuclear spontaneous fission

Recent microscopic studies, based on the theoretical framework of nuclear energy density functionals, have analyzed dynamic (least action) and static (minimum energy) fission paths, and it has been shown that in addition to the important role played by nonaxial and/or octupole collective degrees of freedom, fission paths crucially depend on the approximations adopted in calculating the collective inertia. The dynamics of spontaneous fission of $^{264}$Fm and $^{250}$Fm is explored. The fission paths, action integrals and the corresponding half-lives predicted by the functionals PC-PK1 and DD-PC1 are compared and, in the case of $^{264}$Fm, discussed in relation with recent results obtained using the HFB model based on the Skyrme functional SkM$^*$ and a density dependent mixed pairing interaction. Deformation energy surfaces, collective potentials, and perturbative and nonperturbative cranking collective inertia tensors are calculated using the multidimensionally-constrained relativistic Hartree-Bogoliubov (MDC-RHB) model, with the energy density functionals PC-PK1 and DD-PC1. Pairing correlations are treated in the Bogoliubov approximation using a separable pairing force of finite range. The least-action principle is employed to determine dynamic spontaneous fission paths.

preprint2015arXiv

Multidimensionally-constrained relativistic mean-field study of triple-humped barriers in actinides

Potential energy surfaces (PES's) of actinide nuclei are characterized by a two-humped barrier structure. At large deformations beyond the second barrier the occurrence of a third one was predicted by Mic-Mac model calculations in the 1970s, but contradictory results were later reported. In this paper, triple-humped barriers in actinide nuclei are investigated with covariant density functional theory (CDFT). Calculations are performed using the multidimensionally-constrained relativistic mean field (MDC-RMF) model, with functionals PC-PK1 and DD-ME2. Pairing correlations are treated in the BCS approximation with a separable pairing force of finite range. Two-dimensional PES's of $^{226,228,230,232}$Th and $^{232,234,236,238}$U are mapped and the third minima on these surfaces are located. Then one-dimensional potential energy curves along the fission path are analyzed in detail and the energies of the second barrier, the third minimum, and the third barrier are determined. DD-ME2 predicts the occurrence of a third barrier in all Th nuclei and $^{238}$U. The third minima in $^{230,232}$Th are very shallow, whereas those in $^{226,228}$Th and $^{238}$U are quite prominent. With PC-PK1 a third barrier is found only in $^{226,228,230}$Th. Single-nucleon levels around the Fermi surface are analyzed in $^{226}$Th, and it is found that the formation of the third minimum is mainly due to the $Z=90$ proton energy gap at $β_{20} \approx 1.5$ and $β_{30} \approx 0.7$. The possible occurrence of a third barrier in actinide nuclei depends on the effective interaction used in multidimensional CDFT calculations. More pronounced minima are predicted by the DD-ME2 functional, as compared to the functional PC-PK1. The depth of the third well in Th isotopes decreases with increasing neutron number. The origin of the third minimum is due to the proton $Z=90$ shell gap at relevant deformations.

preprint2014arXiv

Breaking and restoration of rotational symmetry on the lattice for bound state multiplets

We explore the breaking of rotational symmetry on the lattice for bound state energies and practical methods for suppressing this breaking. We demonstrate the general problems associated with lattice discretization errors and finite-volume errors using an $α$ cluster model for $^8$Be and $^{12}$C. We consider the two and three $α$-particle systems and focus on the lowest states with non-zero angular momentum which split into multiplets corresponding to different irreducible representations of the cubic group. We examine the dependence of such splittings on the lattice spacing and box size. We find that lattice spacing errors are closely related to the commensurability of the lattice with the intrinsic length scales of the system. We also show that rotational symmetry breaking effects can be significantly reduced by using improved lattice actions, and that the physical energy levels are accurately reproduced by the weighted average of a given spin multiplets.

preprint2014arXiv

Microscopic description of octupole shape-phase transitions in light actinides and rare-earth nuclei

A systematic analysis of low-lying quadrupole and octupole collective states is presented, based on the microscopic energy density functional framework. By mapping the deformation constrained self-consistent axially symmetric mean-field energy surfaces onto the equivalent Hamiltonian of the $sdf$ interacting boson model (IBM), that is, onto the energy expectation value in the boson condensate state, the Hamiltonian parameters are determined. The study is based on the global relativistic energy density functional DD-PC1. The resulting IBM Hamiltonian is used to calculate excitation spectra and transition rates for the positive- and negative-parity collective states in four isotopic chains characteristic for two regions of octupole deformation and collectivity: Th, Ra, Sm and Ba. Consistent with the empirical trend, the microscopic calculation based on the systematics of $β_{2}$-$β_{3}$ energy maps, the resulting low-lying negative-parity bands and transition rates show evidence of a shape transition between stable octupole deformation and octupole vibrations characteristic for $β_{3}$-soft potentials.

preprint2014arXiv

Multidimensionally-constrained relativistic mean field models and potential energy surfaces of actinide nuclei

By breaking both the axial and the spatial reflection symmetries, we develop multidimensionally constrained relativistic mean field (MDC-RMF) models. The nuclear shape is assumed to be invariant under the reversion of $x$ and $y$ axes, i.e., the intrinsic symmetry group is $V_{4}$ and all shape degrees of freedom $β_{λμ}$ with even $μ$, such as $β_{20}$, $β_{22}$, $β_{30}$, $β_{32}$, $β_{40}$, $\dots$, are included self-consistently. The Dirac wave functions are expanded in an axially deformed harmonic oscillator basis. The RMF functional can be one of the following four forms: the meson exchange or point-coupling nucleon interactions combined with the nonlinear or density-dependent couplings. The pairing effects are taken into account with the BCS approach. The 1-, 2-, and 3-D potential energy surfaces (PES's) of $^{240}$Pu are illustrated for numerical checks and for the study of the effect of the triaxiality on the fission barriers. Potential energy curves of even-even actinide nuclei around the first and second fission barriers are studied systematically. Besides the first ones, the second fission barriers in these nuclei are also lowered considerably by the triaxial deformation. This lowering effect is independent of the effective interactions used in the RMF functionals. Further discussions are made about different predictions on the effect of the triaxiality between the macroscopic-microscopic and MDC-RMF models, possible discontinuities on PES's from self-consistent approaches, and the restoration of broken symmetries. MDC-RMF models give reasonably good description of fission barriers of even-even actinide nuclei. It is important to include both the nonaxial and the reflection asymmetric shapes simultaneously for the study of PES's and fission barriers of actinide nuclei and of those in unknown mass regions such as, e.g., superheavy nuclei.

preprint2014arXiv

Superdeformed $Λ$ hypernuclei from relativistic mean field models

We study the superdeformed (SD) states and corresponding SD hypernuclei of Ar isotopes with the multidimensionally-constrained relativistic mean field (MDC-RMF) models which can accommodate various shape degree of freedom. %The hyperon-nucleon $σ$-$Λ$ and $ω$-$Λ$ interactions are introduced %in the same covariant framework. We found that the density profiles of SD states in Ar isotopes show a strong localization with a ring structure near the surface, while the central part of the density is dilute showing a hole structure. This localization of SD density induces an appreciable deformation in the hyperon wave function and results in a large overlap between the core and the hyperon in the SD hypernuclei of Ar isotopes. Then the $Λ$ separation energy of SD state becomes larger than that of normally deformed or spherical ground state. This feature is different from that found in other nuclei such as $^{32}$S, $^{56}$Ni, and $^{60}$Zn in which the $Λ$ separation energy of larger deformed state is smaller. In this context, the measurement of the $Λ$ separation energy may provide an important information on the localization of the density profile of SD states.

preprint2013arXiv

Exact conservation and breaking of pseudospin symmetry in single particle resonant states

In this contribution we present some results on the study of pseudospin symmetry (PSS) in single particle resonant states. The PSS is a relativistic dynamical symmetry connected with the small component of the nucleon Dirac wave function. Many efforts have been made to study this symmetry in bound states. We recently gave a rigorous justification of the PSS in single particle resonant states by examining the zeros of Jost functions corresponding to the small components of the radial Dirac wave functions and phase shifts of continuum states [1, Phys. Rev. Lett. 109 (2012) 072501]. We have shown that the PSS in single particle resonant states in nuclei is conserved when the attractive scalar and repulsive vector potentials have the same magnitude but opposite sign. Examples of exact conservation and breaking of this symmetry in single particle resonances are given for spherical square-well and Woods-Saxon potentials.

preprint2013arXiv

Multi-dimensional constraint relativistic mean field model and applications in actinide and transfermium nuclei

In this contribution we present some results of potential energy surfaces of actinide and transfermium nuclei from multi-dimensional constrained relativistic mean field (MDC-RMF) models. Recently we developed multi-dimensional constrained covariant density functional theories (MDC-CDFT) in which all shape degrees of freedom $β_{λμ}$ with even $μ$ are allowed and the functional can be one of the following four forms: the meson exchange or point-coupling nucleon interactions combined with the non-linear or density-dependent couplings. In MDC-RMF models, the pairing correlations are treated with the BCS method. With MDC-RMF models, the potential energy surfaces of even-even actinide nuclei were investigated and the effect of triaxiality on the fission barriers in these nuclei was discussed. The non-axial reflection-asymmetric $β_{32}$ shape in some transfermium nuclei with $N=150$, namely $^{246}$Cm, $^{248}$Cf, $^{250}$Fm, and $^{252}$No were also studied.

preprint2013arXiv

Multi-dimensional potential energy surfaces and non-axial octupole correlations in actinide and transfermium nuclei from relativistic mean field models

We have developed multi-dimensional constrained covariant density functional theories (MDC-CDFT) for finite nuclei in which the shape degrees of freedom β_{λμ} with even μ, e.g., β_{20}, β_{22}, β_{30}, β_{32}, β_{40}, etc., can be described simultaneously. The functional can be one of the following four forms: the meson exchange or point-coupling nucleon interactions combined with the non-linear or density-dependent couplings. For the pp channel, either the BCS approach or the Bogoliubov transformation is implemented. The MDC-CDFTs with the BCS approach for the pairing (in the following labelled as MDC-RMF models with RMF standing for "relativistic mean field") have been applied to investigate multi-dimensional potential energy surfaces and the non-axial octupole $Y_{32}$-correlations in N=150 isotones. In this contribution we present briefly the formalism of MDC-RMF models and some results from these models. The potential energy surfaces with and without triaxial deformations are compared and it is found that the triaxiality plays an important role upon the second fission barriers of actinide nuclei. In the study of Y_{32}-correlations in N=150 isotones, it is found that, for 248Cf and 250Fm, β_{32} > 0.03 and the energy is lowered by the β_{32} distortion by more than 300 keV; while for 246Cm and 252No, the pocket with respect to β_{32} is quite shallow.

preprint2013arXiv

Potential energy surfaces of actinide and transfermium nuclei from multi-dimensional constraint covariant density functional theories

Multi-dimensional constrained covariant density functional theories were developed recently. In these theories, all shape degrees of freedom β_{λμ} deformations with even μare allowed, e.g., β_{20}, β_{22}, β_{30}, β_{32}, β_{40}, β_{42}, β_{44}, and so on and the CDFT functional can be one of the following four forms: the meson exchange or point-coupling nucleon interactions combined with the non-linear or density-dependent couplings. In this contribution, some applications of these theories are presented. The potential energy surfaces of actinide nuclei in the (β_{20}, β_{22}, β_{30}) deformation space are investigated. It is found that besides the octupole deformation, the triaxiality also plays an important role upon the second fission barriers. The non-axial reflection-asymmetric β_{32} shape in some transfermium nuclei with N = 150, namely 246Cm, 248Cf, 250Fm, and 252No are studied.

preprint2013arXiv

Pseudospin symmetry in single particle resonances in spherical square wells

The pseudospin symmetry (PSS) has been studied extensively for bound states. Recently we justified rigorously that the PSS in single particle resonant states is exactly conserved when the attractive scalar and repulsive vector potentials of the Dirac Hamiltonian have the same magnitude but opposite sign [PRL 109, 072501 (2012)]. To understand more deeply the PSS, we focus on several issues related to the exact conservation and breaking mechanism of the PSS in single particle resonances. In particular, we are interested in how the energy and width splittings of PS partners depend on the depth of the scalar and vector potentials. We investigate the asymptotic behaviors of radial Dirac wave functions. Spherical square well potentials are employed in which the PSS breaking part in the Jost function can be well isolated. By examining the zeros of Jost functions corresponding to small components of the radial Dirac wave functions, general properties of the PSS are analyzed. By examining the Jost function, the occurrence of intruder orbitals is explained and it is possible to trace continuously the PSS partners from the PSS limit to the case with a finite potential depth. The dependence of the PSS in resonances as well as in bound states on the potential depth is investigated systematically. We find a threshold effect in the energy splitting and an anomaly in the width splitting of pseudospin partners when the depth of the single particle potential varies from zero to a finite value. The conservation and the breaking of the PSS in resonant states and bound states share some similar properties. The appearance of intruder states can be explained by examining the zeros of Jost functions. Origins of the threshold effect in the energy splitting and the anomaly in the width splitting of PS partners, together with many other problems, are still open and should be further investigated.

preprint2013arXiv

Recent progresses on the pseudospin symmetry in single particle resonant states

The pseudospin symmetry (PSS) is a relativistic dynamical symmetry directly connected with the small component of the nucleon Dirac wave function. Much effort has been made to study this symmetry in bound states. Recently, a rigorous justification of the PSS in single particle resonant states was achieved by examining the asymptotic behaviors of the radial Dirac wave functions: The PSS in single particle resonant states in nuclei is conserved exactly when the attractive scalar and repulsive vector potentials have the same magnitude but opposite sign. Several issues related to the exact conservation and breaking mechanism of the PSS in single particle resonances were investigated by employing spherical square well potentials in which the PSS breaking part can be well isolated in the Jost function. A threshold effect in the energy splitting and an anomaly in the width splitting of pseudospin partners were found when the depth of the square well potential varies from zero to a finite value.

preprint2012arXiv

Nonaxial-octupole Y_{32} correlations in N = 150 isotones from multidimensional constrained covariant density functional theories

The non-axial reflection-asymmetric $β_{32}$ shape in some transfermium nuclei with N=150, namely $^{246}$Cm, $^{248}$Cf, $^{250}$Fm, and $^{252}$No are investigated with multidimensional constrained covariant density functional theories. By using the density-dependent point coupling covariant density functional theory with the parameter set DD-PC1 in the particle-hole channel, it is found that, for the ground states of $^{248}$Cf and $^{250}$Fm, the non-axial octupole deformation parameter $β_{32} > 0.03$ and the energy gain due to the $β_{32}$ distortion is larger than 300 keV. In $^{246}$Cm and $^{252}$No, shallow $β_{32}$ minima are found. The occurrence of the non-axial octupole $β_{32}$ correlations is mainly from a pair of neutron orbitals $[734]9/2$ ($νj_{15/2}$) and $[622]5/2$ ($νg_{9/2}$) which are close to the neutron Fermi surface and a pair of proton orbitals $[521]3/2$ ($πf_{7/2}$) and $[633]7/2$ ($πi_{13/2}$) which are close to the proton Fermi surface. The dependence of the non-axial octupole effects on the form of energy density functional and on the parameter set is also studied.

preprint2012arXiv

Potential Energy Surfaces of Actinide Nuclei from a Multi-dimensional Constraint Covariant Density Functional Theory: Barrier Heights and Saddle Point Shapes

For the first time the potential energy surfaces of actinide nuclei in the $(β_{20}, β_{22}, β_{30})$ deformation space are obtained from a multi-dimensional constrained covariant density functional theory. With this newly developed theory we are able to explore the importance of the triaxial and octupole shapes simultaneously along the whole fission path. It is found that besides the octupole deformation, the triaxiality also plays an important role upon the second fission barriers. The outer barrier as well as the inner barrier are lowered by the triaxial deformation compared with axially symmetric results. This lowering effect for the reflection asymmetric outer barrier is 0.5 $\sim$ 1 MeV, accounting for $10 \sim 20%$ of the barrier height. With the inclusion of the triaxial deformation, a good agreement with the data for the outer barriers of actinide nuclei is achieved.

preprint2012arXiv

Pseudospin symmetry in single particle resonant states

The pseudospin symmetry is a relativistic dynamical symmetry connected with the small component of the Dirac spinor. The origin of pseudospin symmetry in single particle bound states in atomic nuclei has been revealed and studied extensively. By examining the zeros of Jost functions corresponding to the small components of Dirac wave functions and phase shifts of continuum states, we show that the pseudospin symmetry in single particle resonant states in nuclei is conserved when the attractive scalar and repulsive vector potentials have the same magnitude but opposite sign. The exact conservation and the breaking of pseudospin symmetry are illustrated for single particle resonances in spherical square-well and Woods-Saxon potentials.

preprint2011arXiv

Quadrupole deformation $(β,γ)$ of light $Λ$ hypernuclei in constrained relativistic mean field model: shape evolution and shape polarization effect of $Λ$ hyperon

The shapes of light normal nuclei and $Λ$ hypernuclei are investigated in the $(β, γ)$ deformation plane by using a newly developed constrained relativistic mean field (RMF) model. As examples, the results of some C, Mg, and Si nuclei are presented and discussed in details. We found that for normal nuclei the present RMF calculations and previous Skyrme-Hartree-Fock models predict similar trends of the shape evolution with the neutron number increasing. But some quantitative aspects from these two approaches, such as the depth of the minimum and the softness in the $γ$ direction, differ a lot for several nuclei. For $Λ$ hypernuclei, in most cases, the addition of a $Λ$ hyperon alters slightly the location of the ground state minimum towards the direction of smaller $β$ and softer $γ$ in the potential energy surface $E \sim (β, γ)$. There are three exceptions, namely, $^{13}_Λ$C, $^{23}_Λ$C, and $^{31}_Λ$Si in which the polarization effect of the additional $Λ$ is so strong that the shapes of these three hypernuclei are drastically different from their corresponding core nuclei.