Linear Response Formula Research Articles

On convergence of linear response formulae in some piecewise hyperbolic maps

Abstract When high-dimensional non-uniformly hyperbolic chaotic systems undergo dynamical perturbations, their long-time statistics are generally observed to respond differentiably with respect to the perturbation. Although important in applications, this differentiability, which is thought to be connected to the dimensionality of the system, has remained resistant to rigorous study outside of the one-dimensional setting. To model non-uniformly hyperbolic systems in multiple dimensions, we consider a family of the mathematically tractable class of piecewise smooth hyperbolic maps, the Lozi maps. For these maps, we prove that the existence of a formal derivative of the response reduces to conditional mixing of the SRB measure on the singularity set. Heuristically, this suggests that Lozi maps, and piecewise uniformly expanding maps more generally, should have linear response.

Regularity and linear response formula of the SRB measures for solenoidal attractors

Abstract We show that a class of higher-dimensional hyperbolic endomorphisms admit absolutely continuous invariant probabilities whose densities are regular and vary differentiably with respect to the dynamical system. The maps we consider are skew-products given by $T(x,y) = (E (x), C(x,y))$ , where E is an expanding map of $\mathbb {T}^u$ and C is a contracting map on each fiber. If $\inf |\!\det DT| \inf \| (D_yC)^{-1}\| ^{-2s}>1$ for some ${s<r-(({u+d})/{2}+1)}$ , $r \geq 2$ , and T satisfies a transversality condition between overlaps of iterates of T (a condition which we prove to be $C^r$ -generic under mild assumptions), then the SRB measure $\mu _T$ of T is absolutely continuous and its density $h_T$ belongs to the Sobolev space $H^s({\mathbb {T}}^u\times {\mathbb {R}}^d)$ . When $s> {u}/{2}$ , it is also valid that the density $h_T$ is differentiable with respect to T. Similar results are proved for thermodynamical quantities for potentials close to the geometric potential.

Classical gravitational observables from the Eikonal operator

We propose two possible eikonal operators encoding the effects of classical radiation as coherent states of gravitons and show how to compute from them different classical observables. In the first proposal, only genuinely propagating gravitons are included, while, in the second, zero-frequency modes are added in order to recover the effects of a static gravitational field. We first calculate the radiated energy-momentum and the change in each particle's momentum, or impulse, to 3PM order finding agreement with the literature. We then calculate the angular momentum of the gravitational field after the collision. In order to do so, we adapt the method of reverse unitarity to the presence of derivatives in the operators describing the angular momentum and reproduce the result of [1] obtained by resumming the small-velocity expansion. As a new application, we derive also the variation in each particle's angular momentum up to 3PM: calculating separately field and particle contributions allows us to check the balance laws explicitly. We also show how the eikonal operator encodes the linear-response formula of Bini-Damour by deriving the linear radiation-reaction contribution to the transverse impulse at 4PM.

Dynamical fluctuations of a tracer coupled to active and passive particles

We study the induced dynamics of an inertial tracer particle elastically coupled to passive or active Brownian particles. We integrate out the environment degrees of freedom to obtain the exact effective equation of the tracer—a generalized Langevin equation in both cases. In particular, we find the exact form of the dissipation kernel and effective noise experienced by the tracer and compare it with the phenomenological modeling of active baths used in previous studies. We show that the second fluctuation-dissipation relation (FDR) does not hold at early times for both cases. However, at finite times, the tracer dynamics violate (obeys) the FDR for the active (passive) environment. We calculate the linear response formulas in this regime for both cases and show that the passive medium satisfies an equilibrium fluctuation response relation, while the active medium does not—we quantify the extent of this violation explicitly. We show that though the active medium generally renders a nonequilibrium description of the tracer, an effective equilibrium picture emerges asymptotically in the small activity limit of the medium. We also calculate the mean squared velocity and mean squared displacement of the tracer and report how they vary with time.

Linear response, Hamiltonian, and radiative spinning two-body dynamics

Using the spinning, supersymmetric worldline quantum field theory formalism we compute the momentum impulse and spin kick from a scattering of two spinning black holes or neutron stars up to quadratic order in spin at third post-Minkowskian (PM) order, including radiation-reaction effects and with arbitrarily misaligned spin directions. Parts of these observables, both conservative and radiative, are also inferred from lower-PM scattering data by extending Bini and Damour's linear response formula to include misaligned spins. By solving Hamilton's equations of motion we also use a conservative scattering angle to infer a complete 3PM two-body Hamiltonian including finite-size corrections and misaligned spin-spin interactions. Finally, we describe mappings to the bound two-body dynamics for aligned spin vectors: including a numerical plot of the binding energy for circular orbits compared with numerical relativity, analytic confirmation of the NNLO PN binding energy, and the energy loss over successive orbits.

Stress Tensor of Single Rigid Dumbbell by Virtual Work Method

We derive the stress tensor of a rigid dumbbell by using the virtual work method. In the virtual work method, we virtually apply a small deformation to the system, and relate the change of the energy to the work done by the stress tensor. A rigid dumbbell consists of two particles connected by a rigid bond of which length is constant (the rigid constraint). The energy of the rigid dumbbell consists only on the kinetic energy. Also, only the deformations which do not violate the rigid constraint are allowed. Thus we need the dynamic equations which are consistent with the rigid constraint to apply the virtual deformation. We rewrite the dynamic equations for the underdamped SLLOD-type dynamic equations into the forms which are consistent with the rigid constraint. Then we apply the virtual deformation to a rigid dumbbell based on the obtained dynamic equations. We derive the stress tensor for the rigid dumbbell model from the change of the kinetic energy. Finally, we take the overdamped limit and derive the stress tensor and the dynamic equation for the overdamped rigid dumbbell. We show that the Green-Kubo type linear response formula can be reproduced by combining the stress tensor and the dynamic equation at the overdamped limit.

Angular momentum loss in gravitational scattering, radiation reaction, and the Bondi gauge ambiguity

Recently, Damour computed the radiation reaction on gravitational scattering as the (linear) response to the angular momentum loss which he found to be of O(G2) in the gravitational constant. This is a puzzle because any amplitude calculation would predict both radiated energy and radiated angular momentum to start only at O(G3). Another puzzle is that the resultant radiation reaction, of O(G3), is nevertheless correct and confirmed by a number of direct calculations. We ascribe these puzzles to the BMS ambiguity in defining angular momentum. The loss of angular momentum is to be counted out from the ADM value and, therefore, should be calculated in the so-called canonical gauge under the BMS transformations in which the remote-past limit of the Bondi angular momentum coincides with the ADM angular momentum. This calculation correctly gives the O(G3) radiative loss. On the other hand, we introduce a gauge in which the Bondi light cones tend asymptotically to those emanating from the center of mass world line. We find that the angular momentum loss in this gauge is precisely the one used by Damour for his radiation reaction result. We call this new gauge “intrinsic” and argue that, although the radiated angular momentum is to be computed in the canonical gauge, any mechanical calculation of gauge-dependent quantities – such as angular momentum – gives the result in the intrinsic gauge. Therefore, it is this gauge that should be used in the linear response formula. This solves the puzzles and establishes the correspondence between the intrinsic mechanical calculations and the Bondi formalism.

Markovian repeated interaction quantum systems

We study a class of dynamical semigroups [Formula: see text] that emerge, by a Feynman–Kac type formalism, from a random quantum dynamical system [Formula: see text] driven by a Markov chain [Formula: see text]. We show that the almost sure large time behavior of the system can be extracted from the large [Formula: see text] asymptotics of the semigroup, which is in turn directly related to the spectral properties of the generator [Formula: see text]. As a physical application, we consider the case where the [Formula: see text]’s are the reduced dynamical maps describing the repeated interactions of a system [Formula: see text] with thermal probes [Formula: see text]. We study the full statistics of the entropy in this system and derive a fluctuation theorem for the heat exchanges and the associated linear response formulas.

Efficient Computation of Linear Response of Chaotic Attractors with One-Dimensional Unstable Manifolds

This paper presents the space-split sensitivity or the S3 algorithm to transform Ruelle's linear response formula into a well-conditioned ergodic-averaging computation. We prove a decomposition of Ruelle's formula that is differentiable on the unstable manifold, which we assume to be one-dimensional. This decomposition of Ruelle's formula ensures that one of the resulting terms, the stable contribution, can be computed using a regularized tangent equation, similar to in a non-chaotic system. The remaining term, known as the unstable contribution, is regularized and converted into an efficiently computable ergodic average. In this process, we develop new algorithms, which may be useful beyond linear response, to compute i) a fundamental statistical quantity we introduce called the density gradient, and ii) the unstable derivatives of the regularized tangent vector field and the unstable direction. We prove that the S3 algorithm, which combines these computational ingredients that enter the stable and unstable contribution, converges like a Monte Carlo approximation of Ruelle's formula. The algorithm presented here is hence a first step toward full-fledged applications of sensitivity analysis in chaotic systems, wherever such applications have been limited due to lack of availability of long-term sensitivities.

Optimal transportation and stationary measures for iterated function systems

AbstractIn this paper we show how ideas, methods and results from optimal transportation can be used to study various aspects of the stationary measures of Iterated Function Systems equipped with a probability distribution. We recover a classical existence and uniqueness result under a contraction-on-average assumption, prove generalised moment bounds from which tail estimates can be deduced, consider the convergence of the empirical measure of an associated Markov chain, and prove in many cases the Lipschitz continuity of the stationary measure when the system is perturbed, with as a consequence a “linear response formula” at almost every parameter of the perturbation.

An ergodic-averaging method to differentiate covariant Lyapunov vectors

Covariant Lyapunov vectors or CLVs span the expanding and contracting directions of perturbations along trajectories in a chaotic dynamical system. Due to efficient algorithms to compute them that only utilize trajectory information, they have been widely applied across scientific disciplines, principally for sensitivity analysis and predictions under uncertainty. In this paper, we develop a numerical method to compute the directional derivatives of CLVs along their own directions. Similar to the computation of CLVs, the present method for their derivatives is iterative and analogously uses the second-order derivative of the chaotic map along trajectories, in addition to the Jacobian. We validate the new method on a super-contracting Smale-Williams Solenoid attractor. We also demonstrate the algorithm on several other examples including smoothly perturbed Arnold Cat maps, and the Lorenz attractor, obtaining visualizations of the curvature of each attractor. Furthermore, we reveal a fundamental connection of the CLV self-derivatives with a statistical linear response formula.

Regularity of Characteristic Exponents and Linear Response for Transfer Operator Cocycles

We consider cocycles obtained by composing sequences of transfer operators with positive weights, associated with uniformly expanding maps (possibly having countably many branches) and depending upon parameters. Assuming $C^k$ regularity with respect to coordinates and parameters, we show that when the sequence is picked within a certain uniform family the top characteristic exponent and generator of top Oseledets space of the cocycle are $C^{k-1}$ in parameters. As applications, we obtain a linear response formula for the equivariant measure associated with random products of uniformly expanding maps, and we study the regularity of the Hausdorff dimension of a repeller associated with random compositions of one-dimensional cookie-cutters.

Linear Response for a Family of Self-consistent Transfer Operators

We study a system of all-to-all weakly coupled uniformly expanding circle maps in the thermodynamic limit. The state of the system is described by a probability measure and its evolution is given by the action of a nonlinear operator, also called a self-consistent transfer operator. We prove that when the coupling is sufficiently small, the system has a unique stable state that satisfies a linear response formula when varying the coupling strength.

Symmetrized decomposition of the Kubo-Bastin formula

The Smrcka-Streda version of Kubo's linear response formula is widely used in the literature to compute non-equilibrium transport properties of heterostructures. It is particularly useful for the evaluation of intrinsic transport properties associated with the Berry curvature of the Bloch states, such as anomalous and spin Hall currents as well as the damping-like component of the spin-orbit torque. Here, we demonstrate in a very general way that the widely used decomposition of the Kubo-Bastin formula introduced by Smrcka and Streda contains an overlap, which has lead to widespread confusion in the literature regarding the Fermi surface and Fermi sea contributions. To remedy this pathology, we propose a new decomposition of the Kubo-Bastin formula based on the permutation properties of the correlation function and derive a new set of formulas, without an overlap, that provides direct access to the transport effects of interest. We apply these new formulas to selected cases and demonstrate that the Fermi sea and Fermi surface contributions can be uniquely addressed with our symmetrized approach.

Linear response for random dynamical systems

We study for the first time linear response for random compositions of maps, chosen independently according to a distribution P. We are interested in the following question: how does an absolutely continuous stationary measure (acsm) of a random system change when P changes smoothly to Pε? For a wide class of one dimensional random maps, we prove differentiability of acsm with respect to ε; moreover, we obtain a linear response formula. Our results cover random maps whose transfer operator does not necessarily admit a spectral gap. We apply our results to iid compositions, with respect to various distributions Pε, of uniformly expanding circle maps, Gauss-Rényi maps (random continued fractions) and Pomeau-Manneville maps. Our results yield an exact formula for the invariant density of random continued fractions; while for Pomeau-Manneville maps our results provide a precise relation between their linear response under certain random perturbations and their linear response under deterministic perturbations.

Analytical techniques for linear response formula of equilibrium states

In this paper, we expose some of the analytical framework used to prove the differentiability of equilibrium states probabilities with respect to the dynamics. We explore the relationship between Birkhoff metrics in cones and their corresponding anisotropic spaces and prove some useful folklore theorems. We apply that framework to revisit several examples recently studied in the research works about equilibrium states.

A modified time domain subspace method for nonlinear identification based on nonlinear separation strategy

Nonlinear factors existing in engineering structures have drawn considerable attention, and nonlinear identification is a competent technique to understand the dynamic characteristics of nonlinear structures. Therefore, in this paper, a novel nonlinear separation subspace identification (NSSI) algorithm based on subspace algorithm and nonlinear separation strategy is proposed to conduct nonlinear parameter identification of nonlinear structures. For the proposed NSSI algorithm, the low-level excitation test is firstly conducted to obtain the transfer matrix in the linear response formula. Then, the obtained transfer matrix is used in the high-level excitation test to calculate the nonlinear response part by the proposed nonlinear separation strategy, and the subspace algorithm is utilized to identify the nonlinear parameter on the modified state-space model including only the nonlinear part. The proposed NSSI algorithm can reduce the coupling error caused by simultaneously processing both the large number part (corresponding to the linear part) and small number part (corresponding to the nonlinear part) in the traditional nonlinear subspace identification (NSI) algorithm. At last, two numerical experiments are given to validate the effectiveness of the developed novel nonlinear identification method. Furthermore, some influence factors are discussed to show the stability of the identification algorithm, and some comparisons between the proposed NSSI method and traditional NSI method are also conducted to demonstrate the advantages of the novel method.

Polaron mobility obtained by a variational approach for lattice Fröhlich models

Charge carrier mobility for a class of lattice models with long-range electron–phonon interaction was investigated. The approach for mobility calculation is based on a suitably chosen unitary transformation of the model Hamiltonian which transforms it into the form where the remaining interaction part can be treated as a perturbation. Relevant spectral functions were then obtained using Matsubara Green’s functions technique and charge carrier mobility was evaluated using Kubo’s linear response formula. Numerical results were presented for a wide range of electron–phonon interaction strengths and temperatures in the case of one-dimensional version of the model. The results indicate that the mobility decreases with increasing temperature for all electron–phonon interaction strengths in the investigated range, while longer interaction range leads to more mobile carriers.

Adiabatic Theorem for Quantum Spin Systems.

The first proof of the quantum adiabatic theorem was given as early as 1928. Today, this theorem is increasingly applied in a many-body context, e.g., in quantum annealing and in studies of topological properties of matter. In this setup, the rate of variation ϵ of local terms is indeed small compared to the gap, but the rate of variation of the total, extensive Hamiltonian, is not. Therefore, applications to many-body systems are not covered by the proofs and arguments in the literature. In this Letter, we prove a version of the adiabatic theorem for gapped ground states of interacting quantum spin systems, under assumptions that remain valid in the thermodynamic limit. As an application, we give a mathematical proof of Kubo's linear response formula for a broad class of gapped interacting systems. We predict that the density of nonadiabatic excitations is exponentially small in the driving rate and the scaling of the exponent depends on the dimension.

Linear and nonlinear response of the Vlasov system with nonintegrable Hamiltonian.

Linear and nonlinear response formulas taking into account all Casimir invariants are derived without use of angle-action variables of a single-particle (mean-field) Hamiltonian. This article deals mainly with the Vlasov system in a spatially inhomogeneous quasistationary state whose associating single-particle Hamiltonian is not integrable and has only one integral of the motion, the Hamiltonian itself. The basic strategy is to restrict the form of perturbation so that it keeps Casimir invariants within a linear order, and the single particle's probabilistic density function is smooth with respect to the single particle's Hamiltonian. The theory is applied for a spatially two-dimensional system and is confirmed by numerical simulations. A nonlinear response formula is also derived in a similar manner.