-
We report on a new all-sky search for periodic gravitational waves in the
frequency band 475-2000 Hz and with a frequency time derivative in the range of
[-1.0e-8, +1e-9] Hz/s. Potential signals could be produced by a nearby spinning
and slightly non-axisymmetric isolated neutron star in our galaxy.
This search uses the data from Advanced LIGO's first observational run O1. No
gravitational wave signals were observed, and upper limits were placed on their
strengths. For completeness, results from the separately published low
frequency search 20-475 Hz are included as well.
Our lowest upper limit on worst-case (linearly polarized) strain amplitude
h_0 is 4e-25 near 170 Hz, while at the high end of our frequency range we
achieve a worst-case upper limit of 1.3e-24. For a circularly polarized source
(most favorable orientation), the smallest upper limit obtained is ~1.5e-25.
-
We report on an all-sky search for periodic gravitational waves in the
frequency band 20-475 Hz and with a frequency time derivative in the range of
[-1.0, +0.1]e-8 Hz/s. Such a signal could be produced by a nearby spinning and
slightly non-axisymmetric isolated neutron star in our galaxy. This search uses
the data from Advanced LIGO's first observational run, O1. No periodic
gravitational wave signals were observed, and upper limits were placed on their
strengths. The lowest upper limits on worst-case (linearly polarized) strain
amplitude h0 are 4e-25 near 170 Hz. For a circularly polarized source (most
favorable orientation), the smallest upper limits obtained are 1.5e-25. These
upper limits refer to all sky locations and the entire range of frequency
derivative values. For a population-averaged ensemble of sky locations and
stellar orientations, the lowest upper limits obtained for the strain amplitude
are 2.5e-25.
-
Contradicting results have been reported in the literature with respect to
the performance of the numerical techniques employed for the study of
supersonic turbulence. We aim at characterising the performance of different
particle-based and grid-based techniques on the modelling of decaying
supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and
three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally
analysed two calculations denoted as PHANTOM A and PHANTOM B using two
different implementations of artificial viscosity. Our analysis indicates that
grid codes tend to be less dissipative than SPH codes, though details of the
techniques used can make large differences in both cases. For example, the
Morris & Monaghan viscosity implementation for SPH results in less dissipation
(PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a
smaller diffusion parameter leads to less dissipation, but results in a larger
bottleneck effect (our ENZO versus FLASH runs). As a general result, we find
that by using a similar number of resolution elements N for each spatial
direction means that all codes (both grid-based and particle-based) show
encouraging similarity of all statistical quantities for isotropic supersonic
turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid
cells), while scales smaller than that are significantly affected by the
specific implementation of the algorithm for solving the equations of
hydrodynamics. At comparable numerical resolution, the SPH runs were on average
about ten times more computationally intensive than the grid runs, although
with variations of up to a factor of ten between the different SPH runs and
between the different grid runs. (abridged)
-
We investigate both analytically and numerically the motion of massless
particles orbiting primary star in a close circular binary system with
particular focus on the gas drag effects. These are the first calculations with
particles ranging in size from 1m to 10km, which account for the presence of a
tidally perturbed gaseous disk. We have found numerically that the radial mass
transport by the tidal waves plays a crucial role in the orbital evolution of
particles. Numerical results are confirmed analytical calculations that do not
assume anything about origin of the radial flow in the disk. We demonstrate
that the migration rate of a particle in a disk out of radial equilibrium is
enhanced due to the enhanced mass flux of gas colliding with the particle and
the migration is always directed inward regardless of the sign of the radial
gas velocity. Within the framework of the perturbation theory we derive
general, formulae for short-term variations of the particle semi-major axis,
eccentricity and inclination in such disk. The formulae account for departures
from axial symmetry by introducing effective components of the gas velocity.
They agree with numerical results within several percent. We have also found in
numerical simulations that the tidal waves introduce coherence in periastron
longitude and eccentricity for particles on neighbouring orbits. The degree of
the coherence depends on the particle size and on the distance from the primary
star, being most prominent for particles with 10m radius. The results are
important mainly in the context of planetesimal formation and, to a lesser
degree, during the early planetesimal accretion stage.
-
We perform numerical simulations of a disc-planet system using various
grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run
for a simple setup where Jupiter and Neptune mass planets on a circular orbit
open a gap in a protoplanetary disc during a few hundred orbital periods. We
compare the surface density contours, potential vorticity and smoothed radial
profiles at several times. The disc mass and gravitational torque time
evolution are analyzed with high temporal resolution. There is overall
consistency between the codes. The density profiles agree within about 5% for
the Eulerian simulations while the SPH results predict the correct shape of the
gap although have less resolution in the low density regions and weaker
planetary wakes. The disc masses after 200 orbital periods agree within 10%.
The spread is larger in the tidal torques acting on the planet which agree
within a factor 2 at the end of the simulation. In the Neptune case the
dispersion in the torques is greater than for Jupiter, possibly owing to the
contribution from the not completely cleared region close to the planet.
-
We reexamine likelihood analyzes of the Local Group (LG) acceleration, paying
particular attention to nonlinear effects. Under the approximation that the
joint distribution of the LG acceleration and velocity is Gaussian, two
quantities describing nonlinear effects enter these analyzes. The first one is
the coherence function, i.e. the cross-correlation coefficient of the Fourier
modes of gravity and velocity fields. The second one is the ratio of velocity
power spectrum to gravity power spectrum. To date, in all analyzes of the LG
acceleration the second quantity was not accounted for. Extending our previous
work, we study both the coherence function and the ratio of the power spectra.
With the aid of numerical simulations we obtain expressions for the two as
functions of wavevector and \sigma_8. Adopting WMAP's best determination of
\sigma_8, we estimate the most likely value of the parameter \beta and its
errors. As the observed values of the LG velocity and gravity, we adopt
respectively a CMB-based estimate of the LG velocity, and Schmoldt et al.'s
(1999) estimate of the LG acceleration from the PSCz catalog. We obtain \beta =
0.66^{+0.21}_{-0.07}; thus our errorbars are significantly smaller than those
of Schmoldt et al. This is not surprising, because the coherence function they
used greatly overestimates actual decoherence between nonlinear gravity and
velocity.
-
We use numerical simulations to study the relation between the velocity of
the Local Group and its gravitational acceleration. This relation serves as a
test for the kinematic origin of the CMB dipole and as a method for estimating
beta. We calculate the misalignment angle between the two vectors and compare
it to the observed value for the PSCz survey. The latter value is near the
upper limit of the 95% confidence interval for the angle; therefore, the
nonlinear effects are unlikely to be responsible for the whole observed
misalignment. We also study the relation between the amplitudes of the LG
velocity and gravity vectors. In an Omega=1 Universe, the smoothed gravity of
the LG turns out to be a biased low estimator of the LG (unsmoothed) velocity.
In an Omega=0.3 Universe, the estimator is biased high. The discussed biases
are, however, only a few per cent, thus the linear theory works to good
accuracy. The gravity-based estimator of the LG velocity has also a scatter
that limits the precision of the estimate of beta in the LG velocity-gravity
comparisons. The random error of beta due to nonlinear effects amounts to
several per cent.
-
A problem of mass flow in the immediate vicinity of a planet embedded in a
protoplanetary disk is studied numerically in two dimensions. Large differences
in temporal and spatial scales involved suggest that a specialized
discretization method for solution of hydrodynamical equations may offer great
savings in computational resources, and can make extensive parameter studies
feasible. Preliminary results obtained with help of Adaptive Mesh Refinement
technique and high-order explicit Eulerian solver are presented. This
combination of numerical techniques appears to be an excellent tool which
allows for direct simulations of mass flow in vicinity of the accretor at
moderate computational cost. In particular, it is possible to resolve the
surface of the planet and to model the process of planet growth with minimal
set of assumptions. Some issues related to visualization of the results and
future prospects are discussed briefly.
We report on an all-sky search for periodic gravitational waves in the frequency band 20-475 Hz and with a frequency time derivative in the range of [-1.0, +0.1]e-8 Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run, O1. No periodic gravitational wave signals were observed, and upper limits were placed on their strengths. The lowest upper limits on worst-case (linearly polarized) strain amplitude h0 are 4e-25 near 170 Hz. For a circularly polarized source (most favorable orientation), the smallest upper limits obtained are 1.5e-25. These upper limits refer to all sky locations and the entire range of frequency derivative values. For a population-averaged ensemble of sky locations and stellar orientations, the lowest upper limits obtained for the strain amplitude are 2.5e-25.
Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)
We investigate both analytically and numerically the motion of massless particles orbiting primary star in a close circular binary system with particular focus on the gas drag effects. These are the first calculations with particles ranging in size from 1m to 10km, which account for the presence of a tidally perturbed gaseous disk. We have found numerically that the radial mass transport by the tidal waves plays a crucial role in the orbital evolution of particles. Numerical results are confirmed analytical calculations that do not assume anything about origin of the radial flow in the disk. We demonstrate that the migration rate of a particle in a disk out of radial equilibrium is enhanced due to the enhanced mass flux of gas colliding with the particle and the migration is always directed inward regardless of the sign of the radial gas velocity. Within the framework of the perturbation theory we derive general, formulae for short-term variations of the particle semi-major axis, eccentricity and inclination in such disk. The formulae account for departures from axial symmetry by introducing effective components of the gas velocity. They agree with numerical results within several percent. We have also found in numerical simulations that the tidal waves introduce coherence in periastron longitude and eccentricity for particles on neighbouring orbits. The degree of the coherence depends on the particle size and on the distance from the primary star, being most prominent for particles with 10m radius. The results are important mainly in the context of planetesimal formation and, to a lesser degree, during the early planetesimal accretion stage.