Effect of turbulence on collisions of dust particles with planetesimals in protoplanetary disks

Planetesimals in gaseous protoplanetary disks may grow by collecting dust particles. Hydrodynamical studies show that small particles generally avoid collisions with the planetesimals because they are entrained by the flow around them. This occurs when $St$, the Stokes number, defined as the ratio of the dust stopping time to the planetesimal crossing time, becomes much smaller than unity. However, these studies have been limited to the laminar case, whereas these disks are believed to be turbulent. We want to estimate the influence of gas turbulence on the dust-planetesimal collision rate and on the impact speeds. We used three-dimensional direct numerical simulations of a fixed sphere (planetesimal) facing a laminar and turbulent flow seeded with small inertial particles (dust) subject to a Stokes drag. A no-slip boundary condition on the planetesimal surface is modeled via a penalty method. We find that turbulence can significantly increase the collision rate of dust particles with planetesimals. For a high turbulence case (when the amplitude of turbulent fluctuations is similar to the headwind velocity), we find that the collision probability remains equal to the geometrical rate or even higher for $St\geq 0.1$, i.e., for dust sizes an order of magnitude smaller than in the laminar case. We derive expressions to calculate impact probabilities as a function of dust and planetesimal size and turbulent intensity

Ion structure in warm dense matter: benchmarking solutions of hypernetted-chain equations by first-principle simulations

We investigate the microscopic structure of strongly coupled ions in warm dense matter using ab initio simulations and hypernetted chain (HNC) equations. We demonstrate that an approximate treatment of quantum effects by weak pseudopotentials fails to describe the highly degenerate electrons in warm dense matter correctly. However, one-component HNC calculations for the ions agree well with first-principles simulations if a linearly screened Coulomb potential is used. These HNC results can be further improved by adding a short-range repulsion that accounts for bound electrons. Examples are given for recently studied light elements, lithium and beryllium, and for aluminum where the extra short-range repulsion is essential

Twisting algebras using non-commutative torsors

Non-commutative torsors (equivalently, two-cocycles) for a Hopf algebra can be used to twist comodule algebras. After surveying and extending the literature on the subject, we prove a theorem that affords a presentation by generators and relations for the algebras obtained by such twisting. We give a number of examples, including new constructions of the quantum affine spaces and the quantum tori.Comment: 27 pages. Masuoka is a new coauthor. Introduction was revised. Sections 1 and 2 were thoroughly restructured. The presentation theorem in Section 3 is now put in a more general framework and has a more general formulation. Section 4 was shortened. All examples (quantum affine spaces and tori, twisting of SL(2), twisting of the enveloping algebra of sl(2)) are left unchange

Modeling Pressure-Ionization of Hydrogen in the Context of Astrophysics

The recent development of techniques for laser-driven shock compression of hydrogen has opened the door to the experimental determination of its behavior under conditions characteristic of stellar and planetary interiors. The new data probe the equation of state (EOS) of dense hydrogen in the complex regime of pressure ionization. The structure and evolution of dense astrophysical bodies depend on whether the pressure ionization of hydrogen occurs continuously or through a ``plasma phase transition'' (PPT) between a molecular state and a plasma state. For the first time, the new experiments constrain predictions for the PPT. We show here that the EOS model developed by Saumon and Chabrier can successfully account for the data, and we propose an experiment that should provide a definitive test of the predicted PPT of hydrogen. The usefulness of the chemical picture for computing astrophysical EOS and in modeling pressure ionization is discussed.Comment: 16 pages + 4 figures, to appear in High Pressure Researc

Coalescence in the 1D Cahn-Hilliard model

We present an approximate analytical solution of the Cahn-Hilliard equation describing the coalescence during a first order phase transition. We have identified all the intermediate profiles, stationary solutions of the noiseless Cahn-Hilliard equation. Using properties of the soliton lattices, periodic solutions of the Ginzburg-Landau equation, we have construct a family of ansatz describing continuously the processus of destabilization and period doubling predicted in Langer's self similar scenario

A protosolar nebula origin for the ices agglomerated by Comet 67P/Churyumov-Gerasimenko

The nature of the icy material accreted by comets during their formation in the outer regions of the protosolar nebula is a major open question in planetary science. Some scenarios of comet formation predict that these bodies agglomerated from crystalline ices condensed in the protosolar nebula. Concurrently, alternative scenarios suggest that comets accreted amorphous ice originating from the interstellar cloud or from the very distant regions of the protosolar nebula. On the basis of existing laboratory and modeling data, we find that the N$_2$/CO and Ar/CO ratios measured in the coma of the Jupiter family comet 67P/Churyumov-Gerasimenko by the ROSINA instrument aboard the European Space Agency's Rosetta spacecraft match those predicted for gases trapped in clathrates. If these measurements are representative of the bulk N$_2$/CO and Ar/CO ratios in 67P/Churyumov-Gerasimenko, it implies that the ices accreted by the comet formed in the nebula and do not originate from the interstellar medium, supporting the idea that the building blocks of outer solar system bodies have been formed from clathrates and possibly from pure crystalline ices. Moreover, because 67P/Churyumov-Gerasimenko is impoverished in Ar and N$_2$, the volatile enrichments observed in Jupiter's atmosphere cannot be explained solely via the accretion of building blocks with similar compositions and require an additional delivery source. A potential source may be the accretion of gas from the nebula that has been progressively enriched in heavy elements due to photoevaporation.Comment: The Astrophysical Journal Letters, in pres

Elemental abundances and minimum mass of heavy elements in the envelope of HD 189733b

Oxygen (O) and carbon (C) have been inferred recently to be subsolar in abundance from spectra of the atmosphere of the transiting hot Jupiter HD 189733b. Yet, the mass and radius of the planet coupled with structure models indicate a strongly supersolar abundance of heavy elements in the interior of this object. Here we explore the discrepancy between the large amount of heavy elements suspected in the planet's interior and the paucity of volatiles measured in its atmosphere. We describe the formation sequence of the icy planetesimals formed beyond the snow line of the protoplanetary disk and calculate the composition of ices ultimately accreted in the envelope of HD 189733b on its migration pathway. This allows us to reproduce the observed volatile abundances by adjusting the mass of ices vaporized in the envelope. The predicted elemental mixing ratios should be 0.15--0.3 times solar in the envelope of HD 189733b if they are fitted to the recent O and C determinations. However, our fit to the minimum mass of heavy elements predicted by internal structure models gives elemental abundances that are 1.2--2.4 times oversolar in the envelope of HD189733b. We propose that the most likely cause of this discrepancy is irradiation from the central star leading to development of a radiative zone in the planet's outer envelope which would induce gravitational settling of elements. Hence, all strongly irradiated extrasolar planets should present subsolar abundances of volatiles. We finally predict that the abundances of nitrogen (N), sulfur (S) and phosphorus (P) are of $\sim$ $2.8 \times 10^{-5}$, $5.3 \times 10^{-6}$ and $1.8 \times 10^{-7}$ relative to H$_2$, respectively in the atmosphere of HD 189733b.Comment: Accepted for publication in Astronomy & Astrophysic

A simple method for measuring plasma power in rf-GDOES instruments

A method for determining plasma power in rf-GDOES is presented. It is based on an effective resistance located in the inductive coil of the impedance matching. The amount of electrical power consumed in the matching system depends on the capacitive current flowing through the matching system, which depends on the applied voltage, the stray capacity, and the frequency. This correction method is experimentally evaluated and compared with the integral plasma power calculatio

High frequency haplotypes are expected events, not historical figures

Cultural transmission of reproductive success states that successful men have more children and pass this raised fecundity to their offspring. Balaresque and colleagues found high frequency haplotypes in a Central Asian Y chromosome dataset, which they attribute to cultural transmission of reproductive success by prominent historical men, including Genghis Khan. Using coalescent simulation, we show that these high frequency haplotypes are consistent with a neutral model, where they commonly appear simply by chance. Hence, explanations invoking cultural transmission of reproductive success are statistically unnecessary