Sharp two- and three-dimensional phase transitional magnetization curves are obtained by an iterative renormalization-group coupling of Ising chains, which are solved exactly. The chains by themselves do not have a phase transition or nonzero magnetization, but the method reflects crossover from temperaturelike to fieldlike renormalization-group flows as the mechanism for the higher-dimensional phase transitions. The magnetization of each chain acts, via the interaction constant, as a magnetic field on its neighboring chains, thus entering its renormalization-group calculation. The method is highly flexible for wide application.

Spontaneous pattern formation is a fundamental scientific problem that has received much attention since the seminal theoretical work of Turing on reaction-diffusion systems. In molecular biophysics, this phenomenon often takes place under the influence of large fluctuations. It is then natural to inquire about the precision of such pattern. In particular, spontaneous pattern formation is a nonequilibrium phenomenon, and the relation between the precision of a pattern and the thermodynamic cost associated with it remains largely unexplored. Here, we analyze this relation with a paradigmatic stochastic reaction-diffusion model, i.e., the Brusselator in one spatial dimension. We find that the precision of the pattern is maximized for an intermediate thermodynamic cost, i.e., increasing the thermodynamic cost beyond this value makes the pattern less precise. Even though fluctuations get less pronounced with an increase in thermodynamic cost, we argue that larger fluctuations can also have a positive effect on the precision of the pattern.

We study the double ionization of atoms subjected to circularly polarized (CP) laser pulses. We analyze two fundamental ionization processes: the sequential (SDI) and nonsequential (NSDI) double ionization in the light of the rotating frame (RF) which naturally embeds nonadiabatic effects in CP pulses. We use and compare two adiabatic approximations: The adiabatic approximation in the laboratory frame (LF) and the adiabatic approximation in the RF. The adiabatic approximation in the RF encapsulates the energy variations of the electrons on subcycle timescales happening in the LF and this, by fully taking into account the ion-electron interaction. This allows us to identify two nonadiabatic effects including the lowering of the threshold intensity at which over-the-barrier ionization happens and the lowering of the ionization time of the electrons. As a consequence, these nonadiabatic effects facilitate over-the-barrier ionization and recollision-induced ionizations. We analyze the outcomes of these nonadiabatic effects on the recollision mechanism. We show that the laser envelope plays an instrumental role in a recollision channel in CP pulses at the heart of NSDI.

Chaotic dynamics of a dynamical system is not necessarily persistent. If there is (without any active intervention from outside) a transition towards a (possibly nonchaotic) attractor, this phenomenon is called transient chaos, which can be observed in a variety of systems, e.g., in chemical reactions, population dynamics, neuronal activity, or cardiac dynamics. Also, chimera states, which show coherent and incoherent dynamics in spatially distinct regions of the system, are often chaotic transients. In many practical cases, the control of the chaotic dynamics (either the termination or the preservation of the chaotic dynamics) is desired. Although the self-termination typically occurs quite abruptly and can so far in general not be properly predicted, previous studies showed that in many systems a ‘terminal transient phase” (TTP) prior to the self-termination existed, where the system was less susceptible against small but finite perturbations in different directions in state space. In this study, we show that, in the specific case of chimera states, these susceptible directions can be related to the structure of the chimera, which we divide into the coherent part, the incoherent part and the boundary in between. That means, in practice, if self-termination is close we can identify the direction of perturbation which is likely to maintain the chaotic dynamics (the chimera state). This finding improves the general understanding of the state space structure during the TTP, and could contribute also to practical applications like future control strategies of epileptic seizures which have been recently related to the collapse of chimera states.

The origins of the large differences observed in the rates at which diverse particles are conveyed along axonal microtubules are still a matter of debate in the literature. There is evidence that certain neurodegenerative diseases may be triggered by disturbances in the related transport processes. Motivated by this, we employ a model to investigate mobility properties of certain cargoes whose dynamics are coupled with that of molecular motors on crowded microtubules. For certain initial and boundary conditions, we use the method of characteristics to resolve perturbatively the pair of equations of Burgers type resulting from a mean-field approach to the original microscopic stochastic model. Extensions to nonperturbative limits are explored numerically. In this context, we are able to figure out conditions under which the cargoes' average velocities may differ up to orders of magnitude just by changing the number of motors on the considered track. We then discuss possibilities to connect these theoretical predictions with available experimental data about axon transport.

Spiral waves of excitation are common in many physical, chemical, and biological systems. In physiological systems like the heart, such waves can lead to cardiac arrhythmias and need to be eliminated. Spiral waves anchor to heterogeneities in the excitable medium, and to eliminate them they need to be unpinned first. Several groups focused on developing strategies to unpin such pinned waves using electric shocks, pulsed electric fields, and recently, circularly polarized electric fields (CPEF). It was shown that in many situations, CPEF is more efficient at unpinning the wave compared to other existing methods. Here, we study how the circularly polarized field acts on the pinned spiral waves and unpins it. We show that the termination always happens within the first rotation of the electric field. For a given obstacle size, there exists a threshold time period of the CPEF below which the spiral can always be terminated. Our analytical formulation accurately predicts this threshold and explains the absence of the traditional unpinning window with the CPEF. We hope our theoretical work will stimulate further experimental studies about CPEF and low energy methods to eliminate spiral waves.

Two dynamical systems unidirectionally coupled in a sender-receiver configuration can synchronize with a nonzero phase lag. In particular, the system can exhibit anticipated synchronization (AS), which is characterized by a negative phase lag, if the receiver also receives a delayed negative self-feedback. Recently, AS was shown to occur between cortical-like neuronal populations in which the self-feedback is mediated by inhibitory synapses. In this biologically plausible scenario, a transition from the usual delayed synchronization (with positive phase lag) to AS can be mediated by the inhibitory conductances in the receiver population. Here we show that depending on the relation between excitatory and inhibitory synaptic conductances the system can also exhibit phase bistability between anticipated and delayed synchronization. Furthermore, we show that the amount of noise at the receiver and the synaptic conductances can mediate the transition from stable phase locking to a bistable regime and eventually to a phase drift. We suggest that our spiking neuronal populations model could be potentially useful to study phase bistability in cortical regions related to bistable perception.

The dynamics of magnetization relaxation in ferrofluids are studied with statistical-mechanical theory and Brownian dynamics simulations. The particle dipole moments are initially perfectly aligned, and the magnetization is equal to its saturation value. The magnetization is then allowed to decay under zero-field conditions toward its equilibrium value of zero. The time dependence is predicted by solving the Fokker-Planck equation for the one-particle orientational distribution function. Interactions between particles are included by introducing an effective magnetic field acting on a given particle and arising from all of the other particles. Two different approximations are proposed and tested against simulations: a first-order modified mean-field theory and a modified Weiss model. The theory predicts that the short-time decay is characterized by the Brownian rotation time ${\tau }_B$, independent of the interaction strength. At times much longer than ${\tau }_B$, the asymptotic decay time is predicted to grow with increasing interaction strength. These predictions are borne out by the simulations. The modified Weiss model gives the best agreement with simulation, and its range of validity is limited to moderate, but realistic, values of the dipolar coupling constant.

We investigate the dissipative mechanisms exhibited by creased material sheets when subjected to mechanical loading, which comes in the form of plasticity and relaxation phenomena within the creases. After demonstrating that plasticity mostly affects the rest angle of the creases, we devise a mapping between this quantity and the macroscopic state of the system that allows us to track its reference configuration along an arbitrary loading path, resulting in a powerful monitoring and design tool for crease-based metamaterials. Furthermore, we show that complex relaxation phenomena, in particular memory effects, can give rise to a nonmonotonic response at the crease level, possibly relating to the similar behavior reported for crumpled sheets. We describe our observations through a classical double-logarithmic time evolution and obtain a constitutive behavior compatible with that of the underlying material. Thus the lever effect provided by the crease allows magnified access to the material's rheology.

In this paper we compare different theoretical approaches to describe the dispersion of collective modes in Yukawa fluids when the interparticle coupling is relatively weak, so that the kinetic and potential contributions to the dispersion relation compete with each other. A thorough comparison with the results from molecular dynamics simulation allows us to conclude that, in the investigated regime, the best description is provided by the sum of the generalized excess bulk modulus and the Bohm-Gross kinetic term.

Heat flow management at the nanoscale is of great importance for emergent quantum technologies. For instance, a thermal sink that can be activated on-demand is a highly desirable tool that may accommodate the need to evacuate excess heat at chosen times, e.g., to maintain cryogenic temperatures or reset a quantum system to ground, and the possibility of controlled unitary evolution otherwise. Here we propose a design of such heat switch based on a single coherently driven qubit. We show that the heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. The complete suppression of the heat flow is a quantum effect occurring for specific driving parameters that we express and we analyze the role of the coherences in the free-qubit energy eigenbasis. We finally study the feasibility of this quantum heat switch in a circuit QED setup involving a charge qubit coupled to thermal resistances. We demonstrate robustness to experimental imperfections such as additional decoherence, paving the road towards experimental verification of this effect.

We investigate the phase transition of the dodecahedron model on the square lattice. The model is a discrete analog of the classical Heisenberg model, which has continuous $O\left(3\right)$ symmetry. In order to treat the large on-site degree of freedom $q=20$, we develop a massively parallelized numerical algorithm for the corner transfer matrix renormalization group method, incorporating EigenExa, the high-performance parallelized eigensolver. The scaling analysis with respect to the cutoff dimension reveals that there is a second-order phase transition at $T_{\mathrm{c}}^{}=0.4398\left(8\right)$ with the critical exponents $\nu =2.88\left(8\right)$ and $\beta =0.21\left(1\right)$. The central charge of the system is estimated as $c=1.99\left(6\right)$.

The statistics of chiral matrix ensembles with uncorrelated but multivariate Gaussian distributed elements is intuitively expected to be driven by many parameters. Contrary to intuition, however, our theoretical analysis reveals the existence of a single parameter, a function of all ensemble parameters, which governs the dynamics of spectral statistics. The analysis not only extends the formulation (known as complexity parameter formulation) for Hermitian ensembles without chirality to those with it but also reveals the underlying connection between chiral complex systems with seemingly different system conditions as well as between other complex systems, e.g., multiparametric Wishart ensembles as well as generalized Calogero-Sutherland Hamiltonians.

The statistics of self-avoiding random walks have been used to model polymer physics for decades. A self-avoiding walk that grows one step at a time on a lattice will eventually trap itself, which occurs after an average of 71 steps on a square lattice. Here, we consider the effect of nearest-neighbor attractive interactions on isolated growing self-avoiding walks, and we examine the effect that self-attraction has both on the statistics of trapping as well as on chain statistics through the transition between expanded and collapsed walks at the theta point. We find that the trapping length increases exponentially with the nearest-neighbor contact energy, but that there is a local minimum in trapping length for weakly self-attractive walks. While it has been controversial whether growing self-avoiding walks have the same asymptotic behavior as traditional self-avoiding walks, we find that the theta point is not at the same location for growing self-avoiding walks, and that the persistence length converges much more rapidly to a smaller value.

Understanding the functional connectivity of the brain has become a major goal of neuroscience. In many situations the relative phase difference, together with coherence patterns, has been employed to infer the direction of the information flow. However, it has been recently shown in local field potential data from monkeys the existence of a synchronized regime in which unidirectionally coupled areas can present both positive and negative phase differences. During the counterintuitive regime, called anticipated synchronization (AS), the phase difference does not reflect the causality. Here we investigate coherence and causality at the alpha frequency band ($f\sim 10$ Hz) between pairs of electroencephalogram (EEG) electrodes in humans during a GO/NO-GO task. We show that human EEG signals can exhibit anticipated synchronization, which is characterized by a unidirectional influence from an electrode A to an electrode B, but the electrode B leads the electrode A in time. To the best of our knowledge, this is the first verification of AS in EEG signals and in the human brain. The usual delayed synchronization (DS) regime is also present between many pairs. DS is characterized by a unidirectional influence from an electrode A to an electrode B and a positive phase difference between A and B which indicates that the electrode A leads the electrode B in time. Moreover we show that EEG signals exhibit diversity in the phase relations: the pairs of electrodes can present in-phase, antiphase, or out-of-phase synchronization with a similar distribution of positive and negative phase differences.

The dynamics of pseudo-spin-1/2 Bose-Einstein condensates with weak spin-orbit coupling through a moving obstacle potential are studied numerically. Four types of wakes are observed and the phase diagrams are determined for different spin-orbit coupling strengths. The conditions to form Bénard–von Kármán vortex street are rather rigorous, and we investigate in detail the dynamical characteristics of the vortex streets. The two point vortices in a pair rotate around their center, and the angular velocity and their distance oscillate periodically. The oscillation intensifies with increasing spin-orbit coupling strengths, and it makes part of the vortex pairs dissociate into separate vortices or combine into single ones and destroys the vortex street in the end. The width $b$ of the street and the distance $l$ between two consecutive vortex pairs of the same circulation are determined by the potential radius and its moving velocity, respectively. The $b/l$ ratios are independent of the spin-orbit coupling strength and fall in the range $0.19–0.27$, which is a little smaller than the stability criterion 0.28 for classical fluids. Proper $b/l$ ratios are necessary to form Bénard–von Kármán vortex street, but the spin-orbit coupling strength affects the stability of the street patterns. Finally, we propose a protocol to experimentally realize the vortex street in ${}^{87}\mathrm{Rb}$ spin-orbit-coupling Bose-Einstein condensates.

We propose a tractable epidemic model that includes containment measures. In the absence of containment measures, the epidemics spread exponentially fast whenever the infectivity rate is positive $\lambda >0$. The containment measures are modeled by considering a time-dependent modulation of the bare infectivity $\lambda$ leading to effective infectivity that decays in time for each infected individual, mimicking, for instance, the combined effect of the asymptomatic onset of the disease, testing policies, and quarantine. We consider a wide range of temporal kernels for effective infectivity, and we investigate the effect of the considered containment measures. We find that not all kernels are able to push the epidemic dynamics below the epidemic threshold with some containment measures only able to reduce the rate of the exponential growth of newly infected individuals. We also propose a pandemic model caused by a growing number of separated foci.

The purpose of this work is to derive a small turbulent Péclet–small turbulent Mach number approximation for hydroradiative turbulent mixing zones encountered in stellar interiors where the radiative conductivity can overwhelms the turbulent transport. To this end, we proceed to an asymptotic analysis and determine orders of magnitude for the fluctuating temperature and pressure, as well as closed expressions for the fluctuating conduction and velocity divergence. The latter is used to extend a Reynolds stress model to the small-Péclet regime. Three-dimensional direct numerical simulations of radiative Rayleigh-Taylor turbulent mixing zones are performed, first, to validate the asymptotic predictions and, second, to validate their use in the Reynolds stress model.

Nowadays hydrodynamic equations coupled with external equation of states provided by quantum mechanical calculations is a widely used approach for simulations of macroscopic degenerate plasmas. Although such an approach is proven to be efficient and shows many good features, especially for large scale simulations, it encounters intrinsic challenges when involving kinetic effects. As a complement, here we have invented a fully kinetic numerical approach for macroscopic degenerate plasmas. This approach is based on first principle Boltzmann-Uhling-Uhlenbeck equations coupled with Maxwell's equation, and is eventually achieved via an existing particle-in-cell simulation code named LAPINS. In this approach, degenerate particles obey Fermi-Dirac statistics and nondegenerate particles follow the typical Maxwell-Boltzmann statistics. The equation of motion of both degenerate and nondegenerate particles are governed by long range collective electromagnetic fields and close particle-particle collisions. Especially, Boltzmann-Uhling-Uhlenbeck collisions ensure that evolution of degenerate particles is enforced by the Pauli exclusion principle. The code is applied to several benchmark simulations, including electronic conductivity for aluminium with varying temperatures from 2 eV to 50 eV, thermalization of alpha particles in a cold fuel shell in inertial confinement fusion, and rapid heating of solid sample by short and intense laser pulses.

In this paper, we study quantum phase transitions and magnetic properties of a one-dimensional spin-1/2 Gamma model, which describes the off-diagonal exchange interactions between edge-shared octahedra with strong spin-orbit couplings along the sawtooth chain. The competing exchange interactions between the nearest neighbors and the second neighbors stabilize the semimetallic ground state in terms of spinless fermions, and give rise to a rich phase diagram, which consists of three gapless phases. We find distinct phases are characterized by the number of Weyl nodes in the momentum space, and such changes in the topology of the Fermi surface without symmetry breaking produce a variety of Lifshitz transitions, in which the Weyl nodes situating at $k=\pi$ change from type I to type II. A coexistence of type-I and type-II Weyl nodes is found in phase II. The information measures including concurrence, entanglement entropy, and relative entropy can effectively signal the second-order transitions. The results indicate that the Gamma model can act as an exactly solvable model to describe Lifshitz phase transitions in correlated electron systems.