Motivated by recent developments in the calibration of stochastic volatility models (SVMs for short), we study continuous-time formulations of martingale optimal transport and martingale Schrödinger bridge problems. We establish duality formulas and also provide alternative proofs, via different techniques, of duality results previously established in the mathematical finance literature. Applications include calibration of SVMs to SPX options, as well as joint calibration to both SPX and VIX options.
논문 ID : 2510.10860제목 : Martingale Optimal Transport and Martingale Schrödinger Bridges for Calibration of Stochastic Volatility Models저자 : Antonios Zitridis분류 : math.OC (최적화 및 제어), math.PR (확률론)발표 시간 : 2025년 10월 12일 (arXiv 사전인쇄본)논문 링크 : https://arxiv.org/abs/2510.10860 본 논문은 확률적 변동성 모델(SVM) 보정 분야의 최근 발전에 영감을 받아, 마팅게일 최적 수송과 마팅게일 슈뢰딩거 다리 문제의 연속시간 표현을 연구한다. 본 논문은 쌍대 공식을 수립하고, 다양한 기법을 통해 금융수학 문헌에서 이전에 수립된 쌍대 결과에 대한 대체 증명을 제공한다. 응용 분야에는 SVM을 SPX 옵션에 보정하는 것과 SPX 및 VIX 옵션에 동시에 보정하는 것이 포함된다.
본 논문이 해결하고자 하는 핵심 문제는 확률적 변동성 모델의 보정 문제이다. Black-Scholes 모델 도입 이후, 학계는 시장 동학을 정확하게 포착할 수 있는 복잡한 변동성 모델 개발에 노력해왔다. 체계적 오류를 수정하고 정확한 예측을 보장하기 위해, 이러한 모델들을 바닐라 옵션(예: S&P 500 옵션 및 VIX 옵션)에 적절히 조정해야 한다.
실제 응용 필요성 : 금융시장의 옵션 가격 책정 정확성은 거래 전략 및 위험 관리에 직접적인 영향을 미친다이론적 완성 : 금융 이론의 일관성을 보장하기 위해 무차익 모델 보정 방법이 필요하다계산 효율성 : 기존 방법은 수치 구현에서 어려움이 있으며, 수치 근사를 지원하기 위한 더 나은 쌍대 공식이 필요하다연속시간 마팅게일 최적 수송 문제와 그 이산 버전 간의 연결이 명확하지 않다 SPX 및 VIX 옵션을 동시에 보정하는 문제를 다루기 위한 통일된 이론 틀이 부족하다 기존 증명 기법은 주로 Fenchel-Rockafellar 정리에 의존하며, 대체 방법이 부족하다 저자는 평균장 게임 이론의 기법, 특히 Von-Neumann 극소극대 정리를 채택하여 마팅게일 최적 수송과 마팅게일 슈뢰딩거 다리 문제에 대한 새로운 이론적 분석 방법을 제공한다.
연속시간 마팅게일 최적 수송 문제의 쌍대 공식을 수립 하고, Von-Neumann 극소극대 정리에 기반한 새로운 증명 방법을 제공한다마팅게일 슈뢰딩거 다리 문제에 대한 새로운 쌍대 공식 증명을 제공 하며, 기존 문헌과 다른 기법 경로를 사용한다SPX-VIX 연합 보정 문제를 해결 하고, VIX 제약이 있는 마팅게일 슈뢰딩거 다리에 대한 엄밀한 쌍대 공식을 제공한다점성해 이론의 응용을 제공 하여 Dirac 델타 항을 포함하는 Hamilton-Jacobi 방정식을 다룬다최적해의 명시적 표현을 제시 하며, 최적 확률 측도의 Radon-Nikodym 도함수를 포함한다주어진 확률 측도 μ 0 , μ 1 , μ 2 ∈ P ( R ) \mu_0, \mu_1, \mu_2 \in P(\mathbb{R}) μ 0 , μ 1 , μ 2 ∈ P ( R ) 가 볼록 순서 관계 μ 0 ≤ c μ 1 ≤ c μ 2 \mu_0 \leq_c \mu_1 \leq_c \mu_2 μ 0 ≤ c μ 1 ≤ c μ 2 를 만족하고, 시간 매개변수 0 < T 1 < T 2 0 < T_1 < T_2 0 < T 1 < T 2 및 T 0 ∈ [ 0 , T 2 ] T_0 \in [0, T_2] T 0 ∈ [ 0 , T 2 ] 일 때, 마팅게일 최적 수송 문제는 다음과 같이 정의된다:
V ( T 0 , μ 0 , μ 1 , μ 2 ) : = inf P ∈ P m { E P [ ∫ T 0 T 2 L ( σ t 2 ) d t ] } V(T_0, \mu_0, \mu_1, \mu_2) := \inf_{P \in P_m} \left\{ E^P \left[ \int_{T_0}^{T_2} L(\sigma_t^2) dt \right] \right\} V ( T 0 , μ 0 , μ 1 , μ 2 ) := inf P ∈ P m { E P [ ∫ T 0 T 2 L ( σ t 2 ) d t ] }
여기서 P m P_m P m 은 마팅게일 조건과 주변 분포 제약을 만족하는 확률 측도의 집합이다.
참조 측도 P 0 P_0 P 0 와 제약 집합 C t 0 ( μ 1 , μ 2 ) C_{t_0}(\mu_1, \mu_2) C t 0 ( μ 1 , μ 2 ) 가 주어질 때, 마팅게일 슈뢰딩거 다리 문제는 다음과 같다:
D P 0 : = inf P ∈ C t 0 ( μ 1 , μ 2 ) H ( P ∣ P 0 ) D_{P_0} := \inf_{P \in C_{t_0}(\mu_1, \mu_2)} H(P|P_0) D P 0 := inf P ∈ C t 0 ( μ 1 , μ 2 ) H ( P ∣ P 0 )
여기서 H ( P ∣ P 0 ) H(P|P_0) H ( P ∣ P 0 ) 는 상대 엔트로피이다.
연속시간 확률적 변동성 모델을 고려한다:
{ d X t = σ ( X t , Y t ) d W t d Y t = b ( X t , Y t ) d t + τ 1 ( X t , Y t ) d W t + τ 2 ( X t , Y t ) d W t ⊥ \begin{cases}
dX_t = \sigma(X_t, Y_t) dW_t \\
dY_t = b(X_t, Y_t) dt + \tau_1(X_t, Y_t) dW_t + \tau_2(X_t, Y_t) dW_t^{\perp}
\end{cases} { d X t = σ ( X t , Y t ) d W t d Y t = b ( X t , Y t ) d t + τ 1 ( X t , Y t ) d W t + τ 2 ( X t , Y t ) d W t ⊥
여기서 W , W ⊥ W, W^{\perp} W , W ⊥ 는 독립적인 브라운 운동이고, σ ( x , y ) = x σ ~ ( y ) \sigma(x,y) = x\tilde{\sigma}(y) σ ( x , y ) = x σ ~ ( y ) 이다.
쌍대 문제의 해는 다음 Hamilton-Jacobi 방정식으로 특성화된다:
{ − ∂ t u + H ( D 2 u 2 ) = δ T 1 ( t ) u 1 u ( T 2 , x ) = u 2 ( x ) \begin{cases}
-\partial_t u + H\left(\frac{D^2u}{2}\right) = \delta_{T_1}(t)u_1 \\
u(T_2, x) = u_2(x)
\end{cases} { − ∂ t u + H ( 2 D 2 u ) = δ T 1 ( t ) u 1 u ( T 2 , x ) = u 2 ( x )
여기서 H ( a ) = sup b ≥ 0 { − a b − L ( b ) } H(a) = \sup_{b \geq 0}\{-ab - L(b)\} H ( a ) = sup b ≥ 0 { − ab − L ( b )} 는 Hamiltonian이다.
기존의 Fenchel-Rockafellar 방법과 달리, 본 논문은 Von-Neumann 극소극대 정리를 채택하여 쌍대 관계를 수립한다. 이는 다음을 검증해야 한다:
본 논문은 Dirac 델타 함수를 포함하는 Hamilton-Jacobi 방정식을 다루기 위한 점성해 이론을 발전시키며, 점프 조건을 정의한다:
u ( T 1 − , x ) − u ( T 1 + , x ) = u 1 ( x ) u(T_1^-, x) - u(T_1^+, x) = u_1(x) u ( T 1 − , x ) − u ( T 1 + , x ) = u 1 ( x )
연합 보정 문제의 경우, 원래의 VIX 등식 제약:
L P ( V P ) = μ 3 L^P(V_P) = \mu_3 L P ( V P ) = μ 3
을 볼록 순서 제약으로 약화시킨다:
L P ( V P ) ≤ c , l μ 3 L^P(V_P) \leq_{c,l} \mu_3 L P ( V P ) ≤ c , l μ 3
이는 문제를 수학적으로 더 다루기 쉽게 만든다.
본 논문은 주로 이론 연구이며, 실험 설정은 정리의 증명과 수학적 검증에 집중한다:
(A1) : L : R → R L: \mathbb{R} \to \mathbb{R} L : R → R 는 볼록 함수이고 p p p -강제적이다(A2) : Hamiltonian H H H 는 C 1 C^1 C 1 함수이다(A3) : 계수 함수의 Lipschitz 성질 및 증가 조건(A4) : 유일성을 위한 추가 정칙성 조건모방 정리를 사용하여 실행 가능성 수립 Von-Neumann 정리 적용하여 inf와 sup 교환 확률적 제어 이론을 활용하여 값 함수 특성화 점성해 이론을 통해 점프 조건 처리 가정 (A1), (A2) 하에서, 볼록 순서의 확률 측도 μ 0 ≤ c μ 1 ≤ c μ 2 \mu_0 \leq_c \mu_1 \leq_c \mu_2 μ 0 ≤ c μ 1 ≤ c μ 2 에 대해, 쌍대 공식이 성립한다:
V ( T 0 , μ 0 , μ 1 , μ 2 ) = sup u 1 , u 2 ∈ Lip { ∫ u ( T 0 , x ) d μ 0 ( x ) − ∫ u 2 ( x ) d μ 2 ( x ) − 1 [ 0 , T 1 ] ( T 0 ) ∫ u 1 ( x ) d μ 1 ( x ) } V(T_0, \mu_0, \mu_1, \mu_2) = \sup_{u_1, u_2 \in \text{Lip}} \left\{ \int u(T_0, x) d\mu_0(x) - \int u_2(x) d\mu_2(x) - \mathbf{1}_{[0,T_1]}(T_0) \int u_1(x) d\mu_1(x) \right\} V ( T 0 , μ 0 , μ 1 , μ 2 ) = sup u 1 , u 2 ∈ Lip { ∫ u ( T 0 , x ) d μ 0 ( x ) − ∫ u 2 ( x ) d μ 2 ( x ) − 1 [ 0 , T 1 ] ( T 0 ) ∫ u 1 ( x ) d μ 1 ( x ) }
여기서 u u u 는 Hamilton-Jacobi 방정식의 점성해이다.
가정 (A3) 하에서, 다음이 성립한다:
inf P ∈ C t 0 ( μ 1 , μ 2 ) H ( P ∣ P 0 ) = sup u 1 , u 2 ∈ Lip { − 1 [ 0 , T 1 ] ( t 0 ) ∫ u 1 d μ 1 − ∫ u 2 d μ 2 + u ( t 0 , X t 0 , Y t 0 ) } \inf_{P \in C_{t_0}(\mu_1, \mu_2)} H(P|P_0) = \sup_{u_1, u_2 \in \text{Lip}} \left\{ -\mathbf{1}_{[0,T_1]}(t_0) \int u_1 d\mu_1 - \int u_2 d\mu_2 + u(t_0, X_{t_0}, Y_{t_0}) \right\} inf P ∈ C t 0 ( μ 1 , μ 2 ) H ( P ∣ P 0 ) = sup u 1 , u 2 ∈ Lip { − 1 [ 0 , T 1 ] ( t 0 ) ∫ u 1 d μ 1 − ∫ u 2 d μ 2 + u ( t 0 , X t 0 , Y t 0 ) }
그리고 최적 측도의 명시적 표현을 제시한다:
d P ∗ d P 0 ∣ F T 2 = e − ∫ t 0 T 2 τ 2 ( X t , Y t ) ∂ y u ∗ ( t , X t , Y t ) d W t ⊥ − 1 2 ∫ t 0 T 2 τ 2 2 ( X t , Y t ) ∣ ∂ y u ∗ ( t , X t , Y t ) ∣ 2 d t \frac{dP^*}{dP_0}\bigg|_{F_{T_2}} = e^{-\int_{t_0}^{T_2} \tau_2(X_t, Y_t) \partial_y u^*(t, X_t, Y_t) dW_t^{\perp} - \frac{1}{2} \int_{t_0}^{T_2} \tau_2^2(X_t, Y_t) |\partial_y u^*(t, X_t, Y_t)|^2 dt} d P 0 d P ∗ F T 2 = e − ∫ t 0 T 2 τ 2 ( X t , Y t ) ∂ y u ∗ ( t , X t , Y t ) d W t ⊥ − 2 1 ∫ t 0 T 2 τ 2 2 ( X t , Y t ) ∣ ∂ y u ∗ ( t , X t , Y t ) ∣ 2 d t
약화된 VIX 제약에 대해, 쌍대 공식을 수립한다:
D P 0 ′ = sup u 1 , u 2 , u 3 { − ∫ u 1 d μ 1 − ∫ u 2 d μ 2 − ∫ u 3 d μ 3 + u ( t 0 , X t 0 , Y t 0 ) } D'_{P_0} = \sup_{u_1, u_2, u_3} \left\{ -\int u_1 d\mu_1 - \int u_2 d\mu_2 - \int u_3 d\mu_3 + u(t_0, X_{t_0}, Y_{t_0}) \right\} D P 0 ′ = sup u 1 , u 2 , u 3 { − ∫ u 1 d μ 1 − ∫ u 2 d μ 2 − ∫ u 3 d μ 3 + u ( t 0 , X t 0 , Y t 0 ) }
여기서 u 3 u_3 u 3 는 볼록이고 하한이 있는 함수이며, u u u 는 최적 제어 문제를 통해 특성화된다.
Strassen 정리와 Brenier 정리를 통해 문제의 실행 가능성을 증명한다 점성해 이론의 표준 결과를 활용하여 Hamilton-Jacobi 방정식 해의 존재 유일성을 보장한다 상태 제약이 있는 경우, DLL06 의 결과를 사용한다 추론 C.5 : 적절한 가정 하에서, 점성해 u ( t , ⋅ ) u(t, \cdot) u ( t , ⋅ ) 는 시간에 대해 균일하게 Lipschitz 연속이다.
보조정리 B.3 : τ 2 \tau_2 τ 2 가 상수일 때, 함수 Φ ( x , y ) \Phi(x,y) Φ ( x , y ) 는 연속이며, 이는 값 함수의 정칙성 수립에 중요하다.
고전 이론 : Monge-Kantorovich 문제 및 그 Benamou-Brenier 연속시간 표현마팅게일 제약 : Mikami-Thieullen과 Tan-Touzi의 반마팅게일 최적 수송연속시간 : Henry-Labordère 등의 연속시간 마팅게일 최적 수송고전 슈뢰딩거 다리 : 두 확률 측도 간의 상대 엔트로피 최소화마팅게일 제약 확장 : Labordère가 도입한 마팅게일 슈뢰딩거 다리금융 응용 : Guyon 등의 VIX 보정 응용본 논문이 이론에 중점을 두고 있지만, 쌍대 공식은 수치 방법의 기초를 제공한다:
GLW17, GLOW22: 쌍대 공식에 기반한 수치 시뮬레이션 Guy22: SPX-VIX 연합 보정의 계산 구현 통일된 이론 틀 : 마팅게일 최적 수송과 마팅게일 슈뢰딩거 다리에 대한 통일된 쌍대 이론을 제공한다새로운 증명 기법 : Von-Neumann 정리는 Fenchel-Rockafellar 정리보다 더 직접적인 증명 경로를 제공한다실제 응용 : 금융에서 중요한 연합 보정 문제를 해결한다VIX 제약 약화 : 등식 제약이 아닌 볼록 순서 제약만 다룬다정칙성 가정 : 계수의 강한 정칙성 가정이 필요하다수치 구현 : 구체적인 수치 알고리즘과 효율성 분석이 제공되지 않는다완전한 VIX 제약 : 원래의 등식 제약 처리 방법 연구수치 알고리즘 : 쌍대 공식에 기반한 효율적인 수치 방법 개발모델 확장 : 더 일반적인 확률적 변동성 모델로의 일반화이론적 혁신 : Von-Neumann 정리를 도입하여 고전 문제에 새로운 관점을 제공한다수학적 엄밀성 : 완전한 증명과 상세한 기술적 세부사항실용적 가치 : 금융의 실제 보정 문제를 해결한다명확한 표현 : 우수한 수학적 표현과 논리적 구조응용 범위 : VIX 제약의 약화는 실제 응용을 제한한다계산 복잡성 : 알고리즘의 계산 복잡도 분석이 없다실증 검증 : 실제 시장 데이터를 이용한 검증이 부족하다이론적 기여 : 최적 수송 이론에 새로운 증명 기법을 제공한다실용적 가치 : 금융공학의 모델 보정에 중요한 의미를 갖는다방법론 : Von-Neumann 방법이 다른 관련 문제에 적용될 수 있다금융기관의 옵션 가격 책정 모델 보정 위험 관리에서의 변동성 모델링 학술 연구의 이론적 분석 도구 정의 C.1 : 함수 u u u 는 방정식 (1.10)의 점성해이다. 만약 그것이 구간 ( T 1 , T 2 ] (T_1, T_2] ( T 1 , T 2 ] 와 [ 0 , T 1 ) [0, T_1) [ 0 , T 1 ) 에서 각각 해당하는 Hamilton-Jacobi 방정식을 만족하고, T 1 T_1 T 1 에서 점프 조건을 만족한다면.
보조정리 2.4 : 확률 밀도의 약 수렴 성질로, Von-Neumann 정리의 조건 수립에 중요하다.
명제 2.2, 2.3 : Fokker-Planck 방정식 해의 컴팩트성 결과로, 최적화 문제에서 수열의 수렴성을 보장한다.
본 논문은 이론 수학과 금융 응용의 교차 분야에서 중요한 기여를 하였으며, 확률적 변동성 모델의 보정을 위한 새로운 수학적 도구를 제공한다. 일부 기술적 세부사항에서 제한이 있지만, 그 이론적 가치와 잠재적 응용 전망은 주목할 만하다.