Three-dimensional radiative MHD simulations of near-surface convection in main sequence cool stars

  

Andrea Perdomo García has recently completed her PhD thesis, Three-dimensional radiative MHD simulations of near-surface convection in main sequence cool stars, supervised by Manuel Collados and myself. She successfully defended her work in front of international committee (Artemio Herrero (IAC), Oskar Steiner (KIS, Freiburg, in photo asking a question) and Maria Bergemann (MPI, Heidelberg)) and graduated cum laude.

Andrea’s thesis explored the atmospheres of cool main-sequence stars through realistic three-dimensional radiative MHD simulations with the MANCHA code. A major part of her work was devoted to one of the central challenges of such simulations: treating radiative transfer accurately while keeping it computationally feasible. She developed and tested opacity-binning strategies for stars ranging from F to M spectral types, investigated the increasingly important role of molecular opacity toward cooler stars, and studied how different opacity treatments affect atmospheric structure and radiative energy exchange.

She then applied these methods to 3D simulations of G2V, K0V, and M2V stars, exploring their convection and photospheric structure and extending the calculations to magnetic simulations in which fields generated through the Biermann battery were amplified by small-scale dynamo action.

It has been a true pleasure to work with Andrea over the years and to follow her development as a researcher. I wish her all the best in her future career and in her new position at the Max Planck Institute in Heidelberg. I am sure there are many interesting problems, simulations, and probably quite a few opacity tables still ahead!

Hydrodynamic simulations of cool stellar atmospheres with MANCHA

 

 

Perdomo García, A., Vitas, N., Khomenko, E., Collados, M., A&A, 688, A27 (2024)

arxiv

Three-dimensional simulations do much more than reproduce what we observe on stellar surfaces. They give us a unique laboratory in which we can experiment with stellar atmospheres: change the physical ingredients, follow their interaction, and understand how convection, radiation and atmospheric structure are connected.

In this paper, we used the MANCHA code to simulate the atmospheres of three cool main-sequence stars, G2V, K0V and M2V, and explored one particularly important ingredient: radiative transfer and opacity. We tested how different approximations to the enormously complex stellar opacity affect radiative energy exchange and, ultimately, the atmospheric structure.

The results show how increasingly important a realistic treatment of opacity becomes as we move toward cooler stars, where molecules start to dominate the radiative properties of the atmosphere.

Another important piece of Andrea’s PhD work, and another step toward realistic 3D modelling across the cool-star sequence. 

Moving to the dark side!

 

I have some big news to share. From March 1, I’ll be starting a new chapter and moving to the dark side of astrophysics: the substellar world of ultracool dwarfs!

This seems like a natural moment to make the move. With the Euclid ESA mission now opening an enormous new window on the faint and cool populations of the Milky Way, we are entering an era in which the number of known ultracool and substellar objects will increase substantially. Turning this wealth of new observations into physical understanding will require careful work on the modelling side.

In some ways, though, I am not moving very far at all. My main interests will remain what they have been for years: developing numerical agoriths for radiative transfer, equations of state, opacities, and numerical modelling of atmospheres. The difference is that I will now be applying them much further down the temperature scale. At the temperatures of ultracool dwarfs, molecules dominate the opacity, chemistry becomes increasingly tricky, condensates and clouds appear, and the coupling between chemistry, radiation and atmospheric structure gets the central stage in modelling problem. 

I’ll be joining the SUBSTELLAR project led by Eduardo Martín, one of the co-discoverers of Teide 1, the first confirmed brown dwarf, announced in 1995. The work is supported by the European Research Council through the ERC Advanced Grant SUBSTELLAR, devoted to pushing the frontier of substellar science with the Euclid mission.

A new wavelength regime, a new class of objects, and plenty of new physics — but, fortunately, still lots of opacities, radiative transfer and numerical challenges. I hope that this move will also allow me to bring some of my experience in stellar astrophysics across the boundary between these fields. I’m looking forward to seeing where it leads.

May the (gravity and Lorentz) Force be with me!