Circularizing Planet Nine through dynamical friction with an extended, cold planetesimal belt

Unexpected clustering in the orbital elements of minor bodies beyond the Kuiper belt has led to speculations that our solar system actually hosts nine planets, the eight established plus a hypothetical "Planet Nine". Several recent studies have shown that a planet with a mass of about 10 Earth masses on a distant eccentric orbit with perihelion far beyond the Kuiper belt could create and maintain this clustering. The evolutionary path resulting in an orbit such as the one suggested for Planet Nine is nevertheless not easily explained. Here we investigate whether a planet scattered away from the giant-planet region could be lifted to an orbit similar to the one suggested for Planet Nine through dynamical friction with a cold, distant planetesimal belt. Recent simulations of planetesimal formation via the streaming instability suggest that planetesimals can readily form beyond 100au. We explore this circularisation by dynamical friction with a set of numerical simulations. We find that a planet that is scattered from the region close to Neptune onto an eccentric orbit has a 20-30% chance of obtaining an orbit similar to that of Planet Nine after 4.6Gyr. Our simulations also result in strong or partial clustering of the planetesimals; however, whether or not this clustering is observable depends on the location of the inner edge of the planetesimal belt. If the inner edge is located at 200au the degree of clustering amongst observable objects is significant.Comment: Accepted to MNRA

Numerical study of resonant spin relaxation in quasi-1D channels

Recent experiments demonstrate that a ballistic version of spin resonance, mediated by spin-orbit interaction, can be induced in narrow channels of a high-mobility GaAs two-dimensional electron gas by matching the spin precession frequency with the frequency of bouncing trajectories in the channel. Contrary to the typical suppression of Dyakonov-Perel' spin relaxation in confined geometries, the spin relaxation rate increases by orders of magnitude on resonance. Here, we present Monte Carlo simulations of this effect to explore the roles of varying degrees of disorder and strength of spin-orbit interaction. These simulations help to extract quantitative spin-orbit parameters from experimental measurements of ballistic spin resonance, and may guide the development of future spintronic devices

Simulations of the tidal interaction and mass transfer of a star in an eccentric orbit around an intermediate-mass black hole: the case of HLX-1

The X-ray source HLX-1 near the spiral galaxy ESO 243-49 is currently the best intermediate-mass black hole candidate. It has a peak bolometric luminosity of $10^{42}$ erg s$^{-1}$, which implies a mass inflow rate of $\sim10^{-4}$ MSun yr$^{-1}$, but the origin of this mass is unknown. It has been proposed that there is a star on an eccentric orbit around the black hole which transfers mass at pericentre. To investigate the orbital evolution of this system, we perform stellar evolution simulations using mesa and SPH simulations of a stellar orbit around an intermediate-mass black hole using fi. We run and couple these simulations using the amuse framework. We find that mass is lost through both the first and second Lagrange points and that there is a delay of up to 10 days between the pericentre passage and the peak mass loss event. The orbital evolution timescales we find in our simulations are larger than what is predicted by analytical models, but these models fall within the errors of our results. Despite the fast orbital evolution, we are unable to reproduce the observed change in outburst period. We conclude that the change in the stellar orbit with the system parameters investigated here is unable to account for all observed features of HLX-1.Comment: accepted for publication in mnra

Simulating the Phases of the Moon Shortly After Its Formation

The leading theory for the origin of the Moon is the giant impact hypothesis, in which the Moon was formed out of the debris left over from the collision of a Mars-sized body with the Earth. Soon after its formation, the orbit of the Moon may have been very different than it is today. We have simulated the phases of the Moon in a model for its formation wherein the Moon develops a highly elliptical orbit with its major axis tangential to the Earth's orbit. This note describes these simulations and their pedagogical value.Comment: 4 pages, 2 figure

The Pairing of Spin-orbit Coupled Fermi Gas in Optical Lattice

We investigate Rashba spin-orbit coupled Fermi gases in square optical lattice by using the determinant quantum Monte Carlo (DQMC) simulations which is free of the sign-problem. We show that the Berezinskii-Kosterlitz-Thoules phase transition temperature is firstly enhanced and then suppressed by spin-orbit coupling in the strong attraction region. In the intermediate attraction region, spin-orbit coupling always suppresses the transition temperature. We also show that the spin susceptibility becomes anisotropic and retains finite values at zero temperature.Comment: 10 pages, 9 figure

A unified framework for the orbital structure of bars and triaxial ellipsoids

We examine a large random sample of orbits in two self-consistent simulations of N-body bars. Orbits in these bars are classified both visually and with a new automated orbit classification method based on frequency analysis. The well-known prograde x1 orbit family originates from the same parent orbit as the box orbits in stationary and rotating triaxial ellipsoids. However, only a small fraction of bar orbits (~4%) have predominately prograde motion like their periodic parent orbit. Most bar orbits arising from the x1 orbit have little net angular momentum in the bar frame, making them equivalent to box orbits in rotating triaxial potentials. In these simulations a small fraction of bar orbits (~7%) are long-axis tubes that behave exactly like those in triaxial ellipsoids: they are tipped about the intermediate axis owing to the Coriolis force, with the sense of tipping determined by the sign of their angular momentum about the long axis. No orbits parented by prograde periodic x2 orbits are found in the pure bar model, but a tiny population (~2%) of short-axis tube orbits parented by retrograde x4 orbits are found. When a central point mass representing a supermassive black hole (SMBH) is grown adiabatically at the center of the bar, those orbits that lie in the immediate vicinity of the SMBH are transformed into precessing Keplerian orbits that belong to the same major families (short-axis tubes, long-axis tubes and boxes) occupying the bar at larger radii. During the growth of an SMBH, the inflow of mass and outward transport of angular momentum transform some x1 and long-axis tube orbits into prograde short-axis tubes. This study has important implications for future attempts to constrain the masses of SMBHs in barred galaxies using orbit-based methods like the Schwarzschild orbit superposition scheme and for understanding the observed features in barred galaxies

Tunable superconducting critical temperature in ballistic hybrid structures with strong spin-orbit coupling

We present a theoretical description and numerical simulations of the superconducting transition in hybrid structures including strong spin-orbit interactions. The spin-orbit coupling is taken to be of Rashba type for concreteness, and we allow for an arbitrary magnitude of the spin-orbit strength as well as an arbitrary thickness of the spin-orbit coupled layer. This allows us to make contact with the experimentally relevant case of enhanced interfacial spin-orbit coupling via atomically thin heavy metal layers. We consider both interfacial spin-orbit coupling induced by inversion asymmetry in an S/F-junction, as well as in-plane spin-orbit coupling in the ferromagnetic region of an S/F/S- and an S/F-structure. Both the pair amplitudes, local density of states and critical temperature show dependency on the Rashba strength and, importantly, the orientation of the exchange field. In general, spin-orbit coupling increases the critical temperature of a proximity system where a magnetic field is present, and enhances the superconducting gap in the density of states. We perform a theoretical derivation which explains these results by the appearance of long-ranged singlet correlations. Our results suggest that $T_c$ in ballistic spin-orbit coupled superconducting structures may be tuned by using only a single ferromagnetic layer.Comment: 14 pages, 14 figures. Published in PR

An Analytic Model of Angular Momentum Transport by Gravitational Torques: From Galaxies to Massive Black Holes

We present analytic calculations of angular momentum transport and gas inflow in galaxies, from scales of ~kpc to deep in the potential of a central black hole (BH). We compare these analytic calculations to numerical simulations and use them to develop a sub-grid model of BH growth that can be incorporated into semi-analytic models or cosmological simulations. Both analytic calculations and simulations argue that the strongest torque on gas arises when non-axisymmetric perturbations to the stellar gravitational potential produces orbit crossings and shocks in the gas. This is true both at large radii ~0.01-1 kpc, where bar-like modes dominate the non-axisymmetric potential, and at smaller radii <10 pc, where a lopsided/eccentric disk dominates. The traditional orbit crossing criterion is not always adequate to predict the locations of, and inflow due to, shocks in gas+stellar disks with finite sound speeds. We derive a modified criterion that predicts the presence of shocks in stellar dominated systems even absent formal orbit crossing. We then derive analytic expressions for the loss of angular momentum and the resulting gas inflow rates in the presence of such shocks. We test our analytic predictions using hydrodynamic simulations at a range of galactic scales, and show that they successfully predict the mass inflow rates and quasi-steady gas surface densities with small scatter (0.3 dex). We use our analytic results to construct a new estimate of the BH accretion rate given galaxy properties at larger radii. This captures the key scalings in the numerical simulations. Alternate estimates such as the local viscous accretion rate or the spherical Bondi rate fail systematically to reproduce the simulations.Comment: 23 pages, 10 figures, minor revisions to match version accepted to MNRA