The Bragg reflection band width and optical rotatory dispersion of liquid crystalline cholesteric blue phases (BPs) I and II are compared by numerical simulations. Attention is paid to the wavelength regions for which the reflection bands with lowest photon energies appear, i.e., the $\left[110\right]$ direction for BP I and the $\left[100\right]$ direction for BP II. Finite difference time domain and $4\times 4$ matrix calculations performed on the theoretical director tensor distribution of BPs with the same material parameters show that BP II, which has simple cubic symmetry, has a wider photonic band gap than BP I, which has body centered cubic symmetry, possibly due to the fact that the density of the double-twist cylinders in BP II are twice that in BP I. The theoretical results on the Bragg reflection band width are supported by reflectance measurements performed on BPs I and II for light incident along the $\left[110\right]$ and $\left[100\right]$ directions, respectively.

We present studies of chiral nematic liquid crystals composed of flexible dimer molecules subject to large dc magnetic fields between 0 and 31 T. We observe that these fields lead to selective reflection of light depending on temperature and magnetic field. The band of reflected wavelengths can be tuned from ultraviolet to beyond the IR-$C$ band. A similar effect induced by electric fields has been presented previously, and was explained by a field-induced oblique-heliconical director deformation in accordance with early theoretical predictions. The use of magnetic field here instead of electric field allows precise measurements of some material constants and holds promise for wireless tuning of selective reflection.

We investigate the stability of thin liquid curtains with respect to two-dimensional perturbations. The dynamics of perturbations with wavelengths exceeding (or comparable to) the curtain's thickness are examined using the lubrication approximation (or a kind of geometric optics). It is shown that, contrary to the previous theoretical results, but in agreement with the experimental ones, all curtains are stable with respect to small perturbations. Large perturbations can still be unstable, however, but only if they propagate upstream and, thus, disrupt the curtain at its outlet. This circumstance enables us to obtain an effective stability criterion by deriving an existence condition for upstream propagating perturbations.

The dynamical scaling analysis for the Kosterlitz-Thouless transition in the nonequilibrium relaxation method is improved by the use of Bayesian statistics and the kernel method. This allows data to be fitted to a scaling function without using any parametric model function, which makes the results more reliable and reproducible and enables automatic and faster parameter estimation. Applying this method, the bootstrap method is introduced and a numerical discrimination for the transition type is proposed.

In the lattice Boltzmann (LB) method, the forcing scheme, which is used to incorporate an external or internal force into the LB equation, plays an important role. It determines whether the force of the system is correctly implemented in an LB model and affects the numerical accuracy. In this paper we aim to clarify a critical issue about the Chapman-Enskog analysis for a class of forcing schemes in the LB method in which the velocity in the equilibrium density distribution function is given by $\mathbf{u}={\sum }_{\alpha }{\mathbf{e}}_{\alpha }{f}_{\alpha }/\rho$, while the actual fluid velocity is defined as $\widehat{\mathbf{u}}=\mathbf{u}+{\delta }_t\mathbf{F}/\left(2\rho \right)$. It is shown that the usual Chapman-Enskog analysis for this class of forcing schemes should be revised so as to derive the actual macroscopic equations recovered from these forcing schemes. Three forcing schemes belonging to the above class are analyzed, among which Wagner's forcing scheme [A. J. Wagner, Phys. Rev. E 74, 056703 (2006)] is shown to be capable of reproducing the correct macroscopic equations. The theoretical analyses are examined and demonstrated with two numerical tests, including the simulation of Womersley flow and the modeling of flat and circular interfaces by the pseudopotential multiphase LB model.

Crook's fluctuation theorem (CFT) and Jarzynski equality (JE) are effective tools for obtaining free-energy difference $\mathrm{\Delta }F({\lambda }_A\to {\lambda }_B,T_0)$ through a set of finite-time protocol driven nonequilibrium transitions between two equilibrium states $A$ and $B$ [parametrized by the time-varying protocol $\lambda \left(t\right)$] at the same temperature $T_0$. Using the generalized dimensionless work function $\mathrm{\Delta }{W}_G$, we extend CFT to transitions between two nonequilibrium steady states (NESSs) created by a thermal gradient. We show that it is possible, provided the period over which the transitions occur is sufficiently long, to obtain $\mathrm{\Delta }F({\lambda }_A\to {\lambda }_B,T_0)$ for different values of $T_0$, using the same set of finite-time transitions between these two NESSs. Our approach thus completely eliminates the need to make new samples for each new $T_0$. The generalized form of JE arises naturally as the average of the exponentiated $\mathrm{\Delta }{W}_G$. The results are demonstrated on two test cases: (i) a single particle quartic oscillator having a known closed form $\mathrm{\Delta }F$, and (ii) a one-dimensional ${\varphi }^4$ chain. Each system is sampled from the canonical distribution at an arbitrary $T^{\prime }$ with $\lambda ={\lambda }_A$, then subjected to a temperature gradient between its ends, and after steady state is reached, the protocol change ${\lambda }_A\to {\lambda }_B$ is effected in time $\tau$, following which $\mathrm{\Delta }{W}_G$ is computed. The reverse path likewise initiates in equilibrium at $T^{\prime }$ with $\lambda ={\lambda }_B$ and the protocol is time reversed leading to $\lambda ={\lambda }_A$ and the reverse $\mathrm{\Delta }{W}_G$. Our method is found to be more efficient than either JE or CFT when free-energy differences at multiple $T_0$'s are required for the same system.

Network reliability is the probability that a dynamical system composed of discrete elements interacting on a network will be found in a configuration that satisfies a particular property. We introduce a reliability property, Ising feasibility, for which the network reliability is the Ising model's partition function. As shown by Moore and Shannon, the network reliability can be separated into two factors: structural, solely determined by the network topology, and dynamical, determined by the underlying dynamics. In this case, the structural factor is known as the joint density of states. Using methods developed to approximate the structural factor for other reliability properties, we simulate the joint density of states, yielding an approximation for the partition function. Based on a detailed examination of why naïve Monte Carlo sampling gives a poor approximation, we introduce a parallel scheme for estimating the joint density of states using a Markov-chain Monte Carlo method with a spin-exchange random walk. This parallel scheme makes simulating the Ising model in the presence of an external field practical on small computer clusters for networks with arbitrary topology with $\sim {10}^6$ energy levels and more than ${10}^{308}$ microstates.

Earthquakes are quite hard to predict. One of the possible reasons can be the fact that the existing catalogs of past earthquakes are limited at most to the order of 100 years, while their characteristic time scale is sometimes greater than that time span. Here we rather use these limitations positively and characterize some large earthquake events as abnormal events that are not included there. When we constructed probabilistic forecasts for large earthquakes in Japan based on similarity and difference to their past patterns—which we call known and unknown abnormalities, respectively—our forecast achieved probabilistic gains of 5.7 and 2.4 against a time-independent model for main shocks with the magnitudes of 7 or above. Moreover, the two abnormal conditions covered 70% of days whose maximum magnitude was 7 or above.

The highly variable dynamics of neocortical circuits observed in vivo have been hypothesized to represent a signature of ongoing stochastic inference but stand in apparent contrast to the deterministic response of neurons measured in vitro. Based on a propagation of the membrane autocorrelation across spike bursts, we provide an analytical derivation of the neural activation function that holds for a large parameter space, including the high-conductance state. On this basis, we show how an ensemble of leaky integrate-and-fire neurons with conductance-based synapses embedded in a spiking environment can attain the correct firing statistics for sampling from a well-defined target distribution. For recurrent networks, we examine convergence toward stationarity in computer simulations and demonstrate sample-based Bayesian inference in a mixed graphical model. This points to a new computational role of high-conductance states and establishes a rigorous link between deterministic neuron models and functional stochastic dynamics on the network level.

We measure the vibrational modes and particle dynamics of quasi-two-dimensional colloidal glasses as a function of interparticle interaction strength. The interparticle attractions are controlled via a temperature-tunable depletion interaction. Specifically, the interparticle attraction energy is increased gradually from a very small value (nearly hard-sphere) to moderate strength $(\sim 4k_{B}T)$, and the variation of colloidal particle dynamics and vibrations are concurrently probed. The particle dynamics slow monotonically with increasing attraction strength, and the particle motions saturate for strengths greater than $\sim 2k_{B}T$, i.e., as the system evolves from a nearly repulsive glass to an attractive glass. The shape of the phonon density of states is revealed to change with increasing attraction strength, and the number of low-frequency modes exhibits a crossover for glasses with weak compared to strong interparticle attraction at a threshold of $\sim 2k_{B}T$. This variation in the properties of the low-frequency vibrational modes suggests a new means for distinguishing between repulsive and attractive glass states.

In this work, with focus on the energy-transport properties in quantum, low-dimensional, graded materials, we address the investigation of the energy (and spin) current in $XXZ$ open chains with graded inner structures and driven out of equilibrium by magnetization pumping applied at the ends. We study several types of graded structures in different situations in order to show a ubiquitous occurrence of energy rectification, even for the system under a homogeneous magnetic field. Due to technical difficulties, we carry out the computation for small chains, but we present arguments that indicate the extension of some results to larger systems. Recalling the generic existence of energy rectification in classical, graded materials, which are described by anharmonic chains of oscillators, and recalling also the anharmonicity of these $XXZ$ models, which involve quartic terms in more transparent representation in terms of fermionic creation and annihilation operators, we may say that our results extend the ubiquity of energy rectification occurrence in classical graded materials to the case of quantum systems.

Recently it was stated that temporal disorder constitutes a relevant perturbation in absorbing phase transitions for all dimensions. However, its effect on systems other than the standard contact process (CP), its competition with other ingredients (e.g., particle diffusion), and other kinds of disorder (besides the standard types) are unknown. In order to shed some light on the above-mentioned points, we investigate a variant of the usual CP, namely, the triplet annihilation model, in which the competition between triplet annihilation and single particle diffusion leads to an unusual phase diagram behavior, with reentrant shape and endless activity for sufficiently large diffusion rates. Two kinds of time-dependent disorder have been considered. In the former, it is introduced in the creation-annihilation parameters (as commonly considered in recent studies), whereas in the latter, the diffusion rate $D$ is allowed to be time dependent. In all cases, the disorder follows a uniform distribution with fixed mean and width $\sigma$. Two values of $\sigma$ have been considered in order to exemplify the regime of “weaker” and “stronger” temporal disorder strengths. Our results show that in the former approach, the disorder suppresses the reentrant phase diagram with a critical behavior deviating from the directed percolation (DP) universality class in the regime of low diffusion rates, while they strongly suggest that the DP class is recovered for larger hopping rates. An opposite scenario is found in the latter disorder approach, with a substantial increase of reentrant shape and the maximum diffusion, in which the reentrant shape also displays a critical behavior consistent with the DP universality class (in similarity with the pure model). In order to compare with very recent claims, the results from taking a bimodal distribution and critical behavior in the limit of strong disorder are presented. Also, the results derived from the mean-field theory are performed, presenting partial agreement with numerical results. Lastly, a comparison with the diffusive disordered CP is undertaken.

Grand-canonical fluctuations of Bose-Einstein condensates of light are accessible to state-of-the-art experiments [J. Schmitt et al., Phys. Rev. Lett. 112, 030401 (2014).]. We phenomenologically describe these fluctuations by using the grand-canonical ensemble for a weakly interacting Bose gas at thermal equilibrium. For a two-dimensional harmonic trap, we use two models for which the canonical partition functions of the weakly interacting Bose gas are given by exact recurrence relations. We find that the grand-canonical condensate fluctuations for weakly interacting Bose gases vanish at zero temperature, thus behaving qualitatively similarly to an ideal gas in the canonical ensemble (or microcanonical ensemble) rather than the grand-canonical ensemble. For low but finite temperatures, the fluctuations remain considerably higher than for the canonical ensemble, as predicted by the ideal gas in the grand-canonical ensemble, thus clearly showing that we are not in a regime in which the ensembles are equivalent.

Under external mechanical loading, glassy materials, ranging from soft matter systems to metallic alloys, often respond via formation of inhomogeneous flow patterns, during yielding. These inhomogeneities can be precursors to catastrophic failure, implying that a better understanding of their underlying mechanisms could lead to the design of smarter materials. Here, extensive molecular dynamics simulations are used to reveal the emergence of heterogeneous dynamics in a binary Lennard-Jones glass, subjected to a constant strain rate. At a critical strain, this system exhibits for all considered strain rates a transition towards the formation of a percolating cluster of mobile regions. We give evidence that this transition belongs to the universality class of directed percolation. Only at low shear rates, the percolating cluster subsequently evolves into a transient (but long-lived) shear band with a diffusive growth of its width. Finally, the steady state with a homogeneous flow pattern is reached. In the steady state, percolation transitions also do occur constantly, albeit over smaller strain intervals, to maintain the stationary plastic flow in the system.

We study the behavior of high-energy positrons emitted from a radioactive source in a magnetospheric dipole field configuration. Because the conservation of the first and second adiabatic invariants is easily destroyed in a strongly inhomogeneous dipole field for high-energy charged particles, the positron orbits are nonintegrable, resulting in chaotic motions. In the geometry of a typical magnetospheric levitated dipole experiment, it is shown that a considerable ratio of positrons from a ${}^{22}\mathrm{Na}$ source, located at the edge of the confinement region, has chaotic long orbit lengths before annihilation. These particles make multiple toroidal circulations and form a hollow toroidal positron cloud. Experiments with a small ${}^{22}\mathrm{Na}$ source in the Ring Trap 1 (RT-1) device demonstrated the existence of such long-lived positrons in a dipole field. Such a chaotic behavior of high-energy particles is potentially applicable to the formation of a dense toroidal positron cloud in the strong-field region of the dipole field in future studies.

We show that solutions to fermion sign problems that are found in the formulation where the path integral is expanded in powers of the interaction in continuous time can be extended to systems involving fermions interacting with dynamical quantum spins. While these sign problems seem unsolvable in the auxiliary field approach, solutions emerge in the world-line representation of quantum spins. Combining this idea with meron-cluster methods, we are able to further extend the class of models that are solvable. We demonstrate these solutions to sign problems by considering several examples of strongly correlated systems that contain the physics of semimetals, insulators, superfluidity, and antiferromagnetism.

Most self-sustained oscillations in biological systems and in technical applications are based on a feedback loop, and it is usually important to know how they will react when an external oscillatory force is applied. Here we investigate the effects that the introduction of a time delay in the feedback can have in the entrainment properties of self-sustained oscillators. To do this, we derive analytic expressions for the periodic trajectories and their asymptotic stability, for a generic external oscillatory force. This allows us to show that, for large quality factors, the resonance frequency does not depend on the feedback delay. When the external force is harmonic, it is shown that the largest entrainment range does not correspond to the time delay that gives the maximal response of the unforced oscillator. In fact, that delay gives the shortest entrainment range.

The heterogeneity of coupled oscillators is important for the degree of their synchronization. According to the classical Kuramoto model, larger heterogeneity reduces synchronization. Here, we show that in a model for coupled pancreatic $\beta$-cells, higher diversity of the cells induces higher synchrony. We find that any system of coupled oscillators that oscillates on two time scales and in which heterogeneity causes a transition from chaotic to damped oscillations on the fast time scale exhibits this property. Thus, synchronization of a subset of oscillating systems can be enhanced by increasing the heterogeneity of the system constituents.

State-space reconstruction is investigated for evaluating the randomness generated by an optically injected semiconductor laser in chaos. The reconstruction of the attractor requires only the emission intensity time series, allowing both experimental and numerical evaluations with good qualitative agreement. The randomness generation is evaluated by the divergence of neighboring states, which is quantified by the time-dependent exponents (TDEs) as well as the associated entropies. Averaged over the entire attractor, the mean TDE is observed to be positive as it increases with the evolution time through chaotic mixing. At a constant laser noise strength, the mean TDE for chaos is observed to be greater than that for periodic dynamics, as attributed to the effect of noise amplification by chaos. After discretization, the Shannon entropies continually generated by the laser for the output bits are estimated in providing a fundamental basis for random bit generation, where a combined output bit rate reaching 200 Gb/s is illustrated using practical tests. Overall, based on the reconstructed states, the TDEs and entropies offer a direct experimental verification of the randomness generated in the chaotic laser.

The motion of a self-propelled particle is affected by its surroundings, such as boundaries or external fields. In this paper, we investigated the bifurcation of the motion of a camphor grain, as a simple actual self-propelled system, confined in a one-dimensional finite region. A camphor grain exhibits oscillatory motion or remains at rest around the center position in a one-dimensional finite water channel, depending on the length of the water channel and the resistance coefficient. A mathematical model including the boundary effect is analytically reduced to an ordinary differential equation. Linear stability analysis reveals that the Hopf bifurcation occurs, reflecting the symmetry of the system.