2024-11-06 | | Total: 6

Shell quenching phenomena in nuclear charge radii are typically observed at the well-established neutron magic numbers. However, the recent discovery of potential new magic numbers at the neutron numbers $N = 32$ and $N = 34$ has sparked renewed interest in this mass region. This work further inspects into the charge radii of nuclei around the $N = 28$ shell closure using the relativistic Hartree-Bogoliubov model. We incorporate meson-exchange and point-coupling effective nucleon-nucleon interactions alongside the Bogoliubov transformation for pairing corrections. To accurately capture the odd-even staggering and shell closure effects observed in charge radii, neutron-proton correlations around Fermi surface are explicitly considered. The charge radii of Ca and Ni isotopes are used to test the theoretical model and show an improvement with neutron-proton pairing corrections, in particular for neutron-rich isotopes.Our calculations reveal a inverted parabolic-like trend in the charge radii along the $N = 28$ isotones for proton numbers $Z$ between 20 and 28. Additionally, the shell closure effect of $Z = 28$ persists across the $N = 28$, $30$, $32$, and $34$ isotonic chains, albeit with a gradual weakening trend. Notably, the significantly abrupt changes in charge radii are observed across $Z = 22$ along both the $N = 32$ and $N = 34$ isotonic chains. This kink at $Z = 22$ comes from the sudden decrease of the neuron-proton correlation around Fermi surfaces across $Z = 22$ for $N = 30$, 32, and 34 isotones, and might provide a signature for identifying the emergence of neutron magic numbers $N=32$ and 34.Furthermore, the calculated charge radii for these isotonic chains ($N = 28$, 30, 32, and 34) can serve as reliable guidelines for future experimental measurements.

This study investigates the tetrahedral structure in $^{80}$Zr and Lambda ($\Lambda$) impurity effect in $^{81}_{~\Lambda}$Zr using the multidimensionally constrained relativistic Hartree-Bogoliubov model. The ground states of both $^{80}$Zr and $^{81}_{~\Lambda}$Zr exhibit a tetrahedral configuration, accompanied by prolate and axial-octupole shape isomers. Our calculations reveal there are changes in the deformation parameters $\beta_{20}$, $\beta_{30}$, and $\beta_{32}$ upon $\Lambda$ binding to $^{80}$Zr, except for $\beta_{32}$ when $\Lambda$ occupies $p$-orbits. Compared to the two shape isomers, the $\Lambda$ particle exhibits weaker binding energy in the tetrahedral state when occupying the $1/2^+[000](\Lambda_s)$ or $1/2^-[110]$ single-particle states. In contrast, the strongest binding occurs for the $\Lambda$ particle in the $1/2^-[101]$ state with tetrahedral shape. Besides, a large $\Lambda$ separation energy may not necessarily correlate with a significant overlap between the density distributions of the $\Lambda$ particle and the nuclear core, particularly for tetrahedral hypernuclei.

The equation of state (EOS) of extremely dense matter is crucial for understanding the properties of rotating neutron stars. Starting from the widely used realistic Bonn potentials rooted in a relativistic framework, we derive EOSs by performing the state-of-the-art relativistic Brueckner-Hartree-Fock (RBHF) calculations in the full Dirac space. The self-consistent and simultaneous consideration of both positive- and negative-energy states (NESs) of the Dirac equation allows us to avoid the uncertainties present in calculations where NESs are treated using approximations. To manifest the impact of rotational dynamics, several structural properties of neutron stars across a wide range of rotation frequencies and up to the Keplerian limit are obtained, including the gravitational and baryonic masses, the polar and equatorial radii, and the moments of inertia. Our theoretical predictions align well with the latest astrophysical constraints from the observations on massive neutron stars and joint mass-radius measurements. The maximum mass for rotating configurations can reach up to $2.93M_{\odot}$ for Bonn A potential, while the radius of a $1.4M_\odot$ neutron star for non-rotating case can be extended to around 17 km through the constant baryonic mass sequences. Relations with good universalities between the Keplerian frequency and static mass as well as radius are obtained, from which the radius of the black widow PSR J0952-0607 is predicted to be less than 19.58 km. Furthermore, to understand how rotation deforms the equilibrium shape of a neutron star, the eccentricity is also calculated. The approximate universality between the eccentricity at the Keplerian frequency and the gravitational mass is found.

The subtraction function plays a pivotal role in calculations involving the forward Compton amplitude, which is crucial for predicting the Lamb shift in muonic atom, as well as the proton-neutron mass difference. In this work, we present a lattice QCD calculation of the subtraction function using two domain wall fermion gauge ensembles at the physical pion mass. We utilize a recently proposed subtraction point, demonstrating its advantage in mitigating statistical and systematic uncertainties by eliminating the need for ground-state subtraction. Our results reveal significant contributions from $N\pi$ intermediate states to the subtraction function. Incorporating these contributions, we compute the proton, neutron and nucleon isovector subtraction functions at photon momentum transfer $Q^2\in[0,2]$ GeV$^2$. For the proton subtraction function, we compare our lattice results with chiral perturbation theory prediction at low $Q^2$ and with the results from the perturbative operator-product expansion at high $Q^2$. Finally, using these subtraction functions as input, we determine their contribution to two-photon exchange effects in the Lamb shift and isovector nucleon electromagnetic self-energy.

We propose a mechanism which explains the masses of $\eta$ and $\eta'$ mesons without invoking the explicit violation of $U(1)_A$ symmetry by the chiral anomaly. It is shown that the U(1) problem, the problem for which the prediction of $\eta$ and $\eta'$ masses in the simple chiral perturbation theory largely deviates from the experimental values, is actually resolved by considering the first order contribution of the disconnected meson correlator with respect to the quark mass. The bound of Weinberg $m_\eta^2 \le 3 m_\pi^2$ is fulfilled by considering the negative squared mass of $\eta$ or $\eta'$ which is just the saddle point of the QCD effective potential, and 20% level agreements with experimental data are obtained by just fitting one low energy constant. We provide the leading chiral Lagrangian due to the disconnected contribution in 3-flavor QCD, and also discuss the 2- and 4-flavor cases as well as the consistency of our mechanism with the chiral restoration at high temperature found in lattice calculations.

Simulations of collisions of fundamental particles on a quantum computer are expected to have an exponential advantage over classical methods and promise to enhance searches for new physics. Furthermore, scattering in scalar field theory has been shown to be BQP-complete, making it a representative problem for which quantum computation is efficient. As a step toward large-scale quantum simulations of collision processes, scattering of wavepackets in one-dimensional scalar field theory is simulated using 120 qubits of IBM's Heron superconducting quantum computer ibm_fez. Variational circuits compressing vacuum preparation, wavepacket initialization, and time evolution are determined using classical resources. By leveraging physical properties of states in the theory, such as symmetries and locality, the variational quantum algorithm constructs scalable circuits that can be used to simulate arbitrarily-large system sizes. A new strategy is introduced to mitigate errors in quantum simulations, which enables the extraction of meaningful results from circuits with up to 4924 two-qubit gates and two-qubit gate depths of 103. The effect of interactions is clearly seen, and is found to be in agreement with classical Matrix Product State simulations. The developments that will be necessary to simulate high-energy inelastic collisions on a quantum computer are discussed.