-
Transmission grid stability with large scale power flows due to renewables
Authors:
María Martínez-Barbeito,
Damià Gomila,
Julian Fritzsch,
Philippe Jacquod,
Pere Colet
Abstract:
We propose a general methodology for identifying critical lines in the long-distance transmission of power across large electrical networks. When the system is pushed to its limits of operation due to large power imbalances or contingencies, the network may lose synchrony. We analyze the Continental Europe grid as a case study and find that instabilities emerge due to topological constraints by wh…
▽ More
We propose a general methodology for identifying critical lines in the long-distance transmission of power across large electrical networks. When the system is pushed to its limits of operation due to large power imbalances or contingencies, the network may lose synchrony. We analyze the Continental Europe grid as a case study and find that instabilities emerge due to topological constraints by which the system loses its stable fixed point (synchronized state), causing the grid to split into two asynchronous zones. We discuss also how the modes of the system provide information on which areas are more susceptible to lose synchrony.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
A Novel Adaptive Inertia Strategy in Large-Scale Electric Power Grids
Authors:
Julian Fritzsch,
Philippe Jacquod
Abstract:
The increasing penetration of new renewable sources of energy in today's power grids is accompanied by a decrease in available electromechanical inertia. This leads to a reduced dynamical stability. To counterbalance this effect, virtual synchronous generators have been proposed to emulate conventional generators and provide inertia to power systems. The high flexibility of these devices makes it…
▽ More
The increasing penetration of new renewable sources of energy in today's power grids is accompanied by a decrease in available electromechanical inertia. This leads to a reduced dynamical stability. To counterbalance this effect, virtual synchronous generators have been proposed to emulate conventional generators and provide inertia to power systems. The high flexibility of these devices makes it possible to control the synthetic inertia they provide and to have them operate even more efficiently than the electromechanical inertia they replace. Here, we propose a novel control scheme for virtual synchronous generators, where the amount of inertia provided is large at short times - thereby absorbing local faults and disturbances as efficiently as conventional generators - but decreases over a tunable time interval to prevent long-time coherent oscillations from setting in. This new model is used to investigate the effect of adaptive inertia on large-scale power grids. Our model outperforms conventional constant inertia in all scenarios and for all performance measures considered. We show how an optimized geographical distribution of adaptive inertia devices not only effectively absorbs local faults, but also significantly improves the damping of inter-area oscillations.
△ Less
Submitted 2 November, 2023;
originally announced November 2023.
-
On the Robustness of Democratic Electoral Processes to Computational Propaganda
Authors:
Glory M. Givi,
Robin Delabays,
Matthieu Jacquemet,
Philippe Jacquod
Abstract:
There is growing evidence of systematic attempts to influence democratic elections by controlled and digitally organized dissemination of fake news. This raises the question of the intrinsic robustness of democratic electoral processes against external influences. Particularly interesting is to identify the social characteristics of a voter population that renders it more resilient against opinion…
▽ More
There is growing evidence of systematic attempts to influence democratic elections by controlled and digitally organized dissemination of fake news. This raises the question of the intrinsic robustness of democratic electoral processes against external influences. Particularly interesting is to identify the social characteristics of a voter population that renders it more resilient against opinion manipulation. Equally important is to determine which of the existing democratic electoral systems is more robust to external influences. Here we construct a mathematical electoral model to address these two questions. We find that, not unexpectedly, biased electorates with clear-cut elections are overall quite resilient against opinion manipulations, because inverting the election outcome requires to change the opinion of many voters. More interesting are unbiased or weakly biased electorates with close elections. We find that such populations are more resilient against opinion manipulations (i) if they are less polarized and (ii) when voters interact more with each other, regardless of their opinion differences, and that (iii) electoral systems based on proportional representation are generally the most robust. Our model qualitatively captures the volatility of the US House of Representatives elections. We take this as a solid validation of our approach.
△ Less
Submitted 2 January, 2024; v1 submitted 22 August, 2023;
originally announced August 2023.
-
Finite-time Correlations Boost Large Voltage-Angle Fluctuations in Electric Power Grids
Authors:
Melvyn Tyloo,
Jason Hindes,
Philippe Jacquod
Abstract:
Decarbonization in the energy sector has been accompanied by an increased penetration of new renewable energy sources in electric power systems. Such sources differ from traditional productions in that, first, they induce larger, undispatchable fluctuations in power generation and second, they lack inertia. Therefore, substituting new renewables for traditional generation induces stronger and more…
▽ More
Decarbonization in the energy sector has been accompanied by an increased penetration of new renewable energy sources in electric power systems. Such sources differ from traditional productions in that, first, they induce larger, undispatchable fluctuations in power generation and second, they lack inertia. Therefore, substituting new renewables for traditional generation induces stronger and more frequent disturbances and modifies the way disturbances propagate across AC electric power grids. Recent measurements have indeed reported long, non-Gaussian tails in the distribution of local grid-frequency data. Large frequency deviations may induce grid instabilities, leading in worst-case scenarios to cascading failures and large-scale blackouts. In this manuscript, we investigate how correlated noise disturbances, characterized by the cumulants of their distribution, propagate through meshed, high-voltage power grids. We show that for a single source of fluctuations, non-Gaussianities in the form of finite skewness and positive kurtosis of the noise distribution propagate over the entire network when the noise correlation time is larger than the network's intrinsic time scales, but that they vanish over short distances if the noise fluctuates rapidly. We furthermore show that a Berry-Esseen theorem leads to the vanishing of non-Gaussianities as the number of uncorrelated noise sources increases. Our results show that the persistence of non-Gaussian fluctuations of feed-in power have a global impact on power-grid dynamics when they fluctuate over time scales larger than the intrinsic time scales of the system, which, we argue, is the relevant regime in real power grids. Our predictions are corroborated by numerical simulations on realistic models of power grids.
△ Less
Submitted 1 March, 2022;
originally announced March 2022.
-
Long Wavelength Coherency in Well Connected Electric Power Networks
Authors:
Julian Fritzsch,
Philippe Jacquod
Abstract:
We investigate coherent oscillations in large scale transmission power grids, where large groups of generators respond in unison to a distant disturbance. Such long wavelength coherent phenomena are known as inter-area oscillations. Their existence in networks of weakly connected areas is well captured by singular perturbation theory. However, they are also observed in strongly connected networks…
▽ More
We investigate coherent oscillations in large scale transmission power grids, where large groups of generators respond in unison to a distant disturbance. Such long wavelength coherent phenomena are known as inter-area oscillations. Their existence in networks of weakly connected areas is well captured by singular perturbation theory. However, they are also observed in strongly connected networks without time-scale separation, where applying singular perturbation theory is not justified. We show that the occurrence of these oscillations is actually generic. Applying matrix perturbation theory, we show that, because these modes have the lowest oscillation frequencies of the system, they are only moderately sensitive to increased network connectivity between well chosen, initially weakly connected areas, and that their general structure remains the same, regardless of the strength of the inter-area coupling. This is qualitatively understood by bringing together the standard singular perturbation theory and Courant's nodal domain theorem.
△ Less
Submitted 1 March, 2022; v1 submitted 9 December, 2021;
originally announced December 2021.
-
Towards Model Reduction for Power System Transients with Physics-Informed PDE
Authors:
Laurent Pagnier,
Michael Chertkov,
Julian Fritzsch,
Philippe Jacquod
Abstract:
This manuscript reports the first step towards building a robust and efficient model reduction methodology to capture transient dynamics in a transmission level electric power system. Such dynamics is normally modeled on seconds-to-tens-of-seconds time scales by the so-called swing equations, which are ordinary differential equations defined on a spatially discrete model of the power grid. Followi…
▽ More
This manuscript reports the first step towards building a robust and efficient model reduction methodology to capture transient dynamics in a transmission level electric power system. Such dynamics is normally modeled on seconds-to-tens-of-seconds time scales by the so-called swing equations, which are ordinary differential equations defined on a spatially discrete model of the power grid. Following Seymlyen (1974) and Thorpe, Seyler, and Phadke (1999), we suggest to map the swing equations onto a linear, inhomogeneous Partial Differential Equation (PDE) of parabolic type in two space and one time dimensions with time-independent coefficients and properly defined boundary conditions. We illustrate our method on the synchronous transmission grid of continental Europe. We show that, when properly coarse-grained, i.e., with the PDE coefficients and source terms extracted from a spatial convolution procedure of the respective discrete coefficients in the swing equations, the resulting PDE reproduces faithfully and efficiently the original swing dynamics. We finally discuss future extensions of this work, where the presented PDE-based modeling will initialize a physics-informed machine learning approach for real-time modeling, $n-1$ feasibility assessment and transient stability analysis of power systems.
△ Less
Submitted 11 July, 2022; v1 submitted 26 October, 2021;
originally announced October 2021.
-
Matrix Perturbation Theory of Inter-Area Oscillations
Authors:
J. Fritzsch,
M. Tyloo,
Ph. Jacquod
Abstract:
Interconnecting power systems has a number of advantages such as better electric power quality, increased reliability of power supply, economies of scales through production and reserve pooling and so forth. Simultaneously, it may jeopardize the overall system stability with the emergence of so-called inter-area oscillations, which are coherent oscillations involving groups of rotating machines se…
▽ More
Interconnecting power systems has a number of advantages such as better electric power quality, increased reliability of power supply, economies of scales through production and reserve pooling and so forth. Simultaneously, it may jeopardize the overall system stability with the emergence of so-called inter-area oscillations, which are coherent oscillations involving groups of rotating machines separated by large distances up to thousands of kilometers. These often weakly damped modes may have harmful consequences for grid operation, yet despite decades of investigations, the mechanisms that generate them are still poorly understood, and the existing theories are based on assumptions that are not satisfied in real power grids where such modes are observed. Here we construct a matrix perturbation theory of large interconnected power systems that clarifies the origin and the conditions for the emergence of inter-area oscillations. We show that coherent inter-area oscillations emerge from the zero-modes of a multi-area network Laplacian matrix, which hybridize only weakly with other modes, even under significant capacity of the inter-area tie-lines, i.e. even when the standard assumption of area partitioning is not satisfied. The general theory is illustrated on a two-area system, and numerically applied to the well-connected PanTaGruEl model of the synchronous grid of continental Europe.
△ Less
Submitted 20 August, 2021;
originally announced August 2021.
-
Reconstructing Network Structures from Partial Measurements
Authors:
Melvyn Tyloo,
Robin Delabays,
Philippe Jacquod
Abstract:
The dynamics of systems of interacting agents is determined by the structure of their coupling network. The knowledge of the latter is, therefore, highly desirable, for instance, to develop efficient control schemes, to accurately predict the dynamics, or to better understand inter-agent processes. In many important and interesting situations, the network structure is not known, however, and previ…
▽ More
The dynamics of systems of interacting agents is determined by the structure of their coupling network. The knowledge of the latter is, therefore, highly desirable, for instance, to develop efficient control schemes, to accurately predict the dynamics, or to better understand inter-agent processes. In many important and interesting situations, the network structure is not known, however, and previous investigations have shown how it may be inferred from complete measurement time series on each and every agent. These methods implicitly presuppose that, even though the network is not known, all its nodes are. Here, we investigate the different problem of inferring network structures within the observed/measured agents. For symmetrically coupled dynamical systems close to a stable equilibrium, we establish analytically and illustrate numerically that velocity signal correlators encode not only direct couplings, but also geodesic distances in the coupling network within the subset of measurable agents. When dynamical data are accessible for all agents, our method is furthermore algorithmically more efficient than the traditional ones because it does not rely on matrix inversion.
△ Less
Submitted 12 October, 2021; v1 submitted 31 July, 2020;
originally announced July 2020.
-
Clusterization and phase diagram of the bimodal Kuramoto model with bounded confidence
Authors:
André Reggio,
Robin Delabays,
Philippe Jacquod
Abstract:
Inspired by the Deffuant and Hegselmann-Krause models of opinion dynamics, we extend the Kuramoto model to account for confidence bounds, i.e., vanishing interactions between pairs of oscillators when their phases differ by more than a certain value. We focus on Kuramoto oscillators with peaked, bimodal distribution of natural frequencies. We show that, in this case, the fixed-points for the exten…
▽ More
Inspired by the Deffuant and Hegselmann-Krause models of opinion dynamics, we extend the Kuramoto model to account for confidence bounds, i.e., vanishing interactions between pairs of oscillators when their phases differ by more than a certain value. We focus on Kuramoto oscillators with peaked, bimodal distribution of natural frequencies. We show that, in this case, the fixed-points for the extended model are made of certain numbers of independent clusters of oscillators, depending on the length of the confidence bound -- the interaction range -- and the distance between the two peaks of the bimodal distribution of natural frequencies. This allows us to construct the phase diagram of attractive fixed-points for the bimodal Kuramoto model with bounded confidence and to analytically explain clusterization in dynamical systems with bounded confidence.
△ Less
Submitted 27 September, 2020; v1 submitted 30 June, 2020;
originally announced July 2020.
-
Primary Control Effort under Fluctuating Power Generation in Realistic High-Voltage Power Networks
Authors:
Melvyn Tyloo,
Philippe Jacquod
Abstract:
Many recent works in control of electric power systems have investigated their synchronization through global performance metrics under external disturbances. The approach is motivated by fundamental changes in the operation of power grids, in particular by the substitution of conventional power plants with new renewable sources of electrical energy. This substitution will simultaneously increase…
▽ More
Many recent works in control of electric power systems have investigated their synchronization through global performance metrics under external disturbances. The approach is motivated by fundamental changes in the operation of power grids, in particular by the substitution of conventional power plants with new renewable sources of electrical energy. This substitution will simultaneously increase fluctuations in power generation and reduce the available mechanical inertia. It is crucial to understand how strongly these two evolutions will impact grid stability. With very few, mostly numerical exceptions, earlier works on performance metrics had to rely on unrealistic assumptions of grid homogeneity. Here we show that a modified spectral decomposition can tackle that issue in inhomogeneous power grids in cases where disturbances occur on time scales that are long compared to the intrinsic time scales of the grid. We find in particular that the magnitude of the transient excursion generated by disturbances with long characteristic times does not depend on inertia. For continental-size, high-voltage power grids, this corresponds to power fluctuations that are correlated on time scales of few seconds or more. We conclude that power fluctuations arising from new renewables will not require per se the deployment of additional rotational inertia. We numerically illustrate our results on the IEEE 118-Bus test case and a model of the synchronous grid of continental Europe.
△ Less
Submitted 30 June, 2020;
originally announced June 2020.
-
Eigenstate thermalization and ensemble equivalence in few-body fermionic systems
Authors:
Philippe Jacquod
Abstract:
We investigate eigenstate thermalization from the point of view of vanishing particle and heat currents between a few-body fermionic Hamiltonian prepared in one of its eigenstates and an external, weakly coupled Fermi-Dirac gas. The latter acts as a thermometric probe, with its temperature and chemical potential set so that there is neither particle nor heat current between the two subsystems. We…
▽ More
We investigate eigenstate thermalization from the point of view of vanishing particle and heat currents between a few-body fermionic Hamiltonian prepared in one of its eigenstates and an external, weakly coupled Fermi-Dirac gas. The latter acts as a thermometric probe, with its temperature and chemical potential set so that there is neither particle nor heat current between the two subsystems. We argue that the probe temperature can be attributed to the few-fermion eigenstate in the sense that (i) it varies smoothly with energy from eigenstate to eigenstate, (ii) it is equal to the temperature obtained from a thermodynamic relation in a wide energy range, (iii) it is independent of details of the coupling between the two systems in a finite parameter range, (iv) it satisfies the transitivity condition underlying the zeroth law of thermodynamics and (v) it is consistent with Carnot's theorem. These conditions are essentially independent of the interaction between the few fermions. When the interaction strength is weak, however, orbital occupancies in the few fermion system differ from the Fermi-Dirac distribution so that partial currents from or to the probe will eventually change its state. We find that (vi) above a certain critical interaction strength, orbital occupancies become close to the Fermi-Dirac distribution, leading to a true equilibrium between the few-fermion system and the probe. We conclude that for few-body systems with sufficiently strong interaction, the eigenstate thermalization hypothesis of Deutsch and Srednicki is complemented by ensemble equivalence: individual many-body eigenstates define a microcanonical ensemble that is equivalent to a canonical ensemble.
△ Less
Submitted 16 September, 2019;
originally announced September 2019.
-
Optimal placement of inertia and primary control : a matrix perturbation theory approach
Authors:
Laurent Pagnier,
Philippe Jacquod
Abstract:
The increasing penetration of inertialess new renewable energy sources reduces the overall mechanical inertia available in power grids and accordingly raises a number of issues of grid stability over short to medium time scales. It has been suggested that this reduction of overall inertia can be compensated to some extent by the deployment of substitution inertia - synthetic inertia, flywheels or…
▽ More
The increasing penetration of inertialess new renewable energy sources reduces the overall mechanical inertia available in power grids and accordingly raises a number of issues of grid stability over short to medium time scales. It has been suggested that this reduction of overall inertia can be compensated to some extent by the deployment of substitution inertia - synthetic inertia, flywheels or synchronous condensers. Of particular importance is to optimize the placement of the limited available substitution inertia, to mitigate voltage angle and frequency disturbances following a fault such as an abrupt power loss. Performance measures in the form of H2-norms have been recently introduced to evaluate the overall magnitude of such disturbances on an electric power grid. However, despite the mathematical conveniance of these measures, analytical results can be obtained only under rather restrictive assumptions of uniform damping ratio, or homogeneous distribution of inertia and/or primary control in the system. Here, we introduce matrix perturbation theory to obtain analytical results for optimal inertia and primary control placement where both are heterogeneous. Armed with that efficient tool, we construct two simple algorithms that independently determine the optimal geographical distribution of inertia and primary control. These algorithms are then implemented on a model of the synchronous transmission grid of continental Europe. We find that the optimal distribution of inertia is geographically homogeneous but that primary control should be mainly located on the slow modes of the network, where the intrinsic grid dynamics takes more time to damp frequency disturbances.
△ Less
Submitted 17 June, 2019;
originally announced June 2019.
-
Rate of change of frequency under line contingencies in high voltage electric power networks with uncertainties
Authors:
Robin Delabays,
Melvyn Tyloo,
Philippe Jacquod
Abstract:
In modern electric power networks with fast evolving operational conditions, assessing the impact of contingencies is becoming more and more crucial. Contingencies of interest can be roughly classified into nodal power disturbances and line faults. Despite their higher relevance, line contingencies have been significantly less investigated analytically than nodal disturbances. The main reason for…
▽ More
In modern electric power networks with fast evolving operational conditions, assessing the impact of contingencies is becoming more and more crucial. Contingencies of interest can be roughly classified into nodal power disturbances and line faults. Despite their higher relevance, line contingencies have been significantly less investigated analytically than nodal disturbances. The main reason for this is that nodal power disturbances are additive perturbations, while line contingencies are multiplicative perturbations, which modify the interaction graph of the network. They are, therefore, significantly more challenging to tackle analytically. Here, we assess the direct impact of a line loss by means of the maximal Rate of Change of Frequency (RoCoF) incurred by the system. We show that the RoCoF depends on the initial power flow on the removed line and on the inertia of the bus where it is measured. We further derive analytical expressions for the expectation and variance of the maximal RoCoF, in terms of the expectations and variances of the power profile in the case of power systems with power uncertainties. This gives analytical tools to identify the most critical lines in an electric power grid.
△ Less
Submitted 21 October, 2019; v1 submitted 13 June, 2019;
originally announced June 2019.
-
Global Robustness vs. Local Vulnerabilities in Complex Synchronous Networks
Authors:
Melvyn Tyloo,
Philippe Jacquod
Abstract:
In complex network-coupled dynamical systems, two questions of central importance are how to identify the most vulnerable components and how to devise a network making the overall system more robust to external perturbations. To address these two questions, we investigate the response of complex networks of coupled oscillators to local perturbations. We quantify the magnitude of the resulting excu…
▽ More
In complex network-coupled dynamical systems, two questions of central importance are how to identify the most vulnerable components and how to devise a network making the overall system more robust to external perturbations. To address these two questions, we investigate the response of complex networks of coupled oscillators to local perturbations. We quantify the magnitude of the resulting excursion away from the unperturbed synchronous state through quadratic performance measures in the angle or frequency deviations. We find that the most fragile oscillators in a given network are identified by centralities constructed from network resistance distances. Further defining the global robustness of the system from the average response over ensembles of homogeneously distributed perturbations, we find that it is given by a family of topological indices known as generalized Kirchhoff indices. Both resistance centralities and Kirchhoff indices are obtained from a spectral decomposition of the stability matrix of the unperturbed dynamics and can be expressed in terms of resistance distances. We investigate the properties of these topological indices in small-world and regular networks. In the case of oscillators with homogeneous inertia and damping coefficients, we find that inertia only has small effects on robustness of coupled oscillators. Numerical results illustrate the validity of the theory.
△ Less
Submitted 9 May, 2019;
originally announced May 2019.
-
Network desynchronization by non-Gaussian fluctuations
Authors:
Jason Hindes,
Philippe Jacquod,
Ira B. Schwartz
Abstract:
Many networks must maintain synchrony despite the fact that they operate in noisy environments. Important examples are stochastic inertial oscillators, which are known to exhibit fluctuations with broad tails in many applications, including electric power networks with renewable energy sources. Such non-Gaussian fluctuations can result in rare network desynchronization. Here we build a general the…
▽ More
Many networks must maintain synchrony despite the fact that they operate in noisy environments. Important examples are stochastic inertial oscillators, which are known to exhibit fluctuations with broad tails in many applications, including electric power networks with renewable energy sources. Such non-Gaussian fluctuations can result in rare network desynchronization. Here we build a general theory for inertial oscillator network desynchronization by non-Gaussian noise. We compute the rate of desynchronization and show that higher-moments of noise enter at specific powers of coupling: either speeding up or slowing down the rate exponentially depending on how noise statistics match the statistics of a network's slowest mode. Finally, we use our theory to introduce a technique that drastically reduces the effective description of network desynchronization. Most interestingly, when instability is associated with a single edge, the reduction is to one stochastic oscillator.
△ Less
Submitted 12 November, 2019; v1 submitted 27 April, 2019;
originally announced April 2019.
-
Noise-Induced Desynchronization and Stochastic Escape from Equilibrium in Complex Networks
Authors:
Melvyn Tyloo,
Robin Delabays,
Philippe Jacquod
Abstract:
Complex physical systems are unavoidably subjected to external environments not accounted for in the set of differential equations that models them. The resulting perturbations are standardly represented by noise terms. We derive conditions under which such noise terms perturb the dynamics strongly enough that they lead to stochastic escape from the initial basin of attraction of an initial stable…
▽ More
Complex physical systems are unavoidably subjected to external environments not accounted for in the set of differential equations that models them. The resulting perturbations are standardly represented by noise terms. We derive conditions under which such noise terms perturb the dynamics strongly enough that they lead to stochastic escape from the initial basin of attraction of an initial stable equilibrium state of the unperturbed system. Focusing on Kuramoto-like models we find in particular that, quite counterintuitively, systems with inertia leave their initial basin faster than or at the same time as systems without inertia, except for strong white-noise perturbations.
△ Less
Submitted 22 December, 2018;
originally announced December 2018.
-
The Key Player Problem in Complex Oscillator Networks and Electric Power Grids: Resistance Centralities Identify Local Vulnerabilities
Authors:
Melvyn Tyloo,
Laurent Pagnier,
Philippe Jacquod
Abstract:
Identifying key players in a set of coupled individual systems is a fundamental problem in network theory. Its origin can be traced back to social sciences and led to ranking algorithms based on graph theoretic centralities. Coupled dynamical systems differ from social networks in that, they are characterized by degrees of freedom with a deterministic dynamics and the coupling between individual u…
▽ More
Identifying key players in a set of coupled individual systems is a fundamental problem in network theory. Its origin can be traced back to social sciences and led to ranking algorithms based on graph theoretic centralities. Coupled dynamical systems differ from social networks in that, they are characterized by degrees of freedom with a deterministic dynamics and the coupling between individual units is a well-defined function of those degrees of freedom. One therefore expects the resulting coupled dynamics, and not only the network topology, to also determine the key players. Here, we investigate synchronizable network-coupled dynamical systems such as high voltage electric power grids and coupled oscillators. We search for nodes which, once perturbed by a local noisy disturbance, generate the largest overall transient excursion away from synchrony. A spectral decomposition of the coupling matrix leads to an elegant, concise, yet accurate solution to this identification problem. We show that, when the internodal coupling matrix is Laplacian, these key players are peripheral in the sense of a centrality measure defined from effective resistance distances. For linearly coupled dynamical systems such as weakly loaded power grids or consensus algorithms, the nodal ranking is efficiently obtained through a single Laplacian matrix inversion, regardless of the operational synchronous state. We call the resulting ranking index LRank. For heavily loaded power grids or coupled oscillators systems closer to the transition to synchrony, nonlinearities render the nodal ranking dependent on the synchronous state. In this case a weighted Laplacian matrix inversion gives another ranking index, which we call WLRank. Quite surprisingly, we find that LRank provides a faithful ranking even for well developed coupling nonlinearities, corresponding to oscillator angle differences up to $40^o$ approximately.
△ Less
Submitted 23 October, 2018;
originally announced October 2018.
-
Disturbance propagation, inertia location and slow modes in large-scale high voltage power grids
Authors:
Laurent Pagnier,
Philippe Jacquod
Abstract:
Conventional generators in power grids are steadily substituted with new renewable sources of electric power. The latter are connected to the grid via inverters and as such have little, if any rotational inertia. The resulting reduction of total inertia raises important issues of power grid stability, especially over short-time scales. We have constructed a model of the synchronous grid of contine…
▽ More
Conventional generators in power grids are steadily substituted with new renewable sources of electric power. The latter are connected to the grid via inverters and as such have little, if any rotational inertia. The resulting reduction of total inertia raises important issues of power grid stability, especially over short-time scales. We have constructed a model of the synchronous grid of continental Europe with which we numerically investigate frequency deviations as well as rates of change of frequency (RoCoF) following abrupt power losses. The magnitude of RoCoF's and frequency deviations strongly depend on the fault location, and we find the largest effects for faults located on the slowest mode - the Fiedler mode - of the network Laplacian matrix. This mode essentially vanishes over Belgium, Eastern France, Western Germany, northern Italy and Switzerland. Buses inside these regions are only weakly affected by faults occuring outside. Conversely, faults inside these regions have only a local effect and disturb only weakly outside buses. Following this observation, we reduce rotational inertia through three different procedures by either (i) reducing inertia on the Fiedler mode, (ii) reducing inertia homogeneously and (iii) reducing inertia outside the Fiedler mode. We find that procedure (iii) has little effect on disturbance propagation, while procedure (i) leads to the strongest increase of RoCoF and frequency deviations. These results for our model of the European transmission grid are corroborated by numerical investigations on the ERCOT transmission grid.
△ Less
Submitted 11 October, 2018;
originally announced October 2018.
-
The Kuramoto model on oriented and signed graphs
Authors:
Robin Delabays,
Philippe Jacquod,
Florian Dörfler
Abstract:
Many real-world systems of coupled agents exhibit directed interactions, meaning that the influence of an agent on another is not reciprocal. Furthermore, interactions usually do not have identical amplitude and/or sign. To describe synchronization phenomena in such systems, we use a generalized Kuramoto model with oriented, weighted and signed interactions. Taking a bottom-up approach, we investi…
▽ More
Many real-world systems of coupled agents exhibit directed interactions, meaning that the influence of an agent on another is not reciprocal. Furthermore, interactions usually do not have identical amplitude and/or sign. To describe synchronization phenomena in such systems, we use a generalized Kuramoto model with oriented, weighted and signed interactions. Taking a bottom-up approach, we investigate the simplest possible oriented networks, namely acyclic oriented networks and oriented cycles. These two types of networks are fundamental building blocks from which many general oriented networks can be constructed. For acyclic, weighted and signed networks, we are able to completely characterize synchronization properties through necessary and sufficient conditions, which we show are optimal. Additionally, we prove that if it exists, a stable synchronous state is unique. In oriented, weighted and signed cycles with identical natural frequencies, we show that the system globally synchronizes and that the number of stable synchronous states is finite.
△ Less
Submitted 22 January, 2019; v1 submitted 30 July, 2018;
originally announced July 2018.
-
Transient Performance of Electric Power Networks under Colored Noise
Authors:
T. Coletta,
B. Bamieh,
Ph. Jacquod
Abstract:
New classes of performance measures have been recently introduced to quantify the transient response to external disturbances of coupled dynamical systems on complex networks. These performance measures are time-integrated quadratic forms in the system's coordinates or their time derivative. So far, investigations of these performance measures have been restricted to Dirac-$δ$ impulse disturbances…
▽ More
New classes of performance measures have been recently introduced to quantify the transient response to external disturbances of coupled dynamical systems on complex networks. These performance measures are time-integrated quadratic forms in the system's coordinates or their time derivative. So far, investigations of these performance measures have been restricted to Dirac-$δ$ impulse disturbances, in which case they can be alternatively interpreted as giving the long time output variances for stochastic white noise power demand/generation fluctuations. Strictly speaking, the approach is therefore restricted to power fluctuating on time scales shorter than the shortest time scales in the swing equations. To account for power productions from new renewable energy sources, we extend these earlier works to the relevant case of colored noise power fluctuations, with a finite correlation time $τ> 0$. We calculate a closed-form expression for generic quadratic performance measures. Applied to specific cases, this leads to a spectral representation of performance measures as a sum over the non-zero modes of the network Laplacian. Our results emphasize the competition between inertia, damping and the Laplacian modes, whose balance is determined to a large extent by the noise correlation time scale $τ$.
△ Less
Submitted 31 January, 2019; v1 submitted 24 July, 2018;
originally announced July 2018.
-
Performance Measures in Electric Power Networks under Line Contingencies
Authors:
Tommaso Coletta,
Philippe Jacquod
Abstract:
Classes of performance measures expressed in terms of ${\cal H}_2$-norms have been recently introduced to quantify the response of coupled dynamical systems to external perturbations. So far, investigations of these performance measures have been restricted to nodal perturbations. Here, we go beyond these earlier works and consider the equally important, but so far neglected case of line perturbat…
▽ More
Classes of performance measures expressed in terms of ${\cal H}_2$-norms have been recently introduced to quantify the response of coupled dynamical systems to external perturbations. So far, investigations of these performance measures have been restricted to nodal perturbations. Here, we go beyond these earlier works and consider the equally important, but so far neglected case of line perturbations. We consider a network-reduced power system, where a Kron reduction has eliminated passive buses. Identifying the effect that a line fault in the physical network has on the Kron-reduced network, we find that performance measures depend on whether the faulted line connects two passive, two active buses or one active to one passive bus. In all cases, performance measures depend quadratically on the original load on the faulted line times a topology dependent factor. Our theoretical formalism being restricted to Dirac-$δ$ perturbations, we investigate numerically the validity of our results for finite-time line faults. We find good agreement with theoretical predictions for longer fault durations in systems with more inertia.
△ Less
Submitted 14 May, 2019; v1 submitted 28 November, 2017;
originally announced November 2017.
-
Robustness of Synchrony in Complex Networks and Generalized Kirchhoff Indices
Authors:
Melvyn Tyloo,
Tommaso Coletta,
Philippe Jacquod
Abstract:
In network theory, a question of prime importance is how to assess network vulnerability in a fast and reliable manner. With this issue in mind, we investigate the response to parameter changes of coupled dynamical systems on complex networks. We find that for specific, non-averaged perturbations, the response of synchronous states critically depends on the overlap between the perturbation vector…
▽ More
In network theory, a question of prime importance is how to assess network vulnerability in a fast and reliable manner. With this issue in mind, we investigate the response to parameter changes of coupled dynamical systems on complex networks. We find that for specific, non-averaged perturbations, the response of synchronous states critically depends on the overlap between the perturbation vector and the eigenmodes of the stability matrix of the unperturbed dynamics. Once averaged over properly defined ensembles of such perturbations, the response is given by new graph topological indices, which we introduce as generalized Kirchhoff indices. These findings allow for a fast and reliable method for assessing the specific or average vulnerability of a network against changing operational conditions, faults or external attacks.
△ Less
Submitted 26 February, 2018; v1 submitted 20 October, 2017;
originally announced October 2017.
-
Resistance distance criterion for optimal slack bus selection
Authors:
Tommaso Coletta,
Philippe Jacquod
Abstract:
We investigate the dependence of transmission losses on the choice of a slack bus in high voltage AC transmission networks. We formulate a transmission loss minimization problem in terms of slack variables representing the additional power injection that each generator provides to compensate the transmission losses. We show analytically that for transmission lines having small, homogeneous resista…
▽ More
We investigate the dependence of transmission losses on the choice of a slack bus in high voltage AC transmission networks. We formulate a transmission loss minimization problem in terms of slack variables representing the additional power injection that each generator provides to compensate the transmission losses. We show analytically that for transmission lines having small, homogeneous resistance over reactance ratios ${r/x\ll1}$, transmission losses are generically minimal in the case of a unique \textit{slack bus} instead of a distributed slack bus. For the unique slack bus scenario, to lowest order in ${r/x}$, transmission losses depend linearly on a resistance distance based indicator measuring the separation of the slack bus candidate from the rest of the network. We confirm these results numerically for several IEEE and Pegase testcases, and show that our predictions qualitatively hold also in the case of lines having inhomogeneous ${r/x}$ ratios, with optimal slack bus choices reducing transmission losses by ${10}\%$ typically.
△ Less
Submitted 7 July, 2017;
originally announced July 2017.
-
A predictive pan-European economic and production dispatch model for the energy transition in the electricity sector
Authors:
Laurent Pagnier,
Philippe Jacquod
Abstract:
The energy transition is well underway in most European countries. It has a growing impact on electric power systems as it dramatically modifies the way electricity is produced. To ensure a safe and smooth transition towards a pan-European electricity production dominated by renewable sources, it is of paramount importance to anticipate how production dispatches will evolve, to understand how incr…
▽ More
The energy transition is well underway in most European countries. It has a growing impact on electric power systems as it dramatically modifies the way electricity is produced. To ensure a safe and smooth transition towards a pan-European electricity production dominated by renewable sources, it is of paramount importance to anticipate how production dispatches will evolve, to understand how increased fluctuations in power generations can be absorbed at the pan-European level and to evaluate where the resulting changes in power flows will require significant grid upgrades. To address these issues, we construct an aggregated model of the pan-European transmission network which we couple to an optimized, few-parameter dispatch algorithm to obtain time- and geographically-resolved production profiles. We demonstrate the validity of our dispatch algorithm by reproducing historical production time series for all power productions in fifteen different European countries. Having calibrated our model in this way, we investigate future production profiles at later stages of the energy transition - determined by planned future production capacities - and the resulting interregional power flows. We find that large power fluctuations from increasing penetrations of renewable sources can be absorbed at the pan-European level via significantly increased electricity exchanges between different countries. We identify where these increased exchanges will require additional power transfer capacities. We finally introduce a physically-based economic indicator which allows to predict future financial conditions in the electricity market. We anticipate new economic opportunities for dam hydroelectricity and pumped-storage plants.
△ Less
Submitted 6 June, 2017;
originally announced June 2017.
-
The Size of the Sync Basin Revisited
Authors:
Robin Delabays,
Melvyn Tyloo,
Philippe Jacquod
Abstract:
In dynamical systems, the full stability of fixed point solutions is determined by their basin of attraction. Characterizing the structure of these basins is, in general, a complicated task, especially in high dimensionality. Recent works have advocated to quantify the non-linear stability of fixed points of dynamical systems through the relative volumes of the associated basins of attraction [D.…
▽ More
In dynamical systems, the full stability of fixed point solutions is determined by their basin of attraction. Characterizing the structure of these basins is, in general, a complicated task, especially in high dimensionality. Recent works have advocated to quantify the non-linear stability of fixed points of dynamical systems through the relative volumes of the associated basins of attraction [D. A. Wiley {\it et al.} Chaos {\bf 16}, 015103 (2006), P. J. Menck {\it et al.} Nat. Phys. {\bf 9}, 89 (2013)]. Here we revisit this issue and propose an efficient numerical method to estimate these volumes. The algorithm first identifies stable fixed points. Second, a set of initial conditions is considered that are randomly distributed at the surface of hypercubes centered on each fixed point. These initial conditions are dynamically evolved. The linear size of each basin of attraction is finally determined by the proportion of initial conditions which converge back to the fixed point. Armed with this algorithm, we revisit the problem considered by Wiley et al. in a seminal paper [D. A. Wiley {\it et al.} Chaos {\bf 16}, 015103 (2006)] that inspired the title of the present manuscript, and consider the equal-frequency Kuramoto model on a cycle. Fixed points of this model are characterized by an integer winding number $q$ and the number $n$ of oscillators. We find that the basin volumes scale as $(1-4q/n)^n$, contrasting with the Gaussian behavior postulated in Wiley et al.'s paper. Finally, we show the applicability of our method to complex models of coupled oscillators with different natural frequencies and on meshed networks.
△ Less
Submitted 14 November, 2017; v1 submitted 1 June, 2017;
originally announced June 2017.
-
How fast can one overcome the paradox of the energy transition? A physico-economic model for the European power grid
Authors:
Laurent Pagnier,
Philippe Jacquod
Abstract:
The paradox of the energy transition is that the low marginal costs of new renewable energy sources (RES) drag electricity prices down and discourage investments in flexible productions that are needed to compensate for the lack of dispatchability of the new RES. The energy transition thus discourages the investments that are required for its own harmonious expansion. To investigate how this parad…
▽ More
The paradox of the energy transition is that the low marginal costs of new renewable energy sources (RES) drag electricity prices down and discourage investments in flexible productions that are needed to compensate for the lack of dispatchability of the new RES. The energy transition thus discourages the investments that are required for its own harmonious expansion. To investigate how this paradox can be overcome, we argue that, under certain assumptions, future electricity prices are rather accurately modeled from the residual load obtained by subtracting non-flexible productions from the load. Armed with the resulting economic indicator, we investigate future revenues for European power plants with various degree of flexibility. We find that, if neither carbon taxes nor fuel prices change, flexible productions would be financially rewarded better and sooner if the energy transition proceeds faster but at more or less constant total production, i.e. by reducing the production of thermal power plants at the same rate as the RES production increases. Less flexible productions, on the other hand, would see their revenue grow more moderately. Our results indicate that a faster energy transition with a quicker withdrawal of thermal power plants would reward flexible productions faster.
△ Less
Submitted 15 June, 2018; v1 submitted 1 June, 2017;
originally announced June 2017.
-
Finite Size Scaling in the Kuramoto Model
Authors:
Tommaso Coletta,
Robin Delabays,
Philippe Jacquod
Abstract:
We investigate the scaling properties of the order parameter and the largest nonvanishing Lyapunov exponent for the fully locked state in the Kuramoto model with a finite number $N$ of oscillators. We show that, for any finite value of $N$, both quantities scale as $(K-K_L)^{1/2}$ with the coupling strength $K$ sufficiently close to the locking threshold $K_L$. We confirm numerically these predict…
▽ More
We investigate the scaling properties of the order parameter and the largest nonvanishing Lyapunov exponent for the fully locked state in the Kuramoto model with a finite number $N$ of oscillators. We show that, for any finite value of $N$, both quantities scale as $(K-K_L)^{1/2}$ with the coupling strength $K$ sufficiently close to the locking threshold $K_L$. We confirm numerically these predictions for oscillator frequencies evenly spaced in the interval $[-1, 1]$ and additionally find that the coupling range $δK$ over which this scaling is valid shrinks like $δK \sim N^{-α}$ with $α\approx1.5$ as $N \rightarrow \infty$. Away from this interval, the order parameter exhibits the infinite-$N$ behavior $r-r_L \sim (K-K_L)^{2/3}$ proposed by Pazó [Phys. Rev. E 72, 046211 (2005)]. We argue that the crossover between the two behaviors occurs because at the locking threshold, the upper bound of the continuous part of the spectrum of the fully locked state approaches zero as $N$ increases. Our results clarify the convergence to the $N \rightarrow \infty$ limit in the Kuramoto model.
△ Less
Submitted 27 April, 2017; v1 submitted 21 December, 2016;
originally announced December 2016.
-
Multistability of Phase-Locking in Equal-Frequency Kuramoto Models on Planar Graphs
Authors:
Robin Delabays,
Tommaso Coletta,
Philippe Jacquod
Abstract:
The number $\mathcal{N}$ of stable fixed points of locally coupled Kuramoto models depends on the topology of the network on which the model is defined. It has been shown that cycles in meshed networks play a crucial role in determining $\mathcal{N}$, because any two different stable fixed points differ by a collection of loop flows on those cycles. Since the number of different loop flows increas…
▽ More
The number $\mathcal{N}$ of stable fixed points of locally coupled Kuramoto models depends on the topology of the network on which the model is defined. It has been shown that cycles in meshed networks play a crucial role in determining $\mathcal{N}$, because any two different stable fixed points differ by a collection of loop flows on those cycles. Since the number of different loop flows increases with the length of the cycle that carries them, one expects $\mathcal{N}$ to be larger in meshed networks with longer cycles. Simultaneously, the existence of more cycles in a network means more freedom to choose the location of loop flows differentiating between two stable fixed points. Therefore, $\mathcal{N}$ should also be larger in networks with more cycles. We derive an algebraic upper bound for the number of stable fixed points of the Kuramoto model with identical frequencies, under the assumption that angle differences between connected nodes do not exceed $π/2$. We obtain $\mathcal{N}\leq\prod_{k=1}^c\left[2\cdot{\rm Int}(n_k/4)+1\right]$, which depends both on the number $c$ of cycles and on the spectrum of their lengths $\{n_k\}$. We further identify network topologies carrying stable fixed points with angle differences larger than $π/2$, which leads us to conjecture an upper bound for the number of stable fixed points for Kuramoto models on any planar network. Compared to earlier approaches that give exponential upper bounds in the total number of vertices, our bounds are much lower and therefore much closer to the true number of stable fixed points.
△ Less
Submitted 24 March, 2017; v1 submitted 8 September, 2016;
originally announced September 2016.
-
In memoriam: Markus Büttiker (1950 - 2013)
Authors:
Carlo Beenakker,
Philippe Jacquod,
Andrew Jordan,
Peter Samuelsson
Abstract:
Introduction for a special issue of Physica E dedicated to the memory of Markus Büttiker.
Introduction for a special issue of Physica E dedicated to the memory of Markus Büttiker.
△ Less
Submitted 1 September, 2017; v1 submitted 7 June, 2016;
originally announced June 2016.
-
Topologically Protected Loop Flows in High Voltage AC Power Grids
Authors:
T. Coletta,
R. Delabays,
I. Adagideli,
Ph. Jacquod
Abstract:
Geographical features such as mountain ranges or big lakes and inland seas often result in large closed loops in high voltage AC power grids. Sizable circulating power flows have been recorded around such loops, which take up transmission line capacity and dissipate but do not deliver electric power. Power flows in high voltage AC transmission grids are dominantly governed by voltage angle differe…
▽ More
Geographical features such as mountain ranges or big lakes and inland seas often result in large closed loops in high voltage AC power grids. Sizable circulating power flows have been recorded around such loops, which take up transmission line capacity and dissipate but do not deliver electric power. Power flows in high voltage AC transmission grids are dominantly governed by voltage angle differences between connected buses, much in the same way as Josephson currents depend on phase differences between tunnel-coupled superconductors. From this previously overlooked similarity we argue here that circulating power flows in AC power grids are analogous to supercurrents flowing in superconducting rings and in rings of Josephson junctions. We investigate how circulating power flows can be created and how they behave in the presence of ohmic dissipation. We show how changing operating conditions may generate them, how significantly more power is ohmically dissipated in their presence and how they are topologically protected, even in the presence of dissipation, so that they persist when operating conditions are returned to their original values. We identify three mechanisms for creating circulating power flows, (i) by loss of stability of the equilibrium state carrying no circulating loop flow, (ii) by tripping of a line traversing a large loop in the network and (iii) by reclosing a loop that tripped or was open earlier. Because voltage angles are uniquely defined, circulating power flows can take on only discrete values, much in the same way as circulation around vortices is quantized in superfluids.
△ Less
Submitted 30 September, 2016; v1 submitted 25 May, 2016;
originally announced May 2016.
-
Multistability of Phase-Locking and Topological Winding Numbers in Locally Coupled Kuramoto Models on Single-Loop Networks
Authors:
Robin Delabays,
Tommaso Coletta,
Philippe Jacquod
Abstract:
Determining the number of stable phase-locked solutions for locally coupled Kuramoto models is a long-standing mathematical problem with important implications in biology, condensed matter physics and electrical engineering among others. We investigate Kuramoto models on networks with various topologies and show that different phase-locked solutions are related to one another by loop currents. The…
▽ More
Determining the number of stable phase-locked solutions for locally coupled Kuramoto models is a long-standing mathematical problem with important implications in biology, condensed matter physics and electrical engineering among others. We investigate Kuramoto models on networks with various topologies and show that different phase-locked solutions are related to one another by loop currents. The latter take only discrete values, as they are characterized by topological winding numbers. This result is generically valid for any network, and also applies beyond the Kuramoto model, as long as the coupling between oscillators is antisymmetric in the oscillators' coordinates. Motivated by these results we further investigate loop currents in Kuramoto-like models. We consider loop currents in nonoriented $n$-node cycle networks with nearest-neighbor coupling. Amplifying on earlier works, we give an algebraic upper bound $\mathcal{N} \le 2 \, {\rm Int}[n/4]+1$ for the number $\cal N$ of different, linearly stable phase-locked solutions. We show that the number of different stable solutions monotonically decreases as the coupling strength is decreased. Furthermore stable solutions with a single angle difference exceeding $π/2$ emerge as the coupling constant $K$ is reduced, as smooth continuations of solutions with all angle differences smaller than $π/2$ at higher $K$. In a cycle network with nearest-neighbor coupling we further show that phase-locked solutions with two or more angle differences larger than $π/2$ are all linearly unstable. We point out similarities between loop currents and vortices in superfluids and superconductors as well as persistent currents in superconducting rings and two-dimensional Josephson junction arrays.
△ Less
Submitted 10 March, 2016; v1 submitted 14 December, 2015;
originally announced December 2015.
-
Linear Stability and the Braess Paradox in Coupled Oscillators Networks and Electric Power Grids
Authors:
Tommaso Coletta,
Philippe Jacquod
Abstract:
We investigate the influence that adding a new coupling has on the linear stability of the synchronous state in coupled oscillators networks. Using a simple model we show that, depending on its location, the new coupling can lead to enhanced or reduced stability. We extend these results to electric power grids where a new line can lead to four different scenarios corresponding to enhanced or reduc…
▽ More
We investigate the influence that adding a new coupling has on the linear stability of the synchronous state in coupled oscillators networks. Using a simple model we show that, depending on its location, the new coupling can lead to enhanced or reduced stability. We extend these results to electric power grids where a new line can lead to four different scenarios corresponding to enhanced or reduced grid stability as well as increased or decreased power flows. Our analysis shows that the Braess paradox may occur in any complex coupled system, where the synchronous state may be weakened and sometimes even destroyed by additional couplings.
△ Less
Submitted 5 April, 2016; v1 submitted 29 May, 2015;
originally announced May 2015.
-
Generation and Detection of Spin Currents in Semiconductor Nanostructures with Strong Spin-Orbit Interaction
Authors:
Fabrizio Nichele,
Szymon Hennel,
Patrick Pietsch,
Werner Wegscheider,
Peter Stano,
Philippe Jacquod,
Thomas Ihn,
Klaus Ensslin
Abstract:
Storing, transmitting, and manipulating information using the electron spin resides at the heart of spintronics. Fundamental for future spintronics applications is the ability to control spin currents in solid state systems. Among the different platforms proposed so far, semiconductors with strong spin-orbit interaction are especially attractive as they promise fast and scalable spin control with…
▽ More
Storing, transmitting, and manipulating information using the electron spin resides at the heart of spintronics. Fundamental for future spintronics applications is the ability to control spin currents in solid state systems. Among the different platforms proposed so far, semiconductors with strong spin-orbit interaction are especially attractive as they promise fast and scalable spin control with all-electrical protocols. Here we demonstrate both the generation and measurement of pure spin currents in semiconductor nanostructures. Generation is purely electrical and mediated by the spin dynamics in materials with a strong spin-orbit field. Measurement is accomplished using a spin-to-charge conversion technique, based on the magnetic field symmetry of easily measurable electrical quantities. Calibrating the spin-to-charge conversion via the conductance of a quantum point contact, we quantitatively measure the mesoscopic spin Hall effect in a multiterminal GaAs dot. We report spin currents of 174 pA, corresponding to a spin Hall angle of 34%.
△ Less
Submitted 18 May, 2015; v1 submitted 14 September, 2014;
originally announced September 2014.
-
Electric Field Control of Soliton Motion and Stacking in Trilayer Graphene
Authors:
Matthew Yankowitz,
Joel I-Jan Wang,
A. Glen Birdwell,
Yu-An Chen,
K. Watanabe,
T. Taniguchi,
Philippe Jacquod,
Pablo San-Jose,
Pablo Jarillo-Herrero,
Brian J. LeRoy
Abstract:
The crystal structure of a material plays an important role in determining its electronic properties. Changing from one crystal structure to another involves a phase transition which is usually controlled by a state variable such as temperature or pressure. In the case of trilayer graphene, there are two common stacking configurations (Bernal and rhombohedral) which exhibit very different electron…
▽ More
The crystal structure of a material plays an important role in determining its electronic properties. Changing from one crystal structure to another involves a phase transition which is usually controlled by a state variable such as temperature or pressure. In the case of trilayer graphene, there are two common stacking configurations (Bernal and rhombohedral) which exhibit very different electronic properties. In graphene flakes with both stacking configurations, the region between them consists of a localized strain soliton where the carbon atoms of one graphene layer shift by the carbon-carbon bond distance. Here we show the ability to move this strain soliton with a perpendicular electric field and hence control the stacking configuration of trilayer graphene with only an external voltage. Moreover, we find that the free energy difference between the two stacking configurations scales quadratically with electric field, and thus rhombohedral stacking is favored as the electric field increases. This ability to control the stacking order in graphene opens the way to novel devices which combine structural and electrical properties.
△ Less
Submitted 29 January, 2014;
originally announced January 2014.
-
Local Temperature of Out-of-Equilibrium Quantum Electron Systems
Authors:
J. Meair,
J. P. Bergfield,
C. A. Stafford,
Ph. Jacquod
Abstract:
We show how the local temperature of out-of-equilibrium, quantum electron systems can be consistently defined with the help of an external voltage and temperature probe. We determine sufficient conditions under which the temperature measured by the probe (i) is independent of details of the system-probe coupling, (ii) is equal to the temperature obtained from an independent current-noise measureme…
▽ More
We show how the local temperature of out-of-equilibrium, quantum electron systems can be consistently defined with the help of an external voltage and temperature probe. We determine sufficient conditions under which the temperature measured by the probe (i) is independent of details of the system-probe coupling, (ii) is equal to the temperature obtained from an independent current-noise measurement, (iii) satisfies the transitivity condition expressed by the zeroth law of thermodynamics, and (iv) is consistent with Carnot's theorem. This local temperature therefore characterizes the system in the common sense of equilibrium thermodynamics, but remains well defined even in out-of-equilibrium situations with no local equilibrium.
△ Less
Submitted 26 June, 2013;
originally announced June 2013.
-
Signatures of Majorana fermions in hybrid normal-superconducting rings
Authors:
Philippe Jacquod,
Markus Buttiker
Abstract:
We investigate persistent currents in metallic rings interrupted by a Coulomb blockaded topological superconducting segment. We show that the presence of Majorana bound states in the superconductor is reflected in the emergence of an h/e harmonics in the persistent current, whose sign is determined by the fermion parity of the superconductor. The Majorana bound states further render the current fi…
▽ More
We investigate persistent currents in metallic rings interrupted by a Coulomb blockaded topological superconducting segment. We show that the presence of Majorana bound states in the superconductor is reflected in the emergence of an h/e harmonics in the persistent current, whose sign is determined by the fermion parity of the superconductor. The Majorana bound states further render the current finite at zero flux, nevertheless the resulting peculiar symmetry of the persistent current is compatible with a free energy that is even in time-reversal symmetry breaking fields. These unique features of the persistent currents are robust against disorder and provide unambiguous signatures of the presence of Majorana fermions.
△ Less
Submitted 26 June, 2013;
originally announced June 2013.
-
Universal Features of Spin Transport and Breaking of Unitary Symmetries
Authors:
Ph. Jacquod,
I. Adagideli
Abstract:
When time-reversal symmetry is broken, quantum coherent systems with and without spin rotational symmetry exhibit the same universal behavior in their electric transport properties. We show that spin transport discriminates between these two cases. In systems with large charge conductance, spin transport is essentially insensitive to the breaking of time-reversal symmetry. However, in the opposite…
▽ More
When time-reversal symmetry is broken, quantum coherent systems with and without spin rotational symmetry exhibit the same universal behavior in their electric transport properties. We show that spin transport discriminates between these two cases. In systems with large charge conductance, spin transport is essentially insensitive to the breaking of time-reversal symmetry. However, in the opposite limit of a single exit channel, spin currents vanish identically in the presence of time-reversal symmetry, but are turned on by breaking it with an orbital magnetic field.
△ Less
Submitted 29 July, 2013; v1 submitted 22 January, 2013;
originally announced January 2013.
-
Mesoscopic behavior of the transmission phase through confined correlated electronic systems
Authors:
Rafael A. Molina,
Peter Schmitteckert,
Dietmar Weinmann,
Rodolfo A. Jalabert,
Philippe Jacquod
Abstract:
We investigate the effect of electronic correlations on the transmission phase of quantum coherent scatterers, considering quantum dots in the Coulomb blockade regime connected to two single-channel leads. We focus on transmission zeros and the associated π-phase lapses that have been observed in interferometric experiments. We numerically explore two types of models for quantum dots: (i) lattice…
▽ More
We investigate the effect of electronic correlations on the transmission phase of quantum coherent scatterers, considering quantum dots in the Coulomb blockade regime connected to two single-channel leads. We focus on transmission zeros and the associated π-phase lapses that have been observed in interferometric experiments. We numerically explore two types of models for quantum dots: (i) lattice models with up to eight sites, and (ii) resonant level models with up to six levels. We identify different regimes of parameters where the presence of electronic correlations is responsible for the increase or the decrease of the number of transmission zeros vs. electrochemical potential on the dot. However, we show that interaction effects cannot reproduce the universal behavior of alternating resonances and phase lapses, experimentally observed in many-electron Coulomb blockaded dots. Our numerical results strongly suggest that the main experimentally observed features are captured by the theory for chaotic ballistic dots of Molina et al., [Phys. Rev. Lett. 108, 076803 (2012)] incorporating one-particle wave-function correlations but ignoring many-particle electronic correlations.
△ Less
Submitted 24 July, 2013; v1 submitted 10 December, 2012;
originally announced December 2012.
-
Scattering Theory of Nonlinear Thermoelectricity in Quantum Coherent Conductors
Authors:
Jonathan Meair,
Philippe Jacquod
Abstract:
We construct a scattering theory of weakly nonlinear thermoelectric transport through sub-micron scale conductors. The theory incorporates the leading nonlinear contributions in temperature and voltage biases to the charge and heat currents. Because of the finite capacitances of sub-micron scale conducting circuits, fundamental conservation laws such as gauge invariance and current conservation re…
▽ More
We construct a scattering theory of weakly nonlinear thermoelectric transport through sub-micron scale conductors. The theory incorporates the leading nonlinear contributions in temperature and voltage biases to the charge and heat currents. Because of the finite capacitances of sub-micron scale conducting circuits, fundamental conservation laws such as gauge invariance and current conservation require special care to be preserved. We do this by extending the approach of Christen and Büttiker [Europhys. Lett. 35, 523 (1996)] to coupled charge and heat transport. In this way we write relations connecting nonlinear transport coefficients in a manner similar to Mott's relation between the linear thermopower and the linear conductance. We derive sum rules that nonlinear transport coefficients must satisfy to preserve gauge invariance and current conservation. We illustrate our theory by calculating the efficiency of heat engines and the coefficient of performance of thermoelectric refrigerators based on quantum point contacts and resonant tunneling barriers. We identify in particular rectification effects that increase device performance.
△ Less
Submitted 24 January, 2013; v1 submitted 23 November, 2012;
originally announced November 2012.
-
Suppression of Interactions in Multimode Random Lasers in the Anderson Localized Regime
Authors:
Peter Stano,
Philippe Jacquod
Abstract:
Understanding random lasing is a formidable theoretical challenge. Unlike conventional lasers, random lasers have no resonator to trap light, they are highly multimode with potentially strong modal interactions and they are based on disordered gain media, where photons undergo random multiple scattering. Interference effects notoriously modify the propagation of waves in such random media, but the…
▽ More
Understanding random lasing is a formidable theoretical challenge. Unlike conventional lasers, random lasers have no resonator to trap light, they are highly multimode with potentially strong modal interactions and they are based on disordered gain media, where photons undergo random multiple scattering. Interference effects notoriously modify the propagation of waves in such random media, but their fate in the presence of nonlinearity and interactions is poorly understood. Here, we present a semiclassical theory for multimode random lasing in the strongly scattering regime. We show that Anderson localization, a wave-interference effect, is not affected by the presence of nonlinearities. To the contrary, its presence suppresses interactions between simultaneously lasing modes. Using a recently constructed theory for complex multimode lasers, we show analytically how Anderson localization justifies a noninteracting, single-pole approximation. Consequently, lasing modes in a strongly scattering random laser are given by long-lived, Anderson localized modes of the passive cavity, whose frequency and wave profile does not vary with pumping, even in the multi-mode regime when mode overlap spatially.
△ Less
Submitted 24 October, 2012;
originally announced October 2012.
-
Onsager Relations in Coupled Electric, Thermoelectric and Spin Transport: The Ten-Fold Way
Authors:
Philippe Jacquod,
Robert S. Whitney,
Jonathan Meair,
Markus Büttiker
Abstract:
Hamiltonian systems can be classified into ten classes, in terms of the presence or absence of time-reversal symmetry, particle-hole symmetry and sublattice/chiral symmetry. We construct a quantum coherent scattering theory of linear transport for coupled electric, heat and spin transport; including the effect of Andreev reflection from superconductors. We derive a complete list of the Onsager rec…
▽ More
Hamiltonian systems can be classified into ten classes, in terms of the presence or absence of time-reversal symmetry, particle-hole symmetry and sublattice/chiral symmetry. We construct a quantum coherent scattering theory of linear transport for coupled electric, heat and spin transport; including the effect of Andreev reflection from superconductors. We derive a complete list of the Onsager reciprocity relations between transport coefficients for coupled electric, spin, thermoelectric and spin caloritronic effects. We apply these to all ten symmetry classes, paying special attention to specific additional relations that follow from the combination of symmetries, beyond microreversibility. We discuss these relations in several illustrative situations. We show the reciprocity between spin-Hall and inverse spin-Hall effects, and the reciprocity between spin-injection and magnetoelectric spin currents. We discuss the symmetry and reciprocity relations of Seebeck, Peltier, spin-Seebeck and spin-Peltier effects in systems with and without coupling to superconductors.
△ Less
Submitted 19 September, 2012; v1 submitted 6 July, 2012;
originally announced July 2012.
-
Macroscopic Coherent Rectification in Andreev Interferometers
Authors:
Jonathan Meair,
Philippe Jacquod
Abstract:
We investigate nonlinear transport through quantum coherent metallic conductors contacted to superconducting components. We find that in certain geometries, the presence of superconductivity generates a large, finite-average rectification effect. Specializing to Andreev interferometers, we show that the direction and magnitude of rectification can be controlled by a magnetic flux tuning the superc…
▽ More
We investigate nonlinear transport through quantum coherent metallic conductors contacted to superconducting components. We find that in certain geometries, the presence of superconductivity generates a large, finite-average rectification effect. Specializing to Andreev interferometers, we show that the direction and magnitude of rectification can be controlled by a magnetic flux tuning the superconducting phase difference at two contacts. In particular, this results in the breakdown of an Onsager reciprocity relation at finite bias. The rectification current is macroscopic in that it scales with the linear conductance, and we find that it exceeds 5% of the linear current at sub-gap biases of few tens of μeV's.
△ Less
Submitted 11 May, 2012;
originally announced May 2012.
-
Emergence of Superlattice Dirac Points in Graphene on Hexagonal Boron Nitride
Authors:
Matthew Yankowitz,
Jiamin Xue,
Daniel Cormode,
Javier D. Sanchez-Yamagishi,
K. Watanabe,
T. Taniguchi,
Pablo Jarillo-Herrero,
Philippe Jacquod,
Brian J. LeRoy
Abstract:
The Schrödinger equation dictates that the propagation of nearly free electrons through a weak periodic potential results in the opening of band gaps near points of the reciprocal lattice known as Brillouin zone boundaries. However, in the case of massless Dirac fermions, it has been predicted that the chirality of the charge carriers prevents the opening of a band gap and instead new Dirac points…
▽ More
The Schrödinger equation dictates that the propagation of nearly free electrons through a weak periodic potential results in the opening of band gaps near points of the reciprocal lattice known as Brillouin zone boundaries. However, in the case of massless Dirac fermions, it has been predicted that the chirality of the charge carriers prevents the opening of a band gap and instead new Dirac points appear in the electronic structure of the material. Graphene on hexagonal boron nitride (hBN) exhibits a rotation dependent Moiré pattern. In this letter, we show experimentally and theoretically that this Moiré pattern acts as a weak periodic potential and thereby leads to the emergence of a new set of Dirac points at an energy determined by its wavelength. The new massless Dirac fermions generated at these superlattice Dirac points are characterized by a significantly reduced Fermi velocity. The local density of states near these Dirac cones exhibits hexagonal modulations indicating an anisotropic Fermi velocity.
△ Less
Submitted 13 February, 2012;
originally announced February 2012.
-
Embedding method for the scattering phase in strongly correlated quantum dots
Authors:
Rafael A. Molina,
Peter Schmitteckert,
Dietmar Weinmann,
Rodolfo A. Jalabert,
Philippe Jacquod
Abstract:
The embedding method for the calculation of the conductance through interacting systems connected to single channel leads is generalized to obtain the full complex transmission amplitude that completely characterizes the effective scattering matrix of the system at the Fermi energy. We calculate the transmission amplitude as a function of the gate potential for simple diamond-shaped lattice models…
▽ More
The embedding method for the calculation of the conductance through interacting systems connected to single channel leads is generalized to obtain the full complex transmission amplitude that completely characterizes the effective scattering matrix of the system at the Fermi energy. We calculate the transmission amplitude as a function of the gate potential for simple diamond-shaped lattice models of quantum dots with nearest neighbor interactions. In our simple models we do not generally observe an interaction dependent change in the number of zeroes or phase lapses that depend only on the symmetry properties of the underlying lattice. Strong correlations separate and reduce the widths of the resonant peaks while preserving the qualitative properites of the scattering phase.
△ Less
Submitted 6 February, 2012;
originally announced February 2012.
-
Spin transistor action from Onsager reciprocity and SU(2) gauge theory
Authors:
I. Adagideli,
V. Lutsker,
M. Scheid,
Ph. Jacquod,
K. Richter
Abstract:
We construct a local gauge transformation to show how, in confined systems, a generic, weak nonhomogeneous SU(2) spin-orbit Hamiltonian reduces to two U(1) Hamiltonians for spinless fermions at opposite magnetic fields, to leading order in the spin-orbit strength. Using an Onsager relation, we further show how the resulting spin conductance vanishes in a two-terminal setup, and how it is turned on…
▽ More
We construct a local gauge transformation to show how, in confined systems, a generic, weak nonhomogeneous SU(2) spin-orbit Hamiltonian reduces to two U(1) Hamiltonians for spinless fermions at opposite magnetic fields, to leading order in the spin-orbit strength. Using an Onsager relation, we further show how the resulting spin conductance vanishes in a two-terminal setup, and how it is turned on by either weakly breaking time-reversal symmetry or opening additional transport terminals. We numerically check our theory for mesoscopic cavities as well as Aharonov-Bohm rings.
△ Less
Submitted 2 January, 2012;
originally announced January 2012.
-
Non-linear spin to charge conversion in mesoscopic structures
Authors:
Peter Stano,
Jaroslav Fabian,
Philippe Jacquod
Abstract:
Motivated by recent experiments [Vera-Marun et al., arXiv:1109.5969], we formulate a non-linear theory of spin transport in quantum coherent conductors. We show how a mesoscopic constriction with energy-dependent transmission can convert a spin current injected by a spin accumulation into an electric signal, relying neither on magnetic nor exchange fields. When the transmission through the constri…
▽ More
Motivated by recent experiments [Vera-Marun et al., arXiv:1109.5969], we formulate a non-linear theory of spin transport in quantum coherent conductors. We show how a mesoscopic constriction with energy-dependent transmission can convert a spin current injected by a spin accumulation into an electric signal, relying neither on magnetic nor exchange fields. When the transmission through the constriction is spin-independent, the spin-charge coupling is non-linear, with an electric signal that is quadratic in the accumulation. We estimate that gated mesoscopic constrictions have a sensitivity that allows to detect accumulations much smaller than a percent of the Fermi energy.
△ Less
Submitted 31 December, 2011;
originally announced January 2012.
-
Scattering phase of quantum dots: Emergence of universal behavior
Authors:
Rafael A. Molina,
Rodolfo A. Jalabert,
Dietmar Weinmann,
Philippe Jacquod
Abstract:
We investigate scattering through chaotic ballistic quantum dots in the Coulomb blockade regime. Focusing on the scattering phase, we show that large universal sequences emerge in the short wavelength limit, where phase lapses of $π$ systematically occur between two consecutive resonances. Our results are corroborated by numerics and are in qualitative agreement with existing experiments.
We investigate scattering through chaotic ballistic quantum dots in the Coulomb blockade regime. Focusing on the scattering phase, we show that large universal sequences emerge in the short wavelength limit, where phase lapses of $π$ systematically occur between two consecutive resonances. Our results are corroborated by numerics and are in qualitative agreement with existing experiments.
△ Less
Submitted 12 December, 2012; v1 submitted 2 November, 2011;
originally announced November 2011.
-
Conductance fluctuations in chaotic systems with tunnel barriers
Authors:
Daniel Waltner,
Jack Kuipers,
Philippe Jacquod,
Klaus Richter
Abstract:
Quantum effects are expected to disappear in the short-wavelength, semiclassical limit. As a matter of fact, recent investigations of transport through quantum chaotic systems have demonstrated the exponential suppression of the weak localization corrections to the conductance and of the Fano factor for shot-noise when the Ehrenfest time exceeds the electronic dwell time. On the other hand, conduc…
▽ More
Quantum effects are expected to disappear in the short-wavelength, semiclassical limit. As a matter of fact, recent investigations of transport through quantum chaotic systems have demonstrated the exponential suppression of the weak localization corrections to the conductance and of the Fano factor for shot-noise when the Ehrenfest time exceeds the electronic dwell time. On the other hand, conductance fluctuations, an effect of quantum coherence, retain their universal value in the limit of the ratio of Ehrenfest time over dwell time to infinity, when the system is ideally coupled to external leads. Motivated by this intriguing result we investigate conductance fluctuations through quantum chaotic cavities coupled to external leads via (tunnel) barriers of arbitrary transparency. Using the trajectory-based semiclassical theory of transport, we find a linear Ehrenfest time-dependence of the conductance variance showing a nonmonotonous, sinusoidal behavior as a function of the transperancy. Most notably, we find an increase of the conductance fluctuations with the Ehrenfest time, above their universal value, for the transparency less than 0.5. These results, confirmed by numerical simulations, show that, contrarily to the common wisdom, effects of quantum coherence may increase in the semiclassical limit, under special circumstances.
△ Less
Submitted 25 August, 2011;
originally announced August 2011.
-
Measuring Spin Accumulations with Current Noise
Authors:
Jonathan Meair,
Peter Stano,
Philippe Jacquod
Abstract:
We investigate the time-dependent fluctuations of the electric current injected from a reservoir with a non-equilibrium spin accumulation into a mesoscopic conductor. We show how the current noise power directly reflects the magnitude of the spin accumulation in two easily noticeable ways. First, as the temperature is lowered, the small-bias noise saturates at a value determined by the spin accumu…
▽ More
We investigate the time-dependent fluctuations of the electric current injected from a reservoir with a non-equilibrium spin accumulation into a mesoscopic conductor. We show how the current noise power directly reflects the magnitude of the spin accumulation in two easily noticeable ways. First, as the temperature is lowered, the small-bias noise saturates at a value determined by the spin accumulation. Second, in the presence of spin-orbit interactions in the conductor, the current noise exhibits a sample-dependent mesoscopic asymmetry under reversal of the electric current direction. These features provide for a purely electric protocol for measuring spin accumulations.
△ Less
Submitted 26 August, 2011; v1 submitted 12 April, 2011;
originally announced April 2011.
-
STM Spectroscopy of ultra-flat graphene on hexagonal boron nitride
Authors:
Jiamin Xue,
Javier Sanchez-Yamagishi,
D. Bulmash,
Philippe Jacquod,
A. Deshpande,
K. Watanabe,
T. Taniguchi,
Pablo Jarillo-Herrero,
B. J. LeRoy
Abstract:
Graphene has demonstrated great promise for future electronics technology as well as fundamental physics applications because of its linear energy-momentum dispersion relations which cross at the Dirac point. However, accessing the physics of the low density region at the Dirac point has been difficult because of the presence of disorder which leaves the graphene with local microscopic electron an…
▽ More
Graphene has demonstrated great promise for future electronics technology as well as fundamental physics applications because of its linear energy-momentum dispersion relations which cross at the Dirac point. However, accessing the physics of the low density region at the Dirac point has been difficult because of the presence of disorder which leaves the graphene with local microscopic electron and hole puddles, resulting in a finite density of carriers even at the charge neutrality point. Efforts have been made to reduce the disorder by suspending graphene, leading to fabrication challenges and delicate devices which make local spectroscopic measurements difficult. Recently, it has been shown that placing graphene on hexagonal boron nitride (hBN) yields improved device performance. In this letter, we use scanning tunneling microscopy to show that graphene conforms to hBN, as evidenced by the presence of Moire patterns in the topographic images. However, contrary to recent predictions, this conformation does not lead to a sizable band gap due to the misalignment of the lattices. Moreover, local spectroscopy measurements demonstrate that the electron-hole charge fluctuations are reduced by two orders of magnitude as compared to those on silicon oxide. This leads to charge fluctuations which are as small as in suspended graphene, opening up Dirac point physics to more diverse experiments than are possible on freestanding devices.
△ Less
Submitted 13 February, 2011;
originally announced February 2011.