Showing 1–50 of 189 results for author: Shapiro, A

Analyzing Near-Field Intensity Distribution in Subwavelength Gratings through Cylindrical Wave Decomposition

Authors: A. S. Bereza, A. E. Chernyavsky, S. V. Perminov, D. A. Shapiro

Abstract: The investigation into the scattering of plane waves by a periodic array of parallel cylinders utilizes the method of cylindrical wave decomposition, thereby reducing the problem complexity to a series of linear algebraic equations. This methodology proves particularly efficacious when the diameter of cylinders is significantly less than the wavelength of incident wave, resulting in a rapid diminu… ▽ More The investigation into the scattering of plane waves by a periodic array of parallel cylinders utilizes the method of cylindrical wave decomposition, thereby reducing the problem complexity to a series of linear algebraic equations. This methodology proves particularly efficacious when the diameter of cylinders is significantly less than the wavelength of incident wave, resulting in a rapid diminution of the solution coefficients as a function of azimuth numbers. Such a reductionist approach facilitates the computation of scattered radiation intensity in near field. Subsequent cross-validation with numerical results corroborates the theoretical findings, showcasing a qualitative concordance between the two. This study underscores the efficacy of cylindrical wave decomposition in simplifying and accurately modeling wave scattering phenomena in structured media. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 14 pages, 8 figures, bibliography 31 item

Distributionally robust stochastic optimal control

Authors: Alexander Shapiro, Yan Li

Abstract: The main goal of this paper is to discuss the construction of distributionally robust counterparts of stochastic optimal control problems. Randomized and non-randomized policies are considered. In particular, necessary and sufficient conditions for the existence of non-randomized policies are given. The main goal of this paper is to discuss the construction of distributionally robust counterparts of stochastic optimal control problems. Randomized and non-randomized policies are considered. In particular, necessary and sufficient conditions for the existence of non-randomized policies are given. △ Less

Submitted 9 June, 2024; originally announced June 2024.

The PLATO Mission

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

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

Submitted 8 June, 2024; originally announced June 2024.

Stellar Surface Magnetic Fields Impact Limb Darkening

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

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

Submitted 29 February, 2024; originally announced March 2024.

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

Cluster structure on genus 2 spherical DAHA: seven-colored flower

Authors: Semeon Arthamonov, Leonid Chekhov, Philippe Di Francesco, Rinat Kedem, Gus Schrader, Alexander Shapiro, Michael Shapiro

Abstract: We construct an embedding of the Arthamonov-Shakirov algebra of genus 2 knot operators into the quantized coordinate ring of the cluster Poisson variety of exceptional finite mutation type $X_7$. The embedding is equivariant with respect to the action of the mapping class group of the closed surface of genus 2. The cluster realization of the mapping class group action leads to a formula for the co… ▽ More We construct an embedding of the Arthamonov-Shakirov algebra of genus 2 knot operators into the quantized coordinate ring of the cluster Poisson variety of exceptional finite mutation type $X_7$. The embedding is equivariant with respect to the action of the mapping class group of the closed surface of genus 2. The cluster realization of the mapping class group action leads to a formula for the coefficient of each monomial in the genus 2 Macdonald polynomial of type $A_1$ as sum over lattice points in a convex polyhedron in 7-dimensional space. △ Less

Submitted 15 June, 2024; v1 submitted 25 February, 2024; originally announced February 2024.

Ruijsenaars wavefunctions as modular group matrix coefficients

Authors: Philippe Di Francesco, Rinat Kedem, Sergey Khoroshkin, Gus Schrader, Alexander Shapiro

Abstract: We give a description of the Hallnäs--Ruijsenaars eigenfunctions of the 2-particle hyperbolic Ruijsenaars system as matrix coefficients for the order 4 element $S\in SL(2,\mathbb{Z})$ acting on the Hilbert space of $GL(2)$ quantum Teichmüller theory on the punctured torus. The $GL(2)$ Macdonald polynomials are then obtained as special values of the analytic continuation of these matrix coefficient… ▽ More We give a description of the Hallnäs--Ruijsenaars eigenfunctions of the 2-particle hyperbolic Ruijsenaars system as matrix coefficients for the order 4 element $S\in SL(2,\mathbb{Z})$ acting on the Hilbert space of $GL(2)$ quantum Teichmüller theory on the punctured torus. The $GL(2)$ Macdonald polynomials are then obtained as special values of the analytic continuation of these matrix coefficients. The main tool used in the proof is the cluster structure on the moduli space of framed $GL(2)$-local systems on the punctured torus, and an $SL(2,\mathbb{Z})$-equivariant embedding of the $GL(2)$ spherical DAHA into the quantized coordinate ring of the corresponding cluster Poisson variety. △ Less

Submitted 15 June, 2024; v1 submitted 21 February, 2024; originally announced February 2024.

Visualizing Magnetic Order in Self-Assembly of Superparamagnetic Nanoparticles

Authors: Xingyuan Lu, Ji Zou, Minh Pham, Arjun Rana, Chen-Ting Liao, Emma Cating Subramanian, Xuefei Wu, Yuan Hung Lo, Charles S. Bevis, Robert M. Karl Jr, Serban Lepadatu, Young-Sang Yu, Yaroslav Tserkovnyak, Thomas P. Russell, David A. Shapiro, Henry C. Kapteyn, Margaret M. Murnane, Robert Streubel, Jianwei Miao

Abstract: We use soft x-ray vector-ptychographic tomography to determine the three-dimensional magnetization field in superparamagnetic nanoparticles self-assembled at the liquid-liquid interface and reveal the magnetic order induced by layered structure. The spins in individual nanoparticles become more aligned with increasing number of layers, resulting in a larger net magnetization. Our experimental resu… ▽ More We use soft x-ray vector-ptychographic tomography to determine the three-dimensional magnetization field in superparamagnetic nanoparticles self-assembled at the liquid-liquid interface and reveal the magnetic order induced by layered structure. The spins in individual nanoparticles become more aligned with increasing number of layers, resulting in a larger net magnetization. Our experimental results show a magnetic short-range order in the monolayer due to the proliferation of thermally induced magnetic vortices and a magnetic long-range order in the bilayer and trilayer, stemming from the strengthened dipolar interactions that effectively suppress thermal fluctuations. We also observe a screening effect of magnetic vortices and the attractive interaction between the magnetic vortices with opposite topological charges. Our work demonstrates the crucial role of layered structure in shaping the magnetization of nanoparticle assemblies, providing new opportunities to modulate these properties through strategic layer engineering. △ Less

Submitted 2 January, 2024; originally announced January 2024.

BGLS: A Python Package for the Gate-by-Gate Sampling Algorithm to Simulate Quantum Circuits

Authors: Alex Shapiro, Ryan LaRose

Abstract: The classical simulation of quantum computers is in general a computationally hard problem. To emulate the behavior of realistic devices, it is sufficient to sample bitstrings from circuits. Recently, arXiv:2112.08499 introduced the so-called gate-by-gate sampling algorithm to sample bitstrings and showed it to be computationally favorable in many cases. Here we present bgls, a Python package whic… ▽ More The classical simulation of quantum computers is in general a computationally hard problem. To emulate the behavior of realistic devices, it is sufficient to sample bitstrings from circuits. Recently, arXiv:2112.08499 introduced the so-called gate-by-gate sampling algorithm to sample bitstrings and showed it to be computationally favorable in many cases. Here we present bgls, a Python package which implements this sampling algorithm. bgls has native support for several states and is highly flexible for use with additional states. We show how to install and use bgls, discuss optimizations in the algorithm, and demonstrate its utility on several problems. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 7 pages, 9 figures, included in Workshops of The International Conference on High Performance Computing, Network, Storage, and Analysis (SC-W 2023)

Journal ref: SC-W '23: Proceedings of the SC '23 Workshops of The International Conference on High Performance Computing, Network, Storage, and Analysis, November 2023, Pages 1436-1442

A roadmap for the atmospheric characterization of terrestrial exoplanets with JWST

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

Abstract: Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a bet… ▽ More Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets. △ Less

Submitted 22 July, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Journal ref: Nature Astronomy (2024) 8, 810-818

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)

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

Rectangularity and duality of distributionally robust Markov Decision Processes

Authors: Yan Li, Alexander Shapiro

Abstract: The main goal of this paper is to discuss several approaches to formulation of distributionally robust counterparts of Markov Decision Processes, where the transition kernels are not specified exactly but rather are assumed to be elements of the corresponding ambiguity sets. The intent is to clarify some connections between the game and static formulations of distributionally robust MDPs, and deli… ▽ More The main goal of this paper is to discuss several approaches to formulation of distributionally robust counterparts of Markov Decision Processes, where the transition kernels are not specified exactly but rather are assumed to be elements of the corresponding ambiguity sets. The intent is to clarify some connections between the game and static formulations of distributionally robust MDPs, and delineate the role of rectangularity associated with ambiguity sets in determining these connections. △ Less

Submitted 5 May, 2024; v1 submitted 21 August, 2023; originally announced August 2023.

Episodic Bayesian Optimal Control with Unknown Randomness Distributions

Authors: Alexander Shapiro, Enlu Zhou, Yifan Lin, Yuhao Wang

Abstract: Stochastic optimal control with unknown randomness distributions has been studied for a long time, encompassing robust control, distributionally robust control, and adaptive control. We propose a new episodic Bayesian approach that incorporates Bayesian learning with optimal control. In each episode, the approach learns the randomness distribution with a Bayesian posterior and subsequently solves… ▽ More Stochastic optimal control with unknown randomness distributions has been studied for a long time, encompassing robust control, distributionally robust control, and adaptive control. We propose a new episodic Bayesian approach that incorporates Bayesian learning with optimal control. In each episode, the approach learns the randomness distribution with a Bayesian posterior and subsequently solves the corresponding Bayesian average estimate of the true problem. The resulting policy is exercised during the episode, while additional data/observations of the randomness are collected to update the Bayesian posterior for the next episode. We show that the resulting episodic value functions and policies converge almost surely to their optimal counterparts of the true problem if the parametrized model of the randomness distribution is correctly specified. We further show that the asymptotic convergence rate of the episodic value functions is of the order $O(N^{-1/2})$. We develop an efficient computational method based on stochastic dual dynamic programming for a class of problems that have convex value functions. Our numerical results on a classical inventory control problem verify the theoretical convergence results and demonstrate the effectiveness of the proposed computational method. △ Less

Submitted 16 August, 2023; originally announced August 2023.

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

Numerical Methods for Convex Multistage Stochastic Optimization

Authors: Guanghui Lan, Alexander Shapiro

Abstract: Optimization problems involving sequential decisions in a stochastic environment were studied in Stochastic Programming (SP), Stochastic Optimal Control (SOC) and Markov Decision Processes (MDP). In this paper we mainly concentrate on SP and SOC modelling approaches. In these frameworks there are natural situations when the considered problems are convex. Classical approach to sequential optimizat… ▽ More Optimization problems involving sequential decisions in a stochastic environment were studied in Stochastic Programming (SP), Stochastic Optimal Control (SOC) and Markov Decision Processes (MDP). In this paper we mainly concentrate on SP and SOC modelling approaches. In these frameworks there are natural situations when the considered problems are convex. Classical approach to sequential optimization is based on dynamic programming. It has the problem of the so-called ``Curse of Dimensionality", in that its computational complexity increases exponentially with increase of dimension of state variables. Recent progress in solving convex multistage stochastic problems is based on cutting planes approximations of the cost-to-go (value) functions of dynamic programming equations. Cutting planes type algorithms in dynamical settings is one of the main topics of this paper. We also discuss Stochastic Approximation type methods applied to multistage stochastic optimization problems. From the computational complexity point of view, these two types of methods seem to be complimentary to each other. Cutting plane type methods can handle multistage problems with a large number of stages, but a relatively smaller number of state (decision) variables. On the other hand, stochastic approximation type methods can only deal with a small number of stages, but a large number of decision variables. △ Less

Submitted 27 March, 2023; originally announced March 2023.

MSC Class: 65K05; 90C15; 90C39; 90C40

Statistical Limit Theorems in Distributionally Robust Optimization

Authors: Jose Blanchet, Alexander Shapiro

Abstract: The goal of this paper is to develop methodology for the systematic analysis of asymptotic statistical properties of data driven DRO formulations based on their corresponding non-DRO counterparts. We illustrate our approach in various settings, including both phi-divergence and Wasserstein uncertainty sets. Different types of asymptotic behaviors are obtained depending on the rate at which the unc… ▽ More The goal of this paper is to develop methodology for the systematic analysis of asymptotic statistical properties of data driven DRO formulations based on their corresponding non-DRO counterparts. We illustrate our approach in various settings, including both phi-divergence and Wasserstein uncertainty sets. Different types of asymptotic behaviors are obtained depending on the rate at which the uncertainty radius decreases to zero as a function of the sample size and the geometry of the uncertainty sets. △ Less

Submitted 26 March, 2023; originally announced March 2023.

MSC Class: 90C15

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

Uncovering Multiple Intrinsic Chiral Phases in (PEA)2PbI4 halide Perovskites

Authors: Shahar Zuri, Arthur Shapiro, Leeor Kronik, Efrat Lifshitz

Abstract: Two-dimensional (2D) halide perovskites offer a unique platform to investigate the ground-state of materials possessing significant anharmonicity. In contrast to three-dimensional perovskites, their 2D counterparts offer substantially fewer degrees of freedom, resulting in multiple well-defined crystal structures. In this work, we thoroughly investigate the anharmonic ground-state of the benchmark… ▽ More Two-dimensional (2D) halide perovskites offer a unique platform to investigate the ground-state of materials possessing significant anharmonicity. In contrast to three-dimensional perovskites, their 2D counterparts offer substantially fewer degrees of freedom, resulting in multiple well-defined crystal structures. In this work, we thoroughly investigate the anharmonic ground-state of the benchmark (PEA)2PbI4 compound, using complementary information from low-temperature x-ray diffraction (XRD) and photoluminescence spectroscopy, supported by density functional theory (DFT) calculations. We extrapolate four crystallographic configurations from the low-temperature XRD. These configurations imply that the ground-state has an intrinsic disorder stemming from two coexisting chiral sub-lattices, each with a bi-oriented organic spacer molecule. We further show evidence that these chiral structures form unevenly populated ground-states, portraying uneven anharmonicity, where the state population may be tuned by surface effects. Our results uncover a disordered ground-state that may induce intrinsic grain boundaries, which cannot be ignored in practical applications. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: manuscripts: 16 pages, 3 figures, 1 table. SI: 14 pages, 5 figures, 9 tables

On Localization of Tight Closure in Line-$S_4$ Quartics

Authors: Levi Borevitz, Naima Nader, Theodore J. Sandstrom, Amelia Shapiro, Austyn Simpson, Jenna Zomback

Abstract: Building on work of Brenner and Monsky from 2010 and on a Hilbert-Kunz calculation of Monsky from 1998, we exhibit a novel example of a hypersurface over $\overline{\mathbb{F}_2}$ in which tight closure does not commute with localization. Our methods involve a surprising tiling argument using Sierpiński triangles, as well as an inspection of a certain dynamical system in characteristic two. Building on work of Brenner and Monsky from 2010 and on a Hilbert-Kunz calculation of Monsky from 1998, we exhibit a novel example of a hypersurface over $\overline{\mathbb{F}_2}$ in which tight closure does not commute with localization. Our methods involve a surprising tiling argument using Sierpiński triangles, as well as an inspection of a certain dynamical system in characteristic two. △ Less

Submitted 6 November, 2022; originally announced November 2022.

Comments: Written during the 2022 SMALL REU at Williams College. Comments welcome

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)

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

Study of detecting behavioral signatures within DeepFake videos

Authors: Qiaomu Miao, Sinhwa Kang, Stacy Marsella, Steve DiPaola, Chao Wang, Ari Shapiro

Abstract: There is strong interest in the generation of synthetic video imagery of people talking for various purposes, including entertainment, communication, training, and advertisement. With the development of deep fake generation models, synthetic video imagery will soon be visually indistinguishable to the naked eye from a naturally capture video. In addition, many methods are continuing to improve to… ▽ More There is strong interest in the generation of synthetic video imagery of people talking for various purposes, including entertainment, communication, training, and advertisement. With the development of deep fake generation models, synthetic video imagery will soon be visually indistinguishable to the naked eye from a naturally capture video. In addition, many methods are continuing to improve to avoid more careful, forensic visual analysis. Some deep fake videos are produced through the use of facial puppetry, which directly controls the head and face of the synthetic image through the movements of the actor, allow the actor to 'puppet' the image of another. In this paper, we address the question of whether one person's movements can be distinguished from the original speaker by controlling the visual appearance of the speaker but transferring the behavior signals from another source. We conduct a study by comparing synthetic imagery that: 1) originates from a different person speaking a different utterance, 2) originates from the same person speaking a different utterance, and 3) originates from a different person speaking the same utterance. Our study shows that synthetic videos in all three cases are seen as less real and less engaging than the original source video. Our results indicate that there could be a behavioral signature that is detectable from a person's movements that is separate from their visual appearance, and that this behavioral signature could be used to distinguish a deep fake from a properly captured video. △ Less

Submitted 8 August, 2024; v1 submitted 6 August, 2022; originally announced August 2022.

Comments: 9 pages

AlexU-AIC at Arabic Hate Speech 2022: Contrast to Classify

Authors: Ahmad Shapiro, Ayman Khalafallah, Marwan Torki

Abstract: Online presence on social media platforms such as Facebook and Twitter has become a daily habit for internet users. Despite the vast amount of services the platforms offer for their users, users suffer from cyber-bullying, which further leads to mental abuse and may escalate to cause physical harm to individuals or targeted groups. In this paper, we present our submission to the Arabic Hate Speech… ▽ More Online presence on social media platforms such as Facebook and Twitter has become a daily habit for internet users. Despite the vast amount of services the platforms offer for their users, users suffer from cyber-bullying, which further leads to mental abuse and may escalate to cause physical harm to individuals or targeted groups. In this paper, we present our submission to the Arabic Hate Speech 2022 Shared Task Workshop (OSACT5 2022) using the associated Arabic Twitter dataset. The shared task consists of 3 sub-tasks, sub-task A focuses on detecting whether the tweet is offensive or not. Then, For offensive Tweets, sub-task B focuses on detecting whether the tweet is hate speech or not. Finally, For hate speech Tweets, sub-task C focuses on detecting the fine-grained type of hate speech among six different classes. Transformer models proved their efficiency in classification tasks, but with the problem of over-fitting when fine-tuned on a small or an imbalanced dataset. We overcome this limitation by investigating multiple training paradigms such as Contrastive learning and Multi-task learning along with Classification fine-tuning and an ensemble of our top 5 performers. Our proposed solution achieved 0.841, 0.817, and 0.476 macro F1-average in sub-tasks A, B, and C respectively. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Journal ref: Proceedings of the OSACT 2022 Workshop, LREC2022, June 2022, 200-208

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)

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

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

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

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

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

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

Chromospheric Extension of the MURaM Code

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

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

Submitted 6 April, 2022; originally announced April 2022.

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

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

Leaders or Followers? A Temporal Analysis of Tweets from IRA Trolls

Authors: Siva K. Balasubramanian, Mustafa Bilgic, Aron Culotta, Libby Hemphill, Anita Nikolich, Matthew A. Shapiro

Abstract: The Internet Research Agency (IRA) influences online political conversations in the United States, exacerbating existing partisan divides and sowing discord. In this paper we investigate the IRA's communication strategies by analyzing trending terms on Twitter to identify cases in which the IRA leads or follows other users. Our analysis focuses on over 38M tweets posted between 2016 and 2017 from… ▽ More The Internet Research Agency (IRA) influences online political conversations in the United States, exacerbating existing partisan divides and sowing discord. In this paper we investigate the IRA's communication strategies by analyzing trending terms on Twitter to identify cases in which the IRA leads or follows other users. Our analysis focuses on over 38M tweets posted between 2016 and 2017 from IRA users (n=3,613), journalists (n=976), members of Congress (n=526), and politically engaged users from the general public (n=71,128). We find that the IRA tends to lead on topics related to the 2016 election, race, and entertainment, suggesting that these are areas both of strategic importance as well having the highest potential impact. Furthermore, we identify topics where the IRA has been relatively ineffective, such as tweets on military, political scandals, and violent attacks. Despite many tweets on these topics, the IRA rarely leads the conversation and thus has little opportunity to influence it. We offer our proposed methodology as a way to track the strategic choices of future influence operations in real-time. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: ICWSM 2022

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

Risk-Averse Stochastic Optimal Control: an efficiently computable statistical upper bound

Authors: Vincent Guigues, Alexander Shapiro, Yi Cheng

Abstract: In this paper, we discuss an application of the SDDP type algorithm to nested risk-averse formulations of Stochastic Optimal Control (SOC) problems. We propose a construction of a statistical upper bound for the optimal value of risk-averse SOC problems. This outlines an approach to a solution of a long standing problem in that area of research. The bound holds for a large class of convex and mono… ▽ More In this paper, we discuss an application of the SDDP type algorithm to nested risk-averse formulations of Stochastic Optimal Control (SOC) problems. We propose a construction of a statistical upper bound for the optimal value of risk-averse SOC problems. This outlines an approach to a solution of a long standing problem in that area of research. The bound holds for a large class of convex and monotone conditional risk mappings. Finally, we show the validity of the statistical upper bound to solve a real-life stochastic hydro-thermal planning problem. △ Less

Submitted 2 May, 2023; v1 submitted 17 December, 2021; originally announced December 2021.

Bayesian Distributionally Robust Optimization

Authors: Alexander Shapiro, Enlu Zhou, Yifan Lin

Abstract: We introduce a new framework, Bayesian Distributionally Robust Optimization (Bayesian-DRO), for data-driven stochastic optimization where the underlying distribution is unknown. Bayesian-DRO contrasts with most of the existing DRO approaches in the use of Bayesian estimation of the unknown distribution. To make computation of Bayesian updating tractable, Bayesian-DRO first assumes the underlying d… ▽ More We introduce a new framework, Bayesian Distributionally Robust Optimization (Bayesian-DRO), for data-driven stochastic optimization where the underlying distribution is unknown. Bayesian-DRO contrasts with most of the existing DRO approaches in the use of Bayesian estimation of the unknown distribution. To make computation of Bayesian updating tractable, Bayesian-DRO first assumes the underlying distribution takes a parametric form with unknown parameter and then computes the posterior distribution of the parameter. To address the model uncertainty brought by the assumed parametric distribution, Bayesian-DRO constructs an ambiguity set of distributions with the assumed parametric distribution as the reference distribution and then optimizes with respect to the worst case in the ambiguity set. We show the consistency of the Bayesian posterior distribution and subsequently the convergence of objective functions and optimal solutions of Bayesian-DRO. Our consistency result of the Bayesian posterior requires simpler assumptions than the classical literature on Bayesian consistency. We also consider several approaches for selecting the ambiguity set size in Bayesian-DRO and compare them numerically. Our numerical experiments demonstrate the out-of-sample performance of Bayesian-DRO in comparison with Kullback-Leibler-based (KL-) and Wasserstein-based empirical DRO as well as risk-neutral Bayesian Risk Optimization. Our numerical results shed light on how to choose the modeling framework (Bayesian-DRO, KL-DRO, Wasserstein-DRO) for specific problems, but the choice for general problems still remains an important and open question. △ Less

Submitted 8 February, 2023; v1 submitted 15 December, 2021; originally announced December 2021.

Multiwavelength Mitigation of Stellar Activity in Astrometric Planet Detection

Authors: Avi Kaplan-Lipkin, Bruce Macintosh, Alexander Madurowicz, Krishnamurthy Sowmya, Alexander Shapiro, Natalie Krivova, Sami K. Solanki

Abstract: Astrometry has long been a promising technique for exoplanet detection. At the theoretical limits, astrometry would allow for the detection of smaller planets than previously seen by current exoplanet search methods, but stellar activity may make these theoretical limits unreachable. Astrometric jitter of a Sun-like star due to magnetic activity in its photosphere induces apparent variability in t… ▽ More Astrometry has long been a promising technique for exoplanet detection. At the theoretical limits, astrometry would allow for the detection of smaller planets than previously seen by current exoplanet search methods, but stellar activity may make these theoretical limits unreachable. Astrometric jitter of a Sun-like star due to magnetic activity in its photosphere induces apparent variability in the photocenter of order $0.5\ \textrm{m}R_\odot$. This jitter creates a fundamental astrophysical noise floor preventing detection of lower-mass planets in a single spectral band. By injecting planet orbits into simulated solar data at five different passbands, we investigate mitigation of this fundamental astrometric noise using correlations across passbands. For a true solar analog and a planet at 1 au semimajor axis, the $6σ$ detection limit set by stellar activity for an ideal telescope at the best single passband is $0.01$ Earth masses. We found that pairs of passbands with highly correlated astrometric jitter due to stellar activity, but with less motion in the redder band, enable higher-precision measurements of the common signal from the planet. Using this method improves detectable planet masses at 1 au by up to a factor of $10$, corresponding to at best $0.005$ Earth masses for a Sun-like star with a perfect telescope. Given these results, we recommend that future astrometry missions consider proceeding with two or more passbands to reduce noise due to stellar activity. △ Less

Submitted 19 April, 2022; v1 submitted 12 December, 2021; originally announced December 2021.

Comments: 11 pages, 9 figures. Published in AJ

Journal ref: 2022 AJ 163 205

High-field magnetic structure of the triangular antiferromagnet RbFe(MoO4)2

Authors: Yu. A. Sakhratov, O. Prokhnenko, A. Ya. Shapiro, H. D. Zhou, L. E. Svistov, A. P. Reyes, O. A. Petrenko

Abstract: The magnetic H - T phase diagram of a quasi-two-dimensional antiferromagnet RbFe(MoO4)2 with an equilateral triangular lattice structure is studied with 87Rb NMR and neutron diffraction techniques. This combination of experimental techniques allows us to determine the ordered components of the magnetic moments on the Fe3+ ions within various high-field phases - the Y, UUD, V, and fan structures, s… ▽ More The magnetic H - T phase diagram of a quasi-two-dimensional antiferromagnet RbFe(MoO4)2 with an equilateral triangular lattice structure is studied with 87Rb NMR and neutron diffraction techniques. This combination of experimental techniques allows us to determine the ordered components of the magnetic moments on the Fe3+ ions within various high-field phases - the Y, UUD, V, and fan structures, stabilized in the compound by the in-plane magnetic field. It is also established that the transition from the V to the fan phase is of first-order, whereas the transition from the fan phase to the polarized paramagnetic phase is continuous. An analysis of the NMR spectra shows that the high-field fan phase of RbFe(MoO4)2 can be successfully described by a periodic commensurate oscillation of the magnetic moments around the field direction in each Fe layer combined with an incommensurate modulation of the magnetic structure perpendicular to the layers. △ Less

Submitted 12 December, 2021; originally announced December 2021.

Comments: 13 pages, 15 figures

Correlative image learning of chemo-mechanics in phase-transforming solids

Authors: Haitao D. Deng, Hongbo Zhao, Norman L. Jin, Lauren Hughes, Benjamin Savitzky, Colin Ophus, Dimitrios Fraggedakis, András Borbély, Young-Sang Yu, Eder Lomeli, Rui Yan, Jueyi Liu, David A. Shapiro, Wei Cai, Martin Z. Bazant, Andrew M. Minor, William C. Chueh

Abstract: Constitutive laws underlie most physical processes in nature. However, learning such equations in heterogeneous solids (e.g., due to phase separation) is challenging. One such relationship is between composition and eigenstrain, which governs the chemo-mechanical expansion in solids. In this work, we developed a generalizable, physically-constrained image-learning framework to algorithmically lear… ▽ More Constitutive laws underlie most physical processes in nature. However, learning such equations in heterogeneous solids (e.g., due to phase separation) is challenging. One such relationship is between composition and eigenstrain, which governs the chemo-mechanical expansion in solids. In this work, we developed a generalizable, physically-constrained image-learning framework to algorithmically learn the chemo-mechanical constitutive law at the nanoscale from correlative four-dimensional scanning transmission electron microscopy and X-ray spectro-ptychography images. We demonstrated this approach on Li$_X$FePO$_4$, a technologically-relevant battery positive electrode material. We uncovered the functional form of composition-eigenstrain relation in this two-phase binary solid across the entire composition range (0 $\leq$ X $\leq$ 1), including inside the thermodynamically-unstable miscibility gap. The learned relation directly validates Vegard's law of linear response at the nanoscale. Our physics-constrained data-driven approach directly visualizes the residual strain field (by removing the compositional and coherency strain), which is otherwise impossible to quantify. Heterogeneities in the residual strain arise from misfit dislocations and were independently verified by X-ray diffraction line profile analysis. Our work provides the means to simultaneously quantify chemical expansion, coherency strain and dislocations in battery electrodes, which has implications on rate capabilities and lifetime. Broadly, this work also highlights the potential of integrating correlative microscopy and image learning for extracting material properties and physics. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Correlative analysis of structure and chemistry of LixFePO4 platelets using 4D-STEM and X-ray ptychography

Authors: L. A. Hughes, Benjamin H. Savitzky, Haitao D. Deng, Norman L. Jin, Eder G. Lomeli, Young-Sang Yu, David A. Shapiro, Patrick Herring, Abraham Anapolsky, William C. Chueh, Colin Ophus, Andrew M. Minor

Abstract: Lithium iron phosphate (LixFePO4), a cathode material used in rechargeable Li-ion batteries, phase separates upon de/lithiation under equilibrium. The interfacial structure and chemistry within these cathode materials affects Li-ion transport, and therefore battery performance. Correlative imaging of LixFePO4 was performed using four-dimensional scanning transmission electron microscopy (4D-STEM),… ▽ More Lithium iron phosphate (LixFePO4), a cathode material used in rechargeable Li-ion batteries, phase separates upon de/lithiation under equilibrium. The interfacial structure and chemistry within these cathode materials affects Li-ion transport, and therefore battery performance. Correlative imaging of LixFePO4 was performed using four-dimensional scanning transmission electron microscopy (4D-STEM), scanning transmission X-ray microscopy (STXM), and X-ray ptychography in order to analyze the local structure and chemistry of the same particle set. Over 50,000 diffraction patterns from 10 particles provided measurements of both structure and chemistry at a nanoscale spatial resolution (16.6-49.5 nm) over wide (several micron) fields-of-view with statistical robustness.LixFePO4 particles at varying stages of delithiation were measured to examine the evolution of structure and chemistry as a function of delithiation. In lithiated and delithiated particles, local variations were observed in the degree of lithiation even while local lattice structures remained comparatively constant, and calculation of linear coefficients of chemical expansion suggest pinning of the lattice structures in these populations. Partially delithiated particles displayed broadly core-shell-like structures, however, with highly variable behavior both locally and per individual particle that exhibited distinctive intermediate regions at the interface between phases, and pockets within the lithiated core that correspond to FePO4 in structure and chemistry.The results provide insight into the LixFePO4 system, subtleties in the scope and applicability of Vegards law (linear lattice parameter-composition behavior) under local versus global measurements, and demonstrate a powerful new combination of experimental and analytical modalities for bridging the crucial gap between local and statistical characterization. △ Less

Submitted 9 July, 2021; originally announced July 2021.

Comments: 17 pages, 4 figures

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

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

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

Submitted 3 July, 2021; originally announced July 2021.

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

Journal ref: 2021, ApJ, 919, 94

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

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

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

Submitted 28 May, 2021; originally announced May 2021.

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

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

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

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

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

Submitted 28 April, 2021; originally announced April 2021.

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

Direct observation of 3D topological spin textures and their interactions using soft x-ray vector ptychography

Authors: Arjun Rana, Chen-Ting Liao, Ezio Iacocca, Ji Zou, Minh Pham, Emma-Elizabeth Cating Subramanian, Yuan Hung Lo, Sinéad A. Ryan, Xingyuan Lu, Charles S. Bevis, Robert M. Karl Jr, Andrew J. Glaid, Young-Sang Yu, Pratibha Mahale, David A. Shapiro, Sadegh Yazdi, Thomas E. Mallouk, Stanley J. Osher, Henry C. Kapteyn, Vincent H. Crespi, John V. Badding, Yaroslav Tserkovnyak, Margaret M. Murnane, Jianwei Miao

Abstract: Magnetic topological defects are energetically stable spin configurations characterized by symmetry breaking. Vortices and skyrmions are two well-known examples of 2D spin textures that have been actively studied for both fundamental interest and practical applications. However, experimental evidence of the 3D spin textures has been largely indirect or qualitative to date, due to the difficulty of… ▽ More Magnetic topological defects are energetically stable spin configurations characterized by symmetry breaking. Vortices and skyrmions are two well-known examples of 2D spin textures that have been actively studied for both fundamental interest and practical applications. However, experimental evidence of the 3D spin textures has been largely indirect or qualitative to date, due to the difficulty of quantitively characterizing them within nanoscale volumes. Here, we develop soft x-ray vector ptychography to quantitatively image the 3D magnetization vector field in a frustrated superlattice with 10 nm spatial resolution. By applying homotopy theory to the experimental data, we quantify the topological charge of hedgehogs and anti-hedgehogs as emergent magnetic monopoles and probe their interactions inside the frustrated superlattice. We also directly observe virtual hedgehogs and anti-hedgehogs created by magnetically inert voids. We expect that this new quantitative imaging method will open the door to study 3D topological spin textures in a broad class of magnetic materials. Our work also demonstrates that magnetically frustrated superlattices could be used as a new platform to investigate hedgehog interactions and dynamics and to exploit optimized geometries for information storage and transport applications. △ Less

Submitted 26 April, 2021; originally announced April 2021.

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

Volume 78, 1 October 2014, Pages 56-64

Authors: A. Lupinacci, J. Kacher, A. A. Shapiro, P. Hosemann, A. M. Minor

Abstract: Characterizing plasticity mechanisms below the ductile-to-brittle transition temperature is traditionally difficult to accomplish in asystematic fashion. Here, we use a new experimental setup to perform in situ cryogenic mechanical testing of pure Sn micropillars at room temperature and at 142°C. Subsequent electron microscopy characterization of the micropillars shows a clear difference in the de… ▽ More Characterizing plasticity mechanisms below the ductile-to-brittle transition temperature is traditionally difficult to accomplish in asystematic fashion. Here, we use a new experimental setup to perform in situ cryogenic mechanical testing of pure Sn micropillars at room temperature and at 142°C. Subsequent electron microscopy characterization of the micropillars shows a clear difference in the deformation mechanisms at room temperature and at cryogenic temperatures. At room temperature, the Sn micropillars deformed through dislocation plasticity, while at142°C they exhibited both higher strength and deformation twinning. Two different orientations were tested, a symmetric (100) orientation and a non-symmetric (451) orientation. The deformation mechanisms were found to be the same for both orientations △ Less

Submitted 1 April, 2021; originally announced April 2021.

Journal ref: Acta Materialia Volume 78, 1 October 2014, Pages 56-64

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

Quantum decorated character stacks

Authors: David Jordan, Ian Le, Gus Schrader, Alexander Shapiro

Abstract: We initiate the study of decorated character stacks and their quantizations using the framework of stratified factorization homology. We thereby extend the construction by Fock and Goncharov of (quantum) decorated character varieties to encompass also the stacky points, in a way that is both compatible with cutting and gluing and equivariant with respect to canonical actions of the modular group o… ▽ More We initiate the study of decorated character stacks and their quantizations using the framework of stratified factorization homology. We thereby extend the construction by Fock and Goncharov of (quantum) decorated character varieties to encompass also the stacky points, in a way that is both compatible with cutting and gluing and equivariant with respect to canonical actions of the modular group of the surface. In the cases $G=SL_2,PGL_2$ we construct a system of categorical charts and flips on the quantum decorated character stacks which generalize the well--known cluster structures on the Fock--Goncharov moduli spaces. △ Less

Submitted 24 February, 2021; originally announced February 2021.

MSC Class: 13F60 81T45

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)

Irradiance Variations of the Sun and Sun-Like Stars -- Overview of Topical Collection

Authors: Greg Kopp, Alexander Shapiro

Abstract: This topical collection summarizes recent advances in observing and modeling irradiance variations of the Sun and Sun-like stars, emphasizing the links between surface magnetic fields and the resulting solar and stellar variability. In particular, the articles composing this collection summarize recent progress in i) solar-irradiance measurements; ii) modeling of solar- and stellar-irradiance vari… ▽ More This topical collection summarizes recent advances in observing and modeling irradiance variations of the Sun and Sun-like stars, emphasizing the links between surface magnetic fields and the resulting solar and stellar variability. In particular, the articles composing this collection summarize recent progress in i) solar-irradiance measurements; ii) modeling of solar- and stellar-irradiance variability; and iii) understanding of the effects of such variability on Earth's climate and exoplanet environments. This topical-collection overview article gives background and more details on these aspects of variability. △ Less

Submitted 13 February, 2021; originally announced February 2021.

Comments: accepted for Solar Physics