Chain motion and viscoelasticity in highly entangled solutions of semiflexible rods

Brownian dynamics simulations are used to study highly entangled solutions of semiflexible polymers. Bending fluctuations of semiflexible rods are signficantly affected by entanglement only above a concentration $c^{**}$, where $c^{**}\sim 10^{3}L^{-3}$ for chains of similar length $L$ and persistence length. For $c > c^{**}$, the tube radius $R_{e}$ approaches a dependence $R_{e} \propto c^{-3/5}$, and the linear viscoelastic response develops an elastic contribution that is absent for $c < c^{**}$. Experiments on isotropic solutions of $F$-actin span concentrations near $c^{**}$ for which the predicted asymptotic scaling of the plateau modulus $G \propto c^{7/5}$ is not yet valid.Comment: 4 pages, 5 figures, submitted to PR

Response of Single Polymers to Localized Step Strains

In this paper, the response of single three-dimensional phantom and self-avoiding polymers to localized step strains are studied for two cases in the absence of hydrodynamic interactions: (i) polymers tethered at one end with the strain created at the point of tether, and (ii) free polymers with the strain created in the middle of the polymer. The polymers are assumed to be in their equilibrium state before the step strain is created. It is shown that the strain relaxes as a power-law in time $t$ as $t^{-\eta}$. While the strain relaxes as $1/t$ for the phantom polymer in both cases; the self-avoiding polymer relaxes its strain differently in case (i) than in case (ii): as $t^{-(1+\nu)/(1+2\nu)}$ and as $t^{-2/(1+2\nu)}$ respectively. Here $\nu$ is the Flory exponent for the polymer, with value $\approx0.588$ in three dimensions. Using the mode expansion method, exact derivations are provided for the $1/t$ strain relaxation behavior for the phantom polymer. However, since the mode expansion method for self-avoiding polymers is nonlinear, similar theoretical derivations for the self-avoiding polymer proves difficult to provide. Only simulation data are therefore presented in support of the $t^{-(1+\nu)/(1+2\nu)}$ and the $t^{-2/(1+2\nu)}$ behavior. The relevance of these exponents for the anomalous dynamics of polymers are also discussed.Comment: 10 pages, 1 figure; minor errors corrected, introduction slightly modified and references expanded; to appear in Phys. Rev.

On a modification method of Lefschetz thimbles

The QCD at finite density is not well understood yet, where standard Monte Carlo simulation suffers from the sign problem. In order to overcome the sign problem, the method of Lefschetz thimble has been explored. Basically, the original sign problem can be less severe in a complexified theory due to the constancy of the imaginary part of an action on each thimble. However, global phase factors assigned on each thimble still remain. Their interference is not negligible in a situation where a large number of thimbles contribute to the partition function, and this could also lead to a sign problem.In this study, we propose a method to resolve this problem by modifying the structure of Lefschetz thimbles such that only a single thimble is relevant to the partition function. It can be shown that observables measured in the original and modified theories are connected by a simple identity. We exemplify that our method works well in a toy model.Comment: 7 pages, 4 figures, talk presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spai

A non-monotonic constitutive model is not necessary to obtain shear banding phenomena in entangled polymer solutions

In 1975 Doi and Edwards predicted that entangled polymer melts and solutions can have a constitutive instability, signified by a decreasing stress for shear rates greater than the inverse of the reptation time. Experiments did not support this, and more sophisticated theories incorporated Marrucci's idea (1996) of removing constraints by advection; this produced a monotonically increasing stress and thus stable constitutive behavior. Recent experiments have suggested that entangled polymer solutions may possess a constitutive instability after all, and have led some workers to question the validity of existing constitutive models. In this Letter we use a simple modern constitutive model for entangled polymers, the non-stretching Rolie-Poly model with an added solvent viscosity, and show that (1) instability and shear banding is captured within this simple class of models; (2) shear banding phenomena is observable for weakly stable fluids in flow geometries that impose a sufficiently inhomogeneous total shear stress; (3) transient phenomena can possess inhomogeneities that resemble shear banding, even for weakly stable fluids. Many of these results are model-independent.Comment: 5 figure

Relation between Confinement and Chiral Symmetry Breaking in Temporally Odd-number Lattice QCD

In the lattice QCD formalism, we investigate the relation between confinement and chiral symmetry breaking. A gauge-invariant analytical relation connecting the Polyakov loop and the Dirac modes is derived on a temporally odd-number lattice, where the temporal lattice size is odd, with the normal (nontwisted) periodic boundary condition for link-variables. This analytical relation indicates that low-lying Dirac modes have little contribution to the Polyakov loop, and it is numerically confirmed at the quenched level in both confinement and deconfinement phases. This fact indicates no direct one-to-one correspondence between confinement and chiral symmetry breaking in QCD. Using the relation, we also investigate the contribution from each Dirac mode to the Polyakov loop. In the confinement phase, we find a new "positive/negative symmetry" of the Dirac-mode matrix element of the link-variable operator, and this symmetry leads to the zero value of the Polyakov loop. In the deconfinement phase, there is no such symmetry and the Polyakov loop is nonzero. Also, we develop a new method for spin-diagonalizing the Dirac operator on the temporally odd-number lattice modifying the Kogut-Susskind formalism.Comment: 15pages, 9 figure

Kinetic Regimes and Cross-Over Times in Many-Particle Reacting Systems

We study kinetics of single species reactions ("A+A -> 0") for general local reactivity Q and dynamical exponent z (rms displacement x_t ~ t^{1/z}.) For small molecules z=2, whilst z=4,8 for certain polymer systems. For dimensions d above the critical value d_c=z, kinetics are always mean field (MF). Below d_c, the density n_t initially follows MF decay, n_0 - n_t ~ n_0^2 Q t. A 2-body diffusion-controlled regime follows for strongly reactive systems (Q>Qstar ~ n_0^{(z-d)/d}) with n_0 - n_t ~ n_0^2 x_t^d. For Q<Qstar, MF kinetics persist, with n_t ~ 1/Qt. In all cases n_t ~ 1/x_t^d at the longest times. Our analysis avoids decoupling approximations by instead postulating weak physically motivated bounds on correlation functions.Comment: 10 pages, 1 figure, uses bulk2.sty, minor changes, submitted to Europhysics Letter

Low-Dimensional Fluctuations and Pseudogap in Gaudin-Yang Fermi Gases

Pseudogap is a ubiquitous phenomenon in strongly correlated systems such as high-$T_{\rm c}$ superconductors, ultracold atoms and nuclear physics. While pairing fluctuations inducing the pseudogap are known to be enhanced in low-dimensional systems, such effects have not been explored well in one of the most fundamental 1D models, that is, Gaudin-Yang model. In this work, we show that the pseudogap effect can be visible in the single-particle excitation in this system using a diagrammatic approach. Fermionic single-particle spectra exhibit a unique crossover from the double-particle dispersion to pseudogap state with increasing the attractive interaction and the number density at finite temperature. Surprisingly, our results of thermodynamic quantities in unpolarized and polarized gases show an excellent agreement with the recent quantum Monte Carlo and complex Langevin results, even in the region where the pseudogap appears.Comment: 6 pages, 5 figure

Cluster mean-field approximations with the coherent-anomaly-method analysis for the driven pair contact process with diffusion

The cluster mean-field approximations are performed, up to 13 cluster sizes, to study the critical behavior of the driven pair contact process with diffusion (DPCPD) and its precedent, the PCPD in one dimension. Critical points are estimated by extrapolating our data to the infinite cluster size limit, which are in good accordance with recent simulation results. Within the cluster mean-field approximation scheme, the PCPD and the DPCPD share the same mean-field critical behavior. The application of the coherent anomaly method, however, shows that the two models develop different coherent anomalies, which lead to different true critical scaling. The values of the critical exponents for the particle density, the pair density, the correlation length, and the relaxation time are fairly well estimated for the DPCPD. These results support and complement our recent simulation results for the DPCPD

On the origin of the unusual behavior in the stretching of single-stranded DNA

Force extension curves (FECs), which quantify the response of a variety of biomolecules subject to mechanical force ($f$), are often quantitatively fit using worm-like chain (WLC) or freely-jointed chain (FJC) models. These models predict that the chain extension, $x$, normalized by the contour length increases linearly at small $f$ and at high forces scale as $x \sim (1 - f^{-\alpha})$ where $\alpha$= 0.5 for WLC and unity for FJC. In contrast, experiments on ssDNA show that over a range of $f$ and ionic concentration, $x$ scales as $x\sim\ln f$, which cannot be explained using WLC or FJC models. Using theory and simulations we show that this unusual behavior in FEC in ssDNA is due to sequence-independent polyelectrolyte effects. We show that the $x\sim \ln f$ arises because in the absence of force the tangent correlation function, quantifying chain persistence, decays algebraically on length scales on the order of the Debye length. Our theory, which is most appropriate for monovalent salts, quantitatively fits the experimental data and further predicts that such a regime is not discernible in double stranded DNA.Comment: Accepted for publication in JC