
We propose a method to extract the properties of the isobaric mass parabola
based on the total double $\beta$ decay energies of isobaric nuclei. Two
important parameters of the mass parabola, the location of the most
$\beta$stable nuclei $Z_{A}$ and the curvature parameter $b_{A}$, are obtained
for 251 A values based on the total double $\beta$ decay energies of nuclei
compiled in AME2016 database.The advantage of this approach is that we can
remove the pairing energy term $P_{A}$ caused by oddeven variation, and the
mass excess $M(A,Z_{A})$ of the most stable nuclide for mass number $A$ in the
performance process, which are used in the mass parabolic fitting method. The
Coulomb energy coefficient $a_{c}=0.6910$ MeV is determined by the mass
difference relation of mirror nuclei, and the symmetry energy coefficient is
also studied by the relation $a_{\rm sym}(A)=0.25b_{A}Z_{A}$.

Recent highprecision measurements of nuclear deep inelastic scattering at
high x and moderate 6 < Q$^2$ < 9GeV$^2$ give a rare opportunity to reach the
quark distributions in the {\it superfast} region, in which the momentum
fraction of the nucleon carried by its constituent quark is larger than the
total fraction of the nucleon at rest, x>1. We derive the leadingorder QCD
evolution equation for such quarks with the goal of relating the moderateQ$^2$
data to the two earlier measurements of superfast quark distributions at large
60 < Q$^2$ < 200~GeV$^2$. Since the highQ$^2$ measurements gave strongly
contradictory estimates of the nuclear effects that generate superfast quarks,
relating them to the highprecision, moderateQ$^2$ data through QCD evolution
allows us to clarify this longstanding issue. Our calculations indicate that
the moderateQ$^2$ data at $x\lesssim 1.05$ are in better agreement with the
highQ$^2$ data measured in (anti)neutrinonuclear reactions which require
substantial highmomentum nuclear effects in the generation of superfast
quarks. Our prediction for the highQ$^2$ and x>1.1 region is somewhat in the
middle of the neutrinonuclear and muonnuclear scattering data.

The nuclear optical model potential (OMP) is generally assumed to be
independent of the orbital angular momentum, $l$, of the interacting nuclei.
Nucleonnucleus and nucleusnucleus interactions are customarily $l$
independent in calculations of nuclear elastic scattering and in standard
reaction codes. The evidence for various forms of $l$ dependence of OMPs is
reviewed and the importance of implementing these forms is evaluated. Existing
arguments and evidence for $l$ dependence are reviewed and new arguments and
calculations are introduced. The relationship is examined between (i) $l$
dependence, and, (ii) the undularity (waviness) of $l$independent potentials
that are $S$matrix equivalent to $l$dependent potentials. Such undularity is
a property of the dynamic polarisation potential (DPP) generated by the
coupling to reaction channels, or by coupling to excited states of the target
or projectile nuclei. Various examples, particularly involving weakly bound
projectile nuclei, are presented. Such undularity also occurs in
$l$independent potentials that have been found in model independent fits to
precise wide angular range elastic scattering angular distributions. Cases of
such phenomenological undularity, for both light and heavy ions, are referenced
and shown to be related to undulatory properties of the dynamic polarisation
potentials (DPPs) arising from channel coupling. Other forms of $l$ dependence,
that could be standard options in direct reaction codes, are noted. Of
particular importance are $l$ dependencies arising from full
antisymmetrization. The case is made that reactioninduced $l$ dependence is a
general property of nucleonnucleus and nucleusnucleus interactions and
represents a valid extension of the nuclear optical model.

Within a Boltzmann transport model, we demonstrate correlation between
stopping observables and shear viscosity in central nuclear collisions at
intermediate energies (on the order of 10 to 1000 MeV/nucleon). The correlation
allows us to assess the viscosity of nuclear matter, by tuning the inmedium
nucleonnucleon cross section in our transport model to agree with nuclear
stopping data. We also calculate the ratio of shear viscosity to entropy
density to determine how close the system is to the universal quantum lower
limit proposed in the context of ultrarelativistic heavy ion collisions.

The distribution of flow harmonics in heavy ion experiment can be
characterized by standardized cumulants. We first model the ellipticity and
power parameters of the ellipticpower distribution by employing MCGlauber
model. Then we use the ellipticpower distribution together with the
hydrodynamic linear response approximation to study the two dimensional
standardized cumulants of elliptic and triangular flow ($v_2$ and $v_3$)
distribution. For the second harmonic, it turns out that finding two
dimensional cumulants in terms of $2q$particle correlation functions
$c_2\{2q\}$ is limited to the skewness. We also show that $c_3\{2\}$,
$c_3\{4\}$, and $c_3\{6\}$, are related to the second, fourth, and sixth
standardized cumulants of the $v_3$ distribution, respectively. The cumulant
$c_{n}\{2q\}$ can be also written in terms of $v_n\{2q\}$. Specifically,
$(v_3\{4\}/v_3\{2\})^4$ turns out to be the kurtosis of the $v_3$
eventbyevent fluctuation distribution. We introduce a new parametrization for
the distribution $p(v_3)$ with $v_3\{2\}$, kurtosis and sixthorder
standardized cumulant being its free parameters. Compared to the Gaussian
distribution, it indicates a more accurate fit with experimental results.
Finally, we compare the kurtosis obtained from simulation with that of
extracted from experimental data for the $v_3$ distribution.

Topological gluon configurations in quantum chromodynamics induce quark
chirality imbalance in local domains, which can result in the chiral magnetic
effect (CME)an electric charge separation along a strong magnetic field.
Experimental searches for the CME in relativistic heavy ion collisions via the
chargedependent azimuthal correlator ($\Delta\gamma$) suffer from large
backgrounds arising from particle correlations (e.g. due to resonance decays)
coupled with the elliptic anisotropy. We propose differential measurements of
the $\Delta\gamma$ as a function of the pair invariant mass ($m_{\rm inv}$), by
restricting to high $m_{\rm inv}$ thus relatively background free, and by
studying the $m_{\rm inv}$ dependence to separate the possible CME signal from
backgrounds. We demonstrate by model studies the feasibility and effectiveness
of such measurements for the CME search.

We introduce two new integral transforms of the quantum mechanical transition
kernel that represent physical information about the path integral. These
transforms can be interpreted as probability distributions on particle
trajectories measuring respectively the relative contribution to the path
integral from paths crossing a given spatial point (the hit function) and the
likelihood of values of the line integral of the potential along a path in the
ensemble (the path averaged potential).

The quadratic form of the isobaric multiplet mass equation (IMME), which was
originally suggested by Wigner and has been generally regarded as valid, is
seriously questioned by recent highprecision nuclear mass measurements. The
usual resolution to this problem is to add empirically the cubic and quartic
$T_z$terms to characterize the deviations from the IMME, but finding the
origin of these terms remains an unsolved difficulty. Based on a strategy
beyond the Wigner's firstorder perturbation, we derive explicitly the cubic
and quartic $T_z$terms. These terms are shown to be generated by the effective
chargesymmetry breaking and chargeindependent breaking interactions in
nuclear medium combined with the Coulomb polarization effect. Calculations for
the $sd$ and lower $fp$shells explore a systematical emergence of the cubic
$T_z$term, suggesting a general deviation from the original IMME.
Intriguingly, the magnitude of the deviation exhibits an oscillationlike
behavior with mass number, modulated by the shell effect.

We investigate probing the hot and dense nuclear matter with strange vector
mesons ($K^*, \bar{K}^*$). Our analysis is based on PHSD which incorporates
partonic and hadronic dof and describes the full dynamics of HICs. This allows
to study the $K^*$ and $\bar{K}^*$ meson formation from the QGP and the
inmedium effects related to the modification of their properties during the
propagation in dense and hot matter. We employ relativistic BreitWigner
spectral functions for the $K^*,\bar{K}^*$ mesons with selfenergies obtained
from a Gmatrix approach to study the role of inmedium effects on the $K^*$
and $\bar{K}^*$ meson dynamics in HIC from FAIR/NICA to LHC energies. According
to our analysis most of the final $K^*/\bar{K}^*$s, that can be observed
experimentally, are produced during the late hadronic phase and stem dominantly
from the $K (\bar{K}) + \pi \to K^*(\bar{K}^*)$ formation channel. The amount
of $K^*/\bar{K}^*$s originating from the QGP channel is comparatively small
even at LHC energies and such $K^*/\bar{K}^*$s can hardly be reconstructed
experimentally due to the rescattering of final pions and (anti)kaons. This
mirrors the results from our previous study on the strange vectormeson
production in HICs at RHIC energies. The influence of the inmedium effects on
the dynamics of the $K^*/\bar{K}^*$ is rather small since they are mostly
produced at low baryon densities. Additional cuts on the shape of the observed
signal and the range of the invariant mass region of the $K^*/\bar{K}^*$ also
affect the final spectra. We demonstrate that the $K^*/\bar{K}^*$ inmedium
effects are more visible at lower beam energy, e.g. FAIR/NICA and BES RHIC
energies, where the production of $K^*/\bar{K}^*$s occurs at larger baryon
densities. Finally, we present the experimental procedures to extract
information on the inmedium masses and widths by fitting final mass spectra at
LHC energies.

We calculate the $\beta$decay of tritium at nexttoleading order in
pionless effective field theory. At this order, a lowenergy parameter $L_{1,
A}$ enters the calculation that is also relevant for a highaccuracy prediction
of the solar protonproton fusion rate. We use the tritium halflife to
determine this parameter and provide uncertainty estimates. We show proper
renormalization of our calculation analyzing the residual cutoff dependence of
observables. We find that nexttoleading order corrections contribute about
$4\%$ to the triton decay GamowTeller strength. We show that these conclusions
are insensitive to different arrangements of the effective range expansion.

The contribution of collective or reaction channels to a local optical model
potential, OMP, can be readily calculated as a dynamical polarization
potential, DPP. The resulting local DPPs commonly have undulatory (`wavy')
features, often including local regions of emissivity in the imaginary
component. We show here that this undularity arises from $l$dependence of the
underlying formal nonlocal and $l$dependent DPP. The $l$independent proton
OMPs, that have the same $S$matrix $S_{lj}$ as phenomenological $l$dependent
potentials, exhibit undulations that are qualitatively similar to undulations
of local DPPs generated by channel coupling. The $l$dependent phenomenological
potentials studied are the potentials that give the best existing fits to the
relevant elastic scattering data and the undulatory potentials presented here,
being $S$matrix equivalent (i.e.\ having the same $S$matrix, $S_{lj}$) give
exactly the same scattering. In addition, we present calculations strongly
suggesting that undularity (`waviness') is a generic property of
$l$independent potentials that are $S$matrix equivalent to $l$dependent
potentials. Implications for the validity of folding models based on a local
density model are noted.

We study the transverse momentum spectrum of hadrons in jets. By measuring
the transverse momentum with respect to a judiciously chosen axis, we find that
this observable is insensitive to (the recoil of) soft radiation. Furthermore,
for small transverse momenta we show that the effects of the jet boundary
factorize, leading to a new transversemomentumdependent (TMD) fragmentation
function. In contrast to the usual TMD fragmentation functions, it does not
involve rapidity divergences and is universal in the sense that it is
independent of the type of process and number of jets. These results directly
apply to subjets instead of hadrons. We discuss potential applications, which
include studying nuclear modification effects in heavyion collisions and
identifying boosted heavy resonances.

We present a lattice calculation of the nucleon isovector axial and induced
pseudoscalar form factors on the CLS ensembles using $N_{\rm f}=2$ dynamical
flavours of nonperturbatively $\mathcal{O}(a)$improved Wilson fermions and an
$\mathcal{O}(a)$improved axial current together with the pseudoscalar density.
Excitedstate effects in the extraction of the form factors are treated using a
variety of methods, with a detailed discussion of their respective merits. The
chiral and continuum extrapolation of the results is performed both using
formulae inspired by Heavy Baryon Chiral Perturbation Theory (HBChPT) and a
global approach to the form factors based on a chiral effective theory (EFT)
including axial vector mesons. Our results indicate that careful treatment of
excitedstate effects is important in order to obtain reliable results for the
axial form factors of the nucleon, and that the main remaining error stems from
the systematic uncertainties of the chiral extrapolation. As final results, we
quote $g_{\rm A} = 1.278 \pm 0.068\genfrac{}{}{0pt}{1}{+0.000}{0.087}$,
$\langle r_{\rm A}^2\rangle = 0.360 \pm
0.036\genfrac{}{}{0pt}{1}{+0.080}{0.088}~\mathrm{fm}^2$, and $g_{\rm P} = 7.7
\pm 1.8 \genfrac{}{}{0pt}{1}{+0.8}{2.0}$ for the axial charge, axial charge
radius and induced pseudoscalar charge, respectively, where the first error is
statistical and the second is systematic.

We show that the production of J/psi or Upsilon pairs in unpolarised pp
collisions is currently the best process to measure the momentumdistribution
of linearlypolarised gluons inside unpolarised protons through the study of
azimuthal asymmetries. Not only the shortdistance coefficients for such
reactions induce the largest possible cos 4 phi modulations, but analysed data
are already available. Among the various final states previously studied in
unpolarised pp collisions within the TMD factorisation approach, diJ/psi
production exhibits by far the largest asymmetries, up to 50% in the region
studied by the ATLAS and CMS experiments. In addition, we use the very recent
LHCb data at 13 TeV to perform the first fit of the unpolarised
transversemomentumdependent gluon distribution.

In this work, we have extracted the initial tempearture from the transverse
momentum spectra of charged particles in $Au+Au$ collisions using STAR data at
RHIC energies from $\sqrt{s_{NN}}$ = 7.7 to 200 GeV. The initial energy density
($\varepsilon$), shear viscosity to entropy density ratio ($\eta/s$), trace
anomaly ($\Delta$), the squared speed of sound ($C_s^2$), entropy density, and
bulk viscosity to entropy density ratio ($\zeta/s$) are obtained and compared
with the lattice QCD calculations for (2+1) flavor. The initial temperatures
obtained are compared with various hadronization and chemical freezeout
temperatures. The analysis of the data shows that the deconfinement to
confinement transition possibly takes place between $\sqrt{s_{NN}}$ = 11.5 and
19.6 GeV.

This short paper is an addendum to a recent publication on charged current
neutrinoinduced pion production (Phys.\ Rev.\ C96 (2017) no.1, 015503). It
presents comparisons of pion production cross sections measured at the T2K near
detector for a CH target.

We investigate photoproduction of $\omega$ off a nucleon target $\gamma
N\to\omega N$ by using a Reggeized model where
$\pi(135)+\sigma(500)+f_1(1285)+f_2(1270)$+Pomeron exchanges are included in
the $t$channel for the reaction at forward angles. The reaction mechanism at
low energy is featured by the dominance of the $\pi$ exchange with the
absorptive cuts introduced to modulate the pion contribution to both $\gamma
p\to \omega p$ and $\gamma n\to \omega n$ reactions. Necessity of $\sigma$
exchange is illustrated in the analysis of the cross section from natural
parity exchanges. Cross sections for differential, total, and spin density
matrix with beam polarization are reproduced and compared with existing data on
$\gamma p\to \omega p$. Scaled differential cross sections of Jefferson Lab
data on $\gamma p\to \omega p$ are investigated at the production angle
$\theta=90^\circ$ with the nonlinear trajectories for a saturation.
Differential and total cross sections for the $\gamma n\to\omega n$ reaction
are analyzed to compare with recent experimental data at the CBELSA/TAPS
Collaboration. An application to the $\gamma p\to\omega\Delta^+$ reaction is
presented to demonstrate the significance of the $\pi$ exchange in the total
and differential cross sections.

The spinindependent and transversity generalised form factors (GFFs) of the
$\phi$ meson are studied using lattice QCD calculations with light quark masses
corresponding to a pion mass $m_\pi\sim450(5)$ MeV. One transversity and three
spinindependent GFFs related to the lowest moments of leadingtwist
spinindependent and transversity gluon distributions are obtained at six
nonzero values of the momentum transfer up to 1.2 GeV$^2$. These quantities
are compared with the analogous spinindependent quark GFFs and the
electromagnetic form factors determined on the same lattice ensemble. The
results show quantitative distinction between the spatial distribution of
transversely polarised gluons, unpolarised gluons, and quarks, and point the
way towards further investigations of the gluon structure of nucleons and
nuclei.

Glauber models are used to calculate geometric quantities in the initial
state of heavy ion collisions, such as impact parameter, number of
participating nucleons and initial eccentricity. Experimental heavyion
collaboration, in particular at RHIC and LHC, use Glauber Model calculations
for various geometric observables. In this document, we describe the
assumptions inherent to the approach, and provide an updated implementation
(v2) of the Monte Carlo based Glauber Model calculation, which originally was
used by the PHOBOS collaboration. The main improvement w.r.t. the earlier
version (arXiv:0805.4411) are the inclusion of Tritium, Helium3, and Uranium,
as well as the treatment of deformed nuclei and GlauberGribov fluctuations of
the proton in p+A collisions. A users' guide (updated to reflect changes in v2)
is provided for running various calculations.

We present a generatorcoordinate method for realistic shellmodel
Hamiltonians that closely approximates the full shell model calculations of the
matrix elements for the neutrinoless double$\beta$ decay of $^{124}$Sn,
$^{130}$Te, and $^{136}$Xe. We treat not only quadrupole deformations but also
the protonneutron pairing amplitudes as generator coordinates. We validate
this method by calculating and comparing spectroscopic quantities with the
exact shell model results and experimental data. Our Hamiltonianbased
generatorcoordinate method produces $0\nu\beta\beta$ matrix elements much
closer to the shell model ones, compared to the existing
energydensityfunctionalbased generatorcoordinate approaches. The remaining
overestimation of $0\nu\beta\beta$ nuclear matrix element suggests that
additional correlations may be needed to be taken into account for $^{124}$Sn,
$^{130}$Te, and $^{136}$Xe when calculating with the Hamiltonianbased
generatorcoordinate method. The validation of this method may open the
possibility of calculating $0\nu\beta\beta$ matrix element of $^{150}$Nd in a
large shellmodel space.

We discuss the ground state of a quark matter containing heavy quarks as
impurities in a simple model which exhibits the QCD Kondo effect. The model
includes a currentcurrent interaction with the color exchange between a light
quark ($\psi$) and a heavy quark ($\Psi$). We introduce a gap function $\Delta
\sim \langle \bar \psi \Psi \rangle $ which represents the correlation between
$\psi$ and $\Psi$, and perform the meanfield approximation assuming that heavy
quarks are uniformly distributed. Values of the gap $\Delta$ measure the
strength of mixing between $\psi$ and $\Psi$. The gap equation obtained from
the minimum of the thermodynamical potential together with the condition for
the heavyquark number conservation turns out to allow for nonzero values of
the gap as the most stable state. We draw a phase diagram in $\mu$ (the
lightquark chemical potential) and $\lambda$ (an analog of the heavyquark
chemical potential) plane, and identify the region where the QCD Kondo effect
occurs.

The $DD^{*}$ potentials are studied within the framework of heavy meson
chiral effective field theory. We have obtained the effective potentials of the
$DD^{*}$ system up to $O(\epsilon^2)$ at one loop level. In addition to the
onepion exchange contribution, the contact and twopion exchange interactions
are also investigated in detail. Furthermore, we have searched for the possible
molecular states by solving Schr\"odinger equation with the potentials. We
notice that the contact and twopion exchange potentials are nonnegligible
numerically and important for the existence of a bound state. In our results,
no bound state is founded in the $I=0$ channel within a wide range of cutoff
parameter, while there exists a bound state in the $I=1$ channel as cutoff is
near $m_\rho$ in our approach.

We argue that the existence of an inhomogeneous external magnetic field can
lead to radial flow in transverse plane. Our aim is to show how the
introduction of a magnetic field generalizes the Bjorken flow. We investigate
the effect of an inhomogeneous weak external magnetic field on the transverse
expansion of inviscid fluid created in high energy nuclear collisions. In
order to simplify our calculation and compare with Gubser model, we consider
the fluid under investigation to be produced in central collisions, at small
impact parameter; azimuthal symmetry has been considered. In our model, we
assume an inhomogeneous external magnetic field following the powerlaw decay
in proper time and having radial inhomogeneity perpendicular to the radial
velocity of the inviscid fluid in the transverse plane; then the space time
evolution of the transverse expansion of the fluid is obtained. We also show
how the existence of an inhomogeneous external magnetic field modifies the
energy density. Finally we use the solutions for the transverse velocity and
energy density in the presence of a weak magnetic field, to estimate the
transverse momentum spectrum of protons and pions emerging from the
Magnetohydrodynamic solutions.

We have analyzed the data on $\psi'$ production in protonnucleus ($p+A$)
collisions, available from the NA50 Collaboration in the SPS energy domain. The
investigated data sets include the absolute production cross sections as well
as $\psi'$toDrell Yan (DY) cross section ratios. An adapted version of two
component QVZ model has been employed to calculate $\psi'$ production cross
sections. For both $\psi'$ and DY production, nuclear modifications to the free
nucleon structure functions are taken into account. For $\psi'$, final state
interaction of the produced $c\bar{c}$ pairs with the nuclear medium is also
taken into account, in accordance with the previously analyzed $J/\psi$ data. A
satisfactory description of the data in $p+A$ collisions is obtained. Model
calculations are extrapolated to low energy regime to predict the $\psi'$
suppression in proton induced collisions at near threshold energies.

The neutron skin and halo distributions in medium and heavy nuclei are
calculated within the extended ThomasFermi approximation. Calculations are
carried out for the effective Skyrmelike forces using the direct variational
method. The analytical expression for the isovector shift of the rms radii
$\Delta r_{np}$ as a sum of skin and halolike terms is obtained. The
contribution of halo and skin terms to $\Delta r_{np}$ are found to be
approximately equal.