A Cubesat Payload for Exoplanet Detection
Next Article in Journal
Performance Analysis of Different Backoff Algorithms for WBAN-Based Emerging Sensor Networks
Previous Article in Journal
An Adaptive Damping Network Designed for Strapdown Fiber Optic Gyrocompass System for Ships
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Cubesat Payload for Exoplanet Detection

1
Department of Industrial Engineering—Aerospace Division, University of Naples Federico II, 80125 Naples, Italy
2
INAF—Astrophysical Observatory of Arcetri, Largo E. Fermi 5, 50125 Florence, Italy
3
INAF—Astrophysical Observatory of Capodimonte, Salita Moiariello, 16, 80131 Naples, Italy
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(3), 493; https://doi.org/10.3390/s17030493
Submission received: 5 January 2017 / Revised: 18 February 2017 / Accepted: 23 February 2017 / Published: 2 March 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
The search for undiscovered planets outside the solar system is a scientific topic that is rapidly spreading into the astrophysical and engineering communities. In this framework, the design of an innovative payload to detect exoplanets from a nano-sized space platform, like a 3U cubesat, is presented. The selected detection method is photometric transit, and the payload aims to detect flux decrements down to ~0.01% with a precision of 12 ppm. The payload design is also aimed at false positive recognition. The solution consists of a four-facets pyramid on the top of the payload, to allow for measurement redundancy and low-resolution spectral dispersion of the star images. The innovative concept is the use of a small and cheap platform for a relevant astronomical mission. The faintest observable target star has V-magnitude equal to 3.38. Despite missions aimed at ultra-precise photometry from microsatellites (e.g., MOST, BRITE), the transit of exoplanets orbiting very bright stars has not yet been surveyed photometrically from space, since any observation from a small/medium sized (30 cm optical aperture) telescope would saturate the detector. This cubesat mission can provide these missing measurements. This work is set up as a demonstrative project to verify the feasibility of the payload concept.

1. Introduction

An exoplanet is a planet that is orbiting around a star other than the Sun. An Earth-like planet is a habitable planet of approximately 1 Earth mass and 1 Earth radius, and in an Earth-like orbit around a Sun-like star at a distance of roughly 1 AU [1].
By now (February 2016) the number of total confirmed exoplanets (From http://exoplanetarchive.ipac.caltech.edu/. The criteria to list an exoplanet as confirmed are as follows: the mass is equal or less than 30 Jupiter masses; the planet is associated with a host star (not free floating); sufficient follow-up observations and validation have been undertaken to deem that the object being a false positive is unlikely; the above information along with further orbital and/or physical properties are available in peer-reviewed publications) is 3449. Among these, the number of planets discovered by the photometric transit method is 2718 (From http://exoplanetarchive.ipac.caltech.edu/), showing that this method is very effective. Of those 2718 planets, the number of terrestrial planets, with planet mass between 0.6 Earth mass and 10 Earth masses, is 80 (From http://exoplanetarchive.ipac.caltech.edu/).
The photometric transit technique consists of measuring the reduction of the star photometric flux when the planet passes in front of the host star as seen from the Earth. In order to realistically observe the exoplanet transit, the exoplanet orbit around the host star must be aligned with the observer. The measurement of the photometric star flux as a function of time is depicted in the so-called light curve of the target star.
The objective of detecting exoplanet transit is becoming a priority among the current astronomical space missions. The following is a brief list of the more relevant past/present/future missions regarding exoplanet transit detection.
CoRoT is the first mission ever designed to detect exoplanet transit from space, led by ESA and by CNES (Centre national d’études spatiales) from 2006 to 2013 [2]. The 27 cm diameter telescope observed in the band from 370 nm to 950 nm, and a one facet prism was installed before the CCDs used for transit detection [3]. Kepler is the NASA’s first mission (from 2009 to 2013) able to find Earth-sized planets around Sun-like stars. The 0.95 m diameter telescope has a FOV (field of view) of 105 deg2, an array of 42 CCDs observing in the wavelength band from 400 nm to 850 nm, with instrument precision of 10−5. It continuously points at a single star field in the Cygnus-Lyra region, with a large number of stars but off the galactic plane to reduce field confusion [4]. TESS is the first spaceborne all-sky exoplanet transit survey from NASA-MIT (Massachusetts Institute of Technology), scheduled in 2017, observing in the band from 600 nm to 1000 nm [5]. CHEOPS is the first small class mission from ESA (European Space Agency), scheduled in 2017, to detect exoplanet transit on stars with V-magnitude between six and twelve anywhere in the sky, observing in the band from 0.4 to 1.1 µm [6]. PLATO is a medium class exoplanet transit detection mission from ESA (scheduled in 2024) in orbit around L2, looking at stars with a V-magnitude between 4 and 16 [7].
There are also cubesat projects that aim at astronomical missions regarding exoplanet detection; a brief project list is provided in the following.
The most relevant cubesat project is the pioneer Exoplanetsat, a 3U cubesat from MIT and DraperLab. As a pathfinder mission, the mission objective is to detect exoplanet transit around bright stars with V-magnitude less than four. The key aspect of this project is a fine target pointing achieved by a closed control loop on the detector position to compensate for the satellite jitter [8]. The project evolved in ASTERIA, a 6U cubesat with improved photometric capabilities. MDOT is a 6U cubesat project from Stanford University in a HEO (High Earth Orbit) orbit. The project aims at taking a direct image of exo-zodiacal dust and transiting exoplanets using an occulter. The exo-zodiacal light is a portion of the star light that is scattered by the micrometric dust grains in the planetary system plane plus a dust thermal emission [9]. To detect the planet transit the central star light must be suppressed through the occulter. CANYVAL-X is a NASA mission and a formation flying demonstrator. The mission consists of 2U cubesat, used for the exoplanet detection, plus a 1U cubesat, used as occulter [10]. DeMi is a MIT study about a Cubesat Deformable Mirror Demonstrator in LEO (Low Earth Orbit) to enable space-based coronagraphic direct imaging of exoplanets from a 3U cubesat platform [11]. Centaur is a pathfinder mission to directly image exoplanets, using a new kilo-deformable mirror (1024 actuators), looking at Alpha Centauri A&B. PicSat is a 3U cubesat with the purpose of monitoring Beta Pictoris, a A6V star of V = 3.86 magnitude. The payload consists of a 35 mm effective aperture objective and a single pixel avalanche photodiode. A single-mode fiber guides the star from the focal plane to the photodiode. As for an exoplanet project, a closed control loop is used to achieve a fine pointing [12].
In a mission devoted to exoplanet detection through photometric transit there are several cases in which the target star brightness decreases for reasons other than a planet transit. The false positive signal sources belong to two possible groups: hardware faults and astronomical scenarios. The second group refers to astronomical objects other than planets that transit in front of the target star. The distinctive element of a planet transit is that the star light curve decreases uniformly in the entire visible band. Unlike stars, the rocky planets temperature produces an approximately uniform spectral emission or absorption in the visible band.
The current approaches to discriminate the true signal from false positives are here briefly described.
Ground based complementary observations are often used to confirm a measurement, like photometry, high resolution spectroscopy, and radial velocity measurements [13].
Another strategy is to evaluate the pixel level data and to identify which pixels in and around the target contain the transit signal. If the event is a false positive, the pixel location of the transit signal does not coincide with the pixel location of the target flux. It works well just in cases of an eclipsing object close enough to the observed target, so that the transit signal is not too much diffused [14].
A quite demanding strategy in terms of database size is ephemerides matching, which consists of first creating catalogs of transiting planets, eclipsing binaries, and other variable stars in the instrument field of view; secondly defining a criterion to compare objects, according to their period and epochs; and then looking for objects in the catalogs that fulfill the matching criterion according to their ephemerides. Found objects indicate at least one false positive event [14].
A more recent approach is the probabilistic method, which consists of computing in a Bayesian framework the probability of the planet transit scenario against an exhaustive set of false positive scenarios. If the planet scenario is the highest probable one, then the planet is validated [13].
The estimated false detection rate of Kepler measurements is in the order of 50% (e.g., giant planet false positive rate indicated in [15]). Most of Kepler false positives are produced by background eclipsing binaries and planets transiting a star that is physically bound with the target star [13].
The false positive rate of CoRoT mission is 83%. In this mission the identification of false positives is realized only by the light-curves analysis, evaluating the transit depth and duration, the curve shape and the presence of color signature [16]. A CoRoT astronomical field is denser than Kepler, so the probability of observing a false positive transit is higher. Moreover, the CoRoT observation approach neglects the centroid follow up during the transit, i.e., measuring the centroid shift during the transit helps in rejecting possible background eclipsing binaries [13,17]. However, the main reason for this rate in CoRoT is the large PSF size of each star, such that the detector saturation was avoided but the overlapping rate of star PSFs was increased.
This project combines the growing interest in the exoplanet search with the increasing success of the cubesat platform. The project aims at detecting exoplanet photometric transit, focusing on very bright stars, and discriminating a false positive detection.

2. The Satellite Project

There are two key aspects when observing an exoplanet transit from space. Space-based observation ensures a high level of precision and consistency for the transit measurement, as it is required due to the tiny flux decrement (0.01%) of a Sun-like star during the transit of an Earth sized planet [1]. The reason is that observations from space do not suffer from atmospheric disturbances (seeing variations or scintillation) and limits, day-time cycle variation, moonlight, weather factors. A drawback of ground telescope observations is the limited scheduled observing time, while a space mission can be devoted even to the observation of a single target over the project lifespan. Even for ground-based instruments that are devoted to exoplanets’ search, the smallest detectable planets are Neptune-sized.

2.1. Mission Analysis

The project objective is to demonstrate the feasibility of a high performance astronomical mission based on a low-cost space platform, and also to show a proof of concept for a valid technique against false positive signals.
The project’s key concepts are the following.
(1)
Exoplanet detection (goal Earth-like planets around Sun-like stars)
(2)
Photometric transit method
(3)
Identification of false positive signals
(4)
3U cubesat platform
Following a system engineering approach [18], each key concept defines a system driving requirement, which in turn determines one or more satellite subsystem requirements. The system driving requirements are shown below with the nomenclature R.x, where x is the key concept number. Table 1 summarizes the satellite subsystem requirements.
R.1
The amplitude of the signal to be detected is 84 ppm. During the transit the ratio between the combined flux from star and planet (F(t)) and the unobstructed star flux (Fs) is approximated by the transit depth expression δ. The transit depth δ is proportional to the square ratio between the planet radius (Rplanet) and the star radius (Rstar) [1], Equation (1) (neglecting at this time the influence of the star limb darkening on the transit depth). The assumed star is the Sun and the assumed planet is the Earth, Table 2.
F ( t ) F s ( R planet R star ) 2
The assumed planet transit duration is 6.5 h [19]. To get a 7 σ certainty of the measured transit [20], the maximum level of tolerable noise is 12 ppm. To be conservative, this requirement is applied to the observation time equal to the duty cycle time, without taking into account the multiple transit observation. A thermal analysis will be included in the error budget, thus reducing the requirement margin. The precision of the measurement requires high-stability pointing along the entire observation. As first choice (in analogy with ExoplanetSat example) the value of 5 arcsec is needed for pointing stability.
R.2
The required payload components are the objective and the detector (photometric payload)
R.3
To discriminate the false positive signal, the photometric signal is slightly dispersed. The observation band is the visible. The required spectral resolution is low (3 bands from 400 nm to 850 nm), since the measurement purpose is to monitor the centroid location of the dispersed signal.
R.4
The standard cubesat structural limits in terms of dimensions and weight must be fulfilled. The 3U-cubesat dimensions are 10 cm × 10 cm × 30 cm (34 cm in 3U+ configuration) with a mass up to 4 kg [21], Figure 1. These small satellites require an orbital deployer (e.g., P-POD from CalPoly) to ensure that the cubesat is safely stored and correctly launched from the launcher. The advantage of choosing a cubesat platform is its straightforward design, realization, and test with a cheap budget. The drawbacks are limited size and mass, and then limited resources for power, computing, and attitude control. The technical challenge is then to fulfill the mission objective through the limited capabilities of the cubesat space platform. Moreover, the requirement of 3U size comes from the consideration that further mass increment (and thus unit increment, till about 8U) would not provide free space for an objective aperture wider than 100 mm (assuming to use a circular aperture shape). The minimum estimated size for the payload is 1 unit and at least 2 other units are required for the satellite’s subsystems. Three units is then considered the most suitable number, as seen in Figure 2.

2.2. Satellite Orbit and Target

Usually the selection of the orbit is influenced by multiple aspects, like the available options from the launch providers, the target position, analogy with existing missions, and considerations on the space radiation environment.
For exoplanet detection missions the orbit options depend on the platform size; for large and medium sized missions, L2 orbit and heliocentric orbits are suitable, while for medium and small platforms, a LEO orbit is common and feasible. If the satellite orbit is far away from Earth, it is possible to avoid the day–night cycle of observation. In the case of a low altitude (400 km/600 km) LEO orbit, instead, the best observing time is the orbit night, and any transit observation is possible only for an observer that lies in the penumbra of the transiting exoplanet.
The launch options for cubesats are the followings.
  • as a piggyback payload from standard launchers,
  • through the cubesat launch service from the ISS (Kibo module),
  • through a private service from a dedicated satellite (e.g., GAUSS Unisat),
  • as a primary payload from a launcher for small satellites (future option).
The launch cost range goes from 40,000 € for a piggyback launch to 150,000 € for a private satellite, or even for free in case of launch campaign for university satellites (QB50 program). Moreover, many start-up companies are building new concept rockets to launch only small satellites with a low cost and easier application process.
The assumed orbit is circular sun-synchronous at 600 km altitude and 98° inclination, as one of the available options from the launch providers (GAUSS private company, DNEPR rocket). From this assumption the computed eclipse time is 31 min over 97 min orbital period duration. Further analyses will be conducted to identify the requirements about the straylight during the transit observation. These future considerations could affect the orbit altitude and platform orientation options.
The number of known dwarf stars (luminosity class V) of V-magnitude less than four is 115 (91 dwarfs have V-magnitude less than 3.8) (From SIMBAD astronomical database, http://simbad.u-strasbg.fr/simbad/). Currently (February 2017) there are four confirmed exoplanets orbiting dwarf stars of V-magnitude less than 4, and six unconfirmed exoplanets orbiting around the same type of star (From http://exoplanet.eu/). None of them has been ever observed with the photometric method from space. One of the confirmed exoplanet host stars is Alpha Centauri B (planet Alpha Centauri Bb), and further investigations are required for the star companion Alpha Centauri A. The first option target is Alpha Centauri AB, which is the closest-to-the-Sun binary star system that could host a planet in the habitable zone, and it is the most feasible target to be observed from a small platform [22]. Alpha Centauri A’s V-magnitude is 0.01, and Alpha Centauri B’s V-magnitude is 1.33. It is still to be defined the observing strategy, e.g., switching from one star to the companion, or choosing one of the two as a permanent target.

2.3. Cubesat Design

The cubesat platform consists of several subsystems. The payload subsystem will be the main focus of this paper, and a preliminary design for the other subsystems will be described.
At the present design stage, the payload consists of the following parts, as seen in Figure 3.
  • a glass pyramid, with a round base and four facets, as seen in Figure 4,
  • a commercial objective,
  • a scientific detector, cmos (1),
  • a second detector used to close the attitude control loop, cmos (2),
  • a two-axis piezo-stage behind the detectors to compensate for the spacecraft jitter in a closed control loop.
The pyramid is located before the objective. The four pyramid facets form four images, and the four sky fields are overlapped in the sky area including the target star. The four images of the target star and surrounding sky field identify four detector areas, in which the photometric measurements take place. These are the so-called detector photometric windows, whose size can be adjusted according to the available data rate.
The pyramid allows for four simultaneous recordings of the transit, so that the overall transit measurement is free from instrumental fault events due to the detector. The probability (P4) that an instrumental fault event happens in all of the four windows is equal to the fourth power of the probability (P1) of an instrumental fault in one window; thus, P4 is greatly lower than P1.
The pyramid also has a dispersive power, and its spectral resolution is determined by the pyramid base angle. The dispersion direction in each photometric window is rotated 90° with respect to the neighboring photometric windows. Before and during the transit the centroid of the photometric signal is computed in each window. If a planet transit occurs, the signal reduction is spectrally uniform and the relative distance between the windows centroids is kept constant. In case of an astronomical false positive event, the transit could be made by the star companion of a binary system, and the signal reduction is not spectrally uniform, due to the transiting star spectral emission and absorption. In this case the relative distances between the window centroids change, as seen in Figure 5.
The pyramid material first option is the commercial BK7 glass; a radiation resistant glass is also considered, i.e., BK7 G18 SCHOOT, and no pyramid design variation is required.
The design of the pyramid consists of the definition of the pyramid base angle according to the required rays deviation and image dispersion. To avoid saturation and to have a broad detector dynamic range, the star spot is defocused. The star defocus is assumed to be 10 × 10 pixels; a simplified sketch of the optical system has been realized in ZEMAX, thanks to the custom zemax library PAM2R [23], as seen in Figure 6. Other assumptions are: the photometric window size is 40 × 40 pixels, and the ray deviation from the optical axis is 300 pixels, in the two perpendicular directions on the detector plane. From geometric considerations the angular deviation α is 3.64°. Considering the material reference refractive index n equal to 1.52, and according to the simplified Equation (2) [24], the pyramid base angle θ is 7°.
Considering the refractive index variation in the band from 410 nm to 850 nm, each star spot is dispersed over 11 pixels, so that a bandwidth of 150 nm corresponds to 3.5 pixels. Adding in quadrature the defocus and the dispersion contribution, the total elongated star dimension is 15 pixels.
θ = α   ( n 1 )
The piezoelectric stage is located behind the focal plane. It is the key element to reduce the spacecraft jitter to few ppm of noise. The stage works in a closed loop that keeps the target image on the same pixels during the scientific exposure. The detector involved in the closed loop is an auxiliary sensor with a pixel size smaller than the scientific detector pixel size, and thus it reaches a higher frame rate. According to [25] this technique (together with the MIT driver code) let to reach the pointing accuracy of 2.3 arcsec.
The first stage option is a custom Physik Instrumente (PI) stage model, with custom changes to fit the cubesat platform. The two axes nanometric PI stage, PI P-733.2, has a travel range of 100 µm (for each axis) and a resolution of 0.1 nm. Another option is XY200M, a two axes piezo stage from CEDRAT Technologies, with a travel range of 200 µm (no load value) and a resolution of 20 nm.
The considered commercial objective is the ZEISS Planar T 1.4/85, with a focal length of 85 mm and the aperture diameter of 60.7 mm. The projected image on the focal plane has a diameter of 43 mm [26]. The objective has been tested for space stresses through a scientific research work from ZEISS laboratory [27].
The selected scientific detector is the HAS2 image sensor. It is the High Accuracy Star tracker CMOS image sensor from ON Semiconductor. A second detector is used to achieve a fine attitude control (e.g., e2V EV76C454 CMOS sensor). The star tracking activity is performed from the HAS2 sensor on the image resulting from the pyramid and the objective combination, and thus the adoption of a cross-correlation algorithm is foreseen. The HAS2 features are listed in Table 3. The HAS temperature features are: dark current doubling for sensor temperature increment of 5.8 °C (average), and voltage-temperature variation of −4.64 mV/°C. If needed, temperature control based on a Peltier module will be considered. The HAS2 sensor has been tested for functionatility up to 300 krad, and up to 42 krad the functionality is guaranteed. Additional considerations and detailed studies on this topic will be carried out in the next steps of the project. However, even at this stage, it can be seen that this sensor is widely used in space applications such as star sensors for attitude determination whose lifetime is several years. The latter performance, even scaled down at the level of a cubesat system, will be in excess of the cubesat satellite’s nominal lifetime.
The payload subsystem should fit in one cubesat unit, that is in a 100 × 100 × 120 mm3 cube, and should weight no more than 1330 g. The pyramid base height is assumed as 10 mm and the pyramid height is 3.8 mm. Table 4 lists the sizes and weights of each payload component; each X- and Y- size must be less than or equal to 100 mm, the sum of Z sizes must be less than or equal to 120 mm, and the sum of the weights must be less than 1330 g.
The attitude determination and control consists of an off-the-shelf module with at least 3 reaction wheels and 3 axis magnetometers, for instance MAI-400 from Maryland Aerospace. The attitude determination is mainly achieved through the star tracker camera HAS2. The fine target pointing is achieved through a closed loop control system realized with a nanometric 2-axis stage and a second auxiliary detector (e.g., PI stage and e2V EV76C454 CMOS, as described above).
The on-board scientific and attitude data handling will be realized through a commercial off the shelf single board computer suitable for cubesats.
The electrical power is supplied by deployable solar panels, and the likely configuration option is to deploy the panels and to arrange them in one single plane. The thermal control will be passive.

3. Star Visual Magnitude Limit

The most demanding requirement refers to the capability of detecting the signal with amplitude equal to 84 ppm and with a maximum noise level of 12 ppm. There is a limit in the target star V-magnitude, below which this requirement is fulfilled. This limit is computed in the following steps.
The first step is the evaluation of the photon flux from typical target stars. The considered V-magnitude (vmag) range is from 0 to 5. The photon flux (sflux) is given in photons per seconds, scaling the reference photon flux of Vega in the visible band. The assumed optical aperture is the one from Zeiss objective (2893 mm2).
The overall signal level is then computed in terms of detected photons, considering a range of possible exposure time (texp) from 0.01 s to 10 s, and of a number of co-added frames (#Frames) from 1 to 10,000, Equation (3).
The second step is the noise computation. The noise sources taken into account are: photonic noise (Nphot), dark current noise (Ndark), and read out noise (Nron) [28], Equations (4)–(6). The detector choice sets the dark current value (Idark), expressed in electrons per pixel per seconds, and the read out characteristic (RON), expressed in electrons per pixel.
The third step is the S/N ratio computation, Equation (7). The dominant noise source is the photon noise, and the S/N can be approximated as the square root of the signal, Equation (8). The requirement of maximum noise level of 12 ppm corresponds to a S/N of 105, and requires a star signal level of 1010 photons. The star spot is intentionally defocused, and the star photons are spread on 10 × 10 pixels, leading to 108 photons per pixel required.
s = s flux × t exp × # Frames
N phot 2 = s
N dark 2 = I dark × t exp × # Frames × # pixel
N ron 2 = RON 2 × # pixel × # Frames
S N =   s N phot 2 + N dark 2 + N ron 2
S N   s ( If readout noise and dark noise are negligible )
The constraints related to the achievable S/N are the following, shown in Table 5.
  • The observation time on each orbit (i.e., the product between the exposure time of a single frame and the number of co-added frames) must be less than or equal to the duty cycle time. The total transit observations should include one full transit time before and after the transit, to better estimate the transit baseline.
  • The exposure time of each frame must be less than the time required to saturate the detector. This constrain is expressed as the number of photons corresponding to the saturation voltage, taking into account the detector conversion factor and the quantum efficiency.
The next computation steps, from number 4 to number 8, identify the star magnitude limit according to the constraints (see Table 6).
The results of the computation show that star V-magnitude equal to 3 is the limit, since the observation of a star of fourth V-magnitude would require an exposure time longer than 4 s. A finer analysis in the magnitude range from 3 to 4 shows that the limiting V-magnitude is 3.38 (Figure 7).
The data rate computation is relative to four 40 × 40 photometric windows and a 12 bit resolution detector; for a target star of V-magnitude 3 (i.e., 453 frames per orbit) the data rate is 34.8 Mbit per orbit. The data will be stored, waiting for the ground station to be in view.

4. Discussion and Conclusions

Payload design for exoplanet transit detection from a 3U cubesat platform is a challenging mission, but it is consistent with the trend of using cheaper space platforms for highly demanding scientific missions. The objective of the project is to analyze the feasibility of this kind of mission, and also of using the pyramid to provide measurements free of false positives.
The 3U cubesat transversal size limits the maximum aperture of the objective, and this affects the V-magnitude of the fainter observable star. From the approximated computations illustrated above, the possible targets are very bright stars with V-magnitudes equal to or lower than 3.38. Even though this is a tight range of possible targets, photometric transit observation from space has never been done for these bright targets.
The attitude subsystem and the target pointing are essential in all the demanding applications for cubesats, especially for astronomical missions. The drawback of aiming at these high-performance objectives is that much customization could be required, which cancels out the low-cost advantage of the cubesat project. The attitude and pointing subsystem could be the subject of future research activities.

Acknowledgments

The funds to publish this article are provided by the annual grant given from University of Naples Federico II to each Ph.D. student attending the Ph.D. program in industrial engineering.

Author Contributions

Ernesto Oliva and Andrea Tozzi contributed to the optical layout of the payload, in particular, Oliva suggested the idea of the pyramid, and they also provided astronomical insights. Schipani contributed in the training and support for the use of PAM2R library and in the pyramid design, and as overall adviser. Accardo and Rufino contributed in guiding the research activity, with tutorship for the entire project, and providing extensive support and insights in the engineering design. Marcella Iuzzolino investigated and developed the topic, matching engineering and astronomical fields.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

FOV Field of view
MITMassachusetts Institute of Technology
NASANational Aeronautics and Space Administration
ESAEuropean Space Agency
CNESCentre national d’études spatiales
TESSTransiting Exoplanet Survey Satellite
CHEOPSCHaracterizing ExOPlanets Satellite
PLATOPLAnetary Transits and Oscillations of stars
CoRoTConvection, Rotation and planetary Transits
HEOHigh Earth Orbit
LEO Low Earth Orbit
CCDCharge Coupled Device
PIPhysik Instrumente
ppmparts per million
S/NSignal to Noise ratio
vmagvisual magnitude

References

  1. Seager, S. Introduction to exoplanets. In Exoplanets, 1st ed.; University of Arizona Press: Tucson, AZ, USA, 2011. [Google Scholar]
  2. Aigrain, S.; Barge, P. The CoRoT Exoplanet Programme: Exploring the Gas-Giant/Terrestrial Planet Transition. In Proceedings of the 14th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, Leicester, UK, 7–11 August 2007; Volume 384, p. 270.
  3. Auvergne, M.; Bodin, P. The CoRoT Satellite in Flight: Description and Performance. Astron. Astrophys. 2009, 506, 411–424. [Google Scholar] [CrossRef]
  4. Borucki, W.; Koch, D. KEPLER: Search for Earth-Size Planets in the Habitable Zone. In Proceedings of the International Astronomical Union Symposium, Ioannina, Greece, 15–19 September 2008; Cambridge University Press: Cambridge, UK, 2008; Volume 4, pp. 289–299. [Google Scholar]
  5. Ricker, R.; Joshua, N.W. Transiting Exoplanet Survey Satellite. J. Astron. Telesc. Instrum. Syst. 2014. [Google Scholar] [CrossRef]
  6. Cessa, V.; Beck, T. Cheops: A space telescope for ultra-high precision photometry of exoplanet transits. In Proceedings of the International Conference on Space Optics, Tenerife, Canary Island, Spain, 7–10 October 2014.
  7. Rauer, H.; Catala, C.; Aerts, C.; Appourchaux, T.; Benz, W.; Brandeker, A.; Christensen-Dalsgaard, J.; Deleuil, M.; Gizon, L.; Goupil, M.J.; et al. The PLATO 2.0 Mission. Exp. Astron. 2014, 38, 249–330. [Google Scholar] [CrossRef]
  8. Smith, M.; Seager, S. The ExoplanetSat mission to detect transiting exoplanets with a cubesat space telescope. In Proceedings of the 25th Annual AIAA/USU Conference on Small Satellites, Logan, UT, USA, 8–11 August 2011.
  9. Grün, E.; Gustafson, B.A.; Dermott, S.; Fechtig, H. Interplanetary Dust, Astronomy and Astrophysics Library; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  10. Park, J.; Park, S. Mission Analysis and CubeSat Design for CANYVAL-X Mission. In Proceedings of the Space Operations 2016 Conference, Deajeon, Korea, 16–20 May 2016.
  11. Cahoy, K.; Marinan, A. CubeSat deformable mirror demonstration. In Proceedings of the Space Telescopes and Instrumentation 2012: Optical, Infrared, and Millimeter Wave, Amsterdam, The Netherlands, 1–6 July 2012.
  12. Nowak, M.; Lacour, S. Reaching Sub-Millimag Photometric Precision on Beta Pictoris with a Nanosat: The PicSat Mission. In Proceedings of the Space Telescopes and Instrumentation, Edinburgh, UK, 26–29 June 2016.
  13. Santerne, A.; Fressin, F. The Contribution of Secondary Eclipses as Astrophysical False Positives to Exoplanet Transit Surveys. Astron. Astrophys. 2013, 557. [Google Scholar] [CrossRef]
  14. Coughlin, J.; Thompson, S. Contamination in the Kepler Field. Identification of 685 KOIs as False Positives via Ephemeris Matching Based on Q1–Q12 Data. Astron. J. 2014, 147, 119. [Google Scholar] [CrossRef]
  15. Santerne, A.; Moutou, C.; Tsantaki, M.; Bouchy, F.; Hébrard, G.; Adibekyan, V.; Almenara, J.M.; Amard, L.; Barros, S.C.; Boisse, I.; et al. SOPHIE velocimetry of Kepler transit candidates. Astron. Astrophys. 2016, 587, A64. [Google Scholar] [CrossRef]
  16. Almenara, J.M.; Deeg, H.J. Rate and Nature of False Positives in the CoRoT Exoplanet Search. Astron. Astrophys. 2009, 506, 337–341. [Google Scholar] [CrossRef]
  17. Bryson, S.; Jenkins, J.M. Identification of Background False Positives from Kepler Data. Publ. Astron. Soc. Pac. 2013, 125, 889. [Google Scholar] [CrossRef]
  18. Wertz, J.; Larson, W. Space Mission Analysis and Design, 3rd ed.; Microcosm Press: El Segundo, CA, USA, 1999; pp. 1–94. [Google Scholar]
  19. Kepler Mission NASA Website. Available online: https://kepler.nasa.gov/science/about/characteristicsOfTransits/transitDetectionAndSimulation/ (accessed on 4 January 2017).
  20. Jenkins, J.M.; Caldwell, D.A. Some Tests to Establish Confidence in Planets Discovered by Transit Photometry. Astrophys. J. 2002, 564, 495. [Google Scholar] [CrossRef]
  21. CubeSat Design Specification. Available online: http://www.cubesat.org/resources/ (accessed on 4 January 2017).
  22. Bendek, E.; Belikov, R.; Thomas, S. Centaur, a scientific and technology pathfinder for direct imaging exoplanets missions. In Proceedings of the 4th Interplanetary Cubesat Workshop 2015, London, UK, 4–25 May 2015.
  23. Antichi, J.; Munari, M. Modeling Pyramidal Sensors in Ray-Tracing Software by a Suitable User-Defined Surface. J. Astron. Telesc. Instrum. Syst. 2016, 2. [Google Scholar] [CrossRef]
  24. Hecht, E. Optics, 4th ed.; Addison-Wesley: San Francisco, CA, USA, 2002. [Google Scholar]
  25. Pong, C. High-Precision Pointing and Attitude Estimation and Control Algorithms for Hardware-Constrained Spacecraft. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2014. [Google Scholar]
  26. Zeiss Website. Available online: http://www.zeiss.it/camera-lenses/it_it/camera_lenses/slr-lenses/planart1485.html (accessed on 4 January 2017).
  27. Zeiss Camera Lens News. Available online: http://lenspire.zeiss.com/en/wp-content/uploads/sites/2/2009/10/cln33_en_web.pdf (accessed on 4 January 2017).
  28. Howell, S.B. Handbook of CCD Astronomy, 2nd ed.; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
Figure 1. Layout of a 3U+ cubesat platform, according to standard cubesat specification.
Figure 1. Layout of a 3U+ cubesat platform, according to standard cubesat specification.
Sensors 17 00493 g001
Figure 2. Cubesat unit tradeoff. The red line represents the available objective aperture size (D) of the payload as a function of the number of units. The blue dash-dot line is a first approximation cost (C) as a function of the unit number. The cost trend follows the linear mass increment trend as the unit number grows. Since the aperture shape in question is circular, D is unchanged until 8 units, which is the minimum number of units to build a square cubesat with a face size of 2 units. The black dashed line corresponds to the 3U cubesat size.
Figure 2. Cubesat unit tradeoff. The red line represents the available objective aperture size (D) of the payload as a function of the number of units. The blue dash-dot line is a first approximation cost (C) as a function of the unit number. The cost trend follows the linear mass increment trend as the unit number grows. Since the aperture shape in question is circular, D is unchanged until 8 units, which is the minimum number of units to build a square cubesat with a face size of 2 units. The black dashed line corresponds to the 3U cubesat size.
Sensors 17 00493 g002
Figure 3. Payload concept.
Figure 3. Payload concept.
Sensors 17 00493 g003
Figure 4. The four-facet round-base pyramid concept.
Figure 4. The four-facet round-base pyramid concept.
Sensors 17 00493 g004
Figure 5. Principle of measurement: (a) shows the ideal spot image of a star, as seen through the payload without the pyramid; (b) shows the ideal spot image of a star as seen through the payload with the pyramid; the four spots are dispersed and the centroids (white circles) are joined by white lines to show their relative distances; (c) shows the case of image (b) during the transit of a planet, a uniform decrement of luminosity determines no shift of centroids; (d) shows the case of image (b) during a false positive transit, where the non-uniform decrement of luminosity determines the shift of centroids and then the change of their relative distances.
Figure 5. Principle of measurement: (a) shows the ideal spot image of a star, as seen through the payload without the pyramid; (b) shows the ideal spot image of a star as seen through the payload with the pyramid; the four spots are dispersed and the centroids (white circles) are joined by white lines to show their relative distances; (c) shows the case of image (b) during the transit of a planet, a uniform decrement of luminosity determines no shift of centroids; (d) shows the case of image (b) during a false positive transit, where the non-uniform decrement of luminosity determines the shift of centroids and then the change of their relative distances.
Sensors 17 00493 g005
Figure 6. First pyramid design in zemax. The red line is the equivalent lens corresponding to the objective.
Figure 6. First pyramid design in zemax. The red line is the equivalent lens corresponding to the objective.
Sensors 17 00493 g006
Figure 7. Integration time as function of target star V-magnitude. The SNR is 105. The red line is the maximum integration time limit, i.e., 4 s.
Figure 7. Integration time as function of target star V-magnitude. The SNR is 105. The red line is the maximum integration time limit, i.e., 4 s.
Sensors 17 00493 g007
Table 1. Matrix of subsystem requirements.
Table 1. Matrix of subsystem requirements.
SubsystemR.1R.2R.3R.4
Cubesat Payload SubsystemThe signal amplitude is 84 ppm, the tolerable noise level is 12 ppm in a duty cycle time observation. The detector resolution should be greater than or equal to 12 bit.The payload is a photometric payload (at least objective + detector) in the visible band.The transit measurement signal must be redundant and spectrally dispersed.The optical aperture diameter should be less than 10 cm.
Cubesat Attitude determination and ControlHigh-stability pointing should be provided to detect the signal decrement (5″). Total stored chemical energy will not exceed 100 Watt-hours.
Cubesat Command and data handlingThe onboard processing should be as limited as possible, to avoid any data alteration. Select the scientific data to be stored according to the available data rate.Total stored chemical energy will not exceed 100 Watt-hours.
Cubesat Structure and Mechanism No more than one unit free space is available for the payload subsystem. Platform dimensions are 10 cm × 10 cm × 34 cm. Platform maximum mass is 4 kg.
Mission OperationsThe minimum mission lifetime is four years, to measure the transit of an exoplanet with an orbital period of one year at least three times.The scientific operations consists in looking at one target star continuously to measure the star flux along the entire orbital period of the exoplanet.
Table 2. Sun and Earth radius, expected transit depth of an Earth-like planet around a Sun-like star.
Table 2. Sun and Earth radius, expected transit depth of an Earth-like planet around a Sun-like star.
ParameterValue
Rstar, Sun volumetric mean radius (km)695,700
Rplanet, Earth volumetric mean radius (km)6371
δ, transit depth 84 × 10−6
Table 3. List of main characteristics of HAS2 star tracker and ZEISS Planar T 1.4/85 objective.
Table 3. List of main characteristics of HAS2 star tracker and ZEISS Planar T 1.4/85 objective.
HAS2ZEISS Objective
ParameterValueParameterValue
Overall dimensions 30 mm × 30 mm f85 mm
Image sensor format1024 × 1024 pixelsF/#1.4
Pixel size 18 μmDiameter60.7 mm
ADC resolution 12 bit Image diameter43 mm
Saturation voltage output 1.49 V FOV12.42 deg
Full well capacity 105 ePixel scale44.15 arcsec/pix
Quantum Efficiency 45% (500–650 nm)
Spectral response 33.3% (400–900 nm )
Conversion factor 14.8 µV/e
Dark current 12.5 e/pix/s
RON 2 e/pix
Table 4. List of payload components and budget of dimensions and weight for the payload unit.
Table 4. List of payload components and budget of dimensions and weight for the payload unit.
ComponentNameLength X (mm)Length Y (mm)Length Z (mm)Weight (g)
Pyramid custom, BK7 glass626213.885.38
ObjectivePlanar T 1.4/85 ZF777762570
DetectorHAS230304.58
Detectore2V101025
Imager boardcustom PCB 3
StagePI P-733.2 CL10010025580
Stage controllercustom (on PCB) 3
Check/Tot each ≤ 100each ≤ 100107.31254.38
Table 5. Constraints, requirements, and assumptions for the limiting V-magnitude computation.
Table 5. Constraints, requirements, and assumptions for the limiting V-magnitude computation.
ConstraintsRequirementsAssumptions
Duty cycle time [min]31S/N ≥105window [pix]40 × 40
Photons to saturate [phot/pix]2.237 × 105Photons/pix108system efficiency70%
star defocus [pix]10 × 10
Table 6. Computation steps to determine the limiting V-magnitude.
Table 6. Computation steps to determine the limiting V-magnitude.
Computation Steps
StepStep Title Step DescriptionValues
1Compute the signal level (see the text above)
2Compute the noise level(see the text above)
3Compute the S/N ratio(see the text above)
4Min number of frames Photon per pixel/Photon to saturate108/(2.237 × 105) = 447 frames
5Max exposure time per frameEclipse time/Min number of frames4 s
6Exposure time limit (tlim) to avoid saturation Compute the detector voltage corresponding to the incoming photon flux for one frame.
Find the exposure time to get the higher voltage value without saturation. The star flux is spread on 10 × 10 pix.
vmag = 0, tlim = 0.17 s;
vmag = 1, tlim = 0.43 s;
vmag = 2, tlim = 1.09 s;
vmag = 3, tlim = 2.74 s;
vmag = 4, tlim = 6.89 s;
vmag = 5, tlim > 10 s;
7Number of co-added frames Compute the number of co-added frames required to get a S/N of 105 (star signal of 1010).vmag = 0, frames = 451;
vmag = 1, frames = 451;
vmag = 2, frames = 451;
vmag = 3, frames = 453;
vmag = 4, frames = 461;
vmag = 5, frames = 823 (10 s);
8Verify the compliance with the duty cycle time constraint. t = frames × tlimvmag = 0, t = 77.75 s;
vmag = 1, t = 194.94 s;
vmag = 2, t = 492.43 s;
vmag = 3, t = 1243 s;
vmag = 4, t = 3175 s;
vmag = 5, t = 8236 s.

Share and Cite

MDPI and ACS Style

Iuzzolino, M.; Accardo, D.; Rufino, G.; Oliva, E.; Tozzi, A.; Schipani, P. A Cubesat Payload for Exoplanet Detection. Sensors 2017, 17, 493. https://doi.org/10.3390/s17030493

AMA Style

Iuzzolino M, Accardo D, Rufino G, Oliva E, Tozzi A, Schipani P. A Cubesat Payload for Exoplanet Detection. Sensors. 2017; 17(3):493. https://doi.org/10.3390/s17030493

Chicago/Turabian Style

Iuzzolino, Marcella, Domenico Accardo, Giancarlo Rufino, Ernesto Oliva, Andrea Tozzi, and Pietro Schipani. 2017. "A Cubesat Payload for Exoplanet Detection" Sensors 17, no. 3: 493. https://doi.org/10.3390/s17030493

APA Style

Iuzzolino, M., Accardo, D., Rufino, G., Oliva, E., Tozzi, A., & Schipani, P. (2017). A Cubesat Payload for Exoplanet Detection. Sensors, 17(3), 493. https://doi.org/10.3390/s17030493

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop