-
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.
-
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.
-
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.
-
Semiclassical theory of the magnetization process of the triangular lattice Heisenberg model
Authors:
Tommaso Coletta,
Tamás Tóth,
Karlo Penc,
Frédéric Mila
Abstract:
Motivated by the numerous examples of 1/3 magnetization plateaux in the triangular lattice Heisenberg an- tiferromagnet with spins ranging from 1/2 to 5/2, we revisit the semiclassical calculation of the magnetization curve of that model, with the aim of coming up with a simple method that allows one to calculate the full mag- netization curve, and not just the critical fields of the 1/3 plateau.…
▽ More
Motivated by the numerous examples of 1/3 magnetization plateaux in the triangular lattice Heisenberg an- tiferromagnet with spins ranging from 1/2 to 5/2, we revisit the semiclassical calculation of the magnetization curve of that model, with the aim of coming up with a simple method that allows one to calculate the full mag- netization curve, and not just the critical fields of the 1/3 plateau. We show that it is actually possible to calculate the magnetization curve including the first quantum corrections and the appearance of the 1/3 plateau entirely within linear spin-wave theory, with predictions for the critical fields that agree to order 1/S with those derived a long-time ago on the basis of arguments that required to go beyond linear spin-wave theory. This calculation relies on the central observation that there is a kink in the semiclassical energy at the field where the classical ground state is the collinear up-up-down structure, and that this kink gives rise to a locally linear behavior of the energy with the field when all semiclassical ground states are compared to each other for all fields. The magnetization curves calculated in this way for spin 1/2, 1 and 5/2 are shown to be in good agreement with available experimental data.
△ Less
Submitted 16 September, 2016; v1 submitted 20 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.
-
Dynamical structure factors and excitation modes of the bilayer Heisenberg model
Authors:
M. Lohöfer,
T. Coletta,
D. G. Joshi,
F. F. Assaad,
M. Vojta,
S. Wessel,
F. Mila
Abstract:
Using quantum Monte Carlo simulations along with higher-order spin-wave theory, bond-operator and strong-coupling expansions, we analyse the dynamical spin structure factor of the spin-half Heisenberg model on the square-lattice bilayer. We identify distinct contributions from the low-energy Goldstone modes in the magnetically ordered phase and the gapped triplon modes in the quantum disordered ph…
▽ More
Using quantum Monte Carlo simulations along with higher-order spin-wave theory, bond-operator and strong-coupling expansions, we analyse the dynamical spin structure factor of the spin-half Heisenberg model on the square-lattice bilayer. We identify distinct contributions from the low-energy Goldstone modes in the magnetically ordered phase and the gapped triplon modes in the quantum disordered phase. In the antisymmetric (with respect to layer inversion) channel, the dynamical spin structure factor exhibits a continuous evolution of spectral features across the quantum phase transition, connecting the two types of modes. Instead, in the symmetric channel we find a depletion of the spectral weight when moving from the ordered to the disordered phase. While the dynamical spin structure factor does not exhibit a well-defined distinct contribution from the amplitude (or Higgs) mode in the ordered phase, we identify an only marginally-damped amplitude mode in the dynamical singlet structure factor, obtained from interlayer bond correlations, in the vicinity of the quantum critical point. These findings provide quantitative information in direct relation to possible neutron or light scattering experiments in a fundamental two-dimensional quantum-critical spin system.
△ Less
Submitted 31 August, 2015;
originally announced August 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.
-
Semiclassical evidence of columnar order in the fully frustrated transverse field Ising model on the square lattice
Authors:
Tommaso Coletta,
Sergey E. Korshunov,
Frédéric Mila
Abstract:
We investigate the zero-temperature phase diagram of the fully frustrated transverse field Ising model on the square lattice both in the classical limit and in the presence of quantum fluctuations. At the classical level (the limit of infinite spin $S$), we find that upon decreasing the transverse field $Γ$ this model exhibits a phase transition from the fully polarized state into an eight-fold de…
▽ More
We investigate the zero-temperature phase diagram of the fully frustrated transverse field Ising model on the square lattice both in the classical limit and in the presence of quantum fluctuations. At the classical level (the limit of infinite spin $S$), we find that upon decreasing the transverse field $Γ$ this model exhibits a phase transition from the fully polarized state into an eight-fold degenerate translational symmetry breaking state. This phase can be identified to correspond to plaquette order in the dimer language and remains the lowest-energy state in the entire range of fields below the critical one, $Γ_c$. The eight-fold degenerate solution which corresponds to columnar order in the dimer language is a saddle point of the classical energy. It is degenerate with the plaquette solution at $Γ=0$ and is only slightly higher in energy in the whole interval $0<Γ<Γ_c$. The effect of quantum fluctuations is investigated in the context of a large S expansion both for the plaquette and columnar structures. For this purpose we employ an approximate method allowing to estimate from above the fluctuation-induced correction to the energy of a configuration which at the classical level is a saddle point of the energy, \textit{not} a local minimum. Although the convergence of the $1/S$ expansion in the $Γ/J\rightarrow 0$ limit remains an open question, harmonic quantum fluctuations show a clear tendency to overcome the energy difference between the two states and to change the classical picture favoring the columnar order over the plaquette one in a wide parameter range.
△ Less
Submitted 21 November, 2014; v1 submitted 15 August, 2014;
originally announced August 2014.
-
Quantum stabilization of classically unstable plateau structures
Authors:
T. Coletta,
M. E. Zhitomirsky,
F. Mila
Abstract:
Motivated by the intriguing report, in some frustrated quantum antiferromagnets, of magnetization plateaus whose simple collinear structure is {\it not} stabilized by an external magnetic field in the classical limit, we develop a semiclassical method to estimate the zero-point energy of collinear configurations even when they do not correspond to a local minimum of the classical energy. For the s…
▽ More
Motivated by the intriguing report, in some frustrated quantum antiferromagnets, of magnetization plateaus whose simple collinear structure is {\it not} stabilized by an external magnetic field in the classical limit, we develop a semiclassical method to estimate the zero-point energy of collinear configurations even when they do not correspond to a local minimum of the classical energy. For the spin-1/2 frustrated square-lattice antiferromagnet, this approach leads to the stabilization of a large 1/2 plateau with "up-up-up-down" structure for J_2/J_1>1/2, in agreement with exact diagonalization results, while for the spin-1/2 anisotropic triangular antiferromagnet, it predicts that the 1/3 plateau with "up-up-down" structure is stable far from the isotropic point, in agreement with the properties of Cs_2CuBr_4.
△ Less
Submitted 13 December, 2012;
originally announced December 2012.
-
Evidence of columnar order in the fully frustrated transverse field Ising model on the square lattice
Authors:
Sandro Wenzel,
Tommaso Coletta,
Sergey E. Korshunov,
Frédéric Mila
Abstract:
Using extensive classical and quantum Monte Carlo simulations, we investigate the ground-state phase diagram of the fully frustrated transverse field Ising model on the square lattice. We show that pure columnar order develops in the low-field phase above a surprisingly large length scale, below which an effective U(1) symmetry is present. The same conclusion applies to the Quantum Dimer Model wit…
▽ More
Using extensive classical and quantum Monte Carlo simulations, we investigate the ground-state phase diagram of the fully frustrated transverse field Ising model on the square lattice. We show that pure columnar order develops in the low-field phase above a surprisingly large length scale, below which an effective U(1) symmetry is present. The same conclusion applies to the Quantum Dimer Model with purely kinetic energy, to which the model reduces in the zero-field limit, as well as to the stacked classical version of the model. By contrast, the 2D classical version of the model is shown to develop plaquette order. Semiclassical arguments show that the transition from plaquette to columnar order is a consequence of quantum fluctuations.
△ Less
Submitted 27 November, 2012; v1 submitted 6 July, 2012;
originally announced July 2012.
-
Semiclassical approach to ground-state properties of hard-core bosons in two dimensions
Authors:
Tommaso Coletta,
Nicolas Laflorencie,
Frédéric Mila
Abstract:
Motivated by some inconsistencies in the way quantum fluctuations are included beyond the classical treatment of hard-core bosons on a lattice in the recent literature, we revisit the large-S semi-classical approach to hard-core bosons on the square lattice at T=0. First of all, we show that, if one stays at the purely harmonic level, the only correct way to get the 1/S correction to the density i…
▽ More
Motivated by some inconsistencies in the way quantum fluctuations are included beyond the classical treatment of hard-core bosons on a lattice in the recent literature, we revisit the large-S semi-classical approach to hard-core bosons on the square lattice at T=0. First of all, we show that, if one stays at the purely harmonic level, the only correct way to get the 1/S correction to the density is to extract it from the derivative of the ground state energy with respect to the chemical potential, and that to extract it from a calculation of the ground state expectation value of the particle number operator, it is necessary to include 1/\sqrt{S} corrections to the harmonic ground state. Building on this alternative approach to get 1/S corrections, we provide the first semiclassical derivation of the momentum distribution, and we revisit the calculation of the condensate density. The results of these as well as other physically relevant quantities such as the superfluid density are systematically compared to quantum Monte Carlo simulations. This comparison shows that the logarithmic corrections in the dilute Bose gas limit are only captured by the semi-classical approach if the 1/S corrections are properly calculated, and that the semi-classical approach is able to reproduce the 1/k divergence of the momentum distribution at k=0. Finally, the effect of 1/S^2 corrections is briefly discussed.
△ Less
Submitted 24 July, 2012; v1 submitted 23 December, 2011;
originally announced December 2011.
-
Phase diagram of the fully frustrated transverse-field Ising model on the honeycomb lattice
Authors:
T. Coletta,
J. -D. Picon,
S. E. Korshunov,
F. Mila
Abstract:
Motivated by the current interest in the quantum dimer model on the triangular lattice, we investigate the phase diagram of the closely related fully-frustrated transverse field Ising model on the honeycomb lattice using classical and semi-classical approximations. We show that, in addition to the fully polarized phase at large field, the classical model possesses a multitude of phases that break…
▽ More
Motivated by the current interest in the quantum dimer model on the triangular lattice, we investigate the phase diagram of the closely related fully-frustrated transverse field Ising model on the honeycomb lattice using classical and semi-classical approximations. We show that, in addition to the fully polarized phase at large field, the classical model possesses a multitude of phases that break the translational symmetry which in the dimer language, correspond to a plaquette phase and a columnar phase separated by an infinite cascade of mixed phases. The modification of the phase diagram by quantum fluctuations has been investigated in the context of linear spin-wave theory. The extrapolation of the semiclassical energies suggests that the plaquette phase extends down to zero field for spin 1/2, in agreement with the $\sqrt{12}\times\sqrt{12}$ phase of the quantum dimer model on the triangular lattice with only kinetic energy.
△ Less
Submitted 8 December, 2010; v1 submitted 28 June, 2010;
originally announced June 2010.
-
Frustration induced plateaux in S>=1/2 Heisenberg spin ladders
Authors:
F. Michaud,
T. Coletta,
S. R. Manmana,
J. -D. Picon,
F. Mila
Abstract:
We study the T = 0 magnetization of frustrated two-leg spin ladders with arbitrary value of the spin S. In the strong rung limit, we use degenerate perturbation theory to prove that frustration leads to magnetization plateaux at fractional values of the magnetization for all spins S, and to determine the critical ratios of parallel to diagonal inter-rung couplings for the appearance of these pla…
▽ More
We study the T = 0 magnetization of frustrated two-leg spin ladders with arbitrary value of the spin S. In the strong rung limit, we use degenerate perturbation theory to prove that frustration leads to magnetization plateaux at fractional values of the magnetization for all spins S, and to determine the critical ratios of parallel to diagonal inter-rung couplings for the appearance of these plateaux. These ratios depend both on the plateau and on the spin. To confirm these results, and to investigate the properties of these ladders away from the strong coupling limit, we have performed extensive density matrix renormalization group (DMRG) calculations for S <= 2. For large enough inter-rung couplings, all plateaux simply disappear, leading to a magnetization curve typical of integer-spin chains in a magnetic field. The intermediate region turns out to be surprisingly rich however, with, upon increasing the inter-rung couplings, the development of magnetization jumps and, in some cases, the appearance of one or more phase transitions inside a given plateau.
△ Less
Submitted 7 December, 2009; v1 submitted 9 July, 2009;
originally announced July 2009.