• ### Scale-free Resilience of Real Traffic Jams(1804.11047)

April 30, 2018 physics.soc-ph, nlin.AO
The concept of resilience can be realized in daily traffic, representing the ability of transportation system to adapt and recover from traffic jams. Although resilience is a critical property needed for understanding and managing the organization balance in traffic networks, a systematic study of resilience as well as its accepted definition in real urban traffic is still missing. Here we define a city traffic resilience as the size of the spatio-temporal clusters of congestion in real traffic, and find that the resilience follows a scale free distribution in both, city road networks and in highways, with different exponents, but similar exponents in different days and cities. The traffic resilience is also revealed to have a novel scaling relation between the cluster size of the spatio-temporal jam and its recovery duration, independent of microscopic details. Our findings provide insight towards universal properties of traffic resilience for better forecasting and mitigation of congestions.
• ### Scale-free Resilience of Real Traffic Jams(1804.11047)

April 30, 2018 physics.soc-ph, nlin.AO
The concept of resilience can be realized in daily traffic, representing the ability of transportation system to adapt and recover from traffic jams. Although resilience is a critical property needed for understanding and managing the organization balance in traffic networks, a systematic study of resilience as well as its accepted definition in real urban traffic is still missing. Here we define a city traffic resilience as the size of the spatio-temporal clusters of congestion in real traffic, and find that the resilience follows a scale free distribution in both, city road networks and in highways, with different exponents, but similar exponents in different days and cities. The traffic resilience is also revealed to have a novel scaling relation between the cluster size of the spatio-temporal jam and its recovery duration, independent of microscopic details. Our findings provide insight towards universal properties of traffic resilience for better forecasting and mitigation of congestions.
• ### Scale-free Resilience of Real Traffic Jams(1804.11047)

April 30, 2018 physics.soc-ph, nlin.AO
The concept of resilience can be realized in daily traffic, representing the ability of transportation system to adapt and recover from traffic jams. Although resilience is a critical property needed for understanding and managing the organization balance in traffic networks, a systematic study of resilience as well as its accepted definition in real urban traffic is still missing. Here we define a city traffic resilience as the size of the spatio-temporal clusters of congestion in real traffic, and find that the resilience follows a scale free distribution in both, city road networks and in highways, with different exponents, but similar exponents in different days and cities. The traffic resilience is also revealed to have a novel scaling relation between the cluster size of the spatio-temporal jam and its recovery duration, independent of microscopic details. Our findings provide insight towards universal properties of traffic resilience for better forecasting and mitigation of congestions.
• ### Response of gadolinium doped liquid scintillator to charged particles: measurement based on intrinsic U/Th contamination(1801.03620)

March 6, 2018 hep-ex, physics.ins-det
A measurement is reported for the response to charged particles of a liquid scintillator named EJ-335 doped with 0.5% gadolinium by weight. This liquid scintillator was used as the detection medium in a neutron detector. The measurement is based on the in-situ $\alpha$-particles from the intrinsic Uranium and Thorium contamination in the scintillator. The $\beta$-$\alpha$ and the $\alpha$-$\alpha$ cascade decays from the U/Th decay chains were used to select $\alpha$-particles. The contamination levels of U/Th were consequently measured to be $(5.54\pm0.15)\times 10^{-11}$ g/g, $(1.45\pm0.01)\times 10^{-10}$ g/g and $(1.07\pm0.01)\times 10^{-11}$ g/g for $^{232}$Th, $^{238}$U and $^{235}$U, respectively, assuming secular equilibrium. The stopping power of $\alpha$-particles in the liquid scintillator was simulated by the TRIM software. Then the Birks constant, $kB$, of the scintillator for $\alpha$-particles was determined to be $(7.28\pm0.23)$ mg/(cm$^{2}\cdot$MeV) by Birks' formulation. The response for protons is also presented assuming the $kB$ constant is the same as for $\alpha$-particles.
• ### Artifacts in magnetic spirals retrieved by transport of intensity equation (TIE)(1709.07160)

Dec. 18, 2017 cond-mat.mtrl-sci
The artifacts in the magnetic structures reconstructed from Lorentz transmission electron microscopy (LTEM) images with TIE method have been analyzed in detail. The processing for the simulated images of Bloch and Neel spirals indicated that the improper parameters in TIE may overestimate the high frequency information and induce some false features in the retrieved images. The specimen tilting will further complicate the analysis of the images because the LTEM image contrast is not the result of the magnetization distribution within the specimen but the integral projection pattern of the magnetic induction filling the entire space including the specimen.
• ### Spin transport in two-layer-CVD-hBN/graphene/hBN heterostructures(1712.00815)

We study room temperature spin transport in graphene devices encapsulated between a layer-by-layer-stacked two-layer-thick chemical vapour deposition (CVD) grown hexagonal boron nitride (hBN) tunnel barrier, and a few-layer-thick exfoliated-hBN substrate. We find mobilities and spin-relaxation times comparable to that of SiO$_2$ substrate based graphene devices, and obtain a similar order of magnitude of spin relaxation rates for both the Elliott-Yafet and D'Yakonov-Perel' mechanisms. The behaviour of ferromagnet/two-layer-CVD-hBN/graphene/hBN contacts ranges from transparent to tunneling due to inhomogeneities in the CVD-hBN barriers. Surprisingly, we find both positive and negative spin polarizations for high-resistance two-layer-CVD-hBN barrier contacts with respect to the low-resistance contacts. Furthermore, we find that the differential spin injection polarization of the high-resistance contacts can be modulated by DC bias from -0.3 V to +0.3 V with no change in its sign, while its magnitude increases at higher negative bias. These features mark a distinctive spin injection nature of the two-layer-CVD-hBN compared to the bilayer-exfoliated-hBN tunnel barriers.
• ### Measurement of the fast neutron background at the China Jinping Underground Laboratory(1706.06831)

Oct. 17, 2017 hep-ex, nucl-ex, physics.ins-det
We report on the measurements of the fluxes and spectra of the environmental fast neutron background at the China Jinping Underground Laboratory (CJPL) with a rock overburden of about 6700 meters water equivalent, using a liquid scintillator detector doped with 0.5% gadolinium. The signature of a prompt nuclear recoil followed by a delayed high energy $\gamma$-ray cascade is used to identify neutron events. The large energy deposition of the delayed $\gamma$-rays from the $(n, \gamma)$ reaction on gadolinium, together with the excellent n-$\gamma$ discrimination capability provides a powerful background suppression which allows the measurement of a low intensity neutron flux. The neutron flux of $(1.51\pm0.03(stat.)\pm0.10(syst.))\times10^{-7}$ cm$^{-2}$s$^{-1}$ in the energy range of 1 -- 10 MeV in the Hall A of CJPL was measured based on 356 days of data. In the same energy region, measurement with the same detector placed in a one meter thick polyethylene room gives a significantly lower flux of $(4.9\pm0.9(stat.)\pm0.5(syst.))\times10^{-9}$ cm$^{-2}$s$^{-1}$ with 174 days of data. This represents a measurement of the lowest environmental fast neutron background among the underground laboratories in the world, prior to additional experiment-specific attenuation. Additionally, the fast neutron spectra both in the Hall A and the polyethylene room were reconstructed with the help of GEANT4 simulation.
• ### Constraining the Mass and Radius of Neutron Stars in Globular Clusters(1709.05013)

Sept. 15, 2017 nucl-th, astro-ph.SR, astro-ph.HE
We analyze observations of eight quiescent low-mass X-ray binaries in globular clusters and combine them to determine the neutron star mass-radius curve and the equation of state of dense matter. We determine the effect that several uncertainties may have on our results, including uncertainties in the distance, the atmosphere composition, the neutron star maximum mass, the neutron star mass distribution, the possible presence of a hotspot on the neutron star surface, and the prior choice for the equation of state of dense matter. We find that the radius of a 1.4 solar mass neutron star is most likely from 10 to 14 km and that tighter constraints are only possible with stronger assumptions about the nature of the neutron stars, the systematics of the observations, or the nature of dense matter. Strong phase transitions are preferred over other models and interpretations of the data with a Bayes factor of 8 or more, and in this case, the radius is likely smaller than 12 km. However, radii larger than 12 km are preferred if the neutron stars have uneven temperature distributions.
• ### T-DB: Toward Fully Functional Transparent Encrypted Databases in DBaaS Framework(1708.08191)

Aug. 28, 2017 cs.CR, cs.DB
Individuals and organizations tend to migrate their data to clouds, especially in a DataBase as a Service (DBaaS) pattern. The major obstacle is the conflict between secrecy and utilization of the relational database to be outsourced. We address this obstacle with a Transparent DataBase (T-DB) system strictly following the unmodified DBaaS framework. A database owner outsources an encrypted database to a cloud platform, needing only to store the secret keys for encryption and an empty table header for the database; the database users can make almost all types of queries on the encrypted database as usual; and the cloud can process ciphertext queries as if the database were not encrypted. Experimentations in realistic cloud environments demonstrate that T-DB has perfect query answer precision and outstanding performance.
• ### T-DB: Toward Fully Functional Transparent Encrypted Databases in DBaaS Framework(1708.08191)

Aug. 28, 2017 cs.CR, cs.DB
Individuals and organizations tend to migrate their data to clouds, especially in a DataBase as a Service (DBaaS) pattern. The major obstacle is the conflict between secrecy and utilization of the relational database to be outsourced. We address this obstacle with a Transparent DataBase (T-DB) system strictly following the unmodified DBaaS framework. A database owner outsources an encrypted database to a cloud platform, needing only to store the secret keys for encryption and an empty table header for the database; the database users can make almost all types of queries on the encrypted database as usual; and the cloud can process ciphertext queries as if the database were not encrypted. Experimentations in realistic cloud environments demonstrate that T-DB has perfect query answer precision and outstanding performance.
• ### Investigation of the Relationship Between the Air Pollution and Solar Activity(1706.06864)

June 20, 2017 physics.ao-ph
How did the Sun affect the air pollution on the Earth? There are few papers about this question. This work investigates the relationship between the air pollution and solar activity by using the geophysical and environmental data during the period of 2000-2016. It is quite certain that the solar activity may impact on the air pollution, but the relationship is very weak and indirect. The Pearson correlation, Spearman rank correlation, Kendalls rank correlation, and conditional probability were adopted to analyze the air pollution index (API), air quality index (AQI), sunspot number (SSN), radio flux at wavelength of 10.7 cm (F10.7), and total solar irradiance (TSI). The analysis implies that the correlation coefficient between API and SSN is weak ($-0.17<r<0.32$) with complex variation. The main results are: (1) For cities with higher air pollution, the probability of high API will be increased along with SSN, then reach to a maximum, and then decreased; (2) For cities with lower air pollution, the API has lower correlation with SSN; (3) The relationship between API and F10.7, or API and TSI are also similar as API and SSN. The solar activities take direct effect on TSI and the energetic particle flux, and indirect and long-term effect on lower atmosphere and weather near the Earth. All of these factors contribute to the air pollution on the Earth.
• ### Investigation of the Relationship Between the Air Pollution and Solar Activity(1706.06864)

June 20, 2017 physics.ao-ph
How did the Sun affect the air pollution on the Earth? There are few papers about this question. This work investigates the relationship between the air pollution and solar activity by using the geophysical and environmental data during the period of 2000-2016. It is quite certain that the solar activity may impact on the air pollution, but the relationship is very weak and indirect. The Pearson correlation, Spearman rank correlation, Kendalls rank correlation, and conditional probability were adopted to analyze the air pollution index (API), air quality index (AQI), sunspot number (SSN), radio flux at wavelength of 10.7 cm (F10.7), and total solar irradiance (TSI). The analysis implies that the correlation coefficient between API and SSN is weak ($-0.17<r<0.32$) with complex variation. The main results are: (1) For cities with higher air pollution, the probability of high API will be increased along with SSN, then reach to a maximum, and then decreased; (2) For cities with lower air pollution, the API has lower correlation with SSN; (3) The relationship between API and F10.7, or API and TSI are also similar as API and SSN. The solar activities take direct effect on TSI and the energetic particle flux, and indirect and long-term effect on lower atmosphere and weather near the Earth. All of these factors contribute to the air pollution on the Earth.
• ### Investigation of the Relationship Between the Air Pollution and Solar Activity(1706.06864)

June 20, 2017 physics.ao-ph
How did the Sun affect the air pollution on the Earth? There are few papers about this question. This work investigates the relationship between the air pollution and solar activity by using the geophysical and environmental data during the period of 2000-2016. It is quite certain that the solar activity may impact on the air pollution, but the relationship is very weak and indirect. The Pearson correlation, Spearman rank correlation, Kendalls rank correlation, and conditional probability were adopted to analyze the air pollution index (API), air quality index (AQI), sunspot number (SSN), radio flux at wavelength of 10.7 cm (F10.7), and total solar irradiance (TSI). The analysis implies that the correlation coefficient between API and SSN is weak ($-0.17<r<0.32$) with complex variation. The main results are: (1) For cities with higher air pollution, the probability of high API will be increased along with SSN, then reach to a maximum, and then decreased; (2) For cities with lower air pollution, the API has lower correlation with SSN; (3) The relationship between API and F10.7, or API and TSI are also similar as API and SSN. The solar activities take direct effect on TSI and the energetic particle flux, and indirect and long-term effect on lower atmosphere and weather near the Earth. All of these factors contribute to the air pollution on the Earth.
• ### Investigation of the Relationship Between the Air Pollution and Solar Activity(1706.06864)

June 20, 2017 physics.ao-ph
How did the Sun affect the air pollution on the Earth? There are few papers about this question. This work investigates the relationship between the air pollution and solar activity by using the geophysical and environmental data during the period of 2000-2016. It is quite certain that the solar activity may impact on the air pollution, but the relationship is very weak and indirect. The Pearson correlation, Spearman rank correlation, Kendalls rank correlation, and conditional probability were adopted to analyze the air pollution index (API), air quality index (AQI), sunspot number (SSN), radio flux at wavelength of 10.7 cm (F10.7), and total solar irradiance (TSI). The analysis implies that the correlation coefficient between API and SSN is weak ($-0.17<r<0.32$) with complex variation. The main results are: (1) For cities with higher air pollution, the probability of high API will be increased along with SSN, then reach to a maximum, and then decreased; (2) For cities with lower air pollution, the API has lower correlation with SSN; (3) The relationship between API and F10.7, or API and TSI are also similar as API and SSN. The solar activities take direct effect on TSI and the energetic particle flux, and indirect and long-term effect on lower atmosphere and weather near the Earth. All of these factors contribute to the air pollution on the Earth.
• ### Modelling the Type Ic SN 2004aw: a Moderately Energetic Explosion of a Massive C+O Star without a GRB(1705.10249)

May 29, 2017 astro-ph.HE
An analysis of the Type Ic supernova (SN) 2004aw is performed by means of models of the photospheric and nebular spectra and of the bolometric light curve. SN2004aw is shown not to be broad-lined'', contrary to previous claims, but rather a fast-lined'' SN Ic. The spectral resemblance to the narrow-lined Type Ic SN1994I, combined with the strong nebular [O I] emission and the broad light curve, point to a moderately energetic explosion of a massive C+O star. The ejected 56Ni mass is ~0.2 Msun. The ejecta mass as constrained by the models is ~3-5 Msun, while the kinetic energy is estimated as KE ~3-6 e51 ergs. The ratio KE/Mej, the specific energy which influences the shape of the spectrum, is therefore ~1. The corresponding zero-age main-sequence mass of the progenitor star may have been ~23-28 Msun. Tests show that a flatter outer density structure may have caused a broad-lined spectrum at epochs before those observed without affecting the later epochs when data are available, implying that our estimate of KE is a lower limit. SN2004aw may have been powered by either a collapsar or a magnetar, both of which have been proposed for gamma-ray burst-supernovae. Evidence for this is seen in the innermost layers, which appear to be highly aspherical as suggested by the nebular line profiles. However, any engine was not extremely powerful, as the outer ejecta are more consistent with a spherical explosion and no gamma-ray burst was detected in coincidence with SN2004aw.
• We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($\tau^{+}\tau^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.
• ### Speeding up $N$-body simulations of modified gravity: Chameleon screening models(1611.09375)

March 12, 2017 astro-ph.CO
We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of $f(R)$ simulations. For example, a test simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is now 5 times faster than before, while a Millennium-resolution simulation for $f(R)$ gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.
• The first observation of the rare decay$B^0_s \to \phi\pi^+\pi^-$ and evidence for $B^0 \to \phi\pi^+\pi^-$ are reported, using $pp$ collision data recorded by the LHCb detector at centre-of-mass energies $\sqrt{s} = 7$ and 8~TeV, corresponding to an integrated luminosity of $3{\mbox{\,fb}^{-1}}$. The branching fractions in the $\pi^+\pi^-$ invariant mass range $400<m(\pi^+\pi^-)<1600{\mathrm{\,Me\kern -0.1em V\!/}c^2}$ are $[3.48\pm 0.23\pm 0.17\pm 0.35]\times 10^{-6}$ and $[1.82\pm 0.25\pm 0.41\pm 0.14]\times 10^{-7}$ for $B^0_s \to \phi\pi^+\pi^-$ and $B^0 \to \phi\pi^+\pi^-$ respectively, where the uncertainties are statistical, systematic and from the normalisation mode $B^0_s \to \phi\phi$. A combined analysis of the $\pi^+\pi^-$ mass spectrum and the decay angles of the final-state particles identifies the exclusive decays $B^0_s \to \phi f_0(980)$, $B_s^0 \to \phi f_2(1270)$, and $B^0_s \to \phi\rho^0$ with branching fractions of $[1.12\pm 0.16^{+0.09}_{-0.08}\pm 0.11]\times 10^{-6}$, $[0.61\pm 0.13^{+0.12}_{-0.05}\pm 0.06]\times 10^{-6}$ and $[2.7\pm 0.7\pm 0.2\pm 0.2]\times 10^{-7}$, respectively.
• The first observation of the rare decay$B^0_s \to \phi\pi^+\pi^-$ and evidence for $B^0 \to \phi\pi^+\pi^-$ are reported, using $pp$ collision data recorded by the LHCb detector at centre-of-mass energies $\sqrt{s} = 7$ and 8~TeV, corresponding to an integrated luminosity of $3{\mbox{\,fb}^{-1}}$. The branching fractions in the $\pi^+\pi^-$ invariant mass range $400<m(\pi^+\pi^-)<1600{\mathrm{\,Me\kern -0.1em V\!/}c^2}$ are $[3.48\pm 0.23\pm 0.17\pm 0.35]\times 10^{-6}$ and $[1.82\pm 0.25\pm 0.41\pm 0.14]\times 10^{-7}$ for $B^0_s \to \phi\pi^+\pi^-$ and $B^0 \to \phi\pi^+\pi^-$ respectively, where the uncertainties are statistical, systematic and from the normalisation mode $B^0_s \to \phi\phi$. A combined analysis of the $\pi^+\pi^-$ mass spectrum and the decay angles of the final-state particles identifies the exclusive decays $B^0_s \to \phi f_0(980)$, $B_s^0 \to \phi f_2(1270)$, and $B^0_s \to \phi\rho^0$ with branching fractions of $[1.12\pm 0.16^{+0.09}_{-0.08}\pm 0.11]\times 10^{-6}$, $[0.61\pm 0.13^{+0.12}_{-0.05}\pm 0.06]\times 10^{-6}$ and $[2.7\pm 0.7\pm 0.2\pm 0.2]\times 10^{-7}$, respectively.
• We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.
• ### Phosphorylation potential and chemical fluxes govern the biological performance of multiple PdP cycles(1612.01233)

Dec. 5, 2016 q-bio.MN
Fission yeast G2/M transition is regulated by a biochemical reaction networks which contains four components: Cdc13, Cdc2, Wee1, and Cdc25. This circuit is characterized by the ultrasensitive responses of Wee1 or Cdc25 to Cdc13/Cdc2 activity, and the bistability of Cdc2 activation. Previous work has shown that this bistability is governed by phosphorylation energy. In this article, we developed the kinetic model of this circuit and conducted further thermodynamic analysis on the role of phosphorylation energy (&[Delta]G). We showed that level &[Delta]G shapes the response curves of Wee1 or Cdc25 to Cdc2 and governs the intrinsic noise level of Cdc2 activity. More importantly, the mutually antagonistic chemical fluxes around the PdP cycles in G2/M circuit were shown to act as a stabilizer of Cdc2 activity against &[Delta]G fluctuations. These results suggests the fundamental role of free energy and chemical fluxes on the sensitivity, bistability and robustness of G2/M transition.
• ### Phosphorylation potential and chemical fluxes govern the biological performance of multiple PdP cycles(1612.01233)

Dec. 5, 2016 q-bio.MN
Fission yeast G2/M transition is regulated by a biochemical reaction networks which contains four components: Cdc13, Cdc2, Wee1, and Cdc25. This circuit is characterized by the ultrasensitive responses of Wee1 or Cdc25 to Cdc13/Cdc2 activity, and the bistability of Cdc2 activation. Previous work has shown that this bistability is governed by phosphorylation energy. In this article, we developed the kinetic model of this circuit and conducted further thermodynamic analysis on the role of phosphorylation energy (&[Delta]G). We showed that level &[Delta]G shapes the response curves of Wee1 or Cdc25 to Cdc2 and governs the intrinsic noise level of Cdc2 activity. More importantly, the mutually antagonistic chemical fluxes around the PdP cycles in G2/M circuit were shown to act as a stabilizer of Cdc2 activity against &[Delta]G fluctuations. These results suggests the fundamental role of free energy and chemical fluxes on the sensitivity, bistability and robustness of G2/M transition.
• ### Phosphorylation potential and chemical fluxes govern the biological performance of multiple PdP cycles(1612.01233)

Dec. 5, 2016 q-bio.MN
Fission yeast G2/M transition is regulated by a biochemical reaction networks which contains four components: Cdc13, Cdc2, Wee1, and Cdc25. This circuit is characterized by the ultrasensitive responses of Wee1 or Cdc25 to Cdc13/Cdc2 activity, and the bistability of Cdc2 activation. Previous work has shown that this bistability is governed by phosphorylation energy. In this article, we developed the kinetic model of this circuit and conducted further thermodynamic analysis on the role of phosphorylation energy (&[Delta]G). We showed that level &[Delta]G shapes the response curves of Wee1 or Cdc25 to Cdc2 and governs the intrinsic noise level of Cdc2 activity. More importantly, the mutually antagonistic chemical fluxes around the PdP cycles in G2/M circuit were shown to act as a stabilizer of Cdc2 activity against &[Delta]G fluctuations. These results suggests the fundamental role of free energy and chemical fluxes on the sensitivity, bistability and robustness of G2/M transition.
• ### Resilience of epidemics on networks(1610.06064)

Epidemic propagation on complex networks has been widely investigated, mostly with invariant parameters. However, the process of epidemic propagation is not always constant. Epidemics can be affected by various perturbations, and may bounce back to its original state, which is considered resilient. Here, we study the resilience of epidemics on networks, by introducing a different infection rate ${\lambda_{2}}$ during SIS (susceptible-infected-susceptible) epidemic propagation to model perturbations (control state), whereas the infection rate is ${\lambda_{1}}$ in the rest of time. Through simulations and theoretical analysis, we find that even for ${\lambda_{2}<\lambda_{c}}$, epidemics eventually could bounce back if control duration is below a threshold. This critical control time for epidemic resilience, i.e., ${cd_{max}}$ can be predicted by the diameter (${d}$) of the underlying network, with the quantitative relation ${cd_{max}\sim d^{\alpha}}$. Our findings can help to design a better mitigation strategy for epidemics.
• ### Resilience of epidemics on networks(1610.06064)

Epidemic propagation on complex networks has been widely investigated, mostly with invariant parameters. However, the process of epidemic propagation is not always constant. Epidemics can be affected by various perturbations, and may bounce back to its original state, which is considered resilient. Here, we study the resilience of epidemics on networks, by introducing a different infection rate ${\lambda_{2}}$ during SIS (susceptible-infected-susceptible) epidemic propagation to model perturbations (control state), whereas the infection rate is ${\lambda_{1}}$ in the rest of time. Through simulations and theoretical analysis, we find that even for ${\lambda_{2}<\lambda_{c}}$, epidemics eventually could bounce back if control duration is below a threshold. This critical control time for epidemic resilience, i.e., ${cd_{max}}$ can be predicted by the diameter (${d}$) of the underlying network, with the quantitative relation ${cd_{max}\sim d^{\alpha}}$. Our findings can help to design a better mitigation strategy for epidemics.