We investigate spike-timing dependent plasticity (STPD) in the case of a synapse connecting two neuronal cells. We develop a theoretical analysis of several STDP rules using Markovian theory. In this context there are two different timescales, fast neuronal activity and slower synaptic weight updates. Exploiting this timescale separation, we derive the long-time limits of a single synaptic weight subject to STDP. We show that the pairing model of presynaptic and postsynaptic spikes controls the synaptic weight dynamics for small external input on an excitatory synapse. This result implies in particular that mean-field analysis of plasticity may miss some important properties of STDP. Anti-Hebbian STDP favors the emergence of a stable synaptic weight. In the case of an inhibitory synapse the pairing schemes matter less, and we observe convergence of the synaptic weight to a nonnull value only for Hebbian STDP. We extensively study different asymptotic regimes for STDP rules, raising interesting questions for future work on adaptative neuronal networks and, more generally, on adaptative systems.

We have used kinetic Monte Carlo (kMC) simulations of a lattice gas to study front fluctuations in the spreading of a nonvolatile liquid droplet onto a solid substrate. Our results are consistent with a diffusive growth law for the radius of the precursor layer, $R\sim t^{\delta }$, with $\delta \approx 1/2$ in all the conditions considered for temperature and substrate wettability, in good agreement with previous studies. The fluctuations of the front exhibit kinetic roughening properties with exponent values which depend on temperature $T$, but become $T$ independent for sufficiently high $T$. Moreover, strong evidence of intrinsic anomalous scaling has been found, characterized by different values of the roughness exponent at short and large length scales. Although such a behavior differs from the scaling properties of the one-dimensional Kardar-Parisi-Zhang (KPZ) universality class, the front covariance and the probability distribution function of front fluctuations found in our kMC simulations do display KPZ behavior, agreeing with simulations of a continuum height equation proposed in this context. However, this equation does not feature intrinsic anomalous scaling, at variance with the discrete model.

We investigate the effects of physicochemical attractions on the transport of finite-sized particles in three-dimensional ordered nanopost arrays using Stokesian dynamics simulations. We find that weak particle-nanopost attractions negligibly affect diffusion due to the dominance of Brownian fluctuations. Strong attractions, however, significantly hinder particle diffusion due to localization of particles around the nanoposts. Conversely, under flow, attractions significantly enhance longitudinal dispersion at low to moderate Péclet number (Pe). At high Pe, by contrast, advection becomes dominant and attractions weakly enhance dispersion. Moreover, attractions frustrate directional locking at moderate flow rates, and shift the onset of this behavior to higher Pe.

The impact of magnetic field (MF) on the parallel resistivity ${\eta }_{\parallel }$ is studied for strongly magnetized plasmas with the electron thermal gyroradius ${\rho }_{\text{the}}$ smaller than the Debye length ${\lambda }_D$ but much larger than the Landau length ${\lambda }_L$. Two previous papers [P. Ghendrih et al., Phys. Lett. A 119, 354 (1987); S. D. Baalrud and T. Lafleur, Phys. Plasmas 28, 102107 (2021)] found ${\eta }_{\parallel }$ to increase monotonically with MF. Unfortunately, both works used predetermined electron distribution functions and are thus not self-consistent. In this paper, we analyze the MF dependence of ${\eta }_{\parallel }$ self-consistently by solving the electron magnetized kinetic equation in a Lorentz gaslike approximation. It is found ${\eta }_{\parallel }$ decreases monotonically with MF, with ${\lambda }_D$ in the usual Coulomb logarithm $ln\mathrm{\Lambda }=ln({\lambda }_D/{\lambda }_L)$ being replaced by ${\rho }_{\text{the}}$. The underlying physics is that the electrons affected only by the collisions with impact parameters between ${\lambda }_L$ and ${\rho }_{\text{the}}$ carry almost all the parallel current.

Nuclear magnetic resonance (NMR) responses of fluids saturating porous media arise from complex relaxation-diffusion dynamics of polarized spins. These constitute a sensitive probe of the microstructure and are described by the Bloch-Torrey equations. An NMR simulation framework based on an augmented lattice Boltzmann method aimed at the fine-scale resolution of nuclear polarization density is presented. The approach encapsulates the time evolution of the full magnetization vector and naturally incorporates the mechanisms of diffusional transport. Spin dephasing mechanisms are fully resolved at tomogram voxel scale to account for magnetic field inhomogeneity. The approach is validated against analytical solutions of spin-echo decays for simple pore geometries. An application to a nano-computed-tomography image of chalk with inhomogeneous internal fields yields $T_2$ spectral measures in good agreement with experiment and illustrates the spatial pore-scale dynamics of net magnetization. Findings establish the feasibility of the framework for pure diffusion and present an approach vector to modeling the evolution of magnetization under flow conditions.

We analyze continuous-time Markovian $\varepsilon \text{−}\mathrm{SIS}$ epidemics with self-infections on the complete graph. The majority of the graphs are analytically intractable, but many physical features of the $\varepsilon \text{−}\mathrm{SIS}$ process observed in the complete graph can occur in any other graph. In this work, we illustrate that the timescales of the $\varepsilon \text{−}\mathrm{SIS}$ process are related to the eigenvalues of the tridiagonal matrix of the SIS Markov chain. We provide a detailed analysis of all eigenvalues and illustrate that the eigenvalues show staircases, which are caused by the nearly degenerate (but strictly distinct) pairs of eigenvalues. We also illustrate that the ratio between the second-largest and third-largest eigenvalue is a good indicator of metastability in the $\varepsilon \text{−}\mathrm{SIS}$ process. Additionally, we show that the epidemic threshold of the Markovian $\varepsilon \text{−}\mathrm{SIS}$ process can be accurately approximated by the effective infection rate for which the third-largest eigenvalue of the transition matrix is the smallest. Finally, we derive the exact mean-field solution for the $\varepsilon \text{−}\mathrm{SIS}$ process on the complete graph, and we show that the mean-field approximation does not correctly represent the metastable behavior of Markovian $\varepsilon \text{−}\mathrm{SIS}$ epidemics.

Stochastic excitability of spiking oscillatory regimes in the calcium kinetics is studied on the basis of the Li-Rinzel conceptual model. The probabilistic mechanisms of the noise-induced generation of large-amplitude oscillations in parametric zones, where the original deterministic model has only stable equilibria, are investigated numerically and analytically. A parametric statistical description of interspike intervals is curried out and the phenomenon of coherence resonance is discussed. For the analytical study of the stochastic excitement, the confidence domain method using a stochastic sensitivity technique is applied. In this analysis, a key role of mutual arrangement of the confidence ellipses and separatrices detaching the sub- and supercritical regions is demonstrated. It is shown that in the Li-Rinzel model such separatrices are the stable manifolds of the saddle equilibria and the transient semiattractors.

Chirality is a design feature of a number of biomolecules (e.g., collagen). In these molecules, cholesteric (chiral-nematic) behavior emerges from a combination of the tendency for the biopolymers to align (nematic interactions) and for the alignment direction to change with position, rotating around an axis normal to the alignment direction. This paper presents self-consistent field theory (SCFT) of chiral-nematic polymers, which takes into account polymer flexibility and the orientational degrees of freedom of polymer segments. Using the resulting SCFT, we construct a phase diagram showing regions of stability for isotropic, nematic, and cholesteric phases. Furthermore, we find that nematic interactions can stabilize the cholesteric phase, pushing the isotropic-cholesteric phase transition to lower cholesteric interaction strength, until the isotropic-nematic-cholesteric triple point is reached.

Polymerized ionic liquids (PILs) are typically single-ion conductors, where one kind of ionic species is either placed as the pendant group to the chain (pendant PILs) or directly incorporated into the polymeric backbone (backbone PILs). This paper compares the thermodynamics, ionic dynamics, and mechanical properties of pendant and backbone PILs. The results indicate that near the glass transition, the energy barrier for ion hopping is much lower for pendant PIL while the backbone PIL shows a much stronger sensitivity to pressure. At the same time, a free-volume based model was proposed here to understand the ion dynamics of both studied PILs at high-pressure conditions. The determined critical volume, quantifying the minimal volume required for ion hopping, of the pendant PIL is significantly reduced compared to the backbone PIL, which is most likely the reason for the enhanced ionic conductivity of the pendant PIL near the glass transition. We found that the proposed model is equivalent to the commonly used pressure counterpart of the Vogel-Fulcher-Tammann equation.

The interaction of relativistic short-pulse lasers with matter produces fast electrons with over megaampere currents, which supposedly heats a solid target isochorically and forms a hot dense plasma. In a picosecond timescale, however, thermal diffusion from hot preformed plasma turns out to be the dominant process of isochoric heating. We describe a heating process, fast thermal diffusion, launched from the preformed plasma heated resistively by the fast electron current. We demonstrate the fast thermal diffusion in the keV range in a solid density plasma by a series of one-dimensional particle-in-cell simulations. A theoretical model of the fast thermal diffusion is developed and we derive the diffusion speed as a function of the laser amplitude and target density. Under continuous laser irradiation, the diffusion front propagates at a constant speed in uniform plasma. Our model can provide a guideline for fast isochoric heating using future kilojoule petawatt lasers.

We develop a three-dimensional (3D) fully self-consistent model for analysis of an ultrashort THz pulse propagation and amplification in a nonequilibrium plasma channel formed in xenon by a femtosecond UV laser pulse. The model is based on the self-consistent solution of a second order wave equation in the cylindrical geometry and the kinetic Boltzmann equation for the electron velocity distribution function (EVDF) at different points of the spatially inhomogeneous nonequilibrium plasma channel. We analyze the wide range of plasma and seed pulse parameters and reveal the optimal regimes for producing high intensity outgoing THz fields as well as highly unipolar THz pulses within the proposed mechanism. It is demonstrated that the process of EVDF relaxation in plasma limits the amplification of THz pulses at the level of $\sim {10}^7\mathrm{W}/\mathrm{c}{\mathrm{m}}^2$. Both focusing features of nonequilibrium plasma and the possibility of producing THz pulses with a high degree of unipolarity are confirmed for the case of 3D geometry.

We study classical Ising spin-$\frac{1}{2}$ models on a two-dimensional (2D) square lattice with ferromagnetic or antiferromagnetic nearest-neighbor interactions, under the effect of a pure imaginary magnetic field. The complex Boltzmann weights of spin configurations cannot be interpreted as a probability distribution, which prevents application of standard statistical algorithms. In this work, the mapping of the Ising spin models under consideration onto symmetric vertex models leads to real (positive or negative) Boltzmann weights. This enables us to apply accurate numerical methods based on the renormalization of the density matrix, namely, the corner transfer matrix renormalization group and the higher-order tensor renormalization group. For the 2D antiferromagnet, varying the imaginary magnetic field, we calculate with high accuracy the curve of critical points related to the symmetry breaking of magnetizations on the interwoven sublattices. The critical exponent $\beta$ and the anomaly number $c$ are shown to be constant along the critical line, equal to their values $\beta =\frac{1}{8}$ and $c=\frac{1}{2}$ for the 2D Ising model in a zero magnetic field. The 2D ferromagnets behave in analogy with their 1D counterparts defined on a chain of sites, namely, there exists a transient temperature which splits the temperature range into its high-temperature and low-temperature parts. The free energy and the magnetization are well defined in the high-temperature region. In the low-temperature region, the free energy exhibits singularities at the Yang-Lee zeros of the partition function and the magnetization is also ill-defined: It varies chaotically with the size of the system. The transient temperature is determined as a function of the imaginary magnetic field by using the fact that from the high-temperature side both the first derivative of the free energy with respect to the temperature and the magnetization diverge at this temperature.

We analyze energy spreading for a system that features mixed chaotic phase space, whose control parameters (or slow degrees of freedom) vary quasistatically. For demonstration purpose we consider the restricted three-body problem, where the distance between the two central stars is modulated due to their Kepler motion. If the system featured hard chaos, one would expect diffusive spreading with coefficient that can be estimated using linear-response (Kubo) theory. But for mixed phase space the chaotic sea is multilayered. Consequently, it becomes a challenge to find a robust procedure that translates the sticky dynamics into a stochastic model. We propose a Poincaré-sequencing method that reduces the multidimensional motion into a one-dimensional random walk in impact space. We test the implied relation between stickiness and the rate of spreading.

Previous works on the flight dynamic stability of insects have focused on relatively large insects. Here, the longitudinal flight dynamic stability of two hovering miniature insects was computed. With the stability properties of the miniature insects from the present work and those of large insects from previous works, we studied the effects of insect size on the stability properties in the full range of insect sizes. The following results were obtained. Although the insects considered have a 30 000-fold difference in mass, their modal structure of flight stability is the same: an unstable oscillatory mode, a stable fast subsidence mode, and a stable slow subsidence mode; because of the unstable mode, the flight is unstable. An approximate analytical expression on the growth rate of the unstable mode as a function of insect mass $\left(m\right)$ was derived. It shows that the time to double the initial values of disturbances $\left(t_d\right)$ is proportional to the 0.17 power of the insect mass $\left(m\right)$. That is, as $m$ becomes smaller, $t_d$ decreases (i.e., the instability becomes faster). This means that miniature insects need a faster nervous system to control the instability than larger insects. For example, the response time (represented by $t_d$) of a miniature insect, the gall midge $(m\approx 0.05\mathrm{mg})$, needs to be faster by about 7 times than that of a larger insect, the hawk moth $(m\approx 1600\mathrm{mg})$.

This article presents the use of artificial neural networks (ANN) to predict nonlocal heat flux transport within hydrodynamic simulations. Several cases of laser driven ablation of a plastic target are considered. The database for the ANN training phase is built using the transport module of the hydrodynamic code CHIC. It covers a range of parameters characteristic of laser experiments in the context of high-energy-density physics. Results show that an ANN can efficiently replace a module of nonlocal transport in one- and two-dimensional hydrodynamic simulations, with an error less than 3% in a radius of $0.5\mu \mathrm{m}$ and an average computation gain of a factor 433 in two dimensions.

A thermal gradient-induced circular motion of particles placed on ferroelectric nematic liquid crystal sessile drops is demonstrated and explained. Unlike hurricanes and tornadoes that are the prime examples for thermal motors and where turbulent flows are apparent, here the texture without tracer particles appears completely steady indicating laminar flow. We provide a simple model showing that the tangential arrangement of the ferroelectric polarization combined with the vertical thermal gradient and the pyroelectricity of the fluid drives the rotation of the tracer particles that become electrically charged in the fluid. These observations provide a fascinating example of the unique nature of fluid ferroelectric liquid crystals.

Epstein's model for a civil disorder is an agent-based model that simulates a social protest process where the central authority uses the police force to dissuade it. The interactions of police officers and citizens produce dynamics that do not yet have any analysis from the sociophysics approach. We present numerical simulations to characterize the properties of the one-dimensional civil disorder model on stationary state. To do this, we consider interactions on a Moore neighborhood and a random neighborhood with two different visions. We introduce a Potts-like energy function and construct the phase diagram using the agent state concentration. We find order-disorder phases and reveal the principle of minimum grievance as the underlying principle of the model's dynamics. Besides, we identify when the system can reach stable or an instability conditions based on the agents' interactions. Finally, we identified the most relevant role of the police based on their capacity to dissuade a protest and their effect on facilitating a stable scenario.

We study the ground-state properties of a quantum sunburst model, composed of a quantum Ising spin ring in a transverse field, symmetrically coupled to a set of ancillary isolated qubits, to maintain a residual translation invariance and also a ${\mathbb{Z}}_2$ symmetry. The large-size limit is taken in two different ways: either by keeping the distance between any two neighboring ancillary qubits fixed or by fixing their number while increasing the ring size. Substantially different regimes emerge, depending on the various Hamiltonian parameters: For small energy scale $\delta$ of the ancillary subsystem and small ring-qubit interaction $\kappa$, we observe rapid and nonanalytic changes in proximity to the quantum transitions of the Ising ring, both first order and continuous, which can be carefully controlled by exploiting renormalization-group and finite-size scaling frameworks. Smoother behaviors are instead observed when keeping $\delta >0$ fixed and in the Ising disordered phase. The effect of an increasing number $n$ of ancillary spins turns out to scale proportionally to $\sqrt{n}$ for sufficiently large values of $n$.

We explore a model system consisting of a particle confined to move along a toroidal helix while being exposed to a static potential as well as a driving force due to a harmonically oscillating electric field. It is shown that in the limit of a vanishing helix radius, the governing equations of motion coincide with those of the well-known Kapitza pendulum—a classical pendulum with oscillating pivot—implying that the driven toroidal helix represents a corresponding generalization. It is shown that the two dominant static fixed points present in the Kapitza pendulum are also present for a finite helix radius. The dependence of the stability of these two fixed points on the helix radius, the driving amplitude, and the static potential are analyzed analytically. These analytical results are subsequently compared to results corresponding of numerical simulations. Additionally, the most prominent deviations of the driven helix from the Kapitza pendulum with respect to the resulting phase space are investigated and analyzed in some detail. These effects include an unusual transition to chaos and an effective directed transport due to the simultaneous presence of multiple chaotic phase space regions.