Showing 1–28 of 28 results for author: Witzke, V

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Stellar Surface Magnetic Fields Impact Limb Darkening

Authors: Nadiia M. Kostogryz, Alexander I. Shapiro, Veronika Witzke, Robert H. Cameron, Laurent Gizon, Natalie A. Krivova, Hans-G. Ludwig, Pierre F. L. Maxted, Sara Seager, Sami K. Solanki, Jeff Valenti

Abstract: Stars appear darker at their limbs than at their disk centers because at the limb we are viewing the higher and cooler layers of stellar photospheres. Limb darkening derived from state-of-the-art stellar atmosphere models systematically fails to reproduce recent transiting exoplanet light curves from the Kepler, TESS, and JWST telescopes -- stellar brightness obtained from measurements drops less… ▽ More Stars appear darker at their limbs than at their disk centers because at the limb we are viewing the higher and cooler layers of stellar photospheres. Limb darkening derived from state-of-the-art stellar atmosphere models systematically fails to reproduce recent transiting exoplanet light curves from the Kepler, TESS, and JWST telescopes -- stellar brightness obtained from measurements drops less steeply towards the limb than predicted by models. All previous models assumed atmosphere devoid of magnetic fields. Here we use our new stellar atmosphere models computed with the 3D radiative magneto-hydrodynamic code MURaM to show that small-scale concentration of magnetic fields on the stellar surface affect limb darkening at a level that allows us to explain the observations. Our findings provide a way forward to improve the determination of exoplanet radii and especially the transmission spectroscopy analysis for transiting planets, which relies on a very accurate description of stellar limb darkening from the visible through the infrared. Furthermore, our findings imply that limb darkening allows measuring the small-scale magnetic field on stars with transiting planets. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 19 pages, 6 figures This is the submission version of manuscript to Nature Astronomy

A roadmap to the efficient and robust characterization of temperate terrestrial planet atmospheres with JWST

Authors: TRAPPIST-1 JWST Community Initiative, :, Julien de Wit, René Doyon, Benjamin V. Rackham, Olivia Lim, Elsa Ducrot, Laura Kreidberg, Björn Benneke, Ignasi Ribas, David Berardo, Prajwal Niraula, Aishwarya Iyer, Alexander Shapiro, Nadiia Kostogryz, Veronika Witzke, Michaël Gillon, Eric Agol, Victoria Meadows, Adam J. Burgasser, James E. Owen, Jonathan J. Fortney, Franck Selsis, Aaron Bello-Arufe, Emeline Bolmont , et al. (57 additional authors not shown)

Abstract: Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets, which have been the favored targets of eight JWST Cycle 1 programs. While Cycle 1 observations have started to yield preliminary insights into the planets, they have… ▽ More Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets, which have been the favored targets of eight JWST Cycle 1 programs. While Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to constrain stellar heterogeneities and gather up to 2$\times$ more transits per JWST hour spent. We conclude that in such systems planets cannot be studied in isolation by small programs, thus large-scale community-supported programs should be supported to enable the efficient and robust exploration of terrestrial exoplanets in the JWST era. △ Less

Submitted 24 October, 2023; v1 submitted 24 October, 2023; originally announced October 2023.

Testing MURaM and MPS-ATLAS against the quiet solar spectrum

Authors: Veronika Witzke, Alexander I. Shapiro, Nadiia M. Kostogryz, Lucien Mauviard, Tanayveer S. Bhatia, Robert Cameron, Laurent Gizon, Damien Przybylski, Sami K. Solanki, Yvonne C. Unruh, Li Yue

Abstract: Three-dimensional (3D) radiative magnetohydrodynamics (MHD) simulations are the only way to model stellar atmospheres without any ad hoc parameterisations. Several 3D radiative MHD codes have achieved good quantitative agreement with observables for our Sun. We aim to validate the most up-to-date version of the MURaM code against well established quiet Sun measurements, in particular spatially ave… ▽ More Three-dimensional (3D) radiative magnetohydrodynamics (MHD) simulations are the only way to model stellar atmospheres without any ad hoc parameterisations. Several 3D radiative MHD codes have achieved good quantitative agreement with observables for our Sun. We aim to validate the most up-to-date version of the MURaM code against well established quiet Sun measurements, in particular spatially averaged measurements that are relevant for stellar studies. This validation extends the number of solar observables that MURaM can reproduce with high precision. It is also an essential condition for using MURaM to accurately calculate spectra of other cool stars. We simulate the solar upper convection zone and photosphere harbouring a small-scale-dynamo. Using time series of 3D snapshots we calculate the spectral irradiance, limb darkening and selected spectral lines, which we compare to observations. The computed observables agree well with the observations, in particular the limb darkening of the quiet Sun is reproduced remarkably well. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 9 Page, 7 figures, accepted in A&A

Journal ref: A&A 681, A81 (2024)

Spectral variability of photospheric radiation due to faculae II: Facular contrasts for cool main-sequence stars

Authors: Charlotte M. Norris, Yvonne C. Unruh, Veronika Witzke, Sami K. Solanki, Natalie A. Krivova, Alexander I. Shapiro, Robert Cameron, Benjamin Beeck

Abstract: Magnetic features on the surface of stars, such as spots and faculae, cause stellar spectral variability on time-scales of days and longer. For stars other than the Sun, the spectral signatures of faculae are poorly understood, limiting our ability to account for stellar pollution in exoplanet transit observations. Here we present the first facular contrasts derived from magnetoconvection simulati… ▽ More Magnetic features on the surface of stars, such as spots and faculae, cause stellar spectral variability on time-scales of days and longer. For stars other than the Sun, the spectral signatures of faculae are poorly understood, limiting our ability to account for stellar pollution in exoplanet transit observations. Here we present the first facular contrasts derived from magnetoconvection simulations for K0, M0 and M2 main-sequence stars and compare them to previous calculations for G2 main-sequence stars. We simulate photospheres and immediate subsurface layers of main-sequence spectral types between K0 and M2, with different injected vertical magnetic fields (0 G, 100 G, 300 G and 500 G) using MURaM, a 3D radiation-magnetohydrodynamics code. We show synthetic spectra and contrasts from the UV (300 nm) to the IR (10000 nm) calculated using the ATLAS9 radiative transfer code. The calculations are performed for nine viewing angles to characterise the facular radiation across the disc. The brightness contrasts of magnetic regions are found to change significantly across spectral type, wavelength and magnetic field strength, leading to the conclusion that accurate contrasts cannot be found by scaling solar values. This is due to features of different size, apparent structure and spectral brightness emerging in the presence of a given magnetic field for different spectral types. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 17 pages, 12 figures, to be published in MNRAS

MPS-ATLAS library of stellar model atmospheres and spectra

Authors: N. Kostogryz, A. I Shapiro, V. Witzke, D. Grant, H. R. Wakeford, K. B. Stevenson, S. K. Solanki, L. Gizon

Abstract: Stellar spectra contain a large amount of information about the conditions in stellar atmospheres. However, extracting this information is challenging and demands comprehensive numerical modelling. Here, we present stellar spectra synthesized using the recently developed state-of-the-art MPS-ATLAS code on a fine grid of stellar fundamental parameters. These calculations have been extensively valid… ▽ More Stellar spectra contain a large amount of information about the conditions in stellar atmospheres. However, extracting this information is challenging and demands comprehensive numerical modelling. Here, we present stellar spectra synthesized using the recently developed state-of-the-art MPS-ATLAS code on a fine grid of stellar fundamental parameters. These calculations have been extensively validated against solar and stellar observations and can be used for various astrophysical applications. The spectra are available at the Max Planck Digital Library (MPDL, https://edmond.mpdl.mpg.de/dataset.xhtml?persistentId=doi:10.17617/3.NJ56TR) and have been also incorporated into the ExoTiC-LD python package (https://github.com/Exo-TiC/ExoTiC-LD/) which returns stellar limb darkening coefficients used for the software package Exoplanet Timeseries Characterisation (ExoTic). △ Less

Submitted 5 March, 2023; originally announced March 2023.

Comments: 3 pages

Can 1D radiative equilibrium models of faculae be used for calculating contamination of transmission spectra?

Authors: Veronika Witzke, Alexander I. Shapiro, Nadiia M. Kostogryz, Robert Cameron, Benjamin V. Rackham, Sara Seager, Sami K. Solanki, Yvonne C. Unruh

Abstract: The reliable characterization of planetary atmospheres with transmission spectroscopy requires realistic modeling of stellar magnetic features, since features that are attributable to an exoplanet atmosphere could instead stem from the host star's magnetic activity. Current retrieval algorithms for analysing transmission spectra rely on intensity contrasts of magnetic features from 1D radiative-co… ▽ More The reliable characterization of planetary atmospheres with transmission spectroscopy requires realistic modeling of stellar magnetic features, since features that are attributable to an exoplanet atmosphere could instead stem from the host star's magnetic activity. Current retrieval algorithms for analysing transmission spectra rely on intensity contrasts of magnetic features from 1D radiative-convective models. However, magnetic features, especially faculae, are not fully captured by such simplified models. Here we investigate how well such 1D models can reproduce 3D facular contrasts, taking a G2V star as an example. We employ the well established radiative magnetohydrodynamic code MURaM to obtain three-dimensional simulations of the magneto-convection and photosphere harboring a local small-scale-dynamo. Simulations without additional vertical magnetic fields are taken to describe the quiet solar regions, while simulations with initially 100 G, 200 G and 300 G vertical magnetic fields are used to represent different magnetic activity levels. Subsequently, the spectra emergent from the MURaM cubes are calculated with the MPS-ATLAS radiative transfer code. We find that the wavelength dependence of facular contrast from 1D radiative-convective models cannot reproduce facular contrasts obtained from 3D modeling. This has far reaching consequences for exoplanet characterization using transmission spectroscopy, where accurate knowledge of the host star is essential for unbiased inferences of the planetary atmospheric properties. △ Less

Submitted 5 November, 2022; originally announced November 2022.

Comments: 7 pages, 2 figures, submitted to APJL

Small-scale dynamo in cool main sequence stars. II. The effect of metallicity

Authors: V. Witzke, H. B. Duehnen, A. I. Shapiro, D. Przybylski, T. S. Bhatia, R. Cameron, S. K. Solanki

Abstract: All cool main sequence stars including our Sun are thought to have magnetic fields. Observations of the Sun revealed that even in quiet regions small-scale turbulent magnetic fields are present. Simulations further showed that such magnetic fields affect the subsurface and photospheric structure, and thus the radiative transfer and emergent flux. Since small-scale turbulent magnetic fields on othe… ▽ More All cool main sequence stars including our Sun are thought to have magnetic fields. Observations of the Sun revealed that even in quiet regions small-scale turbulent magnetic fields are present. Simulations further showed that such magnetic fields affect the subsurface and photospheric structure, and thus the radiative transfer and emergent flux. Since small-scale turbulent magnetic fields on other stars cannot be directly observed, it is imperative to study their effects on the near surface layers numerically. Until recently comprehensive three-dimensional simulations capturing the effect of small-scale turbulent magnetic fields only exists for the solar case. A series of investigations extending SSD simulations for other stars has been started. Here we aim to examine small-scale turbulent magnetic fields in stars of solar effective temperature but different metallicity. We investigate the properties of three-dimensional simulations of the magneto-convection in boxes covering the upper convection zone and photosphere carried out with the MURaM code for metallicity values of $ \rm M/H = \{-1.0, 0.0, 0.5\}$ with and without a small-scale-dynamo. We find that small-scale turbulent magnetic fields enhanced by a small-scale turbulent dynamo noticeably affect the subsurface dynamics and significantly change the flow velocities in the photosphere. Moreover, significantly stronger magnetic field strengths are present in the convection zone for low metallicity. Whereas, at the optical surface the averaged vertical magnetic field ranges from 64G for M/H = 0.5 to 85G for M/H = -1.0. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: 13 pages, 18 figures, submitted to A&A

Journal ref: A&A 669, A157 (2023)

Stellar limb darkening. A new MPS-ATLAS library for Kepler, TESS, CHEOPS, and PLATO passbands

Authors: N. M. Kostogryz, V. Witzke, A. I. Shapiro, S. K. Solanki, P. F. L. Maxted, R. L. Kurucz, L. Gizon

Abstract: The detection of the first exoplanet paved the way into the era of transit photometry space missions with a revolutionary photometric precision that aim at discovering new exoplanetary systems around different types of stars. With this high precision, it is possible to derive very accurately the radii of exoplanets which is crucial for constraining their type and composition. However, it requires… ▽ More The detection of the first exoplanet paved the way into the era of transit photometry space missions with a revolutionary photometric precision that aim at discovering new exoplanetary systems around different types of stars. With this high precision, it is possible to derive very accurately the radii of exoplanets which is crucial for constraining their type and composition. However, it requires an accurate description of host stars, especially their center-to-limb variation of intensities (so called limb darkening) as it affects the planet-to-star radius ratio determination. We aim at improving the accuracy of limb darkening calculations for stars with a wide range of fundamental parameters. We used the recently developed 1D MPS-ATLAS code to compute model atmosphere structures and to synthesize stellar limb darkening on a very fine grid of stellar parameters. For the computations we utilized the most accurate information on chemical element abundances and mixing length parameters including convective overshoot. The stellar limb darkening was fitted using the two most accurate limb darkening laws: the power-2 and 4-parameters non-linear laws. We present a new extensive library of stellar model atmospheric structures, the synthesized stellar limb darkening curves, and the coefficients of parameterized limb-darkening laws on a very fine grid of stellar parameters in the Kepler, TESS, CHEOPS, and PLATO passbands. The fine grid allows overcoming the sizable errors introduced by the need to interpolate. Our computations of solar limb darkening are in a good agreement with available solar measurements at different view angles and wavelengths. Our computations of stellar limb darkening agree well with available measurements of Kepler stars. A new grid of stellar model structures, limb darkening and their fitted coefficients in different broad filters is provided in CDS. △ Less

Submitted 1 July, 2022; v1 submitted 14 June, 2022; originally announced June 2022.

Comments: 19 pages, 17 figures, accepted to Astronomy and Astrophysics. All tables are available in CDS

Journal ref: A&A 666, A60 (2022)

Small-scale Dynamo in Cool Stars: I. Changes in stratification and near-surface convection for main-sequence spectral types

Authors: Tanayveer S. Bhatia, Robert H. Cameron, Sami K. Solanki, Hardi Peter, Damien Przybylski, Veronika Witzke, Alexander Shapiro

Abstract: Some of the small-scale solar magnetic flux can be attributed to a small-scale dynamo (SSD) operating in the near-surface convection. The SSD fields have consequences for solar granular convection, basal flux, as well as chromospheric heating. A similar SSD mechanism is expected to be active in the near-surface convection of other cool main-sequence stars, but this has never been investigated. We… ▽ More Some of the small-scale solar magnetic flux can be attributed to a small-scale dynamo (SSD) operating in the near-surface convection. The SSD fields have consequences for solar granular convection, basal flux, as well as chromospheric heating. A similar SSD mechanism is expected to be active in the near-surface convection of other cool main-sequence stars, but this has never been investigated. We aim to investigate changes in stratification and convection due to inclusion of SSD fields for F3V, G2V, K0V and M0V spectral types in the near-surface convection. 3D magnetohydrodynamic (MHD) models of the four stellar boxes, covering the subsurface convection zone up to the lower photosphere in a small cartesian box, are studied using the \textit{MURaM} radiative-MHD simulation code. The SSD runs are compared against reference hydrodynamic runs. An SSD is found to efficiently produce magnetic field with energies ranging between 5\% to 80\% of the plasma kinetic energy at different depths. This ratio tends to be larger for larger $T_{\mathrm{eff}}$. The relative change in density and gas pressure stratification for the deeper convective layers due to SSD magnetic fields is negligible, except for the F-star. For the F-star, there is a substantial reduction in convective velocities due to Lorentz force feedback from magnetic fields, which, in turn, reduces the turbulent pressure. SSD in near-surface convection for cool main-sequence stars introduces small but significant changes in thermodynamic stratification (especially for the F-star) due to reduction in convective velocities. △ Less

Submitted 20 July, 2023; v1 submitted 31 May, 2022; originally announced June 2022.

Comments: Update (20.07.2023): an errata (appendix D) has been included to correct a small error in the derivation of turbulent pressure in appendix B. In addition, the text has been updated with the published version. 10 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 663, A166 (2022)

Chromospheric Extension of the MURaM Code

Authors: D. Przybylski, R. Cameron, S. K. Solanki, M. Rempel, J. Leenaarts, L. S. Anusha, V. Witzke, A. I. Shapiro

Abstract: Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombi… ▽ More Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 24 pages, 12 figures, to be published in A&A

Journal ref: A&A 664, A91 (2022)

Predictions of Astrometric Jitter for Sun-like Stars. II. Dependence on Inclination, Metallicity, and Active-Region Nesting

Authors: K. Sowmya, N. -E. Nèmec, A. I. Shapiro, E. Işık, V. Witzke, A. Mints, N. A. Krivova, S. K. Solanki

Abstract: Ultra-precise astrometry from the Gaia mission is expected to lead to astrometric detections of more than 20,000 exoplanets in our Galaxy. One of the factors that could hamper such detections is the astrometric jitter caused by the magnetic activity of the planet host stars. In our previous study, we modeled astrometric jitter for the Sun observed equator-on. In this work, we generalize our model… ▽ More Ultra-precise astrometry from the Gaia mission is expected to lead to astrometric detections of more than 20,000 exoplanets in our Galaxy. One of the factors that could hamper such detections is the astrometric jitter caused by the magnetic activity of the planet host stars. In our previous study, we modeled astrometric jitter for the Sun observed equator-on. In this work, we generalize our model and calculate the photocenter jitter as it would be measured by the Gaia and Small-JASMINE missions for stars with solar rotation rate and effective temperature, but with various values of the inclination angle of the stellar rotation axis. In addition, we consider the effect of metallicity and of nesting of active regions (i.e. the tendency of active regions to emerge in the vicinity of each other). We find that, while the jitter of stars observed equator-on does not have any long-term trends and can be easily filtered out, the photocenters of stars observed out of their equatorial planes experience systematic shifts over the course of the activity cycle. Such trends allow the jitter to be detected with continuous measurements, in which case it can interfere with planet detectability. An increase in the metallicity is found to increase the jitter caused by stellar activity. Active-region nesting can further enhance the peak-to-peak amplitude of the photocenter jitter to a level that could be detected by Gaia. △ Less

Submitted 3 July, 2021; originally announced July 2021.

Comments: 24 pages, 20 figures in the main body and 4 in the appendix, accepted for publication in ApJ

Journal ref: 2021, ApJ, 919, 94

MPS-ATLAS: A fast all-in-one code for synthesising stellar spectra

Authors: V. Witzke, A. I. Shapiro, M. Cernetic, R. V. Tagirov, N. M. Kostogryz, L. S. Anusha, Y. C. Unruh, S. K. Solanki, R. L. Kurucz

Abstract: Context. Stellar spectral synthesis is essential for various applications, ranging from determining stellar parameters to comprehensive stellar variability calculations. New observational resources as well as advanced stellar atmosphere modelling, taking three dimensional (3D) effects from radiative magnetohydrodynamics calculations into account, require a more efficient radiative transfer. Aims… ▽ More Context. Stellar spectral synthesis is essential for various applications, ranging from determining stellar parameters to comprehensive stellar variability calculations. New observational resources as well as advanced stellar atmosphere modelling, taking three dimensional (3D) effects from radiative magnetohydrodynamics calculations into account, require a more efficient radiative transfer. Aims. For accurate, fast and flexible calculations of opacity distribution functions (ODFs), stellar atmospheres and stellar spectra we developed an efficient code building on the well-established ATLAS9 code. The new code also paves the way for an easy and fast access to different elemental compositions in stellar calculations. Methods. For the generation of ODF tables we further developed the well-established DFSYNTHE code by implementing additional functionality, and a speed-up by employing a parallel computation scheme. In addition, the line lists used can be changed from Kurucz's recent lists. In particular, we implemented the VALD3 line list. Results. A new code, the Merged Parallelised Simplified ATLAS is presented. It combines the efficient generation of ODF, atmosphere modelling and spectral synthesis in local thermodynamic equilibrium, therefore being an all-in-one code. This all-in-one code provides more numerical functionality and is substantially faster compared to other available codes. The fully portable MPS-ATLAS code is validated against previous ATLAS9 calculations, the PHOENIX code calculations, and high quality observations. △ Less

Submitted 28 May, 2021; originally announced May 2021.

Comments: 22 pages, 22 figures, accepted for publication in A&A

Journal ref: A&A 653, A65 (2021)

Radiative transfer with opacity distribution functions: Application to narrow band filters

Authors: L. S. Anusha, A. I. Shapiro, V. Witzke, M. Cernetic, S. K. Solanki, L. Gizon

Abstract: Modelling of stellar radiative intensities in various spectral pass-bands plays an important role in stellar physics. At the same time the direct calculations of the high-resolution spectrum and then integrating it over the given spectral pass-band is computationally demanding due to the vast number of atomic and molecular lines. This is particularly so when employing three-dimensional (3D) models… ▽ More Modelling of stellar radiative intensities in various spectral pass-bands plays an important role in stellar physics. At the same time the direct calculations of the high-resolution spectrum and then integrating it over the given spectral pass-band is computationally demanding due to the vast number of atomic and molecular lines. This is particularly so when employing three-dimensional (3D) models of stellar atmospheres. To accelerate the calculations, one can employ approximate methods, e.g., the use of Opacity Distribution Functions (ODFs). Generally, ODFs provide a good approximation of traditional spectral synthesis i.e., computation of intensities through filters with strictly rectangular transmission function. However, their performance strongly deteriorates when the filter transmission noticeably changes within its pass-band, which is the case for almost all filters routinely used in stellar physics. In this context, the aims of this paper are a) to generalize the ODFs method for calculating intensities through filters with arbitrary transmission functions; b) to study the performance of the standard and generalized ODFs methods for calculating intensities emergent from 3D models of stellar atmosphere. For this purpose we use the newly-developed MPS-ATLAS radiative transfer code to compute intensities emergent 3D cubes simulated with the radiative magnetohydrodynamics code MURaM. The calculations are performed in the 1.5D regime, i.e., along many parallel rays passing through the simulated cube. We demonstrate that generalized ODFs method allows accurate and fast syntheses of spectral intensities and their centre-to-limb variations. △ Less

Submitted 28 April, 2021; originally announced April 2021.

Comments: 28 pages, 11 figures, accepted for publication in ApJS

Forward modelling of Kepler-band variability due to faculae and spots

Authors: Luke J. Johnson, Charlotte N. Norris, Yvonne C. Unruh, Sami K. Solanki, Natalie Krivova, Veronika Witzke, Alexander I. Shapiro

Abstract: Variability observed in photometric lightcurves of late-type stars (on timescales longer than a day) is a dominant noise source in exoplanet surveys and results predominantly from surface manifestations of stellar magnetic activity, namely faculae and spots. The implementation of faculae in lightcurve models is an open problem, with scaling typically based on spectra equivalent to hot stellar atmo… ▽ More Variability observed in photometric lightcurves of late-type stars (on timescales longer than a day) is a dominant noise source in exoplanet surveys and results predominantly from surface manifestations of stellar magnetic activity, namely faculae and spots. The implementation of faculae in lightcurve models is an open problem, with scaling typically based on spectra equivalent to hot stellar atmospheres or assuming a solar-derived facular contrast. We modelled rotational (single period) lightcurves of active G2, K0, M0 and M2 stars, with Sun-like surface distributions and realistic limb-dependent contrasts for faculae and spots. The sensitivity of lightcurve variability to changes in model parameters such as stellar inclination, feature area coverage, spot temperature, facular region magnetic flux density and active band latitudes is explored. For our lightcurve modelling approach we used actress, a geometrically accurate model for stellar variability. actress generates 2-sphere maps representing stellar surfaces and populates them with user-prescribed spot and facular region distributions. From this, lightcurves can be calculated at any inclination. Quiet star limb darkening and limb-dependent facular contrasts were derived from MURaM 3D magnetoconvection simulations using ATLAS9. 1D stellar atmosphere models were used for the spot contrasts. We applied actress in Monte-Carlo simulations, calculating lightcurve variability amplitudes in the Kepler band. We found that, for a given spectral type and stellar inclination, spot temperature and spot area coverage have the largest effect on variability of all simulation parameters. For a spot coverage of 1%, the typical variability of a solar-type star is around 2 parts-per-thousand. The presence of faculae clearly affects the mean brightness and lightcurve shape, but has relatively little influence on the variability. △ Less

Submitted 23 April, 2021; originally announced April 2021.

Comments: 17 pages, 16 figures. Accepted for publication in MNRAS 22 April 2021

Modeling stellar Ca II H & K emission variations. I. Effect of inclination on the S-index

Authors: K. Sowmya, A. I. Shapiro, V. Witzke, N. -E. Nèmec, T. Chatzistergos, K. L. Yeo, N. A. Krivova, S. K. Solanki

Abstract: The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to… ▽ More The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to be affected not only by the stellar intrinsic properties but also by the inclination angle of the stellar rotation axis. Until now such an inclination effect on S-index has remained largely unexplored. To fill this gap, we develop a physics-based model to calculate S-index, focusing on the Sun. Using the distributions of solar magnetic features derived from observations together with Ca II H & K spectra synthesized in non-local thermodynamic equilibrium, we validate our model by successfully reconstructing the observed variations of solar S-index over four activity cycles. Further, using the distribution of magnetic features over the visible solar disk obtained from surface flux transport simulations, we obtain S-index time series dating back to 1700 and investigate the effect of inclination on S-index variability, both on the magnetic activity cycle and the rotational timescales. We find that when going from an equatorial to a pole-on view, the amplitude of S-index variations decreases weakly on the activity cycle timescale and strongly on the rotational timescale (by about 22% and 81%, respectively, for a cycle of intermediate strength). The absolute value of S-index depends only weakly on the inclination. We provide analytical expressions that model such dependencies. △ Less

Submitted 25 March, 2021; originally announced March 2021.

Comments: 23 pages, 13 figures in the main body and 2 in the appendix, accepted for publication in ApJ

The Dimmest State of the Sun

Authors: K. L. Yeo, S. K. Solanki, N. A. Krivova, M. Rempel, L. S. Anusha, A. I. Shapiro, R. V. Tagirov, V. Witzke

Abstract: How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm… ▽ More How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm-2 to a rise of 6.3 Wm-2. Consequently, the exact contribution by solar forcing to the rise in global temperatures over the past centuries remains inconclusive. Adopting a novel approach based on state-of-the-art solar imagery and numerical simulations, we establish the TSI level of the Sun when it is in its least-active state to be 2.0 +/- 0.7 Wm-2 below the 2019 level. This means TSI could not have risen since the Maunder minimum by more than this amount, thus restricting the possible role of solar forcing in global warming. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: Supporting information: https://agupubs.onlinelibrary.wiley.com/action/downloadSupplement?doi=10.1029%2F2020GL090243&file=grl61316-sup-0001-Text_SI-S01.pdf

Journal ref: Geophysical Research Letters, 47, e2020GL090243 (2020)

Where have all the solar-like stars gone? Rotation period detectability at various inclinations and metallicities

Authors: Timo Reinhold, Alexander I. Shapiro, Veronika Witzke, Nina-E. Nèmec, Emre Işık, Sami K. Solanki

Abstract: The plethora of photometric data collected by the Kepler space telescope has promoted the detection of tens of thousands of stellar rotation periods. However, these periods are not found to an equal extent among different spectral types. Interestingly, early G-type stars with near-solar rotation periods are strongly underrepresented among those stars with known rotation periods. In this study we i… ▽ More The plethora of photometric data collected by the Kepler space telescope has promoted the detection of tens of thousands of stellar rotation periods. However, these periods are not found to an equal extent among different spectral types. Interestingly, early G-type stars with near-solar rotation periods are strongly underrepresented among those stars with known rotation periods. In this study we investigate whether the small number of such stars can be explained by difficulties in the period determination from photometric time series. For that purpose, we generate model light curves of early G-type stars with solar rotation periods for different inclination angles, metallicities and (magnitude-dependent) noise levels. We find that the detectability is determined by the predominant type of activity (i.e. spot or faculae domination) on the surface, which defines the degree of irregularity of the light curve, and further depends on the level of photometric noise. These two effects significantly complicate the period detection and explain the lack of solar-like stars with known near-solar rotation periods. We conclude that the rotation periods of the majority of solar-like stars with near-solar rotation periods remain undetected to date. Finally, we promote the use of new techniques to recover more periods of near-solar rotators. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: 14 pages, 7 figures, accepted for publication in ApJL

Power spectrum of turbulent convection in the solar photosphere

Authors: L. Yelles Chaouche, R. H. Cameron, S. K. Solanki, T. L. Riethmüller, L. S. Anusha, V. Witzke, A. I. Shapiro, P. Barthol, A. Gandorfer, L. Gizon, J. Hirzberger, M. van Noort, J. Blanco Rodríguez, J. C. Del Toro Iniesta, D. Orozco Suárez, W. Schmidt, V. Martínez Pillet, M. Knölker

Abstract: The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical… ▽ More The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law. △ Less

Submitted 18 October, 2020; originally announced October 2020.

Comments: 10 pages, 7 figures, Accepted in A and A

Journal ref: A&A 644, A44 (2020)

Effect of metallicity on the detectability of rotational periods in solar-like stars

Authors: V. Witzke, T. Reinhold, A. I. Shapiro, N. A. Krivova, S. K. Solanki

Abstract: Accurate determination of stellar rotation periods is important for estimating stellar ages as well as for understanding stellar activity and evolution. While for about thirty thousand stars in the Kepler field rotation periods can be determined, there are over hundred thousand stars, especially with low photometric variability and irregular pattern of variations, for which rotational periods are… ▽ More Accurate determination of stellar rotation periods is important for estimating stellar ages as well as for understanding stellar activity and evolution. While for about thirty thousand stars in the Kepler field rotation periods can be determined, there are over hundred thousand stars, especially with low photometric variability and irregular pattern of variations, for which rotational periods are unknown. Here, we investigate the effect of metallicity on the detectability of rotation periods. This is done by synthesising light curves of hypothetical stars, which are identical to our Sun, with the exception of the metallicity. These light curves are then used as an input to the period determination algorithms. We find that the success rate for recovering the rotation signal has a minimum close to the solar metallicity value. This can be explained by the compensation effect of facular and spot contributions. In addition, selecting solar-like stars with near-solar effective temperature, near solar photometric variability, and with metallicity between M/H = -0.35 and M/H = 0.35 from the Kepler sample, we analyse the fraction of stars for which rotational periods have been detected as a function of metallicity. In agreement with our theoretical estimate we found a local minimum for the detection fraction close to the solar metallicity. We further report rotation periods of 87~solar-like Kepler stars for the first time. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 9 pages, 7 figures, accepted in A&A

Journal ref: A&A 634, L9 (2020)

Opacity distribution functions for stellar spectra synthesis

Authors: M. Cernetic, A. I. Shapiro, V. Witzke, N. A. Krivova, S. K. Solanki, R. V. Tagirov

Abstract: Contex. Stellar spectra synthesis is essential for the characterization of potential planetary hosts. In addition, comprehensive stellar variability calculations with fast radiative transfer are needed to disentangle planetary transits from stellar magnetically-driven variability. The planet-hunting space telescopes, such as CoRoT, Kepler, and TESS will bring vast quantities of data, rekindling th… ▽ More Contex. Stellar spectra synthesis is essential for the characterization of potential planetary hosts. In addition, comprehensive stellar variability calculations with fast radiative transfer are needed to disentangle planetary transits from stellar magnetically-driven variability. The planet-hunting space telescopes, such as CoRoT, Kepler, and TESS will bring vast quantities of data, rekindling the interest in fast calculations of the radiative transfer. Aims. We revisit the Opacity Distribution Functions (ODF) approach routinely applied to speedup stellar spectral synthesis. To achieve a considerable speedup relative to the current state-of-the-art, we further optimize the approach and search for the best ODF configuration. Furthermore, we generalize the ODF approach for fast calculations of flux in various filters often used in stellar observations. Methods. In a parameter-sweep-fashion, we generated ODF in the spectral range from UV to IR with different setups. The most accurate ODF configuration for each spectral interval was determined. For calculations of the radiative fluxes through filters we adapted the wavelength grid based on the transmission curve, whereafter the normal ODF procedure was performed. Results. Our optimum ODF configuration allows for a three fold speedup, compared to the previously used ODF configurations. The ODF generalization to calculate fluxes through filters results in a speedup of more than two orders of magnitude. △ Less

Submitted 24 June, 2019; v1 submitted 7 June, 2019; originally announced June 2019.

Comments: 8 pages, 11 figures

Journal ref: A&A 627, A157 (2019)

Evolution and characteristics of forced shear flows in polytropic atmospheres: Large and small Péclet number regimes

Authors: V. Witzke, L. J. Silvers, B. Favier

Abstract: Complex mixing and magnetic field generation occurs within stellar interiors particularly where there is a strong shear flow. To obtain a comprehensive understanding of these processes, it is necessary to study the complex dynamics of shear regions. Due to current observational limitations, it is necessary to investigate the inevitable small-scale dynamics via numerical calculations. Here, we exam… ▽ More Complex mixing and magnetic field generation occurs within stellar interiors particularly where there is a strong shear flow. To obtain a comprehensive understanding of these processes, it is necessary to study the complex dynamics of shear regions. Due to current observational limitations, it is necessary to investigate the inevitable small-scale dynamics via numerical calculations. Here, we examine direct numerical calculations of a local model of unstable shear flows in a compressible polytropic fluid primarily in a two-dimensional domain, where we focus on determining how key parameters affect the global properties and characteristics of the resulting saturated turbulent phase. We consider the effect of varying both the viscosity and the thermal diffusivity on the non-linear evolution. Moreover, our main focus is to understand the global properties of the saturated phase, in particular estimating for the first time the spread of the shear region from an initially hyperbolic tangent velocity profile. We find that the vertical extent of the mixing region in the saturated regime is generally determined by the initial Richardson number of the system. Further, the characteristic quantities of the turbulence, i.e. typical length-scale and the root-mean-square velocity are found to depend on both the Richardson number, and the thermal diffusivity. Finally, we present our findings of our investigation into saturated flows of a `secular' shear instability in the low Péclet number regime with large Richardson numbers. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 14 pager, 9 figures, accepted in MNRAS

From Solar to Stellar Brightness Variations: The Effect of Metallicity

Authors: V. Witzke, A. I. Shapiro, S. K. Solanki, N. A. Krivova, W. Schmutz

Abstract: Context. Comparison studies of Sun-like stars with the Sun suggest an anomalously low photometric variability of the Sun compared to Sun-like stars with similar magnetic activity. Comprehensive understanding of stellar variability is needed, to find a physical reasoning for this observation. Aims. We investigate the effect of metallicity and effective temperature on the photometric brightness chan… ▽ More Context. Comparison studies of Sun-like stars with the Sun suggest an anomalously low photometric variability of the Sun compared to Sun-like stars with similar magnetic activity. Comprehensive understanding of stellar variability is needed, to find a physical reasoning for this observation. Aims. We investigate the effect of metallicity and effective temperature on the photometric brightness change of Sun-like stars seen at different inclinations. The considered range of fundamental stellar parameters is sufficiently small so the stars, investigated here, still count as Sun-like or even as solar twins. Methods. To model the brightness change of stars with solar magnetic activity, we extend a well established model of solar brightness variations, SATIRE (which stands for Spectral And Total Irradiance Reconstruction), which is based on solar spectra, to stars with different fundamental parameters. For that we calculate stellar spectra for different metallicities and effective temperature using the radiative transfer code ATLAS9. Results. We show that even a small change (e.g. within the observational error range) of metallicity or effective temperature significantly affects the photometric brightness change compared to the Sun. We find that for Sun-like stars, the amplitude of the brightness variations obtained for Strömgren (b + y)/2 reaches a local minimum for fundamental stellar parameters close to the solar metallicity and effective temperature. Moreover, our results show that the effect of inclination decreases for metallicity values greater than the solar metallicity. Overall, we find that an exact determination of fundamental stellar parameters is crucially important for understanding stellar brightness changes. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 12 pages, 12 figures, accepted in A&A

Journal ref: A&A 619, A146 (2018)

Response of solar irradiance to sunspot area variations

Authors: T. Dudok de Wit, G. Kopp, A. Shapiro, V. Witzke, M. Kretzschmar

Abstract: One of the important open questions in solar irradiance studies is whether long-term variability (i.e. on timescales of years and beyond) can be reconstructed by means of models that describe short-term variability (i.e. days) using solar proxies as inputs. Preminger and Walton (2005, GRL, 32, 14109) showed that the relationship between spectral solar irradiance and proxies of magnetic-flux emerge… ▽ More One of the important open questions in solar irradiance studies is whether long-term variability (i.e. on timescales of years and beyond) can be reconstructed by means of models that describe short-term variability (i.e. days) using solar proxies as inputs. Preminger and Walton (2005, GRL, 32, 14109) showed that the relationship between spectral solar irradiance and proxies of magnetic-flux emergence, such as the daily sunspot area, can be described in the framework of linear system theory by means of the impulse response. We significantly refine that empirical model by removing spurious solar-rotational effects and by including an additional term that captures long-term variations. Our results show that long-term variability cannot be reconstructed from the short-term response of the spectral irradiance, which cautions the extension of solar proxy models to these timescales. In addition, we find that the solar response is nonlinear in such a way that cannot be corrected simply by applying a rescaling to sunspot area. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Comments: 18 pages. Updated version of an article published in The Astrophysical Journal

Journal ref: The Astrophysical Journal (2018) 853:197

The influence of metallicity on stellar differential rotation and magnetic activity

Authors: Christoffer Karoff, Travis S. Metcalfe, Angela R. G. Santos, Benjamin T. Montet, Howard Isaacson, Veronika Witzke, Alexander I. Shapiro, Savita Mathur, Guy R. Davies, Mikkel N. Lund, Rafael A. Garcia, Allan S. Brun, David Salabert, Pedro P. Avelino, Jennifer van Saders, Ricky Egeland, Margarida S. Cunha, Tiago L. Campante, William J. Chaplin, Natalie Krivova, Sami K. Solanki, Maximilian Stritzinger, Mads F. Knudsen

Abstract: Observations of Sun-like stars over the last half-century have improved our understanding of how magnetic dynamos, like that responsible for the 11-year solar cycle, change with rotation, mass and age. Here we show for the first time how metallicity can affect a stellar dynamo. Using the most complete set of observations of a stellar cycle ever obtained for a Sun-like star, we show how the solar a… ▽ More Observations of Sun-like stars over the last half-century have improved our understanding of how magnetic dynamos, like that responsible for the 11-year solar cycle, change with rotation, mass and age. Here we show for the first time how metallicity can affect a stellar dynamo. Using the most complete set of observations of a stellar cycle ever obtained for a Sun-like star, we show how the solar analog HD 173701 exhibits solar-like differential rotation and a 7.4-year activity cycle. While the duration of the cycle is comparable to that generated by the solar dynamo, the amplitude of the brightness variability is substantially stronger. The only significant difference between HD 173701 and the Sun is its metallicity, which is twice the solar value. Therefore, this provides a unique opportunity to study the effect of the higher metallicity on the dynamo acting in this star and to obtain a comprehensive understanding of the physical mechanisms responsible for the observed photometric variability. The observations can be explained by the higher metallicity of the star, which is predicted to foster a deeper outer convection zone and a higher facular contrast, resulting in stronger variability. △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: Submitted to ApJ

Mean Flow Evolution of Saturated Forced Shear Flows in Polytropic Atmospheres

Authors: V. Witzke, L. J. Silvers

Abstract: In stellar interiors shear flows play an important role in many physical processes. So far helioseismology provides only large-scale measurements, and so the small-scale dynamics remains insufficiently understood. To draw a connection between observations and three-dimensional DNS of shear driven turbulence, we investigate horizontally averaged profiles of the numerically obtained mean state. We f… ▽ More In stellar interiors shear flows play an important role in many physical processes. So far helioseismology provides only large-scale measurements, and so the small-scale dynamics remains insufficiently understood. To draw a connection between observations and three-dimensional DNS of shear driven turbulence, we investigate horizontally averaged profiles of the numerically obtained mean state. We focus here on just one of the possible methods that can maintain a shear flow, namely the average relaxation method. We show that although some systems saturate by restoring linear marginal stability this is not a general trend. Finally, we discuss the reason that the results are more complex than expected. △ Less

Submitted 30 November, 2016; originally announced November 2016.

Comments: 5 pages, 1 figure, Astrofluid Conference Proceedings

Evolution of forced shear flows in polytropic atmospheres: A comparison of forcing methods and energetics

Authors: V. Witzke, L. J. Silvers, B. Favier

Abstract: Shear flows are ubiquitous in astrophysical objects including planetary and stellar interiors, where their dynamics can have significant impact on thermo-chemical processes. Investigating the complex dynamics of shear flows requires numerical calculations that provide a long time evolution of the system. To achieve a sufficiently long lifetime in a local numerical model the system has to be forced… ▽ More Shear flows are ubiquitous in astrophysical objects including planetary and stellar interiors, where their dynamics can have significant impact on thermo-chemical processes. Investigating the complex dynamics of shear flows requires numerical calculations that provide a long time evolution of the system. To achieve a sufficiently long lifetime in a local numerical model the system has to be forced externally. However, at present, there exist several different forcing methods to sustain large-scale shear flows in local models. In this paper we examine and compare various methods used in the literature in order to resolve their respective applicability and limitations. These techniques are compared during the exponential growth phase of a shear flow instability, such as the Kelvin-Helmholtz (KH) instability, and some are examined during the subsequent non-linear evolution. A linear stability analysis provides reference for the growth rate of the most unstable modes in the system and a detailed analysis of the energetics provides a comprehensive understanding of the energy exchange during the system's evolution. Finally, we discuss the pros and cons of each forcing method and their relation with natural mechanisms generating shear flows. △ Less

Submitted 4 August, 2016; originally announced August 2016.

Comments: 14 pages, 9 figures, accepted in MNRAS, http://mnras.oxfordjournals.org/content/early/2016/08/03/mnras.stw1925.abstract?keytype=ref&ijkey=ExcZTzRyI0Gztdd

Shear instabilities in a fully compressible polytropic atmosphere

Authors: V. Witzke, L. J. Silvers, B. Favier

Abstract: Shear flows have an important impact on the dynamics in an assortment of different astrophysical objects including accreditation discs and stellar interiors. Investigating shear flow instabilities in a polytropic atmosphere provides a fundamental understanding of the motion in stellar interiors where turbulent motions, mixing processes, as well as magnetic field generation takes place. Here, a lin… ▽ More Shear flows have an important impact on the dynamics in an assortment of different astrophysical objects including accreditation discs and stellar interiors. Investigating shear flow instabilities in a polytropic atmosphere provides a fundamental understanding of the motion in stellar interiors where turbulent motions, mixing processes, as well as magnetic field generation takes place. Here, a linear stability analysis for a fully compressible fluid in a two-dimensional Cartesian geometry is carried out. Our study focuses on determining the critical Richardson number for different Mach numbers and the destabilising effects of high thermal diffusion. We find that there is a deviation of the predicted stability threshold for moderate Mach number flows along with a significant effect on the growth rate of the linear instability for small Péclet numbers. We show that in addition to a Kelvin-Helmholtz instability a Holmboe instability can appear and we discuss the implication of this in stellar interiors. △ Less

Submitted 10 March, 2015; originally announced March 2015.

Comments: 8 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 577, A76 (2015)