1,257 research outputs found
Strong Lensing Reconstruction
We present a general linear algorithm for measuring the surface mass density
1-\kappa from the observable reduced shear g=\gamma/(1-\kappa) in the strong
lensing regime. We show that in general, the observed polarization field can be
decomposed into ``electric'' and ``magnetic'' components, which have
independent and redundant solutions, but perfectly orthogonal noise properties.
By combining these solutions, one can increase the signal-to-noise ratio by
\sqrt{2}. The solutions allow dynamic optimization of signal and noise, both in
real and Fourier space (using arbitrary smoothing windows). Boundary conditions
have no effect on the reconstructions, apart from its effect on the
signal-to-noise. Many existing reconstruction techniques are recovered as
special cases of this framework. The magnetic solution has the added benefit of
yielding the global and local parity of the reconstruction in a single step.Comment: final accepted version for ApJ
Efficient approximations of neutrino physics for three-dimensional simulations of stellar core collapse
Neutrino transport in spherically symmetric models of stellar core collapse
and bounce has achieved a technically complete level, rewarded by the agreement
among independent groups that a multi-dimensional treatment of the
fluid-instabilities in the post-bounce phase is indispensable to model
supernova explosions. While much effort is required to develop a reliable
neutrino transport technique in axisymmetry, we explore neutrino physics
approximations and parameterizations for an efficient three-dimensional
simulation of the fluid-instabilities in the shock-heated matter that
accumulates between the accretion shock and the protoneutron star. We
demonstrate the reliability of a simple parameterization scheme in the collapse
phase and extend our 3D magneto-hydrodynamical collapse simulations to a
preliminary postbounce evolution. The growth of magnetic fields is
investigated.Comment: 5 pages, 4 figures, in Proceedings of "Nuclei in the Cosmos IX,
Geneva, Jun 25-30", associated movies are displayed at
http://www.physik.unibas.ch/~liebend/displa
Approaching the dynamics of hot nucleons in supernovae
All recent numerical simulations agree that stars in the main sequence mass
range of 9-40 solar masses do not produce a prompt hydrodynamic ejection of the
outer layers after core collapse and bounce. Rather they suggest that stellar
core collapse and supernova explosion are dynamically distinct astrophysical
events, separated by an unspectacular accretion phase of at least ~40 ms
duration. As long as the neutrinospheres remain convectively stable, the
explosion dynamics is determined by the neutrons, protons, electrons and
neutrinos in the layer of impact-heated matter piling up on the protoneutron
star. The crucial role of neutrino transport in this regime has been emphasized
in many previous investigations. Here, we search for efficient means to address
the role of magnetic fields and fluid instabilities in stellar core collapse
and the postbounce phase.Comment: 4 pages, contribution to Nuclei in the Cosmos VIII, Jul. 19-23,
submitted to Nucl. Phys.
Performance improvement of robots using a learning control scheme
Many applications of robots require that the same task be repeated a number of times. In such applications, the errors associated with one cycle are also repeated every cycle of the operation. An off-line learning control scheme is used here to modify the command function which would result in smaller errors in the next operation. The learning scheme is based on a knowledge of the errors and error rates associated with each cycle. Necessary conditions for the iterative scheme to converge to zero errors are derived analytically considering a second order servosystem model. Computer simulations show that the errors are reduced at a faster rate if the error rate is included in the iteration scheme. The results also indicate that the scheme may increase the magnitude of errors if the rate information is not included in the iteration scheme. Modification of the command input using a phase and gain adjustment is also proposed to reduce the errors with one attempt. The scheme is then applied to a computer model of a robot system similar to PUMA 560. Improved performance of the robot is shown by considering various cases of trajectory tracing. The scheme can be successfully used to improve the performance of actual robots within the limitations of the repeatability and noise characteristics of the robot
FISH: A 3D parallel MHD code for astrophysical applications
FISH is a fast and simple ideal magneto-hydrodynamics code that scales to ~10
000 processes for a Cartesian computational domain of ~1000^3 cells. The
simplicity of FISH has been achieved by the rigorous application of the
operator splitting technique, while second order accuracy is maintained by the
symmetric ordering of the operators. Between directional sweeps, the
three-dimensional data is rotated in memory so that the sweep is always
performed in a cache-efficient way along the direction of contiguous memory.
Hence, the code only requires a one-dimensional description of the conservation
equations to be solved. This approach also enable an elegant novel
parallelisation of the code that is based on persistent communications with MPI
for cubic domain decomposition on machines with distributed memory. This scheme
is then combined with an additional OpenMP parallelisation of different sweeps
that can take advantage of clusters of shared memory. We document the detailed
implementation of a second order TVD advection scheme based on flux
reconstruction. The magnetic fields are evolved by a constrained transport
scheme. We show that the subtraction of a simple estimate of the hydrostatic
gradient from the total gradients can significantly reduce the dissipation of
the advection scheme in simulations of gravitationally bound hydrostatic
objects. Through its simplicity and efficiency, FISH is as well-suited for
hydrodynamics classes as for large-scale astrophysical simulations on
high-performance computer clusters. In preparation for the release of a public
version, we demonstrate the performance of FISH in a suite of astrophysically
orientated test cases.Comment: 27 pages, 11 figure
Near term measurements with 21 cm intensity mapping: neutral hydrogen fraction and BAO at z<2
It is shown that 21 cm intensity mapping could be used in the near term to
make cosmologically useful measurements. Large scale structure could be
detected using existing radio telescopes, or using prototypes for dedicated
redshift survey telescopes. This would provide a measure of the mean neutral
hydrogen density, using redshift space distortions to break the degeneracy with
the linear bias. We find that with only 200 hours of observing time on the
Green Bank Telescope, the neutral hydrogen density could be measured to 25%
precision at redshift 0.54<z<1.09. This compares favourably to current
measurements, uses independent techniques, and would settle the controversy
over an important parameter which impacts galaxy formation studies. In
addition, a 4000 hour survey would allow for the detection of baryon acoustic
oscillations, giving a cosmological distance measure at 3.5% precision. These
observation time requirements could be greatly reduced with the construction of
multiple pixel receivers. Similar results are possible using prototypes for
dedicated cylindrical telescopes on month time scales, or SKA pathfinder
aperture arrays on day time scales. Such measurements promise to improve our
understanding of these quantities while beating a path for future generations
of hydrogen surveys.Comment: 6 pages, 5 figures. Submitted to Phys. Rev. D. Addressed reviewer
comments. Changed figure format, added more detailed technical discussion,
and added forecasts for aperture arrays. Added references
Pulsar timing arrays as imaging gravitational wave telescopes: angular resolution and source (de)confusion
Pulsar timing arrays (PTAs) will be sensitive to a finite number of
gravitational wave (GW) "point" sources (e.g. supermassive black hole
binaries). N quiet pulsars with accurately known distances d_{pulsar} can
characterize up to 2N/7 distant chirping sources per frequency bin \Delta
f_{gw}=1/T, and localize them with "diffraction limited" precision \delta\theta
\gtrsim (1/SNR)(\lambda_{gw}/d_{pulsar}). Even if the pulsar distances are
poorly known, a PTA with F frequency bins can still characterize up to
(2N/7)[1-(1/2F)] sources per bin, and the quasi-singular pattern of timing
residuals in the vicinity of a GW source still allows the source to be
localized quasi-topologically within roughly the smallest quadrilateral of
quiet pulsars that encircles it on the sky, down to a limiting resolution
\delta\theta \gtrsim (1/SNR) \sqrt{\lambda_{gw}/d_{pulsar}}. PTAs may be
unconfused, even at the lowest frequencies, with matched filtering always
appropriate.Comment: 7 pages, 1 figure, matches Phys.Rev.D versio
- …