[논문 리뷰] Stress intrepretation of graphene E-2g and A-1g vibrational modes: theoretical analysis

D. Sfyris, Georgios I. Sfyris외 1명|arXiv (Cornell University)|2017. 06. 12.

Graphene research and applications참고 문헌 22인용 수 12

한 줄 요약

이 논문은 IKN 절차로부터 유도된 코시 응력 텐서를 사용하여 그래핀의 E₂g 및 A₁g 진동 모드와 관련된 응력 텐서 성분을 이론적으로 해석하기 위한 프레임워크를 개발한다. 이는 축 방향 응력과 비틀림 응력 성분이 이 모드들과 본질적으로 결합되어 있음을 드러내며, 45° 인장 하중에서 비틀림 응력(σ₁₂)이 지배적이 되어, 변형 분석에서 축 방향 우세성의 가정을 도전한다.

ABSTRACT

We here focus on only one graphene ring and examine to which stress tensor components the E2g and the A1g vibration mode of graphene correspond. These modes are typically related with the G-peak and the D-peak, respectively, and are strongly related to the stress distribution along the specimen. We adopt the theoretical framework of Admal and Tadmor ([1]) for the macroscopic definition of the Cauchy stress tensor and we introduce into this framework the E2g and the A1g as appropriate perturbations. We use these perturbations to the stress tensor expression and evaluate which stress tensor components are related to each vibrational mode. This approach, though qualitative in nature, incorporates all the main physics and reveals that E2g and A1g vibration modes should be related to shear as well as axial stress components when graphene is at rest (i.e. no external applied loading). To bring our framework closer to more concrete results, we evaluate the instantaneous Hardy stress tensor for a pair potential which correspond to the E2g and A1g modes at rest. Our analysis expands to take into account an applied external tensile field. Taking the armchair direction to be along the x-axis, when tension applies along the armchair direction, it is the axial $σ$11 stress component which dominates over $σ$12, $σ$22, which are of smaller order. When tension is along the zig-zag direction, it is the axial $σ$22 stress component that dominates over $σ$12, $σ$11. When tension is at an arbitrary direction between the armchair and the zig-zag direction, all stress components are of the same order and should all be taken into account even at small strains.

연구 동기 및 목표

• 그래핀의 E₂g 및 A₁g 진동 모드와 코시 응력 텐서 성분 간의 관계를 명확히 하기.
• 라만 스펙트로스코피 데이터를 2차원 물질의 매크로스코픽 응력 상태로 변환하는 데 도전하는 문제를 해결하기.
• 다양한 하중 방향이 축 응력과 비틀림 응력 성분의 상대 기여도에 어떻게 영향을 미치는지 조사하기.
• IKN 기반 응력 텐서 프레임워크를 그래핀의 진동 모드 변화를 포함하도록 확장하기.
• 하디 응력 텐서가 원자 척도 역학과 연속체 응력 장 사이의 관계에서 수행하는 역할 평가하기.

제안 방법

• 다체계에서 코시 응력 텐서를 정의하기 위해 IKN 절차를 채택하여, 저자들은 E₂g 및 A₁g 진동 모드를 외부 페르터베이션으로 포함하도록 이를 일반화한다.
• 응력 텐서는 동적 및 운동 에너지 성분으로 분해되며, 잠재 에너지의 연장에 의해 비유일성이 발생하지만, 공간 평균화에서는 이를 상쇄한다.
• 프레임워크는 정지 상태 및 외부 인장장 하에서 응력 텐서 표현에 E₂g 및 A₁g 모드를 조화 페르터베이션으로 포함한다.
• 원자 척도 역학과 연속체 응력 성분을 연결하기 위해, 쌍 위치 에너지 모델에 대해 하디 응력 텐서를 평가한다.
• 분석은 축방향(기준 축 x), 자이그재그(기준 축 y), 그리고 임의의 각도 θ에서의 인장 하중을 고려하며, θ = 45°는 중요한 사례로 간주된다.
• 편미분 이론과 삼각함수 전개를 사용하여, 적용된 인장 각도 θ에 따른 σ₁₁, σ₂₂, σ₁₂의 상대 크기를 평가한다.

실험 결과

연구 질문

• RQ1 정지 상태에서 그래핀의 E₂g 및 A₁g 진동 모드와 관련된 코시 응력 텐서 성분은 무엇인가?
• RQ2 다양한 방향에서 축 인장 하중을 가했을 때 축 응력과 비틀림 응력 성분의 상대 크기는 어떻게 변화하는가?
• RQ3 하디 응력 텐서는 그래핀에서 특정 진동 모드에 의해 유도된 응력 상태를 어느 정도 정확하게 캡처할 수 있는가?
• RQ4 결정학적 축으로부터 45° 각도로 하중을 가했을 때, 비틀림 응력(σ₁₂)은 라만으로 감지 가능한 변형에서 어떤 역할을 하는가?
• RQ5 특히 θ = 45° 근처에서 응력 텐서 성분은 적용된 변형과 하중 각도에 따라 어떻게 스케일링되는가?

주요 결과

• 정지 상태에서 E₂g 및 A₁g 진동 모드는 축 응력 외에도 비틀림 응력 성분을 유도하며, 이는 비틀림 변형과의 본질적 결합을 시사한다.
• 기준 축 방향(θ = 0°)으로 인장 하중을 가할 경우, 축 응력 성분 σ₁₁가 σ₁₂ 및 σ₂₂보다 지배적이다.
• 자이그재그 방향(θ = 90°)으로 인장 하중을 가할 경우, 축 응력 성분 σ₂₂가 σ₁₂ 및 σ₁₁보다 지배적이다.
• θ = 45°에서 비틀림 응력 성분 σ₁₂가 최대에 도달하며, 이는 축 성분 σ₁₁ 및 σ₂₂와 같은 주기의 크기를 가지게 되어, 소변형 영역에서도 무시할 수 없다.
• 임의의 인장 하중 각도에 대해, 특히 45° 근처에서는 축 응력과 비틀림 응력 성분 모두를 동시에 고려해야 하며, 이는 축 및 비틀림 항의 크기가 동등하게 스케일링되기 때문이다.
• 분석 결과, 라만 변형 맵핑에서 축 응력 성분만을 가정하는 것은 정확도를 떨어뜨릴 수 있으며, 특히 축 외부 하중일 경우 비틀림 기여를 포함시켜야 정량적 정확도를 확보할 수 있다.