목차
안녕하세요. 인실리코 DX 사업본부입니다.
이전 글에서 브래그 법칙과 구조 인자를 통해 XRD(X-ray Diffraction) 분석의 이론적 원리를 살펴보았습니다. 시뮬레이션 환경에서 계산된 회절 패턴은 수학적으로 완벽합니다. 그러나 실험실에서 합성한 하얀 분말 샘플을 XRD 장비에 넣고 찍어보면, 노이즈·피크 브로드닝·불순물 혼입 등으로 인해 이론과 현실 사이의 간극을 마주하게 됩니다.
“내 실험 데이터는 왜 교과서(시뮬레이션)와 다를까”
“A성분과 B성분이 정확히 몇 퍼센트씩 섞인 걸까?”
이 난해한 질문에 답하기 위해, 오늘 우리는 Materials Studio의 ‘Reflex’ 모듈을 켭니다. 단일 물질의 회절 패턴 피팅부터 복잡한 유기 혼합물의 정량적 분석(QPA)까지, 원리와 실전을 오가며 튜토리얼을 마스터해 보겠습니다.
첫 단추: “어떤 빛으로 비추었는가?” (Radiation Source)
실험실에서 얻은 데이터와 시뮬레이션 데이터(.xsd)를 겹쳐보았을 때(Overlay), 피크의 위치가 완전히 어긋나 있다면 어떻게 해야할까요? 구조가 틀렸다고 실망하기 전에, 가장 먼저 ‘광원(Radiation)’을 확인해야 합니다.
동일한 결정 구조라도 사용한 광원(Cu-Kα, Mo, Neutron, Electron 등)에 따라 회절 각도(2θ)는 달라집니다. 예를 들어, 구리(Cu) 타겟을 사용한 X-ray의 파장은 약 1.540562 Å입니다. 이 환경에서 특정 피크가 10.7°에서 나타났다고 가정해 봅시다. 만약 파장이 1.0 Å인 중성자(Neutron)나 전자빔(Electron)으로 광원을 바꾸면, 같은 물질이라도 회절 피크는 6.9°로 완전히 이동해 버립니다.
실전 (How-to): 광원 맞추기
1. 예제로 제공된 청바지 염료 물질인 indigo_a.xsd 구조 파일을 엽니다.
2. Reflex Powder Diffraction 다이얼로그의 [Radiation] 탭으로 이동합니다.
3. Type을 X-ray (기본값)에서 Neutron이나 Electron으로 변경한 후 Calculate를 눌러 파장에 따라 패턴이 어떻게 급변하는지 직접 확인해 봅니다.
4. 비교를 위해 다시 X-ray (Copper)로 설정하여 실험 조건과 시뮬레이션 조건을 일치시킵니다.
동질이상(Polymorph) 판별과 피크 살찌우기 (Broadening)
물질의 화학식은 같아도 원자 배열이 다른 동질이상(Polymorph)이 존재할 수 있습니다. Materials Studio 라이브러리 내에 있는 실험 데이터 indigo_1.3cam을 불러와 indigo_b 구조의 시뮬레이션 패턴과 겹쳐보면 엉뚱한 위치에 피크가 뜨거나 빈 곳이 생깁니다. 하지만 indigo_a 시뮬레이션 결과를 겹쳐보면 피크의 위치가 마법처럼 일치합니다. “아, 내가 합성한 물질은 a형 결정이구나!”를 알게 되는 순간입니다.
하지만 위치만 맞았을 뿐, 모양은 여전히 다릅니다. 시뮬레이션 데이터는 바늘처럼 날카로운데, 실험 데이터는 넓게 퍼져(Broadening) 있습니다. 현실의 결정(Crystallite)은 무한하지 않고 유한한 크기를 가지며, 입자가 작아질수록 피크 폭(FWHM)은 뚱뚱해지기 때문입니다. (
)
게다가 장비 자체의 해상도 한계에 의한 오차도 발생합니다.
실전 (How-to): 현실의 크기와 장비 오차 반영하기
1. 입자 크기 반영: Reflex Powder Diffraction의 [Sample] 탭에서 Crystallite size 영역의 상자를 체크합니다. 
의 기본값인 
을 으로 줄여봅니다. 날카롭던 피크가 넓어지며 실험 데이터의 굴곡에 점차 맞아떨어지는 것을 볼 수 있습니다.
2. 장비 오차 보정 (Caglioti Equation): 이번엔 [Profiles] 탭으로 이동합니다. 장비에 의한 FWHM 확장을 보정하기 위해 
파라미터를 조절합니다. 예제에서는 해상도 파라미터인 값을 0.24로 늘렸을 때, 
이상의 고각 영역 피크들이 실험 데이터와 완벽하게 피팅(Fitting)되는 것을 확인할 수 있습니다.
혼합물의 레시피를 밝혀라: 리트벨트 정량 분석 (Rietveld QPA)
단일 물질의 피팅이 끝났다면, 이제 난이도를 대폭 높여 다상(Multi-phase) 혼합물 분석으로 넘어갑니다. 정체를 알 수 없는 하얀 가루에 2형 당뇨병 치료제 성분인 Acetohexamide Form A와 Form B, 그리고 기준 물질인 Corundum(산화알루미늄)이 섞여 있다고 가정합시다. “각 성분이 정확히 몇 퍼센트씩 들어있는가?”를 수치로 밝히는 것이 QPA(Quantitative Phase Analysis)의 목적입니다.
과거에는 특정 피크 하나의 높이만 비교하는 RIR(Reference Intensity Ratio) 방식을 썼지만, 피크가 겹치거나 배향성이 있으면 심각한 오류가 났습니다. Reflex 모듈은 개별 피크가 아닌 패턴 전체(Full Pattern)를 사용하여 이론 값과 실험 값의 오차(R_wp)를 최소화하는 가장 정밀한 Rietveld QPA 방식을 사용합니다.
실전 (How-to): Rietveld QPA 초기 세팅
1. 실험 스펙트럼 데이터인 Acetohexamide_Mixture.xcd를 엽니다.
2. Reflex Powder QPA 다이얼로그의 [Phases] 탭에서 Add new phases를 클릭하여 혼합물 후보군인 Form A, Form B, Al2O3(Corundum) 구조 파일(.xsd) 3개를 모두 로드합니다.
3. 구조적 자유도(DOF)를 주지 않고 초기 Refine(정밀화)을 실행합니다.
4. 분석 결과, 오차율을 나타내는 
값이 약 34.8%(구조 보정 적용에 따라 23% 내외)로 꽤 높게 나옵니다. 아직 실험값과 시뮬레이션 사이의 간극(Difference line)이 큽니다.
분자들의 ‘줄 서기’ 보정 (Preferred Orientation)
오차가 큰 이유는 현실의 분자들이 우리 예상대로 얌전하게 무작위로 쌓이지 않기 때문입니다. 특히 Corundum처럼 특정 축(축)이 다른 축에 비해 3배 이상 긴 형태의 결정들은 바닥에 떨어질 때 마치 동전처럼 특정 방향으로 눕거나 서려는 경향을 보입니다. 이를 우선 배향(Preferred Orientation)이라고 합니다. 배향성이 심하면 특정 각도의 피크만 비정상적으로 치솟게 됩니다.
Reflex에서는 이를 March-Dollase 함수를 통해 수학적으로 보정합니다. 배향 정도를 나타내는 파라미터 가 1이면 무작위, 1보다 작으면 축에 수직, 크면 수평으로 쏠렸다는 의미입니다.
실전 (How-to): 우선 배향(March-Dollase) 보정
1. Study table에서 각 상(Phase)을 선택하고 Edit parameters를 누릅니다.
2. [Sample] 탭에서 Preferred orientation을 찾아 보정 방식을 March-Dollase로 변경합니다.
3. Corundum의 경우 축 방향 계수(
)와 
체크박스를 활성화합니다 (예: 
, 
활성화).
4. 분자 내부의 Torsion(비틀림) 등 구조적 자유도(DOF)와 배향성 파라미터를 활성화하고 다시 Refine을 돌립니다.
5. 결과: 오차율()이 11.08%까지 극적으로 떨어집니다! 도출된 Form A, Form B, Corundum의 질량비도 각각 약 40.68%, 41.86%, 17.46%로 실제 조성(40:40:20)과 매우 유사하게 도출됩니다.
보이지 않는 진실, 절대 질량비 추적 (Internal Standard)
지금까지의 Rietveld 방식은 우리가 입력한 구조물(A, B, Corundum)의 총합을 무조건 100%로 가정합니다(정규화). 하지만 만약 샘플에 비정질(Amorphous) 성분이나 X-ray에 잡히지 않는 미지의 불순물이 섞여 있다면 어떻게 될까요? 실제 A성분은 30%밖에 없는데, 소프트웨어는 50%라고 착각할 수 있습니다.
이를 해결하기 위해 사용하는 끝판왕 기법이 바로 내부 표준(Internal Standard)법입니다. 질량을 정확히 아는 물질(Corundum)을 미리 샘플에 섞어서 찍고, 이 표준물질의 피크 강도를 절대 척도(Scale Factor)로 삼아 나머지 성분의 ‘절대 질량’을 역산하는 방식입니다.
실전 (How-to): Internal Standard 스케일링
1. 단일 Phase A 스펙트럼에 섞인 Corundum 피크의 통합 강도(Integrated intensity)를 계산합니다 (예: 8.80e+3).
2. 1000이라는 임의의 목표값을 강도로 나누어 스케일링 팩터(Scale factor)를 구합니다 (예: 1000 / 8.80e+3 = 0.1137).
3. [Pattern Processing] 탭에서 Background 피크를 직접 제거(Subtract)하고, 방금 구한 스케일링 팩터를 전체 데이터에 곱해 스케일을 보정합니다.
4. 스케일이 맞춰진 실험 패턴들을 혼합물 데이터와 함께 정밀화(Refine)하면, 결과값이 100%로 정규화되지 않는 절대 질량 비율(Absolute weight fraction)을 얻게 됩니다. 실제 예제에서 총 질량합은 92.14%로 측정되며, 이는 Micro-absorption 등의 이유로 놓친 질량이 7.86% 존재함을 의미합니다.
Materials Studio 도입이 필요한 이유
BIOVIA Materials Studio
소재 연구에서 XRD 분석은 단순한 패턴 비교를 넘어, 결정 구조 규명·다상 정량 분석·비정질 함량 추적까지 포함하는 고도화된 데이터 해석 과정입니다. Materials Studio는 이러한 전 과정을 하나의 통합 플랫폼에서 수행할 수 있도록 설계된 계산 기반 소재 분석 솔루션입니다. 특히 Reflex 모듈은 XRD 시뮬레이션, Rietveld 정밀화, 배향성 보정, 결정 크기 및 장비 확장 효과 반영, Internal Standard 기반 절대 정량 분석까지 지원합니다.
이는 실험 데이터와 이론 모델 사이의 오차를 체계적으로 줄이며, 분석 결과의 재현성과 신뢰도를 동시에 확보하는 데 핵심적인 역할을 합니다. 반복적이고 경험에 의존하던 수동 해석을 정량적·물리 기반 모델링으로 전환하는 것이 바로 Materials Studio 도입의 가장 큰 의미입니다.
마치며: 단순한 선(Line) 너머의 물리적 실체를 찾아서
여러분, XRD 데이터 분석은 단순히 컴퓨터 화면에 뜬 두 개의 그래프를 눈대중으로 겹쳐보는 단순 작업이 아닙니다. 광원을 맞추고(Radiation), 입자의 크기를 가늠하며(Broadening), 무작위가 아닌 분자들의 배향성(Preferred Orientation)을 보정하고, 비정질의 함량까지 파악(Internal Standard)해 나가는 이 모든 과정은 역동적인 추리 과정과 같습니다.
Materials Studio의 Reflex 모듈은 이론과 현실 사이의 오차를 논리적으로 줄여 나가며, 거친 1D 그래프를 ‘입자 크기 300 Å’, ‘A성분 40.68%’와 같은 살아있는 인사이트로 번역해 줍니다. 이 튜토리얼을 바탕으로 여러분의 서랍 속에 잠자고 있는 복잡한 실험 데이터를 Reflex에 올려보세요. 그동안 보이지 않던 혼합물 속의 진실이 숫자로 선명하게 나타날 것입니다.
자주 묻는 질문 (Q&A)
Q1. Rietveld 정밀화에서 R-factor는 어느 정도면 신뢰 가능한가요?
단순히 수치만 낮추는 것보다 구조 모델의 물리적 타당성이 중요합니다. 낮은 Rwp 결과가 반드시 정확한 상 비율을 보장하지 않으므로, Rwp 수치와 더불어 잔차(Residual) 그래프가 평탄한지, 그리고 도출된 구조가 물리적으로 타당한지 반드시 교차 검증해야 합니다.
Q2. 결정 크기와 strain 효과를 동시에 고려해야 하나요?
그렇습니다. 결정 크기는 피크의 폭(Broadening)에 영향을 주며, 균일 변형(Homogeneous strain)은 피크 위치 이동(Shift)을, 불균일 변형은 폭의 변화를 일으킵니다. Reflex에서는 Scherrer Equation과 Caglioti Equation을 통해 이를 단계적으로 최적화할 수 있습니다.
Q3. Internal Standard 없이도 정량 분석이 가능한가요?
가능하지만 주의가 필요합니다. Rietveld QPA는 내부 표준 시료 없이도 Scale Factor를 통해 결정질 상 사이의 상대적 Weight %를 산출할 수 있습니다. 하지만 비정질(Amorphous) 함량까지 포함한 절대 정량을 위해서는 Corundum 같은 표준 물질을 섞어 스케일을 조정해야 정확합니다.
Q4. 배향성이 강한 시료는 어떻게 준비하는 것이 좋습니까?
시료 회전, 미세 분쇄, 바인더 혼합 등을 통해 물리적 배향을 최소화할 수 있습니다. 또한 물리적인 처리(미세 분쇄 등) 외에도 S/W 내적으로 Reflex의 March-Dollase 보정 기능을 반드시 활용해야 합니다. 특정 축 방향으로 입자가 배열된 상태를 수치화하여(R0), 배향성 때문에 비정상적으로 높아진 피크 강도를 수학적으로 보정해 주어야 정확한 정량 분석이 가능합니다.
Q5. 시뮬레이션 패턴과 실험 데이터가 계속 불일치한다면?
가장 먼저 광원(Radiation) 설정이 실험 파장과 일치하는지 확인하십시오. 만약 특정 면에서 피크가 나타나지 않는다면 구조 인자(Structure Factor)에 따른 선택 규칙(Selection Rule) 혹은 소멸 법칙에 의해 피크가 사라진 것인지 검토해야 합니다. 또한, 격자 내 원자들의 점유율이나 위치 자유도(DOF)가 적절히 부여되었는지도 단계적으로 점검해야 합니다.
BIOVIA Materials Studio 프로그램 및 사용법이 궁금하시다면 언제든 인실리코로 연락 주세요!
아래의 블로그 구독과Youtube 구독 버튼을 눌러주시면 더 많은 영상을 확인하실 수 있습니다. 감사합니다 :)
