Source author record

Ofer Cohen

Ofer Cohen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

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

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.

preprint2021arXiv

Coronal Mass Ejections and Exoplanets: A Numerical Perspective

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.

preprint2020arXiv

Atmospheric Escape Processes and Planetary Atmospheric Evolution

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.

preprint2020arXiv

The HST PanCET Program: An Optical to Infrared Transmission Spectrum of HAT-P-32Ab

We present a 0.3-5 micron transmission spectrum of the hot Jupiter HAT-P-32Ab observed with the Space Telescope Imaging Spectrograph (STIS) and Wide Field Camera 3 (WFC3) instruments mounted on the Hubble Space Telescope, combined with Spitzer Infrared Array Camera (IRAC) photometry. The spectrum is composed of 51 spectrophotometric bins with widths ranging between 150 and 400 Å, measured to a median precision of 215 ppm. Comparisons of the observed transmission spectrum to a grid of 1D radiative-convective equilibrium models indicate the presence of clouds/hazes, consistent with previous transit observations and secondary eclipse measurements. To provide more robust constraints on the planet's atmospheric properties, we perform the first full optical to infrared retrieval analysis for this planet. The retrieved spectrum is consistent with a limb temperature of 1248$\pm$92 K, a thick cloud deck, enhanced Rayleigh scattering, and $\sim$10x solar H2O abundance. We find log($Z/Z_{\odot}$) = 2.41$_{-0.07}^{+0.06}$, in agreement with the mass-metallicity relation derived for the Solar System.

preprint2020arXiv

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

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.

preprint2020arXiv

WASP-52b. The effect of starspot correction on atmospheric retrievals

We perform atmospheric retrievals on the full optical to infrared ($0.3-5 \, μ\mathrm{m}$) transmission spectrum of the inflated hot Jupiter WASP-52b by combining HST/STIS, WFC3 IR, and Spitzer/IRAC observations. As WASP-52 is an active star which shows both out-of-transit photometric variability and starspot crossings during transits, we account for the contribution of non-occulted active regions in the retrieval. We recover a $0.1-10\times$ solar atmospheric composition, in agreement with core accretion predictions for giant planets, and a weak contribution of aerosols. We also obtain a $<3000$ K temperature for the starspots, a measure which is likely affected by the models used to fit instrumental effects in the transits, and a 5% starspot fractional coverage, compatible with expectations for the host star's spectral type. Such constraints on the planetary atmosphere and on the activity of its host star will inform future JWST GTO observations of this target.

preprint2016arXiv

Stellar flares and the dark energy of CMEs

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.

preprint2016arXiv

The Missing Magnetic Morphology Term in Stellar Rotation Evolution

This study examines the relationship between magnetic field complexity and mass and angular momentum losses. Observations of open clusters have revealed a bimodal distribution of the rotation periods of solar-like stars that has proven difficult to explain under the existing rubric of magnetic braking. Recent studies suggest that magnetic complexity can play an important role in controlling stellar spin-down rates. However, magnetic morphology is still neglected in most rotation evolution models due to the difficulty of properly accounting for its effects on wind driving and angular momentum loss. Using state-of-the-art magnetohydrodynamical magnetized wind simulations we study the effect that different distributions of the magnetic flux at different levels of geometrical complexity have on mass and angular momentum loss rates. Angular momentum loss rates depend strongly on the level of complexity of the field but are independent of the way this complexity is distributed. We deduce the analytical terms representing the magnetic field morphology dependence of mass and angular momentum loss rates. We also define a parameter that best represents complexity for real stars. As a test, we use these analytical methods to estimate mass and angular momentum loss rates for $8$ stars with observed magnetograms and compare them to the simulated results. Magnetic field complexity provides a natural physical basis for stellar rotation evolution models requiring a rapid transition between weak and strong spin-down modes.

preprint2016arXiv

The Space Weather of Proxima Centauri b

A planet orbiting in the "habitable zone" of our closest neighboring star, Proxima Centauri, has recently been discovered, and the next natural question is whether or not Proxima b is "habitable". Stellar winds are likely a source of atmospheric erosion that could be particularly severe in the case of M dwarf habitable zone planets that reside close to their parent star. Here we study the stellar wind conditions that Proxima b experiences over its orbit. We construct 3-D MHD models of the wind and magnetic field around Proxima Centauri using a surface magnetic field map for a star of the same spectral type and scaled to match the observed ~ 600 G surface magnetic field strength of Proxima. We examine the wind conditions and dynamic pressure over different plausible orbits that sample the constrained parameters of the orbit of Proxima b. For all the parameter space explored, the planet is subject to stellar wind pressures of more than 2000 times those experienced by Earth from the solar wind. During an orbit, Proxima b is also subject to pressure changes of 1 to 3 orders of magnitude on timescales of a day. Its magnetopause standoff distance consequently undergoes sudden and periodic changes by a factor of 2 to 5. Proxima b will traverse the interplanetary current sheet twice each orbit, and likely crosses into regions of subsonic wind quite frequently. These effects should be taken into account in any physically realistic assessment or prediction of its atmospheric reservoir, characteristics and loss.

preprint2015arXiv

Detecting Stars at the Galactic Centre via Synchrotron Emission

Stars orbiting within 1$\arcsec$ of the supermassive black hole in the Galactic Centre, Sgr A*, are notoriously difficult to detect due to obscuration by gas and dust. We show that some stars orbiting this region may be detectable via synchrotron emission. In such instances, a bow shock forms around the star and accelerates the electrons. We calculate that around the 10 GHz band (radio) and at 10$^{14}$ Hz (infrared) the luminosity of a star orbiting the black hole is comparable to the luminosity of Sgr A*. The strength of the synchrotron emission depends on a number of factors including the star's orbital velocity. Thus, the ideal time to observe the synchrotron flux is when the star is at pericenter. The star S2 will be $\sim 0.015\arcsec$ from Sgr A* in 2018, and is an excellent target to test our predictions.

preprint2015arXiv

Magnetic complexity as an explanation for bimodal rotation populations among young stars

Observations of young open clusters have revealed a bimodal distribution of fast and slower rotation rates that has proven difficult to explain with predictive models of spin down that depend on rotation rates alone. The Metastable Dynamo Model proposed recently by Brown, employing a stochastic transition probability from slow to more rapid spin down regimes, appears to be more successful but lacks a physical basis for such duality. Using detailed 3D MHD wind models computed for idealized multipole magnetic fields, we show that surface magnetic field complexity can provide this basis. Both mass and angular momentum losses decline sharply with increasing field complexity. Combined with observation evidence for complex field morphologies in magnetically active stars, our results support a picture in which young, rapid rotators lose angular momentum in an inefficient way because of field complexity. During this slow spin-down phase, magnetic complexity is eroded, precipitating a rapid transition from weak to strong wind coupling.

preprint2015arXiv

Quantifying the Difference Between the Flux-Tube Expansion Factor at the Source Surface and at the Alfvén Surface Using A Global MHD Model for the Solar Wind

The potential field approximation has been providing a fast, and computationally inexpensive estimation for the solar corona's global magnetic field geometry for several decades. In contrast, more physics-based global magnetohydrodynamic (MHD) models have been used for a similar purpose, while being much more computationally expensive. Here, we investigate the difference in the field geometry between a global MHD model and the potential field source surface model (PFSSM) by tracing individual magnetic field lines in the MHD model from the Alfven surface (AS), through the source surface (SS), all the way to the field line footpoint, and then back to the source surface in the PFSSM. We also compare the flux-tube expansion at two points at the SS and the AS along the same radial line. We study the effect of solar cycle variations, the order of the potential field harmonic expansion, and different magnetogram sources. We find that the flux-tube expansion factor is consistently smaller at the AS than at the SS for solar minimum and the fast solar wind, but it is consistently larger for solar maximum and the slow solar wind. We use the Wang--Sheeley--Arge (WSA) model to calculate the associated wind speed for each field line, and propagate these solar-wind speeds to 1AU. We find a more than five hours deviation in the arrival time between the two models for 20% of the field lines in the solar minimum case, and for 40% of the field lines in the solar maximum case.

preprint2015arXiv

The dependence of stellar mass and angular momentum losses on latitude and on active region and dipolar magnetic fields

Rotation evolution of late-type stars is dominated by magnetic braking and the underlying factors that control this angular momentum loss are important for the study of stellar spin-down. In this work, we study angular momentum loss as a function of two different aspects of magnetic activity using a calibrated Alfvén wave-driven magnetohydrodynamic wind model: the strengths of magnetic spots and their distribution in latitude. By driving the model using solar and modified solar surface magnetograms, we show that the topology of the field arising from the net interaction of both small-scale and large-scale field is important for spin-down rates and that angular momentum loss is not a simple function of large scale magnetic field strength. We find that changing the latitude of magnetic spots can modify mass and angular momentum loss rates by a factor of two. The general effect that causes these differences is the closing down of large-scale open field at mid- and high-latitudes by the addition of the small-scale field. These effects might give rise to modulation of mass and angular momentum loss through stellar cycles, and present a problem for ab initio attempts to predict stellar spin-down based on wind models. For all the magnetogram cases considered here, from dipoles to various spotted distributions, we find that angular momentum loss is dominated by the mass loss at mid-latitudes. The spin-down torque applied by magnetized winds therefore acts at specific latitudes and is not evenly distributed over the stellar surface, though this aspect is unlikely to be important for understanding spin-down and surface flows on stars.

preprint2014arXiv

A Grid of MHD Models for Stellar Mass Loss and Spin-down Rates of Solar Analogs

Stellar winds are believed to be the dominant factor in spin down of stars over time. However, stellar winds of solar analogs are poorly constrained due to the challenges in observing them. A great improvement has been made in the last decade in our understanding of the mechanisms responsible for the acceleration of the solar wind and in the development of numerical models for solar and stellar winds. In this paper, we present a grid of Magnetohydrodynamic (MHD) models to study and quantify the values of stellar mass-loss and angular momentum loss rates as a function of the stellar rotation period, magnetic dipole component, and coronal base density. We derive simple scaling laws for the loss rates as a function of these parameters, and constrain the possible mass-loss rate of stars with thermally-driven winds. Despite the success of our scaling law in matching the results of the model, we find a deviation between the "solar dipole" case and a real case based on solar observations that overestimates the actual solar mass-loss rate by a factor of 3. This implies that the model for stellar fields might require a further investigation with higher complexity which might include the use of a filling factor for active regions, as well as the distribution of the strength of the small-scale fields. Mass loss rates in general are largely controlled by the magnetic field strength, with the wind density varying in proportion to the confining magnetic pressure $B^2$. We also find that the mass-loss rates obtained using our grid models drop much faster with the increase in rotation period than scaling laws derived using observed stellar activity. For main-sequence solar-like stars, our scaling law for angular momentum loss vs. poloidal magnetic field strength retrieves the well-known Skumanich decline of angular velocity with time.

preprint2014arXiv

Magnetic Modulation of Stellar Angular Momentum Loss

Angular Momentum Loss is important for understanding astrophysical phenomena such as stellar rotation, magnetic activity, close binaries, and cataclysmic variables. Magnetic breaking is the dominant mechanism in the spin down of young late-type stars. We have studied angular momentum loss as a function of stellar magnetic activity. We argue that the complexity of the field and its latitudinal distribution are crucial for angular momentum loss rates. In this work we discuss how angular momentum is modulated by magnetic cycles, and how stellar spin down is not just a simple function of large scale magnetic field strength.

preprint2013arXiv

Implications of mass and energy loss due to coronal mass ejections on magnetically-active stars

Analysis of a database of solar coronal mass ejections (CMEs) and associated flares over the period 1996-2007 finds well-behaved power law relationships between the 1-8 AA flare X-ray fluence and CME mass and kinetic energy. We extrapolate these relationships to lower and higher flare energies to estimate the mass and energy loss due to CMEs from stellar coronae, assuming that the observed X-ray emission of the latter is dominated by flares with a frequency as a function of energy dn/dE=kE^-alpha. For solar-like stars at saturated levels of X-ray activity, the implied losses depend fairly weakly on the assumed value of alpha and are very large: M_dot ~ 5x10^-10 M_sun/yr and E_dot ~ 0.1L_sun. In order to avoid such large energy requirements, either the relationships between CME mass and speed and flare energy must flatten for X-ray fluence >~ 10^31 erg, or the flare-CME association must drop significantly below 1 for more energetic events. If active coronae are dominated by flares, then the total coronal energy budget is likely to be up to an order of magnitude larger than the canonical 10^-3 L_bol X-ray saturation threshold. This raises the question of what is the maximum energy a magnetic dynamo can extract from a star? For an energy budget of 1% of L_bol, the CME mass loss rate is about 5x10^-11 M_sun/yr.

preprint2011arXiv

Modeling the Young Sun's Solar Wind and its Interaction with Earth's Paleomagnetosphere

We present a focused parameter study of solar wind - magnetosphere interaction for the young Sun and Earth, $~3.5$ Ga ago, that relies on magnetohydrodynamic (MHD) simulations for both the solar wind and the magnetosphere. By simulating the quiescent young Sun and its wind we are able to propagate the MHD simulations up to Earth's magnetosphere and obtain a physically realistic solar forcing of it. We assess how sensitive the young solar wind is to changes in the coronal base density, sunspot placement and magnetic field strength, dipole magnetic field strength and the Sun's rotation period. From this analysis we obtain a range of plausible solar wind conditions the paleomagnetosphere may have been subject to. Scaling relationships from the literature suggest that a young Sun would have had a mass flux different from the present Sun. We evaluate how the mass flux changes with the aforementioned factors and determine the importance of this and several other key solar and magnetospheric variables with respect to their impact on the paleomagnetosphere. We vary the solar wind speed, density, interplanetary magnetic field strength and orientation as well as Earth's dipole magnetic field strength and tilt in a number of steady-state scenarios that are representative of young Sun-Earth interaction. This study is done as a first step of a more comprehensive effort towards understanding the implications of Sun-Earth interaction for planetary atmospheric evolution.

preprint2011arXiv

The Independency of Stellar Mass Loss Rates on stellar X-ray Luminosity and Activity Level based on Solar X-ray flux and Solar Wind Observations

Stellar mass loss rates are an important input ingredient for stellar evolution models since they determine stellar evolution parameters such as stellar spin-down and increase in stellar luminosity through the lifetime of a star. Due to the lack of direct observations of stellar winds from Sun-like stars stellar X-ray luminosity and stellar level of activity are commonly used as a proxy for estimating stellar mass loss rates. However, such an intuitive activity --- mass loss rate relation is not well defined. In this paper, I study the mass loss rate of the Sun as a function of its activity level. I compare in-situ solar wind measurements with the solar activity level represented by the solar X-ray flux. I find no clear dependency of the solar mass flux on solar X-ray flux. Instead, the solar mass loss rate is scattered around an average value of $2\cdot10^{-14}\;M_\odot\;yr^{-1}$. This independency of the mass loss rate on level of activity can be explained by the fact that the activity level is governed by the large modulations in the solar close magnetic flux, while the mass loss rate is governed by the rather constant open magnetic flux. I derive a simple expression for stellar mass loss rates as a function of the stellar ambient weak magnetic field, the stellar radius, the stellar escape velocity, and the average height of the Alfvén surface. This expression predicts stellar mass loss rates of $10^{-15}-10^{-12}\;M_\odot\;yr^{-1}$ for Sun-like stars.

preprint2009arXiv

Numerical Simulation of an EUV Coronal Wave Based on the February 13, 2009 CME Event Observed by STEREO

On 13 February 2009, a coronal wave -- CME -- dimming event was observed in quadrature by the STEREO spacecraft. We analyze this event using a three-dimensional, global magnetohydrodynamic (MHD) model for the solar corona. The numerical simulation is driven and constrained by the observations, and indicates where magnetic reconnection occurs between the expanding CME core and surrounding environment. We focus primarily on the lower corona, extending out to $3R_{\odot}$; this range allows simultaneous comparison with both EUVI and COR1 data. Our simulation produces a diffuse coronal bright front remarkably similar to that observed by STEREO/EUVI at 195 Å. It is made up of \emph{two} components, and is the result of a combination of both wave and non-wave mechanisms. The CME becomes large-scale quite low ($<$ 200 Mm) in the corona. It is not, however, an inherently large-scale event; rather, the expansion is facilitated by magnetic reconnection between the expanding CME core and the surrounding magnetic environment. In support of this, we also find numerous secondary dimmings, many far from the initial CME source region. Relating such dimmings to reconnecting field lines within the simulation provides further evidence that CME expansion leads to the "opening" of coronal field lines on a global scale. Throughout the CME expansion, the coronal wave maps directly to the CME footprint. Our results suggest that the ongoing debate over the "true" nature of diffuse coronal waves may be mischaracterized. It appears that \emph{both} wave and non-wave models are required to explain the observations and understand the complex nature of these events.