기본 콘텐츠로 건너뛰기

광속 이야기 (5)






 6. OPERA 중성미자 이상 (Neutrino anomaly)


 중성미자 (Neutrino) 는 보이지도 만질 수도 없고 무엇보다 지구 만한 물체도 그대로 투과해 버리는 도깨비 같은 입자로 우주에서 아주 흔한 입자로 생각되고 있습니다. 1930 년대 W. 파울리에 의해 제안된 이 입자는 경입자의 하나로 스핀은 1/2,  전하 0 인 입자입니다. 방사성 물질의 베타 붕괴때 에너지 보존 법칙이 설명되지 않는 이유를 설명하기 위해 도입한 중성미자는 이후 그 존재가 실제 확인되었습니다. 역사적으로 중성미자는 본래 물리학에서 중요한 입자였지만 특히 질량을 소량이라도 가지고 있는지 아니면 질량이 0 인지가 초미의 관심사였습니다. 


 만약 중성미자가 날아가는 동안 다른 종류로 변하는 중성미자 진동 (Neutrino Oscillation) 을 검출할 수 있다면 중성미자도 질량이 있다는 간접적인 증거가 될 수 있습니다. 따라서 이를 검출하기 위한 여러 실험이 진행되었습니다. 


 지금 부터 설명할 OPERA 실험 역시 마찬가지 입니다. OPERA 라는 단어는 Oscillation Project with Emulsion - tRacking Apparatus 의 약자로 중성미자 진동을 검출하기 위해 계획된 실험입니다. 중성미자는 원자로나 태양, 그리고 항성 및 초신성 등 여러군데에서 발생할 수 있지만 인공적으로도 만들 수 있습니다. 


 스위스 - 프랑스 국경 제네바 교외에 건설된 CERN 연구소에서 중성미자를 730 km 떨어진 이탈리아 중부 그란사소 (Gran Sasso) 에 있는 검출 장치로 발사하면 여기서 중성미자 진동을 발견하는 것이 실험의 주된 목표였습니다. 중성미자는 지구 만한 물질과도 거의 상호 작용없이 그대로 통과가 가능하기 때문에 지하로 통과하는 건 문제가 되지 않지만 검출하기는 역시 꽤 어려워 거대한 검출 장치가 그란 사소에 건설되었습니다. 




(CERN 지하의 뉴트리노 생성기에서 생성된 중성미자 (뮤온 중성미자) 가 지하를 통해 그란 사소로 날아간다  http://en.wikipedia.org/wiki/File:CNGS_layout_%28OPERA_experiment%29.jpg ) 



(CERN 에서 발사된 중성미자는 지구 지각을 통과 그란사소에 도달. )



 OPERA 연구진들은 2011 년 3월 데이터를 분석하는 과정에서 처음 이상한 결과를 발견합니다. 그것은 중성미자들이 대략 60.7 나노초 (nanosecond) 정도 빛보다 빨리 목표지점에 도달했다는 것입니다. 이후 6개월간 재검토와 검사가 이이진 끝에 그들은 9월 23일 자신들이 빛보다 더 빠른 뉴트리노를 검출했다고 보고했습니다. 


 물론 이 결과에 대해서 연구진들은 여러 과학자들의 겸허한 비판과 의견을 수용하기 위해서 이 결과를 숨기지 않고 공개했다고 했는데 솔직히 공개후에 쏟아질 많은 비판을 생각해 보면 이와 같은 자세는 당연한 것이었습니다. 왜 그런지는 차차 설명드립니다. 


 아무튼 이 OPERA 중성미자 이상 (OPERA neutrino anomaly) 은 곧 언론에 대서 특필 됩니다. 국내에서는 아인슈타인이 틀렸다. 혹은 타임 머신 가능 등으로 기사가 났는데 물론 이는 실험결과를 이해하지 못한 데서 나온 무지입니다. 그리고 진짜 빛보다 빠른 중성미자가 존재해도 역시 타임머신은 가능하지 않습니다. 이는 나중에 다시 설명합니다. 


 OPERA 연구팀은 실험 결과를 매우 소상히 공개했는데 우선 속도를 측정하는 방법부터 설명하면 이렇습니다. 속도를 측정하는 가장 기본적인 방법은 이동한 거리 / 걸린 시간 입니다. 지구상의 두 지점인 CERN 의 중성 미자 발생기와 그란사소의 검출기 사이의 거리와 중성미자가 날아간 시간을 검출하기 위해 GPS 에 기반한 매우 정밀한 측정이 이루어졌습니다.   



(중성미자가 날아간 시간을 측정하는 방식을 도표로 그린 것. 이해하기 힘든 그림이지만 아무튼 원자시계를 이용 시간 오차를 10억분의 1초 수준인 1 나노초로 맞추었다. 거리 오차는 20 cm 정도였다.  OPERA time measurement schematic. The PolaRx is the GPS receiver. The BCT, on the left, is associated with the proton beam, the proton spill, which creates neutrinos. This beam hits the "target" creating an intermediate particle, the muon, which then decays inside the "decay tunnel" into a neutrino. The straight line from the "decay tunnel" to "OPERA" actually passes through the solid rock of the earth's crust between CERN and LNGS. The waveform on the left is the measured proton distribution, and on the right, that of the detected neutrinos at OPERA.  ) 




 연구팀은 워낙 정밀한 시간 및 거리 측정이 이루어졌기 때문에 이 결과값에는 오류가 없을 것으로 자신했습니다. 사실 그렇지 않다면 엄청난 비판이 쏟아질 발표를 할 용기를 얻기 힘들었을 것입니다. 연구진에 속했던 과학자들은 발표에는 엄청난 용기가 필요했다고 설명했습니다. 


 연구진은 10월 21일 부터 11월 7 일까지 20개의 중성미자 검출 실험의 결과를 토대로 뉴트리노가 평균 62.1 나노초 광속보다 빨리 도달했다고 발표했습니다. 이는 물론 같은 거리를 진공에서 빛이 이동했을 때 보다 62.1 나노초 정도 더 빨리 목표에 도달했다는 의미입니다. 정규 분포 곡선으로 생각하면 6.2 sigma (표준 편차) 만큼 차이가 났기 때문에 이는 단순 오류로 생각하기 어려운 결과 였습니다. 


(연구진들이 결과를 종합해 제출한 Journal of High Energy Physics 에 2011 년 11월 17일 승인된 논문을 보려면 http://arxiv.org/ftp/arxiv/papers/1109/1109.4897.pdf  로 ) 


 OPERA 연구진들은 중성미자 발생기를 변경해 보는 등 다른 세팅에서도 검증 실험을 진행했지만 여전히 결과는 비슷했습니다. 하지만 이런 여러번의 검증 실험과 정밀한 관측에도 불구하고 이 결과는 역시 과학자 그룹에 비판에 직면했습니다. 


 그 이유는 간단했습니다. 이 실험 결과가 이전의 관측 결과와 완전히 상충되었기 때문입니다.    




댓글

이 블로그의 인기 게시물

통계 공부는 어떻게 하는 것이 좋을까?

 사실 저도 통계 전문가가 아니기 때문에 이런 주제로 글을 쓰기가 다소 애매하지만, 그래도 누군가에게 도움이 될 수 있다고 생각해서 글을 올려봅니다. 통계학, 특히 수학적인 의미에서의 통계학을 공부하게 되는 계기는 사람마다 다르긴 하겠지만, 아마도 비교적 흔하고 난감한 경우는 논문을 써야 하는 경우일 것입니다. 오늘날의 학문적 연구는 집단간 혹은 방법간의 차이가 있다는 것을 객관적으로 보여줘야 하는데, 그려면 불가피하게 통계적인 방법을 쓸 수 밖에 없게 됩니다. 이런 이유로 분야와 주제에 따라서는 아닌 경우도 있겠지만, 상당수 논문에서는 통계학이 들어가게 됩니다.   문제는 데이터를 처리하고 분석하는 방법을 익히는 데도 상당한 시간과 노력이 필요하다는 점입니다. 물론 대부분의 학과에서 통계 수업이 들어가기는 하지만, 그것만으로는 충분하지 않은 경우가 많습니다. 대학 학부 과정에서는 대부분 논문 제출이 필요없거나 필요하다고 해도 그렇게 높은 수준을 요구하지 않지만, 대학원 이상 과정에서는 SCI/SCIE 급 논문이 필요하게 되어 처음 논문을 작성하는 입장에서는 상당히 부담되는 상황에 놓이게 됩니다.  그리고 이후 논문을 계속해서 쓰게 될 경우 통계 문제는 항상 나를 따라다니면서 괴롭히게 될 것입니다.  사정이 이렇다보니 간혹 통계 공부를 어떻게 하는 것이 좋겠냐는 질문이 들어옵니다. 사실 저는 통계 전문가라고 하기에는 실력은 모자라지만, 대신 앞서서 삽질을 한 경험이 있기 때문에 몇 가지 조언을 해줄 수 있을 것 같습니다.  1. 입문자를 위한 책을 추천해달라  사실 예습을 위해서 미리 공부하는 것은 추천하지 않습니다. 기본적인 통계는 학과별로 다르지 않더라도 주로 쓰는 분석방법은 분야별로 상당한 차이가 있을 수 있어 결국은 자신이 주로 하는 부분을 잘 해야 하기 때문입니다. 그러기 위해서는 학과 커리큘럼에 들어있는 통계 수업을 듣는 것이 더 유리합니다. 잘 쓰지도 않을 방법을 열심히 공부하는 것은 아무래도 효율

R 스튜디오 설치 및 업데이트

 R을 설치한 후 기본으로 제공되는 R 콘솔창에서 코드를 입력해 작업을 수행할 수도 있지만, 보통은 그렇게 하기 보다는 가장 널리 사용되는 R 개발환경인 R 스튜디오가 널리 사용됩니다. 오픈 소스 무료 버전의 R 스튜디오는 누구나 설치가 가능하며 편리한 작업 환경을 제공하기 때문에 R을 위한 IDE에서 가장 널리 사용되어 있습니다. 아래 링크에서 다운로드 받습니다.    https://www.rstudio.com/  다운로드 R 이나 혹은 Powerful IDE for R로 들어가 일반 사용자 버전을 받습니다. 오픈 소스 버전과 상업용 버전, 그리고 데스크탑 버전과 서버 버전이 있는데, 일반적으로는 오픈 소스 버전에 데스크탑 버전을 다운로드 받습니다. 상업 버전의 경우 데스크탑 버전의 경우 년간 995달러, 서버 버전은 9995달러를 받고 여러 가지 기술 지원 및 자문을 해주는 기능이 있습니다.   데스크탑 버전을 설치하는 과정은 매우 쉽기 때문에 별도의 설명이 필요하지 않을 것 같습니다. 인스톨은 윈도우, 맥, 리눅스 (우분투/페도라)에 따라 설치 파일이 나뉘지만 설치가 어렵지는 않을 것입니다. 한 가지 주의할 점이라면 R은 사전에 반드시 따로 설치해야 한다는 점입니다. R 스튜디오만 단독 설치하면 아무것도 할 수 없습니다. 뭐 당연한 이야기죠.   설치된 R 스튜디오는 자동으로 업데이틀 체크하지 않습니다. 따라서 업데이트를 위해서는 R 스튜디오에서 Help 로 들어가 업데이트를 확인해야 합니다.     만약 업데이트 할 내용이 없다면 최신 버전이라고 알려줄 것이고 업데이트가 있다면 업데이트를 진행할 수 있도록 도와주게 됩니다. R의 업데이트와 R 스튜디오의 업데이트는 모두 개별적이며 앞서 설명했듯이 R 업데이트는 사실 기존 버전과 병행해서 새로운 버전을 새롭게 설치하는 것입니다. R 스튜디오는 실제로 업데이트가 이뤄지기 때문에 구버전을 지워줄 필요는

R 패키지 설치 및 업데이트 오류 (1)

 R 패키지를 설치하거나 업데이트 하다보면 여러 가지 문제가 생기는 경우들이 있습니다. 이 경우 아예 R을 재설치하는 것도 방법이지만, 어떤 경우에는 이렇게해도 해결이 안되고 계속해서 사용자는 괴롭히는 경우도 있습니다. 이런 경우 중 하나를 소개합니다.  새로운 패키지를 설치, 혹은 업데이트 하는 과정에서 같이 설치하는 패키지 중 하나가 설치가 되지 않는다는 메세지가 계속 나왔는데, 사실은 백신 프로그램 때문이었던 경우입니다.   dplyr 패키지를 업데이트 하려고 했는데, 제대로 되지 않아 다시 설치를 진행했습니다. 그런데 일부 패키지가 제대로 설치되지 않는다는 메세지가 나왔습니다.  > install.packages("dplyr") Error in install.packages : Updating loaded packages > install.packages("dplyr") Installing package into ‘C:/Users/jjy05_000/Documents/R/win-library/3.4’ (as ‘lib’ is unspecified) also installing the dependencies ‘bindr’, ‘bindrcpp’, ‘Rcpp’, ‘rlang’, ‘plogr’ trying URL ' https://cran.rstudio.com/bin/windows/contrib/3.4/bindr_0.1.1.zip ' Content type 'application/zip' length 15285 bytes (14 KB) downloaded 14 KB trying URL ' https://cran.rstudio.com/bin/windows/contrib/3.4/bindrcpp_0.2.2.zip ' Content type 'application/zip' length 620344 b