N-body simulations of oligarchic growth of Mars: Implications for Hf-W chronology

Dauphas and Pourmand (2011) [Nature 473, 489--492] estimated the accretion timescale of Mars to be 1.8 $^{+0.9}_{-1.0}$ Myr from the W isotopes of martian meteorites. This timescale was derived assuming perfect metal-silicate equilibration between the impactor and the target's mantle. However, in the case of a small impactor most likely only a fraction of the target's mantle is involved in the equilibration, while only a small part of the impactor's core equilibrates in the case of a giant impact. We examined the effects of imperfect equilibration using results of high-resolution $N$-body simulations for the oligarchic growth stage. These effects were found to be small as long as a planetary embryo has a deep liquid magma ocean during its accretion. The effect due to partial involvement of the target's mantle in equilibration is small due to the low metal-silicate partition coefficient for W suggested from the low Hf/W ratio of the martian mantle. The effect due to partial involvement of the impactor's core is also small because a large fraction of the embryo mass is delivered from small planetesimals, which are likely to fully equilibrate in the deep magma ocean on the embryo. The accretion timescale of Mars estimated by the Hf-W chronology is shorter than that expected for the minimum mass solar nebula model as long as more than 10% of each impactor's core re-equilibrates with the martian mantle and the final stages of accretion are prolonged. This probably indicates that accretion of Mars proceeded rapidly due to solid and gas surface densities significantly larger than those for the minimum mass solar nebula or due to accretion of small fragments or pebbles.Comment: 37 pages, 10 figures, accepted for publication in EPS

Actinide collisions for QED and superheavy elements with the time-dependent Hartree-Fock theory and the Balian-V\'en\'eroni variational principle

Collisions of actinide nuclei form, during very short times of few zs ($10^{-21}$ s), the heaviest ensembles of interacting nucleons available on Earth. Such collisions are used to produce super-strong electric fields by the huge number of interacting protons to test spontaneous positron-electron pair emission (vacuum decay) predicted by the quantum electrodynamics (QED) theory. Multi-nucleon transfer in actinide collisions could also be used as an alternative way to fusion in order to produce neutron-rich heavy and superheavy elements thanks to inverse quasifission mechanisms. Actinide collisions are studied in a dynamical quantum microscopic approach. The three-dimensional time-dependent Hartree-Fock (TDHF) code {\textsc{tdhf3d}} is used with a full Skyrme energy density functional to investigate the time evolution of expectation values of one-body operators, such as fragment position and particle number. This code is also used to compute the dispersion of the particle numbers (e.g., widths of fragment mass and charge distributions) from TDHF transfer probabilities, on the one hand, and using the Balian-Veneroni variational principle, on the other hand. A first application to test QED is discussed. Collision times in $^{238}$U+$^{238}$U are computed to determine the optimum energy for the observation of the vacuum decay. It is shown that the initial orientation strongly affects the collision times and reaction mechanism. The highest collision times predicted by TDHF in this reaction are of the order of $\sim4$ zs at a center of mass energy of 1200 MeV. According to modern calculations based on the Dirac equation, the collision times at $E_{cm}>1$ GeV are sufficient to allow spontaneous electron-positron pair emission from QED vacuum decay, in case of bare uranium ion collision. A second application of actinide collisions to produce neutron-rich transfermiums is discussed. A new inverse quasifission mechanism associated to a specific orientation of the nuclei is proposed to produce transfermium nuclei ($Z>100$) in the collision of prolate deformed actinides such as $^{232}$Th+$^{250}$Cf. The collision of the tip of one nucleus with the side of the other results in a nucleon flux toward the latter. The probability distributions for transfermium production in such a collision are computed. The produced nuclei are more neutron-rich than those formed in fusion reactions, thus, leading to more stable isotopes closer to the predicted superheavy island of stability. In addition to mass and charge dispersion, the Balian-Veneroni variational principle is used to compute correlations between $Z$ and $N$ distributions, which are zero in standard TDHF calculations.Comment: Proceeding of the FUSION11 conferenc

The effects of short-lived radionuclides and porosity on the early thermo-mechanical evolution of planetesimals

The thermal history and internal structure of chondritic planetesimals, assembled before the giant impact phase of chaotic growth, potentially yield important implications for the final composition and evolution of terrestrial planets. These parameters critically depend on the internal balance of heating versus cooling, which is mostly determined by the presence of short-lived radionuclides (SLRs), such as aluminum-26 and iron-60, as well as the heat conductivity of the material. The heating by SLRs depends on their initial abundances, the formation time of the planetesimal and its size. It has been argued that the cooling history is determined by the porosity of the granular material, which undergoes dramatic changes via compaction processes and tends to decrease with time. In this study we assess the influence of these parameters on the thermo-mechanical evolution of young planetesimals with both 2D and 3D simulations. Using the code family I2ELVIS/I3ELVIS we have run numerous 2D and 3D numerical finite-difference fluid dynamic models with varying planetesimal radius, formation time and initial porosity. Our results indicate that powdery materials lowered the threshold for melting and convection in planetesimals, depending on the amount of SLRs present. A subset of planetesimals retained a powdery surface layer which lowered the thermal conductivity and hindered cooling. The effect of initial porosity was small, however, compared to those of planetesimal size and formation time, which dominated the thermo-mechanical evolution and were the primary factors for the onset of melting and differentiation. We comment on the implications of this work concerning the structure and evolution of these planetesimals, as well as their behavior as possible building blocks of terrestrial planets.Comment: 19 pages, 11 figures, 5 tables; accepted for publication in Icarus; for associated video files, see http://timlichtenberg.net/2015_porosity.html or http://dx.doi.org/10.1016/j.icarus.2016.03.00

Microscopic description of heavy ion collisions around the barrier

7 pages, 9 figuresInternational audienceA microscopic mean-field description of heavy ion collisions is performed in the framework of the time dependent Hartree-Fock theory using a Skyrme energy density functional. A good agreement with experiments is obtained on the position of the fusion barriers for various total masses and mass asymmetries. The excitation function of the 16O+208Pb is overestimated by about 16% above the barrier. Transfer below the barrier is studied in 16O+208Pb central collisions. Heavier systems are considered to study their fusion hindrance. We also compute collision times of the 238U+238U system. The latter has been used to produce super strong electric fields and to test non perturbative quantum electrodynamics theory. Indeed, if the life time of such giant system is of the order of few 10^-21 s, its electric field should induce spontaneous electron-positron pair emissions from vacuum. In our calculations, highest collision times are reached in the 238U+238U reaction for center of mass energies between 1000 and 1300 MeV

Coupling SPH and thermochemical models of planets: Methodology and example of a Mars-sized body

Giant impacts have been suggested to explain various characteristics of terrestrial planets and their moons. However, so far in most models only the immediate effects of the collisions have been considered, while the long-term interior evolution of the impacted planets was not studied. Here we present a new approach, combining 3-D shock physics collision calculations with 3-D thermochemical interior evolution models. We apply the combined methods to a demonstration example of a giant impact on a Mars-sized body, using typical collisional parameters from previous studies. While the material parameters (equation of state, rheology model) used in the impact simulations can have some effect on the long-term evolution, we find that the impact angle is the most crucial parameter for the resulting spatial distribution of the newly formed crust. The results indicate that a dichotomous crustal pattern can form after a head-on collision, while this is not the case when considering a more likely grazing collision. Our results underline that end-to-end 3-D calculations of the entire process are required to study in the future the effects of large-scale impacts on the evolution of planetary interiors.Comment: 29 pages, 10 figures, accepted for publication in Icaru

Late metal-silicate separation on the IAB parent asteroid: Constraints from combined W and Pt isotopes and thermal modelling

The short-lived $^{182}$Hf-$^{182}$W decay system is a powerful chronometer for constraining the timing of metal-silicate separation and core formation in planetesimals and planets. Neutron capture effects on W isotopes, however, significantly hamper the application of this tool. In order to correct for neutron capture effects, Pt isotopes have emerged as a reliable in-situ neutron dosimeter. This study applies this method to IAB iron meteorites, in order to constrain the timing of metal segregation on the IAB parent body. The $\epsilon^{182}$W values obtained for the IAB iron meteorites range from -3.61 $\pm$ 0.10 to -2.73 $\pm$ 0.09. Correlating $\epsilon^{\mathrm{i}}$Pt with $^{182}$W data yields a pre-neutron capture $^{182}$W of -2.90 $\pm$ 0.06. This corresponds to a metal-silicate separation age of 6.0 $\pm$ 0.8 Ma after CAI for the IAB parent body, and is interpreted to represent a body-wide melting event. Later, between 10 and 14 Ma after CAI, an impact led to a catastrophic break-up and subsequent reassembly of the parent body. Thermal models of the interior evolution that are consistent with these estimates suggest that the IAB parent body underwent metal-silicate separation as a result of internal heating by short-lived radionuclides and accreted at around 1.4 $\pm$ 0.1 Ma after CAIs with a radius of greater than 60 km.Comment: 11 pages, 8 figures, 2 tables; open access article under the CC BY-NC-ND license (see http://creativecommons.org/licenses/by-nc-nd/4.0/

Biomechanics of predator–prey arms race in lion, zebra, cheetah and impala

The fastest and most manoeuvrable terrestrial animals are found in savannah habitats, where predators chase and capture running prey. Hunt outcome and success rate are critical to survival, so both predator and prey should evolve to be faster and/or more manoeuvrable. Here we compare locomotor characteristics in two pursuit predator–prey pairs, lion–zebra and cheetah–impala, in their natural savannah habitat in Botswana. We show that although cheetahs and impalas were universally more athletic than lions and zebras in terms of speed, acceleration and turning, within each predator–prey pair, the predators had 20% higher muscle fibre power than prey, 37% greater acceleration and 72% greater deceleration capacity than their prey. We simulated hunt dynamics with these data and showed that hunts at lower speeds enable prey to use their maximum manoeuvring capacity and favour prey survival, and that the predator needs to be more athletic than its prey to sustain a viable success rate

MIMAC : A micro-tpc matrix for directional detection of dark matter

Directional detection of non-baryonic Dark Matter is a promising search strategy for discriminating WIMP events from background. However, this strategy requires both a precise measurement of the energy down to a few keV and 3D reconstruction of tracks down to a few mm. To achieve this goal, the MIMAC project has been developed. It is based on a gaseous micro-TPC matrix, filled with CF4 and CHF3. The first results on low energy nuclear recoils (H, F) obtained with a low mono-energetic neutron field are presented. The discovery potential of this search strategy is discussed and illustrated by a realistic case accessible to MIMAC.Comment: 6 pages, Proc. of the fifth international symposium on large TPCs for low energy rare event detection, Paris, France, Dec. 2010. To appear in Journal of Physic

Henri Temianka Correspondence; (golabek)

This collection contains material pertaining to the life, career, and activities of Henri Temianka, violin virtuoso, conductor, music teacher, and author. Materials include correspondence, concert programs and flyers, music scores, photographs, and books.https://digitalcommons.chapman.edu/temianka_correspondence/3516/thumbnail.jp