Kinetic Monte Carlo method for simulating astrochemical kinetics: Test calculations of molecular hydrogen formation on interstellar dust particles

Molecular hydrogen formation on the surfaces of interstellar dust grains has been studied numerically. Different stochastic methods (analogue Monte Carlo methods, direct solutions of the chemical master equation, continuous-time random-walk simulations, etc.) have been recently applied to the astrochemistry of hydrogen formation; however, the computational efficiency of these approaches is usually low because they are extremely time consuming (Herbst and Shematovich, 2003). A kinetic version of the Monte Carlo method based on splitting by physical processes is presented in the study. Each of the basic physical processes—adsorption of atomic hydrogen on the surface, thermal diffusion, formation of molecular hydrogen, and hydrogen desorption from the surface-are considered as independent random Markovian-type processes and are simulated using stochastic algorithms. The suggested numerical model is computationally efficient allowing calculations to be held on model surfaces with up to 10⁶ active centers, which corresponds to interstellar grains of up to 0.4 μm in diameter. Test calculations of the efficiency of molecular hydrogen formation were held for the two models of the surfaces of interplanetary grains, namely, for the homogeneous surfaces of olivine and amorphous carbon. It was confirmed that the effective formation of molecular hydrogen in diffuse molecular clouds is possible only within a narrow range of dust-particle temperatures. A comparison with the numerical results of other authors is presented.     

The annual period of the Earth’s seismic activity for 1964–2007

The catalog of the United States National Earthquake Information Center (NEIC, 2007) was used for a Fourier analysis of planetary seismic activity from 1964 to 2007 (401219 earthquakes with M ≥ 3 and hypocenter depths H ≥ 1 km) for the Northern Hemisphere (248291 events) and for the Southern Hemisphere (152928 events). The annual periodicity of weak earthquakes (M < 5.0) was verified with a high degree of reliability. All regularities (depending on the geographic latitude, hypocenter depths, and north-south asymmetry) revealed earlier (in 1964–1990) for this period are shown to exist for the period of 1964–2007.     

Melting and Sublimation of Planetary Ices Under Low Pressure Conditions: Laboratory Experiments with a Melting Probe Prototype

One possibility to explore the subsurface layers of icy bodies is to use a probe with a “hot tip", which is able to penetrate ice layers by melting. Such probes have been built and used in the past for the exploration of terrestrial polar ice sheets and may also become useful tools to explore other icy layers in the Solar System. Examples for such layers are the polar areas of Mars or the icy crust of Jupiter’s moon Europa. However, while on Earth a heated probe launched into an ice sheet always causes melting with subsequent refreezing, the behaviour of such a probe in a low pressure environment is quite different. We report on the results of some experiments with a simple “melting probe" prototype with two different kinds of hot tips in a vacuum environment. For one of the tips the probe moved into two types of ice samples: (i) compact water ice and (ii) porous water ice with a snow (firn) like texture. It was also found that the penetration behaviour was basically different for the two sample types even when the same kind of tip was used. While in the porous sample the ice was only subliming, the phase changes occurring during the interaction of the tip with the compact ice are much more complex. Here alternating phases of melting and sublimation occur. The absence of the liquid phase has severe consequences on the performance of a “melting probe" under vacuum conditions: In this environment we find a high thermal resistance between the probe surface and the underlying ice. Therefore, only a low percentage of the heat that is generated in the tip is used to melt or sublime the ice, the bulk of the power is transferred towards the rear end of the probe. This is particularly a problem in the initial phases of an ice penetration experiment, when the probe has not yet penetrated the ice over its whole length. In the compact ice sample, phases could be observed, where a high enough gas pressure had built up locally underneath the probe, so that melting becomes possible. Only during these melting periods the thermal contact between the probe and the ice is good and in consequence the melting probe works effectively.     

Evolution of the photospheric magnetic field and coronal null points before solar flares

Based on a topological model for the magnetic field of a solar active region (AR), we suggest a criterion for the existence of magnetic null points on the separators in the corona. With the problem of predicting solar flares in mind, we have revealed a model parameter whose decrease means that the AR evolves toward a major eruptive flare. We analyze the magnetic field evolution for AR 9077 within two days before the Bastille Day flare on July 14, 2000. The coronal conditions are shown to have become more favorable for magnetic reconnection, which led to a 3B/X5.7 eruptive flare.     

Can the known millisecond pulsars help in the detection of intermediate-mass black holes at the centers of globular clusters?

We consider the possibility of detecting intermediate-mass (10³–10⁴ M _⊙) black holes, whose existence at the centers of globular clusters is expected from optical and infrared observations, using precise pulse arrival timing for the millisecond pulsars in globular clusters known to date. For some of these pulsars closest to the cluster centers, we have calculated the expected delay times of pulses as they pass in the gravitational field of the central black hole. The detection of such a time delay by currently available instruments for the known pulsars is shown to be impossible at a black hole mass of 10³ M _⊙ and very problematic at a black hole mass of 10⁴ M _⊙. In addition, the signal delay will have a negligible effect on the pulsar periods and their first derivatives compared to the current accuracy of their measurements.     

The DynaMICCS perspective

The DynaMICCS mission is designed to probe and understand the dynamics of crucial regions of the Sun that determine solar variability, including the previously unexplored inner core, the radiative/convective zone interface layers, the photosphere/chromosphere layers and the low corona. The mission delivers data and knowledge that no other known mission provides for understanding space weather and space climate and for advancing stellar physics (internal dynamics) and fundamental physics (neutrino properties, atomic physics, gravitational moments...). The science objectives are achieved using Doppler and magnetic measurements of the solar surface, helioseismic and coronographic measurements, solar irradiance at different wavelengths and in-situ measurements of plasma/energetic particles/magnetic fields. The DynaMICCS payload uses an original concept studied by Thalès Alenia Space in the framework of the CNES call for formation flying missions: an external occultation of the solar light is obtained by putting an occulter spacecraft 150 m (or more) in front of a second spacecraft. The occulter spacecraft, a LEO platform of the mini sat class, e.g. PROTEUS, type carries the helioseismic and irradiance instruments and the formation flying technologies. The latter spacecraft of the same type carries a visible and infrared coronagraph for a unique observation of the solar corona and instrumentation for the study of the solar wind and imagers. This mission must guarantee long (one 11-year solar cycle) and continuous observations (duty cycle > 94%) of signals that can be very weak (the gravity mode detection supposes the measurement of velocity smaller than 1 mm/s). This assumes no interruption in observation and very stable thermal conditions. The preferred orbit therefore is the L1 orbit, which fits these requirements very well and is also an attractive environment for the spacecraft due to its low radiation and low perturbation (solar pressure) environment. This mission is secured by instrumental R and D activities during the present and coming years. Some prototypes of different instruments are already built (GOLFNG, SDM) and the performances will be checked before launch on the ground or in space through planned missions of CNES and PROBA ESA missions (PICARD, LYRA, maybe ASPIICS).     

Accretion disc in the massive V448 Cygni system

The UBV light curves of the early-type eclipsing binary V448 Cygni, obtained at the Abastumani Astrophysical Observatory from 1964 to 1967, are re-analysed here. The analysis was made assuming the presence of an accretion disc in the system, as inferred from the light-curve shape and spectroscopic characteristics of the system. The Roche model of a binary was used, containing a geometrically and optically thick accretion disc around the hotter and more massive star. By solving the inverse problem, the orbital elements and the physical parameters of the system components and of the accretion disc were estimated. This result is important for understanding the star formation and evolution processes in the systems with massive components.     

Merging time-scales of stellar subclumps in young star-forming regions

Recent observations and hydrodynamical simulations of star formation inside a giant molecular cloud have revealed that, within a star-forming region, stars do not form evenly distributed throughout this region, but rather in small subclumps. It is generally believed that these subclumps merge and form a young star cluster. The time-scale of this merging process is crucial for the evolution and the possible survival of the final star cluster. The key issue is whether this merging process happens faster than the time needed to remove the residual gas of the cloud. A merging time-scale shorter than the gas-removal time would enhance the survival chances of the resulting star cluster. In this paper, we show by means of numerical simulations that the time-scale of the merging is indeed very fast. Depending on the details of the initial subclump distribution, the merging may occur before the gas is expelled from the newly formed cluster via either supernovae or the winds from massive stars. Our simulations further show that the resulting merger objects have a higher effective star formation efficiency than the overall star-forming region and confirm the results that mass-segregated subclumps form mass-segregated merger objects.     

A dark matter disc in three cosmological simulations of Milky Way mass galaxies

Making robust predictions for the phase-space distribution of dark matter at the solar neighbourhood is vital for dark matter direct-detection experiments. To date, almost all such predictions have been based on simulations that model the dark matter alone. Here, we use three cosmological hydrodynamic simulations of bright, disc-dominated galaxies to include the effects of baryonic matter self-consistently for the first time. We find that the addition of baryonic physics drastically alters the dark matter profile in the vicinity of the solar neighbourhood. A stellar/gas disc, already in place at high redshift, causes merging satellites to be dragged preferentially towards the disc plane where they are torn apart by tides. This results in an accreted dark matter disc that contributes ∼0.25–1.5 times the non-rotating halo density at the solar position. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts Weakly Interacting Massive Particle (WIMP) capture in the Earth and Sun, boosts the annual modulation signal and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined.     

Multi-zone fusion crust formation and classification of the 2004 Auckland meteorite (L6, S5, and W0)

On June 12, 2004, a meteorite passed through Earth's atmosphere and landed under the television in the living room of a house in Auckland, New Zealand. Textural characteristics, the chemistry of olivine (Fa_23–24) and orthopyroxene (Fs_20.7), and the bulk rock triple oxygen isotopes (δ¹⁷O + 3.1; δ¹⁸O + 4.2‰) from the interior of the completely unweathered (W0) 1.3 kg meteorite, hereafter referred to as Auckland, suggest it to be a strongly metamorphosed fragment from the interior of a low iron ordinary chondrite (L6) parent asteroid. The occurrence of maskelynite but shock fracturing of olivine and pyroxene indicates Auckland experienced extreme shock metamorphism (S5), likely during Ordovician fragmentation of the asteroid parent. The fusion crust consists of three zones: (1) an innermost zone containing narrow Fe-Ni-S-bearing veins that migrated along pre-existing shock fractures in olivine and pyroxene; (2) a middle zone in which the meteorite partially melted to form a silicate glass and immiscible blebs of metal and troilite, and is accompanied by unmelted silicate minerals; and (3) an approximately 0.1 mm wide vesicular-rich outermost layer that largely melted, volatilizing sulfides, before quenching to form glass and olivine. Oxygen isotope values of the bulk rock and/or maskelynite of melted rim and modified substrate are 2–3‰ greater than the meteorite interior and indicate that up to 19% of terrestrial atmospheric O₂ was incorporated into the fusion crust during the formation. The fusion crust migrated inwards as ablation occurred, enabling melting, migration, and re-precipitation ± loss of sulfide and metal components, with the prominent glassy rim therefore forming from an already chemically modified zone.