
Searching for the Neutrinoless Double Beta Decay (NLDBD) is now regarded as
the topmost promising technique to explore the nature of neutrinos after the
discovery of neutrino masses in oscillation experiments. PandaXIII (Particle
And Astrophysical Xenon Experiment III) will search for the NLDBD of $^{136}$Xe
at the China Jin Ping underground Laboratory (CJPL). In the first phase of the
experiment, a high pressure gas Time Projection Chamber (TPC) will contain 200
kg, 90% $^{136}$Xe enriched gas operated at 10 bar. Fine pitch micropattern
gas detector (Microbulk Micromegas) will be used at both ends of the TPC for
the charge readout with a cathode in the middle. Charge signals can be used to
reconstruct tracks of NLDBD events and provide good energy and spatial
resolution. The detector will be immersed in a large water tank to ensure
$\sim$5 m of water shielding in all directions. The second phase, a tonscale
experiment, will consist of five TPCs in the same water tank, with improved
energy resolution and better control over backgrounds.

In PYTHIA6 (PYTHIA8) once the transverse momentum $p_T$ of a generated
particle is randomly sampled, $p_x$ and $p_y$ are set on the circle with radius
of $p_T$ randomly. This may largely suppress the development of the final
hadronic state transverse momentum anisotropy from the initial state spatial
asymmetry. We modify PYTHIA6.4 by randomly setting $p_x$ and $p_y$ on the
circumference of an ellipse with the half major and minor axes being
$p_T(1+\delta_p)$ and $p_T(1\delta_p)$, respectively. The modified PYTHIA6.4
is then employed to systematically study the charged particle transverse
momentum asymmetry in the minimum bias pp collisions at $\sqrt s$=0.9, 7, and
14 TeV. The ALICE data on the transverse sphericity as a function of charged
multiplicity, $<S_T>(N_{ch})$, are well reproduced with the modified PYTHIA6.4.
It is found that the predicted charged particle $v_2$ upper limit is a
measurable value, $\sim 0.12$, in the minimum bias pp collisions at $\sqrt s$=7
TeV. We suggest a systematic measurement for the particle transverse momentum
sphericity, eccentricity (ellipticity), and elliptic flow parameter.

We calculate the real and nonstatistical higher moment excitation functions
($\sqrt{s_{NN}}$ = 11.5 to 200 GeV) for the netproton, netbaryon, and the
netcharge number event distributions in the relativistic Au+Au collisions with
the parton and hadron cascade model PACIAE. It turned out that because of the
statistical fluctuation dominance it is very hard to see signature of the CP
singularity in the real higher moment excitation functions. It is found that
the property of higher moment excitation functions are significantly dependent
on the window size, and hence the CP signatures may show only in a definite
window for a given conserved observable nonstatistical higher moments. But for
a given widow size, the CP singularity may show only in the nonstatistical
higher moment excitation functions of a definite conserved observable.

We calculate the heavy quark evolution in heavy ion collisions and show
results for the elliptic flow $v_2$ as well as the nuclear modification factor
$R_{AA}$ at RHIC and LHC energies. For the calculation we implement a Langevin
approach for the transport of heavy quarks in the UrQMD (hydrodynamics +
Boltzmann) hybrid model. As drag and diffusion coefficients we use a Resonance
approach for elastic heavyquark scattering and assume a decoupling temperature
of the charm quarks from the hot medium of $130\,\text{MeV}$. At RHIC energies
we use a coalescence approach at the decoupling temperature for the
hadronization of the heavy quarks to Dmesons and Bmesons and a subfollowing
decay to heavy flavor electrons using PYTHIA. At LHC we use an additional
fragmentation mechanism to account for the higher transverse momenta reached at
higher collision energies.

We use the nonstatistical fluctuation instead of the full one to explore the
higher moment singularities of net proton event distributions in the
relativistic Au+Au collisions at $\sqrt{s_{NN}}$ from 11.5 to 200 GeV
calculated by the parton and hadron cascade model PACIAE. The PACIAE results of
mean ($M$), variance ($\sigma^2$), skewness ($S$), and kurtosis ($\kappa$) are
consistent with the corresponding STAR data. Nonstatistical moments are
calculated as the difference between the moments derived from real events and
the ones from mixed events, which are constructed by combining particles
randomly selected from different real events. An evidence of singularity at
$\sqrt{s_{NN}}\sim$ 60 GeV is first seen in the energy dependent
nonstatistical $S$ and $S\sigma$.

The influence of the chiral mean field on the $K^+$ transverse flow in heavy
ion collisions at SIS energy is investigated within covariant kaon dynamics.
For the kaon mesons inside the nuclear medium a quasiparticle picture
including scalar and vector fields is adopted and compared to the standard
treatment with a static potential. It is confirmed that a Lorentz force from
spatial component of the vector field provides an important contribution to the
inmedium kaon dynamics and strongly counterbalances the influence of the
vector potential on the $K^+$ inplane flow. The FOPI data can be reasonably
described using inmedium kaon potentials based on effective chiral models. The
information on the inmedium $K^+$ potential extracted from kaon flow is
consistent with the knowledge from other sources.

The $\phi$ meson productions in $Au+Au$ and/or $Pb+Pb$ collisions at AGS,
SPS, RHIC, and LHC energies have been studied systematically with a hadron and
string cascade model LUCIAE. After considering the energy dependence of the
model parameter $\alpha$ in string fragmentation function and adjusting it to
the experimental data of charged multiplicity to a certain extent, the model
predictions for $\phi$ meson yield, rapidity, and/or transverse mass
distributions are compatible with the experimental data at AGS, SPS and RHIC
energies. A calculation for $Pb+Pb$ collisions at LHC energy is given as well.
The obtained fractional variable in string fragmentation function shows a
saturation in energy dependence. It is discussed that the saturation of
fractional variable in string fragmentation function might be a qualitative
representation of the energy dependence of nuclear transparency.