• ### Analytic sensitivity analysis for models with correlated input variables(1708.06748)

Aug. 11, 2017 physics.soc-ph
An analytic formula is proposed to characterize the variance propagation from correlated input variables to the model response, by using multi-variate 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.
• ### Harmonic State Space Modeling of a Three-Phase Modular Multilevel Converter(1706.09925)

June 29, 2017 cs.SY
This paper presents the harmonic state space (HSS) modeling of a three-phase modular multilevel converter (MMC). MMC is a converter system with a typical multi-frequency 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 two-level voltage-source converters (VSCs), where only the fundamental-frequency 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 state-space form, is introduced in this paper to model the MMC in order to capture all the harmonics and the frequency couplings. The steady-state and small-signal dynamic HSS models of a three-phase MMC are developed, respectively. The validity of the developed HSS model of a three-phase 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.
• ### A simple method for shifting local dq impedance models to a global reference frame for stability analysis(1706.08313)

June 26, 2017 cs.SY
Impedance-based stability analysis in the dq-domain 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 sub-models 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.
• ### Impedance Analysis of Modular Multilevel Converter Based on Harmonic State-Space Modeling Method(1705.01030)

May 2, 2017 cs.SY
The small-signal impedance modeling of modular multilevel converter (MMC) is the key for analyzing resonance and stability of MMC-based ac power electronics systems. MMC is a converter system with a typical multi-frequency response due to its significant steady-state harmonic components in the arm currents, capacitor voltages, and control signals. Therefore, traditional small-signal modeling methods for 2-level voltage-source converters (VSCs) cannot be directly applied to the MMC. In this paper, the harmonic state-space (HSS) modeling approach is introduced to characterize the harmonic coupling behavior of the MMC. On this basis, the small-signal impedance models of the MMC are developed according to the harmonic linearization principle, which can include all the steady-state harmonic effects of the state variables, leading to the accurate impedance models. Furthermore, in order to reveal the impact of the internal dynamics and closed-loop control on the small-signal impedance of the MMC, three cases are considered in this paper, i.e., open-loop control, ac voltage closed-loop control, and circulating current closed-loop control. Finally, the analytical impedance models are verified by both simulation and experimental results.
• ### Accurate and reduced SISO Sequence Impedance Models of Grid-tied Voltage Source Converter for Small Signal Stability Analysis(1704.04157)

April 13, 2017 cs.SY
Impedance models are widely used in assessing small signal stability of grid-tied voltage source converter (VSC) systems. Recent research has proven that impedance models of grid-tied VSC in both dq and sequence domains are generally Multi-Input Multi-Output (MIMO) systems, and the generalized Nyquist criterion has to be applied for stability analysis to these MIMO systems. However, finding Single-Input and Single-Output (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 grid-tied VSC system, one is derived from the strong grid assumption, and the other is from the closed-loop 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 (matrix-based) 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 net-baryon densities. There is evidence that a Quark-Gluon-Plasma (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 net-baryon densities and lower temperatures, it is expected that the QCD phase diagram exhibits a rich structure, such as a first-order 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 high-energy heavy-ion 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 net-baryon densities, because it is designed to run at unprecedented interaction rates. High-rate operation is the key prerequisite for high-precision measurements of multi-differential 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 baryon-chemical potentials (mu_B > 500 MeV), effects of chiral symmetry, and the equation-of-state 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 high-density QCD matter.
• ### Augmented Sequence Impedance Networks of Grid-tied Voltage Source Converter for Stability Analysis(1703.03977)

March 11, 2017 cs.SY
Impedance-based stability analysis is appealing in the case of Single-Input Single-Output (SISO) systems. However, in the case of grid-tied voltage source converter (VSC) systems, dq impedances of source and load (VSC) subsystems are typically Multi-Input Multi-Output (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 non-coupled 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.
• ### Coupled and decoupled impedance models compared in power electronics systems(1610.04988)

Oct. 17, 2016 cs.SY
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 Single-Input-Single-Output (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 semi-decoupled 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 dq-impedance 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 semi-decoupled models. This can also be viewed as the error in the semi-decoupled Nyquist plot. An example case study consisting of a grid-connected VSC with current controller and PLL is used to compare the different methods. It is found that decoupled and semi-decoupled models in the dq-domain are only applicable in grids with very low X/R-ratio. Furthermore, it is concluded that the decoupled model in the sequence domain gives close to equal results as the semi-decoupled model.
• ### The effects of overtaking strategy in the Nagel-Schreckenberg model(1506.04990)

Sept. 8, 2016 nlin.CG
Based on the Nagel-Schreckenberg (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$
• ### Occurrence of synchronized flow due to overtaking strategy in the Nagel-Schreckenberg model(1608.07948)

Aug. 29, 2016 nlin.CG
The Nagel-Schreckenberg model with overtaking strategy (NSOS) is proposed, and numerical simulations are performed for both closed and open boundary conditions. The fundamental diagram, space-time diagram, and spatial-temporal distribution of speed are investigated. In order to identify the synchronized flow state, both the correlation functions (autocorrelation and cross-correlation) and the one-minute average flow rate vs. density diagram are studied. All the results verify that synchronized flow does occur in our model.
• ### An analytic method for sensitivity analysis of complex systems(1608.07794)

Aug. 28, 2016 physics.soc-ph
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 input-output 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.
• ### Chiral electric field in relativistic heavy-ion collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider(1507.07743)

June 2, 2016 hep-ph, nucl-th
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 off-central nuclear-nuclear collision. We used the Wood-Saxon nucleon distribution to calculate the electric field distributions of the off-central collisions. The chiral electro field spatial distribution at Relativistic Heavy-Ion 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 heavy-ion 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.
• ### A modified sequence domain impedance definition and its equivalence to the dq-domain impedance definition for the stability analysis of AC power electronic systems(1605.00526)

April 28, 2016 cs.SY
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 DC-power 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 non-MFD 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.
• ### Multifractal and Network Analysis of Phase Transition(1601.07715)

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 (MF-DFA) and visibility graph method have been employed to investigate the fluctuation and geometrical structures of magnetization time series of two-dimensional 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 Long-term 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.
• ### The specific charged hadron multiplicity in ep and eD semi-inclusive deep-inelastic scattering in the PYTHIA and PACIAE models(1511.04844)

Dec. 25, 2015 nucl-th
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 semi-inclusive deep-inelastic 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.
• ### The spatial distributions of chiral magnetic field in the RHIC and LHC energy regions(1410.6349)

March 21, 2015 hep-ph, nucl-th
Relativistic heavy-ion collisions can produce extremely strong magnetic field in the collision regions. The spatial variation features of the magnetic fields are analyzed in detail for non-central Pb - Pb collisions at LHC $\sqrt{s_{NN}}$= 900, 2760 and 7000 GeV and Au-Au 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 off-central relativistic heavy-ion 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.

Dec. 24, 2014 hep-ph, nucl-th
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 lepton-nucleon and lepton-nucleus (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 semi-inclusive deep-inelastic 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.
• ### On the lepton-nucleon neutral and charged current deep inelastic scattering cross sections(1412.5268)

Dec. 17, 2014 hep-ph, nucl-th
Based on the requirement in the simulation of lepton-nucleus deep inelastic scattering (DIS), we construct a fortran program LDCS 1.0 calculating the differential and total cross sections for the unpolarized charged lepton-unpolarized nucleon and neutrino-unpolarized 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.
• ### A systematic study of magnetic field in Relativistic Heavy-ion Collisions in the RHIC and LHC energy regions(1408.5694)

Aug. 25, 2014 hep-ph, nucl-th
The features of magnetic field in relativistic heavy-ion 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.
• ### Systematic study of elliptic flow parameter in the relativistic nuclear collisions at RHIC and LHC energies(1307.1241)

April 17, 2014 nucl-th
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.
• ### Optimal Transport in Worldwide Metro Networks(1403.7844)

April 14, 2014 physics.soc-ph
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 power-law 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.
• ### PACIAE 2.1: An updated issue of the parton and hadron cascade model PACIAE 2.0(1206.4795)

Dec. 23, 2012 hep-ph, nucl-th
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.
• ### Comparative study for non-statistical fluctuation of net- proton,baryon, and charge multiplicities(1212.2283)

Dec. 11, 2012 nucl-th
We calculate the real and non-statistical higher moment excitation functions ($\sqrt{s_{NN}}$ = 11.5 to 200 GeV) for the net-proton, net-baryon, and the net-charge 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 non-statistical higher moments. But for a given widow size, the CP singularity may show only in the non-statistical higher moment excitation functions of a definite conserved observable.
• ### Higher moment singularities explored by the net proton non-statistical fluctuations(1205.5634)

May 25, 2012 nucl-th
We use the non-statistical 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. Non-statistical 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 non-statistical $S$ and $S\sigma$.
• ### Universal scaling in sports ranking(1111.2919)

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, top-paid 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 higher-ranked player will top the lower-ranked 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.