We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

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.

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.5 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*.

The formation of the first planetesimals and the final growth of planetary cores likely 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 cut-off for the maximal pebble sizes. We find the dependence of the critical solids size cut-off on the radial separation, disk age, turbulence strength, and the grain-size composition of the boulders in the disk. This size cut-off shapes the size distribution of small pebbles and thereby helps catalyze the early onset of SI formation of planetesimals, and the later core accretion growth through pebble accretion.

The dynamical evolution of triple stellar systems could give rise to the formation of compact binaries and induce binary mergers. Common envelope (CE) evolution, which plays a major role in the evolution of compact binary systems, can similarly play a key role in the evolution of triples. Here we use hydrodynamical simulations coupled with few-body dynamics to provide the first detailed models of triple common envelope (TCE) evolution. We focus on the circumstellar case, where the envelope of an evolved giant engulfs a compact binary orbiting the giant (inner-binary), which then in-spirals into the core of the evolved star. Through our exploratory modeling we find several possible outcomes of such TCE: (1) The merger of the binary inside the third star's envelope; (2) The binary disruption of the in-spiraling binary following its plunge, leading to a chaotic triple dynamics of the stellar-core and the two components of the former disrupted binary. The chaotic evolution typically leads to the in-spiral and merger of at least one of the former binary components with the core, and sometimes to the ejection of the second, or alternatively its further now-binary common-envelope evolution. The in-spiral in TCE leads to overall slower in-spiral, larger mass ejection and the production of more aspherical remnant, compared with a corresponding binary case of similar masses, due to the energy/momentum extraction from the inner-binary. We expect TCE to play a key role in producing various types of stellar-mergers and unique compact binary systems, and potentially induce transient electromagnetic and gravitational-wave sources.

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.

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.

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.

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.

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.

The coalescence of compact binary stars is expected to produce a stochastic background of gravitational waves (GW) observable with future GW detectors. Such backgrounds are usually characterized by their power spectrum as a function of frequency. Here, we present a method to calculate the full 1-point distribution of strain fluctuations. We focus on time series data, but our approach generalizes to the frequency domain. We illustrate how this probability distribution can be evaluated numerically. In addition, we derive accurate analytical asymptotic expressions for the large strain tail, which demonstrate that it is dominated by the nearest source. As an application, we also calculate the distribution of strain fluctuations for the astrophysical GW background produced by binary mergers of compact stars in the Universe, and quantify the extent to which it deviates from a Gaussian distribution. Our approach could be useful for the spectral shape reconstruction of stochastic GW backgrounds.

We introduce a new hybrid method to perform high-resolution tidal disruption simulations, at arbitrary orbits. An SPH code is used to simulate tidal disruptions only in the immediate spatial domain of the star, namely, where the tidal forces dominate over gravity, and then during the fragmentation phase in which the emerging tidal stream may collapse under its own gravity to form fragments. Following each hydrodynamical simulation, an analytical treatment is then applied to instantaneously transfer each fragment back to the tidal sphere for its subsequent disruption, in an iterative process. We validate the hybrid model by comparing it to both an analytical impulse approximation model of single tidal disruptions, as well as full-scale SPH simulations spanning the entire disc formation. The hybrid simulations are essentially indistinguishable from the full-scale SPH simulations, while computationally outperforming their counterparts by orders of magnitude. Thereby our new hybrid approach uniquely enables us to follow the long-term formation and continuous tidal disruption of the planet/planetesimal debris, without the resolution and orbital configuration limitation of previous studies. In addition, we describe a variety of future directions and applications for our hybrid model, which is in principle applicable to any star, not merely white dwarfs

We make use of a new hybrid method to simulate the long-term, multiple-orbit disc formation through tidal disruptions of rocky bodies by white dwarfs, at high-resolution and realistic semi-major axis. We perform the largest-yet suite of simulations for dwarf and terrestrial planets, spanning four orders of magnitude in mass, various pericentre distances and semi-major axes between 3 AU and 150 AU. This large phase space of tidal disruption conditions has not been accessible through the use of previous codes. We analyse the statistical and structural properties of the emerging debris discs, as well as the ejected unbound debris contributing to the population of interstellar asteroids. Unlike previous tidal disruption studies of small asteroids which form ring-like structures on the original orbit, we find that the tidal disruption of larger bodies usually forms dispersed structures of interlaced elliptic eccentric annuli on tighter orbits. We characterize the (typically power-law) size-distribution of the ejected interstellar bodies as well as their composition, rotation velocities and ejection velocities. We find them to be sensitive to the depth (impact parameter) of the tidal disruption. Finally, we briefly discuss possible implications of our results in explaining the peculiar variability of Tabby's star, the origin of the transit events of ZTF J0139+5245 and the formation of a planetary core around SDSS J1228+1040.

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.

Studies of the Galactic Centre suggest that in-situ star formation may have given rise to the observed stellar population near the central supermassive black hole (SMBH). Direct evidence for a recent starburst is provided by the currently observed young stellar disc (2-7 Myr) in the central 0.5 pc of the Galaxy. This result suggests that star formation in galactic nuclei may occur close to the SMBH and produce initially flattened stellar discs. Here we explore the possible build-up and evolution of nuclear stellar clusters near SMBHs through in-situ star formation producing stellar discs similar to those observed in the Galactic Centre and other nuclei. We make use of N-body simulations to model the evolution of multiple young stellar discs and explore the potential observable signatures imprinted by such processes. Each of the five simulated discs is evolved for 100 Myr before the next one is introduced in the system. We find that populations born at different epochs show different morphologies and kinematics. Older and presumably more metal poor populations are more relaxed and extended, while younger populations show a larger amount of rotation and flattening. We conclude that star formation in central discs can reproduce the observed properties of multiple stellar populations in galactic nuclei differing in age, metallicity and kinematic properties.

Type Ia supernovae (SNe) are thought to originate from the thermonuclear explosions of carbon-oxygen (CO) white dwarfs (WDs). The proposed progenitors of standard type Ia SNe have been studied for decades and can be, generally, divided into explosions of CO WDs accreting material from stellar non-degenerate companions (single-degenerate; SD models), and those arising from the explosive interaction of two CO WDs (double-degenerate; DD models). However, current models for the progenitors of such SNe fail to reproduce the diverse properties of the observed explosions, nor do they explain the inferred rates and the characteristics of the observed populations of type Ia SNe and their expected progenitors. Here we show that the little-studied mergers of CO-WDs with hybrid Helium-CO (He-CO) WDs can provide for a significant fraction of the normal type Ia SNe. Here we use detailed thermonuclear-hydrodynamical and radiative-transfer models to show that a wide range of mergers of CO WDs with hybrid He-CO WDs can give rise to normal type Ia SNe. We find that such He-enriched mergers give rise to explosions for which the synthetic light-curves and spectra resemble those of observed type Ia SNe, and in particular, they can produce a wide range of peak-luminosities, MB(MR)~ 18.4 to 19.2 (~ 18.5 to 19:45), consistent with those observed for normal type Ia SNe. Moreover, our population synthesis models show that, together with the contribution from mergers of massive double CO-WDs (producing the more luminous SNe), they can potentially reproduce the full range of type Ia SNe, their rate and delay-time distribution.

Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies, but for the first time perform an extensive statistical study of BH-BH, NS-NS and BH-NS binary mergers by means of direct high-precision regularized $N$-body simulations, including Post-Newtonian (PN) terms up to order PN2.5. We consider different SMBH masses, slopes for the BH mass function, binary semi-major axis and eccentricity distributions, and different spatial distributions for the binaries. We find that the merger rates are a decreasing function of the SMBH mass and are in the ranges $\sim 0.17$-$0.52 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, $\sim 0.06$-$0.10 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ and $\sim 0.04$-$0.16 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ for BH-BH, BH-NS and NS-NS binaries, respectively. However, the rate estimate from this channel remains highly uncertain and depends on the specific assumptions regarding the star-formation history in galactic nuclei and the supply rate of compact objects. We find that $\sim 10\%$--$20\%$ of the mergers enter the LIGO band with eccentricities $\gtrsim 0.1$. We also compare our results to the secular approximation, and show that $N$-body simulations generally predict a larger number of mergers. Finally, these events can also be observable via their electromagnetic counterparts, thus making these compact object mergers especially valuable for cosmological and astrophysical purposes.

The origin of ultra-wide massive binaries (orbital separations $10^3-2\times 10^5$~AU) and their properties are not well characterized nor understood. Here we use the second Gaia data release to search for wide astrometric companions to Galactic O-B5 stars which share similar parallax and proper motion with the primaries. Using the data we characterize the frequency and properties of such binaries. We find an ultra-wide multiplicity fraction of $4.4\pm0.5$ per cent, to our completeness limit (up to $\approx 17$~mag; down to G-stars at distances of 0.3-2~kpc, excluding stars in clusters). The secondary mass-function is generally consistent with a Kroupa initial stellar function; if extrapolated to lower mass companion stars we then might expect a wide-binary fraction of $\sim 27\pm5\%$. In addition we use these data as a verification sample to test the existence of ultra-wide binaries among neutron stars (NSs) and black holes (BHs). We propose that the discovery of such binary can provide unique constraints on the weakest natal kicks possible for NSs/BHs. If a compact object is formed in an ultra-wide binary and receives a very-low natal kick, such a binary should survive as a common proper motion pair. We therefore use Gaia data to search for ultra-wide companions to pulsars (normal and millisecond ones) and X-ray binaries. We find no reliable pairs. Future data could potentially provide stringent constraints through this method.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

The hypothesis of lunar origin by a single giant impact can explain some aspects of the Earth-Moon system. However, it is difficult to reconcile giant impact models with the compositional similarity of the Earth and Moon without violating angular momentum constraints. Furthermore, successful giant impact scenarios require very specific conditions such that they have a low probability of occurring. Here we present numerical simulations suggesting that the Moon could instead be the product of a succession of a variety of smaller collisions. In this scenario, each collision forms a debris disk around the proto-Earth that then accretes to form a moonlet. The moonlets tidally advance outward, and may coalesce to form the Moon. We find that sub-lunar moonlets are a common result of impacts expected onto the proto-Earth in the early solar system and find that the planetary rotation is limited by impact angular momentum drain. We conclude that, assuming efficient merger of moonlets, a multiple impact scenario can account for the formation of the Earth-Moon system with its present properties.

Planet formation begins with collisional growth of small planetesimals accumulating into larger ones. Such growth occurs while planetesimals are embedded in a gaseous protoplanetary disc. However, small-planetesimals experience collisions and gas-drag that lead to their destruction on short timescales, not allowing, or requiring fine tuned conditions for the efficient growth of $\sim$ metre-size objects. Here we show that $\sim10^{4}$ interstellar objects such as the recently detected 1I/2017-U1 ('Oumuamua) could have been captured, and become part of the young Solar System, together with up to hundreds of $\sim$ km sized ones. The capture rates are robust even for conservative assumptions on the protoplanetary disc structure, local stellar environment and planetesimal ISM density. 'Seeding' of such planetesimals then catalyze further planetary growth into planetary embryos, and potentially alleviate the main-challenges with the meter-size growth-"barrier". The capture model is in synergy with the current leading planet formation theories, providing the missing link to the first planetesimals. Moreover, planetesimal capture provides a far more efficient route for lithopanspermia than previously thought.

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.

Neutron-star (NS) - white-dwarf (WD) mergers may give rise to observable explosive transients, but have been little explored. We use 2D coupled hydrodynamical-thermonuclear FLASH-code simulations to study the evolution of WD debris-disks formed following WD-disruptions by NSs. We use a 19-elements nuclear-network and a detailed equation-of-state to follow the evolution, complemented by a post-process analysis using a larger 125-isotopes nuclear-network. We consider a wide range of initial conditions and study the dependence of the results on the NS/WD masses ($1.4-2{\rm M_{\odot}}$;$\,{\rm 0.375-0.7\,M_{\odot}}$, respectively), WD-composition (CO/He/hybrid-He-CO) and the accretion-disk structure. We find that viscous inflow in the disk gives rise to continuous wind-outflow of mostly C/O material mixed with nuclear-burning products arising from a weak detonation occurring in the inner-region of the disk. We find that such transients are energetically weak ($10^{48}-10^{49}$ergs) compared with thermonuclear-supernovae (SNe), and are dominated by the (gravitational) accretion-energy. Although thermonuclear-detonations occur robustly in all of our simulations (besides the He-WD) they produce only little energy $(1-10\%$ of the kinetic energy) and $^{56}{\rm Ni}$ ejecta (few$\times10^{-4}-10^{-3}{\rm M_{\odot}})$, with overall low ejecta masses of $\sim0.01-0.1{\rm M_{\odot}}$. Such explosions may produce rapidly-evolving transients, much shorter and fainter than regular type-Ia SNe. The composition and demographics of such SNe appear to be inconsistent with those of Ca-rich type Ib SNe. Though they might be related to the various classes of rapidly evolving SNe observed in recent years, they are likely to be fainter than the typical ones, and may therefore give rise a different class of potentially observable transients.

Common envelope (CE) is an important phase in the evolution of interacting evolved binary systems. The interaction of the binary components during the CE evolution (CEE) stage gives rise to orbital inspiral and the formation of a short-period binary or a merger, on the expense of extending and/or ejecting the envelope. CEE is not well understood, as hydrodynamical simulations show that only a fraction of the CE-mass is ejected during the dynamical inspiral, in contrast with observations of post-CE binaries. Different CE models suggest different timescales are involved in the CE-ejection, and hence a measurement of the CE-ejection timescale could provide direct constraints on the CEE-process. Here we propose a novel method for constraining the mass-loss timescale of the CE, using post-CE binaries which are part of wide-orbit triple systems. The orbit/existence of a third companion constrains the CE mass-loss timescale, since rapid CE mass-loss may disrupt the triple system, while slower CE mass-loss may change the orbit of the third companion without disrupting it. As first test-cases we examine two observed post-CE binaries in wide triples, Wolf-1130 and GD-319. We follow their evolution due to mass-loss using analytic and numerical tools, and consider different mass-loss functions. We calculate a wide grid of binary parameters and mass-loss timescales in order to determine the most probable mass-loss timescale leading to the observed properties of the systems. We find that mass-loss timescales of the order of $10^{3}-10^{5}{\rm yr}$ are the most likely to explain these systems. Such long timescales are in tension with most of the CE mass-loss models, which predict shorter, dynamical timescales, but are potentially consistent with the longer timescales expected from the dust-driven winds model for CE ejection.

Important clues on the formation and evolution of planetary systems can be inferred from the stellar obliquity $\psi$. We study the distribution of obliquities using the California-Kepler Survey and the TEPCat Catalog of Rossiter-MacLaughlin (RM) measurements, from which we extract, respectively, 275 and 118 targets. We infer a "best fit" obliquity distribution in $\psi$ with a single parameter $\kappa$. Large values of $\kappa$ imply that $\psi$ is distributed narrowly around zero, while small values imply approximate isotropy. Our findings are: (1) The distribution of $\psi$ in Kepler systems is narrower than found by previous studies and consistent with $\kappa\sim15$ (mean $\langle\psi\rangle\sim19^\circ$ and spread $\psi\sim10^\circ$). (2) The value of $\kappa$ in Kepler systems does not depend, at a statistically significant level, on planet multiplicity, stellar multiplicity or stellar age; on the other hand, metal-rich hosts, small planet hosts and long-period planet hosts tend to be more oblique than the general sample (at a $\sim$2.5-$\sigma$ significance level). (3) Hot Jupiter (HJ) systems with RM measurements are consistent with $\kappa\sim2$, more broadly distributed than the general Kepler population. (4) A separation of the RM sample into cooler ($T_{\rm eff}\lesssim$6250 K) and hotter ($T_{\rm eff}\gtrsim$6250 K) HJ hosts results in two distinct distributions, $\kappa_{\rm cooler}\sim4$ and $\kappa_{\rm hotter}\sim1$ (4-$\sigma$ significance), both more oblique than the Kepler sample. We hypothesize that the total mass in planets may be behind the increasing obliquity with metallicity and planet radius, and that the period dependence could be due to primordial disk alignment rather than realignment of stellar spin.

The mergers of neutron stars (NSs) and white dwarfs (WDs) could give rise to explosive transients, potentially observable with current and future transient surveys. However, the expected properties and distribution of such events is not well understood. Here we characterize the rates of such events, their delay time distribution, their progenitors and the distribution of their properties. We use binary populations synthesis models and consider a wide range of initial conditions and physical processes. In particular we consider different common-envelope evolution models and different NS natal kick distributions. We provide detailed predictions arising from each of the models considered. We find that the majority of NS-WD mergers are born in systems in which mass-transfer played an important role, and the WD formed before the NS. For the majority of the mergers the WDs have a carbon-oxygen composition (60-80%) and most of the rest are with oxygen-neon WDs. The rates of NS-WD mergers are in the range of 3-15% of the type Ia supernovae (SNe) rate. Their delay time distribution is very similar to that of type Ia SNe, but slightly biased towards earlier times. They typically explode in young 0.1-1Gyr environments, but have a tail distribution extending to long, Gyrs-timescales. Models including significant kicks give rise to relatively wide offset distribution extending to hundreds of kpcs. The demographic and physical properties of NS-WD mergers suggest they are likely to be peculiar type Ic-like SNe, mostly exploding in late type galaxies. Their overall properties could be related to a class of rapidly evolving SNe recently observed, while they are less likely to be related to the class of Ca-rich SNe.

In hierarchical triple systems, the inner binary is perturbed by a distant companion. For large mutual inclinations, the Lidov-Kozai mechanism secularly excites large eccentricity and inclination oscillations of the inner binary. The maximal eccentricity attained, $e_{\rm max}$ is simply derived and widely used. However, for mildly hierarchical systems (i.e. the companion is relatively close and massive), non-secular perturbations affect the evolution. Here we account for fast non-secular variations and find new analytic formula for $e_{{\rm max}}$, in terms of the system's hierarchy level, correcting previous work and reproducing the orbital flip criteria. We find that $e_{{\rm max}}$ is generally enhanced, allowing closer encounters between the inner binary components, thus significantly changing their interaction and its final outcome. We then extend our approach to include additional relativistic and tidal forces. Using our results, we show that the merger time of gravitational-wave (GW) sources orbiting massive black-holes in galactic nuclei is enhanced compared with previous analysis accounting only for the secular regime. Consequently, this affects the distribution and rates of such GW sources in the relevant mild-hierarchy regime. We test and confirm our predictions with direct N-body and 2.5-level Post-Newtonian codes. Finally, we calculate the formation and disruption rates of hot-Jupiters (HJ) in planetary systems using a statistical approach, which incorporates our novel results for $e_{{\rm max}}$. We find that more HJ migrate from further out, but they are also tidally disrupted more frequently. Remarkably, the overall formation rate of HJs remains similar to that found in previous studies. Nevertheless, the different rates could manifest in different underlying distribution of observed warm-Jupiters.

The Earth-Moon system is suggested to have formed through a single giant collision, in which the Moon accreted from the impact-generated debris disk. However, such giant impacts are rare, and during its evolution the Earth experienced many more smaller impacts, producing smaller satellites that potentially coevolved. In the multiple-impact hypothesis of lunar formation, the current Moon was produced from the mergers of several smaller satellites (moonlets), each formed from debris disks produced by successive large impacts. In the Myrs between impacts, a pre-existing moonlet tidally evolves outward until a subsequent impact forms a new moonlet, at which point both moonlets will tidally evolve until a merger or system disruption. In this work, we examine the likelihood that pre-existing moonlets survive subsequent impact events, and explore the dynamics of Earth-moonlet systems that contain two moonlets generated Myrs apart. We demonstrate that pre-existing moonlets can tidally migrate outward, remain stable during subsequent impacts, and later merge with newly created moonlets (or re-collide with the Earth). Formation of the Moon from the mergers of several moonlets could therefore be a natural byproduct of the Earth's growth through multiple impacts. More generally, we examine the likelihood and consequences of Earth having prior moons, and find that the stability of moonlets against disruption by subsequent impacts implies that several large impacts could post-date Moon formation.

Pebble accretion is a promising process for decreasing growth timescales of planetary cores, allowing gas giants to form at wide orbital separations. However, nebular turbulence can reduce the efficiency of this gas-assisted growth. We present an order of magnitude model of pebble accretion, which calculates the impact of turbulence on the average velocity of small bodies, the radius for binary capture, and the sizes of the small bodies that can be accreted. We also include the effect of turbulence on the particle scale height, which has been studied in previous works. We find that turbulence does not prevent rapid growth in the high-mass regime: the last doubling time to the critical mass to trigger runaway gas accretion ($M \sim 10 \, M_\oplus$) is well within the disk lifetime even for strong ($\alpha \gtrsim 10^{-2}$) turbulence. We find that while the growth timescale is quite sensitive to the local properties of the protoplanetary disk, over large regimes of parameter space large cores grow in less than the disk lifetime if appropriately-sized small bodies are present. Instead, the effects of turbulence are most pronounced for low planetary masses. For strong turbulence the growth timescale is longer than the gas disk lifetime until the core reaches masses $\gtrsim 10^{-2}-10^{-1} M_\oplus$. A "Flow Isolation Mass", at which binary capture ceases, emerges naturally from our model framework. We comment that the dependence of this mass on orbital separation is similar to the semi-major axis distribution of solar-system cores.

During the last stages of the terrestrial planet formation, planets grow mainly through giant-impacts with large planetary embryos. The Earth's Moon was suggested to form through one of these impacts. However, since the proto-Earth has experienced many giant-impacts, several moons are naturally expected to form through a sequence of multiple (including smaller scale) impacts. Each impact potentially forms a sub-Lunar mass moonlet that interacts gravitationally with the proto-Earth and possibly with previously-formed moonlets. Such interactions result in either moonlet-moonlet mergers, moonlet ejections or infall of moonlets on the Earth. The latter possibility, leading to low-velocity moonlet-Earth collisions is explored here for the first time. We make use of SPH simulations and consider a range of moonlet masses, collision impact-angles and initial proto-Earth rotation rates. We find that grazing/tidal-collisions are the most frequent and produce comparable fractions of accreted-material and debris. The latter typically clump in smaller moonlets that can potentially later interact with other moonlets. Other collision geometries are more rare. Head-on collisions do not produce much debris and are effectively perfect mergers. Intermediate impact angles result in debris mass-fractions in the range of 2-25% where most of the material is unbound. Retrograde collisions produce more debris than prograde collisions, whose fractions depend on the proto-Earth initial rotation rate. Moonfalls can slightly change the rotation-rate of the proto-Earth. Accreted moonfall material is highly localized, potentially explaining the isotopic heterogeneities in highly siderophile elements in terrestrial rocks, and possibly forming primordial super-continent topographic features. Our results can be used for simple scaling laws and applied to n-body studies of the formation of the Earth and Moon.

Common-envelope evolution (CEE) is the short-lived phase in the life of an interacting binary-system during which two stars orbit inside a single shared envelope. Such evolution is thought to lead to the inspiral of the binary, the ejection of the extended envelope and the formation of a remnant short-period binary. However, detailed hydrodynamical models of CEE encounter major difficulties. They show that following the inspiral most of the envelope is not ejected; though it expands to larger separations, it remains bound to the binary. Here we propose that dust-driven winds can be produced following the CEE. These can evaporate the envelope following similar processes operating in the ejection of the envelopes of AGB stars. Pulsations in an AGB-star drives the expansion of its envelope, allowing the material to cool down to low temperatures thus enabling dust condensation. Radiation pressure on the dust accelerates it, and through its coupling to the gas it drives winds which eventually completely erode the envelope. We show that the inspiral phase in CE-binaries can effectively replace the role of stellar pulsation and drive the CE expansion to scales comparable with those of AGB stars, and give rise to efficient mass-loss through dust-driven winds.

White dwarfs (WDs) are the stellar core remnants of low mass stars. They are typically divided into three main composition groups: Oxygen Neon (ONe), Carbon Oxygen (CO) and Helium (He) WDs. The evolution of binary systems can significantly change the evolution of the binary stellar components. In particular, striping the envelope of an evolved star can give rise to a core remnant, which can later evolve into a WD with significantly different composition. Here we focus on the formation and evolution of hybrid HeCO WDs. We follow the formation and stellar evolution of such WDs for a range of initial conditions and provide their detailed structure, mass-radius relation and luminosity temperature evolution. We find that both low mass WDs (< 0.45M , typically thought to be He WDs) and intermediate-mass WDs (0.45 < MWD < 0.7, typically thought to be CO WDs) could in fact be hybrid HeCO WDs, with 5-25 (75 -95)% of their mass in He (CO). We use population synthesis calculations to infer the birth rate and properties of such WDs. We find that hybrid HeCO WD comprise the majority of young (< 2Gyr) WDs in binaries, but are more rare among older WDs in binaries. The high frequency and large He content of such WDs could have an important role in WD WD mergers, and may give rise to sub Chandrasekhar thermonuclear supernova explosions.

The Galactic center of the Milky-Way harbors a massive black hole (BH) orbited by a diverse population of young and old stars. A significant fraction of the youngest stars ($\sim4-7$ Myr) reside in a thin stellar disk with puzzling properties; the disk appears to be warped, shows asymmetries, and contains one or more clumpy structures (e.g. IRS 13). Models explaining the clumping invoked the existence of an intermediate-mass BH of $10^{3}-10^{4}$ M$_{\odot}$, but no kinematic evidence for such a BH has been found. Here we use extended $N$-body simulations and hybrid self-consistent field method models to show that naturally formed residual temporal asphericity of the hosting nuclear star cluster gives rise to torques on the disk, which lead to changes in its orientation over time, and to recurrent formation and dissolution of single spiral arm ($m=1$ modes) structures. The changing orientation leads to a flapping-like behavior of the disk and to the formation of a warped disk structure. The spiral arms may explain the over-densities in the disk (clumping) and its observed asymmetry, without invoking the existence of an intermediate-mass BH. The spiral arms are also important for the overall disk evolution and can be used to constrain the structure and composition of the nuclear stellar cluster.

UV and optically selected candidates for stellar tidal disruption events (TDE) often exhibit broad spectral features (HeII emission, H$\alpha$ emission, or absorption lines) on a blackbody-like continuum (1e4K<T<1e5K). The lines presumably emit from TDE debris or circumnuclear clouds photoionized by the flare. Line velocities however are much lower than expected from a stellar disruption by supermassive black hole (SMBH), and are somewhat faster than expected for the broad line region (BLR) clouds of a persistently active galactic nucleus (AGN). The distinctive spectral states are not strongly related to observed luminosity and velocity, nor to SMBH mass estimates. We use exhaustive photoionization modelling to map the domain of fluxes and cloud properties that yield (e.g.) a He-overbright state where a large HeII(4686A)/H$\alpha$ line-ratio creates an illusion of helium enrichment. Although observed line ratios occur in a plausible minority of cases, AGN-like illumination can not reproduce the observed equivalent widths. We therefore propose to explain these properties by a light-echo photoionization model: the initial flash of a hot blackbody (detonation) excites BLR clouds, which are then seen superimposed on continuum from a later, expanded, cooled stage of the central luminous source. The implied cloud mass is substellar, which may be inconsistent with a TDE. Given these and other inconsistencies with TDE models (e.g. host-galaxies distribution) we suggest to also consider alternative origins for these nuclear flares, which we briefly discuss (e.g. nuclear supernovae and starved/subluminous AGNs).

Binary black-holes (BHs) and binary neutron-stars (NSs) mergers had been recently detected through gravitational-wave (GW) emission, with the latter followed by post-merger electromagnetic counterparts, appearing seconds up to weeks after the merger. While post-merger electromagnetic counterparts had been anticipated theoretically, very little electromagnetic precursors to GW-sources had been proposed, and non observed yet. Here we show that a fraction of ${\rm a}\,{\rm few\times}10^{-4}-10^{-1}$ of LIGO GW-sources and short-GRBs, could be preceded by supernovae-explosions years up to decades before the merger. Each of the BH/NS-progenitors in GW-sources are thought to form following a supernova, likely accompanied by a natal velocity-kick to the newly born compact object. The evolution and natal-kicks determine the orbits of surviving binaries, and hence the delay-time between the birth of the compact-binary and its final merger through GW-emission. We use data from binary evolution population-synthesis models to show that the delay-time distribution has a non-negligible tail of ultra-short delay-times between 1-100 yrs, thereby giving rise to potentially observable supernovae precursors to GW-sources. Moreover, future LISA/DECIGO GW space-detectors will enable the detection of GW-inspirals in the pre-mergers stage weeks to decades before the final merger. These ultra-short delay-time sources could therefore produce a unique type of promptly appearing LISA/DECIGO-GW-sources accompanied by \emph{coincident} supernovae. The archival (and/or direct) detection of precursor (coincident) SNe with GW and/or short-GRBs will provide unprecedented characterizations of the merging-binaries, and their prior evolution through supernovae and natal kicks, otherwise inaccessible through other means.

We study the chaotic and secular evolution of hierarchical quadruple systems in the $3+1$ configuration, focusing on the evolution of mutual inclination of the inner binaries as the system undergoes coupled Lidov-Kozai (LK) oscillations. We include short-range forces (SRF; such as those due to tidal and rotational distortions) that control the eccentricity excitation of the inner binary. The evolution of mutual inclination is described, a priori, by two dimensionless parameters, $\pazocal{R}_0$, the ratio between the inner and outer LK time-scales and $\epsilon_{SRF}$, the ratio between the SRF precession and the inner LK precession rates. We find that the chaotic zones for the mutual inclination depend mainly on $\pazocal{R}_0$, while $\epsilon_{SRF}$ controls mainly the range of eccentricity excitation. The mutual inclination evolves chaotically for $1\lesssim \pazocal{R}_0\lesssim 10$, leading to large misalignments. For $0.4 \lesssim \pazocal{R}_0 \lesssim 0.8$, the system could be weakly excited and produce bimodal distribution of mutual inclination angles. Our results can be applied to exomoons-planets in stellar binaries and Warm/Hot Jupiters in stellar triples. Such systems could develop large mutual inclination angles if the inner binary is tight enough, and also high eccentricities, depending of the strength of the short-range forces. Future detections of tilted Warm/Hot Jupiters and exomoons could put our mechanism under observational tests.

Type Ia supernovae (Ia-SNe) are thought to arise from the thermonuclear explosions of white dwarfs (WDs). The progenitors of such explosions are still highly debated; in particular the conditions leading to detonations in WDs are not well understood in most of the suggested progenitor models. Nevertheless, direct head-on collisions of two WDs were shown to give rise to detonations and produce Ia-SNe - like explosions, and were suggested as possible progenitors. The rates of such collisions in dense globular clusters are far below the observed rates of type Ia SNe, but it was suggested that quasi-secular evolution of hierarchical triples could produce a high rate of such collisions. Here we used detailed triple stellar evolution populations synthesis models coupled with dynamical secular evolution to calculate the rates of WD-WD collisions in triples and their properties. We explored a range of models with different realistic initial conditions and derived the expected SNe total mass, mass-ratio and delay time distributions for each of the models. We find that the SNe rate from WD-WD collisions is of the order of 0.1% of the observed Ia-SNe rate across all our models, and the delay-time distribution is almost uniform in time, and is inconsistent with observations. We conclude that SNe from WD-WD collisions in isolated triples can at most provide for a small fraction of Ia-SNe, and can not serve as the main progenitors of such explosions.

Studies suggest that the pollution of white dwarf (WD) atmospheres arises from the accretion of minor planets, but the exact properties of polluting material, and in particular the evidence for water in some cases are not yet understood. Previous works studied the water retention in minor planets around main-sequence and evolving host stars, in order to evaluate the possibility that water survives inside minor planets around WDs. However, all of these studies focused on small, comet-sized to moonlet-sized minor planets, when the inferred mass inside the convection zones of He-dominated WDs could actually also be compatible with much more massive minor planets. In this study we therefore explore for the first time, the water retention inside exo-planetary dwarf planets, or moderate-sized moons, with radii of the order of hundreds of kilometres. We now cover nearly the entire potential mass range of minor planets. The rest of the parameter space considered in this study is identical to that of our previous study, and also includes multiple WD progenitor star masses. We find that water retention in more massive minor planets is still affected by the mass of the WD progenitor, however not as much as when small minor planets were considered. We also find that water retention is now almost always greater than zero. On average, the detected water fraction in He-dominated WD atmospheres should be at least 5\%, irrespective of the assumed initial water composition, if it came from a single accretion event of an icy dwarf planet or moon. This finding also strengthens the possibility of WD habitability. To finalize our previous and current findings, we provide a code which may be freely used as a service to the community. The code calculates ice and water retention by interpolation, spanning the full mass range of both minor planets and their host stars.

Nuclear spiral arms are small-scale transient spiral structures found in the centers of galaxies. Similarly to their galactic-scale counterparts, nuclear spiral arms can perturb the orbits of stars. In the case of the Galactic Center (GC), these perturbations can affect the orbits of stars and binaries in a region extending to several hundred parsecs around the supermassive black hole (MBH), causing diffusion in orbital energy and angular momentum. This diffusion process can drive stars and binaries to close approaches with the MBH, disrupting single stars in tidal disruption events (TDEs), or disrupting binaries, leaving a star tightly bound to the MBH, and an unbound star escaping the galaxy, i.e., a hypervelocity star (HVS). Here, we consider diffusion by nuclear spiral arms in galactic nuclei, specifying to the Milky Way GC. We determine nuclear spiral arm-driven diffusion rates using test-particle integrations, and compute disruption rates. Our TDE rates are up to 20% higher compared to relaxation by single stars. For binaries, the enhancement is up to a factor of ~100, and our rates are comparable to the observed numbers of HVSs and S-stars. Our scenario is complementary to relaxation driven by massive perturbers. In addition, our rates depend on the inclination of the binary with respect to the Galactic plane. Therefore, our scenario provides a novel potential source for the observed anisotropic distribution of HVSs. Nuclear spiral arms may also be important for accelerating the coalescence of binary MBHs, and for supplying nuclear star clusters with stars and gas.

Most observations of polluted white dwarf atmospheres are consistent with accretion of water depleted planetary material. Among tens of known cases, merely two cases involve accretion of objects that contain a considerable mass fraction of water. The purpose of this study is to investigate the relative scarcity of these detections. Based on a new and highly detailed model, we evaluate the retention of water inside icy minor planets during the high luminosity stellar evolution that follows the main sequence. Our model fully considers the thermal, physical, and chemical evolution of icy bodies, following their internal differentiation as well as water depletion, from the moment of their birth and through all stellar evolution phases preceding the formation of the white dwarf. We also account for different initial compositions and formation times. Our results show that previous studies have either underestimated or overestimated water retention. We also reaffirm that water can survive in a variety of circumstances and in great quantities, and therefore other possibilities are discussed in order to explain the infrequency of water detections. We find that the most likely explanation is that water does not survive the dissociation of tidally disrupted minor planets, and instead a direct impact onto the white dwarf is required. We show that this possibility is statistically plausible, and predict that if this hypothesis is correct, future observations of extremely water rich atmospheres will involve only helium dominated DB white dwarfs.

Most studies suggest the pollution of white dwarf (WD) atmospheres arises from accretion of minor planets, but the exact properties of polluting material, and in particular the evidence for water in some cases are not yet understood. Here we study the water retention of small icy bodies in exo-solar planetary systems, as their respective host stars evolve through and off the main sequence and eventually become WDs. We explore for the first time a wide range of star masses and metallicities. We find that the mass of the WD progenitor star is of crucial importance for the retention of water, while its metallicity is relatively unimportant. We predict that minor planets around lower-mass WD progenitors would retain more water in general, and would do so at closer distances from the WD, compared with high-mass progenitors. The dependence of water retention on progenitor mass and other parameters has direct implications for the origin of observed WD pollution, and we discuss how our results and predictions might be tested in the future as more observations of WDs with long cooling-ages become available.

The New Horizons mission to the Kuiper Belt has recently revealed intriguing features on the surface of Charon, including a network of chasmata, cutting across or around a series of high topography features, conjoining to form a belt. It is proposed that this tectonic belt is a consequence of contraction/expansion episodes in the moon's evolution associated particularly with compaction, differentiation and geochemical reactions of the interior. The proposed scenario involves no need for solidification of a vast subsurface ocean and/or a warm initial state. This scenario is based on a new, detailed thermo-physical evolution model of Charon that includes multiple processes. According to the model, Charon experiences two contraction/expansion episodes in its history that may provide the proper environment for the formation of the tectonic belt. This outcome remains qualitatively the same, for several different initial conditions and parameter variations. The precise orientation of Charon's tectonic belt, and the cryovolcanic features observed south of the tectonic belt may have involved a planetary-scale impact, that occurred only after the belt had already formed.

[abridged] In the typical giant-impact scenario for the Moon formation most of the Moon's material originates from the impactor. Any Earth-impactor composition difference should, therefore, correspond to a comparable Earth-Moon composition difference. Analysis of Moon rocks shows a close Earth-Moon composition similarity, posing a challenge for the giant-impact scenario, given that impactors were thought to significantly differ in composition from the planets they impact. Here we use a large set of 140 simulations to show that the composition of impactors could be very similar to that of the planets they impact; in $4.9\%$-$18.2\%$ ($1.9\%$-$6.7\%$) of the cases the resulting composition of the Moon is consistent with the observations of $\Delta^{17}O<15$ ($\Delta^{17}O<6$ ppm). These findings suggest that the Earth-Moon composition similarity could be resolved as to arise from the primordial Earth-impactor composition similarity. Note that although we find the likelihood for the suggested competing model of very high mass-ratio impacts (producing significant Earth-impactor composition mixing) is comparable ($<6.7\%$), this scenario also requires additional fine-tuned requirements of a very fast spinning Earth. Using the same simulations we also explore the composition of giant-impact formed Mars-moons as well as Vesta-like asteroids. We find that the Mars-moon composition difference should be large, but smaller than expected if the moons are captured asteroids. Finally, we find that the left-over planetesimals ('asteroids') in our simulations are frequently scattered far away from their initial positions, thus potentially explaining the mismatch between the current position and composition of the Vesta asteroid.

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.

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.

The 2175 \AA\ UV extinction feature was discovered in the mid-1960s, yet its physical origin remains poorly understood. One suggestion is absorption by Polycyclic Aromatic Hydrocarbons (PAH) molecules, which is supported by theoretical molecular structure computations and by laboratory experiments. PAHs are positively detected by their 3.3, 6.2, 7.7 8.6, 11.3 & 12.7 $\mu$m IR emission bands, which are specified by their modes of vibration. A definitive empirical link between the 2175 \AA\ UV extinction, and the PAH IR emission bands, however, is still missing. We present a new sample of hot stars that have both 2175 \AA\ absorption and PAH IR emission. We find significant shifts of the central wavelength of the UV absorption feature, up to 2350 \AA, but predominantly in stars that also have IR PAH emission. These UV shifts depend on stellar temperature in a fashion that is similar to the shifts of the 6.2 and 7.7$\mu$m PAH IR bands, namely the features are increasingly more red-shifted as the stellar temperature decreases, but only below $\sim 15$ kK. Above 15 kK both UV and IR features retain their nominal values. Moreover, we find a suggestive correlation between the UV and IR shifts. We hypothesize that these similar dependences of both the UV and IR features on stellar temperature hint to a common origin of the two in PAH molecules and may establish the missing link between the UV and IR observations. We further suggest that the shifts depend on molecular size, and that the critical temperature of $\sim 15$ kK above which no shifts are observed is related to the onset of UV driven hot-star winds and their associated shocks.

Nuclear Star Clusters (NSCs) are commonly observed in the centres of most galactic nuclei, including our own Milky Way. While their study can reveal important information about the build-up of the innermost regions of galaxies, the physical processes that regulate their formation are still poorly understood. NSCs might have been formed through gas infall and subsequent in situ star formation, and/or through the infall and merging of multiple star clusters into the centre of the galaxy. Here, we investigate the viability of the latter, by studying direct N-body simulations of inspiralling clusters to the centre of a Milky-Way-like nuclear bulge that hosts a massive black hole. We find that the NSC that forms through this process can show both morphological and kinematical properties that make it comparable with observations of the Milky Way NSC, including significant rotation- a fact that has been so far attributed mainly to gas infall. We explore its kinematic evolution, to see if and how the merger history can imprint fossil records on its dynamical structure. Moreover, we study the effect of stellar foreground contamination in the line-of-sight kinematics of the NSC. Our study shows that no fine tuning of the orientation of the infalling globular clusters is necessary to result in a rotating NSC. We suggest that cluster-inspiral is a viable mechanism for the formation of rotating NSCs.

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; $\mu$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 $\mu$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.

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.

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.

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.

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).

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.

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.

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.

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.

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.

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.

The dwarf planet Pluto is known to host an extended system of five co-planar satellites. Previous studies have explored the formation and evolution of the system in isolation, neglecting perturbative effects by the Sun. Here we show that secular evolution due to the Sun can strongly affect the evolution of outer satellites and rings in the system, if such exist. Although precession due to extended gravitational potential from the inner Pluto-Charon binary quench such secular evolution up to $a_{crit}\sim0.0035$ AU ($\sim0.09$ $R_{Hill}$ the Hill radius; including all of the currently known satellites), outer orbits can be significantly altered. In particular, we find that \emph{co-planar} rings and satellites should not exist beyond $a_{crit}$; rather, satellites and dust particles in these regions secularly evolve on timescales ranging between $10^{4}-10^{6}$ yrs, and quasi-periodically change their inclinations and eccentricities through secular evolution (Lidov-Kozai oscillations). Such oscillations can lead to high inclinations and eccentricities, constraining the range where such satellites (and dust particles) can exist without crossing the orbits of the inner satellites, or crossing the outer Hill stability range. Outer satellites, if such exist are therefore likely to be \emph{irregular} satellites, with orbits limited to be non-circular and/or highly inclined. These could be potentially detected and probed by the New-Horizon mission, possibly providing direct evidence for the secular evolution of the Pluto satellite system, and shedding new light on its origins.

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.

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.

Nuclear stellar cluster (NSCs) are known to exist around massive black holes (MBHs) in galactic nuclei. Two formation scenarios were suggested for their origin: Build-up of NSCs and Continuous in-situ star-formation. Here we 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 show that continuous star-formation 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We present extensive multi-wavelength observations of the extremely rapidly declining Type Ic supernova, SN 2005ek. Reaching a peak magnitude of M_R = -17.3 and decaying by ~3 mag in the first 15 days post-maximum, SN 2005ek is among the fastest Type I supernovae observed to date. The spectra of SN 2005ek closely resemble those of normal SN Ic, but with an accelerated evolution. There is evidence for the onset of nebular features at only nine days post-maximum. Spectroscopic modeling reveals an ejecta mass of ~0.3 Msun that is dominated by oxygen (~80%), while the pseudo-bolometric light curve is consistent with an explosion powered by ~0.03 Msun of radioactive Ni-56. Although previous rapidly evolving events (e.g., SN 1885A, SN 1939B, SN 2002bj, SN 2010X) were hypothesized to be produced by the detonation of a helium shell on a white dwarf, oxygen-dominated ejecta are difficult to reconcile with this proposed mechanism. We find that the properties of SN 2005ek are consistent with either the edge-lit double detonation of a low-mass white dwarf or the iron-core collapse of a massive star, stripped by binary interaction. However, if we assume that the strong spectroscopic similarity of SN 2005ek to other SN Ic is an indication of a similar progenitor channel, then a white-dwarf progenitor becomes very improbable. SN 2005ek may be one of the lowest mass stripped-envelope core-collapse explosions ever observed. We find that the rate of such rapidly declining Type I events is at least 1-3% of the normal SN Ia rate.

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.

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.

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.

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.

(abridged) We consider the effect of stellar remnants on the interstellar medium of a massive star cluster following the initial burst of star formation. We argue that accretion onto stellar-mass black holes (BHs) is an effective mechanism for rapid gas depletion in clusters of all masses, as long as they contain progenitor stars more massive than \gtrsim 50\msun. This scenario is attractive for the progenitor systems of present-day massive globular clusters (GCs) which likely had masses M \gtrsim 10^7\msun. In such clusters, supernovae and stellar winds cannot provide a plausible explanation for the sudden removal of the primordial gas reservoir that is required to explain their complex chemical enrichment history. In order to consider different regimes in the gas accretion rate onto stellar-mass BHs, we consider both the Bondi-Hoyle and Eddington approximations. For either model, our results show that the gas can be significantly depleted within only a few tens of Myrs. This process will affect the distribution of BH masses, and may accelerate the dynamical decoupling of the BH population and, ultimately, their dynamical ejection. Moreover, the timescales for gas depletion are sufficiently short that the accreting BHs could significantly affect the chemistry of subsequent star formation episodes. The gas depletion times and final mass in BHs are sensitive to the assumed model for the accretion rate, and to the initial mass of the most massive BH which, in turn, is determined by the upper mass cut-off of the stellar IMF. Our results imply that the remnant accretion history can have an important bearing on the observed present-day cluster mass-to-light ratio. In particular, we show that an increase of the upper mass cut-off with decreasing metallicity could contribute to the observed anti-correlation between the mass-to-light ratio and the metallicity of GCs.

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).

(abridged) We use newly available empirical binary fractions for globular clusters to carry out a direct test of the binary evolution hypothesis, and of collisional channels that involve binary stars. More specifically, using the previously reported correlation between blue straggler numbers and core mass as a benchmark, we test for correlations with the number of binary stars, as well as with the rates of single-single, single-binary, and binary-binary encounters. Surprisingly, we find that the simple correlation with core mass remains by far the strongest predictor of blue straggler population size, even in our joint models. This is despite the fact that the binary fractions themselves strongly anti-correlate with core mass, just as expected in the binary evolution model. At first sight, these results do not fit neatly with either binary evolution or collisional models in their simplest forms. Arguably the simplest and most intriguing possibility to explain this unexpected result is that observational errors on the core binary fractions are larger than the true intrinsic dispersion associated with their dependence on core mass. In the context of the binary evolution model, this would explain why the combination of binary fraction and core mass is a poorer predictor of blue straggler numbers than core mass alone. It would also imply that core mass is a remarkably clean predictor of core binary fractions. This would be of considerable importance for the dynamical evolution of globular clusters, and provides an important benchmark for models attempting to understand their present-day properties.

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.

The environment near super massive black holes (SMBHs) in galactic nuclei contain a large number of stars and compact objects. A fraction of these are likely to be members of binaries. Here we discuss the binary population of stellar black holes and neutron stars near SMBHs and focus on the secular evolution of such binaries, due to the perturbation by the SMBH. Binaries with highly inclined orbits in respect to their orbit around the SMBH are strongly affected by secular Kozai processes, which periodically change their eccentricities and inclinations (Kozai-cycles). During periapsis approach, at the highest eccentricities during the Kozai-cycles, gravitational wave emission becomes highly efficient. Some binaries in this environment can inspiral and coalesce at timescales much shorter than a Hubble time and much shorter than similar binaries which do not reside near a SMBH. The close environment of SMBHs could therefore serve as catalyst for the inspiral and coalescence of binaries, and strongly affect their orbital properties. Such compact binaries would be detectable as gravitational wave (GW) sources by the next generation of GW detectors (e.g. advanced- LIGO). About 0.5% of such nuclear merging binaries will enter the LIGO observational window while on orbit that are still very eccentric (e>~0.5). The efficient gravitational wave analysis for such systems would therefore require the use of eccentric templates. We also find that binaries very close to the MBH could evolve through a complex dynamical (non-secular) evolution leading to emission of several GW pulses during only a few yrs (though these are likely to be rare). Finally, we note that the formation of close stellar binaries, X-ray binaries and their merger products could be induced by similar secular processes, combined with tidal friction rather than GW emission as in the case of compact object binaries.

Here we propose a mechanism for efficiently growing intermediate mass black holes (IMBH) in disks around supermassive black holes. Stellar mass objects can efficiently agglomerate when facilitated by the gas disk. Stars, compact objects and binaries can migrate, accrete and merge within disks around supermassive black holes. While dynamical heating by cusp stars excites the velocity dispersion of nuclear cluster objects (NCOs) in the disk, gas in the disk damps NCO orbits. If gas damping dominates, NCOs remain in the disk with circularized orbits and large collision cross-sections. IMBH seeds can grow extremely rapidly by collisions with disk NCOs at low relative velocities, allowing for super-Eddington growth rates. Once an IMBH seed has cleared out its feeding zone of disk NCOs, growth of IMBH seeds can become dominated by gas accretion from the AGN disk. However, the IMBH can migrate in the disk and expand its feeding zone, permitting a super-Eddington accretion rate to continue. Growth of IMBH seeds via NCO collisions is enhanced by a pile-up of migrators. We highlight the remarkable parallel between the growth of IMBH in AGN disks with models of giant planet growth in protoplanetary disks. If an IMBH becomes massive enough it can open a gap in the AGN disk. IMBH migration in AGN disks may stall, allowing them to survive the end of the AGN phase and remain in galactic nuclei. Our proposed mechanisms should be more efficient at growing IMBH in AGN disks than the standard model of IMBH growth in stellar clusters. Dynamical heating of disk NCOs by cusp stars is transferred to the gas in a AGN disk helping to maintain the outer disk against gravitational instability. Model predictions, observational constraints and implications are discussed in a companion paper (Paper II).

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.

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.

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.

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.

Galactic nuclei should contain a cluster of stars and compact objects in the vicinity of the central supermassive black hole due to stellar evolution, minor mergers and gravitational dynamical friction. By analogy with protoplanetary migration, nuclear cluster objects (NCOs) can migrate in the accretion disks that power active galactic nuclei by exchanging angular momentum with disk gas. Here we show that an individual NCO undergoing runaway outward migration comparable to Type III protoplanetary migration can generate an accretion rate corresponding to Seyfert AGN or quasar luminosities. Multiple migrating NCOs in an AGN disk can dominate traditional viscous disk accretion and at large disk radii, ensemble NCO migration and accretion could provide sufficient heating to prevent the gravitational instability from consuming disk gas in star formation. The magnitude and energy of the X-ray soft excess observed at ~0.1-1keV in Seyfert AGN could be explained by a small population of ~10^{2}-10^{3} accreting stellar mass black holes or a few ULXs. NCO migration and accretion in AGN disks are therefore extremely important mechanisms to add to realistic models of AGN disks.

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.

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.

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.

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.

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.

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.

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.

Recently several type Ib supernovae (SNe; with the prototypical SN 2005E) have been shown to have atypical properties. These SNe are faint (absolute peak magnitude of ~ -15) and fast SNe that show unique composition. They are inferred to have low ejecta mass (a few tenths of a solar mass) and to be highly enriched in calcium, but poor in silicon elements and nickel. These SNe were therefore suggested to belong to a new class of calcium-rich faint SNe explosions. Their properties were proposed to be the result of helium detonations that may occur on helium accreting white dwarfs. In this paper we theoretically study the scenario of helium detonations, and focus on the results of detonations in accreted helium layers on low mass carbon-oxygen (CO) cores. We present new results from one dimensional simulations of such explosions, including their light curves and spectra. We find that when the density of the helium layer is low enough the helium detonation produces large amounts of intermediate elements, such as calcium and titanium, together with a large amount of unburnt helium. Our results suggest that the properties of calcium-rich faint SNe could indeed be consistent with the helium-detonation scenario on small CO cores. Above a certain density (larger CO cores) the detonation leaves mainly 56Ni and unburnt helium, and the predicted spectrum will unlikely fit the unique features of this class of SNe. Finally, none of our studied models reproduces the bright, fast evolving light curves of another type of peculiar SNe suggested to originate in helium detonations (SNe 1885A, 1939B and 2002bj).

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.

Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white dwarf star composed mainly of carbon and oxygen, whose mass approaches the Chandrasekhar limit. Such SNe are observed in both young and old stellar environments. Here we report our discovery of the faint type Ib SN 2005E in the halo of the nearby isolated galaxy, NGC 1032. The lack of any trace of recent star formation near the SN location (Fig. 1), and the very low derived ejected mass (~0.3 M_sun), argue strongly against a core-collapse origin for this event. Spectroscopic observations and the derived nucleosynthetic output show that the SN ejecta have high velocities and are dominated by helium-burning products, indicating that SN 2005E was neither a subluminous nor a regular SN Ia (Fig. 2). We have therefore found a new type of stellar explosion, arising from a low-mass, old stellar system, likely involving a binary with a primary white dwarf and a helium-rich secondary. The SN ejecta contain more calcium than observed in any known type of SN and likely additional large amounts of radioactive 44Ti. Such SNe may thus help resolve fundamental physical puzzles, extending from the composition of the primitive solar system and that of the oldest stars, to the Galactic production of positrons.

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.

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.

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.

We present new observations of the nuclear star cluster in the central parsec of the Galaxy with the adaptive optics assisted, integral field spectrograph SINFONI on the ESO/VLT. Our work allows the spectroscopic detection of early and late type stars to m_K >= 16, more than 2 magnitudes deeper than our previous data sets. Our observations result in a total sample of 177 bona fide early-type stars. We find that most of these Wolf Rayet (WR), O- and B- stars reside in two strongly warped disks between 0.8" and 12" from SgrA*, as well as a central compact concentration (the S-star cluster) centered on SgrA*. The later type B stars (m_K>15) in the radial interval between 0.8" and 12" seem to be in a more isotropic distribution outside the disks. The observed dearth of late type stars in the central few arcseconds is puzzling, even when allowing for stellar collisions. The stellar mass function of the disk stars is extremely top heavy with a best fit power law of dN/dm ~ m^(-0.45+/-0.3). Since at least the WR/O-stars were formed in situ in a single star formation event ~6 Myrs ago, this mass function probably reflects the initial mass function (IMF). The mass functions of the S-stars inside 0.8" and of the early-type stars at distances beyond 12" are compatible with a standard Salpeter/Kroupa IMF (best fit power law of dN/dm ~ m^(-2.15+/-0.3)).

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.

Young stars observed in the distant Galactic halo are usually thought to have formed elsewhere, either in the Galactic disk or perhaps the Galactic center, and subsequently ejected at high velocities to their current position. However, some of these stars have apparent lifetimes shorter the required flight time from the Galactic disk/center. We suggest that such stars have evolved in close runaway or hypervelocity binaries. Stellar evolution of such binaries can drive them into mass transfer configurations and even mergers. Such evolution could then rejuvenate them (e.g. blue stragglers) and extend their lifetime after their ejection. The extended lifetimes of such stars could then be reconciled with their flight times to the Galactic halo. We study the possibilities of binary runaway and hypervelocity stars and show that such binaries could have been ejected in triple disruptions and other dynamical interactions with stars or with massive black holes. We show that currently observed "too young" star in the halo could have been ejected from the Galactic disk or the Galactic center and be observable in their current position if they were ejected as binaries. Specifically it is shown that the hypervelocity star HE 0437-5439 could be such a rejuvenated star. Other suggestions for its ejection from the LMC are found to be highly unlikely. Moreover, it is shown that its observed metallicity is most consistent with a Galactic origin and a Galactic center origin can not currently rule out. In addition, we suggest that triple disruptions by the massive black hole in the Galactic center could also capture binaries in close orbits near the MBH, some of which may later evolve to become more massive rejuvenated 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.

We use N-body simulations to study the evolution of the orbital eccentricities of stars deposited near (<0.05 pc) the Milky Way massive black hole (MBH), starting from initial conditions motivated by two competing models for their origin: formation in a disk followed by inward migration; and exchange interactions involving a binary star. The first model predicts modest eccentricities, lower than those observed in the S-star cluster, while the second model predicts higher eccentricities than observed. The N-body simulations include a dense cluster of 10 M_sun stellar black holes (SBHs), expected to accumulate near the MBH by mass segregation. Perturbations from the SBHs tend to randomize the stellar orbits, partially erasing the dynamical signatures of their origin. The eccentricities of the initially highly eccentric stars evolve, in 20 Myr (the S-star lifespan), to a distribution that is consistent at the ~95 % level with the observed eccentricity distribution. In contrast, the eccentricities of the initially more circular orbits fail to evolve to the observed values in 20 Myr, arguing against the disk migration scenario. We find that 20 % - 30 % of the S-stars are tidally disrupted by the MBH over their lifetimes, and that the S-stars are not likely to be ejected as hypervelocity stars outside the central 0.05 pc by close encounters with stellar black holes.

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.

The central parsec around the super-massive black hole in the Galactic Center hosts more than 100 young and massive stars. Outside the central cusp (R~1") the majority of these O and Wolf-Rayet (WR) stars reside in a main clockwise system, plus a second, less prominent disk or streamer system at large angles with respect to the main system. Here we present the results from new observations of the Galactic Center with the AO-assisted near-infrared imager NACO and the integral field spectrograph SINFONI on the ESO/VLT. These include the detection of 27 new reliably measured WR/O stars in the central 12" and improved measurements of 63 previously detected stars, with proper motion uncertainties reduced by a factor of four compared to our earlier work. We develop a detailed statistical analysis of their orbital properties and orientations. Half of the WR/O stars are compatible with being members of a clockwise rotating system. The rotation axis of this system shows a strong transition as a function of the projected distance from SgrA*. The main clockwise system either is either a strongly warped single disk with a thickness of about 10 degrees, or consists of a series of streamers with significant radial variation in their orbital planes. 11 out of 61 clockwise moving stars have an angular separation of more than 30 degrees from the clockwise system. The mean eccentricity of the clockwise system is 0.36+/-0.06. The distribution of the counter-clockwise WR/O star is not isotropic at the 98% confidence level. It is compatible with a coherent structure such as stellar filaments, streams, small clusters or possibly a disk in a dissolving state. The observed disk warp and the steep surface density distribution favor in situ star formation in gaseous accretion disks as the origin of the young stars.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.

The study of the formation of molecular hydrogen on low temperature surfaces is of interest both because it allows to explore elementary steps in the heterogeneous catalysis of a simple molecule and because of the applications in astrochemistry. Here we report results of experiments of molecular hydrogen formation on amorphous silicate surfaces using temperature-programmed desorption (TPD). In these experiments beams of H and D atoms are irradiated on the surface of an amorphous silicate sample. The desorption rate of HD molecules is monitored using a mass spectrometer during a subsequent TPD run. The results are analyzed using rate equations and the activation energies of the processes leading to molecular hydrogen formation are obtained from the TPD data. We show that a model based on a single isotope provides the correct results for the activation energies for diffusion and desorption of H atoms. These results can thus be used to evaluate the formation rate of H_2 on dust grains under the actual conditions present in interstellar clouds.

In recent years several hypervelocity stars (HVSs) have been observed in the halo of our Galaxy. Such stars are thought to be ejected through dynamical interactions near the massive black hole (MBH) in the Galactic center. Three scenarios have been suggested for their ejection; binary disruption by a MBH, scattering by inspiraling IMBH and scattering by stellar BHs close to MBH. These scenarios involve different stellar populations in the Galactic center. Here we use observations of the Galactic center stellar population together with dynamical and evolutionary arguments to obtain strong constraints on the nature and origin of HVSs. We show that the IMBH inspiral scenario requires too many (O(10^3)) main sequence B stars to exist close to the MBH (<0.01 pc) at the time of inspiral, where current observations show O(10) such stars. Scattering by SBHs also require too many B stars to be observed in the GC, but it may contribute a small fraction of the currently observed HVSs. The binary disruption scenario is still consistent with current observations. In addition it is shown that recently suggested signatures for HVSs origin such as hypervelocity binaries and slow rotating HVSs are much weaker than suggested and require too large statistics. In addition, we show that due to the conditions close to the MBH most binary star systems are not expected to survive for long in this region. Consequently, unique stellar populations that require long evolution in binaries are not expected to be ejected as HVSs in the BHs scattering mechanisms (this may also be related to to the recently observed asymmetry in the velocity distribution of HVSs).

Recent observations of the Galactic center revealed a nuclear disk of young OB stars near the massive black hole (MBH), in addition to many similar outlying stars with higher eccentricities and/or high inclinations relative to the disk (some of them possibly belonging to a second disk). In addition, observations show the existence of young B stars (the 'S-cluster') in an isotropic distribution in the close vicinity of the MBH ($<0.04$ pc). We use extended N-body simulations to probe the dynamical evolution of these two populations. We show that the stellar disk could have evolved to its currently observed state from a thin disk of stars formed in a gaseous disk, and that the dominant component in its evolution is the interaction with stars in the cusp around the MBH. We also show that the currently observed distribution of the S-stars could be consistent with a capture origin through 3-body binary-MBH interactions. In this scenario the stars are captured at highly eccentric orbits, but scattering by stellar black holes could change their eccentricity distribution to be consistent with current observations.

We compute the Milky Way potential in different cold dark matter (CDM) based models, and compare these with the modified Newtonian dynamics (MOND) framework. We calculate the axis ratio of the potential in various models, and find that isopotentials are less spherical in MOND than in CDM potentials. As an application of these models, we predict the escape velocity as a function of the position in the Galaxy. This could be useful in comparing with future data from planned or already-underway kinematic surveys (RAVE, SDSS, SEGUE, SIM, GAIA or the hypervelocity stars survey). In addition, the predicted escape velocity is compared with the recently measured high proper motion velocity of the Large Magellanic Cloud (LMC). To bind the LMC to the Galaxy in a MOND model, while still being compatible with the RAVE-measured local escape speed at the Sun's position, we show that an external field modulus of less than $0.03 a_0$ is needed.

We show that dynamical relaxation in the aftermath of a galactic merger and the ensuing formation and decay of a binary massive black hole (MBH), are dominated by massive perturbers (MPs) such as giant molecular clouds or clusters. MPs accelerate relaxation by orders of magnitude relative to 2-body stellar relaxation alone, and efficiently scatter stars into the binary MBH's orbit. The 3-body star-binary MBH interactions shrink the binary MBH to the point where energy losses from the emission of gravitational waves (GW) lead to rapid coalescence. We model this process based on observed and simulated MP distributions and take into account the decreased efficiency of the star-binary MBH interaction due to acceleration in the galactic potential. We show that mergers of gas-rich galactic nuclei lead to binary MBH coalescence well within the Hubble time. Moreover, lower-mass binary MBHs (<10^8 Msun) require only a few percent of the typical gas mass in a post-merger nucleus to coalesce in a Hubble time. The fate of a binary MBH in a gas poor galactic merger is less certain, although massive stellar structures (e.g. clusters, stellar rings) could likewise lead to efficient coalescence. These coalescence events are observable by their strong GW emission. MPs thus increase the cosmic rate of such GW events, lead to a higher mass deficit in the merged galactic core and suppress the formation of triple MBH systems and the resulting ejection of MBHs into intergalactic space.

Recent observations of the Galactic center revealed a nuclear disk of young OB stars, in addition to many similar outlying stars with higher eccentricities and/or high inclinations relative to the disk (some of them possibly belonging to a second disk). Binaries in such nuclear disks, if they exist in non-negligible fractions, could have a major role in the evolution of the disks through binary heating of this stellar system. We suggest that interactions with/in binaries may explain some (or all) of the observed outlying young stars in the Galactic center. Such stars could have been formed in a disk, and later on kicked out from it through binary related interactions, similar to ejection of high velocity runaway OB stars in young clusters throughout the galaxy.

The temperatures of interstellar dust grains are analyzed using stochastic simulations, taking into account the grain composition and size and the discreteness of the photon flux. [...] The distribution of grain temperatures is calculated for a broad range of grain sizes and for different intensities of the interstellar radiation field, relevant to diffuse clouds and to PDRs. The dependence of the average grain temperature on its size is shown for different irradiation intensities. It is found that the average temperatures of grains with radii smaller than about 0.02 $\mu$m are reduced due to the fluctuations. The average temperatures of grains of radii larger than about 0.35 $\mu$m are also slightly reduced due to their more efficient emission of infrared radiation, particularly when exposed to high irradiation intensities. The average temperatures <T> of silicate and carbonaceous grains are found to depend on the radiation field intensity X_MMP according to <T>~X_MMP^gamma, where the exponent gamma depends on the grain size and composition. This fitting formula is expected to be useful in simulations of interstellar processes, and can replace commonly used approximations which do not account for the grain temperature fluctuations and for the detailed properties of interstellar dust particles. The implications of the results on molecular hydrogen formation are also discussed. It is concluded that grain-temperature fluctuations tend to reduce the formation rate of molecular hydrogen, and cannot account for the observations of H_2 in photon dominated regions, even in the presence of chemisorption sites.

Experimental results on the formation of molecular hydrogen on amorphous silicate surfaces are presented 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 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 in diffuse interstellar clouds.

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.

We investigate the role of massive perturbers, such as giant molecular clouds or stellar clusters, in supplying low-angular momentum stars that pass very close to the central massive black hole (MBH) or fall into it. We show that massive pe rturbers can play an important role in supplying both binaries and single stars to the vicinity of the MBH. We discuss possible implications for the ejection of high velocity stars; for the capture of stars on tight orbits around the MBH; for the emission of gravitational waves from low-eccentricity inspiraling stars; and for the origin of the young main sequence B stars observed very near the Gal actic MBH. Massive perturbers may also enhance the the growth rate of MBHs, and may accelerate the dynamical orbital decay of coalescing binary MBHs.

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.

Recent laboratory experiments on interstellar dust analogues have shown that H_2 formation on dust grain surfaces is efficient in a range of grain temperatures below 20 K. These results indicate that surface processes may account for the observed H_2 abundance in cold diffuse and dense clouds. However, high abundances of H_2 have also been observed in warmer clouds, including photon-dominated regions (PDRs), where grain temperatures may reach 50 K, making the surface processes extremely inefficient. It was suggested that this apparent discrepancy can be resolved by chemisorption sites. However, recent experiments indicate that chemisorption processes may not be efficient at PDR temperatures. Here we consider the effect of grain porosity on H_2 formation. It is found that porosity extends the efficiency of the recombination process to higher temperatures. This is because H atoms that desorb from the internal surfaces of the pores may re-adsorb many times and thus stay longer on the surface. However, this porosity-driven extension may enable efficient H_2 formation in PDRs only if porosity also contributes to significant cooling of the grains, compared to non-porous grains.

The formation of H2 and HD molecules on interstellar dust grains is studied using rate equation and master equation models. Rate equations are used in the analysis of laboratory experiments which examine the formation of molecular hydrogen on astrophysically relevant surfaces. However, under interstellar conditions, rate equations are not suitable for the calculation of reaction rates on dust-grain surfaces. Due to the low flux and the sub-micron size of the grains, the populations of H and D atoms on a single grain are likely to be small. In this case the reaction rates are dominated by fluctuations and should be calculated using stochastic methods. The rate of molecular hydrogen formation in interstellar clouds is evaluated using the master equation, taking into account the distribution of grain sizes.

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.

The results of experiments on the formation of molecular hydrogen on low density and high density amorphous ice surfaces are analyzed using a rate equation model. The activation energy barriers for the relevant diffusion and desorption processes are obtained. The more porous morphology of the low density ice gives rise to a broader spectrum of energy barriers compared to the high density ice. Inserting these parameters into the rate equation model under steady state conditions we evaluate the production rate of molecular hydrogen on ice-coated interstellar dust grains.