Daisuke Kawata
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Daisuke Kawata has a total of 397 co-authors. Here are the top co-authors by number of shared publications:
- Brad K. Gibson
E.A. Milne Centre for Astrophysics, University of Hull, Hull, HU6 7RX
57 shared publications - Chris B. Brook
Instituto de Astrofisica de Canarias
27 shared publications - Jason A. S. Hunt
University of Surrey
18 shared publications - Robert J. J. Grand
Astrophysics Research Institute, Liverpool John Moores University
17 shared publications - Mark Cropper
Mullard Space Science Laboratory, University College London
14 shared publications Robert J. J. Grand
Heidelberger Institut für Theoretische Studien Schloss-Wolfsbrunnenweg 35 69118 Heidelberg Germany
10 shared publications
Latest Publications
- Open Access
- 24 July 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 551 (1)
Orbital migration and heating history of the Galactic disc: a transition between the bimodal discs
Stellar orbits in the Galactic disc evolve from their birth to the current shape through both radial migration and dynamical heating. The history of their secular evolution is imprinted in the current kinematics and age–metallicity distribution. We construct a chrono-chemo-dynamical model of the disc, incorporating inside–out growth, metallicity evolution, radial migration and heating to fit the observed age–metallicity–kinematics distribution of LAMOST subgiant stars in both the low and high-$\alpha$ disc. By modelling all distribution parameters with spline fitting, we present the first non-parametric stellar migration and heating history of the Galaxy. We determine the heating-to-migration ratio, the ratio of the root-mean-square changes in radial/vertical and azimuthal actions, to be ${\approx} 0.075$ for radial to azimuthal actions and ${\approx} 0.015$ for vertical to azimuthal actions, implying a highly anisotropic diffusion in the action space. Furthermore, we identify a transition in radial migration efficiency coinciding with the transition moment of the bimodal disc, for which the radial migration was more efficient for the high-$\alpha$ disc than for the low-$\alpha$ disc. This transition may be attributed to two correlated scenarios: (1) a bar formation epoch accompanied by violent outward migration, and (2) a drop in the gas mass fraction in the disc when the low-$\alpha$ disc began to form. These findings offer further constraints on the formation mechanisms of bimodal discs, favouring the downsizing scenario. We also briefly discuss the connection between our results and recent high-redshift observations. In addition to the secular evolution history, our model maps the Milky Way interstellar medium metallicity gradient at different lookback times, which we find has only varied a little (in the range of $-0.07~\rm to~-0.10~dex\,kpc^{-1}$) since disc formation.- 24 July 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 551 (1)
Open AccessOrbital migration and heating history of the Galactic disc: a transition between the bimodal discs
Stellar orbits in the Galactic disc evolve from their birth to the current shape through both radial migration and dynamical heating. The history of their secular evolution is imprinted in the current kinematics and age–metallicity distribution. We construct a chrono-chemo-dynamical model of the disc, incorporating inside–out growth, metallicity evolution, radial migration and heating to fit the observed age–metallicity–kinematics distribution of LAMOST subgiant stars in both the low and high-$\alpha$ disc. By modelling all distribution parameters with spline fitting, we present the first non-parametric stellar migration and heating history of the Galaxy. We determine the heating-to-migration ratio, the ratio of the root-mean-square changes in radial/vertical and azimuthal actions, to be ${\approx} 0.075$ for radial to azimuthal actions and ${\approx} 0.015$ for vertical to azimuthal actions, implying a highly anisotropic diffusion in the action space. Furthermore, we identify a transition in radial migration efficiency coinciding with the transition moment of the bimodal disc, for which the radial migration was more efficient for the high-$\alpha$ disc than for the low-$\alpha$ disc. This transition may be attributed to two correlated scenarios: (1) a bar formation epoch accompanied by violent outward migration, and (2) a drop in the gas mass fraction in the disc when the low-$\alpha$ disc began to form. These findings offer further constraints on the formation mechanisms of bimodal discs, favouring the downsizing scenario. We also briefly discuss the connection between our results and recent high-redshift observations. In addition to the secular evolution history, our model maps the Milky Way interstellar medium metallicity gradient at different lookback times, which we find has only varied a little (in the range of $-0.07~\rm to~-0.10~dex\,kpc^{-1}$) since disc formation. - Open Access
- 19 June 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 550 (2)
A simulation-based inference of the Milky Way merger history
Accreted stars in the Milky Way (MW) preserve information about their galaxies of origin in their chemo-dynamical properties. In this study, we use the chemo-dynamical signatures in the merger debris to approximate the posterior distribution of the properties of disrupted satellites at their accretion. Adopting a simulation-based inference framework, we train an ensemble of normalizing flows using samples of merger debris from the Auriga simulations. Applying this methodology to a local sample of accreted stars in the MW, we infer the lookback times, stellar and halo masses, and halo mass merger ratios of several known accretion events: Gaia-Enceladus-Sausage (GES), Helmi streams, Heracles, I’itoi, LMS-1/Wukong, Sagittarius, Sequoia, and Thamnos. Our predictions for accretion time and mass align with previous estimates and follow the expected stellar mass–metallicity relation across redshifts. The total accreted stellar mass from these events is predicted to be $2.2^{+1.1}_{-0.6}\times 10^{9}~\rm {M_{\odot }}$, with the GES and Sagittarius dwarfs being the largest contributors. The total accreted stellar mass from all disrupted progenitors is $1.4^{+1.1}_{-0.4}\times 10^{9}~\rm {M_{\odot }}$, which is consistent with previous mass measurements of this component. We provide a prediction for the evolution of the MW halo mass until the accretion of Sagittarius ($z\approx 1$): We find that the halo mass growth from the time of its first merger ($z\approx 5$) to $z\approx 2$ significantly exceeds the total mass of the known progenitors accreted during that interval, suggesting the presence of unidentified substructures. Our estimate of the Galaxy halo mass after the Sagittarius merger, but prior to the accretion of the Magellanic Clouds, is $5.9^{+1.4}_{-1.1}\times 10^{11}~\rm {M_{\odot }}$.- 19 June 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 550 (2)
Open AccessA simulation-based inference of the Milky Way merger history
Accreted stars in the Milky Way (MW) preserve information about their galaxies of origin in their chemo-dynamical properties. In this study, we use the chemo-dynamical signatures in the merger debris to approximate the posterior distribution of the properties of disrupted satellites at their accretion. Adopting a simulation-based inference framework, we train an ensemble of normalizing flows using samples of merger debris from the Auriga simulations. Applying this methodology to a local sample of accreted stars in the MW, we infer the lookback times, stellar and halo masses, and halo mass merger ratios of several known accretion events: Gaia-Enceladus-Sausage (GES), Helmi streams, Heracles, I’itoi, LMS-1/Wukong, Sagittarius, Sequoia, and Thamnos. Our predictions for accretion time and mass align with previous estimates and follow the expected stellar mass–metallicity relation across redshifts. The total accreted stellar mass from these events is predicted to be $2.2^{+1.1}_{-0.6}\times 10^{9}~\rm {M_{\odot }}$, with the GES and Sagittarius dwarfs being the largest contributors. The total accreted stellar mass from all disrupted progenitors is $1.4^{+1.1}_{-0.4}\times 10^{9}~\rm {M_{\odot }}$, which is consistent with previous mass measurements of this component. We provide a prediction for the evolution of the MW halo mass until the accretion of Sagittarius ($z\approx 1$): We find that the halo mass growth from the time of its first merger ($z\approx 5$) to $z\approx 2$ significantly exceeds the total mass of the known progenitors accreted during that interval, suggesting the presence of unidentified substructures. Our estimate of the Galaxy halo mass after the Sagittarius merger, but prior to the accretion of the Magellanic Clouds, is $5.9^{+1.4}_{-1.1}\times 10^{11}~\rm {M_{\odot }}$. - Open Access
- 16 June 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 550 (2)
Milky Way’s warped disc traced by AGB stars
While the presence of the Galactic warp has long been established from observations of H I gas, the Gaia measurements of over 1 billion stars with parallaxes have enabled much more detailed studies using stellar populations. Here, we demonstrate that asymptotic giant branch (AGB) stars, an evolved phase of low- and intermediate-mass stars, can serve as an effective tracer of the Galactic warp. We use two distinct AGB populations: C-rich AGB stars, representing stars of about 1 Gyr in age with main-sequence masses of 2–2.5 M⊙, and intermediate-mass (3–5 M⊙) O-rich AGB stars, corresponding to ages of 100–300 Myr. The downward warp traced by O-rich AGB stars is consistent with that found from Cepheids, which is expected given their similar ages. The more numerous C-rich AGB stars clearly reveal the Galactic warp over a wide range of azimuthal angles. Their warp appears to reach larger amplitudes than that of Cepheids across azimuthal angles. Our results show that C-rich AGB stars, together with intermediate-mass O-rich AGB stars, provide new constraints on the Galactic warp at intermediate stellar ages, offering a new insight into the stellar age and warp amplitude relation.
- 16 June 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 550 (2)
Open AccessMilky Way’s warped disc traced by AGB stars
While the presence of the Galactic warp has long been established from observations of H I gas, the Gaia measurements of over 1 billion stars with parallaxes have enabled much more detailed studies using stellar populations. Here, we demonstrate that asymptotic giant branch (AGB) stars, an evolved phase of low- and intermediate-mass stars, can serve as an effective tracer of the Galactic warp. We use two distinct AGB populations: C-rich AGB stars, representing stars of about 1 Gyr in age with main-sequence masses of 2–2.5 M⊙, and intermediate-mass (3–5 M⊙) O-rich AGB stars, corresponding to ages of 100–300 Myr. The downward warp traced by O-rich AGB stars is consistent with that found from Cepheids, which is expected given their similar ages. The more numerous C-rich AGB stars clearly reveal the Galactic warp over a wide range of azimuthal angles. Their warp appears to reach larger amplitudes than that of Cepheids across azimuthal angles. Our results show that C-rich AGB stars, together with intermediate-mass O-rich AGB stars, provide new constraints on the Galactic warp at intermediate stellar ages, offering a new insight into the stellar age and warp amplitude relation.
- Open Access
- 13 March 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 548 (1)
Correction to: Towards Galactic archaeology with inferred ages of giant stars from Gaia spectra
- 13 March 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 548 (1)
Open AccessCorrection to: Towards Galactic archaeology with inferred ages of giant stars from Gaia spectra
- Open Access
- 2 January 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 546 (1)
Towards Galactic archaeology with inferred ages of giant stars from Gaia spectra
In the era of Gaia, the accurate determination of stellar ages is transforming Galactic archaeology. We demonstrate the feasibility of inferring stellar ages from Gaia’s RVS spectra and the BP/RP (XP) spectrophotometric data, specifically for red giant branch and high-mass red clump stars. We successfully train two machine learning models, dubbed SIDRA: Stellar age Inference Derived from Gaia spectRA to predict the age. The SIDRA-RVS model uses the RVS spectra and SIDRA-XP the stellar parameters obtained from the XP spectra. Both models use BINGO, an APOGEE-derived stellar age as the training data. SIDRA-RVS estimates ages of stars whose age is around $\tau _\mathrm{BINGO}=10$Gyr with a standard deviation of residuals of $\sim$ 0.12 dex in the unseen test data set, while SIDRA-XP achieves higher precision with residuals $\sim$ 0.064 dex for stars around $\tau _\mathrm{BINGO}=10$ Gyr. Since SIDRA-XP outperforms SIDRA-RVS, we apply SIDRA-XP to analyse the ages for 2218 154 stars. This allowed us to map the chronological and chemical properties of Galactic disc stars, reproducing the known distinct features such as the Gaia-Sausage-Enceladus merger and a potential gas-rich interaction event linked to the first infall of the Sagittarius dwarf galaxy. This study demonstrates that machine learning techniques applied to Gaia’s spectra can provide valuable individual age information, particularly for giant stars, thereby enhancing our understanding of the Milky Way’s formation and evolution.- 2 January 2026
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 546 (1)
Open AccessTowards Galactic archaeology with inferred ages of giant stars from Gaia spectra
In the era of Gaia, the accurate determination of stellar ages is transforming Galactic archaeology. We demonstrate the feasibility of inferring stellar ages from Gaia’s RVS spectra and the BP/RP (XP) spectrophotometric data, specifically for red giant branch and high-mass red clump stars. We successfully train two machine learning models, dubbed SIDRA: Stellar age Inference Derived from Gaia spectRA to predict the age. The SIDRA-RVS model uses the RVS spectra and SIDRA-XP the stellar parameters obtained from the XP spectra. Both models use BINGO, an APOGEE-derived stellar age as the training data. SIDRA-RVS estimates ages of stars whose age is around $\tau _\mathrm{BINGO}=10$Gyr with a standard deviation of residuals of $\sim$ 0.12 dex in the unseen test data set, while SIDRA-XP achieves higher precision with residuals $\sim$ 0.064 dex for stars around $\tau _\mathrm{BINGO}=10$ Gyr. Since SIDRA-XP outperforms SIDRA-RVS, we apply SIDRA-XP to analyse the ages for 2218 154 stars. This allowed us to map the chronological and chemical properties of Galactic disc stars, reproducing the known distinct features such as the Gaia-Sausage-Enceladus merger and a potential gas-rich interaction event linked to the first infall of the Sagittarius dwarf galaxy. This study demonstrates that machine learning techniques applied to Gaia’s spectra can provide valuable individual age information, particularly for giant stars, thereby enhancing our understanding of the Milky Way’s formation and evolution. - Open Access
Orbital migration and heating history of the Galactic disc: a transition between the bimodal discs
Stellar orbits in the Galactic disc evolve from their birth to the current shape through both radial migration and dynamical heating. The history of their secular evolution is imprinted in the current kinematics and age-metallicity distribution. We construct a chrono-chemo-dynamical model of the disc, incorporating inside-out growth, metallicity evolution, radial migration, and heating to fit the observed age-metallicity-kinematics distribution of LAMOST subgiant stars in both the low and high-$α$ disc. By modelling all distribution parameters with spline fitting, we present the first non-parametric stellar migration and heating history of the Galaxy. We determine the heating-to-migration ratio, the ratio of the root-mean-square changes in radial/vertical and azimuthal actions, to be $\approx0.075$ for radial to azimuthal actions and $\approx0.015$ for vertical to azimuthal actions, implying a highly anisotropic diffusion in the action space. Furthermore, we identify a transition in radial migration efficiency coinciding with the transition moment of the bimodal disc, for which the radial migration was more efficient for the high-$α$ disc than for the low-$α$ disc. This transition may be attributed to two correlated scenarios: 1) a bar formation epoch accompanied by violent outward migration, and 2) a drop in the gas mass fraction in the disc when the low-$α$ disc began to form. These findings offer further constraints on the formation mechanisms of bimodal discs, favouring the downsizing scenario. We also briefly discuss the connection between our results and recent high-redshift observations. In addition to the secular evolution history, our model maps the Milky Way ISM metallicity gradient at different lookback times, which we find has only varied a little (in the range of $-0.07~\rm to~-0.10~dex/kpc$) since disc formation.Open AccessOrbital migration and heating history of the Galactic disc: a transition between the bimodal discs
Stellar orbits in the Galactic disc evolve from their birth to the current shape through both radial migration and dynamical heating. The history of their secular evolution is imprinted in the current kinematics and age-metallicity distribution. We construct a chrono-chemo-dynamical model of the disc, incorporating inside-out growth, metallicity evolution, radial migration, and heating to fit the observed age-metallicity-kinematics distribution of LAMOST subgiant stars in both the low and high-$α$ disc. By modelling all distribution parameters with spline fitting, we present the first non-parametric stellar migration and heating history of the Galaxy. We determine the heating-to-migration ratio, the ratio of the root-mean-square changes in radial/vertical and azimuthal actions, to be $\approx0.075$ for radial to azimuthal actions and $\approx0.015$ for vertical to azimuthal actions, implying a highly anisotropic diffusion in the action space. Furthermore, we identify a transition in radial migration efficiency coinciding with the transition moment of the bimodal disc, for which the radial migration was more efficient for the high-$α$ disc than for the low-$α$ disc. This transition may be attributed to two correlated scenarios: 1) a bar formation epoch accompanied by violent outward migration, and 2) a drop in the gas mass fraction in the disc when the low-$α$ disc began to form. These findings offer further constraints on the formation mechanisms of bimodal discs, favouring the downsizing scenario. We also briefly discuss the connection between our results and recent high-redshift observations. In addition to the secular evolution history, our model maps the Milky Way ISM metallicity gradient at different lookback times, which we find has only varied a little (in the range of $-0.07~\rm to~-0.10~dex/kpc$) since disc formation. - Open Access
Modelling the Milky Way's exoplanet population based on cosmological galaxy simulations
In this paper we aim to simulate realistic exoplanet populations across different regions of the MW by combining state-of-the-art cosmological simulations of our Galaxy with exoplanet formation models and observations. We model the exoplanet populations around single stars, using planet occurrence rates and multiplicity depending on stellar mass, metallicity, and planet type, and assign them physical parameters such as mass and orbital period. Focussing first on the solar vicinity, we find mostly metallicity-driven differences in the distributions of non-hosting and planet-hosting single stars. In our simulated solar neighbourhood, 52.5% of all planets are Earth-like (23% of them located in the Habitable Zone), 44% are super-Earths/Neptunes, and 3.5% are giant planets. A comparison with the census of Kepler exoplanets and candidates shows that, when taking into account the most relevant selection effects, we obtain a similar distribution of exoplanets compared to the observed population. However, we also detect significant differences in the exoplanet and host star distributions (e.g. more planets around F-type and red-giant stars compared to observations) that we attribute mostly to a too strong recent star formation and a too large disc scale height in the simulation, as well as to some caveats in our exoplanet population synthesis that will be addressed in future work. Extending our analysis to other regions of the simulated MW and to other simulated galaxies, we find that the relative percentages of planet types remain largely consistent as long as the simulated galaxy matches the morphology and mass of the MW. We have created a fast and flexible framework to produce exoplanet populations based on MW-like simulations that can easily be adapted to produce predictions for the yields of future exoplanet detection missions. (abridged)
Open AccessModelling the Milky Way's exoplanet population based on cosmological galaxy simulations
In this paper we aim to simulate realistic exoplanet populations across different regions of the MW by combining state-of-the-art cosmological simulations of our Galaxy with exoplanet formation models and observations. We model the exoplanet populations around single stars, using planet occurrence rates and multiplicity depending on stellar mass, metallicity, and planet type, and assign them physical parameters such as mass and orbital period. Focussing first on the solar vicinity, we find mostly metallicity-driven differences in the distributions of non-hosting and planet-hosting single stars. In our simulated solar neighbourhood, 52.5% of all planets are Earth-like (23% of them located in the Habitable Zone), 44% are super-Earths/Neptunes, and 3.5% are giant planets. A comparison with the census of Kepler exoplanets and candidates shows that, when taking into account the most relevant selection effects, we obtain a similar distribution of exoplanets compared to the observed population. However, we also detect significant differences in the exoplanet and host star distributions (e.g. more planets around F-type and red-giant stars compared to observations) that we attribute mostly to a too strong recent star formation and a too large disc scale height in the simulation, as well as to some caveats in our exoplanet population synthesis that will be addressed in future work. Extending our analysis to other regions of the simulated MW and to other simulated galaxies, we find that the relative percentages of planet types remain largely consistent as long as the simulated galaxy matches the morphology and mass of the MW. We have created a fast and flexible framework to produce exoplanet populations based on MW-like simulations that can easily be adapted to produce predictions for the yields of future exoplanet detection missions. (abridged)
- Open Access
- 25 September 2025
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 543 (3) , 2275-2291
Thick-to-thin disc transition and gas disc shrinking induced by the Gaia-Sausage-Enceladus merger
Understanding the Milky Way disc formation requires characterizing its structural and kinematic properties as functions of stellar age. Using red giant stars from the Apache Point Observatory Galactic Evolution Experiment DR17 and Gaia DR3, we model the age-dependent stellar kinematics with a quasi-isothermal distribution function and fit disc parameters as a function of age using non-parametric splines. We identify a transition from thick to thin disc populations around 10 Gyr ago. Stars older than this have short scale lengths ($\sim$1.7 kpc), typical of the thick disc, while younger stars exhibit increasing scale length with decreasing age, consistent with inside–out formation of the thin disc. This transition possibly coincides with the end of the starburst triggered by the Gaia-Sausage-Enceladus (GSE) merger. Stars formed around 10 Gyr ago exhibit a dip in scale length, even shorter than that of the thick disc. Comparison with an Auriga simulation suggests that this scale length dip reflects gas disc shrinking caused by the transition from a cold to hot gas accretion mode. We propose the following disc formation scenario: (1) the thick disc formed under cold-mode accretion; (2) the GSE merger triggered a starburst and increased the total mass of the Galaxy, causing the transition to hot-mode accretion; (3) rapid gas consumption led to temporary shrinking of the star-forming gas disc; and then (4) thin disc grows in an inside–out fashion, as the size of the star-forming gas disc grows via hot-mode smooth gas accretion.
- 25 September 2025
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 543 (3) , 2275-2291
Open AccessThick-to-thin disc transition and gas disc shrinking induced by the Gaia-Sausage-Enceladus merger
Understanding the Milky Way disc formation requires characterizing its structural and kinematic properties as functions of stellar age. Using red giant stars from the Apache Point Observatory Galactic Evolution Experiment DR17 and Gaia DR3, we model the age-dependent stellar kinematics with a quasi-isothermal distribution function and fit disc parameters as a function of age using non-parametric splines. We identify a transition from thick to thin disc populations around 10 Gyr ago. Stars older than this have short scale lengths ($\sim$1.7 kpc), typical of the thick disc, while younger stars exhibit increasing scale length with decreasing age, consistent with inside–out formation of the thin disc. This transition possibly coincides with the end of the starburst triggered by the Gaia-Sausage-Enceladus (GSE) merger. Stars formed around 10 Gyr ago exhibit a dip in scale length, even shorter than that of the thick disc. Comparison with an Auriga simulation suggests that this scale length dip reflects gas disc shrinking caused by the transition from a cold to hot gas accretion mode. We propose the following disc formation scenario: (1) the thick disc formed under cold-mode accretion; (2) the GSE merger triggered a starburst and increased the total mass of the Galaxy, causing the transition to hot-mode accretion; (3) rapid gas consumption led to temporary shrinking of the star-forming gas disc; and then (4) thin disc grows in an inside–out fashion, as the size of the star-forming gas disc grows via hot-mode smooth gas accretion.
- Open Access
- 14 August 2025
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 542 (3) , 1776-1790
GalactiKit: reconstructing mergers from z = 0 debris using simulation-based inference in Auriga
We present GalactiKit, a data-driven methodology for estimating the lookback infall time, stellar mass, halo mass, and mass ratio of the disrupted progenitors of Milky Way-like galaxies at the time of infall. GalactiKit uses simulation-based inference to extract the information on galaxy formation processes encoded in the Auriga cosmological magnetohydrodynamic (MHD) simulations of Milky Way-mass haloes to create a model that relates the properties of mergers to those of the corresponding merger debris at $z=0$. We investigate how well GalactiKit can reconstruct the merger properties given the dynamical, chemical, and the combined chemodynamical information of debris. For this purpose, three models were implemented considering the following properties of merger debris: (a) total energy and angular momentum, (b) iron-to-hydrogen and alpha-to-iron abundance ratios, and (c) a combination of all of these. We find that the kinematics of the debris can be used to trace the lookback time at which the progenitor was first accreted into the main halo. However, chemical information is necessary for inferring the stellar and halo masses of the progenitors. In both models (b) and (c), the stellar masses are predicted more accurately than the halo masses, which could be related to the scatter in the stellar mass–halo mass relation. Model (c) provides the most accurate predictions for the merger parameters, which suggests that combining chemical and dynamical data of debris can significantly improve the reconstruction of the Milky Way’s assembly history.- 14 August 2025
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 542 (3) , 1776-1790
Open AccessGalactiKit: reconstructing mergers from z = 0 debris using simulation-based inference in Auriga
We present GalactiKit, a data-driven methodology for estimating the lookback infall time, stellar mass, halo mass, and mass ratio of the disrupted progenitors of Milky Way-like galaxies at the time of infall. GalactiKit uses simulation-based inference to extract the information on galaxy formation processes encoded in the Auriga cosmological magnetohydrodynamic (MHD) simulations of Milky Way-mass haloes to create a model that relates the properties of mergers to those of the corresponding merger debris at $z=0$. We investigate how well GalactiKit can reconstruct the merger properties given the dynamical, chemical, and the combined chemodynamical information of debris. For this purpose, three models were implemented considering the following properties of merger debris: (a) total energy and angular momentum, (b) iron-to-hydrogen and alpha-to-iron abundance ratios, and (c) a combination of all of these. We find that the kinematics of the debris can be used to trace the lookback time at which the progenitor was first accreted into the main halo. However, chemical information is necessary for inferring the stellar and halo masses of the progenitors. In both models (b) and (c), the stellar masses are predicted more accurately than the halo masses, which could be related to the scatter in the stellar mass–halo mass relation. Model (c) provides the most accurate predictions for the merger parameters, which suggests that combining chemical and dynamical data of debris can significantly improve the reconstruction of the Milky Way’s assembly history. - Open Access
Investigating the eccentricity distribution of transiting, long-period giant planets
Eccentric giant planets are predicted to have acquired their eccentricity through two major mechanisms: the Kozai-Lidov effect or planet-planet scattering, but it is normally difficult to separate the two mechanisms and determine the true eccentricity origin for a given system. In this work, we focus on a sample of 92 transiting, long-period giant planets (TLGs) as part of an eccentricity distribution study for this planet population in order to understand their eccentricity origin. Using archival high-contrast imaging observations, public stellar catalogs, precise Gaia astrometry, and the NASA Exoplanet Archive database, we explored the eccentricity distribution correlation with different planet and host-star properties of our sample. We also homogeneously characterized the basic stellar properties for all 86 host-stars in our sample, including stellar age and metallicity. We found a correlation between eccentricity and stellar metallicity, where lower-metallicity stars ([Fe/H] <= 0.1) did not host any planets beyond e > 0.4, while higher-metallicity stars hosted planets across the entire eccentricity range. Interestingly, we found no correlation between the eccentricity distribution and the presence of stellar companions, indicating that planet-planet scattering is likely a more dominant mechanism than the Kozai-Lidov effect for TLGs. This is further supported by an anti-correlation trend found between planet multiplicity and eccentricity, as well as a lack of strong tidal dissipation effects for planets in our sample, which favor planet-planet scattering scenarios for the eccentricity origin.
Open AccessInvestigating the eccentricity distribution of transiting, long-period giant planets
Eccentric giant planets are predicted to have acquired their eccentricity through two major mechanisms: the Kozai-Lidov effect or planet-planet scattering, but it is normally difficult to separate the two mechanisms and determine the true eccentricity origin for a given system. In this work, we focus on a sample of 92 transiting, long-period giant planets (TLGs) as part of an eccentricity distribution study for this planet population in order to understand their eccentricity origin. Using archival high-contrast imaging observations, public stellar catalogs, precise Gaia astrometry, and the NASA Exoplanet Archive database, we explored the eccentricity distribution correlation with different planet and host-star properties of our sample. We also homogeneously characterized the basic stellar properties for all 86 host-stars in our sample, including stellar age and metallicity. We found a correlation between eccentricity and stellar metallicity, where lower-metallicity stars ([Fe/H] <= 0.1) did not host any planets beyond e > 0.4, while higher-metallicity stars hosted planets across the entire eccentricity range. Interestingly, we found no correlation between the eccentricity distribution and the presence of stellar companions, indicating that planet-planet scattering is likely a more dominant mechanism than the Kozai-Lidov effect for TLGs. This is further supported by an anti-correlation trend found between planet multiplicity and eccentricity, as well as a lack of strong tidal dissipation effects for planets in our sample, which favor planet-planet scattering scenarios for the eccentricity origin.