Dark energy,co-evolution of massive black holes with galaxies,and ASTROD-GW

The detection of low frequency band (100 nHz–100 mHz) and very low frequency band (300 pHz–100 nHz) gravitational waves (GWs) is important for exploration of the equation of state of dark energy and the co-evolution of massive black holes (MBHs) with galaxies. Most galaxies are believed to have a massive black hole in the galactic core. In the formation of these black holes, merging and accretion are the two main processes. Merging of massive black holes generate GWs which could be detected by space GW detectors and Pulsar Timing Arrays (PTAs) to cosmological distances. LISA (Laser-Interferometric Space Antenna) is most sensitive to the frequency band 1 mHz–100 mHz, ASTROD-GW (ASTROD [Astrodynamical Space Test of Relativity using Optical Devices] optimized for Gravitational Wave detection) is most sensitive to the frequency band 100 nHz–1 mHz and PTAs are most sensitive to the frequency band 300 pHz–100 nHz. In this paper, we discuss the sensitivities and outlooks of detection of GWs from binary massive black holes in these frequency bands with an emphasis on ASTROD-GW. The GWs generated by the inspirals, merging and subsequent ringdowns of binary black holes are standard sirens to the cosmological distance. Using GW observations, we discuss the methods for determining the equation of state of dark energy and for testing the co-evolution models of massive black holes. ASTROD-GW is an optimization of ASTROD to focus on the goal of detection of GWs. The mission orbits of the 3 spacecraft forming a nearly equilateral triangular array are chosen to be near the Sun-Earth Lagrange points L3, L4 and L5. The 3 spacecraft range interferometrically with one another with arm length about 260 million kilometers. With 52 times longer in arm length compared to that of LISA, the strain detection sensitivity is 52 times better toward larger wavelength. The scientific aim is focused for gravitational wave detection at low frequency. The science goals include detection of GWs from MBHs, and Extreme-Mass-Ratio Black Hole Inspirals (EMRI), and using these observations to find the evolution of the equation of state of dark energy and to explore the co-evolution of massive black holes with galaxies.     

LISA: Heliocentric formation design for the laser interferometer space antenna mission

The LISA (Laser Interferometer Space Antenna) mission has been selected by the European Space Agency’s Science Programme Committee as the third large-class mission of the Cosmic Vision Programme, addressing the science theme of the Gravitational Universe. With a planned launch date in 2034, LISA will be the first ever space-borne Gravitational Wave observatory, relying on laser interferometry between three spacecraft orbiting the Sun in a triangular formation. Airbus is currently leading an industrial Phase A system study on behalf of the European Space Agency. The paper will address the astrodynamics challenges associated with the LISA constellation design, driven by tight requirements on the geometric quality metrics of the near equilateral formation.     

The fundamental physics explorer: An ESA technology reference study

ESA technology reference studies are used as a process to identify key technologies and technical challenges of potential future missions not yet in the science programme. This paper reports on the study of the Fundamental Physics Explorer (FPE), a re-usable platform targeted to small missions testing fundamental laws of physics in space. The study addresses three specific areas of interest: special and general relativity tests based on atomic clocks, experiments on the Weak Equivalence Principle (WEP), and studies of Bose–Einstein condensates under microgravity conditions. Starting from preliminary science objectives and payload requirements, three reference missions in the small/medium class range are discussed, based on a re-adaptation of the LISA Pathfinder spacecraft. A 700/3600 km elliptic orbit has been selected to conduct clock tests of special and general relativity, a 700 km circular orbit to perform experiments on the Weak Equivalence Principle and to study Bose–Einstein condensates, each mission being based on a three-axis stabilised spacecraft. It was determined that adaptation of LISA Pathfinder would be required in order to meet the demands of the FPE missions. Moreover it was established that specific payload and spacecraft technology development would be required to realise such a programme.     

Testing Modified Newtonian Dynamics with LISA Pathfinder

We suggest that LISA Pathfinder, a technology demonstrator for the future gravitational wave observatory LISA, could be used to carry out a direct experimental test of Modified Newtonian Dynamics (MOND). The LISA Pathfinder spacecraft is currently being built and the launch date is just a few years away. No modifications of the spacecraft are required, nor any interference with its nominal mission. The basic concept is to fly LISA Pathfinder through the region around the Sun-Earth saddle point, in an extended mission phase, once the original mission goals are achieved. We examine various strategies to reach the saddle point, and find that the preferred strategy, yielding relatively short transfer times of just over 1 year, probably involves a lunar fly-by. LISA Pathfinder will be able to probe the intermediate MOND regime, i.e. the transition between deep MOND and Newtonian gravity. We present robust estimates of the anomalous gravity gradients that LISA Pathfinder should be exposed to, based on MONDian effects as derived from the Tensor-Vector-Scalar (TeVeS) theory. The spacecraft speed and spatial scale of the MOND signal combine in a way that the spectral signature of the signal falls precisely into LISA Pathfinder’s measurement bandwidth. We find that if the gravity gradiometer on-board the spacecraft achieves its currently predicted sensitivity, these anomalous gradients could not just be detected, but measured in some detail.     

Dynamics of a spacecraft and normalization around Lagrangian points in the Neptune–Triton system

The problem of a spacecraft orbiting the Neptune–Triton system is presented. The new ingredients in this restricted three body problem are the Neptune oblateness and the high inclined and retrograde motion of Triton. First we present some interesting simulations showing the role played by the oblateness on a Neptune’s satellite, disturbed by Triton. We also give an extensive numerical exploration in the case when the spacecraft orbits Triton, considering Sun, Neptune and its planetary oblateness as disturbers. In the plane a × I (a = semi-major axis, I = inclination), we give a plot of the stable regions where the massless body can survive for thousand of years. Retrograde and direct orbits were considered and as usual, the region of stability is much more significant for the case of direct orbit of the spacecraft (Triton’s orbit is retrograde). Next we explore the dynamics in a vicinity of the Lagrangian points. The Birkhoff normalization is constructed around L₂, followed by its reduction to the center manifold. In this reduced dynamics, a convenient Poincaré section shows the interplay of the Lyapunov and halo periodic orbits, Lissajous and quasi-halo tori as well as the stable and unstable manifolds of the planar Lyapunov orbit. To show the effect of the oblateness, the planar Lyapunov family emanating from the Lagrangian points and three-dimensional halo orbits are obtained by the numerical continuation method.     

Synergistic approach of asteroid exploitation and planetary protection

The asteroid and cometary impact hazard has long been recognised as an important issue requiring risk assessment and contingency planning. At the same time asteroids have also been acknowledged as possible sources of raw materials for future large-scale space engineering ventures. This paper explores possible synergies between these two apparently opposed views; planetary protection and space resource exploitation. In particular, the paper assumes a 5 tonne low-thrust spacecraft as a baseline for asteroid deflection and capture (or resource transport) missions. The system is assumed to land on the asteroid and provide a continuous thrust able to modify the orbit of the asteroid according to the mission objective. The paper analyses the capability of such a near-term system to provide both planetary protection and asteroid resources to Earth. Results show that a 5 tonne spacecraft could provide a high level of protection for modest impact hazards: airburst and local damage events (caused by 15–170 m diameter objects). At the same time, the same spacecraft could also be used to transport to bound Earth orbits significant quantities of material through judicious use of orbital dynamics and passively safe aero-capture manoeuvres or low energy ballistic capture. As will be shown, a 5 tonne low-thrust spacecraft could potentially transport between 12 and 350 times its own mass of asteroid resources by means of ballistic capture or aero-capture trajectories that pose very low dynamical pressures on the object.     

Simulation of ASTROD I test mass charging due to solar energetic particles and interplanetary electrons

As ASTROD I travels through space, its test mass will accrue charge due to exposure of the spacecraft to high-energy particles. This test mass charge will result in Coulomb forces between the test mass and the surrounding electrodes. In earlier work, we have used the GEANT 4 toolkit to simulate charging of the ASTROD test mass due to cosmic-ray protons of energies between 0.1 and 1000 GeV at solar maximum and at solar minimum. Here we use GEANT 4 to simulate the charging process due to solar energetic particle events and interplanetary electrons. We then estimate the test mass acceleration noise due to these fluxes. The predicted charging rates range from 2247 e⁺/s to 47,055 e⁺/s, at peak intensity, for the four largest SEP events in September and October 1989. Although the noise due to charging exceeds the ASTROD I budget for the two larger events, it can be suppressed through continuous discharging. The acceleration noise during the two small events is well below the design target. The charging rate of the ASTROD I test mass due to interplanetary electrons in this simulation is about −11% of the cosmic-ray protons at solar minimum, and over −37% at solar maximum. In addition to the Monte Carlo uncertainty, an error of ±30% in the net charging rates should be added to account for uncertainties in the spectra, physics models and geometry implementations.     

Finding the suitable drag-free acceleration noise level for future low-low satellite-to-satellite tracking geodesy missions

This paper evaluates the impact of residual acceleration noise on the estimation of the Earth’s time-varying gravity field for future low-low satellite-to-satellite tracking missions. The goal is to determine the maximum level of residual acceleration noise that does not adversely affect the estimation error. The Gravity Recovery And Climate Experiment (GRACE) has provided monthly average gravity field solutions in spherical harmonic coefficients for more than a decade. It provides information about land and ocean mass variations with a spatial resolution of ～350?km and with an accuracy within 2?cm throughout the entire Earth. GRACE Follow-on was launched in May 2018 to advance the work of GRACE and to test a new laser ranging interferometer, which measures the range between the two satellites with higher precision than the K-Band ranging system used in GRACE. Moreover, there have been simulation studies that show, an additional pair of satellites in an inclined orbit increases the sampling frequency and reduces temporal aliasing errors. Given the fact that future missions will likely continue to use the low-low satellite-to-satellite tracking formation with laser ranging interferometry, it is expected that the residual acceleration noise will become one of the largest error contributor for the time-variable gravity field solution. We evaluate three different levels of residual acceleration noise based on demonstrated drag-free systems to find a suitable drag-free performance target for upcoming geodesy missions. We analyze both a single collinear polar pair and the optimal double collinear pair of drag-free satellites and assume the use of a laser ranging interferometer. A partitioned best linear unbiased estimator that was developed, incorporating several novel features from the ground up is used to compute the solutions in terms of spherical harmonics. It was found that the suitable residual acceleration noise level is around 2?×?10^?12?ms^?2?Hz^?1/2. Decreasing the acceleration noise below this level did not result in more accurate gravity field solutions for the chosen mission architecture.     

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Small satellite survey mission to the second Earth moon

This paper presents an innovative space mission devoted to the survey of the small Earth companion asteroid by means of nano platforms. Also known as the second Earth moon, Cruithne, is the target identified for the mission. Both the trajectory to reach the target and a preliminary spacecraft budget are here detailed. The idea is to exploit high efficient ion thrusters to reduce the propellant mass fraction in such a high total impulse mission (of the order of 1e6 Ns). This approach allows for a 100 kg class spacecraft with a very small Earth escape energy (5 km²/s²) to reach the destination in about 320 days. The 31% propellant mass fraction allows for a payload mass fraction of the order of 8% and this is sufficient to embark on such a small spacecraft a couple of nano-satellites deployed once at the target to carry out a complete survey of the asteroid. Two 2U Cubesats are here considered as representative payload, but also other scientific payloads or different platforms might be considered according with the specific mission needs. The small spacecraft used to transfer these to the target guarantees the manoeuvre capabilities during the interplanetary journey, the protection against radiations along the path and the telecommunication relay functions for the data transmission with Earth stations. The approach outlined in the paper offers reliable solutions to the main issues associated with a deep space nano-satellite mission thus allowing the exploitation of distant targets by means of these tiny spacecraft. The study presents an innovative general strategy for the NEO observation and Cruithne is chosen as test bench. This target, however, mainly for its relevant inclination, requires a relatively large propellant mass fraction that can be reduced if low inclination asteroids are of interest. This might increase the payload mass fraction (e.g. additional Cubesats and/or additional scientific payloads on the main bus) for the same 100 kg class mission.     

Gravitational wave detection on the Moon and the moons of Mars

The Moon and the moons of Mars should be extremely quiet seismically and could therefore become sensitive gravitational wave detectors, if instrumented properly. Highly sensitive displacement sensors could be deployed on these planetary bodies to monitor the motion induced by gravitational waves. A superconducting displacement sensor with a 10-kg test mass cooled to 2 K will have an intrinsic instrument noise of 10⁻¹⁶ m Hz^−1/2. These sensors could be tuned to the lowest two quadrupole modes of the body or operated as a wideband detector below its fundamental mode. An interesting frequency range is 0.1–1 Hz, which will be missed by both the ground detectors on the Earth and LISA and would be the best window for searching for stochastic background gravitational waves. Phobos and Deimos have their lowest quadrupole modes at 0.2–0.3 Hz and could offer a sensitivity h_min ? 10⁻²² Hz^−1/2 within their resonance peaks, which is within two orders of magnitude from the goal of the Big Bang Observer (BBO). The lunar and Martian moon detectors would detect many interesting foreground sources in a new frequency window and could serve as a valuable precursor for BBO.     

Towards development of a consistent orbit series for TOPEX,Jason-1, and Jason-2

The TOPEX/Poseidon, Jason-1 and Jason-2 set of altimeter data now provide a time series of synoptic observations of the ocean that span nearly 17 years from the launch of TOPEX in 1992. The analysis of the altimeter data including the use of altimetry to monitor the global change in mean sea level requires a stable, accurate, and consistent orbit reference over the entire time span. In this paper, we describe the recomputation of a time series of orbits that rely on a consistent set of reference frames and geophysical models. The recomputed orbits adhere to the IERS 2003 standards for ocean and earth tides, use updates to the ITRF2005 reference frame for both the SLR and DORIS stations, apply GRACE-derived models for modeling of the static and time-variable gravity, implement the University College London (UCL) radiation pressure model for Jason-1, use improved troposphere modeling for the DORIS data, and apply the GOT4.7 ocean tide model for both dynamical ocean tide modeling and for ocean loading. The new TOPEX orbits have a mean SLR fit of 1.79 cm compared to 2.21 cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70 cm RMS, and the orbit accuracy is estimated at 1.5–2.0 cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08 cm RMS. The GSFC SLR/DORIS dynamic and reduced-dynamic orbits for Jason-2 agree radially with independent orbits from the CNES and JPL at 0.70–1.06 cm RMS. Applying these new orbits, and using the latest altimeter corrections for TOPEX, Jason-1, and Jason-2 from September 1992 to May 2009, we find a global rate in mean sea level of 3.0 ± 0.4 mm/yr.     

VHF radio response of the near Earth space during solar activity growth in October, 2003

The analysis of observations of very high frequency radio noise intensity at the middle latitude on a frequency f = 500 MHz from 14th till 26th of October, 2003 is presented. These data are compared with the solar radio bursts in the range of frequencies 1–14 MHz registered by RAD2 receiver of the WAVES device installed on board the WIND spacecraft.     

Photon counting detector for space debris laser tracking and lunar laser ranging

We are reporting on a design, construction and performance of solid state photon counting detector package which has been designed for laser tracking of space debris. The detector has been optimized for top photon detection efficiency and detection delay stability. The active area of the commercially available avalanche photodiode manufactured on Si (SAP500 supplied by Laser Components, Inc.) is circular with a diameter of 500 μm. The newly designed control circuit enables to operate the detection sensor at a broad range of biases 5–50 V above its breakdown voltage of 125 V. This permits to select a right trade-off between photon detection efficiency, timing resolution and dark count rate. The photon detection efficiency exceeds 70% at the wavelength of 532 nm. This is the highest photon detection efficiency ever reported for such a device. The timing properties of the detector have been investigated in detail. The timing resolution is better than 80 ps r.m.s, the detection delay is stable within units of picoseconds over several hours of operation. The detection delay stability in a sense of time deviation of 800 fs has been achieved. The temperature change of the detection delay is 0.5 ps/K. The detector has been tested as an echo signal detector in laser tracking of space debris at the satellite laser station in Graz, Austria. Its application in lunar laser ranging is under consideration by several laser stations.     

Periodic attitude motions along planar orbits in the elliptic restricted three-body problem

A way to improve the accuracy of the three-body problem model is taking into account the eccentricity of primary attractors. Elliptic Restricted Three-Body Problem (ER3BP) is a model for studying spacecraft trajectory within the three-body problem such that the orbital eccentricity of primaries is reflected in it. As the principal cause of perturbation in the employed dynamical model, the primaries eccentricity changes the structure of orbits compared to the ideal Circular Restricted Three-Body Problem (CR3BP). It also changes the attitude behavior of a spacecraft revolving along periodic orbits in this regime. In this paper, the coupled orbit-attitude dynamics of a spacecraft in the ER3BP are exploited to find precise periodic solutions as the spacecraft is considered to be in planar orbits around Lagrangian points and Distant Retrograde Orbits (DRO). Periodic solutions are repetitious behaviors in which spacecraft whole dynamics are repeated periodically, these periodic behaviors are the main interest of this study because they are beneficial for future mission designs and allow delineation of the system’s governing dynamics. Previous studies laid the foundation for spacecraft stability analysis or studying pitch motion of spacecraft in the ER3BP regime. While in this paper, at first, initial guesses for correction algorithms were derived through verified search methods, then correction algorithms were used to refine calculated orbit-attitude periodic behaviors. Periodic orbits and full periodic solutions are portrayed and compared to previous studies and simpler models. Natural periodic solutions are valuable information eventuate in the longer functional lifetime of spacecraft. Since the problem assumption considered in this paper is much closer to real mission conditions, these results may be the means to use natural bounded motions in the actual operational environment.     

Response of the mesopause temperatures to solar activity over Yakutia in 1999–2013

OH(6-2) rotational temperature trends and solar cycle effects are studied. Observations were carried out at the Maimaga station (63.04°N, 129.51°E) for the period August 1999 to March 2013. Measurements were conducted with an infrared spectrograph. Temperatures were determined from intensity ratios in the P branch of the OH band. The monthly average residuals of temperature after the subtraction of the mean seasonal variation were used for a search for the solar component of temperature response. The dependence of temperatures on solar activity has been investigated using the Ottawa 10.7 cm flux as a proxy. A linear regression fitting on residual temperatures yields a solar cycle coefficient of 4.24 ± 1.39 K/100 solar flux units (SFU). The cross-correlation analyses showed that changes of the residual temperature follow changes of solar activity with a quasi-two year delay (25 months). The temperature response at the delay of 25 months reaches 7 K/100 SFU. The possible reason of the observed delay can be an influence of quasi-biennial oscillations (QBO) of the atmosphere on the relation of temperature and solar activity. The value of the temperature trend after the subtraction of seasonal and solar components is not statistically significant.     

Dynamics and control of spacecraft hovering using the geomagnetic Lorentz force

To achieve hovering, a spacecraft thrusts continuously to induce an equilibrium state at a desired position. Due to the constraints on the quantity of propellant onboard, long-time hovering around low-Earth orbits (LEO) is hardly achievable using traditional chemical propulsion. The Lorentz force, acting on an electrostatically charged spacecraft as it moves through a planetary magnetic field, provides a new propellantless method for orbital maneuvers. This paper investigates the feasibility of using the induced Lorentz force as an auxiliary means of propulsion for spacecraft hovering. Assuming that the Earth’s magnetic field is a dipole that rotates with the Earth, a dynamical model that characterizes the relative motion of Lorentz spacecraft is derived to analyze the required open-loop control acceleration for hovering. Based on this dynamical model, we first present the hovering configurations that could achieve propellantless hovering and the corresponding required specific charge of a Lorentz spacecraft. For other configurations, optimal open-loop control laws that minimize the control energy consumption are designed. Likewise, the optimal trajectories of required specific charge and control acceleration are both presented. The effect of orbital inclination on the expenditure of control energy is also analyzed. Further, we also develop a closed-loop control approach for propellantless hovering. Numerical results prove the validity of proposed control methods for hovering and show that hovering around low-Earth orbits would be achievable if the required specific charge of a Lorentz spacecraft becomes feasible in the future. Typically, hovering radially several kilometers above a target in LEO requires specific charges on the order of 0.1 C/kg.     

Intracellular molecular distributions in spacecraft experiments in orbit around Earth

It is possible that the nucleolous inside the cell plays the role of a “gravity receptor”. Furthermore, cells up to 10 μm in diameter can demonstrate some effect due to the redistribution of mitochondria or nucleolous. Effects of gravity should be present in various cell systems where larger objects such as the ribosomes move from cell to cell. In this paper we study the effects of gravity on cells. In particular, we examine the resulting intracellular molecular distribution due to Brownian motion and the ordered distribution of molecules under the action of gravity, where n₀ is the number per unit volume at certain level, and n is the number per unit volume above that level. This is an experiment that takes place at a certain orbital altitude in a spacecraft in orbit around Earth, where the acceleration due to the central field is corrected for the oblateness and also the rotation of the Earth. We found that equatorial circular and elliptical orbits have the highest n/n₀ ratios. This experiment takes place in circular and elliptical orbits, with eccentricities e = 0, 0.1 and involves a bacterial cell at an orbital altitude of 300 km. We found that n/n₀ = 1.00299 and 1.0037 respectively, which is still a 0.6–0.7 % higher than n/n₀ = 0.0996685 calculated on the surface of the Earth. Examining mitochondria in similar orbital experiments we found that equatorial orbits result to higher n/n₀ ratios. In particular, we found that n/n₀ = 8.38119, where an elliptical orbit of eccentricity e = 0.1 results to n/n₀ = 13.8525. Both are high above 100%, signifying the importance of Brownian motion over gravity. Our results are of interest to biomedical applications. Molecular concentrations are important for various processes such as the embryogenesis, positional homeostasis and its relation to cell energy expenditure, cell torque, cell deformation, and more. These results indicate that statistical molecular distributions play an important role for the recognition of a particular environment by the cell, in biological space experiment to come.