[논문 리뷰] User Guide for the Discrete Dipole Approximation Code DDSCAT 6.1

B. T. Draine, Piotr J. Flatau|ArXiv.org|2004. 09. 10.

ZnO doping and properties참고 문헌 1인용 수 31

한 줄 요약

이 논문은 이산 두극자 근사(DDA)를 사용하여 임의의 형상을 가진 입자 및 주기적 나노구조물의 전자기 산란과 흡수를 계산하는 무료로 이용 가능한 오픈소스 Fortran-90 코드인 DDSCAT 7.0을 제시한다. 이 소프트웨어는 크기 파라미터가 25 이하이고 굴절률이 1에서 멀리 떨어지지 않은 목표물에 대해 정확한 시뮬레이션을 가능하게 하며, MPI, OpenMP, MKL 및 DDfield 후처리기를 통한 근접장 E/B장 계산과 같은 고급 기능을 지원한다.

ABSTRACT

DDSCAT 6.1 is a software package which applies the discrete dipole approximation (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. DDSCAT 6.1 allows accurate calculations of electromagnetic scattering from targets with size parameters 2 pi a_eff/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 2). DDSCAT 6.1 includes support for MPI and FFTW. We also make available a "plain" distribution of DDSCAT 6.1 that does not include support for MPI, FFTW, or netCDF, but is much simpler to install than the full distribution. The DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions.

연구 동기 및 목표

• 이산 두극자 근사(DDA)를 구현한 오픈소스 소프트웨어 패키지인 DDSCAT 7.0에 대한 포괄적인 사용자 가이드를 제공하는 것.
• 크기 파라미터 ≤25인 경우에 대해 고정밀도로 임의의 형상을 가진 입자, 특히 비균일 및 이방성 재료의 산란 및 흡수를 모델링할 수 있도록 하는 것.
• 1차원 및 2차원 주기적 배열과 근접장 전자기장 계산을 지원하기 위해 MPI, OpenMP, Intel MKL와 같은 고급 계산 기능을 활용한 시뮬레이션을 가능하게 하는 것.
• 방향 평균화 및 장 후처리를 통해 산란 단면적, Mueller 행렬 원소, 복사력의 정확한 계산을 가능하게 하는 것.
• 사용자가 사용자 정의 기하학적 형상을 임포트하고 새로운 목표물 생성 루틴을 추가할 수 있도록 표준화되고 확장 가능한 프레임워크를 제공하는 것.

제안 방법

• 이산 두극자 근사(DDA)는 목표물을 다이폴 모멘트로 이루어진 격자로 이산화하여, 공액 기울기 방법과 같은 반복적 해법기를 사용해 결합된 다이폴 방정식을 푸는 방식이다.
• 코드는 반복적 해법기에서 효율적인 행렬-벡터 곱셈을 위해 빠른 푸리에 변환(FFT)을 사용하며, FFTW 및 Intel MKL와 같은 여러 FFT 라이브러리도 지원한다.
• 주기적 구조물의 경우, 장거리 다이폴-다이폴 상호작용을 처리하기 위해 Ewald 합산 기법을 사용하여 주기적 경계 조건(PBC)을 적용한다.
• 사용자는 매개변수 파일(ddscat.par)을 통해 목표물 기하학 및 구성 요소를 정의하며, 다이폴 위치, 재료, 정렬 방향, 입사파 매개변수 등을 지정한다.
• 코드는 원거리장 및 근접장 계산의 조합을 통해 산란 및 흡수 단면적, Mueller 행렬 원소, 복사력을 계산한다.
• DDfield 코드를 통한 후처리 기능을 통해, 목표물 근처의 사용자 정의 위치에서 전기장 및 자기장을 계산할 수 있으며, 이는 세밀한 장 분석을 가능하게 한다.

실험 결과

연구 질문

• RQ1 임의의 형상 입자 기하학을 위한 효율적이고 병렬 처리가 가능하며 확장 가능한 소프트웨어 패키지에서 이산 두극자 근사를 어떻게 효과적으로 구현할 수 있는가?
• RQ2 현대 HPC 아키텍처에서 대규모 다이폴 상호작용 시스템을 고정밀도와 높은 성능으로 해결하기 위한 최적의 계산 전략은 무엇인가?
• RQ3 DDA는 어떻게 주기적 나노구조 배열을 모델링하고 근접장 전자기 강화를 계산하는 데 확장될 수 있는가?
• RQ4 복잡한 굴절률을 가지며 이방성 재료를 갖는 목표물에 대해 DDA의 정확도와 적용 범위의 한계는 무엇인가?
• RQ5 사용자는 어떻게 신뢰성 있게 복잡한 목표물 기하학, 특히 다층 및 복합 구조를 시뮬레이션 프레임워크에 생성하고 임포트할 수 있는가?

주요 결과

• DDSCAT 7.0은 크기 파라미터가 25 이하이고 굴절률이 |m−1| ≤ 2를 초과하지 않는 목표물에 대해 정확한 전자기 산란 및 흡수 계산을 가능하게 한다.
• 코드는 MPI 및 OpenMP를 통한 병렬 실행을 지원하며, Intel MKL을 활용해 선형 대수 연산을 최적화하여 다중 코어 및 클러스터 시스템에서 성능을 크게 향상시킨다.
• DDfield 후처리기의 통합으로 사용자는 목표물 근처의 임의의 사용자 정의 점에서 전기장 및 자기장을 계산할 수 있어 세밀한 근접장 분석이 가능하다.
• 내장된 생성기와 사용자 정의 입력 파일을 통해 타원체, 실린더, 피라미드, 다층판, 복합 복합체 등 다양한 목표물 기하학을 지원한다.
• 매개변수 파일 ddscat.par는 입사파 편광, 산란 방향, 출력 옵션을 포함한 민감한 설정을 가능하게 하며, 무작위 정렬된 입자를 위한 방향 평균화 기능도 포함되어 있다.
• 코드는 해석적 해와 기존 결과와의 비교를 통해 검증되었으며, 알려진 문제들은 공식 DDSCAT 웹사이트에 문서화되고 업데이트되어 있다.