Exoplanet Code Download

Summer 2007:

MIT Department of Earth, Atmospheric and Planetary Science

Summer UROP with Prof. Sara Seager on Mass-Radius relation of Earth-like exoplanet.

★We Built Computer Models for the interior structure and the atmosphere of extra-solar planets to understand the Mass-Radius relation and atmospheric spectral lines of exoplanets. Our basic assumption for the interior of Earth-like exoplanets is that they all have an iron core, a silicate mantle and a water crust. We also assumed that the temperature dependence of the Equation Of State (EOS) is negligiable, which is a very good approximation.

Based on these assumptions, we developed two computer codes in matlab to interpret the bulk composition of solid exoplanets based on their mass and radius measurements.

For details, please refer to our paper published on Publications of the Astronomical Society of the Pacific (PASP).

Code download instructions:

Frist of all, download the code zip files and decompress them.

Please put the decompressed files into your Matlab work directory so that the Matlab can recognize them. Use addpath command if necessary.

Based in Matlab if using this computer code please cite Li Zeng & Prof. Sara Seager.

ExoterDE(M, Munc, R, Runc); Take a lot of time, but more accurate. This program can deal with the a case of zero uncertainties in planet mass and radius. It plots 1-σ,2-σ,3-σ contours of given Mass, Radius and the uncertainties associated with them.

example: for Mass=10 Earth Mass , Mass uncertainty=0.5 Earth Mass;

Radius= 2 Earth Radius, Radius uncertainty= 0.1 Earth Radius:

just type ExoterDE(10, 0.5, 2, 0.1); in your Matlab command window. Then hit Enter.

ExoterDB(M, Munc, R, Runc); based upon database interpolation. It plots a colormap showing the possible proportions of iron, silicate and water for continuous range of σ from 0 up to 3.

example: for Mass=10 Earth Mass , Mass uncertainty=0.5 Earth Mass;

Radius= 2 Earth Radius, Radius uncertainty= 0.1 Earth Radius:

just type ExoterDB(10, 0.5, 2, 0.1); in your Matlab command window. Then hit Enter.

Home