포노닉 결정

고체역학 계산물리 물리 마지막 수정: 2026-08-17 04:34:07

1. 개요[편집]

포노닉 결정 Phononic Crystal
정의	탄성 물성이 주기적인 인공 구조 → 탄성파 밴드갭
예측	시갈라스·에코노무 (1992) · 쿠슈와하 등 (1993)
첫 실측	마드리드 셈페레 조각상 (마르티네스-살라 등, 1995)
기구 1	브래그 산란 — 격자상수 ≈ 반파장
기구 2	국소 공진 (류 등, 2000) — 격자상수 ≪ 파장
표준 계산	블로흐-플로케 경계조건 + FEM 고유값 해석

진동을 막는 고전적인 방법은 무겁게 만드는 것이다. 주기적으로 만드는 것도 방법이다.

포노닉 결정(phononic crystal)은 밀도와 탄성계수가 공간적으로 주기 배열된 인공 구조로, 특정 주파수 대역의 탄성파·음파가 어느 방향으로도 전파하지 못하게 만든 것이다. 그 금지 대역이 탄성파 밴드갭(elastic/phononic band gap)이다.

원리적 위치는 명확하다. 주기 퍼텐셜 속 전자가 블로흐 정리를 따라 밴드 구조를 갖고, 주기 유전율 속 광자가 같은 이유로 광결정이 되듯, 주기 탄성 매질 속 포논도 같은 논리로 밴드갭을 갖는다. 지배 방정식만 갈아 끼우면 된다.

\nabla\cdot\left(\mathbf{C}(\mathbf{r}):\nabla_{s}\mathbf{u}\right) + \rho(\mathbf{r})\,\omega^{2}\mathbf{u} = 0

$\mathbf{C}$ 는 탄성 텐서, $\nabla_s$ 는 대칭 기울기다. 다만 광자와 결정적으로 다른 점이 하나 있다 — 탄성파는 편광이 셋(종파 1 + 횡파 2)이고 계면에서 서로 변환된다. 이 결합이 완전 밴드갭 설계를 광결정보다 훨씬 까다롭게 만든다.

용어를 정리해 두면, 고체 속 탄성파를 다루면 포노닉 결정, 공기·물 속 음파를 다루면 소닉 결정(sonic crystal)이라 부르는 관행이 있다. 그리고 격자상수가 파장 수준이면 포노닉 결정, 파장보다 훨씬 작으면 음향 메타물질이라는 것이 두 분야를 가르는 기준이다. 아래에서 보듯 이 경계는 생각만큼 깔끔하지 않다.

2. 두 종류의 밴드갭[편집]

2.1. 브래그 산란형[편집]

정통 기구다. 주기 $a$ 인 격자에서 산란된 파들이 상쇄 간섭하려면 브릴루앙 영역 경계에서 조건이 맞아야 하고, 결과적으로 첫 갭의 중심 주파수는 대략 $f \approx c/(2a)$ , 즉 격자상수가 반파장쯤 될 때 열린다. 갭이 열리는 조건은 광결정과 판박이다.

임피던스 대비가 커야 한다. 밀도비와 음속비 둘 다. 강철 원기둥을 공기 중에 배열하면 대비가 $10^4$ 수준이라 갭이 쉽게 열리고, 에폭시 속 유리 구슬처럼 대비가 작으면 잘 안 열린다.
채움률에 최적점이 있다. 너무 성기거나 너무 빽빽하면 닫힌다.
격자 대칭성이 높을수록 유리하다. 브릴루앙 영역이 원에 가까운 삼각격자·벌집격자가 정사각격자보다 완전 갭을 잘 낸다.

역사적으로 시갈라스·에코노무(1992)와 쿠슈와하 등(1993)이 주기 탄성 복합체의 밴드갭을 이론적으로 예측했고, 실험적 확인은 뜻밖의 곳에서 나왔다. 1995년 마르티네스-살라 등이 마드리드에 서 있던 에우세비오 셈페레의 강철 봉 조각상 주위에서 음압을 재 보니 특정 주파수 대역이 죽어 있었다. 아무도 방음을 의도하지 않고 세운 미술품이 첫 소닉 결정 실험 시료가 된 셈이다.¹

브래그형의 치명적 한계는 산수에서 바로 나온다. 100 Hz 소음을 공기 중에서 막으려면 $\lambda\approx 3.4$ m 이라 격자상수가 1.7 m 여야 한다. 저주파일수록 구조가 비현실적으로 커진다.

2.2. 국소 공진형[편집]

2000년 류 등이 Science 에 낸 구조가 이 벽을 뚫었다. 납 구슬을 실리콘 고무로 감싸 에폭시에 격자로 박은 것인데, 단위 셀은 1.55 cm 인데도 400 Hz 근처에서 갭이 열렸다. 그 주파수의 파장은 격자상수의 두 자릿수 이상 크다 — 브래그 조건과 완전히 무관한 곳에서 갭이 생긴 것이다.

기구는 단순하다. 각 산란체가 질량-스프링 공명자이고, 그 고유진동수 바로 위에서 유효 질량밀도가 음수가 된다. 갭 위치를 정하는 것은 격자상수가 아니라 공명자 자체의 고유진동수이므로, 구조를 작게 유지한 채 갭을 저주파로 끌어내릴 수 있다. 이 논문이 사실상 음향 메타물질 분야의 출발점이고, 유효 파라미터 관점의 서술은 그 문서에 있으니 여기서는 반복하지 않는다.

대가는 대역폭이다. 공진 기반이라 갭이 좁고, 고무의 점탄성 손실이 크며, 인과율이 요구하는 두께 하한을 피해 갈 수 없다. “얇고 넓은 대역”은 물리가 금지한다.

항목	브래그형	국소 공진형
갭 위치 결정	격자상수 $a$	공명자 고유진동수
필요 크기	$a \sim \lambda/2$	$a \ll \lambda$
대역폭	상대적으로 넓음	좁음
손실 민감도	낮음	높음 (고무·점탄성)
주 무대	초음파·MHz·구조진동	저주파 차음·방진

실제 소자는 둘을 섞는다. 저주파는 공진으로, 고주파는 브래그로 막는 다중 갭 설계가 흔하다.

3. 밴드구조 계산 — 실무 절차[편집]

이 바닥의 표준 워크플로는 단위 셀 하나 + 블로흐-플로케 경계조건 + 고유값 해석이다. 유한요소법으로 하는 절차를 순서대로 적으면 이렇다.

1. 단위 셀을 모델링한다. 물성이 반복되는 최소 단위를 잡고 메시를 만든다. 여기서 첫 함정 — 마주 보는 경계면의 메시가 절점 단위로 일치해야 한다. 한쪽 면을 메시한 뒤 복사(sweep/copy)하는 것이 정석이고, 자동 메시 생성기에 맡기면 짝이 안 맞아 주기 조건이 걸리지 않는다. 메시 생성 단계에서 이미 승부가 난다.

2. 블로흐-플로케 조건을 건다. 격자 벡터 $\mathbf{a}_i$ 에 대해

\mathbf{u}(\mathbf{r}+\mathbf{a}_{i}) = \mathbf{u}(\mathbf{r})\,e^{\,i\mathbf{k}\cdot\mathbf{a}_{i}}

를 짝지어진 절점 자유도 사이의 선형 구속으로 부과한다. 단순 주기경계조건과의 차이는 위상 인자 $e^{i\mathbf{k}\cdot\mathbf{a}}$ 하나뿐이지만, 이것 때문에 계가 복소 에르미트가 된다. 상미분방정식 쪽의 플로케 이론이 공간 방향으로 옮겨온 것이 정확히 이 조건이다.

3. 파수 $\mathbf{k}$ 를 고정하고 고유값 문제를 푼다.

\left[\mathbf{K}(\mathbf{k}) - \omega^{2}\mathbf{M}\right]\mathbf{u} = \mathbf{0}

각 $\mathbf{k}$ 마다 일반화 고유값 문제이며, 최저 $10\sim20$ 개 밴드만 필요하므로 이동-역변환 란초스(ARPACK/SLEPc)가 정석이다.

4. 기약 브릴루앙 영역 경계를 따라 $\mathbf{k}$ 를 훑는다. 정사각격자면 $\Gamma\!-\!X\!-\!M\!-\!\Gamma$ , 삼각격자면 $\Gamma\!-\!M\!-\!K\!-\!\Gamma$ . 결과를 $\omega$ vs $k$ 로 그린 것이 밴드 다이어그램이고, 모든 밴드가 비어 있는 가로띠가 밴드갭이다.

5. 성능 지표를 뽑는다. 갭 폭을 중심 주파수로 나눈 갭-중심 비(gap-to-midgap ratio) $\Delta\omega/\omega_c$ 가 표준 지표다. 10%를 넘으면 쓸 만하다고 본다.²

6. 유한 구조로 검증한다. 밴드 다이어그램은 무한 주기 구조의 성질이다. 실제 소자는 셀이 5~10개뿐이므로, 유한 개 셀을 늘어놓고 투과 스펙트럼을 계산해 예측한 갭 위치에서 실제로 몇 dB가 떨어지는지 확인해야 한다. 이 단계를 건너뛴 논문은 대체로 현장에서 배신당한다.

4. 완전 밴드갭 vs 방향 갭[편집]

밴드 다이어그램을 보고 흔히 저지르는 오독이다.

방향 갭(directional gap, stop band): $\Gamma\!-\!X$ 처럼 특정 방향으로만 전파가 금지된 경우. 다른 방향으로는 멀쩡히 지나간다. 만들기 쉽고, 1차원 적층 구조(브래그 미러)만으로도 얻는다.
완전 밴드갭(complete band gap): 모든 방향, 모든 편광에 대해 동시에 닫힌 경우. 브릴루앙 영역 경계 전체를 훑었을 때 어떤 밴드도 그 주파수대에 없어야 한다.

고체에서 완전 갭이 어려운 이유가 편광이다. 곧은 원기둥이 박힌 2차원 구조에서는 면외 전단(SH) 모드와 면내(종파+SV 결합) 모드가 분리되므로, 두 계열 각각의 갭이 주파수대에서 겹쳐야 완전 갭이 된다. 한쪽만 보고 “갭이 열렸다”고 하는 것이 흔한 실수다. 반대로 유체를 모재로 쓴 소닉 결정은 유체가 전단을 지지하지 않아 편광 문제가 사라지고, 그래서 강철 봉을 공기 중에 세우기만 해도 갭이 나온다.

5. 계산의 함정[편집]

$\Gamma$ 점의 특이성. $\mathbf{k}=0$ 에서 강성행렬이 강체 운동 3개(2D면 2개)에 대해 특이하다. 음향 분기가 $\omega=0$ 에서 출발하기 때문인데, 이동을 0으로 두면 인수분해가 터진다. 작은 음의 이동을 주는 것이 표준 대응이다.
밴드 극값이 영역 경계에 없을 수 있다. 기약 브릴루앙 영역 경계만 훑는 관행은 편의일 뿐 정리가 아니다. 밴드의 최대·최소가 영역 내부에 있는 반례가 실제로 존재하고, 그러면 경계 스윕은 있지도 않은 갭을 보고한다. 갭이 아슬아슬할 때는 2차원 $\mathbf{k}$ 격자 전체를 훑어 확인해야 한다.³
$\omega(\mathbf{k})$ 와 $\mathbf{k}(\omega)$ 는 다른 문제다. 표준 절차는 실수 $\mathbf{k}$ 를 넣고 $\omega$ 를 얻지만, 실무가 알고 싶은 것은 갭 안에서 셀당 몇 dB로 죽느냐, 즉 복소 파수다. 이건 $\omega$ 를 고정하고 $\mathbf{k}$ 를 미지수로 두는 이차 고유값 문제가 되고, 감쇠상수 $\mathrm{Im}\,k$ 가 직접 나온다. 확장 평면파 전개(EPWE)나 상응하는 FEM 정식화가 여기 쓰인다.
손실을 빼면 국소 공진형 예측이 크게 빗나간다. 고무의 점탄성 손실계수는 무시할 수준이 아니다. 복소 탄성계수를 넣으면 고유값도 복소가 되고 밴드 개념 자체가 흐려지므로, 손실이 큰 구조는 애초에 유한 구조 투과 해석으로 판단하는 편이 낫다. 감쇠 모형 선택이 결과를 좌우한다.
메시 해상도. 최상단 밴드의 파장이 기준이다. 2차 요소 기준 파장당 6~10개는 되어야 하고, 이걸 안 지키면 고주파 밴드가 위로 밀려 갭이 실제보다 넓게 나온다. 메시 독립성 연구를 생략할 수 없다.
평면파 전개(PWE)는 물성 대비가 크면 느리다. $\rho$ 와 $\mathbf{C}$ 가 불연속이라 푸리에 급수가 깁스 현상을 겪고, 고체-유체 혼합계에서는 특히 나쁘다. 그래서 광결정 쪽과 달리 포노닉 결정은 일찍부터 FEM이 주력이 됐다. 1차원이면 전달행렬법, 원기둥·구 배열이면 다중산란 이론(MST)이 정확한 대안이다.

6. 응용[편집]

진동 차단과 방진. 정밀 장비 지지구조, 위성 페이로드 어댑터, 공작기계 베드처럼 특정 대역을 죽이면 되는 곳이 1순위 고객이다. 모드 해석으로 잡은 문제 주파수에 갭 중심을 맞추는 설계가 위상 최적화와 결합해 자동화되고 있다.

방진 트렌치와 지반 구조물. 지반에 주기적으로 파낸 시추공·트렌치·말뚝 배열로 표면파를 막으려는 시도가 있다. 2014년 브륄레 등이 그르노블 인근에서 5 m 깊이 시추공 배열을 실제로 파고 50 Hz 가진을 걸어 표면파가 되돌려지는 것을 측정한 실험이 유명하다. 다만 지진 공학의 관심 대역(0.5~5 Hz) 은 파장이 수십~수백 m 라 브래그형으로는 격자가 비현실적이고, 실용화 논의는 공진형이나 지반 자체를 공명자로 쓰는 쪽으로 간다. 이 분야는 아직 개념 실증 단계이며, 재래식 차수벽·방진 트렌치 대비 우위가 확립됐다고 말하기는 이르다.

MEMS·RF 공진기. 실질적으로 가장 상용에 가까운 분야다. 램파 공진기나 FBAR의 $Q$ 를 깎아먹는 주범이 지지대(anchor)를 통한 음향 에너지 누출인데, 지지대 경로에 포노닉 결정을 심어 동작 주파수를 밴드갭 안에 두면 누출이 급감한다. MEMS 공진기 설계에서 “음향 반사경”이라는 이름으로 이미 쓰이고 있고, 같은 원리로 SAW 소자의 반사기·지연선도 만든다.⁴

열포논 공학. 격자상수를 수십 nm 로 줄이면 열을 나르는 포논 대역에 손댈 수 있다는 발상이다. 실리콘 나노메시로 열전도도를 실제로 한 자릿수 낮춘 결과들이 있고, 열전 소자의 $ZT$ 를 올리는 경로로 주목받았다. 여기서는 과장을 특히 조심해야 한다.

상온 열포논은 파장이 1~10 nm 수준이고 스펙트럼이 매우 넓다. 이 전 대역에 밴드갭을 여는 것은 원리적으로 불가능하다. 노릴 수 있는 것은 일부 대역의 군속도 저하와 산란 증가다.
상온에서 관측된 열전도도 감소의 상당 부분은 밴드 접힘 같은 결맞음 효과가 아니라 표면·계면에서의 비결맞음 산란으로 설명된다는 것이 현재의 중론이다. 결맞음 효과가 분명하게 확인된 것은 주로 극저온 영역이다.
즉 “포노닉 결정이 열을 밴드갭으로 막았다”는 서술은 대개 부정확하다. 표면적을 늘려 산란시킨 것과 구분하려면 온도 의존성과 특성 길이 실험이 필요하다.

7. 관련 문서[편집]

조각가 본인은 강철 봉을 격자로 세운 이유가 순전히 조형적이었을 것이다. 물리학자들이 논문 도판에 미술품 사진을 넣고 Nature 에 실은 몇 안 되는 사례. ↩
갭-중심 비를 안 쓰고 “밴드갭 폭 3 kHz”라고만 적힌 결과는 비교가 불가능하다. 스케일 불변성 때문에 구조를 두 배 키우면 모든 주파수가 절반이 되므로, 절대 폭은 아무 정보가 없다. 이 실수는 놀라울 정도로 자주 반복된다. ↩
주기 연산자의 스펙트럼 가장자리가 브릴루앙 영역 내부에 놓일 수 있다는 것은 수학 쪽에서 반례와 함께 정리된 사실이다. 밴드 다이어그램의 관행적 “경계만 훑기”는 계산 시간을 아끼는 휴리스틱이지 보증이 아니다. 갭 폭이 1% 대인 결과를 들고 오는 사람에게 먼저 물어봐야 할 것. ↩
완전 밴드갭이 열려도 실제 부품은 어딘가에서 지지되어야 하고, 그 지지대가 대개 갭을 우회하는 경로가 된다. “구조는 완벽한데 나사가 다 망쳤다”가 이 분야 실험의 국룰 결말이다. ↩