引力波天文学与理论模型引力波理论建模课件

数值相对论与引力波建模的一些新进展ZhoujianCao(曹周键)DepartmentofAstronomy,BNU2018-5-23第二届引力波天体物理学术研讨会数值相对论与引力波建模的一些新进展ZhoujianCaoContentCurrentstateofGWtheoreticalmodelFiniteelementmethodtoNRTheoreticalmodeltoeccentricbinarysystemSummaryandprospectContentCurrentstateofGWthe引力波探测数据是典型的强噪声、弱信号匹配滤波方法是引力波数据处理的理想方法引力波理论模型可以在既定硬件灵敏度的基础上很好地提高引力波探测的能力GW151226?13,5sigma引力波理论模型是波源参数反演的基础双黑洞?双中子星?星体质量、自旋？……引力波天文学与理论模型引力波探测数据是典型的强噪声、弱信号GW151226引力波理引力波理论建模针对实际的波源天体，求解爱因斯坦方程没有对称性，解析解基本没有希望在极端条件下，近似方法失效数值计算成为主要依赖的方法(数值相对论NR)

引力波理论建模针对实际的波源天体，求解爱因斯坦方程CurrentproblemsfacedbyNRCurrentproblemsfacedbyNRCurrentNRcodesCurrentNRcodesPDEformusedbyAMSS-NCKUPDEformusedbyAMSS-NCKUPDEformforFEPDEformforFEFEnumericalschemeCao,Fu,Ji,andXia,2018FEnumericalschemeCao,Fu,Ji引力波天文学与理论模型引力波理论建模课件Shockforms,butreallystable!ConsistenttoPDEpropertyofnon-linearlydegenerateMoreflexiblethanPSShockforms,butreallystableConvergencebehaviorbetterthanFD,4thvs2ndMoreaccuratethanFD,2ordersbetterConvergencebehaviorbetterthComputationcostofNR10orbits,256CPUs,7daysGW150914:4orbitsGW151226:28orbits3daysx250000=2055years(256CPUs)GW150914ComputationcostofNR10orbit引力波天文学与理论模型引力波理论建模课件StateofartforEOBNRTotalmass:0<M<∞Massratio:1<q<20GW150914(1.24),LVT151012(1.77),GW151226(1.89),GW170104(1.61),GW170608(1.71),GW170814(1.21)Spin:-1<s<0.9GW150914(<0.24),LVT151012(<0.3),GW151226(<0.35),GW170104(<0.42),GW170608(<0.3),GW170814(<0.18)Precessionsimplified(BCV)GW150914,LVT151012&GW151226assumenoprecession,GW170104,GW170608&GW170814gainlittleconstrainonprecessione=0AllannouncedBBHGWeventsassumedStateofartforEOBNRTotalmaLSC,ApJ818,L22(2016)LSC,ApJ818,L22(2016)BinaryBlackHolewitheccentricorbitPlausiblesourceforLIGOManysourcesforLISABinaryBlackHolewitheccentrSEOBNREmodelAssumesmalleccentricityexpansion(c)parteffectivelymorethan4PNElippartabout2PNCao&Han,PRD96,044028(2017)SEOBNREmodelAssumesmallecceNotlikeadiabaticapproximatedmodel,consistentcheckfore=0caseNotlikeadiabaticapproximateCompareto‘standardanswer’NRm1:m2=1:1e=0.02Compareto‘standardanswer’Nm1:m2=1:5e=0.02m1:m2=1:5m1:m2=1:1e=0.2m1:m2=1:1m1:m2=1:1000e=0.3a=0.9a=0.5a=-0.5a=-0.9m1:m2=1:1000a=0.9a=0.5aax-modelforeccentricBBHadiabaticapproximationmodelax-modelforeccentricBBHadiam1:m2=1:1e=0.3Moreaccuratethanadiabaticmodelforlaterinspiralandmergerm1:m2=1:1Moreaccuratethan初步发展有限元数值相对论算法，希望实现大质量比双星系统和中子星双星系统准确建模描述椭圆轨道双黑洞系统的引力波理论模型S