@article{ShpritsKellermanAseevetal.2017, author = {Shprits, Yuri and Kellerman, Adam C . and Aseev, Nikita and Drozdov, Alexander and Michaelis, Ingo}, title = {Multi-MeV electron loss in the heart of the radiation belts}, series = {Geophysical research letters}, volume = {44}, journal = {Geophysical research letters}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2016GL072258}, pages = {1204 -- 1209}, year = {2017}, abstract = {Significant progress has been made in recent years in understanding acceleration mechanisms in the Earth's radiation belts. In particular, a number of studies demonstrated the importance of the local acceleration by analyzing the radial profiles of phase space density (PSD) and observing building up peaks in PSD. In this study, we focus on understanding of the local loss using very similar tools. The profiles of PSD for various values of the first adiabatic invariants during the previously studied 17 January 2013 storm are presented and discussed. The profiles of PSD show clear deepening minimums consistent with the scattering by electromagnetic ion cyclotron waves. Long-term evolution shows that local minimums in PSD can persist for relatively long times. During considered interval of time the deepening minimums were observed around L* = 4 during 17 January 2013 storm and around L* = 3.5 during 1 March 2013 storm. This study shows a new method that can help identify the location, magnitude, and time of the local loss and will help quantify local loss in the future. This study also provides additional clear and definitive evidence that local loss plays a major role for the dynamics of the multi-MeV electrons.}, language = {en} } @phdthesis{Maiti2023, author = {Maiti, Snehanshu}, title = {Magnetohydrodynamic turbulence and cosmic ray transport}, doi = {10.25932/publishup-58903}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-589030}, school = {Universit{\"a}t Potsdam}, pages = {vii, 81}, year = {2023}, abstract = {The first part of the thesis studies the properties of fast mode in magneto hydro-dynamic (MHD) turbulence. 1D and 3D numerical simulations are carried out to generate decaying fast mode MHD turbulence. The injection of waves are carried out in a collinear and isotropic fashion to generate fast mode turbulence. The properties of fast mode turbulence are analyzed by studying their energy spectral density, 2D structure functions and energy decay/cascade time. The injection wave vector is varied to study the dependence of the above properties on the injection wave vectors. The 1D energy spectrum obtained for the velocity and magnetic fields has 𝐸 (𝑘) ∝ 𝑘-2. The 2D energy spectrum and 2D structure functions in parallel and perpendicular directions shows that fast mode turbulence generated is isotropic in nature. The cascade/decay rate of fast mode MHD turbulence is proportional to 𝑘-0.5 for different kinds of wave vector injection. Simulations are also carried out in 1D and 3D to compare balanced and imbalanced turbulence. The results obtained shows that while 1D imbalanced turbulence decays faster than 1D balanced turbulence, there is no difference in the decay of 3D balanced and imbalanced turbulence for the current resolution of 512 grid points. "The second part of the thesis studies cosmic ray (CR) transport in driven MHD turbulence and is strongly dependent on it's properties. Test particle simulations are carried out to study CR interaction with both total MHD turbulence and decomposed MHD modes. The spatial diffusion coefficients and the pitch angle scattering diffusion coefficients are calculated from the test particle trajectories in turbulence. The results confirms that the fast modes dominate the CR propagation, whereas Alfv{\´e}n, slow modes are much less efficient with similar pitch angle scattering rates. The cross field transport on large and small scales are investigated next. On large/global scales, normal diffusion is observed and the diffusion coefficient is suppressed by 𝑀𝜁𝐴 compared to the parallel diffusion coefficients, with 𝜁 closer to 4 in Alfv{\´e}n modes than that in total turbulence as theoretically expected. For the CR transport on scales smaller than the turbulence injection scale 𝐿, both the local and global magnetic reference frames are adopted. Super diffusion is observed on such small scales in all the cases. Particularly, CR transport in Alfv{\´e}n modes show clear Richardson diffusion in the local reference frame. The diffusion transition smoothly from the Richardson's one with index 1.5 to normal diffusion as particle's mean free path decreases from 𝜆∥ ≫ 𝐿 to 𝜆∥ ≪ 𝐿. These results have broad applications to CRs in various astrophysical environments".}, language = {en} } @article{AdolfsHoqueShprits2022, author = {Adolfs, Marjolijn and Hoque, Mohammed Mainul and Shprits, Yuri}, title = {Storm-time relative total electron content modelling using machine learning techniques}, series = {Remote sensing}, volume = {14}, journal = {Remote sensing}, number = {23}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs14236155}, pages = {17}, year = {2022}, abstract = {Accurately predicting total electron content (TEC) during geomagnetic storms is still a challenging task for ionospheric models. In this work, a neural-network (NN)-based model is proposed which predicts relative TEC with respect to the preceding 27-day median TEC, during storm time for the European region (with longitudes 30 degrees W-50 degrees E and latitudes 32.5 degrees N-70 degrees N). The 27-day median TEC (referred to as median TEC), latitude, longitude, universal time, storm time, solar radio flux index F10.7, global storm index SYM-H and geomagnetic activity index Hp30 are used as inputs and the output of the network is the relative TEC. The relative TEC can be converted to the actual TEC knowing the median TEC. The median TEC is calculated at each grid point over the European region considering data from the last 27 days before the storm using global ionosphere maps (GIMs) from international GNSS service (IGS) sources. A storm event is defined when the storm time disturbance index Dst drops below 50 nanotesla. The model was trained with storm-time relative TEC data from the time period of 1998 until 2019 (2015 is excluded) and contains 365 storms. Unseen storm data from 33 storm events during 2015 and 2020 were used to test the model. The UQRG GIMs were used because of their high temporal resolution (15 min) compared to other products from different analysis centers. The NN-based model predictions show the seasonal behavior of the storms including positive and negative storm phases during winter and summer, respectively, and show a mixture of both phases during equinoxes. The model's performance was also compared with the Neustrelitz TEC model (NTCM) and the NN-based quiet-time TEC model, both developed at the German Aerospace Agency (DLR). The storm model has a root mean squared error (RMSE) of 3.38 TEC units (TECU), which is an improvement by 1.87 TECU compared to the NTCM, where an RMSE of 5.25 TECU was found. This improvement corresponds to a performance increase by 35.6\%. The storm-time model outperforms the quiet-time model by 1.34 TECU, which corresponds to a performance increase by 28.4\% from 4.72 to 3.38 TECU. The quiet-time model was trained with Carrington averaged TEC and, therefore, is ideal to be used as an input instead of the GIM derived 27-day median. We found an improvement by 0.8 TECU which corresponds to a performance increase by 17\% from 4.72 to 3.92 TECU for the storm-time model using the quiet-time-model predicted TEC as an input compared to solely using the quiet-time model.}, language = {en} } @article{SmirnovBerrendorfShpritsetal.2020, author = {Smirnov, Artem and Berrendorf, Max and Shprits, Yuri and Kronberg, Elena A. and Allison, Hayley J. and Aseev, Nikita and Zhelavskaya, Irina and Morley, Steven K. and Reeves, Geoffrey D. and Carver, Matthew R. and Effenberger, Frederic}, title = {Medium energy electron flux in earth's outer radiation belt (MERLIN)}, series = {Space weather : the international journal of research and applications}, volume = {18}, journal = {Space weather : the international journal of research and applications}, number = {11}, publisher = {American geophysical union, AGU}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2020SW002532}, pages = {20}, year = {2020}, abstract = {The radiation belts of the Earth, filled with energetic electrons, comprise complex and dynamic systems that pose a significant threat to satellite operation. While various models of electron flux both for low and relativistic energies have been developed, the behavior of medium energy (120-600 keV) electrons, especially in the MEO region, remains poorly quantified. At these energies, electrons are driven by both convective and diffusive transport, and their prediction usually requires sophisticated 4D modeling codes. In this paper, we present an alternative approach using the Light Gradient Boosting (LightGBM) machine learning algorithm. The Medium Energy electRon fLux In Earth's outer radiatioN belt (MERLIN) model takes as input the satellite position, a combination of geomagnetic indices and solar wind parameters including the time history of velocity, and does not use persistence. MERLIN is trained on >15 years of the GPS electron flux data and tested on more than 1.5 years of measurements. Tenfold cross validation yields that the model predicts the MEO radiation environment well, both in terms of dynamics and amplitudes o f flux. Evaluation on the test set shows high correlation between the predicted and observed electron flux (0.8) and low values of absolute error. The MERLIN model can have wide space weather applications, providing information for the scientific community in the form of radiation belts reconstructions, as well as industry for satellite mission design, nowcast of the MEO environment, and surface charging analysis.}, language = {en} } @article{CervantesVillaShpritsAseevetal.2019, author = {Cervantes Villa, Juan Sebastian and Shprits, Yuri and Aseev, Nikita and Drozdov, Alexander and Castillo Tibocha, Angelica Maria and Stolle, Claudia}, title = {Identifying radiation belt electron source and loss processes by assimilating spacecraft data in a three-dimensional diffusion model}, series = {Journal of geophysical research : Space physics}, volume = {125}, journal = {Journal of geophysical research : Space physics}, number = {1}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2019JA027514}, pages = {16}, year = {2019}, abstract = {Data assimilation aims to blend incomplete and inaccurate data with physics-based dynamical models. In the Earth's radiation belts, it is used to reconstruct electron phase space density, and it has become an increasingly important tool in validating our current understanding of radiation belt dynamics, identifying new physical processes, and predicting the near-Earth hazardous radiation environment. In this study, we perform reanalysis of the sparse measurements from four spacecraft using the three-dimensional Versatile Electron Radiation Belt diffusion model and a split-operator Kalman filter over a 6-month period from 1 October 2012 to 1 April 2013. In comparison to previous works, our 3-D model accounts for more physical processes, namely, mixed pitch angle-energy diffusion, scattering by Electromagnetic Ion Cyclotron waves, and magnetopause shadowing. We describe how data assimilation, by means of the innovation vector, can be used to account for missing physics in the model. We use this method to identify the radial distances from the Earth and the geomagnetic conditions where our model is inconsistent with the measured phase space density for different values of the invariants mu and K. As a result, the Kalman filter adjusts the predictions in order to match the observations, and we interpret this as evidence of where and when additional source or loss processes are active. The current work demonstrates that 3-D data assimilation provides a comprehensive picture of the radiation belt electrons and is a crucial step toward performing reanalysis using measurements from ongoing and future missions.}, language = {en} } @article{SmirnovShpritsAllisonetal.2022, author = {Smirnov, Artem and Shprits, Yuri and Allison, Hayley and Aseev, Nikita and Drozdov, Alexander and Kollmann, Peter and Wang, Dedong and Saikin, Anthony}, title = {An empirical model of the equatorial electron pitch angle distributions in earth's outer radiation belt}, series = {Space Weather: the International Journal of Research and Applications}, volume = {20}, journal = {Space Weather: the International Journal of Research and Applications}, number = {9}, publisher = {American Geophysical Union}, address = {Washington, DC}, issn = {1542-7390}, doi = {10.1029/2022SW003053}, pages = {17}, year = {2022}, abstract = {In this study, we present an empirical model of the equatorial electron pitch angle distributions (PADs) in the outer radiation belt based on the full data set collected by the Magnetic Electron Ion Spectrometer (MagEIS) instrument onboard the Van Allen Probes in 2012-2019. The PADs are fitted with a combination of the first, third and fifth sine harmonics. The resulting equation resolves all PAD types found in the outer radiation belt (pancake, flat-top, butterfly and cap PADs) and can be analytically integrated to derive omnidirectional flux. We introduce a two-step modeling procedure that for the first time ensures a continuous dependence on L, magnetic local time and activity, parametrized by the solar wind dynamic pressure. We propose two methods to reconstruct equatorial electron flux using the model. The first approach requires two uni-directional flux observations and is applicable to low-PA data. The second method can be used to reconstruct the full equatorial PADs from a single uni- or omnidirectional measurement at off-equatorial latitudes. The model can be used for converting the long-term data sets of electron fluxes to phase space density in terms of adiabatic invariants, for physics-based modeling in the form of boundary conditions, and for data assimilation purposes.}, language = {en} } @article{WalkerBoyntonShpritsetal.2022, author = {Walker, Simon N. and Boynton, Richard J. and Shprits, Yuri and Balikhin, Michael A. and Drozdov, Alexander}, title = {Forecast of the energetic electron environment of the radiation belts}, series = {Space Weather: The International Journal of Research and Applications}, volume = {20}, journal = {Space Weather: The International Journal of Research and Applications}, number = {12}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2022SW003124}, pages = {21}, year = {2022}, abstract = {Different modeling methodologies possess different strengths and weakness. For instance, data based models may provide superior accuracy but have a limited spatial coverage while physics based models may provide lower accuracy but provide greater spatial coverage. This study investigates the coupling of a data based model of the electron fluxes at geostationary orbit (GEO) with a numerical model of the radiation belt region to improve the resulting forecasts/pastcasts of electron fluxes over the whole radiation belt region. In particular, two coupling methods are investigated. The first assumes an average value for L* for GEO, namely LGEO* L-GEO* = 6.2. The second uses a value of L* that varies with geomagnetic activity, quantified using the Kp index. As the terrestrial magnetic field responds to variations in geomagnetic activity, the value of L* will vary for a specific location. In this coupling method, the value of L* is calculated using the Kp driven Tsyganenko 89c magnetic field model for field line tracing. It is shown that this addition can result in changes in the initialization of the parameters at the Versatile Electron Radiation Belt model outer boundary. Model outputs are compared to Van Allen Probes MagEIS measurements of the electron fluxes in the inner magnetosphere for the March 2015 geomagnetic storm. It is found that the fixed LGEO* L-GEO* coupling method produces a more realistic forecast.}, language = {en} } @article{ProlSmirnovHoqueetal.2022, author = {Prol, Fabricio S. and Smirnov, Artem G. and Hoque, M. Mainul and Shprits, Yuri}, title = {Combined model of topside ionosphere and plasmasphere derived from radio-occultation and Van Allen Probes data}, series = {Scientific reports}, volume = {12}, journal = {Scientific reports}, number = {1}, publisher = {Macmillan Publishers Limited, part of Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-022-13302-1}, pages = {11}, year = {2022}, abstract = {In the last years, electron density profile functions characterized by a linear dependence on the scale height showed good results when approximating the topside ionosphere. The performance above 800 km, however, is not yet well investigated. This study investigates the capability of the semi-Epstein functions to represent electron density profiles from the peak height up to 20,000 km. Electron density observations recorded by the Van Allen Probes were used to resolve the scale height dependence in the plasmasphere. It was found that the linear dependence of the scale height in the topside ionosphere cannot be directly used to extrapolate profiles above 800 km. We find that the dependence of scale heights on altitude is quadratic in the plasmasphere. A statistical model of the scale heights is therefore proposed. After combining the topside ionosphere and plasmasphere by a unified model, we have obtained good estimations not only in the profile shapes, but also in the Total Electron Content magnitude and distributions when compared to actual measurements from 2013, 2014, 2016 and 2017. Our investigation shows that Van Allen Probes can be merged to radio-occultation data to properly represent the upper ionosphere and plasmasphere by means of a semi-Epstein function.}, language = {en} } @article{HaasShpritsAllisonetal.2022, author = {Haas, Bernhard and Shprits, Yuri and Allison, Hayley and Wutzig, Michael and Wang, Dedong}, title = {Which parameter controls ring current electron dynamics}, series = {Frontiers in astronomy and space sciences}, volume = {9}, journal = {Frontiers in astronomy and space sciences}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-987X}, doi = {10.3389/fspas.2022.911002}, pages = {11}, year = {2022}, abstract = {Predicting the electron population of Earth's ring current during geomagnetic storms still remains a challenging task. In this work, we investigate the sensitivity of 10 keV ring current electrons to different driving processes, parameterised by the Kp index, during several moderate and intense storms. Results are validated against measurements from the Van Allen Probes satellites. Perturbing the Kp index allows us to identify the most dominant processes for moderate and intense storms respectively. We find that during moderate storms (Kp < 6) the drift velocities mostly control the behaviour of low energy electrons, while loss from wave-particle interactions is the most critical parameter for quantifying the evolution of intense storms (Kp > 6). Perturbations of the Kp index used to drive the boundary conditions at GEO and set the plasmapause location only show a minimal effect on simulation results over a limited L range. It is further shown that the flux at L \& SIM; 3 is more sensitive to changes in the Kp index compared to higher L shells, making it a good proxy for validating the source-loss balance of a ring current model.}, language = {en} } @phdthesis{Aseev2020, author = {Aseev, Nikita}, title = {Modeling and understanding dynamics of charged particles in the Earth's inner magnetosphere}, doi = {10.25932/publishup-47921}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-479211}, school = {Universit{\"a}t Potsdam}, pages = {xxii, 154}, year = {2020}, abstract = {The Earth's inner magnetosphere is a very dynamic system, mostly driven by the external solar wind forcing exerted upon the magnetic field of our planet. Disturbances in the solar wind, such as coronal mass ejections and co-rotating interaction regions, cause geomagnetic storms, which lead to prominent changes in charged particle populations of the inner magnetosphere - the plasmasphere, ring current, and radiation belts. Satellites operating in the regions of elevated energetic and relativistic electron fluxes can be damaged by deep dielectric or surface charging during severe space weather events. Predicting the dynamics of the charged particles and mitigating their effects on the infrastructure is of particular importance, due to our increasing reliance on space technologies. The dynamics of particles in the plasmasphere, ring current, and radiation belts are strongly coupled by means of collisions and collisionless interactions with electromagnetic fields induced by the motion of charged particles. Multidimensional numerical models simplify the treatment of transport, acceleration, and loss processes of these particles, and allow us to predict how the near-Earth space environment responds to solar storms. The models inevitably rely on a number of simplifications and assumptions that affect model accuracy and complicate the interpretation of the results. In this dissertation, we quantify the processes that control electron dynamics in the inner magnetosphere, paying particular attention to the uncertainties of the employed numerical codes and tools. We use a set of convenient analytical solutions for advection and diffusion equations to test the accuracy and stability of the four-dimensional Versatile Electron Radiation Belt (VERB-4D) code. We show that numerical schemes implemented in the code converge to the analytical solutions and that the VERB-4D code demonstrates stable behavior independent of the assumed time step. The order of the numerical scheme for the convection equation is demonstrated to affect results of ring current and radiation belt simulations, and it is crucially important to use high-order numerical schemes to decrease numerical errors in the model. Using the thoroughly tested VERB-4D code, we model the dynamics of the ring current electrons during the 17 March 2013 storm. The discrepancies between the model and observations above 4.5 Earth's radii can be explained by uncertainties in the outer boundary conditions. Simulation results indicate that the electrons were transported from the geostationary orbit towards the Earth by the global-scale electric and magnetic fields. We investigate how simulation results depend on the input models and parameters. The model is shown to be particularly sensitive to the global electric field and electron lifetimes below 4.5 Earth's radii. The effects of radial diffusion and subauroral polarization streams are also quantified. We developed a data-assimilative code that blends together a convection model of energetic electron transport and loss and Van Allen Probes satellite data by means of the Kalman filter. We show that the Kalman filter can correct model uncertainties in the convection electric field, electron lifetimes, and boundary conditions. It is also demonstrated how the innovation vector - the difference between observations and model prediction - can be used to identify physical processes missing in the model of energetic electron dynamics. We computed radial profiles of phase space density of ultrarelativistic electrons, using Van Allen Probes measurements. We analyze the shape of the profiles during geomagnetically quiet and disturbed times and show that the formation of new local minimums in the radial profiles coincides with the ground observations of electromagnetic ion-cyclotron (EMIC) waves. This correlation indicates that EMIC waves are responsible for the loss of ultrarelativistic electrons from the heart of the outer radiation belt into the Earth's atmosphere.}, language = {en} } @article{KimShprits2019, author = {Kim, Kyung-Chan and Shprits, Yuri}, title = {Statistical Analysis of Hiss Waves in Plasmaspheric Plumes Using Van Allen Probe Observations}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA026458}, pages = {1904 -- 1915}, year = {2019}, abstract = {Plasmaspheric hiss waves commonly observed in high-density regions in the Earth's magnetosphere are known to be one of the main contributors to the loss of radiation belt electrons. There has been a lot of effort to investigate the distributions of hiss waves in the plasmasphere, while relatively little attention has been given to those in the plasmaspheric plume. In this study, we present for the first time a statistical analysis of the occurrence and the spatial distribution of wave amplitudes and wave normal angles for hiss waves in plumes using Van Allen Probes observations during the period of October 2012 to December 2016. Statistical results show that a wide range of hiss wave amplitudes in plumes from a few picotesla to >100 pT is observed, but a modest (<20 pT) wave amplitude is more commonly observed regardless of geomagnetic activity in both the midnight-to-dawn and dusk sector. By contrast, stronger amplitude hiss occurs preferentially during geomagnetically active times in the dusk sector. The wave normal angles are distributed over a broad range from 0° to 90° with a bimodal distribution: a quasi-field-aligned population (<20°) with an occurrence rate of <60\% and an oblique one (>50°) with a relative low occurrence rate of ≲20\%. Therefore, from a statistical point of view, we confirm that the hiss intensity (a few tens of picotesla) and field-aligned hiss wave adopted in previous simulation studies are a reasonable assumption but stress that the activity dependence of the wave amplitude should be considered.}, language = {en} } @article{ZhuShpritsSpasojevicetal.2019, author = {Zhu, Hui and Shprits, Yuri and Spasojevic, M. and Drozdov, Alexander}, title = {New hiss and chorus waves diffusion coefficient parameterizations from the Van Allen Probes and their effect on long-term relativistic electron radiation-belt VERB simulations}, series = {Journal of Atmospheric and Solar-Terrestrial Physics}, volume = {193}, journal = {Journal of Atmospheric and Solar-Terrestrial Physics}, publisher = {Elsevier}, address = {Oxford}, issn = {1364-6826}, doi = {10.1016/j.jastp.2019.105090}, pages = {13}, year = {2019}, abstract = {New wave frequency and amplitude models for the nightside and dayside chorus waves are built based on measurements from the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) instrument onboard the Van Allen Probes. The corresponding 3D diffusion coefficients are systematically obtained. Compared with previous commonly-used (typical) parameterizations, the new parameterizations result in differences in diffusion rates that depend on the energy and pitch angle. Furthermore, one-year 3D diffusive simulations are performed using the Versatile Electron Radiation Belt (VERB) code. Both typical and new wave parameterizations simulation results are in a good agreement with observations at 0.9 MeV. However, the new parameterizations for nightside chorus better reproduce the observed electron fluxes. These parameterizations will be incorporated into future modeling efforts.}, language = {en} } @article{DrozdovAseevEffenbergeretal.2019, author = {Drozdov, Alexander and Aseev, Nikita and Effenberger, Frederic and Turner, Drew L. and Saikin, Anthony and Shprits, Yuri}, title = {Storm Time Depletions of Multi-MeV Radiation Belt Electrons Observed at Different Pitch Angles}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {11}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2019JA027332}, pages = {8943 -- 8953}, year = {2019}, abstract = {During geomagnetic storms, the rapid depletion of the high-energy (several MeV) outer radiation belt electrons is the result of loss to the interplanetary medium through the magnetopause, outward radial diffusion, and loss to the atmosphere due to wave-particle interactions. We have performed a statistical study of 110 storms using pitch angle resolved electron flux measurements from the Van Allen Probes mission and found that inside of the radiation belt (L* = 3 - 5) the number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 30 degrees is higher than number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 75 degrees. We conclude that this result is consistent with electron scattering by whistler and electromagnetic ion cyclotron waves. At the outer edge of the radiation belt (L* >= 5.2) the number of storms that result in depletion is also large (similar to 40-50\%), emphasizing the significance of the magnetopause shadowing effect and outward radial transport.}, language = {en} } @article{CastilloShpritsGanushkinaetal.2019, author = {Castillo, Angelica M. and Shprits, Yuri and Ganushkina, Natalia and Drozdov, Alexander and Aseev, Nikita and Wang, Dedong and Dubyagin, Stepan}, title = {Simulations of the inner magnetospheric energetic electrons using the IMPTAM-VERB coupled model}, series = {Journal of Atmospheric and Solar-Terrestrial Physics}, volume = {191}, journal = {Journal of Atmospheric and Solar-Terrestrial Physics}, publisher = {Elsevier}, address = {Oxford}, issn = {1364-6826}, doi = {10.1016/j.jastp.2019.05.014}, pages = {17}, year = {2019}, abstract = {In this study, we present initial results of the coupling between the Inner Magnetospheric Particle Transport and Acceleration Model (IMPTAM) and the Versatile Electron Radiation Belt (VERB-3D) code. IMPTAM traces electrons of 10-100 keV energies from the plasma sheet (L = 9 Re) to inner L-shell regions. The flux evolution modeled by IMPTAM is used at the low energy and outer L* computational boundaries of the VERB code (assuming a dipole approximation) to perform radiation belt simulations of energetic electrons. The model was tested on the March 17th, 2013 storm, for a six-day period. Four different simulations were performed and their results compared to satellites observations from Van Allen probes and GOES. The coupled IMPTAM-VERB model reproduces evolution and storm-time features of electron fluxes throughout the studied storm in agreement with the satellite data (within similar to 0.5 orders of magnitude). Including dynamics of the low energy population at L* = 6.6 increases fluxes closer to the heart of the belt and has a strong impact in the VERB simulations at all energies. However, inclusion of magnetopause losses leads to drastic flux decreases even below L* = 3. The dynamics of low energy electrons (max. 10s of keV) do not affect electron fluxes at energies >= 900 keV. Since the IMPTAM-VERB coupled model is only driven by solar wind parameters and the Dst and Kp indexes, it is suitable as a forecasting tool. In this study, we demonstrate that the estimation of electron dynamics with satellite-data-independent models is possible and very accurate.}, language = {en} } @article{WoodfieldGlauertMeniettietal.2019, author = {Woodfield, Emma E. and Glauert, Saraha A. and Menietti, J. Douglas and Averkamp, Terrance F. and Horne, Richard B. and Shprits, Yuri}, title = {Rapid Electron Acceleration in Low-Density Regions of Saturn's Radiation Belt by Whistler Mode Chorus Waves}, series = {Geophysical research letters}, volume = {46}, journal = {Geophysical research letters}, number = {13}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2019GL083071}, pages = {7191 -- 7198}, year = {2019}, abstract = {Electron acceleration at Saturn due to whistler mode chorus waves has previously been assumed to be ineffective; new data closer to the planet show it can be very rapid (factor of 104 flux increase at 1 MeV in 10 days compared to factor of 2). A full survey of chorus waves at Saturn is combined with an improved plasma density model to show that where the plasma frequency falls below the gyrofrequency additional strong resonances are observed favoring electron acceleration. This results in strong chorus acceleration between approximately 2.5 R-S and 5.5 R-S outside which adiabatic transport may dominate. Strong pitch angle dependence results in butterfly pitch angle distributions that flatten over a few days at 100s keV, tens of days at MeV energies which may explain observations of butterfly distributions of MeV electrons near L = 3. Including cross terms in the simulations increases the tendency toward butterfly distributions. Plain Language Summary Radiation belts are hazardous regions found around several of the planets in our Solar System. They consist of very hot, electrically charged particles trapped in the magnetic field of the planet. At Saturn the most important way to heat these particles has for many years been thought to involve the particles drifting closer toward the planet. This paper adds to the emerging idea at Saturn that a different way to heat the particles is also possible where the heating is done by waves, in a similar way to what we find at the Earth. We use recent information from the Cassini spacecraft on the number and location of particles and also of the waves strength and location combined with computer simulations to show that a particular wave called chorus is excellent at heating the particles where the surrounding number of cold particles is low.}, language = {en} } @article{DentonOfmanShpritsetal.2019, author = {Denton, Richard E. and Ofman, L. and Shprits, Yuri and Bortnik, J. and Millan, R. M. and Rodger, C. J. and da Silva, C. L. and Rogers, B. N. and Hudson, M. K. and Liu, K. and Min, K. and Glocer, A. and Komar, C.}, title = {Pitch Angle Scattering of Sub-MeV Relativistic Electrons by Electromagnetic Ion Cyclotron Waves}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {7}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9402}, doi = {10.1029/2018JA026384}, pages = {5610 -- 5626}, year = {2019}, abstract = {Electromagnetic ion cyclotron (EMIC) waves have long been considered to be a significant loss mechanism for relativistic electrons. This has most often been attributed to resonant interactions with the highest amplitude waves. But recent observations have suggested that the dominant energy of electrons precipitated to the atmosphere may often be relatively low, less than 1 MeV, whereas the minimum resonant energy of the highest amplitude waves is often greater than 2 MeV. Here we use relativistic electron test particle simulations in the wavefields of a hybrid code simulation of EMIC waves in dipole geometry in order to show that significant pitch angle scattering can occur due to interaction with low-amplitude short-wavelength EMIC waves. In the case we examined, these waves are in the H band (at frequencies above the He+ gyrofrequency), even though the highest amplitude waves were in the He band frequency range (below the He+ gyrofrequency). We also present wave power distributions for 29 EMIC simulations in straight magnetic field line geometry that show that the high wave number portion of the spectrum is in every case mostly due to the H band waves. Though He band waves are often associated with relativistic electron precipitation, it is possible that the He band waves do not directly scatter the sub-megaelectron volts (sub-MeV) electrons, but that the presence of He band waves is associated with high plasma density which lowers the minimum resonant energy so that these electrons can more easily resonate with the H band waves.}, language = {en} } @article{QinHudsonLietal.2019, author = {Qin, Murong and Hudson, Mary and Li, Zhao and Millan, Robyn and Shen, Xiaochen and Shprits, Yuri and Woodger, Leslie and Jaynes, Allison and Kletzing, Craig}, title = {Investigating loss of relativistic electrons associated with EMIC Waves at low L values on 22 June 2015}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {6}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA025726}, pages = {4022 -- 4036}, year = {2019}, abstract = {In this study, rapid loss of relativistic radiation belt electrons at low L* values (2.4-3.2) during a strong geomagnetic storm on 22 June 2015 is investigated along with five possible loss mechanisms. Both the particle and wave data are obtained from the Van Allen Probes. Duskside H+ band electromagnetic ion cyclotron (EMIC) waves were observed during a rapid decrease of relativistic electrons with energy above 5.2 MeV occurring outside the plasma sphere during extreme magnetopause compression. Lower He+ composition and enriched O+ composition are found compared to typical values assumed in other studies of cyclotron resonant scattering of relativistic electrons by EMIC waves. Quantitative analysis demonstrates that even with the existence of He+ band EMIC waves, it is the H+ band EMIC waves that are likely to cause the depletion at small pitch angles and strong gradients in pitch angle distributions of relativistic electrons with energy above 5.2 MeV at low L values for this event. Very low frequency wave activity at other magnetic local time can be favorable for the loss of relativistic electrons at higher pitch angles. An illustrative calculation that combines the nominal pitch angle scattering rate due to whistler mode chorus at high pitch angles with the H+ band EMIC wave loss rate at low pitch angles produces loss on time scale observed at L = 2.4-3.2. At high L values and lower energies, radial loss to the magnetopause is a viable explanation.}, language = {en} } @article{ShpritsVasileZhelayskaya2019, author = {Shprits, Yuri and Vasile, Ruggero and Zhelayskaya, Irina S.}, title = {Nowcasting and Predicting the Kp Index Using Historical Values and Real-Time Observations}, series = {Space Weather: The International Journal of Research and Applications}, volume = {17}, journal = {Space Weather: The International Journal of Research and Applications}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2018SW002141}, pages = {1219 -- 1229}, year = {2019}, abstract = {Current algorithms for the real-time prediction of the Kp index use a combination of models empirically driven by solar wind measurements at the L1 Lagrange point and historical values of the index. In this study, we explore the limitations of this approach, examining the forecast for short and long lead times using measurements at L1 and Kp time series as input to artificial neural networks. We explore the relative efficiency of the solar wind-based predictions, predictions based on recurrence, and predictions based on persistence. Our modeling results show that for short-term forecasts of approximately half a day, the addition of the historical values of Kp to the measured solar wind values provides a barely noticeable improvement. For a longer-term forecast of more than 2 days, predictions can be made using recurrence only, while solar wind measurements provide very little improvement for a forecast with long horizon times. We also examine predictions for disturbed and quiet geomagnetic activity conditions. Our results show that the paucity of historical measurements of the solar wind for high Kp results in a lower accuracy of predictions during disturbed conditions. Rebalancing of input data can help tailor the predictions for more disturbed conditions.}, language = {en} } @article{PostnovOskinovaTorrejon2017, author = {Postnov, K. and Oskinova, Lidia M. and Torrejon, J. M.}, title = {A propelling neutron star in the enigmatic Be-star gamma Cassiopeia}, series = {Monthly notices of the Royal Astronomical Society}, volume = {465}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnrasl/slw223}, pages = {L119 -- L123}, year = {2017}, abstract = {gamma Cassiopeia (gamma Cas), is known to be a binary system consisting of a Be-type star and a low-mass (M similar to 1M(circle dot)) companion of unknown nature orbiting in the Be-disc plane. Here, we apply the quasi-spherical accretion theory on to a compact magnetized star and show that if the low-mass companion of gamma Cas is a fast spinning neutron star, the key observational signatures of. Cas are remarkably well reproduced. Direct accretion on to this fast rotating neutron star is impeded by the propeller mechanism. In this case, around the neutron star magnetosphere a hot shell is formed which emits thermal X-rays in qualitative and quantitative agreement with observed properties of the X-ray emission from gamma Cas. We suggest that gamma Cas and its analogues constitute a new subclass of Be-type X-ray binaries hosting rapidly rotating neutron stars formed in supernova explosions with small kicks. The subsequent evolutionary stage of gamma Cas and its analogues should be the X Per-type binaries comprising low-luminosity slowly rotating X-ray pulsars. The model explains the enigmatic X-ray emission from gamma Cas, and also establishes evolutionary connections between various types of rotating magnetized neutron stars in Be-binaries.}, language = {en} } @article{MarkoetterSintschukBritzkeetal.2022, author = {Mark{\"o}tter, Henning and Sintschuk, Michael and Britzke, Ricardo and Dayani, Shahabeddin and Bruno, Giovanni}, title = {Upgraded imaging capabilities at the BAMline (BESSY II)}, series = {Journal of synchrotron radiation}, volume = {29}, journal = {Journal of synchrotron radiation}, number = {5}, publisher = {International Union of Crystallography}, address = {Chester}, issn = {1600-5775}, doi = {10.1107/S1600577522007342}, pages = {1292 -- 1298}, year = {2022}, abstract = {The BAMline at the BESSY II synchrotron X-ray source has enabled research for more than 20 years in widely spread research fields such as materials science, biology, cultural heritage and medicine. As a nondestructive characterization method, synchrotron X-ray imaging, especially tomography, plays a particularly important role in structural characterization. A recent upgrade of key equipment of the BAMline widens its imaging capabilities: shorter scan acquisition times are now possible, in situ and op erando studies can now be routinely performed, and different energy spectra can easily be set up. In fact, the upgraded double-multilayer monochromator brings full flexibility by yielding different energy spectra to optimize flux and energy resolution as desired. The upgraded detector (based on an sCMOS camera) also allows exploiting the higher flux with reduced readout times. Furthermore, an installed slip ring allows the sample stage to continuously rotate. The latter feature enables tomographic observation of processes occurring in the time scale of a few seconds.}, language = {en} } @article{PruefertBeitzReichetal.2022, author = {Pr{\"u}fert, Christian and Beitz, Toralf and Reich, Olaf and L{\"o}hmannsr{\"o}ben, Hans-Gerd}, title = {Inline process analysis of copper-bearing aerosols using laser-induced breakdown spectroscopy, laser-induced incandescence and optical imaging}, series = {Spectrochimica acta, Part B, Atomic spectroscopy}, volume = {197}, journal = {Spectrochimica acta, Part B, Atomic spectroscopy}, publisher = {Elsevier}, address = {Amsterdam [u.a.]}, issn = {0584-8547}, doi = {10.1016/j.sab.2022.106527}, pages = {11}, year = {2022}, abstract = {The quantification and identification of aerosols in industry plays a key role in process monitoring and control and lays the foundation for process automation aspired by the industry 4.0 initiative. However, measuring particulate matter's mass and number concentrations in harsh environments poses great analytical constraints. The presented approach comprises a comprehensive set of light-and imaging-based techniques, all contactless, in-line, and real-time. It includes, but is not limited to, stroboscopic imaging, laser-induced breakdown spectroscopy (LIBS) and laser-induced incandescence (LII). Stroboscopic imaging confirmed the particles sphericity and was used to measure the particle number density. A phase-selective LIBS setup with low fluence and 500 Hz repetition rate was used to classify each particle with a single-pulse and in real time. Simultaneously, the created plasma was captured by CCD imaging to determine the detection volume and hit rate of the LIBS setup. Both data sets combined were converted to a particle number density, which was consistent with the particle number density of the stroboscopic measurements. Furthermore, using a photodiode and microphone in parallel to the LIBS setup allowed for the photoacoustic normalization of the spectral line intensity at the laser repetition rate of 500 Hz. This was done as a partial photoacoustic normalization method with the cut-off based on the coefficient of variation (CV), reducing it by 25\%. Aside from that photodiode and microphone were proven to be valuable event counting with the advantage of the less spatially constricted. A second laser setup was used for laser -induced incandescence (LII) making it possible to classify the particles based on their incandescence tendency. Given its larger probing volume, LII could be employed at very low particle number densities. With respect to the current literature, this is the first approach of using LII as an in-line, real-time analytical technique for the compositional classification of metal-bearing aerosols.}, language = {en} } @article{ShenarSablowskiHainichetal.2019, author = {Shenar, Tomer and Sablowski, D. P. and Hainich, Rainer and Todt, Helge Tobias and Moffat, Anthony F. J. and Oskinova, Lidia M. and Ramachandran, Varsha and Sana, Hugues and Sander, Andreas Alexander Christoph and Schnurr, O. and St-Louis, N. and Vanbeveren, D. and Gotberg, Y. and Hamann, Wolf-Rainer}, title = {The Wolf-Rayet binaries of the nitrogen sequence in the Large Magellanic Cloud Spectroscopy, orbital analysis, formation, and evolution}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {627}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201935684}, pages = {68}, year = {2019}, abstract = {Context. Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core collapse. It is not known whether core He-burning WR stars (classical WR; cWR) form predominantly through wind stripping (w-WR) or binary stripping (b-WR). Whereas spectroscopy of WR binaries has so-far largely been avoided because of its complexity, our study focuses on the 44 WR binaries and binary candidates of the Large Magellanic Cloud (LMC; metallicity Z approximate to 0.5 Z(circle dot)), which were identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Aims. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at subsolar metallicity and constraining the impact of binary interaction in forming these stars. Methods. Spectroscopy was performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement was performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status was interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically homogeneous evolution. Results. Among our sample, 28/44 objects show composite spectra and are analyzed as such. An additional five targets show periodically moving WR primaries but no detected companions (SB1); two (BAT99 99 and 112) are potential WR + compact-object candidates owing to their high X-ray luminosities. We cannot confirm the binary nature of the remaining 11 candidates. About two-thirds of the WN components in binaries are identified as cWR, and one-third as hydrogen-burning WR stars. We establish metallicity-dependent mass-loss recipes, which broadly agree with those recently derived for single WN stars, and in which so-called WN3/O3 stars are clear outliers. We estimate that 45 +/- 30\% of the cWR stars in our sample have interacted with a companion via mass transfer. However, only approximate to 12 +/- 7\% of the cWR stars in our sample naively appear to have formed purely owing to stripping via a companion (12\% b-WR). Assuming that apparently single WR stars truly formed as single stars, this comprises approximate to 4\% of the whole LMC WN population, which is about ten times less than expected. No obvious differences in the properties of single and binary WN stars, whose luminosities extend down to log L approximate to 5.2 [L-circle dot], are apparent. With the exception of a few systems (BAT99 19, 49, and 103), the equatorial rotational velocities of the OB-type companions are moderate (v(eq) less than or similar to 250 km s(-1)) and challenge standard formalisms of angular-momentum accretion. For most objects, chemically homogeneous evolution can be rejected for the secondary, but not for the WR progenitor. Conclusions. No obvious dichotomy in the locations of apparently single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.}, language = {en} } @article{ShenarHainichTodtetal.2018, author = {Shenar, Tomer and Hainich, Rainer and Todt, Helge Tobias and Moffat, Anthony F. J. and Sander, Andreas Alexander Christoph and Oskinova, Lidia M. and Ramachandran, Varsha and Munoz, M. and Pablo, H. and Sana, Hugues and Hamann, Wolf-Rainer}, title = {The shortest-period Wolf-Rayet binary in the small magellanic cloud}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {616}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201833006}, pages = {15}, year = {2018}, abstract = {Context. SMC AB6 is the shortest-period (P = 6.5 d) Wolf-Rayet (WR) binary in the Small Magellanic Cloud. This binary is therefore a key system in the study of binary interaction and formation of WR stars at low metallicity. The WR component in AB6 was previously found to be very luminous (log L = 6.3 [L-circle dot]) compared to its reported orbital mass (approximate to 8 M-circle dot), placing it significantly above the Eddington limit. Aims. Through spectroscopy and orbital analysis of newly acquired optical data taken with the Ultraviolet and Visual Echelle Spectrograph (UVES), we aim to understand the peculiar results reported for this system and explore its evolutionary history. Methods. We measured radial velocities via cross-correlation and performed a spectral analysis using the Potsdam Wolf-Rayet model atmosphere code. The evolution of the system was analyzed using the Binary Population and Spectral Synthesis evolution code. Results. AB6 contains at least four stars. The 6.5 d period WR binary comprises the WR primary (WN3:h, star A) and a rather rapidly rotating (v(eq) = 265 km s(-1)) early O-type companion (O5.5 V, star B). Static N III and N IV emission lines and absorption signatures in He lines suggest the presence of an early-type emission line star (O5.5 I(f), star C). Finally, narrow absorption lines portraying a long-term radial velocity variation show the existence of a fourth star (O7.5 V, star D). Star D appears to form a second 140 d period binary together with a fifth stellar member, which is a B-type dwarf or a black hole. It is not clear that these additional components are bound to the WR binary. We derive a mass ratio of M-O/M-WR = 2.2 +/- 0.1. The WR star is found to be less luminous than previously thought (log L = 5.9 [L-circle dot]) and, adopting M-O = 41 M-circle dot for star B, more massive (M-WR = 18 M-circle dot). Correspondingly, the WR star does not exceed the Eddington limit. We derive the initial masses of M-i,M-WR = 60 M-circle dot and M-i,M-O = 40 M-circle dot and an age of 3.9 Myr for the system. The WR binary likely experienced nonconservative mass transfer in the past supported by the relatively rapid rotation of star B. Conclusions. Our study shows that AB6 is a multiple - probably quintuple - system. This finding resolves the previously reported puzzle of the WR primary exceeding the Eddington limit and suggests that the WR star exchanged mass with its companion in the past.}, language = {en} } @article{ShenarRichardsonSablowskietal.2017, author = {Shenar, Tomer and Richardson, N. D. and Sablowski, Daniel P. and Hainich, Rainer and Sana, H. and Moffat, A. F. J. and Todt, Helge Tobias and Hamann, Wolf-Rainer and Oskinova, Lidia M. and Sander, Andreas Alexander Christoph and Tramper, Frank and Langer, Norbert and Bonanos, Alceste Z. and de Mink, Selma E. and Gr{\"a}fener, G. and Crowther, Paul and Vink, J. S. and Almeida, Leonardo A. and de Koter, A. and Barb{\´a}, Rodolfo and Herrero, A. and Ulaczyk, Krzysztof}, title = {The tarantula massive binary monitoring}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {598}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201629621}, pages = {16}, year = {2017}, abstract = {We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159 d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300 M-circle dot, making it a candidate for being the most massive star known to date. While the primary is a known late-type, H-rich Wolf-Rayet star (WN6h), the secondary has so far not been unambiguously detected. Using moderate-resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0 : 78 and minimum masses of M-1 sin(3) i approximate to M-2 sin(3) i = 13 +/- 2 M-circle dot, with q = M-2/M-1 = 1.01 +/- 0.07. An analysis of emission excess stemming from a wind-wind collision yields an inclination similar to that obtained from polarimetry (i = 39 +/- 6 degrees). Our analysis thus implies M-1 = 53(-20)(+40) and M2 = 54(-20)(+40) M-circle dot, excluding M-1 > 300 M-circle dot. A detailed comparison with evolution tracks calculated for single and binary stars together with the high eccentricity suggests that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of approximate to 80 M-circle dot and initial masses of M-i,M-1 approximate to 10(5) and M-i,M-2 approximate to 90 M-circle dot, consistent with the upper limits of our derived orbital masses, and would imply an age of approximate to 2.2 Myr.}, language = {en} } @article{SanderFuerstKretschmaretal.2018, author = {Sander, Andreas Alexander Christoph and F{\"u}rst, F. and Kretschmar, P. and Oskinova, Lidia M. and Todt, Helge Tobias and Hainich, Rainer and Shenar, Tomer and Hamann, Wolf-Rainer}, title = {Coupling hydrodynamics with comoving frame radiative transfer}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {610}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201731575}, pages = {19}, year = {2018}, abstract = {Aims. To gain a realistic picture of the donor star in Vela X-1, we constructed a hydrodynamically consistent atmosphere model describing the wind stratification while properly reproducing the observed donor spectrum. To investigate how X-ray illumination affects the stellar wind, we calculated additional models for different X-ray luminosity regimes. Methods. We used the recently updated version of the Potsdam Wolf-Rayet code to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer. Results. The wind flow in Vela X-1 is driven by ions from various elements, with Fe III and S III leading in the outer wind. The model-predicted mass-loss rate is in line with earlier empirical studies. The mass-loss rate is almost unaffected by the presence of the accreting NS in the wind. The terminal wind velocity is confirmed at u(infinity) approximate to 600 km s(-1). On the other hand, the wind velocity in the inner region where the NS is located is only approximate to 100 km s(-1), which is not expected on the basis of a standard beta-velocity law. In models with an enhanced level of X-rays, the velocity field in the outer wind can be altered. If the X-ray flux is too high, the acceleration breaks down because the ionization increases. Conclusions. Accounting for radiation hydrodynamics, our Vela X-1 donor atmosphere model reveals a low wind speed at the NS location, and it provides quantitative information on wind driving in this important HMXB.}, language = {en} } @phdthesis{NovakovicMarinkovic2024, author = {Novakovic-Marinkovic, Nina}, title = {Optical control of bubble domains and skyrmions in thin films}, doi = {10.25932/publishup-64706}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-647069}, school = {Universit{\"a}t Potsdam}, pages = {ii, 106}, year = {2024}, abstract = {Laser induced switching offers an attractive possibility to manipulate small magnetic domains for prospective memory and logic devices on ultrashort time scales. Moreover, optical control of magnetization without high applied magnetic fields allows manipulation of magnetic domains individually and locally, without expensive heat dissipation. One of the major challenges for developing novel optically controlled magnetic memory and logic devices is reliable formation and annihilation of non-volatile magnetic domains that can serve as memory bits in ambient conditions. Magnetic skyrmions, topologically nontrivial spin textures, have been studied intensively since their discovery due to their stability and scalability in potential spintronic devices. However, skyrmion formation and, especially, annihilation processes are still not completely understood and further investigation on such mechanisms are needed. The aim of this thesis is to contribute to better understanding of the physical processes behind the optical control of magnetism in thin films, with the goal of optimizing material parameters and methods for their potential use in next generation memory and logic devices. First part of the thesis is dedicated to investigation of all-optical helicity-dependent switching (AO-HDS) as a method for magnetization manipulation. AO-HDS in Co/Pt multilayer and CoFeB alloys with and without the presence of Dzyaloshinskii-Moriya interaction (DMI), which is a type of exchange interaction, have been investigated by magnetic imaging using photo-emission electron microscopy (PEEM) in combination with X-ray magnetic circular dichroism (XMCD). The results show that in a narrow range of the laser fluence, circularly polarized laser light induces a drag on domain walls. This enables a local deterministic transformation of the magnetic domain pattern from stripes to bubbles in out-of-plane magnetized Co/Pt multilayers, only controlled by the helicity of ultrashort laser pulses. The temperature and characteristic fields at which the stripe-bubble transformation occurs has been calculated using theory for isolated magnetic bubbles, using as parameters experimentally determined average size of stripe domains and the magnetic layer thickness. The second part of the work aims at purely optical formation and annihilation of magnetic skyrmions by a single laser pulse. The presence of a skyrmion phase in the investigated CoFeB alloys was first confirmed using a Kerr microscope. Then the helicity-dependent skyrmion manipulation was studied using AO-HDS at different laser fluences. It was found that formation or annihilation individual skyrmions using AO-HDS is possible, but not always reliable, as fluctuations in the laser fluence or position can easily overwrite the helicity-dependent effect of AO-HDS. However, the experimental results and magnetic simulations showed that the threshold values for the laser fluence for the formation and annihilation of skyrmions are different. A higher fluence is required for skyrmion formation, and existing skyrmions can be annihilated by pulses with a slightly lower fluence. This provides a further option for controlling formation and annihilation of skyrmions using the laser fluence. Micromagnetic simulations provide additional insights into the formation and annihilation mechanism. The ability to manipulate the magnetic state of individual skyrmions is of fundamental importance for magnetic data storage technologies. Our results show for the first time that the optical formation and annihilation of skyrmions is possible without changing the external field. These results enable further investigations to optimise the magnetic layer to maximise the energy gap between the formation and annihilation barrier. As a result, unwanted switching due to small laser fluctuations can be avoided and fully deterministic optical switching can be achieved.}, language = {en} } @article{WoodfieldHorneGlauertetal.2018, author = {Woodfield, Emma E. and Horne, Richard B. and Glauert, S. A. and Menietti, J. D. and Shprits, Yuri and Kurth, William S.}, title = {Formation of electron radiation belts at Saturn by Z-mode wave acceleration}, series = {Nature Communications}, volume = {9}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-018-07549-4}, pages = {7}, year = {2018}, abstract = {At Saturn electrons are trapped in the planet's magnetic field and accelerated to relativistic energies to form the radiation belts, but how this dramatic increase in electron energy occurs is still unknown. Until now the mechanism of radial diffusion has been assumed but we show here that in-situ acceleration through wave particle interactions, which initial studies dismissed as ineffectual at Saturn, is in fact a vital part of the energetic particle dynamics there. We present evidence from numerical simulations based on Cassini spacecraft data that a particular plasma wave, known as Z-mode, accelerates electrons to MeV energies inside 4 RS (1 RS = 60,330 km) through a Doppler shifted cyclotron resonant interaction. Our results show that the Z-mode waves observed are not oblique as previously assumed and are much better accelerators than O-mode waves, resulting in an electron energy spectrum that closely approaches observed values without any transport effects included.}, language = {en} } @article{CaoNiSummersetal.2019, author = {Cao, Xing and Ni, Binbin and Summers, Danny and Shprits, Yuri and Gu, Xudong and Fu, Song and Lou, Yuequn and Zhang, Yang and Ma, Xin and Zhang, Wenxun and Huang, He and Yi, Juan}, title = {Sensitivity of EMIC wave-driven scattering loss of ring current protons to wave normal angle distribution}, series = {Geophysical research letters}, volume = {46}, journal = {Geophysical research letters}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2018GL081550}, pages = {590 -- 598}, year = {2019}, abstract = {Electromagnetic ion cyclotron waves have long been recognized to play a crucial role in the dynamic loss of ring current protons. While the field-aligned propagation approximation of electromagnetic ion cyclotron waves was widely used to quantify the scattering loss of ring current protons, in this study, we find that the wave normal distribution strongly affects the pitch angle scattering efficiency of protons. Increase of peak normal angle or angular width can considerably reduce the scattering rates of <= 10 keV protons. For >10 keV protons, the field-aligned propagation approximation results in a pronounced underestimate of the scattering of intermediate equatorial pitch angle protons and overestimates the scattering of high equatorial pitch angle protons by orders of magnitude. Our results suggest that the wave normal distribution of electromagnetic ion cyclotron waves plays an important role in the pitch angle evolution and scattering loss of ring current protons and should be incorporated in future global modeling of ring current dynamics.}, language = {en} } @article{MassaOskinovaPrinjaetal.2019, author = {Massa, Derck and Oskinova, Lidia M. and Prinja, Raman and Ignace, Richard}, title = {Coordinated UV and X-Ray Spectroscopic Observations of the O-type Giant xi Per}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {873}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/ab0283}, pages = {12}, year = {2019}, abstract = {We present new, contemporaneous Hubble Space Telescope STIS and XMM-Newton observations of the O7. III(n) ((f)) star xi Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of xi Per. The variable wind of this star is known to have a 2.086-day periodicity. We use a simple, heuristic spot model that fits the low-velocity (near-surface) IUE wind line variability very well, to demonstrate that the low-velocity absorption in the new STIS spectra of N IV lambda 1718 and Si IV lambda 1402 vary with the same 2.086-day period. It is remarkable that the period and amplitude of the STIS data agree with those of the IUE spectra obtained 22 yr earlier. We also show that the time variability of the new XMM-Newton fluxes is also consistent with the 2.086-day period. Thus, our new, multiwavelength coordinated observations demonstrate that the mechanism that causes the UV wind line variability is also responsible for a significant fraction of the X-rays in single O stars. The sequence of events for the multiwavelength light-curve minima is Si IV lambda 1402, N IV lambda 1718, and X-ray flux, each separated by a phase of about 0.06 relative to the 2.086-day period. Analysis of the X-ray fluxes shows that they become softer as they weaken. This is contrary to expectations if the variability is caused by periodic excess absorption. Furthermore, the high-resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blueshifted. If we interpret the low-velocity wind line light curves in terms of our model, it implies that there are two bright regions, i.e., regions with less absorption, separated by 180 degrees, on the surface of the star. We note that the presence and persistence of two spots separated by 180 degrees suggest that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in xi Per.}, language = {en} } @article{OskinovaFeldmeierHamann2006, author = {Oskinova, Lidia M. and Feldmeier, Achim and Hamann, Wolf-Rainer}, title = {High-resolution X-ray spectroscopy of bright O-type stars}, series = {Monthly notices of the Royal Astronomical Society}, volume = {372}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1111/j.1365-2966.2006.10858.x}, pages = {313 -- 326}, year = {2006}, abstract = {Archival X-ray spectra of the four prominent single, non-magnetic O stars zeta Pup, zeta Ori, xi Per and zeta Oph, obtained in high resolution with Chandra HETGS/MEG have been studied. The resolved X-ray emission line profiles provide information about the shocked, hot gas which emits the X-radiation, and about the bulk of comparably cool stellar wind material which partly absorbs this radiation. In this paper, we synthesize X-ray line profiles with a model of a clumpy stellar wind. We find that the geometrical shape of the wind inhomogeneities is important: better agreement with the observations can be achieved with radially compressed clumps than with spherical clumps. The parameters of the model, i.e. chemical abundances, stellar radius, mass-loss rate and terminal wind velocity, are taken from existing analyses of UV and optical spectra of the programme stars. On this basis, we also calculate the continuum-absorption coefficient of the cool-wind material, using the Potsdam Wolf-Rayet (POWR) model atmosphere code. The radial location of X-ray emitting gas is restricted from analysing the FIR line ratios of helium-like ions. The only remaining free parameter of our model is the typical distance between the clumps; here, we assume that at any point in the wind there is one clump passing by per one dynamical time-scale of the wind. The total emission in a model line is scaled to the observation. There is a good agreement between synthetic and observed line profiles. We conclude that the X-ray emission line profiles in O stars can be explained by hot plasma embedded in a cool wind which is highly clumped in the form of radially compressed shell fragments.}, language = {en} } @misc{KubatovaHamannKubatetal.2019, author = {Kubatova, Brankica and Hamann, Wolf-Rainer and Kubat, Jiri and Oskinova, Lidia M.}, title = {3D Monte Carlo Radiative Transfer in Inhomogeneous Massive Star Winds}, series = {Radiative signatures from the cosmos}, volume = {519}, journal = {Radiative signatures from the cosmos}, publisher = {Astronomical soc pacific}, address = {San Fransisco}, isbn = {978-1-58381-925-8}, issn = {1050-3390}, pages = {209 -- 212}, year = {2019}, abstract = {Already for decades it has been known that the winds of massive stars are inhomogeneous (i.e. clumped). To properly model observed spectra of massive star winds it is necessary to incorporate the 3-D nature of clumping into radiative transfer calculations. In this paper we present our full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. We use a set of parameters to describe dense as well as the rarefied wind components. At the same time, we account for non-monotonic velocity fields. We show how the 3-D density and velocity wind inhomogeneities strongly affect the resonance line formation. We also show how wind clumping can solve the discrepancy between P v and H alpha mass-loss rate diagnostics.}, language = {en} } @article{RobradeOskinovaSchmittetal.2018, author = {Robrade, Jan and Oskinova, Lidia M. and Schmitt, J. H. M. M. and Leto, Paolo and Trigilio, C.}, title = {Outstanding X-ray emission from the stellar radio pulsar CU Virginis}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {619}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201833492}, pages = {6}, year = {2018}, abstract = {Context. Among the intermediate-mass magnetic chemically peculiar (MCP) stars, CU Vir is one of the most intriguing objects. Its 100\% circularly polarized beams of radio emission sweep the Earth as the star rotates, thereby making this strongly magnetic star the prototype of a class of nondegenerate stellar radio pulsars. While CU Vir is well studied in radio, its high-energy properties are not known. Yet, X-ray emission is expected from stellar magnetospheres and confined stellar winds. Aims. Using X-ray data we aim to test CU Vir for intrinsic X-ray emission and investigate mechanisms responsible for its generation. Methods. We present X-ray observations performed with XMM-Newton and Chandra and study obtained X-ray images, light curves, and spectra. Basic X-ray properties are derived from spectral modelling and are compared with model predictions. In this context we investigate potential thermal and nonthermal X-ray emission scenarios. Results. We detect an X-ray source at the position of CU Vir. With LX approximate to 3 x 10(28) erg s(-1) it is moderately X-ray bright, but the spectrum is extremely hard compared to other Ap stars. Spectral modelling requires multi-component models with predominant hot plasma at temperatures of about T-X = 25MK or, alternatively, a nonthermal spectral component. Both types of model provide a virtually equivalent description of the X-ray spectra. The Chandra observation was performed six years later than those by XMM-Newton, yet the source has similar X-ray flux and spectrum, suggesting a steady and persistent X-ray emission. This is further confirmed by the X-ray light curves that show only mild X-ray variability. Conclusions. CU Vir is also an exceptional star at X-ray energies. To explain its full X-ray properties, a generating mechanism beyond standard explanations, like the presence of a low-mass companion or magnetically confined wind-shocks, is required. Magnetospheric activity might be present or, as proposed for fast-rotating strongly magnetic Bp stars, the X-ray emission of CU Vir is predominantly auroral in nature.}, language = {en} } @article{LetoTrigilioOskinovaetal.2018, author = {Leto, Paolo and Trigilio, C. and Oskinova, Lidia M. and Ignace, R. and Buemi, C. S. and Umana, G. and Ingallinera, A. and Leone, Francesco and Phillips, N. M. and Agliozzo, Claudia and Todt, Helge Tobias and Cerrigone, L.}, title = {A combined multiwavelength VLA/ALMA/Chandra study unveils the complex magnetosphere of the B-type star HR5907}, series = {Monthly notices of the Royal Astronomical Society}, volume = {476}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/sty244}, pages = {562 -- 579}, year = {2018}, abstract = {We present new radio/millimeter measurements of the hot magnetic star HR5907 obtained with the VLA and ALMA interferometers. We find that HR5907 is the most radio luminous early type star in the cm-mm band among those presently known. Its multi-wavelength radio light curves are strongly variable with an amplitude that increases with radio frequency. The radio emission can be explained by the populations of the non-thermal electrons accelerated in the current sheets on the outer border of the magnetosphere of this fast-rotating magnetic star. We classify HR5907 as another member of the growing class of strongly magnetic fast-rotating hot stars where the gyro-synchrotron emission mechanism efficiently operates in their magnetospheres. The new radio observations of HR5907 are combined with archival X-ray data to study the physical condition of its magnetosphere. The X-ray spectra of HR5907 show tentative evidence for the presence of non-thermal spectral component. We suggest that non-thermal X-rays originate a stellar X-ray aurora due to streams of non-thermal electrons impacting on the stellar surface. Taking advantage of the relation between the spectral indices of the X-ray power-law spectrum and the non-thermal electron energy distributions, we perform 3-D modelling of the radio emission for HR5907. The wavelength-dependent radio light curves probe magnetospheric layers at different heights above the stellar surface. A detailed comparison between simulated and observed radio light curves leads us to conclude that the stellar magnetic field of HR 5907 is likely non-dipolar, providing further indirect evidence of the complex magnetic field topology of HR5907.}, language = {en} } @article{OskinovaBulikGomezMoran2018, author = {Oskinova, Lidia M. and Bulik, Tomasz and Gomez-Moran, Ada Nebot}, title = {Infrared outbursts as potential tracers of common-envelope events in high-mass X-ray binary formation}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {613}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201832925}, pages = {7}, year = {2018}, abstract = {Context. Classic massive binary evolutionary scenarios predict that a transitional common-envelope (CE) phase could be preceded as well as succeeded by the evolutionary stage when a binary consists of a compact object and a massive star, that is, a high-mass X-ray binary (HMXB). The observational manifestations of common envelope are poorly constrained. We speculate that its ejection might be observed in some cases as a transient event at mid-infrared (IR) wavelengths. Aims. We estimate the expected numbers of CE ejection events and HMXBs per star formation unit rate, and compare these theoretical estimates with observations. Methods. We compiled a list of 85 mid-IR transients of uncertain nature detected by the Spitzer Infrared Intensive Transients Survey and searched for their associations with X-ray, optical, and UV sources. Results. Confirming our theoretical estimates, we find that only one potential HMXB may be plausibly associated with an IR-transient and tentatively propose that X-ray source NGC4490-X40 could be a precursor to the SPIRITS 16az event. Among other interesting sources, we suggest that the supernova remnant candidate [BWL2012] 063 might be associated with SPIRITS 16ajc. We also find that two SPIRITS events are likely associated with novae, and seven have potential optical counterparts. Conclusions. The massive binary evolutionary scenarios that involve CE events do not contradict currently available observations of IR transients and HMXBs in star-forming galaxies.}, language = {en} } @article{GonzalezGalanOskinovaPopovetal.2018, author = {Gonz{\´a}lez-Gal{\´a}n, Ana and Oskinova, Lidia M. and Popov, Sergei B. and Haberl, F. and K{\"u}hnel, M. and Gallagher, John S. and Schurch, Matthew and Guerrero, Mart{\´i}n A.}, title = {A multiwavelength study of SXP 1062, the long-period X-ray pulsar associated with a supernova remnant}, series = {Monthly notices of the Royal Astronomical Society}, volume = {475}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stx3127}, pages = {2809 -- 2821}, year = {2018}, abstract = {SXP 1062 is a Be X-ray binary (BeXB) located in the Small Magellanic Cloud. It hosts a long-period X-ray pulsar and is likely associated with the supernova remnant MCSNR J0127-7332. In this work we present a multiwavelength view on SXP 1062 in different luminosity regimes. We consider monitoring campaigns in optical (OGLE survey) and X-ray (Swift telescope). During these campaigns a tight coincidence of X-ray and optical outbursts is observed. We interpret this as typical Type I outbursts as often detected in BeXBs at periastron passage of the neutron star (NS). To study different X-ray luminosity regimes in depth, during the source quiescence we observed it with XMM-Newton while Chandra observations followed an X-ray outburst. Nearly simultaneously with Chandra observations in X-rays, in optical the RSS/SALT telescope obtained spectra of SXP 1062. On the basis of our multiwavelength campaign we propose a simple scenario where the disc of the Be star is observed face-on, while the orbit of the NS is inclined with respect to the disc. According to the model of quasi-spherical settling accretion our estimation of the magnetic field of the pulsar in SXP 1062 does not require an extremely strong magnetic field at the present time.}, language = {en} } @article{MartinezChicharroTorrejonOskinovaetal.2018, author = {Martinez-Chicharro, M. and Torrejon, J. M. and Oskinova, Lidia M. and Furst, F. and Postnov, K. and Rodes-Roca, J. J. and Hainich, Rainer and Bodaghee, A.}, title = {Evidence of Compton cooling during an X-ray flare supports a neutron star nature of the compact object in 4U1700-37}, series = {Monthly notices of the Royal Astronomical Society}, volume = {473}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnrasl/slx165}, pages = {L74 -- L78}, year = {2018}, abstract = {Based on new Chandra X-ray telescope data, we present empirical evidence of plasma Compton cooling during a flare in the non-pulsating massive X-ray binary 4U1700-37. This behaviour might be explained by quasi-spherical accretion on to a slowly rotating magnetized neutron star (NS). In quiescence, the NS in 4U1700-37 is surrounded by a hot radiatively cooling shell. Its presence is supported by the detection of mHz quasi-periodic oscillations likely produced by its convection cells. The high plasma temperature and the relatively low X-ray luminosity observed during the quiescence, point to a small emitting area similar to 1 km, compatible with a hotspot on an NS surface. The sudden transition from a radiative to a significantly more efficient Compton cooling regime triggers an episode of enhanced accretion resulting in a flare. During the flare, the plasma temperature drops quickly. The predicted luminosity for such transitions, similar to 3 x 10(35) erg s(-1), is very close to the luminosity of 4U1700-37 during quiescence. The transition may be caused by the accretion of a clump in the stellar wind of the donor star. Thus, a magnetized NS nature of the compact object is strongly favoured.}, language = {en} } @article{TorrejonReigFuerstetal.2018, author = {Torrejon, J. M. and Reig, Pablo and F{\"u}rst, F. and Martinez-Chicharro, M. and Postnov, K. and Oskinova, Lidia M.}, title = {NuSTAR rules out a cyclotron line in the accreting magnetar candidate 4U2206+54}, series = {Monthly notices of the Royal Astronomical Society}, volume = {479}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/sty1628}, pages = {3366 -- 3372}, year = {2018}, abstract = {Based on our new NuSTAR X-ray telescope data, we rule out any cyclotron line up to 60 keV in the spectra of the high-mass X-ray binary 4U2206+54. In particular, we do not find any evidence of the previously claimed line around 30 keV, independently of the source flux, along the spin pulse. The spin period has increased significantly, since the last observation, up to 5750 +/- 10 s, confirming the rapid spin-down rate (nu)over dot = -1.8 x 10(-14) Hz s(-1). This behaviour might be explained by the presence of a strongly magnetized neutron star (B-s > several times 10(13) G) accreting from the slow wind of its main-sequence O9.5 companion.}, language = {en} } @article{GrunerHainichSanderetal.2018, author = {Gruner, David and Hainich, Rainer and Sander, Andreas Alexander Christoph and Shenar, Tomer and Todt, Helge Tobias and Oskinova, Lidia M. and Ramachandran, Varsha and Ayres, T. and Hamann, Wolf-Rainer}, title = {The extreme O-type spectroscopic binary HD 93129A A quantitative, multiwavelength analysis}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {621}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201833178}, pages = {16}, year = {2018}, abstract = {Context. HD 93129A was classified as the earliest O-type star in the Galaxy (O2 If*) and is considered as the prototype of its spectral class. However, interferometry shows that this object is a binary system, while recent observations even suggest a triple configuration. None of the previous spectral analyses of this object accounted for its multiplicity. With new high-resolution UV and optical spectra, we have the possibility to reanalyze this key object, taking its binary nature into account for the first time. Aims. We aim to derive the fundamental parameters and the evolutionary status of HD 93129A, identifying the contributions of both components to the composite spectrum Results. Despite the similar spectral types of the two components, we are able to find signatures from each of the components in the combined spectrum, which allows us to estimate the parameters of both stars. We derive log(L/L-circle dot) = 6.15, T-eff = 52 kK, and log (M)over dot = -4.7[M-circle dot yr(-1)] for the primary Aa, and log(L/L-circle dot) = 5.58, T-eff = 45 kK, and log (M)over dot = -5.8 [M(circle dot)yr(-1)] for the secondary Ab. Conclusions. Even when accounting for the binary nature, the primary of HD 93129A is found to be one of the hottest and most luminous O stars in our Galaxy. Based on the theoretical decomposition of the spectra, we assign spectral types O2 If* and O3 III(f*) to components Aa and Ab, respectively. While we achieve a good fit for a wide spectral range, specific spectral features are not fully reproduced. The data are not sufficient to identify contributions from a hypothetical third component in the system.}, language = {en} } @article{HainichOskinovaShenaretal.2018, author = {Hainich, Rainer and Oskinova, Lidia M. and Shenar, Tomer and Marchant Campos, Pablo and Eldridge, J. J. and Sander, Andreas Alexander Christoph and Hamann, Wolf-Rainer and Langer, Norbert and Todt, Helge Tobias}, title = {Observational properties of massive black hole binary progenitors}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {609}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201731449}, pages = {62}, year = {2018}, abstract = {Context: The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ≈ 36 M⊙ and ≈ 29 M⊙. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models. Aims: Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate. Methods: Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (Potsdam Wolf-Rayet), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. Results: The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. Based on our atmosphere models, we provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Conclusions: While the predicted parameter space for massive BH binary progenitors is partly realized in nature, none of the known massive binaries match our synthetic spectra of massive BH binary progenitors exactly. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.}, language = {en} } @article{FulmerGallagherHamannetal.2020, author = {Fulmer, Leah M. and Gallagher, John S. and Hamann, Wolf-Rainer and Oskinova, Lidia M. and Ramachandran, Varsha}, title = {Testing massive star evolution, star-formation history, and feedback at low metallicity}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {633}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201834314}, pages = {9}, year = {2020}, abstract = {Context. The supergiant ionized shell SMC-SGS 1 (DEM 167), which is located in the outer Wing of the Small Magellanic Cloud (SMC), resembles structures that originate from an energetic star-formation event and later stimulate star formation as they expand into the ambient medium. However, stellar populations within and surrounding SMC-SGS 1 tell a different story. Aims. We present a photometric study of the stellar population encompassed by SMC-SGS 1 in order to trace the history of such a large structure and its potential influence on star formation within the low-density, low-metallicity environment of the SMC. Methods. For a stellar population that is physically associated with SMC-SGS 1, we combined near-ultraviolet (NUV) photometry from the Galaxy Evolution Explorer with archival optical (V-band) photometry from the ESO Danish 1.54 m Telescope. Given their colors and luminosities, we estimated stellar ages and masses by matching observed photometry to theoretical stellar isochrone models. Results. We find that the investigated region supports an active, extended star-formation event spanning similar to 25-40 Myr ago, as well as continued star formation into the present. Using a standard initial mass function, we infer a lower bound on the stellar mass from this period of similar to 3 x 10(4) M-circle dot, corresponding to a star-formation intensity of similar to 6 x 10(-3) M-circle dot kpc(-2) yr(-1). Conclusions. The spatial and temporal distributions of young stars encompassed by SMC-SGS 1 imply a slow, consistent progression of star formation over millions of years. Ongoing star formation, both along the edge and interior to SMC-SGS 1, suggests a combined stimulated and stochastic mode of star formation within the SMC Wing. We note that a slow expansion of the shell within this low-density environment may preserve molecular clouds within the volume of the shell, leaving them to form stars even after nearby stellar feedback expels local gas and dust.}, language = {en} } @article{ToalaBowmanVanReethetal.2022, author = {Toal{\´a}, Jes{\´u}s Alberto and Bowman, Dominic and Van Reeth, Timothy and Todt, Helge Tobias and Dsilva, Karan and Shenar, Tomer and Koenigsberger, Gloria Suzanne and Estrada-Dorado, Sandino and Oskinova, Lidia M. and Hamann, Wolf-Rainer}, title = {Multiple variability time-scales of the early nitrogen-rich Wolf-Rayet star WR 7}, series = {Monthly notices of the Royal Astronomical Society}, volume = {514}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stac1455}, pages = {2269 -- 2277}, year = {2022}, abstract = {We present the analysis of the optical variability of the early, nitrogen-rich Wolf-Rayet (WR) star WR 7. The analysis of multisector Transiting Exoplanet Survey Satellite (TESS) light curves and high-resolution spectroscopic observations confirm multiperiodic variability that is modulated on time-scales of years. We detect a dominant period of 2.6433 +/- 0.0005 d in the TESS sectors 33 and 34 light curves in addition to the previously reported high-frequency features from sector 7. We discuss the plausible mechanisms that may be responsible for such variability in WR 7, including pulsations, binarity, co-rotating interaction regions (CIRs), and clumpy winds. Given the lack of strong evidence for the presence of a stellar or compact companion, we suggest that WR 7 may pulsate in quasi-coherent modes in addition to wind variability likely caused by CIRs on top of stochastic low-frequency variability. WR 7 is certainly a worthy target for future monitoring in both spectroscopy and photometry to sample both the short (less than or similar to 1 d) and long (greater than or similar to 1000 d) variability time-scales.}, language = {en} } @article{OskinovaSchaerer2022, author = {Oskinova, Lidia M. and Schaerer, Daniel}, title = {Ionization of He II in star-forming galaxies by X-rays from cluster winds and superbubbles}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {661}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/202142520}, pages = {6}, year = {2022}, abstract = {The nature of the sources powering nebular He II emission in star-forming galaxies remains debated, and various types of objects have been considered, including Wolf-Rayet stars, X-ray binaries, and Population III stars. Modern X-ray observations show the ubiquitous presence of hot gas filling star-forming galaxies. We use a collisional ionization plasma code to compute the specific He II ionizing flux produced by hot gas and show that if its temperature is not too high (less than or similar to 2.5 MK), then the observed levels of soft diffuse X-ray radiation could explain He II ionization in galaxies. To gain a physical understanding of this result, we propose a model that combines the hydrodynamics of cluster winds and hot superbubbles with observed populations of young massive clusters in galaxies. We find that in low-metallicity galaxies, the temperature of hot gas is lower and the production rate of He II ionizing photons is higher compared to high-metallicity galaxies. The reason is that the slower stellar winds of massive stars in lower-metallicity galaxies input less mechanical energy in the ambient medium. Furthermore, we show that ensembles of star clusters up to similar to 10-20 Myr old in galaxies can produce enough soft X-rays to induce nebular He II emission. We discuss observations of the template low-metallicity galaxy I Zw 18 and suggest that the He II nebula in this galaxy is powered by a hot superbubble. Finally, appreciating the complex nature of stellar feedback, we suggest that soft X-rays from hot superbubbles are among the dominant sources of He II ionizing flux in low-metallicity star-forming galaxies.}, language = {en} }