Showing 1–50 of 64 results for author: Garraffo, C

Unraveling the role of merger histories in the population of Insitu stars: linking IllustrisTNG cosmological simulation to H3 survey

Authors: Razieh Emami, Lars Hernquist, Randall Smith, James F. Steiner, Grant Tremblay, Douglas Finkbeiner, Mark Vogelsberger, Josh Grindlay, Federico Marinacci, Kung-Yi Su, Cecilia Garraffo, Yuan-Sen Ting, Phillip A. Cargile, Rebecca L. Davies, Chloë E. Benton, Yijia Li, Letizia Bugiani, Amir H. Khoram, Sownak Bose

Abstract: We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to… ▽ More We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to define insitu stars. To uncover the underlying factors contributing to deviations between TNG50 and H3 distributions, we scrutinize correlation coefficients among internal drivers(e.g., virial radius, star formation rate [SFR]) and merger-related parameters (such as the effective mass-ratio, mean distance, average redshift, total number of mergers, average spin-ratio and maximum spin alignment between merging galaxies). Notably, we identify significant correlations between deviations from observational data and key parameters such as the median slope of virial radius, mean SFR values, and the rate of SFR change across different redshift scans. Furthermore, positive correlations emerge between deviations from observational data and parameters related to galaxy mergers. We validate these correlations using the Random Forest Regression method. Our findings underscore the invaluable insights provided by the H3 survey in unraveling the cosmic history of galaxies akin to the Milky Way, thereby advancing our understanding of galactic evolution and shedding light on the formation and evolution of Milky Way-like galaxies in cosmological simulations. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 20 pages, 13 figures

Dracula's Chivito: discovery of a large edge-on protoplanetary disk with Pan-STARRS

Authors: Ciprian T. Berghea, Ammar Bayyari, Michael L. Sitko, Jeremy J. Drake, Ana Mosquera, Cecilia Garraffo, Thomas Petit, Ray W. Russell, Korash D. Assani

Abstract: We report the serendipitous discovery of a large edge-on protoplanetary disk associated with the infrared source IRAS 23077+6707. The disk's apparent size in the Pan-STARRS (PS1) images is ~11", making this one of the largest known disks on the sky. It is likely a young system, still surrounded by the envelope which is very faint but still visible in the PS1 images in the northern part (alternativ… ▽ More We report the serendipitous discovery of a large edge-on protoplanetary disk associated with the infrared source IRAS 23077+6707. The disk's apparent size in the Pan-STARRS (PS1) images is ~11", making this one of the largest known disks on the sky. It is likely a young system, still surrounded by the envelope which is very faint but still visible in the PS1 images in the northern part (alternatively this structure could be filaments from the disk itself). We use the PS1 magnitudes and other available photometric data to construct the spectral energy distribution (SED) of the disk. An optical spectrum indicates that the obscured star is hot, most likely late A. We adopt a distance of 300 pc for this object based on Gaia DR3 extinctions. We model the system using the HOCHUNK3D radiative transfer software and find that the system is consistent with a hot star of effective temperature 8000 K surrounded by a disk of size 1650 AU and mass 0.2 M_solar at inclination 82 degrees. △ Less

Submitted 23 April, 2024; v1 submitted 1 February, 2024; originally announced February 2024.

Comments: 11 pages, 4 figures, accepted in ApJL

Heating of the Atmospheres of Short-orbit Exoplanets by Their Rapid Orbital Motion Through an Extreme Space Environment

Authors: Ofer Cohen, Alex Glocer, Cecilia Garraffo, Julian Alvarado-Gomez, Jeremy Drake, Kristina Monsch, Farah Fauth Puigdomenech

Abstract: Exoplanets with short orbit period reside very close to their host stars. They transition very rapidly between different sectors of the circumstellar space environment along their orbit, leading to large variations of the magnetic field in the vicinity of the planet on short timescales. This rapid change of the magnetic flux through the conducting and resistive layer of the planetary upper atmosph… ▽ More Exoplanets with short orbit period reside very close to their host stars. They transition very rapidly between different sectors of the circumstellar space environment along their orbit, leading to large variations of the magnetic field in the vicinity of the planet on short timescales. This rapid change of the magnetic flux through the conducting and resistive layer of the planetary upper atmosphere may drive currents that dissipate in the form of Joule Heating. Here, we estimate the amount of Joule Heating dissipation in the upper atmosphere of Trappist-1e, and two hypothetical planets orbiting the Sun in close-in orbits. We find that the rapid orbital motion could drive a significant amount of atmospheric heating and could significantly affect the planetary atmosphere escape rate. Thus, the process should be accounted for when studying the long-term evolution of exoplanetary atmospheres. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 8 pages, 4 figures, accepted to ApJ

Hyperspectral shadow removal with Iterative Logistic Regression and latent Parametric Linear Combination of Gaussians

Authors: Core Francisco Park, Maya Nasr, Manuel Pérez-Carrasco, Eleanor Walker, Douglas Finkbeiner, Cecilia Garraffo

Abstract: Shadow detection and removal is a challenging problem in the analysis of hyperspectral images. Yet, this step is crucial for analyzing data for remote sensing applications like methane detection. In this work, we develop a shadow detection and removal method only based on the spectrum of each pixel and the overall distribution of spectral values. We first introduce Iterative Logistic Regression (I… ▽ More Shadow detection and removal is a challenging problem in the analysis of hyperspectral images. Yet, this step is crucial for analyzing data for remote sensing applications like methane detection. In this work, we develop a shadow detection and removal method only based on the spectrum of each pixel and the overall distribution of spectral values. We first introduce Iterative Logistic Regression (ILR) to learn a spectral basis in which shadows can be linearly classified. We then model the joint distribution of the mean radiance and the projection coefficients of the spectra onto the above basis as a parametric linear combination of Gaussians. We can then extract the maximum likelihood mixing parameter of the Gaussians to estimate the shadow coverage and to correct the shadowed spectra. Our correction scheme reduces correction artefacts at shadow borders. The shadow detection and removal method is applied to hyperspectral images from MethaneAIR, a precursor to the satellite MethaneSAT. △ Less

Submitted 23 December, 2023; originally announced December 2023.

Generating Images of the M87* Black Hole Using GANs

Authors: Arya Mohan, Pavlos Protopapas, Keerthi Kunnumkai, Cecilia Garraffo, Lindy Blackburn, Koushik Chatterjee, Sheperd S. Doeleman, Razieh Emami, Christian M. Fromm, Yosuke Mizuno, Angelo Ricarte

Abstract: In this paper, we introduce a novel data augmentation methodology based on Conditional Progressive Generative Adversarial Networks (CPGAN) to generate diverse black hole (BH) images, accounting for variations in spin and electron temperature prescriptions. These generated images are valuable resources for training deep learning algorithms to accurately estimate black hole parameters from observati… ▽ More In this paper, we introduce a novel data augmentation methodology based on Conditional Progressive Generative Adversarial Networks (CPGAN) to generate diverse black hole (BH) images, accounting for variations in spin and electron temperature prescriptions. These generated images are valuable resources for training deep learning algorithms to accurately estimate black hole parameters from observational data. Our model can generate BH images for any spin value within the range of [-1, 1], given an electron temperature distribution. To validate the effectiveness of our approach, we employ a convolutional neural network to predict the BH spin using both the GRMHD images and the images generated by our proposed model. Our results demonstrate a significant performance improvement when training is conducted with the augmented dataset while testing is performed using GRMHD simulated data, as indicated by the high R2 score. Consequently, we propose that GANs can be employed as cost effective models for black hole image generation and reliably augment training datasets for other parameterization algorithms. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 11 pages, 7 figures. Accepted by Monthly Notices of the Royal Astronomical Society Journal

Faster Bayesian inference with neural network bundles and new results for $f(R)$ models

Authors: Augusto T. Chantada, Susana J. Landau, Pavlos Protopapas, Claudia G. Scóccola, Cecilia Garraffo

Abstract: In the last few years, there has been significant progress in the development of machine learning methods tailored to astrophysics and cosmology. We have recently applied one of these, namely, the neural network bundle method, to the cosmological scenario. Moreover, we showed that in some cases the computational times of the Bayesian inference process can be reduced. In this paper, we present an i… ▽ More In the last few years, there has been significant progress in the development of machine learning methods tailored to astrophysics and cosmology. We have recently applied one of these, namely, the neural network bundle method, to the cosmological scenario. Moreover, we showed that in some cases the computational times of the Bayesian inference process can be reduced. In this paper, we present an improvement to the neural network bundle method that results in a significant reduction of the computational times of the statistical analysis. The novel aspect consists of the use of the neural network bundle method to calculate the luminosity distance of type Ia supernovae, which is usually computed through an integral with numerical methods. In this work, we have applied this improvement to the Hu-Sawicki and Starobinsky $f(R)$ models. We also performed a statistical analysis with data from type Ia supernovae of the Pantheon+ compilation and cosmic chronometers. Another original aspect of this work is the different treatment we provide for the absolute magnitude of type Ia supernovae during the inference process, which results in different estimates of the distortion parameter than the ones obtained in the literature. We show that the statistical analyses carried out with our new method require lower computational times than the ones performed with both the numerical and the neural network method from our previous work. This reduction in time is more significant in the case of a difficult computational problem such as the ones addressed in this work. △ Less

Submitted 7 June, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 14 pages, 5 figures, 3 tables

Journal ref: Phys. Rev. D 109, 123514 (2024)

Linking circumstellar disk lifetimes to the rotational evolution of low-mass stars

Authors: Kristina Monsch, Jeremy J. Drake, Cecilia Garraffo, Giovanni Picogna, Barbara Ercolano

Abstract: The high-energy radiation emitted by young stars can have a strong influence on their rotational evolution at later stages. This is because internal photoevaporation is one of the major drivers of the dispersal of circumstellar disks, which surround all newly born low-mass stars during the first few million years of their evolution. Employing an internal EUV/X-ray photoevaporation model, we have d… ▽ More The high-energy radiation emitted by young stars can have a strong influence on their rotational evolution at later stages. This is because internal photoevaporation is one of the major drivers of the dispersal of circumstellar disks, which surround all newly born low-mass stars during the first few million years of their evolution. Employing an internal EUV/X-ray photoevaporation model, we have derived a simple recipe for calculating realistic inner disk lifetimes of protoplanetary disks. This prescription was implemented into a magnetic morphology-driven rotational evolution model and is used to investigate the impact of disk-locking on the spin evolution of low-mass stars. We find that the length of the disk-locking phase has a profound impact on the subsequent rotational evolution of a young star, and the implementation of realistic disk lifetimes leads to an improved agreement of model outcomes with observed rotation period distributions for open clusters of various ages. However, for both young star-forming regions tested in our model, the strong bimodality in rotation periods that is observed in hPer could not be recovered. hPer is only successfully recovered, if the model is started from a double-peaked distribution with an initial disk fraction of $65\,\%$. However, at an age of only $\sim 1\,\mathrm{Myr}$, such a low disk fraction can only be achieved if an additional disk dispersal process, such as external photoevaporation, is invoked. These results therefore highlight the importance of including realistic disk dispersal mechanisms in rotational evolution models of young stars. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: accepted for publication in ApJ

Breakthroughs in Cool Star Physics with the Line Emission Mapper X-ray Probe

Authors: Jeremy J. Drake, Julián Alvarado Gomez, Costanza Argiroffi, Ettore Flaccomio, Cecilia Garraffo, Nicolas Grosso, Nazma Islam, Margarita Karovska, Vinay L. Kashyap, Kristina Monsch, Jan-Uwe Ness, Salvatore Sciortino, Bradford Wargelin

Abstract: We outline some of the highlights of the scientific case for the advancement of stellar high energy physics using the Line Emission Mapper X-ray Probe ({\it LEM}). The key to advancements with LEM lie in its large effective area -- up to 100 times that of the {\it Chandra} MEG -- and 1~eV spectral resolution. The large effective area opens up for the first time the ability to study time-dependent… ▽ More We outline some of the highlights of the scientific case for the advancement of stellar high energy physics using the Line Emission Mapper X-ray Probe ({\it LEM}). The key to advancements with LEM lie in its large effective area -- up to 100 times that of the {\it Chandra} MEG -- and 1~eV spectral resolution. The large effective area opens up for the first time the ability to study time-dependent phenomena on their natural timescales at high resolution, such as flares and coronal mass ejections, and also opens the sky to much fainter targets than available to {\it Chandra} or {\it XMM-Newton}. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: A Line Emission Mapper X-ray Probe White Paper

Rapid Spectral Parameter Prediction for Black Hole X-Ray Binaries using Physicalised Autoencoders

Authors: Ethan Tregidga, James Steiner, Cecilia Garraffo, Carter Rhea, Mayeul Aubin

Abstract: Black hole X-ray binaries (BHBs) offer insights into extreme gravitational environments and the testing of general relativity. The X-ray spectrum collected by NICER offers valuable information on the properties and behaviour of BHBs through spectral fitting. However, traditional spectral fitting methods are slow and scale poorly with model complexity. This paper presents a new semi-supervised auto… ▽ More Black hole X-ray binaries (BHBs) offer insights into extreme gravitational environments and the testing of general relativity. The X-ray spectrum collected by NICER offers valuable information on the properties and behaviour of BHBs through spectral fitting. However, traditional spectral fitting methods are slow and scale poorly with model complexity. This paper presents a new semi-supervised autoencoder neural network for parameter prediction and spectral reconstruction of BHBs, showing an improvement of up to a factor of 2,700 in speed while maintaining comparable accuracy. The approach maps the spectral features from the numerous outbursts catalogued by NICER and generalises them to new systems for efficient and accurate spectral fitting. The effectiveness of this approach is demonstrated in the spectral fitting of BHBs and holds promise for use in other areas of astronomy and physics for categorising large datasets. △ Less

Submitted 2 April, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 13 pages, 13 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 529, Issue 2, April 2024, Pages 1654-1666

Simulation-based Inference for Exoplanet Atmospheric Retrieval: Insights from winning the Ariel Data Challenge 2023 using Normalizing Flows

Authors: Mayeul Aubin, Carolina Cuesta-Lazaro, Ethan Tregidga, Javier Viaña, Cecilia Garraffo, Iouli E. Gordon, Mercedes López-Morales, Robert J. Hargreaves, Vladimir Yu. Makhnev, Jeremy J. Drake, Douglas P. Finkbeiner, Phillip Cargile

Abstract: Advancements in space telescopes have opened new avenues for gathering vast amounts of data on exoplanet atmosphere spectra. However, accurately extracting chemical and physical properties from these spectra poses significant challenges due to the non-linear nature of the underlying physics. This paper presents novel machine learning models developed by the AstroAI team for the Ariel Data Challe… ▽ More Advancements in space telescopes have opened new avenues for gathering vast amounts of data on exoplanet atmosphere spectra. However, accurately extracting chemical and physical properties from these spectra poses significant challenges due to the non-linear nature of the underlying physics. This paper presents novel machine learning models developed by the AstroAI team for the Ariel Data Challenge 2023, where one of the models secured the top position among 293 competitors. Leveraging Normalizing Flows, our models predict the posterior probability distribution of atmospheric parameters under different atmospheric assumptions. Moreover, we introduce an alternative model that exhibits higher performance potential than the winning model, despite scoring lower in the challenge. These findings highlight the need to reevaluate the evaluation metric and prompt further exploration of more efficient and accurate approaches for exoplanet atmosphere spectra analysis. Finally, we present recommendations to enhance the challenge and models, providing valuable insights for future applications on real observational data. These advancements pave the way for more effective and timely analysis of exoplanet atmospheric properties, advancing our understanding of these distant worlds. △ Less

Submitted 17 September, 2023; originally announced September 2023.

Comments: Conference proceeding for the ECML PKDD 2023

Numerical quantification of the wind properties of cool main sequence stars

Authors: Judy Chebly, Julián D. Alvarado-Gómez, Katja Poppenhäger, Cecilia Garraffo

Abstract: As a cool star evolves, it loses mass and angular momentum due to magnetized stellar winds which affect its rotational evolution. This change has consequences that range from the alteration of its activity to influences over the atmosphere of any orbiting planet. Despite their importance, observations constraining the properties of stellar winds in cool stars are extremely limited. Therefore, nume… ▽ More As a cool star evolves, it loses mass and angular momentum due to magnetized stellar winds which affect its rotational evolution. This change has consequences that range from the alteration of its activity to influences over the atmosphere of any orbiting planet. Despite their importance, observations constraining the properties of stellar winds in cool stars are extremely limited. Therefore, numerical simulations provide a valuable way to understand the structure and properties of these winds. In this work, we simulate the magnetized winds of 21 cool main-sequence stars (F-type to M-dwarfs), using a state-of-the-art 3D MHD code driven by observed large-scale magnetic field distributions. We perform a qualitative and quantitative characterization of our solutions, analyzing the dependencies between the driving conditions (e.g., spectral type, rotation, magnetic field strength) and the resulting stellar wind parameters (e.g., Alfvén surface size, mass loss rate, angular momentum loss rate, stellar wind speeds). We compare our models with the current observational knowledge on stellar winds in cool stars and explore the behaviour of the mass loss rate as a function of the Rossby number. Furthermore, our 3D models encompass the entire classical Habitable Zones (HZ) of all the stars in our sample. This allows us to provide the stellar wind dynamic pressure at both edges of the HZ and analyze the variations of this parameter across spectral type and orbital inclination. The results here presented could serve to inform future studies of stellar wind-magnetosphere interactions and stellar wind erosion of planetary atmospheres via ion escape processe. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 21 pages, 10 figures

QPOML: A Machine Learning Approach to Detect and Characterize Quasi-Periodic Oscillations in X-ray Binaries

Authors: Thaddaeus J. Kiker, James F. Steiner, Cecilia Garraffo, Mariano Mendez, Liang Zhang

Abstract: Astronomy is presently experiencing profound growth in the deployment of machine learning to explore large datasets. However, transient quasi-periodic oscillations (QPOs) which appear in power density spectra of many X-ray binary system observations are an intriguing phenomena heretofore not explored with machine learning. In light of this, we propose and experiment with novel methodologies for pr… ▽ More Astronomy is presently experiencing profound growth in the deployment of machine learning to explore large datasets. However, transient quasi-periodic oscillations (QPOs) which appear in power density spectra of many X-ray binary system observations are an intriguing phenomena heretofore not explored with machine learning. In light of this, we propose and experiment with novel methodologies for predicting the presence and properties of QPOs to make the first ever detections and characterizations of QPOs with machine learning models. We base our findings on raw energy spectra and processed features derived from energy spectra using an abundance of data from the NICER and RXTE space telescope archives for two black hole low mass X-ray binary sources, GRS 1915+105 and MAXI J1535-571. We advance these non-traditional methods as a foundation for using machine learning to discover global inter-object generalizations between - and provide unique insights about - energy and timing phenomena to assist with the ongoing challenge of unambiguously understanding the nature and origin of QPOs. Additionally, we have developed a publicly available Python machine learning library, QPOML, to enable further Machine Learning aided investigations into QPOs. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 18 pages, 12 figures, accepted by MNRAS

Three-dimensional, Time-dependent MHD Simulation of Disk-Magnetosphere-Stellar Wind Interaction in a T Tauri, Protoplanetary System

Authors: Ofer Cohen, Cecilia Garraffo, Jeremy Drake, Kristina Monsch, Igor Sokolov, Julian Alvarado-Gomez, Federico Fraschetti

Abstract: We present a three-dimensional, time-dependent, MHD simulation of the short-term interaction between a protoplanetary disk and the stellar corona in a T Tauri system. The simulation includes the stellar magnetic field, self-consistent coronal heating and stellar wind acceleration, and a disk rotating at sub-Keplerian velocity to induce accretion. We find that initially, as the system relaxes from… ▽ More We present a three-dimensional, time-dependent, MHD simulation of the short-term interaction between a protoplanetary disk and the stellar corona in a T Tauri system. The simulation includes the stellar magnetic field, self-consistent coronal heating and stellar wind acceleration, and a disk rotating at sub-Keplerian velocity to induce accretion. We find that initially, as the system relaxes from the assumed initial conditions, the inner part of the disk winds around and moves inward and close to the star as expected. However, the self-consistent coronal heating and stellar wind acceleration build up the original state after some time, significantly pushing the disk out beyond $10R_\star$. After this initial relaxation period, we do not find clear evidence of a strong, steady accretion flow funneled along coronal field lines, but only weak, sporadic accretion. We produce synthetic coronal X-ray line emission light curves which show flare-like increases that are not correlated with accretion events nor with heating events. These variations in the line emission flux are the result of compression and expansion due to disk-corona pressure variations. Vertical disk evaporation evolves above and below the disk. However, the disk - stellar wind boundary stays quite stable, and any disk material that reaches the stellar wind region is advected out by the stellar wind. △ Less

Submitted 24 March, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: Accepted to ApJ, 12 pages, 11 figures

Mapping the X-ray variability of GRS1915+105 with machine learning

Authors: Benjamin J. Ricketts, James F. Steiner, Cecilia Garraffo, Ronald A. Remillard, Daniela Huppenkothen

Abstract: Black hole X-ray binary systems (BHBs) contain a close companion star accreting onto a stellar-mass black hole. A typical BHB undergoes transient outbursts during which it exhibits a sequence of long-lived spectral states, each of which is relatively stable. GRS 1915+105 is a unique BHB that exhibits an unequaled number and variety of distinct variability patterns in X-rays. Many of these patterns… ▽ More Black hole X-ray binary systems (BHBs) contain a close companion star accreting onto a stellar-mass black hole. A typical BHB undergoes transient outbursts during which it exhibits a sequence of long-lived spectral states, each of which is relatively stable. GRS 1915+105 is a unique BHB that exhibits an unequaled number and variety of distinct variability patterns in X-rays. Many of these patterns contain unusual behaviour not seen in other sources. These variability patterns have been sorted into different classes based on count rate and color characteristics by Belloni et al (2000). In order to remove human decision-making from the pattern-recognition process, we employ an unsupervised machine learning algorithm called an auto-encoder to learn what classifications are naturally distinct by allowing the algorithm to cluster observations. We focus on observations taken by the Rossi X-ray Timing Explorer's Proportional Counter Array. We find that the auto-encoder closely groups observations together that are classified as similar under the Belloni et al (2000) system, but that there is reasonable grounds for defining each class as made up of components from 3 groups of distinct behaviour. △ Less

Submitted 16 May, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 20 pages, 27 figures. For associated projection code to view the interactive 3D plot, see https://github.com/bjricketts/grs1915-auto-encoder

Revisiting the space weather environment of Proxima Centauri b

Authors: Cecilia Garraffo, Julián D. Alvarado-Gómez, Ofer Cohen, Jeremy J. Drake

Abstract: Close-in planets orbiting around low-mass stars are exposed to intense energetic photon and particle radiation and harsh space weather. We have modeled such conditions for Proxima Centauri b (Garraffo et al. 2016b), a rocky planet orbiting in the habitable-zone of our closest neighboring star, finding a stellar wind pressure three orders of magnitude higher than the solar wind pressure on Earth. A… ▽ More Close-in planets orbiting around low-mass stars are exposed to intense energetic photon and particle radiation and harsh space weather. We have modeled such conditions for Proxima Centauri b (Garraffo et al. 2016b), a rocky planet orbiting in the habitable-zone of our closest neighboring star, finding a stellar wind pressure three orders of magnitude higher than the solar wind pressure on Earth. At that time, no Zeeman-Doppler observations of the surface magnetic field distribution of Proxima Cen were available and a proxy from a star with similar Rossby number to Proxima was used to drive the MHD model. Recently, the first ZDI observation of Proxima Cen became available (Klein et al. 2021). We have modeled Proxima b's space weather using this map and compared it with the results from the proxy magnetogram. We also computed models for a high-resolution synthetic magnetogram for Proxima b generated by a state-of-the-art dynamo model. The resulting space weather conditions for these three scenarios are similar with only small differences found between the models based on the ZDI observed magnetogram and the proxy. We conclude that our proxy magnetogram prescription based on Rossby number is valid, and provides a simple way to estimate stellar magnetic flux distributions when no direct observations are available. Comparison with models based on the synthetic magnetogram show that the exact magnetogram details are not important for predicting global space weather conditions of planets, reinforcing earlier conclusions that the large-scale (low-order) field dominates, and that the small-scale field does not have much influence on the ambient stellar wind. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: Accepted to ApJL, 11 pages, 3 figures

Stellar Energetic Particle Transport in the Turbulent and CME-disrupted Stellar Wind of AU~Microscopii

Authors: F. Fraschetti, J. D. Alvarado-Gómez, J. J. Drake, O. CoheN, C. Garraffo

Abstract: Energetic particles emitted by active stars are likely to propagate in astrospheric magnetized plasma turbulent and disrupted by the prior passage of energetic Coronal Mass Ejections (CMEs). We carried out test-particle simulations of $\sim$ GeV protons produced at a variety of distances from the M1Ve star AU~Microscopii by coronal flares or travelling shocks. Particles are propagated within the l… ▽ More Energetic particles emitted by active stars are likely to propagate in astrospheric magnetized plasma turbulent and disrupted by the prior passage of energetic Coronal Mass Ejections (CMEs). We carried out test-particle simulations of $\sim$ GeV protons produced at a variety of distances from the M1Ve star AU~Microscopii by coronal flares or travelling shocks. Particles are propagated within the large-scale quiescent three-dimensional magnetic field and stellar wind reconstructed from measured magnetograms, and { within the same stellar environment following passage of a $10^{36}$~erg kinetic energy CME}. In both cases, magnetic fluctuations with an isotropic power spectrum are overlayed onto the large scale stellar magnetic field and particle propagation out to the two innnermost confirmed planets is examined. In the quiescent case, the magnetic field concentrates the particles onto two regions near the ecliptic plane. After the passage of the CME, the closed field lines remain inflated and the re-shuffled magnetic field remains highly compressed, shrinking the scattering mean free path of the particles. In the direction of propagation of the CME-lobes the subsequent EP flux is suppressed. Even for a CME front propagating out of the ecliptic plane, the EP flux along the planetary orbits highly fluctuates and peaks at $\sim 2 -3$ orders of magnitude higher than the average solar value at Earth, both in the quiescent and the post-CME cases. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 19 pages, 14 figures, submitted

Space Weather-driven Variations in Ly$α$ Absorption Signatures of Exoplanet Atmospheric Escape: MHD Simulations and the Case of AU Mic b

Authors: O. Cohen, J. D. Alvarado-Gomez, J. J. Drake, L. M. Harbach, C. Garraffo, F. Fraschetti

Abstract: We simulate the space environment around AU Microscopii b and the interaction between the magnetized stellar wind with a planetary atmospheric outflow for ambient stellar wind conditions and Coronal Mass Ejection (CME) conditions. We also calculate synthetic Ly$α$ absorption due to neutral hydrogen in the ambient and the escaping planetary atmosphere affected by this interaction. We find that the… ▽ More We simulate the space environment around AU Microscopii b and the interaction between the magnetized stellar wind with a planetary atmospheric outflow for ambient stellar wind conditions and Coronal Mass Ejection (CME) conditions. We also calculate synthetic Ly$α$ absorption due to neutral hydrogen in the ambient and the escaping planetary atmosphere affected by this interaction. We find that the Ly$α$ absorption is highly variable due to the highly-varying stellar wind conditions. A strong Doppler blue-shift component is observed in the Ly$α$ profile, in contradiction to the actual escape velocity observed in the simulations themselves. This result suggest that the strong Doppler blue-shift is likely attributed to the stellar wind, not the escaping neutral atmosphere, either through its advection of neutral planetary gas, or through the creation of a fast neutral flow via charge exchange between the stellar wind ions and the planetary neutrals. Indeed, our CME simulations indicate a strong stripping of magnetospheric material from the planet, including some of the neutral escaping atmosphere. Our simulations show that the pressure around close-in exoplanets is not much lower, and may be even higher, than the pressure at the top of the planetary atmosphere. Thus, the neutral atmosphere is hydrodynamically escaping with a very small velocity ($<15~km~s^{-1}$). Moreover, our simulations show that an MHD treatment is essential in order to properly capture the coupled magnetized stellar wind and the escaping atmosphere, despite of the atmosphere being neutral. This coupling should be considered when interpreting Ly$α$observations in the context of exoplanets atmospheric escape. △ Less

Submitted 13 June, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 21 pages, accepted to ApJ

Cosmology-informed neural networks to solve the background dynamics of the Universe

Authors: Augusto T. Chantada, Susana J. Landau, Pavlos Protopapas, Claudia G. Scóccola, Cecilia Garraffo

Abstract: The field of machine learning has drawn increasing interest from various other fields due to the success of its methods at solving a plethora of different problems. An application of these has been to train artificial neural networks to solve differential equations without the need of a numerical solver. This particular application offers an alternative to conventional numerical methods, with adva… ▽ More The field of machine learning has drawn increasing interest from various other fields due to the success of its methods at solving a plethora of different problems. An application of these has been to train artificial neural networks to solve differential equations without the need of a numerical solver. This particular application offers an alternative to conventional numerical methods, with advantages such as lower memory required to store solutions, parallelization, and, in some cases, a lower overall computational cost than its numerical counterparts. In this work, we train artificial neural networks to represent a bundle of solutions of the differential equations that govern the background dynamics of the Universe for four different models. The models we have chosen are $Λ\mathrm{CDM}$, the Chevallier-Polarski-Linder parametric dark energy model, a quintessence model with an exponential potential, and the Hu-Sawicki $f(R)$ model. We use the solutions that the networks provide to perform statistical analyses to estimate the values of each model's parameters with observational data; namely, estimates of the Hubble parameter from cosmic chronometers, type Ia supernovae data from the Pantheon compilation, and measurements from baryon acoustic oscillations. The results we obtain for all models match similar estimations done in the literature using numerical solvers. In addition, we estimate the error of the solutions that the trained networks provide by comparing them with the analytical solution when there is one, or to a high-precision numerical solution when there is not. Through these estimations we find that the error of the solutions is at most $\sim1\%$ in the region of the parameter space that concerns the $95\%$ confidence regions that we find using the data, for all models and all statistical analyses performed in this work. △ Less

Submitted 20 March, 2023; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 28 pages, 10 figures, 10 tables, supplemental material available at https://github.com/at-chantada/Supplemental-Materials

Journal ref: Phys. Rev. D 107, 063523 (2023)

Simulating the Space Weather in the AU Mic System: Stellar Winds and Extreme Coronal Mass Ejections

Authors: Julián D. Alvarado-Gómez, Ofer Cohen, Jeremy J. Drake, Federico Fraschetti, Katja Poppenhäger, Cecilia Garraffo, Judy Chebly, Ekaterina Ilin, Laura Harbach, Oleg Kochukhov

Abstract: Two close-in planets have been recently found around the M-dwarf flare star AU Microscopii (AU Mic). These Neptune-sized planets (AU Mic b and c) seem to be located very close to the so-called "evaporation valley" in the exoplanet population, making this system an important target for studying atmospheric loss on exoplanets. This process, while mainly driven by the high-energy stellar radiation, w… ▽ More Two close-in planets have been recently found around the M-dwarf flare star AU Microscopii (AU Mic). These Neptune-sized planets (AU Mic b and c) seem to be located very close to the so-called "evaporation valley" in the exoplanet population, making this system an important target for studying atmospheric loss on exoplanets. This process, while mainly driven by the high-energy stellar radiation, will be strongly mediated by the space environment surrounding the planets. Here we present an investigation on this last area, performing 3D numerical modeling of the quiescent stellar wind from AU Mic, as well as time-dependent simulations describing the evolution of a highly energetic Coronal Mass Ejection (CME) event in this system. Observational constraints on the stellar magnetic field and properties of the eruption are incorporated in our models. We carry out qualitative and quantitative characterizations of the stellar wind, the emerging CMEs, as well as the expected steady and transient conditions along the orbit of both exoplanets. Our results predict an extreme space weather for AU Mic and its planets. This includes sub-Alfvénic regions for the large majority of the exoplanet orbits, very high dynamic and magnetic pressure values in quiescence (varying within $10^{2} - 10^{5}$ times the dynamic pressure experienced by the Earth), and an even harsher environment during the passage of any escaping CME associated with the frequent flaring observed in AU Mic. These space weather conditions alone pose an immense challenge for the survival of the exoplanetary atmospheres (if any) in this system. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: 15 Pages, 7 Figures, 2 Tables, Accepted for publication in The Astrophysical Journal

Coronal Mass Ejections and Exoplanets: A Numerical Perspective

Authors: Julián D. Alvarado-Gómez, Jeremy J. Drake, Ofer Cohen, Federico Fraschetti, Cecilia Garraffo, Katja Poppenhäger

Abstract: Coronal mass ejections (CMEs) are more energetic than any other class of solar phenomena. They arise from the rapid release of up to $10^{33}$ erg of magnetic energy mainly in the form of particle acceleration and bulk plasma motion. Their stellar counterparts, presumably involving much larger energies, are expected to play a fundamental role in shaping the environmental conditions around low-mass… ▽ More Coronal mass ejections (CMEs) are more energetic than any other class of solar phenomena. They arise from the rapid release of up to $10^{33}$ erg of magnetic energy mainly in the form of particle acceleration and bulk plasma motion. Their stellar counterparts, presumably involving much larger energies, are expected to play a fundamental role in shaping the environmental conditions around low-mass stars, in some cases perhaps with catastrophic consequences for planetary systems due to processes such as atmospheric erosion and depletion. Despite their importance, the direct observational evidence for stellar CMEs is almost non-existent. In this way, numerical simulations constitute extremely valuable tools to shed some light on eruptive behavior in the stellar regime. Here we review recent results obtained from realistic modeling of CMEs in active stars, highlighting their key role in the interpretation of currently available observational constraints. We include studies performed on M-dwarf stars, focusing on how emerging signatures in different wavelengths related to these events vary as a function of the magnetic properties of the star. Finally, the implications and relevance of these numerical results are discussed in the context of future characterization of host star-exoplanet systems. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: 8 pages, 3 figures, Accepted for publication in Astronomical Notes (Astronomische Nachrichten). Based on an invited review talk at the XMM-Newton Science Workshop 2021: "A High-Energy View of Exoplanets and Their Environments"

StelNet: Hierarchical Neural Network for Automatic Inference in Stellar Characterization

Authors: Cecilia Garraffo, Pavlos Protopapas, Jeremy J. Drake, Ignacio Becker, Phillip Cargile

Abstract: Characterizing the fundamental parameters of stars from observations is crucial for studying the stars themselves, their planets, and the galaxy as a whole. Stellar evolution theory predicting the properties of stars as a function of stellar age and mass enables translating observables into physical stellar parameters by fitting the observed data to synthetic isochrones. However, the complexity of… ▽ More Characterizing the fundamental parameters of stars from observations is crucial for studying the stars themselves, their planets, and the galaxy as a whole. Stellar evolution theory predicting the properties of stars as a function of stellar age and mass enables translating observables into physical stellar parameters by fitting the observed data to synthetic isochrones. However, the complexity of overlapping evolutionary tracks often makes this task numerically challenging, and with a precision that can be highly variable, depending on the area of the parameter space the observation lies in. This work presents StelNet, a Deep Neural Network trained on stellar evolutionary tracks that quickly and accurately predicts mass and age from absolute luminosity and effective temperature for stars with close to solar metallicity. The underlying model makes no assumption on the evolutionary stage and includes the pre-main sequence phase. We use bootstrapping and train many models to quantify the uncertainty of the model. To break the model's intrinsic degeneracy resulting from overlapping evolutionary paths, we also built a hierarchical model that retrieves realistic posterior probability distributions of the stellar mass and age. We further test and train StelNet using a sample of stars with well-determined masses and ages from the literature. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 17 pages, 15 figures, accepted to AJ

Stellar Winds Drive Strong Variations in Exoplanet Evaporative Outflows and Transit Absorption Signatures

Authors: Laura M. Harbach, Sofia P. Moschou, Cecilia Garraffo, Jeremy J. Drake, Julián D. Alvarado-Gómez, Ofer Cohen, Federico Fraschetti

Abstract: Stellar wind and photon radiation interactions with a planet can cause atmospheric depletion, which may have a potentially catastrophic impact on a planet's habitability. While the implications of photoevaporation on atmospheric erosion have been researched to some degree, studies of the influence of the stellar wind on atmospheric loss are in their infancy. Here, we use three-dimensional magnetoh… ▽ More Stellar wind and photon radiation interactions with a planet can cause atmospheric depletion, which may have a potentially catastrophic impact on a planet's habitability. While the implications of photoevaporation on atmospheric erosion have been researched to some degree, studies of the influence of the stellar wind on atmospheric loss are in their infancy. Here, we use three-dimensional magnetohydrodynamic simulations to model the effect of the stellar wind on the magnetosphere and outflow of a hypothetical planet, modeled to have an H-rich evaporating envelope with a pre-defined mass loss rate, orbiting in the habitable zone close to a low-mass M dwarf. We take the TRAPPIST-1 system as a prototype, with our simulated planet situated at the orbit of TRAPPIST-1e. We show that the atmospheric outflow is dragged and accelerated upon interaction with the wind, resulting in a diverse range of planetary magnetosphere morphologies and plasma distributions as local stellar wind conditions change. We consider the implications of the wind-outflow interaction on potential hydrogen Lyman-alpha (Lya) observations of the planetary atmosphere during transits. The Lya observational signatures depend strongly on the local wind conditions at the time of the observation and can be subject to considerable variation on timescales as short as an hour. Our results indicate that observed variations in exoplanet Lya transit signatures could be explained by wind-outflow interaction. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: Submitted to ApJ. 18 pages, 7 figures

MESA models with magnetic braking

Authors: Seth Gossage, Aaron Dotter, Cecilia Garraffo, Jeremy J. Drake, Stephanie Douglas, Charlie Conroy

Abstract: Two magnetic braking models are implemented in MESA for use in the MIST stellar model grids. Stars less than about 1.3 $M_{\odot}$ are observed to spin down over time through interaction with their magnetized stellar winds (i.e., magnetic braking). This is the basis for gyrochronology, and fundamental to the evolution of lower mass stars. The detailed physics behind magnetic braking are uncertain,… ▽ More Two magnetic braking models are implemented in MESA for use in the MIST stellar model grids. Stars less than about 1.3 $M_{\odot}$ are observed to spin down over time through interaction with their magnetized stellar winds (i.e., magnetic braking). This is the basis for gyrochronology, and fundamental to the evolution of lower mass stars. The detailed physics behind magnetic braking are uncertain, as are 1D stellar evolution models. Thus, we calibrate our models and compare to data from open clusters. Each braking model tested here is capable of reproducing the data, albeit with some important distinctions. The Matt et al. (2015) prescription matches the slowly rotating stars observed in open clusters, but tends to overestimate the presence of rapidly rotating stars. The Garraffo et al. (2018) prescription often produces too much angular momentum loss to accurately match the observed slow sequence for lower mass stars, but reproduces the bimodal nature of slow and rapidly rotating stars observed in open clusters fairly well. We find additional evidence that some level of mass dependency may be missing in these braking models to match the rotation periods observed in clusters older than 1 Gyr better. △ Less

Submitted 22 September, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 18 pages, 6 figures, submitted to ApJ (accepted)

Journal ref: Volume 912, 2021, id.65, 16 pp

An Earth-like stellar wind environment for Proxima Centauri c

Authors: Julián D. Alvarado-Gómez, Jeremy J. Drake, Cecilia Garraffo, Ofer Cohen, Katja Poppenhäger, Rakesh K. Yadav, Sofia P. Moschou

Abstract: A new planet has been recently discovered around Proxima Centauri. With an orbital separation of $\sim$$1.44$ au and a minimum mass of about $7$ $M_{\oplus}$, Proxima c is a prime direct imaging target for atmospheric characterization. The latter can only be performed with a good understanding of the space environment of the planet, as multiple processes can have profound effects on the atmospheri… ▽ More A new planet has been recently discovered around Proxima Centauri. With an orbital separation of $\sim$$1.44$ au and a minimum mass of about $7$ $M_{\oplus}$, Proxima c is a prime direct imaging target for atmospheric characterization. The latter can only be performed with a good understanding of the space environment of the planet, as multiple processes can have profound effects on the atmospheric structure and evolution. Here, we take one step in this direction by generating physically-realistic numerical simulations of Proxima's stellar wind, coupled to a magnetosphere and ionosphere model around Proxima c. We evaluate their expected variation due to the magnetic cycle of the host star, as well as for plausible inclination angles for the exoplanet orbit. Our results indicate stellar wind dynamic pressures comparable to present-day Earth, with a slight increase (by a factor of 2) during high activity periods of the star. A relatively weak interplanetary magnetic field at the distance of Proxima c leads to negligible stellar wind Joule heating of the upper atmosphere (about $10\%$ of the solar wind contribution on Earth) for an Earth-like planetary magnetic field ($0.3$ G). Finally, we provide an assessment of the likely extreme conditions experienced by the exoplanet candidate Proxima d, tentatively located at $0.029$ au with a minimum mass of $0.29$ $M_{\oplus}$. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 9 Pages, 4 Figures, 1 Table. Accepted for publication in The Astrophysical Journal Letters (ApJL)

The High-Energy Radiation Environment Around a 10 Gyr M Dwarf: Habitable at Last?

Authors: Kevin France, Girish Duvvuri, Hilary Egan, Tommi Koskinen, David J. Wilson, Allison Youngblood, Cynthia S. Froning, Alexander Brown, Julian D. Alvarado-Gomez, Zachory K. Berta-Thompson, Jeremy J. Drake, Cecilia Garraffo, Lisa Kaltenegger, Adam F. Kowalski, Jeffrey L. Linsky, R. O. Parke Loyd, Pablo J. D. Mauas, Yamila Miguel, J. Sebastian Pineda, Sarah Rugheimer, P. Christian Schneider, Feng Tian, Mariela Vieytes

Abstract: High levels of X-ray and UV activity on young M dwarfs may drive rapid atmospheric escape on temperate, terrestrial planets orbiting within the liquid water habitable zone. However, secondary atmospheres on planets orbiting older, less active M dwarfs may be stable and present more promising candidates for biomarker searches. We present new HST and Chandra observations of Barnard's Star (GJ 699),… ▽ More High levels of X-ray and UV activity on young M dwarfs may drive rapid atmospheric escape on temperate, terrestrial planets orbiting within the liquid water habitable zone. However, secondary atmospheres on planets orbiting older, less active M dwarfs may be stable and present more promising candidates for biomarker searches. We present new HST and Chandra observations of Barnard's Star (GJ 699), a 10 Gyr old M3.5 dwarf, acquired as part of the Mega-MUSCLES program. Despite the old age and long rotation period of Barnard's star, we observe two FUV ($δ_{130}$ $\approx$ 5000s; $E_{130}$ $\approx$ 10$^{29.5}$ erg each) and one X-ray ($E_{X}$ $\approx$ 10$^{29.2}$ erg) flares, and estimate a high-energy flare duty cycle (defined here as the fraction of the time the star is in a flare state) of $\sim$ 25\%. A 5 A - 10 $μ$m SED of GJ 699 is created and used to evaluate the atmospheric stability of a hypothetical, unmagnetized terrestrial planet in the habitable zone ($r_{HZ}$ $\sim$ 0.1 AU). Both thermal and non-thermal escape modeling indicate (1) the $quiescent$ stellar XUV flux does not lead to strong atmospheric escape: atmospheric heating rates are comparable to periods of high solar activity on modern Earth, and (2) the $flare$ environment could drive the atmosphere into a hydrodynamic loss regime at the observed flare duty cycle: sustained exposure to the flare environment of GJ 699 results in the loss of $\approx$ 87 Earth atmospheres Gyr$^{-1}$ through thermal processes and $\approx$ 3 Earth atmospheres Gyr$^{-1}$ through ion loss processes, respectively. These results suggest that if rocky planet atmospheres can survive the initial $\sim$ 5 Gyr of high stellar activity, or if a second generation atmosphere can be formed or acquired, the flare duty cycle may be the controlling stellar parameter for the stability of Earth-like atmospheres around old M stars. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Comments: Accepted to AJ

The Space Environment and Atmospheric Joule Heating of the Habitable Zone Exoplanet TOI700-d

Authors: O. Cohen, C. Garraffo, S. Moschou, J. Drake, J. Alvarado-Gomez, A. Glocer, F. Fraschetti

Abstract: We investigate the space environment conditions near the Earth-size planet TOI~700~d using a set of numerical models for the stellar corona and wind, the planetary magnetosphere, and the planetary ionosphere. We drive our simulations using a scaled-down stellar input and a scaled-up solar input in order to obtain two independent solutions. We find that for the particular parameters used in our stu… ▽ More We investigate the space environment conditions near the Earth-size planet TOI~700~d using a set of numerical models for the stellar corona and wind, the planetary magnetosphere, and the planetary ionosphere. We drive our simulations using a scaled-down stellar input and a scaled-up solar input in order to obtain two independent solutions. We find that for the particular parameters used in our study, the stellar wind conditions near the planet are not very extreme -- slightly stronger than that near the Earth in terms of the stellar wind ram pressure and the intensity of the interplanetary magnetic field. Thus, the space environment near TOI700-d may not be extremely harmful to the planetary atmosphere, assuming the planet resembles the Earth. Nevertheless, we stress that the stellar input parameters and the actual planetary parameters are unconstrained, and different parameters may result in a much greater effect on the atmosphere of TOI700-d. Finally, we compare our results to solar wind measurements in the solar system and stress that modest stellar wind conditions may not guarantee atmospheric retention of exoplanets. △ Less

Submitted 23 May, 2020; originally announced May 2020.

Comments: accepted to ApJ

Tuning the Exo-Space Weather Radio for Stellar Coronal Mass Ejections

Authors: Julián D. Alvarado-Gómez, Jeremy J. Drake, Federico Fraschetti, Cecilia Garraffo, Ofer Cohen, Christian Vocks, Katja Poppenhäger, Sofia P. Moschou, Rakesh K. Yadav, Ward B. Manchester IV

Abstract: Coronal mass ejections (CMEs) on stars other than the Sun have proven very difficult to detect. One promising pathway lies in the detection of type II radio bursts. Their appearance and distinctive properties are associated with the development of an outward propagating CME-driven shock. However, dedicated radio searches have not been able to identify these transient features in other stars. Large… ▽ More Coronal mass ejections (CMEs) on stars other than the Sun have proven very difficult to detect. One promising pathway lies in the detection of type II radio bursts. Their appearance and distinctive properties are associated with the development of an outward propagating CME-driven shock. However, dedicated radio searches have not been able to identify these transient features in other stars. Large Alfvén speeds and the magnetic suppression of CMEs in active stars have been proposed to render stellar eruptions "radio-quiet". Employing 3D magnetohydrodynamic simulations, we study here the distribution of the coronal Alfvén speed, focusing on two cases representative of a young Sun-like star and a mid-activity M-dwarf (Proxima Centauri). These results are compared with a standard solar simulation and used to characterize the shock-prone regions in the stellar corona and wind. Furthermore, using a flux-rope eruption model, we drive realistic CME events within our M-dwarf simulation. We consider eruptions with different energies to probe the regimes of weak and partial CME magnetic confinement. While these CMEs are able to generate shocks in the corona, those are pushed much farther out compared to their solar counterparts. This drastically reduces the resulting type II radio burst frequencies down to the ionospheric cutoff, which impedes their detection with ground-based instrumentation. △ Less

Submitted 11 April, 2020; originally announced April 2020.

Comments: 13 Pages, 6 Figures, 2 Tables. Accepted for publication in The Astrophysical Journal

Atmospheric Escape Processes and Planetary Atmospheric Evolution

Authors: Guillaume Gronoff, Phil Arras, Suleiman M Baraka, Jared M Bell, Gaël Cessateur, Ofer Cohen, Shannon M. Curry, Jeremy J Drake, Meredith K Elrod, Justin T. Erwin, Katherine Garcia-Sage, Cecilia Garraffo, Alex Glocer, Nicholas Gray Heavens, Kylie Lovato, Romain Maggiolo, Christopher D. Parkinson, Cyril L. Simon Wedlund, Daniel R Weimer, William B. Moore

Abstract: The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to und… ▽ More The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox. △ Less

Submitted 6 March, 2020; originally announced March 2020.

(Simulating) Coronal Mass Ejections in Active Stars

Authors: Julián D. Alvarado-Gómez, Jeremy J. Drake, Cecilia Garraffo, Sofia P. Moschou, Ofer Cohen, Rakesh K. Yadav, Federico Fraschetti

Abstract: The stellar magnetic field completely dominates the environment around late-type stars. It is responsible for driving the coronal high-energy radiation (e.g. EUV/X-rays), the development of stellar winds, and the generation transient events such as flares and coronal mass ejections (CMEs). While progress has been made for the first two processes, our understanding of the eruptive behavior in late-… ▽ More The stellar magnetic field completely dominates the environment around late-type stars. It is responsible for driving the coronal high-energy radiation (e.g. EUV/X-rays), the development of stellar winds, and the generation transient events such as flares and coronal mass ejections (CMEs). While progress has been made for the first two processes, our understanding of the eruptive behavior in late-type stars is still very limited. One example of this is the fact that despite the frequent and highly energetic flaring observed in active stars, direct evidence for stellar CMEs is almost non-existent. Here we discuss realistic 3D simulations of stellar CMEs, analyzing their resulting properties in contrast with solar eruptions, and use them to provide a common framework to interpret the available stellar observations. Additionally, we present results from the first 3D CME simulations in M-dwarf stars, with emphasis on possible observable signatures imprinted in the stellar corona. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: 7 pages, 2 figures, to appear in the proceedings of the IAU Symposium 354, Solar and Stellar Magnetic Fields: Origins and Manifestations. Eds: A. Kosovichev, K. Strassmeier, M. Jardine; based on an invited talk

Coupled MHD -- Hybrid Simulations of Space Plasmas

Authors: S. P. Moschou, I. V. Sokolov, O. Cohen, G. Toth, J. J. Drake, Z. Huang, C. Garraffo, J. D. Alvarado-Gómez, T. Gombosi

Abstract: Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of space plasmas. Hybrid codes, on the other hand, capture important physical processes, such as electric currents and effects of finite Larmor radius, but they can… ▽ More Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of space plasmas. Hybrid codes, on the other hand, capture important physical processes, such as electric currents and effects of finite Larmor radius, but they can be used locally only, since the limitations in available computational resources do not allow for their use throughout a global computational domain. In the present work, we present a new coupled scheme which allows to switch blocks in the block-adaptive grids from fluid MHD to hybrid simulations, without modifying the self-consistent computation of the electromagnetic fields acting on fluids (in MHD simulation) or charged ion macroparticles (in hybrid simulation). In this way, the hybrid scheme can refine the description in specified regions of interest without compromising the efficiency of the global MHD code. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: 13 pages, 3 figures, ASTRONUM 2019 refereed proceedings paper (in press)

Coronal response to magnetically-suppressed CME events in M-dwarf stars

Authors: Julián D. Alvarado-Gómez, Jeremy J. Drake, Sofia P. Moschou, Cecilia Garraffo, Ofer Cohen, Rakesh K. Yadav, Federico Fraschetti

Abstract: We report the results of the first state-of-the-art numerical simulations of Coronal Mass Ejections (CMEs) taking place in realistic magnetic field configurations of moderately active M-dwarf stars. Our analysis indicates that a clear, novel, and observable, coronal response is generated due to the collapse of the eruption and its eventual release into the stellar wind. Escaping CME events, weakly… ▽ More We report the results of the first state-of-the-art numerical simulations of Coronal Mass Ejections (CMEs) taking place in realistic magnetic field configurations of moderately active M-dwarf stars. Our analysis indicates that a clear, novel, and observable, coronal response is generated due to the collapse of the eruption and its eventual release into the stellar wind. Escaping CME events, weakly suppressed by the large-scale field, induce a flare-like signature in the emission from coronal material at different temperatures due to compression and associated heating. Such flare-like profiles display a distinctive temporal evolution in their Doppler shift signal (from red to blue), as the eruption first collapses towards the star and then perturbs the ambient magnetized plasma on its way outwards. For stellar fields providing partial confinement, CME fragmentation takes place, leading to rise and fall flow patterns which resemble the solar coronal rain cycle. In strongly suppressed events, the response is better described as a gradual brightening, in which the failed CME is deposited in the form of a coronal rain cloud leading to a much slower rise in the ambient high-energy flux by relatively small factors ($\sim2-3$). In all the considered cases (escaping/confined) a fractional decrease in the emission from mid-range coronal temperature plasma occurs, similar to the coronal dimming events observed on the Sun. Detection of the observational signatures of these CME-induced features requires a sensitive next generation X-ray space telescope. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: 10 Pages, 5 Figures, Accepted to The Astrophysical Journal Letters

The Stellar CME-flare relation: What do historic observations reveal?

Authors: Sofia-Paraskevi Moschou, Jeremy J. Drake, Ofer Cohen, Julián D. Alvarado-Gómez, Cecilia Garraffo, Federico Fraschetti

Abstract: Solar CMEs and flares have a statistically well defined relation, with more energetic X-ray flares corresponding to faster and more massive CMEs. How this relation extends to more magnetically active stars is a subject of open research. Here, we study the most probable stellar CME candidates associated with flares captured in the literature to date, all of which were observed on magnetically activ… ▽ More Solar CMEs and flares have a statistically well defined relation, with more energetic X-ray flares corresponding to faster and more massive CMEs. How this relation extends to more magnetically active stars is a subject of open research. Here, we study the most probable stellar CME candidates associated with flares captured in the literature to date, all of which were observed on magnetically active stars. We use a simple CME model to derive masses and kinetic energies from observed quantities, and transform associated flare data to the GOES 1--8~Å band. Derived CME masses range from $\sim 10^{15}$ to $10^{22}$~g. Associated flare X-ray energies range from $10^{31}$ to $10^{37}$~erg. Stellar CME masses as a function of associated flare energy generally lie along or below the extrapolated mean for solar events. In contrast, CME kinetic energies lie below the analogous solar extrapolation by roughly two orders of magnitude, indicating approximate parity between flare X-ray and CME kinetic energies. These results suggest that the CMEs associated with very energetic flares on active stars are more limited in terms of the ejecta velocity than the ejecta mass, possibly because of the restraining influence of strong overlying magnetic fields and stellar wind drag. Lower CME kinetic energies and velocities present a more optimistic scenario for the effects of CME impacts on exoplanets in close proximity to active stellar hosts. △ Less

Submitted 21 April, 2019; originally announced April 2019.

Comments: 23 pages, 3 tables, 4 figures, accepted by ApJ

High-Energy Photon and Particle Effects onExoplanet Atmospheres and Habitability

Authors: Jeremy J. Drake, Julián D. Alvarado-Gómez, Vladimir Airapetian, P. Wilson Cauley, Costanza Argiroffi, Matthew K. Browning, Damian J. Christian, Ofer Cohen, Lia Corrales, William Danchi, Miguel de Val-Borro, Chuanfei Dong, William Forman, Kevin France, Elena Gallo, Katherine Garcia-Sage, Cecilia Garraffo, Dawn M. Gelino, Guillaume Gronoff, H. Moritz Günther, Graham M. Harper, Raphaëlle D. Haywood, Margarita Karovska, Vinay Kashyap, Joel Kastner , et al. (19 additional authors not shown)

Abstract: It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current inst… ▽ More It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: Astro2020 Decadal Survey Science White Paper

EUV observations of cool dwarf stars

Authors: Allison Youngblood, Jeremy Drake, James Mason, Rachel Osten, Meng Jin, Adam Kowalski, Kevin France, Brian Fleming, Joel Allred, Ute Amerstorfer, Zachory Berta-Thompson, Vincent Bourrier, Luca Fossati, Cynthia Froning, Cecilia Garraffo, Guillaume Gronoff, Tommi Koskinen, Herbert Lichtenegger

Abstract: The EUV (100-912 Å) is a spectral region notoriously difficult to observe due to attenuation by neutral hydrogen gas in the interstellar medium. Despite this, hundreds to thousands of nearby stars of different spectral types and magnetic activity levels are accessible in the EUV range. The EUV probes interesting and complicated regions in the stellar atmosphere like the lower corona and transition… ▽ More The EUV (100-912 Å) is a spectral region notoriously difficult to observe due to attenuation by neutral hydrogen gas in the interstellar medium. Despite this, hundreds to thousands of nearby stars of different spectral types and magnetic activity levels are accessible in the EUV range. The EUV probes interesting and complicated regions in the stellar atmosphere like the lower corona and transition region that are inaccessible from other spectral regions. In this white paper we describe how direct EUV observations, which require a dedicated grazing-incidence observatory, cannot yet be accurately substituted with models and theory. Exploring EUV emission from cool dwarf stars in the time domain can make a major contribution to understanding stellar outer atmospheres and magnetism, and offers the clearest path toward detecting coronal mass ejections on stars other than the Sun. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey. 8 pages, 2 figures

EUV influences on exoplanet atmospheric stability and evolution

Authors: Allison Youngblood, Kevin France, Tommi Koskinen, Luca Fossati, Ute Amerstorfer, Herbert Lichtenegger, Jeremy Drake, James Mason, Brian Fleming, Joel Allred, Zachory Berta-Thompson, Vincent Bourrier, Cynthia Froning, Cecilia Garraffo, Guillaume Gronoff, Meng Jin, Adam Kowalski, Rachel Osten

Abstract: The planetary effective surface temperature alone is insufficient to characterize exoplanet atmospheres and their stability or evolution. Considering the star-planet system as a whole is necessary, and a critical component of the system is the photoionizing stellar extreme ultraviolet emission (EUV; 100-912 Å). EUV photons drive atmospheric mass loss through thermal and nonthermal processes, and a… ▽ More The planetary effective surface temperature alone is insufficient to characterize exoplanet atmospheres and their stability or evolution. Considering the star-planet system as a whole is necessary, and a critical component of the system is the photoionizing stellar extreme ultraviolet emission (EUV; 100-912 Å). EUV photons drive atmospheric mass loss through thermal and nonthermal processes, and an accurate accounting of the EUV energy deposition in a planet's energy budget is essential, especially for terrestrial habitable zone planets and close-in gaseous planets. Direct EUV observations of exoplanet host stars would require a new, dedicated observatory. Archival observations from the $\textit{EUVE}$ satellite, models, and theory alone are insufficient to accurately characterize EUV spectra of the majority of exoplanet host stars, especially for low-mass stars. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey. 7 pages, 2 figures

Stellar energetic particles in the magnetically turbulent habitable zones of TRAPPIST-1-like planetary systems

Authors: F. Fraschetti, J. J. Drake, J. D. Alvarado-Gomez, S. P. Moschou, C. Garraffo, O. Cohen

Abstract: Planets in close proximity to their parent star, such as those in the habitable zones around M dwarfs, could be subject to particularly high doses of particle radiation. We have carried out test-particle simulations of ~GeV protons to investigate the propagation of energetic particles accelerated by flares or travelling shock waves within the stellar wind and magnetic field of a TRAPPIST-1-like sy… ▽ More Planets in close proximity to their parent star, such as those in the habitable zones around M dwarfs, could be subject to particularly high doses of particle radiation. We have carried out test-particle simulations of ~GeV protons to investigate the propagation of energetic particles accelerated by flares or travelling shock waves within the stellar wind and magnetic field of a TRAPPIST-1-like system. Turbulence was simulated with small-scale magnetostatic perturbations with an isotropic power spectrum. We find that only a few percent of particles injected within half a stellar radius from the stellar surface escape, and that the escaping fraction increases strongly with increasing injection radius. Escaping particles are increasingly deflected and focused by the ambient spiralling magnetic field as the superimposed turbulence amplitude is increased. In our TRAPPIST-1-like simulations, regardless of the angular region of injection, particles are strongly focused onto two caps within the fast wind regions and centered on the equatorial planetary orbital plane. Based on a scaling relation between far-UV emission and energetic protons for solar flares applied to M dwarfs, the innermost putative habitable planet, TRAPPIST-1e, is bombarded by a proton flux up to 6 orders of magnitude larger than experienced by the present-day Earth. We note two mechanisms that could strongly limit EP fluxes from active stars: EPs from flares are contained by the stellar magnetic field; and potential CMEs that might generate EPs at larger distances also fail to escape. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 17 pages, 12 figures, ApJ in press

Breezing through the space environment of Barnard's Star b

Authors: Julián D. Alvarado-Gómez, Cecilia Garraffo, Jeremy J. Drake, Benjamin P. Brown, Jeffrey S. Oishi, Sofia P. Moschou, Ofer Cohen

Abstract: A physically realistic stellar wind model based on Alfvén wave dissipation has been used to simulate the wind from Barnard's Star and to estimate the conditions at the location of its recently discovered planetary companion. Such models require knowledge of the stellar surface magnetic field that is currently unknown for Barnard's Star. We circumvent this by considering the observed field distribu… ▽ More A physically realistic stellar wind model based on Alfvén wave dissipation has been used to simulate the wind from Barnard's Star and to estimate the conditions at the location of its recently discovered planetary companion. Such models require knowledge of the stellar surface magnetic field that is currently unknown for Barnard's Star. We circumvent this by considering the observed field distributions of three different stars that constitute admissible magnetic proxies of this object. Under these considerations, Barnard's Star b experiences less intense wind pressure than the much more close-in planet Proxima~b and the planets of the TRAPPIST-1 system. The milder wind conditions are more a result of its much greater orbital distance rather than in differences in the surface magnetic field strengths of the host stars. The dynamic pressure experienced by the planet is comparable to present-day Earth values, but it can undergo variations by factors of several during current sheet crossings in each orbit. The magnetospause standoff distance would be $\sim$\,$20 - 40$\,\% smaller than that of the Earth for an equivalent planetary magnetic field strength. △ Less

Submitted 29 March, 2019; v1 submitted 1 January, 2019; originally announced January 2019.

Comments: J. D. Alvarado-Gómez & C. Garraffo contributed equally to this work. 8 Pages, 4 Figures. Accepted to The Astrophysical Journal Letters

Synthetic Radio Imaging for Quiescent and CME-flare Scenarios

Authors: Sofia-Paraskevi Moschou, Igor Sokolov, Ofer Cohen, Jeremy J. Drake, Dmitry Borovikov, Justin C. Kasper, Julian D. Alvarado-Gomez, Cecilia Garraffo

Abstract: Radio observations grant access to a wide range of physical processes through different emission mechanisms. These processes range from thermal and quiescent to eruptive phenomena, such as shock waves and particle beams. We present a new synthetic radio imaging tool that calculates and visualizes the Bremsstrahlung radio emission. This tool works concurrently with state-of-the-art Magnetohydrodyna… ▽ More Radio observations grant access to a wide range of physical processes through different emission mechanisms. These processes range from thermal and quiescent to eruptive phenomena, such as shock waves and particle beams. We present a new synthetic radio imaging tool that calculates and visualizes the Bremsstrahlung radio emission. This tool works concurrently with state-of-the-art Magnetohydrodynamic (MHD) simulations of the solar corona using the code BATS-R-US. Our model produces results that are in good agreement with both high and low frequency observations of the solar disk. In this study, a ray-tracing algorithm is used and the radio intensity is computed along the actual curved ray trajectories. We illustrate the importance of refraction in locating the radio emitting source by comparison of the radio imaging illustrations when the line-of-sight instead of the refracted paths are considered. We are planning to incorporate non-thermal radio emission mechanisms in a future version of the radio imaging tool. △ Less

Submitted 25 September, 2018; originally announced September 2018.

Comments: 15 pages, 8 figures, Accepted by ApJ (NOTE: The current version includes extra material compared to the accepted to ApJ version.)

Exoplanet Modulation of Stellar Coronal Radio Emission

Authors: Ofer Cohen, Sofia-Paraskevi Moschou, Alex Glocer, Igor Sokolov, Tsevi Mazeh, Jeremy Drake, Cecilia Garraffo, Julian Alvarado-Gomez

Abstract: The search for exoplanets in the radio bands has been focused on detecting radio emissions produced by the interaction between magnetized planets and the stellar wind (auroral emission). Here we introduce a new tool, which is part of our MHD stellar corona model, to predict the ambient coronal radio emission and its modulations induced by a close planet. For simplicity, the present work assumes th… ▽ More The search for exoplanets in the radio bands has been focused on detecting radio emissions produced by the interaction between magnetized planets and the stellar wind (auroral emission). Here we introduce a new tool, which is part of our MHD stellar corona model, to predict the ambient coronal radio emission and its modulations induced by a close planet. For simplicity, the present work assumes that the exoplanet is stationary in the frame rotating with the stellar rotation. We explore the radio flux modulations using a limited parameter space of idealized cases by changing the magnitude of the planetary field, its polarity, the planetary orbital separation, and the strength of the stellar field. We find that the modulations induced by the planet could be significant and observable in the case of hot Jupiter planets --- above 100% modulation with respect to the ambient flux in the $10-100~MHz$ range in some cases, and 2-10% in the frequency bands above $250~MHz$ for some cases. Thus, our work indicates that radio signature of exoplanets might not be limited to low-frequency radio range. We find that the intensity modulations are sensitive to the planetary magnetic field polarity for short-orbit planets, and to the stellar magnetic field strength for all cases. The new radio tool, when applied to real systems, could provide predictions for the frequency range at which the modulations can be observed by current facilities. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: Accepted to AJ

The Magnetic Nature of the Cataclysmic Variable Period Gap

Authors: C. Garraffo, J. J. Drake, J. D. Alvarado-Gomez, S. P. Moschou, O. Cohen

Abstract: One of the most important problems in the context of cataclysmic variables (CVs) is the lack of observations of systems with periods between 2 and 3.12 hours, known as the period gap. The orbital evolution of CVs with periods shorter than those in the gap is dominated by gravitational radiation while for periods exceeding those of the gap it is dominated by magnetic braking of the secondary star.… ▽ More One of the most important problems in the context of cataclysmic variables (CVs) is the lack of observations of systems with periods between 2 and 3.12 hours, known as the period gap. The orbital evolution of CVs with periods shorter than those in the gap is dominated by gravitational radiation while for periods exceeding those of the gap it is dominated by magnetic braking of the secondary star. Spruit & Ritter (1983) showed that as periods approach 3 hours and secondary stars become fully convective a sharp decline in magnetic dynamo and braking efficiency would result in such a gap. Recent X-ray observations finding coronal magnetic energy dissipation is similar in fully convective and partly radiative M dwarfs cast this theory into doubt. In this work, we use Zeeman-Doppler imaging observations culled from the literature to show that the complexity of the surface magnetic fields of rapidly rotating M dwarfs increases with decreasing rotation period. Garraffo et al. (2018) have shown that the efficiency of angular momentum loss of cool stars declines strongly with increasing complexity of their surface magnetic field. We explore the idea of Taam & Spruit (1989) that magnetic complexity might then explain the period gap. By generating synthetic CV populations using a schematic CV evolutionary approach, we show that the CV period gap can naturally arise as a consequence of a rise in secondary star magnetic complexity near the long period edge of the gap that renders a sharp decline in their angular momentum loss rate. △ Less

Submitted 1 October, 2018; v1 submitted 28 August, 2018; originally announced August 2018.

Comments: 7 pages, 3 figures, accepted to ApJ

Suppression of Coronal Mass Ejections in active stars by an overlying large-scale magnetic field: A numerical study

Authors: Julián D. Alvarado-Gómez, Jeremy J. Drake, Ofer Cohen, Sofia P. Moschou, Cecilia Garraffo

Abstract: We present results from a set of numerical simulations aimed at exploring the mechanism of coronal mass ejection (CME) suppression in active stars by an overlying large-scale magnetic field. We use a state-of-the-art 3D magnetohydrodynamic (MHD) code which considers a self-consistent coupling between an Alfvén wave-driven stellar wind solution, and a first-principles CME model based on the eruptio… ▽ More We present results from a set of numerical simulations aimed at exploring the mechanism of coronal mass ejection (CME) suppression in active stars by an overlying large-scale magnetic field. We use a state-of-the-art 3D magnetohydrodynamic (MHD) code which considers a self-consistent coupling between an Alfvén wave-driven stellar wind solution, and a first-principles CME model based on the eruption of a flux-rope anchored to a mixed polarity region. By replicating the driving conditions used in simulations of strong solar CMEs, we show that a large-scale dipolar magnetic field of $75$ G is able to fully confine eruptions within the stellar corona. Our simulations also consider CMEs exceeding the magnetic energy used in solar studies, which are able to escape the large-scale magnetic field confinement. The analysis includes a qualitative and quantitative description of the simulated CMEs and their dynamics, which reveals a drastic reduction of the radial speed caused by the overlying magnetic field. With the aid of recent observational studies, we place our numerical results in the context of solar and stellar flaring events. In this way, we find that this particular large-scale magnetic field configuration establishes a suppression threshold around $\sim$$3 \times 10^{32}$~erg in the CME kinetic energy. Extending the solar flare-CME relations to other stars, such CME kinetic energies could be typically achieved during erupting flaring events with total energies larger than $6 \times 10^{32}$ erg (GOES class $\sim$X70). △ Less

Submitted 7 June, 2018; originally announced June 2018.

Comments: 16 Pages, 8 figures, Accepted for publication in The Astrophysical Journal

The Revolution Revolution: magnetic morphology driven spin-down

Authors: C. Garraffo, J. J. Drake, A. Dotter, J. Choi, D. J. Burke, S. P. Moschou, J. D. Alvarado-Gomez, V. L. Kashyap, O. Cohen

Abstract: Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magn… ▽ More Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magnetic flux in higher-order multipolar components of the magnetic field. In this work, we present a new model for stellar spin-down that, for the first time, accounts for the stellar surface magnetic field configuration. We show how a magnetic complexity that evolves from complex toward simple configurations as a star spins down can explain the salient features of stellar rotation evolution, including the bimodal distribution of both slow and fast rotators seen in young open clusters. △ Less

Submitted 15 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 8 pages, 4 figures, to appear in ApJ

Energy Dissipation in the Upper Atmospheres of Trappist-1 Planets

Authors: O. Cohen, A. Glocer, C. Garraffo, J. D. Drake, J. M. Bell

Abstract: We present a method to quantify the upper-limit of the energy transmitted from the intense stellar wind to the upper atmospheres of three of the Trappist-1 planets (e, f, and g). We use a formalism that treats the system as two electromagnetic regions, where the efficiency of the energy transmission between one region (the stellar wind at the planetary orbits) to the other (the planetary ionospher… ▽ More We present a method to quantify the upper-limit of the energy transmitted from the intense stellar wind to the upper atmospheres of three of the Trappist-1 planets (e, f, and g). We use a formalism that treats the system as two electromagnetic regions, where the efficiency of the energy transmission between one region (the stellar wind at the planetary orbits) to the other (the planetary ionospheres) depends on the relation between the conductances and impedances of the two regions. Since the energy flux of the stellar wind is very high at these planetary orbits, we find that for the case of high transmission efficiency (when the conductances and impedances are close in magnitude), the energy dissipation in the upper planetary atmospheres is also very large. On average, the Ohmic energy can reach $0.5-1~W/m^2$, about 1\% of the stellar irradiance and 5-15 times the EUV irradiance. Here, using constant values for the ionospheric conductance, we demonstrate that the stellar wind energy could potentially drive large atmospheric heating in terrestrial planets, as well as in hot jupiters. More detailed calculations are needed to assess the ionospheric conductance and to determine more accurately the amount of heating the stellar wind can drive in close-orbit planets. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 6 pages, 1 table, 2 figures, accepted to ApJ Letters

The Solar Wind Environment in Time

Authors: Quentin Pognan, Cecilia Garraffo, Ofer Cohen, Jeremy J. Drake

Abstract: We use magnetograms of 8 solar analogues of ages 30~Myr to 3.6~Gyr obtained from Zeeman Doppler Imaging (ZDI) and taken from the literature, together with two solar magnetograms, to drive magnetohydrodynamical (MHD) wind simulations and construct an evolutionary scenario of the solar wind environment and its angular momentum loss rate. With observed magnetograms of the radial field strength as the… ▽ More We use magnetograms of 8 solar analogues of ages 30~Myr to 3.6~Gyr obtained from Zeeman Doppler Imaging (ZDI) and taken from the literature, together with two solar magnetograms, to drive magnetohydrodynamical (MHD) wind simulations and construct an evolutionary scenario of the solar wind environment and its angular momentum loss rate. With observed magnetograms of the radial field strength as the only variant in the wind model, we find that power law model fitted to the derived angular momentum loss rate against time, $t$, results in a spin down relation $Ω\propto t^{-0.51}$, for angular speed $Ω$, which is remarkably consistent with the well-established Skumanich law $Ω\propto t^{-0.5}$. We use the model wind conditions to estimate the magnetospheric standoff distances for an Earth-like test planet situated at 1~AU for each of the stellar cases, and to obtain trends of minimum and maximum wind ram pressure and average ram pressure in the solar system through time. The wind ram pressure declines with time as $\overline{P_{ram}}\propto t^{2/3}$, amounting to a factor of 50 or so over the present lifetime of the solar system. △ Less

Submitted 16 February, 2018; v1 submitted 14 February, 2018; originally announced February 2018.

Comments: 21 pages, 11 figures, accepted for publication in ApJ, 8th February 2018

A Monster CME Obscuring A Demon Star Flare

Authors: Sofia-Paraskevi Moschou, Jeremy J. Drake, Ofer Cohen, Julian D. Alvarado-Gomez, Cecilia Garraffo

Abstract: We explore the scenario of a Coronal Mass Ejection (CME) being the cause of the observed continuous X-ray absorption of the August 30 1997 superflare on the eclipsing binary Algol (the Demon Star). The temporal decay of the absorption is consistent with absorption by a CME undergoing self-similar evolution with uniform expansion velocity. We investigate the kinematic and energetic properties of th… ▽ More We explore the scenario of a Coronal Mass Ejection (CME) being the cause of the observed continuous X-ray absorption of the August 30 1997 superflare on the eclipsing binary Algol (the Demon Star). The temporal decay of the absorption is consistent with absorption by a CME undergoing self-similar evolution with uniform expansion velocity. We investigate the kinematic and energetic properties of the CME using the ice-cream cone model for its three-dimensional structure in combination with the observed profile of the hydrogen column density decline with time. Different physically justified length scales were used that allowed us to estimate lower and upper limits of the possible CME characteristics. Further consideration of the maximum available magnetic energy in starspots leads us to quantify its mass as likely lying in the range $2\times 10^{21}$-$2\times 10^{22}$ g and kinetic energy in the range $7\times 10^{35}$-$3\times 10^{38}$ erg. The results are in reasonable agreement with extrapolated relations between flare X-ray fluence and CME mass and kinetic energy derived for solar CMEs. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: 10 pages, 5 figures, accepted for publication in ApJ

Far beyond the Sun: I. The beating magnetic heart in Horologium

Authors: Julián D. Alvarado-Gómez, Gaitee A. J. Hussain, Jeremy J. Drake, Jean-François Donati, Jorge Sanz-Forcada, Beate Stelzer, Ofer Cohen, Eliana M. Amazo-Gómez, Jason H. Grunhut, Cecilia Garraffo, Sofia P. Moschou, James Silvester, Mary E. Oksala

Abstract: A former member of the Hyades cluster, $ι$ Horologii ($ι$ Hor) is a planet-hosting Sun-like star which displays the shortest coronal activity cycle known to date ($P_{\rm cyc}\sim$ 1.6 yr). With an age of $\sim$$625$ Myr, $ι$ Hor is also the youngest star with a detected activity cycle. The study of its magnetic properties holds the potential to provide fundamental information to understand the or… ▽ More A former member of the Hyades cluster, $ι$ Horologii ($ι$ Hor) is a planet-hosting Sun-like star which displays the shortest coronal activity cycle known to date ($P_{\rm cyc}\sim$ 1.6 yr). With an age of $\sim$$625$ Myr, $ι$ Hor is also the youngest star with a detected activity cycle. The study of its magnetic properties holds the potential to provide fundamental information to understand the origin of cyclic activity and stellar magnetism in late-type stars. In this series of articles, we present the results of a comprehensive project aimed at studying the evolving magnetic field in this star and how this evolution influences its circumstellar environment. This paper summarizes the first stage of this investigation, with results from a long-term observing campaign of $ι$ Hor using ground-based high-resolution spectropolarimetry. The analysis includes precise measurements of the magnetic activity and radial velocity of the star, and their multiple time-scales of variability. In combination with values reported in the literature, we show that the long-term chromospheric activity evolution of $ι$ Hor follows a beating pattern, caused by the superposition of two periodic signals of similar amplitude at $P_{1} \simeq 1.97 \pm 0.02$ yr and $P_{2} \simeq 1.41 \pm 0.01$ yr. Additionally, using the most recent parameters for $ι$ Hor b in combination with our activity and radial velocity measurements, we find that stellar activity dominates the radial velocity residuals, making the detection of additional planets in this system challenging. Finally, we report here the first measurements of the surface longitudinal magnetic field strength of $ι$ Hor, which displays varying amplitudes within $\pm4$ G and served to estimate the rotation period of the star ($P_{\rm rot} =~7.70^{+0.18}_{-0.67}$ d). △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: 14 pages, 7 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Mottled protoplanetary disk ionization by magnetically-channeled T Tauri star energetic particles

Authors: Federico Fraschetti, Jeremy J. Drake, Ofer Cohen, Cecilia Garraffo

Abstract: The evolution of protoplanetary disks is believed to be driven largely by angular momentum transport resulting from magnetized disk winds and turbulent viscosity. The ionization of the disk that is essential for these processes has been thought due to host star coronal X-rays but could also arise from energetic particles produced by coronal flares or by travelling shock waves and advected by the s… ▽ More The evolution of protoplanetary disks is believed to be driven largely by angular momentum transport resulting from magnetized disk winds and turbulent viscosity. The ionization of the disk that is essential for these processes has been thought due to host star coronal X-rays but could also arise from energetic particles produced by coronal flares or by travelling shock waves and advected by the stellar wind. We have performed test-particle numerical simulations of energetic protons propagating into a realistic T~Tauri stellar wind, including a superposed small-scale magnetostatic turbulence. The isotropic (Kolmogorov power spectrum) turbulent component is synthesised along the individual particle trajectories. We have investigated the energy range $[0.1 - 10]$ GeV, consistent with expectations from {\it Chandra} X-ray observations of large flares on T~Tauri stars and with recent indications by the {\it Herschel} Space Observatory of a significant contribution of energetic particles to the disk ionization of young stars. In contrast with a previous theoretical study finding dominance of energetic particles over X-ray in the ionization throughout the disk, we find that the disk ionization is likely dominated by X-rays over much of its area except within narrow regions where particles are channeled onto the disk by the strongly-tangled and turbulent magnetic field. The radial thickness of such regions is $\sim 5$ stellar radii close to the star and broadens with increasing radial distance. This likely continues out to large distances from the star ($10$ AU or greater) where particles can be copiously advected and diffused by the turbulent wind. △ Less

Submitted 29 December, 2017; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 18 pages, 14 figures, ApJ, in press

The Threatening Environment of the TRAPPIST-1 Planets

Authors: Cecilia Garraffo, Jeremy J. Drake, Ofer Cohen, Julian D. Alvarado-Gomez, Sofia P. Moschou

Abstract: Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind… ▽ More Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind and energetic electromagnetic and particle radiation. Here, we use state-of-the-art 3D magnetohydrodynamic models to simulate the wind around TRAPPIST-1 and study the conditions at each planetary orbit. All planets experience a stellar wind pressure between $10^3$ and $10^5$ times the solar wind pressure on Earth. All orbits pass through wind pressure changes of an order of magnitude and most planets spend a large fraction of their orbital period in the sub-Alfvénic regime. For plausible planetary magnetic field strengths, all magnetospheres are greatly compressed and undergo much more dynamic change than that of the Earth. The planetary magnetic fields connect with the stellar radial field over much of the planetary surface, allowing direct flow of stellar wind particles onto the planetary atmosphere. These conditions could result in strong atmospheric stripping and evaporation and should be taken into account for any realistic assessment of the evolution and habitability of the TRAPPIST-1 planets. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: 12 pages, 5 figures, to appear in ApJ Letters

Giant coronal loops dominate the quiescent X-ray emission in rapidly rotating M stars

Authors: O. Cohen, R. Yadav, C. Garraffo, S. H. Saar, S. J. Wolk, V. L. Kashyap, J. J. Drake, I. Pillitteri

Abstract: Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several ki… ▽ More Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several kilo Gauss large-scale magnetic fields at high latitudes, the corona and its X-ray emission are dominated by star-size large hot loops, while the smaller, underlying colder loops are not visible much in the X-ray. Based on this result we propose that, in rapidly rotating stars, emission from such coronal structures dominates the quiescent, cooler but saturated X-ray emission. △ Less

Submitted 7 November, 2016; originally announced November 2016.

Comments: 9 pages, 4 figures, accepted to ApJ

Stellar flares and the dark energy of CMEs

Authors: Jeremy J. Drake, Ofer Cohen, Cecilia Garraffo, Vinay Kashyap

Abstract: Flares we observe on stars in white light, UV or soft X-rays are probably harbingers of coronal mass ejections (CMEs). If we use the Sun as a guide, large stellar flares will dissipate two orders of magnitude less X-ray radiative energy than the kinetic energy in the associated CME. Since coronal emission on active stars appears to be dominated by flare activity, CMEs pose a quandary for understan… ▽ More Flares we observe on stars in white light, UV or soft X-rays are probably harbingers of coronal mass ejections (CMEs). If we use the Sun as a guide, large stellar flares will dissipate two orders of magnitude less X-ray radiative energy than the kinetic energy in the associated CME. Since coronal emission on active stars appears to be dominated by flare activity, CMEs pose a quandary for understanding the fraction of their energy budget stars can spend on magnetic activity. One answer is magnetic suppression of CMEs, in which the strong large-scale fields of active stars entrap and prevent CMEs unless their free energy exceeds a critical value. The CME-less flaring active region NOAA 2192 presents a possible solar analogue of this. Monster CMEs will still exist, and have the potential to ravage planetary atmospheres. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: 6 pages, to appear in Proceedings of IAU Symposium No. 320