Source author record

Hagai B. Perets

Hagai B. Perets appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

82works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

82 published item(s)

preprint2022arXiv

3D Hydrodynamical Simulations of Helium-Ignited Double-degenerate White Dwarf Mergers

The origins of type Ia supernovae (SNe Ia) are still debated. Some of the leading scenarios involve a double detonation in double white dwarf (WD) systems. In these scenarios, helium shell detonation occurs on top of a carbon-oxygen (CO) WD, which then drives the detonation of the CO-core, producing a SN Ia. Extensive studies have been done on the possibility of a double helium detonation, following a dynamical helium mass-transfer phase onto a CO-WD. However, 3D self-consistent modeling of the double-WD system, the mass transfer, and the helium shell detonation have been little studied. Here we use 3D hydrodynamical simulations to explore this case in which a helium detonation occurs near the point of Roche lobe overflow of the donor WD and may lead to an SN Ia through the dynamically driven double-degenerate double-detonation (D6) mechanism. We find that the helium layer of the accreting primary WD does undergo a detonation, while the underlying carbon-oxygen core does not, leading to an extremely rapid and faint nova-like transient instead of a luminous SN Ia event. This failed core detonation suggests that D6 SNe Ia may be restricted to the most massive carbon-oxygen primary WDs. We highlight the nucleosynthesis of the long-lived radioisotope $^{44}$Ti during explosive helium burning, which may serve as a hallmark both of successful as well as failed D6 events which subsequently detonate as classical double-degenerate mergers.

preprint2022arXiv

A Measurement of Stellar Surface Gravity Hidden in Radial Velocity Differences of Co-moving Stars

The gravitational redshift induced by stellar surface gravity is notoriously difficult to measure for non-degenerate stars, since its amplitude is small in comparison with the typical Doppler shift induced by stellar radial velocity. In this study, we make use of the large observational data set of the Gaia mission to achieve a significant reduction of noise caused by these random stellar motions. By measuring the differences in velocities between the components of pairs of co-moving stars and wide binaries, we are able to statistically measure the combined effects of gravitational redshift and convective blueshifting of spectral lines, and nullify the effect of the peculiar motions of the stars. For the subset of stars considered in this study, we find a positive correlation between the observed differences in Gaia radial velocities and the differences in surface gravity and convective blueshift inferred from effective temperature and luminosity measurements. The results rule out a null signal at the $5σ$ level for our full data-set. Additionally, we study the sub-dominant effects of binary motion, and possible systematic errors in radial velocity measurements within Gaia. Results from the technique presented in this study are expected to improve significantly with data from the next Gaia data release. Such improvements could be used to constrain the mass-luminosity relation and stellar models which predict the magnitude of convective blueshift.

preprint2022arXiv

Are there any pristine comets? Constraints from pebble structure

We show that if comets (or any small icy planetesimals such as Kuiper belt objects) are composed of pebble piles, their internal radiogenic as well as geochemical heating results in considerably different evolutionary outcomes compared to similar past studies. We utilize a 1D thermo-physical evolution code, modified to include state-of-the-art empirical measurements of pebble thermal conductivity and compression, the latter obtained through a new laboratory experiment presented here for the first time. Results indicate that due to the low pebble thermal conductivity, the peak temperatures attained during evolution are much higher than in any previous study given the same formation time. Assuming meteoritic radiogenic abundances, we find that only extremely small, sub-kilometre comets have the potential to retain the primordial, uniform and thermally unprocessed composition from which they formed. Comets with radii in excess of about 20 km are typically swept by rapid and energetically powerful aqueous hydration reactions. Across the full range of comet sizes and formation times, evolutions result in the processing and differentiation of various volatile species, and a radially heterogeneous nucleus stucture. Our computations however also indicate that the assumed fraction of radionuclides is a pivotal free parameter, because isotopic analyses of the only available cometary samples suggest that no 26Al was ever present in comet 81P/Wild 2. We show that if comets formed early in the protoplanetary disc (within 1-3 Myr), the radionuclide abundances indeed must be much smaller than those typically assumed based on meteoritic samples. We discuss the importance of our findings for the formation, present-day attributes and future research of comets.

preprint2022arXiv

Binary evolution, gravitational-wave mergers and explosive transients in multiple-populations gas-enriched globular-clusters

Most globular clusters (GCs) show evidence for multiple stellar populations, suggesting the occurrence of several distinct star-formation episodes. The large fraction of second population (2P) stars observed requires a very large 2P gaseous mass to have accumulated in the cluster core to form these stars. Hence the first population of stars (1P) in the cluster core has had to become embedded in 2P gas, just prior to the formation of later populations. Here we explore the evolution of binaries in ambient 2P gaseous media of multiple-population GCs. We mostly focus on black hole binaries and follow their evolution as they evolve from wide binaries towards short periods through interaction with ambient gas, followed by gravitational-wave (GW) dominated inspiral and merger. We show this novel GW-merger channel could provide a major contribution to the production of GW-sources. We consider various assumptions and initial conditions and calculate the resulting gas-mediated change in the population of binaries and the expected merger rates due to gas-catalyzed GW-inspirals. For plausible conditions and assumptions, we find an expected GW merger rate observable by aLIGO of the order of up to a few tens of $\rm{Gpc^{-3} yr^{-1}}$, and an overall range for our various models of $0.08-25.51 \ \rm{Gpc^{-3} yr^{-1}}$. Finally, our results suggest that the conditions and binary properties in the early stage of GCs could be critically affected by gas-interactions and may require a major revision in the current modeling of the evolution of GCs.

preprint2022arXiv

Chaotic dynamics of wide triples induced by galactic tides: a novel channel for producing compact binaries, mergers, and collisions

Recent surveys show that wide ($>10^4$ AU) binaries and triples are abundant in the field. We study the long-term evolution of wide hierarchical triple systems and the role played by the Galactic tidal (GT) field. We find that when the timescales of the secular von-Ziepel-Lidov-Kozai and the GT oscillations are comparable, triple evolution becomes chaotic which leads to extreme eccentricities. Consequently, the close pericentre approaches of the inner-binary components lead to strong interactions, mergers and collisions. We use a novel secular evolution code to quantify the key parameters and carry out a population-synthesis study of low and intermediate-mass wide-orbit triples. We find that in $\sim9\%$ of low-mass wide-triples the inner main-sequence binaries collide or tidally-inspiral within $10\ \rm Gyr$, with direct collisions are $6$ times more likely to occur. For the intermediate-mass sample, $\sim7.6\%$ of the systems merge or inspiral with roughly equal probabilities. We discuss the relative fractions of different stellar merger/inspiral outcomes as a function of their evolutionary stage (Main-Sequence, MS; Red-Giant, RG; or White-Dwarf, WD), their transient electromagnetic signatures and the final products of the merger/inspiral. In particular, the rate of WD-WD direct-collisions that lead to type-Ia Supernovae is comparable to other dynamical channels and accounts for at most $0.1\%$ of the observed rate. RG inspirals provide a novel channel for the formation of eccentric common-envelope-evolution binaries. The catalysis of mergers/collisions in triples due to GT could explain a significant fraction, or even the vast majority, of blue-stragglers in the field, produce progenitors for cataclysmic-variables, and give-rise to mergers/collisions of double-RG binaries.

preprint2022arXiv

Compact objects formation, retention, and growth through accretion onto gas-embedded white-dwarfs/neutron-stars in gas-enriched globular-clusters

Observations of pulsars in globular clusters (GCs) give evidence that more >10-20% of neutron stars (NSs) ever formed in GCs were retained there. However, the velocity distribution of field pulsars peaks at 5-10 times the escape velocities of GCs. Consequently, only a small fraction of GC-NSs should have been retained, even accounting for low-velocity NSs formed through electron-capture supernovae. Thus, too few low-velocity NSs should have been retained in GCs, giving rise to the NS retention problem in GCs. Here we suggest a novel solution, in which the progenitors of most GC-NSs were ONe white-dwarfs (WDs) that accreted ambient intra-cluster gas and formed low-velocity NSs through accretion induced collapse (AIC). The existence of an early gas-enriched environment in GCs is supported by observations of GC multiple stellar populations. It is thought that 10s-100s of Myrs after the formation of the first generation of stars, and after ONe-WDs were already formed, GCs were replenished with gas which formed a second generation of stars. Accretion of such replenished gas onto the ONe-WDs catalyzed the AIC processes. The number of AIC-formed NSs is then sufficient to explain the large number of NSs retained in GCs. Similar processes might also drive CO-WDs to produce type-Ia SNe or to merge and form NSs, and similarly drive NSs to AIC and mergers producing BHs. Moreover, the wide variety of gas-catalyzed binary mergers and explosive transients suggest to occur in the gas-rich environments of AGN disk could similarly, and even more efficiently, occur in second-generation gas in GCs.

preprint2022arXiv

Constraints on the origins of hypervelocity stars: velocity distribution, mergers and star-formation history

In recent years surveys have identified several dozen B stars in the Milky Way halo moving faster than the local escape speed. The origin of most of these hypervelocity stars (HVSs) is still poorly constrained. Here we show that the velocity distribution, and in particular the deficiency in >700 km/s HVSs is inconsistent with binary disruptions by the massive black hole (MBH) in the Galactic Centre. This conclusion holds in the full and empty loss cone regime, and for secular instabilities in eccentric disks. Accounting for multiple close encounters between binaries and the MBH, does not qualitatively change the results. Moreover, there is no observed counterpart population in the Galactic Centre that is consistent with the HVSs. The star-formation history could be tuned explain the HVS velocity distribution, but this tuning would produce a mismatch with the observed HVS flight times. Frequent stellar collisions of the binary components due to interactions with the MBH do not significantly impact the velocity distribution in the Galactic halo. Such collisions, however, can leave observable remnants in the Galactic Centre, and potentially explain the origins of G2-like dust clouds.

preprint2022arXiv

Inflated Eccentric Migration of evolving gas giants II: Numerical methodology and basic concepts

Hot and Warm Jupiters (HJs&WJs) are gas-giant planets orbiting their host stars at short orbital periods, posing a challenge to their efficient in-situ formation. Therefore, most of the HJs&WJs are thought to have migrated from an initially farther-out birth locations. Current migration models, i.e disc-migration (gas-dissipation driven) and eccentric-migration (tidal evolution driven), fail to produce the occurrence rate and orbital properties of HJs&WJs. Here we study the role of the thermal evolution and its coupling to tidal evolution. We use the AMUSE, numerical environment, and MESA, planetary evolution modeling, to model in detail the coupled internal and orbital evolution of gas-giants during their eccentric-migration. In a companion paper, we use a simple semi-analytic model, validated by our numerical model, and run a population-synthesis study. We consider the initially inflated radii of gas-giants (expected following their formation), as well study the effects of the potential slowed contraction and even re-inflation of gas-giants (due to tidal and radiative heating) on the eccentric-migration. Tidal forces that drive eccentric-migration are highly sensitive to the planetary structure and radius. Consequently, we find that this form of inflated eccentric-migration operates on significantly (up to an order of magnitude) shorter timescales than previously studied eccentric-migration models. Thereby, inflated eccentric-migration gives rise to more rapid formation of HJs&WJs, higher occurrence rates of WJs, and higher rates of tidal disruptions, compared with previous eccentric migration models which consider constant ~Jupiter radii for HJ&WJ progenitors. Coupled thermal-dynamical evolution of eccentric gas-giants can therefore play a key-role in their evolution.

preprint2022arXiv

Inflated Eccentric Migration of Evolving Gas-Giants I: Accelerated Formation and Destruction of Hot and Warm Jupiters

Hot and warm Jupiters (HJs and WJs correspondingly) are gas-giants orbiting their host stars at very short orbital periods ($P_{HJ}<10$ days; $10<P_{WJ}<200$ days). HJs and a significant fraction of WJs are thought to have migrated from an initially farther-out birth locations. While such migration processes have been extensively studied, the thermal evolution of gas-giants and its coupling to the the migration processes are usually overlooked. In particular, gas-giants end their core-accretion phase with large radii and then contract slowly to their final radii. Moreover, intensive heating can slow the contraction at various evolutionary stages. The initial inflated large radii lead to faster tidal migration due to the strong dependence of tides on the radius. Here we explore this accelerated migration channel, which we term inflated eccentric migration, using a semi-analytical self-consistent modeling of the thermal-dynamical evolution of the migrating gas-giants, later validated by our numerical model (see a companion paper, paper II). We demonstrate our model for specific examples and carry a population synthesis study. Our results provide a general picture of the properties of the formed HJs\&WJs via inflated migration, and the dependence on the initial parameters/distributions. We show that tidal migration of gas-giants could occur far more rapidly then previously thought and lead to accelerated destruction and formation of HJs and enhanced formation rate of WJs. Accounting for the coupled thermal-dynamical evolution is therefore critical to the understanding of HJs/WJs formation, evolution and final properties of the population and play a key role in their migration process.

preprint2021arXiv

Extreme mass-ratio gravitational-wave sources: Mass segregation and post binary tidal-disruption captures

The gravitational-wave (GW) inspirals of stellar-mass compact objects onto a supermassive black hole (MBH), are some of the most promising GW sources detectable by next-generation space-born GW-detectors. The rates and characteristics of such extreme mass ratio inspirals (EMRIs) sources are highly uncertain. They are determined by the dynamics of stars near MBHs, and the rate at which compacts objects are driven to the close proximity of the MBH. Here we consider weakly and strongly mass-segregated nuclear clusters, and the evolution of stars captured into highly eccentric orbits following binary disruptions by the MBH. We make use of a Monte-Carlo approach to model the diffusion of both captured objects, and compact-objects brought through two-body relaxation processes. We calculate the rates of GW-inspirals resulting from relaxation-driven objects, and characterize EMRIs properties. We correct previous studies and show that relaxation-driven sources produce GW-sources with lower-eccentricity than previously found, and provide the detailed EMRI eccentricity distribution in the weak and strong mass-segregation regimes. We also show that binary-disruption captured-stars could introduce low-eccentricity GW-sources of stellar black-hole EMRIs in mass-segregated clusters. The eccentricities of the GW-sources from the capture channel, however, are strongly affected by relaxation processes, and are significantly higher than previously suggested. We find that both the rate and eccentricity distribution of EMRIs could probe the dynamics near MBHs, and the contribution of captured stars, characterize the mass-function of stellar compact objects, and verify whether weak or strong mass-segregation processes take place near MBHs.

preprint2021arXiv

Rapid destruction of planetary debris around white dwarfs through aeolian erosion

The discovery of numerous debris disks around white dwarfs (WDs), gave rise to extensive study of such disks and their role in polluting WDs, but the formation and evolution of these disks is not yet well understood. Here we study the role of aeolian (wind) erosion in the evolution of solids in WD debris disks. Aeolian erosion is a destructive process that plays a key role in shaping the properties and size-distribution of planetesimals, boulders and pebbles in gaseous protoplanetary disks. Our analysis of aeolian erosion in WD debris disks shows it can also play an important role in these environments. We study the effects of aeolian erosion under different conditions of the disk, and its erosive effect on planetesimals and boulders of different sizes. We find that solid bodies smaller than $\sim 5 \rm{km}$ will be eroded within the short disk lifetime. We compare the role of aeolian erosion in respect to other destructive processes such as collisional fragmentation and thermal ablation. We find that aeolian erosion is the dominant destructive process for objects with radius $\lesssim 10^3 \rm{cm}$ and at distances $\lesssim 0.6 R_\odot$ from the WD. Thereby, aeolian erosion constitutes the main destructive pathway linking fragmentational collisions operating on large objects with sublimation of the smallest objects and Poynting-Robertson drag, which leads to the accretion of the smallest particles onto the photosphere of WDs, and the production of polluted WDs.

preprint2021arXiv

Return of the TEDI: revisiting the Triple Evolution Dynamical Instability channel in triple stars

Triple-star systems exhibit a phenomenon known as the Triple Evolution Dynamical Instability (TEDI), in which mass loss in evolving triples triggers short-term dynamical instabilities, potentially leading to collisions of stars, exchanges, and ejections. Previous work has shown that the TEDI is an important pathway to head-on stellar collisions in the Galaxy, significantly exceeding the rate of collisions due to random encounters in globular clusters. Here, we revisit the TEDI evolutionary pathway using state-of-the-art population synthesis methods that self-consistently take into account stellar evolution and binary interactions, as well as gravitational dynamics and perturbations from passing stars in the field. We find Galactic TEDI-induced collision rates on the order of 1e-4/yr, consistent with previous studies which were based on more simplified methods. The majority of TEDI-induced collisions involve main sequence stars, potentially producing blue straggler stars. Collisions are also possible involving more evolved stars, potentially producing eccentric post-common-envelope systems, and white dwarfs collisions leading to Type Ia supernovae (although the latter with a negligible contribution to the Galactic rate). In our simulations, the TEDI is not only triggered by adiabatic wind mass loss, but also by Roche lobe overflow in the inner binary: when the donor star becomes less massive than the accretor, the inner binary orbit widens, triggering triple dynamical instability. We find that collision rates are increased by ~17% when fly-bys in the field are taken into account. In addition to collisions, we find that the TEDI produces ~1e-4/yr of unbound stars, although none with escape speeds in excess of 1e3 km/s.

preprint2021arXiv

Spatially Resolved X-ray Study of Supernova Remnant G306.3$-$0.9 with Unusually High Calcium Abundance

G306.3$-$0.9 is an asymmetric Galactic supernova remnant (SNR), whose progenitor has been thought to be a Type Ia supernova (SN), but its high Ca abundance appears inconsistent with the Type Ia origin. Hoping to uncover the reason for its asymmetry and the origin of this SNR, we performed a spatially resolved X-ray spectroscopic analysis of XMM-Newton and Chandra observation data. We divided the SNR into 13 regions and analyzed the spectra using two-temperature models (0.2 keV + 1 keV). Compared to the southwestern regions, the northeastern regions have higher metal abundances and a lower gas density. This suggests that the asymmetric morphology results from the non-uniform ambient environment. We found that neither Type Ia nor core-collapse SN models can account for the abnormally high abundance ratios of Ar/Si, Ca/Si, or the shape of the abundance curve. A comparison with the Ca-rich transient models shows that G306.3$-$0.9 is likely to be the first identified Galactic "Ca-rich transient" remnant, although the theoretical production of element S is lower. We also note that the conclusion for the SNR's origin relies on the measured abundance ratios and existing nucleosynthesis models. Between two groups of Ca-rich transient explosion models, we prefer the He shell detonation for an accreting WD, rather than the merger of a white dwarf and a neutron star.

preprint2021arXiv

The former companion of the hyper-velocity star S5-HVS1

The hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.2 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance in the range 1-10 AU and semimajor axis about $10^3$ AU. We then explore the fate of the captured star. We find that the heat deposited by tidally excited stellar oscillation modes leads to runaway disruption if the pericenter distance is smaller than about 3 AU. Over the past 5 Myr, its angular momentum has been significantly modified by orbital relaxation, which may stochastically drive the pericenter inwards below 3 AU and cause tidal disruption. We find an overall survival probability in the range 5% to 50%, depending on the local relaxation time in the close environment of the captured star, and the initial pericenter at capture. The pericenter distance of the surviving star has migrated to 10-100 AU, making it potentially the most extreme member of the S-star cluster. From the ejection rate of S5-HVS1-like stars, we estimate that there may currently be a few stars in such highly eccentric orbits. They should be detectable (typically Ks < 18.5 mag) by the GRAVITY instrument and by future Extremely Large Telescopes and hence provide an extraordinary probe of the spin of SgrA*.

preprint2020arXiv

A preceding low-virulence strain pandemic inducing immunity against COVID-19

Countries highly exposed to incoming traffic from China were expected to be at the highest risk of COVID-19 spread. However, COVID-19 case numbers (infection levels) are negatively correlated with incoming traffic-level. Moreover, infection levels are positively correlated with population-size, while the latter should only affect infection-level once herd immunity is reached. These could be explained if a low-virulence strain (LVS) began spreading a few months earlier from China, providing immunity from the later emerging known SARS-CoV-2 high-virulence strain (HVS). We find that the dynamics of the COVID-19 pandemic depend on the LVS and HVS spread doubling-times and the delay between their initial onsets. We find that LVS doubling-time to be $T_L\sim1.59\pm0.17$ times slower than the HVS ($T_H$), but its earlier onset allowed its global wide-spread to the levels required for herd-immunity. In countries exposed earlier to the LVS and/or having smaller population-size, the LVS achieved herd-immunity earlier, allowing less time for the spread of the HVS, and giving rise to lower HVS-infection levels. Such model accurately predicts a country's infection-level ({\rm R^{2}=0.74}; p-value of {\rm 5.2\times10^{-13}}), given only its population-size and incoming-traffic from China. It explains the negative correlation with incoming-traffic ($c_{exp}$), the positive correlation with the population size (n_{pop}) and their specific relations (${\rm N}_{\rm cases}\propto n_{pop}^{{\rm T_{L}/{\rm T_{H}}}}\times c_{exp}^{{\rm T_{L}/{\rm T_{H}-1}}}$). We find that most countries should have already achieved herd-immunity. Further COVID-19-spread in these countries is limited and is not expected to rise by more than a factor of 2-3. We suggest tests/predictions to further verify the model and biologically identify the LVS, and discuss the implications.

preprint2020arXiv

Collisional formation of massive exomoons of super-terrestrial exoplanets

Exomoons orbiting terrestrial or super-terrestrial exoplanets have not yet been discovered; their possible existence and properties are therefore still an unresolved question. Here we explore the collisional formation of exomoons through giant planetary impacts. We make use of smooth particle hydrodynamical (SPH) collision simulations and survey a large phase-space of terrestrial/super-terrestrial planetary collisions. We characterize the properties of such collisions, finding one rare case in which an exomoon forms through a graze&capture scenario, in addition to a few graze&merge or hit&run scenarios. Typically however, our collisions form massive circumplanetary discs, for which we use follow-up N-body simulations in order to derive lower-limit mass estimates for the ensuing exomoons. We investigate the mass, long-term tidal-stability, composition and origin of material in both the discs and the exomoons. Our giant-impact models often generate relatively iron-rich moons, that form beyond the synchronous radius of the planet, and would thus tidally evolve outward with stable orbits, rather than be destroyed. Our results suggest that it is extremely difficult to collisionally form currently-detectable exomoons orbiting super-terrestrial planets, through single giant impacts. It might be possible to form massive, detectable exomoons through several mergers of smaller exomoons, formed by multiple impacts, however more studies are required in order to reach a conclusion. Given the current observational initiatives, the search should focus primarily on more massive planet categories. However, about a quarter of the exomoons predicted by our models are approximately Mercury-mass or more, and are much more likely to be detectable given a factor 2 improvement in the detection capability of future instruments, providing further motivation for their development.

preprint2020arXiv

Erosion driven size-redistribution of protoplanetary disk solids and the onset of streaming Instability and Pebble Accretion

The formation of the first planetesimals and the final growth of planetary cores relies on the abundance of small pebbles. The efficiencies of both the streaming instability (SI) process, suggested to catalyze the early growth of planetesimals, and the pebble-accretion process, suggested to accelerate the growth of planetary cores, depend on the sizes of solids residing in the disk. In particular, these processes were found to be sensitive to size distribution of solids, and efficient planetesimal formation and growth through these channels require a limited pebble size distribution. Here we show that aeolian erosion, a process that efficiently grinds down boulders into a mono-sized distribution of pebbles, provides a natural upper limit for the maximal pebble sizes (in terms of their Stokes number). We find the dependence of this upper limit on the radial separation, disk age, turbulence strength, and the grain-size composition of the boulders in the disk. SI is favorable in areas with a Stokes number less than 0.1, which is found in the inner sub-astronomical-unit regions of the disk. This upper limit shapes the size distribution of small pebbles and thereby catalyzes the early onset of planetesimal formation due to SI, and the later core accretion growth through pebble accretion.

preprint2020arXiv

Faint rapid red transients from Neutron star -- CO white-dwarf mergers

Mergers of neutron stars (NS) and white dwarfs (WD) may give rise to observable explosive transient events. We use 3D hydrodynamical (SPH) simulations, as well as 2D hydrodynamical-thermonuclear simulations (using the FLASH AMR code) to model the disruption of CO-WDs by NSs, which produce faint transient events. We post-process the simulations using a large nuclear network and make use of the SuperNu radiation-transfer code to predict the observational signatures and detailed properties of these transients. We calculate the light-curves (LC) and spectra for five models of NS - CO-WD mergers. The small yields of Ni56 (few x 0.001Msun) results in faint, rapidly-evolving reddened transients (RRTs) with B (R) - peak magnitudes of ~ -12 (-13) to ~ -13 (-15), much shorter and fainter than both regular and faint/peculiar type-Ia SNe. We show that the spectra of RRTs share some similarities with rapidly - evolving transients such as SN2010x, though RRTs are significantly fainter, especially in the I/R bands, and show far stronger Si lines. We estimate that the upcoming Large Synoptic Survey Telescope could detect RRTs at a rate of ~ 10 - 70 yr^-1, through observations in the R/I bands.

preprint2020arXiv

Inferred timescales for common envelope ejection using wide astrometric companions

Evolution of close binaries often proceeds through the common envelope stage. The physics of the envelope ejection (CEE) is not yet understood, and several mechanisms were suggested to be involved. These could give rise to different timescales for the CEE mass-loss. In order to probe the CEE-timescales we study wide companions to post-CE binaries. Faster mass-loss timescales give rise to higher disruption rates of wide binaries and result in larger average separations. We make use of data from Gaia DR2 to search for ultra-wide companions (projected separations $10^3$-$2\times 10^5$ a.u. and $M_2 > 0.4$ M$_\odot$) to several types of post-CEE systems, including sdBs, white-dwarf post-common binaries, and cataclysmic variables. We find a (wide-orbit) multiplicity fraction of $1.4\pm 0.2$ per cent for sdBs to be compared with a multiplicity fraction of $5.0\pm 0.2$ per cent for late-B/A/F stars which are possible sdB progenitors. The distribution of projected separations of ultra-wide pairs to main sequence stars and sdBs differs significantly and is compatible with prompt mass loss (upper limit on common envelope ejection timescale of $10^2$ years). The smaller statistics of ultra-wide companions to cataclysmic variables and post-CEE binaries provide weaker constraints. Nevertheless, the survival rate of ultra-wide pairs to the cataclysmic variables suggest much longer, $\sim10^4$ years timescales for the CEE in these systems, possibly suggesting non-dynamical CEE in this regime.

preprint2020arXiv

Origin of (2014) MU69-like Kuiper-belt contact binaries from wide binaries

Following its flyby and first imaging the Pluto-Charon binary, the New Horizons spacecraft visited the Kuiper-Belt-Object (KBO) (486958) 2014 MU69 (Arrokoth). Imaging showed MU69 to be a contact-binary, made of two individual lobes connected by a narrow neck, rotating at low spin period (15.92 h), and having high obliquity (~98 deg), similar to other KBO contact-binaries inferred through photometric observations. The origin of such peculiar configurations is puzzling, and all scenarios suggested for the origins of contact-binaries fail to reproduce such properties and their likely high frequency. Here we show that semi-secular perturbations operating only on ultra-wide (~0.1-0.4 Hill-radius) KBO-binaries can robustly lead to gentle, slow-speed binary mergers at arbitrarily high obliquities, but low rotational velocities, that can reproduce MU69's (and similar oblique contact binaries) characteristics. Using N-body simulations, we find that ~15% of all ultra-wide binaries with cosine-uniform inclination distribution are likely to merge through this process. Moreover, we find that such mergers are sufficiently gentle as to only slightly deform the KBO shape, and can produce the measured rotation speed of MU69. The semi-secular contact-binary formation channel not only explains the observed properties of MU69, but could also apply for other Kuiper/asteroid belt binaries, and for Solar/extra-solar moon systems.

preprint2020arXiv

Rates of Stellar Tidal Disruption

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

preprint2020arXiv

The aeolian-erosion barrier for the growth of metre-size objects in protoplanetary-discs

Aeolian-erosion is a destructive process which can erode small-size planetary objects through their interaction with a gaseous environment. Aeolian-erosion operates in a wide range of environments and under various conditions. Aeolian-erosion has been extensively explored in the context of geophysics in terrestrial planets. Here we show that aeolian-erosion of cobbles, boulders and small planetesimals in protoplanetary-discs can constitute a significant barrier for the early stages of planet formation. We use analytic calculations to show that under the conditions prevailing in protoplanetary-discs small bodies ($10-10^4 \rm{m}$) are highly susceptible to gas-drag aeolian-erosion. At this size-range aeolian-erosion can efficiently erode the planetesimals down to tens-cm size and quench any further growth of such small bodies. It thereby raises potential difficulties for channels suggested to alleviate the metre-size barrier. Nevertheless, the population of $\sim$decimetre-size cobbles resulting from aeolian-erosion might boost the growth of larger (>km size) planetesimals and planetary embryos through increasing the efficiency of pebble-accretion, once/if such large planetesimals and planetary embryos exist in the disc.

preprint2020arXiv

The Wide-Binary Origin of The Pluto-Charon System

The Pluto-Charon binary system is the best-studied representative of the binary Kuiper-belt population. Its origins are vital to understanding the formation of other Kupier-belt objects (KBO) and binaries, and the evolution of the outer solar-system. The Pluto-Charon system is believed to form following a giant impact between two massive KBOs at relatively low velocities. However, the likelihood of a random direct collision between two of the most massive KBOs is low, and is further constrained by the requirement of a low-velocity collision, making this a potentially fine-tuned scenario. Here we expand our previous studies and suggest that the proto-Pluto-Charon system was formed as a highly inclined wide-binary, which was then driven through secular/quasi-secular evolution into a direct impact. Since wide-binaries are ubiquitous in the Kuiper-belt with many expected to be highly inclined, our scenario is expected to be robust. We use analytic tools and few-body simulations of the triple Sun-(proto-)Pluto-Charon system to show that a large parameter-space of initial conditions leads to such collisions. The velocity of such an impact is the escape velocity of a bound system, which naturally explains the low-velocity impact. The dynamical evolution and the origins of the Pluto-Charon system could therefore be traced to similar secular origins as those of other binaries and contact-binaries (e.g. Arrokoth), and suggest they play a key role in the evolution of KBOs.

preprint2019arXiv

Gravitational waves from in-spirals of compact objects in binary common-envelope evolution

Detection of gravitational-wave (GW) sources enables the characterisation of binary compact objects and of their in-spiral. However, other dissipative processes can affect the in-spiral. Here we show that the in-spiral of compact objects through a gaseous common-envelope (CE) arising from an evolved stellar companion produces a novel type of GW-sources, whose evolution is dominated by the dissipative gas dynamical friction effects from the CE, rather than the GW-emission itself. The evolution and properties of the GW-signals differ from those of isolated gas-poor mergers significantly. We find characteristic strains of $\sim10^{-23}$-$10^{-21}$ ($10{\rm kpc}/{D}$) for such sources -- observable by next-generation space-based GW-detectors. The evolution of the GW-signal can serve as a probe of the interior regions of the evolved star, and the final stages of CE-evolution, otherwise inaccessible through other observational means. Moreover, such CE-mergers are frequently followed by observable explosive electromagnetic counterparts and/or the formation of exotic stars.

preprint2019arXiv

Gravitational-wave sources from mergers of binary black-holes catalyzed by fly-bys interactions in the field

Several scenarios were suggested for the origins of gravitational-wave (GW) sources from mergers of stellar binary black holes (BBHs). Here we propose a novel origin through catalyzed formation of GW-sources from ultra-wide binaries in the field. Such binaries experience perturbations from random stellar fly-bys which excite their eccentricities. Once a wide-binary is driven to a sufficiently small peri-center approach, GW-emission becomes significant, and the binary inspirals and merges. We derive an analytic model and verify it with numerical calculation to compute the merger rate to be $\sim10{\rm \times f_{wide}\,{\rm Gpc}^{-3}yr^{-1}}$ ($f_{\rm wide}$ is the fraction of wide BH-binaries), which is comparable to the observationally inferred rate. The observational signatures from this channel include spin-orbit misalignment; preference for high mass-ratio BBH; preference for high velocity-dispersion host-galaxies; and a uniform delay-time distribution.

preprint2016arXiv

A triple origin for the lack of tight coplanar circumbinary planets around short-period binaries

Transiting circumbinary planets are more easily detected around short-period than long-period binaries, but none have yet been observed by {\it Kepler} orbiting binaries with periods shorter than seven days. In triple systems, secular Kozai-Lidov cycles and tidal friction (KLCTF) have been shown to reduce the inner orbital period from $\sim 10^4$ to a few days. Indeed, the majority of short-period binaries are observed to possess a third stellar companion. Using secular evolution analysis and population synthesis, we show that KLCTF makes it unlikely for circumbinary transiting planets to exist around short-period binaries. We find the following outcomes. (1) Sufficiently massive planets in tight and/or coplanar orbits around the inner binary can quench the KL evolution because they induce precession in the inner binary. The KLCTF process does not take place, preventing the formation of a short-period binary. (2) Secular evolution is not quenched and it drives the planetary orbit into a high eccentricity, giving rise to an unstable configuration, in which the planet is most likely ejected from the system. (3) Secular evolution is not quenched but the planet survives the KLCTF evolution. Its orbit is likely to be much wider than the currently observed inner binary orbit, and is likely to be eccentric and inclined with respect to the inner binary. These outcomes lead to two main conclusions: (1) it is unlikely to find a massive planet on a tight and coplanar orbit around a short-period binary, and (2) the properties of circumbinary planets in short-period binaries are constrained by secular evolution.

preprint2016arXiv

Application of Gas Dynamical Friction for Planetesimals: II. Evolution of Binary Planetesimals

One of first the stages of planet formation is the growth of small planetesimals and their accumulation into large planetesimals and planetary embryos. This early stage occurs much before the dispersal of most of the gas from the protoplanetary disk. At this stage gas-planetesimal interactions play a key role in the dynamical evolution of \emph{single} intermediate-mass planetesimals ($m_{p}\sim10^{21}-10^{25}g$) \emph{through gas dynamical friction} (GDF). A significant fraction of all Solar system planetesimals (asteroids and Kuiper-belt objects) are known to be binary planetesimals (BPs). Here, we explore the effects of GDF on the evolution of \emph{binary} planetesimals embedded in a gaseous disk using an N-body code with a fiducial external force accounting for GDF. We find that GDF can induce binary mergers on timescales shorter than the disk lifetime for masses above $m_{p}\gtrsim10^{22}g$ at 1AU, independent of the binary initial separation and eccentricity. Such mergers can affect the structure of merger-formed planetesimals, and the GDF-induced binary inspiral can play a role in the evolution of the planetesimal disk. In addition, binaries on eccentric orbits around the star may evolve in the supersonic regime, where the torque reverses and the binary expands, which would enhance the cross section for planetesimal encounters with the binary. Highly inclined binaries with small mass ratios, evolve due to the combined effects of Kozai-Lidov cycles with GDF which lead to chaotic evolution. Prograde binaries go through semi-regular Kozai-Lidov evolution, while retrograde binaries frequently flip their inclination and $\sim50\%$ of them are destroyed.

preprint2016arXiv

Generalized Hill-Stability Criteria for Hierarchical Three-Body Systems at Arbitrary Inclinations

A fundamental aspect of the three-body problem is its stability. Most stability studies have focused on the co-planar three-body problem, deriving analytic criteria for the dynamical stability of such pro/retrograde systems. Numerical studies of inclined systems phenomenologically mapped their stability regions, but neither complement it by theoretical framework, nor provided satisfactory fit for their dependence on mutual inclinations. Here we present a novel approach to study the stability of hierarchical three-body systems at arbitrary inclinations, which accounts not only for the instantaneous stability of such systems, but also for the secular stability and evolution through Lidov-Kozai cycles and evection. We generalize the Hill-stability criteria to arbitrarily inclined triple systems, explain the existence of quasi-stable regimes and characterize the inclination dependence of their stability. We complement the analytic treatment with an extensive numerical study, to test our analytic results. We find excellent correspondence up to high inclinations $(\sim120^{\circ}$), beyond which the agreement is marginal. At such high inclinations the stability radius is larger, the ratio between the outer and inner periods becomes comparable, and our secular averaging approach is no longer strictly valid. We therefore combine our analytic results with polynomial fits to the numerical results to obtain a generalized stability formula for triple systems at arbitrary inclinations. Besides providing a generalized secular-based physical explanation for the stability of non co-planar systems, our results have direct implications for any triple systems, and in particular binary planets and moon/satellite systems; we briefly discuss the latter as a test case for our models.

preprint2016arXiv

Micro - tidal disruption events by stellar compact objects and the production of ultra-long GRBs

We explore full/partial tidal disruption events (TDEs) of stars/planets by stellar compact objects (Black holes; BHs; or neutron stars; NSs), which we term micro-TDEs. Disruption of a star/planet with mass $M_{\star}$ may lead to the formation of a debris disk around the BH/NS. Efficient accretion of a fraction $(f_{acc}=0.1$ of the debris may then give rise to bright energetic long ($10^{3}-10^{4}\,s$), X-ray/Gamma-ray flares, with total energies of up to $(f_{acc}/0.1)\times10^{52}\,(M_{\star}/{0.6\,M_{\odot}})\,$ergs, possibly resembling ultra-long GRBs/XRFs. The energy of such flares depends on the poorly constrained accretion processes. Significantly fainter flares might be produced if most of the disk mass is blown away through strong outflows. We suggest three dynamical origins for such disruptions. In the first, a star/planet is tidally disrupted following a close random encounter with a BH/NS in a dense cluster. We estimate the BH (NS) micro-TDE rates from this scenario to be few$\times10^{-6}$ (few$\times10^{-7}$) ${\rm yr}^{-1}$ per Milky-Way galaxy. Another scenario involves the interaction of wide companions due to perturbations by stars in the field, likely producing comparable but lower rates. Finally, a third scenario involves a BH/NS which gain a natal velocity kick at birth, leading to a close encounter with a binary companion and the tidal disruption of that companion. Such events could be associated with a supernova, or even with a preceding GRB/XRF event, and would likely occur hours to days after the prompt explosion; the rates of such events could be larger than those obtained from the other scenarios, depending on the preceding complex binary stellar evolution.

preprint2016arXiv

Neutron star natal kicks: Collisions, $μ$TDEs, faint SNe, GRBs and GW sources with preceding electromagnetic counterparts

Based on the observed high velocity of pulsars it is thought that neutron stars (NSs) receive a significant velocity kick at birth. Such natal kicks are considered to play an important role in the the evolution of binary-NS systems. The kick given to the NS (together with the effect of mass loss due to the supernova explosion of the NS progenitor) may result in the binary disruption or lead to a significant change of the binary orbital properties. Here we explore in detail the dynamical aftermath of natal kicks in binary systems, determine their possible outcomes and characterize their relative frequency, making use of analytic arguments and detailed population synthesis models. In a fraction of the cases the kick may cast the NS in such a trajectory as to collide with the binary companion, or pass sufficiently close to it as to disrupt it (micro tidal disruption event; $μ$TDE), or alternatively it could be tidally-captured into a close orbit, eventually forming an X-ray binary. We calculate the rates of direct post-kick physical collisions and the possible potential production of Thorne-Zytkow objects or long-GRBs through this process, estimate the rates X-ray binaries formation and determine the rates of $μ$TDEs and faint supernovae from white dwarf disruptions by NSs. Finally we suggest that natal kicks can produce BH-NS binaries with very short gravitational-wave merger time, possibly giving rise to a new type of promptly appearing eLISA gravitational wave (GW) sources, as well as producing aLIGO binary-merger GW sources with a unique (likely type Ib/c) supernova electromagnetic counterpart which precedes the GW merger.

preprint2016arXiv

Second generation stellar disks in Dense Star Clusters and cluster ellipticities

Globular Clusters (GCs) and Nuclear Star Clusters (NSCs) are typically composed by several stellar populations, characterized by different chemical compositions. Different populations show different ages in NSCs but not necessarily in GCs. The youngest populations in NSCs appear to reside in disk-like structures, as observed in our Galaxy and in M31. Gas infall followed by formation of second generation (SG) stars in GCs may similarly form disk-like structures in the clusters nuclei. Here we explore this possibility and follow the long term evolution of stellar disks embedded in GCs, and study their effects on the evolution of the clusters. We study disks with different masses by means of detailed N-body simulations and explore their morphological and kinematic signatures on the GC structures. We find that as a second generation disk relaxes, the old, first generation, stellar population flattens and becomes more radially anisotropic, making the GC structure become more elliptical. The second generation stellar population is characterized by a lower velocity dispersion, and a higher rotational velocity, compared with the primordial older population. The strength of these kinematic signatures depends both on the relaxation time of the system and on the fractional mass of the second generation disk. We therefore conclude that SG populations formed in flattened configurations will give rise to two systematic trends: (1) Positive correlation between GC ellipticity and fraction of SG population (2) Positive correlation between GC relaxation time and ellipticity. Thereby GC ellipticities and rotation could be related to the formation of SG stars and their initial configuration.

preprint2016arXiv

Secular dynamics of hierarchical quadruple systems: the case of a triple system orbited by a fourth body

We study the secular gravitational dynamics of quadruple systems consisting of a hierarchical triple system orbited by a fourth body. These systems can be decomposed into three binary systems with increasing semimajor axes, binaries A, B and C. The Hamiltonian of the system is expanded in ratios of the three binary separations, and orbit-averaged. Subsequently, we numerically solve the equations of motion. We study highly hierarchical systems that are well described by the lowest-order terms in the Hamiltonian. We find that the qualitative behaviour is determined by the ratio $\mathcal{R}_0$ of the initial Kozai-Lidov (KL) time-scales of the binary pairs AB and BC. If $\mathcal{R}_0\ll 1$, binaries AB remain coplanar if this is initially the case, and KL eccentricity oscillations in binary B are efficiently quenched. If $\mathcal{R}_0\gg 1$, binaries AB become inclined, even if initially coplanar. However, there are no induced KL eccentricity oscillations in binary A. Lastly, if $\mathcal{R}_0\sim 1$, complex KL eccentricity oscillations can occur in binary A that are coupled with the KL eccentricity oscillations in B. Even if binaries A and B are initially coplanar, the induced inclination can result in very high eccentricity oscillations in binary A. These extreme eccentricities could have significant implications for strong interactions such as tidal interactions, gravitational wave dissipation, and collisions and mergers of stars and compact objects. As an example, we apply our results to a planet+moon system orbiting a central star, which in turn is orbited by a distant and inclined stellar companion or planet, and to observed stellar quadruples.

preprint2016arXiv

Secular dynamics of multiplanet systems: implications for the formation of hot and warm Jupiters via high-eccentricity migration

Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.

preprint2016arXiv

Supernovae from direct collisions of white dwarfs and the role of helium shell ignition

Models for supernovae (SNe) arising from thermonuclear explosions of white dwarfs (WDs) have been extensively studied over the last few decades, mostly focusing on the single degenerate (accretion of material of a WD) and double degenerate (WD-WD merger) scenarios. In recent years it was suggested that WD-WD direct collisions provide an additional channel for such explosions. Here we extend the studies of such explosions, and explore the role of Helium-shells in affecting the thermonuclear explosions. We study both the impact of low-mass helium ($\sim0.01$ M$_{\odot})$ shells, as well as high mass shells ($\ge0.1$ M$_{\odot}$). We find that detonation of the massive helium layers precede the detonation of the WD Carbon-Oxygen (CO) bulk during the collision and can change the explosive evolution and outcomes for the cases of high mass He-shells. In particular, the He-shell detonation propagates on the WD surface and inefficiently burns material prior to the CO detonation that later follows in the central parts of the WD. Such evolution leads to larger production of intermediate elements, producing large yields of $^{44}{\rm Ti}$ and $^{48}{\rm Cr}$ relative to the pure CO-CO WD collisions. Collisions of WDs with a low-mass He-shell do not give rise to helium detonation, but helium burning does precede the CO bulk detonation. Such collisions produce a high velocity, low-mass of ejected burned material enriched with intermediate elements, with smaller changes to the overall explosion outcomes. The various effects arising from the contribution of low/high mass He layers change the kinematics and the morphological structure of collision-induced SNe and may thereby provide unique observational signatures for such SNe, and play a role in the chemical enrichment of galaxies and the production of intermediate elements and positrons from their longer-term decay.

preprint2016arXiv

The History of Tidal Disruption Events in Galactic Nuclei

The tidal disruption of a star by a massive black hole (MBH) is thought to produce a transient luminous event. Such tidal disruption events (TDEs) may play an important role in the detection and characterization of MBHs and probe the properties and dynamics of their nuclear stellar clusters (NSCs) hosts. Previous studies estimated the recent rates of TDEs in the local universe. However, the long-term evolution of the rates throughout the history of the universe has been hardly explored. Here we consider the TDE history, using evolutionary models for the evolution of galactic nuclei. We use a 1D Fokker-Planck approach to explore the evolution of MBH-hosting NSCs, and obtain the disruption rates of stars during their evolution. We complement these with an analysis of TDEs history based on N-body simulation data, and find them to be comparable. We consider NSCs that are built-up from close-in star-formation (SF) or from SF/clusters-dispersal far-out, a few pc from the MBH. We also explore cases where primordial NSCs exist and later evolve through additional star-formation/cluster-dispersal processes. We study the dependence of the TDE history on the type of galaxy, as well as the dependence on the MBH mass. These provide several scenarios, with a continuous increase of the TDE rates over time for cases of far-out SF and a more complex behavior for the close-in SF cases. Finally, we integrate the TDE histories of the various scenarios to provide a total TDE history of the universe, which can be potentially probed with future large surveys (e.g. LSST).

preprint2016arXiv

The impact of mass segregation and star-formation on the rates of gravitational-wave sources from extreme mass ratio inspirals

Compact stellar objects inspiralling into massive black holes (MBHs) in galactic nuclei are some of the most promising gravitational wave (GWs) sources for next generation GW-detectors. The rates of such extreme mass ratio inspirals (EMRIs) depend on the dynamics and distribution of compact objects around the MBH. Here we study the impact of mass-segregation processes on EMRI rates. In particular, we provide the expected mass function of EMRIs, given an initial mass function of stellar BHs (SBHs), and relate it to the mass-dependent detection rate of EMRIs. We then consider the role of star formation on the distribution of compact objects and its implication on EMRI rates. We find that the existence of a wide spectrum of SBH masses lead to the overall increase of EMRI rates, and to high rates of the EMRIs from the most-massive SBHs. However, it also leads to a relative quenching of EMRI rates from lower-mass SBHs, and together produces a steep dependence of the EMRI mass function on the highest-mass SBHs. Star-formation history plays a relatively small role in determining the EMRI rates of SBHs, since most of them migrate close to the MBH through mass-segregation rather than forming in-situ. However, the EMRI rate of neutron stars can be significantly increased when they form in-situ close to the MBH, as they can inspiral before relaxation processes significantly segregates them outwards. A reverse but weaker effect of decreasing the EMRI rates from neutron stars and white dwarfs occurs when star-formation proceeds far from the MBH.

preprint2015arXiv

A primordial origin for the composition similarity between the Earth and the Moon

Most of the properties of the Earth-Moon system can be explained by a collision between a planetary embryo and the growing Earth late in the accretion process. Simulations show that most of the material that eventually aggregates to form the Moon originates from the impactor. However, analysis of the terrestrial and lunar isotopic composition show them to be highly similar. In contrast, the compositions of other solar system bodies are significantly different than the Earth and Moon. This poses a major challenge to the giant impact scenario since the Moon-forming impactor is then thought to also have differed in composition from the proto-Earth. Here we track the feeding zones of growing planets in a suite of simulations of planetary accretion, in order to measure the composition of Moon-forming impactors. We find that different planets formed in the same simulation have distinct compositions, but the compositions of giant impactors are systematically more similar to the planets they impact. A significant fraction of planet-impactor pairs have virtually identical compositions. Thus, the similarity in composition between the Earth and Moon could be a natural consequence of a late giant impact.

preprint2015arXiv

Application of gas dynamical friction for planetesimals: I. Evolution of single planetesimals

The growth of small planetesimals into large planetary embryos occurs much before the dispersal of the gas from the protoplanetary disk. The planetesimal - gaseous-disk interactions give rise to migration and orbital evolution of the planetesimals/planets. Small planetesimals are dominated by aerodynamic gas drag. Large protoplanets, $m\sim0.1M_{\oplus}$, are dominated by type I migration \emph{differential} torque. There is an additional mass range, $m\sim10^{21}-10^{25}g$ of \emph{intermediate mass} planetesimals (IMPs), where gravitational interactions with the disk dominate over aerodynamic gas drag, but for which such interactions were typically neglected. Here we model these interactions using the \emph{gas dynamical friction} (GDF) approach, previously used to study the disk-planet interactions at the planetary mass range. We find the critical size where GDF dominates over gas drag, and then study the implications of GDF on single IMPs. We find that planetesimals with small inclinations rapidly become co-planar. Eccentric orbits circularize within a few Myrs, provided the the planetesimal mass is large, $m\gtrsim10^{23}g$ and that the initial eccentricity is low, $e\lesssim0.1$. Planetesimals of higher masses, $m\sim10^{24}-10^{25}g$ inspiral on a time-scale of a few Myrs, leading to \emph{an embryonic migration} to the inner disk. This can lead to an over-abundance of rocky material (in the form of IMPs) in the inner protoplanetary disk ($<1$AU) and induce rapid planetary growth. This can explain the origin of super-Earth planets. In addition, GDF damps the velocities of IMPs, thereby cooling the planetesimal disk and affecting its collisional evolution through quenching the effects of viscous stirring by the large bodies.

preprint2015arXiv

Photometric Amplitude Distribution of Stellar Rotation of Kepler KOIs-Indication for Spin-Orbit Alignment of Cool Stars and High Obliquity for Hot Stars

The observed amplitude of the rotational photometric modulation of a star with spots should depend on the inclination of its rotational axis relative to our line of sight. Therefore, the distribution of observed rotational amplitudes of a large sample of stars depends on the distribution of their projected axes of rotation. Thus, comparison of the stellar rotational amplitudes of the Kepler KOIs with those of Kepler single stars can provide a measure to indirectly infer the properties of the spin-orbit obliquity of Kepler planets. We apply this technique to the large samples of 993 KOIs and 33,614 single Kepler stars in temperature range of 3500-6500 K. We find with high significance that the amplitudes of cool KOIs are larger, on the order of 10%, than those of the single stars. In contrast, the amplitudes of hot KOIs are systematically lower. After correcting for an observational bias, we estimate that the amplitudes of the hot KOIs are smaller than the single stars by about the same factor of 10%. The border line between the relatively larger and smaller amplitudes, relative to the amplitudes of the single stars, occurs at about 6000K. Our results suggest that the cool stars have their planets aligned with their stellar rotation, while the planets around hot stars have large obliquities, consistent with the findings of Winn et al. (2010) and Albrecht et al. (2012). We show that the low obliquity of the planets around cool stars extends up to at least 50 days, a feature that is not expected in the framework of a model that assumes the low obliquity is due to planet-star tidal realignment.

preprint2015arXiv

Short and long term evolution of a stellar disk around a massive black hole: The role of binaries, the cusp and stellar evolution

We study the dynamical evolution of a stellar disk orbiting a massive black hole. We explore the role of two-body relaxation, mass segregation, stellar evolution and binary heating in affecting the disk evolution, and consider the impact of the nuclear cluster structure and the stellar-disk mass-function. We use analytic arguments and numerical calculations, and apply them to study the evolution of a stellar disk (similar to that observed in the Galactic center; GC), both on the short (few Myr) and longer (100 Myr) evolutionary timescales. We find the dominant processes affecting the disk evolution are two-body relaxation and mass segregation where as binary heating have only a little contribution. Massive stars play a dominant role in kinematically heating low mass stars, and driving them to high eccentricities/inclinations. Multi-mass models with realistic mass-functions for the disk stars show the disk structure to be mass stratified, with the most massive stars residing in thinner structures. Stellar evolution plays an important role in decreasing the number of massive stars with time, thereby leading to slower relaxation, where the remnant compact objects of these stars are excited to higher eccentricities/inclinations. At these later evolutionary stages dynamical heating by the nuclear cluster plays a progressively more important role. We conclude that the high eccentricities of the disk-stars in the Galactic Center suggest that the disk formed with initially high eccentricities, or that collective or secular processes dominate the disk evolution. Finally, we find that the disk structure is expected to keep a thin structure even after 100 Myrs. It therefore suggests earlier disks now containing only older, lower mass stars might still be observed in the Galactic center, unless destroyed/smeared by other non-two-body relaxation processes.

preprint2015arXiv

Tidal capture formation of Low Mass X-Ray Binaries from wide binaries in the field

We present a potentially efficient dynamical formation scenario for Low Mass X-ray Binaries (LMXBs) in the field, focusing on black-hole (BH) LMXBs. In this formation channel LMXBs are formed from wide binaries $(>1000$ AU) with a BH component and a stellar companion. The wide binary is perturbed by fly-by's of field stars and its orbit random-walks and changes over time. This diffusion process can drive the binary into a sufficiently eccentric orbit such that the binary components tidally interact at peri-center and the binary evolves to become a short period binary, which eventually evolves into an LMXB. The formation rate of LMXBs through this channel mostly depends on the number of such BH wide binaries progenitors, which in turn depends on the velocity kicks imparted to BHs (or NSs) at birth. We consider several models for the formation and survival of such wide binaries, and calculate the LMXB formation rates for each model. We find that models where BHs form through direct collapse with no/little natal kicks can give rise to high formation rates comparable with those inferred from observations. This formation scenario had several observational signatures: (1) the number density of LMXBs generally follows the background stellar density and (2) the mass function of the BH stellar companion should be comparable to the mass function of the background stellar population, likely peaking at $0.4-0.6$ M$_{\odot}$. The latter aspect, in particular, is unique to this model compared with previously suggested LMXB formation models following common envelope binary stellar evolution. We note that NS LMXBs can similarly form from wide binaries, but their formation rate through this channel is likely significantly smaller due to their much higher natal kicks.

preprint2014arXiv

Age and mass segregation of multiple stellar populations in galactic nuclei and their observational signatures

Nuclear stellar cluster (NSCs) are known to exist around massive black holes (MBHs) in galactic nuclei. They are thought to have formed through in-situ star formation following gas inflow to the nucleus of the galaxy and/or through the infall of multiple stellar clusters. Here we study the latter, and explore the composite structure of the NSC, and its relation to the various stellar populations originating from its progenitor infalling clusters. We use N-body simulations of clusters infall, and show that this scenario may produce observational signatures in the form of age segregation: the distribution of the stellar properties (e.g. stellar age and/or metallicity) in the NSCs reflect the infall history of the different clusters. The stellar populations of clusters infalling at different times (dynamical ages), are differentially segregated in the NSC, and are not fully mixed even after few Gyrs of evolution. Moreover, the radial properties of stellar populations in the progenitor cluster are mapped to their radial distribution in the final NSC, potentially leading to efficient mass segregation in NSCs, even those where relaxation times are longer than a Hubble time. Finally, the overall structures of the stellar populations present non-spherical configurations and show significant cluster to cluster population differences.

preprint2014arXiv

Effects of Intermediate Mass Black Holes on Nuclear Star Clusters

Nuclear star clusters (NSCs) are dense stellar clusters observed in galactic nuclei, typically hosting a central massive black hole. Here we study the possible formation and evolution of NSCs through the inspiral of multiple star clusters hosting intermediate mass black holes (IMBHs). Using an N-body code we examine the dynamics of the IMBHs and their effects on the NSC. We find that IMBHs inspiral to the core of the newly formed NSC and segregate there. Although the IMBHs scatter each other and the stars, none of them is ejected from the NSC. The IMBHs are excited to high eccentricities and their radial density profile develops a steep power-law cusp. The stars also develop a power-law cusp (instead of the central core that forms in their absence), but with a shallower slope. The relaxation rate of the NSC is accelerated due to the presence of IMBHs, which act as massive-perturbers. This in turn fills the loss-cone and boosts the tidal disruption rate of stars both by the MBH and the IMBHs to a value excluded by rate estimates based on current observations. Rate estimates of tidal disruptions can therefore provide a cumulative constraint on the existence of IMBHs in NSCs.

preprint2014arXiv

Formation and evolution of nuclear star clusters with in-situ star-formation: Nuclear cores and age segregation

Nuclear stellar cluster (NSCs) are known to exist around massive black holes (MBHs) in galactic nuclei. Two formation scenarios were suggested for their origin: (1) Build-up of NSCs from consecutive infall of stellar cluster and (2) Continuous in-situ star-formation. Though the cluster-infall scenario has been extensively studied in recent years, the in-situ formation scenario have been hardly explored. Here we use Fokker-Planck (FP) calculations to study the effects of star formation on the build-up of NSCs and its implications for their long term evolution and their resulting structure. We use the FP equation to describe the evolution of several stellar populations, and add appropriate source terms to account for the effects of newly formed stars. We show that continuous star-formation even 1-2 pc away from the MBH can lead to the build-up of an NSC with properties similar to those of the Milky-way NSC. We also find that the general structure of the old stellar population in the NSC with in-situ star-formation could be very similar to the steady-state Bahcall-Wolf cuspy structure. However, its younger stellar population do not yet achieve a steady state. In particular, formed/evolved NSCs with in-situ star-formation contain differential age-segregated stellar populations which are not yet fully mixed. Younger stellar populations formed in the outer regions of the NSC have a cuspy structure towards the NSC outskirts, while showing a core-like distribution inwards; with younger populations having larger core sizes. In principal, such a structure can give rise to an apparent core-like radial distribution of younger (up to 2-3 Gyrs) stars, as observed in the Galactic center. Such an NSC still preserves an underlying stellar cusp of older stars, that can be potentially be missed by current observations of red-giants.

preprint2014arXiv

Formation of irregular and runaway moons/exomoons through moon-moon scattering

Gas giant planets in the Solar system host large satellite systems with multiple regular and irregular moons. Regular moons revolve around their host planet in circular, low inclination short period orbits, and are thought to form in-situ through coagulation processes. In contrast, irregular moons have highly inclined (and even retrograde), typically eccentric and long period orbits around their host planet. Irregular moons are therefore often thought to be unbound objects in helio-centric orbits that were later captured to their current orbits around the planet. Here we study the a different possibility in which regular moons form in-situ, and dynamically evolve through mutual moon-moon scattering. We find that such evolution can excite the satellites into high eccentricities and inclinations. We find that moons are either ejected from the host planet to become runaway moons, or stay bound and become prograde orbiting irregular moons with inclined and eccentric orbits around their host planet. Ejected moons, unbound to the planet, can later be temporarily re-captured by the host planet even at retrograde orbits. Such moons are eventually re-ejected from the system or collide with the planet, at least in the absence of dissipative processes (e.g. collisions with existing bound moons, a debris disk or through tidal interactions with the host planet), not currently modeled. Uncaptured runaway moons may eventually be ejected from the Solar system, or be captured into stable helio-centric orbits and contribute to the populations of asteroidal or trans-Neptunian objects. Such scenarios are potentially relevant both for the gas-giant satellites in the Solar system and for the dynamical evolution of exomoons.

preprint2014arXiv

On the origin of the B-stars in the Galactic center

We present a new directly-observable statistic which uses sky position and proper motion of stars near the Galactic center massive black hole to identify populations with high orbital eccentricities. It is most useful for stars with large orbital periods for which dynamical accelerations are difficult to determine. We apply this statistic to a data set of B-stars with projected radii 0."1 < p < 25" (~0.004 - 1 pc) from the massive black hole in the Galactic center. We compare the results with those from N-body simulations to distinguish between scenarios for their formation. We find that the scenarios favored by the data correlate strongly with particular K-magnitude intervals, corresponding to different zero-age main-sequence (MS) masses and lifetimes. Stars with 14 < mK < 15 (15 - 20 solar masses, t_{MS} = 8-13 Myr) match well to a disk formation origin, while those with mK > 15 (<15 solar masses, t_{MS} >13 Myr), if isotropically distributed, form a population that is more eccentric than thermal, which suggests a Hills binary-disruption origin.

preprint2014arXiv

On the Origin of Young Stars at the Galactic Center

The center of our galaxy is home to a massive black hole, SgrA*, and a nuclear star cluster containing stellar populations of various ages. While the late type stars may be too old to have retained memory of their initial orbital configuration, and hence formation mechanism, the kinematics of the early type stars should reflect their original distribution. In this contribution we present a new statistic which uses directly-observable kinematical stellar data to infer orbital parameters for stellar populations, and is capable of distinguishing between different origin scenarios. We use it on a population of B-stars in the Galactic center that extends out to large radii (0.5 pc) from the massive black hole. We find that the high K-magnitude population form an eccentric distribution, suggestive of a Hills binary-disruption origin.

preprint2014arXiv

Origin of the Galactic 511 keV emission from positrons produced in irregular supernovae

Gamma ray emission of 511 keV lines arising from electron-positron annihilation has been detected from the Galaxy since the 70s. Spatially resolved observations using the INTEGRAL satellite have shown its full sky distribution to be strongly concentrated in the Galactic bulge, with a smaller contribution from the disk, unlike the situation in any other wavelength. The puzzling distribution of the positrons gave rise to various suggestions, including stellar nucleosynthesis in core-collapse (CC) and type Ia thermonuclear supernovae (SNe), accreting compact objects, and more "exotic" explanations of annihilation of dark-matter particles. However, such models encounter difficulties in reproducing the total Galactic 511 keV flux as well as its peculiar bulge-centered distribution. Theoretical models of SNe from thermonuclear Helium detonations on white dwarfs (WDs) were suggested as potential additional sources of positrons, contributing to the Galactic Gamma-ray emission. Here, we show that the recently discovered class of faint, calcium-rich type Ib SNe (with the prototype SN 2005E), thought to arise from such explosions, also called ".Ia" SNe, can explain both the 511 keV flux and its distribution, and can eliminate the need for non-astrophysical (e.g., dark matter annihilation) sources. Such SNe comprise a fraction of only ~2% of all SNe, but they currently inject hundreds of times more positrons (from $^{44}$Ti decay) to the ISM than injected by CC SNe, enough to reproduce the total Galactic 511 keV line emission. They exclusively explode in old environments, and their contribution to the 511 keV emission therefore follows the old (>10 Gyrs) Galactic stellar population, dominated by the bulge, thereby reproducing the observed large bulge-to-disk ratio.

preprint2014arXiv

Secular dynamics in hierarchical three-body systems with mass loss and mass transfer

Recent studies have shown that secular evolution of triple systems can play a major role in the evolution and interaction of their inner binaries. Very few studies explored the stellar evolution of triple systems, and in particular the mass loss phase of the evolving stellar components. Here we study the dynamical secular evolution of hierarchical triple systems undergoing mass loss. We use the secular evolution equations and include the effects of mass-loss and mass-transfer, as well as general relativistic effects. We present various evolutionary channels taking place in such evolving triples, and discuss both the effects of mass-loss and mass-transfer in the inner binary system, as well as the effects of mass-loss/transfer from an outer third companion. We discuss several distinct types/regimes of triple secular evolution, where the specific behavior of a triple system can sensitively depend on its hierarchy and the relative importance of classical and general relativistic effects. We show that the orbital changes due to mass-loss and/or mass-transfer processes can effectively transfer a triple system from one dynamical regime to another. In particular, mass loss/transfer can both induce and quench high amplitude (Lidov-Kozai) variations in the eccentricity and inclination of the inner binaries of evolving triples. They can also change the system dynamics from an orderly periodic behavior to a chaotic one, and vice versa.

preprint2014arXiv

Secular Evolution Of Binaries Near Massive Black Holes: Formation of compact binaries, merger/collision products and G2-like objects

Here we discuss the evolution of binaries around MBH in nuclear stellar clusters. We focus on their secular evolution due to the perturbation by the MBH, while simplistically accounting for their collisional evolution. Binaries with highly inclined orbits in respect to their orbit around the MBH are strongly affected by secular processes, which periodically change their eccentricities and inclinations (e.g. Kozai--Lidov cycles). During periapsis approach, dissipative processes such as tidal friction may become highly efficient, and may lead to shrinkage of a binary orbit and even to its merger. Binaries in this environment can therefore significantly change their orbital evolution due to the MBH third-body perturbative effects. Such orbital evolution may impinge on their later stellar evolution. Here we follow the secular dynamics of such binaries and its coupling to tidal evolution, as well as the stellar evolution of such binaries on longer time-scales. We find that stellar binaries in the central parts of NSCs are highly likely to evolve into eccentric and/or short period binaries, and become strongly interacting binaries either on the main sequence (at which point they may even merge), or through their later binary stellar evolution. The central parts of NSCs therefore catalyze the formation and evolution of strongly interacting binaries, and lead to the enhanced formation of blue stragglers, X-ray binaries, gravitational wave sources and possible supernova progenitors. Induced mergers/collisions may also lead to the formation of G2-like cloud-like objects such as the one recently observed in the Galactic center.

preprint2014arXiv

Stellar dynamics in gas: The role of gas damping

In this paper, we consider how gas damping affects the dynamical evolution of gas-embedded star clusters. Using a simple three-component (i.e. one gas and two stellar components) model, we compare the rates of mass segregation due to two-body relaxation, accretion from the interstellar medium, and gas dynamical friction in both the supersonic and subsonic regimes. Using observational data in the literature, we apply our analytic predictions to two different astrophysical environments, namely galactic nuclei and young open star clusters. Our analytic results are then tested using numerical simulations performed with the NBSymple code, modified by an additional deceleration term to model the damping effects of the gas. The results of our simulations are in reasonable agreement with our analytic predictions, and demonstrate that gas damping can significantly accelerate the rate of mass segregation. A stable state of approximate energy equilibrium cannot be achieved in our model if gas damping is present, even if Spitzer's Criterion is satisfied. This instability drives the continued dynamical decoupling and subsequent ejection (and/or collisions) of the more massive population. Unlike two-body relaxation, gas damping causes overall cluster contraction, reducing both the core and half-mass radii. If the cluster is mass segregated (and/or the gas density is highest at the cluster centre), the latter contracts faster than the former, accelerating the rate of core collapse.

preprint2014arXiv

The Multiple Origin of Blue Straggler Stars: Theory vs. Observations

In this chapter we review the various suggested channels for the formation and evolution of blue straggler stars (BSSs) in different environments and their observational predictions. These include mass transfer during binary stellar evolution - case A/B/C and D (wind Roche-lobe overflow) mass transfer, stellar collisions during single and binary encounters in dense stellar cluster, and coupled dynamical and stellar evolution of triple systems. We also explore the importance of the BSS and binary dynamics in stellar clusters. We review the various observed properties of BSSs in different environments (halo and bulge BSSs, BSSs in globular clusters and BSSs in old open clusters), and compare the current observations with the theoretical predictions for BSS formation. We try to constrain the likely progenitors and processes that play a role in the formation of BSSs and their evolution. We find that multiple channels of BSS formation are likely to take part in producing the observed BSSs, and we point out the strengths and weaknesses of each the formation channel in respect to the observational constraints. Finally we point out directions to further explore the origin of BSS, and highlight eclipsing binary BSSs as important observational tool.

preprint2013arXiv

Environment-derived constraints on the progenitors of low-luminosity type I supernovae

We present a study of the properties of the host galaxies of unusual transient objects of two types, both being sub-luminous compared with the major classes of supernovae. Those of one type exhibit unusually strong calcium features, and have been termed 'Ca-rich'. Those of the second type, with SN2002cx as the prototype and SN2008ha as the most extreme example to date, have some properties in common with the first, but show typically lower ejecta velocities, and different early spectra. We confirm important differences in the environments of the two types, with the Ca-rich transients preferentially occurring in galaxies dominated by old stellar populations. Quantitatively, the association of the the Ca-rich transients with regions of ongoing star formation is well matched to that of type Ia supernovae. The SN2002cx-like transients are very different, with none of the present sample occurring in an early-type host, and a statistical association with star-formation regions similar to that of type II-P supernovae, and therefore a delay time of 30-50 Myrs.

preprint2013arXiv

Evolution of second generation stars in stellar disks of globular and nuclear clusters: Omega Centauri as a test case

Globular clusters (GCs) and many nuclear clusters (NCs) show evidence for hosting multiple generations of stellar populations. Younger stellar populations in NCs appear to reside in disk like structures, including the nuclear cluster in our own Galactic center as well as in M31. Kinematic studies of the anomalous globular cluster Omega Centauri, thought to possibly be a former dwarf galaxy (or a galactic nucleus), show evidence for its hosting of a central, kinematically cold disk component. These observations suggest that formation of second (or multiple) generation stars may occur in flattened disk like structures. Here we use detailed N-body simulations to explore the possible evolution of such stellar disks, embedded in globular clusters. We follow the long term evolution of a disk like structure similar to that observed in Omega Centauri and study its properties. We find that a stellar-disk like origin for second generation stellar populations can leave behind significant kinematic signatures in properties of the clusters, including an anisotropic distribution, and lower velocity dispersions, which can be used to constrain the origin of second generations stars and their dynamical evolution.

preprint2013arXiv

Stability of Satellites in Closely Packed Planetary Systems

We perform numerical integrations of four-body (star, planet, planet, satellite) systems to investigate the stability of satellites in planetary Systems with Tightly-packed Inner Planets (STIPs). We find that the majority of closely-spaced stable two-planet systems can stably support satellites across a range of parameter-space which is only slightly decreased compared to that seen for the single-planet case. In particular, circular prograde satellites remain stable out to $\sim 0.4 R_H$ (where $R_H$ is the Hill Radius) as opposed to $\sim 0.5 R_H$ in the single-planet case. A similarly small restriction in the stable parameter-space for retrograde satellites is observed, where planetary close approaches in the range 2.5 to 4.5 mutual Hill radii destabilize most satellites orbits only if $a\sim 0.65 R_H$. In very close planetary pairs (e.g. the 12:11 resonance) the addition of a satellite frequently destabilizes the entire system, causing extreme close-approaches and the loss of satellites over a range of circumplanetary semi-major axes. The majority of systems investigated stably harbored satellites over a wide parameter-space, suggesting that STIPs can generally offer a dynamically stable home for satellites, albeit with a slightly smaller stable parameter-space than the single-planet case. As we demonstrate that multi-planet systems are not a priori poor candidates for hosting satellites, future measurements of satellite occurrence rates in multi-planet systems versus single-planet systems could be used to constrain either satellite formation or past periods of strong dynamical interaction between planets.

preprint2012arXiv

Eccentric disc instability in stellar discs formed from inspiraling gas clouds in the Galactic Centre

The inspiral of a turbulent molecular cloud in the Galactic Centre may result in the formation of a small, dense and moderately eccentric gas disc around the supermassive black hole (SMBH). Such a disc is unstable to fragmentation and may lead to the formation of young massive stars in the central parsec of the Galaxy. Here we perform high-accuracy direct summation N-body simulations of a ring of massive stars (with initial semi-major axes 0.1 < a/pc < 0.4 and eccentricities 0.2 < e < 0.4), subject to the potential of the SMBH, a stellar cusp, and the parent gas disc, to study how the orbital elements of the ring evolve in time. The initial conditions for the stellar ring are drawn from the results of previous simulations of molecular cloud infall and disruption in the SMBH potential. While semi-major axes do not evolve significantly, the distribution of eccentricities spreads out very fast (~1 Myr) as a consequence of cusp precession. In particular, stellar orbits with initial eccentricity e>0.3 (e<0.3) tend to become even more (less) eccentric, resulting in a bimodal eccentricity distribution. The distribution is qualitatively consistent with that of the massive stars observed in the Galactic Centre's clockwise disc.

preprint2012arXiv

Failed-Detonation Supernovae: Sub-Luminous Low-Velocity Ia Supernovae and Their Kicked Remnant White Dwarfs with Iron-Rich Cores

Type Ia supernovae (SNe Ia) originate from the thermonuclear explosions of carbon-oxygen (C-O) white dwarfs (WDs). The single-degenerate scenario is a well-explored model of SNe Ia where unstable thermonuclear burning initiates in an accreting, Chandrasekhar-mass WD and forms an advancing flame. By several proposed physical processes the rising, burning material triggers a detonation, which subsequently consumes and unbinds the WD. However, if a detonation is not triggered and the deflagration is too weak to unbind the star, a completely different scenario unfolds. We explore the failure of the Gravitationally-Confined Detonation (GCD) mechanism of SNe Ia, and demonstrate through 2D and 3D simulations the properties of failed-detonation SNe. We show that failed-detonation SNe expel a few 0.1 solar masses of burned and partially-burned material and that a fraction of the material falls back onto the WD, polluting the remnant WD with intermediate-mass and iron-group elements, that likely segregate to the core forming an WD whose core is iron rich. The remaining material is asymmetrically ejected at velocities comparable to the escape velocity from the WD, and in response, the WD is kicked to velocities of a few hundred km/s. These kicks may unbind the binary and eject a runaway/hyper-velocity WD. Although the energy and ejected mass of the failed-detonation SN are a fraction of typical thermonuclear SNe, they are likely to appear as sub-luminous low-velocity SNe Ia. Such failed detonations might therefore explain or are related to the observed branch of peculiar SNe Ia, such as the family of low-velocity sub-luminous SNe (SN 2002cx/SN 2008ha-like SNe).

preprint2012arXiv

On the origin of planets at very wide orbits from the re-capture of free floating planets

In recent years several planets have been discovered at wide orbits (>100 AU) around their host stars. Theoretical studies encounter difficulties in explaining their formation and origin. Here we propose a novel scenario for the production of planetary systems at such orbits, through the dynamical recapture of free floating planets (FFPs) in dispersing stellar clusters. This process is a natural extension of the recently suggested scenario for the formation of wide stellar binaries. We use N-body simulations of dispersing clusters with 10-1000 stars and comparable numbers of FFPs to study this process. We find that planets are captured into wide orbits in the typical range ~100-10^6 AU, and have a wide range of eccentricities (thermal distribution). Typically, 3-6 x (f_FFP/1) % of all stars capture a planetary companion with such properties (where f_FFP is the number of FFP per star). The planetary capture efficiency is comparable to that of capture-formed stellar-binaries, and shows a similar dependence on the cluster size and structure. It is almost independent of the specific planetary mass; planets as well as sub-stellar companions of any mass can be captured. The capture efficiency decreases with increasing cluster size, and for a given cluster size the it increases with the host/primary mass. More than one planet can be captured around the same host and planets can be captured into binary systems. Planets can also be captured into pre-existing planetary and into orbits around black holes and massive white dwarfs, if these formed early enough before the cluster dispersal. In particular, stellar black holes have a high capture efficiency (>50 % and 5-10 x (f_FFP/1) % for capture of stars and planetary companions, respectively) due to their large mass. Finally, although rare, two FFPs or brown dwarfs can become bound and form a FFP-binary system with no stellar host.

preprint2012arXiv

Star Hoppers: Planet Instability and Capture in Evolving Binary Systems

Many planets are observed in stellar binary systems, and their frequency may be comparable to that of planetary systems around single stars. Binary stellar evolution in such systems influences the dynamical evolution of the resident planets. Here we study the evolution of a single planet orbiting one star in an evolving binary system. We find that stellar evolution can trigger dynamical instabilities that drive planets into chaotic orbits. This instability leads to planet-star collisions, exchange of the planet between the binary stars ("star-hoppers"), and ejection of the planet from the system. The means by which planets can be recaptured is similar to the pull-down capture mechanism for irregular solar system satellites. Because planets often suffer close encounters with the primary on the asymptotic giant branch, captures during a collision with the stellar envelope are also possible. Such capture could populate the habitable zone around white dwarfs.

preprint2012arXiv

The properties of dynamically ejected runaway and hyper-runaway stars

Runaway stars are stars observed to have large peculiar velocities. Two mechanisms are thought to contribute to the ejection of runaway stars, both involve binarity (or higher multiplicity). In the binary supernova scenario a runaway star receives its velocity when its binary massive companion explodes as a supernova (SN). In the alternative dynamical ejection scenario, runaway stars are formed through gravitational interactions between stars and binaries in dense, compact clusters or cluster cores. Here we study the ejection scenario. We make use of extensive N-body simulations of massive clusters, as well as analytic arguments, in order to to characterize the expected ejection velocity distribution of runaways stars. We find the ejection velocity distribution of the fastest runaways (>~80 km s^-1) depends on the binary distribution in the cluster, consistent with our analytic toy model, whereas the distribution of lower velocity runaways appears independent of the binaries properties. For a realistic log constant distribution of binary separations, we find the velocity distribution to follow a simple power law; Gamma(v) goes like v^(-8/3) for the high velocity runaways and v^(-3/2) for the low velocity ones. We calculate the total expected ejection rates of runaway stars from our simulated massive clusters and explore their mass function and their binarity. The mass function of runaway stars is biased towards high masses, and depends strongly on their velocity. The binarity of runaways is a decreasing function of their ejection velocity, with no binaries expected to be ejected with v>150 km s^-1. We also find that hyper-runaways with velocities of hundreds of km s^-1 can be dynamically ejected from stellar clusters, but only at very low rates, which cannot account for a significant fraction of the observed population of hyper-velocity stars in the Galactic halo.

preprint2012arXiv

The triple evolution dynamical instability: Stellar collisions in the field and the formation of exotic binaries

Physical collisions and close approaches between stars play an important role in the formation of exotic stellar systems. Standard theories suggest that collisions are rare, occurring only via random encounters between stars in dense clusters. We present a different formation pathway, the triple evolution dynamical instability (TEDI), in which mass loss in an evolving triple star system causes orbital instability. The subsequent chaotic orbital evolution of the stars triggers close encounters, collisions, exchanges between the stellar components, and the dynamical formation of eccentric compact binaries (including Sirius like binaries). We demonstrate that the rate of stellar collisions due to the TEDI is approximately 10^{-4} yr^{-1} per Milky-Way Galaxy, which is nearly 30 times higher than the total collision rate due to random encounters in the Galactic globular clusters. Moreover, we find that the dominant type of stellar collisions is qualitatively different; most collisions involve asymptotic giant branch stars, rather than main sequence, or slightly evolved stars, which dominate collisions in globular clusters. The TEDI mechanism should lead us to revise our understanding of collisions and the formation of compact, eccentric binaries in the field.

preprint2012arXiv

Wind-accretion disks in wide binaries, second generation protoplanetary disks and accretion onto white dwarfs

Mass transfer from an evolved donor star to its binary companion is a standard feature of stellar evolution in binaries. In wide binaries, the companion star captures some of the mass ejected in a wind by the primary star. The captured material forms an accretion disk. Here, we study the evolution of wind-accretion disks, using a numerical approach which allows us to follow the long term evolution. For a broad range of initial conditions, we derive the radial density and temperature profiles of the disk. In most cases, wind-accretion leads to long-lived stable disks over the lifetime of the AGB donor star. The disks have masses of a few times 10^{-5}-10^{-3} M_sun, with surface density and temperature profiles that follow broken power-laws. The total mass in the disk scales approximately linearly with the viscosity parameter used. Roughly 50% to 80% of the mass falling into the disk accretes onto the central star; the rest flows out through the outer edge of the disk into the stellar wind of the primary. For systems with large accretion rates, the secondary accretes as much as 0.1 M_sun. When the secondary is a white dwarf, accretion naturally leads to nova and supernova eruptions. For all types of secondary star, the surface density and temperature profiles of massive disks resemble structures observed in protoplanetary disks, suggesting that coordinated observational programs might improve our understanding of uncertain disk physics.

preprint2011arXiv

An Emerging Class of Bright, Fast-evolving Supernovae with Low-mass Ejecta

A recent analysis of supernova (SN) 2002bj revealed that it was an apparently unique type Ib SN. It showed a high peak luminosity, with absolute magnitude M_R -18.5, but an extremely fast-evolving light curve. It had a rise time of <7 days followed by a decline of 0.25 mag per day in B-band, and showed evidence for very low mass of ejecta (<0.15 M_Sun). Here we discuss two additional historical events, SN 1885A and SN 1939B, showing similarly fast light curves and low ejected masses. We discuss the low mass of ejecta inferred from our analysis of the SN 1885A remnant in M31, and present for the first time the spectrum of SN 1939B. The old environments of both SN 1885A (in the bulge of M31) and SN 1939B (in an elliptical galaxy with no traces of star formation activity), strongly support old white dwarf progenitors for these SNe. We find no clear evidence for helium in the spectrum of SN 1939B, as might be expected from a helium-shell detonation on a white dwarf, suggested to be the origin of SN 2002bj. Finally, the discovery of all the observed fast-evolving SNe in nearby galaxies suggests that the rate of these peculiar SNe is at least 1-2 % of all SNe.

preprint2011arXiv

Calcium-rich gap transients in the remote outskirts of galaxies

From the first two seasons of the Palomar Transient Factory, we identify three peculiar transients (PTF09dav, PTF10iuv, PTF11bij) with five distinguishing characteristics: peak luminosity in the gap between novae and supernovae (M_R = 15.5 to -16.5), rapid photometric evolution (rise-time ~12--15 days), large photospheric velocities (~6000 to 11000 km/s), early spectroscopic evolution into nebular phase (~1 to 3 months) and peculiar nebular spectra dominated by Calcium. We also culled the extensive decade-long Lick Observatory Supernova Search database and identified an additional member of this group, SN 2007ke. Our choice of photometric and spectroscopic properties was motivated by SN 2005E (Perets et al. 2010). To our surprise, as in the case of SN 2005E, all four members of this group are also clearly offset from the bulk of their host galaxy. Given the well-sampled early and late-time light curves, we derive ejecta masses in the range of 0.4--0.7 Msun. Spectroscopically, we find that there may be a diversity in the photospheric phase, but the commonality is in the unusual nebular spectra. Our extensive follow-up observations rule out standard thermonuclear and standard core-collapse explosions for this class of "Calcium-rich gap" transients. If the progenitor is a white dwarf, we are likely seeing a detonation of the white dwarf core and perhaps, even shockfront interaction with a previously ejected nova shell. In the less likely scenario of a massive star progenitor, a very non-standard channel specific to a low-metallicity environment needs to be invoked (e.g., ejecta fallback leading to black hole formation). Detection (or lack thereof) of a faint underlying host (dwarf galaxy, cluster) will provide a crucial and decisive diagnostic to choose between these alternatives.

preprint2011arXiv

Planets in evolved binary systems

Exoplanets are typically thought to form in protoplanetary disks left over from protostellar disk of their newly formed host star. However, additional planetary formation and evolution routes may exist in old evolved binary systems. Here we discuss the implications of binary stellar evolution on planetary systems. In these binary systems stellar evolution could lead to the formation of symbiotic stars, where mass is lost from one star and could be transferred to its binary companion, and may form an accretion disk around it. This raises the possibility that such a disk could provide the necessary environment for the formation of a new, second generation of planets in both circumstellar or circumbinary configurations. Pre-existing first generation planets surviving the post-MS evolution of such systems would be dynamically effected by the mass loss in the systems and may also interact with the newly formed disk. Second generation planetary systems should be typically found in white dwarf binary systems, and may show various observational signatures. Most notably, second generation planets could form in environment which are inaccessible, or less favorable, for first generation planets. The orbital phase space available for the second generation planets could be forbidden (in terms of the system stability) to first generation planets in the pre-evolved progenitor binaries. Observations of exo-planets in such forbidden or unfavorable regions could possibly serve to uniquely identify their second generation character. Finally, we point out a few observed candidate second generation planetary systems, including Gl 86, HD 27442 and all of the currently observed circumbinary planet candidates. A second generation origin for these systems could naturally explain their unique configurations.

preprint2011arXiv

The old environment of the faint calcium-rich supernova SN 2005cz

The supernova SN 2005cz has recently attracted some attention, due to the fact that it was spectroscopically similar to type Ib supernovae (SNe), a class that is presumed to result from core-collapse of massive stars, yet it occurred in an elliptical galaxy, where one expects very few massive stars to exist. Two explanations for this remarkable event were put forward. Perets et al. (2010) associate SN 2005cz with the class of Ca-rich, faint SNe Ib, which likely result from old double-white-dwarf systems with a He-rich secondary. On the other hand, Kawabata et al. (2010) suggest that SN 2005cz is indeed a core-collapse event (in a binary system), albeit of a star at the lower end of the mass range, 10-12 M_Sun. The existence of this star in its elliptical host is explained as resulting from low-level star formation (SF) activity in that galaxy. Here we present extensive observations of the location of SN 2005cz, sensitive to a variety of SF tracers, including optical spectroscopy, H_alpha emission, UV emission and HST photometry. We show that NGC 4589, the host galaxy of SN 2005cz, does not show any signatures of a young stellar population or recent SF activity either close to or far from the location of SN 2005cz.

preprint2011arXiv

Triple Disruptions in The Galactic Centre: Captured and Ejected Binaries, Rejuvenated Stars, and Correlated Orbits

The disruption of a binary star by a massive black hole (MBH) typically leads to the capture of one component around the MBH and the ejection of its companion at a high velocity, possibly producing a hypervelocity star. The high fraction of observed triples ($\sim10$% for F/G/K stars and $\sim50$% for OB stars) give rise to the possibility of the disruption of triples by a MBH. Here we study this scenario, and use direct $N$-body integrations to follow the orbits of thousands of triples, during and following their disruption by a MBH (of $4\times10^6$ M$_\odot$, similar to the MBH existing in the Galactic Centre; SgrA$^*$). We find that triple disruption can lead to several outcomes and we discuss their relative frequency. Beside the ejection/capture of single stars, similar to the binary disruption case, the outcomes of triple disruption include the ejection of hypervelocity binaries; capture of binaries around the MBH; collisions between two or all of the triple components (with low enough velocities that could lead to their merger); and the capture of two or even three stars at close orbits around the MBH. The orbits of single stars captured in a single disruption event are found to be correlated. The eccentricity of the mutual orbits of captured/ejected binaries is typically excited to higher values. Stellar evolution of captured/ejected binaries may later result in their coalescence/strong interaction and the formation of hypervelocity blue stragglers or merger remnants in orbits around SgrA*. Finally, the capture of binaries close to the MBH can replenish and increase the binary frequency near the MBH, which is otherwise very low.

preprint2011arXiv

Wind-shearing in gaseous protoplanetary disks and the evolution of binary planetesimals

One of the first stages of planet formation is the growth of small planetesimals. This early stage occurs much before the dispersal of most of the gas from the protoplanetary disk. Due to their different aerodynamic properties, planetesimals of different sizes and shapes experience different drag forces from the gas during this time. Such differential forces produce a wind-shearing (WISH) effect between close by, different size planetesimals. For any two planetesimals, a WISH radius can be considered, at which the differential acceleration due to the wind becomes greater than the mutual gravitational pull between the planetesimals. We find that the WISH radius could be much smaller than the Hill radius, i.e. WISH could play a more important role than tidal perturbations by the star. Here we study the WISH radii for planetesimal pairs of different sizes and compare the effects of wind and gravitational shearing (drag force vs. gravitational tidal force). We then discuss the role of WISH for the stability and survival of binary planetesimals. Binaries are sheared apart by the wind if they are wider than their WISH radius. WISH-stable binaries can inspiral and possibly coalesce due to gas drag. Here, we calculate the WISH radius and the gas drag-induced merger timescale, providing stability and survival criteria for gas-embedded binary planetesimals. Our results suggest that even WISH-stable binaries may merge in times shorter than the lifetime of the gaseous disk. This may constrain currently observed binary planetesimals to have formed far from the star or at a late stage after the dispersal of most of the disk gas. We note that the WISH radius may also be important for other processes such as planetesimal erosion and planetesimal encounters and collisions in a gaseous environment.

preprint2010arXiv

Binary planetesimals and their role in planet formation

One of the main evolutionary stages of planet formation is the dynamical evolution of planetesimal disks. These disks are thought to evolve through gravitational encounters and physical collisions between single planetesimals. In recent years, many binary planetesimals have been observed in the Solar system, indicating that the binarity of planetesimals is high. However, current studies of planetesimal disks formation and evolution do not account for the role of binaries. Here we point out that gravitational encounters of binary planetesimals can have an important role in the evolution of planetesimal disks. Binary planetesimals catalyze close encounters between planetesimals, and can strongly enhance their collision rate. Binaries may also serve as additional heating source of the planetesimal disk, through the exchange of the binaries gravitational potential energy into the kinetic energy of planetesimals in the disk.

preprint2010arXiv

Dynamical constraints on the origin of the young B-stars in the Galactic center

Regular star formation is thought to be inhibited close to the massive black hole (MBH) in the Galactic center. Nevertheless, tens of young main sequence B stars have been observed in an isotropic distribution close to it. Various models have been suggested for the formation of the B-stars closest to the MBH (<0.05 pc; the S-stars), typically involving the migration of these stars from their original birthplace to their currently observed position. Here we explore the orbital phase space distribution of the B-stars throughout the central pc expected from the various suggested models for the origin of the B-stars. We find that most of these models have difficulties in explaining, by themselves, both the population of the S-stars (<0.05 pc), and the population of the young B-stars further away (up to 0.5 pc). Most models grossly over-predict the number of B-stars up to 0.5 pc, given the observed number of S-stars. Such models include the intermediate-mass black hole assisted cluster inspiral scenario, Kozai-like perturbations by two disks, spiral density waves migration in a gaseous disk, and some of the eccentric disk instability models. We focus on one of the other models, the massive perturber induced binary disruption, which is consistent with both the S-stars and the extended population of B-stars further away. For this model we use analytical arguments and N-body simulations to provide further observational predictions. These could be compared with future observations to further support this model, constrain it or refute it. These predictions include the radial distribution of the young B-stars, their eccentricity distribution and its dependence on distance from the MBH (higher eccentricities at larger distances from the MBH), as well as less specific expectations regarding their mass function.

preprint2010arXiv

Dynamics and origins of the young stars in the Galactic center

The environment near the massive black hole (MBH) in the Galactic center is very hostile for star formation. Nevertheless, many young stars (both O and B stars) are observed close the MBH. The B-stars seems to have an isotropic, continuous distribution between 0.01 pc and up to a pc. The O (and Wolf-Rayet;WR) stars, in contrast, seem to be distributed in a coherent disk like configuration, extending only between ~0.04 pc to ~0.5 pc. Our current understanding favors an in-situ formation origin for the more massive (O and WR) stars, in gaseous disk and/or streams from an in-falling gas clump. The B-stars seem to have a different origin, more likely through a dynamical capture, following a binary disruption by the MBH. This scenario could also be able to explain the origin of hypervelocity stars in the Galactic halo. These and other possible origins of the young stars in the Galactic center are briefly reviewed and their possible observational signatures and constraints are detailed.

preprint2010arXiv

Interaction of Atomic and Molecular Hydrogen with Tholin Surfaces at Low Temperatures

We study the interaction of atomic and molecular hydrogen with a surface of tholin, a man-made polymer considered to be an analogue of aerosol particles present in Titan's atmosphere, using thermal programmed desorption at low temperatures below 30 K. The results are fitted and analyzed using a fine-grained rate equation model that describes the diffusion, reaction and desorption processes. We obtain the energy barriers for diffusion and desorption of atomic and molecular hydrogen. These barriers are found to be in the range of 30 to 60 meV, indicating that atom/molecule-surface interactions in this temperature range are dominated by weak adsorption forces. The implications of these results for the understanding of the atmospheric chemistry of Titan are discussed.

preprint2010arXiv

Orbital properties of binary minor planets

Many binary minor planets (BMPs; both binary asteroids and binary Trans-Neptunians objects; TNOs) are known to exist in the Solar system. The currently observed orbital and physical properties of BMPs hold essential information and clues about their origin, their evolution and the conditions underwhich they evolved. Here we study the orbital properties of BMPs with currently known orbital solutions (the BMPs mutual orbits and not the BMP orbits around the sun). We find that BMPs are typically highly inclined relative to their orbit around the sun, with a distribution consistent with an isotropic distribution. BMPs not affected by tidal forces are found to have high eccentricities with non-thermal eccentricity distribution peaking at intermediate eccentricities (typically 0.4-0.6). The high inclinations and eccentricities of the BMPs suggest that BMPs evolved in a dense collisional environment, in which gravitational encounters in addition to tidal and secular Kozai affects played an important role in their orbital evolution.

preprint2010arXiv

Second generation planets

Planets are typically thought to form in protoplanetary disks left over from protostellar disk of their newly formed host star. However, an additional planetary formation route may exist in old evolved binary systems. In such systems stellar evolution could lead to the formation of symbiotic stars, where mass transferred from the expanding evolved star to its binary companion could form an accretion disk around it. Such a disk could provide the necessary environment for the formation of a second generation of planets in both circumstellar or circumbinary configurations. Pre-existing first generation planets surviving the post-MS evolution of such systems may serve as seeds for, and/or interact with, the second generation planets, possibly forming atypical planetary systems. Second generation planetary systems should be typically found in white dwarf binary systems, and may show various observational signatures. Most notably, second generation planets could form in environment which are inaccessible, or less favorable, for first generation planets. The orbital phase space available for the second generation planets could be forbidden (in terms of the system stability) to first generation planets in the pre-evolved progenitor binaries. In addition planets could form in metal poor environments such as globular clusters and/or in double compact object binaries. Observations of planets in such forbidden or unfavorable regions may serve to uniquely identify their second generation character. Finally, we point out a few observed candidate second generation planetary systems, including PSR B1620-26 (in a globular cluster), Gl 86, HD 27442 and all of the currently observed circumbinary planet candidates. A second generation origin for these systems could naturally explain their unique configurations.

preprint2010arXiv

Wind-shearing in gaseous protoplanetary disks

One of the first stages of planet formation is the growth of small planetesimals and their accumulation into large planetesimals and planetary embryos. This early stage occurs much before the dispersal of most of the gas from the protoplanetary disk. Due to their different aerodynamic properties, planetesimals of different sizes/shapes experience different drag forces from the gas at these stage. Such differential forces produce a wind-shearing effect between close by, different size planetesimals. For any two planetesimals, a wind-shearing radius can be considered, at which the differential acceleration due to the wind becomes greater than the mutual gravitational pull between the planetesimals. We find that the wind-shearing radius could be much smaller than the gravitational shearing radius by the Sun (the Hill radius), i.e. during the gas-phase of the disk wind-shearing could play a more important role than tidal perturbations by the Sun. Here we study the wind-shearing radii for planetesimal pairs of different sizes and compare it with gravitational shearing (drag force vs. gravitational tidal forces). We then discuss the role of wind-shearing for the stability and survival of binary planetesimals, and provide stability criteria for binary planetesimals embedded in a gaseous disk.

preprint2009arXiv

Classical Diffusion of a quantum particle in a noisy environment

We study the spreading of a quantum-mechanical wavepacket in a one-dimensional tight-binding model with a noisy potential, and analyze the emergence of classical diffusion from the quantum dynamics due to decoherence. We consider a finite correlation time of the noisy environment, and treat the system by utilizing the separation of fast (dephasing) and slow (diffusion) processes. We show that classical diffusive behavior emerges at long times, and we calculate analytically the dependence of the classical diffusion coefficient on the noise magnitude and correlation time. This method provides a general solution to this problem for arbitrary conditions of the noisy environment. The results are relevant to a large variety of physical systems, from electronic transport in solid state physics, to light transmission in optical devices, diffusion of excitons, and quantum computation.

preprint2009arXiv

Kozai cycles, tidal friction and the dynamical evolution of binary minor planets

In recent years many binary minor planets (BMPs) have been discovered in the Solar system. Many models have been suggested for their formation, but these encounter difficulties explaining their observed characteristics. Here we show that secular perturbations by the Sun (Kozai mechanism) fundamentally change the evolution and the initial distribution of BMPs predicted by such models and lead to unique observational signatures. The Kozai mechanism can lead to a large periodic oscillations in the eccentricity and inclination of highly inclined BMP orbits, where we predict such effects to be observable with current accuracy within a few years (e.g. for the binary asteroid Huenna). In addition, the combined effects of the Kozai mechanism and tidal friction (KCTF) drives BMPs into short period circular orbits. We predict a specific inclination dependent distribution of the separation and eccentricity of BMPs, due to these effects, including a zone of avoidance at the highest inclinations. Specifically the Kozai evolution could explain the recently observed peculiar orbit of the Kuiper belt binary 2001 QW322 . Additionally, the KCTF process could lead to BMPs coalescence and serve as an important route for the formation of irregular shaped single minor planets with large axial tilts.

preprint2009arXiv

On the triple origin of blue stragglers

Blue straggler stars (BSSs) are stars observed to be hotter and bluer than other stars with the same luminosity in their environment. As such they appear to be much younger than the rest of the stellar population. Two main channels have been suggested to produce such stars: (1) collisions between stars in clusters or (2) mass transfer between, or merger of, the components of primordial short-period binaries. Here we suggest a third scenario, in which the progenitor of BSSs are formed in primordial (or dynamically formed) hierarchical triple stars. In such configurations the dynamical evolution of the triples through the Kozai mechanism and tidal friction can induce the formation of very close inner binaries. Angular momentum loss in a magnetized wind or stellar evolution could then lead to the merger of these binaries (or to mass transfer between them) and produce BSSs in binary (or triple) systems. We study this mechanism and its implications and show that it could naturally explain many of the characteristics of the BSS population in clusters, most notably the large binary fraction of long period BSS binaries; their unique period-eccentricity distribution (with typical periods > 700 days); and the typical location of these BSSs in the color-magnitude diagram, far from the cluster turn-off point of their host clusters. We suggest that this scenario has a major (possibly dominant) role in the formation of BSSs in open clusters and give specific predictions for the the BSSs population formed in this manner. We also note that triple systems may be the progenitors of the brightest planetary nebulae in old elliptical galaxies, which possibly evolved from BSSs.

preprint2009arXiv

The Galactic potential and the asymmetric distribution of hypervelocity stars

In recent years several hypervelocity stars (HVSs) have been observed in the halo of our Galaxy. Such HVSs have possibly been ejected from the Galactic center and then propagated in the Galactic potential up to their current position. The recent survey for candidate HVSs show an asymmetry in the kinematics of candidate HVSs (position and velocity vectors), where more outgoing stars than ingoing stars (i.e. positive Galactocentric velocities vs. negative ones) are observed. We show that such kinematic asymmetry, which is likely due to the finite lifetime of the stars and Galactic potential structure, could be used in a novel method to probe and constrain the Galactic potential, identify the stellar type of the stars in the survey and estimate the number of HVSs. Kinematics-independent identification of the stellar types of the stars in such surveys (e.g. spectroscopic identification) could further improve these results. We find that the observed asymmetry between ingoing and outgoing stars favors specific Galactic potential models. It also implies a lower limit of ~54+-8 main sequence HVSs in the survey sample (>=648+-96 in the Galaxy), assuming that all of the main sequence stars in the survey originate from the Galactic center. The other stars in the survey are likely to be hot blue horizontal branch stars born in the halo rather than stars ejected from the Galactic center.

preprint2007arXiv

Molecular Hydrogen Formation on Amorphous Silicates Under Interstellar Conditions

Experimental results on the formation of molecular hydrogen on amorphous silicate surfaces are presented for the first time and analyzed using a rate equation model. The energy barriers for the relevant diffusion and desorption processes are obtained. They turn out to be significantly higher than those obtained earlier for polycrystalline silicates, demonstrating the importance of grain morphology. Using these barriers we evaluate the efficiency of molecular hydrogen formation on amorphous silicate grains under interstellar conditions. It is found that unlike polycrystalline silicates, amorphous silicate grains are efficient catalysts of H$_{2}$ formation within a temperature range which is relevant to diffuse interstellar clouds. The results also indicate that the hydrogen molecules are thermalized with the surface and desorb with low kinetic energy. Thus, they are unlikely to occupy highly excited states.

preprint2006arXiv

Massive perturber-driven interactions of stars with a massive black hole

We study the role of massive perturbers (MPs) in deflecting stars and binaries to almost radial ("loss-cone") orbits, where they pass near the central massive black hole (MBH), interact with it at periapse, and are ultimately destroyed. MPs dominate dynamical relaxation when the ratio of the 2nd moments of the MP and star mass distributions, mu_2=(N_p*<M_p^2>)/(N_s*<M_s^2>), satisfies mu_2>>1. We compile the MP mass function from published observations, and show that MPs in the nucleus of the Galaxy (mainly giant molecular clouds), and plausibly in late type galaxies generally, have 10^2<mu_2<10^8. MPs thus shorten the relaxation timescale by 10^1-10^7 relative to 2-body relaxation by stars alone. We show this increases by 10-1000 the rate of \emph{large}-periapse interactions with the MBH, where loss-cone refilling by stellar 2-body relaxation is inefficient. We extend the Fokker-Planck loss-cone formalism to approximately account for relaxation by rare encounters with MPs. We show that binary stars--MBH exchanges driven by MPs can explain the origin of the young main sequence B stars that are observed very near the Galactic MBH, and increase by orders of magnitude the ejection rate of hyper-velocity stars. In contrast, the rate of \emph{small}-periapse interactions of single stars with the MBH, such as tidal disruption, is only increased by a factor of a few. We suggest that MP-driven relaxation plays an important role in the 3-body exchange capture of single stars on very tight orbits around the MBH. These captured stars may later be disrupted by the MBH via tidal orbital decay or direct scattering into the loss cone; captured compact objects may inspiral into the MBH by the emission of gravitational waves from zero-eccentricity orbits.

preprint2005arXiv

Formation of molecular hydrogen on analogues of interstellar dust grains: experiments and modelling

Molecular hydrogen has an important role in the early stages of star formation as well as in the production of many other molecules that have been detected in the interstellar medium. In this review we show that it is now possible to study the formation of molecular hydrogen in simulated astrophysical environments. Since the formation of molecular hydrogen is believed to take place on dust grains, we show that surface science techniques such as thermal desorption and time-of-flight can be used to measure the recombination efficiency, the kinetics of reaction and the dynamics of desorption. The analysis of the experimental results using rate equations gives useful insight on the mechanisms of reaction and yields values of parameters that are used in theoretical models of interstellar cloud chemistry.