내가 쓴 가정이 다음 질문이 되었다 — 집속초음파와 미세기포 클러스터링

어떤 질문에서 시작했나

물리학 수업에서 운동량 보존 법칙을 배우던 중이었습니다. 두 물체가 충돌할 때 운동량은 보존되면서 에너지는 형태를 바꾸어 전달된다는 대목에서, 초음파가 생체 조직을 통과할 때에도 같은 방식으로 에너지가 전달되지 않을까 하는 의문이 생겼습니다.

이 의문이 향한 곳은 혈액뇌장벽이었습니다. 뇌 혈관에는 약물이 함부로 통과하지 못하게 막는 촘촘한 장벽이 있고, 2023년 승인된 알츠하이머 치료제 레카네맙은 이 장벽 때문에 임상 효과가 제한된다고 알려져 있었습니다. 저자는 레카네맙의 분자량이 약 148 kDa로 리핀스키 5법칙을 크게 벗어난다는 점을 확인했습니다. 즉 이 약은 확산으로 스며들 수 있는 크기가 아니고, 통과시키려면 장벽을 물리적으로 열어야 한다는 뜻입니다. 집속초음파가 미세기포를 진동·붕괴시켜 장벽을 여는 방식이 후보로 떠오른 이유가 여기에 있습니다.

1차 탐구 — 임계 음압을 직접 계산했다

첫 번째 탐구의 질문은 단순했습니다. 장벽을 열려면 음압이 얼마여야 하는가.

저자는 기포의 팽창과 수축을 기술하는 레일리-플레셋 방정식을 4차 룽게-쿠타 수치적분으로 구현한 웹 시뮬레이터를 직접 만들었습니다. 구동 주파수 1.0 MHz, 초기 반경 2.0 마이크로미터 조건에서 음압 세 값을 비교했습니다.

• 0.3 MPa: 최대 붕괴 압력 약 6 MPa. 개방 임계에 미달
• 0.5 MPa: 약 23 MPa. 안정형 캐비테이션 — 기포가 터지지 않고 규칙적으로 진동하는 상태 — 의 안전 개방 영역
• 0.8 MPa: 약 80 MPa. 관성형 캐비테이션으로 넘어가 조직 손상 위험

여기서 쓰인 기계지수는 음압을 주파수의 제곱근으로 나눈 값으로, 초음파가 얼마나 세게 기포를 흔드는지를 나타내는 지표입니다. 0.5 MPa는 기계지수 0.5에 해당했고, 이 값은 실제 전임상에서 보고된 0.5~0.6 MPa 범위와 맞아떨어졌습니다. 교과서 원리에서 출발한 계산이 실측치와 만난 것입니다. 논문의 결과를 자기 계산으로 다시 확인하는 이 방식은 논문 한 편에서 시작하기와 재현도 연구다에서 다루는 접근 그대로입니다.

전환점 — 자기가 쓴 한 문장

이 사례가 공유할 만한 지점은 계산 결과가 아니라 그다음입니다.

1차 탐구를 마무리하면서 저자는 한계를 적었습니다. "기포 하나가 무한한 유체 속에서 홀로 진동한다"는 이상화된 가정 위에 모든 계산이 서 있다는 것이었습니다. 그리고 실제 임상에서 쓰는 미세기포 조영제는 혈관 안에 수백에서 수천 개가 동시에 존재한다는 사실을 알게 됩니다.

보통 한계는 보고서 끝에 적고 잊는 문장입니다. 이 저자는 그 문장을 다음 탐구의 질문으로 바꿨습니다. 기포가 여럿이면 무슨 일이 벌어지는가. 한계를 쓰는 법이 말하는 한계 서술의 쓸모가 바로 이것입니다. 한계를 정확히 적어 두면 다음 사람이 이어받을 수 있고, 그 다음 사람이 자기 자신일 수도 있습니다. 같은 구조로 이어진 다른 사례는 TREM2 사슬에서도 볼 수 있습니다.

2차 탐구 — 기포들은 서로 뭉쳤다

진동하는 기포는 주변 유체에 압력장을 만들고, 그 안에 놓인 다른 기포는 힘을 받습니다. 이 힘을 이차 브제르크네스 힘이라고 부르며, 두 기포가 같은 위상으로 진동하면 끌어당기고 반대 위상이면 밀어내는 상호작용력입니다.

저자는 반경 1, 2, 3, 5, 8 마이크로미터의 다섯 종류 기포를 두고, 간격은 20 마이크로미터로 고정한 뒤 여섯 가지 기포 쌍의 상호작용을 계산했습니다. 구동 조건은 주파수 1.0 MHz, 음압 0.2 MPa로 잡았고 이는 기계지수 0.2에 해당합니다. 1차에서 도출한 0.5보다 낮게 둔 것은, 붕괴가 격렬하지 않은 조건에서 위상차 효과만 분리해 보기 위해서였습니다. 반경이 서로 다른 기포를 같은 규모에서 비교하려고 방정식을 무차원화했고, 붕괴 말기의 급격한 강성 때문에 수치가 발산하는 문제는 강성 방정식용 해법으로 처리했습니다.

결과는 여섯 조합 전부가 인력이었습니다. 부피 진동의 상관계수는 반경이 비슷할수록 1에 가까웠고(3-3 조합 1.000, 3-5 조합 0.963), 차이가 가장 큰 1-8 조합은 0.156으로 훨씬 약했지만 여전히 인력 쪽이었습니다. 뉴턴 제2법칙과 부가질량으로 추정한 접근 시간은 3-3이 약 1.08 마이크로초, 3-5가 약 1.01 마이크로초, 1-8은 약 11.80 마이크로초였습니다. 집속초음파 한 펄스가 수 마이크로초에서 수십 마이크로초라는 점을 생각하면, 크기가 비슷한 기포들은 펄스 하나 안에서도 뭉칠 수 있다는 뜻입니다.

여기서 결론이 1차로 되돌아옵니다. 뭉친 자리는 국소 기포 밀도가 높아지고, 밀도가 높아지면 유효 캐비테이션 강도도 커집니다. 개별 기포 하나를 기준으로 계산한 안전 마진이 실제로는 더 좁을 수 있다는 것입니다.

예상과 어긋난 결과를 그대로 실었다

선형 이론대로라면 공명주파수 — 기포가 가장 잘 흔들리는 고유 진동수 — 가 구동 주파수를 기준으로 서로 반대쪽에 있는 두 기포는 위상이 어긋나 척력을 받아야 합니다. 3-5 조합과 3-8 조합이 그 경우였습니다. 그런데 둘 다 인력이 나왔습니다.

저자는 이 불일치를 지우거나 조건을 바꿔 다시 돌리지 않았습니다. 기계지수 0.2 수준의 진동에서도 이미 비선형성이 커서 선형 공명 이론만으로는 위상 관계를 설명할 수 없다고 적고, 큰 진폭에서 부피 파형에 섞여 드는 배음 성분들의 상관관계가 기본 진동수의 위상차를 압도했을 가능성을 해석으로 제시한 뒤, 배음 성분별 분리 분석을 후속 과제로 남겼습니다. 예상과 다른 결과를 결과로 내놓는 태도에 대해서는 예상과 다른 결과에 따로 정리되어 있습니다.

스스로 밝힌 한계

이 탐구도 한계를 네 가지 적었습니다. 각 기포가 상대의 존재와 무관하게 진동한다고 두고 그 결과로 힘을 계산한 일방향 근사라는 점, 기포 간격을 20 마이크로미터 한 값으로 고정했다는 점, 혈관벽이라는 기하학적 구속을 무시하고 무한한 유체를 가정했다는 점, 그리고 모든 조합이 인력으로 나온 이유를 추가 계산 없이 정성적으로만 설명했다는 점입니다.

다음 사람에게

한계는 감점 요인이 아닙니다. 이 사례에서 한계는 두 번째 탐구를 통째로 만들어 낸 재료였습니다. 한 학생이 감당할 수 있는 크기로 질문을 자르고(적당한 크기의 주제), 자른 자리에서 무엇을 못 다뤘는지 정확히 적어 두면, 그 문장이 다음 탐구의 첫 문장이 됩니다. 보고서 마지막에 한계를 쓸 때 이 점을 기억하면 좋겠습니다.