Showing 1–50 of 307 results for author: Shapiro, I

Low-latitude magnetic flux emergence on rapidly rotating solar-type stars

Authors: Emre Işık, Sami K. Solanki, Robert H. Cameron, Alexander I. Shapiro

Abstract: Besides a dense coverage of their high latitudes by starspots, rapidly rotating cool stars also display low-latitude spots in Doppler images, although generally with a lower coverage. In contrast, flux emergence models of fast-rotating stars predict strong poleward deflection of radially rising magnetic flux as the Coriolis effect dominates over buoyancy, leaving a spot-free band around the equato… ▽ More Besides a dense coverage of their high latitudes by starspots, rapidly rotating cool stars also display low-latitude spots in Doppler images, although generally with a lower coverage. In contrast, flux emergence models of fast-rotating stars predict strong poleward deflection of radially rising magnetic flux as the Coriolis effect dominates over buoyancy, leaving a spot-free band around the equator. To resolve this discrepancy, we consider a flux tube near the base of the convection zone in a solar-type star rotating eight times faster than the Sun, considering field intensification by weak-tube explosions. For the intensification to continue into to the buoyancy-dominated regime, the upper convection zone must have a significantly steeper temperature gradient than in the Sun, by a factor that is comparable with that found in 3D simulations of rotating convection. Within the hypothesis that stellar active regions stem from the base of the convection zone, flux emergence between 1-20 degree latitudes requires highly supercritical field strengths of up to 500 kG in rapidly rotating stars. These field strengths require explosions of 100-kG tubes within the convection zone, compatible with reasonable values of the superadiabatic temperature gradient associated with the more rapid rotation. △ Less

Submitted 18 October, 2024; originally announced October 2024.

Comments: 12 pages, 9 figures, accepted for publication in ApJ

On the renormalization of massive vector field theory coupled to scalar in curved space-time

Authors: Ioseph L. Buchbinder, Públio Rwany B. R. do Vale, Guilherme Y. Oyadomari, Ilya L. Shapiro

Abstract: We consider the renormalization of massive vector field interacting with charged scalar field in curved spacetime. Starting with the theory minimally coupled to external gravity and using the formulations with and without Stückelberg fields, we show that the longitudinal mode of vector field is completely decoupled and the remaining theory of transverse vector field is renormalizable by power coun… ▽ More We consider the renormalization of massive vector field interacting with charged scalar field in curved spacetime. Starting with the theory minimally coupled to external gravity and using the formulations with and without Stückelberg fields, we show that the longitudinal mode of vector field is completely decoupled and the remaining theory of transverse vector field is renormalizable by power counting. The formal arguments based on the covariance and power counting indicate that multiplicative renormalizability of the interacting theory may require introducing two non-minimal terms linear in Ricci tensor in the vector sector. Nevertheless, a more detailed analysis shows that these non-minimal terms violate the decoupling of the longitudinal mode and are prohibited. As a verification of general arguments, we derive the one-loop divergences in the minimal massive scalar QED, using Stückelberg procedure and the heat-kernel technique. The theory without non-minimal terms proves one-loop renormalizable and admits the renormalization group equations for all the running parameters in the scalar and vector sectors. One-loop beta functions do not depend on the gauge fixing and can be used to derive the effective potential. △ Less

Submitted 1 October, 2024; originally announced October 2024.

Comments: 26 pages, 2 figures

Bound states of massive complex ghosts in superrenormalizable quantum gravity theories\

Authors: Manuel Asorey, Gastao Krein, Ilya L. Shapiro

Abstract: One of the remarkable differences between renormalizable quantum gravity with four-derivative action and its superrenormalizable polynomial generalizations is that the latter admit a more sophisticated particle mass spectrum. Already in the simplest superrenormalizable case, the theory has a six-derivative Lagrangian, admitting either a real or complex spectrum of masses. In the case of a real spe… ▽ More One of the remarkable differences between renormalizable quantum gravity with four-derivative action and its superrenormalizable polynomial generalizations is that the latter admit a more sophisticated particle mass spectrum. Already in the simplest superrenormalizable case, the theory has a six-derivative Lagrangian, admitting either a real or complex spectrum of masses. In the case of a real spectrum, there are the graviton, massive unphysical ghosts, and normal particles with masses exceeding the ones of the ghosts. It is also possible to have pairs of complex conjugate massive ghost-like particles. We show that in both cases, these theories do not admit a Källén-Lehmann representation and do not satisfy the positivity criterium of consistency in terms of the fields associated to those particles. In the main part of the work, using a relatively simple Euclidean scalar toy model, we show that the theory with complex spectrum forms bound states confining unphysical massive excitations into a normal composite particle. Finally, we discuss the cosmological implications of such a ghost confinement. △ Less

Submitted 7 September, 2024; v1 submitted 29 August, 2024; originally announced August 2024.

Comments: 20 pages, 1 figure. Replaced version: revised text and references, change in the title

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

Semiclassical bounce with strong minimal assumptions

Authors: Wagno Cesar e Silva, Ilya L. Shapiro

Abstract: We explore the possibility of avoiding cosmological singularity with a bounce solution in the early Universe. The main finding is that simple and well-known semiclassical correction, which describes the mixing of radiation and gravity in the effective action, may provide an analytic solution with a bounce. The solution requires a positive beta function for the total radiation term and the contract… ▽ More We explore the possibility of avoiding cosmological singularity with a bounce solution in the early Universe. The main finding is that simple and well-known semiclassical correction, which describes the mixing of radiation and gravity in the effective action, may provide an analytic solution with a bounce. The solution requires a positive beta function for the total radiation term and the contraction of the Universe at the initial instant. The numerical estimate shows that the bounce may occur in an acceptable range of energies, but only under strong assumptions about the particle physics beyond the Standard Model. △ Less

Submitted 3 October, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: 15 pages, 3 figures. Several detailed explanations, discussions and references added to improve presentation, and change in the title. Fits the published version

MSC Class: 81T20; 83C47; 83F05

Journal ref: Phys. Rev. D 110 (2024) 043540

Scale-dependent cosmology from effective quantum gravity in the invariant framework

Authors: Nicolas R. Bertini, Davi C. Rodrigues, Ilya L. Shapiro

Abstract: We explore the possibility of a consistent cosmology based on the gauge-fixing independent running of the gravitational and cosmological constants ($G$ and $Λ$) in the framework of effective quantum gravity. In particular, their running in this framework was found to satisfy $G \propto Λ^4$. In the cosmological setting, the covariance of the theory provides energy conservation relations, which are… ▽ More We explore the possibility of a consistent cosmology based on the gauge-fixing independent running of the gravitational and cosmological constants ($G$ and $Λ$) in the framework of effective quantum gravity. In particular, their running in this framework was found to satisfy $G \propto Λ^4$. In the cosmological setting, the covariance of the theory provides energy conservation relations, which are impossible to satisfy with the unique scale parameter. However, by introducing the second sub-dominant scale corresponding to the higher-loop corrections and higher-derivative terms, one can close the system of equations for the running of parameters and arrive at the consistent cosmological solutions. This approach yields a change in the cosmological expansion history that affects the ratio of the Hubble parameter today to the Hubble parameter at high redshift. △ Less

Submitted 17 April, 2024; v1 submitted 21 January, 2024; originally announced January 2024.

Comments: 18 pages, 3 figures. v2: Added Cosmic Chronometers analysis. Version accepted in PDU

Large-scale Training of Foundation Models for Wearable Biosignals

Authors: Salar Abbaspourazad, Oussama Elachqar, Andrew C. Miller, Saba Emrani, Udhyakumar Nallasamy, Ian Shapiro

Abstract: Tracking biosignals is crucial for monitoring wellness and preempting the development of severe medical conditions. Today, wearable devices can conveniently record various biosignals, creating the opportunity to monitor health status without disruption to one's daily routine. Despite widespread use of wearable devices and existing digital biomarkers, the absence of curated data with annotated medi… ▽ More Tracking biosignals is crucial for monitoring wellness and preempting the development of severe medical conditions. Today, wearable devices can conveniently record various biosignals, creating the opportunity to monitor health status without disruption to one's daily routine. Despite widespread use of wearable devices and existing digital biomarkers, the absence of curated data with annotated medical labels hinders the development of new biomarkers to measure common health conditions. In fact, medical datasets are usually small in comparison to other domains, which is an obstacle for developing neural network models for biosignals. To address this challenge, we have employed self-supervised learning using the unlabeled sensor data collected under informed consent from the large longitudinal Apple Heart and Movement Study (AHMS) to train foundation models for two common biosignals: photoplethysmography (PPG) and electrocardiogram (ECG) recorded on Apple Watch. We curated PPG and ECG datasets from AHMS that include data from ~141K participants spanning ~3 years. Our self-supervised learning framework includes participant level positive pair selection, stochastic augmentation module and a regularized contrastive loss optimized with momentum training, and generalizes well to both PPG and ECG modalities. We show that the pre-trained foundation models readily encode information regarding participants' demographics and health conditions. To the best of our knowledge, this is the first study that builds foundation models using large-scale PPG and ECG data collected via wearable consumer devices $\unicode{x2013}$ prior works have commonly used smaller-size datasets collected in clinical and experimental settings. We believe PPG and ECG foundation models can enhance future wearable devices by reducing the reliance on labeled data and hold the potential to help the users improve their health. △ Less

Submitted 6 March, 2024; v1 submitted 8 December, 2023; originally announced December 2023.

Comments: Camera ready version for ICLR 2024

Kantoswki-Sachs model with a running cosmological constant and radiation

Authors: Vinicius G. Oliveira, Gil de Oliveira-Neto, Ilya L. Shapiro

Abstract: The simplest anisotropic model of the early Universe is the one with two conformal factors, which can be identified as Kantowski-Sachs metric, or the reduced version of the Bianchi-I metric. To fit the existing observational data, it is important that the anisotropy is washed out in the early stage of the evolution. We explore the possible effect of the running cosmological constant (RCC) on the d… ▽ More The simplest anisotropic model of the early Universe is the one with two conformal factors, which can be identified as Kantowski-Sachs metric, or the reduced version of the Bianchi-I metric. To fit the existing observational data, it is important that the anisotropy is washed out in the early stage of the evolution. We explore the possible effect of the running cosmological constant (RCC) on the dynamics of isotropy, in the case of the space filled by radiation. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: The paper has 16 pages and 16 figures

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)

Antisymmetric Tensor Field and Cheshire Cat Smile of the Local Conformal Symmetry

Authors: Ilya L. Shapiro

Abstract: The conformal version of the antisymmetric second-order tensor field in four spacetime dimensions does not have gauge invariance extensively discussed in the literature for more than half a century. Our first observation is that, when coupled to fermions, only the conformal version provides renormalizability of the theory at the one-loop level. General considerations are supported by the derivatio… ▽ More The conformal version of the antisymmetric second-order tensor field in four spacetime dimensions does not have gauge invariance extensively discussed in the literature for more than half a century. Our first observation is that, when coupled to fermions, only the conformal version provides renormalizability of the theory at the one-loop level. General considerations are supported by the derivation of one-loop divergences in the fermionic sector, indicating good chances for asymptotic freedom. The arguments concerning one-loop renormalizability remain valid in the presence of self-interactions and the masses for both fermion and antisymmetric tensor fields. In the flat spacetime limit, regardless of the conformal symmetry has gone, there is an expectation to meet renormalizability in all loop orders. △ Less

Submitted 20 October, 2023; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: In v2 discussions were extended and added new references, especially to the previous works about conformal operators acting on antisymmetric tensor fields. Version to be submitted for publication. 20 pages, no figures

MSC Class: 81T10; 81T15; 81T17; 81T20

Modeling Stellar CaII H & K Emission Variations: Spot Contribution to the S-index

Authors: K. Sowmya, A. I. Shapiro, L. H. M. Rouppe van der Voort, N. A. Krivova, S. K. Solanki

Abstract: The S-index is a measure of emission in the CaII H & K lines and is a widely used proxy of stellar magnetic activity. It has been assumed until now that the S-index is mainly affected by bright plage regions in the chromosphere. In particular, the effect of starspots on the S-index has been neglected. In this study we revisit this assumption. For this we analyze high-resolution observations of sun… ▽ More The S-index is a measure of emission in the CaII H & K lines and is a widely used proxy of stellar magnetic activity. It has been assumed until now that the S-index is mainly affected by bright plage regions in the chromosphere. In particular, the effect of starspots on the S-index has been neglected. In this study we revisit this assumption. For this we analyze high-resolution observations of sunspots recorded in the CaII H spectral line at the Swedish 1-m Solar Telescope and determine the contrast of spots with respect to the quiet surroundings. We find that the CaII H line core averaged over whole sunspots (including superpenumbrae) is brighter than in the quiet surroundings and that the spot contrast in the line core is comparable to the facular contrast. This allows us to get a first estimate of the influence of spots on the S-index. We show that spots increase the S-index. While this increase is quite small for the Sun, it becomes significantly larger for more active stars. Further, we show that the inclusion of the contribution of spots to the S-index strongly affects the relationship between the S-index and stellar disk area coverages by spots and faculae, and present the new relations. △ Less

Submitted 23 September, 2023; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: 11 pages, 5 figures, accepted for publication in ApJ Letters

On the local term in the anomaly-induced action of Weyl quantum gravity

Authors: Andrei O. Barvinsky, Guilherme H. S. Camargo, Alexei E. Kalugin, Nobuyoshi Ohta, Ilya L. Shapiro

Abstract: The finite local conformally non-invariant $R^2$-term emerges in the one-loop effective action of the model of quantum gravity based on the Weyl-squared classical action. This term is related to the $\Box R$ contribution to the conformal anomaly, which in a wide class of regularization schemes is determined by the second Schwinger-DeWitt (or Gilkey-Seeley) coefficient of the heat kernel expansion… ▽ More The finite local conformally non-invariant $R^2$-term emerges in the one-loop effective action of the model of quantum gravity based on the Weyl-squared classical action. This term is related to the $\Box R$ contribution to the conformal anomaly, which in a wide class of regularization schemes is determined by the second Schwinger-DeWitt (or Gilkey-Seeley) coefficient of the heat kernel expansion for inverse propagators of the theory. The calculation of this term requires evaluating the contributions of the fourth-order derivative minimal and of the second-order nonminimal operators in the tensor and vector sectors of the theory, corresponding to metric, ghost and gauge-fixing operators. To ensure the correctness of existing formulas, we derived (and confirmed) the result using a special technique of calculations, based on the heat-kernel representation of the Euclidean Green's function and the method of universal functional traces. △ Less

Submitted 28 September, 2023; v1 submitted 9 August, 2023; originally announced August 2023.

Comments: 30 pages

New rotation period measurements of 67,163 Kepler stars

Authors: Timo Reinhold, Alexander I. Shapiro, Sami K. Solanki, Gibor Basri

Abstract: The Kepler space telescope leaves a legacy of tens of thousands of stellar rotation period measurements. While many of these stars show strong periodicity, there exists an even bigger fraction of stars with irregular variability for which rotation periods are unknown. As a consequence, many stellar activity studies might be strongly biased toward the behavior of more active stars with measured rot… ▽ More The Kepler space telescope leaves a legacy of tens of thousands of stellar rotation period measurements. While many of these stars show strong periodicity, there exists an even bigger fraction of stars with irregular variability for which rotation periods are unknown. As a consequence, many stellar activity studies might be strongly biased toward the behavior of more active stars with measured rotation periods. To at least partially lift this bias, we apply a new method based on the Gradient of the Power Spectrum (GPS). The maximum of the gradient corresponds to the position of the inflection point (IP). It was shown previously that the stellar rotation period $P_{rot}$ is linked to the inflection point period $P_{IP}$ by the simple equation $P_{rot} = P_{IP}/α$, where $α$ is a calibration factor. The GPS method is superior to classical methods (such as auto-correlation functions (ACF)) because it does not require a repeatable variability pattern in the time series. From the initial sample of 142,168 stars with effective temperature $T_{eff}\leq6500K$ and surface gravity $log g\geq4.0$ in the Kepler archive, we could measure rotation periods for 67,163 stars by combining the GPS and the ACF method. We further report the first determination of a rotation period for 20,397 stars. The GPS periods show good agreement with previous period measurements using classical methods, where these are available. Furthermore, we show that the scaling factor $α$ increases for very cool stars with effective temperatures below 4000K, which we interpret as spots located at higher latitudes. We conclude that new techniques (such as the GPS method) must be applied to detect rotation periods of stars with small and more irregular variabilities. Ignoring these stars will distort the overall picture of stellar activity and, in particular, solar-stellar comparison studies. △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: 15 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 678, A24 (2023)

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

Hopf-cyclic coefficients in the braided setting

Authors: Ilya Shapiro

Abstract: Considering the monoidal category $\mathcal{C}$ obtained as modules over a Hopf algebra $H$ in a rigid braided category $\mathcal{B}$, we prove decomposition results for the Hochschild and cyclic homology categories $HH(\mathcal{C})$ and $HC(\mathcal{C})$ of $\mathcal{C}$. This is accomplished by defining a notion of a (stable) anti-Yetter-Drinfeld module with coefficients in a (stable) braided mo… ▽ More Considering the monoidal category $\mathcal{C}$ obtained as modules over a Hopf algebra $H$ in a rigid braided category $\mathcal{B}$, we prove decomposition results for the Hochschild and cyclic homology categories $HH(\mathcal{C})$ and $HC(\mathcal{C})$ of $\mathcal{C}$. This is accomplished by defining a notion of a (stable) anti-Yetter-Drinfeld module with coefficients in a (stable) braided module over $\mathcal{B}$. When the stable braided module is $HH(\mathcal{B})$, we recover $HH(\mathcal{C})$ and $HC(\mathcal{C})$. The decomposition of $HC(\mathcal{C})$ now follows from that of $HH(\mathcal{B})$. △ Less

Submitted 31 May, 2023; v1 submitted 12 May, 2023; originally announced May 2023.

Comments: Typo in Remark 5.7 fixed, added Lemma 4.26 with a definition and remark, and Corollary 4.44 with a definition, 44 pages

MSC Class: 16E35; 16E40; 16T05; 18D10; 18G90; 18M15; 18N25; 19D55

Forward modelling of brightness variations in Sun-like stars -- II. Light curves and variability

Authors: N. -E. Nèmec, A. I. Shapiro, E. Işik, S. K. Solanki, T. Reinhold

Abstract: The amplitude and morphology of light curves of solar-like stars change substantially with increasing rotation rate: brightness variations get amplified and become more regular, which has so far not been explained. We develop a modelling approach for calculating brightness variations of stars with various rotation rates and use it to explain observed trends in stellar photometric variability. We c… ▽ More The amplitude and morphology of light curves of solar-like stars change substantially with increasing rotation rate: brightness variations get amplified and become more regular, which has so far not been explained. We develop a modelling approach for calculating brightness variations of stars with various rotation rates and use it to explain observed trends in stellar photometric variability. We combine numerical simulations of magnetic Flux Emergence And Transport (FEAT) with a model for stellar brightness variability to calculate synthetic light curves of stars as observed by the Kepler telescope. We compute the distribution of magnetic flux on the stellar surface for various rotation rates and degrees of active-region nesting (i.e., the tendency of active regions to emerge in the vicinity of recently emerged ones). Using the resulting maps of the magnetic flux, we compute the rotational variability of our simulated stellar light curves as a function of rotation rate and nesting of magnetic features and compare our calculations to Kepler observations. We show that both rotation rate and degree of nesting have a strong impact on the amplitude and morphology of stellar light curves. In order to explain the variability of the bulk of \K{} targets with known rotation rates, we need to increase the degree of nesting to values much larger than on the Sun. The suggested increase of nesting with the rotation rate can provide clues to the flux emergence process for high levels of stellar activity. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 10 pages, 15 figures

Journal ref: A&A 672, A138 (2023)

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

Effective approach to the Antoniadis-Mottola model: quantum decoupling of the higher derivative terms

Authors: Wagno Cesar e Silva, Ilya L. Shapiro

Abstract: We explore the decoupling of massive ghost mode in the $4D$ (four-dimensional) theory of the conformal factor of the metric. The model was introduced by Antoniadis and Mottola in [1] and can be regarded as a close analog of the fourth-derivative quantum gravity. The analysis of the derived one-loop nonlocal form factors includes their asymptotic behavior in the UV and IR limits. In the UV (high en… ▽ More We explore the decoupling of massive ghost mode in the $4D$ (four-dimensional) theory of the conformal factor of the metric. The model was introduced by Antoniadis and Mottola in [1] and can be regarded as a close analog of the fourth-derivative quantum gravity. The analysis of the derived one-loop nonlocal form factors includes their asymptotic behavior in the UV and IR limits. In the UV (high energy) domain, our results reproduce the Minimal Subtraction scheme-based beta functions of [1]. In the IR (i.e., at low energies), the diagrams with massive ghost internal lines collapse into tadpole-type graphs without nonlocal contributions and become irrelevant. On the other hand, those structures that contribute to the running of parameters of the action and survive in the IR, are well-correlated with the divergent part (or the leading in UV contributions to the form factors), coming from the effective low-energy theory of the conformal factor. This effective theory describes only the light propagating mode. Finally, we discuss whether these results may shed light on the possible running of the cosmological constant at low energies. △ Less

Submitted 12 July, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 35 pages, 20 figures. Several detailed explanations, discussions, references and Appendices added to improve presentation. Fits the published version

Journal ref: J. High Energ. Phys. 2023 (2023) 97

Pauli equation and charged spin-1/2 particle in a weak gravitational field

Authors: Samuel W. P. Oliveira, Guilherme Y. Oyadomari, Ilya L. Shapiro

Abstract: Using the nonrelativistic approximation in the curved-space Dirac equation, the analog of the Pauli equation is derived for a weak gravitational field with a gauge fixing condition related to the synchronous gauge, in the presence of an electromagnetic field. Different from the previous works which were employing either the exact or conventional Foldy-Wouthuysen transformations, here we perform ca… ▽ More Using the nonrelativistic approximation in the curved-space Dirac equation, the analog of the Pauli equation is derived for a weak gravitational field with a gauge fixing condition related to the synchronous gauge, in the presence of an electromagnetic field. Different from the previous works which were employing either the exact or conventional Foldy-Wouthuysen transformations, here we perform calculations by directly performing nonrelativistic approximation which reduced in the power series expansion in the inverse mass of the spinning particle. On top of that, the equations of motion for the massive spin-$1/2$ charged particle are obtained. The two particular cases of the previously explored backgrounds, namely a) plane gravitational wave and b) homogeneous static gravitational field are considered for control. In the case a) we meet correspondence with the previous results. On the other hand, in case b), there is no correspondence with neither perturbative nor with exact Foldy-Wouthuysen transformations, which we also recalculate and agree with the previous works. The disagreement is a kind of a theoretical challenge and most likely occurs because the potential energy, in the particular case of Newtonian approximation, is proportional to the mass of the particle. △ Less

Submitted 13 October, 2024; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 24 pages, no figures

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)

The background information about perturbative quantum gravity

Authors: Ilya L. Shapiro

Abstract: The purpose of this Chapter is to give a general introduction and status review on the perturbative approach to quantum gravity (QG). This text is a modified version of the corresponding chapters of Part II of the recent textbook on quantum field theory (QFT) and QG, co-authored with I.L. Buchbinder and published in Oxford University Press. We discuss the choice of the starting action in the QG mo… ▽ More The purpose of this Chapter is to give a general introduction and status review on the perturbative approach to quantum gravity (QG). This text is a modified version of the corresponding chapters of Part II of the recent textbook on quantum field theory (QFT) and QG, co-authored with I.L. Buchbinder and published in Oxford University Press. We discuss the choice of the starting action in the QG models, degrees of freedom and propagator of metric perturbations, power counting and renormalizability of these models, the problems related to higher derivative theories and ghosts, such as quantum unitarity and the stability of classical solutions in general relativity; and the perspective to overcome these problems. The gauge fixing and parametrization dependencies are discussed in detail using the corresponding general QFT theorems developed in gauge theories. On top of that, we present a basic example of deriving the one-loop divergences and discuss an important example of the renormalization group in QG. The gauge invariant renormalizability of QG is considered in another Chapter of the Handbook, written together with P.M. Lavrov. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 65 pages, no figures. The chapter for the Section "Effective Quantum Gravity" of the "Handbook of Quantum Gravity" (Eds. C. Bambi, L. Modesto and I.L. Shapiro, Springer Singapore, expected in 2023). Based on sections 18-21 of the textbook "Introduction to Quantum Field Theory with Applications to Quantum Gravity" by I.L. Buchbinder and I.L. Shapiro, with permission from Oxford University Press

Report number: ISBN: 978-981-19-3079-7 MSC Class: 83C45; 81T15

Journal ref: Chapter for the Section "Perturbative Quantum Gravity'' of the "Handbook of Quantum Gravity''. Eds. C. Bambi, L. Modesto and I.L. Shapiro, Springer Singapore, 2023 - link.springer.com/referencework/10.1007/978-981-19-3079-9

Gauge invariant renormalizability of quantum gravity

Authors: P. M. Lavrov, I. L. Shapiro

Abstract: The current understanding of renormalization in quantum gravity (QG) is based on the fact that UV divergences of effective actions in the covariant QG models are covariant local expressions. This fundamental statement plays a central role in QG and, therefore, it is important to prove it for the widest possible range of the QG theories. Using the Batalin-Vilkovisky technique and the background fie… ▽ More The current understanding of renormalization in quantum gravity (QG) is based on the fact that UV divergences of effective actions in the covariant QG models are covariant local expressions. This fundamental statement plays a central role in QG and, therefore, it is important to prove it for the widest possible range of the QG theories. Using the Batalin-Vilkovisky technique and the background field method, we elaborate the proof of gauge invariant renormalizability for a generic model of quantum gravity that is diffeomorphism invariant and does not have additional, potentially anomalous, symmetries. △ Less

Submitted 20 October, 2022; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: References added. 39 pages, no figures. Invited chapter for the Section "Effective Quantum Gravity" of the "Handbook of Quantum Gravity" (Eds. C. Bambi, L. Modesto and I.L. Shapiro, Springer Singapore, expected in 2023). Review paper is mainly based on arXiv:1902.04687, with permission from Physical Review D

Report number: ISBN: 978-981-19-3079-9 MSC Class: 83C45; 81T13; 81T15; 81T70

Journal ref: Handbook of Quantum Gravity''. Eds. C. Bambi, L. Modesto and I.L. Shapiro, Springer Singapore, 2023

A comparison of stochastic and deterministic downscaling in eddy resolving ocean modelling: the Lakshadweep Sea case study

Authors: Georgy I. Shapiro, Jose M. Gonzalez-Ondina, Mohammed Salim, Jiada Tu

Abstract: This study compares the skills of two numerical models having the same horizontal resolution but based on different principles in representing meso- and submesoscale features of ocean dynamics in the Lakshadweep Sea (North Indian Ocean). The first model, titled LD20-NEMO, is based on solving primitive equations using the NEMO (Nucleus for European Modelling of the Ocean) modelling engine. The seco… ▽ More This study compares the skills of two numerical models having the same horizontal resolution but based on different principles in representing meso- and submesoscale features of ocean dynamics in the Lakshadweep Sea (North Indian Ocean). The first model, titled LD20-NEMO, is based on solving primitive equations using the NEMO (Nucleus for European Modelling of the Ocean) modelling engine. The second one, titled LD20-SDD, uses a newer Stochastic-Deterministic Downscaling method. Both models have 1/20o resolution and use the outputs from a Global Ocean Physics Analysis and Forecast model at 1/12o resolution available from Copernicus Marine Service (CMEMS). The LD20-NEMO uses only a 2D set of data from CMEMS as lateral boundary conditions. The LD20-SDD consumes the full 3D set of data from CMEMS and exploits the stochastic properties of these data to generate the downscaled field variables at higher resolution than the parent model. The skills of the three models, CMEMS, LD20-NEMO and LD20-SDD are assessed against remotely sensed and in-situ observations for the four-year period 2015-2018. All models show similar skills in reproducing temperature and salinity, however the SDD version performs slightly better than the NEMO version. This difference in resolution is particularly significant in simulation of vorticity and computation of the share of the sea occupied by highly non-linear processes. While the NEMO and SDD model show similar skill, the SDD model is more computationally efficient than the NEMO model by a large margin. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: 23 pages , 9 figures

Superflares on solar-like stars: A new method for identifying the true flare sources in photometric surveys

Authors: Valeriy Vasilyev, Timo Reinhold, Alexander I. Shapiro, Natalie A. Krivova, Ilya Usoskin, Benjamin T. Montet, Sami K. Solanki, Laurent Gizon

Abstract: Over the past years, thousands of stellar flares have been detected by harvesting data from large photometric surveys. These detections, however, do not account for potential sources of contamination such as background stars appearing in the same aperture as the primary target. We present a new method for identifying the true flare sources in large photometric surveys using data from the Kepler mi… ▽ More Over the past years, thousands of stellar flares have been detected by harvesting data from large photometric surveys. These detections, however, do not account for potential sources of contamination such as background stars appearing in the same aperture as the primary target. We present a new method for identifying the true flare sources in large photometric surveys using data from the Kepler mission. Potential flares are identified in two steps: first, we search the light curves for at least two subsequent data points exceeding a 5σ threshold above the running mean. For these two cadences, we subtract the "quiet" stellar flux from the Kepler pixel data to obtain new images where the potential flare is the main light source. In the second step, we use a Bayesian approach to fit the point spread function of the instrument to determine the most likely location of the flux excess on the detector. We applied our method to 5862 main-sequence stars with near-solar effective temperatures. We found 2274 events exceeding the 5-sigma in at least two consecutive points in the light curves. Applying the second step reduced this number to 342 superflares. Of these, 283 flares happened on 178 target stars, 47 events are associated with fainter background stars, and in 10 cases, the flare location cannot be distinguished between the target and a background star. We also present cases where flares have been reported previously but our technique could not attribute them to the target star. We conclude that 1) identifying outliers in the light curves alone is insufficient to attribute them to stellar flares and 2) flares can only be uniquely attributed to a certain star when the instrument pixel-level data together with the point spread function are taken into account. As a consequence, previous flare statistics are likely contaminated by instrumental effects and unresolved astrophysical sources. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 668, A167 (2022)

Measuring periods in aperiodic light curves -- Applying the GPS method to infer rotation periods of solar-like stars

Authors: Timo Reinhold, Alexander I. Shapiro, Sami K. Solanki, Gibor Basri

Abstract: Light curves of solar-like stars are known to show highly irregular variability. As a consequence, standard frequency analysis methods often fail to detect the correct rotation period. Recently, Shapiro et al. (2020) showed that the periods of such stars could still be measured by considering the Gradient of the Power Spectrum (GPS) instead of the power spectrum itself. In this study, the GPS meth… ▽ More Light curves of solar-like stars are known to show highly irregular variability. As a consequence, standard frequency analysis methods often fail to detect the correct rotation period. Recently, Shapiro et al. (2020) showed that the periods of such stars could still be measured by considering the Gradient of the Power Spectrum (GPS) instead of the power spectrum itself. In this study, the GPS method is applied to model light curves of solar-like stars covering all possible inclination angles and a large range of metallicities and observational noise levels. The model parameters are chosen such that they resemble those of many stars in the Kepler field. We show that the GPS method is able to detect the correct rotation period in 40% of all considered cases, which is more than ten times higher than the detection rate of standard techniques. Thus, we conclude that the GPS method is ideally suited to measure periods of those Kepler stars lacking such a measurement so far. We also show that the GPS method is significantly superior to auto-correlation methods when starspot lifetimes are shorter than a few rotation periods. GPS begins to yield rotation periods that are too short when dominant spot lifetimes are shorter than one rotation period. We conclude that new methods are generally needed to reliably detect rotation periods from sufficiently aperiodic time series -- these periods will otherwise remain undetected. △ Less

Submitted 26 September, 2022; originally announced September 2022.

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

Hopf Cyclic Cohomology and Beyond

Authors: Masoud Khalkhali, Ilya Shapiro

Abstract: This paper is an introduction to Hopf cyclic cohomology with an emphasis on its most recent developments. We cover three major areas: the original definition of Hopf cyclic cohomology by Connes and Moscovici as an outgrowth of their study of transverse index theory on foliated manifolds, the introduction of Hopf cyclic cohomology with coefficients by Hajac-Khalkhali-Rangipour-Sommerhauser, and fin… ▽ More This paper is an introduction to Hopf cyclic cohomology with an emphasis on its most recent developments. We cover three major areas: the original definition of Hopf cyclic cohomology by Connes and Moscovici as an outgrowth of their study of transverse index theory on foliated manifolds, the introduction of Hopf cyclic cohomology with coefficients by Hajac-Khalkhali-Rangipour-Sommerhauser, and finally the latest episode on unifying the coefficients as well as extending the notion to more general settings beyond Hopf algebras. In particular, the last section discusses the relative Hopf cyclic theory that arises in the braided monoidal category settings. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 27 pages, Proceedings of Symposia in Pure Mathematics

MSC Class: 58B34; 19D55; 16E40; 16T05; 16T10; 18D05; 18D10; 18D15

Faculae cancel out on the surfaces of active Suns

Authors: N. -E. Nèmec, A. I. Shapiro, E. Işık, K. Sowmya, S. K. Solanki, N. A. Krivova, R. H. Cameron, L. Gizon

Abstract: Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they… ▽ More Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 Figures, accepted for publication in ApJL

Predictions of Astrometric Jitter for Sun-like Stars. III. Fast Rotators

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

Abstract: A breakthrough in exoplanet detections is foreseen with the unprecedented astrometric measurement capabilities offered by instrumentation aboard Gaia space observatory. Besides, astrometric discoveries of exoplanets are expected from the planned space mission, Small-JASMINE. In this setting, the present series of papers focuses on estimating the effect of magnetic activity of G2V-type host stars o… ▽ More A breakthrough in exoplanet detections is foreseen with the unprecedented astrometric measurement capabilities offered by instrumentation aboard Gaia space observatory. Besides, astrometric discoveries of exoplanets are expected from the planned space mission, Small-JASMINE. In this setting, the present series of papers focuses on estimating the effect of magnetic activity of G2V-type host stars on the astrometric signal. This effect interferes with the astrometric detections of Earth-mass planets. While the first two papers considered stars rotating at the solar rotation rate, this paper focuses on stars having solar effective temperature and metallicity but rotating faster than the Sun, and consequently more active. By simulating the distribution of active regions on such stars using the Flux Emergence And Transport model, we show that the contribution of magnetic activity to the astrometric measurements becomes increasingly significant with increasing rotation rates. We further show that the jitter for the most variable periodic Kepler stars is high enough to be detected by Gaia. Furthermore, due to a decrease in the facula-to-spot area ratio for more active stars, the magnetic jitter is found to be spot-dominated for rapid rotators. Our simulations of the astrometric jitter has the potential to aid the interpretation of data from Gaia and upcoming space astrometry missions. △ Less

Submitted 28 June, 2022; v1 submitted 15 June, 2022; originally announced June 2022.

Comments: 21 pages, 20 figures, accepted for publication in ApJ

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)

Anomaly-induced vacuum effective action with torsion: covariant solution and ambiguities

Authors: Guilherme H. S. Camargo, Ilya L. Shapiro

Abstract: The conformal anomaly in curved spacetime with antisymmetric torsion is reconsidered, taking into account new important details. We formulate, for the first time, the covariant solution of the anomaly-induced effective action. The covariant effective action includes local terms corresponding to total derivatives in the conformal anomaly. The contribution of massless fermions to these terms is char… ▽ More The conformal anomaly in curved spacetime with antisymmetric torsion is reconsidered, taking into account new important details. We formulate, for the first time, the covariant solution of the anomaly-induced effective action. The covariant effective action includes local terms corresponding to total derivatives in the conformal anomaly. The contribution of massless fermions to these terms is characterized by multiplicative anomaly, coming from two different choices of doubling for the spinor operator. On the other hand, the nonlocal part of anomaly-induced action is free of ambiguities and admits a low-energy limit, similar to the effective potential in the metric-scalar theory. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: 21 pages, LaTeX, 1 figure

MSC Class: 81T20; 83D05; 81T17

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)

Nonlocality of quantum matter corrections and cosmological constant running

Authors: E. V. Gorbar, Ilya L. Shapiro

Abstract: Semiclassical contributions to the gravitational action include various terms with zero, two, and four derivatives of the metric as well as nonlocal form factors for these terms. Contributions to some of these terms could be confused with others on a specific metric background or for a particular gauge fixing. We present a critical analysis of the recent works where the tensor structure and the nu… ▽ More Semiclassical contributions to the gravitational action include various terms with zero, two, and four derivatives of the metric as well as nonlocal form factors for these terms. Contributions to some of these terms could be confused with others on a specific metric background or for a particular gauge fixing. We present a critical analysis of the recent works where the tensor structure and the number of derivatives in the action of gravity were not properly taken into account. Taking these relevant aspects into account, we show that although some contributions owing to the quantum fluctuations of massive or massless scalar as well as fermion and vector fields may be attributed to the "running" of the cosmological constant, in reality they correspond to the fourth derivative terms of the action. △ Less

Submitted 9 August, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: Explanations extended, results unchanged, fits the published version. 15 pages, 1 figure

MSC Class: 81T20; 81T17

Trace anomaly and induced action for a metric-scalar background

Authors: Manuel Asorey, Wagno Cesar e Silva, Ilya L. Shapiro, Públio R. B. do Vale

Abstract: The conformal anomaly and anomaly-induced effective action represent useful and economic ways to describe semiclassical contributions to the action of gravity. We discuss the anomaly in the case when the background is formed by metric and scalar fields and formulate the induced action in two standard covariant forms. The analysis of induced action at low energies reveals existing connection to the… ▽ More The conformal anomaly and anomaly-induced effective action represent useful and economic ways to describe semiclassical contributions to the action of gravity. We discuss the anomaly in the case when the background is formed by metric and scalar fields and formulate the induced action in two standard covariant forms. The analysis of induced action at low energies reveals existing connection to the renormalization group and effective potential. The classification of anomalous terms is extended to the scalar background and ambiguities in the total derivative terms in the anomaly are considered using Pauli-Villars regularization. △ Less

Submitted 7 February, 2023; v1 submitted 31 January, 2022; originally announced February 2022.

Comments: Explanations and references added. Fits the version to be published in EPJC. 22 pages, no figures

MSC Class: 81T50; 81T20; 81T15; 83C45

Journal ref: Eur. Phys. J. C 83 (2023) 157

The Effect of Stellar Contamination on Low-resolution Transmission Spectroscopy: Needs Identified by NASA's Exoplanet Exploration Program Study Analysis Group 21

Authors: Benjamin V. Rackham, Néstor Espinoza, Svetlana V. Berdyugina, Heidi Korhonen, Ryan J. MacDonald, Benjamin T. Montet, Brett M. Morris, Mahmoudreza Oshagh, Alexander I. Shapiro, Yvonne C. Unruh, Elisa V. Quintana, Robert T. Zellem, Dániel Apai, Thomas Barclay, Joanna K. Barstow, Giovanni Bruno, Ludmila Carone, Sarah L. Casewell, Heather M. Cegla, Serena Criscuoli, Catherine Fischer, Damien Fournier, Mark S. Giampapa, Helen Giles, Aishwarya Iyer , et al. (36 additional authors not shown)

Abstract: Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectru… ▽ More Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectrum of the star being occulted. However, stars are not homogeneous, constant light sources but have temporally evolving photospheres and chromospheres with inhomogeneities like spots, faculae, plages, granules, and flares. This SAG brought together an interdisciplinary team of more than 100 scientists, with observers and theorists from the heliophysics, stellar astrophysics, planetary science, and exoplanetary atmosphere research communities, to study the current research needs that can be addressed in this context to make the most of transit studies from current NASA facilities like HST and JWST. The analysis produced 14 findings, which fall into three Science Themes encompassing (1) how the Sun is used as our best laboratory to calibrate our understanding of stellar heterogeneities ("The Sun as the Stellar Benchmark"), (2) how stars other than the Sun extend our knowledge of heterogeneities ("Surface Heterogeneities of Other Stars") and (3) how to incorporate information gathered for the Sun and other stars into transit studies ("Mapping Stellar Knowledge to Transit Studies"). In this invited review, we largely reproduce the final report of SAG21 as a contribution to the peer-reviewed literature. △ Less

Submitted 17 March, 2023; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: Invited review in press at RASTI. Based on the ExoPAG SAG21 report (arXiv:2201.09905v1) and refined via feedback from three reviewers. 75 pages, 30 figures, 5 tables

Effective quantum gravity, cosmological constant and the Standard Model of particle physics

Authors: Breno L. Giacchini, Tibério de Paula Netto, Ilya L. Shapiro

Abstract: The renormalization group in effective quantum gravity can be consistently formulated using the Vilkovisky and DeWitt version of effective action and assuming a non-zero cosmological constant. Taking into account that the vacuum counterpart of the cosmological constant is dramatically different from the observed energy density of vacuum, the running of the last quantity in the late cosmology indic… ▽ More The renormalization group in effective quantum gravity can be consistently formulated using the Vilkovisky and DeWitt version of effective action and assuming a non-zero cosmological constant. Taking into account that the vacuum counterpart of the cosmological constant is dramatically different from the observed energy density of vacuum, the running of the last quantity in the late cosmology indicates strong constraints on the physics beyond the minimal Standard Model of particle physics. △ Less

Submitted 16 June, 2022; v1 submitted 12 December, 2021; originally announced December 2021.

Comments: 15 pages, no figures. Added discussions concerning the general status of the Vilkovisky-DeWitt effective action, corresponding references, and the new Appendix, illustrating the choice of the metric. Conclusions are unchanged. Fits the version accepted in PRD

MSC Class: 81T17; 83C45

Journal ref: Phys. Rev. D 105, 126025 (2022)

The running vacuum in effective quantum gravity

Authors: Breno L. Giacchini, Tibério de Paula Netto, Ilya L. Shapiro

Abstract: We briefly review the previous works on the renormalization group in quantum general relativity with the cosmological constant, based on the Vilkovisky and DeWitt version of effective action. On top of that, we discuss the prospects of the applications of this version of renormalization group to the cosmological models with a running Newton constant and vacuum energy density. We briefly review the previous works on the renormalization group in quantum general relativity with the cosmological constant, based on the Vilkovisky and DeWitt version of effective action. On top of that, we discuss the prospects of the applications of this version of renormalization group to the cosmological models with a running Newton constant and vacuum energy density. △ Less

Submitted 2 December, 2021; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: Partially based on the talk presented by I.Sh. at the on-line workshop on quantum gravity, higher derivatives and nonlocality, organized by Luca Buoninfante and Sravan Kumar in TIT/Tokyo. v2: small misprints corrected

MSC Class: 83C45; 81T17

Journal ref: Nuovo Cim.C 45 (2022) 34

Universal leading quantum correction to the Newton potential

Authors: Tibério de Paula Netto, Leonardo Modesto, Ilya L. Shapiro

Abstract: The derivation of effective quantum gravity corrections to Newton's potential is an important step in the whole effective quantum field theory approach. We hereby add new strong arguments in favor of omitting all the diagrams with internal lines of the massive sources, and we also recalculate the corrections to the Newtonian potential using functional methods in an arbitrary parametrization of the… ▽ More The derivation of effective quantum gravity corrections to Newton's potential is an important step in the whole effective quantum field theory approach. We hereby add new strong arguments in favor of omitting all the diagrams with internal lines of the massive sources, and we also recalculate the corrections to the Newtonian potential using functional methods in an arbitrary parametrization of the quantum fluctuations of the metric. The general proof of the gauge- and parametrization-independence within this approach is also explicitly given. On top of that, we argue that the universality of the result holds regardless of the details of the ultraviolet completion of quantum gravity theory. Indeed, it turns out that the logarithm quantum correction depends only on the low energy spectrum of the theory that is responsible for the analytic properties of loop's amplitudes. △ Less

Submitted 16 February, 2022; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: 25 pages, slightly extended discussions. Fits the version to appear in EPJC

On the vector conformal models in an arbitrary dimension

Authors: Manuel Asorey, Lesław Rachwał, Ilya L. Shapiro, Wagno Cesar e Silva

Abstract: The conventional model of the gauge vector field is invariant under the local conformal symmetry only in the four-dimensional space ($4d$). Conformal generalization to an arbitrary dimension $d$ is impossible even for the free theory, differently from scalar and fermion fields. We discuss how to overcome this restriction and eventually construct four vector conformal actions. One of these models i… ▽ More The conventional model of the gauge vector field is invariant under the local conformal symmetry only in the four-dimensional space ($4d$). Conformal generalization to an arbitrary dimension $d$ is impossible even for the free theory, differently from scalar and fermion fields. We discuss how to overcome this restriction and eventually construct four vector conformal actions. One of these models is the particular case of the previously known conformal theory of $n$-forms and others are new, up to our knowledge. In some of these models the gauge invariance is preserved, two of the new models are described by local actions with auxiliary compensating scalar fields, and the extended version of one of these models is on shell equivalent to the last, non-analytic, purely metric version. △ Less

Submitted 3 October, 2021; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 15 pages. A few formulations made more precise, fits the version accepted in EPJP

Journal ref: Eur. Phys. J. Plus 136 (2021) 1043

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)

Renormalization Group in Six-derivative Quantum Gravity

Authors: Leslaw Rachwal, Leonardo Modesto, Aleksandr Pinzul, Ilya L. Shapiro

Abstract: The exact one-loop beta functions for the four-derivative terms (Weyl tensor squared, Ricci scalar squared and the Gauss-Bonnet) are derived for the minimal six-derivative quantum gravity (QG) theory in four spacetime dimensions. The calculation is performed by means of the Barvinsky and Vilkovisky generalized Schwinger-DeWitt technique. With this result we gain, for the first time, the full set o… ▽ More The exact one-loop beta functions for the four-derivative terms (Weyl tensor squared, Ricci scalar squared and the Gauss-Bonnet) are derived for the minimal six-derivative quantum gravity (QG) theory in four spacetime dimensions. The calculation is performed by means of the Barvinsky and Vilkovisky generalized Schwinger-DeWitt technique. With this result we gain, for the first time, the full set of the relevant beta functions in a super-renormalizable model of QG. The complete set of renormalization group (RG) equations, including also those for the Newton and the cosmological constant, is solved explicitly in the general case and for the six-derivative Lee-Wick (LW) quantum gravity proposed in a previous paper by two of the authors. In the ultraviolet regime, the minimal theory is shown to be asymptotically free and describes free gravitons in Minkowski or (anti-) de Sitter ((A)dS) backgrounds, depending on the initial conditions for the RG equations. The ghostlike states appear in complex conjugate pairs at any energy scale consistently with the LW prescription. However, owing to the running, these ghosts may become tachyons. We argue that an extension of the theory that involves operators cubic in Riemann tensor may change the beta functions and hence be capable of overcoming this problem. △ Less

Submitted 27 October, 2021; v1 submitted 28 April, 2021; originally announced April 2021.

Comments: 27 pages, discussion and some comments added, citations added, spelling mistakes corrected; coincides with the journal version

Journal ref: Phys. Rev. D 104 (2021) 085018

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

The Planetary Ephemeris Program: Capability, Comparison, and Open Source Availability

Authors: John F. Chandler, James B. R. Battat, Thomas W. Murphy, Jr., Daniel Reardon, Robert D. Reasenberg, Irwin I. Shapiro

Abstract: We describe for the first time in the scientific literature the Planetary Ephemeris Program (PEP), an open-source general-purpose astrometric data analysis program. We discuss, in particular, the implementation of pulsar timing analysis, which was recently upgraded in PEP to handle more options. This implementation was done independently of other pulsar programs, with minor exceptions that we disc… ▽ More We describe for the first time in the scientific literature the Planetary Ephemeris Program (PEP), an open-source general-purpose astrometric data analysis program. We discuss, in particular, the implementation of pulsar timing analysis, which was recently upgraded in PEP to handle more options. This implementation was done independently of other pulsar programs, with minor exceptions that we discuss. We illustrate the implementation of this capability by comparing the post-fit residuals from the analyses of time-of-arrival observations by both PEP and Tempo2. The comparison shows substantial agreement: 22 ns rms differences for 1,065 pulse time-of-arrival measurements for the millisecond pulsar in a binary system, PSR J1909-3744 (pulse period 2.947108 ms; full-width half-maximum of pulse 43 $μ$s) for epochs in the interval from December 2002 to February 2011. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 9 pages, 1 figure, submitted to Astronomical Journal (AJ)

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

Predictions of astrometric jitter for Sun-like stars. I. The model and its application to the Sun as seen from the ecliptic

Authors: Alexander I. Shapiro, Sami K. Solanki, Natalie A. Krivova

Abstract: The advent of Gaia, capable of measuring stellar wobbles caused by orbiting planets, raised an interest to the astrometric detection of exoplanets. Another source of such wobbles (often also called jitter) is stellar magnetic activity. A quantitative assessment of the stellar astrometric jitter is important for a more reliable astrometric detection and characterisation of exoplanets. We calculate… ▽ More The advent of Gaia, capable of measuring stellar wobbles caused by orbiting planets, raised an interest to the astrometric detection of exoplanets. Another source of such wobbles (often also called jitter) is stellar magnetic activity. A quantitative assessment of the stellar astrometric jitter is important for a more reliable astrometric detection and characterisation of exoplanets. We calculate the displacement of the solar photocentre due to the magnetic activity for an almost 16-year period (February 2, 1999 - August 1, 2014). We also investigate how the displacement depends on the spectral passband chosen for observations, including the wavelength range to be covered by the upcoming Small-JASMINE mission of JAXA. This is done by extending the SATIRE-S model for solar irradiance variability to calculating the displacement of the solar photocentre caused by the magnetic features on the surface of the Sun. We found that the peak to peak amplitude of the solar photocentre displacement would reach 0.5 mas if the Sun were located 10 pc away from the observer and observed in the Gaia G filter. This is by far too small to be detected by the Gaia mission. However, the Sun is a relatively inactive star so that one can expect significantly larger signals for younger, and, consequently, more active stars. The model developed in this study can be combined with the simulations of emergence and surface transport of magnetic flux which have recently became available to model the astrometric jitter over the broad range of magnetic activities. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: Accepted for publication in ApJ

Bounce and stability in the early cosmology with anomaly-induced corrections

Authors: Wagno Cesar e Silva, Ilya L. Shapiro

Abstract: An extremely fast exponential expansion of the Universe is typical for the stable version of the inflationary model, based on the anomaly-induced action of gravity. The total amount of exponential $e$-folds could be very large, before the transition to the unstable version and the beginning of the Starobinsky inflation. Thus, the stable exponential expansion can be seen as a pre-inflationary semic… ▽ More An extremely fast exponential expansion of the Universe is typical for the stable version of the inflationary model, based on the anomaly-induced action of gravity. The total amount of exponential $e$-folds could be very large, before the transition to the unstable version and the beginning of the Starobinsky inflation. Thus, the stable exponential expansion can be seen as a pre-inflationary semiclassical cosmological solution. We explore whether this stable phase could follow after the bounce, subsequent to the contraction of the Universe. Extending the previous consideration of the bounce, we explore both stable expansion and the bounce solutions in the models with non-zero cosmological constant and the presence of background radiation. The critical part of the analysis concerns stability for small perturbations of the Hubble parameter. It is shown that the stability is possible for the variations in the bounce region, but not in the sufficiently distant past in the contraction phase. △ Less

Submitted 25 December, 2020; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 22 pages, 12 figures. Extended version, added more details concerning phase diagrams, many new references and a few comments. Article written by invitation for the Special Issue "Physics and Mathematics of the Dark Universe" in Symmetry (Editor Sergei Ketov)

Journal ref: Symmetry 13 (2021) 50