Як ви вже здогадалися, це питання цікавить команду LIGO. Одночасно з публікацією згаданого вами доповіді про оголошення відкриття, команда LIGO представила ряд супутніх документів з подальшими подробицями про відкриття та прогнози. Один із таких питань стосується вашого питання:
Швидкість злиття двійкових чорних отворів, отриманих на основі розширених спостережень LIGO навколо GW150914
Метод їх оцінки швидкості подій враховує як GW150914, так і ще одну значно слабку (і менш статистично значущу) подію. Вони розглядають ряд моделей того, як може залежати швидкість подій від властивостей системи, і запитують, що спостереження GW150914 та інших кандидатів події означають для загальної швидкості. Результати різняться від моделі до моделі, але вони обрали моделі, які, як вони вважають, можуть приблизно відповідати астрофізично правдоподібній поведінці. Як узагальнено в рефераті:
2 - 53G p c- 3y r- 16 - 400G p c- 3y r- 12 - 400G p c- 3y r- 1
Зауважте, що стаття подається, а не публікується, тобто все ще знаходиться під експертним переглядом. Виступаючи як хтось, хто має досвід таких розрахунків, деякі аспекти методу мені здаються непосидючими, тому я думаю, що варто перевірити цю статтю за кілька тижнів на предмет доопрацювання. Це не вимагає фантазії методології, щоб побачити, що порядок величини тут (від декількох до ~ 100 на кубічний гігапарсек на рік) знаходиться в правильному напрямку. Але в роботі представлена методологія, яка могла б зробити більш детальні та точні оцінки та прогнози по мірі накопичення даних, тому важливо переконатися, що методологія є надійною.