본 논문에서는 비평형 분자동역학 시뮬레이션 기법을 사용하여 알루미늄 박막과 실리콘 웨이퍼 간 열경계저항을 예측하였다. 실리콘의 끝 단 고온부에 열을 공급하고, 같은 양의 열을 알루미늄 끝 단 저온부에서 제거하여 경계면을 통한 열전달이 일어나도록 하였으며, 실리콘 내부와 알루미늄 내부의 선형 온도 변화를 계산함으로써 경계면에서의 온도 차이에 따른 열저항 값을 구하였다. 300K 온도에서 5.13±0.17m2·K/GW의 결과를 얻었으며, 이는 열유속 조건의 변화와 무관함을 확인하였다. 아울러, 펨토초 레이저 기반의 시간영역 열반사율 기법을 사용하여 열경계저항 값을 실험적으로 구하였으며, 시뮬레이션 결과와 비교·검증하였다. 전자빔 증착기를 사용하여 90nm 두께의 알루미늄 박막을 실리콘(100) 웨이퍼 표면에 증착하였으며, 유한차분법을 이용한 수치해석을 통해 열전도 방정식의 해를 구해 실험결과와 곡선맞춤 함으로써 열경계저항을 정량적으로 평가하고 나노스케일에서의 열전달 현상에 관한 특징을 살펴보았다.

탄소나노튜브(CNT) 기반의 멤브레인은 높은 물 전달률과 직경에 따른 이온 배제율로 해수담수, 물질 정화 등을 위한 분리막으로써의 가능성을 보여 주었다. 이온 선택성은 CNT 기반 멤브레인의 응용 분야를 확대하기 위한 중요한 요소이며, 기능기를 이용하여 이온 선택성의 조절이 가능함이 보고되었다. 다양한 원자가/크기의 이온이 혼합될 경우, 이온-기능기간 작용력 뿐만 아니라 이온-이온간의 작용력, 이온의 크기에 의한 반발력 등이 복합적으로 작용한다. 이에 본 연구에서는 분자동역학 전산모사를 통하여, 상이한 원자가/크기를 가진 이온의 혼합이 기능화된 CNT의 이온 선택성에 미치는 영향을 연구 하였다. Potential of Mean Force 계산을 통하여 이온 투과에 대한 자유 에너지 장벽을 계산하였으며, CNT 크기 변화, 전하량 변화를 통하여 이온 선택성과 배제에 영향을 미치는 요소를 분석하였다. 본 연구는 CNT 멤브레인을 이용한 분리막 설계, 생체 이온 전달 채널 모사 등에 유용할 것으로 기대한다.

고분자는 많은 수의 원자로 구성 되어 있기 때문에, 분자동역학 연구 시 다양한 변수를 고려 해야 하는 어려운 연구 대상이다. 구성 원자 개수가 너무 많은 경우 계산 시간이 오래 걸리는 문제점이 있으며, 반대로 원자의 개수를 너무 줄일 경우에는 실제 모사 대상을 정확하게 반영하지 못할 수가 있다. 특히, 기체 분리와 같은 분야에서는 기체 분자가 고분자 내에서 투과하는 것을 반영한 계산이 필요하기 때문에 오랜 시간동안 계산을 하게 되며, 이로 인하여 모델의 크기를 정하는데 있어서 앞서 언급한 모든 요소를 적절히 고려하는 것이 필요하다. 본 연구에서는 분자동역학 전산모사에 있어서 고분자를 구성하는 원자의 개수와 그에 따른 영향을 연구하기 위하여, 고분자 주쇄 길이를 조절하면서 고분자 모델의 특성을 분석하였다.

고분자 분리막의 분자동역학 연구에서는 구성 원자 개수가 많고 투과 거동 계산시 긴 시간을 필요로 하기 때문에 적절한 고분자 주쇄의 길이를 선택하는 것이 매우 중요하다. 본 연구에서는 이러한 고분자 주쇄 길이와 투과 거동 간의 상관 관계가 실제 분자동역학에서 어떻게 나타나는지 조사를 하고자 하였다. 널리 알려진 상용 고분자 Kapton® 폴리이미드 구조 를 이용하여 분자동역학을 수행하였고 기체 투과 거동을 분석하였다. 고분자 주쇄의 움직임은 그 길이와 큰 연관성이 없었으 며, 일반적인 인식과 달리 짧은 주쇄 길이를 갖고 있다고 해서 더 활발하게 움직이는 것은 아니라는 것을 확인할 수 있었다. 또한, 고분자 주쇄의 말단기는 상대적으로 움직이기 쉬울 것이라는 예상과 달리, 말단에 위치하지 않은 경우라도 말단기에 위치한 원자보다 더 높은 움직임을 보이는 경우도 많았다. 최종적으로 기체 분자의 투과 성능에서도 고분자 주쇄의 길이 및 말 단기에 따른 영향은 관찰되지 않았다. 이는 기체 투과 전산모사에서 많이 언급되는 말단기 효과를 실제 전산모사에 적용할 경우, 각각의 모델 특성에 따라 제한적으로 적용을 하고 이에 대한 검증 과정을 반드시 수행하여야 한다는 것을 의미한다.

연료전지용 전해질막의 성능에 있어서 가장 중요한 요소는 수소이온이 전해질막 내부에 형성된 수화채널을 따라 서 얼마나 빨리 전달될 수 있느냐이다. 여기에는 수화채널의 모폴로지 및 수소이온의 확산도 등이 매우 중요한 요소가 되는 데, 이를 규명하기 위하여 다양한 분자동역학 전산모사 연구가 진행되고 있다. 분자동역학 계산에 있어서 각 원자의 움직임 및 상호작용을 미리 변수화 시켜 놓은 force-field는 필수 요소 중 하나로서, 본 연구에서는 이러한 force-field의 종류가 전해 질막 전산모사에 미치는 영향을 분석하기 위하여, 다양한 force-field를 이용하여 연료전지용 전해질막의 수소이온 확산도를 계산하였다. 이 과정에서 non-bonding interaction을 결정하는 전하 값이 수화채널 모폴로지 형성에 매우 중요한 역할을 한다 는 것이 밝혀졌으며, COMPASS force-field가 가장 정확한 수소이온 확산도 값을 얻음으로써 연료전지용 전해질막의 전산모 사에 있어서 가장 적절한 force-field일 것으로 판단된다. 이러한 force-field의 적절한 선정은 최종 분자 구조 뿐만 아니라 수 소이온 확산도에도 큰 영향을 주는 것을 알 수 있었으며, 연료전지용 전해질막 전산모사 수행 시에는 이러한 부분을 충분히 감안하여 force-field를 선택하여야 할 것이다.

기존의 전기 탈이온 공정(EDI)은 전기투석법과 이온교환수지법을 혼합한 공정으로 알려져있다. 여기서의 모듈은 전기투석을 위한 양이온교환막과 음이온교환 막이 있고 이 막 사이에 이온교환수지를 채워 넣는 형태로 제조된다. 이에 이용 되는 양이온교환막으로 술폰화시켜 이온교환능을 가지게한 sulfonated SEBS triblock copolymer를 이용했다. 게다가 술폰화시킨 SEBS, 물 분자, 히드로늄 이온을 모델링하고 술폰화시킨 SEBS, 히드로늄 이온과 물분자의 의 동역학을 분자동역학시뮬레이션을 이용해서 살펴보았다. 그리고 실제 실험을 통한 비교분석 을 진행했다.

수소이온 교환막(PEM; Proton Exchange Membrane)은 연료전지 막-전극 복합체(MEA; Membrane-electrode Assembly)를 구성하는 핵심 소재 중 하나로서, 촉매와 함께 연료전지 성능을 결정하는 중요한 역할을 한다. 이러한 수소이온 교환막의 성능은 내부에 존재하는 수소이온 전달 통로인 수화 채널의 구조에 큰 영향을 받는 것으로 알려져 있다. 분자 동역 학(MD; Molecular Dynamics) 전산모사 기술은 이러한 소재 내부의 분자 및 원자구조를 파악하기 위한 유용한 도구로서, 수 소이온 교환막의 구조 및 특성에 관한 많은 관련 연구가 진행되고 있다. 본 총설에서는 분자동역학 전산모사 관련 연구에 대 한 동향을 정리하고, 이를 통해 어떤 구조적 특징들을 분석할 수 있는지 제시하여, 수소이온 교환막 연구자들과 분리막 연구 자들에게 분자동역학 전산모사 기술의 유용성에 대하여 소개하고자 한다.

전기 탈이온 공정(EDI)은 전기투석법과 이온교환수지법을 혼합한 공정으로 알려져 있고, 모듈은 전기투석을 위한 양이온교환막과 음이온교환막 사이에 이온 교환수지를 채워 넣는 형태로 제조된다. 모듈 제조에 사용되는 양이온교환막을 고무상의 고분자 SEBS triblock copolymer를 이용하여 sulfonation을 진행하여 S-SEBS 양이온교환막을 제조하였다. 제조된 양이온교환막의 특성평가를 진행하였고, 술폰화시킨 SEBS, 물 분자, 히드로늄 이온을 모델링하여 술폰화시킨 SEBS, 히드로늄 이온과 물분자의 의 동역학을 분자동역학시뮬레이션을 이용해서 살펴 본 뒤 비교분석을 진행했다.

기존의 전기 탈이온 공정(EDI)은 전기투석법과 이온교환수지법을 혼합한 공정으로 알려져있다. 여기서의 모듈은 전기투석을 위한 양이온교환막과 음이온교환막이 있고 이 막 사이에 이온교환수지를 채워 넣는 형태로 제조된다. 이것을 기초로 요즘 전기 흡착 탈이온(EAD) 공정은 바이폴라막과 이온교환수지를 이용해서 모듈을 제조하고 있다. 그래서 모듈에 이용되는 양이온교환막으로 술폰화시킨 SEBS triblock copolymer를 이용했다. 게다가 술폰화시킨 SEBS, 물분자, 히드로늄 이온을 모델링하고 술폰화시킨 SEBS, 히드로늄 이온과 물분자의 의 동역학을 분자동역학시뮬레이션을 이용해서 살펴보았다. 그리고 실제 실험을 통한 비교분석을 진행했다.

기존의 전기 탈이온 공정(EDI)은 전기투석법과 이온교환수지법을 혼합한 공정으로 알려져있다. 여기서의 모듈은 전기투석을 위한 양이온교환막과 음이온교환막이 있고 이 막 사이에 이온교환수지를 채워 넣는 형태로 제조된다. 이것을 기초로 요즘 전기 흡착 탈이온(EAD) 공정은 바이폴라막과 이온교환수지를 이용해서 모듈을 제조하고 있다. 그래서 모듈에 이용되는 양이온교환막으로 술폰화시킨 SEBS triblock copolymer를 이용했다. 게다가 술폰화시킨 SEBS, 물 분자, 히드로늄 이온을 모델링하고 술폰화시킨 SEBS, 히드로늄 이온과 물분자의 의 동역학을 분자동역학시뮬레이션을 이용해서 살펴보았다. 그리고 실제 실험을 통한 비교분석을 진행했다.

본 논문에서는 브리징 스케일 분해를 기반으로 멀티스케일 문제에 대한 설계민감도 해석법을 개발하였다. 나노 기술의 급속한 발전으로 인해 나노 수준의 해석의 필요성이 지속적으로 증가하고 있다. 최근 분자동역학과 연속체역학의 연성문제에서 많은 해석 방법들이 개발되었다. 본 논문에서는 연성시스템 해석을 위해 브리징 스케일 기법을 사용한다. 전체 영역의 분자동역학 시스템의 해석은 많은 양의 계산 비용이 들기때문에 분자동역학과 연속체 시뮬레이션의 연성시스템을 선호한다. 분자동역학과 연속체 수준 사이의 정보 교환은 분자동역학과 연속체의 경계에서 일어난다. 브리징 스케일 법에서 일반화된 랑지벵 방정식은 축소된 영역의 분자동역학 시스템 해석을 위하여 요구되고, 시간이력 커널을 사용하여 구한 GLE 힘은 분자동역학 시스템에서 경계에 있는 원자들에 작용한다. 그러므로 분자동역학과 연속체 수준의 시뮬레이션을 분리해서 해석할 수 있으며 계산 과정을 가속시킬 수 있다. 연성문제의 시뮬레이션 이후에는 설계의 최적화를 위해 설계민감도 해석의 필요성이 자연스럽게 나타나며 전체 시스템의 성능은 나노 스케일의 효과를 고려해서 최적화된다. 설계구배 기반 최적화에서 설계민감도가 요구되지만 유한차분법으로 구한 민감도는 문제가 대형화될 때 계산비용의 제한때문에 비실용적이나 해석적 설계민감도는 효율적인 강점을 갖는다. 본 연구에서는 연성된 분자동역학-연속체 멀티스케일 문제에서 해석적 설계민감도를 유도하여 정확성과 향후 최적설계로의 활용 가능성을 확인하였다.

본 논문은 그래핀의 모드 I 균열 진전에 대한 분자동역학 해석과 수치보조장을 사용하는 영역 투영 방법의 역문제 해석방법을 결합하여 균열 선단 응집 법칙을 평가하는 효율적인 방법을 제시하고 있다. 그래핀의 균열 선단 응집 법칙을 결정하는 것은 균열 선단에서 멀리 떨어진 영역의 변위를 사용하여 균열 면에서 미지의 응집 트랙션과 열림 변위를 구하는 역문제를 해석해야 하는데 상호 J-적분과 M-적분의 경로 보존성과 효율적인 수치보조장을 사용하는 방법을 적용하였다. 분자동역학 해석에서 원자 변위를 유한요소 절점 변위로 이동최소자승법을 사용하여 근사하였으며 안정적인 역문제 해석을 통하여 원자 단위의 거동을 연속체 해석으로 연결시킬 수 있는 새로운 방법을 보여주었다.

Free surface를 가진 레나트-존스 이원계의 가열, 냉각에 따른 상전이를 molecular dynamics simulation을 이용하여 연구하였다. 용매-용질우너자간의 size misfit과 용질우너자의 농도변화에 따른 상전이 거동을 살펴본 결과 size misfit과 용질원자의 농도가 증가함에 따라 가열시 용융시 낮은 온도에서 쉽게 발생하였고, 냉각시 임계냉각속도가 감소하여 결정화는 어려워졌다.

정압분자동역학 시뮬레이션에 의하여 주기경계조건을 지닌 L-J 입자들로 구성된 계의 액체-유리-결정 전이를 연구하였다. 원자체적과 엔탈피는 가열 및 냉각과정에서 온도의 함수로 계산되었다. 반경분포함수로부터 유도된 Wendt-Abraham비와 단거리규칙도를 나타내는 각도분포함수를 분석하여 액체, 유리 및 결정상태를 구분하였다. 초기 fcc 결정을 가열하여 얻은 액체상은 급냉시에 비정질화하나 서냉시엔 결정화하였다. 급냉으로 생긴 유리는 다시 서서히 가열하면 fcc로 결정화하였다. 자유표면을 지닌결정은 표면에서부터 용해가 시작되어 벌크에 비하여 낮은 온도에서 녹고 냉각시에는 빠른 냉각속도에서도 결정화가 쉽게 일어나는 경향을 보였다.

버네사이트(birnessite)는 약 7Å의 d-spacing을 가지는 대표적인 층상형 산화망간광물로 높은 양이온 교환능력을 가지기 때문에 지하수나 퇴적물 공극 유체의 화학조성을 결정짓는 중요한 역할을 한다. 버네사이트의 양이온 교환 반응 기작을 규명하기 위해서는 층간 내 양이온의 배위 환경과 결정구조에 대한 원자 수준의 이해가 매우 중요하다. 이번 연구에서는 원자 수준의 계산광물학 방법인 고전 분자동역학(classical molecular dynamics; MD) 시뮬레이션을 이용하여 기존 실험에서 보고된 화학조성을 가지는 삼사정계 N a-와 K-버네사이트의 결정구조, 층간 양이온의 배위 환경 및 적층 구조를 계산하였다. 계산 결과는 기존 X-선 실험에서 보고된 격자 상수와 층간 배위 환경을 잘 재현하여 시뮬레이션 방법의 신뢰성을 보여주었으며, X-선 실험만으로는 구분하기 어려운 층간의 양이온과 물 분자 위치를 구별한 원자 수준의 정보를 제공하였다. 망간 팔면체 층의 적층 순서는 동일하지만 층간 내 N a+와 K+의 위치가 서로 상이하고, 층간 양이온의 배위 환 경과 결정구조 간의 상관관계를 보인다. 원자 수준의 분자동역학 시뮬레이션은 버네사이트의 양이온 교환 반응 기작 규명에 크게 기여할 것으로 기대한다.

란시아이트(ranciéite)는 수화된 Ca2+ 양이온이 망간 원자 빈자리를 아래위로 덮고 층간을 채우고 는 육방정계 층상형 산화망간광물(phyllomanganates)이다. 망간 원자 빈자리를 Mn2+ 양이온이 더 우세하게 채 우는 경우, 다카네라이트(takanelite)라는 광물로 구분하며, 란시아이트와 다카네라이트는 서로 고용체를 이룬 다. 이 광물들은 입자크기가 매우 작고 다른 광물과 함께 산출되기 때문에 실험만으로 정확한 결정구조를 규 명하기 어렵다. 이번 연구에서는 층간 Mn2+/Ca2+ 양이온 비율에 따른 란시아이트-다카네라이트의 결정구조와 층간 구조를 규명하기 위해 고전분자동역학 시뮬레이션(molecular dynamics simulations; MD)을 수행하였다. 연구방법의 적합성을 판단하기 위해 결정구조가 잘 알려진 칼코파나이트 군(chalcophanite group) 광물들에 대해 시뮬레이션 계산을 수행 후 실험 결과와 비교하였다. 이후 층간 양이온 비율에 따른 란시아이트 및 다카 네라이트 모델에 대한 MD 시뮬레이션을 수행하여 양이온 함량에 따른 양이온과 물 분자의 분포 및 (001) 면간거리를 제시한다.

점토광물은 지하수 바닥부터 산림에 이르는 지구의 얇은 표면에 해당하는 ‘크리티컬존(critical zone)’에 존재하는 금속원소의 지구화학적 순환에 중요한 역할을 한다. 입자 크기가 매우 작은 점토광 물에 대한 원자 수준(atomistic scale)의 연구는 지구화학적 순환 과정에 대한 정확한 기작(mechanism) 을 규명할 수 있을 뿐만 아니라 재료개발과 같은 산업분야에도 응용될 수 있다. 원자 간의 페어 퍼텐 셜(pair potential)을 파라미터화한 힘 장(force field)을 사용하는 분자동역학(molecular dynamics) 컴퓨 터 시뮬레이션은 원자 수준의 정보를 제공할 수 있기 때문에 실험과 함께 점토광물의 결정구조와 반 응도 연구에 사용된다. 점토광물 시뮬레이션을 위한 힘 장으로는 이팔면체(dioctahedral) 광물을 기반 으로 만들어진 ClayFF 힘 장이 보편적으로 사용된다. 삼팔면체(trioctahedral) 광물 시뮬레이션에도 ClayFF를 사용하는 연구가 보고되고 있으나, 같은 광물을 계산하더라도 각 연구마다 다른 파라미터 값을 사용하고 있기 때문에 파리미터 선택이 시뮬레이션의 정확도에 어떤 영향을 미치는지 체계적인 테스트가 필요하다. 이번 연구에서는 삼팔면체 광물인 수활석, 리자다이트, 활석을 대상으로 팔면체 마그네슘(Mg)의 원자간 페어 퍼텐셜을 나타내는 파라미터 ‘mgo’와 ‘mgh’를 각각 사용하여 분자동역 학 시뮬레이션 계산결과를 비교하였다. 격자상수, 원자 간의 거리 등 삼팔면체 점토광물의 결정구조는 주어진 두 가지 파라미터에 관계없이 거의 일정한 결과를 보여주었지만, 진동 파워 스펙트럼(vibrational power spectrum)으로 계산한 수산기의 진동수는 파라미터에 따라 상대적으로 뚜렷한 차이를 보였다.

점토광물은 대기부터 지하수에 이르는 크리티컬존(critical zone) 영역에서 금속과 탄소 순환을 결정짓는 역할을 한다. 계산광물학 연구방법 중에 하나인 분자동역학(molecular dynamics) 시뮬레이션 은 지구물질을 원자단위로 계산하기 때문에, 점토광물의 물리화학적 현상들에 대해 원자수준의 자세 한 정보를 제공할 수 있다. 이번 연구에서는 clayFF 힘 장(force field)을 사용한 분자동역학 시뮬레이 션을 이팔면체 점토광물인 깁사이트(gibbsite, Al(OH)3), 카올리나이트(kaolinite, Al2Si2O5(OH)4), 파이 로필라이트(pyrophyllite, Al2Si4O10(OH)2)에 적용하여 300 K, 1기압조건에서 각 광물이 가지는 격자상 수와 원자간 거리를 계산하고 실험결과와 비교하였다. 더불어 수산기의 방향성 및 수소결합의 양상 그리고 파워스펙트럼(power spectrum)을 추가적으로 계산하였다. 계산 결과, 격자상수는 기존의 실험 결과와 약 0.1~3.7% 미만의 오차율을 보였으며, 원자간 거리는 실험결과와 약 5% 미만의 차이를 가 졌다. 깁사이트나 카올리나이트의 팔면체층 표면에 존재하는 수산기가 가지는 신축진동파수(stretching vibrational wavenumber)는 실험값 보다 약 200-300 cm-1 높게 계산되었지만, 팔면체층 표면에 존재하 는 수산기들의 상대적 크기의 경향은 기존 실험 결과와 일치하였다. 팔면체층 표면의 수산기가 (001) 면과 이루는 각도도 기존 실험결과와 상당부분 일치한 반면에 내부 수산기의 경우는 다소 차이를 보 였다. ClayFF를 사용한 분자동역학 시뮬레이션 연구 방법은 이팔면체 점토광물 표면 내(층간) 금속이 온 흡착에 대한 수산기의 역할을 규명하는데 유용한 연구방법이 될 수 있음을 시사한다.

점토광물은 대기부터 지하수에 이르는 크리티컬존(critical zone) 영역에서 금속과 탄소 순환을 결정짓는 역할을 한다. 계산광물학 연구방법 중에 하나인 분자동역학(molecular dynamics) 시뮬레이션 은 지구물질을 원자단위로 계산하기 때문에, 점토광물의 물리화학적 현상들에 대해 원자수준의 자세 한 정보를 제공할 수 있다. 이번 연구에서는 clayFF 힘 장(force field)을 사용한 분자동역학 시뮬레이 션을 이팔면체 점토광물인 깁사이트(gibbsite, Al(OH)3), 카올리나이트(kaolinite, Al2Si2O5(OH)4), 파이 로필라이트(pyrophyllite, Al2Si4O10(OH)2)에 적용하여 300 K, 1기압조건에서 각 광물이 가지는 격자상 수와 원자간 거리를 계산하고 실험결과와 비교하였다. 더불어 수산기의 방향성 및 수소결합의 양상 그리고 파워스펙트럼(power spectrum)을 추가적으로 계산하였다. 계산 결과, 격자상수는 기존의 실험 결과와 약 0.1~3.7% 미만의 오차율을 보였으며, 원자간 거리는 실험결과와 약 5% 미만의 차이를 가 졌다. 깁사이트나 카올리나이트의 팔면체층 표면에 존재하는 수산기가 가지는 신축진동파수(stretching vibrational wavenumber)는 실험값 보다 약 200-300 cm-1 높게 계산되었지만, 팔면체층 표면에 존재하 는 수산기들의 상대적 크기의 경향은 기존 실험 결과와 일치하였다. 팔면체층 표면의 수산기가 (001) 면과 이루는 각도도 기존 실험결과와 상당부분 일치한 반면에 내부 수산기의 경우는 다소 차이를 보 였다. ClayFF를 사용한 분자동역학 시뮬레이션 연구 방법은 이팔면체 점토광물 표면 내(층간) 금속이 온 흡착에 대한 수산기의 역할을 규명하는데 유용한 연구방법이 될 수 있음을 시사한다.

핵폐기물을 고화시키는 재료로 사용하는 붕규산염(borosilicate) 유리의 용해는 지층 처분장에 처리된 고준위 방사성 폐기물의 생태계 유출을 결정할 수 있는 중요한 화학반응이다. 습식 실험에서 유리의 용해속도(dissolution rate)는 유리 화학조성에 의해 크게 좌우되는 것이 관찰된다. 유리의 bulk 구조를 규명한 분광분석 실험에 의하면 유리의 화학조성과 분자수준(molecular-level) 구조(예: SiO4 사면체의 연결구조와 B 원소의 배위구조) 사이의 상관관계가 존재한다. 따라서 화학조성에 따른 유리 용해도의 차이는 조성에 따른 bulk 내부구조의 변화로 이해되어 왔다. 그런데 유리 표면은 수용액과 계면을 이루면서 용해 과정에서 가장 직접적으로 반응하는 부분이기 때문에, 화학조성에 따른 표면구조 변화에 대한 지식 또한 필요하다. 본 논문에서는 분자 동역학(molecular dynamics, MD) 시뮬레이션을 사용하여 4가지의 다른 화학조성을 가지는 소듐붕규산염 유리(xNa2O·B2O3·ySiO2 화학조성)에 대하여 bulk 구조와 실험으로 얻기 어려운 표면(surface) 구조를 연구하였다. MD 시뮬레이션은 유리 표면의 화학조성과 분자수준 구조가 bulk의 것과 매우 상이한 결과를 보여준다. 본 연구의 MD 시뮬레이션 결과는 화학조성에 따른 유리 용해도(특히 초기 용해과정)는 bulk 구조의 변화보다 유리 표면구조의 변화에 의해 크게 좌우될 수 있다는 표면구조에 대한 이해의 중요성을 역설한다.