과학, 중력 렌즈, 시뮬레이션

LensLab

점질량 플럭스 오차 < 10⁻⁶ (독립 적분 대비)

점질량과 특이 등온구(SIS) 렌즈를 다루는 얇은 렌즈 광선 추적 시뮬레이터입니다. 해석적 상 위치와 독립 플럭스 계산으로 검증했습니다.

직접 해 보기

데모는 앱의 원래 언어인 영어로 실행됩니다.

Critical curve (Einstein radius θE; for SIS it also marks the source-plane cut)Unlensed source (2σ)Analytic point-source imagesRuns LensLab's thin-lens equations in the browser; output checked pixel-by-pixel against a CPU reference. Brightness is display-stretched. Field of view 4″, θE = 1″.

b/θE = 0.95, θE = 1 arcsec인 SIS 렌즈의 세 패널: 점 코스틱과 유사 코스틱 옆의 가우스 광원, 임계 곡선과 해석적 +·− 상, 밝은 바깥 상과 희미한 안쪽 상이 보이는 렌더링한 렌즈 밝기. — SIS 시연, b/θE = 0.95 · 광원면, 해석적 점광원 상을 표시한 상면, 렌더링한 유한 광원 · SIS 시연 실행에서 프로젝트의 그림 스크립트로 생성.

질문

단순한 픽셀 광선 추적 렌더러가 렌즈 상과 플럭스를 얼마나 정확히 재현할까요? 그리고 그것을 어떻게 확인할 수 있을까요?

만든 것

점질량 모형과 특이 등온구(SIS) 모형이 함께 쓰는 렌즈 인터페이스를 만들었습니다. 렌즈 방정식, 굴절각, 야코비 행렬의 고윳값, 배율을 담고 있습니다.

해석적 점광원 해가 기준 위치를 줍니다. 타원 가우스 광원은 픽셀 중심에서 표본을 뽑아 표면 밝기를 보존합니다.

설정 파일로 움직이는 CLI가 픽셀 CSV, 그림, 출처 메타데이터를 만듭니다.

평가 방법

상 위치를 렌즈 방정식으로 되돌리면 광원 위치가 소수 12자리까지 재현됩니다. SIS에서는 두 상 사이의 거리가 소수 13자리까지 2θE와 같습니다.

유한 광원의 플럭스는 모델의 배율 함수를 전혀 쓰지 않는 독립 광원면 적분과 비교합니다.

SIS 격자 수렴은 Richardson 외삽으로 확인합니다.

주요 결과

시험한 오프셋 β/θE = 0, 0.5, 1에서 점질량 유한 광원 플럭스는 독립 적분과 상대 오차 10⁻⁶ 미만으로 일치합니다.

SIS 렌더러는 2차로 수렴하며(관측 차수 2 ± 0.2), N = 1024에서 전체 오차가 5×10⁻³ 미만입니다.

간단한 절단 밝기 추정치가 측정한 안쪽 상의 오차를 두 배 이내로 따라갑니다.

제 역할과 에이전트 활용

모델의 가정과 검증 전략을 검토한 뒤, 핵심 결과를 해석해와 독립 수치 계산에 대조했습니다.

Claude와 Codex가 범위가 정해진 조사와 구현을 돕습니다. 결과물을 검토하고, 어떤 주장을 프로젝트에 담을지는 제가 책임집니다.

한계

얇은 렌즈 하나와 원형 모형(점질량, SIS)만 다룹니다. 광원의 표면 밝기는 픽셀 중심에서 표본을 뽑습니다.

PSF, 잡음, 시간 지연, 다중 렌즈 평면은 모형에 없으며, 질량이나 거리를 추정하지 않습니다.

근거와 재현성

렌즈 방정식은 Narayan & Bartelmann의 Lectures on Gravitational Lensing(1996)을 따릅니다. 코드와 테스트 29개가 공개되어 있습니다.

Narayan & Bartelmann, Lectures on Gravitational Lensing (1996) ↗

모든 프로젝트로 돌아가기