• ### HI 21cm Cosmology and the Bi-spectrum: Closure Diagnostics in Massively Redundant Interferometric Arrays(1805.00953)

May 2, 2018 astro-ph.IM
New massively redundant low frequency arrays allow for a novel investigation of closure relations in interferometry. We employ commissioning data from the Hydrogen Epoch of Reionization Array to investigate closure quantities in this densely packed grid array of 14m antennas operating at 100 MHz to 200 MHz. We investigate techniques that utilize closure phase spectra for redundant triads to estimate departures from redundancy for redundant baseline visibilities. We find a median absolute deviation from redundancy in closure phase across the observed frequency range of about 4.5deg. This value translates into a non-redundancy per visibility phase of about 2.6deg, using prototype electronics. The median absolute deviations from redundancy decrease with longer baselines. We show that closure phase spectra can be used to identify ill-behaved antennas in the array, independent of calibration. We investigate the temporal behavior of closure spectra. The Allan variance increases after a one minute stride time, due to passage of the sky through the primary beam of the transit telescope. However, the closure spectra repeat to well within the noise per measurement at corresponding local sidereal times (LST) from day to day. In future papers in this series we will develop the technique of using closure phase spectra in the search for the HI 21cm signal from cosmic reionization.
• Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.
• ### A New Data Source for Inverse Dynamics Learning(1710.02513)

Oct. 6, 2017 cs.RO
Modern robotics is gravitating toward increasingly collaborative human robot interaction. Tools such as acceleration policies can naturally support the realization of reactive, adaptive, and compliant robots. These tools require us to model the system dynamics accurately -- a difficult task. The fundamental problem remains that simulation and reality diverge--we do not know how to accurately change a robot's state. Thus, recent research on improving inverse dynamics models has been focused on making use of machine learning techniques. Traditional learning techniques train on the actual realized accelerations, instead of the policy's desired accelerations, which is an indirect data source. Here we show how an additional training signal -- measured at the desired accelerations -- can be derived from a feedback control signal. This effectively creates a second data source for learning inverse dynamics models. Furthermore, we show how both the traditional and this new data source, can be used to train task-specific models of the inverse dynamics, when used independently or combined. We analyze the use of both data sources in simulation and demonstrate its effectiveness on a real-world robotic platform. We show that our system incrementally improves the learned inverse dynamics model, and when using both data sources combined converges more consistently and faster.
• ### Real-time Perception meets Reactive Motion Generation(1703.03512)

Oct. 6, 2017 cs.RO
We address the challenging problem of robotic grasping and manipulation in the presence of uncertainty. This uncertainty is due to noisy sensing, inaccurate models and hard-to-predict environment dynamics. We quantify the importance of continuous, real-time perception and its tight integration with reactive motion generation methods in dynamic manipulation scenarios. We compare three different systems that are instantiations of the most common architectures in the field: (i) a traditional sense-plan-act approach that is still widely used, (ii) a myopic controller that only reacts to local environment dynamics and (iii) a reactive planner that integrates feedback control and motion optimization. All architectures rely on the same components for real-time perception and reactive motion generation to allow a quantitative evaluation. We extensively evaluate the systems on a real robotic platform in four scenarios that exhibit either a challenging workspace geometry or a dynamic environment. In 333 experiments, we quantify the robustness and accuracy that is due to integrating real-time feedback at different time scales in a reactive motion generation system. We also report on the lessons learned for system building.
• ### SDSS IV MaNGA - Rotation Velocity Lags in the Extraplanar Ionized Gas from MaNGA Observations of Edge-on Galaxies(1704.02582)

April 9, 2017 astro-ph.GA
We present a study of the kinematics of the extraplanar ionized gas around several dozen galaxies observed by the Mapping of Nearby Galaxies at the Apache Point Observatory (MaNGA) survey. We considered a sample of 67 edge-on galaxies out of more than 1400 extragalactic targets observed by MaNGA, in which we found 25 galaxies (or 37%) with regular lagging of the rotation curve at large distances from the galactic midplane. We model the observed $H\alpha$ emission velocity fields in the galaxies, taking projection effects and a simple model for the dust extinction into the account. We show that the vertical lag of the rotation curve is necessary in the modeling, and estimate the lag amplitude in the galaxies. We find no correlation between the lag and the star formation rate in the galaxies. At the same time, we report a correlation between the lag and the galactic stellar mass, central stellar velocity dispersion, and axial ratio of the light distribution. These correlations suggest a possible higher ratio of infalling-to-local gas in early-type disk galaxies or a connection between lags and the possible presence of hot gaseous halos, which may be more prevalent in more massive galaxies. These results again demonstrate that observations of extraplanar gas can serve as a potential probe for accretion of gas.
• SNO+ is a large liquid scintillator-based experiment located 2km underground at SNOLAB, Sudbury, Canada. It reuses the Sudbury Neutrino Observatory detector, consisting of a 12m diameter acrylic vessel which will be filled with about 780 tonnes of ultra-pure liquid scintillator. Designed as a multipurpose neutrino experiment, the primary goal of SNO+ is a search for the neutrinoless double-beta decay (0$\nu\beta\beta$) of 130Te. In Phase I, the detector will be loaded with 0.3% natural tellurium, corresponding to nearly 800 kg of 130Te, with an expected effective Majorana neutrino mass sensitivity in the region of 55-133 meV, just above the inverted mass hierarchy. Recently, the possibility of deploying up to ten times more natural tellurium has been investigated, which would enable SNO+ to achieve sensitivity deep into the parameter space for the inverted neutrino mass hierarchy in the future. Additionally, SNO+ aims to measure reactor antineutrino oscillations, low-energy solar neutrinos, and geoneutrinos, to be sensitive to supernova neutrinos, and to search for exotic physics. A first phase with the detector filled with water will begin soon, with the scintillator phase expected to start after a few months of water data taking. The 0$\nu\beta\beta$ Phase I is foreseen for 2017.
• ### Rapid Evolution of the Photosystem II Electronic Structure during Water Splitting(1506.08862)

Photosynthetic water oxidation is a fundamental process that sustains the biosphere. A Mn$_{4}$Ca cluster embedded in the photosystem II protein environment is responsible for the production of atmospheric oxygen. Here, time-resolved x-ray emission spectroscopy (XES) was used to observe the process of oxygen formation in real time. These experiments reveal that the oxygen evolution step, initiated by three sequential laser flashes, is accompanied by rapid (within 50 $\mu$s) changes to the Mn K$\beta$ XES spectrum. However, no oxidation of the Mn$_{4}$Ca core above the all Mn$^{\text{IV}}$ state was detected to precede O-O bond formation. A new mechanism featuring Mn$^{\text{IV}}$=O formation in the S$_{3}$ state is proposed to explain the spectroscopic results. This chemical formulation is consistent with the unique reactivity of the S$_{3}$ state and explains facilitation of the following S$_{3}$ to S$_{0}$ transition, resolving in part the kinetic limitations associated with O-O bond formation. In the proposed mechanism, O-O bond formation precedes transfer of the final (4$^{\text{th}}$) electron from the Mn$_{4}$Ca cluster, in agreement with experiment.
• The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.
• Super-Kamiokande (SK) can search for weakly interacting massive particles (WIMPs) by detecting neutrinos produced from WIMP annihilations occurring inside the Sun. In this analysis, we include neutrino events with interaction vertices in the detector in addition to upward-going muons produced in the surrounding rock. Compared to the previous result, which used the upward-going muons only, the signal acceptances for light (few-GeV/$c^2$ $\sim$ 200-GeV/$c^2$) WIMPs are significantly increased. We fit 3903 days of SK data to search for the contribution of neutrinos from WIMP annihilation in the Sun. We found no significant excess over expected atmospheric-neutrino background and the result is interpreted in terms of upper limits on WIMP-nucleon elastic scattering cross sections under different assumptions about the annihilation channel. We set the current best limits on the spin-dependent (SD) WIMP-proton cross section for WIMP masses below 200 GeV/$c^2$ (at 10 GeV/$c^2$, 1.49$\times 10^{-39}$ cm$^2$ for $\chi\chi\rightarrow b \bar{b}$ and 1.31$\times 10^{-40}$ cm$^2$ for $\chi\chi\rightarrow\tau^+\tau^-$ annihilation channels), also ruling out some fraction of WIMP candidates with spin-independent (SI) coupling in the few-GeV/$c^2$ mass range.
• ### Study of Redshifted HI from the Epoch of Reionization with Drift scan(1407.4620)

July 17, 2014 astro-ph.CO, astro-ph.IM
The detection of the Epoch of Reionization (EoR) in the redshifted 21-cm line is a challenging task. Here we formulate the detection of the EoR signal using the drift scan strategy. This method potentially has better instrumental stability as compared to the case where a single patch of sky is tracked. We demonstrate that the correlation time between measured visibilities could extend up to 1-2 hr for an interferometer array such as the Murchison Widefield Array (MWA), which has a wide primary beam. We estimate the EoR power based on cross-correlation of visibilities across time and show that the drift scan strategy is capable of the detection of the EoR signal with comparable/better signal-to-noise as compared to the tracking case. We also estimate the visibility correlation for a set of bright point sources and argue that the statistical inhomogeneity of bright point sources might allow their separation from the EoR signal.
• Novalike cataclysmic variables have persistently high mass transfer rates and prominent steady state accretion disks. We present an analysis of infrared observations of twelve novalikes obtained from the Two Micron All Sky Survey, the Spitzer Space Telescope, and the Wide-field Infrared Survey Explorer All Sky Survey. The presence of an infrared excess at >3-5 microns over the expectation of a theoretical steady state accretion disk is ubiquitous in our sample. The strength of the infrared excess is not correlated with orbital period, but shows a statistically significant correlation (but shallow trend) with system inclination that might be partially (but not completely) linked to the increasing view of the cooler outer accretion disk and disk rim at higher inclinations. We discuss the possible origin of the infrared excess in terms of emission from bremsstrahlung or circumbinary dust, with either mechanism facilitated by the mass outflows (e.g., disk wind/corona, accretion stream overflow, and so on) present in novalikes. Our comparison of the relative advantages and disadvantages of either mechanism for explaining the observations suggests that the situation is rather ambiguous, largely circumstantial, and in need of stricter observational constraints.
• ### Stellar Superfluids(1302.6626)

Neutron stars provide a fertile environment for exploring superfluidity under extreme conditions. It is not surprising that Cooper pairing occurs in dense matter since nucleon pairing is observed in nuclei as energy differences between even-even and odd-even nuclei. Since superfluids and superconductors in neutron stars profoundly affect neutrino emissivities and specific heats, their presence can be observed in the thermal evolution of neutron stars. An ever-growing number of cooling neutron stars, now amounting to 13 thermal sources, and several additional objects from which upper limits to temperatures can be ascertained, can now be used to discriminate among theoretical scenarios and even to dramatically restrict properties of nucleon pairing at high densities. In addition, observations of pulsars, including their spin-downs and glitch histories, additionally support the conjecture that superfluidity and superconductivity are ubiquitous within, and important to our understanding of, neutron stars.
• ### Exomoon habitability constrained by illumination and tidal heating(1209.5323)

Oct. 30, 2013 astro-ph.EP
The detection of moons orbiting extrasolar planets ("exomoons") has now become feasible. Once they are discovered in the circumstellar habitable zone, questions about their habitability will emerge. Exomoons are likely to be tidally locked to their planet and hence experience days much shorter than their orbital period around the star and have seasons, all of which works in favor of habitability. These satellites can receive more illumination per area than their host planets, as the planet reflects stellar light and emits thermal photons. On the contrary, eclipses can significantly alter local climates on exomoons by reducing stellar illumination. In addition to radiative heating, tidal heating can be very large on exomoons, possibly even large enough for sterilization. We identify combinations of physical and orbital parameters for which radiative and tidal heating are strong enough to trigger a runaway greenhouse. By analogy with the circumstellar habitable zone, these constraints define a circumplanetary "habitable edge". We apply our model to hypothetical moons around the recently discovered exoplanet Kepler-22b and the giant planet candidate KOI211.01 and describe, for the first time, the orbits of habitable exomoons. If either planet hosted a satellite at a distance greater than 10 planetary radii, then this could indicate the presence of a habitable moon.
• ### Hot Moons and Cool Stars(1301.0235)

Jan. 2, 2013 astro-ph.EP
The exquisite photometric precision of the Kepler space telescope now puts the detection of extrasolar moons at the horizon. Here, we firstly review observational and analytical techniques that have recently been proposed to find exomoons. Secondly, we discuss the prospects of characterizing potentially habitable extrasolar satellites. With moons being much more numerous than planets in the solar system and with most exoplanets found in the stellar habitable zone being gas giants, habitable moons could be as abundant as habitable planets. However, satellites orbiting planets in the habitable zones of cool stars will encounter strong tidal heating and likely appear as hot moons.
• ### Constraints on the habitability of extrasolar moons(1210.5172)

Oct. 18, 2012 astro-ph.EP
Detections of massive extrasolar moons are shown feasible with the Kepler space telescope. Kepler's findings of about 50 exoplanets in the stellar habitable zone naturally make us wonder about the habitability of their hypothetical moons. Illumination from the planet, eclipses, tidal heating, and tidal locking distinguish remote characterization of exomoons from that of exoplanets. We show how evaluation of an exomoon's habitability is possible based on the parameters accessible by current and near-future technology.
• ### Transit Timing Observations from Kepler: VII. Confirmation of 27 planets in 13 multiplanet systems via Transit Timing Variations and orbital stability(1208.3499)

Aug. 16, 2012 astro-ph.EP
We confirm 27 planets in 13 planetary systems by showing the existence of statistically significant anti-correlated transit timing variations (TTVs), which demonstrates that the planet candidates are in the same system, and long-term dynamical stability, which places limits on the masses of the candidates---showing that they are planetary. %This overall method of planet confirmation was first applied to \kepler systems 23 through 32. All of these newly confirmed planetary systems have orbital periods that place them near first-order mean motion resonances (MMRs), including 6 systems near the 2:1 MMR, 5 near 3:2, and one each near 4:3, 5:4, and 6:5. In addition, several unconfirmed planet candidates exist in some systems (that cannot be confirmed with this method at this time). A few of these candidates would also be near first order MMRs with either the confirmed planets or with other candidates. One system of particular interest, Kepler-56 (KOI-1241), is a pair of planets orbiting a 12th magnitude, giant star with radius over three times that of the Sun and effective temperature of 4900 K---among the largest stars known to host a transiting exoplanetary system.
• ### Production of oriented nitrogen-vacancy color centers in synthetic diamond(1112.5757)

Dec. 24, 2011 cond-mat.mtrl-sci
The negatively charged nitrogen-vacancy (NV-) center in diamond is an attractive candidate for applications that range from magnetometry to quantum information processing. Here we show that only a fraction of the nitrogen (typically < 0.5 %) incorporated during homoepitaxial diamond growth by Chemical Vapor Deposition (CVD) is in the form of undecorated NV- centers. Furthermore, studies on CVD diamond grown on (110) oriented substrates show a near 100% preferential orientation of NV- centers along only the [111] and [-1-11] directions, rather than the four possible orientations. The results indicate that NV centers grow in as units, as the diamond is deposited, rather than by migration and association of their components. The NV unit of the NVH- is similarly preferentially oriented, but it is not possible to determine whether this defect was formed by H capture at a preferentially aligned NV center or as a complete unit. Reducing the number of NV orientations from 4 orientations to 2 orientations should lead to increased optically-detected magnetic resonance contrast and thus improved magnetic sensitivity in ensemble-based magnetometry.
• ### A New Limit on Time-Reversal Violation in Beta Decay(1104.2778)

April 26, 2011 nucl-ex
We report the results of an improved determination of the triple correlation $D P \cdot(p_{e}\times p_{\nu})$ that can be used to limit possible time-reversal invariance in the beta decay of polarized neutrons and constrain extensions to the Standard Model. Our result is $D=(-0.96\pm 1.89 (stat)\pm 1.01 (sys))\times 10^{-4}$. The corresponding phase between g_A and g_V is $\phi_{AV} = 180.013^\circ\pm0.028^\circ$ (68 % confidence level). This result represents the most sensitive measurement of D in beta decay.
• ### Adjusting for Network Size and Composition Effects in Exponential-Family Random Graph Models(1004.5328)

Dec. 27, 2010 stat.ME
Exponential-family random graph models (ERGMs) provide a principled way to model and simulate features common in human social networks, such as propensities for homophily and friend-of-a-friend triad closure. We show that, without adjustment, ERGMs preserve density as network size increases. Density invariance is often not appropriate for social networks. We suggest a simple modification based on an offset which instead preserves the mean degree and accommodates changes in network composition asymptotically. We demonstrate that this approach allows ERGMs to be applied to the important situation of egocentrically sampled data. We analyze data from the National Health and Social Life Survey (NHSLS).
• ### Imaging of a Transitional Disk Gap in Reflected Light: Indications of Planet Formation Around the Young Solar Analog LkCa 15(1005.5162)

July 1, 2010 astro-ph.SR, astro-ph.EP
We present H- and Ks-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward-scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the spectral energy distributions (SED) of such systems, comprising an optically thick outer disk with an inner truncation radius of ~46 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary models, our high-contrast imagery imposes an upper limit of 21 Jupiter masses on companions at separations outside of 0.1" and of 13 Jupiter masses outside of 0.2". Thus, we find that a planetary system around LkCa 15 is the most likely explanation for the disk architecture.
• ### Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars. Detailed abundance analysis of OGLE-2008-BLG-209S(0903.3044)

March 17, 2009 astro-ph.GA, astro-ph.SR
AIMS. Our aims are twofold. First we aim to evaluate the robustness and accuracy of stellar parameters and detailed elemental abundances that can be derived from high-resolution spectroscopic observations of microlensed dwarf and subgiant stars. We then aim to use microlensed dwarf and subgiant stars to investigate the abundance structure and chemical evolution of the Milky Way Bulge. [ABRIDGED] METHODS. We present a detailed elemental abundance analysis of OGLE-2008-BLG-209S, the source star of a new microlensing event towards the Bulge, for which we obtained a high-resolution spectrum with the MIKE spectrograph on the Magellan Clay telescope. We have performed four different analyses of OGLE-2008-BLG-209S. [ABRIDGED] We have also re-analysed three previous microlensed dwarf stars OGLE-2006-BLG-265S, MOA-2006-BLG-099S, and OGLE-2007-BLG-349S with the same method. This homogeneous data set, although small, enables a direct comparison between the different stellar populations. RESULTS. We find that OGLE-2008-BLG-209S is a subgiant star that has a metallicity of [Fe/H] ~-0.33. It possesses [alpha/Fe] enhancements similar to what is found for Bulge giant stars at the same metallicity, and what also is found for nearby thick disc stars at the same metallicity. In contrast, the previous three microlensing dwarf stars have very high metallicities, [Fe/H]>+0.4, and more solar-like abundance ratios, i.e. [alpha/Fe]~0. The decrease in the [alpha/Fe] ratio with [Fe/H] is the typical signature of enrichment from low and intermediate mass stars. We furthermore find that the results for the four Bulge stars, in combination with results from studies of giant stars in the Bulge, seem to favour a secular formation scenario for the Bulge.
• ### Electromagnetic Field Creation During EWPT Nucleation With MSSM(hep-ph/0503247)

July 31, 2008 hep-ph
We derive the equations of motion for electroweak MSSM with a right-handed Stop, from which we derive the equations for the electromagnetic field that arises from bubble nucleation and collisions during the first order electroweak phase transition that can occur in this MSSM. Introducing an isospin ansatz we derive e.o.m. for the electrically charged W fields uncoupled from all other fields. These serve as the current for the Maxwell-like e.o.m. for the em field. The resulting electromagnetic field arising during EWPT bubble nucleation is found. This electromagnetic structure, along with that arising from bubble collisions, could seed galactic and extra-galactic magnetic fields.
• ### The bar pattern speed of NGC 4431(astro-ph/0703332)

March 13, 2007 astro-ph
We present surface photometry and stellar kinematics of NGC 4431, a barred dwarf galaxy in the Virgo cluster undergoing a tidal interaction with one of its neighbors, NGC 4436. We measured its bar pattern speed using the Tremaine-Weinberg method, and derived the ratio of the corotation radius, D_L, to the bar semi-major axis, a_B. We found D_L/a_B=0.6^{+1.2}_{-0.4} at 99% confidence level. Albeit with large uncertainty, the probability that the bar ends close to its corotation radius (i.e., 1.0 \leq D_L/a_B \leq 1.4) is about twice as likely as that the bar is much shorter than corotation radius (i.e., D_L/a_B > 1.4).
• ### High-resolution simulations of the final assembly of Earth-like planets 2: water delivery and planetary habitability(astro-ph/0510285)

Jan. 4, 2007 astro-ph
The water content and habitability of terrestrial planets are determined during their final assembly, from perhaps a hundred 1000-km "planetary embryos" and a swarm of billions of 1-10 km "planetesimals." During this process, we assume that water-rich material is accreted by terrestrial planets via impacts of water-rich bodies that originate in the outer asteroid region. We present analysis of water delivery and planetary habitability in five high-resolution simulations containing about ten times more particles than in previous simulations (Raymond et al 2006a, Icarus, 183, 265-282). These simulations formed 15 terrestrial planets from 0.4 to 2.6 Earth masses, including five planets in the habitable zone. Every planet from each simulation accreted at least the Earth's current water budget; most accreted several times that amount (assuming no impact depletion). Each planet accreted at least five water-rich embryos and planetesimals from past 2.5 AU; most accreted 10-20 water-rich bodies. We present a new model for water delivery to terrestrial planets in dynamically calm systems, with low-eccentricity or low-mass giant planets -- such systems may be very common in the Galaxy. We suggest that water is accreted in comparable amounts from a few planetary embryos in a "hit or miss" way and from millions of planetesimals in a statistically robust process. Variations in water content are likely to be caused by fluctuations in the number of water-rich embryos accreted, as well as from systematic effects such as planetary mass and location, and giant planet properties.
• ### High-resolution simulations of the final assembly of Earth-like planets 1: terrestrial accretion and dynamics(astro-ph/0510284)

Oct. 10, 2005 astro-ph
The final stage in the formation of terrestrial planets consists of the accumulation of ~1000-km planetary embryos'' and a swarm of billions of 1-10 km planetesimals.'' During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ~0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: 1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. 2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. 3) If planetary embryos form relatively slowly, following the models of Kokubo & Ida, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. 4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.