Polarization Modeling and Predictions for DKIST Part 2: Application of the Berreman Calculus to Spectral Polarization Fringes of Beamsplitters and Crystal Retarders

We outline polarization fringe predictions derived from a new application of the Berreman calculus for the Daniel K. Inouye Solar Telescope (DKIST) retarder optics. The DKIST retarder baseline design used 6 crystals, single-layer anti-reflection coatings, thick cover windows and oil between all optical interfaces. This new tool estimates polarization fringes and optic Mueller matrices as functions of all optical design choices. The amplitude and period of polarized fringes under design changes, manufacturing errors, tolerances and several physical factors can now be estimated. This tool compares well with observations of fringes for data collected with the SPINOR spectropolarimeter at the Dunn Solar Telescope using bi-crystalline achromatic retarders as well as laboratory tests. With this new tool, we show impacts of design decisions on polarization fringes as impacted by anti-reflection coatings, oil refractive indices, cover window presence and part thicknesses. This tool helped DKIST decide to remove retarder cover windows and also recommends reconsideration of coating strategies for DKIST. We anticipate this tool to be essential in designing future retarders for mitigation of polarization and intensity fringe errors in other high spectral resolution astronomical systems.Comment: Accepted for publication in JATI

Representations of world coordinates in FITS

The initial descriptions of the FITS format provided a simplified method for describing the physical coordinate values of the image pixels, but deliberately did not specify any of the detailed conventions required to convey the complexities of actual image coordinates. Building on conventions in wide use within astronomy, this paper proposes general extensions to the original methods for describing the world coordinates of FITS data. In subsequent papers, we apply these general conventions to the methods by which spherical coordinates may be projected onto a two-dimensional plane and to frequency/wavelength/velocity coordinates.Comment: 15 Pages, 1 figure, LaTex with Astronomy & Astrophysics macro package, submitted to A&A, related papers at http://www.aoc.nrao.edu/~egreise

An open and parallel multiresolution framework using block-based adaptive grids

A numerical approach for solving evolutionary partial differential equations in two and three space dimensions on block-based adaptive grids is presented. The numerical discretization is based on high-order, central finite-differences and explicit time integration. Grid refinement and coarsening are triggered by multiresolution analysis, i.e. thresholding of wavelet coefficients, which allow controlling the precision of the adaptive approximation of the solution with respect to uniform grid computations. The implementation of the scheme is fully parallel using MPI with a hybrid data structure. Load balancing relies on space filling curves techniques. Validation tests for 2D advection equations allow to assess the precision and performance of the developed code. Computations of the compressible Navier-Stokes equations for a temporally developing 2D mixing layer illustrate the properties of the code for nonlinear multi-scale problems. The code is open source

A rarefaction-tracking method for hyperbolic conservation laws

We present a numerical method for scalar conservation laws in one space dimension. The solution is approximated by local similarity solutions. While many commonly used approaches are based on shocks, the presented method uses rarefaction and compression waves. The solution is represented by particles that carry function values and move according to the method of characteristics. Between two neighboring particles, an interpolation is defined by an analytical similarity solution of the conservation law. An interaction of particles represents a collision of characteristics. The resulting shock is resolved by merging particles so that the total area under the function is conserved. The method is variation diminishing, nevertheless, it has no numerical dissipation away from shocks. Although shocks are not explicitly tracked, they can be located accurately. We present numerical examples, and outline specific applications and extensions of the approach.Comment: 21 pages, 7 figures. Similarity 2008 conference proceeding

The impact of organisational external peer review on colorectal cancer treatment and survival in the Netherlands

Background: Organisational external peer review was introduced in 1994 in the Netherlands to improve multidisciplinary cancer care. We examined the clinical impact of this programme on colorectal cancer care. Methods: Patients with primary colorectal cancer were included from 23 participating hospitals and 7 controls. Hospitals from the intervention group were dichotomised by their implementation proportion (IP) of the recommendations from each peer review (high IP vs low IP). Outcome measures were the introduction of new multidisciplinary therapies and survival. Results: In total, 45 705 patients were included (1990-2010). Patients from intervention hospitals more frequently received adjuvant chemotherapy for stage III colon cancer. T2-3/M0 rectal cancer patients from hospitals with a high IP had a higher chance of receiving preoperative radiotherapy (OR 1.31, 95% CI 1.11-1.55) compared with the controls and low IP group (OR 0.75, 95% CI 0.63-0.88). There were no differences in the use of preoperative chemoradiation for T4/M0 rectal cancer. Survival was slightly higher in colon cancer patients from intervention hospitals but unrelated to the phase of the programme in which the hospital was at the time of diagnosis. Conclusions: Some positive effects of external peer review on cancer care were found, but the results need to be interpreted cautiously due to the ambiguity of the outcomes and possible confounding factors

On the validity of mean-field amplitude equations for counterpropagating wavetrains

We rigorously establish the validity of the equations describing the evolution of one-dimensional long wavelength modulations of counterpropagating wavetrains for a hyperbolic model equation, namely the sine-Gordon equation. We consider both periodic amplitude functions and localized wavepackets. For the localized case, the wavetrains are completely decoupled at leading order, while in the periodic case the amplitude equations take the form of mean-field (nonlocal) Schr\"odinger equations rather than locally coupled partial differential equations. The origin of this weakened coupling is traced to a hidden translation symmetry in the linear problem, which is related to the existence of a characteristic frame traveling at the group velocity of each wavetrain. It is proved that solutions to the amplitude equations dominate the dynamics of the governing equations on asymptotically long time scales. While the details of the discussion are restricted to the class of model equations having a leading cubic nonlinearity, the results strongly indicate that mean-field evolution equations are generic for bimodal disturbances in dispersive systems with \O(1) group velocity.Comment: 16 pages, uuencoded, tar-compressed Postscript fil

Kink instabilities in jets from rotating magnetic fields

We have performed 2.5D and 3D simulations of conical jets driven by the rotation of an ordered, large-scale magnetic field in a stratified atmosphere. The simulations cover about three orders of magnitude in distance to capture the centrifugal acceleration as well as the evolution past the Alfven surface. We find that the jets develop kink instabilities, the characteristics of which depend on the velocity profile imposed at the base of the flow. The instabilities are especially pronounced with a rigid rotation profile, which induces a shearless magnetic field. The jet's expansion appears to be limiting the growth of Alfven mode instabilities.Comment: 10 pages, 13 figures, accepted for publication in A&

Relaxation and reconstruction on (111) surfaces of Au, Pt, and Cu

We have theoretically studied the stability and reconstruction of (111) surfaces of Au, Pt, and Cu. We have calculated the surface energy, surface stress, interatomic force constants, and other relevant quantities by ab initio electronic structure calculations using the density functional theory (DFT), in a slab geometry with periodic boundary conditions. We have estimated the stability towards a quasi-one-dimensional reconstruction by using the calculated quantities as parameters in a one-dimensional Frenkel-Kontorova model. On all surfaces we have found an intrinsic tensile stress. This stress is large enough on Au and Pt surfaces to lead to a reconstruction in which a denser surface layer is formed, in agreement with experiment. The experimentally observed differences between the dense reconstruction pattern on Au(111) and a sparse structure of stripes on Pt(111) are attributed to the details of the interaction potential between the first layer of atoms and the substrate.Comment: 8 pages, 3 figures, submitted to Physical Review

Phase Slips and the Eckhaus Instability

We consider the Ginzburg-Landau equation, $\partial_t u= \partial_x^2 u + u - u|u|^2$, with complex amplitude $u(x,t)$. We first analyze the phenomenon of phase slips as a consequence of the {\it local} shape of $u$. We next prove a {\it global} theorem about evolution from an Eckhaus unstable state, all the way to the limiting stable finite state, for periodic perturbations of Eckhaus unstable periodic initial data. Equipped with these results, we proceed to prove the corresponding phenomena for the fourth order Swift-Hohenberg equation, of which the Ginzburg-Landau equation is the amplitude approximation. This sheds light on how one should deal with local and global aspects of phase slips for this and many other similar systems.Comment: 22 pages, Postscript, A

Scalar field induced oscillations of neutron stars and gravitational collapse

We study the interaction of massless scalar fields with self-gravitating neutron stars by means of fully dynamic numerical simulations of the Einstein-Klein-Gordon perfect fluid system. Our investigation is restricted to spherical symmetry and the neutron stars are approximated by relativistic polytropes. Studying the nonlinear dynamics of isolated neutron stars is very effectively performed within the characteristic formulation of general relativity, in which the spacetime is foliated by a family of outgoing light cones. We are able to compactify the entire spacetime on a computational grid and simultaneously impose natural radiative boundary conditions and extract accurate radiative signals. We study the transfer of energy from the scalar field to the fluid star. We find, in particular, that depending on the compactness of the neutron star model, the scalar wave forces the neutron star either to oscillate in its radial modes of pulsation or to undergo gravitational collapse to a black hole on a dynamical timescale. The radiative signal, read off at future null infinity, shows quasi-normal oscillations before the setting of a late time power-law tail.Comment: 12 pages, 13 figures, submitted to Phys. Rev.