The Physics of Wind-Fed Accretion

We provide a brief review of the physical processes behind the radiative driving of the winds of OB stars and the Bondi-Hoyle-Lyttleton capture and accretion of a fraction of the stellar wind by a compact object, typically a neutron star, in detached high-mass X-ray binaries (HMXBs). In addition, we describe a program to develop global models of the radiatively-driven photoionized winds and accretion flows of HMXBs, with particular attention to the prototypical system Vela X-1. The models combine XSTAR photoionization calculations, HULLAC emission models appropriate to X-ray photoionized plasmas, improved models of the radiative driving of photoionized winds, FLASH time-dependent adaptive-mesh hydrodynamics calculations, and Monte Carlo radiation transport. We present two- and three-dimensional maps of the density, temperature, velocity, ionization parameter, and emissivity distributions of representative X-ray emission lines, as well as synthetic global Monte Carlo X-ray spectra. Such models help to better constrain the properties of the winds of HMXBs, which bear on such fundamental questions as the long-term evolution of these binaries and the chemical enrichment of the interstellar medium.Comment: 9 pages including 5 color encapsulated postscript figures; accepted for inclusion in the proceedings of "Cool Discs, Hot Flows: The Varying Faces of Accreting Compact Objects," ed. M. Axelsson (New York: AIP); minor revision which addresses the referee's comments; added Fig. 1 and removed Fig. 3 and the associated tex

Extremely narrow spectrum of GRB110920A: further evidence for localised, subphotospheric dissipation

Much evidence points towards that the photosphere in the relativistic outflow in GRBs plays an important role in shaping the observed MeV spectrum. However, it is unclear whether the spectrum is fully produced by the photosphere or whether a substantial part of the spectrum is added by processes far above the photosphere. Here we make a detailed study of the $\gamma-$ray emission from single pulse GRB110920A which has a spectrum that becomes extremely narrow towards the end of the burst. We show that the emission can be interpreted as Comptonisation of thermal photons by cold electrons in an unmagnetised outflow at an optical depth of $\tau \sim 20$. The electrons receive their energy by a local dissipation occurring close to the saturation radius. The main spectral component of GRB110920A and its evolution is thus, in this interpretation, fully explained by the emission from the photosphere including localised dissipation at high optical depths.Comment: 14 pages, 11 figures, accepted to MNRA

On stoichiometry and intermixing at the spinel/perovskite interface in CoFe2O4/BaTiO3 thin films

The performance of complex oxide heterostructures depends primarily on the interfacial coupling of the two component structures. This interface character inherently varies with the synthesis method and conditions used since even small composition variations can alter the electronic, ferroelectric, or magnetic functional properties of the system. The focus of this article is placed on the interface character of a pulsed laser deposited CoFe2O4/BaTiO3 thin film. Using a range of state-of-the-art transmission electron microscopy methodologies, the roles of substrate morphology, interface stoichiometry, and cation intermixing are determined on the atomic level. The results reveal a surprisingly uneven BaTiO3 substrate surface formed after the film deposition and Fe atom incorporation in the top few monolayers inside the unit cell of the BaTiO3 crystal. Towards the CoFe2O4 side, a disordered region extending several nanometers from the interface was revealed and both Ba and Ti from the substrate were found to diffuse into the spinel layer. The analysis also shows that within this somehow incompatible composite interface, a different phase is formed corresponding to the compound Ba2Fe3Ti5O15, which belongs to the ilmenite crystal structure of FeTiO3 type. The results suggest a chemical activity between these two oxides, which could lead to the synthesis of complex engineered interfaces

Weighted maximal regularity estimates and solvability of non-smooth elliptic systems II

We continue the development, by reduction to a first order system for the conormal gradient, of $L^2$ \textit{a priori} estimates and solvability for boundary value problems of Dirichlet, regularity, Neumann type for divergence form second order, complex, elliptic systems. We work here on the unit ball and more generally its bi-Lipschitz images, assuming a Carleson condition as introduced by Dahlberg which measures the discrepancy of the coefficients to their boundary trace near the boundary. We sharpen our estimates by proving a general result concerning \textit{a priori} almost everywhere non-tangential convergence at the boundary. Also, compactness of the boundary yields more solvability results using Fredholm theory. Comparison between classes of solutions and uniqueness issues are discussed. As a consequence, we are able to solve a long standing regularity problem for real equations, which may not be true on the upper half-space, justifying \textit{a posteriori} a separate work on bounded domains.Comment: 76 pages, new abstract and few typos corrected. The second author has changed nam

One step multiderivative methods for first order ordinary differential equations

A family of one-step multiderivative methods based on Padé approximants to the exponential function is developed. The methods are extrapolated and analysed for use in PECE mode. Error constants and stability intervals are calculated and the combinations compared with well known linear multi-step combinations and combinations using high accuracy Newton-Cotes quadrature formulas as correctors. w926020

Parallel Unsmoothed Aggregation Algebraic Multigrid Algorithms on GPUs

We design and implement a parallel algebraic multigrid method for isotropic graph Laplacian problems on multicore Graphical Processing Units (GPUs). The proposed AMG method is based on the aggregation framework. The setup phase of the algorithm uses a parallel maximal independent set algorithm in forming aggregates and the resulting coarse level hierarchy is then used in a K-cycle iteration solve phase with a $\ell^1$-Jacobi smoother. Numerical tests of a parallel implementation of the method for graphics processors are presented to demonstrate its effectiveness.Comment: 18 pages, 3 figure

Visibilia ex invisibilibus: seeing at the nanoscale for improved preservation of parchment

This paper describes the application of atomic force microscopy (AFM) for the imaging of collagen denaturation as observed in parchment. Parchment is prepared from processed animal skin and collagen is the main component. Large collections in national archives, libraries and religious institutions contain numerous documents written on parchment. Their preservation presents an unsolved problem for conservators. The main challenge is to assess the state of collagen and to detect what conservators refer to as the pre-gelatinised state, which can cause surface cracking resulting in a loss of text and can increase the vulnerability of parchment to aqueous cleaning agents. Atomic force microscopy (AFM) was first used within the Improved Damage Assessment of Parchment (IDAP) project, enabling the characterisation of the collagen structure within parchment at the nanoscale. Damage categories were also established based on the extent of the ordered collagen structure that was observed in the AFM images. This paper describes the work following the IDAP project, where morphological changes in the fibres due to both artificial and natural ageing were observed and linked to observations made by AFM. It also explores the merits and drawbacks of different approaches used for sample preparation and the possibility of using a portable AFM for imaging directly on the surface of documents. A case study on a manuscript from the 18th century is presented

A multidomain spectral method for solving elliptic equations

We present a new solver for coupled nonlinear elliptic partial differential equations (PDEs). The solver is based on pseudo-spectral collocation with domain decomposition and can handle one- to three-dimensional problems. It has three distinct features. First, the combined problem of solving the PDE, satisfying the boundary conditions, and matching between different subdomains is cast into one set of equations readily accessible to standard linear and nonlinear solvers. Second, touching as well as overlapping subdomains are supported; both rectangular blocks with Chebyshev basis functions as well as spherical shells with an expansion in spherical harmonics are implemented. Third, the code is very flexible: The domain decomposition as well as the distribution of collocation points in each domain can be chosen at run time, and the solver is easily adaptable to new PDEs. The code has been used to solve the equations of the initial value problem of general relativity and should be useful in many other problems. We compare the new method to finite difference codes and find it superior in both runtime and accuracy, at least for the smooth problems considered here.Comment: 31 pages, 8 figure

Searching for hidden mirror symmetries in CMB fluctuations from WMAP 7 year maps

We search for hidden mirror symmetries at large angular scales in the WMAP 7 year Internal Linear Combination map of CMB temperature anisotropies using global pixel based estimators introduced for this aim. Two different axes are found for which the CMB intensity pattern is anomalously symmetric (or anti-symmetric) under reflection with respect to orthogonal planes at the 99.84(99.96)% CL (confidence level), if compared to a result for an arbitrary axis in simulations without the symmetry. We have verified that our results are robust to the introduction of the galactic mask. The direction of such axes is close to the CMB kinematic dipole and nearly orthogonal to the ecliptic plane, respectively. If instead the real data are compared to those in simulations taken with respect to planes for which the maximal mirror symmetry is generated by chance, the confidence level decreases to 92.39 (76.65)%. But when the effect in question translates into the anomalous alignment between normals to planes of maximal mirror (anti)-symmetry and these natural axes mentioned. We also introduce the representation of the above estimators in the harmonic domain, confirming the results obtained in the pixel one. The symmetry anomaly is shown to be almost entirely due to low multipoles, so it may have a cosmological and even primordial origin. Contrary, the anti-symmetry one is mainly due to intermediate multipoles that probably suggests its non-fundamental nature. We have demonstrated that these anomalies are not connected to the known issue of the low variance in WMAP observations and we have checked that axially symmetric parts of these anomalies are small, so that the axes are not the symmetry ones.Comment: 18 pages, 10 figures, 2 tables. Consideration and discussion expanded, 5 figures and 1 table added, main conclusions unchange