Three-dimensional radiative transfer calculations on an SIMD machine applied to accretion disks

We have developed a tool to solve the radiative transfer equation for a three-dimensional astrophysical object on the SIMD computer MasPar MP-1. With this tool we can rapidly calculate the image of such an object as seen from an arbitrary direction and at an arbitrary wavelength. Such images and spectra can then be used to directly compare observations with the model. This tool can be applied to many different areas in astrophysics, e.g., HI disks of galaxies and polarized radiative transfer of accretion columns onto white dwarfs. Here we use this tool to calculate the image and spectrum of a simple model of an accretion disk

The accretion flow in the discless intermediate polar V2400 Ophiuchi

RXTE observations confirm that the X-ray lightcurve of V2400 Oph is pulsed at the beat cycle, as expected in a discless intermediate polar. There are no X-ray modulations at the orbital or spin cycles, but optical line profiles vary with all three cycles. We construct a model for line-profile variations in a discless accretor, based on the idea that the accretion stream flips from one magnetic pole to the other, and show that this accounts for the observed behaviour over the spin and beat cycles. The minimal variability over the orbital cycle implies that 1) V2400 Oph is at an inclination of only ~10 deg, and 2) much of the accretion flow is not in a coherent stream, but is circling the white dwarf, possibly as a ring of denser, diamagnetic blobs. We discuss the light this sheds on disc formation in intermediate polars.Comment: 10 pages, 12 figures, To appear in MNRAS, includes low-res figures to reduce siz

A 3D radiative transfer framework: I. non-local operator splitting and continuum scattering problems

We describe a highly flexible framework to solve 3D radiation transfer problems in scattering dominated environments based on a long characteristics piece-wise parabolic formal solution and an operator splitting method. We find that the linear systems are efficiently solved with iterative solvers such as Gauss-Seidel and Jordan techniques. We use a sphere-in-a-box test model to compare the 3D results to 1D solutions in order to assess the accuracy of the method. We have implemented the method for static media, however, it can be used to solve problems in the Eulerian-frame for media with low velocity fields.Comment: A&A, in press. 14 pages, 19 figures. Full resolution figures available at ftp://phoenix.hs.uni-hamburg.de/preprints/3DRT_paper1.pdf HTML version (low res figures) at http://hobbes.hs.uni-hamburg.de/~yeti/PAPERS/3drt_paper1/index.htm

An Optimal Execution Problem with Market Impact

We study an optimal execution problem in a continuous-time market model that considers market impact. We formulate the problem as a stochastic control problem and investigate properties of the corresponding value function. We find that right-continuity at the time origin is associated with the strength of market impact for large sales, otherwise the value function is continuous. Moreover, we show the semi-group property (Bellman principle) and characterise the value function as a viscosity solution of the corresponding Hamilton-Jacobi-Bellman equation. We introduce some examples where the forms of the optimal strategies change completely, depending on the amount of the trader's security holdings and where optimal strategies in the Black-Scholes type market with nonlinear market impact are not block liquidation but gradual liquidation, even when the trader is risk-neutral.Comment: 36 pages, 8 figures, a modified version of the article "An optimal execution problem with market impact" in Finance and Stochastics (2014

V405 Aurigae: A High Magnetic Field Intermediate Polar

Our simultaneous multicolor (UBVRI) circular polarimetry has revealed nearly sinusoidal variation over the WD spin cycle, and almost symmetric positive and negative polarization excursions. Maximum amplitudes are observed in the B and V bands (+-3 %). This is the first time that polarization peaking in the blue has been discovered in an IP, and suggests that V405 Aur is the highest magnetic field IP found so far. The polarized flux spectrum is similar to those found in polars with magnetic fields in the range B ~ 25-50 MG. Our low resolution circular spectropolarimetry has given evidence of transient features which can be fitted by cyclotron harmonics n = 6, 7, and 8, at a field of B = 31.5 +- 0.8 MG, consistent with the broad-band polarized flux spectrum. Timings of the circular polarization zero crossovers put strict upper limits on WD spin period changes and indicate that the WD in V405 Aur is currently accreting closely at the spin equilibrium rate, with very long synchronization timescales, T_s > 10^9 yr. For the observed spin to orbital period ratio, P_{spin}/P_{orb} = 0.0365, and P_{orb} ~ 4.15 hr, existing numerical accretion models predict spin equilibrium condition with B ~ 30 MG if the mass ratio of the binary components is q_1 ~ 0.4. The high magnetic field makes V405 Aur a likely candidate as a progenitor of a polar.Comment: To appear in The Astrophysical Journal, September 1 Issue (2008), 9 pages, 10 figure

A Physical Model of Warped Galaxy Disks

Warped H I gas layers in the outer regions of spiral galaxies usually display a noticeably twisted structure. This structure almost certainly arises primarily as a result of differential precession in the H I disk as it settles toward a preferred orientation in an underlying dark halo potential well that is not spherically symmetric. In an attempt to better understand the structure and evolution of these twisted, warped disk structures, we have adopted the "twist-equation" formalism originally developed by Petterson (1977) to study accretion onto compact objects. Utilizing more recent treatments of this formalism, we have generalized the twist-equation to allow for the treatment of non-Keplerian disks and from it have derived a steady-state structure of twisted disks that develops from free precession in a nonspherical, logarithmic halo potential. We have used this steady-state solution to produce H I maps of five galaxies (M83, NGC 300, NGC 2841, NGC 5033, NGC5055), which match the general features of the observed maps of these galaxies quite well. In addition, the model provides an avenue through which the kinematical viscosity of the H I disk and the quadrupole distortion of the dark halo in each galaxy can be quantified. This generalized equation can also be used to examine the time-evolutionary behavior of warped galaxy disks.Comment: 22 pages, 5 figures, to be published in The Astrophysical Journa

Top-down contingent feature-specific orienting with and without awareness of the visual input

In the present article, the role of endogenous feature-specific orienting for conscious and unconscious vision is reviewed. We start with an overview of orienting. We proceed with a review of masking research, and the definition of the criteria of experimental protocols that demonstrate endogenous and exogenous orienting, respectively. Against this background of criteria, we assess studies of unconscious orienting and come to the conclusion that so far studies of unconscious orienting demonstrated endogenous feature-specific orienting. The review closes with a discussion of the role of unconscious orienting in action control

Dos and don’ts in response priming research

Response priming is a well-understood but sparsely employed paradigm in cognitive science. The method is powerful and well-suited for exploring early visuomotor processing in a wide range of tasks and research fields. Moreover, response priming can be dissociated from visual awareness, possibly because it is based on the first sweep of feedforward processing of primes and targets. This makes it a theoretically interesting device for separating conscious and unconscious vision. We discuss the major opportunities of the paradigm and give specific recommendations (e.g., tracing the time-course of priming in parametric experiments). Also, we point out typical confounds, design flaws, and data processing artifacts