Multi-wavelength simulations of atmospheric radiation from Io with a 3-D spherical-shell backward Monte Carlo radiative transfer model

Conflicting observations regarding the dominance of either sublimation or volcanism as the source of the atmosphere on Io and disparate reports on the extent of its spatial distribution and the absolute column abundance invite the development of detailed computational models capable of improving our understanding of Io’s unique atmospheric structure and origin. Improving upon previous models, Walker et al. (Walker, A.C., Gratiy, S.L., Levin, D.A., Goldstein, D.B., Varghese, P.L., Trafton, L.M., Moore, C.H., Stewart, B. [2009]. Icarus) developed a fully 3-D global rarefied gas dynamics model of Io’s atmosphere including both sublimation and volcanic sources of SO₂ gas. The fidelity of the model is tested by simulating remote observations at selected wavelength bands and comparing them to the corresponding astronomical observations of Io’s atmosphere. The simulations are performed with a new 3-D spherical-shell radiative transfer code utilizing a backward Monte Carlo method. We present: (1) simulations of the mid-infrared disk-integrated spectra of Io’s sunlit hemisphere at 19 μm, obtained with TEXES during 2001-2004; (2) simulations of disk-resolved images at Lyman-α obtained with the Hubble Space Telescope (HST), Space Telescope Imaging Spectrograph (STIS) during 1997-2001; and (3) disk-integrated simulations of emission line profiles in the millimeter wavelength range obtained with the IRAM-30 m telescope in October-November 1999. We found that the atmospheric model generally reproduces the longitudinal variation in band depth from the mid-infrared data; however, the best match is obtained when our simulation results are shifted ∼30° toward lower orbital longitudes. The simulations of Lyman-α images do not reproduce the mid-to-high latitude bright patches seen in the observations, suggesting that the model atmosphere sustains columns that are too high at those latitudes. The simulations of emission line profiles in the millimeter spectral region support the hypothesis that the atmospheric dynamics favorably explains the observed line widths, which are too wide to be formed by thermal Doppler broadening alone.     

Microwave observations of Saturn's rings: anisotropy in directly transmitted and scattered saturnian thermal emission

We present a new Very Large Array (VLA) image of Saturn, made from data taken in October 1998 at a wavelength of λ3.6 cm. The moderate ring opening angle (B≈15°) allows us to explore direct transmission of microwave photons through the A and C rings. We find a strong asymmetry of photons transmitted through the A ring, but not in the C ring, a new diagnostic of wake structure in the ring particles. We also find a weak asymmetry between east and west for the far side of the ansae. To facilitate quantitative comparison between dynamic models of the A ring and radio observations, we extend our Monte Carlo radiative transfer code (described in Dunn et al., 2002, Icarus 160, 132-160) to include idealized wakes. We show the idealized model can reproduce the properties of dynamic simulations in directly transmitted light. We examine the model behavior in directly transmitted and scattered light over a range of physical and geometric wake parameters. Finally, we present a wake model with a plausible set of physical parameters that quantitatively reproduces the observed intensity and asymmetry of the A ring both across the planet and in the ansae.     

Monte Carlo simulations of radio pulses in atmospheric showers using ZHAireS

We present predictions for the radio pulses emitted by extensive air showers using ZHAireS, an AIRES-based Monte Carlo code that takes into account the full complexity of ultra-high energy cosmic-ray induced shower development in the atmosphere, and allows the calculation of the electric field in both the time and frequency domains. We do not presuppose any emission mechanism and our results are compatible with a superposition of geomagnetic and charge excess radio emission effects. We investigate the polarization of the electric field as well as the effects of the refractive index n and shower geometry on the radio pulses. We show that geometry, coupled to the relativistic effects that appear when using a realistic refractive index n > 1, play a prominent role on the radio emission of air showers.     

Multidimensional simulations of radiative transfer in Type Ia supernovae

A three-dimensional Monte Carlo code for modelling radiation transport in Type Ia supernovae is described. In addition to tracking Monte Carlo quanta to follow the emission, scattering and deposition of radiative energy, a scheme involving volume-based Monte Carlo estimators is used to allow properties of the emergent radiation field to be extracted for specific viewing angles in a multidimensional structure. This eliminates the need to compute spectra or light curves by angular binning of emergent quanta. The code is applied to two test problems to illustrate consequences of multidimensional structure on the modelling of light curves. First, elliptical models are used to quantify how large-scale asphericity can introduce angular dependence to light curves. Secondly, a model which incorporates complex structural inhomogeneity, as predicted by modern explosion models, is used to investigate how such structure may affect light-curve properties.     

Examining the wake structure in Saturn's rings from microwave observations over varying ring opening angles and wavelengths

Over the last 15 to 20 years several high quality, high resolution data have been taken with the very large array (VLA). These data exhibit a wide range of ring opening angles (|B|=0 to 26°) and wavelengths (λ=0.7 to 20 cm). At these wavelengths the primary flux from the rings is scattered saturnian thermal emission, with a small contribution coming from the ring particles' own thermal emission. Much of the data do show signs of asymmetries due to wakes either on the ansae or the portion of the rings which occult the planet. As in previous work, we use our Monte Carlo radiative transfer code including idealized wakes [Dunn, D.E., Molnar, L.A., Fix, J.D., 2002. Icarus 160, 132-160; Dunn, D.E., Molnar, L.A., Niehof, J.T., de Pater, I., Lissauer, J.L., 2004. Icarus 171, 183-198] to model the relative contributions of the scattered and thermal radiation emanating from the rings and compare the results to that seen in the data. Although the models do give satisfactory fits to all of our data, we find that no single model can simulate the data at all different |B| and λ. We find that one model works best for moderate and low |B| and another one at higher |B|. The main difference between these models is the ratio of the wake width to their separation. We similarly find that the 2 cm data require higher density wakes than the longer wavelength data, perhaps caused by a preponderance of somewhat smaller ring material in the wakes. We further find evidence for an increase in the physical temperature of the rings with increasing |B|. Continuous observations are required to determine whether the above results regarding variations in wake parameters with |B| and λ are indeed caused by these parameters, or instead by changes over time.     

Energetic Neutral Atoms (ENA) at Mars: Properties of the hydrogen atoms produced upstream of the martian bow shock and implications for ENA sounding technique around non-magnetized planets

We have studied the interaction of fast solar wind hydrogen atoms with the martian atmosphere by a three-dimensional Monte Carlo simulation. These energetic neutral hydrogen atoms, H-ENAs, are formed upstream of the martian bow shock. Both H-ENAs scattered and non-scattered from the martian atmosphere/exosphere were studied. The colliding H-ENAs were found to scatter both to the dayside and nightside. On the dayside they contribute to the so-called H-ENA albedo. On the nightside the heated and scattered hydrogen atoms were found also in the martian wake. The density, the energy distribution function and the direction of the velocity of H-ENAs on the nightside are presented. The present study describes a novel “ENA sounding” technique in which energetic neutral atoms are used to derive information of the properties of planetary exosphere and atmosphere in a similar manner as the solar wind photons are used to derive atmospheric densities by measuring the scattered UV light. A detailed study of the direction and energy of the scattered and non-scattered H-ENAs suggest that the ENA sounding is a method to study the interaction between the planetary atmosphere and the solar wind and to monitor the density, and likely also the magnetization, of the planetary upper atmosphere. Already present-day ENA instrument should be capable to detect the analyzed particle fluxes.     

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.     

The composition of liquid methane-nitrogen aerosols in Titan’s lower atmosphere from Monte Carlo simulations

Molecular level Monte Carlo simulations have been performed with various model potentials for the CH₄-N₂ vapor-liquid equilibrium at conditions prevalent in the atmosphere of Saturn’s moon Titan. With a single potential parameter adjustment to reproduce the vapor-liquid equilibrium at a higher temperature, Monte Carlo simulations are in excellent agreement with available laboratory measurements. The results demonstrate the ability of simple pair potential models to describe phase equilibria with the requisite accuracy for atmospheric modeling, while keeping the number of adjustable parameters at a minimum. This allows for stable extrapolation beyond the range of available laboratory measurements into the supercooled region of the phase diagram, so that Monte Carlo simulations can serve as a reference to validate phenomenological models commonly used in atmospheric modeling. This is most important when the relevant region of the phase diagram lies outside the range of laboratory measurements as in the case of Titan. The present Monte Carlo simulations confirm the validity of phenomenological thermodynamic equations of state specifically designed for application to Titan. The validity extends well into the supercooled region of the phase diagram. The possible range of saturation levels of Titan’s troposphere above altitudes of 7 km is found to be completely determined by the remaining uncertainty of the most recent revision of the Cassini-Huygens data, yielding a saturation of 100 ± 6% with respect to CH₄-N₂ condensation up to an altitude of about 20 km.     

A new look at photometry of the Moon

We use ROLO photometry (Kieffer, H.H., Stone, T.C. [2005]. Astron. J. 129, 2887-2901) to characterize the before and after full Moon radiance variation for a typical highlands site and a typical mare site. Focusing on the phase angle range 45° < α < 50°, we test two different physical models, macroscopic roughness and multiple scattering between regolith particles, for their ability to quantitatively reproduce the measured radiance difference. Our method for estimating the rms slope angle is unique and model-independent in the sense that the measured radiance factor I/F at small incidence angles (high Sun) is used as an estimate of I/F for zero roughness regolith. The roughness is determined from the change in I/F at larger incidence angles. We determine the roughness for 23 wavelengths from 350 to 939 nm. There is no significant wavelength dependence. The average rms slope angle is 22.2° ± 1.3° for the mare site and 34.1° ± 2.6° for the highland site. These large slopes, which are similar to previous “photometric roughness” estimates, require that sub-mm scale “micro-topography” dominates roughness measurements based on photometry, consistent with the conclusions of Helfenstein and Shepard (Helfenstein, P., Shepard, M.K. [1999]. Icarus 141, 107-131). We then tested an alternative and very different model for the before and after full Moon I/F variation: multiple scattering within a flat layer of realistic regolith particles. This model consists of a log normal size distribution of spheres that match the measured distribution of particles in a typical mature lunar soil 72141,1 (McKay, D.S., Fruland, R.M., Heiken, G.H. [1974]. Proc. Lunar Sci. Conf. 5, Geochim. Cosmochim. Acta 1 (5), 887-906). The model particles have a complex index of refraction 1.65-0.003i, where 1.65 is typical of impact-generated lunar glasses. Of the four model parameters, three were fixed at values determined from Apollo lunar soils: the mean radius and width of the log normal size distribution and the real part of the refraction index. We used FORTRAN programs from Mishchenko et al. (Mishchenko, M.I., Dlugach, J.M., Yanovitskij, E.G., Zakharova, N.T. [1999]. J. Quant. Spectrosc. Radiat. Trans. 63, 409-432; Mishchenko, M.I., Travis, L.D., Lacis, A.A. [2002]. Scattering, Absorption and Emission of Light by Small Particles. Cambridge Univ. Press, New York. <http://www.giss.nasa.gov/staff/mmishchenko/books.html>) to calculate the scattering matrix and solve the radiative transfer equation for I/F. The mean single scattering albedo is ω = 0.808, the asymmetry parameter is 〈cos Θ〉 = 0.77 and the phase function is very strongly peaked in both the forward and backward scattering directions. The fit to the observations for the highland site is excellent and multiply scattered photons contribute ?80% of I/F. We conclude that either model, roughness or multiple scattering, can match the observations, but that the strongly anisotropic phase functions of realistic particles require rigorous calculation of many orders of scattering or spurious photometric roughness estimates are guaranteed. Our multiple scattering calculation is the first to combine: (1) a regolith model matched to the measured particle size distribution and index of refraction of the lunar soil, (2) a rigorous calculation of the particle phase function and solution of the radiative transfer equation, and (3) application to lunar photometry with absolute radiance calibration.     

Scattering of Light by Composite Particles in a Planetary Surface

It has been proposed that composite particles containing internal scatterers may provide the explanation for the fact that most photometric studies of planetary surfaces based on Hapke's model of bidirectional reflectance have found the planetary particles to exhibit moderately backscattering phase functions. However, an implicit assumption made in this explanation is that the scattering by composite particles containing multiple internal inclusions in a planetary surface can still be adequately computed using standard radiative transfer theory assuming the composite particles to be the fundamental individual scatterers even though such particles are necessarily in close proximity to each other. In this paper, this assumption is explored by examining the effects of close packing on the light scattering by spherical particles containing isotropic internal scatterers using a Monte Carlo routine. As expected, classical radiative transfer (assuming a random distribution of scattering particles) coupled with the assumption that the composite particle is the fundamental scatterer provides a good approximation in the high porosity limit. However, even for porosities as high as 90% the effects of close packing are clearly seen with the radiative transfer calculation underestimating the scattering by ∼10% at high incidence, emission, and phase angles. As the porosity is lowered further, the discrepancy becomes more severe and can reach 50% or more. In contrast, assuming the individual scatterer properties in the radiative transfer calculation leads to a substantial overestimate of the scattering even for porosities as low as 27.5%. This suggests that parameters derived using the classical radiative transfer theory will yield results intermediate between those of the composite as a whole and those of the internal scatterers. Thus, one should exercise caution in interpreting the results of models based on classical radiative transfer theory in terms of the physical properties of the surface particles and, where possible, the bidirectional reflectance of densely packed composite particles should be computed using more accurate methods such as the stochastic radiative transfer theory.     

The Dust in Comet C/1999 S4 (LINEAR) during Its Disintegration: Narrow-Band Images, Color Maps, and Dynamical Models

Comet C/1999 S4 was observed with the 2m-telescopes of the Bulgarian National Observatory and Pik Terskol Observatory, Northern Caucasus, Russia, at the time of its disintegration. Maps of the dust brightness and color were constructed from images obtained in red and blue continuum windows, free from cometary molecular emissions. We analyze the dust environment of Comet C/1999 S4 (LINEAR) taking into account the observed changes apparent in the brightness images and in plots of Afρ profiles as function of the projected distance ρ from the nucleus. We also make use of the syndyne-synchrone formalism and of a Monte Carlo model based on the Finson-Probstein theory of dusty comets. The brightness and color of individual dust particles, which is needed to derive theoretical brightness and color maps of the cometary dust coma from the Monte Carlo model, is determined from calculations of the light scattering properties of randomly oriented oblate spheroids. In general, the dust of Comet C/1999 S4 (LINEAR) is strongly reddened, with reddening values up to 30%/1000 Å in some locations. Often the reddening is higher in envelopes further away from the nucleus. We observed two outbursts of the comet with brightness peaks on July 14 and just before July 24, 2000, when the final disintegration of the comet started. During both outbursts an excess of small particles was released. Shortly after both outbursts the dust coma “turns blue.” After the first outburst, the whole coma was affected; after the second one only a narrow band of reduced color close to the tail axis was formed. This difference is explained by different terminal ejection speeds, which were much lower than normal in case of the second outburst. In particular in the second, final outburst the excess small particles could originate from fragmentation of “fresh” larger particles.     

Concentric crater fill in Utopia Planitia: History and interaction between glacial “brain terrain” and periglacial mantle processes

At martian mid-to-high latitudes, the surfaces of potentially ice-rich features, including concentric crater fill, lobate debris aprons, and lineated valley fill, typically display a complex texture known as “brain terrain,” due to its resemblance to the complex patterns on brain surfaces. In order to determine the structure and developmental history of concentric crater fill and overlying latitude-dependent mantle (LDM) material, “brain terrain” and polygonally-patterned LDM surfaces are analyzed using HiRISE images from four craters in Utopia Planitia containing concentric crater fill. “Brain terrain” and mantle surface textures are classified based on morphological characteristics: (1) closed-cell “brain terrain,” (2) open-cell “brain terrain,” (3) high-center mantle polygons, and (4) low-center mantle polygons. A combined glacial and thermal-contraction cracking model is proposed for the formation and modification of the “brain terrain” texture of concentric crater fill. A similar model, related to thermal contraction cracking and differential sublimation of underlying ice, is proposed for the formation and development of polygonally patterned mantle material. Both models require atmospheric deposition of ice, likely during periods of high obliquity, but do not require wet active layer processes. Crater dating of “brain terrain” and mantled surfaces suggests a transition at martian mid-latitudes from peak “glacial” conditions occurring within the past ∼10-100 My to a quiescent period followed by a cold-desert “periglacial” period during the past ∼1-2 My.     

Optimising a balloon-borne polarimeter in the hard X-ray domain: From the PoGOLite Pathfinder to PoGO+

PoGOLite is a balloon-borne hard X-ray polarimeter dedicated to the study of point sources. Compton scattered events are registered using an array of plastic scintillator units to determine the polarisation of incident X-rays in the energy range 20–240 keV. In 2013, a near circumpolar balloon flight of 14 days duration was completed after launch from Esrange, Sweden, resulting in a measurement of the linear polarisation of the Crab emission. Building on the experience gained from this Pathfinder flight, the polarimeter is being modified to improve performance for a second flight in 2016. Such optimisations, based on Geant4 Monte Carlo simulations, take into account the source characteristics, the instrument response and the background environment which is dominated by atmospheric neutrons. This paper describes the optimisation of the polarimeter and details the associated increase in performance. The resulting design, PoGO+, is expected to improve the Minimum Detectable Polarisation (MDP) for the Crab from 19.8% to 11.1% for a 5 day flight. Assuming the same Crab polarisation fraction as measured during the 2013 flight, this improvement in MDP will allow a 5σ constrained result. It will also allow the study of the nebula emission only (Crab off-pulse) and Cygnus X-1 if in the hard state.     

Aeolian bedforms, yardangs, and indurated surfaces in the Tharsis Montes as seen by the HiRISE Camera: Evidence for dust aggregates

HiRISE images of Mars with ground sampling down to 25 cm/pixel show that the dust-rich mantle covering the surfaces of the Tharsis Montes is organized into ridges whose form and distribution are consistent with formation by aeolian saltation. Other dusty areas near the volcanoes and elsewhere on the planet exhibit a similar morphology. The material composing these “reticulate” bedforms is constrained by their remote sensing properties and the threshold curve combined with the saltation/suspension boundary, both of which vary as a function of elevation (atmospheric pressure), particle size, and particle composition. Considering all of these factors, dust aggregates are the most likely material composing these bedforms. We propose that airfall dust on and near the volcanoes aggregates in situ over time, maybe due to electrostatic charging followed by cementation by salts. The aggregates eventually reach a particle size at which saltation is possible. Aggregates on the flanks are transported downslope by katabatic winds and form linear and “accordion” morphologies. Materials within the calderas and other depressions remain trapped and are subjected to multidirectional winds, forming an interlinked “honeycomb” texture. In many places on and near the volcanoes, light-toned, low thermal inertia yardangs and indurated surfaces are present. These may represent “duststone” formed when aggregates reach a particle size below the threshold curve, such that they become stabilized and subsequently undergo cementation.     

Photometric and spectral analysis of the distribution of crystalline and amorphous ices on Enceladus as seen by Cassini

Photometric and spectral analysis of data from the Cassini Visual and Infrared Mapping Spectrometer (VIMS) has yielded significant results regarding the properties and composition of the surface of Saturn's satellite Enceladus. We have obtained spectral cubes of this satellite, containing both spatial and spectral information, with a wavelength distribution in the infrared far more extensive than from any previous observations and at much higher spatial resolution. Using a composite mosaic of the satellite, we map the distribution of crystalline and amorphous ices on the surface of Enceladus according to a “crystallinity factor” and also the depth of the temperature- and structure-dependent 1.65 micron water-ice band. These maps show the surface of Enceladus to be mostly crystalline, with a higher degree of crystallinity at the “tiger-stripe” cracks and a larger amorphous signature between these stripes. These results suggest recent geological activity at the “tiger stripe” cracks and an intriguing atmospheric environment over the south pole where amorphous ice is produced either through intense radiative bombardment, flash-freezing of cryovolcanic liquid, or rapid condensation of water vapor particles on icy microspherules or on the surface of Enceladus.     

Scale-dependent infrared radiative damping rates on Mars and their role in the deposition of gravity-wave momentum flux

Using a Curtis-matrix model of 15 μm CO₂ radiative cooling rates for the martian atmosphere, we have computed vertical scale-dependent IR radiative damping rates from 0 to 200 km altitude over a broad band of vertical wavenumbers ∣m∣ = 2π(1-500 km)⁻¹ for representative meteorological conditions at 40°N and average levels of solar activity and dust loading. In the middle atmosphere, infrared (IR) radiative damping rates increase with decreasing vertical scale and peak in excess of 30 days⁻¹ at ∼50-80 km altitude, before gradually transitioning to scale-independent rates above ∼100 km due to breakdown of local thermodynamic equilibrium. We incorporate these computed IR radiative damping rates into a linear anelastic gravity-wave model to assess the impact of IR radiative damping, relative to wave breaking and molecular viscosity, in the dissipation of gravity-wave momentum flux. The model results indicate that IR radiative damping is the dominant process in dissipating gravity-wave momentum fluxes at ∼0-50 km altitude, and is the dominant process at all altitudes for gravity waves with vertical wavelengths ?10-15 km. Wave breaking becomes dominant at higher altitudes only for “fast” waves of short horizontal and long vertical wavelengths. Molecular viscosity plays a negligible role in overall momentum flux deposition. Our results provide compelling evidence that IR radiative damping is a major, and often dominant physical process controlling the dissipation of gravity-wave momentum fluxes on Mars, and therefore should be incorporated into future parameterizations of gravity-wave drag within Mars GCMs. Lookup tables for doing so, based on the current computations, are provided.     

Adaptive optics observations of small moons of saturn

We report near-infrared observations of Prometheus and Janus taken on 9 and 13 November 2000 (UT) with the Palomar Adaptive Optics System on the 5-m Hale telescope at Palomar Observatory. Dione, Rhea, and Tethys were used as guide “stars” for the adaptive optics system, and, though they were outside the isoplanatic patch of the region of interest, they allowed significant correction of the atmospheric turbulence.Prometheus, which is usually impossible to observe from the ground due to scattered light from the A ring, was imaged at superior conjunction with Saturn. At the time of the observations, the rings of Saturn were blocked by the southern limb of the planet while the moon passed just 0.35″ below the planet’s south pole. A K filter, in a methane absorption band, was used to suppress light from the disk of the planet, and template subtraction removed much of the scattered light from the A ring. Prometheus was found to be 21.9 ± 0.1° of mean longitude behind the position predicted by Voyager-era ephemerides, consistent with the orbital lag discovered during the 1995 ring-plane crossing.     

Monte Carlo simulations of supernova light curves and the hypernova SN 1998bw

We propose a method for computing the bolometric light curves of type Ia and Ib/c supernovae based on Monte Carlo simulations of unsteady-state radiative transfer. The method is used to analyze the bolometric light curve of the unusual type Ib/c supernova SN 1998bw associated with GRB 980 425. We show that the previously noted inconsistency in the behavior of simulated light curves at early and late stages, which is attributable to asymmetry effects, can be overcome in a spherically symmetric model. Agreement with observations requires complete ⁵⁶Ni mixing and a higher matter density in the central part of the envelope in the velocity range v<5000 km s^?1 compared to standard models.     

Zero secular torque on asteroids from impinging solar photons in the YORP effect: A simple proof

YORP torques, where “YORP” stands for “Yarokovsky-O’Keefe-Radzievskii-Paddack,” arise mainly from sunlight reflected off a Solar System object and the infrared radiation emitted by it. We show here, through the most elementary demonstration that we can devise, that secular torques from impinging solar photons are generally negligible and thus cause little secular evolution of an asteroid’s obliquity or spin rate.     

The impact and rotational light curves of Comet 9P/Tempel 1

UVES and HIRES high-resolution spectra of Comet 9P/Tempel 1 are used to investigate the impact and rotational light curves of various species with a view toward building a simple model of the distribution and activity of the sources. The emission by OH, NH, CN, C₃, CH, C₂, NH₂, and OI, are analyzed, as well as the light scattered by the dust. It is found that a simple model reproduces fairly well the impact light curves of all species combining the production of the observed molecules and the expansion of the material throughout the slit. The impact light curves are consistent with velocities of 400-600 m/s. Their modeling requires a three-step dissociation sequence “Grand-Parent → Parent → Daughter” to produce the observed molecules. The rotational light curve for each species is explained in terms of a single model with three sources. The dust component can however not easily be explained that way.