Source author record

Yuri Levin

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

38works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

38 published item(s)

preprint2022arXiv

Magnetar bursts due to Alfvén wave nonlinear breakout

The most common form of magnetar activity is short X-ray bursts, with durations from milliseconds to seconds, and luminosities ranging from $10^{36}$ to $10^{43}\ {\rm erg}\,{\rm s}^{-1}$. Recently, an X-ray burst from the galactic magnetar SGR 1935+2154 was detected to be coincident with two fast radio burst (FRB) like events from the same source, providing evidence that FRBs may be linked to magnetar bursts. Using fully 3D force-free electrodynamics simulations, we show that such magnetar bursts may be produced by Alfvén waves launched from localized magnetar quakes: a wave packet propagates to the outer magnetosphere, becomes nonlinear, and escapes the magnetosphere, forming an ultra-relativistic ejecta. The ejecta pushes open the magnetospheric field lines, creating current sheets behind it. Magnetic reconnection can happen at these current sheets, leading to plasma energization and X-ray emission. The angular size of the ejecta can be compact, $\lesssim 1$ sr if the quake launching region is small, $\lesssim 0.01$ sr at the stellar surface. We discuss implications for the FRBs and the coincident X-ray burst from SGR 1935+2154.

preprint2022arXiv

Revisiting Thermal Charge Carrier Refractive Noise in Semiconductor Optics for Gravitational-Wave Interferometers

The test masses in next-generation gravitational-wave interferometers may have a semiconductor substrate, most likely silicon. The stochastic motion of charge carriers within the semiconductor will cause random fluctuations in the material's index of refraction, introducing a noise source called Thermal Charge Carrier Refractive (TCCR) noise. TCCR noise was previously studied in 2020 by Bruns et al., using a Langevin force approach. Here we compute the power spectral density of TCCR noise by both using the Fluctuation-Dissipation theorem (FDT) and accounting for previously neglected effects of the standing wave of laser light which is produced inside the input test mass by its high-reflecting coatings. We quantify our results with parameters from Einstein Telescope, and show that at temperatures of 10 K the amplitude of TCCR noise is up to a factor of $\sqrt{2}$ times greater than what was previously claimed, and from 77 K to 300 K the amplitude is around 5 to 7 orders of magnitude lower than previously claimed when we choose to neglect the standing wave, and is up to a factor of 6 times lower if the standing wave is included. Despite these differences, we still conclude like Bruns et al. that TCCR noise should not be a limiting noise source for next-generation gravitational-wave interferometers.

preprint2020arXiv

Alfvén wave mode conversion in pulsar magnetospheres

The radio emission anomaly coincident with the 2016 glitch of the Vela pulsar may be caused by a star quake that launches Alfvén waves into the magnetosphere, disturbing the original radio emitting region. To quantify the lifetime of the Alfvén waves, we investigate a possible energy loss mechanism, the conversion of Alfvénwaves into fast magnetosonic waves. Using axisymmetric force-free simulations, we follow the propagation of Alfvén waves launched from the stellar surface with small amplitude into the closed zone of a force-free dipolar pulsar magnetosphere. We observe mode conversion happening in the ideal force-free regime. The conversion efficiency during the first passage of the Alfvén wave through the equator can be large, for waves that reach large amplitudes as they travel away from the star, or propagate on the field lines passing close to the Y-point. However, the conversion efficiency is reduced due to dephasing on subsequent passages and considerable Alfvén power on the closed field lines remains. Thus while some leakage into the fast mode happens, we need detailed understanding of the original quenching in order to say whether mode conversion alone can lead to reactivation of the pulsar on a short timescale.

preprint2020arXiv

Constraining temperature distribution inside LIGO test masses from frequencies of their vibrational modes

Thermal distortion of test masses, as well as thermal drift of their vibrational mode frequencies, present a major challenge for operation of the Advanced LIGO and Advanced VIRGO interferometers, reducing optical efficiency, which limits sensitivity and potentially causing instabilities which reduce duty-cycle. In this paper, we demonstrate that test-mass vibrational mode frequency data can be used to overcome some of these difficulties. First, we derive a general expression for the change in a mode frequency as a function of temperature distribution inside the test mass. Then we show how the mode frequency dependence on temperature distribution can be used to identify the wavefunction of observed vibrational modes. We then show how monitoring the frequencies of multiple vibrational modes allows the temperature distribution inside the test mass to be strongly constrained. Finally, we demonstrate using simulations, the potential to improve the thermal model of the test mass, providing independent and improved estimates of important parameters such as the coating absorption coefficient and the location of point absorbers.

preprint2020arXiv

Plasmoid ejection by Alfvén waves and the fast radio bursts from SGR 1935+2154

Using numerical simulations we show that low-amplitude Alfvén waves from a magnetar quake propagate to the outer magnetosphere and convert to "plasmoids" (closed magnetic loops) which accelerate from the star, driving blast waves into the magnetar wind. Quickly after its formation, the plasmoid becomes a thin relativistic pancake. It pushes out the magnetospheric field lines, and they gradually reconnect behind the pancake, generating a variable wind far stronger than the normal spindown wind of the magnetar. Repeating ejections drive blast waves in the amplified wind. We suggest that these ejections generate the simultaneous X-ray and radio bursts detected from SGR 1935+2154. A modest energy budget of the magnetospheric perturbation $\sim 10^{40}$ erg is sufficient to produce the observed bursts. Our simulation predicts a narrow (a few ms) X-ray spike from the magnetosphere, arriving almost simultaneously with the radio burst emitted far outside the magnetosphere. This timing is caused by the extreme relativistic motion of the ejecta.

preprint2020arXiv

Precessing flaring magnetar as a source of repeating FRB 180916.J0158+65

Recently, CHIME detected periodicity in the bursting rate of the repeating FRB 180916.J0158+65. In a popular class of models, the fast radio bursts (FRBs) are created by giant magnetic flares of a hyper-active magnetar driven by fast ambipolar diffusion in the core. We point out that in this scenario the magnetar is expected to precess freely with a period of hours to weeks. The internal magnetic field $B\sim 10^{16}$G deforms the star, and magnetic flares induce sudden changes in magnetic stresses. The resulting torques and displacements of the principal axes of inertia are capable of pumping a significant amplitude of precession. The anisotropy of the flaring FRB activity, combined with precession, implies a strong periodic modulation of the visible bursting rate. The ultra-strong field invoked in the magnetar model provides: (1) energy for the frequent giant flares, (2) the high rate of ambipolar diffusion, releasing the magnetic energy on the timescale $\sim 10^9$s, (3) the core temperature $T\approx 10^9$K, likely above the critical temperature for neutron superfluidity, (4) strong magnetospheric torques, which efficiently spin down the star, and (5) deformation with ellipticity $ε> 10^{-6}$, much greater than the rotational deformation. These conditions result in a precession with negligible viscous damping, and can explain the observed 16 day period in FRB 180916.J0158+65.

preprint2020arXiv

Simulation of a Compact Object with Outflows Moving Through a Gaseous Background

A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e., net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering AGN discs.

preprint2019arXiv

Negative Dynamical Friction on compact objects moving through dense gas

An over-dense wake is created by a gravitating object moving through a gaseous medium, and this wake pulls back on the object and slows it down. This is conventional dynamical friction in a gaseous medium. We argue that if the object drives a sufficiently powerful outflow, the wake is destroyed and instead an extended under-dense region is created behind the object. In this case the overall gravitational force is applied in the direction of the object's motion, producing a negative dynamical friction (NDF). Black holes in dense gas drive powerful outflows and may experience the NDF, although extensive numerical work is probably needed to demonstrate or refute this conclusively. NDF may be important for stellar-mass black holes and neutron stars inside `common envelopes' in binary systems, for stellar mass black holes inside AGN discs, or for massive black holes growing through super-Eddington accretion in early Universe.

preprint2016arXiv

Detecting gravitational-wave memory with LIGO: implications of GW150914

It may soon be possible for Advanced LIGO to detect hundreds of binary black hole mergers per year. We show how the accumulation of many such measurements will allow for the detection of gravitational-wave memory: a permanent displacement of spacetime that comes from strong-field, general relativistic effects. We estimate that Advanced LIGO operating at design sensitivity may be able to make a signal-to-noise ratio 3(5) detection of memory with ~35 (90) events with masses and distance similar to GW150914. We highlight the importance of incorporating higher-order gravitational-wave modes for parameter estimation of binary black hole mergers, and describe how our methods can also be used to detect higher-order modes themselves before Advanced LIGO reaches design sensitivity.

preprint2016arXiv

Fast and slow magnetic deflagration fronts in Type I X-ray bursts

Type I X-ray bursts are produced by thermonuclear runaways that develop on accreting neutron stars. Once one location ignites, the flame propagates across the surface of the star. Flame propagation is fundamental in order to understand burst properties like rise time and burst oscillations. Previous work quantified the effects of rotation on the front, showing that the flame propagates as a deflagration and that the front strongly resembles a hurricane. However the effect of magnetic fields was not investigated, despite the fact that magnetic fields strong enough to have an effect on the propagating flame are expected to be present on many bursters. In this paper we show how the coupling between fluid layers introduced by an initially vertical magnetic field plays a decisive role in determining the character of the fronts that are responsible for the Type I bursts. In particular, on a star spinning at 450 Hz (typical among the bursters) we test seed magnetic fields of $10^{7} - 10^{10}$ G and find that for the medium fields the magnetic stresses that develop during the burst can speed up the velocity of the burning front, bringing the simulated burst rise time close to the observed values. By contrast, in a magnetic slow rotator like IGR J17480--2446, spinning at 11 Hz, a seed field $\gtrsim 10^9$ G is required to allow localized ignition and the magnetic field plays an integral role in generating the burst oscillations observed during the bursts.

preprint2016arXiv

Gravitational-wave cosmology across 29 decades in frequency

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $Ω_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

preprint2016arXiv

Magnetar Outbursts from Avalanches of Hall Waves and Crustal Failures

We explore the interaction between Hall waves and mechanical failures inside a magnetar crust, using detailed one-dimentional models that consider temperature-sensitive plastic flow, heat transport and cooling by neutrino emission, as well as the coupling of the crustal motion to the magnetosphere. We find that the dynamics is enriched and accelerated by the fast, short-wavelength Hall waves that are emitted by each failure. The waves propagate and cause failures elsewhere, triggering avalanches. We argue that these avalanches are the likely sources of outbursts in transient magnetars.

preprint2015arXiv

Double tidal disruptions in galactic nuclei

A star on a nearly radial trajectory approaching a massive black hole (MBH) gets tidally disrupted if it comes sufficiently close to the MBH. Here we explore what happens to binary stars whose centers of mass approach the MBH on nearly radial orbits. The interaction with the MBH often leads to both stars being disrupted in sequence. We argue that such events could produce light curves that are substantially different from those of the single disruptions, with possible features such as two local maxima. Tidal forces from the MBH can also lead the binary components to collide; these merger products can form highly magnetized stars, whose subsequent tidal disruption may enable prompt jet formation.

preprint2015arXiv

Modelling Thermoelastic Distortion of Optics Using Elastodynamic Reciprocity

Thermoelastic distortion resulting from optical absorption by transmissive and reflective optics can cause unacceptable changes in optical systems that employ high power beams. In advanced-generation laser-interferometric gravitational wave detectors for example, optical absorption is expected to result in wavefront distortions that would compromise the sensitivity of the detector; thus necessitating the use of adaptive thermal compensation. Unfortunately, these systems have long thermal time constants and so predictive feed-forward control systems could be required - but the finite-element analysis is computationally expensive. We describe here the use of the Betti-Maxwell elastodynamic reciprocity theorem to calculate the response of linear elastic bodies (optics) to heating that has arbitrary spatial distribution. We demonstrate using a simple example, that it can yield accurate results in computational times that are significantly less than those required for finite-element analyses.

preprint2015arXiv

Rotational effects in thermonuclear Type I Bursts: equatorial crossing and directionality of flame spreading

In a previous study on thermonuclear (type I) nursts on accreting neutron stars we addressed and demonstrated the importance of the effects of rotation, through the Coriolis force, on the propagation of the burning flame. However, that study only analysed cases of longitudinal propagation, where the Coriolis force coefficient $2Ω\cosθ$ was constant. In this paper, we study the effects of rotation on propagation in the meridional (latitudinal) direction, where the Coriolis force changes from its maximum at the poles to zero at the equator. We find that the zero Coriolis force at the equator, while affecting the structure of the flame, does not prevent its propagation from one hemisphere to another. We also observe structural differences between the flame propagating towards the equator and that propagating towards the pole, the second being faster. In the light of the recent discovery of the low spin frequency of burster IGR~J17480-2446 rotating at 11 Hz (for which Coriolis effects should be negligible) we also extend our simulations to slow rotation.

preprint2014arXiv

Intermittency and Lifetime of the 625 Hz QPO in the 2004 Hyperflare from the Magnetar SGR 1806-20 as evidence for magnetic coupling between the crust and the core

Quasi-periodic oscillations (QPOs) detected in the 2004 giant flare from SGR 1806-20 are often interpreted as global magneto-elastic oscillations of the neutron star. There is, however, a large discrepancy between theoretical models, which predict that the highest frequency oscillations should die out rapidly, and the observations, which suggested that the highest-frequency signals persisted for ~100s in X-ray data from two different spacecraft. This discrepancy is particularly important for the high-frequency QPO at ~625 Hz. However, previous analyses did not systematically test whether the signal could also be there in much shorter data segments, more consistent with the theoretical predictions. Here, we test for the presence of the high-frequency QPO at 625 Hz in data from both the Rossi X-ray Timing Explorer (RXTE) and the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) systematically both in individual rotational cycles of the neutron star, as well as averaged over multiple successive rotational cycles at the same phase. We find that the QPO in the RXTE data is consistent with being only present in a single cycle, for a short duration of ~0.5s, whereas the RHESSI data are as consistent with a short-lived signal that appears and disappears as with a long-lived QPO. Taken together, this data provides evidence for strong magnetic interaction between the crust and the core.

preprint2014arXiv

Oblique Shock Breakout in Supernovae and Gamma-Ray Bursts: II. Numerical Solutions For Non-Relativistic Pattern Speeds

Non-spherical explosions develop non-radial flows as the pattern of shock emergence progresses across the stellar surface. In supernovae these flows can limit ejecta speeds, stifle shock breakout emission, and cause collisions outside the star. Similar phenomena occur in stellar and planetary collisions, tidal disruption events, accretion-induced collapses, and propagating detonations. We present two-dimensional, nested-grid Athena simulations of non-radial shock emergence in a frame comoving with the breakout pattern, focusing on the adiabatic, non-relativistic limit in a plane stratified envelope. We set boundary conditions using a known self-similar solution and explore the role of box size and resolution on the result. The shock front curves toward the stellar surface, and exhibits a kink from which weak discontinuities originate. Flow around the point of shock emergence is neither perfectly steady nor self-similar. Waves and vortices, which are not predominantly due to grid effects, emanate from this region. The post-shock flow is deflected along the stellar surface, and its pressure disturbs the stellar atmosphere upstream of the emerging shock. We use the numerical results and their analytical limits to predict the effects of radiation transfer and gravity, which are not included in our simulations.

preprint2014arXiv

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

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

preprint2014arXiv

On the Origin of Young Stars at the Galactic Center

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

preprint2014arXiv

Pulsar scintillations from corrugated reconnection sheets in the ISM

We show that surface waves along interstellar current sheets closely aligned with the line of sight lead to pulsar scintillation properties consistent with those observed. This mechanism naturally produces the length and density scales of the ISM scattering lenses that are required to explain the magnitude and dynamical spectrum of the scintillations. In this picture, the parts of warm ionized interstellar medium that are responsible for the scintillations are relatively quiescent, with scintillation and scattering resulting from weak waves propagating along magnetic domain boundary current sheets, which are both expected from helicity conservation and have been observed in numerical simulations. The model statistically predicts the spacing and amplitudes of inverted parabolic arcs seen in Fourier-transformed dynamical spectra of strongly scintillating pulsars with only 3 parameters. Multi-frequency, multi-epoch low frequency VLBI observations can quantitatively test this picture. If successful, in addition to mapping the ISM, this may open the door to precise nanoarcsecond pulsar astrometry, distance measurements, and emission studies using these 10AU interferometers in the sky.

preprint2014arXiv

Thermoplastic waves in magnetars

Magnetar activity is generated by shear motions of the neutron star surface, which relieve internal magnetic stresses. An analogy with earthquakes and faults is problematic, as the crust is permeated by strong magnetic fields, which greatly constrain crustal displacements. We describe a new deformation mechanism that is specific to strongly magnetized neutron stars. The magnetically stressed crust begins to move because of a thermoplastic instability, which launches a wave that shears the crust and burns its magnetic energy. The propagating wave front resembles the deflagration front in combustion physics. We describe the conditions for the instability, the front structure and velocity, and discuss implications for observed magnetar activity.

preprint2013arXiv

Oblique Shock Breakout in Supernovae and Gamma-Ray Bursts: I. Dynamics and Observational Implications

In a non-spherical stellar explosion, non-radial motions become important near the stellar surface. For realistic deviations from spherical symmetry, non-radial flow dramatically alters the dynamics and emission of shock emergence on a significant fraction of the surface. The breakout flash is stifled, ejecta speeds are limited, and matter is cast sideways. Non-radial ejection allows for collisions outside the star, which may engender a new type of transient. Strongly oblique breakouts are most easily produced in compact stellar progenitors, such as white dwarfs and stripped-envelope core collapse supernovae. We study the shock structure and post-shock acceleration using conservation laws, a similarity analysis, and an approximate theory for oblique shocks. The shock is likely to extend vertically from the stellar surface, then kink before joining a deep asymptotic solution. Outflow from the region crossed by an oblique shock is probably unsteady and may affect the surface ahead of the main shock. We comment on the implications for several notable explosions in which the non-spherical dynamics described in this paper are likely to play an important role. We also briefly consider relativistic and superluminal pattern speeds.

preprint2012arXiv

Creep events and creep noise in gravitational-wave interferometers: basic formalism and stationary limit

In gravitational-wave interferometers, test masses are suspended on thin fibers which experience considerable tension stress. Sudden microscopic stress release in a suspension fiber, which I call a 'creep event', would excite motion of the test mass that would be coupled to the interferometer's readout. The random test-mass motion due to a time-sequence of creep events is referred to as 'creep noise'. In this paper I present an elasto-dynamic calculation for the test-mass motion due to a creep event. I show that within a simple suspension model, the main coupling to the optical readout occurs via a combination of a "dc" horizontal displacement of the test mass, and excitation of the violin and pendulum modes, and not, as was thought previously, via lengthening of the fiber. When the creep events occur sufficiently frequently and their statistics is time-independent, the creep noise can be well-approximated by a stationary Gaussian random process. I derive the functional form of the creep noise spectral density in this limit, with the restrictive assumption that the creep events are statistically independent from each other.

preprint2012arXiv

Magnetically-levitating disks around supermassive black holes

In this paper we report on the formation of magnetically-levitating accretion disks around supermassive black holes. The structure of these disks is calculated by numerically modelling tidal disruption of magnetized interstellar gas clouds. We find that the resulting disks are entirely supported by the pressure of the magnetic fields against the component of gravitational force directed perpendicular to the disks. The magnetic field shows ordered large-scale geometry that remains stable for the duration of our numerical experiments extending over 10% of the disk lifetime. Strong magnetic pressure allows high accretion rate and inhibits disk fragmentation. This in combination with the repeated feeding of manetized molecular clouds to a supermassive black hole yields a possible solution to the long-standing puzzle of black hole growth in the centres of galaxies.

preprint2012arXiv

On the dynamics of mechanical failures in magnetized neutron-star crusts

We consider the dynamics of a mechanical failure induced by a shear stress in a strongly magnetized neutron-star crust. We show that even if the elastic properties of the crust allow the creation of a shear crack, the strongly sheared magnetic field around the crack leads to a back-reaction from the Lorentz force which does not allow large relative displacement of the crack surfaces. Instead, the global evolution of the crack proceeds on a slow resistive time scale, and is unable to release any substantial mechanical energy. Our calculations demostrate that for {\it some} magnetic-field configurations, the magnetic forces cause, effectively, a plastic deformation of the crust when the resulting elastic shear stress exceeds the critical value for mechanical failure.

preprint2012arXiv

Secular Dynamical Anti-Friction in Galactic Nuclei

We identify a gravitational-dynamical process in near-Keplerian potentials of galactic nuclei that occurs when an intermediate-mass black hole (IMBH) is migrating on an eccentric orbit through the stellar cluster towards the central supermassive black hole (SMBH). We find that, apart from conventional dynamical friction, the IMBH experiences an often much stronger systematic torque due to the secular (i.e., orbit-averaged) interactions with the cluster's stars. The force which results in this torque is applied, counterintuitively, in the same direction as the IMBH's precession and we refer to its action as "secular-dynamical anti-friction" (SDAF). We argue that SDAF, and not the gravitational ejection of stars, is responsible for the IMBH's eccentricity increase seen in the initial stages of previous N-body simulations. Our numerical experiments, supported by qualitative arguments, demonstrate that (1) when the IMBH's precession direction is artificially reversed, the torque changes sign as well, which decreases the orbital eccentricity, (2) the rate of eccentricity growth is sensitive to the IMBH migration rate, with zero systematic eccentricity growth for an IMBH whose orbit is artificially prevented from inward migration, and (3) SDAF is the strongest when the central star cluster is rapidly rotating. This leads to eccentricity growth/decrease for the clusters rotating in the opposite/same direction relative to the IMBH's orbital motion.

preprint2012arXiv

Understanding and analysing time-correlated stochastic signals in pulsar timing

Although it is widely understood that pulsar timing observations generally contain time-correlated stochastic signals (TCSSs; red timing noise is of this type), most data analysis techniques that have been developed make an assumption that the stochastic uncertainties in the data are uncorrelated, i.e. "white". Recent work has pointed out that this can introduce severe bias in determination of timing-model parameters, and that better analysis methods should be used. This paper presents a detailed investigation of timing-model fitting in the presence of TCSSs, and gives closed expressions for the post-fit signals in the data. This results in a Bayesian technique to obtain timing-model parameter estimates in the presence of TCSSs, as well as computationally more efficient expressions of their marginalised posterior distribution. A new method to analyse hundreds of mock dataset realisations simultaneously without significant computational overhead is presented, as well as a statistically rigorous method to check the internal consistency of the results. As a by-product of the analysis, closed expressions of the rms introduced by a stochastic background of gravitational-waves in timing-residuals are obtained. Using $T$ as the length of the dataset, and $h_c(1\rm{yr}^{-1})$ as the characteristic strain, this is: $σ_{\rm GWB}^2 = h_{c}(1\rm{yr}^{-1})^2 (9\sqrt[3]{2π^4}Γ(-10/3) / 8008) \rm{yr}^{-4/3} T^{10/3}$.

preprint2011arXiv

Magnetar Oscillations II: spectral method

The seismological dynamics of magnetars is largely determined by a strong hydro-magnetic coupling between the solid crust and the fluid core. In this paper we set up a "spectral" computational framework in which the magnetar's motion is decomposed into a series of basis functions which are associated with the crust and core vibrational eigenmodes. A general-relativistic formalism is presented for evaluation of the core Alfven modes in the magnetic-flux coordinates, as well for eigenmode computation of a strongly magnetized crust of finite thickness. By considering coupling of the crustal modes to the continuum of Alfven modes in the core, we construct a fully relativistic dynamical model of the magnetar which allows: i) Fast and long simulations without numerical dissipation. ii) Very fine sampling of the stellar structure. We find that the presence of strong magnetic field in the crust results in localizing of some high-frequency crustal elasto-magnetic modes with the radial number n>1 to the regions of the crust where the field is nearly horizontal. While the hydro-magnetic coupling of these localized modes to the Alfven continuum in the core is reduced, their energy is drained on a time-scale much less than 1 second. Therefore the puzzle of the observed QPOs with frequencies larger than 600 Hz still stands.

preprint2011arXiv

On the excitation of f-modes and torsional modes by magnetar giant flares

Magnetar giant flares may excite vibrational modes of neutron stars. Here we compute an estimate of initial post-flare amplitudes of both the torsional modes in the magnetar's crust and of the global f-modes. We show that while the torsional crustal modes can be strongly excited, only a small fraction of the flare's energy is converted directly into the lowest-order f-modes. For a conventional model of a magnetar, with the external magnetic field of about 10^{15} Gauss, the gravitational-wave detection of these f-modes with advanced LIGO is unlikely.

preprint2011arXiv

Secular Stellar Dynamics near a Massive Black Hole

The angular momentum evolution of stars close to massive black holes (MBHs) is driven by secular torques. In contrast to two-body relaxation, where interactions between stars are incoherent, the resulting resonant relaxation (RR) process is characterized by coherence times of hundreds of orbital periods. In this paper, we show that all the statistical properties of RR can be reproduced in an autoregressive moving average (ARMA) model. We use the ARMA model, calibrated with extensive N-body simulations, to analyze the long-term evolution of stellar systems around MBHs with Monte Carlo simulations. We show that for a single-mass system in steady-state, a depression is carved out near an MBH as a result of tidal disruptions. Using Galactic center parameters, the extent of the depression is about 0.1 pc, of similar order to but less than the size of the observed "hole" in the distribution of bright late-type stars. We also find that the velocity vectors of stars around an MBH are locally not isotropic. In a second application, we evolve the highly eccentric orbits that result from the tidal disruption of binary stars, which are considered to be plausible precursors of the "S-stars" in the Galactic center. We find that RR predicts more highly eccentric (e > 0.9) S-star orbits than have been observed to date.

preprint2010arXiv

Global Gravitational Instabilities in Discs with Infall

Gravitational instability plays an important role in driving gas accretion in massive protostellar discs. Particularly strong is the global gravitational instability, which arises when the disc mass is of order 0.1 of the mass of the central star and has a characteristic spatial scale much greater than the disc's vertical scale-height. In this paper we use three-dimensional numerical hydrodynamics to study the development of gravitational instabilities in a disc which is embedded in a dense, gaseous envelope. We find that global gravitational instabilities are the dominant mode of angular momentum transport in the disc with infall, in contrast to otherwise identical isolated discs. The accretion torques created by low-order, global modes of the gravitational instability in a disc subject to infall are larger by a factor of several than an isolated disc of the same mass. We show that this global gravitational instability is driven by the strong vertical shear at the interface between the disc and the envelope, and suggest that this process may be an important means of driving accretion on to young stars.

preprint2010arXiv

Magnetar Oscillations I: strongly coupled dynamics of the crust and the core

Quasi-Periodic Oscillations (QPOs) observed during Soft Gamma Repeaters giant flares are commonly interpreted as the torsional oscillations of magnetars. The oscillatory motion is influenced by the strong interaction between the shear modes of the crust and Alfven-like modes in the core. We study the dynamics which arises through this interaction, and present several new results: (1) We show that global {\it edge modes} frequently reside near the edges of the core Alfven continuum. (2) We compute the magnetar's oscillatory motion for realistic axisymmetric magnetic field configurations and core density profiles, but with a simplified model of the elastic crust. We show that one may generically get multiple gaps in the Alfven continuum. One obtains discrete global {\it gap modes} if the crustal frequencies belong to the gaps. (3) We show that field tangling in the core enhances the role of the core discrete Alfven modes and reduces the role of the core Alfven continuum in the overall oscillatory dynamics of the magnetar. (4) We demonstrate that the system displays transient and/or drifting QPOs when parts of the spectrum of the core Alfven modes contain discrete modes which are densely and regularly spaced in frequency. (5) We show that if the neutrons are coupled into the core Alfven motion, then the post-flare crustal motion is strongly damped and has a very weak amplitude. Thus magnetar QPOs give evidence that the proton and neutron components in the core are dynamically decoupled and that at least one of them is a quantum fluid. (6) We show that it is difficult to identify the high-frequency 625 Hz QPO as being due to the physical oscillatory mode of the magnetar, if the latter's fluid core consists of the standard proton-neutron-electron mixture and is magnetised to the same extent as the crust. (Abstract abridged)

preprint2010arXiv

The European Pulsar Timing Array: current efforts and a LEAP toward the future

The European Pulsar Timing Array (EPTA) is a multi-institutional, multi-telescope collaboration, with the goal of using high-precision pulsar timing to directly detect gravitational waves. In this article we discuss the EPTA member telescopes, current achieved timing precision, and near-future goals. We report a preliminary upper limit to the amplitude of a gravitational wave background. We also discuss the Large European Array for Pulsars, in which the five major European telescopes involved in pulsar timing will be combined to provide a coherent array that will give similar sensitivity to the Arecibo radio telescope, and larger sky coverage.

preprint2009arXiv

A New Secular Instability of Eccentric Stellar Disks Around Supermassive Black Holes, with Application to the Galactic Center

We identify a new secular instability of eccentric stellar disks around supermassive black holes. We show that retrograde precession of the stellar orbits, due to the presence of a stellar cusp, induces coherent torques that amplify deviations of individual orbital eccentricities from the average, and thus drive all eccentricities away from their initial value. We investigate the instability using N-body simulations, and show that it can drive individual orbital eccentricities to significantly higher or lower values on the order of a precession time-scale. This physics is relevant for the Galactic center, where massive stars are likely to form in eccentric disks around the SgrA* black hole. We show that the observed bimodal eccentricity distribution of disk stars in the Galactic center is in good agreement with the distribution resulting from the eccentricity instability and demonstrate how the dynamical evolution of such a disk results in several of its stars acquiring high (1-e < 0.1) orbital eccentricity. Binary stars on such highly eccentric orbits would get tidally disrupted by the SgrA* black hole, possibly producing both S-stars near the black hole and high-velocity stars in the Galactic halo.

preprint2009arXiv

Gravitational-wave memory and pulsar timing arrays

Pulsar timing arrays (PTAs) are designed to detect gravitational waves with periods from several months to several years, e.g. those produced by by wide supermassive black-hole binaries in the centers of distant galaxies. Here we show that PTAs are also sensitive to mergers of supermassive black holes. While these mergers occur on a timescale too short to be resolvable by a PTA, they generate a change of metric due to non-linear gravitational-wave memory which persists for the duration of the experiment and could be detected. We develop the theory of the single-source detection by PTAs, and derive the sensitivity of PTAs to the gravitational-wave memory jumps. We show that mergers of $10^8M_{\odot}$ black holes are $2-σ$-detectable (in a direction, polarization, and time-dependent way) out to co-moving distances of $\sim 1$ billion light years. Modern prediction for black-hole merger rates imply marginal to modest chance of an individual jump detection by currently developed PTAs. The sensitivity is expected to be somewhat higher for futuristic PTA experiments with SKA.

preprint2009arXiv

Thermo-refractive and thermo-chemical noise in the beamsplitter of GEO600 gravitational-wave interferometer

Braginsky, Gorodetsky, and Vyatchanin have shown that thermo-refractive fluctuations are an important source of noise in interferometric gravitational-wave detectors. In particular, the thermo-refractive noise in the GEO600 beamsplitter is expected to make a substantial contribution to the interferometer's total noise budget. Here we present a new computation of the GEO600 thermo-refractive noise which takes into account the beam's elliptical profile and, more importantly, the fact that the laser beam induces a standing electromagnetic wave in the beamsplitter. The use of updated parameters results in the overall reduction of the calculated noise amplitude by a factor of about 5 in the low-frequency part of the GEO600 band, compared to the previous estimates. We also find, by contrast with previous calculations, that thermo-refractive fluctuations result in white noise between 600 Hz and 39 MHz, at a level of $8.5\cdot 10^{-24}$Hz$^{-1/2}$. Finally, we describe a new type of thermal noise, which we call the thermo-chemical noise. This is caused by a random motion of optically-active chemical impurities or structural defects in the direction along a steep intensity gradient of the standing wave. We discuss the potential relevance of the thermo-chemical noise for GEO600.

preprint2003arXiv

Stellar disk in the galactic center -- a remnant of a dense accretion disk?

Observations of the galactic center revealed a population of young massive stars within 0.4 pc from Sgr A* -- the presumed location of a supermassive black hole. The origin of these stars is a puzzle as their formation in citu should be suppressed by the black hole's tidal field. We find that out of 13 stars whose 3-dimensional velocities have been measured by Genzel et. al. (2000), 10 lie in a thin disk. The half-opening angle of the disk is consistent with zero within the measurement errors, and does not exceed 10 degrees. We propose that a recent burst of star formation has occurred in a dense gaseous disk around Sgr A*. Such a disk is no longer present because, most likely, it has been accreted by the central black hole. The three-dimensional orbit of S2, the young star closest to Sgr A*, has been recently mapped out with high precision. It is inclined to the stellar disk by 75 degrees. We find that the orbit should undergo Lense-Thirring precession with the period of (5/a) Myr, where a<1 is the dimensionless spin of the black hole. Therefore it is possible that originally S2 orbit lay in the disk plane. If so, we can constrain the black hole spin $a$ be greater than 0.2(t_{S2}/5 Myr), where t_{S2} is the age of S2.

preprint1999arXiv

Nonlinear r-modes in a spherical shell: issues of principle

We use a simple physical model to study the nonlinear behaviour of the r-mode instability. We assume that r-modes (Rossby waves) are excited in a thin spherical shell of rotating incompressible fluid. For this case, exact Rossby wave solutions of arbitrary amplitude are known. We find that: (a) These nonlinear Rossby waves carry ZERO physical angular momentum and positive physical energy, which is contrary to the folklore belief that the r-mode angular momentum and energy are negative. (b) Within our model, we confirm the differential drift reported by Rezzolla, Lamb and Shapiro (1999). Radiation reaction is introduced into the model by assuming that the fluid is electrically charged; r-modes are coupled to electromagnetic radiation through current (magnetic) multipole moments. We find that: (c) To linear order in the mode amplitude, r-modes are subject to the CFS instability, as expected. (d) Radiation reaction decreases the angular velocity of the shell and causes differential rotation (which is distinct from but similar in magnitude to the differential drift reported by Rezzolla et al.) prior to saturation of the r-mode growth. This is contrary to the phenomenological treatments to date, which assume that the loss of stellar angular momentum is accounted for by the r-mode growth. We demonstrate, for the first time, that r-mode radiation reaction leads to differential rotation. (e) We show that for l=2 r-mode electromagnetic radiation reaction is equivalent to gravitational radiation reaction in the lowest post-Newtonian order.