SnSe2 결함 도입으로 인한 SnSe의 고온 열전성능 증대 메커니즘
The Mechanism Behind the High zT of SnSe2 Added SnSe at High Temperatures
Article information
Trans Abstract
SnSe is a promising thermoelectric material due to its low toxicity, low thermal conductivity, and multiple valence band structures, which are ideal for high electronic transport properties. The multiple valence band structure has attracted many attempts to engineer the carrier concentration of the SnSe via doping, to place its fermi level at a position where the maximum number of valence bands can participate in the electronic transport. Up until now, ~5 × 1019 cm-3 was the highest carrier concentration achieved in SnSe via doping. Recently, introducing SnSe2 into SnSe was found to effectively increase the carrier concentration as high as ~6.5 × 1019 cm-3 (at 300 K) due to the generated Sn vacancies. This high carrier concentration at 300 K, combined with the reduction in lattice thermal conductivity due to SnSe2 micro-domains formed within the SnSe lattice, improved the thermoelectric performance (zT) of SnSe – xSnSe2 as high as ~2.2 at 773 K. Here, we analyzed the changes in the electronic band parameters of SnSe as a function of temperature with varying SnSe2 content using the Single Parabolic Band (SPB) model. According to the SPB model, the calculated density-of-states effective mass and the fermi level are changed with temperature in such a way that the Hall carrier concentration (nH) of the SnSe – xSnSe2 samples at 773 K coincides with the optimum nH where the theoretically maximum zT is predicted. To optimize the nH at high temperatures for the highest zT, it is essential to tune the 300 K nH and the rate of nH change with increasing temperature via doping.
1. 서 론
에너지는 우리의 삶을 편리하고 풍요롭게 만드는 데 필수적인 요소이다[1]. 지구 온난화와 같은 환경 문제로 인해 기존의 화석연료를 사용하는 발전 방법 대신 친환경적인 방법으로 에너지를 생산할 수 있는 기술에 많은 관심이 쏠리고 있다. 현재 가장 대중적으로 활용되고 있는 전기 발전 방식인 화석연료를 연소하여 발전하는 화력발전은 생산되는 에너지의 약 60 % 이상을 폐열로 낭비하고 있다[2]. 때문에 폐열로부터 유용한 전기 에너지를 발전할 수 있는 유일한 친환경 발전 방법으로 열전 발전 기술이 더욱 주목을 받고 있다. 열전 발전에 사용되는 열전 소자의 효율은 소자 내 사용된 소재의 열전성능과 밀접하게 연관되어 있다. 따라서, 열전 소자의 발전효율 증대를 위해 다양한 열전소재의 열전성능을 높이는 것이 필수적이며 소재의 열전성능을 높이기 위한 많은 연구가 진행되고 있다. 소재의 열전 성능은 무차원 열전성능지수인 zT로 표현할 수 있다. zT (=S2σT / (κe + κl)) 수식 내 S는 제백계수, σ는 전기전도도, T는 절대 온도, κe는 전체 열전도도 중 전자/홀 등 전하의 기여분, 그리고 κl는 격자열전도도를 나타낸다[3]. 소재의 zT가 높을수록 동일한 열 에너지(온도 차)를 더 큰 전기 에너지로 변환할 수 있다. zT와 밀접한 관련이 있는 S2σ(파워팩터) 내 S와 σ가 상충 관계(trade-off)에 있어 파워팩터의 증가에 의한 zT 증가는 간단하지 않다[4]. 따라서, 많은 연구자들은 resonant state를 형성하거나, 에너지 차이가 작은 두 개의 밴드들을 수렴함으로써(band convergence) S와 σ 간의 상충 관계를 극복하기 위한 다양한 시도들을 해왔다[5-8].
파워팩터 증대 대비 상대적으로 직관적인 zT 증대 방법으로 κl를 감소시키는 전략이 있다. κl는 격자 내 포논의 전달과 직접적으로 연결되어 있는데, 다양한 방법(defect engineering, nanostructuring)으로 포논을 산란하여 열전소재의 κl을 감소시키면 zT를 확실히 개선할 수 있다[9-11]. κl는 zT 수식 내 다른 변수들과 어떤 연관도 없기 때문에 κl의 감소는 zT 증대를 보장한다.
여러 열전소재 중, SnSe는 높은 열전 성능을 가질 수 있는 소재임이 다양한 연구를 통해 이미 보고되었다[12,13]. 또한, SnSe는 열전 성능이 높을 뿐만 아니라 독성이 적어서 실제 열전소자에 적용하기도 적합한 물질이다[14]. SnSe의 열전 성능을 높이기 위해서 다양한 연구가 계속해서 진행되고 있다. SnSe의 나노 구조화를 통해 κl을 크게 감소시킴으로써 우수한 열전 성능을 얻은 연구가 보고된 바 있다[15,16]. SPS(spark plasma sintering)와 같은 가압 소결 방법을 통해 소재 내 이방성을 강화하여 열전 성능을 향상시킨 사례도 보고되었다[17,18]. 또한, band convergence 전략을 이용하여 파워팩터를 향상시킨 결과도 문헌에 보고되었다[19]. SnSe에 Bi, Na, 혹은 Ag를 도핑함으로써 높은 zT를 얻은 연구도 보고되었다[20-22]. 그 외에도, SnSe 2D 나노시트를 사용한 경우, 벌크 SnSe 분말을 사용한 경우 대비 열전 성능이 개선될 수 있다[23]. 최근 Qin et al.은 SnSe에 SnSe2를 주입하여 773 K에서 zT ~ 2.2를 보고하였다[24]. SnSe는 제일원리법을 통해 본질적인 결함이 존재함이 알려져 있다[25]. 이러한 결함은 전하 농도 조절을 어렵게 하기 때문에 기존 도핑 혹은 결정 구조 제어를 통해 도달할 수 있었던 최대 전하농도는 5 × 1019 cm-3 이었다[26]. 하지만 Qin et al. 은 SnSe2를 SnSe에 주입함으로써 6.5 × 1019 cm-3 이상의 전하 농도를 달성하였다. SnSe2를 주입하면서 생기는 추가적인 Sn 공공이 이러한 문제를 해결할 수 있도록 했기 때문이다. 이를 통해 최대 파워팩터가 기대되는 전하농도까지 전하농도를 증가시킬 수 있었다. 그 결과 파워팩터가 순수 SnSe 대비 약 35 % 증가하는 것을 보고하였다 (40 → 54 μW cm-1 K-2). 그뿐만 아니라, SnSe 격자 내 생성된 SnSe2 마이크로 도메인은 κl을 감소시키는데 중요한 역할을 하였다[24,27]. 하지만, SnSe2 주입에 따른 SnSe의 밴드 구조 변화에 대해서는 고찰된 바가 없다. 여러 밴드 인자 중, density-of-states effective mass (md*)에 대한 분석 결과를 보고하였지만, Single Parabolic Band (SPB) 모델을 이용하지 않고 축퇴 영역(degenerate regime)에서만 정확한 것으로 알려져 있는 아래의 근사식을 사용하였기 때문에 계산된 md* 내 오류가 있을 것으로 예상된다.
따라서 본 연구에서는 SPB 모델을 이용하여 SnSe 내 SnSe2의 주입량 변화에 따른 주요 밴드 인자의 변화를 탐구하였다. SPB 모델을 사용하여 md* [28], non-degenerate mobility (μ0), weighted mobility (μw) [29], B-factor와 같은 밴드 인자를 SnSe2 주입량(x, SnSe + xSnSe2 (x = 1,2,3,5))과 온도에 따라 계산하였다[30-32]. SPB 모델 계산 결과, md*, μ0, μw, 그리고 B-factor에서 모두 x의 증가에 따라 증가하다가 다시 감소하는 경향을 나타내었다. 온도가 증가함에 따라 파워 팩터와 비례하는 μw는 감소하지만 zT와 비례하는 B-factor는 증가하는 경향을 나타내었다. 이는 전기적 특성이 감소하였지만, 열적 특성인 κl의 급격한 감소로 인한 보상이 원인이다. B-factor는 zT와 비례하기 때문에 전하농도를 조절할 경우 높은 zT를 얻을 수 있는 지점을 추론할 수 있다. SPB 모델 계산 결과에 따르면 300 K의 경우, x = 2 일 때, 약 14%의 zT 증대가 예상된다 (0.93→ 1.08). 하지만, 고온으로 갈수록 실험 zT와 SPB 모델 예측 최대 zT가 일치하는 것을 볼 수 있다. 이는 고온으로 갈수록 전하농도가 SPB 모델 예측 최적의 전하농도로 수렴하기 때문이다.
2. 실험 방법
md* 계산을 위해 SPB 모델 대신 SPB 모델과 최대 10 % 미만의 오차 내 동일한 결과를 계산하는 아래의 간단한 수식을 사용하였다 (식 (1)) [28].
위 식 (1)에서 n은 chemical 전하 농도이다. 하지만, n을 직접적으로 측정하기 어렵다. 대신 홀 전하농도 (Hall carrier concentration, nH)를 n으로 변환하여 식 (1)에 대입해야 한다.
n과 nH 간의 비율은 rH (Hall factor)로 표현될 수 있는데 이를 S 측정값을 이용하여 간단하게 계산할 수 있다 (식 (2)) [33].
rH는 reduced fermi level(η)에 대한 함수이다. SPB 모델에서 S 또한 수학적으로 η 에만 의존하여 결정 된다. 그러므로 rH는 η 대신에 S에 대한 함수로 표현할 수 있다.
논문에서 추출한 실험값인 nH와 S, 온도를 이용하면 md* 를 계산할 수 있다. md*는 전자의 정지 질량인 m0을 단위 로 표현한다. nH를 식(1)에 넣어서 계산하면 부정확하기에 nH를 n으로 바꾸는 것은 필수적이다[34].
온도에 따른 μ0는 각각의 온도에서 SPB 모델로 얻을 수 있다. nH에 의존하는 μH도 SPB모델을 이용해 계산할 수 있다.
위 식에서 h와 Fn은 플랑크 상수와 n차수의 Fermi integral을 나타낸다 (아래 식 (5)).
nH는 식 (4)로 표현될 수 있고, 이는 md*와 η에 의존한다. μ0를 계산하기 위해서는 식 (3-5)를 사용해야 하고 nH 에 의존하는 μH가 필요하다. SPB 모델에 따르면 식 (1)로 얻은 md*를 이용할 수 있다. nH에 의존하는 μH를 계산하기 위해선, μ0와 md*가 필요하다.
S, σ, 및 T에 의존하는 μw는 식 (6)으로 표현할 수 있다[35]. 식 (6)에서 e, kB는 각각 전자의 전하량, 그리고 볼츠만 상수를 의미한다.
물리학적으로 μw는 μ0와 (md*/me)3/2의 곱으로 나타낼 수 있다.
B-factor는 식 (8)로 나타낼 수 있다. 식 (6)에서 얻은 μw와 실험에서 측정된 열전도도로부터 얻은 κl을 이용해 계산할 수 있다[36,37].
3. 결과 및 고찰
3.1 상태밀도 유효질량(md*) 계산
그림 1(a)는 SnSe + xSnSe2 (x = 1, 2, 3, 5)에서 x 변화에 따른 S 및 log10(n)의 변화를 보여주고 있다. 300 K 에서의 S 및 nH 측정값을 식 (1)을 이용하여 log10(n) 대 S의 그래프에 filled 심볼로 표시하였다. 식 (1) 내 md* 피팅을 통해 filled 심볼로 나타낸 측정값을 관통하는 계산값 또한 실선으로 표시하였다. 기존의 Seebeck pisarenko 그래프에서 흔히 사용하는 S 대 nH 그래프 대신 log10(n) 대 S 그래프를 사용한 이유는 md* 계산에 기존에 사용하던 SPB 모델이 아닌 식 (1)를 사용하였기 때문이다. 모든 degeneracy 영역에서 SPB 모델 계산결과와 동일한 결과를 낼 수 있는 식 (1)는 SPB 모델 내 복잡한 Fermi integral 계산을 거치지 않아도 된다는 장점이 있다. 또한 x-축에 S를 표현하는 경우, 샘플의 fermi level (η)을 더욱 더 직관적으로 파악할 수 있다는 장점이 있다. 300 K일 때, S가 가장 높은 경우는 SnSe+SnSe2 (x = 1)일 때이다 (192 μV K-1). x가 증가함에 따라, 전반적으로 S가 감소하는데 (x = 3 일 때 제외), 이는 x의 증가가 η의 증가를 동반하기 때문이다. 그림 1(a)에 따르면 x의 증가는 n의 증가 또한 야기하는데 n의 증가는 md* 및 η 변화를 동시 고려해야 한다.
그림 1(b)는 그림 1(a)에서 도출한 x 변화에 따른 md* 변화를 보여준다. md*는 x = 3일 때 최댓값(2.02 m0)을 갖고 x > 3 영역에서 다시 감소하는 것으로 계산되었다. x = 3 샘플의 md*(2.02 m0)가 x = 5 샘플의 md*(1.98 m0) 대비 높음에도 불구하고 x = 3 샘플의 n(7.6 × 1019 cm-3)이 x = 5 샘플의 n(7.7 × 1019 cm-3) 보다 낮은 이유는 x = 3 샘플의 η가 x = 5 샘플의 η 보다 작기 때문이다.
그림 1(c)는 온도에 따른 md*의 변화를 보여주고 있다. 모든 x 샘플들에 대해서 고온 영역(T > 650 K)에서 md* 가 크게 감소하는 것을 관찰할 수 있다. Qin et al.은 상온에서부터 700 K 부근까지 증가하던 S는 700 K 이상의 온도에서 감소하고, 700 K까지 감소하던 nH는 700 K 이상의 온도에서 다시 증가하는 결과를 모든 x 샘플에서 관찰할 수 있다고 보고하였다[24]. 그러면서 그 이유를 고온에서 더 많은 밴드(valence band)가 전기적 수송특성에 기여하기 때문이라고 하였다. 동일하게 그림 1(c)에서 고온에서의 md* 감소 또한 전기적 수송특성에 기여하는 밴드 수의 변화로 이해될 수 있다.
Qin et al.에 따르면 x = 2 샘플의 상온 md*는 약 1.9 m0로 계산되었다. 이 때 사용된 근사식은 아래와 같다.
위 수식은 축퇴 영역에서만 정확한 md*를 계산할 수 있는 수식으로 알려져 있다. S 측정값으로 얘기한다면 S가 150 μV K-1 미만일 때 SPB 모델 계산 결과 대비 오류가 10 % 미만인 md*를 계산할 수 있다. 하지만, x와 상관없이 상온 S가 최소 180 μV K-1 이상인 경우 (그림 1(a)), 식 (9)를 이용하여 계산한 md*와 SPB 모델을 이용하여 계산한 md* 간 차이가 클 수 있다. 그림 1(b)에서와 같이 SPB 모델 결과와 모든 축퇴 영역에서 유사한 md*(오류 10 % 미만)를 계산할 수 있는 식 (1)를 사용하여 계산한 결과, x = 2 샘플의 상온 md*는 약 2.02 m0로 식 (9)로 계산한 1.9 m0 대비 약 6 % 더 크게 계산되었다. 고온으로 갈수록 S 측정값이 커지면서 식 (9)를 사용하여 계산한 md*의 오류는 더욱 커지게 된다. 따라서 식 (9)를 이용하여 계산한 고온 x = 2 샘플의 md*는 약 700 K까지 md* ~ 1.9 m0를 유지하는 반면, 실제 식 (1)를 이용하여 정확하게 계산한 x = 2 샘플의 md*는 525 K에서 최대 md*인 ~2.0 m0에 도달한 후 773 K에서 1.6 m0까지 감소하는 것으로 계산되었다. 따라서, 보다 정확한 식 (1)로 계산하였을 때, 고온에서의 md*가 Qin et al.이 보고한 md* 대비 가벼운 것을 알 수 있다.
3.2 비축퇴 이동도(μ0) 계산
그림 2(a)는 SnSe + xSnSe2 (x = 1, 2, 3, 5)에서 x 변화에 따른 μH 및 nH의 변화를 나타낸다. Filled 심볼의 경우, Qin et al.이 문헌에서 보고한 측정값을 나타낸다[24]. 실선의 경우, 식 (3,4) 내 md*에 그림 2(b)의 결과를 대입하고, 나머지 μ0를 측정값(filled 심볼)에 피팅하여 얻었다. 각 x에 대응하는 실선들이 측정값(filled 심볼) 위에 정확하게 위치한다는 것은 피팅하여 얻은 μ0가 실제 소재의 밴드 구조를 잘 반영한다는 것을 의미한다. 우선, x의 증가에 따라 nH가 증가하는 것을 filled 심볼의 위치 변화를 통해 알 수 있다. μH의 경우, x = 2 샘플까지는 증가하지만 x > 2의 경우, SnSe2 도입량이 증가할수록 μH가 감소하는 것을 관찰할 수 있다. 예를 들어, x가 1에서 2로 증가하면서 μH는 4% 증가하지만 (156 → 162 cm2 V-1 s-1), x가 2에서 5로 증가할 때 μH는 22 % 감소하였다 (162 → 126 cm2 V-1 s-1). SnSe 내 도입된 SnSe2는 다양한 스케일(나노 스케일부터 마이크로 스케일)의 결함을 만들기 때문에 x가 증가할수록 μH는 감소하게 된다.
그림 2(b)는 그림 2(a)에서 실선과 filled 심볼 간의 피팅 시 사용된 μ0 계산값을 x에 대해 표시한 것이다. μ0는 비축퇴 상태에서 전하의 이동도를 의미한다. 즉, 전하 간 산란이 최소일 때 소재 내 결함과 전하간의 본질적인 상호 작용을 반영하는 밴드 인자이다. 300K에서 x = 2 샘플의 μ0가 다른 x 대비 가장 높은 것을 확인할 수 있다 (216.4 cm2 V-1 s-1). μH와 마찬가지로 x가 커질수록 μ0가 증가하다가 다시 감소하는 경향을 보인다. 최고점인 x = 2 샘플과 비교하면 최저점인 x = 5 샘플 일때 μ0가 약 25 % 감소하였다 (216.4 → 169.9 cm2 V-1 s-1). μ0는 전하와 포논 간의 상호작용을 나타내는 deformation 포텐셜(Ξ)에 반비례하는데 x ≥ 2 영역에서 μ0가 급격하게 감소하는 것은 Ξ의 급격한 증가와 연관되어 있다. Ξ가 클 경우, 전하와 포논간의 상호작용이 커서 전하의 이동을 포논이 방해하는 것이기 때문에 μ0는 감소하게 된다. 따라서, SnSe2의 주입량이 증가할수록 전하-포논 상호작용이 강해지는 것을 확인할 수 있다.
그림 2(c)는 온도 증가에 따른 μ0의 변화를 보여준다. 온도 증가에 따라 μ0가 감소하는 이유는 μ0가 T3/2에 반비례하기 때문이다. 773 K에서도 x = 2 샘플의 μ0가 가장 높고 (31.9cm2 V-1 s-1), x = 5 샘플의 μ0가 가장 낮은 것을 확인할 수 있다 (24.7cm2 V-1 s-1). 이를 통해 고온에서 Ξ 또한 상온에서의 Ξ와 유사하게 x = 5 샘플의 Ξ가 다른 샘플 대비 가장 낮음을 유추할 수 있다.
3.3 가중 이동도(μw) 계산
그림 3(a)는 SnSe + xSnSe2 (x = 1, 2, 3, 5)에서 x 변화에 따른 상온 μw의 변화를 나타낸다. μw는 두 가지 방법으로 계산할 수 있다. 첫번째 방법은 md*와 μ0를 각각 구한 후, 식 (7)을 이용하여 계산하는 방법이다. 이 경우, md*와 μ0를 각각 구하고 이를 다시 식 (7)에 대입해야 하는 번거로움이 있다. 두 번째 방법은 S 및 σ 측정값을 식 (6)에 대입하여 직접적으로 μw를 계산하는 방법이다. 측정값으로부터 μw를 바로 계산할 수 있기 때문에 첫번째 방법 대비 편리하다는 이점이 있다. 정확도 측면에서는 첫번 째 방법으로 계산한 μw가 더욱 정확하겠지만, 두 번째 방법으로 도출한 μw 또한 첫 번째 방법 대비 오류가 크지 않은 것으로 알려져 있다[35].
그림 3(a)에 나타낸 μw는 위에서 언급한 두번째 방법을 이용하여 계산하였다. Qin et al.이 보고한 온도에 따른 S 및 σ 측정값을 식 (6)에 대입하여 μw를 계산하였다. 그 결과, x = 2 일 때 다른 x 대비 가장 높은 μw를 보였다(603 cm2 V-1 s-1). x = 2일 때 μw가 최고점에 도달한 후 x가 증가함에 따라 μw가 감소하는 경향은 그림 2(b)에서 관찰한 x증가에 따른 μ0 변화 개형과 유사하다. 식 (7)에서 와 같이 μw는 μ0와 md*3/2의 곱으로 정의될 수 있는데, x ≥ 2 구간에서 x 증가에 따른 md*의 변화가 μ0 변화 대비 크지 않기 때문에 μw는 μ0 경향에 더욱 밀접하게 영향을 받는 것으로 분석된다. μw는 이론 파워팩터 최댓값에 직접적으로 비례하기 때문에 nH가 최적으로 맞춰진 상황에서 x = 2의 이론 최대 파워팩터가 다른 샘플 대비 가장 높을 것으로 예측된다.
그림 3(b)는 그림 3(a)에서 계산한 μw를 이용하여 nH에 대해 계산한 이론 파워팩터(실선) 및 Qin et al.이 보고한 파워팩터 측정값(filled 심볼)을 나타낸다[24]. 그림 3(a)에서 보인 x에 따른 μw 변화 개형과 동일하게 그림 3(b) 내 파워팩터 측정값 역시 x = 2 샘플에서 가장 높은 파워팩터 측정값(53.7 μW cm-1K-2)을 보였다. 그 후 x > 2 영역에서 x가 증가할수록 파워팩터 측정값이 감소하는 것을 관찰할 수 있다. 앞서 언급한 것처럼, μw는 이론 최대 파워팩터와 비례한다. 따라서 가장 높은 μw를 갖는 x = 2 샘플( 그림 3(a) 참조)의 이론 최대 파워팩터(56.5 μW cm-1K-2) 또한 다른 샘플들의 이론 최대 파워팩터 대비 높은 것을 확인할 수 있다. 그림 3(b)에 따르면 현재 x = 2 샘플의 nH를 7.2×1019 cm-3로 제어할 경우, 약 5 %의 파워팩터 개선효과를 기대할 수 있다 (53.7 → 56.5 μW cm-1K-2). 파워팩터는 열전소재의 전기적 특성을 대변하는데 SnSe 소재의 경우, x = 2 만큼의 SnSe2를 도입할 경우, md*와 μ0 의 극대화로 인한 파워팩터의 증가가 관찰되었다. SnSe2가 x = 2보다 많이 도입된 경우, 여러 스케일의 결함으로 인해 μ0가 감소했고 이는 다시 파워팩터의 감소로 연결되었다.
3.4 B-factor 계산
그림 4(a)는 그림 3(a)에서 계산한 SnSe + xSnSe2 (x = 1, 2, 3, 5) 샘플들의 μw 계산값 및 300 K κl 측정값(그림 4(b))을 식 (8)에 대입하여 도출한 B-factor를 나타낸다. 식 (8)에 의하면, B-factor는 μw/κl 비율에 비례하는 것을 알 수 있다. 그림 4(b)는 Qin et al.이 보고한 온도에 따른 κl 측정값을 나타낸다[24]. 300 K의 경우, x가 증가할수록 SnSe2 도입에 의한 다양한 스케일의 결함 증가로 κl가 감소하는 것을 관찰할 수 있다. 300 K에서는 x = 1(1.02 W m-1 K-1) x = 2 샘플(0.97W m-1 K-1) 간의 κl 차이가 크지만 300 K 보다 높은 온도에서 x = 1과 x = 2 샘플들 간의 κl 차이는 크지 않는 것을 볼 수 있다. 유사하게 x = 3과 x = 5 샘플들 간의 κl 차이 또한 전 측정 온도 영역에서 크지 않았다. 하지만 x = 1, 2 샘플들과 x = 3, 5 샘플들 간의 κl 차이는 고온(T > 650 K)을 제외하고 모든 온도영역에서 유의미한 차이를 보였다. 이는 x > 2 샘플에서 x ≤ 2 샘플에서 관찰된 Sn 결함 및 SnSe2 마이크로 도메인 농도가 크게 증가하기 때문인 것으로 예상된[24]. 그림 4(b)에서와 같이 300 K에서 x 증가에 따라 κl가 계속적으로 감소하지만, x 증가에 따른 κl 감소 분 대비 x 증가에 따른 μw의 감소분이 상대적으로 더 크기 때문에 B-factor의 x 변화에 따른 개형(그림 4(a)) 또한 μw의 x 변화에 따른 개형(그림 3(a))과 유사한 것을 확인할 수 있다. 예를 들어, 가장 높은 B-factor 값을 갖는 x = 2 샘플의 경우, B-factor는 1.25로 다른 샘플들 대비 가장 높았다.
그림 4(c)는 그림 4(a)에서 계산한 B-factor를 기반으로 예측한 300 K에서의 nH에 따른 이론 zT 변화를 보여준다 (실선). 더불어 Qin et al.이 보고한 zT 측정값을 filled 심볼로 표시하였다. 그림 4(c)에 다르면, x = 2 샘플의 zT 측정값(~0.93) 과 x = 3 샘플의 zT 측정값(~0.92)이 유사한 것을 볼 수 있다. 하지만, SPB 모델을 이용하여 계산한 이론 최대 zT는 x = 2 샘플의 이론 최대 zT(~1.08) 가 x = 3 샘플의 이론 최대 zT(~1.02) 대비 약 5.9 % 높은 것으로 계산되었다. 따라서, x = 3 샘플의 nH를 3.2 × 1019 cm-3로 제어한다면 약 5.9 %의 zT 증대가 가능할 것으로 예상된다.
3.5 고온 zT 예측
그림 5(a)는 온도에 따른 SnSe + xSnSe2 (x = 1, 2, 3, 5) 샘플들의 B-factor를 보여준다. B-factor는 온도가 증가함에 따라 증가하는 경향을 보이는데 이는 B-factor가 T5/2 에 비례하기 때문이다 (식 (8)). 773 K에서는 300 K에서와 달리 x = 3 샘플의 B-factor(4.08)가 x = 2 샘플(3.85) 대비 높은 것을 볼 수 있다. 온도 증가에 따른 B-factor 증가는 온도 증가에 따른 이론 최대 zT 증가로 이어진다. 그림 5(b-e)는 x = 1 - 5 샘플들의 T 및 nH에 따른 이론 zT(실선)를 나타낸다 (측정값은 심볼로 표시하였다). 그림 5(a)의 온도 증가에 따라 B-factor 증가와 마찬가지로 그림 5(b-e)의 이론 최대 zT가 x와 무관하게 T 증가에 따라 항상 증가하는 것을 볼 수 있다. 773 K에서 가장 높은 Bfactor를 갖는 x = 3 샘플(4.08)의 773 K에서의 이론 최대 zT는 ~2.2 다. 773 K에서의 실험 zT값 역시 ~2.2 인것을 볼 때, 추가 도핑 없이 온도 증가에 따른 md* 및 η 변화로 최대 zT를 구현하는 최적의 nH로 nH가 변화되었음을 확인할 수 있다. 예를 들어, x = 2 샘플의 온도에 따른 실험값의 nH 변화를 살펴보면 300 K에서 6.3 × 1019 cm-3이었던 nH가 773 K에서는 4.1 × 1019 cm-3로 약 35% 감소한다. 그런데 이 때, 이론 최대 zT를 얻을 수 있는 최적의 nH 또한 773 K에서 3.7 × 1019 cm-3로 예측되었기 때문에 최대 이론 zT와 유사한 실험 zT를 추가적인 도핑 없이 얻을 수 있었다.
Qin et al.은 SnSe2 도입을 통해 Sn 결함의 농도를 증가시킴으로써 SnSe의 300 K nH를 ~6.5 × 1019 cm-3 수준까지 증가시킨 것이 SnSe + xSnSe2 소재에서 zT ~ 2.2 (773 K)의 높은 열전성능을 얻는 중요한 이유 중 하나라고 주장하였다. 하지만, SPB 모델 계산 결과, 773 K 에서 높은 zT를 얻을 수 있었던 이유 중 하나는 773 K 에서 이론 최대 zT를 구현하는 최적 nH를 실험적으로 얻었기 때문이다. 300 K에서 773 K으로 온도가 증가하면서 변화하는 nH의 온도 증가에 따른 감소분은 도핑 원소 혹은 도입하는 결함의 종류에 따라 달라질 것으로 예상한다 . Qin et al.의 경우, SnSe2 결함 도입에 따른 300 K에서 의 nH 증가, 밴드 인자의 변화 및 SnSe2 결함이 도입된 상태에서 결정된 온도 증가에 따른 md* 및 η 변화 정도(= nH 변화)가 773 K에서 이론 최대 zT를 실험적으로 얻을 수 있었던 요인으로 분석된다.
5. 결 론
SnSe 내 SnSe2를 도핑할 경우, Sn 결함 및 SnSe2 마이크로 결함 생성으로 인해 전하농도가 높아지고, 격자 열 전도도(κl)가 감소됨으로써 773 K에서 최대 zT ~ 2.2를 얻었다는 문헌의 결과가 있었다. 본 연구에서는 SnSe2 도입에 따른 SnSe의 상태밀도 유효질량(md*), 비축퇴 이동도 (μ0), 가중 이동도(μw) 및 B-factor 변화를 SPB 모델과 SPB 모델에 기반한 간단한 수식을 이용하여 계산하였다. 열전소재의 전기적 특성을 대변하는 μw 및 이론 최대 zT 와 밀접한 연관이 있는 B-factor 모두 300 K에서 SnSe – 2SnSe2 (x = 2) 샘플이 가장 높았다. 하지만, 773 K에서는 고온 μw 및 κl 영향으로 SnSe – 3SnSe2 (x = 3) 샘플의 B-factor가 가장 높은 것으로 계산되었다 (4.08). x = 3 샘플의 B-factor로 계산된 이론 최대 zT와 실험 zT는 모두 ~ 2.2 이다. 300 K의 경우, x = 3 샘플의 이론 최대 zT와 실험 zT는 각각 ~1.02 및 ~0.92로 전하농도(nH)를 3.2 × 1019 cm-3로 감소시킬 경우, 약 11 %의 zT 증대가 가능할 것으로 계산되었다. 하지만 773 K의 경우, x = 3 샘플의 nH 측정값이 이미 이론 최대 zT가 구현되는 nH와 일치하기 때문에 최대 zT ~ 2.2를 실험적으로 얻을 수 있었다. 실험적으로 온도가 올라갈수록 nH가 감소하는 정도는 SnSe 소재 내 도핑 원소/소재의 종류 및 도핑량에 따라 달라질 것으로 예상된다. SnSe2 도입에 따른 300 K에서의 nH 증가 및 SnSe2 도입으로 인한 온도에 따른 nH 변화 정도 (온도에 따른 md* 및 페르미 준위 변화 정도)가 SnSe – xSnSe2 소재가 773 K에서 높은 zT를 가질 수 있었던 원인으로 분석된다.
Acknowledgements
This work was supported by the Basic Study and Interdisciplinary R&D Foundation Fund of the University of Seoul (2022) for Hyun-Sik Kim. This work was also financially supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2019R1A6A1A11055660, NRF-2015R1A6A1A03031833) for Won-Seon Seo and Heesun Yang. This work was also supported by the Technology Innovation Program RS-2022-00154720 funded by the Ministry of Trade, Industry & Energy (MOTIE, Korea) for Weon Ho Shin.