⚛️ quantum physics

Variational Quantum Algorithms for Hyperelasticity: Incorporating Nonlinear Constitutive Behavior

본 논문은 보조 변수와 페널티 제약 조건을 통해 신장(stretch)의 유리수 거듭제곱을 포함하는 비선형 구성 방정식이 포함된 1차원 비압축성 초탄성 문제를 해결하기 위해 변분 양자 알고리즘 프레임워크를 확장하며, 반복적인 보정 전략을 통해 솔루션의 정확도를 더욱 향상시킨다.

원저자: Uditnarayan Kouskiya, Caglar Oskay

게시일 2026-08-20

📖 4 분 읽기🧠 심층 분석

CC BY 4.0

원저자: Uditnarayan Kouskiya, Caglar Oskay

원본 논문은 CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

고무줄, 자동차 타이어, 우리 몸의 부드러운 조직과 같이 우리 주변의 재료들을 컴퓨터로 완벽하게 정밀하게 모델링할 수 있는 세상을 상상해 보십시오. 엔지니어와 과학자들은 재료가 압력을 받을 때 어떻게 늘어나고, 찌그러지고, 다시 원래대로 돌아오는지 예측하기 위해 오랫동안 고전 컴퓨터에 의존해 왔습니다. 그러나 이러한 재료들은 종종 복잡하고 비선형적인 방식으로 행동하는데, 이는 작은 충격이 항상 비례적인 움직임으로 이어지지는 않는다는 것을 의미합니다. 이러한 거동 뒤에 숨겨진 수학을 푸는 것은 매우 어려운 일이며, 특히 차세대 컴퓨터인 양자 컴퓨터에서 이를 시뮬레이션하려고 할 때 더욱 그러합니다. 양자 컴퓨터는 오늘날의 어떤 기계보다 특정 문제를 더 빠르게 해결할 수 있을 것으로 기대되지만, 실제 세계의 재료를 흥미롭게 만드는 바로 그 비선형적 동작에는 어려움을 겪습니다. 양자 컴퓨터는 직선과 단순한 관계를 다루는 데는 탁월하지만, 현실의 곡선과 뒤틀림에 직면하면 비틀거립니다.

이것이 바로 밴더빌트 대학교의 연구진인 우디타르나얀 코우스키야(Uditnarjan Kouskiya)와 카글라 오스카이(Caglar Oskay)가 해결하고자 했던 과제입니다. 이들은 토목 공학과 양자 컴퓨팅의 접점에서 연구하고 있으며, 이 분야는 구조 설계의 문제를 해결하기 위해 양자 역학의 기묘한 법칙을 사용하고자 합니다. 그들의 목표는 고무와 같은 재료의 복잡한 수학적 처리를 다루기 위해 이전에 개발했던 새로운 방법을 확장하는 것이었습니다. 구체적으로, 그들은 자신들의 접근 방식이 고무나 생체 조직의 모델에서 흔히 나타나는 분수 승(fractional powers)을 포함하는 복잡한 규칙을 따르는 재료에도 작동할 수 있는지 확인하고자 했습니다. 연구진은 이 연구를 위해 물리적인 양자 컴퓨터를 구축한 것이 아니라, 양자 기계가 어떻게 작동하는지를 모사하는 정교한 고전 컴퓨터 상에서 시뮬레이션을 실행했습니다. 그들의 작업은 특정 전략이 양자 알고리즘의 경직된 논리와 초탄성(hyperelastic) 재료의 유연한 현실 사이의 간극을 메울 수 있음을 보여주는 개념 증명 역할을 합니다.

연구진 작업의 핵심은 변분 양자 알고리즘(Variational Quantum Algorithm)이라 불리는 기술입니다. 간단히 말해, 이는 시스템의 가장 낮은 에너지 상태를 찾기 위해 양자 컴퓨터를 사용하는 방법으로, 이는 주어진 하중 하에서 재료가 취하게 될 가장 안정적인 형태에 해당합니다. 연구진은 재료의 에너지에 대한 수학적 설명을 '비용 함수(cost function)', 즉 양자 컴퓨터가 최소화하려고 노력하는 점수로 변환했습니다. 문제는 이러한 고무와 같은 재료의 방정식에는 분수나 음수의 거듭제곱이 포함되는 경우가 많은데, 현재의 양자 회로는 이를 직접 처리할 수 없다는 점입니다. 이를 해결하기 위해 팀은 영리한 우회 방법을 도입했습니다. 그들은 방정식의 까다로운 부분들을 대신할 일종의 도우미 숫자인 '보조 변수(auxiliary variables)'를 추가했습니다. 그런 다음 시스템에 '벌칙(penalty)'을 추가했습니다. 만약 도우미 숫자가 방정식의 까다로운 부분과 완벽하게 일치하지 않으면 비용 점수가 올라가도록 설정한 것입니다. 이를 통해 양자 컴퓨터가 보조 변수와 까다로운 수학적 부분이 서로 일치하도록 하는 해답을 찾게 함으로써, 복잡한 문제를 양자 기계가 이해할 수 있는 형태로 효과적으로 번역했습니다.

하지만 이 번역은 근사치이며, 모든 근사치가 그렇듯 약간의 오차를 발생시킵니다. 이를 해결하기 위해 연구진은 프로세스의 두 번째 단계인 반복 교정 전략을 개발했습니다. 양자 컴퓨터가 일단 초기적인 거친 해답을 찾아내면, 연구진은 그 결과를 사용하여 '교정(correction)' 필드를 계산합니다. 그런 다음 알고리즘을 다시 실행하여, 첫 번째 답을 개선하는 데 필요한 미세한 조정을 찾도록 요청합니다. 이 과정을 반복하면서 각 단계마다 해답을 정교하게 다듬어 가며 재료의 실제 물리적 거동에 더 가깝게 만듭니다. 이는 라디오를 튜닝하는 것과 비슷합니다. 먼저 일반적인 채널을 찾은 다음, 잡음이 사라지고 신호가 선명해질 때까지 다이얼을 미세하게 조정하는 과정입니다.

팀은 오그던(Ogden) 모델과 무니-리블린(Mooney–Rivlin) 모델을 포함하여 고무와 같은 재료의 세 가지 서로 다른 모델에 대해 이 2단계 접근 방식을 테스트했습니다. 이 모델들은 엔지니어들이 이러한 물질의 변형을 설명할 때 사용하는 표준적인 방식들입니다. 그들은 특정 힘과 경계 조건이 가해진 1밀리미터 길이의 1차원 재료 스트립을 시뮬레이션했습니다. 3개의 큐비트(양자 정보의 기본 단위)를 사용하는 시뮬레이션된 양자 시스템을 통해 알고리즘을 실행했습니다. 결과에 따르면, 초기 벌칙 기반 방식은 괜찮은 추정치를 제공했지만, 반복 교정 과정을 거치면서 정확도가 크게 향상되었습니다. 무니-리블린 모델의 경우, 교정 과정을 몇 차례 거친 후 오차가 4% 이상에서 0.5% 미만으로 떨어졌습니다. 연구진은 첫 번째 추측이 종종 근접하긴 했지만, 미세한 디테일을 포착하는 데 있어서는 후속 교정 작업이 필수적이라는 것을 발견했습니다.

이 연구는 이 방법이 테스트된 특정 유형의 비선형성에 대해서는 실행 가능하다는 결론을 내렸지만, 동시에 현재 기술의 한계도 강조했습니다. 연구진은 자신들이 선택한 단순한 1차원 사례에서는 이 방법이 잘 작동하지만, 다차원 및 복잡한 재료 거동을 포함하는 더 복잡한 실제 시나리오에는 더 많은 큐비트와 더 정교한 보조 변수 처리가 필요할 것이라고 언급했습니다. 그들은 이 방법으로 모든 재료를 양자 컴퓨터로 시뮬레이션하는 문제를 해결했다고 주장하거나, 이 방법이 즉시 산업 현장에 적용될 준비가 되었다고 제안한 것이 아닙니다. 대신, 벌칙 제약과 반복적 정교화의 조합이 이전에는 접근할 수 없었던 어려운 공학 문제들에 대해 양자 알고리즘을 유용하게 만들 수 있다는 명확하고 단계적인 시연을 제시했습니다. 이 연구는 양자 컴퓨팅이 언젠가 엔지니어들이 사용하는 재료의 복잡한 물리학을 이해함으로써 더 안전하고 효율적인 구조물을 설계하는 데 어떻게 도움을 줄 수 있는지 보여주는 중요한 진전입니다.

기술 요약: 비선형 구성 모델을 갖는 초탄성 문제를 위한 변분 양자 알고리즘

문제 정의
본 논문은 근시점(near-term) 양자 하드웨어를 사용하여 비선형 탄성 정적 경계값 문제를 해결하는 과제를 다룹니다. 최근의 발전들이 미분 방정식을 위한 양자 알고리즘에 동기를 부여해 왔으나, 두 가지 주요 장애물이 지속되고 있습니다. 하나는 강력한 이론적 보증을 필요로 하는 많은 알고리즘이 결함 허용(fault-tolerant) 하드웨어를 요구한다는 점이며, 다른 하나는 양자 상태 진화의 선형 구조가 비선형 미분 방정식을 처리하는 것을 어렵게 만든다는 점입니다. 구체적으로, 본 연구는 Ogden 및 Mooney–Rivlin 공식과 같이 유리수 거듭제곱(예: 분수 및 역수 거듭제곱)을 포함하는 구성 모델을 따르는 초탄성 재료를 대상으로 합니다. 이러한 비선형성은 이전에 선형 또는 로그 비선형성을 위해 제안된 얕은 회로(shallow-circuit)의 변분 양자 알고리즘(VQA) 프레임워크 내에서 직접 표현하기 어렵습니다.

방법론
저자들은 유리수 거듭제곱을 포함하는 구성 비선형성을 처리할 수 있도록 이전에 제안된 1D 초탄성 VQA 프레임워크를 확장합니다. 방법론은 세 가지 핵심 구성 요소로 이루어집니다:

비용 함수로의 재구성: 탄성 정적 문제를 위치 에너지 범함수 $E[u]$ 를 이산화함으로써 재구성하고, 변위 자유도를 양자 상태 $|u\rangle = \theta_0 \hat{V}(\theta)|0\rangle$ 로 인코딩합니다. 여기서 $\hat{V}(\theta)$ 는 매개변수화된 양자 회로(ansatz)이며, 위치 에너지는 최소화해야 할 비용 함수 $C(\theta_0, \theta)$ 역할을 합니다.
페널티 방법을 통한 다항식 근사: 현재의 양자 프리미티브(특히 양자 비선형 처리 장치 또는 QNPU)는 임의의 유리수 거듭제곱이 아닌 다항식 연산을 효율적으로 평가할 수 있으므로, 변형 에너지 밀도 $W(\lambda)$ 를 거듭제곱 항들의 합으로 근사합니다. 분수 또는 역수 거듭제곱(예: $p \in \mathbb{Q}$ 인 $\lambda^{p}$ )을 제거하기 위해, $y = \lambda^{p}$ 와 같은 제약 조건을 강제하는 보조 변수 $y$ 를 도입합니다. 이러한 제약 조건은 다음과 같은 페널티 항을 통해 전체 위치 에너지에 통합됩니다:
$\tilde{E}[u, y] = E[u, y] + \frac{\beta}{2} \int_{\Omega_0} (y^b - \lambda^a)^2 dX$
이를 통해 문제는 VQA와 호환되는 다항식 비용 함수로 변환됩니다.
반복 수정 전략: 페널티 공식이 근사 오류(다항식 근사와 페널티 파라미터 $\beta$ 로 인한)를 유발한다는 점을 인식하여, 저자들은 반복 수정 전략을 도입합니다. 근사해 $\tilde{u}$ 로부터부터 수정 필드 $\delta u$ 를 찾아내어 $u = \tilde{u} + \delta u$ 가 되도록 합니다. 현재 근사값에 대해 변형 에너지 밀도의 테일러 전개를 수행함으로써, 수정 항을 위한 새로운 비용 함수를 구축합니다. 이 수정 문제는 일련의 VQA를 통해 순차적으로 해결되어, 솔루션을 반복적으로 업데이트( $u^{(k+1)} = u^{(k)} + \delta u^{(k)}$ )하여 정확도를 향상시킵니다.

주요 기여

VQA 프레임워크의 확장: 본 연구는 선형 및 로그 비선형성을 넘어 유리수 거듭제곱을 포함하는 더 넓은 범위의 구성 모델(특히서 오그던(Ogden) 및 무니-리블린(Mooney–Rivlin) 모델의 적용 가능성을 입증)에 대한 VQA 프레임워크를 확장했습니다.
보조 변수 및 페널티 공식화: 비다항식 구성 항을 처리하기 위해 보조 변수와 페널티 제약을 도입하는 방법을 제안하여, 에너지 범함수를 가용한 양자 알고리즘 프리미티브와 호환되도록 만들었습니다.
반복 수정 알고리즘: 페널티 기반 근사에 내재된 오류를 완화하기 위해, 일련의 VQA 해결 과정을 통해 솔루션을 체계적으로 정밀하게 만드는 새로운 반복 수정 전략을 도입했습니다.

결과
제안된 접근 방식은 Qiskit을 이용한 상태 벡터 시뮬레이션을 통해 1차원 비압축성 오그던 모델( $\alpha=2$ 및 $\alpha=3$ 파라미터)과 무니-리블린 모델에 대해 검증되었습니다.

정확도: 반복 수정 절차는 테스트된 모든 모델에 대해 페널티 기반 솔루션의 정확도를 유의미하게 향상시켰습니다. 해석적 해에 대한 상대적 $L_2$ 오차는 5번의 반복 후에 상당히 감소했습니다. 예를 들어, 무니-리블린 모델은 오차가 4.27%(페널티만 사용 시)에서 0.49%(반복 적용 시)로 감소했습니다.
수렴성: 솔로션이 해석적 응답에 가까워질수록 최적화가 점점 더 어려워지는 현상이 관찰되었으며, 저자들은 상당한 정확도 이득이 이미 달성된 시점인 5회 반복 후에 프로세스를 종료했습니다.
구현: 시뮬레이션에는 3-큐비트 시스템( $N_q=8$ )과 회로 깊이 $d=4$ 가 사용되었습니다.

의의 및 주장
본 논문은 VQA가 이전에는 근시점 양자 알고리즘으로 다루기 어려웠던 유리수 거듭제곱 법칙을 포함하는 비선형 초탄성 문제를 해결하는 데 적응될 수 있음을 보여주는 "개념 증명(proof of concept)"으로서 이 작업을 위치시킵니다. 저자들은 페널티 기반 근사와 반복 수정의 결 조합이 결함 허용 하드웨어를 요구하지 않고도 솔루션 정확도를 개선할 수 있는 실행 가능한 경로를 제공한다고 주장합니다.

저자들은 현재 예시들이 단 하나의 보조 변수만을 필요로 한다는 점을 언급하며, 범위에 대해 겸허한 태도를 유지합니다. 즉, 여러 개의 분수 거듭제곱 항을 포함하는 더 복잡한 구성 모델은 추가적인 변수와 더 비용이 많이 드는 양자 구현을 필요로 할 것입니다. 또한, 본 연구는 현재 1차원 문제로 제한되어 있으며, 저자들은 유한 요소법(finite element method)의 능력을 양자 환경에서 완전히 활용하기 위해 고차원 문제로의 확장과 비정형 격도(unstructured grids)의 처리가 중요한 향로 방향임을 밝히고 있습니다. 본 논문은 이러한 특정 1D 예시에 대해 지수적 이득을 주장하지 않으나, 잠재적 이득에 관한 이전 연구와 유사한 조건에 따라 복잡도 분석이 따른다고 명시했습니다.

유사한 논문