A Robust Determination of the size of quasar accretion disks using gravitational microlensing

Using microlensing measurements from a sample of 27 image-pairs of 19 lensed quasars we determine a maximum likelihood estimate for the accretion disk size of an {{\em}average} quasar of $r_s=4.0^{+2.4}_{-3.1}$ light days at rest frame $=1736$\AA\ for microlenses with a mean mass of $=0.3M_\odot$. This value, in good agreement with previous results from smaller samples, is roughly a factor of 5 greater than the predictions of the standard thin disk model. The individual size estimates for the 19 quasars in our sample are also in excellent agreement with the results of the joint maximum likelihood analysis.Comment: 6 pages, 3 figures, submitted to Ap

A GPU-Enabled, High-Resolution Cosmological Microlensing Parameter Survey

In the era of synoptic surveys, the number of known gravitationally lensed quasars is set to increase by over an order of magnitude. These new discoveries will enable a move from single-quasar studies to investigations of statistical samples, presenting new opportunities to test theoretical models for the structure of quasar accretion discs and broad emission line regions (BELRs). As one crucial step in preparing for this influx of new lensed systems, a large-scale exploration of microlensing convergence-shear parameter space is warranted, requiring the computation of O(10^5) high resolution magnification maps. Based on properties of known lensed quasars, and expectations from accretion disc/BELR modelling, we identify regions of convergence-shear parameter space, map sizes, smooth matter fractions, and pixel resolutions that should be covered. We describe how the computationally time-consuming task of producing ~290000 magnification maps with sufficient resolution (10000^2 pixel/map) to probe scales from the inner edge of the accretion disc to the BELR can be achieved in ~400 days on a 100 teraflop/s high performance computing facility, where the processing performance is achieved with graphics processing units. We illustrate a use-case for the parameter survey by investigating the effects of varying the lens macro-model on accretion disc constraints in the lensed quasar Q2237+0305. We find that although all constraints are consistent within their current error bars, models with more densely packed microlenses tend to predict shallower accretion disc radial temperature profiles. With a large parameter survey such as the one described here, such systematics on microlensing measurements could be fully explored.Comment: 30 pages, 3 figures, accepted for publication in Ap

Continuous and discontinuous modelling of ductile fracture

In many metal forming processes (e.g. blanking, trimming, clipping, machining, cutting) fracture is triggered in order to induce material separation along with a desired product geometry. This type of fracture is preceded by a certain amount of plastic deformation and requires additional energy to be supplied in order for the crack to propagate. It is known as ductile fracture, as opposed to brittle fracture (as in ceramics, concrete, etc). Ductile fracture originates at a microscopic level, as the result of voids initiated at inclusions in the material matrix. These microscopic degradation processes lead to the degradation of the macroscopic mechanical properties, causing softening, strain localisation and finally the formation of macroscopic cracks. The initiation and propagation of cracks has traditionally been studied by fracture mechanics. Yet, the application of this theory to ductile fracture, where highly nonlinear degradation processes (material and geometrical) take place in the fracture process zone, raises many questions. To model these processes, continuum models can be used, either in the form of softening plasticity or continuum damage mechanics. Yet, continuous models can not be applied to model crack propagation, because displacements are no longer continuous across the crack. Hence, a proper model for ductile fracture requires a different approach, one that combines a continuous softening model with a strategy to represent cracks, i.e. displacement discontinuities. This has been the main goal of the present work. In a combined approach, the direction of crack propagation is automatically determined by the localisation pattern, and its rate strongly depends on the evolution of damage (or other internal variables responsible for the strain softening). This contrasts with fracture mechanics, where the material behaviour is not directly linked to the crack propagation criteria. Softening materials have to be supplied with an internal length, which acts as a localisation limiter, thereby ensuring the well-posedness of the governing partial differential equations and mesh independent results. For this purpose, a nonlocal gradient enhancement has been used in this work, which gives similar results to nonlocal models of an integral form. A number of numerical methods are available to model displacement discontinuities in a continuum. In the present context, we have used a remeshing strategy, since it has additional advantages when used with large strain localising material models: it prevents excessive element distortions and allows to optimise the element distriviii bution through mesh adaptivity. As a first step towards a continuum-discontinuum approach, an uncoupled damage model is used first, in which damage merely acts as a crack initiation-propagation indicator, without causing material softening. Since uncoupled models do not lead to material localisation, no regularisation is needed. Yet, uncoupled approaches can not capture the actual failure mechanisms and therefore, in general, can give reliable results only when the size of the fracture process zone is so small that its effect can be neglected. When the size of the fracture process zone is large enough, a truly combined model must be used, which is developed in the second part of this study. Due to softening, the transition from the continuous damage material to the discrete crack occurs gradually, with little stress redistribution, in contrast with the previous uncoupled approach. The gradient regularised softening behaviour is introduced in the yield behaviour of an elastoplastic material. The combined model has been applied satisfactorily to the prediction of ductile failure under shear loading conditions. Third, to be able to apply the model to more general loading conditions, the material description has been improved by introducing the influence of stress triaxiality in the damage evolution of a gradient regularised elastoplastic damage model. The model has been obtained using the continuum damage mechanics concept of effective stress. Results show how compressive (tensile) states of triaxiality may increase (decrease) the material ductility. Finally, the combined approach is applied to the modelling of actual metal forming processes, e.g. blanking, fine blanking, score forming. The gradient regularisation has been implemented in an operator-split manner, which can be very appealing for engineering purposes. To capture the large strain gradients in the localisation zones, a new mesh adaptivity criterion has been proposed. The results of the simulations are in good agreement with experimental data from literature

Detection of chromatic microlensing in Q 2237+0305 A

We present narrowband images of the gravitational lens system Q~2237+0305 made with the Nordic Optical Telescope in eight different filters covering the wavelength interval 3510-8130 \AA. Using point-spread function photometry fitting we have derived the difference in magnitude versus wavelength between the four images of Q~2237+0305. At $\lambda=4110$ \AA, the wavelength range covered by the Str\"omgren-v filter coincides with the position and width of the CIV emission line. This allows us to determine the existence of microlensing in the continuum and not in the emission lines for two images of the quasar. Moreover, the brightness of image A shows a significant variation with wavelength which can only be explained as consequence of chromatic microlensing. To perform a complete analysis of this chromatic event our observations were used together with Optical Gravitational Lensing Experiment light curves. Both data sets cannot be reproduced by the simple phenomenology described under the caustic crossing approximation; using more realistic representations of microlensing at high optical depth, we found solutions consistent with simple thin disk models ($r_{s}\varpropto \lambda^{4/3}$); however, other accretion disk size-wavelength relationships also lead to good solutions. New chromatic events from the ongoing narrow band photometric monitoring of Q~2237+0305 are needed to accurately constrain the physical properties of the accretion disk for this system.Comment: 9 pages, 9 figures, 2 tables. Matches ApJ published version. Some references adde