광집게

물리 계산물리 전자기학 마지막 수정: 2026-09-04 04:47:52

1. 개요[편집]

광집게 Optical Tweezers
고안	1970 방사압 가둠 · 1986 단일빔 기울기힘 트랩 — 아서 애슈킨
노벨 물리학상	2018 (애슈킨)
가두는 힘	기울기 힘 $\mathbf F_{\rm grad}\propto\nabla I$ (보존력)
방해하는 힘	산란력 $\mathbf F_{\rm scat}\propto I\hat k$ (비보존력)
전형적 강성	0.01 ~ 1 pN/nm
측정 가능 힘	0.1 ~ 100 pN
교정법	등분배 · 로렌치안 코너주파수 · 드래그 · 볼츠만

빛으로 물건을 민다는 것까지는 다들 받아들인다. 문제는 어떻게 당기느냐다.

광집게(optical tweezers)는 강하게 집속된 레이저 빔의 세기 기울기가 만드는 힘으로 마이크로미터 크기의 입자·세포·원자를 붙잡아 조작하는 장치다. 1970년 아서 애슈킨이 두 개의 마주 보는 빔으로 라텍스 구슬을 가두는 데 성공했고, 1986년 애슈킨·지에지치·비요크홀름·추가 빔 하나만으로 3차원 가둠을 구현하면서 지금의 형태가 완성됐다. 2018년 노벨 물리학상이 여기에 갔다.¹

이 장치가 물리학에서 차지하는 위치는 좀 특이하다. 광집게는 힘 발생기이자 동시에 위치 검출기라서, 나노미터·피코뉴턴 규모에서 힘과 변위를 동시에 재는 몇 안 되는 도구다. 그 덕에 분자 하나가 내는 힘을 직접 측정하는 단분자 생물물리학이 통째로 이 장치 위에 세워졌고, 반대편 극한에서는 레이저 냉각으로 식힌 중성 원자를 한 개씩 집어 배열하는 양자 시뮬레이터의 손가락 역할을 한다.

2. 두 개의 힘, 그리고 그 경쟁[편집]

빛이 물체에 주는 힘은 성격이 완전히 다른 둘로 갈린다.

산란력(scattering force). 광자의 운동량이 흡수·산란으로 전달되는 힘으로, 빔 진행 방향 $\hat k$ 를 향한다. 세기 $I$ 에 비례하는 비보존력이라 퍼텐셜이 없다. 광집게 입장에서는 입자를 초점 밖으로 밀어내는 방해물이다.
기울기 힘(gradient force). 유도된 전기 쌍극자가 비균일한 장에서 받는 힘으로, 세기가 센 쪽을 향한다( $n_p>n_m$ 일 때). $U=-\tfrac12\alpha\langle E^2\rangle$ 라는 퍼텐셜을 가진 보존력이고, 이것이 실제로 입자를 가둔다.

광집게 설계의 전부는 한 문장으로 요약된다 — 초점의 축방향 세기 기울기를 산란력이 이기지 못할 만큼 가파르게 만들어라. 가로 방향은 빔 단면의 기울기만 있으면 어차피 잡히지만, 축방향은 두 힘이 정면으로 맞선다. 그래서 광집게에는 개구수 NA 1.2 이상의 고배율 대물렌즈가 필수이고, 이 조건을 못 맞추면 입자는 그냥 빔을 타고 날아간다. 실제 평형점은 초점보다 약간 하류에 잡히는데, 거기가 기울기 힘의 되돌림 성분이 산란력과 정확히 균형을 이루는 자리다.

3. 레일리 영역 — 입자가 파장보다 훨씬 작을 때[편집]

$a\ll\lambda$ 면 입자를 점 쌍극자로 근사할 수 있다. 클라우지우스-모소티 분극률

\alpha = 4\pi\varepsilon_0 n_m^2 a^3\,\frac{m^2-1}{m^2+2}, \qquad m=\frac{n_p}{n_m}

를 쓰면 두 힘이 이렇게 나온다( $n_m$ 은 매질, $n_p$ 는 입자의 굴절률).

\mathbf F_{\rm grad}=\frac{2\pi n_m a^3}{c}\left(\frac{m^2-1}{m^2+2}\right)\nabla I, \qquad \mathbf F_{\rm scat}=\frac{n_m\sigma}{c}\,I\,\hat k

\sigma=\frac{128\pi^5a^6}{3\lambda^4}\left(\frac{m^2-1}{m^2+2}\right)^{2}

$\sigma$ 는 레일리 산란 단면적 그대로다. 여기서 지수 두 개가 모든 것을 결정한다.

F_{\rm grad}\propto a^3,\qquad F_{\rm scat}\propto a^6 \quad\Longrightarrow\quad \frac{F_{\rm grad}}{F_{\rm scat}}\propto a^{-3}

입자가 작을수록 축방향 안정성은 좋아진다. 대신 가둠 퍼텐셜의 깊이도 $a^3$ 으로 줄어드니, 열운동에 지지 않으려면 $U_{\max}\gtrsim10\,k_BT$ 를 만족해야 한다. 이 두 조건이 위아래에서 조여 오는 것이 광집게로 잡을 수 있는 입자 크기 범위를 정한다. $m<1$ 인 입자(주변보다 굴절률이 낮은 기포, 흡수체)는 기울기 힘의 부호가 뒤집혀 세기가 약한 쪽으로 밀리므로, 가운데가 어두운 도넛 모드나 광학 병 빔 같은 별도 설계가 필요하다.

4. 미 영역 — 입자가 파장보다 클 때[편집]

세포나 지름 1~5 μm 폴리스티렌 구슬은 $a\gg\lambda$ 라 쌍극자 근사가 안 통한다. 애슈킨(1992)이 정리한 표준 해석은 광선 광학이다. 빔을 광선 다발로 쪼개고, 각 광선이 구를 통과할 때 굴절되며 바뀌는 운동량을 세어 반작용을 더한다.

굴절이 기울기 힘을 만든다. 구가 초점에서 옆으로 밀려나면 세기가 강한 쪽 광선이 더 크게 꺾이므로, 알짜 반작용이 구를 빔 축으로 되돌린다. 구슬이 스스로 렌즈 역할을 해서 자기를 붙잡는 셈이다.
반사가 산란력을 만든다. 표면 반사분은 되돌려 보내지 못하고 그대로 진행 방향으로 민다.
축방향 되돌림은 크게 벌어지는 광선(즉 높은 NA)이 전담한다. 광선 광학으로 계산해 봐도 결론은 같다 — NA 가 낮으면 축방향 트랩이 성립하지 않는다.

문제는 그 사이다. $a\sim\lambda$ 인 중간 영역에서는 레일리도 광선 광학도 둘 다 틀린다. 이 영역은 미 산란 이론을 집속 빔으로 일반화한 일반화 로렌츠-미 이론이나 $T$ -행렬법으로 풀고, 형상이 구가 아니면 이산 쌍극자 근사나 FDTD로 맥스웰 방정식을 직접 푼 뒤 맥스웰 응력 텐서를 표면 적분해 힘을 얻는다. 이때 흔한 사고가 하나 있다. 광학 정리 문서에 정리된 대로, 평면파에 대해 성립하는 전방 산란 진폭 공식은 집속 빔에서 그대로 쓰면 힘이 틀린다 — 입사장과 산란장의 중첩 적분으로 일반화한 형태를 써야 한다.

5. 강성 교정 — 여기서 브라운 운동이 들어온다[편집]

평형점 근처에서 퍼텐셜을 2차까지 전개하면 $U\simeq\tfrac12 k x^2$ 이고, 광집게는 스프링 상수 $k$ (트랩 강성)를 아는 순간 힘 측정기가 된다. 변위를 재서 $F=-kx$ 로 읽으면 되기 때문이다. 그런데 $k$ 는 레이저 출력·입자 크기·굴절률·수차에 다 의존해서 이론으로 예측하기 어렵고, 실험마다 그 자리에서 교정해야 한다.² 그 교정법이 전부 브라운 운동 통계에서 나온다.

물속의 마이크로 구슬은 관성 시간이 $m/\gamma\sim100$ ns 규모라 사실상 과감쇠다(레이놀즈수가 $10^{-5}$ 급인 스토크스 유동 영역이다). 따라서 지배 방정식은 관성을 버린 랑주뱅 동역학이고, 이것은 정확히 오른슈타인-울렌벡 과정이다.

\gamma\dot x = -kx + \sqrt{2\gamma k_BT}\,\eta(t), \qquad \gamma = 6\pi\eta_{\rm vis} a

여기서 나오는 두 가지 통계가 두 가지 교정법이 된다.

1. 등분배 정리. 정상 분포가 볼츠만이므로 등분배 정리가 그대로 적용된다.

\tfrac12 k\langle x^2\rangle=\tfrac12 k_BT \quad\Longrightarrow\quad k=\frac{k_BT}{\langle x^2\rangle}

한 줄이라 매력적이지만 함정이 크다. 검출기의 nm/V 환산이 먼저 정확해야 하고, 무엇보다 $\langle x^2\rangle$ 은 모든 잡음을 그대로 삼킨다. 진동·드리프트·전자잡음이 섞이면 $\langle x^2\rangle$ 이 부풀어 $k$ 를 과소평가하고, 반대로 검출 대역이 좁아 고주파 성분을 잘라 먹으면 $k$ 를 과대평가한다. 스칼라 하나로 요약하는 방법의 숙명이다.

2. 파워 스펙트럼 밀도. 위 확률미분방정식의 파워 스펙트럼은 로렌치안이다.

S_x(f)=\frac{k_BT}{\pi^2\gamma\left(f_c^2+f^2\right)}, \qquad f_c=\frac{k}{2\pi\gamma}

$f\ll f_c$ 에서 평평하고(트랩이 붙잡고 있다), $f\gg f_c$ 에서 $f^{-2}$ 로 떨어진다(트랩이 아직 못 느낀 자유 확산). 꺾이는 코너주파수 $f_c$ 하나가 강성을 준다.

k = 2\pi\gamma f_c

이쪽이 표준이 된 이유는 명확하다. $k$ 를 얻는 데 위치 검출기의 절대 교정이 필요 없다 — $f_c$ 는 가로축(주파수)에서 읽는 값이고, 세로축 눈금이 틀려도 꺾이는 위치는 안 변한다. 게다가 잡음이 섞이면 스펙트럼에 봉우리나 $1/f$ 꼬리로 눈에 보이게 나타나므로, 등분배처럼 조용히 답을 오염시키지 않는다. 저주파 평탄부의 높이를 이론값과 맞추면 검출기 교정 계수까지 덤으로 나온다.

대신 실무의 보정 항목이 길다.

별칭. 표본화 주파수의 절반 위 성분이 접혀 들어온다. 로렌치안을 그냥 맞추면 꼬리가 들려 보이므로 별칭된 로렌치안으로 맞춰야 한다.
검출기의 기생 필터링. 1064 nm에서 실리콘은 거의 투명해 광전자가 깊은 곳에서 느리게 확산해 나온다. 그 결과 수십 kHz 급의 의도치 않은 저역통과가 걸리고, 보정 없이 맞추면 $f_c$ 가 낮게 나온다.
파센 보정. 벽에서 거리 $h$ 안쪽이면 유체 저항 $\gamma$ 가 커진다. 표면에 붙은 세포를 다룰 때 $\gamma$ 를 벌크 값으로 쓰면 $k$ 가 그만큼 틀린다.
유체 기억. 엄밀히는 스펙트럼이 순수 로렌치안이 아니라 $\sqrt f$ 보정(유체 관성)이 붙는다. 정밀 측정에서는 이것까지 넣는다.
적합 방법. 주기도(periodogram)의 오차 분포는 가우시안이 아니라 지수분포다. 그래서 단순 최소제곱이 아니라 블록 평균 후 최대우도추정으로 맞추는 것이 정석이고, 이 처방을 정리한 베르그-쇠렌센·플리우비에르 논문이 사실상 이 분야의 표준 절차서다.

3. 드래그 힘 방법. 시료대를 알려진 속도 $v$ 로 흘리거나 정현파로 흔들어 $F=\gamma v$ 라는 알려진 힘을 주고 변위를 잰다. $k=\gamma v/\Delta x$ . 온도나 등분배 가정에 의존하지 않고, 정현파 구동이면 그 주파수에 뾰족한 봉우리가 서서 검출기 교정까지 동시에 된다. 다만 $\gamma$ 를 알아야 하는 것은 마찬가지다.

4. 볼츠만(퍼텐셜) 방법. 위치 히스토그램에서 $U(x)=-k_BT\ln P(x)$ 로 퍼텐셜 자체를 복원한다. 강성 하나가 아니라 퍼텐셜의 모양을 주므로 비조화성(큰 변위에서 트랩이 무르게 풀리는 구간)을 직접 볼 수 있다. 대신 꼬리까지 통계를 채워야 해서 측정 시간이 길고, 그동안 드리프트가 히스토그램을 뭉갠다.

네 방법이 서로 다른 가정을 쓰므로, 진지한 실험은 최소 둘을 돌려 교차 검증한다. 값이 20% 이상 어긋나면 잡음·드리프트·검출기 비선형 중 하나가 살아 있다는 신호다.

6. 수치적으로 굴려 볼 때[편집]

트랩 안의 구슬을 시뮬레이션하는 일은 확률미분방정식 적분의 교과서 예제다. 다만 이 문제에는 공짜 점심이 하나 있다.

정확한 갱신식이 존재한다. 오른슈타인-울렌벡 과정은 이산 전이 분포가 닫힌 형태로 알려져 있어서, 오일러-마루야마로 근사할 이유가 없다.

x_{n+1}=x_n e^{-\Delta t/\tau}+\sqrt{\frac{k_BT}{k}\left(1-e^{-2\Delta t/\tau}\right)}\;\xi_n, \qquad \tau=\frac{\gamma}{k}

$\xi_n\sim\mathcal N(0,1)$ 이고 이 식은 $\Delta t$ 가 아무리 커도 정확하다. 정상 분산과 자기상관이 구성상 정확히 맞으므로 교정 파이프라인 검증용 합성 데이터를 만들 때는 반드시 이걸 쓴다. 반대로 트랩이 비조화이거나 힘이 위치에 복잡하게 의존하면 정확해가 없어지고, 그때는 $\Delta t\ll\tau$ 를 지키는 일반 적분기로 돌아간다.

검증 지표. 합성 궤적의 PSD를 구해 코너주파수가 입력한 $k/2\pi\gamma$ 로 돌아오는지, $\langle x^2\rangle$ 이 $k_BT/k$ 인지 본다. 두 개가 동시에 맞으면 요동-소산 균형이 제대로 구현된 것이다. 하나만 맞으면 잡음 진폭이나 시간 척도 중 하나가 틀렸다.
자주 나오는 버그. 잡음 진폭에 $\sqrt{\Delta t}$ 를 빠뜨리는 것(그러면 $\Delta t$ 를 바꿀 때마다 온도가 바뀐다), 그리고 표본화율을 코너주파수의 10배 이하로 잡아 놓고 별칭 없이 로렌치안을 맞추려는 것.
힘 계산 쪽. 트랩 자체를 물리적으로 계산하려면 위 미 영역 이야기로 돌아가 $T$ -행렬이나 FDTD를 돌리고, 개방 영역이므로 완전정합층으로 경계를 닫는다. 실무에서는 힘을 매번 푸는 대신 한 번 계산해 3차원 힘 지도로 표를 떠 두고 브라운 동역학에서 보간해 쓰는 것이 보통이다.

7. 무엇을 하는 도구인가[편집]

단분자 생물물리. 이 분야를 통째로 만들어 낸 응용이다. 구슬을 손잡이 삼아 분자에 붙이고 힘을 걸면서 길이를 읽는 단분자 힘 분광이 표준 실험이 됐다.

키네신. 1993년 블록 그룹이 광집게 간섭 검출로 미세소관 위를 걷는 키네신 한 분자의 8 nm 계단을 직접 봤다. 정지력(stall force)은 5~7 pN 규모. 분자 모터가 “걷는다”는 것을 그림이 아니라 시계열로 확정한 실험이다.
DNA. 이중나선을 양끝에서 당기면 지속 길이 50 nm짜리 벌레형 사슬(WLC) 모형으로 잘 맞다가,³ 65 pN 부근에서 갑자기 1.7배로 늘어나는 과신장 전이가 나타난다. RNA 중합효소의 전사 진행, 리보솜의 번역, 단백질 접힘·펼침도 같은 장치로 잰다.
힘 클램프. 변위를 되먹임해 힘을 일정하게 유지하는 모드가 있는데, 되먹임 루프의 설계가 그대로 PID 제어 문제다. 루프 대역폭이 $f_c$ 근처면 트랩 자체의 동역학과 간섭해 발진하므로, 제어기 설계에서 코너주파수를 반드시 고려한다.

원자물리. 여기서는 같은 힘을 다른 이름으로 부른다. 공명에서 멀리 떨어뜨린 레이저의 기울기 힘이 광 쌍극자 트랩이고, 이탈량 $\Delta$ 에 대해 퍼텐셜 깊이는 $I/\Delta$ 로, 원치 않는 광자 산란률은 $I/\Delta^2$ 으로 스케일한다. 멀리 이탈시키고 세기를 올리는 것이 국룰인 이유가 이 비율이다. 마주 오는 빔의 정상파로 주기 퍼텐셜을 만들면 광격자가 되고, 초점을 여러 개 만들어 원자를 하나씩 배치하면 중성원자 양자 시뮬레이터의 트위저 배열이 된다. 보스-아인슈타인 응축의 마지막 증발 단계도 상당수가 자기 트랩이 아니라 광 쌍극자 트랩에서 진행된다.

언쇼와의 관계. 광집게가 성립하는 이유는 언쇼 정리의 전제를 비껴가기 때문이다. 정리가 금지하는 것은 고정된 극자의 정적 가둠인데, 광집게가 쓰는 것은 빛이 그 자리에서 유도한 쌍극자다. 유도 극자의 에너지는 $-\tfrac12\alpha\langle E^2\rangle$ 이라 세기의 극대점에서 안정할 수 있고, 라플라스 방정식이 금지하는 것은 퍼텐셜의 극값이지 $\langle E^2\rangle$ 의 극값이 아니다. 다만 방사압까지 넣으면 이야기가 한 번 더 꼬이는데(애슈킨-고든의 광학 언쇼 정리), 그 논증은 레이저 냉각 문서에 정리돼 있다.

8. 여담[편집]

애슈킨은 1980년대에 이미 살아 있는 박테리아와 바이러스를 가두는 데 성공했는데, 초기 실험에서 가시광 레이저를 쓰다가 세포를 익혀 죽이는 일이 반복됐다. 그가 “옵티큐션”(opticution, optical + electrocution)이라는 단어를 만들어 부른 그 사고다. 근적외선으로 파장을 옮긴 것이 해법이었고, 물의 흡수와 세포 광손상이 동시에 낮은 창이 대략 800~1100 nm 라는 사실이 지금까지 이 분야의 표준 파장을 정하고 있다.
노벨상까지 48년이 걸렸다. 애슈킨은 수상 당시 96세로 역대 최고령 수상자였고, 인터뷰 요청에 “지금 태양광 논문을 쓰는 중이라 바쁘다”고 답한 일화가 유명하다.
광집게의 힘 눈금이 pN 이라는 것은, 뒤집어 보면 모기 한 마리 무게의 100만분의 1 규모를 재고 있다는 뜻이다. 그런데 이 장치의 교정은 결국 “구슬이 얼마나 떨리는가”를 보는 것이므로, 페랭이 1908년에 현미경으로 세던 그 떨림과 원리적으로 같은 신호를 쓴다. 열잡음이 방해물이 아니라 눈금자로 쓰이는 드문 사례다.

9. 관련 문서[편집]

“빛이 물체를 민다”는 사실 자체는 케플러가 혜성 꼬리가 태양 반대편을 향하는 것을 보고 1619년에 추측했고, 맥스웰이 이론적으로 정리했으며, 레베데프가 1901년에 실측했다. 300년 가까이 걸린 셈인데, 이유는 힘이 너무 약해서 미세한 잔류 기체의 열적 흐름(라디오미터 효과)에 완전히 묻혀 버리기 때문이었다. 문방구에서 파는 라디오미터 바람개비가 도는 것은 광압이 아니다. ↩
트랩 강성이 레이저 출력에 비례한다는 것은 맞지만, 그 비례상수가 대물렌즈 수차·구슬 크기·구슬과 시료대 사이 굴절률 불일치에 심하게 흔들린다. 그래서 “같은 장비, 같은 출력, 다른 날”에 강성이 30% 달라지는 일이 예사이고, 실험마다 교정을 다시 하는 것이 이 바닥의 국룰이다. 논문에 강성 값만 적고 교정법을 안 적으면 리뷰어에게 반드시 지적당한다. ↩
지속 길이(persistence length)는 사슬이 방향을 기억하는 거리다. 이중나선 DNA는 50 nm(약 150 염기쌍) 규모라 그 이상 길이에서는 유연한 실처럼 굴고, 단일가닥은 1 nm 정도라 훨씬 흐물흐물하다. 이 차이 덕분에 힘-신장 곡선만 보고도 사슬이 이중가닥인지 단일가닥인지 구분할 수 있어서, 헬리케이스가 나선을 풀어내는 순간을 실시간으로 잡아낼 수 있다. ↩

광집게