We investigate experimentally the formation of focal conic domains of the ferroelectric phase of a liquid crystal, chiral smectic $C$ ($\mathrm{Sm}{C}^{*}$), in the meniscus geometry. The meniscus geometry is formed in the gap between two glass plates which are placed on a common substrate. This gap is called here a physical cavity. Focal conic domains (FCDs) in the physical cavity with dimensions of micrometer scale are investigated under an optical polarizing microscope which enables us to extract the information on the helical structure formation in the constraint and gradient topological meniscus interface. The helical pitch in the FCD is observed to be shorter than in planar confined geometry. A crucial phenomenon of unwrapping and wrapping of helical structure from one FCD to another is also observed. In-plane application of an electric field on a FCD revealed the asymmetric helical unwinding process whereas an increase in temperature has shown symmetrical unwinding. The helical structure based observation is significant for understanding the ferroelectric phase in focal conic domains and their application in microlenses and optical components.

We investigate numerically the structure of twin boundaries of cholesteric blue phases. Our study is based on the Landau–de Gennes continuum theory describing the orientational order of the liquid crystal by a second-rank tensor. We pay particular attention to blue phase I (BP I) with body-centered-cubic symmetry and consider twin boundaries between BP I lattices in which their (110) planes are shared and the ($\overline{1}12$) plane of one lattice is parallel to the ($1\overline{1}2$) plane of the other as observed in previous experiments [Jin et al., Sci. Adv. 6, eaay5986 (2020); Zhang et al., ACS Appl. Mater. Interf. 13, 36130 (2021)]. We discuss two plausible cases in which the twin boundaries are parallel to the $\left\{112\right\}$ planes or the $\left\{111\right\}$ planes. In the former, disclination lines of obtusely bent form penetrate the twin boundaries, and in the latter straight disclination lines as well as bent ones are found at the twin boundaries. The former twin boundaries are energetically less costly, consistent with previous experimental identifications. From our numerical results the free energy of a twin boundary per unit area is estimated to be $\simeq 4\times {10}^{-6}{\mathrm{J}\mathrm{m}}^{-2}$, which indeed indicates that the formation of twin boundaries is not prohibitively costly.

In amorphous solids, plastic flow is prone to localization into shear bands via an avalanche of shear-transformation (ST) rearrangements of constituent atoms or particles. However, such banding instability still remains a lack of direct experimental evidence. Using a real 3D colloidal glass under shear as proof of principle, we study STs’ avalanches into shear banding that is controlled by strain rates. We demonstrate that, accompanying the emergent shear banding, the elastic response fields of the system, typical of a quadrupole for shear and a centrosymmetry for dilatation, lose the Eshelby-type spatial symmetry; instead, a strong correlation appears preferentially along the banding direction. By quantifying the fields’ spatial decay, we identify an elastic criterion for the shear-banding instability, that is, the strongly correlated length of dilatation is smaller than the full length of shear correlation. Specifically, ST-induced free volume has to be confined within the elastic shear domain of ST so that those STs can self-organize to trigger shear banding. This physical picture is directly visualized by tracing the real-space evolution of local dilatation and ST particles. The present work unites the two classical mechanisms: free volume and STs, for the fundamental understanding of shear banding in amorphous solids.

A computational technique based on the lattice Boltzmann method (LBM) is developed to simulate the wettable particles adsorbed to a liquid-vapor interface under gravity. The proposed technique combines the improved smoothed-profile LBM for the treatment of moving solid particles in a fluid and the free-energy LBM for the description of a liquid-vapor system. Five benchmark two-dimensional problems are examined: (A) a stationary liquid drop in the vapor phase; a wettable particle adsorbed to a liquid-vapor interface in (B) the absence and (C) the presence of gravity; (D) two freely moving particles at a liquid-vapor interface in the presence of gravity (i.e., capillary flotation forces); and (E) two vertically constrained particles at a liquid-vapor interface (i.e., capillary immersion forces). The simulation results are in good quantitative agreement with theoretical estimations, demonstrating that the proposed technique can reproduce the capillary interactions between wettable particles at a liquid-vapor interface under gravity.

Entropy production plays a fundamental role in the study of nonequilibrium systems by offering a quantitative handle on the degree of time-reversal symmetry breaking. It depends crucially on the degree of freedom considered as well as on the scale of description. How the entropy production at one resolution of the degrees of freedom is related to the entropy production at another resolution is a fundamental question which has recently attracted interest. This relationship is of particular relevance to coarse-grained and continuum descriptions of a given phenomenon. In this work, we derive the scaling of the entropy production under iterative coarse graining on the basis of the correlations of the underlying microscopic transition rates for noninteracting particles in active disordered media. Our approach unveils a natural criterion to distinguish equilibrium-like and genuinely nonequilibrium macroscopic phenomena based on the sign of the scaling exponent of the entropy production per mesostate.

Discretization of continuous stochastic processes is needed to numerically simulate them or to infer models from experimental time series. However, depending on the nature of the process, the same discretization scheme may perform very differently for the two tasks, if it is not accurate enough. Exact discretizations, which work equally well at any scale, are characterized by the property of invariance under coarse-graining. Motivated by this observation, we build an explicit renormalization group (RG) approach for Gaussian time series generated by autoregressive models. We show that the RG fixed points correspond to discretizations of linear SDEs, and only come in the form of first order Markov processes or non-Markovian ones. This fact provides an alternative explanation of why standard delay-vector embedding procedures fail in reconstructing partially observed noise-driven systems. We also suggest a possible effective Markovian discretization for the inference of partially observed underdamped equilibrium processes based on the exploitation of the Einstein relation.

Diffusion and first passage in the presence of stochastic resetting and potential bias have been of recent interest. We study a few models, systematically progressing in their complexity, to understand the usefulness of resetting. In the parameter space of the models, there are multiple continuous and discontinuous transitions where the advantage of resetting vanishes. We show these results analytically exactly for a tent potential, and numerically accurately for a quartic potential relevant to a magnetic system at low temperatures. We find that the spatial asymmetry of the potential across the barrier, and the number of absorbing boundaries, play a crucial role in determining the type of transition.

In the present work, we revisit a recently proposed and experimentally realized topological two-dimensional lattice with periodically time-dependent interactions. We identify the fundamental solitons, previously observed in experiments and direct numerical simulations, as exact, exponentially localized, periodic in time solutions. This is done for a variety of phase-shift angles of the central nodes upon an oscillation period of the coupling strength. Subsequently, we perform a systematic Floquet stability analysis of the relevant structures. We analyze both their point and their continuous spectrum and find that the solutions are generically stable, aside from the possible emergence of complex quartets due to the collision of bands of continuous spectrum. The relevant instabilities become weaker as the lattice size gets larger. Finally, we also consider multisoliton analogs of these Floquet states, inspired by the corresponding discrete nonlinear Schrödinger (DNLS) lattice. When exciting initially multiple sites in phase, we find that the solutions reflect the instability of their DNLS multi-soliton counterparts, while for configurations with multiple excited sites in alternating phases, the Floquet states are spectrally stable, again analogously to their DNLS counterparts.

Centrality metrics aim to identify the most relevant nodes in a network. In the literature, a broad set of metrics exists, measuring either local or global centrality characteristics. Nevertheless, when networks exhibit a high spectral gap, the usual global centrality measures typically do not add significant information with respect to the degree, i.e., the simplest local metric. To extract different information from this class of networks, we propose the use of the Generalized Economic Complexity index (GENEPY). Despite its original definition within the economic field, the GENEPY can be easily applied and interpreted on a wide range of networks, characterized by high spectral gap, including monopartite and bipartite network systems. Tests on synthetic and real-world networks show that the GENEPY can shed light about the node centrality, carrying information generally poorly correlated with the node number of direct connections (node degree).

Recent publications on spectroscopy of water layers in water bridge structures revealed a significant enhancement of the proton mobility and the dielectric contribution of translational vibrations of water molecules in the interfacial layers compared to bulk water. Herewith, the results of long-term studies of proton dynamics in solid-state acids have shown that proton mobility increases significantly with the predominance of hydronium, but not Zundel, cations in the aqueous phase. In the present work, in the light of these data, we reanalyzed our previously published results on broadband dielectric spectroscopy of bovine heart cytochrome c, bovine serum albumin, and the extracellular matrix and filaments of Shewanella oneidensis MR-1. We revealed that, just as in water bridges, an increase in electrical conductivity in these systems correlates with an increase in the dielectric contribution of water molecular translational vibrations. In addition, the appearance of spectral signatures of the hydronium cations was observed only in those cases when the system revealed noticeable electrical conductivity due to delocalized charge carriers.

Path integrals with complex actions are encountered for many physical systems ranging from spin- or mass-imbalanced atomic gases and graphene to quantum chromodynamics at finite density to the nonequilibrium evolution of quantum systems. Many computational approaches have been developed for tackling the sign problem emerging for complex actions. Among these, complex Langevin dynamics has the appeal of general applicability. One of its key challenges is the potential convergence of the dynamics to unphysical fixed points. The statistical sampling process at such a fixed point is not based on the physical action and hence leads to wrong predictions. Moreover, its unphysical nature is hard to detect due to the implicit nature of the process. In the present work we set up a general approach based on a Markov chain Monte Carlo scheme in an extended state space. In this approach we derive an explicit real sampling process for generalized complex Langevin dynamics. Subject to a set of constraints, this sampling process is the physical one. These constraints originate from the detailed-balance equations satisfied by the Monte Carlo scheme. This allows us to rederive complex Langevin dynamics from a new perspective and establishes a framework for the explicit construction of new sampling schemes for complex actions.

We mimic nanorod-based transparent electrodes as random resistor networks (RRNs) produced by the homogeneous, isotropic, and random deposition of conductive zero-width sticks onto an insulating substrate. We suppose that the number density (the number of objects per unit area of the surface) of these sticks exceeds the percolation threshold, i.e., the system under consideration is a conductor. We computed the electrical conductivity of random resistor networks versus the number density of conductive fillers for the wire-resistance-dominated case, for the junction-resistance-dominated case, and for an intermediate case. We also offer a consistent continuous variant of the mean-field approach. The results of the RRN computations were compared with this mean-field approach. Our computations suggest that, for a qualitative description of the behavior of the electrical conductivity in relation to the number density of conductive wires, the mean-field approximation can be successfully applied when the number density of the fillers $n>2n_c$, where $n_c$ is the percolation threshold. However, note the mean-field approach slightly overestimates the electrical conductivity. We demonstrate that this overestimate is caused by the junction potential distribution.

Schreiber's transfer entropy is an important index for investigating the causal relationship between random variables. The Liang-Kleeman information flow is another analysis to demonstrate the causality within dynamical systems. Horowitz's information flow is introduced through multicomponent stochastic thermodynamics. In this study, I elucidate the relationship between Schreiber's transfer entropy and the Liang-Kleeman information flow through Horowitz's information flow. I consider the case in which the system changes according to the stochastic differential equation. I find that the Liang-Kleeman and Horowitz information flows differ by a term derived from the stochastic fluctuation. I also show that Schreiber's transfer entropy is not less than Horowitz's information flow. This study helps unify various indexes that determine the causal relationship between variables.

The 1D random Boltzmann-Lorentz equation has been connected with a set of stochastic hyperbolic equations. Therefore, the study of the Boltzmann-Lorentz gas with disordered scattering centers has been transformed into the analysis of a set of stochastic telegrapher's equations. For global binary disorder (Markovian and non-Markovian) exact analytical results for the second moment, the velocity autocorrelation function, and the self-diffusion coefficient are presented. We have demonstrated that time-fluctuations in the lost of energy in the telegrapher's equation, can delay the entrance to the diffusive regime, this issue has been characterized by a timescale $t_c$ which is a function of disorder parameters. Indeed, producing a longer ballistic dynamics in the transport process. In addition, fluctuations of the space probability distribution have been studied, showing that the mean value of a stochastic telegrapher's Fourier mode is a good statistical object to characterize the solution of the random Boltzmann-Lorentz gas. In a different context, the stochastic telegrapher's equation has also been related to the run-and-tumble model in Biophysics. Then a discussion devoted to the potential applications when swimmers' speed and tumbling rate have time fluctuations has been pointed out.

Resonant drift of nonlinear spiral waves occurs in various types of excitable media under periodic stimulation. Recently a novel methodology of optogenetics has emerged, which allows to affect excitability of cardiac cells and tissues by optical stimuli. In this paper we study if resonant drift of spiral waves in the heart can be induced by a homogeneous weak periodic optical stimulation of cardiac tissue. We use a two-variable and a detailed model of cardiac tissue and add description of depolarizing and hyperpolarizing optogenetic ionic currents. We show that weak periodic optical stimulation at a frequency equal to the natural rotation frequency of the spiral wave induces resonant drift for both depolarizing and hyperpolarizing optogenetic currents. We quantify these effects and study how the speed of the drift and its direction depend on the initial conditions. We also derive analytical formulas based on the response function theory which correctly predict the drift velocity and its trajectory. We conclude that optogenetic methodology can be used for control of spiral waves in cardiac tissue and discuss its possible applications.

How self-organization leads to the emergence of structure in social populations remains a fascinating and open question in the study of complex systems. One frequently observed structure that emerges again and again across systems is that of self-similar community, i.e., homophily. We use a game theoretic perspective to explore a case in which individuals choose affiliation partnerships based on only two factors: the value they place on having social contacts, and their risk tolerance for exposure to threat derived from social contact (e.g., infectious disease, threatening ideas, etc.). We show how diversity along just these two influences is sufficient to cause the emergence of self-organizing homophily in the population. We further consider a case in which extrinsic social factors influence the desire to maintain particular social ties, and show the robustness of emergent homophilic patterns to these additional influences. These results demonstrate how observable population-level homophily may arise out of individual behaviors that balance the value of social contacts against the potential risks associated with those contacts. We present and discuss these results in the context of outbreaks of infectious disease in human populations. Complementing the standard narrative about how social division alters epidemiological risk, we here show how epidemiological risk may deepen social divisions in human populations.

Probabilistic message-passing algorithms are developed for routing transmissions in multiwavelength optical communication networks, under node- and edge-disjoint routing constraints and for various objective functions. Global routing optimization is a hard computational task on its own but is made much more difficult under the node- and edge-disjoint constraints and in the presence of multiple wavelengths, a problem which dominates routing efficiency in real optical communication networks that carry most of the world's internet traffic. The scalable principled method we have developed is exact on trees but provides good approximate solutions on locally treelike graphs. It accommodates a variety of objective functions that correspond to low latency, load balancing, and consolidation of routes and can be easily extended to include heterogeneous signal-to-noise values on edges and a restriction on the available wavelengths per edge. It can be used for routing and managing transmissions on existing topologies as well as for designing and modifying optical communication networks. Additionally, it provides the tool for settling an open and much-debated question on the merit of wavelength-switching nodes and the added capabilities they provide. The methods have been tested on generated networks such as random-regular, Erdős Rényi, and power-law graphs, as well as on optical communication networks in the United Kingdom and United States. They show excellent performance with respect to existing methodology on small networks and have been scaled up to network sizes that are beyond the reach of most existing algorithms.

A simple coarse-grained model of DNA which includes both Watson-Crick and Hoogsteen base pairing has been used to study the melting and unzipping of triplex DNA. Using Langevin dynamics simulations, we reproduce the qualitative features of one-step and two-step thermal melting of triplex as seen in experiments. The thermal melting phase diagram shows the existence of a stable interchain three-strand complex (bubble-bound state). Our studies based on the mechanical unzipping of a triplex revealed that it is mechanically more stable compared to an isolated duplex-DNA.

Biological membrane is a complex self-assembly of lipids, sterols, and proteins organized as a fluid bilayer of two closely stacked lipid leaflets. Differential molecular interactions among its diverse constituents give rise to heterogeneities in the membrane lateral organization. Under certain conditions, heterogeneities in the two leaflets can be spatially synchronized and exist as registered domains across the bilayer. Several contrasting theories behind mechanisms that induce registration of nanoscale domains have been suggested. Following a recent study showing the effect of position of lipid tail unsaturation on domain registration behavior, we decided to develop an analytical theory to elucidate the driving forces that create and maintain domain registry across leaflets. Towards this, we formulated a Hamiltonian for a stacked lattice system where site variables capture the lipid molecular properties such as the position of unsaturation and various other interactions that could drive phase separation and interleaflet coupling. We solve the Hamiltonian using Monte Carlo simulations and create a complete phase diagram that reports the presence or absence of registered domains as a function of various Hamiltonian parameters. We find that the interleaflet coupling should be described as a competing enthalpic contribution due to interaction of lipid tail termini, primarily due to saturated-saturated interactions, and an interleaflet entropic contribution from overlap of unsaturated tail termini. A higher position of unsaturation is seen to provide weaker interleaflet coupling. Thermodynamically stable nanodomains could also be observed for certain points in the parameter space in our bilayer model, which were further verified by carrying out extended Monte Carlo simulations. These persistent noncoalescing registered nanodomains close to the lower end of the accepted nanodomain size range also point towards a possible “nanoscale” emulsion description of lateral heterogeneities in biological membrane leaflets.