GW170817

Згідно з астрономічними спостереженнями й теоретичними моделями, центральна густина, маса та радіус нейтронних зір становлять 10¹⁵г/см³, 1,3 - 1,4 M_☉ і 10 - 15 км, відповідно, хоча типові значення цих фундаментальних величин досі невідомі[22][23]. Однією з найважливіших проблем є те, що до цих пір невідомий типовий радіус, оскільки: а) коректне визначення радіусу нейтронної зорі використовуючи лише астрономічні спостереження досить складно, а також б) досі не відоме рівняння стану для ядерної матерії з високою густиною як у нейтронній зорі, що ускладнює визначення відношення між її густиною і масою, тим самим вносячи невизначеність в типове значення радіусу [24][23][25].

Нейтронні зорі складаються головним чином із нейтронів, з часткою протонів, значно меншою (~11%), ніж нейтронів[26][27]. Це пов'язано з тим, що густина нейтронної зорі настільки висока, що енергія Фермі вироджених електронів значно перевершує різницю мас між нейтроном і протоном (m_n — _p)c² = 1,293 МеВ: енергія поверхні Фермі електронів ~100 МеВ[19]. Відтак, протягом гравітаційного колапсу масивного зоряного ядра нейтрони можуть утворитись через зворотній процес бета-розпаду вільних протонів та електронів, що сформують нейтрони:

${\displaystyle p+e^{-}\rightarrow n+v_{e}}$

допоки енергія Фермі нейтронів не стане такою ж високою, як електронів[28]. Тож, більша частина нейтронних зір (особливо масивних) складається із нейтронів, що перебувають у стані виродженого фермі-газу, а також невеликого домішку інших частинок. Від подальшого колапсу нейтронну зорю утримує тиск виродженого нейтронного газу. Менш масивні нейтронні зорі можуть складатись з кварків[28].

Більше половини з усіх зір формують подвійні пари, і чим масивнішою є зоря, тим більша ймовірність, що вона подвійна. Подвійні зорі еволюціюють разом і коли колапсують, то можуть залишити по собі нейтронні зорі, що будуть обертатись одна довкола одної й поступово наближатимуться по спіралі. Такі нейтронні зорі обертаються надзвичайно швидко, оскільки колапс масивної наднової до об'єкту в кілька десятків кілометрів збільшує кутову швидкість їх обертання внаслідок збереження кутового моменту[21].

Падіння по спіралі двох нейтронних зір може тривати десятки й сотні мільйонів років, воно поступово пришвидшується і врешті-решт відбувається їх майже миттєве злиття. Таке злиття компактних об'єктів призводить до викидання газу, збагаченого нейтронами. Речовина проходить через процес швидкого захоплення нейтронів (т.з. r-процес, від англ. rapid), утворюючи важкі елементи шляхом захоплення ядрами вільних нейтронів. Ці елементи нестабільні й зазнають радіоактивного розпаду, випромінюючи світло в оптичному та ближньому інфрачервоному діапазоні електромагнітного спектру[21]. Випромінювана енергія в 1000 разів більша, ніж нова зоря, звідси і назва «кілонова». Новоутворені важкі елементи (важчі заліза) викидаються в міжзоряний простір, накопичуються в газі та хмарах пилу, з яких потім утворюються нові зорі та планети. Таке злиття нейтронних зір найчастіше має відбуватись у старих галактиках, де зорі пройшли тривалу еволюцію — від масивних зір до наднових та нейтронних зір[28].

Система подвійних нейтронних зір формується із подвійної системи масивних зір, кожна з яких має масу ${\displaystyle M\gtrsim 8-10}$M_☉[29]. Коли перша з цих зір колапсує, утворюється нейтронна зоря, яка гравітаційно притягує матерію із зорі-компаньйона, прискорюючи її обертання. Коли друга зоря збільшиться і стане червоним гігантом, то нейтронну зорю поглине спільна оболонка, спричиняючи її рух по спіралі [30] Так формується дуже близька подвійна система зір. Вивільнена в цьому процесі (внаслідок акреції та тертя) енергія призведе до втрати гідрогенної оболонки, залишивши подвійну систему нейтронної зорі та гелієвої зорі[31]. Завдяки сильним припливним ефектами ця подвійна система матиме кругову орбіту.

Коли зоря-компаньйон вибухне надновою і утворить другу нейтронну зору, то вона втратить значну частку її маси в результаті викидів. Це призведе до відокремлення щільної подвійної системи, хіба що вибух не повідомить ударної швидкості новій нейтронній зорі[29]. Відтак, можливі два різні шляхи формування систем подвійних нейтронних зір: а) механізм з високою ударною швидкістю та викинутою масою, що пов'язаний із звичайною надновою і утворює системи з високими періодами обертання пульсару та високими ексцентриситетами орбіт; б) механізм з низькими ударними швидкостями та викидами, що пов'язані з захоплюючими електрони надновими, які формують системи з більш швидко обертаючимися пульсарами та меншими ексцентриситетами[32][30].

Після формування бінарної системи нейтронних зір, орбітальне розділення спадає протягом довго часу через випромінювання гравітаційних хвиль[33]. Ця фаза займає майже весь час існування системи подвійних нейтронних зір, за винятком останніх кількох мілісекунд їх життя. Із зменшенням орбітального розділення та наближенням двох нейтронних зір, вони швидко стають нестабільними. В результаті цієї динамічної нестабільності зорі починають швидко обертатись одна довкола одної, входячи в останню фазу злиття. Якщо маси двох нейтронних зір є майже однаковими, то злиття відбуватиметься як повільне зіткнення. Однак, у випадку, коли первинна зоря більш масивна, ніж друга, то остання зазнає припливного руйнування протягом зближення і врешті поєднається з первинною[34]. В цій фазі постають гравітаційні хвилі з максимальною амплітудою, які переносять інформацію про рівняння стану нейтронної зорі. Саме ж злиття може породити теплову енергію, яка врешті стане джерелом короткого гамма-спалаху. Гамма-спалах виникає тоді, коли нагріта ударною хвилею матерія генерує нейтрино і анти-нейтрино, які анігілюються довкола залишків злиття з утворенням високо-енергетичних фотонів[29].

Кінцевою фазою життя подвійної системи нейтронних зір є нова, динамічно стабільна конфігурація. Якщо маса залишку після злиття перевершує граничну масу, M_crit ~ 2,6 - 3,9 M_☉ (покриває, відповідно, діапазон «м'якого» та «жорсткого» рівнянь стану ядерної матерії) то він буде гравітаційно нестабільним і протягом декількох мілісекунд одразу колапсує, утворивши обертову чорну діру[30]. В іншому випадку, після злиття нейтронних зір залишок, залежно від його загальної маси, може пройти один з трьох можливих шляхів:[29]

1. У випадку, якщо загальна маса початкової подвійної системи нейтронних зір буде меншою, аніж максимальна маса необертової сферично симетричної нейтронної зорі, M_max (Ω = 0), то злиття такої подвійної маломасивної системи призведе до утворення стабільного залишку, який ніколи не колапсує в чорну діру;
2. Залишки з масами <M_crit будуть підтримуватись (протягом короткого часу) від колапсу через їх швидке обертання. Якщо маса залишку після злиття буде більшою від граничної маси для ізольованої системи, M_iso, але меншою від маси нейтронної зорі з однорідним обертанням (<1,2 M_iso), то може залишитись т.з. «надмасивна» нейтронна зоря, що буде підтримуватись обертанням твердого тіла і залишатись стабільною протягом кількох хвилин чи довше;
3. У випадку, якщо залишок матиме масу вище надмасивної межі, він може деякий час залишатись стабільним, як «гіпермасивна нейтронна зоря» (ГМНС), що утримується від гравітаційного колапсу швидким диференціальним обертанням. ГМНС може існувати лише від кількох десятків до сотень мілісекунд після злиття, проходячи серію різноманітних коливальних мод, а потім колапсує в обертову чорну діру внаслідок втрати власного диференціального обертання через гравітаційнохвильове випромінювання та гідромагнітні крутні моменти [35]. Енергія, що вивільняється протягом колапсу ГМНС може бути причиною затримки короткого гамма-спалаху. Тоді пік гравітаційно-хвильового випромінювання припадає на початковий етап злиття системи подвійностей, а випромінювання гамма-променів, що викликане колапсом ГМНС в чорну діру, відбувається значно пізніше[29].

У статті 1989 р. Ейхлер та ін. теоретично описали сценарій злиття двох нейтронних зір[36]. У цій статті було зроблено кілька оригінальних передбачень, більшість з яких було підтверджено з відкриттям GW170817 та електромагнітного двійника. Зокрема, Ейхлер та ін. передбачали, що злиття двох нейтронних зір

а) породить спалах гравітаційної радіації;

б) буде місцем нуклеосинтезу важких елементів;

в) джерелом гамма-спалаху

г) частота такого злиття становить 10 подій/рік у межах 100 мегапарсек.

Виділяють два загальні джерела викидів речовини до та після злиття нейтронних зір. Одне з джерел пов'язане з викидом матерії протягом динамічного часу (декілька мілісекунд) або припливними силами, або через викликане гравітаційним стисненням нагрівання на стику між двома зорями, що зливаються[37]. Тому такі викиди називаються «динамічними викидами», загальна маса яких для нейтронних зір, що зливаються, лежить в діапазоні 10^-4 - 10^-2 M_☉[38]. Другим джерелом викидів є витікання (протягом декількох секунд) із диску залишків після злиття: уламки, які не були одразу роз'єднані чи включені до центрального компактного об'єкту, можуть мати достатньо кутового моменту для циркуляції в акреційному диску довкола цього об'єкту[37]. Властивості викидів залежать від того, що саме відбудеться із масивним залишком після злиття нейтронних зір (стане надмасивною нейтронною зорею, гіпермасивною чи колапсує в чорну діру), що, в свою чергу, обумовлюється загальною масою початкової бінарної системи[39].

Теоретичні обчислення та моделювання припускають існування двох головних механізмів викиду матерії в злитті нейтронних зір[37][40]. Перший пов'язаний з гідродинамічними силами, які стискають матерію на стику між двома зорями і вона виштовхується квазі-радіальними пульсаціями залишку, врешті, будучи нагріта ударними хвилями, викидається в різних напрямках[41]. До другого механізму викидів залучені спіральні рукави, утворені припливними взаємодіями протягом злиття двох зір. Внаслідок в'язкого переносу кутового моменту в спіральних рукавах, частка матерії розширяється назовні в екваторіальній площині[38].

Результатом злиття нейтронних зір є утворення акреційного диску довкола центрального залишку. За підрахунками, типова маса диску складатиме ~0,01 - 0,3 M_☉[42]. Викиди матерії вітрами диску залишків (протягом секунд і довше) є другим джерелом викидів, що може конкурувати чи навіть домінувати над динамічними викидами [40]. Маса, викинута вітрами з диску може бути однаковою чи навіть більшою, аніж маса в динамічних викидах.[43]

Одразу після формування диску, він має високу швидкість акреції, будучи рясним джерелом теплових нейтрино, нагрівання яких приводить до втрати маси поверхнею диска [44]. Моделі торів залишків передбачають, що у випадку швидкого формування чорної діри після злиття, маса викинутої матерії із диску буде незначною, через швидке убування світності нейтрино[45][46]. У випадку, якщо залишки після злиття деякий час (довше ~50 мсек) проіснують, як гіпермасивна чи надмасивна нейтронна зоря, то більш велика нейтринна світність із такого компактного залишку викине значну масу речовини, ~10^-3 M_☉[47][48]. Двовимірні гідродинамічні моделі еволюції диску показали, що у випадку швидкого утворення чорної діри після злиття, частка електронів у нейтринних вітрах із диску коливається в діапазоні Y_e ~ 0,2 - 0,4[45]. Цього достатньо для утворення цілого асортименту ядер r-процесу[43].

Загальна маса акреційного диску, вивільнена вітрами під дією в'язких сил, коливається в діапазоні від ~5% (для чорної діри з повільним обертанням) до 30% (для чорної діри зі швидким обертанням)[46]. Проте, якщо центральним залишком після злиття є відносно довготривала гіпермасивна чи надмасивна нейтронна зоря, то завдяки наявності твердої поверхні та вищого рівня нейтринного випромінювання від такого залишку, може вивільнитись значно більша частка маси акреційного диску - до ~90%[49]. Крім того, електронна частка, Y_e викидів монотонно зростає з часом життя гіпермасивної нейтронної зорі, і якщо після злиття залишок нейтронної зорі проіснує довше 300 мс, то більшість вітрових викидів з диску не будуть містити лантаноїдів (Y_e ≥ 0,3)[50].

Окрім динамічних викидів та викидів дисковими вітрами, із довготривалого залишку в процесі його т.з. «скорочення Кельвіна-Гельмгольца» до кінцевого холодного стану, передбачаються додаткові витікання матерії, викликані нейтрино чи магнітним полем[51]. Такі витікання можуть бути додатковим джерелом викидів після злиття і мати значний ефект, якщо залишком стане надмасивна нейтронна зоря.

Наявні у всесвіті гідроген та гелій утворились під час Великого вибуху 13,8 млрд років тому. Важчі елементи, як кисень та вуглець, утворились в ядрах зір внаслідок термоядерного синтезу гелію та гідрогену[52].

Однак, для утворення елементів важчих заліза (наприклад, лантаноїдів), потрібні особливі умови, коли ядра атомів бомбардуються вільними нейтронами[53].

Ядро атома гідрогену, ¹H, складається з одного протона. Ядра ж усіх інших елементів, включаючи важчі ізотопи гідрогену, містять як протони, так і нейтрони (разом відомі, як «нуклони»). Ці нуклони утримуються разом сильною ядерною взаємодією. Маса спокою такого ядра є меншою суми мас спокою вільних протонів та нейтронів. Тож дефект маси перетворюється на енергію, відому як енергія зв'язку ядра. Більша енергія зв'язку (у розрахунку на один нуклон) означає, що нуклони зв'язані в ядрі сильніше. Енергія зв'язку значно зростає від протона до ядра вуглецю, далі поступово росте до заліза, де сягає максимуму, пояснюючи істотну поширеність заліза у всесвіті. Хоча синтез легких елементів зазвичай вивільняє енергію, проте електромагнітне (кулонівське) відштовхування між ядрами перешкоджає наближенню ядер на малу відстань одне до одного, на якій можлива сильна взаємодія. Необхідна для подолання сили електромагнітного відштовхування енергія називається кулонівським бар'єром[52].

Після заліза енергія зв'язку (у розрахунку на один нуклон) зменшується зі збільшенням атомного номера, тому нуклеосинтез елементів, важчих заліза, ускладнений, бо двом позитивно зарядженим ядрам треба подолати кулонівське відштовхування щоб наблизитись достатньо близько для синтезу. Саме тому елементи, важчі від заліза, утворюються шляхом захоплення нейтронів , які не мають заряду, а відтак, можуть наближатися до заряджених ядер, захоплюватися ними та збільшувати масу ядра. Однак із захопленням нейтрона в ядро, те стає радіоактивним і зазнає бета-розпаду. Це пояснює те, чому процес утворення важких елементів є надзвичайно складним [54].

Елементи до заліза мають приблизно однакову кількість нейтронів та протонів. Але після заліза кількість нейтронів перевищує кількість протонів і стабільність ядра досягається лише за рахунок значного надлишку нейтронів. Існують два процеси захоплення нейтронів, що призводять до формування нових елементів. Перший - повільне захоплення нейтронів (так званий s-процес (від англ. slow)) полягає в додаванні нейтронів до ядра по одному. Якщо утворюється нестабільне ядро, воно зазвичай зазнає бета-розпаду, перш ніж захопить наступний нейтрон. В s-процесі швидкість захоплення нейтронів є меншою, ніж бета-розпад, й утворюються лише порівняно стабільні ядра. S-процес відбувається в зорях масою 0,6 - 10 маси Сонця й зупиняється на наймасивніших стабільних ядрах бісмуту[54]. Для утворення більш важких стабільних елементів потрібне захоплення нейтронів швидше від бета-розпаду.

Довгий час загальнопоширеним поміж астрономів було уявлення про наднові, як панівне місце r-процесу. Зокрема, вважалось, що r-процес нуклеосинтезу відбувається у високо-ентропійних, нейтринних вітрах з прото-нейтронних зір, які утворились протягом кількох секунд після вибуху наднових типу ІІ[55][56][57][58]. Втім, таке припущення мало ряд теоретичних недоліків, а нещодавні результати чисельних моделювань показали, що нейтринні вітри позбавлені необхідних фізичних умов для пояснення поширеності всіх важких елементів в галактиці[59].

Накопичені протягом останніх років дані астрономічних спостережень та чисельні моделювання вказують, що джерелом r-процесу має буде щось більш рідкісне, аніж наднова[60][61][62][63].

1957 р. у статті, відомій як B²FH чотири фізики передбачили і пояснили механізм швидкого захоплення нейтронів[64]. Цей механізм вони назвали r-процесом (від англ. rapid). Захоплення нейтронів в r-процесі робить ядра надзвичайно радіоактивними, що потребує надзвичайно швидкого захоплення. У цій публікації припускалося, що r-процес має відбуватись у середовищах із надзвичайно великою густиною нейтронів, де кілька нейтронів може бути захоплено ядром підряд, до того, як нестабільне ядро зазнає бета-розпаду. У такий спосіб утворюються важкі елементи.

Того ж 1957 р. канадсько-американський фізик Аластар Кемерон незалежно від B2FH висунув припущення, що близько половини важчих від заліза елементів у галактиці генеруються лише в середовищах з такою високою густиною нейтронів, що захоплення нейтронів атомними ядрами відбувається значно швидше, аніж бета-розпад[65].

1974 р. Джеймс Латтімер та Девід Шрамм вперше висунули гіпотезу про злиття нейтронної зорі та чорної діри як джерело викидів матерії, багатої вільними нейтронами, що сприятиме r-процесу з дуже низькою часткою електронів, Y_e (відношення протонів до нуклеонів (протонів + нейтронів): ${\displaystyle Y_{e}=n_{p}/(n_{n}+n_{p})}$[66]. 1982 р. Симбалісти та Шрамм вперше висунули ідею, що схожий механізм викиду нейтронної матерії може відбуватись внаслідок злиття подвійної системи нейтронних зір, що і є джерелом r-процесу[67].

На основі перших чисельних моделювань злиття двох нейтронних зір було показано, що результатом злиття подвійної системи має бути припливний викид нейтронно-багатої матерії (${\displaystyle Y_{e}\lesssim 0.1-0.2}$) в орбітальній площині зі швидкістю ~0,2-0,3 c і масою ~10^-4 - 10^-2 M_☉[68][69]. Наслідком такого викиду буде поява важких елементів, передбачувана поширеність яких в цілому узгоджується із спостережуваною поширеністю в сонячній системі[70][71].

Пізніші чисельні моделювання показали, що окрім припливного викиду матерії, на межі контакту між двома зорями що зливаються, повинне формуватись окреме джерело викиду в полярному напрямку - викиди матерії, нагрітої ударними хвилями[72][38]. Нагрівання ударними хвилями та опроміненням нейтрино сприяє слабким взаємодіям, завдяки чому частка електронів динамічного полярного викиду є значно більшою (${\displaystyle Y_{e}\gtrsim 0.25}$) від його первинного значення всередині нейтронної зорі[73][74][75]. Ударно нагріті викиди менш багаті нейтронами, аніж викиди від припливних взаємодій, оскільки їх вища температура (>1 МеВ) уможливлює народження електрон-позитронних пар, що захоплюються відповідно протонами та нейтронами, випромінюючи електронні нейтрино та антинейтрино.[53]. За рахунок того, що частка нейтронів перевершує частку протонів, то в таких викидах значно більше захоплень позитронів, залишаючи матерію з більшою часткою електронів[37].

Ударно нагріті викиди відіграють важливішу роль для «м'якого рівняння стану» (нейтронні зорі менших радіусів), тоді як припливні викиди є панівними в злитті нейтронних зір з більшими асиметріями мас [37]. В останньому випадку злиття буде менш руйнівним, оскільки нейтронна зоря з меншою масою зазнає деформації припливними силами ще на ранній стадії злиття, тож ударне нагрівання буде менш сильним. Наслідком цього будуть більш збагачені на нейтрони припливні викиди[53].

Іншим джерелом викиду матерії r-процесу результати моделювання передбачають уламки від злиття нейтронних зір [76]. Ці уламки мають достатньо кутового моменту для циркуляції в акреційний диск довкола центрального залишку. Цей диск залишків може бути джерелом ультрарелятивістського джету гамма-спалаху[77]. Повільні відтоки матерії (із швидкістю ~0,03 - 0,1 с) з такого диску залишків, що здатні тривати протягом декількох секунд після злиття, можуть бути новим джерелом викиду речовини r-процесу[76][51]. Середнє значення Y_e відтоку акреційного диску зростає із збільшенням часу існування гіпермасивної нейтронної зорі перед тим, як вона колапсує в чорну діру[49].

Динамічні викиди із злиття нейтронних зір мають великі передбачувані маси, в діапазоні від 10^-3 до 10^-2 M_☉[78]. Тож, викинута в злитті двох нейтронних зір матерія має потрібну густину нейтронів та швидкість для запуску ядерних реакцій r-процесу [46].

Обчислення на основі астрономічних спостережень передбачають, що події злиття нейтронних зір є в 100-1000 разів рідшими, аніж гравітаційний колапс наднових типу ІІ[63]. Так, карликові галактики "забруднені" лише кількома подіями злиття бінарних систем нейтронних зір[79][80], чи взагалі однією подією, як карликова галактика Reticulum II[81]. В Reticulum II велика частка зір є високозбагачена елементами r-процесу[82]. Це свідчить про те, що ця галактика була "забруднена" ще на початку її історії єдиною подією r-процесу, яка породила більше важких елементів, аніж здатен нейтринний вітер однієї наднової[83].

Рідкісність злиття нейтронних зір узгоджується з обчисленнями неоднорідності хімічної еволюції, що дозволяє прослідкувати локальні коливання поширеності елементів r-процесу, зумовлені внеском одиничних подій злиття. Так, спостерігаючи хімічну поширеність в близьких тьмяних карликових сфероїдних галактиках, астрономам вдалось встановити джерело r-процесу через розрізнення внеску подій наднових типу ІІ та злиття нейтронних подвійностей[84][85]. Оскільки злиття нейтронних зір є більш рідкісним, ніж вибух наднових типу ІІ, то менш масивні тьмяні карликові сфероїдні галактики не можуть залишити багато нейтронних зір після вибуху масивних зір для формування подвійної системи і їх злиття. Обчислення припускають, що в менш масивних карликових сфероїдних галактиках з масою 10⁵ M_☉, в цілому буде ~500 наднових типу ІІ. Виходячи з отриманої частоти злиття нейтронних зір (1 подія злиття на 1000-2000 подій наднових типу ІІ) на таку галактику припадатиме <1 події злиття нейтронних зір[84]. Тому в менш масивних карликових сфероїдних галактиках (як галактики Дракон, Каріна і Скульптор) було виявлено стале значення [Eu/H] ~ −1,3 безвідносно до значень [ Fe/H ][84]. Тоді як в масивних карликових сфероїдних галактиках було виявлено збільшення [ Eu/H ] із збільшенням [Fe/H]. Це свідчить про те, що r-процес нуклеосинтезу не відбувається в тьмяних, менш масивних карликових сфероїдних галактиках, попри високу частоту подій наднових в цих галактиках[84]. Тож, ці дані є прямим свідченням, що лише такі рідкісні події, як злиття нейтронних зір можуть бути головним місцем r-процесу в галактиці, особливо для нуклідів з A > 130[85].

Результати чисельних моделювань показують, що нуклеосинтез у надзвичайно багатих нейтронами викидах із злиття нейтронних зір може відтворити поширеність важких елементів в сонячній системі незалежно від фізичних властивостей (зоряних мас, співвідношення мас та рівняння стану) подвійної системи[86]. Так, нещодавно були отримані переконливі свідчення про значні події r-процесу відносно нещодавно в Чумацькому шляху, зокрема в Сонячній системі. Таким свідченням нещодавнього збагачення подіями r-процесу стали довгоживучі радіоактивні елементи. Для вимірювання слідів таких подій в глибоководних відкладах на Землі були використані два ізотопи - заліза⁶⁰Fe та плутонію ²²⁴Pu. Ізотоп ⁶⁰Fe утворюється протягом еволюції та вибуху масивних зір, що призводить до появи наднових[87]. Період напіврозпаду ⁶⁰Fe складає 2,6 x 10⁶ років і тому цей ізотоп може бути свідченням недавнього додавання матерії із подій r-процесу, що відбулись кілька мільйонів років тому. Нещодавно ізотоп ⁶⁰Fe був знайдений в земних глибоководних відкладах, що включали зоряні залишки із відносно близького вибуху, що мав місце 2 млн років тому [88][89].

²²⁴Pu має період напіврозпаду 8,1 x 10⁷ років і потребує багатьох подій вибуху наднових. Виявлена в глибоководних відкладах поширеність ²²⁴Pu на два порядки нижче, аніж передбачається у випадку, коли б джерелом були такі часті події, як звичайні наднові з невеликим внеском[60]. З цього відкриття випливало, що нуклеосинтез актиноїдів є дуже рідкісним (уможливлюючи значний розпад ²²⁴Pu з часу останньої події r-процесу) і протягом останніх кількох сотень мільйонів років регулярні наднові не мали значного внеску до їх поширеності в сонячній системі[53]. Подальший аналіз цих експериментальних даних та теоретичні розрахунки вказували на те, що джерелом елементів r-процесу має бути значно більш рідкісна подія, аніж вибух наднових[61]. Втім, такі рідкісні події, як злиття нейтронних зір може пояснити як існування ²²⁴Pu в ранній Сонячній системі, так і низьку поширеність відносно нещодавніх внесків ізотопів ²²⁴Pu в глибоководних відкладах[61].

Однак, на відміну від наднових, злиття нейтронних зір не конденсує зерна пилу, як носіїв синтезованих нуклідів r-процесу (досонячні зерна)[90]. Тому як саме синтезовані в злитті нейтронних зір елементи r-процесу були включені до Сонячної системи наразі залишається незрозумілим.

Схематичне подання ймовірних електромагнітних двійників злиття двох нейтронних зір та кілонова

Щойно синтезовані у викидах матерії із злиття нейтронних зір атомні ядра r-процесу є радіоактивними[69][78]. Із подальшим розширенням цієї матерії відбувається її бета-розпад назад до стабільного стану. Енергія виділена через бета-розпад та поділ ядра може породити і підтримувати тепловий транзієнт - «кілонову», що триватиме від днів до тижнів[91][92].

Кілонова є ключовим електромагнітним двійником злиття нейтронних зір, оскільки порівняно із спрямованим джетом гамма-спалаху, її випромінювання майже ізотропне і може досягнути піку в оптичному діапазоні спектру, будучи відносно легко доступна для пошуку. Яскравість, кольори та тривалість кілонових є свідченням фізичних процесів, що відбуваються під час злиття подвійних нейтронних зір. Крім того, кілонові дозволяють прямо спостерігати та виміряти утворення ядер r-процесу в злитті, відкриваючи унікальну можливість встановити місце нуклеосинтезу важких елементів[92].

1998 р. Л.-С. Лі та Б.Пачінські вперше висунули припущення, що радіоактивний викид із злиття двох нейтронних зір (чи нейтронної зорі та чорної діри) може бути джерелом живлення теплового транзієнта, за аналогією з надновими типу Іа[93]. На основі побудованої простої моделі, вони передбачали, що внаслідок малої маси та високої швидкості (~0,1 c) викиду матерії із злиття нейтронних зір, цей викид швидко стане прозорим для власного фотонного випромінювання[93]. Таке випромінювання повинне досягнути піку протягом одного дня, що значно швидше ніж для регулярної наднової, чиї випромінювання досягають піку протягом тижня і довше.

Не маючи фізичної моделі нуклеосинтезу, Лі та Пачінські параметризували швидкість радіоактивного нагрівання (розпад ядер r-процесу) викиду ${\displaystyle {\dot {Q}}}$ в момент часу t після злиття, як ${\displaystyle {\dot {Q}}\propto t^{-1}}$, залишивши нормування швидкості нагрівання як вільний параметр[93]. Однак, оскільки пікова світність транзієнта пропорційна швидкості нагрівання, то модель Лі-Пачінські передбачала надзвичайно високі значення пікової світності, в діапазоні 10⁴² - 10⁴⁴ ерг сек[93], що більше від світності найяскравішої наднової [92]. Втім подальші астрономічні спостереження, які намагались виявити подібні яскраві транзієнти після коротких гамма-спалахів, таких не виявили[94][95].

2010 р. Б.Метцгер та ін. вперше визначили реальний масштаб світності радіоактивно-підтримуваних транзієнтів злиття нейтронних зір[91]. На основі обчислення мережі ядерних реакцій r-процесу, Метцгер та ін. отримали швидкість радіоактивного нагрівання (передбачаючи ${\displaystyle {\dot {Q}}\propto t^{-1,3}}$ від кількох годин до днів), яку включили в моделі кривих блиску. Дослідники використали більш фізично реалістичну модель непрозорості викидів r-процесу (непрозорість постає із змішання десятків мільйонів атомних спектральних ліній зв'язано-зв'язаних переходів - переходів електронів в атомі з одного енергетичного рівня на інший[96]). В основі цієї моделі було припущення, що непрозорість викидів r-процесу обумовлена непрозорістю в лініях заліза, а не як в простіших моделях - непрозорістю електронного розсіювання. Відтак, для викидів масою 1^-2 M_☉ і швидкістю v ~ 0,1 c, модель Метцгера та ін. передбачала пікову світність в ~3 x 10⁴¹ ерг/сек та спектральний пік у видимому діапазоні[91]. Оскільки така світність приблизно в тисячу разів яскравіша аніж класичні нові (чия пікова світність близька до критичної світності, ~10³⁸ ерг/сек), то електромагнітні транзієнти злиття нейтронних зір, спричинені розпадом ядер r-процесу, були названі Метцгером та ін. «кілонова»[91]. Тож, модель Метцгера та ін. вперше чітко вказала на зв'язок між кілоновими, короткими гамма-спалахами, гравітаційними хвилями та походженням елементів r-процесу (актиноїдів та лантаноїдів).

Різні компоненти викидів із злиття двох нейтронних зір та залежність випромінювань їх кілонових від кута променя зору спостерігача, θ_obs, відносно осі подвійної системи, для двох різних сценаріїв: (ліворуч) швидкого утворення чорної діри та (праворуч) довготривалого магнетарного залишку. Динамічні викиди матерії в орбітальній площині, в обох сценаріях, є багаті нейтронами (з часткою електронів Y_e < 0,1), тому можуть генерувати лантаноїди та пов'язане з ними випромінювання "червоної" кілонови з піком в ближньо-інфрачервоному діапазоні. Тоді як динамічні викиди в полярних напрямках через ударне нагрівання можуть бути досить бідними на вільні нейтрино (Y_e > 0,3), унеможливлюючи синтез лантаноїдів, натомість генеруючи випромінювання "синьої" кілонової в оптичному діапазоні. Найвіддаленіші шари полярних викидів можуть містити вільні нейтрони, розпад яких підтримує ультрафіолетовий транзієнт, що триватиме кілька годин після злиття.

Криві блиску наднових визначаються головним чином радіоактивним розпадом одного типу ядер (ізотопу ⁵⁶Ni), тому їх болометричні криві блиску зазнають експоненціального спаду [30]. В основі ж кілонових є радіоактивний розпад широкого діапазону ядер r-процесу з відмінним періодом напіврозпаду, що обумовлює степеневий спад їх кривих блиску [91]. В перші кілька секунд захоплення нейтронів протягом r-процесу енергія генерується з майже сталою швидкістю, але згодом, коли r-процес "заморожується" і ядра повертаються до стабільного стану, то швидкість генерування енергії наближається до степеневого спаду, ${\displaystyle \propto t^{-\alpha }}$, де α ≈ 1,3 [40].

Різним складовим викинутої в злитті матерії відповідає різна непрозорість, κ. Непрозорість кілоновї визначається головним чином часткою лантаноїдів у викидах матерії: непрозорість викидів з малою часткою лантаноїдів зазвичай становить κ ≈ 0,5 см²/г, тоді як непрозорість багатих лантаноїдами викидів κ ≈ 10 см²/г [97]. Наслідком більшої непрозорості є повільніша часова еволюція кривої блиску та зміщення піку спектрального розподілу енергії до червоних довжин хвиль[96]

За відсутності експериментальних даних про непрозорість в лініях важких елементів r-процесу, в моделі Метцгера та ін. була використана непрозорість багатих на Fe викидів в наднових типу Іа[91]. Згодом Кейсен та ін.[98], Бернес і Кейсен[99] та незалежно Танака і Готокезака [100] вперше визначили світність кілонової, включивши атомні дані про непрозорість ліній для викидів важких елементів r-процесу. Зокрема, в цих обчисленнях було показано, що якщо викиди містять ядра актиноїдів чи лантаноїдів з частково заповненими зовнішніми електронними оболонками f-орбіталі (як відбувається при синтезі ядер r-процесу з масовим числом A≥130), то фотонна непрозорість таких викидів в ультрафіолетовому та оптичному діапазонах має бути в 10-100 разів більшою, аніж для викидів, складених з ядер залізної групи із частково заповненими валентними електронами d-орбіталей[98][100]. Це пов'язано з тим, що частково-заповнена f-орбіталь, з більш щільно розміщеними енергетичними рівнями, уможливлює більше число можливих способів розподілу валентних електронів в цій орбіталі і на порядок більше лінійчатих переходів, аніж в d-орбіталі [98]. Висока непрозорість елементів r-процесу повинна затримувати час еволюції кривої блиску від ~1 дня до ~1 тижня та змістити спектральний пік від видимого (передбаченого моделлю Мецгера та ін.) до ближньо-інфрачервоного діапазону[99][100], призводячи до появи «червоної кілонової» (з непрозорістю κ ≈ 10 см²/г).

Не всі частини викидів зі злиття нейтронних зір обов'язково утворюватимуть важкі ядра r-процесу. Викиди з ${\displaystyle Y_{e}\gtrsim 0.25}$ не матимуть достатньої кількості нейтронів для реакцій захоплення нейтронів, щоб проштовхнути потік нуклонів через другий пік (A≈130) r-процесу [49]. У цьому випадку, як передбачає модель Метцгера та ін., не буде утворення лантаноїдів і викиди із злиття нейтронних зір генеруватимуть яскравіше та синювате оптичне випромінювання, що швидко еволюціонує, т.з. «синю кілонову»(κ ≈ 0,5 см²/г)[91][49]. Синє оптичне випромінювання викидів після злиття є свідченням утворення легших ядер r-процесу, оскільки їх непрозорість лиш трохи вища від непрозорості заліза[97].

Нові теоретичні обчислення та моделювання показали, як включення ефектів переносу нейтрино може призвести до утворення легших ядер r-процесу у викидах після злиття нейтронних зір[73][101][102]. Нагрівання нейтрино відіграє важливу роль в зміні хімічного складу (Y_e) викидів із злиття. Електронні нейтрино та антинейтрино, випромінювані після злиття, характеризуються досить високою світністю (>10⁵³ ерг/cек)[101]. Внаслідок поглинання нейтрино з такою високою світністю викиди стають більш багатшими на протони, оскільки електронні нейтрино перетворюють деяку частку нейтронів в протони через реакції ${\displaystyle n+v_{e}\leftrightarrow p+e^{-}}$[73]. Тому нагріті ударними хвилями полярні динамічні викиди після злиття будуть мати відносно високу частку електронів, ${\displaystyle Y_{e}\gtrsim 0.25}$, і позбавлені лантаноїдів, тим самим роблячи внесок до раннього, синього випромінювання кілонової [101]. Разом з тим, синя кілонова буде видимою лише для променів зору, що не блокуються викидами матерії з високою непрозорістю - «червоною кілоновою»[50].

У порівнянні з колімованим і релятивістсько спрямованим гамма-спалахом, Метцгер та Бергер припустили, що ізотропні випромінювання кілонової роблять їх найбільш перспективним електромагнітним двійником для типового злиття подвійних нейтронних зір на відстані 200 Мпк, що знаходиться в діапазоні модернізованих детекторів LIGO/Virgo [103]. Пізніші обчислення і перші запуски оновленого LIGO дозволили астрономам висунути припущення, що злиття можуть відбуватись значно ближче ніж 200 Мпк, завдяки чому кілонові можна виявити навіть з 1-м телескопами[104].

2013 р. вперше було виявлено інфрачервоне випромінювання після короткого гамма-спалаху GRB 130603B протягом тижня[105][106]. Це відкриття було першим свідченням про прямий зв'язок між злиттям нейтронних зір та короткими гамма-спалахами, а відтак, підтвердженням злиття нейтронних зір як місця утворення важких ядер r-процесу у всесвіті[105].

Існування гравітаційних хвиль вперше було передбачене Альбертом Ейнштейном в 1916 р., через кілька місяців після його публікації рівнянь гравітаційного поля[110]. Відповідно до загальної теорії відносності (ЗВТ), гравітаційні хвилі є збуренням метрики простору-часу і постають як наслідок накладення релятивістської природи на гравітаційні взаємодії. Гравітаційні хвилі генеруються внаслідок асиметричного прискорення системи мас, квадрупольний момент яких змінюється з часом, і поширюються зі швидкістю світла, а їх амплітуда спадає обернено пропорційно відстані від джерела[111]. Гравітаційне випромінювання, генероване орбітальним прискоренням системи мас, призводить до втрати ними енергії, внаслідок чого орбіти зменшуються і врешті решт два масивні тіла зливаються.

Гравітаційні хвилі є хвилями припливної сили. Згідно з принципом еквівалентності, одиничні ізольовані частинки не можуть бути використані для вимірювання гравітаційних хвиль (в силу їх вільного падіння в будь-якому гравітаційному полі і відсутності впливу від проходячої хвилі). Для такого вимірювання потрібні неоднорідності в гравітаційному полі, якими є припливні сили, що переносяться гравітаційними хвилями і які можуть бути виміряні через співставлення положень та взаємодій двох чи більше частинок[110].

Гравітаційне випромінювання в ЗВТ представлене через безслідовий симетричний тензор другого рангу. В загальній системі координат такий тензор має десять незалежних компонент. Однак, подібно до електромагнітного випромінювання, гравітаційне випромінювання в ЗВТ має тільки два незалежні стани поляризації: «+»-поляризацію та «x»-поляризацію (назви пов'язані з формою еквівалентного силового поля, яке вони генерують), що позначаються h₊ і h_x[111]. На відміну від електромагнітних хвиль, кут між двома поляризованими станами складає не ${\displaystyle \pi /2}$, а ${\displaystyle \pi /4}$ (45°). Поляризація гравітаційної хвилі випромінюваного таким джерелом, як подвійна система мас, залежить від орієнтації динаміки в середині цього джерела відносно спостерігача. Відтак, вимірювання поляризації дозволяє визначити орієнтацію подвійної системи[30].

Вплив «+»-поляризованої гравітаційної хвилі на кільце частинок.

Вплив «×»-поляризованої гравітаційної хвилі на кільце частинок.

Гравітаційно-хвильові антени є лінійно-поляризованими квадруполярними детекторами і не чутливі до напрямку хвилі. Тому на основі лише однієї антени не можна визначити напрямок до джерела хвилі. Для цього потрібно одночасне спостереження з використанням трьох і більше детекторів, завдяки чому джерело може бути тріангульоване на небі через вимірювання різниці в часі надходженні сигналу до різних детекторів[111].

Оскільки тканина простору-часу є надто "жорсткою", то амплітуда її деформації дуже мала. Відтак, щоб гравітаційно-хвильовий сигнал міг бути зареєстрований на Землі потрібні зіткнення дуже масивних тіл. Але навіть тоді гравітаційна хвиля, реєстрована земними детекторами, матиме амплітуду лише h~ 10^-21 [30].

На відміну від двох детекторів LIGO, детектор Virgo не виявив GW170817 сигналу, який потрапив у його "сліпу ділянку"[1]. Але цей факт дозволив звузити пошук джерела сигналу на небі до 28 градусів² [112].

Детектори LIGO-Virgo можуть спостерігати гравітаційні хвилі від подвійної нейтронної зорі протягом декількох хвилин. У випадку GW170817, за 100 секунд до зіткнення, нейтронні зорі були на відстані приблизно 400 кілометрів одна від одної й оберталися приблизно 12 разів за секунду. З кожним обертом нейтронні зорі випромінювали гравітаційні хвилі, втрачали енергію й наближалися одна до одної. Зі скороченням орбіти (так званого «падіння по спіралі») збільшувалася швидкість зір, призводячи до збільшення частоти (відомого як чирп, англ. chirp) та амплітуди гравітаційних хвиль. Процес зближення прискорювався, доки дві зорі не злились, утворивши один об'єкт. Випромінена енергія гравітаційної хвилі становила >0,025 M_☉ c²[1].

Гравітаційно-хвильовий сигнал GW170817 є найгучнішим з досі спостережуваних — комбіноване співвідношення сигнал/шум (SNR) становило 32,4 (LIGO-Livingston SNR - 26,4; Ligo Hanford SNR - 18,8; Virgo SNR - 2,0), тоді як SNR для GW150914 лише 24.[2].

Із зареєстрованого GW170817-сигналу, найкраще вимірюваним параметром маси є маса чирпа, ${\displaystyle {\mathcal {M}}}$, (поєднання мас компонент подвійної системи, що обумовлює еволюцію частоти гравітаційного випромінювання і є панівною складовою смуги чутливості детектора), яка становила ${\displaystyle {\mathcal {M}}}$=7000118800000000000♠1.188+0.004
−0.002 M_☉. Тоді як маси зір подвійної системи, внаслідок кореляції їх невизначеностей коливались у широкому діапазоні від 0,86 до 2,26 M_☉.[2].

У цілому, гравітаційно-хвильовий сигнал GW170817, - як і відкриті LIGO перед тим гравітаційні хвилі із зіткнень чорних дір, - цілком узгоджувався із загальною теорією відносності[113][114]:

1. гравітаційні хвилі мають тензорну поляризацію
2. гравітаційні хвилі від зіткнення двох нейтронних зір рухались із тією ж швидкістю, що й світло;
3. гравітаційні хвилі та гамма-промені йшли до Землі 130 млн світлових років і прийшли в межах ~2 секунд;
4. гравітони та фотони із GW170817 події падали до гравітаційного поля Чумацького Шляху в один і той же час, що свідчило про їх падіння з однаковою частотою згідно з принципом еквівалентності.

На основі лише однієї події GW170817 науковцям LIGO-Virgo вдалось встановити частоту злиття подвійних нейтронних зір в рік на середній об'єм простору в кубічний гігапарсек: ${\displaystyle {\mathcal {R}}_{BNS}}$=7003154000000000000♠1540+3200
−1220 Гпк^-3 рік^-1, що відповідає 6 - 120 злиттям подвійних нейтронних зір на рік, коли обсерваторії LIGO-Virgo досягнуть запланованої чутливості детекторів в 2020 р.[1] Така частота злиття узгоджується з частотою, отриманою із попередніх спостережень систем подвійних зір[115][116].

Майже одночасна реєстрація гравітаційно-хвильового сигналу GW170817 та пов'язаного з ним гамма-спалаху GRB 170817A детекторами LIGO та космічними обсерваторіями Fermi/GBM та INTEGRAL. Вгорі: три часові ряди, що відповідають числу фотонів за секунду для двох супутників. Верхні дві панелі - часові ряди Fermi. Ці дві панелі відповідають двом різним діапазонам енергії їх детекторів. Третя панель - часовий ряд INTEGRAL. Зареєстрований Fermi та INTEGRAL гамма-спалах GRB 170817A прийшов через 1,7 сек після гравітаційно-хвильового сигналу GW170817, зареєстрованого LIGO-Livingston (внизу).

Перше повідомлення про GRB170817A було автоматично генероване датчиком гамма-спалахів GBM (англ. Gamma-ray Burst Monitor) на орбітальній обсерваторії Fermi[119] о 12:41:20 UTC, всього через 14 сек після реєстрації ним СГП о 12:41:06 UTC[107]. Згодом, через онлайновий пошук ініційований повідомленнями LIGO-Virgo та Fermi-GBM, GRB170817A був виявлений орбітальною гамма-обсерваторією INTEGRAL використовуючи антиспівпадальний захист (англ. Anti-Coincidence Shield) германієвого гамма-спектрометру (SPI) на його борту [117]. Різниця між GW170817 та GRB170817A становила T0 - t_c = 1,734 ± 0,054 сек[2]

Подальший аналіз даних Fermi-GBM визначив тривалість GRB170817A в T₉₀ = 2,0 ± 0,5 сек, де T₉₀ - інтервал, в межах якого 90% флюенс СГП накопичується в діапазоні енергії 50 - 300 кеВ. Відтак, GRB170817A був класифікований, як короткий спалах гамма-променів із [2]. Піковий потік фотонів GRB170817A, вимірюваний протягом 64 мсек, становив 3,7 ± 0,9 фотонів сек см^-2, а флюенс в межах T₉₀ - (2,8 ± 0,2) × 10^-7 ерг/см^-2 (10 - 1000 кеВ)[107]. GRB170817A є найближчим кСГП з виміряним червоним зміщенням[2].

Затримка електромагнітного сигналу на ~1,7 сек порівняно з гравітаційно-хвильовим може бути пов'язана з ефектом Шапіро (передбачає, що час поширення безмасових часток у викривленому просторі-часі, тобто через гравітаційні поля, дещо збільшується порівняно з гладким простором-часом), хоча не виключені й інші пояснення[1]. Іншими поясненнями затримки гамма-променів може бути[120][121][122][123]:

a) затримка із злиття, допоки не утворилась чорна діра;

б) затримка із злиття, допоки гамма-спалах не був підсилений магнітним полем;

в) час надходження світла від місця випромінювання;

г) уповільнення, зумовлене міжзоряним середовищем

д) екзотична фізика

Така затримка в часі двох сигналів дозволила пов'язати короткі гамма-спалахи зі злиттям нейтронних зір[114]. Крім того, вона накладає обмеження на швидкість гравітаційних хвиль, які за передбаченнями загальної теорії відносності мають ту ж саму швидкість, що й швидкість світла [113]. Тож, відкриття GW170817 та GRB 170817A вперше однозначно показали, що швидкість світла дорівнює швидкості гравітаційних хвиль із точністю до 10^-15[2][1]:

${\displaystyle -3\times 10^{-15}\leq {\frac {v_{gw}-c}{c}}\leq +7\times 10^{-16}}$

GRB170817A характеризувався винятково низьким ізотропним еквівалентом енергії гамма-променів (E_{γ, iso} ≈ 5 × 10⁴⁶ ерг) [114] та піковою фотонною енергією (E_p ~ 40 - 185 кеВ), що разом з пізнім надходженням рентгенівського (через 9 днів)[6] та радіо (через 17 днів)[124][125] випромінювання, є переконливим свідченням випромінювання із відносно вузького релятивістського джету, спостережуваного під великим кутом (θ_obs ≈ 37◦ −42◦) до осі його початкової апертури (θ_obs > θ₀)[126].[127].

Аналіз максимальної правдоподібності ключових параметрів післясвітіння на основі даних радіо- та рентгенівського випромінювання з GRB170817A визначив діапазон енергії для істинної енергії джету кСГП в ⟨E⟩ ≈ 10⁴⁸ - 10⁴⁹[126]. Це, в цілому, узгоджується з енергіями для джетів кСГП, що спрямовані вздовж променя зору спостерігача і які характерні для злиття нейтронних зір[18][95].

Більш детальний аналіз даних Fermi-GRB для GRB170817A виявив два окремі компоненти спалаху: 1) первинний сильний пік, що тривав <0,5 сек, нетепловий комптонівський спектр якого (E_p ~ 185 ± 62 кеВ) був в цілому схожий до регулярних СГПів, ; потім слідував 2) слабший компонент випромінювання, з тривалістю ~2 сек і мав тепловий спектр чорного тіла з kT = 10,3 ± 1,5 кеВ[107]. Перший компонент міг бути позаосьовим проявом більш потужного джету короткого СГП, світність якого була ослаблена релятивістським випромінюванням[128], або ж є свідченням відхилення кута нахилу осі джету СГП від променя зору[129]. В свою чергу джерелом теплової компоненти GRB170817A міг бути гарячий кокон[130][131] чи вихід ударної хвилі внаслідок прориву ультрарелятивістського СГП крізь хмару полярних викидів[132][133].

Оскільки випромінювання червоної кілонови в GW170817 припускають, що після злиття нейтронних зір утворився масивний акреційний диск і формування чорної діри було відносно швидким[120], то система "чорна діра - тор" забезпечує природний механізм для генерування і живлення релятивістського джету СГП[77][134]. Ця система також пояснює часову еволюцію структури джету. Акреційний диск, утворений після злиття, еволюціонує протягом в'язкого часу кілька секунд. За цей час диск втрачає значну частку своєї маси через акрецію та витікання [45]. Якщо колімація джету СГП обумовлена середовищем, сформованим вітрами диску та динамічними викидами, тоді (з послабленням густини довколишньої хмари викидів та сили джету) кут розкриття джету може також збільшитись впродовж кількох секунд, подібно до спосетержуваної затримки GRB170817A в ~1,7 сек[120].

Вірогідність того, що GRB170817A спостерігався під кутом до осі ядра джету кСГП узгоджується з відносно великим кутом нахилу подвійної системи відносно променя зору (θ_obs): θ_obs ≈ 0,2 - 0,6[135]. Іншим свідченням існування більш потужного позаосьового джету в GW170817 було відкриття нетеплових радіо та рентгенівських випромінювань, що слідували після злиття із затримкою в кілька тижнів[124][125] [136][6]. Таке випромінювання характерне для позаосьового післясвітіння від «сирітського» (без реєстрації самого гамма-випромінювання) гамма-спалаху[137].

Тільки через 15 днів після виявлення гравітаційно-хвильового сигналу космічний телескоп «Чандра» зареєстрував рентгеінвське джерело, а опісля був виявлений і радіо сигнал[124][6]. Моделювання показали, що ці рентгенівські та радіо-випромінювання утворилися внаслідок релятивістських джетів з енергією, близькою до космологічних коротких гамма-спалахів, але спрямованих під великим кутом до променя зору, що мало пояснити слабкий короткий гамма-спалах [124].

За отриманими даними, гамма-спалах GRB170817A характеризується наступними рисами [1][114]:

1. розпочався через ~2 сек після злиття подвійної нейтронної зорі;
2. тривав ~2 сек, що близько до тривалості більшості коротких гамма-спалахів;
3. загальна гамма-ізотропна енергія (E_{γ, iso}) ~5 × 10⁴⁶ ерг;
4. пік (E_peak) ~185 кеВ

Тож, GRB170817A набагато слабший, ніж найслабші із зареєстрованих гамма-спалахів, попри те, що джерело гамма-спалаху відносно набагато ближче до Землі і він мав бути набагато більш яскравим [114].

Низька гамма-променева світність GRB170817A не узгоджується із жодним з досі виявлених випромінювань кСГПів [112][138][114][135]. Існує кілька можливих сценаріїв, що пояснюють таку низьку світність кСГП, серед яких найбільш вірогідними є: випромінювання із структурованого джету з ширококутним розподілом [128][139] ; проривне випромінювання із помірно релятивістського кокона [140][141][135][142]; осьове випромінювання множини кСГП з низькою світністю[143][144].

Одним з найбільш обґрунтованих пояснень низької світності GRB170817a є те, що він постав із нагрітого кокона і вісь джету була спрямована під великим кутом до променя зору спостерігача[140][141][142][112]. Злиття двох нейтронних зір супроводжується викидом матерії в довколишнє середовище. З одного боку ці викиди, через лобовий тиск, сприяють колімінації (утворенню паралельних) релятивістських джетів[133]. З іншого боку, коли релятивістський джет проривається через повільні викиди, то значна частка викинутої матерії нагрівається та прискорюється, утворюючи кокон[140]. Розповсюджуючись в поперечному напрямку, кокон розширюється майже ізотропно (до кута ~50°). Стаючи прозорим, такий кокон може підтримувати транзієнт коротких гамма-променів[138][135].

Запропонована Лаццаті та ін. модель передбачає, що поширення джету СГП через баріонні уламки, викинуті довкола місця злиття, призведе до утворення нагрітого ізотропного кокона, через який ранні (протягом перших хвилин після злиття) гамма-промені і будуть спостерігатись[140]. Згідно з цією моделлю, джерелом миттєвої світності GRB 170817A міг бути кСГП, спостережуваний під кутом ~10-20° до осі джету. Однак в оновленій моделі Лаццаті та ін., в основі якої був вже анізотропний кокон, рання світність кСГП припускала спрямованість променя зору спостерігача до осі джету під кутом 40-50°[130]. Однак ця ж сама модель передбачає для реєстрованої Fermi-GBM пікової енергії фотонів GRB170817a (124±52.6 кеВ[107]) кут нахилу 10°. Тож дана модель не повністю пояснює спостережуване випромінювання GRB170817a.

Іншим передбаченням моделей нагрітого кокона є пізня (від кількох тижнів до місяців після злиття) поява рентгенівського випромінювання у випадку позаосьового випромінювання кСГП[130][141]. Причиною такої затримки є те, що по мірі уповільнення трансрелятивістського кокона (з Лоренц-фактором Г~2-3) зовнішнім середовищем та розширення в поперечному напрямку, він починає випромінюватись як регулярне післясвітіння осьового СГП, освітлюючи все більшу частину неба і врешті досягаючи променя зору спостерігача [141]. Передбачається, що випромінювання такого позаосьового "сирітського" (без гамма-променів) післясвітіння повинно бути значно тьмянішим, аніж раннє післясвітіння[145]. Дані спостережень транзієнта GW170817a космічною рентгенівською обсерваторією "Чандра" в цілому узгоджуються з цими моделями нагрітого кокона. На 9 день після злиття нейтронних зір, телескопом Чандра в GW170817 було виявлено лише тьмяне джерело рентгенівського випромінювання [146]. Тоді як на 15 день рентгенівська світність становила L_X ≈ 9 × 10³⁸ ерг/сек [147][6].

Моделі класичного, позаосьового кСГП, - із променем зору спостерігача під кутом до осі сильного, ультрарелятивістського джету, - не в змозі пояснити спостережувані властивості гамма-променів GRB170817a. Тим паче, що у випадку такого сценарію транзієнт GW170817 повинен проявляти яскраве післясвітіння вже приблизно через день після злиття нейтронних зір, що не узгоджується з довгою (від кількох днів до тижнів) затримкою радіо- та рентгенівського випромінювань. Тому для пояснення слабких гамма-променів в події GW170817, Каслівал та ін. запропонували модель ширококутного, помірно-релятивістського (Г≈ 2-3) кокона, що огортає джет кСГП[112]. В цій моделі релятивістський джет запускається після короткої затримки, обумовленої колапсом гіпермасивної нейтронної зорі в чорну діру. По мірі проходження джету через викиди матерії в злитті нейтронних зір, викинута матерія роздувається та утворює герметичний кокон, що розширюється назовні з помірно-релятивістською швидкістю. Згідно з цією моделлю кокона, існують два можливі сценарії для джету: а) у випадку ширококутного джету (кут ≈30°) він виявиться закупореним і не зможе прорватись через викиди матерії; б) у випадку тривалого та вузькокутного джету (≈10°), то він зможе прорватись через викиди та постане для спостерігача як стандартний, осьовий кСГП. Модель Каслівал та ін. передбачає, що всі властивості спостережуваних гамма-променів GRB170817a є наслідком прориву джету через помірно-релятивістський кокон, за умови що Лоренц-фактор джета Г≈2–3, а радіус прориву - ~3 × 10¹¹ см[112].

Виконане Каслівал та ін. релятивістське гідродинамічне моделювання джету всередині розширюючихся викидів матерії нейтронних зір показало, що навіть якщо незначна кількість викинутої матерії (з масою ≈ 3 × 10⁻⁹ M_☉) рухається зі швидкістю 0,8c, то радіус прориву та швидкість будуть достатніми для генерування спостережуваних гамма-променів GRB170817a[112]. В цілому, модель Каслівал та ін. досить добре пояснює багатохвильові властивості транзієнта GW170817, від гамма- до радіо випромінювань.

Фрагмент бесіди в Slack двох астрономів з команди "1M2H", яка першою виявила оптичне джерело після події GW170817 — SSS17a (Swope Supernova Survey 2017a, через 11 годин після відкриття LIGO-Virgo). На другому зображенні помітно спалах у галактиці NGC 4993, де утворився гравітаційно-хвильовий сигнал GW170817

Одразу після реєстрації GW170817 LIGO-Virgo, шість колаборацій астрономів незалежно одна від одної протягом 42 хвилини виявили і почали спостерігати оптичне джерело гравітаційно-хвильового сигналу, згодом назване AT 2017gfo. Ці колаборації включали: 1) колаборацію "Swope Supernova Survey (SSS)" та "One-Meter Two-Hemisphere " (1M2H)[148]; 2) "Distance Less Than 40 Mpc " (DLT40)[149]; 3) групу Dark Energy Survey (DES)[150]; 4) команду відстеження післясвітіння гравітаційних хвиль в "Las Cumbres Observatory" [4]; 5) колаборацію "Mobile Astronomical System of Telescope-Robots" (MASTER)[151] та 6) колаборацію "VIsta Near-infraRed Observations Unveiling Gravitational wave Events " (VINROUGE)[152].

Команда астрономів, використовуючи відносно маленький 1-метровий (в діаметрі) телескоп Swope обсерваторії Лас-Кампанас у Чилі першою виявила яскраве оптичне джерело (SSS17a) в галактиці NGC 4993 (сузір'я Гідри), яка належить до лінзоподібних галактик [153][154], що згодом було підтверджено спостереженнями з багатьох інших телескопів[155][156][148].

Спочатку яскрава світність та синій, невиразний оптичний спектр джерела GW170817 узгоджувався з вибухом молодої наднової. Тому спершу оптичний двійник було зареєстровано як наднову SSS17a (англ. Swope Supernova Survey 2017a). Однак наступної ночі спостереження виявили, що джерело GW170817 істотно потьмяніло в оптичному спектрі, але стало яскравішим в інфрачервоному. Пізніше Центральне бюро астрономічних телеграм Міжнародного астрономічного союзу) перейменувало подію в «AT 2017gfo».

Оптичний транзієнт AT 2017gfo з'явився через ~1 день після злиття нейтронних зір і швидко згасав, з частотою ~2 mag на день в g-смузі, ~1 mag/день в r-смузі і ~0,8 mag/день в i-смузі [4]. Отримані із SOAR) і Магеланових телескопів оптичні спектри та із телескопу Габбла УФ-спектри джерела GW170817 тривали 1,5 - 9,5 днів[157].

Оптичні спектри GW170817 були не схожими на жодні з відомих для наднових[4] і дуже швидко еволюціювали від синього (~6400 K) до червоного (~3500 K) протягом перших трьох днів після злиття[158]. Так, отриманий за допомогою SOAR спектографа англ. Goodman High Throughput Spectrograph найраніший спектр AT 2017gfo, на 1,5 день після злиття, характеризувався панівним синім компонентом (пік ~5000 Å; світність λL_λ≈2×10⁴¹ ерг сек^-1), що вже на 2,5 день змінився червонуватим компонентом (з піком 7000 Å), і цілком змістився з оптичного діапазону на 7,5 день[157].

Злиття двох нейтронних зір, GW170817, супроводжувалось відповідним коротким гамма-спалахом, GRB 170817A. Близько половини із спостережуваних коротких гамма-спалахів мають відповідні післясвітіння[95]. Тому не виключена можливість, що в ранньому, синьому оптичному спектрі злиття може бути присутній домішок такого післясвітіння у випромінюванні із кілонової [158]. У перші дні, AT 2017gfo мав кольори, схожі до раніше спостережуваних післясвітінь коротких гамма-спалахів [157]. Однак, до прикладу, післясвітіння гамма-спалаху GRB 130603B, було на ~4 mag яскравішим, аніж AT 2017 gfo[106][4]. Тому, якщо було післясвітіння пов'язане GW170817/GRB 170817 A, воно було значно слабшим, аніж післясвітіння від GRB 130603B[158].

Оптичні криві блиску кілонови в NGC 4993 Графік показує яскравість кілонової, виявленої в галактиці NGC 4993, виміряну через різні світлофільтри. У синьому світлі оптичне джерело швидко загасло. Але в ближньому інфрачервоному діапазоні об'єкт залишався яскравим деякий час і загасав повільніше. Протягом чотирьох тижнів кілонова змінила колір із яскраво синього на яскраво червоний, що пояснюється r-процесом нуклеосинтезу важких елементів.

Разом з тим, властивості оптичного і УФ спектрів GW170817 не можуть бути пояснені винятково післясвітінням від короткого гамма-спалаху, що зазвичай характеризується більш синім оптичним спектром і значно повільнішою зміною кольору його компонент[159]. Це узгоджуються з даними, отриманими в радіо- та рентгенівському спектрі, що також не виявили значного внеску післясвітіння гамма-спалаху в цей період, вказуючи на відхилення кута нахилу осі джету[136][125].

Світність та швидка зміна кольору AT 2017gfo не можуть бути пояснені простою моделлю випромінюючого абсолютного чорного тіла[157][158][4], однак досить добре узгоджуються з моделями кілонової[91][103][160][100] [8] - оптичного транзієнта, спричиненого радіоактивним розпадом матеріалу r-процесу, викиненого злиттям двох нейтронних зір. Виявлений початковий пік в оптичному спектрі та його швидке згасання через день свідчать про "синю" кілонову, яка, згідно з моделюванням, постає лиш тоді, коли частка електронів (${\displaystyle Y_{e}}$; відношення числа електронів до числа нуклеонів) в значній кількості викинутого матеріалу становить ${\displaystyle Y_{e}\gtrsim 0.3}$, стримуючи утворення лантаноїдів[49].

Оптичні спектри вказують на те, що синя кілонова постала як результат полярних викидів матерії, що складалась з легких ядер r-процесу із масовим числом ${\displaystyle A\lesssim 140}$[157]. Тоді як з утворенням лантаноїдів пов'язаний викид матерії з ${\displaystyle Y_{e}\lesssim 0.2}$. Через більшу непрозорість, випромінювання від багатої лантаноїдами матерії досягають піку в інфрачервоному спектрі протягом тижня після злиття нейтронних зір, формуючи "червону кілонову"[37]. Реєстрований в оптичному діапазоні розподіл спектральної енергії на 2,5-3,5 день був не схожий на спектр жодного із коли-небудь спостережуваних оптичних транзієнтів. Тому жодне з існуючих передбачень обчислювальних моделей чи спостережуваних астрономічних транзієнтів не узгоджується з червоним кольором в оптичному спектрі через 2,5 дні після злиття, за винятком "червоної" кілонової[157].

Еволюція розподілу спектральної енергії оптичного двійника GW170817 (AT2017gfo) протягом перших 12 днів. Помітний перехід від синього до червоного кольору.

Для оптичної світності потрібна низька частка лантаноїдів, оскільки r-процес надважких елементів придушить оптичний потік через непрозорість їх спектральних ліній [160]. Оптичні дані припускають, що такі викиди матерії з незначною часткою лантаноїдів можуть бути видимі лише в межах ${\displaystyle \theta _{obs}\lesssim 45^{\circ }}$ орбітальної осі[157].

На основі порівняння оптичних спектрів з модельними передбаченнями[160], маса викидів "синьої" кілонової повинна становити ~0,003 M_☉[157][158][4]. Відсутність роздільних спектральних ліній в оптичних даних найкраще узгоджується з тими моделями, які передбачають швидкість викидів "синьої" кілонової v = ~0,3c[8]. Згідно з теоретичними розрахунками, така велика швидкість є свідченням того, що "синя" кілонова сформувалась в результаті динамічних полярних викидів матеріалу, нагрітого ударною хвилею в точці зіткнення[38][72]. З цього випливає, що обидва компоненти подвійної системи повинні були бути нейтронними зорями, а не нейтронною зорею і чорною дірою, оскільки в цьому випадку (за відсутності точки зіткнення) єдиним джерелом викидів з високою часткою електронів ${\displaystyle Y_{e}}$ буде лише вітер з акреційного диску[157]. Тож, виявлена в оптичному спектрі значна маса швидких викидів, пов'язаних з «синьою» кілоновою, відкидає можливість того, що джерелом GW170817 є злиття чорної діри та нейтронної зорі. Відтак, оптичні дані суттєво доповнюють гравітаційно-хвильовий сигнал, спираючись лише на який не можна розрізнити складові компоненти подвійної системи.

У випадку, якщо «синя кілонова» дійсно пов'язана з динамічними полярними викидами, то загальна маса викинутого матеріалу дозволяє накласти обмеження на радіус нейтронної зорі подвійної системи[157]. Що компактнішою є нейтронна зоря, тим ближче подвійні зорі можуть наблизитись одна до одної і тим більша орбітальна швидкість при злитті, як наслідок - більша ударна хвиля, яка нагріє та викине більше матеріалу [38]. Моделювання показали, що для нейтронної зорі з малим радіусом (<11 км) маса викидів становить ~10^-2 M_☉, і є на порядок меншою для зір з більшим радіусом (>13 км)[72]. Тому у випадку динамічних полярних викидів, їх більша маса припускає невеликий радіус нейтронної зорі, <12 км [157].

Отож, спостережуваний оптичний двійник GW170817, AT 2017gfo, є першим спектральним свідченням кілонової, підтверджуючи теоретичні припущення, що злиття двох нейтронних зір є місцем r-процесу нуклеосинтезу.

AT2017gfo - джерело GW170817-сигналу. Головна панель — інфрачервоне зображення оптичного транзієнта AT2017gfo (виділено рисками) і його положення в галактиці NGC 4993 (Телескоп Габбла, ІЧ-фільтр F110W). Панелі праворуч показують швидку еволюцію кольору оптичного двійника GW170817 від синього до червоного. Вгорі — зображення до злиття, 2014. Посередині — відкриття оптичного транзієнта. Внизу — через 8,5 діб після відкриття транзієнт згасає й стає червоним.

Тоді як в ультрафіолетовому та оптичному діапазонах спалах згас через кілька днів після реєстрації, в червоному та інфрачервоному діапазонах він тривав ще кілька тижнів. Тож, поряд з оптичними та УФ спостереженнями джерела GW170817, незалежні групи астрономів також спостерігали за ним в інфрачервоному діапазоні[161][4][7][162].

Після виявлення оптичного транзієнта GW170817, група астрономів, використовуючи спектрограф та інфрачервону камеру FLAMINGOS-2 на телескопі Джеміні-підвень в Чилійських Андах, отримала послідовність із семи спектрів ближнього інфрачервоного діапазону (БІЧ) (від 700 до 2500 нм) у період від 1,5 до 10,5 діб після реєстрації гравітаційної хвилі LIGO[161].

Первинна фаза БІЧ спектрів, через півтори доби після злиття нейтронних зір, була дуже гладкою й синьою. Однак, вже на 2,5 день після злиття в БІЧ спектрах з'явився широкий пік (біля 1,05 мкм), який протягом наступних ночей став червоним. А на 4,5 день з'явився другий пік — біля 1,55 мкм, вказуючи на зміну панівних джерел непрозорості[161].

Теоретичні обчислення та моделі припускають, що джерелом яскравого інфрачервоного спалаху після злиття подвійних нейтронних зір є "червона" кілонова — радіоактивний розпад лантаноїдів щойно синтезованих у викинутій після злиття речовині[103][99][160] [163]. Припускається, що на відміну від динамічного викиду в полярному напрямку на межі зіткнення двох нейтронних зір[49], речовина, яка викидається в орбітальній площині припливними силами, має меншу частку електронів ${\displaystyle Y_{e}\lesssim 0,2}$, а відтак має бути багатою на лантаноїди й формувати «червону» кілонову, видиму в БІЧ спектрах[164] [161].

Зіставлення спектрів джерела GW170817 з однокомпонентною моделлю «червоної» кілонової (модель, побудована на основі лише одного елемента — неодиму)[8] у ближньому інфрачервоному діапазоні виявило досить добру узгодженість між модельними й спостережуваними БІЧ спектрами, навіть без додаткового коригування параметрів моделі[161]. Така узгодженість дозволила встановити, що для відтворення БІЧ спектрів через 4,5 доби (спектральні піки біля 1,05 мкм і 1,55 мкм) викид кілонової мав мати масу ~0,04 M_☉, швидкість v = 0,1c і досить великий вміст лантаноїдів , ${\displaystyle X_{lan}\approx 10^{-}2}$. Ці дані БІЧ спектрів, що свідчать про «червону» кілонову, також узгоджуються з комбінованими кривими світності в оптичному та БІЧ діапазонах, отриманими іншою групою астрофізиків за допомогою камери темної енергії (англ. Dark Energy Camera, DECam), телескопів Джеміні-південь та Габбла[162].

Відтак, широкі спектральні піки в ближньому інфрачервоному діапазоні через 2,5 доби після злиття нейтронних зір відповідають значній поширеності лантаноїдів, узгоджуються з даними оптичних та ультрафіолетових спектрів, теоретичними обчисленнями й модельні передбаченнями, що злиття подвійних зір є одним з головних місць r-процесу нуклеосинтезу важких хімічних елементів.

Отримана на основі поєднання ультрафіолетового, оптичного та інфрачервоного спектрів Болометрична світність AT 2017gfo змінювалася від 10⁴² ерг/сек через півдоби після спалаху до 3×10⁴⁰ ерг/сек на 10 добу. На основі оцінки еволюції його ефективної температури, було встановлено, що джерело GW170817 швидко охолоджувалось, від ≈11 000 K до ≈5000 K через добу і до ≈1400 K через 10 діб. Швидкість розширення фотосфери джерела була в діапазоні від 0,3c до 0,1с. Крім цього, було виявлено широкі піки в інфрачервоному спектрі, не схожі на жодні інші виявлені астрономічні транзієнти [112]. Тож, Джерело|поєднання високої швидкості розширення оболонки, швидкого спадання оптичної світності та широких інфрачервоних піків у спектрі дозволили астрономам встановити, що AT 2017gfo пов'язаний саме з GW170817[161][162].

У цілому, криві блиску оптичного джерела GW170817 дуже добре узгоджувались з теоретичними передбаченнями та обчислювальними моделями радіоактивного розпаду важких елементів, утворених під час злиття двох нейтронних зір та появи кілонової[160][91][107].

Разом із тим, пошук високо-енергетичних нейтрино (які мають випромінюватись внаслідок злиття) поблизу джерела сигналу GW170817, не дав результатів. Після ідентифікації галактики NGC 4993, де відбулась подія, астрономи спостерігали її два тижні, прагнучи зареєструвати нейтрино, але не виявили помітного випромінювання[2].

Механізм перетворення енергії в гамма-спалах наразі не достатньо зрозумілий. Більш зрозумілою є природа довгохвильових (від рентгенівських до радіо) післясвітінь. Енергія, вивільнена в результаті вибуху і не випромінена гамма-спалахом, має форму речовини або енергії, що рухається назовні майже зі швидкістю світла. Зіткнення цієї матерії з довколишнім міжзоряним газом утворює релятивістську ударну хвилю, що вільно поширюється міжзоряним простором[165]. Друга, зворотна ударна хвиля може поширитись назад до викинутої вибухом матерії (стадія Сєдова). У цьому процесі речовина в ударних хвилях може підсилити локальні магнітні поля, які, у свою чергу, прискорять заряджені частинки, які випромінюватимуть синхротронне випромінювання в більшій частині електромагнітного спектру[166][167].

Галактика NGC 4993 та післясвітіння гамма-спалаху GRB170817A в оптичному діапазоні (вкладка), зареєстровані. космічним телескопом «Габбл» через 6 днів після реєстрації GW170817

Моделі злиття подвійних нейтронних зір передбачають, що після випромінювань на інших електромагнітних хвилях, має з'явитися також радіо-післясвітіння внаслідок припливного викиду 0,01 - 0,05 сонячних мас речовини на субрелятивітських швидкостях[165][163]. Відповідно до цих моделей, синхротронне випромінювання, яке утворюється внаслідок злиття нейтронних зір, має тривати від кількох місяців до років після злиття з максимум у радіодіапазоні[168]. Оскільки час і яскравість радіо джерела є чутливими до маси та швидкості ударної сили виверження та до густини міжзоряного середовища, то радіо сигнал дозволяє встановити як енергію вибуху внаслідок злиття, так і характеристики навколишнього середовища.

Моделювання також передбачає, що властивості гамма-спалахів та їх післясвітінь зумовлені релятивітським струмінем (джетом), взаємодія якого з середовищем генеруватиме радіо-випромінювання[18]. Однак, у цьому випадку крива блиску радіо сигналу істотно залежатиме від кута між променем зору та напрямком струменя[137].

Лише 2 вересня 2017 р., через 16 днів після реєстрації LIGO-Virgo події GW170817, в радіоастрономічній обсерваторії Дуже Великий Масив (ДВМ) було зареєстровано радіохвилі[124][6]. 5 вересня "Австралійський Компактний Масив Телескопів" (АТCА) також виявив післясвітіння GW170817 у радіодіапазоні. Ці радіоспостереження вказують на два можливі механізми походження радіо хвиль: а) релятивістський викид вибуху спрямований під великим кутом до променя зору; б) після вибуху речовина могла утворити «кокон», що поглинув частину джету й призвів до більшого радіо-випромінювання [124]. Однак, перший механізм не узгоджується з реєстрованим радіопіслясвітінням, яке в такому випадку мало б бути набагато яскравішим, ніж спостерігалося. Слабкість радіопіслясвітіння вказує на користь моделі «кокона» й може свідчити про те, що джерелом радіосигналу міг бути не один гамма-спалах, а спалах гамма-променів[124].

Отриманий різними колабораціями астрономів тепловий спектр оптичного транзієнта GW170817[169][157] [161] є першим прямим підтвердженням моделі кілонови [91]. Форма болометричної кривої (загальне випромінювання у всьому діапазоні електромагнітних хвиль) в цілому узгоджується з передбаченою швидкістю радіоактивного нагрівання (∝ t^-1,3) щойно синтезованих важких ядер r-процесу[91].

Всі астрономічні спостереження транзієнта GW170817 показують, що протягом декількох перших днів електромагнітний двійник характеризувався відмітним синім кольором, що швидко еволюціонував і з спектральним піком в оптичному діапазоні[7][158][6][157][150][170][5] [171][138][162]. Пізніше транзієнт став значно червоним і еволюціонував більш повільно від кількох днів до тижня, з максимальною спектральною інтенсивністю на довжині хвилі ~1,5 мкм [5][112][152][161]. Однак транзієнт GW170817 не мав чітко-визначних спектральних рис, що припускає злиття спектральних ліній внаслідок швидкого (до кількох десятих часток швидкості світла) розширення фотосфери [157]. Разом з тим, в спектрах ближнього інфрачервоного випромінювання спостерігались широкі нерівності[6][161], передбачувані для поглинання лантаноїдів[98].

В цілому, спостереження в оптичних та ближньо-інфрачервоних спектральних діапазонах транзієнта GW170817 узгоджувались з двокомпонентною (синьою та червоною) кілоновою [8][112][4].

Викиди матерії, що породили широкий діапазон спостережуваних електромагнітних хвиль події GW170817, могли постати як із динамічних викидів, так із вітру акреційного диску. Виходячи з підгонки спостережуваних кривих блиску до моделей кілонови та спектрів до обчисленого перенесення випромінювання, "сині" викиди (без лантаноїдів) мали масу ≈ 1 - 2 × 10^-2M_☉, а їх середня швидкість - v ≈ 0,2 c [162][157]. Чисельні моделювання вказують на те, що джерелом такої високої швидкості повинні були бути нагріті ударною хвилею динамічні викиди, а не вітри з акреційного диску[120]. В цьому випадку велика кількість динамічних викидів припускає відносно малий радіус нейтронних зір, що зливались - <11 км[157]. Підтвердження цього результату подальшими чисельними моделюваннями та теоретичними обчисленнями матиме ключове значення для рівняння стану нейтронної зорі[23].

Загальна маса "червоних" викидів (багатих лантаноїдами) становила ≈ 4 × 10^-2 M_☉, з меншою швидкістю розширення, ніж в "синіх" викидах - v ≈ 0,1 c[162][161]. Така велика кількість викидів при такій відносно малій швидкості найкраще узгоджується з вітрами акреційного диску як їх джерела, із великою масою тору > 1,0 M_☉ [120]. Такі вітри із акреційного диску, що мав відносно симетричну геометрію, узгоджуються з результатами спостережень транзієнта GW170817, що не виявили лінійної оптичної поляризації від пізнього випромінювання червоної кілонови[170].

Масивний акреційний тор, що був джерелом викидів червоної кілонової в GW170817, можливий лише у випадку, якщо після злиття нейтронних зір утворилась тимчасово стабільна гіпермасивна нейтронна зоря, а не чорна діра [172][120].

Для спостережуваного із GW170817 випромінювання червоної кілонови викинена дисковими вітрами матерія мала мати частку електронів Y_e <0,25, що згідно з числовим моделюванням еволюції такого диску, передбачає відносно короткий час існування гіпермасивної нейтронної зорі - <100 мсек.[173][174]. Існування такої короткотривалої гіпермасивної нейтронної зорі після злиття узгоджується з передбачуваною помірною кінетичною енергією кілонови ≈ 10⁵¹ ерг[175][176].

Після публікації спостережень електромагітного двійника GW170817 всіма залученими колабораціями астрономів, стало очевидним, що двокомпонента модель не повністю узгоджується з цими даними. Об'єднання та гомогенізація (перехресне калібрування) всього набору ультрафіолетових, оптичних та ближньо-інфрачервоних даних (16 наборів даних, 714 одиничних вимірювань із 38 різних інструментів від 0,45 до 29,4 дня після злиття) змусила до перегляду первинних моделей кілонової в GW170817[121].

Виявилось, що двокомпонента модель кілонової[40][162] не в змозі повністю пояснити всю складність гомогенізованого набору даних. Відтак, для пояснення цих даних двокомпонентна модель була розширена до сферично-симетричної трикомпонентної моделі, до якої був доданий третій, "фіолетовий" компонент (з непрозорістю κ ≈ 3 см²/г)[121]. Цьому компоненту відповідають викиди матерії з часткою лантаноїдів меншою аніж для червоного компоненту, але значно більшою від синього компоненту, тому він з'являється пізніше від синьої кілонової.

При порівнянні двокомпонентної та трикомпонентної моделей з повним набором даних, тільки трикомпонента модель найкраще узгоджувалась з гомогенізованими даними. Згідно зі сферично-симетричною трикомпонентною моделлю, маса і швидкість викинутої радіоактивної матерії r-процесу для кожної з трьох компонент становили[121]:

— Синя кілонова: M = 0,020 M_☉; v = 0,266c

— Фіолетова кілонова: M = 0,047 M_☉; v = 0,152c

— Червона кілонова: M = 0,011 M_☉; v = 0,137c

Отримана загальна маса викидів матерії із злиття нейтронних зір становила ≈ 0,078 M_☉, що є дещо більшою (на ≈ 0,02 - 0,06 M_☉), аніж значення, отримані кількома іншими групами астрономів [112][154]. Параметри в трикомпонентній моделі, в цілому, схожі до двокомпонентної моделі за масами та швидкостями викидів синіх і червоних компонент. Однак в трикомпонентній моделі викиди в червоній компоненті розподілені поміж червоними та фіолетовими компонентами[121].

Співставлення трикомпонентної моделі з повним набором даних виявило, що синя компонента була панівною в усіх діапазонах електромагнітного спектру протягом перших 2-3 днів, тоді як фіолетова компонента була панівною у пізніші часи. Червона компонента, маючи меншу масу викидів, залишалась субдомінуючою протягом всього часу, однак в пізніші часи вносила потрібне випромінювання до червоного діапазону [121]. Початок панування фіолетової компоненти в трикомпонентній моделі на <2-3 день після злиття узгоджується з іншими моделями оптичних та ближньо-інфрачервоних спектрів транзієнта GW170817, які вказують на те, що раннє синє випромінювання найкраще описується викинутою речовиною з градієнтом лантаноїдної частки, яка зростає з часом[161][157].

Згідно з трикомпонентною моделлю, найбільш природним поясненням високої швидкості синіх викидів є відносно багаті протонами (висока частка електронів, Y_e) полярні динамічні викиди, спричинені ударними хвилями із зіткнення двох нейтронних зір[41][72]. В цьому випадку, отримана висока маса викидів ( ≈ 0,02 M_☉) свідчить про малий радіус нейтронної зорі - <12 км[121]. Тоді як найбільш природним джерелом червоних та фіолетових компонент викидів матерії є затримані відтоки (з швидкістю ≈ 0,03 - 0,1c) із акреційного диску, утвореного в злитті [45][46][164]. Відносно висока частка вільних нейтронів в цій матерії (з часткою електронів Y_e < 0,25 - 0,3, необхідною для синтезу ядер лантаноїдів) не узгоджується з довготривалим (>100 мсек) залишком гіпермасивної чи надмасивної нейтронної зорі[50][173], однак досить добре узгоджується з помірною кількістю нейтринного опромінення відтоків із акреційного диску чорної діри[46]. Тож, властивості червоних та фіолетових викидів свідчать про відносно швидке формування чорної діри, як залишку від злиття нейтронних зір в GW170817[121].

Отож, наразі сферично-симетрична, трикомпонентна модель, - із загальною масою викидів ≈ 0,078 M_☉, панівною легкою матерією r-процесу, A<140, та помірними швидкостями ≈ 0,15c, - найкраще узгоджується з кривими блиску всіх діапазонів кілонової в GW170817[121].

Відкриття GW170817 є надзвичайно важливим для космології, оскільки дозволяє безпосередньо виміряти сталу Габбла, H₀ (яка визначає швидкість розширення Всесвіту).

Найбільш поширений спосіб вимірювання H₀ полягає у використанні шкали космічних відстаней (т.зв. «космічної драбини»): поєднання даних про відстані сусідніх зір для визначення відстаней до віддаленіших і застосування отриманих значень для оцінки швидкості розбігання галактик. Однак цей метод не є точним, тому навіть найкращі з наявних сьогодні оцінок сталої Габбла різняться між собою. Значення сталої Габбла, отримане зі спостережень цефеїд і наднових типу Іа[178] є на ~8% більшим, аніж значення, отримане зі спостережень реліктового випромінювання[179]

Виявлення LIGO гравітаційних хвиль відкриває можливість прямо вимірювати сталу Габбла шляхом застосування методу запропонованого ще 30 років тому[180] і пізніше вдосконаленого[181][182]. Цей метод оцінює відстань до галактики, застосовуючи гравітаційно-хвильове спостереження т.зв. «стандартної сирени», яка є гравітаційним аналогом астрономічної стандартної свічки (напр., наднової), знаючи світність якої, можна обчислити відстань.

Стандартні сирени є подвійними системами компактних об'єктів (нейтронних зір або чорних дір), які в міру наближення один до одного (по спіралі) й подальшого злиття, випромінюють частотно модульований гравітаційний сигнал (чирп), що несе інформацію про масу компактних об'єктів. Швидкість, із якою змінюється частота подвійної системи, прямо пов'язана з потужністю генерованих нею гравітаційних хвиль, тобто, наскільки «голосним» є гравітаційно-хвильовий сигнал (звідси й назва — «сирена»). Відстань до джерела визначається шляхом вимірювання амплітуди сигналу. Подібно до того, як спостережувана яскравість зорі залежить від її абсолютної світності та відстані, амплітуда зареєстрованих гравітаційних хвиль залежить як від «гучності» джерела, так від відстані до нього. Аналізуючи частоту такого гравітаційно-хвильового сигналу можна визначити «гучність» (амплітуду) випромінених хвиль. Порівнюючи її з реєстрованою гучністю (амплітудою) можна безпосередньо визначити відстань до джерела.

Реєстрація GW170817 дозволила вперше застосувати метод стандартних сирен для вимірювання сталої Габбла[183]. Через аналіз зареєстрованої амплітуди GW170817 сигналу й моделювання його амплітуди в джерелі вдалось оцінити наскільки вона зменшилась, а отже — визначити відстань до джерела. Поєднавши інформацію про відстань джерела GW170817 з червоним зміщенням його галактики, астрофізики змогли безпосередньо (без космічної шкали відстаней й попередніх вимірювань) виміряти відстань до галактики. Застосувавши баєсовий аналіз, що бере до уваги невизначеність у швидкості та вимірюванні відстані, дослідники отримали нове значення сталої Габбла: ${\displaystyle {\mathcal {H}}_{0}}$=7001700000000000000♠70.0+12.0
−8.0 км сек / Мпк[183]. Попри таку велику невизначеність, отримане лише з однієї події (GW170817) значення H₀ узгоджується з попередніми спостереженнями в електромагнітному спектрі (67 і 72 км сек/Мпк), і цілком незалежне від них. Наступні спостереження злиття нейтронних зір мають накласти додаткові обмеження на це значення й усунути розбіжності між різними експериментами.

Відкриття GW170817, GRB 170817 та кілонови виключило деякі альтернативні моделі гравітації, що намагались пояснити обертання галактик і прискорення всесвіту без залучення темної енергії і темної матерії. Найвідомішими серед таких моделей гравітації є "модифікована ньютонівська динаміка" (MOND) та "ентропійна гравітація".

Виявлення GW170817 вперше підтвердило, що швидкість світла й швидкість гравітаційної хвилі збігаються з точністю до 10^-15. Тоді як багато з моделей модифікованої гравітації передбачають відмінність у швидкостях гравітаційної хвилі та світла (згідно з MOND, гравітаційні хвилі мали були бути зареєстровані на ~445 днів раніше від гамма-променів[184]) чи взагалі заперечують фундаментальність гравітації.

Тож реєстрація лише однієї GW170817 події виключила ці моделі, одночасно наклавши обмеження на інші модифіковані моделі гравітації, які виключають темну матерію та темну енергію і в яких фотони та гравітони рухаються по різних геодезичних лініях[185][186][184][187]. [188]

• Александров Ю. В. Астрофізика : навчальний посібник для студентів напряму «Фізика» класичних університетів – Харків. : ХНУ імені В. Н. Каразіна, 2014. – 216 с.
• Александров Ю. В., Шевченко В. Г. Астрофізика: підручник / Ю. В. Александров, В. Г. Шевченко. – Х.: ХНУ імені В. Н.Каразіна, 2016. – 252 с.
• Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 2.: Теория поля - Рипол Классик, 1958
• Choquet-Bruhat, Y. Introduction to General Relativity, Black Holes, and Cosmology - OUP Oxford, 2014. ISBN 0191644536 - 320 p.
• Schutz B. A First Course in General Relativity. 2 ed. - Cambridge University Press, 2009.
• Boyd R.N. An Introduction to Nuclear Astrophysics - University of Chicago Press, 2007 - 422 p.
• Poisson E., Clifford M.W. Gravity: Newtonian, Post-Newtonian, Relativistic - Cambridge University Press, 2014 - 792 p.
• Yamamoto S. Introduction to Astrochemistry: Chemical Evolution from Interstellar Clouds to Star and Planet Formation - Springer, 2017 - 286 p.

• Попов C. Суперобъекты: Звезды размером с город - Альпина Паблишер, 2015 - 349 с.
• Neutron Stars and Pulsars. - Becker W. ed. - Springer, 2009 - 697 p.
• Haensel P., Potekhin A.Y., Yakovlev D.G.. Neutron Stars 1: Equation of State and Structure - Springer, 2007 - 620 p.
• A. G. W. Cameron, David M. Kahl. Stellar Evolution, Nuclear Astrophysics, and Nucleogenesis - Courier Corporation, 2013 - 208 p.
• Vedrenne G., Atteia J.-L. Gamma-Ray Bursts: The brightest explosions in the Universe - Springer, 2009 - 580 p.
• Rezzolla L., Zanotti O. Relativistic Hydrodynamics - OUP Oxford, 2013 - 735 p.
• Masaru S. Numerical Relativity - World Scientific, 2015 - 844 p.