Spectroscopy of drums and quantum billiards: perturbative and non-perturbative results

We develop powerful numerical and analytical techniques for the solution of the Helmholtz equation on general domains. We prove two theorems: the first theorem provides an exact formula for the ground state of an arbirtrary membrane, while the second theorem generalizes this result to any excited state of the membrane. We also develop a systematic perturbative scheme which can be used to study the small deformations of a membrane of circular or square shapes. We discuss several applications, obtaining numerical and analytical results.Comment: 29 pages, 12 figures, 7 tabl

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

Scale-Free topologies and Activatory-Inhibitory interactions

A simple model of activatory-inhibitory interactions controlling the activity of agents (substrates) through a "saturated response" dynamical rule in a scale-free network is thoroughly studied. After discussing the most remarkable dynamical features of the model, namely fragmentation and multistability, we present a characterization of the temporal (periodic and chaotic) fluctuations of the quasi-stasis asymptotic states of network activity. The double (both structural and dynamical) source of entangled complexity of the system temporal fluctuations, as an important partial aspect of the Correlation Structure-Function problem, is further discussed to the light of the numerical results, with a view on potential applications of these general results.Comment: Revtex style, 12 pages and 12 figures. Enlarged manuscript with major revision and new results incorporated. To appear in Chaos (2006

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

Square Patterns and Quasi-patterns in Weakly Damped Faraday Waves

Pattern formation in parametric surface waves is studied in the limit of weak viscous dissipation. A set of quasi-potential equations (QPEs) is introduced that admits a closed representation in terms of surface variables alone. A multiscale expansion of the QPEs reveals the importance of triad resonant interactions, and the saturating effect of the driving force leading to a gradient amplitude equation. Minimization of the associated Lyapunov function yields standing wave patterns of square symmetry for capillary waves, and hexagonal patterns and a sequence of quasi-patterns for mixed capillary-gravity waves. Numerical integration of the QPEs reveals a quasi-pattern of eight-fold symmetry in the range of parameters predicted by the multiscale expansion.Comment: RevTeX, 11 pages, 8 figure

Longitudinal response functions of 3H and 3He

Trinucleon longitudinal response functions R_L(q,omega) are calculated for q values up to 500 MeV/c. These are the first calculations beyond the threshold region in which both three-nucleon (3N) and Coulomb forces are fully included. We employ two realistic NN potentials (configuration space BonnA, AV18) and two 3N potentials (UrbanaIX, Tucson-Melbourne). Complete final state interactions are taken into account via the Lorentz integral transform technique. We study relativistic corrections arising from first order corrections to the nuclear charge operator. In addition the reference frame dependence due to our non-relativistic framework is investigated. For q less equal 350 MeV/c we find a 3N force effect between 5 and 15 %, while the dependence on other theoretical ingredients is small. At q greater equal 400 MeV/c relativistic corrections to the charge operator and effects of frame dependence, especially for large omega, become more important. In comparison with experimental data there is generally a rather good agreement. Exceptions are the responses at excitation energies close to threshold, where there exists a large discrepancy with experiment at higher q. Concerning the effect of 3N forces there are a few cases, in particular for the R_L of 3He, where one finds a much improved agreement with experiment if 3N forces are included.Comment: 26 pages, 9 figure

Systematic derivation of a rotationally covariant extension of the 2-dimensional Newell-Whitehead-Segel equation

An extension of the Newell-Whitehead-Segel amplitude equation covariant under abritrary rotations is derived systematically by the renormalization group method.Comment: 8 pages, to appear in Phys. Rev. Letters, March 18, 199

Selection of the scaling solution in a cluster coalescence model

The scaling properties of the cluster size distribution of a system of diffusing clusters is studied in terms of a simple kinetic mean field model. It is shown that a one parameter family of mathematically valid scaling solutions exists. Despite this, the kinetics reaches a unique scaling solution independent of initial conditions. This selected scaling solution is marginally physical; i.e., it is the borderline solution between the unphysical and physical branches of the family of solutions.Comment: 4 pages, 5 figure

Evaluation of rate law approximations in bottom-up kinetic models of metabolism.

BackgroundThe mechanistic description of enzyme kinetics in a dynamic model of metabolism requires specifying the numerical values of a large number of kinetic parameters. The parameterization challenge is often addressed through the use of simplifying approximations to form reaction rate laws with reduced numbers of parameters. Whether such simplified models can reproduce dynamic characteristics of the full system is an important question.ResultsIn this work, we compared the local transient response properties of dynamic models constructed using rate laws with varying levels of approximation. These approximate rate laws were: 1) a Michaelis-Menten rate law with measured enzyme parameters, 2) a Michaelis-Menten rate law with approximated parameters, using the convenience kinetics convention, 3) a thermodynamic rate law resulting from a metabolite saturation assumption, and 4) a pure chemical reaction mass action rate law that removes the role of the enzyme from the reaction kinetics. We utilized in vivo data for the human red blood cell to compare the effect of rate law choices against the backdrop of physiological flux and concentration differences. We found that the Michaelis-Menten rate law with measured enzyme parameters yields an excellent approximation of the full system dynamics, while other assumptions cause greater discrepancies in system dynamic behavior. However, iteratively replacing mechanistic rate laws with approximations resulted in a model that retains a high correlation with the true model behavior. Investigating this consistency, we determined that the order of magnitude differences among fluxes and concentrations in the network were greatly influential on the network dynamics. We further identified reaction features such as thermodynamic reversibility, high substrate concentration, and lack of allosteric regulation, which make certain reactions more suitable for rate law approximations.ConclusionsOverall, our work generally supports the use of approximate rate laws when building large scale kinetic models, due to the key role that physiologically meaningful flux and concentration ranges play in determining network dynamics. However, we also showed that detailed mechanistic models show a clear benefit in prediction accuracy when data is available. The work here should help to provide guidance to future kinetic modeling efforts on the choice of rate law and parameterization approaches