A revision of soft proton scattering at grazing incidence and its implementation in the <Emphasis Type="SmallCaps">geant</Emphasis>4 toolkit

The scattering of soft protons inside the Wolter-type optics of X-ray observatories has been proven to concentrate these particles onto the focal plane instruments. The funneling of these protons increases the instrumental background and can also contribute to the degrading of the detectors. The instrumental background and degradation of the detector’s performance experienced by Chandra and XMM-Newton is significantly larger than what was expected on the basis of previous Monte Carlo simulations. For Chandra the main issue is the degradation of the energy resolution due to lattice displacements in the detectors. For XMM the contribution to the instrumental background is more significant. In between, new laboratory measurements as well as a revision of the theory are needed to correctly assess the impact of the environmental radiation for future missions. In this publication we present a Geant4 class that will allow future users to select between either theoretical models or measured data to simulate the scattering of soft protons at grazing angles. To develop this method, we revisit the theory of elastic scattering of protons on polished surfaces and implement these approaches into Geant4. We also implemented recently performed measurements using parts of eROSITA (extended ROentgen Survey with an Imaging Telescope Array) mirror shells as scattering targets as another scattering model to be used within the Geant4 toolkit.     

The spin state of 433 Eros and its possible implications

In this paper, we show that Asteroid (433) Eros is currently residing in a spin-orbit resonance, with its spin axis undergoing a small-amplitude libration about the Cassini state 2 of the proper mode in the nonsingular orbital element sinI/2exp(?Ω), where I the orbital inclination and Ω the longitude of the node. The period of this libration is ?53.4 kyr. By excluding these libration wiggles, we find that Eros' pole precesses with the proper orbital plane in inertial space with a period of ?61.4 kyr. Eros' resonant state forces its obliquity to oscillate with a period of ?53.4 kyr between ?76° and ?89.5°. The observed value of ?89° places it near the latter extreme of this cycle. We have used these results to probe Eros' past orbit and spin evolution. Our computations suggest that Eros is unlikely to have achieved its current spin state by solar and planetary gravitational perturbations alone. We hypothesize that some dissipative process such as thermal torques (e.g., the so-called YORP effect) may be needed in our model to obtain a more satisfactory match with data. A detailed study of this problem is left for future work.     

Elemental Abundances in the Solar Corona as Measured by the X-ray Solar Monitor Onboard Chandrayaan-1

The X-ray Solar Monitor (XSM) on the Indian lunar mission Chandrayaan-1 was flown to complement lunar elemental abundance studies by the X-ray fluorescence experiment C1XS. XSM measured the ≈?1.8?–?20 keV solar X-ray spectrum during its nine months of operation in lunar orbit. The soft X-ray spectra can be used to estimate absolute coronal abundances using intensities of emission-line complexes and the plasma temperature derived from the continuum. The best estimates are obtained from the brightest flare observed by XSM: a C2.8-class flare. The well-known first-ionization potential (FIP) effect is observed; abundances are enhanced for the low-FIP elements Fe, Ca, and Si, while the intermediate-FIP element S shows values close to the photospheric abundance. The derived coronal abundances show a quasi-mass-dependent pattern of fractionation.     

SOVAP/Picard,a Spaceborne Radiometer to Measure the Total Solar Irradiance

The Picard spacecraft was successfully launched on 15 June 2010, into a Sun-synchronous orbit. The mission represents one of the European contributions to solar observations and Essential Climate Variables (ECVs) measurements. The payload is composed of a Solar Diameter Imager and Surface Mapper (SODISM) and two radiometers: SOlar VAriability Picard (SOVAP) and PREcision MOnitor Sensor (PREMOS). SOVAP, a dual side-by-side cavity radiometer, measures the total solar irradiance (TSI). It is the sixth of a series of differential absolute-radiometer-type instruments developed and operated in space by the Royal Meteorological Institute of Belgium. The measurements of SOVAP in the summer of 2010 yielded a TSI value of 1362.1 W?m^?2 with an uncertainty of ±?2.4 W?m^?2 (k=1). During the periods of November 2010 and January 2013, the amplitude of the changes in TSI has been on the order of 0.18 %, corresponding to a range of about 2.4 W?m^?2.     

Resolving Differences in Absolute Irradiance Measurements Between the SOHO/CELIAS/SEM and the SDO/EVE

The Solar EUV Monitor (SEM) onboard SOHO has measured absolute extreme ultraviolet (EUV) and soft X-ray solar irradiance nearly continuously since January 1996. The EUV Variability Experiment (EVE) on SDO, in operation since April of 2010, measures solar irradiance in a wide spectral range that encompasses the band passes (26?–?34 nm and 0.1?–?50 nm) measured by SOHO/SEM. However, throughout the mission overlap, irradiance values from these two instruments have differed by more than the combined stated uncertainties of the measurements. In an effort to identify the sources of these differences and eliminate them, we investigate in this work the effect of reprocessing the SEM data using a more accurate SEM response function (obtained from synchrotron measurements with a SEM sounding-rocket clone instrument taken after SOHO was already in orbit) and time-dependent, measured solar spectral distributions – i.e., solar reference spectra that were unavailable prior to the launch of the SDO. We find that recalculating the SEM data with these improved parameters reduces mean differences with the EVE measurements from about 20 % to less than 5 % in the 26?–?34 nm band, and from about 35 % to about 15 % for irradiances in the 0.1?–?7 nm band extracted from the SEM 0.1?–?50 nm channel.     

The LYRA Instrument Onboard PROBA2: Description and In-Flight Performance

The Large Yield Radiometer (LYRA) is an XUV–EUV–MUV (soft X-ray to mid-ultraviolet) solar radiometer onboard the European Space Agency Project for On-Board Autonomy 2 (PROBA2) mission, which was launched in November 2009. LYRA acquires solar-irradiance measurements at a high cadence (nominally 20?Hz) in four broad spectral channels, from soft X-ray to MUV, which have been chosen for their relevance to solar physics, space weather, and aeronomy. We briefly review the design of the instrument, give an overview of the data products distributed through the instrument website, and describe how the data are calibrated. We also briefly present a summary of the main fields of research currently under investigation by the LYRA consortium.     

Relation between the X-ray and optical luminosities in binary systems with accreting nonmagnetic white dwarfs

We investigate the relation between the optical (g-band) and X-ray (0.5–10 keV) luminosities of accreting nonmagnetic white dwarfs. According to the present-day counts of the populations of star systems in our Galaxy, these systems have the highest space density among the close binary systems with white dwarfs. We show that the dependence of the optical luminosity of accreting white dwarfs on their X-ray luminosity forms a fairly narrow one-parameter curve. The typical half-width of this curve does not exceed 0.2–0.3 dex in optical and X-ray luminosities, which is essentially consistent with the amplitude of the aperiodic flux variability for these objects. At X-ray luminosities L _x ～ 10³² erg s^?1 or lower, the optical g-band luminosity of the accretion flow is shown to be related to its X-ray luminosity by a factor ～2–3. At even lower X-ray luminosities (L _x ? 10³⁰ erg s^?1), the contribution from the photosphere of the white dwarf begins to dominate in the optical spectrum of the binary system and its optical brightness does not drop below M _g ～ 13–14. Using the latter fact, we show that in current and planned X-ray sky surveys, the family of accreting nonmagnetic white dwarfs can be completely identified to the distance determined by the sensitivity of an optical sky survey in this region. For the Sloan Digital Sky Survey (SDSS) with a limiting sensitivity m _g ～ 22.5, this distance is ～400–600 pc.     

First Optical Satellite Tracking Station (OSTS) at NRIAG-Egypt

In March 2019 the first dedicated of optical observation of space debris and artificial satellites (Optical Satellites Tracking Station (OSTS)) in Egypt has been performed by National Research Institute of Astronomy and Geophysics (NRIAG) at Kottamia Observatory. The 0.28 m Telescope is used for tracking and surveying debris and operational satellites at Low Earth Orbit (LEO), High ellipse orbit (HEO) and the Geosynchronous Earth Orbit (GEO). OSTS has also collaborated with International Scientific Optical Network (ISON) for optical observation. The necessary programs to control the telescope, camera and monitoring process are configured. Coordinate corrected metric data is provided with the time information. The system has validated and calibrated processing. The first results of the observations with image processing using Apex-2 software are presented. The optical observations using OSTS are being used to help characterize the debris environment in LEO, HEO and GEO to assist in the modeling projections for space debris database, real population and distribution, detection and orbital determination, and conjunction analysis between operational satellites and/or dangers debris.     

ChemiX: a Bragg crystal spectrometer for the <Emphasis Type="Italic">Interhelioprobe</Emphasis> interplanetary mission

Interhelioprobe (IHP), an analogue to the ESA Solar Orbiter, is the prospective Russian space solar observatory intended for in-situ and remote sensing investigations of the Sun and the inner heliosphere from a heliocentric orbit with the perihelion of about 60 solar radii. One of several instruments on board will be the Bragg crystal spectrometer ChemiX which will measure X-ray spectra from solar corona structures. Analysis of the spectra will allow the determination of the elemental composition of plasma in hot coronal sources like flares and active regions. ChemiX is under development at the Wroc?aw Solar Physics Division of the Polish Academy of Sciences Space Research Centre in collaboration with an international team (see the co-author list). This paper gives an overview of the ChemiX scientific goals and design preparatory to phase B of the instrument development.     

X-ray Luminosity Function of Quasars at 3 &lt; <Emphasis Type="Italic">z</Emphasis> &lt; 5 from XMM-Newton Serendipitous Survey Data

The X-ray luminosity function of distant (3 < z < 5.1) type 1 quasars has been measured. A sample of distant high-luminosity (10⁴⁵ erg s^?1 ≤ L_X,2?10 < 7.5×10⁴⁵ erg s^?1 in the 2–10 keV energy band) quasars from the catalog by Khorunzhev et al. (2016) compiled from the data of the 3XMM-DR4 catalog of the XMM-Newton serendipitous survey and the Sloan Digital Sky Survey (SDSS) has been used. This sample consists of 101 sources. Most of them (90) have spectroscopic redshifts z_spec ? 3; the remaining ones are quasar candidates with photometric redshift estimates z_phot ? 3. The spectroscopic redshifts of eight sources have been measured with the BTA and AZT-33IK telescopes. Owing to the record sky coverage area (?250 sq. deg at X-ray fluxes ~10^?14 erg s^?1 cm^?2 in the 0.5–2 keVband) from which the sample was drawn, we have managed to obtain reliable estimates of the space density of distant X-ray quasars with luminosities L_X,2?10 > 2×10⁴⁵ erg s^?1 for the first time. Their comoving space density remains constant as the redshift increases from z = 3 to 5 to within a factor of 2. The power-law slope of the X-ray luminosity function of distant quasars at its bright end (above the break) has been reliably constrained for the first time. The range of possible slopes for the quasar luminosity and density evolution model is γ₂ = 2.72 _?0.12 ^+0.19 ± 0.21, where initially the lower and upper boundaries of γ₂ with the remaining uncertainty in the detection completeness of X-ray sources in SDSS and subsequently the statistical error of the slope are specified.     

Assessing the Accuracy of CME Speed and Trajectory Estimates from STEREO Observations Through a Comparison of Independent Methods

We have estimated the speed and direction of propagation of a number of Coronal Mass Ejections (CMEs) using single-spacecraft data from the STEREO Heliospheric Imager (HI) wide-field cameras. In general, these values are in good agreement with those predicted by Thernisien, Vourlidas, and Howard in Solar Phys. 256, 111?–?130 (2009) using a forward modelling method to fit CMEs imaged by the STEREO COR2 coronagraphs. The directions of the CMEs predicted by both techniques are in good agreement despite the fact that many of the CMEs under study travel in directions that cause them to fade rapidly in the HI images. The velocities estimated from both techniques are in general agreement although there are some interesting differences that may provide evidence for the influence of the ambient solar wind on the speed of CMEs. The majority of CMEs with a velocity estimated to be below 400 km?s^?1 in the COR2 field of view have higher estimated velocities in the HI field of view, while, conversely, those with COR2 velocities estimated to be above 400 km?s^?1 have lower estimated HI velocities. We interpret this as evidence for the deceleration of fast CMEs and the acceleration of slower CMEs by interaction with the ambient solar wind beyond the COR2 field of view. We also show that the uncertainties in our derived parameters are influenced by the range of elongations over which each CME can be tracked. In order to reduce the uncertainty in the predicted arrival time of a CME at 1 Astronomical Unit (AU) to within six hours, the CME needs to be tracked out to at least 30 degrees elongation. This is in good agreement with predictions of the accuracy of our technique based on Monte Carlo simulations. Within the set of studied CMEs, there are two clear events that were predicted from the HI data to travel over another spacecraft; in-situ measurements at these other spacecraft confirm the accuracy of these predictions. The ability of the HI cameras to image Corotating Interaction Region (CIR)-entrained transients as well as CMEs can result in some ambiguity when trying to distinguishing individual signatures.     

The State of Self-organized Criticality of the Sun During the Last Three Solar Cycles. I. Observations

We analyze the occurrence-frequency distributions of peak fluxes [P], total fluxes [E], and durations [T] of solar flares over the last three solar cycles (during 1980??C?2010) from SMM/HXRBS, CGRO/BATSE, and RHESSI hard X-ray data. From the synthesized data we find powerlaw slopes with mean values of ?? _P =1.73±0.07 for the peak flux, ?? _E =1.62±0.12 for the total flux, and ?? _T =1.99±0.35 for flare durations. We find a tendency of an anti-correlation of the powerlaw slope of peak fluxes with the flare rate or sunspot number as a function of the solar cycle. The occurrence powerlaw slope is always steeper by ??????0.1 during a solar-cycle minimum compared with the previous solar-cycle maximum, but the relative amplitude varies for each cycle or instrument. Since each solar cycle has been observed with a different instrument, part of the variation could be attributed to instrumental characteristics and different event selection criteria used in generating the event catalogs. The relatively flatter powerlaw slopes during solar maxima could indicate more energetic flares with harder electron-energy spectra, probably due to a higher magnetic complexity of the solar corona. This would imply a non-stationarity (or solar-cycle dependence) of the coronal state of self-organized criticality.     

Obliquity evolution of extrasolar terrestrial planets

We have investigated the obliquity evolution of terrestrial planets in habitable zones (at ∼1 AU) in extrasolar planetary systems, due to tidal interactions with their satellite and host star with wide varieties of satellite-to-planet mass ratio (m/Mp) and initial obliquity (γ₀), through numerical calculations and analytical arguments. The obliquity, the angle between planetary spin axis and its orbit normal, of a terrestrial planet is one of the key factors in determining the planetary surface environments. A recent scenario of terrestrial planet accretion implies that giant impacts of Mars-sized or larger bodies determine the planetary spin and form satellites. Since the giant impacts would be isotropic, tilted spins (sinγ₀∼1) are more likely to be produced than straight ones (sinγ₀∼0). The ratio m/Mp is dependent on the impact parameters and impactors' mass. However, most of previous studies on tidal evolution of the planet-satellite systems have focused on a particular case of the Earth-Moon systems in which m/Mp?0.0125 and γ₀∼10° or the two-body planar problem in which γ₀=0° and stellar torque is neglected. We numerically integrated the evolution of planetary spin and a satellite orbit with various m/Mp (from 0.0025 to 0.05) and γ₀ (from 0° to 180°), taking into account the stellar torques and precessional motions of the spin and the orbit. We start with the spin axis that almost coincides with the satellite orbit normal, assuming that the spin and the satellite are formed by one dominant impact. With initially straight spins, the evolution is similar to that of the Earth-Moon system. The satellite monotonically recedes from the planet until synchronous state between the spin period and the satellite orbital period is realized. The obliquity gradually increases initially but it starts decreasing down to zero as approaching the synchronous state. However, we have found that the evolution with initially tiled spins is completely different. The satellite's orbit migrates outward with almost constant obliquity until the orbit reaches the critical radius ∼10-20 planetary radii, but then the migration is reversed to inward one. At the reversal, the obliquity starts oscillation with large amplitude. The oscillation gradually ceases and the obliquity is reduced to ∼0° during the inward migration. The satellite eventually falls onto the planetary surface or it is captured at the synchronous state at several planetary radii. We found that the character change of precession about total angular momentum vector into that about the planetary orbit normal is responsible for the oscillation with large amplitude and the reversal of migration. With the results of numerical integration and analytical arguments, we divided the m/Mp-γ₀ space into the regions of the qualitatively different evolution. The peculiar tidal evolution with initially tiled spins give deep insights into dynamics of extrasolar planet-satellite systems and discussions of surface environments of the planets.     

Experimental study of parameters of X-ray radiation from solar flares using the PENGUIN-M instrument aboard the CORONAS-PHOTON spacecraft

The main characteristics of the PENGUIN-M instrument are given. The instrument has been operating aboard the CORONAS-PHOTON spacecraft (SC) launched into orbit on January 30, 2009. The instrument includes the PENGUIN-MD detector unit (PMD) and the PENGUIN-ME electronic unit (PMD). The purpose of the experiment is to measure the degree of linear polarization of X-ray radiation from solar flares in the energy range of 20–150 keV and to obtain energy spectra of X-ray radiation from solar flares in the energy range of 2–500 keV. The paper describes the instrument, calibration procedure, and in-flight adjustment, and contains the first results of measurements.     

Specifics of the orbital evolution of the CORONAS-F satellite at the final stage of its flight

We report the results of our analysis of the variations in the orbit of the CORONAS-F satellite (August 31, 2001–December 6, 2005) during the low-orbit flight stage. We show that extreme solar events and the accompanying phenomena in the near-Earth space and in Earth’s upper atmosphere during the 2004–2005 period (the declining phase of 23rd solar-activity cycle) significantly affected the duration of the “low-orbit” flight stage of the CORONAS-F satellite.     

A multi-objective approach to the design of low thrust space trajectories using optimal control

In this article, we introduce a novel three-step approach for solving optimal control problems in space mission design. We demonstrate its potential by the example task of sending a group of spacecraft to a specific Earth L ₂ halo orbit. In each of the three steps we make use of recently developed optimization methods and the result of one step serves as input data for the subsequent one. Firstly, we perform a global and multi-objective optimization on a restricted class of control functions. The solutions of this problem are (Pareto-)optimal with respect to ΔV and flight time. Based on the solution set, a compromise trajectory can be chosen suited to the mission goals. In the second step, this selected trajectory serves as initial guess for a direct local optimization. We construct a trajectory using a more flexible control law and, hence, the obtained solutions are improved with respect to control effort. Finally, we consider the improved result as a reference trajectory for a formation flight task and compute trajectories for several spacecraft such that these arrive at the halo orbit in a prescribed relative configuration. The strong points of our three-step approach are that the challenging design of good initial guesses is handled numerically by the global optimization tool and afterwards, the last two steps only have to be performed for one reference trajectory.     

Orbital and physical characteristics of micrometeoroids in the inner solar system as observed by Helios 1

The Helios 1 spacecraft was launched in December 1974 into a heliocentric orbit of 0.3 AU perihelion distance. Helios 2 followed one year later on a similar orbit. Both spaceprobes carry on board micrometeoroid experiments each of which contains two sensors with a total sensitive area of 121 cm². To date, only preliminary data are available from Helios 2. Therefore the results presented here mainly apply to data from Helios 1. The ecliptic sensor of Helios 1 measures dust particles which have trajectories with elevations from ?45° to + 55° with respect to the ecliptic plane. The south sensor detects dust particles with trajectory elevations from ?90° (ecliptic south-pole) to ?4°. The ecliptic sensor is covered by a thin film (3000 Å parylene coated with 750 Å aluminium) as protection against solar radiation. The other sensor is shielded by the spacecraft rim from direct sunlight and has an open aperture. Micrometeoroids are detected by the electric charge produced upon impact. During the first 6 orbits of Helios 1 around the sun the experiment registered a total of 168 meteoroids, 52 particles were detected by the ecliptic sensor and 116 particles by the south sensor. This excess of impacts on the south sensor with regard to the impacts on the ecliptic sensor is due predominantly to small impacts which are characterized by small pulse heights of the charge signals. But also large impacts were statistically significantly more abundant on the south sensor than on the ecliptic sensor. Most impacts on the ecliptic sensor were observed when it was pointing in the direction of motion of Helios (apex direction). In contrast to that the south sensor detected most impacts when it was facing in between the solar and antapex direction. Orbit analysis showed that the “apex” particles which are predominantly detected by the ecliptic sensor have eccentricities e < 0.4 or semi-major axes a ? 0.5 AU. From a comparison with corresponding data from the south sensor it is concluded that the average inclination f of “apex” particles is -i < 30°. The excess of impacts on the south sensor, called “eccentric” particles, have orbit eccentricities e > 0.4 and semimajor axes a > 0.5AU. β-meteoroids leaving the solar system on hyperbolic orbits are directly identified by the observed imbalance of outgoing (away from the sun) and ingoing particles. It is shown that “eccentric” particles, due to their orbital characteristics, should be observable also by the ecliptic sensor. Since they have not been detected by this sensor it is concluded that the only instrumental difference between both sensors, i.e. the entrance film in front of the ecliptic sensor, prevented them from entering it. A comparison with penetration studies proved that particles which do not penetrate the entrance film must have bulk densities ρ(g/cm³) below an upper density limit ρ_max. It is shown that approximately 30% of the “eccentric” particles have densities below ρ_max = 1 g/cm³.     

The dust trail of Comet 67P/Churyumov-Gerasimenko

We report the detection of Comet 67P/Churyumov-Gerasimenko's dust trail and nucleus in 24 μm Spitzer Space Telescope images taken February 2004. The dust trail is not found in optical Palomar images taken June 2003. Both the optical and infrared images show a distinct neck-line tail structure, offset from the projected orbit of the comet. We compare our observations to simulated images using a Monte Carlo approach and a dynamical model for comet dust. We estimate the trail to be at least one orbit old (6.6 years) and consist of particles of size ?100 μm. The neck-line is composed of similar sized particles, but younger in age. Together, our observations and simulations suggest grains 100 μm and larger in size dominate the total mass ejected from the comet. The radiometric effective radius of the nucleus is 1.87±0.08 km, derived from the Spitzer observation. The Rosetta spacecraft is expected to arrive at and orbit this comet in 2014. Assuming the trail is comprised solely of 1 mm radius grains, we compute a low probability (∼¹⁰⁻³) of a trail grain impacting with Rosetta during approach and orbit insertion.     

Observational Tracking of the 2D Structure of Coronal Mass Ejections Between the Sun and 1 AU

The Solar TErrestrial RElations Observatory (STEREO) provides high cadence and high resolution images of the structure and morphology of coronal mass ejections (CMEs) in the inner heliosphere. CME directions and propagation speeds have often been estimated through the use of time-elongation maps obtained from the STEREO Heliospheric Imager (HI) data. Many of these CMEs have been identified by citizen scientists working within the SolarStormWatch project ( www.solarstormwatch.com ) as they work towards providing robust real-time identification of Earth-directed CMEs. The wide field of view of HI allows scientists to directly observe the two-dimensional (2D) structures, while the relative simplicity of time-elongation analysis means that it can be easily applied to many such events, thereby enabling a much deeper understanding of how CMEs evolve between the Sun and the Earth. For events with certain orientations, both the rear and front edges of the CME can be monitored at varying heliocentric distances (R) between the Sun and 1?AU. Here we take four example events with measurable position angle widths and identified by the citizen scientists. These events were chosen for the clarity of their structure within the HI cameras and their long track lengths in the time-elongation maps. We show a linear dependency with R for the growth of the radial width (W) and the 2D aspect ratio (??) of these CMEs, which are measured out to ???0.7?AU. We estimated the radial width from a linear best fit for the average of the four CMEs. We obtained the relationships W=0.14R+0.04 for the width and ??=2.5R+0.86 for the aspect ratio (W and R in units of?AU).