High-Resolution Infrared Imaging of Neptune from the Keck Telescope

We present results of infrared observations of Neptune from the 10-m W. M. Keck I Telescope, using both high-resolution (0.04 arcsecond) broadband speckle imaging and conventional imaging with narrowband filters (0.6 arcsec resolution). The speckle data enable us to track the size and shape of infrared-bright features (“storms”) as they move across the disk and to determine rotation periods for latitudes −30 and −45°. The narrowband data are input to a model that allows us to make estimates of Neptune's stratospheric haze abundance and the size of storm features. We find a haze column density of ∼10⁶ cm⁻² for a haze layer located in the stratosphere, and a lower limit of 10⁷ cm⁻² and an upper limit of 10⁹ cm⁻² for a layer of 0.2 μm particles in the troposphere. We also calculate a lower limit of 7×10⁶ km² for the size of a “storm” feature observed on 13 October 1997.     

A reanalysis of the Apollo light scattering observations, and implications for lunar exospheric dust

Conspicuous excess brightness, exceeding that expected from coronal and zodiacal light (CZL), was observed above the lunar horizon in the Apollo 15 coronal photographic sequence acquired immediately after orbital sunset (surface sunrise). This excess brightness systematically faded as the Command Module moved farther into shadow, eventually becoming indistinguishable from the CZL background. These observations have previously been attributed to scattering by ultrafine dust grains (radius ∼0.1 microns) in the lunar exosphere, and used to obtain coarse estimates of dust concentration at several altitudes and an order-of-magnitude estimate of ∼10⁻⁹ g cm⁻² for the column mass of dust near the terminator, collectively referred to as model “0”.We have reanalyzed the Apollo 15 orbital sunset sequence by incorporating the known sightline geometries in a Mie-scattering simulation code, and then inverting the measured intensities to retrieve exospheric dust concentration as a function of altitude and distance from the terminator. Results are presented in terms of monodisperse (single grain size) dust distributions. For a grain radius of 0.10 microns, our retrieved dust concentration near the terminator (∼0.010 cm⁻³) is in agreement with model “0” at z=10 km, as is the dust column mass (∼3–6×10⁻¹⁰ g cm⁻²), but the present results indicate generally larger dust scale heights, and much lower concentrations near 1 km (<0.08 cm⁻³ vs. a few times 0.1 cm⁻³ for model “0"). The concentration of dust at high altitudes (z>50 km) is virtually unconstrained by the measurements. The dust exosphere extends into shadow a distance somewhere between 100 and 200 km from the terminator, depending on the uncertain contribution of CZL to the total brightness. These refined estimates of the distribution and concentration of exospheric dust above the lunar sunrise terminator should place new and more rigorous constraints on exospheric dust transport models, as well as provide valuable support for upcoming missions such as the Lunar Atmosphere and Dust Environment Explorer (LADEE).     

Improved Determination of Ethane (C2H6) Abundance in Titan's Stratosphere

Ethane spectral lines were observed in emission from Titan in August 1993, October 1995, and September 1996, at a spectral resolution of λ/Δλ≈10⁶, at wavelength 11.7−11.9 μm using the Goddard Infrared Heterodyne Spectrometer at the NASA Infrared Telescope Facility on Mauna Kea, Hawaii. The ethane mole fraction is determined to be (8.8±2.2)×10⁻⁶ (68.3% confidence limits, “1σ”), averaging the retrievals from each observing run obtained using the “recommended” thermal profile of R. V. Yelle, D. Strobel, E. Lellouch, and D. Gautier (1997, in Huygens: Science, Payload, and Mission (J.-P. Lebreton, Ed.), pp. 243-256, European Space Agency SP-1177).     

Global Oscillations at Low Frequency from the SOHO mission (GOLF)

The GOLF experiment on the SOHO mission aims to study the internal structure of the sun by measuring the spectrum of global oscillations in the frequency range 10^–7 to 10^–2 Hz. Bothp andg mode oscillations will be investigated, with the emphasis on the low order long period waves which penetrate the solar core. The instrument employs an extension to space of the proven ground-based technique for measuring the mean line-of-sight velocity of the viewed solar surface. By avoiding the atmospheric disturbances experienced from the ground, and choosing a non-eclipsing orbit, GOLF aims to improve the instrumental sensitivity limit by an order of magnitude to 1 mm s^–1 over 20 days for frequencies higher than 2.10^–4 Hz. A sodium vapour resonance cell is used in a longitudinal magnetic field to sample the two wings of the solar absorption line. The addition of a small modulating field component enables the slope of the wings to be measured. This provides not only an internal calibration of the instrument sensitivity, but also offers a further possibility to recognise, and correct for, the solar background signal produced by the effects of solar magnetically active regions. The use of an additional rotating polariser enables measurement of the mean solar line-of-sight magnetic field, as a secondary objective.     

Decreased values of cosmic dust number density estimates in the Solar System

Experiments to investigate the effect of impacts on side-walls of dust detectors such as the present NASA/ESA Galileo/Ulysses instrument are reported. Side walls constitute 27% of the internal area of these instruments, and increase field of view from 140° to 180°. Impact of cosmic dust particles onto Galileo/Ulysses Al side walls was simulated by firing Fe particles, 0.5-5 μm diameter, 2-50 km s⁻¹, onto an Al plate, simulating the targets of Galileo and Ulysses dust instruments. Since side wall impacts affect the rise time of the target ionization signal, the degree to which particle fluxes are overestimated varies with velocity. Side-wall impacts at particle velocities of 2-20 km s⁻¹ yield rise times 10-30% longer than for direct impacts, so that derived impact velocity is reduced by a factor of ∼2. Impacts on side wall at 20-50 km s⁻¹ reduced rise times by a factor of ∼10 relative to direct impact data. This would result in serious overestimates of flux of particles intersecting the dust instrument at velocities of 20-50 km s⁻¹. Taking into account differences in laboratory calibration geometry we obtain the following percentages for previous overestimates of incident particle number density values from the Galileo instrument [Grün et al., 1992. The Galileo dust detector. Space Sci. Rev. 60, 317-340]: 55% for 2 km s⁻¹ impacts, 27% at 10 km s⁻¹ and 400% at 70 km s⁻¹. We predict that individual particle masses are overestimated by ∼10-90% when side-wall impacts occur at 2-20 km s⁻¹, and underestimated by ∼10-¹⁰² at 20-50 km s⁻¹. We predict that wall impacts at 20-50 km s⁻¹ can be identified in Galileo instrument data on account of their unusually short target rise times. The side-wall calibration is used to obtain new revised values [Krüger et al., 2000. A dust cloud of Ganymede maintained by hypervelocity impacts of interplanetary micrometeoroids. Planet. Space Sci. 48, 1457-1471; 2003. Impact-generated dust clouds surrounding the Galilean moons. Icarus 164, 170-187] of the Galilean satellite dust number densities of 9.4×¹⁰⁻⁵, 9.9×¹⁰⁻⁵, 4.1×¹⁰⁻⁵, and 6.8×¹⁰⁻⁵ m⁻³ at 1 satellite radius from Io, Europa, Ganymede, and Callisto, respectively. Additionally, interplanetary particle number densities detected by the Galileo mission are found to be 1.6×¹⁰⁻⁴, 7.9×¹⁰⁻⁴, 3.2×¹⁰⁻⁵, 3.2×¹⁰⁻⁵, and 7.9×¹⁰⁻⁴ m⁻³ at heliocentric distances of 0.7, 1, 2, 3, and 5 AU, respectively. Work by Burchell et al. [1999b. Acceleration of conducting polymer-coated latex particles as projectiles in hypervelocity impact experiments. J. Phys. D: Appl. Phys. 32, 1719-1728] suggests that low-density “fluffy” particles encountered by Ulysses will not significantly affect our results—further calibration would be useful to confirm this.     

Temporally varying ethylene emission on Jupiter

Ethylene (C₂H₄) emission has been measured in the poles and equator of Jupiter. The 949 cm⁻¹ spectra were recorded with a high resolution spectrometer at the McMath-Pierce telescope at Kitt Peak in October-November 1998 and at the Infrared Telescope Facility at Mauna Kea in June 2000. C₂H₄ is an important product of methane chemistry in the outer planets. Knowledge of its abundance can help discriminate among the various proposed sets of CH₄ photolysis branching ratios at Ly-α, and determine the relative importance of the reaction pathways that produce C₂H₂ and C₂H₆. In the equatorial region the C₂H₄ emission is weak, and we were only able to detect it at high air-mass, near the limb. We derive a peak equatorial molar abundance of C₂H₄ of 4.5×¹⁰⁻⁷-1.7×¹⁰⁻⁶ near 2.2×¹⁰⁻³ mbar, with a total column of 5.7×¹⁰¹⁴-2.2×¹⁰¹⁵ molecules cm⁻² above 10 mbar depending upon choice of thermal profile. We observed enhanced C₂H₄ emission from the poles in the regions where auroras are seen in X-ray, UV, and near infrared images. In 2000 we measured a short-term change in the distribution of polar C₂H₄ emission; the emission in the north IR auroral “hot spot” decreased by a factor of three over a two-day interval. This transient behavior and the sensitivity of C₂H₄ emission to temperature changes near its contribution peak at 5-10 microbar suggests that the polar enhancement is primarily a thermal effect coupled with vertical transport. Comparing our observations from Kitt Peak and Mauna Kea shows that the C₂H₄ emission of the northern non-“hot spot” auroral regions did not change over the three-year period while that in the southern polar regions decreased.     

The DARWIN project

Angel and co-workers have proposed to detect exoplanets around nearby stars in the infrared (6–17 m) and to analyze their spectra, searching for H₂O, CO₂, CH₄, NH₃, and O₃ spectral features. The presence or absence of CO₂ would indicate either a strong similarity or difference with Solar telluric planet atmospheres. Water would indicate a habitable planet, and O₃ would reveal significant photosynthesis activity, due to the presence of carbon chemistry based life. Like these authors, we suggest an infrared nulling interferometer pointing to the star and working as a coronograph. Our main contribution is to propose an observatory made of four to five 1-meter class telescopes observing from about 4 to 5 AU to avoid the Solar Zodiacal Light (ZL) background at 10m instead of four 8-meter ones observing from the Earth vicinity. This allows the mission to be feasible in thenear future. The concept, named DARWIN, is under consideration by the European Space Agency for its Horizon 2000 Plus program.     

The Dual Sources of Io's Sodium Clouds

Io's sodium clouds result mostly from a combination of two atmospheric escape processes at Io. Neutralization of Na⁺ and/or NaX⁺ pickup ions produces the “stream” and the “jet” and results in a rectangular-shaped sodium nebula around Jupiter. Atmospheric sputtering of Na by plasma torus ions produces the “banana cloud” near Io and a diamond-shaped sodium nebula. Charge exchange of thermal Na⁺ with Na in Io's atmosphere does not appear to be a major atmospheric ejection process. The total ejection rate of sodium from Io varied from 3×10²⁶ to 25×10²⁶ atoms/s over seven years of observations. Our results provide further evidence that Io's atmospheric escape is driven from collisionally thick regions of the atmosphere rather than from the exosphere.     

Optimal dynamos in the cores of terrestrial exoplanets: Magnetic field generation and detectability

A potentially promising way to gain knowledge about the internal dynamics of extrasolar planets is by remote measurement of an intrinsic magnetic field. Strong planetary magnetic fields, maintained by internal dynamo action in an electrically conducting fluid layer, are helpful for shielding the upper atmosphere from stellar wind induced mass loss and retaining water over long (Gyr) time scales. Here we present a whole planet dynamo model that consists of three main components: an internal structure model with composition and layers similar to the Earth, an optimal mantle convection model that is designed to maximize the heat flow available to drive convective dynamo action in the core, and a scaling law to estimate the magnetic field intensity at the surface of a terrestrial exoplanet. We find that the magnetic field intensity at the core surface can be up to twice the present-day geomagnetic field intensity, while the magnetic moment varies by a factor of 20 over the models considered. Assuming electron cyclotron emission is produced from the interaction between the stellar wind and the exoplanet magnetic field we estimate the cyclotron frequencies around the ionospheric cutoff at 10 MHz with emission fluxes in the range 10⁻⁴-10⁻⁷ Jy, below the current detection threshold of radio telescopes. However, we propose that anomalous boosts and modulations to the magnetic field intensity and cyclotron emission may allow for their detection in the future.     

Estimating the significance of a signal in a multi-dimensional search

In experiments that are aimed at detecting astrophysical sources such as neutrino telescopes, one usually performs a search over a continuous parameter space (e.g. the angular coordinates of the sky, and possibly time), looking for the most significant deviation from the background hypothesis. Such a procedure inherently involves a “look elsewhere effect”, namely, the possibility for a signal-like fluctuation to appear anywhere within the search range. Correctly estimating the p-value of a given observation thus requires repeated simulations of the entire search, a procedure that may be prohibitively expansive in terms of CPU resources. Recent results from the theory of random fields provide powerful tools which may be used to alleviate this difficulty, in a wide range of applications. We review those results and discuss their implementation, with a detailed example applied for neutrino point source analysis in the IceCube experiment.     

Two-frequency timing of the pulsar B1937+21 in Kalyazin and Kashima in 1997–2002

We present the results from our timing of the millisecond pulsar B1937+21, performed jointly since 1997 on two radio telescopes: the RT-64 in Kalyazin (Russia) at a frequency of 0.6GHz and RT-34 in Kashima (Japan) at a frequency of 2.15 GHz. The rms value of the pulse time of arrival (TOA) residuals for the pulsar at the barycenter of the Solar system is 1.8 μs (the relative variation is ≈10^?14 over the observing period). The TOA residuals are shown to be dominated by white phase noise, which allows this pulsar to be used as an independent time scale keeper. The upper limit for the gravitational background energy density Ω_gh² at frequencies ≈6.5 × 10^?9 Hz is estimated to be no higher than 10^?6. Based on the long-term timing of the pulsar, we have improved its parameters and accurately determined the dispersion measure and its time variation over the period 1984–2002, which was, on average, ?0.00114(3) pc cm^?3 yr^?1.     

Fluctuations of energetic particle flux during solar cycle based on measurements in the solar wind,in the magnetosphere,and at Earth

We present the results of our studies of the cosmic-ray fluctuations in the frequency range 10⁻⁴−1.67 × 10⁻³ Hz based on energetic particle flux measurements on spacecraft in the solar wind, in the magnetosphere, and at Earth in the 11-year solar cycle. The cosmic-ray fluctuation spectrum is shown to have an 11-year modulation related to the solar cycle. A different behavior of the level of energetic particle fluctuations measured in different regions of space is observed for cosmic rays of different origins. We conclude that the new, previously unknown phenomenon of 11-year modulation of the cosmic-ray fluctuation spectrum has been established. A possible explanation of this phenomenon is given.     

Neutrino Dominated Accretion Flow Around a Strange Star

The “strange star - NDAF” model (NDAF: Neutrino Dominated Accretion Flow) is proposed as an alternative central engine of gamma-ray bursts for unifying the interpretation of the prompt emission and postburst activities of gamma-ray bursts. The structure of NDAF around a strange star is calculated. Different from other central compact objects, the strange star will feed back the phase transition energy of strangization on the accretion flow, with neutrinos as energy carriers. The friction between NDAF and strange star is ignored in this paper. The results indicate: firstly, the structure of NDAF around a strange star is sensitive to accretion rate; secondly, if accretion rate is larger than 0.18 M_? s^-1, the “strange star - NDAF” model can unify the explanation on the prompt emission and postburst activities of gamma-ray bursts, and the range of allowable accretion rates is wider than that in frictionless “neutron star - NDAF” models; thirdly, the range of annihilation energy of “strange star - NDAF” model is very wide, when the accretion rate is higher than 0.3 M_? s^-1, the annihilation energy is greater than 10⁵¹ erg; finally, if the accretion rate is greater than 0.3 M_? s^-1, the annihilation energy of “strange star - NDAF” model is larger than what of “black hole - NDAF” model at the same accretion rate by more than one order of magnitude, it is favorable to explaining some extremely energetic gamma-ray bursts.     

MESSENGER observations of the plasma environment near Mercury

The MESSENGER Fast Imaging Plasma Spectrometer (FIPS) measured the bulk plasma characteristics of Mercury's magnetosphere and solar wind environment during the spacecraft's first two flybys of the planet on 14 January 2008 (M1) and 6 October 2008 (M2), producing the first measurements of thermal ions in Mercury's magnetosphere. In this work, we identify major features of the Mercury magnetosphere in the FIPS proton data and describe the data analysis process used for recovery of proton density (n_p) and temperature (T_p) with a forward modeling technique, required because of limitations in measurement geometry. We focus on three regions where the magnetospheric flow speed is likely to be low and meets our criteria for the recovery process: the M1 plasma sheet and the M1 and M2 dayside and nightside boundary-layer regions. Interplanetary magnetic field (IMF) conditions were substantially different between the two flybys, with intense reconnection signatures observed by the Magnetometer during M2 versus a relatively quiet magnetosphere during M1. The recovered ion density and temperature values for the M1 quiet-time plasma sheet yielded n_p∼1–10 cm⁻³, T_p∼2×10⁶ K, and plasma β∼2. The nightside boundary-layer proton densities during M1 and M2 were similar, at n_p∼4–5 cm⁻³, but the temperature during M1 (T_p∼4–8×10⁶ K) was 50% less than during M2 (T_p∼8×10⁶ K), presumably due to reconnection in the tail. The dayside boundary layer observed during M1 had a density of ∼16 cm⁻³ and temperature of 2×10⁶ K, whereas during M2 this region was less dense and hotter (n_p∼8 cm⁻³ and T_p∼10×10⁶ K), again, most likely due to magnetopause reconnection. Overall, the southward interplanetary magnetic field during M2 clearly produced higher T_p in the dayside and nightside magnetosphere, as well as higher plasma β in the nightside boundary, ∼20 during M2 compared with ∼2 during M1. The proton plasma pressure accounts for only a fraction (24% for M1 and 64% for M2) of the drop in magnetic pressure upon entry into the dayside boundary layer. This result suggests that heavy ions of planetary origin, not considered in this analysis, may provide the “missing” pressure. If these planetary ions were hot due to “pickup” in the magnetosheath, the required density for pressure balance would be an ion density of ∼1 cm⁻³ for an ion temperature of ∼10⁸ K.     

Changes in Jupiter’s zonal velocity between 1979 and 2008

We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s⁻¹) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s⁻¹ on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s⁻¹, so velocity fluctuations of 10 m s⁻¹ may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s⁻¹ in 2000 and ∼130 m s⁻¹ in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s⁻¹. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s⁻¹ at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s⁻¹ rather than ∼105 m s⁻¹, then the vertical zonal wind shear is much less than the currently accepted value.     

Under the shadow of the Magellanic Bridge: A measurement of the extragalactic background at ∼ 0.7 keV

We measured the extragalactic 0.7 keV X-ray background by observing the X-ray shadow of a neutral gas cloud in the Magellanic Bridge region. Two ROSAT PSPC observations of total 104 ks were complemented by a detailed H I mapping of the cloud with both the Parkes 64 m telescope and the Australia Telescope Compact Array. From the detected anti-correlation between the observed background intensity and the H I column density of the cloud, we derived the unabsorbed extragalactic background intensity as ∼ 28 keV s⁻¹ cm⁻² keV⁻¹ sr⁻¹ at ∼ 0.7 keV. The 95% confidence lower limit 18 keV s⁻¹ cm⁻² keV⁻¹ sr⁻¹ is greater than the expected point-like source contribution ? 14 keV s⁻¹ cm⁻² keV⁻¹ sr⁻¹, constrained by the mean source spectrum together with the total background intensity in the 1-2 keV band. A significant fraction of the 0.7 keV background likely arises in a diffuse hot intergalactic medium of a few million degrees, as has been predicted in hydrodynamic simulations of cosmological structure formation.Richard McCray     

Radon induced background processes in the KATRIN pre-spectrometer

The KArlsruhe TRItium Neutrino (KATRIN) experiment is a next generation, model independent, large scale tritium β-decay experiment to determine the effective electron anti-neutrino mass by investigating the kinematics of tritium β-decay with a sensitivity of 200 meV/c² using the MAC-E filter technique. In order to reach this sensitivity, a low background level of 10⁻² counts per second (cps) is required. This paper describes how the decay of radon in a MAC-E filter generates background events, based on measurements performed at the KATRIN pre-spectrometer test setup. Radon (Rn) atoms, which emanate from materials inside the vacuum region of the KATRIN spectrometers, are able to penetrate deep into the magnetic flux tube so that the α-decay of Rn contributes to the background. Of particular importance are electrons emitted in processes accompanying the Rn α-decay, such as shake-off, internal conversion of excited levels in the Rn daughter atoms and Auger electrons. While low-energy electrons (<100 eV) directly contribute to the background in the signal region, higher energy electrons can be stored magnetically inside the volume of the spectrometer. Depending on their initial energy, they are able to create thousands of secondary electrons via subsequent ionization processes with residual gas molecules and, since the detector is not able to distinguish these secondary electrons from the signal electrons, an increased background rate over an extended period of time is generated.     

Spectroscopic Observations of the Star Formation Regions in Nearby Galaxies

In recent years the number of worldwide 8∼10 m-class ground-based telescopes is continually increased, the 4 m-diameter or smaller telescopes have become the small and medium-sized telescopes. In order to obtain some noticeable scientific results by using these existing small and medium-sized telescopes, we have to consider very carefully what we can do, and what we can not. For this reason, the Time Allocation Committee of the 2.16 m telescope of the National Astronomical observatories of China (NAOC) has decided to support some key projects since 2013. The long-term project “Spectroscopic Observations of the Star Formation Regions in Nearby Galaxies” proposed by us is one of three key projects, it is supported by the committee with 30 dark/grey nights in each of three years.     

Making other earths: dynamical simulations of terrestrial planet formation and water delivery

We present results from 44 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary “embryos” (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 44 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 “habitable” planets between 0.9 and 1.1 AU. These planets range from dry worlds to “water worlds” with 100+oceans of water (1 ocean=1.5×10²⁴ g), and vary in mass between 0.23M_⊕ and 3.85M_⊕. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call “super embryos,” 1-2M_⊕ protoplanets which can serve as the accretion seeds for 2+M_⊕ planets with large water contents.     

Dark red debris from three short-period comets: 2P/Encke, 22P/Kopff, and 65P/Gunn

We present observations of the extended dust structures near the orbits of three short-period comets: 2P/Encke, 22P/Kopff, and 65P/Gunn. The dust trails were originally discovered by the Infrared Astronomical Satellite (IRAS). Our observations were made using wide-field optical CCD cameras on the University of Hawaii 2.24-m telescope, the Canada-France-Hawaii 3.6-m telescope, and the Kiso 1.05-m Schmidt telescope. We compared the observed images with models and found that the extended structures seen around 2P/Encke and 22P/Kopff before perihelion passage were most likely “dust trails,” whereas images taken after perihelion passage show a high contamination by recently released particles (i.e., particles in Neck-Line structures are visible). We could not confirm the existence of a dust trail from 65P/Gunn within the field of view of the camera used. The effective sizes of the particles responsible for the scattered light were estimated at 1-100 mm (2P/Encke), 1-10 mm (22P/Kopff), and 100 μm-1 mm (65P/Gunn), respectively, which is consistent with previous studies of dust trails made with infrared space telescopes and optical telescopes. We evaluated the mass loss rates of these comets, averaged over their orbits, as reaching (2P/Encke), (22P/Kopff), and (65P/Gunn). These values are consistent with previous work. Therefore, the total amount of material ejected from these three comets is , which would contribute a considerable fraction of the lost within 1 AU that needs to be replaced if the zodiacal cloud is to be maintained in a steady state. We also found that the particles in the dust structures are significantly redder than the Sun and the zodiacal light, and might be redder than the average short-period comet nuclei. Specifically, the reflectivity gradients of 2P/Encke, 22P/Kopff, and 65P/Gunn are 13±7 (% ¹⁰³ Å⁻¹), 20±5 (% ¹⁰³ Å⁻¹), and 15±4 (% ¹⁰³ Å⁻¹), respectively. We examined the change in color with distance from the nucleus. No clear correlation was detected for 2P/Encke or 22P/Kopff to an accuracy of 3-11%, while the 65P/Gunn tail did show color variation, becoming redder with increasing distance from the nucleus. This dark red material, consisting of particles of sand-cobble size, has marginally escaped from the nuclei and will evolve into finer-grained interplanetary dust particles after subsequent collisions.