На логарифмическом графике отношения силы тяжести на поверхности к массе планеты, что означает точка пересечения с осью y?

Я играю с данными с сайта exoplanets.org и меня интересует зависимость гравитации поверхности от массы планеты. Я воспроизвел этот график этот сюжетпосле загрузки их данных и выполнил модель нелинейной регрессии, чтобы соответствовать кривой. К сожалению, ковариационная матрица моей подгонки имеет бесконечные значения, поэтому сейчас я пробую линейную подгонку на логарифмическом графике, показанном ниже. Моя подгонка, для у знак равно а Икс + б , является а знак равно 0,9511 , б знак равно 0,8631 .

Я сейчас думаю о том, что я рисую. Подозреваю, что ничего интересного в нем нет. журнал грамм М р 2 против. журнал М , но, тем не менее, я пытаюсь понять, есть ли смысл в y-перехвате.

Каковы возможные объяснения значения б ?

Учитывая значения для вашего графика, я бы подумал, что ошибка с определением наклона и отсутствие точек между 0,5 и 2 массами Земли ответственны за большое отклонение объекта почти без массы, имеющего значительную поверхностную гравитацию. Было бы лучше построить поверхностную гравитацию для объектов массой от 0,5 до 2 масс Земли, чтобы увидеть, есть ли статистическая ошибка.
Кстати, а откуда данные?
Спасибо за ваши комментарии. Данные взяты с сайта exoplanets.org. Теперь, когда я думаю об этом, log(1)=0, так что точка пересечения с осью y представляет поверхностную гравитацию 1 М Дж планета, нет?
Ты прав. Думаю, мне стоило перепроверить свои расчеты. Итак, поверхностная гравитация планеты с л о грамм ( М О ) знак равно 0 является л о грамм ( а ) знак равно 0,8631 , или а знак равно 7.3 с м / с 2 . Вы уверены, что у вас есть правильные единицы измерения? С 9,98 м / с 2 знак равно 998 с м / с 2 и л о грамм ( 998 с м / с 2 ) знак равно 2,999 .
Да, это сбивает с толку... Прямо сейчас я думаю, что, поскольку гравитация на поверхности также зависит от радиуса, эта линия дает информацию о радиусе. Но ничего особо полезного, т.к. журнал грамм М р 2 против журнал М идет как 1 + к журнал р журнал М , я думаю, после использования некоторых правил логарифмирования. Пожалуйста, дайте мне знать, если это не звучит правильно для вас, я очень новичок в этом типе вещей!
Я просмотрел исходные данные и заметил, что вы нанесли на карту только газовых гигантов (таких как Юпитер). Теперь становится более логичным, что газообразная планета с массой Земли (которая не существовала бы из-за скорости газов; газы просто рассеялись бы) имела поверхностную гравитацию примерно в 100 раз меньше, чем на Земле. Кроме того, поверхностная гравитация каменистых планет увеличивается по мере их уменьшения, поэтому значение для Земли (2,998) находится в середине данных.
Спасибо за этот последний комментарий. Я пытался избавиться от нулевых значений и вырезать значения меньше 0,2. М Дж , который, как я теперь понимаю, слишком велик.
b выглядит как хороший пример опасности экстраполяции кривых наилучшего соответствия за пределы диапазона фактических данных. Планета с нулевой массой не имеет гравитации, любой другой ответ — иллюзия.

Ответы (2)

Я думаю, то, что вы установили здесь, это как раз то, что р имеет тенденцию к увеличению с массой. Плотность планет непостоянна.

Позволять р знак равно р 0 ( М / М е а р т час ) α , так что М знак равно ( 4 / 3 ) π р 3 р 0 ( М / М е а р т час ) α

Затем

грамм знак равно грамм М р 2 знак равно 4 π грамм 3 р р

Заменять р с ( 3 М / 4 π р ) 1 / 3 так что

грамм знак равно 4 π грамм 3 ( 3 М 4 π р ) 1 / 3 р
грамм знак равно ( 4 π 3 ) 2 / 3 грамм М 1 / 3 р 0 2 / 3 ( М / М е а р т час ) 2 α / 3
грамм знак равно ( 4 π 3 ) 2 / 3 грамм М е а р т час 1 / 3 р 0 2 / 3 ( М / М е а р т час ) ( 2 α + 1 ) / 3

Таким образом, за исключением (весьма возможной) алгебраической оговорки, если вы рисуете журнал грамм против журнал М , градиент ( 2 α + 1 ) / 3 , который из вашего сюжета дает α 0,92 - т.е. средняя плотность планеты увеличивается почти линейно с массой.

Тогда перехват

б знак равно журнал [ ( 4 π 3 ) 2 / 3 грамм М е а р т час 1 / 3 р 0 2 / 3 ] ,
что дает р 0 3,5 кг/м 3 (NB: я вычла 2 из вашего б сделать это СИ; давая плотность около 814 кг/м 3 при массе Юпитера).

Тот факт, что плотность почти пропорциональна массе, можно найти из того же набора данных. например см. ниже. Ниже 0,1 массы Юпитера связь, по-видимому, нарушается, хотя на самом деле очень немногие из плотностей таких планет измерены очень точно (поскольку для этого требуется радиус от транзита), но это работает достаточно хорошо в диапазоне, который вы начертили. . Физика здесь заключается в том, что газовые гиганты управляются частично (электронно) вырожденным уравнением состояния, в результате чего все они имеют одинаковый радиус от примерно одной десятой массы Юпитера до примерно 50 масс Юпитера (хотя и со значительным и в значительной степени необъяснимым разбросом) . Таким образом, плотность пропорциональна массе. Это соотношение не работает для небольших каменистых планет, у которых радиус уменьшается для меньших масс.

Плотность планеты в зависимости от массы

Поверхностная гравитация против массы

Для планет с постоянной средней плотностью у вас есть:

М знак равно р × 4 π р 3
и поверхностное значение грамм является:
грамм ( р ) знак равно грамм М р 2 знак равно грамм × р × 4 π × р
Таким образом, для тел постоянной плотности сила тяжести на поверхности пропорциональна радиусу, а наклон как р 0 говорит вам плотность. Так, для тел одинаковой плотности журнал ( грамм ( р ) ) в качестве р 0

По массе:

грамм ( М ) знак равно грамм ( 4 π р ) 2 / 3 М 1 / 3
Таким образом М 0 у нас есть грамм ( М ) 0 и снова журнал ( грамм ( М ) ) и перехват журнал журнал сюжет как М стремится к нулю, позволяет рассчитать среднюю плотность.

Разве это не должно быть грамм знак равно ( 2 3 6 π 2 3 ) ( грамм р 2 М 3 ) ? Обратите внимание, что вопрос ОП на самом деле в основном касается журнал М 0 , но в остальном, я согласен, самая простая физическая вещь, которой он соответствует, - это средняя плотность, подходящая к данным.
Разве проблема с этой идеей не в том, что журнал грамм знак равно ( 1 / 3 ) журнал М + с о н с т , что не то, что получает ОП?
@RobJeffries ОП далеко экстраполирует газовых гигантов, так что р также сильно зависит от М -- но это не должно изменить соответствие средней плотности для любого конкретного фиксированного М , что на самом деле было вопросом ОП, как указано. Правда, экстраполяция ОП физически неуместна, а для легких экзопланет ( М 4 в единицах р знак равно М знак равно 1 ), Сигер и др . (2007)
журнал 10 р знак равно к 1 + 1 3 журнал 10 М к 2 М к 3 ,
и так журнал 10 грамм знак равно 1 3 журнал 10 М 2 к 1 + 4 к 2 М к 3 .
Который 4 к 2 должно быть 2 к 2 в моем комментарии (опечатка).