Coarse graining of master equations with fast and slow states

We propose a general method for simplifying master equations by eliminating from the description rapidly evolving states. The physical recipe we impose is the suppression of these states and a renormalization of the rates of all the surviving states. In some cases, this decimation procedure can be analytically carried out and is consistent with other analytical approaches, like in the problem of the random walk in a double-well potential. We discuss the application of our method to nontrivial examples: diffusion in a lattice with defects and a model of an enzymatic reaction outside the steady state regime.Comment: 9 pages, 9 figures, final version (new subsection and many minor improvements

Rigorous elimination of fast stochastic variables from the linear noise approximation using projection operators

The linear noise approximation (LNA) offers a simple means by which one can study intrinsic noise in monostable biochemical networks. Using simple physical arguments, we have recently introduced the slow-scale LNA (ssLNA) which is a reduced version of the LNA under conditions of timescale separation. In this paper, we present the first rigorous derivation of the ssLNA using the projection operator technique and show that the ssLNA follows uniquely from the standard LNA under the same conditions of timescale separation as those required for the deterministic quasi-steady state approximation. We also show that the large molecule number limit of several common stochastic model reduction techniques under timescale separation conditions constitutes a special case of the ssLNA.Comment: 10 pages, 1 figure, submitted to Physical Review E; see also BMC Systems Biology 6, 39 (2012

Generalized Haldane Equation and Fluctuation Theorem in the Steady State Cycle Kinetics of Single Enzymes

Enyzme kinetics are cyclic. We study a Markov renewal process model of single-enzyme turnover in nonequilibrium steady-state (NESS) with sustained concentrations for substrates and products. We show that the forward and backward cycle times have idential non-exponential distributions: \QQ_+(t)=\QQ_-(t). This equation generalizes the Haldane relation in reversible enzyme kinetics. In terms of the probabilities for the forward ($p_+$) and backward ($p_-$) cycles, $k_BT\ln(p_+/p_-)$ is shown to be the chemical driving force of the NESS, $\Delta\mu$. More interestingly, the moment generating function of the stochastic number of substrate cycle $\nu(t)$, follows the fluctuation theorem in the form of Kurchan-Lebowitz-Spohn-type symmetry. When $\lambda$ = $\Delta\mu/k_BT$, we obtain the Jarzynski-Hatano-Sasa-type equality: $\equiv$ 1 for all $t$, where $\nu\Delta\mu$ is the fluctuating chemical work done for sustaining the NESS. This theory suggests possible methods to experimentally determine the nonequilibrium driving force {\it in situ} from turnover data via single-molecule enzymology.Comment: 4 pages, 3 figure

Shear induced grain boundary motion for lamellar phases in the weakly nonlinear regime

We study the effect of an externally imposed oscillatory shear on the motion of a grain boundary that separates differently oriented domains of the lamellar phase of a diblock copolymer. A direct numerical solution of the Swift-Hohenberg equation in shear flow is used for the case of a transverse/parallel grain boundary in the limits of weak nonlinearity and low shear frequency. We focus on the region of parameters in which both transverse and parallel lamellae are linearly stable. Shearing leads to excess free energy in the transverse region relative to the parallel region, which is in turn dissipated by net motion of the boundary toward the transverse region. The observed boundary motion is a combination of rigid advection by the flow and order parameter diffusion. The latter includes break up and reconnection of lamellae, as well as a weak Eckhaus instability in the boundary region for sufficiently large strain amplitude that leads to slow wavenumber readjustment. The net average velocity is seen to increase with frequency and strain amplitude, and can be obtained by a multiple scale expansion of the governing equations

A Study of the Quasi-elastic (e,e'p) Reaction on 12^{12}C, 56^{56}Fe and 97^{97}Au

We report the results from a systematic study of the quasi-elastic (e,e'p) reaction on $^{12}$C, $^{56}$Fe and $^{197}$Au performed at Jefferson Lab. We have measured nuclear transparency and extracted spectral functions (corrected for radiation) over a Q$^2$ range of 0.64 - 3.25 (GeV/c)$^2$ for all three nuclei. In addition we have extracted separated longitudinal and transverse spectral functions at Q$^2$ of 0.64 and 1.8 (GeV/c)$^2$ for these three nuclei (except for $^{197}$Au at the higher Q$^2$). The spectral functions are compared to a number of theoretical calculations. The measured spectral functions differ in detail but not in overall shape from most of the theoretical models. In all three targets the measured spectral functions show considerable excess transverse strength at Q$^2$ = 0.64 (GeV/c)$^2$, which is much reduced at 1.8 (GeV/c)$^2$.Comment: For JLab E91013 Collaboration, 19 pages, 20 figures, 3 table

Determinants of impact : towards a better understanding of encounters with the arts

The article argues that current methods for assessing the impact of the arts are largely based on a fragmented and incomplete understanding of the cognitive, psychological and socio-cultural dynamics that govern the aesthetic experience. It postulates that a better grasp of the interaction between the individual and the work of art is the necessary foundation for a genuine understanding of how the arts can affect people. Through a critique of philosophical and empirical attempts to capture the main features of the aesthetic encounter, the article draws attention to the gaps in our current understanding of the responses to art. It proposes a classification and exploration of the factors—social, cultural and psychological—that contribute to shaping the aesthetic experience, thus determining the possibility of impact. The ‘determinants of impact’ identified are distinguished into three groups: those that are inherent to the individual who interacts with the artwork; those that are inherent to the artwork; and ‘environmental factors’, which are extrinsic to both the individual and the artwork. The article concludes that any meaningful attempt to assess the impact of the arts would need to take these ‘determinants of impact’ into account, in order to capture the multidimensional and subjective nature of the aesthetic experience

Spin Structure of the Proton from Polarized Inclusive Deep-Inelastic Muon-Proton Scattering

We have measured the spin-dependent structure function $g_1^p$ in inclusive deep-inelastic scattering of polarized muons off polarized protons, in the kinematic range $0.003 < x < 0.7$ and $1 GeV^2 < Q^2 < 60 GeV^2$. A next-to-leading order QCD analysis is used to evolve the measured $g_1^p(x,Q^2)$ to a fixed $Q^2_0$. The first moment of $g_1^p$ at $Q^2_0 = 10 GeV^2$ is $\Gamma^p = 0.136\pm 0.013(stat.) \pm 0.009(syst.)\pm 0.005(evol.)$. This result is below the prediction of the Ellis-Jaffe sum rule by more than two standard deviations. The singlet axial charge $a_0$ is found to be $0.28 \pm 0.16$. In the Adler-Bardeen factorization scheme, $\Delta g \simeq 2$ is required to bring $\Delta \Sigma$ in agreement with the Quark-Parton Model. A combined analysis of all available proton and deuteron data confirms the Bjorken sum rule.Comment: 33 pages, 22 figures, uses ReVTex and smc.sty. submitted to Physical Review

Moment Closure - A Brief Review

Moment closure methods appear in myriad scientific disciplines in the modelling of complex systems. The goal is to achieve a closed form of a large, usually even infinite, set of coupled differential (or difference) equations. Each equation describes the evolution of one "moment", a suitable coarse-grained quantity computable from the full state space. If the system is too large for analytical and/or numerical methods, then one aims to reduce it by finding a moment closure relation expressing "higher-order moments" in terms of "lower-order moments". In this brief review, we focus on highlighting how moment closure methods occur in different contexts. We also conjecture via a geometric explanation why it has been difficult to rigorously justify many moment closure approximations although they work very well in practice.Comment: short survey paper (max 20 pages) for a broad audience in mathematics, physics, chemistry and quantitative biolog