Journal Information
Journal ID (publisher-id): chemical
Title: 대한화학회지
Translated Title (en): Journal of the Korean Chemical Society
ISSN (print): 1017-2548
ISSN (electronic): 2234-8530
Publisher: 대한화학회Korean Chemical Society
다고리 방향족 탄화수소(polycyclic aromatic hydrocarbons, PAHs)는 인간의 활동 등에 의하여 방출되는 환경오염 물질 중의 하나이다. 이 PAHs는 돌연변이 유발성 및 발암성으로 인해 많은 연구가 이루어지고 있다.1-5 PAHs 중 페난트렌(phenanthrene, Phe)은 대기 중에 가장 많이 존재하는 PAH로 알려져 있다.6,7 이 Phe 기체는 대기 중에서 하이드록시 라디칼(·OH), 질산 라디칼(·NO3), 오존(O3)에 의하여 분해되는 것으로 알려져 있으며, Phe의 주된 분해 과정은 하이드록시 라디칼(·OH)에 의한 것으로 제시되고 있다.8-10
·OH의 방향족 고리에 대한 반응은 경쟁적인 두 가지 경로를 통해 일어나는 것으로 여겨지고 있는데,11-17 하나는 방향족 고리의 H 원자가 ·OH에 의하여 분리되어 H2O 분자를 생성하고, PAH는 라디칼이 되는 경로를 거쳐 일어나는 것이다.
그리고 다른 하나는 ·OH가 PAH에 첨가되어 하이드록실-PAH 라디칼을 거쳐 일어나는 것이다.
이러한 반응 경로는 반응속도론적 실험을 통해서 벤젠,11-13 나프탈렌,14,15 안트라센16,17의 경우 낮은 온도에서는 ·OH가 첨가되는 과정이 지배적이고, 높은 온도에서는 H가 분리되는 과정이 중요해지는 것으로 측정되었다. 그러나 PAHs 중 대기 환경 오염 물질로 가장 중요하다고 볼 수 있는 Phe에 대한 ·OH의 반응 경로에 실험적인 측정은 아직 이루어지지 않고 있다.
한편 실험적으로 측정한 IR 스펙트럼과 이론적인 계산으로 구한 IR 스펙트럼을 결합한 IR 스펙트럼 해석은 PAH들을 동정하는데 매우 유용한 기법으로 사용되어 오고 있다.18,19 그러므로 Phe에 대한 ·OH의 반응에서 나타나는 화학종들에 대하여 이론적인 계산을 통해 얻은 IR 스펙트럼은 반응 경로를 실험적으로 결정하는데 매우 중요한 자료 중 하나이다. 따라서 이 연구에서는 ·OH에 의하여 개시되는 Phe 분해 반응 과정을 동정하는데 필요한 모든 화학종에 대한 IR 스펙트럼을 DFT 계산을 사용하여 통하여 구하였고, 각 화학종들을 동정하는데 필요한 IR 스펙트럼의 특징을 고찰하였다.
Phe의 ·OH에 의한 반응에 의하여 생성되는 1차 생성물인 5종의 페난트렌올(phenanthrenols, Phe-ol)은 Fig. 1과 같다. 이 5종의 Phe-ol은 JiYi Lee 등에 의하여 기체 상태에서 측정되었다.20 그러므로 Phe의 ·OH에 의한 반응에서는 Fig. 1에 주어진 5종의 Phe-ol 이성체들이 모두 생성될 수 있음을 알 수 있다.
따라서 대기 중에서 ·OH에 의한 Phe의 분해 경로를 확인하는데 있어서는 Phe와 Phe-ol, 페난트레닐 라디칼(phenanthrenyl radical, Phe 라디칼), 하이드록실 페난트렌 라디칼(hydroxyl-phenanthrene radical, OH-Phe 라디칼) 각 5종에 대한 기체 IR 스펙트럼이 필요하다. 그러므로 이들 모든 화학종을 동정하는 데 필요한 IR 스펙트럼을 DFT 계산을 통하여 구하였다. 이 IR 스펙트럼들을 얻기 위한 모든 계산은 Gaussian 09 프로그램21을 사용하여 수행하였다. DFT 계산을 위한 각 화학종에 대한 초기 구조는 GaussView 5.0 프로그램22을 사용하여 그린 후 높은 대칭 점군으로 고정시켜 Gaussian 09 프로그램의 입력 파일을 작성하였다. 이 입력 파일을 사용하여 구조 최적화 계산은 PM6 → BLYP/6-31g(d,p) → BLYP/6-311++G(d,p)로 단계적으로 계산 수준을 높이며 수행하였다. Phe 라디칼 5종을 제외한 IR 스펙트럼은 BLYP/6-311++G(d,p) 계산에 의하여 얻어진 최적화 구조를 사용하여 진동수 계산(명령어 Freq)을 통해서 얻었다. 반면에 Phe 라디칼 5종의 경우에는 Gaussian 09 프로그램에서 BLYP/6-311++G(d,p) 수준에서 Freq 계산이 수행되지 않아 Opt=CalcAll 명령어를 사용한 계산을 통해서 얻었다.
Phe에서 각 탄소원자의 위치는 Fig. 2(a)에 나타낸 바와 같이 대한화학회에서 권장하는 방법에 따라 번호를 매겼다. 여기에서 Phe 분자는 C2v 대칭을 가지므로 C1, C2, C3, C4, C9는 각각 C8, C7, C6, C5, C10과 동등하고, C4a, C8a는 각각 C4b, C10a와 동등하다.
Phe는 Fig. 3과 같이 ·OH에 의하여 5개의 비동등 수소원자(C1–H, C2–H, C3–H, C4–H, C9–H)의 탈수소화 또는 5개의 비동등 탄소원자(C1, C2, C3, C4, C9)에 대한 ·OH 첨가에 의한 라디칼 생성 과정을 거쳐 5종의 Phe-ol이 생성될 수 있다. 이때 기존의 PAHs에 대한 실험적 연구11-17로부터 탈수소화에 의한 Phe 라디칼을 거치는 과정은 상대적으로 높은 온도에서 진행되고, ·OH 첨가에 의한 OH-Phe 라디칼을 거치는 과정은 상대적으로 낮은 온도에서 진행될 것으로 여겨진다. 이는 ·OH에 의한 Phe의 산화에 대한 이론적 연구에서 각 과정에 대한 활성화 에너지 Ea가 탈수소화 과정은 48.74~50.88 kJ mol−1, ·OH 첨가 과정은 32.09~43.44 kJ mol−1로 얻어지는 것으로부터도 뒷받침된다.23
Phe 라디칼은 Phe보다 단지 H 원자 하나가 적은 물질이므로 Phe 라디칼의 IR 스펙트럼은 Phe의 스펙트럼과 매우 유사한 특징을 가질 수 있다. 그러므로 반응 과정에서 Phe 라디칼의 생성 여부를 판단하기 위해서는 Phe 라디칼의 스펙트럼을 Phe에 의한 스펙트럼과 구별되어야 한다.
Phe와 Phe 라디칼의 IR 스펙트럼을 구하기 위하여 얻은 Phe의 최적화 구조는 Fig. 2(b), Phe에서 H가 하나 제거된 Phe 라디칼들에 대한 최적화 구조는 Fig. 4와 같다. 그리고 이들 최적화 구조를 가지고 진동수 계산에 의해 얻어진 IR 스펙트럼은 Fig. 5에 주어졌다.
Fig. 5의 IR 스펙트럼을 보면 CH들의 out-of-plane 굽힘 진동 영역인 650~850 cm−1에서 강한 세기의 피크들이 나타나고 있다. 이 CH out-of-plane 굽힘 진동 피크들은 Fig. 5에서 보는 바와 같이 Phe 라디칼의 피크들이 Phe의 것들과 뚜렷하게 다름을 알 수 있다.
Fig. 5의 650~850 cm−1 영역에서 Phe는 두 개의 매우 강한 피크(~720, ~800 cm−1)와 하나의 약한 피크(~850 cm−1)가 나타난다. 그러나 9-, 1-Phe 라디칼에서는 ~720 cm−1에서 매우 강한 피크와 비교적 강한 피크의 2개의 피크로 분리되어 나타남을 볼 수 있다. 그리고 2-Phe 라디칼은 ~730 cm−1에서 매우 강한 피크와 780~850 cm−1 영역에서 중간 세기의 3개의 피크가 나타나고 있다.
3-, 4-Phe 라디칼은 650~850 cm−1 영역에서 5개 이상의 피크가 나타나고 있다. 650~850 cm−1 영역에서 Phe는 ~720 cm−1에서 가장 강한 피크, ~810 cm−1에서 두 번째로 강한 피크가 나타나나 3-Phe 라디칼은 이들 피크는 ~810 cm−1에서 가장 강한 피크, ~720 cm−1에서 두 번째로 강한 피크가 나타나고 있다. 그리고 4-Phe 라디칼은 ~720 cm−1과 ~810 cm−1에서 서로 세기가 비슷한 세기의 피크가 나타나고 있음을 볼 수 있다.
따라서 IR 스펙트럼이 ~720 cm−1에서 두 개로 갈라진 강한 피크(9-, 1-Phe 라디칼), 780~850 cm−1에서 중간 세기의 3개의 피크(2-Phe 라디칼), 또는 Phe의 두 강한 세기의 피크 간에 세기 역전(3-, 4-Phe 라디칼)이 일어나면, Phe에서 H 원자가 제거되어 Phe 라디칼이 생성되는 경로의 반응이 일어나는 것을 확인할 수 있을 것이다.
OH-Phe 라디칼은 Phe-ol보다 H 원자가 하나 더 많은 물질이므로 OH-Phe 라디칼의 IR 스펙트럼은 Phe-ol의 스펙트럼과 매우 유사한 형태를 가질 수 있다. 그러므로 반응 과정에서 OH-Phe 라디칼의 생성 여부를 보기 위해서는 OH-Phe 라디칼의 스펙트럼을 Phe-ol에 의한 스펙트럼과 구분하여야 한다.
Phe-ol들의 최적화 구조는 Fig. 1과 같고, Phe에 ·OH가 첨가된 OH-Phe 라디칼들에 대한 최적화 구조는 Fig. 6과 같다. Phe-ol들과 OH-Phe 라디칼들에 대한 최적화된 구조를 가지고 진동수 계산에 의해 얻어진 IR 스펙트럼은 각각 Fig. 7과 8에 주어졌다.
Fig. 7의 Phe-ol들의 스펙트럼에는 CH 신축 진동 피크들이 ~3100 cm−1에만 나타나고 있다. 그러나 Fig. 8의 OH-Phe 라디칼들의 스펙트럼에는 OH가 결합된 탄소에 결합된 CH의 신축 진동 피크가 중간 정도의 세기로 ~3100 cm−1에 나타나는 피크들보다 ~200 cm−1 정도 적색 이동된 ~2900 cm−1에 나타나고 있다. 그리고 Fig. 7의 Phe-ol들의 스펙트럼들에서는 CH in-plane 굽힘 진동 영역인 1000~1300 cm−1에서 강한 세기의 피크들이 나타나고 있으나, Fig. 8의 Phe-ol들의 스펙트럼들에서는 CH out-of-plane 굽힘 진동 영역인 650~850 cm−1에서 강한 세기의 피크들이 나타나고 있다.
그러므로 IR 스펙트럼에서 다른 CH 신축 진동 피크들보다 ~200 cm−1 정도 적색 이동된 피크가 나타나고 강한 세기의 피크들이 650~850 cm−1에서 나타나면, Phe에 ·OH가 첨가되어 OH-Phe 라디칼이 생성되는 경로의 반응이 일어나는 것을 확인할 수 있을 것이다.
이 연구에서 저자들은 BLYP/6-311++G(d,p) 수준에서 Phe, 5종의 Phe 라디칼, 5종의 Phe-ol, 5종의 OH-Phe 라디칼의 진동수 계산을 통한 IR 스펙트럼을 구하였다. 이 결과로부터 페난트렌에 대한 두 가지 ·OH 반응 경로를 IR 분광법을 사용하여 동정하는 방법을 제시하였다.
고온에서 진행되는 것으로 알려진 Phe 라디칼 생성과정은 IR 스펙트럼의 CH out-of-plane 굽힘 진동 영역(650~850 cm−1)에서 Phe와 구분되는 피크들이 나타남을 볼 수 있었다. 그러므로 고온에서 Phe가 H 원자 제거 경로를 통해 진행되는 것은 IR 스펙트럼의 650~850 cm−1 영역에서 Phe의 스펙트럼과 구분되는 피크들이 나타나는지 여부로 확인할 수 있을 것이다.
반면에 저온에서 Phe에 ·OH 첨가되는 경로를 통해 진행되는 것은 Phe의 CH 신축 진동 피크 그룹에서 200 cm−1 정도 적색 이동한 영역에 CH 신축 진동에 의한 피크와 OH 신축 진동 피크(~3600 cm−1)가 나타나는 것으로 동정할 수 있을 것이다.
또한 ~3100 cm−1 영역에 CH 신축 진동 피크 그룹과 ~3600 cm−1에 OH 신축 진동 피크가 나타나고, CH in-plane 굽힘 진동 영역(1000~1300 cm−1)의 피크들의 세기가 CH out-of-plane 굽힘 진동 영역(650~850 cm−1)의 피크들보다 강하게 나타나는 것으로는 Phe와 ·OH 반응의 1차 생성물인 Phe-ol들이 생성된 것을 입증할 수 있을 것이다.
[1] Supplementary material Supporting Information. Additional supporting information is available in the online version of this article.
S. J. Wallace S. R. de Solla P. J. Thomas Harner T. A. Eng V. S. Langlois Ecotox. Environ. Safe.2018150176 [CrossRef]
T. Shimada Toxicol. Res.20173379 [CrossRef]
K. H. Lui B. A. M. Bandowe L. Tian C.-S. Chan J.-J. Cao Z. Ning S. C. Lee K. F. Ho Chemosphere2017169660 [CrossRef]
Y. Kamiya T. Kameda T. Ohura S. Thono Polycycl. Aromat. Comp.201737128 [CrossRef]
M.-J. Lee B.-D. Lee Tetrahedron Lett.2010513782 [CrossRef]
D. Helming W. P. Harger Sci. Total Environ.199414811 [CrossRef]
L. Wang R. Atkinson J. Arey Atmos. Environ.2007412025 [CrossRef]
Y. Zhang B. Yang J. Meng S. Gao X. Dong J. Shu Atmos. Environ.201044697 [CrossRef]
K. Lorenz R. Zellner Ber. Bunsenges. Physik. Chem.198387629 [CrossRef]
R. Atkinson J. Arey B. Zielinska S. M. Aschmann Environ. Sci. Technol.1987211014 [CrossRef]
I. V. Tokmakov M. C. Lin J. Phys. Chem. A200210611309 [CrossRef]
F. Goulay C. Rebrion-Rowe J. L. Le Garrec S. D. Le Picard A. Canosa B. R. Rowe J. Chem. Phys.2005122104308 [CrossRef]
R. Ananthula T. Yamada P. H. Taylor J. Phys. Chem. A20061103559 [CrossRef]
M. Tommasini A. Lucotti M. Alfe A. Ciajolo G. Zerbi Spectrochim. Acta A: Mol. Biomol. Spectrosc.2016152134 [CrossRef]
J. Lee D. A. Lane Atmos. Environ.2010442469 [CrossRef]