Numerical analysis for a unified 2 factor model of structural and reduced form types for corporate bonds with fixed discrete coupon

기술 요약: 고정 이산 쿠폰을 가진 기업 채권의 통합 2-요인 모델에 대한 수치 해석

문제 정의
본 논문은 고정 이산 쿠폰을 지급하는 부도 가능 기업 채권의 가격 결정을 다룬다. 제로 쿠폰 채권(Merton, 1977)과 쿠폰 채권을 위한 구조적 모델(Geske, 1977; Agliardi, 2011)에 대해서는 상당한 이론적 연구가 존재하지만, 고정 이산 쿠폰에 특화되어 구조적 접근법과 축약형(reduced-form) 접근법을 결합한 통합 모델에 대한 해석적 해법에는 공백이 존재한다.

저자들은 통합된 1-요인 모델(고정 이산 쿠폰용)은 해석적 가격 공식이 존재하고, 할인된(discounted) 이산 쿠폰을 위한 통합 2-요인 모델은 해결된 반면, 고정된(fixed) 이산 쿠폰을 위한 통합 2-요인 모델은 해석적 해법이 없는 것으로 보인다고 언급한다. 이는 편미분 방정식(PDE)의 복잡성, 즉 세 개의 변수(기업 가치, 단기 이자율, 시간)와 혼합 도함수를 포함하는 포물선형 방정식 때문인 것으로 보인다. 또한, 모델은 두 가지 유형의 부도를 모두 고려해야 한다:

1. 기대 부도(Expected Default): 기업의 자본이 부채와 쿠폰을 충당하기에 불충분할 때, 오직 이산적인 쿠폰 지급일에만 발생한다.
2. 예상치 못한 부도(Unexpected Default): 일정한 부도 강도(hazard rate)에 의해 지배되며, 언제든 발생할 수 있다.

방법론
이 문제를 해결하기 위해 저자들은 다음을 포함하는 리스크 중립 마팅게일 측도 하의 통합 2-요인 모델을 제안한다:

• 바시체크(Vasicek) 모델: 확률적 단기 이자율($r$)을 위해 사용된다.
• 기하 브라운 운동(Geometric Brownian Motion): 일정한 배당 수익률을 가진 기업 가치($V$)를 위해 사용된다.
• 상관관계: 단기 이자율과 기업 가치를 구동하는 위너 프로세스(Wiener processes) 사이의 상관계수($\rho$).

가격 결정 문제는 PDE의 초기값 문제 집합(식 2.4)으로 정식화된다. 해법은 쿠폰 날짜 사이의 하위 구간들에 대해 정의된 PDE 시퀀스를 포함하며, 경계 조건은 기업 가치가 부채를 충당하는지 여부에 따라 '채권 가격 + 다음 쿠폰'과 '회수 가치' 사이를 전환된다.

일반적인 경우($i < N-1$)에 대해 해석적 공식이 구하기 어렵다고 판단되므로, 저자들은 수치 해석을 위해 **명시적 유한 차분 스킴(explicit finite difference scheme)**을 채택한다. 주요 방법론적 단계는 다음과 같다:

1. 변수 변환: 도메인 $V > 0$를 $x \in \mathbb{R}$로 처리하기 위해 기업 가치 변수를 $x = \ln V$로 변환한다.
2. 이산화: 시간 도함수에는 전방 차분(forward differences)을, 공간 도함수에는 중앙 차분(central differences)을 사용한다. 혼합 도함수($\partial^2 B / \partial x \partial r$)의 경우, 상관계수 $\rho$의 부호에 따라 중앙, 전방, 또는 혼합 전방-후방 차분 스킴을 탐색한다.
3. 안정성 분석: 명시적 스킴에 대한 엄격한 안정성 조건 유도. 저자들은 혼합 도함수로 인해 표준 명시적 스킴이 일반적으로 불안정하다는 것을 보여준다. 이들은 안정성을 보장하기 위해 확산 계수 및 상관계수 $\rho$에 대한 시간 단계($\Delta t$)와 공간 단계($\Delta x, \Delta r$)의 특정 제약 조건을 유도한다.

주요 기여

1. 통합 모델 정식화: 고정 이산 쿠폰 채권을 위해 Agliardi의 구조적 모델을 통합 프레임워크로 일반화한 통합 2-요인 모델을 정식화하였다.
2. 안정성 조건: 이 특정 3-변수 포물선형 PDE(혼합 도함수 포함)에 적용되는 명시적 유한 차분 스킴을 가능하게 하는 안정성 조건을 유도하였다.
   • $\rho = 0$인 경우, 안정성을 위해서는 $\mu_x + \mu_r < 1$ 및 드리프트 항에 대한 특정 경계가 필요하다.
   • $0 < \rho < 1/2$인 경우, 저자들은 안정성을 위해 $0 < \rho < 1/2$와 공간 단계 크기 간의 특정 비율($\sqrt{\mu_r/\mu_x}$)이 필요함을 보여준다.
   • $-1 < \rho < 0$인 경우, 혼합 도함수를 위해 혼합 전방-후방 차분을 사용하는 유사한 제약 조건이 유도된다.
3. 수치적 구현: 해석적 공식이 존재하지 않는 경우를 대비하여 채권 가격, 신용 스프레드, 듀레이션을 계산하기 위한 계산 프레임워크를 제공한다.

결과
$\rho = 0$인 유도된 명시적 스킴을 사용하여, 저자들은 특정 파라미터 세트(예: $N=2$ 쿠폰, 바시체크 파라미터, 특정 부도 강도 $\lambda_0, \lambda_1$)를 바탕으로 수치 결과를 제시한다. 분석 결과는 다음과 같다:

• 채권 가격:
  • 단기 이자율($r$)이 증가함에 따라 가격이 감소한다.
  • 기업 가치($V$)가 증가함에 따라 가격이 증가하며, $V$에 대해 볼록성(convexity)을 보인다.
  • 쿠폰 지급일에 점프(불연속성)를 보이는데, 이는 쿠폰 지급과 부도 위험 조건의 재설정을 반영한다.
• 신용 스프레드 (CS):
  • 기업 가치가 증가함에 따라 CS가 감소한다.
  • 기업 가치가 낮을 경우 만기가 가까워질수록 CS가 급격히 증가한다.
  • 단기 이자율이 증가함에 따라 CS가 감소한다(무위험 구성 요소에 대한 할인 효과 때문).
  • 더 높은 부도 강도($\lambda$)와 함께 CS가 증가한다.
  • 쿠폰 날짜에 점프를 보이며, 점프의 크기는 쿠폰 규모에 따라 커진다.
• 듀레이션:
  • 기업 가치와 단기 이자율이 증가함에 따라 듀레이션이 증가한다.
  • 쿠폰 지급 후 듀레이션이 점프를 보이는데, 이는 현금 흐름에 대한 가중 평균 시간이라는 듀레이션의 정의와 일치한다.

의의 및 주장
본 논문은 Agliardi의 구조적 모델을 고정 이산 쿠폰 채권을 위한 통합 프레임워크로 일반화하는 데 필요한 수치적 도구를 제공한다는 점을 겸허히 주장한다. 저자들은 이 특정 통합 2-요인 모델에 대한 해석적 공식을 얻는 것이 "어렵고" 일반적인 경우 불가능할 것으로 보인다고 밝힌다. 따라서 본 연구의 의의는 다음과 같다:

1. 이 복잡한 PDE 시스템을 적용하는 데 있어 명시적 유한 차분 스킴을 실행 가능하게 만드는 안정성 조건을 확립한 것.
2. 수치적 결과가 실제 금융 상황과 호환됨을 입증한 것. 이는 다양한 시장 조건(이자율, 기업 가치, 부도 강도) 하에서 채권 가격, 신용 스프레드 및 듀레이션의 직관적인 움직임을 통해 증명된다.

본 연구는 해석적 해법이 존재하지 않는 현실적인 채권 구조(고정 이서 쿠폰)를 위한 이론적 구조적/축약형 통합과 실무적 가격 결정 사이의 가교 역할을 한다.