Posts Tagged mass ratio

Recent Postings from mass ratio

A survey of dual active galactic nuclei in simulations of galaxy mergers: frequency and properties

We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger remnant has formed, especially at high luminosities. For typical minor and major mergers, dual activity lasts ~20-70 and ~100-160 Myr, respectively. We also explore the effects of X-ray obscuration from gas, finding that the dual-AGN time decreases at most by a factor of ~2, and of contamination from star formation. Using projected separations and velocity differences rather than three-dimensional quantities can decrease the dual-AGN time-scales by up to ~4, and we apply filters which mimic current observational-resolution limitations. In agreement with observations, we find that, for a sample of major mergers hosting at least one AGN, ~20 per cent of them should harbour dual AGN. We quantify the effects of merger mass ratio (0.1 to 1), geometry (coplanar, prograde, retrograde, and inclined), disc gas fraction, and BH properties, finding that the mass ratio is the most important factor, with the difference between minor and major mergers varying between factors of a few to orders of magnitude, depending on the luminosity and filter used. We also find that a deep imaging survey does not need very high angular resolution, whereas a shallow survey requires it.

A low-mass-ratio and deep contact binary as the progenitor of the merger V1309 Sco

Nova Sco 2008 (=V1309 Sco) is an example of a V838 Mon type eruption rather than a typical classical nova. This enigmatic object was recently shown to have resulted from the merger of two stars in a contact binary. It is the first stellar merger that was identified to be undergoing a common envelope transient. To understand the properties of its binary progenitor, the pre-outburst light curves were analyzed by using the W-D method. The photometric solution of the 2002 light curve shows that it is a deep contact binary (f = 89.5(+-40.5)%) with a mass ratio of 0.094. The asymmetry of the light curve is explained by the presence of a dark spot on the more massive component. The extremely high fill-out factor suggests that the merging of the contact binary is driven by dynamical mass loss from the outer Lagrange point. However, the analysis of the 2004 light curve indicates that no solutions were obtained even at an extremely low mass ratio of q = 0.03. This suggests that the common convective envelope of the binary system disappeared and the secondary component spiraled into the envelope of the primary in 2004. Finally, the ejection of the envelope of the primary produced the outburst.

Fundamental frequencies and resonances from eccentric and precessing binary black hole inspirals

Binary black holes which are both eccentric and undergo precession remain unexplored in numerical simulations. We present simulations of such systems which cover about 50 orbits at comparatively high mass ratios 5 and 7. The configurations correspond to the generic motion of a nonspinning body in a Kerr spacetime, and are chosen to study the transition from finite mass-ratio inspirals to point particle motion in Kerr. We develop techniques to extract analogs of the three fundamental frequencies of Kerr geodesics, compare our frequencies to those of Kerr, and show that the differences are consistent with self-force corrections entering at first order in mass ratio. This analysis also locates orbital resonances where the ratios of our frequencies take rational values. At the considered mass ratios, the binaries pass through resonances in one to two resonant cycles, and we find no discernible effects on the orbital evolution. We also compute the decay of eccentricity during the inspiral and find good agreement with the leading order post-Newtonian prediction.

Hierarchical data-driven approach to fitting numerical relativity data for non-precessing binary black holes, with an application to final spin and radiated energy

Numerical relativity is an essential tool in studying the coalescence of binary black holes (BBHs). It is still computationally prohibitive to cover the BBH parameter space exhaustively, making phenomenological fitting formulas for BBH waveforms and final-state properties important for practical applications. We describe a general hierarchical bottom-up fitting methodology to design and calibrate fits to numerical relativity simulations for the three-dimensional parameter space of quasi-circular non-precessing merging BBHs, spanned by mass ratio and by the individual spin components orthogonal to the orbital plane. Particular attention is paid to incorporating the extreme-mass-ratio limit and to the subdominant unequal-spin effects. As an illustration of the method, we provide two applications, to the final spin and final mass (or equivalently: radiated energy) of the remnant black hole. We obtain results broadly consistent with previously published fits, but with improvements in the approach to extremal limits and for unequal-spin configurations. We also discuss the importance of data quality studies when combining simulations from diverse sources, how detailed error budgets will be necessary for further improvements of these already highly-accurate fits, and how this first detailed study of unequal-spin effects helps in choosing the most informative parameters for future NR runs.

The dust-to-stellar mass ratio as a valuable tool to probe the evolution of local and distant star forming galaxies

The survival of dust grains in galaxies depends on various processes. Dust can be produced in stars, it can grow in the interstellar medium and be destroyed by astration and interstellar shocks. In this paper, we assemble a few data samples of local and distant star-forming galaxies to analyse various dust-related quantities in low and high redshift galaxies, to study how the relations linking the dust mass to the stellar mass and star formation rate evolve with redshift. We interpret the available data by means of chemical evolution models for discs and proto-spheroid (PSPH) starburst galaxies. In particular, we focus on the dust-to-stellar mass (DTS) ratio, as this quantity represents a true measure of how much dust per unit stellar mass survives the various destruction processes in galaxies and is observable. The theoretical models outline the strong dependence of this quantity on the underlying star formation history. Spiral galaxies are characterised by a nearly constant DTS as a function of the stellar mass and cosmic time, whereas PSPHs present an early steep increase of the DTS, which stops at a maximal value and decreases in the latest stages. In their late starburst phase, these models show a decrease of the DTS with their mass, which allows us to explain the observed anti-correlation between the DTS and the stellar mass. The observed redshift evolution of the DTS ratio shows an increase from z~0 to z~1, followed by a roughly constant behaviour at 1<z<2.5. Our models indicate a steep decrease of the global DTS at early times, which implies an expected decrease of the DTS at larger redshift.

Superoutburst of WZ Sge-type Dwarf Nova Below the Period Minimum: ASASSN-15po

We report on a superoutburst of a WZ Sge-type dwarf nova (DN), ASASSN-15po. The light curve showed the main superoutburst and multiple rebrightenings. In this outburst, we observed early superhumps and growing (stage A) superhumps with periods of 0.050454(2) and 0.051809(13) d, respectively. We estimated that the mass ratio of secondary to primary ($q$) is 0.0699(8) by using $P_{\rm orb}$ and a superhump period $P_{\rm SH}$ of stage A. ASASSN-15po [$P_{\rm orb} \sim$ 72.6 min] is the first DN with the orbital period between 67--76 min. Although the theoretical predicted period minimum $P_{\rm min}$ of hydrogen-rich cataclysmic variables (CVs) is about 65--70 min, the observational cut-off of the orbital period distribution at 80 min implies that the period minimum is about 82 min, and the value is widely accepted. We suggest the following four possibilities: the object is (1) a theoretical period minimum object (2) a binary with a evolved secondary (3) a binary with a metal-poor (Popullation II) seconday (4) a binary which was born with a brown-dwarf donor below the period minimum.

A possible formation channel for blue hook stars in globular cluster - II. Effects of metallicity, mass ratio, tidal enhancement efficiency and helium abundance

Employing tidally enhanced stellar wind, we studied in binaries the effects of metallicity, mass ratio of primary to secondary, tidal enhancement efficiency and helium abundance on the formation of blue hook (BHk) stars in globular clusters (GCs). A total of 28 sets of binary models combined with different input parameters are studied. For each set of binary model, we presented a range of initial orbital periods that is needed to produce BHk stars in binaries. All the binary models could produce BHk stars within different range of initial orbital periods. We also compared our results with the observation in the Teff-logg diagram of GC NGC 2808 and {\omega} Cen. Most of the BHk stars in these two GCs locate well in the region predicted by our theoretical models, especially when C/N-enhanced model atmospheres are considered. We found that mass ratio of primary to secondary and tidal enhancement efficiency have little effects on the formation of BHk stars in binaries, while metallicity and helium abundance would play important roles, especially for helium abundance. Specifically, with helium abundance increasing in binary models, the space range of initial orbital periods needed to produce BHk stars becomes obviously wider, regardless of other input parameters adopted. Our results were discussed with recent observations and other theoretical models.

RZ Leonis Minoris Bridging between ER Ursae Majoris-Type Dwarf Nova and Novalike System

We observed RZ LMi, which is renowned for the extremely (~19d) short supercycle and is a member of a small, unusual class of cataclysmic variables called ER UMa-type dwarf novae, in 2013 and 2016. In 2016, the supercycles of this object substantially lengthened in comparison to the previous measurements to 35, 32, 60d for three consecutive superoutbursts. We consider that the object virtually experienced a transition to the novalike state (permanent superhumper). This observed behavior extremely well reproduced the prediction of the thermal-tidal instability model. We detected a precursor in the 2016 superoutburst and detected growing (stage A) superhumps with a mean period of 0.0602(1)d in 2016 and in 2013. Combined with the period of superhumps immediately after the superoutburst, the mass ratio is not as small as in WZ Sge-type dwarf novae, having orbital periods similar to RZ LMi. By using least absolute shrinkage and selection operator (Lasso) two-dimensional power spectra, we detected possible negative superhumps with a period of 0.05710(1)d. We estimated the orbital period of 0.05792d, which suggests a mass ratio of 0.105(5). This relatively large mass ratio is even above ordinary SU UMa-type dwarf novae, and it is also possible that the exceptionally high mass-transfer rate in RZ LMi may be a result of a stripped core evolved secondary which are evolving toward an AM CVn-type object.

WR 148: Identifying the companion of an extreme runaway massive binary

WR 148 (HD 197406) is an extreme runaway system considered to be a potential candidate for a short-period (4.3173 d) rare WR + compact object binary. Provided with new high resolution, high signal-to-noise spectra from the Keck observatory, we determine the orbital parameters for both the primary WR and the secondary, yielding respective projected orbital velocity amplitudes of $88.1\pm3.8$ km s$^{-1}$ and $79.2\pm3.1$ km s$^{-1}$ and implying a mass ratio of $1.1\pm0.1$. We then apply the shift-and-add technique to disentangle the spectra and obtain spectra compatible with a WN7ha and an O4-6 star. Considering an orbital inclination of $\sim67^\circ$, derived from previous polarimetry observations, the system's total mass would be a mere 2-3 M$_{\odot}$ , an unprecedented result for a putative massive binary system. However, a system comprising a 37 M$_{\odot}$ secondary (typical mass of an O5V star) and a 33 M$_{\odot}$ primary (given the mass ratio) would infer an inclination of $\sim18^\circ$. We therefore reconsider the previous methods of deriving the orbital inclination based on time-dependent polarimetry and photometry. While the polarimetric results are inconclusive requiring better data, the photometric results favour low inclinations. Finally, we compute WR 148's space velocity and retrace the runaway's trajectory back to the Galactic plane (GP). With an ejection velocity of $198\pm27$ km s$^{-1}$ and a travel time of $4.7\pm0.8$ Myr to reach its current location, WR 148 was most likely ejected via dynamical interactions in a young cluster.

Deriving analytic solutions for compact binary inspirals without recourse to adiabatic approximations [Cross-Listing]

We utilize the dynamical renormalization group formalism to calculate the real space trajectory of a compact binary inspiral for long times via a systematic resummation of secularly growing terms. This method generates closed form solutions without orbit averaging, and the accuracy can be systematically improved. The expansion parameter is $v^5 \nu \Omega(t-t_0)$ where $t_0$ is the initial time, $t$ is the time elapsed, and $\Omega$ and $v$ are the angular orbital frequency and initial speed, respectively, and $\nu$ is the binary's symmetric mass ratio. We demonstrate how to apply the renormalization group method to resum solutions beyond leading order in two ways. First, we calculate the second order corrections of the leading radiation reaction force, which involves highly non-trivial checks of the formalism (i.e. its renormalizability). Second, we show how to systematically include post-Newtonian corrections to the radiation reaction force. By avoiding orbit averaging we gain predictive power and eliminate ambiguities in the initial conditions. Finally, we discuss how this methodology can be used to find analytic solutions to the spin equations of motion that are valid over long times.

Deriving analytic solutions for compact binary inspirals without recourse to adiabatic approximations

We utilize the dynamical renormalization group formalism to calculate the real space trajectory of a compact binary inspiral for long times via a systematic resummation of secularly growing terms. This method generates closed form solutions without orbit averaging, and the accuracy can be systematically improved. The expansion parameter is $v^5 \nu \Omega(t-t_0)$ where $t_0$ is the initial time, $t$ is the time elapsed, and $\Omega$ and $v$ are the angular orbital frequency and initial speed, respectively, and $\nu$ is the binary's symmetric mass ratio. We demonstrate how to apply the renormalization group method to resum solutions beyond leading order in two ways. First, we calculate the second order corrections of the leading radiation reaction force, which involves highly non-trivial checks of the formalism (i.e. its renormalizability). Second, we show how to systematically include post-Newtonian corrections to the radiation reaction force. By avoiding orbit averaging we gain predictive power and eliminate ambiguities in the initial conditions. Finally, we discuss how this methodology can be used to find analytic solutions to the spin equations of motion that are valid over long times.

Deriving analytic solutions for compact binary inspirals without recourse to adiabatic approximations [Cross-Listing]

We utilize the dynamical renormalization group formalism to calculate the real space trajectory of a compact binary inspiral for long times via a systematic resummation of secularly growing terms. This method generates closed form solutions without orbit averaging, and the accuracy can be systematically improved. The expansion parameter is $v^5 \nu \Omega(t-t_0)$ where $t_0$ is the initial time, $t$ is the time elapsed, and $\Omega$ and $v$ are the angular orbital frequency and initial speed, respectively, and $\nu$ is the binary's symmetric mass ratio. We demonstrate how to apply the renormalization group method to resum solutions beyond leading order in two ways. First, we calculate the second order corrections of the leading radiation reaction force, which involves highly non-trivial checks of the formalism (i.e. its renormalizability). Second, we show how to systematically include post-Newtonian corrections to the radiation reaction force. By avoiding orbit averaging we gain predictive power and eliminate ambiguities in the initial conditions. Finally, we discuss how this methodology can be used to find analytic solutions to the spin equations of motion that are valid over long times.

Relic abundance of MeV millicharged particles [Replacement]

The relic abundance of light millicharged particles (MCP) with the electric charge $e' = 5\cdot 10^{-5} e$ and with the mass slightly below or above the electron mass is calculated. The abundance depends on the mass ratio $\eta=m_X/m_e$ and for $\eta<1$ can be high enough to allow MCP to be the cosmological dark matter or to make a noticeable contribution to it. On the other hand, for $\eta \gtrsim 1$ the cosmological energy density of MCPs can be quite low, $\Omega_X h_0^2 \simeq 0.02$ for scalar MCPs, and $\Omega_X h_0^2 \simeq 0.001$ for spin 1/2 fermions. But even the lowest value of $\Omega_X h_0^2$ is in tension with several existing limits on the MCP abundances and parameters. However, these limits have been derived under some natural or reasonable assumptions on the properties of MCPs. If these assumptions are relaxed, a patch in the mass-charge plot of MCPs may appear, permitting them to be dark matter particles.

Relic abundance of MeV millicharged particles

The relic abundance of light millicharged particles (MCP) with the electric charge $e' = 5\cdot 10^{-5} e$ and with the mass slightly below or above the electron mass is calculated. The abundance depends on the mass ratio $\eta=m_X/m_e$ and for $\eta<1$ can be high enough to allow MCP to be the cosmological dark matter or to make a noticeable contribution to it. On the other hand, for $\eta \gtrsim 1$ the cosmological energy density of MCPs can be quite low, $\Omega_X h_0^2 \simeq 0.02$ for scalar MCPs, and $\Omega_X h_0^2 \simeq 0.001$ for spin 1/2 fermions. But even the lowest value of $\Omega_X h_0^2$ is in tension with several existing limits on the MCP abundances and parameters. However, these limits have been derived under some natural or reasonable assumptions on the properties of MCPs. If these assumptions are relaxed, a patch in the mass-charge plot of MCPs may appear, permitting them to be dark matter particles.

The First Circumbinary Planet Found by Microlensing: OGLE-2007-BLG-349L(AB)c [Replacement]

We present the analysis of the first circumbinary planet microlensing event, OGLE-2007-BLG-349. This event has a strong planetary signal that is best fit with a mass ratio of $q \approx 3.4\times10^{-4}$, but there is an additional signal due to an additional lens mass, either another planet or another star. We find acceptable light curve fits with two classes of models: 2-planet models (with a single host star) and circumbinary planet models. The light curve also reveals a significant microlensing parallax effect, which constrains the mass of the lens system to be $M_L \approx 0.7 M_\odot$. Hubble Space Telescope images resolve the lens and source stars from their neighbors and indicate excess flux due to the star(s) in the lens system. This is consistent with the predicted flux from the circumbinary models, where the lens mass is shared between two stars, but there is not enough flux to be consistent with the 2-planet, 1-star models. So, only the circumbinary models are consistent with the HST data. They indicate a planet of mass $m_c = 80\pm 13\,M_\oplus$, orbiting a pair of M-dwarfs with masses of $M_A = 0.41\pm 0.07 M_\odot$ and $M_B = 0.30\pm 0.07 M_\oplus$, which makes this the lowest mass circumbinary planet system known. The ratio of the separation between the planet and the center-of-mass to the separations of the two stars is $\sim 40$, so unlike most of the circumbinary planets found by Kepler, the planet does not orbit near the stability limit.

The First Circumbinary Planet Found by Microlensing: OGLE-2007-BLG-349L(AB)c

We present the analysis of the first circumbinary planet microlensing event, OGLE-2007-BLG-349. This event has a strong planetary signal that is best fit with a mass ratio of $q \approx 3.4\times 10^{-4}$, but there is an additional signal due to an additional lens mass, either another planet or another star. We find acceptable light curve fits with two classes of models: 2-planet models (with a single host star) and circumbinary planet models. The light curve also reveals a significant microlensing parallax effect, which constraints the mass of the lens system to be $M_L \approx 0.7 M_\odot$. Hubble Space Telescope images resolve the lens and source stars from their neighbors, and indicate excess flux due to the star(s) in the lens system. This is consistent with the predicted flux from the circumbinary models, where the lens mass is shared between two stars, but there is not enough flux to be consistent with the 2-planet, 1-star models. So, only the circumbinary models are consistent with the HST data. They indicate a planet of mass $m_c = 80\pm 13 M_\oplus$, orbiting a pair of M-dwarfs with masses of $M_A = 0.41\pm 0.07 M_\odot$ and $M_B = 0.30\pm 0.07 M_\odot$, which makes this the lowest mass circumbinary planet system known. The ratio of the planet:center-of-mass separation to the separations of the two stars is ~40, so unlike most of the circumbinary planets found by Kepler, the planet does not orbit near the stability limit.

Rotating systems, universal features in dragging and anti-dragging effects, and bounds onto angular momentum

We consider stationary, axially symmetric toroids rotating around spinless black holes, assuming the general-relativistic Keplerian rotation law, in the first post-Newtonian approximation. Numerical investigation shows that the angular momentum accumulates almost exclusively within toroids. It appears that various types of dragging (anti-dragging) effects are positively correlated with the ratio $M_\mathrm{D}/m$ ($M_\mathrm{D}$ is the mass of a toroid and $m$ is the mass of the black hole) - moreover, their maxima are proportional to $M_\mathrm{D}/m$. The horizontal sizes of investigated toroids range from c. 50 to c. 450 of Schwarzschild radii $R_\mathrm{S}$ of the central black hole; their mass $M_\mathrm{D} \in (10^{-4}m, 40m)$ and the radial size of the system is c. 500 $R_\mathrm{S}$. We found that the relative strength of various dragging (anti-dragging) effects does not change with the mass ratio, but it depends on the size of toroids. Several isoperimetric inequalities involving angular momentum are shown to hold true.

Why are pulsar planets rare?

Pulsar timing observations have revealed planets around only a few pulsars. We suggest that the rarity of these planets is due mainly to two effects. First, we show that the most likely formation mechanism requires the destruction of a companion star. Only pulsars with a suitable companion (with an extreme mass ratio) are able to form planets. Second, while a dead zone (a region of low turbulence) in the disk is generally thought to be essential for planet formation, it is most probably rare in disks around pulsars because of the irradiation from the pulsar. The irradiation strongly heats the inner parts of the disk pushing the inner boundary of the dead zone out. We suggest that the rarity of pulsar planets can be explained by the low probability for these two requirements - a very low-mass companion and a dead zone - to be satisfied.

$\gamma$ Doradus Pulsations in the Eclipsing Binary Star KIC 6048106

We present the ${\it Kepler}$ photometry of KIC 6048106 exhibiting O'Connell effect and multiperiodic pulsations. Including a starspot on either of the components, light-curve synthesis indicates that this system is a semi-detached Algol with a mass ratio of 0.211, an orbital inclination of 73.9 deg, and a large temperature difference of 2,534 K. To examine in detail both spot variations and pulsations, we separately analyzed the {\it Kepler} time-series data at the interval of an orbital period by an iterative way. The results reveal that the variable asymmetries of the light maxima can be interpreted as the changes of a magnetic cool spot on the secondary component with time. Multiple frequency analyses were performed in the outside-eclipse light residuals after removal of the binarity effects from the observed {\it Kepler} data. We detected 30 frequencies with signal to noise amplitude ratios larger than 4.0, of which six ($f_2$--$f_6$ and $f_{10}$) can be identified as high-order (17 $\le n \le$ 25) low-degree ($\ell$ = 2) gravity-mode pulsations that were stable during the observing run of 200 d. In contrast, the other frequencies may be harmonic and combination terms. For the six frequencies, the pulsation periods and pulsation constants are in the ranges of 0.352$-$0.506 d and 0.232$-$0.333 d, respectively. These values and the position on the HR diagram demonstrate that the primary star is a $\gamma$ Dor variable. The evolutionary status and the pulsation nature of KIC 6048106 are discussed.

The Post-Starburst Evolution of Tidal Disruption Event Host Galaxies

Tidal Disruption Events (TDEs) favor quiescent host galaxies with strong Balmer absorption lines. Here we study eight hosts of optical/UV-detected TDEs to determine the duration of the recent star formation episode, the time elapsed since it ended, and the fraction of stellar mass produced. Most hosts (6/8) have had short recent starbursts of <200 Myr as opposed to a slower decline in star formation. TDE host galaxies span a wide range of post-starburst ages (60-600 Myr for 6/8 galaxies), indicating that TDEs are not limited to a specific time in their hosts' post-starburst evolution. If the disrupted star was a main sequence star that formed in the burst or before, the post-burst ages provide an independent constraint on its likely mass, excluding O, B and the most massive A stars. If the starburst arose from a galaxy merger, the time elapsed since the starburst began constrains the coalescence timescale and thus limits the merger mass ratio to more equal than 12:1 in most (7/8) TDE hosts. This uncommon ratio, if it also reflects that of the central SMBH binary, disfavors the scenario in which the TDE rate is boosted by the binary but is insensitive to its mass ratio. The fraction of stellar mass created in the burst is 0.5 - 10% for most (7/8) of the TDE hosts, not large enough to explain the increased TDE rate. If more stars are required to boost the TDE rate, the stellar concentration in the core must be more important. TDE host galaxies have stellar masses 10^9.4 - 10^10.3 M$_\odot$, consistent with the SDSS volume-corrected comparison sample and implying central black hole masses of 10^5.5 - 10^7.5 M$_\odot$. Subtracting the absorption line spectra, we uncover hidden emission lines; at least 5 of 8 hosts have ionization sources inconsistent with star formation. These ionization sources may be related to circumnuclear gas, merger shocks, or post-AGB stars.

A numerical investigation of wind accretion in persistent Supergiant X-ray Binaries I - Structure of the flow at the orbital scale

Classical Supergiant X-ray Binaries host a neutron star orbiting a supergiant OB star and display persistent X-ray luminosities of 10$^{35}$ to 10$^{37}$ erg/s. The stellar wind from the massive companion is believed to be the main source of matter accreted by the compact object. With this first paper, we introduce a ballistic model to characterize the structure of the wind at the orbital scale as it accelerates, from the stellar surface to the vicinity of the accretor. Thanks to the parametrization we retained and the numerical pipeline we designed, we can investigate the supersonic flow and the subsequent observables as a function of a reduced set of characteristic numbers and scales. We show that the shape of the permanent flow is entirely determined by the mass ratio, the filling factor, the Eddington factor and the $\alpha$-force multiplier which drives the stellar wind acceleration. Provided scales such as the orbital period are known, we can trace back the observables to evaluate the mass accretion rates, the accretion mechanism, the shearing of the inflow and the stellar parameters. We confront our model to three persistent Supergiant X-ray Binaries (Vela X-1, IGR J18027-2016, XTE J1855-026) and discuss the likelihood of wind-formed accretion discs around the accretors in each case, further investigated in a following paper.

Fractal basins of attraction in the planar circular restricted three-body problem with oblateness and radiation pressure

In this paper we use the planar circular restricted three-body problem where one of the primary bodies is an oblate spheroid or an emitter of radiation in order to determine the basins of attraction associated with the equilibrium points. The evolution of the position of the five Lagrange points is monitored when the values of the mass ratio $\mu$, the oblateness coefficient $A_1$, and the radiation pressure factor $q$ vary in predefined intervals. The regions on the configuration $(x,y)$ plane occupied by the basins of attraction are revealed using the multivariate version of the Newton-Raphson method. The correlations between the basins of convergence of the equilibrium points and the corresponding number of iterations needed in order to obtain the desired accuracy are also illustrated. We conduct a thorough and systematic numerical investigation demonstrating how the dynamical quantities $\mu$, $A_1$, and $q$ influence the basins of attractions. Our results suggest that the mass ratio and the radiation pressure factor are the most influential parameters, while on the other hand the structure of the basins of convergence are much less affected by the oblateness coefficient.

Spin-orbit precession for eccentric black hole binaries at first order in the mass ratio [Replacement]

We consider spin-orbit ("geodetic") precession for a compact binary in strong-field gravity. Specifically, we compute $\psi$, the ratio of the accumulated spin-precession and orbital angles over one radial period, for a spinning compact body orbiting a non-rotating black hole. We show that $\psi$ can be computed for eccentric orbits in both the gravitational self-force and post-Newtonian frameworks, and that the results appear to be consistent. We present a post-Newtonian expansion for $\psi$ at next-to-next-to-leading order, and a Lorenz-gauge gravitational self-force calculation for $\psi$ at first order in the mass ratio. The latter provides new numerical data in the strong-field regime to inform the Effective One-Body model of the gravitational two-body problem. We conclude that $\psi$ complements the Detweiler redshift $z$ as a key invariant quantity characterizing eccentric orbits in the gravitational two-body problem.

Spin-orbit precession for eccentric black hole binaries at first order in the mass ratio

We consider spin-orbit ("geodetic") precession for a compact binary in strong-field gravity. Specifically, we compute $\psi$, the ratio of the accumulated spin-precession and orbital angles over one radial period, for a spinning compact body orbiting a non-rotating black hole. We show that $\psi$ can be computed for eccentric orbits in both the gravitational self-force and post-Newtonian frameworks, and that the results appear to be consistent. We present a post-Newtonian expansion for $\psi$ at next-to-next-to-leading order, and a Lorenz-gauge gravitational self-force calculation for $\psi$ at first order in the mass ratio. The latter provides new numerical data in the strong-field regime to inform the Effective One-Body model of the gravitational two-body problem. We conclude that $\psi$ complements the Detweiler redshift $z$ as a key invariant quantity characterizing eccentric orbits in the gravitational two-body problem.

Improved next-to-leading order tidal heating and torquing of a Kerr black hole

We calculate the energy and angular momentum fluxes across the event horizon of a tidally deformed, rapidly rotating black hole to next-to-leading order in the curvature of the external spacetime. These are expressed in terms of tidal quadrupole moments and their time derivatives, which provide a characterization of a generic tidal environment. As an application of our results, we provide an expression for the energy and angular-momentum fluxes across the horizon when the black hole is a member of a binary system on a slowly-moving, quasi-circular orbit. Our expressions are accurate to 1.5 post-Newtonian order beyond the leading-order fluxes, but they are valid for arbitrary mass ratios. We compare our results to those previously obtained in the case of an extreme mass ratio binary, and find that they do not agree at the 1.5 post-Newtonian order. We investigate a number of possible sources for this discrepancy, but are ultimately unable to resolve it.

Modeling Mergers of Known Galactic Systems of Binary Neutron Stars [Cross-Listing]

We present a study of the merger of six different known galactic systems of binary neutron stars (BNS) of unequal mass with a mass ratio between $0.75$ and $0.99$. Specifically, these systems are J1756-2251, J0737-3039A, J1906+0746, B1534+12, J0453+1559 and B1913+16. We follow the dynamics of the merger from the late stage of the inspiral process up to $\sim$ 20 ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH), using a semi-realistic equation of state (EOS), namely the seven-segment piece-wise polytropic SLy with a thermal component. For the most extreme of these systems ($q=0.75$, J0453+1559), we also investigate the effects of different EOSs: APR4, H4, and MS1. Our numerical simulations are performed using only publicly available open source code such as, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. We show results on the gravitational wave signals, spectrogram and frequencies of the BNS after the merger and the BH properties in the two cases in which the system collapse within the simulated time.

Modeling Mergers of Known Galactic Systems of Binary Neutron Stars

We present a study of the merger of six different known galactic systems of binary neutron stars (BNS) of unequal mass with a mass ratio between $0.75$ and $0.99$. Specifically, these systems are J1756-2251, J0737-3039A, J1906+0746, B1534+12, J0453+1559 and B1913+16. We follow the dynamics of the merger from the late stage of the inspiral process up to $\sim$ 20 ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH), using a semi-realistic equation of state (EOS), namely the seven-segment piece-wise polytropic SLy with a thermal component. For the most extreme of these systems ($q=0.75$, J0453+1559), we also investigate the effects of different EOSs: APR4, H4, and MS1. Our numerical simulations are performed using only publicly available open source code such as, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. We show results on the gravitational wave signals, spectrogram and frequencies of the BNS after the merger and the BH properties in the two cases in which the system collapse within the simulated time.

Modeling Mergers of Known Galactic Systems of Binary Neutron Stars [Cross-Listing]

We present a study of the merger of six different known galactic systems of binary neutron stars (BNS) of unequal mass with a mass ratio between $0.75$ and $0.99$. Specifically, these systems are J1756-2251, J0737-3039A, J1906+0746, B1534+12, J0453+1559 and B1913+16. We follow the dynamics of the merger from the late stage of the inspiral process up to $\sim$ 20 ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH), using a semi-realistic equation of state (EOS), namely the seven-segment piece-wise polytropic SLy with a thermal component. For the most extreme of these systems ($q=0.75$, J0453+1559), we also investigate the effects of different EOSs: APR4, H4, and MS1. Our numerical simulations are performed using only publicly available open source code such as, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. We show results on the gravitational wave signals, spectrogram and frequencies of the BNS after the merger and the BH properties in the two cases in which the system collapse within the simulated time.

Parameter estimates in binary black hole collisions using neural networks

We present an algorithm based on artificial neural networks (ANNs), that estimates the mass ratio in a binary black hole collision out of given Gravitational Wave (GW) strains. In this analysis, the ANN is trained with a sample of GW signals generated with numerical simulations. The effectiveness of the algorithm is evaluated with GWs generated also with simulations for given mass ratios unknown to the ANN. We measure the accuracy of the algorithm in the interpolation and extrapolation regimes. We present the results for noise free signals and signals contaminated with Gaussian noise, in order to foresee the dependence of the method accuracy in terms of the signal to noise ratio.

Parameter estimates in binary black hole collisions using neural networks [Cross-Listing]

We present an algorithm based on artificial neural networks (ANNs), that estimates the mass ratio in a binary black hole collision out of given Gravitational Wave (GW) strains. In this analysis, the ANN is trained with a sample of GW signals generated with numerical simulations. The effectiveness of the algorithm is evaluated with GWs generated also with simulations for given mass ratios unknown to the ANN. We measure the accuracy of the algorithm in the interpolation and extrapolation regimes. We present the results for noise free signals and signals contaminated with Gaussian noise, in order to foresee the dependence of the method accuracy in terms of the signal to noise ratio.

Whimper of a Bang: Documenting the Final Days of the Nearby Type Ia Supernova 2011fe

Using the Hubble Space Telescope (HST) and the Large Binocular Telescope, we followed the evolution of the Type Ia supernova (SN Ia) 2011fe for an unprecedented 1622 days past B-band maximum light and over a factor of 5 million in flux. At 1622 days, the 4000 - 17000 \AA{} quasi-bolometric luminosity is just ($710 \pm 30$ $L_{\odot}$). By measuring the late-time quasi-bolometric light curve we present the first confident detection of 57Co decay in a SN Ia light curve and estimate a mass ratio of log(57Co/56Co) = -1.62+0.08. We do not have a clean detection of 55Fe, but find a limit of 55Fe/57Co $< 0.3$ with 90$\%$ confidence. These abundance ratios provide unique constraints on the progenitor system because the central density of the exploding white dwarf(s) dictates these nucleosynthetic yields. The observed ratios strongly prefer the lower central densities of double-degenerate models (55Fe/57Co = 0.27) over the higher central densities of near Chandrasekhar-mass single-degenerate models (55Fe/57Co = 0.68). We will continue to observe SN 2011fe for another $\sim$900 days with HST and possibly beyond.

The photometric investigation of the newly discovered W UMa type binary system GSC 03122-02426

The $B$ $V$ $R_c$ $I_c$ bands light curves of the newly discovered binary system \astrobj{GSC 03122-02426} are obtained and analyzed using the Wilson-Devinney (W-D) code. The solutions suggest that the mass ratio of the binary system is $q = 2.70$ and the less massive component is $422K$ hotter than the more massive one. We conclude that \astrobj{GSC 03122-02426} is a W-subtype shallow contact (with a contact degree of $f = 15.3\,\%$) binary system. It may be a newly formed contact binary system which is just under geometrical contact and will evolve to be a thermal contact binary system. The high orbital inclination ($i = 81.6^{\circ}$) implies that \astrobj{GSC 03122-02426} is a total eclipsing binary system and the photometric parameters obtained by us are quite reliable. We also estimate the absolute physical parameters of the two components in \astrobj{GSC 03122-02426}, which will provide fundamental information for the research of contact binary systems. The formation and evolutionary scenario of \astrobj{GSC 03122-02426} is discussed.

The first photometric analysis of the W-subtype contact binary UCAC4 436-062932 with O'Connell effect

Two sets of light curves in $V$ $R_c$ $I_c$ bands for a newly discovered binary system UCAC4 436-062932 are obtained and analyzed using the Wilson-Devinney (W-D) code. The two sets of light curves get almost consistent results. The determined mass ratio is about $q = 2.7$ and the less massive component is nearly $250K$ hotter than the more massive one. The solutions conclude that UCAC4 436-062932 is a W-subtype shallow contact (with a contact degree of $f = 20\,\%$) binary system. Since the O'Connell effect appears on one set of the light curves, theories proposed to explain the effect are discussed. We assume that spot model may be the more plausible one to the O'Connell effect appeared on the asymmetric light curves of the binary system UCAC4 436-062932. Therefore, we add a cool spot on the surface of the more massive star (component with lower effective temperature) and get a quite approving results for the light curve fitting. It will provide evidence to support the spot model in the explanatory mechanism of O'Connell effect.

A Solar-type Stellar Companion to a Deep Contact Binary in a Quadruple System

The four-color ($B$ $V$ $R_c$ $I_c$) light curves of V776 Cas are presented and analyzed using the Wilson-Devinney (W-D) method. It is discovered that V776 Cas is an early F-type (F2V) overcontact binary with a very high contact degree ($ f=64.6\,\%$) and an extremely low mass ratio ($q=0.130$), which indicate that it is at the final evolutionary stage of cool short-period binaries. The mass of the primary and secondary stars are calculated to be $M_1 = 1.55(\pm0.04)M_\odot$, $M_2 = 0.20(\pm0.01)M_\odot$. V776 Cas is supposed to be formed from an initially detached binary system via the loss of angular momentum due to the magnetic wind. The initial mass of the present primary and secondary components are calculated to be $M_{1i} = 0.86(\pm0.10)M_\odot$ and $M_{2i} = 2.13(\pm0.04)M_\odot$. The observed-calculated ($O$-$C$) curve exhibits a cyclic period variation, which is due to the light-travel time effect (LTTE) caused by the presence of a third component with a period of 23.7 years. The mass of the third component is estimated to be $M_3 = 1.04(\pm0.03)M_\odot$ and the orbital inclination of the third component is calculated to be $i' = 33.1^{\circ}$. The distance of the binary system to the mass center of the triple system is calculated to be $a'_{12} = 3.45AU$. The presence of the close-in tertiary component may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system.

Photometric and period investigation of the late F-type overcontact binary II UMa

II UMa is a late F-type (F5) contact binary with a close-in tertiary and a distant visual companion. According to the four-color ($B$ $V$ $R_c$ $I_c$) light curves' solutions of II UMa, it is a high fill-out (f=$86.6\,\%$) and low mass ratio ($q = 0.172$) contact binary system, which indicate that it is at the late evolutionary stage of late-type tidal-locked binary stars. The mass of the primary star and secondary one are calculated to be $M_1 = 1.99M_\odot$, $M_2 = 0.34M_\odot$. The primary star has evolved from ZAMS, but it is still before TAMS, and the secondary star is even more evolved. Considering the mass ratio ($M_3/M_1 = 0.67$) obtained by spectroscopic observations, the mass of the close-in tertiary is estimated to be $M_3 = 1.34M_\odot$. The period variations of the binary system is investigated for the first time. According to the observed-calculated ($O$-$C$) curve analysis, a continuous period increase at a rate of $dP/dt=4.88\times{10^{-7}}day\cdot year^{-1}$ is determined. It may be just a part of a cyclic period change, or the combinational period change of a parabolic variation and a cyclic one. More times of minimum light are needed to confirm this. The presence of the tertiary component may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system during the pre-contact stage.

Gravitational-wave tail effects to quartic non-linear order

Gravitational-wave tails are due to the backscattering of linear waves onto the space-time curvature generated by the total mass of the matter source. The dominant tails correspond to quadratic non-linear interactions and arise at the one-and-a-half post-Newtonian (1.5PN) order in the gravitational waveform. Also known are the "tails-of-tails", which are cubically non-linear effects appearing at the 3PN order in the waveform. Here we derive still higher non-linear tail effects, namely those associated with quartic non-linear interactions or "tails-of-tails-of-tails", which are shown to arise at the 4.5PN order. As an application we obtain at that order the complete coefficient in the total gravitational-wave energy flux of compact binary systems moving on circular orbits. Our result perfectly agrees with black-hole perturbation calculations in the limit of extreme mass ratio of the two compact objects.

Evolution of binary seeds in collapsing protostellar gas clouds

We perform three dimensional smoothed particle hydrodynamics (SPH) simulations of gas accretion onto the seeds of binary stars to investigate their short-term evolution. Our simulation setup is more realistic compared to the previous works by taking into account of dynamically evolving envelope with non-uniform distribution of gas density and angular momentum of accreting flow. Our initial condition includes a seed binary and a surrounding gas envelope, modelling the phase of core collapse of gas cloud when the fragmentation has already occurred. We assume that the seed binary has no eccentricity and no growth by gas accretion. The envelope is assumed to be an isothermal gas with no self-gravity. We run multiple simulations with different values of initial mass ratio $q_0$ (the ratio of secondary over primary mass) and gas temperature, and find a critical value of $q_{\rm c} = 0.25$ which distinguishes the later evolution of mass ratio $q$ as a function of time. If $q_0 \ga q_{\rm c}$, the secondary seed grows faster and $q$ increases monotonically towards unity. If $q_0 \la q_{\rm c}$, on the other hand, the primary seed grows faster and $q$ is lower than $q_0$ at the end of the simulation. Based on our numerical results, we analytically calculate the long-term evolution of the seed binary including the growth of binary by gas accretion. We find that the seed binary with $q_0 \ga q_{\rm c}$ evolves towards an equal-mass binary star, and that with $q_0 \la q_{\rm c}$ evolves to a binary with an extreme value of $q$. Binary separation is a monotonically increasing function of time for any $q_0$, suggesting that the binary growth by accretion does not lead to the formation of close binaries.

Evolution of binary seeds in collapsing protostellar gas clouds [Replacement]

We perform three dimensional smoothed particle hydrodynamics (SPH) simulations of gas accretion onto the seeds of binary stars to investigate their short-term evolution. Our simulation setup is more realistic compared to the previous works by taking into account of dynamically evolving envelope with non-uniform distribution of gas density and angular momentum of accreting flow. Our initial condition includes a seed binary and a surrounding gas envelope, modelling the phase of core collapse of gas cloud when the fragmentation has already occurred. We assume that the seed binary has no eccentricity and no growth by gas accretion. The envelope is assumed to be an isothermal gas with no self-gravity. We run multiple simulations with different values of initial mass ratio $q_0$ (the ratio of secondary over primary mass) and gas temperature, and find a critical value of $q_{\rm c} = 0.25$ which distinguishes the later evolution of mass ratio $q$ as a function of time. If $q_0$ >~ $q_{\rm c}$, the secondary seed grows faster and $q$ increases monotonically towards unity. If $q_0$ <~ $q_{\rm c}$, on the other hand, the primary seed grows faster and $q$ is lower than $q_0$ at the end of the simulation. Based on our numerical results, we analytically calculate the long-term evolution of the seed binary including the growth of binary by gas accretion. We find that the seed binary with $q_0$ >~ $q_{\rm c}$ evolves towards an equal-mass binary star, and that with $q_0$ <~ $q_{\rm c}$ evolves to a binary with an extreme value of $q$. Binary separation is a monotonically increasing function of time for any $q_0$, suggesting that the binary growth by accretion does not lead to the formation of close binaries.

OGLE-2016-BLG-0596Lb: High-Mass Planet From High-Magnification Pure-Survey Microlensing Event

We report the discovery of a high mass-ratio planet $q=0.012$, i.e., 13 times higher than the Jupiter/Sun ratio. The host mass is not presently measured but can be determined or strongly constrained from adaptive optics imaging. The planet was discovered in a small archival study of high-magnification events in pure-survey microlensing data, which was unbiased by the presence of anomalies. The fact that it was previously unnoticed may indicate that more such planets lie in archival data and could be discovered by similar systematic study. In order to understand the transition from predominantly survey+followup to predominately survey-only planet detections, we conduct the first analysis of these detections in the observational $(s,q)$ plane. Here $s$ is projected separation in units of the Einstein radius. We find some evidence that survey+followup is relatively more sensitive to planets near the Einstein ring, but that there is no statistical difference in sensitivity by mass ratio.

Binary Neutron Star Mergers and Short Gamma-Ray Bursts: Effects of Magnetic Field Orientation, Equation of State, and Mass Ratio [Cross-Listing]

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

Binary Neutron Star Mergers and Short Gamma-Ray Bursts: Effects of Magnetic Field Orientation, Equation of State, and Mass Ratio [Replacement]

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

Binary Neutron Star Mergers and Short Gamma-Ray Bursts: Effects of Magnetic Field Orientation, Equation of State, and Mass Ratio [Replacement]

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

Binary Neutron Star Mergers and Short Gamma-Ray Bursts: Effects of Magnetic Field Orientation, Equation of State, and Mass Ratio

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

Photometric and H-alpha studies of two extreme mass-ratio short period contact binaries in the direction of open cluster Praesepe

We present the high precession photometric studies in V band and spectroscopic studies centered around H-alpha line for two extreme mass-ratio short period contact binaries ASAS J082243+1927.0 (V1) and ASAS J085710+1856.8 (V2). The variables in study are in the direction of open cluster Praesepe. From the light curve analysis V1 was found to be a low mass-ratio over-contact binary of A-type, with a mass ratio q 0.121 and fill-out factor f 72%, and V2 was found to be a low mass contact binary of W-type, with a mass-ratio q 1.29 and fill-out factor f 11% (marginal contact). The study of H-alpha absorption line profile of both the variables shows variation in equivalent widths (EWs) with orbital phases. The mean EWs of the H-alpha line were obtained as 1.6+or-0.13 A and 1.18+or-0.12 A for V1 and V2 respectively. The variation of H-alpha absorption with respect to phase is explained to be due to chromospheric activity in V1, as evident from the O Connell effect and that due to chromospheric flares, winds or photospheric source in V2. The parameters for both the binaries are also studied with respect to the large sample of well-studied contact binaries and possible mechanisms for merger in contact binaries of both low and high mass ratios emphasized.

An analytic family of post-merger template waveforms [Replacement]

Building on the analytical description of the post-merger (ringdown) waveform of coalescing, non-precessing, spinning, (BBHs) introduced in Phys.~Rev.~D~90, 024054 (2014), we propose an analytic, closed form, time-domain, representation of the $\ell=m=2$ gravitational radiation mode emitted after merger. This expression is given as a function of the component masses and dimensionless spins $(m_{1,2},\chi_{1,2})$ of the two inspiralling objects, as well as of the mass $M_{\rm BH}$ and (complex) frequency $\sigma_{1}$ of the fundamental quasi-normal mode of the remnant black hole. Our proposed template is obtained by fitting the post-merger waveform part of several publicly available numerical relativity simulations from the Simulating eXtreme Spacetimes (SXS) catalog and then suitably interpolating over (symmetric) mass ratio and spins. We show that this analytic expression accurately reproduces ($\sim$~0.01 rad) the phasing of the post-merger data of other datasets not used in its construction. This is notably the case of the spin-aligned run SXS:BBH:0305, whose intrinsic parameters are consistent with the 90\% credible intervals reported by the parameter-estimation followup of GW150914 in Phys. Rev. Lett. 116 (2016) no.24, 241102. Using SXS waveforms as "experimental" data, we further show that our template could be used on the actual GW150914 data to perform a new measure the complex frequency of the fundamental quasi-normal mode so to exploit the complete (high signal-to-noise-ratio) post-merger waveform. The proposed template could also allow the measure of the damping time of the second quasi-normal mode, thus helping to test the general-relativistic no-hair theorem.

Post-merger analytic templates for GW150914

Following the new analytic description of the postmerger (ringdown) waveform of coalescing, nonprecessing, spinning, black hole binaries (BBHs) introduced in Phys.~Rev.~D~90, 024054 (2014), we propose an analytic, closed form, time-domain, representation of the $\ell=m=2$ gravitational radiation mode emitted after merger. This expression is given as a function of the masses and dimensionless spins $(m_{1,2},\chi_{1,2})$ of the two inspiralling objects, as well as of the mass $M_{\rm BH}$ and (complex) frequency $\sigma_{1}$ of the fundamental quasi-normal mode of the final black hole. It is obtained by first fitting the postmerger waveform part of several numerical relativity simulations of the Simulating eXtreme Spacetimes (SXS) catalog and then suitably interpolating over (symmetric) mass ratio and spins. It is demonstrated that this analytic expression reproduces rather accurately the post-merger waveform of other SXS datasets not used to build it and, notably, the dataset SXS:BBH:0305, corresponding to a BBH with mass ratio $m_{1}/m_{2}=1.221$, $\chi_{1}=+0.33$ and $\chi_{2}=-0.4399$, which gives the closest match to the observed signal of GW150914. It would be interesting to use the postmerger analytic waveform template introduced here to perform a new measure of the frequency and damping time of the fundamental quasi-normal mode (and possibly also of individual masses and spins) of GW150914 that exploits the complete post-merger waveform signal and that is not restricted to only the late-time ringdown tail. In addition, the same analytic template could also be used to measure the damping time of the second quasi-normal mode, helping then to test the general-relativistic no-hair theorem.

Turbulence, Transport and Waves in Ohmic Dead Zones

We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change. We find that the residual Reynolds stress near the mid-plane decreases with increasing dead zone thickness, becoming negligible in cases where the active to dead mass ratio is less than a few percent. This implies that purely Ohmic dead zones would be vulnerable to episodic accretion outbursts via the mechanism of Martin & Lubow (2011). We show that even thick dead zones support a large amount of kinetic energy, but this energy is largely in fluid motions that are inefficient at angular momentum transport. Confirming results from Oishi & Mac Low (2009), the perturbed velocity field in the dead zone is dominated by an oscillatory, vertically extended circulation pattern with a low frequency compared to the orbital frequency. This disturbance has the properties predicted for the lowest order r mode in a hydrodynamic disk. We suggest that in a global disk similar excitations would lead to propagating waves, whose properties would vary with the thickness of the dead zone and the nature of the perturbations (isothermal or adiabatic). Flows with similar amplitudes would buckle settled particle layers and could reduce the efficiency of pebble accretion.

Simulating the dust content of galaxies: successes and failures

We present full volume cosmological simulations using the moving-mesh code AREPO to study the coevolution of dust and galaxies. We extend the dust model in AREPO to include thermal sputtering of grains and investigate the evolution of the dust mass function, the cosmic distribution of dust beyond the interstellar medium, and the dependence of dust-to-stellar mass ratio on galactic properties. The simulated dust mass function is well-described by a Schechter fit and lies closest to observations at $z = 0$. The radial scaling of projected dust surface density out to distances of $10 \, \text{Mpc}$ around galaxies with magnitudes $17 < i < 21$ is similar to that seen in Sloan Digital Sky Survey data. At $z = 0$, the predicted dust density of $\Omega_\text{dust} \approx 1.9 \times 10^{-6}$ lies in the range of $\Omega_\text{dust}$ values seen in low-redshift observations. We find that dust-to-stellar mass ratio anti-correlates with stellar mass for galaxies living along the star formation main sequence. Moreover, we estimate the $850 \, \mu\text{m}$ and $1.1 \, \text{mm}$ number density functions for simulated galaxies at $z = 1$ and analyse the relation between dust-to-stellar flux and mass ratios at $z = 0$. At high redshift, our model fails to produce enough dust-rich galaxies, and this tension is not alleviated by adopting a top-heavy initial mass function. We do not capture a decline in $\Omega_\text{dust}$ from $z = 2$ to $z = 0$, which suggests that dust production mechanisms more strongly dependent on star formation may help to produce the observed number of dusty galaxies near the peak of cosmic star formation.

Be discs in binary systems I. Coplanar orbits

Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio, eccentricity) and disc (viscosity) parameters affect the disc structure in coplanar systems. We simulate such binaries with the use of a smoothed particle hydrodynamics code. The main effects of the secondary on the disc are its truncation and the accumulation of material inwards of truncation. In circular or nearly circular prograde orbits, the disc maintains a rotating, constant in shape, configuration, which is locked to the orbital phase. The disc is smaller in size, more elongated and more massive for low viscosity parameter, small orbital separation and/or high mass ratio. Highly eccentric orbits are more complex, with the disc structure and total mass strongly dependent on the orbital phase. We also studied the effects of binarity in the disc continuum emission. Since the infrared and radio SED are sensitive to the disc size and density slope, the truncation and matter accumulation result in considerable modifications in the emergent spectrum. The decretion disc in circular binaries remains constant in shape and phase-locked. Eccentric binaries exhibit strong dependence on the distance to the secondary. We conclude that binarity can serve as an explanation for the variability exhibited in observations of Be stars, and that our model can be used to detect invisible companions.

Be discs in binary systems I. Coplanar orbits [Replacement]

Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio and eccentricity) and disc (viscosity) parameters affect the disc structure in coplanar systems. We simulate such binaries with the use of a smoothed particle hydrodynamics code. The main effects of the secondary on the disc are its truncation and the accumulation of material inwards of truncation. We find two cases with respect to the effects of eccentricity: (i) In circular or nearly circular prograde orbits, the disc maintains a rotating, constant in shape, configuration, which is locked to the orbital phase. The disc is smaller in size, more elongated and more massive for low viscosity parameter, small orbital separation and/or high mass ratio. (ii) Highly eccentric orbits are more complex, with the disc structure and total mass strongly dependent on the orbital phase and the distance to the secondary. We also study the effects of binarity in the disc continuum emission. Since the infrared and radio SED are sensitive to the disc size and density slope, the truncation and matter accumulation result in considerable modifications in the emergent spectrum. We conclude that binarity can serve as an explanation for the variability exhibited in observations of Be stars, and that our model can be used to detect invisible companions.

 

You need to log in to vote

The blog owner requires users to be logged in to be able to vote for this post.

Alternatively, if you do not have an account yet you can create one here.

Powered by Vote It Up

^ Return to the top of page ^