
An analytic formula is proposed to characterize the variance propagation from
correlated input variables to the model response, by using multivariate Taylor
series. With the formula, partial variance contributions to the model response
are then straightforwardly evaluated in the presence of input correlations.
Additionally, an arbitrary variable is represented as the sum of independent
and correlated parts. Universal expressions of the coefficients that specify
the correlated and independent sections of a single variable are derived by
employing linear correlation model. Based on the coefficients, it is nature to
quantify the independent, correlated and coupling contributions to the total
variance of model response. Numerical examples suggest the effectiveness and
validation of our analytic framework for general models. A practical
application of the analytic framework is also proposed to the sensitivity
analysis of a deterministic HIV model.

This paper presents the harmonic state space (HSS) modeling of a threephase
modular multilevel converter (MMC). MMC is a converter system with a typical
multifrequency response due to its significant harmonics in the arm currents,
capacitor voltages, and control signals. These internal harmonic dynamics can
have a great influence on the operation characteristics of MMC. However, the
conventional modeling methods commonly used in twolevel voltagesource
converters (VSCs), where only the fundamentalfrequency dynamic is considered,
will lead to an inaccurate model that cannot accurately reflect the real
dynamic characteristics of MMC. Therefore, the HSS modeling method, in which
harmonics of state variables, inputs, and outputs are posed separately in a
statespace form, is introduced in this paper to model the MMC in order to
capture all the harmonics and the frequency couplings. The steadystate and
smallsignal dynamic HSS models of a threephase MMC are developed,
respectively. The validity of the developed HSS model of a threephase MMC has
been verified by the results from both the nonlinear time domain simulation
model in MATLAB/Simulink and the laboratory prototype with 12 submodules per
arm.

Impedancebased stability analysis in the dqdomain is a widely applied
method for power electronic dominated systems. An inconvenient property with
this method is that impedance models are normally referred to their own local
reference frame, and need to be recalculated when referring to a global
reference frame in a given network. This letter presents a simple method for
translating impedance submodels within a complex network, from their own
reference frames to any given point in the network. What distinguishes this
method is that by using a simple rotational matrix, it only needs impedance
models in their own local reference frames, to be translated to a global
reference in the network. By way of this method, standard circuit analysis
rules for series and parallel connection are applicable, as proven in the
letter. The method is defined and validated for impedances in the dq and
modified sequence domains, and it is shown that the dependency on reference
frame is marginal in the latter. An additional finding from the application of
this method is that components or subsystems with a certain symmetry property
called Mirror Frequency Decoupling are invariant to the choice of reference
frame.

The smallsignal impedance modeling of modular multilevel converter (MMC) is
the key for analyzing resonance and stability of MMCbased ac power electronics
systems. MMC is a converter system with a typical multifrequency response due
to its significant steadystate harmonic components in the arm currents,
capacitor voltages, and control signals. Therefore, traditional smallsignal
modeling methods for 2level voltagesource converters (VSCs) cannot be
directly applied to the MMC. In this paper, the harmonic statespace (HSS)
modeling approach is introduced to characterize the harmonic coupling behavior
of the MMC. On this basis, the smallsignal impedance models of the MMC are
developed according to the harmonic linearization principle, which can include
all the steadystate harmonic effects of the state variables, leading to the
accurate impedance models. Furthermore, in order to reveal the impact of the
internal dynamics and closedloop control on the smallsignal impedance of the
MMC, three cases are considered in this paper, i.e., openloop control, ac
voltage closedloop control, and circulating current closedloop control.
Finally, the analytical impedance models are verified by both simulation and
experimental results.

Impedance models are widely used in assessing small signal stability of
gridtied voltage source converter (VSC) systems. Recent research has proven
that impedance models of gridtied VSC in both dq and sequence domains are
generally MultiInput MultiOutput (MIMO) systems, and the generalized Nyquist
criterion has to be applied for stability analysis to these MIMO systems.
However, finding SingleInput and SingleOutput (SISO) equivalents for this
system is always appealing because of the simplicity and the convenience for
physical interpretation when assessing the stability, compared to MIMO systems.
This paper presents two types of SISO impedance models of gridtied VSC system,
one is derived from the strong grid assumption, and the other is from the
closedloop equivalence. The accuracy of these models is assessed with respect
to the measured impedances in PSCAD/EMTDC, and their effects on the stability
assessment were analyzed as well. It is proven that the accurate SISO model
gives identical result as the MIMO (matrixbased) impedance model with respect
to stability analysis. However, the reduced SISO model may lead to wrong
results if the bandwidth of phase locked loop is large.

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.

Impedancebased stability analysis is appealing in the case of SingleInput
SingleOutput (SISO) systems. However, in the case of gridtied voltage source
converter (VSC) systems, dq impedances of source and load (VSC) subsystems are
typically MultiInput MultiOutput (MIMO) systems in which case the Generalized
Nyquist Criterion (GNC) is required for analyzing the closed loop stability,
which increases the complexity of the analysis. This paper explores further the
coupling between positive and negative sequence impedances, in particular the
dependency and bindings between them. It shows that the couplings in each
subsystem can be compounded into two noncoupled sequence impedances if the
source and load subsystems are viewed as an integrated system instead of as two
separate subsystems. Therefore, two decoupled SISO systems are obtained which
are defined as Augmented Sequence Impedance Networks (ASIN). The stability
analysis of the closed loop system is performed directly on the ASIN with the
principal of Argument. Both numerical and time domain verifications are
presented.

This paper provides a comparative analysis of impedance models for power
electronic converters and systems for the purpose of stability investigations.
Such models can be divided into either decoupled models or matrix models. A
decoupled impedance model is highly appealing since the
SingleInputSingleOutput (SISO) structure makes the analysis and result
interpretation very simple. On the other hand, matrix impedance models are more
accurate, and in some cases necessary. Previous works have applied various
approximations to obtain decoupled models, and both the dq and sequence
domains have been used. This paper introduces the terms decoupled and
semidecoupled impedance models in order to have a clear classification of the
available approximations.
The accuracy of 4 decoupled impedance models are discussed based on the
concept of Mirror Frequency Coupling (MFC). By definition the decoupled models
based on sequence domain impedances will be exact for systems without MFC. In
the general case, they are expected to be more accurate than the decoupled
dqimpedance models. The paper defines a norm $\epsilon$ to measure the degree
of coupling in the impedance matrices. This norm equals the error in the
eigenvalue loci between the matrix and semidecoupled models. This can also be
viewed as the error in the semidecoupled Nyquist plot.
An example case study consisting of a gridconnected VSC with current
controller and PLL is used to compare the different methods. It is found that
decoupled and semidecoupled models in the dqdomain are only applicable in
grids with very low X/Rratio. Furthermore, it is concluded that the decoupled
model in the sequence domain gives close to equal results as the semidecoupled
model.

Based on the NagelSchreckenberg (NS) model with periodic boundary
conditions, we proposed the NSOS model by adding the overtaking strategy (OS).
In our model, overtaking vehicles are randomly selected with probability $q$ at
each time step, and the successful overtaking is determined by their
velocities. We observed that (i) traffic jams still occur in the NSOS model;
(ii) OS increases the traffic flow in the regime where the densities exceed the
maximum flow density. We also studied the phase transition (from free flow
phase to jammed phase) of the NSOS model by analyzing the overtaking success
rate, order parameter, relaxation time and correlation function, respectively.
It was shown that the NSOS model differs from the NS model mainly in the jammed
regime, and the influence of OS on the transition density is dominated by the
braking probability $p$

The NagelSchreckenberg model with overtaking strategy (NSOS) is proposed,
and numerical simulations are performed for both closed and open boundary
conditions. The fundamental diagram, spacetime diagram, and spatialtemporal
distribution of speed are investigated. In order to identify the synchronized
flow state, both the correlation functions (autocorrelation and
crosscorrelation) and the oneminute average flow rate vs. density diagram are
studied. All the results verify that synchronized flow does occur in our model.

Sensitivity analysis is concerned with understanding how the model output
depends on uncertainties (variances) in inputs and then identifies which inputs
are important in contributing to the prediction imprecision. Uncertainty
determination in output is the most crucial step in sensitivity analysis. In
the present paper, an analytic expression, which can exactly evaluate the
uncertainty in output as a function of the output's derivatives and inputs'
central moments, is firstly deduced for general multivariate models with given
relationship between output and inputs in terms of Taylor series expansion. A
$\gamma$order relative uncertainty for output, denoted by
$\mathrm{R^{\gamma}_v}$, is introduced to quantify the contributions of input
uncertainty of different orders. On this basis, it is shown that the widely
used approximation considering the first order contribution from the variance
of input variable can satisfactorily express the output uncertainty only when
the input variance is very small or the inputoutput function is almost linear.
Two applications of the analytic formula are performed to the power grid and
economic systems where the sensitivity of both actual power output and Economic
Order Quantity models are analyzed. The importance of each input variable in
response to the model outputs is quantified by the analytic formula.

It was proposed that the electric fields may lead to chiral separation in
QGP, which is called the chiral electric separation effect. The strong
electromagnetic field and the QCD vacuum can both completely be produced in the
offcentral nuclearnuclear collision. We used the WoodSaxon nucleon
distribution to calculate the electric field distributions of the offcentral
collisions. The chiral electro field spatial distribution at Relativistic
HeavyIon Collider (RHIC) and Large Hadron Collider (LHC) energy regions are
systematically studied in this paper. Compared with magnetic field spatial
distribution, electric field shows some different features in relativistic
heavyion collisions. The dependence of the electric field produced by the
thermal quark in the central position with different impact parameters on the
proper time with different collision energies in the RHIC and LHC energy region
are studied in this paper.

Representations of AC power systems by frequency dependent impedance
equivalents is an emerging technique in the dynamic analysis of power systems
including power electronic converters. The technique has been applied for
decades in DCpower systems, and it was recently adopted to map the impedances
in AC systems. Most of the work on AC systems can be categorized in two
approaches. One is the analysis of the system in the \textit{dq}domain,
whereas the other applies harmonic linearization in the phase domain through
symmetric components. Impedance models based on analytical calculations,
numerical simulation and experimental studies have been previously developed
and verified in both domains independently. The authors of previous studies
discuss the advantages and disadvantages of each domain separately, but neither
a rigorous comparison nor an attempt to bridge them has been conducted. The
present paper attempts to close this gap by deriving the mathematical
formulation that shows the equivalence between the \textit{dq}domain and the
sequence domain impedances. A modified form of the sequence domain impedance
matrix is proposed, and with this definition the stability estimates obtained
with the Generalized Nyquist Criterion (GNC) become equivalent in both domains.
The second contribution of the paper is the definition of a \textit{Mirror
Frequency Decoupled} (MFD) system. The analysis of MFD systems is less complex
than that of nonMFD systems because the positive and negative sequences are
decoupled. This paper shows that if a system is incorrectly assumed to be MFD,
this will lead to an erroneous or ambiguous estimation of the equivalent
impedance.

Many models and real complex systems possess critical thresholds at which the
systems shift from one sate to another. The discovery of the early warnings of
the systems in the vicinity of critical point are of great importance to
estimate how far a system is from a critical threshold. Multifractal Detrended
Fluctuation analysis (MFDFA) and visibility graph method have been employed to
investigate the fluctuation and geometrical structures of magnetization time
series of twodimensional Ising model around critical point. The Hurst exponent
has been confirmed to be a good indicator of phase transition. Increase of the
multifractality of the time series have been observed from generalized Hurst
exponents and singularity spectrum. Both Longterm correlation and broad
probability density function are identified to be the sources of
multifractality of time series near critical regime. Heterogeneous nature of
the networks constructed from magnetization time series have validated the
fractal properties of magnetization time series from complex network
perspective. Evolution of the topology quantities such as clustering
coefficient, average degree, average shortest path length, density,
assortativity and heterogeneity serve as early warnings of phase transition.
Those methods and results can provide new insights about analysis of phase
transition problems and can be used as early warnings for various complex
systems.

We employed the PYTHIA 6.4 model and the extended parton and hadron cascade
model PACIAE 2.2 to comparatively investigate the DIS normalized specific
charged hadron multiplicity in the 27.6 GeV electron semiinclusive
deepinelastic scattering off proton and deuteron. The PYTHIA and PACIAE
results calculated with default model parameters not well and fairly well
reproduce the corresponding HERMES data, respectively. In addition, we have
discussed the effects of the differences between the PYTHIA and PACIAE models.

Relativistic heavyion collisions can produce extremely strong magnetic field
in the collision regions. The spatial variation features of the magnetic fields
are analyzed in detail for noncentral Pb  Pb collisions at LHC
$\sqrt{s_{NN}}$= 900, 2760 and 7000 GeV and AuAu collisions at RHIC
$\sqrt{s_{NN}}$ = 62.4, 130 and 200 GeV. The dependencies of magnetic field on
proper time, collision energies and impact parameters are investigated in this
paper. It is shown that a enormous with highly inhomogeneous spatial
distribution magnetic field can indeed be created in offcentral relativistic
heavyion collisions in RHIC and LHC energy regions. The enormous magnetic
field is quite large, especially just after the collision, and then decreases
rapidly with time. We are surprised to find that the magnetic effect decreases
with increasing energy from RHIC to LHC energy region. It is found that the
magnitude of magnetic field in the LHC energy region is far less than that of
magnetic field in the RHIC energy region.

The parton and hadron cascade model PACIAE 2.1 (cf. Comput. Phys. Commun.184
(2013) 1476) has been upgraded to the new issue of PACIAE 2.2. By this new
issue the leptonnucleon and leptonnucleus (inclusive) deep inelastic
scatterings can also be investigated. As an example, the PACIAE 2.2 model is
enabled to calculate the specific charged hadron multiplicity in the $e^$+p
and $e^$+D semiinclusive deepinelastic scattering at 27.6 GeV electron beam
energy. The calculated results are well comparing with the corresponding HERMES
data. Additionally, the effect of model parameters \alpha and \beta in the Lund
string fragmentation function on the multiplicity is studied.

Based on the requirement in the simulation of leptonnucleus deep inelastic
scattering (DIS), we construct a fortran program LDCS 1.0 calculating the
differential and total cross sections for the unpolarized charged
leptonunpolarized nucleon and neutrinounpolarized nucleon neutral current
(charged current) DIS at leading order. Any set of the experimentally fitted
parton distribution functions could be employed directly. The mass of incident
and scattered leptons is taken into account and the boundary conditions
calculating the single differential and total cross section are studied. The
calculated results well agree with the corresponding experimental data which
indicating the LDCS 1.0 program is good. It is also turned out that the effect
of tauon mass is not negligible in the GeV energy level.

The features of magnetic field in relativistic heavyion collisions are
systematically studied by using a modified magnetic field model in this paper.
The features of magnetic field distributions in the central point are studied
in the RHIC and LHC energy regions. We also predict the feature of magnetic
fields at LHC $\sqrt{s_{NN}}$= 900, 2760 and 7000 GeV based on the detailed
study at RHIC $\sqrt{s_{NN}}$ = 62.4, 130 and 200 GeV. The dependencies of the
features of magnetic fields on the collision energies, centralities and
collision time are systematically investigated, respectively.

We employed the new issue of a parton and hadron cascade model PACIAE 2.1 to
systematically investigate the charged particle elliptic flow parameter $v_2$
in the relativistic nuclear collisions at RHIC and LHC energies. With randomly
sampling the transverse momentum $x$ and $y$ components of the particles
generated in string fragmentation on the circumference of an ellipse instead of
circle originally, the calculated charged particle $v_2(\eta)$ and $v_2(p_T)$
fairly reproduce the corresponding experimental data in the Au+Au/Pb+Pb
collisions at $\sqrt{s_{NN}}$=0.2/2.76 TeV. In addition, the charged particle
$v_2(\eta)$ and $v_2(p_T)$ in the p+p collisions at $\sqrt s$=7 TeV as well as
in the p+Au/p+Pb collisions at $\sqrt{s_{NN}}$=0.2/5.02 TeV are predicted.

Metro networks serve as good examples of traffic systems for understanding
the relations between geometric structures and transport properties.We study
and compare 28 world major metro networks in terms of the Wasserstein distance,
the key metric for optimal transport, and measures geometry related, e.g.
fractal dimension, graph energy and graph spectral distance. The finding of
powerlaw relationships between rescaled graph energy and fractal dimension for
both unweighted and weighted metro networks indicates the energy costs per unit
area are lower for higher dimensioned metros. In L space, the mean Wasserstein
distance between any pair of connected stations is proportional to the fractal
dimension, which is in the vicinity of our theoretical calculations treated on
special regular tree graphs. This finding reveals the geometry of metro
networks and tree graphs are in close proximity to one another. In P space, the
mean Wasserstein distance between any pair of stations relates closely to the
average number of transfers. By ranking several key quantities transport
concerned, we obtain several ranking lists in which New York metro and Berlin
metro consistently top the first two spots.

We have updated the parton and hadron cascade model PACIAE 2.0 for
relativistic nuclear collisions (cf. Comput. Phys. Comm. 183(2012)333) to the
new issue of PACIAE 2.1. The PACIAE model is based on PYTHIA. In the PYTHIA
model, once the generated particle/parton transverse momentum $p_T$ is randomly
selected, its component of $[p_x,p_y]$ was originally put on the circle with
radius $p_T$ randomly. It is now replaced by put on the circumference of
ellipse with half major and minor axes of $p_T(1+\delta_p)$ and
$p_T(1\delta_p)$, respectively, in order to investigate better the transverse
momentum anisotropy in final hadronic state.

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 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$.

Ranking is a ubiquitous phenomenon in the human society. By clicking the web
pages of Forbes, you may find all kinds of rankings, such as world's most
powerful people, world's richest people, toppaid tennis stars, and so on and
so forth. Herewith, we study a specific kind, sports ranking systems in which
players' scores and prize money are calculated based on their performances in
attending various tournaments. A typical example is tennis. It is found that
the distributions of both scores and prize money follow universal power laws,
with exponents nearly identical for most sports fields. In order to understand
the origin of this universal scaling we focus on the tennis ranking systems. By
checking the data we find that, for any pair of players, the probability that
the higherranked player will top the lowerranked opponent is proportional to
the rank difference between the pair. Such a dependence can be well fitted to a
sigmoidal function. By using this feature, we propose a simple toy model which
can simulate the competition of players in different tournaments. The
simulations yield results consistent with the empirical findings. Extensive
studies indicate the model is robust with respect to the modifications of the
minor parts.