이중 불리언 회로의 동기식 동역학과 흡인자 수식

초록

본 논문은 두 개의 사이드 회로가 하나의 노드를 공유하는 이중 불리언 자동회로(dbac)의 동기식(병렬) 업데이트에서 나타나는 주기적 궤도, 즉 attractor의 개수를 정확히 계산하는 식을 제시한다. 회로 길이와 부호(양·음 피드백)의 조합에 따라 유도된 수식은 임의의 주기에 대한 attractor 수와 전체 attractor 수를 모두 구할 수 있게 한다.

상세 분석

이 논문은 불리언 자동회로 네트워크 중에서도 특히 구조가 두 개의 순환 서브그래프(사이드 회로)와 이를 교차하는 하나의 공통 노드로 이루어진 ‘이중 불리언 자동회로(dbac)’에 초점을 맞춘다. 기존 연구에서는 단일 순환 회로의 동기식 동역학을 분석해 주기와 attractor 수를 구했지만, 두 회로가 공유 노드를 통해 결합될 경우 상호작용이 복합적으로 나타나 기존 방법을 그대로 적용하기 어렵다. 저자들은 먼저 dbac의 상호작용 그래프를 (L₁, L₂)라는 두 정수 쌍으로 파라미터화한다. 여기서 L₁, L₂는 각각 사이드 회로의 길이를 의미하고, 각 회로는 양성(+) 또는 음성(–) 피드백을 가질 수 있다. 회로의 부호는 노드 업데이트 규칙에 직접적인 영향을 미치며, 이는 전체 시스템의 상태 전이 그래프에서 허용되는 전이 패턴을 제한한다.

동기식 업데이트를 가정하면 전체 시스템은 2^{L₁+L₂−1}개의 가능한 상태 공간을 갖는다(공통 노드가 중복 계산되지 않음). 저자들은 이 상태 공간을 순환군(Cyclic group) 구조와 동형인 De Bruijn 그래프와 연결시켜, 각 주기 p에 대해 p‑주기 상태들의 존재 여부를 수론적 조건으로 환원한다. 구체적으로, 회로 길이와 부호에 따라 정의된 ‘유효 전이 행렬’의 고유값이 1이 되는 경우에만 p‑주기 attractor가 존재한다는 점을 이용한다. 이를 바탕으로 Möbius 역전(Möbius inversion) 공식을 적용해 p‑주기 attractor의 정확한 개수를

A(p) = (1/p) ∑_{d|p} μ(d) · N(p/d)

와 같이 표현한다. 여기서 μ는 Möbius 함수이며, N(k)는 길이 k에 대한 전이 행렬의 트레이스(또는 고정점 수)로, L₁, L₂, 부호에 따라 명시적인 다항식 형태로 계산된다. 특히, 양성 피드백 회로는 짝수 주기만 허용하고, 음성 피드백 회로는 홀수 주기를 강제하는 등, 부호 조합에 따라 N(k)의 값이 크게 달라진다.

전체 attractor 수는 모든 가능한 주기 p에 대해 A(p)를 합산한 값이며, 저자들은 이를 폐쇄형식으로 정리한다. 예를 들어, 두 회로가 모두 양성이고 길이가 서로소인 경우, 전체 attractor 수는 피보나치 수열과 유사한 재귀식으로 표현될 수 있다. 반면, 하나 이상의 회로가 음성 피드백을 포함하면, attractor 수는 짝·홀수 구분에 따라 두 개의 별도 합으로 나뉘어 계산된다. 이러한 결과는 기존의 단일 회로 분석 결과를 일반화한 것으로, dbac가 복합적인 피드백 구조를 가질 때도 정확히 attractor 수를 예측할 수 있음을 보여준다.

또한, 저자들은 수식의 복잡도를 분석하여, 주어진 (L₁, L₂)와 부호 조합에 대해 O(L₁ + L₂) 시간 안에 전체 attractor 수를 계산할 수 있는 알고리즘을 제시한다. 이는 대규모 네트워크 모델링에서 실용적인 의미를 가진다. 마지막으로, 몇 가지 구체적인 예시와 시뮬레이션 결과를 통해 이론적 수식이 실제 동기식 업데이트 시뮬레이션과 완벽히 일치함을 검증한다. 전체적으로, 이 연구는 복합 순환 구조를 갖는 불리언 네트워크의 동역학을 정량적으로 이해하는 데 중요한 이정표를 제공한다.