The anatomy of a quadruply imaged gravitational lens system

The key to using a strong gravitational lens system to measure the Hubble constant is to obtain an accurate model of the lens potential. In this paper, we investigate the properties of gravitational lens B1608+656, a quadruply-imaged lens system with an extended source intensity distribution. Our analysis is valid for generic quadruply-lensed systems. Limit curves and isophotal separatrices are defined for such systems, and we show that the isophotal separatrices must intersect at the critical curves and the satellite isophotes must be tangent to the limit curves. The current model of B1608+656 (Koopmans et al. 2003) satisfies these criteria for some, but not all, of the isophotal separatrices within the observational uncertainty. We study a non-parametric method of potential reconstruction proposed by Blandford, Surpi & Kundic (2001) and demonstrate that although the method works in principle and elucidates image formation, the initial potential only converges to the true model when it is within ~ 1 percent of the true model.Comment: 12 pages, 12 figures. Minor revisions based on referee's comments after initial submission to MNRA

Computational modelling and experimental characterisation of heterogeneous materials

Heterogeneous materials can exhibit behaviour under load that cannot be described by classical continuum elasticity. Beams in bending can show a relative stiffening as the beam depth tends to zero, a size effect. Size effects are recognised in higher order continuum elastic theories such as micropolar elasticity. The drawback of higher order theories is the requirement of addition constitutive relations and associated properties that are often difficult to establish experimentally. Furthermore the finite element method, of great benefit in classical elasticity, has shown limitations when applied to micropolar elasticity. The determination of additional constitutive properties and the computational modelling of micropolar elasticity will be discussed in the context of a model heterogeneous material loaded in simple 3 point bending. The model material was created by drilling holes in aluminium bar in a regular pattern, with the hole axis normal to the plane of bending. The bending tests show that a size effect is present. These results are compared against modelling the detailed beam geometries in the finite element package ANSYS, which again shows the size effect. These two bending test are used to extract the additional micropolar elastic material properties. A comparison is then made against analytical solutions,numerical solutions using a micropolar beam finite element and a micropolar plane stress control volume method.It will be shown that the need for extensive experimental testing to determine the additional constitutive properties may not be necessary with the appropriate use of numerical methods

A Bayesian analysis of regularised source inversions in gravitational lensing

Strong gravitational lens systems with extended sources are of special interest because they provide additional constraints on the models of the lens systems. To use a gravitational lens system for measuring the Hubble constant, one would need to determine the lens potential and the source intensity distribution simultaneously. A linear inversion method to reconstruct a pixellated source brightness distribution of a given lens potential model was introduced by Warren & Dye. In the inversion process, a regularisation on the source intensity is often needed to ensure a successful inversion with a faithful resulting source. In this paper, we use Bayesian analysis to determine the optimal regularisation constant (strength of regularisation) of a given form of regularisation and to objectively choose the optimal form of regularisation given a selection of regularisations. We consider and compare quantitatively three different forms of regularisation previously described in the literature for source inversions in gravitational lensing: zeroth-order, gradient and curvature. We use simulated data with the exact lens potential to demonstrate the method. We find that the preferred form of regularisation depends on the nature of the source distribution.Comment: 18 pages, 10 figures; Revisions based on referee's comments after initial submission to MNRA

Effects of exercise vs experimental osteoarthritis on imaging outcomes

SummaryObjectiveTo identify changes in imaging outcomes in a controlled model of osteoarthritis (OA) vs exercise.MethodSixteen 2-year-old horses were randomly assigned to an exercise control (n=8) or an exercise OA (n=8) group. All horses had middle carpal joints arthroscopically explored and an osteochondral fragment was induced in one middle carpal joint of the OA group. All horses were treadmill exercised for the duration of the study (91 days). Clinical, radiographic, nuclear scintigraphic, computed tomographic and magnetic resonance imaging (MRI) examinations were performed and outcomes of these were compared between groups. Imaging results were correlated to clinical, biomarker and gross pathologic results.ResultsThe OA group had significant increases in clinical outcomes and most imaging parameters. Specifically, the OA group showed significant increases in radiographic lysis and nuclear scintigraphic uptake. There was very little change in subchondral bone density, but a significant change in subchondral bone edema. Radiographic lysis, radial carpal bone edema and nuclear scintigraphy were strongly correlated with clinical changes and radial carpal bone edema was strongly correlated with changes in Type I and Type II collagen found in the synovial fluid.ConclusionsOA induced significant changes in imaging parameters beyond the adaptation seen with exercise. Bone edema detected with MRI was closely correlated with collagen biomarkers detected in the synovial fluid

First-principles study of (BiScO3){1-x}-(PbTiO3){x} piezoelectric alloys

We report a first-principles study of a class of (BiScO3)_{1-x}-(PbTiO3)_x (BS-PT) alloys recently proposed by Eitel et al. as promising materials for piezoelectric actuator applications. We show that (i) BS-PT displays very large structural distortions and polarizations at the morphotropic phase boundary (MPB) (we obtain a c/a of ~1.05-1.08 and P_tet of ~1.1 C/m^2); (ii) the ferroelectric and piezoelectric properties of BS-PT are dominated by the onset of hybridization between Bi/Pb-6p and O-2p orbitals, a mechanism that is enhanced upon substitution of Pb by Bi; and (iii) the piezoelectric responses of BS-PT and Pb(Zr_{1-x}Ti_x)O3 (PZT) at the MPB are comparable, at least as far as the computed values of the piezoelectric coefficient d_15 are concerned. While our results are generally consistent with experiment, they also suggest that certain intrinsic properties of BS-PT may be even better than has been indicated by experiments to date. We also discuss results for PZT that demonstrate the prominent role played by Pb displacements in its piezoelectric properties.Comment: 6 pages, with 3 postscript figures embedded. Uses REVTEX and epsf macros. Also available at http://www.physics.rutgers.edu/~dhv/preprints/ji_bi/index.htm

Current reversal with type-I intermittency in deterministic inertia ratchets

The intermittency is investigated when the current reversal occurs in a deterministic inertia ratchet system. To determine which type the intermittency belongs to, we obtain the return map of velocities of particle using stroboscopic recording, and numerically calculate the distribution of average laminar length ${}$. The distribution follows the scaling law of ${} \propto {\epsilon}^{-1/2}$, the characteristic relation of type-I intermittency.Comment: 4 pages, 7 figure

Porting Decision Tree Algorithms to Multicore using FastFlow

The whole computer hardware industry embraced multicores. For these machines, the extreme optimisation of sequential algorithms is no longer sufficient to squeeze the real machine power, which can be only exploited via thread-level parallelism. Decision tree algorithms exhibit natural concurrency that makes them suitable to be parallelised. This paper presents an approach for easy-yet-efficient porting of an implementation of the C4.5 algorithm on multicores. The parallel porting requires minimal changes to the original sequential code, and it is able to exploit up to 7X speedup on an Intel dual-quad core machine.Comment: 18 pages + cove

Anomalous enhancement of tetragonality in PbTiO3 induced by negative pressure

Using a first-principles approach based on density-functional theory, we find that a large tetragonal strain can be induced in PbTiO3 by application of a negative hydrostatic pressure. The structural parameters and the dielectric and dynamical properties are found to change abruptly near a crossover pressure, displaying a ``kinky'' behavior suggestive of proximity to a phase transition. Analogous calculations for BaTiO3 show that the same effect is also present there, but at much higher negative pressure. We investigate this unexpected behavior of PbTiO3 and discuss an interpretation involving a phenomenological description in terms of a reduced set of relevant degrees of freedom.Comment: 9 pages, with 9 postscript figures embedded. Uses REVTEX and epsf macros. Also available at http://www.physics.rutgers.edu/~dhv/preprints/st_pbti/index.htm

Synchronization in a System of Globally Coupled Oscillators with Time Delay

We study the synchronization phenomena in a system of globally coupled oscillators with time delay in the coupling. The self-consistency equations for the order parameter are derived, which depend explicitly on the amount of delay. Analysis of these equations reveals that the system in general exhibits discontinuous transitions in addition to the usual continuous transition, between the incoherent state and a multitude of coherent states with different synchronization frequencies. In particular, the phase diagram is obtained on the plane of the coupling strength and the delay time, and ubiquity of multistability as well as suppression of the synchronization frequency is manifested. Numerical simulations are also performed to give consistent results

Epitaxially strained [001]-(PbTiO3_3)1_1(PbZrO3_3)1_1 superlattice and PbTiO3_3 from first principles

The effect of layer-by-layer heterostructuring and epitaxial strain on lattice instabilities and related ferroelectric properties is investigated from first principles for the [001]-(PbTiO$_3$)$_1$(PbZrO$_3$)$_1$ superlattice and pure PbTiO$_3$ on a cubic substrate. The results for the superlattice show an enhancement of the stability of the monoclinic r-phase with respect to pure PbTiO$_3$. Analysis of the lattice instabilities of the relaxed centrosymmetric reference structure computed within density functional perturbation theory suggests that this results from the presence of two unstable zone-center modes, one confined in the PbTiO$_3$ layer and one in the PbZrO$_3$ layer, which produce in-plane and normal components of the polarization, respectively. The zero-temperature dielectric response is computed and shown to be enhanced not only near the phase boundaries, but throughout the r-phase. Analysis of the analogous calculation for pure PbTiO$_3$ is consistent with this interpretation, and suggests useful approaches to engineering the dielectric properties of artificially structured perovskite oxides.Comment: 8 pages, 5 figure