An algorithm for series expansions based on hierarchical rate equations

We propose a computational method to obtain series expansions in powers of time for general dynamical systems described by a set of hierarchical rate equations. The method is generally applicable to problems in both equilibrium and nonequilibrium statistical mechanics such as random sequential adsorption, diffusion-reaction dynamics, and Ising dynamics. New result of random sequential adsorption of dimers on a square lattice is presented.Comment: LaTeX, 9 pages including 1 figur

The Phase Structure of Supersymmetric Sp(2N_c) Gauge Theories with an Adjoint

We study the phase structure of N = 1 supersymmetric Sp(2N_c) gauge theories with 2N_f fundamentals, an adjoint, and vanishing superpotential. Using a-maximization, we derive analytic expressions for the values of N_f below which the first several gauge-invariant operators in the chiral ring violate the unitarity bound and become free fields. In doing so we are able to explicitly check previous conjectures about the behavior of this theory made by Luty, Schmaltz, and Terning. We then compare this to an analysis of the first two 'deconfined' dual descriptions based on the gauge groups Sp(2N_f+2) x SO(2N_c+5) and Sp(2N_f+2) x SO(4N_f+4) x Sp(2N_c+2), finding precise agreement. In particular, we find no evidence for non-obvious accidental symmetries or the appearance of a mixed phase in which one of the dual gauge groups becomes free.Comment: 18 pages, 2 figures; v2: added references to match JHEP versio

Exons, introns and DNA thermodynamics

The genes of eukaryotes are characterized by protein coding fragments, the exons, interrupted by introns, i.e. stretches of DNA which do not carry any useful information for the protein synthesis. We have analyzed the melting behavior of randomly selected human cDNA sequences obtained from the genomic DNA by removing all introns. A clear correspondence is observed between exons and melting domains. This finding may provide new insights in the physical mechanisms underlying the evolution of genes.Comment: 4 pages, 8 figures - Final version as published. See also Phys. Rev. Focus 15, story 1

Comment on "Why is the DNA denaturation transition first order?"

In this comment we argue that while the conclusions in the original paper (Y. Kafri, D. Mukamel and L. Peliti, Phys. Rev. Lett. 85, 4988 (2000)) are correct for asymptotically long DNA chains, they do not apply to the chains used in typical experiments. In the added last paragraph, we point out that for real DNA the average distance between denatured loops is not of the order of the persistence length of a single-stranded chain but much larger. This corroborates our reasoning that the double helix between loops is quite rigid, and thereby our conclusion.Comment: 1 page, REVTeX. Last paragraph adde

Higgs Boson Decays to Neutralinos in Low-Scale Gauge Mediation

We study the decays of a standard model-like MSSM Higgs boson to pairs of neutralinos, each of which subsequently decays promptly to a photon and a gravitino. Such decays can arise in supersymmetric scenarios where supersymmetry breaking is mediated to us by gauge interactions with a relatively light gauge messenger sector (M_{mess} < 100 TeV). This process gives rise to a collider signal consisting of a pair of photons and missing energy. In the present work we investigate the bounds on this scenario within the minimal supersymmetric standard model from existing collider data. We also study the prospects for discovering the Higgs boson through this decay mode with upcoming data from the Tevatron and the LHC.Comment: 18 pages, 5 figures, added references and discussion of neutralino couplings, same as journal versio

Monte Carlo simulation of melting transition on DNA nanocompartment

DNA nanocompartment is a typical DNA-based machine whose function is dependent of molecular collective effect. Fundamental properties of the device have been addressed via electrochemical analysis, fluorescent microscopy, and atomic force microscopy. Interesting and novel phenomena emerged during the switching of the device. We have found that DNAs in this system exhibit a much steep melting transition compared to ones in bulk solution or conventional DNA array. To achieve an understanding to this discrepancy, we introduced DNA-DNA interaction potential to the conventional Ising-like Zimm-Bragg theory and Peyrard-Bishop model of DNA melting. To avoid unrealistic numerical calculation caused by modification of the Peyrard-Bishop nonlinear Hamiltonian with the DNA-DNA interaction, we established coarse-gained Monte Carlo recursion relations by elucidation of five components of energy change during melting transition. The result suggests that DNA-DNA interaction potential accounts for the observed steep transition.Comment: 12 pages, 5 figure

Phase diagram for unzipping DNA with long-range interactions

We present a critique and extension of the mean-field approach to the mechanical pulling transition in bound polymer systems. Our model is motivated by the theoretically and experimentally important examples of adsorbed polymers and double-stranded DNA, and we focus on the case in which quenched disorder in the sequence of monomers is unimportant for the statistical mechanics. We show how including excluded volume interactions in the model affects the phase diagram for the critical pulling force, and we predict a re-entrancy phase at low temperatures which has not been previously discussed. We also consider the case of non-equilibrium pulling, in which the external force probes the local, rather than the global structure of the dsDNA or adsorbed polymer. The dynamics of the pulling transition in such experiments could illuminate the polymer's loop structure, which depends on the nature of excluded volume interactions.Comment: 4 pages, 2 figures; this version clarifies Eq. 8, and corrects errors in Fig.

A length-dynamic Tonks gas theory of histone isotherms

We find exact solutions to a new one-dimensional (1D) interacting particle theory and apply the results to the adsorption and wrapping of polymers (such as DNA) around protein particles (such as histones). Each adsorbed protein is represented by a Tonks gas particle. The length of each particle is a degree of freedom that represents the degree of DNA wrapping around each histone. Thermodynamic quantities are computed as functions of wrapping energy, adsorbed histone density, and bulk histone concentration (or chemical potential); their experimental signatures are also discussed. Histone density is found to undergo a two-stage adsorption process as a function of chemical potential, while the mean coverage by high affinity proteins exhibits a maximum as a function of the chemical potential. However, {\it fluctuations} in the coverage are concurrently maximal. Histone-histone correlation functions are also computed and exhibit rich two length scale behavior.Comment: 5 pp, 3 fig

DNA bubble dynamics as a quantum Coulomb problem

We study the dynamics of denaturation bubbles in double-stranded DNA on the basis of the Poland-Scheraga model. We demonstrate that the associated Fokker-Planck equation is equivalent to a Coulomb problem. Below the melting temperature the bubble lifetime is associated with the continuum of scattering states of the repulsive Coulomb potential, at the melting temperature the Coulomb potential vanishes and the underlying first exit dynamics exhibits a long time power law tail, above the melting temperature, corresponding to an attractive Coulomb potential, the long time dynamics is controlled by the lowest bound state. Correlations and finite size effects are discussed.Comment: 4 pages, 3 figures, revte

Thermal denaturation of fluctuating finite DNA chains: the role of bending rigidity in bubble nucleation

Statistical DNA models available in the literature are often effective models where the base-pair state only (unbroken or broken) is considered. Because of a decrease by a factor of 30 of the effective bending rigidity of a sequence of broken bonds, or bubble, compared to the double stranded state, the inclusion of the molecular conformational degrees of freedom in a more general mesoscopic model is needed. In this paper we do so by presenting a 1D Ising model, which describes the internal base pair states, coupled to a discrete worm like chain model describing the chain configurations [J. Palmeri, M. Manghi, and N. Destainville, Phys. Rev. Lett. 99, 088103 (2007)]. This coupled model is exactly solved using a transfer matrix technique that presents an analogy with the path integral treatment of a quantum two-state diatomic molecule. When the chain fluctuations are integrated out, the denaturation transition temperature and width emerge naturally as an explicit function of the model parameters of a well defined Hamiltonian, revealing that the transition is driven by the difference in bending (entropy dominated) free energy between bubble and double-stranded segments. The calculated melting curve (fraction of open base pairs) is in good agreement with the experimental melting profile of polydA-polydT. The predicted variation of the mean-square-radius as a function of temperature leads to a coherent novel explanation for the experimentally observed thermal viscosity transition. Finally, the influence of the DNA strand length is studied in detail, underlining the importance of finite size effects, even for DNA made of several thousand base pairs.Comment: Latex, 28 pages pdf, 9 figure