Normalized to: Ciecielag, P.
[5]
oai:arXiv.org:1907.06527 [pdf] - 2065185
Are most Cataclysmic Variables in Globular Clusters dynamically formed?
Submitted: 2019-07-15
We have been investigating populations of cataclysmic variables (CVs) in a
set of more than 300 globular cluster (GC) models evolved with the MOCCA code.
One of the main questions we have intended to answer is whether most CVs in GCs
are dynamically formed or not. Contrary to what has been argued for a long
time, we found that dynamical destruction of primordial CV progenitors is much
stronger in GCs than dynamical formation of CVs. In particular, we found that,
on average, the detectable CV population is predominantly composed of CVs
formed via a typical common envelope phase (> 70 per cent). However,
core-collapsed models tend to have higher fractions of bright CVs than
non-core-collapsed ones, which suggests then that the formation of CVs is
indeed slightly favoured through strong dynamical interactions in
core-collapsed GCs, due to the high stellar densities in their cores.
[7]
oai:arXiv.org:1702.05176 [pdf] - 1535915
AstroGrid-PL
Submitted: 2017-02-16
We summarise the achievements of the AstroGrid-PL project, which aims to
provide an infrastructure grid computing, distributed storage and Virtual
Observatory services to the Polish astronomical community. It was developed
from 2011-2015 as a domain grid component within the larger PLGrid Plus project
for scientific computing in Poland.
[8]
oai:arXiv.org:0810.4599 [pdf] - 1001016
Algorithmic comparisons of decaying, isothermal, supersonic turbulence
Kitsionas, S.;
Federrath, C.;
Klessen, R.;
Schmidt, W.;
Price, D.;
Dursi, J.;
Gritschneder, M.;
Walch, S.;
Piontek, R.;
Kim, J.;
Jappsen, A. -K.;
Ciecielag, P.;
Mac Low, M. -M.
Submitted: 2008-10-25, last modified: 2009-09-29
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)
[9]
oai:arXiv.org:0706.0619 [pdf] - 1903
Towards the optimal window for the 2MASS dipole
Submitted: 2007-06-05, last modified: 2008-04-18
A comparison of the 2MASS flux dipole to the CMB dipole can serve as a method
to constrain a combination of the cosmological parameter Omega_m and the
luminosity bias of the 2MASS survey. For this constraint to be as tight as
possible, it is necessary to maximize the correlation between the two dipoles.
This can be achieved by optimizing the survey window through which the flux
dipole is measured. Here we explicitly construct such a window for the 2MASS
survey. The optimization in essence reduces to excluding from the calculation
of the flux dipole galaxies brighter than some limiting magnitude K_min of the
near-infrared K_s band. This exclusion mitigates nonlinear effects and shot
noise from small scales, which decorrelate the 2MASS dipole from the CMB
dipole. Under the assumption of negligible shot noise we find that the optimal
value of K_min is about five. Inclusion of shot noise shifts the optimal K_min
to larger values. We present an analytical formula for shot noise for the 2MASS
flux dipole, to be used in follow-up work with 2MASS data.
The misalignment angle between the two dipoles is a sensitive measure of
their correlation: the higher the correlation, the smaller the expectation
value of the angle. A minimum of the misalignment is thus a sign of the optimal
gravity window. We model analytically the distribution function for the
misalignment angle and show that the misalignment estimated by Maller et al. is
consistent with the assumed underlying model (though it is greater than the
expectation value). We predict with about 90% confidence that the misalignment
will decrease if 2MASS galaxies brighter than K_min = 5 mag are excluded from
the calculation of the flux dipole. This prediction has been indirectly
confirmed by the results of Erdogdu et al. (ABRIDGED)
[10]
oai:arXiv.org:0707.0410 [pdf] - 2761
On gas drag in a circular binary system
Submitted: 2007-07-03
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.
[11]
oai:arXiv.org:astro-ph/0605237 [pdf] - 81925
A comparative study of disc-planet interaction
de Val-Borro, M.;
Edgar, R. G.;
Artymowicz, P.;
Ciecielag, P.;
Cresswell, P.;
D'Angelo, G.;
Delgado-Donate, E. J.;
Dirksen, G.;
Fromang, S.;
Gawryszczak, A.;
Klahr, H.;
Kley, W.;
Lyra, W.;
Masset, F.;
Mellema, G.;
Nelson, R. P.;
Paardekooper, S. -J.;
Peplinski, A.;
Pierens, A.;
Plewa, T.;
Rice, K.;
Schaefer, C.;
Speith, R.
Submitted: 2006-05-09, last modified: 2006-06-01
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.
[12]
oai:arXiv.org:astro-ph/0401195 [pdf] - 62065
Likelihood analysis of the Local Group acceleration revisited
Submitted: 2004-01-12
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.
[13]
oai:arXiv.org:astro-ph/0310895 [pdf] - 60558
Precision Analysis of the Local Group Acceleration
Submitted: 2003-10-31
We reexamine likelihood analyses 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 analyses. 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 analyses 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.
[14]
oai:arXiv.org:astro-ph/0010364 [pdf] - 38711
Gaussianity of Cosmic Velocity Fields and Linearity of the
Velocity-Gravity Relation
Submitted: 2000-10-18, last modified: 2003-02-21
We present a numerical study of the relation between the cosmic peculiar
velocity field and the gravitational acceleration field. We show that on mildly
non-linear scales (4-10 Mpc Gaussian smoothing), the distribution of the
Cartesian coordinates of each of these fields is well approximated by a
Gaussian. In particular, their kurtoses and negentropies are small compared to
those of the velocity divergence and density fields. We find that at these
scales the relation between the velocity and gravity field follows linear
theory to good accuracy. Specifically, the systematic errors in
velocity-velocity comparisons due to assuming the linear model do not exceed 6%
in beta. To correct for them, we test various nonlinear estimators of velocity
from density. We show that a slight modification of the alpha-formula proposed
by Kudlicki et al. yields an estimator which is essentially unbiased and has a
small variance.
[15]
oai:arXiv.org:astro-ph/0202296 [pdf] - 47790
New era in likelihood analyses of the Local Group acceleration
Submitted: 2002-02-15
In maximum-likelihood analyses of the Local Group (LG) acceleration, the
object describing nonlinear effects is the coherence function (CF), i.e. the
cross-correlation coefficient of the Fourier modes of the velocity and gravity
fields. We study the CF both analytically, using perturbation theory, and
numerically, using a hydrodynamic code. The dependence of the function on
Omega_m and the shape of the power spectrum is very weak. The only cosmological
parameter that the CF is strongly sensitive to is the normalization sigma_8 of
the underlying density field. Perturbative approximation for the function turns
out to be accurate as long as sigma_8 is smaller than about 0.3. For higher
normalizations we provide an analytical fit for the CF as a function of sigma_8
and the wavevector. The characteristic decoherence scale which our formula
predicts is an order of magnitude smaller than that determined by Strauss et
al. This implies that present likelihood constraints on cosmological parameters
from analyses of the LG acceleration are significantly tighter than hitherto
reported.
[16]
oai:arXiv.org:astro-ph/0109291 [pdf] - 44790
Local Group velocity versus gravity: the coherence function
Submitted: 2001-09-18, last modified: 2001-11-19
In maximum-likelihood analyses of the Local Group (LG) acceleration, the
object describing nonlinear effects is the coherence function (CF), i.e. the
cross-correlation coefficient of the Fourier modes of the velocity and gravity
fields. We study the CF both analytically, using perturbation theory, and
numerically, using a hydrodynamic code. The dependence of the function on
Omega_m and the shape of the power spectrum is very weak. The only cosmological
parameter that the CF is strongly sensitive to is the normalization sigma_8 of
the underlying density field. Perturbative approximation for the function turns
out to be accurate as long as sigma_8 is smaller than about 0.3. For higher
normalizations we provide an analytical fit for the CF as a function of sigma_8
and the wavevector. The characteristic decoherence scale which our formula
predicts is an order of magnitude smaller than that determined by Strauss et
al. This implies that present likelihood constraints on cosmological parameters
from analyses of the LG acceleration are significantly tighter than hitherto
reported.
[17]
oai:arXiv.org:astro-ph/0101078 [pdf] - 40216
Local Group Velocity Versus Gravity: Nonlinear Effects
Submitted: 2001-01-05
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.
[18]
oai:arXiv.org:astro-ph/0002029 [pdf] - 34372
High-resolution simulations and visualization of protoplanetary disks
Submitted: 2000-02-01, last modified: 2000-02-02
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.
