[논문 리뷰] Extinction calculations of multi-sphere polycrystalline graphitic clusters - A comparison with the 2175 AA peak and between a rigorous solution and discrete-dipole approximations

Anja C. Andersen, J. A. Sotelo외 2명|ArXiv.org|2001. 11. 14.

Theoretical and Computational Physics참고 문헌 2인용 수 4

한 줄 요약

이 연구는 Gérardy & Ausloos의 엄밀한 다극자 방법과 두 가지 이산 두극자 근사법(DDSCAT 및 MarCODES)을 사용하여 다결정 구조의 그래피트 클러스터의 투과도를 계산하고, 분수형과 밀집형 형태를 비교한다. 작은 분수형 클러스터(N ≤ 50, 10 nm 입자)는 4.6 μm⁻¹에서 은하간 2175 Å 특징과 2% 이내의 투과도 피크를 생성하며, 반면 밀집형 클러스터는 4% 이상의 편차를 보여, 분수형 그래피트 클러스터가 2175 Å 피크의 강력한 후보임을 시사한다.

ABSTRACT

Certain dust particles in space are expected to appear as clusters of individual grains. The morphology of these clusters could be fractal or compact. In this paper we study the light scattering by compact and fractal polycrystalline graphitic clusters consisting of touching identical spheres. We compare three general methods for computing the extinction of the clusters in the wavelength range 0.1 - 100 micron, namely, a rigorous solution (Gerardy & Ausloos 1982) and two different discrete-dipole approximation methods -- MarCODES (Markel 1998) and DDSCAT (Draine & Flatau 1994). We consider clusters of N = 4, 7, 8, 27,32, 49, 108 and 343 particles of radii either 10 nm or 50 nm, arranged in three different geometries: open fractal (dimension D = 1.77), simple cubic and face-centred cubic. The rigorous solution shows that the extinction of the fractal clusters, with N < 50 and particle radii 10 nm, displays a peak within 2% of the location of the observed interstellar extinction peak at ~4.6 inverse micron; the smaller the cluster, the closer its peak gets to this value. By contrast, the peak in the extinction of the more compact clusters lie more than 4% from 4.6 inverse micron. At short wavelengths (0.1 - 0.5 micron), all the methods show that fractal clusters have markedly different extinction from those of non-fractal clusters. At wavelengths > 5 micron, the rigorous solution indicates that the extinction from fractal and compact clusters are of the same order of magnitude. It was only possible to compute fully converged results of the rigorous solution for the smaller clusters, due to computational limitations, however, we find that both discrete-dipole approximation methods overestimate the computed extinction of the smaller fractal clusters.

연구 동기 및 목표

• 다결정 구조의 그래피트 클러스터가 은하간 2175 Å 투과도 피크를 재현할 수 있는지 조사하기 위해.
• 엄밀한 전자기 해법(Gérardy & Ausloos)과 이산 두극자 근사법(DDSCAT 및 MarCODES) 간의 투과도 계산 정확도를 비교하기 위해.
• 클러스터 형태—분수형(D = 1.77) 대비 밀집형(단순 입방, 면심 입방)—이 투과도 특성에 미치는 영향을 평가하기 위해.
• 작은 및 중간 크기의 클러스터에 대해 DDA 방법의 수렴성과 신뢰성 평가하기 위해.
• 단일 결정 구조 근사 대비 다결정 구조가 2175 Å 특징의 피크 위치와 너비에 영향을 미치는지 평가하기 위해.

제안 방법

• 1982년 Gérardy & Ausloos의 엄밀한 다극자 이론을 사용하여, Draine & Lee(1984)의 그래프트 유전율 함수를 갖는 접촉하는 구형 클러스터의 투과도를 계산하였다.
• 두 가지 이산 두극자 근사(DDA) 코드인 DDSCAT(Draine & Flatau, 1994)와 MarCODES(Markel, 1998)를 사용하여 비교 분석을 수행하였다.
• 반지름 10 nm 또는 50 nm인 N = 4, 7, 8, 27, 32, 49, 108, 343개의 구를 포함하는 클러스터를 분수형(D = 1.77), 단순 입방, 면심 입방의 세 가지 기하형태로 구성하였다.
• 0.1–100 μm 범위에서 투과도 효율(Q_ext)을 계산하였으며, 작은 클러스터의 수렴성을 위해 최대 11차까지 다극자 수를 증가시켜 검증하였다.
• 투과도 피크의 위치와 반폭(FWHM)을 계산하여 관측된 2175 Å 피크(4.6 μm⁻¹)와 비교하였다.
• 비대칭적인 밀집형 클러스터의 방향성 의존성과 입자 크기 및 클러스터 크기에 대한 민감도를 평가하였다.

실험 결과

연구 질문

• RQ1 작은, 밀집형, 또는 분수형 다결정 구조의 그래피트 클러스터가 관측된 은하간 2175 Å 투과도 피크를 재현할 수 있는가?
• RQ2 0.1–100 μm 범위에서 분수형 클러스터의 투과도 특성은 밀집형 클러스터와 어떻게 비교되는가?
• RQ3 DDSCAT 및 MarCODES DDA 코드가 엄밀한 Gérardy & Ausloos 해법에 비해 이러한 클러스터의 투과도를 예측하는 데 얼마나 정확한가?
• RQ4 단일 결정 구조 근사(1/3–2/3) 대비 다결정 구조를 가정할 경우 2175 Å 특징의 피크 위치와 너비에 영향을 미치는가?
• RQ5 클러스터 크기(N)와 형태(분수형 대비 밀집형)가 투과도 피크의 안정성과 위치에 어떤 영향을 미치는가?

주요 결과

• 작은 분수형 클러스터(N ≤ 50, 10 nm 입자)는 관측된 4.6 μm⁻¹ 은하간 피크와 2% 이내의 투과도 피크를 생성하며, 더 작은 클러스터일수록 더 우수한 일치를 보였다.
• 밀집형 클러스터(단순 입방 및 면심 입방)는 4.6 μm⁻¹에서 4% 이상의 편차를 보이며 관측 결과와의 일치도가 열악함을 시사한다.
• 단파장 영역(0.1–0.5 μm)에서 분수형 클러스터는 밀집형 클러스터와 현저히 다른 투과도를 보이며 형태에 민감한 특성을 나타낸다.
• 파장이 5 μm를 초과하는 영역에서는 분수형 클러스터와 밀집형 클러스터의 투과도가 근사적으로 동일한 주기의 크기를 가지며, 장파장 영역에서 유사한 거동을 보임을 나타낸다.
• 두 DDA 방법(DRSCAT 및 MarCODES)은 엄밀한 Gérardy & Ausloos 해법에 비해 작은 분수형 클러스터의 투과도를 과대평가하며, 특히 장파장 영역에서 그러한 경향이 두드러진다.
• |m| > 2 인 파장에서 흡수를 과대평가하는 경향을 보이며, Draine & Goodman(1993)의 예측과 일치하며, 정확도는 두극자 수와 클러스터 크기에 따라 달라진다.