노드가 서로 다를 때 네트워크가 더 안정해지는 조건

Science ’Disorder-promoted stability’를 ML 연구자의 눈으로 읽기

Disorder-promoted stabilityScience 393(6817), 1241-1249

사이언스 브릿지
복잡계 물리
이질성이 집단 상태를 안정시키는 조건은 야코비안의 비에르미트성에 있다. 경사 흐름과 모멘텀의 차이로 옮겨 읽으면 ML에도 쓸 만한 교훈이 나온다.
공개

2026년 9월 25일

한 줄 요약과 선정 근거

노드의 내부 동역학이 2차원 이상이면, 노드끼리 서로 다를 때 네트워크의 집단 상태가 더 안정할 수 있다. 무작위로 흐트러뜨린 파라미터도 자주 안정성을 높인다.

Science는 같은 호(393권 6817호)에 레이사 드수자(Raissa D’Souza)의 해설 Perspective “Noise and diversity can boost stability”를 함께 실었다. 동기화, 전력망, 생태계에서 따로 보고되던 “이질성이 오히려 안정성을 높인다”는 관찰을 수학적 조건 하나로 묶어 냈다는 점에서 골랐다.

배경지식 세 가지

  1. 야코비안(Jacobian)과 최대 리아푸노프 지수(largest Lyapunov exponent): 평형점 근처의 작은 교란은 \(\delta\dot{x} = J\,\delta x\)를 따른다. \(J\)의 고윳값 실수부 중 최댓값 \(\Lambda_{\max}\)가 음수면 안정하고, 더 음수일수록 교란에서 빨리 회복한다. 논문이 쓰는 안정성 척도다.
  2. 에르미트(Hermitian) 행렬: 실수 행렬에서는 대칭 행렬을 말한다. 대칭 행렬의 고유벡터는 서로 직교한다. 비에르미트 행렬에서는 고유벡터가 기울어지다가 평행해질 수도 있다.
  3. 1차와 2차 노드 동역학: 쿠라모토(Kuramoto) 위상 진동자처럼 노드마다 변수 하나가 1차 방정식을 따르면 1차(1D)다. 전력망의 스윙 방정식처럼 \(\ddot{y}_i + b_i\dot{y}_i + \cdots\) 꼴의 뉴턴 방정식을 따르면 2차(2D)이고, \(b_i\)는 감쇠(damping)다.

무엇을 발견했나

쿠라모토 모형 이후의 통념은 “진동자가 비슷할수록 동기화가 쉽다”였다. 그런데 전력망, 레이저 배열, 전자 회로, 뉴런 회로에서는 이질성이 동기화 안정성을 높인다는 보고가 쌓여 왔다. 논문은 이 모순을 다음처럼 정리한다.

  • 정리: 야코비안이 파라미터 \(b\)에 대해 아핀(affine)이고 최적화 문제 \(\min_{b}\Lambda_{\max}(J(b))\)가 볼록(convex)이면, 전역 최적해 중 적어도 하나는 네트워크의 대칭을 모두 지킨다. 모든 노드가 구조적으로 같은 네트워크라면 그 해는 균일하다. 거꾸로 말해 이질적인 설정만이 최선이려면 문제가 비볼록이어야 한다.
  • 비볼록의 조건: 비볼록은 야코비안이 비에르미트일 때만 생기고, 넓은 부류의 비에르미트 시스템에서는 거의 모든 네트워크에서 생긴다. 2차 노드는 결합이 무방향이어도 야코비안이 비에르미트다.
  • 무질서의 효과: 방향 네트워크에서는 균일 최적점이 안장점(saddle point)이 될 수 있어, 무작위 섭동만으로 안정성이 좋아진다. 작은 섭동 가운데 안정성을 높인 비율은 방향 스몰월드(small-world) 네트워크에서 약 25%, 방향 순환(circulant) 네트워크에서 100%에 가까웠다.
  • 메커니즘: 균일할 때 서로 분리되어 있던 네트워크 모드가 이질성 때문에 섞인다(mode mixing). 무방향 네트워크에서 최적점에 도달하려면 야코비안의 고유벡터 두 개 이상이 평행해져야 한다.
  • 생태계: 일반화 로트카 볼테라(Lotka-Volterra) 모형에서 무질서는 경쟁 네트워크를 불안정하게, 상리공생(mutualistic) 네트워크를 안정하게 만들었다(분석한 네트워크에서 \(\Lambda_{\max}\) 다섯 배 개선). 상리공생이 우세하면 종 수와 연결 밀도가 클수록 이 효과가 잦아진다. “복잡할수록 불안정하다”는 May의 역설에 새 관점을 주는 결과다.
노드마다 상태 하나: 위상 φφ₁φ₂결합 kJ =−b₁−kkk−b₂−k대각선 양쪽의두 칸이 같다대칭 행렬(에르미트)

(가) 1차 노드: 야코비안이 대칭

노드마다 상태 둘: 위치와 속도y₁y₂b₁b₂결합 k속도속도J =00100001−kk−b₁0k−k0−b₂대각선 양쪽의두 칸이 다르다(1과 −k)비대칭 행렬(비에르미트)

(나) 2차 노드: 야코비안이 비대칭

높이는 Λmax, 낮을수록 안정균일 b₁ = b₂(b₁, b₂)(b₂, b₁)가운데균일이 최선

(다) 대칭이면 볼록: 균일한 설정이 최선

높이는 Λmax, 낮을수록 안정균일 b₁ = b₂더 안정한 이질적 설정균일 (안장점)

(라) 비볼록: 균일한 설정이 안장점일 수 있다

야코비안이 대칭이면 가장 안정한 설정 가운데 균일한 것이 늘 있고, 대칭이 아니면 없을 수 있다. 노드마다 상태가 하나인 1차 모형(가)은 결합이 양방향이면 야코비안이 대칭, 곧 에르미트다(b는 노드마다 따로 정하는 파라미터). 이때 Λmax는 (b1, b2)에 대해 볼록하다(다). 노드 번호만 바꾼 두 설정 (b1, b2)와 (b2, b1)은 높이가 같고, 둘을 잇는 점선의 한가운데 바로 아래가 균일한 설정이다. 볼록하면 이 점이 두 설정보다 높을 수 없다. 상태가 위치와 속도 둘인 2차 모형(나)은 결합이 양방향이어도 야코비안이 비대칭이고, 결합이 한 방향이면 1차 모형도 비대칭이 된다. 이때는 Λmax가 볼록하지 않을 수 있다(라). 방향 네트워크에서는 균일한 설정 가운데 가장 좋은 점이 대각선을 따라서는 가장 낮지만 가로질러서는 가장 높은 안장점이 되기도 한다. 그러면 대각선을 벗어나게 조금만 흔들어도 더 안정한 이질적 설정으로 내려가고(주황 화살표), 대각선을 따라 움직이면 올라간다(회색 화살표). 2차 모형을 무방향으로 이으면 최적점이 뾰족한 곳, 곧 미분할 수 없는 점에 놓여 경사를 따라 찾기 어렵다. 지형은 2차원으로 단순화한 개념도다.

어떻게 했나

  • 이론: 볼록 경우의 정리는 젠센 부등식(Jensen’s inequality)으로 증명하고, 비볼록성, 미분 가능성, 모드 섞임에 관한 명제는 보조 자료에서 증명한다.
  • 수치 실험 1: 2차 쿠라모토, 피츠휴 나구모(FitzHugh-Nagumo) 뉴런, 위상 진폭 진동자 등 여러 모형을 노드 100개짜리 무방향 전연결 네트워크(가중치 \(U[0,1]\))에서 최적화했다. 1차 모형만 균일한 설정에서 최적이었다.
  • 수치 실험 2: 방향 스몰월드, 척도 없는(scale-free) 네트워크, 실제 네트워크(예쁜꼬마선충 신경망, 먹이그물, 대장균 전사 조절망, 논리 회로 등)에서 균일 최적점 주변에 무작위 섭동을 1,000번씩 넣고 안정성이 좋아지는 비율을 셌다.
  • 수치 실험 3: 생태계 모형에서 무질서 세기와 상호작용 부호 비율을 바꿔 가며 네트워크 1,000개씩을 분석했다.
  • 새 물리 실험은 없다. 평형 상태 하나를 정해 두고, 그 선형 안정성을 이론과 시뮬레이션으로 따진 논문이다.

분야에서의 의미

  • 동기화 분석의 기본 도구인 마스터 안정성 함수(master stability function)는 노드가 같다는 전제에서 모드를 분리한다. 논문은 이 전제가 지우는 모드 섞임이 바로 안정화의 원천일 수 있음을 보였다.
  • 가장 넓게 적용될 결론은 방법론적 경고다. 다루기 쉽게 차원을 줄인 모형(1차 위상 모형 등)은 비에르미트 구조를 지워 이질성의 안정화 효과를 보지 못한다.
  • 설계 관점에서 이질성은 없앨 잡음이 아니라 쓸 수 있는 자원이 된다. 전력망의 감쇠 배분, 드론 편대, 레이저 배열이 응용 후보로 꼽힌다.
  • 계보: 같은 그룹의 “대칭 상태가 비대칭 시스템을 요구한다”(Nishikawa, Motter, PRL 2016)와 이를 실험으로 보인 연구(Molnar 외, Nature Physics 2020)의 연장선이다. 덧붙이자면, 균일한 감쇠가 최선이라는 추측의 반례는 파동 방정식 감쇠 최적화에서도 나온 적이 있다(Freitas, SIAM J. Control Optim. 1999). 이번 논문의 기여는 이런 관찰을 네트워크 동역학 일반의 조건으로 정리하고 무작위 이질성의 효과까지 보인 데 있다.

ML 연구자를 위한 인사이트

1. 경사 흐름과 모멘텀은 다른 부류다 (방법론 교훈, 직접 적용)

좌표별 가중치 감쇠 \(b_i\)를 둔 경사 흐름을 최솟점 근처에서 선형화하면 \(J = -(H + B)\)로, 논문의 1차 모형과 같은 꼴이다(\(H\)는 헤시안). 모멘텀(heavy ball) \(\ddot{\theta} + B\dot{\theta} + \nabla L(\theta) = 0\)을 선형화하면 아래처럼 논문의 2차 모형과 같은 꼴이 된다.

\[J = \begin{bmatrix} 0 & I \\ -H & -B \end{bmatrix}\]

1차 근사에서 얻은 “균일한 설정이 최선”이라는 직관을 모멘텀 동역학에 그대로 옮기면 안 되는 이유다. 다만 논문의 정리가 라플라시안이 아닌 헤시안에서도 그대로 성립하는지, 실제 학습에서 차이가 얼마나 되는지는 확인되지 않았다.

2. 뉴런의 이질성 (원리 차용, 영감 수준)

스파이킹 신경망에서 뉴런마다 시간 상수가 다르면 학습이 더 안정하고 성능이 좋아진다는 보고가 있다(Perez-Nieves 외, Nature Communications 2021). 이 논문은 그런 현상에 일반 조건 하나를 제시하는 셈이다. 유닛 내부 상태가 2차원 이상이거나 결합이 비대칭일 때 이질성이 안정화에 쓰일 수 있다. 2차 진동자 유닛을 쓰는 순환망(coRNN 등)은 감쇠를 모든 유닛이 공유하는데, 이 설계를 다시 볼 이유가 된다.

3. 대칭점은 안장점일 수 있다 (방법론 교훈, 영감 수준)

균일한 최적점이 안장점이면 작은 무작위 섭동이 개선을 준다. 은닉 유닛을 모두 같은 값으로 초기화하면 학습이 대칭에 갇히고 무작위 초기화가 그 대칭을 깨는 것과 구조가 닮았다. 또 무방향 네트워크의 최적점은 고유벡터가 평행해지는 미분 불가능한 점에 놓인다. 스펙트럼 반경 같은 안정성 지표를 직접 최적화하는 문제에서 경사법이 같은 벽을 만날 수 있다.

4. 생태계와 순환망 초기화 (새 문제와 데이터, 연결 약함)

큰 무작위 시스템의 안정성을 무작위 행렬로 따지는 May의 논증은 순환망 초기화의 “혼돈의 경계(edge of chaos)” 이론과 수학적 뿌리가 같다. 하지만 이번 생태계 결과가 ML에 무엇을 더하는지는 아직 분명하지 않다.

한계와 남은 질문

  • 평형 상태 하나를 정해 두고 그 선형 안정성만 주로 본다. 큰 교란이 들어올 때 어디까지 버티는지, 곧 끌개 영역(basin of attraction)은 보조 자료의 몇몇 예시에서만 다룬다.
  • 아무 이질성이나 돕지는 않는다. 논문의 2차 쿠라모토 진동자 쌍 예에서 이질성이 돕는 곳은 결합이 약한 영역에 한정됐고, 감쇠 차이가 커지면 오히려 균일할 때보다 나빠졌다. 한쪽 감쇠가 음수가 될 만큼 벌어지면 불안정해졌다.
  • 무작위 섭동이 도울 확률은 구조에 크게 좌우된다. 허브가 있는 척도 없는 네트워크에서는 더 낮았다.
  • 실제 네트워크 실험에서도 연결 구조만 가져오고 가중치는 무작위로 새로 매겼다. 실제 상호작용 세기에서도 같은 결과가 나오는지는 열려 있다.
−0.3−0.2−0.100.1−1−0.500.51감쇠 차이 Δb균일 −0.17가장 안정 −0.28불안정Λmax (낮을수록 안정)

(가) 결합 세기 K = 0.05일 때의 단면

00.050.1−1−0.500.51감쇠 차이 Δb0.07이질이 낫다균일이 낫다K = 0.05불안정불안정K

(나) 이질성이 돕는 곳과 해치는 곳

감쇠 차이가 두 진동자의 안정성을 어떻게 바꾸는지 논문 그림 1C의 예로 다시 계산했다. 고유 진동수가 같은 2차 쿠라모토 진동자 두 개를 양방향으로 잇고 위상이 같은 동기 상태 근처에서 선형화하면 ÿ1 + b1ẏ1 + K(y1 − y2) = 0, ÿ2 + b2ẏ2 + K(y2 − y1) = 0이다. 감쇠 평균은 0.75로 두고 b1 = 0.75 + Δb, b2 = 0.75 − Δb로 차이만 바꿨다. Λmax는 야코비안 고윳값 실수부의 최댓값이다. 두 위상을 똑같이 돌리는 방향은 늘 고윳값이 0이어서 뺐다. (가) K = 0.05일 때 Λmax는 균일하면 −0.17이고 Δb = ±0.23에서 −0.28로 가장 낮다. 이 바닥에서 고윳값 두 개가 겹치기 때문에 곡선이 뾰족하다. |Δb|가 0.44를 넘으면 균일할 때보다 나빠지고, 0.78를 넘으면 불안정하다. 이때 b2는 이미 음수다. (나) 이질성이 돕는 곳(주황)은 결합이 약해 균일한 쌍이 과감쇠인 곳(K < 0.75²/8 ≈ 0.07)에만 있고, 결합이 셀수록 좁아진다. 점선은 가장 좋은 Δb, 가로 파선은 (가)의 단면이다.

열린 질문

  • 모멘텀 계열 최적화기에서 파라미터 그룹마다 감쇠를 다르게 두면 실제 신경망 학습의 수렴 안정성은 어떻게 달라지는가?
  • 유닛 내부 상태가 2차원 이상인 순환망에서 유닛 파라미터의 이질성은 1차 유닛과 다르게 작동하는가?
  • 미분 불가능한 최적점과 안장점이 섞인 안정성 지표를 학습 목표로 쓸 때, 무작위 섭동이 경사법보다 나은 조건은 무엇인가?

이해 확인

1. 결합이 무방향이어도 2차 노드 네트워크에서는 이질성이 안정성을 높일 수 있다. 까닭은?

  1. 2차 노드는 변수가 둘이라 계산이 더 정확해서다
  2. 야코비안이 [[0, I], [−L, −B]] 꼴이라 비에르미트가 되고, 그러면 최적화가 비볼록해져 대칭을 깨는 설정이 최선일 수 있어서다
  3. 2차 노드에는 감쇠가 없어서다
  4. 무방향 결합은 늘 불안정해서다

1차 노드에 무방향 결합이면 야코비안이 대칭이라 안정성 최적화가 볼록하고, 구조가 같은 노드들에서는 균일한 설정이 최선이다. 2차 노드에서는 결합이 무방향이어도 야코비안이 대칭이 아니라 문제가 비볼록해질 수 있고, 그때 이질적인 설정이 균일한 설정을 이길 수 있다.

2. “작은 무작위 섭동이 안정성을 높인다”는 효과가 나타나려면 균일 최적점은 어떤 점이어야 하나?

  1. 전역 최솟점
  2. 미분 가능한 안장점
  3. 미분 불가능한 최솟점
  4. 불안정한 평형점

균일 최적점에서 Λmax가 미분 가능하고 그 점이 안장점이면, 방향을 고르지 않은 작은 섭동으로도 안정성이 나아진다. 방향 네트워크에서 이런 경우가 흔했다. 무방향 네트워크의 최적점은 고유벡터가 평행해지는 미분 불가능한 점에 놓여 경사법이나 작은 섭동으로 닿기 어렵다.

3. 생태계 모형에서 무질서의 효과로 맞는 것은?

  1. 경쟁 네트워크와 상리공생 네트워크 모두 불안정해졌다
  2. 경쟁 네트워크는 불안정해지고 상리공생 네트워크는 안정해졌다
  3. 경쟁 네트워크는 안정해지고 상리공생 네트워크는 불안정해졌다
  4. 종 수가 많을수록 효과가 사라졌다

무질서는 경쟁 네트워크를 불안정하게, 상리공생 네트워크를 안정하게 만들었다(분석한 사례에서 Λmax 다섯 배 개선). 상리공생이 우세하면 종 수와 연결 밀도가 클수록 이 효과가 잦아진다. 복잡할수록 불안정하다는 May의 예측과 달리, 복잡하면서도 안정한 시스템이 무질서한 네트워크 가운데 있을 수 있다.

논문 정보

  • 논문: Arthur N. Montanari, Pietro Zanin, Adilson E. Motter, “Disorder-promoted stability”, Science 393(6817), 1241-1249, 2026년 9월 17일. doi:10.1126/science.aeg3946