
We present a detailed analysis of the beam energy dependence of the
mechanisms for generation of directed and elliptic flows in Au+Au collisions
focusing on the role of hadronic rescattering and spectator shadowing within a
microscopic transport model JAM with different equation of state. A systematic
study of the beam energy dependence is performed for Au+Au collisions at $2
\leq \sqrt{s_{NN}}\leq 62.4$ GeV. Transition of the dynamical origin of
directed flow is observed. We found that initial Galuber type nucleonnucleon
collisions generate negative $v_1$ for nucleons at midrapidity due to the
presence of spectator matter, and this negative nucleon $v_1$ is turned to be
positive by the mesonbaryon interactions at the beam energy region of
$\sqrt{s_{NN}} < 30$ GeV. In contrast, above 30 GeV, there is no spectator
shadowing at midrapidity, and initial nucleonnucleon collisions does not
generate directed flow, but subsequent rescatterings among produced particles
generate negative $v_1$ for nucleons. It is demonstrated that negative pion
directed flows are mostly generated by the interaction with the spectator
matter. It is also shown that squeezeout effect is largely suppressed in the
case of softening, which leads to the enhancement of elliptic flow around
$\sqrt{s_{NN}}=57$ GeV. Elliptic flow at midrapidity above 10 GeV is not
influenced by the squeezeout due to spectator matter, while its effect is seen
at the forward rapidity range of $y/y_\mathrm{cm}>0.5$, which decreases as beam
energy increases.

We performed systematic studies for the effects of eventbyevent efficiency
fluctuations on efficiency correction for cumulant analysis in relativistic
heavyion collision experiments. Experimentally, particle efficiencies of
events measured with different experimental conditions should be different. For
fluctuation measurements, the final eventbyevent multiplicity distributions
should be the superposed distributions of various type of events measured under
different conditions.We demonstrate the efficiency fluctuation effects in a
numerical simulation, in which we construct an event ensemble consists of
events with two different efficiencies. By using the mean particle
efficiencies, we find the efficiency corrected cumulants show large deviations
from the original inputs when the discrepancy of the two efficiencies are
large. We further studied the effects of efficiency fluctuations for the
cumulants of netproton distributions by implementing the UrQMD events of Au+Au
collisions at $\sqrt{s_{NN}}=7.7$ GeV in a realistic STAR detector acceptance.
We consider the unequal efficiency in two sides of the Time Projection Chamber
(TPC), multiplicity dependent efficiency and the eventbyevent variations of
the collision vertex position along the longitudinal direction
($V_\mathrm{z}$). When the efficiencies fluctuate dramatically within the
studied event sample, the effects of efficiency fluctuations have significant
impacts on the efficiency corrections of cumulants with the mean efficiencies.
We find that this effect can be effectively suppressed by binning the entire
event ensemble into various subevent samples, in which the efficiency
variations are relative small. The final efficiency corrected cumulants can be
calculated from the weighted average of the corrected factorial moments of the
subevent samples with the mean efficiencies.

We perform a systematic study of elliptic flow ($v_2$) in Au+Au collision at
$\sqrt{s_{NN}}=5$ GeV by using a microscopic transport model JAM. The
centrality, pseudorapidity, transverse momentum and beam energy dependence of
$v_2$ for charged as well as identified hadrons are studied. We investigate the
effects of both the hadronic meanfield and the softening of equation of state
(EoS) on elliptic flow. The softening of EoS is realized by imposing attractive
orbits in two body scattering, which can reduce the pressure of the system. We
found that the softening of EoS leads to the enhancement of $v_2$, while the
hadronic meanfield suppresses $v_2$ relative to the cascade mode. It indicates
that elliptic flow at high baryon density regions is highly sensitive to the
EoS and the enhancement of $v_2$ may probe the signature of a firstorder phase
transition in heavyion collisions at beam energies of a strong baryon stopping
region.

We studied the acceptance dependence of proton cumulants (up to fourth order)
and correlation functions in 05\% most central Au+Au collisions at
$\sqrt{s_\text{NN}}$=7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV from UrQMD
model. We found that high order proton cumulants show suppressions at large
acceptance. By decomposing the proton cumulants into linear combination of
multiproton correlation functions, we observed the twoproton correlation
functions always show negative values due to the effects of baryon number
conservations. The three and fourproton correlation functions are close to
zero and show negligible acceptance dependence. We further observed that the
proton cumulants and correlation functions follow similar trends and show a
scaling behavior when plotting the results versus mean number of protons. The
comparisons between experimental data and the UrQMD calculations show that the
nonmonotonic energy dependence of proton correlation functions measured by
STAR experiment cannot be described by the UrQMD model. The UrQMD calculations
can provide us baselines for the experimental studies of the proton cumulants
and correlation functions. Finally, we propose to measure the rapidity
dependence of the reduced proton correlation functions to search for the QCD
critical point in heavyion collisions.

Fluctuations of conserved quantities, such as baryon, electric charge and
strangeness number, are sensitive observables in relativistic heavyion
collisions to probe the QCD phase transition and search for the QCD critical
point. In this paper, we review the experimental measurements of the cumulants
(up to fourth order) of eventbyevent netproton (proxy for netbaryon),
netcharge and netkaon (proxy for netstrangeness) multiplicity distributions
in Au+Au collisions at $\sqrt{s_{NN}}=7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, 200$
GeV from the first phase of beam energy scan program at the Relativistic
HeavyIon Collider (RHIC). We also summarize the data analysis methods of
suppressing the volume fluctuations, autocorrelations and the unified
description of efficiency correction and error estimation. Based on theoretical
and model calculations, we will discuss the characteristic signatures of
critical point as well as backgrounds for the fluctuation observables in
heavyion collisions. The physics implications and the future second phase of
the beam energy scan (20192020) at RHIC will be also discussed.

We provide a unified description of efficiency correction and error
estimation for moments of conserved quantifies in heavyion collisions. Moments
and cumulants are expressed in terms of the factorial moments, which can be
easily corrected for the efficiency effect. By deriving the covariance between
factorial moments, one can obtain the general error formula for the efficiency
corrected moments based on the error propagation derived from the Delta
theorem. The skellam distribution based Monto Carlo simulation is used to test
the Delta theorem and Bootstrap error estimation methods. The statistical
errors calculated from the two methods can well reflect the statistical
fluctuations of the efficiency corrected moments.

The beam energy dependence of the elliptic flow,$v_2$, is studied in
midcentral Au+Au collisions in the energy range of $3\leq \sqrt{s_{NN}} \leq
30$ GeV within the microscopic transport model JAM. The results of three
different modes of JAM are compared; cascade,hadronic mean field, and a new
mode with modified equations of state, with a first order phase transition
(1.O.P.T.) and with a crossover transition. The standard hadronic mean field
suppresses $v_2$, while the inclusion of the effects of a 1.O.P.T. (and also of
a crossover transition) does enhance $v_2$ at $\sqrt{s_{NN}}<30$ GeV. The
enhancement or suppression of the scaled energy flow, dubbed "elliptic flow"is
understood as being due to out of plane flow, i.e. $v_2<0$, dubbed out of
plane  "squeezeout", which occurs predominantly in the early, compression
stage. Subsequently, the inplane flow dominates, in the expansion stage, $v_2
> 0$. The directed flow, dubbed "bounce off", is an independent measure of the
pressure, which quickly builds up the transverse momentum transfer in the
reaction plane. When the spectator matter leaves the participant fireball
region, where the highest compression occurs, a hard expansion leads to larger
$v_2$. A combined analysis of the three transverse flow coefficients, radial
$v_0$, directed $v_1$ and elliptic $v_2$ flow, in the beam energy range of
$3\leq\sqrt{s_{NN}}\leq10$ GeV, distinguishes the different compression and
expansion scenarios: a characteristic dependence on the early stage equation of
state is observed. The enhancement of both the elliptic and the transverse
radial flow and the simultaneous collapse of the directed flow of nucleons
offers a clear signature if 1.O.P.T. is realized at the highest baryon
densities created in high energy heavyion collisions.

Fluctuations of conserved quantities, such as baryon, electric charge and
strangeness number, are sensitive observables in heavyion collisions to search
for the QCD phase transition and critical point. In this paper, we performed a
systematical analysis on the various cumulants and cumulant ratios of
eventbyevent netstrangeness distributions in Au+Au collisions at
$\sqrt{s_{NN}}$=7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV from UrQMD model. We
performed a systematical study on the contributions from various strange
baryons and mesons to the netstrangeness fluctuations. The results demonstrate
that the cumulants and cumulant ratios of netstrangeness distributions
extracted from different strange particles show very different centrality and
energy dependence behavior. By comparing with the netkaon fluctuations, we
found that the strange baryons play an important role in the fluctuations of
netstrangeness. This study can provide useful baselines to study the QCD phase
transition and search for the QCD critical point by using the fluctuations of
netstrangeness in heavyion collisions experiment. It can help us to
understand noncritical physics contributions to the fluctuations of
netstrangeness.

We discuss properties and applications of factorial cumulants of various
particle numbers and for their mixed channels measured by the eventbyevent
analysis in relativistic heavyion collisions. After defining the factorial
cumulants for systems with multiparticle species, their properties are
elucidated. The uses of the factorial cumulants in the study of critical
fluctuations are discussed. We point out that factorial cumulants play useful
roles in understanding fluctuation observables when they have underlying
physics approximately described by the binomial distribution. As examples, we
suggest novel utilization methods of the factorial cumulants in the study of
the momentum cut and rapidity window dependences of fluctuation observables.

Baryon number fluctuations are sensitive to the QCD phase transition and QCD
critical point. According to the Feynman rules of finitetemperature field
theory, we calculated various order moments and cumulants of the baryon number
distributions in the quasiparticle model of quark gluon plasma. Furthermore,
we compared our results with the experimental data measured by the STAR
experiment at RHIC. It is found that the experimental data can be well
described by the model for the colliding energies above 30 GeV and show large
discrepancies at low energies. It can put new constraint on qQGP model and also
provide a baseline for the QCD critical point search in heavyion collisions at
low energies.

Substantial experimental and theoretical efforts worldwide are devoted to
explore the phase diagram of strongly interacting matter. At LHC and top RHIC
energies, QCD matter is studied at very high temperatures and nearly vanishing
netbaryon densities. There is evidence that a QuarkGluonPlasma (QGP) was
created at experiments at RHIC and LHC. The transition from the QGP back to the
hadron gas is found to be a smooth cross over. For larger netbaryon densities
and lower temperatures, it is expected that the QCD phase diagram exhibits a
rich structure, such as a firstorder phase transition between hadronic and
partonic matter which terminates in a critical point, or exotic phases like
quarkyonic matter. The discovery of these landmarks would be a breakthrough in
our understanding of the strong interaction and is therefore in the focus of
various highenergy heavyion research programs. The Compressed Baryonic Matter
(CBM) experiment at FAIR will play a unique role in the exploration of the QCD
phase diagram in the region of high netbaryon densities, because it is
designed to run at unprecedented interaction rates. Highrate operation is the
key prerequisite for highprecision measurements of multidifferential
observables and of rare diagnostic probes which are sensitive to the dense
phase of the nuclear fireball. The goal of the CBM experiment at SIS100
(sqrt(s_NN) = 2.7  4.9 GeV) is to discover fundamental properties of QCD
matter: the phase structure at large baryonchemical potentials (mu_B > 500
MeV), effects of chiral symmetry, and the equationofstate at high density as
it is expected to occur in the core of neutron stars. In this article, we
review the motivation for and the physics programme of CBM, including
activities before the start of data taking in 2022, in the context of the
worldwide efforts to explore highdensity QCD matter.

For the first time, we investigate susceptibilities of dense quark matter up
to $8$th order using an effective model. Generally higher order
susceptibilities will have more sign changes and larger magnitude, thus should
give more information about the presence and location of the conjectured QCD
critical end point (CEP). Two cases are studied, one with the CEP and one being
crossover phase transition throughout the QCD phase diagram. It is found that a
rapid crossover transition can also give similar sign pattern as the case with
the CEP and yield large signal. We propose that using several kinds of ions for
collision to search for a wider range on the QCD phase diagram may help.

Fluctuations and correlations of conserved charges are sensitive observables
for studying the QCD phase transition and critical point in highenergy
heavyion collisions. We have studied the centrality and energy dependence of
mixedcumulants (up to fourth order) between netbaryon and netstrangeness in
Au+Au collisions at $\sqrt{s_{NN}}$= 7.7, 11.5, 19.6, 27, 39, 62.4, 200 GeV
from UrQMD model. To compare with other theoretical calculations, we normalize
these mixedcumulants by various order cumulants of netstrangeness
distributions. We found that the results obtained from UrQMD calculations are
comparable with the results from Lattice QCD at low temperature and hadron
resonance gas model. The ratios of mixedcumulants
($R_{11}^{BS},R_{13}^{BS},R_{22}^{BS}$,~$R_{31}^{BS}$) from UrQMD calculations
show weak centrality dependence. However, the mixedcumulant ratios
$R_{11}^{BS}$ and $R_{31}^{BS}$ show strong increase at low energy, while the
$R_{13}^{BS}$ snd $R_{22}^{BS}$ are similar at different energies. Furthermore,
we have also studied the correlations between different hadron species and
their contributions to the netbaryon and netstrangeness correlations. These
models studies can provide baselines for searching for the signals of QCD phase
transition and critical point in heavyion collisions.

We evaluate the $2$nd to $4$th order baryon, charge and strangeness
susceptibilities near a chiral critical point using the NambuJonaLasinio
model under different temperature and baryon chemical potential. Baryon number
susceptibilities are found to be of the greatest magnitude, offering the
strongest signal. Whereas the strangeness susceptibilities have the smallest
divergence dominating area, owing to the large strange quark mass. We also make
an attempt to compare our results with experiment data. The trend at high
collision energy are found to be consistent between theory and experiment. The
model calculation predicts more complex behavior at low collision energies,
near the postulated critical end point.

We analyze the rapidity and transverse momentum dependence for the cumulants
of the netproton and netbaryon distributions in Au+Au collisions at
$\sqrt{s_{\text{NN}}} = 5\,\text{GeV}$ with a microscopic hadronic transport
(JAM) model. To study the effects of mean field potential and softening of
equation of state (EoS) on the fluctuations of netproton (baryon) in heavyion
collisions, the calculations are performed with two different modes. The
softening of EoS is realized in the model by implementing the attractive orbit
in the twobody scattering to introduce a reduction pressure of the system. By
comparing the results from the two modes with the results from default cascade,
we find the mean field potential and softening of EoS have strong impacts on
the rapidity distributions ($\text{d}N/\text{d}y$) and the shape of the
netproton (baryon) multiplicity distributions. The netproton (baryon)
cumulants and their ratios calculated from all of the three modes are with
similar trends and show significant suppression with respect to unity, which
can be explained by the presence of baryon number conservations. It indicates
that the effects of mean field potential and softening of EoS might be not the
ingredients that are responsible to the observed strong enhancement in the most
central Au+Au collisions at 7.7 GeV measured by the STAR experiment at RHIC.

Fluctuations of conserved quantities are predicted to be sensitive to the
correlation length and connected to the thermodynamic susceptibility. Thus,
moments of netbaryon, netcharge and netstrangeness have been extensively
studied theoretically and experimentally to explore phase structure and bulk
properties of QCD matters created in heavy ion collision experiment. As the
moments analysis is statistics hungry study, the error estimation is crucial to
extract physics information from the limited experimental data. In this paper,
we will derive the limit distributions and error formula based on Delta theorem
in statistics for various order moments used in the experimental data analysis.
The Monte Carlo simulation is also applied to test the error formula.

Fluctuations of conserved quantities are sensitive observables to probe the
signature of QCD phase transition and critical point in heavyion collisions.
With the UrQMD model, we have studied the centrality and energy dependence of
various order cumulants and cumulant ratios (up to fourth order) of
netproton,netcharge and netkaon multiplicity distributions in Au+Au
collisions at $\sqrt{s_{NN}}$= 7.7, 11.5, 19.6, 27, 39, 62.4, 200 GeV. The
model results show that the production mechanism of the particles and
antiparticles have significant impacts on the cumulants of netparticles
multiplicity distributions and show strong energy dependence. We also made
comparisons between model calculations and experimental data measured in the
first phase of the beam energy scan (BES) program by the STAR experiment at
RHIC. The comparisons indicate that the baryon conservation effect strongly
suppress the cumulants of netproton distributions at low energies and the
nonmonotonic energy dependence for the netproton {\KV} at the most central
Au+Au collisions measured by the STAR experiment can not be described by the
UrQMD model. Since there has no physics of QCD phase transition and QCD
critical point implemented in the UrQMD, the model results provide us baselines
and qualitative estimates about the noncritical background contributions to
the fluctuations observables in heavyion collisions.

We show the importance of separated efficiency corrections between positively
and negatively charged particles for cumulant calculations by Monte Carlo toy
models and analytical calculations. Our results indicate that S{\sigma} in
published netproton results from the STAR experiment will be suppressed about
5 to 10% in central collisions, and 10 to 20% in peripheral collisions at the
beam energy of \sqrt s_{NN} = 62.4 and 200 GeV if the separated efficiencies
are used to efficiency correction.

Beam energy scan programs in heavyion collisions aim to explore the QCD
phase structure at high baryon density. Sensitive observables are applied to
probe the signatures of the QCD phase transition and critical point in
heavyion collisions at RHIC and SPS. Intriguing structures, such as dip, peak
and oscillation, have been observed in the energy dependence of various
observables. In this paper, an overview is given and corresponding physics
implications will be discussed for the experimental highlights from the beam
energy scan programs at the STAR, PHENIX and NA61/SHINE experiments.
Furthermore, the beam energy scan phase II at RHIC (20192020) and other future
experimental facilities for studying the physics at low energies will be also
discussed.

We present the energy dependence of moments of netproton and netcharge
multiplicity distributions in Au+Au collisions measured by the STAR experiment
in the first phase of the Beam Energy Scan (BES) at the Relativistic Heavy Ion
Collider (RHIC). By using the time of flight detector for particle
identification, the upper transverse momentum ($p_{T}$) limit for proton and
antiproton is extended from 0.8 GeV/c up to 2 GeV/c. The $p_{T}$ and rapidity
acceptance dependence study for the moments of netproton distribution show
that the larger the acceptance is, the greater the deviation from unity. The
most pronounced structure is found in the energy dependence of {\KV} of
netproton distributions from the $0\sim5%$ most central collisions within
$0.4<p_{T}<2$ GeV/c and at midrapidity $y<0.5$. At energies above 39 GeV,
the values of {\KV} are close to unity and for energies below 39 GeV, it shows
significant deviation below unity around 19.6 and 27 GeV, then a large increase
above unity is observed at 7.7 GeV.

Fluctuations of conserved quantities in heavyion collisions are used to
probe the phase transition and the QCD critical point for the strongly
interacting hot and dense nuclear matter. The STAR experiment has carried out
moment analysis of netproton (proxy for netbaryon (B)), netkaon (proxy for
netstrangeness (S)), and netcharge (Q). These measurements are important for
understanding the quantum chromodynamics phase diagram. We present the analysis
techniques used in the moment analysis by the STAR experiment and discuss the
moments of netproton and netcharge distributions from the first phase of the
Beam Energy Scan program at the Relativistic Heavy Ion Collider.

We present a systematic comparison between the recently measured cumulants of
the netproton distributions by STAR for 05% central Au+Au collisions at
$\sqrt{s_{NN}}$=7.7200 GeV and two kinds of possible baseline measures, the
Poisson and Binomial baselines. These baseline measures are assuming that the
proton and antiproton distributions independently follow Poisson statistics or
Binomial statistics. The higher order cumulant netproton data are observed to
deviate from all the baseline measures studied at 19.6 and 27 GeV. We also
compare the netproton with netbaryon fluctuations in UrQMD and AMPT model,
and convert the netproton fluctuations to netbaryon fluctuations in AMPT
model by using a set of formula.

Moments (Variance ($\sigma^2$), Skewness($S$), Kurtosis($\kappa$)) of
multiplicity distributions of conserved quantities, such as
netbaryon,netcharge and netstrangeness, are predicted to be sensitive to the
correlation length of the system and connected to the thermodynamic
susceptibilities computed in Lattice QCD and Hadron Resonance Gas (HRG) model.
In this paper, we present several measurement artifacts that could lead to
volume fluctuation and autocorrelation effects in the moment analysis of
netproton multiplicity distributions in heavyion collisions using the UrQMD
model. We discuss methods to overcome these artifacts so that the extracted
moments could be used to obtain physical conclusions. In addition we present
methods to properly estimate the statistical errors in moment analysis.

We present various order moments (Variance ($\sigma^2$), Skewness($S$),
Kurtosis($\kappa$)) of netproton multiplicity distributions in Au+Au
collisions from the first phase of RHIC beam energy scan ({\sNN}=7.7$$200
GeV). The measurements are carried out by the STAR detector at midrapidity
($y<0.5$) and within transverse momentum range $0.4<p_{T}<0.8$ {\gevc}. The
product of moments ({\SD} and \KV) are predicted to be sensitive to the
correlation length and connected to ratios of baryon number susceptibility. We
observe deviations below Poisson expectations in the most central collisions
for all of the energies. These results are compared with a transport model to
understand the effects not related to critical physics. We also discuss the
error estimation methods and the techniques to suppress centrality resolution
effect in the moment analysis.

We present the measurements of the higher moments of the netproton
multiplicity distributions for Au+Au collisions from the first phase of the
beam energy scan program at RHIC. The measurements are carried out at
midrapidity ($y<0.5$) and within the transverse momentum range
$0.4<p_{T}<0.8$ {\gevc} by STAR. The moment products, $\KV$ and $\SD$, which
are related to the ratios of various order baryon number susceptibilities and
sensitive to the correlation length, are applied to search for the QCD critical
point. We report the centrality and energy dependence of the moment products
({\KV} and {\SD}), results are compared with Poisson expectations. The
experiment data are also compared with transport model calculation to
understand the noncritical physics effects.