scispace - formally typeset
Search or ask a question

Showing papers in "Physical Review E in 2022"


Journal ArticleDOI
TL;DR: In this article , the authors propose the notion of "energetic rigidity", in which all nontrivial deformations raise the energy of a structure, which is a more useful notion of rigidity in practice than two commonly used rigidity tests: first-order rigidity and second-order rigidness.
Abstract: Rigidity regulates the integrity and function of many physical and biological systems. This is the first of two papers on the origin of rigidity, wherein we propose that "energetic rigidity," in which all nontrivial deformations raise the energy of a structure, is a more useful notion of rigidity in practice than two more commonly used rigidity tests: Maxwell-Calladine constraint counting (first-order rigidity) and second-order rigidity. We find that constraint counting robustly predicts energetic rigidity only when the system has no states of self-stress. When the system has states of self-stress, we show that second-order rigidity can imply energetic rigidity in systems that are not considered rigid based on constraint counting, and is even more reliable than shear modulus. We also show that there may be systems for which neither first- nor second-order rigidity imply energetic rigidity. The formalism of energetic rigidity unifies our understanding of mechanical stability and also suggests new avenues for material design.

18 citations


Journal ArticleDOI
TL;DR: In this paper , the authors numerically test the complexity growth of unitarily evolving operators for a variety of exemplary systems, including one-dimensional and two-dimensional Ising models as well as one dimensional Heisenberg models.
Abstract: Recently, a hypothesis on the complexity growth of unitarily evolving operators was presented. This hypothesis states that in generic, nonintegrable many-body systems, the so-called Lanczos coefficients associated with an autocorrelation function grow asymptotically linear, with a logarithmic correction in one-dimensional systems. In contrast, the growth is expected to be slower in integrable or free models. In this paper, we numerically test this hypothesis for a variety of exemplary systems, including one-dimensional and two-dimensional Ising models as well as one-dimensional Heisenberg models. While we find the hypothesis to be practically fulfilled for all considered Ising models, the onset of the hypothesized universal behavior could not be observed in the attainable numerical data for the Heisenberg model. The proposed linear bound on operator growth associated with the hypothesis eventually stems from geometric arguments involving the locality of the Hamiltonian as well as the lattice configuration. We derive and investigate a related geometric bound, and we find that while the bound itself is not sharply achieved for any considered model, the hypothesis is nonetheless fulfilled in most cases.

15 citations


Journal ArticleDOI
TL;DR: This work shows that weighted simplicial complexes can represent higher-order networks without loss of information, allowing one at the same time to capture the weighted topology of the data.
Abstract: Hypergraphs and simplical complexes both capture the higher-order interactions of complex systems, ranging from higher-order collaboration networks to brain networks. One open problem in the field is what should drive the choice of the adopted mathematical framework to describe higher-order networks starting from data of higher-order interactions. Unweighted simplicial complexes typically involve a loss of information of the data, though having the benefit to capture the higher-order topology of the data. In this work we show that weighted simplicial complexes allow one to circumvent all the limitations of unweighted simplicial complexes to represent higher-order interactions. In particular, weighted simplicial complexes can represent higher-order networks without loss of information, allowing one at the same time to capture the weighted topology of the data. The higher-order topology is probed by studying the spectral properties of suitably defined weighted Hodge Laplacians displaying a normalized spectrum. The higher-order spectrum of (weighted) normalized Hodge Laplacians is studied combining cohomology theory with information theory. In the proposed framework we quantify and compare the information content of higher-order spectra of different dimension using higher-order spectral entropies and spectral relative entropies. The proposed methodology is tested on real higher-order collaboration networks and on the weighted version of the simplicial complex model "Network Geometry with Flavor."

13 citations


Journal ArticleDOI
Jie Lin, Kai Li, Jizhou He, Jie Ren, Jianhui Wang 
TL;DR: In this paper , the authors investigate heat and work in a Markovian state transition system with cyclic switching hot-cold temperatures, which operates as an Otto heat engine working in long but finite time.
Abstract: The Mpemba effect is a counterintuitive relaxation phenomenon whereby a system with a higher initial temperature may cool down to the thermal state faster than an identical system that was initially prepared at a lower temperature. Here, we investigate heat and work in a Markovian state transition system with cyclic switching hot-cold temperatures, which operates as an Otto heat engine working in long but finite time, either with or without the Mpemba effect. Under the condition of the periodic steady state having been reached, the time durations of the heating and cooling relaxation processes are determined by exploring a distance-from-equilibrium equivalent to the Kullback-Leibler divergence. We then numerically evaluate and compare the averages and variances of both the work and the power output of two scenarios with and without the Mpemba effect. The results show that the Markovian Mpemba effect can enhance the machine performance by significantly increasing the power output for a given efficiency without sacrificing the stability.

13 citations


Journal ArticleDOI
TL;DR: In this paper , the authors apply their formalism to examples in two dimensions: underconstrained random regular spring networks, vertex models, and jammed packings of soft particles, and show that constraint counting is equivalent to energetic rigidity as long as prestresses in the system are sufficiently small.
Abstract: This is the second paper devoted to energetic rigidity, in which we apply our formalism to examples in two dimensions: Underconstrained random regular spring networks, vertex models, and jammed packings of soft particles. Spring networks and vertex models are both highly underconstrained, and first-order constraint counting does not predict their rigidity, but second-order rigidity does. In contrast, spherical jammed packings are overconstrained and thus first-order rigid, meaning that constraint counting is equivalent to energetic rigidity as long as prestresses in the system are sufficiently small. Aspherical jammed packings on the other hand have been shown to be jammed at hypostaticity, which we use to argue for a modified constraint counting for systems that are energetically rigid at quartic order.

12 citations


Journal ArticleDOI
TL;DR: In this paper , a thermal gradient-induced circular motion of particles placed on ferroelectric nematic liquid crystal sessile drops is demonstrated and explained, where the texture without tracer particles appears completely steady indicating laminar flow.
Abstract: 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.

12 citations


Journal ArticleDOI
TL;DR: In this paper , a collection of methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models is presented. But the authors focus on the Holstein model on the square lattice.
Abstract: We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and we report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but it incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

11 citations


Journal ArticleDOI
TL;DR: In this paper , a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions, is presented, which gives a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force.
Abstract: When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force, which plays an essential role in determining both the equilibrium and dynamic properties of interacting fluids. We present a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions. By consideration of local Noether-invariance against spatial distortion of the system we demonstrate the fundamental status of the Yvon-Born-Green (YBG) equation as a local force-balance within the fluid. Using the inhomogeneous Ornstein-Zernike equation we show that the two-body correlations are density functionals and, thus, that the average interparticle force entering the YBG equation is also a functional of the one-body density. The force-based theory we develop provides an alternative to standard density functional theory for the study of inhomogeneous systems both in- and out-of-equilibrium. We compare force-based density profiles to the results of the standard potential-based (dynamical) density functional theory. In-equilibrium, we confirm both analytically and numerically that the standard approach yields profiles that are consistent with the compressibility pressure, whereas the force-density functional gives profiles consistent with the virial pressure. For both approaches we explicitly prove the hard-wall contact theorem that connects the value of the density profile at the hard-wall with the bulk pressure. The structure of the theory offers deep insights into the nature of correlation in dense and inhomogeneous systems.

11 citations


Journal ArticleDOI
TL;DR: In this article , the performance of three-level quantum absorption refrigerators, paradigmatic autonomous quantum thermal machines, and reveal central impacts of strong couplings between the working system and the thermal baths.
Abstract: We study the performance of three-level quantum absorption refrigerators, paradigmatic autonomous quantum thermal machines, and reveal central impacts of strong couplings between the working system and the thermal baths. Using the reaction coordinate quantum master equation method, which treats system-bath interactions beyond weak coupling, we demonstrate that in a broad range of parameters the cooling window at strong coupling can be captured by a weak-coupling theory, albeit with parameters renormalized by the system-bath coupling energy. As a result, at strong system-bath couplings the window of cooling is significantly reshaped compared to predictions of weak-coupling treatments. We further show that strong coupling admits direct transport pathways between the thermal reservoirs. Such beyond-second-order transport mechanisms are typically detrimental to the performance of quantum thermal machines. Our study reveals that it is inadequate to claim for either a suppression or an enhancement of the cooling performance as one increases system-bath coupling-when analyzed against a single parameter and in a limited domain. Rather, a comprehensive approach should be adopted so as to uncover the reshaping of the operational window.

11 citations


Journal ArticleDOI
TL;DR: In this paper , an overview of the ferroelectric nematic liquid crystalline phases is presented, linking history and theoretical predictions to a general outlook of the development and properties of the materials exhibiting Ferroelectric nano-nematic phases.
Abstract: Recent experimental realization of ferroelectric nematic liquid crystalline phases stimulated material development and numerous experimental studies of these phases, guided by their fundamental and applicative interest. In this Perspective, we give an overview of this emerging field by linking history and theoretical predictions to a general outlook of the development and properties of the materials exhibiting ferroelectric nematic phases. We will highlight the most relevant observations to date, e.g., giant dielectric permittivity values, polarization values an order of magnitude larger than in classical ferroelectric liquid crystals, and nonlinear optical coefficients comparable with several ferroelectric solid materials. Key observations of anchoring and electro-optic behavior will also be examined. The collected contributions lead to a final discussion on open challenges in materials development, theoretical description, experimental explorations, and possible applications of the ferroelectric phases.

11 citations


Journal ArticleDOI
TL;DR: In this paper , the authors combine ab initio path integral Monte Carlo (PIMC) simulations with fixed ion configurations from density functional theory molecular dynamics (DFT-MD) simulations to solve the electronic problem for hydrogen under warm dense matter conditions.
Abstract: We combine ab initio path integral Monte Carlo (PIMC) simulations with fixed ion configurations from density functional theory molecular dynamics (DFT-MD) simulations to solve the electronic problem for hydrogen under warm dense matter conditions [Böhme et al., Phys. Rev. Lett. 129, 066402 (2022)0031-900710.1103/PhysRevLett.129.066402]. The problem of path collapse due to the Coulomb attraction is avoided by utilizing the pair approximation, which is compared against the simpler Kelbg pair potential. We find very favorable convergence behavior towards the former. Since we do not impose any nodal restrictions, our PIMC simulations are afflicted with the notorious fermion sign problem, which we analyze in detail. While computationally demanding, our results constitute an exact benchmark for other methods and approximations within DFT. Our setup gives us the unique capability to study important properties of warm dense hydrogen such as the electronic static density response and exchange-correlation kernel without any model assumptions, which will be very valuable for a variety of applications such as the interpretation of experiments and the development of new XC functionals.

Journal ArticleDOI
TL;DR: In this paper, the authors apply Selberg's trace formula to solve problems in hyperbolic band theory, a recently developed extension of Bloch theory to model band structures on experimentally realized Hyperbolic lattices.
Abstract: We apply Selberg's trace formula to solve problems in hyperbolic band theory, a recently developed extension of Bloch theory to model band structures on experimentally realized hyperbolic lattices. For this purpose we incorporate the higher-dimensional crystal momentum into the trace formula and evaluate the summation for periodic orbits on the Bolza surface of genus two. We apply the technique to compute partition functions on the Bolza surface and propose an approximate relation between the lowest bands on the Bolza surface and on the {8,3} hyperbolic lattice. We discuss the role of automorphism symmetry and its manifestation in the trace formula.

Journal ArticleDOI
TL;DR: In this article , the average angle an oxygen atom makes with its neighbors to describe the structural environment of a water molecule was used to estimate the fraction of the two liquid forms in water.
Abstract: The existence of two structural forms in liquid water has been a point of discussion for a long time. A phase transition between these two forms of liquid water has been proposed based on evidence from molecular simulations, and experiments have also been very recently able to track the proposed transition of the low-density liquid form to the high-density liquid form. We propose to use the average angle an oxygen atom makes with its neighbors to describe the structural environment of a water molecule. The distribution of this order parameter is observed to have two peaks with one peak at ∼109.5^{∘}, corresponding to the internal angle of a regular tetrahedron, indicating tetrahedral arrangement. The other peak corresponds to an environment with a tighter arrangement of neighboring molecules. The distribution of O-O-O angles is decomposed into two skewed distributions to estimate the fractions of the two liquid forms in water. A good similarity is observed between the temperature and pressure trends of fractions of locally favored tetrahedral structure (LFTS) form estimated using the new order parameter and the reports in the literature, over a range of temperatures and pressures. We also compare the structural environments indicated by different order parameters and find that the order parameter proposed in this paper captures the structure of first solvation shell of the LFTS accurately.

Journal ArticleDOI
TL;DR: In this paper , the authors investigated the transition to synchronization of oscillator ensembles encoded by simplicial complexes in which pairwise and higher-order coupling weights alter with time through a rate-based adaptive mechanism inspired by the Hebbian learning rule.
Abstract: This Letter investigates the transition to synchronization of oscillator ensembles encoded by simplicial complexes in which pairwise and higher-order coupling weights alter with time through a rate-based adaptive mechanism inspired by the Hebbian learning rule. These simultaneously evolving disparate adaptive coupling weights lead to a phenomenon in that the in-phase synchronization is completely obliterated; instead, the antiphase synchronization is originated. In addition, the onsets of antiphase synchronization and desynchronization are manageable through both dyadic and triadic learning rates. The theoretical validation of these numerical assessments is delineated thoroughly by employing Ott-Antonsen dimensionality reduction. The framework and results of the Letter would help understand the underlying synchronization behavior of a range of real-world systems, such as the brain functions and social systems where interactions evolve with time.

Journal ArticleDOI
TL;DR: In this paper , the authors report a unique phenomenon in which liquid metal droplets (LMDs) under a pure ac electric field pump fluid can be switched by simply shifting the position of the LMD within the cylindrical chamber.
Abstract: We report a unique phenomenon in which liquid metal droplets (LMDs) under a pure ac electric field pump fluid. Unlike the directional pumping that occurs upon reversing the electric field polarity under a dc signal, this phenomenon allows the direction of fluid motion to be switched by simply shifting the position of the LMD within the cylindrical chamber. The physical mechanism behind this phenomenon has been termed Marangoni flow, caused by nonlinear electrocapillary stress. Under the influence of a localized, asymmetric ac electric field, the polarizable surface of the position-offset LMD produces a net time-averaged interfacial tension gradient that scales with twice the field strength, and thus pumps fluid unidirectionally. However, the traditional linear RC circuit polarization model of the LMD/electrolyte interface fails to capture the correct pump-flow direction when the thickness of the LMD oxide skin is non-negligible compared to the Debye length. Therefore, we developed a physical description by treating the oxide layer as a distributed capacitance with variable thickness and connected with the electric double layer. The flow profile is visualized via microparticle imaging velocimetry, and excellent consistency is found with simulation results obtained from the proposed nonlinear model. Furthermore, we investigate the effects of relevant parameters on fluid pumping and discuss a special phenomenon that does not exist in dc control systems. To our knowledge, no previous work addresses LMDs in this manner and uses a zero-mean ac electric field to achieve stable, adjustable directional pumping of a low-conductivity solution.

Journal ArticleDOI
TL;DR: In this paper , a microscopic theory of the CA engine realized with a highly underdamped Brownian particle in a class of nonharmonic potentials was proposed. And the exact expression of the protocol associated with the maximum power for any given efficiency, and analytical results of the power and efficiency statistics for the Brownian CA engine were derived.
Abstract: The Curzon-Ahlborn (CA) efficiency, as the efficiency at the maximum power (EMP) of the endoreversible Carnot engine, has significant impact on finite-time thermodynamics. However, the CA engine is based on many assumptions. In the past few decades, although a lot of efforts have been made, a microscopic theory of the CA engine is still lacking. By adopting the method of the stochastic differential equation of energy, we formulate a microscopic theory of the CA engine realized with a highly underdamped Brownian particle in a class of nonharmonic potentials. This theory gives microscopic interpretation of all assumptions made by Curzon and Ahlborn. In other words, we find a microscopic counterpart of the CA engine in stochastic thermodynamics. Also, based on this theory, we derive the explicit expression of the protocol associated with the maximum power for any given efficiency, and we obtain analytical results of the power and the efficiency statistics for the Brownian CA engine. Our research brings new perspectives to experimental studies of finite-time microscopic heat engines featured with fluctuations.

Journal ArticleDOI
TL;DR: In this article , the performance of slowly driven meso-and microscale refrigerators and heat engines that operate between two thermal baths with a small temperature difference was analyzed and it was shown that such devices can work arbitrarily close to their Carnot limit only if heat leaks between the baths are fully suppressed.
Abstract: We analyze the performance of slowly driven meso- and microscale refrigerators and heat engines that operate between two thermal baths with a small temperature difference. Using a general scaling argument, we show that such devices can work arbitrarily close to their Carnot limit only if heat leaks between the baths are fully suppressed. Their power output is then subject to a universal geometric bound that decays quadratically to zero at the Carnot limit. This bound can be asymptotically saturated in the quasistatic limit if the driving protocols are suitably optimized and the temperature difference between the baths goes to zero with the driving frequency. These results hold under generic conditions for any thermodynamically consistent dynamics admitting a well-defined adiabatic-response regime and a generalized Onsager symmetry. For illustration, we work out models of a qubit-refrigerator and a coherent charge pump operating as a cooling device.

Journal ArticleDOI
TL;DR: In this article , a generalized (2+1)-dimensional short-wave asymptotic model is derived for the nonlinear propagation of a magnetic soliton in a ferromagnetic film.
Abstract: In this paper, we study in detail the nonlinear propagation of a magnetic soliton in a ferromagnetic film. The sample is magnetized to saturation by an external field perpendicular to film plane. A generalized (2+1)-dimensional short-wave asymptotic model is derived. The bilinearlike forms of this equation are constructed and exact magnetic line soliton solutions are exhibited. It is observed that a series of stable lumps can be generated by an unstable magnetic soliton under Gaussian disturbance. Such magnetic lumps are highly stable and can maintain their shapes and velocities during evolution or collision. The interaction between lump and magnetic solitons, as well as the interaction between two lumps, are numerically investigated. We further discuss the nonlinear motion of lumps in ferrites with Gilbert damping and inhomogeneous exchange effects. The results show that the Gilbert-damping effects make the amplitude and velocity of the magnetic lump decay exponentially during propagation. And the shock waves are generated from a lump when quenching the strength of inhomogeneous exchange.

Journal ArticleDOI
TL;DR: In this paper, the authors considered a network of non-identical and globally coupled Stuart-Landau oscillators with additional perturbation and found that the observed dynamical states do not preserve the parity (P) symmetry in the absence of perturbations.
Abstract: Frequency plays a crucial role in exhibiting various collective dynamics in the coexisting corotating and counter-rotating systems. To illustrate the impact of counter-rotating frequencies, we consider a network of nonidentical and globally coupled Stuart-Landau oscillators with additional perturbation. Primarily, we investigate the dynamical transitions in the absence of perturbation, demonstrating that the transition from desynchronized state to cluster oscillatory state occurs through an interesting partial synchronization state in the oscillatory regime. Following this, the system dynamics transits to amplitude death and oscillation death states. Importantly, we find that the observed dynamical states do not preserve the parity (P) symmetry in the absence of perturbation. When the perturbation is increased one can note that the system dynamics exhibits a kind of transition which corresponds to a change from incoherent mixed synchronization to coherent mixed synchronization through a chimera state. In particular, incoherent mixed synchronization and coherent mixed synchronization states completely preserve the P symmetry, whereas the chimera state preserves the P symmetry only partially. To demonstrate the occurrence of such partial symmetry-breaking (chimera) state, we use basin stability analysis and discover that partial symmetry breaking exists as a result of the coexistence of symmetry-preserving and symmetry-breaking behavior in the initial state space. Further, a measure of the strength of P symmetry is established to quantify the P symmetry in the observed dynamical states. Subsequently, the dynamical transitions are investigated in the parametric spaces. Finally, by increasing the network size, the robustness of the chimera state is also inspected, and we find that the chimera state is robust even in networks of larger sizes. We also show the generality of the above results in the related reduced phase. model as well as in other coupled models such as the globally coupled van der Pol and Rössler oscillators.

Journal ArticleDOI
TL;DR: In this article, the underdamped motion of a passive particle in an active environment is studied using phase space path integral method and the probability distribution function of position and velocity for a free and a harmonically bound particle.
Abstract: We study the underdamped motion of a passive particle in an active environment. Using the phase space path integral method we find the probability distribution function of position and velocity for a free and a harmonically bound particle. The environment is characterized by an active noise which is described as the Ornstein-Uhlenbeck process (OUP). Taking two similar, yet slightly different OUP models, it is shown how inertia along with other relevant parameters affect the dynamics of the particle. Further we investigate the work fluctuations of a harmonically trapped particle by considering the trap center being pulled at a constant speed. Finally, the fluctuation theorem of work is validated with an effective temperature in the steady-state limit.

Journal ArticleDOI
TL;DR: In this article , the authors study a few models, systematically progressing in their complexity, to understand the usefulness of resetting and 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.
Abstract: 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.

Journal ArticleDOI
TL;DR: In this paper , the work-minimizing protocol problem in the standard form of an optimal control theory problem is shown to be equivalent to solving a system of Hamiltonian partial differential equations, which admit efficiently calculable numerical solutions.
Abstract: Recent studies have explored finite-time dissipation-minimizing protocols for stochastic thermodynamic systems driven arbitrarily far from equilibrium, when granted full external control to drive the system. However, in both simulation and experimental contexts, systems often may only be controlled with a limited set of degrees of freedom. Here, going beyond slow- and fast-driving approximations employed in previous studies, we obtain exact finite-time optimal protocols for this limited-control setting. By working with deterministic Fokker-Planck probability density time evolution, we can frame the work-minimizing protocol problem in the standard form of an optimal control theory problem. We demonstrate that finding the exact optimal protocol is equivalent to solving a system of Hamiltonian partial differential equations, which in many cases admit efficiently calculable numerical solutions. Within this framework, we reproduce analytical results for the optimal control of harmonic potentials and numerically devise optimal protocols for two anharmonic examples: varying the stiffness of a quartic potential and linearly biasing a double-well potential. We confirm that these optimal protocols outperform other protocols produced through previous methods, in some cases by a substantial amount. We find that for the linearly biased double-well problem, the mean position under the optimal protocol travels at a near-constant velocity. Surprisingly, for a certain timescale and barrier height regime, the optimal protocol is also nonmonotonic in time.

Journal ArticleDOI
TL;DR: It is found that group interactions amplify small initial opinion biases, accelerate the formation of consensus and lead to a drift of the average magnetization.
Abstract: We study the effect of group interactions on the emergence of consensus in a spin system. Agents with discrete opinions {0,1} form groups. They can change their opinion based on their group's influence (voter dynamics), but groups can also split and merge (adaptation). In a hypergraph, these groups are represented by hyperedges of different sizes. The heterogeneity of group sizes is controlled by a parameter β. To study the impact of β on reaching consensus, we provide extensive computer simulations and compare them with an analytic approach for the dynamics of the average magnetization. We find that group interactions amplify small initial opinion biases, accelerate the formation of consensus, and lead to a drift of the average magnetization. The conservation of the initial magnetization, known for basic voter models, is no longer obtained.

Journal ArticleDOI
TL;DR: This paper analyzes the three-dimensional narrow capture problem in the more realistic case of partially reactive target boundaries and derives an asymptotic expansion of the joint probability density for particle position and the so-called boundary local time, which characterizes the amount of time that a Brownian particle spends on a totally reflecting boundary.
Abstract: A general topic of current interest is the analysis of diffusion problems in singularly perturbed domains with small interior targets or traps (the narrow capture problem). One major application is to intracellular diffusion, where the targets typically represent some form of reactive biochemical substrate. Most studies of the narrow capture problem treat the target boundaries as totally absorbing (Dirichlet), that is, the chemical reaction occurs immediately on first encounter between particle and target surface. In this paper, we analyze the three-dimensional narrow capture problem in the more realistic case of partially reactive target boundaries. We begin by considering classical Robin boundary conditions. Matching inner and outer solutions of the single-particle probability density, we derive an asymptotic expansion of the Laplace transformed flux into each reactive surface in powers of ε, where ερ is a given target size. In turn, the fluxes determine the splitting probabilities for target absorption. We then extend our analysis to more general types of reactive targets by combining matched asymptotic analysis with an encounter-based formulation of diffusion-mediated surface reactions. That is, we derive an asymptotic expansion of the joint probability density for particle position and the so-called boundary local time, which characterizes the amount of time that a Brownian particle spends in the neighborhood of a point on a totally reflecting boundary. The effects of surface reactions are then incorporated via an appropriate stopping condition for the boundary local time. Robin boundary conditions are recovered in the special case of an exponential law for the stopping local times. Finally, we illustrate the theory by exploring how the leading-order contributions to the splitting probabilities depend on the choice of surface reactions. In particular, we show that there is an effective renormalization of the target radius of the form ρ→ρ-Ψ[over ̃](1/ρ), where Ψ[over ̃] is the Laplace transform of the stopping local time distribution.

Journal ArticleDOI
TL;DR: In this paper , the authors analyze and compare two strategies currently used to obtain (approximately) steady states in strong-coupling regime, based on perturbative expansions and reaction coordinate mapping.
Abstract: Motivated by the growing importance of strong system-bath coupling in several branches of quantum information and related technological applications, we analyze and compare two strategies currently used to obtain (approximately) steady states in strong-coupling regime. The first strategy is based on perturbative expansions while the second one uses reaction coordinate mapping. Focusing on the widely used spin-boson model, we show that the predictions of these two strategies coincide in many situations. This confirms and strengthens the relevance of both techniques. Beyond that, it is also crucial to know precisely their respective range of validity. In that perspective, thanks to their different limitations, we use one to benchmark the other. We introduce and successfully test some very simple validity criteria for both strategies, bringing some answers to the question of the validity range.

Journal ArticleDOI
TL;DR: In this article , a generic overdamped Langevin dynamics driven out of equilibrium by both time-dependent and non-conservative forces is decomposed into two positive terms, termed excess and housekeeping entropy.
Abstract: For a generic overdamped Langevin dynamics driven out of equilibrium by both time-dependent and nonconservative forces, the entropy production rate can be decomposed into two positive terms, termed excess and housekeeping entropy. However, this decomposition is not unique: There are two distinct decompositions, one due to Hatano and Sasa, the other one due to Maes and Netočný. Here we establish the connection between these two decompositions and provide a simple, geometric interpretation. We show that this leads to a decomposition of the entropy production rate into three positive terms, which we call the excess, housekeeping, and coupling part, respectively. The coupling part characterizes the interplay between the time-dependent and nonconservative forces. We also derive thermodynamic uncertainty relations for the excess and housekeeping entropy in both the Hatano-Sasa and Maes-Netočný decomposition and show that all quantities obey integral fluctuation theorems. We illustrate the decomposition into three terms using a solvable example of a dragged particle in a nonconservative force field.

Journal ArticleDOI
TL;DR: Comparing deep random feature models to deep networks in which all layers are trained provides a detailed characterization of the interplay between width, depth, data density, and prior mismatch and begins to elucidate how architectural details affect generalization performance in this simple class of deep regression models.
Abstract: Understanding how feature learning affects generalization is among the foremost goals of modern deep learning theory. Here, we study how the ability to learn representations affects the generalization performance of a simple class of models: deep Bayesian linear neural networks trained on unstructured Gaussian data. By comparing deep random feature models to deep networks in which all layers are trained, we provide a detailed characterization of the interplay between width, depth, data density, and prior mismatch. We show that both models display samplewise double-descent behavior in the presence of label noise. Random feature models can also display modelwise double descent if there are narrow bottleneck layers, while deep networks do not show these divergences. Random feature models can have particular widths that are optimal for generalization at a given data density, while making neural networks as wide or as narrow as possible is always optimal. Moreover, we show that the leading-order correction to the kernel-limit learning curve cannot distinguish between random feature models and deep networks in which all layers are trained. Taken together, our findings begin to elucidate how architectural details affect generalization performance in this simple class of deep regression models.

Journal ArticleDOI
TL;DR: In this paper, the authors present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number, which imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system.
Abstract: We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. We perform a variety of tests, and in all cases find rapid convergence of the chemical potential-inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

Journal ArticleDOI
TL;DR: In this paper , the authors apply periodic control to realize a quantum thermal transistor, which they termed as the Floquet Quantum thermal Transistor, and achieve large amplification factors even for fixed bath temperatures.
Abstract: We apply periodic control to realize a quantum thermal transistor, which we term as the Floquet Quantum thermal Transistor. Periodic modulation allows us to control the heat flows and achieve large amplification factors even for fixed bath temperatures. Importantly, this transistor effect persists in the cut-off region, where traditional quantum thermal transistors operating in absence of periodic modulation, fail to act as viable heat modulation devices.

Journal ArticleDOI
TL;DR: In this article , the authors studied the dynamics of an interesting class of vector solitons in the long-wave-short-wave resonance interaction (LSRI) system and derived the fundamental and multisoliton solutions in a compact way using Gram-determinants.
Abstract: In this paper, we study the dynamics of an interesting class of vector solitons in the long-wave-short-wave resonance interaction (LSRI) system. The model that we consider here describes the nonlinear interaction of long wave and two short waves and it generically appears in several physical settings. To derive this class of nondegenerate vector soliton solutions we adopt the Hirota bilinear method with the more general form of admissible seed solutions with nonidentical distinct propagation constants. We express the resultant fundamental as well as multisoliton solutions in a compact way using Gram-determinants. The general fundamental vector soliton solution possesses several interesting properties. For instance, the double-hump or a single-hump profile structure including a special flattop profile form results in when the soliton propagates in all the components with identical velocities. Interestingly, in the case of nonidentical velocities, the soliton number is increased to two in the long-wave component, while a single-humped soliton propagates in the two short-wave components. We establish through a detailed analysis that the nondegenerate multisolitons in contrast to the already known vector solitons (with identical wave numbers) can undergo three types of elastic collision scenarios: (i) shape-preserving, (ii) shape-altering, and (iii) a shape-changing collision, depending on the choice of the soliton parameters. Here, by shape-altering we mean that the structure of the nondegenerate soliton gets modified slightly during the collision process, whereas if the changes occur appreciably then we call such a collision as shape-changing collision. We distinguish each of the collision scenarios, by deriving a zero phase shift criterion with the help of phase constants. Very importantly, the shape-changing behavior of the nondegenerate vector solitons is observed in the long-wave mode also, along with corresponding changes in the short-wave modes, and this nonlinear phenomenon has not been observed in the already known vector solitons. In addition, we point out the coexistence of nondegenerate and degenerate solitons simultaneously along with the associated physical consequences. We also indicate the physical realizations of these general vector solitons in nonlinear optics, hydrodynamics, and Bose-Einstein condensates. Our results are generic and they will be useful in these physical systems and other closely related systems including plasma physics when the long-wave-short-wave resonance interaction is taken into account.