Physical review. E最新文献

Robert May famously used random matrix theory to predict that large, complex systems cannot admit stable fixed points. However, this general conclusion is not always supported by empirical observation: from cells to biomes, biological systems are large, complex, and often stable. In this paper, we revisit May's argument in light of recent developments in both ecology and random matrix theory. We focus on competitive systems, and, using a nonlinear generalization of the competitive Lotka-Volterra model, we show that there are, in fact, two kinds of complexity-stability relationships in disordered dynamical systems: if self-interactions grow faster with density than cross-interactions, complexity is destabilizing; but if cross-interactions grow faster than self-interactions, complexity is stabilizing.

We elaborate and validate a generalization of the renowned transition-path-sampling algorithm for a paradigmatic model of active particles, namely, the run-and-tumble particles. Notwithstanding the nonequilibrium character of these particles, we show how the consequent lack of the microscopical reversibility property, which is usually required by transition-path sampling, can be circumvented by identifying reasonable backward dynamics with a well-defined path-probability density. Our method is then applied to characterize the structure and kinetics of rare transition pathways undergone by run-and-tumble particles having to cross a potential barrier in order to find a target.

We develop convergence acceleration procedures that enable a gradient descent-type iteration method to efficiently simulate Hartree-Fock equations for many particles interacting both with each other and with an external potential. Our development focuses on three aspects: (i) optimization of a parameter in the preconditioning operator; (ii) adoption of a technique that eliminates the slowest-decaying mode to the case of many equations (describing many particles); and (iii) a novel extension of the above technique that allows one to eliminate multiple modes simultaneously. We illustrate performance of the numerical method for the two-dimensional model of the first layer of helium atoms above a graphene sheet. We demonstrate that incorporation of aspects (i) and (ii) above into the "plain" gradient descent method accelerates it by at least two orders of magnitude, and often by much more. Aspect (iii)-the multiple-mode elimination-may bring further improvement to the convergence rate compared to aspect (ii), the single-mode elimination. Both single- and multiple-mode elimination techniques are shown to significantly outperform the well-known Anderson Acceleration. We believe that our acceleration techniques can also be employed by other iterative methods, especially those handling hard-core-type interaction of many particles.

The ultimate goal of this paper is to develop a systematic method for deriving quantum master equations that satisfy the requirements of a completely positive and trace-preserving (CPTP) map, further describing thermal relaxation processes. In this paper, we assume that the quantum master equation is obtained through the canonical quantization of the generalized Brownian motion proposed in our recent paper [T. Koide and F. Nicacio, Phys. Lett. A 494, 129277 (2024)0375-960110.1016/j.physleta.2023.129277]. At least classically, this dynamics describes the thermal relaxation process regardless of the choice of the system Hamiltonian. The remaining task is to identify the parameters ensuring that the quantum master equation meets complete positivity. We limit our discussion to many-body quadratic Hamiltonians and establish a CPTP criterion for our quantum master equation. This criterion is useful for applying our quantum master equation to models with interaction such as a network model, which has been used to investigate how quantum effects modify heat conduction.

We propose a reaction-diffusion system that converts topological information of an active nematic into chemical signals. We show that a curvature-activated reaction dipole is sufficient for creating a system that dynamically senses topology by producing a concentration field possessing local extrema coinciding with ±1/2 defects. The enabling term is analogous to polarization charge density seen in dielectric materials. We demonstrate the ability of this system to identify defects in both passive and active nematics. Our results illustrate that a relatively simple feedback scheme, expressed as a system of partial differential equations, is capable of producing chemical signals in response to inherently nonlocal structures in anisotropic media. We posit that such coarse-grained systems can help generate testable hypotheses for regulated processes in biological systems, such as morphogenesis, and motivate the creation of bio-inspired materials that utilize dynamic coupling between nematic structure and biochemistry.

The bacterial flagellar motor is both chemo- and mechanosensitive. It is sensitive to the intracellular concentration of the chemotaxis response regulator CheY-P and to external load conditions. The motor's dose-response curve, which represents the probability of the motor rotating clockwise (CW bias) as a function of CheY-P concentration, characterizes its chemical sensitivity. However, it remains unclear how this dose-response curve depends on the load conditions. Here, we measured the dose-response curves of the motor under various load conditions. Surprisingly, we found that the dose-response curve exhibited minimal changes with load at low CW biases, but shifted leftward with higher sensitivity to CheY-P concentration at high CW biases when the load increased. This observation contradicts previous model predictions that incorporated the effect of stator-rotor interaction on motor switching. Through the development of an Ising-type model for the coupled chemo- and mechanosensitivity of the flagellar switch, we revealed that the mechanism underlying the mechanosensitive dose response is the synergistic interplay between the adaptive remodeling of the motor switch complex and the nonequilibrium effect of the stator-rotor interaction.

We investigate the dynamical phases and phase transitions arising in a classical two-dimensional anisotropic XY model under the influence of a periodically driven temporal external magnetic field in the form of a symmetric square wave. We use a combination of finite temperature classical Monte Carlo simulation, implemented within a CPU+GPU paradigm, utilizing local dynamics provided by the Glauber algorithm and a phenomenological equation-of-motion approach based on relaxational dynamics governed by the time-dependent free energy within a mean-field approximation to study the model. We investigate several parameter regimes of the variables (magnetic field, anisotropy, and the external drive frequency) that influence the anisotropic XY system. We identify four possible dynamical phases: Ising-SBO, Ising-SRO, XY-SBO, and XY-SRO. Both techniques indicate that only three of them (Ising-SRO, Ising-SBO, and XY-SRO) are stable dynamical phases in the thermodynamic sense. Within the Monte Carlo framework, a finite-size scaling analysis, shows that XY-SBO does not survive in the thermodynamic limit giving way to either an Ising-SBO or a XY-SRO regime. The finite-size scaling analysis further shows that the transitions between the three remaining dynamical phases either belong to the two-dimensional Ising universality class or are first-order in nature. Within the mean-field calculations yield three stable dynamical phases, i.e., Ising-SRO, Ising-SBO and XY-SRO, where the final steady state is independent of the initial condition chosen to evolve the equations of motion, as well as a region of bistability where the system flows to either Ising-SBO or XY-SRO (Ising-SRO) depending on the initial condition. Unlike the stable dynamical phases, the XY-SBO represents a transient feature that is eventually lost to either Ising-SBO or XY-SRO. Our mean-field analysis highlights the importance of the competition between switching of the stationary point(s) of the free energy after each half cycle of the external field and the two-dimensional nature of the phase space for the equations of motion.

We explore how the interplay of finite availability, carrying capacity of particles at different parts of a spatially extended system, and particle diffusion between them control the steady-state currents and density profiles in a one-dimensional current-carrying channel connecting the different parts of the system. To study this, we construct a minimal model consisting of two particle reservoirs of finite carrying capacities connected by a totally asymmetric simple exclusion process (TASEP). In addition to particle transport via TASEP between the reservoirs, the latter can also directly exchange particles via Langmuir kinetics-like processes, modeling particle diffusion between them that can maintain a steady current in the system. We calculate the steady-state density profiles and the associated particle currents in the TASEP lane. The resulting phases and the phase diagrams are quite different from an open TASEP, and are characterized by the model parameters defining particle exchanges between the TASEP and the reservoirs, direct particle exchanges between the reservoirs, and the filling fraction of the particles that determines the total resources available. These parameters can be tuned to make the density on the TASEP lane globally uniform or piecewise continuous, and can make the two reservoirs preferentially populated or depopulated.

Buckling instabilities driven by tissue growth underpin key developmental events such as the folding of the brain. Tissue growth is disordered due to cell-to-cell variability, but the effects of this variability on buckling are unknown. Here, we analyze what is perhaps the simplest setup of this problem: the buckling of an elastic rod with fixed ends driven by spatially varying, yet highly symmetric growth. Combining analytical calculations for simple growth fields and numerical sampling of random growth fields, we show that variability can increase as well as decrease the growth threshold for buckling, even when growth variability does not cause any residual stresses. For random growth, we find numerically that the shift of the buckling threshold correlates with spatial moments of the growth field. Our results imply that biological systems can either trigger or avoid buckling by exploiting the spatial arrangement of growth variability.

By means of analytical calculations and numerical simulations, we study the diffusion properties in quasi-two-dimensional structures with two exciton subsystems with an exchange between them. The experimental realization is possible in systems where fast and slow exciton subsystems appear. For substantially different diffusion coefficients of the species, the negative diffusion can be observed if one measures the transport properties of only a single subsystem, just as was obtained in experimental studies for quasi-two-dimensional semiconductor systems. The initial transition from a fast subsystem to a slow one results in a delayed release of fast excitons in the area close to the original excitation spot. Hence, the signal from the fast excitons alone includes the delayed replenishment from the transition of slow quasiparticles. This is seen as the narrowing of the exciton density profile or decrease of mean-squared displacement which is then interpreted as a negative diffusion. We show that the analytical theory matches the available experimental data for negative diffusion quite well. The average diffusion coefficients for the combined population are analytically expressed through the diffusion coefficients of fast and slow excitons. Simple analytical expressions for effective diffusion coefficients obtained in limiting cases are of interest both for theoretical and experimental analysis of not only the exciton transport, but also for a variety of systems, where fast and slow moving subsystems are present.