Нумеричке симулације судара галаксија и настанка морфолошких подструктура
др Станислав Милошевић — у време одбране асистент на Катедри за астрономију Математичког факултета Универзитета у Београду (према реферату комисије и биографији у раду)
Универзитет у Београду, Математички факултет; 2023 (извештај комисије од 17. јула 2023; у НаРДуС-у рад носи датум 30. септембар 2023); 110 страна, основни текст на странама 1–97, 160 библиографских јединица; рад је писан латиницом: „Numeričke simulacije sudara galaksija i nastanka morfoloških podstruktura“; ментор др Мирослав Мићић; комисија: М. Мићић, Б. Вукотић, Д. Илић
Пун текст рада и извештај комисије у Националном репозиторијуму дисертација у Србији (НаРДуС). Осврт: „Стручни критички осврт на дисертацију „Нумеричке симулације судара галаксија и настанка морфолошких подструктура““, 21. септембар 2026., 18 страна, с пописом литературе.
Дисертација симулацијама N-тела испитује настанак звезданих токова и љуски при судару велике спиралне галаксије и патуљастог сателита, и показује да се посматрана расподела металичности дуж Великог звезданог тока у халоу галаксије M31, као и металичност две звездане љуске, може објаснити почетним, линеарно опадајућим градијентом металичности у патуљку. Осврт налази да је реч о солидном и добро изведеном нумеричком раду, чији су резултати објављени у рецензираним часописима и који одликује ретка методолошка уздржаност. Слабости су тројаке: мера слагања с посматрањима погрешно је описана и није пондерисана грешкама, испитан је само један облик почетне расподеле, а рад не приказује тест стабилности модела ни анализу конвергенције и не даје верзије кодова ни доступност података; тврдња о првенству изнета је шире него што грађа носи. Ниједна не руши носеће тезе.
Основне примедбе
| # | Примедба | Шта рад каже или ради | У чему је проблем | Стране | Мерило |
|---|---|---|---|---|---|
| 1 | Мера слагања погрешно описана и непондерисана | Носећа процена да је најбоље слагање расподеле металичности на 2,9 милијарди година почива на мери коју текст зове збиром квадрата („Najpre je izračunata suma kvadrata razlike ove dve vrednosti“), док једначина даје збир апсолутних одступања; мера није пондерисана грешкама посматрања и нема праг прихватања. | Процена параметра без исказане неодређености остаје без мере поузданости. | 74 | Andrae 2010; Andrae, Schulze-Hartung, Melchior 2010; Hogg, Bovy, Lang 2010 |
| 2 | Неодређеност градијента је распон претраге, не грешка | Почетни градијент дат је као −0,3 ± 0,2, при чему је ±0,2 полураспон унапред испитаних вредности од −0,1 до −0,5, а не изведена грешка процене. | Вредност се лако чита као мерна неодређеност иако то није; треба је изричито назвати распоном допустивих градијената. | 70–74 | Andrae 2010 |
| 3 | Испитан само један облик почетне расподеле | Рад сам наводи да „Pokazano je (posmatrački i teorijski) da ti gradijenti mogu biti i pozitivni i negativni“, а испитује само линеарно опадајући градијент, коме варира магнитуду: „Jedino je variran gradijent inicijalne raspodele metaličnosti u patuljastoj galaksiji.“ | Слагање подупире негативан линеаран градијент као доследан модел, али не искључује други облик исте средње вредности; закључак о прогенитору важи под том претпоставком. | 64, 70 | аргумент из грађе самог рада |
| 4 | Монте Карло узоркује само металичност, не и динамику | Хиљаду понављања варира само доделу металичности честицама при фиксираној динамици: једној реализацији почетних услова, орбите, масе и профила халоа. | Наведене грешке по пољу одражавају расипање доделе металичности, а не највеће изворе неодређености. | 70–74 | JCGM 101:2008 |
| 5 | Тврдња о првенству изнета шире него што грађа носи | Закључак каже: „Prvi put smo dali objašnjenje raspodele metaličnosti u GSS i zvezdanim ljuskama.“ Концепт објашњења преко градијента прогенитора увели су раније радови које дисертација сама цитира. | Оригиналан допринос који стоји јесте прва квантитативна репродукција расподеле дуж тока по деветнаест поља, с два пика приписана двема кинематичким групама; тврдњу у закључку треба тако разграничити. | 96 | аргумент из грађе самог рада (Miki и др. 2016; Kirihara и др. 2017; Fardal и др. 2008 које рад цитира) |
| 6 | Без теста стабилности и анализе конвергенције | Рад не приказује развој модела у изолацији пре судара ни анализу конвергенције по броју честица, ублажавању и временском кораку; за све честице користи јединствену дужину ублажавања од 100 pc, иако се масе честица разликују за ред величине. | Равнотежни модели коришћеног пакета по речима његових аутора јесу „susceptible to non-axisymmetric instabilities“, а неравнотежни почетни услови воде вештачкој промени профила и убрзаном плимском распаду сателита. | 38, поглавље 2 | Widrow, Pym, Dubinski 2008; Kazantzidis, Magorrian, Moore 2004; Power и др. 2003; Athanassoula и др. 2000; Merritt 1996 |
| 7 | Поновљивост: верзије кодова и доступност података | Кодови Gadget2 и GalactICs нису наведени с бројем верзије ни адресом, нема изјаве о доступности података, а Табела 2 за осам сценарија патуљка не даје ни број честица ни масу честице. | Масена резолуција тих сценарија не може се реконструисати из рада, а часопис у коме је аутор објавио изричито тражи изјаву о доступности података. | поглавље 2, Табела 2 | Smith, Katz, Niemeyer 2016; Sandve и др. 2013; Wilkinson и др. 2016; упутство ауторима часописа MNRAS (снимак од 19. маја 2023) |
| 8 | Морфолошки закључци из визуелног очитавања мапа | Број токова и љуски, њихове удаљености и време распада остатка изведени су визуелним очитавањем мапа честица, без бројчане мере; правило „Ukoliko se formiraju dva zvezdana toka, možemo pouzdano pretpostaviti da je roditeljska galaksija patuljasta sa diskom.“ преоштро је, јер је број токова временски променљив. | Закључци почивају на једној маси, једном растојању и једној реализацији по сценарију, па их треба читати као демонстрацију типова структура, не као статистички поткрепљено уопштавање; физички распад остатка није разграничен од могућег нумеричког. | 55, поглавље 3 | van den Bosch, Ogiya 2018 |
| 9 | Формалне омашке у једначинама | Поасонова једначина написана је с погрешним предзнаком, безсударна Болцманова с погрешним предзнаком члана с потенцијалом, запис профила централног овала одступа од каноничног облика пакета, а иста ознака за радијус халоа мења значење у суседним једначинама. | Све су формалне омашке: потенцијал и силу рачунају пакети у исправном облику, па се погрешно сложене једначине не преносе ни на један резултат. | поглавља 1 и 2 | Binney, Tremaine 2008; Widrow, Pym, Dubinski 2008 |
| 10 | Уводни преглед, ознаке и научни апарат | Опсег звезданих маса патуљастих сфероидних галаксија дат је као 10⁷–10¹⁰ сунчевих маса, што прецењује стварне вредности за неколико редова величине; симулација Миленијум II приписана је раду за изворну симулацију; рад Miki и сар. (2016) сврстан је у интервал судара од 2–3 милијарде година иако даје знатно краћу скалу; две јединице у попису односе се на исти рад; јединица Хаблове константе одштампана је непотпуно. | Подаци не улазе у моделе, па је реч о чињеничким и редакцијским исправкама у уводном прегледу и апарату. | 8, 16, 24, 37 | McConnachie 2012; Mateo 1998; SI брошура, 9. изд. 2019; ISO 690:2021; Правопис српскога језика 2010 |
Шта рад доноси
- Квантитативна репродукција посматране расподеле металичности дуж Великог звезданог тока по деветнаест посматрачких поља, с објашњењем два пика двема кинематички различитим групама честица.
- Јединствен сударни сценарио: истим негативним градијентом металичности у патуљку истовремено се репродукују ток и обе звездане љуске, чиме металичност љуски постаје независан тест исте почетне расподеле.
- Контролисан оглед о утицају морфологије сателита, нагиба орбите и смера ротације на настале структуре, при истој маси и истом почетном растојању.
- Теоријски увод повезан с методом: теорема виријала изведена је у облику непосредно употребљивом за проверу равнотеже модела.
- Методолошка уздржаност: све је фиксирано осим онога што се испитује, квалитативна поређења изричито су тако означена, а домет метода признат је на правом месту.
Образложење примедби
- Мера слагања погрешно описана и непондерисанаЗбир апсолутних одступања ваљано рангира тренутке, а минимум на 2,9 милијарди година сагласан је с независним динамичким показатељима за интервал 2,4–2,9.
- Неодређеност градијента је распон претраге, не грешкаРеч је о смисленој мери распона у коме модел описује посматрања.
- Испитан само један облик почетне расподелеИзбор није произвољан: циљ је постављен условно, а облик је посматрачки мотивисан особинама патуљастих сфероидних галаксија и репродукује оба посматрана градијента дуж тока.
- Монте Карло узоркује само металичност, не и динамикуОрбита и маса нису слободни параметри него су ограничени посматраним даљинама и брзинама дуж тока.
- Тврдња о првенству изнета шире него што грађа носиОбјашњење два пика привременим структурама при постојању најмање две кинематичке групе звезда проверљив је и оригиналан налаз, што и извештај комисије тачно издваја.
- Без теста стабилности и анализе конвергенцијеПакет изводи почетне услове из равнотежне дистрибуционе функције, па су модели по конструкцији приближно равнотежни, а дужина ублажавања мала је у односу на скале дужине компоненти.
- Поновљивост: верзије кодова и доступност податакаОба кода наведена су својим дефинишућим рецензираним радовима, а модели за ток потпуно су задати масама и бројем честица у Табелама 1 и 3.
- Морфолошки закључци из визуелног очитавања мапаПоглавље је вредно као контролисан оглед у коме се мења само један чинилац, а резултати су објављени у рецензираном раду.
- Формалне омашке у једначинамаТеорема виријала изведена је у облику непосредно употребљивом за проверу равнотеже, чиме је увод повезан с оним што рад стварно ради.
- Уводни преглед, ознаке и научни апаратУсклађеност позивања и пописа литературе потпуна је у оба смера, посредног цитирања нема, а темељни радови цитирани су из изворника.
Осврт не доводи у питање оригиналност рада и не тврди повреду академске честитости. Реч је о стручној критици тачности, метода и закључивања, с дословним наводима и странама дисертације. Као мерило су коришћени само литература, стандарди и прописи који су постојали пре предаје рукописа комисији (јун 2023).