수소결합

계산화학 분자동역학 마지막 수정: 2026-08-16 04:47:52

1. 개요[편집]

수소결합 Hydrogen Bond
분야	계산화학 × 분자동역학 × 물리화학
표기	D–H···A (주개–수소···받개)
정의	IUPAC 2011 권고 (Arunan 외)
에너지 범위	약 1~40 kcal/mol
대표 사례	물, 얼음, DNA 염기쌍, 단백질 이차구조
계산상 지위	힘장에 명시적 항이 없다

공유결합보다 스무 배 약하고 반데르발스보다 열 배 강한, 딱 그 사이의 힘. 그런데 생명은 정확히 그 창문 안에서 굴러간다.

수소결합(hydrogen bond)은 전기음성도가 큰 원자 D에 공유결합된 수소가, 다른 원자 또는 원자단 A의 전자 밀도와 방향성 있는 인력 상호작용을 형성하는 현상이며 관례적으로 D–H···A로 적는다. D를 주개(donor), A를 받개(acceptor)라 부른다. 상호작용 에너지가 대략 1~40 kcal/mol 구간에 퍼져 있어, 열 요동( $k_BT \approx 0.6$ kcal/mol, 상온)으로 끊었다 붙였다 할 수 있으면서도 구조를 유지할 만큼은 강하다. 끊어질 수 있는데 안 끊어지는 이 절묘한 세기가 물의 이상성부터 DNA 이중나선의 지퍼 구조까지를 모두 설명한다.

계산 쪽에서 수소결합이 특별한 이유는 따로 있다. 이 결합은 어떤 시뮬레이션 코드에도 “수소결합 항”으로 존재하지 않는다. 힘장에서는 부분전하 사이의 쿨롱 인력과 레너드-존스 퍼텐셜 반발이 우연히 균형을 이룬 결과로 나타나고, 양자화학에서는 그냥 전자 구조 계산의 부산물이다. 그래서 “수소결합이 몇 개인가”라는 질문에 답하려면 먼저 판정 기준을 정의해야 하는데, 그 기준이 하나로 통일되어 있지 않다.

2. IUPAC 2011 정의와 기하[편집]

2011년 IUPAC 권고(Arunan 외)는 전기음성도 목록으로 후보 원자를 못 박던 낡은 정의를 버리고, **“결합 형성의 증거가 있는가”**를 기준으로 삼았다. 요지는 이렇다.

D–H 결합에서 D가 H보다 전기음성이어야 하고, 상호작용의 증거가 실험적·이론적으로 제시되어야 한다.
인정되는 증거로는 결합 길이 단축, D–H 신축 진동의 적색이동과 세기 증가, NMR 화학이동 및 스핀-스핀 결합 상수, 전자밀도 위상 해석(AIM)의 결합 임계점 등이 열거된다.
이 정의 아래에서 C–H···O, D–H··· $\pi$ , 심지어 M–H···H–X 형태의 이수소 결합까지 수소결합의 가족으로 들어온다. 반대로 단순히 거리가 가깝다는 것만으로는 자격이 없다.

기하학적으로는 세 개의 숫자가 쓰인다.

양	전형적 값	비고
H···A 거리	vdW 반지름 합보다 짧음, 대개 1.5~2.5 Å	짧을수록 강함
D···A 거리	물에서 2.7~3.0 Å	결정학에서 흔히 보고되는 값
D–H···A 각도	110° 이상, 강할수록 180°에 접근	방향성의 근원

각도 조건이 수소결합을 반데르발스 인력과 구별 짓는 핵심이다. 반데르발스 힘은 방향을 가리지 않지만 수소결합은 선형을 선호하고, 그래서 물이 사면체 그물을 짜고 단백질이 $\alpha$ -나선을 감는다. 에너지가 아니라 방향성이 구조를 만든다.

3. 성분 분해 — 무엇이 붙잡고 있는가[편집]

“수소결합은 정전기 상호작용인가 공유결합인가”는 이 바닥의 100년 묵은 논쟁이고, 대칭적응 섭동론(SAPT)이나 에너지 분해 분석(EDA)이 그 답을 성분별 숫자로 제시한다. 물 이량체(CCSD(T)/CBS 기준 결합 에너지 $D_e \approx 5.0$ kcal/mol, 영점 보정 후 $D_0 \approx 3.2$ kcal/mol)의 전형적인 분해는 대략 이렇다.

성분	대략적 크기 (kcal/mol)	성격
정전기	−8 내외	쌍극자·다중극 인력, 가장 큼
교환-반발	+7 내외	파울리 반발, 정전기를 거의 상쇄
유도(분극 + 전하이동)	−2 내외	비가법적, 협동성의 원천
분산	−2 내외	항상 인력, 무시하면 안 됨
합계	−5 내외	큰 수들의 작은 차이

읽는 법이 중요하다. 총합 5 kcal/mol은 8과 7이라는 큰 값의 차이다. 그러니 어느 한 성분을 10%만 틀려도 총합은 20% 틀린다. 수소결합 계산이 유난히 까다로운 근본 이유가 이 상쇄 구조에 있다.

성분 이름에도 함정이 있다. 전하이동을 분극에서 분리하는 것은 기저 의존적이라, 기저를 키우면 전하이동 성분이 줄어들고 분극이 커지는 일이 흔하다. “수소결합의 X%는 전하이동”이라는 문장을 볼 때는 어떤 분해 스킴을 썼는지부터 확인해야 한다. 마찬가지로 분산 성분이 20%나 된다는 사실은, 밀도범함수이론으로 수소결합을 다룰 때 분산 보정을 반드시 켜야 하는 이유가 된다.

4. 협동성 — 덧셈이 안 되는 힘[편집]

수소결합의 가장 성가신 성질은 비가법성이다. 물 분자 세 개를 사슬로 이으면 총 에너지가 이량체 두 개의 합보다 눈에 띄게 크다. 앞의 결합이 만든 분극이 다음 결합을 강화하기 때문이다. 액체 물의 응집 에너지에서 3체 이상의 다체 항이 차지하는 비중은 대략 15~25%로 알려져 있다.

이 협동성은 두 가지 실무적 결과를 낳는다. 첫째, 고정전하 힘장은 원리적으로 이걸 담을 수 없다. 부분전하는 환경에 반응하지 않기 때문이다. TIP3P·SPC/E 같은 고정전하 물 모형은 그 다체 효과를 이체 파라미터에 미리 녹여 넣는 방식으로 우회하며, 그래서 기체상 이량체 에너지는 틀리는 대신 액체 물성은 맞는다. 이런 모형을 기체상 벤치마크로 채점하면 억울한 점수가 나온다. 둘째, 분극성 힘장(AMOEBA, Drude 진동자)이 개발되는 가장 큰 명분이 바로 이 항목이다.

5. 물이 이상한 이유[편집]

물의 기이한 거동은 거의 전부 “분자 하나가 수소결합을 최대 네 개(주개 2 + 받개 2) 만들 수 있고, 그 방향이 사면체로 정해져 있다”는 단 한 문장에서 나온다.

4 °C 밀도 최대. 정확히는 3.98 °C. 얼음 Ih의 사면체 그물은 성기게 뚫려 있어 액체보다 밀도가 낮고(약 9% 팽창), 녹으면 그물이 부분적으로 무너지며 밀도가 오른다. 온도를 더 올리면 열팽창이 이겨서 다시 내려간다. 두 효과가 교차하는 지점이 4 °C 근처다. 호수가 위에서부터 어는 덕에 물고기가 겨울을 난다는 그 이야기.
비정상적으로 높은 끓는점. 같은 족의 H $_2$ S는 −60 °C에서 끓는데 물은 100 °C다. 분자량으로 외삽하면 물은 상온에서 기체여야 한다.
큰 비열과 기화열. 몰 비열 약 75 J/(mol·K), 기화열 약 40.7 kJ/mol. 기화열이 단일 수소결합 에너지(약 5 kcal/mol = 21 kJ/mol)의 두 배쯤인 것은 우연이 아니다 — 분자 하나가 평균 두 개 몫의 결합을 끊고 나가야 하기 때문.
얼음의 잔류 엔트로피. 얼음의 산소 골격은 규칙적이지만 수소 위치는 베르날-파울러 얼음 규칙(산소당 가까운 수소 2개)을 지키는 한 무질서할 수 있다. 파울링이 이 배열 수를 세어 예측한 $R\ln(3/2) \approx 3.4$ J/(mol·K)이 열량 측정과 맞아떨어진 것은 통계역학의 고전적 승리다.¹

6. 대칭 수소결합과 저장벽[편집]

보통의 수소결합에서 양성자는 명백히 D 쪽에 붙어 있고, A 쪽으로 넘어가려면 장벽을 넘어야 한다. 그런데 D와 A가 화학적으로 동등하고 아주 가까워지면 이 이중우물이 얕아지다 못해 하나로 합쳐진다.

이플루오린화 이온 $\mathrm{HF_2^-}$ 가 교과서적 사례다. F···F 거리가 2.28 Å로 극단적으로 짧고 양성자가 정확히 중앙에 앉으며, 결합 세기는 약 39 kcal/mol — 이미 약한 공유결합의 영역이다. 수소결합 에너지 스펙트럼의 최상단이 여기다.
얼음 X. 얼음을 60 GPa 이상으로 압축하면 O···O 거리가 줄어 양성자가 두 산소의 중간에 자리 잡는다. “수소결합과 공유결합의 구분이 사라진 얼음”으로, 고압 물리의 유명한 상이다.
저장벽 수소결합(LBHB). 효소 활성부위에서 이런 구조가 촉매 능력을 설명한다는 가설이 1990년대에 크게 유행했다. 다만 단백질처럼 유전율이 높고 불균일한 환경에서 정말 형성되는지, 형성되더라도 촉매 기여가 유의미한지는 여전히 논쟁 중이다. 계산으로 판정하려면 양성자의 핵 양자 효과까지 다뤄야 해서 난도가 급상승한다.²

7. 시뮬레이션에서 다루기[편집]

7.1. 판정 기준부터 정하라[편집]

분자동역학 궤적에서 “수소결합 개수”를 뽑으려면 정의가 필요하고, 실무에서 쓰이는 것은 크게 둘이다.

기하 기준. 가장 널리 쓰이는 루자르-챈들러(Luzar–Chandler) 기준은 O···O 거리가 3.5 Å 이하(동경분포함수 $g_{\mathrm{OO}}$ 의 첫 최소점)이고 H–O···O 각도가 30° 이내일 때 결합으로 센다. 상온 액체 물에서 분자당 평균 3.5개 안팎이 나온다.
에너지 기준. 쌍 상호작용 에너지가 임계값(예: −10 kJ/mol)보다 낮으면 결합으로 센다.

두 기준이 주는 개수는 다르다. 그러니 “물 분자당 수소결합 3.6개”라는 숫자를 다른 논문과 비교하려면 기준이 같은지부터 확인해야 한다. 임계값 근처에서 깜빡이는 결합을 어떻게 셀 것인가(수명 정의, 연속 vs 간헐 상관함수)도 같은 부류의 문제다.

7.2. 힘장에서의 지위[편집]

앞서 말했듯 현대 생체분자 힘장에 수소결합 항은 없다. 1980년대 AMBER·CHARMM 초기 판본에는 방향성을 강제하는 10-12 퍼텐셜 항이 따로 있었으나, 부분전하와 반발 항만으로도 충분히 재현된다는 결론이 나면서 1990년대에 제거됐다. 그래서 오늘날 힘장의 수소결합 품질은 부분전하 배치와 물 모형 선택이 사실상 전부를 결정한다.

7.3. 양자화학으로 다룰 때의 지뢰밭[편집]

기저중첩오차. 5 kcal/mol짜리 결합을 다루는데 작은 기저의 BSSE가 1 kcal/mol을 넘을 수 있다. 카운터포이즈 보정과 완전기저 극한 외삽은 선택이 아니라 필수다.
분산 누락. 순수 GGA는 분산을 못 담아 이량체를 과소결합시키고, 반대로 어떤 범함수는 자체 상호작용 오차로 양성자 전달 장벽을 과소평가한다. 분산 보정(D3/D4류)이 표준 처방이 된 배경.
핵 양자 효과. 수소는 가볍다. 영점 진동과 터널링이 결합 길이와 세기를 눈에 띄게 바꾸며, D $_2$ O의 녹는점이 3.82 °C로 올라가는 동위원소 효과가 그 직접적 증거다. 고전 분자동역학은 이걸 통째로 놓치므로, 정밀 계산에는 경로적분 분자동역학이 동원된다.
진동수 검증. 수소결합이 형성되면 D–H 신축 진동수가 수십~수백 cm $^{-1}$ 적색이동한다. 계산이 이 이동량을 재현하는지 보는 것이 가장 값싼 검증 및 확인 절차다. 참고로 일부 C–H···O 계에서는 오히려 청색이동이 관측되는데, 이는 수소결합이 아니라는 뜻이 아니라 진동 결합의 다른 메커니즘이 작동한다는 뜻이다.³

8. 관련 문서[편집]

파울링의 계산은 “산소 하나에 붙을 수 있는 수소 배열 6가지 중 이웃과 충돌하지 않는 비율”을 세는 대단히 소박한 조합론인데, 실험값과 소수점 첫째 자리까지 맞는다. 물리학에서 냅킨 계산이 이렇게까지 이긴 사례도 흔치 않다. ↩
그래서 LBHB 논문은 대개 “형성된다”파와 “그냥 짧은 수소결합일 뿐”파로 갈린다. 결정 구조의 전자밀도만으로는 수소 위치를 확정하기 어렵다는 점이 논쟁을 20년째 살려두고 있다. 중성자 회절이 결정적 증거를 줄 수 있는데, 그러려면 큰 결정을 키워야 한다는 현실의 벽이 있다. ↩
청색이동 수소결합은 발견 당시 “수소결합의 정의를 위협하는 반례”로 취급됐지만, 결국 IUPAC 2011 정의가 신축 이동의 부호가 아니라 결합 형성의 증거 전반을 보라고 정리하면서 가족 안으로 편입됐다. 정의를 넓혀 반례를 흡수한 셈인데, 이런 일은 과학사에서 꽤 자주 일어난다. ↩